1Institute of Biotechnology, Shiraz University, Shiraz, Iran
2Department of Genetics and Animal Breeding, Faculty of Animal Science and Fishery, Sari Agricultural Sciences and Natural Resources University, Sari, Iran
3Department of Animal Science, Faculty of Agriculture, Shiraz University, Shiraz, Iran
4Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, Kerman, Iran
Receive Date: 10 February 2016,
Revise Date: 11 June 2016,
Accept Date: 15 June 2016
The aim of this study was to investigate effects of DGAT1, OPN and PPARGC1A candidate genes on milk production traits in Iranian Holstein cattle. Several papers have studied single nucleotide polymorphisms (SNPs) and their association with economic traits in dairy cows, but the combined effect of these genes has not been examined in Iranian Holstein cattle population. Blood samples were collected from 398 registered Holstein cows. Total DNA was extracted using the salting out protocol. The PCR-RFLP technique was used for SNPs genotyping. The largest genotype frequency was estimated as 0.65 for PPARGC1A (c.1892)CT and the least frequency was estimated as 0.09 for DGAT1KK genotype. The allele frequencies were in the range 0.36 to 0.64 for PPARGC1A (c.3359)A and Calleles,respectively. The allelic substitution effects were estimated using a multiple regression model. The effects of allelic substitution for DGAT1K and PPARGC1A (c.1892)T were significant on estimated breeding values for fat percentage (EBVFP) (P<0.01). In addition, the results of multivariate analysis indicated the significant effect of DGAT1 and PPARGC1A (c.1892) on EBVFP (P<0.05). However, there were no association between OPN and PPARGC1A (c.3359) polymorphisms and studied traits.
Candidate gene approach and whole genome scans are two main strategies for QTL identification (Andersson, 2001). The candidate gene approach studies the relationship between the traits and known genes that may be associated with the physiological pathways underlying the trait (Liu et al. 2008). This approach has been successful to some extent. For example, several studies have identified QTLs for milk composition on chromosomes 6 and 14 (Riquet et al. 1999;Farnir et al. 2002; Olsen et al. 2005). The economic traits are polygenic. It means that they are controlled by many loci. Several studies indicated genetic variation in milk production traits cannot be explained by few candidate genes (Kaupe et al. 2007). Therefore, the effects of all candidate loci should be explored together in the same statistical model diacyl glycerol acyltransferase 1(DGAT1) is located near the centromeric region of Bos taurus autosome 14 (BTA14). The first evidence for the effect of DGAT1 variation on milk yield and composition in Holstein cattle was reported by Grisart et al. (2002). DGAT1 is considered as the key enzyme in controlling the synthesis rate of triglycerides in adipocytes. A non-conservative K232A substitution (conservation of alanine to lysine) in DGAT1 was associated with milk production and compositions in Holstein cattle (Thaller et al. 2003). Spelman et al. (2002) and Banos et al. (2008) reported that the K232A substitution in exon 8 of the DGAT1 gene was associated with increasing of milk fat yield and decreasing of milk production and protein yield. Some studies showed that there are significant associations between DGAT1 and milk, fat yield and protein yield. Bovine chromosome six (BTA6) harbors at least six QTLs influencing milk production traits of dairy cattle. The osteopontin(OPN) and peroxisome proliferator activated receptor gamma co-activator 1 Alpha(PPARGC1A) are about 6 Mb apart, which is about 12 cM for this region of chromosome 6 (Olsen etal. 2005). OPN is a strong functional candidate for milk production and it is a highly phosphorylated glycoprotein (Leonard et al. 2005). Schnabel et al. (2005) reported an association between OPN and milk protein percentage in the North American Holstein population. PPARGC1A has main role in fat and glucose metabolism and plays a critical role in the activation of nuclear hormone receptors and transcription factors regulating energy homeostasis (Liang and Ward, 2006; Kowalewska-Luczak et al. 2010). Structure of PPARGC1A gene is made from 13 exons and expressed at different levels in a great number of tissues (Liang and Ward, 2006). Khatib et al. (2007) showed significant associations between PPARGC1A (c.3359)gene, milk yield, milk protein percentage, and somatic cell score in the North American Holstein population. The aim of this study was to investigate the joint effects of OPN, PPARGC1A and DGAT1 candidate genes on milk production traits in Iranian Holstein cattle population.
MATERIALS AND METHODS
Animals and traits
Totally 398 blood samples were collected from Holstein-Friesian cows of Iran, which were distributed in ten dairy herds in two provinces of Iran. The cows were under official milk recording of AnimalBreedingCenter (Karaj-Iran). Finally 372 records for estimated breeding values for milk production adjusted for 305 days (EBVM), fat yield (EBVF) (kg) and fat percentage (EBVFP) were obtained from the AnimalBreedingCenter for analyzing association between genotypes and economics traits. The EBVs were estimated by random regression test day model.
DNA extraction, PCR amplification and SNPs genotyping
DNA extractions were performed using standard salting out protocol (Miller et al. 1988). PCR reactions were performed using standard PCR (Thermo cycler, Biometra, Germany). More details about primers are shown in Table 1 (Kaupe et al. 2004; Weikard et al. 2005; Khatib et al. 2007). PCR reaction for DGAT1 (GenBank: EU077528), OPN (GenBank: NW_255516) and PPARGC1A (GenBank: AY321517) loci were performed in a 25 µL volume using 100 ng genomic DNA, PCR buffer (1X), 1.5 mM MgCl2, 0.2 mM dNTPs, 0.6 pmol of each primer and Taq polymerase enzyme (2U). All accession numbers are available in NCBI site. For DGAT1 gene, the addition of DMSO to the PCR reactions allowed an equal amplification of both alleles. The annealing temperature for DGAT1, OPN and PPARGC1A are considered as 60, 53 and 55 centigrade and finally 411, 290, 195 and 357 base pairs fragments were amplified for DGAT1, OPN, PPARGC1A (c.1892) and PPARGC1A (c.3359). The PCR products (5 µL) were digested using 2 units of the restriction enzymes (FERMENTAS, Lithuania) and separated on a 2% agarose gel. The gels were stained with ethidium bromide and visualized under UV light. Finally, SNPs were genotyped by PCR-RFLP technique. Table 2 shows more detail about restriction enzyme conditions.
Gene and genotype frequencies
The population genetic parameters including gene and genotypic frequencies, Hardy-Weinberg equilibrium (Chi-square test), Indices of genetic diversity in population (Nei (H) and Shannon (I)) were estimated using the PopGene software version 3.1d (Nei, 1977).
Joint analysis of DGAT1, OPN and PPARGC1A variants
The effects of genotypes were tested using the GLM procedure (Pillais trace test) of SPSS (2010) implementing the following fixed model:
yijkmn= µ + Pi + Aj + Dk + Om + eijkmn
y: observation for each trait.
µ: overall mean.
Pi, Aj, Dk andOm: fixed effects of genotypes of PPARGC1A (c.1892), PPARGC1A (c.3359), DGAT1 and OPN genes. eijkmn: residual effect.
The mean comparisons were performed using the Tukey test for significant genotypes.
Table 1 Primer sequence for PCR reaction
Table 2 Restriction enzymes and the digestion conditions
The effects of allele substitutions on milk production traits were tested using the following multiple linear regression models (Knott et al. 1996):
yijkl= µ + bixi + bjxj + bkxk + blxl + eijkl
y: observation for EBVM, EBVF and EBVFP traits.
µ: overall mean.
bi, bj, bk, bl: regression coefficients representing the allelic substitutions for (DGATK, OPNT, PPARGC1A (c.3359)A, PPARGC1A (c.1892)T.
xi, xj, xk, xl: indicator variables for genotypes of DGAT1, OPN, PPARGC1A (c.3359), PPARGC1A (c.1892)loci.
eijkl: residual effect.
RESULTS AND DISCUSSION
The most extreme genotypes frequencies were estimated as 0.65 and 0.09 for PPARGC1A (c.1892)CT and DGAT1KK loci, respectively. Similar results were obtained about genotype frequencies in Holstein cattle population by Khatib et al. (2007), Thaller et al. (2003) and Komisarek and Dorynek (2009). In addition, the most and the least allele frequencies were calculated as 0.64 and 0.36 for A and C alleles of PPARGC1A (c.3359). The joint testing of all of the candidate loci in the population under study indicated significant deviation from the Hardy-Weinberg equilibrium (P<0.05), which seems expectable due to the long time selection for milk production. The genetic diversity indices showed that the population has desirable genetic variation. More detail about values and frequencies are shown in Table 3. Figure 1 gives more information about the results of digestion. Similar results were obtained using the multivariate analysis and between subject test Table 4. The results indicated that the PPARGC1A (c.1892) and DGAT1polymorphismshad significant association with EBVFP. The summaries of statistical analysis are illustrated in Tables 5 and 6. The results obtained are supported from other studies. Weikard et al. (2005) reported significant association between SNP in intron 9 of the PPARGC1A (c.1892) gene and fat yield, which means that the PPARGC1A gene might be involved in genetic variation underlying the QTL for milk fat synthesis on BTA6. Schennink et al. (2009) showed that two SNPs in PPARGC1A (c.3359 and c.1892) had significant effects on fat yield. However, we found that the effect of PPARGC1A (c.3359) and OPN polymorphism were not significant, but different results reported by Zhang et al. (1998) and Mosig et al. (2001), who identified candidate gene affecting milk production traits close to OPN location.
Table 3 Summery of frequencies, H-W equilibrium and genetic variation indices
* (P<0.05) and ** (P<0.01).
Figure 1 Electrophoretic separation of DGAT1 (A), OPN (B), PPARGC1A (c.1892)(C) and PPARGC1A (c.3359)(D)genes PCR products
Figures a, c and d have a common ladder (PUC mix 8)
Figure b: ladder (gene ruler DNA ladder)
Table 4 The results of multivariate analysis (F-statistics) for EBVFP
Table 5 The analysis results of between subjects effects (F-values)
Table 6 The results of mean comparisons for significant genes
The means within the same column with at least one common letter, do not have significant difference (P>0.05).
Table 7 The effect of allele substitution for candidate genes
In addition, Leonard et al. (2005) showed significant association of OPN gene with milk protein percentage. Therefore, we suggest the further studies need to clarify the association between OPN, PPARGC1A (c.3359) and milk production traits in other populations. There are several possible reasons for different results of studies, including differences in allele frequency, the statistical models used to undertake the association analysis and genetic background of the animals in the study (Berry et al. 2010) and environmental circumstances where the animals were producing. The association between DGAT1 gene and EBVFP was significant which may be due to critical role of DGAT1 gene in the synthesis rate of triglycerides (Grisart et al. 2002). Kadlecova et al. (2014) reported significant association between DGAT1 genotypes and fat percentage in primiparous Holstein cows. Anton et al. (2012) indicated the significant effects of the DGAT1 K232A polymorphism on milk yield, fat and protein percentage, as well. In addition, Fontanesi et al. (2015) illustrated that DGAT1 polymorphism was highly associated with fat yield and fat percentage in Reggiania dairy cows (local breed in north of Italy). The results of mean comparisons illustrated that the genotypes of DGAT1KK and PPARGC1A(c.1892)TT had highest EBVFP. More results of mean comparisons are shown in Table 6. Table 7 gives the additive effects (allele substitutions) of the alleles. The result indicated that and DGATK allele increased the EBVFP by +0.04 ± 0.01. Some standard errors were estimated more than their allele substitution effect. It can be due to low number of data, which are used in this study. The results of allele substitutions were confirmed by other studies. Winter et al. (2002) and Strzałkowska et al. (2005) showed that the DGAT1K allele has a positive effect on milk fat content in different cattle breeds. Naslund et al. (2008) reported that DGAT1K variant was associated with an increase in milk fat and protein percentages but decrease milk yield compared with the DGATA variant. Similar results showed that the DGATK allele exceeds of DGATA allele, by (+0.34) percentage unit in fat (Grisart et al. 2002). The DGAT1K allele increases milk fat yield, whereas the DGAT1A allele increases both milk and protein yield (Kaupe et al. 2007; Thaller et al. 2003). According to our finding PPARGC1A (c.1892)T allele decreased the EBVFP by -0.37 ± 0.01. In addition, an association of the PPARGC1A (c.1892)T allele with higher fat yield has been suggested in German Holsteins (Weikard et al. 2005). Alim et al. (2012) indicated that PPARGC1A (c.1892)T allele increased protein yield and protein concentration but there was no association between PPARGC1A (c.1892)T allele and fat yield (%, kg).
Milk and its products are regarded as the most important nutritional resource, meeting the energy requirements and offering high quality protein and various vitamins and minerals. Earlier, most genetic improving programs of agriculturally important livestock population have been carried out through complete phenotypic and pedigree information. However, applying molecular genetic information in breeding stock may lead to a better understanding of quantitative traits. Briefly, the results show that there is significant association between PPARGC1A, DGAT1 and EBVFP trait. Generally, detection and estimation of associations of identified genes and genetic markers with economic traits are the basis of a successful application of marker-assisted selection (MAS) in breeding programs. The MAS strategies can be used for pre-selection of young bulls prior to progeny test.
This study was financially supported by Agricultural Biotechnology Research Institute of Iran.
Alim M.A., Fan Y., Xie Y., Wu X., Sun D., Zhang Y., Zhang S., Zhang Y., Zhang Q. and Lin L. (2012). Single nucleotide polymorphism (SNP) in PPARGC1A gene associates milk production traits in Chinese Holstein cattle. Pakistan Vet. J.32, 609-612.
Andersson L. (2001). Genetic dissection of phenotypic diversity in farm animals. Nat. Genet. 2, 130-138.
Anton I., Kovacs K., Hollo G., Farkas V., Szabo F., Egerszegi I., Ratky J., Zsolnai A. and Brüssow K.P. (2012). Effect of DGAT1, leptin and TG gene polymorphisms on some milk production traits in different dairy cattle breeds in Hungary. Arch. Tierzucht.55(4), 307-314.
Banos G., Woolliams J.A., Woodward B.W., Forbes A.B. and Coffey M.P. (2008). Impact of single nucleotide polymorphisms in leptin, leptin receptor, growth hormone receptor, and diacylglycerol acyltransferase (DGAT1) gene loci on milk production, feed and body energy traits of UK dairy cows. J. Dairy Sci.91, 3190-3200.
Berry D.P., Howard D., Boyle P.O., Water S., Kearney J.F. and McCabe M. (2010). Associations between the K232A polymorphism in the diacylglycerol-O-transferase 1 (DGAT1) gene and performance in Irish Holstein-Friesian dairy cattle. Irish J. Agric. Food Res. 49, 1-9.
Farnir F., Grisart B., Coppieters W., Riquet J., Berzi P., Cambisano N., Karim L., Mni M., Simon P. and Wagenaar D. (2002). Simultaneous mining of linkage and LD to fine map QTL in outbred half-sib pedigrees: revisiting the location of a QTL with major effect on milk production on bovine chromosome 14. Genetics. 161, 275-287.
Fontanesi L., Scotti L., Samore A.B., Bagnato A. and Russo V. (2015). Association of 20 candidate gene markers with milk production and composition traits in sires of reggiana breed, a local dairy cattle population. Livest.Sci. 176, 14-21.
Grisart B., Coppieters W., Farnir F., Karim L., Ford C., Berzi P., Cambisano N., Mni M., Reid S., Simon P., Spelman R., Georges M. and Snell R. (2002). Positional candidate cloning of a QTL in dairy cattle: identification of a missense mutation in the bovine DGAT1 gene with major effect on milk yield and composition. Genome. Res.12, 222-231.
Kadlecova V., Nemeckova D., Katerina Jecminkova K. and Stadni K.L. (2014). Association of bovine DGAT1 and leptin genes polymorphism with milk production traits and energy balance indicators in primiparous Holstein cows. Mljekarstvo. 64(1), 19-26.
Kaupe B., Winter A., Fries R. and Erhardt G. (2004). DGAT1 polymorphism in Bos indicus and Bos taursus cattle breeds.J. Dairy Res. 71, 182-187.
Kaupe B., Brandt H., Prinzenberg E.M. and Erhardt G. (2007). Joint analysis of the influence of CYP11B1 and DGAT1 genetic variation on milk production, somatic cell score, conformation, reproduction and productive lifespan in German Holstein cattle. J. Anim. Sci. 85, 11-21.
Khatib H., Zaitoun I., Wiebelhaus-Finger J., Chang Y.M. and Rosa G.J.M. (2007). The association of bovine PPARGC1A and OPN genes with milk composition in two independent Holstein cattle populations. J. Dairy Sci. 90, 2966-2970.
Knott S.A., Elsen J.M. and Haley C.S. (1996). Methods for multiple-marker mapping of quantitative trait loci in half-sib populations. Theor. Appl. Genet.93, 71-80.
Komisarek J. and Dorynek Z. (2009). Effect of ABCG2, PPARGC1A, OLR1 and SCD1 gene polymorphism on estimated breeding values for functional and production traits in Polish Holstein-Friesian bulls. J. Appl. Genet.50, 125-132.
Kowalewska-LuczakI., Kulig H. and Kmiec M. (2010). Associations between the bovine PPARGC1A gene and milk production traits. Czech J. Anim. Sci.55, 195-199.
Leonard S., Khatib H., Schutzkus V., Chang Y.M. and Maltecca C. (2005). Effects of the osteopontin gene variants on milk production traits in dairy cattle. J. Dairy Sci.88, 4083-4086.
Liang H. and Ward W. (2006). PGC-1α: a key regulator of energy metabolism. Adv. Physiol. Educ.30, 145-151.
Liu X., Zhang H., li H., li N., Zhang Y., Zhang Q., Wang S., Wang Q. and Wang H. (2008). Fine mapping quantitative trait loci for body weight and abdominal fat traits: effects of marker density and sample size. Poult. Sci.87, 1314-1319.
Miller S.A., Dykes D.D. and Polesky H.F. (1988). A simple salting-out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res.16, 1215-1219.
Mosig M.O., Lipkin E., Khutoreskaya G., Tchourzyna E., Soller M. and Friedmann A. (2001). A whole genome scan for quantitative trait loci affecting milk protein percentage in Israeli-Holstein cattle, by means of selective milk DNA pooling in a daughter design, using an adjusted false discovery rate criterion. Genetics.157, 1683-1698.
Naslund J., Fikse W.F., Pielberg G.R. and Lunden A. (2008). Frequency and effect of the bovine acyl-CoA: diacylglycerolAcyltransferase 1 (DGAT1) K232A polymorphism in Swedish dairy cattle. J. Dairy Sci.91, 2127-2134.
Nei M. (1977). F-statistics and analysis of gene diversity in subdivided populations. Ann. Hum. Genet.41, 225-233.
Olsen H.G., Lien S., Gautier M., Nilsen H., Roseth A., Berg P.R., Sundsaasen K.K., Svendsen M. and Meuwissen T.H. (2005). Mapping of a milk production quantitative trait locus to a 420-kb region on bovine chromosome 6. Genetics.169, 275-283.
Riquet J., Coppieters W., Cambisano N., Arranz J.J., Berzi P., Davis S.K., Grisart B., Farnir F., Karim L., Mni M., Simon P., Taylor J.F., Vanmanshoven P., Wagenaar D., Womack J.E. and Georges M. (1999). Fine-mapping of quantitative trait loci by identity by descent in outbred populations: application to milk production in dairy cattle. Proc. Nat. Acad. Sci. USA. 16, 9252-9257.
Schennink A., Bovenhuis H., Leon-Kloosterziel K.M., Arendonk A.M. and Visker M.H.P. (2009). Effect of polymorphisms in the FASN, OLR1, PPARGC1A, PRL and STAT5A genes on bovine milk-fat composition. Anim. Genet. 40, 909-916.
Schnabel R.D., Kim J.J., Ashwell M.S., Sonstegard T.S., Van Tassell C.P., Connor E.E. and Taylor J.F. (2005). Fine-mapping milk production quantitative trait loci on BTA6: analysis of the bovine osteopontin gene. Proc. Nat. Acad. Sci. USA. 102, 6896-6901.
Spelman R.J., Ford C.A., Mcelhinney P., Gregory G.C. and Snell R.G. (2002). Characterization of the DGAT1 gene in the New Zealand dairy population. J. Dairy Sci.85, 3514-3517.
SPSS Inc. (2010). Statistical Package for Social Sciences Study. SPSS for Windows, Version 20. Chicago SPSS Inc.
Strzałkowska N., Siadkowska E., Słoniewski K., Krzyżewski J. and Zwierzchowski L. (2005). Effect of the DGAT1 gene polymorphism on milk production traits in Black-and-White (Friesian) cows. Anim.Sci. Pap. Rep. 23, 189-197.
Thaller G., Kramer W., Winter A., Kaupe B., Erhardt G. and Fries R. (2003). Effects of DGAT1 variants on milk production traits in German cattle breeds. J. Anim. Sci.81, 1911-1918.
Weikard R., Kuhn C., Goldammer T., Freyer G. and Schwerin M. (2005). The bovine PPARGC1A gene: molecular characterization and association of an SNP with variation of milk fat synthesis. Physiol. Genomics. 21, 1-13.
Winter A., Krämer W., Werner F.A., Kollers S., Kata S., Durstewitz G., Buitkamp J., Womack J.E., Thaller G. and Fries R. (2002). Association of a lysine-232/alanine polymorphism in a bovine gene encoding acyl-CoA: diacylglycerol acyltransferase (DGAT1) with variation at a quantitative trait locus for milk fat content. Proc. Nat. Acad. Sci. USA.99, 9300-9305.
Zhang Q., Boichard D., Hoeschele I., Ernst C., Eggen A., Murkve B., Pfistergenskow M., Witte L.A., Grignola F.E., Uimari P., Thaller G. and Bishop M.D. (1998). Mapping QTL for milk production and health of dairy cattle in a large outbred pedigree. Genetics. 149, 1959-1973.