Whole-genome comparison of Mycobacterium tuberculosis clinical and laboratory strains.

Virulence and immunity are poorly understood in Mycobacterium tuberculosis. We sequenced the complete genome of the M. tuberculosis clinical strain CDC1551 and performed a whole-genome comparison with the laboratory strain H37Rv in order to identify polymorphic sequences with potential relevance to disease pathogenesis, immunity, and evolution. We found large-sequence and single-nucleotide polymorphisms in numerous genes. Polymorphic loci included a phospholipase C, a membrane lipoprotein, members of an adenylate cyclase gene family, and members of the PE/PPE gene family, some of which have been implicated in virulence or the host immune response. Several gene families, including the PE/PPE gene family, also had significantly higher synonymous and nonsynonymous substitution frequencies compared to the genome as a whole. We tested a large sample of M. tuberculosis clinical isolates for a subset of the large-sequence and single-nucleotide polymorphisms and found widespread genetic variability at many of these loci. We performed phylogenetic and epidemiological analysis to investigate the evolutionary relationships among isolates and the origins of specific polymorphic loci. A number of these polymorphisms appear to have occurred multiple times as independent events, suggesting that these changes may be under selective pressure. Together, these results demonstrate that polymorphisms among M. tuberculosis strains are more extensive than initially anticipated, and genetic variation may have an important role in disease pathogenesis and immunity.

Comparative genomics and understanding of microbial biology.

The sequences of close to 30 microbial genomes have been completed during the past 5 years, and the sequences of more than 100 genomes should be completed in the next 2 to 4 years. Soon, completed microbial genome sequences will represent a collection of >200,000 predicted coding sequences. While analysis of a single genome provides tremendous biological insights on any given organism, comparative analysis of multiple genomes provides substantially more information on the physiology and evolution of microbial species and expands our ability to better assign putative function to predicted coding sequences.

Characterization of Porphyromonas gingivalis insertion sequence-like element ISPg5.

Porphyromonas gingivalis, a black-pigmented, gram-negative anaerobe, is found in periodontitis lesions, and its presence in subgingival plaque significantly increases the risk for periodontitis. In contrast to many bacterial pathogens, P. gingivalis strains display considerable variability, which is likely due to genetic exchange and intragenomic changes. To explore the latter possibility, we have studied the occurrence of insertion sequence (IS)-like elements in P. gingivalis W83 by utilizing a convenient and rapid method of capturing IS-like sequences and through analysis of the genome sequence of P. gingivalis strain W83. We adapted the method of Matsutani et al. (S. Matsutani, H. Ohtsubo, Y. Maeda, and E. Ohtsubo, J. Mol. Biol. 196:445-455, 1987) to isolate and clone rapidly annealing DNA sequences characteristic of repetitive regions within a genome. We show that in P. gingivalis strain W83, such sequences include (i) nucleotide sequence with homology to tRNA genes, (ii) a previously described IS element, and (iii) a novel IS-like element. Analysis of the P. gingivalis genome sequence for the distribution of the least used tetranucleotide, CTAG, identified regions in many of the initial 218 contigs which contained CTAG clusters. Examination of these CTAG clusters led to the discovery of 11 copies of the same novel IS-like element identified by the repeated sequence capture method of Matsutani et al. This new 1,512-bp IS-like element, designated ISPg5, has features of the IS3 family of IS elements. When a recombinant plasmid containing much of ISPg5 was used in Southern analysis of several P. gingivalis strains, including clinical isolates, diversity among strains was apparent. This suggests that ISPg5 and other IS elements may contribute to strain diversity and can be used for strain fingerprinting.

Genome data: what do we learn?

Genome sequence information has continued to accumulate at a spectacular pace during the past year. Details of the sequence and gene content of human chromosome 22 were published. The sequencing and annotation of the first two Arabidopsis thaliana chromosomes was completed. The sequence of chromosome 3 from Plasmodium falciparum, the second sequenced malaria chromosome, was reported, as was that of chromosome 1 from Leishmania major. The complete genomic sequences of five microbes were reported. Approaches to using data from completely sequenced microbial genomes in phylogenetic studies are being explored, as is the application of microarrays to whole genome expression analysis.

Complete genome sequence of Neisseria meningitidis serogroup B strain MC58.

The 2,272,351-base pair genome of Neisseria meningitidis strain MC58 (serogroup B), a causative agent of meningitis and septicemia, contains 2158 predicted coding regions, 1158 (53.7%) of which were assigned a biological role. Three major islands of horizontal DNA transfer were identified; two of these contain genes encoding proteins involved in pathogenicity, and the third island contains coding sequences only for hypothetical proteins. Insights into the commensal and virulence behavior of N. meningitidis can be gleaned from the genome, in which sequences for structural proteins of the pilus are clustered and several coding regions unique to serogroup B capsular polysaccharide synthesis can be identified. Finally, N. meningitidis contains more genes that undergo phase variation than any pathogen studied to date, a mechanism that controls their expression and contributes to the evasion of the host immune system.

Genome sequence of the radioresistant bacterium Deinococcus radiodurans R1.

The complete genome sequence of the radiation-resistant bacterium Deinococcus radiodurans R1 is composed of two chromosomes (2,648,638 and 412,348 base pairs), a megaplasmid (177,466 base pairs), and a small plasmid (45,704 base pairs), yielding a total genome of 3,284, 156 base pairs. Multiple components distributed on the chromosomes and megaplasmid that contribute to the ability of D. radiodurans to survive under conditions of starvation, oxidative stress, and high amounts of DNA damage were identified. Deinococcus radiodurans represents an organism in which all systems for DNA repair, DNA damage export, desiccation and starvation recovery, and genetic redundancy are present in one cell.

Genomic loci of the Porphyromonas gingivalis insertion element IS1126.

The Porphyromonas gingivalis genome contains multiple copies of insertion element IS1126. When chromosomal DNA digests of different strains were probed with IS1126, between 25 and 35 hybridizing fragments per genome were detected, depending on the strain. Unrelated strains had very different restriction fragment length polymorphism (RFLP) patterns. When different laboratory copies of a specific strain were examined, the IS1126 RFLP patterns were very similar but small differences were observed, indicating that element-associated changes had occurred during laboratory passage. Within the next year, genome sequencing, assembly, and annotation for P. gingivalis W83 will be completed. Because repetitive elements complicate the assembly of randomly sequenced DNA fragments, we isolated and sequenced the flanking regions of IS1126 copies in strain W83. We also isolated and sequenced the flanking regions of IS1126 copies in strain ATCC 33277 in order to compare insertion sites in phylogenetically divergent strains. We identified 37 new sequences flanking IS1126 from strain ATCC 33277 and 30 from strain W83. The insertion element was found between genes except where it transposed into another insertion element. Examination of identifiable flanking genes or open reading frames indicated that the insertion sites were different in the two strains, except that both strains possess an insertion adjacent to the Lys-gingipain gene (J. P. Lewis and F. L. Macrina, Infect. Immun. 66:3035-3042, 1998). Most of the genes or sequences flanking IS1126 in ATCC 33277 were present in W83 but were contiguous and not insertion element associated. Thus, where genes were identified in both strains, their order was maintained, indicating that the two genomes are organized similarly, but the loci of IS1126 are different. In both strains, insertion element-associated duplicated target sites were lost from several copies of IS1126, providing evidence of homologous recombination between elements. Larger organizational differences between the genomes, such as deletions and inversions, may result from insertion element-mediated recombination events.

Evidence for lateral gene transfer between Archaea and bacteria from genome sequence of Thermotoga maritima.

The 1,860,725-base-pair genome of Thermotoga maritima MSB8 contains 1,877 predicted coding regions, 1,014 (54%) of which have functional assignments and 863 (46%) of which are of unknown function. Genome analysis reveals numerous pathways involved in degradation of sugars and plant polysaccharides, and 108 genes that have orthologues only in the genomes of other thermophilic Eubacteria and Archaea. Of the Eubacteria sequenced to date, T. maritima has the highest percentage (24%) of genes that are most similar to archaeal genes. Eighty-one archaeal-like genes are clustered in 15 regions of the T. maritima genome that range in size from 4 to 20 kilobases. Conservation of gene order between T. maritima and Archaea in many of the clustered regions suggests that lateral gene transfer may have occurred between thermophilic Eubacteria and Archaea.

Alignment of whole genomes.

A new system for aligning whole genome sequences is described. Using an efficient data structure called a suffix tree, the system is able to rapidly align sequences containing millions of nucleotides. Its use is demonstrated on two strains of Mycoplasma tuberculosis, on two less similar species of Mycoplasma bacteria and on two syntenic sequences from human chromosome 12 and mouse chromosome 6. In each case it found an alignment of the input sequences, using between 30 s and 2 min of computation time. From the system output, information on single nucleotide changes, translocations and homologous genes can easily be extracted. Use of the algorithm should facilitate analysis of syntenic chromosomal regions, strain-to-strain comparisons, evolutionary comparisons and genomic duplications.

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.

DNA-PKcs: a T-cell tumour suppressor encoded at the mouse scid locus.

Severe combined immunodeficiency (SCID) mice are defective in their ability to rearrange their variable (V), diversity (D) and joining (J) genetic elements to generate functional immunoglobulin (Ig) and T-cell receptor (TCR) molecules; as a result, they lack mature B and T cells. These mice are highly sensitive to ionizing radiation, suggesting that the product of the scid gene plays a critical role in both V(D)J recombination and DNA double-strand break repair. Recent studies suggest that the SCID defect lies in the gene encoding the catalytic subunit of DNA-dependent protein kinase (DNA-PK; refs 6-8), a nuclear protein made up of the Ku 70 and Ku 86 subunits as well as the large catalytic subunit, DNA-PKcs. Other reports have implied that the SCID phenotype correlates with nonsense mutations at the extreme 3' end of Prkdc, the DNA-PKcs gene. The identity of the gene remains in doubt, however, because the consequences of genetic inactivation of Prkdc have not been determined. This study shows that complete inactivation of Prkdc in a novel insertional mouse mutant recapitulates the SCID phenotype and that Prkdc and scid are alleic. Significantly, DNA-PKcs null mice demonstrate complete penetrance of thymic lymphoblastic lymphomas, strongly suggesting that Prkdc functions in mice as a T-cell tumour suppressor and, by virtue of its association with DNA repair and recombination, belongs to the 'caretaker' class of tumour-suppressor genes that includes ATM, BRCA1 and BRCA2 (ref. 15).

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

Strategies for whole microbial genome sequencing and analysis.

The introduction of methods for automated DNA sequence analysis nearly a decade ago, together with more recent advances in the field of bioinformatics, have revolutionized biology and medicine and have ushered in a new era of genomic science, the study of genes and genomes. These new technologies have had an impact on many areas of research, including the association between genes and disease, in DNA-based diagnostics, and in the sequencing of genomes from human and other model organisms. The demonstration in 1995, that automated DNA sequencing methods could be used to decipher the entire genome sequence of a free-living organism, Haemophilus influenzae, was a milestone in both the genomics and microbial fields [1]. Since the first report of the complete sequence of H. influenzae, these methodologies have been adopted by laboratories around the world. The complete genomic sequence of five eubacterial species [1-5], one archaea [6], and the eukaryote, Saccharomyces cerevisiae [7], have been reported in the last 18 months. At the beginning of 1997 more than a dozen microbial genome projects are at or near completion, with many others in progress. It is likely that in the next few years we will see the complete sequence of perhaps as many as 30-40 microbial genomes. In this article, we will review methods for whole genome sequencing and analysis and examine how this information can be exploited to better understand microbial physiology and evolution.

DNA repeats identify novel virulence genes in Haemophilus influenzae.

The whole genome sequence (1.83 Mbp) of Haemophilus influenzae strain Rd was searched to identify tandem oligonucleotide repeat sequences. Loss or gain of one or more nucleotide repeats through a recombination-independent slippage mechanism is known to mediate phase variation of surface molecules of pathogenic bacteria, including H. influenzae. This facilitates evasion of host defenses and adaptation to the varying microenvironments of the host. We reasoned that iterative nucleotides could identify novel genes relevant to microbe-host interactions. Our search of the Rd genome sequence identified 9 novel loci with multiple (range 6-36, mean 22) tandem tetranucleotide repeats. All were found to be located within putative open reading frames and included homologues of hemoglobin-binding proteins of Neisseria, a glycosyltransferase (lgtC gene product) of Neisseria, and an adhesin of Yersinia. These tetranucleotide repeat sequences were also shown to be present in two other epidemiologically different H. influenzae type b strains, although the number and distribution of repeats was different. Further characterization of the lgtC gene showed that it was involved in phenotypic switching of a lipopolysaccharide epitope and that this variable expression was associated with changes in the number of tetranucleotide repeats. Mutation of lgtC resulted in attenuated virulence of H. influenzae in an infant rat model of invasive infection. These data indicate the rapidity, economy, and completeness with which whole genome sequences can be used to investigate the biology of pathogenic bacteria.

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.

Whole-genome random sequencing and assembly of Haemophilus influenzae Rd.

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.

Frequency and distribution of DNA uptake signal sequences in the Haemophilus influenzae Rd genome.

The naturally transformable, Gram-negative bacterium Haemophilus influenzae Rd preferentially takes up DNA of its own species by recognizing a 9-base pair sequence, 5'-AAGTGCGGT, carried in multiple copies in its chromosome. With the availability of the complete genome sequence, 1465 copies of the 9-base pair uptake site have been identified. Alignment of these sites unexpectedly reveals an extended consensus region of 29 base pairs containing the core 9-base pair region and two downstream 6-base pair A/T-rich regions, each spaced about one helix turn apart. Seventeen percent of the sites are in inverted repeat pairs, many of which are located downstream to gene termini and are capable of forming stem-loop structures in messenger RNA that might function as signals for transcription termination.