Core genome multilocus sequence typing scheme for Bacillus cereus group bacteria

Antimicrobial and Phylogenomic Characterization of Bacillus cereus Group Strains Isolated from Different Food Sources in Italy

Background:Bacillus cereus is a widespread environmental Gram-positive bacterium which is especially common in soil and dust. It produces two types of toxins that cause vomiting and diarrhea. At present, foodborne outbreaks due to Bacillus cereus group bacteria (especially Bacillus cereus sensu stricto) are rising, representing a serious problem in the agri-food supply chain. Methods: In this work, we analyzed 118 strains belonging to the Bacillus cereus group, isolated from several food sources, for which in vitro and in silico antibiotic resistance assessments were performed. Results: Many strains showed intermediate susceptibility to clindamycin, erythromycin, and tetracycline, suggesting an evolving acquisition of resistance against these antibiotics. Moreover, one strain showed intermediate resistance to meropenem, an antibiotic currently used to treat infections caused by Bacillus cereus. In addition to the phenotypic antimicrobial resistance profile, all strains were screened for the presence/absence of antimicrobial genes via whole-genome sequencing. There was inconsistency between the in vitro and in silico analyses, such as in the case of vancomycin, for which different isolates harbored resistance genes but, phenotypically, the same strains were sensitive. Conclusions: This would suggest that antibiotic resistance is a complex phenomenon due to a variety of genetic, epigenetic, and biochemical mechanisms.

Development and implementation of a core genome multilocus sequence typing scheme for Haemophilus influenzae

Haemophilus influenzae is part of the human nasopharyngeal microbiota and a pathogen causing invasive disease. The extensive genetic diversity observed in H. influenzae necessitates discriminatory analytical approaches to evaluate its population structure. This study developed a core genome multilocus sequence typing (cgMLST) scheme for H. influenzae using pangenome analysis tools and validated the cgMLST scheme using datasets consisting of complete reference genomes (N = 14) and high-quality draft H. influenzae genomes (N = 2297). The draft genome dataset was divided into a development dataset (N = 921) and a validation dataset (N = 1376). The development dataset was used to identify potential core genes, and the validation dataset was used to refine the final core gene list to ensure the reliability of the proposed cgMLST scheme. Functional classifications were made for all the resulting core genes. Phylogenetic analyses were performed using both allelic profiles and nucleotide sequence alignments of the core genome to test congruence, as assessed by Spearman's correlation and ordinary least square linear regression tests. Preliminary analyses using the development dataset identified 1067 core genes, which were refined to 1037 with the validation dataset. More than 70% of core genes were predicted to encode proteins essential for metabolism or genetic information processing. Phylogenetic and statistical analyses indicated that the core genome allelic profile accurately represented phylogenetic relatedness among the isolates (R 2 = 0.945). We used this cgMLST scheme to define a high-resolution population structure for H. influenzae, which enhances the genomic analysis of this clinically relevant human pathogen.

Antimicrobial Resistance Profile, Whole-Genome Sequencing and Core Genome Multilocus Sequence Typing of B. anthracis Isolates in Croatia from 2001 to 2022

Bacillus anthracis, the causative agent of anthrax disease, is a worldwide threat to livestock, wildlife and public health. It is also considered one of the most important pathogens of bioterrorism. Rapid and reliable diagnosis and administration of antimicrobials are essential for effective anthrax treatment. In this study, we determined the in vitro susceptibilities of 40 isolates of B. anthracis isolated in Croatia over the recent two decades to 18 antimicrobials. Whole-genome sequencing was performed, and bioinformatics tools were used to determine virulence factors and antimicrobial resistance genes. Core genome-based multilocus sequence typing was used for isolate comparison and phylogenetic analysis. All isolates were susceptible to all antimicrobials recommended for post-exposure prophylaxis or anthrax therapy. Susceptibility was found to all other tested antimicrobials that are an alternative for primary therapy. We found two beta-lactamase genes, but their expression is not sufficient to confer resistance. In all isolates used in this study, we found 21 virulence genes, 8 of which are responsible for toxin and capsule production. As far as phylogenetic analysis is concerned, the B. anthracis isolates from Croatia are categorised into two clades. The first is clade A, subclade Trans Eurasia, and the other is clade B, subclade B2.

Cereulide production capacities and genetic properties of 31 emetic Bacillus cereus group strains

The highly potent toxin cereulide is a frequent cause of foodborne intoxications. This extremely resistant toxin is produced by Bacillus cereus group strains carrying the plasmid encoded cesHPTABCD gene cluster. It is known that the capacities to produce cereulide vary greatly between different strains but the genetic background of these variations is not clear. In this study, cereulide production capacities were associated with genetic characteristics. For this, cereulide levels in cultures of 31 strains were determined after incubation in tryptic soy broth for 24 h at 24 °C, 30 °C and 37 °C. Whole genome sequencing based data were used for an in-depth characterization of gene sequences related to cereulide production. The taxonomy, population structure and phylogenetic relationships of the strains were evaluated based on average nucleotide identity, multi-locus sequence typing (MLST), core genome MLST and single nucleotide polymorphism analyses. Despite a limited strain number, the approach of a genome wide association study (GWAS) was tested to link genetic variation with cereulide quantities. Our study confirms strain-dependent differences in cereulide production. For most strains, these differences were not explainable by sequence variations in the cesHPTABCD gene cluster or the regulatory genes abrB, spo0A, codY and pagRBc. Likewise, the population structure and phylogeny of the tested strains did not comprehensively reflect the cereulide production capacities. GWAS yielded first hints for associated proteins, while their possible effect on cereulide synthesis remains to be further investigated.

Prevalence and genomic characterization of the Bacillus cereus group strains contamination in food products in Southern China

The Bacillus cereus group, as one of the important opportunistic foodborne pathogens, is considered a risk to public health due to foodborne diseases and an important cause of economic losses to food industries. This study aimed to gain essential information on the prevalence, phenotype, and genotype of B. cereus group strains isolated from various food products in China. A total of 890 strains of B. cereus group bacteria from 1181 food samples from 2020 to 2023 were identified using the standardized detection method. These strains were found to be prevalent in various food types, with the highest contamination rates observed in cereal flour (55.8 %) and wheat/rice noodles (45.7 %). The tested strains exhibited high resistance rates against penicillin (98.5 %) and ampicillin (98.9 %). Strains isolated from cereal flour had the highest rate of meropenem resistance (7.8 %), while strains from sausages were most resistant to vancomycin (16.8 %). A total of 234 out of the 891 B. cereus group strains were randomly selected for WGS analysis, 18.4 % of which displayed multidrug resistance. The species identification by WGS analysis revealed the presence of 10 distinct species within the B. cereus group, with B. cereus species being the most prevalent. The highest level of species diversity was observed in sausages. Notably, B. anthracis strains lacking the anthrax toxin genes were detected in flour-based food products and sausages. A total of 20 antibiotic resistance genes have been identified, with β-lactam resistance genes (bla1, bla2, BcI, BcII, and blaTEM-116) being the most common. The B. tropicus strains exhibit the highest average number of virulence genes (23.4). The diarrheal virulence genes nheABC, hblACD, and cytK were found in numerous strains. Only 4 of the 234 (1.7 %) sequenced strains contain the ces gene cluster linked to emetic symptoms. These data offer valuable insights for public health policymakers on addressing foodborne B. cereus group infections and ensuring food safety.

Co-inoculation of the endophytes Bacillus thuringiensis CAPE95 and Paenibacillus polymyxa CAPE238 promotes Tropaeolum majus L. growth and enhances its root bacterial diversity

Tropaeolum majus L. is a versatile edible plant that is widely explored due to its medicinal properties and as a key element in intercropping systems. Its growth could be improved by the use of biofertilizers that can enhance nutrient uptake by the plant or provide tolerance to different abiotic and biotic stresses. In a previous study, 101 endophytes isolated from T. majus roots showed more than three plant growth-promoting (PGP) features in vitro, such as phosphate mineralization/solubilization, production of siderophores, antimicrobial substances and indole-related compounds, and presence of the nifH gene. To provide sustainable alternatives for biofertilization, the genomes of two promising endophytes-CAPE95 and CAPE238-were sequenced to uncover metabolic pathways related to biofertilization. Greenhouse experiments were conducted with 216 seeds and 60 seedlings, half co-inoculated with the endophytes (treatment) and half inoculated with 1X PBS (control), and the impact of the co-inoculation on the plant's bacteriome was accessed through 16S rRNA gene metabarcoding. The strains CAPE95 and CAPE238 were taxonomically assigned as Bacillus thuringiensis and Paenibacillus polymyxa, respectively. Metabolic pathways related to the enhancement of nutrient availability (nitrogen fixation, sulfate-sulfur assimilation), biosynthesis of phytohormones (indole-3-acetic acid precursors) and antimicrobial substances (bacilysin, paenibacillin) were found in their genomes. The in vivo experiments showed that treated seeds exhibited faster germination, with a 20.3% higher germination index than the control on the eleventh day of the experiment. Additionally, treated seedlings showed significantly higher plant height and leaf diameters (p < 0.05). The bacterial community of the treated plants was significantly different from that of the control plants (p < 0.001) and showed a higher richness and diversity of species (Chao and Shannon indexes, p < 0.001). A higher relative abundance of potential synergistic PGP bacteria was also shown in the bacteriome of the treated plants, such as Lysinibacillus and Geobacter. For the first time, co-inoculation of B. thuringiensis and P. polymyxa was shown to have great potential for application as a biofertilizer to T. majus plants. The bacterial consortium used here could also be explored in other plant species in the future.

Investigating environmental transmission to resolve a Bacillus cereus group outbreak in a neonatal intensive care unit using core genome multilocus sequence typing

BACKGROUND

We analyzed an outbreak of Bacillus cereus group (Bcg) at a single-center neonatal intensive care unit level IV by conducting comprehensive sampling of both patients and the environment.

METHODS

Between 06/2020 and 10/2021, all Bcg isolates identified by both regular colonization screening and additional sampling of the environment were subjected to whole-genome sequencing, followed by in vitro extraction of MLST ST, resistance genes and virulence factors. Using publicly available genome sequences, we defined an ad hoc core genome multilocus sequence typing (cgMLST) scheme comprising 2759 target genes for Bcg typing, which we applied to the detected isolates. We have compared the results with a stable cgMLST that was published in the meantime and completed the investigation with a SNP analysis.

RESULTS

We analyzed 28 Bcg isolates from patient and environmental samples using MLST and cgMLST. This revealed multiple sequence types, with ST127 being the most common (n = 13). Both cgMLST schemes grouped ten of the 13 ST127 isolates into a cluster, including two invasive isolates from two different patients and several environmental samples. SNP analysis postulated a screen from a ventilation machine as a possible reservoir.

CONCLUSION

In sensitive settings such as neonatal intensive care units, considering the environment in outbreak analyses is crucial, especially when investigating potential transmission routes through shared devices. When dealing with widespread bacteria such as Bcg, high-resolution typing techniques are necessary. In this study, we successfully resolved an outbreak of Bcg infections using a custom cgMLST scheme combined with a SNP analysis.

The toxicological spectrum of the Bacillus cereus toxin cereulide points towards niche-specific specialisation

Most microbes share their environmental niches with very different forms of life thereby engaging in specialised relationships to enable their persistence. The bacterium Bacillus cereus occurs ubiquitously in the environment with certain strain backgrounds causing foodborne and opportunistic infections in humans. The emetic lineage of B. cereus is capable of producing the toxin cereulide, which evokes emetic illnesses. Although food products favouring the accumulation of cereulide are known, the ecological role of cereulide and the environmental niche of emetic B. cereus remain elusive. To better understand the ecology of cereulide-producing B. cereus, we systematically assayed the toxicological spectrum of cereulide on a variety of organisms belonging to different kingdoms. As cereulide is a potassium ionophore, we further tested the effect of environmental potassium levels on the action of cereulide. We found that adverse effects of cereulide exposure are species-specific, which can be exacerbated with increased environmental potassium. Additionally, we demonstrate that cereulide is produced within an insect cadaver indicating its potential ecological function for a saprophytic lifestyle. Collectively, distinct cereulide susceptibilities of other organisms may reflect its role in enabling competitive niche specialization of emetic B. cereus.