Comparative transcriptional profiling of Bacillus cereus sensu lato strains during growth in CO2-bicarbonate and aerobic atmospheres

Effect of CO2 Concentrations on Entomopathogen Fitness and Insect-Pathogen Interactions

Numerous insect species and their associated microbial pathogens are exposed to elevated CO2 concentrations in both artificial and natural environments. However, the impacts of elevated CO2 on the fitness of these pathogens and the susceptibility of insects to pathogen infections are not well understood. The yellow mealworm, Tenebrio molitor, is commonly produced for food and feed purposes in mass-rearing systems, which increases risk of pathogen infections. Additionally, entomopathogens are used to control T. molitor, which is also a pest of stored grains. It is therefore important to understand how elevated CO2 may affect both the pathogen directly and impact on host-pathogen interactions. We demonstrate that elevated CO2 concentrations reduced the viability and persistence of the spores of the bacterial pathogen Bacillus thuringiensis. In contrast, conidia of the fungal pathogen Metarhizium brunneum germinated faster under elevated CO2. Pre-exposure of the two pathogens to elevated CO2 prior to host infection did not affect the survival probability of T. molitor larvae. However, larvae reared at elevated CO2 concentrations were less susceptible to both pathogens compared to larvae reared at ambient CO2 concentrations. Our findings indicate that whilst elevated CO2 concentrations may be beneficial in reducing host susceptibility in mass-rearing systems, they may potentially reduce the efficacy of the tested entomopathogens when used as biological control agents of T. molitor larvae. We conclude that CO2 concentrations should be carefully selected and monitored as an additional environmental factor in laboratory experiments investigating insect-pathogen interactions.

Effects of carbon dioxide and oxygen on the growth rate of various food spoilage bacteria

The growth of six bacterial species (Carnobacterium maltaromaticum, Bacillus weihenstephanensis, Bacillus cereus, Paenibacillus spp., Leuconostoc mesenteroides and Pseudomonas fragi) was studied in various gas compositions. Growth curves were obtained at various oxygen concentrations (between 0.1 and 21%), or various carbon dioxide concentrations (between 0 and 100%). Decreasing the O₂ concentration from 21% to about 3-5% has no effect on the bacterial growth rates, which are only affected by low oxygen levels. For each strain studied, the growth rate decreased linearly with carbon dioxide concentration, except for L. mesenteroides which remained insensible to this gas. Conversely, the most sensitive strain was totally inhibited by 50% of carbon dioxide in the gas phase at 8 °C. Predictive models were fitted, and the parameters characterizing the inhibitory effect of these two gases were estimated. This study provides new tools to help the food industry design suitable packaging for MAP storage.

From cereus to anthrax and back again: Assessment of the temperature-dependent phenotypic switching in the "cross-over" strain Bacillus cereus G9241

Bacillus cereus G9241 was isolated from a Louisiana welder suffering from an anthrax-like infection. The organism carries two transcriptional regulators that have previously been proposed to be incompatible with each other in Bacillus anthracis: the pleiotropic transcriptional regulator PlcR found in most members of the Bacillus cereus group but truncated in all B. anthracis isolates, and the anthrax toxin regulator AtxA found in all B. anthracis strains and a few B. cereus sensu stricto strains. Here we report cytotoxic and hemolytic activity of cell free B. cereus G9241 culture supernatants cultured at 25°C to various eukaryotic cells. However, this is not observed at the mammalian infection relevant temperature 37°C, behaving much like the supernatants generated by B. anthracis. Using a combination of genetic and proteomic approaches to understand this unique phenotype, we identified several PlcR-regulated toxins to be secreted highly at 25°C compared to 37°C. Furthermore, results suggest that differential expression of the protease involved in processing the PlcR quorum sensing activator molecule PapR appears to be the limiting step for the production of PlcR-regulated toxins at 37°C, giving rise to the temperature-dependent hemolytic and cytotoxic activity of the culture supernatants. This study provides an insight on how B. cereus G9241 is able to "switch" between B. cereus and B. anthracis-like phenotypes in a temperature-dependent manner, potentially accommodating the activities of both PlcR and AtxA.

Direct Regulons of AtxA, the Master Virulence Regulator of Bacillus anthracis

AtxA, the master virulence regulator of Bacillus anthracis, regulates the expression of three toxins and genes for capsule formation that are required for the pathogenicity of B. anthracis. Recent transcriptome analyses showed that AtxA affects a large number of genes on the chromosome and plasmids, suggesting a role as a global regulator. However, information on genes directly regulated by AtxA is scarce. In this work, we conducted genome-wide analyses and cataloged the binding sites of AtxA in vivo and transcription start sites on the B. anthracis genome. By integrating these results, we detected eight genes as direct regulons of AtxA. These consisted of five protein-coding genes, including two of the three toxin genes, and three genes encoding the small RNAs XrrA and XrrB and a newly discovered 95-nucleotide small RNA, XrrC. Transcriptomes from single-knockout mutants of these small RNAs revealed changes in the transcription levels of genes related to the aerobic electron transport chain, heme biosynthesis, and amino acid metabolism, suggesting their function for the control of cell physiology. These results reveal the first layer of the gene regulatory network for the pathogenicity of B. anthracis and provide a data set for the further study of the genomics and genetics of B. anthracis.

IMPORTANCEBacillus anthracis is the Gram-positive bacterial species that causes anthrax. Anthrax is still prevalent in countries mainly in Asia and Africa, where it causes economic damage and remains a public health issue. The mechanism of pathogenicity is mainly explained by the three toxin proteins expressed from the pXO1 plasmid and by proteins involved in capsule formation expressed from the pXO2 plasmid. AtxA is a protein expressed from the pXO1 plasmid that is known to upregulate genes involved in toxin production and capsule formation and is thus considered the master virulence regulator of B. anthracis. Therefore, understanding the detailed mechanism of gene regulation is important for the control of anthrax. The significance of this work lies in the identification of genes that are directly regulated by AtxA via genome-wide analyses. The results reveal the first layer of the gene regulatory network for the pathogenicity of B. anthracis and provide useful resources for a further understanding of B. anthracis.

You Can't B. cereus - A Review of Bacillus cereus Strains That Cause Anthrax-Like Disease

Emerging strains of Bacillus cereus, traditionally considered a self-limiting foodborne pathogen, have been associated with anthrax-like disease in mammals, including humans. The strains have emerged by divergent evolution and, as exchange of genetic material in the Bacillus genus occurs naturally, it is possible that further isolates will be identified in the future. The strains vary in their genotypes and phenotypes, combining traits of both B. cereus and B. anthracis species. Cases of anthrax-like disease associated with these strains result in similar symptoms and mortality rates as those caused by B. anthracis. The strains are susceptible to frontline antibiotics used in the treatment of anthrax and existing vaccines provide protection in animal models. The emergence of these strains has reignited the debate surrounding classification of the B. cereus sensu lato group and serves as a reminder that the field of medical microbiology is constantly changing and remains an important and ongoing area of research.

The Bacillus cereus Group: Bacillus Species with Pathogenic Potential

The Bacillus cereus group includes several Bacillus species with closely related phylogeny. The most well-studied members of the group, B. anthracis, B. cereus, and B. thuringiensis, are known for their pathogenic potential. Here, we present the historical rationale for speciation and discuss shared and unique features of these bacteria. Aspects of cell morphology and physiology, and genome sequence similarity and gene synteny support close evolutionary relationships for these three species. For many strains, distinct differences in virulence factor synthesis provide facile means for species assignment. B. anthracis is the causative agent of anthrax. Some B. cereus strains are commonly recognized as food poisoning agents, but strains can also cause localized wound and eye infections as well as systemic disease. Certain B. thuringiensis strains are entomopathogens and have been commercialized for use as biopesticides, while some strains have been reported to cause infection in immunocompromised individuals. In this article we compare and contrast B. anthracis, B. cereus, and B. thuringiensis, including ecology, cell structure and development, virulence attributes, gene regulation and genetic exchange systems, and experimental models of disease.

Link between Heterotrophic Carbon Fixation and Virulence in the Porcine Lung Pathogen Actinobacillus pleuropneumoniae

Actinobacillus pleuropneumoniae is a capnophilic pathogen of the porcine respiratory tract lacking enzymes of the oxidative branch of the tricarboxylic acid (TCA) cycle. We previously claimed that A. pleuropneumoniae instead uses the reductive branch in order to generate energy and metabolites. Here, we show that bicarbonate and oxaloacetate supported anaerobic growth of A. pleuropneumoniae Isotope mass spectrometry revealed heterotrophic fixation of carbon from stable isotope-labeled bicarbonate by A. pleuropneumoniae, which was confirmed by nano-scale secondary ion mass spectrometry at a single-cell level. By gas chromatography-combustion-isotope ratio mass spectrometry we could further show that the labeled carbon atom is mainly incorporated into the amino acids aspartate and lysine, which are derived from the TCA metabolite oxaloacetate. We therefore suggest that carbon fixation occurs at the interface of glycolysis and the reductive branch of the TCA cycle. The heme precursor δ-aminolevulinic acid supported growth of A. pleuropneumoniae, similar to bicarbonate, implying that anaplerotic carbon fixation is needed for heme synthesis. However, deletion of potential carbon-fixing enzymes, including PEP-carboxylase (PEPC), PEP-carboxykinase (PEPCK), malic enzyme, and oxaloacetate decarboxylase, as well as various combinations thereof, did not affect carbon fixation. Interestingly, generation of a deletion mutant lacking all four enzymes was not possible, suggesting that carbon fixation in A. pleuropneumoniae is an essential metabolic pathway controlled by a redundant set of enzymes. A double deletion mutant lacking PEPC and PEPCK was not impaired in carbon fixation in vitro but showed reduction of virulence in a pig infection model.

High Levels of CO2 Induce Spoilage by Leuconostoc mesenteroides by Upregulating Dextran Synthesis Genes

During nonventilated storage of carrots, CO2 gradually accumulates to high levels and causes modifications in the carrot's microbiome toward dominance of Lactobacillales and Enterobacteriales The lactic acid bacterium Leuconostoc mesenteroides secretes a slimy exudate over the surface of the carrots. The objective of this study was to characterize the slime components and the potential cause for its secretion under high CO2 levels. A proteomic analysis of the exudate revealed bacterial glucosyltransferases as the main proteins, specifically, dextransucrase. A chemical analysis of the exudate revealed high levels of dextran and several simple sugars. The exudate volume and dextran amount were significantly higher when L. mesenteroides was incubated under high CO2 levels than when incubated in an aerated environment. The treatment of carrot medium plates with commercial dextransucrase or exudate protein extract resulted in similar sugar profiles and dextran production. Transcriptome analysis demonstrated that dextran production is related to the upregulation of the L. mesenteroides dextransucrase-encoding genes dsrD and dsrT during the first 4 to 8 h of exposure to high CO2 levels compared to aerated conditions. A phylogenetic analysis of L. mesenteroides YL48 dsrD revealed a high similarity to other dsr genes harbored by different Leuconostoc species. The ecological benefit of dextran production under elevated CO2 requires further investigation. However, this study implies an overlooked role of CO2 in the physiology and fitness of L. mesenteroides in stored carrots, and perhaps in other food items, during storage under nonventilated conditions.IMPORTANCE The bacterium Leuconostoc mesenteroides is known to cause spoilage of different types of foods by secreting a slimy fluid that damages the quality and appearance of the produce. Here, we identified a potential mechanism by which high levels of CO2 affect the spoilage caused by this bacterium by upregulating dextran synthesis genes. These results have broader implications for the study of the physiology, degradation ability, and potential biotechnological applications of Leuconostoc.

Certhrax Is an Antivirulence Factor for the Anthrax-Like Organism Bacillus cereus Strain G9241

Bacillus cereus G9241 caused a life-threatening anthrax-like lung infection in a previously healthy human. This strain harbors two large virulence plasmids, pBCXO1 and pBC210, that are absent from typical B. cereus isolates. The pBCXO1 plasmid is nearly identical to pXO1 from Bacillus anthracis and carries genes (pagA1, lef, and cya) for anthrax toxin components (protective antigen [called PA1 in G9241], lethal factor [LF], and edema factor [EF], respectively). The plasmid also has an intact hyaluronic acid capsule locus. The pBC210 plasmid has a tetrasaccharide capsule locus, a gene for a PA1 homolog called PA2 (pagA2), and a gene (cer) for Certhrax, an ADP-ribosyltransferase toxin that inactivates vinculin. LF, EF, and Certhrax require PA for entry into cells. In this study, we asked what role PA1, PA2, LF, and Certhrax play in the pathogenicity of G9241. To answer this, we generated isogenic deletion mutations in the targeted toxin gene components and then assessed the strains for virulence in highly G9241-susceptible (A/J) and moderately G9241-sensitive (C57BL/6) mice. We found that full virulence of G9241 required PA1 and LF, while PA2 contributed minimally to pathogenesis of G9241 but could not functionally replace PA1 as a toxin-binding subunit in vivo Surprisingly, we discovered that Certhrax attenuated the virulence of G9241; i.e., a Δcer Δlef mutant strain was more virulent than a Δlef mutant strain following subcutaneous inoculation of A/J mice. Moreover, the enzymatic activity of Certhrax contributed to this phenotype. We concluded that Certhrax acts as an antivirulence factor in the anthrax-like organism B. cereus G9241.

RNA Sequencing and Analysis in Microorganisms for Metabolic Network Reconstruction

There is a strict interplay between metabolic networks and transcriptional regulation in bacteria; indeed, the transcriptome regulation, affecting the expression of large gene sets, can be used to predict the likely "on" or "off" state of metabolic genes as a function of environmental factors. Up to date, many bacterial transcriptomes have been studied by RNAseq, hundreds of experiments have been performed, and Giga bases of sequences have been produced. All this transcriptional information could potentially be integrated into metabolic networks in order to obtain a more comprehensive view of their regulation and to increase their prediction power.To get high-quality transcriptomic data, to be integrated into metabolic networks, it is paramount to clearly know how to produce highly informative RNA sequencing libraries and how to manage RNA sequencing data.In this chapter, we will get across the main steps of an RNAseq experiment: from removal of ribosomal RNAs, to strand-specific library preparation, till data analysis and integration. We will try to share our experience and know-how, to give you a precise protocol to follow, and some useful recommendations or tips and tricks to adopt in order to go straightforward toward a successful RNAseq experiment.

Severe hypercapnia and outcome of mechanically ventilated patients with moderate or severe acute respiratory distress syndrome

PURPOSE

To analyze the relationship between hypercapnia developing within the first 48 h after the start of mechanical ventilation and outcome in patients with acute respiratory distress syndrome (ARDS).

PATIENTS AND METHODS

We performed a secondary analysis of three prospective non-interventional cohort studies focusing on ARDS patients from 927 intensive care units (ICUs) in 40 countries. These patients received mechanical ventilation for more than 12 h during 1-month periods in 1998, 2004, and 2010. We used multivariable logistic regression and a propensity score analysis to examine the association between hypercapnia and ICU mortality.

MAIN OUTCOMES

We included 1899 patients with ARDS in this study. The relationship between maximum PaCO₂ in the first 48 h and mortality suggests higher mortality at or above PaCO₂ of ≥50 mmHg. Patients with severe hypercapnia (PaCO₂ ≥50 mmHg) had higher complication rates, more organ failures, and worse outcomes. After adjusting for age, SAPS II score, respiratory rate, positive end-expiratory pressure, PaO₂/FiO₂ ratio, driving pressure, pressure/volume limitation strategy (PLS), corrected minute ventilation, and presence of acidosis, severe hypercapnia was associated with increased risk of ICU mortality [odds ratio (OR) 1.93, 95% confidence interval (CI) 1.32 to 2.81; p = 0.001]. In patients with severe hypercapnia matched for all other variables, ventilation with PLS was associated with higher ICU mortality (OR 1.58, CI 95% 1.04-2.41; p = 0.032).

CONCLUSIONS

Severe hypercapnia appears to be independently associated with higher ICU mortality in patients with ARDS.

TRIAL REGISTRATION

Clinicaltrials.gov identifier, NCT01093482.

Changes in Bacillus anthracis CodY regulation under host-specific environmental factor deprived conditions

Background

Host-specific environmental factors induce changes in Bacillus anthracis gene transcription during infection. A global transcription regulator, CodY, plays a pivotal role in regulating central metabolism, biosynthesis, and virulence in B. anthracis. In this study, we utilized RNA-sequencing to assess changes in the transcriptional patterns of CodY-regulated B. anthracis genes in response to three conditions of environmental starvation: iron, CO₂, or glucose deprivation. In addition, we performed chromatin immunoprecipitation on newly identified CodY-mediated genes.

Results

Environmental deprivation induced transcriptional changes in CodY-regulated genes in both wild-type and codY null strains, and both CodY-specific and environment-specific patterns were observed. In the iron-depleted condition, overexpression of iron homeostasis genes was observed independent of codY deletion; however, transcription of siderophore and amino acid biosynthesis genes was CodY dependent. Although CodY has a significant regulatory role in central metabolism and the carbon overflow pathway, metabolism-associated genes exhibited CodY-independent expression patterns under glucose starvation. Genes that were differentially expressed in response to CO₂ availability showed CodY-dependent regulation, though their maximal expression did require a supply of CO₂/bicarbonate.

Conclusions

We speculate that CodY regulates the expression of environmental-responsive genes in a hierarchical manner and is likely associated with other transcription regulators that are specific for a particular environmental change.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3004-8) contains supplementary material, which is available to authorized users.

Capsules, toxins and AtxA as virulence factors of emerging Bacillus cereus biovar anthracis

Emerging B. cereus strains that cause anthrax-like disease have been isolated in Cameroon (CA strain) and Côte d’Ivoire (CI strain). These strains are unusual, because their genomic characterisation shows that they belong to the B. cereus species, although they harbour two plasmids, pBCXO1 and pBCXO2, that are highly similar to the pXO1 and pXO2 plasmids of B. anthracis that encode the toxins and the polyglutamate capsule respectively. The virulence factors implicated in the pathogenicity of these B. cereus bv anthracis strains remain to be characterised. We tested their virulence by cutaneous and intranasal delivery in mice and guinea pigs; they were as virulent as wild-type B. anthracis. Unlike as described for pXO2-cured B. anthracis, the CA strain cured of the pBCXO2 plasmid was still highly virulent, showing the existence of other virulence factors. Indeed, these strains concomitantly expressed a hyaluronic acid (HA) capsule and the B. anthracis polyglutamate (PDGA) capsule. The HA capsule was encoded by the hasACB operon on pBCXO1, and its expression was regulated by the global transcription regulator AtxA, which controls anthrax toxins and PDGA capsule in B. anthracis. Thus, the HA and PDGA capsules and toxins were co-regulated by AtxA. We explored the respective effect of the virulence factors on colonisation and dissemination of CA within its host by constructing bioluminescent mutants. Expression of the HA capsule by itself led to local multiplication and, during intranasal infection, to local dissemination to the adjacent brain tissue. Co-expression of either toxins or PDGA capsule with HA capsule enabled systemic dissemination, thus providing a clear evolutionary advantage. Protection against infection by B. cereus bv anthracis required the same vaccination formulation as that used against B. anthracis. Thus, these strains, at the frontier between B. anthracis and B. cereus, provide insight into how the monomorphic B. anthracis may have emerged.

Anthrax is caused by the bacterium Bacillus anthracis that affects all mammals worldwide. It emerged more than 10,000 years ago from a Bacillus cereus precursor. In the past decade, B. cereus bacteria were isolated in the USA from anthrax-like pneumonia cases. They harbour one virulence plasmid very similar to the toxin–encoding plasmid of B. anthracis. Recently, an anthrax-like disease in great apes in Africa was caused by emerging B. cereus strains, named B. cereus biovar anthracis. These strains are atypical as they possess both plasmids coding for toxin and capsule similar to those so far found only in B. anthracis. These unusual pathogenic B. cereus are currently neglected. We explored the virulence of these pathogens and their colonisation and dissemination capacity within the murine host. We found that these toxinogenic strains harbour two capsules, the classical B. anthracis capsule and an additional polysaccharidic capsule. This latter capsule confers virulence alone or in combination with toxins. Both capsules are concomitantly expressed, under the control of a common global regulator and host signals. Our results show that acquisition of new genetic information by these B. cereus clearly gives them a selective advantage, favouring their dissemination within infected hosts and the environment.

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.

Identification of pBC218/pBC210 Genes of Bacillus cereus G9241 in Five Florida Soils Using qPCR

The distribution of the virulent plasmid pBC210 of B. cereus that carries several B. anthracis genes and has been implicated in lethal anthrax-like pulmonary disease is unknown. We screened our collection of 103 B. cereus isolates and 256 soil samples using a quantitative PCR (qPCR) assay that targeted three open reading frames putatively unique to pBC210. When tested with DNA from 2 B. cereus strains carrying pBC210, and 64 Gram-positive and 55 Gram-negative bacterial species, the assay had 100% sensitivity and specificity. None of the DNA from the B. cereus isolates yielded positive amplicons but DNA extracted from five soils collected in Florida gave positive results for all three target sequences of pBC210. While screening confirms that pBC210 is uncommon in B. cereus, this study is the first to report that pBC210 is present in Florida soils. This study improves our knowledge of the distribution of pBC210 in soils and, of public health importance, the potential threat of B. cereus isolates carrying the toxin-carrying plasmid. We demonstrated that sequences of pBC210 can be found in a larger geographical area than previously thought and that finding more B. cereus carrying the virulent plasmid is a possibility in the future.

Transcriptome analysis identifies Bacillus anthracis genes that respond to CO2 through an AtxA-dependent mechanism

Background

Upon infection of a mammalian host, Bacillus anthracis responds to host cues, and particularly to elevated temperature (37°C) and bicarbonate/CO₂ concentrations, with increased expression of virulence factors that include the anthrax toxins and extracellular capsular layer. This response requires the presence of the pXO1 virulence plasmid-encoded pleiotropic regulator AtxA. To better understand the genetic basis of this response, we utilized a controlled in vitro system and Next Generation sequencing to determine and compare RNA expression profiles of the parental strain and an isogenic AtxA-deficient strain in a 2 × 2 factorial design with growth environments containing or lacking carbon dioxide.

Results

We found 15 pXO1-encoded genes and 3 chromosomal genes that were strongly regulated by the separate or synergistic actions of AtxA and carbon dioxide. The majority of the regulated genes responded to both AtxA and carbon dioxide rather than to just one of these factors. Interestingly, we identified two previously unrecognized small RNAs that are highly expressed under physiological carbon dioxide concentrations in an AtxA-dependent manner. Expression levels of the two small RNAs were found to be higher than that of any other gene differentially expressed in response to these conditions. Secondary structure and small RNA-mRNA binding predictions for the two small RNAs suggest that they may perform important functions in regulating B. anthracis virulence.

Conclusions

A majority of genes on the virulence plasmid pXO1 that are regulated by the presence of either CO₂ or AtxA separately are also regulated synergistically in the presence of both. These results also elucidate novel pXO1-encoded small RNAs that are associated with virulence conditions.

Exoproteome analysis of a novel strain of Bacillus cereus implicated in disease resembling cutaneous anthrax

Bacillus cereus belongs to B. cereus sensu lato group, shared by six other related species including Bacillus anthracis. B. anthracis is the causative agent for serious illness affecting a wide range of animals as well as humans and is a category A Biological and Toxin Warfare (BTW) agent. Recent studies indicate that a Bacillus species other than B. anthracis can cause anthrax-like disease and role of anthrax virulence plasmids (pXO1 and pXO2) on the pathogenicity of B. cereus has been documented. B. cereus strain TF5 was isolated from the tissue fluid of cutaneous anthrax-like skin lesions of a human patient from an anthrax endemic area in India. The strain harboured a PA gene, however, presence of pXO1 or pXO2-like plasmids could not be ascertained using reported primers. Abundant exoproteome of the strain in the early stationary phase was elucidated using a 2-DE MS approach and compared with that from a reference B. cereus strain. Analysis of proteins showing qualitative and quantitative differences between the two strains indicated an altered regulatory mechanism and putative role of S-layer protein and sphingomyelinase in the pathogenesis of strain TF5. Phylogenetic analysis of the S-layer protein indicated close affiliation of the strain with anthracis-like B. cereus strains such as B. cereus var. anthracis strain CI; whereas sphingomyelinase exhibited specific relationship with all the strains of B. anthracis apart from that with anthracis-like B. cereus strains.

Delineating the effect of host environmental signals on a fully virulent strain of Bacillus anthracis using an integrated transcriptomics and proteomics approach

Pathogenic bacteria sense the host environment and regulate expression of virulence-related genes. Environmental signals like temperature, bicarbonate/CO2 and glucose induce toxin production in Bacillus anthracis, but the mechanisms by which these signals contribute to virulence and overall physiological adaptation remains elusive. An integrated, systems level investigation using transcriptomics and iTRAQ-based proteomics was done to assess the effect of temperature, bicarbonate/CO2 and glucose on B. anthracis. Significant changes observed in amino acid, carbohydrate, energy and nucleotide metabolism indicates events of metabolic readjustments by environmental factors. Directed induction of genes involved in polyamine biosynthesis and iron metabolism revealed the redirection of cellular metabolite pool towards iron uptake. Protein levels of glycolytic enzymes, ptsH and Ldh along with transcripts involved in immune evasion (mprF, bNOS, Phospholipases and asnA), cell surface remodeling (rfbABCD, antABCD, and cls) and utilization of lactate (lutABC) and inositol showed constant repression under environmental perturbations. Discrepancies observed in mRNA/protein level of genes involved in glycolysis, protein synthesis, stress response and nucleotide metabolism hinted at the existence of additional regulatory layers and illustrated the utility of an integrated approach. The above findings might assist in the identification of novel adaptive strategies of B. anthracis during host associated survival and pathogenesis.

BIOLOGICAL SIGNIFICANCE

In this study, the changes observed at both transcript and protein level were quantified and integrated to understand the effect of host environmental factors (host temperature, bicarbonate and glucose) in shaping the physiology and adaptive strategies of a fully virulent strain of B. anthracis for efficient survival and virulence in its host. Perturbations affecting toxin production were found to concordantly affect vital metabolic pathways and several known as well as novel virulence factors. These changes act as a valuable asset for generating testable hypotheses that can be further verified by detailed molecular and mutant studies to identify novel adaptive strategies of B. anthracis during infection. Adaptation of an integrated transcriptomics and proteomics approach also led to the identification of discrepancies between mRNA/protein levels among genes across major functional categories. Few of these discrepancies have been previously reported in literature for model organisms. However their existence in B. anthracis and that too as a result of growth perturbations have not been reported till date. These findings demonstrate a substantial role of regulatory processes post mRNA synthesis via post transcriptional, translational or protein degradation mechanisms. This article is part of a Special Issue entitled: Proteomics of non-model organisms.

Carbon dioxide-sensing in organisms and its implications for human disease

The capacity of organisms to sense changes in the levels of internal and external gases and to respond accordingly is central to a range of physiologic and pathophysiologic processes. Carbon dioxide, a primary product of oxidative metabolism is one such gas that can be sensed by both prokaryotic and eukaryotic cells and in response to altered levels, elicit the activation of multiple adaptive pathways. The outcomes of activating CO2-sensitive pathways in various species include increased virulence of fungal and bacterial pathogens, prey-seeking behavior in insects as well as taste perception, lung function, and the control of immunity in mammals. In this review, we discuss what is known about the mechanisms underpinning CO2 sensing across a range of species and consider the implications of this for physiology, disease progression, and the possibility of developing new therapeutics for inflammatory and infectious disease.

An efficient rRNA removal method for RNA sequencing in GC-rich bacteria

Background

Next generation sequencing (NGS) technologies have revolutionized gene expression studies and functional genomics analysis. However, further improvement of RNA sequencing protocols is still desirable, in order to reduce NGS costs and to increase its accuracy. In bacteria, a major problem in RNA sequencing is the abundance of ribosomal RNA (rRNA), which accounts for 95-98% of total RNA and can therefore hinder sufficient coverage of mRNA, the main focus of transcriptomic studies. Thus, efficient removal of rRNA is necessary to achieve optimal coverage, good detection sensitivity and reliable results. An additional challenge is presented by microorganisms with GC-rich genomes, in which rRNA removal is less efficient.

Results

In this work, we tested two commercial kits for rRNA removal, either alone or in combination, on Burkholderia thailandensis. This bacterium, chosen as representative of the important Burkholderia genus, which includes both pathogenic and environmental bacteria, has a rather large (6.72 Mb) and GC-rich (67.7%) genome. Each enriched mRNA sample was sequenced through paired-end Illumina GAIIx run in duplicate, yielding between 10 and 40 million reads. We show that combined treatment with both kits allows an mRNA enrichment of more than 238-fold, enabling the sequencing of almost all (more than 90%) B. thailandensis transcripts from less than 10 million reads, without introducing any bias in mRNA relative abundance, thus preserving differential expression profile.

Conclusions

The mRNA enrichment protocol presented in this work leads to an increase in detection sensitivity up to 770% compared to total RNA; such increased sensitivity allows for a corresponding reduction in the number of sequencing reads necessary for the complete analysis of whole transcriptome expression profiling. Thus we can conclude that the MICROBExpress/Ovation combined rRNA removal method could be suitable for RNA sequencing of whole transcriptomes of microorganisms with high GC content and complex genomes enabling at the same time an important scaling down of sequencing costs.

Divergence of the SigB regulon and pathogenesis of the Bacillus cereus sensu lato group

Background

The Bacillus cereus sensu lato group currently includes seven species (B. cereus, B. anthracis, B. mycoides, B. pseudomycoides, B. thuringiensis, B. weihenstephanensis and B. cytotoxicus) that recent phylogenetic and phylogenomic analyses suggest are likely a single species, despite their varied phenotypes. Although horizontal gene transfer and insertion-deletion events are clearly important for promoting divergence among these genomes, recent studies have demonstrated that a major basis for phenotypic diversity in these organisms may be differential regulation of the highly similar gene content shared by these organisms. To explore this hypothesis, we used an in silico approach to evaluate the relationship of pathogenic potential and the divergence of the SigB-dependent general stress response within the B. cereus sensu lato group, since SigB has been demonstrated to support pathogenesis in Bacillus, Listeria and Staphylococcus species.

Results

During the divergence of these organisms from a common “SigB-less” ancestor, the placement of SigB promoters at varied locations in the B. cereus sensu lato genomes predict alternative structures for the SigB regulon in different organisms. Predicted promoter changes suggesting differential transcriptional control of a common gene pool predominate over evidence of indels or horizontal gene transfer for explaining SigB regulon divergence.

Conclusions

Four lineages of the SigB regulon have arisen that encompass different gene contents and suggest different strategies for supporting pathogenesis. This is consistent with the hypothesis that divergence within the B. cereus sensu lato group rests in part on alternative strategies for regulation of a common gene pool.

Strand-specific RNA-seq reveals ordered patterns of sense and antisense transcription in Bacillus anthracis

BACKGROUND

Although genome-wide transcriptional analysis has been used for many years to study bacterial gene expression, many aspects of the bacterial transcriptome remain undefined. One example is antisense transcription, which has been observed in a number of bacteria, though the function of antisense transcripts, and their distribution across the bacterial genome, is still unclear.

METHODOLOGY/PRINCIPAL FINDINGS

Single-stranded RNA-seq results revealed a widespread and non-random pattern of antisense transcription covering more than two thirds of the B. anthracis genome. Our analysis revealed a variety of antisense structural patterns, suggesting multiple mechanisms of antisense transcription. The data revealed several instances of sense and antisense expression changes in different growth conditions, suggesting that antisense transcription may play a role in the ways in which B. anthracis responds to its environment. Significantly, genome-wide antisense expression occurred at consistently higher levels on the lagging strand, while the leading strand showed very little antisense activity. Intrasample gene expression comparisons revealed a gene dosage effect in all growth conditions, where genes farthest from the origin showed the lowest overall range of expression for both sense and antisense directed transcription. Additionally, transcription from both strands was verified using a novel strand-specific assay. The variety of structural patterns we observed in antisense transcription suggests multiple mechanisms for this phenomenon, suggesting that some antisense transcription may play a role in regulating the expression of key genes, while some may be due to chromosome replication dynamics and transcriptional noise.

CONCLUSIONS/SIGNIFICANCE

Although the variety of structural patterns we observed in antisense transcription suggest multiple mechanisms for antisense expression, our data also clearly indicate that antisense transcription may play a genome-wide role in regulating the expression of key genes in Bacillus species. This study illustrates the surprising complexity of prokaryotic RNA abundance for both strands of a bacterial chromosome.

Methylation-dependent DNA restriction in Bacillus anthracis

Bacillus anthracis, the causative agent of anthrax, is poorly transformed with DNA that is methylated on adenine or cytosine. Here we characterize three genetic loci encoding type IV methylation-dependent restriction enzymes that target DNA containing C5-methylcytosine (m5C). Strains in which these genes were inactivated, either singly or collectively, showed increased transformation by methylated DNA. Additionally, a triple mutant with an ~30-kb genomic deletion could be transformed by DNA obtained from Dam(+)Dcm(+)E. coli, although at a low frequency of ~10(-3) transformants/10(6)cfu. This strain of B. anthracis can potentially serve as a preferred host for shuttle vectors that express recombinant proteins, including proteins to be used in vaccines. The gene(s) responsible for the restriction of m6A-containing DNA in B. anthracis remain unidentified, and we suggest that poor transformation by such DNA could in part be a consequence of the inefficient replication of hemimethylated DNA in B. anthracis.

Whole-genome phylogenies of the family Bacillaceae and expansion of the sigma factor gene family in the Bacillus cereus species-group

BACKGROUND

The Bacillus cereus sensu lato group consists of six species (B. anthracis, B. cereus, B. mycoides, B. pseudomycoides, B. thuringiensis, and B. weihenstephanensis). While classical microbial taxonomy proposed these organisms as distinct species, newer molecular phylogenies and comparative genome sequencing suggests that these organisms should be classified as a single species (thus, we will refer to these organisms collectively as the Bc species-group). How do we account for the underlying similarity of these phenotypically diverse microbes? It has been established for some time that the most rapidly evolving and evolutionarily flexible portions of the bacterial genome are regulatory sequences and transcriptional networks. Other studies have suggested that the sigma factor gene family of these organisms has diverged and expanded significantly relative to their ancestors; sigma factors are those portions of the bacterial transcriptional apparatus that control RNA polymerase recognition for promoter selection. Thus, examining sigma factor divergence in these organisms would concurrently examine both regulatory sequences and transcriptional networks important for divergence. We began this examination by comparison to the sigma factor gene set of B. subtilis.

RESULTS

Phylogenetic analysis of the Bc species-group utilizing 157 single-copy genes of the family Bacillaceae suggests that several taxonomic revisions of the genus Bacillus should be considered. Within the Bc species-group there is little indication that the currently recognized species form related sub-groupings, suggesting that they are members of the same species. The sigma factor gene family encoded by the Bc species-group appears to be the result of a dynamic gene-duplication and gene-loss process that in previous analyses underestimated the true heterogeneity of the sigma factor content in the Bc species-group.

CONCLUSIONS

Expansion of the sigma factor gene family appears to have preferentially occurred within the extracytoplasmic function (ECF) sigma factor genes, while the primary alternative (PA) sigma factor genes are, in general, highly conserved with those found in B. subtilis. Divergence of the sigma-controlled transcriptional regulons among various members of the Bc species-group likely has a major role in explaining the diversity of phenotypic characteristics seen in members of the Bc species-group.

Bacillus cereus G9241 makes anthrax toxin and capsule like highly virulent B. anthracis Ames but behaves like attenuated toxigenic nonencapsulated B. anthracis Sterne in rabbits and mice

Bacillus cereus G9241 was isolated from a welder with a pulmonary anthrax-like illness. The organism contains two megaplasmids, pBCXO1 and pBC218. These plasmids are analogous to the Bacillus anthracis Ames plasmids pXO1 and pXO2 that encode anthrax toxins and capsule, respectively. Here we evaluated the virulence of B. cereus G9241 as well as the contributions of pBCXO1 and pBC218 to virulence. B. cereus G9241 was avirulent in New Zealand rabbits after subcutaneous inoculation and attenuated 100-fold compared to the published 50% lethal dose (LD(50)) values for B. anthracis Ames after aerosol inoculation. A/J and C57BL/6J mice were comparably susceptible to B. cereus G9241 by both subcutaneous and intranasal routes of infection. However, the LD(50)s for B. cereus G9241 in both mouse strains were markedly higher than those reported for B. anthracis Ames and more like those of the toxigenic but nonencapsulated B. anthracis Sterne. Furthermore, B. cereus G9241 spores could germinate and disseminate after intranasal inoculation into A/J mice, as indicated by the presence of vegetative cells in the spleen and blood of animals 48 h after infection. Lastly, B. cereus G9241 derivatives cured of one or both megaplasmids were highly attenuated in A/J mice. We conclude that the presence of the toxin- and capsule-encoding plasmids pBCXO1 and pBC218 in B. cereus G9241 alone is insufficient to render the strain as virulent as B. anthracis Ames. However, like B. anthracis, full virulence of B. cereus G9241 for mice requires the presence of both plasmids.

AtxA, a Bacillus anthracis global virulence regulator

Fifteen years ago, AtxA was isolated as a toxin gene activator and five years later it was shown to be a Bacillus anthracis master regulator. AtxA controls the expression of more than a hundred genes belonging to all genetic elements, the chromosome and both virulence plasmids, including those encoding the major virulence factors. AtxA can activate or repress gene expression. The mechanism by which AtxA exerts its control is unknown; it is indirect on some genes but may be direct on others. The expression of many AtxA-controlled genes is induced by the presence of bicarbonate/CO(2). AtxA links the metabolic state and virulence gene expression.

Activation of the latent PlcR regulon in Bacillus anthracis

Many genes in Bacillus cereus and Bacillus thuringiensis are under the control of the transcriptional regulator PlcR and its regulatory peptide, PapR. In Bacillus anthracis, the causative agent of anthrax, PlcR is inactivated by truncation, and consequently genes having PlcR binding sites are expressed at very low levels when compared with B. cereus. We found that activation of the PlcR regulon in B. anthracis by expression of a PlcR–PapR fusion protein does not alter sporulation in strains containing the virulence plasmid pXO1 and thereby the global regulator AtxA. Using comparative 2D gel electrophoresis, we showed that activation of the PlcR regulon in B. anthracis leads to upregulation of many proteins found in the secretome of B. cereus, including phospholipases and proteases, such as the putative protease BA1995. Transcriptional analysis demonstrated expression of BA1995 to be dependent on PlcR–PapR, even though the putative PlcR recognition site of the BA1995 gene does not exactly match the PlcR consensus sequence, explaining why this protein had escaped recognition as belonging to the PlcR regulon. Additionally, while transcription of major PlcR-dependent haemolysins, sphingomyelinase and anthrolysin O is enhanced in response to PlcR activation in B. anthracis, only anthrolysin O contributes significantly to lysis of human erythrocytes. In contrast, the toxicity of bacterial culture supernatants from a PlcR-positive strain towards murine macrophages occurred independently of anthrolysin O expression in vitro and in vivo.

Proteomic studies of Bacillus anthracis

Bacillus anthracis is a Gram-positive, spore-forming bacterium representing the etiological cause of anthrax, a rare lethal disease of animals and humans. Development of anthrax countermeasures has gained increasing attention owing to the potential use of B. anthracis spores as a bioterror weapon. The various forms of infection by B. anthracis are characterized both by toxemia and septicemia, both of which are the result of spore entry into the host followed by their germination into rapidly multiplying, toxin-producing bacilli. Following the publication of the bacterial genome, proteomic studies were carried out to determine the protein composition of the spore and identify exposed vegetative (membrane-located or secreted) proteins. These studies included comparison of strains differing in their virulence, cultured under different conditions and, in some cases, were complemented by serological inspection, which addressed expression during infection of proteomically identified proteins and their immunogenicity. The proteomic approach emerged as a valuable strategy for the generation of a pool of potential B. anthracis protein targets for further evaluation in detection, diagnostics, therapy and prophylaxis, and contributed to the elucidation of some aspects of the pathogenesis of the disease.

Transcriptional profiling of Bacillus anthracis Sterne (34F2) during iron starvation

Lack of available iron is one of many environmental challenges that a bacterium encounters during infection and adaptation to iron starvation is important for the pathogen to efficiently replicate within the host. Here we define the transcriptional response of B. anthracis Sterne (34F₂) to iron depleted conditions. Genome-wide transcript analysis showed that B. anthracis undergoes considerable changes in gene expression during growth in iron-depleted media, including the regulation of known and candidate virulence factors. Two genes encoding putative internalin proteins were chosen for further study. Deletion of either gene (GBAA0552 or GBAA1340) resulted in attenuation in a murine model of infection. This attenuation was amplified in a double mutant strain. These data define the transcriptional changes induced during growth in low iron conditions and illustrate the potential of this dataset in the identification of putative virulence determinants for future study.

Structure and complexity of a bacterial transcriptome

Although gene expression has been studied in bacteria for decades, many aspects of the bacterial transcriptome remain poorly understood. Transcript structure, operon linkages, and information on absolute abundance all provide valuable insights into gene function and regulation, but none has ever been determined on a genome-wide scale for any bacterium. Indeed, these aspects of the prokaryotic transcriptome have been explored on a large scale in only a few instances, and consequently little is known about the absolute composition of the mRNA population within a bacterial cell. Here we report the use of a high-throughput sequencing-based approach in assembling the first comprehensive, single-nucleotide resolution view of a bacterial transcriptome. We sampled the Bacillus anthracis transcriptome under a variety of growth conditions and showed that the data provide an accurate and high-resolution map of transcript start sites and operon structure throughout the genome. Further, the sequence data identified previously nonannotated regions with significant transcriptional activity and enhanced the accuracy of existing genome annotations. Finally, our data provide estimates of absolute transcript abundance and suggest that there is significant transcriptional heterogeneity within a clonal, synchronized bacterial population. Overall, our results offer an unprecedented view of gene expression and regulation in a bacterial cell.