The DNA sequence and biological annotation of human chromosome 1

The reference sequence for each human chromosome provides the framework for understanding genome function, variation and evolution. Here we report the finished sequence and biological annotation of human chromosome 1. Chromosome 1 is gene-dense, with 3,141 genes and 991 pseudogenes, and many coding sequences overlap. Rearrangements and mutations of chromosome 1 are prevalent in cancer and many other diseases. Patterns of sequence variation reveal signals of recent selection in specific genes that may contribute to human fitness, and also in regions where no function is evident. Fine-scale recombination occurs in hotspots of varying intensity along the sequence, and is enriched near genes. These and other studies of human biology and disease encoded within chromosome 1 are made possible with the highly accurate annotated sequence, as part of the completed set of chromosome sequences that comprise the reference human genome.     

Functional genetic analysis of mouse chromosome 11

Now that the mouse and human genome sequences are complete, biologists need systematic approaches to determine the function of each gene. A powerful way to discover gene function is to determine the consequence of mutations in living organisms. Large-scale production of mouse mutations with the point mutagen N-ethyl-N-nitrosourea (ENU) is a key strategy for analysing the human genome because mouse mutants will reveal functions unique to mammals, and many may model human diseases. To examine genes conserved between human and mouse, we performed a recessive ENU mutagenesis screen that uses a balancer chromosome, inversion chromosome 11 (refs 4, 5). Initially identified in the fruitfly, balancer chromosomes are valuable genetic tools that allow the easy isolation of mutations on selected chromosomes. Here we show the isolation of 230 new recessive mouse mutations, 88 of which are on chromosome 11. This genetic strategy efficiently generates and maps mutations on a single chromosome, even as mutations throughout the genome are discovered. The mutations reveal new defects in haematopoiesis, craniofacial and cardiovascular development, and fertility.     

The dynamics of chromosome evolution in birds and mammals

Comparative mapping, which compares the location of homologous genes in different species, is a powerful tool for studying genome evolution. Comparative maps suggest that rates of chromosomal change in mammals can vary from one to ten rearrangements per million years. On the basis of these rates we would expect 84 to 600 conserved segments in a chicken comparison with human or mouse. Here we build comparative maps between these species and estimate that numbers of conserved segments are in the lower part of this range. We conclude that the organization of the human genome is closer to that of the chicken than the mouse and by adding comparative mapping results from a range of vertebrates, we identify three possible phases of chromosome evolution. The relative stability of genomes such as those of the chicken and human will enable the reconstruction of maps of ancestral vertebrates.     

Comparative analyses of multi-species sequences from targeted genomic regions

The systematic comparison of genomic sequences from different organisms represents a central focus of contemporary genome analysis. Comparative analyses of vertebrate sequences can identify coding and conserved non-coding regions, including regulatory elements, and provide insight into the forces that have rendered modern-day genomes. As a complement to whole-genome sequencing efforts, we are sequencing and comparing targeted genomic regions in multiple, evolutionarily diverse vertebrates. Here we report the generation and analysis of over 12 megabases (Mb) of sequence from 12 species, all derived from the genomic region orthologous to a segment of about 1.8 Mb on human chromosome 7 containing ten genes, including the gene mutated in cystic fibrosis. These sequences show conservation reflecting both functional constraints and the neutral mutational events that shaped this genomic region. In particular, we identify substantial numbers of conserved non-coding segments beyond those previously identified experimentally, most of which are not detectable by pair-wise sequence comparisons alone. Analysis of transposable element insertions highlights the variation in genome dynamics among these species and confirms the placement of rodents as a sister group to the primates.     

DNA sequence and analysis of human chromosome 18

Chromosome 18 appears to have the lowest gene density of any human chromosome and is one of only three chromosomes for which trisomic individuals survive to term. There are also a number of genetic disorders stemming from chromosome 18 trisomy and aneuploidy. Here we report the finished sequence and gene annotation of human chromosome 18, which will allow a better understanding of the normal and disease biology of this chromosome. Despite the low density of protein-coding genes on chromosome 18, we find that the proportion of non-protein-coding sequences evolutionarily conserved among mammals is close to the genome-wide average. Extending this analysis to the entire human genome, we find that the density of conserved non-protein-coding sequences is largely uncorrelated with gene density. This has important implications for the nature and roles of non-protein-coding sequence elements.     

Comparative chromosome mapping of a conserved homoeo box region in mouse and human

Specific genes are assumed to regulate pattern formation in the mammalian embryo, but as yet none has been identified unequivocally. It is possible that such genes in mammals may be identified by virtue of a conserved coding sequence, because many of the Drosophila melanogaster homoeotic and segmentation genes, which have crucial roles in the regulation of segmental pattern formation during embryonic development, contain a 180-base pair (bp) DNA sequence, the homoeo box, and that sequences homologous to the Drosophila homoeo box are also present in 6-10 copies in higher animals, including mammals. Although the assumption that the homoeo box identifies genes responsible for pattern formation in mammals remains to be validated, it is a particularly attractive hypothesis given the strong conservation of homoeo boxes over vast evolutionary distances. Here we report the localization of a human homoeo box region, previously cloned and shown to contain two homoeo boxes within a sequence of 5-kilobases (kb), to the long arm of chromosome 17. We show that two single-copy homoeo box-flanking probes derived from this region strongly hybridize to single-copy restriction fragments in mouse genomic DNA and that these conserved homoeo box-flanking sequences map to mouse chromosome 11. This may be significant as several genes that map to chromosome 17 in human also map to chromosome 11 in the mouse, implying that a segment of mouse chromosome 11 is homologous to a region of human chromosome 17. Taken together, these data suggest that the homoeo box region detected with our probes is highly conserved in human and mouse.     

A conserved sequence in the immunoglobulin J kappa-C kappa intron: possible enhancer element

Several functionally important genetic elements (such as the TATA box, mRNA splice sequences, poly(A) addition signal) were first detected as short segments of unexplained sequence homology within non-coding regions of different genes. A short region of unknown sequence in the intron between the human J kappa and C kappa immunoglobulin coding regions was found to be sufficiently homologous to the corresponding segment of the mouse gene to form stable heteroduplexes. Although no specific function has yet been definitely ascribed to this region (which we call the kappa intron conserved region, or KICR), some functional significance has been inferred from the findings that (1) activation of B lymphocytes induces a DNase hypersensitivity site in this region and (2) deletions including this region reduce expression of kappa genes introduced into lymphoid cells. To delineate the KICR more precisely and to test the generality of the sequence conservation in a third species, we have sequenced this region of the human and mouse genes and have examined the corresponding region of a recently cloned rabbit kappa gene. We find a segment of about 130 base pairs (bp) that shows striking conservation in all three species, demonstrating homology significantly higher than within the C kappa coding region itself. Two short sequences from the J kappa-C kappa intron that were noted by other investigators to be homologous to proposed 'enhancer' sequences both lie within the conserved region.     

In vivo enhancer analysis of human conserved non-coding sequences

Identifying the sequences that direct the spatial and temporal expression of genes and defining their function in vivo remains a significant challenge in the annotation of vertebrate genomes. One major obstacle is the lack of experimentally validated training sets. In this study, we made use of extreme evolutionary sequence conservation as a filter to identify putative gene regulatory elements, and characterized the in vivo enhancer activity of a large group of non-coding elements in the human genome that are conserved in human-pufferfish, Takifugu (Fugu) rubripes, or ultraconserved in human-mouse-rat. We tested 167 of these extremely conserved sequences in a transgenic mouse enhancer assay. Here we report that 45% of these sequences functioned reproducibly as tissue-specific enhancers of gene expression at embryonic day 11.5. While directing expression in a broad range of anatomical structures in the embryo, the majority of the 75 enhancers directed expression to various regions of the developing nervous system. We identified sequence signatures enriched in a subset of these elements that targeted forebrain expression, and used these features to rank all approximately 3,100 non-coding elements in the human genome that are conserved between human and Fugu. The testing of the top predictions in transgenic mice resulted in a threefold enrichment for sequences with forebrain enhancer activity. These data dramatically expand the catalogue of human gene enhancers that have been characterized in vivo, and illustrate the utility of such training sets for a variety of biological applications, including decoding the regulatory vocabulary of the human genome.     

Integration of telomere sequences with the draft human genome sequence

Telomeres are the ends of linear eukaryotic chromosomes. To ensure that no large stretches of uncharacterized DNA remain between the ends of the human working draft sequence and the ends of each chromosome, we would need to connect the sequences of the telomeres to the working draft sequence. But telomeres have an unusual DNA sequence composition and organization that makes them particularly difficult to isolate and analyse. Here we use specialized linear yeast artificial chromosome clones, each carrying a large telomere-terminal fragment of human DNA, to integrate most human telomeres with the working draft sequence. Subtelomeric sequence structure appears to vary widely, mainly as a result of large differences in subtelomeric repeat sequence abundance and organization at individual telomeres. Many subtelomeric regions appear to be gene-rich, matching both known and unknown expressed genes. This indicates that human subtelomeric regions are not simply buffers of nonfunctional 'junk DNA' next to the molecular telomere, but are instead functional parts of the expressed genome.     

Sequencing of Aspergillus nidulans and comparative analysis with A. fumigatus and A. oryzae

The aspergilli comprise a diverse group of filamentous fungi spanning over 200 million years of evolution. Here we report the genome sequence of the model organism Aspergillus nidulans, and a comparative study with Aspergillus fumigatus, a serious human pathogen, and Aspergillus oryzae, used in the production of sake, miso and soy sauce. Our analysis of genome structure provided a quantitative evaluation of forces driving long-term eukaryotic genome evolution. It also led to an experimentally validated model of mating-type locus evolution, suggesting the potential for sexual reproduction in A. fumigatus and A. oryzae. Our analysis of sequence conservation revealed over 5,000 non-coding regions actively conserved across all three species. Within these regions, we identified potential functional elements including a previously uncharacterized TPP riboswitch and motifs suggesting regulation in filamentous fungi by Puf family genes. We further obtained comparative and experimental evidence indicating widespread translational regulation by upstream open reading frames. These results enhance our understanding of these widely studied fungi as well as provide new insight into eukaryotic genome evolution and gene regulation.     

Genome sequence and comparative analysis of the model rodent malaria parasite Plasmodium yoelii yoelii

Species of malaria parasite that infect rodents have long been used as models for malaria disease research. Here we report the whole-genome shotgun sequence of one species, Plasmodium yoelii yoelii, and comparative studies with the genome of the human malaria parasite Plasmodium falciparum clone 3D7. A synteny map of 2,212 P. y. yoelii contiguous DNA sequences (contigs) aligned to 14 P. falciparum chromosomes reveals marked conservation of gene synteny within the body of each chromosome. Of about 5,300 P. falciparum genes, more than 3,300 P. y. yoelii orthologues of predominantly metabolic function were identified. Over 800 copies of a variant antigen gene located in subtelomeric regions were found. This is the first genome sequence of a model eukaryotic parasite, and it provides insight into the use of such systems in the modelling of Plasmodium biology and disease.     

DNA sequence and comparative analysis of chimpanzee chromosome 22

Human-chimpanzee comparative genome research is essential for narrowing down genetic changes involved in the acquisition of unique human features, such as highly developed cognitive functions, bipedalism or the use of complex language. Here, we report the high-quality DNA sequence of 33.3 megabases of chimpanzee chromosome 22. By comparing the whole sequence with the human counterpart, chromosome 21, we found that 1.44% of the chromosome consists of single-base substitutions in addition to nearly 68,000 insertions or deletions. These differences are sufficient to generate changes in most of the proteins. Indeed, 83% of the 231 coding sequences, including functionally important genes, show differences at the amino acid sequence level. Furthermore, we demonstrate different expansion of particular subfamilies of retrotransposons between the lineages, suggesting different impacts of retrotranspositions on human and chimpanzee evolution. The genomic changes after speciation and their biological consequences seem more complex than originally hypothesized.     

Conservation of Y-linked genes during human evolution revealed by comparative sequencing in chimpanzee

The human Y chromosome, transmitted clonally through males, contains far fewer genes than the sexually recombining autosome from which it evolved. The enormity of this evolutionary decline has led to predictions that the Y chromosome will be completely bereft of functional genes within ten million years. Although recent evidence of gene conversion within massive Y-linked palindromes runs counter to this hypothesis, most unique Y-linked genes are not situated in palindromes and have no gene conversion partners. The 'impending demise' hypothesis thus rests on understanding the degree of conservation of these genes. Here we find, by systematically comparing the DNA sequences of unique, Y-linked genes in chimpanzee and human, which diverged about six million years ago, evidence that in the human lineage, all such genes were conserved through purifying selection. In the chimpanzee lineage, by contrast, several genes have sustained inactivating mutations. Gene decay in the chimpanzee lineage might be a consequence of positive selection focused elsewhere on the Y chromosome and driven by sperm competition.     

Topological domains in mammalian genomes identified by analysis of chromatin interactions

The spatial organization of the genome is intimately linked to its biological function, yet our understanding of higher order genomic structure is coarse, fragmented and incomplete. In the nucleus of eukaryotic cells, interphase chromosomes occupy distinct chromosome territories, and numerous models have been proposed for how chromosomes fold within chromosome territories. These models, however, provide only few mechanistic details about the relationship between higher order chromatin structure and genome function. Recent advances in genomic technologies have led to rapid advances in the study of three-dimensional genome organization. In particular, Hi-C has been introduced as a method for identifying higher order chromatin interactions genome wide. Here we investigate the three-dimensional organization of the human and mouse genomes in embryonic stem cells and terminally differentiated cell types at unprecedented resolution. We identify large, megabase-sized local chromatin interaction domains, which we term 'topological domains', as a pervasive structural feature of the genome organization. These domains correlate with regions of the genome that constrain the spread of heterochromatin. The domains are stable across different cell types and highly conserved across species, indicating that topological domains are an inherent property of mammalian genomes. Finally, we find that the boundaries of topological domains are enriched for the insulator binding protein CTCF, housekeeping genes, transfer RNAs and short interspersed element (SINE) retrotransposons, indicating that these factors may have a role in establishing the topological domain structure of the genome.     

Functional genomic analysis of C. elegans chromosome I by systematic RNA interference

Complete genomic sequence is known for two multicellular eukaryotes, the nematode Caenorhabditis elegans and the fruit fly Drosophila melanogaster, and it will soon be known for humans. However, biological function has been assigned to only a small proportion of the predicted genes in any animal. Here we have used RNA-mediated interference (RNAi) to target nearly 90% of predicted genes on C. elegans chromosome I by feeding worms with bacteria that express double-stranded RNA. We have assigned function to 13.9% of the genes analysed, increasing the number of sequenced genes with known phenotypes on chromosome I from 70 to 378. Although most genes with sterile or embryonic lethal RNAi phenotypes are involved in basal cell metabolism, many genes giving post-embryonic phenotypes have conserved sequences but unknown function. In addition, conserved genes are significantly more likely to have an RNAi phenotype than are genes with no conservation. We have constructed a reusable library of bacterial clones that will permit unlimited RNAi screens in the future; this should help develop a more complete view of the relationships between the genome, gene function and the environment.     

Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution

We present here a draft genome sequence of the red jungle fowl, Gallus gallus. Because the chicken is a modern descendant of the dinosaurs and the first non-mammalian amniote to have its genome sequenced, the draft sequence of its genome--composed of approximately one billion base pairs of sequence and an estimated 20,000-23,000 genes--provides a new perspective on vertebrate genome evolution, while also improving the annotation of mammalian genomes. For example, the evolutionary distance between chicken and human provides high specificity in detecting functional elements, both non-coding and coding. Notably, many conserved non-coding sequences are far from genes and cannot be assigned to defined functional classes. In coding regions the evolutionary dynamics of protein domains and orthologous groups illustrate processes that distinguish the lineages leading to birds and mammals. The distinctive properties of avian microchromosomes, together with the inferred patterns of conserved synteny, provide additional insights into vertebrate chromosome architecture.