Metabolomics comparison of rumen fluid and milk in dairy cattle using proton nuclear magnetic resonance spectroscopy
Article information
Abstract
Objective
The metabolites that constitute the rumen fluid and milk in dairy cattle were analyzed using proton nuclear magnetic resonance (1H-NMR) spectroscopy and compared with the results obtain for other dairy cattle herds worldwide. The aim was to provide basic dataset for facilitating research on metabolites in rumen fluid and milk.
Methods
Six dairy cattle were used in this study. Rumen fluid was collected using a stomach tube, and milk was collected using a pipeline milking system. The metabolites were determined by 1H-NMR spectroscopy, and the obtained data were statistically analyzed by principal component analysis, partial least squares discriminant analysis, variable importance in projection scores, and metabolic pathway data using Metaboanalyst 4.0.
Results
The total numbers of metabolites in rumen fluid and milk were measured to be 186 and 184, and quantified as 72 and 109, respectively. Organic acid and carbohydrate metabolites exhibited the highest concentrations in rumen fluid and milk, respectively. Some metabolites that have been associated with metabolic diseases (acidosis and ketosis) in cows were identified in rumen fluid, and metabolites associated with ketosis, somatic cell production, and coagulation properties were identified in milk.
Conclusion
The metabolites measured in rumen fluid and milk could potentially be used to detect metabolic diseases and evaluate milk quality. The results could also be useful for metabolomic research on the biofluids of ruminants in Korea, while facilitating their metabolic research.
INTRODUCTION
Metabolomic studies have been actively conducted to not only investigate metabolic diseases found in animals, but also related human diseases, drug toxicity, and altered gene functions [1–3]. Metabolomics is the study of metabolic mechanisms that are achieved by measuring metabolic products in biological specimens as a whole and by elucidating the underlying chemical processes [2]. The most common techniques used for metabolomics are nuclear magnetic resonance (NMR) spectroscopy, liquid chromatography-mass spectrometry (LC/MS) and gas chromatography-MS (GC/MS). NMR spectroscopy is less sensitive than LC/MS and GC/MS, while providing at relatively simple sample preparation and highly reproducible molecule quantification [3,4]. In addition, NMR is a powerful analytical technique for the identification of metabolic biomarkers [5]. Therefore, NMR would be suitable for the analysis of numerous metabolites in the biofluids of ruminants.
Numerous researches have previously analyzed the metabolites in ruminant biofluids (e.g., rumen fluid, serum, milk, and urine) using NMR spectroscopy, and by the identification and quantification of metabolites through concentration comparison [1,6–8]. Moreover, metabolic pathways have been statistically analyzed by principal component analysis (PCA), partial least square discriminant analysis (PLS-DA), determination of variable importance in projection (VIP) scores, and hierarchical clustering analysis (with heatmaps) [1,6–8]. In particular, research is being conducted on finding new biomarkers for metabolic diseases found in ruminants (acidosis, ketosis, etc.) as well as for the early diagnosis of these diseases that may occur in the livestock industry; this has been achieved by identifying specific metabolites in ruminant biofluids [9,10]. Therefore, metabolites in ruminant biofluids must be studied to prevent and minimize the damages caused by metabolic diseases.
In Korea, NMR metabolomics studies have been conducted on samples such as natural food (potatoes and Astragalus root) products [11,12], human urine [13], liver, serum, urine of piglets [14], and obese mice [15]. Recently, metabolomic research has been conducted by comparing the composition of rumen fluid in Hanwoo’s cattle with those of the volatile fatty acids and monosaccharides in metabolites by using NMR spectroscopy, high performance liquid chromatography and high-performance anion-exchange chromatography [16]. A study has also investigated the effects of different roughage diets that were fed at concentrated ratios on the milk productivity and milk metabolites in dairy cows [17]. However, extremely few studies have been conducted on metabolites that constitute ruminant biofluids in Korea by using proton nuclear magnetic resonance (1H-NMR).
Therefore, this study aimed to measure the metabolites in dairy cow rumen fluid and milk by using 1H-NMR spectroscopy; these metabolites were quantified and classified to construct a database for each sample, wherein the concentration of each metabolite was provided. The role of metabolites in rumen fluid and milk observed in this study were compared to that of the metabolites in ruminant biofluids reported in a previous study. Moreover, the results of this study provide a useful a database for the analysis of metabolites in ruminant biofluids.
MATERIALS AND METHODS
All experimental protocols used in this study were approved by the National Institute of Animal Science Department of Animal Resources Development Dairy Science Division (Cheonan, Chungcheongnam-do, Korea; NIAS-201908).
Animals and sample collection
Six dairy cattle were used in this study, and all animals were fed the same diet that was composed of total mixed ration (TMR); the fed amounts were defined by their voluntary intake. The chemical composition of the TMR have been presented in Table 1. Rumen fluid samples were collected using stomach tube from dairy cows. All rumen fluid collected using conical tube (30 mL each). Subsequently, rumen fluid samples centrifuged at 806 g for 15 min to remove feed particles and the supernatant was stored at −80°C until analyzed for metabolites using 1H-NMR. Milk samples were collected by using pipeline milking system and then transferred to conical tubes (30 mL each). Samples were stored at −80°C until analyzing for metabolites by 1H-NMR.
NMR spectroscopy
The rumen fluid sample was recentrifuged at 12,902 g for 10 min and the supernatant was collected 300 μL. Standard buffer solution (2,2,3,3-d(4)-3-(trimethylsilyl)propionic acid [TSP] sodium salt) was added to 300 μL of supernatant in deuterium oxide (D2O) solvent/standard buffer solution (300 μL). The supernatants (600 μL) were transferred to 5 mm NMR tubes for NMR analysis [6].
The collected milk sample was centrifuged at 4,000 g for 15 min to remove the lipid layer in supernatant. Thereafter, the mixture of milk (250 μL) and D2O (300 μL) were transferred to 5 mm NMR tubes for NMR analysis [6].
1H-NMR spectra of rumen fluid and milk samples were obtained on a SPE-800 MHz NMR-MS Spectrometer (Bruker BioSpin AG, Billerica, MA, USA) at 64 K using a 5 mm triple-resonance inverse cryoprobe with Z-gradients (Bruker BioSpin CO., USA). The pulse sequence used for the rumen fluid and milk were a NOESY and Carr-Purcell-Meiboom-Gill pulse sequence collecting 64,000 data points with 128 transients, a spectral width of 16,025.641 Hz, a relaxation delay of 4.0 s, and an acquisition time of 2.0 s.
NMR metabolite measurement, quantification, and statistical analysis
The processed spectra were imported into the Chenomx NMR suite 8.4 software (Chenomx, Edmonton, Canada) for identification and quantification. The baseline and phase were matched for comparison between samples using the NMR processor. The following procedure was employed for qualitative and quantitative analysis of the metabolites in samples: the spectral width was 10 ppm and was referenced to the TSP signal at 0 ppm. The resources used were the Livestock Metabolite Database (http://www.lmdb.ca), Bovine Metabolite Database (http://www.bmdb.ca), and Chenomx library. Moreover, metabolite qualitative and quantitative were performed using the Chenomx profiler program.
Statistical analyses of the metabolite data were conducted using MetaboAnalyst version 4.0 (http://www.metaboanalyst.ca), an open source R-based program for metabolomics. During the statistical analysis the resulting metabolites were subjected to sample normalization, data transformation, and data scaling by using the “sum”, “log”, “pareto” functions, respectively. Univariate Student’s t-tests were used to identify differences between the metabolite profiles of the rumen fluid and milk samples. PCA and PLS-DA were used as multivariate data analysis techniques to generate a classification model and provide quantitative information for discriminating the metabolites. The different rumen fluid and milk metabolites were determined on the basis of a statistically significant threshold of VIP scores. Metabolites with VIP scores higher than 1.5 were obtained through PLS-DA. Metabolic pathways were quantified and, common metabolites in rumen fluid as well as milk metabolites of the other studied animals were statistically analyzed by MetaboAnalyst 4.0 for performing a metabolic pathway analysis, based on a database source by Kyoto encyclopedia of genes and genomes (http://www.kegg.com).
RESULTS
Measured quantification of metabolites by 1H-NMR spectroscopy
The results in shown Figure 1 and Supplementary Table S1, S2, and S3 reveal the measured and quantified compound in rumen fluid (A) and milk (B) that were measured by 1H-NMR. In the rumen fluid, 186 metabolites were measured and classified into 13 chemical classes. The classes with the most metabolites were others (34), carboxylic acids (27), and carbohydrates (27); the classes with the three highest concentrations were organic acids (22.759 mM), carbohydrates (4.871 mM) and lipids (0.961 mM). In addition, 72 metabolites were quantified (n≥4) in the rumen fluid. In milk, 184 metabolites were measured and classified into 14 chemical classes. The classes with the most metabolites were others (34), carboxylic acids (32), and carbohydrates (28); the classes with the three highest concentrations were carbohydrates (115 mM), carboxylic acids (9.703 mM), and lipids (6.270 mM). In addition, 109 metabolites were quantified (n≥4) in milk.

The classification of measured metabolites according to chemical class in rumen fluid (A) and, milk (B) by 1H-NMR analysis. Each square box color indicates the classification of metabolites, the numbers represents the measured metabolites, and the numbers in parentheses indicate the sum of the total concentrations of the measured metabolites. 1H-NMR, proton nuclear magnetic resonance.
1H-NMR spectroscopy analysis
Representative 1H-NMR spectra of 58 and 63 metabolites measured as single, doublet, and triplet peaks, as well as the reference substance TSP, in the rumen fluid and milk samples are shown in Supplementary Figures S1 and S2, respectively. Our data showed that the metabolites were different between rumen fluid and milk. To visualize the differences among the metabolites data, we performed PCA and PLS-DA (Figures 2, 3). Both score plots revealed differences in rumen fluid and milk, which were well separated in PC1 (39.8%) and PC2 (9.2%) for PCA and PC1 (39.8%) and PC2 (7.3%) PLS-DA. These results indicated significant variation among the different classes and concentrations of metabolites in the rumen fluid and milk. As shown in Figure 4, the two biofluids exhibited completely different metabolite profiles. VIP scores were also utilized to identify the metabolites that affected the differentiation in the PLS-DA score plot, and 21 metabolites were found to be significantly different (VIP score >1.5) in the rumen fluid and milk (Figure 4). In the rumen fluid, 14 metabolites (propionate, butyrate, valerate, acetate, caprate, isovalerate, glucose, isobutyrate, N-phenylacetylglycine, acetamide, 3-phenylpropionate, 4-methylhistidine, N-acetylglycine, and 2-aminobutyrate) had significantly higher concentrations than those present in milk. In milk, 7 metabolites (guanidinoacetate, glycine, lactose, glycolate, tartrate, galactose, and sn-glycero-3-phosphocholine) had significantly higher concentrations than those present in the rumen fluid.

Principal components analysis (PCA) score plot based on metabolites data in rumen fluid and milk by 1H-NMR analysis. On the score plot, each point represents an individual sample, with the green dot representing the rumen fluid group (n = 6), and the red dot representing the milk group (n = 6). The abscissa and ordinate represent the variance associated with PC 1 and 2, respectively. 1H-NMR, proton nuclear magnetic resonance.

Partial least square discriminant analysis (PLS-DA) score plot of rumen fluid and milk, by 1H-NMR analysis. The shaded ellipses represent the 95% confidence interval estimated from the score. On the score plot, each point represents an individual sample, with the green dot representing the rumen fluid group (n = 6), and the red dot representing the milk group (n = 6). The abscissa and ordinate represent the variance associated with components 1 and 2, respectively. 1H-NMR, proton nuclear magnetic resonance.

Variable importance in projection (VIP) scores of metabolites in rumen fluid and milk by 1H-NMR analysis. The selected metabolites were those with VIP score>1.5. Heat map with red or green boxes on the right indicates high and low abundance ratio, respectively, of the corresponding metabolite in rumen fluid and milk. VIP score was based on the PLS-DA model. 1H-NMR, proton nuclear magnetic resonance; PLS-DA, partial least square discriminant analysis; N-PAG, N-phenylacetylglycine; 3-PP, 3-phenylpropionate; 4-MH, 4-methylhistidine; N-agly, N-acetylglycine; AABA, 2-aminobutyrate; sn-glycero-3-P, sn-glycero-3-phosphocholine. VIP score value: propionate, 2.6698; butyrate, 2.6299; valerate, 2.3674; guanidoacetate, 2.2495; acetate, 2.1503; caprate, 2.1438; glycine, 2.1405; isovalerate, 2.1325; lactose, 2.0742; glucose, 1.9032; isobutyrate, 1.8975; N-PAG, 1.7465; glycolate, 1.7320; acetamide, 1.7051; tartrate, 1.6819; 3-PP, 1.6819; 4-MH, 1.6332; galactose, 1.6281; N-Agly, 1.5878; AABA, 1.5771; sn-glycero-3-P, 1.5696.
Top 30 average concentrations of rumen fluid and milk metabolites
The top 30 average concentrations of metabolites in the rumen fluid and milk are shown in Table 2 and 3. Among the metabolites measured in the rumen fluid, acetate, propionate, and butyrate (classified as organic acids) had the highest concentrations. In contrast, uracil (classified as a nucleoside and nucleotide), 3-hydroxyphenylacetate (classified as a carboxylic acid), and imidazole (classified as an imidazolinone) had the lowest concentrations. Among the metabolites measured in milk, lactose (classified as a carbohydrate), guanidinoacetate (classified as a carboxylic acid), and ethylene glycol (classified as a lipid) had the highest concentrations. In contrast, tartrate (classified as a benzoic acid), N-carbamoyl-β-alanine (classified as an organic acid), and acetone (classified as a others compound) had the lowest concentrations.

Average concentrations (mean±standard deviation) of top 30 metabolites in rumen fluid by 1H-NMR analysis (n≥4)
Common metabolites and metabolic pathways in the rumen fluid and milk
Fourth-three common metabolites were quantified (n≥4) in the rumen fluid and milk (Supplementary Table S1–3). An analysis showed that the common metabolites (Table 4, Figure 5) followed 16 metabolic pathways, among which the top five pathways were pyruvate metabolism, glycolysis/gluconeogenesis, glyoxylate and dicarboxylate metabolism, galactose metabolism, and glycerophospholipid metabolism. The metabolites in the relevant pathways were mainly carbohydrates and amino acids.

Metabolic pathway mapping of common quantified (n≥4) metabolites between rumen fluid and milk. The pathway impact analysis was performed using Metaboanalyst 4.0 software. The x-axis represents the pathway impact, and y-axis represents the pathway enrichment. The results are presented graphically as a bubble plot. The darker color and larger the size represent higher p-value from enrichment analysis and greater impact from the pathway topology analysis, respectively. Metabolic pathway name: 1, pyruvate metabolism; 2, glycolysis/glucoconeogenesis; 3, glyoxylate and dicarboxylate metabolism; 4, galactose metabolism; 5, glycerophospholipid metabolism; 6, ether lipid metabolism; 7, starch and sucrose metabolism; 8, tryptophan metabolism; 9, glycine, serine and threonine metabolism; 10, vitamin B6 metabolism; 11, arginine and proline metabolism; 12, amino sugar and nucleotide sugar metabolism; 13, histidine metabolism; 14, phenylalanine metabolism; 15, citrate cycle (tricarboxylic acid cycle); 16, propanoate metabolism.
DISCUSSION
Metabolites in ruminant biofluids can be affected by the species of ruminant, type of feed consumed, breeding environment and season, and the adopted analytical technique. Saleem et al [18] analyzed a total of 246 metabolites using 1H-NMR, GC/MS, LCMS/MS, and ICP-MS to study bovine rumination according to feeding concentrate ratio, 50 of which were measured by 1H-NMR. In this study, a total of 186 metabolites were measured by 1H-NMR, 72 of which were quantified in rumen fluid. O’Callaghan et al [6] studied bovine milk under different types of irradiation type by 1H-NMR and quantified a total of 49 metabolites. In this study, a total of 184 metabolites were measured by 1H-NMR, among which 109 were quantified in milk. It is possible that the higher NMR frequency and the Chenomx programs used in this study resulted in the quantification of more metabolites.
PCA is a multivariate method that can be used to classify and analyze large amounts of data with numerous variables and provide information on metabolite changes [19]. PLS-DA visualizes metabolic changes that are likely to affect classification for two or more individuals. VIP scores and heatmaps indicate the differences in the contributions of metabolites when the two groups are divided according to PLS-DA [20]. In Korea, multivariate analysis has been widely used for food origin determination [11,12] and human diseases [13]. Multivariate analysis is also an essential method for metabolite research; thus, studies on the metabolites found in ruminant biofluids are also be necessary. The PCA of the metabolites in rumen fluid and milk showed a separation between the groups. The PLS-DA was influenced by propionate, butyrate, valerate in rumen fluid and guanidinoacetate, glycine, and lactose metabolites in milk, which also demonstrated a separation between the groups.
Previous studies have shown that acetate, propionate, and butyrate in the rumen account for approximately 70% of ruminant energy and the majority of metabolite concentrations [18], which in accordance with this study. In terms of metabolic diseases, a higher concentrate ratio of the agricultural feed can cause acidosis [21]. According to Ametaj et al [22], a diet with a high concentrate ratio increases methylamine, glucose, alanine, maltose, propionate, uracil, valerate, xanthine, ethanol, and phenylacetate concentrations and decreases 3-phenylpropionate concentrations. In this study, the feeding concentrate ratio was measured in relation to the metabolites. As reported by Wang et al [23], the concentrations of metabolites based on biogenic amines (tyramine, putrescine, histamine, methylamine, and tryptamine) in the rumen fluid were observed to be higher in bovines suffering from acidosis. Sato and Shiogama [24] reported that changes in the concentrations of acetone and isopropanol in rumen fluid might be an indicator of ketosis. Compared with the metabolites detected in the above-mentioned previous studies [22–24], a majority of these were also detected in the rumen fluid in this study, including the diuretic-related biogenic amine metabolite histamine, methylamine, acetone, and isopropanol. Therefore, the metabolites measured in rumen fluid in this study can be used to verify the occurrence of diseases such as acidosis and ketosis in dairy cattle.
Numerous studies have been conducted on metabolites present in milk. Among them, a few have determined the chemical composition and nutritional value of milk along with the bioactive compounds present in milk; further, they have elucidated, the potential use of biomarkers in milk as a diagnostic tool using NMR spectroscopy [25]. According to Klein et al [9], the milk samples from cows exhibiting ketosis showed higher concentrations of ketone body metabolites (3-hydroxybutyrate [BHBA], acetoacetate, and acetone) and a higher glycerophosphocholine to phosphocholine ratio. Moreover, the presence of somatic cells was found to affect milk quality. Sundekilde et al [25] classified BHBA, lactate, lactose, hippurate, acetate, butyrate, and fumarate in milk as metabolites associated with somatic cell production, wherein the concentrations of hippurate and fumarate decreased with increasing number of somatic cells [26]. In addition, choline, carnitine, citrate, and lactose in milk are known to be associated with coagulation properties [27], and citirate, N-acetylcarbohydrates, trimethylamine, lecithin, and lactose in milk are known to be closely associated with the milk quality [28–30]. All metabolites associated with somatic cells were identified in this study, except for lactate and hippurate; further, the metabolites associated with coagulation properties were also detected. In addition, all metabolites associated with milk quality were detected except N-acetylcarbohydrates and lecithin, and all the ketone body metabolites were observed. Therefore, milk quality testing using 1H-NMR spectroscopy might be a viable tool for detecting the metabolic disease of ketosis in dairy cattle.
IMPLICATIONS
Proton nuclear magnetic resonance spectroscopy and statistical analyses were employed to analyze the metabolites in dairy cattle rumen fluid and milk. The metabolites measured in the rumen fluid and milk were mostly consistent with those reported in studies conducted abroad, and may be useful for predicting metabolic diseases, and milk quality. Furthermore, this report on metabolites in ruminant biofluids, which was achieved using the proton nuclear magnetic resonance analysis, will contribute to all the future ruminal metabolism studies in Korea.
Supplementary Information
Notes
CONFLICT OF INTEREST
We certify that there is no conflict of interest with any financial organization regarding the material discussed in the manuscript.
ACKNOWLEDGMENTS
This work was carried out with the support of “Cooperative Research Program for Agriculture Science and Technology Development (Project No. PJ01503902)” Rural Development Administration, Republic of Korea.