- Research Article
- Open access
- Published:
Breeding for resilience in finishing pigs can decrease tail biting, lameness and mortality
Genetics Selection Evolution volume 56, Article number: 48 (2024)
Abstract
Background
Previous research showed that deviations in longitudinal data are heritable and can be used as a proxy for pigs’ general resilience. However, only a few studies investigated the relationship between these resilience traits and other traits related to resilience and welfare. Therefore, this study investigated the relationship between resilience traits derived from deviations in longitudinal data and traits related to animal resilience, health and welfare, such as tail and ear biting wounds, lameness and mortality.
Results
In our experiment, 1919 finishing pigs with known pedigree (133 Piétrain sires and 266 crossbred dams) were weighed every 2 weeks and scored for physical abnormalities, such as lameness and ear and tail biting wounds (17,066 records). Resilience was assessed via deviations in body weight, deviations in weighing order and deviations in observed activity during weighing. The association between these resilience traits and physical abnormality traits was investigated and genetic parameters were estimated. Deviations in body weight had moderate heritability estimates (h2 = 25.2 to 36.3%), whereas deviations in weighing order (h2 = 4.2%) and deviations in activity during weighing (h2 = 12.0%) had low heritability estimates. Moreover, deviations in body weight were positively associated and genetically correlated with tail biting wounds (rg = 0.22 to 0.30), lameness (rg = 0.15 to 0.31) and mortality (rg = 0.19 to 0.33). These results indicate that events of tail biting, lameness and mortality are associated with deviations in pigs’ body weight evolution. This relationship was not found for deviations in weighing order and activity during weighing. Furthermore, individual body weight deviations were positively correlated with uniformity at the pen level, providing evidence that breeding for these resilience traits might increase both pigs’ resilience and within-family uniformity.
Conclusions
In summary, our findings show that breeding for resilience traits based on deviations in longitudinal weight data can decrease pigs’ tail biting wounds, lameness and mortality while improving uniformity at the pen level. These findings are valuable for pig breeders, as they offer evidence that these resilience traits are an indication of animals’ general health, welfare and resilience. Moreover, these results will stimulate the quantification of resilience via longitudinal body weights in other species.
Background
Resilience is a multifaceted concept overlapping with other related concepts such as robustness, tolerance, and resistance [1]. In this study, we define resilience based on Berghof et al. [2], building on Colditz and Hine [3], as “the capacity of the animal to be minimally affected by disturbances or to rapidly return to the state pertained before exposure to a disturbance”. This definition implies that disturbances can affect animals, albeit with individual variation in the magnitude of the impact exerted by these disturbances. As a result, deviations in longitudinal data can be used as a proxy for overall resilience as more resilient animals will show less deviations from their optimal production level with a quick recovery, after a challenge [4]. Testing this hypothesis of resilience necessitates specific conditions and assumptions. First, the animals tested must encounter at least one challenge since evaluating resilience in an ‘optimal’ environment is unfeasible. Second, collection of longitudinal data is needed, as it is impossible to quantify temporal variability within an individual from a single data point [2]. These longitudinal data may encompass various aspects such as (re)production traits (e.g., weight or feed intake), immunological traits (e.g., viral load, antibody levels), or physiological traits (e.g., body temperature, heart rate) [2]. However, to quantify resilience effectively, these traits must be impacted by the challenges under consideration. And third, inferring the optimal (expected) performance of an animal from longitudinal data necessitates proper statistical modeling. Only then can deviations from this “optimal” performance be quantified.
Breeding resilient livestock is becoming increasingly important [4, 5]. The increasing intensification in pig farming, especially in the European Union (EU), results in more pigs to be managed per farmer [6]. As a result, improving pigs’ resilience is desired as it would reduce costs of labor and treatment while enhancing animal welfare and maintaining optimal production levels [5, 7, 8]. In addition, the goal of the EU Green Deals’ farm to fork strategy is to reduce antimicrobials by 50% in livestock and aquaculture by 2030 [9]. Besides, the focus on animal welfare is increasing in the EU, framed by a stringent welfare legislation. A recent example is the European citizens’ initiative ‘End the cage age’ [10], which could potentially lead to the prohibition of all cage systems in pig husbandry. The call to raise pigs in a more natural setting implies exposing them to a broader spectrum of challenging environments. In addition, the escalating impact of climate change is resulting in a significant rise in the frequency of heat waves [11, 12]. Finally, there is mounting pressure to incorporate more by-products into pig feed to alleviate food-feed competition. However, augmenting the proportion of by-products in pigs’ diets may present a dietary challenge, certainly since most breeding pigs are selected on high-quality, nutrient rich diets [13].
Scheffer et al. [4] indicated that a holistic approach might be optimal for improving general resilience. Indeed, previous studies have shown that breeding against specific diseases can be effective, but might increase susceptibility to other diseases [1, 14] and as such, not improve general resilience. Moreover, resilience indicators based on deviations from longitudinal data are more informative, as they reflect more continuous variation in phenotypical differences compared to classical scoring of resilience indicators, such as mortality (”dead” vs “alive”) or disease status (”not”, ”mildly” or “severely” affected) [15, 16]. For these reasons, Scheffer et al. [4] suggested that these new resilience traits based on deviations in longitudinal data could have substantial impact.
In the past, collecting longitudinal data was often practically challenging and expensive. However, due to technological developments the collection of longitudinal data will probably become standard in animal breeding [4]. Automatic feeding stations can already record individual pigs’ body weight, feed intake and feeding behaviour at every visit. Moreover, computer vision and wearable devices can phenotype pigs almost continuously [4], for example to monitor body composition [17], tail biting events [18] heart rate, respiration rate and body temperature [19].
Previous research showed that deviations in longitudinal data are heritable and these studies suggested that these resilience traits can be used as a proxy for pigs’ general resilience. Deviations in daily feed intake were estimated to be moderately heritable, with heritability (h2) estimates ranging from 8 to 26% [20], 7 to 11% [21], 9 to 23% [22]. Deviations in feeding behaviour, such as number of daily visits or daily visit duration, generally had higher heritabilities ranging from 16 to 20% [21], 16 to 28% [22] and 36 to 40% [23]. Heritability estimates of pigs’ body weight deviations range from 3 to 4% [24], 3 to 20% [22] and 31% [25], depending on the study. However, only a few studies looked at the relationship of these “indirect” traits with specific resilience-related traits in pigs, such as disease, mortality and social stressors. Putz et al. [20] and Cheng et al. [26] estimated moderate to high genetic correlations (rg = 0.37 to 0.85) between deviations in daily feed intake and mortality and antibiotic treatments, indicating that more deviations were associated with an increased mortality and an increased number of antibiotic treatments. In another study, a phenotypical association was found between deviations in activity levels and resilience in pigs infected with porcine reproductive and respiratory syndrome virus (PRRSV) [27]. Here, the odds of being non-resilient increased by 1.42 with a one-unit difference in root mean squared error (RMSE) of activity three days post-challenge.
Relationships between traits related to resilience and welfare, on the one hand, and deviations in longitudinal data, on the other hand, have already been observed in other species. In chicken, body weight deviations were positively related with lesion scores after infection [28]. In dairy cattle, a decrease in milk yield variations has been shown to be associated with improved longevity (rg = − 0.28 to − 0.34), udder health (rg = − 0.22 to − 0.32) and ketosis (rg = − 0.27 to − 0.33) [29]. In humans, general or systemic resilience based on deviations from natural time series, such as recovery time of blood pressure or heart rate, are related with overall mortality risk [4].
In the current study, resilience traits were estimated as deviations in body weight, deviations in weighing order and deviations in activity levels during weighing. First, resilience traits from deviations in longitudinal data were associated with physical abnormality scores and mortality. Second, variance components and heritabilities were estimated for these traits. And third, genetic correlations were estimated between resilience traits and traits related to resilience and welfare (e.g. mortality and tail biting) to evaluate if deviation traits are informative for breeders to improve pigs’ general resilience.
Methods
Animals and data collection
In total before quality control, 3696 finishing pigs with 20,833 body weight recordings from Vlaamse Piétrain Fokkerij (Belgium) were examined in this study. All pigs were crossbreds originating from Piétrain sires (N = 141) and Pig Improvement Company (PIC)-based dams bred in a criss-cross system (N = 352) with known pedigree. Piétrain sires were not preselected for this study and were a representative sample for the population. Most Piétrain sires from this study were homozygous positive for the halothane sensitivity (RYR1) genotype (93.3%), while 4.2% of the boars were heterozygous and 2.5% of the boars were homozygous negative.
Before weaning, piglets were vaccinated at the sow farm against Mycoplasma hyopneumoniae, as routinely done in Belgian pig production. At a mean age of 73.3 days (sd = 3.0) and weight of 23.2 kg (sd = 4.5), pigs were transferred from the sow farm to the experimental farm and upon arrival, pigs received anthelminthic medication and vaccination against Lawsonia intracellularis. Pigs were raised between July 2020 and July 2021 in a single pig building with 17 identical compartments composed of eight fully slatted pens (2.5 m × 4.0 m) per compartment. New pigs entered the stable in batches of 200 to 320 pigs every 2 weeks, filling two or three compartments with the progeny of 8 to 12 Piétrain sires. Per sire, a median of 26 crossbred piglets (half-sibs and full-sibs from the same sire) were equally distributed over two mixed-sex pens (sows and barrows) in the same compartment, mostly 13 piglets per pen. Once the stable was filled, it took about 2 months before a new round of pigs could enter, as the previous pigs were not sent to slaughter yet. During this experiment, we were able to evaluate 17 batches of pigs, distributed over 317 pens. Food and water were available to the animals ad libitum and were supplied via one trough and one nipple drinker per pen throughout the duration of the experiment. Over the course of the finishing period, which lasted about 120 days, observations were made until the pigs reached slaughter weight. The mean slaughter age and weight were 191.0 days (sd = 25.1) and 121.8 kg (sd = 10.4). The specific challenges presented to the pigs in our experiment were not recorded. However, pigs were kept in commercial conditions, automatically creating at least some challenges of various origin over the finishing period. Examples of such challenges are transportation stress upon arrival, disease pressure, heat stress in summer, social stress and stress related to a feed transition halfway the finishing phase.
Data collection
The data collection and experimental setup used in this study was previously described in detail by Gorssen et al. [30]. In short, pigs were weighed individually using radio frequency identification (RFID) ear tags and a RFID reader connected to a ground scale weighing platform. In total, 3358 pigs were weighed three days and 13 days after arrival on the farm. Hereafter, due to practical limitations, a subset of one out of two pens per sire was selected for subsequent weight recordings, at least at bi-weekly intervals (1919 pigs with 17,066 body weight recordings). For all pigs, slaughter weight and lean meat percentage were individually collected at the slaughterhouse of the Belgian Pork Group (site Meer, Belgium) using AutoFom III™ (Frontmatec, Smoerum A/S, Denmark) [31]. Moreover, cumulative feed intake was recorded weekly at the pen level for all pens. As explained by Gorssen et al. [32], weighing order on the pen level was recorded per weighing. Within a pen, the pig that was weighed first was scored with ‘1’, the second with ‘2’, etc.
In addition, during each weighing, pigs were visually scored individually on a set of physical abnormalities by the first author: presence of umbilical hernia (0 = not present, 1 = present), ear swellings or hematomas (0 = none, 1 = one ear, 2 = both ears), lameness (0 = not lame, 1 = lame), ear biting wounds (0 = none, 1 = one ear, 2 = both ears) and tail biting wounds (0 = none, 1 = small scratches, 3 = bloody and/or infected tail; (see Additional file 1 Figure S1). These scores were treated as continuous variables in further analyses. For tail biting wounds, we choose to not score bloody/infected tails with a score of ‘2’ but with a score of ‘3’, to give more weight to severely wounded tails. We would like to note that the used scoring method for these physical abnormalities is imperfect, as we assume that these scores are equidistant and can be treated as continuous variables. Although these measurements are imperfect and partly subjective, they will at least partly be related to pigs’ health, welfare and/or resilience status. Furthermore, although medical treatments were often started upon onset of problems, these treatments were not recorded in this study.
Estimates of activity levels during weighing using video analysis were only available for a subset of the total dataset and were recorded between January 2021 and July 2022 (1556 pigs, 7428 records), as described in Gorssen et al. [30]. In short, the activity level was estimated based on a linear combination of the mean speed, straightness index and sinuosity index of a pig’s tail base coordinates over time during weighing, with a 1/3 weight per trait after rescaling [31].
Ear tissue was collected for at least one finishing pig per sire-dam combination using the ‘Tissue Sampling Unit for DNA testing’ from Allflex® (New Zealand). In total, 413 finishing pigs were genotyped on medium-density single nucleotide polymorphism (SNP) chips: 122 finishing pigs were genotyped on the IMAGE porcine array (10,107 SNPs), whereas 291 pigs were genotyped on the GGP Porcine 50K chip from Neogen® (USA). Two SNP genotyping arrays were used, because these pigs were genotyped in different projects. Furthermore, the Piétrain sires were previously genotyped on the GGP Porcine 50K chip and the data were available for this research. Collection of ear tissues were approved by the Animal Ethics Committee of KU Leuven (P004/2020).
Data handling and quality control of weight records and genotypes
Before quality control (QC), 20,833 individual body weight records were available on 3696 pigs from 317 pens. First, we checked for duplicated records (same ID appearing in different pens or multiple recordings on the same day) and we removed duplicates if needed. Next, the number of unique sires per pen was checked and pens with more than one unique sire were removed (N = 1 pen, 24 records). Outlying weights were detected for pigs with at least four weight records using second order polynomial regression and comparing predicted with observed weights. Based on the histogram of deviations from predicted vs observed weights, deviations of more than three standard deviations were visually inspected. However, all weight trajectories seemed valid and no outlying weights were set to missing.
Genotype QC was done with the PLINK v1.9 software [33] (PLINK commands in brackets). The same QC settings were used for both arrays. Individual QC was first done on outlying missing (–missing) and heterozygosity (–het) values, with a minimum individual call-rate of 85% and heterozygosity maximum three standard deviations from the mean. Eight pigs did not pass this QC step due to 15 to 24% missing SNPs. Hereafter, QC of SNPs was performed by retaining only autosomal SNPs (–autosome) and setting the SNP call rate at a minimum of 95% (–geno 0.05). Lastly, duplications were checked (–genome), but all samples had relatedness (PI_HAT) < 0.90.
Finishing pigs genotyped on the IMAGE array [34] had 9072 SNPs and 122 pigs were kept after QC whereas for the finishing pigs genotyped on the GGP porcine array 45,824 SNPs and 283 pigs remained. After QC, arrays were merged (–bmerge) resulting in a dataset of 3769 overlapping SNPs and 405 finishing pigs. This dataset was subsequently merged with genotypes from Piétrain sire line (100% overlap, 2805 genotypes after merging) for parentage verification using the seekparentf90 software from the blupf90 suite of programs [35]. Parentage was not confirmed for three pigs. Although parentage was confirmed for their pedigree (half)sibs, it was decided to remove the complete pens with uncertain paternity (N = 3 pens, 313 records).
After QC for both phenotypic records and genotypes, the dataset contained 20,496 records on 3649 pigs (Fig. 1a). Resilience traits were constructed based on weight deviations, which were only estimated for pigs with at least five records. Consequently, a dataset of 17,066 records on 1919 pigs were left for analysing longitudinal resilience traits. These pigs, housed in 151 pens on paternal full- and half-sib basis, descended from 133 sires and 266 dams.
Trait definition
After QC, new traits were defined at the individual pig level. First, average daily gain (ADG) was determined as weight gain (in kg) divided by period of time in days. ADG from birth until first weight recording upon entering the finishing farm was named ADGyouth, whereas ADG from birth until last weight recording was named ADGlife. As birth weights were not recorded, weight at first/last recording was simply divided by age in days to calculate ADGyouth and ADGlife. For physical abnormalities, scored at each weighing event, the maximum score per pig was retained. For example, if a pig had severe tail biting wounds (score ‘3’), but this wound healed over time (score ‘0’ or ‘1’), the maximum score of ‘3’ was retained for statistical analysis. The sum of scores per individual was also calculated to estimate the duration of an event. However, the (genetic) correlations with a maximum score were very high. Therefore, only the results for maximum score are shown in this paper (detailed results not shown). Average feed intake (FI; kg/d) and feed conversion ratio (FCR; kg/kg) were estimated at the pen level. FI was estimated as the total amount of feed intake of the pen in kg, divided by the number of pigs and duration in days. FCR was estimated by dividing total feed intake (kg) by cumulative weight of a pen (kg), including pigs that died. Moreover, the coefficients of variation (CV) of pigs’ weight within a pen were estimated for weighing of pigs between 70 and 89 days of age (CVstart; mean age of 81.8 days, sd = 3.7), weighings between 115 and 130 days of age (CVmiddle; mean age of 120.8 days, sd = 3.0), and of weighings between 170 and 185 days of age (CVend; mean age of 176.2 days, sd = 2.8), as \(\text{CV}=\frac{\text{Standard deviation of weight at pen level}}{\text{Mean weight at pen level}}\).
Resilience traits were constructed based on deviations in longitudinal weight data. Theoretical or predicted weights were estimated using linear modelling and Gompertz growth curve modelling [36] in R [37]. As done by [9], a linear regression of weight as a function of age was performed (lmer function in R) and then the root mean squared error (RMSE) was estimated of observed versus predicted weights. This approach seems justified, as finishing pigs are more or less in their linear growth phase [38]. Then, we transformed RMSE using the natural logarithm (ln function in R) of the MSE (lnMSEweight). Gompertz modelling was performed with the nls function and the formula:
where the weight of an individual \(\text{i}\) at a given age in days \(\text{j}\) (weightij) can be predicted based on a combination of three Gompertz growth curve parameters (\({\text{A}}_{\text{i}}\), \({\text{B}}_{\text{i}}\) and \({\text{k}}_{\text{i}}\)) with a given time or age \({\text{t}}_{\text{ij}}\), whereas \({\upvarepsilon }_{\text{ij}}\) is residual error. Per individual and per weighing, this approach led to an observed weight and a predicted Gompertz weight (Fig. 2a and d). Hereafter, we calculated the difference between observed and predicted weights (Fig. 2b and e). Finally, the natural logarithm of the variance of these differences in observed versus expected weights was estimated (lnvarweight) according to [2]. Next, we also calculated the natural logarithm of the variance of a pigs’ weight after standardizing (lnvarstandardized), as explained in [28]. All weights were standardized by age with a bin of 3 days (for example, ages 75–77, 78–80,…), to ensure a sufficiently large population per age group. After standardization, weights were distributed with a mean of zero and a standard deviation of 1 per age bin (Fig. 1b). Weight standardization might be needed to correct for a scaling effect: as the mean weight increases over time in finishing pigs, the variance tends to increase as well [39]. Figure 2 shows the methodology for estimation of lnvarweight and lnvarstandardized for a pig with little deviations in body weight (Fig. 2a–c) and a pig with strong deviations in body weight (Fig. 2d–f).
Finally, deviations in pigs’ weighing order (lnMSEorder) and deviations in pigs’ mean activity levels during weighing (lnMSEactivity) were estimated similar to lnMSEweight. These traits were linearly regressed to age, and then the natural logarithm of MSE of observed versus predicted values was estimated. However, we were only able to record lnMSEorder and lnMSEactivity for a subset of 898 pigs because video recordings started halfway through our experiment, as explained in Gorssen et al. [30].
Statistical analysis of resilience traits and resilience-related traits
The statistical associations between the proposed resilience traits lnvarweight, lnvarstandardized, lnMSEweight and lnMSEactivity and resilience-related traits, such as mortality and tail and ear biting wounds were investigated. In first instance, boxplots and pairwise correlation plots were visually inspected for resilience traits according to the scores for resilience-related traits (Fig. 3) (see Additional file 2 Figure S2, Additional file 3 Figure S3, Additional file 4 Figure S4, Additional file 5 Figure S5, Additional file 6 Figure S6 and Additional file 7 Figure S7). Hereafter, linear mixed model regression analysis was applied (lmer package in R) as follows:
where, \(\mathbf{y}\) is the vector of phenotypes for a specific resilience trait (lnvarweight, lnvarstandardized, lnMSEweight or lnMSEactivity); \(\mathbf{b}\) is a vector of the fixed effects (sex, 2 levels) and covariables (specific resilience related trait score (tail wounds, ear wounds, hematomas, umbilical hernia, lameness, mortality); age at last recording and ADGlife); \(\mathbf{c}\) is the vector of common environmental pen effects (151 levels), following a normal distribution \(\mathbf{c}\sim N(\boldsymbol{0},{\mathbf{I}\upsigma }_{\text{c}}^{2})\), where \(\mathbf{I}\) is the identity matrix; \(\mathbf{e}\) is the vector of residual effects assumed to follow a normal distribution \(\mathbf{e}\sim N(\boldsymbol{0},{\mathbf{I}\upsigma }_{\text{e}}^{2})\); \(\mathbf{X}\) and \(\mathbf{W}\) are incidence matrices for, respectively, fixed effects, and the random pen effects (including both animal effects and common environmental effects).
Genetic modelling
Genetic parameters were estimated with the remlf90 software [35] using single-step genomic prediction (OPTION SNP_file). Heritability was calculated as the ratio of additive genetic variance to total variance, and common environmental pen effect (c2) as the ratio of variance explained by common environmental pen effects (c) to total variance. As described by [23, 40, 41] for an exponential model, the genetic coefficient of variation (GCV) was estimated for lnvarweight, lnvarstandardized, lnMSEweight and lnMSEactivity as: \(\text{GCV}=\sqrt{{\upsigma }_{\text{a}}^{2}}\). Next, we estimated genetic correlations (rg) via bivariate animal models for all trait combinations. The estimated animal models were as follows:
where, \(\mathbf{y}\) is the vector with phenotypes; \(\mathbf{b}\) is a vector with the fixed effects (sex, 2 levels; parity, 5 levels) and covariates (maximum age, mean time between consecutive records); \(\mathbf{a}\) is a vector of the additive genetic effects (5037 animals in pedigree, 2709 with genotype information). The assumption is that \(\mathbf{a}\) follows a normal distribution for the \(\mathbf{H}\) matrix following [42,43,44], using single-step genomic evaluation with both pedigree (\(\mathbf{A}\)) and genomic (\(\mathbf{G}\)) relationship matrices: \(\mathbf{a}\sim \text{N}(\boldsymbol{0},{\mathbf{H}\upsigma }_{\text{a}}^{2})\). \(\mathbf{c}\) is the vector of common environmental pen effects (151 levels), following a normal distribution \(\mathbf{c}\sim N(\boldsymbol{0},{\mathbf{I}\upsigma }_{\text{c}}^{2})\), where \(\mathbf{I}\) is the identity matrix; \(\mathbf{e}\) is the vector of residual effects assumed to follow a normal distribution \(\mathbf{e}\sim N(\boldsymbol{0},{\mathbf{I}\upsigma }_{\text{e}}^{2})\); \(\mathbf{X}\), \(\mathbf{Z}\) and \(\mathbf{W}\) are incidence matrices for, respectively, fixed effects, random animal effects and random common environmental pen effects. Within our experimental design, family groups are confounded with pens. The dataset contained 151 pens in total with an average of 12.7 pigs per pen (sd = 1.0; minimum 9.0; maximum 15.0) and every pen was sired by the same sire. Of these 151 pens, 7 pens consisted of full-sibs only, 118 pens consisted of a mixture of full-sibs and half-sibs originating from 2 dams per pen and 26 pens consisted of a mixture of full-sibs and half-sibs originating from 3 dams per pen. The average genetic relatedness across individuals within a single pen was 42.9% (sd = 10.8%; estimated from H-matrix). Within the complete set of 1919 studied pigs, the average genetic relatedness was 7.7% with a standard deviation of 5.3%.
Parity of the sow was divided in five classes: parity ‘1’ (N = 500), parity ‘2 to 3’ (N = 668), parity ‘4 to 5’ (N = 426), parity ‘6 to 7’ (N = 190) and parity ‘8+’ (N = 96).
Similarly, we estimated genetic correlations (rg) between traits using bivariate animal models as follows:
Similar to the single-trait animal model, \(\mathbf{y}\boldsymbol{1}\) and \(\mathbf{y}\boldsymbol{2}\) represent vectors with phenotypes for the studied traits; \(\mathbf{b}\boldsymbol{1}\) and \(\mathbf{b}\boldsymbol{2}\) are vectors of the fixed effects and covariates; \(\mathbf{a}\boldsymbol{1}\) and \(\mathbf{a}\boldsymbol{2}\) are vectors of additive genetic effects, assumed to follow a normal distribution for the \(\mathbf{H}\) matrix using single-step genomic evaluation:
\(\mathbf{c}\boldsymbol{1}\) and \(\mathbf{c}\boldsymbol{1}\) are vectors of common environmental pen effects, assumed to follow a normal distribution \(\left[\begin{array}{c}\mathbf{c}\boldsymbol{1}\\ \mathbf{c}\boldsymbol{2}\end{array}\right]\sim \text{N}(\left[\begin{array}{c}\boldsymbol{0}\\ \boldsymbol{0}\end{array}\right],\left[\begin{array}{cc}{\upsigma }_{\text{c}1}^{2}& {\upsigma }_{\text{c}1,\text{c}2}\\ {\upsigma }_{\text{c}1,\text{c}2}& {\upsigma }_{\text{c}2}^{2}\end{array}\right]\otimes \mathbf{I})\); \(\mathbf{e}\boldsymbol{1}\) and \(\mathbf{e}\boldsymbol{2}\) are vectors of residual effects, assumed to follow a normal distribution \(\left[\begin{array}{c}\mathbf{e}\boldsymbol{1}\\ \mathbf{e}\boldsymbol{2}\end{array}\right]\sim \text{N}(\left[\begin{array}{c}\boldsymbol{0}\\ \boldsymbol{0}\end{array}\right],\left[\begin{array}{cc}{\upsigma }_{\text{e}1}^{2}& {\upsigma }_{\text{e}1,\text{e}2}\\ {\upsigma }_{\text{e}1,\text{e}2}& {\upsigma }_{\text{e}2}^{2}\end{array}\right])\); \(\mathbf{X}\boldsymbol{1}\), \(\mathbf{X}\boldsymbol{2}\), \(\mathbf{Z}\boldsymbol{1}\), \(\mathbf{Z}\boldsymbol{2}\), \(\mathbf{W}\boldsymbol{1}\) and \(\mathbf{W}\boldsymbol{2}\) are incidence matrices for fixed effects, random animal effects and random common environmental pen effects, respectively.
Results
Association analysis of resilience traits with physical abnormalities and mortality
The evolution of (standardized) weights according to age (in Fig. 1) shows that finishing pigs are more or less in their linear growth trajectory and most animals are within three standard deviations of the mean. The sizes of the effect of the association analysis of resilience traits according to mortality and physical abnormality scores are in Table 1. Based on these results, lnvarweight, lnvarstandardized and lnMSEweight were strongly associated with mortality, lameness and tail biting wound scores. On average, pigs challenged by at least one of these conditions had significantly (p < 0.001) more deviations in body weight over time. This was also visible in the boxplots of resilience traits according to mortality and physical abnormality scores, which are shown in Fig. 3 for lnvarweight. For the other resilience traits, these boxplots are provided in Additional file 2 Figure S2, Additional file 3 Fig. 3, Additional file 4 Figure S4 and Additional file 5 Figure S5. These results show that tail biting wounds, lameness and mortality were associated with an increase in the deviations of predicted versus observed weight. An example of the association between tail biting wounds and lnvarweight is shown in Fig. 4. In this specific pen, tail biting wounds started from an age of 140 days on all pigs, indicating an outbreak of tail biting. For two pigs, this led to a drastic decrease in body weight (Fig. 4a), while most pigs recovered after an initial decrease in body weight gain. At the pen level (Fig. 4b), standardized weights decreased after the tail biting event, meaning that the pigs’ weight relative to the population average decreased.
Moreover, ear biting wounds were significantly associated with lnvarstandardized (p < 0.05), but not with lnvarweight (p > 0.10) and lnMSEweight (0.05 < p < 0.10). In contrast, deviations in weighing order (lnMSEorder) and activity during weighing (lnMSEactivity) were not associated with mortality and physical abnormalities (p > 0.10), except for a negative association between umbilical hernia and lnMSEactivity (p < 0.05).
The traits CVstart, CVmiddle, CVend, FI and FCR were recorded at the pen level. Therefore, phenotypical correlations of these traits with individually recorded traits were estimated using the pen averages. The pairwise correlation matrix of the mean of the resilience traits at the pen level compared with traits recorded at the pen level is shown in Additional file 8 Figure S8. Most correlations were significantly different from zero and significance values of these correlations are also in Additional file 8 Figure S8. Interestingly, body weight deviations (lnvarweight, lnvarstandardized, lnMSEweight) were phenotypically lowly correlated with CV at the pen level, especially at the end of the finishing phase (CVend) (rp = 0.14 to 0.29). Moreover, lnvarweight, lnvarstandardized and lnMSEweight had a negative correlation with FI (rp = − 0.23 to − 0.14) and a positive correlation with FCR (rp = 0.05 to 0.21), indicating that fewer deviations in body weight were lowly correlated with higher FI and lower FCR at the pen level. Correlations between body weight deviation traits and ADGlife were not consistent: lnvarweight had a positive correlation (rp = 0.10) with ADGlife, whereas this correlation was negative after standardizing weights (rp = − 0.10 between lnvarstandardized and ADGlife) [Table 3 and Additional file 9 Figure S9].
Genetic parameters
An overview of the estimated heritability and variance components is in Table 2. Gompertz curve parameters B and k had a moderate to high heritability (h2 = 40.3 and 47.3%). In contrast, the A parameter was estimated to be lowly heritable, but with a very high common environmental pen effect estimate (c2 = 91.2). Resilience traits from body weight deviations (lnvarweight, lnvarstandardized and lnMSEweight) were moderately heritable (h2 = 25.2 to 36.3%), but with very high GCV estimates of 44.5 to 59.5% (estimated from σa). Moreover, when standardized weights were used the common environmental pen effect estimate was lower (c2 = 7.7% for lnvarstandardized) compared to deviations in observed weights with Gompertz modeling (c2 = 8.4% for lnvarweight) and linear modeling (c2 = 14.3% for lnMSEweight). This indicates that scaling effects in observed weight over time might bias the estimates of weight deviations. Deviations in pigs’ weighing order (h2 = 4.2%) and activity during weighing over time (h2 = 12.0%) were lowly heritable and did not significantly differ from zero (using estimate ± 1.96*se). The heritability estimates for physical abnormality scores varied widely. Mortality (h2 = 6.1%), lameness (h2 = 8.3%) and umbilical hernia (h2 = 10.1%) had low heritability estimates. Tail wound score was moderately heritable (h2 = 16.4%), whereas hematomas (h2 = 41.6%) and ear wound scores (h2 = 46.3%) had high heritability estimates.
An overview of the estimated phenotypic and genetic correlations is in Table 3. Genetic correlations of lnMSEorder and lnMSEactivity were not taken into account, as standard errors were very large. The resilience traits based on body deviations (lnvarweight, lnMSEweight and lnvarstandardized) were internally lowly to moderately correlated (rp = 0.37 to 0.57; rg = 0.11 to 0.66, se = 0.19 to 0.26). Interestingly, body weight deviations were genetically positively correlated with tail biting wounds (rp = 0.09 to 0.21; rg = 0.22 to 0.39, se = 0.25 to 0.36), lameness (rp = 0.19 to 0.24; rg = 0.15 to 0.44, se = 0.05 to 0.70) and mortality (rp = 0.21 to 0.28; rg = 0.19 to 0.43, se = 0.07 to 0.09). It should be noted that the estimated standard errors of these genetic correlations were quite high, with an average standard error of 0.24 and a range from 0.00 to 1.47 (see Additional file 10 Table S1). ADGlife was positively correlated with lnvarweight (rp = 0.10, rg = 0.35, se = 0.18), but this positive correlation was not present after standardization of weights in lnvarstandardized (rp = − 0.10, rg = 0.02, se = 0.18). However, lnvarstandardized was positively correlated with meat percentage of the carcass (rp = 0.09, rg = 0.33, se = 0.18). Ear biting wound scores were positively correlated with tail biting wound scores (rp = 0.19; rg = 0.32, se = 0.33). Moreover, tail biting wound scores were also positively correlated with lameness (rp = 0.12; rg = 0.38, se = 1.07) and mortality (rp = 0.12; rg = 0.30, se = 0.12), and mortality and lameness were internally moderately to highly correlated (rp = 0.32; rg = 0.79, se = 0.13).
Discussion
This study investigated the (genetic) relationship between resilience traits quantified from deviations in longitudinal data in pigs with traits related to pig welfare, health and resilience, such as lameness, tail biting wounds and mortality. These resilience traits are promising proxies for a general resilience trait in pigs, but currently evidence is scarce for this putative relationship. Our results showed a highly significant (p < 0.001) association and a positive genetic correlation between resilience traits derived from body weight deviations (lnvarweight, lnvarstandardized and lnMSEweight) and tail biting wounds, lameness and mortality, although standard errors of genetic correlations were high. Moreover, these resilience traits were moderately heritable and phenotypically lowly favourably correlated with FI and FCR at the pen level. Furthermore, individual body weight deviations were correlated with uniformity at the pen level just before slaughter (CVend), providing evidence that breeding for these resilience traits might increase pigs’ within-individual uniformity and resilience as well as within-family uniformity at the pen level. In summary, our results suggest that breeding for our proposed resilience traits based on body weight deviations would increase pigs’ general resilience, uniformity and FCR at the pen level. These findings are highly valuable for pig breeders, and might be transferable to other (livestock) species.
Heritability estimates of lnvarweight, lnvarresilience and lnMSEweight were moderate ranging from 25.2 to 36.3%. These estimates are on the high end compared to previous estimates for resilience parameters based on deviations in body weight (h2 = 3 to 31%) [22, 24, 25], deviations in feed intake (h2 = 7 to 31%) [20,21,22,23] and deviations in feeding behaviour (h2 = 16 to 40%) [21,22,23]. In addition, GCV estimates for these traits were very high in this study, ranging from 44.5 to 59.5%. In contrast, Homma et al. and Gorssen et al. [22, 23], respectively, reported GCV estimates ranging from 22 to 39% and 21 to 33% for resilience parameters based on deviations in body weight, feed intake and feeding behaviour. We argue that these high heritability estimates and genetic coefficient of variation estimates might be due to a combination of high standardization and a reasonably challenging environment in our study, as explained in a previous study [30]. All studied pigs were kept in the same pig building during 1 year, hence, environmental variation might have been relatively low across pens, possibly reducing common environmental and/or residual variance. Moreover, our finishing pigs were kept under commercial conditions, with a higher pig density and disease pressure compared to purebred pigs in high health breeding farms as was done in Revilla et al. [24].
A significant (p < 0.01) and positive relationship between resilience traits based on body weight deviations (lnvarweight, lnvarresilience and lnMSEweight) and tail biting, lameness and mortality was found (Fig. 2 and Table 1) (Additional file 2 Figure S2, Additional file 3 Figure S3, Additional file 4 Figure S4, Additional file 5 Figure S5 and Additional file 6 Figure S6). Moreover, phenotypic and genetic correlations were also positive (rp = 0.09 to 0.28; rg = 0.15 to 0.44). These results indicate that breeding pigs for decreased deviations in body weight will also reduce tail biting, lameness and mortality. As such, these resilience traits might prove a valuable proxy for pig general health, welfare and resilience. Our positive correlation between deviations in body weight and mortality are consistent with Putz et al. [20] and Cheng et al. [26], who estimated moderate to high genetic correlations (rg = 0.37 to 0.85) between deviations in daily feed intake and mortality and antibiotic treatments. To our knowledge, our study is the first to report a favourable association between resilience parameters derived from deviations in body weight and tail biting wounds and lameness. However, several studies previously mentioned that average daily weight gain can be temporarily seriously affected in tail-bitten pigs [45,46,47,48] and in pigs suffering lameness [49].
We have to note that in our study, the CVstart, CVmiddle, CVend, FI and FCR traits were recorded at the pen level. Therefore, phenotypical correlations of these traits with individually recorded traits were estimated using the pen averages. Using this methodology, we found a positive correlation between resilience indicators from body weight deviations and coefficient of variation at the end of the finishing phase (rp = 0.14 to 0.29), showing that deviations on the individual level are also lowly to moderately correlated with uniformity at the pen level. These findings demonstrate that breeding pigs for an increased resilience based on individual body weight deviations (i.e., within-individual uniformity) could also improve the uniformity of pigs’ weight at the pen level (i.e., within-family uniformity). However, we were not able to estimate genetic correlations between within-individual uniformity and within-family uniformity due to the limited number of pens (N = 151). Future research could investigate this putative genetic relationship more in detail. In addition, lnvarweight, lnvarstandardized and lnMSEweight were also favourably correlated with FI (rp = − 0.23 to − 0.14) and FCR (rp = 0.05 to 0.21), implying that an increase in body weight deviations (lower resilience) is associated with a decreased FI and increased FCR at the pen level. Gorssen et al. [22] reported a similar favorable relationship between FCR and lnvarweight (rp = 0.17 and rg = 0.37) at the individual level, although they reported an unfavorable genetic correlation between FI and lnvarweight (rp = − 0.03 and rg = 0.26). As FI and FCR were recorded at the pen level, genetic correlations were not estimated due to the limited amount of pens (N = 151).
In contrast to deviations in body weight, lnMSEorder and lnMSEactivity had low heritability estimates (h2 = 4.2% and 12.0%), which did not significantly differ from zero. No clear association was found between lnMSEactivity and resilience-related traits, such as mortality. This is somewhat in contrast with the results of van der Zande et al. [27], who found a phenotypical association between deviations in activity levels and disease status and morbidity in pigs three days after infection with the porcine reproductive and respiratory syndrome virus (PRRSV).
Heritability estimates for tail wound scores (h2 = 16.4%) were low to moderate, and comparable to other studies reporting heritabilities for tail lesions, such as Martinsen et al. (h2 = 5%) [50], Hermesch et al. (h2 = 9 to 25%) [51] and a study on clinical tail biters by Breuer et al. (h2 = 0 to 27%) [52]. In our study, tail wound scores were also positively correlated with ear biting wound scores (rp = 0.19; rg = 0.32), which intuitively makes sense and was also found in previous studies [53,54,55,56]. Furthermore, ear biting wound scores and hematomas were estimated to be moderately to highly heritable (h2 = 46.3 and 41.6%). These results are interesting, as little is known about the heritability of these traits in pigs. Estimating heritabilities for ear and tail biting is challenging, as it is difficult to distinguish victims from biters and neutral pigs. In our study, however, every pen (~ 13 pigs) consisted of a mixture of half-sibs and full-sibs that originated from the same Piétrain sire and one to three crossbred dams. Therefore, pens with an incidence of tail and/or ear biting wounds had a similar genetic background, making it possible to estimate genetic variances for the occurence of biting behaviour, in general. However, this experimental design also complicated genetic modeling and can affect the accuracy of the estimates. We were quite rigid in our statistical modelling and used a random common environmental pen effect, although this pen-effect might partly overlap with the genetic component. Nevertheless, heritability estimates were significantly higher than zero. Litter effects were not included in the genetic model, as we believe this would overlap with genetic effects. A model including a random effect of the litter seemed to overestimate heritabilities for tail biting wounds (h2 = 77.3%) and ear biting wounds (h2 = 89.9%), whereas a model with both litter and pen as random effects experienced convergence issues (detailed results not shown). However, as we used an animal model and ssGBLUP while correcting for common pen effects, litter effects will be at least partially captured via the H-matrix. Within pens, the average genetic relatedness estimated from the H-matrix across individuals was 42.9%, with a mean standard deviation of 10.8%. Moreover, within the studied set of 1919 finishing pigs the average genetic relatedness was 7.6% with a standard deviation of 5.3%. Therefore, we argue that there was sufficient variability in relationships among pigs to accurately estimate genetic parameters.
No clear association was found with our resilience traits and ear biting wounds, although lnvarstandardized was significantly associated with ear biting wounds (p < 0.05) and had a positive genetic correlation (rp = 0.06; rg = 0.38). In our study, we observed that pigs usually already showed ear biting lesions upon arrival in the finishing farm. Only a few pigs showed new ear biting wounds after arrival, suggesting that this behaviour originated mostly in the farrowing and/or nursery farm. Although no data were available on the precise age of onset of ear biting in this study, previous studies [56] reported that ear biting behaviour often starts during the weaner stage (5–12 weeks of age) and mostly within a few weeks post-weaning. We observed that the prevalence of tail biting wounds increased at an age of approximately 140–160 days in our study. This might also explain why ear biting wounds do not show a relationship with all deviation in body weight traits in our dataset: as the ear biting challenge was mostly initiated before weight recording, the weight deviations might have occured before arrival at the finishing farm. The same reasoning might be applicable to ear hematomas: most hematomas originated early in life. This is in line with [57], who reported the highest incidence of hematomas 4 to 5 weeks after weaning. Therefore, we hypothesise that ear biting wounds and ear hematomas in pigs might be related with deviations in early-life body weight. For example, [57, 58] indicated that ear hematomas can temporarily negatively impact average daily gain up to 2 weeks after incidence, followed by compensatory growth. However, future research is necessary to test this hypothesis. This hypothesis also implies that our proposed resilience traits might be sensitive to the observation period, as indicated in Berghof et al. [2] and shown in Gorssen et al. [22]. Hence, deviations in body weight during the finishing phase might genetically differ from deviations in body weight during early-life in pigs. Furthermore, no significant association was found between our resilience parameters and umbilical hernias. Possibly, this is because umbilical hernias are a partly heritable congenital effect with a gradual development over a pigs’ lifetime. As a result, the presented challenge of an umbilical hernia might lower the average weight gain of pigs as observed in Searcy-Bernal et al. [59] and Straw et al. [60], but it might not cause substantial short-term deviations in observed versus expected weight. In addition, piglets with a umbilical hernia are generally excluded for progeny testing, and this preselection effect might also have influenced this relationship.
Although our proposed resilience traits lnvarweight and lnMSEweight are promising, it is important to note that these traits are based on deviations of observed versus expected weights using a model (Gompertz versus linear modeling). As indicated in a previous paper by Gorssen et al. [22], modeling theoretical weights based on observations is difficult, certainly for severely challenged animals (as can be seen in Fig. 2d, for example). Therefore, lnvarstandardized may be preferred as a resilience trait, as this trait does not rely upon a theoretical model, but reflects longitudinal changes of an individual animal’s weight compared to the population average. Moreover, Gorssen et al. [22] showed that lnvarstandardized was more robust to a low data density and standardizes the natural increase of body weight variance over time, as can be seen in Fig. 1a and b.
The proposed resilience traits have the advantage over recording physical abnormality traits only that they offer a way to refine phenotypes and allow a holistic approach [4]. Many resilience-related traits are recorded as binary (e.g. ‘lame’ vs not ‘lame’) or ordinal variables (e.g. tail wound score), although the underlying biological mechanisms are most probably complex and continuous [61]. As Fig. 3 shows, much of the variation in lnvarweight is still present in pigs without any visible physical abnormality. Pigs with a high lnvarweight score, but without visible deficiencies might be presented with other challenges, such as infections and/or disease, thermal discomfort, social stressors and/or bullying. However, we did not record these traits in the present study. Challenged pigs also showed great variation in lnvarweight scores. Some pigs with severe tail biting wounds, for example (Figs. 3 and 4), did not show much body weight deviations, indicating they might be more resilient to this stressor, whereas others were severely affected. Future studies should focus on this complex relationship.
For pig breeders, manually obtaining repeated records of pigs’ body weight records as was done in this study might be too costly and practically infeasible. However, new precision livestock technologies are developed with unprecedented speed, facilitating the collection of longitudinal data. For example, automated feeding stations are already used by most pig breeding companies. Furthermore, computer vision and wearable devices would allow to monitor almost continuously phenotypes such as body composition [17], tail biting events [18], heart rate, respiration rate and body temperature [19]. As a result, our findings will be valuable for pig breeders, as they show that resilience indicators based on deviations in longitudinal body weight in finishing pigs are heritable and associated with pig general health, welfare and resilience.
Conclusions
This study showed that resilience traits derived from deviations in body weight are promising proxies for a pig’s general health, welfare and resilience. Our results showed a highly significant (p < 0.001) association between body weight deviation traits lnvarweight, lnvarstandardized and lnMSEweight and tail biting wounds, lameness and mortality. Moreover, these resilience traits were moderately heritable and also genetically correlated with tail biting wounds, lameness and mortality. Furthermore, individual body weight deviations were positively correlated with uniformity at the pen level, providing evidence that breeding for these resilience traits might also increase pigs’ within-family uniformity. Our findings will be valuable for pig breeders, as they provide evidence that resilience indicators from deviations in body weight in finishing pigs are heritable and associated with pig general health, welfare and resilience.
Availability of data and materials
The datasets generated and/or analysed during the current study are not publicly available due to data restriction from the Vlaamse Piétrain Fokkerij vzw but are available from the corresponding author on reasonable request and with permission of the Vlaamse Piétrain Fokkerij vzw.
References
Putz A. Quantifying resilience in sows and wean-to-finish pigs. Ames: Iowa State University; 2019.
Berghof TVL, Poppe M, Mulder HA. Opportunities to improve resilience in animal breeding programs. Front Genet. 2019;9:692.
Colditz IG, Hine BC. Resilience in farm animals: biology, management, breeding and implications for animal welfare. Anim Prod Sci. 2016;56:1961–83.
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.
Knap PW. The scientific development that we need in the animal breeding industry. J Anim Breed Genet. 2020;137:343–4.
Bellini S, et al. The pig sector in the European Union. In: Lacolina L, Penrith ML, Bellini S, Chanais E, Jori F, Montoya M, et al., editors. Understanding and combatting African Swine Fever: A European perspective. Wageningen: Wageningen Academic Publishers; 2021. p. 183–94.
Kasper C, Ribeiro D, de Almeida AM, Larzul C, Liaubet L, Murani E. Omics application in animal science—a special emphasis on stress response and damaging behaviour in pigs. Genes. 2020;11:920.
de SouzaIung LH, Carvalheiro R, de Neves HHR, Mulder HA. Genetics and genomics of uniformity and resilience in livestock and aquaculture species: a review. J Anim Breed Genet. 2020;137:263–80.
Food and Agriculture Organization. Moving towards a more healthy and sustainable EU food system, a corner stone of the European Green Deal. Rome: Food and Agriculture Organization; 2020.
Rodenburg TB, van Gerwen M, Meijer E, Tobias TJ, Giersberg MF, Goerlich VC, et al. End the cage age: looking for alternatives: Overview of alternatives to cage systems and the impact on animal welfare and other aspects of sustainability. 2020. https://www.europarl.europa.eu/RegData/etudes/STUD/2020/658539/IPOL_STU(2020)658539_EN.pdf. Accessed 14 May 2024.
Mayorga EJ, Renaudeau D, Ramirez BC, Ross JW, Baumgard LH. Heat stress adaptations in pigs. Anim Front. 2019;9:54–61.
Rousi E, Kornhuber K, Beobide-Arsuaga G, Luo F, Coumou D. Accelerated western European heatwave trends linked to more-persistent double jets over Eurasia. Nat Commun. 2022;13:3851.
Van Zanten HHE, Van Ittersum MK, De Boer IJM. The role of farm animals in a circular food system. Glob Food Sec. 2019;21:18–22.
Gross WG, Siegel PB, Hall RW, Domermuth CH, DuBoise RT. Production and persistence of antibodies in chickens to sheep erythrocytes. 2. Resistance to infectious diseases. Poult Sci. 1980;59:205–10.
Gorssen W, Maes D, Meyermans R, Depuydt J, Janssens S, Buys N. High heritabilities for antibiotic usage show potential to breed for disease resistance in finishing pigs. Antibiotics. 2021;10:829.
Rydhmer L. Advances in understanding the genetics of pig behaviour. Cambridge: Burleigh Dodds Science Publishing; 2021.
Fernandes AFA, Dórea JRR, Valente BD, Fitzgerald R, Herring W, Rosa GJM. Comparison of data analytics strategies in computer vision systems to predict pig body composition traits from 3D images. J Anim Sci. 2020;98: skaa250.
Liu D, Oczak M, Maschat K, Baumgartner J, Pletzer B, He D, et al. A computer vision-based method for spatial-temporal action recognition of tail-biting behaviour in group-housed pigs. Biosyst Eng. 2020;195:27–41.
Wen H, Johnson JS, Henrique P, Freitas PHF, Maskal JM, Burd MK, et al. Genetic parameter estimation of various body temperature and respiration rate indicators in maternal line pigs under heat stress conditions. In: Proceedings of the 12th World Congress on Genetics Applied to Livestock Production: 3–8 August 2022; Rotterdam. 2022.
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.
Kavlak AT, Uimari P. Improving welfare of pigs through selection for resilience. In: Proceedings of the 73rd Annual Meeting of the European Federation of Animal Science: 5–8 September 2022; Porto. 2022.
Gorssen W, Winters C, Meyermans R, Chapard L, Hooyberghs K, Janssens S, et al. A promising resilience parameter for breeding: the use of weight and feed trajectories in growing pigs. J Anim Sci Biotechnol. 2023;14:101.
Homma C, Hirose K, Ito T, Kamikawa M, Toma S, Nikaido S, et al. Estimation of genetic parameter for feed efficiency and resilience traits in three pig breeds. Animal. 2021;15: 100384.
Revilla M, Lenoir G, Flatres-Grall L, Munoz-Tamayo R, Friggens NC. Quantifying growth perturbations over the fattening period in swine via mathematical modelling. Peer Community J. 2022;2: e9.
Laghouaouta H, Pena RN, Ros-Freixedes R, Reixach J, Díaz M, Estany J, et al. A methodology to quantify resilience in growing pigs. Animals. 2021;11:2970.
Cheng J, Putz AM, Harding JCS, Dyck MK, Fortin F, Plastow GS, et al. Genetic analysis of disease resilience in wean-to-finish pigs from a natural disease challenge model. J Anim Sci. 2020;98: skaa244.
van der Zande LE, Dunkelberger JR, Rodenburg TB, Bolhuis JE, Mathur PK, Cairns WJ, et al. Quantifying individual response to PRRSV using dynamic indicators of resilience based on activity. Front Vet Sci. 2020;7:325.
Berghof TVL, Bovenhuis H, Mulder HA. Body weight deviations as indicator for resilience in layer chickens. Front Genet. 2019;10:1216.
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.
Gorssen W, Winters C, Meyermans R, D’Hooge R, Janssens S, Buys N. Estimating genetics of body dimensions and activity levels in pigs using automated pose estimation. Sci Rep. 2022;12:15384.
Frontmatec. Fully automatic ultrasonic carcass grading AutoFom IIITM. Accessed 21 Sept 2023.
Gorssen W, Meyermans R, Janssens S, Buys N. Pigs’ weighing order stability is favourably linked with tail biting and production: a new proxy for social structure? In: Proceedings of the 12th World Congress on Genetics Applied to Livestock Production: 3–8 August 2022; Rotterdam. 2022.
Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience. 2015;4:7.
Tixier-Boichard M, Peynot N, Weigend S, Monteagudo L, Martinez Martinez A, Delgado JV, et al. Mapping genetic diversity in European gene banks: preliminary results on chickens for the validation of IMAGE001 array. In: Proceedings of the 12th World Congress on Genetics Applied to Livestock Production: 3–8 August 2022; Rotterdam. 2022.
Misztal I, Tsuruta S, Lourenco D, Aguilar I, Legarra A, Vitezica Z. Manual for BLUPF90 family of programs. Athens: University of Georgia; 2014.
Winsor CP. The Gompertz curve as a growth curve. Proc Natl Acad Sci USA. 1932;18:1–8.
R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2020.
Koivula M, Sevón-Aimonen ML, Strandén I, Matilainen K, Serenius T, Stalder KJ, et al. Genetic (co)variances and breeding value estimation of Gompertz growth curve parameters in Finnish Yorkshire boars, gilts and barrows. J Anim Breed Genet. 2008;125:168–75.
Lande R. On comparing coefficients of variation. Syst Zool. 1977;26:214–7.
Mulder HA, Bijma P, Hill WG. Prediction of breeding values and selection responses with genetic heterogeneity of environmental variance. Genetics. 2007;175:1895–910.
Hill WG, Mulder HA. Genetic analysis of environmental variation. Genet Res. 2010;92:381–95.
Legarra A, Aguilar I, Misztal I. A relationship matrix including full pedigree and genomic information. J Dairy Sci. 2009;92:4656–63.
Christensen OF, Lund MS. Genomic prediction when some animals are not genotyped. Genet Sel Evol. 2010;42:2.
Aguilar I, Misztal I, Johnson DL, Legarra A, Tsuruta S, Lawlor TJ. Hot topic: a unified approach to utilize phenotypic, full pedigree, and genomic information for genetic evaluation of Holstein final score. J Dairy Sci. 2010;93:743–52.
Leite NG, Knol EF, Nuphaus S, Vogelzang R, Tsuruta S, Wittmann M, et al. The genetic basis of swine inflammation and necrosis syndrome and its genetic association with post-weaning skin damage and production traits. J Anim Sci. 2023;101: skad067.
Wallgren P, Lindahl E. The influence of tail biting on performance of fattening pigs. Acta Vet Scand. 1996;37:453–60.
Sinisalo A, Niemi JK, Heinonen M, Valros A. Tail biting and production performance in fattening pigs. Livest Sci. 2012;143:220–5.
Schrøder-Petersen DL, Simonsen HB. Tail biting in pigs. Vet J. 2001;162:196–210.
Jensen TB, Kristensen HH, Toft N. Quantifying the impact of lameness on welfare and profitability of finisher pigs using expert opinions. Livest Sci. 2012;149:209–14.
Martinsen KH, Nordbø Ø, Enger EG, Olsen D, Grindflek E. Genetic analysis of feeding behaviour and tail lesions. In: Proceedings of the 12th World Congress on Genetics Applied to Livestock Production: 3–8 August 2022; Rotterdam. 2022.
Hermesch S, Guy SZY. Genetic analysis of tail-biting victims in pigs. In: Proceedings of the 23rd conference of the Association for the Advancement of Animal Breeding and Genetics: 27 Oct-1 Nov 2019; Armidale; 2019.
Breuer K, Sutcliffe MEM, Mercer JT, Rance KA, O’Connell NE, Sneddon IA, et al. Heritability of clinical tail-biting and its relation to performance traits. Livest Prod Sci. 2005;93:87–94.
Camerlink I, Ursinus WW, Bijma P, Kemp B, Bolhuis JE. Indirect genetic effects for growth rate in domestic pigs alter aggressive and manipulative biting behaviour. Behav Genet. 2015;45:117–26.
Ursinus WW, Van Reenen CG, Kemp B, Bolhuis JE. Tail biting behaviour and tail damage in pigs and the relationship with general behaviour: predicting the inevitable? Appl Anim Behav Sci. 2014;156:22–36.
Leite NL, Knol EF, Nuphaus S, Tsuruta S, Vogelzang R, Lourenco D. Swine inflammation and necrosis syndrome (SINS) and its association with biting behavior after weaning. In: Proceedings of the 12th World Congress on Genetics Applied to Livestock Production: 3–8 August 2022; Rotterdam. 2022.
Diana A, Carpentier L, Piette D, Boyle LA, Berckmans D, Norton T. An ethogram of biter and bitten pigs during an ear biting event: first step in the development of a precision livestock farming tool. Appl Anim Behav Sci. 2019;215:26–36.
Van Den Hof J, Beek J, Chiers K, Maes D. Risk factors for auricular hematomas in piglets. Vlaams Diergeneeskd Tijdschr. 2014;83:113–8.
Dewey C, Sunstrum J, Richardson K. Management of ear hematomas in pigs. J Swine Health Prod. 2018;26:137–41.
Searcy-Bernal R, Gardner IA, Hird DW. Effects of and factors associated with umbilical hernias in a swine herd. J Am Vet Med Assoc. 1994;204:1660–4.
Straw B, Bates R, May G. Anatomical abnormalities in a group of finishing pigs: prevalence and pig performance. J Swine Health Prod. 2009;17:28–31.
Guy S, Thomson PC, Hermesch S. Selection of pigs for improved coping with health and environmental challenges: Breeding for resistance or tolerance? Front Genet. 2012;3:281.
Acknowledgements
We would like to acknowledge Steven Broekx for his help in constructing the experimental setup. Furthermore, we would like to thank Nathalie Cenens, Simon Cardinaels and Hannelore Vermoesen for their help in weighing the pigs.
Funding
This study was partially funded by an FR PhD fellowship (1104320N; WG) and two SB PhD fellowships (1S05818N (CW) and 1S37119N (RM)) of the Research Foundation Flanders (FWO). Moreover, RM and LC were also funded by a KU Leuven C2 project, grant number C24/18/036. RM was also partially funded by Internal Funds KU Leuven (PDMT2/23/035). KH was funded by the UNIPIG project of VLAIO (HBC2019.2866). The funding bodies played no role in the design of the study,collection analysis, interpretation of data and in writing the manuscript.
Author information
Authors and Affiliations
Contributions
WG analysed the data and wrote the manuscript. WG, SJ. and NB designed and conceived this study. WG, CW, LC, KH, RM, JD, SJ, HM and NB critically reviewed the analyses and the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
All experimental procedures were approved by the Animal Ethics Committee of KU Leuven (P004/2020) in accordance with European Community Council Directive 86/609/EEC, the ARRIVE guidelines and the ILAR Guide to the Care and Use of Experimental Animals.
Consent for publication
Not applicable.
Competing interests
JD is an employee of Vlaamse Piétrain Fokkerij vzw and provided the data for this study, although Vlaamse Piétrain Fokkerij vzw did not fund this study. Prevalence of mortality and physical abnormality scores were not provided in this manuscript due to restrictions from Vlaamse Piétrain Fokkerij vzw. Apart from this restriction, all authors declare that the results are presented in full and as such present no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Figure S1.
Scoring sheet for physical abnormalities. Graphic representation of the scoring method of physical abnormalities in finishing pigs during the experiment. The scored abnormalities were tail biting wounds, ear biting wounds, ear hematomas or ear swellings and umbilical hernia.
Additional file 2: Figure S2.
Boxplots of lnvarstandardized according to mortality and physical abnormality scores. Boxplots of lnvarstandardized according to mortality and physical abnormality scores. Small dots indicate scores for individual pigs. Effect sizes and significance of these differences were statistically tested and are in Table 1.
Additional file 3: Figure S3.
Boxplots of lnMSEweight according to mortality and physical abnormality scores. Boxplots of lnMSEweight according to mortality and physical abnormality scores. Small dots indicate scores for individual pigs. Effect sizes and significance of these differences were statistically tested and are in Table 1.
Additional file 4: Figure S4.
Boxplots of lnMSEorder according to mortality and physical abnormality scores. Boxplots of lnMSEorder according to mortality and physical abnormality scores. Small dots indicate scores for individual pigs. Effect sizes and significance of these differences were statistically tested and are in Table 1.
Additional file 5: Figure S5.
Boxplots of lnMSEactivity according to mortality and physical abnormality scores. Boxplots of lnMSEactivity according to mortality and physical abnormality scores. Small dots indicate scores for individual pigs. Effect sizes and significance of these differences were statistically tested and are in Table 1.
Additional file 6: Figure S6.
Pairwise correlation plots of resilience traits according to physical abnormality scores at the individual level. Pairwise correlation plots for all evaluated resilience traits according to mortality and physical abnormality scores, based on individual scores. Below the diagonal the pairwise correlation plots are shown. Above the diagonal Pearson correlations are shown. °: correlation is significantly different from zero with p < 0.10. *: correlation is significantly different from zero with p < 0.05. **: correlation is significantly different from zero with p < 0.01. ***: correlation is significantly different from zero with p < 0.001.
Additional file 7: Figure S7.
Pairwise correlation plots of resilience traits according to physical abnormality scores at the pen level. Pairwise correlation plots for all evaluated resilience traits according to mortality and physical abnormality scores, based on the mean score at the pen level. Below the diagonal the pairwise correlation plots are shown. Above the diagonal Pearson correlations are shown. °: correlation is significantly different from zero with p < 0.10. *: correlation is significantly different from zero with p < 0.05. **: correlation is significantly different from zero with p < 0.01. ***: correlation is significantly different from zero with p < 0.001.
Additional file 8: Figure S8.
Pairwise correlation plots of resilience traits according to uniformity, feed intake and feed conversion ratio at the pen level. Pairwise correlation plots for all evaluated resilience traits according to coefficient of variation (CV) at the start, middle and end of the finishing period, as well as the feed conversion ratio. As these traits were recorded at the pen level, the pen-mean of the resilience trait was calculated. Below the diagonal the pairwise correlation plots are shown. Above the diagonal Pearson correlations are shown. °: correlation is significantly different from zero with p < 0.10. *: correlation is significantly different from zero with p < 0.05. **: correlation is significantly different from zero with p < 0.01. ***: correlation is significantly different from zero with p < 0.001.
Additional file 9: Figure S9.
Pairwise correlation plots of resilience traits according to production traits at the individual level. Pairwise correlation plots within all evaluated resilience traits and between resilience traits and production traits, such as the Gompertz growth curve parameters (A, B, k), average daily gain at start (ADGyouth) and over the pigs life (ADGlife), the recorded meat percentage at slaughterhouse, mean feed intake at pen level (FI) and mean feed conversion ratio at the pen level (FCR). Note that the correlations with FI and FCR might be inflated here, as these were recorded at the pen level, but correlated with traits recorded at the individual level. The correlations at the pen level, without possible inflation of the correlation were also constructed (see Additional file 8 Figure S8). Below the diagonal the pairwise correlation plots are shown. Above the diagonal Pearson correlations are shown. °: correlation is significantly different from zero with p < 0.10. *: correlation is significantly different from zero with p < 0.05. **: correlation is significantly different from zero with p < 0.01. ***: correlation is significantly different from zero with p < 0.001.
Additional file 10: Table S1.
Genetic correlations between all trait combinations using bivariate models. Genetic correlations between all trait combinations using bivariate models. Columns A (Trait1) and B (Trait2) give the trait names used in the bivariate animal model. Columns C (h2_trait1) and D (h2_trait2) show the estimated heritability for trait 1 and trait 2. Column E (genetic_correlation) shows the estimated genetic correlation, with the estimated standard error in column F (SE_genetic_correlation). Genetic correlations of lnMSEorder and lnMSEactivity were not taken into account, as standard errors were very high.
Rights and permissions
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 Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Gorssen, W., Winters, C., Meyermans, R. et al. Breeding for resilience in finishing pigs can decrease tail biting, lameness and mortality. Genet Sel Evol 56, 48 (2024). https://doi.org/10.1186/s12711-024-00919-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12711-024-00919-1