|
|
||||||||
GENETICS |
,2
ski*
* Department of Mathematical and Statistical Methods, and
Department of Genetics and Animal Breeding, August Cieszkowski Agricultural University of Poznan, PL60-637 Poznan, Poland
2 Corresponding author: tomasz{at}jay.au.poznan.pl
| ABSTRACT |
|---|
|
|
|---|
Key Words: Gibbs sampling laying hen reproductive trait single gene threshold model
| INTRODUCTION |
|---|
|
|
|---|
Selection intensity may increase when information on the existence of single loci is included into genetic improvement programs (Lahav et al., 2006). At the moment, many major genes affecting reproductive ability are identified for livestock, mainly for sheep (Davis, 2005), cattle (Cobanoglu et al., 2005), and pigs (Rathje et al., 1997). Druyan and Cahaner (2007) suggested that 2 complementary dominant genes affect ascites resistance and ascites susceptibility in broilers.
Bayesian marker-free segregation analysis seems to be a very useful tool for detection of single loci (Kadarmideen and Janss, 2005). This approach supplies information on hypothetical genotypic effects and their frequencies. Both fertility and hatchability are recorded as binary traits. These traits are influenced by many genetic and environmental factors, so the threshold animal model seems to be the most adequate for their statistical analysis (Sørensen et al., 1995; Bennewitz et al., 2007; Skotarczak et al., 2007). In this model, it is assumed that a categorical distribution of observations is determined by an unobserved continuous variable called liability. For the binary recorded traits, the threshold between 2 categories is set in point 0 and the variance component for liability is 1.
The objectives of this study were to verify the hypothesis on segregation of single genes (1 vs. 2) affecting fertility and hatchability and to estimate heritability of these traits.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Data of 2,040 and 2,015 dams from full-pedigreed strains of Rhode Island Red (R33) and New Hampshire (N88) from 1 pedigree farm located in Poland were analyzed.
The birds were naturally mated and kept on litter. The environmental conditions (e.g., feeding level) did not change considerably over time. Both populations were kept under long-term selection. The within-strain generation selection was based on the classical selection index described by W
yk (1978). The five following production traits were included in the selection: initial egg production (until the 35th week), average egg weight between the 33rd and 35th week of age, body weight at the 18th week of age, and age at sexual maturity. Fertility was checked by candling on the eighth day of incubation. The percentage of fertilized eggs and the percentage of the eggs hatched of fertilized eggs were registered for dams only. These reproductive eggs were collected between the 45th and 54th week of age. In a binomial phenotypic scale, 10 eggs per dam were included into the analysis. The number of reproductive eggs studied was restricted by the algorithm applied.
The data were classified according to generation (8 levels) and hatch period (4 levels). Averages for studied traits fluctuated negligible over generations and hatch periods. Therefore, they are not shown in the present study. The brief description of the data sets is given in Table 1
.
|
| Methods |
|---|
|
|
|---|
![]() |
![]() |
where u = a (10s x 1) vector of unobserved liability; X = a (10s x b) design matrix of nongenetic effects; β = a (b x 1) vector of fixed effects of hatch period (4 levels) and generation (8 levels); Z = a (10s x q) design matrix relating polygenic and single locus effects to observations; W and V = (q x 3) random matrices containing information of genotype of each individual, each row of these matrices has 1 of the following forms: [1, 0, 0], [0, 1, 0], or [0, 0, 1] corresponding to genotypes AA, Aa, or aa and BB, Bb, or bb, respectively; µ = [µAA, 0, –µAA]' and
= [
BB, 0, –
BB]' = the vectors in which the first element is the effect of dominant homozygote (e.g., AA), the effect of heterozygote (e.g., Aa) is assumed to be equal to 0 and the effect of recessive homozygote (e.g., aa) is the opposite of the effect of dominant homozygote; a = a (q x 1) vector of random additive polygenic effects; S = a (10s x s) design matrix relating environmental effects to observations; p = a (s x 1) vector of permanent environmental effects; and e = a (10s x 1) vector of random errors effects. Moreover, s denotes the number of recorded hens, and q is equal to the number of individuals included in the pedigree.
The Bayesian methods with Gibbs sampling algorithm have been used to the estimation of unknown parameters in the above introduced models. The known formulas for the ordinary threshold animal model have been adapted for the case when repeated observations were collected for 1 individual. Improper prior uniform distributions were assumed for vectors β, µ, and
. The following multivariate normal distributions were assumed for the random vectors: a
N(0, A
a2), where A = a (q x q) relationship matrix, p
N(0, Is
p2), e
N(0, I10s). The inverted
2 distribution was taken as a prior distribution of variance components
a2 and
p2.
Because it is well known, the conditional posterior distributions for vectors β, µ,
, a, and p in the presented models are normal with means equal to the solutions of the appropriate mixed model equations (Sørensen and Gianola, 2002). In the considered case, the number of the mixed model equations is 10 times higher than the number of observed individuals. Therefore, the following formulas were used to calculate the expected values and variances in every step of Gibbs sampling procedure; for the vector of fixed effects β:
![]() |
where xi = the ith column of X; X–i = matrix X without the ith column; β–i = vector β without the ith element, i = 1,...,b; for the vector of additive genetic effects a:
![]() |
where zi = the ith column of Z; Ai,i = the element of A–1 in the ith row and ith column; Ai,–i = the ith row of matrix A–1 without the ith element; a–i is vector a without the ith element, i = 1,...,q. In both formulas, the element ZV
should be used in model 2 only.
Further, for the elements of vectors µ and
, the following assumptions were made:
![]() |
where w1' and w3' = the 1st and 3rd column of matrix w1 and (zw)1' and (zw)3' = the 1st and 3rd column of matrix ZW, respectively.
![]() |
where v1 and v3 = the 1st and 3rd column of matrix V and (zv)1 and (zv) = the 1st and 3rd column of matrix ZV. In every step of Gibbs sampling procedure, it was assumed that µAA > 0 and
BB > 0.
The following conditional posterior normal distribution was assumed for the unobservable liability: ui
N[xiR β + (zw)iR µ + (zv)iR
+ ziRa + siR p; 1], where xiR = the ith row of matrix X; (zw)iR = the ith row of matrix ZW; (zv)iR = the ith row of matrix ZV; ziR = the ith row of matrix Z, i = 1,...,10s. The expression (zv)iR
should be omitted in model 1.
Moreover, the following inverted
2 distributions were assumed for variance components:
![]() |
and
![]() |
where va,Sa2,vp,Sp2 = the hyperparameters (
a and
p are degrees of freedom, Sa2 and Sp2 are scale parameters).
According to the suggestions of Guo and Thompson (1994), the elements of the unknown genotypes table G were generated from the following formula:
![]() |
where Gi = the genotype of the ith individual; G–i = a table of the genotypes of all individuals excluding the ith individual; Goi = the genotype of the progeny of the ith individual; Gmi = the genotype of the ith individuals mate; Gsi, Gdi are the genotypes of the ith individuals parents, i = 1,...,q. When the individual is not observed, the last term will be substituted by 1. Moreover, the element ZV
was used in model 2 only. It should be stressed that in model 2, the table G was generated for each matrix W and V separately. In both models in the first step of Gibbs sampling, it was assumed that matrices W and V have the following form: W = V = [0:1:0]. Further, for the major gene, Mendelian transmission probabilities were assumed.
To estimate the genotypes for the individuals, the frequencies of alleles among the founders groups are required. The frequencies were generated from the beta distribution according to the following formulas (Kadarmideen and Janss, 2005):
![]() |
![]() |
where nA and na = the number of alleles A and a and nB and nb = the number of alleles B and b, respectively, in the group of founders.
The significance of the major gene effects was checked on the basis of the highest posterior density regions (HPDR; Scott, 1992), which were built for the major gene variance components
GA2 = 2fA(1 – fA)µAA2 and
GB2 = 2fB(1 – fB)µBB2. For example, if the 95% HPDR for
GA2 included value 0, it was stated that the effect of AA genotype is not significant.
Moreover, additional analyses in a submodel were conducted. The submodel did not contain the major genes effect (i.e., it had the following form: u = Xβ + Za + Sp + e). The estimation of parameter β, a, p, u,
a2, and
p2 in the submodel was carried out in a similar way as in models 1 and 2, but in the formulas for conditional posterior distributions, the components containing µ and
were omitted.
On the basis of the estimated variance components, the following genetic parameters were calculated:
![]() |
![]() |
![]() |
![]() |
In the above formulas,
GB2 is omitted in model 1, and
GA2 and
GB2 are omitted in the submodel.
The simulation studies were performed previously to verify the properties of the implemented algorithm (Skotarczak et al., 2007).
In each analyzed case, 1,000,000 rounds of Gibbs sampling were conducted. The first 400,000 steps were discarded as a burn-in period. The important results were collected from every 10th iteration. The means of the posterior distributions were calculated as the point estimators of the unknown parameters.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
|
|
|
|
As mentioned above, the segregation of 1 major gene is suggested for hatchability in strain R33. This trait has a complex background influenced by many genetic and environmental factors. Hatchability is a sum of embryo survivals in a certain period of time—the present study from the eighth to the 21st day of incubation. In general, the ratio of embryonic mortality for chicken is relatively high (Liptoi and Hidas, 2006). It can be strongly affected by recessive genotypes and chromosomal abnormalities. A number of genes determining mortality of embryos were presented by Liptoi and Hidas (2006). It should be noted that the majority of them were identified in first decades of the 20th century, for instance the so-called creeper gene. Hence, they were usually rejected from populations. However, in the last decades, new mutations have been also identified. Delany et al. (1994) found the variation of rRNA gene copy number, which affected early embryonic mortality in chicken. A number of lethal and sublethal genes were described for other animal species. So, the arguments presented above seem to confirm the existence of single loci for hatchability.
It should be recalled that fertility was examined on the eighth day of incubation. Therefore, the percentage of fertilized eggs is underestimated by early embryo mortality. Although the 95% HPDR for single gene variance of fertility of strain R33 included zero, this criterion for genotypic mean (for 1 locus) suggested a possibility of segregation of a single gene. Szwaczkowski et al. (2006), based on the Bayesian linear model, found a single locus responsible for the percentage of fertilized eggs in the population studied (R33). In the present paper, estimated allele frequencies were 0.308 for fertility and 0.547 for hatchability, with relatively large standard deviations. However, Beaumont et al. (1997) estimated the positive genetic correlations between susceptibilities to the different stages of embryonic death. Is this the same locus with various effects over time? The hypothesis can be verified by detailed molecular segregation and linkage analysis.
The results correspond with estimates of additive polygenic variances for strain R33 (see Tables 4
and 5
). When single locus effects were omitted, the estimates of polygenic variance increased strongly, especially for hatch-ability. In the case of this trait, a similar tendency was also registered for permanent environmental variance estimates. Of course, it considerably influenced heritability estimates. Relatively large participation of major gene variance in total variance compared with polygenic heritability should be stressed. This confirms a segregation of 1 biallelic single locus for hatchability in this strain.
For the second population (N88), the estimates of variance components, and in consequence, the genetic parameters, were balanced for the 3 models used. Contrary to strain R33, it indicates no segregation of single genes determining fertility and hatchability.
Although heritability varies across populations, models, and methods, a number of authors (Foerster, 1993; Brah et al., 1999; Bennewitz et al., 2007) concluded that reproductive traits are lowly heritable. Moreover, the heritability for fertility is usually lower than for hatchability. Despite different measurements of these characters, the results correspond with reports by many authors (Sapp et al., 2005; Bennewitz et al., 2007).
Finally, further analysis based on the molecular study should be suggested to localize the locus responsible for hatchability.
| FOOTNOTES |
|---|
Received for publication November 22, 2007. Accepted for publication February 5, 2008.
| REFERENCES |
|---|
|
|
|---|
Bennewitz, J., O. Morgades, R. Preisinger, G. Thaller, and E. Kalm. 2007. Variance component and breeding value estimation for reproductive traits in laying hens using a Bayesian threshold model. Poult. Sci. 86:823–828.
Brah, G. S., S. S. Jagtar, and M. L. Chaudhary. 1999. Variance and covariance component analysis of incubational mortality in chickens. Arch. Tierz. 42:295–302.
Bumstead, N. 1998. Genomic mapping of resistance to Mareks disease. Avian Pathol. 27:278–281.
Cobanoglu, O., P. J. Berger, and B. W. Kirkpatrick. 2005. Genome screen for twinning rate QTL in four North American Holstein families. Anim. Genet. 36:303–308.[CrossRef][Web of Science][Medline]
Davis, G. H. 2005. Major genes affecting ovulation rate in sheep. Genet. Sel. Evol. 37:11–23.[CrossRef]
Delany, M. E., D. E. Muscarella, and S. E. Bloom. 1994. Effects of rRNA gene copy number and nuclear variation on early development: Inhibition of gastrulation in rDNA-deficient chick embryos. J. Hered. 85:211–217.
Druyan, S., and A. Cahaner. 2007. Segregation among test-cross progeny suggests that two complementary dominant genes explain the difference between ascites-resistant and ascites-susceptible broiler lines. Poult. Sci. 86:2295–2300.
Foerster, A. 1993. Züchterische Möglichkeiten einer Verbesserung der Schlupfrate in Reinzuchtlinien eines Zuchtprogrammes für braune Legehybriden. Institute für Tierzucht und Tierhaltung, Christian-Albrechts Universität, Kiel, Germany.
Gowe, R. S., R. W. Fairfull, I. McMillan, and G. S. Schmidt. 1993. A strategy for maintaining high fertility and hatchability in multiple-trait egg stock selection program. Poult. Sci. 72:1433–1448.[Web of Science]
Guo, S. W., and E. A. Thompson. 1994. Monte Carlo estimation of mixed model for large complex pedigrees. Biometrics 50:417–432.[CrossRef][Web of Science][Medline]
Hocking, P. 2005. Review of QTL mapping results in chickens. Worlds Poult. Sci. J. 61:215–226.[CrossRef][Web of Science]
Janss, L. L. G., R. Thompson, and J. A. M. van Arendonk. 1995. Application of Gibbs sampling for inference in a mixed major gene-polygenic inheritance model in animal populations. Theor. Appl. Genet. 91:1137–1147.[Web of Science]
Jassim, E. W., M. Grossman, W. J. Koops, and R. A. J. Luykx. 1996. Multiphasic analysis of embryonic mortality in chicken. Poult. Sci. 75:464–471.[Web of Science][Medline]
Kadarmideen, H. N., and L. L. G. Janss. 2005. Evidence of a major gene from Bayesian segregation analyses of liability to osteochondral diseases in pigs. Genetics 171:1195–1206.
Lahav, T., G. Atzmon, S. Blum, G. Ben-Ari, S. Weigend, A. Cahaner, U. Lavi, and J. Hillel. 2006. Marker-assisted selection based on multi-trait economic index in chicken: Experimental results and simulation. Anim. Genet. 37:482–488.[CrossRef][Web of Science][Medline]
Liptoi, K., and A. Hidas. 2006. Investigation of possible genetic background of early embryonic mortality in poultry. Worlds Poult. Sci. J. 62:326–337.[CrossRef][Web of Science]
Rao, Y., X. Shen, M. Xia, C. Luo, Q. Nie, D. Zhang, and X. Zhang. 2007. SNP mapping of QTL affecting growth and fatness on chicken GGA1. Genet. Sel. Evol. 39:569–583.[CrossRef][Web of Science][Medline]
Rathje, T. A., G. A. Rohrer, and R. K. Johnson. 1997. Evidence for quantitative trait loci affecting ovulation rate in pigs. J. Anim. Sci. 75:1486–1494.
Sapp, R. L., R. Rekaya, I. Misztal, and T. Wing. 2005. Longitudinal multiple-trait versus cumulative single-trait analysis of male and female fertility and hatchability in chickens. Poult. Sci. 84:1010–1014.
Schmidt, G. S., E. A. P. Figueiredo, A. P. Ledur, and M. C. Fairfull. 1994. Selection for reproductive traits in a white egg stock breeding programme. Pages 17–20 in Proc. 5th World Congr. Genet. Appl. Livest. Prod., Guelph, Canada.
Scott, D. W. 1992. Multivariate density estimation: Theory, practice and visualization. John Wiley and Sons Inc., New York, NY.
Skotarczak, E., K. Moli
ski, T. Szwaczkowski, and A. Dobek. 2007. Bayesian analysis of genetic backgrounds of twinning rate of Thoroughbred horses. J. Anim. Feed Sci. 16:527–538.
Sørensen, D. A., S. Andersen, D. Gianola, and I. Kørsgaard. 1995. Bayesian inference in threshold models using Gibbs sampling. Genet. Sel. Evol. 27:229–249.[CrossRef][Web of Science]
Sørensen, D., and D. Gianola. 2002. Likelihood, Bayesian, and MCMC Methods in Quantitative Genetics. Springer-Verlag New York Inc., New York, NY.
Szwaczkowski, T. 2003. Use of mixed model methodology in poultry breeding: Estimation of genetic parameters. Pages 165–201 in Poultry Genetics, Breeding and Biotechnology. W. M. Muir and S. E. Aggrey, ed. CAB Int., Wallingford, Oxfordshire, UK.
Szwaczkowski, T., E. Skotarczak, A. Dobek, and K. Moli
ski. 2006. Preliminary studies on genetic backgrounds for fertility and hatchability in laying hens. Page 233 in XII Eur. Poult. Congr., Verona, Italy.
Tuiskula-Haavisto, M., M. Honkaturia, J. Vikki, D. J. de Koning, N. F. Schulman, and A. Maki-Tanila. 2002. Mapping of quantitative trait loci affecting quality and production traits in laying hens. Poult. Sci. 81:919–927.
Van Kaam, J. B. C. H. M., M. A. M. Groenen, H. Bovenhuis, A. L. J. Veenendaal, and J. A. M. Van Arendonk. 1999. Whole genome scan in chickens for quantitative trait loci affecting growth and feed efficiency. Poult. Sci. 78:15–23.
W
yk, S. 1978. System SELEKT dla stad zarodowych drobiu. Wyniki Prac Badawczych Zakladu Hodowli Drobiu Instytutu Zootechniki 7:7–22. (in Polish)
Zhou, H., N. Deeb, C. M. Evock-Clover, A. D. Mitchell, C. M. Ashwell, and S. J. Lamont. 2007. Genome-wide linkage analysis to identify chromosomal regions affecting phenotypic traits in the chicken. III. Skeletal integrity. Poult. Sci. 86:255–266.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |