Advertisement
Research| Volume 102, ISSUE 5, P3851-3867, May 2019

Metabolomic analysis of significant changes in Lactobacillus casei Zhang during culturing to generation 4,000 under conditions of glucose restriction

  • Lin Pan
    Affiliations
    Key Laboratory of Dairy Biotechnology and Engineering, Education Ministry of China, Huhhot, 010018, Inner Mongolia, China

    Key Laboratory of Dairy Products Processing, Ministry of Agricultural, Department of Food Science and Engineering, Inner Mongolia Agricultural University, Hohhot 010018, China
    Search for articles by this author
  • Jie Yu
    Affiliations
    Key Laboratory of Dairy Biotechnology and Engineering, Education Ministry of China, Huhhot, 010018, Inner Mongolia, China

    Key Laboratory of Dairy Products Processing, Ministry of Agricultural, Department of Food Science and Engineering, Inner Mongolia Agricultural University, Hohhot 010018, China
    Search for articles by this author
  • Dongyan Ren
    Affiliations
    Key Laboratory of Dairy Biotechnology and Engineering, Education Ministry of China, Huhhot, 010018, Inner Mongolia, China

    Key Laboratory of Dairy Products Processing, Ministry of Agricultural, Department of Food Science and Engineering, Inner Mongolia Agricultural University, Hohhot 010018, China
    Search for articles by this author
  • Caiqing Yao
    Affiliations
    Key Laboratory of Dairy Biotechnology and Engineering, Education Ministry of China, Huhhot, 010018, Inner Mongolia, China

    Key Laboratory of Dairy Products Processing, Ministry of Agricultural, Department of Food Science and Engineering, Inner Mongolia Agricultural University, Hohhot 010018, China
    Search for articles by this author
  • Yongfu Chen
    Affiliations
    Key Laboratory of Dairy Biotechnology and Engineering, Education Ministry of China, Huhhot, 010018, Inner Mongolia, China

    Key Laboratory of Dairy Products Processing, Ministry of Agricultural, Department of Food Science and Engineering, Inner Mongolia Agricultural University, Hohhot 010018, China
    Search for articles by this author
  • Bilige Menghe
    Correspondence
    Corresponding author
    Affiliations
    Key Laboratory of Dairy Biotechnology and Engineering, Education Ministry of China, Huhhot, 010018, Inner Mongolia, China

    Key Laboratory of Dairy Products Processing, Ministry of Agricultural, Department of Food Science and Engineering, Inner Mongolia Agricultural University, Hohhot 010018, China
    Search for articles by this author
Open ArchivePublished:March 14, 2019DOI:https://doi.org/10.3168/jds.2018-15702

      ABSTRACT

      Lactic acid bacteria are being consumed more frequently as awareness of their health benefits has increased. The industrial production of lactic acid bacteria requires a comprehensive understanding of their survival stress, especially regarding changes in metabolic substances in a glucose-limited environment. In the present study, a metabolomic approach was applied to investigate Lactobacillus casei Zhang using cultures from a common ancestor that were permitted to evolve under conditions with normal or glucose-restricted media for up to 4,000 generations. Metabolomic analyses of intracellular and extracellular differential metabolites under De Man, Rogosa and Sharpe broth (2% vol/vol glucose; Oxoid Ltd., Basingstoke, UK) and glucose-restricted (0.02% vol/vol glucose in De Man, Rogosa and Sharpe broth) conditions were performed at generations 0, 2,000, and 4,000 and revealed 23 different metabolites. Myristic acid, ergothioneine, Lys-Thr, and palmitamide contents exhibited significant reductions between 0 and 4,000 generations, whereas nicotinate, histidine, palmitic acid, l-lysine, urocanate, thymine, and other substances increased. The dynamics of the pathways involved in AA metabolism, including glycine, serine, and threonine metabolism, histidine metabolism, lysine degradation, and arginine and proline metabolism, were also a focus of the present study. There were also changes in several other metabolic pathways, including vitamin B6, thiamine, nicotinate, and nicotinamide, according to generation time. Additionally, in the present study we screened for key metabolites involved in the glucose-restricted response and provided a theoretical basis for comprehensively revealing the regulatory mechanisms associated with L. casei Zhang glucose restriction at the metabolic level. These findings also provide novel ideas and methods for analyzing the glucose-restricted stress response at the metabolic level.

      Key words

      INTRODUCTION

      Lactic acid bacteria (LAB) can ferment lactose to produce lactic acid and are widely distributed in nature. Lactic acid bacteria in the form of probiotics play a vital role in human health and as an important industrial microorganism. Lactic acid bacteria and their metabolites have been widely used in food, medicine, feed, and so on (
      • Burgess C.
      • O'Connell-Motherway M.
      • Sybesma W.
      • Hugenholtz J.
      • van Sinderen D.
      Riboflavin production in Lactococcus lactis: Potential for in situ production of vitamin-enriched foods.
      ;
      • De Angelis M.
      • Gobbetti M.
      Environmental stress responses in Lactobacillus: A review.
      ;
      • Parvez S.
      • Malik K.A.
      • Ah Kang S.
      • Kim H.Y.
      Probiotics and their fermented food products are beneficial for health.
      ).
      Lactic acid bacteria inevitably face various environmental stresses during their production and application, such as osmotic pressure, acid environments, oxidative stress, cold, heat, pH, and starvation. These stresses can affect the physiological state and characteristics of cells and may limit the growth of LAB in other manufacturing processes (
      • De Angelis M.
      • Gobbetti M.
      Environmental stress responses in Lactobacillus: A review.
      ). Among these stresses in particular, glucose restriction has been reported as one of the most important limitations to growth. During the process of reproduction, if LAB cannot replenish glucose in time, the bacteria will enter a state of starvation and exhibit survival stress; it is worth noting that LAB can also indirectly enter a nutrient-deficient state through other stresses; LAB regulate genetic regulators of different responses associated with other stresses when faced with carbohydrate restriction (
      • Hussain M.A.
      • Knight M.I.
      • Britz M.L.
      Proteomic analysis of lactose-starved Lactobacillus casei during stationary growth phase.
      ).
      • Al-Naseri A.
      • Bowman J.P.
      • Wilson R.
      • Nilsson R.E.
      • Britz M.L.
      Impact of lactose starvation on the physiology of Lactobacillus casei GCRL163 in the presence or absence of Tween 80.
      demonstrated that Lactobacillus casei GCRL163 inhibits phosphotransferase abundance, lactose metabolism, and galactose metabolism under starvation conditions. The same process occurs during synthesis of nucleotides and proteins (
      • Al-Naseri A.
      • Bowman J.P.
      • Wilson R.
      • Nilsson R.E.
      • Britz M.L.
      Impact of lactose starvation on the physiology of Lactobacillus casei GCRL163 in the presence or absence of Tween 80.
      ). Thus, for an experimenter, it is important to understand the conditions that are beneficial or harmful to LAB and the mechanisms that allow them to survive under stress and carry out metabolic activities.
      Lactobacillus casei Zhang is a probiotic bacteria strain isolated from the traditional homemade milk products of herdsmen in Zhenglan Banner, Xilinhot, in the Inner Mongolia region of China. Lactobacillus casei Zhang is highly tolerant to acid and bile salt stress and has probiotic functions that lower blood lipids, improve immune regulation, enhance antioxidation, and antagonize the growth of intestinal pathogens (
      • Ya T.
      • Zhang Q.
      • Chu F.
      • Merritt J.
      • Bilige M.
      • Sun T.
      • Du R.
      • Zhang H.
      Immunological evaluation of Lactobacillus casei Zhang: A newly isolated strain from koumiss in Inner Mongolia, China.
      ;
      • Wu R.
      • Wang L.
      • Wang J.
      • Li H.
      • Menghe B.
      • Wu J.
      • Guo M.
      • Zhang H.
      Isolation and preliminary probiotic selection of lactobacilli from koumiss in Inner Mongolia.
      ;
      • Zhang Y.
      • Du R.
      • Wang L.
      • Zhang H.
      The antioxidative effects of probiotic Lactobacillus casei Zhang on the hyperlipidemic rats.
      ). To assess the effects of L. casei Zhang on upper respiratory tract infections and intestinal diseases, researchers divided 137 volunteers into adult and elderly groups to perform randomized double-blind controlled trials using a daily oral dose of L. casei Zhang (9 log cfu). The administration of L. casei Zhang reduced plasma proinflammatory cytokines (IL-1), increased anti-inflammatory cytokines (IL-4 and IL-10), and upregulated the whole-blood expression levels of cluster of differentiation (CD) 4, CD8, CD44, CD27, and C-X-C motif chemokine receptor 5 (CXCR5). Additionally, the elderly group had an increased number of bowel movements per week, which alleviated red blood cell abnormalities, ameliorated the symptoms of upper respiratory tract infection in adults, and significantly shortened the duration of symptoms. These findings demonstrated the potential of L. casei Zhang for reducing the symptoms of upper respiratory tract infections and maintaining gastrointestinal health (
      • Hor Y.Y.
      • Lew L.C.
      • Lau A.S.Y.
      • Ong J.S.
      • Chuah L.O.
      • Lee Y.Y.
      • Choi S.B.
      • Rashid F.
      • Wahid N.
      • Sun Z.H.
      • Kwok L.Y.
      • Zhang H.P.
      • Liong M.T.
      Probiotic Lactobacillus casei Zhang (LCZ) alleviates respiratory, gastrointestinal and RBC abnormality via immuno-modulatory, anti-inflammatory and anti-oxidative actions.
      ).
      Metabolomics can be used to qualitatively and quantitatively analyze low molecular weight metabolites in cells that are produced by microbial metabolism, including peptides, AA, carbohydrates, nucleic metabolites, vitamins, organic acids, and minerals (
      • Jewett M.C.
      • Hofmann G.
      • Nielsen J.
      Fungal metabolite analysis in genomics and phenomics.
      ;
      • Wishart D.S.
      Metabolomics: Applications to food science and nutrition research.
      ).
      • Foschi C.
      • Laghi L.
      • Parolin C.
      • Giordani B.
      • Compri M.
      • Cevenini R.
      • Marangoni A.
      • Vitali B.
      Novel approaches for the taxonomic and metabolic characterization of lactobacilli: Integration of 16S rRNA gene sequencing with MALDI-TOF MS and 1H-NMR.
      combined 16S rRNA gene sequencing with proton nuclear magnetic resonance to determine the classification and metabolic properties of LAB and found that metabolomics can be used to correlate specific biological activities with taxonomy as well as to understand the mechanisms underlying bacteriostasis in some LAB species. Under glucose-restricted conditions, Lactobacillus plantarum H4 produces acetate in addition to acetic acid in the pH range of 3.8 to 6, and l-lactate accumulation was greater than that of d-lactate (
      • Bobillo M.
      • Marshall V.M.
      Effect of acidic pH and salt on acid end-products by Lactobacillus plantarum in aerated, glucose-limited continuous culture.
      ). These findings demonstrate that L. plantarum H4 adjusts to changes in environmental conditions by altering its metabolic pathways. There is also evidence showing that heterofermentative L. casei can degrade arginine during culture and lead to the formation of ammonia and citrulline (
      • Toh H.
      • Oshima K.
      • Nakano A.
      • Takahata M.
      • Murakami M.
      • Takaki T.
      • Nishiyama H.
      • Igimi S.
      • Hattori M.
      • Morita H.
      Genomic adaptation of the Lactobacillus casei group.
      ). Under stress, Lactobacillus will become deficient in essential nutrients, which limits the exponential phase of microbial growth, allowing it to enter the stationary phase and engage in carbohydrate starvation, resulting in the consumption of cellular energy (
      • De Angelis M.
      • Gobbetti M.
      Environmental stress responses in Lactobacillus: A review.
      ). At present, LAB, as an important fermenting microorganism, are faced with widespread carbon source stress in production applications, which will seriously affect their physiological functions and restrict the development of related products.
      Ultraperformance liquid chromatography quadrupole time-of-flight MS in the MSE mode (UPLC-Q-TOF MSE) was used to analyze changes in metabolites during the different evolutionary processes between the original probiotic L. casei Zhang and other bacteria that were subcultured to generations 2,000 and 4,000 under glucose-restricted conditions. Understanding metabolic changes in L. casei Zhang under nutrient stress will provide practical guidance for the industrial production of probiotics.

      MATERIALS AND METHODS

      Bacterial Strain and Growth Conditions

      The original strain of L. casei Zhang was deposited in the Key Laboratory of Dairy Biotechnology and Engineering, Education Ministry of China. The obtained ampoule strain was activated in De Man, Rogosa and Sharpe (MRS) broth medium (CM0359; Oxoid Ltd., Basingstoke, UK) and then subcultured twice at 37°C in an incubator (MIR-1620; Sanyo, Osaka, Japan).

      Experimental Evolution

      After 2 generations of cultivation (24 h/generation), L. casei Zhang was streaked on MRS agar medium (CM0361; Oxoid Ltd.). After incubation at 37°C for 72 h, 3 colonies were randomly selected and cultured in a normal 2% (vol/vol) glucose MRS broth medium and a 0.02% (vol/vol) glucose-restricted medium. Every 24 h, L. casei Zhang strains cultured in MRS broth or a glucose-restricted environment were inoculated into MRS broth and MRS containing 0.02% glucose sugar at an inoculation amount of 1% (vol/vol, 50 μL); 100-fold bacterial growth per day represents almost 6.6 generations per subculture (
      • Wang J.
      • Dong X.
      • Shao Y.
      • Guo H.
      • Pan L.
      • Hui W.
      • Kwok L.Y.
      • Zhang H.
      • Zhang W.
      Genome adaptive evolution of Lactobacillus casei under long-term antibiotic selection pressures.
      ). The 4,000-generation evolution in the laboratory took approximately 603 d to complete. Strains of different generations were obtained for vacuum freeze-drying preservation during the long-term cultivation (−50°C, 15.0 Pa for 48–72 h).

      Sample Preparation for UPLC-Q-TOF MSE Analysis

      Cultivation of L. casei Zhang

      The original strain of L. casei Zhang, and the L. casei Zhang grown in the glucose-restricted medium deposited in the ampoule tube, were cultured in the MRS. The strains growing in normal or glucose-restricted environments were subcultured in MRS broth twice and the bacterial solution was transferred to fresh medium at 1% volume of inoculation every 24 h.

      Sample Quenching and Metabolite Extraction

      Samples were collected at the mid-exponential phase (∼8 h). The activated samples were dispensed into 2-mL centrifuge tubes and centrifuged at 7,000 × g for 5 min at 4°C; supernatant retention was performed to collect extracellular metabolic samples. Briefly, 200 μL of extracellular sample was obtained, and 800 μL of methanol/acetonitrile at the ratio of 1:1 (vol/vol) was added to each extracellular sample. Then, the sample was vortexed for 30 s and subjected to ultrasound for 10 min (4°C cold bath). Following the ultrasound, samples were allowed to stand at −20°C for 60 min and were then centrifuged at 13,000 × g for 15 min at 4°C. Next, 500 μL of the supernatant was collected, vacuum-dried at 35°C, reconstituted with 200 μL of 50% acetonitrile, vortexed for 30 s, and subjected to ultrasound for 10 min (4°C cold bath). Next, the sample was centrifuged at 13,000 × g for 15 min at 4°C and kept on the machine for testing.
      The bacteria sludge cells were quickly added to a precooled 0.9% NaCl, vortexed, centrifuged at 0°C for 5 min at 5,000 × g, and washed. This process was repeated twice, and the supernatant was discarded; then, the remaining bacterial sludge was collected and −40°C precooled 60% methanol/0.9% NaCl was quickly added to the centrifuge tube for quenching using 4 times the volume of the sludge (
      • de Koning W.
      • Vandam K.
      A method for the determination of changes of glycolytic metabolites in yeast on a subsecond time scale using extraction at neutral pH.
      ;
      • Jensen N.B.
      • Jokumsen K.V.
      • Villadsen J.
      Determination of the phosphorylated sugars of the Embden-Meyerhoff-Parnas pathway in Lactococcus lactis using a fast sampling technique and solid phase extraction.
      ;
      • Oldiges M.
      • Takors R.
      Applying metabolic profiling techniques for stimulus-response experiments: Chances and pitfalls.
      ;
      • Mashego M.R.
      • Rumbold K.
      • De Mey M.
      • Vandamme E.
      • Soetaert W.
      • Heijnen J.J.
      Microbial metabolomics: Past, present and future methodologies.
      ;
      • Canelas A.B.
      • Ras C.
      • Pierick A.T.
      • van Dam J.C.
      • Heijnen J.J.
      • van Gulik W.M.
      Leakage-free rapid quenching technique for yeast metabolomics.
      ). The cells were then quenched for 5 min and centrifuged at 10,000 × g for 5 min at 0°C; the supernatant was then discarded. The quenched cells were subjected to liquid nitrogen milling for 20 min (a mortar was precooled at −80°C) and ground to a fine powder; 50 mg of powder was placed in a 1.5-mL centrifuge tube. Then, 1 mL of −40°C 70% methanol/water mixture was added to the tube, vortexed, fully mixed, and extracted to suspend the cells. Liquid nitrogen was used to repeatedly freeze and thaw the cells (4 times; liquid nitrogen freezing for 3 min, −20°C refrigerator melting), followed by centrifugation at 10,000 × g for 5 min at −20°C; the supernatant was then collected. The cell pellet was continuously suspended in 1 mL of −40°C 70% methanol/water, vortexed for 30 s, and centrifuged at 8,000 × g for 5 min at 0°C; the supernatants obtained by the 2 extractions were then combined (
      • Chen M.M.
      • Li A.L.
      • Sun M.C.
      • Feng Z.
      • Meng X.C.
      • Wang Y.
      Optimization of the quenching method for metabolomics analysis of Lactobacillus bulgaricus..
      ). The supernatant was stored at −80°C until analysis; this fraction was labeled as the intracellular sample. Each sample (2 replicates per sample) was processed in parallel to yield 6 replicates.

      UPLC-Q-TOF MSE Analysis

      An Acquity UPLC/Xevo G2 Q TOF-MSE instrument (Waters Corp., Milford, MA) was used to extract the MS spectral peaks of the samples. Samples were filtered and placed in a sample vial for testing. When the samples in this experiment were analyzed by metabolomics, the quadrupole time-of-flight mass spectrometer carried out data scanning using positive electrospray ionization (ESI+) and negative electrospray ionization (ESI) ion modes. The cone voltage was 40 V, the ion energy was 1, the collision energy was 6, the sample infusion flow rate was 10 µL/min, the sample fill volume was 250 µL, the LockSpray capillary was at 3 kV, and the acquisition scan range was 100 to 1,000 mass-to-charge ratio (m/z) at a scan frequency of 10.00 s.
      To ensure the accuracy of data acquisition, leucine-cephalin was used as the lock mass for both positive ion mode (556.2771 [M + H]+) and negative ion mode (554.2615 [M − H]). Pretreated samples were chromatographed using a C18-UPLC system. The column was an HSS T3 (2.1 mm × 100 mm, 1.8 µm; Waters Corp.) with a sample infusion flow rate of 0.3 µL/min; the column was maintained at 35°C. The injection volume was 4 µL.
      For ESI+ ion mode, mobile phase A was an aqueous phase containing 0.1% formic acid (Sigma-Aldrich, St. Louis, MO), and mobile phase B was an organic phase containing 0.1% formic acid in acetonitrile (Optima LC/MS-grade; Fisher Scientific, Grand Island, NY). The elution gradient is shown in Table 1. For scanning in ESI ion mode, the mobile phase consisted of water (Millipore, Bedford, MA) containing 10% acetonitrile (solution A) and 50% acetonitrile (solution B). The elution gradient was as shown in Table 2.
      Table 1Ultraperformance liquid chromatography quadrupole time-of-flight MS in the MSE mode mobile phase gradient program for the analysis of samples in positive ion mode
      Time (min)Flow rate (mL/min)Mobile phase A (%)Mobile phase B (%)Curve
      Initial0.3955Initial
      20.39556
      70.380206
      120.335656
      140.320806
      170.320806
      180.39556
      Table 2Ultraperformance liquid chromatography quadrupole time-of-flight MS in the MSE mode mobile phase gradient program for the analysis of samples in negative ion mode
      Time (min)Flow rate (mL/min)Mobile phase A (%)Mobile phase B (%)Curve
      Initial0.39010Initial
      20.390106
      70.380206
      120.330706
      140.330706
      150.390106

      UPLC-Q-TOF MSE Data Processing and Multivariate Statistical Analyses

      The data processing method has been described in detail previously (
      • Pan L.
      • Yu J.
      • Mi Z.
      • Mo L.
      • Jin H.
      • Yao C.
      • Ren D.
      • Menghe B.
      A metabolomics approach uncovers differences between traditional and commercial dairy products in Buryatia (Russian Federation).
      ). Briefly, the data were preprocessed in MassLynx version 4.1 (Waters Corp.) and then imported into the EZinfo software (Waters Corp.) for centralization and Pareto homogenization to reduce deviation due to large-scale peaks.
      Simultaneous multivariate statistical analysis was performed using SIMCA-P + 11 (Umetrics AB, Umeå, Sweden) and MetaboAnalyst (www.metaboanalyst.ca) software. For each treatment, 3 comparisons were made: 0 versus 2,000 generations, 0 versus 4,000 generations, and 2,000 versus 4,000 generations. To determine changes in metabolites between the different generations of L. casei Zhang, variables with variable importance in the projection values greater than 1 (fold change ≥2 and P ≤ 0.5) were selected as potential metabolic discriminant variables (
      • Chong I.G.
      • Jun C.H.
      Performance of some variable selection methods when multicollinearity is present.
      ;
      • Xi X.
      • Kwok L.Y.
      • Wang Y.
      • Ma C.
      • Mi Z.
      • Zhang H.
      Ultra-performance liquid chromatography-quadrupole-time of flight mass spectrometry MSE-based untargeted milk metabolomics in dairy cows with subclinical or clinical mastitis.
      ). Based on these variables, the different metabolic substances were further analyzed and identified.

      Identification of Differential Metabolites

      The significantly different substances obtained by screening were automatically searched for in several databases, including METLIN (http://www.metlin.scripps.edu/), KEGG (www.genome.jp/kegg), ChemSpider (www.chemspider.com), and the Human Metabolome Database (www.hmdb.ca) based on their mass-to-charge ratio and retention time. Finally, correlation and metabolic pathway analyses were performed on the various retrieved metabolites.

      RESULTS

      Sample Principal Component Analysis Results

      To compare the differences between L. casei Zhang generation 0 and glucose-restricted generation 2,000, multivariate analyses were used to segment the effective spectral area, and principal component analysis (PCA) was performed. In the PCA model, the abscissa is the first principal component and the ordinate is the second principal component. In each score graph, the addition of the 2 principal components explained more than 95% of the variance, as shown in Figure 1.
      Figure thumbnail gr1
      Figure 1Principal component (PC) analysis score plots of intracellular (IC) samples. Data were obtained in positive electrospray ionization (ESI+: A, B, C) and negative electrospray ionization (ESI: D, E, F) ion modes. A and D: generation 0 versus generation 2,000; B and E: generation 0 versus generation 4,000; C and F: generation 2,000 versus generation 4,000.
      In Figure 1, it can be seen that all PCA score plots have clear separation between different generations, indicating that the chemical compositions of samples differed among L. casei Zhang under stress conditions. Taking Figure 1A as an example, the 2 sets of samples are distributed in different regions, without intersection or overlap, and there is an obvious discrete phenomenon. Six parallel samples of each type are closely aggregated.
      Samples between generation 0 and glucose-restricted generation 2,000 had values of 78.90 and 19.90% in the principal component 1 and principal component 2 directions, respectively, and the 2 principal components accounted for most of the variance among the different treatment samples. Principal component 1 had the largest contribution rate; this indicates that the metabolites of generations 0 and 2,000 are clearly dissociable and confirms the effectiveness of the sampling method and the reproducibility of the test. The extracellular analysis method is the same as the intracellular method; the PCA results are shown in Figure 2.
      Figure thumbnail gr2
      Figure 2Principal component (PC) analysis score plots of extracellular (EC) samples. Data were obtained in positive electrospray ionization (ESI+; A, B, C) and negative electrospray ionization (ESI; D, E, F) ion modes. A and D: generation 0 versus generation 2,000; B and E: generation 0 versus generation 4,000; C and F: generation 2,000 versus generation 4,000.
      Figures 2A and D and Figures 2B and E are well separated by PCA analysis; only Figures 2C and F overlap, which indicates that the PCA model cannot adequately distinguish between the 2 groups. In addition, the overall dispersion of the 2 types of sample is large, indicating that the individual differences among the samples are also relatively large. However, as can be seen from Figure 1, Figure 2, the PCA plots of both intracellular and extracellular metabolites are well separated, indicating that the metabolites of different generations will vary significantly during the laboratory evolution of L. casei Zhang in a glucose-restricted environment.

      Sample Orthogonal Partial Least Squares Discriminant Analysis Results

      The differences between the extracellular and intracellular metabolite samples were further examined by orthogonal partial least squares discriminant analysis (OPLS-DA). In OPLS-DA, The R2Y parameter describes the percentage of variation explained by the model and Q2 describes the predictive ability of the model. In Figure 3A, obtained through OPLS-DA in positive ion mode, the R2Y and Q2 values of the generation 0 versus generation 2,000 intracellular metabolite samples are equal to 0.996 and 0.936, respectively. This indicates that the 2 models are reliable and that the metabolic samples had significant differences. The R2Y and Q2 values of all models were relatively high, indicating that the models were reliable and that both sets of samples in each OPLS-DA score plot had metabolic differences.
      Figure thumbnail gr3
      Figure 3Orthogonal partial least squares discriminant analysis score plots and S-plots of the data for intracellular (IC) samples obtained in positive electrospray ionization (ESI+) ion mode. A and D: generation 0 versus 2,000; B and E: generation 0 versus 4,000; C and F: generation 2,000 versus 4,000.
      In Figure 3D, each point represents a variable, and the greater the distance of a variable from the origin point, the higher the confidence level regarding its contribution to the clustering observed in the OPLS-DA score plots. The variables in each S-plot map were screened for final identification of the substance based on variable importance in the projection and fold-change values. The OPLS-DA score plots and S-plots from generation 0 to 4,000 in both positive and negative ion modes are shown in Supplemental Figures S1 through S3 (https://doi.org/10.3168/jds.2018-15702).

      Identification of Different Metabolites

      According to the data on retention time and mass-to-charge ratio, a total of 8,376 differentially expressed metabolites were found between generations 0 and 2,000, including intracellular and extracellular metabolites; among these, 5,298 were in ESI+ ion mode. A total of 10,234 substances differed between generations 0 and 4,000, including 7,374 in ESI+ ion mode, whereas 3,952 differentially expressed substances were found between generations 2,000 and 4,000 (of which 2,205 were in ESI+ ion mode). The differential substances ultimately conforming to the screening principle described above are discussed below, and the differential metabolites of intracellular generations 0 and 2,000 are summarized in Table 3.
      Table 3Generation 0 versus generation 2,000: differential intracellular metabolites at positive and negative ions
      RT_Exact MS
      RT_Exact MS = retention time and exact mass.
      MetabolitesMolecular formulaVIP
      Variable importance in the projection.
      Fold changeP-value
      Increased in positive ion mode
       1.13_119.0473l-ThreonineC4H9NO31.041720.136428.49E-05
       2.13_120.0790l-ErythroseC4H8O43.859270.111672.49E-09
       1.12_123.0424NicotinateC6H4NO21.360290.0902435.56E-08
       4.89_126.0533ThymineC5H6N2O21.647260.277251.05E-05
       8.52_135.0793AdenineC5H5N51.894820.285610.010734
       4.66_136.0599HypoxanthineC5H4N4O1.479240.0942835.87E-06
       3.32_138.0532UrocanateC6H6N6O21.366210.0395060.00082435
       4.89_144.0645Octanoic acidC8H16O21.089490.0728310.0016261
       4.20_146.0588l-LysC6H14N2O21.615730.0260494.72E-07
       4.90_150.0542d-RibuloseC5H10O51.271490.185163.58E-06
       1.13_165.0556l-PheC9H11NO21.035630.0664469.08E-06
       2.13_166.0857d-Arabinonic acidC5H10O61.148510.081385.05E-08
       4.58_168.0648HomogentisateC8H7O41.272140.277230.0067467
       6.77_169.0758PyridoxineC8H11NO31.224690.225460.0003405
       4.88_186.0762Ala–ProC8H14N2O32.294690.292181.05E-07
       4.19_188.0713Homo-cis-aconitateC7H8O61.712210.0464982.83E-08
       1.99_189.0875N-Acetyl-glutamateC7H11NO52.306440.28412.35E-06
       4.07_204.0869TryptophanC11H12N2O22.669640.188144.20E-07
       4.07_210.0767d-Glucaric acidC6H10O81.701030.139641.78E-07
       11.45_228.2328Myristic acidC14H28O21.039950.259290.00015464
       8.62_243.0886CytidineC9H13N3O52.343930.254641.18E-06
       9.31_255.0661DihydroneopterinC9H13N5O42.622420.227831.92E-05
       13.32_256.2646Palmitic acidC16H32O21.260380.431520.010902
       9.49_285.0758Leu–Gly–ProC13H23N3O41.871680.158031.83E-06
       4.89_294.0978Ile–TyrC15H22N2O41.222010.307065.19E-05
       13.39_300.2901Retinoic acidC20H28O21.140880.154030.00024178
       9.64_309.0878Tyr–GlnC14H19N3O51.551370.205513.47E-05
       4.89_372.1299BiocytinC16H28N4O4S1.088220.0615361.61E-05
      Reduced in positive ion mode
       7.65_131.0685CreatineC4H9N3O21.059280.00911710.0091171
       8.25_133.0845Aspartic acidC4H7NO42.669560.00200910.0020091
       10.07_139.1100HistidinalC6H9N3O1.483950.00183360.0018336
       7.65_175.0966CitrullineC6H13N3O31.105330.00367250.0036725
       7.65_177.1121N-Formyl-l-methionineC6H11NO3S2.672140.00429310.0042931
       7.32_221.1387N-Acetyl-d-glucosamineC8H15NO61.524760.0025520.002552
       1.06_230.1872ErgothioneineC9H15N3O2S1.676430.000505110.00050511
       7.64_283.1750Gln–HisC11H17N5O41.248460.000434460.00043446
       7.32_309.1909Tyr–LysC15H23N3O41.401440.00862280.0086228
      Increased in negative ion mode
       1.05_100.98411-Aminocyclopropane 1-carboxylateC4H7NO21.591470.00218282.93E-08
       1.47_112.9843CreatinineC4H7N3O2.384910.0904310.00021513
       2.16_122.9663NicotinateC6H5NO21.513720.10686.97E-05
       1.84_144.97372-OxoglutaramateC5H7NO41.429250.169860.004854
       1.04_147.0443O-Acetyl-l-serineC5H9NO41.584740.00870793.48E-05
       1.20_150.9617GuanineC5H5N51.037730.00516945.16E-09
       1.99_162.9586N-Acetyl-l-cysteineC5H9NO3S1.88440.138530.0001293
       9.67_243.1951Lys–ProC11H21N3O31.773310.00141160.0018486
       12.62_246.9924Pyridoxal-5-phosphateC8H10NO6P1.024090.122210.0057913
       11.54_283.2632StearamideC18H37NO4.00320.000296920.0004265
       1.81_116.9579l-Aspartate-4-semialdehydeC4H7NO31.070030.00460037.43E-06
      Reduced in negative ion mode
       9.09_199.1696LauramideC12H25NO1.35532536.060.0060526
       9.46_243.1958Gln–ProC10H17N3O41.76496525.263.21E-05
       10.01_265.1458l-Aspartic acid l-ornithineC9H19N3O62.306611,107.92.84E-05
       10.67_293.1779Phe–GlnC14H19N3O42.14739503.540.0021298
       10.37_309.1712Tyr–GlnC14H19N3O51.2100339.550.00059251
      1 RT_Exact MS = retention time and exact mass.
      2 Variable importance in the projection.
      In Table 3 we can see that, via metabolomic analysis, 53 significantly different substances between generations 0 and 2,000 were identified based on the screening principle, containing AA, organic acids, vitamins, pyrimidines, and short peptides. In total, 39 were upregulated and are involved in histidine metabolism; tryptophan metabolism; biotin metabolism; pyruvate metabolism; alanine, aspartate, and glutamate metabolism; pyrimidine metabolism; purine metabolism; valine, leucine, and isoleucine biosynthesis; and vitamin digestion and absorption. There were 14 downregulated substances involved in arginine and proline metabolism; glycine, serine, and threonine metabolism; histidine metabolism; arginine biosynthesis; and biosynthesis of AA.
      In the extracellular metabolomic analysis of generations 0 and 2,000, 81 significantly different metabolites were identified under ESI+ and ESI modes (Table 4). These 81 substances are involved in metabolic pathways including carbohydrate, energy, lipid, nucleotide, and AA metabolism. The number of significantly different metabolites was larger than that of intracellular substances, and more metabolic pathways were involved.
      Table 4Generation 0 versus generation 2,000: differential extracellular metabolites at positive and negative ions
      RT_Exact MS
      RT_Exact MS = retention time and exact mass.
      MetabolitesMolecular formulaVIP
      Variable importance in the projection.
      Fold changeP-value
      Increased in positive ion mode
       13.04_106.0836GlycerateC3H6O41.376960.0647540.00010451
       9.04_111.0422HistamineC5H9N31.323710.0601766.25E-06
       9.05_113.0576CreatinineC4H7N3O1.12740.244789.76E-07
       11.62_121.0276l-CysteineC3H7NO2S1.334410.036470.00068041
       3.07_126.0533ThymineC5H6N2O21.380250.377882.22E-06
       3.54_129.1007PipecolateC6H10NO21.345340.109384.83E-07
       9.01_131.0691l-IsoleucineC6H13NO21.363240.34641.11E-06
       6.93_133.0848Aspartic acidC4H7NO42.627250.0924831.01E-05
       4.64_136.0608HypoxanthineC5H4N4O1.121880.409622.02E-07
       8.99_137.0588TyramineC8H11NO1.44780.258622.27E-08
       5.45_138.054UrocanateC6H6N2O21.1920.109810.0033048
       11.62_149.0229l-MethionineC5H11NO2S4.240850.228810.0011476
       4.86_150.0543d-RibuloseC5H10O51.773830.208332.52E-09
       4.51_150.0548GuanineC5H5N5O51.157840.0997190.0020137
       8.84_155.0701HistidineC6H9N3O21.563220.364951.90E-05
       5.45_168.0654PhosphoenolpyruvateC3H7O6P1.377410.0847760.0017744
       6.71_169.0765PyridoxineC8H11NO31.409960.318772.53E-06
       9.00_175.0967CitrullineC6H13N3O31.476530.363464.00E-06
       7.95_177.1123N-Formyl-l-methionineC6H11NO3S3.511790.0326431.56E-07
       10.47_181.0767TyrosineC9H11NO31.102250.120716.00E-07
       1.78_189.0879N-Acetyl-glutamateC7H11NO52.386580.000530120.015579
       4.16_195.1137DopaquinoneC9H9NO41.834360.240378.53E-08
       9.03_199.0972O-Phospho-l-homoserineC4H10NO6P1.599520.0986391.18E-05
       9.82_200.2018Lauric acidC12H24O21.317530.181461.20E-05
       5.41_204.0876TryptophanC11H12N2O21.851810.198210.0045673
       5.56_210.0773PhosphocreatineC4H10N3O5P1.353990.187120.0090698
       4.48_218.1504Thr–ValC9H18N2O41.373570.119983.61E-06
       3.85_219.1164Asn–SerC7H13N3O51.25030.247153.86E-07
       1.08_223.0049N-Acetyl-l-tyrosineC11H13NO41.145260.186512.46E-05
       1.13_226.0459CarnosineC9H14N4O31.079020.240061.05E-07
       11.39_228.234Myristic acidC14H28O21.726470.11293.59E-07
       11.78_240.2334PalmitaldehydeC16H32O1.129810.315352.46E-06
       9.03_243.1236CytidineC9H13N3O51.918260.133772.27E-06
       2.44_252.0747UbiquinolC19H28O41.024640.200190.00018987
       13.15_256.265Palmitic acidC16H32O26.071640.0948824.06E-11
       4.46_258.1459Cys–HisC9H14N4O3S1.11130.0896580.018617
       9.03_261.1327Gly–TrpC13H15N3O31.1110.180450.0031832
       9.61_263.0826Asn–MetC9H17N3O4S1.444260.269986.02E-09
       10.47_274.2172Lys–GlnC11H22N4O41.449440.0917743.14E-06
       4.59_295.0962Asn–TyrC13H17N3O51.314430.00251369.77E-08
       13.22_300.29069-Cis-retinoic acidC20H28O22.597460.0743042.66E-08
       4.24_303.1382Val–TrpC16H21N3O31.160460.450853.01E-09
       9.04_305.1589Trp–ThrC15H19N3O41.076580.272430.0062877
       9.61_309.0875Tyr–GlnC14H19N3O51.18214.76764.58E-05
       4.85_332.1941Gln–TrpC16H20N4O41.731750.0971696.31E-08
      Reduced in positive ion mode
       4.16_146.0595GlutamineC5H10N2O31.607152.62544.15E-08
       1.87_166.0473d-Arabinonic acidC5H10O61.129643.49926.09E-05
       5.06_199.1808LauramideC12H25NO1.1325732.1681.43E-09
       4.15_188.0715Homo-cis-aconitateC7H8O62.09792990.717.71E-12
       1.39_212.1041Pro–ProC10H16N2O31.087733.84160.001635
       1.07_230.1873ErgothioneineC9H15N3O2S2.882496.990.00089026
       4.38_235.1089Asn–CysC7H13N3O4S1.3173.33165.77E-05
       2.80_247.1302Lys–ThrC10H21N3O41.20217.7827.83E-09
       1.60_254.1625His–ValC11H18N4O31.073876.27527.48E-07
       10.47_255.0778N-Ribosyl-nicotinamideC11H16N2O51.331822.61371.13E-09
       1.24_276.1448SaccharopineC11H20N2O61.3969911.1011.49E-06
       1.97_310.1292PorphyrinC20H14N41.1275215.6959.83E-07
       3.87_343.15His–TrpC17H19N5O31.123323.50891.95E-05
       8.97_360.2057AldosteroneC21H28O51.23831349.48.39E-09
      Increased in negative ion mode
       1.14_100.0299IsohexanalC6H12O1.925320.123980.014571
       2.02_100.97431-Aminocyclopropane 1-carboxylateC4H7NO22.11696−1.82446.06E-05
       1.12_101.01315-AminopentanalC5H11NO2.11987−1.71470.012176
       1.20_112.97371-Pyrroline 2-carboxylateC5H7NO21.50273−1.17260.011679
       1.70_126.9932Piperideine-2-carboxylateC6H9NO21.11773−8.36040.0061875
       2.21_144.96252-OxoglutaramateC5H7NO41.88478−1.94460.014777
      Reduced in negative ion mode
       1.15_116.04162-OxoisovalerateC5H8O31.193535.91190.0010377
       1.12_130.02281-Pyrroline 4-hydroxy-2- carboxylateC5H8O3N1.41212.9160.0017742
       8.01_144.036Octanoic acidC8H16O21.0472234.8866.55E-05
       9.14_199.161LauramideC12H25NO1.62047630.160.0014492
       8.85_201.1129HeteropyrithiamineC11H13N41.13337279.329.77E-05
       1.13_203.0665Lys–AlaC9H19N3O31.2869511.9540.0031665
       6.03_220.0714Asp–SerC7H12N2O61.231983.67680.016655
       10.57_233.0657Lys–SerC9H19N3O41.12395290.40.00056947
       4.94_243.1724CytidineC9H13N3O51.33214419.50.0010608
       8.62_243.1911Lys–ProC11H21N3O32.24671990.021.09E-09
       1.12_251.0834DeoxyadenosineC10H13N5O32.0697839.6290.00071983
       5.16_255.0814DihydroneopterinC9H13O4N51.128514.25260.01321
       10.60_255.2246PalmitamideC16H33NO7.816271,963.71.01E-05
       11.31_283.2548StearamideC18H37NO7.99357131223.67E-06
       9.35_299.1956DehydrosphinganinC18H37NO21.1018665.7540.0013811
       4.40_314.1949ProgesteroneC21H30O21.5214591.050.0041312
      1 RT_Exact MS = retention time and exact mass.
      2 Variable importance in the projection.
      In addition, L. casei Zhang was subcultured to generation 4,000. The differential metabolites of generations 0 and 4,000 are shown in Supplemental Tables S1 and S2 (https://doi.org/10.3168/jds.2018-15702). Forty-seven different substances, including the following were found in generation 4,000 but not in generation 2,000: the AA l-ornithine, l-leucine, and l-homoserine; the organic acids salicylic acid, succinic, acid, and linolenic acid; the short peptides valine–alanine, isoleucine–proline, leucine–serine, histidine–cysteine, and tartrate; and pyridoxamine, 2-phospho-glycerate, N-acetylcytidine, guanosine, and indole. Meanwhile, 49 metabolites in generation 2,000 were not detected in generation 4,000 and hence were not screened.
      Finally, we conducted a comparative analysis of the metabolites that differed between L. casei Zhang generations 2,000 and 4,000. The results are shown in Supplemental Tables S3 and S4 (https://doi.org/10.3168/jds.2018-15702). There were 66 metabolomic differences between generations 2,000 and 4,000, including both intracellular and extracellular metabolites, of which 37 were upregulated and 29 were downregulated. Figure 4 shows that 23 metabolites were shared by all 3 groups; generation 0 versus generation 2,000, generation 0 versus generation 4,000 and generation 2,000 versus generation 4,000.
      Figure thumbnail gr4
      Figure 4The distribution of significant metabolic differences among the 3 generations (generation 0 vs. 2,000; generation 0 vs. 4,000; generation 2,000 vs. 4,000).
      The content changes of the 23 different substances in the 3 groups from generation 0 to 4,000 are shown in Figure 5. In generation 0, the contents of myristic acid, ergothioneine, lysine–threonine, and palmitamide were higher than the contents of other metabolites. As the experiment progressed and they were subcultured to generation 4,000, the levels of the substances with a high content in generation 0 decreased, which may have been due to adaptability during long-term generation. It was also obvious that the color of substances such as phosphocreatine, hypoxantine, octanoic acid, l-lysine, O-acetyl-l-serine, urocanate, l-phenylalanine, stearamide, lysine–serine, thymine, and d-ribulose changed from blue to green, which indicates that they were slightly more abundant in generation 4,000 than generation 0.
      Figure thumbnail gr5
      Figure 5Changes in the contents of 23 substances differentially expressed between generations 0 and 4,000.

      Metabolite Enrichment Pathway Analysis

      All 175 significantly different metabolites (including intracellular and extracellular metabolites from generation 0 vs. 2,000, generation 0 vs. 4,000, and generation 2,000 vs. 4,000) were analyzed by a pathway analysis (Figure 6). As shown in Figure 6, metabolites from the 0 versus 2,000, 0 versus 4,000, and 2,000 versus 4,000 generation comparisons are involved in various metabolic processes, such as nucleotide metabolism and AA metabolism.
      Figure thumbnail gr6
      Figure 6Pathway analysis of significantly differentiated metabolites. (A) Generation 0 versus 2,000. 1 = alanine, aspartate, and glutamate metabolism; 2 = histidine metabolism; 3 = glycine, serine, and threonine metabolism; 4 = vitamin B6 metabolism; 5 = cysteine and methionine metabolism; 6 = lysine biosynthesis; 7 = aminoacyl-tRNA biosynthesis; and 8 = arginine and proline metabolism. (B) Generation 0 versus 4,000. 1 = histidine metabolism; 2 = alanine, aspartate, and glutamate metabolism; 3 = glycine, serine, and threonine metabolism; 4 = arginine and proline metabolism; 5 = thiamine metabolism; 6 = phenylalanine metabolism; and 7 = aminoacyl-tRNA biosynthesis. (C) Generation 2,000 versus 4,000. 1 = histidine metabolism; 2 = nicotinate and nicotinamide metabolism; 3 = lysine degradation; 4 = arginine and proline metabolism; 5 = cysteine and methionine metabolism; 6 = alanine, aspartate, and glutamate metabolism; and 7 = aminoacyl-tRNA biosynthesis.
      Supplemental Figure S4 (https://doi.org/10.3168/jds.2018-15702) shows a metabolite enrichment analysis chart. In total, 32 metabolic subpathways were shared by all 3 generations, including histidine metabolism, fatty acid metabolism, arginine and proline metabolism, the urea cycle, β-alanine metabolism, glycine and serine metabolism, aspartate metabolism, purine metabolism, nicotinate and nicotinamide metabolism, glycerolipid metabolism, lysine degradation, tyrosine metabolism, and vitamin B6 metabolism.
      For the generation 0 versus 2,000 comparison, 3 metabolic pathways for metabolite enrichment, histidine metabolism, arginine and proline metabolism, and vitamin B6 metabolism were observed, as shown in Figure 6. Similarly, arginine and proline metabolism as well as thiamine metabolism were observed in the 0 versus 4,000 generation comparison. For the 2,000 versus 4,000 generation comparison, lysine degradation, nicotinate and nicotinamide metabolism, and histidine metabolism were observed in both the enrichment and pathway analyses.
      All of the differential substances identified in the present experiment were analyzed in terms of their metabolic pathways (Table 5). The metabolic pathway analysis revealed that the metabolites that significantly changed during subculturing of the original strain from generation 0 to 4,000 were mainly responsible for the changes in alanine, aspartate, and glutamate metabolism; glycine, serine, and threonine metabolism; histidine metabolism; cysteine and methionine metabolism; valine, leucine, and isoleucine biosynthesis; α-linolenic acid metabolism; and other metabolic pathways (Table 5).
      Table 5Significant metabolic pathways of metabolites from generation 0 to 2,000 and 4,000
      Pathway nameRaw P-value
      The original P-value calculated from the enrichment analysis.
      −log (P-value)Holm P-valueFDR
      False discovery rate.
      Impact
      The pathway impact value calculated from pathway topology analysis.
      Alanine, aspartate, and glutamate metabolism0.00231716.06740.171460.0231710.54606
      Glycine, serine, and threonine metabolism1.43E-0511.1580.00112655.70E-040.34468
      Histidine metabolism0.0016616.40030.124570.0221460.3032
      Cysteine and methionine metabolism0.00159566.44050.121270.0221460.20644
      Valine, leucine, and isoleucine biosynthesis0.0039945.5230.283570.0319520.20403
      α-Linolenic acid metabolism0.702590.35298110.20335
      Arginine and proline metabolism0.00342955.67530.246930.0304850.19684
      Tryptophan metabolism0.216831.528710.691820.17943
      Retinol metabolism0.224841.492410.691820.17514
      Lysine degradation0.299061.207110.854460.17297
      Thiamine metabolism0.00231716.06740.171460.0231710.16019
      Nicotinate and nicotinamide metabolism0.265241.327110.785890.14392
      Tyrosine metabolism0.194821.635710.691820.14174
      Lysine biosynthesis0.00847944.77010.593560.0616680.13361
      Phenylalanine metabolism0.0339433.383110.210720.13206
      Folate biosynthesis0.516290.66109110.13123
      Biotin metabolism0.0710422.644510.378890.12195
      Pyrimidine metabolism0.0945342.358810.44490.11351
      Aminoacyl-tRNA biosynthesis9.77E-0713.8397.82E-057.82E-050.11268
      Steroid hormone biosynthesis0.37720.9749810.914420.11043
      Sulfur metabolism0.0342423.374310.210720.10985
      Vitamin B6 metabolism0.377040.9754110.914420.10376
      Glycolysis or gluconeogenesis0.72660.31937110.1035
      Ascorbate and aldarate metabolism0.848650.16411110.07571
      Pantothenate and CoA biosynthesis0.0945412.358710.44490.07366
      Ubiquinone and other terpenoid-quinone biosynthesis0.435020.83237110.0709
      Sphingolipid metabolism0.648130.43366110.06438
      Arachidonic acid metabolism0.926540.076294110.05672
      Fatty acid metabolism0.612130.4908110.04482
      Valine, leucine, and isoleucine degradation0.221061.509310.691820.03889
      Purine metabolism0.168941.778210.675750.03469
      Glyoxylate and dicarboxylate metabolism0.333131.099210.888360.03291
      d-Glutamine and d-glutamate metabolism0.367581.000810.914420.02674
      β-Alanine metabolism0.102892.274110.457290.02328
      Pentose phosphate pathway0.737880.30397110.02181
      Glycerolipid metabolism0.737880.30397110.0206
      Butanoate metabolism0.812990.20704110.01774
      Amino sugar and nucleotide sugar metabolism0.882540.12495110.01487
      Citrate cycle (tricarboxylic acid cycle)0.194531.637110.691820.01446
      Nitrogen metabolism7.88E-047.1460.060680.0157610.0083
      Phenylalanine, tyrosine, and tryptophan biosynthesis5.77E-047.45720.0450270.0153940.008
      Pentose and glucuronate interconversions0.892220.11404110.00638
      Selenoamino acid metabolism0.600910.50931110.00482
      Glutathione metabolism0.199591.611510.691820.0019
      Propanoate metabolism0.769020.26264110.00134
      d-Arginine and d-ornithine metabolism0.0391443.240510.223680
      Cyanoamino acid metabolism0.13611.994310.573070
      Fatty acid biosynthesis0.321771.133910.887630
      Taurine and hypotaurine metabolism0.5660.56917110
      Fatty acid elongation in mitochondria0.676490.39083110
      Pyruvate metabolism0.737880.30397110
      Porphyrin and chlorophyll metabolism0.987970.012104110
      1 The original P-value calculated from the enrichment analysis.
      2 False discovery rate.
      3 The pathway impact value calculated from pathway topology analysis.
      The metabolites of generation 0 versus those of generation 4,000 were analyzed using iPath (https://pathways.embl.de/ipath3.cgi; Figure 7). The primary metabolic pathways of the laboratory evolution of L. casei Zhang affected by the differential metabolites were AA metabolism, carbohydrate metabolism, nucleotide metabolism and cofactors, and vitamin metabolism.
      Figure thumbnail gr7
      Figure 7Analysis of significant metabolites using iPath (https://pathways.embl.de/ipath3.cgi) metabolic pathways between generation 0 and 4,000. Red node = metabolites; red line = interaction.

      DISCUSSION

      Effects of Glucose Restriction on Carbohydrate Metabolism

      A glucose-restricted environment generally occurs during the stationary phase of microbial growth. Under normal culture conditions, metabolic diversity of probiotic L. casei Zhang will occur. A significant scarcity of carbon sources can cause lactobacillus to stop growing (to ensure survival), such that it enters a stable period. Because glucose is depleted during later growth stages, the cells begin to use other sugars (
      • Wang J.
      • Zhang W.
      • Zhong Z.
      • Wei A.
      • Bao Q.
      • Zhang Y.
      • Sun T.
      • Postnikoff A.
      • Meng H.
      • Zhang H.
      Gene expression profile of probiotic Lactobacillus casei Zhang during the late stage of milk fermentation.
      ).
      • Borch E.
      • Berg H.
      • Holst O.
      Heterolactic fermentation by a homofermentative Lactobacillus sp. during glucose limitation in anaerobic continuous culture with complete cell recycle.
      studied Lactobacillus sp. under limited glucose conditions in a continuous anaerobic culture with complete cell recycling and found that homofermentation transitioned into heterolactic fermentation when glucose was continuously depleted within the cell cycle such that the cells were starved. Under conditions of limited glucose, especially during cell cycle reculture, AA are widely produced and utilized, and some sulfides are produced (
      • Borch E.
      • Berg H.
      • Holst O.
      Heterolactic fermentation by a homofermentative Lactobacillus sp. during glucose limitation in anaerobic continuous culture with complete cell recycle.
      ). This is consistent with the present experimental results.

      Effects of Glucose Restriction on AA Metabolism

      Amino acids are involved in many physiological functions of cells. Under stress conditions, AA can accumulate in significant amounts and function as cell-protective agents (
      • Devantier R.
      • Scheithauer B.
      • Villas-Bôas S.G.
      • Pedersen S.
      • Olsson L.
      Metabolite profiling for analysis of yeast stress response during very high gravity ethanol fermentations.
      ;
      • Jozefczuk S.
      • Klie S.
      • Catchpole G.
      • Szymanski J.
      • Cuadros-Inostroza A.
      • Steinhauser D.
      • Selbig J.
      • Willmitzer L.
      Metabolomic and transcriptomic stress response of Escherichia coli.
      ;
      • Kang H.J.
      • Yang H.J.
      • Kim M.J.
      • Han E.S.
      • Kim H.J.
      • Kwon D.Y.
      Metabolomic analysis of meju during fermentation by ultra performance liquid chromatography-quadrupole-time of flight mass spectrometry (UPLC-Q-TOF MS).
      ;
      • Li H.
      • Ma M.
      • Luo S.
      • Zhang R.
      • Han P.
      • Hu W.
      Metabolic responses to ethanol in Saccharomyces cerevisiae using a gas chromatography tandem mass spectrometry-based metabolomics approach.
      ). Of the 23 substances that were detected at the intersection of the 3 experimental groups, the contents of substances involved in the AA metabolic pathway were elevated. For example, urocanic acid was detected in generation 0, 2,000, and 4,000, and the relative content increased with generation time. Urocanic acid is believed to be useful as a sun-screening agent in cosmetics and pharmaceuticals (
      • Yamamoto K.
      • Sato T.
      • Tosa T.
      • Chibata I.
      Continuous production of urocanic acid by immobilized Achromobacter liquidum cells.
      ), and
      • Hogan J.A.
      • Maniatis A.
      • Moloney W.C.
      The serum Lactobacillus casei folate clearance test in various hematologic disorders.
      reported that elevated levels of urocanic acid may be due to reduced folate levels. Urocanic acid is the final product of histidine metabolism, which breaks histidine down into glutamic acid in the presence of formiminoglutamic acid (
      • Chanarin I.
      Urocanic acid and formimino-glutamic acid excretion in megaloblastic anaemia and other conditions: The effect of specific therapy.
      ). Humans metabolize histidine to uric acid via the activity of l-histidine ammonium lyase, which then enters the tricarboxylic acid cycle (
      • Bentley S.
      • Bottarelli A.
      • Bonardi S.
      Histamine as biogenic amines in food.
      ). Several AA such as glycine, serine, threonine, and alanine in Saccharomyces cerevisiae cells increased significantly (
      • Li H.
      • Ma M.
      • Luo S.
      • Zhang R.
      • Han P.
      • Hu W.
      Metabolic responses to ethanol in Saccharomyces cerevisiae using a gas chromatography tandem mass spectrometry-based metabolomics approach.
      ). Similarly, changes in AA caused by starvation can also lead to changes in proteins. According to
      • Kunji E.R.S.
      • Ubbink T.
      • Matin A.
      • Poolman B.
      • Konings W.N.
      Physiological responses to Lactococcus lactis ML3 to altering conditions of growth and starvation.
      , arginine starvation in Lactococcus lactis ssp. lactis led to the induction of 15 proteins, whereas 14 proteins were induced after galactose starvation. Bacteria adapt to nutrition limitations by assuming a physiological state that is characterized by downregulation of nucleic acid and protein synthesis and simultaneous upregulation of protein degradation and AA synthesis (
      • Chatterji D.
      • Ojha A.K.
      Revisiting the stringent response, ppGpp and starvation signaling.
      ;
      • Hussain M.A.
      • Knight M.I.
      • Britz M.L.
      Proteomic analysis of lactose-starved Lactobacillus casei during stationary growth phase.
      ). Some researchers also reported that the size and fatty acid composition of microbial cells will change under growth conditions adapted to nutrient deficiencies (
      • Rice S.A.
      • Oliver J.D.
      Starvation response of the marine barophile CNPT-3.
      ). At the same time, the rate of cellular protein synthesis decreases, and the protein expression associated with long-term survival or environmental stress increases significantly (
      • Schultz J.E.
      • Matin A.
      Molecular and functional characterization of a carbon starvation gene of Escherichia coli..
      ). In the present study, 23 unique metabolites were identified between generations 2,000 and 4,000. For example, O-phospho-l-serine was involved in glycine, serine, and threonine metabolism, glyceraldehyde-3-phosphate was involved in glycolysis and gluconeogenesis and vitamin B6 metabolism, and 2-oxo-5-aminovalerate was involved in arginine and proline metabolism. The levels of these 3 substances were elevated when subcultured in a glucose-limited environment, which corresponds to the present finding that most of the significant metabolites participated in the AA metabolic pathway. In a study investigating the restricted carbon source culturing of L. lactis IL-1403,
      • Redon E.
      • Loubiere P.
      • Cocaign-Bousquet M.
      Transcriptome analysis of the progressive adaptation of Lactococcus lactis to carbon starvation.
      reported that citric acid is used preferentially when a carbon is scarce. To maintain the supply of cellular energy, the cells hydrolyze arginine to produce ornithine through the arginine deiminase pathway during the growth phase and then produce ATP to supply energy (
      • Redon E.
      • Loubiere P.
      • Cocaign-Bousquet M.
      Transcriptome analysis of the progressive adaptation of Lactococcus lactis to carbon starvation.
      . The accumulation of these AA represents a stress response of the cells to carbon stress to preserve the integrity of the cells.

      Effects of Glucose Restriction on Nucleotide Metabolism

      Thymine was the base of deoxyribonucleic acid in the present experiment, and the thymine content increased during the subcultuing of L. casei Zhang from generation 0 to 4,000.
      • Soska J.
      Growth of Lactobacillus acidophilus in the absence of folic acid.
      cultured Lactobacillus acidophilus R-26 in a medium without folic acid and found that the organism required thymine in addition to deoxynucleosides, purines, pyrimidines, and most AA. When thymine was present in the medium, the bacteria grew exponentially and without limitation. The presence of thymine can meet the demand for folic acid by Lactobacillus leichmannii (
      • Thompson R.B.
      Observations on the effects of vitamin B12, liver extracts, folic acid and thymine on the maturation of megaloblasts in culture.
      ) and, in the absence of thymine, there is a 150% increase in DNA content that induces an extra cycle of chromosome replication (
      • Elli M.
      • Zink R.
      • Reniero R.
      • Morelli L.
      Growth requirements of Lactobacillus johnsonii in skim and UHT milk.
      ;
      • Soska J.
      Growth of Lactobacillus acidophilus in the absence of folic acid.
      ). Uridine is a nucleoside that serves as a major component of ribonucleic acid. In the present study, uridine increased from generation 2,000 to 4,000, whereas guanosine decreased from generation 0 to 4,000. Guanosine is a proven taste enhancer when fermenting bread with Lactobacillus plantarum (
      • McCabe C.
      • Rolls E.T.
      Umami: A delicious flavor formed by convergence of taste and olfactory pathways in the human brain.
      ) and can provide significant growth stimulation for the Lactobacillus johnsonii NCC533 strain; however, nucleotide supplementation is required for L. plantarum growth in reconstituted skim milk (
      • Elli M.
      • Zink R.
      • Reniero R.
      • Morelli L.
      Growth requirements of Lactobacillus johnsonii in skim and UHT milk.
      ). Furthermore, Lactobacillus gasseri PA-3 shows greater proliferation in the presence of guanosine during culture (
      • Yamada N.
      • Saito C.
      • Murayama-Chiba Y.
      • Kano H.
      • Asami Y.
      Lactobacillus gasseri PA-3 utilizes the purines GMP and guanosine and decreases their absorption in rats.
      ). Metabolomic analyses were used to assess Streptococcus intermedius under aerobic and anaerobic conditions, and it was revealed that the intracellular guanosine concentration was higher during anaerobic growth and that the cells absorbed more under anaerobic conditions (
      • Fan F.
      • Michelle L.M.
      • Brian E.M.
      • Dawn M.E.B.
      • Michael G.M.
      Metabolic and transcriptomic profiling of Streptococcus intermedius during aerobic and anaerobic growth.
      ). Additionally, analyses of the presence of several typical molecules associated with the flavor and taste of foods during the fermentation of bread with Lactobacillus plantarum have detected uridine and guanosine (
      • Valerio F.
      • Conte A.
      • Di Biase M.
      • Lattanzio V.M.T.
      • Lonigro S.L.
      • Padalino L.
      • Pontonio E.
      • Lavermicocca P.
      Formulation of yeast-leavened bread with reduced salt content by using a Lactobacillus plantarum fermentation product.
      ).

      Effects of Glucose Restriction on the Metabolism of Cofactors and Vitamins

      9-Cis-retinoic acid, which is associated with retinol metabolism, is the active substance in vitamin A (
      • Priyamvada S.
      • Anbazhagan A.N.
      • Gujral T.
      • Borthakur A.
      • Saksena S.
      • Gill R.K.
      • Alrefai W.A.
      • Dudeja P.K.
      All-trans-retinoic acid increases SLC 26A3 DRA (down-regulated in adenoma) expression in intestinal epithelial cells via HNF-1β..
      ) and vitamin B6 and is associated with the metabolism of pyridoxine and pyridoxamine during the metabolism of cofactors and vitamins. It has been reported that L. casei can be cultured on a medium containing riboflavin, pantothenic acid, nicotinic acid, and pyridine as AA sources. It has also been reported that L. casei can be cultured on a medium containing riboflavin, pantothenic acid, nicotinic acid, and pyridine as AA sources (
      • Hutchings B.
      • Bohonos N.
      • Peterson W.H.
      Growth factors for bacteria. 13. Purification and properties of an eluate factor required by certain lactic acid bacteria.
      ).
      In the present study, there was a significant decrease in nicotinate, which is associated with nicotinate and nicotinamide metabolism, from generation 2,000 to 4,000. The essential growth factors needed by Lactobacillus include riboflavin, pantothenic acid, nicotinic acid, pyridoxine, biotin, and the eluate factor such as folic acid (
      • Hutchings B.
      • Bohonos N.
      • Peterson W.H.
      Growth factors for bacteria. 13. Purification and properties of an eluate factor required by certain lactic acid bacteria.
      ). Lactobacillus gasseri can produce NAD using nicotinate and nicotinamide in the presence of nicotinamidase and nicotinate phosphoribosyltransferase (
      • Foster J.W.
      • Moat A.G.
      Nicotinamide adenine dinucleotide biosynthesis and pyridine nucleotide cycle metabolism in microbial systems.
      ), whereas L. johnsonii and L. gasseri lack the corresponding enzymes to synthesize thiamine, biotin, pyridoxine, and other substances involved in vitamin digestion and absorption (
      • Azcarate-Peril M.A.
      • Altermann E.
      • Goh Y.J.
      • Tallon R.
      • Sanozky-Dawes R.B.
      • Pfeiler E.A.
      • O'Flaherty S.
      • Buck B.L.
      • Dobson A.
      • Duong T.
      • Miller M.J.
      • Barrangou R.
      • Klaenhammer T.R.
      Analysis of the genome sequence of Lactobacillus gasseri ATCC 33323 reveals the molecular basis of an autochthonous intestinal organism.
      ). For L. casei LC2W, nicotinate is an essential nutrient under anaerobic and aerobic conditions and can increase the yield of L. casei LC2W exopolysaccharide (
      • Xu N.
      • Liu J.
      • Ai L.Z.
      • Liu L.M.
      Reconstruction and analysis of the genome-scale metabolic model of Lactobacillus casei LC2W.
      ).
      Although nicotinate content was reduced by the influence of glucose restriction during the process of generation in the present study, L. casei Zhang still produced beneficial substances. Therefore, when L. casei Zhang is used as a starter, it may lead to a decrease in beneficial substances when long-term glucose sources are insufficient.
      Due to the current practice of application of L. casei Zhang in food and feedstuffs, it is important to study the metabolic kinetics of L. casei Zhang under glucose-restricted conditions. Metabolomic analyses of the metabolites of L. casei Zhang during long-term passage under conditions of glucose restriction will provide indispensable supporting data for further understanding the properties of L. casei Zhang.

      CONCLUSIONS

      We studied metabolomics-related differences between an original strain of L. casei Zhang and strains cultured in a glucose-restricted medium for 2,000 and 4,000 generations using ultraperformance liquid chromatography quadrupole time-of-flight MS. Significant changes in metabolites were observed under the glucose-restricted conditions. Changes in the metabolites, including AA, nucleotides, and short peptides, resulted in changes in metabolic pathways; the changes were more pronounced when passaged to generation 4,000. The metabolically significant metabolites significantly affected several metabolic pathways, including histidine metabolism, arginine and proline metabolism, vitamin B6 metabolism, thiamine metabolism, lysine degradation, and nicotinate and nicotinamide metabolism. Metabolomics was used to study changes in the metabolites associated with L. casei Zhang under conditions of glucose restriction and to screen for key metabolites involved in the glucose-restricted response. The present findings may provide a theoretical basis for a comprehensive understanding of the metabolic regulatory mechanisms of L. casei Zhang during glucose restriction. These findings may also lead to novel ideas and methods for assessing the glucose-restricted stress response at a metabolic level.

      ACKNOWLEDGMENTS

      This research was supported by National Natural Science Foundation of China (Beijing, China; Grant No. 31660450), Program Funded by University for Fostering Outstanding Young Scholars (Grant No. 2017XYQ-4), the Natural Science Foundation of Inner Mongolia, China (Grant No. 2018MS03041), and the Natural Science Foundation of Inner Mongolia, China (Grant No. 2017ZD07).

      Supplementary Material

      REFERENCES

        • Al-Naseri A.
        • Bowman J.P.
        • Wilson R.
        • Nilsson R.E.
        • Britz M.L.
        Impact of lactose starvation on the physiology of Lactobacillus casei GCRL163 in the presence or absence of Tween 80.
        J. Proteome Res. 2013; 12 (24066708): 5313-5322
        • Azcarate-Peril M.A.
        • Altermann E.
        • Goh Y.J.
        • Tallon R.
        • Sanozky-Dawes R.B.
        • Pfeiler E.A.
        • O'Flaherty S.
        • Buck B.L.
        • Dobson A.
        • Duong T.
        • Miller M.J.
        • Barrangou R.
        • Klaenhammer T.R.
        Analysis of the genome sequence of Lactobacillus gasseri ATCC 33323 reveals the molecular basis of an autochthonous intestinal organism.
        Appl. Environ. Microbiol. 2008; 74 (18539810): 4610-4625
        • Bentley S.
        • Bottarelli A.
        • Bonardi S.
        Histamine as biogenic amines in food.
        Ingeniería Alimentos. 1995; 5: 35
        • Bobillo M.
        • Marshall V.M.
        Effect of acidic pH and salt on acid end-products by Lactobacillus plantarum in aerated, glucose-limited continuous culture.
        J. Appl. Bacteriol. 1992; 73: 67-70
        • Borch E.
        • Berg H.
        • Holst O.
        Heterolactic fermentation by a homofermentative Lactobacillus sp. during glucose limitation in anaerobic continuous culture with complete cell recycle.
        J. Appl. Microbiol. 1991; 71: 265-269
        • Burgess C.
        • O'Connell-Motherway M.
        • Sybesma W.
        • Hugenholtz J.
        • van Sinderen D.
        Riboflavin production in Lactococcus lactis: Potential for in situ production of vitamin-enriched foods.
        Appl. Environ. Microbiol. 2004; 70 (15466513): 5769-5777
        • Canelas A.B.
        • Ras C.
        • Pierick A.T.
        • van Dam J.C.
        • Heijnen J.J.
        • van Gulik W.M.
        Leakage-free rapid quenching technique for yeast metabolomics.
        Metabolomics. 2008; 4: 226-239
        • Chanarin I.
        Urocanic acid and formimino-glutamic acid excretion in megaloblastic anaemia and other conditions: The effect of specific therapy.
        Br. J. Haematol. 1963; 9 (14019995): 141-157
        • Chatterji D.
        • Ojha A.K.
        Revisiting the stringent response, ppGpp and starvation signaling.
        Curr. Opin. Microbiol. 2001; 4 (11282471): 160-165
        • Chen M.M.
        • Li A.L.
        • Sun M.C.
        • Feng Z.
        • Meng X.C.
        • Wang Y.
        Optimization of the quenching method for metabolomics analysis of Lactobacillus bulgaricus..
        J. Zhejiang Univ. Sci. B. 2014; 15 (24711354): 333-342
        • Chong I.G.
        • Jun C.H.
        Performance of some variable selection methods when multicollinearity is present.
        Chemom. Intell. Lab. Syst. 2005; 78: 103-112
        • De Angelis M.
        • Gobbetti M.
        Environmental stress responses in Lactobacillus: A review.
        Proteomics. 2004; 4 (14730676): 106-122
        • de Koning W.
        • Vandam K.
        A method for the determination of changes of glycolytic metabolites in yeast on a subsecond time scale using extraction at neutral pH.
        Anal. Biochem. 1992; 204 (1514678): 118-123
        • Devantier R.
        • Scheithauer B.
        • Villas-Bôas S.G.
        • Pedersen S.
        • Olsson L.
        Metabolite profiling for analysis of yeast stress response during very high gravity ethanol fermentations.
        Biotechnol. Bioeng. 2005; 90 (15812801): 703-714
        • Elli M.
        • Zink R.
        • Reniero R.
        • Morelli L.
        Growth requirements of Lactobacillus johnsonii in skim and UHT milk.
        Int. Dairy J. 1999; 9: 507-513
        • Fan F.
        • Michelle L.M.
        • Brian E.M.
        • Dawn M.E.B.
        • Michael G.M.
        Metabolic and transcriptomic profiling of Streptococcus intermedius during aerobic and anaerobic growth.
        Metabolomics. 2016; 12: 46
        • Foschi C.
        • Laghi L.
        • Parolin C.
        • Giordani B.
        • Compri M.
        • Cevenini R.
        • Marangoni A.
        • Vitali B.
        Novel approaches for the taxonomic and metabolic characterization of lactobacilli: Integration of 16S rRNA gene sequencing with MALDI-TOF MS and 1H-NMR.
        PLoS One. 2017; 12 (28207855): e0172483
        • Foster J.W.
        • Moat A.G.
        Nicotinamide adenine dinucleotide biosynthesis and pyridine nucleotide cycle metabolism in microbial systems.
        Microbiol. Rev. 1980; 44 (6997723): 83-105
        • Hogan J.A.
        • Maniatis A.
        • Moloney W.C.
        The serum Lactobacillus casei folate clearance test in various hematologic disorders.
        Blood. 1964; 24 (14200254): 187-197
        • Hor Y.Y.
        • Lew L.C.
        • Lau A.S.Y.
        • Ong J.S.
        • Chuah L.O.
        • Lee Y.Y.
        • Choi S.B.
        • Rashid F.
        • Wahid N.
        • Sun Z.H.
        • Kwok L.Y.
        • Zhang H.P.
        • Liong M.T.
        Probiotic Lactobacillus casei Zhang (LCZ) alleviates respiratory, gastrointestinal and RBC abnormality via immuno-modulatory, anti-inflammatory and anti-oxidative actions.
        J. Funct. Foods. 2018; 44: 235-245
        • Hussain M.A.
        • Knight M.I.
        • Britz M.L.
        Proteomic analysis of lactose-starved Lactobacillus casei during stationary growth phase.
        J. Appl. Microbiol. 2009; 106 (19302099): 764-773
        • Hutchings B.
        • Bohonos N.
        • Peterson W.H.
        Growth factors for bacteria. 13. Purification and properties of an eluate factor required by certain lactic acid bacteria.
        J. Biol. Chem. 1941; 141: 521-528
        • Jensen N.B.
        • Jokumsen K.V.
        • Villadsen J.
        Determination of the phosphorylated sugars of the Embden-Meyerhoff-Parnas pathway in Lactococcus lactis using a fast sampling technique and solid phase extraction.
        Biotechnol. Bioeng. 1999; 63 (10099615): 356-362
        • Jewett M.C.
        • Hofmann G.
        • Nielsen J.
        Fungal metabolite analysis in genomics and phenomics.
        Curr. Opin. Biotechnol. 2006; 17 (16488600): 191-197
        • Jozefczuk S.
        • Klie S.
        • Catchpole G.
        • Szymanski J.
        • Cuadros-Inostroza A.
        • Steinhauser D.
        • Selbig J.
        • Willmitzer L.
        Metabolomic and transcriptomic stress response of Escherichia coli.
        Mol. Syst. Biol. 2010; 6 (20461071): 364
        • Kang H.J.
        • Yang H.J.
        • Kim M.J.
        • Han E.S.
        • Kim H.J.
        • Kwon D.Y.
        Metabolomic analysis of meju during fermentation by ultra performance liquid chromatography-quadrupole-time of flight mass spectrometry (UPLC-Q-TOF MS).
        Food Chem. 2011; 127 (25214096): 1056-1064
        • Kunji E.R.S.
        • Ubbink T.
        • Matin A.
        • Poolman B.
        • Konings W.N.
        Physiological responses to Lactococcus lactis ML3 to altering conditions of growth and starvation.
        Arch. Microbiol. 1993; 159: 372-379
        • Li H.
        • Ma M.
        • Luo S.
        • Zhang R.
        • Han P.
        • Hu W.
        Metabolic responses to ethanol in Saccharomyces cerevisiae using a gas chromatography tandem mass spectrometry-based metabolomics approach.
        Int. J. Biochem. Cell Biol. 2012; 44 (22504284): 1087-1096
        • Mashego M.R.
        • Rumbold K.
        • De Mey M.
        • Vandamme E.
        • Soetaert W.
        • Heijnen J.J.
        Microbial metabolomics: Past, present and future methodologies.
        Biotechnol. Lett. 2007; 29 (17091378): 1-16
        • McCabe C.
        • Rolls E.T.
        Umami: A delicious flavor formed by convergence of taste and olfactory pathways in the human brain.
        Eur. J. Neurosci. 2007; 25 (17432971): 1855-1864
        • Oldiges M.
        • Takors R.
        Applying metabolic profiling techniques for stimulus-response experiments: Chances and pitfalls.
        Adv. Biochem. Eng. Biotechnol. 2005; 92 (15791937): 173-196
        • Pan L.
        • Yu J.
        • Mi Z.
        • Mo L.
        • Jin H.
        • Yao C.
        • Ren D.
        • Menghe B.
        A metabolomics approach uncovers differences between traditional and commercial dairy products in Buryatia (Russian Federation).
        Molecules. 2018; 23 (29565828): E735
        • Parvez S.
        • Malik K.A.
        • Ah Kang S.
        • Kim H.Y.
        Probiotics and their fermented food products are beneficial for health.
        J. Appl. Microbiol. 2006; 100 (16696665): 1171-1185
        • Priyamvada S.
        • Anbazhagan A.N.
        • Gujral T.
        • Borthakur A.
        • Saksena S.
        • Gill R.K.
        • Alrefai W.A.
        • Dudeja P.K.
        All-trans-retinoic acid increases SLC 26A3 DRA (down-regulated in adenoma) expression in intestinal epithelial cells via HNF-1β..
        J. Biol. Chem. 2015; 290 (25887398): 15066-15077
        • Redon E.
        • Loubiere P.
        • Cocaign-Bousquet M.
        Transcriptome analysis of the progressive adaptation of Lactococcus lactis to carbon starvation.
        J. Bacteriol. 2005; 187 (15866950): 3589-3592
        • Rice S.A.
        • Oliver J.D.
        Starvation response of the marine barophile CNPT-3.
        Appl. Environ. Microbiol. 1992; 58 (16348747): 2432-2437
        • Schultz J.E.
        • Matin A.
        Molecular and functional characterization of a carbon starvation gene of Escherichia coli..
        J. Mol. Biol. 1991; 218: 129-140
        • Soska J.
        Growth of Lactobacillus acidophilus in the absence of folic acid.
        J. Bacteriol. 1966; 91 (4957024): 1840-1847
        • Thompson R.B.
        Observations on the effects of vitamin B12, liver extracts, folic acid and thymine on the maturation of megaloblasts in culture.
        Blood. 1952; 7 (14925145): 522-525
        • Toh H.
        • Oshima K.
        • Nakano A.
        • Takahata M.
        • Murakami M.
        • Takaki T.
        • Nishiyama H.
        • Igimi S.
        • Hattori M.
        • Morita H.
        Genomic adaptation of the Lactobacillus casei group.
        PLoS One. 2013; 8 (24116025): e75073
        • Valerio F.
        • Conte A.
        • Di Biase M.
        • Lattanzio V.M.T.
        • Lonigro S.L.
        • Padalino L.
        • Pontonio E.
        • Lavermicocca P.
        Formulation of yeast-leavened bread with reduced salt content by using a Lactobacillus plantarum fermentation product.
        Food Chem. 2017; 221 (27979245): 582-589
        • Wang J.
        • Dong X.
        • Shao Y.
        • Guo H.
        • Pan L.
        • Hui W.
        • Kwok L.Y.
        • Zhang H.
        • Zhang W.
        Genome adaptive evolution of Lactobacillus casei under long-term antibiotic selection pressures.
        BMC Genomics. 2017; 18 (28438179): 320
        • Wang J.
        • Zhang W.
        • Zhong Z.
        • Wei A.
        • Bao Q.
        • Zhang Y.
        • Sun T.
        • Postnikoff A.
        • Meng H.
        • Zhang H.
        Gene expression profile of probiotic Lactobacillus casei Zhang during the late stage of milk fermentation.
        Food Control. 2012; 25: 321-327
        • Wishart D.S.
        Metabolomics: Applications to food science and nutrition research.
        Trends Food Sci. Technol. 2008; 19: 482-493
        • Wu R.
        • Wang L.
        • Wang J.
        • Li H.
        • Menghe B.
        • Wu J.
        • Guo M.
        • Zhang H.
        Isolation and preliminary probiotic selection of lactobacilli from koumiss in Inner Mongolia.
        J. Basic Microbiol. 2009; 49 (19219898): 318-326
        • Xi X.
        • Kwok L.Y.
        • Wang Y.
        • Ma C.
        • Mi Z.
        • Zhang H.
        Ultra-performance liquid chromatography-quadrupole-time of flight mass spectrometry MSE-based untargeted milk metabolomics in dairy cows with subclinical or clinical mastitis.
        J. Dairy Sci. 2017; 100 (28342601): 4884-4896
        • Xu N.
        • Liu J.
        • Ai L.Z.
        • Liu L.M.
        Reconstruction and analysis of the genome-scale metabolic model of Lactobacillus casei LC2W.
        Gene. 2015; 554 (25452194): 140-147
        • Ya T.
        • Zhang Q.
        • Chu F.
        • Merritt J.
        • Bilige M.
        • Sun T.
        • Du R.
        • Zhang H.
        Immunological evaluation of Lactobacillus casei Zhang: A newly isolated strain from koumiss in Inner Mongolia, China.
        BMC Immunol. 2008; 9 (19019236): 68
        • Yamada N.
        • Saito C.
        • Murayama-Chiba Y.
        • Kano H.
        • Asami Y.
        Lactobacillus gasseri PA-3 utilizes the purines GMP and guanosine and decreases their absorption in rats.
        Nucleosides Nucleotides Nucleic Acids. 2018; 3 (29723107): 1-9
        • Yamamoto K.
        • Sato T.
        • Tosa T.
        • Chibata I.
        Continuous production of urocanic acid by immobilized Achromobacter liquidum cells.
        Biotechnol. Bioeng. 1974; 16 (4441634): 1601-1610
        • Zhang Y.
        • Du R.
        • Wang L.
        • Zhang H.
        The antioxidative effects of probiotic Lactobacillus casei Zhang on the hyperlipidemic rats.
        Eur. Food Res. Technol. 2010; 231: 151-158