Integrated genome-based assessment of safety and probiotic characteristics of Lactiplantibacillus plantarum PMO 08 isolated from kimchi

Abstract

Lactiplantibacillus plantarum PMO 08 has been used as a probiotic starter culture for plant-based fermented beverages, with various health-promoting effects such as cholesterol-lowering and anti-inflammatory activities. This study aimed to analyze the genome sequence of Lp. plantarum PMO 08 and identify its safety and probiotic characteristics at the genomic level. For this, complete genome sequencing was conducted to investigate the genes associated with risk and probiotic characteristics by using Pacbio combined with Illumina HiSeq. This bacterial strain has one circular chromosome of 3,247,789 bp with 44.5% G + C content and two plasmids of 50,296 bp with 39.0% G + C content and 19,592 bp with 40.5% G + C content. Orthologous average nucleotide identity analysis showed that PMO 08 belongs to the Lp. plantarum group with 99.14% similarity to Lp. plantarum WCFS1. No deleterious genes were determined in the virulence factor analysis, and no hemolysin activity or secondary bile salt synthesis were detected in vitro test. In the case of antibiotic resistance analysis, PMO 08 was resistant to ampicillin in vitro test, but these genes were not transferable. In addition, the strain showed same carbohydrate utilization with Lp. plantarum WCFS1, except for mannopyranoside, which only our strain can metabolize. The strain also harbors a gene for inositol monophosphatase family protein related with phytate hydrolysis and have several genes for metabolizing various carbohydrate which were rich in plant environment. Furthermore, in probiotic characteristics several genes involved in phenotypes such as acid/bile tolerance, adhesion ability, and oxidative stress response were detected in genome analysis. This study demonstrates that Lp. plantarum PMO 08 harbors several probiotic-related genes (with no deleterious genes) and is a suitable probiotic starter for plant-based fermentation.

Abstract

Lactiplantibacillus plantarum PMO 08 has been used as a probiotic starter culture for plant-based fermented beverages, with various health-promoting effects such as cholesterol-lowering and anti-inflammatory activities. This study aimed to analyze the genome sequence of Lp. plantarum PMO 08 and identify its safety and probiotic characteristics at the genomic level. For this, complete genome sequencing was conducted to investigate the genes associated with risk and probiotic characteristics by using Pacbio combined with Illumina HiSeq. This bacterial strain has one circular chromosome of 3,247,789 bp with 44.5% G + C content and two plasmids of 50,296 bp with 39.0% G + C content and 19,592 bp with 40.5% G + C content. Orthologous average nucleotide identity analysis showed that PMO 08 belongs to the Lp. plantarum group with 99.14% similarity to Lp. plantarum WCFS1. No deleterious genes were determined in the virulence factor analysis, and no hemolysin activity or secondary bile salt synthesis were detected in vitro test. In the case of antibiotic resistance analysis, PMO 08 was resistant to ampicillin in vitro test, but these genes were not transferable. In addition, the strain showed same carbohydrate utilization with Lp. plantarum WCFS1, except for mannopyranoside, which only our strain can metabolize. The strain also harbors a gene for inositol monophosphatase family protein related with phytate hydrolysis and have several genes for metabolizing various carbohydrate which were rich in plant environment. Furthermore, in probiotic characteristics several genes involved in phenotypes such as acid/bile tolerance, adhesion ability, and oxidative stress response were detected in genome analysis. This study demonstrates that Lp. plantarum PMO 08 harbors several probiotic-related genes (with no deleterious genes) and is a suitable probiotic starter for plant-based fermentation.

Citation: Oh YJ, Kim S-A, Yang SH, Kim DH, Cheng Y-Y, Kang JI, et al. (2022) Integrated genome-based assessment of safety and probiotic characteristics of Lactiplantibacillus plantarum PMO 08 isolated from kimchi. PLoS ONE 17(10): e0273986. https://doi.org/10.1371/journal.pone.0273986

Editor: Guadalupe Virginia Nevárez-Moorillón, Universidad Autonoma de Chihuahua, MEXICO

Received: March 2, 2022; Accepted: August 19, 2022; Published: October 3, 2022

Copyright: © 2022 Oh et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: The complete sequences of chromosome and plasmids of L. plantarum PMO 08 were deposited in GenBank under accession number, CP062059.1 (chromosome), CP062060.1 (plasmid 1), and CP062061.1 (plasmid 2).

Funding: NSH designed experiments. SYL provided experimental material and discussion. YJO, SAK, SHY, DHK, YYC, and JIK performed the analysis and statistical analysis. YJO and SAK drafted the manuscript. All authors read and approved the final manuscript. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: Authors YJO, JIK, and SYL are employee of Pulmuone. The other authors have declared that no competing interests exist.

1. Introduction

Probiotics are live organisms, often bacteria, with beneficial health effects such as lowering cholesterol and modulating of gut-microbiota when consumed in sufficient amounts [1]. In the last few decades, many probiotics have been developed from dairy products and used as starter cultures in the fermentation industry. Probiotics from plant-based foods are attracting attention as alternatives due to the increasing limitations of milk-based foods as a result of lactose intolerance and milk allergy for instance [2, 3]. Meanwhile, plant-based products contain anti-nutritional factor, such as phytates, which interfere with mineral absorption [4]. However, fermentation can improve nutritional properties by degrading phytates to improve mineral absorption [5]. For these reasons, scientists have focused on identifying novel food-grade strains with functional properties for plant-based fermentation [6].

Lactiplantibacillus plantarum is found in various environmental niches including the gastrointestinal and vaginal tracts, plants, and many food products [7]. This bacterial strain is generally recognized as safe (GRAS) and has been consumed as a fermented food with documented safety, offering diverse applications [8]. Owing to their high tolerance to gastrointestinal conditions and health-promoting activities, many Lp. plantarum species have been developed as probiotics [9]. Similarly, Lp. plantarum PMO 08 has been developed as a probiotic starter for plant-based fermented beverages [10]. This strain was isolated from kimchi, a Korean fermented vegetable, with high stability under gastrointestinal conditions [11]. In addition, several studies reported its health-promoting activities such as anti-inflammatory in vitro [12], cholesterol-lowering in mice [13], and gut microbiome shaping in human [14]. As many bioinformatics tools have been developed over the past decades, a genetic information of certain bacteria can accelerate probiotic research at the genomic level [15, 16]. Genome analysis can be used for safety assessments, such as the determination of genes for antimicrobial resistance and virulence factors [17]. In addition, probiotic features, such as tolerance to acid and bile salts, adhesion ability to epithelial cells, and antioxidant capacity, can be determined at the genomic level [18]. Furthermore, the availability of complete genome sequences has significantly expanded our understanding of the biology of these microorganisms [19].

The goal of this study was to analyze the whole genome sequence and investigate the safety and probiotic traits of Lp. plantarum PMO 08 at the genomic level. To achieve this, we conducted whole-genome sequencing of this strain using Pacbio combined with Illumina HiSeq. In addition, we analyzed its COG functional categories, average nucleotide identity (ANI)-based phylogenic taxonomy, and the presence of safety-related genes that encode virulence factors and toxins. Furthermore, we demonstrated the biochemical and probiotic properties of the PMO 08 strain in genomic levels and compared with results obtained from various in vitro studies.

2. Materials and methods

2.1. Microorganism and culture conditions

Lactiplantibacillus plantarum PMO 08 was isolated from kimchi [11] and Lp. plantarum WCFS1 was obtained from American Type Culture Collection (ATCC). The genomic DNA of Lp. plantarum PMO 08 was extracted from overnight culture in MRS medium (Difco, Detroit, MI, USA) at 37°C using a bacterial genomic DNA prep kit (Solgent, Daejeon, Korea) according to the manufacturer’s instructions. The prepared genomic DNA was kept at -20°C for future use. The quality and quantity of the isolated DNA were determined by SPECTRA MAX 190 (Molecular Devices, CA, USA).

2.2. Whole-genome sequencing and genome assembly

Genome sequencing based on hybrid technologies was conducted at Macrogen (Seoul, Korea). A PacBio RS II sequencer (Pacific Biosciences, Menlo Park, CA, USA) was used for long-read sequencing to construct the contigs, and a HiSeq 2500 instrument (Illumina, San Diego, CA, USA) was used for short reads to improve the accuracy of the contigs. Sequencing libraries were generated with a SMRTbell template prep kit v3.0 (Pacific Biosciences) for the PacBio RS II and a TruSeq Nano DNA kit (Illumina) for the HiSeq 2500 following the manufacturer’s protocols. The quality of the raw data sequences was checked using FastQC v0.11.9 [20]. The generated reads were assembled de novo using Hierarchical Genome Assembly Process (HGAP v3.0), with a default minimum length of 6 kb [21]. When both ends of the contig overlapped, the contig was regarded as circular, and overlapped regions were manually trimmed. To make the contigs more accurate, a final polishing step was performed using Pilon v1.21 [22]. The completeness of the final assemblies was evaluated using the Benchmarking Universal Single-Copy Orthologs (BUSCO v3.02) software with parameters set to default [23]. Finally, the genome sequence of Lp. plantarum PMO 08 was annotated using the NCBI Prokaryotic Genome Annotation Pipeline (PGAP v4.11) [24]. All software were run with default settings unless otherwise specified. Further COG functional categories were analyzed using the Eggnog-mapper (v5.0) [25].

2.3. ANI-based taxonomy analysis

Species designation of the strain was first determined using 16S rRNA gene using the Ribosomal Database Project Classifier (http://rdp.cme.msu.edu/). Further analysis, Orthologous (Ortho) ANI was calculated using CJ Bioscience’s OrthoANI Tool (OAT), available on the EzBioCloud server (https://www.ezbiocloud.net/tools/orthoani) [26].

2.4. Determination of virulence factors and undesirable genes

The presence of genes encoding virulence factors and toxins in the PMO 08 genome was determined using the virulence factor database (VFDB) available at http://www.mgc.ac.cn/cgi-bin/VFs/v5/main.cgi [27]. Two search criteria that is a stringent search using the cut-off values at >80% identity and >60% coverage and a less stringent search with cut-off values >60% similarity, >60% coverage, and E-values <1e-10 were used to identify the possible virulence genes. In addition, the BlastKoala search tool in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (Release 90.1) available at https://www.kegg.jp/ was used to identify virulence factors and undesirable genes [28].

2.5. Determination of antibiotic resistance

The antibiotic resistance of PMO 08 was investigated using both phenotypic and genotypic methods. The resistance phenotype of the strain was tested as recommended by the European Food Safety Authority (EFSA) [29]. The susceptibility of the strain to seven antimicrobial drugs (ampicillin, gentamicin, kanamycin, erythromycin, clindamycin, tetracycline, and chloramphenicol) was determined. The minimum inhibitory concentration (MIC) for each antimicrobial was evaluated using the microdilution method as described in the international standard ISO 10932:201020. Briefly, the strain was grown on MRS agar plates containing each antibiotic for 16–24 h. The antimicrobial solutions were prepared by dissolving each antimicrobial powder in an appropriate solvent and adjusting its potency, as suggested in the ISO standard. Lacticaseibacillus rhamnosus GG was used as a quality control strain to ensure the performance of the prepared antimicrobial solutions. The MIC of each antimicrobial was determined in triplicate after incubation for 48 h.

Antibiotic resistance genes were searched using the publicly available databases such as BlastKOALA search tool under “Brite ko01504: Antimicrobial resistance genes” [28], the Comprehensive Antibiotic Resistance Database (CARD) [30], and ResFinder 4.0 [31]. The transferability of the antibiotic resistance genes found in the genome was investigated by their locations in two mobile elements: plasmids and bacteriophages. The existence of prophages in the genome was determined using the PHASTER tool available at http://phaster.ca/ [32] and mobile elements were investigated by using Mobile Element Finder available at https://cge.cbs.dtu.dk/services/MobileElementFinder/ [33].

2.6. Biochemical profiles

Carbohydrate utilization patterns were analyzed using the API CHL kit (BioMériux Co., Marcy-l’Étoile, France), according to the manufacturer’s instructions. Briefly, bacterial cells were harvested after cultivation and resuspended in the API 50CH medium. A 120 μL aliquot of the bacterial resuspension was pipetted into each tube of the test strips, and then mineral oil was dropped in the cupules to cover the tubes, followed by incubation at 37°C for 48 h. Fermentation patterns were recorded as positive if the blue indicator in the medium changed to yellow (with the exception of tube #25, which changed to black).

2.7. Probiotic-related genes

The sequence information for different probiotic genes [34] of related strains was obtained from the NCBI database and used to identify the probiotic genes present in Lp. plantarum strains by a sequence similarity search using BLASTP and sequence similarity dropped below 70%.

3. Results and discussion

3.1. Genomic features

To analyze the genetic characteristics of Lp. plantarum PMO 08, genomic DNA was sequenced. In PacBio RS II sequencer, total of 152,212 high-quality reads and 1,601,406,015 bp were produced, with an average read length of 10,520 bp. In case of Illumina Hiseq sequencer, total of 8,257,502 reads and 834,007,702 bp were generated. The assembly process generated one circular chromosome of 3,247,789 bp with 44.51% G + C content as well as two plasmids of 50,296 bp with 39.05% G + C content and 19,592 bp with 40.49% G + C content (Fig 1). In addition, genome analysis revealed that Lp. plantarum PMO 08 contained 2,972 protein-coding genes, 71 tRNA genes, 16 rRNA genes, and four ncRNA genes (Table 1). Notably, no CRISPR sequences were detected. Further COG functional categories are presented in S1 Fig.

[Figure omitted. See PDF.]

Circles represent the following characteristics from the outermost circle to the center: (1) contig information, (2) coding sequences on forward strand, (3) coding sequences on reverse strand, (4) tRNAs and rRNAs, (5) GC skew, and (6) GC ratio. G, guanine; C, cytosine.

[Figure omitted. See PDF.]

3.2. ANI-based taxonomic analysis

ANI is a measure of nucleotide-level genomic similarity between the coding regions of two genomes, and OrthoANI is an improved version that compares orthologous fragment pairs [26]. The taxonomy of the PMO 08 strain was analyzed by comparing its OrthoANI value with that of other lactic acid bacteria (Fig 2). All Lp. plantarum strains showed more than 99% OrthoANI values with PMO 08, with the type strain DSM 20174 showing the highest value (99.79%). In general, two genomes are considered the same species when the ANI value is higher than 95–96% [35]. With this, PMO 08 showed 99.14% OrthoANI values with WCFS1 and 65.96% with LGG, which are commercial probiotics. These results indicate that strain PMO 08 belongs to the species Lp. plantarum. Meanwhile, Lp. argentoratensis which showed 95.48% OrthoANI with Lp. plantarum PMO 08 was belonging to Lactobacillus plantarum subsp. argentoratensis, but recently reclassified as Lp. argentoratensis sp. nov. [36].

[Figure omitted. See PDF.]

Each cell represents the OrthoANI values between the row and the corresponding genomes of the column. The phylogenetic tree on the left represents a phylogenetic distance constructed by CJ Bioscience’s OAT, available on the EzBioCloud server.

3.3 Determination of virulence factor and toxin-encoding genes

The virulence factor database (VFDB) and BlastKoala were used to detect genes that encode virulence factors and toxins that may exist in the PMO 08 genome. No virulence genes were detected under stringent criteria of >80% identity and >60% coverage in VFDB. However, three categories related to virulence factors and undesirable metabolites were detected in BlastKoala (Table 2).

[Figure omitted. See PDF.]

The first category is the bacterial toxin, hemolysin III (hlyIII), which is classified as a membrane-damaging toxin. When the hemolysin activity of PMO 08 was tested using the supernatant and agar plate method, the strain did not show any erythrocyte lysis (Fig 3). In addition, the sequence of hlyIII in PMO 08 was compared with other lactic acid bacteria such as Lp. plantarum WCFS1, JDM1, and 299v, which are commercial probiotics known to have non-hemolytic activity. As result, hlyIII in PMO 08 showed 100% similarity with the genes of these strains. Furthermore, several studies have demonstrated that hlyIII is widespread in lactobacilli, and that these microbes do not show any hemolytic activity [37]. Therefore, we concluded that hlyIII in PMO 08 is not a safety concern.

[Figure omitted. See PDF.]

(a) Hemolysin activity of Lp. plantarum PMO 08 analyzed on a 5% sheep blood agar plate. (b) Hemolysin activity of culture supernatant from Lp. plantarum PMO 08 treated with erythrocytes. Triton X-100 and Listeria monocytogenes used as a positive control and Lp. rhamnosus GG used as a negative control.

The second category includes genes that produce D-lactate, which may cause metabolic stress in patients [38]. In the KEGG database, two genes were identified as D-lactate producers in PMO 08 and their enzymes were D-lactate dehydrogenase and lactate racemase (Table 2). Recently, several studies have reported that both L- and D-lactate can be metabolized in the human body and are regarded as safe for healthy individuals [39]. However, Monroe (2019) [40] reported that D-lactate can cause acidosis in infants and patients with short bowel syndrome [40]. Therefore, a general precaution regarding D-lactate-producing bacteria is necessary for infants or patients with short bowel syndrome [41].

The third category includes genes that produce secondary bile salts that can promote colorectal cancer [42]. Secondary bile salts, such as deoxycholic acid, can be synthesized by bile salt hydrolase (EC 3.5.1.24), which hydrolyzes the amide bond of 7α-dehydroxylated primary bile salt [43]. In the KEGG database, three choloylglycine hydrolases (bile salt hydrolases) were identified as being involved in primary or secondary bile acid biosynthesis in the PMO 08 strain. However, when the bile salt deconjugation activity of PMO 08 was tested on MRS agar plates containing 0.5% taurodeoxycholic acid, the strain showed no activity (Fig 4). These results indicated that they do not contribute to the synthesis of secondary bile salts. According to the guidelines for the evaluation of probiotics in food issued by the FAO/WHO, bile salt hydrolase is one of the desirable properties for probiotics because strains with this enzyme are resistant to gastric conditions [44]. Therefore, we concluded that bile salt hydrolase and its gene in PMO 08 are not a safety concern.

[Figure omitted. See PDF.]

Plates were incubated anaerobically at 37°C for 48 h. Enterococcus faecium ATCC 19434 used as a positive control and Lacticaseibacillus rhamnosus GG used as a negative control.

The last category includes genes that produce biogenic amine such as histamine through amino acid metabolism which can cause hypotension, flushing, and headache [45]. In the KEGG database, no biogenic producing genes were detected.

3.4 Antibiotic resistance

To analyze the antibiotic resistance of Lp. plantarum PMO 08, a MIC test was performed according to the EFSA guidelines (Table 3). Among the eight antibiotics tested, ampicillin exceeded the cutoff value, with a MIC of 4 μg/mL. Ampicillin belongs to the β-lactam antibiotic group and acts as an inhibitor of transpeptidase, which is essential for cell wall assembly [46].

[Figure omitted. See PDF.]

The primary concern regarding antibiotic resistance genes in lactic acid bacteria is the possibility of transferring such genes to other bacteria via mobile elements, such as plasmids or bacteriophages. Therefore, when the MIC value is above the cut-off value, it is necessary to determine whether these genes are intrinsic or acquired resistance through genetic analysis.

The presence of antibiotic resistance genes was analyzed using BlastKoala, CARD database, and ResFinder 4.0, and their transfer possibilities were investigated by PHASTER database and Mobile Element Finder (Tables 4 and 5). In BlastKoala, 14 genes related to antibiotic resistance were detected in the chromosome, but not in the plasmids. In CARD and ResFinder, no antibiotic resistance genes were detected in criteria with strict or perfect (S1 Table). In phage analysis, six phage regions were detected: three in the chromosome, two in plasmid 1, and one in plasmid 2 (Table 5). In addition, in mobile element analysis, nine regions were determined: five ISP2, one ISP1, and three ISLhe30, and all of them were insertion sequence (Table 6). However, antibiotic resistance genes were not present in these phage regions or mobile element sequence indicating a low possibility of transfer to other bacteria.

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

[Figure omitted. See PDF.]

Intrinsic antibiotic resistance is chromosomally encoded and not transferable to other bacteria [47]. In case of vancomycin, the resistance or sensitivity is from intrinsic nature of bacteria such as Enterococcus, Lactobacillus, and Bifidobacterium; their peptidoglycan cell walls have D-Ala-D-Lac termini that have low binding affinity to vancomycin [48, 49].

In contrast, acquired antibiotic resistance located on mobile genetic elements, such as plasmids or phage regions, can be transferred from one bacterium to another, contributing to a greater and more widespread risk than intrinsic resistance [47]. When we analyzed the presence of antibiotic resistance genes in the genome of the PMO 08 strain, we found that they were located on the chromosome, but not in the plasmid or phage regions. Therefore, we concluded that PMO 08 has a low possibility of transferring antibiotic-resistance genes to other bacteria.

3.5 Carbohydrate utilization pathway and genes

To investigate the metabolic characteristics of Lp. plantarum PMO 08, carbohydrate utilization capacity was analyzed using 49 substrates (Table 6) and matched at the genome level using KEGG pathway analysis (Fig 5 and S2 Table). PMO 08 and WCFS1 commonly utilized L-arabinose, D-ribose, D-galactose, D-glucose, D-fructose, D-mannose, mannitol, sorbitol, N-acetyl-glucosamine, amygdalin, arbutin, esculin, salicin, cellobiose, maltose, lactose, melibiose, sucrose, trehalose, melizitose, raffinose, gentibiose, turanose, and tagatose. Notably, PMO 08 could utilize both α-methyl-D-mannoside and gluconate, whereas WCFS1 could not. Accordingly, this result indicates that PMO 08 has different biochemical attributes than WCFS1 in terms of carbohydrate utilization patterns. In the glycolysis/gluconeogenesis pathway, PMO 08 possesses a PTS system and glucokinase to utilize glucose. In addition, this strain possesses a fructose PTS system, fructose kinase, mannose PTS system, and galactose kinase for sugar intake. These results match the carbohydrate utilization patterns seen in Table 7.

[Figure omitted. See PDF.]

PTS, phosphotransferase; PP, phosphorylase.

[Figure omitted. See PDF.]

However, PMO 08 could not consume L-rhamnose, even though this strain shows a complete pathway to utilize this deoxy sugar. Moreover, this strain has an incomplete pathway to utilize tagatose, but can metabolize that carbon source. In the pentose phosphate pathway, PMO 08 harbors ribokinase to metabolize ribose. However, it does not carry the complete pathway to consume the D/L- form of arabinose or xylose, and these results correspond to those in Table 7, except for L-arabinose. PMO 08 has an incomplete L-arabinose pathway but can utilize that carbon source in vitro. For disaccharide and oligosaccharide metabolism, PMO 08 uses a sucrose PTS system, cellobiose PTS system, maltose phosphorylase, and trehalose phosphorylase to utilize these carbon sources. In addition, this strain possesses α- or β-galactosidase to utilize lactose, melibiose, and raffinose. These results correspond with those of the in vitro tests shown in Table 7.

However, in the case of polysaccharides, PMO 08 has a complete pathway for utilizing starch and glycogen, but the strain cannot metabolize these polysaccharides. In addition, this strain does not carry the genes to consume inulin, and this result is consistent with the experimental data in Table 7. In the case of sugar alcohols, PMO 08 has a mannitol and sorbitol PTS system to utilize these carbon sources. However, it possesses an incomplete pathway to use ribitol, galactitol, inositol, xylitol, and D/L- form of arabitol, and these results were matched with the in vitro test in Table 7, except for galactitol. Although the strain can metabolize galactitol in the experiment, the related genes were not found in the KEGG analysis. In amino sugar and nucleotide sugar metabolism, PMO 08 has a sugar PTS system to metabolize N-acetylglucosamine, which corresponds with the in vitro test. Lastly, PMO 08 carries several β-glucosidase enzymes that use glucosides such as amygdalin, arbutin, esculin, and salicin, which is consistent with the results shown in Table 7. PMO 08 could utilize α-methyl-D-mannoside, whereas WCFS1 could not. Considering the origin of PMO 08 (kimchi, [10]), which is rich in hemicellulose, compared to that of WCFS1 (human saliva, [50]). PMO 08 has multiple copies of mannose PTS system EII components while WCFS1 has one copy of that gene (S2 Table). In addition, the PMO 08 harbors an inositol monophosphatase family protein, which is related with phytate hydrolysis [51]. Taken together, the experimental results were consistent with genome analysis in many carbon sources, but some (L-arabinose, rhamnose, galactitol, starch, glycogen, and D-tagatose) were not matched; thus, further studies are required.

3.6 Probiotic characteristics

To verify the probiotic characteristics of Lp. plantarum PMO 08 at the genomic level, genes related to probiotic properties were annotated using NCBI Prokaryotic Genome Annotation Pipeline (PGAP v4.11). As shown in Table 8, several genes related to acid or bile salt tolerance, adhesion to epithelial cells, and stress scavenging were found on the chromosome of PMO 08.

[Figure omitted. See PDF.]

Na+/H+ antiporters, cation:proton antiporters, F0F1 ATP synthase subunits, and chaperonin GroEL were annotated as contributing resistance to low pH in PMO 08. These genes are associated with the regulation of cytoplasmic pH by hydrolyzing ATP to pump H+ from cells, which allows for the maintenance of pH homeostasis [52]. Generally, F0F1 ATP synthase subunits are highly expressed under low pH conditions [53]. Choloylglycine hydrolase (bile salt hydrolase) and GCN5-N-acetyltransferase (GNAT) proteins were found to be related to bile salt tolerance in PMO 08. Bile salt hydrolase catalyzes the hydrolysis of bile salts and is thought to be involved in the bile salt resistance of bacteria [54]. In addition, studies have reported that bile salt hydrolase can lower cholesterol levels by facilitating the integration of cholesterol or bile into the bacterial membrane [55]. GNAT family proteins contribute to the maintenance of resistance to bile salts through acetylation of the outer membrane glycolipid enterobacterial common antigen [56].

Another essential probiotic characteristic is the ability to adhere to the epithelial membrane. In the genome analysis of PMO 08, elongation factor, LPXTG cell wall anchor domain-containing protein, tyrosine protein phosphatase, CpsD/CapB family tyrosine-protein kinase, and flagellar biosynthesis proteins were annotated as adhesion-related proteins. Among them, the elongation factor Tu acts as an aminoacyl tRNA binding to the active site of ribosomes and this indicates the adhesive ability of PMO 08 to mucin [57]. LPXTG cell wall anchor domain-containing proteins or tyrosine-protein phosphatases act as mucus adhesins to mediate binding to human epithelial cells or adhesive signaling functions [58, 59]. Finally, several antioxidant-related genes, such as triose-phosphate isomerase, thiol peroxidase, glutathione peroxidase, and thioredoxin family proteins were detected in the PMO 08 strain. In case of triose-phosphate isomerase, it functions as an antioxidant via glutathionylation [60]. Furthermore, the thioredoxin-thioredoxin reductase and GSH- glutaredoxin systems, which maintain intracellular dithiol/disulfide homeostasis in lactic acid bacteria, play an important role in the defense against oxidative stress. The thioredoxin system, comprising NADPH, thioredoxin reductase, and thioredoxin, shuttles electrons to thiol-dependent peroxidases to maintain redox homeostasis and protect probiotic bacteria from ROS and RNS damage [61]. This system controls the thiol-disulfide balance and thus plays an essential role in DNA and protein repair by reducing ribonucleotide reductase, methionine sulfoxide reductases, and regulating the activity of numerous redox-sensitive transcription factors [62]. Overall, these results suggest that Lp. plantarum PMO 08 has probiotic characteristics and defense mechanism against oxidative stress that can be verified using in silico genomic analysis.

4. Conclusion

Lp. plantarum PMO 08 has been used as a plant-based probiotic starter, but further investigation is limited because of the absence of information regarding whole genome sequences. In this study, we analyzed the whole genome sequence of the PMO 08 strain and found that it contains one chromosome and two plasmids with a high homology (≥ 99%) of ANI value to Lp. plantarum WCFS1. We conducted safety analysis both in vitro and in silico, and no deleterious genes were detected. In the case of antibiotic resistance genes, PMO 08 was resistant to ampicillin, but these genes were not transferable. In addition, PMO 08 can utilize various carbohydrates that are favorable for adaptation to the plant environment, and these results are consistent with carbohydrate utilization enzymes at the genomic level. The strain also harbors a gene for inositol monophosphatase family protein related with phytate hydrolysis. Furthermore, we verified probiotic characteristics, such as acid or bile salt tolerance, adhesion to epithelial cells, and stress scavenging of PMO 08 at the genomic level. This study demonstrates that Lp. plantarum PMO 08 contains several probiotic-related genes (while harboring no deleterious genes) and is a suitable probiotic starter for plant-based fermentation. This information will provide considerable insight into the physiology of this organism and make it possible to understand their genetic content, predict important capabilities, and manipulate it for improved beneficial activities.

Supporting information

S1 Fig. COG analysis of annotated genes for Lactiplantibacillus plantarum PMO 08.

Functional categories were presented on the right part of the figure with blanket which shows number of genes and percentage.

https://doi.org/10.1371/journal.pone.0273986.s001

(DOCX)

S1 Table. Antimicrobial resistance genes determined by CARD database and ResFinder sever 4.1.

https://doi.org/10.1371/journal.pone.0273986.s002

(XLSX)

S2 Table. Comparison of enzyme related with carbohydrate utilization between Lactiplantibacillus plantarum PMO 08 and Lp. plantarum WCFS1.

https://doi.org/10.1371/journal.pone.0273986.s003

(XLSX)

Acknowledgments

This study was supported by Pulmuone Co.

Citation: Oh YJ, Kim S-A, Yang SH, Kim DH, Cheng Y-Y, Kang JI, et al. (2022) Integrated genome-based assessment of safety and probiotic characteristics of Lactiplantibacillus plantarum PMO 08 isolated from kimchi. PLoS ONE 17(10): e0273986. https://doi.org/10.1371/journal.pone.0273986

About the Authors:

Young Joo Oh

Contributed equally to this work with: Young Joo Oh, Seul-Ah Kim

Roles: Investigation, Writing – original draft

Affiliation: Pulmuone Co., Ltd., Cheongju, Republic of Korea

Seul-Ah Kim

Contributed equally to this work with: Young Joo Oh, Seul-Ah Kim

Roles: Formal analysis, Investigation, Writing – original draft

Affiliation: Brain Korea 21 Center for Bio-Health Industry, Division of Animal, Horticultural, and Food Sciences, Chungbuk National University, Cheongju, Chungbuk, Republic of Korea

Soo Hwi Yang

Roles: Formal analysis, Methodology

Affiliation: Brain Korea 21 Center for Bio-Health Industry, Division of Animal, Horticultural, and Food Sciences, Chungbuk National University, Cheongju, Chungbuk, Republic of Korea

Da Hye Kim

Roles: Formal analysis

Affiliation: Brain Korea 21 Center for Bio-Health Industry, Division of Animal, Horticultural, and Food Sciences, Chungbuk National University, Cheongju, Chungbuk, Republic of Korea

Ya-Yun Cheng

Roles: Investigation, Methodology, Software

Affiliation: Brain Korea 21 Center for Bio-Health Industry, Division of Animal, Horticultural, and Food Sciences, Chungbuk National University, Cheongju, Chungbuk, Republic of Korea

Jung Il Kang

Roles: Project administration, Resources

Affiliation: Pulmuone Co., Ltd., Cheongju, Republic of Korea

Sang Yun Lee

Roles: Investigation, Resources

Affiliation: Pulmuone Co., Ltd., Cheongju, Republic of Korea

Nam Soo Han

Roles: Conceptualization, Writing – review & editing

E-mail: [email protected]

Affiliation: Brain Korea 21 Center for Bio-Health Industry, Division of Animal, Horticultural, and Food Sciences, Chungbuk National University, Cheongju, Chungbuk, Republic of Korea

https://orcid.org/0000-0002-8460-8165

References

1. Hill C, Guarner F, Reid G, Gibson GR, Merenstein DJ, Pot B, et al. Expert consensus document: the international scientific association for probiotics and prebiotics consensus statement on the scope and appropriate use of the term probiotic. Nat Rev Gastro Hepat. 2014.

2. Choi AR, Patra JK, Kim WJ, Kang SS. Antagonistic activities and probiotic potential of lactic acid bacteria derived from a plant-based fermented food. Front Microbiol. 2018;9:1963. pmid:30197633

3. Peres CM, Peres C, Hernández-Mendoza A, Malcata FX. Review on fermented plant materials as carriers and sources of potentially probiotic lactic acid bacteria–with an emphasis on table olives. Trends Food Sci Technol. 2012;26(1):31–42.

4. Jo YM, Kim GY, Kim SA, Cheon SW, Kang CH, Han NS. Limosilactobacillus fermentum MG7011: An amylase and phytase producing starter for the preparation of rice-based probiotic beverages. Front Microbiol. 2021;2895.

5. Rasika DM, Vidanarachchi JK, Rocha RS, Balthazar CF, Cruz AG, Sant’Ana AS, et al. Plant-based milk substitutes as emerging probiotic carriers. Curr Opin Food Sci. 2021;38:8–20.

6. Lin MN, Chiu TH, Chang CE, Lin MN. The impact of a plant-based dietary pattern on dementia risk: a prospective cohort study. Innov Aging. 2019;3(Suppl 1):S734.

7. Liu CJ, Wang R, Gong FM, Liu XF, Zheng HJ, Luo YY, et al. Complete genome sequences and comparative genome analysis of Lactobacillus plantarum strain 5–2 isolated from fermented soybean. Genomics. 2015;106(6):404–11.

8. Seddik HA, Bendali F, Gancel F, Fliss I, Spano G, Drider D. Lactobacillus plantarum and its probiotic and food potentialities. Probiotics Antimicrob Proteins. 2017;9(2):111–22.

9. Behera SS, Ray RC, Zdolec N. Lactobacillus plantarum with functional properties: an approach to increase safety and shelf-life of fermented foods. Biomed Res Int. 2018.

10. Oh YJ, Kim TS, Moon HW, Lee SY, Lee SY, Ji GE, et al. Lactobacillus plantarum PMO 08 as a Probiotic starter culture for plant-based fermented beverages. Molecules. 2020;25(21):5056. pmid:33143293

11. Jang SE, Hyun YJ, Trinh HT, Han MJ, Kim DH. Anti-scratching behavioral effect of Lactobacillus plantarum PM008 isolated from kimchi in mice. Immunopharmacol Immunotoxicol. 2011;33(3):539–44. pmid:21275795

12. Jang SE, Hyun YJ, Oh YJ, Choi KB, Kim T, Yeo IH, et al. Adhesion activity of Lactobacillus plantarum PM 008 isolated from kimchi on the intestine of mice. J Bacteriol Virol. 2011;41(2):83–90.

13. Oh YJ, Kim HJ, Kim TS, Yeo IH, Ji GE. Effects of Lactobacillus plantarum PMO 08 alone and combined with chia seeds on metabolic syndrome and parameters related to gut health in high-fat diet-induced obese mice. J Med Food. 2019;22(12):1199–207.

14. Oh YJ, Nam K, Kim Y, Lee SY, Kim HS, Kang JI, et al. Effect of a nutritionally balanced diet comprising whole grains and vegetables alone or in combination with probiotic supplementation on the gut microbiota. Prev Nutr Food Sci. 2021;26(2):121. pmid:34316477

15. Ventura M, Turroni F, Zomer A, Foroni E, Giubellini V, Bottacini F, et al. The Bifidobacterium dentium Bd1 genome sequence reflects its genetic adaptation to the human oral cavity. PLoS Genet. 2009;5(12):e1000785.

16. Castro-López C, Santiago-López L, Vallejo-Cordoba B, González-Córdova AF, Liceaga AM, García HS, et al. An insight to fermented edible insects: a global perspective and prospective. Food Res Int. 2020;137:109750. pmid:33233312

17. Chokesajjawatee N, Santiyanont P, Chantarasakha K, Kocharin K, Thammarongtham C, Lertampaiporn S, et al. Safety assessment of a nham starter culture Lactobacillus plantarum BCC9546 via whole-genome analysis. Sci Rep. 2020;10(1):1–12.

18. Kapse N, Engineer A, Gowdaman V, Wagh S, Dhakephalkar P. Functional annotation of the genome unravels probiotic potential of Bacillus coagulans HS243. Genomics. 2019;111(4):921–9.

19. Ventura M, Turroni F, van Sinderen D. Probiogenomics as a tool to obtain genetic insights into adaptation of probiotic bacteria to the human gut. Bioengineered. 2012;3(2):73–9.

20. Andrews S, Krueger F, Seconds-Pichon A, Biggins F, Wingett S. Fast QC: a quality control tool for high throughput sequence data. Babraham Bioinformatics. 2015

21. Hellström D, Saghir M. Packaging and logistics interactions in retail supply chains. Packag Technol Sci. 2007;20(3), 197–216.

22. Walker BJ, Abeel T, Shea T, Priest M, Abouelliel A, Sakthikumar S, et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PloS One. 2014; 9(11), e112963. pmid:25409509

23. Simão FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinform. 2015. 31(19), 3210–3212. pmid:26059717

24. Tatusova T, DiCuccio M, Badretdin A, Chetvernin V, Nawrocki EP, Zaslavsky L, et al. NCBI Prokaryotic genome annotation pipeline. Nucleic Acids Res. 2016;44:6614–6624. pmid:27342282

25. Huerta-Cepas J, Szklarczyk D, Heller D, Hernández-Plaza A, Forslund SK, Cook H. EggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 2019;47:D309–14. pmid:30418610

26. Lee I, Kim YO, Park SC, Chun J. OrthoANI: an improved algorithm and software for calculating average nucleotide identity. Int J Syst Evol Microbiol. 2016;66(2):1100–3. pmid:26585518

27. Chen L, Yang J, Yu J, Yao Z, Sun L, Shen Y, et al. VFDB: a reference database for bacterial virulence factors. Nucleic Acids Res. 2005;33(suppl_1):D325–D8. pmid:15608208

28. Kanehisa M, Sato Y, Morishima K. BlastKOALA and GhostKOALA: KEGG tools for functional characterization of genome and metagenome sequences. J Mol Biol. 2016;428(4):726–31. pmid:26585406

29. Leuschner RG, Robinson TP, Hugas M, Cocconcelli PS, Richard-Forget F, Klein G, et al. Qualified presumption of safety (QPS): a generic risk assessment approach for biological agents notified to the European Food Safety Authority (EFSA). Trends Food Sci Technol. 2010;21(9):425–35.

30. Alcock BP, Raphenya AR, Lau TT, Tsang KK, Bouchard M, Edalatmand A, et al. CARD 2020: antibiotic resistome surveillance with the comprehensive antibiotic resistance database. Nucleic Acids Res. 2020;48(D1), D517–D525. pmid:31665441

31. Bortolaia V, Kaas RS, Ruppe E, Roberts MC, Schwarz S, Cattoir V, et al. ResFinder 4.0 for predictions of phenotypes from genotypes. J Antimicrob Chemother. 2020;75(12), 3491–3500. pmid:32780112

32. Arndt D, Marcu A, Liang Y, Wishart DS. PHAST, PHASTER and PHASTEST: tools for finding prophage in bacterial genomes. Brief Bioinform. 2019;20(4):1560–7. pmid:29028989

33. Johansson MH, Bortolaia V, Tansirichaiya S, Aarestrup FM, Roberts AP, Petersen TN. Detection of mobile genetic elements associated with antibiotic resistance in Salmonella enterica using a newly developed web tool: MobileElementFinder. J Antimicrob Chemother. 2021;76(1), 101–109. pmid:33009809

34. Alayande KA, Aiyegoro OA, Nengwekhulu TM, Katata-Seru L, Ateba CN. Integrated genome-based probiotic relevance and safety evaluation of Lactobacillus reuteri PNW1. PLoS One. 2020;15(7), e0235873. pmid:32687505

35. Yoon SH, Ha SM, Lim J, Kwon S, Chun J. A large-scale evaluation of algorithms to calculate average nucleotide identity. Antonie Van Leeuwenhoek. 2017;110(10):1281–6. pmid:28204908

36. Li TT, Liu DD, Fu ML, Gu CT. Proposal of Lactobacillus kosoi Chiou et al. 2018 as a later heterotypic synonym of Lactobacillus micheneri McFrederick et al. 2018, elevation of Lactobacillus plantarum subsp. argentoratensis to the species level as Lactobacillus argentoratensis sp. nov., and Lactobacillus zhaodongensis sp. nov., isolated from traditional Chinese pickle and the intestinal tract of a honey bee (Apis mellifera). Int J Syst Evol Microb. 2020;70(5), 3123–3133.

37. Zhang L, Ma H, Kulyar MFEA, Pan H, Li K, Li A, et al. Complete genome analysis of Lactobacillus fermentum YLF016 and its probiotic characteristics. Microb Pathog. 2021:105212.

38. Jin Q, Jung JY, Kim YJ, Eom H-J, Kim S-Y, Kim T-J, et al. Production of L-lactate in Leuconostoc citreum via heterologous expression of L-lactate dehydrogenase gene. J Biotechnol. 2009;144(2):160–4.

39. Levitt MD, Levitt DG. Quantitative evaluation of D-lactate pathophysiology: new insights into the mechanisms involved and the many areas in need of further Investigation. Clin Exp Gastroenterol. 2020;13:321. pmid:32982363

40. Monroe GR, van Eerde AM, Tessadori F, Duran KJ, Savelberg SM, van Alfen JC, et al. Identification of human D-lactate dehydrogenase deficiency. Nat Commun. 2019;10(1):1–8.

41. Bianchetti DG, Amelio GS, Lava SA, Bianchetti MG, Simonetti GD, Agostoni C, et al. D-lactic acidosis in humans: systematic literature review. Pediatr Nephrol. 2018;33(4):673–681. pmid:29218437

42. Bustos AY, de Valdez GF, Fadda S, Taranto MP. New insights into bacterial bile resistance mechanisms: the role of bile salt hydrolase and its impact on human health. Food Res Int. 2018;112:250–62. pmid:30131136

43. Bachmann V, Kostiuk B, Unterweger D, Diaz-Satizabal L, Ogg S, Pukatzki S. Bile salts modulate the mucin-activated type VI secretion system of pandemic Vibrio cholerae. PLoS Negl Trop Dis. 2015;9(8):e0004031. pmid:26317760

44. Ramakrishna BS. FAO/WHO guidelines to assess probiotic efficacy for human consumption–would they suffice? In: Nair GB, Takeda Y, editors. ECAB health impact of probiotics: vision & opportunities-E-Book. 2014;105.

45. Karovičová J, Kohajdová Z, Šimko P, Lukáčová D. Using capillary isotachophoresis for the determination of biogenic amines and D‐isocitric acid in food products. Food/Nahrung. 2003;47(3), 188–190. pmid:12866622

46. Cao C, Wang J, Liu Y, Kwok LY, Zhang H, Zhang W. Adaptation of Lactobacillus plantarum to ampicillin involves mechanisms that maintain protein homeostasis. mSystems. 2020;5(1):e00853–19. pmid:31992633

47. Li T, Teng D, Mao R, Hao Y, Wang X, Wang J. A critical review of antibiotic resistance in probiotic bacteria. Food Res Int. 2020:109571. pmid:32846610

48. Zhang S, Oh JH, Alexander LM, Özçam M, van Pijkeren JP. D-Alanyl- D-alanine ligase as a broad-host-range counterselection marker in vancomycin-resistant lactic acid bacteria. J Bacteriol. 2018;200(13), e00607–17. pmid:29686137

49. Jose B, Antoci V Jr, Zeiger AR, Wickstrom E, Hickok NJ. Vancomycin covalently bonded to titanium beads kills Staphylococcus aureus. Chem Biol. 2005;12(9):1041–8.

50. Kleerebezem M, Boekhorst J, van Kranenburg R, Molenaar D, Kuipers OP, Leer R, et al. Complete genome sequence of Lactobacillus plantarum WCFS1. PNAS. 2003;100(4):1990–5. pmid:12566566

51. Fu B, Huang X, Ma J, Chen Q, Zhang Q, Yu P. Characterization of an inositol-producing Lactobacillus plantarum strain and the assessment of its probiotic potential and antibacterial activity. LWT. 2022;153:112553.

52. Stasic AJ, Chasen NM, Dykes EJ, Vella SA, Asady B, Starai VJ, et al. The toxoplasma vacuolar H+-ATPase regulates intracellular pH and impacts the maturation of essential secretory proteins. Cell Rep. 2019;27(7):2132–46. e7. pmid:31091451

53. Šeme H, Gjuračić K, Kos B, Fujs Š, Štempelj M, Petković H, et al. Acid resistance and response to pH-induced stress in two Lactobacillus plantarum strains with probiotic potential. Benef Microbes. 2015;6(3):369–79. pmid:25380802

54. Wang G, Yu H, Feng X, Tang H, Xiong Z, Xia Y, et al. Specific bile salt hydrolase genes in Lactobacillus plantarum AR113 and relationship with bile salt resistance. LWT. 2021;145:111208.

55. Begley M, Hill C, Gahan CG. Bile salt hydrolase activity in probiotics. Appl Environ Microbiol. 2006;72(3):1729–38. pmid:16517616

56. Burckhardt RM, Escalante-Semerena JC. Small-molecule acetylation by GCN5-related N-acetyltransferases in bacteria. Microbiol Mol Biol Rev. 2020;84(2):e00090–19. pmid:32295819

57. Cai G, Wu D, Li X, Lu J. Levan from Bacillus amyloliquefaciens JN4 acts as a prebiotic for enhancing the intestinal adhesion capacity of Lactobacillus reuteri JN101. Int J Biol Macromol. 2020;146:482–7.

58. Glenting J, Beck HC, Vrang A, Riemann H, Ravn P, Hansen AM, et al. Anchorless surface associated glycolytic enzymes from Lactobacillus plantarum 299v bind to epithelial cells and extracellular matrix proteins. Microbiol Res. 2013;168(5):245–53.

59. Aricescu AR, Hon WC, Siebold C, Lu W, van der Merwe PA, Jones EY. Molecular analysis of receptor protein tyrosine phosphatase μ‐mediated cell adhesion. EMBO J. 2006;25(4):701–12.

60. Ito H, Iwabuchi M, Ogawa KI. The sugar-metabolic enzymes aldolase and triose-phosphate isomerase are targets of glutathionylation in Arabidopsis thaliana: detection using biotinylated glutathione. Plant Cell Physiol. 2003;44(7), 655–660.

61. Feng T, Wang J. Oxidative stress tolerance and antioxidant capacity of lactic acid bacteria as probiotic: A systematic review. Gut Microb. 2020;12(1), 1801944. pmid:32795116

62. Lu J, Holmgren A. The thioredoxin antioxidant system. Free Radic Biol Med. 2014;66, 75–87. pmid:23899494