<figure><img src="image_1.png" /></figure>

HAL
open science

# A novel modelling approach to quantify the response of dairy goats to a high-concentrate diet

Masoomeh Taghipoor, Maud Delattre, Sylvie Giger-Reverdin

► To cite this version:

Masoomeh Taghipoor, Maud Delattre, Sylvie Giger-Reverdin. A novel modelling approach to quantify the response of dairy goats to a high-concentrate diet. Scientific Reports, 2020, 10 (1), pp.20376. 10.1038/s41598-020-77353-y. hal-03132771

HAL Id: hal-03132771

https://hal.inrae.fr/hal-03132771

Submitted on 24 Feb 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# A novel modelling approach to quantify the response of dairy goats to a high-concentrate diet

Masoomeh Taghipoor, Maud Delattre, Sylvie Giger-Reverdin

► To cite this version:

Masoomeh Taghipoor, Maud Delattre, Sylvie Giger-Reverdin. A novel modelling approach to quantify the response of dairy goats to a high-concentrate diet. Scientific Reports, 2020, 10 (1), pp.20376. 10.1038/s41598-020-77353-y. hal-03132771

HAL Id: hal-03132771

https://hal.inrae.fr/hal-03132771v1

Submitted on 24 Feb 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

OPEN

# A novel modelling approach to quantify the response of dairy goats to a high-concentrate diet

Masoomeh Taghipoor<sup>D1✉</sup>, Maud Delattre<sup>2</sup> & Sylvie Giger-Reverdin<sup>1</sup>

High-producing ruminants need high-concentrate diets to satisfy their nutrient requirements and meet performance objectives. However, such diets induce sub-acute ruminal acidosis (SARA), which will adversely affect dry matter intake and lead to lower production performance. This work develops a novel modelling approach to quantify the capacity of dairy goats to adapt to a high-concentrate diet challenge at the individual level. The animal model used was dairy goats (from Saanen or Alpine breed), and rumen pH was used as the indicator of the response. A three-step modelling procedure was developed to quantify daily scores and produce a single global index for animals' adaptive response to the new diet. The first step summarizes the post-prandial kinetics of rumen acid status using three synthetic variables. In the second step, the effect of time on the response of goats is described, in the short and long terms. In the last step, a metric based on phase trajectories ranks goats for their resilience capacity. This modelling procedure showed a high variability among the goats in response to the new diet, highlighting in particular their daily and general strategies to buffer the effect of the diet change. Two main categories of adaptive strategies were observed: (i) acid status increased, but the goats tried to minimize its variations, and (ii) acid status oscillated between increases and decreases. Such phenotyping, alongside other behavioral, digestive, and metabolic measures, can help to determine biomarkers of goats' capacity to adapt to diets of higher nutritive value and to increase production performance without compromising their health status. Quantifying the capacity of goats to buffer the effect of highly fermentable diets helps to better adapt feed to animals in precision livestock farming. This procedure is generic and can be adapted to any indicator of animal health and performance. In particular, several indicators can be combined to assess multi-performance, which is of major interest in the context of selection for robust animals.

Although ruminants can valorise fibrous low-quality diets<fnref n="1" />, high-producing ruminants need high-concentrate diets to satisfy their nutrient requirements. Such diets are often rich in highly and rapidly fermentable carbohydrates, and so increase the production of volatile fatty acids, which in turn decreases rumen pH to non-physiological levels². When rumen pH falls below a threshold often set at around 6.0, cellulolysis is inhibited and therefore dry matter intake is depressed³. Animals are then considered to be in a condition of sub-acute rumen acidosis (SARA), some physiological and pathological consequences of which have been extensively reviewed by Owens et al.⁴ The main effects of SARA are decreased feed efficiency and a modification of various production traits. Daily rumen acid status pattern varies considerably between diets and between animals fed the same diet⁵. It is well known that within a herd, some animals are susceptible to SARA when fed a high-concentrate diet, while others are tolerant⁶. Several indicators have been proposed to evaluate this pattern and the occurrence of SARA: initial and final pH values, amplitude of pH variations, area under a given pH threshold and the pH curve, time under a given threshold, etc.⁷ However, to our knowledge, no quantification of individuals' response in terms of rumen acid status has yet been carried out. Assigning an index to each individual is a step forward in the search for non-invasive biomarkers of the animal's capacity to buffer the effect of a high-concentrate diet. The comprehension of mechanisms underlying this capacity will be assisted by associating this global index with other behavioural and physiological measures in response to perturbation⁴,⁶,⁸. Ongoing environmental, economic and societal changes have meanwhile led animal scientists and selection companies to consider new strategies to select more robust animals⁹⁻¹¹. The index described here allows ranking of ruminants for their capacity to buffer the effect of a highly fermentable diet, and subsequently refining strategies to select for animals more resistant to

<footnote n="1">¹Université Paris-Saclay, INRAE, AgroParisTech, UMR Modélisation Systémique Appliquée aux Ruminants, 75005 Paris, France. ²Université Paris-Saclay, INRAE, MalAGE, 78350 Jouy-en-Josas, France. ³email: Masoomeh.Taghipoor@inrae.fr</footnote>
<figure><img src="image_2.png" /></figure>

OPEN

# A novel modelling approach to quantify the response of dairy goats to a high-concentrate diet

Masoomeh Taghipoor<sup>D1✉</sup>, Maud Delattre<sup>2</sup> & Sylvie Giger-Reverdin<sup>1</sup>

<figure><img src="image_3.png" /></figure>

High-producing ruminants need high-concentrate diets to satisfy their nutrient requirements and meet performance objectives. However, such diets induce sub-acute ruminal acidosis (SARA), which will adversely affect dry matter intake and lead to lower production performance. This work develops a novel modelling approach to quantify the capacity of dairy goats to adapt to a high-concentrate diet challenge at the individual level. The animal model used was dairy goats (from Saanen or Alpine breed), and rumen pH was used as the indicator of the response. A three-step modelling procedure was developed to quantify daily scores and produce a single global index for animals' adaptive response to the new diet. The first step summarizes the post-prandial kinetics of rumen acid status using three synthetic variables. In the second step, the effect of time on the response of goats is described, in the short and long terms. In the last step, a metric based on phase trajectories ranks goats for their resilience capacity. This modelling procedure showed a high variability among the goats in response to the new diet, highlighting in particular their daily and general strategies to buffer the effect of the diet change. Two main categories of adaptive strategies were observed: (i) acid status increased, but the goats tried to minimize its variations, and (ii) acid status oscillated between increases and decreases. Such phenotyping, alongside other behavioral, digestive, and metabolic measures, can help to determine biomarkers of goats' capacity to adapt to diets of higher nutritive value and to increase production performance without compromising their health status. Quantifying the capacity of goats to buffer the effect of highly fermentable diets helps to better adapt feed to animals in precision livestock farming. This procedure is generic and can be adapted to any indicator of animal health and performance. In particular, several indicators can be combined to assess multi-performance, which is of major interest in the context of selection for robust animals.

Although ruminants can valorise fibrous low-quality diets<fnref n="1" />, high-producing ruminants need high-concentrate diets to satisfy their nutrient requirements. Such diets are often rich in highly and rapidly fermentable carbohydrates, and so increase the production of volatile fatty acids, which in turn decreases rumen pH to non-physiological levels². When rumen pH falls below a threshold often set at around 6.0, cellulolysis is inhibited and therefore dry matter intake is depressed³. Animals are then considered to be in a condition of sub-acute rumen acidosis (SARA), some physiological and pathological consequences of which have been extensively reviewed by Owens et al.⁴ The main effects of SARA are decreased feed efficiency and a modification of various production traits. Daily rumen acid status pattern varies considerably between diets and between animals fed the same diet⁵. It is well known that within a herd, some animals are susceptible to SARA when fed a high-concentrate diet, while others are tolerant⁶. Several indicators have been proposed to evaluate this pattern and the occurrence of SARA: initial and final pH values, amplitude of pH variations, area under a given pH threshold and the pH curve, time under a given threshold, etc.⁷ However, to our knowledge, no quantification of individuals' response in terms of rumen acid status has yet been carried out. Assigning an index to each individual is a step forward in the search for non-invasive biomarkers of the animal's capacity to buffer the effect of a high-concentrate diet. The comprehension of mechanisms underlying this capacity will be assisted by associating this global index with other behavioural and physiological measures in response to perturbation⁴,⁶,⁸. Ongoing environmental, economic and societal changes have meanwhile led animal scientists and selection companies to consider new strategies to select more robust animals⁹⁻¹¹. The index described here allows ranking of ruminants for their capacity to buffer the effect of a highly fermentable diet, and subsequently refining strategies to select for animals more resistant to

<footnote n="1">¹Université Paris-Saclay, INRAE, AgroParisTech, UMR Modélisation Systémique Appliquée aux Ruminants, 75005 Paris, France. ²Université Paris-Saclay, INRAE, MalAGE, 78350 Jouy-en-Josas, France. ³email: Masoomeh.Taghipoor@inrae.fr</footnote>
SARA. Finally, in livestock precision farming, this index helps to adapt feed composition to an animal's rumen capacity and thereby efficiently valorise food composition¹².

The main objective of this study was to develop a modelling approach to quantify animals' capacity to regulate rumen acid status in response to a high-concentrate diet challenge, at the individual level. The animal model used was the dairy goat (from Alpine or Saanen breeds). The indicator of rumen acid status was pH, the post-prandial kinetics of which were measured throughout the 9 days of sampling. Various authors have quantified animals' response to different types of perturbations, such as feed restriction, weaning, response to mycotoxin challenges or general perturbations with unknown origins¹³⁻¹⁶. In these studies, the indicator of the animal response (body weight, feed intake, biochemical samples in plasma, milk yield, etc.) was measured at one given time scale (daily, weekly, etc.). The main difficulty in quantifying goats' response in terms of pH in the present trial was the presence of two different and imbricated sampling time scales, namely the hourly scale of post-prandial pH and the daily scale of the experimentation. To address this problem, an original three-step modelling procedure was developed. Firstly, the variations of the post-prandial acid status pattern were modelled to obtain some biologically relevant synthetic variables describing the post-prandial kinetics. Secondly, a mixed model was developed for each of the selected synthetic variables to study the short- and long-term effects of the new diet. It was hypothesised that goats would be able to regulate their rumen pH in the short term, whereas in the long term it might become more and more difficult to overcome the challenge of the high-concentrate diet. Finally, the variables with significant long-term effects were used to develop a modelling approach drawing on health trajectories¹⁷ to quantify individual goats' capacity to meet the challenge of a new diet. Using this approach, it was possible to assign a single numerical index to each goat describing its performance throughout the period of experimentation.

Since working with the small values of hydrogen ion concentration ([H⁺], moles per litre) was inconvenient, Sørensen¹⁸ advocated a pH scale, where pH = -log₁₀[H⁺]. Although this is a useful scale to measure rumen acidity, it may lead to errors in statistical tests and mathematical calculations¹⁹. In this work, values of [H⁺] were therefore used for model development and statistical tests. A comparison between results obtained using [H⁺] and pH is made in "Results" section.

## Materials and methods

The goats were cared for and handled in accordance with the French legislation on animal experimentation and in line with the European Convention for the Protection of Vertebrates Used for Experimental and Other Scientific Purposes (European Directive 86/609). All experimental procedures followed the guidelines for the care and use of experimental animals and were approved by the local ethics committee (Comité d'Ethique en Expérimentation Animale, COMETHEA 45, registered as 15-04). This study was conducted, from mid-April to mid-May 2015, at the experimental farm of the Inrae-AgroParisTech MoSAR Research Unit (Thiverval-Grignon, France, 48°51' N 1°55' E, 70 m a.s.l.).

**Experimental design.** Eight rumen-cannulated dairy goats (four Alpine and four Saanen) were habituated to a total mixed ration with a low concentrate level (20% on a dry matter basis). This was abruptly switched to a total mixed ration with a high concentrate level (50%). The two diets were formulated to be iso-protein (iso-PDI or truly digestible protein) according to the renewed French Inrae System for ruminants²⁰. Goats were fed in accordance with their nutritional requirements and in average 13% of refusal was observed (13% ± 8%). This means that goats were able to express their natural feeding behaviour. Moreover, the diet was a TMR (Total Mixed Ratio) to prevent from sorting.

Rumen fluid was sampled through the cannula before and 1, 2, 4 and 6 h after the morning feed delivery on 2 days before the change of diet, *d*₁ and *d*₂ (where *d*₂ = *d*₁ + 3), then on four consecutive days following the diet change, *d*₃, *d*₄, *d*₅ and *d*₆, and finally once a week for 3 weeks, *d*₇, *d*₈ and *d*₉, where *d*₇ = *d*₆ + 4, *d*₈ = *d*₇ + 6, and *d*₉ = *d*₈ + 7. Rumen pH was measured immediately after sampling. Figure 1 illustrates the experimental design.

To allow for the time interval between milking, feed was delivered shortly after milking (7 a.m. and 3 p.m.), with 1/3 of the total ration in the morning and 2/3 in the afternoon. Diet samples were collected for 9 days corresponding to rumen sampling. They were analysed separately for dry matter (DM)²¹ and cell wall components estimated by the neutral detergent fibre method (NDF) of Van Soest and Wine²² modified by Giger et al.²³ These last authors suggest using a heat-stable alpha-amylase without sodium sulphite and decalin as proposed by Robertson and Van Soest²⁴. Lignocellulose (ADF) and acid detergent lignin were obtained using a sequential approach on the NDF residue²³. Starch content was analysed according to AFNOR (1997) method²⁵. Total N was determined by the Dumas technique²⁶, and crude protein was estimated as 6.25 N. The nutritive value of the diets was obtained according to the renewed French INRA System for ruminants²⁰. Details are provided in supplementary Table S1.

**Modelling procedure.** Values of pH were first transformed by applying the logarithmic function to obtain the values of hydrogen ion concentration [H⁺](pH = -log₁₀[H⁺]).

The main difficulty in developing a model to describe such data was the presence of two different imbricated time scales, namely several weeks of experimentation versus the post-prandial kinetics scale (hours, Fig. 1). Accordingly, the post-prandial kinetics of [H⁺] were first summarised with small numbers of synthetic variables, which were then integrated into the weekly data analysis.

The modelling procedure is detailed in three steps. The first step describes the procedure to determine the synthetic variables, the second step explains the development of the mixed model to study the effect of time on goats' response, and in the third step, an original method is introduced to quantify animals' capacity to adapt to the new diet.
SARA. Finally, in livestock precision farming, this index helps to adapt feed composition to an animal's rumen capacity and thereby efficiently valorise food composition¹².

The main objective of this study was to develop a modelling approach to quantify animals' capacity to regulate rumen acid status in response to a high-concentrate diet challenge, at the individual level. The animal model used was the dairy goat (from Alpine or Saanen breeds). The indicator of rumen acid status was pH, the post-prandial kinetics of which were measured throughout the 9 days of sampling. Various authors have quantified animals' response to different types of perturbations, such as feed restriction, weaning, response to mycotoxin challenges or general perturbations with unknown origins¹³⁻¹⁶. In these studies, the indicator of the animal response (body weight, feed intake, biochemical samples in plasma, milk yield, etc.) was measured at one given time scale (daily, weekly, etc.). The main difficulty in quantifying goats' response in terms of pH in the present trial was the presence of two different and imbricated sampling time scales, namely the hourly scale of post-prandial pH and the daily scale of the experimentation. To address this problem, an original three-step modelling procedure was developed. Firstly, the variations of the post-prandial acid status pattern were modelled to obtain some biologically relevant synthetic variables describing the post-prandial kinetics. Secondly, a mixed model was developed for each of the selected synthetic variables to study the short- and long-term effects of the new diet. It was hypothesised that goats would be able to regulate their rumen pH in the short term, whereas in the long term it might become more and more difficult to overcome the challenge of the high-concentrate diet. Finally, the variables with significant long-term effects were used to develop a modelling approach drawing on health trajectories¹⁷ to quantify individual goats' capacity to meet the challenge of a new diet. Using this approach, it was possible to assign a single numerical index to each goat describing its performance throughout the period of experimentation.

Since working with the small values of hydrogen ion concentration ([H⁺], moles per litre) was inconvenient, Sørensen¹⁸ advocated a pH scale, where pH = -log₁₀[H⁺]. Although this is a useful scale to measure rumen acidity, it may lead to errors in statistical tests and mathematical calculations¹⁹. In this work, values of [H⁺] were therefore used for model development and statistical tests. A comparison between results obtained using [H⁺] and pH is made in "Results" section.

## Materials and methods

The goats were cared for and handled in accordance with the French legislation on animal experimentation and in line with the European Convention for the Protection of Vertebrates Used for Experimental and Other Scientific Purposes (European Directive 86/609). All experimental procedures followed the guidelines for the care and use of experimental animals and were approved by the local ethics committee (Comité d'Ethique en Expérimentation Animale, COMETHEA 45, registered as 15-04). This study was conducted, from mid-April to mid-May 2015, at the experimental farm of the Inrae-AgroParisTech MoSAR Research Unit (Thiverval-Grignon, France, 48°51' N 1°55' E, 70 m a.s.l.).

**Experimental design.** Eight rumen-cannulated dairy goats (four Alpine and four Saanen) were habituated to a total mixed ration with a low concentrate level (20% on a dry matter basis). This was abruptly switched to a total mixed ration with a high concentrate level (50%). The two diets were formulated to be iso-protein (iso-PDI or truly digestible protein) according to the renewed French Inrae System for ruminants²⁰. Goats were fed in accordance with their nutritional requirements and in average 13% of refusal was observed (13% ± 8%). This means that goats were able to express their natural feeding behaviour. Moreover, the diet was a TMR (Total Mixed Ratio) to prevent from sorting.

Rumen fluid was sampled through the cannula before and 1, 2, 4 and 6 h after the morning feed delivery on 2 days before the change of diet, *d*₁ and *d*₂ (where *d*₂ = *d*₁ + 3), then on four consecutive days following the diet change, *d*₃, *d*₄, *d*₅ and *d*₆, and finally once a week for 3 weeks, *d*₇, *d*₈ and *d*₉, where *d*₇ = *d*₀ + 4, *d*₈ = *d*₇ + 6, and *d*₉ = *d*₈ + 7. Rumen pH was measured immediately after sampling. Figure 1 illustrates the experimental design.

To allow for the time interval between milking, feed was delivered shortly after milking (7 a.m. and 3 p.m.), with 1/3 of the total ration in the morning and 2/3 in the afternoon. Diet samples were collected for 9 days corresponding to rumen sampling. They were analysed separately for dry matter (DM)²¹ and cell wall components estimated by the neutral detergent fibre method (NDF) of Van Soest and Wine²² modified by Giger et al.²³ These last authors suggest using a heat-stable alpha-amylase without sodium sulphite and decalin as proposed by Robertson and Van Soest²⁴. Lignocellulose (ADF) and acid detergent lignin were obtained using a sequential approach on the NDF residue²³. Starch content was analysed according to AFNOR (1997) method²⁵. Total N was determined by the Dumas technique²⁶, and crude protein was estimated as 6.25 N. The nutritive value of the diets was obtained according to the renewed French INRA System for ruminants²⁰. Details are provided in supplementary Table S1.

**Modelling procedure.** Values of pH were first transformed by applying the logarithmic function to obtain the values of hydrogen ion concentration [H⁺](pH = -log₁₀[H⁺]).

The main difficulty in developing a model to describe such data was the presence of two different imbricated time scales, namely several weeks of experimentation versus the post-prandial kinetics scale (hours, Fig. 1). Accordingly, the post-prandial kinetics of [H⁺] were first summarised with small numbers of synthetic variables, which were then integrated into the weekly data analysis.

The modelling procedure is detailed in three steps. The first step describes the procedure to determine the synthetic variables, the second step explains the development of the mixed model to study the effect of time on goats' response, and in the third step, an original method is introduced to quantify animals' capacity to adapt to the new diet.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Experimental design. Goats were fed the standard diet in the first week (the first 2 days). Their diet was changed to the high-concentrate diet in the next 4 weeks of experimentation (Days 3 to 9). At each sampling day, pH was recorded before and 1, 2, 4 and 6 h after feed delivery. The presence of two embedded dynamics can be seen in the scheme: one for post-prandial kinetics of pH, the second for the sampling days.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Illustration of the parameters of the quadratic function used to fit post-prandial [H⁺] curves. Parameter v₀ is the initial value of [H⁺], A is the maximum deviation from the initial value and represents the difference between max[H⁺] and v₀. Parameter R is not illustrated in the figure and is the percentage of recovery after 6 h. θ is the threshold of acidity.</figcaption></figure>

**Step 1: Synthetic variables to describe the post-prandial kinetics of [H⁺].** To determine the synthetic variables, the post-prandial kinetics of [H⁺] were described with a quadratic function,

$$[H^+](t) = at^2 + bt + c, \quad \text{Model (1.1)}$$

where t is time of sampling and a, b and c are the quadratic function coefficients. This function was re-parametrised to obtain a function with biologically meaningful parameters²⁷,

$$[H^+](t) = f(t, v_0, A, R), \quad \text{Model (1.2)}$$

where v₀, A and R are the new parameters of the model, v₀ is the initial value of [H⁺], A is the amplitude of acidosis or the extent of the deviation from v₀, and R represents the recovery capacity of the animal 6 h after feed delivery (see supplementary Eq. 1 for the re-parametrisation).

The threshold of acidosis θ is the value of [H⁺] above which the animal is assumed to be in a condition of rumen acidosis (Fig. 2). This model was then fitted to all the individual post-prandial kinetics of [H⁺] to estimate the triplet v₀, A, R. Besides obtaining biologically meaningful parameters, the re-parametrisation of the quadratic function helped to avoid the statistical difficulties arising from the initiation of the model parameters (a, b, c) for 9 days of post-prandial kinetics of eight goats.

Besides v₀, A, R, other synthetic variables can be calculated from Model 1.2. The variable duration (dur) stands for the duration of acidosis when [H⁺] exceeds θ, amplitude of acidosis (AmpAc) is associated with the difference between the maximum of [H⁺] and θ, and the variable varlast is the estimation of [H⁺] 6 h after the feed delivery. Variables AmpAc and dur are dependent on the values of θ.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Experimental design. Goats were fed the standard diet in the first week (the first 2 days). Their diet was changed to the high-concentrate diet in the next 4 weeks of experimentation (Days 3 to 9). At each sampling day, pH was recorded before and 1, 2, 4 and 6 h after feed delivery. The presence of two embedded dynamics can be seen in the scheme: one for post-prandial kinetics of pH, the second for the sampling days.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Illustration of the parameters of the quadratic function used to fit post-prandial [H⁺] curves. Parameter v₀ is the initial value of [H⁺], A is the maximum deviation from the initial value and represents the difference between max[H⁺] and v₀. Parameter R is not illustrated in the figure and is the percentage of recovery after 6 h. θ is the threshold of acidity.</figcaption></figure>

**Step 1: Synthetic variables to describe the post-prandial kinetics of [H⁺].** To determine the synthetic variables, the post-prandial kinetics of [H⁺] were described with a quadratic function,

$$[H^+](t) = at^2 + bt + c, \quad \text{Model (1.1)}$$

where t is time of sampling and a, b and c are the quadratic function coefficients. This function was re-parametrised to obtain a function with biologically meaningful parameters²⁷,

$$[H^+](t) = f(t, v_0, A, R), \quad \text{Model (1.2)}$$

where v₀, A and R are the new parameters of the model, v₀ is the initial value of [H⁺], A is the amplitude of acidosis or the extent of the deviation from v₀, and R represents the recovery capacity of the animal 6 h after feed delivery (see supplementary Eq. 1 for the re-parametrisation).

The threshold of acidosis θ is the value of [H⁺] above which the animal is assumed to be in a condition of rumen acidosis (Fig. 2). This model was then fitted to all the individual post-prandial kinetics of [H⁺] to estimate the triplet v₀, A, R. Besides obtaining biologically meaningful parameters, the re-parametrisation of the quadratic function helped to avoid the statistical difficulties arising from the initiation of the model parameters (a, b, c) for 9 days of post-prandial kinetics of eight goats.

Besides v₀, A, R, other synthetic variables can be calculated from Model 1.2. The variable duration (dur) stands for the duration of acidosis when [H⁺] exceeds θ, amplitude of acidosis (AmpAc) is associated with the difference between the maximum of [H⁺] and θ, and the variable varlast is the estimation of [H⁺] 6 h after the feed delivery. Variables AmpAc and dur are dependent on the values of θ.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Illustration of Model 2. “Week” is a factor and varies from w1 to w5, “days” is a numerical variable associated with the days of sampling. As seen in the graph, $\alpha_j$ determines the intercept associated with each week and $\beta$ is the slope of daily variation of the synthetic variable during the first week after the diet change. Slope $\beta$ is applied only during week 2 because of the use of the indicator function.</figcaption></figure>

**Step 2: Statistical analysis on synthetic variables.** To study the effect of time on the goats’ response to the dietary change during the 5 weeks of experimentation, a mixed model was developed for each of the synthetic variables (Model 2), with individual goats as random effect. Fixed effects were the qualitative variable “week” for the weeks of sampling, associated with the response of the animal in the long term. The numerical variable “days”, which represented days of sampling, was introduced to describe the short-term response of goats during the first week. The variable “days” was taken to be positive after the delivery of high-concentrate diet (Day 1 was the first day of delivery of the high-concentrate diet). The random effect “goats” helped to quantify the contribution of individual variability among the goats to the total variance of the model$^{28}$.

$$v_{ijk} = \alpha_j + G_i + \beta \cdot (\text{days}_{k,1j=2}) + \epsilon_{ijk} \quad (\text{Model 2})$$

$$G_i \in \mathcal{N}(0, \sigma_A^2), \quad i = 1, \dots, 8$$

$$\epsilon_{ijk} \in \mathcal{N}(0, \sigma^2), \quad j = 1, \dots, 5, \quad k = 1, \dots, 9.$$

where i stands for the number of animals, j the number of weeks, and k the number of days of sampling, $v_{ijk}$ is the synthetic variable, $G_i$ is the random intercept for goat i, $\alpha_j$ is the fixed effect of "week" and $\beta$ is the fixed effect of "days". The indicator function ($1_{j=2}$) indicates that the effect of $\beta$ is applied only during the second week of sampling (first week after the delivery of high-concentrate diet),

$$1_{j=2} = \begin{cases} 1, & j = 2 \\ 0, & j \neq 2 \end{cases}$$

The model has a random intercept for the difference among goats at the beginning of experimentation. Figure 3 illustrates the fixed effect of the model.

**Step 3: Quantification of goats’ response.** To quantify individual capacity to adapt to the new diet, a method based on health trajectories$^{17,29}$ was developed and applied on synthetic variables with a significant long-term effect. In this method, for a given animal, the phase space for synthetic variables was plotted. Let us consider the example of the two-dimensional phase space (phase plane) for variables $v_0$ and A (Fig. 4). Based on daily variations of both synthetic variables under study for $i \in \{1, \dots, 8\}$, a daily score $s_{i \to i+1}$ from -2 to 2 is associated with an animal’s response from day $d_i$ to $d_{i+1}$. For both variables, a decrease is associated with an expected adaptive behaviour, and conversely, an increase is associated with a non-capacity to adapt (see Table 1 for the definition of each synthetic variable). Therefore, if both variables decrease, the score is 2, and if both variables increase, the score is -2. An increase in $v_0$ and a decrease in A is associated with the effort of the animal to adapt, i.e., despite the increase in the initial value $v_0$, the maximum has not changed or is decreased, which results in a decrease in A: therefore, to record the effort of adaptation, the score is 1. Finally, the score of -1 is assigned if despite the decrease in $v_0$, the value of A increases (maximum post-prandial $H^+$ increases). To consider the variations of synthetic variables and the number of days between two sampling days $d_i$ and $d_{i+1}$, the relative Euclidean distance in the phase plane was calculated and used as a weighting associated with daily scores (Model 3).

$$\omega_{i \to i+1} = \frac{1}{d_{i+1} - d_i} \cdot \sqrt{(A(d_{i+1}) - A(d_i))^2 + (v_0(d_{i+1}) - v_0(d_i))^2} \quad (\text{Model 3})$$

Using this method, a daily index and a global index of animal adaptive capacity are defined: daily index is the daily score associated with each animal, and global index is calculated as the sum of daily scores from day $d_2$ to $d_9$ (first day of diet change is $d_3$),

$$d_{i \to i+1} = s_{i \to i+1} \cdot \omega_{i \to i+1}$$

$$GI = \sum_{2}^{8} d_{i \to i+1},$$
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Illustration of Model 2. “Week” is a factor and varies from w1 to w5, “days” is a numerical variable associated with the days of sampling. As seen in the graph, $\alpha_j$ determines the intercept associated with each week and $\beta$ is the slope of daily variation of the synthetic variable during the first week after the diet change. Slope $\beta$ is applied only during week 2 because of the use of the indicator function.</figcaption></figure>

**Step 2: Statistical analysis on synthetic variables.** To study the effect of time on the goats' response to the dietary change during the 5 weeks of experimentation, a mixed model was developed for each of the synthetic variables (Model 2), with individual goats as random effect. Fixed effects were the qualitative variable “week” for the weeks of sampling, associated with the response of the animal in the long term. The numerical variable “days”, which represented days of sampling, was introduced to describe the short-term response of goats during the first week. The variable “days” was taken to be positive after the delivery of high-concentrate diet (Day 1 was the first day of delivery of the high-concentrate diet). The random effect “goats” helped to quantify the contribution of individual variability among the goats to the total variance of the model$^{28}$.

$$v_{ijk} = \alpha_j + G_i + \beta \cdot (\text{days}_{k,1j=2}) + \epsilon_{ijk} \quad (\text{Model 2})$$

$$G_i \in \mathcal{N}(0, \sigma_A^2), \quad i = 1, \dots, 8$$

$$\epsilon_{ijk} \in \mathcal{N}(0, \sigma^2), \quad j = 1, \dots, 5, \quad k = 1, \dots, 9.$$

where $i$ stands for the number of animals, $j$ the number of weeks, and $k$ the number of days of sampling, $v_{ijk}$ is the synthetic variable, $G_i$ is the random intercept for goat $i$, $\alpha_j$ is the fixed effect of "week" and $\beta$ is the fixed effect of "days". The indicator function $(1_{j=2})$ indicates that the effect of $\beta$ is applied only during the second week of sampling (first week after the delivery of high-concentrate diet),

$$1_{j=2} = \begin{cases} 1, & j = 2 \\ 0, & j \neq 2 \end{cases}$$

The model has a random intercept for the difference among goats at the beginning of experimentation. Figure 3 illustrates the fixed effect of the model.

**Step 3: Quantification of goats' response.** To quantify individual capacity to adapt to the new diet, a method based on health trajectories$^{17,29}$ was developed and applied on synthetic variables with a significant long-term effect. In this method, for a given animal, the phase space for synthetic variables was plotted. Let us consider the example of the two-dimensional phase space (phase plane) for variables $v_0$ and $A$ (Fig. 4). Based on daily variations of both synthetic variables under study for $i \in \{1, \dots, 8\}$, a daily score $s_{i \to i+1}$ from -2 to 2 is associated with an animal's response from day $d_i$ to $d_{i+1}$. For both variables, a decrease is associated with an expected adaptive behaviour, and conversely, an increase is associated with a non-capacity to adapt (see Table 1 for the definition of each synthetic variable). Therefore, if both variables decrease, the score is 2, and if both variables increase, the score is -2. An increase in $v_0$ and a decrease in $A$ is associated with the effort of the animal to adapt, i.e., despite the increase in the initial value $v_0$, the maximum has not changed or is decreased, which results in a decrease in $A$: therefore, to record the effort of adaptation, the score is 1. Finally, the score of -1 is assigned if despite the decrease in $v_0$, the value of $A$ increases (maximum post-prandial $H^+$ increases). To consider the variations of synthetic variables and the number of days between two sampling days $d_i$ and $d_{i+1}$, the relative Euclidean distance in the phase plane was calculated and used as a weighting associated with daily scores (Model 3).

$$\omega_{i \to i+1} = \frac{1}{d_{i+1} - d_i} \cdot \sqrt{(A(d_{i+1}) - A(d_i))^2 + (v_0(d_{i+1}) - v_0(d_i))^2} \quad (\text{Model 3})$$

Using this method, a daily index and a global index of animal adaptive capacity are defined: daily index is the daily score associated with each animal, and global index is calculated as the sum of daily scores from day $d_2$ to $d_9$ (first day of diet change is $d_3$),

$$d_{i \to i+1} = s_{i \to i+1} \cdot \omega_{i \to i+1}$$

$$GI = \sum_{2}^{8} d_{i \to i+1},$$
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Example of a phase plane of $v_0$ and A, representing variations of A as a function of $v_0$, rescaled by $10^7$ for clarity. Numbers (1 to 9) on the curve are days of sampling. A score $s_{i \to i+1}$ is assigned to the variation from day $d_i$ to day $d_{i+1}$, for $i \in \{1, \dots, 8\}$. The small plane on the right-hand side of the graph illustrates the assignment of scores. The score of -2 is assigned to $v_0^+ A^+$ describing an increase in both variables, a score of -1 is assigned to $v_0^+ A^+$, i.e. decreasing $v_0$ and increasing A, a score of 1 is associated with $v_0^+ A^-$ describing increasing $v_0$ and decreasing A, and finally 2 for the decrease in both variables. Each score is then weighted by $\omega_{i \to i+1}$ the relative Euclidian distance between days $d_i$ and $d_{i+1}$, calculated as the Euclidian distance between the 2 days divided by $(d_{i+1} - d_i)$. For this curve the set of scores can be written $\{-1 \omega_{1 \to 2}, -1 \omega_{2 \to 3}, -1 \omega_{3 \to 4}, 1 \omega_{4 \to 5}, -1 \omega_{5 \to 6}, 2 \omega_{6 \to 7}, -1 \omega_{7 \to 8}, 1 \omega_{8 \to 9}\}$. The adaptive capacity index associated with this animal can then be defined as the sum of the scores of the 9 days of sampling.</figcaption></figure>

## Group 1. parameters of the model

<table><thead><tr><th>Synthetic variables</th><th>Initial values</th><th>Definition</th></tr></thead><tbody><tr><td>$v_0$ (moles/litre)</td><td>$[H^+](t_0)$</td><td>Initial value of $[H^+]$ as estimated by the model</td></tr><tr><td>A (moles/litre)</td><td>$\max([H^+]) - v_0$</td><td>Difference between maximum of $[H^+_{estim}]$ and $v_0$</td></tr><tr><td>R (dimensionless)</td><td>$\frac{\max([H^+])-v_0}{\max([H^+])-v_{last}} \times 100$</td><td>The capacity of a goat to reach its initial $[H^+_{estim}]$ at $t=6$ h</td></tr></tbody></table>

## Group 2. calculated from the quadratic function for $\theta = 5.5$

<table><caption>**Table 1.** Two groups of candidate synthetic variables to characterise the post-prandial kinetics of $[H^+]$. $[H^+_{estim}]$ is the estimation of $[H^+]$ by the model. Variable dur is the time interval of acidosis where $[H^+_{estim}] < \theta$, and AmpAc is the extent of the deviation from the threshold $\theta = 5.5$. Variable $v_{last}$ is the value of $H^+$ at $t=6$ h, as estimated by the model. In contrast to Group 2, the synthetic variables of Group 1 are independent of the value of $\theta$.</caption><thead><tr><th>Synthetic variables</th><th>Formula</th><th>Definition</th></tr></thead><tbody><tr><td>dur (hours)</td><td>Time ($[H^+_{estim}] < \theta$)</td><td>Time interval of acidosis measured by the time that $[H^+_{estim}]$ is under the threshold $\theta$ of acidity</td></tr><tr><td>AmpAc (moles/litre)</td><td>($\max([H^+_{estim}]) - \theta$)</td><td>Amplitude of acidity, the difference between $\theta$ and $\max([H^+_{estim}]$</td></tr><tr><td>$v_{last}$</td><td>$[H^+](t_{last})$</td><td>The value of $[H^+]$ at $t=6$ h as estimated by the model</td></tr></tbody></table>

where, $d_{i \to i+1}$ stands for the daily score assigned to the daily response of a goat, and GI is the global index for each goat.

**Statistical analysis.** Software R version 3.5.3 (R core team 2019; https://www.R-project.org/)³⁰ was used for all statistical analyses. For correlation analysis, because the distribution of synthetic variables was not Gaussian, the Spearman correlation test was used. To study the effect of fixed variables of the model, the likelihood ratio test (ANOVA function in R) was used. When significant effects of weeks were found, a post hoc comparison test with FDR correction was applied (emeans package³¹). The function nls of package stat was used to estimate the parameters of the quadratic function (Model 1.2). For each of the synthetic variables, the Gaussian distribution was checked with a Shapiro-Wilk test³². In the absence of normality, a box-cox transformation was applied to determine the adequate power transformation, using the boxcox function of R. Finally, the function lme4 was used for the mixed models³³.

## Results

Figure 5 shows the postprandial kinetics of $[H^+]$ along the 9 sampling days of the trial, for all the goats. No apparent health problem was observed, even though the range of rumen pH observed varied from 5.29 to 7.07 (see Supplementary Fig. S1 online).
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Example of a phase plane of $v_0$ and A, representing variations of A as a function of $v_0$, rescaled by $10^7$ for clarity. Numbers (1 to 9) on the curve are days of sampling. A score $s_{i \to i+1}$ is assigned to the variation from day $d_i$ to day $d_{i+1}$, for $i \in \{1, \dots, 8\}$. The small plane on the right-hand side of the graph illustrates the assignment of scores. The score of -2 is assigned to $v_0^+ A^+$ describing an increase in both variables, a score of -1 is assigned to $v_0^+ A^+$, i.e. decreasing $v_0$ and increasing A, a score of 1 is associated with $v_0^+ A^-$ describing increasing $v_0$ and decreasing A, and finally 2 for the decrease in both variables. Each score is then weighted by $\omega_{i \to i+1}$ the relative Euclidian distance between days $d_i$ and $d_{i+1}$, calculated as the Euclidian distance between the 2 days divided by $(d_{i+1} - d_i)$. For this curve the set of scores can be written $\{-1 \omega_{1 \to 2}, -1 \omega_{2 \to 3}, -1 \omega_{3 \to 4}, 1 \omega_{4 \to 5}, -1 \omega_{5 \to 6}, 2 \omega_{6 \to 7}, -1 \omega_{7 \to 8}, 1 \omega_{8 \to 9}\}$. The adaptive capacity index associated with this animal can then be defined as the sum of the scores of the 9 days of sampling.</figcaption></figure>

<table><caption>**Table 1.** Two groups of candidate synthetic variables to characterise the post-prandial kinetics of $[H^+]$. $[H^+_{estim}]$ is the estimation of $[H^+]$ by the model. Variable dur is the time interval of acidosis where $[H^+_{estim}] < \theta$, and AmpAc is the extent of the deviation from the threshold $\theta = 5.5$. Variable $v_{last}$ is the value of $H^+$ at $t=6$ h, as estimated by the model. In contrast to Group 2, the synthetic variables of Group 1 are independent of the value of $\theta$.</caption><thead><tr><th colspan="3">Group 1. parameters of the model</th></tr><tr><th>Synthetic variables</th><th>Initial values</th><th>Definition</th></tr></thead><tbody><tr><td>$v_0$ (moles/litre)</td><td>$[H^+](t_0)$</td><td>Initial value of $[H^+]$ as estimated by the model</td></tr><tr><td>A (moles/litre)</td><td>$\max([H^+]) - v_0$</td><td>Difference between maximum of $[H^+_{estim}]$ and $v_0$</td></tr><tr><td>R (dimensionless)</td><td>$\frac{\max([H^+])-v_0}{\max([H^+])-v_{last}} \times 100$</td><td>The capacity of a goat to reach its initial $[H^+_{estim}]$ at $t=6$ h</td></tr></tbody></table><table><thead><tr><th colspan="3">Group 2. calculated from the quadratic function for $\theta = 5.5$</th></tr><tr><th>Synthetic variables</th><th>Formula</th><th>Definition</th></tr></thead><tbody><tr><td>dur (hours)</td><td>Time ($[H^+_{estim}] < \theta$)</td><td>Time interval of acidosis measured by the time that $[H^+_{estim}]$ is under the threshold $\theta$ of acidity</td></tr><tr><td>AmpAc (moles/litre)</td><td>($\max([H^+_{estim}]) - \theta$)</td><td>Amplitude of acidity, the difference between $\theta$ and $\max([H^+_{estim}]$</td></tr><tr><td>$v_{last}$</td><td>$[H^+](t_{last})$</td><td>The value of $[H^+]$ at $t=6$ h as estimated by the model</td></tr></tbody></table>

where, $d_{i \to i+1}$ stands for the daily score assigned to the daily response of a goat, and GI is the global index for each goat.

**Statistical analysis.** Software R version 3.5.3 (R core team 2019; https://www.R-project.org/)³⁰ was used for all statistical analyses. For correlation analysis, because the distribution of synthetic variables was not Gaussian, the Spearman correlation test was used. To study the effect of fixed variables of the model, the likelihood ratio test (ANOVA function in R) was used. When significant effects of weeks were found, a post hoc comparison test with FDR correction was applied (emeans package³¹). The function nls of package stat was used to estimate the parameters of the quadratic function (Model 1.2). For each of the synthetic variables, the Gaussian distribution was checked with a Shapiro-Wilk test³². In the absence of normality, a box-cox transformation was applied to determine the adequate power transformation, using the boxcox function of R. Finally, the function lme4 was used for the mixed models³³.

## Results

Figure 5 shows the postprandial kinetics of $[H^+]$ along the 9 sampling days of the trial, for all the goats. No apparent health problem was observed, even though the range of rumen pH observed varied from 5.29 to 7.07 (see Supplementary Fig. S1 online).
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Individuals' daily variations of [H⁺]. Each panel represents the post-prandial kinetics of individual goats along the 9 days of sampling. Different colours are associated with different time points of post-prandial kinetics of [H⁺]. For clarity, values of [H⁺] are rescaled by 10⁷.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Graph of correlation among synthetic variables. Correlation values and their significance are noted. Symbols *, ** and *** mean *p* values smaller than 0.05, 0.01, and 0.001, respectively. Because the distribution of synthetic variables was not Gaussian, the Spearman correlation test was used. *a*, *b* and *c* are the coefficients of the original quadratic function. *A*, *v*₀ and *R* are parameters of the new re-parametrised function. AmpAc and dur are parameters calculated after fitting the model to individual [H⁺] daily kinetics.</figcaption></figure>

**Step 1: Synthetic variables.** Figure 6 shows the correlation among different synthetic variables introduced in Table 1 (see supplementary Table S1 online for descriptive statistics of the synthetic variables). There was a high positive correlation between *A* and AmpAc; both indicators of the maximum post-prandial [H⁺] deviation (*r* = 0.89, *n* = 72). The variable *a* represents the curvature (second derivative) of the quadratic function, and was highly correlated with AmpAc and *A*, explaining the fact that the larger *a*, the greater the deviation of the function from the origin (with respectively, *r* = -0.90 and *r* = -0.96, *n* = 72). Different thresholds of acidity θ are reported in the literature³⁴, and the evolution of the synthetic variables should be considered in relation to θ. For example, AmpAc and dur are dependent on the values of θ. Increasing θ will down-shift both AmpAc and dur. This was expected, given the definitions of these variables in Table 1. Instead of using θ-dependent variables, three synthetic variables *v*₀, *A* and *R* were used to describe the post-prandial [H⁺] kinetics for each goat.

**Step 2: Statistical analysis on synthetic variables.** The Shapiro test showed that none of the three variables *v*₀, *A* and *R* were Gaussian (*p* < 0.05). A power transformation was applied (box-cox transformation),
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Individuals' daily variations of [H⁺]. Each panel represents the post-prandial kinetics of individual goats along the 9 days of sampling. Different colours are associated with different time points of post-prandial kinetics of [H⁺]. For clarity, values of [H⁺] are rescaled by 10⁷.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Graph of correlation among synthetic variables. Correlation values and their significance are noted. Symbols *, ** and *** mean *p* values smaller than 0.05, 0.01, and 0.001, respectively. Because the distribution of synthetic variables was not Gaussian, the Spearman correlation test was used. *a*, *b* and *c* are the coefficients of the original quadratic function. *A*, *v*₀ and *R* are parameters of the new re-parametrised function. AmpAc and dur are parameters calculated after fitting the model to individual [H⁺] daily kinetics.</figcaption></figure>

**Step 1: Synthetic variables.** Figure 6 shows the correlation among different synthetic variables introduced in Table 1 (see supplementary Table S1 online for descriptive statistics of the synthetic variables). There was a high positive correlation between *A* and AmpAc; both indicators of the maximum post-prandial [H⁺] deviation (*r* = 0.89, *n* = 72). The variable *a* represents the curvature (second derivative) of the quadratic function, and was highly correlated with AmpAc and *A*, explaining the fact that the larger *a*, the greater the deviation of the function from the origin (with respectively, *r* = -0.90 and *r* = -0.96, *n* = 72). Different thresholds of acidity θ are reported in the literature³⁴, and the evolution of the synthetic variables should be considered in relation to θ. For example, AmpAc and dur are dependent on the values of θ. Increasing θ will down-shift both AmpAc and dur. This was expected, given the definitions of these variables in Table 1. Instead of using θ-dependent variables, three synthetic variables *v*₀, *A* and *R* were used to describe the post-prandial [H⁺] kinetics for each goat.

**Step 2: Statistical analysis on synthetic variables.** The Shapiro test showed that none of the three variables *v*₀, *A* and *R* were Gaussian (*p* < 0.05). A power transformation was applied (box-cox transformation),
<table><thead><tr><th rowspan="2"></th><th colspan="5">Weeks</th><th colspan="2">p-values</th></tr><tr><th>W1</th><th>W2</th><th>W3</th><th>W4</th><th>W5</th><th>Week</th><th>Days</th></tr></thead><tbody><tr><th>v<sub>0</sub><sup>λ<sub>1</sub></sup></th><td>0.052<sup>a</sup></td><td>0.046<sup>c</sup></td><td>0.049<sup>bc</sup></td><td>0.047<sup>c</sup></td><td>0.054<sup>b</sup></td><td>&lt; 0.001</td><td>0.11</td></tr><tr><th>A<sup>λ<sub>2</sub></sup></th><td>0.016<sup>a</sup></td><td>0.017<sup>ab</sup></td><td>0.018<sup>a</sup></td><td>0.020<sup>bc</sup></td><td>0.023<sup>c</sup></td><td>&lt; 0.001</td><td>0.11</td></tr><tr><th>R<sup>λ<sub>3</sub></sup></th><td>4.40</td><td>3.98</td><td>4.03</td><td>4.08</td><td>4.29</td><td>0.36</td><td>0.37</td></tr></tbody></table>

**Table 2.** Effect of the shift from the standard to a high-concentrate diet on synthetic variables describing the post-prandial kinetics of rumen hydrogen ion concentration [H⁺] along the 5 weeks of experimentation. Means without a common superscript letter differ significantly (*p* < 0.05, a-d within a row). A<sup>λ<sub>1</sub></sup>, v<sub>0</sub><sup>λ<sub>2</sub></sup> and R<sup>λ<sub>3</sub></sup> are the power transformed values of A and v<sub>0</sub> and R. W1 is associated with the time interval of distribution of the standard diet, and W2 to W5 with the distribution of the high-concentrate diet. Week and days describe the significance of fixed effects of the mixed model.

<table><thead><tr><th>Goat #</th><th>d<sub>1→2</sub></th><th>d<sub>2→3</sub></th><th>d<sub>3→4</sub></th><th>d<sub>4→5</sub></th><th>d<sub>5→6</sub></th><th>d<sub>6→7</sub></th><th>d<sub>7→8</sub></th><th>d<sub>8→9</sub></th><th>Index</th></tr></thead><tbody><tr><td>1</td><td>0.89</td><td>-5.64</td><td>-17.74</td><td>-2.22</td><td>-0.64</td><td>1.1</td><td>2.9</td><td>-4.19</td><td>-26.42</td></tr><tr><td>2</td><td>-1.09</td><td>-2.14</td><td>-6.36</td><td>-1.34</td><td>6.2</td><td>-1.15</td><td>0.82</td><td>-2.69</td><td>-6.66</td></tr><tr><td>3</td><td>-0.8</td><td>-1.24</td><td>-13.91</td><td>13.12</td><td>-4.39</td><td>1.41</td><td>-2.38</td><td>0.88</td><td>-6.51</td></tr><tr><td>4</td><td>-0.72</td><td>-9.5</td><td>-1.3</td><td>10.84</td><td>4.06</td><td>-3.87</td><td>-1.53</td><td>-2.71</td><td>-4.01</td></tr><tr><td>5</td><td>-0.27</td><td>-2.38</td><td>-6.83</td><td>5.64</td><td>-3.77</td><td>2.98</td><td>-1.53</td><td>1.19</td><td>-4.7</td></tr><tr><td>6</td><td>-1.25</td><td>5.99</td><td>-5.78</td><td>0.6</td><td>-3.97</td><td>1.54</td><td>-0.99</td><td>-1.7</td><td>-4.31</td></tr><tr><td>7</td><td>3.08</td><td>7.8</td><td>-1.11</td><td>-16.19</td><td>-14.88</td><td>4.89</td><td>-0.96</td><td>-4.08</td><td>-24.53</td></tr><tr><td>8</td><td>0.7</td><td>-4.6</td><td>2.57</td><td>-3.76</td><td>5.77</td><td>-0.54</td><td>0.46</td><td>-3.67</td><td>-3.76</td></tr></tbody></table>

**Table 3.** Daily scores and global index for goats' response to high-concentrate diet challenge. *d*<sub>*i*→*i*+1</sub> is the value assigned to the quality of the response of each goat from day *d*<sub>*i*</sub> to *d*<sub>*i*+1</sub> for *i* ∈ {1, .., 8}. Decreasing values of A from day *d*<sub>*i*</sub> to *d*<sub>*i*+1</sub> are associated with the capacity to adapt to the new diet, which is recognised with a positive daily score. Conversely, increasing values of A from day *d*<sub>*i*</sub> to *d*<sub>*i*+1</sub> show that goats had difficulties adapting to the new diet, and negative daily scores are assigned. The value of the score depends on the variation of *v*<sub>0</sub> from day *d*<sub>*i*</sub> to *d*<sub>*i*+1</sub> and also the weighting *ω*<sub>*j*</sub> associated with the relative Euclidian distance between *d*<sub>*i*</sub> and *d*<sub>*i*+1</sub>.

and then the mixed model (Model 2) was used on power transformed variables v<sub>0</sub><sup>λ<sub>1</sub></sup>, A<sup>λ<sub>2</sub></sup> and R<sup>λ<sub>3</sub></sup>. When significant differences were detected, differences among means were tested using Tukey's comparison test with FDR correction (false discovery rate, emmeans package of R). For all variables, the Gaussian distribution of residuals and the heteroscedasticity of the model were checked. Results are presented in Table 2.

Variable v<sub>0</sub><sup>λ<sub>1</sub></sup> for λ<sub>1</sub> = 0.20 was not affected by the fixed effect of “days” (χ<sup>2</sup>(1) = 38.7, *p* = 0.11), but was significantly affected by the fixed effect of “week” (χ<sup>2</sup>(1) = 2.51, *p* < 0.001). The synthetic variable A<sup>λ<sub>2</sub></sup> for λ<sub>2</sub> = 0.28 was not affected by the effect of “days” (χ<sup>2</sup>(1) = 2.52, *p* = 0.11), but was significantly affected by the fixed effect of “week” (χ<sup>2</sup>(1) = 28.4, *p* < 0.001). The anova test showed that the effects of “days” and “week” were not significant (*p* > 0.05) on post-prandial recovery capacity R.

The contribution of the individual variability to the total variance of each of the synthetic variables can be presented as the coefficient ρ = (σ<sup>2</sup><sub>A</sub> - 28) / (σ<sup>2</sup><sub>A</sub> + σ<sup>2</sup><sub>R</sub>). σ<sup>2</sup><sub>A</sub> is the variance associated with the goat effect (between-goats variation) and σ<sup>2</sup> is the variance of residuals (within-goats variation). For the synthetic variable A, 17% of the variation was due to between-goats variability (*ρ* = 0.17), more than 50% of the total variance of v<sub>0</sub> was described by the goat effect (*ρ* = 0.56), and only 10% of total variance of R was due to individual variability.

**Step 3: Quantification of goats' response.** Only synthetic variables with a significant long-term effect (v<sub>0</sub> and A) were used to calculate goats' capacity to adapt to the diet change. Table 3 gives daily scores and global index for each goat. For a better presentation of the results, all scores were rescaled by 10<sup>7</sup>: the higher the global index, the better the adaptive capacity of the animal. The global index for all the goats was negative. The results in Table 3 show six goats with larger global indices (> -10), while two others had rather small indices (< -20). Goat 8 had the best adaptive response (least negative), and Goat 1 the poorest response in terms of capacity to adapt to the new diet.

Figure 7 shows two main strategies of response in goats. The first strategy involved keeping variations of [H⁺] as small as possible, resulting in smaller Euclidian distances (Goats 8, 2 and 6). The second strategy consisted in an oscillating response, mainly in variations of A (Goats 3 and 5).

**Analysis of pH data.** The effect of the diet change was studied on the synthetic variables calculated from pH, for the purpose of comparing with [H⁺] results.
<table><caption>**Table 2.** Effect of the shift from the standard to a high-concentrate diet on synthetic variables describing the post-prandial kinetics of rumen hydrogen ion concentration [H⁺] along the 5 weeks of experimentation. Means without a common superscript letter differ significantly (*p* < 0.05, a-d within a row). A<sup>λ<sub>1</sub></sup>, v<sub>0</sub><sup>λ<sub>2</sub></sup> and R<sup>λ<sub>3</sub></sup> are the power transformed values of A and v<sub>0</sub> and R. W1 is associated with the time interval of distribution of the standard diet, and W2 to W5 with the distribution of the high-concentrate diet. Week and days describe the significance of fixed effects of the mixed model.</caption><thead><tr><th rowspan="2"></th><th colspan="5">Weeks</th><th colspan="2">p-values</th></tr><tr><th>W1</th><th>W2</th><th>W3</th><th>W4</th><th>W5</th><th>Week</th><th>Days</th></tr></thead><tbody><tr><th>v<sub>0</sub><sup>λ<sub>1</sub></sup></th><td>0.052<sup>a</sup></td><td>0.046<sup>c</sup></td><td>0.049<sup>bc</sup></td><td>0.047<sup>c</sup></td><td>0.054<sup>b</sup></td><td>&lt; 0.001</td><td>0.11</td></tr><tr><th>A<sup>λ<sub>2</sub></sup></th><td>0.016<sup>a</sup></td><td>0.017<sup>ab</sup></td><td>0.018<sup>a</sup></td><td>0.020<sup>bc</sup></td><td>0.023<sup>c</sup></td><td>&lt; 0.001</td><td>0.11</td></tr><tr><th>R<sup>λ<sub>3</sub></sup></th><td>4.40</td><td>3.98</td><td>4.03</td><td>4.08</td><td>4.29</td><td>0.36</td><td>0.37</td></tr></tbody></table>

<table><caption>**Table 3.** Daily scores and global index for goats' response to high-concentrate diet challenge. *d*<sub>*i*→*i*+1</sub> is the value assigned to the quality of the response of each goat from day *d*<sub>*i*</sub> to *d*<sub>*i*+1</sub> for *i* ∈ {1, .., 8}. Decreasing values of A from day *d*<sub>*i*</sub> to *d*<sub>*i*+1</sub> are associated with the capacity to adapt to the new diet, which is recognised with a positive daily score. Conversely, increasing values of A from day *d*<sub>*i*</sub> to *d*<sub>*i*+1</sub> show that goats had difficulties adapting to the new diet, and negative daily scores are assigned. The value of the score depends on the variation of *v*<sub>0</sub> from day *d*<sub>*i*</sub> to *d*<sub>*i*+1</sub> and also the weighting *ω*<sub>*j*</sub> associated with the relative Euclidian distance between *d*<sub>*i*</sub> and *d*<sub>*i*+1</sub>.</caption><thead><tr><th>Goat #</th><th>d<sub>1→2</sub></th><th>d<sub>2→3</sub></th><th>d<sub>3→4</sub></th><th>d<sub>4→5</sub></th><th>d<sub>5→6</sub></th><th>d<sub>6→7</sub></th><th>d<sub>7→8</sub></th><th>d<sub>8→9</sub></th><th>Index</th></tr></thead><tbody><tr><td>1</td><td>0.89</td><td>-5.64</td><td>-17.74</td><td>-2.22</td><td>-0.64</td><td>1.1</td><td>2.9</td><td>-4.19</td><td>-26.42</td></tr><tr><td>2</td><td>-1.09</td><td>-2.14</td><td>-6.36</td><td>-1.34</td><td>6.2</td><td>-1.15</td><td>0.82</td><td>-2.69</td><td>-6.66</td></tr><tr><td>3</td><td>-0.8</td><td>-1.24</td><td>-13.91</td><td>13.12</td><td>-4.39</td><td>1.41</td><td>-2.38</td><td>0.88</td><td>-6.51</td></tr><tr><td>4</td><td>-0.72</td><td>-9.5</td><td>-1.3</td><td>10.84</td><td>4.06</td><td>-3.87</td><td>-1.53</td><td>-2.71</td><td>-4.01</td></tr><tr><td>5</td><td>-0.27</td><td>-2.38</td><td>-6.83</td><td>5.64</td><td>-3.77</td><td>2.98</td><td>-1.53</td><td>1.19</td><td>-4.7</td></tr><tr><td>6</td><td>-1.25</td><td>5.99</td><td>-5.78</td><td>0.6</td><td>-3.97</td><td>1.54</td><td>-0.99</td><td>-1.7</td><td>-4.31</td></tr><tr><td>7</td><td>3.08</td><td>7.8</td><td>-1.11</td><td>-16.19</td><td>-14.88</td><td>4.89</td><td>-0.96</td><td>-4.08</td><td>-24.53</td></tr><tr><td>8</td><td>0.7</td><td>-4.6</td><td>2.57</td><td>-3.76</td><td>5.77</td><td>-0.54</td><td>0.46</td><td>-3.67</td><td>-3.76</td></tr></tbody></table>

and then the mixed model (Model 2) was used on power transformed variables v<sub>0</sub><sup>λ<sub>1</sub></sup>, A<sup>λ<sub>2</sub></sup> and R<sup>λ<sub>3</sub></sup>. When significant differences were detected, differences among means were tested using Tukey's comparison test with FDR correction (false discovery rate, emmeans package of R). For all variables, the Gaussian distribution of residuals and the heteroscedasticity of the model were checked. Results are presented in Table 2.

Variable v<sub>0</sub><sup>λ<sub>1</sub></sup> for λ<sub>1</sub> = 0.20 was not affected by the fixed effect of “days” (χ<sup>2</sup>(1) = 38.7, *p* = 0.11), but was significantly affected by the fixed effect of “week” (χ<sup>2</sup>(1) = 2.51, *p* < 0.001). The synthetic variable A<sup>λ<sub>2</sub></sup> for λ<sub>2</sub> = 0.28 was not affected by the effect of “days” (χ<sup>2</sup>(1) = 2.52, *p* = 0.11), but was significantly affected by the fixed effect of “week” (χ<sup>2</sup>(1) = 28.4, *p* < 0.001). The anova test showed that the effects of “days” and “week” were not significant (*p* > 0.05) on post-prandial recovery capacity R.

The contribution of the individual variability to the total variance of each of the synthetic variables can be presented as the coefficient $\rho = \frac{\sigma_A^2 - 28}{\sigma_A^2 + \sigma_A^2}$. σ<sup>2</sup><sub>A</sub> is the variance associated with the goat effect (between-goats variation) and σ<sup>2</sup> is the variance of residuals (within-goats variation). For the synthetic variable A, 17% of the variation was due to between-goats variability (*ρ* = 0.17), more than 50% of the total variance of v<sub>0</sub> was described by the goat effect (*ρ* = 0.56), and only 10% of total variance of R was due to individual variability.

**Step 3: Quantification of goats' response.** Only synthetic variables with a significant long-term effect (v<sub>0</sub> and A) were used to calculate goats' capacity to adapt to the diet change. Table 3 gives daily scores and global index for each goat. For a better presentation of the results, all scores were rescaled by 10<sup>7</sup>: the higher the global index, the better the adaptive capacity of the animal. The global index for all the goats was negative. The results in Table 3 show six goats with larger global indices (> -10), while two others had rather small indices (< -20). Goat 8 had the best adaptive response (least negative), and Goat 1 the poorest response in terms of capacity to adapt to the new diet.

Figure 7 shows two main strategies of response in goats. The first strategy involved keeping variations of [H⁺] as small as possible, resulting in smaller Euclidian distances (Goats 8, 2 and 6). The second strategy consisted in an oscillating response, mainly in variations of A (Goats 3 and 5).

**Analysis of pH data.** The effect of the diet change was studied on the synthetic variables calculated from pH, for the purpose of comparing with [H⁺] results.
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Phase trajectories of goats in response to the high-concentrate diet challenge. Each graph shows the trajectories of each individual goat in terms of A and $v_0$, represented by dashed lines. Goat numbers from 1 to 8 are stated above each graph. Sampling days (from 1 to 9) are in blue inside each graph. Both variables $v_0$ and A are rescaled by $10^7$.</figcaption></figure>

The initial pH (variable $v_0$) was significantly affected by the fixed effects of both “week” ($\chi^2(1) = 72.70, p < 0.001$) and “days” ($\chi^2(1) = 11.7, p < 0.001$). In all, 81% of the variation of $v_0$ was related to the individual variability. The deviation from the initial pH (variable A) was significantly affected by the fixed effect of “week” ($\chi^2(1) = 17.1, p < 0.01$), while the effect of “days” was not significant. Results showed that 38% of the variation of A was related to the individual variability of goats. No effect of diet on the recovery capacity was observed, and only 9% of the total variation was related to individual variability (see supplementary Table S2 online).

Discussion

An original modelling procedure was developed to quantify goats' capacity to adapt to the challenge of a high-concentrate diet, at the individual level. Using the metric defined in Model 3, a global index was assigned to each goat, which allowed ranking the goats for their capacity to buffer the effect of the new diet. A daily score was also assigned, which made it possible to compare the evolution of goats' strategy of response. Results in Table 3 show that Goat 8 had the best adaptive capacity. Although several negative daily scores for this goat could be observed, small variations in $[H^+]$ (small Euclidian distance) resulted in a high global index. Despite the close scores of Goats 1 and 7 (the poorest responses), a large difference between their daily scores was found. This indicates that these goats exhibited different underlying mechanisms and daily strategies to adapt to the effect of the high-concentrate diet. The use of daily scores in Table 3, together with other systemic measures (e.g. plasma metabolites, feed intake, rumination, etc.) will help shed light on the underlying mechanisms of animal adaptive strategies.

Differences in daily strategies of response (daily scores) justifies the large between-goats variability of $v_0$ and A, that has been presented in Result section. The between-goats variability of $v_0$ is linked to the animal's rumen microbial community, which ferments the diet, to the buffering effect of saliva, and to the absorptive capacity of the animal's rumen wall³⁵. A large between-goats variability for variable A, can be explained by the large variation in rumen acid status for some goats, while it was more stable for others. The rumen acid status has an impact on the rumen microbial community, as a marked drop in pH (a large increase in A), which decreases its richness and its diversity, as observed by Zhang et al.³⁶ in goats fed high-grain diets. It has been shown that for ruminants fed the same diet, some develop acidosis while others keep their neutral pH as observed in cows³⁷ or goats³⁸. It seems that feeding behaviour³⁹ and rumen adaptive capacity⁴⁰ are among the main factors regulating animals' response to a change of diet. This link needs to be studied in greater depth, and in particular the link between feeding behaviour and evolution of microbiota, and the link between feeding behaviour and feed efficiency⁴¹.

The global index was developed for two indicators of the animals' adaptive capacity. However, this method can be transposed to more than two indicators. In the case of *n* indicators of performance, an *n*-dimensional
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Phase trajectories of goats in response to the high-concentrate diet challenge. Each graph shows the trajectories of each individual goat in terms of A and $v_0$, represented by dashed lines. Goat numbers from 1 to 8 are stated above each graph. Sampling days (from 1 to 9) are in blue inside each graph. Both variables $v_0$ and A are rescaled by $10^7$.</figcaption></figure>

The initial pH (variable $v_0$) was significantly affected by the fixed effects of both “week” ($\chi^2(1) = 72.70, p < 0.001$) and “days” ($\chi^2(1) = 11.7, p < 0.001$). In all, 81% of the variation of $v_0$ was related to the individual variability. The deviation from the initial pH (variable A) was significantly affected by the fixed effect of “week” ($\chi^2(1) = 17.1, p < 0.01$), while the effect of “days” was not significant. Results showed that 38% of the variation of A was related to the individual variability of goats. No effect of diet on the recovery capacity was observed, and only 9% of the total variation was related to individual variability (see supplementary Table S2 online).

Discussion

An original modelling procedure was developed to quantify goats' capacity to adapt to the challenge of a high-concentrate diet, at the individual level. Using the metric defined in Model 3, a global index was assigned to each goat, which allowed ranking the goats for their capacity to buffer the effect of the new diet. A daily score was also assigned, which made it possible to compare the evolution of goats' strategy of response. Results in Table 3 show that Goat 8 had the best adaptive capacity. Although several negative daily scores for this goat could be observed, small variations in $[H^+]$ (small Euclidian distance) resulted in a high global index. Despite the close scores of Goats 1 and 7 (the poorest responses), a large difference between their daily scores was found. This indicates that these goats exhibited different underlying mechanisms and daily strategies to adapt to the effect of the high-concentrate diet. The use of daily scores in Table 3, together with other systemic measures (e.g. plasma metabolites, feed intake, rumination, etc.) will help shed light on the underlying mechanisms of animal adaptive strategies.

Differences in daily strategies of response (daily scores) justifies the large between-goats variability of $v_0$ and A, that has been presented in Result section. The between-goats variability of $v_0$ is linked to the animal's rumen microbial community, which ferments the diet, to the buffering effect of saliva, and to the absorptive capacity of the animal's rumen wall³⁵. A large between-goats variability for variable A, can be explained by the large variation in rumen acid status for some goats, while it was more stable for others. The rumen acid status has an impact on the rumen microbial community, as a marked drop in pH (a large increase in A), which decreases its richness and its diversity, as observed by Zhang et al.³⁶ in goats fed high-grain diets. It has been shown that for ruminants fed the same diet, some develop acidosis while others keep their neutral pH as observed in cows³⁷ or goats³⁸. It seems that feeding behaviour³⁹ and rumen adaptive capacity⁴⁰ are among the main factors regulating animals' response to a change of diet. This link needs to be studied in greater depth, and in particular the link between feeding behaviour and evolution of microbiota, and the link between feeding behaviour and feed efficiency⁴¹.

The global index was developed for two indicators of the animals' adaptive capacity. However, this method can be transposed to more than two indicators. In the case of *n* indicators of performance, an *n*-dimensional
phase space is needed instead of a phase plane (in two dimensions). The definition of the metric in Model 3 should thus be slightly modified. In addition, the metric used in this work was weighted by the number of days between two sampling days $d_i$ and $d_{i+1}$, which was a strong assumption made to cancel the effect of variable intervals between different recording days (Fig. 1). The optimal use of this approach would be when the indicator of performance is recorded with the same time intervals. Comparing results of weighted and unweighted metrics (see supplementary Table S3 online), the least adaptive animals kept their rank, but the other animals with better capacity to adapt changed rank. This suggests to be cautious in the use of weighted or unweighted metrics and for adapting them to the traits under study.

The method proposed here draws on the method of health trajectories first introduced by Schneider²⁹, and later applied by other authors to distinguish between resistance and tolerance capacities of animals in response to pathogens¹⁷,⁴². This is of major interest in refining selection strategies for an improved adaptive capacity in the face of pathogens. In this approach, the phase plane of pathogen load versus an indicator of performance such as body weight, feed intake or other health indicators was plotted. A score was then assigned to each time step depending on the health status of the animal. Using this method, the authors concluded on the dynamic interrelationship between resistance and tolerance in the animal's response. In our study, the phase plane was presented for two indicators of a goat's performance, i.e., the phase plane represented a response-response plot, which differed slightly from the strategy of presenting a cause-response plot (e.g. pathogen load versus performance).

Other models have been developed to quantify animals' response to perturbations, using time series data¹³⁻¹⁶. In all these models, data were available at a high frequency, and on a given time scale, whereas in the present work, the model had to consider two embedded time scales. One originality of this work was the development of a modelling procedure in three steps that enabled these two time scales to be combined to produce one global index. A mathematical model was developed to transform the post-prandial kinetics of [H⁺] into three synthetic variables. Other authors have also suggested using synthetic variables to describe and summarise the post-prandial kinetics of pH. Molina-Alcaide et al.⁴³ used the average pH to study the influence of diet on ruminal fermentation. Despite the utility of average pH, it does not consider the dynamic characteristics of post-prandial pH. Other synthetic variables were proposed by Dragomir et al.⁷ to analyse the post-prandial pH evolution during the 8 h after delivery of the diet. In their study, two main types of synthetic variables were proposed, the first calculated directly from the initial data and the second corresponding to the parameters of a cubic model ($pH_t = at^3 + bt^2 + ct + d$), or derived therefrom. They finally used a PCA analysis to study the relevant descriptors of the post-prandial evolution of pH. In our approach, given the five post-prandial measures of pH, a quadratic function with three parameters was used to avoid over-parametrisation of the model⁴⁴. The re-parametrisation of the model allowed coefficients *a*, *b* and *c* to be replaced by a set of variables with biological meaningful parameters. The synthetic variable $v_0$ corresponds to the rumen status of the animal just before the morning feed allowance. Since previous feed delivery was performed 16 h before the morning sampling, $v_0$ was associated with the basal level, where there is a balance between the production of volatile fatty acids from fermented feeds in the rumen, elimination by absorption through the rumen wall, and travel to the lower digestive tract⁴⁵. The synthetic variable *A* represents the maximum within-day variation of the acid status of the animal. Large variation of *A* is therefore associated with large variations in the rumen microflora in its environment, which makes it less efficient, entailing lower feed efficiency⁵. A large variation of *A* might damage the rumen wall, with serious pathological consequences⁴⁶.

The results confirmed the hypothesis that variables *A* and $v_0$ were not affected by the short effect of dietary change (fix effect of “days”), illustrating the capacity of goats to overcome the high-concentrate diet challenge in the short term. These variables were little affected by the fixed effect of “week” before Week 4. This suggests that during the first weeks after the diet change, goats succeeded in regulating [H⁺]. These results can be explained by the ability of goats to modify their feeding behaviour when fed a high-concentrate diet. For example, they could decrease their feed intake, their intake rate and increase the time spent ruminating per unit of intake³⁸,⁴⁷. In some cases, they can also look for fibres, which are less fermentable than starch³⁹. Other authors reported some changes in behaviour of cows in cases of disease⁴⁸. At Week 5, a significant increase in *A* and $v_0$ was observed. However, the authors suggest that the response of only 1 week is not enough to conclude on adaptive capacity in the long term. For a better understanding of the mechanisms underlying these results, they must be analysed together with feed intake and time spent ruminating. This analysis lay outside the scope of this work, which focused more narrowly on the modelling aspects of merging the time scales and developing the approach to quantifying the global index.

The reported threshold of acidosis $\theta$ ranges from 5.5 to 6.0⁴⁹,⁵⁰. This should be considered in the definition of synthetic variables. In the model developed here, all three synthetic variables for [H⁺] are threshold-independent. This means that the analysis is valid whatever the value of the threshold for acidosis. This is of major interest when comparing small and large ruminants. The literature shows that goats have a higher pH (+0.4 pH) than cattle when fed diets with similar NDF (neutral detergent fibre) or concentrate percentage³⁴, which implies that the threshold under which animals develop acidosis might be different for large and small ruminants.

Importantly, the greater the number of animals being studied, the more meaningful is the between-goats variability and the statistical significance of the results, whence the utility of pooling several trials and performing a meta-analysis on different data sets⁵¹. Given the small number of goats studied, the statistical power of the tests was estimated. Assuming that the true parameter values were those estimated from the sample, the statistical power for the test of the day effect was around 40% with eight goats. With 20 goats, it reached 70%, while for the week effect, the power was almost 100% even with eight goats. This means a sufficient power for the week effect, but a rather low power for the days effect. This was expected, because of the smaller number of data to estimate the days effect (four observations per goat) than to estimate the week effect.
phase space is needed instead of a phase plane (in two dimensions). The definition of the metric in Model 3 should thus be slightly modified. In addition, the metric used in this work was weighted by the number of days between two sampling days $d_i$ and $d_{i+1}$, which was a strong assumption made to cancel the effect of variable intervals between different recording days (Fig. 1). The optimal use of this approach would be when the indicator of performance is recorded with the same time intervals. Comparing results of weighted and unweighted metrics (see supplementary Table S3 online), the least adaptive animals kept their rank, but the other animals with better capacity to adapt changed rank. This suggests to be cautious in the use of weighted or unweighted metrics and for adapting them to the traits under study.

The method proposed here draws on the method of health trajectories first introduced by Schneider²⁹, and later applied by other authors to distinguish between resistance and tolerance capacities of animals in response to pathogens¹⁷,⁴². This is of major interest in refining selection strategies for an improved adaptive capacity in the face of pathogens. In this approach, the phase plane of pathogen load versus an indicator of performance such as body weight, feed intake or other health indicators was plotted. A score was then assigned to each time step depending on the health status of the animal. Using this method, the authors concluded on the dynamic interrelationship between resistance and tolerance in the animal's response. In our study, the phase plane was presented for two indicators of a goat's performance, i.e., the phase plane represented a response-response plot, which differed slightly from the strategy of presenting a cause-response plot (e.g. pathogen load versus performance).

Other models have been developed to quantify animals' response to perturbations, using time series data¹³⁻¹⁶. In all these models, data were available at a high frequency, and on a given time scale, whereas in the present work, the model had to consider two embedded time scales. One originality of this work was the development of a modelling procedure in three steps that enabled these two time scales to be combined to produce one global index. A mathematical model was developed to transform the post-prandial kinetics of [H⁺] into three synthetic variables. Other authors have also suggested using synthetic variables to describe and summarise the post-prandial kinetics of pH. Molina-Alcaide et al.⁴³ used the average pH to study the influence of diet on ruminal fermentation. Despite the utility of average pH, it does not consider the dynamic characteristics of post-prandial pH. Other synthetic variables were proposed by Dragomir et al.⁷ to analyse the post-prandial pH evolution during the 8 h after delivery of the diet. In their study, two main types of synthetic variables were proposed, the first calculated directly from the initial data and the second corresponding to the parameters of a cubic model ($pH_t = at^3 + bt^2 + ct + d$), or derived therefrom. They finally used a PCA analysis to study the relevant descriptors of the post-prandial evolution of pH. In our approach, given the five post-prandial measures of pH, a quadratic function with three parameters was used to avoid over-parametrisation of the model⁴⁴. The re-parametrisation of the model allowed coefficients *a*, *b* and *c* to be replaced by a set of variables with biological meaningful parameters. The synthetic variable $v_0$ corresponds to the rumen status of the animal just before the morning feed allowance. Since previous feed delivery was performed 16 h before the morning sampling, $v_0$ was associated with the basal level, where there is a balance between the production of volatile fatty acids from fermented feeds in the rumen, elimination by absorption through the rumen wall, and travel to the lower digestive tract⁴⁵. The synthetic variable *A* represents the maximum within-day variation of the acid status of the animal. Large variation of *A* is therefore associated with large variations in the rumen microflora in its environment, which makes it less efficient, entailing lower feed efficiency⁵. A large variation of *A* might damage the rumen wall, with serious pathological consequences⁴⁶.

The results confirmed the hypothesis that variables *A* and $v_0$ were not affected by the short effect of dietary change (fix effect of “days”), illustrating the capacity of goats to overcome the high-concentrate diet challenge in the short term. These variables were little affected by the fixed effect of “week” before Week 4. This suggests that during the first weeks after the diet change, goats succeeded in regulating [H⁺]. These results can be explained by the ability of goats to modify their feeding behaviour when fed a high-concentrate diet. For example, they could decrease their feed intake, their intake rate and increase the time spent ruminating per unit of intake³⁸,⁴⁷. In some cases, they can also look for fibres, which are less fermentable than starch³⁹. Other authors reported some changes in behaviour of cows in cases of disease⁴⁸. At Week 5, a significant increase in *A* and $v_0$ was observed. However, the authors suggest that the response of only 1 week is not enough to conclude on adaptive capacity in the long term. For a better understanding of the mechanisms underlying these results, they must be analysed together with feed intake and time spent ruminating. This analysis lay outside the scope of this work, which focused more narrowly on the modelling aspects of merging the time scales and developing the approach to quantifying the global index.

The reported threshold of acidosis $\theta$ ranges from 5.5 to 6.0⁴⁹,⁵⁰. This should be considered in the definition of synthetic variables. In the model developed here, all three synthetic variables for [H⁺] are threshold-independent. This means that the analysis is valid whatever the value of the threshold for acidosis. This is of major interest when comparing small and large ruminants. The literature shows that goats have a higher pH (+0.4 pH) than cattle when fed diets with similar NDF (neutral detergent fibre) or concentrate percentage³⁴, which implies that the threshold under which animals develop acidosis might be different for large and small ruminants.

Importantly, the greater the number of animals being studied, the more meaningful is the between-goats variability and the statistical significance of the results, whence the utility of pooling several trials and performing a meta-analysis on different data sets⁵¹. Given the small number of goats studied, the statistical power of the tests was estimated. Assuming that the true parameter values were those estimated from the sample, the statistical power for the test of the day effect was around 40% with eight goats. With 20 goats, it reached 70%, while for the week effect, the power was almost 100% even with eight goats. This means a sufficient power for the week effect, but a rather low power for the days effect. This was expected, because of the smaller number of data to estimate the days effect (four observations per goat) than to estimate the week effect.
The differences between results for [H⁺] and pH emphasise that the systematic use of pH as an indicator of rumen acidosis can lead to errors, in particular when using statistical tests or mathematical calculations, as explained also by Murphy¹⁹.

In conclusion, this work proposes a novel modelling approach to quantify the adaptive response of goats to a high-concentrate diet in terms of [H⁺] variations, sampling at two different time scales. This model assigns daily scores and a single global index to each goat's response. The method is generic and can be adapted to several indicators of performance, and to other types of perturbations with several sampling time scales.

## Data availability

The R script of the model and data are available in public repository Zenodo at https://doi.org/10.5281/zenodo.0.3741774.

Received: 10 April 2020; Accepted: 2 November 2020

Published online: 23 November 2020

## References

1. Hofmann, R. R. Evolutionary steps of ecophysiological adaptation and diversification of ruminants—A comparative view of their digestive-system. *Oecologia* **78**, 443–457 (1989).

2. Krause, K. M. & Oetzel, G. R. Understanding and preventing subacute ruminal acidosis in dairy herds: A review. *Anim. Feed Sci. Technol.* **126**, 215–236 (2006).

3. Mould, F. L. & Ørskov, E. R. Manipulation of rumen fluid pH and its influence on cellulolysis in sacco, dry matter degradation and the rumen microflora of sheep offered either hay or concentrate. *Anim. Feed Sci. Technol.* **10**, 1–14 (1983).

4. Owens, F. N., Secrist, D. S., Hill, W. J. & Gill, D. R. Acidosis in cattle: A review. *J. Anim. Sci.* **76**, 275–286 (1998).

5. Krause, K. M. & Oetzel, G. R. Inducing subacute ruminal acidosis in lactating dairy cows. *J. Dairy Sci.* **88**, 3633–3639 (2005).

6. Gao, X. & Oba, M. Relationship of severity of subacute ruminal acidosis to rumen fermentation, chewing activities, sorting behavior, and milk production in lactating dairy cows fed a high-grain diet. *J. Dairy Sci.* **97**, 3006–3016 (2014).

7. Dragomir, C., Sauvant, D., Peyraud, J.-L., Giger-Reverdin, S. & Michalet-Doreau, B. Meta-analysis of 0 to 8 h post-prandial evolution of ruminal pH. *Animal* **2**, 1437–1448 (2008).

8. Fievez, V., Colman, E., Castro-Montoya, J. M., Stefanov, I. & Vlaeminck, B. Milk odd- and branched-chain fatty acids as biomarkers of rumen function—An update. *Anim. Feed Sci. Technol.* **172**, 51–65 (2012).

9. Friggens, N. C., Blanc, F., Berry, D. P. & Puillet, L. Review: Deciphering animal robustness. A synthesis to facilitate its use in livestock breeding and management. *Animal* https://doi.org/10.1017/S175173111700088X (2017).

10. Ingrisch, J. & Bahn, M. Towards a comparable quantification of resilience. *Trends Ecol. Evol.* **33**, 251–259 (2018).

11. Kitano, H. Grand challenges in systems physiology. *Front. Physiol.* **1**, 3 (2010).

12. Giger-Reverdin, S., Duvaux-Ponter, C., Sauvant, D. & Friggens, N. C. Repeatability of traits for characterizing feed intake patterns in dairy goats: A basis for phenotyping in the precision farming context. *Animal* **14**, 1083–1092 (2020).

13. Ben Abdelkrim, A., Puillet, L., Gomes, P. & Martin, O. Lactation curve model with explicit representation of perturbations as a phenotyping tool for dairy livestock precision farming. *PCI Animal*, bioRxiv, 661249, ver. 4 (2019).

14. Friggens, N. C., Etienne, M. P. & Schmidely, P. Characterizing individual differences in animal responses to a nutritional challenge: Toward improved robustness measures. *J. Dairy Sci.* **99**, 2704–2718 (2016).

15. Nguyen-Ba, H., van Milgen, J. & Taghipoor, M. A procedure to quantify the feed intake response of growing pigs to perturbations. *Animal* **14**, 253–260 (2020).

16. Revilla, M. *et al*. Towards the quantitative characterisation of piglets' robustness to weaning: A modelling approach. *Animal* https://doi.org/10.1017/S1751731119000843 (2019).

17. Lough, G., Lengeling, A., Bergmann, S., Doeschl-Wilson, A. B. & Kyriazakis, I. Health trajectories reveal the dynamic contributions of host genetic resistance and tolerance to infection outcome. *Proc. R. Soc. B Biol. Sci.* **282**, 20152151 (2015).

18. Sørensen, S. Enzymstudien II: Über die Messung und die Bedeutung der Wasserstoffionenkonzentration bei enzymatischen Prozessen, cited by Jensen, W. B. The symbol for pH. *Biochem. Z.* **21**, 131–200 (1909).

19. Murphy, M. R. Analysing and presenting pH data. *J. Dairy Sci.* **65**, 161–163 (1982).

20. INRA. INRA feeding system for ruminants (Wageningen Academic Publishers, Wageningen, 2018). https://doi.org/10.3920/978-90-8686-872-8.

21. ISO. Animal feeding stuffs - Determination of moisture and other volatile matter content. ISO 64961999 (1999).

22. Van Soest, P. J. & Wine, R. H. Use of detergents in the analysis of fibrous feeds. IV. Determination of plant cell-wall constituents. *J. Assoc. Off. Anal. Chem.* **50**, 50–55 (1967).

23. Giger, S., Thivend, P., Sauvant, D., Dorléans, M. & Journaix, P. Etude de l'influence préalable de différents traitements amylolytiques sur la teneur en résidu NDF d'aliments du bétail. (Effect of different amylolytic pretreatments on NDF content in feedstuffs). *Ann. Zootech.* **36**, 39–48 (1987).

24. Robertson, J. B. & Van Soest, P. J. Chapter 8. The detergent system of analysis and its application to human foods. In *The Analysis of Dietary Fiber in Food* (eds James, W. P. T. & Theander, O.) 123–158 (Marcel Dekker, Inc., New York, 1981).

25. AFNOR. Aliments des animaux. Dosage de l'amidon. Méthode enzymatique. (Animal Feedingstuffs. Determination of starch content by an enzymatic method). (1997).

26. Sweeney, R. A. & Rexroad, P. R. Comparison of Leco-FP-228 nitrogen determinator with AOAC copper catalyst Kjeldahl method for crude protein. *J. Assoc. Off. Anal. Chem.* **70**, 1028–1030 (1987).

27. Ratkowsky, D. A. *Applied statistics handbook of nonlinear regression models*. (1990).

28. Demidenko, E. *Mixed Models: Theory and Applications* (Wiley Series in Probability and Statistics) (Wiley-Interscience, New York, 2004).

29. Schneider, D. S. Tracing personalized health curves during infections. *PLoS Biol* **9**, e1001158 (2011).

30. R Core Team. R: A Language and Environment for Statistical Computing. (2019).

31. Searle, S. R., Speed, F. M. & Milliken, G. A. Population marginal means in the linear model: An alternative to least squares means. *Am. Stat.* **34**, 216–221 (1980).

32. Royston, P. Remark AS R94: A remark on algorithm AS 181: The W-test for normality. *Appl. Stat.* **44**, 547 (1995).

33. Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting Linear Mixed-Effects Models Using lme4. *J. Stat. Softw.* **67**, (2015).

34. Sauvant, D., Giger-Reverdin, S. & Peyraud, J.-L. 15. Digestive welfare and rumen acidosis. In *INRA Feeding System for Ruminants* 213–218 (Wageningen Academic Publishers, Wageningen, 2018). https://doi.org/10.3920/978-90-8686-872-8

35. Aschenbach, J. R., Penner, G. B., Stumpff, F. & Gäbel, G. Role of fermentation acid absorption in the regulation of ruminal pH. *J. Anim. Sci.* **89**, 1092–1107 (2011).
The differences between results for [H⁺] and pH emphasise that the systematic use of pH as an indicator of rumen acidosis can lead to errors, in particular when using statistical tests or mathematical calculations, as explained also by Murphy¹⁹.

In conclusion, this work proposes a novel modelling approach to quantify the adaptive response of goats to a high-concentrate diet in terms of [H⁺] variations, sampling at two different time scales. This model assigns daily scores and a single global index to each goat's response. The method is generic and can be adapted to several indicators of performance, and to other types of perturbations with several sampling time scales.

## Data availability

The R script of the model and data are available in public repository Zenodo at https://doi.org/10.5281/zenodo.0.3741774.

Received: 10 April 2020; Accepted: 2 November 2020

Published online: 23 November 2020

## References

1. Hofmann, R. R. Evolutionary steps of ecophysiological adaptation and diversification of ruminants—A comparative view of their digestive-system. *Oecologia* **78**, 443–457 (1989).

2. Krause, K. M. & Oetzel, G. R. Understanding and preventing subacute ruminal acidosis in dairy herds: A review. *Anim. Feed Sci. Technol.* **126**, 215–236 (2006).

3. Mould, F. L. & Ørskov, E. R. Manipulation of rumen fluid pH and its influence on cellulolysis in sacco, dry matter degradation and the rumen microflora of sheep offered either hay or concentrate. *Anim. Feed Sci. Technol.* **10**, 1–14 (1983).

4. Owens, F. N., Secrist, D. S., Hill, W. J. & Gill, D. R. Acidosis in cattle: A review. *J. Anim. Sci.* **76**, 275–286 (1998).

5. Krause, K. M. & Oetzel, G. R. Inducing subacute ruminal acidosis in lactating dairy cows. *J. Dairy Sci.* **88**, 3633–3639 (2005).

6. Gao, X. & Oba, M. Relationship of severity of subacute ruminal acidosis to rumen fermentation, chewing activities, sorting behavior, and milk production in lactating dairy cows fed a high-grain diet. *J. Dairy Sci.* **97**, 3006–3016 (2014).

7. Dragomir, C., Sauvant, D., Peyraud, J.-L., Giger-Reverdin, S. & Michalet-Doreau, B. Meta-analysis of 0 to 8 h post-prandial evolution of ruminal pH. *Animal* **2**, 1437–1448 (2008).

8. Fievez, V., Colman, E., Castro-Montoya, J. M., Stefanov, I. & Vlaeminck, B. Milk odd- and branched-chain fatty acids as biomarkers of rumen function—An update. *Anim. Feed Sci. Technol.* **172**, 51–65 (2012).

9. Friggens, N. C., Blanc, F., Berry, D. P. & Puillet, L. Review: Deciphering animal robustness. A synthesis to facilitate its use in livestock breeding and management. *Animal* https://doi.org/10.1017/S175173111700088X (2017).

10. Ingrisch, J. & Bahn, M. Towards a comparable quantification of resilience. *Trends Ecol. Evol.* **33**, 251–259 (2018).

11. Kitano, H. Grand challenges in systems physiology. *Front. Physiol.* **1**, 3 (2010).

12. Giger-Reverdin, S., Duvaux-Ponter, C., Sauvant, D. & Friggens, N. C. Repeatability of traits for characterizing feed intake patterns in dairy goats: A basis for phenotyping in the precision farming context. *Animal* **14**, 1083–1092 (2020).

13. Ben Abdelkrim, A., Puillet, L., Gomes, P. & Martin, O. Lactation curve model with explicit representation of perturbations as a phenotyping tool for dairy livestock precision farming. *PCI Animal*, bioRxiv, 661249, ver. 4 (2019).

14. Friggens, N. C., Etienne, M. P. & Schmidely, P. Characterizing individual differences in animal responses to a nutritional challenge: Toward improved robustness measures. *J. Dairy Sci.* **99**, 2704–2718 (2016).

15. Nguyen-Ba, H., van Milgen, J. & Taghipoor, M. A procedure to quantify the feed intake response of growing pigs to perturbations. *Animal* **14**, 253–260 (2020).

16. Revilla, M. *et al*. Towards the quantitative characterisation of piglets' robustness to weaning: A modelling approach. *Animal* https://doi.org/10.1017/S1751731119000843 (2019).

17. Lough, G., Lengeling, A., Bergmann, S., Doeschl-Wilson, A. B. & Kyriazakis, I. Health trajectories reveal the dynamic contributions of host genetic resistance and tolerance to infection outcome. *Proc. R. Soc. B Biol. Sci.* **282**, 20152151 (2015).

18. Sørensen, S. Enzymstudien II: Über die Messung und die Bedeutung der Wasserstoffionenkonzentration bei enzymatischen Prozessen, cited by Jensen, W. B. The symbol for pH. *Biochem. Z.* **21**, 131–200 (1909).

19. Murphy, M. R. Analysing and presenting pH data. *J. Dairy Sci.* **65**, 161–163 (1982).

20. INRA. INRA feeding system for ruminants (Wageningen Academic Publishers, Wageningen, 2018). https://doi.org/10.3920/978-90-8686-872-8.

21. ISO. Animal feeding stuffs - Determination of moisture and other volatile matter content. ISO 64961:1999 (1999).

22. Van Soest, P. J. & Wine, R. H. Use of detergents in the analysis of fibrous feeds. IV. Determination of plant cell-wall constituents. *J. Assoc. Off. Anal. Chem.* **50**, 50–55 (1967).

23. Giger, S., Thivend, P., Sauvant, D., Dorléans, M. & Journaix, P. Etude de l'influence préalable de différents traitements amylolytiques sur la teneur en résidu NDF d'aliments du bétail. (Effect of different amylolytic pretreatments on NDF content in feedstuffs). *Ann. Zootech.* **36**, 39–48 (1987).

24. Robertson, J. B. & Van Soest, P. J. Chapter 8. The detergent system of analysis and its application to human foods. In *The Analysis of Dietary Fiber in Food* (eds James, W. P. T. & Theander, O.) 123–158 (Marcel Dekker, Inc., New York, 1981).

25. AFNOR. Aliments des animaux. Dosage de l'amidon. Méthode enzymatique. (Animal Feedingstuffs. Determination of starch content by an enzymatic method). (1997).

26. Sweeney, R. A. & Rexroad, P. R. Comparison of Leco-FP-228 nitrogen determinator with AOAC copper catalyst Kjeldahl method for crude protein. *J. Assoc. Off. Anal. Chem.* **70**, 1028–1030 (1987).

27. Ratkowsky, D. A. *Applied statistics handbook of nonlinear regression models*. (1990).

28. Demidenko, E. *Mixed Models: Theory and Applications* (Wiley Series in Probability and Statistics) (Wiley-Interscience, New York, 2004).

29. Schneider, D. S. Tracing personalized health curves during infections. *PLoS Biol* **9**, e1001158 (2011).

30. R Core Team. R: A Language and Environment for Statistical Computing. (2019).

31. Searle, S. R., Speed, F. M. & Milliken, G. A. Population marginal means in the linear model: An alternative to least squares means. *Am. Stat.* **34**, 216–221 (1980).

32. Royston, P. Remark AS R94: A remark on algorithm AS 181: The W-test for normality. *Appl. Stat.* **44**, 547 (1995).

33. Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting Linear Mixed-Effects Models Using lme4. *J. Stat. Softw.* **67**, (2015).

34. Sauvant, D., Giger-Reverdin, S. & Peyraud, J.-L. 15. Digestive welfare and rumen acidosis. In *INRA Feeding System for Ruminants* 213–218 (Wageningen Academic Publishers, Wageningen, 2018). https://doi.org/10.3920/978-90-8686-872-8

35. Aschenbach, J. R., Penner, G. B., Stumpff, F. & Gäbel, G. Role of fermentation acid absorption in the regulation of ruminal pH. *J. Anim. Sci.* **89**, 1092–1107 (2011).
36. Zhang, R. Y. et al. Response of rumen microbiota, and metabolic profiles of rumen fluid, liver and serum of goats to high-grain diets. *Animal* **13**, 1855–1864 (2019).

37. Oetzel, G. R. Monitoring and testing dairy herds for metabolic disease. *Vet. Clin. N. Am. Food Anim. Pract.* **20**, 651–674 (2004).

38. Desnoyers, M., Giger-Reverdin, S., Sauvant, D. & Duvaux-Ponter, C. The use of a multivariate analysis to study between-goat variability in feeding behavior and associated rumen pH patterns. *J. Dairy Sci.* **94**, 842–852 (2011).

39. Giger-Reverdin, S. Recent advances in the understanding of subacute ruminal acidosis (SARA) in goats, with focus on the link to feeding behaviour. *Small Rumin. Res.* **163**, 24–28 (2018).

40. Giger-Reverdin, S., Domange, C., Broudiscou, L. P., Sauvant, D. & Berthelot, V. Rumen function in goats, an example of adaptive capacity. *J. Dairy Res.* **87**, 45–51 (2020).

41. Schären, M. et al. Interrelations between the rumen microbiota and production, behavioral, rumen fermentation, metabolic, and immunological attributes of dairy cows. *J. Dairy Sci.* **101**, 4615–4637 (2018).

42. Doeschl-wilson, A. B., Bishop, S. C., Kyriazakis, I., Villanueva, B. & Banos, G. Novel methods for quantifying individual host response to infectious pathogens for genetic analyses. **3**, 1–9 (2012).

43. Molina-Alcaide, E., Pascual, M. R., Cantalapiedra-Hijar, G., Morales-Garcia, E. Y. & Martin-Garcia, A. I. Effects of concentrate replacement by feed blocks on ruminal fermentation and microbial growth in goats and single-flow continuous-culture fermenters. *J. Anim. Sci.* **87**, 1321–1333 (2009).

44. Muñoz-Tamayo, R. et al. Review: To be or not to be an identifiable model. Is this a relevant question in animal science modelling?. *Animal* **12**, 701–712 (2018).

45. Allen, M. S. Relationship between fermentation acid production in the rumen and the requirement for physically effective fiber. *J. Dairy Sci.* **80**, 1447–1462 (1997).

46. Das, S. K., Misra, S. K. & Basak, D. K. Pathological changes in experimental rumen acidosis in goats. *Indian Vet. J.* **69**, 495–497 (1992).

47. Desnoyers, M., Giger-Reverdin, S., Duvaux-Ponter, C. & Sauvant, D. Modeling of off-feed periods caused by subacute acidosis in intensive lactating ruminants: Application to goats. *J. Dairy Sci.* **92**, 3894–3906 (2009).

48. von Keyserlingk, M. A. G. & Weary, D. M. Review: Feeding behaviour of dairy cattle: Meaures and applications. *Can. J. Anim. Sci.* **90**, 303–309 (2010).

49. Calsamiglia, S., Blanch, M., Ferret, A. & Moya, D. Is subacute ruminal acidosis a pH related problem? Causes and tools for its control. *Anim. Feed Sci. Technol.* **172**, 42–50 (2012).

50. Nocek, J. E. Bovine acidosis: Implications on laminitis. *J. Dairy Sci.* **80**, 1005–1028 (1997).

51. Sauvant, D., Schmidely, P., Daudin, J. J. & St-Pierre, N. R. Meta-analyses of experimental data in animal nutrition. *Animal* **2**, 1203–1214 (2008).

## Acknowledgements

The authors thank Joseph Tessier, Alexandra Eymard and all the technical staff of the MoSAR experimental farm at Thiverval-Grignon, for the care of the goats and their work on the experimental protocol.

## Author contributions

S.G.-R. designed the experiment, M.T. performed analyses of data and the modelling procedure. M.D. performed statistical analyses. S.G.-R. and M.T. wrote the manuscript and prepared the figures and tables in the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41598-020-77353-y.

Correspondence and requests for materials should be addressed to M.T.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

**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/.

© The Author(s) 2020
36. Zhang, R. Y. et al. Response of rumen microbiota, and metabolic profiles of rumen fluid, liver and serum of goats to high-grain diets. *Animal* **13**, 1855–1864 (2019).

37. Oetzel, G. R. Monitoring and testing dairy herds for metabolic disease. *Vet. Clin. N. Am. Food Anim. Pract.* **20**, 651–674 (2004).

38. Desnoyers, M., Giger-Reverdin, S., Sauvant, D. & Duvaux-Ponter, C. The use of a multivariate analysis to study between-goat variability in feeding behavior and associated rumen pH patterns. *J. Dairy Sci.* **94**, 842–852 (2011).

39. Giger-Reverdin, S. Recent advances in the understanding of subacute ruminal acidosis (SARA) in goats, with focus on the link to feeding behaviour. *Small Rumin. Res.* **163**, 24–28 (2018).

40. Giger-Reverdin, S., Domange, C., Broudiscou, L. P., Sauvant, D. & Berthelot, V. Rumen function in goats, an example of adaptive capacity. *J. Dairy Res.* **87**, 45–51 (2020).

41. Schären, M. et al. Interrelations between the rumen microbiota and production, behavioral, rumen fermentation, metabolic, and immunological attributes of dairy cows. *J. Dairy Sci.* **101**, 4615–4637 (2018).

42. Doeschl-wilson, A. B., Bishop, S. C., Kyriazakis, I., Villanueva, B. & Banos, G. Novel methods for quantifying individual host response to infectious pathogens for genetic analyses. **3**, 1–9 (2012).

43. Molina-Alcaide, E., Pascual, M. R., Cantalapiedra-Hijar, G., Morales-Garcia, E. Y. & Martin-Garcia, A. I. Effects of concentrate replacement by feed blocks on ruminal fermentation and microbial growth in goats and single-flow continuous-culture fermenters. *J. Anim. Sci.* **87**, 1321–1333 (2009).

44. Muñoz-Tamayo, R. et al. Review: To be or not to be an identifiable model. Is this a relevant question in animal science modelling?. *Animal* **12**, 701–712 (2018).

45. Allen, M. S. Relationship between fermentation acid production in the rumen and the requirement for physically effective fiber. *J. Dairy Sci.* **80**, 1447–1462 (1997).

46. Das, S. K., Misra, S. K. & Basak, D. K. Pathological changes in experimental rumen acidosis in goats. *Indian Vet. J.* **69**, 495–497 (1992).

47. Desnoyers, M., Giger-Reverdin, S., Duvaux-Ponter, C. & Sauvant, D. Modeling of off-feed periods caused by subacute acidosis in intensive lactating ruminants: Application to goats. *J. Dairy Sci.* **92**, 3894–3906 (2009).

48. von Keyserlingk, M. A. G. & Weary, D. M. Review: Feeding behaviour of dairy cattle: Meaures and applications. *Can. J. Anim. Sci.* **90**, 303–309 (2010).

49. Calsamiglia, S., Blanch, M., Ferret, A. & Moya, D. Is subacute ruminal acidosis a pH related problem? Causes and tools for its control. *Anim. Feed Sci. Technol.* **172**, 42–50 (2012).

50. Nocek, J. E. Bovine acidosis: Implications on laminitis. *J. Dairy Sci.* **80**, 1005–1028 (1997).

51. Sauvant, D., Schmidely, P., Daudin, J. J. & St-Pierre, N. R. Meta-analyses of experimental data in animal nutrition. *Animal* **2**, 1203–1214 (2008).

## Acknowledgements

The authors thank Joseph Tessier, Alexandra Eymard and all the technical staff of the MoSAR experimental farm at Thiverval-Grignon, for the care of the goats and their work on the experimental protocol.

## Author contributions

S.G.-R. designed the experiment, M.T. performed analyses of data and the modelling procedure. M.D. performed statistical analyses. S.G.-R. and M.T. wrote the manuscript and prepared the figures and tables in the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41598-020-77353-y.

Correspondence and requests for materials should be addressed to M.T.

Reprints and permissions information is available at www.nature.com/reprints.

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

<figure><img src="image_11.png" /></figure>

**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/.

© The Author(s) 2020