|
|
||||||||
ANIMAL GENETICS |

* Department of Animal Sciences, University of Illinois, Urbana-Champaign 61801; and
and
ARS, USDA, U.S. Meat Animal Research Center, Clay Center, NE 68933-0166
| Abstract |
|---|
|
|
|---|
Key Words: Heritability Population Mortality Sheep Survival Analysis
| Introduction |
|---|
|
|
|---|
Numerous studies have concluded that there is substantial genetic variation for the various causes of mortality. Mortality differences in sheep breeds associated with different causes have been reported (Nash et al., 1997
; Mukasa-Mugerwa et al., 2000
; Nguti et al., 2003
). Genetic variation in cause of mortality has been reported by Lush et al. (1948)
and Robertson and Lerner (1949)
in poultry, and by Gama et al. (1991)
in sheep. Dürr et al. (1999
and 2002
) and Roxström and Strandberg (2002)
identified genetic variability in dairy cow survival associated with different removal reasons. Therefore, the cause of mortality must be included when modeling mortality.
Competing risks analysis is an approach to account for cause of mortality when modeling mortality. A competing risk analysis considers that the competition among the different causes of mortality results in the observed mortality by one of the different mortality causes. Thus, a lamb that is alive has an associated risk of mortality due to each mortality cause. With the occurrence of mortality due to one cause, the risks of mortality due to the other causes are removed. The objective of this study was to assess the environmental and genetic factors associated with different causes of lamb mortality using complementary models in a competing risks approach.
| Materials and Methods |
|---|
|
|
|---|
Mortality records were obtained on lambs from F3 and advanced generations of a terminal sire composite population (breed composition: 50% Columbia, 25% Hampshire, and 25% Suffolk) at the U.S. Meat Animal Research Center, Clay Center, Nebraska. Leymaster (1991)
and Mousa et al. (1999)
have previously described the breed formation and management of this population. Lamb records were available from 1985 through 1997 and included the date of mortality and cause of mortality determined by necropsy. General descriptions of the different sheep mortality causes can be found in the Sheep Industry Handbook (American Sheep Industry Association, 2003
). Records from lambs that lost their identification were removed because the status of these lambs was unknown. Records from lambs raised in the nursery were removed because the nursery environment was expected to provide a different risk than the dam environment. The data consisted of 8,301 lambs from 299 sires and 2,475 dams.
Causes of mortality from birth to weaning were grouped into the following categories: disease (DIS; excluding pneumonia), dam-related (DAMR; including dystocia and starvation), pneumonia (PNEU), and other (OTHER; including eating disorders, injury, weather related, physical reasons, and unspecified causes). These categories were consistent with biological considerations of individual causes and had a minimum number of records necessary to obtain precise estimates. The overall (OVERALL) category included mortality by any cause and corresponded with the mortality definition used by Southey et al. (2001
, 2003)
. A summary of the number of records by cause of mortality is given in Table 1
.
|
Lamb mortality or failure time accounting for cause of mortality was analyzed using a competing risks approach by assigning a hypothetical potential or latent failure time to each category of mortality. Each category of mortality has an unobserved or latent time of failure that can occur provided that the other category of mortality cannot occur. The observed mortality time is the minimum of these latent failure times and results in the observed category of mortality. After the occurrence of mortality by one category, the nonminimum times associated with the other categories are no longer observable. Thus, the latent failure times associated with categories other than the observed category of mortality are considered to be right-censored at observed mortality because the mortality associated with these categories must be larger than the observed mortality (Crowder, 2001
). A hypothetical illustration of this concept considers that the latent failure times for a lamb from DAMR, DIS, PNEU, and OTHER categories are 3, 15, 40, and 50 d, respectively. The observed failure time or mortality record for this lamb would be 3 d, because this is the minimum of the latent failure times. This observed mortality record would be treated as observed for the DAMR category, but as censored at a value of 3 d for the DIS, PNEU, and OTHER categories.
Given the available information on time, cause of the event, and time-independent explanatory variables, it is not possible to distinguish between an independent competing risk model and any cause-specific hazard model (Crowder, 2001
; Kalbfleisch and Prentice, 2002
). To overcome the identifiability problem inherent to competing risks models, the causes of mortality were assumed to act independently; thus, the latent failure times are treated as independent (Crowder, 2001
; Kalbfleisch and Prentice, 2002
). This assumption was made for this analysis because there is no additional information on the potential differences between causes to justify any specific dependency.
The likelihood function of continuous-time mortality records can be factored into independent cause-specific likelihood functions associated with each competing risk (Allison, 1997
; Crowder, 2001
; Kalbfleisch and Prentice, 2002
). Within cause, the mortality records from other causes are treated as censored observations at the time of the mortality (Allison, 1997
; Crowder, 2001
; Kalbfleisch and Prentice, 2002
). Therefore, the competing risks model for continuous-time mortality can be implemented using standard survival analysis.
The likelihood function of discrete-time mortality records cannot be factored in a manner similar to the continuous-time mortality model, and requires the simultaneous estimation of all cause specific models (Allison, 1997
; Crowder, 2001
). Allison (1997)
described a procedure to analyze each cause separately based on the multinomial analysis approach of Begg and Gray (1984)
. This approach was implemented using a modification of the discrete-time approach described by Southey et al. (2003)
by creating separate animal-time data sets for each cause of mortality. In each data set, the mortality record of every lamb was initially divided into discrete 1-wk intervals. The number of discrete time intervals varied between causes to avoid a potential extreme value problem (Misztal et al., 1989
), with very few or no mortality recorded in an interval. The number of discrete-time intervals was 4 wk for the DAMR category, 6 wk for DIS category, and 7 wk for the PNEU and OTHER categories. Each interval consisted of a binary response variable that indicated either the occurrence or nonoccurrence for the specific cause of mortality. For example, a lamb that died in the third week due to PNEU would have three time intervals in the PNEU category data set, where the first two time intervals would be recorded as nonoccurrence of mortality and the third time interval would be recorded as occurrence of mortality due to PNEU. The same lamb in a data set for another category would consist of the first two time intervals being recorded as nonoccurrences of mortality and the third time interval recorded as censored.
A Weibull model was used for the continuous-time mortality records using the generalized linear mixed model class in Matvec (Kachman and Fernando, 2002
). A complementary log-log link function was used for the analysis of discrete-time mortality records in ASREML (Gilmour et al., 1999
) because the coefficients have the same relative risk interpretation as coefficients from the Cox proportional hazards model (Allison, 1997
; Kalbfleisch and Prentice, 2002
). The explanatory variables for both approaches were sex, age of dam (four levels: 1, 2, 3, and 4+ yr of age), type of birth (three levels: single-, twin-, and triplet or higher-born), and contemporary group (18 levels derived from year of lambing and breeding period within year). In addition, the discrete time interval was included as an explanatory variable in the discrete-time approach. Sire, animal, and maternal effects models were fitted for both approaches. The maternal effects model includes additive genetic variances of direct effects, additive genetic effects for maternal effects, and the possible covariance between the direct additive genetic and maternal effects. Heritabilities were computed as 
/(
+
2/6), where 
is the additive genetic variance and
2/6 is the variance of the Weibull model and complementary log-log function (Fahrmeir and Tutz, 1994
). The formula to compute the standard error of the heritability is provided in the ASREML manual (Gilmour et al., 1999
).
| Results |
|---|
|
|
|---|
|
Triplet- or higher-born lambs had significantly (P < 0.001) greater hazard of mortality than single- or twin-born lambs, and twin-born lambs had significantly (P< 0.001) greater hazard of mortality than single-born lambs in all categories. The highest hazard of mortality associated with type of birth occurred in the DAMR category, with single- and twin-born lambs having 81 and 71% lower hazard of mortality than triplet or higher-born lambs, respectively. The lowest hazard of mortality associated with type of birth occurred in the DIS and PNEU categories, with single-born lambs having 68 and 59% lower hazard of mortality than triplet-or higher-born lambs respectively, and twin-born lambs having 40 and 44% lower hazard of mortality than triplet or higher-born lambs, respectively.
Lambs from 3-yr-old dams had a risk of mortality similar to lambs from dams 4 yr of age and older in all categories. The OVERALL hazard of mortality of lambs from 2-yr-old dams was greater than of lambs from dams 4 yr and older but was not significantly different in any individual category. Lambs from 2-yr-old dams had significantly higher risk of mortality (P < 0.01) than lambs from 3-yr-old dams except in the DIS category, with the differences being marginal (P < 0.10). One-year-old dams provided a significantly (P < 0.001) poorer maternal environment because their lambs had a greater hazard of mortality than lambs from older dams.
The probability of lamb mortality varied depending on the discrete time interval and was greatest in the first week of life for all categories (Figure 1
). The probability of lamb mortality, after adjusting for the terms in the model, during the first week was 9.8% in the OVERALL category, indicating that for the average lamb, there was a 10% probability of mortality in the first week. The probability of mortality in the OVERALL category was under 2% from the second week of age until weaning.
|
The variance and associated heritability estimates for different categories are presented in Table 3
for the continuous-time approach and in Table 4
for the discrete-time approach. Similar estimates of sire variance and heritability were obtained for the OVERALL cause and individual categories for the continuous- and discrete-time approaches using the sire model. In general, the continuous-time approach using the animal model provided greater estimates of additive genetic variances and, hence, greater heritability estimates than the discrete-time approach. The maternal effects model only converged for the continuous-time analysis of OVERALL and PNEU causes, assuming a zero correlation between direct and maternal effects. For the OVERALL category, the estimated direct additive variance was 0.450 ± 0.119 and the estimated maternal variance was 0.206 ±0.073, resulting in an additive heritability estimate of 0.213 ±0.048. For the PNEU category, the estimated direct additive variance was 0.539 ±0.228 and the estimated maternal variance was 0.112 ±0.153, resulting in an additive heritability estimate of 0.235 ±0.080.
|
|
| Discussion |
|---|
|
|
|---|
The influence of the maternal environment in this study was in agreement with Gama et al. (1991)
, who analyzed different definitions of lamb mortality, including a respiratory cause, from birth to weaning. Gama et al. (1991)
reported low heritability estimates (<10%) of mortality using a half-sib model, but higher estimates of heritability for total mortality from birth to weaning (30%) and birth to 24 h (89%) under full-sib model. Gama et al. (1991)
reported similar estimates of heritability associated with respiratory causes between half-sib and full-sib approaches. This suggested that, in agreement with this study, there was minimal maternal influence on respiratory causes of mortality.
Gama et al. (1991)
used a cumulative binary approach, where the response variable was defined as occurrence or nonoccurrence. In contrast, this study used a discrete-time approach that Southey et al. (2003)
had shown to provide similar estimates of variance components to the continuous-time and cumulative binary analysis reported by Southey et al. (2001)
. This cumulative binary approach could lead to biased parameter estimates in survival analysis with a single cause of mortality due to different degrees of censoring within the levels of an explanatory variable. Ingram and Kleinman (1989)
showed that the cumulative binary analysis provided large biases when censoring only occurred in one level of a two level factor but showed no bias when censoring occurred in both levels of the factor. Consequently, the cumulative binary analysis applied to the data sets used in this study requires additional caution compared with the discrete-time analysis due to the increased censoring that occurs by treating mortality from other causes as censored.
The estimates of additive genetic variance suggest that identifying different categories of mortality would permit potential selection against specific causes of interest. Estimates of genetic parameters and selection experiments indicate that selection on mortality due to the factors that contribute to the DAMR category can be effective in decreasing mortality in this category. Cloete et al. (2002)
reported low heritability estimates for lambing and neonatal behavior characteristics, and Hinch (1997)
found substantial breed differences. Selection for improved rearing ability has also improved lamb survival (Atkins, 1980
; Donnelly, 1982
; Haughey, 1983
). Cloete and Scholtz (1998)
observed that lambs from the high rearing ability line had higher survival rates and were quicker to progress from standing to sucking than lambs from the low rearing ability line. Indirect response to selection has also increased pelvic size (Knight et al., 1988
; Kilgour and Haughey, 1993
), thus reducing lambing difficulty. Bennett (2001)
reported decreased calving difficulty in cattle selected for improved calving ease.
Cold resistance as an indicator of exposure is another important determinant of lamb survival in many environments. Although this aspect was not directly considered in this study, there have been reports of genetic variation. Experiments with cold resistance using water bath tests have reported heritability estimates of 0.3 (Slee and Stott, 1986
), 0.36 (Wolff et al., 1987
), and 0.7 (Slee et al., 1991
). Slee and Stott (1986)
reported a realized heritability of 0.27 for increased cold resistance but a realized heritability of zero for decreased cold resistance. In addition, a major gene associated with cold resistance has been identified (Slee and Simpson, 1991
; Simpson and Slee, 1998
).
The competing risks approach used here assumed that the categories of mortality are independent, even though some of the categories may not be strictly independent, such as PNEU and OTHER diseases. The grouping of individual causes into broader categories can result in loss of information when compared with the study of specific individual causes. However, the grouping is expected to increase the power if the grouped causes share the genetic mechanism of control. For example, the Nramp locus in mice is known to provide resistance to a variety of diverse and unrelated pathogens (Vidal et al., 1993
).
The magnitude of the correlations between breeding value estimates from different mortality categories was used as an indicator of the suitability of the independence of categories assumption. The correlations between the individual categories were close to zero (range = 0.12 to 0.19) and smaller than the correlation between each individual category and the OVERALL category (range 0.18 to 0.50). Thus, the response to selection for mortality, regardless of cause, is likely to be limited because few animals would be favorable for several causes. The low magnitudes of the correlations also suggest that the assumption of independent causes may be appropriate since high correlations would be expected if the causes were indicators of the same response variable. However, the use of breeding values to assess independence is challenging due to bias associated with censoring (Vukasinovic et al., 1999
), the lack of identifiability in competing risks, and because only the minimum mortality time of the potential failure times is recorded.
The estimates of the variance components for OVERALL mortality are in agreement with our previous reports (Southey et al., 2001
, 2003
) applying continuous-and discrete-time survival analysis. These similarities between the different approaches were expected from the results of Ingram and Kleinman (1989)
and Doksum and Gasko (1990)
and suggest that a variety of approaches using different software can be applied to study mortality. Our studies have shown that there is genetic variation for lamb survival assuming either a single cause or multiple causes. The estimates of individual categories support the suggestion of Cundiff et al. (1982)
that selection of components of mortality is likely to be more effective than selection solely on mortality, regardless of cause.
| Implications |
|---|
|
|
|---|
| Footnotes |
|---|
2 Correspondencee-mail: southey{at}uiuc.edu.
Received for publication December 17, 2003. Accepted for publication June 14, 2004.
| Literature Cited |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
W. Mekkawy, R. Roehe, R. M. Lewis, M. H. Davies, L. Bunger, G. Simm, and W. Haresign Genetic relationship between longevity and objectively or subjectively assessed performance traits in sheep using linear censored models J Anim Sci, November 1, 2009; 87(11): 3482 - 3489. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Maxa, A. R. Sharifi, J. Pedersen, M. Gauly, H. Simianer, and E. Norberg Genetic parameters and factors influencing survival to twenty-four hours after birth in Danish meat sheep breeds J Anim Sci, June 1, 2009; 87(6): 1888 - 1895. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. M. Dwyer Genetic and physiological determinants of maternal behavior and lamb survival: Implications for low-input sheep management J Anim Sci, April 1, 2008; 86(14_suppl): E246 - E258. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |