|
|
||||||||


,1
* CERSA-Fondazione Parco Tecnologico Padano and
and
Istituto Biologia e Biotecnologia Agraria, Consiglio Nazionale delle Ricerche Segrate, Italy 20090; and
and
Department of Animal Science, University of Tennessee, Knoxville 37996-4574
| Abstract |
|---|
|
|
|---|
Key Words: Gametes Imprinting Litter Size Pigs
| Introduction |
|---|
|
|
|---|
The second approach has used pedigree and phenotypic data to estimate variance components indicative of imprinting. Three different statistical models have been proposed. The first model (Schaeffer et al., 1989
) involves estimating effects of gametes from a given parental origin. Gametic effects are correlated according to the gametic relationship matrix (Smith and Allaire 1985
). This model is the most appropriate method for estimating effects of imprinting (Essl and Voith, 2002b
), but it has not been routinely applied due to a lack of software. Tier and Sölkner (1993)
treated gametes as homozygous diploid individuals and included them in a standard relationship matrix. Twice the variance of these effects was considered equivalent to the gametic variance of Schaeffer et al. (1989)
. This approach has been applied to swine (De Vries et al., 1994
) and beef cattle (Engellandt and Tier, 2002
). Essl and Voith (2002a)
proposed a third test based on the difference between variances of maternal and paternal gametic effects obtained from dam and sire models, respectively. Essl and Voith (2002a)
applied both of these latter two methods to data from dairy cattle.
Studies of imprinting in livestock have generally been applied to production traits. According to certain theories about the evolutionary basis of imprinting (Reik and Walter, 2001
), reproductive traits could also be affected. The objectives of this study were to develop software for and to apply the model proposed by Schaeffer et al. (1989)
to test for effects of imprinting on litter size in swine.
| Materials and Methods |
|---|
|
|
|---|
Methods
The models for evaluation of gametic effects included effects of contemporary group, parity, mate, and permanent environment, in addition to the genetic effects. Contemporary groups were defined on two levels, herds and herd-year-seasons (HYS) of farrowing. Data were from 735 Landrace and 2,455 Yorkshire herds. The HYS were of variable length, ranging from 2 mo to 1 yr, depending on the number of animals per HYS. Each HYS contained at least two contemporaries. The HYS were formed in a stepwise manner. Initially, HYS of 2-mo seasons were created. All HYS with at least two animals were kept. The HYS with a single animal were then merged with adjacent HYS, forming HYS four months in length. This process was repeated in 2-mo steps until reaching a maximum length of one year. Animals with no herdmates within the same year were not included in the analysis. This approach resulted in the definition of 8,440 and 22,702 HYS for Landrace and Yorkshires, respectively. Parity numbers ranged from 1 to 17 (and unknown) in the original datasets and these were defined into six different groups (1st, 2nd, 3rd to 5th, 6th to 9th,
10th, and unknown), based on average litter sizes and number of animals per parity. The effect of mate was included to account for effects of differences in fertility of service sires on litter size. The numbers of mates (service sires) were 5,289 and 14,193 for Landrace and Yorkshires, respectively.
The model for the estimation of gametic effects was defined by the equation:
![]() | [1] |
where yijklmn is the observed value for litter size of the oth record; hi is the fixed effect of the ith herd, pj is the fixed effect of the jth parity group, hysk is the random effect of the kth HYS, ml is the random effect of the lth mate (i.e., service sire), pem is the random effect of the permanent environment of dam m, am is the random additive genetic effect of animal m, gn is the random effect of the nth gamete (of either maternal or paternal origin for animal m), and eijklmno is the random residual for record o.
This model assumes that additive genetic effects are distributed MVN(0,A
2a) and gametic effects are distributed MVN(0,G
2g), where A and G are additive genetic and gametic relationship matrices (Smith and Allaire, 1985
), respectively, and that variability in gametic effects is observed only in maternally or paternally transmitted gametes, not both. The model also assumes that other potential genetic influences such as effects of the fetuses, maternal effects, or cytoplasmic effects are negligible.
Variance components were estimated by a Bayesian procedure that employed Gibbs sampling to obtain posterior distributions of the various parameters (Sorensen, 1999
). Prior distributions were uniform, normal, and inverted
2 for fixed effects, random effects, and variances, respectively. Priors for variances were informative but received very little weight. Fixed and random effects were sampled from normal distributions and variance components were sampled from inverted
2 distributions. For each analysis, three chains of 100,000 cycles were generated (i.e., 300,000 in total). Based on visual inspection, the first 5,000 samples from each chain were discarded as a burn-in period to ensure sampling from the appropriate marginal distributions. Posterior means and SD were estimated from the respective posterior distributions based on saving every twentieth sample. This sampling strategy was adopted to ensure that all posterior means were based on at least 75 independent samples.
Simulated Data
To aid in interpretation of results based on real data, similar analyses were also applied to sets of simulated data. Similar statistical analyses had been previously applied to other sets of real and simulated data (our unpublished observations) and two observations had been made. First, some population structures precluded the detection of imprinting effects. For example, when the model defined by Eq. [1]
was previously applied to carcass data in swine (for which animals with data had no offspring), no significant gametic variance was observed, even when an effect of imprinting was simulated and added to each of the phenotypic records. In that analysis, gametic variance was partitioned into additive genetic variance. Second, in previous analyses of simulated data, estimates of gametic variance were small (<1.0% of total variance), but non-zero in situations where no imprinting effects were simulated. Therefore, the goals of the analyses of simulated data were 1) to determine if the population structure for the litter size data was sufficient for the detection of imprinting effects (i.e., to help avoid type-II errors) and 2) to observe estimates of gametic variance in situations for which no imprinting effects existed (i.e., to help avoid type-I errors).
The simulated data was obtained using the Landrace data as a template. Data were simulated according to Eq. [1]
. The mean of simulated litter size was set at 10.0 and all effects were simulated by random sampling from normal distributions. Standard deviations for effects of parity, herd, HYS, mate, permanent environment, and residual were 0.5, 0.6, 0.13, 0.16, 0.34, and 5.65, respectively, and were based on a preliminary analysis of the data using REML. The underlying variance of single gametic effects was either 0.0 (i.e., no imprinting) or 0.15. The value of 0.15 was arbitrarily chosen to equal one-third of additive genetic variance. Two gametes were simulated for each animal, but only one of the two contributed to phenotypes, because only maternal or paternal imprinting was simulated in any given model. Additive genetic variance was either 0.54 (no imprinting) or 0.45. An infinitesimal model was assumed for both gametic and additive genetic effects and inbreeding of parents was accounted for in Mendelian sampling of genetic effects.
Five sets of data were simulated with no imprinting effect and these data sets were each evaluated twice, with models that included 1) maternal imprinting or 2) paternal imprinting. A sampling chain of 50,000 rounds was generated for each of these analyses. One set of data each was generated with real effects of paternal imprinting and maternal imprinting. Each of these datasets was analyzed twice, once with effects of gametes of paternal origin in the model (maternal imprinting) and once with maternal gametic effects (paternal imprinting). Each analysis included 100,000 rounds of Gibbs sampling. Prior values used for gametic and additive genetic variances were 0.10 and 0.50, respectively, for all models. The true values were used as priors for other sources of variance. Posterior means were calculated for the various parameters.
| Results |
|---|
|
|
|---|
|
|
Two primary conclusions were drawn from the analyses of simulated data. The first was that (at least for these data) a low but non-zero estimate of gametic variance from this type of analysis cannot be interpreted as an indication of a small but real influence of gametic imprinting. Southwood et al. (1989)
and Boettcher et al. (1997)
reported similar results when applying models that included cytoplasmic effects to data for which no such effects were simulated. In those studies, the estimates of cytoplasmic variance were small but non-zero, even in situations where data were simulated without a true cytoplasmic effect, apparently due to partitioning of a small proportion of the additive genetic effect into the cytoplasmic genetic effect. The second conclusion was that the structure of data analyzed in these studies was sufficient to detect effects of gametic imprinting when they existed. Estimates of all variance components were correct when effects of expression of gametes of either paternal or maternal origin were simulated and the correct model was applied. As mentioned previously, some data structures (data from only a single generation, for example) are not sufficient for detecting such effects by applying similar approaches for analysis (our unpublished observations).
Real Data
Landrace Data.
Posterior means and SD of the various components of variance when analyses that assumed variability in effects of certain maternally or paternally transmitted genes were applied to the Landrace data are shown in Table 3
. The first half of the table shows the estimates of variances from the model in which expression of genes of paternal origin (i.e., maternal imprinting) was assumed. None of the individual factors in the model accounted for as much as 10% of the total variance and residual effects accounted for >80% of the variability in litter size. Heritability for litter size was 0.078, which is at the lower end the range of estimates previously reported in the literature (e.g., Lamberson, 1990
; Roehe and Kennedy, 1995
; Chen et al., 2003
). The posterior estimate of gametic variance was different from zero since the mean was 0.053 and the lowest sampled value was 0.018, but these values were very similar to estimates obtained from simulated data for which no effect of imprinting was simulated (Table 1
). These results suggest that maternally imprinted loci have no detectable effect on litter size in the Landrace breed. At best, such effects could explain <1% of the variance, but such an effect is too small to be detected unambiguously by the analysis undertaken.
|
The posterior mean of additive genetic variance from the model with maternal expression (0.471) was less than the estimate with paternal expression (0.539). Consequently, the heritability from the maternal expression model was 1% lower (0.068 vs. 0.078) than from the model with maternal imprinting. The difference in the additive genetic variance between the two analyses was basically equal in magnitude to the difference in gametic variances between the two models. The residual variance and its proportion of total variance were essentially the same from the two analyses.
Yorkshire Data.
Results from analyses of the Yorkshire data are in Table 4
. Estimates from both of the models (for expression of genes from gametes of either paternal or maternal origin) were very similar to each other. Specifically, both the posterior means and SD of variances for HYS and mate effects were equal up to three decimal places. Posterior means for all other sources of variance differed only at the third decimal place (i.e., would round to the same values at the second decimal place). Of particular interest in this study were the estimates of the gametic variances. Posterior means were 0.060 and 0.057 when expression of genes of paternal and maternal origins, respectively, was assumed. These estimates were only somewhat greater than the overall means across replicates from data simulated without an imprinting effect (Table 1
), and were lower than some of the posterior means from the five individual replicates. This result suggests that the true value of the gametic effect is not large in the Yorkshire population from which these data were recorded.
|
| Discussion |
|---|
|
|
|---|
Another interesting aspect of the results was that a significant estimate of gametic variance was observed only for the Landrace breed. To ensure that this result was not due to differences in data structure between the populations, we resimulated data with variable effects of gametes of maternal origin using the Yorkshire population as a template and estimated variances using a model that included these effects. The resulting estimates for gametic variance were correct, indicating the differences observed were not due to data structure.
Given that the differences between Landrace and Yorkshire were not due to differences in the structure of the data, one may be tempted to conclude that no genes associated with litter size are imprinted in the Yorkshire breed. However, nonsignificant gametic variances from this analysis are not necessarily indicative that gametic imprinting does not occur within the population tested. Rather, a significant estimate of gametic variance for a given parental origin suggests that genetic variability exists in imprinted genes that affect litter size. Only a very small proportion of genes is expected to be imprinted (Barlow, 1995
). Perhaps an equally small number of imprinted genes control litter size in both breeds, but polymorphism or variation in the effect of only one (or very few) of these genes is present only in the Landrace breed. Alternative alleles of the estrogen receptor allele have been shown to have different (unimprinted) effects on litter size in certain breeds and lines of swine (Rothschild et al., 1996
). However, in other populations, this gene can either be nonpolymorphic or have no observable difference between alternative alleles (Drogemuller et al., 2001
).
Assuming the variability observed is associated with polymorphism at a single biallelic locus, the substitution affect of the superior allele could be substantial. If alternative alleles each had frequencies of 50%, the gametic variance of 0.121 would correspond to a difference of approximately 0.70 piglets per litter. This substitution effect would increase for greater disparities in gene frequencies.
The analysis applied to the data was based on the assumption of a normal distribution of gametic effects rather than a biallelic locus. Therefore, we used simulation to confirm the ability of the previously described model (Eq. [1]
) to detect gametic variability associated with a single, biallelic loci. We used the Landrace data structure to simulate litter size controlled in part by an imprinted biallelic locus. We simulated two alleles, each with frequencies of 0.5 and with effects differing by 0.75 piglets. A single replication was generated each for a paternally and maternally imprinted gametic effect. Estimates of variances were 0.152 (paternal imprinting) and 0.138 (maternal), neither significantly different from the expected value of 0.141 (posterior SD were approximately 0.03 in each case).
Although variability in paternally imprinted genes that affect litter size may be the most plausible reason for the significant gametic variance observed in the Landrace breed, other possible explanations exist. As was noted earlier in the explanation of the analysis, the model used assumed that other genetic effects, such as fetal, maternal, or cytoplasmic effects, had no influence on litter size. If such factors were important but ignored, they could cause a bias in the estimate of gametic variance. The genotypes of the pigs of the litters were not considered in this study. In theory, differences in genotypes of embryos could cause variability in their survival. However, Gama et al. (1991)
estimated the direct genetic effects of embryonic survival and reported a heritability of <1%. Maternal sisters share the same maternal effect and cytoplasmic effects and would be expected to have correlated effects due to maternally expressed genes. Thus, these factors could affect estimates of the influence of paternally imprinted genes. However, the importance of these factors is questionable. Previous studies have estimated variance for maternal effects on litter size and the influence of this factor has generally been quite small, usually <1% of the total variance in most situations (e.g., Roehe and Kennedy, 1995
; Chen et al., 2003
). With respect to cytoplasmic effects, Robison (1998)
reported differences of approximately 0.75 piglets in mean litter size between lines formed by reciprocal crossing of Yorkshires and Duroc breeds. He attributed the difference to cytoplasmic effects. However, the study of Robison (1998)
dealt with differences in cytoplasmic effects between breeds rather than variability within a given breed. No other evidence for cytoplasmic effects on litter size within a breed has been published. If the results of this study were due to cytoplasmic effects rather than imprinted effects, such a result would be biologically interesting in its own right. Molecular genetic approaches may be needed to precisely define the biological cause for the results observed in this study. The identification of the gene or genes associated with the observed effect could produce information useful in marker assisted selection programs.
In theory, accounting for effect of imprinting (or another non-Mendelian factor) in national genetic evaluations could improve accuracy of EBV. However, such a procedure would add an additional level of complexity to the model. A gametic model such as the one applied in this study includes an additional number of factors equal to twice the number of animals. Given the relatively small influence of this factor on phenotypic variability (less than 2% of the total variance), and the fact that it seems to be maternally transmitted, accuracies of EBV of boars from current genetic evaluations that ignore imprinting may not be improved greatly by adapting the currently used genetic models to account for gametic effects.
| Implications |
|---|
|
|
|---|
1 Correspondence: Palazzo LITA, Via Fratelli Cervi 93, 20090 Segrate (phone: +39 0221013508; fax: +39 0226412135; E-mail: boettch{at}ibba.cnr.it).
Received for publication February 25, 2003. Accepted for publication May 21, 2003.
| Literature Cited |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. M. Patten and D. Haig Reciprocally Imprinted Genes and the Response to Selection on One Sex Genetics, July 1, 2008; 179(3): 1389 - 1394. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |