|
|
||||||||
Departamento de Mejora Genética Animal, INIA, Ctra La Coruña km 7, 28040 Madrid, Spain
2 Correspondence:
phone: 34-91-3476757; fax: 34-91-3572293; E-mail:
silio{at}inia.es.
| Abstract |
|---|
|
|
|---|
G). Inbreeding depression, expressed as the performance decrease relative to the mean, per 10% increase of the inbreeding coefficient, was -3.18, -2.31 and -5.37% for weight at 120 d (a), and -3.16, -2.08 and -6.49% for daily gain (b) in each of the three strains. The results indicate that the inbreeding effects are dependent on the level of previous inbreeding, the more previous inbreeding, the less inbreeding depression.
Key Words: Bayesian Theory Coancestry Growth Inbreeding Depression Pigs
| Introduction |
|---|
|
|
|---|
In livestock species, evidence for the inbreeding effects comes partially from data recorded in progressively inbred populations (Burrow, 1993). This procedure poses different problems that can be avoided when data are available from contemporary unselected individuals inbred to various degrees. The goal of this experiment was to analyze the effect on growth traits of a fast increase of inbreeding in three related strains of pigs with a different past level of inbreeding. To obtain accurate inferences about inbreeding effects, data were analyzed using a Bayesian procedure, which takes into account all the available information and allows a joint analysis of parameters of growth functions, (co)variance components, systematic effects, and breeding values (Varona et al., 1997).
| Material and Methods |
|---|
|
|
|---|
Data analyzed in this paper come from three related strains of Iberian pigs maintained in a conservation program established from four Spanish and Portuguese ancient populations of this Mediterranean breed: a) Torbiscal, which is the result of blending the four ancient founder populations; b) Guadyerbas, which is one of these founder populations, representing the ancient earliest mature black hairless type of the Iberian breed; and c) Gamito, a hairless strain developed exclusively from Guadyerbas but showing a red coat color. These populations of different effective size have previously experienced similar mild selection pressure (Béjar et al., 1993).
An experimental trial was carried out to test the growth performance of contemporary groups of inbred and control pigs of the three strains. A combination of full-sib and minimum coancestry matings were planned for each of the available boars of each strain. Two full-sib piglets (1 male and 1 female) were sampled at weaning from each of the 116 litters born in four testing batches and each pair of full-sib pigs was separately penned and weighed monthly. Animals were fed twice daily with a diet containing 3,450 kcal of digestible energy and 186 g of crude protein per kilogram of DM. The average daily ration per pig increased from 0.75 kg up to a maximum value of 3.25 kg. The number of boars, dams, and progeny tested in the inbred and control groups of each strain are presented in Table 1
. The growth data consisted of 1,615 weights recorded on 229 pigs between 50 and 250 d of age. Pigs were weighed 3 to 8 times. The average pig weight increased from 13.8 to 93.5 kg during the trial.
|
The genetic relationships among all the animals were known back to the founder pigs of the herd and diverse pedigree tools can monitor the loss of genetic variability throughout the history of these populations (Rodrigáñez et al., 2000; Toro et al., 2000). The basic concepts in the pedigree analysis are the coancestry (kinship) coefficient between two individuals fxy (probability of identity by descent of two genes taken at random from each individual at the locus) and the inbreeding coefficient Fx of an individual (probability of identity by descent of the two genes carried by this individual at a given locus). Both concepts are closely related with the concept of genetic contribution (ci) defined by James and MacBride (1958) as the proportion of all pathways leading back from descendants that pass through the individual ancestor i.
Caballero and Toro (2000) have shown that in a pedigree with a total of M individuals, where N0 are founders (animal with unknown parents), the average pairwise coancestry (including reciprocals and self-coancestries) f of a given group of individuals can be partitioned into components of coancestry due to the genetic contributions of founders:
![]() |
and nonfounders:
![]() |
Another parameter related to genetic contribution is the number of discrete generation equivalents
This parameter is equal to the number of generations in a pedigree with no overlapping generations (Woolliams and Mäntysaari, 1995). Under random mating, inbreeding is delayed with respect to coancestry for one generation. The departure from random mating can be measured with the coefficient
that indicates the degree of deviation from Hardy-Weinberg proportions and it is related to the previous coefficients by (1 - F) = (1 - f)(1 -
). A positive value of
means that relatives are preferentially mated.
The Growth Function Approach
Varona et al. (1997) proposed a Bayesian procedure to analyze growth performance data from production functions that has been applied to these data. In the analyzed growth period, the weight of tested pigs is adequately modeled by a linear function. The two parameters of this function, a (intercept) andb (slope), correspond to productive traits: weight at 120 d of age and daily gain, respectively. The kth weight of the jth animal recorded on day xjk - 120 is considered as a sample from the normal distribution:
![]() |
where
is the variance of weight given the production function parameters. The analysis of the growth function parameters requires other variables related to genetic and environmental relationships between the growth functions of different animals. These sources of variation were represented by the following animal model:
![]() |
where
= matrix of values of intercept and slope (a and b) of orderN x 2; ß = [ßa, ßb], u = [ua, ub], c = [ca, cb] and e = [ea, eb] are matrices of systematic, additive genetic, common litter environmental and residual effects for both parameters;G, C, and R are their respective (co)variance matrices; andX, Z1, and Z2 are known incidence matrices. The systematic effects considered in ß were two sexes, four test groups, and six experimental groups (three strains x two degrees of inbreeding).
The joint posterior distribution given the weights can be explained as their likelihood given the growth function parameters multiplied by the prior distribution of these parameters. In this application, the likelihood is equal to
![]() |
![]() |
with N being the number of pigs with records and nj the number of weight records of the jth animal.
The following prior distribution of the linear growth function parameters was assumed:
![]() |
![]() |
The prior distributions of breeding values and common litter environmental effects given the corresponding (co)variances are, respectively, f(u | G) ~ N (0, A
G) andf(c | C) ~ N (0,I116 C), where A is the numerator relationship matrix. The other prior distributions:
f(G), f(C), (R) and f (ß) are taken as flat, as is usual in Bayesian analysis to reflect a vague previous knowledge about these parameters.
As a consequence, the joint posterior distribution given the performances was
![]() |
![]() |
![]() |
Gibbs Sampler
The marginal posterior density of each unknown parameter was determined using a Gibbs sampling algorithm (Wang et al., 1994), which requires independent draws from all the full conditional posterior distributions. Conditional posterior distributions are univariate normal both for the parameters of linear function (a, b) and systematic, genetic, and common litter environmental effects (ßa, ßb, ua, ub, ca, cb). Details about all these distributions can be found in Rodríguez et al. (1996) and Varona et al. (1997). The conditional posterior distributions of the (co)variance matrices (G, C, R) are bivariate inverted Wishart distributions
and
, whereN, and nu, andnc are, respectively, the number of pigs with records, animals in the pedigree, and the number of litters and Su,Sc, and Se parameters are 2 x 2 matrices with the following expressions:
![]() |
![]() |
![]() |
The convergence was assessed by the double-chain method (García Cortés et al., 1998). The length of the Gibbs sampling was 2,050,000, discarding the first 50,000 (warm-up) and saving only one sample from each 200 iterations. The effective number of samples was calculated from the estimated autocorrelations using time series methods (Geyer, 1992) and ranged from 400 to 2,700 for the different parameters resulting in low values for Monte-Carlo standard errors. Marginal densities were estimated from the Gibbs samples using nonparametric kernel estimators (Silverman, 1986) and the main statistics of the posterior distributions were computed from these densities.
Inbreeding Effects
The marginal posterior distributions of the differences between the effects of the inbred and control groups on the growth function parameters (a, b) allow inferences to be made about the inbreeding effects on growth performance in the three strains of Iberian pigs. The differences
and
were calculated for each strain S in each iteration. The following scaled expressions were also calculated:
![]() | [1] |
![]() | [2] |
These expressions scale the differences between average a and b of inbred and selected group of each strain S to a 10% increase in average inbreeding coefficient F [1] and additionally to the different mean value of intercept and slope of each strain [2]
. Finally, the differences of the other two strains from the Torbiscal strain for the expressions [1]
and [2]
were also calculated.
| Results and Discussion |
|---|
|
|
|---|
The average values of coancestry (f) and inbreeding (F) coefficients of the pigs of the six experimental groups and of their respective parents are presented in Table 1
. With an unequal number of males and females, the inbreeding and coancestry coefficients of the inbred and control groups can be calculated as an unweighted or weighted average of the two sexes. Both values were close in the analyzed cases and the unweighted coefficients are shown. The F values are dependent on the mating system. As expected, the inbred groups present positive values of
. The matings of breeding animals in control groups resulted in negative values for
, reflecting the fact that matings between close relatives were avoided, as was done in the development of these lines. The genetic diversity of a population (GD = 1 - f) can be partitioned into diversity between and within individuals. Caballero and Toro (2000) have shown that the proportion of genetic diversity between individuals can be expressed by a coefficient closely related to
:G = (1 +
)/2. The values of G for the inbred progenies of Guadyerbas, Gamito, and Torbiscal (0.581, 0.589, and 0.598) confirm that with full-sib matings, most of the genetic diversity is between individuals. However, the values of G for the correspondent control progenies (0.471, 0.470, and 0.476) indicate a greater proportion of genetic diversity within individuals.
There are clear differences for f and F values among the three strains derived from the different effective sizes: Torbiscal, with a broader origin and higher number of breeding animals, have lower values for f and F than the other two strains. The complete pedigrees of the Guadyerbas, Gamito, and Torbiscal strains include 751, 770, and 2,884 breeding animals, respectively. To evaluate the progression of coancestry in these strains and the respective contribution of founder and nonfounder individuals, mean values for breeding animals born every 3 yr were calculated. The results of the pedigree analysis for the first 36 yr are shown in Figure 1
. The average value of coancestry coefficient f increased in this period from 0.052 to 0.318 in Guadyerbas, from 0.094 to 0.377 in Gamito, and from 0.015 to 0.133 in Torbiscal. The corresponding value of F increased from 0.000 to 0.278 in Guadyerbas, from 0.000 to 0.348 in Gamito, and from 0.016 to 0.125 in Torbiscal.
|
Variance Ratios and Estimates of Linear Growth Parameters
The Bayesian procedure provides estimates of growth function parameters of each animal that take account of their own performances and, through prior distributions, the information from relatives and other unrelated animals sharing the same systematic effects (Varona et al., 1997). Moreover, this approach provides an estimation of the variance of performances given the production function parameters
, that includes errors in measure or physiological variations of the animals along the fattening period. This procedure can be improved by modeling the increase of
with the age of animals, as has been done in the analysis of growth data from chickens by Mignon-Grasteau et al. (2000). This more complex approach was discarded for the present analysis given its higher computational cost and the moderate value of the scale factor between the first and last weights recorded in the trial.
The main statistics of the marginal posterior distributions (m.p.d.) of
and the heritabilities, common environmental coefficients for the intercept (a, weight at 120th d) and slope (b, rate of daily gain) of the growth linear functions of tested pigs are summarized in the Table 2
. Point estimates of heritabilities for a (0.576) and b (0.547) were lower than those obtained in related strains of Iberian pigs using this Bayesian procedure, but the point estimate of genetic correlation (0.989) was very similar (Fernández et al., 2002). However, common litter environmental effects, which accounted for 6 to 8% of the variance in the present study, were not taken into account.
|
Inbreeding Effects on Intercept and Slope
The trial was designed as a comparison of the effects on growth of discrete levels of inbreeding, generated by contemporary control and full-sib matings in three populations of different previous inbreeding. The purpose was to avoid some of the problems of most of the studies of inbreeding depression based on records from mildly inbred populations of domestic animals: 1) environmental effects can be confused with inbreeding effects, when the different classes of inbreeding are assayed in different generations and environments (Lynch and Walsh, 1998) and 2) the use of pedigree-based inbreeding coefficients to estimate inbreeding depression for traits under either direct or indirect selection is potentially biased (Groen et al., 1995).
The inbreeding effects on the linear growth parameters are the main parameters of interest of this study. The m.p.d. of the effect of differences between inbred and control pigs on the a and b in the three strains are presented in Tables 3
and 4
. Absolute and relative values of effects of inbreeding differences between groups have been considered, because inbreeding depression values are frequently scaled to the effect of a 10% increase in inbreeding on the mean trait of noninbred animals. The conventional techniques of animal breeding only allow approximate solutions to the scaled effects. The use of Gibbs sampling permits drawing of m.p.d. of functions of estimated parameters (as the scaled effects) expressing its uncertainty (Blasco, 2001). According to the values of 95% highest posterior density intervals (95%HPD), negative inbreeding effects on the intercept and slope of the growth function are clearly evidenced in Torbiscal but not in Guadyerbas and Gamito strains. In the last population, the 95%HPD confidence intervals of absolute or scaled inbreeding effects on growth parameters included zero and positive values. For a similar trait, Bereskin et al. (1968) have reported a decrease of 4.3% of the mean of the weight at 154 d per 10% of increase of inbreeding. This value is close to the point estimates of scaled inbreeding effects on weight at 120 d obtained in this study (from 2.45 to 5.54% of the mean).
|
|
|
|
| Implications |
|---|
|
|
|---|
| Footnotes |
|---|
Received for publication August 22, 2001. Accepted for publication April 3, 2002.
| Literature Cited |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M.-H. Li, I. Stranden, and J. Kantanen Genetic diversity and pedigree analysis of the Finnsheep breed J Anim Sci, May 1, 2009; 87(5): 1598 - 1605. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. I. Fernandez, E. Alves, A. Fernandez, E. de Pedro, M. A. Lopez-Garcia, C. Ovilo, M. C. Rodriguez, and L. Silio Mitochondrial genome polymorphisms associated with longissimus muscle composition in Iberian pigs J Anim Sci, June 1, 2008; 86(6): 1283 - 1290. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |