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: Though gene prediction will improve with every new protein that is discovered and through improvements in the current set of tools, there is a long way to go before the authors can decipher the precise exonic structure of every gene in the human genome using purely computational methodology.
Abstract: One of the first useful products from the human genome will be a set of predicted genes. Besides its intrinsic scientific interest, the accuracy and completeness of this data set is of considerable importance for human health and medicine. Though progress has been made on computational gene identification in terms of both methods and accuracy evaluation measures, most of the sequence sets in which the programs are tested are short genomic sequences, and there is concern that these accuracy measures may not extrapolate well to larger, more challenging data sets. Given the absence of experimentally verified large genomic data sets, we constructed a semiartificial test set comprising a number of short single-gene genomic sequences with randomly generated intergenic regions. This test set, which should still present an easier problem than real human genomic sequence, mimics the approximately 200kb long BACs being sequenced. In our experiments with these longer genomic sequences, the accuracy of GENSCAN, one of the most accurate ab initio gene prediction programs, dropped significantly, although its sensitivity remained high. Conversely, the accuracy of similarity-based programs, such as GENEWISE, PROCRUSTES, and BLASTX was not affected significantly by the presence of random intergenic sequence, but depended on the strength of the similarity to the protein homolog. As expected, the accuracy dropped if the models were built using more distant homologs, and we were able to quantitatively estimate this decline. However, the specificities of these techniques are still rather good even when the similarity is weak, which is a desirable characteristic for driving expensive follow-up experiments. Our experiments suggest that though gene prediction will improve with every new protein that is discovered and through improvements in the current set of tools, we still have a long way to go before we can decipher the precise exonic structure of every gene in the human genome using purely computational methodology.
204 citations
TL;DR: The breakpoints of the cosmopolitan inversion 2j of D. buzzatii were cloned and sequenced and it appears that the two pairs of target site duplications generated upon insertion were exchanged during the inversion event, and that the inversions arose by ectopic recombination between two copies of the transposon that were in opposite orientations.
Abstract: Although polymorphic inversions in Drosophila are very common, the origin of these chromosomal rearrangements is unclear. The breakpoints of the cosmopolitan inversion 2j of D. buzzatii were cloned and sequenced. Both breakpoints contain large insertions corresponding to a transposable element. It appears that the two pairs of target site duplications generated upon insertion were exchanged during the inversion event, and that the inversion arose by ectopic recombination between two copies of the transposon that were in opposite orientations. This is apparently the mechanism by which transposable elements generate natural inversions in Drosophila.
203 citations
TL;DR: The status of the map of genes whose chromosomal locations have been determined in both mouse and man is assessed, progress towards saturated maps is evaluated, and the manner in which this information is now being used in fields as diverse as medical and evolutionary genetics is illustrated.
200 citations
TL;DR: A set of four energy transfer dyes are synthesized and demonstrated their use in automated DNA sequencing and their values were reduced by 20-25% compared with the dichlororhodamine dyes alone.
Abstract: We have synthesized a set of four energy transfer dyes and demonstrated their use in automated DNA sequencing. The donor dyes are the 5- or 6-carboxy isomers of 4'-aminomethylfluorescein and the acceptor dyes are a novel set of four 4,7-dichloro-substituted rhodamine dyes which have narrower emission spectra than the standard, unsubstituted rhodamines. A rigid amino acid linker, 4-aminomethylbenzoic acid, was used to separate the dyes. The brightness of each dye in an automated sequencing instrument equipped with a dual line argon ion laser (488 and 514 nm excitation) was 2-2.5 times greater than the standard dye-primers with a 2 times reduction in multicomponent noise. The overall improvement in signal-to-noise was 4- to 5-fold. The utility of the new dye set was demonstrated by sequencing of a BAC DNA with an 80 kb insert. Measurement of the extinction coefficients and the relative quantum yields of the dichlororhodamine components of the energy transfer dyes showed their values were reduced by 20-25% compared with the dichlororhodamine dyes alone.
197 citations
TL;DR: The increased levels of Alu RNAs caused by cellular exposure to different stresses regulate protein synthesis by antagonizing PKR activation, which provides a functional role for mammalian short interspersed elements, prototypical junk DNA.
Abstract: Cell stress, viral infection, and translational inhibition increase the abundance of human Alu RNA, suggesting that the level of these transcripts is sensitive to the translational state of the cell To determine whether Alu RNA functions in translational homeostasis, we investigated its role in the regulation of double-stranded RNA-activated kinase PKR We found that overexpression of Alu RNA by cotransient transfection increased the expression of a reporter construct, which is consistent with an inhibitory effect on PKR Alu RNA formed stable, discrete complexes with PKR in vitro, bound PKR in vivo, and antagonized PKR activation both in vitro and in vivo Alu RNAs produced by either overexpression or exposure of cells to heat shock bound PKR, whereas transiently overexpressed Alu RNA antagonized virus-induced activation of PKR in vivo Cycloheximide treatment of cells decreased PKR activity, coincident with an increase in Alu RNA These observations suggest that the increased levels of Alu RNAs caused by cellular exposure to different stresses regulate protein synthesis by antagonizing PKR activation This provides a functional role for mammalian short interspersed elements, prototypical junk DNA
194 citations