Research Article | Open Access
Tara M. Henagan, Laura K. Stewart, Laura A. Forney, Lauren M. Sparks, Neil Johannsen, Timothy S. Church, "PGC1α −1 Nucleosome Position and Splice Variant Expression and Cardiovascular Disease Risk in Overweight and Obese Individuals", PPAR Research, vol. 2014, Article ID 895734, 7 pages, 2014. https://doi.org/10.1155/2014/895734
PGC1α −1 Nucleosome Position and Splice Variant Expression and Cardiovascular Disease Risk in Overweight and Obese Individuals
PGC1, a transcriptional coactivator, interacts with PPARs and others to regulate skeletal muscle metabolism. PGC1 undergoes splicing to produce several mRNA variants, with the NTPGC1 variant having a similar biological function to the full length PGC1 (FLPGC1). CVD is associated with obesity and T2D and a lower percentage of type 1 oxidative fibers and impaired mitochondrial function in skeletal muscle, characteristics determined by PGC1 expression. PGC1 expression is epigenetically regulated in skeletal muscle to determine mitochondrial adaptations, and epigenetic modifications may regulate mRNA splicing. We report in this paper that skeletal muscle PGC1 −1 nucleosome (−1N) position is associated with splice variant NTPGC1 but not FLPGC1 expression. Division of participants based on the −1N position revealed that those individuals with a −1N phased further upstream from the transcriptional start site (UP) expressed lower levels of NTPGC1 than those with the −1N more proximal to TSS (DN). UP showed an increase in body fat percentage and serum total and LDL cholesterol. These findings suggest that the −1N may be a potential epigenetic regulator of NTPGC1 splice variant expression, and −1N position and NTPGC1 variant expression in skeletal muscle are linked to CVD risk. This trial is registered with clinicaltrials.gov, identifier NCT00458133.
Obesity and type 2 diabetes (T2D), well-known risk factors for cardiovascular disease (CVD), are marked by metabolic disturbances resulting partially from skeletal muscle mitochondrial maladaptations [1–3]. CVD itself is associated with a low percentage of type 1 oxidative skeletal muscle fibers and impaired mitochondrial function in skeletal muscle. Associated maladaptations, including decreases in mitochondrial number and function, are highly dependent on controllable risk factors which have the potential to alter the epigenome [4–6]. Peroxisome proliferator-activated receptor gamma coactivator 1 alpha (PGC1) is an important regulator of mitochondrial adaption and metabolism in several tissues, including skeletal muscle, due to its transcriptional coactivator function and binding to PPARs, estrogen related receptor alpha (ERR), nuclear receptor factor 1 (Nrf1), and others. PGC1 has been recently shown to be epigenetically regulated , and its splicing produces the novel, biologically relevant N-terminal truncated mRNA variant (NTPGC1) . In adipocytes, NTPGC1 expression acts in a similar manner to the unspliced PGC1 variant (FLPGC1) to determine mitochondrial adaptations and compensates for loss or downregulation of FLPGC1 [8, 9].
FLPGC1 and NTPGC1 expression link environmental stimuli to mitochondrial adaptations and metabolism. For example, it has been recently shown that PGC1 is hypermethylated via recruitment of DNA methyltransferase 3b (DNMT3b) in T2D muscle and upon treatment with the fatty acids palmitate and oleate . The nucleosome core position determines recruitment of DNMT3b  as well as determining the chromatin structure and access of RNA pol II and other transcription factors to the DNA template for transcription to successfully occur [11, 12]. The nucleosome position itself is partially determined by the combinatorial effects of epigenetic modifications in the cell  and may serve as an identifier of epigenetically regulated genomic loci in addition to its role in chromatin dynamics and gene regulation. Interestingly, new studies have also suggested a novel role for nucleosomes, particularly the −1 nucleosome (−1N) which is the first nucleosome within the promoter region that is directly upstream of the TSS , in determining mRNA splicing and variant expression .
Although NTPGC1 has been shown to be expressed in other tissues such as skeletal muscle, it is unknown if NTPGC1 expression is altered in association with disease risk, if its expression is linked to beneficial metabolic outcomes similar to increased FLPGC1 expression, or if NTPGC1 expression is epigenetically regulated. For this study, we used a subset of baseline muscle samples from the Health Benefits of Aerobic and Resistance Training in Type 2 Diabetes (HART-D) study  to define the relationship between the −1N position in PGC1, FLPGC1, and NTPGC1 splice variant expression in skeletal muscle and cardiometabolic risk factors in overweight/obese individuals with T2D.
2. Materials and Methods
Collection and use of skeletal muscle samples from the HART-D study [12, 14] were approved by the Institutional Review Board of the Pennington Biomedical Research Center. The clinical trial has been registered at clinicaltrials.gov, identifier NCT00458133. From a pool of 80 samples, 15 were randomly selected to represent tissue from overweight/obese (BMI > 30) individuals ( female and 11 male), aged 39–67, with T2D. Researchers involved with the analysis of the muscle tissue and resultant data were blinded to any previously obtained data from all participants. Skeletal muscle samples were collected by biopsy and cryopreserved until use in this study as previously reported .
2.2. Scanning qPCR
Genomic and mononucleosomal DNA, or DNA within one nucleosome, were isolated from 10 mg of skeletal muscle tissue ground under liquid nitrogen as previously described . Briefly, nuclei were extracted from quadriceps muscle in a 0.25 M sucrose buffer (0.25 M sucrose, 10 mM Tris-acetate pH 8.1, 1 mM EDTA, 1 mM DTT, 1 mM sodium orthovanadate, and 1 1X complete protease inhibitor tablet (Roche 11873580001)). After washing, pellets were resuspended in 0.25 sucrose buffer containing 0.1 N CaCl2 and 4 mM MgCl2 and incubated with micrococcal nuclease (MNase; Roche) for mononucleosomal DNA extraction or without MNase for genomic DNA extraction for 15 min at 37°C. EDTA was added to a final concentration of 5 mM. Pellets were formed by centrifugation and lysed with 0.25 M lysis buffer (50 mM Tris, pH 8.1, 10 mM EDTA, 1% SDS, and 1 1X complete protease inhibitor tablet (Roche 11873580001)). All samples were treated with 0.1 mg/mL of proteinase K (Qiagen) overnight at 37°C to remove histone proteins.
Scanning qPCR was performed as previously described [16, 17]. Overlapping primers (sequences presented in Table 1) were designed to cover the PGC1 gene promoter region, ranging from ~−800 nucleotide (nt) to the −100 nt (Figure 1). PCR products for both mononucleosomal and genomic DNA samples were run on a 1.5% agarose gel and visualized on a Molecular Imager Gel Doc XR (Biorad, Hercules, CA). Densitometry was performed using MacBiophotonics ImageJ (Bethesda, MD), with mononucleosomal band intensity being divided by the intensity of the corresponding input genomic DNA.
(a) Scanning PCR primers
Quadriceps muscle was ground under liquid nitrogen with a mortar and pestle, homogenized with Trizol reagent per the manufacturer’s protocol (Life Technologies, Foster City, CA) and column purified with a RNeasy kit (Qiagen, Valencia, CA). cDNA synthesis was carried out using M-MLV reverse transcriptase per the manufacturer’s protocol (Promega, Madison, WI). qRT-PCR was performed on the ABI7900HT platform with SYBR Green PCR Master Mix (Applied Biosystems, Foster City, CA). Previously published primers targeted to FLPGCC1 and NTPGC1 were used . Cyclophilin B was used as the internal control. Data were analyzed using the standard curve method.
2.4. Determining −1N Subject Categories
Subjects were divided into two groups based on the results of our scanning qPCR, which was used to map the −1N position in the PGC1 promoter region. Subjects with a −1N positioned approximately between the −290 and −440 nt were designated to the upstream category (UP). Subjects with a −1 nucleosome positioned approximately between the −170 and −320 nt were designated to the downstream (DN) category.
2.5. Statistical Analyses
Previously measured anthropometric and cardiometabolic risk factors  and FLPGC1 and NTPGC1 mRNA expressions were averaged for each group, UP and DN. All data were analyzed by Student’s -test with being considered significant using GraphPad Prism 4.0 software.
3.1. Nucleosome Position and Splice Variant Expression
Several environmentally induced PGC1 splice variants are expressed in skeletal muscle, including FLPGC1 and NTPGC1 which have both been shown to act similarly to induce beneficial mitochondrial adaptations and improve metabolism in adipocytes [8, 18]. The −1N, which is important in determining gene transcription and expression and has recently been shown to determine splice variant expression , was mapped in the PGC1 promoter to determine its association with splice variant expression. Overlapping primer pairs were designed to span the PGC1 promoter region from approximately the −800 nt to the −100 nt (Figure 1(a)), and scanning qPCR was used to map the −1N position in PGC1. Scanning qPCR gives a high resolution map of nucleosome position and occupancy at a specific genomic locus and is a common method for targeted nucleosome mapping [17, 19]. The −1 nucleosome showed phasing, ranging from the −440 nucleotide (nt) to the −170 nt (Figure 1(b)). Analysis of the phasing showed that individuals could be grouped based on the –1N being shifted either further upstream, away from the TSS (UP, ), or closer downstream, toward the TSS (DN, ) (Figure 1(b)).
After grouping participant data based on the −1N position (Figure 1(b)), mRNA expressions of FLPGC1 and NTPGC1 were measured and analyzed. There was no significant difference in FLPGC1 mRNA expression between groups (; Figure 2(a)). However, we observed a significant decrease in NTPGC1 mRNA expression in UP compared to DN (; Figure 2(b)). These data suggest that nucleosome positioning in PGC1 may play a role in splice variant expression.
3.2. Nucleosome Position and CVD Risk Factors
CVD risk is associated with obesity and T2D and individuals with lower expression of skeletal muscle PGC1 exhibit higher disease risk . When subjects were divided into groups based on −1N position in PGC1, no differences in body weight (Figure 3(a)) or age (UP y; DN y) existed between UP and DN. Interestingly, percent body fat was lower in UP compared with DN (), although percentage of lean mass was not different (Figure 3(a)). BMI was not statistically significant between groups (Figure 3(a)). No significant differences in systolic (SBP) or diastolic blood pressure (DBP) were evident (Figure 3(b)). In UP, total serum cholesterol () and low density lipoprotein (LDL; ) cholesterol were lower, and there was no difference in high density lipoprotein (HDL) cholesterol between groups (Figure 3(c)). Serum triglycerides, free fatty acids (FFA), fasting blood glucose, and insulin levels were not different between groups (Figure 3(d)). These data show that individuals with a −1N positioned proximal to the TSS in the PGC1 promoter and with higher levels of NTPGC1 exhibit increased CVD risk as assessed by adiposity, total cholesterol, and LDL cholesterol.
(a) Body weight and composition
(b) Blood pressure
(d) Triglycerides, free fatty acid, and insulin sensitivity
We mapped the −1N within the PGC1 promoter region to provide insight into differential epigenetic regulation of PGC1, including FLPGC1 and NTPGC1 splice variant expression, in overweight/obese individuals with T2D and to provide evidence that alterations in PGC1 −1N position and splice variant expression are associated with differential CVD risk . Our results provide evidence that proximal positioning of the −1N in PGC1 is associated with increased CVD risk and increased NTPGC1 expression. Interestingly, we found that, in overweight/obese individuals with T2D exhibiting higher adverse CVD risk, the −1N was positioned over a regulatory epigenetic site in the PGC1 promoter , which may be dependent on nucleosome positioning . The −1N was associated with the degree of adiposity but not fasting insulin or glucose levels, with those individuals exhibiting a −1N proximal to the TSS, over the regulatory epigenetic site, being more obese and having higher levels of total and LDL cholesterol. These data suggest that the chromatin structure of PGC1 [7, 14] is related to the degree of overweight/obesity and obesity-associated CVD risk in individuals with T2D. Indeed, others have noted that epigenetic regulation of skeletal muscle PGC1 resulting in decreased gene expression leads to a reduction in skeletal muscle mitochondrial number and decreased expression of PGC1 target genes in association with disease state, specifically insulin resistance and T2D [7, 21].
Although −1N position was not associated with alterations in FLPGC1 expression, those individuals with the −1N more proximal to the TSS showed an increase in NTPGC1 expression in addition to increased CVD risk. Recent research indicates a role of the −1N in regulating transcript processing via mRNA splicing . Here, we found differential −1N in PGC1 predicted NTPGC1 but not FLPGC1 expression in skeletal muscle, with a significant increase in NTPGC1 in the skeletal muscle of DN individuals who had their −1N positioned proximal to the TSS. These observed differences in splice variant expression between UP and DN groups in the present study suggest that the −1N may regulate PGC1 splicing and variant expression. Although the mechanism of this regulation is yet to be explored, it is possible that phasing of the −1N may determine periodicity of downstream nucleosomes or decrease transcription elongation rate, leading to intron 6 inclusion and NTPGC1 expression . Importantly, NTPGC1 has been shown to translocate to the nucleus in adipose tissue, where it acts in a similar manner to FLPGC1 to regulate nuclear-encoded mitochondrial gene expression [8, 22, 23]. Splice variant expression may also be important in determining mitochondrial function and number in skeletal muscle [9, 24–26], and NTPGC1 expression may be increased to compensate for lack of change in FLPGC1 . Indeed, it has recently been reported that upregulation of NTPGC1 in myotubes increases glucose transporter and mitochondrial gene expression, which may account partially for the similarities between FLPGC1 and NTPGC1 in their insulin sensitizing effects . Although it is possible that NTPGC1 has differential functions in various tissues, we speculate that in the present study NTPGC1 is upregulated to compensate for lack of change in FLPGC1 in the obese and diabetic state. However, the magnitude of NTPGC1 increase may not be adequate or sufficient to ameliorate obesity-associated metabolic dysfunction and CVD risk in this population. The lack of data on differential expression of PGC1 target genes in the present study is a limitation and the focus of future studies on these particular differences in obese and diabetic individuals will provide insight into the molecular mechanisms linking PGC1 DNA methylation to −1N positioning, the role of epigenetics in determining splice variant expression, and the associations of splice variants’ expression to PGC1 target gene activation and disease state. Additionally, further analysis of splice variant expression levels should be conducted in lean versus overweight/obese individuals as well as other diseased and nondiseased populations to determine a minimum level of FL- or NTPGC1 that is sufficient to decrease CVD risk.
Our data revealed that downstream −1N in the PGC1 promoter is associated with higher adiposity and adverse health risk, whereas those individuals with an upstream −1N had lower adiposity and obesity-related risk for CVD. Our data suggest that −1N positioning may be a potential epigenetic mechanism that regulates NTPGC1 splice variant expression, and this variant expression is linked to CVD risk in overweight/obese individuals with T2D.
|BMI:||Body mass index|
|DNMT3b:||DNA methyltransferase 3b|
|DBP:||Diastolic blood pressure|
|ERR:||Estrogen related receptor alpha|
|FFA:||Free fatty acids|
|HART-D:||Health Benefits of Aerobic and Resistance Training in Type 2 Diabetes|
|HDL:||High density lipoprotein|
|FLPGC1:||Full length peroxisome proliferator-activated receptor gamma coactivator 1 alpha|
|LDL:||Low density lipoprotein|
|Nrf1:||Nuclear receptor factor 1|
|NTPGC1:||N truncated peroxisome proliferator-activated receptor gamma coactivator 1 alpha|
|PGC1:||Peroxisome proliferator-activated receptor gamma coactivator 1 alpha|
|PPAR:||Peroxisome proliferator-activated receptor|
|SBP:||Systolic blood pressure|
|T2D:||Type 2 diabetes|
|TSS:||Transcriptional start site|
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This work was supported by NIH NIDDK 5T32-DK064584-09 (Tara M. Henagan), NIH COBRE (8P20-GM103528-07), and NIH NIDDK DK068298 (Timothy S. Church).
- R. A. DeFronzo and D. Tripathy, “Skeletal muscle insulin resistance is the primary defect in type 2 diabetes,” Diabetes Care, vol. 32, pp. S157–S163, 2009.
- K. E. Shortreed, M. P. Krause, J. H. Huang et al., “Muscle-specific adaptations, impaired oxidative capacity and maintenance of contractile function characterize diet-induced obese mouse skeletal muscle,” PLoS ONE, vol. 4, no. 10, Article ID e7293, 2009.
- L. M. Sparks, H. Xie, R. A. Koza et al., “A high-fat diet coordinately downregulates genes required for mitochondrial oxidative phosphorylation in skeletal muscle,” Diabetes, vol. 54, no. 7, pp. 1926–1933, 2005.
- Z. Vucetic, J. Kimmel, K. Totoki, E. Hollenbeck, and T. M. Reyes, “Maternal high-fat diet alters methylation and gene expression of dopamine and opioid-related genes,” Endocrinology, vol. 151, no. 10, pp. 4756–4764, 2010.
- C. Gallou-Kabani and C. Junien, “Nutritional epigenomics of metabolic syndrome: new perspective against the epidemic,” Diabetes, vol. 54, no. 7, pp. 1899–1906, 2005.
- R. Feil and M. F. Fraga, “Epigenetics and the environment: emerging patterns and implications,” Nature Reviews Genetics, vol. 13, no. 2, pp. 97–109, 2012.
- R. Barrès, M. E. Osler, J. Yan et al., “Non-CpG methylation of the PGC-1α promoter through DNMT3B controls mitochondrial density,” Cell Metabolism, vol. 10, no. 3, pp. 189–198, 2009.
- Y. Zhang, P. Huypens, A. W. Adamson et al., “Alternative mRNA splicing produces a novel biologically active short isoform of PGC-1α,” The Journal of Biological Chemistry, vol. 284, no. 47, pp. 32813–32826, 2009.
- J. L. Ruas, J. P. White, R. R. Rao et al., “A PGC-1α isoform induced by resistance training regulates skeletal muscle hypertrophy,” Cell, vol. 151, no. 6, pp. 1319–1331, 2012.
- H. Takeshima, I. Suetake, H. Shimahara, K. Ura, S. Tate, and S. Tajima, “Distinct DNA methylation activity of Dnmt3a and Dnmt3b towards naked and nucleosomal DNA,” The Journal of Biochemistry, vol. 139, no. 3, pp. 503–515, 2006.
- R. F. Luco, M. Allo, I. E. Schor, A. R. Kornblihtt, and T. Misteli, “Epigenetics in alternative pre-mRNA splicing,” Cell, vol. 144, no. 1, pp. 16–26, 2011.
- L. M. Sparks, N. M. Johannsen, T. S. Church et al., “Nine months of combined training improves Ex vivo skeletal muscle metabolism in individuals with type 2 diabetes,” The Journal of Clinical Endocrinology and Metabolism, vol. 98, no. 4, pp. 1694–1702, 2013.
- G. E. Zentner and S. Henikoff, “Regulation of nucleosome dynamics by histone modifications,” Nature Structural & Molecular Biology, vol. 20, no. 3, pp. 259–266, 2013.
- T. S. Church, S. N. Blair, S. Cocreham et al., “Effects of aerobic and resistance training on hemoglobin A1c levels in patients with type 2 diabetes: a randomized controlled trial,” The Journal of the American Medical Association, vol. 304, no. 20, pp. 2253–2262, 2010.
- Z. Li, J. Schug, G. Tuteja, P. White, and K. H. Kaestner, “The nucleosome map of the mammalian liver,” Nature Structural and Molecular Biology, vol. 18, no. 6, pp. 742–746, 2011.
- K. M. Creamer, G. Job, S. Shanker et al., “The Mi-2 homolog Mit1 actively positions nucleosomes within heterochromatin to suppress transcription,” Molecular and Cellular Biology, vol. 34, no. 11, pp. 2046–2061, 2014.
- E. A. Sekinger, Z. Moqtaderi, and K. Struhl, “Intrinsic histone-DNA interactions and low nucleosome density are important for preferential accessibility of promoter regions in yeast,” Molecular Cell, vol. 18, no. 6, pp. 735–748, 2005.
- H.-J. Jun, Y. Joshi, Y. Patil, R. C. Noland, and J. S. Chang, “NT-PGC-1α activation attenuates high-fat diet-induced obesity by enhancing brown fat thermogenesis and adipose tissue oxidative metabolism,” Diabetes, vol. 63, no. 11, pp. 3615–3625, 2014.
- J. Infante, G. L. Law, and E. Young, “Analysis of nucleosome positioning using a nucleosome-scanning assay,” in Chromatin Remodeling, R. H. Morse, Ed., pp. 63–87, Humana Press, 2012.
- C. Handschin and B. M. Spiegelman, “The role of exercise and PGC1α in inflammation and chronic disease,” Nature, vol. 454, no. 7203, pp. 463–469, 2008.
- R. Barrès, J. Yan, B. Egan et al., “Acute exercise remodels promoter methylation in human skeletal muscle,” Cell Metabolism, vol. 15, no. 3, pp. 405–411, 2012.
- J. S. Chang, P. Huypens, Y. Zhang, C. Black, A. Kralli, and T. W. Gettys, “Regulation of NT-PGC-1α subcellular localization and function by protein kinase A-dependent modulation of nuclear export by CRM1,” The Journal of Biological Chemistry, vol. 285, no. 23, pp. 18039–18050, 2010.
- J. S. Chang and T. W. Gettys, “Analyzing phosphorylation-dependent regulation of subcellular localization and transcriptional activity of transcriptional coactivator NT-PGC-1α,” Methods in Molecular Biology, vol. 952, pp. 163–173, 2013.
- D. Knutti and A. Kralli, “PGC-1, a versatile coactivator,” Trends in Endocrinology and Metabolism, vol. 12, no. 8, pp. 360–365, 2001.
- H. Liang and W. F. Ward, “PGC-1α: a key regulator of energy metabolism,” American Journal of Physiology: Advances in Physiology Education, vol. 30, no. 4, pp. 145–151, 2006.
- P. J. Fernandez-Marcos and J. Auwerx, “Regulation of PGC-1α, a nodal regulator of mitochondrial biogenesis,” American Journal of Clinical Nutrition, vol. 93, no. 4, pp. 884S–890S, 2011.
- X. Wen, J. Wu, J. S. Chang et al., “Effect of exercise intensity on isoform-specific expressions of NT-PGC-1 mRNA in mouse skeletal muscle,” BioMed Research International, vol. 2014, Article ID 402175, 11 pages, 2014.
Copyright © 2014 Tara M. Henagan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.