|
|
||||||||




* Department of Animal Science, Agricultural University of Norway, N-1432 Aas, Norway
Centre UdL-IRTA, Department of Animal Production, 25198 Lleida, Spain
Department of Morphology, Genetics and Aquatic Biology, The Norwegian School of Veterinary Science, N-0033 Oslo, Norway
Agricultural Research Centre MTT, Animal Production Research, 31600 Jokioinen, Finland
|| Chair of Animal Breeding, Technical University of Munich, D-85380 Freising-Weihenstephan, Germany
Corresponding author:
Hanne Gro Olsen; e-mail:
hanne-gro.olsen{at}ihf.nlh.no.
| ABSTRACT |
|---|
|
|
|---|
Key Words: quantitative trait loci milk cattle
Abbreviation key: BTA = bovine (Bos taurus) chromosome, CI = confidence interval, EBV = estimated breeding value, NCM = Norwegian Dairy Cattle Map
| INTRODUCTION |
|---|
|
|
|---|
Recently, a male genetic marker map based on six paternal half-sib families from the Norwegian Dairy Cattle population has been constructed (Våge et al., 2000). The aim of the present work was to use the genotypic information of these six families in an autosomal genome scan for QTL affecting milk production traits in the Norwegian Dairy Cattle population.
| MATERIALS AND METHODS |
|---|
|
|
|---|
For the present study, six paternal half-sib families were chosen because of their relatively large number of progeny-tested sons. The total number of sons in the study were 284, ranging from 32 sons of sire 2052 to 71 sons of sire 2005 (Table 1
). The average number of effective daughters in the dataset was as large as 843, and only one son had less than 150 daughters.
|
Marker Map
The Norwegian Cattle Map (NCM) covers all 29 autosomal chromosomes and summarizes to a total length of 2682 cM (Våge et al., 2000). The map consists of 288 markers, with an average of 131 informative sons per marker. The average length of marker intervals is 12.5 cM, and approximately 82% of the intervals are shorter than 20 cM. The chromosome assignment, locus order and interval lengths agree with the most comprehensive linkage map available in cattle (Kappes et al., 1997). The marker map was constructed using six half-sib families (Table 1
), with marker information for sires and their sons only. Detailed information about the NCM is available at http://www.nlh.no/Institutt/IHF/Genkartstorfe/.
The marker density of the NCM was increased for chromosome 6 by genotyping four additional microsatellites (FBN9, FBN12, FBN13 and FBN15) with primers and PCR conditions as described by Weikard et al. (1997). Linkage order and map distances among markers were estimated using the CRI-MAP 2.4 program (Green et al., 1990), with map distances based on Haldanes mapping function.
Performance Data
In Norway, progeny testing is based on more than 250 daughters per son, which is about five times higher than commonly used in other European or American populations. Estimated breeding values (EBV) of the sons from the six families were used for performance data in the analyses. EBV for the five milk production traits (milk yield, fat percentage, fat yield, protein percentage and protein yield) were available from the national genetic evaluation in June 2000 carried out by GENO Breeding and A.I. Association, and evaluated using BLUP with a single trait sire-maternal grandsire model.
Statistical Analyses
The multiple marker regression method described by de Koning et al. (2001) was applied for QTL mapping. As the method assumes that all sires and dams are unrelated, and that each dam has only one offspring, information about the extensive relationships between animals was not included in the analysis. A weighted least squares analysis was utilized, where the weights were the number of daughters per son.
First, the most likely linkage phase of each of a sires chromosomes was estimated based on marker information of the sire and his sons. For a given heterozygous marker, the allele inherited from sire to son could be determined unequivocally if the son possessed only one of the sires alleles. If the son was heterozygous for the same alleles as his sire, the most likely paternal allele was estimated from the flanking markers. The most likely haplotype of the sire was taken as the linkage phase minimizing the number of recombinant offspring. One of the haplotypes of each chromosome of the sire was arbitrarily chosen as the first haplotype. Then, for each son, the probability of inheriting the sires first haplotype was calculated at each cM along each chromosome conditional on informative flanking markers.
At each cM, a QTL was fitted using least squares methodology by regressing the trait score of the individual on the probability of receiving the first haplotype at that point, according to the model:
![]() |
where
| Yij | = | EBV of son j within sire family i (i = 1 to 6);
|
| µ | = | overall mean;
|
| ai | = | fixed polygenic effect of sire i;
|
| bi | = | regression coefficient within sire i; (i.e. allele substitution effect for a putative QTL)
|
| Xij | = | probability of son j of inheriting the first haplotype from sire i conditional on the informative flanking markers of the son; and
|
| eij | = | residual effect.
|
This model assumed one QTL with an additive effect and no interaction because dominance cannot be detected using a half-sib structure, and because detection of epistasis would require several QTL. The residual effects were assumed normally distributed.
Regression was nested within sire because the QTL and marker alleles as well as their respective linkage phases might be different across sires. For each sire, a within-family test statistic for the regression coefficients was calculated at each cM along the genome, and the position with the maximum F-ratio was considered to be the most likely QTL location for that family. An across-family F-ratio was then obtained by pooling the regression mean squares over grandsires and dividing it by the residual mean squares. Again, this analysis was repeated at each position along the genome to find an estimate of the overall best position. The across family test statistic was used for determining whether a QTL was segregating in the tested material, and to find the most likely overall position for this QTL. The within-family test statistic revealed which families were heterozygous for the QTL, and could give an indication of whether one of more QTL were segregating if the most likely QTL position differed between families.
Chromosome-wise significance values (Pchromosome) for both the within-family and the across-family analyses were determined empirically applying the permutation test of Churchill and Doerge (1994), with 10,000 permutations per chromosome. As performing a full permutation test for all 29 chromosomes simultaneously was too computer demanding, the genome-wise significance levels, taking into account the multiple testing of the entire genome, were derived by a Bonferroni correction as:
![]() |
where Pchromosome is the significance level of a specific chromosome and r was obtained as the length of the specific chromosome divided by the total length of the autosomal genome (de Koning et al., 1999).
Empirical 95% confidence intervals (CI) for the QTL positions were obtained by applying the bootstrapping approach proposed by Visscher et al. (1996).
In addition, chromosome 6 was tested for the presence of two QTL as described by Velmala et al. (1999). A model testing the presence of two QTL vs. none was fitted as:
![]() |
where the subscripts 1 and 2 refer to the position of the first and the second QTL, respectively, and the other terms are the same as for the one-QTL model. The two-QTL model was fitted by grid searching, i.e. each combination of 5-cM positions was tested with a minimum distance of 10 cM between the two intervals. The most probable positions of the two QTL were taken as the combination giving the highest F-value. Empirical significance values were determined by permutation. If the test statistic for two vs no QTL was significant, the presence of two vs one QTL was tested with an F-test on the mean squares from the two-QTL model compared to that of the best one-QTL model. Whether the two-QTL model explained more of the variation than the one-QTL model was evaluated from a standard F table, with the number of grandsires as nominator and the total number of sons minus three times the number of grandsires as the denominator (Velmala et al, 1999).
| RESULTS |
|---|
|
|
|---|
|
|
Rather large QTL effects were detected for both percentage traits. The allele substitution effect on protein percentage was estimated to be approximately 0.16% and 0.20% (1.6 and 1.9 genetic standard deviations) in family 2402 and family 2463, respectively, whereas the effect on fat percentage was approximately 0.17% and 0.42% (0.8 and 2.1 genetic standard deviations) in the same two families. Segregation of QTL for milk yield was detected only in family 2463, with an allele substitution effect of approximately 350 kg milk, or 0.85 genetic standard deviations.
Confidence intervals for the most likely QTL positions on chromosome 6 were estimated by a bootstrapping approach. The 95% CI for the two percentage traits were almost identical in position and size, covering positions 29 to 45 cM for fat percentage and positions 30 to 46 for protein percentage. An overlapping, but much wider, confidence interval was found for milk yield, encompassing positions 22 through 128.
For the other chromosomes, a putative QTL affecting milk yield on BTA18 was the most significant, with a genome-wise significance level of 8%. Although segregation was detected in one family only, the substitution effect in this family corresponded to approximately 396 kg or 0.96 genetic standard deviations.
| DISCUSSION |
|---|
|
|
|---|
The common haplotype surrounding the BTA6 QTL found in families 2402 and 2463 seems to cause a marked reduction in fat and protein percentages, and a modest increase in milk yield. Such an effect might be explained by gene(s) affecting the lactose yield. As lactose is the major osmotic component of the milk, increased lactose yield would result in increased milk yield, but unaltered yield of fat and protein. Thus, fat and protein percentages would be reduced. However, there are no indications that genes with such effects are present in the relevant area of BTA6.
Unfavorable genetic correlation between increased milk production and enhanced frequency of mastitis is well known. In a recent study of the same families (Klungland et al., 2001), putative QTL affecting clinical mastitis were localized to BTA3, 4, 6, 14 and 27, respectively. Whereas the mastitis BTA3 QTL is positioned on the other end of the chromosome, the BTA6 QTL has similar location for both traits. Based on these findings it could be argued that the mastitis QTL on BTA6 is caused by high milk yield, rather than a change in the immune response or, alternatively, that the relatively wide confidence interval for the mastitis QTL could contain several independent QTL. The fact that the QTL for milk production and clinical mastitis were detected in different families supports the latter hypothesis.
Several studies have reported another QTL on chromosome 6, which is close to the casein gene cluster (e.g. Kühn et al., 1996; Velmala et al., 1999). As a previous study in Norwegian Dairy Cattle also reported an association between protein yield and a specific casein haplotype (Lien et al., 1995), a test for the presence of two QTL on chromosome 6 was performed. However, no support for a second QTL in the casein region was found. A possible explanation for this discrepancy may be, in addition to the use of different statistical approaches, that the animal material was somewhat different in the two studies.
In recent years, a large number of breeds have been tested for the presence of QTL using an array of different methods. From these studies, as many as 20 chromosomes have been suggested to harbor QTL that affect milk traits (Mosig et al. 2001). Typically, only a few of the putative QTL detected in a given study are confirmed by other studies. Of the putative QTL reported in this study, the findings on chromosomes 3 and 20 are strongly supported by the literature. On BTA3, one family was segregating for a putative QTL affecting fat yield near the centromere. Previously, Heyen et al. (1999) found significant marker effects for at least one QTL affecting fat percentage, protein yield and protein percentage in the same region in the US Holstein population. Lipkin et al. (1998) and Zhang et al. (1998) also reported a QTL for protein percentage at the centromeric end of BTA3. On BTA20, a QTL affecting protein yield was segregating in one family in our study, between markers BMS2361 and UWCA26. This confirms the results found in Holstein (Arranz et al., 1998; Georges et al., 1995; Zhang et al., 1998) and in Finnish Ayrshire (Mäki-Tanila et al., 1998). Although the across-family test statistic for BTA3 and BTA20 did not exceed the genome-wise significance level in our study, the presence of similar patterns in different studies strongly support these positions as putative QTL for milk production traits. Other studies have reported the presence of QTL on chromosomes showing no significant effect in this study, e.g. BTA1 (e.g. Georges et al., 1995; Nadesalingam et al., 2001), BTA9 (Georges et al., 1995; Vilkki et al., 1997; Zhang et al., 1998; Wiener et al., 2000) and BTA14 (e.g. Coppieters et al., 1998; Heyen et al., 1999). The BTA14 QTL has later been fine-mapped (Riquet et al., 1999; Looft et al, 2001), and a candidate gene and a causal mutation have been detected (Grisart et al., 2001). This inconsistency across studies can partly be explained by the many different statistical methods and designs used in the studies, and by how the significance thresholds are set to account for the multiple testing performed in each study. For most breeds, family sizes are lower than optimum for QTL detection, thus detection power is low in most studies. As a consequence, many QTL are missed, and the effect of those detected is often overestimated. Additionally, as most QTL are segregating in only a few families, the sires in the studied sample might not be heterozygous for a QTL segregating in the population as a whole. The density and informativity of the surrounding markers could also be insufficient for linkage detection. Confirmation of results in independent studies is therefore needed to reveal which effects are authentic. Importantly, analysis of data from different breeds or populations might also provide additional insight into the genetics controlling the trait of interest.
Most genome scans using similar experimental designs report confidence intervals of at least 20 to 30 cM. Thus, the 16 cM intervals estimated for the percentage traits on BTA6 are surprisingly narrow, a fact that probably reflects the large effect of the QTL. These intervals, however, cannot be substantially reduced by denser marker maps due to the lack of recombinations needed for positioning the QTL. Marker information may be utilized in e.g., marker-assisted selection, separating linked QTL, estimation of pleiotropic effects and positional cloning. However, these implementations require the QTL to be mapped to an interval of 1 to 2 cM to be efficient. Other methods with much higher precision than the genome scan are therefore needed.
A more precise location of the QTL could be achieved by identifying additional families segregating the QTL and increase the number of genotyped sons within these families. This is not a realistic approach in Norwegian Dairy Cattle, however, as the number of sire families with a large number of progeny-tested sons is limited. Alternatively, the granddaughter design could be supplemented by genotyping a very large number of daughters in a daughter design. Recently, Ron et al. (2001) were able to confine a QTL for protein percentage on BTA6 to a confidence interval of only 4 cM for two families of the Israeli Holstein population by using this approach. However, this approach is not easily applicable to the Norwegian Dairy Cattle population as herd size is small, making the collection of DNA samples very expensive. Thus, a more realistic approach for increasing the precision of the estimates is to utilize historical recombinations in linkage disequilibrium or identical by descent mapping (Meuwissen et al., in press; Riquet et al., 1999) to refine the location of the BTA6 QTL in Norwegian Dairy Cattle. When reaching an acceptable refinement, the QTL positions can be compared with comparative regions of the human genome in an effort to identify positional candidate genes.
For the putative QTL on bovine chromosomes 3 and 6, no obvious candidate genes for milk production were identified by examination of the corresponding human genome information. On chromosome 20, however, the genes encoding the prolactine and growth hormone receptors are mapped to the same interval as was reported in this and other studies, and are therefore possible candidates for further research.
| CONCLUSION |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication February 16, 2002. Accepted for publication May 24, 2002.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. Lillehammer, B. J. Hayes, T. H. E. Meuwissen, and M. E. Goddard Gene by environment interactions for production traits in Australian dairy cattle J Dairy Sci, August 1, 2009; 92(8): 4008 - 4017. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Momke, H. Scholz, K. Doll, J. Rehage, and O. Distl Mapping Quantitative Trait Loci for Left-Sided Displacement of the Abomasum in German Holstein Dairy Cows J Dairy Sci, November 1, 2008; 91(11): 4383 - 4392. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. D. Daetwyler, F. S. Schenkel, M. Sargolzaei, and J. A. B. Robinson A Genome Scan to Detect Quantitative Trait Loci for Economically Important Traits in Holstein Cattle Using Two Methods and a Dense Single Nucleotide Polymorphism Map J Dairy Sci, August 1, 2008; 91(8): 3225 - 3236. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Lillehammer, M. E. Goddard, H. Nilsen, E. Sehested, H. G. Olsen, S. Lien, and T. H. E. Meuwissen Quantitative Trait Locus-by-Environment Interaction for Milk Yield Traits on Bos taurus Autosome 6 Genetics, July 1, 2008; 179(3): 1539 - 1546. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Druet, S. Fritz, D. Boichard, and J. J. Colleau Estimation of genetic parameters for quantitative trait Loci for dairy traits in the French holstein population. J Dairy Sci, October 1, 2006; 89(10): 4070 - 4076. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Khatib, S. D. Leonard, V. Schutzkus, W. Luo, and Y. M. Chang Association of the OLR1 Gene with Milk Composition in Holstein Dairy Cattle J Dairy Sci, May 1, 2006; 89(5): 1753 - 1760. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Y. Chen, Q. Zhang, C. C. Yin, C. K. Wang, W. J. Gong, and G. Mei Detection of Quantitative Trait Loci Affecting Milk Production Traits on Bovine Chromosome 6 in a Chinese Holstein Population by the Daughter Design J Dairy Sci, February 1, 2006; 89(2): 782 - 790. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Cohen-Zinder, E. Seroussi, D. M. Larkin, J. J. Loor, A. E.-v. d. Wind, J.-H. Lee, J. K. Drackley, M. R. Band, A.G. Hernandez, M. Shani, et al. Identification of a missense mutation in the bovine ABCG2 gene with a major effect on the QTL on chromosome 6 affecting milk yield and composition in Holstein cattle Genome Res., July 1, 2005; 15(7): 936 - 944. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Weikard, C. Kuhn, T. Goldammer, G. Freyer, and M. Schwerin The bovine PPARGC1A gene: molecular characterization and association of an SNP with variation of milk fat synthesis Physiol Genomics, March 21, 2005; 21(1): 1 - 13. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Szyda, Z. Liu, F. Reinhardt, and R. Reents Estimation of Quantitative Trait Loci Parameters for Milk Production Traits in German Holstein Dairy Cattle Population J Dairy Sci, January 1, 2005; 88(1): 356 - 367. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. G. Olsen, S. Lien, M. Gautier, H. Nilsen, A. Roseth, P. R. Berg, K. K. Sundsaasen, M. Svendsen, and T. H. E. Meuwissen Mapping of a Milk Production Quantitative Trait Locus to a 420-kb Region on Bovine Chromosome 6 Genetics, January 1, 2005; 169(1): 275 - 283. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. J. Boettcher, A. Caroli, A. Stella, S. Chessa, E. Budelli, F. Canavesi, S. Ghiroldi, and G. Pagnacco Effects of Casein Haplotypes on Milk Production Traits in Italian Holstein and Brown Swiss Cattle J Dairy Sci, December 1, 2004; 87(12): 4311 - 4317. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. G. Olsen, S. Lien, M. Svendsen, H. Nilsen, A. Roseth, M. Aasland Opsal, and T. H. E. Meuwissen Fine Mapping of Milk Production QTL on BTA6 by Combined Linkage and Linkage Disequilibrium Analysis J Dairy Sci, March 1, 2004; 87(3): 690 - 698. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |