|
|
||||||||
ANIMAL GENETICS |
,3


* Area de Producció Animal, Centre UdL-IRTA, 25198 Lleida, Spain;
and
Departament de Ciència Animal i dels Aliments, Universitat Autònoma de Barcelona, 08193 Bellaterra (Barcelona), Spain; and
and
Dto. Reproducción Animal, Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria, 28080, Madrid, Spain
| Abstract |
|---|
|
|
|---|
Key Words: Bayes Factor Neonate Metabolism Survival Vitality
| Introduction |
|---|
|
|
|---|
Early physiological or behavioral variables may also represent an alternative approach by predicting the survival ability of the piglets. After birth, newborn pigs must adapt to the new environment quickly because a delay could create a great disadvantage in terms of access to resources (Fraser and Thompson, 1991
). In addition, Herpin et al. (1996
, 1998)
reported that the degree of asphyxiation suffered during delivery had a significant influence on piglet adaptation and posterior viability. This period of hypoxia also initiates drastic changes in various physiological variables that, which in turn can influence piglet survivability (Herpin et al., 1996
, 1998
; Tuchscherer et al., 2000
). Several authors have related preweaning piglet mortality with early metabolic variables such as body temperature or arterial oxygen saturation (Herpin and Le Dividich, 1998
, Herpin et al., 2001
; Casellas et al., 2004a
), or initial suckling behavior (Herpin et al., 1996
; Tuchscherer et al., 2000
). Moreover, heart rate has been also included in the vitality score (Randall, 1971
).
However, knowledge about the additive genetic variation of these traits is practically nonexistent. One important question is whether there is a genetic background in the observed phenotypic variation. Recently, a Bayes factor (BF) procedure for testing the genetic background of quantitative traits was developed (García-Cortés et al., 2001
; Varona et al., 2001
). In this context, this article is an attempt to examine the additive genetic background of several physiological and vitality variables closely related to the first minutes of life after birth.
| Materials and Methods |
|---|
|
|
|---|
One hundred eight litters from Iberian x Meishan F1 gilts mated with eight Iberian x Meishan F1 boars were recorded in Nova Genètica farm of Solsona (Lleida, Spain) between November 2001 and May 2002. These F1 individuals came from three unrelated Iberian boars and 18 Meishan sows, and they provided 1,039 piglets born alive, 34 stillbirths, and 27 mummified fetuses. Matings were designed to avoid possible inbreeding. These data were obtained in the scope of a F2 design for detection of QTL for reproductive traits (Grant AGL2000-1229-C03). Primiparous sows were housed in climate-controlled rooms (24°C) during farrowing and lactation, and penned in standard farrowing crates with heating plates for piglets (38°C). Feeding of sows was restricted during the gestation period (9.2 MJ of NE, 13.5% CP, and 0.48% lysine), and ad libitum during lactation (9.8 MJ of NE, 17.5% CP, and 0.82% lysine). During parturition, veterinary interventions were kept to a minimum; 20 I.U. of oxytocin was administered by i.m. injection when the interval between births exceeded 30 min. Furthermore, sows that did not deliver before d 113 of gestation received 175 µg of cloprostenol to elicit farrowing (intravulvar injection).
Animal procedures followed ethical norms of IRTA.
Physiological and Vitality Trait Recording
Birth weight, time, and presentation (anterior or posterior) were recorded for each piglet, and the vitality score was determined following in part the procedures described by Randall (1971)
. Three variables were scored by subjective evaluation (see Casellas et al., 2004b
): time to the onset of respiration, muscle tone, and attempts to stand. The sum of these three values defined vitality score. Piglets were monitored for heart rate (HRL), arterial oxygen saturation (OS1), and rectal temperature at birth (RT1) and 60 min later (HR2, OS2, RT2), using a Vet/Ox 4404 monitor pulsoxymeter (Heska Corp., Fort Collins, CO). Given that all the piglets had black skin, heart rate and arterial oxygen saturation measurements were taken by tongue, using an Exotic Reflectance Sensor (Heska Corp.). The time interval between birth and the measurement of OS was also recorded, which was used to adjust arterial oxygen saturation (OS) to 90 s after birth (OS1), applying the GLM procedure of SAS (SAS Inst., Inc., Cary, NC) because Herpin et al. (1998)
observed a large increase in OS during the first seconds of life. All piglets were ear tagged for identification purposes, weighed, and placed near the posterior side of their mother. This allowed for recording the time that elapsed between birth and first teat contact, as well as the first suckle (TS).
Bayes Factor Analysis
For each variable the significant effects previously analyzed in Casellas et al. (2004b)
were considered (Table 1
), with the exception of birth weight (BWT), which was analyzed following in part the fixed model proposed by Roehe (1999)
, with sex and litter size as systematic effects, and a random source of variation characterized by the litter. We estimated the BF between the models with and without genetic effects following García-Cortés et al. (2001)
and Varona et al. (2001)
. The first model (Model 1) includes the additive genetic effects:
|
![]() |
where y contains phenotypic records; X, Z1, and Z2 are the incidence matrices relating observations to systematic (ß), additive genetic (u), and random litter (p) effects; and e is the vector of residuals. The vectors u, p, and e are assumed to be normally distributed.
![]() |
![]() |
![]() |
where
2u is the additive genetic variance,
2p is the variance of the random litter effects,
2e is the residual variance, and A is the numerator relationship matrix. Model 1 can be reparameterized as:
![]() |
where
![]() |
Consequentsly,
![]() |
![]() |
where
2 =
2u/
2t is the intraclass correlation, and
2t = (
2u +
2e) is the variance of e*. The joint distribution of all variables in Model 1 is:
![]() |
where
![]() |
![]() |
and 0 otherwise for each level i of ß)
![]() |
![]() |
![]() |
![]() |
where k1, k2, and k3 are three small enough values to ensure a flat distribution over the parametric space.
The alternative model (Model 2) is:
![]() |
where
![]() |
![]() |
Then, the joint distribution of records and parameters is:
![]() |
where we can assume that prior distributions p2(ß), p2 (p |
2p), p2(
2p), and p2(
2t) are identical to the prior distributions of the previous model. And the likelihood of the model is
![]() |
According to Varona et al. (2001)
, only the analysis with the complex model (Model 1) is required. Following García-Cortés et al. (2001)
and Varona et al. (2001)
:
![]() |
because p1(
2 = 0) = 1 or simultaneously:
![]() |
A BF12 (BF21) sss>1 (<1) indicates that the model with additive variance components (heritable) is more suitable. On the contrary, a BF12 (BF21) <1 (>1) indicates that the non heritable model is more probable. From now on, we will refer to BF12 as BF.
According to Jeffreys (1984)
, the BF can be classified according to the levels of evidence:
|
Sampling from the conditional distribution of
2 was performed using a Gibbs sampler (Gelfand and Smith, 1990
), with a Metropolis-Hastings step (Hastings, 1970
). To calculate the posterior distribution of the heritability (h2), the following formula was used every iteration:
![]() |
A total of 12,500 iterations were performed, and 10,000 iterations were used after discarding the first 2,500. All correlated samples were used to calculate the posterior distributions using the ergodic property of the chain (Gilks et al., 1996
). Convergence was checked using the algorithm of Raftery and Lewis (1996)
.
| Results |
|---|
|
|
|---|
Five variables reached average BF larger than 1: heart rate 60 min after birth (48.90), rectal temperature at birth (13.82), birth weight (2.27), heart rate at birth (1.22), and the interval between birth and first teat contact (1.03). Following the classification of the BF by Jeffreys (1984)
, there was very strong evidence for the presence of genetic background of HR2, and strong evidence in RT1. The posterior modes (means in parentheses) of their heritability were 0.13 (0.22), 0.21 (0.38), 0.29 (0.32), 0.40 (0.39), and 0.17 (0.38) for BWT, heart rate at birth and 60 min later, RT1, and time between birth and first teat contact, respectively (Table 2
). The rest of the variables attained BF <1, varying between 0.56 for elapsed time between birth and first suckle and 0.29 for OS1, all with posterior modes for heritability smaller than 0.13 (Table 2
). The distributions of the heritability estimates were clearly asymmetrical, revealing great differences between the posterior mode and the posterior mean for the heritability, with the exception of the HR2 and RT1 variables.
|
| Discussion |
|---|
|
|
|---|
The application of the BF for testing the genetic background described by García-Cortés et al. (2001)
requires the reparameterization of the model in terms of phenotypic variance and intraclass correlation because the BF can only be computed between nested models that differ in a bounded variable using this approach. In this case, the bounded variable is the intraclass correlation, and the comparison between models is calculated between a model with zero intraclass correlation (nonheri-table) against a model with intraclass correlation between 0 and 1 (heritable).
The main disadvantage of the application of the BF is its strong dependence on the assumed prior distributions. However, in this case, the prior distributions for all parameters, with the only exception of the intraclass correlation, are cancelled because they are assumed to be the same in both candidate models (Varona et al., 2001
). In addition, the final formula only includes the marginal prior and posterior distribution for the in-traclass correlation (
2), and, as a consequence, the BF results are robust for modifications of this prior distribution. When the prior distribution of
2 is modified, the posterior distribution changes in the same direction and, as a consequence, the BF remains stable.
Physiological variables and initial feeding behavior may play an important role during adaptation to the external environment and subsequent survivability in newborn pigs (Fraser and Thompson, 1991
; Herpin et al., 1998
; Casellas et al., 2004a
). Variables like heart rate, rectal temperature, or OS characterize the metabolic situation of piglets after birth, the degree of asphyxiation suffered during delivery, and the ability to recover physiological levels (Herpin et al., 1996
, 1998
; Casellas et al., 2004b
). Notwithstanding, our knowledge of these variables is limited, especially concerning their genetic component. To our knowledge, this is the first study that analyzed the genetic background of heart rate, arterial oxygen saturation, or rectal temperature in neonates. Thus, our results open a new research field related to neonate physiology in general and the genetic aspects of newborn pig metabolism in particular.
Rectal temperature at birth or during the first day of life is historically the best-studied variable among those treated in this investigation (see Herpin and Le Dividich, 1998
), perhaps with the exception of birth weight, but our knowledge about the possible heritability of body temperature is restricted to adult individuals. In this context, the evidence of heritable variation for body temperature has been the object of an important discussion in recent decades. Silcock and Parsons (1973)
observed an additive genetic component for body temperature in mice, but subsequent studies were unable to find significant heritabilities (Lacy and Lynch, 1978
; Lynch et al., 1988
). Recently, a response of rectal temperature to selection was demonstrated in laboratory rats (Gordon and Rezvani, 2001
) but, on the other hand, the heritability estimates for rectal temperature in leaf-eared mice (Phyllotis darwini) did not reach the statistical threshold, although other variables closely related to animal thermoregulation, such as thermal conductance, showed significant heritabilities (Nespolo et al., 2003
). Unfortunately, our results did not allow us to clarify this controversy. Although rectal temperature at birth showed a large BF (13.82), with a moderate to high posterior mode for the heritability (0.40), rectal temperature 60 min after birth reached a lower posterior mode (0.11) and BF (0.42). This implies a great change in heritability over a short temporal space, although Casellas et al. (2004b)
found a significant phenotypic correlation between both measures in the same data set. This results are perhaps comparable with the ones obtained in leaf-eared mice, where diurnal temperature was clearly nonheritable, whereas the genetic component of nocturnal temperature was almost significant (P < 0.10; Nespolo et al., 2003
).
Moderate heritabilities were estimated for heart rate variables, with a BF >1, especially for heart rate 60 min after birth (48.90). In this context, an analysis carried out in adult humans also described a moderate heritability for heart rate (Singh et al., 1999
), and Barrey (1992)
reported a high heritability for the maximum heart rate in adult human athletes (0.72). These results indicate the presence of a genetic component for these variables, compatible with the description of a QTL associated with heart rate in an F2 experiment involving CBA/CaJ and BALB/cJ mice strains (Sugiyama et al., 2002
). However, both heart rate variables we measured have to be considered as different traits because Casellas et al. (2004b)
did not find a significant genetic correlation.
Bayes factors for arterial oxygen saturations did not reach values greater than 1, so there is no evidence for the existence of an additive genetic component for these variables in neonates. It is also remarkable that the BF for OS1 was 0.29, or 3.45 in the opposite sense; thus, there is substantial evidence in favor of the model without genetic determinism. Notwithstanding, a major gene for OS has been described in Tibetan highlanders, obtaining a heritability of 0.41 ± 0.14 (Beall et al., 1994
). Unfortunately, the effect of this major gene has never been evaluated in neonates or in swine. Similar values were obtained for the interval between birth and the first teats contact, as well as the interval between birth and the first colostrum intake. Only the interval between birth and the first teat contact reached a BF >1, but very close to 1. This result implies that the time spent searching for teats may be influenced by an additive component, but further studies are necessary to confirm this hypothesis.
Birth weight reached a moderate heritability (0.13) with a BF of 2.27. It did not provide us with substantial evidence regarding the genetic background of this trait; the posterior mode for heritability was similar to the value reported by Kerr and Cameron (1995)
, and slightly greater than the one described by Roehe (1999)
. However, the results of this study are limited by the structure of the data set, given that the experiment was designed for QTL detection in reproductive traits. The F2 design is not an appropriate design for detecting genetic background because in the F2 generation, genetic components related to dominance and linkage disequilibrium between loci usually appear. The additive genetic variance that we were able to detect in this study was composed of half the genetic variance of the parental populations plus the segregation variance. This latter variance is the one that we shall be able to exploit in QTL analysis to detect specific loci related to the differences in these behavioral and physiological variables. Unfortunately, data from F0 and F1 generations were not available to allow for the discrimination between additive variances from the parental populations and the segregation variance (Birchmeier et al., 2002
). In addition, it is expected that half the heterosis is still present in this F2 generation and, in our opinion, its effect will be absorbed into the systematic effects.
| Implications |
|---|
|
|
|---|
| Footnotes |
|---|
3 Both authors contributed equally to this manuscript. ![]()
2 Correspondencephone: 34 973 702637; e-mail: Luis.Varona{at}irta.es.
Received for publication October 8, 2004. Accepted for publication October 27, 2004.
| Literature Cited |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
X. F. de Sevilla, E. Fabrega, J. Tibau, and J. Casellas Genetic background and phenotypic characterization over two farrowings of leg conformation defects in Landrace and Large White sows J Anim Sci, May 1, 2009; 87(5): 1606 - 1612. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |