Search

The genome of the bacterium Borrelia burgdorferi B31, the aetiologic agent of Lyme disease, contains a linear chromosome of 910,725 base pairs and at least 17 linear and circular plasmids with a combined size of more than 533,000 base pairs. The chromosome contains 853 genes encoding a basic set of proteins for DNA replication, transcription, translation, solute transport and energy metabolism, but, like Mycoplasma genitalium, it contains no genes for cellular biosynthetic reactions. Because B. burgdorferi and M. genitalium are distantly related eubacteria, we suggest that their limited metabolic capacities reflect convergent evolution by gene loss from more metabolically competent progenitors. Of 430 genes on 11 plasmids, most have no known biological function; 39% of plasmid genes are paralogues that form 47 gene families. The biological significance of the multiple plasmid-encoded genes is not clear, although they may be involved in antigenic variation or immune evasion.

Nerve growth factor-induced differentiation of adrenal chromaffin PC-12 cells to a neuronal phenotype involves alterations in gene expression and represents a model system to study neuronal differentiation. We have used the expressed-sequence-tag approach to identify approximately 600 differentially expressed mRNAs in untreated and nerve growth factor-treated PC-12 cells that encode proteins with diverse structural and biochemical functions. Many of these mRNAs encode proteins belonging to cellular pathways not previously known to be regulated by nerve growth factor. Comparative expressed-sequence-tag analysis provides a basis for surveying global changes in gene-expression patterns in response to biological signals at an unprecedented scale, is a powerful tool for identifying potential interactions between different cellular pathways, and allows the gene-expression profiles of individual genes belonging to a particular pathway to be followed.

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 ∼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.

Automated DNA sequencing is an extremely valuable technique which requires very high quality DNA templates to be carried out successfully. While it has been possible to readily produce large numbers of such templates from M13 or other single-stranded vectors for several years, the sequencing of double-stranded DNA templates using the ABI 373 DNA Sequencer has had a considerably lower success rate. We describe how the combination of a new fluorescent, dideoxy sequencing method, called cycle-sequencing, coupled with modifications to template isolation procedures based on Qiagen columns, makes fluorescent sequencing of double-stranded templates a reliable procedure. From a single five milliliter culture enough DNA can be isolated (up to 20 micrograms) to do 4-8 sequencing reactions, each of which yields 400-500 bases of high quality sequence data. These procedures make the routine use of double-stranded DNA templates a viable strategy in automated DNA sequencing projects.

The sequence of five non-contiguous genomic fragments encompassing 14.4 kilobases from the NF1 locus have been determined by fluorescence-based automated DNA sequence analysis. These fragments included one kilobase of the NF1 coding region, which resulted in the identification of the intron/exon boundaries of five exons. Based on these sequences, five new NF1 exon-PCR assays have been developed, that could be useful for detecting new NF1 mutations. The genomic sequences were analyzed for the presence of Alu repetitive elements and their classification is described. This analysis may provide some insight into the characterization of genetic rearrangements resulting in disruption of the NF1 gene.

Automated partial DNA sequencing was conducted on more than 600 randomly selected human brain complementary DNA (cDNA) clones to generate expressed sequence tags (ESTs). ESTs have applications in the discovery of new human genes, mapping of the human genome, and identification of coding regions in genomic sequences. Of the sequences generated, 337 represent new genes, including 48 with significant similarity to genes from other organisms, such as a yeast RNA polymerase II subunit; Drosophila kinesin, Notch, and Enhancer of split; and a murine tyrosine kinase receptor. Forty-six ESTs were mapped to chromosomes after amplification by the polymerase chain reaction. This fast approach to cDNA characterization will facilitate the tagging of most human genes in a few years at a fraction of the cost of complete genomic sequencing, provide new genetic markers, and serve as a resource in diverse biological research fields.

The high degree of similarity between the mouse and human genomes is demonstrated through analysis of the sequence of mouse chromosome 16 (Mmu 16), which was obtained as part of a whole-genome shotgun assembly of the mouse genome. The mouse genome is about 10% smaller than the human genome, owing to a lower repetitive DNA content. Comparison of the structure and protein-coding potential of Mmu 16 with that of the homologous segments of the human genome identifies regions of conserved synteny with human chromosomes (Hsa) 3, 8, 12, 16, 21, and 22. Gene content and order are highly conserved between Mmu 16 and the syntenic blocks of the human genome. Of the 731 predicted genes on Mmu 16, 509 align with orthologs on the corresponding portions of the human genome, 44 are likely paralogous to these genes, and 164 genes have homologs elsewhere in the human genome; there are 14 genes for which we could find no human counterpart.

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 ∼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 ∼13,600 genes, somewhat fewer than the smaller Caenorhabditis elegans genome, but with comparable functional diversity.

Archaeoglobus fulgidus is the first sulphur-metabolizing organism to have its genome sequence determined. Its genome of 2,178,400 base pairs contains 2,436 open reading frames (ORFs). The information processing systems and the biosynthetic pathways for essential components (nucleotides, amino acids and cofactors) have extensive correlation with their counterparts in the archaeon Methanococcus jannaschii. The genomes of these two Archaea indicate dramatic differences in the way these organisms sense their environment, perform regulatory and transport functions, and gain energy. In contrast to M. jannaschii, A. fulgidus has fewer restriction-modification systems, and none of its genes appears to contain inteins. A quarter (651 ORFs) of the A. fulgidus genome encodes functionally uncharacterized yet conserved proteins, two-thirds of which are shared with M. jannaschii (428 ORFs). Another quarter of the genome encodes new proteins indicating substantial archaeal gene diversity.

The complete 1.66-megabase pair genome sequence of an autotrophic archaeon, Methanococcus jannaschii, and its 58- and 16-kilobase pair extrachromosomal elements have been determined by whole-genome random sequencing. A total of 1738 predicted protein-coding genes were identified; however, only a minority of these (38 percent) could be assigned a putative cellular role with high confidence. Although the majority of genes related to energy production, cell division, and metabolism in M. jannaschii are most similar to those found in Bacteria, most of the genes involved in transcription, translation, and replication in M. jannaschii are more similar to those found in Eukaryotes.

Large-scale DNA sequencing and analysis is providing important molecular information on complex biological systems.

A total of 116,118 basepairs (bp) derived from three cosmids spanning the ERCC1 locus of human chromosome 19q13.3 have been sequenced with automated fluorescence-based sequencers and analysed by polymerase chain reaction amplification and computer methods. The assembled sequence forms two contigs totalling 105,831 bp, which contain a human fosB proto-oncogene, a gene encoding a protein phosphatase, two genes of unknown function and the previously-characterized ERCC1 DNA repair gene. This light band region has a high average density of 1.4 Alu repeats per kilobase. Human chromosome light bands could therefore contain up to 75,000 genes and 1.5 million Alu repeats.

Genomic DNA that encodes the beta 1 subunit of the human gamma-aminobutyric acidA (GABAA) receptor was cloned and mapped. Exons and flanking introns (greater than 14 kb) were sequenced to determine the structural organization of the gene. The gene was localized on human chromosome 4, in bands p12-13. The beta 1 subunit is encoded by a relatively large gene (greater than 65 kb) on nine exons. In contrast to other conserved regions of the subunit polypeptide, the proposed channel-forming domain (M2) is derived from more than one exon. The organization of exons was compared with that of the genes that code for subunits of nicotinic acetylcholine receptors. There is no evidence for conservation of gene structure between these two members of the proposed gene superfamily. However, intron-exon junctions were found to be conserved precisely between subtypes of GABAA receptor subunits.

Two cDNA clones, lambda RHM-MF and lambda RHB-DAR, encoding the muscarinic cholinergic receptor and the beta-adrenergic receptor, respectively, have been isolated from a rat heart cDNA library. The cDNA clones were characterized by restriction mapping and automated DNA sequence analysis utilizing fluorescent dye primers. The rat heart muscarinic receptor consists of 466 amino acids and has a calculated molecular weight of 51,543. The rat heart beta-adrenergic receptor consists of 418 amino acids and has a calculated molecular weight of 46,890. The two cardiac receptors have substantial amino acid homology (27.2% identity, 50.6% with favored substitutions). The rat cardiac beta receptor has 88.0% homology (92.5% with favored substitutions) with the human brain beta receptor and the rat cardiac muscarinic receptor has 94.6% homology (97.6% with favored substitutions) with the porcine cardiac muscarinic receptor. The muscarinic cholinergic and beta-adrenergic receptors appear to be as conserved as hemoglobin and cytochrome c but less conserved than histones and are clearly members of a multigene family. These data support our hypothesis, based upon biochemical and immunological evidence, that suggests considerable structural homology and evolutionary conservation between adrenergic and muscarinic cholinergic receptors. To our knowledge, this is the first report utilizing automated DNA sequence analysis to determine the structure of a gene.

Two cDNA clones (3.7 kb and 4.8 kb) encoding a Drosophila muscarinic acetylcholine receptor were isolated from a Drosophila head cDNA library and characterized by automated DNA sequence analysis. The Drosophila muscarinic receptor contains 788 amino acids with a calculated Mr of 84 807 and displays greater than 60% homology with mammalian muscarinic receptors. The muscarinic receptor maps to the tip of the right arm of the second chromosome of the Drosophila genome.

We have extensively characterized the structure and evolution of the adrenergic and muscarinic cholinergic receptors. This analysis has included studies with monoclonal antibodies, protein purification, target size analysis, gene cloning and sequencing and gene expression (1–7). Based upon the results of biochemical and immunological studies, we proposed that adrenergic and muscarinic cholinergic receptors are highly conserved proteins. This hypothesis has been confirmed by gene cloning experiments which have provided the primary structures of a number of adrenergic and muscarinic cholinergic receptors.

Two cDNA clones, lambda-CLFV-108 and lambda-CLFV-119, encoding for the beta-adrenergic receptor, have been isolated from a human brain stem cDNA library. One human genomic clone, LCV-517 (20 kb), was characterized by restriction mapping and partial sequencing. The human brain beta-receptor consists of 413 amino acids with a calculated Mr of 46480. The gene contains three potential glucocorticoid receptor-binding sites. The beta-receptor expressed in human brain was homology with rodent (88%) and avian (52%) beta-receptors and with porcine muscarinic cholinergic receptors (31%), supporting our proposal [(1984) Proc. Natl. Acad. Sci. USA 81, 272 276] that adrenergic and muscarinic cholinergic receptors are structurally related. This represents the first cloning of a neurotransmitter receptor gene from human brain.

An approach for genome analysis based on sequencing and assembly of unselected pieces of DNA from the whole chromosome has been applied to obtain the complete nucleotide sequence (1,830,137 base pairs) of the genome from the bacterium Haemophilus influenzae Rd. This approach eliminates the need for initial mapping efforts and is therefore applicable to the vast array of microbial species for which genome maps are unavailable. The H. influenzae Rd genome sequence (Genome Sequence DataBase accession number L42023) represents the only complete genome sequence from a free-living organism.