Responses of selected plasma metabolites to a two-day nutritional challenge of goats divergently selected for functional longevity

Understanding the extent to which genetics × environment plays a role in shaping individual strategies to environmental challenges is of considerable interest for future selection of more resilient animals. Accordingly, the objective of this study was to evaluate the metabolic responses to a nutritional challenge of goats divergently selected for functional longevity based on plasma metabolites and the repeatability of these responses across 2 experimental farms and years. We carried out 6 different experimental trials from years 2018 to 2022 (4 trials on site Bourges (2018-21) and 2 trials (2021-22) on site Grignon) in which 267 first kidding goats, daughters of Alpine bucks divergently selected for functional longevity, longevity plus (n = 137), and longevity minus (n = 130), were exposed to a 2-d nutritional challenge in early lactation. The experiments consisted of a 5 or 7-d control period (pre-challenge) on a standard lactation diet followed by a 2-d nutritional challenge with straw-only feeding and then a 7 or 10-d recovery period on a standard lactation diet, for site Bourges and Grignon, respectively. During the challenge plasma metabolite composition was recorded daily. Linear mixed-effects models were used to analyze all traits, considering the individual as a random effect and the 2x2 treatments (i.e., genetic line and year nested in site) and litter size as fixed effects. The linear mixed-effects model using a piecewise arrangement was used to analyze the response/recovery profiles to the nutritional challenge. Random parameters estimated for each individual, using the mixed-effects models without the fixed effects of genetic line, were used in a Sparse Partial Least Square Discriminant Analysis (sPLS-DA) to compare the goat metabolism response to the challenge on a multivariate scale. The plasma metabolites, glucose, β-hydroxybutyrate (BHB), and nonesterified fatty acids (NEFA), and urea concentrations responded to the 2-d nutritional challenge. Selection for functional longevity did not affect plasma glucose, NEFA, BHB, and urea response/recoveries to a 2-d nutritional challenge. However, site, trial, and litter size affected these responses. Moreover, the plasma metabolites seem not to fully recover to prechallenge levels after the recovery phase. The sPLS-DA analysis did not discriminate between the 2 longevity lines. We observed meaningful between-individuals' variability in plasma BHB, especially on the prechallenge and rate of response and rate of recovery from the 2-d nutritional challenge (CV = 26.2%, 36.1%, and 41.2%, repeatability = 0.749, 0.322, and 0.741, respectively). Plasma NEFA recovery from challenge also demonstrated high between-individuals' variability (CV = 16.4%, repeatability = 0.323). Selection for functional longevity did not affect plasma metabolites responses to a 2-d nutritional challenge in dairy goats. Plasma NEFA and BHB response/recovery presented high between-individuals' variability, indicating individual adaptative characteristics to nutritional challenges not related to the environmental conditions but to inherent individual characteristics.

Multivariate analysis of milk metabolite measures shows potential for deriving new resilience phenotypes

In a context of growing interest in breeding more resilient animals, a noninvasive indicator of resilience would be very valuable. We hypothesized that the time-course of concentrations of several milk metabolites through a short-term underfeeding challenge could reflect the variation of resilience mechanisms to such a challenge. We submitted 138 one-year-old primiparous goats, selected for extreme functional longevity (i.e., productive longevity corrected for milk yield [60 low longevity line goats and 78 high longevity line goats]), to a 2-d underfeeding challenge during early lactation. We measured the concentration of 13 milk metabolites and the activity of 1 enzyme during prechallenge, challenge, and recovery periods. Functional principal component analysis summarized the trends of milk metabolite concentration over time efficiently without preliminary assumptions concerning the shapes of the curves. We first ran a supervised prediction of the longevity line of the goats based on the milk metabolite curves. The partial least square analysis could not predict the longevity line accurately. We thus decided to explore the large overall variability of milk metabolite curves with an unsupervised clustering. The large year × facility effect on the metabolite concentrations was precorrected for. This resulted in 3 clusters of goats defined by different metabolic responses to underfeeding. The cluster that showed higher β-hydroxybutyrate, cholesterol, and triacylglycerols increase during the underfeeding challenge was associated with poorer survival compared with the other 2 clusters. These results suggest that multivariate analysis of noninvasive milk measures show potential for deriving new resilience phenotypes.

Milk metabolites can characterise individual differences in animal resilience to a nutritional challenge in lactating dairy goats

The aim of this study is built in two phases: to quantify the ability of novel milk metabolites to measure between-animal variability in response and recovery profiles to a short-term nutritional challenge, then to derive a resilience index from the relationship between these individual variations. At two different stages of lactation, sixteen lactating dairy goats were exposed to a 2-d underfeeding challenge. The first challenge was in late lactation, and the second was carried out on the same goats early in the following lactation. During the entire experiment period, samples were taken at each milking for milk metabolite measures. For each metabolite, the response profile of each goat was characterised using a piecewise model for describing the dynamic pattern of response and recovery profiles after the challenge relative to the start of the nutritional challenge. Cluster Analysis identified three types of response/recovery profiles per metabolite. Using cluster membership, multiple correspondence analyses (MCAs) were performed to further characterise response profile types across animals and metabolites. This MCA analysis identified three groups of animals. Further, discriminant path analysis was able to separate these groups of multivariate response/recovery profile type based on threshold levels of three milk metabolites: β-hydroxybutyrate, free glucose and uric acid. Further analyses were done to explore the possibility of developing an index of resilience from milk metabolite measures. Different types of performance response to short-term nutritional challenge can be distinguished using multivariate analyses of a panel of milk metabolites.