|
|
||||||||
Area de Producció Animal, Centre UdL-IRTA, Rovira Roure 177, 25198 Lleida, Spain
2 Correspondence:
phone: 34-973-702576; fax: 34-973-238301; E-mail:
joseluis.noguera{at}irta.es.
| Abstract |
|---|
|
|
|---|
Key Words: Bayesian Theory Pigs Prolificacy Selection Responses
| Introduction |
|---|
|
|
|---|
Two alternative approaches have been proposed for increasing response to selection for litter size. One is the "hyperprolific" selection scheme, based on the application of a high selection intensity over a large population of females evaluated with an individual selection index (Legault and Gruand, 1976), and the other is the use of a family selection index in the genetic evaluation of selection candidates (Avalos and Smith, 1987).
Another alternative is to combine both approaches, the use of BLUP (utilizing all family information) together with intense selection. The aim of the present paper is to report on the results of a large-scale selection experiment to increase litter size in a Landrace pig population. The experiment was carried out using both family information and high selection intensity. The first step was to evaluate the efficiency in choosing individuals for selected and control lines using standardized selection differentials. The second step was to quantify the success of the selection experiment by estimating direct genetic responses for number of piglets born alive. Finally, correlated selection responses for production traits (weight and backfat thickness) were estimated.
| Materials and Methods |
|---|
|
|
|---|
Genetic evaluation for number of piglets born alive (NBA) was performed using a repeatability animal model (h2 = 0.06, r = 0.13) that included herd-year-season and parity as fixed effects. Breeding values and permanent environmental effects were considered random. The breed of the service boar (Landrace vs Large White) was taken into account implicitly by the herd-year-season effect. Genetic parameters were previously estimated by Restricted Maximum Likelihood using all previously available data.
Once the 3,034 sows were evaluated, two genetic lines of 160 sows and 25 boars were established from this population (Figure 1
). One selected (H) and one control line (C) were established. From the 3,034 female candidates evaluated, the control line C was established by choosing 160 females at random. The 25 Landrace boars currently used in the population were also used in line C (MC). After this, the 160 sows with the highest breeding value for NBA were used to make the H line. Boars in line H (MH) were selected from the 25 litters (one male per litter) with the highest breeding values for NBA from the 961 Landrace litters available in the selection herds. All these animals in both lines constituted the generation G0. Sows from line H and C of G0 had an average of 3.09 and 2.31 parities, respectively.
|
Response to selection for NBA in the experiment was estimated with the offspring of G0 in both lines, H and C. In line H, 204 daughters and 24 sons were chosen randomly among the offspring of the mating between 117 sows and 24 boars (one son per boar) to constitute G1. In line C, 152 daughters and 23 sons were chosen randomly among the offspring of the mating between 123 sows and 23 boars (one son per boar) to constitute G1. Generation G1 in lines H and C were kept under the same intensive conditions. Gilts were mated at around 9 mo of age. No culling based on litter size was performed, and sows were kept in the herd until they were culled for other reasons (diseases, leg problems, etc.). A random sample of the G1 offspring from each line, H and C, was recorded for weight (WT) and backfat thickness (BF) at around 175 d of age. A new generation (G2) was formed for both lines to add information that is mostly captured in G1 parents. This was decided because of the limited capacity of the avalaible facilities, which made it impossible to increase the size of the experiment at G1. Males and females were selected randomly with the restriction that one son was kept from each boar. The same management as previous generations was followed. A detailed description of the selection experiment data used in this study is given in Tables 1
and 2
.
|
|
The analysis of the experiment was performed using a multivariate Bayesian analysis. Mixed-model techniques do not provide exact inference about genetic change when variances are unknown. In contrast, Bayesian inference, by means of Markov Chain Monte Carlo methods, provides a full description of selection response through its marginal posterior distribution (Sorensen et al., 1994).
Model of Analysis
NBA was analyzed for the first six parities; considering each parity as a different trait. The following multivariate model was assumed:
![]() |
where
Four classes were considered for each herd and year: December to February, March to May, June to August, and September to November summing up to 226 levels per parity. Age of farrowing ranged between 280 and 530 d for the first parity, 433 to 684 d for the second parity, 544 to 846 d for the third, 693 to 990 d for the fourth, 883 to 1,140 d for the fifth, and 975 to 1,288 d for the sixth. In each parity 10 balanced classes were considered, with approximately 10% of the total of farrowings per class.
For WT and BT, the model of analysis was:
![]() |
where
i
Bayesian Analysis
A multivariate implementation of the above animal models was performed for the whole Landrace pig population, as described in Noguera et al. (2002). Joint posterior distributions of data and parameters were constructed by multiplying the multivariate normal likelihood of observed and missing data by the prior distributions. These prior distributions were multivariate normal for breeding values and permanent environmental effects, flat for systematic effect, and inverted Wishart for (co)variance components. A Gibbs sampling algorithm (Geman and Geman, 1984; Gelfand and Smith, 1990) with a data augmentation step (Tanner and Wong, 1987) was carried out. Convergence was checked using the algorithms of Raftery and Lewis (1992) and García-Cortés et al. (1998). Effective sample size was computed following the algorithm of Geyer (1992). Estimates for genetic parameters were reported in a companion paper (Noguera et al., 2002).
Standardized Selection Differentials
The standardized selection differentials were calculated as a posteriori measures of the efficiency of the selection process. They were calculated from the difference between breeding values of selected individuals in G0 and the breeding values of the base population, and its magnitude is expressed in terms of genetic standard deviations.
From the Gibbs sampling output at every iteration, the following statistics were calculated for NBA for the i-th parity:
![]() |
![]() |
![]() |
where
hm(i) and
hf(i) are the means of breeding values of males and females in G0 of line H weighted by the number of offspring of every individual in G1,
pm(i) and
pf(i) are the means of breeding values of males and females that were candidates for selection in generation G0,
cf(i) is the weighted mean of breeding values of females for G0 of line C, and
u(i) is the genetic standard deviation calculated following Sorensen et al. (2001). From that,
Shm(i) and
Shf(i) are random samples from the posterior marginal distribution of the selection differentials for males and females producing progeny in line H in the i-th parity. Moreover,
Scf(i) is a random sample from the posterior marginal distribution of the selection differential for females in the creation of line C in the i-th parity. However, all available males (25) at the moment of creation of lines were used in line C. Thus, the selection differential for males in the i-th parity (
Scm(i)) in line C should be close to zero, except for unequal contributions of males to the next generation caused by chance.
From the samples of selection differentials, density estimation techniques (Silverman, 1986) were used to calculate the posterior marginal densities.
Response to Selection
From every iteration of Gibbs sampling, the following statistics were computed for NBA in each parity:
![]() |
where
h1is the posterior mean of breeding values of individuals of G1 in line H for the i-th parity,
c1(i) is the posterior mean of breeding values of individuals of G1 in line C for the i-th parity, andR(i) is a random sample of the marginal posterior distributions of response to selection for that parity.
For correlated responses in WT and BT, the following statistics were also computed:
![]() |
![]() |
where
h1(WT) and
h1(BT) are the posterior means of breeding values of individuals of G1 in line H for WT and BT, respectively;
c1(WT) and
c1(BT) are the posterior means of breeding values of G1 individuals in line C for WT and BT, respectively. Thus,R(WT) and R (BT) are random samples of the marginal posterior distributions of response to selection for WT and BT, respectively.
From the samples of R, density estimation techniques (Silverman, 1986) were used to calculate the posterior marginal densities.
| Results |
|---|
|
|
|---|
The effective sample size computed with the procedure of Geyer (1992) ranged from 40.25 to 4,231.5. Posterior means and standard deviations were directly estimated from the mean and standard deviation of the samples.
Phenotypic Data
Phenotypic means and standard deviations for NBA for the different parities and generation in each line are presented in Table 1
. Line H showed values higher than line C in all parities and generations except of the sixth parity of G1 and the fifth parity of G2. Ages at farrowing were very similar for both lines in each generation. Phenotypic means and their standard deviations for WT and BT in the G1 are presented in Table 2
. Average weight and backfat thickness were very similar in both lines.
Standardized Selection Differentials
Marginal posterior means and posterior standard deviations (PSD) of standardized selection differentials for G0 in all parities are presented in Table 3
. Posterior means of standardized selection differentials for females in line H ranged from 0.70 (PSD 0.12) in the sixth parity to 0.94 (PSD 0.06) in the fifth. For males in the line H, posterior means ranged from 0.22 (PSD 0.19) in the sixth parity to 0.34 (PSD 0.10) in the fourth. In all cases except the sixth parity in males, highest posterior density regions (HPD) of 95% did not include zero.
|
Direct Genetic Response to Selection for NBA
Marginal posterior distributions of direct genetic responses for NBA are presented in Figure 2
. All distributions were nearly symmetric. Therefore, posterior means, modes, and medians were very similar (Table 4
). Posterior means of direct genetic response for NBA ranged from 0.32 (PSD 0.08) in the first parity to 0.64 (PSD 0.08) in the fourth one.
|
|
|
|
| Discussion |
|---|
|
|
|---|
Our results are in agreement with Le Roy et al. (1987) and with Bolet et al. (1989), who observed a general tendency to increase selection efficiency in later parities. They are also in agreement with Tartar (1981), who showed theoretically that selection over an average of the first five parities results in greater response in later parities. Moreover, average response for the six parities (0.46) is very similar to the expected response (0.41) obtained in a simulation experiment although a repeatability animal model was assumed (Noguera et al., 1994).
The classical "hyperprolific" selection scheme (Legault and Gruand, 1976), based on high selection intensity for litter size applied to a large population of females through individual selection index, resulted in an initial increase of the female side of around one piglet (Bolet and Legault, 1982). After 20 yr of selection, accumulated response was around 1.4 piglets per litter (Bidanel et al., 1994; Herment et al., 1994), but annual genetic gain was smaller (0.07 piglets per year) because of the need for backcrossing boars from the progeny of dams with extreme prolificacy to other sows with extreme prolificacy during repeated cycles of selection.
Higher rates of genetic improvement in litter size when selection is based on information from relatives have been reported. Thus, Wang et al. (1994) reported an experiment to improve total number born per litter with a selection criteria consisting of breeding values predicted with BLUP. The rate of genetic response was about 1.6% per year. More recently, Sorensen et al. (2000) reported a direct genetic response to selection of 0.43 piglets in a large-scale selection experiment for total number of piglets born. The criterion of selection was the average predicted breeding values (BLUP) of the sire and dam of the litter that contributed piglets to the selected line. The average response in litter size in the experiment reported in this paper was in line with these previous experiments.
The calculation of standarized selection differentials allows evaluation of the efficiency of a selection process given the available variance. As observed in Table 3
, selection was effective in both males and females in line H, but greater values were observed for females. Sows were selected using their own information plus their relatives information. However, males were selected based on information of their ancestors only. On the contrary, the HPD of 95% for selection differentials for females in line C always included zero. Although selection was performed at random, culling of individuals resulted in slight positive values of selection differentials. Selection differentials for males in line C should be around zero, because all available males in the population were used to create G1.
Correlated responses in WT and BT were not observed in our experiment and are also consistent with results of other authors (Brien, 1986; Kuhlers and Jungst, 1991; Kerr and Cameron, 1996). It is consistent with estimates of genetic correlations close to zero between litter size and WT or BT reported by Noguera et al. (2002).
The results of the scheme proposed demonstrate an important response of litter size and do not produce undesirable correlated response in production traits. Nevertheless, a disadvantage of "hyperprolific" schemes can be a genetic lag for production traits due to selection pressure on litter size (Bichard and Seidel, 1982; Bidanel et al., 1994).
| Implications |
|---|
|
|
|---|
| Footnotes |
|---|
Received for publication September 7, 2001. Accepted for publication June 6, 2002.
| Literature Cited |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. Rosendo, L. Canario, T. Druet, J. Gogue, and J. P. Bidanel Correlated responses of pre- and postweaning growth and backfat thickness to six generations of selection for ovulation rate or prenatal survival in French Large White pigs J Anim Sci, December 1, 2007; 85(12): 3209 - 3217. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Varona, D. Sorensen, and R. Thompson Analysis of Litter Size and Average Litter Weight in Pigs Using a Recursive Model Genetics, November 1, 2007; 177(3): 1791 - 1799. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Rosendo, T. Druet, J. Gogue, L. Canario, and J. P. Bidanel Correlated responses for litter traits to six generations of selection for ovulation rate or prenatal survival in French Large White pigs J Anim Sci, July 1, 2007; 85(7): 1615 - 1624. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Estany, D. Villalba, J. Tibau, J. Soler, D. Babot, and J. L. Noguera Correlated response to selection for litter size in pigs: I. Growth, fat deposition, and feeding behavior traits J Anim Sci, October 1, 2002; 80(10): 2556 - 2565. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Estany, D. Villalba, M. Tor, D. Cubilo, and J. L. Noguera Correlated response to selection for litter size in pigs: II. Carcass, meat, and fat quality traits J Anim Sci, October 1, 2002; 80(10): 2566 - 2573. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |