- Research article
- Open Access
Lactation-related metabolic mechanism investigated based on mammary gland metabolomics and 4 biofluids’ metabolomics relationships in dairy cows
BMC Genomics volume 18, Article number: 936 (2017)
Lactation is extremely important for dairy cows; however, the understanding of the underlying metabolic mechanisms is very limited. This study was conducted to investigate the inherent metabolic patterns during lactation using the overall biofluid metabolomics and the metabolic differences from non-lactation periods, as determined using partial tissue-metabolomics. We analyzed the metabolomic profiles of four biofluids (rumen fluid, serum, milk and urine) and their relationships in six mid-lactation Holstein cows and compared their mammary gland (MG) metabolomic profiles with those of six non-lactating cows by using gas chromatography-time of flight/mass spectrometry.
In total, 33 metabolites were shared among the four biofluids, and 274 metabolites were identified in the MG tissues. The sub-clusters of the hierarchical clustering analysis revealed that the rumen fluid and serum metabolomics profiles were grouped together and highly correlated but were separate from those for milk. Urine had the most different profile compared to the other three biofluids. Creatine was identified as the most different metabolite among the four biofluids (VIP = 1.537). Five metabolic pathways, including gluconeogenesis, pyruvate metabolism, the tricarboxylic acid cycle (TCA cycle), glycerolipid metabolism, and aspartate metabolism, showed the most functional enrichment among the four biofluids (false discovery rate < 0.05, fold enrichment >2). Clear discriminations were observed in the MG metabolomics profiles between the lactating and non-lactating cows, with 54 metabolites having a significantly higher abundance (P < 0.05, VIP > 1) in the lactation group. Lactobionic acid, citric acid, orotic acid and oxamide were extracted by the S-plot as potential biomarkers of the metabolic difference between lactation and non-lactation. The TCA cycle, glyoxylate and dicarboxylate metabolism, glutamate metabolism and glycine metabolism were determined to be pathways that were significantly impacted (P < 0.01, impact value >0.1) in the lactation group. Among them, the TCA cycle was the most up-regulated pathway (P < 0.0001), with 7 of the 10 related metabolites increased in the MG tissues of the lactating cows.
The overall biofluid and MG tissue metabolic mechanisms in the lactating cows were interpreted in this study. Our findings are the first to provide an integrated insight and a better understanding of the metabolic mechanism of lactation, which is beneficial for developing regulated strategies to improve the metabolic status of lactating dairy cows.
Dairy cows contribute the most important source of dairy foods for humans, and this relies largely on the cows’ lactation performance . The inherent metabolic status during lactation and the difference from the non-lactation period are the most important metabolic biological processes for determining milk production and sustainability, and they depend on the development of the mammary gland (MG) and the overall coordinated metabolism and physiology.
Metabolomics, a vital part of systemic biology, can be utilized to explore the disease, diet or environment-related comprehensive metabolism by analyzing endogenous small molecules in the biofluids or tissues . To date, biofluid metabolomics combined with tissue metabolomics has provided valuable information on the overall or partial metabolic mechanisms . Biofluid metabolomics has been extensively applied in disease diagnosis , biomarker discovery  and novel pathway identification  in human clinical studies. Recently, biofluid metabolomics has been widely used to investigate the relationship between the rumen fluid metabolites and cow health , the effect of heat stress on milk metabolites , and the urinary biomarkers under different quality forage diets . However, most biofluid metabolomics studies in dairy cows have been limited to a single biofluid and/or focused on the changes in the physiological conditions under different treatments. Multiple biofluid metabolomics and their relationships under the same treatment allow new insights into the inherent metabolic pattern and global metabolism, which could contribute to improving the metabolic status for better productivity and sustainability. Tissue metabolomics can capture the subtle metabolic variation of the specific organs or parts of the body and reflect the actual biological processes that affect gene expression, transcription and translation . It has been revealed that MG metabolomics is an effective way to explore the initiation, maintenance and regulation of lactation [11, 12] and can provide direct connections with the milking phenotypes. Therefore, identifying the metabolites in the MG during lactation and comparing their key pathways with non-lactating cows could enhance our understanding of the lactation mechanism.
In this study, we hypothesize about the inherent metabolic characteristics during lactation and their differences from the non-lactation period, which can be addressed by biofluid metabolomics combined with tissue metabolomics. The relationships among four biofluids’ (rumen fluid, serum, milk and urine) metabolomics profiles from the lactating cows were analyzed to investigate the inherent metabolic patterns and possible effects on lactation sustainability. The MG metabolomics profiles were compared between lactating and non-lactating cows to identify the key metabolites and pathways and their potential regulatory roles in lactation.
The experimental procedures were approved by the Animal Care Committee at Zhejiang University (Hangzhou, China) and were in accordance with the university’s guidelines for animal research.
Study design and sample collection
In total, 12 Holstein dairy cows, 6 lactating and 6 non-lactating, were used in this study. Among them, 6 multiparous dairy cows with similar milk yields (30.4 ± 2.29 kg/d, mean ± SD) and at similar lactation stages (days in milk = 164 ± 19.6 d, mean ± SD) were fed a diet with a forage-to-concentrate ratio of 45:55, and they were provided 16.7% (DM basis) crude protein and 1.57 Mcal/kg net energy for lactation, following the protocol described in our previous study . The diet ingredients and nutrient composition are commonly applied in state-of-the-art farms throughout China (Additional file 1: Table S1). Approximately 1000 mg of MG tissues was collected immediately after slaughter according to previously described methods  from each cow in both the lactation and non-lactation groups. The tissues were then stored at −80 °C until metabolite extraction. Biofluid samples, including 50 mL of rumen fluid, 50 mL of milk, 10 mL of blood and 10 mL of urine, were collected from the lactating cows before the morning feeding at 6:00 AM. Rumen fluid was collected using an oral stomach tube following the standard procedures , milk samples were collected using a milk sampling device (Waikato Milking Systems NZ Ltd., Waikato, Hamilton, New Zealand), blood samples were collected from the jugular vein using pro-coagulation 10-mL tubes, and urine samples were collected using the vulval stimulation method. To minimize any possible degradation of the metabolites, rumen fluid, milk and urine samples were infused into a 15-mL sterilized centrifuge tube and immediately placed in liquid nitrogen. Before storage in a − 80 °C freezer, the rumen fluid, milk and urine samples were centrifuged at 6000×g, 3000×g, and 3000×g, respectively, and incubated at 4 °C for 15 min. Blood samples were centrifuged at 4 °C, 3000×g for 15 min within 20 min after sample collection. More details can be obtained from our previous study .
Extraction of compounds
The biofluid metabolomics was examined using an Agilent 7890 gas chromatography system equipped with a Pegasus 4D time of flight mass spectrometer (LECO, St. Joseph, MI, USA) . The tissue metabolomics procedures were performed as follows. First, 100 mg of MG tissue from each sample was added to a 2-mL Eppendorf tube with 0.4 mL of methanol-chloroform (Vmethanol: Vchloroform = 3:1) and 30 μL of L-2-chlorophenylalanine (1 mg/mL, stored in dH2O) and was mixed by vortexing for 10 s. Second, steel balls were placed in the tube and milled for 5 min at 55 Hz. The sample was then centrifuged at 4 °C at 12,000 rpm for 15 min. Third, approximately 0.4 mL of supernatant was transferred into a 2 mL silylated vial. An equal volume (10 μL) of each sample was placed in a new 2 mL silylated vial as a mixed sample for the quality control of the stability of the equipment system, the standard deviation of the beginning, middle and ending retention time of the mixed samples was less than 0.2, which indicates good stability.
The extracts were dried in a vacuum concentrator at 30 °C for 1.5 h. Methoxymethyl amine salt (80 μL; dissolved in pyridine to a final concentration of 20 mg/mL) was added to the dried metabolites, mixed and gently sealed. The solution was then incubated at 37 °C for 2 h in an oven. Then, 100 μL of bis trifluoroacetamide (containing 1% trimethylchlorosilane, v/v) was added to each sample, which was sealed again and incubated at 70 °C for 1 h. In addition, 10 μL of fatty acid methyl esters (a standard mixture of fatty acid methyl esters, C8-C16: 1 mg/mL; C18-C30: 0.5 mg/mL in chloroform) was added to the mixed sample and cooled to room temperature. Finally, the samples were mixed well and subjected to gas chromatography-time of flight/mass spectrometry (GC-TOF/MS) testing.
Metabolite identification by GC-TOF/MS
GC-TOF/MS analysis of the tissue samples was performed on an Agilent 7890 gas chromatograph system in cooperation with a Pegasus HT time-of-flight mass spectrometer (LECO, St. Joseph, MI, USA). The system used a DB-5 MS capillary column coated with 5% cross-linked diphenyl and 95% dimethyl polysiloxane (30 m × 250 μm inner diameter, 0.25 μm film thickness; J&W Scientific, Folsom, CA, USA). One μL aliquot of the analyte was added in splitless mode. Helium was used as the carrier gas. The flow of the front inlet purge was 3 mL min−1, and the gas flow rate through the column was 20 mL min−1. The original temperature was set at 50 °C and was maintained for 1 min. The temperature was increased to 330 °C at a speed of 10 °C min−1 and was maintained for 5 min at 330 °C. Temperatures of 280 °C, 280 °C, and 220 °C were used for the injection, transfer line, and ion source, respectively. The energy was −70 eV in electron influence mode. The full-scan mode of the mass spectrometry data was 85 m/z - 600 m/z at a rate of 20 spectra per second after a solvent delay of 366 s.
Identification of differentially expressed metabolites
Using the interquartile range denoising method, missing values of the raw data were filled by half of the minimum value and valid peaks were detected and only the metabolites remained. Additionally, an internal standard (L-2-chlorobenzene alanine) log normalization method was applied. The resulting three-dimensional data, including the peak number, sample name and normalized peak area, were uploaded to the online analysis platform - Metaboanalyst 3.0 (http://www.metaboanalyst.ca/) for univariate and multivariate analyses. Both one-way ANOVA and post hoc analysis were used to identify the important metabolites of the four biofluids. When the P value was less than 0.05, the metabolite was characterized as differentially expressed. Univariate analyses provided a preliminary overview of the features that were potentially significant. The false discovery rate (FDR) was used to conduct multiple comparisons testing, with an FDR value less than 0.05 indicating significance. Multivariate analyses included hierarchical cluster analysis (HCA), principal component analysis (PCA), partial least squares discriminant analysis (PLS-DA), and orthogonal projections to latent structures (OPLS). PCA was used to visualize the dataset of multitudinous variables in a 2 or 3-dimensional figure and display the general similarity and difference. The first principal component of the PLS (variable importance projection, VIP value) was obtained to identify the differentially expressed metabolites. Variables with VIP values exceeding 1.0 and P values less than 0.05 were selected as differentially expressed metabolites. An S-plot was generated to further identify the statistically significant and potentially biochemically significant metabolites in both their contributions to the module variables and the reliability of the module.
The compound names and relative concentrations of the shared metabolites were imported into Metaboanalyst 3.0 (http://www.metaboanalyst.ca) to perform functional enrichment and impact pathway analyses (integrates pathway enrichment analysis and pathway topology analysis). The functional enrichment analysis was based on several libraries containing approximately 6300 groups of biologically meaningful metabolite sets. Functional enrichment analysis was performed using the metabolite set enrichment analysis (MSEA) to investigate the enrichment of the predefined groups of functionally related metabolites, which are usually associated with biological pathways . Significant functional pathways were selected using a fold enrichment (FE) threshold greater than 2.0 and an FDR value less than 0.05 . Pathway topological analysis was used to calculate the importance of the metabolites based on the out-of-degree centrality and the relative betweenness measures of a metabolite in a given metabolic network . The impact pathway value was calculated as the sum of the matched important metabolites of all the metabolites in each pathway . Impact values greater than 0.10 and P values less than 0.05 were defined as significant impact pathways. Extremely different pathways were defined with P values less than 0.01. Additionally, commercial databases, including the Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.genome.jp/kegg/) and the Small Molecular Pathway Database (SMPDB, http://smpdb.ca), were used to search for metabolites and for the integrated pathway analysis.
The 4% fat-corrected milk yield of the six lactating cows was 26.29 ± 2.66 kg/d (mean ± SD), and the milk protein, milk fat and lactose contents were 3.29 ± 0.16%, 4.17 ± 0.37% and 4.91 ± 0.14% (mean ± SD), respectively. Detailed results can be found in our previous study . No adverse events occurred in any animal.
Metabolites identified in biofluids using GC-TOF/MS
From the metabolomics datasets of the four biofluids in six lactating cows, 33 shared metabolites were identified in the rumen fluid, serum, milk and urine, including noradrenaline, methylmalonic acid, glycine, lyxose, L-malic acid, thymol, creatine, 5-methoxytryptamine, oxoproline, glycerol, L-threose, m-cresol, aminomalonic acid, 2-hydroxybutanoic acid, alanine, hydroxylamine, phosphate, fumaric acid, glucose, prostaglandin E2, isoleucine, lactose, malonic acid, N-methyl-L-glutamic acid, phenylethylamine, 2,4-diaminobutyric acid, oxalic acid, lactic acid, 4-androsten-1-beta-ol-3,17-dione, norleucine, asparagine, conduritol-b-epoxide, and 5-aminovaleric acid. The relative abundance of each metabolite is presented in Additional file 1: Table S2.
When we compared the MG metabolomics and milk metabolomics of the 6 same lactating dairy cows in the lactation group, a total of 118 shared metabolites, 156 MG-specific metabolites and 67 milk-specific metabolites were identified (Fig. 1, Additional file 1: Table S3).
Metabolites identified in the MG tissue using GC-TOF/MS
In total, 394 valid peaks were detected from the analytes of the MG tissues. Based on the LECO/Fiehn Metabolomics Library, 274 metabolites were characterized and quantified (Additional file 1: Table S4). Among them, stigmasterol, sarcosine, palatinitol, oxamide, lyxose, L-threose, D-erythronolactone, cholesterol, carnitine and carbamoyl-aspartic acid were identified in the lactation group only and mainly belong to lipid or lipid-like molecules, carbohydrates and amino acids. Conversely, terephthalic acid, phenyl beta-D-glucopyranoside, octadecanol and alpha-ketoisocaproic acid were found only in the non-lactation dairy cow group.
Multivariate statistical analysis
The 3D–PCA score plot of the metabolic profiles showed significantly separated clusters among the four biofluids (Fig. 2a). All score plots for the rumen fluid, serum, milk, and urine sample were in the Hotelling T2 ellipse with 95% confidence. As shown in the PLS-DA score map, the samples were clearly separated into four parts, indicating the differential metabolic profiling of the four biofluids samples (Fig. 2b). HCA revealed sub-clusters containing varying numbers of metabolites, with different metabolomics profiles from each biofluid being clustered together (Fig. 2c). The metabolites from the rumen and serum samples were clustered together, which then clustered with the milk samples, but the cluster derived from the urine samples was separated from the other three biofluids. Several metabolites with higher relative concentrations were identified in each biofluid: glycerol, noradrenaline, lyxose, 5-methoxytryptamine and glucose in rumen fluid; lactic acid, glycine, N-methyl-L-glutamic acid, phenylethylamine, aminomalonic acid and glucose in the serum; glycerol, oxoproline, oxalic acid, L-malic acid, lyxose and L-threose in the milk; and isoleucine, m-cresol, 2-hydroxybutanoic acid, phosphate, creatine and lactose in the urine. HCA also revealed the presence of one cluster that included an amino acid and its derivatives (e.g., glycine, aminomalonic acid, N-methyl-L-glutamic acid, 5-methoxytryptamine, noradrenaline and phenylethylamine) and an organic acid and its derivatives (e.g., glucose, lactic acid and malonic acid). Another cluster consisted of aliphatic compounds, such as oxoproline, L-threose, oxalic acid, L-malic acid, glycerol and lyxose.
For the MG metabolomics, the heat map revealed the 6 lactating and 6 non-lactating dairy cows were assigned to 2 significantly separated clusters (Fig. 3a). The overview of the heatmap showed that the concentrations of most MG metabolites in the lactation group were higher than those in MG of non-lactating dairy cows. The 3-D PCA score map of GC-TOF/MS metabolic profiles of MG tissue showed significant discrimination between the two animal groups (Fig. 3b). All score plots for the MG samples in the 2 groups were in the 95% Hotelling T2 ellipse in the PLS-DA score map (Fig. 3c). In summary, clear separation in metabolomic profiles was found between the lactating and non-lactating dairy cow groups.
Identification of differentially expressed metabolites
The differentially expressed metabolites identified from the four biofluids are listed in Additional file 1: Table S5. Twenty-six metabolites showed obvious differences among the four biofluids, with four shared metabolites being differentially expressed in each comparison (FDR value <0.01): oxoproline, L-threose, creatine and lactose. As illustrated in Fig. 4, the VIP values of the PLS-DA model >1 included creatine, L-malic acid, oxalic acid, lyxose, lactose, phenylethylamine, glycerol, L-threose, m-cresol and oxoproline, with creatine exhibiting maximum integrated classification performance (VIP = 1.537).
In total, 73 differentially expressed metabolites were identified to in MG tissues between the lactating and non-lactating cows (VIP > 1, P < 0.05 & Q < 0.10), with 54 significantly higher and 19 significantly lower relative concentration metabolites in the lactating cows (Additional file 1: Table S4). The S-plot showed the key metabolites (Fig. 5). On the left-hand side of the S-plot, 4 variables with strong model contribution and high statistical reliability were explored as potential biomarkers to characterize the metabolic discriminations between lactation and non-lactation: lactobionic acid (Pcov = −4.669, Pcorr = −0.806), citric acid (Pcov = −2.983, Pcorr = −0.793), orotic acid (Pcov = −2.786, Pcorr = −0.887) and oxamide (Pcov = −2.491, Pcorr = −0.975).
KEGG pathway analysis
Overall, 20 pathways were obtained when the 33 shared metabolites from the four biofluids were imported into the KEGG analysis. Figure 6 shows the functional enrichment of the different pathways. The most enriched functional pathways included gluconeogenesis (FDR < 0.001, FE = 3.79), pyruvate metabolism (FDR = 0.003, FE = 2.68), the citric acid cycle (FDR = 0.007, FE = 2.63), glycerolipid metabolism (FDR = 0.017, FE = 2.01), and aspartate metabolism (FDR = 0.021, FE = 3.24). For the MG metabolomics, significantly higher abundance metabolites were involved in the 22 pathways in the lactation group (Additional file 1: Table S6). The key functional impact of the pathways is illustrated in Fig. 7. The results of the enrichment and impact pathways demonstrated that there were four significantly up-regulated pathways in the lactating cows: the tricarboxylic acid cycle (TCA), glyoxylate and dicarboxylate metabolism, glutamate metabolism, and glycine biosynthesis and degradation (P < 0.05, impact value >0.10). The integrated overview pathway map combined the above four key cellular pathways with the corresponding most significant and relevant metabolites (Fig. 8). Among these were seven metabolites that are involved in the TCA cycle, glyoxylate and dicarboxylate pathway, three metabolites that are involved in the glutamate metabolism pathway and three metabolites that are involved in the glycine biosynthesis and degradation pathway.
Remarkable progress has been made in our understanding of dairy cow metabolism, especially of the nutrient requirements of the biochemical reactions and functions . However, the complex metabolic mechanisms that occur during lactation are still considered a “black box” due to our lack of in-depth understanding . The increasing application of the -omics technologies (i.e., metagenomics, transcriptomics, proteomics and metabolomics) have become a powerful tool leading to the understanding of this “black box” . Indeed, these -omics technologies are being used to evaluate the physiological and molecular changes of dairy cows under different environments, such as varied diets and management. The inherent metabolic status during lactation and the difference in the non-lactation period play a vital role in the lactation physiology, which determines the milk production, udder health and dairy sustainability . Very little information is available on the metabolomics profiles of these two stages in the dairy cows. To our knowledge, this paper is the first to address the lactation-related metabolic profiles by using tissue and biofluid metabolomics. Both the overall lactation-related metabolism (variation and coordinated metabolic mechanism based on four biofluids’ metabolomics) and partial metabolism (comparison metabolomics in MGs between lactating and non-lactating cows) provide novel insights and a greater understanding of the metabolic mechanism of lactation.
For dairy cows, common biofluids include the rumen fluid, blood, milk and urine. The rumen fluid metabolites are critical for microbial metabolism and growth, host nutrient utilization and health, all of which can influence milk production . Blood is an intermediary biofluid that consists of metabolites that are excreted from the body and links different organs and tissues . The milk metabolomics reflect the product quality, metabolic processes and status of MG . Urine is rich in metabolites and contains many biochemical pathway signals . It is easy to distinguish the rumen fluid, serum, milk and urine based on their physical characteristics, such as the color and odor. However, the metabolic relationships, which may contain subtle metabolic profiling for the whole body, are not well understood. According to the metabolic flow in the dairy cow, the nutrients in the diet are first ingested and digested in the rumen, absorbed or passed through the rumen, transferred into blood in the small intestines, and then allocated to and utilized in the different organs, such as the MG. Waste materials are then extracted from the body through the urine . It is reasonable to consider that the metabolites in the rumen fluid and serum are similar, whereas those in the urine are the most different from the other biofluids. In this study, creatine was characterized as an important metabolite that crosses four biofluids that have potential associations with lactation performance and health in dairy cows. Creatine is an important intermediate metabolite in the energy reactions, and its phosphorylated form is an important metabolite that is part of the energy shuttle . In the reverse reaction, a high-energy phosphate group, such as ATP, is transferred to creatine to form phosphocreatine and ADP. There are two ways in which creatine may accumulate in animals due to the continual replacement of lost creatinine : diet and de novo synthesis. Creatine is mainly synthesized in the liver from three amino acids: arginine, glycine and methionine. It is generally used by tissues, absorbed into blood and excreted in urine . The relative concentration across the four biofluids from high to low were serum > urine > milk > rumen fluid, which represented the basic metabolic distribution of this metabolite. Milk creatine can help infants who suffer from defects in the arginine:glycine amidinotransferase and glycine amidinotransferase . It is reported that more creatine is needed to ease the serious negative energy balance in dairy cows . The creatine content is decreased in cows with ovarian inactivity, which may be related to the increased energy consumption in inactive cows . Additionally, creatine was identified as a potential biomarker in diagnosing the heat stress status in dairy cows, which may be attributed to the phosphocreatine in the muscle tissue that have been mobilized for energy supply .
MSEA is a group-based approach that does not require the pre-selection of metabolites by using an arbitrary threshold. The key idea behind MSEA is to directly investigate the enrichment of pre-defined groups of functionally related metabolites (or metabolites sets) instead of individual metabolites, which have proven to be successful in deriving new information from the untargeted metabolomics studies . In this study, all of the shared metabolites in the four biofluids were evaluated together, and their related biological information was incorporated into the results. For the mid-lactation dairy cows, the most active and important pathways, including gluconeogenesis, pyruvate metabolism, the TCA cycle, glycerolipid metabolism and aspartate metabolism, which were identified, reflected and confirmed by the overall metabolomics profiles of the four biofluids. Among them, gluconeogenesis, pyruvate metabolism and the TCA cycle are the key pathways of energy metabolism. Gluconeogenesis generates glucose from lactate, glycerol and glucogenic amino acids . Glucose is extremely important for milk synthesis. In dairy cows, gluconeogenesis mostly occurs in the liver and provides up to 90% of the glucose required for host maintenance and production . Pyruvate is the starting point of gluconeogenesis and the end product for glycolysis, and it can also be generated by the transamination of alanine . Almost all of the dietary carbohydrates are fermented to volatile fatty acids (acetate, propionate and butyrate) in the rumen of the dairy cows, with propionate as the predominant substrate for gluconeogenesis . Pyruvate is an important intermediate metabolite for the generation of propionate from the succinic pathway or the lactate pathway . It can be converted by the pyruvate dehydrogenase complex into acetyl-CoA, which can then enter the TCA cycle. The TCA cycle plays a central role in cellular respiration and the supply of energy to all living cells , which is of paramount significance to the cell’s metabolic efficiency and, therefore, to the cow’s metabolism and production . Insufficient energy is highly associated with ketosis (a severe metabolic disease) and can lead to decreased milk production . In the glycerolipid metabolism pathway, glycerol can be converted into glucose by the liver and provide energy for cellular metabolism and lactation maintenance . Aspartate is a precursor of many compounds that are involved in cellular signaling, such as N-acetyl-aspartate. It is also a metabolite in the urea cycle and participates in gluconeogenesis in dairy cows . Based on the overview map of the KEGG pathway, these 5 key metabolic pathways were closely integrated together, suggesting that the carbohydrate pathways and energy pathways play vital roles in regulating lactation maintenance for the dairy cow.
Compared to the increase in the human breast metabolomics studies, to our knowledge, only 1 study has reported the goat MG secretory tissue metabolomics using NMR and identified 46 metabolites, with lactose, glutamate, glycine and lactate being the most abundant . In our study, the most abundant metabolites in the lactating dairy cow MGs were lactic acid, lactobionic acid, oxoproline, alanine, diglycerol, glycine, citric acid, creatine and glutamic acid, of which lactobionic acid and oxoproline were first identified following other studies [27, 47]. The milk metabolomics of the dairy cows mainly focused on the composition and phenotype-related analyses, including distinguishing the formula milk from breast milk , comparing the colostrum and milk , discerning Holstein cows and other minor dairy animals , and analyzing the relationship with heat stress  or methane emission . In this study, we also compared the similarities and differences between MG and milk metabolomics from the same lactating dairy cows. However, the four potential biomarkers (citric acid, lactobionic acid, oxamide, and orotic acid), including those metabolites that were found specifically in the MG from lactating cows, were also identified in the milk, which indicates that milk better reflects the physiological status of the MG compared to the non-lactating cows.
In this study, the S-plot combined with the impact pathway analysis was proven to be an effective and easy way to screen most statistically significant metabolites and functional impact pathways. The S-plot indicated that lactobionic acid, citric acid, orotic acid and oxamide (identified only in the lactation group) play vital roles in discriminating between lactating and non-lactating cows. The TCA cycle, glutamate metabolism and glycine metabolism pathways are the most important and work together in the MG for lactation initiation. Milk component synthesis and secretion is the main transition in lactation initiation. Lactobionic acid is the intermediate metabolite in the lactose biosynthesis pathway, is produced by lactose oxidation and has strong mineral complex properties that serve as promising bioactive ingredients in human nutrition . Citric acid is formed in the TCA cycle or from the diet and participates in the intermediate metabolism of carbohydrate oxidation in animal tissues . The level of citric acid is significantly higher in the lactating dairy cows and may enhance energy by participating in the TCA cycle. Orotic acid is found at a high concentration in bovine milk and almost exclusively originates in the MG cells . Orotic acid has been applied to improve athletic performance and body composition, and it has proven to enhance ATP . Oxamide is a diamide that is derived from oxalic acid and participates in the TCA cycle. In this study, oxamide and other secondary metabolites were significantly higher in the lactating dairy cows, indicating that the TCA cycle is more active during lactation initiation. Oxamide was also identified in the milk in the lactation group, suggesting that the TCA cycle is up-regulated in the lactation group.
Moreover, the functional impact pathway analysis showed that the TCA cycle pathway had the highest P value between the groups, which was mainly attributed to the 7 up-regulated metabolites of the 10 main substrates in the TCA cycle. It has been reported that a critical energy output is required during the initiation of lactation to support milk synthesis and secretion of mammary epithelial cells . Figure 8 depicts the integration of the most significant pathways and the most significant and relevant metabolites, which suggests the different metabolic mechanisms in the MG of the lactating cows compared to the non-lactating cows. Therefore, the downstream metabolites in the glutamic acid pathway entered the TCA cycle to generate succinate in the mitochondria through succinate semialdehyde, thereby releasing NADH. It is well known that NADH is used by the oxidative phosphorylation pathway to generate ATP , which can transport chemical energy within the cells for metabolism and provide a larger amount of energy for lactation maintenance and related biological processes. In addition to supplying energy sources, the TCA cycle integrates many other pathways to unify carbohydrate, protein and fat metabolism . For example, citric acid can be transported out of the mitochondrion to produce cytosolic acetyl-CoA for fatty acid and cholesterol (which further synthesizes steroid hormones and vitamin D) synthesis, which is also very important for lactation maintenance . Another intermediate of glutamate metabolism, γ-glutamyl-cysteine, is used to synthesize glutathione, which is vital for preventing damage to important cellular components of the mammary epithelial cells . As a prime metabolic source of glutathione, creatine, purines and serine and a protein constituent in the lactating mammary epithelial cells, glycine plays a significant role in various biological processes . With respect to the Gly pathway, two metabolites were highlighted: sarcosine and glycocyamine. Sarcosine is an intermediate and byproduct in glycine synthesis and degradation . It can be rapidly degraded into glycine in the mitochondria, which serves as a substrate along with acetyl-CoA to generate L-2-aminoacetoacetate and interacts with the pyruvate metabolism pathway. Glycocyamine is a precursor of creatine, which serves as an essential substrate for muscle energy metabolism and has been used in the feed additive industry to improve production quality . The above results suggest that major energy-related metabolic changes occur in the MG of lactating cows at the initiation of lactation to accommodate the demand for increased milk synthesis nutrients.
Using MG tissue and multi-biofluid metabolomics, the metabolic mechanisms of inherent lactation metabolic patterns and differences with non-lactation periods were addressed in this study. Based on the HCA of the 33 shared metabolites, the rumen fluid and serum were highly correlated and grouped together with milk. Creatine was characterized as the key metabolite to explain the biological variation among the four biofluids. For the lactating cows, lactobionic acid, citric acid, orotic acid and oxamide were identified as potential biomarkers among 54 differentially expressed metabolites. Gluconeogenesis, pyruvate metabolism, the TCA cycle, glycerolipid metabolism and the aspartate metabolism pathways were the most functionally enriched pathways. Integrated analysis of the differentially expressed metabolites involved in the TCA cycle, glutamate metabolism and glycine biosynthesis and degradation pathways revealed the probable key metabolic mechanism in the MG during lactation. Overall, our results provide a better physiological understanding of the lactation metabolism of mid-lactation dairy cows, which can help elucidate the regulated metabolic strategies for the lactating dairy cows in the future. More importantly, the combined application of the multi-biofluids and the tissue metabolomics in this study provide new insights into addressing complex biological questions. Further studies are required to validate the potential biomarkers and pathways identified in this study.
False discovery rate
Gas chromatography-time of flight/mass spectrometry
Hierarchical cluster analysis
metabolite set enrichment analysis
Orthogonal projections to latent structures
Principal component analysis
Partial least squares discriminant analysis
Tricarboxylic acid cycle
Variable importance projection
Ferreira AM, Bislev SL, Bendixen E, Almeida AM. The mammary gland in domestic ruminants: a systems biology perspective. J Proteome. 2013;94:110–23.
Sejrsen K, Hvelplund T, Nielsen MO. Ruminant physiology: digestion, metabolism and impact of nutrition on gene expression, immunology and stress. The Netherlands: Wageningen Academic Pub; 2006.
Johnson CH, Ivanisevic J, Siuzdak G. Metabolomics: beyond biomarkers and towards mechanisms. Nat Rev Mol Cell Biol. 2016;17(7):451–9.
Mörén L, Bergenheim AT, Ghasimi S, Brännström T, Johansson M, Antti H. Metabolomic screening of tumor tissue and serum in glioma patients reveals diagnostic and prognostic information. Meta. 2015;5(3):502–20.
Ryan D, Newnham ED, Prenzler PD, Gibson PR. Metabolomics as a tool for diagnosis and monitoring in coeliac disease. Metabolomics. 2015;11(4):980–90.
Sun H, Wang B, Wang J, Liu H, Liu J. Biomarker and pathway analyses of urine metabolomics in dairy cows when corn stover replaces alfalfa hay. J Anim Sci Biotechnol. 2016;7(1):49.
Menni C, Graham D, Kastenmüller G, Alharbi NH, Alsanosi SM, McBride M, Mangino M, Titcombe P, Shin S-Y, Psatha M. Metabolomic identification of a novel pathway of blood pressure regulation involving hexadecanedioate. Hypertension. 2015;66(2):422–9.
Saleem F, Ametaj B, Bouatra S, Mandal R, Zebeli Q, Dunn S, Wishart D. A metabolomics approach to uncover the effects of grain diets on rumen health in dairy cows. J Dairy Sci. 2012;95(11):6606–23.
Tian H, Zheng N, Wang W, Cheng J, Li S, Zhang Y, Wang J. Integrated Metabolomics study of the milk of heat-stressed lactating dairy cows. Sci Rep. 2016;6:24208.
Huang Q, Tan Y, Yin P, Ye G, Gao P, Lu X, Wang H, Xu G. Metabolic characterization of hepatocellular carcinoma using nontargeted tissue metabolomics. Cancer Res. 2013;73(16):4992–5002.
Rawson P, Stockum C, Peng L, Manivannan B, Lehnert K, Ward HE, Berry SD, Davis SR, Snell RG, McLauchlan D. Metabolic proteomics of the liver and mammary gland during lactation. J Proteome. 2012;75(14):4429–35.
Li Z, Liu H, Jin X, Lo L, Liu J. Expression profiles of microRNAs from lactating and non-lactating bovine mammary glands and identification of miRNA related to lactation. BMC Genomics. 2012;13(1):1.
Wang B, Mao S, Yang H, Wu Y, Wang J, Li S, Shen Z, Liu J. Effects of alfalfa and cereal straw as a forage source on nutrient digestibility and lactation performance in lactating dairy cows. J Dairy Sci. 2014;97(12):7706–15.
Apelo SA, Singer L, Lin X, McGilliard M, St-Pierre N, Hanigan M. Isoleucine, leucine, methionine, and threonine effects on mammalian target of rapamycin signaling in mammary tissue. J Dairy Sci. 2014;97(2):1047–56.
Shen J, Song L, Sun H, Wang B, Chai Z, Chacher B, Liu J. Effects of corn and soybean meal types on rumen fermentation, nitrogen metabolism and productivity in dairy cows. Asian austral. J Anim Sci. 2015;28(3):351.
Sun HZ, Wang DM, Wang B, Wang JK, Liu HY, Guan le L, Liu JX. Metabolomics of four biofluids from dairy cows: potential biomarkers for milk production and quality. J Proteome Res. 2015;14(2):1287–98.
Xia J, Sinelnikov IV, Han B, Wishart DS. MetaboAnalyst 3.0—making metabolomics more meaningful. Nucleic Acid Res. 2015;43(W1):W251–W57.
Xia J, Wishart DS. MSEA: a web-based tool to identify biologically meaningful patterns in quantitative metabolomic data. Nucleic Acid Res. 2010;38(suppl 2):W71–7.
Abbasi A, Hossain L, Leydesdorff L. Betweenness centrality as a driver of preferential attachment in the evolution of research collaboration networks. J Inf Secur. 2012;6(3):403–12.
Xia J, Wishart DS. MetPA: a web-based metabolomics tool for pathway analysis and visualization. Bioinformatics. 2010;26(18):2342–4.
Wang D, Liang G, Wang B, Sun H, Liu J, Guan LL. Systematic microRNAome profiling reveals the roles of microRNAs in milk protein metabolism and quality: insights on low-quality forage utilization. Sci Rep. 2016;6:21194.
Drackley J, Donkin S, Reynolds C. Major advances in fundamental dairy cattle nutrition. J Dairy Sci. 2006;89(4):1324–36.
Garnsworthy PC: Nutrition and lactation in the dairy cow: Elsevier; 2013.
Loor JJ, Bionaz M, Invernizzi G. Systems biology and animal nutrition: insights from the dairy cow during growth and the lactation cycle. Syst biol. Livest Sci. 2011:215–46.
Mao SY, Huo WJ, Zhu WY. Microbiome–metabolome analysis reveals unhealthy alterations in the composition and metabolism of ruminal microbiota with increasing dietary grain in a goat model. Environ Microbiol. 2016;18(2):525–41.
Kenéz Á, Dänicke S, Rolle-Kampczyk U, von Bergen M, Huber K. A metabolomics approach to characterize phenotypes of metabolic transition from late pregnancy to early lactation in dairy cows. Metabolomics. 2016;12(11):165.
Yang Y, Zheng N, Zhao X, Zhang Y, Han R, Yang J, Zhao S, Li S, Guo T, Zang C. Metabolomic biomarkers identify differences in milk produced by Holstein cows and other minor dairy animals. J Proteome. 2016;136:174–82.
Wang X, Zhang A, Han Y, Wang P, Sun H, Song G, Dong T, Yuan Y, Yuan X, Zhang M. Urine metabolomics analysis for biomarker discovery and detection of jaundice syndrome in patients with liver disease. Mol Cell Proteomics. 2012;11(8):370–80.
Mercimek-Mahmutoglu S, Stoeckler-Ipsiroglu S, Adami A, Appleton R, Araújo HC, Duran M, Ensenauer R, Fernandez-Alvarez E, Garcia P, Grolik C. GAMT deficiency features, treatment, and outcome in an inborn error of creatine synthesis. Neurol. 2006;67(3):480–4.
Wyss M, Kaddurah-Daouk R. Creatine and creatinine metabolism. Physiol Rev. 2000;80(3):1107–213.
Poortmans JR, Kumps A, Duez P, Fofonka A, Carpentier A, Francaux M. Effect of oral creatine supplementation on urinary methylamine, formaldehyde, and formate. Med Sci Sport Exer. 2005;37(10):1717.
Schulze A. Creatine deficiency syndromes. Mol Cell Biochem. 2003;244(1-2):143–50.
Wang Y, Gao Y, Xia C, Zhang H, Qian W, Cao Y. Pathway analysis of plasma different metabolites for dairy cow ketosis. Ital J Anim Sci. 2016;15(3):545–51.
Xu C, Xia C, Sun Y, Xiao X, Wang G, Fan Z, Shu S, Zhang H, Xu C, Yang W. Metabolic profiles using 1H-NMR spectroscopy in postpartum dairy cows with ovarian inactivity. Theriogenology. 2016;86(6):1475–81.
Tian H, Wang W, Zheng N, Cheng J, Li S, Zhang Y, Wang J. Identification of diagnostic biomarkers and metabolic pathway shifts of heat-stressed lactating dairy cows. J Proteome. 2015;125:17–28.
Nelson DL, Lehninger AL, Cox MM: Lehninger principles of biochemistry: Macmillan; 2008.
Aschenbach JR, Kristensen NB, Donkin SS, Hammon HM, Penner GB. Gluconeogenesis in dairy cows: the secret of making sweet milk from sour dough. IUBMB Life. 2010;62(12):869–77.
Denton R, Halestrap A. Regulation of pyruvate metabolism in mammalian tissues. Essays Biochem. 1978;15:37–77.
Zhang Q, Koser SL, Bequette BJ, Donkin SS. Effect of propionate on mRNA expression of key genes for gluconeogenesis in liver of dairy cattle. J Dairy Sci. 2015;98(12):8698–709.
Jeyanathan J, Martin C, Morgavi D. The use of direct-fed microbials for mitigation of ruminant methane emissions: a review. Animal. 2014;8(02):250–61.
Grassian AR, Parker SJ, Davidson SM, Divakaruni AS, Green CR, Zhang X, Slocum KL, Pu M, Lin F, Vickers C. IDH1 mutations alter citric acid cycle metabolism and increase dependence on oxidative mitochondrial metabolism. Cancer Res. 2014;74(12):3317–31.
Antunes-Fernandes E, van Gastelen S, Dijkstra J, Hettinga K, Vervoort J. Milk metabolome relates enteric methane emission to milk synthesis and energy metabolism pathways. J Dairy Sci. 2016;99(8):6251–62.
Bravo D, Wall E. The rumen and beyond: nutritional physiology of the modern dairy cow1. J Dairy Sci. 2016;99(6):4939–40.
Doran AG, Berry DP, Creevey CJ. Whole genome association study identifies regions of the bovine genome and biological pathways involved in carcass trait performance in Holstein-Friesian cattle. BMC Genomics. 2014;15(1):1.
Piccioli-Cappelli F, Loor J, Seal C, Minuti A, Trevisi E. Effect of dietary starch level and high rumen-undegradable protein on endocrine-metabolic status, milk yield, and milk composition in dairy cows during early and late lactation. J Dairy Sci. 2014;97(12):7788–803.
Palma M, Hernández-Castellano LE, Castro N, Arguëllo A, Capote J, Matzapetakis M, de Almeida AM. NMR-metabolomics profiling of mammary gland secretory tissue and milk serum in two goat breeds with different levels of tolerance to seasonal weight loss. Mol BioSyst. 2016;12:2094–107.
Pisano MB, Scano P, Murgia A, Cosentino S, Caboni P. Metabolomics and microbiological profile of Italian mozzarella cheese produced with buffalo and cow milk. Food Chem. 2016;192:618–24.
Scano P, Murgia A, Demuru M, Consonni R, Caboni P. Metabolite profiles of formula milk compared to breast milk. Food Res Int. 2016;87:76–82.
Curtasu M, Theil P, Hedemann M. Metabolomic profiles of colostrum and milk from lactating sows. J Anim Sci. 2016;94(7supplement3):272–5.
Schaafsma G. Lactose and lactose derivatives as bioactive ingredients in human nutrition. Int Dairy J. 2008;18(5):458–65.
Krebs HA, Johnson W. The role of citric acid in intermediate metabolism in animal tissues. FEBS Lett. 1980;117(S1):2–10.
Saidi B, Warthesen J. Analysis and stability of orotic acid in milk. J Dairy Sci. 1989;72(11):2900–5.
Loeffler M, Carrey EA, Zameitat E. Orotic acid, more than just an intermediate of Pyrimidine de novo synthesis. J Genet Genomics. 2015;42(5):207–19.
Chen C-C, Stairs DB, Boxer RB, Belka GK, Horseman ND, Alvarez JV, Chodosh LA. Autocrine prolactin induced by the Pten–Akt pathway is required for lactation initiation and provides a direct link between the Akt and Stat5 pathways. Genes Dev. 2012;26(19):2154–68.
Malloy CR, Sherry AD, Jeffrey F. Evaluation of carbon flux and substrate selection through alternate pathways involving the citric acid cycle of the heart by 13C NMR spectroscopy. J Biol Chem. 1988;263(15):6964–71.
He W. Miao FJ-P, Lin DC-H, Schwandner RT, Wang Z, Gao J, Chen J-L, Tian H, Ling L: citric acid cycle intermediates as ligands for orphan G-protein-coupled receptors. Nature. 2004;429(6988):188–93.
Akers RM. Lactation and the mammary gland. Ames: Wiley; 2016.
Kannan N, Nguyen LV, Makarem M, Dong Y, Shih K, Eirew P, Raouf A, Emerman JT, Eaves CJ. Glutathione-dependent and-independent oxidative stress-control mechanisms distinguish normal human mammary epithelial cell subsets. P Natl Acad Scd USA. 2014;111(21):7789–94.
Shennan D, McNeillie S, Curran D. The effect of a hyposmotic shock on amino acid efflux from lactating rat mammary tissue: stimulation of taurine and glycine efflux via a pathway distinct from anion exchange and volume-activated anion channels. Exp Physiol. 1994;79(5):797–808.
Khan AP, Rajendiran TM, Bushra A, Asangani IA, Athanikar JN, Yocum AK, Mehra R, Siddiqui J, Palapattu G, Wei JT. The role of sarcosine metabolism in prostate cancer progression. Neoplasia. 2013;15(5):491–IN13.
Ostojic SM. Advanced physiological roles of guanidinoacetic acid. Eur J Nutr. 2015;54(8):1211–5.
The authors gratefully thank the personnel of Hangjiang Dairy Farm (Hangzhou, China) for their assistance in the feeding and care of the animals. We would like to thank B. Wang from the Institute of Dairy Science of Zhejiang University (Hangzhou, China) for assistance in animal trial and sampling.
This study was supported by the National Natural Science Foundation of China (No. 31472121) and the China Agriculture Research System (No. CARS-37). The funding body has not participated in or interfered with the research.
Availability of data and materials
Supplementary materials in a file accompany this paper. The raw.xps files of rumen, serum, milk, urine, liver tissue and MG tissue metabolomics have been submitted to MetaboLights (EMBI-EBI) under accession number MTBLS411.
The procedures of this study were approved by the Animal Care and Use Committee of Zhejiang University, Hangzhou, P. R. China and were in accordance with the university’s guidelines for animal research. All the animals used in this study were selected from Hangjiang Dairy Farm (Hangzhou, China).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Ingredients and nutrient composition of the experimental diet. Table S2. The relative abundance of mutual metabolites in rumen fluid, serum, milk and urine. Table S3. Mutual metabolites identified in the mammary gland and milk in lactation group. Table S4 The metabolites identified in the lactation (N) and non-lactation (NL) groups. Table S5. Differentially expressed metabolites among the 4 biofluids identified by one-way ANOVA and post-hoc analysis. Table S6. Significantly up-regulated pathways in lactation group. (DOCX 70 kb)