Polymorphisms in the promoter region of the Chinese bovine PPARGC1A gene.
The peroxisome proliferator-activated receptor gamma coactivator-1 alpha (PPARGC1A, also known as PGC-1a) protein, encoded by the PPARGC1A gene, plays an important role in energy homeostasis. PGC-1a is a metabolic switch, which convergently regulate metabolic pathways through its pleiotropic interactions with nuclear receptors (NRs) and other non-NR transcription factors. PGC-1a transcriptionally activates a complex pathway of mitochondrial biogenesis, lipid and glucose metabolism (Lin et al., 2005; Puigserver, 2005). The expression of PGC-1a is enriched in tissues where mitochondria are abundant and oxidative metabolism is active, such as brown adipose tissue (BAT), the cardiac and skeletal muscles (Puigserver et al., 1998).
Independent studies in humans have identified a common hPPARGC1A polymorphism, Gly482Ser, associated with type 2 diabetes (Barroso et al., 2006; Bhat et al., 2008; Povel et al., 2010; Geloneze et al., 2012). The most common polymorphism c.1892-19T>C in bovine PPARGC1A intron 9 was studied for the associations with milk production and milk-fat composition traits (Weikard et al., 2005; Khatib et al., 2007; Komisarek and Dorynek, 2009; Schennink et al., 2009; Kowalewska-Luczak et al., 2011). Several other polymorphisms in the bovine PPARGC1A gene have been identified and studied for the association with meat quality and carcass phenotypes (Soria et al., 2009; Ryu et al., 2012).
As an important component of the world's bovine population, the Chinese native population possesses very abundant genetic resources. However, up to now, there was no information about the genetic variation of PPARGC1A gene promoter in Chinese indigenous bovine breeds. Therefore, a scan of the genetic variations within the PPARGC1A gene promoter in Chinese native bovines, should provide important information on bovine genetic resources that could be useful for selection and breeding through marker assisted selection (MAS).
MATERIAL AND METHODS
Bovine populations, genomic DNA isolation and pooling
Genomic DNA samples were obtained from 808 bovines belonging to three cattle breeds and yaks. Three Chinese indigenous cattle breeds, Jiaxian (JX, n = 148), Nanyang (NY, n = 278), Qinchuan (QC, n = 325) were investigated. The Jiaxian cattle were from the breeding farm of Jiaxian cattle (Jiaxian county, Henan Province, China); the Nanyang cattle were from the breeding farm of Nanyang cattle (Nanyang city, Henan Province, China); the Qinchuan cattle were from the reserved farm (Fufeng county, Shaanxi Province, China), the breeding and seed centre of Qinchuan cattle (Yangling District, Shaanxi Province, China). The yaks (n = 57) were from the breeding farm of yaks (Datong county, Qinghai province, China). DNA was isolated from 1 ml whole blood samples according to Sambrook and Russell (2001). DNA of 50 to 100 randomly chosen individuals from the same breed were mixed to one DNA pool (Bansal et al., 2002).
Six primer pairs were designed with the Primer V5.0 software based on NCBI database (GenBank accession number AC_000163) to amplify the promoter region of the bovine PPARGC1A gene. Detailed information about oligonucleotide primers, amplicon sizes and corresponding annealing temperatures is given in Table 1.
Single nucleotide polymorphisms (SNPs) detection
The genomic DNA pool samples were used as templates. Each amplification reaction was carried out in a total volume of 25 [micro]l in the presence of 50 ng bovine genomic DNA pool as template, 0.5 [micro]M of each primer, 1 x buffer (including 1.5 mM [MgCl.sub.2]), 200 [micro]M dNTPs (dATP, dTTP, dCTP and dGTP), and 0.625 U Taq DNA polymerase (MBI, Vilnius, Lithuania). The cycling protocol was 4 min at 95[degrees]C, 35 cycles of denaturing at 94[degrees]C for 30 s, annealing at the optimized temperature for 30 s, extending at 72[degrees]C for 30 s, with a final extension at 72[degrees]C for 10 min. All PCR products were then sequenced using ABI 3730 sequencer (ABI, Foster City, CA, USA), in both the forward and reverse directions, and the sequences were imported into the Megalign software module of DNASTAR (version 4.0) and BioXM software (version 2.6) to search for SNPs.
The genomic DNA samples were used as templates to individually and specifically amplify the promoter region of PPARGC1A. SSCP analysis of PCR fragments obtained above was performed on electrophoresis units (Bio-Rad, Amersham), coupled with a refrigerated system. Aliquots of 5 [micro]l PCR products were mixed with 5 denaturing solution (95% formamide, 25 mM EDTA, 0.025% xylene-cyanole and 0.025% bromophenol blue), heated for 10 min at 98[degrees]C and chilled on ice. Denatured DNA was subjected to PAGE (80 x 73 x 0.75 mm) in 1 x TBE buffer and constant voltage (200 V) for 2.5 to 3.0 h. After electrophoresis, gels were stained with 0.1% silver nitrate and visualized with 2% NaOH solution (supplied with 0.1% formaldehyde) according to Zhang et al. (2007). After the polymorphisms were detected, three PCR products of each different electrophoresis pattern were sequenced in both directions in an ABI 3730 sequencer (ABI, Foster City, CA, USA) and analyzed with BioXM software (version 2.6).
Allele frequencies were determined for each breed by direct counting. The program MatInspector (http://www.genomatrix.de) was performed to identify any known transcription factors that may potentially bind to polymorphic promoter sequences in the TRANSFAC database.
RESULTS AND DISCUSSION
In the present study, genomic DNA of all bovine breeds was successfully amplified using primer pairs for the promoter region of the bovine PPARGC1A gene. A total of 6 SNPs and one 4 bp insertion variation in the promoter region of bovine PPARGC1A gene were identified. In the PCR product primer1, the SNP -259 T>A and -301_298insCTTT were found; in primer3, -915 A>G; in primer4, -1175 T>G; in primer5, -1590 C>T, -1665 C>T and -1690 G>A. All these polymorphisms are numbered relative to the adenine of the translation start codon. The sequencing results and SNPs found are shown in Figure 1. The nature and the distribution of PPARGC1A promoter variations are shown in Table 2.
The SNPs -259 T>A, -1175 T>G, -1590 C>T and -1665 C>T occurred only in the Yak population. Conversely, the SNP -915 A>G occurred only in cattle breeds (JX, NY and QC), while -301_-298insCTTT and SNP -1690 G>A occurred both in Yak and cattle populations. The genetic differences of PPARGC1A promoter in the present study give support to the notion that the yak is an independent genus, Poephagus. This notion was also supported by the reports from Zhang (Zhang et al., 2010; Zhang et al., 2011). The discrepancy in different populations was probably caused by the following two possible reasons. First, the yak was an ancient animal species with many wild bovine characteristics, which did not undergo much artificial selective breeding. Meanwhile, the Jiaxian, Nanyang and Qinchuan cattle, three of the most representative indigenous bovine (Bos taurus) breeds in China, are draft and beef dual-purpose cattle breeds which have been selected for beef production in the past decades. Second, yaks must survive in harsh environments (low temperature, hypoxia with high altitude). So, some variations may only exist in yaks as a result of adaptation their environment.
Polymorphisms in promoter region may affect gene transcription through changing the affinity between the transcription factors and promoter sequences. To investigate the potential effect that some of the polymorphisms could have on gene transcription, the program MatInspector (http://www.genomatrix.de) was performed in the TRANSFAC database to identify any known transcription factors that may potentially bind to polymorphic promoter sequences. The -301_-298insCTTT makes the recognition sequence of sex-determining region Y (SRY) repeated. The SRY locus on the mammalian Y chromosome is the developmental switch responsible for testis determination. Inconsistent with this important function, SRY has been demonstrated to modulate the catecholamine pathway, so it should have functional consequences in the central and peripheral nervous system (Turner et al., 2011). The -915 A>G makes the recognition sequence of myeloid-specific zinc finger-1 (MZF-1) disappeared. The zinc finger gene MZF-1 is preferentially expressed in primitive hematopoietic cells and plays an important role in regulating myelopoiesis (Hromas et al., 1995). In this polymorphic site, the yak population was in homozygous AA status, which preserves the recognition sequence of MZF-1. These traits are crucial for yaks to survive in severe environments. The -1665 C>T is in the core recognition sequence of octamer factor 1 (Oct-1). It is now known that Oct-1 controls not only housekeeping but also numerous tissue-specific genes. The latter include the genes for interleukins (IL) 2, 8, 3, 5; the granulocyte-macrophagal colony-stimulating factor, immunoglobulins a; the light and heavy chains of immunoglobulins; the endocrine-associated Pit-1 gene; the genes for gonadoliberin, prolactin, the thyroid transcription factor 1, and thyrotropin (Sytina and Pankratova, 2003). It is expected that these polymorphisms regulate PPARGC1A gene transcription and might have consequences at a regulatory level, but there are no data available to suggest the interaction between the variations and the PPARGC1A promoter activity until now. The potential biological function of the mutation in bovine PPARGC1A promoter remains to be further studied.
This study was supported by the National Natural Science Foundation of China (No. 30972080, 30901023), Program of National Beef Cattle Industrial Technology System (No. CARS-38), Agricultural Science and Technology Innovation Projects of Shaanxi Province (No.2012NKC01-13).
Bansal, A., D. van den Boom, S. Kammerer, C. Honisch, G. Adam, C. R. Cantor, P. Kleyn and A. Braun. 2002. Association testing by DNA pooling: An effective initial screen. Proc. Natl. Acad. Sci. USA. 99:16871-16874.
Barroso, I., J. Luan, M. S. Sandhu, P. W. Franks, V. Crowley, A. J. Schafer, S. O'Rahilly and N. J. Wareham. 2006. Meta-analysis of the gly482ser variant in ppargc1a in type 2 diabetes and related phenotypes. Diabetologia 49:501-505.
Bhat, A., A. Koul, E. Rai, S. Sharma, M. K. Dhar and R. N. Bamezai. 2007. Pgc-1 alpha thr394thr and gly482ser variants are significantly associated with t2dm in two north indian populations: A replicate case-control study. Hum. Genet. 121:609-614.
Geloneze, S. R., B. Geloneze, J. Morari, J. R. Matos-Souza, M. M. Lima, E. A. Chaim, J. C. Pareja and L. A. Velloso. 2012. Pgc1 alpha gene gly482ser polymorphism predicts improved metabolic, inflammatory and vascular outcomes following bariatric surgery. Int. J. Obes. 36:363-368.
Hromas, R., J. Morris, K. Cornetta, D. Berebitsky, A. Davidson, M. Sha, G. Sledge and F. Rauscher. 1995. Aberrant expression of the myeloid zinc finger gene, mzf-1, is oncogenic. Cancer Res. 55:3610-3614.
Khatib, H., I. Zaitoun, J. Wiebelhaus-Finger, Y. M. Chang and G. J. Rosa. 2007. The association of bovine ppargc1a and opn genes with milk composition in two independent holstein cattle populations. J. Dairy Sci. 90:2966-2970.
Komisarek, J. and Z. Dorynek. 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-Luczak, I., H. Kulig and M. Kosobucki. 2011. Effect of polymorphism in the ppargc1a gene on the milk production traits of polish holstein-friesian cows. Tieraerztl. Umsch. 66:147-150.
Lin, J. D., C. Handschin and B. M. Spiegelman. 2005. Metabolic control through the pgc-1 family of transcription coactivators. Cell Metab. 1:361-370.
Povel, C. M., E. J. Feskens, S. Imholz, E. E. Blaak, J. M. Boer and M. E. Dolle. 2010. Glucose levels and genetic variants across transcriptional pathways: Interaction effects with bmi. Int. J. Obes. 34:840-845.
Puigserver, P. 2005. Tissue-specific regulation of metabolic pathways through the transcriptional coactivator pgc1-alpha. Int. J. Obes. 29:S5-S9.
Puigserver, P., Z. Wu, C. W. Park, R. Graves, M. Wright and B. M. Spiegelman. 1998. A cold-inducible coactivator of nuclear receptors linked to adaptive thermogenesis. Cell 92:829-839.
Ryu, J., Y. Kim, C. Kim, J. Kim and C. Lee. 2012. Association of bovine carcass phenotypes with genes in an adaptive thermogenesis pathway. Mol. Biol. Rep. 39:1441-1445.
Sambrook, J. and D. W. Russell. 2001. Molecular cloning: A laboratory manual (3rd ed.). Cold Spring Harbor Laboratory Press, New York.
Schennink, A., H. Bovenhuis, K. M. Leon-Kloosterziel, J. A. van Arendonk and M. H. Visker. 2009. Effect of polymorphisms in the fasn, olr1, ppargc1a, prl and stat5a genes on bovine milkfat composition. Anim. Genet. 40:909-916.
Soria, L. A., P. M. Corva, A. Branda Sica, E. L. Villarreal, L. M. Melucci, C. A. Mezzadra, J. Papaleo Mazzucco, G. Fernandez Macedo, C. Silvestro, A. Schor and M. C. Miquel. 2009. Association of a novel polymorphism in the bovine ppargc1a gene with growth, slaughter and meat quality traits in brangus steers. Mol. Cell. Probes. 23:304-308.
Sytina, E. V. and E. V. Pankratova. 2003. Transcription factor oct1: Plasticity and multiplicity of functions. Mol. Biol. 37:637-648.
Turner, M. E., D. Ely, J. Prokop and A. Milsted. 2011. Sry, more than testis determination? Am. J. Physiol-Reg. Integr. Comp. Physiol. 301:R561-R571.
Weikard, R., C. Kuhn, T. Goldammer, G. Freyer and M. Schwerin. 2005. The bovine ppargc1a gene: Molecular characterization and association of an snp with variation of milk fat synthesis. Physiol. Genomics 21:1-13.
Zhang, C. L., Y. H. Wang, H. Chen, X. Y. Lan and C. Z. Lei. 2007. Enhance the efficiency of single-strand conformation polymorphism analysis by short polyacrylamide gel and modified silver staining. Anal. Biochem. 365:286-287.
Zhang, Q., H. Chen, S. Zhao, L. Zhang, L. Zhang and X. Wang. 2010. Polymorphisms in the promoter region of bovine prkab1 gene. Mol. Biol. Rep. 37:435-440.
Zhang, Q., S. Zhao, H. Chen, L. Zhang, L. Zhang, F. Li and X. Wang. 2011. Snp discovery and haplotype analysis in the bovine prkaa2 gene. Mol. Biol. Rep. 38:1551-1556.
M. J. Li, M. Liu, D. Liu, X. Y. Lan, C. Z. Lei, D. Y. Yang (1) and H. Chen *
College of Animal Science and Technology, Northwest A&F University, Shaanxi Key Laboratory of Molecular Biology for Agriculture, Yangling, Shaanxi 712100, China
* Corresponding Author: H. Chen. Tel: +86-29-87092102, Fax: +86-29-87092164, E-mail: email@example.com
(1) Department of Biology, Dezhou University, Dezhou, Shandong 253023, China.
Submitted Oct. 8, 2012; Accepted Nov. 19, 2012; Revised Jan. 8, 2013
Table 1. Primer information of bovine PPARGC1A gene Primer Primer sequence (5' to 3') Size (bp) P1 1F: TATTATATCGCCAGGGCTC 432 1R: TCGATGTCACTCCATACAGA P2 2F: GGTATTTTGTTGCTTGTTCTCC 417 2R: AATAACTAAAATCCCCTGCAATC P3 3F: TGCGTTAGCTCTCATTGTCTC 363 3R: AACAAGCAACAAAATACCATCC P4 4F: TGATAAAGGGAAGGACAAT 347 4R: TCGAGCCTGAGACAATGAG P5 5F: TTTTCCTGTTCTTTCTTGAT 352 5R: TCTATTGTCCTTCCCTTTAT P6 6F: TCAGTGAACACGAATAAAA 290 6R: GAAAGAACAGGAAAATACAA Primer Tm Note (ref. AC_000163) P1 49.1 -385 to 47 P2 48.0 -798 to -382 P3 50.5 -1143 to -781 P4 48.0 -1461 to -1115 P5 52.4 -1791 to -1440 P6 46.9 -2072 to -1783 PCR primers were designed according to the bovine PPARGC1A gene (GenBank accession No. AC 000163). Location of primers referred to the translation start site, where +1 corresponding to the adenine of the translation initiation codon ATG. Table 2. The SNPs identified in the promoter region of the bovine PPARGC1A gene in Chinese breeds and allele frequencies of the SNPs Polymorphism site Allelic variant frequency JX (n = 148) NY (n = 278) -259 T>A 0 0 -301_-298insCTTT 0.30 0.25 -915 A>G 0.33 0.40 -1175 T>G 0 0 -1590 C>T 0 0 -1665 C>T 0 0 -1690 G>A 0.51 0.41 Polymorphism site Allelic variant frequency QC (n = 328) Yak (n = 57) -259 T>A 0 1 -301_-298insCTTT 0.22 0.95 -915 A>G 0.26 0 -1175 T>G 0 1 -1590 C>T 0 1 -1665 C>T 0 1 -1690 G>A 0.53 1 Polymorphisms in promoter region numbered relative to the translation start site; Adenine of the start codon ATG is counted as +1 (GenBank accession No. AC 000163). QC = Qinchuan cattle; NY = Nanyang cattle; JX = Jiaxian cattle.
|Printer friendly Cite/link Email Feedback|
|Author:||Li, M.J.; Liu, M.; Liu, D.; Lan, X.Y.; Lei, C.Z.; Yang, D.Y.; Chen, H.|
|Publication:||Asian - Australasian Journal of Animal Sciences|
|Date:||Apr 1, 2013|
|Previous Article:||Comparative analysis of Nkx2-5/GATA4/TBX5 expression in chicken, quail and chicken-quail hybrids during the early stage of cardiac development in...|
|Next Article:||Ovarian follicular dynamics, ovarian follicular growth, oocyte yield, in vitro embryo production and repeated oocyte pick up in Thai native heifers...|