Genomic characterization of the Bacillus cereus sensu lato species: backdrop to the evolution of Bacillus anthracis

In silico and structural analysis of Bacillus licheniformis FAO.CP7 pullulanase isolated from cocoa (Theobroma cacao L.) pod waste

Pullulanase (EC 3.2.1.41) is an important debranching enzyme that plays a critical role in maximizing the abundant energy present in branched polysaccharides. Its unique ability to efficiently degrade branched polysaccharides makes it crucial in industries like biofuels, food, and pharmaceuticals. Therefore, discovering microbes that produce pullulanase and thrive in harsh industrial conditions holds significant potential for optimizing large-scale bioprocessing. This unique property has made pullulanase an important enzyme in the industry. Thus, the search for microbes that have the pullulanase production properties and capacity to withstand harsh industrial conditions will be of high industrial relevance. Therefore, this study aimed to amplify, sequence, and molecularly characterize the pullulanase gene encoding extracellular pullulanase in Bacillus licheniformis strain FAO.CP7 (Accession No: MN150530.1.) which was obtained from cocoa pods using several bioinformatics tools. The amplified PulA gene had a nucleotide sequence of 2247 base pairs encoding a full-length open reading frame (ORF) pullulanase protein of 748 amino-acids residues with molecular weight 82.39 kDa and theoretical isoelectric point of 6.47, respectively. The deduced pullulanase protein had an aliphatic index of 77.66. Using BLASTp, the deduced amino acid sequence of the pullulanase gene showed 85% homologies with those from B. licheniformis strains. Multiple sequence alignment of PulA protein sequence showed that it contains YNWGYNP motif which is also found in all type I pullulanase protein sequences analysed. The restriction mapping of the gene showed that it can be digested with several restriction enzymes. Further analysis revealed that the deduced protein had a hydrophobicity score of - 0.37 without a transmembrane helix. Overall, this study revealed the PulA gene of B. licheniformis strain FAO.CP7 obtained from cocoa pods and its deduced protein show significant potential for enhancing starch bioprocessing. With further optimization, it could offer substantial benefits to starch-based biotechnological industries.

Bacillus bombysepticus, an Entomopathogen in Yellow-Spined Bamboo Locust with Biocontrol Potential

Ceracris kiangsu (Tsai) (Orthoptera: Acrididae) is a severe pest of the bamboo Dendrocalamus barbatus Hsuch & D.Z. Li that is widely grown in plantations in Vietnam. A biocontrol option to manage this pest is a priority for the forest sector. This study evaluates whether there are potential biocontrol candidates within the C. kiangsu population. Bacteria were isolated from C. kiangsu nymphs that were parasitized in the field. In a nursery experiment, three Bacillus bombysepticus Wang isolates (FPRC17, FPRC23 and FPRC30) caused 63-68% mortality of C. kiangsu nymphs 9 days after spraying, and the results were similar to applying a commercial Bacillus thuringiensis product. In a field trial, feeding damage from a C. kiangsu infestation was reduced by over 70% at 15 days after spraying these isolates. This finding reveals the potential for the development of a natural biopesticide for bamboo locust management in Vietnam.

Investigating Anthrax-Associated Virulence Genes among Archival and Contemporary Bacillus cereus Group Genomes

Bacillus anthracis causes anthrax through virulence factors encoded on two plasmids. However, non-B. anthracis organisms within the closely related, environmentally ubiquitous Bacillus cereus group (BCG) may cause an anthrax-like disease in humans through the partial adoption of anthrax-associated virulence genes, challenging the definition of anthrax disease. To elucidate these phenomena and their evolutionary past, we performed whole-genome sequencing on non-anthracis BCG isolates, including 93 archival (1967-2003) and 5 contemporary isolates (2019-2023). We produced annotated genomic assemblies and performed a pan-genome analysis to identify evidence of virulence gene homology and virulence gene acquisition by linear inheritance or horizontal gene transfer. At least one anthrax-associated virulence gene was annotated in ten isolates. Most homologous sequences in archival isolates showed evidence of pseudogenization and subsequent gene loss. The presence or absence of accessory genes, including anthrax-associated virulence genes, aligned with the phylogenetic structure of the BCG core genome. These findings support the hypothesis that anthrax-associated virulence genes were inherited from a common ancestor in the BCG and were retained or lost across different lineages, and contribute to a growing body of work informing public health strategies related to anthrax surveillance and identification.

Rapid evolutionary tuning of endospore quantity versus quality trade-off via a phase-variable contingency locus

The UV resistance of bacterial endospores is an important quality supporting their survival in inhospitable environments and therefore constitutes an essential driver of the ecological success of spore-forming bacteria. Nevertheless, the variability and evolvability of this trait are poorly understood. In this study, directed evolution and genetics approaches revealed that the Bacillus cereus pdaA gene (encoding the endospore-specific peptidoglycan-N-acetylmuramic acid deacetylase) serves as a contingency locus in which the expansion and contraction of short tandem repeats can readily compromise (PdaAOFF) or restore (PdaAON) the pdaA open reading frame. Compared with B. cereus populations in the PdaAON state, populations in the PdaAOFF state produced a lower yield of viable endospores but endowed them with vastly increased UV resistance. Moreover, selection pressures based on either quantity (i.e., yield of viable endospores) or quality (i.e., UV resistance of viable endospores) aspects could readily shift populations between PdaAON and PdaAOFF states, respectively. Bioinformatic analysis also revealed that pdaA homologs within the Bacillus and Clostridium genera are often equipped with several short tandem repeat regions, suggesting a wider implementation of the pdaA-mediated phase variability in other sporeformers as well. These results for the first time reveal (1) pdaA as a phase-variable contingency locus in the adaptive evolution of endospore properties and (2) bet-hedging between what appears to be a quantity versus quality trade-off in endospore crops.

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

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

Whole genome sequence analyses of thermotolerant Bacillus sp. isolates from food

The Bacillus cereus group, also known as B. cereus sensu lato (B. cereus s.l.), is composed of various Bacillus species, some of which can cause diarrheal or emetic food poisoning. Several emerging highly heat-resistant Bacillus species have been identified, these include B. thermoamylovorans, B. sporothermodurans, and B. cytotoxicus NVH 391-98. Herein, we performed whole genome analysis of two thermotolerant Bacillus sp. isolates, Bacillus sp. B48 and Bacillus sp. B140, from an omelet with acacia leaves and fried rice, respectively. Phylogenomic analysis suggested that Bacillus sp. B48 and Bacillus sp. B140 are closely related to B. cereus and B. thuringiensis, respectively. Whole genome alignment of Bacillus sp. B48, Bacillus sp. B140, mesophilic strain B. cereus ATCC14579, and thermophilic strain B. cytotoxicus NVH 391-98 using the Mauve program revealed the presence of numerous homologous regions including genes responsible for heat shock in the dnaK gene cluster. However, the presence of a DUF4253 domain-containing protein was observed only in the genome of B. cereus ATCC14579 while the intracellular protease PfpI family was present only in the chromosome of B. cytotoxicus NVH 391-98. In addition, prophage Clp protease-like proteins were found in the genomes of both Bacillus sp. B48 and Bacillus sp. B140 but not in the genome of B. cereus ATCC14579. The genomic profiles of Bacillus sp. isolates were identified by using whole genome analysis especially those relating to heat-responsive gene clusters. The findings presented in this study lay the foundations for subsequent studies to reveal further insights into the molecular mechanisms of Bacillus species in terms of heat resistance mechanisms.

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.

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.

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

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

Adaptation and phenotypic diversification of Bacillus thuringiensis biofilm are accompanied by fuzzy spreader morphotypes

Bacillus cereus group (Bacillus cereus sensu lato) has a diverse ecology, including various species that produce biofilms on abiotic and biotic surfaces. While genetic and morphological diversification enables the adaptation of multicellular communities, this area remains largely unknown in the Bacillus cereus group. In this work, we dissected the experimental evolution of Bacillus thuringiensis 407 Cry- during continuous recolonization of plastic beads. We observed the evolution of a distinct colony morphotype that we named fuzzy spreader (FS) variant. Most multicellular traits of the FS variant displayed higher competitive ability versus the ancestral strain, suggesting an important role for diversification in the adaptation of B. thuringiensis to the biofilm lifestyle. Further genetic characterization of FS variant revealed the disruption of a guanylyltransferase gene by an insertion sequence (IS) element, which could be similarly observed in the genome of a natural isolate. The evolved FS and the deletion mutant in the guanylyltransferase gene (Bt407ΔrfbM) displayed similarly altered aggregation and hydrophobicity compared to the ancestor strain, suggesting that the adaptation process highly depends on the physical adhesive forces.

Bioplastic accumulates antibiotic and metal resistance genes in coastal marine sediments

The oceans are increasingly polluted with plastic debris, and several studies have implicated plastic as a reservoir for antibiotic resistance genes and a potential vector for antibiotic-resistant bacteria. Bioplastic is widely regarded as an environmentally friendly replacement to conventional petroleum-based plastic, but the effects of bioplastic pollution on marine environments remain largely unknown. Here, we present the first evidence that bioplastic accumulates antibiotic resistance genes (ARGs) and metal resistance genes (MRGs) in marine sediments. Biofilms fouling ceramic, polyethylene terephthalate (PET), and polyhydroxyalkanoate (PHA) were investigated by shotgun metagenomic sequencing. Four ARG groups were more abundant in PHA: trimethoprim resistance (TMP), multidrug resistance (MDR), macrolide-lincosamide-streptogramin resistance (MLS), and polymyxin resistance (PMR). One MRG group was more abundant in PHA: multimetal resistance (MMR). The relative abundance of ARGs and MRGs were strongly correlated based on a Mantel test between the Bray-Curtis dissimilarity matrices (R = 0.97, p < 0.05) and a Pearson's analysis (R = 0.96, p < 0.05). ARGs were detected in more than 40% of the 57 metagenome-assembled genomes (MAGs) while MRGs were detected in more than 90% of the MAGs. Further investigation (e.g., culturing, genome sequencing, antibiotic susceptibility testing) revealed that PHA biofilms were colonized by hemolytic Bacillus cereus group bacteria that were resistant to beta-lactams, vancomycin, and bacitracin. Taken together, our findings indicate that bioplastic, like conventional petroleum-based plastic, is a reservoir for resistance genes and a potential vector for antibiotic-resistant bacteria in coastal marine sediments.

Genotype and virulence gene analyses of Bacillus cereus group clinical isolates from the Chinese softshell turtle (Pelodiscus sinensis) in Taiwan

Chinese softshell turtles (Pelodiscus sinensis) (CST) are susceptible to infections by bacteria belonging to the Bacillus cereus group (Bcg). Bcg includes several closely related species, two of which, B. cereus and B. thuringiensis, are pathogens of aquatic animals or insects. In the present study, we collected 57 Bcg isolates obtained from diseased CST from 2016 to 2019 in Kaohsiung and Pingtung, the areas with the most CST farms in Taiwan. All isolates were divided into four genotypes with two restriction enzymes, SmaI and NotI, by pulsed-field gel electrophoresis and enterobacterial repetitive intergenic consensus-polymerase chain reaction (ERIC-PCR). Representative isolates from each genotype were subjected to phylogenetic tree analysis using 16S rDNA and pyruvate carboxylase genes as phylogenetic markers, and these CST isolates appeared in different clades. PCR was performed targeting six selected virulence genes, four of which were detected in CST isolates, including cytotoxin K (1/57), hblC of the haemolysin BL complex (46/57), nheA of the non-haemolytic enterotoxin complex (52/57) and enterotoxin FM (57/57), whereas cereulide synthetase and cereulide peptide synthase-like genes were not detected in any isolates.

Comparative genome analysis of Bacillus thuringiensis strain HD521 and HS18-1

Bacillus thuringiensis (Bt) is an important biological insecticide used to management of different agricultural pests by producing toxic parasporal crystals proteins. Strain HD521 has an antagonistic effect against Rhizoctonia solani AG1IA, the causal agent of rice sheath blight. This strain with three cry7 genes can the formation of bipyramidal parasporal crystals (BPCs). BPCs are used for insecticidal activities against Henosepilachna vigintioctomaculata larva (Coleoptera). Strain HS18-1 contains different types of BPCs encoding genes and has effective toxicity for Lepidoptera and Diptera insects. Here we report the whole genome sequencing and assembly of HD521 and HS18-1 strains and analyzed the genome constitution covering virulence factors, types of plasmid, insertion sequences, and prophage sequences. The results showed that the genome of strain HD521 contains a circular chromosome and six circular plasmids, encoding eight types of virulence protein factors [Immune Inhibitor A, Hemolytic Enterotoxin, S-layer protein, Phospholipase C, Zwittermicin A-resistance protein, Metalloprotease, Chitinase, and N-acyl homoserine lactonase (AiiA)], four families of insertion sequence, and comprises six pro-phage sequences. The genome of strain HS18-1 contains one circular chromosome and nine circular plasmids, encoding five types of virulence protein factors [Hemolytic Enterotoxin, S-layer protein, Phospholipase C, Chitinase, and N-acyl homoserine lactonase (AiiA)] and four families of insertion sequence, and comprises of three pro-phage sequences. The obtained results will contribute to deeply understand the B. thuringiensis strain HD521 and HS18-1 at the genomic level.

Plant Polysaccharides Modulate Biofilm Formation and Insecticidal Activities of Bacillus thuringiensis Strains

After the biological pesticide Bacillus thuringiensis (Bt) is applied to the field, it has to remain on the surface of plants to have the insecticidal activities against insect pests. Bt can form biofilms on the surface of vegetable leaves, which were rich in polysaccharides. However, the relationship between polysaccharides of the leaves and the biofilm formation as well as the insecticidal activities of Bt is still unknown. Herein, this study focused on the effects of plant polysaccharides pectin and xylan on biofilm formation and the insecticidal activities of Bt strains. By adding pectin, there were 88 Bt strains with strong biofilm formation, 69 strains with weak biofilm formation, and 13 strains without biofilm formation. When xylan was added, 13 Bt strains formed strong biofilms, 98 strains formed weak biofilms, and 59 strains did not form biofilms. This indicated that two plant polysaccharides, especially pectin, modulate the biofilm formation of Bt strains. The ability of pectin to induce biofilm formation was not related to Bt serotypes. Pectin promoted the biofilms formed by Bt cells in the logarithmic growth phase and lysis phase at the air–liquid interface, while it inhibited the biofilms formed by Bt cells in the sporangial phase at the air–liquid interface. The dosage of pectin was positively correlated with the yield of biofilms formed by Bt cells in the logarithmic growth phase or lysis phase at the solid–liquid interfaces. Pectin did not change the free-living growth and the cell motility of Bt strains. Pectin can improve the biocontrol activities of the spore–insecticidal crystal protein mixture of Bt and BtK commercial insecticides, as well as the biofilms formed by the logarithmic growth phase or lysis phase of Bt cells. Our findings confirmed that plant polysaccharides modulate biofilm formation and insecticidal activities of Bt strains and built a foundation for the construction of biofilm-type Bt biopesticides.

A putative exosporium lipoprotein GBAA0190 of Bacillus anthracis as a potential anthrax vaccine candidate

BACKGROUND

Bacillus ancthracis causes cutaneous, pulmonary, or gastrointestinal forms of anthrax. B. anthracis is a pathogenic bacterium that is potentially to be used in bioterrorism because it can be produced in the form of spores. Currently, protective antigen (PA)-based vaccines are being used for the prevention of anthrax, but it is necessary to develop more safe and effective vaccines due to their prolonged immunization schedules and adverse reactions.

METHODS

We selected the lipoprotein GBAA0190, a potent inducer of host immune response, present in anthrax spores as a novel potential vaccine candidate. Then, we evaluated its immune-stimulating activity in the bone marrow-derived macrophages (BMDMs) using enzyme-linked immunosorbent assay (ELISA) and Western blot analysis. Protective efficacy of GBAA0190 was evaluated in the guinea pig (GP) model.

RESULTS

The recombinant GBAA0190 (r0190) protein induced the expression of various inflammatory cytokines including tumor necrosis factor-α (TNF-α), interleukin-6 (IL-6), monocyte chemoattractant protein-1 (MCP-1), and macrophage inflammatory protein-1α (MIP-1α) in the BMDMs. These immune responses were mediated through toll-like receptor 1/2 via activation of mitogen-activated protein (MAP) kinase and Nuclear factor-κB (NF-κB) pathways. We demonstrated that not only immunization of r0190 alone, but also combined immunization with r0190 and recombinant PA showed significant protective efficacy against B. anthracis spore challenges in the GP model.

CONCLUSIONS

Our results suggest that r0190 may be a potential target for anthrax vaccine.

Phylogenetic Methods for Genome-Wide Association Studies in Bacteria

Genome-wide association studies in bacteria have great potential to deliver a better understanding of the genetic basis of many biologically important phenotypes, including antibiotic resistance, pathogenicity, and host adaptation. Such studies need however to account for the specificities of bacterial genomics, especially in terms of population structure, homologous recombination, and genomic plasticity. A powerful way to tackle this challenge is to use a phylogenetic approach, which is based on long-standing methodology for the evolutionary analysis of bacterial genomic data. Here we present both the theoretical and practical aspects involved in the use of phylogenetic methods for bacterial genome-wide association studies.

A fine-scale map of genome-wide recombination in divergent Escherichia coli population

Recombination is one of the most important molecular mechanisms of prokaryotic genome evolution, but its exact roles are still in debate. Here we try to infer genome-wide recombination within a species, utilizing a dataset of 149 complete genomes of Escherichia coli from diverse animal hosts and geographic origins, including 45 in-house sequenced with the single-molecular real-time platform. Two major clades identified based on physiological, clinical and ecological characteristics form distinct genetic lineages based on scarcity of interclade gene exchanges. By defining gene-based syntenies for genomic segments within and between the two clades, we build a fine-scale recombination map for this representative global E. coli population. The map suggests extensive within-clade recombination that often breaks physical linkages among individual genes but seldom interrupts the structure of genome organizational frameworks as well as primary metabolic portfolios supported by the framework integrity, possibly due to strong natural selection for both physiological compatibility and ecological fitness. In contrast, the between-clade recombination declines drastically when phylogenetic distance increases to the extent where a 10-fold reduction can be observed, establishing a firm genetic barrier between clades. Our empirical data suggest a critical role for such recombination events in the early stage of speciation where recombination rate is associated with phylogenetic distance in addition to sequence and gene variations. The extensive intraclade recombination binds sister strains into a quasisexual group and optimizes genes or alleles to streamline physiological activities, whereas the sharply declined interclade recombination split the population into clades adaptive to divergent ecological niches.

Development of a visible loop mediated isothermal amplification assay for rapid detection of Bacillus anthracis

Distressing effects on animal and human health with lethal progression, being used as bioweapon and shared features with non-pathogenic bacteria demands sensitive, specific, safe, cost effective and rapid detection methods for anthrax causing organisms. Conventional microbiology based diagnostics for anthrax are time consuming and need sophisticated equipment, while molecular diagnostics require less time and labor. The Loop mediated isothermal amplification assay (LAMP) is rapid, sensitive and specific assay and requires no specialized equipment. In the present study, we developed a LAMP assay for rapid as well as specific detection of Bacillus anthracis. The optimized assay produced positive results with the Sterne strain and one field isolate of B. anthracis and, negative results with other bacteria of the same and different genera within 2 h. Sensitivity was 500 fg of total DNA of B. anthracis, which was 100 times more sensitive than conventional PCR. The present study also demonstrated that the simple method of total DNA extraction by repeated boiling and freezing will not adversely affect the LAMP results. In conclusion, the optimized LAMP assay is a promising tool for the specific, sensitive, less time-consuming diagnosis for anthrax causing bacteria and also, for detecting the virulence of suspected B. anthracis cultures.

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.

Detection of Bacillus anthracis and Bacillus anthracis-like spores in soil from state of Rio de Janeiro, Brazil

BACKGROUND

Bacillus anthracis is the aetiologic agent of anthrax, a re-emerging, septicaemic, haemorrhagic and lethal disease that affects humans, domestic ruminants and wildlife. Plasmids pXO1 and pXO2 are attributes that confer pathogenicity to B. anthracis strains. This bacterium was used as biological weapon in the World Wars and in the biological attack in the United States of America at 2001. B. anthracis is classified as a Tier 1 bioterrorism agent by the Centers for Diseases Control and Prevention. Anthrax is recognised as a re-emerging disease. Several studies concerning the dynamics of B. anthracis cycle in soil revealed that nonpathogenic B. anthracis strains due to lack of pXO2 plasmid are commonly found in some types of soil.

OBJECTIVES

This study aimed isolation and identification of B. anthracis spores in soil samples of the state of Rio de Janeiro, Brazil.

METHODS

Phenotypic and genotypic approaches were used to identify isolates including MALDI-TOF/MS, motility test, susceptibility to gamma phage and penicillin, survey for pag and cap genes as surrogates of pXO1 and pXO2 plasmids, respectively, and sequencing of 16SrRNA-encoding gene. Physicochemical analysis of the soil samples were carried out to describe soil characteristics.

FINDINGS

We observed the presence of one B. anthracis pXO1+ and pXO2- isolated from clay loam soil; one B. anthracis-like strain pXO1+ and pXO2-isolated from loamy sand; and 10 Bacillus spp. strains sensitive to phage-gamma that need better characterisation to define which their species were recovered from loamy sand.

MAIN CONCLUSIONS

This work showed promising results and it was the first study to report results from an active surveillance for B. anthracis in Brazil.

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.

Cereulide Synthetase Acquisition and Loss Events within the Evolutionary History of Group III Bacillus cereus Sensu Lato Facilitate the Transition between Emetic and Diarrheal Foodborne Pathogens

Cereulide-producing members of Bacillus cereussensu lato group III (also known as emetic B. cereus) possess cereulide synthetase, a plasmid-encoded, nonribosomal peptide synthetase encoded by the ces gene cluster. Despite the documented risks that cereulide-producing strains pose to public health, the level of genomic diversity encompassed by emetic B. cereus has never been evaluated at a whole-genome scale. Here, we employ a phylogenomic approach to characterize group III B. cereussensu lato genomes which possess ces (ces positive) alongside their closely related, ces-negative counterparts (i) to assess the genomic diversity encompassed by emetic B. cereus and (ii) to identify potential ces loss and/or gain events within the evolutionary history of the high-risk and medically relevant sequence type (ST) 26 lineage often associated with emetic foodborne illness. Using all publicly available ces-positive group III B. cereussensu lato genomes and the ces-negative genomes interspersed among them (n = 159), we show that emetic B. cereus is not clonal; rather, multiple lineages within group III harbor cereulide-producing strains, all of which share an ancestor incapable of producing cereulide (posterior probability = 0.86 to 0.89). Members of ST 26 share an ancestor that existed circa 1748 (95% highest posterior density [HPD] interval = 1246.89 to 1915.64) and first acquired the ability to produce cereulide before 1876 (95% HPD = 1641.43 to 1946.70). Within ST 26 alone, two subsequent ces gain events were observed, as well as three ces loss events, including among isolates responsible for B. cereussensu lato toxicoinfection (i.e., "diarrheal" illness).IMPORTANCEB. cereus is responsible for thousands of cases of foodborne disease each year worldwide, causing two distinct forms of illness: (i) intoxication via cereulide (i.e., emetic syndrome) or (ii) toxicoinfection via multiple enterotoxins (i.e., diarrheal syndrome). Here, we show that emetic B. cereus is not a clonal, homogenous unit that resulted from a single cereulide synthetase gain event followed by subsequent proliferation; rather, cereulide synthetase acquisition and loss is a dynamic, ongoing process that occurs across lineages, allowing some group III B. cereussensu lato populations to oscillate between diarrheal and emetic foodborne pathogens over the course of their evolutionary histories. We also highlight the care that must be taken when selecting a reference genome for whole-genome sequencing-based investigation of emetic B. cereussensu lato outbreaks, since some reference genome selections can lead to a confounding loss of resolution and potentially hinder epidemiological investigations.

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

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

Influence of the genus Bacillus bacteria on the content of H2O2 and the activity of hydrolases and their inhibitors in potato plants during Phytophthora infestans Mont. de Bary infection

The authors studied the effect of treatment with bacteria Bacillus subtilis Cohn (strains 26D) and B. thuringiensis Berliner (strain B-6066) on the hydrogen peroxide (H2O2) content, the activity of hydrolytic enzymes and their protein inhibitors in potato plants (Solanum tuberosum L.) in connection with development of resistance to the late blight pathogen - oomycete Phytophthora infestans Mont. de Bary. Studies were carried out on potato plants of the susceptible Early Rose potato cultivar that were treated with a suspension of B. subtilis and B. thuringiensis bacteria (108 cells/ml) and infected with P. infestans (107 spores/ml). A decrease in the degree of leaf damage by oomycete was revealed under the influence of the genus Bacillus bacteria, depending on the strain. The increase in potato resistance to P. infestans infection was mediated by the stimulating effect of the B. subtilis 26D and the B. thuringiensis B-6066 bacteria on the concentration of H2O2, the modulating effect on the activity of hydrolytic enzymes and the enhancement of the transcriptional activity of protease and amylase inhibitor genes in plant tissues. Differences in the degree of activation of the transcriptional activity of hydrolase inhibitor genes by the B. subtilis 26D and the B. thuringiensis B-6066 bacteria were revealed, which suggests differential ways of forming the potato resistance to P. infestans under their influence.

Comparative Genome Analysis of Bacillus sporothermodurans with Its Closest Phylogenetic Neighbor, Bacillus oleronius, and Bacillus cereus and Bacillus subtilis Groups

Bacillus sporothermodurans currently possesses one of the most highly heat-resistant spores (HRS), which can withstand ultra-high temperature (UHT) processing. Determination of multiple whole genome sequences of B. sporothermodurans provided an opportunity to perform the first comparative genome analysis between strains and with B. oleronius, B. cereus, and B. subtilis groups. In this study, five whole genome sequences of B. sporothermodurans strains, including those belonging to the HRS clone (SAD and BR12) normally isolated from UHT milk, were compared with the aforementioned Bacillus species for gene clusters responsible for heat resistance. In the phylogenomic analysis, B. sporothermodurans, with its closest phylogenetic neighbor, B. oleronius, clustered with B. thermoamylovorans and B. thermotolerans. Heat shock proteins GrpE, GroES, GroEL, and DnaK presented identical sequences for all B. sporothermodurans strains, indicating that differences in functional efficiency are not involved in the thermal resistance variations. However, comparing all species evaluated, B. sporothermodurans exhibited a different gene configuration in the chromosomal region of the heat shock protein GrpE. Furthermore, only B. sporothermodurans strains presented the stage II sporulation protein P gene located in this region. Multisequence alignment and phylogenetic analysis of the ClpB protein showed differences for HRS and non-HRS strains. The study identified ClpC, ClpE, and ClpX as the three ATPases putatively involved in protein disaggregation in B. sporothermodurans. Bacillussporothermodurans exhibits high homology with other Bacillus species in the DnaK, DnaJ, GroEL, and GroES cluster of genes involved in heat resistance. The data presented here pave the way to select and evaluate the phenotypic effects of genes putatively involved in heat resistance.

Novel xylanase producing Bacillus strain X2: molecular phylogenetic analysis and its application for production of xylooligosaccharides

A Bacillus strain X2 that produced extracellular endo-xylanase (GH 11) (EC: 3.2.1.8) was isolated from the soil of the Northeast India region. This aerobic culture was Gram positive and endospore forming. Chemotaxonomic characterization showed variance with the fatty acid profile of related species in the Bacillus subtilis group. In Bacillus strain X2, distinct occurrence of iso-C14:0 lipids is absent in other related species. The 16S rRNA gene sequence homology showed 99% similarity with Bacillus subtilis subsp. inaquosorum. The phylogenetic analysis by the multilocus sequence analysis (MLSA) of the nucleotide sequence of six concatenated genes (16S rRNA, groEL, gyrA, polC, purH and rpoB) resolved the taxonomic position of the Bacillus strain X2 in the Bacillus subtilis subsp. group. The MLSA showed that it is a member of a clade that includes Bacillus subtilis subsp. stercoris. In in silico DNA-DNA hybridization (DDH), the highest matching score was obtained with Bacillus subtilis subsp. stercoris (87%). The in silico DDH of the genome (G + C 43.7 mol %) shared 48.5%, with Bacillus subtilis subsp. inaquosorum. The MLSA phylogenetic tree and the highest degree of DNA hybridization, indicating that it belongs to the Bacillus subtilis subspecies stercoris.

A Unique Isolation of a Lytic Bacteriophage Infected Bacillus anthracis Isolate from Pafuri, South Africa

Bacillus anthracis is a soil-borne, Gram-positive endospore-forming bacterium and the causative agent of anthrax. It is enzootic in Pafuri, Kruger National Park in South Africa. The bacterium is amplified in a wild ungulate host, which then becomes a source of infection to the next host upon its death. The exact mechanisms involving the onset (index case) and termination of an outbreak are poorly understood, in part due to a paucity of information about the soil-based component of the bacterium’s lifecycle. In this study, we present the unique isolation of a dsDNA bacteriophage from a wildebeest carcass site suspected of having succumbed to anthrax. The aggressively lytic bacteriophage hampered the initial isolation of B. anthracis from samples collected at the carcass site. Classic bacteriologic methods were used to test the isolated phage on B. anthracis under different conditions to simulate deteriorating carcass conditions. Whole genome sequencing was employed to determine the relationship between the bacterium isolated on site and the bacteriophage-dubbed Bacillus phage Crookii. The 154,012 bp phage belongs to Myoviridae and groups closely with another African anthrax carcass-associated Bacillus phage WPh. Bacillus phage Crookii was lytic against B. cereus sensu lato group members but demonstrated a greater affinity for encapsulated B. anthracis at lower concentrations (<1 × 10⁸ pfu) of bacteriophage. The unusual isolation of this bacteriophage demonstrates the phage’s role in decreasing the inoculum in the environment and impact on the life cycle of B. anthracis at a carcass site.

Molecular and Genomic Characterization of PFAB2: A Non-virulent Bacillus anthracis Strain Isolated from an Indian Hot Spring

Background

Thermophilic bacilli in both aerobic or facultative anaerobic forms have been isolated for over a hundred years from different mesophilic or thermophilic environments as they are potential source of bioactive secondary metabolites. But the taxonomic resolution in the Bacillus genus at species or at strain level is very challenging for the insufficient divergence of the 16S rRNA genes. One such recurring problem is among Bacillus anthracis, B. cereus and B. thuringiensis. The disease-causing B. anthracis strains have their characteristic virulence factors coded in two well-known plasmids, namely pXO1 (toxin genes) and pXO2 (capsule genes).

Objective

The present study aimed at the molecular and genomic characterization of a recently reported thermophilic and environmental isolate of B. anthracis, strain PFAB2.

Methods

We performed comparative genomics between the PFAB2 genome and different strains of B. anthracis, along with closely related B. cereus strains.

Results

The pangenomic analysis suggests that the PFAB2 genome harbors no complete prophage genes. Cluster analysis of Bray-Kurtis similarity resemblance matrix revealed that gene content of PFAB2 is more closely related to other environmental strains of B. anthracis. The secretome analysis and the in vitro and in vivo pathogenesis experiments corroborate the avirulent phenotype of this strain. The most probable explanation for this phenotype is the apparent absence of plasmids harboring genes for capsule biosynthesis and toxins secretion in the draft genome. Additional features of PFAB2 are good spore-forming and germinating capabilities and rapid replication ability.

Conclusion

The high replication rate in a wide range of temperatures and culture media, the non-pathogenicity, the good spore forming capability and its genomic similarity to the Ames strain together make PFAB2 an interesting model strain for the study of the pathogenic evolution of B. anthracis.

Novel Sequence Type in Bacillus cereus Strains Associated with Nosocomial Infections and Bacteremia, Japan

This sequence type was dominant in isolates from bacteremia patients in 3 hospitals.

Bacillus cereus is associated with foodborne illnesses characterized by vomiting and diarrhea. Although some B. cereus strains that cause severe extraintestinal infections and nosocomial infections are recognized as serious public health threats in healthcare settings, the genetic backgrounds of B. cereus strains causing such infections remain unknown. By conducting pulsed-field gel electrophoresis and multilocus sequence typing, we found that a novel sequence type (ST), newly registered as ST1420, was the dominant ST isolated from the cases of nosocomial infections that occurred in 3 locations in Japan in 2006, 2013, and 2016. Phylogenetic analysis showed that ST1420 strains belonged to the Cereus III lineage, which is much closer to the Anthracis lineage than to other Cereus lineages. Our results suggest that ST1420 is a prevalent ST in B. cereus strains that have caused recent nosocomial infections in Japan.

Comparative Genomic and Phylogenomic Analyses Clarify Relationships Within and Between Bacillus cereus and Bacillus thuringiensis: Proposal for the Recognition of Two Bacillus thuringiensis Genomovars

The present study was designed to clarify the taxonomic status of two species classified as Bacillus cereus sensu lato, namely B. cereus sensu stricto and Bacillus thuringiensis. To this end, nearly 900 whole genome sequences of strains assigned to these taxa were the subject of comparative genomic and phylogenomic analyses. A phylogenomic tree based on core gene sequences showed that the type strains of B. cereus and B. thuringiensis formed a well-supported monophyletic clade that was clearly separated from corresponding clades composed of the remaining validly published species classified as B. cereus sensu lato. However, since average nucleotide identity and digital DNA-DNA hybridization similarities between the two types of Bacillus were slightly higher than the thresholds used to distinguish between closely related species we conclude that B. cereus and B. thuringiensis should continue to be recognized as validly published species. The B. thuringiensis strains were assigned to two genomically distinct groups, we propose that these taxa be recognized as genomovars, that is, as B. thuringiensis gv. thuringiensis and B. thuringiensis gv. cytolyticus. The extensive comparative genomic data clearly show that the distribution of pesticidal genes is irregular as strains identified as B. thuringiensis were assigned to several polyphyletic groups/subclades within the B. cereus-B. thuringiensis clade. Consequently, we recommend that genomic or equivalent molecular systematic features should be used to identify B. thuringiensis strains as the presence of pesticidal genes cannot be used as a diagnostic marker for this species. Comparative taxonomic studies are needed to find phenotypic properties that can be used to distinguish between the B. thuringiensis genomovars and between them and B. cereus.

Bacillus cereus cshA Is Expressed during the Lag Phase of Growth and Serves as a Potential Marker of Early Adaptation to Low Temperature and pH

The spore-forming bacterium B. cereus is a major cause of foodborne outbreaks in Europe. Some B. cereus strains can grow at low temperatures and low pH in many processed foods. Modeling of the bacterial lag time is hampered by a lack of knowledge of the timing of events occurring during this phase. In this context, the identification of lag phase markers, not currently available, could be a real advance for the better prediction of lag time duration. Currently, no molecular markers of this phase are available. By determining that cshA was always expressed early during the lag phase, we provide a molecular marker of the early adaptation process of B. cereus cells when exposed to low temperature and pH.

Bacterial adaptation is characterized by a lag phase during which cells do not multiply or modify their physiology to cope with the constraints of their environment. Our aim was to determine a sequence of events during the lag phase of growth at low temperature and pH for three Bacillus cereus strains. The onsets of expression of two genes, one of which is essential for stress adaptation (cshA, coding for a RNA helicase) and one of which is involved in the transition between lag phase and exponential phase (abrB, coding for a transition regulator), were determined using fluorescent transcriptional reporter systems. Regardless of the stressing conditions and the tested strains, the cshA promoter was active very early, while the biomass increased and always did so before the first cell division. At 12°C and pH 7.0, the onset of cshA promoter activity occurred at between 3 h and 7 h, while the bacterial counts started to increase at between 12 h and 13 h. At pH 5.0 and at 20°C or 30°C, the onset of cshA promoter activity occurred before 1 h and earlier than at pH 7.0. In contrast, the onset of abrB promoter activity depended on the strain and the stressing conditions. In the ATCC 14579 strain, the onset of abrB promoter activity always started at between 30 min and 3 h, before biomass increased and cell division occurred. For the other strains, it took place along with the first cell division at 12°C but did so much later during growth under the other tested conditions.

IMPORTANCE The spore-forming bacterium B. cereus is a major cause of foodborne outbreaks in Europe. Some B. cereus strains can grow at low temperatures and low pH in many processed foods. Modeling of the bacterial lag time is hampered by a lack of knowledge of the timing of events occurring during this phase. In this context, the identification of lag phase markers, not currently available, could be a real advance for the better prediction of lag time duration. Currently, no molecular markers of this phase are available. By determining that cshA was always expressed early during the lag phase, we provide a molecular marker of the early adaptation process of B. cereus cells when exposed to low temperature and pH.

In silico prospection of microorganisms to produce polyhydroxyalkanoate from whey: Caulobacter segnis DSM 29236 as a suitable industrial strain

Polyhydroxyalkanoates (PHAs) are polyesters of microbial origin that can be synthesized by prokaryotes from noble sugars or lipids and from complex renewable substrates. They are an attractive alternative to conventional plastics because they are biodegradable and can be produced from renewable resources, such as the surplus of whey from dairy companies. After an in silico screening to search for ß-galactosidase and PHA polymerase genes, several bacteria were identified as potential PHA producers from whey based on their ability to hydrolyse lactose. Among them, Caulobacter segnis DSM 29236 was selected as a suitable strain to develop a process for whey surplus valorization. This microorganism accumulated 31.5% of cell dry weight (CDW) of poly(3-hydroxybutyrate) (PHB) with a titre of 1.5 g l-1 in batch assays. Moreover, the strain accumulated 37% of CDW of PHB and 9.3 g l-1 in fed-batch mode of operation. This study reveals this species as a PHA producer and experimentally validates the in silico bioprospecting strategy for selecting microorganisms for waste re-valorization.

Pangenomic Approach To Understanding Microbial Adaptations within a Model Built Environment, the International Space Station, Relative to Human Hosts and Soil

Understanding underlying mechanisms involved in microbial persistence in the built environment (BE) is essential for strategically mitigating potential health risks. To test the hypothesis that BEs impose selective pressures resulting in characteristic adaptive responses, we performed a pangenomics meta-analysis leveraging 189 genomes (accessed from GenBank) of two epidemiologically important taxa, Bacillus cereus and Staphylococcus aureus, isolated from various origins: the International Space Station (ISS; a model BE), Earth-based BEs, soil, and humans. Our objectives were to (i) identify differences in the pangenomic composition of generalist and host-associated organisms, (ii) characterize genes and functions involved in BE-associated selection, and (iii) identify genomic signatures of ISS-derived strains of potential relevance for astronaut health. The pangenome of B. cereus was more expansive than that of S. aureus, which had a dominant core component. Genomic contents of both taxa significantly correlated with isolate origin, demonstrating an importance for biogeography and potential niche adaptations. ISS/BE-enriched functions were often involved in biosynthesis, catabolism, materials transport, metabolism, and stress response. Multiple origin-enriched functions also overlapped across taxa, suggesting conserved adaptive processes. We further characterized two mobile genetic elements with local neighborhood genes encoding biosynthesis and stress response functions that distinctively associated with B. cereus from the ISS. Although antibiotic resistance genes were present in ISS/BE isolates, they were also common in counterparts elsewhere. Overall, despite differences in microbial lifestyle, some functions appear common to remaining viable in the BE, and those functions are not typically associated with direct impacts on human health. IMPORTANCE The built environment contains a variety of microorganisms, some of which pose critical human health risks (e.g., hospital-acquired infection, antibiotic resistance dissemination). We uncovered a combination of complex biological functions that may play a role in bacterial survival under the presumed selective pressures in a model built environment-the International Space Station-by using an approach to compare pangenomes of bacterial strains from two clinically relevant species (B. cereus and S. aureus) isolated from both built environments and humans. Our findings suggest that the most crucial bacterial functions involved in this potential adaptive response are specific to bacterial lifestyle and do not appear to have direct impacts on human health.

Role of plasmid plasticity and mobile genetic elements in the entomopathogen Bacillus thuringiensis serovar israelensis

Bacillus thuringiensis is a well-known biopesticide that has been used for more than 80 years. This spore-forming bacterium belongs to the group of Bacillus cereus that also includes, among others, emetic and diarrheic pathotypes of B. cereus, the animal pathogen Bacillus anthracis and the psychrotolerant Bacillus weihenstephanensis. Bacillus thuringiensis is rather unique since it has adapted its lifestyle as an efficient pathogen of specific insect larvae. One of the peculiarities of B. thuringiensis strains is the extent of their extrachromosomal pool, with strains harbouring more than 10 distinct plasmid molecules. Among the numerous serovars of B. thuringiensis, 'israelensis' is certainly emblematic since its host spectrum is apparently restricted to dipteran insects like mosquitoes and black flies, vectors of human and animal diseases such as malaria, yellow fever, or river blindness. In this review, the putative role of the mobile gene pool of B. thuringiensis serovar israelensis in its pathogenicity and dedicated lifestyle is reviewed, with specific emphasis on the nature, diversity, and potential mobility of its constituents. Variations among the few related strains of B. thuringiensis serovar israelensis will also be reported and discussed in the scope of this specialised insect pathogen, whose lifestyle in the environment remains largely unknown.

RefSeq database growth influences the accuracy of k-mer-based lowest common ancestor species identification

In order to determine the role of the database in taxonomic sequence classification, we examine the influence of the database over time on k-mer-based lowest common ancestor taxonomic classification. We present three major findings: the number of new species added to the NCBI RefSeq database greatly outpaces the number of new genera; as a result, more reads are classified with newer database versions, but fewer are classified at the species level; and Bayesian-based re-estimation mitigates this effect but struggles with novel genomes. These results suggest a need for new classification approaches specially adapted for large databases.

Environmentally-driven gene content convergence and the Bacillus phylogeny

Background

Members of the Bacillus genus have been isolated from a variety of environments. However, the relationship between potential metabolism and the niche from which bacteria of this genus have been isolated has not been extensively studied. The existence of a monophyletic aquatic Bacillus group, composed of members isolated from both marine and fresh water has been proposed. Here, we present a phylogenetic/phylogenomic analysis to investigate the potential relationship between the environment from which group members have been isolated and their evolutionary origin. We also carried out hierarchical clustering based on functional content to test for potential environmental effects on the genetic content of these bacteria.

Results

The phylogenetic reconstruction showed that Bacillus strains classified as aquatic have evolutionary origins in different lineages. Although we observed the presence of a clade consisting exclusively of aquatic Bacillus, it is not comprised of the same strains previously reported. In contrast to phylogeny, clustering based on the functional categories of the encoded proteomes resulted in groups more compatible with the environments from which the organisms were isolated. This evidence suggests a detectable environmental influence on bacterial genetic content, despite their different evolutionary origins.

Conclusion

Our results suggest that aquatic Bacillus species have polyphyletic origins, but exhibit convergence at the gene content level.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1261-7) contains supplementary material, which is available to authorized users.

Genomic characterization of three novel Basilisk-like phages infecting Bacillus anthracis

Background

In the present study, we sequenced the complete genomes of three novel bacteriophages v_B-Bak1, v_B-Bak6, v_B-Bak10 previously isolated from historical anthrax burial sites in the South Caucasus country of Georgia. We report here major trends in the molecular evolution of these phages, which we designate as “Basilisk-Like-Phages” (BLPs), and illustrate patterns in their evolution, genomic plasticity and core genome architecture.

Results

Comparative whole genome sequence analysis revealed a close evolutionary relationship between our phages and two unclassified Bacillus cereus group phages, phage Basilisk, a broad host range phage (Grose JH et al., J Vir. 2014;88(20):11846-11860) and phage PBC4, a highly host-restricted phage and close relative of Basilisk (Na H. et al. FEMS Microbiol. letters. 2016;363(12)). Genome comparisons of phages v_B-Bak1, v_B-Bak6, and v_B-Bak10 revealed significant similarity in sequence, gene content, and synteny with both Basilisk and PBC4. Transmission electron microscopy (TEM) confirmed the three phages belong to the Siphoviridae family. In contrast to the broad host range of phage Basilisk and the single-strain specificity of PBC4, our three phages displayed host specificity for Bacillus anthracis. Bacillus species including Bacillus cereus, Bacillus subtilis, Bacillus anthracoides, and Bacillus megaterium were refractory to infection.

Conclusions

Data reported here provide further insight into the shared genomic architecture, host range specificity, and molecular evolution of these rare B. cereus group phages. To date, the three phages represent the only known close relatives of the Basilisk and PBC4 phages and their shared genetic attributes and unique host specificity for B. anthracis provides additional insight into candidate host range determinants.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5056-4) contains supplementary material, which is available to authorized users.

Fine-scale differentiation between Bacillus anthracis and Bacillus cereus group signatures in metagenome shotgun data

BACKGROUND

It is possible to detect bacterial species in shotgun metagenome datasets through the presence of only a few sequence reads. However, false positive results can arise, as was the case in the initial findings of a recent New York City subway metagenome project. False positives are especially likely when two closely related are present in the same sample. Bacillus anthracis, the etiologic agent of anthrax, is a high-consequence pathogen that shares >99% average nucleotide identity with Bacillus cereus group (BCerG) genomes. Our goal was to create an analysis tool that used k-mers to detect B. anthracis, incorporating information about the coverage of BCerG in the metagenome sample.

METHODS

Using public complete genome sequence datasets, we identified a set of 31-mer signatures that differentiated B. anthracis from other members of the B. cereus group (BCerG), and another set which differentiated BCerG genomes (including B. anthracis) from other Bacillus strains. We also created a set of 31-mers for detecting the lethal factor gene, the key genetic diagnostic of the presence of anthrax-causing bacteria. We created synthetic sequence datasets based on existing genomes to test the accuracy of a k-mer based detection model.

RESULTS

We found 239,503 B. anthracis-specific 31-mers (the Ba31 set), 10,183 BCerG 31-mers (the BCerG31 set), and 2,617 lethal factor k-mers (the lef31 set). We showed that false positive B. anthracis k-mers-which arise from random sequencing errors-are observable at high genome coverages of B. cereus. We also showed that there is a "gray zone" below 0.184× coverage of the B. anthracis genome sequence, in which we cannot expect with high probability to identify lethal factor k-mers. We created a linear regression model to differentiate the presence of B. anthracis-like chromosomes from sequencing errors given the BCerG background coverage. We showed that while shotgun datasets from the New York City subway metagenome project had no matches to lef31 k-mers and hence were negative for B. anthracis, some samples showed evidence of strains very closely related to the pathogen.

DISCUSSION

This work shows how extensive libraries of complete genomes can be used to create organism-specific signatures to help interpret metagenomes. We contrast "specialist" approaches to metagenome analysis such as this work to "generalist" software that seeks to classify all organisms present in the sample and note the more general utility of a k-mer filter approach when taxonomic boundaries lack clarity or high levels of precision are required.

Genomic Characterization and Copy Number Variation of Bacillus anthracis Plasmids pXO1 and pXO2 in a Historical Collection of 412 Strains

Bacillus anthracis microorganisms are of historical and epidemiological importance and are among the most homogenous bacterial groups known, even though the B. anthracis genome is rich in mobile elements. Mobile elements can trigger the diversification of lineages; therefore, characterizing the extent of genomic variation in a large collection of strains is critical for a complete understanding of the diversity and evolution of the species. Here, we sequenced a large collection of B. anthracis strains (>400) that were recovered from human, animal, and environmental sources around the world. Our results confirmed the remarkable stability of gene content and synteny of the anthrax plasmids and revealed no signal of plasmid exchange between B. anthracis and pathogenic B. cereus isolates but rather predominantly vertical descent. These findings advance our understanding of the biology and pathogenomic evolution of B. anthracis and its plasmids.

Bacillus anthracis plasmids pXO1 and pXO2 carry the main virulence factors responsible for anthrax. However, the extent of copy number variation within the species and how the plasmids are related to pXO1/pXO2-like plasmids in other species of the Bacillus cereussensu lato group remain unclear. To gain new insights into these issues, we sequenced 412 B. anthracis strains representing the total phylogenetic and ecological diversity of the species. Our results revealed that B. anthracis genomes carried, on average, 3.86 and 2.29 copies of pXO1 and pXO2, respectively, and also revealed a positive linear correlation between the copy numbers of pXO1 and pXO2. No correlation between the plasmid copy number and the phylogenetic relatedness of the strains was observed. However, genomes of strains isolated from animal tissues generally maintained a higher plasmid copy number than genomes of strains from environmental sources (P < 0.05 [Welch two-sample t test]). Comparisons against B. cereus genomes carrying complete or partial pXO1-like and pXO2-like plasmids showed that the plasmid-based phylogeny recapitulated that of the main chromosome, indicating limited plasmid horizontal transfer between or within these species. Comparisons of gene content revealed a closed pXO1 and pXO2 pangenome; e.g., plasmids encode <8 unique genes, on average, and a single large fragment deletion of pXO1 in one B. anthracis strain (2000031682) was detected. Collectively, our results provide a more complete view of the genomic diversity of B. anthracis plasmids, their copy number variation, and the virulence potential of other Bacillus species carrying pXO1/pXO2-like plasmids.

IMPORTANCEBacillus anthracis microorganisms are of historical and epidemiological importance and are among the most homogenous bacterial groups known, even though the B. anthracis genome is rich in mobile elements. Mobile elements can trigger the diversification of lineages; therefore, characterizing the extent of genomic variation in a large collection of strains is critical for a complete understanding of the diversity and evolution of the species. Here, we sequenced a large collection of B. anthracis strains (>400) that were recovered from human, animal, and environmental sources around the world. Our results confirmed the remarkable stability of gene content and synteny of the anthrax plasmids and revealed no signal of plasmid exchange between B. anthracis and pathogenic B. cereus isolates but rather predominantly vertical descent. These findings advance our understanding of the biology and pathogenomic evolution of B. anthracis and its plasmids.

Assembly and Function of the Bacillus anthracis S-Layer

Bacillus anthracis, the anthrax agent, is a member of the Bacillus cereus sensu lato group, which includes invasive pathogens of mammals or insects as well as nonpathogenic environmental strains. The genes for anthrax pathogenesis are located on two large virulence plasmids. Similar virulence plasmids have been acquired by other B. cereus strains and enable the pathogenesis of anthrax-like diseases. Among the virulence factors of B. anthracis is the S-layer-associated protein BslA, which endows bacilli with invasive attributes for mammalian hosts. BslA surface display and function are dependent on the bacterial S-layer, whose constituents assemble by binding to the secondary cell wall polysaccharide (SCWP) via S-layer homology (SLH) domains. B. anthracis and other pathogenic B. cereus isolates harbor genes for the secretion of S-layer proteins, for S-layer assembly, and for synthesis of the SCWP. We review here recent insights into the assembly and function of the S-layer and the SCWP.

Spores and soil from six sides: interdisciplinarity and the environmental biology of anthrax (Bacillus anthracis)

Environmentally transmitted diseases are comparatively poorly understood and managed, and their ecology is particularly understudied. Here we identify challenges of studying environmental transmission and persistence with a six-sided interdisciplinary review of the biology of anthrax (Bacillus anthracis). Anthrax is a zoonotic disease capable of maintaining infectious spore banks in soil for decades (or even potentially centuries), and the mechanisms of its environmental persistence have been the topic of significant research and controversy. Where anthrax is endemic, it plays an important ecological role, shaping the dynamics of entire herbivore communities. The complex eco-epidemiology of anthrax, and the mysterious biology of Bacillus anthracis during its environmental stage, have necessitated an interdisciplinary approach to pathogen research. Here, we illustrate different disciplinary perspectives through key advances made by researchers working in Etosha National Park, a long-term ecological research site in Namibia that has exemplified the complexities of the enzootic process of anthrax over decades of surveillance. In Etosha, the role of scavengers and alternative routes (waterborne transmission and flies) has proved unimportant relative to the long-term persistence of anthrax spores in soil and their infection of herbivore hosts. Carcass deposition facilitates green-ups of vegetation to attract herbivores, potentially facilitated by the role of anthrax spores in the rhizosphere. The underlying seasonal pattern of vegetation, and herbivores' immune and behavioural responses to anthrax risk, interact to produce regular 'anthrax seasons' that appear to be a stable feature of the Etosha ecosystem. Through the lens of microbiologists, geneticists, immunologists, ecologists, epidemiologists, and clinicians, we discuss how anthrax dynamics are shaped at the smallest scale by population genetics and interactions within the bacterial communities up to the broadest scales of ecosystem structure. We illustrate the benefits and challenges of this interdisciplinary approach to disease ecology, and suggest ways anthrax might offer insights into the biology of other important pathogens. Bacillus anthracis, and the more recently emerged Bacillus cereus biovar anthracis, share key features with other environmentally transmitted pathogens, including several zoonoses and panzootics of special interest for global health and conservation efforts. Understanding the dynamics of anthrax, and developing interdisciplinary research programs that explore environmental persistence, is a critical step forward for understanding these emerging threats.

Adaptive Strategies of Bacillus thuringiensis Isolated from Acid Mine Drainage Site in Sabah, Malaysia

The adaptive process in bacteria is driven by specific genetic elements which regulate phenotypic characteristics such as tolerance to high metal ion concentrations and the secretion of protective biofilms. Extreme environments such as those associated with heavy metal pollution and extremes of acidity offer opportunities to study the adaptive mechanisms of microorganisms. This study focused on the genome analysis of Bacillus thuringiensis (Bt MCMY1), a gram positive rod shaped bacterium isolated from an acid mine drainage site in Sabah, Malaysia by using a combination of Single Molecule Real Time DNA Sequencing, Scanning Electron Microscopy (SEM) and Fourier Transform Infrared Spectroscopy (FTIR). The genome size of Bt MCMY1 was determined to be 5,458,152 bases which was encoded on a single chromosome. Analysis of the genome revealed genes associated with resistance to Copper, Mercury, Arsenic, Cobalt, Zinc, Cadmium and Aluminum. Evidence from SEM and FTIR indicated that the bacterial colonies form distinct films which bear the signature of polyhydroxyalkanoates (PHA) and this finding was supported by the genome data indicating the presence of a genetic pathway associated with the biosynthesis of PHAs. This is the first report of a Bacillus sp. isolated from an acid mine drainage site in Sabah, Malaysia and the genome sequence will provide insights into the manner in which B. thuringiensis adapts to acid mine drainage.

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.

Intraclade Variability in Toxin Production and Cytotoxicity of Bacillus cereus Group Type Strains and Dairy-Associated Isolates

While some species in the Bacillus cereus group are well-characterized human pathogens (e.g., B. anthracis and B. cereus sensu stricto), the pathogenicity of other species (e.g., B. pseudomycoides) either has not been characterized or is presently not well understood. To provide an updated characterization of the pathogenic potential of species in the B. cereus group, we classified a set of 52 isolates, including 8 type strains and 44 isolates from dairy-associated sources, into 7 phylogenetic clades and characterized them for (i) the presence of toxin genes, (ii) phenotypic characteristics used for identification, and (iii) cytotoxicity to human epithelial cells. Overall, we found that B. cereus toxin genes are broadly distributed but are not consistently present within individual species and/or clades. After growth at 37°C, isolates within a clade did not typically show a consistent cytotoxicity phenotype, except for isolates in clade VI (B. weihenstephanensis/B. mycoides), where none of the isolates were cytotoxic, and isolates in clade I (B. pseudomycoides), which consistently displayed cytotoxic activity. Importantly, our study highlights that B. pseudomycoides is cytotoxic toward human cells. Our results indicate that the detection of toxin genes does not provide a reliable approach to predict the pathogenic potential of B. cereus group isolates, as the presence of toxin genes is not always consistent with cytotoxicity phenotype. Overall, our results suggest that isolates from multiple B. cereus group clades have the potential to cause foodborne illness, although cytotoxicity is not always consistently found among isolates within each clade.IMPORTANCE Despite the importance of the Bacillus cereus group as a foodborne pathogen, characterizations of the pathogenic potential of all B. cereus group species were lacking. We show here that B. pseudomycoides (clade I), which has been considered a harmless environmental microorganism, produces toxins and exhibits a phenotype consistent with the production of pore-forming toxins. Furthermore, B. mycoides/B. weihenstephanensis isolates (clade VI) did not show cytotoxicity when grown at 37°C, despite carrying multiple toxin genes. Overall, we show that the current standard methods to characterize B. cereus group isolates and to detect the presence of toxin genes are not reliable indicators of species, phylogenetic clades, or an isolate's cytotoxic capacity, suggesting that novel methods are still needed for differentiating pathogenic from nonpathogenic species within the B. cereus group. Our results also contribute data that are necessary to facilitate risk assessments and a better understanding as to which B. cereus group species are likely to cause foodborne illness.

Structures of the Peptidoglycan N-Acetylglucosamine Deacetylase Bc1974 and Its Complexes with Zinc Metalloenzyme Inhibitors

The cell wall peptidoglycan is recognized as a primary target of the innate immune system, and usually its disintegration results in bacterial lysis. Bacillus cereus, a close relative of the highly virulent Bacillus anthracis, contains 10 polysaccharide deacetylases. Among these, the peptidoglycan N-acetylglucosamine deacetylase Bc1974 is the highest homologue to the Bacillus anthracis Ba1977 that is required for full virulence and is involved in resistance to the host's lysozyme. These metalloenzymes belong to the carbohydrate esterase family 4 (CE4) and are attractive targets for the development of new anti-infective agents. Herein we report the first X-ray crystal structures of the NodB domain of Bc1974, the conserved catalytic core of CE4s, in the unliganded form and in complex with four known metalloenzyme inhibitors and two amino acid hydroxamates that target the active site metal. These structures revealed the presence of two conformational states of a catalytic loop known as motif-4 (MT4), which were not observed previously for peptidoglycan deacetylases, but were recently shown in the structure of a Vibrio clolerae chitin deacetylase. By employing molecular docking of a substrate model, we describe a catalytic mechanism that probably involves initial binding of the substrate in a receptive, more open state of MT4 and optimal catalytic activity in the closed state of MT4, consistent with the previous observations. The ligand-bound structures presented here, in addition to the five Bc1974 inhibitors identified, provide a valuable basis for the design of antibacterial agents that target the peptidoglycan deacetylase Ba1977.

A minireplicon of plasmid pBMB26 represents a new typical replicon in the megaplasmids of Bacillus cereus group

A new minireplicon (rep26 minireplicon) from pBMB26, the 188 kb indigenous plasmid related to spore-crystal association (SCA) phenotype in Bacillus thuringiensis strain YBT-020, was characterized. A 12 kb EcoRI fragment, which encoded 10 putative open reading frames (ORFs), was capable of supporting replication when cloned in a replication probe vector. Deletion and frame shift mutation analysis showed that a 4.1 kb region encompassing two putative ORFs (orf21 and orf22) was essential for the plasmid replication in B. thuringiensis. Gene orf21 encoding a 49.8 kDa protein (named Rep26) with a helix-turn-helix motif showed no homology with known replication proteins and gene orf22 encoding a protein of 82.6 kDa showed homology to bacterial PcrA helicase. The replication origin of rep26 minireplicon was proved to be located in the coding region of orf21. Plasmid stability experiments indicated that the recombinant plasmid containing rep26 minireplicon has excellent segregational stability. BLASTP analysis revealed that amino acid sequences of ORF21 and ORF22 were well conserved among Bacillus cereus group strains. The rep26 minireplicon was widely distributed and could be defined as a new typical replicon in the megaplasmids of B. cereus group.