|
|
||||||||

* Department of Dairy Science, University of Wisconsin, Madison 53706
Genetic Epidemiology Division, St. Jamess University Hospital, Leeds LS9 7TF, United Kingdom
1 Corresponding author: hkhatib{at}wisc.edu
| ABSTRACT |
|---|
|
|
|---|
Key Words: OPN PPARGC1A candidate gene single nucleotide polymorphism
| INTRODUCTION |
|---|
|
|
|---|
Ron et al. (2001) localized a QTL affecting protein percentage to a confidence interval of 4 cM in the region near 50 cM on chromosome 6. Recently, Olsen et al. (2005) positioned a QTL affecting milk production traits at an interval of 420 kb between the genes ABCG2 and LAP3 on bovine chromosome 6. This narrow region harbors only 6 genes including the osteopontin (OPN) gene. Based on the aforementioned QTL studies, several groups investigated possible associations between genes in this region and milk production traits in dairy cattle. Weikard et al. (2005) reported significant association between a single nucleotide polymorphism (SNP) in intron 9 of the peroxisome proliferator activated receptor gamma coactivator 1 alpha (PPARGC1A) gene (close to microsatellite BM143) and milk fat yield in the German Holstein population. Leonard et al. (2005) and Schnabel et al. (2005) reported an association between OPN and milk protein percentage in the North American Holstein population. Furthermore, Schnabel et al. (2005) presented several lines of evidence for a quantitative trait nucleotide, located upstream of OPN promoter region, causing variation in milk protein percentage. On the other hand, Cohen-Zinder et al. (2005) reported the identification of a causative mutation in ABCG2 affecting milk yield and composition and excluded the OPN mutation reported in the Schnabel study from causality. The results for these conflicting candidate genes for the QTL in the middle of chromosome 6 necessitate additional populations to be tested for these genes (de Koning, 2006). In this study, the association of OPN variants with milk composition was investigated in an independent Holstein population for the validation of results previously obtained in a different granddaughter design Holstein population. In addition, the association between PPARGC1A variants, also located in the middle of chromosome 6, and milk production traits was also investigated in the University of Wisconsin (UW) daughter design and in the Cooperative Dairy DNA Repository (CDDR) granddaughter design resource populations. The OPN and PPARGC1A are about 6 Mb apart, which is about 12 cM for this region of chromosome 6 (http://www.ncbi.nlm.nih.gov/mapview).
| MATERIALS AND METHODS |
|---|
|
|
|---|
In addition, semen samples from 12 Holstein sires and their 581 sons were obtained from the CDDR, which is maintained by the USDA Bovine Functional Genomics Laboratory (Beltsville, MD). Predicted transmitting ability data for sons in the CDDR and yield deviation for daughters in the UW resource populations were obtained for milk, protein, and fat yields (kg), protein and fat percentages, and SCS from the USDA Animal Improvement Programs Laboratory (Beltsville, MD).
OPN and PPARGC1A SNP
It was previously shown that a C/T polymorphism at position 8514 of OPN (GenBank accession number NW_255516) is associated with milk protein and fat percentages in 1,362 bulls obtained from 28 sire families from the CDDR population (Leonard et al., 2005). To validate this association, 891 cows from the UW resource population were genotyped for this polymorphism. The association of PPARGC1A SNP with milk production was investigated in the UW and the CDDR resource populations using 2 SNP reported by Weikard et al. (2005); SNP C/T at position 1892 and SNP A/C at position 3359. The SNP 1892 was genotyped only for the UW resource population; SNP 3359 was genotyped in both populations.
SNP Genotyping
Genomic DNA was extracted from blood samples using GFX Genomic Blood DNA Purification Kit (Amersham Biosciences, Piscataway, NJ) and from semen samples by standard methods using proteinase K and phenol/chloroform. Primer sequences and amplification conditions for OPN were described in Leonard et al. (2005). Primers 1892F 5'-CATAGCCGGCGCCCCAGGTAAGATGCACGTTGGC-3' and 1892R 5'-CTGGTACTCCTCGTAGCTGTC-3' were used to amplify 195 bp in intron 9 to genotype SNP C/T at position 1892 of PPARGC1A. Primers 3359F 5'-GCGAGCACGGTGTTACATTACTAAGGAGAGTTGGCTAG-3' and 3359R 5'-GTTGTGTTGCACTCAATGGAC-3' were used to amplify 357 bp in the 3' untranslated region (3' UTR) of PPARGC1A to genotype SNP A/C at position 3359. Primers 1892R and 3359R were used by Weikard et al. (2005). Amplification was performed in a 15-µL reaction volume, which included 25 ng of genomic DNA, 25 ng of each primer, 100 µM of each dNTP, 1.5 µL of 10x PCR buffer (Promega, Madison, WI), and 0.2 U of Taq DNA polymerase (Promega). The temperature cycles were as follows: 95°C for 5 min, followed by 32 cycles of 94°C for 45 s, 50°C for 45 s, 72°C for 45 s, and a final extension at 72°C for 7 min. For genotyping of SNP C/T at position 1892, PCR products were digested with the restriction enzyme HaeIII that distinguishes alleles C and T of the SNP. For genotyping of SNP A/C at position 3359, PCR products were digested with the restriction enzyme NheI that distinguishes alleles A and C. The digestion products were electrophoresed on a 2.5% agarose gel.
Statistical Analysis
The association of OPN and PPARGC1A SNP with production traits was studied in 2 Holstein populations (CDDR and UW). The CDDR phenotypic data referred to bulls PTA data for milk, protein and fat yield, protein and fat percentages, and SCS. The UW resource population data was composed of yield deviations of cows for the same traits.
CDDR Population.
Associations between OPN variants and the phenotypic traits recorded in the CDDR population were published by Leonard et al. (2005). Therefore, in the present work only the PPARGC1A locus was considered in this population. Because the CDDR data are PTA, effects of an additional positive allele in the sons of a sire are necessarily additive. Hence, for this data set an allele substitution model having the following structure was used for the analysis of each trait:
![]() |
where yij is the PTA value relative to son j of sire i, µ is a general constant (intercept), si is the fixed effect of sire i, ß is the regression coefficient representing half of the allele substitution effect (
/2), xij is the number of A alleles (0, 1, or 2) at the PPARGC1A locus of son j of sire i, and
ij is a residual term. Reliabilities of the sons PTA were incorporated as weights in the model to obtain weighted least squares estimates of the allele substitution effects.
UW Resource Population.
Because the UW resource data are individual cow records, allele effects are not necessarily additive. Hence, data relative to each trait in the UW resource population were analyzed using the following fixed effects model:
![]() |
where yijkl represents the yield deviation of daughter j of sire i;
is an effect associated with M. paratuberculosis infection status; dij is an indicator variable assuming the values 0 and 1 for noninfected and infected cows, respectively; ok is the effect of the OPN genotype (k = CC, CT, or TT); and pl is the effect of the PPARGC1A genotype (l = AA, AC, or CC). The remaining terms in the model are as previously defined. Both OPN and PPARGC1A were analyzed in the same model to control for variation introduced by each locus. All interaction terms were included in the error term. This was not done for CDDR data because only the PPARGC1A locus was considered in this population at this time. Additive genetic effect of each locus was estimated as half of the difference between the 2 homozygous groups, as (ôCC ôTT)/2 for OPN and (
AA
CC)/2 for PPARGC1A. The dominance effect was estimated as the difference between the heterozygous group and the average of the 2 homozygous groups in each locus. The analyses were implemented using the GLM procedure of SAS (SAS Institute, 1999).
| RESULTS |
|---|
|
|
|---|
|
|
| DISCUSSION |
|---|
|
|
|---|
The PPARGC1A gene was chosen for association test with milk fat traits because of several QTL studies that showed effects of PPARGC1A region on milk production and because of possible involvement of this gene in fat metabolism (Weikard et al., 2005). The authors reported a significant association between SNP 1892 in intron 9 and milk yield and milk fat traits in 434 bulls from the German Holstein population. The SNP 3359 in the 3' UTR showed some trend of association with fat traits, but it did not reach the significance level of 5%, so the authors attributed this association to the linkage disequilibrium between the 2 SNP.
In contrast to Weikard et al. (2005), in the current study a significant association was found between SNP 3359 and protein percentage and milk yield, whereas SNP 1892 did not show significant association with any of the examined traits. These conflicting findings between results from this study and results from Weikard et al. (2005) could be explained by the 2 different genetic backgrounds of the studied populations or by the different number of samples included in each study. In the current study, PPARGC1A was genotyped in 2 independent populations with more than 1,400 individuals, whereas in the study by Weikard et al. (2005) 434 bulls were genotyped. Also, different associations of different SNP in the same gene in 2 or more populations could be a result of linkage disequilibrium of these SNP with the causative mutation in the same gene or in other genes. The SNP showing significant association with milk production traits would provide an excellent opportunity for marker-assisted selection (MAS) programs in dairy cattle. The aim of MAS is to substitute selection at the DNA level for selection on the basis of phenotype (Soller, 1994). Also, MAS can increase allele frequencies of the favorable alleles, as has been shown for the calpain and calpastatin genes in beef cattle (Pollak, 2005). However, for MAS to be efficient in genetic improvement programs, it is assumed that the effects of genes used in MAS are known and are consistent in different genetic backgrounds and environments (Dekkers and Hospital, 2002). The inconsistent association results between our data and Weikards necessitate further investigation prior to applications in MAS programs.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication December 4, 2006. Accepted for publication February 4, 2007.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
H. Khatib, W. Huang, X. Wang, A. H. Tran, A. B. Bindrim, V. Schutzkus, R. L. Monson, and B. S. Yandell Single gene and gene interaction effects on fertilization and embryonic survival rates in cattle J Dairy Sci, May 1, 2009; 92(5): 2238 - 2247. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. J. Hausman, M. V. Dodson, K. Ajuwon, M. Azain, K. M. Barnes, L. L. Guan, Z. Jiang, S. P. Poulos, R. D. Sainz, S. Smith, et al. BOARD-INVITED REVIEW: The biology and regulation of preadipocytes and adipocytes in meat animals J Anim Sci, April 1, 2009; 87(4): 1218 - 1246. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Khatib, C. Maltecca, R. L. Monson, V. Schutzkus, X. Wang, and J. J. Rutledge The fibroblast growth factor 2 gene is associated with embryonic mortality in cattle J Anim Sci, September 1, 2008; 86(9): 2063 - 2067. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. P. P. Macciotta, M. Mele, G. Conte, A. Serra, M. Cassandro, R. Dal Zotto, A. Cappio Borlino, G. Pagnacco, and P. Secchiari Association Between a Polymorphism at the Stearoyl CoA Desaturase Locus and Milk Production Traits in Italian Holsteins J Dairy Sci, August 1, 2008; 91(8): 3184 - 3189. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Wang, C. Maltecca, R. Tal-Stein, E. Lipkin, and H. Khatib Association of Bovine Fibroblast Growth Factor 2 (FGF2) Gene with Milk Fat and Productive Life: An Example of the Ability of the Candidate Pathway Strategy to Identify Quantitative Trait Genes J Dairy Sci, June 1, 2008; 91(6): 2475 - 2480. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Khatib, R. L. Monson, V. Schutzkus, D. M. Kohl, G. J. M. Rosa, and J. J. Rutledge Mutations in the STAT5A Gene Are Associated with Embryonic Survival and Milk Composition in Cattle J Dairy Sci, February 1, 2008; 91(2): 784 - 793. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |