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.

Genomic sequence of a Lyme disease spirochaete, Borrelia burgdorferi.

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.

The complete genome sequence of the hyperthermophilic, sulphate-reducing archaeon Archaeoglobus fulgidus.

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 genome sequence of the gastric pathogen Helicobacter pylori.

Helicobacter pylori, strain 26695, has a circular genome of 1,667,867 base pairs and 1,590 predicted coding sequences. Sequence analysis indicates that H. pylori has well-developed systems for motility, for scavenging iron, and for DNA restriction and modification. Many putative adhesins, lipoproteins and other outer membrane proteins were identified, underscoring the potential complexity of host-pathogen interaction. Based on the large number of sequence-related genes encoding outer membrane proteins and the presence of homopolymeric tracts and dinucleotide repeats in coding sequences, H. pylori, like several other mucosal pathogens, probably uses recombination and slipped-strand mispairing within repeats as mechanisms for antigenic variation and adaptive evolution. Consistent with its restricted niche, H. pylori has a few regulatory networks, and a limited metabolic repertoire and biosynthetic capacity. Its survival in acid conditions depends, in part, on its ability to establish a positive inside-membrane potential in low pH.

Complete genome sequence of the methanogenic archaeon, Methanococcus jannaschii.

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.

The minimal gene complement of Mycoplasma genitalium.

The complete nucleotide sequence (580,070 base pairs) of the Mycoplasma genitalium genome, the smallest known genome of any free-living organism, has been determined by whole-genome random sequencing and assembly. A total of only 470 predicted coding regions were identified that include genes required for DNA replication, transcription and translation, DNA repair, cellular transport, and energy metabolism. Comparison of this genome to that of Haemophilus influenzae suggests that differences in genome content are reflected as profound differences in physiology and metabolic capacity between these two organisms.

Initial assessment of human gene diversity and expression patterns based upon 83 million nucleotides of cDNA sequence.

In an effort to identify new genes and analyse their expression patterns, 174,472 partial complementary DNA sequences (expressed sequence tags (ESTs)), totalling more than 52 million nucleotides of human DNA sequence, have been generated from 300 cDNA libraries constructed from 37 distinct organs and tissues. These ESTs have been combined with an additional 118,406 ESTs from the database dbEST, for a total of 83 million nucleotides, and treated as a shotgun sequence assembly project. The assembly process yielded 29,599 distinct tentative human consensus (THC) sequences and 58,384 non-overlapping ESTs. Of these 87,983 distinct sequences, 10,214 further characterize previously known genes based on statistically significant similarity to sequences in the available databases; the remainder identify previously unknown genes. Thirty tissues were sampled by over 1,000 ESTs each; only eight genes were matched by ESTs from all 30 tissues, and 227 genes were represented in 20 or more of the tissues sampled with more than 1,000 ESTs. Approximately 40% of identified human genes appear to be associated with basic energy metabolism, cell structure, homeostasis and cell division, 22% with RNA and protein synthesis and processing, and 12% with cell signalling and communication.

Comparative expressed-sequence-tag analysis of differential gene expression profiles in PC-12 cells before and after nerve growth factor treatment.

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.

Expressed genes, Alu repeats and polymorphisms in cosmids sequenced from chromosome 4p16.3.

The sequences of three cosmids (90 kilobases) from the Huntington's disease region in chromosome 4p16.3 have been determined. A 30,837 base overlap of DNA sequenced from two individuals was found to contain 72 DNA sequence polymorphisms, an average of 2.3 polymorphisms per kilobase (kb). The assembled 58 kb contig contains 62 Alu repeats, and eleven predicted exons representing at least three expressed genes that encode previously unidentified proteins. Each of these genes is associated with a CpG island. The structure of one of the new genes, hda1-1, has been determined by characterizing cDNAs from a placental library. This gene is expressed in a variety of tissues and may encode a novel housekeeping gene.

Automated DNA sequencing and analysis of 106 kilobases from human chromosome 19q13.3.

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.

Rapid and reliable fluorescent cycle sequencing of double-stranded templates.

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.

Sequencing and analysis of genomic fragments from the NF1 locus.

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.

An adipose tissue-specific beta-adrenergic receptor. Molecular cloning and down-regulation in obesity.

Clones encoding an atypical beta-adrenergic receptor were isolated from a rat brown adipose tissue cDNA library. This receptor expressed in Chinese hamster ovary (CHO) cells displays a low affinity for beta-adrenergic antagonists and a high affinity for BRL 37344, an agonist that selectively stimulates lipolysis in adipose tissue. The rank order of potency for agonist-mediated increases in intracellular cAMP in transfected cells correlates with that for agonist-mediated stimulation of lipolysis in brown adipocytes. Northern blot analysis demonstrates that this receptor subtype is expressed only in brown and white adipose tissue where it represents the predominant beta-receptor subtype. The amount of atypical beta-adrenergic receptor present in adipose tissue of obese (fa/fa) Zucker rats is reduced by up to 71% as compared with lean (Fa/Fa) control animals. These findings suggest that a change in the expression of this beta-adrenergic receptor subtype may play a role in obesity.

Isolation, characterization, and localization of human genomic DNA encoding the beta 1 subunit of the GABAA receptor (GABRB1).

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.

Complementary DNA sequencing: expressed sequence tags and human genome project.

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.

Cloning, localization, and permanent expression of a Drosophila octopamine receptor.

A cDNA for a member of the G protein-coupled receptor family was isolated from Drosophila using a probe derived from a human beta 2-adrenergic receptor cDNA. This Drosophila receptor gene is localized at 99A10-B1 on the right arm of chromosome 3 and is preferentially expressed in Drosophila heads. The insect octopamine receptor has been permanently expressed in mammalian cells, where it mediates the attenuation of adenylate cyclase activity and exhibits a pharmacological profile consistent with an octopamine type 1 receptor. Sequence and pharmacological comparisons indicate that the octopamine receptor is unique but closely related to mammalian adrenergic receptors, perhaps as an evolutionary precursor.

Site-directed mutagenesis of m1 muscarinic acetylcholine receptors: conserved aspartic acids play important roles in receptor function.

Muscarinic acetylcholine receptors contain a region encompassing the second and third transmembrane domains that is rich in conserved aspartic acid residues. To investigate the role of four conserved aspartic acids at positions 71, 99, 105, and 122 in muscarinic receptor function, point mutations in the rat m1 muscarinic receptor gene were made that converted each Asp to Asn, and wild type or mutant genes were stably expressed in Chinese hamster ovary cells that normally lack muscarinic receptors. Substitution of Asp71 or Asp122 with Asn produced mutant receptors that displayed high affinity for carbachol but decreased efficacy and potency, respectively, in agonist-induced activation of phosphoinositide hydrolysis, suggesting that these residues may mediate receptor-GTP binding protein interactions. Substitution of Asp99 or Asp105 with Asn produced marked decreases in ligand binding affinities and/or covalent incorporation of [3H] propylbenzilylcholine mustard, suggesting that these residues may be involved in receptor-ligand interactions.

Cloning, sequence analysis and chromosome localization of a Drosophila muscarinic acetylcholine receptor.

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.