The sequence of the human genome

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.

The genome sequence of Drosophila melanogaster

The fly Drosophila melanogaster is one of the most intensively studied organisms in biology and serves as a model system for the investigation of many developmental and cellular processes common to higher eukaryotes, including humans. We have determined the nucleotide sequence of nearly all of the approximately 120-megabase euchromatic portion of the Drosophila genome using a whole-genome shotgun sequencing strategy supported by extensive clone-based sequence and a high-quality bacterial artificial chromosome physical map. Efforts are under way to close the remaining gaps; however, the sequence is of sufficient accuracy and contiguity to be declared substantially complete and to support an initial analysis of genome structure and preliminary gene annotation and interpretation. The genome encodes approximately 13,600 genes, somewhat fewer than the smaller Caenorhabditis elegans genome, but with comparable functional diversity.

Comparative genomics of the eukaryotes

A comparative analysis of the genomes of Drosophila melanogaster, Caenorhabditis elegans, and Saccharomyces cerevisiae-and the proteins they are predicted to encode-was undertaken in the context of cellular, developmental, and evolutionary processes. The nonredundant protein sets of flies and worms are similar in size and are only twice that of yeast, but different gene families are expanded in each genome, and the multidomain proteins and signaling pathways of the fly and worm are far more complex than those of yeast. The fly has orthologs to 177 of the 289 human disease genes examined and provides the foundation for rapid analysis of some of the basic processes involved in human disease.

A general approach to single-nucleotide polymorphism discovery

Single-nucleotide polymorphisms (SNPs) are the most abundant form of human genetic variation and a resource for mapping complex genetic traits. The large volume of data produced by high-throughput sequencing projects is a rich and largely untapped source of SNPs (refs 2, 3, 4, 5). We present here a unified approach to the discovery of variations in genetic sequence data of arbitrary DNA sources. We propose to use the rapidly emerging genomic sequence as a template on which to layer often unmapped, fragmentary sequence data and to use base quality values to discern true allelic variations from sequencing errors. By taking advantage of the genomic sequence we are able to use simpler yet more accurate methods for sequence organization: fragment clustering, paralogue identification and multiple alignment. We analyse these sequences with a novel, Bayesian inference engine, POLYBAYES, to calculate the probability that a given site is polymorphic. Rigorous treatment of base quality permits completely automated evaluation of the full length of all sequences, without limitations on alignment depth. We demonstrate this approach by accurate SNP predictions in human ESTs aligned to finished and working-draft quality genomic sequences, a data set representative of the typical challenges of sequence-based SNP discovery.

A BMP homolog acts as a dose-dependent regulator of body size and male tail patterning in Caenorhabditis elegans

We cloned the dbl-1 gene, a C. elegans homolog of Drosophila decapentaplegic and vertebrate BMP genes. Loss-of-function mutations in dbl-1 cause markedly reduced body size and defective male copulatory structures. Conversely, dbl-1 overexpression causes markedly increased body size and partly complementary male tail phenotypes, indicating that DBL-1 acts as a dose-dependent regulator of these processes. Evidence from genetic interactions indicates that these effects are mediated by a Smad signaling pathway, for which DBL-1 is a previously unidentified ligand. Our study of the dbl-1 expression pattern suggests a role for neuronal cells in global size regulation as well as male tail patterning.

Identification and characterization of a p53 homologue in Drosophila melanogaster

The tumor suppressor gene p53 in mammalian cells plays a critical role in safeguarding the integrity of genome. It functions as a sequence-specific transcription factor. Upon activation by a variety of cellular stresses, p53 transactivates downstream target genes, through which it regulates cell cycle and apoptosis. However, little is known about p53 in invertebrates. Here we report the identification and characterization of a Drosophila p53 homologue gene, dp53. dp53 encodes a 385-amino acid protein with significant homology to human p53 (hp53) in the region of the DNA-binding domain, and to a lesser extent the tetramerization domain. Purified dp53 DNA-binding domain protein was shown to bind to the consensus hp53-binding site by gel mobility analysis. In transient transfection assays, expression of dp53 in Schneider cells transcriptionally activated promoters that contained consensus hp53-responsive elements. Moreover, a mutant dp53 (Arg-155 to His-155), like its hp53 counterpart mutant, exerted a dominant-negative effect on transactivation. Ectopic expression of dp53 in Drosophila eye disk caused cell death and led to a rough eye phenotype. dp53 is expressed throughout the development of Drosophila with highest expression levels in early embryogenesis, which has a maternal component. Consistent with this, dp53 RNA levels were high in the nurse cells of the ovary. It appears that p53 is structurally and functionally conserved from flies to mammals. Drosophila will provide a useful genetic system to the further study of the p53 network.

Age-related mutations and chronic myelomonocytic leukemia

Chronic myelomonocytic leukemia (CMML) is a hematologic malignancy nearly confined to the elderly. Previous studies to determine incidence and prognostic significance of somatic mutations in CMML have relied on candidate gene sequencing, although an unbiased mutational search has not been conducted. As many of the genes commonly mutated in CMML were recently associated with age-related clonal hematopoiesis (ARCH) and aged hematopoiesis is characterized by a myelomonocytic differentiation bias, we hypothesized that CMML and aged hematopoiesis may be closely related. We initially established the somatic mutation landscape of CMML by whole exome sequencing followed by gene-targeted validation. Genes mutated in ⩾10% of patients were SRSF2, TET2, ASXL1, RUNX1, SETBP1, KRAS, EZH2, CBL and NRAS, as well as the novel CMML genes FAT4, ARIH1, DNAH2 and CSMD1. Most CMML patients (71%) had mutations in ⩾2 ARCH genes and 52% had ⩾7 mutations overall. Higher mutation burden was associated with shorter survival. Age-adjusted population incidence and reported ARCH mutation rates are consistent with a model in which clinical CMML ensues when a sufficient number of stochastically acquired age-related mutations has accumulated, suggesting that CMML represents the leukemic conversion of the myelomonocytic-lineage-biased aged hematopoietic system.

Trimethylpsoralen induces small deletion mutations in Caenorhabditis elegans

To examine the mutagenic spectrum of 4,5',8-trimethylpsoralen (TMP) in Caenorhabditis elegans, we isolated mutations in the unc-22 and pal-1 genes following TMP mutagenesis and analyzed them for restriction fragment length polymorphisms by Southern blot. Eleven of 21 unc-22 mutations exhibited restriction fragment length polymorphisms, 8 of which were deletions of between 0.10 and 15 kb in length. Both of two pal-1 mutations were also small deletions within this size range. Comparison of our results with previous studies on mutagenesis by gamma-rays and x-rays suggests that the mutagenic spectrum of TMP may be similar. TMP should be useful in generating mutations that cause complete loss of function of single genes and that are likely to result in allele-specific DNA polymorphisms.

Identification of key concepts in biomedical literature using a modified Markov heuristic

MOTIVATION

The recent explosion of interest in mining the biomedical literature for associations between defined entities such as genes, diseases and drugs has made apparent the need for robust methods of identifying occurrences of these entities in biomedical text. Such concept-based indexing is strongly dependent on the availability of a comprehensive ontology or lexicon of biomedical terms. However, such ontologies are very difficult and expensive to construct, and often require extensive manual curation to render them suitable for use by automatic indexing programs. Furthermore, the use of statistically salient noun phrases as surrogates for curated terminology is not without difficulties, due to the lack of high-quality part-of-speech taggers specific to medical nomenclature.

RESULTS

We describe a method of improving the quality of automatically extracted noun phrases by employing prior knowledge during the HMM training procedure for the tagger. This enhancement, when combined with appropriate training data, can greatly improve the quality and relevance of the extracted phrases, thereby enabling greater accuracy in downstream literature mining tasks.