|
|
||||||||
1 USDA-ARS, Bovine Functional Genomics Laboratory, Beltsville, MD 20705
2 Department of Animal Sciences, University of Illinois, Urbana 61801
3 Department of Animal Science, University of Minnesota, St. Paul 55108
4 USDA-ARS, Animal Improvement Programs Laboratory, Beltsville, MD 20705
5 Institute of Animal Sciences, ARO, The Volcani Center, Bet Dagan 50250 Israel
Corresponding author: Melissa S. Ashwell; e-mail: mashwell{at}anri.barc.usda.gov.
| ABSTRACT |
|---|
|
|
|---|
Key Words: quantitative trait locus pregnancy rate milk production traits
Abbreviation key: AIPL = Animal Improvement Programs Laboratory, BTA = Bos taurus autosome, DBDR = Dairy Bull DNA Repository, PL = productive life, DGAT1 = acylCoA:diacylglycerol acyltransferase gene, GHR = growth hormone receptor gene
| INTRODUCTION |
|---|
|
|
|---|
In many species, especially livestock species, studies designed for identification of QTL are based on crosses of genetically distinct breeds or inbred lines. This approach has been less used in dairy cattle research because development of experimental populations is prohibitively expensive and the generation interval is long. In addition, the population structure that results from the use of elite AI bulls as sires is quite powerful and requires only modest amounts of genotyping for QTL identification. In 1991, the Dairy Bull DNA Repository (DBDR; Da et al., 1994) was initiated as a collection of semen from 35 dairy grandsire families that was available to collaborators for studies to identify QTL affecting economically important traits in Holstein cattle. Two genome-wide scans were conducted in DBDR families, focusing on identification of QTL affecting the milk production traits and SCS. Originally two groups conducted and reported independent research findings (Heyen et al. 1999; Ashwell et al. 2001). These datasets have now been merged for further analysis allowing a more thorough coverage of the bovine genome. Originally the independent genome scans were analyzed one marker at a time using analysis of variance methods, with interval regression analysis conducted only on specific chromosomes. Those studies lacked thorough grandsire-specific chromosomal coverage; the merged results presented here increased the coverage by combining genotypes from the different sets of DNA markers, such that interval regression analysis could be more thoroughly carried out on all 29 autosomes.
Since the two datasets were generated and merged, researchers with the USDA-ARS Animal Improvement Programs Laboratory (AIPL) have calculated genetic evaluations for female fertility (VanRaden et al., 2003). The female fertility evaluations are expressed as the pregnancy rate, which is defined as the percentage of nonpregnant cows that become pregnant during each 21-d period. Heritability for this trait has been estimated at 4% (VanRaden et al., 2003), making this an ideal candidate trait for a marker-assisted selection program.
The objective of this study was to identify QTL affecting milk production traits using the merged set of DBDR genotypes as well as provide the first report of putative QTL affecting the pregnancy rate trait.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Originally the two research groups (University of Illinois, Urbana/ARO, Israel and USDA-ARS, Beltsville) conducted independent genome scans. Each group selected eight families for study, six that were common across both studies (families 1 to 5, 8) and four families that were studied by only one of the two groups (families 6 and 7 by University of Illinois/ARO; families 9 and 12 by USDA-ARS).
Genotyping
DNA was extracted from semen as previously described (Ashwell et al., 1996; Heyen et al., 1999). The PCR was carried out using either fluorescently-labeled primers or by [
-32P] incorporation into PCR products as previously described (Ashwell et al., 1996; Heyen et al., 1999). For each individual genome scan, microsatellite markers were selected at approximately 20-cM intervals from published bovine maps (Bishop et al., 1994; Barendse et al., 1994; Ma et al., 1996; Barendse et al., 1997; Kappes et al., 1997). Heyen et al. (1999) selected 174 markers; Ashwell et al. (2001) selected 232 markers, 38 of which were genotyped by both groups. Duplicate markers were treated as separate loci because allele-calling methods were not consistent across both datasets. When the two datasets were merged, genome coverage was estimated to be 2713.5 cM (90%), assuming a 3000-cM genome. The average number of markers per chromosome was 12.5, with an average spacing of 7.4 cM.
Phenotypic Data
Data for milk yield and composition, SCS, and productive life (PL) collected through November 2001 were processed as part of the routine USDA/DHIA genetic evaluation procedure by the AIPL of USDA-ARS. The female fertility trait, pregnancy rate, is a new genetic evaluation being calculated by AIPL. Evaluations are expressed as the pregnancy rate, which is defined as the percentage of nonpregnant cows that become pregnant during each 21-d period (VanRaden and Tooker, 2003; VanRaden et al., 2003). The PTA values for pregnancy rate trait ranged from -4.9 to 3.3%, with an average of -0.4% and a standard deviation of 1.1 for the sires used in this study. Pregnancy status is determined from the date of last breeding and is verified using the next calving date or veterinary diagnosis when available. Cows that are sold for reproductive problems are assumed not pregnant. Pregnancy status is determined only within the first 250 d of lactation because many cows not pregnant by 250 d are not bred again but instead are culled. Mean pregnancy rate for Holsteins in 1995 was 23% and SD of sire evaluations was 1%. These evaluations for pregnancy rate are obtained from data on days open, and an increase of 4 d open equals a 1% decrease in pregnancy rate. Data since 1960 from the first 5 lactations are included, for a total of 40 million records. Data are adjusted for month of calving, age at first calf, and lactation number and are computed using the same animal model evaluation programs used for yield traits, productive life, and SCS (VanRaden and Tooker, 2003; VanRaden et al., 2003). The present study included 3 mo less data than the first official evaluations that were released in February 2003.
Statistical Analysis
Data from a total of eight traits were analyzed using a regression approach originally described by Haley and Knott (1992). A web-based version of this regression interval mapping method is now available (Seaton et al., 2001). The software, QTL Express (http://qtl.cap.ed.ac.uk), analyzes data from F2, half-sib, and sib-pair families to detect QTL. The software allows the user to fit one or two QTL in the model and includes tools for permutation and bootstrap analyses to calculate chromosome-wise significance thresholds and 95% confidence intervals, respectively. For this study, 1000 permutations were studied for each trait to determine the P < 0.05 and P < 0.01 chromosome-wise significance thresholds, and the regression interval analysis was conducted at 1-cM intervals along the chromosome. Both within- and across-family analyses were conducted, fitting one and two QTL in the model. Data included daughter deviations for milk, fat, and protein yield, fat and protein percentage, SCS, and PL, weighted by their respective reliabilities. At the time this analysis was conducted, no daughter deviations were available for pregnancy rate, therefore, PTA for pregnancy rate were used. Within family analyses correspond to a contrast of two sire marker alleles with one degree of freedom, while the across-family analysis tests for the evidence of different marker effects across families. Marker effects are nested within sire family so that no assumptions are made regarding the phase of QTL allelic effects in the across-family analysis.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
Recently, the same group (Blott et al., 2003) identified a mutation in the growth hormone receptor (GHR) gene that is associated with an effect on milk yield and milk composition. GHR mapped to BTA20 at 43 cM on their map and is associated with a major effect on milk protein percentage. Results from the DBDR merged dataset identified a QTL affecting milk protein percentage at 51 cM, based on the MARC map (Kappes et al., 1997). The QTL identified in our DBDR study most likely corresponds to the QTL identified by Blott et al. (2003).
Marker Effects on SCS and Productive Life
No evidence favoring the two QTL models over the one QTL model was found in the within and across-family analyses for SCS and productive life. Eight significant effects on SCS were identified within families on eight chromosomes when one QTL was modeled (Table 1
). Four significant effects on SCS were identified in the across-family analysis (Table 2
). For SCS, a total of 290 chromosome-family combinations were tested, with three significant tests expected by chance at P < 0.01 and eight were observed. In the across-family analysis, 29 tests were conducted, with less than one significant test expected by chance, and four were observed. All significant effects on this trait were detected in the independent DBDR studies, except for the within-family associations on BTA15 and BTA20. Boichard and Bishop (1997) also detected an association on BTA15 at 36 cM, close to the same location reported here. Although many groups have detected QTL on BTA20, none have been reported affecting SCS.
A putative QTL affecting SCS was detected on BTA26 at 0 cM in families 3 and 7. A SCS QTL on BTA26 was reported by Zhang et al. (1998) and Boichard and Bishop (1997); however, they report the location to be toward the telomere, between TGLA429 and BM804. DBDR families were not genotyped for marker TGLA429; however, families 3 and 7, as well as 4, 5, and 6, were heterozygous at BM804 and all their sons were genotyped. Therefore, one explanation for these findings is that two QTL affecting SCS are located on BTA26 at opposite ends of the chromosome and were detected in different families.
Only one significant effect on PL, modeling one QTL, was detected. The putative QTL was detected within family 1 on BTA7 (Table 1
). Neither of the independent scans detected a significant association affecting PL, however the Heyen et al. (1999; see also http://cagst.animal.uiuc.edu/genemap/WEB/Table1.html) study detected a suggestive association (P < 0.1) in another family. Comparison of this finding with those from other groups is difficult because few groups have tried to identify PL QTL and the phenotype and genetic evaluation systems differ from country to country (van der Linde and de Jong, 2002). A total of 29 significance tests were conducted, with three significant tests expected by chance at P < 0.01; therefore the effect detected on BTA7 may be a false-positive association.
Marker Effects on Pregnancy Rate
The PTA for pregnancy rate (VanRaden et al., 2002) were recently computed, allowing, for the first time, identification of putative QTL affecting this trait. Seven significant effects were identified within families on six chromosomes (Table 3
). A total of 290 significant tests were conducted, with three significant tests expected by chance at P < 0.01 and seven observed. The most significant within-family effect (F statistic = 14.1) was located on BTA18 in DBDR family 4 (Figure 1
). The estimate of the difference of the allelic effect for this QTL was 0.57 (± 0.15) units (Table 3
). The only significant associations previously detected on this chromosome affected rump angle; however, this was in a different DBDR family (Ashwell et al., 2001). No significant effects were detected across families. This is likely if the heterozygosity at the QTL is low due to selection (Heyen et al., 1999) and QTL will only be detected within families and not across the families. No evidence was found favoring the two QTL model over the one QTL model in the within- and across-family analyses for this trait.
|
|
Schrooten et al. (2000) analyzed "nonreturn daughters" and "interval calving to first insemination" fertility traits in their QTL study. Suggestive linkage was found on BTA2 and BTA9 affecting the nonreturn daughters trait. These locations cannot be confirmed with the pregnancy rate findings. Larsson and Andersson-Eklund (2002) analyzed "interval calving to first insemination," fertility treatments, and a combined fertility index in their partial genome scan to identify QTL affecting reproduction traits in Swedish dairy cattle. Putative QTL affecting the interval from calving to first insemination were identified on BTA3, BTA11, and BTA25 and cannot be confirmed by the pregnancy rate findings. They also detected QTL on BTA6 affecting fertility treatments and the combined fertility index. Schrooten et al. (2000) also identified suggestive linkage on BTA6 and BTA17 affecting the interval from calving to first insemination. The putative QTL on BTA6 identified by Schrooten et al. (2000), located at 106 cM in their study, may be supported by the pregnancy rate QTL detected on this chromosome at 122 cM in our study. In addition, BTA6 is also known to harbor QTL effecting traits such as milk production (Kühn et al., 1999; Nadesalingam et al., 2001; Ron et al., 2001), clinical mastitis (Klungland et al., 2001) and calving ease (Schrooten et al., 2000; Larsson and Andersson-Eklund, 2002). Once the genes underlying these QTL are identified, researchers may be able to determine the pleiotropic effects on important traits in dairy breeding.
The putative pregnancy rate QTL on BTA14 may be explained by effects of the DGAT1 gene discussed above. DGAT1 maps to the centromeric end of BTA14 and is included in the region where a pregnancy rate QTL was detected. DGAT1 catalyzes the last step in triglyceride synthesis and Grisart et al. (2002) identified a mutation in this gene that is associated with a major effect on milk fat content. It has been shown that changes in the fat:protein ratio in milk during early lactation has a negative effect on fertility (de Vries and Veerkamp, 2000). Therefore mutations in the DGAT1 gene may have pleiotropic effects on the pregnancy rate.
Previously, no fertility QTL were reported on BTA27. However, a previous study (Ashwell et al., 2001) detected a QTL affecting the conformation trait, dairy form, on the telomeric end of BTA27 in the same region as the pregnancy rate QTL. Casas et al. (2002) also identified a putative QTL affecting marbling in beef cattle in the same region of chromosome 27. Therefore, results suggest at least one gene affecting fat metabolism is located on this chromosome that may also affect fertility. Rogers et al. (1999) studied the genetic correlations between type traits and three groups of diseases: reproductive, foot and leg, and metabolic and digestive diseases. Results from this study showed that the genetic correlation between dairy form and the three disease categories were negative and moderate in magnitude. Metabolic disorders (such as milk fever and ketosis) are indicators of negative energy balance (Rogers et al., 1999), which is known to affect reproductive performance (de Vries and Veerkamp, 2000). Therefore, selection for increased dairy form may lead to cows that are more prone to reproductive and metabolic diseases (Rogers et al., 1999).
It should be emphasized that the results associated with pregnancy rate are preliminary and require validation. However, based on these results, one can speculate that identification of the genes underlying the pregnancy rate QTL may actually be genes affecting body condition and the metabolic state of dairy animals. In several cases there are effects on milk production traits in regions where putative pregnancy rate QTL have been detected, so it may be difficult to improve female fertility without sacrificing some milk production. However, it may be more economically advantageous to sacrifice some gains in milk production in order to be able to breed cows in a shorter period of time.
These data are just the first step in identifying and locating QTL. Work is underway to validate these QTL. Once these results are validated, fine-mapping studies will begin, so the genes responsible for the observed effects can be identified and incorporated into marker-assisted selection programs.
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication December 20, 2002. Accepted for publication June 1, 2003.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
R. M. Demeter, G. C. B. Schopen, A. G. J. M. O. Lansink, M. P. M. Meuwissen, and J. A. M. van Arendonk Effects of milk fat composition, DGAT1, and SCD1 on fertility traits in Dutch Holstein cattle J Dairy Sci, November 1, 2009; 92(11): 5720 - 5729. [Abstract] [Full Text] [PDF] |
||||
![]() |
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] |
||||
![]() |
J. B. Cole, P. M. VanRaden, J. R. O'Connell, C. P. Van Tassell, T. S. Sonstegard, R. D. Schnabel, J. F. Taylor, and G. R. Wiggans Distribution and location of genetic effects for dairy traits J Dairy Sci, June 1, 2009; 92(6): 2931 - 2946. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. K. Hoglund, B. Guldbrandtsen, G. Su, B. Thomsen, and M. S. Lund Genome scan detects quantitative trait loci affecting female fertility traits in Danish and Swedish Holstein cattle J Dairy Sci, May 1, 2009; 92(5): 2136 - 2143. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. M.-R. Arnould, H. Soyeurt, N. Gengler, F. G. Colinet, M. V. Georges, C. Bertozzi, D. Portetelle, and R. Renaville Genetic analysis of lactoferrin content in bovine milk J Dairy Sci, May 1, 2009; 92(5): 2151 - 2158. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. M. VanRaden, C. P. Van Tassell, G. R. Wiggans, T. S. Sonstegard, R. D. Schnabel, J. F. Taylor, and F. S. Schenkel Invited Review: Reliability of genomic predictions for North American Holstein bulls J Dairy Sci, January 1, 2009; 92(1): 16 - 24. [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] |
||||
![]() |
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] |
||||
![]() |
E. Casas and G. D. Snowder A putative quantitative trait locus on chromosome 20 associated with bovine pathogenic disease incidence J Anim Sci, October 1, 2008; 86(10): 2455 - 2460. [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] |
||||
![]() |
J. Conington, G. Cao, A. Stott, and L. Bunger Breeding for resistance to mastitis in United Kingdom sheep, a review and economic appraisal Vet Rec., March 22, 2008; 162(12): 369 - 376. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. I. Weller, M. Golik, S. Reikhav, R. Domochovsky, E. Seroussi, and M. Ron Detection and Analysis of Quantitative Trait Loci Affecting Production and Secondary Traits on Chromosome 7 in Israeli Holsteins J Dairy Sci, February 1, 2008; 91(2): 802 - 813. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. D. Nkrumah, E. L. Sherman, C. Li, E. Marques, D. H. Crews Jr., R. Bartusiak, B. Murdoch, Z. Wang, J. A. Basarab, and S. S. Moore Primary genome scan to identify putative quantitative trait loci for feedlot growth rate, feed intake, and feed efficiency of beef cattle J Anim Sci, December 1, 2007; 85(12): 3170 - 3181. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. S. Lund, G. Sahana, L. Andersson-Eklund, N. Hastings, A. Fernandez, N. Schulman, B. Thomsen, S. Viitala, J. L. Williams, A. Sabry, et al. Joint Analysis of Quantitative Trait Loci for Clinical Mastitis and Somatic Cell Score on Five Chromosomes in Three Nordic Dairy Cattle Breeds J Dairy Sci, November 1, 2007; 90(11): 5282 - 5290. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Gutierrez-Gil, M. F. El-Zarei, Y. Bayon, L. Alvarez, L. F. de la Fuente, F. S. Primitivo, and J. J. Arranz Short Communication: Detection of Quantitative Trait Loci Influencing Somatic Cell Score in Spanish Churra Sheep J Dairy Sci, January 1, 2007; 90(1): 422 - 426. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Casas and R. T. Stone Putative quantitative trait loci associated with the probability of contracting infectious bovine keratoconjunctivitis J Anim Sci, December 1, 2006; 84(12): 3180 - 3184. [Abstract] [Full Text] [PDF] |
||||
![]() |
O. Cobanoglu, I. Zaitoun, Y. M. Chang, G. E. Shook, and H. Khatib Effects of the signal transducer and activator of transcription 1 (STAT1) gene on milk production traits in Holstein dairy cattle. J Dairy Sci, November 1, 2006; 89(11): 4433 - 4437. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. S. Sharma, G. B. Jansen, N. A. Karrow, D. Kelton, and Z. Jiang Detection and characterization of amplified fragment length polymorphism markers for clinical mastitis in canadian holsteins. J Dairy Sci, September 1, 2006; 89(9): 3653 - 3663. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Holmberg and L. Andersson-Eklund Quantitative trait Loci affecting fertility and calving traits in Swedish dairy cattle. J Dairy Sci, September 1, 2006; 89(9): 3664 - 3671. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Longeri, M. Polli, M. G. Strillacci, A. B. Samore, and M. Zanotti Short communication: quantitative trait loci affecting the somatic cell score on chromosomes 4 and 26 in Italian Holstein cattle. J Dairy Sci, August 1, 2006; 89(8): 3175 - 3177. [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] |
||||
![]() |
S. Leonard, H. Khatib, V. Schutzkus, Y. M. Chang, and C. Maltecca Effects of the Osteopontin Gene Variants on Milk Production Traits in Dairy Cattle J Dairy Sci, November 1, 2005; 88(11): 4083 - 4086. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. S. Ashwell, D. W. Heyen, J. I. Weller, M. Ron, T. S. Sonstegard, C. P. Van Tassell, and H. A. Lewin Detection of Quantitative Trait Loci Influencing Conformation Traits and Calving Ease in Holstein-Friesian Cattle J Dairy Sci, November 1, 2005; 88(11): 4111 - 4119. [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] |
||||
![]() |
E E Connor, D L Wood, T S Sonstegard, A F da Mota, G L Bennett, J L Williams, and A V Capuco Chromosomal mapping and quantitative analysis of estrogen-related receptor alpha-1, estrogen receptors alpha and beta and progesterone receptor in the bovine mammary gland J. Endocrinol., June 1, 2005; 185(3): 593 - 603. [Abstract] [Full Text] [PDF] |
||||
![]() |
Z. Madeja, T. Adamowicz, A. Chmurzynska, T. Jankowski, J. Melonek, M. Switonski, and T. Strabel Short Communication: Effect of Leptin Gene Polymorphisms on Breeding Value for Milk Production Traits J Dairy Sci, November 1, 2004; 87(11): 3925 - 3927. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |