Complete Genome Sequencing Revealed the Potential Application of a Novel Weizmannia coagulans PL-W Production with Promising Bacteriocins in Food Preservative

Weizmannia coagulans is an important potential probiotic with dual characteristics of Bacillus and Lactobacillus. This study describes a novel Weizmannia coagulans PL-W with excellent antibacterial activity isolated from Mongolian traditional cheese, in which safety and probiotic potential were evaluated by complete genome sequencing. The crude bacteriocins of W. coagulans PL-W showed antibacterial activity against various foodborne pathogens, including Listeria monocytogenes CMCC 54,004, Bacillus cereus ATCC 14,579, and Staphylococcus aureus ATCC 25,923. Moreover, the crude bacteriocins have outstanding stability against pH, temperature, surfactants, and are sensitive to protease. The complete genome sequencing revealed W. coagulans PL-W consists of 3,666,052-base pair (bp) circular chromosomes with a GC content of 46.24% and 3485 protein-coding genes. It contains 84 tRNA, 10 23S rRNA, 10 16S rRNA, and 10 5S rRNA. In addition, no risk-related genes such as acquired antibiotic resistance genes, virulence, and pathogenic factors were identified, demonstrating that W. coagulans PL-W is safe to use. Furthermore, the presence of gene clusters involved in bacteriocin synthesis, adhesion-related genes, and genes contributing to acid and bile tolerance indicate that W. coagulans PL-W is a potential candidate probiotic. Thus, antimicrobial activity and genome characterization of W. coagulans PL-W demonstrate that it has extensive potential applications as a food protective culture.


Introduction
Foodborne diseases are a significant health problem for all of humanity [1]. Various antibiotics have been used for a long time to inhibit the growth of pathogenic bacteria and prevent their threat to human health. However, the overuse of antibiotics has led to increased multidrug resistance in bacteria, which is a global public health crisis that threatens our ability to treat bacterial infections [2,3]. On the other hand, several studies have shown that the human commensal microbiota can be affected by the continuous use of antibiotics [4,5]. Therefore, it is urgent to search for new antimicrobial substances, such as bacteriocin, to prevent the adverse effects of traditional antimicrobial substances. In this regard, several studies have shown that bacteriocins from probiotics which are generally recognized as safe (GRAS), could be used as a substitute for traditional antibiotics in the future [6,7].
Bacteriocins are antimicrobial peptides synthesized from ribosomes produced by bacteria, which inhibit the growth of related (narrow spectrum) or nonrelated (broad spectrum) microorganisms [8]. At present, the bacteriocins from microbials approved for application in food additives are mainly from Nisin and Pediocin PA-1 produced by lactic acid bacteria (LAB) [9,10]. Bacteriocins predominantly exert their antibacterial activity by influencing gene and protein replication, pore formation, and membrane permeabilization [11,12], which could prevent target bacteria from evolving into corresponding drug-resistant strains, a Diameter of inhibition zone: +++, ≥15 mm; ++, ≥12 mm; +, >8 mm; -, no inhibition zone. The diameter of the hole was 8 mm. The diameter of the circular hole was 8 mm.

Isolation of Antimicrobial Substance-Producing W. coagulans and Crude Antimicrobial Substance Preparation
The method of isolation of W. coagulans was based on the previous method with some modifications [27]. Firstly, Mongolian traditional cheese was homogenized in sterilized water and then heated to 80 • C for 10 min. Subsequently, samples were diluted in a gradient and then spread individually on MRS agar with 2 g/L CaCO 3 and incubated at 45 • C for 48 h. Bacterial colonies that showed clear circles on the plates containing CaCO 3 were individually picked and then inoculated on MRS broth medium and incubated at 45 • C under the shaking condition (200 rpm) for 48 h.
The cell-free supernatant (CFS) was obtained by centrifugation at 10,000× g for 10 min, and then ammonium sulfate was slowly added to the CFS to 80% saturation by stirring and at 4 • C overnight. To collect the crude antimicrobial substance, the mixture was centrifuged at 10,000× g for 20 min at 4 • C, and then the precipitate was resuspended in 1 mL of PBS (pH 6.8). The activity of crude antimicrobial substance against L. monocytogenes CMCC 54,004 was determined by the agar well diffusion method. Strains with strong antimicrobial activity against the L. monocytogenes CMCC 54,004 were selected and identified by Gram staining, lactate production capacity, catalase and oxidase activities, ability to grow at different temperatures (45-60 • C), NaCl concentrations (1-5%), and various sugars. The 16S rDNA sequencing was performed for genotype identification as described previously [28]. Subsequently, the MEGA 7.0 software was used to analyze the phylogenetics of strains.
2.3. Characteristics of Antimicrobial Substance Production in W. coagulans PL-W 2.3.1. Kinetics of Growth and Crude Antimicrobial Substance Production in W. coagulans PL-W W. coagulans PL-W was grown in 30 mL of MRS medium at 45 • C for 48 h and then inoculated into MRS medium at a 2% inoculum (v/v). The cell density at 600 nm was determined every 4 h while CFS was collected to obtain the crude antimicrobial substance and assess antibacterial activity by the method described in Section 2.2 [28].

Physicochemical Properties of Crude Antimicrobial Substance
The stability of the crude antimicrobial substance was examined by detecting the changes in antibacterial activity against L. monocytogenes CMCC 54,004 under different conditions [29]. The enzyme stability of the crude antimicrobial substance was evaluated by incubation with lipase, α-amylase, proteinase K, neutral protease, flavor enzyme, trypsin, and pepsin for 2 h under the optimal reaction conditions for each enzyme, while the final concentration of each enzyme was 1 mg/mL. The pH of the crude antimicrobial substance was adjusted from 2 to 10 with 2 M HCL and 2 M NaOH to detect the change in its antibacterial activity. The thermal stability of the crude antimicrobial substance was tested by incubation at 4-121 • C for 10 and 30 min, respectively. To test the crude antimicrobial substance stability with the surfactant, the crude antimicrobial substance was incubated with various chemical reagents 1% (v/v), including Tween 80, EDTA, urea, and SDS at 37 • C for 2 h. The antimicrobial activity of the untreated antimicrobial substance was measured and used as the positive control.

Antimicrobial Spectrum Assay of Crude Bacteriocins
To investigate the antimicrobial spectrum of W. coagulans PL-W, the antimicrobial activity of the crude bacteriocins of W. coagulans PL-W against a range of indicator strains, including food spoilage bacteria and food-borne pathogens (Table 1), was determined by using the pour plate method described by An et al. [30]. The minimal inhibitory concentration (MIC) and minimal bactericidal concentration (MBC) of the crude bacteriocin were calculated by observing the growth of L. monocytogenes CMCC 54004 mixed with various concentrations of the crude bacteriocins [31].

Purification of Bacteriocin
To identify the bacteriocin, firstly, ultrafiltration tubes of 10 KD and 3 KD were used to isolate the crude bacteriocin. The fractions containing proteins larger than 10 KD, between 3-10 KD, and smaller than 3 KD were tested for antibacterial activity. Then, the active fractions were loaded onto a C18 reverse-phase column (5 µm, 4.6 × 250 mm, Agilent, Santa Clara, CA, USA), connected to a reverse-phase high-performance liquid chromatography (RP-HPLC) system, and eluted at 0.5 mL/min flow rate by a linear gradient elution with 95% water-acetonitrile (5-95%) containing 0.1% trifluoroacetic acid (TFA) for 30 min. The different peaks were collected at an absorbance of 280 nm, and then concentrated using 1 KD ultrafiltration tubes for antibacterial activity evaluation. Using Tricine-SDS-PAGE (16.5% separated and 4% concentrated gel), the range of molecular mass of the collected active fractions was analyzed.

Genome Sequencing, Assembly, Annotation, and Classification
W. coagulans PL-W grown to mid-logarithmic phase were collected and genomic DNA were extracted using the Bacterial Genomic DNA Isolation kit (TianGen, Beijing, China) according to the kit instructions. The PromethION platform was used to sequence the W. coagulans PL-W genome. Subsequently, the sequences were mix assembled with Unicycler and corrected with Pilon (parameter: OFF). After removing the redundant part, the Circlator (parameter: fixStart) was used to move the origin of the sequence to the replication start site of the genome to obtain the final genome sequence. The encoding gene was predicted with prodigal, tRNA, rRNA, and other ncRNAs were predicted using trnascan-SE. Interproscan was used for annotation of genome-encoded proteins. BlastP was used to align the encoded proteins to KEGG, RefSeq, and COG databases, and the best results with alignment coverage greater than 30% were retained as annotation results. The circos was used to draw the nuclear genome circle map.
Using W. coagulans PL-W as a reference strain, digital DNA-DNA hybridization was performed using an online tool (http://ggdc.dsmz.de/, (accessed on 5 September 2022)). The W. coagulans PL-W genome in the form of FASTA was submitted to the Type (Strain) Genome Server (TYGS) and the genome was compared with the TYGS database.

Prediction of the Safety of W. coagulans PL-W
The Comprehensive Antibiotic Research Database (CARD) and Resistance Gene Identification Tool (RGI) were used to analyze the presence of drug-resistance genes in W. coagulans PL-W. The VFDB webserver was used to predict putative virulence factors and the pathogen finder webserver was used to predict bacterial pathogenicity.

Prediction of the Probiotic Characteristics of W. coagulans PL-W
The Hidden Markov model (HMM) was used to detect genes related to acid and bile tolerance in the genome of W. coagulans PL-W and various proteins related to adhesion and aggregation were searched in the genome annotation data. AntiSMASH5 and BAGEL4 were used to predict non-ribosomal synthetic secondary metabolites (NRPS) and bacteriocin synthesis gene clusters in the W. coagulans PL-W genome, respectively.

Statistical Analysis
Three parallel groups were set up for each group of experiments. GraphPad Prism software was used for statistical analysis of the experimental data by a one-way ANOVA, and values of p < 0.05 were considered statistically significant.

Screening and Identification of Antimicrobial Substance-Producing Strains
Bacteriocins are a potential alternative to antibiotics since they are safe and it is not easy to produce antibiotic resistance to indicator strains [32]. Bacteriocins produced by Bacillus have a broad inhibitory spectrum, which makes them of excellent research significance [16]. The present study isolated 243 potential lactate-producing strains from Mongolian traditional cheese. Among these strains, strain PL-W possessed the highest antibacterial activity to indicator strains (L. monocytogenes CMCC 54,004) and was selected as the target strain. It was characterized as Gram-positive, spore-forming, rod-shaped, lactate produced, and positive for catalase, indole, and the Voges-Proskauer test, and grew in medium containing 5% sodium chloride and fermented sucrose, glucose, lactose, and arabinose. The phylogenetic analysis according to the strain 16S rRNA sequence showed that the selected strain belonged to the same evolutionary branch as W. coagulans and shared 100% support degree with W. coagulans 683, W. coagulans DSM1 ATCC 7050, and W. coagulans NBRC 12,583 (Supplementary Materials, Figure S1). Thus, we classified the strain PL-W as W. coagulans PL-W.

Crude Antimicrobial Substance Production Properties of W. coagulans PL-W 3.2.1. Kinetics of Growth and Crude Antimicrobial Substance Production in W. coagulans PL-W
Kinetics of crude antimicrobial substance production of W. coagulans PL-W cultured in MRS broth at 45 • C with shaking are presented in Figure 1. The inhibition zone size of crude antimicrobial substance against L. monocytogenes 54,004 was used to detect its production. The result showed that W. coagulans PL-W entered the stable phase after 24 h of culture. Meanwhile, the inhibition zone against L. monocytogenes 54,004 was detected and reached the largest at 32 h, indicating the production of the crude antimicrobial substance may reach the maximum mass. The antibacterial activity of the crude antimicrobial starts to decrease slightly after 32 h, probably due to degradation by other substances secreted by W. coagulans PL-W. Thus, the crude antimicrobial substance obtained from W. coagulans PL-W was cultured at 32 h for further characterization.

Characterization of Crude Antimicrobial Substance
Bacteriocins are antimicrobial peptides synthesized by ribosomes, produced by LAB, which have drawn wide attention owing to their characterization as food antiseptics. For example, Nisin, produced by Lactococcus lactis, has been approved by the Food and Drug Administration (FDA) and used for food preservation since 1950 [33]. Pediocin, bavaricin, leucocin, and sakacin are other bacteriocins produced by LAB that have been approved for use in food in some countries [34,35]. Although some Bacillus species are involved in

Characterization of Crude Antimicrobial Substance
Bacteriocins are antimicrobial peptides synthesized by ribosomes, produced by LAB, which have drawn wide attention owing to their characterization as food antiseptics. For example, Nisin, produced by Lactococcus lactis, has been approved by the Food and Drug Administration (FDA) and used for food preservation since 1950 [33]. Pediocin, bavaricin, leucocin, and sakacin are other bacteriocins produced by LAB that have been approved for use in food in some countries [34,35]. Although some Bacillus species are involved in various food fermentation processes, none of the Bacillus bacteriocins have been approved as food preservatives [36,37]. W. coagulans has been reported as safe by the FDA and the European Union Food Safety Authority (EFSA), which was classified as Bacillus coagulans due to its Bacillus characteristics. Therefore, there are further opportunities for W. coagulans bacteriocins to prove useful for food preservation. In this study, the antimicrobial substance of W. coagulans PL-W was treated with different enzymes to determine its nature. The antimicrobial activity against L. monocytogenes CMCC 54,004 indicated that the antimicrobial substance remained stable after lipase and α-amylase treatment. On the contrary, the antimicrobial substance was sensitive to proteinase K, indicating that the antimicrobial substance is a bacteriocin (Figure 2A). Furthermore, we examined the stability of the crude bacteriocin of W. coagulans PL-W under different conditions, including pH, temperature, and surfactants, to predict its application in the food industry. In the test of pH stability, the antibacterial activity of crude bacteriocins remained above 90% at pH 2-8, which indicates that the crude bacteriocins still have specific activity under acidic, neutral, and weakly alkaline conditions ( Figure 2B). Therefore, bacteriocins produced by W. coagulans PL-W have the potential to be applied to acidic, neutral, and weakly alkaline foods. The antibacterial activities of the bacteriocins were nearly not altered under 4 • C to 40 • C, indicating that bacteriocins may have storage stability. Although the activity decreased after treatment at 60-100 • C, the crude bacteriocins show more than 70% activity making them valuable for heat-processed foods ( Figure 2C). The stability of bacteriocins against surfactants facilitates their use in emulsifying food. Our study showed that the crude bacteriocins had little influence on antibacterial activity after treatment with different surfactants, suggesting that bacteriocins could be used in emulsified food ( Figure 2D).

Antibacterial Spectrum of Crude Bacteriocins
The antibacterial spectrum of crude bacteriocins was determined by assaying for the antimicrobial activity against several indicator strains. The crude bacteriocins displayed a wide range of the antibacterial spectrum, including Gram-positive and Gram-negative bacteria (Table 1). Significantly, some of them are common spoilage bacteria in the food industry, such as L. monocytogenes CMCC 54,004, S. aureus ATCC 25,923, and B. cereus

Antibacterial Spectrum of Crude Bacteriocins
The antibacterial spectrum of crude bacteriocins was determined by assaying for the antimicrobial activity against several indicator strains. The crude bacteriocins displayed a wide range of the antibacterial spectrum, including Gram-positive and Gram-negative bacteria (Table 1). Significantly, some of them are common spoilage bacteria in the food industry, such as L. monocytogenes CMCC 54,004, S. aureus ATCC 25,923, and B. cereus ATCC 14,579. The result of this study suggests that W. coagulans PL-W and its bacteriocin products have great utilization potential in food preservation.

Purification of Bacteriocin
The crude antimicrobial substances were separated into different fractions using 10 KD and 3 KD ultrafiltration tubes. The result showed that only fractions with molecular mass between 3-10 kDa exhibited antibacterial activity ( Figure 3A). Therefore, it was presumed that the molecular mass of the bacteriocin was among 3-10 kDa. Further purification of the active fraction using a C18 column showed that only the third peak was able to inhibit the growth of the indicator strain ( Figure 3B). Tricine-SDS-PAGE analysis of the molecular mass of peak 3 displayed a single band around 7 kDa ( Figure 3C). Consequently, we conclude that W. coagulans PL-W expressed a bacteriocin with a molecular mass near 7 kDa.

General Genome Features of W. coagulans PL-W
W. coagulans are widely used as probiotics in the food industry, medicine, and an breeding [17]. Genetic analyses assist humans in finding potential probiotic strains genetically encoded properties such as bile acid resistance, epithelial adhesion, and teriocins, but also determine whether a strain is safe to use from a genetic perspe [26,38]. There is currently insufficient information on the probiotic effects of W. coag on the genetic basis. Therefore, the present study aims to elucidate the safety and prob properties of W. coagulans PL-W and explore its potential application in food preserv through genomic analysis. Whole genome sequencing was performed using the Pr thION platform. A total of 3,554,753,206 raw reads were used for genome assembly the unicycler assembler, version 0.4.8 (Ryan R. Wick, Victoria, Australia). The com genome of W. coagulans PL-W consists of 2 contigs of 3,666,052 bp with a GC conte 46.24%. None of the plasmid sequences were validated with the Plasmid Finder 2 The protein-coding genes, ribosomal RNA, and transfer RNA of the genome are show Table 2. The functional classification of protein was analyzed in Table 3 through the database. KEGG and GO databases were also used functionally to annotate protein ing genes (Supplementary Materials, Figures S2 and S3). The complete genome i mation of W. coagulans PL-W is shown in Figure 4.

General Genome Features of W. coagulans PL-W
W. coagulans are widely used as probiotics in the food industry, medicine, and animal breeding [17]. Genetic analyses assist humans in finding potential probiotic strains with genetically encoded properties such as bile acid resistance, epithelial adhesion, and bacteriocins, but also determine whether a strain is safe to use from a genetic perspective [26,38]. There is currently insufficient information on the probiotic effects of W. coagulans on the genetic basis. Therefore, the present study aims to elucidate the safety and probiotic properties of W. coagulans PL-W and explore its potential application in food preservation through genomic analysis. Whole genome sequencing was performed using the Prome-thION platform. A total of 3,554,753,206 raw reads were used for genome assembly with the unicycler assembler, version 0.4.8 (Ryan R. Wick, Victoria, Australia). The complete genome of W. coagulans PL-W consists of 2 contigs of 3,666,052 bp with a GC content of 46.24%. None of the plasmid sequences were validated with the Plasmid Finder 2 tool. The protein-coding genes, ribosomal RNA, and transfer RNA of the genome are shown in Table 2. The functional classification of protein was analyzed in Table 3 through the COG database. KEGG and GO databases were also used functionally to annotate protein-coding genes (Supplementary Materials, Figures S2 and S3). The complete genome information of W. coagulans PL-W is shown in Figure 4. genome of W. coagulans PL-W consists of 2 contigs of 3,666,052 bp with a GC content of 46.24%. None of the plasmid sequences were validated with the Plasmid Finder 2 tool. The protein-coding genes, ribosomal RNA, and transfer RNA of the genome are shown in Table 2. The functional classification of protein was analyzed in Table 3 through the COG database. KEGG and GO databases were also used functionally to annotate protein-coding genes (Supplementary Materials, Figures S2 and S3). The complete genome information of W. coagulans PL-W is shown in Figure 4.

Taxonomic Classification and Phylogeny
Phylogenetic analysis at the genome level has contributed to evaluating the diversity of W. coagulans species and the taxonomic status of the species. Type (Strain) Genome Server (TYGS) is a high-throughput platform that is the most advanced genome-based classification platform [39]. The results of TYGS analysis demonstrated that W. coagulans PL-W is similar to W. coagulans ATCC7050 and W. coagulans DSM1 ( Figure 5). Digital DNA-DNA hybridization (DDH) of the W. coagulans PL-W genome was performed using the publicly available genome sequences of three W. coagulans (DSM1, 2-6, ATCC7050, and XZL9) using the Genome-Genome Distance Calculator (GGDC) 3.0 server [40]. The results indicated that W. coagulans PL-W shares 94.21% and 85.32% similarity with W. coagulans XZL9 and W. coagulans 2-6, proving that W. coagulans PL-W is a new strain of Weizmannia. Using JSpeciesWS to calculate average nucleotide identity (ANI) using the method described previously revealed a 98.23% similarity of W. coagulans PL-W with W. coagulans XZL9 [41]. Thus, W. coagulans PL-W was identified as a member of the W. coagulans species based on >70% similarity in DDH and~95% or higher ANI with a similar reference strain. coagulans species based on >70% similarity in DDH and ~95% or higher ANI with a similar reference strain.

Safety Assessment of W. coagulans PL-W
Safety assessment is an essential process for selecting probiotics. According to the recommendations of Qualified Safety Presumption (QPS) approved by the European

Safety Assessment of W. coagulans PL-W
Safety assessment is an essential process for selecting probiotics. According to the recommendations of Qualified Safety Presumption (QPS) approved by the European Food Safety Authority (EFSA), the presence and potential mobility of antibiotic resistance genes should be considered when selecting new probiotics since the increasing resistance of bacteria to antibiotics poses a great threat to human health [42]. The CARD database contains bacterial drug resistance genes from different environmental sources (such as the gut, domestic wastewater, rivers, etc.) and their annotated information such as their resistance spectrum, mechanism of action, ontology, COG, and CDD [43]. Analysis of the W. coagulans PL-W genome from the CARD database showed that W. coagulans PL-W does not contain any antibiotic resistance genes or acquired antimicrobial genes detected. Similarly, phenotypic studies of W. coagulans PL-W sensitivity to the antibiotics were all below the cut-off values mentioned in the EFSA prescribed standards (Table 4). The Virulence Factor Database (VFDB) collects information on the virulence factors of bacterial pathogens [44]. The existence of any virulence genes was not detected on the W. coagulans PL-W genome using the VFDB service. At the same time, we found that W. coagulans PL-W did not exhibit hemolytic zones on the Columbia Agar plates with sheep blood (Supplementary Materials, Figure S4). Based on the above analysis, we promoted that W. coagulans PL-W was presumed to be a safe strain.

Assessment of Probiotic Properties
The ability to endure harsh conditions in the gastrointestinal tract is one of the criteria for selecting probiotics [45]. A large number of genes involved in acid and bile salt tolerance have been identified in the W. coagulans PL-W genome. ATP synthase (F1-F0-ATPase), a group of proteins mainly involved in acid resistance, are known to effectively pump protons out of the cell by hydrolyzing ATP, maintaining a low proton concentration inside the cell, thereby improving acid tolerance and thus maintaining pH homeostasis. There are eight coding F0F1 ATP synthase genes (atpC, atpD, atpG, atpA, atpH, atpF, atpE, and atpB) that were identified in the W. coagulans PL-W genome. In addition, two (Na+/H+) transport protein genes, one (H+/Cl−) antitransporter gene, ClcA, and one (Ca 2+ /H+) antitransporter gene, chaA, were detected in the W. coagulans PL-W genome, which has been shown to play critical roles in pH and homeostasis of cells [46,47]. Bile salt hydrolases belong to the family of glycine hydrolase, which act by binding bile salts to counteract the harmful effects of bile. A gene-encoding cholylglycine hydrolase was identified in the W. coagulans PL-W genome. Two sodium bile acid symporter family genes were also identified, and the presence of these genes was beneficial for the bile tolerance of W. coagulans PL-W [48]. The survival rate of W. coagulans PL-W at 1% bile salt concentration and pH 2.0 further validated the genomic data.
Adhesion to the intestinal epithelium has been considered an important probiotic property since adhesion can promote intestinal colonization while competitively excluding harmful pathogens [49,50]. We searched the genes of adhesion, colonization, mucin binding, flagellar hook, and fibrinogen/fibronectin binding from annotated data. There were 11 genes found to encode adhesion-related proteins in W. coagulans PL-W, including hook-associated protein FlgK and FlgL, flagellar filament capping protein FliD, flagellar hook-basal body complex protein FliE, flagellar hook-length control protein FliK, flagellar hook assembly protein FlgD, flagellar basal body rod protein FlgF, fibronectin/fibrinogenbinding protein (FbpA), DUF817 domains (fibronectin/fibrinogen binding protein), and segregation and condensation protein (scpA and scpB) (fibronectin-binding protein). Flagella can act directly as adhesins and play a key role in colonization by facilitating bacterial motility [51]. Thus, several proteins related to flagellar formation were detected in W. coagulans PL-W, indicating that W. coagulans PL-W may have good motility. Fibronectin is an important multidomain glycoprotein with various adhesion properties and serves as an essential link between cells and their extracellular matrix [52]. Additionally, some similar types of adhesion proteins have been identified in genome-wide analyses of W. coagulans S-Lac and W. coagulans GBI-30, helping the W. coagulans to colonize in host intestines and act as probiotic [50].

Antimicrobial Compound Gene Prediction and Validation
The characteristics of synthesizing antimicrobial compounds play an important role in the competition of survival exclusion when probiotics are colonized in the gut [11]. Genome analysis is expected to pave the way to finding new antimicrobial compounds and understanding the mechanisms of antimicrobial compound production in bacteria [53]. Therefore, BAGEL 4.0 and AntiSMASH 5.0 were used to predict antimicrobial compounds in W. coagulans PL-W. Polyketides (PKs)-T3PK3 biosynthesis, betalactone, and two RiPPlike compounds were identified by antiSMASH. The two RiPP-like compounds were both identified as Circularin A and Amylocyclicin by BAGEL 4.0 ( Figure 6). Circularin A and Amylocyclicin are typically produced as a propeptide after the leader peptide is cleaved, and then ligation between the N-and C-termini results in cyclic antimicrobial peptides. Circularin A produced by Clostridium beijerinckii ATCC 25,752 has been reported to have broadspectrum antibacterial activity [54], and Amylocyclicin produced by Bacillus amyloliticus FZB42 has a wide range of antibacterial activity against Gram-positive bacteria [55]. It is worth noting that the putative Circularin A and Amylocyclicin synthetic gene clusters of W. coagulans PL-W were not found to be similar to any known cluster by antiSMASH servicer. The ABC transporter may be responsible for the transport of mature peptides to protect itself from bacteriocin attack [56], which was both found on the putative Circularin A and Amylocyclicin synthesis gene cluster, suggesting that these two bacteriocins may be secreted in the form of ABC transport. A gene annotated with CirC protein that had a modification function related to bacteriocin cyclization was detected in the Circularin A synthesis gene cluster [57]. Meanwhile, a histidine kinase gene was detected on the putative Amylocyclicin synthesis gene cluster, which predicted that Amylocyclicin synthesis might be regulated by quorum sensing [58]. A transcriptional regulatory protein and competence regulatory protein may be involved in the transcription and cyclization of Amylocyclicin. peptides. Circularin A produced by Clostridium beijerinckii ATCC 25,752 has been reported to have broad-spectrum antibacterial activity [54], and Amylocyclicin produced by Bacillus amyloliticus FZB42 has a wide range of antibacterial activity against Gram-positive bacteria [55]. It is worth noting that the putative Circularin A and Amylocyclicin synthetic gene clusters of W. coagulans PL-W were not found to be similar to any known cluster by antiSMASH servicer. The ABC transporter may be responsible for the transport of mature peptides to protect itself from bacteriocin attack [56], which was both found on the putative Circularin A and Amylocyclicin synthesis gene cluster, suggesting that these two bacteriocins may be secreted in the form of ABC transport. A gene annotated with CirC protein that had a modification function related to bacteriocin cyclization was detected in the Circularin A synthesis gene cluster [57]. Meanwhile, a histidine kinase gene was detected on the putative Amylocyclicin synthesis gene cluster, which predicted that Amylocyclicin synthesis might be regulated by quorum sensing [58]. A transcriptional regulatory protein and competence regulatory protein may be involved in the transcription and cyclization of Amylocyclicin. Tricine-SDS-PAGE analysis of the bacteriocins indicated that its molecular mass was near 7 kDa ( Figure 3C); combining the genomic information, we concluded the bacteriocin extracted from W. coagulans PL-W was Circular A and the entire amino acid sequence was MGLFHVASKFHVSAGIASGVVTAVLHAGTIASIIGAVTVVMSGGVDAILDMGW-TAFIAEVKHLAKEYGKKRAIAW. Although the genome was confirmed for the presence of the gene of Amylocyclicin, purification of crude bacteriocin did not find it. Thus, we speculate that Amylocyclicin may be expressed only under certain specific growth conditions, or Amylocyclicin did not show inhibitory activity against the L. monocytogenes CMCC 54,004 that we used in the purification process. In this regard, it may be possible to use heterologous expression methods in the future to detect the antimicrobial ability of Amylocyclicin or attempt to use other bacteria rather than L. monocytogenes as the indicator strain when purifying the bacteriocins. In general, the expression of Circularin A may provide a competitive advantage for W. coagulans PL-W.

Conclusions
The current study identifies a bacteriocins-producing W. coagulans PL-W from Mongolian traditional cheese. The crude bacteriocin extracted from W. coagulans PL-W by 80% ammonium sulfate not only exhibited a broad antimicrobial spectrum including a Grampositive bacterium and a Gram-negative bacterium, some of which are common food- Tricine-SDS-PAGE analysis of the bacteriocins indicated that its molecular mass was near 7 kDa ( Figure 3C); combining the genomic information, we concluded the bacteriocin extracted from W. coagulans PL-W was Circular A and the entire amino acid sequence was MGLFHVASKFHVSAGIASGVVTAVLHAGTIASIIGAVTVVMSGGVDAILD-MGWTAFIAEVKHLAKEYGKKRAIAW. Although the genome was confirmed for the presence of the gene of Amylocyclicin, purification of crude bacteriocin did not find it. Thus, we speculate that Amylocyclicin may be expressed only under certain specific growth conditions, or Amylocyclicin did not show inhibitory activity against the L. monocytogenes CMCC 54,004 that we used in the purification process. In this regard, it may be possible to use heterologous expression methods in the future to detect the antimicrobial ability of Amylocyclicin or attempt to use other bacteria rather than L. monocytogenes as the indicator strain when purifying the bacteriocins. In general, the expression of Circularin A may provide a competitive advantage for W. coagulans PL-W.

Conclusions
The current study identifies a bacteriocins-producing W. coagulans PL-W from Mongolian traditional cheese. The crude bacteriocin extracted from W. coagulans PL-W by 80% ammonium sulfate not only exhibited a broad antimicrobial spectrum including a Gram-positive bacterium and a Gram-negative bacterium, some of which are common foodborne pathogens, but has good stability to surfactant, heat, and pH. Genome analysis indicated that W. coagulans PL-W is a safe strain to use. In addition, the strain harbors genes encoding two bacteriocins which might ensure W. coagulans PL-W has excellent antibacterial ability. Research of the crude bacteriocin characteristics and genomic sequencing suggested that W. coagulans PL-W is a safe candidate for controlling foodborne pathogens and is promising for use in the food industry. Moreover, the whole genomic results will contribute to further in vitro and in vivo investigations of W. coagulans PL-W to prospect its application as a probiotic.
Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/foods12010216/s1, Figure S1: Phylogenetic analysis of 16S rRNA sequence of W. coagulans PL-W; Figure S2: KEGG category distribution of functional annotation results; Figure S3: WEGO category distribution of functional annotation results. Figure S4: Evaluation of hemolytic activity of W. coagulans PL-W in Columbia blood agar plates.
Funding: This study was supported by the National Natural Science Foundation of China (grant number 32172172).