# Using egg production longitudinal recording to study the genetic background of resilience in purebred and crossbred laying hens

Nicolas Bédère, Tom Berghof, Katrijn Peeters, Marie-Hélène Pinard-van Der Laan, Jeroen Visscher, Ingrid David, Han Mulder

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

## ► To cite this version:

Nicolas Bédère, Tom Berghof, Katrijn Peeters, Marie-Hélène Pinard-van Der Laan, Jeroen Visscher, et al.. Using egg production longitudinal recording to study the genetic background of resilience in purebred and crossbred laying hens. Genetics Selection Evolution, 2022, 54 (1), 12 p. 10.1186/s12711-022-00716-8. http://hal-03657027

HAL Id: hal-03657027

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

Submitted on 2 May 2022

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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
<figure><img src="image_3.png" /></figure>

# Using egg production longitudinal recording to study the genetic background of resilience in purebred and crossbred laying hens

Nicolas Bedere¹*, Tom V. L. Berghof²,³, Katrijn Peeters⁴, Marie-Hélène Pinard-van der Laan⁵, Jeroen Visscher⁴, Ingrid David⁶ and Han A. Mulder²

**Abstract**

**Background:** There is growing interest in using genetic selection to obtain more resilient farm animals (i.e. that are minimally affected by disturbances or rapidly recover from them). The aims of this study were to: (i) estimate the genetic parameters of resilience indicator traits based on egg production data, (ii) assess whether these traits are genetically correlated in purebreds and crossbreds, and (iii) assess the genetic correlations of these traits with egg production (EP) as total number of eggs between 25 and 83 weeks. Purebred hens (33,825 from a White Leghorn (WA) line and 34,397 from a Rhode Island (BD) line were housed in individual cages, while crossbred hens were housed in collective cages of 6 to 8 paternal half-sibs (12,852 WA and 3898 BD crossbred groups, where the name of the group refers to the line used as the sire). Deviations of a hen's weekly egg production from the average of the corresponding batch were calculated. Resilience indicator traits investigated were the natural logarithm of the variance (LNVAR), the skewness (SKEW), and the lag-one autocorrelation (AUTO-R) of these deviations.

**Results:** In both purebred lines, EP was estimated to be lowly heritable (WA: 0.11 and BD: 0.12). Resilience indicators were also estimated to be lowly heritable in both lines (LNVAR: 0.10 and 0.12, SKEW: 0.04 and 0.02, AUTO-R: 0.06 and 0.08 in WA and BD, respectively). In both crossbred groups, EP, AUTO-R, and SKEW were estimated to be less heritable than in purebreds (EP: $h^2 \le 0.07$; and resilience indicator traits: $h^2 \le 0.03$), while LNVAR had an $h^2$ estimate that was similar to or higher in crossbreds ($h^2$ ranged from 0.13 to 0.21) than in purebreds. In both purebreds and crossbreds, resilience indicator traits were estimated to have favorable genetic correlations with EP and between each other. For all traits and in both lines, estimates of genetic correlations between purebreds and crossbreds ($r_{pc}$) differed from 1 and ranged from 0.16 to 0.63.

**Conclusions:** These results show that selection for resilience based on EP data can be considered in breeding programs for layers. Genetic improvement of resilience in crossbreds can be achieved by using information on purebreds, but would be greatly enhanced by the integration of information on crossbreds in breeding programs.

**Background**

Over the past 60 years, breeding programs have resulted
in very productive farm animals, at the expense of
functional traits related to fitness (longevity, health,

reproduction, etc.) [1, 2]. These trade-offs between bio-
logical functions can be explained by the “resource allo-
cation theory”: when two biological processes share
the same resources, they are competing in situations of
limited resources [3]. Situations of limited resources are
more likely to occur nowadays because of the increasing
competition for natural resources with other animals,
or with humans [4]. Moreover, situations of additional

<footnote>*Correspondence: nicolas-bedere@inrae.fr
¹ PEGASE, INRAE, Institut Agro, 35590 Saint Gilles, France
Full list of author information is available at the end of the article</footnote>
energy expenses are more frequent because of negative social interactions due to group housing [5] and increasing risks of diseases, temperature fluctuations, for instance with outdoor access and reduced use of preventive antibiotics. Therefore, there is growing interest to select for improved resilience in farm animals, although improvement of general resilience is not commonly implemented in breeding programs [6].

Resilience is defined as “the capacity of an animal to be minimally affected by disturbances or to rapidly return to the state pertained before exposure to a disturbance” [5, 6]. One reason why resilience has not yet been implemented in breeding programs is due to the lack of selection criteria and recording tools. Opportunities to construct selection criteria for resilience lie in the collection of repeated information per individual or per family [7]. Indeed, the calculation of resilience indicator traits from longitudinal data is possible in various species [8]. Animals that are exposed to a disturbance are known to show a different biological response than expected without a disturbance (e.g. growth, feed intake, milk production) [5, 9, 10]. Thus, the most studied resilience indicators are derived from deviations between the observed and expected performance (such as the variance of the deviations) or directly from longitudinal analyses (such as the slope of a reaction norm). These indicators are known to be associated with health-status and periods of perturbations [11–14]. Only a few studies have investigated the genetic parameters of resilience indicator traits derived from longitudinal data in farm animals: dairy cows [15, 16], pigs [17], and layers [18]. Resilience indicator traits based on longitudinal recording of milk production [15, 16], feed intake [17], or body weight [18] have been shown to be heritable and have favorable genetic correlations with health-related traits.

Previous studies in laying hens have investigated resilience indicators that were derived from body weight [18], egg weight [19], or eggshell color [20]. Because these traits were monitored only at specific time points, the number of records was very small in these studies (3 to 8 records during the production cycle per hen). However, breeders record the number of eggs laid daily, at the individual level in nucleus farms or at a group cage level in recurrent test farms. Resilience indicators that are derived from daily data on egg production have not been investigated but may hold potential to implement selection on resilience in breeding programs for layers. The breeding goal of poultry breeding companies is to genetically improve multiple traits in commercial hens, often a 3- or 4-way cross. To do that, they use a pyramidal breeding program, i.e., they perform genetic selection on purebreds, to genetically improve performance of crossbreds. Records on individual performance are required in the

nucleus to select birds with the desired characteristics. Consequently, pure line birds are individually housed to enable individual identification and recording of own performance (e.g. individual egg traits). Two-way crossbred animals are placed in recurrent field tests (RT), in collective cages (i.e. housed in groups) to mimic commercial farm conditions. Crossbred phenotypes are recorded on the group level. These crossbred animals are pedigreed to allow the estimation of genetic parameters and breeding values of purebreds. The purebred–crossbred genetic correlation for a trait is often lower than 1 [21, 22], which implies that crossbred information is valuable to estimate breeding values of purebred birds for crossbred performance [23, 24].

Our objectives were: (i) to study the genetic parameters of resilience indicator traits that are derived from longitudinal egg production data, (ii) to assess whether these traits are genetically correlated in purebreds and crossbreds (purebred–crossbred correlation), and (iii) to assess the genetic correlations of these traits with egg production. To our knowledge, this is the first time that egg production longitudinal recording is used to study resilience in laying hens. The purebred–crossbred correlation determines whether selection for improved resilience in purebreds will result in improved resilience in crossbreds. It also determines the added value of collecting crossbred information: the lower the purebred–crossbred correlation, the higher the value of crossbred information when crossbred performance is the breeding goal [23, 24]. Estimates of genetic correlations with other traits will provide knowledge about the overlap between the genetic background of the traits and practical insights about the integration of resilience traits in breeding programs.

## Methods

### Study populations

In this study, two purebred lines from Hendrix Genetics were studied: a White Leghorn line, called WA, and a Rhode Island line, called BD. The WA purebreds were individually housed on a nucleus farm in The Netherlands (33,825 layers) and distributed among 13 batches from 2012 to 2017. The BD purebreds were individually housed on a nucleus farm in France (34,397 layers) and distributed among 14 batches from 2012 to 2017 (further details about the study design are in Additional file 1: Tables S1 and S2). Recurrent field tests were carried out in Canada (CA) and The Netherlands (NL) by crossing purebred sires from the WA or the BD line with dams from five and seven other lines, respectively. The crossbreds were housed in cages, in groups of 6 to 8 paternal half-sibs. All crossbreds were beak-trimmed but the purebreds were not. The final dataset included records on
12,852 groups of WA crossbreds (6154 records on cages of 6 hens and 6698 records on cages of 7 hens) distributed among four CA farms and six NL farms and 17 batches, and records on 3898 groups of BD crossbreds (2428 records on cages of 7 hens and 1470 records on cages of 8 hens) distributed among one CA farm and four NL farms and nine batches. For crossbred RT, 71 and 29% of the sires used in the WA and BD lines, respectively, had offspring in both countries (CA and NL). The purebred and crossbred populations consisted mostly of half-sib families with, on average, 50 offspring per sire in purebreds and 70 offspring per sire in crossbreds. The pedigree for each recorded individual was traced back between 4 and 16 generations (paternal only for crossbreds).

## Egg production recording

Eggs were collected several times a week (with the interval between two collections ranging from 1 to 4 days). Therefore, the raw data consisted of dates and number of eggs collected. Egg production was recorded individually for purebreds and as a pooled average for a cage for crossbreds, i.e., the total number of eggs produced during a particular period was divided by the average number of hens alive in the cage during that period. Observations for individual purebred hens with two eggs per day were not discarded from the analyses, since the first egg could have been laid after collection time on the day before and the second egg could have been laid before collection time on the day of recording. Observations of three or more eggs per day were discarded from the analyses, because this is not biologically possible, even with different collection times in successive days. Daily egg production was estimated for every day between two recording dates by dividing the total number of eggs by the number of days in the period. These daily averages were summed over 7 days to obtain the weekly egg production (total number of eggs laid per week per hen). The trait egg production (EP) was defined as the total number of eggs laid per hen between 25 and 83 weeks of age.

## Resilience indicator traits based on egg production

The longitudinal egg production data enabled the computation of the types of resilience indicators based on deviations from the average performance of a contemporary group that had been used for other species in other studies [6, 16, 18]. The average weekly number of eggs of the batch was calculated as the total number of eggs divided by the number of hens alive in the batch and used to calculate deviations of individual weekly EP records from the batch average. Resilience indicator traits evaluated were the natural logarithm of the variance (LNVAR) of an individual's or group's deviations, the skewness of the distribution of the deviations

(SKEW), and the lag-one autocorrelation (AUTO-R) of the deviations. These three traits are complementary [6]: a higher LNVAR is an indication of poor resilience of an animal to a perturbation, while resilient animals show a smaller range of deviations from the expected performance and therefore a smaller LNVAR; a more negative SKEW means that the animal shows more negative than positive deviations or that the animal experiences severe drops from the expected performance, which is an indication of poor resilience; and the closer AUTO-R is to 1, the more two subsequent deviations are alike, which is interpreted as a slower recovery from the perturbation and thus an indication of poor resilience. Observations for the resilience traits were discarded if they were outside the mean plus or minus four times the standard deviation within populations (combination of lines and purebreds/crossbreds separately). However, only the observation for the particular resilience indicator trait was discarded, i.e., observations for all other resilience indicators and egg production were kept. This resulted in the removal of about 4 and 8% of the raw data for each resilience indicator trait in the WA and BD lines, respectively.

Data handling and graphs were performed in base R [25]. Resilience indicator traits were calculated in R, i.e. LNVAR with the log and var functions, AUTO-R with the acf function, and SKEW with the skewness function of the e1071 package [26].

## Estimation of variance components for the same trait between purebreds and crossbreds

Within sire lines, variance components were estimated using bivariate analyses, considering records of the purebreds and records of the crossbreds from the same sires' line as two different traits, to allow estimation of the purebred-crossbred correlation. The model used for purebred phenotypes was the following animal model:

$$y_{ijl} = \mu + \text{HatchLoc}_i + \text{SurvClass}_j + a_l + e_{ijl},$$

where y is the observed phenotype (EP, LNVAR, SKEW, AUTO-R) of the lth purebred hen; HatchLoc is the fixed effect for the ith combination of hatch week, barn, corridor, and row of the cage; SurvClass is the fixed effect of the jth class of lifespan duration, with classes comprising intervals of four weeks, from the shortest lifespan observed to the longest (24 to 112 weeks of age), which allowed for the potential effect of survival on phenotypes to be accounted for without assuming a linear relationship between survival and the phenotype (as a covariate would do); a is the random additive genetic effect of the lth purebred layer; and e is the random residual effect.

The model used for crossbred phenotypes was the following sire model:
$$y_{ijkl} = \mu + \text{HatchLoc}_i + \text{SurvClass}_j + \text{MaternalLine}_k + s_l + e_{ijkl},$$

where y is the mean observed phenotype (EP, LNVAR, SKEW, or AUTO-R) from a cage of crossbred halfsib hens from the 1th sire; HatchLoc and SurvClass are the same fixed effects as in the animal model; MaternalLine is the fixed effect of the kth maternal line of the crossbred offspring; s is the random sire genetic effect of the lth purebred sire on phenotype of the crossbred offspring; and e is the random residual effect.

The animal effects for the purebred phenotypes and the sire effects for the crossbred phenotypes were assumed to be bivariate normally distributed as:

$$\begin{bmatrix} \mathbf{a}_p \\ \mathbf{s}_c \end{bmatrix} \sim \text{MVN}\left(\begin{pmatrix} 0 \\ 0 \end{pmatrix}, \mathbf{A} \times \begin{bmatrix} \sigma_{a_p}^2 & \sigma_{a_p s_c} \\ \sigma_{a_p s_c} & \sigma_{s_c}^2 \end{bmatrix}\right) \text{ with } \sigma_{s_c}^2 = \frac{1}{4}\sigma_{a_c}^2,$$

where a_p is the vector of breeding values for purebred performance and s_c is the vector of sire genetic effects for crossbred performance (p indicates purebreds and c indicates crossbreds; A is the additive genetic relationship matrix; σ_a_p^2 is the additive genetic variance for purebred performance (animal model); σ_s_c^2 is the sire genetic variance for crossbred performance (sire model); σ_a_c^2 is the additive genetic variance for crossbred performance; and σ_a_p s_c is the genetic covariance between purebred and crossbred performance (note: σ_a_p s_c = r_p c × σ_a_p × σ_s_c = 1/2 × r_p c × σ_a_p × σ_a_c, where r_p c is the purebred-crossbred genetic correlation). Heterogeneity of residual variance was modeled with a single residual variance for individuals for purebred phenotypes and a separate residual variance for each cage size for crossbred phenotypes (cages of 6 or 7 hens in the WA line and cages of 7 or 8 hens in the BD line). We chose this approach because we observed that the residual variance varied with cage size. The residual effects were assumed to be distributed as:

$$\begin{bmatrix} \mathbf{e}_p \\ \mathbf{e}_{c6} \\ \mathbf{e}_{c7} \end{bmatrix} \sim \text{MVN}\left(\begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix}, \mathbf{I} \times \begin{bmatrix} \sigma_{e_p}^2 & 0 & 0 \\ 0 & \sigma_{e_{c6}}^2 & 0 \\ 0 & 0 & \sigma_{e_{c7}}^2 \end{bmatrix}\right) \text{ for the WA line,}$$

$$\begin{bmatrix} \mathbf{e}_p \\ \mathbf{e}_{c7} \\ \mathbf{e}_{c8} \end{bmatrix} \sim \text{MVN}\left(\begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix}, \mathbf{I} \times \begin{bmatrix} \sigma_{e_p}^2 & 0 & 0 \\ 0 & \sigma_{e_{c7}}^2 & 0 \\ 0 & 0 & \sigma_{e_{c8}}^2 \end{bmatrix}\right) \text{ for the BD line,}$$

where I is the identity matrix, σ_e_p^2 is the residual variance for purebreds, σ_e_c^2 is the residual variance for crossbreds in cages of 6, 7, or 8 grouped hens. A residual covariance between purebred and crossbred performance was not

modeled because hens were present in only one of the three environments.

## Estimation of variance components for different traits within purebreds and crossbreds

Bivariate analyses were used to study the genetic relationship between traits, using the same models presented above (animal model for purebred phenotypes and sire model for crossbred phenotypes). Within each purebred population and for each pair of traits, the genetic effects were assumed to be distributed as:

$$\begin{bmatrix} \mathbf{a}_{\text{trait.1}} \\ \mathbf{a}_{\text{trait.2}} \end{bmatrix} \sim \text{MVN}\left(\begin{pmatrix} 0 \\ 0 \end{pmatrix}, \mathbf{A} \times \begin{bmatrix} \sigma_{a_{\text{trait.1}}}^2 & \sigma_{a_{\text{trait.1,trait.2}}} \\ \sigma_{a_{\text{trait.1,trait.2}}} & \sigma_{a_{\text{trait.2}}}^2 \end{bmatrix}\right),$$

where the subscripts trait.1 and trait.2 indicate the two traits analyzed, A is the additive genetic relationship matrix, σ_a_{ trait.1}^2 is the additive genetic variance, σ_a_{ trait.1,trait.2} is the genetic covariance between the two traits.

Within each crossbred population and for each pair of traits, the genetic effects were assumed to be distributed as:

$$\begin{bmatrix} \mathbf{s}_{\text{trait.1}} \\ \mathbf{s}_{\text{trait.2}} \end{bmatrix} \sim \text{MVN}\left(\begin{pmatrix} 0 \\ 0 \end{pmatrix}, \mathbf{A} \times \begin{bmatrix} \sigma_{s_{\text{trait.1}}}^2 & \sigma_{s_{\text{trait.1,trait.2}}} \\ \sigma_{s_{\text{trait.1,trait.2}}} & \sigma_{s_{\text{trait.2}}}^2 \end{bmatrix}\right) \\ \text{with } \sigma_{s_c}^2 = \frac{1}{4}\sigma_{a_c}^2.$$

The residual effects were assumed to be distributed as:

$$\begin{bmatrix} \mathbf{e}_{\text{trait.1}} \\ \mathbf{e}_{\text{trait.2}} \end{bmatrix} \sim \text{MVN}\left(\begin{pmatrix} 0 \\ 0 \end{pmatrix}, \mathbf{I} \times \begin{bmatrix} \sigma_{e_{\text{trait.1}}}^2 & \sigma_{e_{\text{trait.1,trait.2}}} \\ \sigma_{e_{\text{trait.1,trait.2}}} & \sigma_{e_{\text{trait.2}}}^2 \end{bmatrix}\right),$$

where I is the identity matrix, σ_e_p^2 is the residual variance for purebreds, σ_e_c^2 is the residual variance for crossbreds in cages of 6, 7, or 8 grouped hens. A residual covariance between the two traits, because in this case the two traits were measured on the same individuals. Each cage size was analyzed separately (2 sizes each for WA and BD crossbreds), because of the different residual variances.
Variance component estimation was performed with the statistical software ASReml 4.2 [27].

## Heritability estimation

For purebred phenotypes, heritability estimates were calculated from the estimates of variance components obtained with the preceding models as:

$$h_p^2 = \frac{\sigma_{a_p}^2}{\sigma_{a_p}^2 + \sigma_{e_p}^2}$$

For crossbred phenotypes, the calculation of the heritability needs to account for (i) the fact that the phenotypes are the average per cage of the pooled performance from all the layers in the cage and (ii) for the half-sib relationships between the grouped layers. The formula to scale heritability estimated with pooled data to an individual-based heritability was established in previous studies [28–30], as shown in Additional file 2, and was as follows (applied separately to estimates of variance components for each cage size):

$$h_{c·n}^2 = \frac{4\sigma_{s_c}^2}{\sigma_{P_c^*}^2} \text{ with } \sigma_{P_c^*}^2 = \sigma_{s_c}^2 + n\sigma_{e_{c·n}}^2,$$

where *n* is the initial number of hens grouped in the cage (ranging from 6 to 8), $4\sigma_{s_c}^2$ is the additive genetic variance for crossbreds; and $\sigma_{e_c}^2$ is the residual variance for crossbreds in cages of 6, 7, or 8 grouped hens. However, because, mathematically, the group phenotype for LNVAR is not the mean of individual records of LNVAR (see the details of the transformations in Additional file 2), its heritability was estimated as:

$$h_{c·n}^2 = \frac{4\sigma_{s_c}^2}{\sigma_{s_c}^2 + \sigma_{e_{c·n}}^2}.$$

## Significance testing

Likelihood-ratio (LR) tests were performed to test if the estimates of heritability were statistically significantly different from 0 or if the genetic correlations were different from 0 or 1. The LR was calculated as:

$$LR = -2(\ell_{full} - \ell_{constrained}),$$

where $\ell_{full}$ is the likelihood of the univariate model for estimation of heritabilities and of the bivariate model for estimation of genetic correlations between traits, and $\ell_{constrained}$ is the likelihood of the corresponding constrained model. To test whether heritability estimates were statistically different from 0, the constrained models

were identical to the univariate full models, except that they were fitted without the random genetic effect. To test whether the genetic correlations were statistically significantly different from 0 or 1, the constrained models were identical to the bivariate full models, except that the genetic covariance was set to 0 or 1 (using the *rr*1 structure for the genetic effects in ASReml 4.2 [27]). Assuming that $LR \sim \chi_{(n)}^2$, where *n* is the reduction of the number of degrees of freedom between the compared models, the p-value can be retrieved from the chi-squared distribution. Note that when the comparison of models corresponds to a test of a parameter on the boundary of parameter space (test of heritability equal to 0 or correlation equal to 1), the distribution of this test statistic under the null hypothesis is a 50:50 mixture of $\chi_0^2$ and $\chi_1^2$ distributions (i.e. the threshold value for $\alpha=0.05$ is $\chi^2 = 2.706$) [31].

## Results

### Observed egg production and resilience indicator traits

The two purebred lines were high-producing layers. The observed mean EP (between 25 and 83 weeks of age) was 371 and 347 eggs in the WA and BD purebreds, respectively (Table 1). Compared to the purebreds, the crossbreds had a similar EP, a lower LNVAR (i.e. smaller the fluctuations), a higher SKEW (i.e. a less symmetric distribution of the deviations), and a comparable AUTO-R (a higher AUTO-R corresponds to two subsequent deviations being more alike). The standard deviations were large for all traits (phenotypic coefficient of variation of about 10% for most traits), although they were smaller in crossbreds than in purebreds. It is important to keep in mind that crossbreds were housed in groups while purebreds were housed individually and averages typically have a lower variance than individual records.

<table><caption>**Table 1** Phenotypic means and standard deviations (SD) of the four analyzed traits in purebreds and crossbreds from the WA and BD lines</caption><thead><tr><th rowspan="2">Traits</th><th rowspan="2">Type</th><th colspan="2">WA line</th><th colspan="2">BD line</th></tr><tr><th>Mean</th><th>SD</th><th>Mean</th><th>SD</th></tr></thead><tbody><tr><td rowspan="2">EP</td><td>Purebred</td><td>371</td><td>51</td><td>347</td><td>69</td></tr><tr><td>Crossbred</td><td>368</td><td>24</td><td>341</td><td>36</td></tr><tr><td rowspan="2">LNVAR</td><td>Purebred</td><td>-0.81</td><td>0.94</td><td>-0.33</td><td>1.09</td></tr><tr><td>Crossbred</td><td>-2.32</td><td>0.73</td><td>-1.84</td><td>0.64</td></tr><tr><td rowspan="2">SKEW</td><td>Purebred</td><td>-1.51</td><td>1.03</td><td>-1.20</td><td>1.15</td></tr><tr><td>Crossbred</td><td>-0.58</td><td>0.65</td><td>-0.40</td><td>0.63</td></tr><tr><td rowspan="2">AUTO-R</td><td>Purebred</td><td>0.21</td><td>0.24</td><td>0.36</td><td>0.29</td></tr><tr><td>Crossbred</td><td>0.19</td><td>0.27</td><td>0.34</td><td>0.23</td></tr></tbody></table>

EP total number of eggs laid per layer between 25 and 83 weeks of age, LNVAR natural logarithm of the variance of the deviations between the observed phenotype and the average of the batch, SKEW skewness of the distribution of these deviations, AUTO-R lag-one autocorrelation of these deviations
<table><caption>**Table 2** Estimates of genetic parameters with their standard error in parentheses for pure- and crossbred layers of the WA line</caption><thead><tr><th></th><th>EP</th><th>LNVAR</th><th>SKEW</th><th>AUTO-R</th></tr></thead><tbody><tr><td><strong>Genetic variance</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>σ<sup>2</sup><sub>d<sub>p</sub></sub></td><td>107 (10)</td><td>0.065 (0.006)</td><td>0.037 (0.006)</td><td>0.003 (&lt;0.001)</td></tr><tr><td>σ<sup>2</sup><sub>d<sub>c</sub></sub></td><td>91 (9)</td><td>0.046 (0.008)</td><td>0.037 (0.006)</td><td>0.003 (&lt;0.001)</td></tr><tr><td><strong>Residual variance</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>σ<sup>2</sup><sub>e<sub>p</sub></sub></td><td>873 (10)</td><td>0.609 (0.006)</td><td>0.920 (0.008)</td><td>0.047 (&lt;0.001)</td></tr><tr><td>σ<sup>2</sup><sub>e<sub>c-6</sub></sub></td><td>214 (4)</td><td>0.343 (0.007)</td><td>0.399 (0.008)</td><td>0.060 (0.001)</td></tr><tr><td>σ<sup>2</sup><sub>e<sub>c-7</sub></sub></td><td>242 (5)</td><td>0.336 (0.006)</td><td>0.382 (0.007)</td><td>0.065 (0.001)</td></tr><tr><td><strong>Total variance</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>σ<sup>2</sup><sub>P<sub>p</sub></sub></td><td>980 (9)</td><td>0.675 (0.006)</td><td>0.957 (0.008)</td><td>0.050 (&lt;0.001)</td></tr><tr><td>σ<sup>2</sup><sub>P<sub>c-6</sub></sub></td><td>1307 (25)</td><td>0.356 (0.007)</td><td>2.400 (0.045)</td><td>0.361 (0.007)</td></tr><tr><td>σ<sup>2</sup><sub>P<sub>c-7</sub></sub></td><td>1711 (31)</td><td>0.346 (0.006)</td><td>2.681 (0.048)</td><td>0.456 (0.008)</td></tr><tr><td><strong>Genetic parameters</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>h<sup>2</sup><sub>P</sub></td><td>0.11 (0.01)</td><td>0.10 (0.01)</td><td>0.04 (0.01)</td><td>0.06 (0.01)</td></tr><tr><td>h<sup>2</sup><sub>c-6</sub></td><td>0.07 (0.01)</td><td>0.13 (0.02)</td><td>0.006 (0.002)</td><td>0.02 (&lt;0.01)</td></tr><tr><td>h<sup>2</sup><sub>c-7</sub></td><td>0.05 (0.01)</td><td></td><td></td><td>0.01 (&lt;0.01)</td></tr><tr><td>r<sub>pc</sub></td><td>0.31 (0.10)</td><td>0.16 (0.13)</td><td>0.20 (0.23)</td><td>0.63 (0.14)</td></tr></tbody></table>

The variances reported are the variances estimated by the bivariate analyses with the same trait in purebreds and crossbreds analyzed jointly

EP total number of eggs laid per layer between 25 and 83 weeks of age, LNVAR natural logarithm of the variance of the deviations between the observed phenotype and the average of the batch, SKEW skewness of the distribution of these deviations, AUTO-R lag-one autocorrelation of these deviations

σ<sup>2</sup><sub>d<sub>p</sub></sub>: additive genetic variances for purebreds; σ<sup>2</sup><sub>d<sub>c</sub></sub>: additive genetic variances for crossbreds; σ<sup>2</sup><sub>e<sub>p</sub></sub>: residual variances for purebreds; σ<sup>2</sup><sub>e<sub>c</sub></sub>: residual variances for r crossbreds; where n is the number of hens in the cage; σ<sup>2</sup><sub>P<sub>p</sub></sub>: total variance for purebreds; σ<sup>2</sup><sub>P<sub>c-6</sub></sub>: total variance for crossbreds; h<sup>2</sup><sub>P</sub>: heritability estimates for purebreds; h<sup>2</sup><sub>c-6</sub>: heritability estimates for crossbreds; r<sub>pc</sub>: genetic correlation between purebred and crossbreds

## Genetic parameters of egg production and resilience indicator traits

For purebreds, EP were estimated to be moderately heritable (WA: 0.11 in Table 2; BD: 0.12 in Table 3), while resilience indicators traits were estimated to be low to moderately heritable, depending on the trait. Heritability estimates were equal to 0.10 and 0.12 for LNVAR, 0.04 and 0.02 for SKEW, and 0.06 and 0.08 for AUTO-R in the WA and BD lines, respectively.

For crossbred phenotypes, EP and resilience indicator traits were estimated to be less heritable than in purebreds, except for LNVAR. Although estimates of additive genetic variances were of the same magnitude for purebred and crossbred phenotypes, the adjusted phenotypic variances (σ<sup>2</sup><sub>P<sub>p</sub></sub>) were larger in crossbreds than the phenotypic variances in purebreds (σ<sup>2</sup><sub>P<sub>c</sub></sub>), except for LNVAR. Heritability estimates ranged from 0.04 to 0.07 for EP, from 0.13 to 0.21 for LNVAR and from 0.01 to 0.03 for AUTO-R (see Tables 2 and 3). Heritability estimates for SKEW were 0.006 and 0.005 in cages of six and seven WA crossbreds, respectively

<table><caption>**Table 3** Estimates of genetic parameters with their standard error in parentheses for pure- and crossbred layers of the BD line</caption><thead><tr><th></th><th>EP</th><th>LNVAR</th><th>SKEW</th><th>AUTO-R</th></tr></thead><tbody><tr><td><strong>Genetic variance</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>σ<sup>2</sup><sub>d<sub>p</sub></sub></td><td>325 (29)</td><td>0.137 (0.012)</td><td>0.027 (0.006)</td><td>0.006 (0.001)</td></tr><tr><td>σ<sup>2</sup><sub>d<sub>c</sub></sub></td><td>96 (19)</td><td>0.055 (0.015)</td><td>0.016 (0.013)</td><td>0.009 (0.002)</td></tr><tr><td><strong>Residual variance</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>σ<sup>2</sup><sub>e<sub>p</sub></sub></td><td>2400 (27)</td><td>1.022 (0.011)</td><td>1.236 (0.011)</td><td>0.074 (0.001)</td></tr><tr><td>σ<sup>2</sup><sub>e<sub>c-6</sub></sub></td><td>293 (9)</td><td>0.242 (0.008)</td><td>0.369 (0.011)</td><td>0.046 (0.002)</td></tr><tr><td>σ<sup>2</sup><sub>e<sub>c-7</sub></sub></td><td>269 (11)</td><td>0.357 (0.014)</td><td>0.382 (0.015)</td><td>0.049 (0.002)</td></tr><tr><td><strong>Total variance</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>σ<sup>2</sup><sub>P<sub>p</sub></sub></td><td>2725 (24)</td><td>1.159 (0.010)</td><td>1.264 (0.010)</td><td>0.08 (0.001)</td></tr><tr><td>σ<sup>2</sup><sub>P<sub>c-6</sub></sub></td><td>2078 (64)</td><td>0.256 (0.008)</td><td>2.586 (0.079)</td><td>0.327 (0.010)</td></tr><tr><td>σ<sup>2</sup><sub>P<sub>c-7</sub></sub></td><td>2173 (88)</td><td>0.371 (0.014)</td><td>3.063 (0.119)</td><td>0.396 (0.016)</td></tr><tr><td><strong>Genetic parameters</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>h<sup>2</sup><sub>P</sub></td><td>0.12 (0.01)</td><td>0.12 (0.01)</td><td>0.02 (&lt;0.01)</td><td>0.08 (0.01)</td></tr><tr><td>h<sup>2</sup><sub>c-6</sub></td><td>0.05 (0.01)</td><td>0.21 (0.06)</td><td>0.006<sup>a</sup> (0.005)</td><td>0.03 (0.01)</td></tr><tr><td>h<sup>2</sup><sub>c-7</sub></td><td>0.04 (0.01)</td><td>0.15 (0.04)</td><td>0.005<sup>a</sup> (0.004)</td><td>0.02 (0.01)</td></tr><tr><td>r<sub>pc</sub></td><td>0.52 (0.14)</td><td>0.47 (0.16)</td><td>0.37 (0.40)</td><td>0.56 (0.17)</td></tr></tbody></table>

The variances reported are the variances estimated by the bivariate analyses with the same trait in purebreds and crossbreds analyzed jointly

EP total number of eggs laid per layer between 25 and 83 weeks of age, LNVAR natural logarithm of the variance of the deviations between the observed phenotype and the average of the batch, SKEW skewness of the distribution of these deviations, AUTO-R lag-one autocorrelation of these deviations

σ<sup>2</sup><sub>d<sub>p</sub></sub>: additive genetic variances for purebreds; σ<sup>2</sup><sub>d<sub>c</sub></sub>: additive genetic variances for crossbreds; σ<sup>2</sup><sub>e<sub>p</sub></sub>: residual variances for purebreds; σ<sup>2</sup><sub>e<sub>c</sub></sub>: residual variances for r crossbreds; where n is the number of hens in the cage; σ<sup>2</sup><sub>P<sub>p</sub></sub>: total variance for purebreds; σ<sup>2</sup><sub>P<sub>c-6</sub></sub>: total variance for crossbreds; h<sup>2</sup><sub>P</sub>: heritability estimates for purebreds; h<sup>2</sup><sub>c-6</sub>: heritability estimates for crossbreds; r<sub>pc</sub>: genetic correlation between purebred and crossbreds

<sup>a</sup> Indicates h<sup>2</sup> not significantly different from 0 (see "Estimation of variance components for the same trait between purebreds and crossbreds" section about the LRT)

(Table 2), and were not statistically significantly different from 0 in BD crossbreds (Table 3). In summary, the resilience indicator traits showed low heritability estimates in purebreds and even lower estimates in crossbreds.

## Genetic correlations between the resilience traits

Phenotypic correlations are in Additional file 3: Tables S3–S5.

In purebreds, estimates of genetic correlations were moderately negative between LNVAR and SKEW, low between LNVAR and AUTO-R (negative or positive depending on the line), and moderately negative between SKEW and AUTO-R in both lines (Table 4).

In crossbreds, the estimate of the genetic correlation between LNVAR and SKEW was −0.67 in cages of eight BD crossbreds and was not statistically significantly different from 0 for the other crossbreds (Table 5).
<table><caption>**Table 4** Estimates of genetic correlations between traits in purebreds with their standard error in parentheses for the WA line (below the diagonal) and the BD line (above the diagonal)</caption><thead><tr><th>WA\BD</th><th>EP</th><th>LNVAR</th><th>SKEW</th><th>AUTO-R</th></tr></thead><tbody><tr><td>EP</td><td></td><td>-0.83 (0.04)</td><td>0.49 (0.15)</td><td>0.38 (0.11)</td></tr><tr><td>LNVAR</td><td>-0.72 (0.04)</td><td></td><td>-0.67 (0.11)</td><td>0.14 (0.10)</td></tr><tr><td>SKEW</td><td>-0.05 (0.09)</td><td>-0.46 (0.07)</td><td></td><td>-0.55 (0.15)</td></tr><tr><td>AUTO-R</td><td>0.25 (0.09)</td><td>-0.01 (0.08)</td><td>-0.21 (0.09)</td><td></td></tr></tbody></table>

EP total number of eggs laid per layer between 25 and 83 weeks of age, LNVAR natural logarithm of the variance of the deviations between the observed phenotype and the average of the batch, SKEW skewness of the distribution of these deviations, AUTO-R lag-one autocorrelation of these deviations

<table><caption>**Table 5** Estimates of genetic correlations between traits in the WA crossbreds with their standard error in parentheses, for cages of six (below the diagonal) and seven hens (above the diagonal)</caption><thead><tr><th>6\7 hens</th><th>EP</th><th>LNVAR</th><th>SKEW</th><th>AUTO-R</th></tr></thead><tbody><tr><td>EP</td><td></td><td>-0.82 (0.06)</td><td>0.19<sup>a</sup> (0.24)</td><td>-0.22<sup>a</sup> (0.13)</td></tr><tr><td>LNVAR</td><td>-0.74 (0.09)</td><td></td><td>-0.16<sup>a</sup> (0.23)</td><td>0.56 (0.10)</td></tr><tr><td>SKEW</td><td>-0.38 (0.18)</td><td>-0.18<sup>a</sup> (0.28)</td><td></td><td>0.46<sup>a</sup> (0.30)</td></tr><tr><td>AUTO-R</td><td>0.21 (0.18)</td><td>0.21<sup>a</sup> (0.25)</td><td>0.34<sup>a</sup> (0.39)</td><td></td></tr></tbody></table>

EP total number of eggs laid per layer between 25 and 83 weeks of age, LNVAR natural logarithm of the variance of the deviations between the observed phenotype and the average of the batch, SKEW skewness of the distribution of these deviations, AUTO-R lag-one autocorrelation of these deviations

<table><caption>**Table 6** Estimates of genetic correlations between traits in the BD crossbreds with their standard error in parentheses for cages of seven (below the diagonal) and eight hens (above the diagonal)</caption><thead><tr><th>7\8 hens</th><th>EP</th><th>LNVAR</th><th>SKEW</th><th>AUTO-R</th></tr></thead><tbody><tr><td>EP</td><td></td><td>-0.80<sup>b</sup> (0.13)</td><td>0.40 (0.27)</td><td>-0.39<sup>a</sup> (0.44)</td></tr><tr><td>LNVAR</td><td>-0.54 (0.20)</td><td></td><td>-0.67 (0.27)</td><td>0.83<sup>b</sup> (0.33)</td></tr><tr><td>SKEW</td><td>-0.34<sup>a</sup> (0.51)</td><td>-0.70<sup>ab</sup> (0.44)</td><td></td><td>-0.42<sup>a</sup> (0.67)</td></tr><tr><td>AUTO-R</td><td>0.14<sup>a</sup> (0.26)</td><td>0.49 (0.18)</td><td>-0.40<sup>a</sup> (0.45)</td><td></td></tr></tbody></table>

EP total number of eggs laid per layer between 25 and 83 weeks of age, LNVAR natural logarithm of the variance of the deviations between the observed phenotype and the average of the batch, SKEW skewness of the distribution of these deviations; AUTO-R lag-one autocorrelation of these deviations

Estimates of the genetic correlation between LNVAR and AUTO-R ranged from 0.21 to 0.83 and was not statistically significantly different from 0 in cages of six WA crossbreds (Table 6). None of the genetic correlations between SKEW and AUTO-R were statistically significantly different from 0.

In summary, our results indicate synergetic genetic correlations between the resilience indicator traits at least in the purebreds. This suggests that including one trait as a selection criterion will lead to genetic improvements of the other traits.

## Genetic correlations between egg production and resilience indicator traits

In both purebred lines, the genetic correlations were highly negative between EP and LNVAR, from lowly negative to moderately positive between EP and SKEW, and moderately positive between EP and AUTO-R (Table 4).

In the crossbreds, genetic correlation estimates were highly negative between EP and LNVAR, ranged from -0.38 to 0.40 between EP and SKEW (Tables 5 and 6). Although in the same direction as their respective estimates in the purebreds, the magnitude of the estimate of the genetic correlation between EP and SKEW was higher in the crossbreds. Estimates of the genetic correlations of EP with AUTO-R ranged from -0.39 to 0.21.

In summary, estimates of genetic correlations of EP with LNVAR were in all cases negative and therefore favorable. In contrast, estimates of genetic correlations of EP with SKEW and AUTO-R were sometimes positive and sometimes negative, but never very strong.

## Purebred-crossbred genetic correlations

Estimates of genetic correlations for the same trait between purebreds and crossbreds ($r_{pc}$) ranged from 0.16 to 0.63 across traits (Tables 2 and 3). Generally speaking, the corresponding estimates of $r_{pc}$ were higher in the BD than in the WA line. All estimates of $r_{pc}$ in this study were statistically significantly different from 1, except for the estimate for SKEW, which had large standard errors. This indicates that the traits in purebreds and crossbreds (EP and the resilience indicator traits) are genetically different traits but do share part of the same genetic background.

## Discussion

Our goal was to study general resilience of layer chickens, which is expected to be related to the capacity of an animal to cope with different types of challenges: disease, climate change, negative social interactions, etc. [5, 6, 10]. Thus, we explored the genetic background of LNVAR, AUTO-R, and SKEW of the deviations of weekly egg production of an individual hen from batch average to assess their potential as selection criteria for resilience in layers. Our results indicate that these resilience indicator traits are lowly heritable and are genetically different in purebred and crossbred laying hens. Moreover, LNVAR showed favorable genetic correlations with egg production.

<footnote><sup>a</sup> Indicates that the estimate of the genetic correlation was not significantly different from 0 (see “Methods” section on LRT)</footnote>

<footnote><sup>a</sup> Indicates that the estimate of the genetic correlation was not significantly different from 0 (see “Methods” section on LRT)</footnote>

<footnote><sup>b</sup> Indicates that the estimate of the genetic correlation was not significantly different from 1 (see “Methods” section on LRT)</footnote>
## Egg production

Our estimates of heritability for EP in purebreds were low and consistent with those from previous studies that were performed on EP data for similar time frames and environments (ranging from 0.01 to 0.20 across the laying period [32, 33]). In general, reported estimates of heritability for egg production traits are higher for purebreds than for crossbreds [34, 35], which appears unexpected because part of the dominance variance being captured in the additive genetic variance, hence heritability is usually higher in crossbreds than in purebreds [36]. However, estimates of residual variances vary greatly between studies and are mostly higher in crossbreds than in purebreds, which may be due to environmental conditions of crossbreds varying more than those of purebreds and makes comparisons complicated. Because purebred EP is recorded in individual cages and crossbred EP is recorded in group cages, social interactions may increase phenotypic variation for crossbred EP. Differences in heritability estimates between individual and group records have also been reported in the literature [28, 29, 37]. These differences can be explained by differences in the estimation of variance components, either because of incomplete information on the relationships between individuals due to the use of group information instead of individual information (in purebreds the complete pedigree is known while in purebred only the sire branch is known), or modelling choices (e.g. taking social interactions into account or not). Yet, when using individual data to simulate group records, estimates of the variance components based on individual and group records are the same [30].

## Using LNVAR and AUTO-R to improve resilience to disturbances

LNVAR is expected to be an indicator of both the severity and of the duration of a perturbation that the animal is experiencing, i.e. animals with a small variance are less affected by a perturbation than those with a large variance [6]. In our study, LNVAR was the most promising resilience selection criterion, because it had the highest heritability and also had a favorable genetic correlation with EP.

Our estimates of heritability for LNVAR of purebred EP (estimates of 0.11 in WA and 0.12 in BD) were similar to the estimate of heritability for LNVAR of body weight (repeated records) of layers related to the WA line (estimate of 0.10) [18]. Few studies have reported genetic parameters for LNVAR for multiple observations on one individual. LNVAR is by definition the natural logarithm of the variance of the residuals from an animal model of a trait [6]. In layers, heritabilities of 0.01 have been reported for residual variation of eggshell color at the individual egg level and of 0.15 at the level

of 15 to 20 eggs measured per hen [20], and of 0.10 for residual variation of egg weight (also repeated records) [19]. In broilers, low heritabilities (lower than 0.05) have been reported for residual variation of body weight (one record per individual) [38–40]. Other studies in dairy cattle, pigs, and aquatic species, with either one record per individual or repeated records have also reported mostly low heritabilities for residual variance (ranging from 0.01 to 0.14; [16, 41–45]).

In our study, the genetic correlation of LNVAR with EP was favorable and negative (ranging from −0.83 to −0.54). In other species, LNVAR has also been shown to have favorable genetic correlations with other traits. In pigs, lower day-to-day variance in feed intake (root mean square error of prediction of feed intake) was genetically correlated with lower mortality and fewer health treatments [17]. In dairy cattle, LNVAR of milk yield deviations was found to be favorably genetically correlated with health (udder and hoof), longevity (productive longevity, calf and maternal survival at birth), fertility (calving interval, first to last service interval), metabolic status (ketosis resistance), and production traits (milk yield, persistency, body condition score, and feed intake) [15, 16].

To our knowledge, our study is the first to report genetic correlations of LNVAR with EP in layers. These correlations imply that selecting for a reduced LNVAR would increase EP, which is desirable. EP has been a selection criterion for many generations and, therefore, has probably contributed to indirect genetic improvement of LNVAR. Given the low heritability of EP, adding LNVAR in the breeding goal is an opportunity to improve not only resilience but also production, in spite of the low heritability of LNVAR. The next step towards the inclusion of LNVAR in a breeding program is to derive economic values or to set the targeted genetic gain. A simplified simulation study has already shown that selecting for resilience indicators has a beneficial economic impact, and even more if health traits are not included in the breeding program, because of reduced labor and treatments [6].

In summary, although selecting layers on EP already improves LNVAR, including LNVAR in breeding programs may further enhance the improvement of resilience. The absence of trade-offs makes it relatively easy to increase resilience without a loss in selection response in EP. In addition, literature on other species shows that LNVAR is favorably associated with general immunity and health [6, 45].

The resilience measure AUTO-R is expected to be informative regarding the duration of a perturbation and thus to be an indicator of an animal's recovery capacity: close to 0 (subsequent deviations of EP are uncorrelated)
for animals with high recovery rates, close to +1 (subsequent deviations of EP are correlated) for animals with a slow recovery rate; and close to -1 (subsequent deviations of EP are opposite) for animals showing a compensatory response to the perturbation [6]. In our study, AUTO-R had a low heritability estimate in both purebreds and in crossbreds. These heritability estimates are similar to the heritability for AUTO-R of body weight deviations in layers (value of 0.09) [18]. Similarly, in dairy cattle, the heritability for AUTO-R of milk yield deviations was estimated at 0.09 [16]. In our study, AUTO-R was estimated to be unfavorably genetically correlated with EP in purebreds, but this relationship was less clear in crossbreds. Other studies on layers have not shown a clear genetic relationship of AUTO-R with health: a non-significant genetic correlation between AUTO-R of body weight deviations and titers of natural antibody (NAbs) isotypes in serum was reported [18]. However, in dairy cattle, AUTO-R of milk yield deviations was found to be favorably genetically correlated with health, longevity, fertility, metabolic, and production traits, but at a lower level than LNVAR [16].

Before using AUTO-R as a selection criterion, its definition as a trait needs to be further discussed to reach a consensus, and its biological impact on production traits and its relationship with health needs to be understood. Knowing both the LNVAR and AUTO-R of an animal can lead to complicated interpretations: phenotypically, animals with a low LNVAR could display an AUTO-R close to 1 because they are hardly affected by perturbations, and their subsequent deviations from EP are more alike. In the present study, we found that LNVAR and AUTO-R were favorably genetically correlated in all populations but their genetic correlations were less than 1. However, if layers have genetically a high AUTO-R when faced with a severe perturbation, they will show a slow recovery. Recovery from a perturbation may be an interesting trait in the case of severe perturbations (such as a virus outbreak or a heatwave), thus it may be more appropriate to calculate AUTO-R during periods with large deviations instead of over the whole period.

In summary, selecting purebred layers on EP is expected to increase AUTO-R, which could decrease recovery rates from perturbation. Including AUTO-R in breeding programs may further enhance the improvement of resilience because it would put more emphasis on recovery rate. The trade-off between AUTO-R and EP also highlights the need to quantify economic values when considering the inclusion of AUTO-R in breeding programs.

SKEW is expected to be an indicator of both the direction (positive or negative) and severity of a perturbation that an animal is experiencing, i.e. animals with a negative

SKEW experience more negative deviations than others, which is interpreted as poor resilience [6]. However, heritability estimates for SKEW were very low in the purebreds and crossbreds from both lines, similar to previous studies in layers [18] and in cattle [16], which reported that SKEW of body weight deviations and of milk yield deviations are not heritable. SKEW is known to be highly sensitive to outliers and thus to contain more noise [6]. We conclude that SKEW does not seem to be a useful resilience indicator trait when breeding for resilient laying hens. Overall, estimates of the genetic correlation of SKEW with LNVAR and AUTO-R were favorable, which means that genetic selection on LNVAR or AUTO-R is not expected to result in degraded SKEW.

## Resilience indicators in purebreds and crossbreds are genetically different

In our study, all $r_{pc}$ estimates were statistically significantly lower than 1, ranging from 0.16 to 0.63. This suggests that EP and resilience indicators traits are genetically different between purebreds and crossbreds. Our estimates of $r_{pc}$ were rather low compared to those in the literature, which range from 0.45 to 0.87 in pigs [21] and poultry [46]. Three main factors can contribute to $r_{pc}$ being lower than 1: (i) genotype-by-genotype interactions ($G \times G$), (ii) genotype-by-environment ($G \times E$) interactions, and (iii) differences in measurement or definition [20, 21].

First, crossbreds only share half of the genetic characteristics of each purebred parental line. Allele frequencies likely differ between purebreds and crossbreds and such differences may have direct effects on the trait (causal variants) or indirect effects (also known as $G \times G$: dominance or epistasis), depending on the loci involved.

Second, $G \times E$ interactions can contribute to the same trait having a different genetic architecture in different environments. Purebreds and crossbreds were raised in very different environments in our study in terms of climate, management, feed, barns, disease pressure, etc. The farms where purebreds were kept have a higher level of biosecurity than regular farms because the breeding companies cannot take any risk with the population under selection. RT-farms generally belong to farmers who collaborate with breeding companies and, although their health and care practices are usually among the best ones, the layers will experience more challenges on an RT-farm than on a nucleus farm. Commercial farm environments can still be very different from RT-farm environments (cage-free systems, winter garden access, outdoor access, etc.). This diversity of environments and their susceptibility to changes (heat waves, disease outbreak, etc.) may influence the resilience performance of hens due to differences in types, severity, and duration
of the perturbations. In addition, purebreds were housed individually, while crossbreds were in group cages, in which case the presence of cage mates creates positive (e.g. grooming) as well as negative (e.g. pecking) social interactions that can affect the individuals' wellbeing and performance. As a consequence of G × E interactions, re-ranking of sires may occur. In both lines, estimates of r<sub>pc</sub> were lower for LNVAR than for EP or AUTO-R, which may be due to G × E interactions having a greater role in traits based on variances than in traits based on means. In trout, G × E interactions were reported to be larger for the variance of body weight than for body weight itself [47].

Third, differences in the definition or measurement of the traits may contribute to r<sub>pc</sub> being lower than 1. Although the definition of the variable may be the same, values may differ due to differences in the recording machines and protocols used, (e.g. housing conditions such as individual vs. group), the commitment of the collaborators, units, errors, etc., that affect the quality of the raw data. These differences can contribute to reducing the r<sub>pc</sub>. Although the recording of the performance was the same between purebreds and crossbreds, dealing with records from individually-housed hens or group-housed hens may lead to different traits (e.g. implication of social interactions) and thus result in reduced r<sub>pc</sub> (theoretical derivation [20]).

All these elements can contribute to the value of recording the traits in crossbreds in a commercial environment. Breeding programs would benefit from considering both pure- and crossbred data to improve response to selection for crossbred performance in commercial environments [23, 48, 49], because it can increase the accuracy of the estimated breeding values (EBV) of purebred animals for crossbred performance.

## Longitudinal recording to genetically improve general resilience in farm animals

Our study shows that longitudinal recording offers opportunities to develop resilience indicators that are heritable. Here, resilience indicator traits are based on deviations of production performance from the undisturbed phenotype (here defined as the contemporary group average). An alternative is to estimate the undisturbed phenotype of an individual based on various modeling approaches [14, 16, 50]. Both approaches may underestimate the undisturbed phenotype, leading to underestimation of deviations and, thus, overestimation of resilience. On the one hand, if a perturbation affects the majority of the individuals in the contemporary group, the contemporary group average will be lower than the true undisturbed phenotype. On the other hand, if a perturbation affects the phenotype of

an individual permanently (i.e. poor recovery), the individual's trajectory based on observations will be lower than the undisturbed phenotype. Another approach is to estimate breeding values that allow to calculate an undisturbed trajectory (e.g. breeding values for curve parameters) [6]. Nevertheless, when modeled, a correction for the contemporary group effect is needed to account for the difference in disturbance level between contemporary groups. More research is needed to provide insight into the genetic correlations between resilience indicator traits estimated with all three approaches and which approach would best predict resilience [51].

General resilience is likely a combination of resilience of various biological functions, including food intake, digestion, growth, maintenance, reproduction, thermoregulation, immunity, behavior, etc. [5, 6, 10]. In our study, LNVAR and AUTO-R of deviations of weekly EP may capture the resilience of (re)production, i.e. the capacity of EP to be minimally affected by disturbances or to rapidly return to the undisturbed performance. This could lead to resilience of EP but not necessarily to general resilience. Thus, future studies should jointly study LNVAR of EP, LNVAR of bodyweight, and LNVAR of feed intake (e.g. cross-correlations of deviations). Longitudinal data on different traits (egg weight, body weight, and feed intake) could also be used to identify periods with a clear perturbation at the herd level [51], which would enable the study of the capacity to recover from a perturbation (i.e. AUTO-R or slope of reaction norm). Resilience indicators based on different traits could also be combined in a selection index for general resilience.

Tools are available now for phenotype recording in layers, such as electronic nests, feeders, weighing scales, cameras, microphones, RFID chips (radio-frequency identification), other embedded sensors, NIRS or MIRS (near- or mid-infrared spectroscopy), etc., which allow the concomitant high-throughput phenotyping of different traits. Based on our results and those reported in other livestock species that show the potential of resilience indicator traits for genetic improvement of resilience, especially LNVAR, further investigations in chickens are needed (i) to find the best (combination of) production traits to calculate resilience indicators, and (ii) to record available perturbation-related indicators and relate them with the resilience indicator traits.

## Conclusions

This study shows that LNVAR and AUTO-R of weekly EP are lowly heritable in purebreds as well as in crossbreds, while SKEW had very low heritability estimates. We found that these three resilience indicator traits have favorable genetic correlations with each other and with EP. Compared to purebreds, EP and the resilience
indicator traits show lower heritability estimates in crossbreds, although estimates of genetic additive variance were roughly at the same level in crossbreds and purebreds. Studies on LNVAR can go forward towards its inclusion in breeding programs, i.e. by estimating economic values, and evaluating the economic and environmental impacts. AUTO-R may add complementary information, but both its definition and its joint inclusion with LNVAR in breeding programs need to be further studied. Estimates of genetic correlations between purebreds and crossbreds of production and resilience indicator traits were low in both lines, which suggests that resilience indicator traits in crossbreds are genetically different from resilience indicator traits in purebreds, which may be partly due to individual recording in the purebreds and group recording in the crossbreds. This implies that strategies that include information on crossbreds in breeding programs would be beneficial to the genetic improvement of resilience in crossbreds.

## Supplementary Information

The online version contains supplementary material available at https://doi.org/10.1186/s12711-022-00716-8.

**Additional file 1: Table S1.** Distribution of hens from the WA line according to country and farms. **Table S2.** Distribution of hens from the BD line according to country and farms.

**Additional file 2.** Heritability estimated with pooled records of half-sibs, for regular and resilience indicator traits [28, 30, 52].

**Additional file 3: Table S3.** Phenotypic correlations among traits in purebreds, and their standard error below in italic characters, for the WA line (below the diagonal) and the BD line (above the diagonal). **Table S4.** Phenotypic correlations among traits in the WA crossbreds, and their standard error below in italic characters, for cages of 6 hens (below the diagonal) and cages of 7 hens (above the diagonal). **Table S5.** Phenotypic correlations among traits in the BD crossbreds, and their standard error below in italic characters, for cages of 7 hens (below the diagonal) and cages of 8 hens (above the diagonal).

### Acknowledgements

The authors are grateful to WUR colleagues: Marieke Poppe for discussions on resilience indicators; Piter Bijma and Henk Bovenhuis for discussions on estimating heritability with pooled records; and Mario Calus and Joas Bos for sharing their unpublished literature review on genetic correlations between purebred and crossbred broilers and layers.

### Author contributions

NB, TB, MHPVDL, HM designed the research. KP and JV provided the data and participated in discussions about analyses and results. ID derived the equations for scaling heritability based on pooled records. NB and TB analyzed the data. NB wrote the manuscript. All authors participated in discussions about the interpretation of the results, and revised the manuscript. All authors read and approved the final manuscript.

### Funding

This work is financially supported by the Netherlands Organization for Scientific Research Earth and Life Sciences (NWO-ALW; ALWSA.2016.4). This research was made possible by funding from SusAn, an ERA-Net co-funded under the European Union's Horizon 2020 research and innovation program (www.era-susan.eu), under Grant Agreement no. 696231. This work is in kind supported by Hendrix Genetics. NB was visiting WUR for a stay of 3 months

and then 1 week per month during a year with the financial support of the Animal Genetics Division of INRAE, and the financial support of the graduate school Wageningen Institute of Animal Sciences (WIAS). TB is supported by the Alexander von Humboldt Foundation during part of this study.

### Availability of data and materials

The data that support the findings of this study are available from Hendrix Genetics but restrictions apply to the availability of these data, which were used under license for the current study, and thus are not publicly available. Data are however available from the authors upon reasonable request and with permission of Hendrix Genetics (katrijn.peeters@hendrix-genetics.com).

## Declarations

### Ethics approval and consent to participate

Not applicable.

### Consent for publication

Not applicable.

### Competing interests

NB, TB, MHPVDL, ID and HM declare that they have no competing interests. KP and JV are employed by Hendrix Genetics. The datasets are of interest to commercial targets of Hendrix Genetics, but this interest did not influence the results in this manuscript in any matter.

### Author details

¹PEGASE, INRAE, Institut Agro, 35590 Saint Gilles, France. ²Wageningen University & Research Animal Breeding & Genomics, P.O. Box 338, 6700 AH Wageningen, The Netherlands. ³Reproductive Biotechnology, TUM School of Life Sciences, Technical University of Munich, Liesel-Beckmann-Strasse 1, 85354 Freising, Germany. ⁴Hendrix Genetics B.V., P.O. Box 114, 5830 AC Boxmeer, The Netherlands. ⁵Université Paris-Saclay, INRAE, AgroParisTech, GABI, 78350 Jouy-en-Josas, France. ⁶Gen PhySE, Université de Toulouse, INRAE, ENVT, 31320 Castanet Tolosan, France.

Received: 19 March 2021 Accepted: 28 March 2022

Published online: 20 April 2022

## References

1. Rauw WM, Kanis E, Noordhuizen-Stassen EN, Grommers FJ. Undesirable side effects of selection for high production efficiency in farm animals: a review. Livest Prod Sci. 1998;56:15–33.

2. Rauw WM. Immune response from a resource allocation perspective. Front Genet. 2012;3:267.

3. Rendel JM. Correlation between the number of scutellar and abdominal bristles in *Drosophila melanogaster*. Genetics. 1963;48:391–408.

4. Dumont B, Fortun-Lamoth E, Jouven M, Thomas M, Tichit M. Prospects from agroecology and industrial ecology for animal production in the 21st century. Animal. 2013;7:1028–43.

5. Colditz IG, Hine BC. Resilience in farm animals: biology, management, breeding and implications for animal welfare. Anim Prod Sci. 2016;56:1961–83.

6. Berghof TVL, Poppe M, Mulder HA. Opportunities to improve resilience in animal breeding programs. Front Genet. 2019;9:692.

7. Mulder HA. Is GxE a burden or a blessing? Opportunities for genomic selection and big data. J Anim Breed Genet. 2017;134:435–6.

8. Scheffer M, Bolhuis JE, Borsboom D, Buchman TG, Gijzel SMW, Goulson D, et al. Quantifying resilience of humans and other animals. Proc Natl Acad Sci USA. 2018;115:11883–90.

9. Doeschl-Wilson AB, Villanueva B, Kyriazakis I. The first step toward genetic selection for host tolerance to infectious pathogens: obtaining the tolerance phenotype through group estimates. Front Genet. 2012;3:265.

10. Friggens NC, Blanc F, Berry DP, Puillet L. Review: Deciphering animal robustness. A synthesis to facilitate its use in livestock breeding and management. Animal. 2017;11:2237–51.

11. de Haas Y, Veerkamp RF, Barkema HW, Gröhn YT, Schukken YH. Associations between pathogen-specific cases of clinical mastitis and somatic cell count patterns. J Dairy Sci. 2004;87:95–105.
12. Codrea MC, Højsgaard S, Friggens NC. Differential smoothing of time-series measurements to identify disturbances in performance and quantify animal response characteristics: an example using milk yield profiles in dairy cows. *J Anim Sci.* 2011;89:3089–98.

13. Byrd CJ, Johnson JS, Radcliffe JS, Craig BA, Eicher SD, Lay DC. Nonlinear analysis of heart rate variability for evaluating the growing pig stress response to an acute heat episode. *Animal*. 2019;14:379–87.

14. Nguyen-Ba H, van Milgen J, Taghipoor M. A procedure to quantify the feed intake response of growing pigs to perturbations. *Animal*. 2020;14:253–60.

15. Elgersma GG, de Jong G, van der Linde R, Mulder HA. Fluctuations in milk yield are heritable and can be used as a resilience indicator to breed healthy cows. *J Dairy Sci.* 2018;101:1240–50.

16. Poppe M, Veerkamp RF, van Pelt ML, Mulder HA. Exploration of variance, autocorrelation, and skewness of deviations from lactation curves as resilience indicators for breeding. *J Dairy Sci.* 2020;103:1667–84.

17. Putz AM, Harding JCS, Dyck MK, Fortin F, Plastow GS, Dekkers JCM, et al. Novel resilience phenotypes using feed intake data from a natural disease challenge model in wean-to-finish pigs. *Front Genet.* 2019;9:660.

18. Berghof TVL, Bovenhuis H, Mulder HA. Body weight deviations as indicator for resilience in layer chickens. *Front Genet.* 2019;10:1216.

19. Wolc A, Arango J, Settar P, Fulton JE, O'Sullivan NP, Preisinger R, et al. Genome-wide association analysis and genetic architecture of egg weight and egg uniformity in layer chickens: genome-wide association analysis and genetic architecture of egg. *Anim Genet.* 2012;43:S87–96.

20. Mulder HA, Visscher J, Fablet J. Estimating the purebred–crossbred genetic correlation for uniformity of eggshell color in laying hens. *Genet Sel Evol.* 2016;48:39.

21. Wientjes YCJ, Calus MPL. Board Invited Review: The purebred–crossbred correlation in pigs: a review of theory, estimates, and implications. *J Anim Sci.* 2017;95:3467–78.

22. Duenk P, Bijma P, Wientjes YCJ, Calus MPL. Review: optimizing genomic selection for crossbred performance by model improvement and data collection. *J Anim Sci.* 2021;99: skab205.

23. Wei M, van der Werf HJH. Maximizing genetic response in cross-breds using both purebred and crossbred information. *Anim Sci.* 1994;59:401–13.

24. Bijma P, Woollams JA, van Arendonk JAM. Genetic gain of pure line selection and combined crossbred purebred selection with constrained inbreeding. *Anim Sci.* 2001;72:225–32.

25. R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2020.

26. Meyer D, Dimitriadou E, Hornik K, Weingessel A, Leisch F, Chang C-C, et al. e1071: Misc functions of the Department of Statistics, Probability Theory Group. 2020. https://CRAN.R-project.org/package=e1071/. Accessed 4 Mar 2021.

27. Gilmour AR, Gogel BJ, Cullis BR, Welham SJ, Thompson R. ASReml user guide release 4.2 functional specification. Hemel Hempstead: VSN International Ltd; 2021.

28. Biscarini F, Bovenhuis H, van Arendonk JAM. Estimation of variance components and prediction of breeding values using pooled data. *J Anim Sci.* 2008;86:2845–52.

29. Peeters K, Ellen ED, Bijma P. Using pooled data to estimate variance components and breeding values for traits affected by social interactions. *Genet Sel Evol.* 2013;45:27.

30. Peeters K. Genetics of social interactions in laying hens: improving survival and productivity Ph.D. thesis, Wageningen University; 2015.

31. Stram DO, Lee JW. Variance components testing in the longitudinal mixed effects model. *Biometrics*. 1994;50:1171–7.

32. Wolc A, Arango J, Settar P, O'Sullivan NP, Dekkers JCM. Evaluation of egg production in layers using random regression models. *Poult Sci.* 2011;90:30–4.

33. Liu Z, Yang N, Yan Y, Li G, Liu A, Wu G, et al. Genome-wide association analysis of egg production performance in chickens across the whole laying period. *BMC Genet.* 2019;20:67.

34. Besbes B, Ducrocq V, Foulley JL, Protais M, Tavernier A, Tixier-Boichard M, et al. Estimation of genetic parameters of egg production traits of laying hens by restricted maximum likelihood applied to a multiple-trait reduced animal model. *Genet Sel Evol.* 1992;24:539–52.

35. Wei M. Combined crossbred and purebred selection in animal breeding. Ph.D. thesis, Wageningen University; 1992.

36. Falconer DS, Mackay TFC. Introduction to quantitative genetics. 4th ed. Harlow: Longmans Green; 1996.

37. Nurgiartiningsih V, Mielenz N, Preisinger R, Schmutz M, Schüler L. Genetic parameters for egg production and egg weight of laying hens housed in single and group cages. *Arch Anim Breed.* 2002;45:501–8.

38. Rowe SJ, White IM, Avendaño S, Hill WG. Genetic heterogeneity of residual variance in broiler chickens. *Genet Sel Evol.* 2006;38:617–35.

39. Mulder HA, Hill WG, Vereijken A, Veerkamp RF. Estimation of genetic variation in residual variance in female and male broiler chickens. *Animal*. 2009;3:1673–80.

40. Wolc A, White IMS, Avendano S, Hill WG. Genetic variability in residual variation of body weight and conformation scores in broiler chickens. *Poult Sci.* 2009;88:1156–61.

41. Rönnegård L, Felleki M, Fikse WF, Mulder HA, Strandberg E. Variance component and breeding value estimation for genetic heterogeneity of residual variance in Swedish Holstein dairy cattle. *J Dairy Sci.* 2013;96:2627–36.

42. Vandenplas J, Bastin C, Engler N, Mulder HA. Genetic variance in micro-environmental sensitivity for milk and milk quality in Walloon Holstein cattle. *J Dairy Sci.* 2013;96:5977–90.

43. Sell-Kubiak E, Bijma P, Knol EF, Mulder HA. Comparison of methods to study uniformity of traits: application to birth weight in pigs. *J Anim Sci.* 2015;93:900–11.

44. lung LHS, Neves HHR, Mulder HA, Carvalheiro R. Genetic control of residual variance of yearling weight in Nellore beef cattle. *J Anim Sci.* 2017;95:1425–33.

45. lung LHDS, Carvalheiro R, Neves HHDR, Mulder HA. Genetics and genomics of uniformity and resilience in livestock and aquaculture species: a review. *J Anim Breed Genet.* 2020;137:263–80.

46. Calus MPL, Bos JA, Duenk P, Wientjes YCJ. The purebred–crossbred correlation in broilers and layers: a review. In: *Proceedings of the 71st annual meeting of the European Federation of Animal Science (EAAP)*, 1–4 December 2020; Online virtual meeting; 2020.

47. Sae-Lim P, Kause A, Janhunen M, Vehviläinen H, Koskinen H, Gjerde B, et al. Genetic (co)variance of rainbow trout (*Oncorhynchus mykiss*) body weight and its uniformity across production environments. *Genet Sel Evol.* 2015;47:46.

48. Bijma P, van Arendonk JAM. Maximizing genetic gain for the sire line of a crossbreeding scheme utilizing both purebred and crossbred information. *Anim Sci.* 1998;66:529–42.

49. Dekkers JCM. Marker-assisted selection for commercial crossbred performance. *J Anim Sci.* 2007;85:2104–14.

50. Revilla M, Guillaume L, Loïc F-G, Rafael M-T, Friggens NC. Quantifying growth perturbations over the fattening period in swine via mathematical modelling. *Peer Community J.* 2022;2:e9.

51. Poppe M, Mulder HA, Kamphuis C, Veerkamp RF. Between-herd variation in resilience and relations to herd performance. *J Dairy Sci.* 2021;104:616–27.

52. Eriksson M. A simulation method for skewness correction. *Report Uppsala University No. U.U.D.M.* project report 2008:22; 2008. p. 40. https://www.diva-portal.org/smash/get/diva2:302313/FULLTEXT01.pdf. Accessed 18 Feb 2022.

## Publisher's Note

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

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

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