Advertisement
Research Short Communication| Volume 102, ISSUE 12, P10838-10844, December 2019

Short communication: Complete genome sequence of Lactobacillus plantarum J26, a probiotic strain with immunomodulatory activity

  • Author Footnotes
    * These authors contributed equally to this work.
    Ziwei Zhang
    Footnotes
    * These authors contributed equally to this work.
    Affiliations
    Key Laboratory of Dairy Science, Ministry of Education, Department of Food Science, Northeast Agricultural University, Harbin, 150030, China
    Search for articles by this author
  • Author Footnotes
    * These authors contributed equally to this work.
    Chaoxin Man
    Footnotes
    * These authors contributed equally to this work.
    Affiliations
    Key Laboratory of Dairy Science, Ministry of Education, Department of Food Science, Northeast Agricultural University, Harbin, 150030, China
    Search for articles by this author
  • Linlin Sun
    Affiliations
    Key Laboratory of Dairy Science, Ministry of Education, Department of Food Science, Northeast Agricultural University, Harbin, 150030, China
    Search for articles by this author
  • Xinyan Yang
    Affiliations
    Key Laboratory of Dairy Science, Ministry of Education, Department of Food Science, Northeast Agricultural University, Harbin, 150030, China
    Search for articles by this author
  • Mingyu Li
    Affiliations
    Key Laboratory of Dairy Science, Ministry of Education, Department of Food Science, Northeast Agricultural University, Harbin, 150030, China
    Search for articles by this author
  • Wei Zhang
    Affiliations
    Key Laboratory of Dairy Science, Ministry of Education, Department of Food Science, Northeast Agricultural University, Harbin, 150030, China
    Search for articles by this author
  • Yujun Jiang
    Correspondence
    Corresponding author
    Affiliations
    Key Laboratory of Dairy Science, Ministry of Education, Department of Food Science, Northeast Agricultural University, Harbin, 150030, China
    Search for articles by this author
  • Author Footnotes
    * These authors contributed equally to this work.
Open ArchivePublished:September 20, 2019DOI:https://doi.org/10.3168/jds.2019-16593

      ABSTRACT

      Lactobacillus plantarum J26, a significant probiotic isolated from Chinese traditional fermented dairy products, exerts a positive immunomodulatory effect by regulating the expression of immune-related genes. We investigated expression of the cytokines IL-1α, IL-1β, IL-6, and tumor necrosis factor-α in the intestinal tract of mice stimulated by L. plantarum J26. In vivo, these cytokines were upregulated, peaked on d 5, and then decreased to the control level, indicating that L. plantarum J26 could induce expression of the genes encoding these proinflammatory cytokines. Teichoic acids produced by L. plantarum are recognized as key immunomodulatory molecules involved in the regulation of the host immune response. To better understand the genetic basis of this immunomodulatory mechanism, we sequenced and analyzed the whole genome of L. plantarum J26. The genome of L. plantarum J26 contains a circular chromosome and 4 circular plasmids. Lactobacillus plantarum J26 was predicted to synthesize ribitol-type backbones of wall teichoic acid. Furthermore, orthologous average nucleotide identity (OrthoANI) values showed that the genome was highly similar (>98.00%) to other L. plantarum strains, especially to L. plantarum ST-III and JDM1. The genomic data of L. plantarum J26 provide a genetic basis to further elucidate its mechanism of immunoregulation and will facilitate its application in the functional dairy food industry.

      Key words

      Short Communication

      Lactobacillus constitutes a large and diverse genus within the lactic acid bacteria and its genetic diversity is larger than that of a normal family (
      • Sun Z.
      • Harris H.M.
      • McCann A.
      • Guo C.
      • Argimon S.
      • Zhang W.
      • Yang X.
      • Jeffery I.B.
      • Cooney J.C.
      • Kagawa T.F.
      • Liu W.
      • Song Y.
      • Salvetti E.
      • Wrobel A.
      • Rasinkangas P.
      • Parkhill J.
      • Rea M.C.
      • O'Sullivan O.
      • Ritari J.
      • Douillard F.P.
      • Paul Ross R.
      • Yang R.
      • Briner A.E.
      • Felis G.E.
      • de Vos W.M.
      • Barrangou R.
      • Klaenhammer T.R.
      • Caufield P.W.
      • Cui Y.
      • Zhang H.
      • O'Toole P.W.
      Expanding the biotechnology potential of lactobacilli through comparative genomics of 213 strains and associated genera.
      ). Lactobacillus plantarum is a gram-positive bacterium (
      • Hugenholtz P.
      The Genera of Lactic Acid Bacteria.
      ) that inhabits relatively abundant ecological niches and plays a significant role in food microbiology and human nutrition because of its fermentative and probiotic functions (
      • Salvetti E.
      • Torriani S.
      • Felis G.E.
      The genus Lactobacillus: A taxonomic update.
      ). Recently, numerous studies have investigated the health benefits of L. plantarum, such as effects on oxidative stress regulation (
      • Li L.
      • Jiang Y.J.
      • Yang X.Y.
      • Liu Y.
      • Wang J.Y.
      • Man C.X.
      Immunoregulatory effects on Caco-2 cells and mice of exopolysaccharides isolated from Lactobacillus acidophilus NCFM.
      ), reduction of cholesterol level (
      • Barreto F.M.
      • Colado Simao A.N.
      • Morimoto H.K.
      • Batisti Lozovoy M.A.
      • Dichi I.
      • Helena da Silva Miglioranza L.
      Beneficial effects of Lactobacillus plantarum on glycemia and homocysteine levels in postmenopausal women with metabolic syndrome.
      ), and management of bacterial composition in feces (
      • Goossens D.
      • Jonkers D.
      • Russel M.
      • Stobberingh E.
      • Van Den Bogaard A.
      • Stockbrügger R.
      The effect of Lactobacillus plantarum 299 v on the bacterial composition and metabolic activity in faeces of healthy volunteers: A placebo-controlled study on the onset and duration of effects.
      ). Among these properties, the potential of L. plantarum to inhibit inflammation has received extensive attention. Lactobacillus plantarum GB-LP2 has also shown antiviral effects against the influenza virus in mice (
      • Choi S.-W.
      • Youn H.-N.
      • Hong W.
      • Park J.-K.
      • Yuk S.-S.
      • Kwon J.-H.
      • Noh J.-Y.
      • Kang J.-S.
      • Cho K.-J.
      • Ryu J.-J.
      • Lee J.-B.
      • Park S.-Y.
      • Choi I.-S.
      • Lee S.-W.
      • Song C.-S.
      Intranasal administration model for evaluating protection against influenza virus in mice.
      ). Similarly, L. plantarum PS128 reduced LPS-induced proinflammatory cytokine production in a mouse macrophage cell model (
      • Liu W.H.
      • Yang C.H.
      • Lin C.T.
      • Li S.W.
      • Cheng W.S.
      • Jiang Y.P.
      • Wu C.C.
      • Chang C.H.
      • Tsai Y.C.
      Genome architecture of Lactobacillus plantarum PS128, a probiotic strain with potential immunomodulatory activity.
      ).
      Teichoic acid (TA), the major component of cell walls of gram-positive bacteria, is recognized as key immunomodulatory molecule involved in the regulation of host immune response. For example, heat-killed L. plantarum MYL26 and bacterial cell wall extracts are able to reduce LPS-induced inflammation by impairing toll-like receptor 4 (TLR4)-nuclear factor κB (NF-κB) signal transduction in vitro, suggesting that constituents of bacterial cell wall help attenuate inflammation (
      • Chiu Y.-H.
      • Lu Y.-C.
      • Ou C.-C.
      • Lin S.-L.
      • Tsai C.-C.
      • Huang C.-T.
      • Lin M.-Y.
      Lactobacillus plantarum MYL26 induces endotoxin tolerance phenotype in Caco-2 cells.
      ). In addition, an investigation found that d-alanylation of TA abolished the pro- and anti-inflammatory response in vivo (
      • Smelt M.J.
      • de Haan B.J.
      • Bron P.A.
      • van Swam I.
      • Meijerink M.
      • Wells J.M.
      • Kleerebezem M.
      • Faas M.M.
      • de Vos P.
      The impact of Lactobacillus plantarum WCFS1 teichoic acid D-alanylation on the generation of effector and regulatory T-cells in healthy mice.
      ). However, immunomodulatory capability varies in different strains, probably because the types and structures of TA are significantly different (
      • Han S.H.
      • Kim J.H.
      • Martin M.
      • Michalek S.M.
      • Nahm M.H.
      Pneumococcal lipoteichoic acid (lta) is not as potent as staphylococcal lta in stimulating toll-like receptor 2.
      ).
      Lactobacillus plantarum J26 is a common probiotic strain isolated from traditional fermented dairy products. Our recent investigation showed that strain J26 could alleviate oxidative stress by modulating the production of antioxidant enzymes significantly in vitro (
      • Hou Y.
      • Li X.
      • Liu X.
      • Zhang Y.
      • Zhang W.
      • Man C.
      • Jiang Y.
      Transcriptomic responses of Caco-2 cells to Lactobacillus rhamnosus GG and Lactobacillus plantarum J26 against oxidative stress.
      ). In the current study, we confirmed that the strain exhibits potent immunomodulatory activity. To better understand the immunomodulatory mechanism of L. plantarum J26 and obtain detailed insight into the genetic basis, we analyzed the whole genome sequence of L. plantarum J26 combined with biochemical assays.
      Lactobacillus plantarum J26 was isolated from Chinese traditional fermented dairy products in Inner Mongolia (China). It was grown in de Man, Rogosa, and Sharpe (Qingdao Hope Bio-Technology Co. Ltd., Qingdao, China) broth at 37°C for 16 h through 2 propagation steps. Genomic DNA was extracted using the Qiagen DNA Mini kit (Qiagen, Hilden, Germany) following the manufacturer's instructions.
      Quantity and quality determinations of genomic DNA in L. plantarum J26 were conducted by Qubit 2.0 (Thermo Fisher Scientific, Waltham, MA) and Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA). The genomic DNA was sequenced by Illumina Miseq (insert size of 400 bp; Illumina Inc., San Diego, CA) with paired-end sequencing mode and Pacbio sequencing (20,000 bp template library; Pacific Biosciences, Menlo Park, CA) with standard sequencing mode using next-generation sequencing technology and single molecule real-time (SMRT) technology. A total of 3,221,902 reads were obtained. Genome assembly was performed using A5-miseq v20150522 (
      • Coil D.
      • Jospin G.
      • Darling A.E.
      A5-miseq: An updated pipeline to assemble microbial genomes from Illumina MiSeq data.
      ) and CANU (
      • Koren S.
      • Walenz B.P.
      • Berlin K.
      • Miller J.R.
      • Bergman N.H.
      • Phillippy A.M.
      Canu: Scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation.
      ) software for the data obtained by second- and third-generation sequencing platforms. Then, the assembled results were confirmed by MUMmer software (http://mummer.sourceforge.net/) to determine the positional relationship between contigs and to fill gaps (
      • Delcher A.L.
      • Salzberg S.L.
      • Phillippy A.M.
      Using MUMmer to identify similar regions in large sequence sets.
      ). Finally, a complete genome without gaps was constructed after the correction using Pilon software (
      • Walker B.J.
      • Abeel T.
      • Shea T.
      • Priest M.
      • Abouelliel A.
      • Sakthikumar S.
      • Cuomo C.A.
      • Zeng Q.
      • Wortman J.
      • Young S.K.
      • Earl A.M.
      Pilon: An integrated tool for comprehensive microbial variant detection and genome assembly improvement.
      ). Gene prediction was carried out by GeneMarkS (http://topaz.gatech.edu/). The predictions of transfer (t)RNA, rRNA, and small (s)RNA were identified using tRNAscan-SE (
      • Lowe T.M.
      • Eddy S.R.
      tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence.
      ), rRNAmmer (
      • Lagesen K.
      • Hallin P.
      • Rodland E.A.
      • Staerfeldt H.H.
      • Rognes T.
      • Ussery D.W.
      RNAmmer: Consistent and rapid annotation of ribosomal RNA genes.
      ), and Rfam (
      • Gardner P.P.
      • Daub J.
      • Tate J.G.
      • Nawrocki E.P.
      • Kolbe D.L.
      • Lindgreen S.
      • Wilkinson A.C.
      • Finn R.D.
      • Griffiths-Jones S.
      • Eddy S.R.
      • Bateman A.
      Rfam: Updates to the RNA families database.
      ) softwares, respectively. Genome annotation was accomplished by National Center for Biotechnology Information (NCBI) Prokaryotic Genome Annotation Pipeline (
      • Pruitt K.D.
      • Tatusova T.
      • Brown G.R.
      • Maglott D.R.
      NCBI Reference Sequences (RefSeq): Current status, new features and genome annotation policy.
      ) and confirmed by BLAST analysis (www.ncbi.nlm.nih.gov/protein). Using the multiple sequence alignment software MUSCLE (https://www.ebi.ac.uk/Tools/msa/muscle/) to align the core gene sequences of the L. plantarum J26 and other species, we constructed a phylogenetic tree based on the core genes (
      • Kumar S.
      • Stecher G.
      • Tamura K.
      MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets.
      ). Furthermore, the average nucleotide identity of the genome sequence between the L. plantarum J26 and other reference strains was evaluated by orthologous average nucleotide identity (OrthoANI; http://www.ezbiocloud.net/tools/orthoani) tool.
      The expression of immune-related genes was evaluated in vivo. All of the procedures used here were approved by the Institutional Animal Care and Use Committee of Northeast Agricultural University in China, and the animal experiments were performed following the Northeast Agricultural University guidelines for Laboratory Animals Care and Use.
      Forty BALB/cl mice, 8 wk old and 18 to 22 g of BW, were provided by Heilongjiang Province Medical University Laboratory Animal Center (Harbin, China). Mice were randomly divided into 2 groups and acclimated in pathogen-free cages under controlled conditions (22°C ± 2°C, relative humidity 55% ± 5%, 12-h light/dark cycles) with free access to water and a standard diet for 7 d. After this period, the experimental group was given 0.3 mL of skim milk containing L. plantarum J26 (108 cfu/mL) each day. The mice in the control group were given skim milk without L. plantarum J26. The mice were killed at 1, 3, 5, or 7 d after supplementation of bacteria. The colon and ileum tissue were removed, washed with cold saline, and placed in liquid nitrogen (
      • Jiang Y.
      • Li L.
      • Sun H.
      • Shan Y.
      • Liu Y.
      • Li L.
      • Qu B.
      • Man C.
      Induction of cytokines via NF-κB and p38 MAP kinase signalling pathways associated with the immunomodulation by Lactobacillus plantarum NDC 75017 in vitro and in vivo..
      ). To evaluate the mRNA expression levels of IL1A, IL1B, IL6, and TNFA, RNA of colon and ileum was extracted using a commercial kit (RNAprep Tissue Kit, Tiangen Biotech, Shanghai, China) and converted to cDNA using a PrimeScript RT reagent kit (Takara, Dalian, China). The expression of immune-related genes was analyzed by using real-time quantitative PCR (ABI Prism 7500 system; Applied Biosystems, Foster City, CA). Approximately 2 μL of the RNA sample was added to the 20-μL total reaction volume according to the requirements of the SYBR PremicExTaqTMII (PerfectReal-time) kit (Takara). The GAPDH gene was used as the reference gene. To amplify the selected immune-related genes IL1A, IL1B, IL6, and TNFA and GAPDH, specific primers were designed using Primer 5.0 software (Premier Biosoft, USA; http://www.premierbiosoft.com/primerdesign/index.html) as shown in Table 1. The protocol for the amplification reactions was as follows: 95°C for 30 s (denaturation) and 40 cycles of 95°C for 5 s and 60°C for 34 s (amplification and quantification). The results were expressed as relative values after normalization to GAPDH mRNA (
      • Jiang Y.
      • Li L.
      • Sun H.
      • Shan Y.
      • Liu Y.
      • Li L.
      • Qu B.
      • Man C.
      Induction of cytokines via NF-κB and p38 MAP kinase signalling pathways associated with the immunomodulation by Lactobacillus plantarum NDC 75017 in vitro and in vivo..
      ).
      Table 1The base sequence, annealing temperatures, and sizes of PCR products for target gene-specific primers
      Target genePrimer sequence
      F = forward; R = reverse.
      (5′→3′)
      Annealing temperature (°C)Product size (bp)
      GAPDHF: GCCTGGAGAAACCTGCC55.4200
      R: ATACCAGGAAATGAGCTTGACA56.7
      IL1AF: TCTGCCATTGACCATCTC50.7183
      R: AATCTTCCCGTTGCTTG50.8
      IL1BF: AAGTTGACGGACCCCA50.6126
      R: GTGATACTGCCTGCCTGA51.3
      IL6F: GAAGTGATTCTTACGCA42.5247
      R: GTTTAGGTGGAGGTGTC43.0
      TNFAF: CGTAGCAAACCACCAAG49.3149
      R: CCGTGAAGAGAACCTGG50.4
      1 F = forward; R = reverse.
      The complete genome of L. plantarum J26 comprises a circular chromosome (3,096,468 bp) and 4 circular plasmids: pJ26p1 (83,889 bp), pJ26p2 (43,153 bp), pJ26p3 (35,191 bp), and pJ26p4 (8,686 bp), with GC contents of 44.79, 41.14, 40.30, 45.48, and 35.92%, respectively (Figure 1). Among the 3,289 identified protein-coding genes, 3,068 were in the chromosome and 105, 62, 44, and 10 in plasmids pJ26p1 to pJ26p4, respectively. Furthermore, the chromosome included 16 rRNA operons, 71 tRNA, and 1 sRNA (Table 2). Each gene was assigned to a functional category according to the Clusters of Orthologous Groups database (www.ncbi.nlm.nih.gov/COG; Table 3). Among 23 functional categories, the most abundant were assigned to transcription (248 genes); 236 genes were assigned to carbohydrate transport and metabolism; and 137 genes were responsible for cell wall, membrane, and envelope biogenesis.
      Figure thumbnail gr1
      Figure 1Circular genome map of Lactobacillus plantarum J26. From periphery to center: genome sequence (ring 1 and 2), Cluster of Orthologous Groups of proteins (COG) annotated coding sequences (ring 3 and 4), intact prophages in red and clustered regularly interspaced short palindromic repeats (CRISPR) in black (ring 5), noncoding (nc)RNA genes (ring 6), GC content (ring 7), and GC skew [(G+C)/(G+C)] (ring 8), where values >0 are in green, and values <0 are in purple. The smaller maps show the genome map of plasmids in L. plantarum J26. The image was created using the software Circos (http://circos.ca/software/download/circos/).
      Table 2Features of Lactobacillus plantarum J26 genome
      Feature
      CDS = coding sequence; tRNA = transfer RNA; sRNA = small RNA sequences.
      ChromosomePlasmid 1Plasmid 2Plasmid 3Plasmid 4
      Size (bp)3,096,46883,88943,15335,1918,686
      GC content (%)44.7941.1440.3045.4835.92
      Predicted genes3,068105624410
      Protein-coding genes (CDS)3,068105624410
      rRNA operons160000
      tRNA710000
      sRNA12111
      1 CDS = coding sequence; tRNA = transfer RNA; sRNA = small RNA sequences.
      Table 3Cluster of Orthologous Groups of proteins (COG) categories of coding proteins in Lactobacillus plantarum J26
      COG classNameCountProportion (%)
      CEnergy production and conversion1114.21
      DCell cycle control, cell division, chromosome partitioning361.37
      EAmino acid transport and metabolism2198.31
      FNucleotide transport and metabolism1023.87
      GCarbohydrate transport and metabolism2368.96
      HCoenzyme transport and metabolism1264.78
      ILipid transport and metabolism953.61
      JTranslation, ribosomal structure and biogenesis2037.71
      KTranscription2489.42
      LReplication, recombination and repair1144.33
      MCell wall/membrane/envelope biogenesis1375.20
      NCell motility170.65
      OPosttranslational modification, protein, turnover, chaperones863.26
      PInorganic ion transport and metabolism1294.90
      QSecondary metabolites biosynthesis, transport and catabolism281.06
      RGeneral function prediction only2379.00
      SFunction unknown1616.11
      TSignal transduction mechanisms1184.48
      UIntracellular trafficking, secretion, and vesicular transport180.68
      VDefense mechanisms762.89
      WExtracellular structures30.11
      XMobilome: prophages, transposons1345.09
      A phylogenetic tree was constructed based on homology of core genes of L. plantarum J26 with related strains. Phylogenetic tree analysis (Figure 2) showed that L. plantarum J26 was more closely related to L. plantarum JDM1 than to other L. plantarum strains and evolved from the common ancestor compared with ST-III (which contains TA synthesis protein B instead of TA synthesis proteins TagF1 and TagF2 (
      • Liu W.H.
      • Yang C.H.
      • Lin C.T.
      • Li S.W.
      • Cheng W.S.
      • Jiang Y.P.
      • Wu C.C.
      • Chang C.H.
      • Tsai Y.C.
      Genome architecture of Lactobacillus plantarum PS128, a probiotic strain with potential immunomodulatory activity.
      ). The genome of J26 was compared with 10 other L. plantarum strains. The OrthoANI values (Table 4) indicated that the genome of L. plantarum J26 was closest to L. plantarum JDM1 (99.30% OrthoANI), followed by L. plantarum 5-2 (98.73%) and L. plantarum WCFS1 (98.71%). Generally, two genomes are considered the same species when the ANI value is higher than 95 to 96% (
      • Lee I.
      • Ouk Kim Y.
      • Park S.C.
      • Chun J.
      OrthoANI: An improved algorithm and software for calculating average nucleotide identity.
      ). Consequently, L. plantarum J26 was confirmed to belong to the species of L. plantarum.
      Figure thumbnail gr2
      Figure 2Phylogenetic tree based on core gene sequence showing the phylogenetic relationships of 11 Lactobacillus plantarum strains. Phylogenetic analysis was evaluated by MUSCLE software (https://www.ebi.ac.uk/Tools/msa/muscle/).
      Table 4Average nucleotide identity (ANI; %) values between the strains and different Lactobacillus plantarum species
      Strain code165-2JDM1LPL-1LZ206LZ227LZ95WCFS1J26ST-IIIZJ316
      16100.0098.7198.9099.0799.1298.8898.8098.7698.5398.7399.11
      5-298.71100.0098.9498.7398.7698.7199.0799.0098.7399.0098.73
      JDM198.9098.94100.0098.8798.7298.8398.9698.9599.3098.8798.75
      LPL-199.0798.7398.87100.0098.7098.8098.7498.7998.7098.6798.70
      LZ20699.1298.7698.7298.70100.0099.2798.7398.6998.5198.6699.71
      LZ22798.8898.7198.8398.8099.27100.0098.6798.6698.4098.6699.31
      LZ9598.8099.0798.9698.7498.7398.66100.0098.9398.7199.8998.83
      WCFS198.7699.0098.9598.7998.6998.6698.93100.0098.7198.9998.77
      J2698.5398.7399.3098.7098.5198.4098.7198.71100.0098.6298.56
      ST-III98.7399.0098.8798.6798.6698.6699.8998.9998.62100.0098.73
      ZJ31699.1198.7598.7598.7099.7199.3198.8398.7798.5698.73100.00
      Intestinal epithelial cells stimulated by microbes can secrete cytokines and pattern recognition receptors to participate in appropriate immune regulation (
      • Round J.L.
      • Mazmanian S.K.
      The gut microbiota shapes intestinal immune responses during health and disease.
      ). As shown in Figure 3, we investigated mRNA expression of proinflammatory cytokines IL-1α, IL-1β, IL-6, and TNF-α in the intestinal tract of mice stimulated with L. plantarum J26. These cytokines were all upregulated, with expression peaking at d 5, and then decreasing to control levels. Moreover, the expression of IL1A, IL1B, IL6, and TNFA was significantly different from that in the control group (P < 0.01) on d 5, indicating that L. plantarum J26 was able to trigger the mRNA expression of the genes encoding these proinflammatory cytokines in vivo, results consistent with our previous study with another L. plantarum strain (
      • Jiang Y.
      • Li L.
      • Sun H.
      • Shan Y.
      • Liu Y.
      • Li L.
      • Qu B.
      • Man C.
      Induction of cytokines via NF-κB and p38 MAP kinase signalling pathways associated with the immunomodulation by Lactobacillus plantarum NDC 75017 in vitro and in vivo..
      ). Although IL-1α, IL-1β, IL-6, and TNF-α are considered proinflammatory cytokines, the specific mechanism by which L. plantarum J26 regulates their expression is unknown. We found no evidence that the immunomodulatory activity of L. plantarum J26 was directly related to certain host cells; our results indicate that the immune responses induced by J26 were likely linked to the type of immune-related cytokines, which is the focus of our further research.
      Figure thumbnail gr3
      Figure 3Expression of IL1A, IL1B, IL6, and TNFA genes induced by Lactobacillus plantarum J26 in mice. Mice were intragastrically administered 0.3 mL of L. plantarum J26 (gray bars) or an equal amount of sterilized skim milk (white bars) for 7 d, and the expression of IL1A (a), IL1B (b), IL6 (c), and TNFA (d) was determined. Results are expressed as mean ± SD of 3 independent determinations. *P < 0.05, **P < 0.01 compared with control.
      The immunomodulatory mechanism of J26 was further investigated by analysis of the whole genome of L. plantarum J26. Teichoic acids consist of 2 types: wall teichoic acid (WTA) and lipoteichoic acid (LTA). For L. plantarum WCFS1, d-alanylated TA in the cell walls account for several host immunomodulatory effects (
      • Smelt M.J.
      • de Haan B.J.
      • Bron P.A.
      • van Swam I.
      • Meijerink M.
      • Wells J.M.
      • Kleerebezem M.
      • Faas M.M.
      • de Vos P.
      The impact of Lactobacillus plantarum WCFS1 teichoic acid D-alanylation on the generation of effector and regulatory T-cells in healthy mice.
      ). Because TA are recognized as immunoregulatory molecules, the TA-related genes of L. plantarum J26 have been searched by in silico analysis. The synthesis genes encoding the WTA backbones are tag and tar homologs, whereas those of LTA are ltaS. The genes encoding LTA synthase (ltaS), TA glycosylation proteins (gtcA1, gtcA2, gtcA3), and d-alanylation protein (dltX) were found in the L. plantarum J26 genome. Two gene regions, tagD1-tagF2 (tagD1, tagF1, and tagF2) and tarI-tarL (tarI, tarJ, tarK and tarL), were regarded as the synthesis genes of glycerol-type and ribitol-type backbones of WTA. In L. plantarum J26, only a glycerol-3-phosphate cytidylyltransferase gene tagD1 was detected. The genes encoding TA synthesis protein F (tagF1 and tagF2) were absent. However, there are 2 TA synthesis protein B genes (tagB1 and tagB2) instead, and tagB1 and tagB2 found in L. plantarum were supposed as tarK and tarL. The tarI and tarJ genes were annotated as d-ribitol-5-phosphate cytidylyltransferase and ribitol-5-phosphate 2-dehydrogenase. Hence, we propose that L. plantarum J26 contains tarI-tarL region and synthesizes ribitol-type backbones of WTA. Besides, 2 TA transporters (permease protein, tag; and ATP-binding protein, tagH), subunits of the ABC transporter complex, were identified and may be associated with the immunomodulatory effects of L. plantarum J26.
      In the current study, we demonstrated that the L. plantarum J26 is a probiotic strain with immunomodulatory activity (regulating the expression of immune-related genes). The genome data of L. plantarum J26 provide substantial information for understanding the genetic basis of its probiotic properties. In addition, we discovered genes involved in key immunomodulatory molecules (teichoic acids) in the genome of J26. These results contribute to the development of immune-related functional dairy foods containing lactic acid bacteria. Further investigation into the molecular mechanism by which the TA-related gene products affect the host immune response are ongoing. The complete genomic sequence of L. plantarum J26 has been deposited in National Center for Biotechnology Information under GenBank accession nos. CP033616, CP033617, CP033618, CP033619, and CP033620.

      ACKNOWLEDGMENTS

      This work was supported by the National Natural Science Foundation of China (Beijing; project number 31871828) and the National High-Level Talents Special Support Program of China (810098).

      REFERENCES

        • Barreto F.M.
        • Colado Simao A.N.
        • Morimoto H.K.
        • Batisti Lozovoy M.A.
        • Dichi I.
        • Helena da Silva Miglioranza L.
        Beneficial effects of Lactobacillus plantarum on glycemia and homocysteine levels in postmenopausal women with metabolic syndrome.
        Nutrition. 2014; 30 (24613434): 939-942
        • Chiu Y.-H.
        • Lu Y.-C.
        • Ou C.-C.
        • Lin S.-L.
        • Tsai C.-C.
        • Huang C.-T.
        • Lin M.-Y.
        Lactobacillus plantarum MYL26 induces endotoxin tolerance phenotype in Caco-2 cells.
        BMC Microbiol. 2013; 13 (23937116): 190
        • Choi S.-W.
        • Youn H.-N.
        • Hong W.
        • Park J.-K.
        • Yuk S.-S.
        • Kwon J.-H.
        • Noh J.-Y.
        • Kang J.-S.
        • Cho K.-J.
        • Ryu J.-J.
        • Lee J.-B.
        • Park S.-Y.
        • Choi I.-S.
        • Lee S.-W.
        • Song C.-S.
        Intranasal administration model for evaluating protection against influenza virus in mice.
        J. Bacteriol. Virol. 2015; 45: 44
        • Coil D.
        • Jospin G.
        • Darling A.E.
        A5-miseq: An updated pipeline to assemble microbial genomes from Illumina MiSeq data.
        Bioinformatics. 2015; 31 (25338718): 587-589
        • Delcher A.L.
        • Salzberg S.L.
        • Phillippy A.M.
        Using MUMmer to identify similar regions in large sequence sets.
        Curr. Protoc. Bioinformatics. 2003; (10.1002/0471250953.bi1003s00.): 10.3.1-10.3.18
        • Gardner P.P.
        • Daub J.
        • Tate J.G.
        • Nawrocki E.P.
        • Kolbe D.L.
        • Lindgreen S.
        • Wilkinson A.C.
        • Finn R.D.
        • Griffiths-Jones S.
        • Eddy S.R.
        • Bateman A.
        Rfam: Updates to the RNA families database.
        Nucleic Acids Res. 2009; 37: D136-D140
        • Goossens D.
        • Jonkers D.
        • Russel M.
        • Stobberingh E.
        • Van Den Bogaard A.
        • Stockbrügger R.
        The effect of Lactobacillus plantarum 299 v on the bacterial composition and metabolic activity in faeces of healthy volunteers: A placebo-controlled study on the onset and duration of effects.
        Aliment. Pharmacol. Ther. 2003; 18 (12950422): 495-505
        • Han S.H.
        • Kim J.H.
        • Martin M.
        • Michalek S.M.
        • Nahm M.H.
        Pneumococcal lipoteichoic acid (lta) is not as potent as staphylococcal lta in stimulating toll-like receptor 2.
        Infect. Immun. 2003; 71 (14500472): 5541-5548
        • Hou Y.
        • Li X.
        • Liu X.
        • Zhang Y.
        • Zhang W.
        • Man C.
        • Jiang Y.
        Transcriptomic responses of Caco-2 cells to Lactobacillus rhamnosus GG and Lactobacillus plantarum J26 against oxidative stress.
        J. Dairy Sci. 2019; 102 (31255276): 7684-7696
        • Hugenholtz P.
        The Genera of Lactic Acid Bacteria.
        Blackie Academic and Professional, London, UK1998
        • Jiang Y.
        • Li L.
        • Sun H.
        • Shan Y.
        • Liu Y.
        • Li L.
        • Qu B.
        • Man C.
        Induction of cytokines via NF-κB and p38 MAP kinase signalling pathways associated with the immunomodulation by Lactobacillus plantarum NDC 75017 in vitro and in vivo..
        J. Funct. Foods. 2016; 20: 215-225
        • Koren S.
        • Walenz B.P.
        • Berlin K.
        • Miller J.R.
        • Bergman N.H.
        • Phillippy A.M.
        Canu: Scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation.
        Genome Res. 2017; 27: 722-736
        • Kumar S.
        • Stecher G.
        • Tamura K.
        MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets.
        Mol. Biol. Evol. 2016; 33 (27004904): 1870-1874
        • Lagesen K.
        • Hallin P.
        • Rodland E.A.
        • Staerfeldt H.H.
        • Rognes T.
        • Ussery D.W.
        RNAmmer: Consistent and rapid annotation of ribosomal RNA genes.
        Nucleic Acids Res. 2007; 35: 3100-3108
        • Lee I.
        • Ouk Kim Y.
        • Park S.C.
        • Chun J.
        OrthoANI: An improved algorithm and software for calculating average nucleotide identity.
        Int. J. Syst. Evol. Microbiol. 2016; 66 (26585518): 1100-1103
        • Li L.
        • Jiang Y.J.
        • Yang X.Y.
        • Liu Y.
        • Wang J.Y.
        • Man C.X.
        Immunoregulatory effects on Caco-2 cells and mice of exopolysaccharides isolated from Lactobacillus acidophilus NCFM.
        Food Funct. 2014; 5 (25340590): 3261-3268
        • Liu W.H.
        • Yang C.H.
        • Lin C.T.
        • Li S.W.
        • Cheng W.S.
        • Jiang Y.P.
        • Wu C.C.
        • Chang C.H.
        • Tsai Y.C.
        Genome architecture of Lactobacillus plantarum PS128, a probiotic strain with potential immunomodulatory activity.
        Gut Pathog. 2015; 7 (26279684): 22
        • Lowe T.M.
        • Eddy S.R.
        tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence.
        Nucleic Acids Res. 1997; 25: 955-964
        • Pruitt K.D.
        • Tatusova T.
        • Brown G.R.
        • Maglott D.R.
        NCBI Reference Sequences (RefSeq): Current status, new features and genome annotation policy.
        Nucleic Acids Res. 2012; 40 (22121212): D130-D135
        • Round J.L.
        • Mazmanian S.K.
        The gut microbiota shapes intestinal immune responses during health and disease.
        Nat. Rev. Immunol. 2009; 9 (19343057): 313-323
        • Salvetti E.
        • Torriani S.
        • Felis G.E.
        The genus Lactobacillus: A taxonomic update.
        Probiotics Antimicrob. Proteins. 2012; 4 (26782181): 217-226
        • Smelt M.J.
        • de Haan B.J.
        • Bron P.A.
        • van Swam I.
        • Meijerink M.
        • Wells J.M.
        • Kleerebezem M.
        • Faas M.M.
        • de Vos P.
        The impact of Lactobacillus plantarum WCFS1 teichoic acid D-alanylation on the generation of effector and regulatory T-cells in healthy mice.
        PLoS One. 2013; 8 (23646181)e63099
        • Sun Z.
        • Harris H.M.
        • McCann A.
        • Guo C.
        • Argimon S.
        • Zhang W.
        • Yang X.
        • Jeffery I.B.
        • Cooney J.C.
        • Kagawa T.F.
        • Liu W.
        • Song Y.
        • Salvetti E.
        • Wrobel A.
        • Rasinkangas P.
        • Parkhill J.
        • Rea M.C.
        • O'Sullivan O.
        • Ritari J.
        • Douillard F.P.
        • Paul Ross R.
        • Yang R.
        • Briner A.E.
        • Felis G.E.
        • de Vos W.M.
        • Barrangou R.
        • Klaenhammer T.R.
        • Caufield P.W.
        • Cui Y.
        • Zhang H.
        • O'Toole P.W.
        Expanding the biotechnology potential of lactobacilli through comparative genomics of 213 strains and associated genera.
        Nat. Commun. 2015; 6 (26415554): 8322
        • Walker B.J.
        • Abeel T.
        • Shea T.
        • Priest M.
        • Abouelliel A.
        • Sakthikumar S.
        • Cuomo C.A.
        • Zeng Q.
        • Wortman J.
        • Young S.K.
        • Earl A.M.
        Pilon: An integrated tool for comprehensive microbial variant detection and genome assembly improvement.
        PLoS One. 2014; 9 (25409509)e112963