- Research article
- Open Access
Assessment and molecular characterization of Bacillus cereus isolated from edible fungi in China
BMC Microbiology volume 20, Article number: 310 (2020)
Bacillus cereus is a foodborne pathogen commonly found in nature and food and can cause food spoilage and health issues. Although the prevalence of B. cereus in foods has been reported worldwide, the extent of contamination in edible fungi, which has become increasingly popular as traditional or functional food, is largely unknown. Here we investigated the prevalence, toxin genes’ distribution, antibiotic resistance, and genetic diversity of B. cereus isolated from edible fungi in China.
Six hundred and ninety-nine edible fungi samples were collected across China, with 198 (28.3%) samples found to be contaminated by B. cereus, with an average contamination level of 55.4 most probable number (MPN)/g. Two hundred and forty-seven B. cereus strains were isolated from the contaminated samples. Seven enterotoxin genes and one cereulide synthetase gene were detected. The detection frequencies of all enterotoxin genes were ≥ 80%, whereas the positive rate of the cesB gene in B. cereus was 3%. Most isolates were resistant to penicillins, β-lactam/β-lactamase inhibitor combinations, cephems, and ansamycins, but were susceptible to penems, aminoglycosides, macrolides, ketolide, glycopeptides, quinolones, phenylpropanol, tetracyclines, lincosamides, streptogramins, and nitrofurans. Meanwhile, 99.6% of all isolates displayed multiple antimicrobial resistance to three or more classes of antimicrobials. Using genetic diversity analysis, all isolates were defined in 171 sequence types (STs), of which 83 isolates were assigned to 78 new STs.
This study provides large-scale insight into the prevalence and potential risk of B. cereus in edible fungi in China. Approximately one-third of the samples were contaminated with B. cereus, and almost all isolates showed multiple antimicrobial resistance. Detection frequencies of all seven enterotoxin genes were equal to or more than 80%. These new findings may indicate a need for proper pre-/post-processing of edible fungi to eliminate B. cereus, thereby preventing the potential risk to public health.
Foodborne disease is an important public health issue worldwide, which is among the leading causes of morbidity and mortality [1,2,3,4]. Data from the foodborne disease outbreak surveillance and reporting system showed 5021 foodborne disease outbreaks were reported between 2001 to 2010 in China . These outbreaks involved 140,101 cases and 1427 deaths, with fatality rate of 1.0%, and pathogenic microorganisms accounted for most of the incidents and cases . Ninety-eight foodborne disease outbreaks were reported in Shandong, China in 2014, involving 1238 patients and four deaths, mostly caused by pathogenic microorganisms . The top five bacterial pathogens responsible for foodborne outbreaks from 2011 to 2016 in China were Vibrio parahaemolyticus, Salmonella spp., Staphylococcus aureus (enterotoxin), Bacillus cereus, and diarrheagenic Escherichia coli . Notably, there had been 4342 foodborne cases caused by B. cereus during these 6 years.
B. cereus is a Gram-positive, spore-forming bacterium widespread in natural environment [8, 9]. It can cause two types of food poisoning in humans, including diarrheal and emetic syndromes, along with a variety of local and systemic infections, such as endophthalmitis, endocarditis, meningitis, osteomyelitis, wound infections, and septicemia [10, 11]. It causes diarrhea and emesis via the production of diarrheal enterotoxins and emetic toxin, respectively. The enterotoxins include hemolysin BL (Hbl), nonhemolytic enterotoxin (Nhe), and the single-protein enterotoxin, cytotoxin K (CytK) [9, 12, 13]. Hbl and Nhe are tripartite toxins, in which all three components are necessary for maximal cytotoxic activity [14, 15]. The emetic toxin cereulide is a cyclic depsipeptide toxin encoded by the ces gene cluster [16, 17]. Emetic symptoms usually occur within 0.5-6 h after consuming contaminated foods . Cereulide is pre-formed in food and is difficult to inactivate either during food processing or in the gastrointestinal tract, owing to its high resistance to heat treatments, extreme pH conditions, and enzymolysis [9, 16].
To date, antibiotic treatment is one choice in clinical treatment for severe food poisoning and other body or tissue infections caused by B. cereus [18,19,20]. However, studies have reported B. cereus to be resistant to different antibiotic agents, thus constituting a fundamental problem worldwide [21,22,23]. Information is therefore necessary for an antibiotic resistance profile of this bacterium.
B. cereus has been isolated from many foodstuffs, such as rice, pasta, vegetables, meat and meat products, soup, and milk and other dairy products [24,25,26,27,28,29]. However, the actual incidence rate of B. cereus in edible fungi remains largely unknown. From culinary perspective, edible fungi are consumed worldwide. The nutritional properties and health benefits of ‘medicinal’ edible fungi have been known in China for over 2000 years. Besides, edible fungi have been a part of traditional Chinese medicine, and have also been used in Japan and Malaysia . B. cereus contamination could lead to spoilage of edible fungi and cause food poisoning. It is therefore necessary to investigate the prevalence of B. cereus in edible fungi.
The objective of this study was to investigate the contamination level and toxin gene distribution, including seven diarrheal enterotoxin genes and one emetic toxin gene across the 247 B. cereus strains isolated from edible fungi in China. In addition, we characterized the isolates based on their antimicrobial resistance and genotypic diversity. Thereby, we provided a systematic risk assessment for B. cereus isolated from edible fungi in China.
Contamination of edible fungi with B. cereus
The prevalence of B. cereus in 699 edible fungi samples examined in this study is shown in Table 1. B. cereus was detected in 28.3% (198/699) of all samples collected and 247 isolates were isolated from these contaminated samples. Average contamination level for the positive samples was 55.4 MPN/g. 20.6% (45/218) of Flammulina velutipes, 34.5% (39/113) of Pleurotus ostreatus, 57.8% (63/109) of Lentinus edodes, 14.7% (14/95) of Pleurotus eryngii, 16.0% (15/94) of Hypsizygus marmoreus, and 31.4% (22/70) of other species were contaminated with B. cereus, respectively. The contamination levels of 81.8% (162/198) of positive samples ranged from 3 to 1100 MPN/g, with that in five samples (5/198, 2.5%) exceeding 1100 MPN/g.
Distribution of toxin genes
The positive rate of all seven enterotoxin genes was ≥80%. Specifically, 80, 89, 91, 100, 90, 98, and 83% of any of the isolates harbored hblA, hblC, hblD, nheA, nheB, nheC, and cytK, respectively (Fig. 1; Additional file 1: Table S1). In contrast, 3% of the isolates were positive for the cereulide synthetase gene B (cesB). All isolates were positive for at least two toxin genes, and 157 (64%) isolates harbored seven or more toxin genes. Among the 30 different gene profiles, the most commonly observed gene profile was hblA-hblC-hblD-nheA-nheB-nheC-cytK, possessed by 62% of the isolates.
Antibiotic susceptibility of B. cereus in edible fungi
The antibiotic susceptibility results of 247 B. cereus isolates to 20 selected antimicrobials are shown in Fig. 2. The diameters of inhibition zones in the antimicrobial susceptibility test are included in Additional file 2: Table S2. All isolates were resistant to ampicillin (AMP) and penicillin (P). Most isolates were resistant to amoxicillin-clavulanic acid (AMC; 99.2%), cephalothin (KF; 85.4%), cefoxitin (FOX; 91.9%), and rifampin (RD; 91.1%). Over 93% of the isolates showed susceptibility to imipenem (IPM; 99.2%), gentamicin (CN; 98.0%), ciprofloxacin (CIP; 93.1%), and chloramphenicol (C; 93.5%). Moreover, most isolates showed intermediate resistance to clindamycin (DA; 83.8%) and quinupristin-dalfopristin (QD; 72.1%) (Additional file 3: Table S3).
All isolates were resistant to more than two antibiotics. Five isolates (2.0%) were resistant to ≥10 antibiotics, and 33 isolates (13.4%) showed resistance to at least eight antibiotics (Additional file 4: Table S4). For the most common antimicrobial resistant pattern, 33.6% of the isolates harbored the profile AMP-P-AMC-KF-FOX-RD (Additional file 4: Table S4). In addition, 99.6% of all isolates displayed multiple antimicrobial resistance to three and more classes of antimicrobials (Fig. 2). Approximately 91.9 and 39.3% of the isolates showed multiple antimicrobial resistance to more than three or four classes of antimicrobials, respectively (Fig. 2).
Multi-locus sequence typing and clustering of the isolates
Of the 171 STs assigned to all isolates, ST770 was the most prevalent but only included 11 isolates, and 135 STs (78.9%) included only one isolate. Eighty-three (33.6%) of the 247 isolates belonged to one of 78 novel STs. Of all the isolates, 141 singletons and seven clonal complexes (CCs) were assigned. The most prevalent CC was the ST-142 complex, which included 40 isolates (Additional file 4: Table S4). The ST-8, ST-18, ST-23, ST-97, ST-111, and ST-205 complexes included 3, 27, 10, 5, 2, and 19 isolates, respectively (Fig. 3; Additional file 4: Table S4). Regarding the phylogenetic relationship, all isolates were grouped into ten clusters with the cut-off value of 43% similarity (Additional file 5: Figure S1) and the cluster numbers were defined depended on the description of Guinebretière et al. . Cluster III contained the most strains. Cluster IX only contained one isolate identified in this study, and cluster X only possessed Bacillus anthracis ATCC4728, indicating these strains were far away from the others. Isolates in cluster I were phylogenetically map together with Bacillus pseudomycoides DSM 12442, and isolates in cluster III were map together with Bacillus weihenstephanensis WSBC 10204 and Bacillus cytotoxicus NVH 391-98, which indicating that these strains may be phylogenetically close with each other.
B. cereus is an invasive and opportunistic pathogen that can cause gastrointestinal diseases and severe nosocomial infections. To date, prevalence studies of pathogenic B. cereus isolated from edible fungi are limited. We therefore investigated the contamination characteristics of B. cereus from edible fungi collected between 2011 to 2016. B. cereus was present in 28.3% of edible fungi samples collected from major cities of China. This suggested that more attention should be paid to the extent of contamination in edible fungi. Contamination may occur at any point during cultivation, processing, storage, and sale of the edible fungi. Compared to previous surveys with other foodstuffs (8.2% in infant formula, 3.8% in rice flour, and 9% in ready-to-eat foods), the contamination rate of B. cereus isolated from edible fungi was higher [32, 33]. Unlike our previous studies which showed another important foodborne pathogen, Listeria monocytogenes, to be highly associated with F. velutipes [34, 35]; here, B. cereus was frequently detected in L. edodes. Since > 50% of the L. edodes samples were positive for B. cereus, further studies regarding the contamination source would be recommended, considering the B. cereus strains to possibly be the contamination sources for the environment and other foodstuffs. In addition, the close correlation between B. cereus and L. edodes should also be investigated. Although there is no specified standard for the microbiological limits of B. cereus in edible fungi in China, the amounts between 103 and 105 CFU/g of B. cereus in ready-to-eat foods are “acceptable”, but not “satisfactory”, according to the Microbiological Guidelines for Food of Hong Kong, China , hence suggesting a potential for public health risk. Studies have shown that more than 100 CFU/g of B. cereus in food samples may lead to foodborne infection [37, 38]. In this study, contamination levels in 2.5% of positive samples exceeded 1100 MPN/g (Table 1). A B. cereus strain (strain number 49-1), isolated from one of these samples, harbored the cesB gene and was simultaneously identified as a ST26 strain. The reported emetic type strains F4810/72 and NC7401 [39, 40], were also identified as ST26 based on their complete genome sequences in the NCBI genome database. These results indicated a potential hazard of these contaminated edible fungi and the need for adequate cleaning and heating measures before consumption.
According to toxin gene distribution, 30 different toxin gene profiles were identified from the isolates. The detection frequency of nheABC (89%) was similar to the findings in previous studies reporting that over 90% of isolates, from a variety of food samples, contained the nheABC gene cluster [28, 29, 41]. Besides, the detection frequencies of nheC, nheB, and nheA were close to that of B. cereus isolated from different food samples (100, 99, and 96%, respectively) , in line with these three genes being generally prevalent in B. cereus isolates. The detection frequencies of nheA, nheB, and nheC were not equal, which is also consistent with the previous studies [38, 42, 43]. These three genes were not always simultaneously found in the same isolate for unknown reasons, as were the genes hblA, hblC, and hblD. One possible explanation could be that the general primers used for detection cannot cover the polymorphism of target genes in different strains and affected the amplification efficiency in PCR. One hundred and ninety-one (77%) isolates simultaneously carried the three hemolysin genes, hblA, hblC, and hblD, which had higher detection frequencies than in other foodstuffs (55% in raw milk samples, and 45% in pasteurized milk samples) [27, 28]. The most common toxin gene profile among all isolates (62%) was hblA-hblC-hblD-nheA-nheB-nheC-cytK, revealing that B. cereus isolated from edible fungi mainly harbored the diarrheal enterotoxin-encoding genes, in line with previous research . Although the detection frequency of cesB is lower than that of all the enterotoxin genes tested, it is still higher than that in a previous study , indicating the potential virulence of these isolates.
B. cereus has been reported to be resistant to different antimicrobials [45,46,47]. In this study, we found that B. cereus isolates from edible fungi were mainly resistant to β-lactams (except cefotetan and imipenem) and ansamycins, especially ampicillin and penicillin. Previous studies had shown that B. cereus harbors the antibiotic resistance genes encoding β-lactamases [11, 23, 48, 49]. Therefore, B. cereus has the ability to degrade β-lactams and show resistant characteristics [23, 28, 29, 48, 49]. In contrast, B. cereus isolated from edible fungi were still sensitive to aminoglycosides, ketolide, glycopeptides, quinolones, and phenylpropanol antibiotics, in agreement with some previous studies [24, 27,28,29, 50]. Ciprofloxacin, clindamycin, and aminoglycoside antibiotics have been highly effective in the treatment of diseases caused by B. cereus [11, 51,52,53]. However, 83.8% of the isolates in this study showed intermediate resistance to clindamycin, indicating that clindamycin may be less effective at inhibiting these isolates. Lee et al. have reported that B. cereus isolates from Sunsik were resistant to ampicillin, cefoxitin, and penicillin, which were in concordance with our results. However, their isolates were fully sensitive to chloramphenicol, gentamicin, imipenem, sulfamethoxazole, and tetracycline , whereas 22.7% of our isolates showed resistance to sulfamethoxazole. Luna et al.  and Ikram et al.  had reported B. cereus to be resistant to ampicillin, penicillin, and sensitive to chloramphenicol, ciprofloxacin, and gentamicin, similar to our findings; however, 72.1% of our isolates showed intermediate resistance to quinupristin-dalfopristin and 91.1% of our isolates showed resistance to rifampin. These differences may be caused by the differences in sample sources, and sampling sites. Remarkably, we found 99.6% of all the isolates to be simultaneously resistant to antibiotics from three or more antimicrobial classes, with 91.9 and 39.3% of the isolates showing multiple antimicrobial resistance to more than three or four classes, respectively. The emergence of multiple antimicrobial resistant bacteria has led to a reduction in the types and quantities of effective antibacterial drugs available for pathogen elimination, and may further lead to a decline in the effectiveness of clinical treatment for bacterial infection. The potential hazards of multiple antimicrobial resistant isolates may need further investigation.
Of the isolates studied, 33.6% (83/247) were assigned to new STs, indicating the diverse bacterial resources not yet identified in this study. Of the 171 STs, 135 had only one isolate; thus, at least 54.7% of the isolates had a unique ST from each other, revealing the genetic diversity of isolates from edible fungi. The most prevalent ST was ST770, and the dominant CC was the ST-142 clonal complex, contradicting with previous studies [57, 58]. The reason for this phenomenon may be the difference in isolation sources of the strains. From clustering results in evolutionary analysis (Additional file 5: Figure S1), we could find specific characteristics of antibiotic resistance of the isolates. Clusters IX and X possessed only one strain each. Only isolates in cluster VI showed no QD (quinupristin-dalfopristin) resistance. Isolates in cluster I, III, IV and V showed resistance to TEL (telithromycin), and isolates in cluster III, V and VI showed resistance to TEC (teicoplanin). There were only 20, 7, and 8 isolates showing resistance to QD, TEL, and TEC. Three (isolates 49-1, 633-1C, and 1284) of the seven cesB-positive isolates identified in this study were grouped in cluster III, and they were likely to be gathered with the B. cereus emetic type strains NC7401 and F4810/72. Thus, the pathogenic potential of these strains may be higher than the others, and need to be investigated in future. Interestingly, the other four cesB-positive isolates were randomly distributed in different branches, suggesting the phylogenetic relationship between them to be distant, probably evolved from different origins.
The prevalence of B. cereus in about one third of the edible fungi samples in China indicated a source of potential infectious exposure to the consumers. Edible fungi samples with a B. cereus contamination level of ≥1100 MPN/g suggest the importance of performing extensive cleaning and/or processing before consumption. Of the 247 isolates, 171 STs were assigned, of which 135 STs only contained one isolate and 78 new STs were defined. Moreover, 64% of the isolates harbored all the detected (seven) enterotoxins, among which four isolates further harbored the cereulide synthetase gene cesB. Three cesB-positive strains were identified as the same ST as B. cereus clinical emetic type strains and simultaneously phylogenetically clustered together with them. Isolates exhibiting multiple antimicrobial resistance and multiple toxin genes may indicate potential health risks related to B. cereus from edible fungi. Therefore, further studies would be recommended to identify the pathogenic potential of these strains.
Edible fungi sample collection
From July 2011 to January 2016, 699 edible fungi samples (including 218 Flammulina velutipes, 113 Pleurotus ostreatus, 109 Lentinus edodes, 95 Pleurotus eryngii, 94 Hypsizygus marmoreus, and 70 other species) were obtained from 39 cities, covering all the provincial capital cities of China. The details of sampling positions are shown in Additional file 6: Figure S2.
The samples were stored in sealed bags and placed in cooling boxes (1-4 °C). They were transported to the laboratory within 12 h and analyzed immediately according to the general sample collection guidelines of the National Food Safety Standard .
Identification of B. cereus in edible fungi
All samples were subjected to qualitative and quantitative analysis for B. cereus, according to the B. cereus test standards given by the National Food Safety Standard [60, 61] with some modifications, as previously described [28, 29]. Twenty-five grams of the collected samples were cut and mixed with 225 mL phosphate buffered saline (PBS, 0.01 mol/L) under aseptic conditions. Then, 10- and 100-fold dilutions of the homogenized solution were prepared. One milliliter of each dilution was inoculated into nine milliliters of tryptic soy broth (TSB) with polymyxin B (Huankai, Guangzhou, China) in different tubes. Each experiment was repeated in triplicate in parallel. The tubes were incubated at 30 °C for 8-18 h. The samples were streaked onto mannitol egg yolk polymyxin (MYP) agar plates (Huankai, Guangzhou, China) and incubated at 30 °C for 24-48 h.
The presumptive colonies, in pink color with a pink halo, were selected and streaked on the chromogenic B. cereus agar plates (Huankai, Guangzhou, China) and then incubated at 30 °C for 24 h. The blue-green colonies were selected and streaked on the nutrient agar plates (Huankai, Guangzhou, China). The selected colonies were incubated and used for further analysis, including biochemical identification, parasporal crystal observation, hemolysis test, MYP agar plate test, root growth observation, catalase test, motility test, nitrate reduction test, casein decomposition test, lysozyme tolerance test, glucose utilization test, and acetyl methyl alcohol test, described in the National Food Safety Standard [60, 61] and the Bacteriological Nnalytical Manuals of the U.S. Food and Drug Administration . The isolate, which 1) produces large Gram-positive rods with spores that does not swell the sporangium; 2) produces lecithinase and does not ferment mannitol on MYP agar; 3) grows and produces acid from glucose anaerobically; 4) reduces nitrate to nitrite (a few strains may be negative); 5) produces acetylmethylcarbinol (VP-positive); 6) decomposes L-tyrosine; 7) grows in the presence of 0.001% lysozyme; 8) is actively motile and strongly hemolytic; 9) does not produce rhizoid colonies or protein toxin crystals, can be identified as B. cereus. The most probable number (MPN), which represents the number of corresponding bacteria most likely to be present per gram of sample, of B. cereus for the samples was determined according to the National Food Safety Standard [60, 61], the Bacteriological Nnalytical Manuals of the U.S. Food and Drug Administration , and previous studies [29, 63, 64].
Detection of toxin genes
Genomic DNA was extracted from B. cereus using a HiPure DNA Extraction Kit (Magen, Guangzhou, China), according to the manufacturer’s instructions. Eight toxin genes, including seven enterotoxin genes (nheA, nheB, nheC, hblC, hblD, hblA, and cytK) and a cereulide synthetase gene B (cesB) were detected. The primers used for the detection of toxin genes and their annealing temperature are shown in Additional file 7: Table S5.
The polymerase chain reaction (PCR) mixture (25 μL) comprised of 50 ng genomic DNA, 0.5 μL of each primer solution (concentration 10 μM), and 12.5 μL PCR Premix Taq™ (Takara, China). The amplification process included: initial denaturation at 94 °C for 5 min, 30 cycles of 94 °C for 1 min, 55 °C or 58 °C (58 °C for cesB and 55 °C for the other toxin genes) for 1 min, 72 °C for 1 min, and a final 10 min extension at 72 °C [17, 43, 65, 66]. PCR was conducted in a TONE-96G PCR Thermal Cycler (Analytik Jena, Jena, Germany). Amplicons were subjected to electrophoresis using 1.5% agarose gel containing 0.01% Gold View. The gels were visualized by a UV Imaging System (Bio-Rad, Hercules, USA). A 2000-bp DNA ladder (Dongsheng, Guangzhou, China) was used as a molecular weight marker. The images were captured in TIFF file format for further analysis.
Antimicrobial susceptibility test
Antimicrobial susceptibility was evaluated using the Kirby-Bauer disk diffusion method on Mueller-Hinton (MH) agar for 20 selected antimicrobial drugs belonging to 16 different classes [28, 29, 67]. The details of antimicrobial agents (Oxoid, Basingstoke, UK) are shown in Fig. 2, and Additional file 3: Table S3. Following the methods of CLSI  and previous research [23, 68,69,70,71], results were expressed as resistant (R), intermediate (I), and sensitive (S). The zone diameter interpretive criteria are shown in Additional file 8: Table S6, with reference to CLSI . In addition, multiple antimicrobial resistant isolates, showing resistance to antibiotics from three or more antimicrobial classes, were also evaluated.
Multi-locus sequence typing and clustering of different isolates
Seven housekeeping genes (glp, gmk, ilvD, pta, pur, pycA, and tpi) were amplified, sequenced, and analyzed to obtain the multi-locus sequence type of the isolates according to the MLST scheme which was developed by Martin Maiden, Gus Priest, and Maggie Barker  on the PubMLST official website (http://pubmlst.org/bcereus/info/primers.shtml). PHYLOViZ 2.0 software (Instituto de Microbiologia, Portugal) was used to construct a minimum spanning tree of the isolates , to visualize the genetic diversity and inter-strain relationships. The phylogenetic relationship between different isolates was studied depending on the different STs and seven housekeeping genes, using the BioNumerics software (version 7.6; Applied Maths, Belgium) with the arithmetic mean (UPGMA) method. The strains used in the phylogenetic study contained 247 isolates identified in this study, and nine isolates (B. cereus ATCC14579, Bacillus mycoides DSM 2048, Bacillus pseudomycoides DSM 12442, Bacillus weihenstephanensis WSBC 10204, Bacillus anthracis ATCC4728, Bacillus thuringiensis ATCC10792, Bacillus cytotoxicus NVH 391-98 and two clinical emetic type strains B. cereus NC7401 and B. cereus F4810/72) from the B. cereus MLST database (www.mlstoslo.uio.no/index.html) or NCBI genome database (https://www.ncbi.nlm.nih.gov/genome/). The names of isolates, MLST sequence type (ST), MLST clonal complex (CC), antibiotic resistance patterns, and presence of toxin genes were also included.
Availability of data and materials
Genetic polymorphisms data used for Additional file 5: Figure S1 has been list in Additional file 4: Table S4. The sequence type and corresponding DNA sequences of seven housekeeping genes (glp, gmk, ilvD, pta, pur, pycA, and tpi) were all available and can be downloaded from the public MLST database (https://pubmlst.org/bcereus/) according to the STs of the isolates list in Additional file 4: Table S4. All other data generated or analyzed during this study are included in this manuscript and its supplementary information files.
Most probable number
- cesB :
Cereulide synthetase gene B
Multi-locus sequence typing
Phosphate buffered saline
Tryptic soy broth
Mannitol egg yolk polymyxin
Polymerase chain reaction
Multiple antimicrobial resistance
Wu G, Yuan Q, Wang L, Zhao J, Chu Z, Zhuang M, Zhang Y, Wang K, Xiao P, Liu Y, et al. Epidemiology of foodborne disease outbreaks from 2011 to 2016 in Shandong Province, China. Medicine (Baltimore). 2018;97(45):e13142.
Jahan S. Epidemiology of Foodborne Illness. In., vol. 336; 2012.
Helms M, Vastrup P, Gerner-Smidt P, Mølbak K. Short and long term mortality associated with foodborne bacterial gastrointestinal infections: registry based study. BMJ. 2003;326(7385):357.
Alcorn T, Ouyang Y. China's invisible burden of foodborne illness. Lancet. 2012;379(9818):789–90.
Xu J, Zhang J. Analysis of foodborne disease outbreaks in China from 2001 to 2010. Chin Agri Sci Bull. 2012;28:313–6 (in Chinese).
Zong JW, Zhou JY, Zhao JS, Chu ZH, Wang LS, Yan JY. Monitoring data of foodborne disease outbreaks in Shandong Province, 2014. Pract Prev Med. 2018;1(25):122–5 (in Chinese).
Liu J, Bai L, Li W, Han H, Fu P, Ma X, Bi Z, Yang X, Zhang X, Zhen S, et al. Trends of foodborne diseases in China: lessons from laboratory-based surveillance since 2011. Front Med. 2018;12(1):48–57.
European Food Safety Authority (EFSA). Opinion of the scientific panel on biological hazards on Bacillus cereus and other Bacillus spp. in foodstuffs. EFSA J. 2005;175:1–48. https://doi.org/10.2903/j.efsa.2005.175.
Stenfors Arnesen LP, Fagerlund A, Granum PE. From soil to gut: Bacillus cereus and its food poisoning toxins. FEMS Microbiol Rev. 2008;32(4):579–606.
Drobniewski FA. Bacillus cereus and related spcies. Clin Microbiol Rev. 1993;6(4):324–38.
Bottone EJ. Bacillus cereus, a volatile human pathogen. Clin Microbiol Rev. 2010;23(2):382–98.
Lund T, De Buyser ML, Granum PE. A new cytotoxin from Bacillus cereus that may cause necrotic enteritis. Mol Microbiol. 2000;38(2):254–61.
Ehling-Schulz M, Guinebretière MH, Monthán A, Berge O, Fricker M, Svensson B. Toxin gene profiling of enterotoxic and emetic Bacillus cereus. FEMS Microbiol Lett. 2006;260(2):232–40.
Beecher DJ, Macmillan JD. Characterization of the components of hemolysin BL from Bacillus cereus. Infect Immun. 1991;59(5):1778–84.
Lindbäck T, Fagerlund A, Rødland MS, Granum PE. Characterization of the Bacillus cereus Nhe enterotoxin. Microbiology. 2004;150(Pt 12):3959–67.
Agata N, Ohta M, Mori M, Isobe M. A novel dodecadepsipeptide, cereulide, is an emetic toxin of Bacillus cereus. FEMS Microbiol Lett. 1995;129(1):17–20.
Ehling-Schulz M, Vukov N, Schulz A, Shaheen R, Andersson M, Märtlbauer E, Scherer S. Identification and partial characterization of the nonribosomal peptide synthetase gene responsible for cereulide production in emetic Bacillus cereus. Appl Environ Microbiol. 2005;71(1):105–13.
Kervick GN, Flynn HW Jr, Alfonso E, Miller D. Antibiotic therapy for Bacillus species infections. Am J Ophthalmol. 1990;110(6):683–7.
Dierick K, Van Coillie E, Swiecicka I, Meyfroidt G, Devlieger H, Meulemans A, Hoedemaekers G, Fourie L, Heyndrickx M, Mahillon J. Fatal family outbreak of Bacillus cereus-associated food poisoning. J Clin Microbiol. 2005;43(8):4277–9.
Veysseyre F, Fourcade C, Lavigne JP, Sotto A. Bacillus cereus infection: 57 case patients and a literature review. Med Mal Infect. 2015;45(11–12):436–40.
Agersø Y, Jensen LB, Givskov M, Roberts MC. The identification of a tetracycline resistance gene tet(M), on a Tn916-like transposon, in the Bacillus cereus group. FEMS Microbiol Lett. 2002;214(2):251–6.
Brown JR, Gentry D, Becker JA, Ingraham K, Holmes DJ, Stanhope MJ. Horizontal transfer of drug-resistant aminoacyl-transfer-RNA synthetases of anthrax and gram-positive pathogens. EMBO Rep. 2003;4(7):692–8.
Torkar KG, Bedenić B. Antimicrobial susceptibility and characterization of metallo-beta-lactamases, extended-spectrum beta-lactamases, and carbapenemases of Bacillus cereus isolates. Microb Pathog. 2018;118:140–5.
Wijnands LM, Dufrenne JB, Rombouts FM, In't Veld PH, Van Leusden FM. Prevalence of potentially pathogenic Bacillus cereus in food commodities in The Netherlands. J Food Prot. 2006;69(11):2587–94.
Chaves JQ, Pires ES, Vivoni AM. Genetic diversity, antimicrobial resistance and toxigenic profiles of Bacillus cereus isolated from food in Brazil over three decades. Int J Food Microbiol. 2011;147(1):12–6.
Chon JW, Yim JH, Kim HS, Kim DH, Kim H, Oh DH, Kim SK, Seo KH. Quantitative prevalence and toxin gene profile of Bacillus cereus from ready-to-eat vegetables in South Korea. Foodborne Pathog Dis. 2015;12(9):795–9.
Cui Y, Liu X, Dietrich R, Märtlbauer E, Cao J, Ding S, Zhu K. Characterization of Bacillus cereus isolates from local dairy farms in China. FEMS Microbiol. Lett. 2016;363(12):fnw096.
Gao T, Ding Y, Wu Q, Wang J, Zhang J, Yu S, Yu P, Liu C, Kong L, Feng Z. Prevalence, virulence genes, antimicrobial susceptibility, and genetic diversity of Bacillus cereus isolated from pasteurized milk in China. Front Microbiol. 2018;9:533.
Yu P, Yu S, Wang J, Guo H, Zhang Y, Liao X, Zhang J, Wu S, Gu Q, Xue L, et al. Bacillus cereus isolated from vegetables in China: incidence, genetic diversity, virulence genes, and antimicrobial resistance. Front Microbiol. 2019;10:948.
Sadler M. Nutritional properties of edible fungi. Nutr Bull. 2003;28(3):305–8.
Guinebretière MH, Thompson FL, Sorokin A, Normand P, Dawyndt P, Ehling-Schulz M, Svensson B, Sanchis V, Nguyen-The C, Heyndrickx M, et al. Ecological diversification in the Bacillus cereus group. Environ Microbiol. 2008;10(4):851–65.
Zhang Y, Chen J, Feng C, Zhan L, Zhang J, Li Y, Yang Y, Chen H, Zhang Z, Zhang Y, et al. Quantitative prevalence, phenotypic and genotypic characteristics of Bacillus cereus isolated from retail infant foods in China. Foodborne Pathog Dis. 2017;14(10):564–72.
Forero AY, Galindo M, Morales GE. Isolation of Bacillus cereus in school restaurants in Colombia. Biomedica. 2018;38(3):338–44.
Chen M, Wu Q, Zhang J, Guo W, Wu S, Yang X. Prevalence and contamination patterns of Listeria monocytogenes in Flammulina velutipes plants. Foodborne Pathog Dis. 2014;11(8):620–7.
Wu S, Wu Q, Zhang J, Chen M, Yan ZA, Hu H. Listeria monocytogenes prevalence and characteristics in retail raw foods in China. PLoS One. 2015;10(8):e0136682.
Food and environmental hygiene department. Microbiological Guidelines for Food. Hong Kong: Food and environmental hygiene department. 2014.
Ghelardi E, Celandroni F, Salvetti S, Barsotti C, Baggiani A, Senesi S. Identification and characterization of toxigenic Bacillus cereus isolates responsible for two food-poisoning outbreaks. FEMS Microbiol Lett. 2002;208(1):129–34.
Glasset B, Herbin S, Guillier L, Cadel-Six S, Vignaud ML, Grout J, Pairaud S, Michel V, Hennekinne JA, Ramarao N, et al. Bacillus cereus-induced food-borne outbreaks in France, 2007 to 2014: epidemiology and genetic characterisation. Euro Surveill. 2016;21(48):30413.
Häggblom MM, Apetroaie C, Andersson MA, Salkinoja-Salonen MS. Quantitative analysis of cereulide, the emetic toxin of Bacillus cereus, produced under various conditions. Appl Environ Microbiol. 2002;68(5):2479–83.
Takeno A, Okamoto A, Tori K, Oshima K, Hirakawa H, Toh H, Agata N, Yamada K, Ogasawara N, Hayashi T, et al. Complete genome sequence of Bacillus cereus NC7401, which produces high levels of the emetic toxin cereulide. J Bacteriol. 2012;194(17):4767–8.
Kim JB, Kim JM, Cho SH, Oh HS, Choi NJ, Oh DH. Toxin genes profiles and toxin production ability of Bacillus cereus isolated from clinical and food samples. J Food Sci. 2011;76(1):T25–9.
Carter L, Chase HR, Gieseker CM, Hasbrouck NR, Stine CB, Khan A, Ewing-Peeples LJ, Tall BD, Gopinath GR. Analysis of enterotoxigenic Bacillus cereus strains from dried foods using whole genome sequencing, multi-locus sequence analysis and toxin gene prevalence and distribution using endpoint PCR analysis. Int J Food Microbiol. 2018;284:31–9.
Hansen BM, Hendriksen NB. Detection of enterotoxic Bacillus cereus and Bacillus thuringiensis strains by PCR analysis. Appl Environ Microbiol. 2001;67(1):185–9.
Saleh-Lakha S, Leon-Velarde CG, Chen S, Lee S, Shannon K, Fabri M, Downing G, Keown B. A study to assess the numbers and prevalence of Bacillus cereus and its toxins in pasteurized fluid milk. J Food Prot. 2017;80(7):1085–9.
Barbosa TM, Levy SB. The impact of antibiotic use on resistance development and persistence. Drug Resist Updat. 2000;3(5):303–11.
Savić D, Miljković-Selimović B, Lepšanović Z, Tambur Z, Konstantinović S, Stanković N, Ristanović E. Antimicrobial susceptibility and beta-lactamase production in Bacillus cereus isolates from stool of patients, food and environment samples. Vojnosanit Pregl. 2016;73(10):904–9.
Shawish R, Tarabees R. Prevalence and antimicrobial resistance of Bacillus cereus isolated from beef products in Egypt. Open Vet J. 2017;7(4):337–41.
Madonna MJ, Zhu YF, Lampen JO. Nucleotide sequence of the beta-lactamase I gene of Bacillus cereus strains 569/H and 5/B. Nucleic Acids Res. 1987;15(4):1877.
Hussain M, Pastor FI, Lampen JO. Cloning and sequencing of the blaZ gene encoding beta-lactamase III, a lipoprotein of Bacillus cereus 569/H. J Bacteriol. 1987;169(2):579–86.
Park YB, Kim JB, Shin SW, Kim JC, Cho SH, Lee BK, Ahn J, Kim JM, Oh DH. Prevalence, genetic diversity, and antibiotic susceptibility of Bacillus cereus strains isolated from rice and cereals collected in Korea. J Food Prot. 2009;72(3):612–7.
Gascoigne AD, Richards J, Gould K, Gibson GJ. Successful treatment of Bacillus cereus infection with ciprofloxacin. Thorax. 1991;46(3):220–1.
Kemmerly SA, Pankey GA. Oral ciprofloxacin therapy for Bacillus cereus wound infection and bacteremia. Clin Infect Dis. 1993;16(1):189.
Katsuya H, Takata T, Ishikawa T, Sasaki H, Ishitsuka K, Takamatsu Y, Tamura K. A patient with acute myeloid leukemia who developed fatal pneumonia caused by carbapenem-resistant Bacillus cereus. J Infect Chemother. 2009;15(1):39–41.
Lee N, Sun JM, Kwon KY, Kim HJ, Koo M, Chun HS. Genetic diversity, antimicrobial resistance, and toxigenic profiles of Bacillus cereus strains isolated from Sunsik. J Food Prot. 2012;75(2):225–30.
Luna VA, King DS, Gulledge J, Cannons AC, Amuso PT, 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(3):555–67.
Ikram S, Heikal A, Finke S, Hofgaard A, Rehman Y, Sabri AN, Økstad OA. Bacillus cereus biofilm formation on central venous catheters of hospitalised cardiac patients. Biofouling. 2019;35(2):204–16.
Liu Y, Lai Q, Du J, Shao Z. Genetic diversity and population structure of the Bacillus cereus group bacteria from diverse marine environments. Sci Rep. 2017;7(1):689.
Yang Y, Yu X, Zhan L, Chen J, Zhang Y, Zhang J, Chen H, Zhang Z, Zhang Y, Lu Y, et al. Multilocus sequence type profiles of Bacillus cereus isolates from infant formula in China. Food Microbiol. 2017;62:46–50.
The Hygiene Ministry of China. National food safety standard. Food microbiological examination: General guidelines. Beijing: The Hygiene Ministry of China; 2010.
The Hygiene Ministry of China. National food safety standard. Food microbiological examination: Bacillus cereus test. Beijing: The Hygiene Ministry of China. 2003.
The Hygiene Ministry of China. National food safety standard. Food microbiological examination: Bacillus cereus test. Beijing: The Hygiene Ministry of China; 2014.
Tallent, S. M., Rhodehamel, E. J., Harmon, S. M., and Bennett, R. W. (2012). Bacteriological analytical manual; methods for specific pathogens. U. S. Food and Drug Administration. Chapter 14, Bacillus cereus (updated November 2019). 2012. Available at: https://www.fda.gov/food/laboratory-methods-food/bam-bacillus-cereus. Accessed 16 April, 2020.
Hariram U, Labbé R. Spore prevalence and toxigenicity of Bacillus cereus and Bacillus thuringiensis isolates from U.S. retail spices. J. Food Prot. 2015;78(3):590–6.
Harper NM, Getty KJ, Schmidt KA, Nutsch AL, Linton RH. Comparing the mannitol-egg yolk-polymyxin agar plating method with the three-tube most-probable-number method for enumeration of Bacillus cereus spores in raw and high-temperature, short-time pasteurized milk. J Food Prot. 2011;74(3):461–4.
Fagerlund A, Ween O, Lund T, Hardy SP, Granum PE. Genetic and functional analysis of the cytK family of genes in Bacillus cereus. Microbiology. 2004;150:2689–97.
Oltuszak-Walczak E, Walczak P. PCR detection of cytK gene in Bacillus cereus group strains isolated from food samples. J Microbiol Methods. 2013;95(2):295–301.
The Clinical. And laboratory standards institute (CLSI). Performance standards for antimicrobial susceptibility testing; twentieth informational supplement. Approved standard-M100-S20. Wayne, PA: The Clinical and Laboratory Standards Institute; 2010.
Mohammed MJ, Marston CK, Popovic T, Weyant RS, Tenover FC. Antimicrobial susceptibility testing of Bacillus anthracis: comparison of results obtained by using the National Committee for clinical laboratory standards broth microdilution reference and Etest agar. Gradient diffusion methods. J Clin Microbiol. 2002;40(6):1902–7.
Osman KM, Kappell AD, Orabi A, Al-Maary KS, Mubarak AS, Dawoud TM, Hemeg HA, Moussa IMI, Hessain AM, Yousef HMY, et al. Poultry and beef meat as potential seedbeds for antimicrobial resistant enterotoxigenic Bacillus species: a materializing epidemiological and potential severe health hazard. Sci Rep. 2018;8(1):11600.
Park KM, Jeong M, Park KJ, Koo M. Prevalence, enterotoxin genes, and antibiotic resistance of Bacillus cereus isolated from raw vegetables in Korea. J Food Prot. 2018;81(10):1590–7.
Yusuf U, Kotwal SK, Gupta S, Ahmed T. Identification and antibiogram pattern of Bacillus cereus from the milk and milk products in and around Jammu region. Vet World. 2018;11(2):186–91.
Jolley KA, Bray JE, Maiden MCJ. Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications. Wellcome Open Res. 2018; 3:124.
Ribeiro-Gonçalves B, Francisco AP, Vaz C, Ramirez M, Carriço JA. PHYLOViZ online: web-based tool for visualization, phylogenetic inference, analysis and sharing of minimum spanning trees. Nucleic Acids Res. 2016;44(W1):W246–51.
The authors would like to thank the mlstoslo web database (www.mlstoslo.uio.no/) for the assistance of B. cereus MLST information of different strains.
The present work was financially supported by the Guangdong Technological Innovation Strategy of Special Funds (Key-Area of Research and Development Program of Guangdong Province, grant number 2018B020205003), National Key R&D Program of China (grant number 2018YFC1602500), National Natural Science Foundation of China (grant numbers 31730070 and 31701195), Science and Technology Program of Guangdong Province (grant number 2017A070702018), GDAS’ Special Project of Science and Technology Development (grant number 2017GDASCX-0201), and the 1000-Youth Elite Program and the 111 Project. The funding bodies were not involved in the design of the study and collection, analysis, interpretation of data, and writing of the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Prevalence of toxin genes in B. cereus isolated from edible fungi in China.
The diameters of inhibition zones in antimicrobial susceptibility test. Among the zone diameter interpretive criteria for the antimicrobials used, the minimum value for the antimicrobial resistance interpretive criteria is 10 mm. When there is no bacteriostatic zone or the bacteriostatic zone diameter is less than 8 mm, the results are directly recorded as “R”.
Results of antibiotic resistance test for B. cereus isolates in the study.
The MLST typing results, and the profiles of antibiotic resistance and toxin genes of all isolates. “*” represents a new ST.
Phylogenetic analysis of 247 B. cereus isolates from edible fungi and eight type strains based on the MLST typing results. A cut-off value of 43% similarity was applied to define the clusters. ATCC14579: B. cereus ATCC14579; DSM 2048: Bacillus mycoides DSM 2048; DSM 12442: Bacillus pseudomycoides DSM 12442; WSBC 10204: Bacillus weihenstephanensis WSBC 10204; ATCC4728: Bacillus anthracis ATCC4728; ATCC10792: Bacillus thuringiensis ATCC10792; NVH 391-98: Bacillus cytotoxicus NVH 391-98; NC7401: clinical emetic type strain B. cereus NC7401; F4810/72: clinical emetic type strain B. cereus F4810/72.
Sampling cities where the edible fungi were collected. Sampling plan design and this sampling map were all done by ourselves.
Primers used in this study.
Zone diameter interpretive criteria for the antimicrobials used in this study.
About this article
Cite this article
Liu, C., Yu, P., Yu, S. et al. Assessment and molecular characterization of Bacillus cereus isolated from edible fungi in China. BMC Microbiol 20, 310 (2020). https://doi.org/10.1186/s12866-020-01996-0
- Bacillus cereus
- Edible fungi
- Antibiotic resistance
- Genetic diversity