Next Article in Journal
Elevated Seawater Temperature and Infection with Neoparamoeba perurans Exacerbate Complex Gill Disease in Farmed Atlantic Salmon (Salmo salar) in British Columbia, Canada
Previous Article in Journal
An Improved Method for Quantification of Viable Fusarium Cells in Infected Soil Products by Propidium Monoazide Coupled with Real-Time PCR
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 10
Issue 5
10.3390/microorganisms10051038
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Preclinical Safety Assessment of Bacillus subtilis BS50 for Probiotic and Food Applications

by
Laura M. Brutscher
1,
Claudia Borgmeier
2,
Sean M. Garvey
3 and
Jessica L. Spears
1,*
1
BIO-CAT Microbials, LLC, Shakopee, MN 55379, USA
2
BRAIN Biotech AG, 64673 Zwingenberg, Germany
3
BIO-CAT, Inc., Troy, VA 22974, USA
*
Author to whom correspondence should be addressed.
Microorganisms 2022, 10(5), 1038; https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms10051038
Submission received: 19 April 2022 / Revised: 11 May 2022 / Accepted: 13 May 2022 / Published: 17 May 2022
(This article belongs to the Section Food Microbiology)
Browse Figures
Versions Notes

Abstract

:
Despite the commercial rise of probiotics containing Bacillaceae spp., it remains important to assess the safety of each strain before clinical testing. Herein, we performed preclinical analyses to address the safety of Bacillus subtilis BS50. Using in silico analyses, we screened the 4.15 Mbp BS50 genome for genes encoding known Bacillus toxins, secondary metabolites, virulence factors, and antibiotic resistance. We also assessed the effects of BS50 lysates on the viability and permeability of cultured human intestinal epithelial cells (Caco-2). We found that the BS50 genome does not encode any known Bacillus toxins. The BS50 genome contains several gene clusters involved in the biosynthesis of secondary metabolites, but many of these antimicrobial metabolites (e.g., fengycin) are common to Bacillus spp. and may even confer health benefits related to gut microbiota health. BS50 was susceptible to seven of eight commonly prescribed antibiotics, and no antibiotic resistance genes were flanked by the complete mobile genetic elements that could enable a horizontal transfer. In cell culture, BS50 cell lysates did not diminish either Caco-2 viability or monolayer permeability. Altogether, BS50 exhibits a robust preclinical safety profile commensurate with commercial probiotic strains and likely poses no significant health risk to humans.

1. Introduction

Bacillus subtilis is a Gram-positive bacterium with a long history of use in molecular biology, industry, medicine, and fermented foods [1,2]. Bacillus strains are particularly useful for their ability to produce and secrete enzymes in mass and amenability to genetic manipulation. In the past two decades, many strains of Bacillus spp. have been used as human probiotics and direct-fed microbial for animal health. Probiotics are live microorganisms that, when administered in adequate amounts, confer a health benefit on the host [3]. Probiotics may provide health benefits such as supporting digestion, gastrointestinal (GI) health, immune health, beneficial resident gut microbes, and mood and stress response [4,5,6,7,8]. Some of the most commonly used probiotic strains include members of the Lactobacillaceae family (Bacillota phylum, formally known as Firmicutes), including the Lactiplantibacillus, Lacticaseibacillus, and Lactobacillus genera. Common probiotic strains also include Bacillus spp. and Weizmannia coagulans (formally Bacillus coagulans) strains from the Bacillaceae family of the Bacillota phylum and Bifidobacterium spp. from the Actinomycetota (formally Actinobacteria) phylum.
Bacillaceae species are well-suited for probiotic applications because they can be manufactured as spores that persist without refrigeration and resist the acidic and high bile salt conditions that occur throughout the GI tract of humans and monogastric animals [9]. Bacillus subtilis (or B. subtilis), in particular, has a history of safe consumption across the globe. B. subtilis has been used in traditional fermented foods of many East Asian cultures for centuries, including the use of B. subtilis subsp. natto for commercial production of natto, a traditional Japanese dish containing fermented soybean [10]. B. subtilis strains have also been detected in Korean kimchi, Egyptian kishk, and other cultural adaptations of fermented soy, including miso and thua nao [11,12,13,14].
In addition to work utilizing in silico and in vitro studies, several animal toxicity studies have demonstrated the safety of B. subtilis for human use [15,16,17,18,19,20]. Clinical trials of B. subtilis and W. coagulans (formally Bacillus coagulans) strain supplementations have also shown safety and tolerance in humans, as well as digestive and GI health benefits in subjects with inflammatory bowel syndrome [21,22,23,24,25], dyspepsia [26,27], as well as individuals with or without mild symptoms of GI distress [28,29,30,31,32,33,34,35,36,37,38,39,40]. For example, the B. subtilis strain MB40 has been shown to be safe and support GI health in a randomized, double-blind, placebo-controlled trial of 100 healthy adults [17,35]. Additionally, B. inaquosorum DE111 supplementation has been shown to be safe in both adult and pediatric human subjects [41,42,43,44,45,46,47,48]. Altogether, these studies provide a large body of clinical evidence that Bacillaceae spp., including B. subtilis, are safe for human consumption.
In this work, we performed preclinical studies to determine the safety of B. subtilis strain BS50 for probiotic applications. BS50 is a unique Bacillus subtilis strain that was isolated from soil and shows promise as a probiotic; preliminary assays indicate that BS50 exhibits enhanced heat tolerance and survivability in a simulated gastric model (unpublished data). To date, no known serious adverse effects have been reported from B. subtilis doses up to 10 billion colony-forming units (CFU)/day. At least five B. subtilis strains are the subject of “generally regarded as safe” (GRAS) dossiers, for which the Food and Drug Administration (FDA) has issued “no objection letters” for safe use in food [49,50,51,52,53]. Furthermore, the European Food Safety Authority (EFSA) maintains a qualified presumption of safety (QPS) list of biological agents that includes B. subtilis, which allows their use in food with no restrictions on age or exposure limit [54].
It is essential to assess the safety of each individual strain before clinical testing and safe use in dietary supplements, food, and beverages. Importantly, several Bacillus spp., including B. cereus, are capable of producing emetic toxins (e.g., cereulide), hemolytic and non-hemolytic enterotoxins, as well as cytotoxins (e.g., cytotoxin K), all of which can cause serious illness in humans and animals [55,56,57,58,59,60,61]. Another potential concern for probiotic strains is the presence of antibiotic resistance genes with flanking genetic sequences than can enable horizontal transfer to pathogenic bacteria in the GI tract [62,63,64,65,66,67]. To assess if B. subtilis BS50 poses any safety concerns to humans ahead of clinical testing, the BS50 genome was screened for genes encoding virulence factors, Bacillus toxins, and antibiotic resistance. We also performed in vitro antibiotic susceptibility tests and viability and permeability assays in human colon-derived Caco-2 cells.

2. Materials and Methods

2.1. Bacillus subtilis BS50 Isolation

B. subtilis BS50 (ATCC Accession No. PTA-127287, hereafter referred to as “BS50”) is a Gram-positive, spore-forming facultative bacterium that was isolated at BIO-CAT Microbials, LLC (Shakopee, MN, USA) from soil collected from Gallatin County, Montana, USA (collected on 4 July 2015). Isolation was performed by diluting the soil sample in Butterfield’s buffer and heating the sample up to 80 °C for 7 min to enrich for spore-forming bacteria. Serial dilutions of the sample were then plated on tryptic soy agar (TSA) plates that were incubated at 37 °C overnight. BS50 was a product of one of the resulting colonies.

2.2. Genome Sequencing

Genomic DNA was isolated from tryptic soy broth (TSB) shake flask cultures using Genomic Tip 100/G (Qiagen, Hilden, Germany) in accordance with the manufacturer’s instructions. To obtain high purity DNA appropriate for sequencing, DNA was extracted via Genomic Clean and Concentrator columns (Zymo Research, Irvine, CA, USA) and afterward checked for quality and quantity using the deNovix dsDNA Broad Range fluorometric Assay (Wilmington, DE, USA). The sample was multiplexed and pooled with other libraries using SQK-LSK109 chemistry, and Native Barcode Extension packs EXP-NBD104 and EXP-NBD114 from Oxford Nanopore Technologies (Oxford, UK). All necessary clean-up steps were carried out using Clean NA magnetic beads for next-generation sequencing (Clean NA, Waddinxveen, Netherlands). Genome sequencing took place on MinIOn FlowCells FLO-MIN106D over 48–72 h (Oxford Nanopore Technologies, Oxford, UK). The full genome was assembled with Flye [68] using default settings. The BS50 genome comprises a single, circular contig 4,150,844 bp in length. No plasmids were detected. The BS50 genome has a GC content of 43.7%.

2.3. BS50 Taxonomic Classification via Multilocus Sequence Typing

Using BLAST+ command-line software [69], the nucleotide BLAST (BLASTn) algorithm [70] was used to identify nucleotide sequences in the BS50 genome and 20 other Bacillus genomes that aligned with six genes from the genome of B. subtilis subspecies subtilis, strain 168—one of the longest existing and most extensively studied strains of B. subtilis (type strain Marburg derived) [71,72]: rpoB (GeneID: 936335), purH (GeneID: 936053), gyrA (GeneID: 940002), groEL (GeneID: 938045), polC (GeneID: 939620), and 16S rRNA (GeneID: 936895). These genes are standard “housekeeping” genes for Bacillus spp. and are commonly used for phylogenetic analysis of Bacillus spp. [73]. For each strain, the sequences aligning to these six genes were then concatenated into single nucleotide sequences (~19,616 nt). The strains used for comparison were selected based on having a complete genome available in the NCBI (National Center for Biotechnology Information) database or if they were currently used in probiotic supplements or food (i.e., MB40, BEST195, and DE111).
Multiple sequence alignment of the concatenated sequences for each Bacillus strain was performed using MAFFT [74] (accessed 10 June 2021). The multiple sequence alignment file produced by MAFFT was then input into MEGA X [75] for phylogenetic tree construction. Their evolutionary history was inferred using the Maximum Likelihood method and Tamura-Nei model [76]. Data was bootstrapped 50 times. The tree with the highest log likelihood (−30362.06) was chosen. Initial tree(s) for the heuristic search were obtained automatically by applying neighbor-joining, and BioNJ algorithms to a matrix of pairwise distances estimated using the Tamura–Nei model and then selecting the topology with a superior log-likelihood value. This analysis involved 21 nucleotide sequences. Codon positions included were 1st + 2nd + 3rd+Noncoding. There was a total of 15,093 positions in the final dataset.
In order to further characterize the sequence identities between the whole genomes of BS50 and 20 other B. subtilis strains, pairwise BLASTn alignments between BS50 and each Bacillus strain were performed via the NCBI website (accessed 25 January 2022) by uploading BS50 as the query and the other Bacillus genome as the subject. Default settings were used.

2.4. BLASTn Screen for Known Bacillus Toxins

A BLASTn search was completed via the NCBI website (accessed 2 June 2021) to determine the presence or absence of toxin genes commonly associated with the Bacillus genus. A table of the genes that were screened is shown in Table 1. In addition, positive control genes were identified in B. subtilis glutamyl-tRNA(Gln) amidotransferase subunit and B. cereus methionyl-tRNA synthetase. These genes were used as a query against the subject sequence B. subtilis BS50 genome to demonstrate the BLASTn algorithm was able to generate a match both within and across species when one existed. Each toxin gene DNA sequence was identified using NCBI gene or NCBI nucleotide databases. The sequence for the B. cereus cereulide gene cluster (cesHPTABCD) was obtained from the 270 kb plasmid pCER270 sequence (NC_010924.1, location: 15094 to 38668) [77,78]. Finally, each toxin gene DNA sequence was used as a query against the subject sequence BS50 genome. All nucleotide BLASTn alignments were run using default parameters.

2.5. BLASTx Screen for Known Bacillus Toxins

A translated nucleotide BLAST search was completed via the NCBI website (accessed 2 June 2021) to determine the presence or absence of coding sequences that are homologous to toxins commonly associated with the Bacillus genus. Protein sequences related to the control and toxin genes previously included in the BLASTn analysis were identified (http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/protein (accessed on 4 June 2021). These protein sequences were used as subjects against the query B. subtilis BS50 translated genome. All BLASTx alignments were run using default parameters.

2.6. In Silico PCR Amplification of BS50 for Bacillus Toxins

In silico PCR amplification was accessed online (4 June 2021) to search the B. subtilis BS50 genome for toxins via gene primer matches [79]. Ten sets of sequence primers for Bacillus toxin DNA amplification [80,81,82] were used to complete the virtual PCR (Supplemental Table S2). The following parameters were used to closely mimic an actual PCR run: two mismatches allowed, no mismatch allowed in the last nucleotide of the 3′ end, and a maximum band length of 10,000 nucleotides. As a positive control for the primers, the same set of primers was screened against the B. cereus genome, generating matches in all cases. As a control for the virtual PCR protocol, primers for 16S rRNA were used to show that the program would find a match when one was present.

2.7. Secondary Metabolite Screen via AntiSMASH

To determine if BS50 has the capacity to produce secondary metabolites, the BS50 genome was submitted to the online database antiSMASH bacterial version 6.0.1 (accessed 18 January 2022) [83]. Default settings were used; detection strictness was set to relaxed, and the features KnownClusterBlast, ActiveSiteFinder, RREFinder, and SubClusterBlast were turned on.

2.8. Secreted Protein via SignalP 6.0 Analysis

To determine if the BS50 genome encodes secreted proteins, it was uploaded onto the online server PATRIC [84] and annotated and translated via the RAST Tool Kit (RASTtk) [85] (accessed 26 March 2021). The translated amino acid sequences from the annotated BS50 genome were then analyzed for the presence of secreted proteins using the online SignalP 6.0 database [86] by setting the organism as “other” and setting model mode to “fast”. SignalP 6.0 was accessed on 18 January 2022. SignalP utilizes a machine learning model that predicts the presence of signal peptide motifs (i.e., Sec/SPI, Sec/SPII, Sec/SPIII, Tat/SPI, Tat/SPII) and the location of their cleavage sites [86].

2.9. Virulence Factor Screen via VFDB

To assess if the BS50 genome encodes for virulence factors (VF) or proteins involved in VF synthesis, the virulence factor database (VFDB) [87] was accessed online (17 January 2022), and the “full dataset” of VF-associated protein sequences was downloaded. The “full dataset” includes 1381 amino acid sequences for both verified and predicted VF-associated proteins from 954 medically relevant bacterial strains, whereas the “core dataset” only includes sequences of experimentally verified VF-associated proteins. The full dataset includes 36 VF-associated proteins from 164 strains and eight species of Bacillus, including proteins related to adherence (e.g., BslA), antiphagocytosis (e.g., capsule), iron acquisition, enzymes (e.g., InhA), regulation (e.g., AtxA), secretion systems (e.g., T7SS), and toxins (e.g., ALO, anthrax toxin, cereulide, certhrax, CytK, HBL, and Nhe) [87]. Since the dataset was primarily curated from medically relevant Bacillus strains, VF detection in BS50 was potentially limited. Using the BLASTx algorithm [70] with local BLAST+ command-line software [69], the BS50 genome was translated and screened against the VF dataset. Hits with <20% coverage were excluded from analysis, and multiple hits aligned to the same region of the BS50 genome were screened for the hit with the highest bit score.

2.10. Antimicrobial Resistance Gene and Mobile Genetic Element Screen

The BS50 genome was screened for antibiotic resistance factors using the Resistance Gene Identifier (RGI), which is part of the Comprehensive Antibiotic Resistance Database (CARD) [88,89]. RGI is a web-based platform that utilizes BLAST to predict complete “resistomes” from genomic and metagenomic data. The BS50 genome sequence was submitted to the RGI CARD webserver (accessed 24 April 2021) using the following criteria: Perfect, Strict, complete genes only, 95% identity nudge used. Identity nudge allows any loose hit with at least 95% identity to be scored as a strict hit.
To screen the BS50 genome for mobile genetic elements (MGE), the “A CLAssification of Mobile genetic Elements” (ACLAME) [90] database, version 0.4, was downloaded (1 June 2021) and aligned against the BS50 genome using the BLASTn [70] command with local BLAST+ software [69] under default parameters. The database contains 125,190 nucleotide sequences of predicted MGEs from prophages, virus, and bacterial plasmids. The BS50 genome was screened for known insertion sequences using the online program ISfinder [91] (accessed 1 June 2021), which utilizes the BLASTn algorithm [70] to search for nucleotide sequences that match insertion sequences.
To assess if MGEs or insertion sequences present within the BS50 could play a role in antibiotic resistance gene transfer, the loci of the sequences were manually compared to the loci of antibiotic resistance genes. Mobile genetic elements and insertion sequences that were not within five Kb of the loci of antibiotic resistance genes were not considered to play a role in antibiotic resistance gene transfer [92].

2.11. Antibiotic Minimum Inhibitory Concentration (MIC) Evaluation of BS50

MIC evaluation of BS50 against eight commonly prescribed antibiotics (i.e., chloramphenicol, clindamycin, erythromycin, gentamicin, kanamycin, streptomycin, oxytetracycline, and vancomycin) was completed by BioSciences (Bozeman, MT, USA; report number 2105336-202). The MIC of each antibiotic was determined based upon the methodology described in Clinical and Laboratory Sciences Institute (CLSI) Document M07 11th edition [93]. BS50 cells (3.93 × 106 CFU/mL per well) were exposed to each of the 10 different dilutions of each antibiotic in sterile nutrient broth. Following an appropriate incubation period, the MIC of each antibiotic was determined visually and documented. Enterococcus faecalis (ATCC Accession No. 29212) and Staphylococcus aureus (ATCC #29213) (2.96 × 106 and 8.25 × 105 CFU/mL per well, respectively) were tested in tandem with BS50 to verify the methodology performed in this study, and they exhibited MICs within the CLSI quality control range. BS50 was deemed susceptible or resistant to particular antibiotics based on specific MIC thresholds established by the European Food Safety Authority (EFSA) for Bacillus strains [94,95].

2.12. Blood Hemolysis Assay

BS50 was streaked onto sheep blood agar plates to assess its ability to lyse blood cells. After incubation overnight, the agar was inspected for alpha-or beta-hemolysis. Alpha-hemolysis, or incomplete hemolysis, is indicated by a discolored, darkened, or green medium color after test culture growth. Beta-hemolysis, or complete hemolysis, is indicated by a clear and colorless medium after growth. An indiscernible change in the color of the agar indicates that no hemolysis occurred (i.e., gamma-hemolysis).

2.13. Caco-2 Cell Viability Assay

The effects of BS50 cell lysate on Caco-2 cell viability were tested at Charles River Laboratories (Bristol, UK). Caco-2 cells are an immortalized epithelial cell line of human colorectal adenocarcinoma cells. To generate the cell lysate, BS50 cells were harvested from overnight bacterial cultures and washed. The cells were lysed via enzymatic and mechanical bead-based processes. The final lysate was filtered through a 0.2 μM filter to remove any remaining cells. The final sterile-filtered lysate was plated on TSA to ensure it was free of viable cells. A “blank” sample was used as a process control sample for the lysate production method. The blank sample was sterile, uninoculated media that was treated exactly as the lysates were, including all spins, washes, lysing, and filtering steps. To perform the assay, Caco-2 cells were harvested, counted, and plated into 96-well flat-bottomed plates at 1 × 104 cells/well in 100 µL volumes and left to adhere overnight at 37 °C, 5% CO2 in a humidified chamber. Cells were treated with BS50 lysate and incubated for an additional 48 h. Controls included cells that were left untreated and cells that were fully lysed at the time of treatment. Cell treatments were done in technical triplicate. Caco-2 cell viability was assayed using a CellTiter-Glo® intracellular ATP quantification assay (Promega, Madison, WI, USA), alongside an ATP standard curve as per the manufacturer’s guidelines. Luminescence was quantified using a GloMax® Plate reader (Promega). Levels of intracellular ATP in test conditions were quantified using the standard curve. ATP concentrations were tested for statistical significance using the Kruskal–Wallis test followed by a post-hoc Dunn’s test with Bonferroni correction for multiple testing in R Studio (Version 4.0.5). p-values less than 0.05 were considered significant.

2.14. Caco-2 Cell Transepithelial Electrical Resistance (TEER) Assay

The TEER assay was used to determine the effect of BS50 on gut barrier permeability (Charles River Laboratories, Portishead, UK). To generate a Caco-2 monolayer, Caco-2 cells were seeded on Transwell inserts over 14 days. At day 14, the polarized Caco-2 monolayers were pre-treated with a 1:5 dilution of BS50 lysate, sterile media process control, or lipopolysaccharide (LPS) control and left to incubate for 48 h. There was also a non-treatment control. TEER was measured before treatment (0 h), and at 2, 4, 6, 24, and 48 h after treatment. The TEER assays were performed twice on separate dates, with separate cell lysate preparations. Since the starting TEER values (ohm/cm2) at 0 h varied across treatments and trials, the TEER fold-changes were calculated relative to 0 h. Fold-change data from both trials were then combined and statistically analyzed as duplicates via the Kruskal–Wallis test, followed by a post-hoc Dunn’s test with Bonferroni correction for multiple testing in R Studio (Version 4.0.5). p-values less than 0.05 were considered significant.

3. Results

3.1. Taxonomic Classification of BS50

To confirm that BS50 is taxonomically a Bacillus subtilis strain, a phylogenetic tree of BS50 and 20 Bacillus strains was generated using concatenated ~20,000 nt sequences containing six “housekeeping” genes (i.e., rpoB, purH, gyrA, groEL, polC, 16S rRNA) [73]. The phylogenetic tree shows that BS50 aligns closely with other common B. subtilis strains, including the B. subtilis type strain 168 and B. subtilis MB40, a commercial probiotic strain (Figure 1). BS50 also closely aligns with commercial stains previously classified as Bacillus subtilis subsp. such as B. inaquosorum DE111. Pairwise whole genome alignments were performed between BS50 and the other Bacillus strains using BLASTn (Supplemental Table S1). Bacterial genomes sharing at least 95% average nucleotide identity are generally accepted as belonging to the same species [96,97]. The BS50 genome has 98.5% sequence identity to B. subtilis MB40 and 99.0% identity to B. subtilis subsp. natto BEST195, a B. subtilis strain commonly found in Japanese fermented natto beans (Supplemental Table S1). These data further support the classification of BS50 as a bona fide B. subtilis strain.

3.2. BLASTn Screen for Known Bacillus Toxins

To screen the BS50 genome for toxin-encoding genes, the nucleotide sequences of known Bacillus toxins were aligned against the BS50 genome using BLASTn. The control genes, gatA and metG, yielded positive matches of 98% identity with 100% sequence coverage and 71% identity with 95% sequence coverage, respectively. The metG gene from B. cereus was used as a control for cross-species sequence matches to ensure that BLASTn could identify matches within BS50 when a gene from a different species was used as the input. Because B. subtilis and B. cereus are different species, a high identity is not expected. Thus, 71% identity with 95% sequence coverage satisfies its use as a control gene for cross-species matches (Table 1). No significant similarities were found between the query toxin sequences and the BS50 genome. The identified matches, including HblA, entFM, cytK, and NheA, B, C from B. cereus and NheA, B, C from B. mycoides, were the only partial matches that covered less than 25% of the toxin gene sequences.
The B. cereus cereulide gene cluster (cesHPTABCD) from the 270 kb plasmid pCER270 sequence (NC_010924.1, location: 15094 to 38668) was also aligned against the BS50 genome. Only 50% coverage and 79% sequence identity were achieved, suggesting an incomplete cereulide gene cluster in the BS50 genome (Table 1).

3.3. BLASTx Screen for Known Bacillus Toxins

To further account for the ability of BS50 to produce toxin-encoding genes, the translated BS50 genome was aligned against the amino acid sequences of known Bacillus toxins using BLASTx.
The control proteins, GatA and MetG, yielded positive matches of 100% identity and 74.16% identity, respectively. Because B. subtilis and B. cereus are different species, a high identity was not expected, and thus, a 74.16% identity further satisfies its use as a control gene for cross-species matches (Table 2).
No significant similarities were found between the query toxin protein sequences and the translated BS50 genome. The alignment between the translated BS50 genome and EntFM from B. cereus exhibited only 52.21% identity over a span of 113 amino acids. The EntFM protein sequence is 426 amino acids long, and the alignment only covered 26.5% of the EntFM protein sequence, which is insufficient coverage to conclude that BS50 produces the EntFM protein. The BS50 genome was translated and compared to the seven proteins encoded by the B. cereus cereulide gene cluster cesHPTABCD. There were matches between the BS50 genome and the protein sequences of CesA, CesB, CesC, CesH, CesP, and CesT all of which were less than 40% identical. CesH aligned at a locus of the BS50 genome that was roughly 1.3 Mb upstream of the other cereulide biosynthesis protein alignments. There were no significant matches with CesD (Table 2).

3.4. In Silico PCR Amplification of BS50 for Bacillus Toxins

Virtual PCR only yielded matches using the positive control 16S rRNA and spoIVA primers. None of the 11 queried toxin genes were detected in the BS50 genome using virtual PCR (Supplemental Table S2).

3.5. Secondary Metabolite Screen via AntiSMASH

To determine if BS50 has the ability to produce secondary metabolites, the BS50 genome was screened for secondary metabolite biosynthetic gene clusters using the online tool, antiSMASH [83]. Ten unique secondary metabolites (two terpene hits) were predicted in the BS50 genome (Table 3).

3.6. Secreted Protein via SignalP 6.0 Analysis

To determine if the BS50 genome encodes for secreted proteins, the translated BS50 genome was analyzed for the presence of secreted proteins using the online SignalP 6.0 database [86]. As a result, 151 proteins were predicted with a greater than 50% likelihood to have Sec/SPI motifs, 93 proteins were expected to have Sec/SPII motifs, four proteins were predicted to have Tat/SPI motifs, and three proteins were predicted to have Sec/SPIII motifs.

3.7. Virulence Factor Screen via VFDB

To assess if BS50 genome encodes for virulence factors (VF), the virulence factor database (VFDB) [87] was aligned against the BS50 genome using BLASTx. There were 12 hits for VF-associated proteins in the BS50 genome (Table 4).

3.8. Antibiotic Resistance Gene Analysis

The online tool RGI was used to screen the BS50 genome for antibiotic resistance genes. RGI identified one perfect, three strict, and 275 loose hits. Of the 275 loose hits, only 12 hits had at least a 95% identity and were nudged to strict hits (Table 5). Based on the presence of a gene with roughly 98% identity to aadK, an aminoglycoside 6-adenylyltransferase that is part of the ANT6 gene family, BS50 is predicted to be resistant to streptomycin. BS50 is also predicted to be resistant to the macrolides spiramycin and telithromycin due to the presence of mph(K), a macrolide phosphotransferase. Additionally, BS50 is predicted to be resistant to tetracycline due to the presence of a tetracycline efflux pump (Tet(L)). In total, there are 16 potential resistance gene hits including aadK, mphK, and tet (45), but only seven hits that cover more than 90% of the reference gene sequence.

3.9. Insertion Sequences and Mobile Genetic Element Analysis

To assess if the antibiotic resistance genes present within the BS50 genome have the ability to be horizontally transferred to other bacteria, the BS50 genome was screened for insertion sequences using ISfinder and other mobile genetic elements using the ACLAME database (4.0). ISfinder found no matches between the BS50 genome and known insertion sequences with coverages greater than 15%. There were 122 unique loci in the BS50 genome that aligned with known mobile genetic element sequences from the ACLAME database with greater than 50% coverage, e-values less than 1.3 × 10−11, and bit scores greater than 65. To assess if these putative mobile genetic elements could play a role in antibiotic resistance gene transfer, the loci of sequences in the BS50 genome matching mobile genetic elements were then compared to the loci of antibiotic resistance genes identified via RGI. Out of the 122 loci that aligned to mobile genetic elements from the ACLAME database (4.0), one was found within five kb of an antibiotic resistance gene. The nucleotide sequence for the cupin domain-containing protein (NC_006322.1 (1,461,102 to 1,461,695)) was detected 1641 bp upstream of the blt-encoding gene (start position: 3,686,740; stop position 3,687,924). However, the nucleotide sequence for the cupin domain-containing protein only aligned to the BS50 genome with 80.3% similarity and 67% coverage, for which the 174 nt of the 5′ region did not align.

3.10. MIC Evaluation of BS50 against Eight Antibiotics

BS50 sensitivity to eight medically relevant antibiotics, including chloramphenicol, clindamycin, erythromycin, gentamicin, kanamycin, streptomycin, oxytetracycline, and vancomycin was determined by MIC methods [93]. BS50 was susceptible to seven of eight antibiotics and exhibited resistance against streptomycin (Table 6).

3.11. Blood Hemolysis Assay

To characterize any potential hemolytic activity, BS50 cells were streaked onto sheep blood agar plates and incubated overnight. The agar displayed a greenish hue surrounding the streaks where BS50 colonies grew, indicating that BS50 exhibits alpha-hemolysis.

3.12. Caco-2 Cell Viability Assay

Caco-2 cells were treated with BS50 lysate to test for deleterious effects on cell viability. While there was a significant difference in ATP concentrations between the cell lysis control and the untreated control (p = 0.014), the cells exposed to BS50 lysate showed similar ATP concentrations to the untreated control (p = 0.423) (Figure 2). Similarly, there was no significant difference in ATP concentrations between the untreated Caco-2 cell control and the blank sample, nor between the BS50 treatment and blank sample treatment.

3.13. Caco-2 Cell TEER Assay

TEER assays were performed to determine the effect of BS50 on gut barrier permeability (Figure 3). Due to variations in the initial TEER measurements across wells, fold-changes relative to 0 h from both trials were combined into one data set for statistical analysis. There were no significant differences in TEER fold-change values between the untreated control, blank process control, and cells treated with BS50 lysate at both 24 h and 48 h post-treatment (p > 0.2), whereas the LPS control lowered TEER compared to all other treatments at 24 h (p < 0.006).

4. Discussion

Spore-forming bacteria, particularly several Bacillaceae strains, are increasingly used in dietary supplements, food, and beverages due to their resistance to high temperatures and stability during manufacture, storage, and transportation [98]. Furthermore, the European Food Safety Authority (EFSA) has identified 17 Bacillaceae spp. with Qualified Presumption of Safety (QPS) status, including B. subtilis, B. amyloliquefaciens, B. licheniformis, W. coagulans, and P. megaterium, which are used as probiotics for humans and animals [54]. Regardless of the established safety of numerous Bacillaceae species, it is important to assess the safety of each individual strain, as reflected in the QPS qualifications that strains are required to meet (e.g., lack of acquired antimicrobial resistance, lack of cytotoxicity). We show here that B. subtilis strains BS50 show a robust preclinical safety profile. BS50 is a unique B. subtilis strain with at least 98% sequence similarity to commercial probiotic strains such as B. subtilis subsp. natto and B. subtilis MB40 (Supplemental Table S1).
Bacillaceae spp., such as B. anthracis, B. cereus, and B. thuringiensis, are pathogenic in humans and animals [55,56,57,58,59]. B. cereus produces the emetic toxin cereulide, enterotoxins haemolysin BL (Hbl) and non-hemolytic enterotoxin (Nhe), and cytotoxin K (CytK) [60,61]. Other strains such as B. subtilis, B. mojavensis, B. pumilus, and B. fusiformis can produce cytotoxic and emetic toxins [99,100,101]. In order to address if BS50 is capable of producing toxins, we utilized BLASTn and BLASTx to screen the BS50 genome against the nucleotide and amino acid sequences of known Bacillus toxins, including the Bacillus cereus cereulide gene cluster (cesHPTABCD, 24-kb gene cluster belonging to the 270 kb plasmid pCER270) (Table 1 and Table 2). There were matches between the translated BS50 genome and the protein sequences of CesH, CesP, CesT, CesA, CesB, and CesC, but they had less than 40% sequence identity, and they did not contiguously align throughout the genome. Further, while most of these matches aligned with greater than 90% coverage, CesA and CesB aligned with less than 65% coverage, and there were no significant matches with CesD (Table 2). Given the absence of CesD in the BS50 genome, non-contiguous alignment, and the low sequence identity and/or coverage to CesH, CesP, CesT, CesA, CesB, and CesC, there is sufficient evidence to conclude that BS50 does not contain a functioning cereulide synthase cluster.
The BS50 genome was also screened in silico for virulence factors and secondary metabolites. It was found that the BS50 genome contains secondary metabolite biosynthetic gene clusters and encodes several proteins that are associated with virulence in pathogenic organisms. However, the products encoded by these genes are not innately toxic. Contrary to primary metabolites, secondary metabolites are non-essential small organic molecules that may contribute to evolutionary fitness over time, such as improving survival against competing organisms in the same niche [102]. For example, a few secondary metabolites (e.g., bacillibactin and fengycin) that are synthesized by non-ribosomal peptide synthases (NRPS) were predicted to be produced by BS50. Bacillibactin is a catecholate siderophore encoded by the dhb operon (as detected in Table 4) and is involved in the chelation and utilization of ferric iron [103,104].
Due to its ability to bind and remove iron, bacillibactin has been proposed to treat Parkinson’s disease since patients exhibit an accumulation of iron in the brain’s substantia nigra [105]. In silico analysis also predicted that BS50 produced fengycin, an established antimicrobial in preclinical studies and suggested bioactive in a clinical observational trial; The presence of fecal Bacillus spp. was correlated with the reduced fecal occurrence of the pathogen Staphylococcus aureus in a rural Thai population [106]. Preclinical experiments suggest that fengycin production by B. subtilis is required to exert this pathogen exclusion effect [106]. Two antibiotic-encoding genes were also detected in the BS50 genome, including bacilysin and bacillaene. Bacilysin is a non-ribosomally synthesized dipeptide antibiotic that inhibits Gram-negative foodborne pathogens [107,108,109]. Bacillaene is a polyene antibiotic that can accelerate biofilm formation and has activity against a broad spectrum of bacteria, including S. aureus and E. coli [110,111,112,113]. It functions by inhibiting bacterial protein synthesis, but it cannot inhibit eukaryotic protein synthesis. BS50 also encodes for genes involved in capsular polyglutamate synthesis and transport (i.e., CapA, CapB, and CapC). Polyglutamate can enhance the pathogenesis of B. anthracis and S. epidermidis by evading the innate immune response [114,115]. Interestingly, poly-γ-glutamic acid isolated from a novel B. sonorensis strain has been shown to inhibit S. aureus and E. coli growth [116].
Most of the secondary metabolites and VF-associated proteins that were detected in the BS50 genome are also widely present throughout many Bacillus genomes [102]. As mentioned in [102], surfactin, plipastatin/fengycin, bacillibactin, bacillaene, and bacilysin are produced by 99%, 97%, 99%, 77%, and 93% of B. subtilis strains tested. Subtilosin A is also produced by several B. subtilis strains, including Strain 22a, a wild strain of B. subtilis isolated from a fermented soybean product [117,118]. All four strains of B. subtilis and no other species isolated from the mushroom substrate (including Lactococcus lactis, B. lichenimormis, and B. sonorensis) produce subtilomycin [119]. As mentioned prior, BS50 encodes genes involved in the biosynthesis of polyglutamate (Table 4). Polyglutamate is produced by many commensal Bacillus strains and is found in several Bacillus-fermented foods, including natto [120]. In a study examining polyglutamate synthesis in fermented foods, 4.7%, 1.8%, and 3.0% of the Bacillus-like strains isolated from Cheongkukjang, Doenjang, and Kochujang samples, respectively, produced polyglutamate [121]. Because these metabolites/virulence factors predicted to be synthesized by B. subtilis BS50 are produced by other species of B. subtilis, these properties should be considered intrinsic.
BS50 was also screened for the presence of antibiotic resistance encoding genes and susceptibility to antibiotics. The emergence of multidrug resistant pathogens is a major global health concern, and overuse of antibiotics has contributed to a greater incidence of antibiotic-resistant pathogens [62,122,123]. Additionally, antibiotic resistance genes present in plasmids, transposons, and integrons can be transferred from one bacteria to another via horizontal gene transfer [63,64,65,66,67]. The GI tracts of humans and animals contain complex and diverse microbial communities that may contribute to the transfer of antibiotic resistance genes from commensal organisms to potentially pathogenic bacteria [124]. BS50 was predicted to encode 16 antibiotic resistance genes that can provide resistance against multiple types of antibiotics, including fluoroquinolones, aminoglycosides, macrolides, lincosamides, tetracyclines, phenicols, nucleoside antibiotics, and peptide antibiotics (Table 5). BS50 was then tested in vitro for susceptibility/resistance against a comprehensive suite of medically relevant antibiotics as established by EFSA guidelines [94,95]; in vitro susceptibility tests determined that BS50 was resistant to the aminoglycoside streptomycin and susceptible to one phenicol antibiotic, two macrolides/lincosamides, two aminoglycosides, one glycopeptide, and oxytetracycline (Table 6).
Streptomycin resistance is widespread throughout Bacillus species and is most likely a part of their intrinsic genetic makeup rather than having acquired resistance from transferable genetic elements [125]. Regarding antibiotic resistance gene transfer, no plasmids were detected during BS50 genome assembly. While 122 regions of the BS50 genome aligned with mobile genetic elements from the ACLAME database, only one mobile genetic element was within five kb of any antibiotic resistance genes detected via CARD. The mobile genetic element cupin-domain-containing protein was detected 1641 bp upstream of the blt gene, which confers resistance against fluoroquinolone antibiotics and acridine dyes. However, 174 nt of the 5′ region of the sequence encoding for the cupin-domain-containing protein did not align to the BS50 genome, suggesting that this gene is non-functional and/or truncated. Thus, BS50 is at low risk of transferring antibiotic resistance genes to human gut-resident bacteria.
Of note, the BS50 genome encodes for a hemolysin, putative membrane hydrolase (hlyIII) (Table 4). In turn, BS50 was streaked onto sheep blood agar plates to assess its ability to lyse blood cells, and it was determined that BS50 exhibits incomplete hemolysis. Hemolytic activity has been detected throughout several Bacillus strains isolated from commercially available probiotics [126]. While this may present a safety concern if BS50 comes into contact with the bloodstream, the likelihood of an oral probiotic translocating through the intestinal barrier into the bloodstream is small and has only been reported at very low rates in hospitalized patients [127]. Nonetheless, to address potential concerns with gut barrier integrity and translocation, human colon-derived Caco-2 epithelial cell ATP viability and TEER tests were performed. We established that BS50 lysates did not negatively affect Caco-2 cell viability or monolayer permeability. Maintenance of Caco-2 cell viability and monolayer barrier integrity during BS50 lysate exposure, together with the in silico safety profile, suggest that BS50 will not be toxic to enterocytes in the human intestine or affect gut barrier integrity. A clinical trial in healthy adults has been initiated to better understand the safety and tolerability of BS50 in humans (ClinicalTrials.gov (last accessed on 18 April 2022). Identifier: NCT04655352).

5. Conclusions

Based on the results from in silico and in vitro analyses, BS50 is expected to be safe for human consumption. A clinical trial is being conducted to support the safe use of this strain by humans at anticipated rates of consumption from use in food or dietary supplements.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/microorganisms10051038/s1, Supplemental Table S1: Summary of BLASTn pairwise alignments between BS50 and 20 other Bacillus genomes; Supplemental Table S2: Summary of virtual PCR results for the absence of Bacillus toxins in BS50.

Author Contributions

Study conceptualization, J.L.S. and C.B.; methodology and data analysis, L.M.B.; manuscript preparation, L.M.B., J.L.S. and S.M.G.; manuscript review and editing, all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by BIO-CAT Microbials, LLC and BIO-CAT, Inc.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data not presented within the article or Supplementary Materials is available upon request from the corresponding author. The data are not publicly available due to privacy.

Acknowledgments

The authors thank Louise Brackenbury and Nicolas Pionnier at Charles River Discovery Research Services for their assistance in performing the Caco-2 cell ATP and TEER assays.

Conflicts of Interest

L.M.B. and J.L.S. are employees of BIO-CAT Microbials, LLC, which provided funding for this study and manufactured the investigated microbial strain. S.M.G. is an employee of BIO-CAT, Inc., which provided funding for this study. The funders were involved in the design of the study, in the collection, analyses, and interpretation of data, in the writing of the manuscript, and in the decision to publish the results. The authors declare no other conflict of interest.

References

  1. Harwood, C.R. Bacillus subtilis and its relatives: Molecular biological and industrial workhorses. Trends Biotechnol. 1992, 10, 247–256. [Google Scholar] [CrossRef]
  2. Su, Y.; Liu, C.; Fang, H.; Zhang, D. Bacillus subtilis: A universal cell factory for industry, agriculture, biomaterials and medicine. Microb. Cell Fact. 2020, 19, 173. [Google Scholar] [CrossRef]
  3. Hill, C.; Guarner, F.; Reid, G.; Gibson, G.R.; Merenstein, D.J.; Pot, B.; Morelli, L.; Canani, R.B.; Flint, H.J.; Salminen, S.; et al. Expert consensus document: The international scientific association for probiotics and prebiotics consensus statement on the scope and appropriate use of the term probiotic. Nat. Rev. Gastroenterol. Hepatol. 2014, 11, 506–514. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Patterson, E.; Griffin, S.M.; Ibarra, A.; Ellsiepen, E.; Hellhammer, J. Lacticaseibacillus paracasei Lpc-37® improves psychological and physiological markers of stress and anxiety in healthy adults: A randomized, double-blind, placebo-controlled and parallel clinical trial (the Sisu study). Neurobiol. Stress 2020, 13, 100277. [Google Scholar] [CrossRef] [PubMed]
  5. Venkataraman, R.; Madempudi, R.S.; Neelamraju, J.; Ahire, J.J.; Vinay, H.R.; Lal, A.; Thomas, G.; Stephen, S. Effect of Multi-strain Probiotic Formulation on Students Facing Examination Stress: A Double-Blind, Placebo-Controlled Study. Probiotics Antimicrob. Proteins 2021, 13, 12–18. [Google Scholar] [CrossRef]
  6. Wauters, L.; Van Oudenhove, L.; Accarie, A.; Geboers, K.; Geysen, H.; Toth, J.; Luypaerts, A.; Verbeke, K.; Smokvina, T.; Raes, J.; et al. Lactobacillus rhamnosus CNCM I-3690 decreases subjective academic stress in healthy adults: A randomized placebo-controlled trial. Gut Microbes 2022, 14, 2031695. [Google Scholar] [CrossRef]
  7. Plaza-Diaz, J.; Ruiz-Ojeda, F.J.; Gil-Campos, M.; Gil, A. Mechanisms of Action of Probiotics. Adv. Nutr. 2019, 10, S49–S66. [Google Scholar] [CrossRef] [Green Version]
  8. Sanders, M.E.; Merenstein, D.J.; Reid, G.; Gibson, G.R.; Rastall, R.A. Probiotics and prebiotics in intestinal health and disease: From biology to the clinic. Nat. Rev. Gastroenterol. Hepatol. 2019, 16, 605–616. [Google Scholar] [CrossRef]
  9. Cutting, S.M. Bacillus probiotics. Food Microbiol. 2011, 28, 214–220. [Google Scholar] [CrossRef]
  10. Shurtleff, W.; Aoyagi, A. History of Natto and Its Relatives. Soyinfo Center 2012. [Google Scholar]
  11. Jeon, H.L.; Lee, N.K.; Yang, S.J.; Kim, W.S.; Paik, H.D. Probiotic characterization of Bacillus subtilis P223 isolated from kimchi. Food Sci. Biotechnol. 2017, 26, 1641–1648. [Google Scholar] [CrossRef] [PubMed]
  12. Kotb, E. Purification and partial characterization of serine fibrinolytic enzyme from Bacillus megaterium KSK-07 isolated from kishk, a traditional Egyptian fermented food. Appl. Biochem. Microbiol. 2015, 51, 34–43. [Google Scholar] [CrossRef]
  13. Chantawannakul, P.; Oncharoen, A.; Klanbut, K.; Chukeatirote, E.; Lumyong, S. Characterization of proteases of Bacillus subtilis strain 38 isolated from traditionally fermented soybean in Northern Thailand. Sci. Asia 2002, 28, 241. [Google Scholar] [CrossRef]
  14. Inatsu, Y.; Nakamura, N.; Yuriko, Y.; Fushimi, T.; Watanasiritum, L.; Kawamoto, S. Characterization of Bacillus subtilis strains in Thua nao, a traditional fermented soybean food in northern Thailand. Lett. Appl. Microbiol. 2006, 43, 237–242. [Google Scholar] [CrossRef]
  15. Lee, N.K.; Kim, W.S.; Paik, H.D. Bacillus strains as human probiotics: Characterization, safety, microbiome, and probiotic carrier. Food Sci. Biotechnol. 2019, 28, 1297–1305. [Google Scholar] [CrossRef]
  16. Sorokulova, I.B.; Pinchuk, I.V.; Denayrolles, M.; Osipova, I.G.; Huang, J.M.; Cutting, S.M.; Urdaci, M.C. The safety of two Bacillus probiotic strains for human use. Dig. Dis. Sci. 2008, 53, 954–963. [Google Scholar] [CrossRef]
  17. Spears, J.L.; Kramer, R.; Nikiforov, A.I.; Rihner, M.O.; Lambert, E.A. Safety Assessment of Bacillus subtilis MB40 for Use in Foods and Dietary Supplements. Nutrients 2021, 13, 733. [Google Scholar] [CrossRef]
  18. Tompkins, T.A.; Hagen, K.E.; Wallace, T.D.; Fillion-Forté, V. Safety evaluation of two bacterial strains used in asian probiotic products. Can. J. Microbiol. 2008, 54, 391–400. [Google Scholar] [CrossRef]
  19. Hong, H.A.; Huang, J.M.; Khaneja, R.; Hiep, L.V.; Urdaci, M.C.; Cutting, S.M. The safety of Bacillus subtilis and Bacillus indicus as food probiotics. J. Appl. Microbiol. 2008, 105, 510–520. [Google Scholar] [CrossRef]
  20. Li, A.; Jiang, X.; Wang, Y.; Zhang, L.; Zhang, H.; Mehmood, K.; Li, Z.; Waqas, M.; Li, J. The impact of Bacillus subtilis 18 isolated from Tibetan yaks on growth performance and gut microbial community in mice. Microb. Pathog. 2019, 128, 153–161. [Google Scholar] [CrossRef]
  21. Dolin, B.J. Effects of a propietary Bacillus coagulans preparation on symptoms of diarrhea-predominant irritable bowel syndrome. Methods Find. Exp. Clin. Pharmacol. 2009, 31, 655. [Google Scholar] [CrossRef] [PubMed]
  22. Hun, L. Original Research: Bacillus coagulans Significantly Improved Abdominal Pain and Bloating in Patients with IBS. Postgrad. Med. 2009, 121, 119–124. [Google Scholar] [CrossRef] [PubMed]
  23. Majeed, M.; Nagabhushanam, K.; Natarajan, S.; Sivakumar, A.; Ali, F.; Pande, A.; Majeed, S.; Karri, S.K. Bacillus coagulans MTCC 5856 supplementation in the management of diarrhea predominant Irritable Bowel Syndrome: A double blind randomized placebo controlled pilot clinical study. Nutr. J. 2016, 15, 21. [Google Scholar] [CrossRef] [Green Version]
  24. Madempudi, R.S.; Ahire, J.J.; Neelamraju, J.; Tripathi, A.; Nanal, S. Randomized clinical trial: The effect of probiotic Bacillus coagulans Unique IS2 vs. placebo on the symptoms management of irritable bowel syndrome in adults. Sci. Rep. 2019, 9, 12210. [Google Scholar] [CrossRef]
  25. Gupta, A.K.; Maity, C. Efficacy and safety of Bacillus coagulans LBSC in irritable bowel syndrome: A prospective, interventional, randomized, double-blind, placebo-controlled clinical study [CONSORT Compliant]. Medicine 2021, 100, e23641. [Google Scholar] [CrossRef] [PubMed]
  26. Kumar, V.; Sudha, K.; Bennur, S.; Dhanasekar, K. A prospective, randomized, open-label, placebo-controlled comparative study of Bacillus coagulans GBI-30,6086 with digestive enzymes in improving indigestion in geriatric population. J. Fam. Med. Prim. Care 2020, 9, 1108. [Google Scholar] [CrossRef]
  27. Wauters, L.; Slaets, H.; De Paepe, K.; Ceulemans, M.; Wetzels, S.; Geboers, K.; Toth, J.; Thys, W.; Dybajlo, R.; Walgraeve, D.; et al. Efficacy and safety of spore-forming probiotics in the treatment of functional dyspepsia: A pilot randomised, double-blind, placebo-controlled trial. Lancet Gastroenterol. Hepatol. 2021, 6, 784–792. [Google Scholar] [CrossRef]
  28. Tompkins, T.; Xu, X.; Ahmarani, J. A comprehensive review of post-market clinical studies performed in adults with an Asian probiotic formulation. Benef. Microbes 2010, 1, 93–106. [Google Scholar] [CrossRef]
  29. Zeng, J.; Wang, C.T.; Zhang, F.S.; Qi, F.; Wang, S.F.; Ma, S.; Wu, T.J.; Tian, H.; Tian, Z.T.; Zhang, S.L.; et al. Effect of probiotics on the incidence of ventilator-associated pneumonia in critically ill patients: A randomized controlled multicenter trial. Intensive Care Med. 2016, 42, 1018–1028. [Google Scholar] [CrossRef]
  30. Maity, C.; Gupta, A.K. A prospective, interventional, randomized, double-blind, placebo-controlled clinical study to evaluate the efficacy and safety of Bacillus coagulans LBSC in the treatment of acute diarrhea with abdominal discomfort. Eur. J. Clin. Pharmacol. 2019, 75, 21–31. [Google Scholar] [CrossRef]
  31. Stecker, R.A.; Moon, J.M.; Russo, T.J.; Ratliff, K.M.; Mumford, P.W.; Jäger, R.; Purpura, M.; Kerksick, C.M. Bacillus coagulans GBI-30, 6086 improves amino acid absorption from milk protein. Nutr. Metab. 2020, 17, 93. [Google Scholar] [CrossRef] [PubMed]
  32. Tarik, M.; Ramakrishnan, L.; Bhatia, N.; Goswami, R.; Kandasamy, D.; Roy, A.; Chandran, D.S.; Singh, A.; Upadhyay, A.D.; Kalaivani, M.; et al. The effect of Bacillus coagulans Unique IS-2 supplementation on plasma amino acid levels and muscle strength in resistance trained males consuming whey protein: A double-blind, placebo-controlled study. Eur. J. Nutr. 2022. [Google Scholar] [CrossRef] [PubMed]
  33. Hanifi, A.; Culpepper, T.; Mai, V.; Anand, A.; Ford, A.L.; Ukhanova, M.; Christman, M.; Tompkins, T.A.; Dahl, W.J. Evaluation of Bacillus subtilis R0179 on gastrointestinal viability and general wellness: A randomised, double-blind, placebo-controlled trial in healthy adults. Benef. Microbes 2015, 6, 19–27. [Google Scholar] [CrossRef] [PubMed]
  34. Hatanaka, M.; Yamamoto, K.; Suzuki, N.; Iio, S.; Takara, T.; Morita, H.; Takimoto, T.; Nakamura, T. Effect of Bacillus subtilis C-3102 on loose stools in healthy volunteers. Benef. Microbes 2018, 9, 357–365. [Google Scholar] [CrossRef] [PubMed]
  35. Penet, C.; Kramer, R.; Little, R.; Spears, J.L.; Parker, J.; Iyer, J.K.; Guthrie, N.; Evans, M. A Randomized, Double-blind, Placebo-controlled, Parallel Study Evaluating the Efficacy of Bacillus subtilis MB40 to Reduce Abdominal Discomfort, Gas, and Bloating. Altern. Ther. Health Med. 2021, 27, 146–157. [Google Scholar]
  36. Hatanaka, M.; Kanzato, H.; Tsuda, R.; Nadaoka, I.; Yasue, M.; Hoshino, T.; Iio, S.; Takara, T. Safety evaluation of the excessive intake of Bacillus subtilis C-3102 in healthy Japanese adults: A randomized, placebo-controlled, double-blind, parallel-group, comparison trial. Toxicol. Rep. 2020, 7, 46–58. [Google Scholar] [CrossRef]
  37. Lefevre, M.; Racedo, S.M.; Ripert, G.; Housez, B.; Cazaubiel, M.; Maudet, C.; Jüsten, P.; Marteau, P.; Urdaci, M.C. Probiotic strain Bacillus subtilis CU1 stimulates immune system of elderly during common infectious disease period: A randomized, double-blind placebo-controlled study. Immun. Ageing 2015, 12, 24. [Google Scholar] [CrossRef] [Green Version]
  38. Horosheva, T.V.; Vodyanoy, V.; Sorokulova, I. Efficacy of Bacillus probiotics in prevention of antibiotic-associated diarrhoea: A randomized, double-blind, placebo-controlled clinical trial. JMM Case Rep. 2014, 1, e004036. [Google Scholar] [CrossRef] [Green Version]
  39. Lefevre, M.; Racedo, S.M.; Denayrolles, M.; Ripert, G.; Desfougères, T.; Lobach, A.R.; Simon, R.; Pélerin, F.; Jüsten, P.; Urdaci, M.C. Safety assessment of Bacillus subtilis CU1 for use as a probiotic in humans. Regul. Toxicol. Pharmacol. 2017, 83, 54–65. [Google Scholar] [CrossRef] [Green Version]
  40. Kalman, D.S.; Schwartz, H.I.; Alvarez, P.; Feldman, S.; Pezzullo, J.C.; Krieger, D.R. A prospective, randomized, double-blind, placebo-controlled parallel-group dual site trial to evaluate the effects of a Bacillus coagulans-based product on functional intestinal gas symptoms. BMC Gastroenterol. 2009, 9, 85. [Google Scholar] [CrossRef] [Green Version]
  41. Freedman, K.E.; Hill, J.L.; Wei, Y.; Vazquez, A.R.; Grubb, D.S.; Trotter, R.E.; Wrigley, S.D.; Johnson, S.A.; Foster, M.T.; Weir, T.L. Examining the gastrointestinal and immunomodulatory effects of the novel probiotic Bacillus subtilis de111. Int. J. Mol. Sci. 2021, 22, 2453. [Google Scholar] [CrossRef]
  42. Paytuví-Gallart, A.; Sanseverino, W.; Winger, A.M. Daily intake of probiotic strain Bacillus subtilis DE111 supports a healthy microbiome in children attending day-care. Benef. Microbes 2020, 11, 611–620. [Google Scholar] [CrossRef]
  43. Toohey, J.C.; Townsend, J.R.; Johnson, S.B.; Toy, A.M.; Vantrease, W.C.; Bender, D.; Crimi, C.C.; Stowers, K.L.; Ruiz, M.D.; VanDusseldorp, T.A.; et al. Effects of Probiotic (Bacillus subtilis) Supplementation During Offseason Resistance Training in Female Division I Athletes. J. Strength Cond. Res. 2020, 34, 3173–3181. [Google Scholar] [CrossRef] [PubMed]
  44. Townsend, J.R.; Bender, D.; Vantrease, W.C.; Sapp, P.A.; Toy, A.M.; Woods, C.A.; Johnson, K.D. Effects of probiotic (Bacillus subtilis de111) supplementation on immune function, hormonal status, and physical performance in division i baseball players. Sports 2018, 6, 70. [Google Scholar] [CrossRef] [Green Version]
  45. Townsend, J.R.; Vantrease, W.C.; Jones, M.D.; Sapp, P.A.; Johnson, K.D.; Beuning, C.N.; Haase, A.A.; Boot, C.M. Plasma amino acid response to whey protein ingestion following 28 days of probiotic (Bacillus subtilis de111) supplementation in active men and women. J. Funct. Morphol. Kinesiol. 2021, 6, 1. [Google Scholar] [CrossRef] [PubMed]
  46. Trotter, R.E.; Vazquez, A.R.; Grubb, D.S.; Freedman, K.E.; Grabos, L.E.; Jones, S.; Gentile, C.L.; Melby, C.L.; Johnson, S.A.; Weir, T.L. Bacillus subtilis DE111 intake may improve blood lipids and endothelial function in healthy adults. Benef. Microbes 2020, 11, 621–630. [Google Scholar] [CrossRef] [PubMed]
  47. Colom, J.; Freitas, D.; Simon, A.; Brodkorb, A.; Buckley, M.; Deaton, J.; Deaton, J. Presence and Germination of the Probiotic Bacillus subtilis DE111® in the Human Small Intestinal Tract: A Randomized, Crossover, Double- Blind, and Placebo-Controlled Study. Front. Microbiol. 2021, 12, 2189. [Google Scholar] [CrossRef] [PubMed]
  48. Cuentas, A.M.; Deaton, J.; Khan, S.; Davidson, J.; Ardita, C. The Effect of Bacillus subtilis DE111 on the Daily Bowel Movement Profile for People with Occasional Gastrointestinal Irregularity. J. Probiotics Health 2017, 5, 4. [Google Scholar] [CrossRef]
  49. FDA GRN. GRAS Notification 905. Available online: https://www.fda.gov/media/139501/download (accessed on 18 April 2022).
  50. FDA GRN. GRAS Notification 831. Available online: https://www.fda.gov/media/132389/download (accessed on 18 April 2022).
  51. FDA GRN. GRAS Notification 955. Available online: https://www.fda.gov/media/146438/download (accessed on 18 April 2022).
  52. FDA GRN. GRAS Notification 969. Available online: https://fda.report/media/150352/GRAS-Notice-GRN-969-Bacillus-subtilis.pdf (accessed on 18 April 2022).
  53. FDA GRN. GRAS Notification 956. Available online: https://www.fda.gov/media/146998/download (accessed on 18 April 2022).
  54. Koutsoumanis, K.; Allende, A.; Alvarez-Ordóñez, A.; Bolton, D.; Bover-Cid, S.; Chemaly, M.; Davies, R.; De Cesare, A.; Hilbert, F.; Lindqvist, R.; et al. Update of the list of QPS-recommended biological agents intentionally added to food or feed as notified to EFSA 14: Suitability of taxonomic units notified to EFSA until March 2021. EFSA J. 2021, 19, e06689. [Google Scholar] [CrossRef]
  55. Damgaard, P.H.; Granum, P.E.; Bresciani, J.; Torregrossa, M.V.; Eilenberg, J.; Valentino, L. Characterization of Bacillus thuringiensis isolated from infections in burn wounds. FEMS Immunol. Med. Microbiol. 1997, 18, 47–53. [Google Scholar] [CrossRef] [Green Version]
  56. Hernandez, E.; Ramisse, F.; Ducoureau, J.P.; Cruel, T.; Cavallo, J.D. Bacillus thuringiensis subsp. konkukian (Serotype H34) superinfection: Case report and experimental evidence of pathogenicity in immunosuppressed mice. J. Clin. Microbiol. 1998, 36, 2138–2139. [Google Scholar] [CrossRef] [Green Version]
  57. Little, S.F.; Ivins, B.E. Molecular pathogenesis of Bacillus anthracis infection. Microbes Infect. 1999, 1, 131–139. [Google Scholar] [CrossRef]
  58. Kotiranta, A.; Lounatmaa, K.; Haapasalo, M. Epidemiology and pathogenesis of Bacillus cereus infections. Microbes Infect. 2000, 2, 189–198. [Google Scholar] [CrossRef]
  59. Raymond, B.; Johnston, P.R.; Nielsen-LeRoux, C.; Lereclus, D.; Crickmore, N. Bacillus thuringiensis: An impotent pathogen? Trends Microbiol. 2010, 18, 189–194. [Google Scholar] [CrossRef] [PubMed]
  60. Dietrich, R.; Jessberger, N.; Ehling-Schulz, M.; Märtlbauer, E.; Granum, P.E. The Food Poisoning Toxins of Bacillus cereus. Toxins 2021, 13, 98. [Google Scholar] [CrossRef] [PubMed]
  61. Schoeni, J.L.; Lee Wong, A.C. Bacillus cereus Food Poisoning and Its Toxins. J. Food Prot. 2005, 68, 636–648. [Google Scholar] [CrossRef] [PubMed]
  62. Sengupta, S.; Chattopadhyay, M.K.; Grossart, H.P. The multifaceted roles of antibiotics and antibiotic resistance in nature. Front. Microbiol. 2013, 4, 47. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Alekshun, M.N.; Levy, S.B. Molecular Mechanisms of Antibacterial Multidrug Resistance. Cell 2007, 128, 1037–1050. [Google Scholar] [CrossRef] [Green Version]
  64. van Reenen, C.A.; Dicks, L.M.T. Horizontal gene transfer amongst probiotic lactic acid bacteria and other intestinal microbiota: What are the possibilities? A review. Arch. Microbiol. 2011, 193, 157–168. [Google Scholar] [CrossRef]
  65. Santagati, M.; Campanile, F.; Stefani, S. Genomic diversification of enterococci in hosts: The role of the mobilome. Front. Microbiol. 2012, 3, 95. [Google Scholar] [CrossRef] [Green Version]
  66. Blair, J.M.A.; Webber, M.A.; Baylay, A.J.; Ogbolu, D.O.; Piddock, L.J.V. Molecular mechanisms of antibiotic resistance. Nat. Rev. Microbiol. 2015, 13, 42–51. [Google Scholar] [CrossRef] [PubMed]
  67. Wozniak, R.A.F.; Waldor, M.K. Integrative and conjugative elements: Mosaic mobile genetic elements enabling dynamic lateral gene flow. Nat. Rev. Microbiol. 2010, 8, 552–563. [Google Scholar] [CrossRef] [PubMed]
  68. Kolmogorov, M.; Yuan, J.; Lin, Y.; Pevzner, P.A. Assembly of long, error-prone reads using repeat graphs. Nat. Biotechnol. 2019, 37, 540–546. [Google Scholar] [CrossRef] [PubMed]
  69. Camacho, C.; Coulouris, G.; Avagyan, V.; Ma, N.; Papadopoulos, J.; Bealer, K.; Madden, T.L. BLAST+: Architecture and applications. BMC Bioinform. 2009, 10, 421. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Altschul, S.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  71. Burkholder, P.R.; Giles, N.H. Induced biochemical mutations in Bacillus subtilis. Am. J. Bot. 1947, 34, 345–348. [Google Scholar] [CrossRef]
  72. Zeigler, D.R.; Prágai, Z.; Rodriguez, S.; Chevreux, B.; Muffler, A.; Albert, T.; Bai, R.; Wyss, M.; Perkins, J.B. The origins of 168, W23, and other Bacillus subtilis legacy strains. J. Bacteriol. 2008, 190, 6983–6995. [Google Scholar] [CrossRef] [Green Version]
  73. Kubo, Y.; Rooney, A.P.; Tsukakoshi, Y.; Nakagawa, R.; Hasegawa, H.; Kimura, K. Phylogenetic analysis of Bacillus subtilis strains applicable to natto (fermented soybean) production. Appl. Environ. Microbiol. 2011, 77, 6463–6469. [Google Scholar] [CrossRef] [Green Version]
  74. Katoh, K.; Rozewicki, J.; Yamada, K.D. MAFFT online service: Multiple sequence alignment, interactive sequence choice and visualization. Brief. Bioinform. 2018, 20, 1160–1166. [Google Scholar] [CrossRef] [Green Version]
  75. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  76. Tamura, K.; Nei, M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol. Biol. Evol. 1993, 10, 512–526. [Google Scholar] [CrossRef] [PubMed]
  77. Mei, X.; Xu, K.; Yang, L.; Yuan, Z.; Mahillon, J.; Hu, X. The genetic diversity of cereulide biosynthesis gene cluster indicates a composite transposon Tnces in emetic Bacillus weihenstephanensis. BMC Microbiol. 2014, 14, 149. [Google Scholar] [CrossRef] [Green Version]
  78. Ehling-Schulz, M.; Fricker, M.; Grallert, H.; Rieck, P.; Wagner, M.; Scherer, S. Cereulide synthetase gene cluster from emetic Bacillus cereus: Structure and location on a mega virulence plasmid related to Bacillus anthracis toxin plasmid pXO1. BMC Microbiol. 2006, 6, 20. [Google Scholar] [CrossRef] [Green Version]
  79. San Millán, R.M.; Martínez-Ballesteros, I.; Rementeria, A.; Garaizar, J.; Bikandi, J. Online exercise for the design and simulation of PCR and PCR-RFLP experiments. BMC Res. Notes 2013, 6, 2–5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Agata, N.; Ohta, M.; Mori, M.; Isobe, M. A novel dodecadepsipeptide, cereulide, is an emetic toxin of Bacillus cereus. FEMS Microbiol. Lett. 1995, 129, 17–19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Asano, S.I.; Nukumizu, Y.; Bando, H.; Iizuka, T.; Yamamoto, T. Cloning of novel enterotoxin genes from Bacillus cereus and Bacillus thuringiensis. Appl. Environ. Microbiol. 1997, 63, 1054–1057. [Google Scholar] [CrossRef] [Green Version]
  82. Mäntynen, V.; Lindström, K. A rapid PCR-based DNA test for enterotoxic Bacillus cereus. Appl. Environ. Microbiol. 1998, 64, 1634–1639. [Google Scholar] [CrossRef] [Green Version]
  83. Blin, K.; Shaw, S.; Kloosterman, A.M.; Charlop-Powers, Z.; van Wezel, G.P.; Medema, M.H.; Weber, T. antiSMASH 6.0: Improving cluster detection and comparison capabilities. Nucleic Acids Res. 2021, 49, W29–W35. [Google Scholar] [CrossRef]
  84. Wattam, A.R.; Abraham, D.; Dalay, O.; Disz, T.L.; Driscoll, T.; Gabbard, J.L.; Gillespie, J.J.; Gough, R.; Hix, D.; Kenyon, R.; et al. PATRIC, the bacterial bioinformatics database and analysis resource. Nucleic Acids Res. 2014, 42, 581–591. [Google Scholar] [CrossRef] [Green Version]
  85. Brettin, T.; Davis, J.J.; Disz, T.; Edwards, R.A.; Gerdes, S.; Olsen, G.J.; Olson, R.; Overbeek, R.; Parrello, B.; Pusch, G.D.; et al. RASTtk: A modular and extensible implementation of the RAST algorithm for building custom annotation pipelines and annotating batches of genomes. Sci. Rep. 2015, 5, 8365. [Google Scholar] [CrossRef] [Green Version]
  86. Teufel, F.; Almagro Armenteros, J.J.; Johansen, A.R.; Gíslason, M.H.; Pihl, S.I.; Tsirigos, K.D.; Winther, O.; Brunak, S.; von Heijne, G.; Nielsen, H. SignalP 6.0 predicts all five types of signal peptides using protein language models. Nat. Biotechnol. 2022. [Google Scholar] [CrossRef] [PubMed]
  87. Liu, B.; Zheng, D.; Jin, Q.; Chen, L.; Yang, J. VFDB 2019: A comparative pathogenomic platform with an interactive web interface. Nucleic Acids Res. 2019, 47, D687–D692. [Google Scholar] [CrossRef]
  88. McArthur, A.G.; Waglechner, N.; Nizam, F.; Yan, A.; Azad, M.A.; Baylay, A.J.; Bhullar, K.; Canova, M.J.; De Pascale, G.; Ejim, L.; et al. The comprehensive antibiotic resistance database. Antimicrob. Agents Chemother. 2013, 57, 3348–3357. [Google Scholar] [CrossRef] [Green Version]
  89. Alcock, B.P.; Raphenya, A.R.; Lau, T.T.Y.; Tsang, K.K.; Bouchard, M.; Edalatmand, A.; Huynh, W.; Nguyen, A.L.V.; Cheng, A.A.; Liu, S.; et al. CARD 2020: Antibiotic resistome surveillance with the comprehensive antibiotic resistance database. Nucleic Acids Res. 2020, 48, D517–D525. [Google Scholar] [CrossRef] [PubMed]
  90. Leplae, R.; Lima-Mendez, G.; Toussaint, A. ACLAME: A CLAssification of mobile genetic elements, update 2010. Nucleic Acids Res. 2009, 38, 57–61. [Google Scholar] [CrossRef] [Green Version]
  91. Siguier, P.; Perochon, J.; Lestrade, L.; Mahillon, J.; Chandler, M. ISfinder: The reference centre for bacterial insertion sequences. Nucleic Acids Res. 2006, 34, 32–36. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  92. Ben Maamar, S.; Glawe, A.J.; Brown, T.K.; Hellgeth, N.; Hu, J.; Wang, J.P.; Huttenhower, C.; Hartmann, E.M. Mobilizable antibiotic resistance genes are present in dust microbial communities. PLoS Pathog. 2020, 16, e1008211. [Google Scholar] [CrossRef]
  93. Weinstein, M.P.; Patel, J.B.; Burnhman, C.-A.; ZImmer, B.L. M07 Methods for Dilution Antimicrobial Susceptibility Tests for Bacteria That Grow Aerobically; Clinical and Laboratory Standards Institute: Wayne, PA, USA, 2018; p. 91. [Google Scholar]
  94. Rychen, G.; Aquilina, G.; Azimonti, G.; Bampidis, V.; Bastos, M.D.L.; Bories, G.; Chesson, A.; Cocconcelli, P.S.; Flachowsky, G.; Gropp, J.; et al. Guidance on the characterisation of microorganisms used as feed additives or as production organisms. EFSA J. 2018, 16, e05206. [Google Scholar] [CrossRef]
  95. European Food Safety Authority (EFSA) European Food Safety Authority, EFSA Panel on Additives and Products or Substances used in Animal Feed (FEEDAP). Scientific opinion Guidance on the assessment of bacterial susceptibility to antimicrobials of human and veterinary importance. EFSA J. 2012, 10, 2740. [Google Scholar] [CrossRef]
  96. Konstantinidis, K.T.; Tiedje, J.M. Genomic insights that advance the species definition for prokaryotes. Proc. Natl. Acad. Sci. USA 2005, 102, 2567–2572. [Google Scholar] [CrossRef] [Green Version]
  97. Richter, M.; Rosselló-Móra, R. Shifting the genomic gold standard for the prokaryotic species definition. Proc. Natl. Acad. Sci. USA 2009, 106, 19126–19131. [Google Scholar] [CrossRef] [Green Version]
  98. Chaucheyras-Durand, F.; Durand, H. Probiotics in Animal Nutrition and Health. Benef. Microbes 2016, 1, 3–9. [Google Scholar] [CrossRef] [PubMed]
  99. From, C.; Pukall, R.; Schumann, P.; Hormazábal, V.; Granum, P.E. Toxin-producing ability among Bacillus spp. outside the Bacillus cereus group. Appl. Environ. Microbiol. 2005, 71, 1178–1183. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  100. Salkinoja-Salonen, M.S.; Vuorio, R.; Andersson, M.A.; Kämpfer, P.; Andersson, M.C.; Honkanen-Buzalski, T.; Scoging, A.C. Toxigenic strains of Bacillus licheniformis related to food poisoning. Appl. Environ. Microbiol. 1999, 65, 4637–4645. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  101. Taylor, J.M.W.; Sutherland, A.D.; Aidoo, K.E.; Logan, N.A. Heat-stable toxin production by strains of Bacillus cereus, Bacillus firmus, Bacillus megaterium, Bacillus simplex and Bacillus licheniformis. FEMS Microbiol. Lett. 2005, 242, 313–317. [Google Scholar] [CrossRef] [Green Version]
  102. Harwood, C.R.; Mouillon, J.M.; Pohl, S.; Arnau, J. Secondary metabolite production and the safety of industrially important members of the Bacillus subtilis group. FEMS Microbiol. Rev. 2018, 42, 721–738. [Google Scholar] [CrossRef]
  103. Raza, W.; Hussain, Q.; Shen, Q. Production, regulation and transportation of bacillibactin in Bacillus subtilis. J. Chem. Soc. Pak. 2012, 34, 996–1002. [Google Scholar]
  104. May, J.J.; Wendrich, T.M.; Marahiel, M.A. The dhb Operon of Bacillus subtilis Encodes the Biosynthetic Template for the Catecholic Siderophore 2,3-Dihydroxybenzoate-Glycine-Threonine Trimeric Ester Bacillibactin. J. Biol. Chem. 2001, 276, 7209–7217. [Google Scholar] [CrossRef] [Green Version]
  105. Kaur, D.; Yantiri, F.; Rajagopalan, S.; Kumar, J.; Mo, J.Q.; Boonplueang, R.; Viswanath, V.; Jacobs, R.; Yang, L.; Beal, M.F.; et al. Genetic or Pharmacological Iron Chelation Prevents MPTP-Induced Neurotoxicity In Vivo. Neuron 2003, 37, 899–909. [Google Scholar] [CrossRef] [Green Version]
  106. Piewngam, P.; Zheng, Y.; Nguyen, T.H.; Dickey, S.W.; Joo, H.-S.; Villaruz, A.E.; Glose, K.A.; Fisher, E.L.; Hunt, R.L.; Li, B.; et al. Pathogen elimination by probiotic Bacillus via signalling interference. Nature 2018, 562, 532–537. [Google Scholar] [CrossRef]
  107. Kenig, M.; Abraham, E.P. Antimicrobial Activities and Antagonists of Bacilysin and Anticapsin. J. Gen. Microbiol. 1976, 94, 37–45. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  108. Nannan, C.; Vu, H.Q.; Gillis, A.; Caulier, S.; Nguyen, T.T.T.; Mahillon, J. Bacilysin within the Bacillus subtilis group: Gene prevalence versus antagonistic activity against Gram-negative foodborne pathogens. J. Biotechnol. 2021, 327, 28–35. [Google Scholar] [CrossRef] [PubMed]
  109. Steinborn, G.; Hajirezaei, M.R.; Hofemeister, J. bac genes for recombinant bacilysin and anticapsin production in Bacillus host strains. Arch. Microbiol. 2005, 183, 71–79. [Google Scholar] [CrossRef]
  110. Chen, X.H.; Vater, J.; Piel, J.; Franke, P.; Scholz, R.; Schneider, K.; Koumoutsi, A.; Hitzeroth, G.; Grammel, N.; Strittmatter, A.W.; et al. Structural and functional characterization of three polyketide synthase gene clusters in Bacillus amyloliquefaciens FZB 42. J. Bacteriol. 2006, 188, 4024–4036. [Google Scholar] [CrossRef] [Green Version]
  111. Patel, P.S.; Huang, S.; Fisher, S.; Pirnik, D.; Aklonis, C.; Dean, L.; Meyers, E.; Fernandes, P.; Mayerl, F. Bacillaene, a novel inhibitor of procaryotic protein synthesis produced by Bacillus subtilis: Production, taxonomy, isolation, physico-chemical characterization and biological activity. J. Antibiot. 1995, 48, 997–1003. [Google Scholar] [CrossRef] [Green Version]
  112. Li, H.; Han, X.; Dong, Y.; Xu, S.; Chen, C.; Feng, Y.; Cui, Q.; Li, W. Bacillaenes: Decomposition Trigger Point and Biofilm Enhancement in Bacillus. ACS Omega 2021, 6, 1093–1098. [Google Scholar] [CrossRef] [PubMed]
  113. Butcher, R.A.; Schroeder, F.C.; Fischbach, M.A.; Straight, P.D.; Kolter, R.; Walsh, C.T.; Clardy, J. The identification of bacillaene, the product of the PksX megacomplex in Bacillus subtilis. Proc. Natl. Acad. Sci. USA 2007, 104, 1506–1509. [Google Scholar] [CrossRef] [Green Version]
  114. Sharma, S.; Bhatnagar, R.; Gaur, D. Bacillus anthracis Poly-γ-D-Glutamate Capsule Inhibits Opsonic Phagocytosis by Impeding Complement Activation. Front. Immunol. 2020, 11, 462. [Google Scholar] [CrossRef]
  115. Kocianova, S.; Vuong, C.; Yao, Y.; Voyich, J.M.; Fischer, E.R.; DeLeo, F.R.; Otto, M. Key role of poly-γ-DL-glutamic acid in immune evasion and virulence of Staphylococcus epidermidis. J. Clin. Investig. 2005, 115, 688–694. [Google Scholar] [CrossRef]
  116. Ijadi Bajestani, M.; Mousavi, S.M.; Mousavi, S.B.; Jafari, A.; Shojaosadati, S.A. Purification of extra cellular poly-γ-glutamic acid as an antibacterial agent using anion exchange chromatography. Int. J. Biol. Macromol. 2018, 113, 142–149. [Google Scholar] [CrossRef]
  117. Stein, T.; Düsterhus, S.; Stroh, A.; Entian, K.D. Subtilosin Production by Two Bacillus subtilis Subspecies and Variance of the sbo-alb Cluster. Appl. Environ. Microbiol. 2004, 70, 2349–2353. [Google Scholar] [CrossRef] [Green Version]
  118. Zheng, G.; Slavik, M.F. Isolation, partial purification and characterization of a bacteriocin produced by a newly isolated Bacillus subtilis strain. Lett. Appl. Microbiol. 1999, 28, 363–367. [Google Scholar] [CrossRef] [PubMed]
  119. Bolocan, A.S.; Pennone, V.; O’Connor, P.M.; Coffey, A.; Nicolau, A.I.; McAuliffe, O.; Jordan, K. Inhibition of Listeria monocytogenes biofilms by bacteriocin-producing bacteria isolated from mushroom substrate. J. Appl. Microbiol. 2017, 122, 279–293. [Google Scholar] [CrossRef]
  120. Bhat, A.R.; Irorere, V.U.; Bartlett, T.; Hill, D.; Kedia, G.; Morris, M.R.; Charalampopoulos, D.; Radecka, I. Bacillus subtilis natto: A non-toxic source of poly-γ-glutamic acid that could be used as a cryoprotectant for probiotic bacteria. AMB Express 2013, 3, 36. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  121. Kang, S.-E.; Rhee, J.-H.; Park, C.; Sung, M.-H.; Lee, I.-H. Distribution of poly-gamma-glutamate producers in Korean fermented foods, Cheongkukjang, Doenjang, and Kochujang. Food Sci. Biotechnol. 2005, 14, 704–708. [Google Scholar]
  122. Laxminarayan, R.; Duse, A.; Wattal, C.; Zaidi, A.K.M.; Wertheim, H.F.L.; Sumpradit, N.; Vlieghe, E.; Hara, G.L.; Gould, I.M.; Goossens, H.; et al. Antibiotic resistance—The need for global solutions. Lancet Infect. Dis. 2013, 13, 1057–1098. [Google Scholar] [CrossRef] [Green Version]
  123. Davies, J. Origins and evolution of antibiotic resistance. Microbiologia 1996, 12, 9–16. [Google Scholar] [CrossRef] [Green Version]
  124. Aarts, H.; Margolles, A. Antibiotic resistance genes in food and gut (non pathogenic) bacteria. Bad genes in good bugs. Front. Microbiol. 2014, 5, 2014–2015. [Google Scholar] [CrossRef] [Green Version]
  125. Agersø, Y.; Bjerre, K.; Brockmann, E.; Johansen, E.; Nielsen, B.; Siezen, R.; Stuer-Lauridsen, B.; Wels, M.; Zeidan, A.A. Putative antibiotic resistance genes present in extant Bacillus licheniformis and Bacillus paralicheniformis strains are probably intrinsic and part of the ancient resistome. PLoS ONE 2019, 14, e0210363. [Google Scholar] [CrossRef]
  126. Deng, F.; Chen, Y.; Sun, T.; Wu, Y.; Su, Y.; Liu, C.; Zhou, J.; Deng, Y.; Wen, J. Antimicrobial resistance, virulence characteristics and genotypes of Bacillus spp. from probiotic products of diverse origins. Food Res. Int. 2021, 139, 109949. [Google Scholar] [CrossRef]
  127. Yelin, I.; Flett, K.B.; Merakou, C.; Mehrotra, P.; Stam, J.; Snesrud, E.; Hinkle, M.; Lesho, E.; McGann, P.; McAdam, A.J.; et al. Genomic and epidemiological evidence of bacterial transmission from probiotic capsule to blood in ICU patients. Nat. Med. 2019, 25, 1728–1732. [Google Scholar] [CrossRef] [PubMed]
Microorganisms 10 01038 g001 550
Figure 1. Maximum likelihood phylogenetic tree of BS50 and 20 other B. subtilis strains based on a concatenated sequence of the genes rpoB, purH, gyrA, groEL, polC, and 16S rRNA (15,093 nt). The bar indicates the rate of substitutions per nucleotide.
Figure 1. Maximum likelihood phylogenetic tree of BS50 and 20 other B. subtilis strains based on a concatenated sequence of the genes rpoB, purH, gyrA, groEL, polC, and 16S rRNA (15,093 nt). The bar indicates the rate of substitutions per nucleotide.
Microorganisms 10 01038 g001
Microorganisms 10 01038 g002 550
Figure 2. Effect of BS50 lysate treatment on Caco-2 cell viability after 48 h, as determined by ATP concentrations. Assay controls included untreated Caco-2 cells and cells that were fully lysed at the time of treatment. Data are expressed as mean ± standard deviation across technical replicates (n = 3).
Figure 2. Effect of BS50 lysate treatment on Caco-2 cell viability after 48 h, as determined by ATP concentrations. Assay controls included untreated Caco-2 cells and cells that were fully lysed at the time of treatment. Data are expressed as mean ± standard deviation across technical replicates (n = 3).
Microorganisms 10 01038 g002
Microorganisms 10 01038 g003 550
Figure 3. Effects of BS50 lysates on Caco-2 cell monolayer TEER in two separate trials (A,B). TEER was measured before treatment (0 h) and 2, 4, 6, 24, and 48 h after treatment. Square, untreated Caco-2 cells; diamond, “blank” lysate processing control; circle, BS50 lysate treatment; triangle, LPS treatment (TEER reduction control). Data are shown as two separate trials without replication within each trial (n = 1). Values on the y-axis are plotted on a logarithmic scale.
Figure 3. Effects of BS50 lysates on Caco-2 cell monolayer TEER in two separate trials (A,B). TEER was measured before treatment (0 h) and 2, 4, 6, 24, and 48 h after treatment. Square, untreated Caco-2 cells; diamond, “blank” lysate processing control; circle, BS50 lysate treatment; triangle, LPS treatment (TEER reduction control). Data are shown as two separate trials without replication within each trial (n = 1). Values on the y-axis are plotted on a logarithmic scale.
Microorganisms 10 01038 g003
Table
Table 1. Summary of BLASTn screening results for known Bacillus toxin genes in the BS50 genome.
Table 1. Summary of BLASTn screening results for known Bacillus toxin genes in the BS50 genome.
GeneOrganismAccessionMax Score% CoverageE-Value% Identity
gatAB. subtilis9387482405100%098%
metGB. cereus6157831391195%071%
HblAB. licheniformisKM514479.1No significant similarity
HblAB. cereusKF681259.135.612%0.02182%
HblCB. cereusJQ039142.1No significant similarity
HblDB. cereusJQ039158.1No significant similarity
NheA,B,CB. cereusDQ885236.142422%9 × 10−11870%
NheA,B,CB. mycoidesDQ153260.182.83%0.00268%
NheAB. cereusFN825684.1No significant similarity
NheA,B,CB. thuringiensisEU925144.1No significant similarity
entFMB. cereusAY789084.15914%9 × 10−0975%
cytKB. mycoidesAY871809.1No significant similarity
cytKB. licheniformisKM657965.1No significant similarity
cytKB. cereusDQ019311.137.41%0.04492%
HlyIIB. thuringiensis564444080No significant similarity
cesHPTABCDB. cereusNC_010924.110950%1 × 10−2279%
Table
Table 2. Summary of BLASTx screening results for known Bacillus toxin genes in the BS50 genome.
Table 2. Summary of BLASTx screening results for known Bacillus toxin genes in the BS50 genome.
ProteinOrganismAccessionMax ScoreE-Value% Identity
GatAB. subtilisNP_388550.18790100%
MetGB. cereusWP_079994147.1946074.16%
cytKB. mycoidesAAW56196.1No significant similarity found
EntFMB. cereusAAX14641.11214 × 10−2952.21%
cytKB. cereusAAY84864.1No significant similarity found
NheAB. mycoidesAAZ82480.1No significant similarity found
NheBB. mycoidesAAZ82481.1No significant similarity found
NheCB. mycoidesAAZ82482.1No significant similarity found
NheAB. cereusABI52601.1No significant similarity found
NheBB. cereusABI52602.1No significant similarity found
NheCB. cereusABI52603.1No significant similarity found
NheAB. cereusCBL95107.1No significant similarity found
NheA, partialB. thuringiensisACM18211.1No significant similarity found
NheBB. thuringiensisACM18212.1No significant similarity found
NheC, partialB. thuringiensisACM18213.1No significant similarity found
HblDB. cereusAFN08801.1No significant similarity found
HblCB. cereusAFN08807.1No significant similarity found
HblAB. cereusAII31101.1No significant similarity found
HblDB. licheniformisAIR07774.1No significant similarity found
HblAB. licheniformisAIR07775.1No significant similarity found
cytKB. licheniformisAIS75096.1No significant similarity found
CesAB. cereusWP_002081542.11250034.42%
CesBB. cereusWP_000953496.1776036.32%
CesCB. cereusWP_000590108.11446 × 10−3831.51%
CesDB. cereusWP_001008264.1No significant similarity found
CesHB. cereusWP_000291846.1532 × 10−0722.05%
CesPB. cereusWP_000680399.11293 × 10−3331.16%
CesTB. cereusWP_000764755.11164 × 10−2930.22%
Table
Table 3. Summary of secondary metabolite screening results for BS50 using antiSMASH.
Table 3. Summary of secondary metabolite screening results for BS50 using antiSMASH.
Cluster TypeMost Similar Cluster% Identity
NRPS (Non-ribosomal peptide synthases)Surfactin78%
NRPSFengycin100%
NRPSBacillibactin100%
OtherBacilysin100%
Polyketide + NRPBacillaene100%
RiPP: ThiopeptideSubtilosin A100%
RiPP: ThiopeptideSubtilomycin100%
CDPSN/AN/A
TerpeneN/AN/A
T3PKSN/AN/A
Table
Table 4. Summary of BS50 genome screening for virulence factors using VFDB.
Table 4. Summary of BS50 genome screening for virulence factors using VFDB.
GeneCategoryOrganismAccession% Ident% CoverageE
non-ribosomal peptide synthetase, DhbFBacillibactin;
Nutritional/
Metabolic factor		B. sub 168NP_39107699.1990
2,3-dihydroxybenzoate adenylase DhbEBacillibactin;
Nutritional/
Metabolic factor		B. sub 168NP_38972399.41000
isochorismate synthase DhbCBacillibactin;
Nutritional/
Metabolic factor		B. sub 168NP_39107898.51000
isochorismatase, DhbBBacillibactin;
Nutritional/
Metabolic factor		B. sub 168NP_39147199.71000
2,3-dihydroxybenzoate-2,3-
dehydrogenase, DhbA		Bacillibactin;
Nutritional/
Metabolic factor		B. sub 168NP_39107999.21000
gamma-glutamyltranspeptidase,
required for polyglutamate
anchoring to peptidoglycan		Capsule;
Immune
modulation		B. sub 168NP_39146998.91000
CapB, involved in
Poly-gamma-glutamate synthesis		Capsule;
Immune
modulation		B. sub 168NP_39107799.31000
CapA, required for
Poly-gamma-glutamate transport		Capsule;
Immune
modulation		B. sub 168NP_39108099.21000
CapC, involved in
Poly-gamma-glutamate synthesis		Capsule;
Immune
modulation		B. sub 168NP_3900621001000
endopeptidase Clp
ATP-binding chain C		ClpC;
Stress survival		B. sub 168NP_39147098.71000
(tufA) elongation factor TuEF-Tu;
Adherence		Lm EGD-eNP_46376372.6890
(hlyIII) putative membrane
hydrolase		Hemolysin III; ExotoxinFranc.WP_01392240674.7991.39 ×
10−142
Table
Table 5. Summary of antibiotic resistance genes detected in the BS50 genome using CARD.
Table 5. Summary of antibiotic resistance genes detected in the BS50 genome using CARD.
ARO Term (Gene)AMR Gene FamilyDrug Class% Identity% LengthRGI
Criteria
ykkDsmall multidrug resistance (SMR) antibiotic efflux pumpaminoglycoside antibiotic, tetracycline antibiotic, phenicol antibiotic100101.9Strict
lmrBATP-binding cassette (ABC) antibiotic efflux pumplincosamide antibiotic96.7100.42Strict
ykkCsmall multidrug resistance (SMR) antibiotic efflux pumpaminoglycoside antibiotic, tetracycline antibiotic, phenicol antibiotic100100Perfect
tet(45)major facilitator superfamily (MFS) antibiotic efflux pumptetracycline antibiotic75.8100Strict
mphKmacrolide phosphotransferase (MPH)macrolide antibiotic97.7100Strict
bltmajor facilitator superfamily (MFS) antibiotic efflux pumpfluoroquinolone antibiotic,
acridine dye		99.898.5Strict
Bacillus subtilis pgsA with mutation conferring resistance to daptomycindaptomycin resistant pgsApeptide antibiotic99.790.53Strict
Bacillus subtilis mprFdefensin resistant mprFpeptide antibiotic99.776.87Strict
vmlRABC-F ATP-binding cassette ribosomal protection proteinmacrolide antibiotic,
lincosamide antibiotic, streptogramin antibiotic, tetracycline antibiotic, oxazolidinone antibiotic, phenicol antibiotic, pleuromutilin antibiotic		98.575.5Strict
aadKANT(6)aminoglycoside antibiotic97.863.03Strict
bmrmajor facilitator superfamily (MFS) antibiotic efflux pumpfluoroquinolone antibiotic,
nucleoside antibiotic, acridine dye, phenicol antibiotic		10047.3Strict
tmrBtunicamycin resistance
protein		nucleoside antibiotic97.642.13Strict
aadKANT(6)aminoglycoside antibiotic97.239.44Strict
vmlRABC-F ATP-binding cassette ribosomal protection proteinmacrolide antibiotic,
lincosamide antibiotic, streptogramin antibiotic, tetracycline antibiotic, oxazolidinone antibiotic, phenicol antibiotic, pleuromutilin antibiotic		96.427.24Strict
tmrBtunicamycin resistance proteinnucleoside antibiotic10026.9Strict
Bacillus subtilis mprFdefensin resistant mprFpeptide antibiotic10016.36Strict
Table
Table 6. In vitro minimum inhibitory concentrations of antibiotics for BS50. The last column includes EFSA-recommended MIC thresholds for antibiotic resistance in Bacillus strains [94,95].
Table 6. In vitro minimum inhibitory concentrations of antibiotics for BS50. The last column includes EFSA-recommended MIC thresholds for antibiotic resistance in Bacillus strains [94,95].
AntibioticsTypeMIC (µg/mL)EFSA MIC (µg/mL)
Resistance Threshold
ChloramphenicolPhenicol28
ClindamycinMacrolides, lincosamides0.54
ErythromycinMacrolides, lincosamides<0.06254
GentamicinAminoglycosides0.54
KanamycinAminoglycoside28
StreptomycinAminoglycoside1258
OxytetracyclineTetracycline88
VancomycinGlycopeptide0.254
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Brutscher, L.M.; Borgmeier, C.; Garvey, S.M.; Spears, J.L. Preclinical Safety Assessment of Bacillus subtilis BS50 for Probiotic and Food Applications. Microorganisms 2022, 10, 1038. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms10051038

AMA Style

Brutscher LM, Borgmeier C, Garvey SM, Spears JL. Preclinical Safety Assessment of Bacillus subtilis BS50 for Probiotic and Food Applications. Microorganisms. 2022; 10(5):1038. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms10051038

Chicago/Turabian Style

Brutscher, Laura M., Claudia Borgmeier, Sean M. Garvey, and Jessica L. Spears. 2022. "Preclinical Safety Assessment of Bacillus subtilis BS50 for Probiotic and Food Applications" Microorganisms 10, no. 5: 1038. https://0-doi-org.brum.beds.ac.uk/10.3390/microorganisms10051038

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop