- Short Report
- Open Access
Two new splice variants in porcine PPARGC1A
BMC Research Notesvolume 1, Article number: 138 (2008)
Peroxisome proliferator-activated receptor γ coactivator 1α (PPARGC1A) is a coactivator with a vital and central role in fat and energy metabolism. It is considered to be a candidate gene for meat quality in pigs and is involved in the development of obesity and diabetes in humans. How its many functions are regulated, is however still largely unclear. Therefore a transcription profile of PPARGC1A in 32 tissues and 4 embryonic developmental stages in the pig was constructed by screening its cDNA for possible splice variants with exon-spanning primers.
This led to the discovery of 2 new splice variants in the pig, which were subsequently also detected in human tissues. In these variants, exon 8 was either completely or partly (the last 66 bp were conserved) spliced out, potentially coding for a much shorter protein of respectively 337 and 359 amino acids (aa), of which the first 291 aa would be the same compared to the complete protein (796 aa).
Considering the functional domains of the PPARGC1A protein, it is very likely these splice variants considerably affect the function of the protein and alternative splicing could be one of the mechanisms by which the diverse functions of PPARGC1A are regulated.
Peroxisome proliferator-activated receptor γ coactivator 1α (PPARGC1A) is a transcriptional coactivator with many diverse functions and has a pivotal role in fat and energy metabolism. This cold- and exercise-inducible gene is crucial to adaptive thermogenesis and is an essential regulator of adipogenesis, adipocyte differentiation and mitochondrial biogenesis/respiration [1–4]. Recently, it has been shown that it is also involved in angiogenesis . It exerts its function through a whole range of nuclear hormone receptors and other transcription factors, and is primarily expressed in tissues with high energy demands . Besides having an important influence on the regulation and composition of the body weight, it also is an important factor in determining muscle fibre type composition [7–9]. It has been shown that PPARGC1A increases the amount of oxidative muscle fibres, and that it also is expressed at a higher level in these muscle fibres.
For several reasons, porcine PPARGC1A is an interesting candidate gene for meat quality, an economically important and complex characteristic which is composed of many different traits. Associations have been found between mutations in the coding region of PPARGC1A and certain fat characteristics in the pig [10–12]. Other interesting findings are that PPARGC1A is the only candidate gene so far that was located in the QTL region for leaf fat weight and backfat on chromosome 8p21 [11, 13], and that other candidate genes for meat quality, like GLUT4, are regulated by PPARGC1A .
As explained above, PPARGC1A has many functions which can strongly differ between tissues. It has been shown that there is a variation in mRNA expression in the pig, not only between tissues, but also between different locations within the longissimus dorsi muscle . Only very little is known about this multifunctional gene in the pig and its possible use as a selection marker in the pig industry. In order to get a better understanding of the regulation of the many functions of PPARGC1A, exon-spanning primers were used to construct a detailed transcription profile of its presence in 32 different tissues and several embryonic developmental stages in the pig. The aim was to identify possible splice variants, because they could provide an explanation for the regulation of the tissue-dependent functions of PPARGC1A.
Tissue samples were collected from a freshly slaughtered female, commercial, hybrid pig and immediately submerged in RNA later (Sigma-Aldrich, Bornem, Belgium), according to the instructions manual. Testis was collected from a similar male pig. Total RNA was extracted with the Aurum Total RNA Fatty and Fibrous Tissue Kit (Bio-Rad, Nazareth, Belgium), according to the manufacturer's protocol which included an on-column DNase treatment. Ovaries were collected at a local slaughterhouse from pigs at slaughter age, and used for in vitro embryo production as described in Bijttebier et al. . RNA extraction from embryonic samples (for the 2–4 cell, 8 cell, morula and blastocyst stage respectively 15, 12, 8 and 6 pooled embryos were used) was performed with the PicoPure RNA Isolation Kit (Arcturus, Mountain View, USA), according to the instructions manual, after which a DNase treatment was carried out with RQ1 RNase-free DNase (Promega, Leiden, The Netherlands). Both DNase treatments were verified by a minus reverse transcription (RT) control PCR and RNA integrity was checked, as described in Erkens et al. . Also, RNA purity and concentration were measured with the ND-1000 Spectrophotometer (NanoDrop, Wilmington, USA). Next, the iScript cDNA Synthesis Kit (which contains both oligo dT and random primers; Bio-Rad, Nazareth, Belgium) was used to convert approximately 1 μg of RNA from each sample to cDNA, according to the manufacturer's protocol. This RT step was verified by a control PCR , in which a no-template control was included to check for DNA contamination. Ready-to-use human cDNA from kidney and liver was provided by Prof. Vandesompele (Department of Pediatrics and Medical Genetics, Ghent University), to verify whether the detected splice variants also occur in human tissues.
Porcine sequences [GenBank:AH013726] and [GenBank:AY346131], found in the NCBI database , were used for the design of exon-spanning primers with Primer3  (Table 1). This way, possible splice variants for each of the 13 exons (except the outer ones) could be detected. During primer design, Mfold  and Blast  were used to check for possible secondary structures and primer specificity, respectively. PCR conditions for each primer were optimized with FastStart Taq DNA Polymerase (Roche, Vilvoorde, Belgium) and included a no-template control. Also, a genomic DNA control was included to check for possible amplification of pseudogenes. The annealing temperature used for all primer pairs was 60°C.
The amplicons from all primer pairs, except PGC1A+/-Ex7,9, were sequenced by direct sequencing. Because the use of PGC1A+/-Ex7,9 resulted in multiple amplicons, the GENECLEAN II Kit (Qbiogene, Brussels, Belgium) was used to first isolate and purify the separate amplicons from the agarose gel, before sequencing. Sequencing of the amplicons was conducted on an Applied Biosystems 3730xl DNA Analyser with the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Lennik, Belgium), according to the manufacturer's protocol.
Results and discussion
All exon-spanning primer pairs in all tissues resulted in one amplicon of the expected length (Table 1), except for PGC1A+/-Ex7,9. With this last primer pair, 3 different amplicons were detected and the results were tissue-dependent (Table 2). The forward and reverse primer of PGC1A+/-Ex7,9 are located in exon 7 and exon 9, respectively, which means its amplicon completely contains exon 8. Besides the expected fragment of 1017 bp (which incorporates all 916 bp of exon 8), 2 new splice variants of PPARGC1A in the pig were detected. Sequencing revealed that in these splice variants exon 8 was either partly (the last 66 bp of exon 8 were conserved) or completely spliced out, which resulted in amplicons of 167 and 101 bp, respectively (Figure 1, 2). Additional PCRs with longer amplicons (exon 4–9 and exon 4–12) confirmed the existence of both exon 8 splice variants and indicate that the other exons were preserved in the rest of the transcript. Taken together with the fact that no pseudogenes for PPARGC1A are described in any species and the use of a genomic DNA control, this makes it unlikely that our observations are the result of pseudogenes.
As can be seen from Figure 2, both splice variant boundaries have GC-AG splice sites, instead of following the more common GT-AG rule. Table 2 shows that one or both splice variants were found in almost every tissue that was tested. They could also be detected in the pre-implantation embryonic stages, which implies PPARGC1A is involved in early development. An interesting finding however, is the fact that the complete amplicon (1017 bp) was not detected in any of the 4 embryonic stages. This could indicate certain functions of PPARGC1A are switched off or altered during early development, but without a functional analysis of the splice variants it is not possible to discuss the effects on its functionality or to give an explanation for these findings. In testis, a PCR artefact (221 bp) was detected with a 92% identity to human outer dense fibre of sperm tails 2 (ODF2) and was deposited in GenBank as an EST [GenBank:EY122774].
Only very little is known about the existence of splice variants of PPARGC1A in any species, and up until now nothing was known about it in the pig. There have been previous reports suggesting the existence of splice variants in which exon 8 was possibly spliced out, in rat skeletal muscle and brown adipose tissue, but their sequence was not determined [21, 22]. However, this is the first study giving a detailed description of the actual presence and sequence of 2 splice variants in a whole range of porcine tissues.
The 2 newly identified splice variants possibly give rise to a much shorter protein of respectively 359 and 337 aa, depending on whether exon 8 is partly or completely spliced out (Figure 3). This is much shorter than the complete protein with 796 aa. The first 291 aa of both variants are identical to the complete protein, but the rest of the aa sequence is completely different. It can be expected that this will have an important influence on the function of the produced protein (Figure 3). In the study by Baar et al.  an increase of a smaller PPARGC1A protein was detected after exercise in rats. This was consistent with the increase of the smaller cDNA band they detected, although it was not established if that protein originated from the shorter mRNA and if this had any functional significance.
The PPARGC1A protein can generally be divided into 3 regions (Figure 3a). The N-terminal region consists of a transcriptional activation domain which contains an essential LXXLL motif, and is involved in the activation of many transcription factors. A less distinct, central region contains both an inhibitory domain and several interaction domains (PPARγ, NRF-1, MEF2C). At the C-terminal end, an RNA processing domain is located, which comprises 2 serine-arginine-rich domains (SR) and an RNA recognition motif (RRM) [1, 23–25]. Figure 3 shows that in the putative protein from both splice variants the N-terminal activation domain is conserved. The central region however is only partly conserved. The aa in the first part of the inhibitory and NRF-1 domain are conserved, but the aa of the complete PPARγ interaction region are altered. The second part of the inhibitory domain and NRF-1 interaction region is either altered or absent. The RNA processing domain at the C-terminal end of the complete protein is completely absent in the splice variants. These results suggest that the putative proteins from both splice variants show some remarkable alterations and this is likely to have a large impact on the function of PPARGC1A.
Currently, human medicine shows a great interest in PPARGC1A, because of the important role this gene plays in the worldwide problems concerning obesity, insulin resistance and correlated diseases, such as type 2 diabetes mellitus [26, 27]. It has also been shown recently that a lower expression of PPARGC1A is involved in the onset of multiple neurodegenerative diseases, like Parkinson's, Alzheimer's and Huntington's disease . Because of its significance, human kidney and liver tissue were also tested for the presence of the newly detected splice variants in the pig. This showed that both splice variants were also found in human liver and only the shortest splice variant in human kidney. In porcine kidney, both splice variants were detected, indicating the existence of possible species differences. The discovery of these new splice variants could therefore be of importance for the human research regarding PPARGC1A as well.
The results from this study contribute to a better understanding of this complex gene and are of possible use not only for research in the pig industry regarding meat quality and carcass composition, but also for human research. Considering the functional domains of the PPARGC1A protein, it is very likely these splice variants considerably affect the function of the protein and alternative splicing could be one of the mechanisms by which the diverse functions of PPARGC1A are regulated.
Puigserver P, Wu Z, Park CW, Graves R, Wright M, Spiegelman BM: A cold-inducible coactivator of nuclear receptors linked to adaptive thermogenesis. Cell. 1998, 92: 829-839. 10.1016/S0092-8674(00)81410-5.
Spiegelman BM, Puigserver P, Wu Z: Regulation of adipogenesis and energy balance by PPARγ and PGC-1. Int J Obes Relat Metab Disord. 2000, 24 Suppl 4: S8-S10.
Yoon JC, Puigserver P, Chen GX, Donovan J, Wu ZD, Rhee J, Adelmant G, Stafford J, Kahn CR, Granner DK, Newgard CB, Spiegelman BM: Control of hepatic gluconeogenesis through the transcriptional coactivator PGC-1. Nature. 2001, 413: 131-138. 10.1038/35093050.
Medina-Gomez G, Gray S, Vidal-Puig A: Adipogenesis and lipotoxicity: role of peroxisome proliferator-activated receptor γ (PPARγ) and PPARγ coactivator-1 (PGC1). Public Health Nutr. 2007, 10: 1132-1137. 10.1017/S1368980007000614.
Arany Z, Foo SY, Ma Y, Ruas JL, Bommi-Reddy A, Girnun G, Cooper M, Laznik D, Chinsomboon J, Rangwala SM, Baek KH, Rosenzweig A, Spiegelman BM: HIF-independent regulation of VEGF and angiogenesis by the transcriptional coactivator PGC-1α. Nature. 2008, 451: 1008-1012. 10.1038/nature06613.
Larrouy D, Vidal H, Andreelli F, Laville M, Langin D: Cloning and mRNA tissue distribution of human PPARγ coactivator-1. Int J Obesity. 1999, 23: 1327-1332. 10.1038/sj.ijo.0801106.
Dulloo AG, Samec S: Uncoupling proteins: their roles in adaptive thermogenesis and substrate metabolism reconsidered. Br J Nutr. 2001, 86: 123-139. 10.1079/BJN2001412.
Lin J, Wu H, Tarr PT, Zhang C-Y, Wu Z, Boss O, Michael LF, Puigserver P, Isotani E, Olson EN, Lowell BB, Bassel-Duby R, Spiegelman BM: Transcriptional co-activator PGC-1α drives the formation of slow-twitch muscle fibres. Nature. 2002, 418: 797-801. 10.1038/nature00904.
Mortensen OH, Frandsen L, Schjerling P, Nishimura E, Grunnet N: PGC-1α and PGC-1β have both similar and distinct effects on myofiber switching toward an oxidative phenotype. Am J Physiol-Endocrinol Metab. 2006, 291: E807-816. 10.1152/ajpendo.00591.2005.
Kunej T, Wu XL, Berlic TM, Michal JJ, Jiang Z, Dovc P: Frequency distribution of a Cys430Ser polymorphism in peroxisome proliferator-activated receptor-γ coactivator-1 (PPARGC1) gene sequence in Chinese and Western pig breeds. J Anim Breed Genet. 2005, 122: 7-11. 10.1111/j.1439-0388.2004.00498.x.
Jacobs K, Rohrer G, Van Poucke M, Piumi F, Yerle M, Barthenschlager H, Mattheeuws M, Van Zeveren A, Peelman LJ: Porcine PPARGC1A (peroxisome proliferative activated receptor gamma coactivator 1A): coding sequence, genomic organization, polymorphisms and mapping. Cytogenet Genome Res. 2006, 112: 106-113. 10.1159/000087521.
Stachowiak M, Szydlowski M, Cieslak J, Switonski M: SNPs in the porcine PPARGC1a gene: interbreed differences and their phenotypic effects. Cell Mol Biol Lett. 2007, 12: 231-239. 10.2478/s11658-006-0066-7.
Rohrer GA, Keele JW: Identification of quantitative trait loci affecting carcass composition in swine: I. Fat deposition traits. J Anim Sci. 1998, 76: 2247-2254.
Michael LF, Wu Z, Cheatham RB, Puigserver P, Adelmant G, Lehman JJ, Kelly DP, Spiegelman BM: Restoration of insulin-sensitive glucose transporter (GLUT4) gene expression in muscle cells by the transcriptional coactivator PGC-1. Proc Natl Acad Sci USA. 2001, 98: 3820-3825. 10.1073/pnas.061035098.
Erkens T, van Poucke M, Vandesompele J, Goossens K, Van Zeveren A, Peelman LJ: Development of a new set of reference genes for normalization of real-time RT-PCR data of porcine backfat and longissimus dorsi muscle, and evaluation with PPARGC1A. BMC Biotechnol. 2006, 6: 41-10.1186/1472-6750-6-41.
Bijttebier J, Van Soom A, Meyer E, Mateusen B, Maes D: Preovulatory follicular fluid during in vitro maturation decreases polyspermic fertilization of cumulus-intact porcine oocytes in vitro maturation of porcine oocytes. Theriogenology. 2008, 70: 715-724. 10.1016/j.theriogenology.2008.04.046.
National Center for Biotechnology Information. [http://www.ncbi.nlm.nih.gov]
Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. 2000, Totowa: Humana Press, 365-386. [http://frodo.wi.mit.edu/cgi-bin/primer3/primer3_www.cgi]
Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31: 3406-3415. 10.1093/nar/gkg595. [http://frontend.bioinfo.rpi.edu/applications/mfold/cgi-bin/dna-form1.cgi]
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410. [http://www.ncbi.nlm.nih.gov/blast]
Kakuma T, Wang Z-W, Pan W, Unger RH, Zhou Y-T: Role of leptin in peroxisome proliferator-activated receptor gamma coactivator-1 expression. Endocrinology. 2000, 141: 4576-4582. 10.1210/en.141.12.4576.
Baar K, Wende AR, Jones TE, Marison M, Nolte LA, Chen M, Kelly DP, Holloszy JO: Adaptations of skeletal muscle to exercise: rapid increase in the transcriptional coactivator PGC-1. FASEB J. 2002, 16: 1879-1886. 10.1096/fj.02-0367com.
Wu Z, Puigserver P, Andersson U, Zhang C, Adelmant G, Mootha V, Troy A, Cinti S, Lowell B, Scarpulla RC, Spiegelman BM: Mechanisms controlling mitochondrial biogenesis and respiration through the thermogenic coactivator PGC-1. Cell. 1999, 98: 115-124. 10.1016/S0092-8674(00)80611-X.
Puigserver P, Spiegelman BM: Peroxisome proliferator-activated receptor-γ coactivator 1α (PGC-1α): transcriptional coactivator and metabolic regulator. Endocr Rev. 2003, 24: 78-90. 10.1210/er.2002-0012.
Sano M, Tokudome S, Shimizu N, Yoshikawa N, Ogawa C, Shirakawa K, Endo J, Katayama T, Yuasa S, Ieda M, Makino S, Hattori F, Tanaka H, Fukuda K: Intramolecular control of protein stability, subnuclear compartmentalization, and coactivator function of peroxisome proliferator-activated receptor γ coactivator 1α. J Biol Chem. 2007, 282: 25970-25980. 10.1074/jbc.M703634200.
Ek J, Andersen G, Urhammer SA, Gæde PH, Drivshom T, Borch-Johnsen K, Hansen T, Pedersen O: Mutation analysis of peroxisome proliferator activated receptor-γ coactivator-1 (PGC-1) and relationships of identified amino acid polymorphisms to type II diabetes mellitus. Diabetologia. 2001, 44: 2220-2226. 10.1007/s001250100032.
Hammarstedt A, Jansson P-A, Wesslau C, Yang X, Smith U: Reduced expression of PGC-1 and insulin-signaling molecules in adipose tissue is associated with insulin resistance. Biochem Biophys Res Commun. 2003, 301: 578-582. 10.1016/S0006-291X(03)00014-7.
St-Pierre J, Drori S, Uldry M, Silvaggi JM, Rhee J, Jäger S, Handschin C, Zheng K, Lin J, Yang W, Simon DK, Bachoo R, Spiegelman BM: Suppression of reactive oxygen species and neurodegeneration by the PGC-1 transcriptional activators. Cell. 2006, 127: 397-408. 10.1016/j.cell.2006.09.024.
The authors would like to thank Jo Bijttebier and Isabel Lemahieu (Department of Obstetrics, Reproduction and Herd Health, Ghent University), Filip Barbé (Laboratory of Virology, Faculty of Veterinary Medicine, Ghent University), and Prof. Vandesompele and Nurten Yigit (Department of Pediatrics and Medical Genetics, Ghent University) for providing samples. This work was supported by Ghent University BOF-grant 01J02707.
The authors declare that they have no competing interests.
TE participated in the study design, performed part of the experimental procedures and was the primary author of the manuscript. KB performed part of the experimental procedures and helped to draft the manuscript. AVZ and LJP participated in the design of the project, helped to draft the manuscript and supervised the study. All authors read and approved the final manuscript.