Composite genome sequence of Bacillus clausii, a probiotic commercially available as Enterogermina®, and insights into its probiotic properties

Composite genome sequence of Bacillus clausii, a probiotic commercially available as Enterogermina®, and insights into its probiotic properties | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research article Composite genome sequence of Bacillus clausii, a probiotic commercially available as Enterogermina®, and insights into its probiotic properties Indu Khatri, Gaurav Sharma, Srikrishna Subramanian This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.2.15490/v3 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 30 Dec, 2019 Read the published version in BMC Microbiology → Version 3 posted 2 You are reading this latest preprint version Show more versions Abstract Background: Some of the spore-forming strains of Bacillus probiotics are marketed commercially as they survive harsh gastrointestinal conditions and bestow health benefits to the host. Results: We report the composite genome of Bacillus clausii ENTPro from a commercially available probiotic Enterogermina® and compare it with the genomes of other Bacillus probiotics. We find that the members of B. clausii species harbor high heterogeneity at the species as well as genus level. The genes conferring resistance to chloramphenicol, streptomycin, rifampicin, and tetracycline in the B. clausii ENTPro strain could be identified. The genes coding for the bacteriocin gallidermin, which prevents biofilm formation in the pathogens S. aureus and S. epidermidis, were also identified. KEGG Pathway analysis suggested that the folate biosynthesis pathway, which depicts one of the important roles of probiotics in the host, is conserved completely in B. subtilis and minimally in Bacillus clausii and other probiotics. Conclusions: We identified various antibiotic resistance, bacteriocins, stress-related, and adhesion-related domains, and industrially-relevant pathways, in the genomes of these probiotic bacteria that are likely to help them survive in the harsh gastrointestinal tract, facilitating adhesion to host epithelial cells, persistence during antibiotic treatment and combating bacterial infections. General Microbiology Bacteriocins Gastrointestinal-tract Phylogeny Resistome Pathogenicity. Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Background Probiotics are live microbes which when consumed in sufficient amount helps to resume the original gut microflora, distressed by diarrhea or antibiotic intake [1]. Most bacterial probiotics such as Lactobacillus and Bifidobacteria, which are inhabitants of the gut, are available as lyophilized preparations of vegetative cells while some probiotic bacterial preparations that belong to the genus Bacillus are available in the form of spores [2]. Bacterial spores are dormant and resistant to heat, desiccation, dehydration and are extremely stable, which is a desirable property for probiotics [3, 4]. The spores of Bacillus germinates in the gut and the vegetative cells are vital for the human gut health [5]. Bacillus subtilis belongs to one of the most studied and explored family Bacillaceae and many of its strains are being used as probiotics since the 1990s. The role of Bacillus species ranges from probiotic nature of B. subtilis , B. clausii , B. coagulans , B. pumilus and other strains to biological control agents ( B. thuringiensis and B. sphaericus ), and pathogenicity ( B. anthracis and B. cereus ). Several strains are economically important ( B. subtilis ) whereas others have medical importance ( B. licheniformis ) [6]. Some Bacillus spp. are industrially-important and produce proteins such as alkaline proteases, xylanases, amylases, and cellulases [4, 7]. In 2001, some B. subtilis strains, which were used in the probiotics and soap industry, were reclassified as B. clausii [8]. B. clausii spores are marketed as the probiotic Enterogermina ® which consists of four Bacillus strains O/C, N/R, SIN, and T that are resistant to Chloramphenicol, Novobiocin/Rifampicin, Neomycin/Streptomycin, and Tetracycline, respectively [9, 10]. Although these four strains are known to have been derived from a single penicillin-resistant strain, B. subtilis ATCC 9799 [9, 11], secretome analysis have revealed variation in the expression level of some of the secreted proteins [12]. Also, the O/C strain of B. clausii inhibits the cytotoxic effect induced by the Clostridium difficile and B. cereus toxins [13]. The intrinsic antibiotic resistance in probiotics is considered advantageous in cases of antibiotics-probiotics combination prescriptions to restore healthy gut [14, 15]. The mode of action of B. clausii as a probiotic is not clear, but the strains have been reported to secrete some proteins that are involved in the immunomodulatory mechanism, adaptation and their colonization in the human gastrointestinal tract (GIT) [12, 13, 16–18]. Uniquely, B. clausii harbor erm(34) gene that imparts the resistance to erythromycin. The erm(34) gene is not a homolog of erm(A) , erm(B) , erm(C) and erm(TR) genes for MLSB resistance in Gram-positive human pathogens and erm(D) , erm(K) , and erm(J) characterized in B. licheniformis , B. halodurans , and B. anthracis, respectively [19]. Various clinical trials and molecular studies [8, 13, 16, 20–23] have been performed to identify the major features that demarcate B. clausii probiotic strains from other Bacillus spp., but still, the genomic reasons of its probiotic activity have not been reported before. Therefore, we sequenced the composite genome sequence of Bacillus clausii (composite of all four strains of B. clausii used in the probiotic formulation) from Enterogermina ® , an oral probiotic, marketed by Sanofi in India. The composite genome obtained from the sequencing of this probiotics was named as B. clausii ENTPro. We performed an extensive analysis to identify the genomic features known to impart probiotic properties in B. clausii viz., adhesion to gut, withstanding harsh conditions in the gut, antibiotic resistance, and biosynthesis pathways. In addition, to gain insight into the genomic features of different probiotics, we have compared pathways, types of bacteriocins and antibiotic resistance genes in different Bacillus probiotics. Results Genome features of B. clausii ENTPro De novo assembly of PacBio sequencing reads of B. clausii ENTPro gDNA resulted in two contigs: one long circular contig of 4,264,866 base pairs (bp) and one short circular 31,475 bp contig. The long contig represents the composite circular chromosome of Bacillus clausii ENTPro with an average GC content of 44.75% (Fig. 1 ) and the smaller one (GC Content: 39.9%) is likely a plasmid. In addition, Illumina sequencing-based assembly resulted in 4.3 Mbp genome from 36 contigs and N50 of 344,696 bp, which overlaps completely with the genome assembled using PacBio reads. The composite genome obtained from PacBio sequencing reads was submitted to GenBank [NC_006582.1] and further used for all the comparisons in this study. Bacillus clausii ENTPro genome is 99.8% similar to another probiotic strain B. clausii B106 [NFZO01] ( Fig. S1A ), followed by 94.3% similarity to B. clausii KSM-K16 [NC_006582.1] ( Fig. S1B ), whereas other members of the same species are 50-94% similar. This suggests that the members of this species are quite diverse as characterized by their GGDC values ( Table S1 ). Our analysis suggests that probiotic strains within B. clausii such as ENTPro, B106, and UBBC-07 are highly similar to each other as compared to other strains. The plasmid sequence is novel and does not have any close similarity with other plasmids in the NCBI nucleotide database (NT). Most of the proteins encoded by the plasmid sequence are hypothetical and are not functionally characterized. We mapped Illumina reads against the plasmid database downloaded from NCBI to identify if we could obtain hits to any previously known plasmids. Very few reads mapped on to known plasmids and no full plasmid could be retrieved using the Illumina reads. Therefore, we concluded that the identified plasmid sequence harbored by B. clausii ENTPro is novel. Annotation of the B. clausii ENTPro genome revealed the presence of 4,384 protein-coding sequences, which constitute 86.73% of the genome with an average length of 843 bp (ranging from 113 to 9,509 bp) (Table 1) . A total of 1,215 Coding DNA Sequences (CDS) were annotated as hypothetical proteins, accounting for 27.72% of the total proteins. The ENTPro genome has all the three proteins R (restriction), M (modification), and S (specificity) that belongs to the Type I RM system. m6A methylation was observed in >96% of the motifs G A GNNNNNNRTGC and GC A YNNNNNNCTC in the genome at 2 nd and 3 rd positions, respectively. There are 75 tRNA genes and seven complete rRNA operons (>99% identity) in the B. clausii ENTPro genome. 16S rRNAs obtained from the de novo assembly of B. clausii ENTPro genome shows 99.8% similarity with B. clausii Enterogermina strains O/C, T, N/R, and SIN. This is in line to previously known variations in 16S rRNA genes in bacterial genomes [24]. Most of the varying sites were present in the V1 region of the 16S rRNA sequences even in B. clausii KSM-K16 and B. clausii DSM 8716 (Fig. S2) . Amongst the total proteome, ~75% (3,311) proteins could be categorized into Clusters of Orthologous Groups (COGs) functional groups. Among these mapped proteins, ~35% belonged to the metabolism category, ~14% to cellular processes and signaling and ~16% proteins to information storage and processing. According to COG mapping data, 152 proteins are involved in signal transduction mechanisms (COG: T) and 44 proteins were reported to function in secondary metabolites biosynthesis, transport, and catabolism (COG: Q). COG assignments to proteomes of B. clausii members revealed that all the organisms have similar number of proteins assigned to various COG categories ( Fig. 2 ). Phylogenetic position of B. clausii as inferred from housekeeping proteins-based phylogeny Phylogenetically, B. clausii clustered in a separate clade with further grouping within this clade ( Fig. 3 ). The phylogenetic tree reveals that ENTPro strain is closest to the B106 strain of B. clausii . Both these probiotic strains are further similar to another probiotic strain UBBC-07 of B. clausii . All these probiotic strains share a common ancestor with industrial B. clausii KSM-K16 strain. This phylogenetic placement of B. clausii probiotic strains is concordant with the whole genome similarity matrix as obtained by genome-genome distance calculator (GGDC) [25]. Other B. clausii “Heroin” strains form several different groups within the B. clausii clade. Interestingly, the B. clausii proteome matches the proteome of other Bacillus species at <70% identity. This clearly suggests the genomic heterogeneity of B. clausii genome in comparison to other Bacillus species. We also included all Bacillus probiotics genomes in phylogenetic analysis to investigate their position phylogenetically [26]. Bacillus probiotics shared clades with their species members. Interestingly, probiotic strains cluster together e.g. B. clausii , B. coagulans and B. subtilis . B. clausii ENTPro as a derived strain from four different strains B. clausii Enterogermina ® is a mixture of four different strains each of which is supposed to confer resistance against specific antibiotics, namely novobiocin and rifampicin (strain N/R), chloramphenicol (strain O/C), streptomycin and neomycin (strain SIN) and tetracycline (strain T) [12]. The specific genes conferring resistance could not be traced in the literature so different in silico strategies were employed to identify possible genes that could help impart resistance to these antibiotics in B. clausii ENTPro ( Table S2 and S3) . Rifampicin: Rifampicin resistance is acquired by specific mutations at positions 516, 526 and 531 in the rpoB gene in Escherichia coli [27]. These mutations are mapped in the center of the rpoB gene in 3 regions: one cluster covering 507-533 amino-acid (AA); cluster II covering AA 563–572 and cluster III with AA change at position 687, which altogether are referred to as RIF resistance determining region (RRDR) [27]. In order to find the presence of RRDR region in RpoB protein in ENTPro, the RpoB protein sequences from all Bacillus spp. were retrieved and aligned with E. coli RpoB protein sequence [Accession Number: NP_418414.1]. P 524 ->L (corresponding to 567 AA position in E. coli RpoB protein sequence) AA change was observed in B. clausii ENTPro strain that was not observed in other Bacillus spp. ( Fig. S3 ). Chloramphenicol: Chloramphenicol acetyltransferase, involved in conferring resistance against chloramphenicol [28], was identified from the proteome of B. clausii ENTPro [Accession Number: WP_035203840.1]. Streptomycin: Pfam domains, known to impart resistance against streptomycin, were identified in B. clausii ENTPro. Nine proteins in B. clausii had the Pfam domain PF02522, PF01636, PF01909, PF04439, PF04655, PF07091, PF07827, and PF10706 that has core domain aminoglycoside. Two proteins had streptomycin adenylyltransferase domain (PF04439), six proteins have aminoglycoside phosphotransferase [PF01909] domain and one protein has Kanamycin nucleotidyltransferase [PF07827] domain in their sequence ( Table S2 ). This suggests the presence of domains that are involved in imparting resistance to streptomycin. In addition, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis of the organism reveal the presence of complete KEGG pathway for the streptomycin biosynthesis in the B. clausii ENTPro ( Fig. S4 ). Tetracycline: The domains conferring tetracycline resistance [RF0133, RF0134, RF0135, and RF0127] were present in B. clausii ENTPro ( Fig. 4 ). The presence of these genes in the composite genome of B. clausii ENTPro was further confirmed by mapping the Illumina reads to these genes. Probiotic Properties in B. clausii ENTPro Probiotics are beneficial components of microbiota that modulates immunological, respiratory and gastrointestinal functions [29]. For imparting these functions probiotics adhere to the mucosal membrane to interact with the host and have acidic, alkaline and oxidative stress resistance and stress adaptation proteins [6]. Probiotics are believed to have good adherence capacity which promotes gut residence time, pathogen elimination and adhesion to the epithelial layer of host cells and exerting immune modulation. Pfam analysis reveals the presence of three proteins involved in adhesion namely a mucus-binding protein with ‘Gram_pos_anchor’ Pfam domain [PF00746] at the C-terminus, a collagen-binding protein with LPXTG motif at the C-terminus and a fibronectin-binding protein [30] (Table S2) . These adhesion proteins may help facilitate the probiotic bacterium to bind and help in the direct interactions with the intestinal mucosa layer. Probiotic B. clausii has to encounter various harsh environmental conditions during transit in the GIT such as the acidic environment in the stomach, bile juice environment in the small intestine, oxidative stress, and osmotic stress [1]. When a bacterium faces an acidic environment, H + homeostasis is maintained by F0F1 ATP synthase pump, which work by hydrolyzing ATP to pump protons (H + ) from the cytoplasm [1, 31]. We found that this synthase complex is present in ENTPro genome as a full operon [DB29_02342--DB29_02349]. The bacteria have to face the toxicity of bile salts that induce intracellular acidification and act as detergents that disrupt biological membranes [32]. Five proteins were identified that were involved in bile tolerance mechanism; two belong to ornithine decarboxylase [33] and three to sodium bile acid symporter family [34, 35] (Table S2) . B. clausii ENTPro also harbors general stress adaptation proteins. The universal stress protein UspA [PF00582] is important for survival during cellular growth arrest and reprograms the cell towards defense and escape during cellular stress [36, 37]. Molecular chaperones that may impart resistance against environmental stress were obtained through annotation and Pfam domain search such as the chaperonin GroES [PF00166] and GroEL [38, 39] and one heat shock protein 33 [PF01430], two copies of cold shock proteins CspA [PF00313], three Clp protease [PF00574] and HtpX and HrcA-like heat shock proteins. These proteins play an important role in basic cellular functions that includes growth, the stability of DNA and RNA and they also prevent the formation of inclusion bodies [40–42]. For hyperosmotic stress and heat resistance, B. clausii ENTPro harbors one copy each of the chaperone protein DnaJ [PF00226] and nucleotide exchange factor GrpE [PF01025]. Also, two methionine sulfoxide reductase A [43] [PF01625] were present in B. clausii ENTPro that provides resistance in oxidative stress (Table S2) . This suggests that B. clausii ENTPro has proteins to improve adhesion and handling stress and harsh conditions in the human gut. Antibiotic Resistance in Bacillus Probiotics Antibiotic resistance is a common phenomenon in Gram-positive bacteria [44–46]. It is accomplished by genes acquired either horizontally through plasmids, or foreign DNA recombination, or mutations at different chromosomal loci in the bacterial genome [47]. It is preferred that probiotic strains carry few antibiotic resistance genes as possible so that they are not a putative source for transferring these genes to other gut bacteria including pathogens [46]. However, on the other hand since some of these probiotics are administrated alongside antibiotics, some resistance to commonly administrated antibiotics are desirable. Presence of a novel plasmid sequence in B. clausii ENTPro could be a possible source of antibiotic-resistance gene transfer but we could not identify any potential antibiotic-resistance domain(s) in the plasmid. We also searched for the presence of antibiotic resistance genes and efflux pumps in the genomes with multiple methods to avoid false positives. The Chloramphenicol acetyltransferase, that confers resistance against chloramphenicol, is absent in B. amyloliquefaciens and B. coagulans whereas chloramphenicol efflux pump was present in B. amyloliquefaciens ( Fig. 4 ). This would imply the presence of chloramphenicol resistance in all the Bacillus probiotics except B. coagulans . Different classes of beta-lactamase were present in one or the other Bacillus probiotics that clearly suggest the presence of resistance against Penicillin in all the Bacillus probiotics. Multidrug resistance protein, a universal stress protein, EmrB, and its efflux pump, tetracycline resistance protein, and penicillin-binding protein are present in all the Bacillus probiotics. This suggests that most of the Bacillus probiotics are resistant to common antibiotics. Erythromycin resistance was identified by subjecting the erm (34) gene sequence (GenBank Identifier: AY234334) of B. clausii DSM8716 to BLASTn against all Bacillus genomes. This gene was identified in B. clausii ENTPro named as “SSU rRNA (adenine(1518)-N(6)/adenine(1519)-N(6))-dimethyltransferase” (GenBank Identifier: ALA53582). The gene was also identified in all the B. clausii genomes. The gene sequence shared 61% identity to rRNA adenine methyltransferase of B. halodurans and 57% identity to rRNA adenine methyltransferase of B. licheniformis , B. anthracis , B. sonorensis and B. fordii . The rRNA adenine methyltransferase gene from other Bacillus spp. shared 20-50% identity with erm (34) gene. The result reveals that the erm (34) gene is unique to B. clausii and is not present in other members of the Bacillus genus. Vancomycin resistance, as observed from KEGG pathway analysis, ( Fig. S5 ) was identified only in B. toyonensis while absent in other Bacillus probiotics. The accessory proteins of vancomycin resistance operon were present in some of the Bacillus probiotics, but resistance-conferring genes were completely absent. We would like to add an advisory note that previous studies have shown that an organism may exhibit intrinsic resistance to a few antibiotics that could not be related to its genotype [46]. Though we have endeavored to relate the genome-level occurrence of antibiotic resistance proteins or domains to their probable phenotypes, we have not performed any phenotypic studies to substantiate these analyses and/or confirm for intrinsic resistance. Further, the current situation may constitute a safety concern because of the possibility of transfer of antibiotic gene transfer to gut flora [46]. Bacteriocins in Bacillus probiotics Bacteriocins are proteinaceous toxins produced by bacteria that act as narrow-spectrum antibiotics to inhibit the growth of similar or closely related bacterial strains [48, 49]. They can help probiotics to survive the toxins produced by invading bacteria by inhibiting their growth and hence can result in beneficial effects on the hosts. The identified bacteriocins in all the probiotics are represented in a presence-absence binary matrix in Fig. 5 . Several of these bacteriocins are already well utilized in therapeutics [50] and their spectrum against pathogens is well established [9, 51–53]. Gallidermin identified via in silico analysis in B. clausii genomes is known to efficiently prevent biofilm formation in the pathogens S. aureus and S. epidermidis species [26]. This bacteriocin has also been reported to be effective in skin disorders including acne, eczema, folliculitis, and impetigo where the targets organisms are Propionibacteria, Staphylococci, and Streptococci [50]. Lacticin 3147 A2 and Leucocyclin Q as identified in B. amyloliquefaciens are broad-spectrum bacteriocins. Lacticin has been used effectively in the treatment of bacterial mastitis, Staphylococcal and Enterococcal infections including vancomycin-resistant Enterococci [50] and is effective against Listeria infections [51]. Similarly, leucocyclicin Q exhibit bactericidal or bacteriostatic effects on Gram-positive bacteria, including food-borne pathogens, such as Lactococcus , Weissella paramesenteroides , Pediococcus dextrinicus , Enterococcus , Streptococcus , and Leuconostoc [52]. Plantazolicin identified in B. amyloliquefaciens and B. pumilus has nematicidal activity [54]. Cirucularin A produced by B. coagulans has been reported to be the most effective bacteriocin against C. tyrobutyricum NIZOB570, a known cheese-spoilage bacterium [55] and also Lactococci , Enterococci, and some Lactobacillus strains [56]. LichenicidinVK21A2 identified in B. paralicheniformis is considered as self-immunity bacteriocin that exhibits antimicrobial activity against several strains of Listeria monocytogenes , methicillin-resistant S. aureus , and vancomycin-resistant Enterococcus [57]. Zoocin A in B. toyonensis shows antimicrobial activity against several other Streptococci by cleaving the peptidoglycan cross-links of the target cell wall [58]. Subtilosin A produced by B. subtilis is also a broad range bacteriocin that is effective against Listeria monocytogenes , and strains of E. faecalis , P. gingivalis , K. rhizophila , Enterobacter aerogenes , Streptococcus pyogenes, and Shigella sonnei [53]. Sporulation-killing factor skfA produced by B. subtilis induces the lysis of other B. subtilis cells that have not entered the sporulation pathway. This cannibalistic behavior provides a source of nutrients to support those cells that have entered sporulation [59, 60]. At high concentrations, it can also inhibit the growth of other bacteria [61]. The presence of well-characterized bacteriocins in the Bacillus probiotics suggests their important role in fighting against the pathogen in the gut. Folate Biosynthesis Pathways in Bacillus Probiotics The gut microbiota aids the host, playing a crucial role in nutrient digestion and energy recovery. Due to potentially relevant applications, the capacity to yield folate has been investigated in various probiotic strains. Previously, the presence of these pathways was reported in Lactobacillus and Bifidobacterium probiotics but was not explored in Bacillus probiotics [62] except B. subtilis [63]. We performed the identification of key components of folate production pathways in Bacillus probiotics using KEGG Pathway database [64]. The analysis of genome sequences of Bacillus probiotics revealed the presence of complete operon to synthesize para-aminobenzoic acid (PABA) de novo only in B. subtilis probiotics ( Fig. 6). On the other hand, the enzymes, necessary for chorismate conversion into PABA are present in almost all the Bacillus probiotics. Moreover, the shikimate pathway for chorismate production is complete only in B. subtilis , B. pumilus and B. toyonensis , while it is partial in all the other Bacillus probiotics. On the other hand, Bacillus probiotic strains contain the genes of DHPPP de novo biosynthetic pathway, the gene encoding dihydropteroate synthase (EC 2.5.1.15) and gene encoding dihydropteroate synthase (EC 2.5.1.15). Therefore, it is expected that these strains are not auxotrophic for folates or DHP but can produce folate in the presence of PABA supplementation. The presence/absence of the components of the folate biosynthesis pathway is reported based on KEGG pathway analysis. Previous studies have revealed that B. subtilis genome harbor all the pathways components and have been engineered for folate production [63, 65, 66]. Discussion In this study, we report the complete composite genome of Bacillus clausii ENTPro, sequenced from a commercially available probiotic, Enterogermina ® that is a mixture of four closely related strains. The oral probiotics are popularly recommended by physicians as an adjunct to antibiotic therapy to avoid antibiotic-induced diarrhea and/or other gastrointestinal distress. We determined the composite genome of the probiotic Bacillus clausii , to identify the features responsible for its probiotic properties and correlate it to the phenotypic properties mentioned previously in the literature [26, 50, 61, 67–71]. While literature studies also mention the presence of four closely related strains in the probiotic Enterogermina ®, [2, 12, 72] our single chromosome assembly suggests that the different strains are inseparable. However, variations in number of bacteria from each strain in the formulation has also been reported previously [73]. Our investigation revealed the presence of antibiotic resistance genes that B. clausii ENTPro harbors. However, we could not find the genes conferring Novobiocin and Neomycin resistance, possibly because the genome reported here is a composite of the four strains. The possibility of the identified genome to be strain-specific could explain these missed antibiotic resistance genes. Previously reported secretome analysis [12] supports the presence of four related strains, and our genome data suggests that they are very closely related. COG assignments of all the B. clausii strains did not show any differences in probiotic strains versus industrial strains of B. clausii . Apart from these common features within the members of B. clausii genus, the low genome-genome distance revealed that its members are quite diverse. We further used phylogenetic methods to understand the relationship of the members of B. clausii and other Bacillus using representatives of each species. The phylogenetic tree clearly has a separate group where members of B. clausii lies with B. lehensis , B. halodurans whereas B. subtilis, B. coagulans formed separate clades. The B. clausii proteome analysis also supported its distant placement in the phylogenetic tree from other Bacillus members. These results indicate that B. clausii is unique in comparison to other Bacillus species. Using in silico mining approaches, we previously reported the sporulation proteins as well as various other proteins that might play a role in probiotic function such as molecular chaperones, stress proteins, flagellin, and mucin binding protein in two other Bacillus probiotics marketed in India namely B. coagulans S-lac and B. subtilis TOA JPC [67]. We also reported the presence of adhesins, which might aid in adhesion to the mucosal layer of host tissues. All these proteins are present in the currently compared probiotic genomes as well. In addition to these, in the current study, we looked for the antibiotic resistance genes, bacteriocins, and folate biosynthesis pathways. Multiple strategies were used in our analysis to comprehensively catalog these domains, nevertheless, these bacteria may harbor genes that confer intrinsic antibiotic resistance. For example, the presence of chloramphenicol, tetracycline and vancomycin resistance in B. toyonensis has been reported previously [68, 69]. The bile tolerance is reported for B. clausii Enterogermina in a recent study [74] that substantiates our findings related to the genes exhibiting bile tolerance. We have also explored the occurrence of bacteriocins in different Bacillus probiotics. The presence of different bacteriocins makes these organisms unique. For example, B. clausii would be effective in S. aureus infections [26] in skin disorders [50] while B. amyloliquefaciens in food industry especially yogurt beverage and probiotic [70], B. toyonensis shows antimicrobial activity against several other Streptococci [58, 75, 76]. One of the important roles of probiotics is nutrient digestion and energy recovery by producing folate. Folate occurs naturally in food [62] and it is required for the efficiency of DNA replication, repair and methylation process in humans [77, 78]. The folate-producing probiotic strains could possibly confer protection against cancer, inflammation, cardiovascular disease, stress, and depression [62, 63, 77, 78]. This role of probiotics has been extensively studied for their commercial utilization in folate production [62, 79]. We investigated the pathways involved in the production of folic acid in these probiotics. Not surprisingly, we found all the pathway components intact in B. subtilis which is being engineered for the enhanced synthesis of folate [63]. Some core components of the folate synthesis pathways were present in other Bacillus probiotics, suggesting that they may also be potential sources for the de novo synthesis of folate. The composite genome of Bacillus clausii ENTPro and the comparative analysis presented here have helped identify several genes and pathways of interest required for probiotic action and can serve as a starting point for the experimental characterization of these gene products and bacteriocins in order to obtain a deeper understanding into the mechanism of probiotic action of these important microbes. Conclusions The composite circular genome of Bacillus clausii strain ENTPro, isolated from Enterogermina ® , an oral probiotic, marketed by Sanofi in India, is reported. The genomes of different strains are inseparable as suggested by complete circular assembly using long PacBio reads. B. clausii ENTPro shares high similarity with probiotic strains of genus B. clausii whereas it is quite diverse as compared to other Bacillus probiotics. The ML tree based on 25 housekeeping protein sequences, clearly places Bacillus clausii in a separate clade as an outgroup of the Bacillus species. In this study, we report the genes that are responsible for conferring antibiotic resistance in B. clausii . We could identify all the antibiotic resistance genes that are indicative of the presence of all four strains in the assembled composite genome. Also, we compared the presence of antibiotic resistance-conferring genes and related pathways in all the probiotic Bacillus genomes. The most important finding of our study is the identification of bacteriocins in Bacillus probiotic genomes which could be directly related to their usage in food and beverage industry. For e.g. gallidermin bacteriocin identified in Bacillus clausii functions against S. aureus biofilm formation and infections [26, 80]. The bacteriocins in B. amyloliquefaciens fight against foodborne pathogens [81] which clearly indicates and supports its usage in the yoghurt beverage industry. B. paralicheniformis secretes bacteriocins to prevent Listeria , S. aureus, and Enterococcus borne infections [57, 82]. The other important aspect we studied was the presence of genes necessary for the production of folate. We found that B. subtilis can produce folate de novo whereas other Bacillus probiotics depend on supplements viz. pABA to produce the same. Several important components and alternative pathways for folate production were present in other Bacillus probiotics but not complete like B. subtilis . While identifying several genes and pathways of interest is insufficient to explain the concerted probiotic action, we believe our study shed light on several genomic aspects of different Bacillus probiotics. We trust that the comparative genomics analysis presented here will pave the way for experimental characterization of our findings, and to possibly engineer these organisms for enhanced probiotic actions. Methods Isolation and purification of Bacillus clausii genomic DNA: B. clausii spore suspension drug “Enterogermina ® , Sanofi-Aventis” (Batch No. 120965; Mfd. 12/2011 and Exp. 11/2013) was procured from a drugstore in Chandigarh, India and was cultured in March 2013. Bacterial cells were suspended in Milli-Q water, serially diluted, and plated on ATCC medium: 688 nutrient agar plates. The plates were incubated at 25°C for 48 hours. DNA isolation was performed using the ZR Fungal/Bacterial DNA miniprep kit (Zymogen) as per instructions in its user manual. After isolation, the genomic DNA was treated with RNase A (1µl of a 10 µg/mL stock solution for 100µl of a solution containing DNA) and incubated at 37ºC for 30 minutes. Then, 1/10 volume of 3M sodium acetate (pH 5.2) and 2.5 volumes of absolute ethanol was added followed by incubation at -20ºC overnight and centrifugation at 14,000 rpm for 30 minutes at 4ºC. The supernatant was carefully discarded; the pellet was rinsed with 70% ethanol and centrifuged again at 14,000 rpm for 15 minutes at 4ºC. The ratio of OD at 260/280 nm was >1.8 as observed by NanoDropND-1000 spectrophotometer. Genome Sequencing: PacBio Sequencing: The probiotic B. clausii ENTPro was sequenced using PacBio P6C4 chemistry at Genome Quebec Centre, McGill University. DNA samples were sheared and concentrated using AMPure magnetic beads and treated by ExoVII to remove single-stranded ends. SMRTbell libraries were created using the ‘Procedure and Checklist–20 kb Template Preparation Using BluePippin TM Size Selection System protocol. Size Selection was performed to retain longer reads (>10k reads) for sequencing. Blunt ligation reactions were prepared and SMRTbell templates were purified using AMPure magnetic beads. BluePippin TM The Size-selected SMRTbell templates were annealed and polymerase was added for Sequencing. Single SMRT cell was run on the PacBio RS II system using P6C4 chemistry and a 180-minute data collection mode. Illumina Sequencing: B. clausii ENTPro was also sequenced using the Illumina HiSeq PE platform. The library preparation was carried out according to the TruSeq DNA sample preparation protocol (Illumina, Inc., San Diego, CA) at C-CAMP, Bangalore, India. One μg of bacterial DNA was sheared to an average length of 300 to 400 bp. End repair, A-tailing, and adapter ligation (~120 base adapter) procedure was performed according to paired-end DNA sample preparation kit (Illumina, Size selection of adapter-ligated DNA was done in a range of 400 to 550 bases for DNA library). The insert size was taken in a range of 280 to 430 bases for DNA library. PCR enrichment was performed for eight cycles, and the samples were validated on a Bioanalyzer. Libraries were sequenced in a paired-end 100 base run, using TruSeq PE Cluster Kit v3-cBot-HS for cluster generation on C-bot and TruSeq SBS Kit v3-HS (Catalog No.: PE-401-3001) for sequencing on the Illumina HiSeq 1000 platform according to recommended protocols. Genome Assembly and Annotation: The PacBio reads were assembled de novo using Hierarchical Genome Assembly Process (HGAP) v2.0 [71] in SMRT portal using default parameters. Functional annotation was carried out by RAST (Rapid Annotation using Subsystem Technology) [83, 84], tRNA was predicted by tRNAscan-SE 1.23 [85] and rRNA genes by RNAmmer 1.2 [86]. The taxonomic characterization of the contigs was performed by subjecting the contigs to BLASTn [87] against NT database. The methylome was deduced by RS Modification and Motif analysis in SMRT portal ( https://github.com/PacificBiosciences ). The plasmid sequence was confirmed by plasmidSPAdes [88]. Phylogenetic Analysis The Bacillus genus comprises more than 1000 genomes. To be efficiently able to plot the phylogenetic position of Bacillus clausii in Bacillaceae family, we retrieved the Bacillus genomes that have been classified as representative genomes by NCBI (Accession Numbers in Table S4 ). The representative genomes from genus Listeria and Clostridium were selected as outliers. Twenty five housekeeping proteins (ribosomal protein S11, ribosomal protein S13, ribosomal protein S19, ribosomal protein S2, ribosomal protein S3, ribosomal protein S5, ribosomal protein S9, ribosomal protein L11, ribosomal protein L13, ribosomal protein L19, ribosomal protein L2, ribosomal protein L20, ribosomal protein L27, ribosomal protein L3, ribosomal protein L4, ribosomal protein L5, ribosomal protein L6, ribosomal protein L7/L12, CTP synthase, DNA gyrase subunit B, DNA mismatch repair protein MutS, DNA primase, elongation factor Ts, elongation factor Tu, Transcription termination protein NusA) were retrieved from all the genomes and were concatenated in a particular order. The sequences were aligned using Muscle [89] and the phylogenetic inference was drawn using Maximum likelihood [ML] approach based on PROTGAMMA model in RAxML [90] (Bootstrap: 100). Comparative Genomics For comparative analysis, all B. clausii genomes available on March 2018 were downloaded from NCBI. In addition, the genomes with reported probiotics properties were downloaded from the NCBI for comparison. All these genomes were annotated again using the RAST server [83] to remove the bias from different annotation strategies. COGs were identified by subjecting the proteomes of these organisms to BLASTp against COG database [91] at E-value 1e-5. Identification of genome features contributing to probiotic properties of B. clausii ENTPro B. clausii ENTPro was scanned using Hidden Markov Model (HMM) [92, 93] for the presence of specific domains involved in acid tolerance, Adhesion, antibiotic resistance, antimicrobial production, heavy metal resistance, bile resistance, oxidative and universal stress resistance, and riboflavin synthesis. Identification of Antibiotic Resistance genes The Comprehensive Antibiotic Resistance Database (CARD) [94] and Pfam domains were downloaded and hmmscan [93] was run locally against the proteome of all the organisms to identify the domains that could impart antibiotic resistance. Rifampicin resistance was identified based on the mutation in RpoB genes. Chloramphenicol resistance was identified based on the presence of chloramphenicol acetyltransferase gene in the proteome of the respective organisms. The presence of erm (34) was found by subjecting its gene sequence to BLASTn against the NR database. Streptomycin biosynthesis and Vancomycin resistance were identified from the pathway analysis against the KEGG database [64]. Identification of Bacteriocins Bacteriocins were reported as per the identification from the BAGEL3 server [95]. Any bacteriocin is considered present in a species if that bacteriocin is present at least in one more strain of that species. KEGG Pathway Analysis The proteomes of all the organisms were subject to bidirectional best hits to KEGG database [64, 96] to identify the components of folate biosynthesis pathways, streptomycin biosynthesis, and vancomycin resistance pathways. List of Abbreviations AA: amino-acid; bp: base pairs; CDS: Coding DNA Sequences; COGs: Clusters of Orthologous Groups; GGDC: Genome-Genome Distance Calculator; GIT: Gastrointestinal Tract; KEGG: Kyoto Encyclopedia of Genes and Genomes; M: modification; NT: nucleotide database; PABA: para-aminobenzoic acid; R: restriction; RRDR: RIF resistance determining region; S: specificity Declarations Ethics approval and consent to participate Not applicable. Consent for publication Not Applicable Availability of supporting data This Whole Genome Shotgun project for chromosome and plasmid of Bacillus clausii ENTPro has been deposited at DDBJ/EMBL/GenBank under the accession CP012475 and CP012476. Competing interests All the authors declare that they have no conflict of interests. Funding This work was supported by the Council of Scientific Research (CSIR) Network projects on Man as a super-organism: Understanding the Human Microbiome (HUM-CSIR-BSC-0119) and the Department of Biotechnology project BTISNET (GAP001). IK and GS were supported by research fellowships from the University Grants Commission (UGC) and Council of Scientific and Industrial Research (CSIR) and respectively. The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript. Authors’ Contributions SS conceived the idea; GS isolated genomic DNA carried out strain identification and performed genome assembly and annotation. IK and GS carried out the comparative analysis. IK, GS, and SS wrote the manuscript. All authors have read and approved the manuscript. Acknowledgments We thank the next-generation genomics facility at C-CAMP ( http://www.ccamp.res.in/ ), Bangalore and Genome Quebec Centre, McGill University, Canada for help in obtaining NGS data. References Fuller R. Probiotics in man and animals. J Appl Bacteriol. 1989;66:365–78. http://www.ncbi.nlm.nih.gov/pubmed/2666378. Cutting SM. Bacillus probiotics. Food Microbiol. 2011;28:214–20. doi:10.1016/j.fm.2010.03.007. Nicholson WL, Munakata N, Horneck G, Melosh HJ, Setlow P. Resistance of Bacillus endospores to extreme terrestrial and extraterrestrial environments. Microbiol Mol Biol Rev. 2000;64:548–72. http://www.ncbi.nlm.nih.gov/pubmed/10974126. Casula G, Cutting SM. Bacillus probiotics: spore germination in the gastrointestinal tract. Appl Env Microbiol. 2002;68:2344–52. http://www.ncbi.nlm.nih.gov/pubmed/11976107. Kim SG, Becattini S, Moody TU, Shliaha P V., Littmann ER, Seok R, et al. Microbiota-derived lantibiotic restores resistance against vancomycin-resistant Enterococcus. Nature. 2019;572:665–9. doi:10.1038/s41586-019-1501-z. Senok AC, Ismaeel AY, Botta GA. Probiotics: facts and myths. Clin Microbiol Infect. 2005;11:958–66. doi:10.1111/j.1469-0691.2005.01228.x. Porwal S, Lal S, Cheema S, Kalia VC. Phylogeny in aid of the present and novel microbial lineages: diversity in Bacillus . PLoS One. 2009;4:e4438. doi:10.1371/journal.pone.0004438. Senesi S, Celandroni F, Tavanti A, Ghelardi E. Molecular characterization and identification of Bacillus clausii Strains marketed for use in oral bacteriotherapy. Appl Env Microbiol. 2001;67:834–9. doi:10.1128/AEM.67.2.834-839.2001. Ciffo F. Determination of the spectrum of antibiotic resistance of the “Bacillus subtilis” strains of Enterogermina. Chemioterapia. 1984;3:45–52. Mazza P, Zani F, Martelli P. Studies on the antibiotic resistance of Bacillus subtilis strains used in oral bacteriotherapy. Boll Chim Farm. 1992;131:401–8. http://www.ncbi.nlm.nih.gov/pubmed/1299263. Green DH, Wakeley PR, Page A, Barnes A, Baccigalupi L, Ricca E, et al. Characterization of two Bacillus probiotics. Appl Environ Microbiol. 1999;65:4288–91. http://www.ncbi.nlm.nih.gov/pubmed/10473456. Accessed 29 Apr 2018. Lippolis R, Siciliano RA, Mazzeo MF, Abbrescia A, Gnoni A, Sardanelli AM, et al. Comparative secretome analysis of four isogenic Bacillus clausii probiotic strains. Proteome Sci. 2013;11:28. doi:10.1186/1477-5956-11-28. Ripert G, Racedo SM, Elie A-M, Jacquot C, Bressollier P, Urdaci MC. Secreted Compounds of the Probiotic Bacillus clausii Strain O/C Inhibit the Cytotoxic Effects Induced by Clostridium difficile and Bacillus cereus Toxins. Antimicrob Agents Chemother. 2016;60:3445–54. doi:10.1128/AAC.02815-15. Varankovich N V., Nickerson MT, Korber DR. Probiotic-based strategies for therapeutic and prophylactic use against multiple gastrointestinal diseases. Front Microbiol. 2015;6:685. doi:10.3389/fmicb.2015.00685. Hammad AM, Shimamoto T. Towards a compatible probiotic-antibiotic combination therapy: assessment of antimicrobial resistance in the Japanese probiotics. J Appl Microbiol. 2010;109:1349–60. doi:10.1111/j.1365-2672.2010.04762.x. Lopetuso LR, Scaldaferri F, Franceschi F, Gasbarrini A. Bacillus clausii and gut homeostasis: state of the art and future perspectives. Expert Rev Gastroenterol Hepatol. 2016;:1–6. doi:10.1080/17474124.2016.1200465. Pradhan B, Guha D, Ray P, Das D, Aich P. Comparative Analysis of the Effects of Two Probiotic Bacterial Strains on Metabolism and Innate Immunity in the RAW 264.7 Murine Macrophage Cell Line. Probiotics Antimicrob Proteins. 2016;8:73–84. doi:10.1007/s12602-016-9211-4. Patrone V, Molinari P, Morelli L. Microbiological and molecular characterization of commercially available probiotics containing Bacillus clausii from India and Pakistan. Int J Food Microbiol. 2016;237:92–7. doi:10.1016/j.ijfoodmicro.2016.08.012. Bozdogan B, Galopin S, Leclercq R. Characterization of a new erm-related macrolide resistance gene present in probiotic strains of Bacillus clausii . Appl Env Microbiol. 2004;70:280–4. http://www.ncbi.nlm.nih.gov/pubmed/14711653. Marseglia GL, Tosca M, Cirillo I, Licari A, Leone M, Marseglia A, et al. Efficacy of Bacillus clausii spores in the prevention of recurrent respiratory infections in children: a pilot study. Ther Clin Risk Manag. 2007;3:13–7. doi:10.2147/TCRM.2007.3.1.13. Gian Luigi Marseglia, Mariangela Tosca, Ignazio Cirillo, Amelia Licari, Maddalena Leone, Alessia Marseglia, Anna Maria Castellazzi and GC. Efficacy of Bacillus clausii spores in the prevention of recurrent respiratory infections in children: a pilot study. Ther Clin Risk Manag. 2007;3:13–7. https://www.dovepress.com/efficacy-of-bacillus-clausii-spores-in-the-prevention-of-recurrent-res-peer-reviewed-article-TCRM. Accessed 28 Apr 2018. Lakshmi SG, Jayanthi N, Saravanan M, Ratna MS. Safety assesment of Bacillus clausii UBBC07, a spore forming probiotic. Toxicol Reports. 2017;4:62–71. doi:10.1016/J.TOXREP.2016.12.004. Canani RB, Cirillo P, Terrin G, Cesarano L, Spagnuolo MI, De Vincenzo A, et al. Probiotics for treatment of acute diarrhoea in children: randomised clinical trial of five different preparations. BMJ. 2007;335:340. doi:10.1136/bmj.39272.581736.55. Větrovský T, Baldrian P. The Variability of the 16S rRNA Gene in Bacterial Genomes and Its Consequences for Bacterial Community Analyses. PLoS One. 2013;8:e57923. doi:10.1371/journal.pone.0057923. Auch AF, von Jan M, Klenk H-P, Göker M. Digital DNA-DNA hybridization for microbial species delineation by means of genome-to-genome sequence comparison. Stand Genomic Sci. 2010;2:117–34. doi:10.4056/sigs.531120. Saising J, Dube L, Ziebandt A-K, Voravuthikunchai SP, Nega M, Götz F. Activity of gallidermin on Staphylococcus aureus and Staphylococcus epidermidis biofilms. Antimicrob Agents Chemother. 2012;56:5804–10. doi:10.1128/AAC.01296-12. Goldstein BP. Resistance to rifampicin: a review. J Antibiot. 2014;67:625–30. doi:10.1038/ja.2014.107. Schwarz S, Kehrenberg C, It Doublet B, Cloeckaert A. Molecular basis of bacterial resistance to chloramphenicol and florfenicol. 2004. doi:10.1016/j.femsre.2004.04.001. Floch MH, Walker WA. Probiotics are considered nutritional supplements. J Clin Gastroenterol. 2011;45 Suppl:S107. doi:10.1097/MCG.0b013e3182317e74. Navarre WW, Schneewind O. Surface proteins of gram-positive bacteria and mechanisms of their targeting to the cell wall envelope. Microbiol Mol Biol Rev. 1999;63:174–229. http://www.ncbi.nlm.nih.gov/pubmed/10066836. Accessed 20 Jul 2018. Cotter PD, Hill C. Surviving the acid test: responses of gram-positive bacteria to low pH. Microbiol Mol Biol Rev. 2003;67:429–53, table of contents. doi:10.1128/MMBR.67.3.429-453.2003. Begley M, Gahan CG, Hill C. The interaction between bacteria and bile. FEMS Microbiol Rev. 2005;29:625–51. doi:10.1016/j.femsre.2004.09.003. Azcarate-Peril MA, Altermann E, Hoover-Fitzula RL, Cano RJ, Klaenhammer TR. Identification and inactivation of genetic loci involved with Lactobacillus acidophilus acid tolerance. Appl Environ Microbiol. 2004;70:5315–22. doi:10.1128/AEM.70.9.5315-5322.2004. Hagenbuch B, Dawson P. The sodium bile salt cotransport family SLC10. Pflugers Arch Eur J Physiol. 2004;447:566–70. doi:10.1007/s00424-003-1130-z. Wong MH, Oelkers P, Craddock AL, Dawson PA. Expression cloning and characterization of the hamster ileal sodium-dependent bile acid transporter. J Biol Chem. 1994;269:1340–7. http://www.ncbi.nlm.nih.gov/pubmed/8288599. Accessed 20 Jul 2018. Nachin L, Nannmark U, Nystrom T. Differential Roles of the Universal Stress Proteins of Escherichia coli in Oxidative Stress Resistance, Adhesion, and Motility. J Bacteriol. 2005;187:6265–72. doi:10.1128/JB.187.18.6265-6272.2005. Seifart Gomes C, Izar B, Pazan F, Mohamed W, Mraheil MA, Mukherjee K, et al. Universal stress proteins are important for oxidative and acid stress resistance and growth of Listeria monocytogenes EGD-e in vitro and in vivo. PLoS One. 2011;6:e24965. doi:10.1371/journal.pone.0024965. Ventura M, Canchaya C, Zink R, Fitzgerald GF, van Sinderen D. Characterization of the groEL and groES loci in Bifidobacterium breve UCC 2003: genetic, transcriptional, and phylogenetic analyses. Appl Environ Microbiol. 2004;70:6197–209. doi:10.1128/AEM.70.10.6197-6209.2004. Susin MF, Baldini RL, Gueiros-Filho F, Gomes SL. GroES/GroEL and DnaK/DnaJ have distinct roles in stress responses and during cell cycle progression in Caulobacter crescentus . J Bacteriol. 2006;188:8044–53. doi:10.1128/JB.00824-06. Veinger L, Diamant S, Buchner J, Goloubinoff P. The small heat-shock protein IbpB from Escherichia coli stabilizes stress-denatured proteins for subsequent refolding by a multichaperone network. J Biol Chem. 1998;273:11032–7. http://www.ncbi.nlm.nih.gov/pubmed/9556585. Accessed 20 Jul 2018. Narberhaus F. Alpha-crystallin-type heat shock proteins: socializing minichaperones in the context of a multichaperone network. Microbiol Mol Biol Rev. 2002;66:64–93; table of contents. http://www.ncbi.nlm.nih.gov/pubmed/11875128. Accessed 20 Jul 2018. Jakob U, Gaestel M, Engel K, Buchner J. Small heat shock proteins are molecular chaperones. J Biol Chem. 1993;268:1517–20. http://www.ncbi.nlm.nih.gov/pubmed/8093612. Accessed 20 Jul 2018. Fu X, Adams Z, Liu R, Hepowit NL, Wu Y, Bowmann CF, et al. Methionine Sulfoxide Reductase A (MsrA) and Its Function in Ubiquitin-Like Protein Modification in Archaea . MBio. 2017;8. doi:10.1128/mBio.01169-17. Marshall B, Petrowski D, Levy SB, Summers AO. Inter- and intraspecies spread of Escherichia coli in a farm environment in the absence of antibiotic usage. Proc Natl Acad Sci U S A. 1990;87:6609–13. doi:10.1073/pnas.87.17.6609. Gibson MK, Forsberg KJ, Dantas G. Improved annotation of antibiotic resistance determinants reveals microbial resistomes cluster by ecology. ISME J. 2015;9:207–16. doi:10.1038/ismej.2014.106. Gueimonde M, Sánchez B, G de Los Reyes-Gavilán C, Margolles A. Antibiotic resistance in probiotic bacteria. Front Microbiol. 2013;4:202. doi:10.3389/fmicb.2013.00202. Martinez JL, Baquero F. Mutation frequencies and antibiotic resistance. Antimicrob Agents Chemother. 2000;44:1771–7. http://www.ncbi.nlm.nih.gov/pubmed/10858329. Accessed 29 Apr 2018. Cotter PD, Hill C, Ross RP. What’s in a name? Class distinction for bacteriocins - Author reply. Nat Rev Microbiol. 2006;4. doi:DOI 10.1038/nrmicro1273-c2. Heng NCK, Tagg JR. What’s in a name? Class distinction for bacteriocins. Nat Rev Microbiol. 2006;4. doi:DOI 10.1038/nrmicro1273-c1. Field D, Cotter PD, Hill C, Ross RP. Bioengineering Lantibiotics for Therapeutic Success. Front Microbiol. 2015;6:1363. doi:10.3389/fmicb.2015.01363. McAuliffe O, Ryan MP, Ross RP, Hill C, Breeuwer P, Abee T. Lacticin 3147, a broad-spectrum bacteriocin which selectively dissipates the membrane potential. Appl Environ Microbiol. 1998;64:439–45. http://www.ncbi.nlm.nih.gov/pubmed/9464377. Accessed 29 Apr 2018. Masuda Y, Ono H, Kitagawa H, Ito H, Mu F, Sawa N, et al. Identification and characterization of leucocyclicin Q, a novel cyclic bacteriocin produced by Leuconostoc mesenteroides TK41401. Appl Environ Microbiol. 2011;77:8164–70. doi:10.1128/AEM.06348-11. Shelburne CE, An FY, Dholpe V, Ramamoorthy A, Lopatin DE, Lantz MS. The spectrum of antimicrobial activity of the bacteriocin subtilosin A. J Antimicrob Chemother. 2006;59:297–300. doi:10.1093/jac/dkl495. Liu Z, Budiharjo A, Wang P, Shi H, Fang J, Borriss R, et al. The highly modified microcin peptide plantazolicin is associated with nematicidal activity of Bacillus amyloliquefaciens FZB42. Appl Microbiol Biotechnol. 2013;97:10081–90. doi:10.1007/s00253-013-5247-5. Kawai Y, Kemperman R, Kok J, Saito T. The circular bacteriocins gassericin A and circularin A. Curr Protein Pept Sci. 2004;5:393–8. http://www.ncbi.nlm.nih.gov/pubmed/15544534. Accessed 29 Apr 2018. Kemperman R, Kuipers A, Karsens H, Nauta A, Kuipers O, Kok J. Identification and characterization of two novel clostridial bacteriocins, circularin A and closticin 574. Appl Environ Microbiol. 2003;69:1589–97. doi:10.1128/AEM.69.3.1589-1597.2003. Begley M, Cotter PD, Hill C, Ross RP. Identification of a novel two-peptide lantibiotic, lichenicidin, following rational genome mining for LanM proteins. Appl Environ Microbiol. 2009;75:5451–60. doi:10.1128/AEM.00730-09. Simmonds RS, Pearson L, Kennedy RC, Tagg JR. Mode of action of a lysostaphin-like bacteriolytic agent produced by Streptococcus zooepidemicus 4881. Appl Environ Microbiol. 1996;62:4536–41. http://www.ncbi.nlm.nih.gov/pubmed/8953725. Accessed 29 Apr 2018. Liu W-T, Yang Y-L, Xu Y, Lamsa A, Haste NM, Yang JY, et al. Imaging mass spectrometry of intraspecies metabolic exchange revealed the cannibalistic factors of Bacillus subtilis. Proc Natl Acad Sci U S A. 2010;107:16286–90. doi:10.1073/pnas.1008368107. González-Pastor JE, Hobbs EC, Losick R. Cannibalism by sporulating bacteria. Science. 2003;301:510–3. doi:10.1126/science.1086462. Lin D, Qu LJ, Gu H, Chen Z. A 3.1-kb genomic fragment of Bacillus subtilis encodes the protein inhibiting growth of Xanthomonas oryzae pv. oryzae. J Appl Microbiol. 2001;91:1044–50. http://www.ncbi.nlm.nih.gov/pubmed/11851812. Accessed 29 Apr 2018. Rossi M, Amaretti A, Raimondi S. Folate Production by Probiotic Bacteria. Nutrients. 2011;3:118–34. doi:10.3390/nu3010118. Zhu T, Pan Z, Domagalski N, Koepsel R, Ataai MM, Domach MM. Engineering of Bacillus subtilis for enhanced total synthesis of folic acid. Appl Environ Microbiol. 2005;71:7122–9. doi:10.1128/AEM.71.11.7122-7129.2005. Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 1999;27:29–34. http://www.ncbi.nlm.nih.gov/pubmed/9847135. Accessed 30 Jul 2017. de Crécy-Lagard V. Variations in metabolic pathways create challenges for automated metabolic reconstructions: Examples from the tetrahydrofolate synthesis pathway. Comput Struct Biotechnol J. 2014;10:41–50. doi:10.1016/j.csbj.2014.05.008. Salem AR, Pattison JR, Foster MA. Folic acid and the methylation of homocysteine by Bacillus subtilis . Biochem J. 1972;126:993–1004. http://www.ncbi.nlm.nih.gov/pubmed/4627401. Accessed 20 Jul 2018. Khatri I, Sharma S, Ramya TNC, Subramanian S. Complete Genomes of Bacillus coagulans S-lac and Bacillus subtilis TO-A JPC, Two Phylogenetically Distinct Probiotics. PLoS One. 2016;11:e0156745. doi:10.1371/journal.pone.0156745. Casanovas-Massana A, Sala-Comorera L, Blanch AR. Quantification of tetracycline and chloramphenicol resistance in digestive tracts of bulls and piglets fed with Toyocerin®, a feed additive containing Bacillus toyonensis spores. Vet Microbiol. 2014;173:59–65. doi:10.1016/J.VETMIC.2014.07.005. Zhang S, Hu Y, Fan Q, Wang X, He J. Two-component system YvqEC-dependent bacterial resistance against vancomycin in Bacillus thuringiensis . Antonie Van Leeuwenhoek. 2015;108:365–76. doi:10.1007/s10482-015-0489-0. Abriouel H, Franz CMAP, Omar N Ben, Gálvez A. Diversity and applications of Bacillus bacteriocins. FEMS Microbiol Rev. 2011;35:201–32. doi:10.1111/j.1574-6976.2010.00244.x. Chin CS, Alexander DH, Marks P, Klammer AA, Drake J, Heiner C, et al. Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat Methods. 2013;10:563–9. doi:10.1038/nmeth.2474. Urdaci MC, Bressollier P, Pinchuk I. Bacillus clausii Probiotic Strains. J Clin Gastroenterol. 2004;38 Supplement 2:S86–90. doi:10.1097/01.mcg.0000128925.06662.69. Ghelardi E, Celandroni F, Salvetti S, Gueye SA, Lupetti A, Senesi S. Survival and persistence of Bacillus clausii in the human gastrointestinal tract following oral administration as spore-based probiotic formulation. J Appl Microbiol. 2015;119:552–9. doi:10.1111/jam.12848. Vecchione A, Celandroni F, Mazzantini D, Senesi S, Lupetti A, Ghelardi E. Compositional Quality and Potential Gastrointestinal Behavior of Probiotic Products Commercialized in Italy. Front Med. 2018;5:59. doi:10.3389/fmed.2018.00059. Gargis AS, O’Rourke A-LD, Sloan GL, Simmonds RS. Prevalence and acquisition of the genes for zoocin A and zoocin A resistance in Streptococcus equi subsp. zooepidemicus . J Mol Evol. 2009;68:498–505. doi:10.1007/s00239-009-9221-x. Williams LD, Burdock GA, Jiménez G, Castillo M. Literature review on the safety of Toyocerin®, a non-toxigenic and non-pathogenic Bacillus cereus var. toyoi preparation. Regul Toxicol Pharmacol. 2009;55:236–46. doi:10.1016/J.YRTPH.2009.07.009. Fuchs CS, Willett WC, Colditz GA, Hunter DJ, Stampfer MJ, Speizer FE, et al. The influence of folate and multivitamin use on the familial risk of colon cancer in women. Cancer Epidemiol Biomarkers Prev. 2002;11:227–34. http://www.ncbi.nlm.nih.gov/pubmed/11895870. Accessed 19 Jul 2018. White E, Shannon JS, Patterson RE, Hunter DJ, Stampfer MJ, Speizer FE, et al. Relationship between vitamin and calcium supplement use and colon cancer. Cancer Epidemiol Biomarkers Prev. 1997;6:769–74. http://www.ncbi.nlm.nih.gov/pubmed/9332757. Accessed 19 Jul 2018. Schallmey M, Singh A, Ward OP. Developments in the use of Bacillus species for industrial production. Can J Microbiol. 2004;50:1–17. doi:10.1139/w03-076. Kellner R, JUNG G, HORNER T, ZAHNER H, SCHNELL N, ENTIAN K-D, et al. Gallidermin: a new lanthionine-containing polypeptide antibiotic. Eur J Biochem. 1988;177:53–9. doi:10.1111/j.1432-1033.1988.tb14344.x. Kalyon B, Helaly SE, Scholz R, Nachtigall J, Vater J, Borriss R, et al. Plantazolicin A and B: Structure Elucidation of Ribosomally Synthesized Thiazole/Oxazole Peptides from Bacillus amyloliquefaciens FZB42. Org Lett. 2011;13:2996–9. doi:10.1021/ol200809m. Alvarez-Ordóñez A, Begley M, Clifford T, Deasy T, Considine K, O’Connor P, et al. Investigation of the Antimicrobial Activity of Bacillus licheniformis Strains Isolated from Retail Powdered Infant Milk Formulae. Probiotics Antimicrob Proteins. 2014;6:32–40. doi:10.1007/s12602-013-9151-1. Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, et al. The RAST Server: Rapid Annotations using Subsystems Technology. BMC Genomics. 2008;9:75. doi:10.1186/1471-2164-9-75. Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, et al. The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). Nucleic Acids Res. 2014;42 Database issue:D206-14. doi:10.1093/nar/gkt1226. Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25:955–64. http://www.ncbi.nlm.nih.gov/pubmed/9023104. Accessed 30 Jul 2017. Lagesen K, Hallin P, Rodland EA, Staerfeldt HH, Rognes T, Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007;35:3100–8. doi:10.1093/nar/gkm160. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ, Yu Y, et al. Basic local alignment search tool. J Mol Biol. 1990;215:403–10. doi:10.1016/S0022-2836(05)80360-2. Antipov D, Hartwick N, Shen M, Raiko M, Lapidus A, Pevzner PA. plasmidSPAdes: assembling plasmids from whole genome sequencing data. Bioinformatics. 2016;32:btw493. doi:10.1093/bioinformatics/btw493. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7. doi:10.1093/nar/gkh340. Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312–3. doi:10.1093/bioinformatics/btu033. Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin E V, et al. The COG database: an updated version includes eukaryotes. BMC Bioinformatics. 2003;4:41. doi:10.1186/1471-2105-4-41. Eddy SR. Hidden Markov models. Curr Opin Struct Biol. 1996;6:361–5. http://www.ncbi.nlm.nih.gov/pubmed/8804822. Eddy SR. Accelerated Profile HMM Searches. PLoS Comput Biol. 2011;7:e1002195. doi:10.1371/journal.pcbi.1002195. McArthur AG, Waglechner N, Nizam F, Yan A, Azad MA, Baylay AJ, et al. The comprehensive antibiotic resistance database. Antimicrob Agents Chemother. 2013;57:3348–57. doi:10.1128/AAC.00419-13. van Heel AJ, de Jong A, Montalbán-López M, Kok J, Kuipers OP. BAGEL3: Automated identification of genes encoding bacteriocins and (non-)bactericidal posttranslationally modified peptides. Nucleic Acids Res. 2013;41 Web Server issue:W448-53. doi:10.1093/nar/gkt391. Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M. The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004;32:277D – 280. doi:10.1093/nar/gkh063. Tables Table 1: Genome assembly and annotation statistics for Bacillus clausii ENTPro composite genome. Chromosome genome assembly and annotation statistics of Bacillus clausii ENTPro Chromosome Plasmid Sequencing data P6 polymerase and C4 [P6C4] Chemistry based PacBio sequencing Bio Project Number PRJNA242453 NCBI Accession number CP012475 CP012476 Genome size (in bp) 4,264,866 31,475 GC content (%) 44.75 39.9 Chromosome/Contig 1 1 CDS 4,384 40 % Coding sequences 86.73 84.45 CDS from (+) strand 2,254 35 CDS from (-) strand 2,130 5 Max. CDS length 9,509 2,711 Mean CDS length 843 664 Hypothetical proteins 1,215 33 Hypothetical proteins (%) 27.72 82.5 tRNA 76 NA rRNA 7 operons (21 rRNAs) NA Supplementary Files SuppFiguresClausii.pdf SupplementaryTablesBMCmicroGK.xlsx Cite Share Download PDF Status: Published Journal Publication published 30 Dec, 2019 Read the published version in BMC Microbiology → Version 3 posted Submission checks completed at journal 11 Dec, 2019 Editorial decision: Accept 11 Dec, 2019 You are reading this latest preprint version Show more versions Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-6155","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research article","associatedPublications":[],"authors":[{"id":249481,"identity":"b4f06eba-5d28-457a-8132-688df73917ee","order_by":1,"name":"Indu Khatri","email":"","orcid":"https://orcid.org/0000-0002-7993-1953","institution":"Leiden University Medical Center","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Indu","middleName":"","lastName":"Khatri","suffix":""},{"id":249482,"identity":"5fe8deec-fa2a-431b-8b81-e3771294748f","order_by":2,"name":"Gaurav Sharma","email":"","orcid":"","institution":"UC Davis","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Gaurav","middleName":"","lastName":"Sharma","suffix":""},{"id":249483,"identity":"542c57d8-36f6-48ad-9a9d-2728e10effeb","order_by":3,"name":"Srikrishna Subramanian","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA5ElEQVRIie3RPQrCMBTA8RcCcWntJi0BvULFSRC9SksHJ1FnBwOFN/mxKl7EsZLBpeDqqHgBoZO4mOoiCI+ODvkvISG/JBAAm+0P84EpbsamKGdMfVazKqQjgFcmUO6FWH0TsmCn02K67w9XjdPh8tjD2PMypicEkfUY5TZPRigT3l7m0N1uItAbgjQdhtxFXhIhGUIYngG0Q5O0cHE+FDKpPUsyOGU0kQ5T0kUdGSL4+xaIaBIsGBpybJuHdYIl+qF/jhVJ/Lx2Mw+btVq7+Hp/YC/01loXFPk5Air+js1ms9moXoGUQcpb9mZTAAAAAElFTkSuQmCC","orcid":"https://orcid.org/0000-0002-3263-1048","institution":"","correspondingAuthor":true,"submittingAuthor":false,"prefix":"","firstName":"Srikrishna","middleName":"","lastName":"Subramanian","suffix":""}],"badges":[],"createdAt":"2019-09-26 14:40:12","currentVersionCode":3,"declarations":"","doi":"10.21203/rs.2.15490/v3","doiUrl":"https://doi.org/10.21203/rs.2.15490/v3","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1186/s12866-019-1680-7","type":"published","date":"2019-12-30T12:00:00+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":249588,"identity":"d6ba5766-9f60-440c-93ce-f00b960f5a77","added_by":"auto","created_at":"2019-12-11 21:24:20","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":458601,"visible":true,"origin":"","legend":"Circular representation of the B. clausii ENTPro composite genome. Here B. clausii ENTPro was taken as reference genome and circles from inside to outside represents as: Circle 1 and 2 represent GC content and GC skew of B. clausii ENTPro, circle 3 represent encoded RNAs in B. clausii ENTPro; circle 4 represent genes encoded by B. clausii ENTPro; circle 5 depicts B. clausii ENTPro composite chromosome; Circle 6 depicts the mapping of B. clausii ENTPro genome against genome of B. clausii KSM-K16; further the circles 7 to 13 represents the genomes of B. sp. JCM 19045, B. sp. JCM 19046, B. sp. JCM 19047, B. lehensis G1, B. halodurans C-125, B. cellulosilyticus DSM 2522, B. pseudofirmus OF4 mapped on B. clausii ENTPro respectively. BRIG 0.95 was used to build the circular representation. Mapping studies were done using BLASTn with an E-value cut-off 1e-5.","description":"","filename":"1.png","url":"https://assets-eu.researchsquare.com/files/771f7b58-e357-4d78-a9e8-1c7cc017ef1d/v3/1.png"},{"id":249591,"identity":"c86aca2e-26c8-4fc4-8bf0-86356c6da386","added_by":"auto","created_at":"2019-12-11 21:24:21","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":53794,"visible":true,"origin":"","legend":"Cluster of Orthologous groups [COG] categories in B. clausii genomes. The X-axis represents the COG groups and the Y-axis represents the average number of proteins in respective COG groups. The genomes of B. clausii are clustered as per the properties. B. clausii Heroin represents all the organisms isolated from the Heroin samples [PRJNA395369]. B clausii Probiotics represents all the probiotic strains as a single category. The COG categories are identified by capital letters as follows: A, RNA processing and modification; B, Chromatin structure and dynamics; C, energy production and conversion; D, cell cycle control, cell division and chromosome partitioning; E, amino acid transport and metabolism; F, nucleotide transport and metabolism; G, carbohydrate transport and metabolism; H, coenzyme transport and metabolism; I, lipid transport and metabolism; J, translation; K, transcription; L, replication; M, cell wall/membrane/envelope biogenesis; N, cell motility; O, posttranslational modification, protein turnover, chaperones; P, inorganic ion transport and metabolism; Q, secondary metabolites biosynthesis, transport and catabolism; R, general function prediction only; S, function unknown; T, signal transduction mechanisms; U, intracellular trafficking and secretion; V, defense mechanisms; and X, Mobilome: prophages, transposons.","description":"","filename":"2.png","url":"https://assets-eu.researchsquare.com/files/771f7b58-e357-4d78-a9e8-1c7cc017ef1d/v3/2.png"},{"id":249595,"identity":"3f304c9a-5e18-4f91-8a2a-92205b8207da","added_by":"auto","created_at":"2019-12-11 21:24:22","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":409214,"visible":true,"origin":"","legend":"Housekeeping proteins based Maximum Likelihood phylogenetic tree. \nThe outgroups are colored with a yellow background, Bacillus species are colored with light green background and B. clausii members are colored with dark green background. The Bacillus probiotics are written with red labels with double the size of the rest of the organisms. Values on branches represent the bootstrap values.","description":"","filename":"3.png","url":"https://assets-eu.researchsquare.com/files/771f7b58-e357-4d78-a9e8-1c7cc017ef1d/v3/3.png"},{"id":249597,"identity":"b0937385-8e2b-447f-9e0f-410045fb2768","added_by":"auto","created_at":"2019-12-11 21:24:22","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":92297,"visible":true,"origin":"","legend":"The binary matrix of Antibiotic resistance in Bacillus probiotics. Red marks the presence whereas black is absent. The rows represent the name of the Antibiotic resistance categories and the columns are the Bacillus probiotics. The genomes in Bacillus spp. are clustered as per the species they belong to. Multiple strains of same species are clustered together and are represented by the species rather than strains.","description":"","filename":"4.png","url":"https://assets-eu.researchsquare.com/files/771f7b58-e357-4d78-a9e8-1c7cc017ef1d/v3/4.png"},{"id":249599,"identity":"3918de88-5241-4ceb-9d9b-ef8528cec28c","added_by":"auto","created_at":"2019-12-11 21:24:23","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":74267,"visible":true,"origin":"","legend":"The binary matrix of Bacteriocins in Bacillus probiotics. Red marks the presence of the bacteriocin whereas black is absent. The rows represent the name of the bacteriocins, and the columns are the Bacillus probiotics. The genomes in Bacillus spp. are clustered as per the species they belong to. Multiple strains of same species are clustered together and are represented by the species rather than strains.","description":"","filename":"5.png","url":"https://assets-eu.researchsquare.com/files/771f7b58-e357-4d78-a9e8-1c7cc017ef1d/v3/5.png"},{"id":249601,"identity":"170b6bea-e74f-45ab-9ba4-a5086fbf2920","added_by":"auto","created_at":"2019-12-11 21:24:23","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":97208,"visible":true,"origin":"","legend":"The binary matrix of components of folate biosynthesis pathways in Bacillus probiotics. Sea green marks the presence whereas yellow is absent. The rows represent the name of the Bacillus probiotics whereas column represents the pathway components.","description":"","filename":"6.png","url":"https://assets-eu.researchsquare.com/files/771f7b58-e357-4d78-a9e8-1c7cc017ef1d/v3/6.png"},{"id":13481393,"identity":"b5655b9f-8cc4-49a5-8cae-025ed17db75b","added_by":"auto","created_at":"2021-09-16 21:46:57","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":2074566,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-6155/v3/17862440-a918-4645-95cd-0c80c0641d5d.pdf"},{"id":249592,"identity":"fb132279-bdae-4ec5-93c2-dcbd0d896ab6","added_by":"auto","created_at":"2019-12-11 21:24:21","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"supplement","size":4683444,"visible":true,"origin":"","legend":"","description":"","filename":"SuppFiguresClausii.pdf","url":"https://assets-eu.researchsquare.com/files/771f7b58-e357-4d78-a9e8-1c7cc017ef1d/v3/Supp-Figures-Clausii.pdf"},{"id":249589,"identity":"18f6816e-28d1-470e-ba74-2d7f7db7a262","added_by":"auto","created_at":"2019-12-11 21:24:20","extension":"xlsx","order_by":0,"title":"","display":"","copyAsset":false,"role":"supplement","size":23085,"visible":true,"origin":"","legend":"","description":"","filename":"SupplementaryTablesBMCmicroGK.xlsx","url":"https://assets-eu.researchsquare.com/files/771f7b58-e357-4d78-a9e8-1c7cc017ef1d/v3/Supplementary Tables BMCmicro-GK.xlsx"}],"financialInterests":"","formattedTitle":"Composite genome sequence of Bacillus clausii, a probiotic commercially available as Enterogermina®, and insights into its probiotic properties","fulltext":[{"header":"Background","content":"\u003cp\u003eProbiotics are live microbes which when consumed in sufficient amount helps to resume the original gut microflora, distressed by diarrhea or antibiotic intake [1]. Most bacterial probiotics such as \u003cem\u003eLactobacillus\u003c/em\u003e and \u003cem\u003eBifidobacteria,\u003c/em\u003e which are inhabitants of the gut, are available as lyophilized preparations of vegetative cells while some probiotic bacterial preparations that belong to the genus \u003cem\u003eBacillus\u003c/em\u003e are available in the form of spores [2]. Bacterial spores are dormant and resistant to heat, desiccation, dehydration and are extremely stable, which is a desirable property for probiotics [3, 4]. The spores of \u003cem\u003eBacillus\u003c/em\u003e germinates in the gut and the vegetative cells are vital for the human gut health [5]. \u003cem\u003eBacillus subtilis\u003c/em\u003e belongs to one of the most studied and explored family \u003cem\u003eBacillaceae and \u003c/em\u003emany of its strains are being used as probiotics since the 1990s. The role of \u003cem\u003eBacillus\u003c/em\u003e species ranges from probiotic nature of \u003cem\u003eB. subtilis\u003c/em\u003e, \u003cem\u003eB. clausii\u003c/em\u003e, \u003cem\u003eB. coagulans\u003c/em\u003e, \u003cem\u003eB. pumilus\u003c/em\u003e and other strains to biological control agents (\u003cem\u003eB. thuringiensis \u003c/em\u003eand \u003cem\u003eB. sphaericus\u003c/em\u003e), and pathogenicity (\u003cem\u003eB. anthracis\u003c/em\u003e and \u003cem\u003eB. cereus\u003c/em\u003e). Several strains are economically important (\u003cem\u003eB. subtilis\u003c/em\u003e) whereas others have medical importance (\u003cem\u003eB. licheniformis\u003c/em\u003e) [6]. Some \u003cem\u003eBacillus\u003c/em\u003e spp. are industrially-important and produce proteins such as alkaline proteases, xylanases, amylases, and cellulases [4, 7].\u003c/p\u003e\n\u003cp\u003eIn 2001, some \u003cem\u003eB. subtilis\u003c/em\u003e strains, which were used in the probiotics and soap industry, were reclassified as \u003cem\u003eB. clausii\u003c/em\u003e [8]. \u003cem\u003eB. clausii\u003c/em\u003e spores are marketed as the probiotic Enterogermina\u003cstrong\u003e\u0026reg; \u003c/strong\u003ewhich consists of four \u003cem\u003eBacillus\u003c/em\u003e strains O/C, N/R, SIN, and T that are resistant to Chloramphenicol, Novobiocin/Rifampicin, Neomycin/Streptomycin, and Tetracycline, respectively [9, 10]. Although these four strains are known to have been derived from a single penicillin-resistant strain, \u003cem\u003eB. subtilis\u003c/em\u003e ATCC 9799 [9, 11], secretome analysis have revealed variation in the expression level of some of the secreted proteins [12]. Also, the O/C strain of \u003cem\u003eB. clausii\u003c/em\u003e inhibits the cytotoxic effect induced by the \u003cem\u003eClostridium difficile\u003c/em\u003e and \u003cem\u003eB. cereus \u003c/em\u003etoxins [13]. The intrinsic antibiotic resistance in probiotics is considered advantageous in cases of antibiotics-probiotics combination prescriptions to restore healthy gut [14, 15]. The mode of action of \u003cem\u003eB. clausii\u003c/em\u003e as a probiotic is not clear, but the strains have been reported to secrete some proteins that are involved in the immunomodulatory mechanism, adaptation and their colonization in the human gastrointestinal tract (GIT) [12, 13, 16\u0026ndash;18]. Uniquely, \u003cem\u003eB. clausii\u003c/em\u003e harbor \u003cem\u003eerm(34)\u003c/em\u003e gene that imparts the resistance to erythromycin. The \u003cem\u003eerm(34)\u003c/em\u003e gene is not a homolog of \u003cem\u003eerm(A)\u003c/em\u003e, \u003cem\u003eerm(B)\u003c/em\u003e, \u003cem\u003eerm(C)\u003c/em\u003e and \u003cem\u003eerm(TR)\u003c/em\u003e genes for MLSB resistance in Gram-positive human pathogens and \u003cem\u003eerm(D)\u003c/em\u003e, \u003cem\u003eerm(K)\u003c/em\u003e, and \u003cem\u003eerm(J) \u003c/em\u003echaracterized in \u003cem\u003eB. licheniformis\u003c/em\u003e, \u003cem\u003eB. halodurans\u003c/em\u003e, and \u003cem\u003eB. anthracis,\u003c/em\u003e respectively [19].\u003c/p\u003e\n\u003cp\u003eVarious clinical trials and molecular studies [8, 13, 16, 20\u0026ndash;23] have been performed to identify the major features that demarcate \u003cem\u003eB. clausii\u003c/em\u003e probiotic strains from other \u003cem\u003eBacillus\u003c/em\u003e spp., but still, the genomic reasons of its probiotic activity have not been reported before. Therefore, we sequenced the composite genome sequence of \u003cem\u003eBacillus clausii\u003c/em\u003e (composite of all four strains of \u003cem\u003eB. clausii \u003c/em\u003eused in the probiotic formulation) from Enterogermina\u003cstrong\u003e\u0026reg;\u003c/strong\u003e, an oral probiotic, marketed by Sanofi in India. The composite genome obtained from the sequencing of this probiotics was named as \u003cem\u003eB. clausii \u003c/em\u003eENTPro. We performed an extensive analysis to identify the genomic features known to impart probiotic properties in \u003cem\u003eB. clausii\u003c/em\u003e \u003cem\u003eviz.,\u003c/em\u003e adhesion to gut, withstanding harsh conditions in the gut, antibiotic resistance, and biosynthesis pathways. In addition, to gain insight into the genomic features of different probiotics, we have compared pathways, types of bacteriocins and antibiotic resistance genes in different \u003cem\u003eBacillus\u003c/em\u003e probiotics.\u003c/p\u003e"},{"header":"Results","content":"\u003cp\u003e\u003cstrong\u003eGenome features of \u003cem\u003eB. clausii \u003c/em\u003eENTPro\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eDe novo\u003c/em\u003e assembly of PacBio sequencing reads of \u003cem\u003eB. clausii\u003c/em\u003e ENTPro gDNA resulted in two contigs: one long circular contig of 4,264,866 base pairs (bp) and one short circular 31,475 bp contig. The long contig represents the composite circular chromosome of \u003cem\u003eBacillus clausii\u003c/em\u003e ENTPro with an average GC content of 44.75% \u003cstrong\u003e(Fig. \u003c/strong\u003e\u003ca href=\"http://genomebiology.com/2004/5/10/r77/table/T1\"\u003e\u003cstrong\u003e1\u003c/strong\u003e\u003c/a\u003e\u003cstrong\u003e)\u003c/strong\u003e and the smaller one (GC Content: 39.9%) is likely a plasmid. In addition, Illumina sequencing-based assembly resulted in 4.3 Mbp genome from 36 contigs and N50 of 344,696 bp, which overlaps completely with the genome assembled using PacBio reads. The composite genome obtained from PacBio sequencing reads was submitted to GenBank [NC_006582.1] and further used for all the comparisons in this study. \u003cem\u003eBacillus clausii\u003c/em\u003e ENTPro genome is 99.8% similar to another probiotic strain \u003cem\u003eB. clausii\u003c/em\u003e B106 [NFZO01] (\u003cstrong\u003eFig. S1A\u003c/strong\u003e), followed by 94.3% similarity to \u003cem\u003eB. clausii KSM-K16 \u003c/em\u003e[NC_006582.1] (\u003cstrong\u003eFig. S1B\u003c/strong\u003e), whereas other members of the same species are 50-94% similar. This suggests that the members of this species are quite diverse as characterized by their GGDC values (\u003cstrong\u003eTable S1\u003c/strong\u003e). Our analysis suggests that probiotic strains within \u003cem\u003eB. clausii\u003c/em\u003e such as ENTPro, B106, and UBBC-07 are highly similar to each other as compared to other strains.\u003c/p\u003e\n\u003cp\u003eThe plasmid sequence is novel and does not have any close similarity with other plasmids in the NCBI nucleotide database (NT). Most of the proteins encoded by the plasmid sequence are hypothetical and are not functionally characterized. We mapped Illumina reads against the plasmid database downloaded from NCBI to identify if we could obtain hits to any previously known plasmids. Very few reads mapped on to known plasmids and no full plasmid could be retrieved using the Illumina reads. Therefore, we concluded that the identified plasmid sequence harbored by \u003cem\u003eB. clausii\u003c/em\u003e ENTPro is novel.\u003c/p\u003e\n\u003cp\u003eAnnotation of the \u003cem\u003eB. clausii\u003c/em\u003e ENTPro genome revealed the presence of 4,384 protein-coding sequences, which constitute 86.73% of the genome with an average length of 843 bp (ranging from 113 to 9,509 bp) \u003cstrong\u003e(Table 1)\u003c/strong\u003e. A total of 1,215 Coding DNA Sequences (CDS) were annotated as hypothetical proteins, accounting for 27.72% of the total proteins. The ENTPro genome has all the three proteins R (restriction), M (modification), and S (specificity) that belongs to the Type I RM system. m6A methylation was observed in \u0026gt;96% of the motifs G\u003cstrong\u003eA\u003c/strong\u003eGNNNNNNRTGC and GC\u003cstrong\u003eA\u003c/strong\u003eYNNNNNNCTC in the genome at 2\u003csup\u003end\u003c/sup\u003e and 3\u003csup\u003erd\u003c/sup\u003e positions, respectively. There are 75 tRNA genes and seven complete rRNA operons (\u0026gt;99% identity) in the \u003cem\u003eB. clausii\u003c/em\u003e ENTPro genome. 16S rRNAs obtained from the \u003cem\u003ede novo\u003c/em\u003e assembly of \u003cem\u003eB. clausii \u003c/em\u003eENTPro genome shows 99.8% similarity with \u003cem\u003eB. clausii\u003c/em\u003e Enterogermina strains O/C, T, N/R, and SIN. This is in line to previously known variations in 16S rRNA genes in bacterial genomes [24]. \u0026nbsp;Most of the varying sites were present in the V1 region of the 16S rRNA sequences even in \u003cem\u003eB. clausii\u003c/em\u003e KSM-K16 and \u003cem\u003eB. clausii\u003c/em\u003e DSM 8716\u003cstrong\u003e(Fig. S2)\u003c/strong\u003e.\u003c/p\u003e\n\u003cp\u003eAmongst the total proteome, ~75% (3,311) proteins could be categorized into Clusters of Orthologous Groups (COGs) functional groups. Among these mapped proteins, ~35% belonged to the metabolism category, ~14% to cellular processes and signaling and ~16% proteins to information storage and processing. According to COG mapping data, 152 proteins are involved in signal transduction mechanisms (COG: T) and 44 proteins were reported to function in secondary metabolites biosynthesis, transport, and catabolism (COG: Q). COG assignments to proteomes of \u003cem\u003eB. clausii\u003c/em\u003e members revealed that all the organisms have similar number of proteins assigned to various COG categories (\u003cstrong\u003eFig. 2\u003c/strong\u003e).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003ePhylogenetic position of \u003cem\u003eB. clausii\u003c/em\u003e as inferred from housekeeping proteins-based phylogeny\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003ePhylogenetically, \u003cem\u003eB. clausii\u003c/em\u003e clustered in a separate clade with further grouping within this clade (\u003cstrong\u003eFig. 3\u003c/strong\u003e). The phylogenetic tree reveals that ENTPro strain is closest to the B106 strain of \u003cem\u003eB. clausii\u003c/em\u003e. Both these probiotic strains are further similar to another probiotic strain UBBC-07 of \u003cem\u003eB. clausii\u003c/em\u003e. All these probiotic strains share a common ancestor with industrial \u003cem\u003eB. clausii\u003c/em\u003e KSM-K16 strain. This phylogenetic placement of \u003cem\u003eB. clausii\u003c/em\u003e probiotic strains is concordant with the whole genome similarity matrix as obtained by genome-genome distance calculator (GGDC) [25]. Other \u003cem\u003eB. clausii\u003c/em\u003e \u0026ldquo;Heroin\u0026rdquo; strains form several different groups within the \u003cem\u003eB. clausii\u003c/em\u003e clade. Interestingly, the \u003cem\u003eB. clausii\u003c/em\u003e proteome matches the proteome of other \u003cem\u003eBacillus \u003c/em\u003especies at \u0026lt;70% identity. This clearly suggests the genomic heterogeneity of \u003cem\u003eB. clausii\u003c/em\u003e genome in comparison to other \u003cem\u003eBacillus \u003c/em\u003especies. We also included all \u003cem\u003eBacillus\u003c/em\u003e probiotics genomes in phylogenetic analysis to investigate their position phylogenetically [26]. \u003cem\u003eBacillus\u003c/em\u003e probiotics shared clades with their species members. Interestingly, probiotic strains cluster together e.g. \u003cem\u003eB. clausii\u003c/em\u003e, \u003cem\u003eB. coagulans\u003c/em\u003e and \u003cem\u003eB. subtilis\u003c/em\u003e.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003e\u003cem\u003eB. clausii\u003c/em\u003e ENTPro as a derived strain from four different strains\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eB. clausii \u003c/em\u003eEnterogermina\u003cstrong\u003e\u0026reg;\u003c/strong\u003e is a mixture of four different strains each of which is supposed to confer resistance against specific antibiotics, namely novobiocin and rifampicin (strain N/R), chloramphenicol (strain O/C), streptomycin and neomycin (strain SIN) and tetracycline (strain T) [12]. The specific genes conferring resistance could not be traced in the literature so different \u003cem\u003ein silico\u003c/em\u003e strategies were employed to identify possible genes that could help impart resistance to these antibiotics in \u003cem\u003eB. clausii\u003c/em\u003e ENTPro \u003cstrong\u003e(\u003c/strong\u003e\u003cstrong\u003eTable S2 and S3)\u003c/strong\u003e.\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eRifampicin:\u003c/em\u003e Rifampicin resistance is acquired by specific mutations at positions 516, 526 and 531 in the \u003cem\u003erpoB\u003c/em\u003e gene in \u003cem\u003eEscherichia coli\u003c/em\u003e [27]. These mutations are mapped in the center of the \u003cem\u003erpoB\u003c/em\u003e gene in 3 regions: one cluster covering 507-533 amino-acid (AA); cluster II covering AA 563\u0026ndash;572 and cluster III with AA change at position 687, which altogether are referred to as RIF resistance determining region (RRDR) [27]. In order to find the presence of RRDR region in RpoB protein in ENTPro, the RpoB protein sequences from all \u003cem\u003eBacillus \u003c/em\u003espp. were retrieved and aligned with \u003cem\u003eE. coli\u003c/em\u003e RpoB protein sequence [Accession Number: NP_418414.1]. P\u003csub\u003e524\u003c/sub\u003e-\u0026gt;L (corresponding to 567 AA position in \u003cem\u003eE. coli\u003c/em\u003e RpoB protein sequence) AA change was observed in \u003cem\u003eB. clausii\u003c/em\u003e ENTPro strain that was not observed in other \u003cem\u003eBacillus \u003c/em\u003espp. (\u003cstrong\u003eFig. S3\u003c/strong\u003e).\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eChloramphenicol:\u003c/em\u003e Chloramphenicol acetyltransferase, involved in conferring resistance against chloramphenicol [28], was identified from the proteome of \u003cem\u003eB. clausii\u003c/em\u003e ENTPro [Accession Number: WP_035203840.1].\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eStreptomycin:\u003c/em\u003e Pfam domains, known to impart resistance against streptomycin, were identified in \u003cem\u003eB. clausii\u003c/em\u003e ENTPro. Nine proteins in \u003cem\u003eB. clausii\u003c/em\u003e had the Pfam domain PF02522, PF01636, PF01909, PF04439, PF04655, PF07091, PF07827, and PF10706 that has core domain aminoglycoside. Two proteins had streptomycin adenylyltransferase domain (PF04439), six proteins have aminoglycoside phosphotransferase [PF01909] domain and one protein has Kanamycin nucleotidyltransferase [PF07827] domain in their sequence (\u003cstrong\u003eTable S2\u003c/strong\u003e). This suggests the presence of domains that are involved in imparting resistance to streptomycin. In addition, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis of the organism reveal the presence of complete KEGG pathway for the streptomycin biosynthesis in the \u003cem\u003eB. clausii\u003c/em\u003e ENTPro (\u003cstrong\u003eFig. S4\u003c/strong\u003e).\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eTetracycline:\u003c/em\u003e The domains conferring tetracycline resistance [RF0133, RF0134, RF0135, and RF0127] were present in \u003cem\u003eB. clausii\u003c/em\u003e ENTPro (\u003cstrong\u003eFig. 4\u003c/strong\u003e). The presence of these genes in the composite genome of \u003cem\u003eB. clausii\u003c/em\u003e ENTPro was further confirmed by mapping the Illumina reads to these genes.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eProbiotic Properties in \u003cem\u003eB. clausii\u003c/em\u003e ENTPro\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eProbiotics are beneficial components of microbiota that modulates immunological, respiratory and gastrointestinal functions [29]. For imparting these functions probiotics adhere to the mucosal membrane to interact with the host and have acidic, alkaline and oxidative stress resistance and stress adaptation proteins [6]. Probiotics are believed to have good adherence capacity which promotes gut residence time, pathogen elimination and adhesion to the epithelial layer of host cells and exerting immune modulation.\u003c/p\u003e\n\u003cp\u003ePfam analysis reveals the presence of three proteins involved in adhesion namely a mucus-binding protein with \u0026lsquo;Gram_pos_anchor\u0026rsquo; Pfam domain [PF00746] at the C-terminus, a collagen-binding protein with LPXTG motif at the C-terminus and a fibronectin-binding protein [30] \u003cstrong\u003e(Table S2)\u003c/strong\u003e. These adhesion proteins may help facilitate the probiotic bacterium to bind and help in the direct interactions with the intestinal mucosa layer.\u003c/p\u003e\n\u003cp\u003eProbiotic \u003cem\u003eB. clausii \u003c/em\u003ehas to encounter various harsh environmental conditions during transit in the GIT such as the acidic environment in the stomach, bile juice environment in the small intestine, oxidative stress, and osmotic stress [1]. When a bacterium faces an acidic environment, H\u003csup\u003e+\u003c/sup\u003e homeostasis is maintained by F0F1 ATP synthase pump, which work by hydrolyzing ATP to pump protons (H\u003csup\u003e+\u003c/sup\u003e) from the cytoplasm [1, 31]. We found that this synthase complex is present in ENTPro genome as a full operon [DB29_02342--DB29_02349].\u003c/p\u003e\n\u003cp\u003eThe bacteria have to face the toxicity of bile salts that induce intracellular acidification and act as detergents that disrupt biological membranes [32]. Five proteins were identified that were involved in bile tolerance mechanism; two belong to ornithine decarboxylase [33] and three to sodium bile acid symporter family [34, 35] \u003cstrong\u003e(Table S2)\u003c/strong\u003e.\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eB. clausii\u003c/em\u003e ENTPro also harbors general stress adaptation proteins. The universal stress protein UspA [PF00582] is important for survival during cellular growth arrest and reprograms the cell towards defense and escape during cellular stress [36, 37]. Molecular chaperones that may impart resistance against environmental stress were obtained through annotation and Pfam domain search such as the chaperonin GroES [PF00166] and GroEL [38, 39] and one heat shock protein 33 [PF01430], two copies of cold shock proteins CspA [PF00313], three Clp protease [PF00574] and HtpX and HrcA-like heat shock proteins. These proteins play an important role in basic cellular functions that includes growth, the stability of DNA and RNA and they also prevent the formation of inclusion bodies [40\u0026ndash;42].\u003c/p\u003e\n\u003cp\u003eFor hyperosmotic stress and heat resistance, \u003cem\u003eB. clausii\u003c/em\u003e ENTPro harbors one copy each of the chaperone protein DnaJ [PF00226] and nucleotide exchange factor GrpE [PF01025]. Also, two methionine sulfoxide reductase A [43] [PF01625] were present in \u003cem\u003eB. clausii\u003c/em\u003e ENTPro that provides resistance in oxidative stress \u003cstrong\u003e(Table S2)\u003c/strong\u003e. This suggests that \u003cem\u003eB. clausii\u003c/em\u003e ENTPro has proteins to improve adhesion and handling stress and harsh conditions in the human gut.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAntibiotic Resistance in \u003cem\u003eBacillus\u003c/em\u003e Probiotics\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAntibiotic resistance is a common phenomenon in Gram-positive bacteria [44\u0026ndash;46]. It is accomplished by genes acquired either horizontally through plasmids, or foreign DNA recombination, or mutations at different chromosomal loci in the bacterial genome [47]. It is preferred that probiotic strains carry few antibiotic resistance genes as possible so that they are not a putative source for transferring these genes to other gut bacteria including pathogens [46]. However, on the other hand since some of these probiotics are administrated alongside antibiotics, some resistance to commonly administrated antibiotics are desirable.\u003c/p\u003e\n\u003cp\u003ePresence of a novel plasmid sequence in \u003cem\u003eB. clausii\u003c/em\u003e ENTPro could be a possible source of antibiotic-resistance gene transfer but we could not identify any potential antibiotic-resistance domain(s) in the plasmid. We also searched for the presence of antibiotic resistance genes and efflux pumps in the genomes with multiple methods to avoid false positives.\u003c/p\u003e\n\u003cp\u003eThe Chloramphenicol acetyltransferase, that confers resistance against chloramphenicol, is absent in \u003cem\u003eB. amyloliquefaciens\u003c/em\u003e and \u003cem\u003eB. coagulans\u003c/em\u003e whereas chloramphenicol efflux pump was present in \u003cem\u003eB. amyloliquefaciens \u003c/em\u003e(\u003cstrong\u003eFig. 4\u003c/strong\u003e). This would imply the presence of chloramphenicol resistance in all the \u003cem\u003eBacillus\u003c/em\u003e probiotics except \u003cem\u003eB. coagulans\u003c/em\u003e. Different classes of beta-lactamase were present in one or the other \u003cem\u003eBacillus\u003c/em\u003e probiotics that clearly suggest the presence of resistance against Penicillin in all the \u003cem\u003eBacillus\u003c/em\u003e probiotics. Multidrug resistance protein, a universal stress protein, EmrB, and its efflux pump, tetracycline resistance protein, and penicillin-binding protein are present in all the \u003cem\u003eBacillus\u003c/em\u003e probiotics. This suggests that most of the \u003cem\u003eBacillus\u003c/em\u003e probiotics are resistant to common antibiotics.\u003c/p\u003e\n\u003cp\u003eErythromycin resistance was identified by subjecting the \u003cem\u003eerm (34)\u003c/em\u003e gene sequence (GenBank Identifier: AY234334) of \u003cem\u003eB. clausii\u003c/em\u003e DSM8716 to BLASTn against all \u003cem\u003eBacillus\u003c/em\u003e genomes. This gene was identified in \u003cem\u003eB. clausii\u003c/em\u003e ENTPro named as \u0026ldquo;SSU rRNA (adenine(1518)-N(6)/adenine(1519)-N(6))-dimethyltransferase\u0026rdquo; (GenBank Identifier: ALA53582). The gene was also identified in all the \u003cem\u003eB. clausii \u003c/em\u003egenomes. The gene sequence shared 61% identity to rRNA adenine methyltransferase of \u003cem\u003eB. halodurans\u003c/em\u003e and 57% identity to rRNA adenine methyltransferase of \u003cem\u003eB. licheniformis\u003c/em\u003e, \u003cem\u003eB. anthracis\u003c/em\u003e, \u003cem\u003eB. sonorensis \u003c/em\u003eand\u003cem\u003e B. fordii\u003c/em\u003e. The rRNA adenine methyltransferase gene from other \u003cem\u003eBacillus \u003c/em\u003espp. shared 20-50% identity with \u003cem\u003eerm (34)\u003c/em\u003e gene. The result reveals that the \u003cem\u003eerm (34)\u003c/em\u003e gene is unique to \u003cem\u003eB. clausii\u003c/em\u003e and is not present in other members of the \u003cem\u003eBacillus\u003c/em\u003e genus.\u003c/p\u003e\n\u003cp\u003eVancomycin resistance, as observed from KEGG pathway analysis, (\u003cstrong\u003eFig. S5\u003c/strong\u003e) was identified only in \u003cem\u003eB. toyonensis \u003c/em\u003ewhile absent in other \u003cem\u003eBacillus\u003c/em\u003e probiotics. The accessory proteins of vancomycin resistance operon were present in some of the \u003cem\u003eBacillus\u003c/em\u003e probiotics, but resistance-conferring genes were completely absent.\u003c/p\u003e\n\u003cp\u003eWe would like to add an advisory note that previous studies have shown that an organism may exhibit intrinsic resistance to a few antibiotics that could not be related to its genotype [46]. Though we have endeavored to relate the genome-level occurrence of antibiotic resistance proteins or domains to their probable phenotypes, we have not performed any phenotypic studies to substantiate these analyses and/or confirm for intrinsic resistance. Further, the current situation may constitute a safety concern because of the possibility of transfer of antibiotic gene transfer to gut flora [46].\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eBacteriocins in \u003cem\u003eBacillus\u003c/em\u003e probiotics\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eBacteriocins are proteinaceous toxins produced by bacteria that act as narrow-spectrum antibiotics to inhibit the growth of similar or closely related bacterial strains [48, 49].\u0026nbsp; They can help probiotics to survive the toxins produced by invading bacteria by inhibiting their growth and hence can result in beneficial effects on the hosts. The identified bacteriocins in all the probiotics are represented in a presence-absence binary matrix in \u003cstrong\u003eFig. 5\u003c/strong\u003e. Several of these bacteriocins are already well utilized in therapeutics [50] and their spectrum against pathogens is well established [9, 51\u0026ndash;53].\u003c/p\u003e\n\u003cp\u003eGallidermin identified via \u003cem\u003ein silico\u003c/em\u003e analysis in \u003cem\u003eB. clausii\u003c/em\u003e genomes is known to efficiently prevent biofilm formation in the pathogens \u003cem\u003eS. aureus\u003c/em\u003e and \u003cem\u003eS. epidermidis\u003c/em\u003e species [26]. This bacteriocin has also been reported to be effective in skin disorders including acne, eczema, folliculitis, and impetigo where the targets organisms are Propionibacteria, Staphylococci, and Streptococci [50].\u003c/p\u003e\n\u003cp\u003eLacticin 3147 A2 and Leucocyclin Q as identified in \u003cem\u003eB. amyloliquefaciens\u003c/em\u003e are broad-spectrum bacteriocins. Lacticin has been used effectively in the treatment of bacterial mastitis, Staphylococcal and Enterococcal infections including vancomycin-resistant \u003cem\u003eEnterococci\u003c/em\u003e\u0026nbsp; [50] and is effective against \u003cem\u003eListeria\u003c/em\u003e infections [51]. Similarly, leucocyclicin Q exhibit bactericidal or bacteriostatic effects on Gram-positive bacteria, including food-borne pathogens, such as \u003cem\u003eLactococcus\u003c/em\u003e, \u003cem\u003eWeissella\u003c/em\u003e \u003cem\u003eparamesenteroides\u003c/em\u003e, \u003cem\u003ePediococcus\u003c/em\u003e \u003cem\u003edextrinicus\u003c/em\u003e, \u003cem\u003eEnterococcus\u003c/em\u003e, \u003cem\u003eStreptococcus\u003c/em\u003e, and \u003cem\u003eLeuconostoc\u003c/em\u003e [52]. Plantazolicin identified in \u003cem\u003eB. amyloliquefaciens\u003c/em\u003e and \u003cem\u003eB. pumilus\u003c/em\u003e has nematicidal activity [54]. Cirucularin A produced by \u003cem\u003eB. coagulans\u003c/em\u003e has been reported to be the most effective bacteriocin against \u003cem\u003eC.\u0026nbsp; tyrobutyricum\u003c/em\u003e NIZOB570,\u0026nbsp; a\u0026nbsp; known cheese-spoilage bacterium [55] and also Lactococci\u003cem\u003e,\u003c/em\u003e Enterococci, and some \u003cem\u003eLactobacillus\u003c/em\u003e strains [56]. LichenicidinVK21A2 identified in \u003cem\u003eB. paralicheniformis\u003c/em\u003e is considered as self-immunity bacteriocin that exhibits antimicrobial activity against several strains of \u003cem\u003eListeria monocytogenes\u003c/em\u003e, methicillin-resistant \u003cem\u003eS. aureus\u003c/em\u003e, and vancomycin-resistant Enterococcus [57]. Zoocin A in \u003cem\u003eB. toyonensis\u003c/em\u003e shows antimicrobial activity against several other Streptococci by cleaving the peptidoglycan cross-links of the target cell wall [58].\u003c/p\u003e\n\u003cp\u003eSubtilosin A produced by \u003cem\u003eB. subtilis\u003c/em\u003e is also a broad range bacteriocin that is effective against \u003cem\u003eListeria monocytogenes\u003c/em\u003e, and strains of \u003cem\u003eE. faecalis\u003c/em\u003e, \u003cem\u003eP. gingivalis\u003c/em\u003e, \u003cem\u003eK. rhizophila\u003c/em\u003e, \u003cem\u003eEnterobacter\u003c/em\u003e \u003cem\u003eaerogenes\u003c/em\u003e, \u003cem\u003eStreptococcus\u003c/em\u003e \u003cem\u003epyogenes,\u003c/em\u003e and \u003cem\u003eShigella\u003c/em\u003e \u003cem\u003esonnei\u003c/em\u003e [53]. Sporulation-killing factor skfA produced by \u003cem\u003eB. subtilis\u003c/em\u003e induces the lysis of other \u003cem\u003eB. subtilis\u003c/em\u003e cells that have not entered the sporulation pathway. This cannibalistic behavior provides a source of nutrients to support those cells that have entered sporulation [59, 60]. At high concentrations, it can also inhibit the growth of other bacteria [61]. The presence of well-characterized bacteriocins in the \u003cem\u003eBacillus\u003c/em\u003e probiotics suggests their important role in fighting against the pathogen in the gut.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFolate Biosynthesis Pathways in \u003cem\u003eBacillus\u003c/em\u003e Probiotics\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe gut microbiota aids the host, playing a crucial role in nutrient digestion and energy recovery. Due to potentially relevant applications, the capacity to yield folate has been investigated in various probiotic strains. Previously, the presence of these pathways was reported in \u003cem\u003eLactobacillus\u003c/em\u003e and \u003cem\u003eBifidobacterium\u003c/em\u003e probiotics but was not explored in \u003cem\u003eBacillus\u003c/em\u003e probiotics [62] except \u003cem\u003eB. subtilis\u003c/em\u003e [63]. We performed the identification of key components of folate production pathways in \u003cem\u003eBacillus\u003c/em\u003e probiotics using KEGG Pathway database [64]. The analysis of genome sequences of \u003cem\u003eBacillus\u003c/em\u003e probiotics revealed the presence of complete operon to synthesize para-aminobenzoic acid (PABA) \u003cem\u003ede novo\u003c/em\u003e only in \u003cem\u003eB. subtilis\u003c/em\u003e probiotics (\u003cstrong\u003eFig. 6).\u003c/strong\u003e On the other hand, the enzymes, necessary for chorismate conversion into PABA are present in almost all the \u003cem\u003eBacillus\u003c/em\u003e probiotics. Moreover, the shikimate pathway for chorismate production is complete only in \u003cem\u003eB. subtilis\u003c/em\u003e, \u003cem\u003eB. pumilus\u003c/em\u003e and \u003cem\u003eB. toyonensis\u003c/em\u003e, while it is partial in all the other \u003cem\u003eBacillus\u003c/em\u003e probiotics. On the other hand, \u003cem\u003eBacillus\u003c/em\u003e probiotic strains contain the genes of DHPPP \u003cem\u003ede novo\u003c/em\u003e biosynthetic pathway, the gene encoding dihydropteroate synthase (EC 2.5.1.15) and gene encoding dihydropteroate synthase (EC 2.5.1.15). Therefore, it is expected that these strains are not auxotrophic for folates or DHP but can produce folate in the presence of PABA supplementation. The presence/absence of the components of the folate biosynthesis pathway is reported based on KEGG pathway analysis. Previous studies have revealed that \u003cem\u003eB. subtilis\u003c/em\u003e genome harbor all the pathways components and have been engineered for folate production [63, 65, 66].\u003c/p\u003e"},{"header":"Discussion","content":"\u003cp\u003eIn this study, we report the complete composite genome of \u003cem\u003eBacillus clausii \u003c/em\u003eENTPro, sequenced from a commercially available probiotic, Enterogermina\u003cstrong\u003e\u0026reg; \u003c/strong\u003ethat is a mixture of four closely related strains. The oral probiotics are popularly recommended by physicians as an adjunct to antibiotic therapy to avoid antibiotic-induced diarrhea and/or other gastrointestinal distress. We determined the composite genome of the probiotic \u003cem\u003eBacillus clausii\u003c/em\u003e, to identify the features responsible for its probiotic properties and correlate it to the phenotypic properties mentioned previously in the literature [26, 50, 61, 67\u0026ndash;71]. While literature studies also mention the presence of four closely related strains in the probiotic Enterogermina\u003cstrong\u003e\u0026reg;, \u003c/strong\u003e[2, 12, 72] our single chromosome assembly suggests that the different strains are inseparable. However, variations in number of bacteria from each strain in the formulation has also been reported previously [73].\u003c/p\u003e\n\u003cp\u003eOur investigation revealed the presence of antibiotic resistance genes that \u003cem\u003eB. clausii\u003c/em\u003e ENTPro harbors. However, we could not find the genes conferring Novobiocin and Neomycin resistance, possibly because the genome reported here is a composite of the four strains. The possibility of the identified genome to be strain-specific could explain these missed antibiotic resistance genes. Previously reported secretome analysis [12] supports the presence of four related strains, and our genome data suggests that they are very closely related.\u003c/p\u003e\n\u003cp\u003eCOG assignments of all the \u003cem\u003eB. clausii\u003c/em\u003e strains did not show any differences in probiotic strains versus industrial strains of \u003cem\u003eB. clausii\u003c/em\u003e. Apart from these common features within the members of \u003cem\u003eB. clausii\u003c/em\u003e genus, the low genome-genome distance revealed that its members are quite diverse. We further used phylogenetic methods to understand the relationship of the members of \u003cem\u003eB. clausii\u003c/em\u003e and other \u003cem\u003eBacillus \u003c/em\u003eusing representatives of each species. The phylogenetic tree clearly has a separate group where members of \u003cem\u003eB. clausii\u003c/em\u003e lies with \u003cem\u003eB. lehensis\u003c/em\u003e, \u003cem\u003eB. halodurans\u003c/em\u003e whereas \u003cem\u003eB. subtilis, B. coagulans\u003c/em\u003e formed separate clades. The \u003cem\u003eB. clausii\u003c/em\u003e proteome analysis also supported its distant placement in the phylogenetic tree from other \u003cem\u003eBacillus\u003c/em\u003e members. These results indicate that \u003cem\u003eB. clausii \u003c/em\u003eis unique in comparison to other \u003cem\u003eBacillus\u003c/em\u003e species.\u003c/p\u003e\n\u003cp\u003eUsing \u003cem\u003ein silico\u003c/em\u003e mining approaches, we previously reported the sporulation proteins as well as various other proteins that might play a role in probiotic function such as molecular chaperones, stress proteins, flagellin, and mucin binding protein in two other \u003cem\u003eBacillus\u003c/em\u003e probiotics marketed in India namely \u003cem\u003eB. coagulans \u003c/em\u003eS-lac and \u003cem\u003eB. subtilis\u003c/em\u003e TOA JPC [67]. We also reported the presence of adhesins, which might aid in adhesion to the mucosal layer of host tissues. All these proteins are present in the currently compared probiotic genomes as well. In addition to these, in the current study, we looked for the antibiotic resistance genes, bacteriocins, and folate biosynthesis pathways. Multiple strategies were used in our analysis to comprehensively catalog these domains, nevertheless, these bacteria may harbor genes that confer intrinsic antibiotic resistance. For example, the presence of chloramphenicol, tetracycline and vancomycin resistance in \u003cem\u003eB. toyonensis\u003c/em\u003e has been reported previously [68, 69]. The bile tolerance is reported for \u003cem\u003eB. clausii\u003c/em\u003e Enterogermina in a recent study [74] that substantiates our findings related to the genes exhibiting bile tolerance. We have also explored the occurrence of bacteriocins in different \u003cem\u003eBacillus\u003c/em\u003e probiotics. The presence of different bacteriocins makes these organisms unique. For example, \u003cem\u003eB. clausii\u003c/em\u003e would be effective in \u003cem\u003eS. aureus \u003c/em\u003einfections [26] in skin disorders [50] while \u003cem\u003eB.\u003c/em\u003e \u003cem\u003eamyloliquefaciens\u003c/em\u003e in food industry especially yogurt beverage and probiotic [70], \u003cem\u003eB. toyonensis \u003c/em\u003eshows antimicrobial activity against several other Streptococci [58, 75, 76].\u003c/p\u003e\n\u003cp\u003eOne of the important roles of probiotics is nutrient digestion and energy recovery by producing folate. Folate occurs naturally in food [62] and it is required for the efficiency of DNA replication, repair and methylation process in humans [77, 78]. The folate-producing probiotic strains could possibly confer protection against cancer, inflammation, cardiovascular disease, stress, and depression [62, 63, 77, 78]. This role of probiotics has been extensively studied for their commercial utilization in folate production [62, 79]. We investigated the pathways involved in the production of folic acid in these probiotics. Not surprisingly, we found all the pathway components intact in \u003cem\u003eB. subtilis\u003c/em\u003e which is being engineered for the enhanced synthesis of folate [63]. Some core components of the folate synthesis pathways were present in other \u003cem\u003eBacillus\u003c/em\u003e probiotics, suggesting that they may also be potential sources for the \u003cem\u003ede novo\u003c/em\u003e synthesis of folate.\u003c/p\u003e\n\u003cp\u003eThe composite genome of \u003cem\u003eBacillus clausii\u003c/em\u003e ENTPro and the comparative analysis presented here have helped identify several genes and pathways of interest required for probiotic action and can serve as a starting point for the experimental characterization of these gene products and bacteriocins in order to obtain a deeper understanding into the mechanism of probiotic action of these important microbes.\u0026nbsp;\u0026nbsp;\u0026nbsp;\u0026nbsp;\u0026nbsp;\u0026nbsp;\u0026nbsp;\u0026nbsp;\u0026nbsp;\u0026nbsp;\u0026nbsp;\u003c/p\u003e"},{"header":"Conclusions","content":"\u003cp\u003eThe composite circular genome of \u003cem\u003eBacillus clausii\u003c/em\u003e strain ENTPro, isolated from Enterogermina\u003cstrong\u003e\u0026reg;\u003c/strong\u003e, an oral probiotic, marketed by Sanofi in India, is reported. The genomes of different strains are inseparable as suggested by complete circular assembly using long PacBio reads. \u003cem\u003eB. clausii\u003c/em\u003e ENTPro shares high similarity with probiotic strains of genus \u003cem\u003eB. clausii\u003c/em\u003e whereas it is quite diverse as compared to other \u003cem\u003eBacillus\u003c/em\u003e probiotics. The ML tree based on 25 housekeeping protein sequences, clearly places \u003cem\u003eBacillus clausii\u003c/em\u003e in a separate clade as an outgroup of the \u003cem\u003eBacillus\u003c/em\u003e species.\u003c/p\u003e\n\u003cp\u003eIn this study, we report the genes that are responsible for conferring antibiotic resistance in \u003cem\u003eB. clausii\u003c/em\u003e. We could identify all the antibiotic resistance genes that are indicative of the presence of all four strains in the assembled composite genome. Also, we compared the presence of antibiotic resistance-conferring genes and related pathways in all the probiotic \u003cem\u003eBacillus\u003c/em\u003e genomes. The most important finding of our study is the identification of bacteriocins in \u003cem\u003eBacillus\u003c/em\u003e probiotic genomes which could be directly related to their usage in food and beverage industry. For e.g. gallidermin bacteriocin identified in \u003cem\u003eBacillus clausii\u003c/em\u003e functions against \u003cem\u003eS. aureus\u003c/em\u003e biofilm formation and infections [26, 80]. The bacteriocins in \u003cem\u003eB. amyloliquefaciens\u003c/em\u003e fight against foodborne pathogens [81] which clearly indicates and supports its usage in the yoghurt beverage industry. \u003cem\u003eB. paralicheniformis\u003c/em\u003e secretes bacteriocins to prevent \u003cem\u003eListeria\u003c/em\u003e, \u003cem\u003eS. aureus,\u003c/em\u003e and \u003cem\u003eEnterococcus\u003c/em\u003e borne infections [57, 82].\u003c/p\u003e\n\u003cp\u003eThe other important aspect we studied was the presence of genes necessary for the production of folate. We found that \u003cem\u003eB. subtilis\u003c/em\u003e can produce folate \u003cem\u003ede novo\u003c/em\u003e whereas other \u003cem\u003eBacillus\u003c/em\u003e probiotics depend on supplements viz. pABA to produce the same. Several important components and alternative pathways for folate production were present in other \u003cem\u003eBacillus\u003c/em\u003e probiotics but not complete like \u003cem\u003eB. subtilis\u003c/em\u003e. While identifying several genes and pathways of interest is insufficient to explain the concerted probiotic action, we believe our study shed light on several genomic aspects of different \u003cem\u003eBacillus\u003c/em\u003e probiotics. We trust that the comparative genomics analysis presented here will pave the way for experimental characterization of our findings, and to possibly engineer these organisms for enhanced probiotic actions.\u003c/p\u003e"},{"header":"Methods","content":"\u003cp\u003e\u003cstrong\u003eIsolation and purification of \u003cem\u003eBacillus clausii \u003c/em\u003egenomic DNA:\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eB. clausii\u003c/em\u003e spore suspension drug \u0026ldquo;Enterogermina\u003cstrong\u003e\u0026reg;\u003c/strong\u003e, Sanofi-Aventis\u0026rdquo; (Batch No. 120965; Mfd. 12/2011 and Exp. 11/2013) was procured from a drugstore in Chandigarh, India and was cultured in March 2013. Bacterial cells were suspended in Milli-Q water, serially diluted, and plated on ATCC medium: 688 nutrient agar plates. The plates were incubated at 25\u0026deg;C for 48 hours. DNA isolation was performed using the ZR Fungal/Bacterial DNA miniprep kit (Zymogen) as per instructions in its user manual. After isolation, the genomic DNA was treated with RNase A (1\u0026micro;l of a 10 \u0026micro;g/mL stock solution for 100\u0026micro;l of a solution containing DNA) and incubated at 37\u0026ordm;C for 30 minutes. Then, 1/10 volume of 3M sodium acetate (pH 5.2) and 2.5 volumes of absolute ethanol was added followed by incubation at -20\u0026ordm;C overnight and centrifugation at 14,000 rpm for 30 minutes at 4\u0026ordm;C. The supernatant was carefully discarded; the pellet was rinsed with 70% ethanol and centrifuged again at 14,000 rpm for 15 minutes at 4\u0026ordm;C. The ratio of OD at 260/280 nm was \u0026gt;1.8 as observed by NanoDropND-1000 spectrophotometer.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eGenome Sequencing:\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003ePacBio Sequencing:\u003c/strong\u003e The probiotic \u003cem\u003eB. clausii \u003c/em\u003eENTPro was sequenced using PacBio P6C4 chemistry at Genome Quebec Centre, McGill University. DNA samples were sheared and concentrated using AMPure magnetic beads and treated by ExoVII to remove single-stranded ends. SMRTbell libraries were created using the \u0026lsquo;Procedure and Checklist\u0026ndash;20\u0026thinsp;kb Template Preparation Using BluePippin\u003csup\u003eTM\u003c/sup\u003e Size Selection System protocol. Size Selection was performed to retain longer reads (\u0026gt;10k reads) for sequencing. Blunt ligation reactions were prepared and SMRTbell templates were purified using AMPure magnetic beads. BluePippin\u003csup\u003eTM\u003c/sup\u003e The Size-selected SMRTbell templates were annealed and polymerase was added for Sequencing. Single SMRT cell was run on the PacBio RS II system using P6C4 chemistry and a 180-minute data collection mode.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eIllumina Sequencing: \u003c/strong\u003e\u003cem\u003eB. clausii \u003c/em\u003eENTPro was also sequenced using the Illumina HiSeq PE platform. The library preparation was carried out according to the TruSeq DNA sample preparation protocol (Illumina, Inc., San Diego, CA) at C-CAMP, Bangalore, India. One \u0026mu;g of bacterial DNA was sheared to an average length of 300 to 400 bp. End repair, A-tailing, and adapter ligation (~120 base adapter) procedure was performed according to paired-end DNA sample preparation kit (Illumina, Size selection of adapter-ligated DNA was done in a range of 400 to 550 bases for DNA library). The insert size was taken in a range of 280 to 430 bases for DNA library. PCR enrichment was performed for eight cycles, and the samples were validated on a Bioanalyzer. Libraries were sequenced in a paired-end 100 base run, using TruSeq PE Cluster Kit v3-cBot-HS for cluster generation on C-bot and TruSeq SBS Kit v3-HS (Catalog No.: PE-401-3001) for sequencing on the Illumina HiSeq 1000 platform according to recommended protocols.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eGenome Assembly and Annotation:\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe PacBio reads were assembled \u003cem\u003ede novo\u003c/em\u003e using Hierarchical Genome Assembly Process (HGAP) v2.0 [71] in SMRT portal using default parameters. Functional annotation was carried out by RAST (Rapid Annotation using Subsystem Technology) [83, 84], tRNA was predicted by tRNAscan-SE 1.23 [85] and rRNA genes by RNAmmer 1.2 [86]. The taxonomic characterization of the contigs was performed by subjecting the contigs to BLASTn [87] against NT database. The methylome was deduced by RS Modification and Motif analysis in SMRT portal (\u003ca href=\"https://github.com/PacificBiosciences\"\u003ehttps://github.com/PacificBiosciences\u003c/a\u003e). The plasmid sequence was confirmed by plasmidSPAdes [88].\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003ePhylogenetic Analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe \u003cem\u003eBacillus\u003c/em\u003e genus comprises more than 1000 genomes. To be efficiently able to plot the phylogenetic position of \u003cem\u003eBacillus clausii\u003c/em\u003e in \u003cem\u003eBacillaceae\u003c/em\u003e family, we retrieved the \u003cem\u003eBacillus\u003c/em\u003e genomes that have been classified as representative genomes by NCBI (Accession Numbers in \u003cstrong\u003eTable S4\u003c/strong\u003e). The representative genomes from genus \u003cem\u003eListeria\u003c/em\u003e and \u003cem\u003eClostridium\u003c/em\u003e were selected as outliers. Twenty five housekeeping proteins (ribosomal protein S11, ribosomal protein S13, ribosomal protein S19, ribosomal protein S2, ribosomal protein S3, ribosomal protein S5, ribosomal protein S9, ribosomal protein L11, ribosomal protein L13, ribosomal protein L19, ribosomal protein L2, ribosomal protein L20, ribosomal protein L27, ribosomal protein L3, ribosomal protein L4, ribosomal protein L5, ribosomal protein L6, ribosomal protein L7/L12, CTP synthase, DNA gyrase subunit B, DNA mismatch repair protein MutS, DNA primase, elongation factor Ts, elongation factor Tu, Transcription termination protein NusA) were retrieved from all the genomes and were concatenated in a particular order. The sequences were aligned using Muscle [89] and the phylogenetic inference was drawn using Maximum likelihood [ML] approach based on PROTGAMMA model in RAxML [90] (Bootstrap: 100).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eComparative Genomics\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eFor comparative analysis, all \u003cem\u003eB. clausii\u003c/em\u003e genomes available on March 2018 were downloaded from NCBI. In addition, the genomes with reported probiotics properties were downloaded from the NCBI for comparison. All these genomes were annotated again using the RAST server [83] to remove the bias from different annotation strategies. COGs were identified by subjecting the proteomes of these organisms to BLASTp against COG database [91] at E-value 1e-5.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eIdentification of genome features contributing to probiotic properties of \u003cem\u003eB. clausii \u003c/em\u003eENTPro\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eB. clausii\u003c/em\u003e ENTPro was scanned using Hidden Markov Model (HMM) [92, 93] for the presence of specific domains involved in acid tolerance, Adhesion, antibiotic resistance, antimicrobial production, heavy metal resistance, bile resistance, oxidative and universal stress resistance, and riboflavin synthesis.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eIdentification of Antibiotic Resistance genes\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe Comprehensive Antibiotic Resistance Database (CARD) [94] and Pfam domains were downloaded and hmmscan [93] was run locally against the proteome of all the organisms to identify the domains that could impart antibiotic resistance. Rifampicin resistance was identified based on the mutation in \u003cem\u003eRpoB\u003c/em\u003e genes. Chloramphenicol resistance was identified based on the presence of chloramphenicol acetyltransferase gene in the proteome of the respective organisms. The presence of \u003cem\u003eerm \u003c/em\u003e(34) was found by subjecting its gene sequence to BLASTn against the NR database. Streptomycin biosynthesis and Vancomycin resistance were identified from the pathway analysis against the KEGG database [64].\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eIdentification of Bacteriocins\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eBacteriocins were reported as per the identification from the BAGEL3 server [95]. Any bacteriocin is considered present in a species if that bacteriocin is present at least in one more strain of that species.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eKEGG Pathway Analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe proteomes of all the organisms were subject to bidirectional best hits to KEGG database [64, 96] to identify the components of folate biosynthesis pathways, streptomycin biosynthesis, and vancomycin resistance pathways.\u003c/p\u003e"},{"header":"List of Abbreviations","content":"\u003cp\u003eAA: amino-acid; bp: base pairs; CDS: Coding DNA Sequences; COGs: Clusters of Orthologous Groups; GGDC: Genome-Genome Distance Calculator; GIT: Gastrointestinal Tract; KEGG: Kyoto Encyclopedia of Genes and Genomes; M: modification; NT: nucleotide database; PABA: para-aminobenzoic acid; R: restriction; RRDR: RIF resistance determining region; S: specificity\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eEthics approval and consent to participate\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConsent for publication\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot Applicable\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of supporting data\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis Whole Genome Shotgun project for chromosome and plasmid of \u003cem\u003eBacillus clausii \u003c/em\u003eENTPro has been deposited at DDBJ/EMBL/GenBank under the accession CP012475 and CP012476.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAll the authors declare that they have no conflict of interests.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis work was supported by the Council of Scientific Research (CSIR) Network projects on Man as a super-organism: Understanding the\u0026nbsp;Human\u0026nbsp;Microbiome (HUM-CSIR-BSC-0119) and the Department of Biotechnology project BTISNET (GAP001). IK and GS were supported by research fellowships from the University Grants Commission (UGC) and Council of Scientific and Industrial Research (CSIR) and respectively. \u0026nbsp;The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthors\u0026rsquo; Contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eSS conceived the idea; GS isolated genomic DNA carried out strain identification and performed genome assembly and annotation. IK and GS carried out the comparative analysis. IK, GS, and SS wrote the manuscript. All authors have read and approved the manuscript.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgments\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe thank the next-generation genomics facility at C-CAMP (\u003ca href=\"http://www.ccamp.res.in/\"\u003ehttp://www.ccamp.res.in/\u003c/a\u003e), Bangalore and Genome Quebec Centre, McGill University, Canada for help in obtaining NGS data.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eFuller R. Probiotics in man and animals. J Appl Bacteriol. 1989;66:365\u0026ndash;78. http://www.ncbi.nlm.nih.gov/pubmed/2666378.\u003c/li\u003e\n\u003cli\u003eCutting SM. \u003cem\u003eBacillus\u003c/em\u003e probiotics. Food Microbiol. 2011;28:214\u0026ndash;20. doi:10.1016/j.fm.2010.03.007.\u003c/li\u003e\n\u003cli\u003eNicholson WL, Munakata N, Horneck G, Melosh HJ, Setlow P. Resistance of \u003cem\u003eBacillus\u003c/em\u003e endospores to extreme terrestrial and extraterrestrial environments. Microbiol Mol Biol Rev. 2000;64:548\u0026ndash;72. http://www.ncbi.nlm.nih.gov/pubmed/10974126.\u003c/li\u003e\n\u003cli\u003eCasula G, Cutting SM. \u003cem\u003eBacillus\u003c/em\u003e probiotics: spore germination in the gastrointestinal tract. Appl Env Microbiol. 2002;68:2344\u0026ndash;52. http://www.ncbi.nlm.nih.gov/pubmed/11976107.\u003c/li\u003e\n\u003cli\u003eKim SG, Becattini S, Moody TU, Shliaha P V., Littmann ER, Seok R, et al. Microbiota-derived lantibiotic restores resistance against vancomycin-resistant Enterococcus. Nature. 2019;572:665\u0026ndash;9. doi:10.1038/s41586-019-1501-z.\u003c/li\u003e\n\u003cli\u003eSenok AC, Ismaeel AY, Botta GA. Probiotics: facts and myths. Clin Microbiol Infect. 2005;11:958\u0026ndash;66. doi:10.1111/j.1469-0691.2005.01228.x.\u003c/li\u003e\n\u003cli\u003ePorwal S, Lal S, Cheema S, Kalia VC. Phylogeny in aid of the present and novel microbial lineages: diversity in \u003cem\u003eBacillus\u003c/em\u003e. PLoS One. 2009;4:e4438. doi:10.1371/journal.pone.0004438.\u003c/li\u003e\n\u003cli\u003eSenesi S, Celandroni F, Tavanti A, Ghelardi E. Molecular characterization and identification of \u003cem\u003eBacillus clausii\u003c/em\u003e Strains marketed for use in oral bacteriotherapy. Appl Env Microbiol. 2001;67:834\u0026ndash;9. doi:10.1128/AEM.67.2.834-839.2001.\u003c/li\u003e\n\u003cli\u003eCiffo F. Determination of the spectrum of antibiotic resistance of the \u0026ldquo;Bacillus subtilis\u0026rdquo; strains of Enterogermina. Chemioterapia. 1984;3:45\u0026ndash;52.\u003c/li\u003e\n\u003cli\u003eMazza P, Zani F, Martelli P. Studies on the antibiotic resistance of \u003cem\u003eBacillus subtilis\u003c/em\u003e strains used in oral bacteriotherapy. Boll Chim Farm. 1992;131:401\u0026ndash;8. http://www.ncbi.nlm.nih.gov/pubmed/1299263.\u003c/li\u003e\n\u003cli\u003eGreen DH, Wakeley PR, Page A, Barnes A, Baccigalupi L, Ricca E, et al. Characterization of two \u003cem\u003eBacillus\u003c/em\u003e probiotics. Appl Environ Microbiol. 1999;65:4288\u0026ndash;91. http://www.ncbi.nlm.nih.gov/pubmed/10473456. Accessed 29 Apr 2018.\u003c/li\u003e\n\u003cli\u003eLippolis R, Siciliano RA, Mazzeo MF, Abbrescia A, Gnoni A, Sardanelli AM, et al. Comparative secretome analysis of four isogenic \u003cem\u003eBacillus clausii \u003c/em\u003eprobiotic strains. Proteome Sci. 2013;11:28. doi:10.1186/1477-5956-11-28.\u003c/li\u003e\n\u003cli\u003eRipert G, Racedo SM, Elie A-M, Jacquot C, Bressollier P, Urdaci MC. Secreted Compounds of the Probiotic \u003cem\u003eBacillus clausii \u003c/em\u003eStrain O/C Inhibit the Cytotoxic Effects Induced by \u003cem\u003eClostridium difficile\u003c/em\u003e and \u003cem\u003eBacillus cereus\u003c/em\u003e Toxins. Antimicrob Agents Chemother. 2016;60:3445\u0026ndash;54. doi:10.1128/AAC.02815-15.\u003c/li\u003e\n\u003cli\u003eVarankovich N V., Nickerson MT, Korber DR. Probiotic-based strategies for therapeutic and prophylactic use against multiple gastrointestinal diseases. Front Microbiol. 2015;6:685. doi:10.3389/fmicb.2015.00685.\u003c/li\u003e\n\u003cli\u003eHammad AM, Shimamoto T. Towards a compatible probiotic-antibiotic combination therapy: assessment of antimicrobial resistance in the Japanese probiotics. J Appl Microbiol. 2010;109:1349\u0026ndash;60. doi:10.1111/j.1365-2672.2010.04762.x.\u003c/li\u003e\n\u003cli\u003eLopetuso LR, Scaldaferri F, Franceschi F, Gasbarrini A. \u003cem\u003eBacillus clausii \u003c/em\u003eand gut homeostasis: state of the art and future perspectives. Expert Rev Gastroenterol Hepatol. 2016;:1\u0026ndash;6. doi:10.1080/17474124.2016.1200465.\u003c/li\u003e\n\u003cli\u003ePradhan B, Guha D, Ray P, Das D, Aich P. Comparative Analysis of the Effects of Two Probiotic Bacterial Strains on Metabolism and Innate Immunity in the RAW 264.7 Murine Macrophage Cell Line. Probiotics Antimicrob Proteins. 2016;8:73\u0026ndash;84. doi:10.1007/s12602-016-9211-4.\u003c/li\u003e\n\u003cli\u003ePatrone V, Molinari P, Morelli L. Microbiological and molecular characterization of commercially available probiotics containing \u003cem\u003eBacillus clausii \u003c/em\u003efrom India and Pakistan. Int J Food Microbiol. 2016;237:92\u0026ndash;7. doi:10.1016/j.ijfoodmicro.2016.08.012.\u003c/li\u003e\n\u003cli\u003eBozdogan B, Galopin S, Leclercq R. Characterization of a new erm-related macrolide resistance gene present in probiotic strains of \u003cem\u003eBacillus clausii\u003c/em\u003e. Appl Env Microbiol. 2004;70:280\u0026ndash;4. http://www.ncbi.nlm.nih.gov/pubmed/14711653.\u003c/li\u003e\n\u003cli\u003eMarseglia GL, Tosca M, Cirillo I, Licari A, Leone M, Marseglia A, et al. Efficacy of \u003cem\u003eBacillus clausii \u003c/em\u003espores in the prevention of recurrent respiratory infections in children: a pilot study. Ther Clin Risk Manag. 2007;3:13\u0026ndash;7. doi:10.2147/TCRM.2007.3.1.13.\u003c/li\u003e\n\u003cli\u003eGian Luigi Marseglia, Mariangela Tosca, Ignazio Cirillo, Amelia Licari, Maddalena Leone, Alessia Marseglia, Anna Maria Castellazzi and GC. Efficacy of \u003cem\u003eBacillus clausii \u003c/em\u003espores in the prevention of recurrent respiratory infections in children: a pilot study. Ther Clin Risk Manag. 2007;3:13\u0026ndash;7. https://www.dovepress.com/efficacy-of-bacillus-clausii-spores-in-the-prevention-of-recurrent-res-peer-reviewed-article-TCRM. Accessed 28 Apr 2018.\u003c/li\u003e\n\u003cli\u003eLakshmi SG, Jayanthi N, Saravanan M, Ratna MS. Safety assesment of \u003cem\u003eBacillus clausii\u003c/em\u003e UBBC07, a spore forming probiotic. Toxicol Reports. 2017;4:62\u0026ndash;71. doi:10.1016/J.TOXREP.2016.12.004.\u003c/li\u003e\n\u003cli\u003eCanani RB, Cirillo P, Terrin G, Cesarano L, Spagnuolo MI, De Vincenzo A, et al. Probiotics for treatment of acute diarrhoea in children: randomised clinical trial of five different preparations. BMJ. 2007;335:340. doi:10.1136/bmj.39272.581736.55.\u003c/li\u003e\n\u003cli\u003eVětrovsk\u0026yacute; T, Baldrian P. The Variability of the 16S rRNA Gene in Bacterial Genomes and Its Consequences for Bacterial Community Analyses. PLoS One. 2013;8:e57923. doi:10.1371/journal.pone.0057923.\u003c/li\u003e\n\u003cli\u003eAuch AF, von Jan M, Klenk H-P, G\u0026ouml;ker M. Digital DNA-DNA hybridization for microbial species delineation by means of genome-to-genome sequence comparison. Stand Genomic Sci. 2010;2:117\u0026ndash;34. doi:10.4056/sigs.531120.\u003c/li\u003e\n\u003cli\u003eSaising J, Dube L, Ziebandt A-K, Voravuthikunchai SP, Nega M, G\u0026ouml;tz F. Activity of gallidermin on \u003cem\u003eStaphylococcus aureus\u003c/em\u003e and \u003cem\u003eStaphylococcus epidermidis\u003c/em\u003e biofilms. Antimicrob Agents Chemother. 2012;56:5804\u0026ndash;10. doi:10.1128/AAC.01296-12.\u003c/li\u003e\n\u003cli\u003eGoldstein BP. Resistance to rifampicin: a review. J Antibiot. 2014;67:625\u0026ndash;30. doi:10.1038/ja.2014.107.\u003c/li\u003e\n\u003cli\u003eSchwarz S, Kehrenberg C, It Doublet B, Cloeckaert A. Molecular basis of bacterial resistance to chloramphenicol and florfenicol. 2004. doi:10.1016/j.femsre.2004.04.001.\u003c/li\u003e\n\u003cli\u003eFloch MH, Walker WA. Probiotics are considered nutritional supplements. J Clin Gastroenterol. 2011;45 Suppl:S107. doi:10.1097/MCG.0b013e3182317e74.\u003c/li\u003e\n\u003cli\u003eNavarre WW, Schneewind O. Surface proteins of gram-positive bacteria and mechanisms of their targeting to the cell wall envelope. Microbiol Mol Biol Rev. 1999;63:174\u0026ndash;229. http://www.ncbi.nlm.nih.gov/pubmed/10066836. Accessed 20 Jul 2018.\u003c/li\u003e\n\u003cli\u003eCotter PD, Hill C. Surviving the acid test: responses of gram-positive bacteria to low pH. Microbiol Mol Biol Rev. 2003;67:429\u0026ndash;53, table of contents. doi:10.1128/MMBR.67.3.429-453.2003.\u003c/li\u003e\n\u003cli\u003eBegley M, Gahan CG, Hill C. The interaction between bacteria and bile. FEMS Microbiol Rev. 2005;29:625\u0026ndash;51. doi:10.1016/j.femsre.2004.09.003.\u003c/li\u003e\n\u003cli\u003eAzcarate-Peril MA, Altermann E, Hoover-Fitzula RL, Cano RJ, Klaenhammer TR. Identification and inactivation of genetic loci involved with \u003cem\u003eLactobacillus acidophilus\u003c/em\u003e acid tolerance. Appl Environ Microbiol. 2004;70:5315\u0026ndash;22. doi:10.1128/AEM.70.9.5315-5322.2004.\u003c/li\u003e\n\u003cli\u003eHagenbuch B, Dawson P. The sodium bile salt cotransport family SLC10. Pflugers Arch Eur J Physiol. 2004;447:566\u0026ndash;70. doi:10.1007/s00424-003-1130-z.\u003c/li\u003e\n\u003cli\u003eWong MH, Oelkers P, Craddock AL, Dawson PA. Expression cloning and characterization of the hamster ileal sodium-dependent bile acid transporter. J Biol Chem. 1994;269:1340\u0026ndash;7. http://www.ncbi.nlm.nih.gov/pubmed/8288599. Accessed 20 Jul 2018.\u003c/li\u003e\n\u003cli\u003eNachin L, Nannmark U, Nystrom T. Differential Roles of the Universal Stress Proteins of \u003cem\u003eEscherichia coli\u003c/em\u003e in Oxidative Stress Resistance, Adhesion, and Motility. J Bacteriol. 2005;187:6265\u0026ndash;72. doi:10.1128/JB.187.18.6265-6272.2005.\u003c/li\u003e\n\u003cli\u003eSeifart Gomes C, Izar B, Pazan F, Mohamed W, Mraheil MA, Mukherjee K, et al. Universal stress proteins are important for oxidative and acid stress resistance and growth of \u003cem\u003eListeria monocytogenes\u003c/em\u003e EGD-e in vitro and in vivo. PLoS One. 2011;6:e24965. doi:10.1371/journal.pone.0024965.\u003c/li\u003e\n\u003cli\u003eVentura M, Canchaya C, Zink R, Fitzgerald GF, van Sinderen D. Characterization of the groEL and groES loci in \u003cem\u003eBifidobacterium breve\u003c/em\u003e UCC 2003: genetic, transcriptional, and phylogenetic analyses. Appl Environ Microbiol. 2004;70:6197\u0026ndash;209. doi:10.1128/AEM.70.10.6197-6209.2004.\u003c/li\u003e\n\u003cli\u003eSusin MF, Baldini RL, Gueiros-Filho F, Gomes SL. GroES/GroEL and DnaK/DnaJ have distinct roles in stress responses and during cell cycle progression in \u003cem\u003eCaulobacter crescentus\u003c/em\u003e. J Bacteriol. 2006;188:8044\u0026ndash;53. doi:10.1128/JB.00824-06.\u003c/li\u003e\n\u003cli\u003eVeinger L, Diamant S, Buchner J, Goloubinoff P. The small heat-shock protein IbpB from \u003cem\u003eEscherichia coli\u003c/em\u003e stabilizes stress-denatured proteins for subsequent refolding by a multichaperone network. J Biol Chem. 1998;273:11032\u0026ndash;7. http://www.ncbi.nlm.nih.gov/pubmed/9556585. Accessed 20 Jul 2018.\u003c/li\u003e\n\u003cli\u003eNarberhaus F. Alpha-crystallin-type heat shock proteins: socializing minichaperones in the context of a multichaperone network. Microbiol Mol Biol Rev. 2002;66:64\u0026ndash;93; table of contents. http://www.ncbi.nlm.nih.gov/pubmed/11875128. Accessed 20 Jul 2018.\u003c/li\u003e\n\u003cli\u003eJakob U, Gaestel M, Engel K, Buchner J. Small heat shock proteins are molecular chaperones. J Biol Chem. 1993;268:1517\u0026ndash;20. http://www.ncbi.nlm.nih.gov/pubmed/8093612. Accessed 20 Jul 2018.\u003c/li\u003e\n\u003cli\u003eFu X, Adams Z, Liu R, Hepowit NL, Wu Y, Bowmann CF, et al. Methionine Sulfoxide Reductase A (MsrA) and Its Function in Ubiquitin-Like Protein Modification in \u003cem\u003eArchaea\u003c/em\u003e. MBio. 2017;8. doi:10.1128/mBio.01169-17.\u003c/li\u003e\n\u003cli\u003eMarshall B, Petrowski D, Levy SB, Summers AO. Inter- and intraspecies spread of \u003cem\u003eEscherichia coli\u003c/em\u003e in a farm environment in the absence of antibiotic usage. Proc Natl Acad Sci U S A. 1990;87:6609\u0026ndash;13. doi:10.1073/pnas.87.17.6609.\u003c/li\u003e\n\u003cli\u003eGibson MK, Forsberg KJ, Dantas G. Improved annotation of antibiotic resistance determinants reveals microbial resistomes cluster by ecology. ISME J. 2015;9:207\u0026ndash;16. doi:10.1038/ismej.2014.106.\u003c/li\u003e\n\u003cli\u003eGueimonde M, S\u0026aacute;nchez B, G de Los Reyes-Gavil\u0026aacute;n C, Margolles A. Antibiotic resistance in probiotic bacteria. Front Microbiol. 2013;4:202. doi:10.3389/fmicb.2013.00202.\u003c/li\u003e\n\u003cli\u003eMartinez JL, Baquero F. Mutation frequencies and antibiotic resistance. Antimicrob Agents Chemother. 2000;44:1771\u0026ndash;7. http://www.ncbi.nlm.nih.gov/pubmed/10858329. Accessed 29 Apr 2018.\u003c/li\u003e\n\u003cli\u003eCotter PD, Hill C, Ross RP. What\u0026rsquo;s in a name? Class distinction for bacteriocins - Author reply. Nat Rev Microbiol. 2006;4. doi:DOI 10.1038/nrmicro1273-c2.\u003c/li\u003e\n\u003cli\u003eHeng NCK, Tagg JR. What\u0026rsquo;s in a name? Class distinction for bacteriocins. Nat Rev Microbiol. 2006;4. doi:DOI 10.1038/nrmicro1273-c1.\u003c/li\u003e\n\u003cli\u003eField D, Cotter PD, Hill C, Ross RP. Bioengineering Lantibiotics for Therapeutic Success. Front Microbiol. 2015;6:1363. doi:10.3389/fmicb.2015.01363.\u003c/li\u003e\n\u003cli\u003eMcAuliffe O, Ryan MP, Ross RP, Hill C, Breeuwer P, Abee T. Lacticin 3147, a broad-spectrum bacteriocin which selectively dissipates the membrane potential. Appl Environ Microbiol. 1998;64:439\u0026ndash;45. http://www.ncbi.nlm.nih.gov/pubmed/9464377. Accessed 29 Apr 2018.\u003c/li\u003e\n\u003cli\u003eMasuda Y, Ono H, Kitagawa H, Ito H, Mu F, Sawa N, et al. Identification and characterization of leucocyclicin Q, a novel cyclic bacteriocin produced by \u003cem\u003eLeuconostoc mesenteroides\u003c/em\u003e TK41401. Appl Environ Microbiol. 2011;77:8164\u0026ndash;70. doi:10.1128/AEM.06348-11.\u003c/li\u003e\n\u003cli\u003eShelburne CE, An FY, Dholpe V, Ramamoorthy A, Lopatin DE, Lantz MS. The spectrum of antimicrobial activity of the bacteriocin subtilosin A. J Antimicrob Chemother. 2006;59:297\u0026ndash;300. doi:10.1093/jac/dkl495.\u003c/li\u003e\n\u003cli\u003eLiu Z, Budiharjo A, Wang P, Shi H, Fang J, Borriss R, et al. The highly modified microcin peptide plantazolicin is associated with nematicidal activity of \u003cem\u003eBacillus amyloliquefaciens\u003c/em\u003e FZB42. Appl Microbiol Biotechnol. 2013;97:10081\u0026ndash;90. doi:10.1007/s00253-013-5247-5.\u003c/li\u003e\n\u003cli\u003eKawai Y, Kemperman R, Kok J, Saito T. The circular bacteriocins gassericin A and circularin A. Curr Protein Pept Sci. 2004;5:393\u0026ndash;8. http://www.ncbi.nlm.nih.gov/pubmed/15544534. Accessed 29 Apr 2018.\u003c/li\u003e\n\u003cli\u003eKemperman R, Kuipers A, Karsens H, Nauta A, Kuipers O, Kok J. Identification and characterization of two novel clostridial bacteriocins, circularin A and closticin 574. Appl Environ Microbiol. 2003;69:1589\u0026ndash;97. doi:10.1128/AEM.69.3.1589-1597.2003.\u003c/li\u003e\n\u003cli\u003eBegley M, Cotter PD, Hill C, Ross RP. Identification of a novel two-peptide lantibiotic, lichenicidin, following rational genome mining for LanM proteins. Appl Environ Microbiol. 2009;75:5451\u0026ndash;60. doi:10.1128/AEM.00730-09.\u003c/li\u003e\n\u003cli\u003eSimmonds RS, Pearson L, Kennedy RC, Tagg JR. Mode of action of a lysostaphin-like bacteriolytic agent produced by \u003cem\u003eStreptococcus zooepidemicus\u003c/em\u003e 4881. Appl Environ Microbiol. 1996;62:4536\u0026ndash;41. http://www.ncbi.nlm.nih.gov/pubmed/8953725. Accessed 29 Apr 2018.\u003c/li\u003e\n\u003cli\u003eLiu W-T, Yang Y-L, Xu Y, Lamsa A, Haste NM, Yang JY, et al. Imaging mass spectrometry of intraspecies metabolic exchange revealed the cannibalistic factors of \u003cem\u003eBacillus subtilis.\u003c/em\u003e Proc Natl Acad Sci U S A. 2010;107:16286\u0026ndash;90. doi:10.1073/pnas.1008368107.\u003c/li\u003e\n\u003cli\u003eGonz\u0026aacute;lez-Pastor JE, Hobbs EC, Losick R. Cannibalism by sporulating bacteria. Science. 2003;301:510\u0026ndash;3. doi:10.1126/science.1086462.\u003c/li\u003e\n\u003cli\u003eLin D, Qu LJ, Gu H, Chen Z. A 3.1-kb genomic fragment of \u003cem\u003eBacillus subtilis\u003c/em\u003e encodes the protein inhibiting growth of Xanthomonas oryzae pv. oryzae. J Appl Microbiol. 2001;91:1044\u0026ndash;50. http://www.ncbi.nlm.nih.gov/pubmed/11851812. Accessed 29 Apr 2018.\u003c/li\u003e\n\u003cli\u003eRossi M, Amaretti A, Raimondi S. Folate Production by Probiotic Bacteria. Nutrients. 2011;3:118\u0026ndash;34. doi:10.3390/nu3010118.\u003c/li\u003e\n\u003cli\u003eZhu T, Pan Z, Domagalski N, Koepsel R, Ataai MM, Domach MM. Engineering of \u003cem\u003eBacillus subtilis\u003c/em\u003e for enhanced total synthesis of folic acid. Appl Environ Microbiol. 2005;71:7122\u0026ndash;9. doi:10.1128/AEM.71.11.7122-7129.2005.\u003c/li\u003e\n\u003cli\u003eOgata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 1999;27:29\u0026ndash;34. http://www.ncbi.nlm.nih.gov/pubmed/9847135. Accessed 30 Jul 2017.\u003c/li\u003e\n\u003cli\u003ede Cr\u0026eacute;cy-Lagard V. Variations in metabolic pathways create challenges for automated metabolic reconstructions: Examples from the tetrahydrofolate synthesis pathway. Comput Struct Biotechnol J. 2014;10:41\u0026ndash;50. doi:10.1016/j.csbj.2014.05.008.\u003c/li\u003e\n\u003cli\u003eSalem AR, Pattison JR, Foster MA. Folic acid and the methylation of homocysteine by \u003cem\u003eBacillus subtilis\u003c/em\u003e. Biochem J. 1972;126:993\u0026ndash;1004. http://www.ncbi.nlm.nih.gov/pubmed/4627401. Accessed 20 Jul 2018.\u003c/li\u003e\n\u003cli\u003eKhatri I, Sharma S, Ramya TNC, Subramanian S. Complete Genomes of \u003cem\u003eBacillus coagulans\u003c/em\u003e S-lac and \u003cem\u003eBacillus subtilis\u003c/em\u003e TO-A JPC, Two Phylogenetically Distinct Probiotics. PLoS One. 2016;11:e0156745. doi:10.1371/journal.pone.0156745.\u003c/li\u003e\n\u003cli\u003eCasanovas-Massana A, Sala-Comorera L, Blanch AR. Quantification of tetracycline and chloramphenicol resistance in digestive tracts of bulls and piglets fed with Toyocerin\u0026reg;, a feed additive containing \u003cem\u003eBacillus toyonensis\u003c/em\u003e spores. Vet Microbiol. 2014;173:59\u0026ndash;65. doi:10.1016/J.VETMIC.2014.07.005.\u003c/li\u003e\n\u003cli\u003eZhang S, Hu Y, Fan Q, Wang X, He J. Two-component system YvqEC-dependent bacterial resistance against vancomycin in \u003cem\u003eBacillus thuringiensis\u003c/em\u003e. Antonie Van Leeuwenhoek. 2015;108:365\u0026ndash;76. doi:10.1007/s10482-015-0489-0.\u003c/li\u003e\n\u003cli\u003eAbriouel H, Franz CMAP, Omar N Ben, G\u0026aacute;lvez A. Diversity and applications of \u003cem\u003eBacillus\u003c/em\u003e bacteriocins. FEMS Microbiol Rev. 2011;35:201\u0026ndash;32. doi:10.1111/j.1574-6976.2010.00244.x.\u003c/li\u003e\n\u003cli\u003eChin CS, Alexander DH, Marks P, Klammer AA, Drake J, Heiner C, et al. Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat Methods. 2013;10:563\u0026ndash;9. doi:10.1038/nmeth.2474.\u003c/li\u003e\n\u003cli\u003eUrdaci MC, Bressollier P, Pinchuk I. \u003cem\u003eBacillus clausii\u003c/em\u003e Probiotic Strains. J Clin Gastroenterol. 2004;38 Supplement 2:S86\u0026ndash;90. doi:10.1097/01.mcg.0000128925.06662.69.\u003c/li\u003e\n\u003cli\u003eGhelardi E, Celandroni F, Salvetti S, Gueye SA, Lupetti A, Senesi S. Survival and persistence of \u003cem\u003eBacillus clausii\u003c/em\u003e in the human gastrointestinal tract following oral administration as spore-based probiotic formulation. J Appl Microbiol. 2015;119:552\u0026ndash;9. doi:10.1111/jam.12848.\u003c/li\u003e\n\u003cli\u003eVecchione A, Celandroni F, Mazzantini D, Senesi S, Lupetti A, Ghelardi E. Compositional Quality and Potential Gastrointestinal Behavior of Probiotic Products Commercialized in Italy. Front Med. 2018;5:59. doi:10.3389/fmed.2018.00059.\u003c/li\u003e\n\u003cli\u003eGargis AS, O\u0026rsquo;Rourke A-LD, Sloan GL, Simmonds RS. Prevalence and acquisition of the genes for zoocin A and zoocin A resistance in \u003cem\u003eStreptococcus equi\u003c/em\u003e subsp. \u003cem\u003ezooepidemicus\u003c/em\u003e. J Mol Evol. 2009;68:498\u0026ndash;505. doi:10.1007/s00239-009-9221-x.\u003c/li\u003e\n\u003cli\u003eWilliams LD, Burdock GA, Jim\u0026eacute;nez G, Castillo M. Literature review on the safety of Toyocerin\u0026reg;, a non-toxigenic and non-pathogenic Bacillus cereus var. \u003cem\u003etoyoi\u003c/em\u003e preparation. Regul Toxicol Pharmacol. 2009;55:236\u0026ndash;46. doi:10.1016/J.YRTPH.2009.07.009.\u003c/li\u003e\n\u003cli\u003eFuchs CS, Willett WC, Colditz GA, Hunter DJ, Stampfer MJ, Speizer FE, et al. The influence of folate and multivitamin use on the familial risk of colon cancer in women. Cancer Epidemiol Biomarkers Prev. 2002;11:227\u0026ndash;34. http://www.ncbi.nlm.nih.gov/pubmed/11895870. Accessed 19 Jul 2018.\u003c/li\u003e\n\u003cli\u003eWhite E, Shannon JS, Patterson RE, Hunter DJ, Stampfer MJ, Speizer FE, et al. Relationship between vitamin and calcium supplement use and colon cancer. Cancer Epidemiol Biomarkers Prev. 1997;6:769\u0026ndash;74. http://www.ncbi.nlm.nih.gov/pubmed/9332757. Accessed 19 Jul 2018.\u003c/li\u003e\n\u003cli\u003eSchallmey M, Singh A, Ward OP. Developments in the use of \u003cem\u003eBacillus\u003c/em\u003e species for industrial production. Can J Microbiol. 2004;50:1\u0026ndash;17. doi:10.1139/w03-076.\u003c/li\u003e\n\u003cli\u003eKellner R, JUNG G, HORNER T, ZAHNER H, SCHNELL N, ENTIAN K-D, et al. Gallidermin: a new lanthionine-containing polypeptide antibiotic. Eur J Biochem. 1988;177:53\u0026ndash;9. doi:10.1111/j.1432-1033.1988.tb14344.x.\u003c/li\u003e\n\u003cli\u003eKalyon B, Helaly SE, Scholz R, Nachtigall J, Vater J, Borriss R, et al. Plantazolicin A and B: Structure Elucidation of Ribosomally Synthesized Thiazole/Oxazole Peptides from \u003cem\u003eBacillus amyloliquefaciens\u003c/em\u003e FZB42. Org Lett. 2011;13:2996\u0026ndash;9. doi:10.1021/ol200809m.\u003c/li\u003e\n\u003cli\u003eAlvarez-Ord\u0026oacute;\u0026ntilde;ez A, Begley M, Clifford T, Deasy T, Considine K, O\u0026rsquo;Connor P, et al. Investigation of the Antimicrobial Activity of \u003cem\u003eBacillus licheniformis\u003c/em\u003e Strains Isolated from Retail Powdered Infant Milk Formulae. Probiotics Antimicrob Proteins. 2014;6:32\u0026ndash;40. doi:10.1007/s12602-013-9151-1.\u003c/li\u003e\n\u003cli\u003eAziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, et al. The RAST Server: Rapid Annotations using Subsystems Technology. BMC Genomics. 2008;9:75. doi:10.1186/1471-2164-9-75.\u003c/li\u003e\n\u003cli\u003eOverbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, et al. The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). Nucleic Acids Res. 2014;42 Database issue:D206-14. doi:10.1093/nar/gkt1226.\u003c/li\u003e\n\u003cli\u003eLowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25:955\u0026ndash;64. http://www.ncbi.nlm.nih.gov/pubmed/9023104. Accessed 30 Jul 2017.\u003c/li\u003e\n\u003cli\u003eLagesen K, Hallin P, Rodland EA, Staerfeldt HH, Rognes T, Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007;35:3100\u0026ndash;8. doi:10.1093/nar/gkm160.\u003c/li\u003e\n\u003cli\u003eAltschul SF, Gish W, Miller W, Myers EW, Lipman DJ, Yu Y, et al. Basic local alignment search tool. J Mol Biol. 1990;215:403\u0026ndash;10. doi:10.1016/S0022-2836(05)80360-2.\u003c/li\u003e\n\u003cli\u003eAntipov D, Hartwick N, Shen M, Raiko M, Lapidus A, Pevzner PA. plasmidSPAdes: assembling plasmids from whole genome sequencing data. Bioinformatics. 2016;32:btw493. doi:10.1093/bioinformatics/btw493.\u003c/li\u003e\n\u003cli\u003eEdgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792\u0026ndash;7. doi:10.1093/nar/gkh340.\u003c/li\u003e\n\u003cli\u003eStamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312\u0026ndash;3. doi:10.1093/bioinformatics/btu033.\u003c/li\u003e\n\u003cli\u003eTatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin E V, et al. The COG database: an updated version includes eukaryotes. BMC Bioinformatics. 2003;4:41. doi:10.1186/1471-2105-4-41.\u003c/li\u003e\n\u003cli\u003eEddy SR. Hidden Markov models. Curr Opin Struct Biol. 1996;6:361\u0026ndash;5. http://www.ncbi.nlm.nih.gov/pubmed/8804822.\u003c/li\u003e\n\u003cli\u003eEddy SR. Accelerated Profile HMM Searches. PLoS Comput Biol. 2011;7:e1002195. doi:10.1371/journal.pcbi.1002195.\u003c/li\u003e\n\u003cli\u003eMcArthur AG, Waglechner N, Nizam F, Yan A, Azad MA, Baylay AJ, et al. The comprehensive antibiotic resistance database. Antimicrob Agents Chemother. 2013;57:3348\u0026ndash;57. doi:10.1128/AAC.00419-13.\u003c/li\u003e\n\u003cli\u003evan Heel AJ, de Jong A, Montalb\u0026aacute;n-L\u0026oacute;pez M, Kok J, Kuipers OP. BAGEL3: Automated identification of genes encoding bacteriocins and (non-)bactericidal posttranslationally modified peptides. Nucleic Acids Res. 2013;41 Web Server issue:W448-53. doi:10.1093/nar/gkt391.\u003c/li\u003e\n\u003cli\u003eKanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M. The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004;32:277D \u0026ndash; 280. doi:10.1093/nar/gkh063.\u003c/li\u003e\n\u003c/ol\u003e"},{"header":"Tables","content":"\u003cp style=\"text-align: justify; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003eTable 1:\u003c/strong\u003e Genome assembly and annotation statistics for \u003cem\u003eBacillus clausii\u003c/em\u003e ENTPro composite genome.\u003c/p\u003e\n\u003ctable style=\"width: 475.25pt; border-collapse: collapse;\" width=\"634\"\u003e\n\u003ctbody\u003e\n\u003ctr style=\"height: 1.0pt;\"\u003e\n\u003ctd style=\"width: 475.25pt; border: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" colspan=\"4\" width=\"634\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003eChromosome \u003c/span\u003e\u003c/strong\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003egenome assembly and annotation statistics of \u003cem\u003eBacillus clausii \u003c/em\u003eENTPro\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr style=\"height: 1.0pt;\"\u003e\n\u003ctd style=\"width: 148.35pt; border: solid black 1.0pt; border-top: none; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"198\"\u003e\u0026nbsp;\u003c/td\u003e\n\u003ctd style=\"width: 173.2pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"231\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003eChromosome\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 153.0pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"204\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003ePlasmid\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"border: none; padding: 0in 0in 0in 0in;\" width=\"1\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr style=\"height: 1.0pt;\"\u003e\n\u003ctd style=\"width: 148.35pt; border: solid black 1.0pt; border-top: none; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"198\"\u003e\n\u003cp style=\"text-align: justify; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003eSequencing data\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 326.2pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" colspan=\"2\" width=\"435\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003eP6 polymerase and C4 [P6C4] Chemistry based PacBio sequencing\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"border: none; padding: 0in 0in 0in 0in;\" width=\"1\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr style=\"height: 1.0pt;\"\u003e\n\u003ctd style=\"width: 148.35pt; border: solid black 1.0pt; border-top: none; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"198\"\u003e\n\u003cp style=\"text-align: justify; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003eBio Project Number\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 326.2pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" colspan=\"2\" width=\"435\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003ePRJNA242453\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"border: none; padding: 0in 0in 0in 0in;\" width=\"1\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr style=\"height: 1.0pt;\"\u003e\n\u003ctd style=\"width: 148.35pt; border: solid black 1.0pt; border-top: none; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"198\"\u003e\n\u003cp style=\"text-align: justify; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003eNCBI Accession number\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 173.2pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"231\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003eCP012475\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 153.0pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"204\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003eCP012476\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"border: none; padding: 0in 0in 0in 0in;\" width=\"1\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr style=\"height: 1.0pt;\"\u003e\n\u003ctd style=\"width: 148.35pt; border: solid black 1.0pt; border-top: none; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"198\"\u003e\n\u003cp style=\"text-align: justify; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003eGenome size (in bp)\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 173.2pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"231\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e4,264,866\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 153.0pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"204\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e31,475\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"border: none; padding: 0in 0in 0in 0in;\" width=\"1\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr style=\"height: 1.0pt;\"\u003e\n\u003ctd style=\"width: 148.35pt; border: solid black 1.0pt; border-top: none; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"198\"\u003e\n\u003cp style=\"text-align: justify; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003eGC content (%)\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 173.2pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"231\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e44.75\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 153.0pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"204\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e39.9\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"border: none; padding: 0in 0in 0in 0in;\" width=\"1\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr style=\"height: 1.0pt;\"\u003e\n\u003ctd style=\"width: 148.35pt; border: solid black 1.0pt; border-top: none; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"198\"\u003e\n\u003cp style=\"text-align: justify; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003eChromosome/Contig\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 173.2pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"231\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e1\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 153.0pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"204\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e1\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"border: none; padding: 0in 0in 0in 0in;\" width=\"1\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr style=\"height: 1.0pt;\"\u003e\n\u003ctd style=\"width: 148.35pt; border: solid black 1.0pt; border-top: none; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"198\"\u003e\n\u003cp style=\"text-align: justify; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003eCDS\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 173.2pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"231\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e4,384\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 153.0pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"204\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e40\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"border: none; padding: 0in 0in 0in 0in;\" width=\"1\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr style=\"height: 1.0pt;\"\u003e\n\u003ctd style=\"width: 148.35pt; border: solid black 1.0pt; border-top: none; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"198\"\u003e\n\u003cp style=\"text-align: justify; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e% Coding sequences\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 173.2pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"231\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e86.73\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 153.0pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"204\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e84.45\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"border: none; padding: 0in 0in 0in 0in;\" width=\"1\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr style=\"height: 1.0pt;\"\u003e\n\u003ctd style=\"width: 148.35pt; border: solid black 1.0pt; border-top: none; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"198\"\u003e\n\u003cp style=\"text-align: justify; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003eCDS from (+) strand\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 173.2pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"231\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e2,254\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 153.0pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"204\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e35\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"border: none; padding: 0in 0in 0in 0in;\" width=\"1\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr style=\"height: 1.0pt;\"\u003e\n\u003ctd style=\"width: 148.35pt; border: solid black 1.0pt; border-top: none; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"198\"\u003e\n\u003cp style=\"text-align: justify; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003eCDS from (-) strand\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 173.2pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"231\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e2,130\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 153.0pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"204\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e5\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"border: none; padding: 0in 0in 0in 0in;\" width=\"1\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr style=\"height: 1.0pt;\"\u003e\n\u003ctd style=\"width: 148.35pt; border: solid black 1.0pt; border-top: none; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"198\"\u003e\n\u003cp style=\"text-align: justify; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003eMax. CDS length\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 173.2pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"231\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e9,509\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 153.0pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"204\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e2,711\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"border: none; padding: 0in 0in 0in 0in;\" width=\"1\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr style=\"height: 1.0pt;\"\u003e\n\u003ctd style=\"width: 148.35pt; border: solid black 1.0pt; border-top: none; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"198\"\u003e\n\u003cp style=\"text-align: justify; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003eMean CDS length\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 173.2pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"231\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e843\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 153.0pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"204\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e664\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"border: none; padding: 0in 0in 0in 0in;\" width=\"1\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr style=\"height: 1.0pt;\"\u003e\n\u003ctd style=\"width: 148.35pt; border: solid black 1.0pt; border-top: none; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"198\"\u003e\n\u003cp style=\"text-align: justify; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003eHypothetical proteins\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 173.2pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"231\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e1,215\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 153.0pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"204\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e33\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"border: none; padding: 0in 0in 0in 0in;\" width=\"1\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr style=\"height: 1.0pt;\"\u003e\n\u003ctd style=\"width: 148.35pt; border: solid black 1.0pt; border-top: none; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"198\"\u003e\n\u003cp style=\"text-align: justify; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003eHypothetical proteins (%)\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 173.2pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"231\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e27.72\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 153.0pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"204\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e82.5\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"border: none; padding: 0in 0in 0in 0in;\" width=\"1\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr style=\"height: 1.0pt;\"\u003e\n\u003ctd style=\"width: 148.35pt; border: solid black 1.0pt; border-top: none; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"198\"\u003e\n\u003cp style=\"text-align: justify; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003etRNA\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 173.2pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"231\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e76\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 153.0pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"204\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003eNA\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"border: none; padding: 0in 0in 0in 0in;\" width=\"1\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003ctr style=\"height: 1.0pt;\"\u003e\n\u003ctd style=\"width: 148.35pt; border: solid black 1.0pt; border-top: none; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"198\"\u003e\n\u003cp style=\"text-align: justify; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003erRNA\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 173.2pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"231\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003e7 operons (21 rRNAs)\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"width: 153.0pt; border-top: none; border-left: none; border-bottom: solid black 1.0pt; border-right: solid black 1.0pt; padding: .75pt 4.55pt 0in 4.55pt; height: 1.0pt;\" width=\"204\"\u003e\n\u003cp style=\"text-align: center; line-height: 150%; tab-stops: 7.1pt;\"\u003e\u003cstrong\u003e\u003cspan style=\"font-size: 11.0pt; line-height: 150%;\"\u003eNA\u003c/span\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003c/td\u003e\n\u003ctd style=\"border: none; padding: 0in 0in 0in 0in;\" width=\"1\"\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003c/td\u003e\n\u003c/tr\u003e\n\u003c/tbody\u003e\n\u003c/table\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":false,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"bmc-microbiology","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"mcro","sideBox":"Learn more about [BMC Microbiology](http://bmcmicrobiol.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/mcro","title":"BMC Microbiology","twitterHandle":"#bmcmicrobiology","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"Bacteriocins, Gastrointestinal-tract, Phylogeny, Resistome, Pathogenicity.","lastPublishedDoi":"10.21203/rs.2.15490/v3","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.2.15490/v3","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"Background: Some of the spore-forming strains of Bacillus probiotics are marketed commercially as they survive harsh gastrointestinal conditions and bestow health benefits to the host.\nResults: We report the composite genome of Bacillus clausii ENTPro from a commercially available probiotic Enterogermina® and compare it with the genomes of other Bacillus probiotics. We find that the members of B. clausii species harbor high heterogeneity at the species as well as genus level. The genes conferring resistance to chloramphenicol, streptomycin, rifampicin, and tetracycline in the B. clausii ENTPro strain could be identified. The genes coding for the bacteriocin gallidermin, which prevents biofilm formation in the pathogens S. aureus and S. epidermidis, were also identified. KEGG Pathway analysis suggested that the folate biosynthesis pathway, which depicts one of the important roles of probiotics in the host, is conserved completely in B. subtilis and minimally in Bacillus clausii and other probiotics.\nConclusions: We identified various antibiotic resistance, bacteriocins, stress-related, and adhesion-related domains, and industrially-relevant pathways, in the genomes of these probiotic bacteria that are likely to help them survive in the harsh gastrointestinal tract, facilitating adhesion to host epithelial cells, persistence during antibiotic treatment and combating bacterial infections.","manuscriptTitle":"Composite genome sequence of Bacillus clausii, a probiotic commercially available as Enterogermina®, and insights into its probiotic properties","msid":"","msnumber":"","nonDraftVersions":[{"code":3,"date":"2019-12-11 21:24:18","doi":"10.21203/rs.2.15490/v3","editorialEvents":[{"type":"communityComments","content":0},{"type":"checksComplete","content":"","date":"2019-12-11T12:00:00+00:00","index":"","fulltext":""},{"type":"decision","content":"Accept","date":"2019-12-11T12:00:00+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"bmc-microbiology","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"mcro","sideBox":"Learn more about [BMC Microbiology](http://bmcmicrobiol.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/mcro","title":"BMC Microbiology","twitterHandle":"#bmcmicrobiology","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true}},{"code":2,"date":"2019-11-27 16:49:23","doi":"10.21203/rs.2.15490/v2","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Minor revision","date":"2019-12-04T12:00:00+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2019-11-25T12:00:00+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2019-11-24T12:00:00+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2019-11-24T12:00:00+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"bmc-microbiology","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"mcro","sideBox":"Learn more about [BMC Microbiology](http://bmcmicrobiol.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/mcro","title":"BMC Microbiology","twitterHandle":"#bmcmicrobiology","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true}},{"code":1,"date":"2019-10-01 23:25:07","doi":"10.21203/rs.2.15490/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Minor revision","date":"2019-11-06T12:00:00+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2019-10-29T12:00:00+00:00","index":3,"fulltext":"Form responses:\n---\n* Are the methods appropriate and well described?: **Yes**\n* Does the work include the necessary controls?: **Unable to assess**\n* Are the conclusions drawn adequately supported by the data shown?: **Yes**\n* Are you able to assess any statistics in the manuscript or would you recommend an additional statistical review?: **Not relevant to this manuscript**\n* Quality of written English: **Acceptable**\n* Declaration of competing interests: **I declare that I have no competing interests.**\n* Does the work include the necessary controls?\nIf not, please specify which controls are required in your comments to the authors.\tUnable to assess: **Are the conclusions drawn adequately supported by the data shown?\nIf not, please explain in your comments to the authors.\tYes**\n"},{"type":"editorInvitedReview","content":"","date":"2019-10-16T12:00:00+00:00","index":2,"fulltext":"Recommendation: Reviewer's comments unavailable due to the journal's policy.\n"},{"type":"editorInvitedReview","content":"","date":"2019-10-16T12:00:00+00:00","index":1,"fulltext":"Recommendation: Reviewer's comments unavailable due to the journal's policy.\n"},{"type":"reviewerAgreed","content":"","date":"2019-10-09T12:00:00+00:00","index":3,"fulltext":""},{"type":"reviewerAgreed","content":"","date":"2019-10-07T12:00:00+00:00","index":2,"fulltext":""},{"type":"reviewersInvited","content":"","date":"2019-10-01T12:00:00+00:00","index":"","fulltext":""},{"type":"reviewerAgreed","content":"","date":"2019-10-01T12:00:00+00:00","index":1,"fulltext":""},{"type":"checksComplete","content":"","date":"2019-09-26T12:00:00+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2019-09-22T12:00:00+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2019-09-21T12:00:00+00:00","index":"","fulltext":""},{"type":"submitted","content":"","date":"2019-09-16T12:00:00+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"bmc-microbiology","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"mcro","sideBox":"Learn more about [BMC Microbiology](http://bmcmicrobiol.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/mcro","title":"BMC Microbiology","twitterHandle":"#bmcmicrobiology","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"753b7bfd-f327-4a8e-9828-cba07de9c003","owner":[],"postedDate":"December 11th, 2019","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[{"id":38635,"name":"General Microbiology"}],"tags":[],"updatedAt":"","versionOfRecord":{"articleIdentity":"rs-6155","link":"https://doi.org/10.1186/s12866-019-1680-7","journal":{"identity":"bmc-microbiology","isVorOnly":false,"title":"BMC Microbiology"},"publishedOn":"2019-12-30 12:00:00","publishedOnDateReadable":"December 30th, 2019"},"versionCreatedAt":"2019-12-11 21:24:18","video":"","vorDoi":"10.1186/s12866-019-1680-7","vorDoiUrl":"https://doi.org/10.1186/s12866-019-1680-7","workflowStages":[]},"version":"v3","identity":"rs-6155","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"identity":"rs-6155","version":["v3"]},"buildId":"7rjqhiLT3MXkJMwkYKINL","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}