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

Unveiling intraspecific diversity and evolutionary dynamics of the foodborne pathogen Bacillus paranthracis through high-quality pan-genome analysis

Understanding the evolutionary dynamics of foodborne pathogens throughout host-associated habitats is of utmost importance. Bacterial pan-genomes, as dynamic entities, are strongly influenced by ecological lifestyles. As a phenotypically diverse species in the Bacillus cereus group, Bacillus paranthracis is recognized as an emerging foodborne pathogen and a probiotic simultaneously. This poorly understood species is a suitable study model for adaptive pan-genome evolution. In this study, we determined the biogeographic distribution, abundance, genetic diversity, and genotypic profiles of key genetic elements of B. paranthracis. Metagenomic read recruitment analyses demonstrated that B. paranthracis members are globally distributed and abundant in host-associated habitats. A high-quality pan-genome of B. paranthracis was subsequently constructed to analyze the evolutionary dynamics involved in ecological adaptation comprehensively. The open pan-genome indicated a flexible gene repertoire with extensive genetic diversity. Significant divergences in the phylogenetic relationships, functional enrichment, and degree of selective pressure between the different components demonstrated different evolutionary dynamics between the core and accessory genomes driven by ecological forces. Purifying selection and gene loss are the main signatures of evolutionary dynamics in B. paranthracis pan-genome. The plasticity of the accessory genome is characterized by horizontal gene transfer (HGT), massive gene losses, and weak purifying or positive selection, which might contribute to niche-specific adaptation. In contrast, although the core genome dominantly undergoes purifying selection, its association with HGT and positively selected mutations indicates its potential role in ecological diversification. Furthermore, host fitness-related dynamics are characterized by the loss of secondary metabolite biosynthesis gene clusters (BGCs) and CAZyme-encoding genes and the acquisition of antimicrobial resistance (AMR) and virulence genes via HGT. This study offers a case study of pan-genome evolution to investigate the ecological adaptations reflected by biogeographical characteristics, thereby advancing the understanding of intraspecific diversity and evolutionary dynamics of foodborne pathogens.

Development of an Online Genome Sequence Comparison Resource for Bacillus cereus sensu lato Strains Using the Efficient Composition Vector Method

An automated method was developed for differentiating closely related B. cereus sensu lato (s.l.) species, especially biopesticide Bacillus thuringiensis, from other human pathogens, B. anthracis and B. cereus sensu stricto (s.s.). In the current research, four typing methods were initially compared, including multi-locus sequence typing (MLST), single-copy core genes phylogenetic analysis (SCCGPA), dispensable genes content pattern analysis (DGCPA) and composition vector tree (CVTree), to analyze the genomic variability of 23 B. thuringiensis strains from aizawai, kurstaki, israelensis, thuringiensis and morrisoni serovars. The CVTree method was the best option to be used for typing B. thuringiensis strains since it proved to be the fastest method, whilst giving high-resolution data about the strains. In addition, CVTree agrees well with ANI-based method, revealing the relationship between B. thuringiensis and other B. cereus s.l. species. Based on these data, an online genome sequence comparison resource was built for Bacillus strains called the Bacillus Typing Bioinformatics Database to facilitate strain identification and characterization.

Pathogenicity and Genomic Characterization of a Novel Genospecies, Bacillus shihchuchen, of the Bacillus cereus Group Isolated from Chinese Softshell Turtle (Pelodiscus sinensis)

The Chinese softshell turtle (CST; Pelodiscus sinensis) is a freshwater aquaculture species of substantial economic importance that is commercially farmed across Asia, particularly in Taiwan. Although diseases caused by the Bacillus cereus group (Bcg) pose a major threat to commercial CST farming systems, information regarding its pathogenicity and genome remains limited. Here, we investigated the pathogenicity of Bcg strains isolated in a previous study and performed whole-genome sequencing. Pathogenicity analysis indicated that QF108-045 isolated from CSTs caused the highest mortality rate, and whole-genome sequencing revealed that it was an independent group distinct from other known Bcg genospecies. The average nucleotide identity compared to other known Bcg genospecies was below 95%, suggesting that QF108-045 belongs to a new genospecies, which we named Bacillus shihchuchen. Furthermore, genes annotation revealed the presence of anthrax toxins, such as edema factor and protective antigen, in QF108-045. Therefore, the biovar anthracis was assigned, and the full name of QF108-045 was Bacillus shihchuchen biovar anthracis. In addition to possessing multiple drug-resistant genes, QF108-045 demonstrated resistance to various types of antibiotics, including penicillins (amoxicillin and ampicillin), cephalosporins (ceftifour, cephalexin, and cephazolin), and polypeptides, such as vancomycin.

Environmental sporobiota: Occurrence, dissemination, and risks

Spore-forming bacteria known as sporobiota are widespread in diverse environments from terrestrial and aquatic habitats to industrial and healthcare systems. Studies on sporobiota have been mainly focused on food processing and clinical fields, while a large amount of sporobiota exist in natural environments. Due to their persistence and capabilities of transmitting virulence factors and antibiotic resistant genes, environmental sporobiota could pose significant health risks to humans. These risks could increase as global warming and environmental pollution has altered the life cycle of sporobiota. This review summarizes the current knowledge of environmental sporobiota, including their occurrence, characteristics, and functions. An interaction network among clinical-, food-related, and environment-related sporobiota is constructed. Recent and effective methods for detecting and disinfecting environmental sporobiota are also discussed. Key problems and future research needs for better understanding and reducing the risks of environmental sporobiota and sporobiome are proposed.

Selecting for infectivity across metapopulations can increase virulence in the social microbe Bacillus thuringiensis

Passage experiments that sequentially infect hosts with parasites have long been used to manipulate virulence. However, for many invertebrate pathogens, passage has been applied naively without a full theoretical understanding of how best to select for increased virulence and this has led to very mixed results. Understanding the evolution of virulence is complex because selection on parasites occurs across multiple spatial scales with potentially different conflicts operating on parasites with different life histories. For example, in social microbes, strong selection on replication rate within hosts can lead to cheating and loss of virulence, because investment in public goods virulence reduces replication rate. In this study, we tested how varying mutation supply and selection for infectivity or pathogen yield (population size in hosts) affected the evolution of virulence against resistant hosts in the specialist insect pathogen Bacillus thuringiensis, aiming to optimize methods for strain improvement against a difficult to kill insect target. We show that selection for infectivity using competition between subpopulations in a metapopulation prevents social cheating, acts to retain key virulence plasmids, and facilitates increased virulence. Increased virulence was associated with reduced efficiency of sporulation, and possible loss of function in putative regulatory genes but not with altered expression of the primary virulence factors. Selection in a metapopulation provides a broadly applicable tool for improving the efficacy of biocontrol agents. Moreover, a structured host population can facilitate artificial selection on infectivity, while selection on life-history traits such as faster replication or larger population sizes can reduce virulence in social microbes.

Identification of Genetic Markers for the Detection of Bacillus thuringiensis Strains of Interest for Food Safety

Bacillus thuringiensis (Bt), belonging to the Bacillus cereus (Bc) group, is commonly used as a biopesticide worldwide due to its ability to produce insecticidal crystals during sporulation. The use of Bt, especially subspecies aizawai and kurstaki, to control pests such as Lepidoptera, generally involves spraying mixtures containing spores and crystals on crops intended for human consumption. Recent studies have suggested that the consumption of commercial Bt strains may be responsible for foodborne outbreaks (FBOs). However, its genetic proximity to Bc strains has hindered the development of routine tests to discriminate Bt from other Bc, especially Bacillus cereus sensu stricto (Bc ss), well known for its involvement in FBOs. Here, to develop tools for the detection and the discrimination of Bt in food, we carried out a genome-wide association study (GWAS) on 286 complete genomes of Bc group strains to identify and validate in silico new molecular markers specific to different Bt subtypes. The analyses led to the determination and the in silico validation of 128 molecular markers specific to Bt, its subspecies aizawai, kurstaki and four previously described proximity clusters associated with these subspecies. We developed a command line tool based on a 14-marker workflow, to carry out a computational search for Bt-related markers from a putative Bc genome, thereby facilitating the detection of Bt of interest for food safety, especially in the context of FBOs.

Conjugation across Bacillus cereus and kin: A review

Horizontal gene transfer (HGT) is a major driving force in shaping bacterial communities. Key elements responsible for HGT are conjugation-like events and transmissible plasmids. Conjugative plasmids can promote their own transfer as well as that of co-resident plasmids. Bacillus cereus and relatives harbor a plethora of plasmids, including conjugative plasmids, which are at the heart of the group species differentiation and specification. Since the first report of a conjugation-like event between strains of B. cereus sensu lato (s.l.) 40 years ago, many have studied the potential of plasmid transfer across the group, especially for plasmids encoding major toxins. Over the years, more than 20 plasmids from B. cereus isolates have been reported as conjugative. However, with the increasing number of genomic data available, in silico analyses indicate that more plasmids from B. cereus s.l. genomes present self-transfer potential. B. cereus s.l. bacteria occupy diverse environmental niches, which were mimicked in laboratory conditions to study conjugation-related mechanisms. Laboratory mating conditions remain nonetheless simplistic compared to the complex interactions occurring in natural environments. Given the health, economic and ecological importance of strains of B. cereus s.l., it is of prime importance to consider the impact of conjugation within this bacterial group.

Simultaneous Discrimination of Cereulide-Producing Bacillus cereus and Psychrotolerant B. cereus Group by Matrix-Assisted Laser Desorption Ionization-Time-of-Flight Mass Spectrometry

ABSTRACT

Cereulide-producing Bacillus cereus, which causes foodborne illnesses with vomiting, and psychrotolerant B. cereus group strains such as Bacillus mycoides, which can grow at ≥7°C and cause spoilage of refrigerated foods, are significant concerns for the food industry. Rapid and simple methods to discriminate the cereulide-producing B. cereus and psychrotolerant B. cereus group strains from other B. cereus group strains are needed. We developed a novel, rapid, and simple method with matrix-assisted laser desorption ionization-time-of-flight mass spectrometry (MALDI-TOF MS) analysis for simultaneous discrimination of these two groups from other B. cereus group strains. A potassium adduct of cereulide was used to detect cereulide-producing B. cereus, and three ribosomal subunit proteins (L30, S16, and S20) were used to detect psychrotolerant B. cereus group. A total of 51 B. cereus group strains were analyzed by MALDI-TOF MS. The biomarkers allowed successful discrimination of 16 cereulide-producing B. cereus and 15 psychrotolerant B. cereus group strains from other B. cereus group strains. The results showed that this MALDI-TOF MS analysis allows simultaneous discrimination of cereulide-producing B. cereus and psychrotolerant B. cereus group strains from other B. cereus group strains. This efficient method has the potential to be a valuable tool for ensuring food safety.

HIGHLIGHTS

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.

Differentiation of Bacillus cereus and Bacillus thuringiensis Using Genome-Guided MALDI-TOF MS Based on Variations in Ribosomal Proteins

Bacillus cereus and B. thuringiensis are closely related species that are relevant to foodborne diseases and biopesticides, respectively. Unambiguous differentiation of these two species is crucial for bacterial taxonomy. As genome analysis offers an objective but time-consuming classification of B. cereus and B. thuringiensis, in the present study, matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) was used to accelerate this process. By combining in silico genome analysis and MALDI-TOF MS measurements, four species-specific peaks of B. cereus and B. thuringiensis were screened and identified. The species-specific peaks of B. cereus were m/z 3211, 6427, 9188, and 9214, and the species-specific peaks of B. thuringiensis were m/z 3218, 6441, 9160, and 9229. All the above peaks represent ribosomal proteins, which are conserved and consistent with the phylogenetic relationship between B. cereus and B. thuringiensis. The specificity of the peaks was robustly verified using common foodborne pathogens. Thus, we concluded that genome-guided MALDI-TOF MS allows high-throughput differentiation of B. cereus and B. thuringiensis and provides a framework for differentiating other closely related species.

Bacterial cooperation through horizontal gene transfer

Cooperation exists across all scales of biological organization, from genetic elements to complex human societies. Bacteria cooperate by secreting molecules that benefit all individuals in the population (i.e., public goods). Genes associated with cooperation can spread among strains through horizontal gene transfer (HGT). We discuss recent findings on how HGT mediated by mobile genetic elements promotes bacterial cooperation, how cooperation in turn can facilitate more frequent HGT, and how the act of HGT itself may be considered as a form of cooperation. We propose that HGT is an important enforcement mechanism in bacterial populations, thus creating a positive feedback loop that further maintains cooperation. To enforce cooperation, HGT serves as a homogenizing force by transferring the cooperative trait, effectively eliminating cheaters.

Antimicrobial resistance in Bacillus-based biopesticide products

The crisis of antimicrobial resistant bacterial infections is one of the most pressing public health issues. Common agricultural practices have been implicated in the generation of antimicrobial resistant bacteria. Biopesticides, live bacteria used for pest control, are non-pathogenic and considered safe for consumption. Application of bacteria-based pesticides to crops in high concentrations raises the possibility of unintentional contributions to the movement and generation of antimicrobial resistance genes in the environment. However, the presence of clinically relevant antimicrobial resistance genes and their resistance phenotypes are currently unknown. Here we use a combination of multiple bioinformatic and microbiological techniques to define resistomes of widely used biopesticides and determine how the presence of suspected antimicrobial resistance genes translates to observable resistance phenotypes in several biopesticide products. Our results demonstrate that biopesticide products are reservoirs of clinically relevant antimicrobial resistance genes and bear resistance to multiple drug classes.

ORFograph: search for novel insecticidal protein genes in genomic and metagenomic assembly graphs

BACKGROUND

Since the prolonged use of insecticidal proteins has led to toxin resistance, it is important to search for novel insecticidal protein genes (IPGs) that are effective in controlling resistant insect populations. IPGs are usually encoded in the genomes of entomopathogenic bacteria, especially in large plasmids in strains of the ubiquitous soil bacteria, Bacillus thuringiensis (Bt). Since there are often multiple similar IPGs encoded by such plasmids, their assemblies are typically fragmented and many IPGs are scattered through multiple contigs. As a result, existing gene prediction tools (that analyze individual contigs) typically predict partial rather than complete IPGs, making it difficult to conduct downstream IPG engineering efforts in agricultural genomics.

METHODS

Although it is difficult to assemble IPGs in a single contig, the structure of the genome assembly graph often provides clues on how to combine multiple contigs into segments encoding a single IPG.

RESULTS

We describe ORFograph, a pipeline for predicting IPGs in assembly graphs, benchmark it on (meta)genomic datasets, and discover nearly a hundred novel IPGs. This work shows that graph-aware gene prediction tools enable the discovery of greater diversity of IPGs from (meta)genomes.

CONCLUSIONS

We demonstrated that analysis of the assembly graphs reveals novel candidate IPGs. ORFograph identified both already known genes "hidden" in assembly graphs and potential novel IPGs that evaded existing tools for IPG identification. As ORFograph is fast, one could imagine a pipeline that processes many (meta)genomic assembly graphs to identify even more novel IPGs for phenotypic testing than would previously be inaccessible by traditional gene-finding methods. While here we demonstrated the results of ORFograph only for IPGs, the proposed approach can be generalized to any class of genes. Video abstract.

Comparative Genomics of Prophages Sato and Sole Expands the Genetic Diversity Found in the Genus Betatectivirus

Tectiviruses infecting the Bacillus cereus group represent part of the bacterial "plasmid repertoire" as they behave as linear plasmids during their lysogenic cycle. Several novel tectiviruses have been recently found infecting diverse strains belonging the B. cereus lineage. Here, we report and analyze the complete genome sequences of phages Sato and Sole. The linear dsDNA genome of Sato spans 14,852 bp with 32 coding DNA sequences (CDSs), whereas the one of Sole has 14,444 bp comprising 30 CDSs. Both phage genomes contain inverted terminal repeats and no tRNAs. Genomic comparisons and phylogenetic analyses placed these two phages within the genus Betatectivirus in the family Tectiviridae. Additional comparative genomic analyses indicated that the "gene regulation-genome replication" module of phages Sato and Sole is more diverse than previously observed among other fully sequenced betatectiviruses, displaying very low sequence similarities and containing some ORFans. Interestingly, the ssDNA binding protein encoded in this genomic module in phages Sato and Sole has very little amino acid similarity with those of reference betatectiviruses. Phylogenetic analyses showed that both Sato and Sole represent novel tectivirus species, thus we propose to include them as two novel species in the genus Betatectivirus.

Evolution of virulence in a novel family of transmissible mega-plasmids

Some Serratia entomophila isolates have been successfully exploited in biopesticides due to their ability to cause amber disease in larvae of the Aotearoa (New Zealand) endemic pasture pest, Costelytra giveni. Anti-feeding prophage and ABC toxin complex virulence determinants are encoded by a 153-kb single-copy conjugative plasmid (pADAP; amber disease-associated plasmid). Despite growing understanding of the S. entomophila pADAP model plasmid, little is known about the wider plasmid family. Here, we sequence and analyse mega-plasmids from 50 Serratia isolates that induce variable disease phenotypes in the C. giveni insect host. Mega-plasmids are highly conserved within S. entomophila, but show considerable divergence in Serratia proteamaculans with other variants in S. liquefaciens and S. marcescens, likely reflecting niche adaption. In this study to reconstruct ancestral relationships for a complex mega-plasmid system, strong co-evolution between Serratia species and their plasmids were found. We identify 12 distinct mega-plasmid genotypes, all sharing a conserved gene backbone, but encoding highly variable accessory regions including virulence factors, secondary metabolite biosynthesis, Nitrogen fixation genes and toxin-antitoxin systems. We show that the variable pathogenicity of Serratia isolates is largely caused by presence/absence of virulence clusters on the mega-plasmids, but notably, is augmented by external chromosomally encoded factors.

The Distribution of Several Genomic Virulence Determinants Does Not Corroborate the Established Serotyping Classification of Bacillus thuringiensis

Bacillus thuringiensis, commonly referred to as Bt, is an object of the lasting interest of microbiologists due to its highly effective insecticidal properties, which make Bt a prominent source of biologicals. To categorize the exuberance of Bt strains discovered, serotyping assays are utilized in which flagellin serves as a primary seroreactive molecule. Despite its convenience, this approach is not indicative of Bt strains' phenotypes, neither it reflects actual phylogenetic relationships within the species. In this respect, comparative genomic and proteomic techniques appear more informative, but their use in Bt strain classification remains limited. In the present work, we used a bottom-up proteomic approach based on fluorescent two-dimensional difference gel electrophoresis (2D-DIGE) coupled with liquid chromatography/tandem mass spectrometry(LC-MS/MS) protein identification to assess which stage of Bt culture, vegetative or spore, would be more informative for strain characterization. To this end, the proteomic differences for the israelensis-attributed strains were assessed to compare sporulating cultures of the virulent derivative to the avirulent one as well as to the vegetative stage virulent bacteria. Using the same approach, virulent spores of the israelensis strain were also compared to the spores of strains belonging to two other major Bt serovars, namely darmstadiensis and thuringiensis. The identified proteins were analyzed regarding the presence of the respective genes in the 104 Bt genome assemblies available at open access with serovar attributions specified. Of 21 proteins identified, 15 were found to be encoded in all the present assemblies at 67% identity threshold, including several virulence factors. Notable, individual phylogenies of these core genes conferred neither the serotyping nor the flagellin-based phylogeny but corroborated the reconstruction based on phylogenomics approaches in terms of tree topology similarity. In its turn, the distribution of accessory protein genes was not confined to the existing serovars. The obtained results indicate that neither gene presence nor the core gene sequence may serve as distinctive bases for the serovar attribution, undermining the notion that the serotyping system reflects strains' phenotypic or genetic similarity. We also provide a set of loci, which fit in with the phylogenomics data plausibly and thus may serve for draft phylogeny estimation of the novel strains.

Genome evolution and the emergence of pathogenicity in avian Escherichia coli

Chickens are the most common birds on Earth and colibacillosis is among the most common diseases affecting them. This major threat to animal welfare and safe sustainable food production is difficult to combat because the etiological agent, avian pathogenic Escherichia coli (APEC), emerges from ubiquitous commensal gut bacteria, with no single virulence gene present in all disease-causing isolates. Here, we address the underlying evolutionary mechanisms of extraintestinal spread and systemic infection in poultry. Combining population scale comparative genomics and pangenome-wide association studies, we compare E. coli from commensal carriage and systemic infections. We identify phylogroup-specific and species-wide genetic elements that are enriched in APEC, including pathogenicity-associated variation in 143 genes that have diverse functions, including genes involved in metabolism, lipopolysaccharide synthesis, heat shock response, antimicrobial resistance and toxicity. We find that horizontal gene transfer spreads pathogenicity elements, allowing divergent clones to cause infection. Finally, a Random Forest model prediction of disease status (carriage vs. disease) identifies pathogenic strains in the emergent ST-117 poultry-associated lineage with 73% accuracy, demonstrating the potential for early identification of emergent APEC in healthy flocks.

mMGE: a database for human metagenomic extrachromosomal mobile genetic elements

Extrachromosomal mobile genetic elements (eMGEs), including phages and plasmids, that can move across different microbes, play important roles in genome evolution and shaping the structure of microbial communities. However, we still know very little about eMGEs, especially their abundances, distributions and putative functions in microbiomes. Thus, a comprehensive description of eMGEs is of great utility. Here we present mMGE, a comprehensive catalog of 517 251 non-redundant eMGEs, including 92 492 plasmids and 424 759 phages, derived from diverse body sites of 66 425 human metagenomic samples. About half the eMGEs could be further grouped into 70 074 clusters using relaxed criteria (referred as to eMGE clusters below). We provide extensive annotations of the identified eMGEs including sequence characteristics, taxonomy affiliation, gene contents and their prokaryotic hosts. We also calculate the prevalence, both within and across samples for each eMGE and eMGE cluster, enabling users to see putative associations of eMGEs with human phenotypes or their distribution preferences. All eMGE records can be browsed or queried in multiple ways, such as eMGE clusters, metagenomic samples and associated hosts. The mMGE is equipped with a user-friendly interface and a BLAST server, facilitating easy access/queries to all its contents easily. mMGE is freely available for academic use at: https://mgedb.comp-sysbio.org.

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.

Strong Environment-Genotype Interactions Determine the Fitness Costs of Antibiotic Resistance In Vitro and in an Insect Model of Infection

The acquisition of antibiotic resistance commonly imposes fitness costs, a reduction in the fitness of bacteria in the absence of drugs. These costs have been quantified primarily using in vitro experiments and a small number of in vivo studies in mice, and it is commonly assumed that these diverse methods are consistent. Here, we used an insect model of infection to compare the fitness costs of antibiotic resistance in vivo to those in vitro Experiments explored diverse mechanisms of resistance in a Gram-positive pathogen, Bacillus thuringiensis, and a Gram-negative intestinal symbiont, Enterobacter cloacae Rifampin resistance in B. thuringiensis showed fitness costs that were typically elevated in vivo, although these were modulated by genotype-environment interactions. In contrast, resistance to cefotaxime via derepression of AmpC β-lactamase in E. cloacae resulted in no detectable costs in vivo or in vitro, while spontaneous resistance to nalidixic acid, and carriage of the IncP plasmid RP4, imposed costs that increased in vivo Overall, fitness costs in vitro were a poor predictor of fitness costs in vivo because of strong genotype-environment interactions throughout this study. Insect infections provide a cheap and accessible means of assessing the fitness consequences of resistance mutations, data that are important for understanding the evolution and spread of resistance. This study emphasizes that the fitness costs imposed by particular mutations or different modes of resistance are extremely variable and that only a subset of these mutations is likely to be prevalent outside the laboratory.

Completed Genomic Sequence of Bacillus thuringiensis HER1410 Reveals a Cry-Containing Chromosome, Two Megaplasmids, and an Integrative Plasmidial Prophage

Bacillus thuringiensis is the most used biopesticide in agriculture. Its entomopathogenic capacity stems from the possession of plasmid-borne insecticidal crystal genes (cry), traditionally used as discriminant taxonomic feature for that species. As such, crystal and plasmid identification are key to the characterization of this species. To date, about 600 B. thuringiensis genomes have been reported, but less than 5% have been completed, while the other draft genomes are incomplete, hindering full plasmid delineation. Here we present the complete genome of Bacillus thuringiensis HER1410, a strain closely related to B. thuringiensis entomocidus and a known host for a variety of Bacillus phages. The combination of short and long-read techniques allowed fully resolving the genome and delineation of three plasmids. This enabled the accurate detection of an unusual location of a unique cry gene, cry1Ba4, located in a genomic island near the chromosome replication origin. Two megaplasmids, pLUSID1 and pLUSID2 could be delineated: pLUSID1 (368 kb), a likely conjugative plasmid involved in virulence, and pLUSID2 (156 kb) potentially related to the sporulation process. A smaller plasmidial prophage pLUSID3, with a dual lifestyle whose integration within the chromosome causes the disruption of a flagellar key component. Finally, phylogenetic analysis placed this strain within a clade comprising members from the B. thuringiensis serovar thuringiensis and other serovars and with B. cereus s. s. in agreement with the intermingled taxonomy of B. cereus sensu lato group.

Discrimination of psychrotolerant Bacillus cereus group based on MALDI-TOF MS analysis of ribosomal subunit proteins

Psychrotolerant species of the Bacillus cereus group, Bacillus mycoides and Bacillus weihenstephanensis, can grow at ≥ 7 °C and are significant concerns for the food industry due to their ability to cause spoilage of refrigerated food. In addition to that, some strains of B. weihenstephanensis can produce emetic toxin, namely cereulide, which is known to cause vomiting. Therefore, rapid and simple methods to discriminate psychrotolerant B. cereus group species are crucial. Here, matrix assisted laser desorption ionization time of flight mass spectrometry (MALDI-TOF MS) and the S10-spc-alpha operon gene encoded ribosomal protein mass spectrum (S10-GERMS) method were used to discriminate psychrotolerant species of the B. cereus group based on a set of four ribosomal subunit proteins (S10, S16, S20 and L30). A total of 36 strains of B. cereus group were cultured on LB agar, and analyzed by MALDI-TOF MS. The four biomarkers successfully discriminated 12 strains of psychrotolerant species from mesophilic species of the B. cereus group. Furthermore, the four biomarkers also classified some Bacillus thuringiensis strains. MALDI-TOF MS analysis using the S10-GERMS method allowed simple and rapid discrimination of psychrotolerant species of the B. cereus group from other mesophilic species. This method has a possibility to enable manufacturers and distributors of refrigerated foods to control psychrotolerant species of the B. cereus group effectively.

Temperature-Dependent Growth Characteristics of Bacillus thuringiensis in a Ratatouille Food Model

ABSTRACT

In contrast to Bacillus cereus, the role of Bacillus thuringiensis in foodborne illness has been controversially discussed. As B. thuringiensis-based biopesticides containing a mixture of crystal toxins and viable spores are widely used, a current European Food Safety Authority opinion underlines the need for additional data to enable risk assessment. However, it is currently poorly understood if B. thuringiensis is able to multiply in food, which is crucial to sound risk assessment. Therefore, the aim of this study was to investigate growth of selected B. thuringiensis strains from food and insecticides in a ratatouille food model. To this end, the growth parameters of three B. thuringiensis strains were determined: insecticide strain ABTS-351 (CH_119, B. thuringiensis serovar kurstaki), insecticide strain ABTS-1857 (CH_121, B. thuringiensis serovar aizawai), and CH_48 (wild-type B. thuringiensis isolated from rosemary), producing extremely high levels of enterotoxins. After an initial drop in colony counts, we observed a statistically significant growth for the tested B. thuringiensis strains between 6 and 24 h at 22, 30, and 37°C, conditions mimicking prolonged holding times. We were also able to show that the enterotoxin overproducer CH_48 can grow up to 108 CFU/g in the ratatouille matrix within 24 h at 37°C. The two midlevel enterotoxin formers ABTS-351 (CH_119) and ABTS-1857 (CH_121) isolated from biopesticides exhibited growth between 6 and 24 h, with one of the strains growing to 107 CFU/g. To our knowledge, this is the first study providing evidence of B. thuringiensis growth in a food model with intact competitive flora. Our findings suggest strain-specific variation and stress the complexity of assessing the risk related to B. thuringiensis in food, indicating that some strains can represent a risk to consumer health when vegetable-based foods are stored under conditions of prolonged temperature abuse.

HIGHLIGHTS

Proposal of a Taxonomic Nomenclature for the Bacillus cereus Group Which Reconciles Genomic Definitions of Bacterial Species with Clinical and Industrial Phenotypes

The Bacillus cereus group comprises numerous closely related species, including bioterrorism agent B. anthracis, foodborne pathogen B. cereus, and biopesticide B. thuringiensis Differentiating organisms capable of causing illness or death from those used in industry is essential for risk assessment and outbreak preparedness. However, current species definitions facilitate species-phenotype incongruences, particularly when horizontally acquired genes are responsible for a phenotype. Using all publicly available B. cereus group genomes (n = 2,231), we show that current species definitions lead to overlapping genomospecies clusters, in which 66.2% of genomes belong to multiple genomospecies at a conventional 95 average nucleotide identity (ANI) genomospecies threshold. A genomospecies threshold of ≈92.5 ANI is shown to reflect a natural gap in genome similarity for the B. cereus group, and medoid genomes identified at this threshold are shown to yield resolvable genomospecies clusters with minimal overlap (six of 2,231 genomes assigned to multiple genomospecies; 0.269%). We thus propose a nomenclatural framework for the B. cereus group which accounts for (i) genomospecies using resolvable genomospecies clusters obtained at ≈92.5 ANI, (ii) established lineages of medical importance using a formal collection of subspecies names, and (iii) heterogeneity of clinically and industrially important phenotypes using a formalized and extended collection of biovar terms. We anticipate that the proposed nomenclature will remain interpretable to clinicians, without sacrificing genomic species definitions, which can in turn aid in pathogen surveillance; early detection of emerging, high-risk genotypes; and outbreak preparedness.IMPORTANCE Historical species definitions for many prokaryotes, including pathogens, have relied on phenotypic characteristics that are inconsistent with genome evolution. This scenario forces microbiologists and clinicians to face a tradeoff between taxonomic rigor and clinical interpretability. Using the Bacillus cereus group as a model, a conceptual framework for the taxonomic delineation of prokaryotes which reconciles genomic definitions of species with clinically and industrially relevant phenotypes is presented. The nomenclatural framework outlined here serves as a model for genomics-based bacterial taxonomy that moves beyond arbitrarily set genomospecies thresholds while maintaining congruence with phenotypes and historically important species names.

Bacteria from natural populations transfer plasmids mostly towards their kin

Plasmids play a key role in microbial ecology and evolution, yet the determinants of plasmid transfer rates are poorly understood. Particularly, interactions between donor hosts and potential recipients are understudied. Here, we investigate the importance of genetic similarity between naturally co-occurring Escherichia coli isolates in plasmid transfer. We uncover extensive variability, spanning over five orders of magnitude, in the ability of isolates to donate and receive two different plasmids, R1 and RP4. Overall, transfer is strongly biased towards clone-mates, but not correlated to genetic distance when donors and recipients are not clone-mates. Transfer is limited by the presence of a functional restriction-modification system in recipients, suggesting sharing of strain-specific defence systems contributes to bias towards kin. Such restriction of transfer to kin sets the stage for longer-term coevolutionary interactions leading to mutualism between plasmids and bacterial hosts in natural communities.

Diversity of the Rap-Phr quorum-sensing systems in the Bacillus cereus group

Bacteria of the Bacillus cereus group colonize several ecological niches and infect different hosts. Bacillus cereus, a ubiquitous species causing food poisoning, Bacillus thuringiensis, an entomopathogen, and Bacillus anthracis, which is highly pathogenic to mammals, are the most important species of this group. These species are closely related genetically, and their specific toxins are encoded by plasmids. The infectious cycle of B. thuringiensis in its insect host is regulated by quorum-sensing systems from the RNPP family. Among them, the Rap-Phr systems, which are well-described in Bacillus subtilis, regulate essential processes, such as sporulation. Given the importance of these systems, we performed a global in silico analysis to investigate their prevalence, distribution, diversity and their role in sporulation in B. cereus group species. The rap-phr genes were identified in all selected strains with 30% located on plasmids, predominantly in B. thuringiensis. Despite a high variability in their sequences, there is a remarkable association between closely related strains and their Rap-Phr profile. Based on the key residues involved in RapH phosphatase activity, we predicted that 32% of the Rap proteins could regulate sporulation by preventing the phosphorylation of Spo0F. These Rap are preferentially located on plasmids and mostly related to B. thuringiensis. The predictions were partially validated by in vivo sporulation experiments suggesting that the residues linked to the phosphatase function are necessary but not sufficient to predict this activity. The wide distribution and diversity of Rap-Phr systems could strictly control the commitment to sporulation and then improve the adaptation capacities of the bacteria to environmental changes.

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.

Bacillus wiedmannii biovar thuringiensis: A Specialized Mosquitocidal Pathogen with Plasmids from Diverse Origins

Bacillus cereus sensu lato also known as B. cereus group is composed of an ecologically diverse bacterial group with an increasing number of related species, some of which are medically or agriculturally important. Numerous efforts have been undertaken to allow presumptive differentiation of B. cereus group species from one another. FCC41 is a Bacillus sp. strain toxic against mosquito species like Aedes aegypti, Aedes (Ochlerotatus) albifasciatus, Culex pipiens, Culex quinquefasciatus, and Culex apicinus, some of them responsible for the transmission of vector-borne diseases. Here, we report the complete genome sequence of FCC41 strain, which consists of one circular chromosome and eight circular plasmids ranging in size from 8 to 490 kb. This strain harbors six crystal protein genes, including cry24Ca, two cry4-like and two cry52-like, a cry41-like parasporin gene and multiple virulence factors. The phylogenetic analysis of the whole-genome sequence of this strain with molecular approaches places this strain into the Bacillus wiedmannii cluster. However, according with phenotypical characteristics such as the mosquitocidal activity due to the presence of Cry proteins found in the parasporal body and cry genes encoded in plasmids of different sizes, indicate that this strain could be renamed as B. wiedmannii biovar thuringiensis strain FCC41.

Enterotoxin Production of Bacillus thuringiensis Isolates From Biopesticides, Foods, and Outbreaks

While the relevance of Bacillus (B.) cereus as a major cause of gastroenteritis is undisputed, the role of the closely related B. thuringiensis in foodborne disease is unclear. B. thuringiensis strains frequently harbor enterotoxin genes. However, the organism has only very rarely been associated with foodborne outbreaks, possibly due to the fact that during outbreak investigations, B. cereus is routinely not differentiated from B. thuringiensis. A recent EFSA scientific opinion stresses the urgent need for further data allowing for improved risk assessment, in particular as B. thuringiensis is a commonly used biopesticide. Therefore, the aim of this study was to gain further insights into the hazardous potential of B. thuringiensis. To this end, 39 B. thuringiensis isolates obtained from commercially used biopesticides, various food sources, as well as from foodborne outbreaks were characterized by panC typing, panC-based SplitsTree analysis, toxin gene profiling, FTIR spectroscopic analysis, a cytotoxicity assay screening for enterotoxic activity, and a sphingomyelinase assay. The majority of the tested B. thuringiensis isolates exhibited low (23%, n = 9) or mid level enterotoxicity (74%, n = 29), and produced either no (59%, n = 23) or low levels (33%, n = 13) of sphingomyelinase, which is reported to act synergistically with enterotoxins Nhe and Hbl. One strain isolated from rosemary was however classified as highly enterotoxic surpassing the cytotoxic activity of the high-level reference strain by a factor of 1.5. This strain also produced vast amounts of sphingomyelinase. Combining all results obtained in this study into a fingerprint pattern, several enterotoxic biopesticide strains were indistinguishable from those of isolates from foods or collected in association with outbreaks. Our study shows that many B. thuringiensis biopesticide strains exhibit mid-level cytotoxicity in a Vero cell assay and that some of these strains cannot be differentiated from isolates collected from foods or in association with outbreaks. Thus, we demonstrate that the use of B. thuringiensis strains as biopesticides can represent a food safety risk, underpinning the importance of assessing the hazardous potential of each strain and formulation used.

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.