|
|
||||||||



* Institut National de la Recherche Agronomique, Station de Recherches Zootechniques, Domaine Duclos, 97170 Petit Bourg, Guadeloupe (French West Indies);
and
INRA, Station de Génétique Quantitative et Appliquée,78352 Jouy-en-Josas Cedex, France;
and
INRA, Domaine de Gardel, 97129 Moule, Guadeloupe(French West Indies);
and
INRA, Département de Santé Animale, 37380 Nouzilly, France
| Abstract |
|---|
|
|
|---|
Key Words: Genetic Variation Goats Mortality Strongylidae Survival
| Introduction |
|---|
|
|
|---|
This article reports a study of mortality in Creole goat, a highly prolific meat breed of the Caribbean characterized by Alexandre et al. (1999)
and Mandonnet et al. (2002)
. Since genetic variability of resistance to gastrointestinal nematode parasites was assessed (Mandonnet et al., 2001
), breeding for improved resistance in postweaning Creole kids is one way for reducing mortality. However, managing flock effects that increase risk of death could lead to more rapid improvement of tolerance in kids. Survival analysis methodology (Klein and Moeschberger, 1997
) was considered to be the method of choice to assess variability of mortality in our experimental conditions. Since our main objective was to estimate genetic variability of resistance criteria, kids were drenched to reduce mortality rate. Even in situations of low incidence, such as ours, survival models proved to be efficient tools (Ducrocq et al., 2000
). The rate at which the animals died is described without restricting the analysis to predefined cutpoints in time. The models use the available information, regardless of whether it comes from animals alive at the end of the period or from dead animals.
The aims of this study were to identify management factors influencing kids mortality due to infection with gastrointestinal strongyles and to estimate genetic parameters of this trait using survival analysis.
| Materials and Methods |
|---|
|
|
|---|
The data of 837 kids born between 1995 and 1998, sired from 48 bucks and 250 does, were collected from the experimental Creole goat flock of INRA-Gardel, located in the French West Indies. An intensive reproduction system was applied (three parturitions every 2 yr). The kidding periods were the dry season (mid-February through mid-March), the intermediate season (mid-June through mid-July), and the wet season (mid-October through mid-November). The flock produced, on average, a 200-kid cohort every 4 mo. Throughout the year, the animals grazed on Digitaria decumbens irrigated pastures managed in a rotational system. Kids were weaned at an average age of 82 d. After weaning, males and females were grazed in separate paddocks for 8 mo (two seasons). The stocking rates ranged on average from 1.2 t of live weight/ha at the beginning of the period to 1.6 t/ha at the end. Sires pedigrees were traced back to the foundation in 1979. The causes of mortality were recorded. Death due to strongyles was the only source of mortality considered here.
Coccidiosis at weaning and infection by Moniezia during fattening were controlled by regular drenchings. Cowdriosis was completely controlled by twice-monthly accaricide application. Some cases of pneumonia and footrot were diagnosed during the wet season. Affected kids were removed from the experiment. Kids were naturally infected with gastrointestinal strongyles (mainly Haemonchus contortus and Trichostrongylus colubriformis). During the postweaning fattening period, kids were drenched every 8 wk with levamisole (12 mg/kg). Levamisole was effective (Barré et al., 1997
) until the last three cohorts, when worm populations became resistant (G. Aumont, unpublished data).
Fecal samples were collected after 6 and 7 wk of each infection period (the interval between two drenchings) during fattening. Blood samples were collected for each animal every seventh week. Live weights were recorded at weaning, during fattening, at drenching, and in the middle of each infection period.
Fecal egg counts (FEC) were estimated using a modified McMaster method for rapid determination (Aumont et al., 1997a
). In addition, fecal cultures were prepared to assess the composition of nematode burdens. Packed cell volume (PCV) was measured by the capillary microhematocrit method.
Survival Analysis
Theory.
For a complete description of theory, see Kalbfleisch and Prentice (1980)
or Klein and Moeschberger (1997)
. The survival data of infected kids were described with a proportional hazards model, for which the hazard function h(t;xm) at time t of the animal m is written as follows:
![]() | [1] |
where t is time in days since weaning, h0(.) is the baseline hazard function and describes the overall risk of dying of the population, xm(t) is a vector of (possibly time-dependent) fixed covariates or indicator variables with ß as the vector of regression parameters, zm(t) is an incidence vector relating the hazard function to a vector of random (possibly correlated) effects. The function h0(.) can be left completely unspecified (defining the so-called Cox model) or may follow a parametric hazard, the most popular and flexible being the Weibull hazard function, (a two-parameter hazard distribution defined as h0(t)=
(
t)
-1). In order for the model to better fit the data, different baseline hazard functions h0,n(.) can be defined for each level (or stratum) n of a particular factor. The need for different baselines can be assessed using graphical tests. When the plots of the natural logarithm of [-ln
0,n(t)] against the natural logarithm of t (ln t) are parallel lines, a unique baseline can be assumed overall strata. Furthermore, a Weibull proportional hazard model can be assumed if a straight line is obtained.
Model selection.
The fixed effects of sex (male or female), parity of the dam (lactation one or two, lactation three to five, and lactation six and more), cohort (nine levels from 1995-1 to 1997-3), litter size (one, two, or three and more), rearing mode (maternal or artificial), and weight at weaning (as a continuous covariate) were tested using likelihood ratio tests. Only significant effects were included in the final model. The baseline hazard function was initially stratified per level of each effect, except the cohort effect and the effect of weight at weaning. Then the values of ln[-ln
0,n(t)] were plotted against ln t to decide whether these baselines could be grouped together. Treatment application was modeled as a time-dependent variable (level = 1 during the 3 wk following each drenching, level = 2 after the 3 wk following each drenching). The effects of live weight, FEC, and PCV were also tested as time-dependent continuous covariates. Finally, sire effects were added and were assumed to follow a multivariate normal distribution. The sire variance was estimated using a Bayesian approach (Ducrocq and Casella, 1996
). The characteristics (mean, mode, standard deviation) of the approximate marginal posterior density of this sire variance were calculated. All computations were done using Survival Kit version 3.12, a set of FORTRAN programs written with animal breeding applications in mind (Ducrocq and Solkner, 1998
).
| Results |
|---|
|
|
|---|
Model Selection and Fixed Effects Estimates
Baseline survivor functions for the effect of rearing mode were grouped together because the plots of ln[-ln
0,n(t)] against ln t appeared to be roughly parallel. Lines were not parallel for the sex effect, requiring the consideration of two different baselines (Figure 1
). However, these lines are relatively straight, indicating that a Weibull model stratified by sex could be assumed. For each level of sex, the two baseline parameters
and
ln
(= intercept) graphically estimated from the Cox model or directly estimated by maximum likelihood techniques for the Weibull model are presented in Table 1
. They appear to differ when one or the other model is used, but this is partly the consequence of a strong negative correlation between the estimates of
and
ln
. Because fewer a priori assumptions are postulated, only the estimates obtained from the Cox model are presented (Table 2
) and discussed.
|
|
|
> 1) and decreased in females (
< 1). This corresponds to a death rate more than three times greater in males than in females (Figure 2
|
Estimated Survivor Curves
The estimated survivor curves are presented in Figure 2
for different populations characterized by different levels of sex, rearing mode, weaning weight, drenching, and cohort. Whatever the effect, survival rate mainly decreased at three steps: at the beginning of fattening (between 20 and 40 d), around 6 mo of age (between 100 and 120 d), and at the end of fattening (between 220 and 240 d). A proportion of 95% of females were still alive after 166 d vs. 119 d for males. The threshold of 90% of kids still alive was reached after 246 d and 115 d, respectively, for maternal and artificial rearing. After 260 d, the survival rate was about 15% higher for heavy kids (12.5 kg) at weaning than for light ones (6.5 kg). With drenching, 90% of the kids were still alive after 246 d. Without drenching, this mortality rate of 10% was reached 45 d earlier. In the different cohorts, the survival rates 260 d after weaning ranged from 97.1 to 85.5%.
Variance Components and Genetic Evaluation
Genetic parameters obtained with Cox models are reported in Table 3
. As expected, the standard deviations and the skewness of the posterior density of sire effect were large and substantially changed when pedigree information was included. Using the formula of Yazdi et al. (2002)
, the heritability estimated with a sire model, in the unrealistic situation of no censoring, ranged from 0.50 to 0.70. After correction for the high censoring rate, the equivalent heritability (i.e., the heritability that can be used with selection index theory [Yazdi et al., 2002
] to approximate reliabilities of genetic evaluations) was very low (around 0.05), but also very imprecise.
|
| Discussion |
|---|
|
|
|---|
The adverse effect of artificial rearing on survival is probably due to poorer grazing ability and fewer contacts with parasites before weaning. There is no evidence of protection against gastrointestinal nematode infections provided by the ingestion of colostrum from the dam (Dineen et al., 1978
). Therefore, since parity of the dam and litter size were not significant, no maternal influence seems to modify survival of kids infected with gastrointestinal nematode parasites. These results are consistent with the lack of maternal effects reported for resistance criteria in Creole kids after 6 mo of age (Mandonnet et al., 2001
). Postweaning survival and resistance can be evaluated on individual performances only. This suggests a good degree of adaptability to harsh tropical conditions in this local breed.
The inclusion of weight or PCV as a time-dependent covariate in the model eliminated the significance of the other effects. This result indicates that death in infected kids occurs after weakening (anemia) and severe loss of weight. High infection level (egg output) certainly increases the risk of death. However, it is not a direct cause of mortality. Some kids can harbor a high worm burden without dying. Very similar risk, of death per unit of PCV and FEC reported by Nguti et al. (2003)
. They also highlighted the important influence of reduced weight in likelihood of death. This result justifies the consideration of mortality as a resilience criterion compared with egg output, which is a resistance criterion.
Variability of small ruminant mortality is commonly analyzed as a binary trait (Lancelot et al., 2002
; Rege et al., 2002
; Baker et al., 2003
). Currently no data has been published using survival analysis in kids. The major advantages of this methodology are the complete use of the information and the unique ability to incorporate covariates that vary with time, such as treatment, live weight, resistance criteria. Our experimental conditions were very similar to those found in clinical biometrics, where the methodology has primarily been developed, with a low rate of death and small data sets. Finally, it allowed for drawing conclusions on covariates of particular interest. A Cox model was used rather than a Weibull model because fewer assumptions are necessary.
The high censoring rate in our data set had two main consequences. First, it prevented us from correctly testing interactions between effects. Second, the level of mortality following strongyle infection was too low to express significant family differences. No genetic variability on mortality related to strongyles could be formally assessed from our data set on goat. Larger data sets are needed for a more precise assessment of genetic variability, but are not easy to obtain, at least when strongyle infection is to be studied in detail. However heritability estimates of mortality have generally ranged from 0 to 0.1 in sheep (Lopez-Villalobos and Garrick, 1999
; Morris et al., 2000
; Cloete et al., 2001
), even in harsh environmental conditions (Burfening, 1993
; Snyman et al., 1998
). Therefore, little genetic progress should be expected for this critical component of flock productivity. When mortality is studied using survival analysis, higher genetic variability can be assessed (Nguti et al., 2003
). Southey et al. (2001)
reported estimates ranging from 0.1 to 0.2. They promoted selection to improve productivity and profitability of sheep production. Bishop et al. (2002)
also suggested putting emphasis on this trait. Genetic variability in mortality would enable genetic improvement of tolerance to gastrointestinal strongyles infection and would therefore reduce the effects of disease in the flocks. Also, increased tolerance would not place much selective pressure on the pathogen, unlike resistance. At the moment, references on goat mortality are scarce and no genetic parameters have been previously published.
| Implications |
|---|
|
|
|---|
| Footnotes |
|---|
2 Correspondence: INRA-URZ, Prise deau (phone: +33-590-25-54-08; fax: 33-590-25-59-36; E-mail: mandonne{at}antilles.inra.fr).
Received for publication October 17, 2002. Accepted for publication June 3, 2003.
| Literature Cited |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
L. Limea, M. Boval, N. Mandonnet, G. Garcia, H. Archimede, and G. Alexandre Growth performance, carcass quality, and noncarcass components of indigenous Caribbean goats under varying nutritional densities J Anim Sci, November 1, 2009; 87(11): 3770 - 3781. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. C. Bambou, R. Arquet, H. Archimede, G. Alexandre, N. Mandonnet, and E. Gonzalez-Garcia Intake and digestibility of naive kids differing in genetic resistance and experimentally parasitized (indoors) with Haemonchus contortus in two successive challenges J Anim Sci, July 1, 2009; 87(7): 2367 - 2375. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |