Initial sequencing and analysis of the human genome.
TL;DR: The results of an international collaboration to produce and make freely available a draft sequence of the human genome are reported and an initial analysis is presented, describing some of the insights that can be gleaned from the sequence.
Abstract: The human genome holds an extraordinary trove of information about human development, physiology, medicine and evolution. Here we report the results of an international collaboration to produce and make freely available a draft sequence of the human genome. We also present an initial analysis of the data, describing some of the insights that can be gleaned from the sequence.
Citations
More filters
TL;DR: The definition and use of family-specific, manually curated gathering thresholds are explained and some of the features of domains of unknown function (also known as DUFs) are discussed, which constitute a rapidly growing class of families within Pfam.
Abstract: Pfam is a widely used database of protein families and domains. This article describes a set of major updates that we have implemented in the latest release (version 24.0). The most important change is that we now use HMMER3, the latest version of the popular profile hidden Markov model package. This software is approximately 100 times faster than HMMER2 and is more sensitive due to the routine use of the forward algorithm. The move to HMMER3 has necessitated numerous changes to Pfam that are described in detail. Pfam release 24.0 contains 11,912 families, of which a large number have been significantly updated during the past two years. Pfam is available via servers in the UK (http://pfam.sanger.ac.uk/), the USA (http://pfam.janelia.org/) and Sweden (http://pfam.sbc.su.se/).
14,075 citations
TL;DR: Comparative genomic analysis indicates vertebrate expansions of genes associated with neuronal function, with tissue-specific developmental regulation, and with the hemostasis and immune systems are indicated.
Abstract: A 2.91-billion base pair (bp) consensus sequence of the euchromatic portion of the human genome was generated by the whole-genome shotgun sequencing method. The 14.8-billion bp DNA sequence was generated over 9 months from 27,271,853 high-quality sequence reads (5.11-fold coverage of the genome) from both ends of plasmid clones made from the DNA of five individuals. Two assembly strategies-a whole-genome assembly and a regional chromosome assembly-were used, each combining sequence data from Celera and the publicly funded genome effort. The public data were shredded into 550-bp segments to create a 2.9-fold coverage of those genome regions that had been sequenced, without including biases inherent in the cloning and assembly procedure used by the publicly funded group. This brought the effective coverage in the assemblies to eightfold, reducing the number and size of gaps in the final assembly over what would be obtained with 5.11-fold coverage. The two assembly strategies yielded very similar results that largely agree with independent mapping data. The assemblies effectively cover the euchromatic regions of the human chromosomes. More than 90% of the genome is in scaffold assemblies of 100,000 bp or more, and 25% of the genome is in scaffolds of 10 million bp or larger. Analysis of the genome sequence revealed 26,588 protein-encoding transcripts for which there was strong corroborating evidence and an additional approximately 12,000 computationally derived genes with mouse matches or other weak supporting evidence. Although gene-dense clusters are obvious, almost half the genes are dispersed in low G+C sequence separated by large tracts of apparently noncoding sequence. Only 1.1% of the genome is spanned by exons, whereas 24% is in introns, with 75% of the genome being intergenic DNA. Duplications of segmental blocks, ranging in size up to chromosomal lengths, are abundant throughout the genome and reveal a complex evolutionary history. Comparative genomic analysis indicates vertebrate expansions of genes associated with neuronal function, with tissue-specific developmental regulation, and with the hemostasis and immune systems. DNA sequence comparisons between the consensus sequence and publicly funded genome data provided locations of 2.1 million single-nucleotide polymorphisms (SNPs). A random pair of human haploid genomes differed at a rate of 1 bp per 1250 on average, but there was marked heterogeneity in the level of polymorphism across the genome. Less than 1% of all SNPs resulted in variation in proteins, but the task of determining which SNPs have functional consequences remains an open challenge.
12,098 citations
TL;DR: In a four-genome analysis of 3' UTRs, approximately 13,000 regulatory relationships were detected above the estimate of false-positive predictions, thereby implicating as miRNA targets more than 5300 human genes, which represented 30% of the gene set.
11,624 citations
TL;DR: A mature web tool for rapid and reliable display of any requested portion of the genome at any scale, together with several dozen aligned annotation tracks, is provided at http://genome.ucsc.edu.
Abstract: As vertebrate genome sequences near completion and research refocuses to their analysis, the issue of effective genome annotation display becomes critical. A mature web tool for rapid and reliable display of any requested portion of the genome at any scale, together with several dozen aligned annotation tracks, is provided at http://genome.ucsc.edu. This browser displays assembly contigs and gaps, mRNA and expressed sequence tag alignments, multiple gene predictions, cross-species homologies, single nucleotide polymorphisms, sequence-tagged sites, radiation hybrid data, transposon repeats, and more as a stack of coregistered tracks. Text and sequence-based searches provide quick and precise access to any region of specific interest. Secondary links from individual features lead to sequence details and supplementary off-site databases. One-half of the annotation tracks are computed at the University of California, Santa Cruz from publicly available sequence data; collaborators worldwide provide the rest. Users can stably add their own custom tracks to the browser for educational or research purposes. The conceptual and technical framework of the browser, its underlying MYSQL database, and overall use are described. The web site currently serves over 50,000 pages per day to over 3000 different users.
9,605 citations
TL;DR: Velvet represents a new approach to assembly that can leverage very short reads in combination with read pairs to produce useful assemblies and is in close agreement with simulated results without read-pair information.
Abstract: We have developed a new set of algorithms, collectively called "Velvet," to manipulate de Bruijn graphs for genomic sequence assembly. A de Bruijn graph is a compact representation based on short words (k-mers) that is ideal for high coverage, very short read (25-50 bp) data sets. Applying Velvet to very short reads and paired-ends information only, one can produce contigs of significant length, up to 50-kb N50 length in simulations of prokaryotic data and 3-kb N50 on simulated mammalian BACs. When applied to real Solexa data sets without read pairs, Velvet generated contigs of approximately 8 kb in a prokaryote and 2 kb in a mammalian BAC, in close agreement with our simulated results without read-pair information. Velvet represents a new approach to assembly that can leverage very short reads in combination with read pairs to produce useful assemblies.
9,389 citations
References
More filters
TL;DR: This is the second half of a two-part, noncomprehensive review that presents an overview of different types of automation equipment used in genome sequencing.
Abstract: Automation for genomics has enabled a 43-fold increase in the total finished human genomic sequence in the world in the past four years. This is the second half of a two-part, noncomprehensive review that presents an overview of different types of automation equipment used in genome sequencing. The first part of the review, published in the previous issue, focused on automated procedures used to prepare DNA for sequencing or analysis. This second part of the review presents a look at available DNA sequencers and array technology and concludes with a look at future technologies. Alternate sequencing technologies including mass spectrometry, biochips, and single molecule analysis are included in this review.
133 citations
TL;DR: Position effect variegation (PEV), inducing allele-specific transcriptional repression of a gene located more centromeric, has been postulated as the underlying genetic mechanism of FSHD and detection of PEV in F SHD patients may require analysis of more specific cell types at particular developmental stages.
Abstract: Facioscapulohumeral muscular dystrophy (FSHD) is an autosomal dominant, neuromuscular disorder characterized by progressive weakness of muscles in the face, shoulder and upper arm. Deletion of integral copies of a 3.3 kb repeated unit from the subtelomeric region on chromosome 4q35 has been shown to be associated with FSHD. These repeated units which are apparently not transcribed, map very close to the 4q telomere and belong to a 3.3 kb repeat family dispersed over heterochromatic regions of the genome. Hence, position effect variegation (PEV), inducing allele-specific transcriptional repression of a gene located more centromeric, has been postulated as the underlying genetic mechanism of FSHD. This hypothesis has directed the search for the FSHD gene to the region centromeric to the repeated units. A CpG island was identified and found to be associated with the 5' untranslated region of a novel human gene, FRG1 (FSHD Region Gene 1). This evolutionary conserved gene is located about 100 kb proximal to the repeated units and belongs to a multigene family with FRG1 related sequences on multiple chromosomes. The mature chromosome 4 FRG1 transcript is 1042 bp in length and contains nine exons which encode a putative protein of 258 amino acid residues. Transcription of FRG1 was detected in several human tissues including placenta, lymphocytes, brain and muscle. To investigate a possible PEV mechanism, allele-specific FRG1 steady-state transcript levels were determined using RNA-based single-strand conformation polymorphism (SSCP) analysis. A polymorphic fragment contained within the first exon of FRG1 was amplified from reverse transcribed RNA from lymphocytes and muscle biopsies of patients and controls. No evidence for PEV mediated repression of allelic transcription was obtained in these tissues. However, detection of PEV in FSHD patients may require analysis of more specific cell types at particular developmental stages.
133 citations
TL;DR: The striking similarity between the evolutionary patterns of the EAR genes and those of the major histocompatibility complex, immunoglobulin, and T cell receptor genes stands in strong support of the hypothesis that host-defense and generation of diversity are among the primary physiological function of the rodent EARs.
Abstract: The mammalian RNase A superfamily comprises a diverse array of ribonucleolytic proteins that have a variety of biochemical activities and physiological functions Two rapidly evolving RNases of higher primates are of particular interest as they are major secretory proteins of eosinophilic leukocytes and have been found to possess anti-pathogen activities in vitro To understand how these RNases acquired this function during evolution and to develop animal models for the study of their functions in vivo, it is necessary to investigate these genes in many species Here, we report the sequences of 38 functional genes and 23 pseudogenes of the eosinophil-associated RNase (EAR) family from 5 rodent species Our phylogenetic analysis of these genes showed a clear pattern of evolution by a rapid birth-and-death process and gene sorting, a process characterized by rapid gene duplication and deactivation occurring differentially among lineages This process ultimately generates distinct or only partially overlapping inventories of the genes, even in closely related species Positive Darwinian selection also contributed to the diversification of these EAR genes The striking similarity between the evolutionary patterns of the EAR genes and those of the major histocompatibility complex, immunoglobulin, and T cell receptor genes stands in strong support of the hypothesis that host-defense and generation of diversity are among the primary physiological function of the rodent EARs The discovery of a large number of divergent EARs suggests the intriguing possibility that these proteins have been specifically tailored to fight against distinct rodent pathogens
132 citations
131 citations
TL;DR: A 9.7 kb segment encompassing exons 7-10 of the adrenoleukodystrophy locus of the X chromosome has duplicated to specific locations near the pericentromeric regions of human chromosomes 2p11,10p 11, 16p11 and 22q11, indicating that the autosomal ALD paralogs arose relatively recently during the course of higher primate evolution.
Abstract: A 9.7 kb segment encompassing exons 7-10 of the adrenoleukodystrophy (ALD) locus of the X chromosome has duplicated to specific locations near the pericentromeric regions of human chromosomes 2p11,10p11, 16p11 and 22q11. Comparative sequence analysis reveals 92-96% nucleotide identity, indicating that the autosomal ALD paralogs arose relatively recently during the course of higher primate evolution (5-10 million years ago). Analysis of sequences flanking the duplication region identifies the presence of an unusual GCTTTTTGC repeat which may be a sequence-specific integration site for the process of pericentromeric-directed transposition. The breakpoint sequence and phylogenetic analysis predict a two-step transposition model, in which a duplication from Xq28 to pericentromeric 2p11 occurred once, followed by a rapid distribution of a larger duplicon cassette among the pericentromeric regions. In addition to facilitating more effective mutation detection among ALD patients, these findings provide further insight into the molecular basis underlying a pericentromeric-directed mechanism for non-homologous interchromosomal exchange.
130 citations