The fungus Escovopsis (Ascomycota: Hypocreales): a critical review of its biology and parasitism of attine ant colonies

Two biological phenomena that contribute to increasing complexity in biological systems are mutualistic symbiotic interactions and the evolution of sociality. These two phenomena are also of fundamental importance to our understanding of the natural world. An organism that poses a threat to one or both of these is therefore also of great interest as it represents a challenge that mutualistic symbioses and social organisms have to overcome. This is the case with the fungus Escovopsis (Ascomycota: Hypocreales), which attacks the fungus garden of attine ants (Formicidae: Attina) such as the leaf cutters. This parasite has attracted much high-profile scientific interest for considerable time, and its study has been fruitful in understanding evolutionary, ecological and behavioural processes. Despite this, much of the biology and ecology of this organism remains unknown. Here we discuss this fungus and three sister genera (Escovopsioides, Luteomyces and Sympodiorosea) that until recently were considered as a single group. We first describe its position as the most highly specialised microbial symbiont in this system other than the mutualistic fungal cultivar itself and as that of greatest scientific interest. We then review the taxonomic history of the group and its macroevolution and biogeography. We examine what we know of its life cycle in the field - surprisingly little is known of how it is transmitted between colonies, but we explain what is known to date. We then review how it interacts with its host(s), first at the level of its direct interaction with the basidiomycete host fungi wherein we show the evidence for it being a mycoparasite; then at the colony level where empirical evidence points towards it being a parasite with a very low virulence or even merely a opportunist. Finally, we offer directions for future research.

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.

Comparative genomics of Bacillus cereus sensu lato spp. biocontrol strains in correlation to in-vitro phenotypes and plant pathogen antagonistic capacity

Bacillus cereus sensu lato (Bcsl) strains are widely explored due to their capacity to antagonize a broad range of plant pathogens. These include B. cereus sp. UW85, whose antagonistic capacity is attributed to the secondary metabolite Zwittermicin A (ZwA). We recently isolated four soil and root-associated Bcsl strains (MO2, S-10, S-25, LSTW-24) that displayed different growth profiles and in-vitro antagonistic effects against three soilborne plant pathogens models: Pythium aphanidermatum (oomycete) Rhizoctonia solani (basidiomycete), and Fusarium oxysporum (ascomycete). To identify genetic mechanisms potentially responsible for the differences in growth and antagonistic phenotypes of these Bcsl strains, we sequenced and compared their genomes, and that of strain UW85 using a hybrid sequencing pipeline. Despite similarities, specific Bcsl strains had unique secondary metabolite and chitinase-encoding genes that could potentially explain observed differences in in-vitro chitinolytic potential and anti-fungal activity. Strains UW85, S-10 and S-25 contained a (~500 Kbp) mega-plasmid that harbored the ZwA biosynthetic gene cluster. The UW85 mega-plasmid contained more ABC transporters than the other two strains, whereas the S-25 mega-plasmid carried a unique cluster containing cellulose and chitin degrading genes. Collectively, comparative genomics revealed several mechanisms that can potentially explain differences in in-vitro antagonism of Bcsl strains toward fungal plant pathogens.

A small bacteriophage protein determines the hierarchy over co-residential jumbo phage in Bacillus thuringiensis serovar israelensis

Bacillus thuringiensis serovar israelensis is the most widely used biopesticide against insects, including vectors of animal and human diseases. Among several extrachromosomal elements, this endospore-forming entomopathogen harbors two bacteriophages: a linear DNA replicon named GIL01 that does not integrate into the chromosome during lysogeny and a circular-jumbo prophage known as pBtic235. Here, we show that GIL01 hinders the induction of cohabiting prophage pBtic235. The GIL01-encoded small protein, gp7, which interacts with the host LexA repressor, is a global transcription regulator and represses the induction of pBtic235 after DNA damage to presumably allow GIL01 to multiply first. In a complex with host LexA in stressed cells, gp7 down-regulates the expression of more than 250 host and pBtic235 genes, many of which are involved in the cellular functions of genome maintenance, cell-wall transport, and membrane and protein stability. We show that gp7 homologs that are found exclusively in bacteriophages act in a similar fashion to enhance LexA's binding to DNA, while likely also affecting host gene expression. Our results provide evidence that GIL01 influences both its host and its co-resident bacteriophage.

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.

Adaptation of Bacillus thuringiensis to Plant Colonization Affects Differentiation and Toxicity

The Bacillus cereus group (Bacillus cereus sensu lato) has a diverse ecology, including various species that are vertebrate or invertebrate pathogens. Few isolates from the B. cereus group have however been demonstrated to benefit plant growth. Therefore, it is crucial to explore how bacterial development and pathogenesis evolve during plant colonization. Herein, we investigated Bacillus thuringiensis (Cry-) adaptation to the colonization of Arabidopsis thaliana roots and monitored changes in cellular differentiation in experimentally evolved isolates. Isolates from two populations displayed improved iterative ecesis on roots and increased virulence against insect larvae. Molecular dissection and recreation of a causative mutation revealed the importance of a nonsense mutation in the rho transcription terminator gene. Transcriptome analysis revealed how Rho impacts various B. thuringiensis genes involved in carbohydrate metabolism and virulence. Our work suggests that evolved multicellular aggregates have a fitness advantage over single cells when colonizing plants, creating a trade-off between swimming and multicellularity in evolved lineages, in addition to unrelated alterations in pathogenicity. IMPORTANCE Biologicals-based plant protection relies on the use of safe microbial strains. During application of biologicals to the rhizosphere, microbes adapt to the niche, including genetic mutations shaping the physiology of the cells. Here, the experimental evolution of Bacillus thuringiensis lacking the insecticide crystal toxins was examined on the plant root to reveal how adaptation shapes the differentiation of this bacterium. Interestingly, evolution of certain lineages led to increased hemolysis and insect larva pathogenesis in B. thuringiensis driven by transcriptional rewiring. Further, our detailed study reveals how inactivation of the transcription termination protein Rho promotes aggregation on the plant root in addition to altered differentiation and pathogenesis in B. thuringiensis.

The roles of environmental variation and parasite survival in virulence-transmission relationships

Disease outbreaks are a consequence of interactions among the three components of a host-parasite system: the infectious agent, the host and the environment. While virulence and transmission are widely investigated, most studies of parasite life-history trade-offs are conducted with theoretical models or tractable experimental systems where transmission is standardized and the environment controlled. Yet, biotic and abiotic environmental factors can strongly affect disease dynamics, and ultimately, host-parasite coevolution. Here, we review research on how environmental context alters virulence-transmission relationships, focusing on the off-host portion of the parasite life cycle, and how variation in parasite survival affects the evolution of virulence and transmission. We review three inter-related 'approaches' that have dominated the study of the evolution of virulence and transmission for different host-parasite systems: (i) evolutionary trade-off theory, (ii) parasite local adaptation and (iii) parasite phylodynamics. These approaches consider the role of the environment in virulence and transmission evolution from different angles, which entail different advantages and potential biases. We suggest improvements to how to investigate virulence-transmission relationships, through conceptual and methodological developments and taking environmental context into consideration. By combining developments in life-history evolution, phylogenetics, adaptive dynamics and comparative genomics, we can improve our understanding of virulence-transmission relationships across a diversity of host-parasite systems that have eluded experimental study of parasite life history.

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.

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.

Application of Bacillus thuringiensis strains with conjugal and mobilizing capability drives gene transmissibility within Bacillus cereus group populations in confined habitats

Background

Bacillus thuringiensis bacteria share similar genetic, physiological, and biochemical characteristics with other members of the Bacillus cereus group. Their diversity and entomopathogenic origin are related to their mobile genetic elements. However, the effects of wide-spread application of B. thuringiensis-based pesticides on genetically related B. cereus group populations present in the environment remain poorly understood.

Results

We first identified pBMB76 from B. thuringiensis tenebrionis as a new conjugative plasmid. Mixed mating experiments suggested that pBMB76 may compete with pHT73, another known conjugative plasmid. Applications of single (tenebrionis 4AA1 and kurstaki HD73 carrying pBMB76 and pHT73, respectively) and mixed (4AA1 + HD73) B. thuringiensis strains were performed in confined plot habitats (soil and leaf) over two planting seasons. In total, 684 B. cereus group isolates were randomly selected from different treatment sets, and the transmissibility and occurrence rate of potential conjugative plasmids were surveyed. Results showed that the percentage of isolates with plasmid mobility was markedly enhanced in the B. thuringiensis-sprayed groups. Furthermore, we performed multi-locus sequence typing (MLST) for a subset of 291 isolates, which indicated that the dominant sequence types in the treated habitats were identical or related to the corresponding sprayed formulations.

Conclusions

The application of B. thuringiensis strains with conjugal and mobilizing capability drove gene transmissibility within the B. cereus group populations in confined habitats and potentially modified the population structure.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-020-02047-4.

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.

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.

Optimal Response to Quorum-Sensing Signals Varies in Different Host Environments with Different Pathogen Group Size

Quorum sensing describes the ability of microbes to alter gene regulation according to their local population size. Some successful theory suggests that this is a form of cooperation, namely, investment in shared products is only worthwhile if there are sufficient bacteria making the same product. This theory can explain the genetic diversity in these signaling systems in Gram-positive bacteria, such as Bacillus and Staphylococcus sp. The possible advantages gained by rare genotypes (which can exploit the products of their more common neighbors) could explain why different genotypes can coexist. We show that while these social interactions can occur in simple laboratory experiments, they do not occur in naturalistic infections using an invertebrate pathogen, Bacillus thuringiensis. Instead, our results suggest that different genotypes are adapted to differently sized hosts. Overall, social models are not easily applied to this system, implying that a different explanation for this form of quorum sensing is required.

The persistence of genetic variation in master regulators of gene expression, such as quorum-sensing systems, is hard to explain. Here, we investigated two alternative hypotheses for the prevalence of polymorphic quorum sensing in Gram-positive bacteria, i.e., the use of different signal/receptor pairs (‘pherotypes’) to regulate the same functions. First, social interactions between pherotypes or ‘facultative cheating’ may favor rare variants that exploit the signals of others. Second, different pherotypes may increase fitness in different environments. We evaluated these hypotheses in the invertebrate pathogen Bacillus thuringiensis, using three pherotypes expressed in a common genetic background. Facultative cheating could occur in well-mixed host homogenates provided there was minimal cross talk between competing pherotypes. However, facultative cheating did not occur when spatial structure was increased in static cultures or in naturalistic oral infections, where common pherotypes had higher fitness. There was clear support for environment-dependent fitness; pherotypes varied in responsiveness to signals and in mean competitive fitness. Notably, competitive fitness varied with group size. In contrast to typical social evolution models of quorum sensing which predict higher response to signal at larger group size, the pherotype with highest responsiveness to signals performed best in smaller hosts where infections have a lower pathogen group size. In this system, low signal abundance appears to limit fitness in hosts, while the optimal level of response to signals varies in different host environments.

The CRISPR-Cas systems were selectively inactivated during evolution of Bacillus cereus group for adaptation to diverse environments

CRISPR-Cas systems are considered as barriers to horizontal gene transfer (HGT). However, the influence of such systems on HGT within species is unclear. Also, little is known about the impact of CRISPR-Cas systems on bacterial evolution at the population level. Here, using Bacillus cereus sensu lato as model, we investigate the interplay between CRISPR-Cas systems and HGT at the population scale. We found that only a small fraction of the strains have CRISPR-Cas systems (13.9% of 1871), and most of such systems are defective based on their gene content analysis. Comparative genomic analysis revealed that the CRISPR-Cas systems are barriers to HGT within this group, since strains harboring active systems contain less mobile genetic elements (MGEs), have lower fraction of unique genes and also display limited environmental distributions than strains without active CRISPR-Cas systems. The introduction of a functional CRISPR-Cas system into a strain lacking the system resulted in reduced adaptability to various stresses and decreased pathogenicity of the transformant strain, indicating that B. cereus group strains could benefit from inactivating such systems. Our work provides a large-scale case to support that the CRISPR-Cas systems are barriers to HGT within species, and that in the B. cereus group the inactivation of CRISPR-Cas systems correlated with acquisition of MGEs that could result in better adaptation to diverse environments.

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.

Crystal toxins and the volunteer's dilemma in bacteria

The growth and virulence of the bacteria Bacillus thuringiensis depend on the production of Cry toxins, which are used to perforate the gut of its host. Successful invasion of the host relies on producing a threshold amount of toxin, after which there is no benefit from producing more toxin. Consequently, the production of Cry toxin appears to be a different type of social problem compared with the public goods scenarios that bacteria usually encounter. We show that selection for toxin production is a volunteer's dilemma. We make specific predictions that (a) selection for toxin production depends upon an interplay between the number of bacterial cells that each host ingests and the genetic relatedness between those cells; (b) cheats that do not produce toxin gain an advantage when at low frequencies, and at high bacterial density, allowing them to be maintained in a population alongside toxin‐producing cells. More generally, our results emphasize the diversity of the social games that bacteria play.

Interaction between Insects, Toxins, and Bacteria: Have We Been Wrong So Far?

Toxins are a major virulence factor produced by many pathogenic bacteria. In vertebrates, the response of hosts to the bacteria is inseparable from the response to the toxins, allowing a comprehensive understanding of this tripartite host-pathogen-toxin interaction. However, in invertebrates, this interaction has been investigated by two complementary but historically distinct fields of research: toxinology and immunology. In this article, I highlight how such dichotomy between these two fields led to a biased, or even erroneous view of the ecology and evolution of the interaction between insects, toxins, and bacteria. I focus on the reason behind such a dichotomy, on how to bridge the fields together, and on confounding effects that could bias the outcome of the experiments. Finally, I raise four questions at the border of the two fields on the cross-effects between toxins, bacteria, and spores that have been largely underexplored to promote a more comprehensive view of this interaction.

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.

c-di-GMP Regulates Various Phenotypes and Insecticidal Activity of Gram-Positive Bacillus thuringiensis

C-di-GMP has been well investigated to play significant roles in the physiology of many Gram-negative bacteria. However, its effect on Gram-positive bacteria is less known. In order to more understand the c-di-GMP functions in Gram-positive bacteria, we have carried out a detailed study on the c-di-GMP-metabolizing enzymes and their physiological functions in Bacillus thuringiensis, a Gram-positive entomopathogenic bacterium that has been applied as an insecticide successfully. We performed a systematic study on the ten putative c-di-GMP-synthesizing enzyme diguanylate cyclases (DGCs) and c-di-GMP-degrading enzyme phosphodiesterases (PDEs) in B. thuringiensis BMB171, and artificially elevated the intracellular c-di-GMP level in BMB171 by deleting one or more pde genes. We found increasing level of intracellular c-di-GMP exhibits similar activities as those in Gram-negative bacteria, including altered activities in cell motility, biofilm formation, and cell-cell aggregation. Unexpectedly, we additionally found a novel function exhibited by the increasing level of c-di-GMP to promote the insecticidal activity of this bacterium against Helicoverpa armigera. Through whole-genome transcriptome profile analyses, we found that 4.3% of the B. thuringiensis genes were differentially transcribed when c-di-GMP level was increased, and 77.3% of such gene products are involved in some regulatory pathways not reported in other bacteria to date. In summary, our study represents the first comprehensive report on the c-di-GMP-metabolizing enzymes, their effects on phenotypes, and the transcriptome mediated by c-di-GMP in an important Gram-positive bacterium.

An appeal for a more evidence based approach to biopesticide safety in the EU

EFSA responded to our perspective article on the safe use of the insect pathogen Bacillus thrurigiensis (Bt). In doing so they admitted that there is no direct evidence to suggest that B. thuringiensis can cause diarrhoea. They nevertheless continue to repeat the assertion that Bt cannot be distinguished from Bacillus cereus, even though nearly all Bt strains, and certainly all biopesticide strains, can be distinguished from B. cereus using multi-locus sequencing typing. EFSA also continue to repeat the unsupported and speculative hypothesis that Bt strains could be capable of causing cryptic infections in humans. This hypothesis is very much against the weight of all available safety and epidemiological data. Moreover, genotyping schemes of B. cereus group clinical infections also show that biopesticide strains have never been associated with human infections. Our position that Bt biopesticides and Bt isolates from the clade dominated by invertebrate pathogens are incapable of causing infections in humans is well supported by the international community of scientists familiar with the data on the safety of Bt after more than four decades of extensive use in agriculture and forestry.

Spores and vegetative cells of phenotypically and genetically diverse Bacillus cereus sensu lato are common bacteria in fresh water of northeastern Poland

Gram-positive rods Bacillus cereus sensu lato (sl) are common in natural habitats and food products. It is believed that they are restricted to spores; however, their ecology in aquatic habitats is still poorly investigated. Thus, the aim of the study was to assess the rain-dependent fluctuations in the concentration of B. cereus sl vegetative cells and spores, with evaluation of their phylogenetic and population structure in relation to the toxicity and psychrotolerance. We proved that vegetative cells of B. cereus sl are widely distributed in fresh water of rivers and lakes, being as common as spores. Moreover, heavy rain has a huge impact on their concentration in undisturbed environments. The diversity of B. cereus sl reflects the multiple sources of bacteria and the differences between their distinct environments. Next, their diverse genetic structure and phenotypes better fit their ecological properties than their taxonomic affiliation.

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.

Novel motB as a potential predictive tool for identification of B. cereus, B. thuringiensis and differentiation from other Bacillus species by triplex real-time PCR

Quantitative triplex real-time PCR (qPCR) offers an alternative method for detection of bacterial contamination. It provides quantitation of the number of gene copies. In our study, we established a qPCR assay to detect and quantify the specificity towards Bacillus cereus and B. thuringiensis. The assay was designed to detect a 280 bp fragment of motB gene encoding the flagellar motor protein, specific for detection of B. cereus and B. thuringiensis, excluding other group species B. pseudomycoides, B. mycoides and B. weihenstephanensis. Specificity of the assay was confirmed with 111 strains belonging to Bacillus cereus group and performed against 58 B. cereus, 50 B. thuringiensis, 3 other Bacillus bacteria and 9 non-Bacillus bacteria. Detection limit was determined for each assay. Direct analysis of samples revealed the specificity towards identification and characterization of B. cereus group cultured in nutrient media. Based on results, it was observed that motB showed 97% specificity towards B. cereus strains, 98% for B. thuringiensis but other B. cereus group showed less sensitivity (0%), thus, provides an efficient tool to identify B. cereus and B. thuringiensis. Further, environmental and food samples do not require band isolation, re-amplification or sequence identification. Thus, reducing the time and cost of analysis.

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.

Outer Membrane Vesicles (OMVs) of Gram-negative Bacteria: A Perspective Update

Outer Membrane Vesicles (OMVs) of Gram-negative bacteria are spherical membrane-enclosed entities of endocytic origin. Reported in the consortia of different bacterial species, production of OMVs into extracellular milieu seems essential for their survival. Enriched with bioactive proteins, toxins, and virulence factors, OMVs play a critical role in the bacteria-bacteria and bacteria-host interactions. Emergence of OMVs as distinct cellular entities helps bacteria in adaptating to diverse niches, in competing with other bacteria to protect members of producer species and more importantly play a crucial role in host-pathogen interaction. Composition of OMV, their ability to modulate host immune response, along with coordinated secretion of bacterial effector proteins, endows them with the armory, which can withstand hostile environments. Study of the OMV production under natural and diverse stress conditions has broadened the horizons, and also opened new frontiers in delineating the molecular machinery involved in disease pathogenesis. Playing diverse biological and pathophysiological functions, OMVs hold a great promise in enabling resurgence of bacterial diseases, in concomitance with the steep decline in the efficiency of antibiotics. Having multifaceted role, their emergence as a causative agent for a series of infectious diseases increases the probability for their exploitation in the development of effective diagnostic tools and as vaccines against diverse pathogenic species of Gram-negative origin.

The global regulator CodY is required for the fitness of Bacillus cereus in various laboratory media and certain beverages

The impact of gene mutations on the growth of the cells can be studied using pure cultures. However, the importance of certain proteins and pathways can be also examined via co-culturing wild type and its mutant derivative. Here, the relative fitness of a mutant strain that lacks the global nitrogen regulator, CodY, was examined in Bacillus cereus, a food poisoning Gram-positive bacterium. Fitness measurements revealed that the ΔcodY strain was outcompeted when cocultured with the wild-type ATCC 14579 under various rich laboratory medium, and also when inoculated in certain beverages. In nutrient-poor minimal medium, the ΔcodY mutant had comparable fitness to the wild-type strain. Interestingly, the relative fitness of the ΔcodY strain was antagonistic when it was cultivated in apple or orange juices due to unknown properties of these beverages, highlighting the importance of chemical composition of the test medium during the bacterial fitness measurements.

Limiting opportunities for cheating stabilizes virulence in insect parasitic nematodes

Cooperative secretion of virulence factors by pathogens can lead to social conflict when cheating mutants exploit collective secretion, but do not contribute to it. If cheats outcompete cooperators within hosts, this can cause loss of virulence. Insect parasitic nematodes are important biocontrol tools that secrete a range of significant virulence factors. Critically, effective nematodes are hard to maintain without live passage, which can lead to virulence attenuation. Using experimental evolution, we tested whether social cheating might explain unstable virulence in the nematode Heterorhabditis floridensis by manipulating relatedness via multiplicity of infection (MOI), and the scale of competition. Passage at high MOI, which should reduce relatedness, led to loss of fitness: virulence and reproductive rate declined together and all eight independent lines suffered premature extinction. As theory predicts, relatedness treatments had more impact under stronger global competition. In contrast, low MOI passage led to more stable virulence and increased reproduction. Moreover, low MOI lineages showed a trade-off between virulence and reproduction, particularly for lines under stronger between-host competition. Overall, this study indicates that evolution of virulence theory is valuable for the culture of biocontrol agents: effective nematodes can be improved and maintained if passage methods mitigate possible social conflicts.

High instability of a nematicidal Cry toxin plasmid in Bacillus thuringiensis

In bacterial pathogens, virulence factors are often carried on plasmids and other mobile genetic elements, and as such, plasmid evolution is central in understanding pathogenicity. Bacillus thuringiensis is an invertebrate pathogen that uses plasmid-encoded crystal (Cry) toxins to establish infections inside the host. Our study aimed to quantify stability of two Cry toxin-encoding plasmids, BTI_23p and BTI_16p, under standard laboratory culturing conditions. These two plasmids are part of the genome of the B. thuringiensis strain MYBT18679, which is of particular interest because of its high pathogenicity towards nematodes. One of the plasmids, BTI_23p, was found to be highly unstable, with substantial loss occurring within a single growth cycle. Nevertheless, longer term experimental evolution in the absence of a host revealed maintenance of the plasmid at low levels in the bacterial populations. BTI_23p encodes two nematicidal Cry toxins, Cry21Aa2 and Cry14Aa1. Consistent with previous findings, loss of the plasmid abolished pathogenicity towards the nematode Caenorhabditis elegans, which could be rescued by addition of Cry21Aa2-expressing Escherichia coli. These results implicate BTI_23p as a plasmid that is required for successful infection, yet unstable when present at high frequency in the population, consistent with the role of Cry toxins as public goods.

Live to cheat another day: bacterial dormancy facilitates the social exploitation of β-lactamases

The breakdown of antibiotics by β-lactamases may be cooperative, since resistant cells can detoxify their environment and facilitate the growth of susceptible neighbours. However, previous studies of this phenomenon have used artificial bacterial vectors or engineered bacteria to increase the secretion of β-lactamases from cells. Here, we investigated whether a broad-spectrum β-lactamase gene carried by a naturally occurring plasmid (pCT) is cooperative under a range of conditions. In ordinary batch culture on solid media, there was little or no evidence that resistant bacteria could protect susceptible cells from ampicillin, although resistant colonies could locally detoxify this growth medium. However, when susceptible cells were inoculated at high densities, late-appearing phenotypically susceptible bacteria grew in the vicinity of resistant colonies. We infer that persisters, cells that have survived antibiotics by undergoing a period of dormancy, founded these satellite colonies. The number of persister colonies was positively correlated with the density of resistant colonies and increased as antibiotic concentrations decreased. We argue that detoxification can be cooperative under a limited range of conditions: if the toxins are bacteriostatic rather than bacteridical; or if susceptible cells invade communities after resistant bacteria; or if dormancy allows susceptible cells to avoid bactericides. Resistance and tolerance were previously thought to be independent solutions for surviving antibiotics. Here, we show that these are interacting strategies: the presence of bacteria adopting one solution can have substantial effects on the fitness of their neighbours.

Bacterial Cooperation Causes Systematic Errors in Pathogen Risk Assessment due to the Failure of the Independent Action Hypothesis

The Independent Action Hypothesis (IAH) states that pathogenic individuals (cells, spores, virus particles etc.) behave independently of each other, so that each has an independent probability of causing systemic infection or death. The IAH is not just of basic scientific interest; it forms the basis of our current estimates of infectious disease risk in humans. Despite the important role of the IAH in managing disease interventions for food and water-borne pathogens, experimental support for the IAH in bacterial pathogens is indirect at best. Moreover since the IAH was first proposed, cooperative behaviors have been discovered in a wide range of microorganisms, including many pathogens. A fundamental principle of cooperation is that the fitness of individuals is affected by the presence and behaviors of others, which is contrary to the assumption of independent action. In this paper, we test the IAH in Bacillus thuringiensis (B.t), a widely occurring insect pathogen that releases toxins that benefit others in the inoculum, infecting the diamondback moth, Plutella xylostella. By experimentally separating B.t. spores from their toxins, we demonstrate that the IAH fails because there is an interaction between toxin and spore effects on mortality, where the toxin effect is synergistic and cannot be accommodated by independence assumptions. Finally, we show that applying recommended IAH dose-response models to high dose data leads to systematic overestimation of mortality risks at low doses, due to the presence of synergistic pathogen interactions. Our results show that cooperative secretions can easily invalidate the IAH, and that such mechanistic details should be incorporated into pathogen risk analysis.

Biogenesis and multifaceted roles of outer membrane vesicles from Gram-negative bacteria

Outer membrane vesicles (OMVs) released from Gram-negative bacteria consist of lipids, proteins, lipopolysaccharides and other molecules. OMVs are associated with several biological functions such as horizontal gene transfer, intracellular and intercellular communication, transfer of contents to host cells, and eliciting an immune response in host cells. Although hypotheses have been made concerning the mechanism of biogenesis of these vesicles, research on OMV formation is far from complete. The roles of outer membrane components, bacterial quorum sensing molecules and some specific proteins in OMV biogenesis have been studied. This review discusses the different models that have been proposed for OMV biogenesis, along with details of the biological functions of OMVs and the likely scope of future research.

The genomes, proteomes, and structures of three novel phages that infect the Bacillus cereus group and carry putative virulence factors

This article reports the results of studying three novel bacteriophages, JL, Shanette, and Basilisk, which infect the pathogen Bacillus cereus and carry genes that may contribute to its pathogenesis. We analyzed host range and superinfection ability, mapped their genomes, and characterized phage structure by mass spectrometry and transmission electron microscopy (TEM). The JL and Shanette genomes were 96% similar and contained 217 open reading frames (ORFs) and 220 ORFs, respectively, while Basilisk has an unrelated genome containing 138 ORFs. Mass spectrometry revealed 23 phage particle proteins for JL and 15 for Basilisk, while only 11 and 4, respectively, were predicted to be present by sequence analysis. Structural protein homology to well-characterized phages suggested that JL and Shanette were members of the family Myoviridae, which was confirmed by TEM. The third phage, Basilisk, was similar only to uncharacterized phages and is an unrelated siphovirus. Cryogenic electron microscopy of this novel phage revealed a T=9 icosahedral capsid structure with the major capsid protein (MCP) likely having the same fold as bacteriophage HK97 MCP despite the lack of sequence similarity. Several putative virulence factors were encoded by these phage genomes, including TerC and TerD involved in tellurium resistance. Host range analysis of all three phages supports genetic transfer of such factors within the B. cereus group, including B. cereus, B. anthracis, and B. thuringiensis. This study provides a basis for understanding these three phages and other related phages as well as their contributions to the pathogenicity of B. cereus group bacteria. Importance: The Bacillus cereus group of bacteria contains several human and plant pathogens, including B. cereus, B. anthracis, and B. thuringiensis. Phages are intimately linked to the evolution of their bacterial hosts and often provide virulence factors, making the study of B. cereus phages important to understanding the evolution of pathogenic strains. Herein we provide the results of detailed study of three novel B. cereus phages, two highly related myoviruses (JL and Shanette) and an unrelated siphovirus (Basilisk). The detailed characterization of host range and superinfection, together with results of genomic, proteomic, and structural analyses, reveal several putative virulence factors as well as the ability of these phages to infect different pathogenic species.

Division of labour and terminal differentiation in a novel Bacillus thuringiensis strain

A major challenge in bacterial developmental biology has been to understand the mechanisms underlying cell fate decisions. Some differentiated cell types display cooperative behaviour. Cooperation is one of the greatest mysteries of evolutionary biology and microbes have been considered as an excellent system for experimentally testing evolution theories. Bacillus thuringiensis (Bt) is a spore-forming bacterium, which is genetically closely related to B. anthracis, the agent of anthrax, and to B. cereus, an opportunistic human pathogen. The defining feature that distinguishes Bt from its relatives is its ability to produce crystal inclusions in the sporulating cells. These toxins are solubilized after ingestion and are cooperative public goods in insect hosts. In this study, we describe a Bt strain LM1212 that presents the unique ability to terminally differentiate into crystal producers and spore formers. Transcriptional analysis based on lacZ and gfp reporter genes suggested that this phenotype is the consequence of a new type of cell differentiation associated with a novel regulation mode of cry gene expression. The differentiating crystal-producer phenotype has higher spore productivity than a typical Bt strain and is better able to compete with Cry toxin null 'cheaters'. Potentially, this division of labour provides additional fitness benefits in terms of spore viability or durability of Cry toxin.

Regulation of cry gene expression in Bacillus thuringiensis

Bacillus thuringiensis differs from the closely related Bacillus cereus group species by its ability to produce crystalline inclusions. The production of these crystals mainly results from the expression of the cry genes, from the stability of their transcripts and from the synthesis, accumulation and crystallization of large amounts of insecticidal Cry proteins. This process normally coincides with sporulation and is regulated by various factors operating at the transcriptional, post-transcriptional, metabolic and post-translational levels.