|
|
||||||||
1 Animal Breeding and Genetics Group, Wageningen Institute of Animal Sciences, Wageningen University, The Netherlands
2 Holland Genetics, CR Delta, 6802 EB Arnhem, The Netherlands
3 Biometris, 6700 AC Wageningen, The Netherlands
Corresponding author: C. Schrooten; e-mail: chris.schrooten{at}wur.nl.
| ABSTRACT |
|---|
|
|
|---|
= 0.01, chromosomewise). Fifteen MQR were found on Bos taurus autosome (BTA) 14. Four or more MQR were found on BTA 6, 13, 19, 22, 23, and 25. Eight MQR involving udder conformation and udder health and 4 MQR involving production traits and udder health were found. Five MQR were identified for combinations of fertility and udder conformation traits, and another 5 MQR were identified for combinations of fertility and production traits. For 22 MQR, the difference between the correlation attributable to the MQR and the overall genetic correlation was >0.60. Although the false discovery rate was relatively high (0.52), it was considered important to present these results to assess potential consequences of using these MQR for marker-assisted selection.
Key Words: genome scan multiple trait quantitative trait locus dairy cattle
Abbreviation key: BTA = Bos taurus autosome, MAS = marker-assisted selection, MQR = multiple trait quantitative trait region.
| INTRODUCTION |
|---|
|
|
|---|
Selection for a certain trait can lead to genetic changes in other traits because of the genetic correlations between the traits. One of the possible causes for genetic correlation is pleiotropy. Pleiotropic genes influence 2 or more traits. Closely linked QTL can also contribute to genetic correlation because of linkage disequilibrium. Pleiotropic effects of QTL, or closely linked QTL, each affecting a different trait, can affect the value of individual QTL for MAS.
The contribution of identified QTL to the overall genetic correlation can be determined by multiple trait QTL analysis. Multiple trait methods have been developed (e.g., Jiang and Zeng, 1995; Korol et al., 1995; Knott and Haley, 2000; Korol et al., 2001). Weller et al. (1996) developed a method to address the effect of QTL on multiple traits by canonical trait analysis. Schrooten and Bovenhuis (2002) showed that analysis of the covariance between contrasts from single trait analyses for different traits can reveal chromosomal regions affecting multiple traits (multiple trait quantitative trait regions; MQR). They developed a method that utilizes results from single trait QTL analyses to identify MQR in a granddaughter design. The method, however, cannot distinguish between pleiotropic QTL or closely linked QTL.
In this project, the method described by Schrooten and Bovenhuis (2002) is applied to data from a granddaughter design. The goal is to identify chromosomal regions affecting multiple traits (MQR) in dairy cattle.
| MATERIALS AND METHODS |
|---|
|
|
|---|
![]() | ([1]) |
where yij = deregressed breeding value of son j of grandsire i, µ = overall mean, gsi = fixed effect of grandsire i, bik = regression coefficient for grandsire i at location k on the chromosome, Xijk = probability that son j receives a chromosomal segment from grandsire i at position k, and eijk = random residual.
Significant QTL in the single trait analyses were determined by comparing the test statistic at each location with the 10% genomewise threshold, as described by Churchill and Doerge (1994). In this procedure, phenotypes were shuffled within families and reassigned to the sons, thus breaking the association between phenotypes and genotypes. Significance thresholds were obtained from analysis of 10,000 permuted data sets. A QTL was considered suggestive when the threshold would yield one false positive result in a genome scan (Churchill and Doerge, 1994). This threshold corresponds to a chromosomewise type I error (
) of 0.0345.
To detect chromosomal regions affecting multiple traits, the covariance between estimated contrasts for pairs of traits, i.e., the regression coefficient in Equation (1), was determined at each location separately. As shown by Schrooten and Bovenhuis (2002), the expected covariance can be written as
![]() | ([2]) |
where
and
= contrasts for trait v and w for grandsire i; rA = genetic correlation between trait v and trait w; rE = environmental correlation between trait v and trait w;
and
= heritabilities of trait v and trait w, respectively; ns = number of sires per grandsire; ng = number of daughters per sire;
pv and
pw = phenotypic standard deviations of trait v and trait w, respectively; phet = fraction that is heterozygous for the QTL; and av and aw = one-half the difference between the 2 homozygous genotypes (Falconer and Mackay, 1996), i.e., the allele substitution effect when there is no dominance. Note that parameter bik in Equation (1) is an estimate of "a" in Equation (2). Standard deviations, heritabilities, and correlations do not include the effect of the QTL, i.e., these parameters are corrected for the effect of the QTL. A more detailed description of the method and discussion of its properties can be found in Schrooten and Bovenhuis (2002).
When there is no QTL, or the QTL affects only one trait, Equation (2) shows that the covariance does not depend on the values of av and aw, but solely on the polygenic and environmental parameters. Deviations from this expected covariance indicate the presence of an MQR, i.e., a pleiotropic QTL, or closely linked QTL, each affecting one of the traits. It is not possible to distinguish between QTL with a pleiotropic effect and closely linked QTL in a statistical analysis. In this paper, we will therefore consider chromosomal regions affecting multiple traits (MQR), because pleiotropic QTL or closely linked QTL may be involved.
In the analysis, the expected covariance under the null hypothesis, i.e., a situation without a pleiotropic QTL or closely linked QTL, can be determined by calculating the average covariance in 10,000 permuted data sets. Permuted data sets were obtained as described for the single trait analysis. Here, permutation was also used to construct a 99% confidence interval for the average or expected covariance. Covariances outside the confidence interval indicate that this chromosomal region is affecting both traits. Effects can have the same direction of change, i.e., the value of both traits is either increased or decreased, or have opposite direction of change. This can, in some cases, result in a large difference between the correlation because of the MQR and the overall genetic correlation. To detect chromosomal regions where this is the case, trait combinations with large differences between the MQR correlation and the overall genetic correlation were identified. Overall genetic correlations were not based on the data set used in the current project, but on data used for national evaluations (A. Harbers, 2003, personal communication).
| RESULTS |
|---|
|
|
|---|
|
) that was used to construct the confidence interval for the covariance.
|
|
Chromosomes with
4 MQR are presented graphically in Figure 2
and are listed in Table 2
. Results for these chromosomes will be discussed in more detail. Results for other chromosomes are only presented in Table 2
.
|
BTA 13.
On BTA 13, 4 MQR involving production (kg milk and protein and fat percentage) and udder conformation traits (fore udder attachment and udder depth) were located at 0 cM. In the single trait analysis, significant QTL for fore udder attachment and udder depth were found in this region. Overall genetic correlations were close to zero for all 4 trait combinations. The MQR correlations were positive for combinations of udder conformation traits and percentage traits (0.57 to 0.75). The correlation attributable to the MQR affecting fore udder attachment and kilograms of milk was negative (0.63).
BTA 14.
Bos taurus autosome 14 contained 15 MQR. Covariances exceeded the confidence interval in the chromosomal region between 0 and 17 cM, with most extreme covariances between 7 and 13 cM (Table 2
; Figure 2
). In Figure 2
, only this region of BTA 14 is shown. Thirteen of the 15 MQR involved 1 or 2 production traits, and 6 of these involved 2 production traits. Traits other than production that were affected by these MQR were fore udder attachment, front teat placement, udder depth, and milking speed. Most overall genetic correlations between the traits involved were close to zero, except for the correlation between kilograms of milk and kilograms of fat (0.37), among percentage traits (0.78) and between milk yield and percentage traits (0.59 and 0.65). Correlations attributable to the MQR were positive between udder conformation and percentage traits (0.46 to 0.55) and negative between udder conformation traits and kilograms of milk (0.61 to 0.63). The MQR correlations among yield traits were comparable with overall genetic correlations, with the exception of the MQR correlation between kilograms of milk and kilograms of fat (0.24) and kilograms of protein with protein and fat percentages (0.62 to 0.77). The MQR correlations with milking speed were positive for udder conformation traits and fat percentage (0.44 to 0.56) and negative for milk yield (0.46). In the single trait analysis, significant QTL were found for fat percentage, kilograms of milk, and kilograms of protein, whereas suggestive QTL were identified for protein percentage, kilograms of fat, and milking speed in the same region as the MQR.
BTA 19.
Five of the 6 MQR on this chromosome are located at 0 cM and affect one of the udder conformation traits and somatic cell score or 2 udder conformation traits. The sixth MQR was at 39 cM. The correlation attributable to this MQR affecting the 2 fertility traits, interval calvingfirst insemination and nonreturn of daughters at 56 d after insemination, was moderately negative (0.51), whereas the overall genetic correlation was moderately positive (0.36). The MQR correlation between rear udder height and front teat placement was strongly positive (0.82), whereas the overall genetic correlation between these traits was close to zero. For udder depth, fore udder attachment, and front teat placement, the most probable location for QTL, from the single trait analysis, was at 35, 66, and 66 cM, respectively.
BTA 22.
On BTA 22, MQR were found at locations 0 (nonreturn of daughters at 56 d after insemination with both kilograms of fat and kilograms of protein), 18 (front teat placement and protein percentage), and 85 cM (fore udder attachment and nonreturn at 56 d after insemination). The single trait analysis only revealed a suggestive QTL at 0 cM for front teat placement. For the 2 MQR at 0 cM, there was a large difference between the correlation attributable to the MQR (0.58 to 0.64) and the overall genetic correlation (0.19 to 0.28).
BTA 23.
The MQR on BTA 23 affected production and udder conformation traits. These MQR were found in the region between 8 and 14 cM. In the single trait analysis, suggestive QTL for fore udder attachment and kilograms of protein were found in the same region.
BTA 25.
On BTA 25, 5 MQR were located at 80 cM. Four of these MQR involved central ligament. Positive effects on central ligament, caused by the MQR, were associated with positive effects for yield traits, somatic cell score, and rear udder height. The overall genetic correlations for the relevant combinations ranged from around zero to moderately positive. There was no significant or suggestive QTL for central ligament at this location nor for the other trait involved in these 4 MQR. The highest test statistic for central ligament was found at 24 cM. Also for the 5th MQR on this chromosome, involving fat percentage and kilograms of protein, there was no significant or suggestive QTL in the single trait analysis.
MQR with Correlations Deviating from the Overall Genetic Correlation
For a number of MQR, the correlation caused by the MQR differed largely from the overall genetic correlation. In Table 2
, correlations are indicated where the difference between MQR correlation and the overall genetic correlation was
0.60. Also indicated are MQR where, in addition, the MQR correlation and the overall genetic correlation had opposite sign, and the overall genetic correlation had an absolute value of
0.10. This was the case for 12 MQR (Table 2
).
Most of the MQR with a correlation deviating from the overall genetic correlation have been indicated in the previous sections already. In addition, positive MQR correlations were observed between milking speed and kilograms of protein (BTA 4), and protein percentage and interval calving first insemination (BTA 18). Negative MQR correlations were observed between teat length and nonreturn of daughters 56 d after insemination (BTA 1), fore udder attachment and kilograms of fat (BTA 10), teat length and kilograms of protein (BTA 10), and udder depth and fore udder attachment (BTA 20). The correlation caused by the MQR on BTA 20 affecting udder depth and fore udder attachment was 0.19. On BTA 19, another MQR affecting udder depth and fore udder attachment was observed. The correlation attributable to this MQR was +0.91, which differs largely from the correlation attributable to the MQR on BTA 20.
| DISCUSSION |
|---|
|
|
|---|
Permutation test.
In the method used in this paper, a permutation test is used to identify MQR. Covariances between contrasts are determined in the original data set as well as in a number of permuted data sets. The applied method consists of computation of the covariance between contrasts and comparison of this covariance with the covariance in permuted data sets. Permutation takes into account the characteristics of the data set, e.g., pedigree structure, marker density, and trait distribution. It is, however, not possible to obtain permuted data sets with a QTL affecting only one trait, which would probably be the ideal situation. As indicated by Schrooten and Bovenhuis (2002), the standard deviation of the covariance may be higher in a data set containing a QTL affecting only one trait, as compared with the permuted data sets having no QTL. Depending on the characteristics of the QTL, this phenomenon can lead to higher type I errors than expected. Elevated covariances between contrasts normally indicate the presence of MQR, but in some cases may also be caused by one QTL affecting one trait. Figure 3
shows the relationship between QTL effect and type I error for a trait with heritability of 0.60 and for QTL effects up to 0.5
a. These results were obtained by simulating a granddaughter design of 20 grandsires with 50 sons per grandsire. Each son had 100 daughters with records on 2 traits. Heritability of the traits was 0.6 and 0.35, respectively, and the effect of the QTL on trait 1 ranged from 0.1 to 0.5
a. The second trait was not affected by the QTL. One-half of the grandsires were heterozygous for the bi-allelic QTL. Grandsires were informative for markers at 5-cM intervals on a chromosome of 100 cM. The QTL was located at 30 cM. One thousand data sets were generated, and each data set was permuted 1000 times and analyzed to obtain 95 and 99% confidence intervals for the covariance between contrasts. The number of data sets where the covariance between contrasts was outside the confidence interval was counted to obtain the real type I error. Figure 3
shows that type I error increases with increasing effects of the QTL. The impact on type I error rates is small for QTL effects up to 0.2
a.
|
For a large number of MQR presented in Table 2
, QTL for one of the traits involved have been reported in literature. Most studies only reported on a limited number of traits; therefore, only a limited number of MQR found in the current research could be confirmed from literature, especially those on BTA 6 and BTA 14, which affected production traits. These will be mentioned here, and the QTL from literature affecting one of the traits involved in the MQR from Table 2
will not be discussed.
BTA 6.
Numerous researchers have studied QTL on BTA 6, indicating the presence of multiple QTL on this chromosome (reviewed by Mosig et al., 2001). Spelman et al. (1996) found a QTL affecting both protein percentage and milk yield in 2 grandsire families. Using interval mapping and multiple regression analysis in a daughter design, Ron et al. (2001) found results supporting the hypothesis that one QTL is segregating in 2 different families and that this QTL affects 5 production traits: kilograms of milk, fat percentage, protein percentage, kilograms of fat, and kilograms of protein. This QTL was located close to 55 cM, which corresponds to 40 cM on the map used in our study. Freyer et al. (2002) reported a QTL on BTA 6 that affected both protein yield and fat yield. This QTL was located at 70 cM, which corresponds to 52 cM on the map used in the current study. Freyer et al. (2002) also found QTL at the casein cluster that affected several production traits. The casein cluster corresponds to a location of 95 cM on the map used in the current study. Compared with the MQR involving production traits in our study, QTL presented by Freyer et al. (2002) were located more toward the center of the chromosome.
BTA 14.
Recently, a mutation in the DGAT1 gene on BTA 14 was reported, and an effect on production traits was noted (Grisart et al., 2002). This gene has an effect on several milk production traits and is located in the region between the first 2 markers used in this study, close to 0 cM. The substitution effects for milk production traits were kilograms of milk, 158; fat percentage, +0.17%; protein percentage, +0.04%; kilograms of fat, +5.23; and kilograms of protein, 2.82 (Grisart et al., 2002). Effects were expressed in terms of daughter yield deviations. Based on these effects, an increase of kilograms of milk, for example, is expected to result in a decrease of fat percentage and a decrease of kilograms of fat. In the current study, which consists of partly the same data as used by Grisart et al. (2002), a region on BTA 14 affecting multiple milk production traits was reported (Figure 2
). For all combinations of production traits presented in Figure 2
, the effects on the traits involved were in agreement with the effects of the gene as presented by Grisart et al. (2002), which supports the results presented in this study. The MQR for production traits on BTA 14 in the current study showed most extreme covariances at 7, 11, and 12 cM, which agrees with the location of the DGAT1 gene close to 0 cM.
MQR with Correlations Deviating from the Overall Genetic Correlation
In Table 2
, a number of MQR are shown for which the covariance caused by the MQR and the overall genetic covariance differ largely. The overall genetic covariance is the result of the summation of all covariances attributable to individual QTL, either pleiotropic or closely linked. The contribution of each MQR to the overall genetic covariance depends on the size and the sign of the effects on the traits involved. If the genetic correlation is strongly positive, it is expected that the covariance at most of the individual MQR is positive. There might, however, be a few MQR with a negative contribution to the overall genetic covariance, i.e., MQR with correlations that strongly deviate from the overall genetic correlation. As was shown by Schrooten and Bovenhuis (2002), the statistical power to detect these MQR is increased compared with MQR that do not deviate from overall genetic correlations. For pairs of traits with high genetic correlations, it is expected that there are relatively few MQR that strongly deviate from the overall genetic correlation, but the power to find these MQR is higher (Korol et al., 1995; Schrooten and Bovenhuis, 2002). For cases in which the overall genetic correlation between traits is unfavorable, MAS for MQR that strongly deviate from the overall genetic correlation offers excellent opportunities to achieve progress for both traits. This was demonstrated by de Koning and Weller (1994).
Comparison with Results from Single Trait Analysis
The analysis described in this paper resulted in regions affecting multiple traits, in some cases affecting traits for which no QTL were detected in the single trait analysis. As was discussed by Korol et al. (1995), methods that use information from multiple traits have higher power, which can account for a number of QTL not previously identified. However, in the present study, we did not account for multiple testing when identifying MQR. In total, 3045 tests were performed, and MQR were reported when the covariance was outside the 99% chromosomewise confidence interval. The type I error of 0.01 was divided over both sides of the interval when constructing the confidence interval. In the single trait analysis, 435 tests were performed, and QTL were considered suggestive when the chromosomewise type I error was <0.0345. Combining the number of detected QTL and MQR (Table 1
) with the respective number of tests yields a false discovery rate (Benjamini and Hochberg, 1995) of (0.0345 x 435)/52 = 0.29 for the single QTL analysis and (0.01 x 3045)/59 = 0.52 for the MQR analysis. The latter is relatively high. Therefore, MQR reported here should be considered suggestive rather than significant. However, it was considered important to present these results to assess potential consequences of using the given MQR in MAS.
To compare results from single trait analysis and the analysis described in this paper, consider for example the results for BTA 25. Four MQR involving central ligament were detected at 80 cM, whereas, in the single trait analysis for this trait, a suggestive QTL was located at 24 cM. At first glance, these results do not seem to be in agreement. A within-family analysis on this chromosome for central ligament, however, indicated 2 families with suggestive QTL at 80 cM, in addition to other families with suggestive QTL for central ligament at or around 24 cM (results not shown). Within-family contrasts for other traits, such as rear udder height, were significant in these families. As a result, the covariance between contrasts for these traits was outside the confidence interval for the expected covariance under the null hypothesis. Although these QTL were not detected in the across-family analysis, they can be detected in a within-family analysis and also by the method used in this paper.
The MQR identified in the current study include a number of potentially interesting regions from the point of view of selection. For example, the MQR on BTA 25 should be explored further: favorable effects on kilograms of milk and kilograms of protein are associated with favorable effects on udder traits, especially central ligament. Udder traits such as rear udder height, front teat placement, fore udder attachment, and udder depth are affected by MQR on BTA 19, which could be exploited in MAS. Favorable effects of MQR on BTA 13, on udder traits such as udder depth and fore udder attachment, are associated with favorable effects on fat and protein percentages.
| CONCLUSIONS |
|---|
|
|
|---|
| APPENDIX |
|---|
|
|
|---|
|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication June 3, 2003. Accepted for publication March 15, 2004.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. S. Lund, B. Guldbrandtsen, A. J. Buitenhuis, B. Thomsen, and C. Bendixen Detection of Quantitative Trait Loci in Danish Holstein Cattle Affecting Clinical Mastitis, Somatic Cell Score, Udder Conformation Traits, and Assessment of Associated Effects on Milk Yield J Dairy Sci, October 1, 2008; 91(10): 4028 - 4036. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Kolbehdari, Z. Wang, J. R. Grant, B. Murdoch, A. Prasad, Z. Xiu, E. Marques, P. Stothard, and S. S. Moore A Whole-Genome Scan to Map Quantitative Trait Loci for Conformation and Functional Traits in Canadian Holstein Bulls J Dairy Sci, July 1, 2008; 91(7): 2844 - 2856. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |