Articles Service
Research article
Safety Assessment of Lactiplantibacillus (formerly Lactobacillus) plantarum Q180
1Probiotics Research Laboratory, Chong Kun Dang Bio Research Institute (CKDBIO), Gyeonggi 15064, Republic of Korea
2Department of Life Science, Chung-Ang University, Seoul 06974, Republic of Korea
J. Microbiol. Biotechnol. 2021; 31(10): 1420-1429
Published October 28, 2021 https://doi.org/10.4014/jmb.2106.06066
Copyright © The Korean Society for Microbiology and Biotechnology.
Abstract
Keywords
Graphical Abstract
Introduction
Probiotics are defined as “live microorganisms which when administered in adequate amounts confer a health benefit on the host” (FAO/WHO 2002). While probiotics have several known benefits, there have been safety issues related to their use in humans and animals. In 2002, the Food and Agriculture Organization/World Health Organization (FAO/WHO) reported that probiotics may cause side-effect such as systemic infection, deleterious metabolism, and excessive immune response in susceptible subjects or exhibit harmful gene transfer (FAO/WHO 2002). Furthermore, the European Food Safety Authority (EFSA) announced safety assessment guidance for probiotics by focusing on antimicrobial resistance (AMR) [1] and made it mandatory to examine susceptibility of all bacterial strains used as feed additives to the most relevant antibiotics. As a basic requirement, the minimum inhibitory concentration (MIC) should be determined for nine antibiotics (ampicillin, vancomycin, tetracycline, and others). The resistance of a bacterial strain to any specific antibiotic higher than the microbiological cut-off values defined by EFSA is deemed indicative of the presence of acquired resistance. Further, additional information is needed on the genetic basis of the AMR. For the genetic safety evaluation of microorganisms, EFSA recently recommended the taxonomic identification and characterization of their potential functional traits of concern, which may include virulence factors (VFs) and AMR [2, 3].
In our previous studies, strain Q180 isolated from the feces of a healthy Korean adult was found to exhibit blood triglyceride lowering effects in vivo and in clinical trials [21-23]. We intended to register strain Q180 as a probiotic strain with the Korea Food and Drug Administration (KFDA) and commercialize it as a probiotic product in Korea. Therefore, this study was aimed to verify the safety of strain Q180 through bioinformatic analysis and other safety tests. In this direction, we performed MIC test, hemolysis assay, biogenic amine (BA) production, and bioinformatic analyses to identify the genes related to AMR, VFs, or MGEs.
Materials and Methods
Bacterial Strains and Culture Conditions
Strain Q180 was provided by the Korea Food Research Institute (KFRI, Korea) [23], and
Whole Genome Sequencing of Strain Q180
The genomic DNA of strain Q180 was extracted and completely sequenced using the combination of PacBio RS II and Illumina HiSeq 2500 platforms at Macrogen (Korea), as previously described [24]. In brief, the genomic DNA of strain Q180 was sequenced using PacBio RS II with a 10 kb library, and the resulting sequencing reads were
Phylogenetic and Genome-Related Analyses
Taxonomic identification of strain Q180 was conducted through phylogenetic analyses based on the 16S rRNA gene and whole-genome sequences. For the 16S rRNA gene sequence-based phylogenetic analysis, the 16S rRNA gene sequences of strain Q180 and closely related type strains were aligned using the fast secondary-structure-aware infernal aligner available in the ribosomal database project [27]. A maximum-likelihood (ML) tree with bootstrap values (1,000 replications) was constructed using the MEGA7 software [28]. The 16S rRNA gene sequence similarities between strain Q180 and closely related type strains were calculated using the EzTaxon-e server (http://www.ezbiocloud.net/). For the genome-based phylogenetic analysis, 92 housekeeping core genes from the genomes of strain Q180 and closely related type strains were extracted using the UBCG pipeline (www.ezbiocloud.net/tools/ubcg) [29], and an ML tree with bootstrap values (1,000 replications) based on the concatenated nucleotide sequences of the 92 housekeeping core genes was constructed using the MEGA7 software.
The genome relatedness between strain Q180 and closely related type strains was evaluated through average nucleotide identity (ANI) and digital DNA-DNA hybridization (DDH) analyses using the Orthologous Average Nucleotide Identity Tool software (www.ezbiocloud.net/sw/oat) [30] and the server-based Genome-to-Genome Distance Calculator (http://ggdc.dsmz.de/distcalc2.php) [31], respectively.
Genomic Analysis
The circular map of the assembled genome of strain Q180 was visualized using a web-based CGview program [32]. The complete genome of strain Q180 was submitted to GenBank for the gene prediction and functional annotation. Protein-coding sequences in the genomes of strain Q180 and
Bioinformatic Analysis of AMR Genes, VF- and Toxin-Related Genes, and MGEs
AMR genes present in the genomes of all
VF- and toxin-related genes, including those associated with enterotoxin, leucotoxin, cytolysin, cytotoxin K, hemolysis, BA production, hyaluronidase, aggregation, enterococcal surface protein, endocarditis antigen, collagen adhesion, cereulide, sex pheromone, and serine protease were searched in the genomes of strain Q180,
To assess the transferability of AMR/VF genes in strain Q180, prophages and insertion sequences (including transposons), which are representative MGEs, were investigated using PHAge Search Tool Enhanced Release (PHASTER) [39] and ISfinder [40], respectively.
Prophage Induction Test in Strain Q180
The inducibility of the putative prophage sequences as active phages in strain Q180 was evaluated using a mitomycin C approach, as described previously [41]. In brief, strain Q180 was cultured to an optical density of 0.2 (600 nm) at 37°C in 100 ml MRS broth supplemented with 10 mM CaCl2 and mitomycin C was added to be a final concentration of 0.6 μg/ml, followed by additional incubation at 30°C for 18 h. Phage DNA was isolated from the cell culture supernatants, according to the phage DNA extraction protocol described previously [42], and detected on 1.0% (w/v) agarose gel.
Determination of MICs of Antibiotics
Phenotypic resistance of strain Q180 and
Phenotypic Tests for Hemolytic and Gelatinase Activities and BA-producing Ability
Phenotypic tests for hemolytic and gelatinase activities of strain Q180 were performed, and BA-producing ability of strain Q180 was investigated.
To evaluate hemolytic activity, each strain was streaked on blood agar containing 5% (w/v) sheep blood (MB Cell, Korea) and the plates were incubated at 37°C for 48 h. The hemolytic activity (
The BA-producing ability of each strain was assessed according to a previously described procedure [45]. In brief, the test strains were cultivated in MRS broth containing histidine, tyrosine, ornithine, and lysine (each 0.25%) at 37°C for 2 days. The culture broths were syringe-filtered (0.2 μm; Biofact, Korea). For 9-fluorenylmethoxy carbonyl (FMOC) derivatization, 20 μl of the filtrates or standards (5, 10, and 20 μM of histamine, tyramine, putrescine, and cadaverine; Sigma-Aldrich, USA) and 200 μl of 1.5 mM FMOC (in acetone) were added to 200 μl of 0.5 M sodium borate buffer (pH 8.5) containing 20 μM norvaline (internal standard) and vigorously mixed. After 3 min incubation at room temperature (dark), 50 μl of 10 mM glycine in 0.5 M sodium borate buffer was added to the mixture to remove excess FMOC. FMOC-derivatized BA was analyzed by high-performance liquid chromatography (HPLC; Shimadzu, Japan) equipped with a reverse-phase C18 column (250 × 4.6 mm) and a fluorescence detector (RF-10AXL) as described by Brückner
Results and Discussion
Genomic Sequencing and Taxonomic Identification of Strain Q180
The genome of strain Q180 was completely sequenced with high quality and taxonomically identified through 16S rRNA gene- and whole genome-based phylogenetic analyses. Phylogenetic analysis based on 16S rRNA gene sequences revealed that strain Q180 formed a phylogenic lineage with both
-
Fig. 1. Maximum-likelihood phylogenetic trees showing the relationships between strain Q180 and the closely related
Lactiplantibacillus type strains, based on the 16S rRNA gene (A) and concatenated 92 housekeeping core gene (B) sequences. Levilactobacillus tongjiangensis LMG 26013T (JQCL00000000) was used as the outgroup (not shown). The bars, 0.005 and 0.05, represent changes per nucleotide.
-
Fig. 2. Heat-maps showing pair-wise average nucleotide identity (ANI) and digital DNA-DNA hybridization (DDH) values of strain Q180 and the closely related
Lactiplantibacillus type strains.
Genomic Features of Strain Q180
The key genomic features of strain Q180, including GC skew, protein coding sequences (CDSs), COG categories, and G+C contents, are graphically depicted in Fig. 3. The general genomic features of strain Q180 were summarized and compared with those of the type strain of
-
Table 1 . General features of the genomes of strain Q180 and
Lp. plantarum DSM 20174T.Feature Q180 DSM 20174T No. of contigs 1 2 Chromosome size (bp) 3,197,263 3,242,936 Plasmid size (bp) – 7,218 GC content (%) 44.6 44.5 Total genes 3,301 3,060 Protein coding genes 3,049 2,921 rRNA genes 16 16 tRNA genes 68 71 Pseudogenes 168 48 Proteins with function prediction 2,373 2,518 Proteins assigned to COG 2,526 2,525 GenBank acc. no. CP073753 CP039121
-
Fig. 3. A circular map representing the chromosome of strain Q 180.
Forward strand and reverse strand coding sequences on the outermost two circles of the map are differently colored according to the COG categories of the right side. GC skews (GC skew+: green, GC skew-: pink) and G+C content (yellow) are drawn on the third and fourth circles, respectively.
Bioinformatic Analysis of AMR Genes and VF- and Toxin-Related Genes
During the evaluation of the safety of probiotic strains, their phenotypic properties related to AMR are considered important (WHO 2019). As the genomic analysis of microbes provides information on AMR genes as well as their potential transferability, genomic analysis is necessary for the evaluation of the potential risks of new probiotic strains. Previous genomic analysis studies have shown that many
-
Table 2 . The distribution of antimicrobial resistance (AMR) genes in the genomes of
Lactiplantibacillus plantarum strains.ARO category Target antibiotics Antibiotic gene database† CARD ARG-ANNOT ResFinder P S L P S L P S L Antibiotic inactivation enzyme Carbapenem, aminoglycoside, phenicol antibiotic, lincosamide, and streptogramin 14 1 4 11 2 1 2 1 0 Antibiotic target protection protein Tetracycline and lincosamide 0 1 4 1 0 7 1 0 0 Antibiotic efflux pump Tetracycline, aminoglycoside, fluoroquinolone, fosfomycin, and macrolide 0 0 66 0 0 16 0 0 0 ARO, antibiotic resistance ontology
†The AMR genes were searched in a total of 583
Lp. plantarum genomes retrieved from GenBank based on the CARD, ARGANNOT, and ResFinder databases. The results represent the number of genomes containing at least one or more AMR genes. P, perfect hits: 100% similarity sequences with the database sequences; S, strict hits: 90%–100% similarity sequences with the database sequences; L, loose hit: 50%–90% similarity sequences with the database sequences.
The VF- and toxin-related genes were bioinformatically analyzed in the genomes of strain Q180,
-
Table 3 . Bioinformatic analysis for the presence of putative virulence factor- and toxin-related genes in the genomes of strain Q180,
Lp. plantarum DSM 20174T,E. faecium ATCC 19434T, andS. aureus ATCC 6538.Class Gene Lp. plantarum E. faecium S. aureus Q180 DSM 20174T ATCC 19434T ATCC 6538 Enterotoxin selk, selq, set – – – + Leucotoxin lukD – – – + Cytolysin cylA – – – – Cytotoxin K cytK – – – – Hemolysin hbl – – – + Gelatinase gelE – – – + Amino acid decarboxylase hdc1 ,hdc2 – – – – tdc – – + – odc – – – ldc – – – – Hyaluronidase hyl – – + + Aggregation substance asa1 – – – Enterococcal surface protein esp – – – – Endocarditis antigen efaA – – – – Adhesion of collagen ace – – – – Cereulide cesA – – – – Sex pheromones ccf ,cob ,cpd – – – – Serine protease sprE – – – + Transposon-related genes int ,intTN – – – –
To assess the transferability of AMR/VF or toxin genes in strain Q180, the presence of MGEs, including prophages or insertion sequences, was bioinformatically investigated. Two putative intact prophages (regions 2 and 3) and three incomplete phage sequences were identified in the genome of strain Q180 (Fig. 4A), suggesting the possibility of active lateral gene transfer by phage infection. However, the prophage sequence of region 2 did not have an integrase protein, a key protein for prophage induction (Fig. 4B). Although some key proteins for prophage induction were identified from the prophage sequence of region 3, the sequences were more similar to bacterial protein-coding genes than viral genes. In-vitro induction test of the prophages also showed that the prophages were not induced (data not shown). These results suggest that the prophage sequences identified in the genome of strain Q180 may be not inducible and strain Q180 is safe from the transferability of AMR/VF or toxin genes by a prophage.
-
Fig. 4. Bioinformatic analysis of phage sequences in the genome of strain Q180.
A, the locations of putative prophage regions identified in the chromosome of strain Q180; B and C, the genomic maps of putative prophage regions 2 and 3 in panel A. The phage sequences were analyzed using PHASTER.
The bioinformatic analysis also showed that two types of insertion sequences (IS1182 and ISL3) annotated as transposases probably originating from
Antibiotic Susceptibility Assay
Although most strains, including
The antibiotic resistance assays using the E-test strip and broth microdilution susceptibility assay methods showed that the MIC values of the seven antibiotics suggested by the EFSA guideline for strain Q180 were similar to those of
-
Table 4 . Minimum inhibitory concentrations (MICs) of antibiotics for strain Q180 and
Lp. plantarum ATCC 14917T.Antibiotic MIC (mg/l) Strain Q180 Strain ATCC 14917T Cut-off value‡ E-test strip Broth assay† E-test strip Broth assay† Ampicillin 1.0±0.0 0.25±0.0 0.38±0.0 0.5±0.0 2.0 Gentamicin 1.5±0.0 8.0±0.0 1.0±0.0 8.0±0.0 16.0 Kanamycin 24.0±0.0 64.0±0.0 24.0±0.0 64.0±0.0 64.0 Erythromycin 0.25±0.0 0.4±0.1 0.25±0.0 1.0±0.0 1.0 Clindamycin 0.25±0.0 0.125±0.0 0.19±0.0 0.125±0.0 2.0 Tetracycline 32.0±0.0 32.0±0.0 24.0±0.0 32.0±0.0 32.0 Chloramphenicol 8.0±0.0 8.0±0.0 6.0±0.0 8.0±0.0 8.0 †The assay was performed using the standard broth microdilution susceptibility assay protocol recommended by the National Committee for Clinical Laboratory Standards [43].
‡Microbiological cut-off values for antibiotics for
Lp. plantarum , as provided by the EFSA 2012 guideline.
Hemolytic and Gelatinase Activities and BA-Producing Ability
To confirm the safety of a probiotic strain, it is necessary to phenotypically verify its virulence or toxin-related properties through in vitro tests. In this study, the hemolytic and gelatinase activities and BA-producing ability of strain Q180 were tested, together with
-
Fig. 5. Phenotypic tests. A, hemolysis test; B, gelatinase activity; and C, biogenic amine production of 1, strain Q180; 2,
Lp. plantarum ATCC 14917T; 3,E. faecium ATCC 19434T; and 4,S. aureus ATCC 6538.
Gelatinase is a hydrophobic metalloprotease that cleaves insulin, casein, hemoglobin, collagen, and gelatin. Because it has been reported that pathogenic bacteria such as enterococci responsible for endocarditis and bacteremia have the ability to penetrate the tissue through their gelatinase activities [58], gelatinase is considered as one of the VFs. In the present study, no gelatinase activity was detected for strain Q180 as well as
BAs such as histamine, tyramine, cadaverine, and putrescine, mainly formed by the decarboxylation of their corresponding amino acids, act as neurotransmitters in organisms; however, excessive intake of BAs can cause diseases or disorders such as diarrhea, vomiting, sweating, and tachycardia [61]. Here, we found that strain Q180 had no BA-producing ability, as observed for
In this study, the safety of strain Q180 known to exhibit postprandial lipid-lowering effect was evaluated through bioinformatic and phenotypic analyses to test its potential as a probiotic strain. First, strain Q180 was subjected to phylogenetic and genome-related analyses, which showed that it belongs to
Supplemental Materials
Conflicts of Interest
The authors have no financial conflicts of interest to declare.
References
- Additives EPo, Feed PoSuiA. 2012. Guidance on the assessment of bacterial susceptibility to antimicrobials of human and veterinary importance.
EFSA J. 10 : 2740. - Li J, Ren F, Gu H, Li X, Gan B. 2011. Safety evaluation in vitro of
Enterococcus durans from Tibetan traditional fermented yak milk.J. Microbiol. 49 : 721-728. - Lee DY, Seo Y-S, Rayamajhi N, Kang ML, Lee SI, Yoo HS. 2009. Isolation, characterization, and evaluation of wild isolates of
Lactobacillus reuteri from pig feces.J. Microbiol. 47 : 663-672. - De Vries MC, Vaughan EE, Kleerebezem M, de Vos WM. 2006.
Lactobacillus plantarum -survival, functional and potential probiotic properties in the human intestinal tract.Int. Dairy J. 16 : 1018-1028. - Piqué N, Berlanga M, Miñana-Galbis D. 2019. Health benefits of heat-killed (Tyndallized) probiotics: An overview.
Int. J. Mol. Sci. 20 : 2534. - Liu Y-W, Liong M-T, Tsai Y-C. 2018. New perspectives of
Lactobacillus plantarum as a probiotic: The gut-heart-brain axis.J. Microbiol. 56 : 601-613. - Zhang M, Jiang Y, Cai M, Yang Z. 2020. Characterization and ACE inhibitory activity of fermented milk with probiotic
Lactobacillus plantarum K25 as analyzed by GC-MS-based metabolomics approach.J. Microbiol. Biotechnol. 30 : 903-911. - Zheng J, Wittouck S, Salvetti E, Franz CM, Harris HM, Mattarelli P,
et al . 2020. A taxonomic note on the genusLactobacillus : Description of 23 novel genera, emended description of the genusLactobacillus Beijerinck 1901, and union ofLactobacillaceae andLeuconostocaceae .Int. J. Syst. Evol. Microbiol. 70 : 2782-2858. - Mukerji P, Roper JM, Stahl B, Smith AB, Burns F, Rae JC,
et al . 2016. Safety evaluation of AB-LIFE®(Lactobacillus plantarum CECT 7527, 7528 and 7529): Antibiotic resistance and 90-day repeated-dose study in rats.Food Chem. Toxicol. 92 : 117-128. - Olek A, Woynarowski M, Ahrén IL, Kierkuś J, Socha P, Larsson N,
et al . 2017. Efficacy and safety ofLactobacillus plantarum DSM 9843 (LP299V) in the prevention of antibiotic-associated gastrointestinal symptoms in children-randomized, double-blind, placebo-controlled study.J. Pediatr. 186 : 82-86. - Gotteland M, Cires MJ, Carvallo C, Vega N, Ramirez MA, Morales P,
et al . 2014. Probiotic screening and safety evaluation ofLactobacillus strains from plants, artisanal goat cheese, human stools, and breast milk.J. Med. Food 17 : 487-495. - Goel A, Halami PM, Tamang JP. 2020. Genome analysis of
Lactobacillus plantarum isolated from some Indian fermented foods for bacteriocin production and probiotic marker genes.Front. Microbiol. 11 : 40. - Park DM, Bae J-H, Kim MS, Kim H, Kang SD, Shim S,
et al . 2019. Suitability ofLactobacillus plantarum SPC-SNU 72-2 as a probiotic starter for sourdough fermentation.J. Microbiol. Biotechnol. 29 : 1729-1738. - Gueimonde M, Sánchez B, de Los Reyes-Gavilán CG, Margolles A. 2013. Antibiotic resistance in probiotic bacteria.
Front. Microbiol. 4 : 202. - Devirgiliis C, Zinno P, Perozzi G. 2013. Update on antibiotic resistance in foodborne
Lactobacillus andLactococcus species.Front. Microbiol. 4 : 301. - Abriouel H, Muñoz MdCC, Lerma LL, Montoro BP, Bockelmann W, Pichner R,
et al . 2015. New insights in antibiotic resistance ofLactobacillus species from fermented foods.Food Res. Int. 78 : 465-481. - Campedelli I, Mathur H, Salvetti E, Clarke S, Rea MC, Torriani S,
et al . 2019. Genus-wide assessment of antibiotic resistance inLactobacillus spp.Appl. Environ. Microbiol. 85 : e01738-18. - Chokesajjawatee N, Santiyanont P, Chantarasakha K, Kocharin K, Thammarongtham C, Lertampaiporn S,
et al . 2020. Safety assessment of a nham starter cultureLactobacillus plantarum BCC9546 via whole-genome analysis.Sci. Rep. 10 : 10241. - Klarin B, Larsson A, Molin G, Jeppsson B. 2019. Susceptibility to antibiotics in isolates of
Lactobacillus plantarum RAPD‐type Lp299v, harvested from antibiotic treated, critically ill patients after administration of probiotics.Microbiologyopen 8 : e00642. - Shao Y, Zhang W, Guo H, Pan L, Zhang H, Sun T. 2015. Comparative studies on antibiotic resistance in
Lactobacillus casei andLactobacillus plantarum .Food Control 50 : 250-258. - Park S-Y, Cho S-A, Kim S-H, Lim S-D. 2014. Physiological characteristics and anti-obesity effect of
Lactobacillus plantarum Q180 isolated from feces.Korean J. Food Sci. Anim. Resour. 34 : 647-655. - Kwon J, Kim B, Lee C, Joung H, Kim B-K, Choi IS,
et al . 2020. Comprehensive amelioration of high-fat diet-induced metabolic dysfunctions through activation of the PGC-1α pathway by probiotics treatment in mice.PLoS One 15 : e0228932. - Park YE, Kim MS, Shim KW, Kim Y-I, Chu J, Kim B-K,
et al . 2020. Effects ofLactobacillus plantarum Q180 on postprandial lipid levels and intestinal environment: A double-blind, randomized, placebo-controlled, parallel trial.Nutrients 12 : 255. - Baek JH, Kim KH, Moon JY, Yeo S-H, Jeon CO. 2020.
Acetobacter oryzoeni sp. nov., isolated from Korean rice wine vinegar.Int. J. Syst. Evol. Microbiol. 70 : 2026-2033. - Berlin K, Koren S, Chin C-S, Drake JP, Landolin JM, Phillippy AM. 2015. Assembling large genomes with single-molecule sequencing and locality-sensitive hashing.
Nat. Biotechnol. 33 : 623-630. - Walker BJ, Abeel T, Shea T, Priest M, Abouelliel A, Sakthikumar S,
et al . 2014. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement.PLoS One 9 : e112963. - Cole JR, Wang Q, Cardenas E, Fish J, Chai B, Farris RJ,
et al . 2009. The Ribosomal Database Project: improved alignments and new tools for rRNA analysis.Nucleic Acids Res. 37 : 141-145. - Kumar S, Stecher G, Tamura K. 2016. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets.
Mol. Biol. Evol. 33 : 1870-1874. - Na S-I, Kim YO, Yoon S-H, Ha S-m, Baek I, Chun J. 2018. UBCG: up-to-date bacterial core gene set and pipeline for phylogenomic tree reconstruction.
J. Microbiol. 56 : 280-285. - Lee I, Kim YO, Park S-C, Chun J. 2016. OrthoANI: an improved algorithm and software for calculating average nucleotide identity.
Int. J. Syst. Evol. Microbiol. 66 : 1100-1103. - Meier-Kolthoff JP, Auch AF, Klenk H-P, Göker M. 2013. Genome sequence-based species delimitation with confidence intervals and improved distance functions.
BMC Bioinformatics 14 : 60. - Grant JR, Stothard P. 2008. The CGView Server: a comparative genomics tool for circular genomes.
Nucleic Acids Res. 36 : 181-184. - Huerta-Cepas J, Szklarczyk D, Heller D, Hernández-Plaza A, Forslund SK, Cook H,
et al . 2019. eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses.Nucleic Acids Res. 47 : 309-314. - Buchfink B, Xie C, Huson DH. 2015. Fast and sensitive protein alignment using DIAMOND.
Nat. Methods 12 : 59-60. - Alcock BP, Raphenya AR, Lau TT, Tsang KK, Bouchard M, Edalatmand A,
et al . 2020. CARD 2020: antibiotic resistome surveillance with the comprehensive antibiotic resistance database.Nucleic Acids Res. 48 : 517-525. - Gupta SK, Padmanabhan BR, Diene SM, Lopez-Rojas R, Kempf M, Landraud L,
et al . 2014. ARG-ANNOT, a new bioinformatic tool to discover antibiotic resistance genes in bacterial genomes.Antimicrob. Agents Chemother. 58 : 212-220. - Bortolaia V, Kaas RS, Ruppe E, Roberts MC, Schwarz S, Cattoir V,
et al . 2020. ResFinder 4.0 for predictions of phenotypes from genotypes.J. Antimicrob. Chemother. 75 : 3491-3500. - Chen L, Yang J, Yu J, Yao Z, Sun L, Shen Y,
et al . 2005. VFDB: a reference database for bacterial virulence factors.Nucleic Acids Res. 33 : 325-328. - Arndt D, Grant JR, Marcu A, Sajed T, Pon A, Liang Y,
et al . 2016. PHASTER: a better, faster version of the PHAST phage search tool.Nucleic Acids Res. 44 : 16-21. - Siguier P, Pérochon J, Lestrade L, Mahillon J, Chandler M. 2006. ISfinder: the reference centre for bacterial insertion sequences.
Nucleic Acids Res. 34 : 32-36. - Mercanti DJ, Carminati D, Reinheimer JA, Quiberoni A. 2011. Widely distributed lysogeny in probiotic lactobacilli represents a potentially high risk for the fermentative dairy industry.
Int. J. Food Microbiol. 144 : 503-510. - Binetti AG, Del Río B, Martìn MC, Alvarez MA. 2005. Detection and characterization of
Streptococcus thermophilus bacteriophages by use of the antireceptor gene sequence.Appl. Environ. Microbiol. 71 : 6096-6103. - Clinical and laboratory standards institute. 2012. Methods for dilution antimicrobial susceptibility tests for bacteria that grow aerobically; approved standard CLSI.
In: document M07-A9 . Wayne, PA, USA. - Beecher DJ, Schoeni JL, Wong A. 1995. Enterotoxic activity of hemolysin BL from
Bacillus cereus .Infect. Immun. 63 : 4423-4428. - Kim KH, Lee SH, Chun BH, Jeong SE, Jeon CO. 2019.
Tetragenococcus halophilus MJ4 as a starter culture for repressing biogenic amine (cadaverine) formation during saeu-jeot (salted shrimp) fermentation.Food Microbiol. 82 : 465-473. - Brückner H, Flassig S, Kirschbaum J. 2012. Determination of biogenic amines in infusions of tea (
Camellia sinensis ) by HPLC after derivatization with 9-fluorenylmethoxycarbonyl chloride (Fmoc-Cl).Amino Acids 42 : 877-885. - Chun J, Oren A, Ventosa A, Christensen H, Arahal DR, da Costa MS,
et al . 2018. Proposed minimal standards for the use of genome data for the taxonomy of prokaryotes.Int. J. Syst. Evol. Microbiol. 68 : 461-466. - Domig KJ, Zycka-Krzesinska J, Bardowski J, Morelli L. 2008. Molecular assessment of erythromycin and tetracycline resistance genes in lactic acid bacteria and bifidobacteria and their relation to the phenotypic resistance.
Int. J. Probiotics Prebiotics 3 : 271-280. - Hummel AS, Hertel C, Holzapfel WH, Franz CM. 2007. Antibiotic resistances of starter and probiotic strains of lactic acid bacteria.
Appl. Environ. Microbiol. 73 : 730-739. - Mayrhofer S, Van Hoek AH, Mair C, Huys G, Aarts HJ, Kneifel W,
et al . 2010. Antibiotic susceptibility of members of theLactobacillus acidophilus group using broth microdilution and molecular identification of their resistance determinants.Int. J. Food Microbiol. 144 : 81-87. - Lahtinen SJ, Boyle RJ, Margolles A, Frias R, Gueimonde M. 2009. Safety assessment of probiotics.
In: Prebiotics probiotics Sci. Technol . Springer-Verlag press, Berlin. - Cebeci A, Gürakan C. 2003. Properties of potential probiotic
Lactobacillus plantarum strains.Food Microbiol. 20 : 511-518. - Huycke MM, Spiegel CA, Gilmore MS. 1991. Bacteremia caused by hemolytic, high-level gentamicin-resistant
Enterococcus faecalis .Antimicrob. Agents Chemother. 35 : 1626-1634. - Sharma P, Tomar SK, Goswami P, Sangwan V, Singh R. 2014. Antibiotic resistance among commercially available probiotics.
Food Res. Int. 57 : 176-195. - Flórez AB, Egervärn M, Danielsen M, Tosi L, Morelli L, Lindgren S,
et al . 2006. Susceptibility ofLactobacillus plantarum strains to six antibiotics and definition of new susceptibility-resistance cutoff values.Microb. Drug Resist. 12 : 252-256. - Zhang F, Gao J, Wang B, Huo D, Wang Z, Zhang J,
et al . 2018. Whole-genome sequencing reveals the mechanisms for evolution of streptomycin resistance inLactobacillus plantarum .J. Dairy Sci. 101 : 2867-2874. - Facklam RR, Carvalho MdGS, Teixeira LM. 2002. History, taxonomy, biochemical characteristics, and antibiotic susceptibility testing of enterococci.
In: The enterococci: pathogenesis, molecular biology, and antibiotic resistance . ASM press, Washington, D.C., USA. - Zoletti GO, Pereira EM, Schuenck RP, Teixeira LM, Siqueira JF Jr, dos Santos KRN. 2011. Characterization of virulence factors and clonal diversity of
Enterococcus faecalis isolates from treated dental root canals.Res. Microbiol. 162 : 151-158. - Abouloifa H, Rokni Y, Bellaouchi R, Ghabbour N, Karboune S, Brasca M,
et al . 2019. Characterization of probiotic properties of antifungalLactobacillus strains isolated from traditional fermenting green olives.Probiotics Antimicrob. Proteins 12 : 683-696. - Kenfack CHM, Ngoufack FZ, Kaktcham PM, Wang YR, Zhu T, Yin L. 2018. Safety and antioxidant properties of five probiotic
Lactobacillus plantarum strains isolated from the digestive tract of honey bees.Am. J. Microbiol. Res. 6 : 1-8. - Wójcik W, Łukasiewicz M, Puppel K. 2021. Biogenic amines: formation, action and toxicity-a review.
J. Sci. Food Agric. 101 : 2634-2640. - Kim S, Huang E, Park S, Holzapfel W, Lim S-D. 2018. Physiological characteristics and anti-obesity effect of
Lactobacillus plantarum K10.Korean J. Food Sci. Anim. Resour. 38 : 554-569. - Gao Y, Liu Y, Sun M, Zhang H, Mu G, Tuo Y. 2020. Physiological function analysis of
Lactobacillus plantarum Y44 based on genotypic and phenotypic characteristics.J. Dairy Sci. 103 : 5916-5930. - Priyanka V, Ramesha A, Gayathri D, Vasudha M. 2020. Molecular characterization of non-biogenic amines producing
Lactobacillus plantarum GP11 isolated from traditional pickles using HRESI-MS analysis.J. Food Sci. Technol. 58 : 2216-2226.
Related articles in JMB
Article
Research article
J. Microbiol. Biotechnol. 2021; 31(10): 1420-1429
Published online October 28, 2021 https://doi.org/10.4014/jmb.2106.06066
Copyright © The Korean Society for Microbiology and Biotechnology.
Safety Assessment of Lactiplantibacillus (formerly Lactobacillus) plantarum Q180
Yoo Jin Kwon1†, Byung Hee Chun2†, Hye Su Jung2, Jaeryang Chu1, Hyunchae Joung1, Sung Yurb Park1, Byoung Kook Kim1, and Che Ok Jeon2*
1Probiotics Research Laboratory, Chong Kun Dang Bio Research Institute (CKDBIO), Gyeonggi 15064, Republic of Korea
2Department of Life Science, Chung-Ang University, Seoul 06974, Republic of Korea
Correspondence to:Che Ok Jeon, cojeon@cau.ac.kr
Abstract
The safety of the probiotic strain Q180, which exerts postprandial lipid-lowering effects, was bioinformatically and phenotypically evaluated. The genome of strain Q180 was completely sequenced, and single circular chromosome of 3,197,263 bp without any plasmid was generated. Phylogenetic and related analyses using16S rRNA gene and whole-genome sequences revealed that strain Q180 is a member of Lactiplantibacillus (Lp., formerly Lactobacillus) plantarum. Antimicrobial resistance (AMR) genes were bioinformatically analyzed using all Lp. plantarum genomes available in GenBank, which showed that AMR genes are present differently depending on Lp. plantarum strains. Bioinformatic analysis demonstrated that some mobile genetic elements such as prophages and insertion sequences were identified in the genome of strain Q180, but because they did not contain harmful genes such as AMR genes and virulence factor (VF)- and toxin-related genes, it was suggested that there is no transferability of harmful genes. The minimum inhibition concentrations of seven tested antibiotics suggested by the European Food Safety Authority guidelines were slightly lower than or equal to the microbiological cut-off values for Lp. plantarum. Strain Q180 did not show hemolytic and gelatinase activities and biogenic amine-producing ability. Taken together, this study demonstrated the safety of strain Q180 in terms of absence of AMR genes and VF- and toxin-related genes as a probiotic strain.
Keywords: Lactiplantibacillus plantarum Q180, safety, antibiotic resistance, virulence factor, probiotics
Introduction
Probiotics are defined as “live microorganisms which when administered in adequate amounts confer a health benefit on the host” (FAO/WHO 2002). While probiotics have several known benefits, there have been safety issues related to their use in humans and animals. In 2002, the Food and Agriculture Organization/World Health Organization (FAO/WHO) reported that probiotics may cause side-effect such as systemic infection, deleterious metabolism, and excessive immune response in susceptible subjects or exhibit harmful gene transfer (FAO/WHO 2002). Furthermore, the European Food Safety Authority (EFSA) announced safety assessment guidance for probiotics by focusing on antimicrobial resistance (AMR) [1] and made it mandatory to examine susceptibility of all bacterial strains used as feed additives to the most relevant antibiotics. As a basic requirement, the minimum inhibitory concentration (MIC) should be determined for nine antibiotics (ampicillin, vancomycin, tetracycline, and others). The resistance of a bacterial strain to any specific antibiotic higher than the microbiological cut-off values defined by EFSA is deemed indicative of the presence of acquired resistance. Further, additional information is needed on the genetic basis of the AMR. For the genetic safety evaluation of microorganisms, EFSA recently recommended the taxonomic identification and characterization of their potential functional traits of concern, which may include virulence factors (VFs) and AMR [2, 3].
In our previous studies, strain Q180 isolated from the feces of a healthy Korean adult was found to exhibit blood triglyceride lowering effects in vivo and in clinical trials [21-23]. We intended to register strain Q180 as a probiotic strain with the Korea Food and Drug Administration (KFDA) and commercialize it as a probiotic product in Korea. Therefore, this study was aimed to verify the safety of strain Q180 through bioinformatic analysis and other safety tests. In this direction, we performed MIC test, hemolysis assay, biogenic amine (BA) production, and bioinformatic analyses to identify the genes related to AMR, VFs, or MGEs.
Materials and Methods
Bacterial Strains and Culture Conditions
Strain Q180 was provided by the Korea Food Research Institute (KFRI, Korea) [23], and
Whole Genome Sequencing of Strain Q180
The genomic DNA of strain Q180 was extracted and completely sequenced using the combination of PacBio RS II and Illumina HiSeq 2500 platforms at Macrogen (Korea), as previously described [24]. In brief, the genomic DNA of strain Q180 was sequenced using PacBio RS II with a 10 kb library, and the resulting sequencing reads were
Phylogenetic and Genome-Related Analyses
Taxonomic identification of strain Q180 was conducted through phylogenetic analyses based on the 16S rRNA gene and whole-genome sequences. For the 16S rRNA gene sequence-based phylogenetic analysis, the 16S rRNA gene sequences of strain Q180 and closely related type strains were aligned using the fast secondary-structure-aware infernal aligner available in the ribosomal database project [27]. A maximum-likelihood (ML) tree with bootstrap values (1,000 replications) was constructed using the MEGA7 software [28]. The 16S rRNA gene sequence similarities between strain Q180 and closely related type strains were calculated using the EzTaxon-e server (http://www.ezbiocloud.net/). For the genome-based phylogenetic analysis, 92 housekeeping core genes from the genomes of strain Q180 and closely related type strains were extracted using the UBCG pipeline (www.ezbiocloud.net/tools/ubcg) [29], and an ML tree with bootstrap values (1,000 replications) based on the concatenated nucleotide sequences of the 92 housekeeping core genes was constructed using the MEGA7 software.
The genome relatedness between strain Q180 and closely related type strains was evaluated through average nucleotide identity (ANI) and digital DNA-DNA hybridization (DDH) analyses using the Orthologous Average Nucleotide Identity Tool software (www.ezbiocloud.net/sw/oat) [30] and the server-based Genome-to-Genome Distance Calculator (http://ggdc.dsmz.de/distcalc2.php) [31], respectively.
Genomic Analysis
The circular map of the assembled genome of strain Q180 was visualized using a web-based CGview program [32]. The complete genome of strain Q180 was submitted to GenBank for the gene prediction and functional annotation. Protein-coding sequences in the genomes of strain Q180 and
Bioinformatic Analysis of AMR Genes, VF- and Toxin-Related Genes, and MGEs
AMR genes present in the genomes of all
VF- and toxin-related genes, including those associated with enterotoxin, leucotoxin, cytolysin, cytotoxin K, hemolysis, BA production, hyaluronidase, aggregation, enterococcal surface protein, endocarditis antigen, collagen adhesion, cereulide, sex pheromone, and serine protease were searched in the genomes of strain Q180,
To assess the transferability of AMR/VF genes in strain Q180, prophages and insertion sequences (including transposons), which are representative MGEs, were investigated using PHAge Search Tool Enhanced Release (PHASTER) [39] and ISfinder [40], respectively.
Prophage Induction Test in Strain Q180
The inducibility of the putative prophage sequences as active phages in strain Q180 was evaluated using a mitomycin C approach, as described previously [41]. In brief, strain Q180 was cultured to an optical density of 0.2 (600 nm) at 37°C in 100 ml MRS broth supplemented with 10 mM CaCl2 and mitomycin C was added to be a final concentration of 0.6 μg/ml, followed by additional incubation at 30°C for 18 h. Phage DNA was isolated from the cell culture supernatants, according to the phage DNA extraction protocol described previously [42], and detected on 1.0% (w/v) agarose gel.
Determination of MICs of Antibiotics
Phenotypic resistance of strain Q180 and
Phenotypic Tests for Hemolytic and Gelatinase Activities and BA-producing Ability
Phenotypic tests for hemolytic and gelatinase activities of strain Q180 were performed, and BA-producing ability of strain Q180 was investigated.
To evaluate hemolytic activity, each strain was streaked on blood agar containing 5% (w/v) sheep blood (MB Cell, Korea) and the plates were incubated at 37°C for 48 h. The hemolytic activity (
The BA-producing ability of each strain was assessed according to a previously described procedure [45]. In brief, the test strains were cultivated in MRS broth containing histidine, tyrosine, ornithine, and lysine (each 0.25%) at 37°C for 2 days. The culture broths were syringe-filtered (0.2 μm; Biofact, Korea). For 9-fluorenylmethoxy carbonyl (FMOC) derivatization, 20 μl of the filtrates or standards (5, 10, and 20 μM of histamine, tyramine, putrescine, and cadaverine; Sigma-Aldrich, USA) and 200 μl of 1.5 mM FMOC (in acetone) were added to 200 μl of 0.5 M sodium borate buffer (pH 8.5) containing 20 μM norvaline (internal standard) and vigorously mixed. After 3 min incubation at room temperature (dark), 50 μl of 10 mM glycine in 0.5 M sodium borate buffer was added to the mixture to remove excess FMOC. FMOC-derivatized BA was analyzed by high-performance liquid chromatography (HPLC; Shimadzu, Japan) equipped with a reverse-phase C18 column (250 × 4.6 mm) and a fluorescence detector (RF-10AXL) as described by Brückner
Results and Discussion
Genomic Sequencing and Taxonomic Identification of Strain Q180
The genome of strain Q180 was completely sequenced with high quality and taxonomically identified through 16S rRNA gene- and whole genome-based phylogenetic analyses. Phylogenetic analysis based on 16S rRNA gene sequences revealed that strain Q180 formed a phylogenic lineage with both
-
Figure 1. Maximum-likelihood phylogenetic trees showing the relationships between strain Q180 and the closely related
Lactiplantibacillus type strains, based on the 16S rRNA gene (A) and concatenated 92 housekeeping core gene (B) sequences. Levilactobacillus tongjiangensis LMG 26013T (JQCL00000000) was used as the outgroup (not shown). The bars, 0.005 and 0.05, represent changes per nucleotide.
-
Figure 2. Heat-maps showing pair-wise average nucleotide identity (ANI) and digital DNA-DNA hybridization (DDH) values of strain Q180 and the closely related
Lactiplantibacillus type strains.
Genomic Features of Strain Q180
The key genomic features of strain Q180, including GC skew, protein coding sequences (CDSs), COG categories, and G+C contents, are graphically depicted in Fig. 3. The general genomic features of strain Q180 were summarized and compared with those of the type strain of
-
Table 1 . General features of the genomes of strain Q180 and
Lp. plantarum DSM 20174T..Feature Q180 DSM 20174T No. of contigs 1 2 Chromosome size (bp) 3,197,263 3,242,936 Plasmid size (bp) – 7,218 GC content (%) 44.6 44.5 Total genes 3,301 3,060 Protein coding genes 3,049 2,921 rRNA genes 16 16 tRNA genes 68 71 Pseudogenes 168 48 Proteins with function prediction 2,373 2,518 Proteins assigned to COG 2,526 2,525 GenBank acc. no. CP073753 CP039121
-
Figure 3. A circular map representing the chromosome of strain Q 180.
Forward strand and reverse strand coding sequences on the outermost two circles of the map are differently colored according to the COG categories of the right side. GC skews (GC skew+: green, GC skew-: pink) and G+C content (yellow) are drawn on the third and fourth circles, respectively.
Bioinformatic Analysis of AMR Genes and VF- and Toxin-Related Genes
During the evaluation of the safety of probiotic strains, their phenotypic properties related to AMR are considered important (WHO 2019). As the genomic analysis of microbes provides information on AMR genes as well as their potential transferability, genomic analysis is necessary for the evaluation of the potential risks of new probiotic strains. Previous genomic analysis studies have shown that many
-
Table 2 . The distribution of antimicrobial resistance (AMR) genes in the genomes of
Lactiplantibacillus plantarum strains..ARO category Target antibiotics Antibiotic gene database† CARD ARG-ANNOT ResFinder P S L P S L P S L Antibiotic inactivation enzyme Carbapenem, aminoglycoside, phenicol antibiotic, lincosamide, and streptogramin 14 1 4 11 2 1 2 1 0 Antibiotic target protection protein Tetracycline and lincosamide 0 1 4 1 0 7 1 0 0 Antibiotic efflux pump Tetracycline, aminoglycoside, fluoroquinolone, fosfomycin, and macrolide 0 0 66 0 0 16 0 0 0 ARO, antibiotic resistance ontology.
†The AMR genes were searched in a total of 583
Lp. plantarum genomes retrieved from GenBank based on the CARD, ARGANNOT, and ResFinder databases. The results represent the number of genomes containing at least one or more AMR genes. P, perfect hits: 100% similarity sequences with the database sequences; S, strict hits: 90%–100% similarity sequences with the database sequences; L, loose hit: 50%–90% similarity sequences with the database sequences..
The VF- and toxin-related genes were bioinformatically analyzed in the genomes of strain Q180,
-
Table 3 . Bioinformatic analysis for the presence of putative virulence factor- and toxin-related genes in the genomes of strain Q180,
Lp. plantarum DSM 20174T,E. faecium ATCC 19434T, andS. aureus ATCC 6538..Class Gene Lp. plantarum E. faecium S. aureus Q180 DSM 20174T ATCC 19434T ATCC 6538 Enterotoxin selk, selq, set – – – + Leucotoxin lukD – – – + Cytolysin cylA – – – – Cytotoxin K cytK – – – – Hemolysin hbl – – – + Gelatinase gelE – – – + Amino acid decarboxylase hdc1 ,hdc2 – – – – tdc – – + – odc – – – ldc – – – – Hyaluronidase hyl – – + + Aggregation substance asa1 – – – Enterococcal surface protein esp – – – – Endocarditis antigen efaA – – – – Adhesion of collagen ace – – – – Cereulide cesA – – – – Sex pheromones ccf ,cob ,cpd – – – – Serine protease sprE – – – + Transposon-related genes int ,intTN – – – –
To assess the transferability of AMR/VF or toxin genes in strain Q180, the presence of MGEs, including prophages or insertion sequences, was bioinformatically investigated. Two putative intact prophages (regions 2 and 3) and three incomplete phage sequences were identified in the genome of strain Q180 (Fig. 4A), suggesting the possibility of active lateral gene transfer by phage infection. However, the prophage sequence of region 2 did not have an integrase protein, a key protein for prophage induction (Fig. 4B). Although some key proteins for prophage induction were identified from the prophage sequence of region 3, the sequences were more similar to bacterial protein-coding genes than viral genes. In-vitro induction test of the prophages also showed that the prophages were not induced (data not shown). These results suggest that the prophage sequences identified in the genome of strain Q180 may be not inducible and strain Q180 is safe from the transferability of AMR/VF or toxin genes by a prophage.
-
Figure 4. Bioinformatic analysis of phage sequences in the genome of strain Q180.
A, the locations of putative prophage regions identified in the chromosome of strain Q180; B and C, the genomic maps of putative prophage regions 2 and 3 in panel A. The phage sequences were analyzed using PHASTER.
The bioinformatic analysis also showed that two types of insertion sequences (IS1182 and ISL3) annotated as transposases probably originating from
Antibiotic Susceptibility Assay
Although most strains, including
The antibiotic resistance assays using the E-test strip and broth microdilution susceptibility assay methods showed that the MIC values of the seven antibiotics suggested by the EFSA guideline for strain Q180 were similar to those of
-
Table 4 . Minimum inhibitory concentrations (MICs) of antibiotics for strain Q180 and
Lp. plantarum ATCC 14917T..Antibiotic MIC (mg/l) Strain Q180 Strain ATCC 14917T Cut-off value‡ E-test strip Broth assay† E-test strip Broth assay† Ampicillin 1.0±0.0 0.25±0.0 0.38±0.0 0.5±0.0 2.0 Gentamicin 1.5±0.0 8.0±0.0 1.0±0.0 8.0±0.0 16.0 Kanamycin 24.0±0.0 64.0±0.0 24.0±0.0 64.0±0.0 64.0 Erythromycin 0.25±0.0 0.4±0.1 0.25±0.0 1.0±0.0 1.0 Clindamycin 0.25±0.0 0.125±0.0 0.19±0.0 0.125±0.0 2.0 Tetracycline 32.0±0.0 32.0±0.0 24.0±0.0 32.0±0.0 32.0 Chloramphenicol 8.0±0.0 8.0±0.0 6.0±0.0 8.0±0.0 8.0 †The assay was performed using the standard broth microdilution susceptibility assay protocol recommended by the National Committee for Clinical Laboratory Standards [43]..
‡Microbiological cut-off values for antibiotics for
Lp. plantarum , as provided by the EFSA 2012 guideline..
Hemolytic and Gelatinase Activities and BA-Producing Ability
To confirm the safety of a probiotic strain, it is necessary to phenotypically verify its virulence or toxin-related properties through in vitro tests. In this study, the hemolytic and gelatinase activities and BA-producing ability of strain Q180 were tested, together with
-
Figure 5. Phenotypic tests. A, hemolysis test; B, gelatinase activity; and C, biogenic amine production of 1, strain Q180; 2,
Lp. plantarum ATCC 14917T; 3,E. faecium ATCC 19434T; and 4,S. aureus ATCC 6538.
Gelatinase is a hydrophobic metalloprotease that cleaves insulin, casein, hemoglobin, collagen, and gelatin. Because it has been reported that pathogenic bacteria such as enterococci responsible for endocarditis and bacteremia have the ability to penetrate the tissue through their gelatinase activities [58], gelatinase is considered as one of the VFs. In the present study, no gelatinase activity was detected for strain Q180 as well as
BAs such as histamine, tyramine, cadaverine, and putrescine, mainly formed by the decarboxylation of their corresponding amino acids, act as neurotransmitters in organisms; however, excessive intake of BAs can cause diseases or disorders such as diarrhea, vomiting, sweating, and tachycardia [61]. Here, we found that strain Q180 had no BA-producing ability, as observed for
In this study, the safety of strain Q180 known to exhibit postprandial lipid-lowering effect was evaluated through bioinformatic and phenotypic analyses to test its potential as a probiotic strain. First, strain Q180 was subjected to phylogenetic and genome-related analyses, which showed that it belongs to
Supplemental Materials
Conflicts of Interest
The authors have no financial conflicts of interest to declare.
Fig 1.
Fig 2.
Fig 3.
Fig 4.
Fig 5.
-
Table 1 . General features of the genomes of strain Q180 and
Lp. plantarum DSM 20174T..Feature Q180 DSM 20174T No. of contigs 1 2 Chromosome size (bp) 3,197,263 3,242,936 Plasmid size (bp) – 7,218 GC content (%) 44.6 44.5 Total genes 3,301 3,060 Protein coding genes 3,049 2,921 rRNA genes 16 16 tRNA genes 68 71 Pseudogenes 168 48 Proteins with function prediction 2,373 2,518 Proteins assigned to COG 2,526 2,525 GenBank acc. no. CP073753 CP039121
-
Table 2 . The distribution of antimicrobial resistance (AMR) genes in the genomes of
Lactiplantibacillus plantarum strains..ARO category Target antibiotics Antibiotic gene database† CARD ARG-ANNOT ResFinder P S L P S L P S L Antibiotic inactivation enzyme Carbapenem, aminoglycoside, phenicol antibiotic, lincosamide, and streptogramin 14 1 4 11 2 1 2 1 0 Antibiotic target protection protein Tetracycline and lincosamide 0 1 4 1 0 7 1 0 0 Antibiotic efflux pump Tetracycline, aminoglycoside, fluoroquinolone, fosfomycin, and macrolide 0 0 66 0 0 16 0 0 0 ARO, antibiotic resistance ontology.
†The AMR genes were searched in a total of 583
Lp. plantarum genomes retrieved from GenBank based on the CARD, ARGANNOT, and ResFinder databases. The results represent the number of genomes containing at least one or more AMR genes. P, perfect hits: 100% similarity sequences with the database sequences; S, strict hits: 90%–100% similarity sequences with the database sequences; L, loose hit: 50%–90% similarity sequences with the database sequences..
-
Table 3 . Bioinformatic analysis for the presence of putative virulence factor- and toxin-related genes in the genomes of strain Q180,
Lp. plantarum DSM 20174T,E. faecium ATCC 19434T, andS. aureus ATCC 6538..Class Gene Lp. plantarum E. faecium S. aureus Q180 DSM 20174T ATCC 19434T ATCC 6538 Enterotoxin selk, selq, set – – – + Leucotoxin lukD – – – + Cytolysin cylA – – – – Cytotoxin K cytK – – – – Hemolysin hbl – – – + Gelatinase gelE – – – + Amino acid decarboxylase hdc1 ,hdc2 – – – – tdc – – + – odc – – – ldc – – – – Hyaluronidase hyl – – + + Aggregation substance asa1 – – – Enterococcal surface protein esp – – – – Endocarditis antigen efaA – – – – Adhesion of collagen ace – – – – Cereulide cesA – – – – Sex pheromones ccf ,cob ,cpd – – – – Serine protease sprE – – – + Transposon-related genes int ,intTN – – – –
-
Table 4 . Minimum inhibitory concentrations (MICs) of antibiotics for strain Q180 and
Lp. plantarum ATCC 14917T..Antibiotic MIC (mg/l) Strain Q180 Strain ATCC 14917T Cut-off value‡ E-test strip Broth assay† E-test strip Broth assay† Ampicillin 1.0±0.0 0.25±0.0 0.38±0.0 0.5±0.0 2.0 Gentamicin 1.5±0.0 8.0±0.0 1.0±0.0 8.0±0.0 16.0 Kanamycin 24.0±0.0 64.0±0.0 24.0±0.0 64.0±0.0 64.0 Erythromycin 0.25±0.0 0.4±0.1 0.25±0.0 1.0±0.0 1.0 Clindamycin 0.25±0.0 0.125±0.0 0.19±0.0 0.125±0.0 2.0 Tetracycline 32.0±0.0 32.0±0.0 24.0±0.0 32.0±0.0 32.0 Chloramphenicol 8.0±0.0 8.0±0.0 6.0±0.0 8.0±0.0 8.0 †The assay was performed using the standard broth microdilution susceptibility assay protocol recommended by the National Committee for Clinical Laboratory Standards [43]..
‡Microbiological cut-off values for antibiotics for
Lp. plantarum , as provided by the EFSA 2012 guideline..
References
- Additives EPo, Feed PoSuiA. 2012. Guidance on the assessment of bacterial susceptibility to antimicrobials of human and veterinary importance.
EFSA J. 10 : 2740. - Li J, Ren F, Gu H, Li X, Gan B. 2011. Safety evaluation in vitro of
Enterococcus durans from Tibetan traditional fermented yak milk.J. Microbiol. 49 : 721-728. - Lee DY, Seo Y-S, Rayamajhi N, Kang ML, Lee SI, Yoo HS. 2009. Isolation, characterization, and evaluation of wild isolates of
Lactobacillus reuteri from pig feces.J. Microbiol. 47 : 663-672. - De Vries MC, Vaughan EE, Kleerebezem M, de Vos WM. 2006.
Lactobacillus plantarum -survival, functional and potential probiotic properties in the human intestinal tract.Int. Dairy J. 16 : 1018-1028. - Piqué N, Berlanga M, Miñana-Galbis D. 2019. Health benefits of heat-killed (Tyndallized) probiotics: An overview.
Int. J. Mol. Sci. 20 : 2534. - Liu Y-W, Liong M-T, Tsai Y-C. 2018. New perspectives of
Lactobacillus plantarum as a probiotic: The gut-heart-brain axis.J. Microbiol. 56 : 601-613. - Zhang M, Jiang Y, Cai M, Yang Z. 2020. Characterization and ACE inhibitory activity of fermented milk with probiotic
Lactobacillus plantarum K25 as analyzed by GC-MS-based metabolomics approach.J. Microbiol. Biotechnol. 30 : 903-911. - Zheng J, Wittouck S, Salvetti E, Franz CM, Harris HM, Mattarelli P,
et al . 2020. A taxonomic note on the genusLactobacillus : Description of 23 novel genera, emended description of the genusLactobacillus Beijerinck 1901, and union ofLactobacillaceae andLeuconostocaceae .Int. J. Syst. Evol. Microbiol. 70 : 2782-2858. - Mukerji P, Roper JM, Stahl B, Smith AB, Burns F, Rae JC,
et al . 2016. Safety evaluation of AB-LIFE®(Lactobacillus plantarum CECT 7527, 7528 and 7529): Antibiotic resistance and 90-day repeated-dose study in rats.Food Chem. Toxicol. 92 : 117-128. - Olek A, Woynarowski M, Ahrén IL, Kierkuś J, Socha P, Larsson N,
et al . 2017. Efficacy and safety ofLactobacillus plantarum DSM 9843 (LP299V) in the prevention of antibiotic-associated gastrointestinal symptoms in children-randomized, double-blind, placebo-controlled study.J. Pediatr. 186 : 82-86. - Gotteland M, Cires MJ, Carvallo C, Vega N, Ramirez MA, Morales P,
et al . 2014. Probiotic screening and safety evaluation ofLactobacillus strains from plants, artisanal goat cheese, human stools, and breast milk.J. Med. Food 17 : 487-495. - Goel A, Halami PM, Tamang JP. 2020. Genome analysis of
Lactobacillus plantarum isolated from some Indian fermented foods for bacteriocin production and probiotic marker genes.Front. Microbiol. 11 : 40. - Park DM, Bae J-H, Kim MS, Kim H, Kang SD, Shim S,
et al . 2019. Suitability ofLactobacillus plantarum SPC-SNU 72-2 as a probiotic starter for sourdough fermentation.J. Microbiol. Biotechnol. 29 : 1729-1738. - Gueimonde M, Sánchez B, de Los Reyes-Gavilán CG, Margolles A. 2013. Antibiotic resistance in probiotic bacteria.
Front. Microbiol. 4 : 202. - Devirgiliis C, Zinno P, Perozzi G. 2013. Update on antibiotic resistance in foodborne
Lactobacillus andLactococcus species.Front. Microbiol. 4 : 301. - Abriouel H, Muñoz MdCC, Lerma LL, Montoro BP, Bockelmann W, Pichner R,
et al . 2015. New insights in antibiotic resistance ofLactobacillus species from fermented foods.Food Res. Int. 78 : 465-481. - Campedelli I, Mathur H, Salvetti E, Clarke S, Rea MC, Torriani S,
et al . 2019. Genus-wide assessment of antibiotic resistance inLactobacillus spp.Appl. Environ. Microbiol. 85 : e01738-18. - Chokesajjawatee N, Santiyanont P, Chantarasakha K, Kocharin K, Thammarongtham C, Lertampaiporn S,
et al . 2020. Safety assessment of a nham starter cultureLactobacillus plantarum BCC9546 via whole-genome analysis.Sci. Rep. 10 : 10241. - Klarin B, Larsson A, Molin G, Jeppsson B. 2019. Susceptibility to antibiotics in isolates of
Lactobacillus plantarum RAPD‐type Lp299v, harvested from antibiotic treated, critically ill patients after administration of probiotics.Microbiologyopen 8 : e00642. - Shao Y, Zhang W, Guo H, Pan L, Zhang H, Sun T. 2015. Comparative studies on antibiotic resistance in
Lactobacillus casei andLactobacillus plantarum .Food Control 50 : 250-258. - Park S-Y, Cho S-A, Kim S-H, Lim S-D. 2014. Physiological characteristics and anti-obesity effect of
Lactobacillus plantarum Q180 isolated from feces.Korean J. Food Sci. Anim. Resour. 34 : 647-655. - Kwon J, Kim B, Lee C, Joung H, Kim B-K, Choi IS,
et al . 2020. Comprehensive amelioration of high-fat diet-induced metabolic dysfunctions through activation of the PGC-1α pathway by probiotics treatment in mice.PLoS One 15 : e0228932. - Park YE, Kim MS, Shim KW, Kim Y-I, Chu J, Kim B-K,
et al . 2020. Effects ofLactobacillus plantarum Q180 on postprandial lipid levels and intestinal environment: A double-blind, randomized, placebo-controlled, parallel trial.Nutrients 12 : 255. - Baek JH, Kim KH, Moon JY, Yeo S-H, Jeon CO. 2020.
Acetobacter oryzoeni sp. nov., isolated from Korean rice wine vinegar.Int. J. Syst. Evol. Microbiol. 70 : 2026-2033. - Berlin K, Koren S, Chin C-S, Drake JP, Landolin JM, Phillippy AM. 2015. Assembling large genomes with single-molecule sequencing and locality-sensitive hashing.
Nat. Biotechnol. 33 : 623-630. - Walker BJ, Abeel T, Shea T, Priest M, Abouelliel A, Sakthikumar S,
et al . 2014. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement.PLoS One 9 : e112963. - Cole JR, Wang Q, Cardenas E, Fish J, Chai B, Farris RJ,
et al . 2009. The Ribosomal Database Project: improved alignments and new tools for rRNA analysis.Nucleic Acids Res. 37 : 141-145. - Kumar S, Stecher G, Tamura K. 2016. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets.
Mol. Biol. Evol. 33 : 1870-1874. - Na S-I, Kim YO, Yoon S-H, Ha S-m, Baek I, Chun J. 2018. UBCG: up-to-date bacterial core gene set and pipeline for phylogenomic tree reconstruction.
J. Microbiol. 56 : 280-285. - Lee I, Kim YO, Park S-C, Chun J. 2016. OrthoANI: an improved algorithm and software for calculating average nucleotide identity.
Int. J. Syst. Evol. Microbiol. 66 : 1100-1103. - Meier-Kolthoff JP, Auch AF, Klenk H-P, Göker M. 2013. Genome sequence-based species delimitation with confidence intervals and improved distance functions.
BMC Bioinformatics 14 : 60. - Grant JR, Stothard P. 2008. The CGView Server: a comparative genomics tool for circular genomes.
Nucleic Acids Res. 36 : 181-184. - Huerta-Cepas J, Szklarczyk D, Heller D, Hernández-Plaza A, Forslund SK, Cook H,
et al . 2019. eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses.Nucleic Acids Res. 47 : 309-314. - Buchfink B, Xie C, Huson DH. 2015. Fast and sensitive protein alignment using DIAMOND.
Nat. Methods 12 : 59-60. - Alcock BP, Raphenya AR, Lau TT, Tsang KK, Bouchard M, Edalatmand A,
et al . 2020. CARD 2020: antibiotic resistome surveillance with the comprehensive antibiotic resistance database.Nucleic Acids Res. 48 : 517-525. - Gupta SK, Padmanabhan BR, Diene SM, Lopez-Rojas R, Kempf M, Landraud L,
et al . 2014. ARG-ANNOT, a new bioinformatic tool to discover antibiotic resistance genes in bacterial genomes.Antimicrob. Agents Chemother. 58 : 212-220. - Bortolaia V, Kaas RS, Ruppe E, Roberts MC, Schwarz S, Cattoir V,
et al . 2020. ResFinder 4.0 for predictions of phenotypes from genotypes.J. Antimicrob. Chemother. 75 : 3491-3500. - Chen L, Yang J, Yu J, Yao Z, Sun L, Shen Y,
et al . 2005. VFDB: a reference database for bacterial virulence factors.Nucleic Acids Res. 33 : 325-328. - Arndt D, Grant JR, Marcu A, Sajed T, Pon A, Liang Y,
et al . 2016. PHASTER: a better, faster version of the PHAST phage search tool.Nucleic Acids Res. 44 : 16-21. - Siguier P, Pérochon J, Lestrade L, Mahillon J, Chandler M. 2006. ISfinder: the reference centre for bacterial insertion sequences.
Nucleic Acids Res. 34 : 32-36. - Mercanti DJ, Carminati D, Reinheimer JA, Quiberoni A. 2011. Widely distributed lysogeny in probiotic lactobacilli represents a potentially high risk for the fermentative dairy industry.
Int. J. Food Microbiol. 144 : 503-510. - Binetti AG, Del Río B, Martìn MC, Alvarez MA. 2005. Detection and characterization of
Streptococcus thermophilus bacteriophages by use of the antireceptor gene sequence.Appl. Environ. Microbiol. 71 : 6096-6103. - Clinical and laboratory standards institute. 2012. Methods for dilution antimicrobial susceptibility tests for bacteria that grow aerobically; approved standard CLSI.
In: document M07-A9 . Wayne, PA, USA. - Beecher DJ, Schoeni JL, Wong A. 1995. Enterotoxic activity of hemolysin BL from
Bacillus cereus .Infect. Immun. 63 : 4423-4428. - Kim KH, Lee SH, Chun BH, Jeong SE, Jeon CO. 2019.
Tetragenococcus halophilus MJ4 as a starter culture for repressing biogenic amine (cadaverine) formation during saeu-jeot (salted shrimp) fermentation.Food Microbiol. 82 : 465-473. - Brückner H, Flassig S, Kirschbaum J. 2012. Determination of biogenic amines in infusions of tea (
Camellia sinensis ) by HPLC after derivatization with 9-fluorenylmethoxycarbonyl chloride (Fmoc-Cl).Amino Acids 42 : 877-885. - Chun J, Oren A, Ventosa A, Christensen H, Arahal DR, da Costa MS,
et al . 2018. Proposed minimal standards for the use of genome data for the taxonomy of prokaryotes.Int. J. Syst. Evol. Microbiol. 68 : 461-466. - Domig KJ, Zycka-Krzesinska J, Bardowski J, Morelli L. 2008. Molecular assessment of erythromycin and tetracycline resistance genes in lactic acid bacteria and bifidobacteria and their relation to the phenotypic resistance.
Int. J. Probiotics Prebiotics 3 : 271-280. - Hummel AS, Hertel C, Holzapfel WH, Franz CM. 2007. Antibiotic resistances of starter and probiotic strains of lactic acid bacteria.
Appl. Environ. Microbiol. 73 : 730-739. - Mayrhofer S, Van Hoek AH, Mair C, Huys G, Aarts HJ, Kneifel W,
et al . 2010. Antibiotic susceptibility of members of theLactobacillus acidophilus group using broth microdilution and molecular identification of their resistance determinants.Int. J. Food Microbiol. 144 : 81-87. - Lahtinen SJ, Boyle RJ, Margolles A, Frias R, Gueimonde M. 2009. Safety assessment of probiotics.
In: Prebiotics probiotics Sci. Technol . Springer-Verlag press, Berlin. - Cebeci A, Gürakan C. 2003. Properties of potential probiotic
Lactobacillus plantarum strains.Food Microbiol. 20 : 511-518. - Huycke MM, Spiegel CA, Gilmore MS. 1991. Bacteremia caused by hemolytic, high-level gentamicin-resistant
Enterococcus faecalis .Antimicrob. Agents Chemother. 35 : 1626-1634. - Sharma P, Tomar SK, Goswami P, Sangwan V, Singh R. 2014. Antibiotic resistance among commercially available probiotics.
Food Res. Int. 57 : 176-195. - Flórez AB, Egervärn M, Danielsen M, Tosi L, Morelli L, Lindgren S,
et al . 2006. Susceptibility ofLactobacillus plantarum strains to six antibiotics and definition of new susceptibility-resistance cutoff values.Microb. Drug Resist. 12 : 252-256. - Zhang F, Gao J, Wang B, Huo D, Wang Z, Zhang J,
et al . 2018. Whole-genome sequencing reveals the mechanisms for evolution of streptomycin resistance inLactobacillus plantarum .J. Dairy Sci. 101 : 2867-2874. - Facklam RR, Carvalho MdGS, Teixeira LM. 2002. History, taxonomy, biochemical characteristics, and antibiotic susceptibility testing of enterococci.
In: The enterococci: pathogenesis, molecular biology, and antibiotic resistance . ASM press, Washington, D.C., USA. - Zoletti GO, Pereira EM, Schuenck RP, Teixeira LM, Siqueira JF Jr, dos Santos KRN. 2011. Characterization of virulence factors and clonal diversity of
Enterococcus faecalis isolates from treated dental root canals.Res. Microbiol. 162 : 151-158. - Abouloifa H, Rokni Y, Bellaouchi R, Ghabbour N, Karboune S, Brasca M,
et al . 2019. Characterization of probiotic properties of antifungalLactobacillus strains isolated from traditional fermenting green olives.Probiotics Antimicrob. Proteins 12 : 683-696. - Kenfack CHM, Ngoufack FZ, Kaktcham PM, Wang YR, Zhu T, Yin L. 2018. Safety and antioxidant properties of five probiotic
Lactobacillus plantarum strains isolated from the digestive tract of honey bees.Am. J. Microbiol. Res. 6 : 1-8. - Wójcik W, Łukasiewicz M, Puppel K. 2021. Biogenic amines: formation, action and toxicity-a review.
J. Sci. Food Agric. 101 : 2634-2640. - Kim S, Huang E, Park S, Holzapfel W, Lim S-D. 2018. Physiological characteristics and anti-obesity effect of
Lactobacillus plantarum K10.Korean J. Food Sci. Anim. Resour. 38 : 554-569. - Gao Y, Liu Y, Sun M, Zhang H, Mu G, Tuo Y. 2020. Physiological function analysis of
Lactobacillus plantarum Y44 based on genotypic and phenotypic characteristics.J. Dairy Sci. 103 : 5916-5930. - Priyanka V, Ramesha A, Gayathri D, Vasudha M. 2020. Molecular characterization of non-biogenic amines producing
Lactobacillus plantarum GP11 isolated from traditional pickles using HRESI-MS analysis.J. Food Sci. Technol. 58 : 2216-2226.