Evolution of pathogenicity in the Bacillus cereus group

Evaluating the Safety of Bacillus cereus GW-01 Obtained from Sheep Rumen Chyme

Bacillus cereus is responsible for 1.4-12% food poisoning outbreaks worldwide. The safety concerns associated with the applications of B. cereus in health and medicine have been controversial due to its dual role as a pathogen for foodborne diseases and a probiotic in humans and animals. In this study, the pathogenicity of B. cereus GW-01 was assessed by comparative genomic, and transcriptome analysis. Phylogenetic analysis based on a single-copy gene showed clustering of the strain GW-01, and 54 B. cereus strains from the NCBI were classified into six major groups (I-VI), which were then associated with the source region and sequence types (STs). Transcriptome results indicated that the expression of most genes related with toxins secretion in GW-01 was downregulated compared to that in the lag phase. Overall, these findings suggest that GW-01 is not directly associated with pathogenic Bacillus cereus and highlight an insightful strategy for assessing the safety of novel B. cereus strains.

New Approach Methods to Assess the Enteropathogenic Potential of Strains of the Bacillus cereus Group, including Bacillus thuringiensis

Bacillus cereus (Bc) is a wide group of Gram-positive and spore-forming bacteria, known to be the etiological agents of various human infections, primarily food poisoning. The Bc group includes enteropathogenic strains able to germinate in the digestive tract and to produce enterotoxins such as Nhe, Hbl, and CytK. One species of the group, Bacillus thuringiensis (Bt), has the unique feature of producing insecticidal crystals during sporulation, making it an important alternative to chemical pesticides to protect crops from insect pest larvae. Nevertheless, several studies have suggested a link between the ingestion of pesticide strains and human cases of food poisoning, calling their safety into question. Consequently, reliable tools for virulence assessment are worth developing to aid decision making in pesticide regulation. Here, we propose complementary approaches based on two biological models, the human intestinal Caco-2 cell line and the insect Drosophila melanogaster, to assess and rank the enteric virulence potency of Bt strains in comparison with other Bc group members. Using a dataset of 48 Bacillus spp. strains, we showed that some Bc group strains, including Bt, were able to induce cytotoxicity in Caco-2 cells with concomitant release of IL-8 cytokine, a landmark of pro-inflammatory response. In the D. melanogaster model, we were able to sort a panel of 39 strains into four different classes of virulence, ranging from no virulence to strong virulence. Importantly, for the most virulent strains, mortality was associated with a loss of intestinal barrier integrity. Interestingly, although strains can share a common toxinotype, they display different degrees of virulence, suggesting the existence of specific mechanisms of virulence expression in vivo in the intestine.

JMM Profile: Bacillus anthracis

Principal routes of Bacillus anthracis infection and stages of anthrax pathogenesis, consistent with current understandings. Depending on the route of infection, germination of spores may happen in extracellular tissue fluid, or following phagocytosis (a). Successful infection of host cells leads to toxin-associated cell death and release of vegetative cells and toxin (b). Toxin binds and enters other host cells (c), including those of the immune system, disrupting function. In some cases this leads to systemic disease, which typically is fatal.

High Genetic Diversity and Virulence Potential in Bacillus cereus sensu lato Isolated from Milk and Cheeses in Apulia Region, Southern Italy

The Bacillus cereus group includes species that act as food-borne pathogens causing diarrheal and emetic symptoms. They are widely distributed and can be found in various foods. In this study, out of 550 samples of milk and cheeses, 139 (25.3%) were found to be contaminated by B. cereus sensu lato (s.l.). One isolate per positive sample was characterized by Multilocus Sequence Typing (MLST) and for the presence of ten virulence genes. Based on MLST, all isolates were classified into 73 different sequence types (STs), of which 12 isolates were assigned to new STs. Virulence genes detection revealed that 90% and 61% of the isolates harboured the nheABC and the hblCDA gene cluster, respectively. Ninety-four percent of the isolates harboured the enterotoxin genes entS and entFM; 8% of the isolates possessed the ces gene. Thirty-eight different genetic profiles were identified, suggesting a high genetic diversity. Our study clearly shows the widespread diffusion of potentially toxigenic isolates of B. cereus s.l. in milk and cheeses in the Apulia region highlighting the need to adopt GMP and HACCP procedures along every step of the milk and cheese production chain in order to reduce the public health risk linked to the consumption of foods contaminated by B. cereus s.l.

Phylogenetic Analysis of Bacterial Pathogen Genomes

The development of high-throughput sequencing technology has led to a significant reduction in the time and cost of sequencing whole genomes of bacterial pathogens. Studies can sequence and compare hundreds or even thousands of genomes within a given bacterial population. A phylogenetic tree is the most frequently used method of depicting the relationships between these bacterial pathogen genomes. However, the presence of homologous recombination in most bacterial pathogen species can invalidate the application of standard phylogenetic tools. Here we describe a method to produce phylogenetic analyses that accounts for the disruptive effect of recombination. This allows users to investigate the recombination events that have occurred, as well as to produce more meaningful phylogenetic analyses which recover the clonal genealogy representing the clonal relationships between genomes.

A Comparative Analysis of the Core Proteomes within and among the Bacillus subtilis and Bacillus cereus Evolutionary Groups Reveals the Patterns of Lineage- and Species-Specific Adaptations

By integrating phylogenomic and comparative analyses of 1104 high-quality genome sequences, we identify the core proteins and the lineage-specific fingerprint proteins of the various evolutionary clusters (clades/groups/species) of the Bacillus genus. As fingerprints, we denote those core proteins of a certain lineage that are present only in that particular lineage and absent in any other Bacillus lineage. Thus, these lineage-specific fingerprints are expected to be involved in particular adaptations of that lineage. Intriguingly, with a few notable exceptions, the majority of the Bacillus species demonstrate a rather low number of species-specific fingerprints, with the majority of them being of unknown function. Therefore, species-specific adaptations are mostly attributed to highly unstable (in evolutionary terms) accessory proteomes and possibly to changes at the gene regulation level. A series of comparative analyses consistently demonstrated that the progenitor of the Cereus Clade underwent an extensive genomic expansion of chromosomal protein-coding genes. In addition, the majority (76–82%) of the B. subtilis proteins that are essential or play a significant role in sporulation have close homologs in most species of both the Subtilis and the Cereus Clades. Finally, the identification of lineage-specific fingerprints by this study may allow for the future development of highly specific vaccines, therapeutic molecules, or rapid and low-cost molecular tests for species identification.

Phylogenetic and protein prediction analysis reveals the taxonomically diverse distribution of virulence factors in Bacillus cereus strains

Bacillus cereus is a food contaminant with widely varying enterotoxic potential due to its virulence proteins. In this article, phylogenetic analysis of the amino acid sequences from the whole-genomes of 41 strains, evolutionary distance calculation of the amino acid sequences of the virulence genes, and functional and structural predictions of the virulence proteins were performed to reveal the taxonomically diverse distribution of virulence factors. The genome evolution of the strains showed a clustering trend based on the protein-coding virulence genes. The strains of B. cereus have evolved into non-toxic risk and toxic risk clusters with medium-high- and medium-low-risk subclusters. The evolutionary transfer distances of incomplete virulence genes relative to housekeeping genes were greater than those of complete virulence genes, and the distance values of HblACD were higher than those of nheABC and CytK among the complete virulence genes. Cytoplasmic localization was impossible for all the virulence proteins, and NheB, NheC, Hbl-B, and Hbl-L₁ were predicted to be extracellular. Nhe and Hbl proteins except CytK had similar spatial structures. The predicted structures of Nhe and Hbl mainly showed ‘head’ and ‘tail’ domains. The ‘head’ of NheA and Hbl-B, including two α-helices separated by β-tongue strands, might play a special role in the formation of Nhe trimers and Hbl trimers, respectively. The ‘cap’ of CytK, which includes two ‘latches’ with many β-sheets, formed a β-barrel structure with pores, and a ‘rim’ balanced the structure. The evolution of B. cereus strains showed a clustering tendency based on the protein-coding virulence genes, and the complete virulence-gene operon combination had higher relative genetic stability. The beta-tongue or latch associated with β-sheet folding might play an important role in the binding of virulence structures and pore-forming toxins in B. cereus.

Signatures of selection in core and accessory genomes indicate different ecological drivers of diversification among Bacillus cereus clades

Bacterial clades are often ecologically distinct, despite extensive horizontal gene transfer (HGT). How selection works on different parts of bacterial pan‐genomes to drive and maintain the emergence of clades is unclear. Focusing on the three largest clades in the diverse and well‐studied Bacillus cereus sensu lato group, we identified clade‐specific core genes (present in all clade members) and then used clade‐specific allelic diversity to identify genes under purifying and diversifying selection. Clade‐specific accessory genes (present in a subset of strains within a clade) were characterized as being under selection using presence/absence in specific clades. Gene ontology analyses of genes under selection revealed that different gene functions were enriched in different clades. Furthermore, some gene functions were enriched only amongst clade‐specific core or accessory genomes. Genes under purifying selection were often clade‐specific, while genes under diversifying selection showed signs of frequent HGT. These patterns are consistent with different selection pressures acting on both the core and the accessory genomes of different clades and can lead to ecological divergence in both cases. Examining variation in allelic diversity allows us to uncover genes under clade‐specific selection, allowing ready identification of strains and their ecological niche.

Three Novel Bacteriophages, J5a, F16Ba, and z1a, Specific for Bacillus anthracis, Define a New Clade of Historical Wbeta Phage Relatives

Bacillus anthracis is a potent biowarfare agent, able to be highly lethal. The bacteria dwell in the soil of certain regions, as natural flora. Bacteriophages or their lytic enzymes, endolysins, may be an alternative for antibiotics and other antibacterials to fight this pathogen in infections and to minimize environmental contamination with anthrax endospores. Upon screening environmental samples from various regions in Poland, we isolated three new siphophages, J5a, F16Ba, and z1a, specific for B. anthracis. They represent new species related to historical anthrax phages Gamma, Cherry, and Fah, and to phage Wbeta of Wbetavirus genus. We show that the new phages and their closest relatives, phages Tavor_SA, Negev_SA, and Carmel_SA, form a separate clade of the Wbetavirus genus, designated as J5a clade. The most distinctive feature of J5a clade phages is their cell lysis module. While in the historical phages it encodes a canonical endolysin and a class III holin, in J5a clade phages it encodes an endolysin with a signal peptide and two putative holins. We present the basic characteristic of the isolated phages. Their comparative genomic analysis indicates that they encode two receptor-binding proteins, of which one may bind a sugar moiety of B. anthracis cell surface.

Pathogenic Impacts of Bacillus cereus Strains on Crassostrea gigas

Regarding the economic importance of bivalve farming, a great deal of interest has recently been devoted to studying the pathogenesis of infectious diseases of these mollusks to prepare for public health emergencies. Bacillus cereus is one of these pathogens; it is a ubiquitous soil bacterium responsible for many types of gastrointestinal diseases associated with food. This study was conducted to determine the pathogenic effect of B. cereus on Crassostrea gigas. This effect was studied by assessing hemocytes death using flow cytometry analysis. The results showed that only ∼15% of C. gigas were able to survive after B. cereus artificial infection with 10⁸ CFU (colony-forming unit)/oyster. Evenly, the percentage of nonviable hemocytes gradually increased with the concentration of B. cereus, with a peak value of ∼40% after infection. Indeed, findings showed that this strain is harmful to C. gigas.

Microbial evolution and transitions along the parasite-mutualist continuum

Virtually all plants and animals, including humans, are home to symbiotic microorganisms. Symbiotic interactions can be neutral, harmful or have beneficial effects on the host organism. However, growing evidence suggests that microbial symbionts can evolve rapidly, resulting in drastic transitions along the parasite-mutualist continuum. In this Review, we integrate theoretical and empirical findings to discuss the mechanisms underpinning these evolutionary shifts, as well as the ecological drivers and why some host-microorganism interactions may be stuck at the end of the continuum. In addition to having biomedical consequences, understanding the dynamic life of microorganisms reveals how symbioses can shape an organism's biology and the entire community, particularly in a changing world.

Exploiting targeted and untargeted approaches for the analysis of bacterial metabolites under altered growth conditions

In the host, pathogenic microorganisms have developed stress responses to cope with constantly changing environments. Stress responses are directly related to changes in several metabolomic pathways, which could hamper microorganisms' unequivocal identification. We evaluated the effect of various in vitro stress conditions (acidic, basic, oxidative, ethanolic, and saline conditions) on the metabolism of Staphylococcus aureus, Bacillus cereus, and Pseudomonas aeruginosa, which are common lung pathogens. The metabolite profiles of the bacteria were analyzed using liquid chromatography coupled to triple quadrupole and quadrupole time-of-flight mass spectrometry. The advantages of targeted and untargeted analysis combined with univariate and multivariate statistical analysis (principal component analysis, hierarchical cluster analysis, partial least square discriminant analysis, random forest) were combined to unequivocally identify bacterial species. In normal in vitro conditions, the targeted methodology, based on the analysis of primary metabolites, enabled the rapid and efficient discrimination of the three bacteria. In changing in vitro conditions and specifically in presence of the various stressors, the untargeted methodology proved to be more valuable for the global and accurate differentiation of the three bacteria, also considering the type of stress environment within each species. In addition, species-specific metabolites (i.e., fatty acids, polysaccharides, peptides, and nucleotide bases derivatives) were putatively identified. Good intra-day repeatability and inter-day repeatability (< 10% RSD and < 15% RSD, respectively) were obtained for the targeted and the untargeted methods. This untargeted approach highlights its importance in unusual (and less known) bacterial growth environments, being a powerful tool for infectious disease diagnosis, where the accurate classification of microorganisms is sought.

New genetic biomarkers to differentiate non-pathogenic from clinically relevant Bacillus cereus strains

OBJECTIVES

Bacillus cereus is responsible for food poisoning and rare but severe clinical infections. The pathogenicity of strains varies from harmless to lethal strains. However, there are currently no markers, either alone or in combination, to differentiate pathogenic from non-pathogenic strains. The objective of the study was to identify new genetic biomarkers to differentiate non-pathogenic from clinically relevant B. cereus strains.

METHODS

A first set of 15 B. cereus strains were compared by RNAseq. A logistic regression model with lasso penalty was applied to define combination of genes whose expression was associated with strain pathogenicity. The identified markers were checked for their presence/absence in a collection of 95 B. cereus strains with varying pathogenic potential (food-borne outbreaks, clinical and non-pathogenic). Receiver operating characteristic area under the curve (AUC) analysis was used to determine the combination of biomarkers, which best differentiate between the "disease" versus "non-disease" groups.

RESULTS

Seven genes were identified during the RNAseq analysis with a prediction to differentiate between pathogenic and non-pathogenic strains. The validation of the presence/absence of these genes in a larger collection of strains coupled with AUC prediction showed that a combination of four biomarkers was sufficient to accurately discern clinical strains from harmless strains, with an AUC of 0.955, sensitivity of 0.9 and specificity of 0.86.

CONCLUSIONS

These new findings help in the understanding of B. cereus pathogenic potential and complexity and may provide tools for a better assessment of the risks associated with B. cereus contamination to improve patient health and food safety.

The Food Poisoning Toxins of Bacillus cereus

Bacillus cereus is a ubiquitous soil bacterium responsible for two types of food-associated gastrointestinal diseases. While the emetic type, a food intoxication, manifests in nausea and vomiting, food infections with enteropathogenic strains cause diarrhea and abdominal pain. Causative toxins are the cyclic dodecadepsipeptide cereulide, and the proteinaceous enterotoxins hemolysin BL (Hbl), nonhemolytic enterotoxin (Nhe) and cytotoxin K (CytK), respectively. This review covers the current knowledge on distribution and genetic organization of the toxin genes, as well as mechanisms of enterotoxin gene regulation and toxin secretion. In this context, the exceptionally high variability of toxin production between single strains is highlighted. In addition, the mode of action of the pore-forming enterotoxins and their effect on target cells is described in detail. The main focus of this review are the two tripartite enterotoxin complexes Hbl and Nhe, but the latest findings on cereulide and CytK are also presented, as well as methods for toxin detection, and the contribution of further putative virulence factors to the diarrheal disease.

Multi-locus sequence typing and virulence profile in Bacillus cereus sensu lato strains isolated from dairy products

Members of Bacillus cereus group are important food contaminants and they are of relevant interest in food safety and public heath due to their ability to cause two distinct forms of food poisoning, emetic and diarrhoeal syndrome. In the present study, 90 strains of B. cereus isolated from dairy products, have been typed using Multilocus Sequence Typing (MLST) analysis and investigated for the occurrence of 10 enterotoxigenic genes (hblA, hblC, hblD, nheA, nheB, nheC, cytK, entFM, entS and bceT) and one emetogenic gene (ces), to determine their genetic diversity. A total of 58 sequence types were identified and among these 17 were signalled as new profiles. Among the virulence genes, the majority of our strains carried the entS (92%), entFM (86%), nhe (82%) and cytK (72%) genes. All remaining genes were identified in at least one strain with different prevalence, stressing the genetic diversity, how even the different grade of pathogenicity of B. cereus isolated from dairy products.

Divergence in environmental adaptation between terrestrial clades of the Bacillus cereus group

The Bacillus cereus group encompasses beneficial and harmful species in diverse niches and has a much debated taxonomy. Investigating whether selection has led to ecological divergence between phylogenetic clades can help understand the basis of speciation, and has implications for predicting biological safety across this group. Using three most terrestrial species in this group (B. cereus, Bacillus thuringiensis and Bacillus mycoides) we charactererized ecological specialization in terms of resource use, thermal adaptation and fitness in different environmental conditions and tested whether taxonomic species or phylogenetic clade best explained phenotypic variation. All isolates grew vigorously in protein rich media and insect cadavers, but exploitation of soil or plant derived nutrients was similarly weak for all. For B. thuringiensis and B. mycoides, clade and taxonomic species were important predictors of relative fitness in insect infections. Fully psychrotolerant isolates could outcompete B. thuringiensis in insects at low temperature, although psychrotolerance predicted growth in artificial media better than clade. In contrast to predictions, isolates in the Bacillus anthracis clade had sub-optimal growth at 37°C. The common ecological niche in these terrestrial B. cereus species is the ability to exploit protein rich resources such as cadavers. However, selection has led to different phylogenetic groups developing different strategies for accessing this resource. Thus, clades, as well as traditional taxonomic phenotypes, predict biologically important traits.

The assessment of leading traits in the taxonomy of the Bacillus cereus group

Bacillus cereus sensu lato strains (B. cereus group) are widely distributed in nature and have received interest for decades due to their importance in insect pest management, food production and their positive and negative repercussions in human health. Consideration of practical uses such as virulence, physiology, morphology, or ill-defined features have been applied to describe and classify species of the group. However, current comparative studies have exposed inconsistencies between evolutionary relatedness and biological significance among genomospecies of the B. cereus group. Here, the combined analyses of core-based phylogeny and all versus all Average Nucleotide Identity values based on 2116 strains were conducted to update the genomospecies circumscriptions within B. cereus group. These analyses suggested the existence of 57 genomospecies, 37 of which are novel, thus indicating that the taxonomic identities of more than 39% of the analyzed strains should be revised or updated. In addition, we found that whole-genome in silico analyses were suitable to differentiate genomospecies such as B. anthracis, B. cereus and B. thuringiensis. The prevalence of toxin and virulence factors coding genes in each of the genomospecies of the B. cereus group was also examined, using phylogeny-aware methods at wide-genome scale. Remarkably, Cry and emetic toxins, commonly assumed to be associated with B. thuringiensis and emetic B. paranthracis, respectively, did not show a positive correlation with those genomospecies. On the other hand, anthrax-like toxin and capsule-biosynthesis coding genes were positively correlated with B. anthracis genomospecies, despite not being present in all strains, and with presumably non-pathogenic genomospecies. Hence, despite these features have been so far considered relevant for industrial or medical classification of related species of the B. cereus group, they were inappropriate for their circumscription. In this study, genomospecies of the group were accurately affiliated and representative strains defined, generating a rational framework that will allow comparative analysis in epidemiological or ecological studies. Based on this classification the role of specific markers such as Type VII secretion system, cytolysin, bacillolysin, and siderophores such as petrobactin were pointed out for further analysis.

Genotypic Profiling of Bacillus cereus Recovered from Some Retail Foods in Ogun State, Nigeria, and Their Phylogenetic Relationship

Identifying Bacillus cereus with conventional methods is neither specific nor rapid because of the close relatedness of the B. cereus group, hence the need for molecular methods. Genotypic profiling of B. cereus isolates from food was obtained by Random Amplified Polymorphic DNA-polymerase chain reaction (RAPD-PCR) using OPR13 primer. A dendrogram was drawn with the Numerical Taxonomy System of Statistic (NTSYS) software. Thirty of the isolates were subjected to molecular identification by 16S rDNA sequencing. The thirty sequences were deposited in GenBank for accession number. Phylogenetic relationship of the 16S rDNA sequence obtained was carried out with the Multiple Alignment using Fast Fourier Transform (MAFFT) software version 7.0. The evolutionary tree was drawn using the Molecular Evolutionary Genetics Analysis (MEGA 6) software. The dendrogram generated for the RAPD profile showed that all the strains are closely related, with a similarity coefficient of 70%. The isolates were confirmed with 16S rDNA sequencing as B. cereus. The thirty sequences deposited in GenBank were given accession numbers: KX574760–KX574769, KX610811–KX610820, MT757957-MT757963, and MT772282-MT772284. By comparing the phylogenetic relationship, eleven of the strains did not cluster with the reference strains from the GenBank but form distinct clades, which means they are likely to be of different ancestors. Conventional methods rarely differentiate bacteria of the same species into clade, neither can it describe their ancestral lineage. Therefore, it is important to employ molecular methods in identifying bacteria to give detailed information about them.

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.

Bacillus paranthracis Isolate from Blood of Fatal Ebola Virus Disease Case

A Bacillus paranthracis isolate was cultured from the blood of a fatal Ebola virus disease (EVD) case in Liberia and was identified by whole genome sequencing. Although B. paranthracis has only recently been described and is poorly characterized, this case may represent the bacterial co-infection of an EVD patient.

The Population Structure of Pseudomonas aeruginosa Is Characterized by Genetic Isolation of exoU+ and exoS+ Lineages

The diversification of microbial populations may be driven by many factors including adaptation to distinct ecological niches and barriers to recombination. We examined the population structure of the bacterial pathogen Pseudomonas aeruginosa by analyzing whole-genome sequences of 739 isolates from diverse sources. We confirmed that the population structure of P. aeruginosa consists of two major groups (referred to as Groups A and B) and at least two minor groups (Groups C1 and C2). Evidence for frequent intragroup but limited intergroup recombination in the core genome was observed, consistent with sexual isolation of the groups. Likewise, accessory genome analysis demonstrated more gene flow within Groups A and B than between these groups, and a few accessory genomic elements were nearly specific to one or the other group. In particular, the exoS gene was highly overrepresented in Group A compared with Group B isolates (99.4% vs. 1.1%) and the exoU gene was highly overrepresented in Group B compared with Group A isolates (95.2% vs. 1.8%). The exoS and exoU genes encode effector proteins secreted by the P. aeruginosa type III secretion system. Together these results suggest that the major P. aeruginosa groups defined in part by the exoS and exoU genes are divergent from each other, and that these groups are genetically isolated and may be ecologically distinct. Although both groups were globally distributed and caused human infections, certain groups predominated in some clinical contexts.

Comparative genomic survey of Bacillus cereus sensu stricto isolates from the dairy production chain in Brazil

The genomes of 262 Bacillus cereus isolates were analyzed including 69 isolates sampled from equipment, raw milk and dairy products from Brazil. The population structure of isolates showed strains belonging to known phylogenetic groups II, III, IV, V and VI. Almost all the isolates obtained from dairy products belonged to group III. Investigation of specific alleles revealed high numbers of isolates carrying toxin-associated genes including cytK (53.62%), hblA (59.42%), hblC (44.93%), hblD (53.62%), nheA (84.06%), nheB (89.86%) and nheC (84.06%) with isolates belonging to groups IV and V having significant higher prevalence of hblACD and group IV of CytK genes. Strains from dairy products had significantly lower prevalence of CytK and hblACD genes compared to isolates from equipment and raw milk/bulk tanks. Genes related to sucrose metabolism were detected at higher frequency in isolates obtained from raw milk compared to strains from equipment and utensils. The population genomic analysis demonstrated the diversity of strains and variability of putative function among B. cereus group isolates in Brazilian dairy production, with large numbers of strains potentially able to cause foodborne illness. This detailed information will contribute to targeted interventions to reduce milk contamination and spoilage associated with B. cereus in Brazil.

Diversity of Bacillus cereus sensu lato mobilome

BACKGROUND

Bacillus cereus sensu lato s.l.) is a group of bacteria displaying close phylogenetic relationships but a high ecological diversity. The three most studied species are Bacillus anthracis, Bacillus cereus sensu stricto and Bacillus thuringiensis. While some species are pathogenic to mammals or associated with food poisoning, Bacillus thuringiensis is a well-known entomopathogenic bacterium used as biopesticide worldwide. B. cereus s.l. also contains a large variety of mobile genetic elements (MGEs).

RESULTS

In this study, we detail the occurrence and plasmid vs. chromosome distribution of several MGEs in 102 complete and annotated genomes of B. cereus s.l. These MGEs include 16 Insertion Sequence (IS) families, the Tn3 family, 18 different Bacillus cereus repeats (BCRs) and 30 known group II introns.

CONCLUSIONS

Our analysis not only shows the diversity of these MGEs among strains of the same species and between different species within the B. cereus s.l. group, but also highlights the potential impact of these elements on the plasticity of the plasmid pool, and the TEs (Transposable Elements) - species relationship within B. cereus s.l.

Bactéries sporulées et biofilms : un problème récurrent dans les lignes de production de lait reconstitué ou recombiné pasteurisé

In the dairy industry, bacterial contaminants persist on equipment surfaces due to spore and biofilm formation. These are involved in cross-contamination problems that affect the quality of processed products and limit their life. The pasteurization process, in which milk is submitted to moderate heat treatment, is inefficient against bacterial spores. The most prevalent sporulated bacteria belong to Bacillus and related genera. The situation is more complicated in countries where pasteurized milk is derived from imported milk powder originally contaminated by bacterial spores. Studies have shown biofilm formation on dairy equipment by mesophilic strains from the group Bacillus cereus and thermophilic strains from the genus Geobacillus. These biofilms are resistant to cleaning procedures and are sources of chronic contamination of pasteurized milk. This review analyzes the dairy situation in Algeria exposed to sporulated flora and derived biofilm problems, with the aim of proposing efficient solutions in the light of current knowledge. [Journal translation].

Genes Associated With Psychrotolerant Bacillus cereus Group Isolates

The Bacillus cereus group comprises 18 different species, including human pathogens as well as psychrotolerant strains that are an important cause of fluid milk spoilage. To enhance our understanding of the genetic markers associated with psychrotolerance (defined here as > 1 log₁₀ increase in cfu/mL after 21 days incubation at 6°C) among dairy-associated B. cereus group isolates, we used genetic (whole genome sequencing) and phenotypic methods [growth in Skim Milk Broth (SMB) and Brain Heart Infusion (BHI) broth] to characterize 23 genetically-distinct representative isolates from a collection of 503 dairy-associated isolates. Quality threshold clustering identified three categories of psychrotolerance: (i) 14 isolates that were not psychrotolerant in BHI or SMB, (ii) 6 isolates that were psychrotolerant in BHI but not in SMB, and (iii) 2 isolates that were psychrotolerant in BHI and SMB. One isolate, which was psychrotolerant in BHI broth but was just below the cut-off of >1 log₁₀ cfu/mL increase in SMB was not assigned to a cluster. A maximum likelihood phylogeny constructed with core genome single nucleotide polymorphisms classified all psychrotolerant isolates (i.e., psychrotolerant in BHI) into clade VI (representing B. mycoides/weihenstephanensis). Analysis of correlations between gene ortholog presence or absence patterns and psychrotolerance identified 206 orthologous gene clusters that were significantly overrepresented among psychrotolerant strains, including two clusters of cold shock proteins, which were identified in 8/9 and 7/9 psychrotolerant isolates. Gene ontology analyses revealed 36 gene ontology terms that were overrepresented in psychrotolerant isolates, including putrescine catabolic processes and putrescine transmembrane transporter activity. Lastly, Hidden Markov Model searches identified three protein family motifs, including cold shock domain proteins and fatty acid hydroxylases that were significantly associated with psychrotolerance in BHI broth. Analyses of CspA sequences revealed a positive association between psychrotolerant strains and a previously identified “psychrotolerant” CspA sequence. Overall, our data highlight genetic and phenotypic differences in psychrotolerance among B. cereus group dairy-associated isolates and show that psychrotolerance is dependent on the growth medium. We also identified a number of gene targets that could be used for specific detection or control of psychrotolerant B. cereus group isolates.

Pathogenicity, population genetics and dissemination of Bacillus anthracis

Bacillus anthracis, the etiological agent of anthrax, procures its particular virulence by a capsule and two AB type toxins: the lethal factor LF and the edema factor EF. These toxins primarily disable immune cells. Both toxins are translocated to the host cell by the adhesin-internalin subunit called protective antigen PA. PA enables LF to reach intra-luminal vesicles, where it remains active for long periods. Subsequently, LF translocates to non-infected cells, leading to inefficient late therapy of anthrax. B. anthracis undergoes slow evolution because it alternates between vegetative and long spore phases. Full genome sequence analysis of a large number of worldwide strains resulted in a robust evolutionary reconstruction of this bacterium, showing that B. anthracis is split in three main clades: A, B and C. Clade A efficiently disseminated worldwide underpinned by human activities including heavy intercontinental trade of goat and sheep hair. Subclade A.Br.WNA, which is widespread in the Northern American continent, is estimated to have split from clade A reaching the Northern American continent in the late Pleistocene epoch via the former Bering Land Bridge and further spread from Northwest southwards. An alternative hypothesis is that subclade A.Br.WNA. evolved from clade A.Br.TEA tracing it back to strains from Northern France that were assumingly dispatched by European explorers that settled along the St. Lawrence River. Clade B established mostly in Europe along the alpine axis where it evolved in association with local cattle breeds and hence displays specific geographic subclusters. Sequencing technologies are also used for forensic applications to trace unintended or criminal acts of release of B. anthracis. Under natural conditions, B. anthracis generally affects domesticated and wild ruminants in arid ecosystems. The more recently discovered B. cereus biovar anthracis spreads in tropical forests, where it threatens particularly endangered primate populations.

Lineage-specific plasmid acquisition and the evolution of specialized pathogens in Bacillus thuringiensis and the Bacillus cereus group

Bacterial plasmids can vary from small selfish genetic elements to large autonomous replicons that constitute a significant proportion of total cellular DNA. By conferring novel function to the cell, plasmids may facilitate evolution but their mobility may be opposed by co-evolutionary relationships with chromosomes or encouraged via the infectious sharing of genes encoding public goods. Here, we explore these hypotheses through large-scale examination of the association between plasmids and chromosomal DNA in the phenotypically diverse Bacillus cereus group. This complex group is rich in plasmids, many of which encode essential virulence factors (Cry toxins) that are known public goods. We characterized population genomic structure, gene content and plasmid distribution to investigate the role of mobile elements in diversification. We analysed coding sequence within the core and accessory genome of 190 B. cereus group isolates, including 23 novel sequences and genes from 410 reference plasmid genomes. While cry genes were widely distributed, those with invertebrate toxicity were predominantly associated with one sequence cluster (clade 2) and phenotypically defined Bacillus thuringiensis. Cry toxin plasmids in clade 2 showed evidence of recent horizontal transfer and variable gene content, a pattern of plasmid segregation consistent with transfer during infectious cooperation. Nevertheless, comparison between clades suggests that co-evolutionary interactions may drive association between plasmids and chromosomes and limit wider transfer of key virulence traits. Proliferation of successful plasmid and chromosome combinations is a feature of specialized pathogens with characteristic niches (Bacillus anthracis, B. thuringiensis) and has occurred multiple times in the B. cereus group.

The Replacement of five Consecutive Amino Acids in the Cyt1A Protein of Bacillus thuringiensis Enhances its Cytotoxic Activity against Lung Epithelial Cancer Cells

Cyt1A protein is a cytolytic protein encoded by the cyt gene of Bacillus thuringiensis subsp. israelensis (Bti) as part of the parasporal crystal proteins produced during the sporulation. Cyt1A protein is unique compared to the other endotoxins present in these parasporal crystals. Unlike δ-endotoxins, Cyt1A protein does not require receptors to bind to the target cell and activate the toxicity. It has the ability to affect a broad range of cell types and organisms, due to this characteristic. Cyt1A has been recognized to not only target the insect cells directly, but also recruit other endotoxins by acting as receptors. Due to these mode of actions, Cyt1A has been studied for its cytolytic activity against human cancer cell lines, although not extensively. In this study, we report a novel Cyt1A protein produced by a Bti strain QBT229 isolated from Qatar. When tested for its cytotoxicity against lung cancer cells, this local strain showed considerably higher activity compared to that of the reference Bti and other strains tested. The possible reasons for such enhanced activity were explored at the gene and protein levels. It was evidenced that five consecutive amino acid replacements in the β8 sheet of the Cyt1A protein enhanced the cytotoxicity against the lung epithelial cancer cells. Such novel Cyt1A protein with high cytotoxicity against lung cancer cells has been characterized and reported through this study.

Subterranean Mammals: Reservoirs of Infection or Overlooked Sentinels of Anthropogenic Environmental Soiling?

Global reports of emergent pathogens in humans have intensified efforts to identify wildlife reservoirs. Subterranean mammals, such as bathyergid mole rats, are largely overlooked, despite their high-level exposure to soil-dwelling microbes. Initial assessment of bathyergid reservoir potential was determined using a broad-range 16S rRNA PCR approach, which revealed an 83% PCR-positivity for the 234 bathyergid lung samples evaluated. The presence of the Bacillus cereus complex, a ubiquitous bacterial assemblage, containing pathogenic and zoonotic species, was confirmed through nucleotide sequencing, prior to group- and species-specific PCR sequencing. The latter allowed for enhanced placement and prevalence estimations of Bacillus in four bathyergid species sampled across a range of transformed landscapes in the Western Cape Province, South Africa. Two novel Bacillus strains (1 and 2) identified on the basis of the concatenated 16S rRNA-groEL-yeaC data set (2066 nucleotides in length), clustered with B. mycoides (ATCC 6462) and B. weihenstephanensis (WSBC 10204), within a well-supported monophyletic lineage. The levels of co-infection, evaluated with a groEL strain-specific assay, developed specifically for this purpose, were high (71%). The overall Bacillus presence of 17.95% (ranging from 0% for Georychus capensis to 45.35% for Bathyergus suillus) differed significantly between host species (χ² = 69.643; df = 3; P < 0.05), being significantly higher in bathyergids sampled near an urban informal settlement (χ² = 70.245; df = 3; P < 0.05). The results highlight the sentinel potential of soil-dwelling mammals for monitoring anthropogenically introduced, opportunistic pathogens and the threats they pose to vulnerable communities, particularly in the developing world.

Identification and genomic comparison of temperate bacteriophages derived from emetic Bacillus cereus

Cereulide-producing Bacillus cereus isolates can cause serious emetic (vomiting) syndrome and even acute lethality. As mobile genetic elements, the exploration of prophages derived from emetic B. cereus isolates will help in our understanding of the genetic diversity and evolution of these pathogens. In this study, five temperate phages derived from cereulide-producing B. cereus strains were induced, with four of them undergoing genomic sequencing. Sequencing revealed that they all belong to the Siphoviridae family, but presented in different forms in their hosts. PfNC7401 and PfIS075 have typical icosahedral heads, probably existing alone as phagemids in the host with self-replicating capability in the lysogenic state. PfEFR-4, PfEFR-5, and PfATCC7953 have elongated heads, with the genomes of the former two identified as linear dsDNA, which could be integrated into the host genome during the lysogenic state. Genomic comparison of the four phages with others also derived from emetic B. cereus isolates showed similar genome structures and core genes, thus displaying host spectrum specificity. In addition, phylogenic analysis based on the complete genome and conserved tail fiber proteins of 36 Bacillus species-derived phages confirmed that the phages derived from emetic B. cereus strains were highly similar. Furthermore, one endolysin LysPfEFR-4 was cloned and showed lytic activity against all tested emetic B. cereus strains and cross-lytic activity against some other pathogenic bacteria, implying a potential to control bacterial contamination in the food supply.

Comparative Genomics of Bacillus thuringiensis Reveals a Path to Specialized Exploitation of Multiple Invertebrate Hosts

Understanding the genetic basis of host shifts is a key genomic question for pathogen and parasite biology. The Bacillus cereus group, which encompasses Bacillus thuringiensis and Bacillus anthracis, contains pathogens that can infect insects, nematodes, and vertebrates. Since the target range of the essential virulence factors (Cry toxins) and many isolates is well known, this group presents a powerful system for investigating how pathogens can diversify and adapt to phylogenetically distant hosts. Specialization to exploit insects occurs at the level of the major clade and is associated with substantial changes in the core genome, and host switching between insect orders has occurred repeatedly within subclades. The transfer of plasmids with linked cry genes may account for much of the adaptation to particular insect orders, and network analysis implies that host specialization has produced strong associations between key toxin genes with similar targets. Analysis of the distribution of plasmid minireplicons shows that plasmids with orf156 and orf157, which carry genes encoding toxins against Lepidoptera or Diptera, were contained only by B. thuringiensis in the specialized insect clade (clade 2), indicating that tight genome/plasmid associations have been important in adaptation to invertebrate hosts. Moreover, the accumulation of multiple virulence factors on transposable elements suggests that cotransfer of diverse virulence factors is advantageous in terms of expanding the insecticidal spectrum, overcoming insect resistance, or through gains in pathogenicity via synergistic interactions between toxins.IMPORTANCE Population genomics have provided many new insights into the formation, evolution, and dynamics of bacterial pathogens of humans and other higher animals, but these pathogens usually have very narrow host ranges. As a pathogen of insects and nematodes, Bacillus thuringiensis, which produces toxins showing toxicity to many orders of insects and other invertebrates, can be used as a model to study the evolution of pathogens with wide host ranges. Phylogenomic analysis revealed that host specialization and switching occur at the level of the major clade and subclade, respectively. A toxin gene co-occurrence network indicates that multiple toxins with similar targets were accumulated by the same cell in the whole species. This accumulation may be one of the strategies that B. thuringiensis has used to fight against host resistance. This kind of formation and evolution of pathogens represents a different path used against multiple invertebrate hosts from that used against higher animals.

Pan-genome and phylogeny of Bacillus cereus sensu lato

BACKGROUND

Bacillus cereus sensu lato (s. l.) is an ecologically diverse bacterial group of medical and agricultural significance. In this study, I use publicly available genomes and novel bioinformatic workflows to characterize the B. cereus s. l. pan-genome and perform the largest phylogenetic and population genetic analyses of this group to date in terms of the number of genes and taxa included. With these fundamental data in hand, I identify genes associated with particular phenotypic traits (i.e., "pan-GWAS" analysis), and quantify the degree to which taxa sharing common attributes are phylogenetically clustered.

METHODS

A rapid k-mer based approach (Mash) was used to create reduced representations of selected Bacillus genomes, and a fast distance-based phylogenetic analysis of this data (FastME) was performed to determine which species should be included in B. cereus s. l. The complete genomes of eight B. cereus s. l. species were annotated de novo with Prokka, and these annotations were used by Roary to produce the B. cereus s. l. pan-genome. Scoary was used to associate gene presence and absence patterns with various phenotypes. The orthologous protein sequence clusters produced by Roary were filtered and used to build HaMStR databases of gene models that were used in turn to construct phylogenetic data matrices. Phylogenetic analyses used RAxML, DendroPy, ClonalFrameML, PAUP*, and SplitsTree. Bayesian model-based population genetic analysis assigned taxa to clusters using hierBAPS. The genealogical sorting index was used to quantify the phylogenetic clustering of taxa sharing common attributes.

RESULTS

The B. cereus s. l. pan-genome currently consists of ≈60,000 genes, ≈600 of which are "core" (common to at least 99% of taxa sampled). Pan-GWAS analysis revealed genes associated with phenotypes such as isolation source, oxygen requirement, and ability to cause diseases such as anthrax or food poisoning. Extensive phylogenetic analyses using an unprecedented amount of data produced phylogenies that were largely concordant with each other and with previous studies. Phylogenetic support as measured by bootstrap probabilities increased markedly when all suitable pan-genome data was included in phylogenetic analyses, as opposed to when only core genes were used. Bayesian population genetic analysis recommended subdividing the three major clades of B. cereus s. l. into nine clusters. Taxa sharing common traits and species designations exhibited varying degrees of phylogenetic clustering.

CONCLUSIONS

All phylogenetic analyses recapitulated two previously used classification systems, and taxa were consistently assigned to the same major clade and group. By including accessory genes from the pan-genome in the phylogenetic analyses, I produced an exceptionally well-supported phylogeny of 114 complete B. cereus s. l. genomes. The best-performing methods were used to produce a phylogeny of all 498 publicly available B. cereus s. l. genomes, which was in turn used to compare three different classification systems and to test the monophyly status of various B. cereus s. l. species. The majority of the methodology used in this study is generic and could be leveraged to produce pan-genome estimates and similarly robust phylogenetic hypotheses for other bacterial groups.

In defense of Bacillus thuringiensis, the safest and most successful microbial insecticide available to humanity - a response to EFSA

The Bacillus cereus group contains vertebrate pathogens such as B. anthracis and B. cereus and the invertebrate pathogen B. thuringiensis (Bt). Microbial biopesticides based on Bt are widely recognised as being among the safest and least environmentally damaging insecticidal products available. Nevertheless, a recent food-poisoning incident prompted a European Food Safety Authority review which argued that Bt poses a health risk equivalent to B. cereus, a causative agent of diarrhoea. However, a critical examination of available data, and this latest incident, provides no solid evidence that Bt causes diarrhoea. Although relatively high levels of B. cereus-like spores can occur in foods, genotyping demonstrates that these are predominantly naturally occurring strains rather than biopesticides. Moreover, MLST genotyping of >2000 isolates show that biopesticide genotypes have never been isolated from any clinical infection. MLST data demonstrate that B. cereus group is heterogeneous and formed of distinct clades with substantial differences in biology, ecology and host association. The group posing the greatest risk (the anthracis clade) is distantly related to the clade containing all biopesticides. These recent data support the long-held view that Bt and especially the strains used in Bt biopesticides are very safe for humans.

The Bacterial Sequential Markov Coalescent

Bacteria can exchange and acquire new genetic material from other organisms directly and via the environment. This process, known as bacterial recombination, has a strong impact on the evolution of bacteria, for example, leading to the spread of antibiotic resistance across clades and species, and to the avoidance of clonal interference. Recombination hinders phylogenetic and transmission inference because it creates patterns of substitutions (homoplasies) inconsistent with the hypothesis of a single evolutionary tree. Bacterial recombination is typically modeled as statistically akin to gene conversion in eukaryotes, i.e., using the coalescent with gene conversion (CGC). However, this model can be very computationally demanding as it needs to account for the correlations of evolutionary histories of even distant loci. So, with the increasing popularity of whole genome sequencing, the need has emerged for a faster approach to model and simulate bacterial genome evolution. We present a new model that approximates the coalescent with gene conversion: the bacterial sequential Markov coalescent (BSMC). Our approach is based on a similar idea to the sequential Markov coalescent (SMC)-an approximation of the coalescent with crossover recombination. However, bacterial recombination poses hurdles to a sequential Markov approximation, as it leads to strong correlations and linkage disequilibrium across very distant sites in the genome. Our BSMC overcomes these difficulties, and shows a considerable reduction in computational demand compared to the exact CGC, and very similar patterns in simulated data. We implemented our BSMC model within new simulation software FastSimBac. In addition to the decreased computational demand compared to previous bacterial genome evolution simulators, FastSimBac provides more general options for evolutionary scenarios, allowing population structure with migration, speciation, population size changes, and recombination hotspots. FastSimBac is available from https://bitbucket.org/nicofmay/fastsimbac, and is distributed as open source under the terms of the GNU General Public License. Lastly, we use the BSMC within an Approximate Bayesian Computation (ABC) inference scheme, and suggest that parameters simulated under the exact CGC can correctly be recovered, further showcasing the accuracy of the BSMC. With this ABC we infer recombination rate, mutation rate, and recombination tract length of Bacillus cereus from a whole genome alignment.

Simulating Intestinal Growth Conditions Enhances Toxin Production of Enteropathogenic Bacillus cereus

Bacillus cereus is a ubiquitous bacterial pathogen increasingly reported to be the causative agent of foodborne infections and intoxications. Since the enterotoxins linked to the diarrheal form of food poising are foremost produced in the human intestine, the toxic potential of enteropathogenic B. cereus strains is difficult to predict from studies carried out under routine cultivation procedures. In this study, toxigenic properties of a panel of strains (n = 19) of diverse origin were compared using cell culture medium pre-incubated with CaCo-2 cells to mimic intestinal growth conditions. Shortly after contact of the bacteria with the simulated host environment, enterotoxin gene expression was activated and total protein secretion of all strains was accelerated. Although the signal stimulating enterotoxin production still needs to be elucidated, it could be shown that it originated from the CaCo-2 cells. Overall, our study demonstrates that the currently used methods in B. cereus diagnostics, based on standard culture medium, are not allowing a conclusive prediction of the potential health risk related to a certain strain. Thus, these methods should be complemented by cultivation procedures that are simulating intestinal host conditions.

Genetic diversity and population structure of the Bacillus cereus group bacteria from diverse marine environments

The phylogenetic diversity of marine bacteria belonged to the Bacillus cereus group has not been well investigated. Here, we present the genetic diversity and population structure of 71 bacteria from diverse marine environments, using a multilocus sequence typing (MLST) approach and the analyses of digital DNA-DNA hybridization (dDDH) and average nucleotide identity (ANI) based on some representative genomic sequences. The MLST analysis demonstrated that these isolates were highly diverse and a wide distribution in marine environments and some of them showed niche specificity to some extent. They were assigned to 27 sequence types (STs) with 23 novel STs. Phylogenetic analysis of 82 bacteria containing 11 type strains based on MLST discriminated them as 20 clusters including 10 new ones. Both the dDDH and ANI results supported the proposition that each of 20 clusters represented one independent species, including 10 putative novel species. Values of 98.3% of MLST similarity and 96.2% of ANI were proposed as the standard for the species definition of this group. In summary, the first insight into the phylogenetic diversity of the group bacteria from marine environments will contribute to better understanding of their ecological role and evolution in contrast with terrestrial environments.

Is Predominant Clonal Evolution a Common Evolutionary Adaptation to Parasitism in Pathogenic Parasitic Protozoa, Fungi, Bacteria, and Viruses?

We propose that predominant clonal evolution (PCE) in microbial pathogens be defined as restrained recombination on an evolutionary scale, with genetic exchange scarce enough to not break the prevalent pattern of clonal population structure. The main features of PCE are (1) strong linkage disequilibrium, (2) the widespread occurrence of stable genetic clusters blurred by occasional bouts of genetic exchange ('near-clades'), (3) the existence of a "clonality threshold", beyond which recombination is efficiently countered by PCE, and near-clades irreversibly diverge. We hypothesize that the PCE features are not mainly due to natural selection but also chiefly originate from in-built genetic properties of pathogens. We show that the PCE model obtains even in microbes that have been considered as 'highly recombining', such as Neisseria meningitidis, and that some clonality features are observed even in Plasmodium, which has been long described as panmictic. Lastly, we provide evidence that PCE features are also observed in viruses, taking into account their extremely fast genetic turnover. The PCE model provides a convenient population genetic framework for any kind of micropathogen. It makes it possible to describe convenient units of analysis (clones and near-clades) for all applied studies. Due to PCE features, these units of analysis are stable in space and time, and clearly delimited. The PCE model opens up the possibility of revisiting the problem of species definition in these organisms. We hypothesize that PCE constitutes a major evolutionary strategy for protozoa, fungi, bacteria, and viruses to adapt to parasitism.

Genotypic heterogeneity of emetic toxin producing Bacillus cereus isolates from China

Emetic toxin-producing Bacillus cereus (emetic B. cereus) is the third member of B. cereus group whose toxins are encoded by megaplasmids, beside anthrax and insecticidal toxins of B. anthracis and B. thuringiensis, respectively. A total of 18 emetic isolates collected from food poisoning events, clinical and non-random food samples in Zhejiang province of China, were analyzed by plasmid screening, pulse field gel electrophoresis, multilocus sequence typing, and toxic gene identification to investigate their genotypic diversity. In this study, 13 plasmid profile types, 14 pulse types and 6 different STs from emetic isolates were detected, in which ST 1035,1038,1053,1054 and 1065 were first assigned and reported. The toxic gene ces existed on its own, or coexisted with other toxic genes bceT, cytk, entFM and nhe, but never with hbl in emetic isolates. The results demonstrated that the emetic B. cereus strains from China were heterologous at genotypic level.

MALDI-TOF MS portrait of emetic and non-emetic Bacillus cereus group members

The number of foodborne intoxications caused by emetic Bacillus cereus isolates has increased significantly. As such, rapid and reliable methods to identify emetic strains appear to be clinically relevant. In this study, intact cell matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) was used to differentiate emetic and non-emetic bacilli. The phyloproteomic clustering of 34 B. cereus emetic and 88 non-emetic isolates classified as B. cereus, Bacillus thuringiensis, Bacillus weihenstephanensis, and Bacillus mycoides, showed (i) a clear separation of both groups at a similarity level of 43%, and (ii) a high relatedness among the emetic isolates (similarity of 78%). Specifically, 83 mass peak classes were recognized in the spectral window range between m/z 4000 and 12 000 that were tentatively assigned to 41 protein variants based on a bioinformatic approach. Mass variation between the emetic and the non-emetic subsets was recorded for 27 of them, including ten ribosomal subunit proteins, for which inter-strain polymorphism was confirmed by gene sequencing. Additional peaks were assigned to other proteins such as small acid soluble proteins, cold shock proteins and hypothetical proteins, e.g., carbohydrate kinase. Moreover, the results were supported by in silico analysis of the biomarkers in 259 members of B. cereus group, including Bacillus anthracis, based on their whole-genome sequences. In conclusion, the proteomic profiling by MALDI-TOF MS is a promising and rapid method for pre-screening B. cereus to identify medically relevant isolates and for epidemiologic purposes.

Feces Derived Allergens of Tyrophagus putrescentiae Reared on Dried Dog Food and Evidence of the Strong Nutritional Interaction between the Mite and Bacillus cereus Producing Protease Bacillolysins and Exo-chitinases

Tyrophagus putrescentiae (Schrank, 1781) is an emerging source of allergens in stored products and homes. Feces proteases are the major allergens of astigmatid mites (Acari: Acaridida). In addition, the mites are carriers of microorganisms and microbial adjuvant compounds that stimulate innate signaling pathways. We sought to analyze the mite feces proteome, proteolytic activities, and mite-bacterial interaction in dry dog food (DDF). Proteomic methods comprising enzymatic and zymographic analysis of proteases and 2D-E-MS/MS were performed. The highest protease activity was assigned to trypsin-like proteases; lower activity was assigned to chymotrypsin-like proteases, and the cysteine protease cathepsin B-like had very low activity. The 2D-E-MS/MS proteomic analysis identified mite trypsin allergen Tyr p3, fatty acid-binding protein Tyr p13 and putative mite allergens ferritin (Grp 30) and (poly)ubiquitins. Tyr p3 was detected at different positions of the 2D-E. It indicates presence of zymogen at basic pI, and mature-enzyme form and enzyme fragment at acidic pI. Bacillolysins (neutral and alkaline proteases) of Bacillus cereus symbiont can contribute to the protease activity of the mite extract. The bacterial exo-chitinases likely contribute to degradation of mite exuviae, mite bodies or food boluses consisting of chitin, including the peritrophic membrane. Thus, the chitinases disrupt the feces and facilitate release of the allergens. B. cereus was isolated and identified based on amplification and sequencing of 16S rRNA and motB genes. B. cereus was added into high-fat, high-protein (DDF) and low-fat, low-protein (flour) diets to 1 and 5% (w/w), and the diets palatability was evaluated in 21-day population growth test. The supplementation of diet with B. cereus significantly suppressed population growth and the suppressive effect was higher in the high-fat, high-protein diet than in the low-fat, low-protein food. Thus, B. cereus has to coexist with the mite in balance to be beneficial for the mite. The mite-B. cereus symbiosis can be beneficial-suppressive at some level. The results increase the veterinary and medical importance of the allergens detected in feces. The B. cereus enzymes/toxins are important components of mite allergens. The strong symbiotic association of T. putrescentiae with B. cereus in DDF was indicated.

The Phylogeny of Bacillus cereus sensu lato

The three main species of the Bacillus cereus sensu lato , B. cereus , B. thuringiensis , and B. anthracis , were recognized and established by the early 1900s because they each exhibited distinct phenotypic traits. B. thuringiensis isolates and their parasporal crystal proteins have long been established as a natural pesticide and insect pathogen. B. anthracis , the etiological agent for anthrax, was used by Robert Koch in the 19th century as a model to develop the germ theory of disease, and B. cereus , a common soil organism, is also an occasional opportunistic pathogen of humans. In addition to these three historical species designations, are three less-recognized and -understood species: B. mycoides , B. weihenstephanensis , and B. pseudomycoides . All of these “species” combined comprise the Bacillus cereus sensu lato group. Despite these apparently clear phenotypic definitions, early molecular approaches to separate the first three by various DNA hybridization and 16S/23S ribosomal sequence analyses led to some “confusion” because there were limited differences to differentiate between these species. These and other results have led to frequent suggestions that a taxonomic change was warranted to reclassify this group to a single species. But the pathogenic properties of B. anthracis and the biopesticide applications of B. thuringiensis appear to “have outweighed pure taxonomic considerations” and the separate species categories are still being maintained. B. cereus sensu lato represents a classic example of a now common bacterial species taxonomic quandary.

Infection of Tribolium castaneum with Bacillus thuringiensis: quantification of bacterial replication within cadavers, transmission via cannibalism, and inhibition of spore germination

Reproduction within a host and transmission to the next host are crucial for the virulence and fitness of pathogens. Nevertheless, basic knowledge about such parameters is often missing from the literature, even for well-studied bacteria, such as Bacillus thuringiensis, an endospore-forming insect pathogen, which infects its hosts via the oral route. To characterize bacterial replication success, we made use of an experimental oral infection system for the red flour beetle Tribolium castaneum and developed a flow cytometric assay for the quantification of both spore ingestion by the individual beetle larvae and the resulting spore load after bacterial replication and resporulation within cadavers. On average, spore numbers increased 460-fold, showing that Bacillus thuringiensis grows and replicates successfully in insect cadavers. By inoculating cadaver-derived spores and spores from bacterial stock cultures into nutrient medium, we next investigated outgrowth characteristics of vegetative cells and found that cadaver-derived bacteria showed reduced growth compared to bacteria from the stock cultures. Interestingly, this reduced growth was a consequence of inhibited spore germination, probably originating from the host and resulting in reduced host mortality in subsequent infections by cadaver-derived spores. Nevertheless, we further showed that Bacillus thuringiensis transmission was possible via larval cannibalism when no other food was offered. These results contribute to our understanding of the ecology of Bacillus thuringiensis as an insect pathogen.

Food-bacteria interplay: pathometabolism of emetic Bacillus cereus

Bacillus cereus is a Gram-positive endospore forming bacterium known for its wide spectrum of phenotypic traits, enabling it to occupy diverse ecological niches. Although the population structure of B. cereus is highly dynamic and rather panmictic, production of the emetic B. cereus toxin cereulide is restricted to strains with specific genotypic traits, associated with distinct environmental habitats. Cereulide is an ionophoric dodecadepsipeptide that is produced non-ribosomally by an enzyme complex with an unusual modular structure, named cereulide synthetase (Ces non-ribosomal peptide synthetase). The ces gene locus is encoded on a mega virulence plasmid related to the B. anthracis toxin plasmid pXO1. Cereulide, a highly thermo- and pH- resistant molecule, is preformed in food, evokes vomiting a few hours after ingestion, and was shown to be the direct cause of gastroenteritis symptoms; occasionally it is implicated in severe clinical manifestations including acute liver failures. Control of toxin gene expression in emetic B. cereus involves central transcriptional regulators, such as CodY and AbrB, thereby inextricably linking toxin gene expression to life cycle phases and specific conditions, such as the nutrient supply encountered in food matrices. While in recent years considerable progress has been made in the molecular and biochemical characterization of cereulide toxin synthesis, far less is known about the embedment of toxin synthesis in the life cycle of B. cereus. Information about signals acting on toxin production in the food environment is lacking. We summarize the data available on the complex regulatory network controlling cereulide toxin synthesis, discuss the role of intrinsic and extrinsic factors acting on toxin biosynthesis in emetic B. cereus and stress how unraveling these processes can lead to the development of novel effective strategies to prevent toxin synthesis in the food production and processing chain.

Progressive genomic convergence of two Helicobacter pylori strains during mixed infection of a patient with chronic gastritis

OBJECTIVE

To study the detailed nature of genomic microevolution during mixed infection with multiple Helicobacter pylori strains in an individual.

DESIGN

We sampled 18 isolates from a single biopsy from a patient with chronic gastritis and nephritis. Whole-genome sequencing was applied to these isolates, and statistical genetic tools were used to investigate their evolutionary history.

RESULTS

The genomes fall into two clades, reflecting colonisation of the stomach by two distinct strains, and these lineages have accumulated diversity during an estimated 2.8 and 4.2 years of evolution. We detected about 150 clear recombination events between the two clades. Recombination between the lineages is a continuous ongoing process and was detected on both clades, but the effect of recombination in one clade was nearly an order of magnitude higher than in the other. Imputed ancestral sequences also showed evidence of recombination between the two strains prior to their diversification, and we estimate that they have both been infecting the same host for at least 12 years. Recombination tracts between the lineages were, on average, 895 bp in length, and showed evidence for the interspersion of recipient sequences that has been observed in in vitro experiments. The complex evolutionary history of a phage-related protein provided evidence for frequent reinfection of both clades by a single phage lineage during the past 4 years.

CONCLUSIONS

Whole genome sequencing can be used to make detailed conclusions about the mechanisms of genetic change of H. pylori based on sampling bacteria from a single gastric biopsy.

Burkholderia pseudomallei sequencing identifies genomic clades with distinct recombination, accessory, and epigenetic profiles

Burkholderia pseudomallei (Bp) is the causative agent of the infectious disease melioidosis. To investigate population diversity, recombination, and horizontal gene transfer in closely related Bp isolates, we performed whole-genome sequencing (WGS) on 106 clinical, animal, and environmental strains from a restricted Asian locale. Whole-genome phylogenies resolved multiple genomic clades of Bp, largely congruent with multilocus sequence typing (MLST). We discovered widespread recombination in the Bp core genome, involving hundreds of regions associated with multiple haplotypes. Highly recombinant regions exhibited functional enrichments that may contribute to virulence. We observed clade-specific patterns of recombination and accessory gene exchange, and provide evidence that this is likely due to ongoing recombination between clade members. Reciprocally, interclade exchanges were rarely observed, suggesting mechanisms restricting gene flow between clades. Interrogation of accessory elements revealed that each clade harbored a distinct complement of restriction-modification (RM) systems, predicted to cause clade-specific patterns of DNA methylation. Using methylome sequencing, we confirmed that representative strains from separate clades indeed exhibit distinct methylation profiles. Finally, using an E. coli system, we demonstrate that Bp RM systems can inhibit uptake of non-self DNA. Our data suggest that RM systems borne on mobile elements, besides preventing foreign DNA invasion, may also contribute to limiting exchanges of genetic material between individuals of the same species. Genomic clades may thus represent functional units of genetic isolation in Bp, modulating intraspecies genetic diversity.

Structural insights into Bacillus thuringiensis Cry, Cyt and parasporin toxins

Since the first X-ray structure of Cry3Aa was revealed in 1991, numerous structures of B. thuringiensis toxins have been determined and published. In recent years, functional studies on the mode of action and resistance mechanism have been proposed, which notably promoted the developments of biological insecticides and insect-resistant transgenic crops. With the exploration of known pore-forming toxins (PFTs) structures, similarities between PFTs and B. thuringiensis toxins have provided great insights into receptor binding interactions and conformational changes from water-soluble to membrane pore-forming state of B. thuringiensis toxins. This review mainly focuses on the latest discoveries of the toxin working mechanism, with the emphasis on structural related progress. Based on the structural features, B. thuringiensis Cry, Cyt and parasporin toxins could be divided into three categories: three-domain type α-PFTs, Cyt toxin type β-PFTs and aerolysin type β-PFTs. Structures from each group are elucidated and discussed in relation to the latest data, respectively.

Crystalliferous Bacillus cereus group bacteria from a Maryland hardwood forest are dominated by psychrotolerant strains

Crystal-forming bacteria of the Bacillus cereus group were isolated from soil samples collected at different elevations within a mixed hardwood forest in central Maryland, and their phylogenetic relationships determined by multilocus sequence analysis. The vast majority of isolates obtained were associated with two phylogenetic groups known to be psychrotolerant, with very few isolates representing phylogenetic groups more typically associated with Bacillus thuringiensis. Isolates from the psychrotolerant groups were found to grow on solid media at 7°C. Isolates of 11 highly related, novel sequence types (STs) from the psychrotolerant group that includes Bacillus weihenstephanensis were generally found at higher elevations, and were not associated with soils near streams. Isolates of two related STs from the second psychrotolerant group were nearly always found at the bottoms of ravines near streams, in areas abundant in earthworm castings.