Genomic insights into the taxonomic status of the Bacillus cereus group

Whole Genome Sequencing Reveals Antimicrobial Resistance and Virulence Genes of Both Pathogenic and Non-Pathogenic B. cereus Group Isolates from Foodstuffs in Thailand

Members of the Bacillus cereus group are spore-forming Gram-positive bacilli that are commonly associated with diarrheal or emetic food poisoning. They are widespread in nature and frequently present in both raw and processed food products. Here, we genetically characterized 24 B. cereus group isolates from foodstuffs. Whole-genome sequencing (WGS) revealed that most of the isolates were closely related to B. cereus sensu stricto (12 isolates), followed by B. pacificus (5 isolates), B. paranthracis (5 isolates), B. tropicus (1 isolate), and "B. bingmayongensis" (1 isolate). The most detected virulence genes were BAS_RS06430, followed by bacillibactin biosynthesis genes (dhbA, dhbB, dhbC, dhbE, and dhbF), genes encoding the three-component non-hemolytic enterotoxin (nheA, nheB, and nheC), a gene encoding an iron-regulated leucine-rich surface protein (ilsA), and a gene encoding a metalloprotease (inhA). Various biofilm-associated genes were found, with high prevalences of tasA and sipW genes (matrix protein-encoding genes); purA, purC, and purL genes (eDNA synthesis genes); lytR and ugd genes (matrix polysaccharide synthesis genes); and abrB, codY, nprR, plcR, sinR, and spo0A genes (biofilm transcription regulator genes). Genes related to fosfomycin and beta-lactam resistance were identified in most of the isolates. We therefore demonstrated that WGS analysis represents a useful tool for rapidly identifying and characterizing B. cereus group strains. Determining the genetic epidemiology, the presence of virulence and antimicrobial resistance genes, and the pathogenic potential of each strain is crucial for improving the risk assessment of foodborne B. cereus group strains.

Microbiological evaluation of the indigenous fermented condiment okpeye available at various retail markets in the south-eastern region of Nigeria

In Africa, indigenous fermented condiments contribute to food security as a low-cost source of protein. Okpeye is an indigenous fermented condiment produced from Prosopis africana seeds. The reliance on spontaneous fermentation processes and unhygienic practices during production often results in the contamination of the final product with microbial hazards. A microbiological evaluation of 18 commercial samples of okpeye purchased from six markets in two cities in southeastern Nigeria was conducted. Fifty-nine (59) bacteria were isolated and identified at the species level by phenotyping and sequencing the 16S rRNA, gyrB and rpoB genes. Bacillus (47.4 %) and Staphylococcus (42.3 %) were the predominant bacterial genera in okpeye. Overall, B. amyloliquefaciens and S. simulans were the most frequently occurring bacteria and were present in all samples. In addition, B. cereus was isolated in samples obtained from all markets. Other bacterial species included B. velezensis, Oceanobacillus caeni, S. cohnii, Escherichia fergusonni and Vagacoccus lutrae. The B. cereus isolates (10) were screened for the presence of 8 enterotoxin genes (hblA, hblC, hblD, nheA, nheB, nheC, cytK, entFM) and one emetic gene (cesB). The non-haemolytic enterotoxin (nheABC) and haemolytic enterotoxin (hblABD) complexes were present in 70 % and 50 % of B. cereus respectively. The positive rate of cytK and entFM genes was 70 %, while the cesB gene was 30 %. Antibiotic susceptibility assessment showed that most of the isolates were susceptible to gentamicin, tetracycline, streptomycin, and erythromycin but resistant to ciprofloxacin and vancomycin. These findings highlight the need for further controls to reduce contamination with potential pathogenic bacteria in indigenous fermented condiments such as okpeye. There is also a need to educate producers regarding hygienic practices to safeguard public health and food security.

Clostridium strain FAM25158, a unique endospore-forming bacterium related to Clostridium tyrobutyricum and isolated from Emmental cheese shows low tolerance to salt

The genus Clostridium is a large and diverse group of species that can cause food spoilage, including late blowing defect (LBD) in cheese. In this study, we investigated the taxonomic status of strain FAM25158 isolated from Emmental cheese with LBD using a polyphasic taxonomic and comparative genomic approach. A 16S rRNA gene sequence phylogeny suggested affiliation to the Clostridium sensu stricto cluster, with Clostridium tyrobutyricum DSM 2637T being the closest related type strain (99.16% sequence similarity). Average Nucleotide Identity (ANI) analysis revealed that strain FAM25158 is at the species threshold with C. tyrobutyricum, with ANI values ranging from 94.70 to 95.26%, while the digital DNA-DNA hybridization values were below the recommended threshold, suggesting that FAM25158 is significantly different from C. tyrobutyricum at the genomic level. Moreover, comparative genomic analysis between FAM25158 and its four closest C. tyrobutyricum relatives revealed a diversity of metabolic pathways, with FAM25158 differing from other C. tyrobutyricum strains by the presence of genes such as scrA, srcB, and scrK, responsible for sucrose utilization, and the absence of many important functional genes associated with cold and osmolality adaptation, which was further supported by phenotypic analyses. Surprisingly, strain FAM25158 exhibited unique physiologic traits, such as an optimal growth temperature of 30°C, in contrast to its closest relatives, C. tyrobutyricum species with an optimal growth temperature of 37°C. Additionally, the growth of FAM25158 was inhibited at NaCl concentrations higher than 0.5%, a remarkable observation considering its origin from cheese. While the results of this study provide novel information on the genetic content of strain FAM25158, the relationship between its genetic content and the observed phenotype remains a topic requiring further investigation.

De novo assembly and comparative genome analysis for polyhydroxyalkanoates-producing Bacillus sp. BNPI-92 strain

BACKGROUND

Certain Bacillus species play a vital role in polyhydroxyalkanoate (PHA) production. However, most of these isolates did not properly identify to species level when scientifically had been reported.

RESULTS

From NGS analysis, 5719 genes were predicted in the de novo genome assembly. Based on genome annotation using RAST server, 5,527,513 bp sequences were predicted with 5679 bp number of protein-coding sequence. Its genome sequence contains 35.1% and 156 GC content and contigs, respectively. In RAST server analysis, subsystem (43%) and non-subsystem coverage (57%) were generated. Ortho Venn comparative genome analysis indicated that Bacillus sp. BNPI-92 shared 2930 gene cluster (core gene) with B. cereus ATCC 14579 T (AE016877), B. paranthracis Mn5T (MACE01000012), B. thuringiensis ATCC 10792 T (ACNF01000156), and B. antrics Amen T (AE016879) strains. For our strain, the maximum gene cluster (190) was shared with B. cereus ATCC 14579 T (AE016877). For Ortho Venn pair wise analysis, the maximum overlapping gene clusters thresholds have been detected between Bacillus s p.BNPI-92 and Ba. cereus ATCC 14579 T (5414). Average nucleotide identity (ANI) such as OriginalANI and OrthoANI, in silicon digital DND-DNA hybridization (isDDH), Type (Strain) Genome Server (TYGS), and Genome-Genome Distance Calculator (GGDC) were more essentially related Bacillus sp. BNPI-92 with B. cereus ATCC 14579 T strain. Therefore, based on the combination of RAST annotation, OrthoVenn server, ANI and isDDH result Bacillus sp.BNPI-92 strain was strongly confirmed to be a B. cereus type strain. It was designated as B. cereus BNPI-92 strain. In B. cereus BNPI-92 strain whole genome sequence, PHA biosynthesis encoding genes such as phaP, phaQ, phaR (PHA synthesis repressor phaR gene sequence), phaB/phbB, and phaC were predicted on the same operon. These gene clusters were designated as phaPQRBC. However, phaA was located on other operons.

CONCLUSIONS

This newly obtained isolate was found to be new a strain based on comparative genomic analysis and it was also observed as a potential candidate for PHA biosynthesis.

Get to Know Your Neighbors: Characterization of Close Bacillus anthracis Isolates and Toxin Profile Diversity in the Bacillus cereus Group

Unexpected atypical isolates of Bacillus cereus s.l. occasionally challenge conventional microbiology and even the most advanced techniques for anthrax detection. For anticipating and gaining trust, 65 isolates of Bacillus cereus s.l. of diverse origin were sequenced and characterized. The BTyper3 tool was used for assignation to genomospecies B. mosaicus (34), B. cereus s.s (29) and B. toyonensis (2), as well as virulence factors and toxin profiling. None of them carried any capsule or anthrax-toxin genes. All harbored the non-hemolytic toxin nheABC and sphygomyelinase spH genes, whereas 41 (63%), 30 (46%), 11 (17%) and 6 (9%) isolates harbored cytK-2, hblABCD, cesABCD and at least one insecticidal toxin gene, respectively. Matrix-assisted laser desorption ionization-time of flight mass spectrometry confirmed the production of cereulide (ces genes). Phylogeny inferred from single-nucleotide polymorphisms positioned isolates relative to the B. anthracis lineage. One isolate (BC38B) was of particular interest as it appeared to be the closest B. anthracis neighbor described so far. It harbored a large plasmid similar to other previously described B. cereus s.l. megaplasmids and at a lower extent to pXO1. Whereas bacterial collection is enriched, these high-quality public genetic data offer additional knowledge for better risk assessment using future NGS-based technologies of detection.

Plant-Associated Representatives of the Bacillus cereus Group Are a Rich Source of Antimicrobial Compounds

Seventeen bacterial strains able to suppress plant pathogens have been isolated from healthy Vietnamese crop plants and taxonomically assigned as members of the Bacillus cereus group. In order to prove their potential as biocontrol agents, we perform a comprehensive analysis that included the whole-genome sequencing of selected strains and the mining for genes and gene clusters involved in the synthesis of endo- and exotoxins and secondary metabolites, such as antimicrobial peptides (AMPs). Kurstakin, thumolycin, and other AMPs were detected and characterized by different mass spectrometric methods, such as MALDI-TOF-MS and LIFT-MALDI-TOF/TOF fragment analysis. Based on their whole-genome sequences, the plant-associated isolates were assigned to the following species and subspecies: B. cereus subsp. cereus (6), B. cereus subsp. bombysepticus (5), Bacillus tropicus (2), and Bacillus pacificus. These three isolates represent novel genomospecies. Genes encoding entomopathogenic crystal and vegetative proteins were detected in B. cereus subsp. bombysepticus TK1. The in vitro assays revealed that many plant-associated isolates enhanced plant growth and suppressed plant pathogens. Our findings indicate that the plant-associated representatives of the B. cereus group are a rich source of putative antimicrobial compounds with potential in sustainable agriculture. However, the presence of virulence genes might restrict their application as biologicals in agriculture.

Bacillus dicomae sp. nov., a new member of the Bacillus cereus group isolated from medicinal plant Dicoma anomala

A Gram-stain-positive, endospore-forming endophytic bacterial strain designated MHSD28T was isolated from surface-sterilized leaves of Dicoma anomala collected from Eisleben, Botlokwa, Limpopo Province, South Africa. The phenotypic and phylogenetic characteristics of strain MHSD28T were consistent with those of members within the Bacillus cereus group. Comparative analysis between this strain and its relatives confirmed that it belongs to this group and forms a monophyletic branch. The digital DNA-DNA hybridization values between strain MHSD28T and its relatives were lower than the 70 % threshold for species delineation. To further determine its phylogenetic position, multi-locus sequence analysis (MLSA) based on five concatenated housekeeping gene (gyrB, atpD, DnaK, rpoB and rpoD) sequences, phenotypic analysis, matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) biotyper identification, fatty acid and polar lipid profile analyses were carried out. Phenotypic characterization, MLSA, whole genome sequence based analyses and MALDI-TOF results placed strain MHSD28T within the B. cereus group. The major fatty acids were iso-C15 : 0 and summed feature 3 and the main polar lipids were diphosphatidylglycerol, phosphatidylglycerol and phosphatidylethanolamine. The respiratory quinone was menaquinone-7. The cell-wall peptidoglycan structure included meso-diaminopimelic acid. Considering the above results, strain MHSD28T represents a novel species of the B. cereus group, for which the name Bacillus dicomae sp. nov. is proposed. The type strain is MHSD28T (=BD 2262T=LMG 32287T=CECT 30671T).

Elucidation of the Bovine Intramammary Bacteriome and Resistome from healthy cows of Swiss dairy farms in the Canton Tessin

Healthy, untreated cows of nine dairy herds from the Swiss Canton Tessin were analyzed three times within one year to identify the most abundant species of the intramammary bacteriome. Aseptically collected milk samples were cultured and bacteria identified using MALDI-TOF. Of 256 cows analyzed, 96% were bacteriologically positive and 80% of the 1,024 quarters were positive for at least one bacterial species. 84.5% of the quarters were healthy with somatic cell counts (SCC) < 200,000 cells/mL, whereas 15.5% of the quarters showed a subclinical mastitis (SCC ≥ 200,000 cells/mL). We could assign 1,288 isolates to 104 different bacterial species including 23 predominant species. Non-aureus staphylococci and mammaliicocci (NASM) were most prevalent (14 different species; 73.5% quarters). Staphylococcus xylosus and Mammaliicoccus sciuri accounted for 74.7% of all NASM isolates. To describe the intramammary resistome, 350 isolates of the predominant species were selected and subjected to short-read whole genome sequencing (WGS) and phenotypic antibiotic resistance profiling. While complete genomes of eight type strains were available, the remaining 15 were de novo assembled with long reads as a resource for the community. The 23 complete genomes served for reference-based assembly of the Illumina WGS data. Both chromosomes and mobile genetic elements were examined for antibiotic resistance genes (ARGs) using in-house and online software tools. ARGs were then correlated with phenotypic antibiotic resistance data from minimum inhibitory concentration (MIC). Phenotypic and genomic antimicrobial resistance was isolate-specific. Resistance to clindamycin and oxacillin was most frequently observed (65 and 30%) in Staphylococcus xylosus but could not be linked to chromosomal or plasmid-borne ARGs. However, in several cases, the observed antimicrobial resistance could be explained by the presence of mobile genetic elements like tetK carried on small plasmids. This represents a possible mechanism of transfer between non-pathogenic bacteria and pathogens of the mammary gland within and between herds. The-to our knowledge-most extensive bacteriome reported and the first attempt to link it with the resistome promise to profoundly affect veterinary bacteriology in the future and are highly relevant in a One Health context, in particular for mastitis, the treatment of which still heavily relies on antibiotics.

Retrospective surveillance of viable Bacillus cereus group contaminations in commercial food and feed vitamin B2 products sold on the Belgian market using whole-genome sequencing

Bacillus cereus is a spore-forming bacterium that occurs as a contaminant in food and feed, occasionally resulting in food poisoning through the production of various toxins. In this study, we retrospectively characterized viable B. cereus sensu lato (s.l.) isolates originating from commercial vitamin B₂ feed and food additives collected between 2016 and 2022 by the Belgian Federal Agency for the Safety of the Food Chain from products sold on the Belgian market. In total, 75 collected product samples were cultured on a general medium and, in case of bacterial growth, two isolates per product sample were collected and characterized using whole-genome sequencing (WGS) and subsequently characterized in terms of sequence type (ST), virulence gene profile, antimicrobial resistance (AMR) gene profile, plasmid content, and phylogenomic relationships. Viable B. cereus was identified in 18 of the 75 (24%) tested products, resulting in 36 WGS datasets, which were classified into eleven different STs, with ST165 (n = 10) and ST32 (n = 8) being the most common. All isolates carried multiple genes encoding virulence factors, including cytotoxin K-2 (52.78%) and cereulide (22.22%). Most isolates were predicted to be resistant to beta-lactam antibiotics (100%) and fosfomycin (88.89%), and a subset was predicted to be resistant to streptothricin (30.56%). Phylogenomic analysis revealed that some isolates obtained from different products were closely related or even identical indicating a likely common origin, whereas for some products the two isolates obtained did not show any close relationship to each other or other isolates found in other products. This study reveals that potentially pathogenic and drug-resistant B. cereus s.l. can be present in food and feed vitamin B₂ additives that are commercially available, and that more research is warranted to assess whether their presence in these types of products poses a threat to consumers.

Retrospective Screening of Anthrax-like Disease Induced by Bacillus tropicus str. JMT from Chinese Soft-Shell Turtles in Taiwan

Bacillus cereus is ubiquitous in the environment and a well-known causative agent of foodborne disease. Surprisingly, more and more emerging strains of atypical B. cereus have been identified and related to severe disease in humans and mammals such as chimpanzees, apes, and bovine. Recently, the atypical B. cereus isolates, which mainly derive from North America and Africa, have drawn great attention due to the potential risk of zoonosis. The cluster of B. cereus carries several anthrax-like virulent genes that are implicated in lethal disease. However, in non-mammals, the distribution of atypical B. cereus is still unknown. In this study, we conducted a retrospective screening of the 32 isolates of Bacillus spp. from diseased Chinese soft-shelled turtles from 2016 to 2020. To recognize the causative agent, we used various methods, such as sequencing analysis using PCR-amplification of the 16S rRNA gene, multiplex PCR for discriminating, and colony morphology by following previous studies. Furthermore, the digital DNA-DNA hybridization (dDDH) and average nucleotide identity (ANI) values were calculated, respectively, below the 70 and 96% cutoff to define species boundaries. According to the summarized results, the pathogen is taxonomically classified as Bacillus tropicus str. JMT (previous atypical Bacillus cereus). Subsequently, analyses such as targeting the unique genes using PCR and visual observation of the bacteria under various staining techniques were implemented in our study. Our findings show that all (32/32, 100%) isolates in this retrospective screening share similar phenotypical properties and carry the protective antigen (PA), edema factor (EF), hyaluronic acid (HA), and exopolysaccharide (Bps) genes on their plasmids. In this study, the results indicate that the geographic distribution and host range of B. tropicus were previously underestimated.

DNA Nanomachine (DNM) Biplex Assay for Differentiating Bacillus cereus Species

Conventional methods for the detection and differentiation of Bacillus cereus group species have drawbacks mostly due to the complexity of genetic discrimination between the Bacillus cereus species. Here, we describe a simple and straightforward assay based on the detected unamplified bacterial 16S rRNA by DNA nanomachine (DNM). The assay uses a universal fluorescent reporter and four all-DNA binding fragments, three of which are responsible for "opening up" the folded rRNA while the fourth stand is responsible for detecting single nucleotide variation (SNV) with high selectivity. Binding of the DNM to 16S rRNA results in the formation of the 10-23 deoxyribozyme catalytic core that cleaves the fluorescent reporter and produces a signal, which is amplified over time due to catalytic turnover. This developed biplex assay enables the detection of B. thuringiensis 16S rRNA at fluorescein and B. mycoides at Cy5 channels with a limit of detection of 30 × 10³ and 35 × 10³ CFU/mL, respectively, after 1.5 h with a hands-on time of ~10 min. The new assay may simplify the analysis of biological RNA samples and might be useful for environmental monitoring as a simple and inexpensive alternative to amplification-based nucleic acid analysis. The DNM proposed here may become an advantageous tool for detecting SNV in clinically significant DNA or RNA samples and can easily differentiate SNV under broadly variable experimental conditions and without prior amplification.

Elucidating the genomic history of commercially used Bacillus thuringiensis subsp. tenebrionis strain NB176

Bacillus thuringiensis subsp. tenebrionis (Btt) produces a coleopteran-specific crystal protoxin protein (Cry3Aa δ-endotoxin). After its discovery in 1982, the strain NB125 (DSM 5526) was eventually registered in 1990 to control the Colorado potato beetle (Leptinotarsa decemlineata). Gamma-irradiation of NB125 resulted in strain NB176-1 (DSM 5480) that exhibited higher cry3Aa production and became the active ingredient of the plant protection product Novodor^® FC. Here, we report a comparative genome analysis of the parental strain NB125, its derivative NB176-1 and the current commercial production strain NB176. The entire genome sequences of the parental and derivative strains were deciphered by a hybrid de novo approach using short (Illumina) and long (Nanopore) read sequencing techniques. Genome assembly revealed a chromosome of 5.4 to 5.6 Mbp and six plasmids with a size range from 14.9 to 250.5 kbp for each strain. The major differences among the original NB125 and the derivative strains NB176-1 and NB176 were an additional copy of the cry3Aa gene, which translocated to another plasmid as well as a chromosomal deletion (~ 178 kbp) in NB176. The assembled genome sequences were further analyzed in silico for the presence of virulence and antimicrobial resistance (AMR) genes.

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.

Characterization of the pathogenicity of a Bacillus cereus isolate from the Mariana Trench

Bacillus cereus is an important opportunistic pathogen widely distributed in the environment. In this study, we reported the isolation and characterization of a B. cereus isolate, MB1, from the Challenger Deep of the Mariana Trench. MB1 is aerobic, motile, and able to form endospores. It possesses 5966 genes distributed on a circular chromosome and two plasmids. The MB1 genome contains 14 sets of 23S, 5S, and 16S ribosomal RNA operons, 106 tRNA genes, 4 sRNA genes, 12 genomic islands, 13 prophages, and 302 putative virulence genes, including enterotoxins and cytolysins. Infection studies showed that MB1 was able to cause acute and lethal infection in fish and mice, and was highly toxic to mammalian cells. MB1 induced, in a dose-dependent manner, pyroptotic cell death, characterized by activation of caspase-1, cleavage of gasdermin D, and release of IL-1β and IL-18. MB1 spores exhibited swimming and haemolytic capacity, but were severely attenuated in pathogenicity, which, however, was regained to the full extent when the spores germinated under suitable conditions. Taken together, these results provide new insights into the biological and pathogenic mechanism of deep sea B. cereus.

High-Quality Genome Sequences of Six Actinobacterial Strains Isolated from Granite, Granodiorite, and Tourmaline Rock Surfaces Sampled from Tamil Nadu, India, and New England, United States

Here, we announce four contiguous and two high-quality draft genome sequences of six actinobacterial strains (Blastococcus, Georgenia, Nocardioides, Allobranchiibius, Yimella, and Williamsia) that were isolated from rock samples obtained from Indian historical ruins and colonial building stones in New England, United States. These new sequences expand the genome datasets recovered from stone-dwelling microbes and will allow the prediction of their potential role in the stone microbiome.

Draft Genome Sequences for Bacteria Associated with Root Nodules of Alnus incana in New England

Nine bacterial strains isolated from the root nodules of Alnus incana were sequenced to determine their potential roles in plant health. The selected bacterial isolates belonged to the genera Bacillus, Herbaspirillum , Pantoea , Paenibacillus , and Rothia . Here, we report the draft genome sequences.

Conjugation across Bacillus cereus and kin: A review

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

Ecophysiological and genomic analyses of a representative isolate of highly abundant Bacillus cereus strains in contaminated subsurface sediments

Bacillus cereus strain CPT56D-587-MTF (CPTF) was isolated from the highly contaminated Oak Ridge Reservation (ORR) subsurface. This site is contaminated with high levels of nitric acid and multiple heavy metals. Amplicon sequencing of the 16S rRNA genes (V4 region) in sediment from this area revealed an amplicon sequence variant (ASV) with 100% identity to the CPTF 16S rRNA sequence. Notably, this CPTF-matching ASV had the highest relative abundance in this community survey, with a median relative abundance of 3.77% and comprised 20%-40% of reads in some samples. Pangenomic analysis revealed that strain CPTF has expanded genomic content compared to other B. cereus species-largely due to plasmid acquisition and expansion of transposable elements. This suggests that these features are important for rapid adaptation to native environmental stressors. We connected genotype to phenotype in the context of the unique geochemistry of the site. These analyses revealed that certain genes (e.g. nitrate reductase, heavy metal efflux pumps) that allow this strain to successfully occupy the geochemically heterogenous microniches of its native site are characteristic of the B. cereus species while others such as acid tolerance are mobile genetic element associated and are generally unique to strain CPTF.

Identification of Universally Applicable and Species-Specific Marker Peptides for Bacillus anthracis

Anthrax is a zoonotic infection caused by the bacterium Bacillus anthracis (BA). Specific identification of this pathogen often relies on targeting genes located on two extrachromosomal plasmids, which represent the major pathogenicity factors of BA. However, more recent findings show that these plasmids have also been found in other closely related Bacillus species. In this study, we investigated the possibility of identifying species-specific and universally applicable marker peptides for BA. For this purpose, we applied a high-resolution mass spectrometry-based approach for 42 BA isolates. Along with the genomic sequencing data and by developing a bioinformatics data evaluation pipeline, which uses a database containing most of the publicly available protein sequences worldwide (UniParc), we were able to identify eleven universal marker peptides unique to BA. These markers are located on the chromosome and therefore, might overcome known problems, such as observable loss of plasmids in environmental species, plasmid loss during cultivation in the lab, and the fact that the virulence plasmids are not necessarily a unique feature of BA. The identified chromosomally encoded markers in this study could extend the small panel of already existing chromosomal targets and along with targets for the virulence plasmids, may pave the way to an even more reliable identification of BA using genomics- as well as proteomics-based techniques.

Draft Genomes Sequences of 11 Geodermatophilaceae Strains Isolated from Building Stones from New England and Indian Stone Ruins found at historic sites in Tamil Nadu, India

Metagenomic analysis of stone microbiome from samples collected in New England, USA and Tamil Nadu, India identified numerous Actinobacteria including Geodermatphilaceae. A culture-dependent approach was performed as a companion study with this culture-independent metagenomic analysis of these stone samples and resulted in the isolation of eleven Geodermatphilaceae strains (2 Geodermatophilus and 9 Blastococcus strains). The genomes of the 11 Geodermatphilaceae strains were sequenced and analyzed. The genomes for the two Geodermatophilus isolates, DF1-2 and TF2-6, were 4.45 and 4.75 Mb, respectively, while the Blastococcus genomes ranged in size from 3.98 to 5.48 Mb. Phylogenetic analysis, digital DNA:DNA hybridization (dDDH), and comparisons of the average nucleotide identities (ANI) suggest the isolates represent novel Geodermatophilus and Blastococcus species. Functional analysis of the Geodermatphilaceae genomes provides insight on the stone microbiome niche.

The SGNH hydrolase family: a template for carbohydrate diversity

The substitution and de-substitution of carbohydrate materials are important steps in the biosynthesis and/or breakdown of a wide variety of biologically important polymers. The SGNH hydrolase superfamily is a group of related and well-studied proteins with a highly conserved catalytic fold and mechanism composed of 16 member families. SGNH hydrolases can be found in vertebrates, plants, fungi, bacteria, and archaea, and play a variety of important biological roles related to biomass conversion, pathogenesis, and cell signaling. The SGNH hydrolase superfamily is chiefly composed of a diverse range of carbohydrate-modifying enzymes, including but not limited to the carbohydrate esterase families 2, 3, 6, 12 and 17 under the carbohydrate-active enzyme classification system and database (CAZy.org). In this review, we summarize the structural and functional features that delineate these subfamilies of SGNH hydrolases, and which generate the wide variety of substrate preferences and enzymatic activities observed of these proteins to date.

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.

Strains Associated with Two 2020 Welder Anthrax Cases in the United States Belong to Separate Lineages within Bacillus cereus sensu lato

Anthrax-causing members of Bacillus cereus sensu lato (s.l.) pose a serious threat to public health. While most anthrax-causing strains resemble B. anthracis phenotypically, rare cases of anthrax-like illness caused by strains resembling “B. cereus” have been reported. Here, whole-genome sequencing was used to characterize three B. cereus s.l. isolates associated with two 2020 welder anthrax cases in the United States, which resembled “B. cereus” phenotypically. Comparison of the three genomes sequenced here to all publicly available, high-quality B. cereus s.l. genomes (n = 2890 total genomes) demonstrated that genomes associated with each case effectively belonged to separate species at the conventional 95% average nucleotide identity prokaryotic species threshold. Two PubMLST sequence type 78 (ST78) genomes affiliated with a case in Louisiana were most closely related to B. tropicus and possessed genes encoding the Bps exopolysaccharide capsule, as well as hemolysin BL (Hbl) and cytotoxin K (CytK). Comparatively, a ST108 genome associated with a case in Texas was most closely related to B. anthracis; however, like other anthrax-causing strains most closely related to B. anthracis, this genome did not possess Bps-, Hbl-, or CytK-encoding genes. Overall, results presented here provide insights into the evolution of anthrax-causing B. cereus s.l.

Genomic Sequencing of Bacillus cereus Sensu Lato Strains Isolated from Meat and Poultry Products in South Africa Enables Inter- and Intranational Surveillance and Source Tracking

Members of the Bacillus cereus sensu lato species complex, also known as the B. cereus group, vary in their ability to cause illness but are frequently isolated from foods, including meat products; however, food safety surveillance efforts that use whole-genome sequencing (WGS) often neglect these potential pathogens. Here, we evaluate the surveillance and source tracking potential of WGS as applied to B. cereus sensu lato by (i) using WGS to characterize B. cereus sensu lato strains isolated during routine surveillance of meat products across South Africa (n = 25) and (ii) comparing the genomes sequenced here to all publicly available, high-quality B. cereus sensu lato genomes (n = 2,887 total genomes). Strains sequenced here were collected from meat products obtained from (i) retail outlets, processing plants, and butcheries across six South African provinces (n = 23) and (ii) imports held at port of entry (n = 2). The 25 strains sequenced here were partitioned into 15 lineages via in silico seven-gene multilocus sequence typing (MLST). While none of the South African B. cereus sensu lato strains sequenced here were identical to publicly available genomes, six MLST lineages contained multiple strains sequenced in this study, which were identical or nearly identical at the whole-genome scale (≤3 core single nucleotide polymorphisms). Five MLST lineages contained (nearly) identical genomes collected from two or three South African provinces; one MLST lineage contained nearly identical genomes from two countries (South Africa and the Netherlands), indicating that B. cereus sensu lato can spread intra- and internationally via foodstuffs. IMPORTANCE Nationwide foodborne pathogen surveillance programs that use high-resolution genomic methods have been shown to provide vast public health and economic benefits. However, Bacillus cereus sensu lato is often overlooked during large-scale routine WGS efforts. Thus, to our knowledge, no studies to date have evaluated the potential utility of WGS for B. cereus sensu lato surveillance and source tracking in foodstuffs. In this preliminary proof-of-concept study, we applied WGS to B. cereus sensu lato strains collected via South Africa's national surveillance program of domestic and imported meat products, and we provide strong evidence that B. cereus sensu lato can be disseminated intra- and internationally via the agro-food supply chain. Our results showcase that WGS has the potential to be used for source tracking of B. cereus sensu lato in foods, although future WGS and metadata collection efforts are needed to ensure that B. cereus sensu lato surveillance initiatives are on par with those of other foodborne pathogens.

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

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

Genomic insights into the diversity of non-coding RNAs in Bacillus cereus sensu lato

Bacillus cereus sensu lato is a group of bacteria of medical and agricultural importance in different ecological niches and with controversial taxonomic relationships. Studying the composition of non-coding RNAs (ncRNAs) in several bacterial groups has been an important tool for identifying genetic information and better understanding genetic regulation towards environment adaptation. However, to date, no comparative genomics study of ncRNA has been performed in this group. Thus, this study aimed to identify and characterize the set of ncRNAs from 132 strains of Bacillus cereus, Bacillus thuringiensis and Bacillus anthracis to obtain an overview of the diversity and distribution of these genetic elements in these species. We observed that the number of ncRNAs differs in the chromosomes of the three species, but not in the plasmids, when species or phylogenetic clusters were compared. The prevailing functional/structural category was Cis-reg and the most frequent class was Riboswitch. However, in plasmids, the class Group II intron was the most frequent. Also, nine ncRNAs were selected for validation in the strain B. thuringiensis 407 by RT-PCR, which allowed to identify the expression of the ncRNAs. The wide distribution and diversity of ncRNAs in the B. cereus group, and more intensely in B. thuringiensis, may help improve the abilities of these species to adapt to various environmental changes. Further studies should address the expression of these genetic elements in different conditions.

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

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

The Fate of Bacteria of the Bacillus cereus Group in the Amoeba Environment

The Bacillus cereus sensu lato group consists of several closely related species, including B. anthracis, B. cereus sensu stricto, and B. thuringiensis. Spores of these pathogenic bacteria are commonly found in the soil but evidence suggests that they are unable to grow in such a natural environment in the absence of nutrient input. Amoebas have been reported to be an amplifier for several species of pathogenic bacteria and their potential involvement to explain the large amount of B. thuringiensis and B. cereus spores in soil has been frequently proposed. Here, we studied the fate of Bacillus and amoebas when cultured together. We show that the virulence factors produced by B. thuringiensis and B. cereus do not affect the amoeba Acanthamoeba castellanii, which, on the contrary, can phagocytose and effectively digest vegetative Bacillus cells to grow and prevent the formation of cysts. Bacterial spores can germinate in the amoeba environment and the vegetative cells can then form chains or aggregates that appear to be less efficiently phagocyted by the amoeba. The use of transcriptional fusions between fluorescent reporter genes and stationary phase- and sporulation-specific promoters showed that the sporulation process occurs more efficiently in the presence of amoebas than in their absence. Moreover, our results showed the amoeba environment to promote spore germination and allow the bacteria to complete their developmental cycle. Overall, this study suggests that the amoeba-Bacillus interaction creates a virtuous circle in which each protagonist helps the other to develop.

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

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

Draft Genome Sequences of 6 Actinobacterial Strains Isolated from Rock Surfaces Obtained from Indian Stone Ruins in Tamil Nadu, India, and Rocks from New England, United States

Here, we report the draft genome sequences obtained for 6 actinobacterial strains isolated from stone surfaces acquired from New England and Indian ruins. These strains were sequenced to determine their potential functional roles in the stone microbiome. The strains belong to the genera Allobranchiibius, Agrococcus, Dermococcus, Leifsonia, and Mycobacterium.

Development of a high resolution melting method based on a novel molecular target for discrimination between Bacillus cereus and Bacillus thuringiensis

Delimitation within the Bacillus cereus group is confusing due to the highly similar genetic background of its constituent bacteria. This study aimed to develop a rapid and efficient method for the identification of Bacillus cereus and Bacillus thuringiensis, two closely related species within the B. cereus group. Using average nucleotide identity analysis (ANI) and ribosomal multilocus sequence typing (rMLST), the authenticity of the genomes of B. cereus and B. thuringiensis was determined. Emetic B. cereus and Bacillus bombysepticus were also included to provide novel genomic insights into the boundaries within the B. cereus group. Using pan-genome analysis, ispD, a novel core and single-copy molecular target, was identified for the differentiation between B. cereus and B. thuringiensis. Based on the single nucleotide polymorphism within ispD, a high resolution melting (HRM) method for the determination of B. cereus and B. thuringiensis was developed. This method can not only distinguish B. cereus and B. thuringiensis, but can also separate B. cereus from other foodborne pathogenic bacteria. The detection limit of this method could reach 1 pg of pure genomic DNA and 3.7 × 10² cfu/mL of pure culture. Moreover, this new method could effectively differentiate B. cereus and B. thuringiensis in spiked, mixed, and real food samples. Collectively, the established HRM method can provide a new reference paradigm for the sensitive and specific nucleic acid detection of pathogens with identical genomes.

Biofilm-associated heat resistance of Bacillus cereus spores in vitro and in a food model, Cheonggukjang jjigae

Bacillus cereus, a foodborne pathogen, is capable of forming spores and biofilms as methods to withstand environmental stresses. These bacterial structures are an issue for food safety as they aid the bacteria survive heat sterilisation processes of foods and food contact surfaces. This study was conducted to investigate the role of the biofilm structure in providing an extra layer of protection to spores against heat treatments. For this, heat resistance of B. cereus spores in intact biofilms was compared to that of planktonic spores in vitro and in a Cheonggukjang jjigae food model. Using methods developed in this study to measure the wet and dry heat resistance of spores in intact biofilms, it was found that B. cereus spores have significantly higher heat resistances when present in biofilms rather than as planktonic spores, and that dry heat is less effective than wet heat at killing spores in biofilms. In further detail, for wet heat treatments, spores in biofilms of the strain isolated from Cheonggukjang (Korean fermented whole soybean), B. cereus CH3, had generally higher wet heat resistances than the reference strain, B. cereus ATCC 10987, both in vitro and in the Cheonggukjang jjigae food model. However, the spores in biofilms of the two strains showed similar heat resistance to dry heat, with some exceptions, when biofilms were formed in vitro or in Cheonggukjang jjigae broth. Meanwhile, B. cereus ATCC 10987 spores in biofilms had higher or similar wet heat resistances in vitro compared to in Cheonggukjang jjigae broth. Wet heat resistances of B. cereus CH3 spores in biofilms were all statistically similar regardless of biofilm formation media (brain heart infusion and Cheonggukjang jjigae broths). For dry heat, spores in biofilms of both B. cereus strains were more heat resistant when biofilms were formed in the Cheonggukjang jjigae food model rather than in vitro. Altogether, heat resistances of spores in biofilms formed in vitro and in the food environment were found to be different depending on the tested B. cereus strain, but higher than planktonic spores in any case. This is the first study examining the heat resistance of B. cereus spores in intact biofilms matrices attached to the surface, both in vitro and in a food model. Therefore, this research is valuable to understand the protective effects of biofilms formed in food environments and to reduce the food safety risks associated with B. cereus.

Profiling of Bacillus cereus on Canadian grain

Microorganisms that cause foodborne illnesses challenge the food industry; however, environmental studies of these microorganisms on raw grain, prior to food processing, are uncommon. Bacillus cereus sensu lato is a diverse group of bacteria that is common in our everyday environment and occupy a wide array of niches. While some of these bacteria are beneficial to agriculture due to their entomopathogenic properties, others can cause foodborne illness; therefore, characterization of these bacteria is important from both agricultural and food safety standpoints. We performed a survey of wheat and flax grain samples in 2018 (n = 508) and 2017 (n = 636) and discovered that B. cereus was present in the majority of grain samples, as 56.3% and 85.2%, in two years respectively. Whole genome sequencing and comparative genomics of 109 presumptive B. cereus isolates indicates that most of the isolates were closely related and formed two genetically distinct groups. Comparisons to the available genomes of reference strains suggested that the members of these two groups are not closely related to strains previously reported to cause foodborne illness. From the same data set, another, genetically more diverse group of B. cereus was inferred, which had varying levels of similarity to previously reported strains that caused disease. Genomic analysis and PCR amplification of genes linked to toxin production indicated that most of the isolates carry the genes nheA and hbID, while other toxin genes and gene clusters, such as ces, were infrequent. This report of B. cereus on grain from Canada is the first of its kind and demonstrates the value of surveillance of bacteria naturally associated with raw agricultural commodities such as cereal grain and oilseeds.

Trichococcus shcherbakoviae subsp. psychrophilus subsp. nov., a psychrotolerant facultative anaerobe isolated from a cold spring

A psychrotolerant facultative anaerobe, strain SKBG^T, was isolated from the bottom sediments of the cold mineral spring Buxichen (Buryatia, Russia). Gram-positive non-motile cocci with a diameter of 1.75-2.5 µm were observed singly or in long chains. Cells grew in the temperature range from ̶ 5-35 °C. Growth was observed within the pH range of 7.0-9.5, with the optimum growth at pH 7.6 and at a NaCl concentration from 0-1.0 % (optimum 0.1 % (w/v)). Strain SKBG^T was a chemoorganoheterotroph that used sugars and some organic acids as substrates. The predominant fatty acids in cell walls were С_16:1ω9, С_18:1ω9, and С_16 : 0. The 16S rRNA gene sequence of strain SKBG^T shared high similarity (>99 %) with those of the type strains of the genus Trichococcus. Digital DNA-DNA hybridization (dDDH) and average nucleotide identity (ANI) values between strain SKBG^T and Trichococcus shcherbakoviae Art^T (=DSM 107162^T=VKM B-3260^T) were 70.1 and 95.4 %, respectively. The genomic DNA G+C content of strain SKBG^T was 47.1 mol%. Compared with the type strain of T. shcherbakoviae, the new strain was characterized by a temperature optimum for growth (10 °C) significantly lower than that of T. shcherbakoviae DSM 107162^T (20-30 °C). Based on phenotypic and genomic characteristics, the isolate SKBG^T was classified as T. shcherbakoviae subsp. psychrophilus subsp. nov. The type strain is SKBG^T (=VKM B-3241^Т=JCM 33326^T).

Comprehensive Genomic Analysis of the Endophytic Bacillus altitudinis Strain GLB197, a Potential Biocontrol Agent of Grape Downy Mildew

Bacillus has been extensively studied for agricultural application as a biocontrol agent. B. altitudinis GLB197, an endophytic bacterium isolated from grape leaves, exhibits distinctive inhibition to grape downy mildew based on unknown mechanisms. To determine the genetic traits involved in the mechanism of biocontrol and host-interaction traits, the genome sequence of GLB197 was obtained and further analyzed. The genome of B. altitudinis GLB197 consisted of one plasmid and a 3,733,835-bp circular chromosome with 41.56% G + C content, containing 3,770 protein-coding genes. Phylogenetic analysis of 17 Bacillus strains using the concatenated 1,226 single-copy core genes divided into different clusters was conducted. In addition, average nucleotide identity (ANI) values indicate that the current taxonomy of some B. pumilus group strains is incorrect. Comparative analysis of B. altitudinis GLB197 proteins with other B. altitudinis strains identified 3,157 core genes. Furthermore, we found that the pan-genome of B. altitudinis is open. The genome of B. altitudinis GLB197 contains one nonribosomal peptide synthetase (NRPS) gene cluster which was annotated as lichenysin. Interestingly, the cluster in B. altitudinis has two more genes than other Bacillus strains (lgrD and lgrB). The two genes were probably obtained via horizontal gene transfer (HGT) during the evolutionary process from Brevibacillus. Taken together, these observations enable the future application of B. altitudinis GLB197 as a biocontrol agent for control of grape downy mildew and promote our understanding of the beneficial interactions between B. altitudinis GLB197 and plants.

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.

Phenotypic plasticity and evolution of thermal tolerance in bacteria from temperate and hot spring environments

Phenotypic plasticity allows individuals to respond to the selective forces of a new environment, followed by adaptive evolution. We do not know to what extent phenotypic plasticity allows thermal tolerance evolution in bacteria at the border of their physiological limits. We analyzed growth and reaction norms to temperature of strains of two bacterial lineages, Bacillus cereus sensu lato and Bacillus subtilis sensu lato, that evolved in two contrasting environments, a temperate lagoon (T) and a hot spring (H). Our results showed that despite the co-occurrence of members of both lineages in the two contrasting environments, norms of reactions to temperature exhibited a similar pattern only in strains within the lineages, suggesting fixed phenotypic plasticity. Additionally, strains from the H environment showed only two to three degrees centigrade more heat tolerance than strains from the T environment. Their viability decreased at temperatures above their optimal for growth, particularly for the B. cereus lineage. However, sporulation occurred at all temperatures, consistent with the known cell population heterogeneity that allows the Bacillus to anticipate adversity. We suggest that these mesophilic strains survive in the hot-spring as spores and complete their life cycle of germination and growth during intermittent opportunities of moderate temperatures. The limited evolutionary changes towards an increase in heat tolerance in bacteria should alert us of the negative impact of climate change on all biological cycles in the planet, which at its most basic level depends on microorganisms.

Pan-Genome Portrait of Bacillus mycoides Provides Insights into the Species Ecology and Evolution

Bacillus mycoides is poorly known despite its frequent occurrence in a wide variety of environments. To provide direct insight into its ecology and evolutionary history, a comparative investigation of the species pan-genome and the functional gene categorization of 35 isolates obtained from soil samples from northeastern Poland was performed. The pan-genome of these isolates is composed of 20,175 genes and is characterized by a strong predominance of adaptive genes (∼83%), a significant amount of plasmid genes (∼37%), and a great contribution of prophages and insertion sequences. The pan-genome structure and phylodynamic studies had suggested a wide genomic diversity among the isolates, but no correlation between lineages and the bacillus origin was found. Nevertheless, the two B. mycoides populations, one from Białowieża National Park, the last European natural primeval forest with soil classified as organic, and the second from mineral soil samples taken in a farm in Jasienówka, a place with strong anthropogenic pressure, differ significantly in the frequency of genes encoding proteins enabling bacillus adaptation to specific stress conditions and production of a set of compounds, thus facilitating their colonization of various ecological niches. Furthermore, differences in the prevalence of essential stress sigma factors might be an important trail of this process. Due to these numerous adaptive genes, B. mycoides is able to quickly adapt to changing environmental conditions.

IMPORTANCE This research allows deeper understanding of the genetic organization of natural bacterial populations, specifically, Bacillus mycoides, a psychrotrophic member of the Bacillus cereus group that is widely distributed worldwide, especially in areas with continental cold climates. These thorough analyses made it possible to describe, for the first time, the B. mycoides pan-genome, phylogenetic relationship within this species, and the mechanisms behind the species ecology and evolutionary history. Our study indicates a set of functional properties and adaptive genes, in particular, those encoding sigma factors, associated with B. mycoides acclimatization to specific ecological niches and changing environmental conditions.

Bacillus cereus in Packaging Material: Molecular and Phenotypical Diversity Revealed

The Bacillus cereus group has been isolated from soils, water, plants and numerous food products. These species can produce a variety of toxins including several enterotoxins [non-hemolytic enterotoxin (Nhe), hemolysin BL (Hbl), cytotoxin K, and enterotoxin FM], the emetic toxin cereulide and insecticidal Bt toxins. This is the first study evaluating the presence of B. cereus in packaging material. Among 75 different isolates, four phylogenetic groups were detected (II, III, IV, and VI), of which the groups III and IV were the most abundant with 46.7 and 41.3%, respectively. One isolate was affiliated to psychrotolerant group VI. Growth experiments showed a mesophilic predominance. Based on PCR analysis, nhe genes were detectable in 100% of the isolates, while hbl genes were only found in 50.7%. The cereulide encoding gene was found in four out of 75 isolates, no isolate carried a crystal toxin gene. In total, thirteen different toxin gene profiles were identified. We showed that a variety of B. cereus group strains can be found in packaging material. Here, this variety lies in the presence of four phylogenetic groups, thirteen toxin gene profiles, and different growth temperatures. The results suggest that packaging material does not contain significant amounts of highly virulent strains, and the low number of cereulide producing strains is in accordance with other results.

A Whole-Genome-Based Gene-by-Gene Typing System for Standardized High-Resolution Strain Typing of Bacillus anthracis

Whole-genome sequencing (WGS) has been established for bacterial subtyping and is regularly used to study pathogen transmission, to investigate outbreaks, and to perform routine surveillance. Core-genome multilocus sequence typing (cgMLST) is a bacterial subtyping method that uses WGS data to provide a high-resolution strain characterization. This study aimed at developing a novel cgMLST scheme for Bacillus anthracis, a notorious pathogen that causes anthrax in livestock and humans worldwide. The scheme comprises 3,803 genes that were conserved in 57 B. anthracis genomes spanning the whole phylogeny. The scheme has been evaluated and applied to 584 genomes from 50 countries. On average, 99.5% of the cgMLST targets were detected. The cgMLST results confirmed the classical canonical single-nucleotide-polymorphism (SNP) grouping of B. anthracis into major clades and subclades. Genetic distances calculated based on cgMLST were comparable to distances from whole-genome-based SNP analysis with similar phylogenetic topology and comparable discriminatory power. Additionally, the application of the cgMLST scheme to anthrax outbreaks from Germany and Italy led to a definition of a cutoff threshold of five allele differences to trace epidemiologically linked strains for cluster typing and transmission analysis. Finally, the association of two clusters of B. anthracis with human cases of injectional anthrax in four European countries was confirmed using cgMLST. In summary, this study presents a novel cgMLST scheme that provides high-resolution strain genotyping for B. anthracis. This scheme can be used in parallel with SNP typing methods to facilitate rapid and harmonized interlaboratory comparisons, essential for global surveillance and outbreak analysis. The scheme is publicly available for application by users, including those with little bioinformatics knowledge.

Biodegradation of Paracetamol by Some Gram-Positive Bacterial Isolates

Bacterial isolates with the capacity to remove paracetamol were selected from an activated sludge sample collected in an oxidation ditch of a wastewater treatment plant. Among these, twelve bacterial isolates were selected according to their capacity to grow in the presence of paracetamol. They were identified using the colony morphotype procedure and by 16S rRNA gene sequencing analysis, but only four of them showed the ability to utilise paracetamol as the sole carbon source in the presence of a nitrogen supply. Those four bacterial isolates were assigned to species of the genera Bacillus, [Brevibacterium], Corynebacterium and Enterococcus. Bacterial isolates were cultured in liquid mineral salt medium (MSM) spiked with 200 mg/L of paracetamol at 28 °C in the dark. In cultures inoculated with [Brevibacterium] frigoritolerans, Corynebacterium nuruki and Enterococcus faecium, removal of 97 ± 4%, 97 ± 6% and 86.9 ± 0.8% of paracetamol at 200 mg/L were obtained, respectively, while in the presence of a species belonging to Bacillus cereus group removal of the drug below the limits of detection was attained with evidence of mineralisation, after 144 h of incubation. During the degradation process, the metabolites 4-aminophenol, hydroquinone and 2-hexenoic acid were detected. As far as we know, these species are herein first-time described as paracetamol degraders.

Current Classification of the Bacillus pumilus Group Species, the Rubber-Pathogenic Bacteria Causing Trunk Bulges Disease in Malaysia as Assessed by MLSA and Multi rep-PCR Approaches

Bacillus pumilus is the causal agent of trunk bulges disease affecting rubber and rubberwood quality and yield production. In this study, B. pumilus and other closely related species were included in B. pumilus group, as they shared over 99.5% similarity from 16S rRNA analysis. Multilocus sequence analysis (MLSA) of five housekeeping genes and repetitive elements-based polymerase chain reaction (rep-PCR) using REP, ERIC, and BOX primers conducted to analyze the diversity and systematic relationships of 20 isolates of B. pumilus group from four rubber tree plantations in Peninsular Malaysia (Serdang, Tanah Merah, Baling, and Rawang). Multi rep-PCR results revealed the genetic profiling among the B. pumilus group isolates, while MLSA results showed 98-100% similarity across the 20 isolates of B. pumilus group species. These 20 isolates, formerly established as B. pumilus, were found not to be grouped with B. pumilus. However, being distributed within distinctive groups of the B. pumilus group comprising of two clusters, A and B. Cluster A contained of 17 isolates close to B. altitudinis, whereas Cluster B consisted of three isolates attributed to B. safensis. This is the first MLSA and rep-PCR study on B. pumilus group, which provides an in-depth understanding of the diversity of these rubber-pathogenic isolates in Malaysia.

A virulent Bacillus cereus strain from deep-sea cold seep induces pyroptosis in a manner that involves NLRP3 inflammasome, JNK pathway, and lysosomal rupture

Recent studies indicate that the Bacillus species is distributed in deep-sea environments. However, no specific studies on deep-sea Bacillus cereus have been documented. In the present work, we isolated a B. cereus strain, H2, from the deep-sea cold seep in South China Sea. We characterized the pathogenic potential of H2 and investigated H2-induced death of different types of cells. We found that H2 was capable of tissue dissemination and causing acute mortality in mice and fish following intraperitoneal/intramuscular injection. In vitro studies revealed that H2 infection of macrophages induced pyroptosis and activation of the NLRP3 inflammasome pathway that contributed partly to cell death. H2 infection activated p38, JNK, and ERK, but only JNK proved to participate in H2-triggered cell death. Reactive oxygen species (ROS) and intracellular Ca²⁺ were essential to H2-induced activation of JNK and NLRP3 inflammasome. In contrast, lysosomal rupture and cathepsins were required for H2-induced NLRP3 inflammasome activation but not for JNK activation. This study revealed for the first time the virulence characteristics of deep-sea B. cereus and provided new insights into the mechanism of B. cereus infection.

The auxin-producing Bacillus thuringiensis RZ2MS9 promotes the growth and modifies the root architecture of tomato (Solanum lycopersicum cv. Micro-Tom)

Strains of Bacillus thuringiensis (Bt) are commonly commercialized as bioinoculants for insect pest control, but their benefits go beyond their insecticidal property: they can act as plant growth-promoters. Auxins play a major role in the plant growth promotion. However, the mechanism of auxin production by the Bacilli group, and more specifically by Bt strains, is unclear. In previous work, the plant growth-promoting rhizobacterium (PGPR) B. thuringiensis strain RZ2MS9 increased the corn roots. This drew our attention to the strain's auxin production trait, earlier detected in vitro. Here, we demonstrate that in its genome, RZ2MS9 harbours the complete set of genes required in two pathways that are used for Indole acetic acid (IAA) production. We also detected that the strain produces almost five times more IAA during the stationary phase. The bacterial application increased the shoot dry weight of the Micro-Tom (MT) tomato by 24%. The application also modified MT root architecture, with an increase of 26% in the average lateral root length and inhibition of the axial root. At the cellular level, RZ2MS9-treated MT plants presented elongated root cortical cells with intensified mitotic activity. Altogether, these are the best characterized auxin-associated phenotypes. Besides that, no growth alteration was detected in the auxin-insensitive diageotropic (dgt) plants either with or without the RZ2MS9 inoculation. Our results suggest that auxins play an important role in the ability of B. thuringiensis RZ2MS9 to promote MT growth and provide a better understanding of the auxin production mechanism by a Bt strain.

Dissecting the Environmental Consequences of Bacillus thuringiensis Application for Natural Ecosystems

Bacillus thuringiensis (Bt), a natural pathogen of different invertebrates, primarily insects, is widely used as a biological control agent. While Bt-based preparations are claimed to be safe for non-target organisms due to the immense host specificity of the bacterium, the growing evidence witnesses the distant consequences of their application for natural communities. For instance, upon introduction to soil habitats, Bt strains can affect indigenous microorganisms, such as bacteria and fungi, and further establish complex relationships with local plants, ranging from a mostly beneficial demeanor, to pathogenesis-like plant colonization. By exerting a direct effect on target insects, Bt can indirectly affect other organisms in the food chain. Furthermore, they can also exert an off-target activity on various soil and terrestrial invertebrates, and the frequent acquisition of virulence factors unrelated to major insecticidal toxins can extend the Bt host range to vertebrates, including humans. Even in the absence of direct detrimental effects, the exposure to Bt treatment may affect non-target organisms by reducing prey base and its nutritional value, resulting in delayed alleviation of their viability. The immense phenotypic plasticity of Bt strains, coupled with the complexity of ecological relationships they can engage in, indicates that further assessment of future Bt-based pesticides' safety should consider multiple levels of ecosystem organization and extend to a wide variety of their inhabitants.

Bacillus pumilus Group Comparative Genomics: Toward Pangenome Features, Diversity, and Marine Environmental Adaptation

Background

Members of the Bacillus pumilus group (abbreviated as the Bp group) are quite diverse and ubiquitous in marine environments, but little is known about correlation with their terrestrial counterparts. In this study, 16 marine strains that we had isolated before were sequenced and comparative genome analyses were performed with a total of 52 Bp group strains. The analyses included 20 marine isolates (which included the 16 new strains) and 32 terrestrial isolates, and their evolutionary relationships, differentiation, and environmental adaptation.

Results

Phylogenomic analysis revealed that the marine Bp group strains were grouped into three species: B. pumilus, B. altitudinis and B. safensis. All the three share a common ancestor. However, members of B. altitudinis were observed to cluster independently, separating from the other two, thus diverging from the others. Consistent with the universal nature of genes involved in the functioning of the translational machinery, the genes related to translation were enriched in the core genome. Functional genomic analyses revealed that the marine-derived and the terrestrial strains showed differences in certain hypothetical proteins, transcriptional regulators, K⁺ transporter (TrK) and ABC transporters. However, species differences showed the precedence of environmental adaptation discrepancies. In each species, land specific genes were found with possible functions that likely facilitate survival in diverse terrestrial niches, while marine bacteria were enriched with genes of unknown functions and those related to transcription, phage defense, DNA recombination and repair.

Conclusion

Our results indicated that the Bp isolates show distinct genomic features even as they share a common core. The marine and land isolates did not evolve independently; the transition between marine and non-marine habitats might have occurred multiple times. The lineage exhibited a priority effect over the niche in driving their dispersal. Certain intra-species niche specific genes could be related to a strains adaptation to its respective marine or terrestrial environment(s). In summary, this report describes the systematic evolution of 52 Bp group strains and will facilitate future studies toward understanding their ecological role and adaptation to marine and/or terrestrial environments.

Keeping up with the Bacillus cereus group: taxonomy through the genomics era and beyond

The Bacillus cereus group, also known as B. cereus sensu lato (s.l.), is a species complex that contains numerous closely related lineages, which vary in their ability to cause illness in humans and animals. The classification of B. cereus s.l. isolates into species-level taxonomic units is thus essential for informing public health and food safety efforts. However, taxonomic classification of these organisms is challenging. Numerous-often conflicting-taxonomic changes to the group have been proposed over the past two decades, making it difficult to remain up to date. In this review, we discuss the major nomenclatural changes that have accumulated in the B. cereus s.l. taxonomic space prior to 2020, particularly in the genomic sequencing era, and outline the resulting problems. We discuss several contemporary taxonomic frameworks as applied to B. cereus s.l., including (i) phenotypic, (ii) genomic, and (iii) hybrid nomenclatural frameworks, and we discuss the advantages and disadvantages of each. We offer suggestions as to how readers can avoid B. cereus s.l. taxonomic ambiguities, regardless of the nomenclatural framework(s) they choose to employ. Finally, we discuss future directions and open problems in the B. cereus s.l. taxonomic realm, including those that cannot be solved by genomic approaches alone.

Genetic diversity and the presence of circular plasmids in Bacillus cereus isolates of clinical and environmental origin

The diversity of 61 Bacillus cereus strains isolated from different clinical specimens, food including raw milk and milk products, and water was evaluated. PFGE analysis could discriminate 61 distinct pulsotypes with similarity levels from 25 to 82%, which were divided into 13 clonal complexes. The similarity between clonal complexes was at least 40%. Clinical strains were divided into 10 clonal complexes, while the strains, isolated from milk, food and water were included in 9, 6 and 6 clonal complexes, respectively. Three clonal complexes were dominated by clinical isolates, while they were absent in two complexes. Bacterial isolates from food, being a probable source of alimentary toxoinfection, showed low similarity to isolates from stool specimens. The isolates from both sources were classified together in only 4 out of 13 clonal complexes. The large circular and linear plasmids with the sizes between 50 and 200 kb were detected in 24 (39.3%) and 14 (23%) B. cereus strains, respectively. Thirteen (21.3%) strains contained only one plasmid, two plasmids were found in 6 (9.8%) of strains, and three or more plasmids were obtained in 5 (8.2%) of tested strains. The plasmids were confirmed in 30.8% and 40% of isolates from clinical specimens and food and milk samples, respectively. No clear correlation between the PFGE profiles, the source as well as plasmid content among all tested strains was observed.

Antifungal potential against Sclerotinia sclerotiorum (Lib.) de Bary and plant growth promoting abilities of Bacillus isolates from canola (Brassica napus L.) roots

Endophytic bacteria show important abilities in promoting plant growth and suppressing phytopathogens, being largely explored in agriculture as biofertilizers or biocontrol agents. Bacteria from canola roots were isolated and screened for different plant growth promotion (PGP) traits and biocontrol of Sclerotinia sclerotiorum. Thirty isolates belonging to Bacillus, Paenibacillus, Lysinibacillus, and Microbacterium genera were obtained. Several isolates produced auxin, siderophores, hydrolytic enzymes, fixed nitrogen and solubilized phosphate. Five isolates presented antifungal activity against S. sclerotiorum by the dual culture assay and four of them also inhibited fungal growth by volatile organic compounds production. All antagonistic isolates belonged to the Bacillus genus, and had their genomes sequenced for the search of biosynthetic gene clusters (BGC) related to antimicrobial metabolites. These isolates were identified as Bacillus safensis (3), Bacillus pumilus (1), and Bacillus megaterium (1), using the genomic metrics ANI and dDDH. Most strains showed several common BGCs, including bacteriocin, polyketide synthase (PKS), and non-ribosomal peptide synthetase (NRPS), related to pumilacidin, bacillibactin, bacilysin, and other antimicrobial compounds. Pumilacidin-related mass peaks were detected in acid precipitation extracts through MALDI-TOF analysis. The genomic features demonstrated the potential of these isolates in the suppression of plant pathogens; however, some aspects of plant-bacterial interactions remain to be elucidated.