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: Genome compaction and synteny conservation between human and Fugu were studied using data from public databases to gain insights to the extent and mechanisms of divergence.
Abstract: Background: Knowledge of the amount of gene order and synteny conservation between two species gives insights to the extent and mechanisms of divergence. The vertebrate Fugu rubripes (pufferfish) has a small genome with little repetitive sequence which makes it attractive as a model genome. Genome compaction and synteny conservation between human and Fugu were studied using data from public databases.
Methods: Intron length and map positions of human and Fugu orthologues were compared to analyse relative genome compaction and synteny conservation respectivley. The divergence of these two genomes by genome rearrangement was simulated and the results were compared to the real data.
Results: Analysis of 199 introns in 22 orthologous genes showed an eight-fold average size reduction in Fugu, consistent with the ratio of total genome sizes. There was no consistent pattern relating the size reduction in individual introns or genes to gene base composition in either species. For genes that are neighbours in Fugu (genes from the same cosmid or GenBank entry), 40–50% have conserved synteny with a human chromosome. This figure may be underestimated by as much as two-fold, due to problems caused by incomplete human genome sequence data and the existence of dispersed gene families. Some genes that are neighbours in Fugu have human orthologues that are several megabases and tens of genes apart. This is probably caused by small inversions or other intrachromosomal rearrangements.
Conclusions: Comparison of observed data to computer simulations suggests that 4000–16 000 chromosomal rearrangements have occured since Fugu and human shared a common ancestor, implying a faster rate of rearrangement than seen in human/mouse comparisons.
85 citations
TL;DR: If 5'-flanking sequences, including a GC-motif in the -40 region, are removed from the genes, transcription is reduced with a factor 10 or more, suggesting that sequences upstream of the coding region are important for the level of transcription.
Abstract: The human 5S rRNA genes are found in clusters of tandem repeated units. We have cloned and partially characterized six restriction fragments from two clusters of 2.3 kb and 1.6 kb repeats, respectively. Four fragments from the cluster of 2.3 kb repeats contain a 5S rRNA gene and one fragment contains a gene variant with an additional nucleotide in the internal control region. A fragment from the 1.6 kb cluster contains a gene and is highly homologous to the 2.3 kb repeats, except for a large deletion in the 3'-flanking region starting 12 bp downstream of the gene. The number of genes and closely related gene variants is found to be 300-400 per haploid human genome. 100-150 of these are found in 2.3 kb repeats and 5-10 are found in 1.6 kb repeats. The total number of 5S rRNA sequences, including pseudogenes, is 1700-2000 per haploid genome. The genes and the gene variant are transcribed equally efficient in a HeLa cell extract. If 5'-flanking sequences, including a GC-motif in the -40 region, are removed from the genes, transcription is reduced with a factor 10 or more, suggesting that sequences upstream of the coding region are important for the level of transcription.
85 citations
TL;DR: The immunoglobulin kappa locus and its immediate surroundings, which are described in the present report, comprise 3 Mb of DNA, i.e., 0.1% or one per mill of the 3000 Mb of the human genome.
84 citations
TL;DR: Replication asynchrony across the GABRA5 duplication is found, suggesting for the first time that the imprinted part of chromosome 15q extends proximal of the region commonly deleted in Angelman and Prader-Willi syndromes.
Abstract: We report the identification of a partial duplication of GABRA5 , a gene within the imprinted 15q11-q13 region. The duplicated locus maps to the pericentromeric region of 15q, proximal to the large deletions associated with Angelman and Prader-Willi syndromes. We also observed variation in the number of copies of this locus in different individuals, indicating that the duplication is part of a variable repeat. Investigation of the duplication in individuals with a normal karyotype revealed between one and four copies of the repeat on each chromosome 15, whereas from eight to 20 copies were found in individuals possessing a cytogenetically detectable elongation of the 15q region. The variable region is roughly 1 Mb in size and contains two other non-processed duplications, the immunoglobulin heavy chain (IgH) D segment gene and the neurofibromatosis type 1 (NF1) gene. One unit of the pericentromeric repeat is thus composed of duplications of genes from different chromosomal regions. Moreover, we have found replication asynchrony across the GABRA5 duplication, suggesting for the first time that the imprinted part of chromosome 15q extends proximal of the region commonly deleted in Angelman and Prader-Willi syndromes.
84 citations
TL;DR: FIP was found to be present in goldfish brain and to be formed as a labeled product in vivo and the results suggest that FIP may be formed at extracellular and postsynaptic loci where multiple associating inputs interact at CA1.
Abstract: Ependymin, a glycoprotein of the brain ECF, has been implicated in the neurochemistry of memory and neuronal regeneration. Three behavioral experiments (swimming with a float, avoidance conditioning, and classical conditioning) in the goldfish and one in the mouse (T-maze learning) indicate that ependymin has a role in the synaptic changes that take place in the consolidation step of memory formation and the activity-dependent phase of sharpening of goldfish retinotectal connections during neuronal regeneration. The ECF concentration of the protein was found to decrease after the goldfish learned to associate a light stimulus (CS) with the subsequent arrival of a shock (US): paired CS-US gave changes whereas an unpaired presentation of CS-US gave no changes relative to the unstimulated controls. Ependymin is present in ECF as a mixture of three disulfide-linked dimers of two acidic (alpha and beta) polypeptide chains (37 kDa and 31 kDa). Upon removal of its N-linked glycan fragment by N-glycosidase F, the beta chain yields gamma-ependymin (26 kDa). Determinations of the amino acid sequence of gamma-ependymin indicate that it is a unique protein with no long sequence homologies to any known polypeptide. There are, however, small segments (5-7 amino acids long) with homologies to fibronectin, laminin, and tubulin. Ependymin has the capacity to polymerize into FIP (after activation by phosphorylation) in response to events that deplete ECF calcium. FIP is insoluble in 2% SDS in 6 M urea, 10 mM Ca2+Ac2, 100% acetic acid, chloroform/methanol (2/1), saturated KCNS, and even 100% trifluoroacetic acid. FIP was found to be present in goldfish brain and to be formed as a labeled product in vivo. Ependymin's FIP-forming property was used to propose a molecular hypothesis for generating synaptic changes in response to local extracellular depletions of calcium at sites of "associating inputs." The model assumes that, following NMDA receptor stimulation, the translocated PKC that is generated activates extracellular ependymin by converting it to its phosphorylated form using presynaptically released ATP. The hypothesis was tested in studies of LTP of rat hippocampal slices at CA1. After LTP, new sites that stained with antisera to ependymin, visible at 100x, were obtained in its potentiated radiatum in the CA1 region but not in the unpotentiated CA3. Electron microscopic studies showed that the horseradish peroxidase reaction product obtained was localized at synaptic clefts and postsynaptic regions. The results suggest that FIP may be formed at extracellular and postsynaptic loci where multiple associating inputs interact at CA1.
83 citations