Next Article in Journal
Diagnosis and Treatment of Fungal Infections in Lung Transplant Recipients
Next Article in Special Issue
Secreted Extracellular Products of Flavobacterium covae as Potential Immunogenic Factors for Protection against Columnaris Disease in Channel Catfish (Ictalurus punctatus)
Previous Article in Journal
Identification of Trichomonas vaginalis 5-Nitroimidazole Resistance Targets
Previous Article in Special Issue
Development of Bioluminescent Virulent Aeromonas hydrophila for Understanding Pathogenicity
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Pathogens
Volume 12
Issue 5
10.3390/pathogens12050693
pathogens-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

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

by
Jia-Ming Tsai
1,2,
Hsin-Wei Kuo
3 and
Winton Cheng
1,*
1
Department of Aquaculture, National Pingtung University of Science and Technology, Pingtung 91201, Taiwan
2
Fish Doctor Veterinary Clinic, Pingtung 94042, Taiwan
3
General Research Service Center, National Pingtung University of Science and Technology, Pingtung 91201, Taiwan
*
Author to whom correspondence should be addressed.
Pathogens 2023, 12(5), 693; https://doi.org/10.3390/pathogens12050693
Submission received: 16 January 2023 / Revised: 19 April 2023 / Accepted: 8 May 2023 / Published: 10 May 2023
(This article belongs to the Special Issue Emerging Infections in Aquatic Animals)
Browse Figures
Versions Notes

Abstract

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

1. Introduction

The Chinese soft-shell turtle, Tryonyx sinensis is well-adapted to aquatic environments and is favored for commercial production because of its high nutritional and medicinal values. In Taiwan, the main pond-breeding area of turtles is located in Kaohsiung and Pingtung counties. In Taiwan, due to the favorable climate, the period for Chinese soft-shell turtle egg-laying lasts about 9 months. The collection and exportation of eggs to China has produced the most important economic income of this industry. Generally, disease occurs when turtles are cultured at high density as a result of overstocked ponds or when turtles are exposed to poor water quality. In recent years, bacterial disease has increased the difficulty of aquaculture management and has led to limited yields, often resulting in significant economic loss [1,2,3].
In 2011, several cases of Chinese soft-shell turtle disease originating from B. cereus infection were reported for the first time in the pond-cultured farms in Guangdong Province, China [4]. Soon after, similar cases also appeared in Taiwan [5]. In addition to B. cereus, B. thuringiensis also was isolated from diseased Chinese soft-shell turtles [6]. B. cereus group has become a newly emerging bacterial pathogen of Chinese soft-shell turtle in Taiwan. Up to the present date, recurrent epizootic outbreaks have been reported in several pond-cultured farms in Pingtung and Kaohsiung Counties [7]. The emerging epidemic has resulted in severe economic losses to the Chinese soft-shell turtle culture industry in recent years.
Bacteria within the Bacillus cereus group, also known as B. cereus sensu lato, are Gram-positive, spore-forming, facultatively anaerobic, rod-shaped bacteria, consisting of at least eight closely related species that are widespread in natural environments [8]. Moreover, nine novel species of the Bacillus cereus group were proposed, including Bacillus tropicus, Bacillus paranthracis, Bacillus pacificus, Bacillus albus, Bacillus mobilis, Bacillus luti, Bacillus proteolyticus, Bacillus nitratireducens, Bacillus paramycoides [9]. The most prominent members of this group include Bacillus anthracis, Bacillus cereus, Bacillus thuringiensis, and Bacillus mycoides. The group has played a critical role in many aspects of human life, such as terrorism [10], human health [11], agriculture [12], aquaculture [4,6].
B. cereus is an opportunistic pathogen and a well-known causative agent of foodborne bacterial intoxicants that are known to cause an emetic or diarrheal type of food-associated disease, which is generally self-limiting. Most isolates appear to be harmless [13].
Therefore, it was a dramatic turn when B. cereus G9241 was initially isolated from a Louisiana welder in 1994, and subsequently identified as the causative agent of severe anthrax-like pneumonia [14]. In previous study, anthrax lead by B. anthracis can induce the systemic effects of edema and severe bacteremia, with bacilli often noticed within blood vessels of the various vital organs [15]. The marked bacteremia are one of the key feature in histopathology [16]. However, a similar disease, such as anthrax, could be triggered by the emerging pathogens.
Since 2004, more and more newly emerging strains of atypical B. cereus have been shown to be associated with the anthrax-like disease in humans and mammals such as chimpanzee, apes, and bovine [14,17,18,19,20,21,22]. These atypical B. cereus isolates are capable of causing food-borne illnesses and can cause severe, even fatal, infections in non-immunocompromised humans and various animals [23]. As such, B. cereus G9241 was thought to differ from typical B. cereus isolates by the presence of unique plasmids. The strains harbor two large virulence plasmids, pBCXO1 and pBC210 (or pBCXO2).
Two components, anthrax toxin and an extracellular capsule, are required for full virulence. The plasmid called pBCXO1, which is the anthrax toxin source of B. cereus G9241, shares a high similarity with plasmid pXO1 of B. anthracis and consists of the toxin genes pagA, lef, and cya for anthrax toxin components, protective antigen (PA), lethal factor (LF), and edema factor (EF), respectively.
For extracellular capsules, the hyaluronic acid (HA) component of the capsule is encoded by the hasACB operon located on pBCXO1, and the exopolysaccharide (Bps) component is encoded by a nine-gene operon, bpsX-H, located on plasmid pBC210. This suspected zoonotic disease has attracted significant attention due to the atypical B. cereus causing an anthrax-like disease in mammals. Additionally, it is worth mentioning that Bacillus tropicus isolated from animals and humans also induced anthrax-like disease in 1941 [24].
On the basis of its highly analogous traits, many members of the genus Bacillus could be our suspect, causing the disease and high mortality in the Chinese soft-shell turtle industry. In the study, we estimated the analytical methods of the previous studies [6] and provided an insight into genetic diversity. Genome sequence-based classification was employed to determine taxonomic classification, including orthologous average nucleotide identity (OrthoANI) [25] and digital DNA-DNA hybridization (dDDH) [26].
However, as for non-mammals, the distribution of atypical B. cereus, which has been implicated in lethal disease, is unknown. In this study, we showed the phenotypical and molecular analysis of Bacillus species JMT isolates from Chinese soft-shell turtle with anthrax-like disease. Furthermore, we proposed the taxonomic change for several atypical B. cereus isolates.

2. Materials and Methods

2.1. Case Data Collection in Present Study

2.1.1. Case Profile of Chinese Soft-Shell Turtles

In this study, we collected bacterial isolates, temporarily designated as Bacillus species JMT, from diseased Chinese soft-shell turtles. From 2016 to 2020, the epidemic led by Bacillus spp. JMT among turtles cultured in Kaohsiung (K) and Pingtung (PT), two geographic districts in southern Taiwan, caused severe mortality in up to 100% of infected animal. A total of 32 cases infected by Bacillus spp. JMT were provided by breeders via purposive sampling in the subadult/adult populations of pond-cultured Chinese soft-shell turtle from different ponds. The sources of the Bacillus spp. JMT isolates used in the study are listed in Table 1.

2.1.2. Identification of Bacteremia by Histopathology

Isolates JMT107-17 to JMT109-32, were collected for histopathological analysis (Table 1). Vital organs including the heart, lungs, liver, spleen, kidneys, and GastroIntestinal-tract (GI-tract) were sampled and fixed in 10% neutral buffered formalin. Subsequently, organ tissues were processed with standard techniques for dehydration, clearing, embedding, sectioning, staining (both Gram’s stain and hematoxylin and eosin stain), embedding in paraffin, and sectioning. The bacteremia were identified by the presentence of bacilli within blood vessels in various tissues.

2.2. Bacteriology Test

2.2.1. Bacterial Isolation, Culture and Collection

Bacteria were isolated from the liver and spleen of Chinese soft-shell turtles by aseptic operation and cultured on trypticase soy agar supplemented with 5% goat blood (BAP) at room temperature (28 °C). The gray–white colonies with β-hemolytic properties can be seen by the naked eye within 24 h. Furthermore, B. thuringiensis ATCC 10792 and B. cereus ATCC 14579, purchased from Bioresource Collection and Research Center (Hsinchu County, Taiwan), were utilized as the negative positive controls, respectively, to compare with the collected wild strains.

2.2.2. Parasporal Body Visualization with Coomassie Blue Stain

A 0.5 McFarland turbidity standard of the sporulated colony from all isolates was transferred to the microscope slides, and the slides were air-dried before staining. The specimens were stained (0.133% Coomassie Blue stain in 50% acetic acid [27]), rinsed with distilled water, air-dried, and observed by microscope using a 100× oil immersion objective. The appearance of the parasporal body was clearly identified by the presence of abundant dark-blue stains and diamond-shaped objects.

2.2.3. Capsule Visualization with Liu’s Stain and Negative Stain

The spleen from the diseased Chinese soft-shell turtles was imprinted on a microscope slide and dried by air. The tissue stamp smear was completed using Liu’s Solutions including Liu’s A and B solutions (Baso company, Taiwan). Liu’s staining was performed according to the manufacturer’s protocol. Then, the smears were washed with tap water and air-dried prior to microscopic examination.
The cells from a single colony of each bacterial isolates were scattered into 100 μL of sterilized saline solution and 5 μL of the mixture was pipetted on to a microscope slide and subsequently covered with a coverslip. One drop of India ink (Becton, Dickinson and Company, Franklin Lakes, NJ, USA) was applied on the edge of the cover slide. The presence of encapsulated cells on the wet mount slide was examined microscopically at 100× oil-immersion magnification.

2.2.4. Plasmid DNA and Genomic DNA Extraction

Plasmid and genomic DNA templates were prepared from the 16 h culture grown on trypticase soy agar, supplemented with 5% sheep blood. Cells were resuspended in 200 μL of sterile MilliQ water. The DNA of each bacterial strain was extracted and purified using Blood and Tissue Genomic DNA Miniprep System and Mini Plus Plasmid DNA Extraction System (Viogene, New Taipei City, Taiwan) according to the manufacturer’s instructions. Finally, purified DNA was used as a template for subsequent PCR reaction and stored at −20 °C until use. Plasmid DNAwaas used for these primer sets, lef3/lef4, JMlef436f/JMlef436r, JMbpsB208f/JMbpsB208r, JMpagA318f/JMpagA318r, JMhasA700f/JMhasA700r, JMcya882f/JMcya882r; genomic DNA was used for the rest of the primer sets (Table 2).

2.2.5. Oligonucleotide Primer Design

16S rRNA and five sets of primers were used to differentiate between specific species among Bacillus spp. [29]. PCR assays targeting the groEL gene were used to amplify 400 bp of expectant product from Bacillus group species. The rest of the sets targeting the gryB yielded 253 bp amplicon, 299 bp amplicon, 475 bp amplicon, and 604 bp amplicon for B. anthracis, B. thuringiensis, B. cereus, and B. mycoides, respectively. Additionally, we further employed a transcriptional regulator (XRE) as the new biomarker of B. thuringiensis [30].
In this study, we designed several new primer sets for detecting unique virulence genes. Specific oligonucleotide primers were designed for hyaluronan synthase and exopolysaccharide biosynthesis protein based on the hasA and bpsB sequences, respectively (GenBank accession number AAEK01000036 and AAEK01000004); whereas the primers for the PA, EF, and LF were designed from the pagA, cya, and lef sequences, respectively (GenBank accession number CP009592). To confirm the results, another published primer set was undertaken to detect the lef gene [31]. The oligonucleotides of all primer sets in this study were synthesized by Tri-I Biotech, Inc. (New Taipei City, Taiwan). The sequences of the primers used in this study are shown in Table 2.

2.2.6. PCR Condition, Gene Cloning and Sequencing

For each reaction, 50 ng of DNA template and 10 pmole of each primer was added to AccuPower PCR Premix (Bioneer, Daejeon, Korea), consisting of 250 µM of each of the nucleotides dATP, dCTP, dGTP, and dTTP, PCR buffer (10 mM Tris—HCl, 30 mM KCl, 1.5 mM MgCl2), and 1 U of Top DNA polymerase. The reaction was added to a total volume of 20 µL with nuclease-free water. The reaction parameters were an initial denaturation at 94 °C for 5 min, 35 cycles of amplification with denaturation at 94 °C for 30 s, annealing at 55 °C for 30 s, an extension at 72 °C for 60 s, and final extension at 72 °C for 7 min. The PCR conditions could be slightly varied to adapt to each primer set.
The PCR reaction was performed in a FlexCycler2 Thermal Cyclers (Analytik Jena GmbH, Jena, Germany). Finally, the PCR products were analyzed by electrophoresis in 2% agarose gels stained with DNA-View (BIOTOOLS Co., Ltd., Jupiter, FL, USA) and immediately obrserved with UVP GelSolo (Analytik Jena GmbH, Germany) for analysis.
PCR products amplified by each primer set were cloned into the yT&A vector (Yeastern Biotech, Taipei, Taiwan) according to the manufacturer’s instructions and transformed into Escherichia coli DH—5α (Yeastern Biotech, Taipei, Taiwan). The recombinant plasmid DNA from the different clones was extracted separately and sequenced using an ABI3730 automatic sequencer (Tri-I Biotech, Inc., New Taipei City, Taiwan). Both forward and reverse strands were sequenced for each individual. For the purpose of ensuring sequence quality, both strands were sequenced as a crosscheck. Determined sequences were compared with the reference strains in GenBank.
The partial hasA fragment (GenBank accession number OP588641), partial bpsB fragment (GenBank accession number OP588640), partial pagA fragment (GenBank accession number OP583942), and partial cya fragment (GenBank accession number OP584757) served as the positive control. In addition, the lef positive control was kindly provided by the Institute of Medical Science, Tzu-Chi University, Taiwan [32].

2.3. Whole Genome Sequencing

2.3.1. Oxford Nanopore Technologies

Genomic DNA of the representative isolate (JM105-2) was extracted using Pure extraction QIAGEN Genomic-tip 20/G Kit (QIAGEN GmbH, Hilden, Germany). The gDNA extractions were evaluated for quantity and quality using Qubit (Life Technologies Corporation, Carlsbad, CA, USA), and Q6000B UV-Vis spectrophotometer (Quawell Technology Limited, Carlsbad, CA, USA), following the manufacturer’s instructions. Nanopore DNA Libraries were prepared using 2000 ng of gDNA per sample and preprocessed by Short Read Eliminator Kit (Circulomics Inc., Baltimore, MD, USA) to enhance long-read sequencing read lengths by removing short DNA.
Subsequently, libraries construction followed the native barcoding genomic DNA protocol (EXP-NBD114 and SQK-LSK109) (Oxford Nanopore Technologies, Oxford, UK) for combined FFPE repair and end-prep, ligation of barcodes, and the ligation of sequencing adapters. Sequencing was performed with MinKNOW software (version 21.05.8) on the GridION Mk1 using Nanopore sequencing platform by Oxford Nanopore Technologies (ONT) [33]. Subsequently, basecalling reads (Super-accurate basecalling model), and demultiplexing, adapter trimming was analyzed with Guppy (version 4.3.4) (Oxford Nanopore Technologies). De novo assembly was performed using SPAdes (version 3.14.1) [34]. The whole genome sequence of the representative isolate (JMT105-2) was deposited into GenBank (GenBank accession number CP115303, CP115304).
Finally, genome maps of the chromosome and plasmid of JMT105-2 isolate were drawn by Circos software [35] and CLCbio (CLC bio, Aarhus, Denmark), respectively.

2.3.2. Orthologous Average Nucleotide Identity and Digital DNA-DNA Hybridization

Orthologous average nucleotide identity (OrthoANI) between a pair of genomes was calculated by using CJ Bioscience’s Orthologous Average Nucleotide Identity Tool (OAT) [25].
Genome-to-genome distance (GGD) analysis was performed using GGD calculator (GGDC) with the GBDP3.0_BLASTPLUS program provided by German Collection of Microorganisms and Cell Cultures (DSMZ, Braunschweig, Germany) (http://ggdc.dsmz.de, accessed on 1 December 2022) [36]. Distances were calculated using recommended GBDP method, formula d4 [26]. These distances were transformed to values analogous to DNA-DNA hybridization (DDH) using a generalized linear model (GLM) and produced digital DDH (dDDH) values. The cutoff values of ANI and dDDH are, respectively, 95–96% and 70% [25,26]. If the value is above the cut-off point, it is recognized as the same species.
Whole genome sequences of the strains for OrthoANI and dDDH analyses were obtained from the GenBank database (https://www.ncbi.nlm.nih.gov/genome, accessed on 3 December 2022). These accession number included Bacillus tropicus JMT105-2 (GenBank accession number CP115303), Bacillus cereus G9241 (GenBank accession number NZ_CP009590), Bacillus cereus strain 03BB87 (GenBank accession number NZ_CP009941), Bacillus tropicus strain EMB20(GenBank accession number NZ_CP078081), Bacillus cereus ATCC 14579(GenBank accession number NC_004722), Bacillus thuringiensis serovar berliner ATCC 10792 (GenBank accession number NZ_CM000753), Bacillus anthracis str. Ames (GenBank accession number NC_003997), Bacillus mycoides DSM 2048 (GenBank accession number NZ_CM000742).

3. Results

3.1. Identification of Bacillus spp. JMT by Morphology and PCR

At first, the colony of Bacillus spp. JMT isolates was a white–grey color with a double zone of β-hemolysis in the absence of the rhizoid growth by macroscopic observation after a 24 h incubation period (Supplementary Figure S1).
The parasporal bodies were directly and strikingly marked by the presence of numerous dark blue staining objects. B. thuringiensis ATCC 10792 was perfectly demonstrated as a positive control (Figure 1a). Therefore, the dark blue staining bodies were not present in Bacillus spp. JMT in the microscopic field (Figure 1b).
The DNA of Bacillus spp. JMT isolates were analyzed using the 27f/1492r primer set, and the nucleotide sequence identity among the isolates was 100% after sequencing. The 16S rRNA fragment (1516 bp) of the JMT105-2 isolate was deposited in GenBank (GenBank accession numbers OP592263) and BLAST was conducted on NCBI database. The result we received shared the highest nucleotide sequence identity of 99~100% within three Bacillus species: B. anthracis, B. cereus, and B. thuringiensis. Therefore, the published multiplex PCR assay was undertaken to discriminate between the Bacillus relatives in the present study by recognizing specific DNA segments. We processed all bacterial isolates in this study, including B. cereus ATCC 14579, B. thuringiensis ATCC 10792 and all Bacillus spp. JMT isolates. One universal and one specific amplicon of each sample in electrophoresis were visible in the PCR reaction of the former two isolates. The rest of the Bacillus spp. JMT isolates were unanimous in presenting amplicons and agreed with B. thuringiensis ATCC 10792 pattern, which amplified the partial groEL fragment (400 bp) and the partial gyrB fragment (299 bp) (Supplementary Figure S2).
The transcriptional regulator (XRE) was described as a reliable and capable tool for identifying B. thuringiensis in a previous study [30]. The novel biomarker was used to estimate and confirm all bacterial strains. Both the standard strains used as control, B. cereus ATCC 14579, B. thuringiensis ATCC 10792, showed the expected results. Thus, all Bacillus spp. JMT isolates were negative and identical in the use of PCR targeting XRE (Supplementary Figure S2).

3.2. Capability of Capsule Production

The capsule-producing abilities of 32 Bacillus spp. JMT isolates from clinical sources were notably identical (Table 3). We performed a microscopical examination of the wet-mount slides after dropping India ink. One thick and clear layer surrounding the bacilli was observed, as shown in Figure 2a. The capsule was indicated by arrow signs. Negative staining of India ink highlights the location and thickness of the capsule produced by Bacillus spp. JMT. In addition, several heterophils with left-shifting were observed on tissue-imprinted smears with Liu’s stain. The capsule was visible as metachromatic color (magenta) and well-demarcated around the bacilli (Bacillus spp. JMT), as shown in Figure 2b. In comparison with encapsulated Bacillus spp. JMT, only bacilli were observed in B. cereus ATCC 14579 strains.

3.3. Anthrax-like Virulence Genes

3.3.1. PCR Assay

A total of five primer sets were used for screening in this study, specifically targeting each virulence gene of interest. In the capsule-associated genes, bpsB and hasA were found in all the Bacillus spp. JMT (32/32, 100%).
As for anthrax-like toxin-associated genes, pagA and cya were present (32/32, 100%) in all Bacillus spp. JMT isolates. Unexpectedly, the results of amplifying the partial lef gene were negative in all the Bacillus spp. JMT isolates using both newly designed and published primer sets. Whether capable of producing the unique capsule or an anthrax-like toxin, all Bacillus spp. JMT isolates retained identical properties. All results of the PCR assay are shown in Table 3. The following four nucleotide sequence documents obtained in the current study were submitted to GenBank. These sequences, derived from the JMT105-2 isolate, correspond to the partial hasA fragment (GenBank accession number OP588641), partial bpsB fragment (GenBank accession number OP588640), partial pagA fragment (GenBank accession number OP583942), and partial cya fragment (GenBank accession number OP584757).

3.3.2. The Absent Evidence of lef Gene via Nanopore Sequencing

A total of 1,442,375 reads (2,441,281,061 bp) were obtained with an average coverage of 236.8. The average read length was 1692.5 bp and the N50 value was 5,257,674 bp. The high-quality short-read and long-read sequences were De novo assembled into a complete chromosome and a plasmid using SPAdes Genome Assembler (v. 3.14.1) with default settings. The genome size of JMT105-2 is 5,487,300 bp, containing a chromosome (5,257,674 bp) and a plasmid (229,626 bp) (Supplementary Figure S3). No evidence of known lef gene sequences was present in these Nanopore sequencing data.

3.4. Bacteremia

By necropsy examination, the diseased Chinese soft-shell turtle exhibited some morphological changes, such as a generalized edema and darkened and enlarged spleen (Figure 3 and Figure 4). Further, these lesions correspond microscopically to moderate-to-severe edema depending on varying degrees of progress. The vital organs, the lung, liver, spleen, kidney, and GI tract showed diffused pathological changes characterized by inflammatory cells (mainly heterophils) infiltration and bacilli could be noted in blood vessels (Figure 4). Gram-positive bacilli had the same appearance as cytological smears and bacterial cultures. The result of the pathological analysis revealed a highly consistent presence of Bacillus spp. JMT-induced bacteremia among 16 clinical cases (16/16, 100%).

3.5. OrthoANI and GGD Analyses

Among the eight Bacillus spp. strains, the value of seven paired-genome combinations was higher than both thresholds in terms of orthologous average nucleotide identity and digital DNA–DNA hybridization similarities. The calculation recognized as the same species consisted of the following combinations: 1, Bacillus tropicus JMT105-2(CP115303) and Bacillus cereus G9241(NZ_CP009590.1); 2, Bacillus tropicus JMT105-2(CP115303) and Bacillus cereus strain 03BB87(NZ_CP009941); 3, Bacillus cereus G9241(NZ_CP009590.1) and Bacillus cereus strain 03BB87(NZ_CP009941); 4, Bacillus tropicus JMT105-2(CP115303) and Bacillus tropicus strain EMB20(NZ_CP078081.1); 5, Bacillus cereus G9241(NZ_CP009590.1) and Bacillus tropicus strain EMB20(NZ_CP078081.1); 6, Bacillus cereus strain 03BB87(NZ_CP009941) and Bacillus tropicus JMT105-2(CP115303); 7, Bacillus cereus ATCC 14579(NC_004722) and Bacillus thuringiensis serovar berliner ATCC 10792(NZ_CM000753) (listed in Table 4).

4. Discussion

Aquaculture production is essential to the economy and food supply of Taiwan. With the increasing demand for rich protein and high profitability, the aquaculture industry has been motivated to develop in all aspects. The Chinese soft-shell turtle, Tryonyx sinensis, is one of the most common and important commercially cultured aquaculture species raised in Taiwan and Southeast Asia. The disease induced by Bacillus spp. JMT has had a significant adverse impact on the industry. However, we have known so little about it that we struggle with taxonomic classification. Several methods were developed to deal with this situation.

4.1. JMT Isolates Belong to Bacillus tropicus

The results demonstrated that the etiological agent could be B. tropicus. According to phenotypic and molecular analysis, Bacillus spp. JMT isolates belong to B. tropicus and share the same pattern. In the process of identification, we used various methods to rule out and narrow down possible species. The colony morphology of Bacillus spp. JMT indicated a low probability of it being B. anthracis or B. mycoid. B. anthacis, as it is nonhemolytic while B. mycoid forms long hair or root-like colonies as rhizoid growth. Thus, the visualization of the parasporal bodies was critical in differentiating B. thuringiensis from the B. cereus group. The presentation of the crystal protein could be visually discerned by modified Coomassie staining [27]. The assay provided a reliable result in identifying Bacillus spp. JMT isolates and led to a definite result; it did not belong to B. thuringiensis in terms of phenotypic properties. Nonetheless, Bacillus spp. JMT and B. cereus 14,579 showed slightly different colony morphologies in terms of texture. The former colony appeared smoother than the latter.
Following molecular analysis, the 16S rRNA fragment was indistinct and not beneficial in distinguishing between species due to its highly conservative DNA sequence. Furthermore, we used a similar strategy as the previous study [6], employing multiplex PCR assay for identification. The result unexpectedly indicated that Bacillus spp. JMT belonged to B. thuringiensis. After that, the PCR assay targeting the molecular biomarker (XRE) and the multiplex PCR assay results did not match, either.
The major difference between B. cereus and B. thuringiensis is that the latter possesses a unique ability to synthesize parasporal crystalline inclusions. Surprisingly, the multiplex PCR primer sets for discriminating in this study apparently failed to identify Bacillus species and cannot address Bacillus spp. JMT isolates.
Interestingly, the multiplex PCR results of the study by Chen et al. [6] were the same as ours. Moreover, they claimed that B. thuringiensis was the pathogen and highly pathogenic to juvenile Chinese soft-shell turtles. Nevertheless, Bacillus spp. JMT did not possess the essential ability to produce B. thuringiensis crystals in our study.
The conflict between the multiplex PCR results and crystal production ability of the assay requires further high-resolution analysis. Genome-sequence-based analyses provided new insight into the similarity between Bacillus species. The digital DNA–DNA hybridization and orthologous average nucleotide identity consequently became a standard technique for the circumscription of bacterial species; we can obtain relatively reliable information to determine genomic similarities. An advantage of these genome-based methods is that they often form sharp, clear clusters of strains compared to those only circumscribed by phenotypic traits [27]. With in silico analyses, the valid published atypical B. cereus (strain G9241 and 03BB87) and Bacillus spp. JMT105-2 isolate were suggested to be the same species as B. tropicus strain EMB20, not B. cereus ATCC 14579 or Bacillus thuringiensis serovar berliner ATCC 10792. It is worth mentioning that although B. cereus and B. thuringiensis share genomic similarities, they were still considered to be different species [37]. On the basis of the morphology and PCR data, along with the high OrthoANI and dDDH values of Bacillus spp. JMT and B. tropicus strain EMB20, these isolates from diseased Chinese shell-turtle were assigned as Bacillus tropicus JMT.
As previously mentioned, we solved a confusing problem based on the bacterial isolates collected in the study. The members of the genus Bacillus are too close to be distinguished. In previous studies [6,7,38], B. cereus and B. thuringiensis were considered the etiological agents of Chinese soft-shell turtles. Due to the lack of genome sequence-based analyses, it was difficult to address the taxonomic recognition of these close species-level lineages.
In Taiwan, many Chinese soft-shell turtle B. cereus-associated researchers unofficially posit the following theory: the origin of the etiological agent is a biological agricultural pesticide. The typical pond-cultured farms are usually located by country roads and surrounded by agricultural farms in the vicinity. It seems rational to assume that B. thuringiensis might be an agent; therefore, it became the most likely suspect. B. thuringiensis is deadly to certain herbivorous insects through the production of crystal toxins and is the most frequently used biological pesticide globally.
In this study, all Bacillus spp. JMT isolates provided identical results in several tests. Regardless, the etiological agent that triggered the mass death during outbreaks is more likely to be B. tropicus than either B. cereus or B. thuringiensis in Taiwan. According to the summarized results, Bacillus spp. JMT is referred to as Bacillus tropicus JMT in the following section of the study.

4.2. Genome Analysis-Based Reclassification of Atypical Bacillus cereus as Bacillus tropicus

B. tropicus was previously recognized as atypical B. cereus. In this study, we reclarified the taxonomic status of B. tropicus JMT, as well as these atypical B. cereus strains.
B. tropicus JMT retains several anthrax-like virulence genes and displays different characteristics, such as B. cereus G9241 and B. cereus 03BB87, not the non-pathogenic B. tropicus strain [39,40,41]. In previous studies, these atypical B. cereus strains have been widely accepted to overcome host immunity and cause severe symptoms, even death. Since B. cereus G9241 has been identified and characterized, the stereotypes of pathogenic B. cereus isolates not only originate from food-borne incidents but are also related to severe pulmonary disorders. In addition to being an opportunistic pathogen, it is an etiologic agent of emerging disease, requiring vigilance regarding public health issues. These atypical B. cereus strains carry unique virulence toxins encoded on the megaplasmids [14,17,18,20,42,43]. These atypical B. cereus plasmids have an analogous function to the pXO1 and pXO2 plasmids of B. anthracis and have proven to be sources of pathogenicity [44]. Previous studies have emphasized the diversity between typical and atypical strains in terms of pathogenicity and phenotypic traits. Surprisingly, the genome-sequence-based analyses indicated that we should verify the taxonomic status of all atypical B. cereus-associated strains which induce anthrax-like disease.

4.3. The Anthrax-like Property of B. tropicus JMT Isolated from Chinese Soft-Shell Turtle

The infectious mechanism of B. tropicus JMT isolates could be similar to its relative, B. anthracis, which causes a famous zoonotic disease. Anthrax is an acute bacterial disease mainly occurring in wild and domestic herbivores and it occasionally causes acute infection in humans. Like anthrax, atypical B. cereus strains were reported to mainly be isolated from mammals. In the present study, we conducted a retrospective screening of the isolates from 2016 to 2020. Our findings provide solid evidence of the presence of B. tropicus JMT and its critical pathogenic ability against a reptile: Chinese soft-shell turtle. The capsule could escape from phagocytosis via the host immune system, so this has been considered a virulence factor. For instance, the commonly used laboratory strain B. anthracis Sterne harbors plasmid pXO1 only; consequently, its pathogenicity is reduced. Nevertheless, the encapsulated wild B. anthracis strain can inhibit both complement-dependent and complement-independent opsonic phagocytosis. [45]. Our findings demonstrated that the capsule production ability was prevalent among the B. tropicus JMT isolates. The thick and unstainable coat around the bacilli was visualized by microscopy, using India ink staining in all isolates. The capsule was also visible on the smear with Liu’s staining, which has been used for clinical aquatic animal diagnosis in Taiwan. Its application as a capsule stain with excellent discrimination is similar to the standard staining method used for B. anthracis [46]. This characteristic was in agreement with the PCR results of this study and previous atypical strains.

4.4. Anthrax-like Disease

The anthrax-like diseases induced by atypical B. cereus were usually severe to lethal in mammals in previous studies. We also observed a similar situation in the reptile, Chinese soft-shell turtle. Since 2016, the emerging pathogen has led to a great epidemic of diseases. Bacteremia was recognizable in all vital organs sampled for pathological diagnosis in this study. Previous studies have suggested that the disease may be related to emetic toxin enterotoxins secreted by B. cereus [7]. Nevertheless, the classic symptoms we observed in this study were a similar systemic infectious disease rather than an enteric one.
The causative agent, B. tropicus JMT, induced similar severe systemic pathological changes in all sampled specimens in this study. The data showed that the diseased animals in this study would suffer from a severe systemic reaction and the circulatory system would eventually be invaded by B. tropicus JMT in the terminal stage. In the aquaculture field, the individuals engage the pathogen, inducing a severe serial immune response in the diseased animals. The diseased individuals usually showed lethargic swimming behavior on water surfaces with twisted necks and rigid legs. Further, they exhibited reduced alertness to the surrounding environment. With the course of disease development, abnormal behaviors were observed, accompanied by prominent reductions in feed intake reduction and general edema of the tissues (Figure 3.).
The various cumulative mortality rate of each farm was dependent on the window of opportunity for appropriate antibiotic treatment at the onset of the outbreak. With the progressively declining health condition, the prognosis of the animals acquiring B. tropicus JMT infection became increasingly worse. In later stages of outbreak, the cumulative probability of mortality reached 100%. Therefore, the probability of a widespread disaster is much higher in the absence of medical treatment.

4.5. Spore Contamination

The environmental pond was contaminated by the sporulated forms after the disease outbreak. Medical treatment with various antibiotics could be effective on the suffering animals. Unfortunately, the rest of animals in the pond were not exempt from invasion of the disease caused by the spores.
The carnivorous habits of the soft-shell turtles also exacerbate the spread of infectious diseases. For animals in the pond, dying animals and the corpses infected with B. tropicus JMT were hotspot sources of infection during the outbreak and threatened healthy individuals with overwhelming exposure; in the cultured environment, they grew and sporulated. Subsequently, the sporulated forms shed and were spread throughout the cultured pond by dying animals. When the external conditions were appropriate, the spores germinated and multiplied, increasing the risk of transmission to other healthy individuals in the same pond. The spores are relatively resistant to heat, dryness, and disinfectants. Consequently, the contamination will last as long as there is no way to resist the disease. As a result, healthy individuals acquired the disease annually from batch to batch in the pond contaminated by the spores of B. tropicus JMT.
Atypical B. cereus isolates rarely inhabit natural soils in the environment [47]. As there is an increasing number contaminated farming areas in Taiwan, the land used to farm Chinese soft-shell turtles has gradually been lost. This concern is consistent with our clinical observations. Owing to the difficulty of spore clearance, B. tropicus JMT imposes an adverse economic impact by dramatically decreasing production. In the meantime, a useful decontamination program is desperately needed to control the challenges we will face in the future.

4.6. Absence of lef Gene in PCR Assay and Nanopore Sequencing

The lef gene was absent in this study. Five newly designed primer sets and one published primer set for screening virulence toxins were employed to confirm the presence of the megaplasmid. The four primer sets were related to the anthrax tripartite toxin; the rest targeted bpxB and hasA, corresponding to capsule production. Four of the six primer sets amplified the expected PCR products among all the B. tropicus JMT isolates. The PCR results imply that B. tropicus JMT carries the plasmid that plays a critical role in pathogenicity. Unexpectedly, neither the newly designed nor published primer sets for targeting lef could amplify the PCR product. To the best of our knowledge, the genes of lef, pagA, and cya (anthrax tripartite toxin) are encoded on the same plasmid in atypical B. cereus and B. anthracis. Sequence mutations or incomplete plasmids in B. tropicus JMT probably result in amplification failure. The phenomenon of carrying a partial plasmid occurred in B. cereus FL2013 [48].
Furthermore, whole-genome sequencing is a powerful tool to provide a high-resolution view of the entire genome. The data of the de novo assembled genome showed no evidence of the segment of lef gene. The preliminary result in this study reveals the absence of the lef gene.

4.7. Remarkable Capsule-Production Capacity of B. tropicus JMT

B. tropicus JMT can synthesize astonishingly thick capsules. Both staining visualization methods demonstrated that B. tropicus JMT isolates possess an inherent ability to produce capsules.
As for the capsule production, the PCR results correspond to the above-mentioned phenotypic characteristics of B. tropicus JMT. The hasACB genes are encoded on the plasmid of B. anthrax and atypical B. cereus. In contrast to B. anthracis, atypical B. cereus contains non-mutated hasACB genes and is functional for capsule production [49]. The capsule produced by B. anthracis is considered to help avoid opsonization and phagocytosis of the host immune system [50,51]. Similarly, the atypical B. cereus possess this ability due to their two unique natural gifts. Both the virulence plasmids make the required contributions for synthesizing a unique capsule structured with hyaluronic acid (hasACB) and polysaccharide (bpsX-H). Thus, they differ from the poly-γ-d-glutamic acid (PGA) capsule produced by B. anthracis.
In terms of pathogenicity and the thickness of the produced capsule, hyaluronic acid is thought to be a more influential capsule componant than exopolysaccharide in B. cereus G9241 [44]. B. tropicus JMT isolates carry both the hasA and bpsB genes related to the formation of capsules. According to the PCR results, capsule-associated genes were present in all B. tropicus JMT isolates, implying that they could possibly provide ancillary pathogenic factors. However, whether they have an analogous effect on capsule production compared to previous atypical B. cereus isolates remains unknown.

4.8. A Biomarker Tool for Rapid Detecting B. tropicus JMT

The powerful primer sets used in the present study are vital for detecting pathogenic B. tropicus JMT isolates. B. cereus group members are too similar to be recognized as B. tropicus JMT in phenotypic analysis. Despite its high pathogenicity and infectivity against Chinese soft-shell turtle, there is no evidence that any other animals or humans could naturally acquire the disease and no medical cases have been documented to date. The living organisms have varying innate susceptibility to the pathogen. Generally, individual species and genetic background are critical determinants of pathogenicity [52]. Using the One Health basics, we recognize that the health of people, animals, and ecosystems is closely connected. Moreover, the spread of many emerging infectious diseases among various animals is also linked and interdependent. B. tropicus JMT carries anthrax-like virulence genes, which inevitably give rise to concerns about zoonotic infectious diseases. A thorough understanding of the host range and geographic distribution of B. tropicus JMT is essential for the control and prevention of disease outbreaks. Based on previous studies, the distribution of atypical B. cereus, which usually infects mammals, has been restricted to North America and Africa. This study is the first retrospective screening report of anthrax-like virulent genes in B. tropicus (atypical B. cereus) isolated from a reptile and reveals that the atypical strains are pathogens of not only mammals but also reptiles. There is an imperceptible threat that is expanding to affect human and animal health. Additionally, the geographic distribution and host range of B. tropicus were underestimated. Extensive surveillance will allow us to better understand and evaluate the epidemic cross-species risk. Furthermore, it will help raise alertness and develop strategies to prevent potential threats. Likewise, for the uphill challenges faced by the aquaculture industry, precise and rapid identification of the etiologic agent is also increasingly important to recover economic losses. Consequently, a reliable tool is needed to assess potentially contagious hazards.
The amplification of each specific virulence gene by PCR can be a valuable tool to distinguish and identify pathogenic B. tropicus JMT. The powerful PCR assay based on the anthrax-like gene primer sets JMbpsB208f/JMbpsB208r, JMpagA318f/JMpagA318r, JMhasA700f/JMhasA700r, and JMcya882f/JMcya882r displayed a much stronger ability to identify B. tropicus JMT isolates, with the highest detection rate (32/32 = 100%). These results help to confirm that this method is likely to be a viable candidate for screening the causative agents.

5. Conclusions

In this study, the bacterial isolates collected from Chinese soft-shell turtles cultured in southern Taiwan during the epizootics of 2016–2020 were B. tropicus str. JMT in terms of taxonomic classification. They harbor anthrax-like virulence genes, as shown using a PCR assay (32/32 = 100%). Furthermore, genome sequence-based analyses with high resolution can successfully unveil the lurking emerging pathogen in turtle aquaculture.
Our findings demonstrate that the geographic distribution of atypical B. cereus is not only in North America and Africa, but also in Asia. Consequently, the clinical observations imply that B. tropicus JMT is a frank pathogen against its host, Chinese soft-shell turtle. Lastly, further studies are needed to elucidate B. tropicus JMT’s pathogenic ability against other animals.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pathogens12050693/s1. Figure S1. The colony morphology of Bacillus tropicus JMT; Figure S2. Gels of PCR assays in this study; Figure S3. Genome map of JMT105-2 isolate contains a chromosome and a plasmid.

Author Contributions

Conceptualization, W.C. and J.-M.T.; methodology, W.C. and J.-M.T.; resources, W.C. and H.-W.K.; data curation, J.-M.T.; writing—original draft preparation, J.-M.T.; writing—review and editing, W.C. and H.-W.K.; supervision, W.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partly funded by the Taiwan Fisheries Agency of Fund, grant number 111-Adjust-Fishery-01(06).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors greatly acknowledge Der-Shan Sun for sharing precious material generously (Institute of Medical Science, Tzu-Chi University, Taiwan, ROC) and Feng-Chia Hsieh for providing advices (Taiwan Agricultural Chemicals and Toxic Substances Research Institute).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, J.; Zhu, N.; Kong, L.; Bei, Y.; Zheng, T.; Ding, X.; He, Z. First Case of Soft Shell Disease in Chinese Soft-shelled Turtle (Trionyx sinens) Associated with Aeromonas sobria–A. veronii complex. Aquaculture 2013, 406, 62–67. [Google Scholar] [CrossRef]
  2. Hsieh, C.; Chang, T.; Shen, Y.; Chang, C.; Tu, C.; Tung, M.; Chen, L.; Tsai, S. Pathological and PCR Detection of Mycobacteriosis in Pond-cultured Chinese Soft Shell Turtles, Trionyx sinensis. Aquaculture 2006, 261, 10–16. [Google Scholar] [CrossRef]
  3. Jensen, H. Mycobacterium kansasii Infection in Chinese Soft Shell Turtle (Pelodiscus sinensis). Vet. Rec. 2003, 152, 474–476. [Google Scholar]
  4. Tan, A.-P.; Zhao, F.; Jiang, L.; Luo, L.; Wang, W.-L.; Peng, H.-L.; Chen, Y.-L.; Zou, W.-M. Isolation and Identification of Bacillus cereus from Trionyx sinensis. Guangdong Agric. Sci. 2011, 20, 115–119. [Google Scholar]
  5. Yeh, S.-H. Pathological and Molecular Studies on Bacillus spp. of Soft-Shell Turtle (Trionyx sinensis). Master Thesis, National Pingtung University of Science and Technology, Pingtung, Taiwan, 2013. [Google Scholar]
  6. Chen, J.; Zhu, N.; Kong, L.; Bei, Y.; Zheng, T.; Ding, X.; He, Z. First Reported Fatal Bacillus thuringiensis Infections in Chinese Soft-shelled Turtles (Trionyx sinensis). Aquaculture 2014, 428, 16–20. [Google Scholar] [CrossRef]
  7. Cheng, L.W.; Rao, S.; Poudyal, S.; Wang, P.C.; Chen, S.C. Genotype and Virulence Gene Analyses of Bacillus cereus Group Clinical Isolates from the Chinese Softshell Turtle (Pelodiscus sinensis) in Taiwan. J. Fish Dis. 2021, 44, 1515–1529. [Google Scholar] [CrossRef]
  8. Liu, Y.; Lai, Q.; Göker, M.; Meier-Kolthoff, J.P.; Wang, M.; Sun, Y.; Wang, L.; Shao, Z. Genomic Insights into the Taxonomic Status of the Bacillus cereus group. Sci. Rep. 2015, 5, 14082. [Google Scholar] [CrossRef] [PubMed]
  9. 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]
  10. Inglesby, T.V.; O’Toole, T.; Henderson, D.A.; Bartlett, J.G.; Ascher, M.S.; Eitzen, E.; Friedlander, A.M.; Gerberding, J.; Hauer, J.; Hughes, J. Anthrax as a Biological Weapon, 2002: Updated Recommendations for Management. Jama 2002, 287, 2236–2252. [Google Scholar] [CrossRef]
  11. Goepfert, J.; Spira, W.; Kim, H. Bacillus cereus: Food Poisoning Organism. A Review. J. Milk Food Technol. 1972, 35, 213–227. [Google Scholar] [CrossRef]
  12. Höfte, H.; Whiteley, H. Insecticidal Crystal Proteins of Bacillus thuringiensis. Microbiol. Mol. Biol. Rev. 1989, 53, 242–255. [Google Scholar] [CrossRef] [PubMed]
  13. Tewari, A.; Abdullah, S. Bacillus cereus Food Poisoning: International and Indian Perspective. J. Food Sci. Technol. 2015, 52, 2500–2511. [Google Scholar] [CrossRef]
  14. Hoffmaster, A.R.; Ravel, J.; Rasko, D.A.; Chapman, G.D.; Chute, M.D.; Marston, C.K.; De, B.K.; Sacchi, C.T.; Fitzgerald, C.; Mayer, L.W.; et al. Identification of Anthrax Toxin Genes in a Bacillus cereus Associated with an Illness Resembling Inhalation Anthrax. Proc. Natl. Acad. Sci. USA 2004, 101, 8449–8454. [Google Scholar] [CrossRef] [PubMed]
  15. Firoved, A.M.; Miller, G.F.; Moayeri, M.; Kakkar, R.; Shen, Y.; Wiggins, J.F.; McNally, E.M.; Tang, W.-J.; Leppla, S.H. Bacillus anthracis Edema Toxin Causes Extensive Tissue Lesions and Rapid Lethality in Mice. Am. J. Pathol. 2005, 167, 1309–1320. [Google Scholar] [CrossRef]
  16. Abramova, F.A.; Grinberg, L.M.; Yampolskaya, O.V.; Walker, D.H. Pathology of Inhalational Anthrax in 42 Cases from the Sverdlovsk Outbreak of 1979. Proc. Natl. Acad. Sci. USA 1993, 90, 2291–2294. [Google Scholar] [CrossRef]
  17. Avashia, S.B.; Riggins, W.; Lindley, C.; Hoffmaster, A.; Drumgoole, R.; Nekomoto, T.; Jackson, P.J.; Hill, K.K.; Williams, K.; Lehman, L. Fatal Pneumonia among Metalworkers Due to Inhalation Exposure to Bacillus cereus Containing Bacillus anthracis Toxin Genes. Clin. Infect. Dis. 2007, 44, 414–416. [Google Scholar] [CrossRef] [PubMed]
  18. Hoffmaster, A.R.; Hill, K.K.; Gee, J.E.; Marston, C.K.; De, B.K.; Popovic, T.; Sue, D.; Wilkins, P.P.; Avashia, S.B.; Drumgoole, R. Characterization of Bacillus cereus Isolates Associated with Fatal Pneumonias: Strains Are Closely Related to Bacillus anthracis and Harbor B. anthracis Virulence Genes. J. Clin. Microbiol. 2006, 44, 3352–3360. [Google Scholar] [CrossRef]
  19. Klee, S.R.; Brzuszkiewicz, E.B.; Nattermann, H.; Bruggemann, H.; Dupke, S.; Wollherr, A.; Franz, T.; Pauli, G.; Appel, B.; Liebl, W.; et al. The Genome of a Bacillus Isolate Causing Anthrax in Chimpanzees Combines Chromosomal Properties of B. cereus with B. anthracis Virulence Plasmids. PLoS ONE 2010, 5, e10986. [Google Scholar] [CrossRef]
  20. Klee, S.R.; Ozel, M.; Appel, B.; Boesch, C.; Ellerbrok, H.; Jacob, D.; Holland, G.; Leendertz, F.H.; Pauli, G.; Grunow, R.; et al. Characterization of Bacillus anthracis-like Bacteria Isolated from Wild Great Apes from Cote d’Ivoire and Cameroon. J. Bacteriol. 2006, 188, 5333–5344. [Google Scholar] [CrossRef]
  21. Leendertz, F.H.; Ellerbrok, H.; Boesch, C.; Couacy-Hymann, E.; Mätz-Rensing, K.; Hakenbeck, R.; Bergmann, C.; Abaza, P.; Junglen, S.; Moebius, Y. Anthrax Kills Wild Chimpanzees in a Tropical Rainforest. Nature 2004, 430, 451–452. [Google Scholar] [CrossRef]
  22. Pilo, P.; Rossano, A.; Bamamga, H.; Abdoulkadiri, S.; Perreten, V.; Frey, J. Bovine Bacillus anthracis in Cameroon. Appl. Environ. Microbiol. 2011, 77, 5818–5821. [Google Scholar] [CrossRef] [PubMed]
  23. Baldwin, V.M. You Can’t B. cereus–A Review of Bacillus cereus Strains That Cause Anthrax-like Disease. Front. Microbiol. 2020, 11, 1731. [Google Scholar] [CrossRef] [PubMed]
  24. Heaslip, W. Bacillus tropicus, a New Species Isolated from Man and Animals Described and Compared with Other Bacilli Resembling, Bacillus anthracis. Med. J. 1941, 2, 536–540. [Google Scholar] [CrossRef]
  25. Lee, I.; Kim, Y.O.; Park, S.-C.; Chun, J. OrthoANI: An Improved Algorithm and Software for Calculating Average Nucleotide Identity. Int. J. Syst. Evol. Microbiol. 2016, 66, 1100–1103. [Google Scholar] [CrossRef]
  26. Meier-Kolthoff, J.P.; Auch, A.F.; Klenk, H.-P.; Göker, M. Genome Sequence-based Species Delimitation with Confidence Intervals and Improved Distance Functions. BMC Bioinform. 2013, 14, 60. [Google Scholar] [CrossRef]
  27. Ammons, D.; Rampersad, J.; Khan, A. Usefulness of Staining Parasporal Bodies When Screening for Bacillus thuringiensis. J. Invertebr. Pathol. 2002, 79, 203–204. [Google Scholar] [CrossRef]
  28. Lane, D.J. 16S/23S rRNA sequencing. In Nucleic Acid Techniques in Bacterial Systematic; Goodfellow, E.S.A.M., Ed.; John Wiley and Sons: New York, NY, USA, 1991; pp. 115–175. [Google Scholar]
  29. Park, S.-H.; Kim, H.-J.; Kim, J.-H.; Kim, T.-W.; Kim, H.-Y. Simultaneous Detection and Identification of Bacillus cereus Group Bacteria Using Multiplex PCR. J. Microbiol. Biotechnol. 2007, 17, 1177–1182. [Google Scholar]
  30. Wei, S.; Chelliah, R.; Park, B.-J.; Kim, S.-H.; Forghani, F.; Cho, M.S.; Park, D.-S.; Jin, Y.-G.; Oh, D.-H. Differentiation of Bacillus thuringiensis From Bacillus cereus Group Using a Unique Marker Based on Real-Time PCR. Front. Microbiol. 2019, 10, 883. [Google Scholar] [CrossRef]
  31. Ramisse, V.; Patra, G.; Garrigue, H.; Guesdon, J.-L.; Mock, M. Identification and Characterization of Bacillus anthracis by Multiplex PCR Analysis of Sequences on Plasmids pXO1 and pXO2 and Chromosomal DNA. FEMS Microbiol. Lett. 1996, 145, 9–16. [Google Scholar] [CrossRef]
  32. Chang, H.-H.; Tsai, M.-F.; Chung, C.-P.; Chen, P.-K.; Hu, H.-I.; Kau, J.-H.; Huang, H.-H.; Lin, H.-C.; Sun, D.-S. Single-step Purification of Recombinant Anthrax Lethal Factor from Periplasm of Escherichia coli. J. Biotechnol. 2006, 126, 277–285. [Google Scholar] [CrossRef]
  33. Jain, M.; Olsen, H.E.; Paten, B.; Akeson, M. The Oxford Nanopore MinION: Delivery of Nanopore Sequencing to the Genomics Community. Genome Biol. 2016, 17, 239. [Google Scholar] [CrossRef] [PubMed]
  34. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D. SPAdes: A New Genome Assembly Algorithm and Its Applications to Single-cell Sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef] [PubMed]
  35. Krzywinski, M.; Schein, J.; Birol, I.; Connors, J.; Gascoyne, R.; Horsman, D.; Jones, S.J.; Marra, M.A. Circos: An Information Aesthetic for Comparative Genomics. Genome Res. 2009, 19, 1639–1645. [Google Scholar] [CrossRef] [PubMed]
  36. Meier-Kolthoff, J.P.; Carbasse, J.S.; Peinado-Olarte, R.L.; Göker, M. TYGS and LPSN: A Database Tandem for Fast and Reliable Genome-based Classification and Nomenclature of Prokaryotes. Nucleic Acids Res. 2022, 50, D801–D807. [Google Scholar] [CrossRef] [PubMed]
  37. Baek, I.; Lee, K.; Goodfellow, M.; Chun, J. Comparative Genomic and Phylogenomic Analyses Clarify Relationships within and between Bacillus cereus and Bacillus thuringiensis: Proposal for the Recognition of Two Bacillus thuringiensis Genomovars. Front. Microbiol. 2019, 10, 1978. [Google Scholar] [CrossRef] [PubMed]
  38. Meng, Q.; Yin, F.; Fu, C.; Chen, F.; Liu, C.; Yuan, N.; Wang, L.; Zhang, H.; Qian, D. Isolation, Identification and Pathogenicity Analysis of Bacillus cereus from Chinese Soft-shelled Turtles, Pelodiscus sinensis. Acta Hydrobiol. Sin. 2019, 43, 570–578. [Google Scholar]
  39. Samanta, S.; Datta, D.; Halder, G. Biodegradation Efficacy of Soil Inherent Novel sp. Bacillus tropicus (MK318648) onto Low Density Polyethylene Matrix. J. Polym. Res. 2020, 27, 324. [Google Scholar] [CrossRef]
  40. Thakur, P.; Singh, A.K.; Mukherjee, G. Isolation and Characterization of Alkaline Pectinase Productive Bacillus tropicus from Fruit and Vegetable Waste Dump Soil. Braz. Arch. Biol. Technol. 2021, 64, 319. [Google Scholar] [CrossRef]
  41. Uddin, M.; Swathi, K.; Anil, A.; Boopathy, R.; Ramani, K.; Sekaran, G. Biosequestration of Lignin in Municipal Landfill Leachate by Tailored Cationic Lipoprotein Biosurfactant through Bacillus tropicus Valorized Tannery Solid Waste. J. Environ. Manag. 2021, 300, 113755. [Google Scholar] [CrossRef]
  42. Brezillon, C.; Haustant, M.; Dupke, S.; Corre, J.P.; Lander, A.; Franz, T.; Monot, M.; Couture-Tosi, E.; Jouvion, G.; Leendertz, F.H.; et al. Capsules, Toxins and AtxA as Virulence Factors of Emerging Bacillus cereus Biovar anthracis. PLoS Negl. Trop. Dis. 2015, 9, e0003455. [Google Scholar] [CrossRef]
  43. Sue, D.; Hoffmaster, A.R.; Popovic, T.; Wilkins, P.P. Capsule Production in Bacillus cereus Strains Associated with Severe Pneumonia. J. Clin. Microbiol. 2006, 44, 3426–3428. [Google Scholar] [CrossRef] [PubMed]
  44. Oh, S.Y.; Budzik, J.M.; Garufi, G.; Schneewind, O. Two Capsular Polysaccharides Enable Bacillus cereus G9241 to Cause Anthrax-like Disease. Mol. Microbiol. 2011, 80, 455–470. [Google Scholar] [CrossRef] [PubMed]
  45. Sharma, S.; Bhatnagar, R.; Gaur, D. Bacillus anthracis Poly-γ-d-Glutamate Capsule Inhibits Opsonic Phagocytosis by Impeding Complement Activation. Front. Immunol. 2020, 11, 462. [Google Scholar] [CrossRef]
  46. Owen, M.P.; Schauwers, W.; Hugh-Jones, M.E.; Kiernan, J.A.; Turnbull, P.C.; Beyer, W. A Simple, Reliable M’Fadyean Stain for Visualizing the Bacillus anthracis Capsule. J. Microbiol. Methods 2013, 92, 264–269. [Google Scholar] [CrossRef] [PubMed]
  47. Luna, V.A.; Nguyen, K.; Gilling, D.H. Identification of pBC218/pBC210 Genes of Bacillus cereus G9241 in Five Florida Soils Using qPCR. Int. Sch. Res. Not. 2014, 2014, 197234. [Google Scholar] [CrossRef] [PubMed]
  48. Marston, C.K.; Ibrahim, H.; Lee, P.; Churchwell, G.; Gumke, M.; Stanek, D.; Gee, J.E.; Boyer, A.E.; Gallegos-Candela, M.; Barr, J.R.; et al. Anthrax Toxin-Expressing Bacillus cereus Isolated from an Anthrax-like Eschar. PLoS ONE 2016, 11, e0156987. [Google Scholar] [CrossRef]
  49. Pena-Gonzalez, A.; Rodriguez-R, L.M.; Marston, C.K.; Gee, J.E.; Gulvik, C.A.; Kolton, C.B.; Saile, E.; Frace, M.; Hoffmaster, A.R.; Konstantinidis, K.T. Genomic Characterization and Copy Number Variation of Bacillus anthracis Plasmids pXO1 and pXO2 in a Historical Collection of 412 Strains. MSystems 2018, 3, e00065-18. [Google Scholar] [CrossRef]
  50. Scorpio, A.; Chabot, D.J.; Day, W.A.; O’brien, D.K.; Vietri, N.J.; Itoh, Y.; Mohamadzadeh, M.; Friedlander, A.M. Poly-γ-glutamate Capsule-Degrading Enzyme Treatment Enhances Phagocytosis and Killing of Encapsulated Bacillus anthracis. Antimicrob. Agents Chemother. 2007, 51, 215–222. [Google Scholar] [CrossRef]
  51. Scorpio, A.; Chabot, D.J.; Day, W.A.; Hoover, T.A.; Friedlander, A.M. Capsule Depolymerase Overexpression Reduces Bacillus anthracis Virulence. Microbiology 2010, 156, 1459–1467. [Google Scholar] [CrossRef]
  52. Wilson, M.K.; Vergis, J.M.; Alem, F.; Palmer, J.R.; Keane-Myers, A.M.; Brahmbhatt, T.N.; Ventura, C.L.; O’Brien, A.D. Bacillus cereus G9241 Makes Anthrax Toxin and Capsule like Highly Virulent B. anthracis Ames but Behaves Like Attenuated Toxigenic Nonencapsulated B. anthracis Sterne in Rabbits and Mice. Infect. Immun. 2011, 79, 3012–3019. [Google Scholar] [CrossRef]
Pathogens 12 00693 g001 550
Figure 1. The parasporal bodies of B. thuringiensis ATCC 10792 (the arrows) denoted by the presence of dark blue and diamond-shape objects (a). No parasporal bodies were present on the specimen of Bacillus spp. JMT (b).
Figure 1. The parasporal bodies of B. thuringiensis ATCC 10792 (the arrows) denoted by the presence of dark blue and diamond-shape objects (a). No parasporal bodies were present on the specimen of Bacillus spp. JMT (b).
Pathogens 12 00693 g001
Pathogens 12 00693 g002 550
Figure 2. Microscopic images were captured (at 100× oil immersion) of Bacillus spp. JMT was observed with India ink staining. The unstainable and thick halo surrounding the bacilli on the images was designated the capsule of Bacillus spp. JMT (the black arrows) (a). Metachromatic color (magenta) and well-demarcated around the bacilli in Liu’s stain (the black arrows) (b).
Figure 2. Microscopic images were captured (at 100× oil immersion) of Bacillus spp. JMT was observed with India ink staining. The unstainable and thick halo surrounding the bacilli on the images was designated the capsule of Bacillus spp. JMT (the black arrows) (a). Metachromatic color (magenta) and well-demarcated around the bacilli in Liu’s stain (the black arrows) (b).
Pathogens 12 00693 g002
Pathogens 12 00693 g003 550
Figure 3. The edema of the diseased Chinese soft-shell turtle can be observed in the neck (dorsal side) (a) and limbs (ventral side) (the white arrows) (b). The enlarged spleen (the black arrows) was well recognized in the abdominal cavity (b).
Figure 3. The edema of the diseased Chinese soft-shell turtle can be observed in the neck (dorsal side) (a) and limbs (ventral side) (the white arrows) (b). The enlarged spleen (the black arrows) was well recognized in the abdominal cavity (b).
Pathogens 12 00693 g003
Pathogens 12 00693 g004 550
Figure 4. A number of inflammatory cells (mainly heterophil) infiltrating the spleen (the arrows) (a). Bacilli denoted by the blue color stained with Gram’s staining in the lung (the arrows) (b). Edema (* asterisk) showed (c,d).
Figure 4. A number of inflammatory cells (mainly heterophil) infiltrating the spleen (the arrows) (a). Bacilli denoted by the blue color stained with Gram’s staining in the lung (the arrows) (b). Edema (* asterisk) showed (c,d).
Pathogens 12 00693 g004
Table
Table 1. Background information of Bacillus spp. JMT isolates were analyzed in this study.
Table 1. Background information of Bacillus spp. JMT isolates were analyzed in this study.
TimeCase (Isolate) ID.LocationTurtle AgeCumulative
Mortality		Number of SamplesFixed
Tissues a
April 2016JMT105-1Changzhi Township (PT)adult-5−
May 2016JMT105-2Ligang Township (PT)subadult-8−
March 2016JMT105-3Ligang Township (PT)adult-4−
April 2016JMT105-4Jiuru Township (PT)adult-3−
April 2016JMT105-5Meinong District (K)adult-7−
April 2016JMT105-6Jiuru Township (PT)subadult-3−
May 2016JMT105-7Yanpu Township (PT)subadult-5−
May 2017JMT106-8Jiuru Township (PT)subadult800/1000 (80%)3−
March 2017JMT106-9Changzhi Township (PT)adult300/1200 (25%)3−
March 2017JMT106-10Ligang Township (PT)subadult450/1000 (45%)2−
June 2017JMT106-11Changzhi Township (PT)adult900/1300 (69.2%)4−
March 2017JMT106-12Jiuru Township (PT)adult450/1350 (33.3)3−
April 2017JMT106-13Ligang Township (PT)subadult150/700 (21.4)4−
May 2017JMT106-14Ligang Township (PT)subadult750/1400 (53%)4−
April 2018JMT107-15Jiuru Township (PT)adult650/900 (72.2%)3−
June 2018JMT107-16Jiuru Township (PT)subadult380/1100 (34.5%)3−
May 2018JMT107-17Yanpu Township (PT)subadult230/900 (25.6%)5+
March 2018JMT107-18Meinong District (K)subadult300/1500 (20%)7+
March 2018JMT107-19Ligang Township (PT)adult400/1300 (30.8%)3+
April 2018JMT107-20Ligang Township (PT)subadult600/1300 (46.2%)3+
May 2018JMT107-21Ligang Township (PT)adult1200/1200 (100%)2+
March 2018JMT107-22Changzhi Township (PT)adult800/1450 (55.2%)5+
May 2019JMT108-23Changzhi Township (PT)subadult130/1200 (59.1%)4+
May 2019JMT108-24Ligang Township (PT)adult350/1200 (29.2%)3+
May 2019JMT108-25Ligang Township (PT)subadult400/1250 (32%)3+
June 2019JMT108-26Yanpu Township (PT)subadult350/1450 (24.1%)6+
March 2019JMT108-27Pingtung City (PT)subadult500/130 (38.5%)5+
March 2019JMT108-28Ligang Township (PT)adult400/900 (44.4%)5+
April 2020JMT109-29Yanpu Township (PT)subadult700/1250 (56%)4+
May 2020JMT109-30Ligang Township (PT)subadult130/2400 (54.2%)4+
May 2020JMT109-31Ligang Township (PT)subadult2300/2300 (100%)3+
June 2020JMT109-32Jiuru Township (PT)adult900/1300 (69.2%)3+
(K) denotes Kaohsiung City. (PT) denotes Pingtung County. (-) denotes no data in the column of Cumulative Mortality. (−) denote no tissues available in the column of fixed tissue. (+) denotes there were tissues for histopathology analysis. a included the heart, lungs, liver, spleen, kidneys, and GI-tract.
Table
Table 2. Primer sets used in the present study.
Table 2. Primer sets used in the present study.
SpeciesPrimerPrimer Sequence (5′-3′)Amplicon Size (bp)Target GeneReference
Bacteria27fAGAGTTTGATCMTGGCTCAG~150016S rRNA[28]
1492rTACGGYTACCTTGTTACGACTT
B. cereus groupBCGSH-1FGTGCGAACCCAATGGGTCTTC400groEL[29]
BCGSH-1RCCTTGTTGTACCACTTGCTC
B. anthracisBASH-2FGGTAGATTAGCAGATTGCTCTTCAAAAGA253gyrB[29]
BASH-2RACGAGCTTTCTCAATATCAAAATCTCCGC
B. thuringiensisBTJH-1FGCTTACCAGGGAAATTGGCAG299gyrB[29]
BTJH-RATCAACGTCGGCGTCGG
B. cereusBCJH-FTCATGAAGAGCCTGTGTACG475gyrB[29]
BCJH-1RCGACGTGTCAATTCACGCGC
B. mycoidesBMSH-FTTTTAAGACTGCTCTAACACGTGTAAT603gyrB[29]
BMSH-RTTCAATAGCAAAATCCCCACCAAT
B. thuringiensisXRE1AAGATATTGCAAGCGGTAAGAT246XRE[30]
XRE2GTTTTGTTTCAGCATTCCAGTAA
B. anthracislef3CTTTTGCATATTATATCGAGC385lef[31]
lef4GAATCACGAATATCAATTTGTAGC
B. cereus G9241JMlef436fGCGGTCATGGTGATGTAGGT436lefThis study
JMlef436rACGTTCAGTGCCTTTTCAGT
JMpagA318fCTGGGACGGCTCCAATCTAC318pagAThis study
JMpagA318rAGCCTGTATCCACCCTCACT
JMcya882fTGCACCTGACCATAGAACGG882cyaThis study
JMcya882rTCCGGTTTCCTCTCAATTCCA
JMhasA700fGTAATGGCCACAACTGGTCACGTGAAC700hasAThis study
JMhasA700rAACGTGTTCCCCAACCATTAGA
JMbpsB208fAAAACAGACCCCCAATCACC208bpsBThis study
JMbpsB208rGTTTCCGCAGGTCAGCATTT
Table
Table 3. Laboratory findings of the characteristics of B. cereus group in this study.
Table 3. Laboratory findings of the characteristics of B. cereus group in this study.
Bacterial IsolatesColony
Morphology		Microscopical ExaminationPCR Assay
Rhizoid
Growth		Parasporal
Body		Capsule Expression B. c aB. t aXRE bpsBhasApagAcyalef
Bacillus spp. JMT (32)−−32+−32+−32+32+32+32+−
B. cereus ATCC 14579 (1)−−−1+−−−−−−−
B. thuringiensis ATCC 10792 (1)−1+−−1+1+−−−−−
Colony morphology: (+) denotes a positive of rhizoid growth visualization; (−) denote a negative of rhizoid growth visualization. Microscopical examination: (+) denotes a positive of parasporal body and capsule visualization; (−) denote a negative of parasporal body and capsule visualization. PCR assay: (+) and (−) respectively denote a positive and negative amplicon after PCR with different primer sets. B. c: B. cereus; B. t: B. thuringiensis. a The results were based on multiplex PCR assay in this study.
Table
Table 4. The calculation of OrthoANI and dDDH between paired genomes.
Table 4. The calculation of OrthoANI and dDDH between paired genomes.
12345678
1 96.64 #96.56 #96.83 #92.0691.8094.8189.68
270.20 # 99.98 #97.63 #91.9791.7794.8289.65
370.10 #99.80 # 97.59 #91.9191.6194.8389.66
472.20 #78.40 #78.30 # 92.1491.9295.15 #89.58
546.3046.1046.1046.60 96.80 #91.7189.49
645.4045.2045.0045.4071.20 # 91.4589.58
759.0058.7058.6060.9045.0044.10 89.54
838.7038.6038.6038.5038.1038.3038.20
# denotes the value is higher than the OrthoANI or GGD thresholds. dDDH values are in the bottom-left and OrthoANI values in the top-right. Value unit: %. Strains: 1, Bacillus tropicus JMT105-2(CP115303); 2, Bacillus cereus G9241(NZ_CP009590); 3, Bacillus cereus strain 03BB87(NZ_CP009941); 4, Bacillus tropicus strain EMB20(NZ_CP078081); 5, Bacillus cereus ATCC 14579(NC_004722); 6, Bacillus thuringiensis serovar berliner ATCC 10792(NZ_CM000753); 7, Bacillus anthracis str. Ames (NC_003997); 8, Bacillus mycoides DSM 2048 (NZ_CM000742).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tsai, J.-M.; Kuo, H.-W.; Cheng, W. Retrospective Screening of Anthrax-like Disease Induced by Bacillus tropicus str. JMT from Chinese Soft-Shell Turtles in Taiwan. Pathogens 2023, 12, 693. https://doi.org/10.3390/pathogens12050693

AMA Style

Tsai J-M, Kuo H-W, Cheng W. Retrospective Screening of Anthrax-like Disease Induced by Bacillus tropicus str. JMT from Chinese Soft-Shell Turtles in Taiwan. Pathogens. 2023; 12(5):693. https://doi.org/10.3390/pathogens12050693

Chicago/Turabian Style

Tsai, Jia-Ming, Hsin-Wei Kuo, and Winton Cheng. 2023. "Retrospective Screening of Anthrax-like Disease Induced by Bacillus tropicus str. JMT from Chinese Soft-Shell Turtles in Taiwan" Pathogens 12, no. 5: 693. https://doi.org/10.3390/pathogens12050693

APA Style

Tsai, J. -M., Kuo, H. -W., & Cheng, W. (2023). Retrospective Screening of Anthrax-like Disease Induced by Bacillus tropicus str. JMT from Chinese Soft-Shell Turtles in Taiwan. Pathogens, 12(5), 693. https://doi.org/10.3390/pathogens12050693

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