Characterization of a Deep Sea Bacillus toyonensis Isolate: Genomic and Pathogenic Features

Introduction

Deep hydrothermal vents are localized areas of the deep sea with high tectonic activities that cause the escaping of hot, mineral-rich water from the sea floor, which forms chimney-like black or white smokers. Compared to other areas of the deep ocean, hydrothermal regions are rich in various sulfur compounds and have large amounts of small molecules of hydrocarbons, which provide the basis for the development of chemoautotroph microbial communities that can utilize the sulfur compounds and hydrocarbons (Marteinsson et al., 1996; Ruger et al., 2000; Burgaud et al., 2009; Le Calvez et al., 2009; Wang et al., 2020). The food and energy produced by the microorganisms enable the growth of invertebrate animals, such as mussels, tube worms, and crabs, thus forming the deep-sea food web (Dover, 2002; Nishijima et al., 2010; Van Dover, 2014).

Bacillus spp. are present in soil, air, and the deep-sea (Liu et al., 2006; Stenfors Arnesen et al., 2008). In animals, B. cereus can be a transient presence in the guts of mammals, and for some insects, B. cereus is part of the gut microbial communities (Jensen et al., 2003). They have been studied mostly as probiotics for animal feeds and medicines (Cutting, 2011; Elshaghabee et al., 2017). However, in the B. cereus group, which is a species complex within the Bacillus genus, some members, such as Bacillus anthracis, Bacillus thuringiensis, and B. cereus, are well known as pathogens and can cause fatal diseases in humans and animals (Bottone, 2010). Bacillus anthracis causes anthrax mediated by two key virulence factors: a tripartite toxin and a poly-D-glutamic acid capsule (Bourgogne et al., 2003; Liu et al., 2014). Bacillus thuringiensis, which has been used in agriculture as a biopesticide, produces proteins that aggregate into crystals and are toxic to insects (Berry et al., 2002; Bravo et al., 2007; Li et al., 2008). Some members of the Bacillus cereus group can produce several toxins, including the single-protein toxin, cytotoxin K (CytK), and the two protein toxin complexes formed by hemolysin BL (Hbl) and the non-hemolytic enterotoxin (Nhe) (Stenfors Arnesen et al., 2008; Ehling-Schulz et al., 2019).

Bacillus toyonensis is a member of the Bacillus cereus group, with strain BCT-7112^T being the type strain (Jimenez et al., 2013). Strain BCT-7112^T is a soil isolate and has been used for many years as an additive of animal feeds. In a recent study, it was demonstrated that although BCT-7112^T possesses the genes coding for hemolysin BL (Hbl) and the non-hemolytic enterotoxin (Nhe), the expression levels of these genes were undetectable or very low compared to the B. cereus reference strains, suggesting that the ability of BCT-7112^T to produce functional enterotoxins is low or unlikely (Abdulmawjood et al., 2019). Recently, two other members of B. toyonensis were reported as potentially novel agents for biocontrol of plants (Lopes et al., 2017; Rojas-Solis et al., 2020). B. toyonensis-like strains have been shown to exist in neritic regions (Okaiyeto et al., 2015; Tallur et al., 2016), however, to our knowledge, no B. toyonensis isolates from the deep sea (depth >1,000 m) have been documented.

In this study, we reported the isolation of a B. toyonensis strain, P18, from a deep sea hydrothermal field and the analysis of its genomic, biological, and virulence features. We found that P18 possesses unique genetic traits and was able to cause lethal infection in fish and mice upon artificial inoculation. Our results provide the first evidence that pathogenic B. toyonensis exits in the deep-sea environment.

Materials and Methods

Ethics Statement

The experiments involving live animals were approved by the Animal Ethics Committee of Institute of Oceanology, Chinese Academy of Sciences. All methods were carried out in accordance with the relevant guidelines.

Animals

The mice and fish used in this study were purchased and maintained as reported previously (Gu et al., 2019). Briefly, clinically healthy BALB/c mice (female, 8–10 weeks, and 14 ± 2 g) were purchased from Qingdao Daren Fortune Animal Technology Co., Ltd. Clinically healthy Japanese flounder (Paralichthys olivaceus) (15 ± 0.5g) were purchased from a local fish farm in Shandong Province, China and maintained at 20°C in aerated seawater in the laboratory. For tissue collection, the fish and mice were euthanized with an overdose of tricaine methanesulfonate (Sigma, St. Louis, USA) and ketamine (80 mg/kg) (Ketavet, Pfizer, Berlin, Germany), respectively (Gu et al., 2019).

Cell Culture

The cell lines used in this study were cultured as reported previously (Gu et al., 2019; Zhao et al., 2019). Briefly, RAW264.7 cells and Hela cells were cultured in DMEM (Gibco, USA) with 10% fetal bovine serum (Gibco, USA) and 1% penicillin and streptomycin (Beyotime Biotechnology, China) at 37°C and 5% CO2. FG-9307 cells were cultured in L-15 (Genom, China) with 10% fetal bovine serum and 1% penicillin and streptomycin at 24°C.

Bacterial Isolation

P18 was isolated from a sediment sample collected from the hydrothermal field in southern Okinawa Trough (127°04′30.980″ E, 27°16′01.064″ N, depth of 1446 m), where the water temperature was around 4°C. The sediment sample was obtained by push-cores on the Remotely Operated Vehicle (ROV) equipped on the KEXUE vessel (Sun et al., 2015). The sample was plated on Marine agar 2216 medium (Sun et al., 2015) and incubated at 28°C under aerobic conditions for 3 days. One of the colonies that appeared on the plate was named P18, which, after re-purification, was stored at −80°C with 30% (v/v) glycerol.

Morphology and Growth of P18

For morphological observation, P18 was cultured in Marine 2216E medium at 28°C to an OD₆₀₀ of 0.8. The culture was centrifuged at 5,000 rpm for 10 min, and the bacterial pellet was re-suspended in ddH₂O. For spore observation, P18 was cultured at 28°C in Marine 2216E medium containing 5mg/L MnSO₄ for 3 days with shaking (180 rpm/min) (Osadchaia et al., 1997); the culture was then centrifuged and re-suspended in ddH₂O. After fixation with glutaraldehyde and dehydration with acetone, the vegetative cells and spores of P18 were observed with a transmission electron microscope (Hitachi, JEM-2100, Japan). The spores were also stained with a Spore Staining Kit (HaiBo, China) and observed with a light microscope (Ti-S/L100, Nikon, Japan). The temperature range of growth was determined by culturing P18 at 4, 16, 28, 37, and 50°C in Marine 2216E medium for 36 h to 7 days. The pH range was determined by culturing P18 in Marine 2216E medium of different pH (pH 5–10, with an increment of 1 pH unit). The NaCl range was determined by culturing P18 in 2216E medium with different salinities (NaCl 0–8%, w/v, with an increment of 1%). The motility of P18 was determined as reported previously (Gu et al., 2019).

Genome Analysis

Genome sequencing was carried out by Novogene (Beijing, China) as reported previously (Zhao et al., 2019). Briefly, clean data were obtained by controlling the quality of the original data and filtering out the low-quality sequences for subsequent analysis. SMRT software (Link V5.0.1) (Berlin et al., 2015) was used to assemble reads. Plasmid sequences were distinguished from chromosomal sequence based on the method of Gomez-Escribano et al. (2015). Coding genes were predicted with GeneMarkS 4.28 software (http://topaz.gatech.edu/) (Besemer et al., 2001). Genomic islands were predicted using IslandPath-DIOMB 0.2 software (Hsiao et al., 2003). ncRNA sequences were predicted using rRNAmmer 1.2 (Lagesen et al., 2007) and tRNAscan 1.3.1 (Schattner et al., 2005) softwares and the Rfam (Griffiths-Jones et al., 2005) database. Repeated sequences and prophages were predicted using RepeatMasker 4.0.5 (Chen, 2004) and phiSpy 2.3 (Akhter et al., 2012) softwares, respectively. All programs were run with default settings. Functional annotation was performed using the following databases: KEGG (http://www.genome.jp/kegg/) (Kanehisa et al., 2004), COG (http://www.ncbi.nlm.nih.gov/COG/) (Tatusov et al., 2003), SwissProt (http://www.uniprot.org/) (Bairoch and Apweiler, 2000), GO (http://www.geneontology.org/) (Ashburner et al., 2000), and NR (http://www.ncbi.nlm.nih.gov/RefSeq/) (Li et al., 2002). Diamond (v0.7.9.58) blastp was used for the analysis with the COG, KEGG, and NR databases (identity >40% and E-value <1E-05). Blast was used for the analysis with the SwissProt database (identity >40% and E-value <1E-05). Gene mapping was performed using Circos software by combining the assembled genome sequence with the predicted results of the coding genes. Virulence related genes were predicted using the VFDB (Virulence Factors of Pathogenic Bacteria Database; http://www.mgc.ac.cn/VFs/) database (Chen et al., 2005) with Diamond (v0.7.9.58) blastp (identity >40% and E-value <1E-05). The genomes of P18 and BCT-7112^T (Harris, 2007) were compared using MUMmer 3.23 software (Delcher et al., 2003) to determine the collinearity. LASTZ 1.02.00 software was used to confirm Translocation/Trans, Inversion/Inv, and Trans+Inv regions between P18 and BCT-7112^T. Gene clusters were identified using the cd-hit software (Li and Godzik, 2006), with identity >50% and length difference <30%. The average nucleotide identity (ANI) between P18 and other strains of Bacillus (Supplementary Table S1) was calculated using the EzBiocloud web service (https://www.ezbiocloud.net/tools/ani) (Kim et al., 2014). The chromosome and plasmid sequences of P18 have been deposited in GenBank under the accession numbers CP064875 and CP064876, respectively.

Infection Assay

In vivo infection was performed as reported previously (Gu et al., 2019). Briefly, for infection of flounder, the fish were injected intramuscularly with P18 at the dose of 3 × 10⁴ to 3 × 10⁶ CFU/g (with 10-fold interval) and monitored for mortality for 10 days. For the comparison of the mortality-inducing capacity of P18 and Bacillus subtilis 168 [purchased from China General Microbiological Culture Collection Center (CGMCC, http://www.cgmcc.net)], the bacteria were cultured to an OD₆₀₀ of 0.8. The bacteria were collected by centrifugation and resuspended in PBS. The same dose of P18 and B. subtilis 168 were used to inject flounder as above. For infection of mice, BALB/c mice were injected intraperitoneally with P18 at the dose of 6.25 × 10⁵ to 5 × 10⁶ CFU/g (with 2-fold interval) and monitored for mortality for 10 days. For all infections, the control groups were injected with PBS. The median lethal dose (LD50) was determined as reported previously (Gu et al., 2019). For the detection of bacterial dissemination in host tissues, liver, spleen, kidney, and blood of flounder (five fish/time point) were sampled at 12, 24, and 48 hpi; liver, spleen, kidney, and intestine of mice (three mice/time point) were sampled at 2, 4, and 6 hpi. The tissues were homogenized in PBS, and bacterial numbers in the homogenates were determined by plate count as reported previously (Zhao et al., 2019).

Histopathology

For histopathological analysis, the tissues prepared above were immersed in 4% PFA and stored at room temperature. Tissue sectioning and HE staining were performed by Wuhan Servicebio Technology Co., Ltd. (Wuhan, China). After dewaxing, hematoxylin staining, eosin staining, and dehydration, the tissue sections were observed with an optical microscope (Eclipse E100, Nikon, Japan).

Cellular Infection

P18 was cultured in L-15 or DMEM medium at 28°C overnight. The culture was centrifuged at 8,000 rpm, and the bacteria and the supernatant were collected. The bacteria were resuspended in L-15 or DMEM. The supernatant was filtered with a 0.22 Millipore filters (Millex^®-GP, Cork, Ireland). FG-9307, RAW264.7, and Hela cells were culture as described above for overnight in 24-well plates. The cells were infected with P18 at a MOI of 3:1 or treated with the bacterial supernatant (5%, v/v). After incubation at 28°C for 0.5–1.5 h, the cells were observed with a light microscope.

Serum Resistance and Hemolytic Activity

Serum resistance and hemolytic activity were analyzed as reported previously (Gu et al., 2019). Briefly, to examine the hemolytic activity of P18, P18 and B. subtilis 168 were cultured to an OD₆₀₀ of 0.8. The bacteria were collected by centrifugation and resuspended in PBS. Ten microliters of the bacterial suspension were added to the filter paper on a LB agar plate containing 5% sterile sheep blood (Hopebio, China). Triton X-100 and PBS were used as the positive and negative controls, respectively. The plate was incubated at 28°C for 24 h and then observed for hemolysis.

Results

Morphological and Growth Characteristics of P18

Strain P18 was isolated from a deep-sea sediment sample collected from the hydrothermal field in southern Okinawa Trough (127°04′30.980″ E, 27°16′01.064″ N, depth 1446 m). Transmission electron microscopy (TEM) showed that P18 was rod-shaped and exhibited peritrichous flagella (Figure 1A). It is able to produce endospores under specific conditions (Figure 1B). P18 exhibited swimming and swarming abilities in 0.3 and 0.5% agar, respectively (Figure 1C). P18 grew well at pH 5–9, 4–37°C, and in the presence of up to 8.0% NaCl (Supplementary Figure S1).

FIGURE 1

Figure 1 Morphology and motility of P18. (A) P18 cells of spore (upper) and vegetative (lower) growth observed with a transmission electron microscope. (B) P18 cells in vegetative (pinkish red) and spore (bluish green) forms observed with a microscope. (C) Growth of P18 in marine 2216E containing 0.3, 0.5, or 1.5% (w/v) agar.

Genome Features of P18

P18 contains a circular chromosome of 5 250 895 bp and a plasmid (named pBtP18) of 536 892 bp, with an average G+C content of 35.42 and 32.86%, respectively (Figure 2; Table 1). The chromosome and the plasmid encode 5,380 and 523 genes, respectively, with average gene lengths of 820 and 754 bp, respectively. There are 2,229 genes of P18 that encode hypothetical proteins and could not be annotated functionally based on the COG database. P18 has 14 sets of rRNAs (5S, 16S, and 23S rRNA), 105 tRNAs, and 4 sRNAs, all of which are present on the chromosome, and 20 genomic islands (GIs), 16 of which are present on the chromosome (Table 1). In addition, P18 possesses 260 interspersed repeats, 367 tandem repeats, 298 mini-satellites, and 1 microsatellite. P18 exhibits 98.67, 91.04, 91.22, 91.43, and 91.54% ANI values with Bacillus toyonensis BCT-7112^T, Bacillus luti TD41^T, Bacillus wiedmannii FSL W8-0169^T, Bacillus cereus ATCC 14579^T, and Bacillus thringiensis ATCC 10792^T, respectively (Supplementary Table S1). Based on these results, P18 was classified as a member of B. toyonensis.

FIGURE 2

Figure 2 Genomic maps of the P18 chromosome (A) and plasmid (B). (A) From the outside circle to the inside circle are the coordinates of the genome sequence, annotation information of gene function (based on COG, KEGG, and GO databases), ncRNA, GC content of the genome, distribution of the value of GC skew. (B) From the outside circle to the inside circle are COG annotation information, the coordinates of the plasmid sequence, GC content of the plasmid, and distribution of the value of GC skew.

TABLE 1

Table 1 Features of P18 and BCT-7112^T.

Comparative Analysis of the Genomes of P18 and Its Close Homologue BCT-7112^T

P18 is similar to B. toyonensis BCT-7112^T in gene number, average gene length, G+C content, and the number of non-coding RNAs, but has a larger genome size and possesses one, rather than two, plasmids (Table 1). However, the size of pBtP18 is much larger than that of either of the plasmids of BCT-7112 ^T, and is 6.3 times of the combined sizes of the two plasmids of BCT-7112^T. pBtP18 also possesses more GIs than the BCT-7112^T plasmids, which makes the total number of GIs in P18 higher than that in BCT-7112^T. Most of the chromosomal genes of P18 and BCT-7112^T exhibit strong collinearity, but the genes in the plasmids of the two strains exhibit translocation and inversion (Figure 3A). P18 and BCT-7112^T share 4151 common gene clusters, and contain 1401 and 529 specific genes, respectively. Of the 1,401 P18-unique genes, 960 encode hypothetical proteins with unknown function in the COG database, and 441 were classified into 20 COG functional categories, with the majority of these genes belonging to the categories of R (General function), K (Transcription), G (Carbohydrate transport and metabolism), V (Defense mechanisms), M (Cell wall/membrane/envelope biogenesis), T (Signal transduction mechanisms), and X (Mobilome: prophages, transposons) (Figure 3B).

FIGURE 3

Figure 3 Genome comparison between P18 and BCT-7112^T. (A) Chromosome (left) and plasmid (right) alignment of P18 and BCT-7112^T. Upper and lower axes of the linear synteny graph are constructed with the same proportion of size reduction of both sequences. Each pair nucleic acid sequence of the two alignments is marked in the coordinate diagram according to its position information, and the height of the filled color in the block indicates similarity of sequence alignment. The color of the lines between the two axes indicates the type of comparison. (B) COG functional classification of the P18-specific genes. (C) The functional categories of P18-specific virulence genes. The genes of each category are listed on the right.

Virulence-Related Genes of P18

A total of 319 virulence genes were identified in P18 based on a database of known virulence genes (Supplementary File S1). Most of these virulence genes are also present in BCT-7112^T. Twenty-four of the virulence genes are specific to P18 and absent in BCT-7112^T. These genes are classified into the functions associated with capsule, regulator, flagella, adhesion, perforin, and lipase (Figure 3C). P18 carries a number of toxin-related genes, including hemolysins (hlyIII, hblA, hblB, hblC, and hblD), anthrolysin (alo), phospholipase C (SM-PLC, PI-PLC, and PC-PLC), and non-hemolytic enterotoxins (nheA, nheB, and nheC).

Infectivity and Lethality of P18 to Fish and Mice

Live infection study in fish and mice models showed that the LD₅₀ of P18 in Japanese flounder and mice were 1.4 × 10⁵ and 2.5 × 10⁶ CFU/g, respectively, when P18 was inoculated into the animals via i.m. (for fish) and i.p. (for mice) injection. When inoculated into flounder at the dose of 3 × 10⁵ CFU/g, P18 induced 80% accumulated mortality, whereas B. subtilis 168 in no mortality (Supplementary Figure S2). P18-injected fish began to exhibit skin ulceration and swelling at 12 hpi, and the symptom worsened over time (Supplementary Figure S3A). P18 was detected in the blood, spleen, and kidney of the infected fish from 12 to 48 hpi (Figure 4A). Mice infected with 2.5 × 10⁶ CFU/g of P18 exhibited weakness in the limb and vomiting at 4 hpi, and mortality began to occur at 6 hpi (Supplementary Figure S3B). At 6 hpi, bacterial dissemination in the liver, spleen, intestine, and kidney of the mice was observed (Figure 4B).

FIGURE 4

Figure 4 Dissemination of P18 in infected animals. Japanese flounder (A) and mice (B) were inoculated with P18, and bacterial recovery from the tissues was determined at different time points. The results are the means of triplicate experiments and shown as means ± Standard Error of Mean (SEM).

Tissue Lesion Induced by P18

Histopathology showed that in P18-infecetd fish, the arrangement of kidney renal tubules became looser, and the lymphocyte number in the mesenchyme dramatically decreased (Figure 5Aa). In spleen, there was a considerable reduction in the amount of immune cells in the white pulp and mesenchyme, accompanied by slight expansion of splenic sinus (Figure 5Ab). In P18-infecetd mice, kidney exhibited narrowed renal capsules and hemorrhage, and the intestine exhibited abscission of epithelial cells and hyperemia of blood vessels (Figure 5B).

FIGURE 5

Figure 5 Histopathologic examination of P18-infected fish (A) and mice (B). Japanese flounder (A) and mice (B) were infected with or without P18 for 48 h (A) or 6 h (B), and the tissue structures were examined with a microscope. Bar size, 50 (A) and 100 (B) µm.

Cytopathic Effect of P18

As shown in Figure 6, when RAW264.7 cells (mouse monocytes-macrophages) were incubated with P18 at an MOI of 3:1, the cells rapidly rounded out, and cell death was observed at 1.5 h after incubation. When incubated with HeLa cells (human epithelial cells) or the fish epithelial cells FG-9307, P18 induced rapid cellular membrane bubbling and damage that eventually led to cell death. Similar cytopathic effects were observed with the cells treated with the culture supernatant of P18 (Figure 6).

FIGURE 6

Figure 6 The effect of P18 on mammalian and fish cells. RAW264.7 cells (upper), Hela cells (middle), and FG-9307 cells were treated with or without (control) P18 or its culture supernatant, and the cells were observed with a microscope.

Serum Resistance and Hemolytic Activity of P18

Serum resistance analysis showed that P18 exhibited survival rates of 78 and 70% in flounder and mice sera, respectively. Hemolytic analysis indicated that P18 caused lysis of sheep blood cells to the effect comparable to that caused by 5% Triton X-100, but, unlike Triton X-100, P18 induced the formation of double hemolysis halos on the blood agar plate (Figure 7).

FIGURE 7

Figure 7 Hemolytic activity of P18. P18, Bacillus subtilis 168 (a non-hemolytic bacterium), 5% Triton X-100, 0.2% Triton X-100, and PBS were spotted onto the filter discs in sheep blood agar plate, and hemolysis was observed after 24 h incubation.

Discussion

In this study, we examined the biological features of a B. toyonensis strain, P18, isolated from a deep hydrothermal field. We found that P18 was able to form spores and grow at conditions resembling that of the local deep sea environment (low temperature and high salinity). However, we do not know whether P18 exists predominately as vegetative cells or spores or both under the native condition. P18 was classified as a member of B. toyonensis based on its high ANI with the type strain of B. toyonensis, BCT-7112^T. BCT-7112^T is a soil bacterium isolated in Japan in 1966 for use as a feed additive and was re-classified in 2013 as a novel member of the Bacillus cereus group (Jimenez et al., 2013). Compared to BCT-7112^T, the genome of P18 is 762368 bp larger and consequently possesses 664 more coding genes. One distinguishing feature of P18 is the presence of an exceedingly large plasmid, pBtP18, which in part accounts for the larger genomic size of P18. In contrast to the two plasmids of BCT-7112^T, which contain a single or no GI, pBtP18 contains 4 GIs with genes involved in, among others, transposition. Consistently, collinearity analysis showed that the genes of pBtP18 exhibit marked translocation and inversion, which is most likely caused by the transposases in the GIs of pBtP18. These observations are in agreement with a previous report which showed that the GIs in the genome of B. cereus probably play a role in horizontal gene transfer and, as a result, cause gene gain or loss  (Zhang and Zhang, 2003). Since P18 is a spore-forming bacterium, the possibility exists that it may have originated from a non-deep sea environment, from which the spores were carried to the deep sea. Nevertheless, the distinct genetic features of P18, together with its growth property, observed in our study suggest a living adaption to the deep sea environment.

Although a large amount of genes are shared by P18 and BCT-7112^T, there are still 1401 P18-unique genes that are absent in BCT-7112^T. In a previous study of the genomes of 66 deep-sea bacteria, it was found that extracellular peptidase and carbohydrate metabolism genes are widespread in most bacteria (Li et al., 2016). In line with this observation, we found that some of the P18-unique genes with predictable functions are related to carbohydrate transport and metabolism, cell membrane/wall/envelope synthesis, transcription, signal transduction, transposition, and other general functions. However, it is of note that most of the P18-unique genes encode hypothetical proteins that cannot be annotated in function based on the current COG database. It is possible that these hypothetical genes may encode deep sea-unique proteins with functions to facilitate bacterial survival under the specific conditions of the hydrothermal filed.

B. toyonensis has been studied primarily as a probiotic for preventing microbial diseases of crops or improving the immune response of animals (Santos et al., 2018). To date, little is known about the pathogenicity of B. toyonensis. Several studies have provided evidences that suggested a lack of or low cytotoxic potential for some B. toyonensis strains (Guinebretiere et al., 2010; Johler et al., 2018; Abdulmawjood et al., 2019; Carroll et al., 2020; Mendez Acevedo et al., 2020). Two other reports showed that B. toyonensis could exhibit larvacidal effect against beet armyworm (Eski et al., 2018) or induce HeLa cell membrane damage (Miller et al., 2018). In our study, a large account of putative virulence genes were identified in P18, including those encoding capsule, flagella, phospholipase C, and cytotoxins, which are known to facilitate the infection of B. cereus and other bacteria (Dons et al., 2004; van Asten et al., 2004; Ramarao and Lereclus, 2005). In line with this observation, P18 possesses visible flagella and exhibited apparent swimming and swarming capacities. Live animal infection showed that P18 was able to induce acute infection in fish and mice upon im/ip injection and cause high mortality. Following inoculation, P18 rapidly appeared in multiple internal tissues, and evident histological changes were detected in some organs of the infected animals. These results suggest that P18 possesses in vivo invasion capacity that enables the bacteria to disseminate in host tissues and cause pathological lesions.

Genes coding for several types of cytotoxins were found in P18, including hbl, nhe, alo, and plc. Hbl and Nhe are prominent cytotoxins of the B. cereus group. Hbl can activate the NLRP3 inflammatory response and induce lysis of animal cells such as RAW267.4 (Sastalla et al., 2013). Similarly, Nhe triggers osmotic lysis of the target cells following pore formation in the plasma membrane (Fagerlund et al., 2008). Anthrolysin O (ALO) was identified as a hemolysin and a cholesterol-dependent cytolysin produced by Bacillus anthracis and is involved in anthrax pathogenesis (Bourdeau et al., 2009; Tonello and Zornetta, 2012). Phospholipase C is also an inducer of membrane damage by its ability to hydrolyze membrane phospholipids (Goldfine et al., 1993; Goebel and Kuhn, 2000). In accordance with the abundance of cytolytic genes in its genome, P18 displayed strong cytopathic effects on different types of vertebrate cells. The observation of cytotoxicity conferred by the culture supernatant of P18 suggests extracellular toxin production. In addition to fish cells and murine cells, P18 exhibited marked ability to lyse sheep red blood cells. The hemolytic activity of P18 may account for the hemorrhage observed in some of the tissues of P18-infected animals.

In summary, we demonstrated in this study that pathogenic B. toyonensis exists in deep hydrothermal systems. We found that B. toyonensis P18 possesses a large number of virulence genes and genes of unknown function. Consistently, P18 is highly toxic to fish and mammalian cells and can invade into host tissues and cause mortality. These results provide the first insight into the pathogenicity of B. toyonensis associated with deep sea environments.

Data Availability Statement

The data presented in the study are deposited in the GenBank repository, accession number CP064875 and CP064876.

Ethics Statement

The animal study was reviewed and approved by Animal Ethics Committee of Institute of Oceanology, Chinese Academy of Sciences.

Author Contributions

J-cL, HL, and YZ performed the experiments and analyzed the data. JZ obtained the deep sea sample and isolated the bacteria. LS conceived and designed the experiments. J-cL and LS wrote the paper. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the grants from Qingdao National Laboratory for Marine Science and Technology (QNLM2016ORP0309), the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA22050402), the Collaborative Research Grant (KLMVI-OP-202002) of CAS Key Laboratory of Molecular Virology & Immunology, Institut Pasteur of Shanghai, Chinese Academy of Sciences, and the Taishan Scholar Program of Shandong Province.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2021.629116/full#supplementary-material

References

Abdulmawjood A., Herrmann J., Riede S., Jimenez G., Becker A., Breves G. (2019). Evaluation of enterotoxin gene expression and enterotoxin production capacity of the probiotic strain Bacillus toyonensis BCT-7112T. PLoS One 14 (4), e0214536. doi: 10.1371/journal.pone.0214536

PubMed Abstract | CrossRef Full Text | Google Scholar

Akhter S., Aziz R. K., Edwards R. A. (2012). PhiSpy: a novel algorithm for finding prophages in bacterial genomes that combines similarity- and composition-based strategies. Nucleic Acids Res. 40 (16), e126. doi: 10.1093/nar/gks406

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashburner M., Ball C. A., Blake J. A., Botstein D., Butler H., Cherry J. M., et al. (2000). Gene ontology: tool for the unification of biology. Gene Ontology Consortium Nat. Genet. 25 (1), 25–29. doi: 10.1038/75556

CrossRef Full Text | Google Scholar

Bairoch A., Apweiler R. (2000). The SWISS-PROT protein sequence database and its supplement TrEMBL in 2000. Nucleic Acids Res. 28 (1), 45–48. doi: 10.1093/nar/28.1.45

PubMed Abstract | CrossRef Full Text | Google Scholar

Berlin K., Koren S., Chin C. S., Drake J. P., Landolin J. M., Phillippy A. M. (2015). Assembling large genomes with single-molecule sequencing and locality-sensitive hashing. Nat. Biotechnol. 33 (6), 623–630. doi: 10.1038/nbt.3238

PubMed Abstract | CrossRef Full Text | Google Scholar

Berry C., O’Neil S., Ben-Dov E., Jones A. F., Murphy L., Quail M. A., et al. (2002). Complete sequence and organization of pBtoxis, the toxin-coding plasmid of Bacillus thuringiensis subsp. Israelensis Appl. Environ. Microbiol. 68 (10), 5082–5095. doi: 10.1128/aem.68.10.5082-5095.2002

CrossRef Full Text | Google Scholar

Besemer J., Lomsadze A., Borodovsky M. (2001). GeneMarkS: a self-training method for prediction of gene starts in microbial genomes. Implications for finding sequence motifs in regulatory regions. Nucleic Acids Res. 29 (12), 2607–2618. doi: 10.1093/nar/29.12.2607

PubMed Abstract | CrossRef Full Text | Google Scholar

Bottone E. J. (2010). Bacillus cereus, a volatile human pathogen. Clin. Microbiol. Rev. 23 (2), 382–398. doi: 10.1128/CMR.00073-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Bourdeau R. W., Malito E., Chenal A., Bishop B. L., Musch M. W., Villereal M. L., et al. (2009). Cellular functions and X-ray structure of anthrolysin O, a cholesterol-dependent cytolysin secreted by Bacillus anthracis. J. Biol. Chem. 284 (21), 14645–14656. doi: 10.1074/jbc.M807631200

PubMed Abstract | CrossRef Full Text | Google Scholar

Bourgogne A., Drysdale M., Hilsenbeck S. G., Peterson S. N., Koehler T. M. (2003). Global effects of virulence gene regulators in a Bacillus anthracis strain with both virulence plasmids. Infect. Immun. 71 (5), 2736–2743. doi: 10.1128/iai.71.5.2736-2743.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Bravo A., Gill S. S., Soberon M. (2007). Mode of action of Bacillus thuringiensis Cry and Cyt toxins and their potential for insect control. Toxicon 49 (4), 423–435. doi: 10.1016/j.toxicon.2006.11.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Burgaud G., Le Calvez T., Arzur D., Vandenkoornhuyse P., Barbier G. (2009). Diversity of culturable marine filamentous fungi from deep-sea hydrothermal vents. Environ. Microbiol. 11 (6), 1588–1600. doi: 10.1111/j.1462-2920.2009.01886.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Carroll L. M., Cheng R. A., Kovac J. (2020). No Assembly Required: Using btyper3 to assess the congruency of a proposed taxonomic framework for the Bacillus cereus group with historical typing methods. Front. Microbiol. 11, 580691. doi: 10.3389/fmicb.2020.580691

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen L., Yang J., Yu J., Yao Z., Sun L., Shen Y., et al. (2005). VFDB: a reference database for bacterial virulence factors. Nucleic Acids Res. 33 (Database issue), D325–D328. doi: 10.1093/nar/gki008

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen N. (2004). Using RepeatMasker to identify repetitive elements in genomic sequences. Curr. Protoc. Bioinf. Chapter 4, Unit 4 10. doi: 10.1002/0471250953.bi0410s05

CrossRef Full Text | Google Scholar

Cutting S. M. (2011). Bacillus probiotics. Food Microbiol. 28 (2), 214–220. doi: 10.1016/j.fm.2010.03.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Delcher A. L., Salzberg S. L., Phillippy A. M. (2003). Using MUMmer to identify similar regions in large sequence sets. Curr. Protoc. Bioinf. Chapter 10, Unit 10 13. doi: 10.1002/0471250953.bi1003s00

CrossRef Full Text | Google Scholar

Dons L., Eriksson E., Jin Y., Rottenberg M. E., Kristensson K., Larsen C. N., et al. (2004). Role of flagellin and the two-component CheA/CheY system of Listeria monocytogenes in host cell invasion and virulence. Infect. Immun. 72 (6), 3237–3244. doi: 10.1128/IAI.72.6.3237-3244.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Dover C. V. (2002). Trophic relationships among invertebrates at the Kairei hydrothermal vent field (Central Indian Ridge). Marine Biol. 141 (4), 761–772. doi: 10.1007/s00227-002-0865-y

CrossRef Full Text | Google Scholar

Ehling-Schulz M., Lereclus D., Koehler T. M. (2019). The Bacillus cereus group: Bacillus species with pathogenic potential. Microbiol. Spectr. 7 (3), 875–902. doi: 10.1128/microbiolspec.GPP3-0032-2018

CrossRef Full Text | Google Scholar

Elshaghabee F. M. F., Rokana N., Gulhane R. D., Sharma C., Panwar H. (2017). Bacillus as potential probiotics: status, concerns, and future perspectives. Front. Microbiol. 8, 1490. doi: 10.3389/fmicb.2017.01490

PubMed Abstract | CrossRef Full Text | Google Scholar

Eski A., Demir I., Gullu M., Demirbag Z. (2018). Biodiversity and pathogenicity of bacteria associated with the gut microbiota of beet armyworm, Spodoptera exigua Hubner (Lepidoptera: Noctuidae). Microb. Pathog. 121, 350–358. doi: 10.1016/j.micpath.2018.05.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Fagerlund A., Lindback T., Storset A. K., Granum P. E., Hardy S. P. (2008). Bacillus cereus Nhe is a pore-forming toxin with structural and functional properties similar to the ClyA (HlyE, SheA) family of haemolysins, able to induce osmotic lysis in epithelia. Microbiol. (Reading) 154 (Pt 3), 693–704. doi: 10.1099/mic.0.2007/014134-0

CrossRef Full Text | Google Scholar

Goebel W., Kuhn M. (2000). Bacterial replication in the host cell cytosol. Curr. Opin. Microbiol. 3 (1), 49–53. doi: 10.1016/s1369-5274(99)00050-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Goldfine H., Johnston N. C., Knob C. (1993). Nonspecific phospholipase C of Listeria monocytogenes: activity on phospholipids in Triton X-100-mixed micelles and in biological membranes. J. Bacteriol. 175 (14), 4298–4306. doi: 10.1128/jb.175.14.4298-4306.1993

PubMed Abstract | CrossRef Full Text | Google Scholar

Gomez-Escribano J. P., Castro J. F., Razmilic V., Chandra G., Andrews B., Asenjo J. A., et al. (2015). The Streptomyces leeuwenhoekii genome: de novo sequencing and assembly in single contigs of the chromosome, circular plasmid pSLE1 and linear plasmid pSLE2. BMC Genomics 16, 485. doi: 10.1186/s12864-015-1652-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Griffiths-Jones S., Moxon S., Marshall M., Khanna A., Eddy S. R., Bateman A. (2005). Rfam: annotating non-coding RNAs in complete genomes. Nucleic Acids Res. 33, 121–124. doi: 10.1093/nar/gki081

CrossRef Full Text | Google Scholar

Gu H. J., Sun Q. L., Luo J. C., Zhang J., Sun L. (2019). A first study of the virulence potential of a Bacillus subtilis isolate from deep-sea hydrothermal vent. Front. Cell Infect. Microbiol. 9, 183. doi: 10.3389/fcimb.2019.00183

PubMed Abstract | CrossRef Full Text | Google Scholar

Guinebretiere M. H., Velge P., Couvert O., Carlin F., Debuyser M. L., Nguyen-The C. (2010). Ability of Bacillus cereus group strains to cause food poisoning varies according to phylogenetic affiliation (groups I to VII) rather than species affiliation. J. Clin. Microbiol. 48 (9), 3388–3391. doi: 10.1128/JCM.00921-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Harris R. S. (2007). Improved pairwise alignment of genomic DNA. Ph.D. thesis (Penn State Univ).

Google Scholar

Hsiao W., Wan I., Jones S. J., Brinkman F. S. (2003). IslandPath: aiding detection of genomic islands in prokaryotes. Bioinformatics 19 (3), 418–420. doi: 10.1093/bioinformatics/btg004

PubMed Abstract | CrossRef Full Text | Google Scholar

Jensen G. B., Hansen B. M., Eilenberg J., Mahillon J. (2003). The hidden lifestyles of Bacillus cereus and relatives. Environ. Microbiol. 5 (8), 631–640. doi: 10.1046/j.1462-2920.2003.00461.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Jimenez G., Urdiain M., Cifuentes A., Lopez-Lopez A., Blanch A. R., Tamames J., et al. (2013). Description of Bacillus toyonensis sp. nov., a novel species of the Bacillus cereus group, and pairwise genome comparisons of the species of the group by means of ANI calculations. Syst. Appl. Microbiol. 36 (6), 383–391. doi: 10.1016/j.syapm.2013.04.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Johler S., Kalbhenn E. M., Heini N., Brodmann P., Gautsch S., Bagcioglu M., et al. (2018). Enterotoxin production of Bacillus thuringiensis isolates from biopesticides, foods, and outbreaks. Front. Microbiol. 9, 1915. doi: 10.3389/fmicb.2018.01915

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanehisa M., Goto S., Kawashima S., Okuno Y., Hattori M. (2004). The KEGG resource for deciphering the genome. Nucleic Acids Res. 32 (Database issue), D277–D280. doi: 10.1093/nar/gkh063

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim M., Oh H. S., Park S. C., Chun J. (2014). Towards a taxonomic coherence between average nucleotide identity and 16S rRNA gene sequence similarity for species demarcation of prokaryotes. Int. J. Syst. Evol. Microbiol. 64 (Pt 2), 346–351. doi: 10.1099/ijs.0.059774-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Lagesen K., Hallin P., Rodland E. A., Staerfeldt H. H., Rognes T., Ussery D. W. (2007). RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 35 (9), 3100–3108. doi: 10.1093/nar/gkm160

PubMed Abstract | CrossRef Full Text | Google Scholar

Le Calvez T., Burgaud G., Mahe S., Barbier G., Vandenkoornhuyse P. (2009). Fungal diversity in deep-sea hydrothermal ecosystems. Appl. Environ. Microbiol. 75 (20), 6415–6421. doi: 10.1128/AEM.00653-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Li W., Godzik A. (2006). Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22 (13), 1658–1659. doi: 10.1093/bioinformatics/btl158

PubMed Abstract | CrossRef Full Text | Google Scholar

Li W., Jaroszewski L., Godzik A. (2002). Tolerating some redundancy significantly speeds up clustering of large protein databases. Bioinformatics 18 (1), 77–82. doi: 10.1093/bioinformatics/18.1.77

PubMed Abstract | CrossRef Full Text | Google Scholar

Li X. Q., Tan A., Voegtline M., Bekele S., Chen C. S., Aroian R. V. (2008). Expression of Cry5B protein from Bacillus thuringiensis in plant roots confers resistance to root-knot nematode. Biol. Control 47 (1), 97–102. doi: 10.1016/j.biocontrol.2008.06.007

CrossRef Full Text | Google Scholar

Li M., Jain S., Dick G. J. (2016). Genomic and transcriptomic resolution of organic matter utilization among deep-sea bacteria in guaymas basin hydrothermal plumes. Front. Microbiol. 7, 1125. doi: 10.3389/fmicb.2016.01125

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu B., Wang Y. Q., Zhang X. B. (2006). Characterization of a recombinant maltogenic amylase from deep sea thermophilic Bacillus sp WPD616. Enzyme Microbial Technol. 39 (4), 805–810. doi: 10.1016/j.enzmictec.2006.01.003

CrossRef Full Text | Google Scholar

Liu S., Moayeri M., Leppla S. H. (2014). Anthrax lethal and edema toxins in anthrax pathogenesis. Trends Microbiol. 22 (6), 317–325. doi: 10.1016/j.tim.2014.02.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Lopes R., Cerdeira L., Tavares G. S., Ruiz J. C., Blom J., Horacio E. C. A., et al. (2017). Genome analysis reveals insights of the endophytic Bacillus toyonensis BAC3151 as a potentially novel agent for biocontrol of plant pathogens. World J. Microbiol. Biotechnol. 33 (10), 185. doi: 10.1007/s11274-017-2347-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Marteinsson V. T., Birrien J. L., Jeanthon C., Prieur D. (1996). Numerical taxonomic study of thermophilic Bacillus isolated from three geographically separated deep-sea hydrothermal vents. FEMS Microbiol. Ecol. 21 (4), 255–266. doi: 10.1111/j.1574-6941.1996.tb00122.x

CrossRef Full Text | Google Scholar

Mendez Acevedo M., Carroll L. M., Mukherjee M., Mills E., Xiaoli L., Dudley E. G., et al. (2020). Novel effective Bacillus cereus group species “Bacillus clarus“ is represented by antibiotic-producing strain ATCC 21929 isolated from soil. mSphere 5 (6), e00882–20. doi: 10.1128/mSphere.00882-20

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller R. A., Jian J., Beno S. M., Wiedmann M., Kovac J. (2018). Intraclade variability in toxin production and cytotoxicity of Bacillus cereus group type strains and dairy-associated isolates. Appl. Environ. Microbiol. 84 (6), e02479–17 . doi: 10.1128/AEM.02479-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Nishijima M., Lindsay D. J., Hata J., Nakamura A., Kasai H., Ise Y., et al. (2010). Association of thioautotrophic bacteria with deep-sea sponges. Mar. Biotechnol. (NY) 12 (3), 253–260. doi: 10.1007/s10126-009-9253-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Okaiyeto K., Nwodo U. U., Mabinya L. V., Okoh A. I. (2015). Bacillus toyonensis strain AEMREG6, a bacterium isolated from South African marine environment sediment samples produces a glycoprotein bioflocculant. Molecules 20 (3), 5239–5259. doi: 10.3390/molecules20035239

PubMed Abstract | CrossRef Full Text | Google Scholar

Osadchaia A. I., Kudriavtsev V. A., Safronova L. A., Kozachko I. A., Smirnov V. V. (1997). [Stimulation of growth and spore formation of Bacillus subtilis by optimization of carbohydrate nutrition during submerged cultivation]. Prikl. Biokhim. Mikrobiol. 33 (3), 321–324. doi: 10.1002/1098-1004(200006)15:6<577::AID-HUMU10>3.0.CO;2-Z

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramarao N., Lereclus D. (2005). The InhA1 metalloprotease allows spores of the B. cereus group to escape macrophages. Cell Microbiol. 7 (9), 1357–1364. doi: 10.1111/j.1462-5822.2005.00562.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Rojas-Solis D., Vences-Guzman M. A., Sohlenkamp C., Santoyo G. (2020). Bacillus toyonensis COPE52 modifies lipid and fatty acid composition, exhibits antifungal activity, and stimulates growth of tomato plants under saline conditions. Curr. Microbiol. 77 (10), 2735–2744. doi: 10.1007/s00284-020-02069-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruger H. J., Fritze D., Sproer C. (2000). New psychrophilic and psychrotolerant Bacillus marinus strains from tropical and polar deep-sea sediments and emended description of the species. Int. J. Syst. Evol. Microbiol. 50 (Pt 3), 1305–1313. doi: 10.1099/00207713-50-3-1305

PubMed Abstract | CrossRef Full Text | Google Scholar

Santos F. D. S., Menegon Y. A., Piraine R. E. A., Rodrigues P. R. C., Cunha R. C., Leite F. P. L. (2018). Bacillus toyonensis improves immune response in the mice vaccinated with recombinant antigen of bovine herpesvirus type 5. Benef. Microbes 9 (1), 133–142. doi: 10.3920/BM2017.0021

PubMed Abstract | CrossRef Full Text | Google Scholar

Sastalla I., Fattah R., Coppage N., Nandy P., Crown D., Pomerantsev A. P., et al. (2013). The Bacillus cereus Hbl and Nhe tripartite enterotoxin components assemble sequentially on the surface of target cells and are not interchangeable. PLoS One 8 (10), e76955. doi: 10.1371/journal.pone.0076955

PubMed Abstract | CrossRef Full Text | Google Scholar

Schattner P., Brooks A. N., Lowe T. M. (2005). The tRNAscan-SE, snoscan and snoGPS web servers for the detection of tRNAs and snoRNAs. Nucleic Acids Res. 33 (Web Server issue), W686–W689. doi: 10.1093/nar/gki366

PubMed Abstract | CrossRef Full Text | Google Scholar

Stenfors Arnesen L. P., Fagerlund A., Granum P. E. (2008). From soil to gut: Bacillus cereus and its food poisoning toxins. FEMS Microbiol. Rev. 32 (4), 579–606. doi: 10.1111/j.1574-6976.2008.00112.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun Q. L., Wang M. Q., Sun L. (2015). Characteristics of the cultivable bacteria from sediments associated with two deep-sea hydrothermal vents in Okinawa Trough. World J. Microbiol. Biotechnol. 31 (12), 2025–2037. doi: 10.1007/s11274-015-1953-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Tallur P. N., Sajjan D. B., Mulla S. I., Talwar M. P., Pragasam A., Nayak V. M., et al. (2016). Characterization of antibiotic resistant and enzyme producing bacterial strains isolated from the Arabian Sea. 3 Biotech. 6 (1), 28. doi: 10.1007/s13205-015-0332-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Tatusov R. L., Fedorova N. D., Jackson J. D., Jacobs A. R., Kiryutin B., Koonin E. V., et al. (2003). The COG database: an updated version includes eukaryotes. BMC Bioinf. 4, 41. doi: 10.1186/1471-2105-4-41

CrossRef Full Text | Google Scholar

Tonello F., Zornetta I. (2012). Bacillus anthracis factors for phagosomal escape. Toxins (Basel) 4 (7), 536–553. doi: 10.3390/toxins4070536

PubMed Abstract | CrossRef Full Text | Google Scholar

van Asten F. J., Hendriks H. G., Koninkx J. F., van Dijk J. E. (2004). Flagella-mediated bacterial motility accelerates but is not required for Salmonella serotype Enteritidis invasion of differentiated Caco-2 cells. Int. J. Med. Microbiol. 294 (6), 395–399. doi: 10.1016/j.ijmm.2004.07.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Dover C. L. (2014). Impacts of anthropogenic disturbances at deep-sea hydrothermal vent ecosystems: a review. Mar. Environ. Res. 102, 59–72. doi: 10.1016/j.marenvres.2014.03.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang W., Li Z., Zeng L., Dong C., Shao Z. (2020). The oxidation of hydrocarbons by diverse heterotrophic and mixotrophic bacteria that inhabit deep-sea hydrothermal ecosystems. ISME J. 14 (8), 1994–2006. doi: 10.1038/s41396-020-0662-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang R., Zhang C. T. (2003). Identification of genomic islands in the genome of Bacillus cereus by comparative analysis with Bacillus anthracis. Physiol. Genomics 16 (1), 19–23. doi: 10.1152/physiolgenomics.00170.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao Y., Chen C., Gu H. J., Zhang J., Sun L. (2019). Characterization of the genome feature and toxic capacity of a Bacillus wiedmannii isolate from the hydrothermal field in okinawa trough. Front. Cell Infect. Microbiol. 9, 370. doi: 10.3389/fcimb.2019.00370

PubMed Abstract | CrossRef Full Text | Google Scholar