Next Article in Journal
Possibility of Using Different Calcium Compounds for the Manufacture of Fresh Acid Rennet Cheese from Goat’s Milk
Previous Article in Journal
Re-Rolling Treatment in the Fermentation Process Improves the Aroma Quality of Black Tea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Untargeted Metabolomics and Physicochemical Analysis Revealed the Quality Formation Mechanism in Fermented Milk Inoculated with Lactobacillus brevis and Kluyveromyces marxianus Isolated from Traditional Fermented Milk

College of Food Science, Sichuan Agricultural University, Ya’an 625014, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Foods 2023, 12(19), 3704; https://doi.org/10.3390/foods12193704
Submission received: 30 August 2023 / Revised: 24 September 2023 / Accepted: 27 September 2023 / Published: 9 October 2023
(This article belongs to the Section Food Microbiology)

Abstract

:
Traditional fermented milk from the western Sichuan plateau of China has a unique flavor and rich microbial diversity. This study explored the quality formation mechanism in fermented milk inoculated with Lactobacillus brevis NZ4 and Kluyveromyces marxianus SY11 (MFM), the dominant microorganisms isolated from traditional dairy products in western nan. The results indicated that MFM displayed better overall quality than the milk fermented with L. brevis NZ4 (LFM) and K. marxianus SY11 (KFM), respectively. MFM exhibited good sensory quality, more organic acid types, more free amino acids and esters, and moderate acidity and ethanol concentrations. Non-targeted metabolomics showed a total of 885 metabolites annotated in the samples, representing 204 differential metabolites between MFM and LFM and 163 between MFM and KFM. MFM displayed higher levels of N-acetyl-L-glutamic acid, cysteinyl serine, glaucarubin, and other substances. The differential metabolites were mainly enriched in pathways such as glycerophospholipid metabolism, arginine biosynthesis, and beta-alanine metabolism. This study speculated that L. brevis affected K. marxianus growth via its metabolites, while the mixed fermentation of these strains significantly changed the metabolism pathway of flavor-related substances, especially glycerophospholipid metabolism. Furthermore, mixed fermentation modified the flavor and quality of fermented milk by affecting cell growth and metabolic pathways.

1. Introduction

Western Sichuan presents unique microbial resources due to its specific geographical and climatic conditions. Naturally fermented yak yogurt contains abundant lactic acid bacteria (LAB) and yeasts, contributing a unique flavor and texture to fermented dairy products on the western Sichuan plateau of China. Screening and utilizing the dominant strains can help develop fermented dairy products with regional characteristics and provide nutritional enrichment.
Fermented dairy products are popular with consumers due to their rich nutritional value and unique flavor. Commercial fermented milk is mainly produced using LAB as a starter, while traditional products like milk dregs and kefir primarily use yeasts in addition to LAB, such as Kluyveromyces marxianus and Saccharomyces cerevisiae [1,2,3]. Ni [4] studied the yeast diversity in traditionally fermented products in Xinjiang and Qinghai. The results showed that the yeast species and quantity varied in different regions, with K. unisporus (49.43%) and K. marxianus (26.44%) dominating in Xinjiang, and P. fermentans (41.89%), S. cerevisiae (27.02%), and K. marxianus (24.32%) being the most abundant in Qinghai. K. marxianus is widely present in traditional fermented dairy products. Screening and applying these yeasts as starters or co-starters can help develop new dairy products with local characteristics.
In milk, lactose represents the main carbohydrate. When LAB and yeast are simultaneously inoculated into fermented milk, they can utilize lactose metabolites in complementary ways. Therefore, mixed fermentation is widely used in dairy products [5,6]. The interaction between LAB and yeast during the fermentation process is extremely complex and includes the complementation of their metabolites, influence on cell growth, and impact on the population sensing phenomenon [7,8,9,10]. The complementarity between LAB and yeast in metabolic products is the most common. Arakawa [11] found that LAB uses lactose to produce galactose, which is utilized by yeast to produce pyruvate, while LAB and yeast use pyruvate to form lactic acid.
Sun [12] used metabolomics to compare the metabolite changes during single and mixed fermentation with L. casei M2-L2 and S. cerevisiae M4-Y2, revealing that mixed fermentation produced more long-chain fatty acids. Higher methionine, phenylalanine, tyrosine, and arginine levels were found in the mixed fermentation system, confirming that amino acids might promote LAB proliferation.
Therefore, this study uses quality analysis and metabolomics technology to examine small- and medium-sized molecular metabolite production during LAB and yeast fermentation while analyzing the metabolic pathways involved in the fermentation process using the bioinformatics database [13] and assessing the mechanism behind flavor substance formation. This may provide a basis for accurately regulating the flavor of LAB- and yeast-fermented milk.

2. Materials and Methods

2.1. Materials and Equipment

L. brevis NZ4 (OQ660312) and K. marxianus SY11 (MZ144218) were screened from naturally fermented dairy products in western Sichuan, China. Milk was purchased from local pasture, while the UHPLC-Q Exactive HF-X was obtained from Waltham, MA, USA. The pure methanol, acetonitrile, formic acid, and propanol were supplied by Fisher from Waltham, MA, USA.

2.2. Sample Preparation

2.2.1. Preparation of the Strains

The activated L. brevis NZ4 and K. marxianus SY11 strains were inoculated into the milk and fermented at 36 °C and 28 °C, respectively, until set; after which, they were refrigerated at 4 °C as starters.

2.2.2. Preparation of the Fermented Milk

The sterilized milk was inoculated with the starters at a 5% ratio and fermented at 35 °C until set. The fermented milks inoculated with L. brevis NZ4, K. marxianus SY11, and L. brevis NZ4 + K. marxianus SY11 (1:1) were labeled LFM, KFM, and MFM, respectively. The samples were refrigerated at 4 °C for subsequent examination.

2.3. Sensory Evaluation, Physical and Chemical Indicators, and Viable Counts

Ten panelists were selected and trained to evaluate the tissue status, odor, and taste of the fermented milk. The acidity and ethanol content were determined according to the GB 5009.239-2016 standard of China and a method described by Wang [14], using phenolphthalein as the indicator, titration with 0.1 N NaOH, and calculation. The organic acids content was obtained via HPLC according to the GB 5009.239-2016 standard of China, while amino acids were determined using an amino acid analyzer, reference GB 5009.14-2016. The live LAB was counted according to the GB 4789.35-2016 standard of China, using the dilution coating plate method, and the same method applied to the counting of yeasts, but referenced the GB 4789.15-2016 standard of China.

2.4. UHPLC-Q Exactive HF-X Data Acquisition

2.4.1. Sample Pretreatment and Application

Here, 200 mg of the fermented milk and grinding beads was placed in a centrifuge tube; after which, 400 μL of an extraction solution (methanol: water = 4:1 (v:v)) containing 0.02 mg/mL of internal standard (L-2-chlorophenylalanin) was added for metabolite extraction. The extract was frozen and ground for 6 min (−10 °C, 50 Hz), followed by low-temperature ultrasound extraction for 30 min (5 °C, 40 kHz). The sample was allowed to stand at −20 °C for 30 min and centrifuged for 15 min (4 °C, 13,000 r/min); after which, the supernatant was collected for analysis. The sample metabolites were mixed at equal volumes to prepare a quality control sample (QC).

2.4.2. Liquid Chromatography Conditions

A UHPLC-Q Exactive HF-X system with an HSS T3 chromatographic column (100 mm × 2.1 mm i.d., 1.8 µm) was used to separate a 2 μL sample, which was subjected to mass spectrometric detection. Mobile phase A consisted of 95% water +5% acetonitrile (containing 0.1% formic acid), while mobile phase B comprised 47.5% acetonitrile +47.5% isopropanol +5% water (containing 0.1% formic acid). These procedures were performed using methods described by Meng [15].

2.4.3. Mass Spectrometric Conditions

The mass spectrometric data were collected using a Thermo UHPLC-Q Exactive HF-X Mass Spectrometer equipped with an electrospray ionization (ESI) source operating in either positive or negative ion mode. The optimal conditions included a heating temperature of 425 °C, a capillary temperature of 325 °C, a sheath gas flow rate of 50 arb, an aux gas flow rate of 13 arb, an ion-spray voltage floating (ISVF) of −3500 V in negative mode and 3500 V in positive mode, and a normalized collision energy of 20–40–60 V rolling for MS/MS. The full MS resolution was 60,000, while that of MS/MS was 7500. The data were acquired in Data-Dependent Acquisition (DDA) mode, while detection occurred in a 70–1050 m/z mass range. These procedures were performed using methods described by Fang [16].

2.4.4. Data Processing and Metabolite Annotation

The raw LC–MS data were imported into Progenesis QI metabolomics processing software (Version 2.0, Waters Corporation, Milford, CT, USA) for baseline filtering, peak recognition, integration, retention time correction, and peak alignment. The mass spectrometric information was matched with the HMDB and Metlin databases to obtain the metabolite data. After searching the database, the data matrix was filled with missing values after removing them using an 80% rule. The response intensity of the sample mass spectrum peak was subjected to sum normalization to obtain the data matrix. The variables with a relative standard deviation (RSD) exceeding 30% in the QC samples were deleted, and the data matrix was obtained via log10 logarithmization for subsequent analysis.
The preprocessed matrix file was subjected to differential assessment using R software package ropls (Version 1.6.2) for principal component analysis (PCA) and orthogonal partial least squares discriminant analysis (OPLS-DA). Differential metabolite selection was based on the variable importance in projection (VIP) value obtained via OPLS-DA and the p-value of the Student’s t-test. VIP > 1 and p < 0.05 values indicated differential metabolites, which were annotated using the KEGG database to determine the related pathways. Python (Version 3.10) software package scipy.stats was used for pathway enrichment analysis, and Fisher’s exact test was employed to determine the most relevant biological pathways for experimental processing.

3. Results

3.1. The Sensory Description, Physicochemical Indicators, and Viable Counts of the Fermented Milk

The sensory description, acidity, ethanol content, and viable LAB and yeast count of the fermented milk are shown in Table 1. The acidity values of all the fermented milk samples were between 70 and 80 °T, which was within the acceptable range for consumers [17]. The acidity of KFM reached 72.34 °T without the participation of LAB, indicating that K. marxianus SY11 displayed a strong acid-producing ability. The low alcohol level in the fermented milk, ranging between 1 and 4 g/L, indicated that the yeast provided a mellow flavor to the fermentation system without causing health or other issues due to high alcohol content. The MFM exhibited a LAB level of 10.35 lg CFU/mL and a yeast content of 9.86 lg CFU/mL, both at a higher order of magnitude, showing that the two microorganisms exhibited excellent proliferation in the fermentation system. The sensory description results indicated that KFM exhibited an uneven, thin texture with bubbles, whey precipitation, an obvious alcoholic flavor, and a bland, sour taste. Therefore, K. marxianus SY11 alone was unsuitable for milk fermentation. LFM displayed a smooth, thick texture with a distinct sour aroma, showing that L. brevis NZ4 was appropriate as a common yogurt starter. MFM showed a smooth, thick texture with micro-bubbles and the moderate alcohol, ester, and sour aroma characteristic of fermented milk in western Sichuan, China. Therefore, combining K. marxianus SY11 and L. brevis NZ4 as a starter culture is suitable for producing local, uniquely flavored dairy products.
The volatile flavor compounds in the three kinds of fermented milk samples detected via GC–MS are shown in Table 2. A total of thirty-one volatile flavor compounds were identified, including eight acids, five alcohols, nine esters, three alkanes, and six other volatile components. The order of the volatile substance types and quantities in the fermented milk was MFM (16 types) > LMF (14 types) = KFM (14 types). Acids represented the main volatile components in LMF, while KFM contained mostly acids and alcohols, and MFM consisted primarily of acids and esters. These results indicated a higher aroma compound abundance when combining L. brevis and K. marxianus for milk fermentation, while the acids reacted with the alcohol to increase the ester content.
As shown in Table 3, seven common organic acids were detected in the three kinds of fermented milk. The lactic acid content was highest in MFM, followed by acetic acid, citric acid, succinic acid, oxalic acid, malic acid, and tartaric acid. MFM also displayed the highest organic acid levels. The malic acid content in KFM was 0.28 mg/100 g, significantly lower than the other two fermented milk samples. Some studies have shown that yeast degrades malic acid [18]. In addition, malic acid is an intermediate product of the TCA cycle. Other research has revealed that malic acid is completely metabolized during mixed L. brevis and K. marxianus fermentation, while acetic acid is produced in large quantities [19], which is consistent with the malic and acetic acids change trend in this study.
The amino acid determination results are shown in Table 4. A total of 17 free amino acids were detected in the three types of fermented milk. The total free amino acid content was significantly higher in MFM and KFM than in LFM and included seven essential amino acids: threonine, valine, methionine, isoleucine, leucine, phenylalanine, and lysine. These amino acids must be obtained from food due to limited synthesis in the human body. Of these, the lysine levels were the highest, reaching 30.15 mg/100 mL in MFM. Lysine is generally higher in animals than cereal proteins and is vital for promoting human growth and development and regulating metabolic balance [20]. Although arginine can be synthesized in the human body, the synthesis rate is slow. Therefore, it is also known as a semi-essential amino acid [21]. The arginine content in the three fermented milk samples ranged from 0.38 mg/100 mL to 0.51 mg/100 mL.

3.2. Metabolite Profiles and Data Analysis

3.2.1. Data Quality Control Analysis

The PCA is shown in Figure 1. The tight QC sample clustering on the PCA map indicated the stability and high repeatability of the program. The R2X values of the PCA model in positive and negative ion modes were 0.537 and 0.555, respectively, indicating high fitting accuracy. The three experimental samples showed significant clustering in the PCA plots in both positive and negative ion modes, indicating excellent repeatability and minimal differences in the sample group.

3.2.2. Overview of the Metabolites

Metabolomics detected a total of 9384 peaks in both positive and negative ion modes. Of these, 901 metabolites were annotated via the KEGG and HMDB databases. After removing data with missing values > 20% and QC validation RSD values < 30%, 885 metabolites were ultimately annotated, including 627 for positive ions and 258 for negative ions. KEGG annotated a total of 295 metabolic products, mainly including phospholipids (20 types), fatty acids (14 types), carboxylic acids (12 types), and glycerol phosphate choline (12 types). HMDB annotated a total of 667 metabolites, primarily including lipids (282 species), organic acids (136 species), oxygen-containing organic compounds, such as monosaccharides (70 species), and organic heterocyclic compounds, such as nucleic acids, vitamins, antibiotics (69 species), and benzoids (36 species).

3.3. Differential Metabolite Analysis

3.3.1. Overview of the Differences between the Samples

PLS-DA was used to analyze the differences between the samples, as shown in Figure 2A. Positive ion mode revealed significant differences between MFM, LFM, and KFM. However, the distance between LFM and KFM was relatively small, possibly due to the substantial differences between MFM and the other two types. As shown in Figure 2C, considerable differences were also evident between MFM and LFM or KFM in negative ion mode, while an intersection was apparent between LFM and KFM. This may be because the differential metabolites produced by KFM and LFM are mostly volatile substances and cannot be fully reflected by HPLC–MS. A permutation test confirmed the PLS-DA model stability (Figure 3B,D). In positive ion mode, PLS-DA produced a cumulative R2X value of 0.7177, a cumulative R2Y value of 0.991, and a cumulative Q2 value of 0.4027. Negative ion mode produced a cumulative R2X value of 0.754, a cumulative R2Y value of 0.988, and a cumulative Q2 value of 0.3492, indicating that the PLS-DA models displayed excellent fitting and high predictive ability.

3.3.2. Differential Metabolite and Pathway Screening

OPLS-DA was used to obtain the VIP and Student’s test p-value. VIP values > 1 and p < 0.05 indicated differential metabolites, of which 204 were evident between MFM and LFM, 163 between MFM and KFM, and 151 between LFM and KFM. The metabolites related to the fermented milk flavor (including organic acids, amino acids, sugars, and fatty acids) were screened. The 50 most abundant metabolites were selected to prepare the cluster diagram, as shown in Figure 3. The results showed that the flavor-related metabolites of MFM and KFM were similar, while significant differences were evident between those of LFM and KFM.
The relevant metabolites were searched in the KEGG database and enriched in the KEGG pathway. As shown in Figure 4, five metabolic pathways were highly enriched in the test group, including glycerophospholipid metabolism, arginine biosynthesis, beta-alanine metabolism, arginine and proline metabolism, and benzoxazinoid biosynthesis. A glycerophospholipid metabolic pathway map was drawn using Illustrator software (Version CS5) to further analyze the related upstream and downstream metabolites.
As shown in Figure 5, the annotated metabolites annotated in this pathway included 1-Acyl-sn-glyceeo-3-phosphocholine, CDP-diacylglycerol, and phosphatidyl-L-serine. Phosphatidyl-L-serine was significantly upregulated in MFM compared with LFM and KFM, while CDP-diacylglycerol and phosphatidyl-glycerophosphate were substantially higher than in KFM. The 1-Acyl-sn-glyceeo-3-phospholine level in MFM was significantly lower than in LFM and KFM, suggesting it represented a precursor of this metabolic pathway. The phosphate tidylcholine and 1,2-diacyl-sn-glycerol-3P levels were considerably lower in MFM than in KFM. These two substances may play an intermediary role in the metabolic pathway.

4. Discussion

Fermented dairy products are highly popular with consumers worldwide. Yeast can change the quality during the fermentation and ripening of these products. Its high protein decomposition ability can produce flavor-affecting precursors as well as alcohol and carbon dioxide that reduce the pH and inhibit harmful bacterial growth [22,23]. Therefore, this study used mixed L. brevis and K. marxianus fermentation to prepare alcohol-flavored fermented milk.
MFM and LFM displayed high LAB levels at about 1010 CFU/mL. The yeast content was significantly higher in MFM than in KFM, while the ethanol level was lower. This may be because the ethanol in MFM is utilized by LAB or esterified with other substances. Some studies have shown that ethanol and lactic acid esterification can reduce the damage to their own cells to a certain extent [24]. The types and relative content of the esters in the volatile MFM flavor compounds were higher than in KFM and LFM, further indicating a higher number of esterification reactions between the acids and alcohols in MFM. In addition, lactic acid has a chelating effect on cations and can reduce their impact on yeast cells [25,26]. LAB and yeast interact directly or indirectly. Wang [27] revealed that K. marxianus significantly promoted L. casei proliferation, while L. lactis substantially encouraged K. marxianus growth.
The free amino acid content was significantly higher in MFM and KFM than in LFM. The nutritional factors required for LAB growth were relatively abundant, with many amino acids consumed during fermentation. Moreover, the free amino acid content was low in LFM since its protein hydrolase enzyme system was not as rich as MFM, which exhibited a higher amino acid level due to protein hydrolyzation by K. marxianus. Since the amino acids produced by K. marxianus were used by L. brevis to maintain growth and metabolism, the free amino acid content was slightly lower in MFM than in KFM. The analysis of the basic physical and chemical indexes showed that L. brevis and K. marxianus eventually improved the quality of fermented milk by promoting cell growth and causing changes in other indexes.
The flavor substance sources in fermented milk are mostly represented by fat, lactose, and protein as precursors, producing amino acids, sulfur-containing amino acids, peptides, lipids, sugars, and other substances under the action of microorganisms. Since these substances further constitute the flavor components in fermented milk [28], amino acids and their derivatives, organic acids, sugars, and fatty acids are important analytical targets. In this study, the overall abundance of the flavor-related metabolites was higher in KFM than in LFM and included N-acetyltryptophan, methionyl-threonine, and cysteinyl-proline, most of which represented amino acid derivatives. Studies have shown that yeast increases the amino acids in the mixed fermentation system to maintain LAB growth and metabolism [29]. Using transcriptome technology, some studies found significant changes in the biosynthesis or degradation and the transport and metabolism of amino acids and other gene pathways in milk fermented with LAB and yeast [30]. From genotype to phenotype, LAB and yeast may interact using amino acids as a medium during the co-fermentation process. In addition, significant changes in the amino acid metabolism pathway were evident in fermented sourdough, indicating that this phenomenon is widespread and not limited to fermented dairy products [31]. Nopaline acid, neryl rhamnosyl-glucoside, oleyl alcohol, sphingosine, 19-hydroxycinnamyl 19-glucoside blumenol CO-[rhamnosyl-(1->6)-glucoside], and other substance levels, which mostly represented glycosides, were higher in LFM than in MFM. The glycosides produced via L. brevis metabolism provided the carbon source for K. marxianus growth. Kadyan [32] used non-targeted metabolomics to determine the key L. plantarum and K. lactis metabolites during single and mixed fermentation. The results indicated a lower sugar content in the mixed fermentation system, possibly improving the sugar utilization rate. Liu [33] revealed that mixed L. paracasei and K. marxianus fermentation significantly changed sucrose metabolism, glycolysis, and other pathways by stimulating sugar utilization by microorganisms. The flavor metabolites in MFM included N-acetyl-L-glutamic acid, cysteine-serine, glaucarubin, maleic acid, and serotonin, with fatty acid metabolism-related substances being the most abundant at a total of 97 types. Tian [34] tested the lipase activity of five yeast products, revealing that K. maximus lipase activity was the highest, and speculated that yeast participated in protein and fat metabolism, accelerating decomposition and producing various fatty and amino acids. This represents the main source of many flavor substances in fermented foods [35]. Other studies have shown that mixed milk fermentation using LAB and yeast increased the total fatty acid and free fatty acid content and the degree of protein hydrolysis, which is consistent with the results of this study [36].
The selected flavor-related differential metabolites in the KEGG metabolic pathway showed that glycerophospholipid metabolism displayed the highest enrichment. Eight test group substances were enriched into the pathway. Glycerophospholipid metabolism included two sub-pathways, namely phosphatidylcholine and phosphatidylethanolamine biosynthesis, while the substances in the test group were mainly enriched in phosphatidylethanolamine biosynthesis, including four compounds, as shown in Figure 5, of which phosphatidylserine played the most significant role. Phosphatidylserine is formed by phosphatidylserine synthase with L-serine and CDP-diacylglycerol as precursors and is widely found in bacterial, yeast, and mammalian cell membranes. It represents the primary acidic phospholipid in human cell membranes and usually forms part of the cell membrane lobule. It accounts for about 2–20% of the total phospholipid mass in human plasma and cell membranes [37,38]. Some studies have shown that phosphatidylserine displays excellent anti-inflammatory and anti-depression properties and improves cognitive brain function [39]. Zhang [40] screened 94 lipid-related differential metabolites in sheep’s milk using lipomics. In addition, glycerophospholipid metabolism represents a vital fatty acid metabolism pathway and plays a crucial role in product flavor formation. Zhang [41] measured the flavor changes during cheese fermentation in Inner Mongolia, screening 37 different metabolites and six key metabolic pathways. The metabolic pathway was closely related to lipid catabolism, indicating that the catabolic lipid activity during the cheese production process affected its flavor formation. Chen [42] analyzed the flavor precursor changes in purple-leaf tea during different production processes using a lipomics system. The results revealed alpha-linolenic acid metabolism and glycerolipid metabolism as critical lipid metabolism pathways in the raw tea and subsequent products, respectively.
Of the other metabolic pathways noted in this experiment, arginine biosynthesis, beta-alanine metabolism, arginine and proline metabolism, and benzoxazinoid biosynthesis exhibited a higher abundance, most of which were related to amino acid synthesis and metabolism. Yang [43] revealed that the total amino acid content initially increased when mixing Lactobacillus sanfranciscensis (L. sanfranciscensis) and Saccharomyces cerevisiae (S. cerevisiae) to ferment sourdough, followed by a decline. It is speculated that the two microorganisms consume a large number of amino acids during the initial fermentation stage, while the lower pH in the dough during the later fermentation stage increases the protease activity in the flour, promoting amino acid degradation. Furthermore, S. cerevisiae and L. sanfranciscensis were exposed to metabolic activity, such as protein hydrolysis, where the protein was decomposed into amino acids, increasing the amino acid content in the dough [44,45].
In conclusion, mixed L. brevis and K. marxianus fermentation promotes the metabolism of milk proteins, carbohydrates, and fats. The significant enrichment of differential metabolites in the glycerophospholipid metabolism pathway can explain the mechanism behind the changes in fermented milk quality caused by mixed fermentation.

5. Conclusions

This study reveals the quality formation in fermented milk inoculated with Lactobacillus brevis NZ4 and Kluyveromyces marxianus SY11, the dominant microorganisms isolated from traditional dairy products. MFM displays a better overall quality with good sensory characteristics, more organic acid types, higher levels of free amino acids and volatile substance esters, and lower acidity and ethanol concentrations than LFM and KFM. The differential metabolites in MFM are mainly enriched in pathways such as glycerophospholipid metabolism, arginine biosynthesis, beta-alanine metabolism, arginine, and proline metabolism. The mixed fermentation of these two strains significantly changes the flavor-related substance metabolism pathway. The results will pave the way for industrial production and flavor regulation of local special fermented milk.

Author Contributions

Conceptualization, X.-L.A., methodology, X.-L.A., S.-L.L. and Y.Y.; validation, Y.-M.L., H.-Y.K., H.-L.L. and T.H.; formal analysis, X.-L.A. and Y.-M.L.; resources, L.-K.Z.; data curation, X.-Y.L.; writing—original draft preparation, Y.-M.L., H.-Y.K., H.-L.L. and T.H.; writing—review and editing, X.-L.A., S.-L.L., S.-J.C., Y.Y., X.-L.A. and Y.-M.L. contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This study is supported by the Natural Science Foundation of Sichuan Province of China, 2022NSFSC0115.

Data Availability Statement

The data presented in this study are available in the article.

Acknowledgments

This research was financially supported by the Natural Science Foundation of Sichuan Province of China (grant number 2022NSFSC0115). We thank the staff at the Microbiological Team, College of Food Science, Sichuan Agricultural University, for their assistance.

Conflicts of Interest

The authors declare that they have no known competing financial interest or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Şafak, H.; Gün, İ.; Tudor Kalit, M.; Kalit, S. Physico-Chemical, Microbiological and Sensory Properties of Water Kefir Drinks Produced from Demineralized Whey and Dimrit and Shiraz Grape Varieties. Foods 2023, 12, 1851. [Google Scholar] [CrossRef] [PubMed]
  2. Zhong, Z.; Hou, Q.; Kwok, L.; Zheng, Y.; Sun, M.B.; Zhang, H. Bacterial microbiota compositions of naturally fermented milk are shapedby both geographic origin and sample type. J. Dairy Sci. 2016, 99, 7832–7841. [Google Scholar] [CrossRef] [PubMed]
  3. Dalmasso, A.; Civera, T.; Pattono, D.; Cardazzo, B.; Bottero, M.T. Characterization of microbiota in Plaisentif cheese by high-throughput sequencing. LWT-Food Sci. Technol. 2016, 69, 490–496. [Google Scholar] [CrossRef]
  4. Ni, H.J. Identification and Biodiversity of Yeasts Isolated from Traditional Fermented Milk Products in Xinjiang and Qinghai of China; Inner Mongolia Agricultural University: Hohhot, China, 2009. [Google Scholar]
  5. Saleena, L.A.K.; Phing, P.L.; Gan, R.Y.; Al-Nabulsi, A.; Osaili, T.; Kamal-Eldin, A.; Ayyash, M. Fermented dairy products from Middle Eastern and Northern African (MENA) countries: Insight on production and physiochemical characteristics. Int. Dairy J. 2023, 141, 105614. [Google Scholar] [CrossRef]
  6. Mehraban Sangatash, M.; Azadi, H.; Ehtiati, A. Evaluation of physicochemical properties of yogurt fortified with inactive baker’s yeast. J. Food Sci. Technol. 2023, 20, 210–229. [Google Scholar]
  7. Peng, X.; Yue, Q.; Chi, Q.; Liu, Y.; Tian, T.; Dai, S.; Jiang, L. Microbial Diversity and Flavor Regularity of Soy Milk Fermented Using Kombucha. Foods 2023, 12, 884. [Google Scholar] [CrossRef] [PubMed]
  8. Fang, L.; Wang, W.; Dou, Z.; Chen, J.; Meng, Y.; Cai, L.; Li, Y. Effects of mixed fermentation of different lactic acid bacteria and yeast on phytic acid degradation and flavor compounds in sourdough. LWT 2023, 174, 114438. [Google Scholar] [CrossRef]
  9. Proust, L.; Haudebourg, E.; Sourabié, A.; Pedersen, M.; Besançon, I.; Monnet, V.; Juillard, V. Multi-omics approach reveals how yeast extract peptides shape Streptococcus thermophilus metabolism. Appl. Environ. Microbiol. 2020, 86, e01446-20. [Google Scholar] [CrossRef]
  10. Jia, H.; Wu, Z.; Tan, J.; Wu, S.; Yang, C.; Raza, S.H.A.; Wang, M.; Song, G.; Shi, Y.; Zan, L.; et al. Lnc-TRTMFS promotes milk fat synthesis via the miR-132x/RAI14/mTOR pathway in BMECs. J. Anim. Sci. 2023, 101, skad218. [Google Scholar] [CrossRef]
  11. Arakawa, K.; Miyamoto, M.; Miyamoto, T. Interaction between lactic acid bacteria and yeasts in airag, an alcoholic fermented milk. Anim. Sci. J. 2013, 84, 66–74. [Google Scholar]
  12. Sun, Z.H. Characterization of Dominant Bacteria and Volatile Flavors in Koumiss and Study on Its Metabolomics; Zhejiang Gongshang University: Hangzhou, China, 2019. [Google Scholar]
  13. Ralli, E.; Amargianitaki, M.; Manolopoulou, E.; Misiak, M.; Markakis, G.; Tachtalidou, S.; Spyros, A. NMR spectroscopy protocols for food metabolomics applications. In Metabolic Profiling: Methods and Protocols; Humana Press: New York, NY, USA, 2018; pp. 203–211. [Google Scholar]
  14. Wang, X.C. Isolation and Identification of Non-Saccharomyces cerevisiae from East Foot of Helanshan Moutains of Ningxia and Fermentation Characteristics; Ningxia University: Yinchuan, China, 2016. [Google Scholar]
  15. Meng, F.B.; Zhou, L.; Li, J.J.; Li, Y.C.; Wang, M.; Zou, L.H.; Chen, W.J. The combined effect of protein hydrolysis and Lactobacillus plantarum fermentation on antioxidant activity and metabolomic profiles of quinoa beverage. Food Res. Int. 2022, 157, 111416. [Google Scholar] [CrossRef]
  16. Fang, X.; Liu, Y.; Xiao, J.; Ma, C.; Huang, Y. GC–MS and LC-MS/MS metabolomics revealed dynamic changes of volatile and non-volatile compounds during withering process of black tea. Food Chem. 2023, 410, 135396. [Google Scholar] [CrossRef] [PubMed]
  17. Peng, J.; Ma, L.; Kwok, L.Y.; Zhang, W.; Sun, T. Untargeted metabolic footprinting reveals key differences between fermented brown milk and fermented milk metabolomes. J. Dairy Sci. 2022, 105, 2771–2790. [Google Scholar] [CrossRef] [PubMed]
  18. Wang, J.L.; Yan, Y.C.; Li, Q.Y.; Zhang, F.S.; Li, X.; Lu, Z.Q. Research progress of biodegradable citric acid and its influencing factors. Mod. Food Sci. Technol. 2022, 38, 347–357. [Google Scholar]
  19. Herrero, M.; Garcia, L.A.; Diaz, M. Organic Acids in Cider with Simultaneous Inoculation of Yeast and Malolactic Bacteria: Effect of Fermentation Temperature. J. Inst. Brew. 1999, 105, 229–232. [Google Scholar] [CrossRef]
  20. Wang Li Luo, H.X.; Li, S.R.; Chi, Q.B.; Wang, H.H. Amino acid composition and nutritional evaluation of different Hypsizygus marmoreus samples. Sci. Technol. Food Ind. 2016, 37, 338–341. [Google Scholar]
  21. Tian, H.M. Preliminary Optimization of L-arginine Fermentation Technology on the Genetically Engineered Bacteria Corynebacterium glutamicum; South China University of Technology: Guangzhou, China, 2016. [Google Scholar]
  22. Liu, Y.; Cao, Y.; Yohannes Woldemariam, K.; Zhong, S.; Yu, Q.; Wang, J. Antioxidant effect of yeast on lipid oxidation in salami sausage. Front. Microbiol. 2023, 13, 1113848. [Google Scholar] [CrossRef] [PubMed]
  23. Zhang, K.; Zhang, T.T.; Guo, R.R.; Ye, Q.; Zhao, H.L.; Huang, X.H. The Regulation of Key Flavor of Traditional Fermented Food by Microbial Metabolism: A Review. Food Chem. X 2023, 19, 100871. [Google Scholar] [CrossRef]
  24. Sieuwerts, S.; Bron, P.A.; Smid, E.J. Mutually stimulating interactions between lactic acid bacteria and Saccharomyces cerevisiae in sourdough fermentation. LWT 2018, 90, 201–206. [Google Scholar] [CrossRef]
  25. Wang, Y.; Han, J.; Ren, Q.; Liu, Z.; Zhang, X.; Wu, Z. The Involvement of Lactic Acid Bacteria and Their Exopolysaccharides in the Biosorption and Detoxication of Heavy Metals in the Gut. Biol. Trace Elem. Res. 2023, 1–14. [Google Scholar] [CrossRef]
  26. Abduxukur, D.; Tursuntay, A.; Zhu, X.; Wang, X.; Rahman, N. Antioxidant Capacity of Lactic Acid Bacteria and Yeasts from Xinjiang Traditional Fermented Dairy Products. Fermentation 2023, 9, 639. [Google Scholar] [CrossRef]
  27. Wang, X.B.; Wu, Y.; Wang, L.; Miao, S.; Zou, W.; Li, Z.H.; Li, W.W. Interaction between lactic acid bacteria and yeast during koumiss fermentation. Food Sci. 2015, 36, 125–129. [Google Scholar]
  28. Wang, M.S.; Fan, M.; Zheng, A.R.; Wei, C.K.; Liu, D.H.; Thaku, K.; Wei, Z.J. Characterization of a fermented dairy, sour cream: Lipolysis and the release profile of flavor compounds. Food Chem. 2023, 423, 136299. [Google Scholar] [CrossRef] [PubMed]
  29. Ponomarova, O.; Gabrielli, N.; Sévin, D.C.; Mülleder, M.; Zirngibl, K.; Bulyha, K.; Patil, K.R. Yeast creates a niche for symbiotic lactic acid bacteria through nitrogen overflow. Cell Syst. 2017, 5, 345–357. [Google Scholar] [CrossRef] [PubMed]
  30. Abrahams, A.M. Proteomics and transcriptomics and their application in fermented foods. In Indigenous Fermented Foods for the Tropics; Academic Press: Cambridge, MA, USA, 2023; pp. 377–391. [Google Scholar]
  31. Zhang, G.H.; Wang, W.; Tu, J.; Zhang, W.Z.; He, G.Q. Analysis of microbial metabolism mechanism in traditional sourdough fermentation based on mac rotranscriptomics technology. J. Chin. Cereals Oils Assoc. 2019, 34, 10–16. [Google Scholar]
  32. Kadyan, S.; Rashmi, H.M.; Pradhan, D.; Kumari, A.; Chaudhari, A.; Deshwal, G.K. Effect of lactic acid bacteria and yeast fermentation on antimicrobial, antioxidative and metabolomic profile of naturally carbonated probiotic whey drink. LWT 2021, 142, 111059. [Google Scholar] [CrossRef]
  33. Liu, N.; Qin, L.; Mazhar, M.; Miao, S. Integrative transcriptomic-proteomic analysis revealed the flavor formation mechanism and antioxidant activity in rice-acid inoculated with Lactobacillus paracasei and Kluyveromyces marxianus. J. Proteom. 2021, 238, 104158. [Google Scholar] [CrossRef] [PubMed]
  34. Tian, Y.C. Study on the Influence of the Yeast on Lipids Metabolism in Traditional Dairy Products; Dalian Polytechnic University: Dalian, China, 2013. [Google Scholar]
  35. Wang, X.; Li, W.; Mahsa, G.C.; Zhang, C.; Ma, K.; Rui, X.; Li, W. Co-cultivation effects of Lactobacillus helveticus SNA12 and Kluveromyces marxiensis GY1 on the probiotic properties, flavor, and digestion in fermented milk. Food Res. Int. 2023, 169, 112843. [Google Scholar] [CrossRef]
  36. Li, L. The Effects of Kluyveromyces marxianus to Lipid Metabolism in Fermented Milk; Jilin Agricultural University: Changchun, China, 2015. [Google Scholar]
  37. Calianese, D.C.; Birge, R.B. Biology of phosphatidylserine (PS): Basic physiology and implications in immunology, infectious disease, and cancer. Cell Commun. Signal. 2020, 18, 41. [Google Scholar] [CrossRef]
  38. Partoazar, A.; Seyyedian, Z.; Zamanian, G.; Saffari, P.M.; Muhammadnejad, A.; Dehpour, A.R.; Goudarzi, R. Neuroprotective phosphatidylserine liposomes alleviate depressive-like behavior related to stroke through neuroinflammation attenuation in the mouse hippocampus. Psychopharmacology 2021, 238, 1531–1539. [Google Scholar] [CrossRef]
  39. Liu, B.Y. Mechanistic Study of Phosphatidylserine Emulsion to Improve Brain Cognitive Function and Preliminary Investigation of Method Parameters; Liaocheng University: Liaocheng, China, 2022. [Google Scholar]
  40. Zhang, R.; Zhu, Z.; Jia, W. Molecular mechanism associated with the use of magnetic fermentation in modulating the dietary lipid composition and nutritional quality of goat milk. Food Chem. 2022, 366, 130554. [Google Scholar] [CrossRef]
  41. Zhang, X.; Zheng, Y.; Feng, J.; Zhou, R.; Ma, M. Integrated metabolomics and high-throughput sequencing to explore the dynamic correlations between flavor related metabolites and bacterial succession in the process of Mongolian cheese production. Food Res. Int. 2022, 160, 111672. [Google Scholar] [CrossRef] [PubMed]
  42. Chen, X.; Wang, P.; Wei, M.; Lin, X.; Gu, M.; Fang, W.; Ye, N. Lipidomics analysis unravels changes from flavor precursors in different processing treatments of purple-leaf tea. J. Sci. Food Agric. 2022, 102, 3730–3741. [Google Scholar] [CrossRef] [PubMed]
  43. Yang, H.Y. Intraspecific Diversity of Saccharomyces cerevisiae and Lactobacillus sanfranciscensis Isolates from Chinese Traditional Sourdough and Their Interaction; Zhejiang University: Hangzhou, China, 2018. [Google Scholar]
  44. Liu, Y.; Danial, M.; Liu, L.; Sadiq, F.A.; Wei, X.; Zhang, G. Effects of Co-Fermentation of Lactiplantibacillus plantarum and Saccharomyces cerevisiae on Digestive and Quality Properties of Steamed Bread. Foods 2023, 12, 3333. [Google Scholar] [CrossRef] [PubMed]
  45. Fu, W.; Jia, X.; Liu, C.; Meng, X.; Zhang, K.; Tao, S.; Xue, W. Sourdough yeast-bacteria interactions results in reduced immunogenicity by increasing depolymerization and hydrolysis of gluten. Innov. Food Sci. Emerg. Technol. 2023, 84, 103281. [Google Scholar] [CrossRef]
Figure 1. Principal component analysis (PCA) score plots of 3 groups of samples in (A) positive [electrospray ionization (ESI+)] (R2X = 0.537) and (B) negative (ESI−) (R2X = 0.555) ion modes. MFM: milk fermented with L. brevis NZ4 + K. marxianus SY11 (1:1); LFM: milk fermented with L. brevis NZ4; KFM: milk fermented with K. marxianus SY11; and QC: quality control; the shadow represents the confidence interval.
Figure 1. Principal component analysis (PCA) score plots of 3 groups of samples in (A) positive [electrospray ionization (ESI+)] (R2X = 0.537) and (B) negative (ESI−) (R2X = 0.555) ion modes. MFM: milk fermented with L. brevis NZ4 + K. marxianus SY11 (1:1); LFM: milk fermented with L. brevis NZ4; KFM: milk fermented with K. marxianus SY11; and QC: quality control; the shadow represents the confidence interval.
Foods 12 03704 g001
Figure 2. (A,B) Partial least squares discriminant analysis (PLS-DA) of scores and permutation test plots for the fermented milk analyzed in the positive ion mode. (C,D) Partial least squares discriminant analysis (PLS-DA) of scores and permutation test plots for the fermented milk analyzed in the negative ion mode; the shadow represents the confidence interval.
Figure 2. (A,B) Partial least squares discriminant analysis (PLS-DA) of scores and permutation test plots for the fermented milk analyzed in the positive ion mode. (C,D) Partial least squares discriminant analysis (PLS-DA) of scores and permutation test plots for the fermented milk analyzed in the negative ion mode; the shadow represents the confidence interval.
Foods 12 03704 g002
Figure 3. Heatmap of the 50 flavor-related differential metabolites identified in fermented milk samples.
Figure 3. Heatmap of the 50 flavor-related differential metabolites identified in fermented milk samples.
Foods 12 03704 g003
Figure 4. The metabolome view map of the significant metabolic pathways characterized in the fermented milk samples. Circles represent metabolic pathways. The X-axis represents the pathway impact, and the Y-axis represents pathway enrichment. Larger sizes and darker colors indicate greater pathway enrichment and higher pathway impact values, respectively.
Figure 4. The metabolome view map of the significant metabolic pathways characterized in the fermented milk samples. Circles represent metabolic pathways. The X-axis represents the pathway impact, and the Y-axis represents pathway enrichment. Larger sizes and darker colors indicate greater pathway enrichment and higher pathway impact values, respectively.
Foods 12 03704 g004
Figure 5. Partial glycerophospholipid metabolic pathway map; red dots represent a relative increase in the content of metabolites, while blue dots represent a relative decrease in the content of metabolites; the upper arrow represents the increase of the metabolite, the lower arrow represents the decrease of the metabolite, and the flat arrow represents no significant change.
Figure 5. Partial glycerophospholipid metabolic pathway map; red dots represent a relative increase in the content of metabolites, while blue dots represent a relative decrease in the content of metabolites; the upper arrow represents the increase of the metabolite, the lower arrow represents the decrease of the metabolite, and the flat arrow represents no significant change.
Foods 12 03704 g005
Table 1. Sensory description, physicochemical indicators, and microbial count of the fermented milk.
Table 1. Sensory description, physicochemical indicators, and microbial count of the fermented milk.
ItemMFMLFMKFM
Acidity (°T)76.93 ± 0.48 a78.25 ± 0.59 a72.34 ± 0.83 b
Ethanol content (g/L)3.22 ± 0.03 b1.56 ± 0.12 c3.99 ± 0.09 a
LAB (lg CFU/mL)10.35 ± 0.27 a10.6 ± 0.13 a
Yeast (lg CFU/mL)9.86 ± 0.33 a9.09 ± 0.22 b
Sensory descriptionAppearanceThick, no wheyThick, a small amount of whey Thin, a small amount of whey
TextureMicro-bubble, smoothNo bubble, smoothBubble, unsmooth
FlavorModerate alcohol, ester, and sour aroma Obvious sour aromaObvious alcohol flavor and mild acid aroma
Note: MFM: milk fermented with L. brevis NZ4 + K. marxianus SY11 (1:1); LFM: milk fermented with L. brevis NZ4; KFM: milk fermented with K. marxianus SY11; and a, b, and c represent the significant differences in the same column, p < 0.05.
Table 2. Relative content of volatile flavor compounds in the fermented milk.
Table 2. Relative content of volatile flavor compounds in the fermented milk.
Volatile SubstancesRelative Content (%)
MFMLFMKFM
AcidsAcetic acid25.16 ± 2.6517.72 ± 3.4928.28 ± 1.02
Butanoic acid6.51 ± 0.72
Pentanoic acid6.16 ± 0.909.10 ± 0.7412.27 ± 0.68
n-Decanoic acid13.72 ± 0.5111.76 ± 1.4915.19 ± 2.48
Octanoic acid12.21 ± 1.2411.60 ± 1.40
Propanedioic acid6.68 ± 0.52
Hexanoic acid, 4-methyl-5.29 ± 0.28
6-Hydroxy-2-naphthoic acid2.57 ± 0.23
AlcoholsEthanol10.33 ± 2.81
3-Nonanol, 3-methyl-7.27 ± 0.76
3-Pentanol, 3-methyl-4.02 ± 0.20
1-Propanol, 3-chloro-0.28 ± 0.04
4-Methoxy-4-methyl-2-pentanol5.46 ± 0.27
Esters2-Chloroethyl benzoate13.18 ± 2.4310.41 ± 1.1316.30 ± 2.49
Butanoic acid, 3-methyl-2.54 ± 1.76
Oxalic acid, cyclobutyl hexyl ester2.20 ± 0.46
1-Propen-2-ol, formate1.61 ± 0.04
Benzeneacetic acid, α-oxo-, methyl ester5.67 ± 0.57
Ethyl hydrogen oxalate8.40 ± 0.48
Ethyl oxamate10.27 ± 0.57
Ethanol, 2-methoxy-, acetate3.64 ± 0.57
2-Propenoic acid, butyl ester4.07 ± 0.45
TaxanesPentane, 2,3,3,4-tetramethyl-4.59 ± 0.57
Cyclobutane, methoxy-0.70 ± 0.27
3-Pentanamine0.79 ± 0.09
QthersAcetoin6.98 ± 2.5322.77 ± 2.4113.88 ± 2.42
Acetone2.32 ± 0.892.14 ± 0.883.74 ± 0.24
2-Propanamine, N-methyl-0.25 ± 0.18
Toluene0.77 ± 0.17
Formic acid hydrazide2.50 ± 0.07
Acetamide, N-methyl-1.72 ± 0.13
Table 3. Organic acid content of the fermented milk.
Table 3. Organic acid content of the fermented milk.
Type (mg/100 g)MFMLFMKFM
Oxalic acid17.58 ± 0.18 c21.44 ± 0.14 b23.04 ± 0.14 a
Tartaric acid1.89 ± 0.01 c2.02 ± 0.004 b2.52 ± 0.01 a
Malic acid6.57 ± 0.03 b8.51 ± 0.18 a0.28 ± 0.02 c
Lactic acid224.29 ± 1.67 a200.05 ± 1.50 b190.83 ± 0.57 c
Acetic acid88.50 ± 0.73 b91.74 ± 0.52 a92.93 ± 1.53 a
Citric acid55.88 ± 0.37 b60.92 ± 0.97 a59.71 ± 0.41 a
Succinic acid27.76 ± 1.21 a25.26 ± 1.10 a26.82 ± 2.20 a
Total422.46 ± 3.33 a409.94 ± 2.08 b396.13 ± 1.34 c
Note: a, b, and c represent the significant differences in the same column, p < 0.05.
Table 4. Free amino acid content of the fermented milk.
Table 4. Free amino acid content of the fermented milk.
TypeFree Amino Acid Content (mg/100 mL)
MFMLFMKFM
Aspartic acid10.81 ± 0.16 a7.13 ± 0.10 c9.89 ± 0.29 b
Threonine2.49 ± 0.12 c4.10 ± 0.15 a3.14 ± 0.23 b
Serine1.37 ± 0.10 a0.70 ± 0.23 b1.08 ± 0.09 ab
Glutamic acid40.22 ± 0.28 a35.67 ± 0.50 b38.95 ± 0.87 a
Glycine6.63 ± 0.04 b7.03 ± 0.05 a6.79 ± 0.17 ab
Cysteine0.61 ± 0.26 a0.28 ± 0.01 a0.91 ± 0.37 a
Alanine0.30 ± 0.03 b0.64 ± 0.18 ab0.84 ± 0.15 a
Valine0.55 ± 0.21 b0.45 ± 0.01 b0.90 ± 0.05 a
Methionine2.45 ± 0.71 a2.15 ± 1.04 a3.90 ± 1.10 a
Isoleucine4.03 ± 0.15 b4.08 ± 0.20 b5.09 ± 0.03 a
Leucine6.11 ± 0.07 b4.36 ± 0.06 c6.45 ± 0.08 a
Tyrosine3.10 ± 1.18 a3.58 ± 0.21 a3.90 ± 0.82 a
Phenylalanine2.69 ± 0.05 ab2.27 ± 0.05 b2.98 ± 0.28 a
Lysine30.15 ± 1.22 a24.08 ± 0.30 b28.72 ± 0.61 a
Histidine2.21 ± 0.15 a2.15 ± 0.49 a2.41 ± 0.01 a
Arginine0.38 ± 0.04 a0.45 ± 0.17 a0.51 ± 0.13 a
Proline9.10 ± 0.23 a8.89 ± 0.12 a9.65 ± 0.47 a
Total123.19 ± 0.62 a108.00 ± 1.82 b126.12 ± 3.82 a
Note: a, b, and c represent the significant differences in the same column, p < 0.05.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ao, X.-L.; Liao, Y.-M.; Kang, H.-Y.; Li, H.-L.; He, T.; Zou, L.-K.; Liu, S.-L.; Chen, S.-J.; Yang, Y.; Liu, X.-Y. Untargeted Metabolomics and Physicochemical Analysis Revealed the Quality Formation Mechanism in Fermented Milk Inoculated with Lactobacillus brevis and Kluyveromyces marxianus Isolated from Traditional Fermented Milk. Foods 2023, 12, 3704. https://doi.org/10.3390/foods12193704

AMA Style

Ao X-L, Liao Y-M, Kang H-Y, Li H-L, He T, Zou L-K, Liu S-L, Chen S-J, Yang Y, Liu X-Y. Untargeted Metabolomics and Physicochemical Analysis Revealed the Quality Formation Mechanism in Fermented Milk Inoculated with Lactobacillus brevis and Kluyveromyces marxianus Isolated from Traditional Fermented Milk. Foods. 2023; 12(19):3704. https://doi.org/10.3390/foods12193704

Chicago/Turabian Style

Ao, Xiao-Lin, Yi-Mo Liao, Hai-Yan Kang, Hong-Li Li, Tian He, Li-Kou Zou, Shu-Liang Liu, Shu-Juan Chen, Yong Yang, and Xing-Yan Liu. 2023. "Untargeted Metabolomics and Physicochemical Analysis Revealed the Quality Formation Mechanism in Fermented Milk Inoculated with Lactobacillus brevis and Kluyveromyces marxianus Isolated from Traditional Fermented Milk" Foods 12, no. 19: 3704. https://doi.org/10.3390/foods12193704

APA Style

Ao, X. -L., Liao, Y. -M., Kang, H. -Y., Li, H. -L., He, T., Zou, L. -K., Liu, S. -L., Chen, S. -J., Yang, Y., & Liu, X. -Y. (2023). Untargeted Metabolomics and Physicochemical Analysis Revealed the Quality Formation Mechanism in Fermented Milk Inoculated with Lactobacillus brevis and Kluyveromyces marxianus Isolated from Traditional Fermented Milk. Foods, 12(19), 3704. https://doi.org/10.3390/foods12193704

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop