Comparative genome sequencing for discovery of novel polymorphisms in Bacillus anthracis

Cloning, purification and crystallization of Bacillus anthracis class C acid phosphatase

Cloning, expression, purification and crystallization studies of a recombinant class C acid phosphatase from the Category A pathogen Bacillus anthracis are reported. Large diffraction-quality crystals were grown in the presence of HEPES and Jeffamine ED-2001 at pH 7.0. The crystals belong to space group P2(1)2(1)2(1), with unit-cell parameters a = 53.4, b = 90.1, c = 104.2 angstroms. The asymmetric unit is predicted to contain two protein molecules with a solvent content of 38%. Two native data sets were collected from the same crystal before and after flash-annealing. The first data set had a mosaicity of 1.6 degrees and a high-resolution limit of 1.8 angstroms. After flash-annealing, the apparent mosaicity decreased to 0.9 degrees and the high-resolution limit of usable data increased to 1.6 angstroms. This crystal form is currently being used to determine the structure of B. anthracis class C acid phosphatase with experimental phasing techniques.

Growth characteristics of virulent Bacillus anthracis and potential surrogate strains

The objectives of this study were to compare generation and lag times of virulent Bacillus anthracis strains with those of other Bacillus strains, to identify possible surrogates for growth studies, and to determine if the B. cereus module of the U.S. Department of Agriculture Pathogen Modeling Program (PMP) had predictive value for B. anthracis. Growth characteristics of B. anthracis, B. cereus, B. mycoides, and B. subtilis strains in brain heart infusion broth at pH 6.5, 6.0, and 5.5 were determined by absorbance measurements. Growth curves of B. anthracis Sterne and B. cereus strains appeared similar, and the generation times for strain Sterne fell within the PMP's 95% confidence interval for B. cereus. However, the virulent B. anthracis strains Vollum and Pasteur had shorter generation times than the avirulent Sterne strain and most other surrogates and were lower than the PMP's 95% confidence interval for B. cereus. Growth curves of B. cereus ATCC 9818 and B. subtilis ATCC 6633 were more similar to those of virulent B. anthracis strains, but all potential surrogates had significantly different generation times and lag times under some conditions.

PulseNet USA: a five-year update

PulseNet USA is the molecular surveillance network for foodborne infections in the United States. Since its inception in 1996, it has been instrumental in detection, investigation and control of numerous outbreaks caused by Shiga toxin-producing Escherichia coli O157:[H7] (STEC O157), Salmonella enterica, Listeria monocytogenes, Shigella spp., and Campylobacter. This paper describes the current status of the network, including the methodologies used and its future possibilities. The currently preferred subtyping method in the network is pulsed-field gel electrophoresis (PFGE), a proven highly discriminatory molecular subtyping method. New simpler sequencebased subtyping methods are under development and validation to complement and eventually replace PFGE. PulseNet is essentially a cluster detection network, but the data in the system will now also be used in attribution analyses of sporadic infections. The PulseNet platform will also be used as a primary tool in preparedness and response to acts of food bioterrorism.

Single-nucleotide repeat analysis for subtyping Bacillus anthracis isolates

Single-nucleotide repeats (SNRs) are variable-number tandem repeats that display very high mutation rates. In an outbreak situation, the use of a marker system that exploits regions with very high mutation rates, such as SNRs, allows the differentiation of isolates with extremely low levels of genetic diversity. This report describes the identification and analysis of SNR loci of Bacillus anthracis. SNR loci were selected in silico, and the loci with the highest diversity were used to design and test locus-specific primers against a number of B. anthracis strains with the same multilocus variable-number tandem repeat analysis (MLVA) genotype. SNR markers that allowed strains with the same MLVA genotype to be differentiated from each other were identified. The resulting SNR marker system can be used as a molecular epidemiological tool in a natural outbreak or bioterrorism event, offering the best chance of distinguishing very closely related isolates.

Observing local and global properties of metabolic pathways: ‘load points’ and ‘choke points’ in the metabolic networks

MOTIVATION

The local and global aspects of metabolic network analyses allow us to identify enzymes or reactions that are crucial for the survival of the organism(s), therefore directing us towards the discovery of potential drug targets.

RESULTS

We demonstrate a new method ('load points') to rank the enzymes/metabolites in the metabolic network and propose a model to determine and rank the biochemical lethality in metabolic networks (enzymes/metabolites) through 'choke points'. Based on an extended form of the graph theory model of metabolic networks, metabolite structural information was used to calculate the k-shortest paths between metabolites (the presence of more than one competing path between substrate and product). On the basis of these paths and connectivity information, load points were calculated and used to empirically rank the importance of metabolites/enzymes in the metabolic network. The load point analysis emphasizes the role that the biochemical structure of a metabolite, rather than its connectivity (hubs), plays in the conversion pathway. In order to identify potential drug targets (based on the biochemical lethality of metabolic networks), the concept of choke points and load points was used to find enzymes (edges) which uniquely consume or produce a particular metabolite (nodes). A non-pathogenic bacterial strain Bacillus subtilis 168 (lactic acid producing bacteria) and a related pathogenic bacterial strain Bacillus anthracis Sterne (avirulent but toxigenic strain, producing the toxin Anthrax) were selected as model organisms. The choke point strategy was implemented on the pathogen bacterial network of B.anthracis Sterne. Potential drug targets are proposed based on the analysis of the top 10 choke points in the bacterial network. A comparative study between the reported top 10 bacterial choke points and the human metabolic network was performed. Further biological inferences were made on results obtained by performing a homology search against the human genome.

AVAILABILITY

The load and choke point modules are introduced in the Pathway Hunter Tool (PHT), the basic version of which is available on http://www.pht.uni-koeln.de.

A novel FtsZ-like protein is involved in replication of the anthrax toxin-encoding pXO1 plasmid in Bacillus anthracis

Plasmid pXO1 encodes the tripartite anthrax toxin, which is the major virulence factor of Bacillus anthracis. In spite of the important role of pXO1 in anthrax pathogenesis, very little is known about its replication and maintenance in B. anthracis. We cloned a 5-kb region of the pXO1 plasmid into an Escherichia coli vector and showed that this plasmid can replicate when introduced into B. anthracis. Mutational analysis showed that open reading frame 45 (repX) of pXO1 was required for the replication of the miniplasmid in B. anthracis. Interestingly, repX showed limited homology to bacterial FtsZ proteins that are involved in cell division. A mutation in the predicted GTP binding domain of RepX abolished its replication activity. Genes almost identical to repX are contained on several megaplasmids in members of the Bacillus cereus group, including a B. cereus strain that causes an anthrax-like disease. Our results identify a novel group of FtsZ-related initiator proteins that are required for the replication of virulence plasmids in B. anthracis and possibly in related organisms. Such replication proteins may provide novel drug targets for the elimination of plasmids encoding the anthrax toxin and other virulence factors.

Probing genomic diversity and evolution of Escherichia coli O157 by single nucleotide polymorphisms

Infections by Shiga toxin-producing Escherichia coli O157:H7 (STEC O157) are the predominant cause of bloody diarrhea and hemolytic uremic syndrome in the United States. In silico comparison of the two complete STEC O157 genomes (Sakai and EDL933) revealed a strikingly high level of sequence identity in orthologous protein-coding genes, limiting the use of nucleotide sequences to study the evolution and epidemiology of this bacterial pathogen. To systematically examine single nucleotide polymorphisms (SNPs) at a genome scale, we designed comparative genome sequencing microarrays and analyzed 1199 chromosomal genes (a total of 1,167,948 bp) and 92,721 bp of the large virulence plasmid (pO157) of eleven outbreak-associated STEC O157 strains. We discovered 906 SNPs in 523 chromosomal genes and observed a high level of DNA polymorphisms among the pO157 plasmids. Based on a uniform rate of synonymous substitution for Escherichia coli and Salmonella enterica (4.7x10(-9) per site per year), we estimate that the most recent common ancestor of the contemporary beta-glucuronidase-negative, non-sorbitolfermenting STEC O157 strains existed ca. 40 thousand years ago. The phylogeny of the STEC O157 strains based on the informative synonymous SNPs was compared to the maximum parsimony trees inferred from pulsed-field gel electrophoresis and multilocus variable numbers of tandem repeats analysis. The topological discrepancies indicate that, in contrast to the synonymous mutations, parts of STEC O157 genomes have evolved through different mechanisms with highly variable divergence rates. The SNP loci reported here will provide useful genetic markers for developing high-throughput methods for fine-resolution genotyping of STEC O157. Functional characterization of nucleotide polymorphisms should shed new insights on the evolution, epidemiology, and pathogenesis of STEC O157 and related pathogens.

The Bacillus anthracis chromosome contains four conserved, excision-proficient, putative prophages

Background

Bacillus anthracis is considered to be a recently emerged clone within the Bacillus cereus sensu lato group. The B. anthracis genome sequence contains four putative lambdoid prophages. We undertook this study in order to understand whether the four prophages are unique to B. anthracis and whether they produce active phages.

Results

More than 300 geographically and temporally divergent isolates of B. anthracis and its near neighbors were screened by PCR for the presence of specific DNA sequences from each prophage region. Every isolate of B. anthracis screened by PCR was found to produce all four phage-specific amplicons whereas none of the non-B. anthracis isolates, produced more than one phage-specific amplicon. Excision of prophages could be detected by a PCR based assay for attP sites on extra-chromosomal phage circles and for attB sites on phage-excised chromosomes. SYBR-green real-time PCR assays indicated that prophage excision occurs at very low frequencies (2 × 10^-5 - 8 × 10^-8/cell). Induction with mitomycin C increased the frequency of excision of one of the prophages by approximately 250 fold. All four prophages appear to be defective since, mitomycin C induced culture did not release any viable phage particle or lyse the cells or reveal any phage particle under electron microscopic examination.

Conclusion

The retention of all four putative prophage regions across all tested strains of B. anthracis is further evidence of the very recent emergence of this lineage and the prophage regions may be useful for differentiating the B. anthracis chromosome from that of its neighbors. All four prophages can excise at low frequencies, but are apparently defective in phage production.

Diagnostic oligonucleotide microarray fingerprinting of Bacillus isolates

A genome-independent microarray and new statistical techniques were used to genotype Bacillus strains and quantitatively compare DNA fingerprints with the known taxonomy of the genus. A synthetic DNA standard was used to understand process level variability and lead to recommended standard operating procedures for microbial forensics and clinical diagnostics.

Genetic analysis and attribution of microbial forensics evidence

Because of the availability of pathogenic microorganisms and the relatively low cost of preparing and disseminating bioweapons, there is a continuing threat of biocrime and bioterrorism. Thus, enhanced capabilities are needed that enable the full and robust forensic exploitation and interpretation of microbial evidence from acts of bioterrorism or biocrimes. To respond to the need, greater resources and efforts are being applied to the burgeoning field of microbial forensics. Microbial forensics focuses on the characterization, analysis and interpretation of evidence for attributional purposes from a bioterrorism act, biocrime, hoax or inadvertent agent release. To enhance attribution capabilities, a major component of microbial forensics is the analysis of nucleic acids to associate or eliminate putative samples. The degree that attribution can be addressed depends on the context of the case, the available knowledge of the genetics, phylogeny, and ecology of the target microorganism, and technologies applied. The types of genetic markers and features that can impact statistical inferences of microbial forensic evidence include: single nucleotide polymorphisms, repetitive sequences, insertions and deletions, mobile elements, pathogenicity islands, virulence and resistance genes, house keeping genes, structural genes, whole genome sequences, asexual and sexual reproduction, horizontal gene transfer, conjugation, transduction, lysogeny, gene conversion, recombination, gene duplication, rearrangements, and mutational hotspots. Nucleic acid based typing technologies include: PCR, real-time PCR, MLST, MLVA, whole genome sequencing, and microarrays.

Characterization of germination receptors of Bacillus cereus ATCC 14579

Specific amino acids, purine ribonucleosides, or a combination of the two is required for efficient germination of endospores of Bacillus cereus ATCC 14579. A survey including 20 different amino acids showed that l-alanine, l-cysteine, l-threonine, and l-glutamine are capable of initiating the germination of endospores of B. cereus ATCC 14579. In addition, the purine ribonucleosides inosine and adenosine can trigger germination of the spores. Advanced annotation of the B. cereus ATCC 14579 genome revealed the presence of seven putative germination (ger) operons, termed gerG, gerI, gerK, gerL, gerQ, gerR, and gerS. To determine the role of the encoded putative receptors in nutrient-induced germination, disruption mutants were constructed by the insertion of pMUTIN4 into each of the seven operons. Four of the seven mutants were affected in the germination response to amino acids or purine ribonucleosides, whereas no phenotype could be attributed to the mutants with disrupted gerK, gerL, and gerS loci. The strain with a disrupted gerR operon was severely hampered in the ability to germinate: germination occurred in response to l-glutamine but not in the presence of any of the other amino acids tested. The gerG mutant showed significantly reduced l-glutamine-induced germination, which points to a role of this receptor in the l-glutamine germination signaling pathway. gerR, gerI, and gerQ mutants showed reduced germination rates in the presence of inosine, suggesting a role for these operons in ribonucleoside signaling. Efficient germination by the combination of l-glutamine and inosine was shown to involve the gerG and gerI operons, since the germination of mutants lacking either one of these receptors was significantly reduced. Germination triggered by the combination of l-phenylalanine and inosine was lost in the gerI mutant, indicating that both molecules are effective at the GerI receptor.

The impact of comparative genomics on infectious disease research

The past decade has witnessed a revolution in infectious disease research, fuelled by the accumulation of a huge amount of DNA sequence data. The avalanche of genome sequence information has largely promoted the development of comparative genomics, which exploits available genome sequences to perform either inter- or intra-species comparisons of bacterial genome contents, or performs comparisons between the human genome and those of other organisms. This review aims to summarize how comparative genomics is being extensively used in infectious disease research, such as in the studies to identify virulence determinants, antimicrobial drug targets, vaccine candidates and new markers for diagnostics. These applications hold considerable promise for alleviating the burden of infectious diseases in the coming years.

Genetic analysis of transfer-related regions of the vancomycin resistance Enterococcus conjugative plasmid pHTbeta: identification of oriT and a putative relaxase gene

The pHT plasmids pHTalpha (65.9 kbp), pHTbeta (63.7 kbp), and pHTgamma (66.5 kbp) are highly conjugative pheromone-independent pMG1-like plasmids that carry Tn1546-like transposons encoding vancomycin resistance. pHTbeta is the prototype plasmid, and the pHTalpha and pHTgamma plasmids are derivatives of the insertion into pHTbeta of an IS232-like (2.2 kbp) element and a group II intron (2.8 kbp), respectively. The complete nucleotide sequence of the pHTbeta plasmid was determined and, with the exception of the Tn1546-like insertion (10,851 bp), was found to be 52,890 bp. Sixty-one open reading frames (ORFs) having the same transcript orientation were identified. A homology search revealed that 22 of the pHTbeta (pHT) plasmid ORFs showed similarities to the ORFs identified on the pXO2 plasmid (96.2 kbp), which is the virulence plasmid essential for capsule formation by Bacillus anthracis; however, the functions of most of the ORFs remain unknown. Most other ORFs did not show any significant homology to reported genes for which functions have been analyzed. To investigate the highly efficient transfer mechanism of the pHT plasmid, mutations with 174 unique insertions of transposon Tn917-lac insertion mutants of pHTbeta were obtained. Of the 174 derivatives, 92 showed decrease or loss in transfer frequency, and 74 showed normal transfer frequency and LacZ expression. Eight derivatives showed normal transfer and no LacZ expression. Inserts within the 174 derivatives were mapped to 124 different sites on pHTbeta. The Tn917-lac insertions which resulted in altered transfer frequency mapped to three separate regions designated I, II, and III, which were separated by segments in which insertions of Tn917-lac did not affect transfer. There was no region homologous to the previously reported oriT sequences in the pHT plasmid. The oriT was cloned by selection for the ability to mobilize the vector plasmid pAM401. The oriT region resided in a noncoding region (192 bp) between ORF31 and ORF32 and contained three direct repeat sequences and two inverted repeat sequences. ORF34, encoding a 506-amino-acid protein which was located downstream of the oriT region, contains the three conserved motifs (I to III) of the DNA relaxase/nickase of mobile plasmids. The transfer abilities of the Tn917-lac-insertion mutants of ORF34 or a mutant of ORF34 with an in-frame motif III deletion were completely abolished. The sequence of the oriT region and the deduced relaxase/nickase protein of ORF34 showed no significant similarity to the oriT and relaxase/nickase of other conjugative plasmids, respectively. The putative relaxase/nickase protein of ORF34 could be classified as a new member of the MOB(MG) family.

Biocrimes, microbial forensics, and the physician

Schutzer and colleagues give guidance for physicians who believe that one of their patients has been a victim of an act of bioterror or another biocrime.

YrxA is the transcriptional regulator that represses de novo NAD biosynthesis in Bacillus subtilis

The first genetic, in vivo, and in vitro evidences that YrxA is the regulator of NAD de novo biosynthesis in Bacillus subtilis are hereby reported. The protein is essential to the transcription repression of the divergent operons nadBCA and nifS-yrxA in the presence of nicotinic acid and binds to their shared operator-promoter region.

Validation of a pXO2-A PCR assay to explore diversity among Italian isolates of Bacillus anthracis strains closely related to the live, attenuated Carbosap vaccine

Several circulating Bacillus anthracis strains isolated in Italy and belonging to the A1.a cluster, genotype 3 (A1.a-3) are genotypically indistinguishable from Carbosap, a live attenuated vaccine strain, containing both pXO1 and pXO2 plasmids. The genotype was assessed by using eight-locus multilocus variable-number tandem repeat analysis. We describe here the use of a ninth locus able to explore variability among strains that have the same genotype. It is important to be able to genotype the wild isolate of B. anthracis strains from outbreaks of anthrax in areas where Carbosap vaccination of cattle and sheep is common practice. A total of 27 representative field strains isolated in Italy and four vaccinal strains, namely, Carbosap, Sterne, Pasteur I, and Pasteur II, were characterized by a ninth marker, called pXO2-A. Twenty-three field strains were genotype 3 and therefore identical to Carbosap. The marker was in the pXO2 plasmid and is based on the polymorphism of the already-known VX2-3 locus. Detection was obtained by PCR with fluorescence-labeled forward primers in order to produce appropriate fragments for capillary electrophoresis with an ABI 310 genetic analyzer. Genetic relationships showed heterogeneity in all of the examined samples. Interestingly, with respect to genotype 3, samples grouped into eight different subtypes, A to H, and the subtype G, had only two samples indistinguishable from Carbosap. The results of the present study confirm the validity of a hierarchical progressive protocol for discrimination among closely related isolates.

Comparisons of dN/dS are time dependent for closely related bacterial genomes

The ratio of non-synonymous (dN) to synonymous (dS) changes between taxa is frequently computed to assay the strength and direction of selection. Here we note that for comparisons between closely related strains and/or species a second parameter needs to be considered, namely the time since divergence of the two sequences under scrutiny. We demonstrate that a simple time lag model provides a general, parsimonious explanation of the extensive variation in the dN/dS ratio seen when comparing closely related bacterial genomes. We explore this model through simulation and comparative genomics, and suggest a role for hitch-hiking in the accumulation of non-synonymous mutations. We also note taxon-specific differences in the change of dN/dS over time, which may indicate variation in selection, or in population genetics parameters such as population size or the rate of recombination. The effect of comparing intra-species polymorphism and inter-species substitution, and the problems associated with these concepts for asexual prokaryotes, are also discussed. We conclude that, because of the critical effect of time since divergence, inter-taxa comparisons are only possible by comparing trajectories of dN/dS over time and it is not valid to compare taxa on the basis of single time points.

NhaK, a novel monovalent cation/H+ antiporter of Bacillus subtilis

Four Na+/H+ antiporters, Mrp, TetA(L), NhaC, and MleN have so far been described in Bacillus subtilis 168. We identified an additional Na+/H+ antiporter, YvgP, from B. subtilis that exhibits homology to the cation: proton antiporter-1 (CPA-1) family. The yvgP-dependent complementation observed in a Na+(Ca2+)/H+ antiporter-defective Escherichia coli mutant (KNabc) suggested that YvgP effluxed Na+ and Li+. In addition, effects of yvgP expression on a K+ uptake-defective mutant of E. coli indicated that YvgP also supported K+ efflux. In a fluorescence-based assay of everted membrane vesicles prepared from E. coli KNabc transformants, YvgP-dependent Na+ (K+, Li+, Rb+)/H+ antiport activity was demonstrated. Na+ (K+, Li+)/H+ activity was higher at pH 8.5 than at pH 7.5. Mg2+, Ca2+ and Mn2+ did not serve as substrates but they inhibited YvgP antiport activities. Studies of yvgP expression in B. subtilis, using a reporter gene fusion, showed a significant constitutive level of expression that was highest in stationary phase, increasing as stationary phase progressed. In addition, the expression level was significantly increased in the presence of added K+ and Na+.

Adaptation of multilocus sequencing for studying variation within a major clone: evolutionary relationships of Salmonella enterica serovar Typhimurium

Serovar Typhimurium of Salmonella enterica is a model organism for studies of pathogenesis that exhibits phage-type variation and variation in host range and virulence, but in a recent study showed no sequence variation in four genes, indicating the clonal nature of this serovar. We determined the relationships of 46 Typhimurium isolates of nine phage types using mutational changes detected either by matching AFLP (amplified fragment length polymorphism) fragments to computer-modeled LT2 AFLP fragments or by sequencing intergenic regions. Fifty-one polymorphic sites were detected, which gave a single phylogenetic tree. Comparison with genome sequences of five other serovars, Typhi, Paratyphi A, Gallinarum, Enteritidis, and Pullorum, enabled determination of the root of the tree. Only two parallel events were observed, giving high confidence in the tree branching order. The mutation-based tree provided a high level of consistency and a clear lineage for the Typhimurium isolates studied. This enabled us to show that for seven of the nine phage types used, the isolates studied have a single origin, but that two phage types clearly have more than one independent origin. We found that sequencing intergenic regions provides a good strategy for detection of mutational polymorphisms and study of phylogenetic relationships of closely related isolates and would be applicable to many other species.

Essential cysteine residues in Bacillus subtilis spore photoproduct lyase identified by alanine scanning mutagenesis

Endospore-forming bacteria (Bacillus and Clostridium spp.) are highly ultraviolet (UV) resistant and repair UV-induced DNA damage in part using the spore-specific DNA repair enzyme spore photoproduct (SP) lyase. SP lyase in all known sporeformers contains four conserved cysteine residues; three absolutely conserved residues are located at the "Radical SAM" consensus (C91xxxC95xxC98), which presumably participates in [4Fe-4S] cluster formation. A fourth conserved cysteine, the function of which is unknown, is located at C141 in SP lyase from all Bacillus spp. sequenced to date. To probe the function of the fourth cysteine, each conserved cysteine in the B. subtilis SP lyase was systematically altered to alanine by site-directed mutagenesis. UV-visible spectroscopy of wild-type and mutant SP lyases indicated that C141 does not participate in [4Fe-4S] formation and redox chemistry; however, in vivo SP lyase activity was abolished in all mutants, indicating an essential role for C141 in SP lyase activity.

Comparative sequencing of the serine-aspartate repeat-encoding region of the clumping factor B gene (clfB) for resolution within clonal groups of Staphylococcus aureus

Molecular techniques such as spa typing and multilocus sequence typing use DNA sequence data for differentiating Staphylococcus aureus isolates. Although spa typing is capable of detecting both genetic micro- and macrovariation, it has less discriminatory power than the more labor-intensive pulsed-field gel electrophoresis (PFGE) and costly genomic DNA microarray analyses. This limitation hinders strain interrogation for newly emerging clones and outbreak investigations in hospital or community settings where robust clones are endemic. To overcome this constraint, we developed a typing system using DNA sequence analysis of the serine-aspartate (SD) repeat-encoding region within the gene encoding the keratin- and fibrinogen-binding clumping factor B (clfB typing) and tested whether it is capable of discriminating within clonal groups. We analyzed 116 S. aureus strains, and the repeat region was present in all isolates, varying in sequence and in length from 420 to 804 bp. In a sample of 36 well-characterized genetically diverse isolates, clfB typing subdivided identical spa and PFGE clusters which had been discriminated by whole-genome DNA microarray mapping. The combination of spa typing and clfB typing resulted in a discriminatory power (99.5%) substantially higher than that of spa typing alone and closely approached that of the whole-genome microarray (100.0%). clfB typing also successfully resolved genetic differences among isolates differentiated by PFGE that had been collected over short periods of time from single hospitals and that belonged to the most prevalent S. aureus clone in the United States. clfB typing demonstrated in vivo, in vitro, and interpatient transmission stability yet revealed that this locus may be recombinogenic in a primarily clonal population structure. Taken together, these data show that the SD repeat-encoding region of clfB is a highly stable marker of microvariation, that in conjunction with spa typing it may serve as a DNA sequence-based alternative to PFGE for investigating genetically similar strains, and that it is useful for analyzing collections of isolates in both long-term population-based and local epidemiologic studies.

A comparative study ofBacillus cereus,Bacillus thuringiensis andBacillus anthracis extracellular proteomes

Bacillus cereus, Bacillus thuringiensis and Bacillus anthracis are closely related species that share a similar genetic background but occupy different ecological niches. Virulence plasmids bearing genes coding for toxins, may explain, at least partly, this specialization. We have compared by 2-DE in the early stationary phase of growth the extracellular proteomes of three strains of these species that have lost their virulence plasmids. Proteins expected to be secreted or to belong to the cell wall or to the cytosol were found in the three proteomes. For the cell wall and cytosolic proteins located in the extracellular space, the three proteomes were similar. Cytosolic proteins included enolase, GroEL, PdhB, PdhD, SodA and others. Cell surface proteins were mainly autolysins, proteases, nucleotidases and OppAs. In contrast, the secreted proteins profiles of B. cereus and B. thuringiensis were quite different from that of B. anthracis. B. cereus and B. thuringiensis extracellular proteomes both contained large amounts of secreted degradative enzymes and toxins, including nine proteases, three phospholipases, two haemolysins and several enterotoxins. Most of the genes encoding these enzymes and toxins are controlled by the transcriptional activator PlcR. The extracellular proteome of the pXO1-, pXO2- B. anthracis 9131 strain contained only one secreted protein: the metalloprotease InhA1, also found in the proteomes of the two other strains and possibly involved in antibacterial peptide degradation.

Comparative genome analysis ofBacillus cereusgroup genomes withBacillus subtilis

Genome features of the Bacillus cereus group genomes (representative strains of Bacillus cereus, Bacillus anthracis and Bacillus thuringiensis sub spp. israelensis) were analyzed and compared with the Bacillus subtilis genome. A core set of 1381 protein families among the four Bacillus genomes, with an additional set of 933 families common to the B. cereus group, was identified. Differences in signal transduction pathways, membrane transporters, cell surface structures, cell wall, and S-layer proteins suggesting differences in their phenotype were identified. The B. cereus group has signal transduction systems including a tyrosine kinase related to two-component system histidine kinases from B. subtilis. A model for regulation of the stress responsive sigma factor sigmaB in the B. cereus group different from the well studied regulation in B. subtilis has been proposed. Despite a high degree of chromosomal synteny among these genomes, significant differences in cell wall and spore coat proteins that contribute to the survival and adaptation in specific hosts has been identified.

Molecular diversity of Bacillus anthracis in Italy

We used multiple-locus variable-number tandem repeat analysis (MLVA) to type 64 Bacillus anthracis isolates from outbreaks that have occurred during the past 40 years in Italy. MLVA of the 64 isolates revealed 10 unique genotypes; 9 of these genotypes and the majority of isolates (63/64) belonged to the previously described genetic cluster A1.a. Within the A1.a isolates, two previously described genotypes (G1 and G3), which differ by a single mutation in the pX01 locus, account for the majority of isolates in the country (53/63). The low diversity of B. anthracis genotypes in Italy suggests a single, dominant historical introduction, followed by limited localized differentiation.

Unusual group II introns in bacteria of the Bacillus cereus group

A combination of sequence and structure analysis and reverse transcriptase PCR experiments was used to characterize the group II introns in the complete genomes of two strains of the pathogen Bacillus cereus. While B. cereus ATCC 14579 harbors a single intron element in the chromosome, B. cereus ATCC 10987 contains three introns in the chromosome and four in its 208-kb pBc10987 plasmid. The most striking finding is the presence in B. cereus ATCC 10987 of an intron [B.c.I2(a)] located on the reverse strand of a gene encoding a putative cell surface protein which appears to be correlated to strains of clinical origin. Because of the opposite orientation of B.c.I2(a), the gene is disrupted. Even more striking is that B.c.I2(a) splices out of an RNA transcript corresponding to the opposite DNA strand. All other intragenic introns studied here are inserted in the same orientation as their host genes and splice out of the mRNA in vivo, setting the flanking exons in frame. Noticeably, B.c.I3 in B. cereus ATCC 10987 represents the first example of a group II intron entirely included within a conserved replication gene, namely, the alpha subunit of DNA polymerase III. Another striking finding is that the observed 3' splice site of B.c.I4 occurs 56 bp after the predicted end of the intron. This apparently unusual splicing mechanism may be related to structural irregularities in the 3' terminus. Finally, we also show that the intergenic introns of B. cereus ATCC 10987 are transcribed with their upstream genes and do splice in vivo.

Amplified intergenic locus polymorphism as a basis for bacterial typing of Listeria spp. and Escherichia coli

DNA-based methods are increasingly important for bacterial typing. The high number of polymorphic sites present among closely related bacterial genomes is the basis for the presented method. The method identifies multilocus genomic polymorphisms in intergenic regions termed AILP (amplified intergenic locus polymorphism). For each locus, a pair of unique PCR primers was designed to amplify an intergenic sequence from one open reading frame (ORF) to the adjacent ORF. Presence, absence, and size variation of the amplification products were identified and used as genetic markers for rapidly differentiating among strains. Polymorphism was evaluated using 18 AILP sites among 28 strains of Listeria monocytogenes and 6 strains of Listeria spp. and 30 AILP markers among 27 strains of Escherichia coli. Up to four alleles per locus were identified among Listeria strains, and up to six were identified among E. coli strains. In both species, more than half of the AILP sites revealed intraspecies polymorphism. The AILP data were applied to phylogenetic analysis among Listeria and E. coli strains. A clear distinction between L. monocytogenes and Listeria spp. was demonstrated. In addition, the method separated L. monocytogenes into the three known lineages and discriminated the most common virulent serotypic group, 4b. In E. coli, AILP analysis separated the known groups as well as the virulent O157:H7 isolates. These findings for both Listeria and E. coli are in agreement with other phylogenetic studies using molecular markers. The AILP method was found to be rapid, simple, reproducible, and a low-cost method for initial bacterial typing that could serve as a basis for epidemiological investigation.

Molecular characterization of plasmid pBM300 from Bacillus megaterium QM B1551

Strain QM B1551 of Bacillus megaterium contains seven compatible plasmids: two small rolling circle plasmids and five theta-replicating plasmids with cross-hybridizing replicons. To expand our understanding of these plasmids, the replicon region (6.7 kb) from pBM300 was cloned, sequenced, and functionally characterized. Sequence analysis showed that the replication protein (RepM300) was highly homologous to two other plasmid Rep proteins of the same strain but to no other known proteins. Furthermore, the location of the replication origin was within the RepM300 coding region, and the origin contained three 12-base direct repeats. Deletion analysis of the replicon confirmed the role of the Rep protein and showed that open reading frame 2 (ORF2) was required for stability. However, the protein encoded by ORF2 is entirely different from the replicon stability proteins encoded by the other two replicons. The entire plasmid was isolated from the plasmid array by integrating a spectinomycin resistance gene and transforming a plasmidless strain, PV361. Complete sequencing showed that pBM300 was 26,300 bp long, had a G+C content of 35.2%, and contained 20 ORFs, two of which encoded proteins that had no similarity to other proteins in the database. The proteins encoded by the plasmid ORFs had similarity to proteins for mobilization and transfer, an integrase, a rifampin resistance protein, a cell wall hydrolase, glutathione synthase, and a biotin carboxylase. The similarities were to several gram-positive genera and a few gram-negative genera and archaea. oriT and ssoT-like regions were detected near two mob genes. These results suggest that pBM300 is a mobilizable hybrid plasmid that confers increased metabolic and germination ability on its host. Its replicon also helps define a new plasmid family.

Proteomics reveals that proteins expressed during the early stage of Bacillus anthracis infection are potential targets for the development of vaccines and drugs

In this review, we advance a new concept in developing vaccines and/or drugs to target specific proteins expressed during the early stage of Bacillus anthracis (anthrax) infection and address existing challenges to this concept. Three proteins (immune inhibitor A, GPR-like spore protease, and alanine racemase) initially identified by proteomics in our laboratory were found to have differential expressions during anthrax spore germination and early outgrowth. Other studies of different bacillus strains indicate that these three proteins are involved in either germination or cytotoxicity of spores, suggesting that they may serve as potential targets for the design of anti-anthrax vaccines and drugs.

Comparative secretome analyses of three Bacillus anthracis strains with variant plasmid contents

Bacillus anthracis, the causative agent of anthrax, secretes numerous proteins into the extracellular environment during infection. A comparative proteomic approach was employed to elucidate the differences among the extracellular proteomes (secretomes) of three isogenic strains of B. anthracis that differed solely in their plasmid contents. The strains utilized were the wild-type virulent B. anthracis RA3 (pXO1(+) pXO2(+)) and its two nonpathogenic derivative strains: the toxigenic, nonencapsulated RA3R (pXO1(+) pXO2(-)) and the totally cured, nontoxigenic, nonencapsulated RA3:00 (pXO1(-) pXO2(-)). Comparative proteomics using two-dimensional gel electrophoresis followed by computer-assisted gel image analysis was performed to reveal unique, up-regulated, or down-regulated secretome proteins among the strains. In total, 57 protein spots, representing 26 different proteins encoded on the chromosome or pXO1, were identified by peptide mass fingerprinting. S-layer-derived proteins, such as Sap and EA1, were most frequently observed. Many sporulation-associated enzymes were found to be overexpressed in strains containing pXO1(+). This study also provides evidence that pXO2 is necessary for the maximal expression of the pXO1-encoded toxins lethal factor (LF), edema factor (EF), and protective antigen (PA). Several newly identified putative virulence factors were observed; these include enolase, a high-affinity zinc uptake transporter, the peroxide stress-related alkyl hydroperoxide reductase, isocitrate lyase, and the cell surface protein A.

Chips and SNPs, bugs and thugs: a molecular sleuthing perspective

Recent events both here and abroad have focused attention on the need for ensuring a safe and secure food supply. Although much has been written about the potential of particular select agents in bioterrorism, we must consider seriously the more mundane pathogens, especially those that have been implicated previously in foodborne outbreaks of human disease, as possible agents of bioterrorism. Given their evolutionary history, the enteric pathogens are more diverse than agents such as Bacillus anthracis, Francisella tularensis, or Yersinia pestis. This greater diversity, however, is a double-edged sword; although diversity affords the opportunity for unequivocal identification of an organism without the need for whole-genome sequencing, the same diversity can confound definitive forensic identification if boundaries are not well defined. Here, we discuss molecular approaches used for the identification of Salmonella enterica, Escherichia coli, and Shigella spp. and viral pathogens and discuss the utility of these approaches to the field of microbial molecular forensics.

Efflux-mediated antimicrobial resistance

Antibiotic resistance continues to plague antimicrobial chemotherapy of infectious disease. And while true biocide resistance is as yet unrealized, in vitro and in vivo episodes of reduced biocide susceptibility are common and the history of antibiotic resistance should not be ignored in the development and use of biocidal agents. Efflux mechanisms of resistance, both drug specific and multidrug, are important determinants of intrinsic and/or acquired resistance to these antimicrobials, with some accommodating both antibiotics and biocides. This latter raises the spectre (as yet generally unrealized) of biocide selection of multiple antibiotic-resistant organisms. Multidrug efflux mechanisms are broadly conserved in bacteria, are almost invariably chromosome-encoded and their expression in many instances results from mutations in regulatory genes. In contrast, drug-specific efflux mechanisms are generally encoded by plasmids and/or other mobile genetic elements (transposons, integrons) that carry additional resistance genes, and so their ready acquisition is compounded by their association with multidrug resistance. While there is some support for the latter efflux systems arising from efflux determinants of self-protection in antibiotic-producing Streptomyces spp. and, thus, intended as drug exporters, increasingly, chromosomal multidrug efflux determinants, at least in Gram-negative bacteria, appear not to be intended as drug exporters but as exporters with, perhaps, a variety of other roles in bacterial cells. Still, given the clinical significance of multidrug (and drug-specific) exporters, efflux must be considered in formulating strategies/approaches to treating drug-resistant infections, both in the development of new agents, for example, less impacted by efflux and in targeting efflux directly with efflux inhibitors.

Virulence signatures: microarray-based approaches to discovery and analysis

Rapid, accurate, and sensitive detection of biothreat agents requires a broad-spectrum assay capable of discriminating between closely related microbial or viral pathogens. Moreover, in cases where a biological agent release has been identified, forensic analysis demands detailed genetic signature data for accurate strain identification and attribution. To date, nucleic acid sequences have provided the most robust and phylogentically illuminating signature information. Nucleic acid signature sequences are not often linked to genomic or extrachromosomal determinants of virulence, a link that would further facilitate discrimination between pathogens and closely related species. Inextricably coupling genetic determinants of virulence with highly informative nucleic acid signatures would provide a robust means of identifying human, livestock, and agricultural pathogens. By means of example, we present here an overview of two general applications of microarray-based methods for: (1) the identification of candidate virulence factors; and (2) the analysis of genetic polymorphisms that are coupled to Bacillus anthracis virulence factors using an accurate, low cost solid-phase mini-sequencing assay. We show that microarray-based analysis of gene expression can identify potential virulence associated genes for use as candidate signature targets, and, further, that microarray-based single nucleotide polymorphism assays provide a robust platform for the detection and identification of signature sequences in a manner independent of the genetic background in which the signature is embedded. We discuss the strategy as a general approach or pipeline for the discovery of virulence-linked nucleic acid signatures for biothreat agents.

A new Bacillus cereus DNA-binding protein, HlyIIR, negatively regulates expression of B. cereus haemolysin II

Haemolysin II, HlyII, is one of several cytotoxic proteins produced by Bacillus cereus, an opportunistic human pathogen that causes food poisoning. The hlyII gene confers haemolytic activity to Escherichia coli cells. Here a new B. cereus gene, hlyIIR, which is located immediately downstream of hlyII and regulates hlyII expression, is reported. The deduced amino acid sequence of HlyIIR is similar to prokaryotic DNA-binding transcriptional regulators of the TetR/AcrA family. Measurements of haemolytic activity levels and of hlyII promoter activity levels using gene fusions and primer-extension assays demonstrated that, in E. coli, hlyII transcription decreased in the presence of hlyIIR. Recombinant HlyIIR binds to a 22 bp inverted DNA repeat centred 48 bp upstream of the hlyII promoter transcription initiation point. In vitro transcription studies showed that HlyIIR inhibits transcription from the hlyII promoter by binding to the 22 bp repeat and RNA polymerase, and by decreasing the formation of the catalytically competent open promoter complex.

Microbial forensics: the next forensic challenge

Pathogens and toxins can be converted to bioweapons and used to commit bioterrorism and biocrime. Because of the potential and relative ease of an attack using a bioweapon, forensic science needs to be prepared to assist in the investigation to bring perpetrators to justice and to deter future attacks. A new subfield of forensics--microbial forensics--has been created, which is focused on characterization of evidence from a bioterrorism act, biocrime, hoax, or an inadvertent release. Forensic microbiological investigations are essentially the same as any other forensic investigation regarding processing. They involve crime scene(s) investigation, chain of custody practices, evidence collection, handling and preservation, evidence shipping, analysis of evidence, interpretation of results, and court presentation. In addition to collecting and analyzing traditional forensic evidence, the forensic investigation will attempt to determine the etiology and identity of the causal agent, often in a similar fashion as in an epidemiologic investigation. However, for attribution, higher-resolution characterization is needed. The tools for attribution include genetic- and nongenetic-based assays and informatics to attempt to determine the unique source of a sample or at least eliminate some sources. In addition, chemical and physical assays may help determine the process used to prepare, store, or disseminate the bioweapon. An effective microbial forensics program will require development and/or validation of all aspects of the forensic investigative process, from sample collection to interpretation of results. Quality assurance (QA) and QC practices, comparable to those used by the forensic DNA science community, are being implemented. Lastly, partnerships with other laboratories will be requisite, because many of the necessary capabilities for analysis will not reside in the traditional forensic laboratory.

Time dependency of molecular rate estimates and systematic overestimation of recent divergence times

Studies of molecular evolutionary rates have yielded a wide range of rate estimates for various genes and taxa. Recent studies based on population-level and pedigree data have produced remarkably high estimates of mutation rate, which strongly contrast with substitution rates inferred in phylogenetic (species-level) studies. Using Bayesian analysis with a relaxed-clock model, we estimated rates for three groups of mitochondrial data: avian protein-coding genes, primate protein-coding genes, and primate d-loop sequences. In all three cases, we found a measurable transition between the high, short-term (< 1-2 Myr) mutation rate and the low, long-term substitution rate. The relationship between the age of the calibration and the rate of change can be described by a vertically translated exponential decay curve, which may be used for correcting molecular date estimates. The phylogenetic substitution rates in mitochondria are approximately 0.5% per million years for avian protein-coding sequences and 1.5% per million years for primate protein-coding and d-loop sequences. Further analyses showed that purifying selection offers the most convincing explanation for the observed relationship between the estimated rate and the depth of the calibration. We rule out the possibility that it is a spurious result arising from sequence errors, and find it unlikely that the apparent decline in rates over time is caused by mutational saturation. Using a rate curve estimated from the d-loop data, several dates for last common ancestors were calculated: modern humans and Neandertals (354 ka; 222-705 ka), Neandertals (108 ka; 70-156 ka), and modern humans (76 ka; 47-110 ka). If the rate curve for a particular taxonomic group can be accurately estimated, it can be a useful tool for correcting divergence date estimates by taking the rate decay into account. Our results show that it is invalid to extrapolate molecular rates of change across different evolutionary timescales, which has important consequences for studies of populations, domestication, conservation genetics, and human evolution.

SNPsFinder--a web-based application for genome-wide discovery of single nucleotide polymorphisms in microbial genomes

Single nucleotide polymorphisms (SNPs) are the most abundant form of genetic variations in closely related microbial species, strains or isolates. Some SNPs confer selective advantages for microbial pathogens during infection and many others are powerful genetic markers for distinguishing closely related strains or isolates that could not be distinguished otherwise. To facilitate SNP discovery in microbial genomes, we have developed a web-based application, SNPsFinder, for genome-wide identification of SNPs. SNPsFinder takes multiple genome sequences as input to identify SNPs within homologous regions. It can also take contig sequences and sequence quality scores from ongoing sequencing projects for SNP prediction. SNPsFinder will use genome sequence annotation if available and map the predicted SNP regions to known genes or regions to assist further evaluation of the predicted SNPs for their functional significance. SNPsFinder can generate PCR primers for all predicted SNP regions according to user's input parameters to facilitate experimental validation. The results from SNPsFinder analysis are accessible through the World Wide Web.

AVAILABILITY

The SNPsFinder program is available at http://snpsfinder.lanl.gov/.

SUPPLEMENTARY INFORMATION

The user's manual is available at http://snpsfinder.lanl.gov/UsersManual/

Real-time PCR assay for a unique chromosomal sequence of Bacillus anthracis

Real-time PCR has become an important method for the rapid identification of Bacillus anthracis since the 2001 anthrax mailings. Most real-time PCR assays for B. anthracis have been developed to detect virulence genes located on the pXO1 and pXO2 plasmids. In contrast, only two published chromosomal targets exist, the rpoB gene and the gyrA gene. In the present study, subtraction-hybridization with a plasmid-cured B. anthracis tester strain and a Bacillus cereus driver was used to find a unique chromosomal sequence. By targeting this region, a real-time assay was developed with the Ruggedized Advanced Pathogen Identification Device. Further testing has revealed that the assay has 100% sensitivity and 100% specificity, with a limit of detection of 50 fg of DNA. The results of a search for sequences with homology with the BLAST program demonstrated significant alignment to the recently published B. anthracis Ames strain, while an inquiry for protein sequence similarities indicated homology with an abhydrolase from B. anthracis strain A2012. The importance of this chromosomal assay will be to verify the presence of B. anthracis independently of plasmid occurrence.

Sequence, Transcription Activity, and Evolutionary Origin of the R-BodyCoding Plasmid pKAP298 from the Intracellular Parasitic BacteriumCaedibacter taeniospiralis

We isolated the intracellular parasitic bacterium Caedibacter taeniospiralis from cultures of the freshwater ciliate Paramecium tetraurelia strain 298. Plasmid pKAP298 as well as the total RNA were isolated from the bacteria. pKAP298 was totally sequenced (49.1 kb; NCBI accession number AY422720). From southern blots of pKAP-fragments and Digoxigenin-labeled cDNA of the Caedibacter-RNA, we generated transcription maps of pKAP298. The observed transcription activity indicated functions of the plasmid besides the synthesis of the R-body, a complex protein inclusion associated with toxic effects of Caedibacter cells on host paramecia. We identified 63 potential protein coding regions on pKAP298, and a novel transposon as well as known transposons were characterized. A group II intron was identified. Homologies with putative phage genes were detected on pKAP298 that direct to the evolution of pKAP298 from a bacteriophage. This original phage most probably belonged to the Caudovirales. Hints on a toxin coding region of pKAP298 are given: a protein with homology to the Soj-/ParA-family also has homologies to a membrane associated ATPase, which is involved in eukaryotic ATPase dependent ion carriers and may be associated with toxic effects on paramecia ingesting this protein.

IS231-MIC231 elements from Bacillus cereus sensu lato are modular

Summary IS231A was originally discovered in Bacillus thuringiensis as a typical 1.6 kb insertion sequence (IS) displaying 20 bp inverted repeats (IR) flanking a transposase gene. A first major variation of this canonical organization was found in MIC231A1. This mobile insertion cassette (MIC), delineated by IS231A-related extremities, contained an active d-stereospecific endopeptidase (adp) gene instead of a transposase. Interestingly, it was shown that MIC231A1 can be mobilized in trans by the IS231A transposase. In this paper, we show that this family of IS231-MIC231 elements can be extended to a broad range of related entities displaying higher levels of structural complexity. Several IS231A-like elements contained, upstream of their transposase gene, passenger genes coding for putative antibiotic resistances or regulatory factors. Furthermore, the diversity of the MIC231 elements ranged from empty cassettes to structures carrying up to three passenger genes. Among these, MIC231V carried, in addition to an adp gene, an active fosfomycin resistance determinant. In vivo transposition assays showed that MIC231V is also trans-activated by the IS231A transposase. These results lend further support to the potential contribution of these modular mobile elements to the genome plasticity of the Bacillus cereus/B. thuringiensis group.

Biosynthesis of the angiogenesis inhibitor borrelidin by Streptomyces parvulus Tü4055: insights into nitrile formation

The 18-membered polyketide macrolide borrelidin exhibits a number of important biological activities, including potent angiogenesis inhibition. This has prompted two recent total syntheses as well as the cloning of the biosynthetic gene cluster from Streptomyces parvulus Tü4055. Borrelidin possesses some unusual structural characteristics, including a cyclopentane carboxylic acid moiety at C17 and a nitrile moiety at C12 of the macrocyclic ring. Nitrile groups are relatively rare in nature, and little is known of their biosynthesis during secondary metabolism. The nitrile group of borrelidin is shown here to arise from the methyl group of a methylmalonyl-CoA extender unit incorporated during polyketide chain extension. Insertional inactivation of two genes in the borrelidin gene cluster, borI (coding for a cytochrome P450 monooxygenase) and borJ (coding for an aminotransferase), generated borrelidin non-producing mutants. These mutants accumulated different compounds lacking the C12 nitrile moiety, with the product of the borI-minus mutant (12-desnitrile-12-methyl-borrelidin) possessing a methyl group and that of the borJ-minus mutant (12-desnitrile-12-carboxyl-borrelidin) a carboxyl group at C12. The former but not the latter was converted into borrelidin when biotransformed by an S. parvulus mutant that is deficient in the biosynthesis of the borrelidin starter unit. This suggests that 12-desnitrile-12-methyl-borrelidin is a competent biosynthetic intermediate, whereas the carboxylated derivative is a shunt metabolite. Bioconversion of 12-desnitrile-12-methyl-borrelidin into borrelidin was also achieved in a heterologous system co-expressing borI and borJ in Streptomyces albus J1074. This bioconversion was more efficient when borK, which is believed to encode a dehydrogenase, was simultaneously expressed with borI and borJ. On the basis of these findings, a pathway is proposed for the formation of the nitrile moiety during borrelidin biosynthesis.

Functional cloning of Bacillus anthracis dihydrofolate reductase and confirmation of natural resistance to trimethoprim

Bacillus anthracis is reported to be naturally resistant to trimethoprim (TMP), a drug that inhibits dihydrofolate reductase (DHFR), a key enzyme in the folate pathway. A microdilution broth assay established that the MIC of TMP for B. anthracis Sterne is >2,048 but < or =4,096 microg/ml. A putative DHFR sequence was amplified from B. anthracis Sterne genomic DNA. The PCR product was cloned into the Invitrogen pCRT7/CT-TOPO vector, followed by transformation into Escherichia coli TOP10F' chemically competent cells. Plasmid DNA from a clone showing the correct construct with a thrombin cleavage site attached downstream from the terminus of the cloned PCR product was transformed into E. coli BL21 Star (DE3)pLysS competent cells for expression of the six-histidine-tagged fusion protein and purification on a His-Bind resin column. Functionality of the purified Sterne recombinant DHFR (Sterne rDHFR) was confirmed in an established enzyme assay. The 50% inhibitory concentrations of TMP and methotrexate for the Sterne rDHFR were found to be 77,233 and 12.2 nM, respectively. TMP resistance was observed with E. coli BL21 Star (DE3)pLysS competent cells transformed with the Sterne DHFR gene. Alignment of the amino acid sequence of the Sterne DHFR gene revealed 100% homology with various virulent strains of B. anthracis. These results confirm the natural resistance of B. anthracis to TMP and clarify that the resistance is correlated to a lack of selectivity for the chromosomally encoded gene product. These findings will assist in the development of narrow-spectrum antimicrobial agents for treatment of anthrax.

Microarray-based resequencing of multiple Bacillus anthracis isolates

Custom-designed resequencing arrays were used to generate 3.1 Mb of genomic sequence from a panel of 56 Bacillus anthracis strains. Sequence quality was shown to be very high by replication and by comparison to independently generated shotgun sequence

We used custom-designed resequencing arrays to generate 3.1 Mb of genomic sequence from a panel of 56 Bacillus anthracis strains. Sequence quality was shown to be very high by replication (discrepancy rate of 7.4 × 10^-7) and by comparison to independently generated shotgun sequence (discrepancy rate < 2.5 × 10^-6). Population genomics studies of microbial pathogens using rapid resequencing technologies such as resequencing arrays are critical for recognizing newly emerging or genetically engineered strains.

The bcr1 DNA repeat element is specific to the Bacillus cereus group and exhibits mobile element characteristics

Bacillus cereus strains ATCC 10987 and ATCC 14579 harbor an approximately 155-bp repeated element, bcr1, which is conserved in B. cereus, B. anthracis, B. thuringiensis, and B. mycoides but not in B. subtilis and B. licheniformis. In this study, we show by Southern blot hybridizations that bcr1 is present in all 54 B. cereus group strains tested but absent in 11 Bacillus strains outside the group, suggesting that bcr1 may be specific and ubiquitous to the B. cereus group. By comparative analysis of the complete genome sequences of B. cereus ATCC 10987, B. cereus ATCC 14579, and B. anthracis Ames, we show that bcr1 is exclusively present in the chromosome but absent from large plasmids carried by these strains and that the numbers of full-length bcr1 repeats for these strains are 79, 54, and 12, respectively. Numerous copies of partial bcr1 elements are also present in the three genomes (91, 128, and 53, respectively). Furthermore, the genomic localization of bcr1 is not conserved between strains with respect to chromosomal position or organization of gene neighbors, as only six full-length bcr1 loci are common to at least two of the three strains. However, the intergenic sequence surrounding a specific bcr1 repeat in one of the three strains is generally strongly conserved in the other two, even in loci where bcr1 is found exclusively in one strain. This finding indicates that bcr1 either has evolved by differential deletion from a very high number of repeats in a common ancestor to the B. cereus group or is moving around the chromosome. The identification of bcr1 repeats interrupting genes in B. cereus ATCC 10987 and ATCC 14579 and the presence of a flanking TTTAT motif in each end show that bcr1 exhibits features characteristic of a mobile element.

Identification and inactivation of genetic loci involved with Lactobacillus acidophilus acid tolerance

Amino acid decarboxylation-antiporter reactions are one of the most important systems for maintaining intracellular pH between physiological limits under acid stress. We analyzed the Lactobacillus acidophilus NCFM complete genome sequence and selected four open reading frames with similarities to genes involved with decarboxylation reactions involved in acid tolerance in several microorganisms. Putative genes encoding an ornithine decarboxylase, an amino acid permease, a glutamate gamma-aminobutyrate antiporter, and a transcriptional regulator were disrupted by insertional inactivation. The ability of L. acidophilus to survive low-pH conditions, such as those encountered in the stomach or fermented dairy foods, was investigated and compared to the abilities of early- and late-stationary-phase cells of the mutants by challenging them with a variety of acidic conditions. All of the integrants were more sensitive to low pH than the parental strain. Interestingly, each integrant also exhibited an adaptive acid response during logarithmic growth, indicating that multiple mechanisms are present and orchestrated in L. acidophilus in response to acid challenge.

Sulfatases: Structure, Mechanism, Biological Activity, Inhibition, and Synthetic Utility

Sulfatases, which cleave sulfate esters in biological systems, play a key role in regulating the sulfation states that determine the function of many physiological molecules. Sulfatase substrates range from small cytosolic steroids, such as estrogen sulfate, to complex cell-surface carbohydrates, such as the glycosaminoglycans. The transformation of these molecules has been linked with important cellular functions, including hormone regulation, cellular degradation, and modulation of signaling pathways. Sulfatases have also been implicated in the onset of various pathophysiological conditions, including hormone-dependent cancers, lysosomal storage disorders, developmental abnormalities, and bacterial pathogenesis. These findings have increased interest in sulfatases and in targeting them for therapeutic endeavors. Although numerous sulfatases have been identified, the wide scope of their biological activity is only beginning to emerge. Herein, accounts of the diversity and growing biological relevance of sulfatases are provided along with an overview of the current understanding of sulfatase structure, mechanism, and inhibition.