Identification of Bacillus anthracis specific chromosomal sequences by suppressive subtractive hybridization

New bacteriophage-derived lysins, LysJ and LysF, with the potential to control Bacillus anthracis

Bacillus anthracis is an etiological agent of anthrax, a severe zoonotic disease that can be transmitted to people and cause high mortalities. Bacteriophages and their lytic enzymes, endolysins, have potential therapeutic value in treating infections caused by this bacterium as alternatives or complements to antibiotic therapy. They can also be used to identify and detect B. anthracis. Endolysins of two B. anthracis Wbetavirus phages, J5a and F16Ba which were described by us recently, differ significantly from the best-known B. anthracis phage endolysin PlyG from Wbetavirus genus bacteriophage Gamma and a few other Wbetavirus genus phages. They are larger than PlyG (351 vs. 233 amino acid residues), contain a signal peptide at their N-termini, and, by prediction, have a different fold of cell binding domain suggesting different structural basis of cell epitope recognition. We purified in a soluble form the modified versions of these endolysins, designated by us LysJ and LysF, respectively, and depleted of signal peptides. Both modified endolysins could lyse the B. anthracis cell wall in zymogram assays. Their activity against the living cells of B. anthracis and other species of Bacillus genus was tested by spotting on the layers of bacteria in soft agar and by assessing the reduction of optical density of bacterial suspensions. Both methods proved the effectiveness of LysJ and LysF in killing the anthrax bacilli, although the results obtained by each method differed. Additionally, the lytic efficiency of both proteins was different, which apparently correlates with differences in their amino acid sequence. KEY POINTS: • LysJ and LysF are B. anthracis-targeting lysins differing from lysins studied so far • LysJ and LysF could be overproduced in E. coli in soluble and active forms • LysJ and LysF are active in killing cells of B. anthracis virulent strains.

Rapid identification of Bacillus anthracis by real-time PCR with dual hybridization probes in environmental swabs

In the present study, we report the development of a real-time PCR assay for the identification of Bacillus anthracis, based on the amplification of a unique chromosomal marker, the E4 sequence, with dual hybridization probes. The assay was evaluated using a panel of ten B. anthracis strains, two B. anthracis isolates from human clinical samples, 12 B. anthracis environmental swabs and 40 non- B. anthracis strains. All 12 B. anthracis strains and clinical isolates were correctly detected, and the method did not show cross-reactions with other micro-organisms. Likewise, the E4 sequence was not found in those strains of B. thuringiensis and B. cereus closely related (homology > 90%) to B. anthracis by computer analysis. On the other hand, this molecular assay showed a high analytical sensitivity, 3.5 genome equivalents per reaction at 95% probability. Furthermore, the real-time PCR assay allowed sequence-specific detection of the amplicon (melting peak with a Tm of 63.5 °C ± 0.5 °C) without post-amplification procedures, which offers an additional advantage over other qPCR assays for B. anthracis detection. Finally, the performance of the method was successfully evaluated in 12 environmental samples. In summary, we have developed a rapid and specific method for the molecular identification of Bacillus anthracis in environmental samples.

Divergence of protein-coding capacity and regulation in the Bacillus cereus sensu lato group

BACKGROUND

The Bacillus cereus sensu lato group contains ubiquitous facultative anaerobic soil-borne Gram-positive spore-forming bacilli. Molecular phylogeny and comparative genome sequencing have suggested that these organisms should be classified as a single species. While clonal in nature, there do not appear to be species-specific clonal lineages, excepting B. anthracis, in spite of the wide array of phenotypes displayed by these organisms.

RESULTS

We compared the protein-coding content of 201 B. cereus sensu lato genomes to characterize differences and understand the consequences of these differences on biological function. From this larger group we selected a subset consisting of 25 whole genomes for deeper analysis. Cluster analysis of orthologous proteins grouped these genomes into five distinct clades. Each clade could be characterized by unique genes shared among the group, with consequences for the phenotype of each clade. Surprisingly, this population structure recapitulates our recent observations on the divergence of the generalized stress response (SigB) regulons in these organisms. Divergence of the SigB regulon among these organisms is primarily due to the placement of SigB-dependent promoters that bring genes from a common gene pool into/out of the SigB regulon.

CONCLUSIONS

Collectively, our observations suggest the hypothesis that the evolution of these closely related bacteria is a consequence of two distinct processes. Horizontal gene transfer, gene duplication/divergence and deletion dictate the underlying coding capacity in these genomes. Regulatory divergence overlays this protein coding reservoir and shapes the expression of both the unique and shared coding capacity of these organisms, resulting in phenotypic divergence. Data from other organisms suggests that this is likely a common pattern in prokaryotic evolution.

Phages preying on Bacillus anthracis, Bacillus cereus, and Bacillus thuringiensis: past, present and future

Many bacteriophages (phages) have been widely studied due to their major role in virulence evolution of bacterial pathogens. However, less attention has been paid to phages preying on bacteria from the Bacillus cereus group and their contribution to the bacterial genetic pool has been disregarded. Therefore, this review brings together the main information for the B. cereus group phages, from their discovery to their modern biotechnological applications. A special focus is given to phages infecting Bacillus anthracis, B. cereus and Bacillus thuringiensis. These phages belong to the Myoviridae, Siphoviridae, Podoviridae and Tectiviridae families. For the sake of clarity, several phage categories have been made according to significant characteristics such as lifestyles and lysogenic states. The main categories comprise the transducing phages, phages with a chromosomal or plasmidial prophage state, γ-like phages and jumbo-phages. The current genomic characterization of some of these phages is also addressed throughout this work and some promising applications are discussed here.

Existence of separate domains in lysin PlyG for recognizing Bacillus anthracis spores and vegetative cells

As a potential antimicrobial, the bacteriophage lysin PlyG has been reported to specifically recognize Bacillus anthracis vegetative cells only and to kill B. anthracis vegetative cells and its germinating spores. However, how PlyG interacts with B. anthracis spores remains unclear. Herein, a 60-amino-acid domain in PlyG (residues 106 to 165), located mainly in the previously identified catalytic domain, was found able to specifically recognize B. anthracis spores but not vegetative cells. The exosporium of the spores was found to be the most probable binding target of this domain. This is the first time that a lysin for spore-forming bacteria has been found to have separate domains to recognize spores and vegetative cells, which might help in understanding the coevolution of phages with spore-forming bacteria. Besides providing new biomarkers for developing better assays for identifying B. anthracis spores, the newly found domain may be helpful in developing PlyG as a preventive antibiotic to reduce the threat of anthrax in suspected exposures to B. anthracis spores.

Developing an integrated proteo-genomic approach for the characterisation of biomarkers for the identification of Bacillus anthracis

Bacillus anthracis is the causative agent of anthrax, an acute and often fatal disease in humans. Due to the high genomic relatedness within the Bacillus cereus group of species it is a challenge to identify B. anthracis consistently. Alternative strategies such as proteomics coupled with mass spectrometry (MS) provide a powerful approach for biomarker discovery. However, validating and evaluating these markers, particularly for genetically homogeneous species such as B. anthracis are challenging. The objective of this study is to develop a robust biomarker discovery and validation pipeline, using proteomic methodology combined with in silico and molecular approaches, to determine a biomarker list, using B. anthracis as a model. In this exploratory study we profiled the proteome of B. anthracis and genetically related species using GeLC-Liquid Chromatography MS/MS (GeLC-LC MS/MS), identifying peptides that could be used to detect B. anthracis. Peptides were filtered to remove low quality identifications. Using comparative bioinformatic approaches, matching and searching against genomic sequence data a shortlist of peptide biomarkers was determined and validated using DNA sequencing, against a panel of closely related strains, to determine marker specificity. Further validation was performed using MS quantitation methods to assess sensitivity and specificity. A biomarker discovery pipeline was successfully developed in this study, comprising four distinct stages: proteome profiling, comparative bioinformatic validation, DNA sequencing and MS validation. Using the pipeline, 5379 peptides specific for Bacillus species and 36 peptides specific for B. anthracis were identified and validated. The 36 peptides, representing 30 proteins were derived from over 15 different clusters of orthologous group categories, including proteins involved in transcription, energy production/conservation as well as multifunctional proteins. We demonstrated that the peptide biomarkers identified in this study could be detected in a complex background, in which 0.1 μg of protein extract from B. anthracis was spiked into 9.90 μg of B. cereus protein extracts. The integration of both stable non-redundant peptides with molecular methodology for marker discovery and validation, improves the robustness of identifying and characterising candidate biomarkers for the identification of bacteria such as B. anthracis.

Characterization of an environmental strain of Bacillus thuringiensis from a hot spring in Western Himalayas

Bacillus anthracis, the etiological agent of anthrax, is responsible for a serious and often fatal disease of mammalian livestock and humans and is an important biological warfare agent. Bacillus sp. AKG was isolated from a hot spring in western Himalayas and species-specific primers targeting gyrB gene identified the strain as B. anthracis within cereus-group. Cloning, sequencing, and phylogenetic analysis of the partial gyrB sequence from strain AKG indicated a close affiliation with B. anthracis and a few recently isolated strains of B. thuringiensis (e.g., strain Al Hakam and serovar konkukian). Phylogenetic analysis of two other housekeeping genes, clpC and gdpD yielded similar results. This observation is further substantiated by phylogenetic reconstruction using concatenated sequences (1680 bases) of the three genes (gyrB, clpC, and gdpD). Phenotypic features indicated a non-anthracis affiliation for the strain AKG. A novel strategy to distinguish among strains of B. anthracis, B. cereus, and B. thuringiensis based on whole proteome comparison was developed and tested for the identification of this environmental strain. Proteome comparison was used to establish the identity of this unknown environmental strain. Group of replicate 2DE gels for whole cell proteome were generated for each of the three species and strain AKG. Protein spots unique to each group and those showing match between the groups, in a pair-wise comparison, indicated strain AKG as a member of B. thuringiensis. This strategy can be used to assign strains of B. cereus group to their respective species.

Glycosyltransferase: a specific marker for the discrimination of Bacillus anthracis from the Bacillus cereus group

Bacillus anthracis, the aetiological agent of anthrax, has been taxonomically classified with the Bacillus cereus group, which comprises B. cereus, Bacillus thuringiensis, Bacillus mycoides, Bacillus pseudomycoides and Bacillus weihenstephanensis. Although the pathogenesis and ecological manifestations may be different, B. anthracis shares a high degree of DNA sequence similarity with its group member species. As a result, the discrimination of B. anthracis from its close relatives in the B. cereus group is still quite difficult. Suppression subtractive hybridization (SSH) was performed to search for genomic differences between a B. anthracis Korean isolate CR and the most closely related B. cereus type strain KCTC 3624(T). Two-hundred and five B. anthracis CR clones obtained by SSH underwent Southern hybridization, and comparative sequences were analysed using the blast program from the National Center for Biotechnology Information (NCBI). Subsequently, primer sets based on the glycosyltransferase group 1 family protein gene specific to B. anthracis were designed from the sequences of subtracted clones, and their specificities were evaluated using eight B. anthracis, 33 B. cereus, 10 B. thuringiensis, six B. mycoides, one B. pseudomycoides, one B. weihenstephanensis and 19 strains from 11 other representative Bacillus species. PCR primers specific for the glycosyltransferase group 1 family protein gene did not amplify the desired products from any of the Bacillus strains under examination, except B. anthracis alone. These findings may be useful in the future development of efficient diagnostic tools for the rapid identification of B. anthracis from other members of the B. cereus group.

Genetic distribution of 295 Bacillus cereus group members based on adk-screening in combination with MLST (Multilocus Sequence Typing) used for validating a primer targeting a chromosomal locus in B. anthracis

The genetic distribution of 295 Bacillus cereus group members has been investigated by using a modified Multilocus Sequence Typing method (MLST). By comparing the nucleic acid sequence of the adk gene fragment, isolates of B. cereus group members most related to B. anthracis may be easily identified. The genetic distribution, with focus on the B. anthracis close neighbours, was used to evaluate a new primer set for specific identification of B. anthracis. This primer set, BA5510-1/2, targeted the putative B. anthracis specific gene BA5510. Real-time PCR using BA5510-1/2 amplified the target fragment from all B. anthracis strains tested and only two (of 289) non-B. anthracis strains analysed. This is one of the most thoroughly validated chromosomal B. anthracis markers for real-time PCR identification, in which the screened collection contained several very closely related B. anthracis strains.

Proteomic profiling and identification of immunodominant spore antigens of Bacillus anthracis, Bacillus cereus, and Bacillus thuringiensis

Differentially expressed and immunogenic spore proteins of the Bacillus cereus group of bacteria, which includes Bacillus anthracis, Bacillus cereus, and Bacillus thuringiensis, were identified. Comparative proteomic profiling of their spore proteins distinguished the three species from each other as well as the virulent from the avirulent strains. A total of 458 proteins encoded by 232 open reading frames were identified by matrix-assisted laser desorption ionization-time-of-flight mass spectrometry analysis for all the species. A number of highly expressed proteins, including elongation factor Tu (EF-Tu), elongation factor G, 60-kDa chaperonin, enolase, pyruvate dehydrogenase complex, and others exist as charge variants on two-dimensional gels. These charge variants have similar masses but different isoelectric points. The majority of identified proteins have cellular roles associated with energy production, carbohydrate transport and metabolism, amino acid transport and metabolism, posttranslational modifications, and translation. Novel vaccine candidate proteins were identified using B. anthracis polyclonal antisera from humans postinfected with cutaneous anthrax. Fifteen immunoreactive proteins were identified in B. anthracis spores, whereas 7, 14, and 7 immunoreactive proteins were identified for B. cereus and in the virulent and avirulent strains of B. thuringiensis spores, respectively. Some of the immunodominant antigens include charge variants of EF-Tu, glyceraldehyde-3-phosphate dehydrogenase, dihydrolipoamide acetyltransferase, Delta-1-pyrroline-5-carboxylate dehydrogenase, and a dihydrolipoamide dehydrogenase. Alanine racemase and neutral protease were uniquely immunogenic to B. anthracis. Comparative analysis of the spore immunome will be of significance for further nucleic acid- and immuno-based detection systems as well as next-generation vaccine development.

Detailed genomic analysis of the Wbeta and gamma phages infecting Bacillus anthracis: implications for evolution of environmental fitness and antibiotic resistance

Phage-mediated lysis has been an essential laboratory tool for rapidly identifying Bacillus anthracis for more than 40 years, relying on the gamma phage derivative of a Bacillus cereus prophage called W. The complete genomic sequences of the temperate W phage, referred to as Wbeta, and its lytic variant gamma were determined and found to encode 53 open reading frames each, spanning 40,864 bp and 37,373 bp, respectively. Direct comparison of the genomes showed that gamma evolved through mutations at key loci controlling host recognition, lysogenic growth, and possibly host phenotypic modification. Included are a cluster of point mutations at the gp14 tail fiber locus of gamma, encoding a protein that, when fused to green fluorescent protein, binds specifically to B. anthracis. A large 2,003-bp deletion was also identified at the gamma lysogeny module, explaining its shift from a temperate to a lytic lifestyle. Finally, evidence of recombination was observed at a dicistronic Wbeta locus, encoding putative bacterial cell surface-modifying proteins, replaced in gamma with a locus, likely obtained from a B. anthracis prophage, encoding demonstrable fosfomycin resistance. Reverse transcriptase PCR analysis confirmed strong induction at the dicistronic Wbeta locus and at four other phage loci in B. anthracis and/or B. cereus lysogens. In all, this study represents the first genomic and functional description of two historically important phages and is part of a broader investigation into contributions of phage to the B. anthracis life cycle. Initial findings suggest that lysogeny of B. anthracis promotes ecological adaptation, rather than virulence, as with other gram-positive pathogens.

A cDNA macroarray approach to parasite-induced gene expression changes in a songbird host: genetic response of house finches to experimental infection by Mycoplasma gallisepticum

In 1994, the bacterial parasite Mycoplasma gallisepticum expanded its host range and swept through populations of a novel host--eastern US populations of the house finch (Carpodacus mexicanus). This epizootic caused a dramatic decline in finch population numbers, has been shown to have caused strong selection on house finch morphology, and presumably caused evolutionary change at the molecular level as finches evolved enhanced resistance. As a first step toward identifying finch genes that respond to infection by Mycoplasma and which may have experienced natural selection by this parasite, we used suppression subtractive hybridization (SSH) and cDNA macroarray approaches to identify differentially expressed genes regulated by the Mycoplasma parasite. Two subtractive cDNA libraries consisting of 16,512 clones were developed from spleen using an experimentally uninfected bird as the 'tester' and an infected bird as 'driver', and vice versa. Two hundred and twenty cDNA clones corresponding 34 genes with known vertebrate homologues and a large number of novel transcripts were found to be qualitatively up- or down-regulated genes by high-density filter hybridization. These gene expression changes were further confirmed by a high throughout reverse Northern blot approach and in specific cases by targeted Northern analysis. blast searches show that heat shock protein (HSP) 90, MHC II-associated invariant chain (CD74), T-cell immunoglobulin mucin 1 (TIM1), as well as numerous novel expressed genes not found in the databases were up- or down-regulated by the host in response to this parasite. Our results and macroarray resources provide a foundation for molecular co-evolutionary studies of the Mycoplasma parasite and its recently colonized avian host.

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.