|
|
||||||||
ANIMAL GENETICS |

* Department of Animal Science, University of Nebraska, Lincoln 68581-0908; and
and
Department of Animal Science, North Carolina State University, Raleigh 27695
| Abstract |
|---|
|
|
|---|
= 0.01, 0.05, and 0.10). Thirty-one putative QTL for reproductive traits and two QTL for birth weight were identified (P < 0.10). One Mendelian QTL for FF (P < 0.05), one for NBA (P < 0.05), three for NSB (P < 0.05), three for NN (P < 0.05), seven for AP (P < 0.10), five for MUM (P < 0.10), and one for BWT (P < 0.10) were found. Partial imprinting of QTL affecting OR (P < 0.01), BWT (P < 0.05), and MUM (P < 0.05) was detected. There were four paternally expressed QTL for NN (P < 0.10) and one each for AP (P < 0.05) and MUM (P < 0.10). Maternally expressed QTL affecting NSB (P < 0.10), NN (P < 0.10), and MUM (P < 0.10) were detected. No QTL were detected with half-sib analyses. Multiple QTL models with imprinting effects are more appropriate for analyzing F2 data than single Mendelian QTL line-cross models.
Key Words: Imprinting Pigs Reproduction Quantitative Trait Loci Weight
| Introduction |
|---|
|
|
|---|
Haley et al. (1994)
used regression to estimate additive (a) and dominance (d) coefficients in F2 populations. Based on their model and line origin probabilities, de Koning et al. (2000)
defined contrasts to estimate a, d, and imprinted QTL effects in F2 line-cross data.
Rathje et al. (1997)
and Cassady et al. (2001)
identified 16 putative QTL for reproductive traits with single Mendelian QTL models applied to an F2 population from lines selected for increased litter size. An assumption was that marker genotypes associated with multiple QTL were uncorrelated. If genotypes are correlated, QTL effects will be overestimated and the variance explained by QTL will not be additive, in which case single-QTL models may detect only one or two QTL explaining significant portions of the variation (Schork et al., 1993
). Composite interval mapping (Zeng, 1993
, 1994
) with multiple QTL models was developed to increase the power to detect additional QTL.
The objectives of this study were to reanalyze the F2 data reported by Cassady et al. (2001)
with multiple QTL models applied sequentially to identify Mendelian and imprinted QTL affecting reproduction and early growth in pigs.
| Materials and Methods |
|---|
|
|
|---|
Modeling of QTL
Let the quantitative phenotype, y, be a linear function of a single QTL that influences it:
![]() |
where µ is the overall mean,
is the effect of the QTL, and
is the random environmental deviation. Genotypic-specific means of the QTL effect with two alleles, Q and q, are shown in Table 1
. Because QTL genotypes are not exactly known, probabilities associated with QTL effects are estimated with marker information and used in regression models to identify QTL.
|
|
Data for each trait were fitted first to a single Mendelian QTL, line-cross model. The reduced model included fixed effects of replicate, sire-dam combination (included to adjust for polygenic effects), and covariates of number of fully formed pigs in a litter for BWT and number of pigs weaned and age at weaning for WWT. The full model also included regressions on contrasts of line origin probabilities for additive and dominance effects.
Imprinting effects were tested comparing the following models:
![]() |
![]() |
![]() |
![]() |
where yij is the phenotype of the ith F2 offspring; µ is the combined fixed effects of intercept, replicate, polygenic effect defined as sire-dam combination, covariate according to trait analyzed (the number of fully formed pigs in a litter for birth weight and number of pigs weaned, and age at weaning for weaning weight); p, m, and d are the paternal, maternal, and dominance effects for the imprinted QTL, respectively; cpi is the coefficient of the ith individual for the paternal component at the imprinted QTL; cmi is the coefficient of the ith individual for the maternal component at the imprinted QTL; cdi is the coefficient of the ith individual for the dominance component at the imprinted QTL; and eij is the residual error.
The coefficients cpi and cmi were calculated as described in Table 2
. Calculated LOD scores were used to compare F to N to detect imprinting effects. If both Mendelian and imprinting scans indicated a QTL at a position, the imprinting model was tested against the Mendelian model and, when nonsignificant, the Mendelian model was considered appropriate. To determine the mode of action, LOD scores for F vs. P and F vs. M were calculated. Paternal expression was indicated when only F vs. M was significant, maternal expression was indicated if only F vs. P was significant, and partial imprinting was indicated if both or neither of these contrasts were significant. As QTL scans identify chromosomal regions that may harbor a single gene or multiple genes affecting the trait, and imprinted genes have tendencies to cluster together in the genome, it is possible to identify partially imprinted QTL, indicating that both maternal and paternal imprinted genes in that region affect the trait.
A sequential multiple QTL search with forward selection model building procedures was used to find best-fitting models. From results of single Mendelian QTL models for each trait, the position of the QTL with the highest LOD score that exceeded the critical level (
= 0.10) was chosen as a background effect. An additional Mendelian QTL model was fitted with the position of background QTL fixed in the model. The model was:
![]() |
where yijk is the phenotype of the ith F2 offspring; µ is the combined fixed effects of intercept, replicate, polygenic effect defined as sire-dam combination, and covariate according to trait analyzed (the number of fully formed pigs in a litter for birth weight and number of pigs weaned and age at weaning for weaning weight); al, dl, an, and dn are the additive and dominance effects for QTL l and n, respectively; cail is the coefficient of the ith individual for the additive effect at the lth QTL (selected by the largest significant LOD score); cdil is the coefficient of the ith individual for the dominance effect at the lth QTL; cain is the coefficient of the ith individual for the additive effect at the nth QTL; cdin is the coefficient of the ith individual for the dominance effect at the nth QTL; and eijk is the residual error.
Calculated LOD scores were used to compare the full model with the reduced, single-QTL model. Rounds of genome scans with n 1 identified Mendelian and imprinted QTL as background effects in the reduced model were completed until the largest LOD score corresponding to the nth QTL was less than the genome-wide threshold level (
= 0.10). The effects of QTL incorporated into multiple QTL models as background effects were estimated in the full model. The effects of QTL identified but not used as background effects were estimated in the model in which it was last significant.
Data also were fitted to half-sib QTL models. Reduced models included fixed effects of replicate, sire-dam combination, and appropriate covariates for BWT and WWT. Full models also included covariate coefficients of probability of inheriting the control line allele within an F1 sire family.
Thresholds
Genome-wide significance levels (
= 0.01, 0.05, and 0.10) were estimated from 475 permutations of the data, the number required to estimate
0.05 with a SE of 0.01 for y, the number of times that the test statistic exceeded the critical value (i.e., distributed binomially with parameters N and
). When shuffling data, associations between effects in the reduced model were retained. Thresholds for critical value of
= 0.01, 0.05, and 0.10 were LOD scores that exceeded the 99th, 95th, and 90th percentiles, respectively, of the ranked scores.
| Results |
|---|
|
|
|---|
Progression of genome scans with putative QTL, locations, LOD scores, mode of action, and selection of background loci are reported in Table 3
. Estimates of additive, dominance, and imprinting effects are in Table 4
. Additive effects were estimated as the mean of individuals homozygous for the allele inherited from Line C minus the mean of those homozygous for alleles inherited from Line I. Mean genotype specific deviations from the average of the two homozygotes were calculated and are also in Table 4
. Deviations for Mendelian QTL are a, d, d, and a for the CC, CI, IC, and II genotypes, respectively. Mean deviations for imprinted QTL are (p + m), (d + p m), (d p + m), and (p m) for the CC, CI, IC, and II genotypes.
|
|
Mendelian QTL in SSC11 at 52 and 71 cM affected FF and NBA, respectively (Table 3
). Neither of these positions significantly affected NSB, which is the difference between FF and NBA; however, Mendelian QTL on SSC13 (P < 0.05), SSC5 (P < 0.10), and SSC12 (P < 0.10) affecting NSB were identified in single QTL models (Table 3
). After fitting the QTL in SSC13 in the model, the one on SSC5 remained (P < 0.10), but the one on SSC12 was not significant after fitting both the SSC13 and SSC5 positions. Dominance expression of the SSC13 QTL (a = 0.39 ± 0.12, d = 0.53 ± 0.20), overdominance expression of the one in SSC5 (a = 0.06 ± 0.17, d = 1.00 ± 0.29), and additive expression of the one in SSC12 (a = 0.38 ± 0.13, d = 0.40 ± 0.23) were observed (Table 4
).
One Mendelian QTL affecting MUM on SSC12 at 98 cM (P < 0.10) was identified in a single QTL model (Table 3
). The four additional Mendelian QTL identified in multiple-QTL models were on SSC12 at 70 cM (P < 0.05) that entered the model after fitting the QTL at 98 cM, on SSC6 at 64 cM (P < 0.10), which became significant after fitting the two QTL on SSC12, and one each on SSC2 (P < 0.10) and SSC6 at 165 cM (P < 0.10) that were identified in four-QTL models. With the exceptions of the SSC2 QTL, for which the estimate of d was 0.02 ± 0.13, and the SSC6 QTL at 64 cM, for which the estimate of a was 0.06 ± 0.10, these QTL tended to be dominant with absolute values of a and d ranging from 0.16 to 0.45 (Table 4
).
The Mendelian QTL for NN with the greatest LOD score in a single QTL model (Table 3
) was on SSC11 (P < 0.05). The second and third QTL to be added were in SSC8 and SSC7, respectively, (P < 0.05). The LOD score for the SSC7 QTL increased from 2.30 (P < 0.10) in a single QTL model to 2.90 (P < 0.05) after fitting the positions of QTL on SSC11 and SSSC8. Overdominance was found for the SSC11 QTL (a = 0.05 ± 0.10, d = 0.72 ± 0.17), whereas the magnitude of estimates of a and d were similar for the other two QTL, indicating dominance (Table 4
).
Three of the QTL for AP identified in single-QTL models (Table 3
) were on SSC8 at 101 (P < 0.10), 136 (P < 0.10), and 172 cM (P < 0.05); two were on SSC7 at 1 (P < 0.05) and 58 cM (P < 0.10); and one was on SSC12 (P < 0.10). After fitting the position on SSC8 at 172 cM, only the two QTL on SSC7 remained significant; however, a fourth scan in which these three positions were fixed identified another Mendelian QTL on SSC18 (P < 0.10). These QTL tended to be dominant as estimates of d were similar to or greater than estimates of a in each case (Table 4
).
There was evidence for a dominant Mendelian QTL affecting BWT (a = 5.32 ± 1.92 g, d = 8.62 ± 3.65 g, Table 4
) on SSC12 (P < 0.10). No Mendelian QTL for WWT or OR were found.
Imprinted QTL
Six imprinted QTL were identified in single, imprinted-QTL models, one each for OR, NSB, MUM, NN, AP, and BWT (Table 3
). Those for OR (SSC9, P < 0.01), MUM (SSC6 at 81 cM, P < 0.05), and BWT (SSC6, P < 0.05) were partially imprinted (Table 4
). The QTL on SSC14 for NSB was maternally expressed (P < 0.10) and those for AP (SSC15 at 98 cM, P < 0.05) and NN (SSC15 at 109 cM, P < 0.05) were paternally expressed.
Six more imprinted QTL were found in multiple-QTL models (Table 3
). Using three QTL as background effects, a maternally expressed QTL affecting MUM on SSC2 at 6 cM (P < 0.10) and a paternally expressed QTL on SSC6 at 191 cM (P < 0.10) were found. Using two QTL as background effects, there was evidence for four more imprinted QTL affecting NN. These included paternally expressed QTL on SSC6 at 85 cM (P < 0.05) and SSC15 at 64 cM (P < 0.10) and QTL on SSC1 (P < 0.10) and SSC6 (171 cM in single-QTL model, 163 cM in four-QTL model, P < 0.10) that were classified as Mendelian in single- and two-QTL models but were identified as imprinted in three- and four-QTL models.
Half-Sib Analyses
No QTL were identified with half-sib analyses.
| Discussion |
|---|
|
|
|---|
Although the power to identify QTL in segregating populations is decreased due to the combination of differences in gene frequencies and size of genetic effects (de Koning et al., 2002
), the model-building procedures resulted in 15 additional Mendelian or imprinted QTL that were not detected by Cassady et al. (2001)
. With the exception of the QTL on SSC12 affecting MUM at 98 cM and NSB at 37 cM, results of single QTL scans were identical to those of Cassady et al. (2001)
. The discrepancies on SSC12 were due to an error in the marker data in the earlier analysis resulting in mis-identification of QTL regions. By using imprinted QTL models or by fixing background QTL to decrease residual variation, two additional QTL for NSB, three for NN, two for AP, and seven for MUM were detected.
The power of multiple-QTL models compared with a single-QTL model is best illustrated for MUM. A QTL in SSC12 close to marker SWR10 at 98 cM was located with a single QTL model; a nonsignificant local peak existed in the 60- to 75-cM interval near marker S0090 (Figure 1
). A second model in which the QTL near SWR10 was fixed identified another QTL on SSC12 within the interval of the local peak at 70 cM. Linkage between these positions was expected to create a correlation between them such that the QTL at 70 cM may not be detected after fixing the one at 98 cM. Apparently, adequate recombination had occurred so that the variation explained by the QTL near marker S0090 was significant. By fixing two QTL in models for MUM, a Mendelian QTL on SSC6 (64 cM) and an imprinted QTL on SSC6 (81 cM) were identified. However, the QTL at 64 cM was not significant in a four-QTL model in which the imprinted QTL on SSC6 was included as the third background position, probably due to correlation between these positions caused by linkage. Four additional QTL for MUM, two on SSC2 at 6 and 29 cM and two on SSC6 at 165 and 191 cM, were identified in the four-QTL model.
|
Mendelian QTL for MUM on SSC12, SSC2, and SSC6 were found. Fixing the QTL on SSC12 (98 cM, a = 0.32 ± 0.09) and SSC2 (6 cM, a = 0.27 ± 0.08) for Line C alleles was estimated to cause a decrease of 1.18 mummified pigs at birth compared with fixing Line I alleles. However, fixing the Line I alleles for the QTL on SSC12 (70 cM, a = 0.16 ± 0.09) and SSC6 (165 cM, a = 0.21 ± 0.09) is expected to decrease MUM by 0.74 compared with Line C alleles. Dominance was found for the Mendelian QTL in SSC6 at 64 cM (a = 0.06 ± 0.10, d = 0.77 ± 0.23). Individuals heterozygous for Line I and Line C alleles were estimated to have 0.77 fewer mummified piglets than the average of those homozygous for Line I and Line C alleles. Partial imprinting occurred for the QTL on SSC6 (81 cM). Identification of partial imprinting indicated that there may be at least two imprinted QTL, one paternally and one maternally expressed. The joint effect of inheriting a Line C allele from the sire and a Line I allele from the dam (CI genotype) was estimated to result in the least mummified pigs of any genotype at this QTL position, 0.58 mummified pigs less than the reciprocal genotype (Table 4
.). Fixing either Line C or Line I alleles is expected to increase the number of mummified pigs by approximately 0.8 compared with the CI genotype. There was complete paternal expression of the QTL on SSC6 (191 cM) affecting MUM. Fixing Line I alleles was estimated to decrease mummified pigs by at least 0.28 pigs compared with other genotypic means. The CI genotype was estimated to have the largest number of mummified pigs.
Similar application of estimates of QTL effects for other traits produces the mean genotypic differences in Table 4
. Interpretations are straightforward for Mendelian QTL, but are less so for imprinted QTL. For OR and BWT, the joint effect of inheriting a Line I allele from the sire and a Line C allele from the dam was estimated to maximize response compared with other genotypes at these QTL. However, in the case of the OR QTL, the dominance effect exceeded other effects and resulted in both heterozygotes with greater OR than either homozygote. In contrast, the dominance effect for the BWT QTL was small compared with other effects; thus, the IC heterozygote was predicted to have the greatest BWT and the CI heterozygote the least.
When a QTL was imprinted and the magnitude of the dominance effect was at least twice as large as the imprinting effect, as found for maternally expressed QTL for NN, both heterozygotes had means outside the range of the two homozygotes. However, imprinted QTL for which d was small, such as paternally expressed QTL for NN and maternally expressed QTL for MUM, resulted in positive effects in one heterozygote and one homozygote and negative effects in the other heterozygote and homozygote. When the estimated effect was negative, means for CC and CI (IC) genotypes were negative for paternally (maternally) expressed QTL and means for the other genotypes were positive. The four paternally expressed QTL for NN followed this pattern, suggesting that the Line C allele inherited from the sire would decrease the number of nipples. Maternally expressed QTL for MUM indicated a decrease in mummified pigs when the Line C alleles were inherited from the dam.
Composite interval mapping was designed to remove the residual error associated with known markers or QTL affecting the trait to detect additional QTL explaining less variation. The increase in power was readily apparent in scans for QTL for MUM (Figure 1
) and AP (Figure 2
). Some QTL identified in single-QTL models were not significant in full models as illustrated for AP in Figure 3
. For AP, six QTL were identified in the single Mendelian QTL scan; however, three of these QTL were not significant in subsequent scans. If several QTL are acting independently, then single- and multiple-QTL models are expected to identify the same QTL. But if epistasis exists, effects of multiple QTL are not independent (Schork et al., 1993
), and QTL identified in single QTL scans may not be significant in sequential scans in which the one explaining the most variation is fixed. Another possibility is that only one QTL exists, but there is a covariance between marker genotypes at the two positions. Depending on the magnitude of the covariance, a single scan might identify two QTL, but after fixing the one explaining the greatest variation, a sequential scan may not identify additional QTL. Covariation between marker genotypes could be because chromosomal associations that occurred in base generation animals of the selection lines were maintained during the selection experiment, or be due to random sampling.
|
|
Other imprinted effects reported in swine include a QTL on SSC4 affecting abdominal fat (Knott et al., 1998
) and regions on SSC6 affecting i.m. fat and on SSC7 affecting muscle depth (de Koning et al., 2000
). In humans, imprinted genes can be found on every chromosome and tend to cluster together (Tycko and Morison, 2002
; Okita et al., 2003
). Five of the 12 imprinted QTL identified in this study were on SSC6 and three were on SSC15.
Half-sib analyses have the potential to identify QTL when F1 sires differ in QTL genotype. This can occur when founder lines are segregating, even if they do not differ in allele frequencies. Therefore, half-sib analyses may identify QTL not identified in line-cross analyses. However, if founder populations are completely inbred, a half-sib analysis requires four times the number of individuals to obtain the same power of the F2 line-cross analysis (Weller, 2001
). The power of half-sib analyses also increases with increasing family size. The LOD peaks in half-sib analyses did not reach the critical threshold levels and no QTL were indicated. However, LOD peaks for NSB and NN coincided with positions of peaks in the line-cross analyses. Family sizes in this experiment were insufficient to detect QTL with the half-sib analyses.
| Implications |
|---|
|
|
|---|
| Footnotes |
|---|
2 Correspondencephone: 402-472-6404; fax: 402-472-6362; e-mail: rjohnson5{at}unl.edu.
Received for publication May 10, 2004. Accepted for publication September 2, 2004.
| Literature Cited |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
S. Sato, K. Atsuji, N. Saito, M. Okitsu, S. Sato, A. Komatsuda, T. Mitsuhashi, K. Nirasawa, T. Hayashi, Y. Sugimoto, et al. Identification of quantitative trait loci affecting corpora lutea and number of teats in a Meishan x Duroc F2 resource population J Anim Sci, November 1, 2006; 84(11): 2895 - 2901. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |