Front. Microbiol.Frontiers in MicrobiologyFront. Microbiol.1664-302XFrontiers Media S.A.10.3389/fmicb.2024.1347830MicrobiologyOriginal ResearchUnraveling the antimicrobial potential of Lactiplantibacillus plantarum strains TE0907 and TE1809 sourced from Bufo gargarizans: advancing the frontier of probiotic-based therapeuticsHuangFeiyun1ZhaoYanni1HouYusen2YangYu1YueBisong1ZhangXiuyue3*1Key Laboratory of Bio-Resources and Eco-Environment (Ministry of Education), College of Life Sciences, Sichuan University, Chengdu, China2School of Pharmacy, Chengdu University of Traditional Chinese Medicine, Chengdu, China3Sichuan Key Laboratory of Conservation Biology on Endangered Wildlife, College of Life Sciences, Sichuan University, Chengdu, China
Edited by: Yiannis Kourkoutas, Democritus University of Thrace, Greece
Reviewed by: Eugenia Bezirtzoglou, Democritus University of Thrace, Greece
Gregoria Mitropoulou, Democritus University of Thrace, Greece
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Introduction
In an era increasingly defined by the challenge of antibiotic resistance, this study offers groundbreaking insights into the antibacterial properties of two distinct Lactiplantibacillus plantarum strains, TE0907 and TE1809, hailing from the unique ecosystem of Bufo gargarizans. It uniquely focuses on elucidating the intricate components and mechanisms that empower these strains with their notable antibacterial capabilities.
Methods
The research employs a multi-omics approach, including agar diffusion tests to assess antibacterial efficacy and adhesion assays with HT-29 cells to understand the preliminary mechanisms. Additionally, gas chromatography-mass spectrometry (GC-MS) is employed to analyze the production of organic acids, notably acetic acid, and whole-genome sequencing is utilized to identify genes linked to the biosynthesis of antibiotics and bacteriocin-coding domains.
Results
The comparative analysis highlighted the exceptional antibacterial efficacy of strains TE0907 and TE1809, with mean inhibitory zones measured at 14.97 and 15.98 mm, respectively. A pivotal discovery was the significant synthesis of acetic acid in both strains, demonstrated by a robust correlation coefficient (cor ≥ 0.943), linking its abundance to their antimicrobial efficiency. Genomic exploration uncovered a diverse range of elements involved in the biosynthesis of antibiotics similar to tetracycline and vancomycin and potential regions encoding bacteriocins, including Enterolysin and Plantaricin.
Conclusion
This research illuminates the remarkable antibacterial efficacy and mechanisms intrinsic to L. plantarum strains TE0907 and TE1809, sourced from B. gargarizans. The findings underscore the strains' extensive biochemical and enzymatic armamentarium, offering valuable insights into their role in antagonizing enteric pathogens. These results lay down a comprehensive analytical foundation for the potential clinical deployment of these strains in safeguarding animal gut health, thereby enriching our understanding of the role of probiotic bacteria in the realm of antimicrobial interventions.
Gastrointestinal (GI) disorders, predominantly triggered by bacterial pathogens like Pseudomonas aeruginosa, Campylobacter jejuni, Escherichia coli, critically impact global health, with their management complicated by escalating antibiotic resistance worldwide, and it is urgent to develop new green antibiotic products (Greenwood-Van Meerveld et al., 2017; Zhao et al., 2017; Yang et al., 2022). Probiotic lactic acid bacteria (LAB) emerge as potential antibiotic alternatives, effectively inhibiting enteric pathogens for both prevention and treatment. Nonetheless, given the complex and varied composition of gastrointestinal microbiota in humans and animals, coupled with the diversity of probiotic strains and their individualized effects, the study of the interaction mechanisms between probiotics and other microbial communities is still in its nascent phase (Gomaa, 2020). Therefore, there is a critical need for comprehensive research to uncover the inhibitory constituents of probiotics, aimed at facilitating strategic and efficacious colonization of beneficial elements in human and animal gastrointestinal tracts, crucial for maintaining gut microbiota balance and suppressing the growth of intestinal pathogens.
Lactobacillus, a genus within the lactic acid bacteria (LAB) family, is a cornerstone of probiotic research, acclaimed for its vital role in fermentation and significant health benefits, such as enhancing gut flora and boosting immunity (Ballini et al., 2023). With over 220 species, it is the most expansive and diverse LAB group (http://www.bacterio.net/lactobacillus.html). Research shows Lactobacillus species exhibit varied genetic and physiological traits across different environments. Some, like Lactobacillus delbrueckii and Lacticaseibacillus rhamnosus, occupy specialized niches, while Lactiplantibacillus plantarum is remarkably versatile, thriving in diverse habitats from dairy to human tracts, reflecting its ability to adapt metabolically (Valan Arasu et al., 2015; Berbegal et al., 2016; Min Hsiu et al., 2016; Mendoza et al., 2023; Rebaza-Cardenas et al., 2023). Moreover, Endorsed by the European Food Safety Authority (EFSA) and the U.S. Food and Drug Administration (FDA) for its safety, L. plantarum is a versatile probiotic, widely employed in fermenting diverse foods and drinks, and in medical and pharmaceutical applications (Todorov and Franco, 2010; Capozzi et al., 2012). Its health-promoting and antibacterial qualities have not only led to innovative probiotic products but also enhanced food safety through biopreservation (Tang et al., 2023).
Numerous studies have focused on investigating the potential impact of L. plantarum in managing diverse gastrointestinal disorders, including infectious diarrhea, antibiotic-associated diarrhea, irritable bowel syndrome (IBS), and inflammatory bowel disease (Ji et al., 2023; Kumar et al., 2023; Yang et al., 2023). L. plantarum 299v ameliorated Escherichia coli-induced intestinal permeability in rats, mitigating a 53% surge in mannitol permeation through the intestinal barrier, as gauged using Ussing chambers (Mangell et al., 2002). Similarly, Yang et al. (2014) revealed that L. plantarum CGMCC 1258 effectively combats enterotoxigenic E. coli-induced diarrhea in young pigs by bolstering intestinal integrity and regulating tight junction proteins (Yang et al., 2014). Nevertheless, the current body of research on the probiotic mechanism of L. plantarum still needs to be improved due to existing technical limitations. Despite its notable colonization capabilities and ability to effectively target various pathogens, further investigation is required to comprehend its potential fully.
Currently, the integration of multi-omics sequencing methodologies, proficiently addresses the constraints inherent to in vitro assays of LAB, enabling a profound elucidation of its antimicrobial mechanisms and salutary impacts. The availability of complete genome sequences derived from L. plantarum DNA strains permitted us to understand better this LAB species' potential for adaptation and functionality. The complete genome sequencing of L. plantarum strains, starting with the human saliva strain L. plantarum WCFS1, has significantly deepened insights into LAB's adaptability and functionality (Kleerebezem et al., 2003). This genomic revelation, followed by a series of sequenced genomes, has facilitated a comprehensive understanding of LAB's evolutionary journey and environmental interactions through comparative genomics. Siezen et al.'s (2010) study elegantly dissected the genomic landscape of L. plantarum across 42 strains, shedding light on its exceptional adaptability. They identified a core set of 2,049 genes common to all strains, with 121 unique genes distinguishing L. plantarum from other LAB species. Comparative analysis revealed gene variations compared to the WCFS1 genome, with many distinctive genes located in lifestyle adaptation regions, highlighting the bacterium's evolutionary ingenuity (Siezen et al., 2010). Furthermore, in metabolomics, the integration of gas chromatography, gas chromatography-mass spectrometry, and high-performance liquid chromatography is paramount for the precise quantification of short-chain fatty acids, owing to their exceptional volatility and versatility, finding extensive utility in medical and chemical research domains (De Preter et al., 2009; De Baere et al., 2013; Mulat and Feilberg, 2015; Marques et al., 2023).
Residing in the Chordata phylum and Bufonidae family, Bufo gargarizans, a medicinal amphibian indigenous to China, has garnered significant academic interest for its anti-inflammatory and anti-neoplastic properties, particularly in its secretions and dermal extracts (Wang et al., 2012; Qi et al., 2014). This species, characterized by its rich integumentary and digestive microbiota, stands at the forefront of probiotic-related research. While literature on its probiotic associations remains sparse, studies on amphibian-derived probiotics offer valuable perspectives (Kueneman et al., 2016; Antwis and Harrison, 2018).
Our research team, having sourced specimens of B. gargarizans from the Jiuzhaigou Valley Scenic and Historic Interest Area in Sichuan Province, pioneered the isolation and screening of two prominent L. plantarum strains from their intestinal microbiota. In transcending traditional phenotypic analyses of LAB strains, this study meticulous control experiments, genomic profiling, and metabolomic explorations collectively decode the antibacterial dynamics inherent to L. plantarum originating. Notably, our novel undertakings in genomic visualizations and assessments of short-chain fatty acid synthesis chart unprecedented territory in this field. Collectively, our findings constitute a robust multi-omics testament to the clinical significance of these Lactobacillus strains in elucidating molecular mechanisms and bolstering intestinal fortification.
2 Materials and methods2.1 Experimental methods2.1.1 Experimental strains
We selected the L. plantarum strains TE0907 and TE1809 (Supplementary material S1), stemming from preliminary in vitro screenings within our facility, originally derived from the intestinal milieu of mature B. gargarizans specimens. These initial evaluations vouched for their commendable safety profile, manifesting through the absence of hemolysis, susceptibility to a septet of prevalent antibiotics (Clindamycin, Clarithromycin, Gentamicin, Chloromycetin, Tetracycline, Ampicillin, Ceftriaxone) (Supplementary material S2), and marked gastrointestinal resilience. Currently, these strains have been meticulously conserved at the esteemed China General Microbiological Culture Collection Center (CGMCC), with their corresponding study sequences diligently archived in the Sequence Read Archive (SRA) of the National Center for Biotechnology Information (NCBI) website (https://www.ncbi.nlm.nih.gov/), accessible under the BioProject designation PRJNA1066612. The reference strain for the experiment was L. rhamnosus GG ATCC 53103 (LGG). This strain was graciously provided by Yijia Biotechnology Co. Ltd., Chengdu.
2.1.2 Determination of bacterial strains' cell adhesion
HT-29 cells, post-subculturing, were standardized to 5 × 106 cells/mL. In each culture plate well, 1 mL of this cell suspension was seeded and incubated at 37°C in a 5% CO2 environment until monolayer formation. Following DMEM medium removal and a PBS (Aoke, Chengdu) wash, 1 mL bacterial suspension was inoculated into each well in triplicate for every sample. After 2 h at the same conditions, residual bacterial suspensions were discarded. Wells were then thoroughly rinsed with PBS to eliminate unbound bacteria. Cells were detached utilizing 0.7 mL of 0.25% trypsin-EDTA (Aoke, Chengdu) for 10 min, halted by adding 0.3 mL DMEM. Subsequent cell suspensions were quantified using a 10-fold dilution method. Adhesion efficiency was assessed by comparing bacterial concentrations pre- and post-incubation.
2.1.3 Pathogen susceptibility assay
The Oxford cup assay evaluated the antagonistic potential against various pathogens. Bacterial strains of interest, namely E. coli CMCC (B) 44102, S. aureus CMCC (B) 26003, P. aeruginosa CMCC (B), C. jejuni ATCC 33291, C. coli ATCC 43478, S. typhimurium ATCC 14028, S. haemolyticus ATCC29970, Y. enterocolitica ATCC 23715, P. multocida ATCC51689, L. monocytogenes ATCC 19115, and Salmonella paratyphi type B CMCC (B) 50094, were propagated in LB medium (Huankai, Guangdong) until the exponential phase. Post-cultivation, these strains were collected, rinsed with PBS, and adjusted to a 108 CFU/mL final concentration. Petri dishes were prepared with 20 mL of LB agar (Huankai, Guangdong) enriched with 1.5% agarose for the assay. Upon solidification, three symmetrically spaced Oxford cups were anchored onto each plate, dispensing 200 μL of TE0907 and TE1809′s HMV18 cell-free supernatant (CFS). L. rhamnosus GG ATCC 53103 (LGG) was employed as a reference strain for comparative control. The inhibitory zone diameters of each strain were standardized using the Z-Score normalization technique.
2.1.4 Exploration of antibacterial active ingredients
In order to further understand the effective ingredients involved in the excellent antibacterial properties of the two probiotics mentioned above, the conditioned culture supernatant (CFS) of the probiotics was prepared again. Utilizing a cohort of five meticulously selected animal pathogens as indicative models, the prepared cell-free supernatant (CFS) was subjected to a series of distinct, non-sequential interventions designed to elucidate the spectrum of antibacterial components. These interventions comprised the assay of hydrogen peroxide via specialized test strips (Huankai, Guangdong), thermal treatment in a 70°C precisely controlled water bath (Yiheng, Shanghai), pH normalization to 7 using a 1.0 M NaOH solution (Jinbohua, Tianjin), and the incorporation of proteinase K (Suolaibao, Beijing) into the CFS. This strategic approach was employed to dissect the antibacterial constituent profile inherent within these strains.
2.1.5 GC-MS analysis of acid production capacity
We utilized a Shimadzu GCMS-QP2010 Plus gas chromatograph-mass spectrometer with a RESTEK Rtx-5 column (30 m × 0.25 mm × 0.25 μm). The temperature program started at 40°C for 5 min, increased to 150°C at 5°C/min, then to 280°C at 10°C/min, holding for 2 min. High-purity helium (>99.999%) was the carrier gas at 1.0 mL/min. Mass spectrometry parameters included electron impact ionization at 200°C, an interface temperature of 220°C, and a m/z 33–500 scan range. For sample preparation, fermentation broth (4 mL) was mixed with 10 μL of a 2-ethylbutyric acid internal standard (200 μg/mL). Samples were injected in split mode (1:3 ratio, 1 μL volume) at 270 μL with a 0.1-min solvent delay. This method quantified acids such as acetic, butyric, and valeric.
Furthermore, Pearson correlation analysis was employed to investigate the association between the concentrations of various short-chain fatty acids in separate bacterial strains and their average inhibitory effects against 11 tested pathogenic bacteria.
2.1.6 Genome sequencing, genome annotation and gene annotation related to bacteriostasis
The genomes of L. plantarum TE0907 and L. plantarum TE1809 were sequenced using the HiSeqTM 2000 platform by Biotechnology Corporation in Shanghai, China. To ensure precise annotation information for the predicted protein-encoding genes, multiple databases, such as NR (Altschul et al., 1990), Swiss-Prot, GO, KEGG, and COG (Tatusov et al., 2003), were employed for Blast alignment. Furthermore, Circos (Krzywinski et al., 2009) generated a circular genome chromosome map that visualizes chromosome sequences, gene positions, mutation information, and other relevant data. Lastly, the BAGEL4 database (van Heel et al., 2018) was employed for predictive alignment of bacteriocin-related genes within the genome of LAB.
2.1.7 Statistical analysis
Statistical analyses were performed using GraphPad Prism version 8 (California, USA), presenting data as means ± standard deviation (SD) and employing one-way ANOVA with the Student's t-test for comparisons. Significance was established at P < 0.05; all experiments were independently conducted and iterated at least thrice.
3 Results3.1 Evaluation of adhesion rate of intestinal epidermis
TE0907 (55.02 ± 3.18%) exhibits no significant disparity when contrasted with LGG (55.71 ± 41.29%) (P > 0.05), while TE1809 (85.71 ± 3.70%) significantly surpasses LGG (Figure 1). Hence, at this preliminary stage, it can be inferred that the mentioned strains possess notable adhesion capabilities, implying their potential to colonize the intestinal epithelium, forming the basis for their prospective antibacterial functions.
Adhesion rates of LGG, TE0907, and TE1809 to HT-29 cells.
3.2 Antibacterial activity assessment
Prominent disparities were observed in the inhibitory prowess displayed by distinct LAB strains. In conjunction with the LGG strain, the two candidate strains manifested antagonistic capacities against pathogenic bacteria. Notably, a preponderance of these probiotics unveiled heightened antimicrobial efficacy against P. multocida, yet presented a subdued inhibitory response to S. aureus (Figure 2A). Contrastingly, TE0907(average value: 14.97 mm) and TE1809 (15.98 mm) outstripped LGG (12.56 mm) in their antimicrobial potency against an eclectic assortment of animal pathogens, especially underscoring their efficacy against gram-positive bacteria, thereby elucidating their broad-spectrum inhibitory attributes (Figure 2B). Intriguingly, our two strains of L. plantarum, when juxtaposed with the L. plantarum DY1 studied by Mao et al. (2020) (E. coli: 12.89 ± 0.21 mm; Salmonella: 15.93 ± 0.22 mm; S. aureus: 15.70 ± 0.41 mm), exhibited an equally compelling or even superior antibacterial zone diameter using a comparable method for antimicrobial performance testing. This denotes an unequivocal competitive edge in antimicrobial efficacy within the same species.
Antibacterial activity assessment. (A) LGG and the inhibitory diameter distribution chart of L. plantarum TE0907 and L. plantarum TE1809; (B) statistical diagram for bacteriostatic test of LGG and L. plantarum TE0907 and L. plantarum TE1809.
3.3 Exploration of antibacterial active ingredients
CFS showed no evidence of hydrogen peroxide production after evaluating hydrogen peroxide test strips. However, the antibacterial potency of CFS from both TE0907 and TE1809 was lost after pH adjustment, underscoring the primary role of organic acids in its antibacterial effect. Post proteinase K treatment, there was a marginal decline in the antibacterial activity of the strains, suggesting that L. plantarum might produce antibacterial peptides or proteins that play a minor role in their inhibitory action. Notably, following heat exposure to the CFS, the antibacterial efficacy of both strains remained largely intact against most tested bacteria (Table 1).
Statistics of active antibacterial ingredients (a: P < 0.05; b: P < 0.01; c: P < 0.001; d: P < 0.0001).
Strains
Treatment
S. aureus (mm)
P. aeruginosa (mm)
S. haemolyticus (mm)
Salmonella paratyphi type B (mm)
E. coli (mm)
L. rhamnosus GG
Blank
0.00
0.00
0.00
0.00
0.00
CFS
9.20 ± 0.04
14.42 ± 0.20
12.67 ± 0.09
13.45 ± 0.70
13.63 ± 0.77
pH 7.0
0.00d
0.00d
0.00d
0.00d
0.00d
70°C
9.38 ± 0.13
13.10 ± 0.02d
12.25 ± 0.36
12.27 ± 0.22a
12.17 ± 0.23a
Proteinase K
8.22 ± 0.01d
11.42 ± 0.01d
10.38 ± 0.05d
11.33 ± 0.21b
10.36 ± 0.06c
L. plantarum TE0907
Blank
0.00
0.00
0.00
0.00
0.00
CFS
13.75 ± 0.41
16.18 ± 0.39
15.23 ± 0.21
16.27 ± 1.16
13.92 ± 0.23
pH 7.0
0.00d
0.00d
0.00d
0.00d
0.00d
70°C
13.37 ± 0.25
14.56 ± 0.04c
13.94 ± 0.4b
15.94 ± 0.17
13.21 ± 0.01
Proteinase K
11.21 ± 0.03d
12.55 ± 0.27d
12.33 ± 0.28d
14.09 ± 0.05a
10.85 ± 0.37c
L. plantarum TE1809
Blank
0.00
0.00
0.00
0.00
0.00
CFS
13.36 ± 0.28
17.06 ± 0.18
13.80 ± 0.95
14.29 ± 0.11
13.99 ± 0.77
pH 7.0
0.00d
0.00d
0.00d
0.00d
0.00d
70°C
13.51 ± 0.38
15.75 ± 0.07d
13.73 ± 0.16
14.37 ± 0.52
13.99 ± 0.12
Proteinase K
9.26 ± 0.16d
12.42 ± 0.12d
12.39 ± 0.23a
9.78 ± 0.79c
12.32 ± 0.1b
3.4 Acid production capacity and acid production component analysis
GC-MS analysis identified acetic acid (3.576–13.08 μg/mL), N-butyric acid (0.018–0.052 μg/mL), isobutyric acid (0.06–0.18 μg/mL), isovaleric acid (0.061–0.41 μg/mL), and isocaproic acid 2-Methylbutyric acid (0.022–0.23 μg/mL) as the salient SCFAs, a revelation corroborated by Figure 3A total ion chromatogram, which adeptly profiles these SCFAs with commendable separation integrity, negligible tailing, and a harmonized baseline. The juxtaposition of the internal standard, 2-methylbutyric acid, with other pivotal peaks underscores the assay's unparalleled separation prowess, buttressing its robustness, especially in isolating SCFAs, a claim bolstered by the conspicuous absence of peak interference (Figures 3C, D). It is also paramount to note that acetic acid dominates the SCFAs landscape across all strains and exhibits a marked affinity with the levels of other acids. Furthermore, our study uncovered a robust correlation (correlation coefficient ≥ 0.943) between the antimicrobial effectiveness and the presence of acetic acid, a connection far more distinct compared to other short-chain fatty acids (SCFAs), which was significantly relevant in inhibiting S. haemolyticus and Y. enterocolitica (P < 0.05) (Figure 3B). Strains proficient in metabolizing robust concentrations of acetic acid may engender a pronouncedly acidic milieu, detrimental to the proliferation of pathogenic organisms, thereby effectuating a competitive exclusion dynamic.
Study on the production of SCFAs by LGG and 2 LAB strains using GC-MS. (A) SCFAs content of each strain; (B) The correlation between the short-chain fatty acid content of each strain and the inhibitory effects on certain pathogens was analyzed using the R package ggcor (SH, S. haemolyticus; LM, L. monocytogenes; YE, Y. enterocolitica); (C, D) Total ion flow chromatogram for the determination of SCFAs in L. plantarum TE0907 and L. plantarum TE1809 using GC-MS.
3.5 Genomic analysis
After second-generation Illumina HiSeq sequencing and quality control, over 85% of the bases in the genomes of two bacterial strains were of high quality, ensuring reliable subsequent analysis. TE0907′s genome spans 3,154,032 base pairs with a 44.69% GC content and 3,055 protein-coding genes, while TE1809′s genome measures 3,294,124 base pairs, a 44.32% GC content and houses 3,227 protein-coding genes. The sequencing data is diligently archived in the Sequence Read Archive (SRA) of the National Center for Biotechnology Information (NCBI) website (https://www.ncbi.nlm.nih.gov/), and is accessible under the BioProject designation PRJNA1066612.
Both LAB strains exhibited extensive functional gene annotations, with 99.79% (2,910 genes) and 99.55% (3,067 genes) aligning to the NR database, predominantly corresponding to the L. plantarum species, consistent with prior 16S rRNA gene-based identifications (Figure 4A; Supplementary Figure S1A). We assembled the genomes of two bacterial strains and constructed their chromosomal circos plots (Figure 4B; Supplementary Figures S1B). The GO database identified 2,108 and 2,125 genes in these strains, mainly in the molecular function domain, highlighting their advanced biosynthetic aptitude. KEGG annotations revealed their affinity to metabolic pathways, with TE0907 and TE1809 having specific gene counts for glucose and amino acid metabolism and a notable focus on secondary metabolite production. COG analysis emphasized their protein synthesis prowess, with over 150 genes related to translation and ribosomal activities (Figures 4C–E; Supplementary Figures S1C–E). Additionally, the annotation outcomes of bacteriocin and immune-related genes in the KEGG database revealed that the genomes of the two strains exhibited the presence of genes associated with tetracycline biosynthesis (ko00253), vancomycin antibiotic synthesis (ko01055), and cytochrome P450 metabolism of xenobiotics (ko00980), which are known to modulate diseases relevant to humans (Table 2).
Genome functional annotation of L. plantarum TE0907. (A) Chromosome Circos plot of the L. plantarum TE0907 genome [annotations for the positive and negative strands are represented on the circular diagram's outer and inner sides, respectively. The outermost layer marks chromosomal DNA base positions. Layers two and three show protein-coding genes, while layers four to seven represent COG and KEGG functional annotations. Layers eight and nine cover GO annotations, and layers ten and eleven contain rRNA (5S, 16S, 23S), tRNA, and sRNA gene information. The twelfth layer displays relative GC content, and the thirteenth indicates GC skew; Color annotations are provided in Supplementary material S3; (B–E) NR, COG, GO, KEGG database annotation of the L. plantarum TE0907 genome based on NCBI Blast + alignment.
Bacteriocin, immunoregulation and related genes in the genomes of TE0907 and TE1809.
The bacteriocin gene cluster known as Plantaricin was identified in TE0907 and TE1809, respectively. TE0907 exhibited six distinct gene clusters: Plantaricin_E, Plantaricin_F, Plantaricin_J, Plantaricin_K, Plantaricin_NC8-alpha, and Plantaricin_NC8-beta (Figure 5A). Conversely, the L. plantarum TE1809 genome unveiled two salient regions: the first encompassing gene clusters Plantaricin_E, Plantaricin_F, Plantaricin_N, and Plantaricin_J, and the second highlighting the Enterolysin_A gene cluster (Figures 5B, C). Multiple open reading frames (ORF) within these bacteriocins underscore a potential correlation between the observed antibacterial characteristics in both strains and the bacteriocins' existence. Employing a genomic lens, we have corroborated the efficacious antibacterial components metabolized by the strain, as initially discerned in our antecedent controlled variable investigations.
Metabolic system analysis of TE0907 and TE1809. Annotation of bacteriocins in the L. plantarum TE0907 and L. plantarum TE1809 genome. (A) Alignment-based prediction of bacteriocin-related genes encoded within the L. plantarum TE0907 genome using the BAGEL4 database; (B, C) alignment-based prediction of bacteriocin-related genes encoded within the L. plantarum TE1809; (D, E) annotation of the L. plantarum TE0907 and L. plantarum TE1809 genome CAZy database alignment based on HMMER3 (E-value < 1e-5); (F) heat map of glycoside hydrolase (GH) distribution in the genomes of L. plantarum TE0907 and L. plantarum TE1809.
The two strains manifest distinct quantities of genes for carbohydrate-active enzymes (CAZymes), with 136 and 133 identified in each, and notably feature a substantial presence of GHs genes, at 48 (35.3%) and 49 (36.1%), respectively, followed by genes for glycosyl transferases (GTs) at 35 (25.7%) and 34 (25.6%), and genes for carbohydrate esterases (CEs) at 25 (18.4%) and 24 (18.0%), respectively, illustrating a diverse enzymatic composition. Furthermore, it is worth noting that only the genome of L. plantarum TE0907 has been subjected to gene annotation, explicitly focusing on identifying genes encoding polysaccharide lyases (PLs) (Figures 5D, E). Heatmaps were constructed using the top 10 GHs in terms of abundance from two selected genomes to provide a more intuitive representation of GH distribution. The results indicated that GH1, GH5, and GH13 have a relatively higher prevalence in the genomes of each bacterial strain (Figure 5F). Simultaneously, this revelation underscores the potent metabolic and biosynthetic prowess of the L. plantarum strains, thereby elucidating the molecular physiological underpinnings pivotal for their pronounced antibacterial potency.
Melting curve analysis of the genomes of strains TE0907 and TE1809 showcased a coordinated reduction in shared genes and a rise in non-redundant genes (Figure 6A; Supplementary Figure S2A), hinting at an emerging open pan-genome in the L. plantarum genomes under investigation. Notably, both strains consistently shared 1,367 and 1,499 genes, respectively, aligning with the gene profiles of all 15 reference L. plantarum genomes. Furthermore, with 100 and 153 unique genes in their genomes, they surpassed the gene counts in most other reference strains (100 vs. 9–150; 153 vs. 36–115) (Figure 6B; Supplementary Figure S2B). While 1,367 and 1,499 genes were commonly shared across various L. plantarum strains, the distribution of single-copy and multi-copy homologous genes remained consistent: for TE0907, 1267–1285 single-copy (1275 in TE0907) and 72–100 multi-copy (92 in TE0907); for TE1809, 1267–1295 single-copy and 72–100 multi-copy (92 in TE1809) (Figure 6C; Supplementary Figure S2C). These genes chiefly involve cell wall/membrane biogenesis, replication, recombination, repair, and general function prediction (Figures 6E, F; Supplementary Figure S2F). Phylogenetically, amino acid sequences from homologous genes place TE0907 close to plant-associated L. plantarum JDM1, while TE1809 aligns closely with L. plantarum subsp. plantarum ST.III (Figure 6D; Supplementary Figure S1D).
Pan-genome analysis of L. plantarum TE0907. (A) Box plot of Core-Pan gene dilution (green represents the Pan gene dilution curve, red represents the Core gene dilution curve. The x-axis represents the number of samples selected for each statistical analysis, and the y-axis represents the distribution of the number of samples); (B) homologous gene cluster petal diagram (showing the number of shared and unique orthologs clusters); (C) bar chart of homologous gene number statistics [single-copy: number of single-copy homologous genes, multiple-copy: number of multiple-copy homologous genes; core: genes that are present in all species; dispensable: genes that are present in some (more than one) species; specific: genes that are present in only one species]; (D) a phylogenetic tree diagram constructed using the Neighbor-Joining clustering method based on the pan-genome; (E, F) bar chart of COG annotations for shared and unique orthologs clusters.
4 Discussion
L. plantarum strains with superior probiotic properties play a pivotal role in maintaining gut microbial homeostasis, offering targeted treatment against pathogens, and mitigating the overuse of antibiotics (Gurunathan et al., 2023). In this research, through a combination of bacterial genomic sequencing and GC-MS spectrometry analysis, we, for the first time, unveiled the mechanistic actions and associated functional genes of two L. plantarum strains derived from B. gargarizans in terms of their antimicrobial activity and their capability to modulate the gut microbial composition. This investigation furnishes a more holistic and precise approach to evaluating the efficacy of probiotics, providing novel insights into the field.
The colonization proficiency of gastrointestinal tract (GIT) strains, crucial for their probiotic efficacy, involves not only adherence to mucus and epithelial cells for a competitive advantage, but also the cultivation of substantial interactions with the host, a process essential for probiotic recognition and immune response stimulation, thus significantly magnifying their beneficial effects compared to strains that merely pass through the GIT (Bangotra et al., 2023; Dahiya and Nigam, 2023; Han et al., 2023). In addition, the colonization process by probiotic strains serves as a preventive measure against the attachment of enteric pathogens to the cells of the digestive tract (Son et al., 2023). Extensive research has demonstrated that probiotics exert inhibitory effects on the adhesion of gastrointestinal pathogens, including S. typhimurium, Clostridium sporogenes, and Enterococcus faecalis, to Caco-2 cells (Bachtarzi et al., 2024; Jain et al., 2023; Lyu et al., 2023). Such adhesion is a prerequisite for the beneficial role of LAB in countering pathogenic bacteria, balancing the intestinal microbiota, and modulating immune responses (Grigoryan et al., 2018; Van Zyl et al., 2020). A key factor in this adhesion process is the probiotic mucin-binding proteins, with the MUB gene playing a central role in adherence to mucin proteins and being a ubiquitous feature in probiotic genomes (Singh et al., 2018). This gene-centric perspective is instrumental in stratifying probiotic strains, emphasizing those with higher expression of MUB gene proteins. Our experiments have highlighted the exceptional adhesion rates of TE0907 and TE1809 to the HT-29 cell line, aligning with existing literature and indicating a fundamental requirement for effective antibacterial action. The distinctive properties of L. plantarum, including its strong adhesion to intestinal epithelial cells, ability to persistently colonize, and the antagonistic effects of its metabolites on pathogenic bacteria, solidify its standing as a prototypical probiotic (Dell'Anno et al., 2021; Jeng and Yan, 2022).
Subsequent to meeting the prerequisite of intestinal epithelial adherence, the endeavor progresses toward identifying strains with antimicrobial properties, wherein the assessment of LAB's antibacterial activity typically involves utilizing corresponding pathogenic bacteria as indicator organisms to determine the strains' effectiveness. Our selection was based on the antimicrobial capacity of the isolates against E. coli, S. aureus, and C. jejuni, etc., which have been linked to numerous foodborne outbreaks (Gray and Killinger, 1966; Lowy, 1998; Kerr and Snelling, 2009; Gomes et al., 2016; Dar et al., 2017), with criteria emphasizing isolates that exhibited a zone of inhibition surpassing that of a control strain, specifically a commercially recognized strain known for producing antimicrobial metabolites (Ahmed et al., 2023). Our goal was to find new potential antimicrobial metabolites-producing strains and ones that could have a more significant potential than cells already reported, when compared to the L. plantarum DY1 as reported by Mao et al. (2020), our strains not only exhibit enhanced antibacterial activity but also potentially possess a broader spectrum of microbial inhibition. Probiotic bacteria can render the intestinal environment less conducive to developing pathogens, serving as an additional mechanism for pathogen displacement within the gut. For example, LAB suppress pathogenic bacteria by producing organic acids, antimicrobial peptides, and hydrogen peroxide through metabolism, thereby altering the microbial structure in the intestines and protecting the host's health. Lactic acid in its undissociated state acts as a permeabilizer for the outer cell membrane of Gram-negative bacteria. Once inside the bacterial cytoplasm, it dissociates (Zhang et al., 2023). This dissociation leads to a decrease in intracellular pH level due to the accumulation of ionized forms of organic acids and other antimicrobial substances, ultimately resulting in a bactericidal effect (Jyung et al., 2023). Moreover, short-chain fatty acids (SCFAs), which are vital byproducts of the gut microbiome, have diverse roles. They can modify chemotaxis and phagocytosis, induce reactive oxygen species (ROS), regulate cell division and function, and exert anti-inflammatory, antitumorigenic, and antibacterial effects, besides influencing gut integrity (Tan et al., 2014). Among these SCFAs, acetic acid is notable for being the predominant short-chain fatty acid produced by isolated LAB strains in this study. It is not only a critical metabolic product in the fermentation processes of most bacteria but also a primary substrate for cholesterol synthesis. Additionally, it plays a role in the production of butyric acid (Vidra and Németh, 2018). Interestingly, our findings reveal a striking association (cor ≥ 0.943) between antimicrobial efficacy and acetic acid, a correlation significantly more pronounced than those observed with alternative SCFAs. Consequently, we hypothesize that the strain produces an optimal concentration of acetic acid, which is likely to permeate the cell membranes of S. haemolyticus and Y. enterocolitica, subsequently lowering the intracellular pH levels. The acidic environment induced by the acetic acid may lead to protein denaturation, disrupting the physiological activities of the bacteria and resulting in widespread cellular dysfunction. This mechanism suggests a targeted antimicrobial strategy that exploits the disruptive potential of acetic acid on pathogenic bacterial cells (Ewadh et al., 2013; Karabiyikli and Sengun, 2017). This association highlight acetic acid's potential central role in antimicrobial activity, warranting further exploration of its mechanistic contributions within microbial ecology and host-microbe dynamics for devising innovative antimicrobial approaches.
LAB-bacteriocins comprise a heterogeneous group of physicochemically diverse ribosomally-synthesized peptides or proteins showing a narrow or broad antimicrobial activity spectrum against Gram-positive bacteria (Choi et al., 2023; Bertuzzi et al., 2024). Bacteriocins, including nisin derived from L. lactis, plantaricin generated by L. plantarum, and lacticin B obtained from L. acidophilus, have antimicrobial activity against food-borne enteropathogens such as Listeria, Clostridium, Bacillus, methicillin-resistant S. aureus, and vancomycin-resistant enterococci (VRE) (Garg et al., 2014; Mathipa and Thantsha, 2017). Bacteriocins can have either a bacteriostatic or a direct bactericidal impact on pathogens, limiting the cells' ability to establish colonization inside the gastrointestinal tract. The antibacterial properties of bacteriocins enable probiotic microorganisms that produce these compounds to acquire a competitive edge within the intricate gastrointestinal milieu (Oftedal, 2023). We conducted the first genomic sequencing analysis of LAB strains from B. gargarizan. Furthermore, both strains harbor genes linked to antibiotic biosynthesis, cytochrome P450 metabolism of heterocycles, suggesting their antimicrobial and immune-modulating capacities. Through annotation in the BACEL4 database, we discovered six types of plantaricins were detected in L. plantarum TE0907, which have been proven to inhibit pathogenic bacteria such as S. aureus, Salmonella, and E. coli. Among them, Plantaricin_E and Plantaricin_F form a pair of Class IIb bacteriocins, synergistically exhibiting antimicrobial effects (Wang et al., 2022). The bacteriocin paired with Plantaricin_J is Plantaricin_K. The genome of L. plantarum TE1809 contains two potential regions. It is worth noting that these bacteriocins have multiple ORFs, which encode specific immune proteins responsible for bacteriocin secretion, specific proteins involved in modification, and specific proteins involved in transport and leader cleavage. Most ORFs share homology with plnD, plnI, plnH, and plnS, indicating that they may all belong to the pln operon and collectively participate in bacteriocin production (Diep et al., 2009). It is speculated that plantaricin played a significant role in the antimicrobial activity of the selected bacterial strain in this screening. In investigating L. plantarum-derived bacteriocins, a complex array of antimicrobial actions has been revealed, primarily focused on undermining the structural and functional integrity of bacterial membranes. Plantaricin FB-2, for example, disrupts the cell membrane, leading to the leakage of vital cytoplasmic substances and an increase in reactive oxygen species (Li et al., 2023). Other plantaricins weaken the proton motive force or create membrane pores, inducing a bactericidal effect on related microorganisms (Wang et al., 2015; Goel and Halami, 2023). Plantaricin A (PlnA) enhances antibiotic effectiveness by increasing bacterial outer membrane permeability and disrupting membrane integrity through interactions with lipopolysaccharides (Meng et al., 2022). These diverse mechanisms highlight the potential of plantaricins as powerful antimicrobial agents, offering a promising approach to combat pathogens and advance microbial therapies.
Upon meticulous examination of metabolic capabilities, the L. plantarum strains manifest a pronounced functional enrichment in the realms of carbohydrate metabolism and biosynthesis, with a substantial proportion exceeding 8% of protein-coding genes being intricately associated with carbohydrate transport and metabolic processes. Comparative genomics within the same taxonomic clade elucidates that L. fermentum YLF016 and L. pentosus ZFM94 are endowed with 6.00 and 8.28% of their genomic repertoire dedicated to such functionalities, respectively (Ye et al., 2020; Zhang et al., 2022). Notwithstanding the species-specific variance in gene abundance, a juxtaposition of the two strains reveals a comprehensive array of genes implicated in carbohydrate and lactic acid metabolism. This discovery further corroborates prior academic investigations into the probiotic mechanisms attributed to the metabolic derivatives of Lactobacillus strains, reinforcing their role as the fundamental molecular underpinnings for inhibiting pathogens and bolstering host gastrointestinal wellbeing (Sakaridis et al., 2012).
5 Conclusion
We have elucidated, for the first time, the antibacterial efficacy and underlying mechanism of L. plantarum derived from B. gargarizans utilizing multi-omics approaches, enhancing our comprehension of the diverse array of chemicals and enzymes produced by these microorganisms, as a result, the range of potential uses. The conspicuous antimicrobial attributes manifested in strains TE0907 and TE1809 can largely be ascribed to their salient adhesive prowess and metabolic vigor. Of significance, the resultant antimicrobial dynamism is chiefly propelled by the generation of acetic acid coupled with the production of class II and III bacteriocins.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was funded by the National Natural Science Foundation of China (U21A20409).
We sincerely thank Dr. Shengzhi Yang for his invaluable guidance on the experimental methodologies.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2024.1347830/full#supplementary-material
Genome functional annotation of L. plantarum TE1809. (A) Chromosome Circos plot of the L. plantarum TE1809 genome (Annotation as in Figure 4A); (B–E) NR, COG, GO, KEGG database annotation of the L. plantarum TE1809 genome based on NCBI Blast+ alignment.
Pan-genome analysis of L. plantarum TE1809. (A) Box plot of Core-Pan gene dilution (green represents the Pan gene dilution curve, red represents the Core gene dilution curve. The x-axis represents the number of samples selected for each statistical analysis, and the y-axis represents the distribution of the number of samples); (B) homologous gene cluster petal diagram (showing the number of shared and unique orthologs clusters); (C) bar chart of homologous gene number statistics [Single-copy: number of single-copy homologous genes, Multiple-copy: number of multiple-copy homologous genes; Core: genes that are present in all species; Dispensable: genes that are present in some (more than one) species; Specific: genes that are present in only one species]; (D) A phylogenetic tree diagram constructed using the Neighbor-Joining clustering method based on the pan-genome; (E, F) Bar chart of COG annotations for shared and unique orthologs clusters.
Genomic sequences of bacterial strains.
The sensitivity testing results of TE0907 and TE1809 against seven common antibiotics [sensitive (S), intermediate (I), and resistant (R), with antimicrobial disk concentrations in μg per disk and inhibition zone diameter criteria in millimeters] are based on the standards established by the Clinical and Laboratory Standards Institute (CLSI) of the United States.
Genomic chromosomal circos plot annotation file.
ReferencesAhmedN. A.KhattabR. A.RagabY. M.HassanM. (2023). Safety assessment of Enterococcus lactis strains complemented with comparative genomics analysis reveals probiotic and safety characteristics of the entire species. BMC Genom.24:667. 10.1186/s12864-023-09749-937932698AltschulS. F.GishW.MillerW.MyersE. W.LipmanD. J. (1990). Basic local alignment search tool. J. Mol. Biol.215, 403–410. 10.1016/S0022-2836(05)80360-2AntwisR. E.HarrisonX. A. (2018). Probiotic consortia are not uniformly effective against different amphibian chytrid pathogen isolates. Mol. Ecol.27, 577–589. 10.1111/mec.1445629218845BachtarziN.GomriM. A.MeradjiM.Gil-CardosoK.OrtegaN.ChomiciuteG.. (2024). In vitro assessment of biofunctional properties of Lactiplantibacillus plantarum strain Jb21-11 and the characterization of its exopolysaccharide. Int. Microbiol.27, 239–256. 10.1007/s10123-023-00387-537286917BalliniA.CharitosI. A.CantoreS.TopiS.BottalicoL.SantacroceL. (2023). About functional foods: the probiotics and prebiotics state of art. Antibiotics12:635. 10.3390/antibiotics1204063537106999BangotraR.KourB.BajajB. K. (2023). Probiotic functional attributes of lactic acid bacteria from indigenously fermented milk product kalarei. Indian J. Exp. Biol. 61, 7. 10.56042/ijeb.v61i07.3116BerbegalC.PeñaN.RussoP.GriecoF.PardoI.FerrerS.. (2016). Technological properties of Lactobacillus plantarum strains isolated from grape must fermentation. Food Microbiol.57, 187–194. 10.1016/j.fm.2016.03.00227052718BertuzziM. A.GamboniJ. E.SlavutskyA. M.IbargurenC. (2024). “Novel food packaging systems with antimicrobial agents from microbial source,” in Food Packaging and Preservation, eds JaiswalA. K.ShankarS. (Academic Press), 91–111. 10.1016/B978-0-323-90044-7.00006-9CapozziV.RussoP.LaderoV.FernándezM.FioccoD.AAlvarezA.. (2012). Biogenic amines degradation by Lactobacillus plantarum: toward a potential application in wine. Front. Microbiol.3:122. 10.3389/fmicb.2012.0012222485114ChoiG.-H.HolzapfelW. H.TodorovS. D. (2023). Diversity of the bacteriocins, their classification and potential applications in combat of antibiotic resistant and clinically relevant pathogens. Crit. Rev. Microbiol.49, 578–597. 10.1080/1040841X.2022.209022735731254DahiyaD.NigamP. S.-N. (2023). Inclusion of dietary-fibers in nutrition provides prebiotic substrates to probiotics for the synthesis of beneficial metabolites SCFA to sustain gut health minimizing risk of IBS, IBD, CRC. Recent Progr. Nutr.3:2303017. 10.21926/rpn.2303017DarM. A.AhmadS. M.BhatS. A.AhmedR.UrwatU.MumtazP. T.. (2017). Salmonella typhimurium in poultry: a review. Worlds Poult. Sci. J.73, 345–354. 10.1017/S0043933917000204De BaereV.EeckhautS.SteppeM.MaesschalckC. D.BackerP. D.Van ImmerseelF.. (2013). Development of a HPLC–UV method for the quantitative determination of four short-chain fatty acids and lactic acid produced by intestinal bacteria during in vitro fermentation. J. Pharm. Biomed. Anal.80, 107–115. 10.1016/j.jpba.2013.02.03223542733De PreterV.StaeyenG. V.EsserD.RutgeertsP.VerbekeK. (2009). Development of a screening method to determine the pattern of fermentation metabolites in faecal samples using on-line purge-and-trap gas chromatographic–mass spectrometric analysis. J. Chromatogr. A1216, 1476–1483. 10.1016/j.chroma.2008.12.09519167006Dell'AnnoM.GirominiC.ReggiS.CavalleriM.MoscatelliA.OnelliE.. (2021). Evaluation of adhesive characteristics of L. plantarum and L. reuteri Isolated from Weaned Piglets. Microorganisms9:1587. 10.3390/microorganisms908158734442665DiepB. D.StraumeD.KjosM.TorresC.NesI. F. (2009). An overview of the mosaic bacteriocin pln loci from Lactobacillus plantarum. Peptides30, 1562–1574. 10.1016/j.peptides.2009.05.01419465075EwadhM.HasanH.BnyanI.MousaF.SultanJ.EwadhM. (2013). Antibacterial activity of 2-(2-Hydroxy phenylimino) acetic acid. Adv. Life Sci. Technol.7, 15–19.GargN.OmanT. J.WangT.-S. A.GonzaloC. V. G. D.WalkerS.Van Der DonkW. A. (2014). Mode of action and structure–activity relationship studies of geobacillin I. J. Antibiot.67, 133–136. 10.1038/ja.2013.112GoelA.HalamiP. M. (2023). Structural and biosynthetic diversity of plantaricins from. Appl. Microbiol. Biotechnol.107, 5635–5649. 10.1007/s00253-023-12692-037493805GomaaE. Z. (2020). Human gut microbiota/microbiome in health and diseases: a review. Antonie Van Leeuwenhoek113, 2019–2040. 10.1007/s10482-020-01474-7GomesA. T. T.EliasW. P.ScaletskyI. C. A.GuthB. E. C.RodriguesJ. F.PiazzaR. M. F.. (2016). Diarrheagenic Escherichia coli. Braz. J. Microbiol.47, 3–30. 10.1016/j.bjm.2016.10.015GrayM. L.KillingerA. H. (1966). Listeria monocytogenes and listeric infections. Bacteriol. Rev.30, 309–382. 10.1128/br.30.2.309-382.19664956900Greenwood-Van MeerveldB.JohnsonA. C.GrundyD. (2017). Gastrointestinal physiology and function. Gastrointest. Pharmacol.239, 1–16. 10.1007/164_2016_118GrigoryanS.BazukyanI.TrchounianA. (2018). Aggregation and adhesion activity of lactobacilli isolated from fermented products in vitro and in vivo: a potential probiotic strain. Probiot. Antimicrob. Prot.10, 269–276. 10.1007/s12602-017-9283-928455689GurunathanS.ThangarajP.KimJ.-H. (2023). Postbiotics: functional food materials and therapeutic agents for cancer, diabetes, and inflammatory diseases. Foods13:89. 10.3390/foods1301008938201117HanM.LeiW.LiangJ.LiH.HouM.GaoZ. (2023). The single-cell modification strategies for probiotics delivery in inflammatory bowel disease: a review. Carbohydr. Polym.324, 121472. 10.1016/j.carbpol.2023.121472JainM.StittG.SonL.EnioutinaE. Y. (2023). Probiotics and their bioproducts: a promising approach for targeting methicillin-resistant Staphylococcus aureus and vancomycin-resistant enterococcus. Microorganisms11:2393. 10.3390/microorganisms1110239337894051JengH.-S.YanT.-R. (2022). Lactiplantibacillus plantarum E51 protects against Clostridioides difficile-induced damages on Caco-2 intestinal barrier functions. Arch. Microbiol.204:290. 10.1007/s00203-022-02837-635503482JiJ.JinW.LiuS.JiaoZ.LiX. (2023). Probiotics, prebiotics, and postbiotics in health and disease. MedComm4:e420. 10.1002/mco2.42037929014JyungS.KangJ.-W.KangD.-H. (2023). Inactivation of Listeria monocytogenes through the synergistic interaction between plasma-activated water and organic acid. Food Res. Int.167:112687. 10.1016/j.foodres.2023.11268737087257KarabiyikliS.SengunI. (2017). “Beneficial effects of acetic acid bacteria and their food products,” in Acetic Acid Bacteria, eds KarabiyikliE.SengunI. (CRC Press), 321–342.KerrG. K.SnellingA. M. (2009). Pseudomonas aeruginosa: a formidable and ever-present adversary. J. Hosp. Infect.73, 338–344. 10.1016/j.jhin.2009.04.02019699552KleerebezemM.BoekhorstJ.KranenburgR. V.MolenaarD.KuipersO. P.LeerE.. (2003). Complete genome sequence of Lactobacillus plantarum WCFS1. Proc. Nat. Acad. Sci. U. S. A.100, 1990–1995. 10.1073/pnas.033770410012566566KrzywinskiM.ScheinJ.BirolI.ConnorsJ.GascoyneR.HorsmanD.. (2009). Circos: an information aesthetic for comparative genomics. Genome Res.19, 1639–1645. 10.1101/gr.092759.10919541911KuenemanG. J.WoodhamsD. C.HarrisR.ArcherH. M.KnightR.McKenzieV. J. (2016). Probiotic treatment restores protection against lethal fungal infection lost during amphibian captivity. Proc. R. Soc. B Biol. Sci.283:20161553. 10.1098/rspb.2016.155327655769KumarS.AhmadM. F.NathP.RoyR.BhattacharjeeR.ShamaE.. (2023). Controlling intestinal infections and digestive disorders using probiotics. J. Med. Food26, 705–720. 10.1089/jmf.2023.006237646629LiY. J.YuS.WengP. F.WuZ. F.LiuY. N. (2023). Purification and antimicrobial mechanism of a novel bacteriocin produced by Lactiplantibacillus plantarum FB-2. Lwt Food Sci. Technol.185:115123. 10.1016/j.lwt.2023.115123LowyF. D. (1998). Staphylococcus aureus infections. New Engl. J. Med.339, 520–532. 10.1056/NEJM199808203390806LyuM.BaiY.OriharaK.MiyanagaK.YamamotoN. (2023). GAPDH Released from Lactobacillus johnsonii MG enhances barrier function by upregulating genes associated with tight junctions. Microorganisms11:1393. 10.3390/microorganisms1106139337374895MangellP.NejdforsP.WangM.AhrnéS.WeströmB.ThorlaciusH.. (2002). Lactobacillus plantarum 299v inhibits Escherichia coli-induced intestinal permeability. Dig. Dis. Sci.47, 511–516. 10.1023/A:101794753153611911334MaoY.ZhangX.XuZ. (2020). Identification of antibacterial substances of Lactobacillus plantarum DY-6 for bacteriostatic action. Food Sci. Nutr.8, 2854–2863. 10.1002/fsn3.158532566203MarquesC.DinisL.-T.SantosM. J.MotaJ.VilelaA. (2023). Beyond the bottle: exploring health-promoting compounds in wine and wine-related products—extraction, detection, quantification, aroma properties, and terroir effects. Foods12:4277. 10.3390/foods1223427738231704MathipaM. G.ThantshaM. S. (2017). Probiotic engineering: towards development of robust probiotic strains with enhanced functional properties and for targeted control of enteric pathogens. Gut Pathog.9, 1–17. 10.1186/s13099-017-0178-928491143MendozaM. R.KimS. H.VasquezR.HwangI.-C.ParkY.-S.PaikH. D.. (2023). Bioinformatics and its role in the study of the evolution and probiotic potential of lactic acid bacteria. Food Sci. Biotechnol.32, 389–412. 10.1007/s10068-022-01142-836911331MengQ. F.LiuY. N.NieT.TangC.LyuF. X.BieX. M.. (2022). Plantaricin A, derived from Lactiplantibacillus plantarum, reduces the intrinsic resistance of gram-negative bacteria to hydrophobic antibiotics. Appl. Environ. Microbiol.88, e0037122. 10.1128/aem.00371-2235499329Min HsiuC.FengH. S.HuaC. J.FangL. M.ShuhC. C.ChenW. S. (2016). Antibacterial activity Lactobacillus plantarum isolated from fermented vegetables and investigation of the plantaricin genes. Afr. J. Microbiol. Res.10, 796–803. 10.5897/AJMR2016.7922MulatD. G.FeilbergA. (2015). GC/MS method for determining carbon isotope enrichment and concentration of underivatized short-chain fatty acids by direct aqueous solution injection of biogas digester samples. Talanta143, 56–63. 10.1016/j.talanta.2015.04.06526078128OftedalT. F. (2023). Bacteriocins: from discovery to characterization and applications.QiJ.TanC. K.HashimiS. M.ZulfikerA. H. M.GoodD.WeiM. Q. (2014). Toad glandular secretions and skin extractions as anti-inflammatory and anticancer agents. Evid. Based Comp. Altern. Med. 2014:312684. 10.1155/2014/31268424734105Rebaza-CardenasT. D.Silva-CajaleónK.SabaterC.DelgadoS.Montes-VillanuevaN. D.Ruas-MadiedoP. (2023). “Masato de yuca” and “Chicha de Siete Semillas” two traditional vegetable fermented beverages from Peru as source for the isolation of potential probiotic bacteria. Probiot. Antimicrob. Prot.15, 300–311. 10.1007/s12602-021-09836-x34453308SakaridisI.SoultosN.DovasC. I.PapavergouE.AmbrosiadisI.KoidisP. (2012). Lactic acid bacteria from chicken carcasses with inhibitory activity against Salmonella spp. and Listeria monocytogenes. Anaerobe18, 62–66. 10.1016/j.anaerobe.2011.09.00921983108SiezenJ. RTzenevaV. A.CastioniA.WelsM.PhanH. T. K.. (2010). Phenotypic and genomic diversity of Lactobacillus plantarum strains isolated from various environmental niches. Environ. Microbiol.12, 758–773. 10.1111/j.1462-2920.2009.02119.x20002138SinghK. S.KumarS.MohantyA. K.GroverS.KaushikJ. K. (2018). Mechanistic insights into the host-microbe interaction and pathogen exclusion mediated by the Mucus-binding protein of Lactobacillus plantarum. Sci. Rep.8:14198. 10.1038/s41598-018-32417-y30242281SonY. S.KwonM.SonN.KimS.-K.SonM.-Y. (2023). Multilayer coating with red ginseng dietary fiber improves intestinal adhesion and proliferation of probiotics in human intestinal epithelial models. J. Microbiol. Biotechnol.33:1309. 10.4014/jmb.2305.0501337528560TanJ.McKenzieC.PotamitisM.ThorburnA. N.MackayC. R.MaciaL.. (2014). The role of short-chain fatty acids in health and disease. Adv. Immunol.121, 91–119. 10.1016/B978-0-12-800100-4.00003-924388214TangH.HuangW.YaoY.-F. (2023). The metabolites of lactic acid bacteria: classification, biosynthesis and modulation of gut microbiota. Microbial Cell10:49. 10.15698/mic2023.03.79236908281TatusovL. R.FedorovaN. D.JacksonJ. D.JacobsA. R.KiryutinB.KooninE. V.. (2003). The COG database: an updated version includes eukaryotes. BMC Bioinformat.4:41. 10.1186/1471-2105-4-4112969510TodorovS. D.FrancoB. D. (2010). Lactobacillus plantarum: characterization of the species and application in food production. Food Rev. Int.26, 205–229. 10.1080/87559129.2010.48411324499153Valan ArasuM.JungM. W.KimD. H.ParkH. S.IlavenilS.Al-DhabiN. A.. (2015). Identification and phylogenetic characterization of novel Lactobacillus plantarum species and their metabolite profiles in grass silage. Ann. Microbiol.65, 15–25. 10.1007/s13213-014-0830-2van HeelA. J.de JongA.SongC. X.VielJ. H.KokJ.KuipersO. P. (2018). BAGEL4: a user-friendly web server to thoroughly mine RiPPs and bacteriocins. Nucleic Acids Res.46, W278–W81. 10.1093/nar/gky38329788290Van ZylW. F.DeaneS. M.DicksL. M. T. (2020). Molecular insights into probiotic mechanisms of action employed against intestinal pathogenic bacteria. Gut Microbes12:1831339. 10.1080/19490976.2020.183133933112695VidraA.NémethÁ. (2018). Bio-produced acetic acid: a review. Period. Polytech. Chem. Eng.62, 245–256. 10.3311/PPch.11004WangJ.LiL.ZhaoX.ZhouZ. J. (2015). Partial characteristics and antimicrobial mode of pediocin produced by PA003. Ann. Microbiol.65, 1753–1762. 10.1007/s13213-014-1014-9WangJ.-Y.ChenL.ZhengZ.WangQ.GuoJ.XuL.. (2012). Cinobufocini inhibits NF-κB and COX-2 activation induced by TNF-α in lung adenocarcinoma cells. Oncol. Rep.27, 1619–1624. 10.3892/or.2012.164722267101WangY.WeiY.ShangN.LiP. (2022). Synergistic inhibition of plantaricin E/F and lactic acid against Aeromonas hydrophila LPL-1 reveals the novel potential of class IIb bacteriocin. Front. Microbiol.13:774184. 10.3389/fmicb.2022.77418435242114YangJ.MengL.LiY.HuangH. (2023). Strategies for applying probiotics in the antibiotic management of Clostridioides difficile infection. Food Funct.14, 8711–8733. 10.1039/D3FO02110F37725066YangK. M.JiangZ. Y.ZhengC. T.WangL.YangX. F. (2014). Effect of Lactobacillus plantarum on diarrhea and intestinal barrier function of young piglets challenged with enterotoxigenic Escherichia coli K88. J. Anim. Sci.92, 1496–1503. 10.2527/jas.2013-661924492550YangS.LiuY.YangN.LanY.LanW.FengJ.. (2022). The gut microbiome and antibiotic resistome of chronic diarrhea rhesus macaques (Macaca mulatta) and its similarity to the human gut microbiome. Microbiome10, 1–13. 10.1186/s40168-021-01218-335139923YeK.LiP.GuQ. (2020). Complete genome sequence analysis of a strain Lactobacillus pentosus ZFM94 and its probiotic characteristics. Genomics112, 3142–3149. 10.1016/j.ygeno.2020.05.01532450257ZhangL.MaH.KulyarM. F.-A.PanH.LiK.LiA.. (2022). Complete genome analysis of Lactobacillus fermentum YLF016 and its probiotic characteristics. Microb. Pathog.162:105212. 10.1016/j.micpath.2021.10521234597776ZhangY.ZhangY.LiuF.MaoY.ZhangY.ZengH.. (2023). Mechanisms and applications of probiotics in prevention and treatment of swine diseases. Porcine Health Manag.9:5. 10.1186/s40813-022-00295-6ZhaoN.LiM.LuoJ.WangS.LiuS.WangS.. (2017). Impacts of canine distemper virus infection on the giant panda population from the perspective of gut microbiota. Sci. Rep.7:39954. 10.1038/srep3995428051146