Plasma metabolomics reveals major changes in carbohydrate, lipid, and protein metabolism of abruptly weaned beef calves

Scientific Reports volume 13, Article number: 8176 (2023) Cite this article

424 Accesses

1 Altmetric

Metrics details

1H NMR-based metabolomics was used to study the effect of abrupt weaning on the blood metabolome of beef calves. Twenty Angus calves (258 ± 5 kg BW; 5 to 6 months old) were randomly assigned to a non-weaned (NW) group that remained grazing with their dam or a weaned (W) group that underwent abrupt separation from their dam to a separate paddock on d 0 of the study. Body weight, behaviour, and blood samples for cortisol and metabolomics were measured at d 0, 1, 2, 7, and 14 of the study. On d 1 and 2, W calves spent less time grazing and ruminating, and more time vocalising and walking, had a greater concentration of cortisol, NEFA, 3-hydroxybutyrate, betaine, creatine, and phenylalanine, and lesser abundance of tyrosine (P < 0.05) compared to NW calves. Compared to NW calves at d 14, W calves had greater (P < 0.01) relative abundance of acetate, glucose, allantoin, creatinine, creatine, creatine phosphate, glutamate, 3-hydroxybutyrate, 3-hydroxyisobutyrate, and seven AA (alanine, glutamate, leucine, lysine, phenylalanine, threonine and valine) but lesser (P < 0.05) relative abundance of low density and very low-density lipids, and unsaturated lipids. Both PCA and OPLS-DA showed no clustering or discrimination between groups at d 0 and increasing divergence to d 14. Blood metabolomics is a useful tool to quantify the acute effects of stress in calves during the first 2 days after abrupt weaning, and longer-term changes in carbohydrate, lipid and protein metabolism due to nutritional changes from cessation of milk intake and greater reliance on forage intake.

Nutritional requirements of new-born and pre-ruminant beef calves are met exclusively from protein, fat, carbohydrates (mainly lactose), vitamins and minerals absorbed from milk that bypasses the undeveloped rumen1. The transition from the pre-ruminant to the ruminant stage in beef calves raised on pasture occurs gradually as consumption of solid feed and rumen fermentation increases and the animal becomes more reliant on microbial protein and VFA2,3. One of the most outstanding features of the ruminant animal is the high gluconeogenic potential from propionate, lactate and AA1. The abrupt removal of a calf from its dam (weaning) at 3 to 6 months of age is a standard practice in beef cattle production systems4. Weaned calves may experience nutritional, social, and psychological stress which causes behavioural, metabolic, physiological, and immunological changes5,6. Therefore, weaning practices should consider animal welfare, metabolic function, and nutritional requirements of the calves2. Previous studies have focussed on changes in behaviour, stress hormones, and immunology of weaned calves, which are more evident during the first 48 h after weaning5,7,8,9. The catabolic effect of stress hormones such as cortisol can explain the increased concentrations of NEFA in blood and some of the decline in ADG after weaning10,11. However, the cessation of milk intake of calves at pasture is expected to produce longer-term metabolic changes because dietary fat is reduced due to the low concentration in forages whereas fiber and ruminal fermentation increases. Surprisingly, little is known about changes in metabolic function of grazing calves after weaning. This information could be used to develop new strategies to alleviate metabolic stress, improve welfare, and minimize the decline in animal production after weaning.

Metabolic function is reflected in the blood metabolome which comprises metabolites such as lipids, sugars and AA that influence cellular, tissue and organ function12,13. The metabolome is downstream of the genome, transcriptome, and proteome, and is the closest ‘omic’ to the phenotype. This led to the suggestion that the metabolome may be the best indicator of alterations in biological function14,15,16. Another important feature of the blood metabolome is that it represents the integration of external (e.g., diet) and internal (e.g. genotype) factors that influence metabolism15. In pigs, weaning was associated with changes in both metabolic function and the blood metabolome and these changes could be partially countered by supplementation with arginine17. A previous study reported that neonatal calves fed colostrum showed an increase in serum abundance of glutamate, histidine, methionine, phenylalanine, tyrosine, tryptophan, valine, leucine, isoleucine and proline, and a reduction of glutamine13. However, no studies were published assessing changes in the blood metabolome of grazing beef calves after weaning under commercial conditions.

The aim of the present study was to characterize changes in the blood metabolome of Angus beef calves subjected to abrupt weaning. This information could improve our understanding of the metabolic impact of abrupt weaning in beef calves and inform novel strategies to manage the transition. The hypothesis was that abrupt weaning in beef calves was associated with stress and nutritional changes that produce metabolic responses reflected in the blood metabolome.

The study had animal ethics approval from The University of Sydney Animal Ethics Committee: Protocol 767. All methods were performed in accordance with ethics approval and the relevant State and National acts, regulations and codes including Animal Research Act 1985 (https://www.dpi.nsw.gov.au/about-us/legislation/list/animal-research), Animal Research Regulation 2021 (https://legislation.nsw.gov.au/view/pdf/asmade/sl-2021-477) and Australian code for the care and use of animals for scientific purposes 8th Edition 2013 (https://www.nhmrc.gov.au/about-us/publications/australian-code-care-and-use-animals-scientific-purposes).

Twenty multiparous Angus cows and their calves (5 to 6 months old) were maintained under standard grazing conditions and management at The University of Sydney John Pye Farm (Greendale, New South Wales, Australia). The pastures were dominated by Paspalum (Paspalum dilatatum) with relatively low densities of Chicory (Cichorium intybus), Plantain (Plantago lanceolata), Red Clover (Trifolium pratense) and White Clover (Trifolium repens) and contained 10.7% CP, 60.7% NDF, and 20.6% ADF. Calves were spring born and had remained with their dam from birth.

On d 0, calves were weighed, blocked by sex, and randomly assigned to a non-weaned group that remained with their dam (NW; n = 10; 257 ± 6 kg BW) or a weaned group that underwent abrupt separation from their dam (W; n = 10; 258 ± 5 kg BW) and sent to a separate paddock with no visual or auditory contact with their dam. Each group had 7 male (non-castrated) and 3 female calves. Pasture in the two paddocks was similar in composition and quantity. This was verified by pasture sampling and analysis at the start (d 0) and end (d 14) of the study. For this purpose, ten quadrats (0.5 m × 0.5 m) were randomly distributed in each paddock and the forage inside the quadrats was cut at ground level, dried for 48 h at 65 °C, and weighed to calculate available forage biomass. The two groups of calves remained in their respective paddock from d 0 to end (d 14) of the study. Body weight was recorded for all calves at d 0, 2, 7 and 14. Blood samples were collected by jugular venipuncture using evacuated tubes containing EDTA (BD Vacutainer, Becton Dickinson, NJ, USA) at d 0, 1, 2, 7, and 14 between 09:00 AM and 10:00 AM. Samples were refrigerated at 4 °C for approximately 30 min, centrifuged (2000×g for 30 min), and plasma stored at − 80 °C until metabolomics analysis.

Behaviour was measured using scan sampling at 3 min intervals for each animal on day − 2, − 1, 0, 1, 2, 7, and 14 in relation to the day of weaning (d 0). Observations were performed directly in the field using a recording sheet containing animal number, date, and time (at 3-min intervals) whereas the behaviour performed by the animals was recorded by two observers during the periods: 7 A.M. to 9 A.M., 11 A.M. to 1 P.M., 2 P.M. to 4 P.M., and 5 P.M. to 7 P.M. The behaviours recorded included grazing, ruminating, resting, suckling, walking, vocalising (mooing), and drinking as described18,19.

Concentrations of cortisol in plasma were determined in a single assay according to the manufacturer's instructions (ADVIA Centaur® Cortisol assay; Siemens Healthcare Pty Ltd, Bayswater, Victoria, Australia). Assay sensitivity was 2.0 ng/mL, and the intra-assay coefficient of variation was 3.5%. The NEFA were also determined in a single assay (Randox Australia Pty Ltd, Parramatta, Australia) previously used in cattle20. Assay sensitivity was 0.03 mEq/L, and the intra-assay coefficient of variation was 1.4%.

Sample preparation and acquisition methods for 1H NMR were based on published protocols21,22. Plasma was allowed to thaw a room temperature and an aliquot (350 µL) was mixed with 350 µL of aqueous (80% H2O:20% D2O) phosphate buffer solution including 0.075 M NaH2PO4, pH = 7.4 (KOH adjusted), 0.1% sodium azide, 1 mM 3-trimethylsilyl-1-[2,2,3,3,-2H4] propionate (TSP) in Eppendorf tubes. Samples were placed on a vortex for 30 s and then centrifuged at 6000×g for 10 min. Aliquots of the supernatant (600 µL) were transferred into 5 mm NMR tubes for 1H NMR analysis.

1H NMR spectra were acquired with a Bruker Avance III 400 MHz spectrometer operating at 400.13 MHz for 1H at 310 K equipped with a 5 mm broad-band inverse configuration probe. Samples were analysed in random order and automation with a SampleCase 24 sample automation system. Samples were analysed using water suppressed 1D NMR spectrum using the NOESYPRESAT pulse sequence (160 transients) and a Carr-Purcell-Meiboom-Gill (CPMG) spin echo sequence with pre-saturation (160 transients). Irradiation of the solvent (water) resonance was applied during pre-saturation delay (2.0 s) for all spectra and for the water suppressed 1D NMR spectra also during the mixing time (0.1 s). The pulse sequence parameters including the 90° pulse (~ 12 µs), pulse frequency offset (~ 1880 Hz), receiver gain (90.5), and pulse powers were optimised on a representative sample and then set constant for the cohort analysed. The spectral width was 30 ppm for NOESY experiments and 20 ppm for CPMG experiments. The resulting Fourier induction decays were processed with an exponential line broadening of 0.3 Hz prior to Fourier transformation, which were collected with approximately 32 k real data points.

1H NMR raw datasets were automatically phased, baseline corrected and referenced to the α-C1H glucose doublet (5.23 ppm) using MATLAB 7.0 software (MathWorks, Natick, MA). To reduce analytical variation between samples, the residual water signal (4.67–4.98 ppm) was truncated from the data set. Probabilistic quotient normalisation was performed across the cohort23. Assignment of endogenous metabolites was made to a high confidence using Chenomx® (Chenomx Inc., Edmonton, AB, Canada) and by reference to published literature, online resources, and spiking experiments24,25,26. Statistical Recoupling of Variables27 was used to cluster pre-processed NMR spectrum data into features or peaks representing the relative metabolite abundance.

Following the processing of the 1H NMR data, multivariate statistical analysis was performed using both MATLAB 7.0 and SAS (SAS Inc, Cary, NC). Data for behaviour, ADG, cortisol, NEFA and the relative abundance of identified metabolites were analysed using a mixed-effects linear regression model with treatment as a fixed effect, time as a repeated measure subjected to the random effect of calf, and treatment × time interaction. Sex was originally included in the model with the respective interaction, but this proved to be insignificant (P > 0.05) and thus deleted from the model. The spatial power covariance structure was selected based on the lowest Bayesian Information Criterion, which accounts for the uneven distance between repeated measures. All variables and residuals were tested for normality, random distribution and mean of zero. Data for cortisol and NEFA were transformed to log10 before analysis to normalise the distribution. The behavioural data was used to calculate the percentage of observations of each behaviour for each animal and day, and then data transformed using the arcsine square root before statistical analysis. Differences between treatment means were calculated within each point in time using pairwise comparisons.

Principal component analysis (PCA) on the normalised spectra was used to identify potential unusual spectra (outliers) and detect obvious patterns or trends in the metabolite profiles28. No outliers were detected and then a second PCA analysis was done with the relative abundance of the 26 identified metabolites at each time point to detect the clustering of treatments groups after plotting principal components. Orthogonal partial least-squares discriminant analysis (OPLS-DA) was used to discriminate animals between the treatment groups (supervised classification). These OPLS-DA models identify spectrums or abundance of metabolites that allow assigning each animal to a particular treatment group. The OPLS-DA models generated were run through 1000 cross-validation permutation testing to assess the validity of the supervised models28,29,30,31. The OPLS-DA models were generated from both the entire 1H NMR normalised spectra and the 26 identified metabolites to determine if a dataset was better at discriminating animals between weaning groups. The predictive ability of OPLS-DA models was measured by the Q2 value30.

No significant differences (P > 0.05; data not shown) were found between treatments in green or dry pasture biomass at the start (2130 ± 217 and 555 ± 67 kg DM, respectively) or end (1753 ± 105 and 668 ± 47 kg DM, respectively) of the trial. The BW on d 0 and d 14 did not differ (P > 0.05) between treatments (data not shown). Overall ADG was greater (P = 0.014) for NW calves (1.05 ± 0.079 kg/d) compared with W calves (0.74 ± 0.084 kg/d). The treatment × day interaction (P < 0.001) indicated that W calves lost weight between d 0 and d 2 while NW calves gained weight in this period (P = 0.06; Fig. 1). The ADG between d 2 and d 7 was greater in W compared to NW calves (P < 0.05) and the opposite effect was found between d 7 and d 14 (P < 0.05; Fig. 1). All behaviours showed a treatment × day interaction (P < 0.001; Fig. 2). Both W and NW calves spent a similar amount of time for each activity on d − 2 and − 1 (P > 0.10). In contrast, W calves spent a shorter time grazing compared to NW calves on d 0 (P < 0.05), d 1 (P = 0.06) and d 7 (P < 0.05) but more on d 2 (P < 0.05). W calves spent less time ruminating and more time vocalising, and walking compared to NW calves on d 0, d 1 and d 2 (P < 0.05; Fig. 2). There were no behavioural differences (P > 0.05) between W and NW calves at d 14 except for a shorter time spent sucking from the dams and longer walking in W compared to NW calves, respectively (P < 0.05).

Growth rate of weaned and not weaned Angus calves. Weaning was at day 0. W is main effect of weaning; D is main effect of Day relative to weaning. †, *, **, *** values within day differ or tended to differ between treatment groups (P ≤ 0.10, P ≤ 0.05, P ≤ 0.01, P ≤ 0.001).

Behaviour of weaned and not weaned Angus beef calves relative to the time of weaning (day 0 was after weaning). W is main effect of weaning; D is main effect of Day relative to weaning. †, *, **, *** values within day differ between treatment groups (P ≤ 0.10, P ≤ 0.05, P ≤ 0.01, P ≤ 0.001).

A weaning × day interaction was observed for plasma concentrations of cortisol (P < 0.05) and NEFA (P < 0.001; Fig. 3). The W calves tended to have greater concentration of cortisol at d 1 (P < 0.10) and 2 (P < 0.01) compared to NW calves but no differences between groups were observed at d 7 and d 14 (P > 0.10). Plasma concentrations of NEFA were greater (P < 0.01) in W compared to NW calves on d 1, d 2 and d 7 (P < 0.05), with no differences at d 0 and d 14 (P > 0.10; Fig. 3).

Blood concentration of cortisol and non-esterified fatty acids (NEFA) for weaned and not weaned Angus beef calves relative to the time of weaning at day 0 (blood sample before weaning). W is main effect of weaning; D is main effect of Day relative to weaning. †, *, **, *** values within day differ between treatment groups (P ≤ 0.10, P ≤ 0.05, P ≤ 0.01, P ≤ 0.001).

A total of 26 metabolites were identified from the 1H NMR spectra (Figs. 4, 5, and 6). The relative concentration of metabolites for NW calves did not change over time (P > 0.05) whereas all metabolites changed over time for W calves (P < 0.05). Differences between NW and W calves emerged after weaning for 24 of the 26 metabolites. The exceptions were the AA histidine and glutamine which did not show a treatment × day interaction (P > 0.05). Some metabolites showed the greatest difference between weaning groups at d 1 and d 2 and others at d 14 (Figs. 4, 5, and 6). Weaned calves had lesser abundance of tyrosine and greater abundance of isoleucine compared to NW calves on d 1 or d 2 (P < 0.10) with no differences at d 0, d 7 and d 14 (P > 0.10). Betaine was greater in W compared to NW calves on d 2 and d 7 (P < 0.05; Fig. 6). The relative abundance of other metabolites was greater for W compared to NW calves only at d 7 or d 14, or both d 7 and d 14 including a group of AA (valine, alanine, threonine, leucine, lysine, threonine), allantoin, acetate, and glycoprotein; and lower abundance in W compared to NW calves for VLDL/LDL and unsaturated lipids (P < 0.05). Another group of metabolites showed greater relative abundance in W compared to NW calves from d 1 or d 2 to d 14 including phenylalanine, glutamate, creatinine, creatine, creatine phosphate, glucose, 3-hydroxybutyrate, and 3-hydroxyisobutyrate (P < 0.10; Figs. 4, 5, 6).

Relative concentration of plasma metabolites of weaned and not weaned Angus beef calves relative to the time of weaning at day 0 (blood sample before weaning). W is main effect of weaning; D is main effect of Day relative to weaning. †, *, **, *** values within day differ between treatment groups (P ≤ 0.10, P ≤ 0.05, P ≤ 0.01, P ≤ 0.001).

The unsupervised PCA plots with the 26 identified metabolites as predictors in Fig. 7 represent the first principal component (PC1, X-axis) which accounts for as much of the variability in the data as possible plotted against the second principal component (PC2, Y-axis). The PC2 accounts for the greatest possible variance in the data under the constraint that it is orthogonal to the preceding component. The variance explained by the first two PCs was over 50% for any day in relation to weaning (Fig. 7). Before weaning (d 0), the PCA did not display an obvious discrimination between NW and W calves based on their metabolomes as shown by the overlap of the two groups of calves (Fig. 7A). At d 1 after weaning, the PCA showed a separation between NW and W calves (Fig. 7B). This separation remained evident at d 2 (Fig. 7C), d 7 (Fig. 7D) and d 14 (Fig. 7E).

Principal Component Analysis of the 1H NMR spectra of non-weaned (open circles) and weaned (closed circles) Angus calves (A) immediately before abrupt weaning, (B) d 1, (C) d 2, (D) d 7, and (E) d 14 after weaning. The ellipses serve to highlight the lack of clustering at d 0 and the clustering of the non-weaned (dotted outline) and the weaned group (solid outline) on different days after weaning.

A multi-vector supervised OPLS-DA scatter plot using the entire spectrum of each sample also highlighted the differences in the blood metabolome of calves depending on time and weaning status (Fig. 8). The blood metabolome was markedly altered at d 1 after weaning which continued to d 14 after weaning. The increasing separation between W and NW calves with time indicated that the metabolome continued to diverge with increasing time after weaning.

Stability of the metabotypes of Angus calves related to time and weaning treatment group (weaned or non-weaned calves). Two dimensional OPLS model scores plot (2 predictive and 5 orthogonal components) derived from 1H NMR spectra from the plasma of Angus calves (time point, dark blue through to red) that were non-weaned (open circles) and weaned (closed circles). Time points are d 0 (purple), 1 (blue), 2 (green), 7 (orange), and 14 (red) after weaning.

Separate discriminatory models were subsequently created for each day of blood collection and either using the entire 1H NMR spectrum or the 26 identified metabolites (Table 1). At d 0, the Q2 value indicated that the OPLS-DA had no predictive validity of whether calves belonged to the NW or W groups (Table 1). The Q2 values for d 1 to 14 were > 0.40 which confirmed that the OPLS-DA could accurately discriminate between NW and W calves using both the entire spectrum and the relative abundance of metabolites (Table 1). Predictions using the metabolites seemed to perform slightly better than the entire spectrum due the slightly higher Q2. Random permutation testing confirmed all models were able to differentiate between treatment groups after the change in weaning status.

The present study sought to characterize the effect of abrupt weaning on the blood metabolome of beef calves measured using 1H NMR and supported by traditional measurements of stress and energy balance such as ADG, cortisol, NEFA, and behaviour. The PCA at each sampling date was performed to assess potential clustering of calves into either a W or NW group according to the metabolome as an unsupervised classification method. This analysis confirmed that weaning had major effects on metabolic function to allow the separation of groups based purely on the blood metabolome. These results were further confirmed using supervised classification (OPLS-DA), which showed that accurate discrimination can be achieved either using the entire NMR spectrum or the identified metabolites. Furthermore, both PCA and OPLS-DA analysis demonstrated that the metabolome of W and NW calves did not differ before weaning but a clear divergence of groups existed since the time of weaning. This is the first published study to assess the changes in the blood metabolome of beef calves after weaning to the authors’ best knowledge. Differences in the blood metabolome between NW and W calves were not explained by pasture total or green biomass available, which was similar for the two groups. Therefore, the metabolic changes observed in W calves are likely due to a combination of physiological and psychological stressors, type of nutrients consumed, and changes in metabolic function as a result of the discontinuation of milk intake after physical separation of the calf from the dam, cessation of suckling, removal of social support from the dam, and change of social organisation2,4,32. It is important to note that calves of the present study were 5 to 6 months-old and the development of the reticulorumen is expected to be completed with solid feed available at this age33. However, calves still suckled during the present study suggesting that a proportion of the energy and nutrients consumed came from the dam's milk. Milk consumed bypasses the rumen to be digested in the abomasum and absorbed in the lower gastrointestinal tract33. The reduction in plasma relative concentration of VLDL/LDL and unsaturated lipids in W calves from d 0 to d 14 may reflect the abrupt reduction in fat intake from milk and reliance on forage intake. Forages have very low fat content (~ 1% of DM) whereas bovine milk has high fat content of approximately 30% of DM33.

In contrast to the long-term decline of lipids, the marked differences between W and NW calves during the first week after weaning in the present study could reflect physiological changes from stress hormones. Prior research concluded that the stress of weaning was acute and of short duration (up to 48 h) if the conditions are appropriate2. Findings of the present study for ADG, cortisol, and behaviour (vocalising, walking and grazing) support this observation because these variables returned to pre-weaning values by d 7 in agreement with previous studies19. Elevated blood cortisol influences a broad range of physiological processes in cattle including metabolism34,35 and immune function5,6,36. In the present study, the concentration of NEFA, betaine and 3-hydroxybutyrate also peaked at d 2 after weaning but remained elevated until d 7 and 14, respectively. In contrast, cortisol and behaviours indicative of stress had already returned to pre-weaning values by d 7. Elevated NEFA, betaine, and 3-hydroxybutyrate in cattle generally reflect an increase in lipolysis and body fat mobilisation which can be caused by stress-mediated glucocorticoids, high energy demand, or negative energy balance2,37. The W calves lost weight at d 2 after weaning however both NW and W calves gained weight thereafter indicating there was no negative energy balance. Therefore, the finding that W calves had lower ADG compared to NW calves suggests a cortisol-mediated effect on energy utilisation or a nutritional challenge from lower nutrient intake, or both.

The temporal trends in the abundance of NEFA for W calves differed markedly from those of lipids suggesting that these compounds do not share the same metabolic pathway despite both being part of lipid metabolism. This speculation is supported by a lack of effect of increasing fat concentration in milk replacer on blood NEFA and triglycerides of calves38. Furthermore, there was no correlation between NEFA and lipids in the present study (P > 0.05; data not shown). Therefore, lipids measured by 1H NMR may derive from dietary sources whereas NEFA may be derived from fat mobilisation.

Betaine is a methyl donor that participates in protein and lipid metabolism, improves the antioxidant capacity to protect cells from stress, and has osmoregulatory properties in cells39. In transition dairy cows, the production of betaine increased under negative energy balance, leading to an increase in plasma concentration of NEFA and 3-hydroxybutyrate to support energy production via lipolysis when plasma glucose concentration decreased40. Furthermore, betaine supplementation during periods of heat stress improved milk production in dairy cattle41,42 and growth and fat deposition in beef cattle43,44. It is considered likely that W calves of the present study increased betaine production to improve lipolysis and oxidation of triglycerides and thus, supply energy required for the stress response of weaning. However, the present findings differ from lactating dairy cows because calves were not in negative energy balance and glucose concentration increased in W compared to NW calves, unlike transition dairy cows.

The relative abundance of tyrosine decreased in W calves in the first 2 days after weaning in the present study, suggesting tyrosine also played an important role in the stress response to weaning. Tyrosine is a precursor of several neurotransmitters of the catecholamine pathway involved in the stress response including dopamine, norepinephrine and epinephrine45,46,47. It is hypothesized that the acute stress response to weaning increased the demand for and utilisation of tyrosine leading to a reduction of its concentration in blood. However, research on the role of tyrosine in the stress response of cattle is scarce and more research is required to draw more informed conclusions.

Glycoproteins are acute phase proteins that increase with inflammation, infection, stress, and trauma48. The increase in the relative abundance of glycoprotein in W calves compared to NW calves may have reflected the general stress response of W calves to abrupt weaning. However, acute phase proteins are generally not specific, and it is unclear if glycoproteins measured by NMR are the same as those measured by other techniques such as ELISA or chromatography. Nevertheless, these findings are important because demonstrate the potential of blood metabolomics to measure multiple metabolites that indicate different responses and metabolic pathways.

In contrast to the metabolites that seemed to reflect the acute stress of weaning, other metabolites reached their greatest relative concentration at d 14 such as valine, alanine, leucine, lysine, creatinine, creatine, creatine phosphate, acetate, and glucose. It is suggested that these metabolites reflect a change in metabolic function either due to diet change, chronic stress, or both. The fact that no significant changes in the metabolome of NW calves occurred over time in the present study suggest that these metabolites reflect metabolic adjustments of W calves. However, the present study only measured the metabolome to 14 d after weaning and it would be of interest to measure metabolites for a longer period after weaning to determine if these changes are a result of long-term metabolic changes with the cessation of milk intake or, alternatively, short-term metabolic adjustments or a consequence of the acute stress response.

The sharp and sustained increase in circulating glucose of W compared to NW calves could be interpreted to reflect an increase of gluconeogenesis either because of the stress response or metabolic adjustments to produce more energy from ruminal volatile fatty acids and gluconeogenic AA rather than from milk, or both. On one hand, glucocorticoids promote glucogenic metabolism during stress through the mobilisation of glycogen stores in liver and muscle45,49,50,51. However, cortisol and stress behaviours peaked at d 2 after weaning and returned to pre-weaning values by d 7 whereas glucose peaked at d 14 in the present study. On the other hand, the rate of gluconeogenesis increases with the transition from the pre-ruminant to the ruminant stage1. All calves of the present study consumed medium quality pasture and therefore dietary glucose was expected to be minimal. Gluconeogenesis in the liver may be supported by propionate absorbed from the rumen, lactose, lactate, and glucogenic AA such as alanine and glutamine1, and valine, isoleucine, and threonine47. Alanine, valine and threonine appeared to have supported gluconeogenesis in the present study increasing to d 14 like glucose. Propionate was not identified in the NMR spectra although acetate showed a large increase at d 7 and 14. However, lactate, glutamine, and isoleucine showed no differences between W and NW calves at d 7 or d 14.

Acetate is a VFA produced in higher proportion during fermentation of high forage diets in the rumen where it is absorbed into the bloodstream52. Therefore, the greater abundance of acetate in W compared to NW calves at d 7 and d 14 could have been due to higher intake of forage and fibre in W calves (partly substituting milk). After absorption, acetate can be converted to 3-hydroxybutyrate, oxidized via the tricarboxylic acid cycle (TCA) or used for fatty acid synthesis53. However, the concentration of acetate did not peak at the same time as 3-hydroxybutyrate in the present study.

The temporal trend of 3-hydroxybutyrate of W calves in the present study was interesting because it peaked at d 2 but remained elevated until d 14. This metabolite can originate from either incomplete oxidation of fatty acids during lipolysis or the conversion of acetate during absorption from the rumen54. Therefore, it is possible that the trends observed in W calves reflected both metabolic processes, with greater abundance of 3-hydroxybutyrate during the first 2 days after weaning reflecting incomplete hepatic oxidation of fatty acids (NEFA) mobilised during the acute stress response. The relative abundance of 3-hydroxybutyrate could has then been sustained by the conversion of acetate absorbed in the rumen due to increased forage intake. This is supported by results from previous studies which demonstrated that 3-hydroxybutyrate increased during periods of high energy demand or negative energy balance in post-partum dairy cows or calves after transport and weaning37,55,56. It is important to note that the relative contribution of different metabolic pathways to the observed trends and abundance of circulating metabolites cannot be measured with the experimental design of the present study. It is also possible that other metabolites also reflect the interaction between multiple metabolic pathways such as glucose, creatine, phenylalanine, glutamate, and 3-hydroxyisobutyrate.

Creatine, creatine phosphate and creatinine of W calves were maximally different to NW calves at d 14. These metabolites are critical in energy metabolism of the ruminant animal to store and recycle energy. Elevated creatine kinase, creatine and creatinine are associated with excessive muscle breakdown, muscle damage during periods of increased or unaccustomed physical activity and stress45,46,47. Creatine is one of the main energy sources for brain tissue and skeletal muscle during physical activity, which is converted to creatine phosphate by the enzyme creatine kinase and stored as a rich energy source47. This reaction is reversible and thus creatine phosphate can be converted to creatine to produce energy (ATP), and creatine can then be converted to creatinine in an irreversible reaction which is an inactive end product of muscle metabolism excreted in urine47,57. For this reason, creatinine is often used an indicator of proteolysis and body muscle mass, and typically decreases in blood concurrent with increases of 3-hydroxybutyrate and NEFA and with a reduction of glucose in post-partum dairy cows losing BW under negative energy balance58,59. However, these changes observed in dairy cows are in contrast with the results in weaned calves of the present study where all three metabolites increased to d 14 after weaning. The W calves of the present study showed increased physical activity immediately post-weaning at d 1 and d 2 but declined at d 7 and d 14, in agreement with previous studies9,19. Therefore, the increased abundance of creatine and creatinine of W calves to d 14 appears not to be explained by physical activity, behavioural indicators of stress, cortisol, or NEFA. It could be possible that the production of creatine and storage of creatine phosphate continues after a period of acute stress or high physical activity leading to ‘chronic’ elevation of these metabolites. However, this latter hypothesis cannot be confirmed with data from the present study and further research is needed to determine whether such profound and long-term changes in energy and protein metabolism are a result of earlier acute stress resulting in chronic elevation of creatine and creatinine or to metabolic changes due to discontinuation of milk intake. Regardless, W calves seemed to have been in a catabolic situation until d 14 with both proteolysis and lipolysis occurring as reflected by high relative concentrations of creatine, creatinine, 3-hydroxybutyrate, glucose, and some AA.

Allantoin increased at d 7 and d 14 in W compared to NW calves in the present study. Allantoin is one of the purine derivatives which increases with the flow of microbial protein from the rumen and with OM intake60, and it is also involved in the urea cycle and excreted in urine47. In addition, allantoin has been suggested as a biomarker of chronic oxidative stress, illness and senescence in both cattle and humans61,62. Chronic stress was not evident in the present study, and it is hypothesized that the greater allantoin concentration of W compared to NW calves is a result of greater solid feed intake and ruminal microbial activity after discontinuation of milk intake. However, feed intake was not measured in the present study and thus further research is needed to confirm this hypothesis. Similarly, the reasons for the sustained increase in valine, alanine, threonine, leucine, lysine, phenylalanine and glutamate in W calves to d 14 are difficult to explain in the present study because total protein intake was not measured (from both milk and pasture). Faecal N concentration measured on the same days of blood sampling as an indicator of protein intake showed no differences between W and NW calves at any time point with an overall average of 1.12 ± 0.033 and 1.15 ± 0.032% of DM for W and NW calves, respectively (data not shown; P > 0.05). This could imply that total protein flow to the lower gastrointestinal tract was not different between W and NW calves.

Previous research that measured metabolites in blood were undertaken in calves fed milk replacers which did not experience the stress from separation of the dam as in the present study. In addition, milk allowance was gradually reduced in most previous studies with dairy63 and beef calves64 and these studies have been done in small pens where animals cannot fully express their stress behaviours such as walking long distances. For instance, discontinuation of feeding milk replacer showed no effects on glucose or NEFA in blood64,65, and increased glucose, acetate, propionate and butyrate, and a decrease in NEFA63. Interestingly, 3-hydroxybutyrate increased after weaning in all the above studies and there is agreement this metabolite comes from the conversion of acetate and butyrate in the ruminal epithelium as a result of greater solid feed intake fermented in the rumen of weaned calves63. A study with beef calves at foot on pasture that were placed in pens at weaning produced similar results to the present study, where both NEFA and 3-hydroxybutyrate increased after weaning although glucose and creatine kinase were not affected3 unlike the present study. It is important to note that the latter study did not have a control NW group like the present study. A review of the effects of weaning on calves concluded that circulating glucose produced inconsistent results across studies but the causes for these discrepancies were not identified2.

The present study confirmed that blood metabolomics is a strong platform to investigate and understand metabolic changes resulting from complex interactions that arise from stress and dietary changes. The combination of high throughput 1H NMR and advanced statistics is a powerful combination that generates a large amount of biological information. The field of metabolomics can make an important contribution to increasing the understanding of complex metabolic changes. This could facilitate the development of strategies that enhance the management, performance and welfare of production animals. For example, supplements for calves at weaning could be formulated with high concentration of precursors of the metabolites required for gluconeogenesis or the stress response such as tyrosine. Furthermore, the nutritional requirements of weaned calves could be managed using new technologies to deliver targeted amounts and types of supplements66. These approaches can also be integrated with strategies that reduce physical activity such as yard weaning, fenceline contact, or two-stage weaning which are known to lower the stress response and energy expenditure8,9.

In conclusion, several key metabolites showed dramatic changes in calves to d 14 after weaning reflecting marked changes in carbohydrate, lipid, and protein metabolism due to both stress and diet change. Metabolic changes due to stress seemed of short duration being most marked on d 2 after weaning, whereas longer-term metabolic changes were observed until d 14 after weaning. The drastic short-term changes in lipid metabolism of W calves were reflected through a sharp increase of NEFA and betaine as methyl donor to facilitate lipid oxidation. The metabolite 3-hydroxybutyrate seemed to be involved in both short- and long-term changes in lipid metabolism peaking at d 2 and remaining high until d 14 after weaning. Long-term changes in lipid metabolism were shown through lower concentration of LDL/VLDL and unsaturated lipids, which were most evident at d 14 after weaning and could be the result of diet change from discontinuation of milk intake. The most important metabolic change associated with carbohydrate metabolism was the significant increase in glucose concentration in plasma of W calves until d 14, which could have been due to an increase of gluconeogenesis. In addition, a sharp increase in acetate observed on d 7 and 14 could reflect greater fermentation of solid feed in the rumen which may also have contributed to a sustained high abundance of 3-hydroxybutyrate of W calves. Changes in protein metabolism were evidenced through an increase in several AA and allantoin up until d 14. Finally, profound long-term changes in energy metabolism of peripheral muscle and brain tissue were reflected through a large increase in the relative concentration of circulating creatine, creatine phosphate, and creatinine in W calves up to d 14. Further research is needed to ascertain whether the long-term metabolic changes of calves after weaning are a permanent shift in metabolism or a transient change because of acute or chronic stress, and whether metabolite concentrations return to pre-weaning values.

Data is available upon reasonable request to the corresponding author.

Aschenbach, J. R., Kristensen, N. B., Donkin, S. S., Hammon, H. M. & Penner, G. B. Gluconeogenesis in dairy cows: The secret of making sweet milk from sour dough. IUBMB Life 62, 869–877. https://doi.org/10.1002/iub.400 (2010).

Article CAS PubMed Google Scholar

Lynch, E., McGee, M. & Earley, B. Weaning management of beef calves with implications for animal health and welfare. J. Appl. Anim. Res. 47, 167–175. https://doi.org/10.1080/09712119.2019.1594825 (2019).

Article Google Scholar

Lynch, E. M., McGee, M., Doyle, S. & Earley, B. Effect of pre-weaning concentrate supplementation on peripheral distribution of leukocytes, functional activity of neutrophils, acute phase protein and behavioural responses of abruptly weaned and housed beef calves. BMC Vet. Res. 8, 1. https://doi.org/10.1186/1746-6148-8-1 (2012).

Article CAS PubMed PubMed Central Google Scholar

Enriquez, D., Hotzel, M. J. & Ungerfeld, R. Minimising the stress of weaning of beef calves: A review. Acta Vet. Scand. https://doi.org/10.1186/1751-0147-53-28 (2011).

Article PubMed PubMed Central Google Scholar

Lynch, E. M., Earley, B., McGee, M. & Doyle, S. Effect of abrupt weaning at housing on leukocyte distribution, functional activity of neutrophils, and acute phase protein response of beef calves. BMC Vet. Res. https://doi.org/10.1186/1746-6148-6-39 (2010).

Article PubMed PubMed Central Google Scholar

O’Loughlin, A., McGee, M., Doyle, S. & Earley, B. Biomarker responses to weaning stress in beef calves. Res. Vet. Sci. 97, 458–463. https://doi.org/10.1016/j.rvsc.2014.06.003 (2014).

Article CAS PubMed Google Scholar

O’Loughlin, A., McGee, M., Waters, S. M., Doyle, S. & Earley, B. Examination of the bovine leukocyte environment using immunogenetic biomarkers to assess immunocompetence following exposure to weaning stress. BMC Vet. Res. 7, 45. https://doi.org/10.1186/1746-6148-7-45 (2011).

Article CAS PubMed PubMed Central Google Scholar

Haley, D. B., Bailey, D. W. & Stookey, J. M. The effects of weaning beef calves in two stages on their behavior and growth rate. J. Anim. Sci. 83, 2205–2214 (2005).

Article CAS PubMed Google Scholar

Price, E. O., Harris, J. E., Borgwardt, R. E., Sween, M. L. & Connor, J. M. Fenceline contact of beef calves with their dams at weaning reduces the negative effects of separation on behavior and growth rate. J. Anim. Sci. 81, 116–121 (2003).

Article CAS PubMed Google Scholar

Lefcourt, A. M. & Elsasser, T. H. Adrenal responses of angus X hereford cattle to the stress of weaning. J. Anim. Sci. 73, 2669–2676 (1995).

Article CAS PubMed Google Scholar

Qiu, X. et al. Genetic effects on acute phase protein response to the stresses of weaning and transportation in beef calves. J. Anim. Sci. 85, 2367–2374. https://doi.org/10.2527/jas.2006-843 (2007).

Article CAS PubMed Google Scholar

Goldansaz, S. A. et al. Livestock metabolomics and the livestock metabolome: A systematic review. PLoS ONE https://doi.org/10.1371/journal.pone.0177675 (2017).

Article PubMed PubMed Central Google Scholar

Qi, Y. et al. Exploration of the relationship between intestinal colostrum or milk, and serum metabolites in neonatal calves by metabolomics analysis. J. Agric. Food Chem. 66, 7200–7208. https://doi.org/10.1021/acs.jafc.8b01621 (2018).

Article CAS PubMed Google Scholar

Pantophlet, A. J. et al. The use of metabolic profiling to identify insulin resistance in veal calves. PLoS ONE https://doi.org/10.1371/journal.pone.0179612 (2017).

Article PubMed PubMed Central Google Scholar

Fontanesi, L. Metabolomics and livestock genomics: Insights into a phenotyping frontier and its applications in animal breeding. Anim. Front. 6, 73–79. https://doi.org/10.2527/af.2016-0011 (2016).

Article Google Scholar

D’Occhio, M. J., Baruselli, P. S. & Campanile, G. Influence of nutrition, body condition, and metabolic status on reproduction in female beef cattle: A review. Theriogenology 125, 277–284. https://doi.org/10.1016/j.theriogenology.2018.11.010 (2019).

Article CAS PubMed Google Scholar

He, Q. H. et al. Dietary supplementation with l-arginine partially counteracts serum metabonome induced by weaning stress in piglets. J. Proteome Res. 10, 5214–5221. https://doi.org/10.1021/pr200688u (2011).

Article ADS CAS PubMed Google Scholar

Gonzalez, L. A., Bishop-Hurley, G. J., Handcock, R. N. & Crossman, C. Behavioral classification of data from collars containing motion sensors in grazing cattle. Comput. Electron Agr. 110, 91–102. https://doi.org/10.1016/j.compag.2014.10.018 (2015).

Article Google Scholar

Alvez, P., Quintans, G., Hotzel, M. J. & Ungerfeld, R. Two-step weaning in beef calves: Permanence of nose flaps for 7 or 21 days does not influence the behaviour response. Anim. Prod. Sci. 56, 866–870. https://doi.org/10.1071/An14643 (2016).

Article Google Scholar

Wathes, D. C. et al. Differences between primiparous and multiparous dairy cows in the inter-relationships between metabolic traits, milk yield and body condition score in the periparturient period. Domest. Anim. Endocrin. 33, 203–225. https://doi.org/10.1016/j.domaniend.2006.05.004 (2007).

Article CAS Google Scholar

Dona, A. C. et al. Precision high-throughput proton NMR spectroscopy of human urine, serum, and plasma for large-scale metabolic phenotyping. Anal. Chem. 86, 9887–9894. https://doi.org/10.1021/ac5025039 (2014).

Article CAS PubMed Google Scholar

Connolly, S., Dona, A., Hamblin, D., D’Occhio, M. J. & Gonzalez, L. A. Changes in the blood metabolome of Wagyu crossbred steers with time in the feedlot and relationships with marbling. Sci. Rep. 10, 18987. https://doi.org/10.1038/s41598-020-76101-6 (2020).

Article CAS PubMed PubMed Central Google Scholar

Craig, A., Cloareo, O., Holmes, E., Nicholson, J. K. & Lindon, J. C. Scaling and normalization effects in NMR spectroscopic metabonomic data sets. Anal. Chem. 78, 2262–2267. https://doi.org/10.1021/ac0519312 (2006).

Article CAS PubMed Google Scholar

Dona, A. C. et al. A guide to the identification of metabolites in NMR-based metabonomics/metabolomics experiments. Comput. Struct. Biotechnol. 14, 135–153. https://doi.org/10.1016/j.csbj.2016.02.005 (2016).

Article CAS Google Scholar

Nicholson, J. K., Foxall, P. J. D., Spraul, M., Farrant, R. D. & Lindon, J. C. 750-MHz H-1 and H-1-C-13 NMR-spectroscopy of human blood-plasma. Anal. Chem. 67, 793–811. https://doi.org/10.1021/ac00101a004 (1995).

Article CAS PubMed Google Scholar

Psychogios, N. et al. The human serum metabolome. PLoS ONE https://doi.org/10.1371/journal.pone.0016957 (2011).

Article PubMed PubMed Central Google Scholar

Blaise, B. J. et al. Statistical recoupling prior to significance testing in nuclear magnetic resonance based metabonomics. Anal. Chem. 81, 6242–6251. https://doi.org/10.1021/ac9007754 (2009).

Article CAS PubMed Google Scholar

Lindon, J. C., Holmes, E. & Nicholson, J. K. Pattern recognition methods and applications in biomedical magnetic resonance. Prog. Nucl. Mag. Res. Sp. 39, 1–40. https://doi.org/10.1016/S0079-6565(00)00036-4 (2001).

Article CAS Google Scholar

Gromski, P. S. et al. A tutorial review: Metabolomics and partial least squares-discriminant analysis—A marriage of convenience or a shotgun wedding. Anal. Chim. Acta 879, 10–23. https://doi.org/10.1016/j.aca.2015.02.012 (2015).

Article CAS PubMed Google Scholar

Triba, M. N. et al. PLS/OPLS models in metabolomics: The impact of permutation of dataset rows on the K-fold cross-validation quality parameters. Mol. BioSyst. 11, 13–19. https://doi.org/10.1039/c4mb00414k (2015).

Article CAS PubMed Google Scholar

Trygg, J. & Wold, S. Orthogonal projections to latent structures (O-PLS). J. Chemometr. 16, 119–128. https://doi.org/10.1002/cem.695 (2002).

Article CAS Google Scholar

Lyles, J. L., Calvo-Lorenzo, M. S. & Bill, E. Kunkle Interdisciplinary Beef Symposium: Practical developments in managing animal welfare in beef cattle: what does the future hold?. J. Anim. Sci. 92, 5334–5344. https://doi.org/10.2527/jas.2014-8149 (2014).

Article CAS PubMed Google Scholar

Huber, J. T. Development of the digestive and metabolic apparatus of the calf. J. Dairy Sci. 52, 1303–1315. https://doi.org/10.3168/jds.S0022-0302(69)86744-5 (1969).

Article CAS PubMed Google Scholar

Burdick, N. C. et al. Relationships between cortisol concentrations and cow temperament with calf exit velocity from 3 weeks of age through weaning. J. Dairy Sci. 90, 528–528 (2007).

Google Scholar

Saco, Y. et al. Evaluation of serum cortisol, metabolic parameters, acute phase proteins and faecal corticosterone as indicators of stress in cows. Vet. J. 177, 439–441. https://doi.org/10.1016/j.tvjl.2007.05.019 (2008).

Article CAS PubMed Google Scholar

Carroll, J. A. & Forsberg, N. E. Influence of stress and nutrition on cattle immunity. Vet. Clin. N. Am. Food A 23, 105. https://doi.org/10.1016/j.cvfa.2007.01.003 (2007).

Article Google Scholar

Gonzalez, L. A. et al. Effect of the number of concentrate feeding places per pen on performance, behavior, and welfare indicators of Friesian calves during the first month after arrival at the feedlot. J. Anim. Sci. 86, 419–431. https://doi.org/10.2527/jas.2007-0362 (2008).

Article CAS PubMed Google Scholar

Yohe, T. T. et al. Metabolic and performance responses to the replacement of lactose by fat in milk replacer formulations for dairy calves. Animal 15, 100031. https://doi.org/10.1016/j.animal.2020.100031 (2021).

Article CAS PubMed Google Scholar

McFadden, J. W. et al. Symposium review: One-carbon metabolism and methyl donor nutrition in the dairy cow. J. Dairy Sci. 103, 5668–5683. https://doi.org/10.3168/jds.2019-17319 (2020).

Article CAS PubMed Google Scholar

Coleman, D. N. et al. Choline supply during negative nutrient balance alters hepatic cystathionine beta-synthase, intermediates of the methionine cycle and transsulfuration pathway, and liver function in Holstein cows. J. Dairy Sci. 102, 8319–8331. https://doi.org/10.3168/jds.2019-16406 (2019).

Article CAS PubMed Google Scholar

Monteiro, A. P. A. et al. Effects of feeding betaine-containing liquid supplement to transition dairy cows. J. Dairy Sci. 100, 1063–1071. https://doi.org/10.3168/jds.2016-11452 (2017).

Article CAS PubMed Google Scholar

Wang, C., Liu, H., Wang, C., Liu, J. & Liu, H. Effects of dietary rumen-protected betaine on lactation performance and serum metabolites of mid-lactation Holstein dairy cows. J. Agric. Food Chem. 68, 13154–13159. https://doi.org/10.1021/acs.jafc.9b07453 (2020).

Article CAS PubMed Google Scholar

DiGiacomo, K., Warner, R. D., Leury, B. J., Gaughan, J. B. & Dunshea, F. R. Dietary betaine supplementation has energy-sparing effects in feedlot cattle during summer, particularly in those without access to shade. Anim. Prod. Sci. 54, 450–458. https://doi.org/10.1071/AN13418 (2014).

Article CAS Google Scholar

Wang, C. et al. Effects of rumen-protected folic acid and betaine supplementation on growth performance, nutrient digestion, rumen fermentation and blood metabolites in Angus bulls. Br. J. Nutr. 123, 1109–1116. https://doi.org/10.1017/S0007114520000331 (2020).

Article CAS PubMed Google Scholar

Schaefer, A. L., Dubeski, P. L., Aalhus, J. L. & Tong, A. K. W. Role of nutrition in reducing antemortem stress and meat quality aberrations. J. Anim. Sci. 79, E91–E101. https://doi.org/10.2527/jas2001.79E-SupplE91x (2001).

Article Google Scholar

Lieberman, H. R., Georgelis, J. H., Maher, T. J. & Yeghiayan, S. K. Tyrosine prevents effects of hyperthermia on behavior and increases norepinephrine. Physiol. Behav. 84, 33–38. https://doi.org/10.1016/j.physbeh.2004.10.023 (2005).

Article CAS PubMed Google Scholar

Engelking, L. R. Textbook of Veterinary Physiological Chemistry (Academic Press/Elsevier, 2015).

Google Scholar

Cray, C. Acute phase proteins in animals. Prog. Mol. Biol. Transl. Sci. 105, 113–150. https://doi.org/10.1016/B978-0-12-394596-9.00005-6 (2012).

Article CAS PubMed Google Scholar

Bueno, A. R., Rasby, R. & Clemens, E. T. Age at weaning and the endocrine response to stress. Arq. Bras. Med. Vet. Zoo. 55, 1–7. https://doi.org/10.1590/S0102-09352003000100001 (2003).

Article Google Scholar

Gardner, G. E., McGilchrist, P. & Pethick, D. W. Ruminant glycogen metabolism. Anim. Prod. Sci. 54, 1575–1583. https://doi.org/10.1071/An14434 (2014).

Article CAS Google Scholar

Yanchev, I. et al. Effect of weaning on some plasma metabolites in black-and-white calves fed diets with or without supplemental chromium picolinate. Bulg. J. Agric. Sci. 14, 329–334 (2008).

Google Scholar

Sicilianojones, J. & Murphy, M. R. Production of volatile fatty-acids in the rumen and cecum-colon of steers as affected by forage—concentrate and forage physical form. J. Dairy Sci. 72, 485–492. https://doi.org/10.3168/jds.S0022-0302(89)79130-X (1989).

Article CAS PubMed Google Scholar

Church, D. C. The Ruminant Animal: Digestive Physiology and Nutrition (Waveland Press, 1993).

Google Scholar

Hocquette, J. F. et al. Intramuscular fat content in meat-producing animals: Development, genetic and nutritional control, and identification of putative markers. Animal 4, 303–319. https://doi.org/10.1017/S1751731109991091 (2010).

Article CAS PubMed Google Scholar

McArt, J. A. A., Nydam, D. V., Oetzel, G. R., Overton, T. R. & Ospina, P. A. Elevated non-esterified fatty acids and beta-hydroxybutyrate and their association with transition dairy cow performance. Vet. J. 198, 560–570. https://doi.org/10.1016/j.tvjl.2013.08.011 (2013).

Article CAS PubMed Google Scholar

Ospina, P. A., Mcart, J. A., Overton, T. R., Stokol, T. & Nydam, D. V. Using nonesterified fatty acids and beta-hydroxybutyrate concentrations during the transition period for herd-level monitoring of increased risk of disease and decreased reproductive and milking performance. Vet. Clin. N. Am.-Food A 29, 387. https://doi.org/10.1016/j.cvfa.2013.04.003 (2013).

Article Google Scholar

Wyss, M. & Kaddurah-Daouk, R. Creatine and creatinine metabolism. Physiol. Rev. 80, 1107–1213 (2000).

Article CAS PubMed Google Scholar

Spaans, O. K. et al. Temporal profiles describing markers of inflammation and metabolism during the transition period of pasture-based, seasonal-calving dairy cows. J. Dairy Sci. https://doi.org/10.3168/jds.2021-20883 (2022).

Article PubMed Google Scholar

Sadri, H., Ghaffari, M. H. & Sauerwein, H. Invited review: Muscle protein breakdown and its assessment in periparturient dairy cows. J. Dairy Sci. 106, 822–842. https://doi.org/10.3168/jds.2022-22068 (2023).

Article CAS PubMed Google Scholar

Gonda, H. L. & Lindberg, J. E. Effect of diet on milk allantoin and its relationship with urinary allantoin in dairy cows. J. Dairy Sci. 80, 364–373. https://doi.org/10.3168/jds.S0022-0302(97)75946-0 (1997).

Article CAS PubMed Google Scholar

Kand’ar, R. & Zakova, P. Allantoin as a marker of oxidative stress in human erythrocytes. Clin. Chem. Lab. Med. 46, 1270–1274. https://doi.org/10.1515/Cclm.2008.244 (2008).

Article CAS PubMed Google Scholar

Liao, Y. et al. Metabolomics profiling of serum and urine in three beef cattle breeds revealed different levels of tolerance to heat stress. J. Agric. Food Chem. 66, 6926–6935. https://doi.org/10.1021/acs.jafc.8b01794 (2018).

Article CAS PubMed Google Scholar

Klotz, J. L. & Heitmann, R. N. Effects of weaning and ionophore supplementation on selected blood metabolites and growth in dairy calves. J. Dairy Sci. 89, 3587–3598. https://doi.org/10.3168/jds.S0022-0302(06)72398-0 (2006).

Article CAS PubMed Google Scholar

Yuste, S. et al. Effect of preweaning diet on performance, blood metabolites and rumen fermentation around weaning in calves of two beef breeds. Anim. Prod. Sci. 60, 1018–1027 (2020).

Article CAS Google Scholar

Schwarzkopf, S. et al. Weaning Holstein calves at 17 weeks of age enables smooth transition from liquid to solid feed. Animals https://doi.org/10.3390/ani9121132 (2019).

Article PubMed PubMed Central Google Scholar

Gonzalez, L. A., Kyriazakis, I. & Tedeschi, L. O. Review: Precision nutrition of ruminants: Approaches, challenges and potential gains. Animal 12, s246–s261. https://doi.org/10.1017/S1751731118002288 (2018).

Article CAS PubMed Google Scholar

Download references

We thank Dr. Ian Luck for expert advice on 1H NMR and Mr. P. Lipscombe, Ms J. Lipscombe and Ms. S. Steedman for the management and care of animals. We also thank Dr Augusto Imaz for the preparation of the figures, Mr. Allan Lisle for exert statistical advice and Emeritus Professor Alan Bell for constructive comments on an early version of the paper. JGSC was supported by the CNPq Science Without Borders Program (Brazil) and BCK was supported by the CAPES Foundation (Ministry of Education of Brazil), Prof. Luciano A. Gonzalez was supported by the McCaughey Memorial Trust and the Nancy Roma Paech bequest. The study was partly funded by the Nancy Roma Paech Bequest to Prof. Michael J. D’Occhio.

The study was funded by the Nancy Roma Paech bequest, the McCaughey Memorial Trust, and The University of Sydney.

Sydney Institute of Agriculture, and School of Life and Environmental Sciences, Faculty of Science, The University of Sydney, Camden, NSW, 2570, Australia

Luciano A. González, Julia G. S. Carvalho, Bruno C. Kuinchtner & Michael J. D’Occhio

Departamento de Reprodução Animal, Faculdade de Medicina Veterinária e Zootecnia, Universidade de São Paulo, São Paulo, SP, Brazil

Julia G. S. Carvalho & Pietro S. Baruselli

Natural Pasture Ecology Laboratory (LEPAN), Universidade Federal de Santa Maria, Santa Maria, RS, Brazil

Bruno C. Kuinchtner

Kolling Institute of Medical Research, Northern Medical School, University of Sydney, St Leonards, NSW, 2065, Australia

Anthony C. Dona

You can also search for this author in PubMed Google Scholar

J.G.S.C.: conceptualization, data curation, formal analysis, investigation, methodology, writing—original draft; B.C.K., A.C.D., conceptualization, data curation, formal analysis, investigation, methodology, writing—review and editing; P.S.B.: conceptualization, writing—review and editing; M.J.D.: conceptualization, methodology, funding acquisition, investigation, methodology, project administration, resources, supervision, writing—review and editing; L.A.G., conceptualization, data curation, formal analysis, methodology, supervision, writing—original draft.

Correspondence to Luciano A. González.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

González, L.A., Carvalho, J.G.S., Kuinchtner, B.C. et al. Plasma metabolomics reveals major changes in carbohydrate, lipid, and protein metabolism of abruptly weaned beef calves. Sci Rep 13, 8176 (2023). https://doi.org/10.1038/s41598-023-35383-2

Download citation

Received: 07 April 2023

Accepted: 17 May 2023

Published: 20 May 2023

DOI: https://doi.org/10.1038/s41598-023-35383-2

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.