Next Article in Journal
Hepatitis B Virus Integration into Transcriptionally Active Loci and HBV-Associated Hepatocellular Carcinoma
Previous Article in Journal
Genome Characterization and Probiotic Potential of Corynebacterium amycolatum Human Vaginal Isolates
 
 
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 2
10.3390/microorganisms10020252
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

Antimicrobial Susceptibility Profile and Whole-Genome Analysis of a Strong Biofilm-Forming Bacillus Sp. B87 Strain Isolated from Food

by
Phornphan Sornchuer
1,2,*,
Kritsakorn Saninjuk
3,
Parisa Prathaphan
2,
Rattana Tiengtip
4 and
Suphot Wattanaphansak
3
1
Microbiology and Immunology, Department of Preclinical Science, Faculty of Medicine, Thammasat University, Klongluang 12120, Pathum Thani, Thailand
2
Thammasat University Research Unit in Nutraceuticals and Food Safety, Faculty of Medicine, Thammasat University, Klongluang 12120, Pathum Thani, Thailand
3
Departments of Veterinary Medicine, Faculty of Veterinary Science, Chulalongkorn University, Bangkok 10330, Thailand
4
Laboratory Section, Faculty of Medicine, Thammasat University, Klongluang 12120, Pathum Thani, Thailand
*
Author to whom correspondence should be addressed.
Microorganisms 2022, 10(2), 252; https://doi.org/10.3390/microorganisms10020252
Submission received: 12 December 2021 / Revised: 12 January 2022 / Accepted: 20 January 2022 / Published: 23 January 2022
(This article belongs to the Section Biofilm)
Browse Figures
Versions Notes

Abstract

:
Members of the Bacillus cereus group are considered to be foodborne pathogens commonly associated with diarrheal and emetic gastrointestinal syndromes. Biofilm formation is a major virulence determinant of various pathogenic bacteria, including the B. cereus strains, since it can protect the bacteria against antimicrobial agents and the host immune response. Moreover, a biofilm allows the exchange of genetic material, such as antimicrobial resistance genes, among the different bacterial strains inside the matrix. The aim of the current study was to genotypically and phenotypically characterize Bacillus sp. B87, a strain that was isolated from food and which exhibited strong biofilm-forming capacity. Based on the analysis of the phylogenetic relationship, the isolate was phylogenetically mapped close to Bacillus pacificus. Antimicrobial susceptibility testing revealed that the isolate was resistant to tetracycline and β-lactam antimicrobial agents, which corresponded with the genotypic characterization using the whole-genome analysis. The genome of Bacillus sp. B87 carried the three-component non-hemolytic enterotoxin (NHE), which is a type of enterotoxin that causes diarrheal symptoms. In addition, the genome also contained several genes that participate in biofilm formation, including the pelDEADAFG operon. These findings expand our understanding of antimicrobial resistance and virulence in Bacillus species based on the link between genotypic and phenotypic characterization.

1. Introduction

Bacillus cereus are Gram-positive, spore-forming bacteria that inhabit food, soil, and other natural media. B. cereus is a known human pathogen that causes food poisoning with emetic or diarrheal symptoms. Emetic strains of B. cereus can secrete the highly thermo- and pH-resistant toxin cereulide, which is pre-formed in food and causes vomiting a few hours after consumption [1]. Cereulide is an ionophoric dodecadepsipeptide that is produced by cereulide synthetase or Ces non-ribosomal peptide synthetase. For diarrheal strains of B. cereus, spores in contaminated food are consumed by the host, germinate within the small intestine, and the resulting vegetative cells then grow and produce enterotoxins. Three types of enterotoxins are secreted by B. cereus: the three-component enterotoxin hemolysin BL (HBL), the three-component non-hemolytic enterotoxin (NHE), and the single-component enterotoxin cytotoxin K (CytK) [2]. In addition to enterotoxins, B. cereus produces several other toxins and degradative enzymes, with most of these products controlled by the PlcR transcriptional activator [3]. PlcR is one of the B. cereus quorum-sensing systems that helps the bacterium adapt to diverse conditions [4].
B. cereus is typically resistant to β-lactam antimicrobial agents, such as penicillin G, ampicillin, and cefotaxime [5], due to the production of β-lactamase enzymes [6]. Bacterial resistance to commonly used antimicrobial agents, such as erythromycin, tetracycline, and streptomycin, can be a consequence of both nature and nurture [5,7]. Tetracycline is a broad-spectrum antimicrobial agent with activity against a wide range of bacteria, including Gram-positive and Gram-negative isolates. B. cereus is generally susceptible to tetracycline, but the resistance of B. cereus to this antimicrobial agent has been reported in some countries [8].
Biofilm formation is a major virulence determinant of various pathogenic bacteria, especially in the B. cereus group [9]. The formation of biofilms by bacteria can be associated with chronic infections in human and animal hosts. Moreover, biofilm formation allows the development and transfer of antimicrobial resistance through the bacterial interactions that occur within the biofilm [10,11]. The key genes for biofilm formation comprise those encoding biofilm transcriptional regulators, matrix structural genes, potential extracellular DNA synthesis genes, and cyclic-di-GMP metabolism genes [12]. In addition, several gene loci are involved in biofilm formation, including genes encoding the lipopeptide kurstakin, genes encoding the cyclic-di-GMP responsive effector protein BspA, and genes encoding the c-di-GMP synthesizing enzyme [13,14,15,16]. However, there may be strain-dependent variation in the mechanisms of biofilm formation among members of the B. cereus group.
The accessibility of whole-genome sequencing (WGS) has facilitated the assessment of bacterial genomes through bioinformatics analysis for the genetic potential to produce virulence factors and proteins involved in antimicrobial resistance and biofilm formation. In this study, the strong biofilm-forming strain, Bacillus sp. B87, isolated from food, was characterized genotypically and phenotypically, and genomic comparisons with other relevant B. cereus genomes were performed. The study aimed to generate insights into the genetic basis of antimicrobial resistance and virulence of this foodborne pathogen.

2. Materials and Methods

2.1. Bacterial Strains and Culture Conditions

The Bacillus sp. B87 used in this study was isolated from a spicy mussel salad in Pathum Thani province, Thailand [17]. Bacteria were aerobically grown in Luria–Bertani (LB) broth (Difco Laboratories, Detroit, MI, USA), with shaking at 180 rpm at 35 ± 2 °C.

2.2. Biofilm Formation Assay

The biofilm formation assay was based on a previously described protocol with minor modifications [18]. Briefly, Bacillus sp. B87 was grown overnight in LB broth at 35 °C and 180 rpm to generate inoculum cultures. Overnight cultures were adjusted to an optical density at 600 nm (OD600) of 0.01 in LB. Next, 100 µL of the adjusted bacterial culture was added to each well of a pre-sterilized 96-well flat-bottomed polystyrene microtiter plate, followed by incubation at 35 °C and 50 rpm for 24 h. Planktonic bacteria were then removed, and the wells were washed with distilled water and air-dried. Biofilm cells were stained with 200 µL of 0.3% crystal violet for 10 min, washed with distilled water, and air-dried. The crystal violet in the biofilm cells was solubilized with 200 µL of 70% ethanol, and the optical density at 590 nm (OD590) was measured.

2.3. Antimicrobial Susceptibility Tests

The antimicrobial susceptibility of Bacillus sp. B87 was determined using the Kirby–Bauer disk diffusion method according to standard criteria of the Clinical and Laboratory Standards Institute (CLSI) 2010 [19]. Briefly, the isolate was grown overnight at 35 °C on a nutrient agar (NA; Oxoid, United Kingdom) and the culture was compared with 0.5 McFarland turbidity standards. The culture was then applied onto Mueller–Hinton agar (MHA) plates using a sterile cotton swab, and the inoculated plates were dried at room temperature. The antimicrobial agents tested in this study included ampicillin (AMP, 10 µg), amoxicillin–clavulanic acid (AMC, 20 µg/10 mg), penicillin G (PEN, 10 U), gentamicin (GEN, 10 µg), imipenem (IPM, 10 µg), vancomycin (VAN, 30 µg), chloramphenicol (CHL, 30 µg), ciprofloxacin (CIP, 5 µg), tetracycline (TET, 30 µg), trimethoprim–sulfamethoxazole (SXT, 1.25 µg/23.75 µg), and erythromycin (ERY, 15 µg). Based on the zones of inhibition, Bacillus sp. B87 was classified as sensitive (S), intermediate (I), or resistant (R) to each antimicrobial agent according to the interpretative criteria for Staphylococcus spp., following CLSI guidelines [20]. Staphylococcus aureus ATCC 25923 was used as a control strain for the antimicrobial susceptibility tests.

2.4. Whole-Genome Sequencing, Assembly, and Annotation

Genomic DNA was extracted from Bacillus sp. B87 using a GF-1 Bacterial DNA Extraction Kit (Vivantis Technologies, Selangor, Malaysia) according to the manufacturer’s instructions. DNA quality was assessed via spectrophotometry and gel electrophoresis. Purified high molecular weight DNA samples with a 260/280 nm absorbance ratio of 1.8–2.0 and a 260/230 nm absorbance ratio of 2.0–2.2 were used for library construction and sequencing. The DNA sequencing library was prepared using a QIAGEN FX kit (Qiagen, Valencia, CA, USA), which fragments the gDNA using an enzymatic reaction, cleans the fragmented DNA with magnetic beads, and then ligates an adaptor index to the fragmented DNA. The quality and quantity of the indexed libraries were determined using an Agilent 2100 Bioanalyzer and a Denovix fluorometer, and the libraries were then pooled in equimolar quantities. Cluster generation and paired-end 2 × 150 nucleotide read sequencing were performed on an Illumina HiseqXten (Illumina Inc., San Diego, CA, USA).
The quality of the raw sequencing reads was assessed using FASTQC software. Adaptors and poor-quality reads were removed using Fastp, and the filtered reads were used as inputs for the Unicycler genome assembly program. The genome of Bacillus sp. B87 was annotated with the Rapid Annotation using Subsystem Technology tool kit (RASTtk) in PATRIC (Pathosystems Resource Integration Center). Sequences were queried using the BTyper tool, the Virulence Factor Database (VFDB) and Victors resource (for virulence factors), and the Comprehensive Antibiotic Resistance Database (CARD) and the NCBI National Database of Antibiotic Resistant Organisms (NDARO) (for antimicrobial resistance). In addition, the functional annotation of genes in terms of the Kyoto Encyclopedia of Genes and Genomes (KEGG) orthology assignments and predictions of KEGG pathways were performed through the KEGG Automatic Annotation Server (KAAS; https://www.genome.jp/kegg/kaas/ accessed on 18 October 2021) using the bi-directional best hit (BBH) method [21]. The circular genome map was constructed by using the circular viewer of PATRIC. A comparison of syntenic analyses in Bacillus sp. B87 and other species of the genus Bacillus was performed by using the Easyfig program [22].

2.5. Phylogenetic Analysis

A phylogenomic tree based on the core genes was generated with PATRIC Phylogenetic Tree Building Service [23]. The default was set for codon trees, which utilizes both the protein and gene sequences from PATRIC’s global protein families (PGFams). Protein sequences were aligned using MUSCLE, and the nucleotide coding gene sequences were aligned using the Codon_align function of BioPython. A concatenated alignment of all proteins and nucleotides was converted to a phylip formatted file, and then a partitions file for RaxML was constructed. Support values were created using 100 rounds of the “Rapid” bootstrapping option of RaxML. The phylogenomic classification of Bacillus sp. B87 was also performed using the Type (strain) Genome Server (TYGS), a free bioinformatics platform for a whole-genome-based taxonomic analysis [24].

2.6. Nucleotide Sequence Accession Numbers

The draft genomes of Bacillus sp. B87 were deposited in the NCBI database under the accession numbers SRR16129018. The BioProject ID in GenBank is PRJNA7260213.

3. Results

3.1. Biofilm-Forming Ability of Bacillus Sp. B87

The ability of Bacillus sp. B87 to produce biofilm was determined in comparison with B. cereus ATCC 14579. A microtiter plate assay for biofilm formation revealed that the biofilm level of Bacillus sp. B87 (OD590 0.545 ± 0.149) after incubation for 24 h was significantly higher than that of B. cereus ATCC 14579 (OD590 0.035 ± 0.002) (Figure 1).

3.2. Antimicrobial Resistance Profile of Bacillus Sp. B87

Bacillus sp. B87 was tested for susceptibility to 11 selected antimicrobial agents, as shown in Table 1. Bacillus sp. B87 was susceptible to most of the tested antimicrobial agents, including GEN (24.8 mm. ± 1.4), IPM (39.4 mm. ± 0.3), VAN (22.8 mm. ± 0.5), CHL (30.7 mm. ± 0.7), CIP (33.2 mm. ± 1.7), SXT (21.0 mm. ± 1.6), and ERY (28.8 mm. ± 0.6). The isolate was resistant to antimicrobial agents in the β-lactam category, including AMP (14.3 mm. ± 0.8), AMC (8.2 mm. ± 0.7), and PEN (10.6 mm. ± 2.3), and was also resistant to tetracycline (12.5 mm. ± 1.4).

3.3. Genetic Features of Bacillus Sp. B87

Genomic features and annotation information for the genome of Bacillus sp. B87 are summarized in Table 2. The draft genome sequence had an estimated length of 5,448,163 bp, a GC content of 35.18%, and contained 5661 coding sequences. The circular representation of the Bacillus sp. B87 draft genome was generated using the circular viewer of PATRIC and is shown in Figure 2. Potential genes in the draft genome of Bacillus sp. B87 were investigated and annotated based on different biological processes and metabolic pathways using the RAST server (Figure 3). The predicted genes included 711 genes involved in metabolism, 272 genes involved in cellular processes, and 142 genes involved in stress response, defense, and virulence. Genes involved in prophages, transposable elements, and plasmids were also found in the draft genome of Bacillus sp. B87, and were classified as: subclass: pathogenicity islands; subsystem name: Listeria Pathogenicity Island LIPI-1 extended. These pathogenicity islands were also present in the genomes of B. cereus ATCC 14579 and B. anthracis str. Ames.
A phylogenetic tree based on core genes was reconstructed in PATRIC using the whole-genome sequence of Bacillus sp. B87 (Figure 4A). This phylogenomic analysis revealed that Bacillus sp. B87 was closely related to strains of Bacillus pacificus. The 16S rRNA gene (Figure 4B) and whole-genome (Figure 4C) phylogeny reconstructions using TYGS confirmed a close association of Bacillus sp. B87 with B. pacificus.

3.4. Antimicrobial Resistance Genes

Antimicrobial resistance (AMR) genes that are associated with resistance to one or more antimicrobial agents were predicted based on the CARD and NDARO databases, and Bacillus sp. B87 contained nine genes connected with resistance to different antimicrobial agents (Table 3). Moreover, 43 AMR genes were annotated according to the PATRIC database using K-mer Search (Supplementary Table S1). These findings might be associated with the observed antimicrobial resistance phenotypes to β-lactam antimicrobial agents and tetracycline shown in Table 1.

3.5. Biofilm Formation Genes

According to RASTtk, available in PATRIC, the genes involved in the formation of biofilm by Bacillus sp. B87 were identified. Nine genes were detected (Table 4) according to the following hierarchical classification: superclass: cellular process; class: microbial communities; subclass: quorum sensing and biofilm formation. Key biofilm-formation genes were also investigated using the KAAS database, and this analysis identified the genes encoding biofilm transcriptional regulators, matrix protein-encoding genes, putative matrix polysaccharide synthesis genes, and extracellular DNA (eDNA) synthesis genes (Table 5).
The differences in gene synteny between the genomes of Bacillus sp. B87 and B. cereus ATCC 14579 were apparent with the pelDEADAFG operon (Figure 5). Bacillus sp. B87 contained pelA, pelD, pelF, and pelG, and was therefore similar to B. cereus ATCC 10987 and Bacillus sp. EB422 (B. pacificus).

3.6. Virulence Factor Genes

Thirteen genes in the draft genome of Bacillus sp. B87 were classified as virulence factors according to the Victors (nine genes) and VFDB (four genes) databases (Table 6). Among these genes were nheA, nheB, and nheC, which encode the three-component NHE complex, a type of enterotoxin that causes diarrheal symptoms. However, the cytotoxin K virulence factor, which is also associated with diarrheal illness and was present in B. cereus ATCC 14579, was not detected in Bacillus sp. B87 (Supplementary Table S2). Virulence genes were also predicted by using the BTyper tool, and this analysis predicted 17 virulence genes in the genome of Bacillus sp. B87 (Table 6).

4. Discussion

This study aimed to genotypically and phenotypically characterize Bacillus sp. B87, a strain with strong biofilm-forming activity that was previously isolated from a spicy mussel salad in Pathum Thani province, Thailand. The phylogenetic analysis revealed that Bacillus sp. B87 was closely related to B. pacificus, which had been reported to be isolated from the sediment of the Pacific Ocean [25]. Currently, little is known about the phenotype and genotype of this bacterial strain.
The alignment of the Bacillus sp. B87 genome indicated that the isolate possesses the genes nheA, nheB, and nheC, but not hblA, hblC, hblD, or cytK. In B. cereus, the toxins that are associated with diarrheal diseases are HBL, NHE, CytK, and enterotoxin FM [26,27,28,29,30,31,32]. Indeed, isolates that carry the genes encoding HBL might be more virulent [33]. However, a strain lacking the HBL operon—B. cereus ATCC 10987—was reported to exhibit strong cytopathogenic activity in Vero cells, since it could produce a large amount of the NHE mRNA [34]. The genome of Bacillus sp. B87 also contains Listeria Pathogenicity Island LIPI-1 extended. LIP-1 in Listeria monocytogenes harbors several important genes, including prfA, plcA, hly, mpl, actA, and plcB, which are involved in host invasion and cellular proliferation [35]. Orthologs for LIP-1 were also detected in B. cereus ATCC 14579 and B. cereus strain FORC_021, which was isolated from a knife used at a sashimi restaurant in the Republic of Korea [36]. However, the coding sequence (CDS) counts of LIPI-1 genes in each strain were dissimilar in that report. The presence of Listeria Pathogenicity Island LIPI-1 extended in Bacillus sp. B87 might be responsible for the virulence of this strain and requires further investigation.
B. cereus is typically resistant to β-lactam antimicrobial agents, including penicillin G, ampicillin, and amoxicillin–clavulanic acid [17,20,37,38]. Some bacterial strains can produce β-lactamase enzymes that are responsible for the resistance to β-lactam antimicrobial agents. Three different β-lactamases have been classified in strains of B. cereus, namely, β-lactamase I, II, and III [39]. B. cereus β-lactamase II (BcII) is a heat-stable metallo-β-lactamase (MBL) that shares high sequence homology with Bla2 from B. anthracis [40]. MBL catalyzes the hydrolysis of β-lactam antimicrobial agents, including penicillin, cephalosporin, carbapenem, and cephamycin [41]. However, BcII catalyzes the hydrolysis of penicillin at higher rates than cephalosporin and carbapenem [42]. The current study revealed the presence of the BcII gene in the draft genome of Bacillus sp. B87, and this gene might play an important role in the β-lactam resistance phenotype of the isolate.
Bacillus sp. B87 carries the tet(45) gene, which contributes to the tetracycline resistance phenotype in this isolate. Tet45 is a tetracycline efflux pump closely related to TetL [43]. Homologs of Tet45 have been found in the genomes of strains of B. cereus and B. thuringiensis. Furthermore, tetA-carrying B. cereus was reported to be susceptible to tetracycline [44]. Moreover, some isolates of B. cereus were phenotypically resistant to tetracycline, even though they did not carry tetA, tetB, or tetC, and this might be due to the presence of other tetracycline resistance genes, such as tetM and tetL [45]. B. cereus is generally susceptible to tetracycline, but the resistance of B. cereus to tetracycline has been reported in some countries [8]. In the current study, Bacillus sp. B87 was resistant to tetracycline, since it possessed the gene encoding the tetracycline resistance efflux pump. This phenotype is concordant with a previous report on B. cereus strain MS532a, which presented a tetracycline-resistant phenotype and carried the tet(45) tetracycline resistance gene [38]. Horizontal gene transfer may have contributed to the dissemination and persistence of tet(45) in environments such as a poultry litter-impacted soil [43]. The presence of the tet(45) gene in the genome of Bacillus sp. B87, a strain isolated from spicy mussel salad, may be due to the contamination of vegetables, mussels, or other ingredients. Tetracycline-resistant isolates of B. cereus have previously been shown to carry the tet(L) gene on a plasmid, while other species of the genus Bacillus carry either tet(L) or tet(K) on plasmids and/or in the chromosome [46,47,48]. Both genes—tet(L) and tet(K)—can occasionally be mobilized in the presence of conjugative plasmids, but are not themselves able to independently transfer, thus decelerating their spread within the population [45]. The ability to form biofilm in bacteria is an important mechanism that allows the bacteria to resist antimicrobial agents and disinfectants, as well as to evade the host immune system. B. cereus can form biofilms on various surfaces [49,50,51], including plastic, glass wool, and stainless steel, and the biofilm cells were more resistant to sanitizers compared with planktonic cells [52]. Biofilm plays an important role in the exchange of antimicrobial resistance genes [53]. In strains of Escherichia coli and Pseudomonas aeruginosa, sub-inhibitory concentrations of tetracycline and cephradine induce biofilm formation and increase the transfer rate of the pB10 plasmid among the biofilm biomass at higher rates compared with no antimicrobial treatment [54]. The findings from the current study were congruent with a previous report [55] that resistant bacteria could form biofilms. Tetracycline-resistant Bacillus sp. B87, with strong biofilm-formation ability, may develop greater resistance to different antimicrobial agents and other environmental stressors if left unmonitored.
The regulation network controlling B. cereus biofilm formation involves several pathways. PlcR is a pleiotropic regulator that controls the expression of genes encoding several enterotoxins, hemolysins, phospholipases, and proteases in B. cereus [9,56]. The plcR gene is instrumental in biofilm formation, since the deletion of this gene in B. cereus ATCC 14579 resulted in a significant increase in the amount of biofilm [57]. Two roles have been reported for the regulator CodY: the repression of biofilm formation in B. cereus ATCC 14579 [58] and the promotion of biofilm formation in B. cereus UW101C [59]. The codY gene operon has a role in pellicle biofilm formation and swarming motility [60]. An operon including tasA, tapA, and sipW has been proposed to be involved in biofilm formation in B. subtilis [61]. The transcription of these genes is regulated by SinI and repressed by SinR [62]. The transcription of sinI is activated by the master regulator of sporulation Spo0A [63]. In the present study, Bacillus sp. B87 was found to carry the genes tasA and sipW, which may be sufficient for biofilm production. The deletion of the genomic region encoding two orthologs of the amyloid-like protein TasA and SipW signal peptidase inhibited biofilm assembly [64]. However, the genes calY and tapA were not detected in the isolate. Mutations in tasA or calY did not completely prevent biofilm formation. Moreover, the lack of tapA may not completely interrupt biofilm production, which was congruent with a previous report [65].
Several gene loci are involved in biofilm formation, and there may be strain-to-strain variability in matrix component utilization among isolates of B. cereus. For example, the pelDEADAFG operon identified in strain ATCC 10987 is a crucial locus for biofilm formation in B. cereus. However, this locus is not present in strain ATCC 14579 of this species. This operon might be required for biofilm formation in Bacillus cereus strains, at least on polystyrene surfaces, since ATCC 10987 formed biofilms, while ATCC 14579 did not form biofilms on polystyrene 96-well plates [66]. In addition, the deletion of any of pelDEADAFG in ATCC 10987 resulted in the reduction of bacterial adhesion to the wells of the plastic microtiter dish [67]. The bioinformatics analysis and validation of Pel production in Bacillus sp. B87 in the current study suggested that the pelDEADAFG operon was present and potentially functional in biofilm-forming strains. In addition, the presence of the gene recA in the genome of Bacillus sp. B87 may enhance the biofilm-forming capacity of this isolate, since it has been reported that the biofilm formation and swarming motility of B. cereus 905 are promoted by RecA [68]. However, further investigation via a gene deletion approach is required to confirm this theory.

5. Conclusions

This study reported the phenotypic and genotypic characterization of Bacillus sp. B87, an isolate with strong biofilm-forming capacity. The isolate was resistant to β-lactam antimicrobial agents and tetracycline, since it carried BcII and tet(45) genes in the genome. Consequently, the strain can survive in tetracycline-containing environments where it may induce the biofilm-forming capacity of this isolate, and this assumption needs to be further elucidated. Moreover, the presence of the pelDEADAFG operon in the genome of Bacillus sp. B87 might play an important role in the biofilm-forming capacity of this isolate. Biofilm formation would allow the bacteria to exchange genetic material among strains, and could potentially lead to the development of greater resistance to different antimicrobial agents and other environmental stressors. In conclusion, the identification of genes encoding virulence factors and antimicrobial resistance in foodborne bacteria should be considered to be potential key points to assess human health risks from the bacteria. Moreover, the findings from this study suggest that WGS analysis could be an effective tool to elucidate the pathogenic potential of the B. cereus group. However, the molecular basis proposed in this study needs to be further clarified through gene knockout and protein characterization, since several gene homologs in members of the B. cereus group have unique and varied functions.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/microorganisms10020252/s1, Table S1: PATRIC database using K-mer Search, Table S2: PATRIC virulence factors.

Author Contributions

Conceptualization, P.S. and K.S.; methodology, P.S., R.T. and P.P.; software, K.S.; validation, P.S., K.S. and S.W.; formal analysis, P.S.; writing—original draft preparation, P.S.; writing—review and editing, P.S., K.S., P.P., R.T. and S.W.; visualization, P.S.; supervision, K.S.; project administration, P.S.; funding acquisition, P.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the Faculty of Medicine, Thammasat University, Thailand. The grant number is 2-16/2563. This work was also supported by the Research Group in Multidrug Resistant Bacteria and the Antimicrobial Herbal Extracts from Faculty of Medicine, Thammasat University, and Thammasat University Research Unit in Nutraceuticals and Food Safety.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw sequence generated for this study was deposited into the National Centre for Biotechnology Information (NCBI) Sequence Read Archive (SRA) under the accession numbers SRR16129018, BioProject ID: PRJNA726021. All the other data supporting the finding of this study are available in this published article and its supplementary information files.

Acknowledgments

The authors would like to acknowledge Sumet Amonyingcharoen for his assistance with preparing the figures.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ehling-Schulz, M.; Frenzel, E.; Gohar, M. Food-bacteria interplay: Pathometabolism of emetic Bacillus cereus. Front. Microbiol. 2015, 6, 704. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Stenfors Arnesen, L.P.; Fagerlund, A.; Granum, P.E. From soil to gut: Bacillus cereus and its food poisoning toxins. FEMS Microbiol. Rev. 2008, 32, 579–606. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Gohar, M.; Faegri, K.; Perchat, S.; Ravnum, S.; Okstad, O.A.; Gominet, M.; Kolsto, A.B.; Lereclus, D. The PlcR virulence regulon of Bacillus cereus. PLoS ONE 2008, 3, e2793. [Google Scholar] [CrossRef] [PubMed]
  4. de Been, M.; Francke, C.; Moezelaar, R.; Abee, T.; Siezen, R.J. Comparative analysis of two-component signal transduction systems of Bacillus cereus, Bacillus thuringiensis and Bacillus anthracis. Microbiology 2006, 152 Pt 10, 3035–3048. [Google Scholar] [CrossRef] [Green Version]
  5. Citron, D.M.; Appleman, M.D. In Vitro activities of daptomycin, ciprofloxacin, and other antimicrobial agents against the cells and spores of clinical isolates of Bacillus species. J. Clin. Microbiol. 2006, 44, 3814–3818. [Google Scholar] [CrossRef] [Green Version]
  6. Park, Y.B.; Kim, J.B.; Shin, S.W.; Kim, J.C.; Cho, S.H.; Lee, B.K.; Ahn, J.; Kim, J.M.; Oh, D.H. Prevalence, genetic diversity, and antibiotic susceptibility of Bacillus cereus strains isolated from rice and cereals collected in Korea. J. Food. Prot. 2009, 72, 612–617. [Google Scholar] [CrossRef]
  7. Jensen, L.B.; Baloda, S.; Boye, M.; Aarestrup, F.M. Antimicrobial resistance among Pseudomonas spp. and the Bacillus cereus group isolated from Danish agricultural soil. Environ. Int. 2001, 26, 581–587. [Google Scholar] [CrossRef]
  8. Luna, V.A.; King, D.S.; Gulledge, J.; Cannons, A.C.; Amuso, P.T.; Cattani, J. Susceptibility of Bacillus anthracis, Bacillus cereus, Bacillus mycoides, Bacillus pseudomycoides and Bacillus thuringiensis to 24 antimicrobials using Sensititre automated microbroth dilution and Etest agar gradient diffusion methods. J. Antimicrob. Chemother. 2007, 60, 555–567. [Google Scholar] [CrossRef] [Green Version]
  9. Rutherford, S.T.; Bassler, B.L. Bacterial quorum sensing: Its role in virulence and possibilities for its control. Cold Spring Harb Perspect. Med. 2012, 2, a012427. [Google Scholar] [CrossRef]
  10. Gomes, F.; Saavedra, M.J.; Henriques, M. Bovine mastitis disease/pathogenicity: Evidence of the potential role of microbial biofilms. Pathog. Dis. 2016, 74, ftw006. [Google Scholar] [CrossRef] [Green Version]
  11. McCarthy, H.; Rudkin, J.K.; Black, N.S.; Gallagher, L.; O’Neill, E.; O’Gara, J.P. Methicillin resistance and the biofilm phenotype in Staphylococcus aureus. Front. Cell Infect. Microbiol. 2015, 5, 1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Ikram, S.; Heikal, A.; Finke, S.; Hofgaard, A.; Rehman, Y.; Sabri, A.N.; Okstad, O.A. Bacillus cereus biofilm formation on central venous catheters of hospitalised cardiac patients. Biofouling 2019, 35, 204–216. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Fagerlund, A.; Smith, V.; Rohr, A.K.; Lindback, T.; Parmer, M.P.; Andersson, K.K.; Reubsaet, L.; Okstad, O.A. Cyclic diguanylate regulation of Bacillus cereus group biofilm formation. Mol. Microbiol. 2016, 101, 471–494. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Gelis-Jeanvoine, S.; Canette, A.; Gohar, M.; Caradec, T.; Lemy, C.; Gominet, M.; Jacques, P.; Lereclus, D.; Slamti, L. Genetic and functional analyses of krs, a locus encoding kurstakin, a lipopeptide produced by Bacillus thuringiensis. Res. Microbiol. 2017, 168, 356–368. [Google Scholar] [CrossRef] [PubMed]
  15. Sudarsan, N.; Lee, E.R.; Weinberg, Z.; Moy, R.H.; Kim, J.N.; Link, K.H.; Breaker, R.R. Riboswitches in eubacteria sense the second messenger cyclic di-GMP. Science 2008, 321, 411–413. [Google Scholar] [CrossRef] [Green Version]
  16. Tang, Q.; Yin, K.; Qian, H.; Zhao, Y.; Wang, W.; Chou, S.H.; Fu, Y.; He, J. Cyclic di-GMP contributes to adaption and virulence of Bacillus thuringiensis through a riboswitch-regulated collagen adhesion protein. Sci. Rep. 2016, 6, 28807. [Google Scholar] [CrossRef]
  17. Sornchuer, P.; Tiengtip, R. Prevalence, virulence genes, and antimicrobial resistance of Bacillus cereus isolated from foodstuffs in Pathum Thani Province, Thailand. Pharm. Sci. Asia 2021, 48, 194–203. [Google Scholar] [CrossRef]
  18. Hsueh, Y.H.; Somers, E.B.; Lereclus, D.; Wong, A.C. Biofilm formation by Bacillus cereus is influenced by PlcR, a pleiotropic regulator. Appl. Environ. Microbiol. 2006, 72, 5089–5092. [Google Scholar] [CrossRef] [Green Version]
  19. Clinical and Laboratory Standards Institute. Performance Standards for Antimicrobial Susceptibility Testing, 20th Informational Supplement; CLSI Document M100-S20; CLSI: Wayne, PA, USA, 2010. [Google Scholar]
  20. Gao, T.; Ding, Y.; Wu, Q.; Wang, J.; Zhang, J.; Yu, S.; Yu, P.; Liu, C.; Kong, L.; Feng, Z.; et al. Prevalence, virulence genes, antimicrobial susceptibility, and genetic diversity of Bacillus cereus isolated from pasteurized milk in China. Front. Microbiol. 2018, 9, 533. [Google Scholar] [CrossRef] [Green Version]
  21. Moriya, Y.; Itoh, M.; Okuda, S.; Yoshizawa, A.C.; Kanehisa, M. KAAS: An automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007, 35, W182-5. [Google Scholar] [CrossRef] [Green Version]
  22. Sullivan, M.J.; Petty, N.K.; Beatson, S.A. Easyfig: A genome comparison visualizer. Bioinformatics 2011, 27, 1009–1010. [Google Scholar] [CrossRef] [PubMed]
  23. Fiedler, G.; Herbstmann, A.D.; Doll, E.; Wenning, M.; Brinks, E.; Kabisch, J.; Breitenwieser, F.; Lappann, M.; Bohnlein, C.; Franz, C. Taxonomic evaluation of the Heyndrickxia (Basonym Bacillus) sporothermodurans Group (H. sporothermodurans, H. vini, H. oleronia) based on whole genome sequences. Microorganisms 2021, 9, 246. [Google Scholar] [CrossRef] [PubMed]
  24. Meier-Kolthoff, J.P.; Goker, M. TYGS is an automated high-throughput platform for state-of-the-art genome-based taxonomy. Nat. Commun. 2019, 10, 2182. [Google Scholar] [CrossRef] [PubMed]
  25. Liu, Y.; Du, J.; Lai, Q.; Zeng, R.; Ye, D.; Xu, J.; Shao, Z. Proposal of nine novel species of the Bacillus cereus group. Int. J. Syst. Evol. Microbiol. 2017, 67, 2499–2508. [Google Scholar] [CrossRef] [PubMed]
  26. Granum, P.E.; O’Sullivan, K.; Lund, T. The sequence of the non-haemolytic enterotoxin operon from Bacillus cereus. FEMS Microbiol. Lett. 1999, 177, 225–229. [Google Scholar] [CrossRef] [Green Version]
  27. Beecher, D.J.; Wong, A.C.L. Tripartite haemolysin BL: Isolation and characterization of two distinct homologous sets of components from a single Bacillus cereus isolate. Microbiology 2000, 146 Pt 6, 1371–1380. [Google Scholar] [CrossRef] [Green Version]
  28. Lund, T.; De Buyser, M.L.; Granum, P.E. A new cytotoxin from Bacillus cereus that may cause necrotic enteritis. Mol. Microbiol. 2000, 38, 254–261. [Google Scholar] [CrossRef]
  29. Hansen, B.M.; Hendriksen, N.B. Detection of enterotoxic Bacillus cereus and Bacillus thuringiensis strains by PCR analysis. Appl. Environ. Microbiol. 2001, 67, 185–189. [Google Scholar] [CrossRef] [Green Version]
  30. Fagerlund, A.; Ween, O.; Lund, T.; Hardy, S.P.; Granum, P.E. Genetic and functional analysis of the cytK family of genes in Bacillus cereus. Microbiology 2004, 150 Pt 8, 2689–2697. [Google Scholar] [CrossRef] [Green Version]
  31. Sergeev, N.; Distler, M.; Vargas, M.; Chizhikov, V.; Herold, K.E.; Rasooly, A. Microarray analysis of Bacillus cereus group virulence factors. J. Microbiol. Methods 2006, 65, 488–502. [Google Scholar] [CrossRef]
  32. Hernandez, A.G.C.; Ortiz, V.G.; Gomez, J.L.A.; Lopez, M.A.R.; Morales, J.A.R.; Macias, A.F.; Hidalgo, E.A.; Ramirez, J.N.; Gallardo, F.J.F.; Gutierrez, M.C.G.; et al. Detection of Bacillus cereus sensu lato isolates posing potential health risks in Mexican chili powder. Microorganisms 2021, 9, 2226. [Google Scholar] [CrossRef] [PubMed]
  33. Salamitou, S.; Ramisse, F.; Brehelin, M.; Bourguet, D.; Gilois, N.; Gominet, M.; Hernandez, E.; Lereclus, D. The plcR regulon is involved in the opportunistic properties of Bacillus thuringiensis and Bacillus cereus in mice and insects. Microbiology 2000, 146 Pt 11, 2825–2832. [Google Scholar] [CrossRef] [Green Version]
  34. Lindback, T.; Okstad, O.A.; Rishovd, A.L.; Kolsto, A.B. Insertional inactivation of hblC encoding the L2 component of Bacillus cereus ATCC 14579 haemolysin BL strongly reduces enterotoxigenic activity, but not the haemolytic activity against human erythrocytes. Microbiology 1999, 145 Pt 11, 3139–3146. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Chen, Y.; Luo, Y.; Curry, P.; Timme, R.; Melka, D.; Doyle, M.; Parish, M.; Hammack, T.S.; Allard, M.W.; Brown, E.W.; et al. Assessing the genome level diversity of Listeria monocytogenes from contaminated ice cream and environmental samples linked to a listeriosis outbreak in the United States. PLoS ONE 2017, 12, e0171389. [Google Scholar] [CrossRef] [PubMed]
  36. Chung, H.Y.; Lee, K.H.; Ryu, S.; Yoon, H.; Lee, J.H.; Kim, H.B.; Kim, H.; Jeong, H.G.; Choi, S.H.; Kim, B.S. Genome Sequence of Bacillus cereus FORC_021, a Food-Borne Pathogen Isolated from a Knife at a Sashimi Restaurant. J. Microbiol. Biotechnol. 2016, 26, 2030–2035. [Google Scholar] [CrossRef] [Green Version]
  37. Owusu-Kwarteng, J.; Wuni, A.; Akabanda, F.; Tano-Debrah, K.; Jespersen, L. Prevalence, virulence factor genes and antibiotic resistance of Bacillus cereus sensu lato isolated from dairy farms and traditional dairy products. BMC Microbiol. 2017, 17, 65. [Google Scholar] [CrossRef] [Green Version]
  38. Fiedler, G.; Schneider, C.; Igbinosa, E.O.; Kabisch, J.; Brinks, E.; Becker, B.; Stoll, D.A.; Cho, G.S.; Huch, M.; Franz, C. Antibiotics resistance and toxin profiles of Bacillus cereus-group isolates from fresh vegetables from German retail markets. BMC Microbiol. 2019, 19, 250. [Google Scholar] [CrossRef]
  39. Chen, Y.; Succi, J.; Tenover, F.C.; Koehler, T.M. Beta-lactamase genes of the penicillin-susceptible Bacillus anthracis Sterne strain. J. Bacteriol. 2003, 185, 823–830. [Google Scholar] [CrossRef] [Green Version]
  40. Hawk, M.J.; Breece, R.M.; Hajdin, C.E.; Bender, K.M.; Hu, Z.; Costello, A.L.; Bennett, B.; Tierney, D.L.; Crowder, M.W. Differential binding of Co(II) and Zn(II) to metallo-beta-lactamase Bla2 from Bacillus anthracis. J. Am. Chem. Soc. 2009, 131, 10753–10762. [Google Scholar] [CrossRef] [Green Version]
  41. Page, M.I.; Badarau, A. The mechanisms of catalysis by metallo beta-lactamases. Bioinorg. Chem. Appl. 2008, 2008, 576297. [Google Scholar] [CrossRef] [Green Version]
  42. Felici, A.; Amicosante, G.; Oratore, A.; Strom, R.; Ledent, P.; Joris, B.; Fanuel, L.; Frere, J.M. An overview of the kinetic parameters of class B beta-lactamases. Biochem. J. 1993, 291 Pt 1, 151–155. [Google Scholar] [CrossRef] [PubMed]
  43. You, Y.; Hilpert, M.; Ward, M.J. Identification of Tet45, a tetracycline efflux pump, from a poultry-litter-exposed soil isolate and persistence of tet(45) in the soil. J. Antimicrob. Chemother. 2013, 68, 1962–1969. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Rather, M.; Aulakh, R.; Munshi, Z.H.; Hussain, S.; Shah, s.; Wani, N. Antibiogram and detection of tetracycline resistance genes (tetA, tetB) of Bacillus cereus from foods of animal origin. SKUAST J. Res. 2015, 17, 99–105. [Google Scholar]
  45. Agerso, Y.; Jensen, L.B.; Givskov, M.; Roberts, M.C. The identification of a tetracycline resistance gene tet(M), on a Tn916-like transposon, in the Bacillus cereus group. FEMS Microbiol. Lett. 2002, 214, 251–256. [Google Scholar] [CrossRef]
  46. Bernhard, K.; Schrempf, H.; Goebel, W. Bacteriocin and antibiotic resistance plasmids in Bacillus cereus and Bacillus subtilis. J. Bacteriol. 1978, 133, 897–903. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Eccles, S.J.; Chopra, I. Biochemical and genetic characterization of the tet determinant of Bacillus plasmid pAB124. J. Bacteriol. 1984, 158, 134–140. [Google Scholar] [CrossRef] [Green Version]
  48. Ives, C.L.; Bott, K.F. Cloned Bacillus subtilis chromosomal DNA mediates tetracycline resistance when present in multiple copies. J. Bacteriol. 1989, 171, 1801–1810. [Google Scholar] [CrossRef] [Green Version]
  49. Auger, S.; Krin, E.; Aymerich, S.; Gohar, M. Autoinducer 2 affects biofilm formation by Bacillus cereus. Appl. Environ. Microbiol. 2006, 72, 937–941. [Google Scholar] [CrossRef] [Green Version]
  50. Oosthuizen, M.C.; Steyn, B.; Theron, J.; Cosette, P.; Lindsay, D.; Von Holy, A.; Brozel, V.S. Proteomic analysis reveals differential protein expression by Bacillus cereus during biofilm formation. Appl. Environ. Microbiol. 2002, 68, 2770–2780. [Google Scholar] [CrossRef] [Green Version]
  51. Peng, J.S.; Tsai, W.C.; Chou, C.C. Surface characteristics of Bacillus cereus and its adhesion to stainless steel. Int. J. Food Microbiol. 2001, 65, 105–111. [Google Scholar] [CrossRef]
  52. Peng, J.S.; Tsai, W.C.; Chou, C.C. Inactivation and removal of Bacillus cereus by sanitizer and detergent. Int. J. Food Microbiol. 2002, 77, 11–18. [Google Scholar] [CrossRef]
  53. Abe, K.; Nomura, N.; Suzuki, S. Biofilms: Hot spots of horizontal gene transfer (HGT) in aquatic environments, with a focus on a new HGT mechanism. FEMS Microbiol. Ecol. 2020, 96, fiaa031. [Google Scholar] [CrossRef] [PubMed]
  54. Salcedo, D.; Lee, J.; Ha, U.-H.; Kim, S. The effects of antibiotics on the biofilm formation and antibiotic resistance gene transfer. Desalination Water Treat. 2014, 54, 3582–3588. [Google Scholar] [CrossRef]
  55. Stachurová, T.; Malachová, K.; Semerád, J.; Sterniša, M.; Rybková, Z.; Smole Možina, S. Tetracycline induces the formation of biofilm of bacteria from different phases of wastewater treatment. Processes 2020, 8, 989. [Google Scholar] [CrossRef]
  56. Slamti, L.; Lereclus, D. A cell-cell signaling peptide activates the PlcR virulence regulon in bacteria of the Bacillus cereus group. EMBO J. 2002, 21, 4550–4559. [Google Scholar] [CrossRef] [Green Version]
  57. Ryu, J.H.; Beuchat, L.R. Biofilm formation and sporulation by Bacillus cereus on a stainless steel surface and subsequent resistance of vegetative cells and spores to chlorine, chlorine dioxide, and a peroxyacetic acid-based sanitizer. J. Food Prot. 2005, 68, 2614–2622. [Google Scholar] [CrossRef] [PubMed]
  58. Lindback, T.; Mols, M.; Basset, C.; Granum, P.E.; Kuipers, O.P.; Kovacs, A.T. CodY, a pleiotropic regulator, influences multicellular behaviour and efficient production of virulence factors in Bacillus cereus. Environ. Microbiol. 2012, 14, 2233–2246. [Google Scholar] [CrossRef] [Green Version]
  59. Hsueh, Y.H.; Somers, E.B.; Wong, A.C. Characterization of the codY gene and its influence on biofilm formation in Bacillus cereus. Arch. Microbiol. 2008, 189, 557–568. [Google Scholar] [CrossRef]
  60. Yan, F.; Yu, Y.; Gozzi, K.; Chen, Y.; Guo, J.H.; Chai, Y. Genome-wide investigation of biofilm formation in Bacillus cereus. Appl. Environ. Microbiol. 2017, 83, e00561-17. [Google Scholar] [CrossRef] [Green Version]
  61. Candela, T.; Fagerlund, A.; Buisson, C.; Gilois, N.; Kolsto, A.B.; Okstad, O.A.; Aymerich, S.; Nielsen-Leroux, C.; Lereclus, D.; Gohar, M. CalY is a major virulence factor and a biofilm matrix protein. Mol. Microbiol. 2019, 111, 1416–1429. [Google Scholar] [CrossRef] [Green Version]
  62. Kearns, D.B.; Chu, F.; Branda, S.S.; Kolter, R.; Losick, R. A master regulator for biofilm formation by Bacillus subtilis. Mol. Microbiol. 2005, 55, 739–749. [Google Scholar] [CrossRef] [PubMed]
  63. Majed, R.; Faille, C.; Kallassy, M.; Gohar, M. Bacillus cereus biofilms-same, only different. Front. Microbiol. 2016, 7, 1054. [Google Scholar] [CrossRef] [PubMed]
  64. Caro-Astorga, J.; Perez-Garcia, A.; de Vicente, A.; Romero, D. A genomic region involved in the formation of adhesin fibers in Bacillus cereus biofilms. Front. Microbiol. 2014, 5, 745. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Bianco, A.; Capozzi, L.; Monno, M.R.; Del Sambro, L.; Manzulli, V.; Pesole, G.; Loconsole, D.; Parisi, A. Characterization of Bacillus cereus group isolates from human bacteremia by whole-genome sequencing. Front. Microbiol. 2021, 11, 599524. [Google Scholar] [CrossRef] [PubMed]
  66. Karunakaran, E.; Biggs, C.A. Mechanisms of Bacillus cereus biofilm formation: An investigation of the physicochemical characteristics of cell surfaces and extracellular proteins. Appl. Microbiol. Biotechnol. 2011, 89, 1161–1175. [Google Scholar] [CrossRef]
  67. Whitfield, G.B.; Marmont, L.S.; Bundalovic-Torma, C.; Razvi, E.; Roach, E.J.; Khursigara, C.M.; Parkinson, J.; Howell, P.L. Discovery and characterization of a Gram-positive Pel polysaccharide biosynthetic gene cluster. PLoS Pathog. 2020, 16, e1008281. [Google Scholar] [CrossRef] [Green Version]
  68. Gao, T.; Ding, M.; Wang, Q. The recA gene is crucial to mediate colonization of Bacillus cereus 905 on wheat roots. Appl. Microbiol. Biotechnol. 2020, 104, 9251–9265. [Google Scholar] [CrossRef]
Microorganisms 10 00252 g001 550
Figure 1. Biofilm formation assay of Bacillus sp. B87. Submerged biofilms of Bacillus sp. B87 and B. cereus ATCC 14579 were visualized via crystal violet staining. Data represent the mean of three independent experiments, each consisting of six internal replicates. Error bars indicate standard deviations. Statistical analysis includes Student’s paired t-test (* p < 0.01).
Figure 1. Biofilm formation assay of Bacillus sp. B87. Submerged biofilms of Bacillus sp. B87 and B. cereus ATCC 14579 were visualized via crystal violet staining. Data represent the mean of three independent experiments, each consisting of six internal replicates. Error bars indicate standard deviations. Statistical analysis includes Student’s paired t-test (* p < 0.01).
Microorganisms 10 00252 g001
Microorganisms 10 00252 g002 550
Figure 2. Circular representation of the Bacillus sp. B87 draft genome. Circular genome visualization was generated using the circular viewer of PATRIC. Outer to center: contigs, forward CDS, reverse CDS, non-CDS features, AMR genes, VF Genes, transporters, and drug targets. The two inner tracks are GC content and GC skew.
Figure 2. Circular representation of the Bacillus sp. B87 draft genome. Circular genome visualization was generated using the circular viewer of PATRIC. Outer to center: contigs, forward CDS, reverse CDS, non-CDS features, AMR genes, VF Genes, transporters, and drug targets. The two inner tracks are GC content and GC skew.
Microorganisms 10 00252 g002
Microorganisms 10 00252 g003 550
Figure 3. Overview of the subsystem categories of the annotated draft whole genome of Bacillus sp. B87 from the RAST server. The pie chart shows the number of genes related to individual subsystems. The bar graph on the left reveals the subsystem coverage. The ratio of coding sequences annotated in the SEED subsystem (33%) and outside of the SEED subsystem (67%) is indicated.
Figure 3. Overview of the subsystem categories of the annotated draft whole genome of Bacillus sp. B87 from the RAST server. The pie chart shows the number of genes related to individual subsystems. The bar graph on the left reveals the subsystem coverage. The ratio of coding sequences annotated in the SEED subsystem (33%) and outside of the SEED subsystem (67%) is indicated.
Microorganisms 10 00252 g003
Microorganisms 10 00252 g004a 550Microorganisms 10 00252 g004b 550
Figure 4. Phylogenetic analysis of Bacillus sp. B87. (A) The codon tree method selects single-copy PATRIC PGFams and analyzes aligned proteins and coding DNA from single-copy genes using the program RAxML. (B,C) TYGS results of Bacillus sp. B87 based on 16S rRNA gene (B) and whole-genome (C) sequences.
Figure 4. Phylogenetic analysis of Bacillus sp. B87. (A) The codon tree method selects single-copy PATRIC PGFams and analyzes aligned proteins and coding DNA from single-copy genes using the program RAxML. (B,C) TYGS results of Bacillus sp. B87 based on 16S rRNA gene (B) and whole-genome (C) sequences.
Microorganisms 10 00252 g004aMicroorganisms 10 00252 g004b
Microorganisms 10 00252 g005 550
Figure 5. Genetic organization of the pelDEADAFG operon. The pelDEADAFG operons were compared using the Easyfig tool with homologous clusters found in related species of the genus Bacillus. Arrows indicate the transcription direction of each CDS with different colors (red, sodA; green, amiS; light blue, pel; grey, galE; white, unknown function (ND)). The shade of gray indicates the degree of nucleotide sequence homology (%) according to BLASTN.
Figure 5. Genetic organization of the pelDEADAFG operon. The pelDEADAFG operons were compared using the Easyfig tool with homologous clusters found in related species of the genus Bacillus. Arrows indicate the transcription direction of each CDS with different colors (red, sodA; green, amiS; light blue, pel; grey, galE; white, unknown function (ND)). The shade of gray indicates the degree of nucleotide sequence homology (%) according to BLASTN.
Microorganisms 10 00252 g005
Table
Table 1. Antimicrobial susceptibility testing of Bacillus sp. B87.
Table 1. Antimicrobial susceptibility testing of Bacillus sp. B87.
CategoryAntimicrobial AgentInterpretation
β-lactamAmpicillin (10 µg)R
Amoxicillin–clavulanic acid (20 µg/10 µg)R
Penicillin (10 U)R
AminoglycosidesGentamicin (10 µg)S
CarbapenemsImipenem (10 µg)S
GlycopeptidesVancomycin (30 µg)S
PhenicolsChloramphenicol (30 µg)S
FluoroquinolonesCiprofloxacin (5 µg)S
TetracyclinesTetracycline (30 µg)R
Folate pathway inhibitorsTrimethoprim–sulfamethoxazole
(1.25 µg/23.75 µg)		S
MacrolidesErythromycin (15 µg)S
S: Sensitive, R: Resistant.
Table
Table 2. Genomic features and annotation information of the chromosome of Bacillus sp. B87.
Table 2. Genomic features and annotation information of the chromosome of Bacillus sp. B87.
Genome FeaturesChromosome
Genome length (bp)5,448,163
Protein-coding genes5661
GC content (%)35.18
The number of tRNA77
The number of rRNA5
Contigs117
Contig L5011
Contig N50161,893
Table
Table 3. AMR genes prediction of Bacillus sp. B87 based on NDARO and CARD databases.
Table 3. AMR genes prediction of Bacillus sp. B87 based on NDARO and CARD databases.
GenesProductSource IDSource Organism
CARD database
Translation elongation factor TuYP_006374661.1Enterococcus faecium DO
BLA1Class A beta-lactamase (EC 3.5.2.6)AAR20595.1B. anthracis
FosBFosfomycin resistance protein FosBNP_831795.1B. cereus ATCC 14579
BcIISubclass B1 beta-lactamase
(EC 3.5.2.6) => BcII family		AAA22562.1B. cereus
dfrEThymidylate synthase (EC 2.1.1.45)AAD01867.1E. faecalis
NDARO database
Class A beta-lactamase (EC 3.5.2.6)WP_063842248.1B. cereus
Subclass B1 beta-lactamase
(EC 3.5.2.6) => BcII family		WP_000799223.1B. cereus group
Fosfomycin resistance protein FosBWP_000943763.1Bacillus
Tetracycline resistance, MFS efflux pump => Tet(45)WP_063855885.1Bhargavaea cecembensis
Table
Table 4. Prediction of biofilm formation genes of Bacillus sp. B87 based on RASTtk in PATRIC.
Table 4. Prediction of biofilm formation genes of Bacillus sp. B87 based on RASTtk in PATRIC.
Biofilm Formation GenesProduct
lsrRTranscriptional repressor of lsr operon
lsrKAutoinducer 2 (AI-2) kinase LsrK (EC 2.7.1.-)
lsrDAutoinducer 2 (AI-2) ABC transport system, membrane channel protein LsrD
lsrCAutoinducer 2 (AI-2) ABC transport system, membrane channel protein LsrC
Cupin domain protein in autoinducer 2 (AI-2)-related operon
N-acyl homoserine lactone hydrolase
3-hydroxy-5-phosphonooxypentane-2,4-dione thiolase (EC 2.3.1.245)
Autoinducer 2 (AI-2) ABC transporter, dimeric ATP-binding protein
Autoinducer 2 (AI-2) ABC transporter, substrate-binding protein
Table
Table 5. Identification of genes in Bacillus sp. B87 relevant to biofilm formation annotated using the KAAS database.
Table 5. Identification of genes in Bacillus sp. B87 relevant to biofilm formation annotated using the KAAS database.
KEGG OrthologyGenesProtein Product
Biofilm transcriptional regulators
K06284abrBAbrB family transcriptional regulator, transcriptional pleiotropic regulator of transition state genes
K03706codYTranscriptional pleiotropic repressor
K20480nprRHTH-type transcriptional regulator, quorum-sensing regulator NprR
K20391plcRHTH-type transcriptional regulator
K20390papRRegulatory peptide PapR
K06372sinIAntagonist of SinR
K19449sinRXRE family transcriptional regulator, master regulator for biofilm formation
K07699Spo0ATwo-component system, response regulator, stage 0 sporulation protein A
Matrix protein-encoding genes
K06336tasASpore coat-associated protein N
K13280sipWSignal peptidase I
Putative matrix polysaccharide synthesis genes
K07705lytRTwo-component system, LytTR family, response regulator LytT
K00012ugdUDPglucose 6-dehydrogenase
K21006pelAPolysaccharide biosynthesis protein PelA
K21009pelDPolysaccharide biosynthesis protein PelD
K21011pelFPolysaccharide biosynthesis protein PelF
K21012pelGPolysaccharide biosynthesis protein PelG
eDNA synthesis genes
K01939purAAdenylosuccinate synthase
K01923purCPhosphoribosylaminoimidazole-succinocarboxamide synthase
K23269purLPhosphoribosylformylglycinamidine synthase subunit PurL
Table
Table 6. Prediction of genes related to virulence factor of Bacillus sp. B87 according to Victors and VFDB databases, and BTyper tool.
Table 6. Prediction of genes related to virulence factor of Bacillus sp. B87 according to Victors and VFDB databases, and BTyper tool.
GenesProductSource IDSource Organism
Victors database
sodA2Superoxide dismutase [Mn] (EC 1.15.1.1)227818216B. anthracis str. CDC 684
sigBRNA polymerase sigma factor SigB227816152B. anthracis str. CDC 684
nosNitric oxide synthase oxygenase (EC 1.-.-.-)227818215B. anthracis str. CDC 684
codYGTP-sensing transcriptional pleiotropic repressor CodY227813264B. anthracis str. CDC 684
recARecA protein15926868S. aureus subsp. aureus N315
phnXPhosphonoacetaldehyde hydrolase (EC 3.11.1.1)47526609B. anthracis str. ‘Ames Ancestor’
sodCSuperoxide dismutase [Cu-Zn] precursor (EC 1.15.1.1)227817676B. anthracis str. CDC 684
sodA1Superoxide dismutase [Mn] (EC 1.15.1.1)227817051B. anthracis str. CDC 684
clpXATP-dependent Clp protease ATP-binding subunit ClpX227817253B. anthracis str. CDC 684
VFDB database
nheCEnterotoxin CVFG016286B. cereus ATCC 10987
inhAImmune inhibitor A, metalloprotease (EC 3.4.24.-)VFG016338B. anthracis str. Sterne
nheBNon-hemolytic enterotoxin lytic component L1VFG016278B. cereus ATCC 10987
nheANon-hemolytic enterotoxin AVFG016270B. cereus ATCC 10987
BTyper tool
bpsFBacillus cereus exo-polysaccharide operon gene F tyrosine protein kinase [plasmid pBC218] B. cereus str. G9241
entFMEnterotoxin B. cereus ATCC 14579
bceTDiarrheal toxin B. cereus
plcA1-Phosphatidylinositol phosphodiesterase precursor B. cereus ATCC 14579
entAEnterotoxin/cell-wall binding protein B. cereus ATCC 14579
bpsEBacillus cereus exo-polysaccharide operon gene E UTP--glucose-1-phosphate uridylyltransferase [plasmid pBC218] B. cereus str. G9241
inhA2Immune inhibitor A precursor B. cereus ATCC 14579
nheCEnterotoxin C B. cereus ATCC 14579
cerACereolysin A B. cereus
bpsHBacillus cereus exo-polysaccharide operon gene H LytR family transcriptional regulator [plasmid pBC218] B. cereus str. G9241
inhA1Immune inhibitor A precursor B. cereus ATCC 14579
nheANon-hemolytic enterotoxin lytic component L2 B. cereus ATCC 14579
nheBNon-hemolytic enterotoxin lytic component L1 B. cereus ATCC 14579
cerBCereolysin B
plcBPhospholipase C B. cereus ATCC 14579
sphSphingomyelinase C B. anthracis str. Ames
plcRTranscriptional regulator B. anthracis str. Ames
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sornchuer, P.; Saninjuk, K.; Prathaphan, P.; Tiengtip, R.; Wattanaphansak, S. Antimicrobial Susceptibility Profile and Whole-Genome Analysis of a Strong Biofilm-Forming Bacillus Sp. B87 Strain Isolated from Food. Microorganisms 2022, 10, 252. https://doi.org/10.3390/microorganisms10020252

AMA Style

Sornchuer P, Saninjuk K, Prathaphan P, Tiengtip R, Wattanaphansak S. Antimicrobial Susceptibility Profile and Whole-Genome Analysis of a Strong Biofilm-Forming Bacillus Sp. B87 Strain Isolated from Food. Microorganisms. 2022; 10(2):252. https://doi.org/10.3390/microorganisms10020252

Chicago/Turabian Style

Sornchuer, Phornphan, Kritsakorn Saninjuk, Parisa Prathaphan, Rattana Tiengtip, and Suphot Wattanaphansak. 2022. "Antimicrobial Susceptibility Profile and Whole-Genome Analysis of a Strong Biofilm-Forming Bacillus Sp. B87 Strain Isolated from Food" Microorganisms 10, no. 2: 252. https://doi.org/10.3390/microorganisms10020252

APA Style

Sornchuer, P., Saninjuk, K., Prathaphan, P., Tiengtip, R., & Wattanaphansak, S. (2022). Antimicrobial Susceptibility Profile and Whole-Genome Analysis of a Strong Biofilm-Forming Bacillus Sp. B87 Strain Isolated from Food. Microorganisms, 10(2), 252. https://doi.org/10.3390/microorganisms10020252

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