|
|
||||||||
ANIMAL GENETICS |
,3
* MTT Agrifood Research Finland, Animal Production Research, Animal Breeding, Jokioinen, Finland; and
and
Department of Animal Science, Iowa State University, Ames 50011-3150
| Abstract |
|---|
|
|
|---|
Key Words: Genetic Correlation Heritability Lifetime Prolificacy Longevity Sow
| Introduction |
|---|
|
|
|---|
Sow longevity can be defined as "stayability" or length of productive herd life (LPL). Stayability is a binary trait, measuring whether a sow has survived in a herd until some defined fixed parity or time. Similarly, LPL is the number of days between a beginning event, such as date of birth or date of first farrowing and culling of a sow. These traits continue to be breeding objectives for many seedstock suppliers; however, indirect selection for LPL and LTP has been traditionally employed in most breeding programs because productive life can be recorded only after a sow has been culled or death occurs. Indirect selection for this type of trait may be more useful because it can be carried out with a shorter generation interval.
To compare the efficiency of direct and indirect selection for sow longevity, the genetic parameters (heritabilities and genetic correlations) should be known. The effectiveness of indirect and direct selection for swine longevity traits is likely population-dependent and should be evaluated before selection is actually implemented. Information regarding the effectiveness of direct and indirect selection potential for sow longevity is sparse, particularly for the Finnish Landrace and Large White populations. Genetic correlations between sow longevity and other economically important swine traits are needed to determine the importance of their association.
| Materials and Methods |
|---|
|
|
|---|
|
Statistical Analyses
In the analysis of survival data, the main challenge is to use information on animals that are still alive, and to model factors that are time-dependent (e.g., the effect of farm or some disease). Ducrocq and Sölkner (1998)
developed the Survival Kit to analyze this type of data. Using proportional hazard models, it is possible to account for censored records, to model time-dependent factors, and to fit the distribution of longevity data (Ducrocq and Sölkner, 1998
); however, it is possible to run only single-trait analysis with the Survival Kit. Therefore, a single-trait proportional hazard model and a multitrait linear model were fitted to the current data. The single-trait analyses were carried out with the Survival Kit (Ducrocq and Sölkner, 2001
), and multitrait analyses with DMU package (Madsen and Jensen, 2000
).
In survival analysis, the hazard function of a sows LPL, t days after first farrowing, was defined as:
![]() |
where 
(
t)
1 is the Weibull baseline hazard function with location (
) and shape (
) parameters, ß is a vector of fixed and random effects, and x(t) is the corresponding incidence matrix. The suitability of the Weibull distribution was assessed by plotting ln{ln[S(t)]} against ln(t), where S(t) is the Kaplan-Maier survivor function (Figure 1
). The plot of this relationship produced a relatively straight line, which indicates that the Weibull distribution fits the data very well. Moreover, the lines of different stratums produced from the data used in this study were approximately parallel; hence, one baseline hazard function was assumed over the data.
|
where
is the sire variance, and A is the additive relationship matrix. The heritabilities of LPL on the logarithmic scale of survival analysis were calculated as showed by Ducrocq (2001)
![]() |
where
2/6 is the variance of extreme value distribution.
In multitrait analysis, the following linear sire model, in matrix notation, was fitted for all of the traits studied using DMU-package (Madsen and Jensen, 2000
):
![]() |
where y is the vector of observations of the seven traits considered simultaneously (LPL, LTP, NW, AFF, FFI, ADG, and either backfat thickness or leg score), X and Z are the incidence matrices for fixed (b) and random sire (u) effects, and e is the vector of residuals. The censored records (sows sold or still alive) of LPL and LTP were treated as missing. The effect of farm and year interaction was included in vector b for all the traits studied. In addition, farrowing month (NW, AFF, ADG, backfat thickness, leg score), mating type (NW), age of litter at weaning (NW, FFI), and breeding consultant (ADG, fat percent, leg score) were the other effects included in vector b. Moreover, the fixed regression of test weight was included in the statistical models for ADG and fat percent. The (co)variance structure of random effects of u, and e were assumed to be var(u) = A * G0, and var(e) = I * R0, where A is additive relationship matrix among sires, I is the identity matrix, and G0 and R0 are the variance-covariance matrices for additive genetic sire and residual effects, respectively.
Two types of "genetic" correlations between LPL and the other traits studied are presented. From multitrait analysis, the genetic correlations were obtained directly. In addition, the correlations between estimated breeding values for LPL in survival analysis and the breeding values of other traits studied in linear model analysis were estimated.
| Results |
|---|
|
|
|---|
|
|
|
Correlations Between Breeding Values
Correlations of the LPL breeding values obtained from survival analysis with the other traits studied are presented in Tables 3
and 4
. Although the correlations among LPL breeding values estimated using linear model and survival analysis were moderate or high (Pearson correlations ranged between 0.40 and 0.72, and Spearman rank correlations similarly ranged between 0.38 and 0.72), it should be noted that differences in sire ranking between these methods occurred. The correlations were higher in the Large White (Pearson correlations = 0.52 to 0.72; Spearman rank correlations = 0.48 to 0.72) compared with those from the Landrace (Pearson correlations = 0.40 to 0.57; Spearman rank correlations = 0.38 to 0.55) analysis. Similarly, survival analysis breeding values for true LPL (Spearman rank correlations was 0.55 in Landrace and 0.72 in Large White) had higher correlation with linear model breeding values than those for functional LPL (the corresponding correlations being 0.38 and 0.48).
In general, the correlations between the estimated breeding values were lower than the corresponding genetic correlations obtained from a multiple trait analysis (Tables 3
and 4
). However, most of the estimates had the same sign (12 out of 14 estimates), and the tendency for changes in the magnitude of correlation estimates were similar between the two methods.
Phenotypic and Genetic Correlations
Genetic and phenotypic correlations are presented in Tables 3
and 4
. In general, phenotypic correlations were very similar between the breeds, whereas some genetic correlations differed between breeds. Moreover, the phenotypic correlations were commonly very low. However, there was an indication that LPL and LTP are very closely associated, as both phenotypic and genetic correlations were greater than 0.90. Because of that, the correlations with the other traits studied are very similar between LPL and LTP.
All the prolificacy traits were genetically correlated with LPL and LTP and the correlations were generally greater than 0.13 (Tables 3
and 4
). This indicates that selection for more piglets weaned in the first litter and for short first farrowing interval will have a beneficial indirect effect for LTP and for LPL. The absolute values of these genetic correlations ranged between 0.30 and 0.54 and were similar between the Landrace and Large White populations; however, genetic correlation differences for AFF with LTP and LPL were found between the breeds. In Landrace, the correlations were positive (0.17 with PLP and 0.13 with LTP), whereas they were negative in Large White (0.28 with LPL and 0.29 with LTP).
In Landrace, the genetic correlations between overall leg score and longevity (0.32 with LPL and 0.28 with LTP; Table 3
) were moderate but positive. The corresponding correlations in Large White were 0.17 and 0.19 (Table 4
). Although the correlations of overall leg action with LPL and LTP were low to moderate, it may be said that selection for leg conformation, measured when the sow has reached 100 kg, is beneficial for improving sow longevity in an indirect manner.
Among the other traits studied, there were moderate genetic and phenotypic correlations between ADG and backfat thickness in both breeds (rg = 0.32 and rp = 0.40 in Landrace and rg = 0.39 and rp = 0.40 in Large White; Tables 3
and 4
). There were genetic correlations between AFF and FFI (average over breeds = 0.40), between AFF and ADG (0.40), and between backfat thickness and FFI (0.29). In addition, there were a few moderate genetic correlations present only in the Landrace breed (NW and daily gain = 0.27; AFF and backfat = 0.18; AFF and Score = 0.37; Tables 3
and 4
).
| Discussion |
|---|
|
|
|---|
Heritability Estimates for LPL
The heritability estimates for LPL were lower for the linear model (0.05 to 0.10) than for the survival analysis (0.16 to 0.19) in the current study. This is in agreement with the estimates presented in the literature (Tholen et al., 1996a
,b
; López-Serrano et al., 2000
), which reported stayability heritability estimates ranging from 0.02 to 0.11 when analyzed using linear models. The heritability estimates have ranged between 0.11 and 0.31 when analyzed using survival analyses (Yazdi et al., 2000a
,b
). Although heritability estimates from linear model analyses had been generally lower than those from survival analyses, moderate heritabilities estimates were obtained when censored records were accounted in the linear model analyses (Guo et al., 2001
). In their analyses, the estimated heritabilities for LPL and LTP ranged between 0.22 and 0.25, which were substantially higher than the current estimates when a linear model was applied and censored records were not considered (Guo et al., 2001
).
Current heritability estimates indicate that environmental effects may be modeled more precisely in survival analysis compared with the linear model analysis (i.e., it may be assumed that the higher heritability estimates are due to the possibility to model farm-year effect as time-dependent). On the other hand, it should be remembered that the residual variance was not estimated in the survival analysis, and the heritability estimates were calculated assuming that the residual effects were following extreme value distribution with the variance
2/6. Therefore, one might argue that the different heritability estimates are not comparable.
Correlation Estimates
The results of the current study demonstrate that LPL and LTP are very similar, both genetically and phenotypically (rg = 0.96 in Landrace and 0.97 in Large White; rp = 0.94 in both the breeds). Additionally, the correlations of LPL and LTP with the other traits studied were very similar. In general, our correlation estimates between LPL and prolificacy traits are very similar to those reported by Tholen et al. (1996b)
, except that their genetic correlations between litter size at first parity and stayability were negative in one dataset, ranging from 0.04 to 0.25. However, the same study reported correlation estimates between litter size and stayability ranging from 0.25 to 0.45, which were of the same magnitude as the current estimates between number of piglets weaned and LPL (0.39 in Landrace and 0.30 in Large White). Similarly, the genetic correlations reported by Tholen et al. (1996b)
for farrowing interval and stayability ranged between 0.24 and 0.54, compared with current estimates, which ranged between 0.40 and 0.43. Moreover, fixed effects of litter size and farrowing interval also significantly affected LPL in a Swedish pig population (Yazdi et al., 2000a
,b
).
Although the estimated genetic correlations of LPL and LTP with litter size and farrowing interval were very similar between the breeds, the correlations with AFF differed between the breeds (Tables 3
and 4
). In Large White, negative genetic associations between AFF and longevity (0.28 with LPL and 0.29 with LTP) were found, whereas the signs in corresponding correlations were positive in the Landrace (0.17 and 0.13) analyses. Similarly, the number of piglets weaned was positively correlated (0.21) with AFF in the Large White analysis, whereas the correlation in the Land-race population was approximately zero (Table 3
). These results concur with our previously reported studies, where the genetic correlations between AFF and litter size, and piglet survival have been different between the two breeds (Serenius et al., 2004a
,b
). Thus, the genetic correlation differences associated with AFF among the Finnish Landrace and Large White breeds is likely due to differences in the genetic makeup of the two breeds.
In the current analysis, there was no clear genetic association of LPL and LTP with production traits. The only substantial genetic correlation exists between LPL (and LTP) and backfat thickness in the Large White population (0.22). These findings are in agreement with the literature estimates, where both zero and unfavorable genetic correlations of LPL and LTP with backfat thickness and ADG have been reported. In the study of López-Serrano et al. (2000)
, stayability was positively correlated with backfat thickness, ranging from 0.11 to 0.27 and negatively associated with ADG (ranging from 0.06 to 0.32). However, Tholen et al. (1996b)
reported genetic correlation estimates between stayability and backfat thickness that ranged from 0.03 to 0.36, and between stayability and ADG from 0.02 to 0.13. Thus, genetic correlations among longevity and production traits seem to differ depending on the population from which they are estimated. Additionally, the differences could be the result of variation in the standardized weight at which backfat measures were obtained, back-fat measurement procedures, or other testing methodology.
In Landrace, a moderate genetic correlation (0.32) was found between LPL and leg conformation score. The genetic correlation between the same traits in the Large White breed was approximately half (0.17) that found in the Landrace breed. Similarly, López-Serrano et al. (2000)
reported that there was a positive genetic correlation between stayability and leg score in the German Landrace population, and the corresponding correlation was approaching zero in the German Large White population. It may be concluded that leg conformation is genetically correlated with LPL. However, it should be recognized that the reported heritabilities of leg conformation traits are very low (0.06 in the current study); hence, one should not expect large nor rapid genetic improvement in LPL through selection for leg conformation.
Based on heritability estimates, it seems that survival analysis may be the most appropriate method of evaluating swine longevity traits compared with linear models. However, there is one major concern relating to the use of survival analysis in breeding value estimation: because of computational problems, multiple-trait analyses involving longevity and other economically important traits are not currently possible. As stated earlier, LPL and LTP have genetic correlations with litter size, farrowing interval, and leg conformation that are relatively high compared with genetic correlations among longevity and measures of leg soundness.
| Implications |
|---|
|
|
|---|
| Footnotes |
|---|
3 Collaborator, via a fellowship under the OECD Cooperative Research Program: Biological Resource Management for Sustainable Agriculture Systems. ![]()
2 Correspondence: 109 Kildee Hall (phone: 515-294-4103; fax: 515-294-5698; e-mail: timo.serenius{at}mtt.fi).
Received for publication April 24, 2004. Accepted for publication July 16, 2004.
| Literature Cited |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
J. Tarres, J. P. Bidanel, A. Hofer, and V. Ducrocq Analysis of longevity and exterior traits on Large White sows in Switzerland J Anim Sci, November 1, 2006; 84(11): 2914 - 2924. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Serenius, K. J. Stalder, T. J. Baas, J. W. Mabry, R. N. Goodwin, R. K. Johnson, O. W. Robison, M. Tokach, and R. K. Miller National Pork Producers Council Maternal Line National Genetic Evaluation Program: A comparison of sow longevity and trait associations with sow longevity J Anim Sci, September 1, 2006; 84(9): 2590 - 2595. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Serenius and K. J. Stalder Selection for sow longevity J Anim Sci, April 1, 2006; 84(13_suppl): E166 - E. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Arango, I. Misztal, S. Tsuruta, M. Culbertson, and W. Herring Study of codes of disposal at different parities of Large White sows using a linear censored model J Anim Sci, September 1, 2005; 83(9): 2052 - 2057. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |