Michael Fonstein

Genome sequence of Bacillus cereus and comparative analysis with Bacillus anthracis

Bacillus cereus is an opportunistic pathogen causing food poisoning manifested by diarrhoeal or emetic syndromes. It is closely related to the animal and human pathogen Bacillus anthracis and the insect pathogen Bacillus thuringiensis, the former being used as a biological weapon and the latter as a pesticide. B. anthracis and B. thuringiensis are readily distinguished from B. cereus by the presence of plasmid-borne specific toxins (B. anthracis and B. thuringiensis) and capsule (B. anthracis). But phylogenetic studies based on the analysis of chromosomal genes bring controversial results, and it is unclear whether B. cereus, B. anthracis and B. thuringiensis are varieties of the same species or different species. Here we report the sequencing and analysis of the type strain B. cereus ATCC 14579. The complete genome sequence of B. cereus ATCC 14579 together with the gapped genome of B. anthracis A2012 enables us to perform comparative analysis, and hence to identify the genes that are conserved between B. cereus and B. anthracis, and the genes that are unique for each species. We use the former to clarify the phylogeny of the cereus group, and the latter to determine plasmid-independent species-specific markers.

Experimental Determination and System Level Analysis of Essential Genes in Escherichia coli MG1655

Defining the gene products that play an essential role in an organisms functional repertoire is vital to understanding the system level organization of living cells. We used a genetic footprinting technique for a genome-wide assessment of genes required for robust aerobic growth of Escherichia coli in rich media. We identified 620 genes as essential and 3,126 genes as dispensable for growth under these conditions. Functional context analysis of these data allows individual functional assignments to be refined. Evolutionary context analysis demonstrates a significant tendency of essential E. coli genes to be preserved throughout the bacterial kingdom. Projection of these data over metabolic subsystems reveals topologic modules with essential and evolutionarily preserved enzymes with reduced capacity for error tolerance.

WIT: integrated system for high-throughput genome sequence analysis and metabolic reconstruction

The WIT (What Is There) (http://wit.mcs.anl.gov/WIT2/) system has been designed to support comparative analysis of sequenced genomes and to generate metabolic reconstructions based on chromosomal sequences and metabolic modules from the EMP/MPW family of databases. This system contains data derived from about 40 completed or nearly completed genomes. Sequence homologies, various ORF-clustering algorithms, relative gene positions on the chromosome and placement of gene products in metabolic pathways (metabolic reconstruction) can be used for the assignment of gene functions and for development of overviews of genomes within WIT. The integration of a large number of phylogenetically diverse genomes in WIT facilitates the understanding of the physiology of different organisms.

Genome Sequence and Analysis of the Oral Bacterium Fusobacterium nucleatum Strain ATCC 25586

We present a complete DNA sequence and metabolic analysis of the dominant oral bacterium Fusobacterium nucleatum. Although not considered a major dental pathogen on its own, this anaerobe facilitates the aggregation and establishment of several other species including the dental pathogens Porphyromonas gingivalis and Bacteroides forsythus. The F. nucleatum strain ATCC 25586 genome was assembled from shotgun sequences and analyzed using the ERGO bioinformatics suite (http://www.integratedgenomics.com). The genome contains 2.17 Mb encoding 2,067 open reading frames, organized on a single circular chromosome with 27% GC content. Despite its taxonomic position among the gram-negative bacteria, several features of its core metabolism are similar to that of gram-positive Clostridium spp., Enterococcus spp., and Lactococcus spp. The genome analysis has revealed several key aspects of the pathways of organic acid, amino acid, carbohydrate, and lipid metabolism. Nine very-high-molecular-weight outer membrane proteins are predicted from the sequence, none of which has been reported in the literature. More than 137 transporters for the uptake of a variety of substrates such as peptides, sugars, metal ions, and cofactors have been identified. Biosynthetic pathways exist for only three amino acids: glutamate, aspartate, and asparagine. The remaining amino acids are imported as such or as di- or oligopeptides that are subsequently degraded in the cytoplasm. A principal source of energy appears to be the fermentation of glutamate to butyrate. Additionally, desulfuration of cysteine and methionine yields ammonia, H(2)S, methyl mercaptan, and butyrate, which are capable of arresting fibroblast growth, thus preventing wound healing and aiding penetration of the gingival epithelium. The metabolic capabilities of F. nucleatum revealed by its genome are therefore consistent with its specialized niche in the mouth.

The ERGOTM genome analysis and discovery system

The ERGO (http://ergo.integratedgenomics.com/ERGO/) genome analysis and discovery suite is an integration of biological data from genomics, biochemistry, high-throughput expression profiling, genetics and peer-reviewed journals to achieve a comprehensive analysis of genes and genomes. Far beyond any conventional systems that facilitate functional assignments, ERGO combines pattern-based analysis with comparative genomics by visualizing genes within the context of regulation, expression profiling, phylogenetic clusters, fusion events, networked cellular pathways and chromosomal neighborhoods of other functionally related genes. The result of this multifaceted approach is to provide an extensively curated database of the largest available integration of genomes, with a vast collection of reconstructed cellular pathways spanning all domains of life. Although access to ERGO is provided only under subscription, it is already widely used by the academic community. The current version of the system integrates 500 genomes from all domains of life in various levels of completion, 403 of which are available for subscription.

From Genetic Footprinting to Antimicrobial Drug Targets: Examples in Cofactor Biosynthetic Pathways

Novel drug targets are required in order to design new defenses against antibiotic-resistant pathogens. Comparative genomics provides new opportunities for finding optimal targets among previously unexplored cellular functions, based on an understanding of related biological processes in bacterial pathogens and their hosts. We describe an integrated approach to identification and prioritization of broad-spectrum drug targets. Our strategy is based on genetic footprinting in Escherichia coli followed by metabolic context analysis of essential gene orthologs in various species. Genes required for viability of E. coli in rich medium were identified on a whole-genome scale using the genetic footprinting technique. Potential target pathways were deduced from these data and compared with a panel of representative bacterial pathogens by using metabolic reconstructions from genomic data. Conserved and indispensable functions revealed by this analysis potentially represent broad-spectrum antibacterial targets. Further target prioritization involves comparison of the corresponding pathways and individual functions between pathogens and the human host. The most promising targets are validated by direct knockouts in model pathogens. The efficacy of this approach is illustrated using examples from metabolism of adenylate cofactors NAD(P), coenzyme A, and flavin adenine dinucleotide. Several drug targets within these pathways, including three distantly related adenylyltransferases (orthologs of the E. coli genes nadD, coaD, and ribF), are discussed in detail.

The Rhodobacter capsulatus genome

The genome of Rhodobacter capsulatus has been completely sequenced. It consists of a single chromosome containing 3.5 Mb and a circular plasmid of 134 kb. This effort, started in 1992, began with a fine-structure restriction map of an overlapping set of cosmids that covered the genome. Cosmid sequencing led to a gapped genome that was filled by primer walking on the chromosome and by using lambda clones. Methods had to be developed to handle strong stops in the high GC (68%) inserts. Annotation was done with the ERGO system at Integrated Genomics, as was the reconstruction of the cells metabolism. It was possible to recognize 3709 orfs of which functional assignments could be made with high confidence to 2392 (65%). Unusual features include the presence of numerous cryptic phage genomes embedded in the chromosome.

Refinement of the high-resolution physical and genetic map of Rhodobacter capsulatus and genome surveys using blots of the cosmid encyclopedia.

Cosmids from a library containing Rhodobacter capsulatus DNA fragments were previously ordered in two contigs: one corresponding to the chromosome and one to a 134 kb plasmid. This map contained 40 regions connected only by colony hybridization. To confirm the linkage and correct the map, the actual sizes of the overlaps were determined by blot‐hybridization with Rhodobacter chromosomal DNA and by mapping of additional cosmids. Several revisions of the earlier map include single cosmid shifts and inversions. One additional gap in a cosmid contig was also found, raising the possibility that the chromosome is not a contiguous circle. About 2500 additional EcoRI,BamHI and HindIII restriction sites were added to the 560 EcoRV sites previously mapped onto the Rhodobacter chromosome, increasing the resolution of the physical map to the size of individual genes. Twenty‐five new markers were located on the genetic map. The 48 markers now mapped represent nearly 300 genes and ORFs cloned from different species of Rhodobacter. The orientation of transcription of the four rrn operons was established using 16S rRNA‐ and 23S rRNA‐specific probes and digestion with the rare‐cutting enzyme, CeuI. Gel blots of 192 cosmids of the miniset of R.capsulatus digested with EcoRV were prepared. Such a hybridization template represents the whole genome cut into 560 DNA fragments varying in size from 0.4 to 25 kb. This template was used for high‐resolution mapping of single genes, analysis of total genomic DNAs from related Rhodobacter strains and differentially expressed RNAs.

Rhodobacter capsulatus SB1003

DNA fragments from a Sau3a partial digest of Rhodobacter capsulatus chromosomal DNA were size selected and cloned in the BamRI site of Lorist 6 by cosmid arm cloning and clones were individually transferred to eighteen 96-well plates (see also Chapter 28). Using a specially constructed printing device, individual cosmid clones were replicated onto nylon filters as ordered sets (18X coverage). These sets of clones were hybridized with individual XbaI and AseI fragments prepared from PFGE gels. Forty-one restriction sites for AseI and XbaI were mapped onto the 3.7-Mb genome by standard analyses of PFGE blot-hybridization data and of linking cosmids, revealed as clones hybridizing with two restriction fragments at once (Fonstein et al., 1992). Other types of hybridization (fragments to fragments, fragments to cosmid sets and riboprobes generated from the ends of the cosmids to cosmid sets) made it possible to vary the probe/target ratio in the hybridizations, thus distinguishing mapping artifacts connected with different repeated DNA sequences. At this stage, the cosmid clones were grouped in about 80 subcontigs, corresponding to the macro-restriction fragments from the PFG and the regions surrounding these rare sites. They formed two groups, one corresponding to the chromosome of R. capsulatus, the other to its 134-kb plasmid.

Encyclopedias of Bacterial Genomes

Summarizing recent studies of bacterial chromosomes, participants in the First International Symposium on Small Genomes in Paris in 1993 declared the birth of a new discipline, bacterial genomics. The subject of genomics can be defined as the study of integral genome structures, integral genome properties and the evolution of genomes. Genome encyclopedias (ordered sets of overlapping clones containing the entire genome), together with Pulsed/Field gel electrophoresis (PFGE)-related methods (see Chapters 24 to 26), and genome sequencing (see Chapts in section 2 of this monograph) are the specific tools of genomics.