Introduction

The exclusive bamboo diet of giant pandas (Ailuropoda melanoleuca) is unique within the Order Carnivora and well documented in the literature. There are over 1,200 species of bamboo, of which giant pandas consume more than 601, 2. The growth of the bamboo follows a cyclic temporal pattern, which results in seasonal nutritional changes within the different plant parts3. Pandas primarily eat shoots when shoots are available which is only for a few months of the year during the shooting season3. During the shooting season bamboo invests its energy into shoot growth, not leaf development, but at the end of the shooting season bamboo leaves increase in nutrient value as the plants put growth and energy into photosynthesis4. When the shooting season is over, pandas eat the bamboo leaves. In the winter the leaves lose much of their available nutrient value while the bamboo goes dormant. During this time the culms have highest available nutrient composition (compared to leaves, and when no shoots are available), so pandas eat culms. Because of the cyclicity of the bamboo growth, the giant pandas consume different parts accordingly.

In the wild, more than one kind of bamboo species is available for pandas for the most time of a year. Due to such factors as season, transportation, weather, etc., captive giant pandas are generally provided with only one species of bamboo. Thus, pandas in captivity primarily eat shoots in spring, leaves in summer and fall, and culms in winter and early spring before the shoots grow in the next cycle5. Adult giant pandas are provided with several kinds of bamboo species each year6,7,8. At the Chengdu Research Base of Giant Panda Breeding (Chengdu, Sichuan Province, China), three main species fed to pandas at this facility include: Phyllostachys B., of which pandas only consume the culms, Bashania F., of which they only consume the leaves, and the bamboo shoots of Qiongzhuea O 9. Such single-plant part provisioning may continue for several months. Previous studies show that the nutritional value varies greatly in different parts of the bamboo body10,11,12,13. Moreover, the bioavailability of different plant parts varies significantly14, 15. So, if the consumption of single-parts continues for a long period of time, it may potentially lead to nutritional imbalances and negatively impact the health of giant pandas.

The objectives of this study were to determine the affect of consuming only certain parts of bamboo on the metabolism and health of pandas through the use of metabolomics. Nuclear Magnetic Resonance Spectrometry (NMR), Gas Chromatography-Mass Spectrometry (GC-MS) and Liquid Chromatography-Mass Spectrometry (LC-MS) are currently widely adopted methods for metabolomics analysis and detection. For NMR sample preparation, the sample pre-treatment is simple or even unnecessary under certain conditions, thus it is non-destructive and will not affect the structure or property of the samples16, 17, however, NMR has low sensitivity. In comparison, both types of MS have feature high efficiency, great sensitivity, qualitative and quantitative determination, and can detect and identify multiple metabolites in a complex mixture. Because of these advantages, MS is extensively used in metabolomic analysis, however the technique has several disadvantages such as matrix interference and difficulty in data comparison18. Using both GC-MS and NMR metabolomic analysis, we examined the metabolites of the blood serum extracted from captive giant pandas which were fed diets of different composition. The nutrient digestibility, metabolomic data and biochemical indicators were collected for these animals. All of these parameters may provide a biological explanation for the metabolism and general health of pandas in response to consuming different plant parts of bamboo.

Results

Daily food intake, nutrient digestibility and body mass gain

The daily dry matter intake differed significantly when the pandas consumed different bamboo parts ad libitum. The level of daily dry matter intake, crude fiber and energy intake were higher (p < 0.05) in “culm” group than in “shoot” and “leaf” groups, while the amount of ether extract intake and crude protein intake were the lowest in “culm” group (p < 0.05). For the “shoot” group, the level of daily dry matter intake was the lowest (p < 0.05) among the three groups, while the crude protein intake, digestibility and amount digested were the highest among the three groups (p < 0.05). The “shoot” group also had lower level of crude fiber intake and energy intake in comparison with the “culm” (p < 0.05)and “leaf” groups. Compared with “culm” and “leaf” groups, “shoot” group had higher body mass gain (p < 0.05) (Table 1).

Table 1 Daily nutrient intake (dry matter basis), nutrient apparent digestibility and body mass gain of pandas.
Full size table

Biochemical parameters

Clinical chemistry results show that BUN/CRE, BUN, TCHO, TG, VLDL, WBC and neutrophil levels were significantly different between these three groups (p < 0.05). Among these parameters, the BUN/CRE, BUN and WBC levels in the “shoot” group were much higher (p < 0.05) than in other two groups, and the Neutrophil% is higher in “shoot” group than in “leaf” (p < 0.05). And the TCHO, TG, VLDL levels in the “leaf” group were much higher (p < 0.05) than in the other two groups (Table 2).

Table 2 Conventional biochemical serum parameters from giant pandas kept on culm, leaf and shoot feeding regimes.
Full size table

NMR Results

1H NMR spectra of serum samples

Figure 1 shows the representative 1H NMR spectra of serum samples obtained from the three different groups and a total of 35 discrepant metabolites were shown in Table 3.

Figure 1
figure 1

Representative CPMG1H NMR spectra (δ0.7–5.4 and δ5.4-9.0) of serum obtained from “leaf” group, “culm” group and “shoot” group. The region of δ5.4–9.0 (in the dashed box) was magnified 16 times compared with corresponding region of δ0.7–5.4 for the purpose of clarity. Keys: 3-HB: 3-Hydroxybutyrate; 4HPPA, 4-Hydroxyphenylpyruvate; Ace: Acetate; Acet: Acetone; Aceta: Acetoacetate; Ala: Alanine; Asn: Asparagine; Cho: Choline; Cit: Citrate; Cr: Creatine; DMA: Dimethylamine; For: Formate; Fum: Fumarate; G: Glycerol; Glc: Glucose; Gln: Glutamine; Glu: Glutamate; Gly: Clycine; GPC: Glycerolphosphocholine; Ile: Isoleucine; L1: LDL, CH 3-(CH2)n-; L2: VLDL, CH 3-(CH2)n-; L3: LDL, CH3-(CH 2)n-; L4: VLDL, CH3-(CH 2)n-; L5: VLDL, -CH 2-CH2-C=O; L6: Lipid, -CH 2-CH=CH-; L7: Lipid, -CH 2-C=O; L8: Lipid,=CH-CH 2-CH=; L9: Lipid, -CH=CH-; Lac: Lactate; Leu: Leucine; Lys: Lysine; m-I: myo-Inositol; NAG: N-acetyl glycoprotein signals; OAG: O-acetyl glycoprotein signals; PC: Phophocholine; Phe: Phenylalanine; Pry: Pyruvate; Scy: Scyllitol; Tyr: Tyrosine; U1: Unknown1; Val: Valine.

Full size image
Table 3 OPLS-DA coefficients derived from the NMR data of metabolites in serum obtained from different groups.
Full size table

Multivariate data analysis of NMR data

A principal component analysis (PCA) was applied to the serum spectral data to investigate the correlations between the three groups. The PCA results showed a tendency to separate the three groups in the score plot (Fig. 2). PLS-DA was conducted on the serum spectra of the “culm”, “shoot” and “leaf” groups.

Figure 2
figure 2

PCA score plots based on the 1H NMR spectra of serum metabolites obtained from the “shoot”, “culm” and “leaf” group. (A) Total PCA score plots (R2X = 62.7%, Q2 = 0.479); (B) Shoot-Culm (R2X = 62.9%, Q2 = 0.426); (C) Shoot-Leaf (R2X = 80.6%, Q2 = 0.711); (D) Culm-Leaf (R2X = 67.6%, Q2 = 0.438). Two serum samples from “leaf” group were excluded because it positioned outside the Hotelling’s T2 ellipse on the PCA score plot.

Full size image

Furthermore, the metabolic changes were evident from the separation between the three bamboo forage groups plotted in the PLS-DA when expanding to the supervised multivariable statistical method, OPLS-DA (Fig. 3). The S-plots of OPLS-DA indicated that the metabolic profile of the groups deviated from each other, suggesting that significant biochemical changes were induced by the different bamboo diets. The corresponding coefficient analysis showed that the level of Valine, Leucine, Isoleucine, Lysine, Dimethylamine, Asparagine, Phenylalanine, Tyrosine, Creatine, Pyruvate, Lactate, Acetone, Acetoacetate, Scyllitol and Allantoin in the serum of shoot-fed pandas increased compared with the “culm” group (p < 0.05), while the level of Histidine, 4-hydroxyphenylpyruvate, Glutamine, Glycine, Alanine, Fumarate and lipids in the serum decreased (p < 0.05, Fig. 3 and Table 3). Compared with the “leaf” group, the level of Valine, Lysine, Asparagine, Phenylalanine, Tyrosine, Glutamine, Alanine, Creatine, Pyruvate, Citrate, Lactate, Acetoacetate, Scyllitol, O-acetyl-glycoprotein and Allantoin in the serum of the “shoot” group increased (p < 0.05), and most of the lipid compounds in the serum decreased (p < 0.05, Fig. 3 and Table 3). Compared with the “leaf” group, the level of Histidine, Phenylalanine, 4-hydroxyphenylpyruvate, Glutamine, Glycine, Alanine, Choline, Fumarate and Glycerol in the serum of the “culm” group increased (p < 0.05), and the level of Creatine, Acetone, N-acetyl-glycoprotein and most of the lipids compounds in the serum decreased (p < 0.05, Fig. 3 and Table 3).

Figure 3
figure 3

OPLS-DA scores plots (left panel) derived from 1H NMR spectra of serum and corresponding coefficient loading plots (right panel) obtained from different groups. The color map shows the significance of metabolites variations between the two classes. Peaks in the positive direction indicate metabolites that are more abundant in the groups in the positive direction of first principal component. Consequently, metabolites that are more abundant in the groups in the negative direction of first primary component are presented as peaks in the negative direction.

Full size image

GC-MS Results

The GC-MS total ion current chromatograms of the 28 serum samples from “culm”, “leaf” and “shoot” group are shown in Fig. 4. There are obvious changes in the three groups. Based on the LECO-Fiehn Rtx5 database and related references, most of the peaks were identified as metabolites, including the following: amino acids, glucose and fatty acids. These metabolites were mainly involved in glucose metabolism, lipid metabolism and amino acid metabolism.

Figure 4
figure 4

Representative GC-MS total ion chromatograms from the serum samples of “culm”, “leaf” and “shoot” group. Blue: Shoot group; Red: Leaf group; Black: Culm group.

Full size image

Multivariate statistical analysis of metabolites

Missing values of raw data were filled up by half of the minimum value, which means the peak area of all of samples. There are 107 missing values, and then, 438 peaks were detected and 186 metabolites could be left through an interquartile range denoising method. In addition, an internal standard normalization method used D-Phenylalanine, 2,4-dichloro-N-[(1,1-dimethylethoxy)carbonyl] as the interior label was employed in this data analysis. The resulting three-dimensional data set includes involving the peak number, sample name, and normalized peak area and was analyzed using the SIMCA software package (V14, Umetrics AB, Umea, Sweden) for principal component analysis (PCA) and the orthogonal projections to latent structures-discriminant analysis (OPLS). The PCA showed the distribution of the original unsupervised data structure (Fig. 5). In order to obtain a higher level of group separation and a better understanding of the variables responsible for group classification, a supervised partial least squares discriminant analysis (OPLS-DA)was applied. A 7-fold cross validation was used to estimate the optimal component number of the model and avoid the over-fitting. The R2 and Q2 intercept values were R2 = 0.827, Q2 = −0.365 (Culm vs Shoot); R2 = 0.860, Q2 = −0.227 (Leaf vs Shoot); R2 = 0.820, Q2 = −0.266 (Leaf vs Culm) after 200 permutations. The low values of the Q2 intercept indicate the robustness of the models, and thus show a low risk of over fitting and reliability. Based on the OPLS, a loading plot was constructed, which show the contribution of variables to the differences between each pair of groups (Fig. 6). The OPLS loadings also show the most important variables which were situated fararthest from the origin (Fig. 6). To refine this analysis, the first principal component of the variable importance projection (VIP) was obtained. The VIP values exceeding 1.0 were first identified as changed metabolites. In step 2, the remaining variables, that were different between the two groups, were evaluated for their significance to the model as assessed by Student’s T test (T-test), P > 0.05, variables were discarded between two comparison groups. In addition, commercial databases, including KEGG (http://www.genome.jp/kegg/) and NIST (http://www.nist.gov/index.html) were utilized to search for the MS library signals and associated gene expression of metabolites. Table 4 shows the identified differences metabolites in the three groups. Twenty-three metabolites with VIP-values greater than 1.0 and P values less than 0.05 were to be significantly different in the different groups. Compared with the “shoot”group, the level of Beta-alanine, Glycine, Alanine, Aminomalonic acid, Glucose, Aconitic acid, 1,5-anhydroglucitol, Threonic acid, Dl-anabasine, Threitol and Arbutin of the culm-fed pandas increased, and the level of Tryptophan, Lysine, Phenylalanine, Fructose, Mannose, Pyruvate, Lactate, D-glyceric acid and D-(glycerol 1-phosphate) decreased. Compared with the “shoot” group, the level of L-cysteine, 1,5-anhydroglucitol, Threonic acid, Dl-anabasine, and Threitol of the “leaf” group increased, and the level of Alanine and Lactate decreased. Compared with the “culm” group, the level of L-cysteine, Mannose, D-glyceric acid and Elaidic acid of the “leaf” group increased, and the level of Serine, Glycine, Alanine, Aminomalonic acid, Glucose, Aconitic acid, Lactate, and Arbutin decreased.

Figure 5
figure 5

Score plot of PCA model of TOTAL samples.

Full size image
Figure 6
figure 6

Score plot of OPLS model obtained from Culm vs Shoot (A); Leaf vs Shoot (C); Culm vs Leaf (E).Permutation analysis of OPLS-DA models from Culm vs Shoot (B); Leaf vs Shoot (D); Culm vs Leaf (F). Two hundred permutations were performed, and the resulting R2 and Q2 values were plotted. (Green triangle): R2; (blue square): Q2. The green line represents the regression line for R2 and the blue line for Q2.

Full size image
Table 4 The fold change in metabolite abundance among the Culm, Leaf and Shoot groupsa.
Full size table

Discussion

Bamboo is the staple food of giant pandas, however, there are various parts of bamboo for giant pandas to consume. Previous studies on captive giant panda nutrition have not explored the effect of different single- plant bamboo parts on the metabolism and health of pandas19. The result of the digestibility trial shows that the CP digestibility is higher in “shoot” group than “leaf” group (p < 0.05), and “leaf” group is higher than “culm” group (p < 0.05). Interestingly, the CF composition is much higher in “culm” group than in “leaf” group (p < 0.05), and “leaf” group is higher than in “shoot” group (p < 0.05) (Table 5). Based on these results, we speculated that high level of fiber will reduce the efficiency of the CP digestibility, which is similar in other monogastric mammals20, 21. Therefore, although the amount of daily CP intake in “shoot” group was similar to that of the “leaf” group, the digestible CP in “shoot” group was much higher than “leaf” group (p < 0.05).

Table 5 The nutrient composition of the culm in Phyllostachys B., leaf in Bashania F. and shoot in Qiongzhuea O.
Full size table

The daily EE intake from leaves was the highest, while the digestibility of EE from leaves was the lowest. In contrast, the daily EE intake from culms was the lowest, while the digestibility of EE from culms was the highest. This may be related to the low enzymatic activity of the lipase in the intestinal tract of giant pandas22. The EE in the leaves was too high for giant pandas to fully digest, which might result in the lower digestibility of EE from the leaves than that from the culms. These results show that the nutrients of bamboo will affect their gastrointestinal digestion and absorption capacity.

Maintenance of a stable body mass is often used as an indicator of health. Our data showed there was a significant difference in mass gain that correlated with experimental diet groups. During our study “shoot” group had higher body mass gain than “culm” and “leaf” groups. The “culm” pandas digested and absorbed the most calories and fiber compared to pandas in the other groups, and leaf group pandas digested and absorbed more calories than “shoot” pandas, however, “culm” and “leaf” pandas digested and absorbed less protein than shoot pandas. Therefore, we speculated that fiber digestion demands energy, and demands more energy to consume than it gives back in caloric value. So it takes physical effort, and energy to eat fiber, but no calories are returned and is thus negative in net caloric value. To explore the accuracy of this assumption, metabonomic methods were used to analyze the serum metabolites of the different treatment-fed pandas.

Creatine is involved in major energy metabolism in vertebrates23, 24, and contains more energy than ATP. High levels of 1,5-anhydroglucitol content in serum is an indicator of hypoglycemia25, 26. “culm” and “leaf” groups had lower creatine contents and higher 1,5-anhydroglucitol levels in serum than the “shoot” group, which shows that the energy reserve in the “culm” and “leaf” groups is lower than the “shoot” group pandas. Moreover, the “leaf” group had lower tricarboxylic acid cycle activity than “shoot” group, with lower citrate concentrations, because pyruvate is an important substrate of tricarboxylic acid cycle and citrate is involved in tricarboxylic acid cycle27. The “culm” group had a more compromised energy source and tricarboxylic acid cycle activity compared to the “shoot” group with lower concentrations of fructose, pyruvate, mannose, and creatine, and higher levels of aconitic acid, fumarate, 1,5-anhydroglucitol. The “culm” group had higher levels of alanine and aminomalonic acid than the “shoot” group most likely because alanine can react with oxomalonate to create the bi-products aminomalonic acid and pyruvate28. In addition, the fat intake was lower in the “culm” group than in the “shoot” group, however, the serum lipid level was higher in the “culm” animals than in the “shoot” group. These data show that “culm” group pandas were in short supply of energy with depressed tricarboxylic acid cycle activity, and mobilized amino acids, glycogen and lipids to compensatively trigger the tricarboxylic acid cycle to meet their energy requirements29,30,31. Thus, culm pandas were in trend of negative energy balance. Usually, negative energy balance is harmful to health32, 33. Therefore, we conclude that fiber digestion demands more energy to consume than it gives back in caloric value, which lead to lower growth rate of culm- and leaf-fed pandas.

Some studies show that ruminant and monogastric herbivores can use the energy produced through fiber fermentation and promote the body mass growth34, 35. However, the monogastric non-herbivores are different from the ruminant and monogastric herbivores, and a high fiber diet is not beneficial for energy reserves and body mass growth36, 37. Therefore, even though bamboo is the staple food of giant pandas, the bioavailability of fiber is relatively low, and high level fiber might reduce the efficiency of the CP digestibility. These results are in accordance with the conclusion of the previous studies which show that giant pandas’ gut microbiota is distinct from those of red pandas and other herbivores, and the coevolution of the giant panda and its gut microbiota are aberrant38,39,40, even though it has been 4.2 million years since they evolved from carnivory to herbivory41. Moreover, compared with other groups, shoot ingestion improved the most amino acids concentrations in the blood serum of “shoot” group. Combined with higher average body mass gain in “shoot” group, we also found that more digested and absorbed protein in bamboo enhanced body growth. Overall, because of the digestion, absorption and metabolism characteristics of pandas, different bamboo parts can affect the maintenance of body mass for pandas, thus impacting health.

The “leaf” group pandas had the highest triglyceride and total cholesterol concentrations compared to the other groups, and the “culm” group pandas were in negative energy balance. The “shoot” group had a greater increase in mass, however, the “shoot” group had higher levels of ketone bodies, lactate, urea, BUN, BUN/CRE than “culm” and “leaf” fed pandas, and “shoot” group pandas had higher WBC and neutrophils (Table 2) implying a potential inflammation reaction in the body. Because “shoot” fed pandas ingested more protein and amino acids than the other experimental groups, they had higher BUN and BUN/CRE, and metabolic burden aggravates the production of ketone bodies and lactate42,43,44. If pandas consume “shoot” for a long time, they may potentially suffer from nephropathy or other diseases. Because the high levels of BUN and BUN/CRE are the indicators of low glomerular filtration rate and kidney damage. Continuous high level of BUN tends to cause kidney disease. In addition, according to our husbandry observations, the number of female pandas in oestrous will decrease when culms are lacking in the diet (unpublished data). This phenomenon might be the result of lower fiber intake, because culms contain higher fiber than the other bamboo parts tested. In summary, different bamboo parts can affect the health of pandas.

The results of this study are meaningful not only for the management of ex-situ giant pandas but also for the conservation of the in-situ population as well. Bamboo species consumed by giant pandas vary with different habitat and elevations. In the Qinling Mountains nine bamboo species belonging to five genera have been described, however giant pandas prefer to consume Bashania fargesii and Fargesia spathacea 4, 45. These two species are the dominant native bamboos in the region. Schaller (1985) observed that giant pandas in the Wolong Nature Reserve fed mainly on Sinarundinaria fangiana and Fargesia spathacea, while Sinarundinaria nitida was also consumed, but to a much lesser extent3. Older studies reported that giant pandas in the Min Mountains consumed a single bamboo species, whereas animals in the Liang Mountains forage upon five bamboo species46. These data show that the nutrition intake of the wild pandas in the different reserves might vary significantly, this may affect wild giant panda health and reproductive abilities. As a result, it is necessary to estimate whether bamboo species found in the various nature reserves could provide a balanced nutrition for wild giant pandas, and in addition aid conservation managers in selecting the appropriate habitat for giant panda reintroductions.

Conclusion

The “culm” group absorbed the highest amount of calories and fiber, but was in short energy supply with a depressed tricarboxylic acid cycle activity, therefore, digestion of fiber requires energy input and yields low caloric extraction from the culm and leaf. In addition, the high level fiber will reduce the CP digestibility for giant panda. So, the high level of fiber not only make pandas in trend of negative energy balance, but also disturb the digestion of CP, these characteristic is similar with monogastric non-herbivores. The “culm” and “leaf” groups absorbed less amount of protein, and had a lower body mass increase than the “shoot” group, this showed that protein was benefit for body mass gain. However, long-term intake of shoots may have a potentially negative effect on the health. Therefore, the balance of an annual plant parts diverse bamboo diet is important for the health of giant pandas. In addition, we need to protect and conserve the biological diversity of bamboo in giant panda habitat, which is crucial for the survival and health of the in-situ population.

Methods

Study subjects and husbandry methods

A total of 17 adult captive giant pandas from the Chengdu Research Base of Giant Panda Breeding were the subjects of this study. There were three food groups based on bamboo species (Table 6). None of the female pandas were in estrus, pregnant, or nursing during the study. The bamboo species provided were Phyllostachys bissetii (culm group) from which the pandas only consume the culm, Bashania fargesii (leaf group) in which the pandas only consume the leaves and bamboo shoots of Qiongzhuea opienenss in which the pandas only consume the shoot (shoot group). Following the normal feeding protocol of the Chengdu Research Base of Giant Panda Breeding, all of the pandas were fed bamboo ad libitum, with experimental treatments lasting for 20 consecutive days. The supplemented diets were provided at the same mass every day among the three groups. The daily feed intake was recorded. The giant pandas were weighed at the beginning of each dietary treatment and every 15 days after consecutively consuming one specific bamboo species.

Table 6 Information of the adult captive giant pandas.
Full size table

Blood collection

Twenty-eight blood samples in total were collected from three groups after single bamboo plant parts intake for 20 consecutive days: nine blood samples from the “shoot” group, ten from the “culm” group and nine from the “leaf” group. Blood samples were collected from the giant pandas during routine physical examinations by experienced staff veterinarians at the Chengdu Research Base of Giant Panda Breeding. All sample collection and utility protocols in this study were approved by Chengdu Research Base of Giant Panda Breeding. The experimental procedures were fully in compliance with the current laws on animal welfare and research in China.

Blood chemistry measurements

Chemistry analyses were performed with an automatic biochemical analyzer (OLYMPUS AU-2700), which included alkaline phosphatase (ALP), blood urea nitrogen (BUN), creatinine (CRE), total cholesterol (TCHO), triglyceride (TG), very low-density lipoprotein (VLDL), white blood cell (WBC), LYMPH, Neutrophil and Monocytes.

Apparent digestibility trial

The apparent digestibility trial was conducted for consecutive 3 days after one-plant bamboo part intake for 20 days. All of the food provided, orts and feces were weighed. To calculate the amount of bamboo intake, three control groups were made to test the water loss under the same environmental conditions as the study animals were kept in. Samples of bamboo which the pandas ate and feces were collected twice per day, weighed and stored for nutrient composition analysis. The items of the nutrient composition include dry matter (DM), crude protein (CP), ether extract (EE), crude fiber (CF) and energy.

Dry Matter was measured using an oven drying method (DHG-9070A), crude protein was measured by kjeldahl method (K1100, Automatic kjeldahl apparatus), ether extract was measured by Soxhlet extraction method (SOX500), crude fiber was measured by Acid base digestion method (ANKOM A200i), energy concentration in bamboo and feces was measured by oxygen bomb calorimeter calorimetric method (PARR 6400 Calorimeter).

The apparent digestibility macro-nutrients were calculated as follows:

$${\rm{Apparent}}\,{\rm{digestibility}} \% =\frac{{\rm{nutrient}}\,\mathrm{digested}\,(\mathrm{nutrient}\,\mathrm{consumed}\,-\,{\rm{nutrient}}\,{\rm{excreted}})}{{\rm{nutrient}}\,\mathrm{consumed}\,(\mathrm{nutrient}\,\mathrm{offered}\,-\,{\rm{nutrient}}\,{\rm{refused}})}$$

Sample preparation of 1H-NMR metabolomics of serum

serum samples were removed from −80 °C storage and thawed at room temperature. 200 μL of serum and 400 μL of phosphate buffer compounded by D2O (Deuterium Oxide) were pipetted individually and mixed together for field-frequency lock. After centrifugation at 4 °C and 12,000 rpm for 10 min, 560 μL of the supernatants were transferred into NMR tubes (Norell, ST500-7) for 1H-NMR analysis47.

1H-NMR spectra of the serum samples were acquired at 298 K on a 600 MHz high resolution spectrometer (Agilent 600DD2, USA), operating at a 1H frequency of 600.11 MHz with three high resolution resonance probes. A total of 64 transients were collected with a spectral width of 20 ppm. The 90° pulse length was adjusted to 10 μs and the recycle delay is 2 s. A water-presaturated Carr-Purcell-Meiboom-Gill (CPMG) pulse sequence (RD -90sonτ-180-Pτ)n-ACQ) was used to record the NMR spectrum. A total spin-spin relaxation delay(2 nτ) of 200 ms was employed.

Sample preparation of GC-MS metabolomics

Take 100 μL sample into 1.5 mL EP tubes. Add 0.4 mL methanol and 20 μL of L-2-Chlorophenylalanine (0.1 mg/mL stock in dH2O) as an internal standard and vortex for 10 s. Centrifuge for 15 min at 12000 rpm, 4 °C. Transfer the supernatant (0.39 mL) into a fresh 2 ml GC/MS glass vial. Take 13 μL from each sample and pool as QC sample. Dry in a vacuum concentrator without heating. Add 60 μL of methoxyamination reagent (20 mg/mL. in pyridine) and incubate for 20 min at 80 °C. Add 80 μL BSTFA regent (1% TMCS, v/v) to the sample aliquots and incubate for 1 h at 70 °C. Add 10 μL FAMEs (Standard mixture of fatty acid methyl esters, C8-C16:1 mg/mL; C18-C24:0.5 mg/mL in chloroform) to the QC sample after it cools to room temperature. Mix well for GC-MS analysis.

GC-MS analysis. GC-MS analysis was performed using an Agilent 7890 gas chromatograph system coupled with a Pegasus 4D time-of-flight mass spectrometer. The system utilized a DB-5MS capillary column coated with 5% diphenyl cross-linked with 95% dimethylpolysiloxane (30 m × 250 μm inner diameter, 0.25 μm film thickness; J&W Scientific, Folsom, CA, USA). A 1 μL aliquot of the analyte was injected in splitless mode. Helium was used as the carrier gas, the front inlet purge flow was 3 mL min−1, and the gas flow rate through the column was 20 mL min−1. The initial temperature was kept at 80 °C for 1 min, then raised to 295 °C at a rate of 10 °C min−1, then kept for 7.5 min at 295 °C. The injection, transfer line, and ion source temperatures were 280, 280, and 220 °C, respectively. The energy was −70 eV in electron impact mode. The mass spectrometry data were acquired in full-scan mode with the m/z range of 85–600 at a rate of 20 spectra per second after a solvent delay of 465 s.

Data analysis

The mean daily food intake, nutrient digestibility, daily body mass gain, and conventional serum biochemical parameters were determined statistically relevant by one-way ANOVA analysis with SPSS 17.0 software (SPSS Inc., Chicago, IL). The distribution and variance of the data were tested before ANOVA analysis to ensure the parametric tests could be used. Datasets were performed by using post-hoc tests (LSD) for multiple comparisons to determine the statistical differences between groups. Data is presented as mean ± SD. The level of significance used was p < 0.05.

Data of 1H-NMR: The 1H-NMR spectra were processed using MestReNova software (V7.0, Mestrelab Research, Santiago de Compostella, Spain). All spectra were manually corrected for Fourier transformation, phase adjustment, baseline distortions and calibration. The signal-noise ratio was increased by multiplying an exponential function with a line-broadening factor of 1 Hz before Fourier transformation. The spectra of serum were referenced to the peak of the lactate at a chemical shift of δ1.33 ppm. The serum spectra integrating range and the signal integral computed in 0.002 ppm intervals across the region d 0.7–9.5 ppm. The discarded regions include 5.5–6.0 ppm for urea and 4.450–5.214 ppm for H2O.

Multivariate data analysis was performed on the normalized NMR datasets with the software package SIMCA-P+ (V11.0, Umetrics AB, Umea, Sweden). Principal component analysis (PCA) was conducted by using mean center scaling on the dataset to overview the intrinsic similarity/dissimilarity within the dataset.

Projection to latent structure–discriminant analysis (PLS–DA) and orthogonal projection to latent structure–discriminant analysis (OPLS–DA) were further conducted to obtain the metabolites with significant contributions to intergroup differentiation by using unit-variance scaled NMR data as the X-matrix and class information as the Y-matrix48. The quality of the model was assessed by the parameters R2X and Q2, which represent the total explained variations for the X matrix and the model predictability, respectively. The models were certified by using a six-segment cross validation49, 50. The metabolites associated with the group separations were indicated by coefficient coded loading plots calculated by back transformation of the loadings. P values were calculated by the correlation coefficient, and the correlation coefficient was calculated by the SIMCA software package (V14, Umetrics AB, Umea, Sweden). The correlation coefficient of│r│ > 0.632 was used as the cutoff value for the statistical significance based on the discrimination significance at the level of P = 0.05 and df (degree of freedom) = 8.

Data of GC-MS: Chroma TOF4.3X software of LECO Corporation and LECO-Fiehn Rtx5 database were used for raw peaks selection, the data baselines filtering and calibration of the baseline, peak alignment, deconvolution analysis, peak identification and integration of the peak area51. The RI (retention time index) method was used in the peak identification, and the RI tolerance was 5000. GC-MS spectral data were imported to SIMCA-P+ for principal components analysis (PCA), and orthogonal partial least-squares discriminant analysis (OPLS-DA).Variable Influence on Projection (VIP) values, which correspond to the importance of the variables were calculated, and variables with a VIP value larger than 1.0 were considered significant and used for further analysis and identification of the responsible peak(s) within the spectrum.