Comparative genome sequencing for discovery of novel polymorphisms in Bacillus anthracis

Genome differences that distinguish Bacillus anthracis from Bacillus cereus and Bacillus thuringiensis

The three species of the group 1 bacilli, Bacillus anthracis, B. cereus, and B. thuringiensis, are genetically very closely related. All inhabit soil habitats but exhibit different phenotypes. B. anthracis is the causative agent of anthrax and is phylogenetically monomorphic, while B. cereus and B. thuringiensis are genetically more diverse. An amplified fragment length polymorphism analysis described here demonstrates genetic diversity among a collection of non-anthrax-causing Bacillus species, some of which show significant similarity to B. anthracis. Suppression subtractive hybridization was then used to characterize the genomic differences that distinguish three of the non-anthrax-causing bacilli from B. anthracis Ames. Ninety-three DNA sequences that were present in B. anthracis but absent from the non-anthrax-causing Bacillus genomes were isolated. Furthermore, 28 of these sequences were not found in a collection of 10 non-anthrax-causing Bacillus species but were present in all members of a representative collection of B. anthracis strains. These sequences map to distinct loci on the B. anthracis genome and can be assayed simultaneously in multiplex PCR assays for rapid and highly specific DNA-based detection of B. anthracis.

The genome sequence of Bacillus anthracis Ames and comparison to closely related bacteria

Bacillus anthracis is an endospore-forming bacterium that causes inhalational anthrax. Key virulence genes are found on plasmids (extra-chromosomal, circular, double-stranded DNA molecules) pXO1 (ref. 2) and pXO2 (ref. 3). To identify additional genes that might contribute to virulence, we analysed the complete sequence of the chromosome of B. anthracis Ames (about 5.23 megabases). We found several chromosomally encoded proteins that may contribute to pathogenicity--including haemolysins, phospholipases and iron acquisition functions--and identified numerous surface proteins that might be important targets for vaccines and drugs. Almost all these putative chromosomal virulence and surface proteins have homologues in Bacillus cereus, highlighting the similarity of B. anthracis to near-neighbours that are not associated with anthrax. By performing a comparative genome hybridization of 19 B. cereus and Bacillus thuringiensis strains against a B. anthracis DNA microarray, we confirmed the general similarity of chromosomal genes among this group of close relatives. However, we found that the gene sequences of pXO1 and pXO2 were more variable between strains, suggesting plasmid mobility in the group. The complete sequence of B. anthracis is a step towards a better understanding of anthrax pathogenesis.

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.

The sigmaE regulon and the identification of additional sporulation genes in Bacillus subtilis

We report the identification and characterization on a genome-wide basis of genes under the control of the developmental transcription factor sigma(E) in Bacillus subtilis. The sigma(E) factor governs gene expression in the larger of the two cellular compartments (the mother cell) created by polar division during the developmental process of sporulation. Using transcriptional profiling and bioinformatics we show that 253 genes (organized in 157 operons) appear to be controlled by sigma(E). Among these, 181 genes (organized in 121 operons) had not been previously described as members of this regulon. Promoters for many of the newly identified genes were located by transcription start site mapping. To assess the role of these genes in sporulation, we created null mutations in 98 of the newly identified genes and operons. Of the resulting mutants, 12 (in prkA, ybaN, yhbH, ykvV, ylbJ, ypjB, yqfC, yqfD, ytrH, ytrI, ytvI and yunB) exhibited defects in spore formation. In addition, subcellular localization studies were carried out using in-frame fusions of several of the genes to the coding sequence for GFP. A majority of the fusion proteins localized either to the membrane surrounding the developing spore or to specific layers of the spore coat, although some fusions showed a uniform distribution in the mother cell cytoplasm. Finally, we used comparative genomics to determine that 46 of the sigma(E)-controlled genes in B.subtilis were present in all of the Gram-positive endospore-forming bacteria whose genome has been sequenced, but absent from the genome of the closely related but not endospore-forming bacterium Listeria monocytogenes, thereby defining a core of conserved sporulation genes of probable common ancestral origin. Our findings set the stage for a comprehensive understanding of the contribution of a cell-specific transcription factor to development and morphogenesis.

Methionine regeneration and aminotransferases in Bacillus subtilis, Bacillus cereus, and Bacillus anthracis

The conversion of ketomethiobutyrate to methionine has been previously examined in a number of organisms, wherein the aminotransferases responsible for the reaction have been found to be members of the Ia subfamily (L. C. Berger, J. Wilson, P. Wood, and B. J. Berger, J. Bacteriol. 183:4421-4434, 2001). The genome of Bacillus subtilis has been found to contain no subfamily Ia aminotransferase sequences. Instead, the analogous enzymes in B. subtilis were found to be members of the If subfamily. These putative aspartate aminotransferases, the yugH, ywfG, ykrV, aspB, and patA gene products, have been cloned, expressed, and characterized for methionine regeneration activity. Only YkrV was able to convert ketomethiobutyrate to methionine, and it catalyzed the reaction only when glutamine was used as amino donor. In contrast, subcellular homogenates of B. subtilis and Bacillus cereus utilized leucine, isoleucine, valine, alanine, phenylalanine, and tyrosine as effective amino donors. The two putative branched-chain aminotransferase genes in B. subtilis, ybgE and ywaA, were also cloned, expressed, and characterized. Both gene products effectively transaminated branched-chain amino acids and ketoglutarate, but only YbgE converted ketomethiobutyrate to methionine. The amino donor preference for methionine regeneration by YbgE was found to be leucine, isoleucine, valine, phenylalanine, and tyrosine. The B. subtilis ybgE gene is a member of the family III of aminotransferases and falls in a subfamily designated here IIIa. Examination of B. cereus and Bacillus anthracis genome data found that there were no subfamily IIIa homologues in these organisms. In both B. cereus and B. anthracis, two putative branched-chain aminotransferases and two putative D-amino acid aminotransferases were discovered as members of subfamily IIIb. These four sequences were cloned from B. cereus, expressed, and characterized. Only the gene product from the sequence designated Bc-BCAT2 was found to convert ketomethiobutyrate to methionine, with an amino donor preference of leucine, isoleucine, valine, phenylalanine, and tyrosine. The B. anthracis homologue of Bc-BCAT2 was also cloned, expressed, and characterized and was found to be identical in activity. The aminooxy compound canaline was found to be an uncompetitive inhibitor of B. subtilis YbgE and also inhibited growth of B. subtilis and B. cereus in culture.

Identification of the immunodominant protein and other proteins of the Bacillus anthracis exosporium

Spores of Bacillus anthracis, the causative agent of anthrax, are enclosed by a prominent loose-fitting, balloon-like layer called the exosporium. Although the exosporium serves as the source of surface antigens and a primary permeability barrier of the spore, its molecular structure and function are not well characterized. In this study, we identified five major proteins in purified B. anthracis (Sterne strain) exosporia. One protein was the recently identified collagen-like glycoprotein BclA, which appears to be a structural component of the exosporium hair-like nap. Using a large panel of unique antispore monoclonal antibodies, we demonstrated that BclA is the immunodominant antigen on the B. anthracis spore surface. We also showed that the BclA protein and not a carbohydrate constituent directs the dominant immune response. In addition, the length of the central (GXX)(n) repeat region of BclA appears to be strain specific. Two other unique proteins, BxpA and BxpB, were identified. BxpA is unusually rich in Gln and Pro residues and contains several different tandem repeats, which also exhibit strain-specific variation. In addition, BxpA was found to be cleaved approximately in half. BxpB appears to be glycosylated or associated with glycosylated material and is encoded by a gene that (along with bclA) may be part of an exosporium genomic island. The other two proteins identified were alanine racemase and superoxide dismutase, both of which were reported to be associated with the surface of other Bacillus spores. Possible functions of the newly identified proteins are discussed.

Crystal structure of a hypothetical protein, TM841 of Thermotoga maritima, reveals its function as a fatty acid-binding protein

We determined the three-dimensional (3D) crystal structure of protein TM841, a protein product from a hypothetical open-reading frame in the genome of the hyperthermophile bacterium Thermotoga maritima, to 2.0 A resolution. The protein belongs to a large protein family, DegV or COG1307 of unknown function. The 35 kDa protein consists of two separate domains, with low-level structural resemblance to domains from other proteins with known 3D structures. These structural homologies, however, provided no clues for the function of TM841. But the electron density maps revealed clear density for a bound fatty-acid molecule in a pocket between the two protein domains. The structure indicates that TM841 has the molecular function of fatty-acid binding and may play a role in the cellular functions of fatty acid transport or metabolism.

Bacillus anthracis virulence in Guinea pigs vaccinated with anthrax vaccine adsorbed is linked to plasmid quantities and clonality

Bacillus anthracis is a bacterial pathogen of great importance, both historically and in the present. This study presents data collected from several investigations and indicates that B. anthracis virulence is associated with the clonality and virulence of plasmids pXO1 and pXO2. Guinea pigs vaccinated with Anthrax Vaccine Adsorbed were challenged with 20 B. anthracis isolates representative of worldwide genetic diversity. These same isolates were characterized with respect to plasmid copy number by using a novel method of quantitative PCR developed for rapid and efficient detection of B. anthracis from environmental samples. We found that the copy numbers for both pXO1 and pXO2 differed from those in previously published reports. By combining the data on survival, plasmid copy numbers, and clonality, we developed a model predicting virulence. This model was validated by using a randomly chosen set of 12 additional B. anthracis isolates. Results from this study will be helpful in future efforts to elucidate the basis for variation in the virulence of this important pathogen.

Microbe forensics: oxygen and hydrogen stable isotope ratios in Bacillus subtilis cells and spores

Bacillus subtilis, a Gram-positive, endospore-forming soil bacterium, was grown in media made with water of varying oxygen (delta(18)O) and hydrogen (deltaD) stable isotope ratios. Logarithmically growing cells and spores were each harvested from the cultures and their delta(18)O and deltaD values determined. Oxygen and hydrogen stable isotope ratios of organic matter were linearly related with those of the media water. We used the relationships determined in these experiments to calculate the effective whole-cell fractionation factors between water and organic matter for B. subtilis. We then predicted the delta(18)O and deltaD values of spores produced in nutritionally identical media and local water sources for five different locations around the United States. Each of the measured delta(18)O and deltaD values of the spores matched the predicted values within a 95% confidence interval, indicating that stable isotope ratio analyses may be a powerful tool for tracing the geographic point-of-origin for microbial products.

Analysis of the quorum-sensing regulon of the opportunistic pathogen Burkholderia cepacia H111 by proteomics

Burkholderia cepacia H111, an important pathogen for persons suffering from cystic fibrosis, employs a quorum-sensing (QS) system, cep, to control expression of virulence factors as well as the formation of biofilms. The QS system is thought to ensure that pathogenic traits are only expressed when the bacterial population density is high enough to overwhelm the host before it is able to mount an efficient response. In this study, we compared the protein pattern of the intracellular, extracellular, and surface protein fractions of an AHL-deficient cepI mutant with the one of the parent strain H111 by means of two-dimensional gel electrophoresis (2-DE). Our analysis showed that 55 proteins out of 985 detected spots were differentially expressed; these are expected to represent QS-controlled gene products. Addition of the respective signal molecules to the growth medium of the cep mutant fully restored the wild-type protein expression profile. In total about 5% of the B. cepacia proteome was downregulated and 1% upregulated in the cepI mutant, indicating that quorum sensing represents a global regulatory system. Nineteen proteins were identified with high confidence by N-terminal sequence analysis.

NotI passporting to identify species composition of complex microbial systems

We describe here a new method for large-scale scanning of microbial genomes on a quantitative and qualitative basis. To achieve this aim we propose to create NotI passports: databases containing NotI tags. We demonstrated that these tags comprising 19 bp of sequence information could be successfully generated using DNA isolated from intestinal or fecal samples. Such NotI passports allow the discrimination between closely related bacterial species and even strains. This procedure for generating restriction site tagged sequences (RSTS) is called passporting and can be adapted to any other rare cutting restriction enzyme. A comparison of 1312 tags from available sequenced Escherichia coli genomes, generated with the NotI, PmeI and SbfI restriction enzymes, revealed only 219 tags that were not unique. None of these tags matched human or rodent sequences. Therefore the approach allows analysis of complex microbial mixtures such as in human gut and identification with high accuracy of a particular bacterial strain on a quantitative and qualitative basis.

Bacillus anthracis diverges from related clades of the Bacillus cereus group in 16S-23S ribosomal DNA intergenic transcribed spacers containing tRNA genes

Mung bean nuclease treatment of 16S-23S ribosomal DNA intergenic transcribed spacers (ITS) amplified from several strains of the six species of the Bacillus cereus group showed that B. anthracis Davis TE702 and B. mycoides G2 have other intermediate fragments in addition to the 220- and 550-bp homoduplex fragments typical of the B. cereus group. Long and intermediate homoduplex ITS fragments from strains Davis TE702 and G2 and from another 19 strains of the six species were sequenced. Two main types of ITS were found, either with two tRNA genes (tRNA(Ile) and tRNA(Ala)) or without any at all. Strain Davis TE702 harbors an additional ITS with a single tRNA gene, a hybrid between the tRNA(Ile) and tRNA(Ala) genes, suggesting that a recombination event rather than a deletion generated the single tDNA-containing ITS. Strain G2 showed an additional ITS of intermediate length with no tDNA and no similarity to other known sequences. Neighbor-joining analysis of tDNA-containing long ITS indicated that B. cereus and B. thuringiensis represent a single clade. Three signature sequences discriminated B. anthracis from B. cereus and B. thuringiensis, indicating that the anthrax agent started evolving separately from the related clades of the B. cereus group. B. mycoides and B. weienstephanensis were very closely related, while B. pseudomycoides appeared the most distant species.

Genomewide pattern of synonymous nucleotide substitution in two complete genomes of Mycobacterium tuberculosis

Comparison of the pattern of synonymous nucleotide substitution between two complete genomes of Mycobacterium tuberculosis at 3298 putatively orthologous loci showed a mean percent difference per synonymous site of 0.000328 0.000022. Although 80.5% of loci showed no synonymous or nonsynonymous nucleotide differences, the level of polymorphism observed at other loci was greater than suggested by previous studies of a small number of loci. This level of nucleotide difference leads to the conservative estimate that the common ancestor of these two genotypes occurred approximately 35000 ago, which is twice as high as some recent estimates of the time of origin of this species. Our results suggest that a large number of loci should be examined for an accurate assessment of the level of nucleotide diversity in natural populations of pathogenic microorganisms.

Protocol for real-time PCR identification of anthrax spores from nasal swabs after broth enrichment

A mass-screening protocol for the diagnosis of anthrax from nasal swabs based on an enrichment step in liquid medium was devised. Incubation for growth was performed in autoclavable vials and racks which allow real-time PCR analysis of sterilized cultures. A dual-color PCR was set up with primers and probes for the chromosomal marker rpoB and the plasmid marker lef. Specific primer and probe sets were designed for the differentiation of Bacillus anthracis from B. cereus and for the differentiation of the Sterne vaccine strain from field isolates and the Ames strain, which was used in the recent anthrax bioterrorist attack. The present protocol thus combines the high specificity and sensitivity of real-time PCR with excellent biosafety and the low hands-on time necessary for the processing of large numbers of samples, which is extremely important during control programs involving the processing of large numbers of samples.

Molecular subtyping of Bacillus anthracis and the 2001 bioterrorism-associated anthrax outbreak, United States

Molecular subtyping of Bacillus anthracis played an important role in differentiating and identifying strains during the 2001 bioterrorism-associated outbreak. Because B. anthracis has a low level of genetic variability, only a few subtyping methods, with varying reliability, exist. We initially used multiple-locus variable-number tandem repeat analysis (MLVA) to subtype 135 B. anthracis isolates associated with the outbreak. All isolates were determined to be of genotype 62, the same as the Ames strain used in laboratories. We sequenced the protective antigen gene (pagA) from 42 representative outbreak isolates and determined they all had a pagA sequence indistinguishable from the Ames strain (PA genotype I). MLVA and pagA sequencing were also used on DNA from clinical specimens, making subtyping B. anthracis possible without an isolate. Use of high-resolution molecular subtyping determined that all outbreak isolates were indistinguishable by the methods used and probably originated from a single source. In addition, subtyping rapidly identified laboratory contaminants and nonoutbreak-related isolates.