O-GlcNAcylation mediates the control of cytosolic phosphoenolpyruvate carboxykinase activity via Pgc1α

Pedro Latorre; Luis Varona; Carmen Burgos; José A. Carrodeguas; Pascual López-Buesa

doi:10.1371/journal.pone.0179988

Abstract

PGC1α is a coactivator of many transcription factors and cytosolic phosphoenolpyruvate carboxykinase (PCK1) is a key enzyme for gluconeogenesis. PGC1α interacts with the transcription factor PPARγ to stimulate PCK1 expression and thus de novo glucose synthesis. These proteins are not only important for central energy metabolism but also for supplying intermediates for other metabolic pathways, including lipidogenesis and protein synthesis and might therefore be important factors in the ethiopathogenesis of metabolic disorders like diabetes but also in other pathologies like cancer. Since polymorphisms in these proteins have been related to some phenotypic traits in animals like pigs and PGC1α G482S polymorphism increases fat deposition in humans, we have investigated the molecular basis of such effects focusing on a commonly studied polymorphism in pig Pgc1α, which changes a cysteine at position 430 (WT) of the protein to a serine (C430S). Biochemical analyses show that Pgc1α WT stimulates higher expression of human PCK1 in HEK293T and HepG2 cells. Paradoxically, Pgc1α WT is less stable than Pgc1α p.C430S in HEK293T cells. However, the study of different post-translational modifications shows a higher O-GlcNAcylation level of Pgc1α p.C430S. This higher O-GlcNAcylation level significantly decreases the interaction between Pgc1α and PPARγ demonstrating the importance of post-translational glycosylation of PGC1α in the regulation of PCK1 activity. This, furthermore, could explain at least in part the observed epistatic effects between PGC1α and PCK1 in pigs.

Citation: Latorre P, Varona L, Burgos C, Carrodeguas JA, López-Buesa P (2017) O-GlcNAcylation mediates the control of cytosolic phosphoenolpyruvate carboxykinase activity via Pgc1α. PLoS ONE 12(6): e0179988. https://doi.org/10.1371/journal.pone.0179988

Editor: Ramona Natacha PENA i SUBIRÀ, Universitat de Lleida, SPAIN

Received: March 16, 2017; Accepted: June 7, 2017; Published: June 23, 2017

Copyright: © 2017 Latorre et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All relevant data are within the paper and its Supporting Information file.

Funding: This work has been funded by grants AGL2015-66177 to P.L.B. and J.A.C. from Ministerio de Economia, Industria y Competitividad (http://www.mineco.gob.es/portal/site/mineco/), UZ 2014-CIE-03 to P.L.B. and UZ-2015-BIO-01 to J.A.C. from University of Zaragoza (http://www.unizar.es). Funders' role has been limited to financial support. P.L. is supported by a predoctoral fellowship from ‘la Caixa’ Foundation (https://www.fundacionbancarialacaixa.org/corporate/home_es.html).

Competing interests: The authors have declared that no competing interests exist.

Introduction

PGC1α was first discovered as a cold-inducible protein involved in muscle thermogenesis [1]. To coordinate metabolic responses PGC1α interacts with a large number of transcription factors, including PPARγ [1], stimulating the expression of gluconeogenic enzymes such as cytosolic phosphoenolpyruvate carboxykinase (PCK1). PCK1 catalyzes the reversible reaction that decarboxylates OAA (oxaloacetic acid) to obtain PEP (phosphoenolpyruvate). PEP is a precursor of several biomolecules, being a key regulatory enzyme of gluconeogenesis, glyceroneogenesis and other metabolic pathways [2]. We have recently reported the effects of pig Pck1 c.A2456C substitution on fat distribution in pigs [3]. The c.A2456C substitution produces a Met139Leu substitution in Pck1 that modifies the kinetic properties and protein stability in vivo and in vitro. This SNP is also associated, among others, with a decrease in backfat thickness (BT), a change in water-holding capacity and an increase in intramuscular fat content (IMF), traits related to fat metabolism and meat quality.

PGC1α regulates not only metabolic responses, it also contributes to mitochondrial biogenesis [4], interacts with proteins with histone acetyltransferase activity [5] and determines the type of muscle fiber [6]. In humans, the G482S variant of PGC1α has been shown to modify fat metabolism in hepatocytes [7]. Several polymorphisms have been described in pig Pgc1α [8–10]. The most studied one [8], c.T1378A, produces a single amino acid change from cysteine to serine at position 430 (Cys430Ser). The reports about its phenotypic consequences are conflicting, ranging from the finding of no effects [9,10] to changes in muscle pH values [11] or intramuscular fat content [12].

Given the role of PGC1α in regulating PCK1, we hypothesized that the Cys430Ser substitution in Pgc1α modulates PCK1 activity. We have studied whether post-translational modifications of PGC1α could influence expression of PCK1, focusing on differences in glycosylation between Pgc1α wild type (WT) and p.C430S proteins. We have found that although the WT variant is less stable than the p.C430S one, the former is responsible for higher expression levels of PCK1. Higher O-GlcNAcylation of p.C430S decreases its interaction with PPARγ, which explains the observed effect. Furthermore, this molecular mechanism can explain the potential epistatic effect observed between Pgc1α and Pck1.

Materials and methods

Cloning

Human DYK tagged O-GlcNAc transferase (OGT, OHu54262) and pig Pgc1α-myc (OSe00030D; p.C430S) clones were obtained from GenScript. Pgc1α was subcloned into pEGFP-N1 using BamHI and XhoI restriction enzymes. Site-directed mutagenesis was performed using Pfu Ultra II HS DNA polymerase (Agilent) and DpnI (Fermentas). The following primers were used to create p.430C variant in pig Pgc1α (PGC1aSer430Cys-F 5’-ccacagactcagaccagtgctacctgaccgagacgtcggag-3’ and PGC1aSer430Cys-R 5’-ctccgacgtctcggtcaggtagcactggtctgagtctgtgg-3’). The stop codon was eliminated to overexpress Pgc1α-myc-GFP using PGC1aEGFP-F 5-catctcagaagaggatctgttggatccaccggtcgccacc-3 and PGC1aEGFP-R 5-ggtggcgaccggtggatccaacagatcctcttctgagatg-3.

Quantitative PCR assays

To determine differences in Pgc1α activity, mRNA levels of endogenous PCK1 were measured. HEK293T cells (4x10⁴) were grown on 96-well plates using complete DMEM (DMEM supplemented with 10% fetal bovine serum, 2 mM L-glutamine, 0.1 mg/mL streptomycin and 100 U/mL penicillin). After 24 hours, cells were transfected using GeneJuice (Novagen) following manufacturer's instructions. 24 hours later, mRNA was extracted and cDNA was obtained using Cells to Ct 1-Step Taqman kit (Thermo Fisher Scientific). PCK1 cDNA was amplified using a human PCK1 probe set (Thermo Fisher Scientific) in an Mx3005P qPCR System (Agilent). GAPDH was used as control. PCK1 was quantified as ; where Δ_Ct is Pck1_Ct − GAPDH_Ct from each sample.

Endogenous PCK1 detection

HEK293T cells (ATCC, CRL-3216) (7.5x10⁵ per well) were grown on 6-well plates and treated as described above. HepG2 (1x10⁶) cells were maintained in MEM supplemented with 10% FBS, grown on 6-well plates and reverse transfected using GeneJuice: trypsinized cells were directly plated on the transfection mixture, prepared following manufacturer’s instructions. 24 hours post-transfection cells were lysed in PBS containing 1% NP-40 (Sigma) and a protease inhibitor cocktail (1 μM PMSF, 10 μM benzamidine and 0.5 μM leupeptin). Samples were boiled in SDS loading buffer and loaded into a 10% polyacrylamide gel. Proteins were transferred onto a PVDF membrane and blocked with PBS containing 0.5% Tween20 (Sigma) and 10% skimmed milk for 1 h. The membrane was incubated overnight at 4^°C with Anti-PCK1 (Abcam Ab28455, 1:500) and anti-actin (Sigma A2066 1:20000) and then washed and incubated with secondary HRP-conjugated anti-rabbit antibody (Millipore 12–348, 1:10000) for 1h at RT. Western blots were visualized by ECL (Millipore).

Glucose production assay

To analyze metabolic differences between the p.430C and p.430S variants, we measured the ability of HEK293T cells to produce glucose when overexpressing each variant. Glucose production was measured as previously described [13]. Cells were grown on 24-well plates. 24 hours after transfection, medium was removed, cells were washed once with PBS and 0.5 mL of glucose production medium was added (DMEM without glucose and phenol red supplemented with 2 mM sodium pyruvate and 20 mM sodium lactate). After 6 hours, half of the medium was collected and a colorimetric assay was performed (Sigma, GAGO20). Data were normalized to the total protein content from the whole cell lysates measured by Bradford assays.

Subcellular localization of Pgc1α variants

Transfected HEK293T cells were grown on 10 mm round cover slides in a 24-well dish plate and incubated in complete DMEM supplemented with 2 μM Hoechst 33342 (Invitrogen) for 25 minutes to stain nuclei. Then, cells were fixed in 3.7% formaldehyde at RT for 20 minutes, washed with PBS three times and once with water. Samples were mounted on microscope slides using ProLong Gold Antifade Mountant (Thermo Fisher) and observed using a DMI6000 Leica fluorescent microscope.

Protein stability in HEK293T cells

HEK293T cells (7.5x10⁵) were grown on 6-well plates and treated and transfected as described above. Medium was replaced with complete DMEM supplemented with 150 μg/mL cycloheximide (Sigma) to inhibit protein synthesis. Cells were collected at different time points (0, 5, 10, 20, 40 and 60 minutes) in PBS supplemented with the aforementioned protease inhibitor cocktail and 10 μM of proteasome inhibitor MG132 (Sigma) to avoid rapid degradation of Pgc1α. Cells were pelleted and lysed in SDS-PAGE loading buffer. Proteins were separated in 8% polyacrylamide gels and transferred onto a PVDF membrane. Blocked membranes were incubated with anti-myc (Invitrogen R950-25, 1:5000) and anti-actin and visualized as described previously. Western blots were analyzed using ImageJ (NIH).

Immunoprecipitation and co-immunoprecipitation

To determine differences in post-translational modifications (PTMs), each variant was immunoprecipitated from mammalian cells. HEK293T cells were grown on 6-well plates, co-transfected with or without OGT and either Pgc1α WT or Pgc1α p.C430S. When specified, cells were incubated in DMEM supplemented with PUGNAc at 100 μ3ϕoρ 18η to inhibit O-GlcNAcase. Cells were lysed in 50 mM TRIS, pH 8, 150 mM NaCl, 1% NP-40, 0.1 mM EDTA, 2 mM nicotinamide, 10 mM NaF and protease inhibitor cocktail. Clarified cell lysates were incubated with 5 μg of anti-myc (Invitrogen 9E10 clone) at 4^°C overnight. Pgc1α was immunoprecipitated using 25 μL of Dynabeads (Invitrogen) at 4^°C for 2 hours. Non-transfected cells were used as control. Samples were boiled in SDS-PAGE loading buffer, run in 8% gels and transferred to PVDF membranes. These membranes were incubated with anti-myc, anti-O-GlcNAc (Abcam, Ab2739, 1:5000), anti-acetyl-lysine (Cell Signal, 9441, 1:1000), anti-PhosphoSer/Thr (Abcam, Ab17464, 1:1000) and anti-PPARγ (Santa Cruz, sc-7196, 1:500). Clean blot HRP detection kit (Thermo Fisher Scientific) was used to detect PPARγ from immunoprecipitates. Band density was analyzed using ImageJ.

Statistical analysis

Data normality was analyzed using Shapiro-Wilk test. For glucose, qPCR and PCK1 western blotting, ANOVA tests were performed following Tukey's post-hoc tests. In the case of cycloheximide, immunoprecipitation and co-immunoprecipitation assays, Welch's two sample t-tests were performed.

Isolation of genomic DNA

Genomic DNA was extracted from tails of Duroc x Landrace/Large White piglets. DNA was extracted using the Realpure genomic DNA extraction kit (Durviz) following the manufacturer's instructions. DNA concentration in the extract was quantified using a NanoDrop 1000 spectrophotometer (Thermo Fisher Scientific).

Phenotypic recording and analysis of the Pgc1α c.T1378A and Pck1 c.A2456C genotypes

We used a total of 202 pigs that had been phenotypically characterized and also genotyped with respect to the c.A2456C SNP in Pck1 in a previous work [3]. DNA amplification was carried out with AccuPrime^TMTaq DNA Polymerase High Fidelity (Invitrogen) using the primers described by Kunej et al. [8] to amplify the Pgc1α SNP region. The PCR was performed in an UNO-Thermoblock apparatus (Biometra). The Pgc1α PCR product was purified using the NucleoSpin Extract II kit (Macherey-Nagel) following manufacturer′s instructions, quantified using a NanoDrop 1000 spectrophotometer and sequenced at Sistemas Genómicos (Paterna, Valencia, Spain) to identify the Pgc1α genotype.

Epistatic analysis

The statistical analysis for detection of additive and additive x additive epistatic analysis was performed with a linear model under the Bayesian framework. The model of analysis was:

y = wb + Xβ + Zu + Tp + λ_Pgc1αa_Pgc1α + λ_Pck1a_Pck1 + λ_axai_axa + e, where y is the vector of phenotypic data, b is the covariate with live weight at slaughter, β is the vector of batch effects, u is the vector of polygenic additive genetic effects, p is the vector of litter effects, a_Pgc1α, a_Pck1 and are the additive effects of Pgc1α, Pck1, respectively, i_axa is the additive x additive epistatic interaction effect and e is the vector of residual effects. Further, w is the vector of live weights at slaughter, X, Z and T are the corresponding incidence matrices that links systematic and polygenic additive effects with data, λ_Pgc1a, λ_Pck1 and λ_axa are the vectors of the additive and additive x additive coefficients calculated with the orthogonal approach proposed by Álvarez-Castro and Carlborg [14]. We assumed bounded uniform distributions for the covariate (b), systematic effects (β) and the additive and additive x additive effects (a_Pgc1α, a_Pck1 and i_axa). Further, the following multivariate prior Gaussian distributions were assumed for polygenic and litter effects: and , where A is the numerator relationship matrix calculated from the pedigree information and I is the identity matrix. Finally, prior distributions for variance components (, and ) were bounded uniform. The analysis was performed using a Gibbs Sampler [15] using a single chain of 100,000 iterations after discarding the first 10,000 for the calculation of the posterior distributions. Relevance of additive or additive x additive effects was determined when the posterior distribution had a probability over (or below) zero greater than 0.9 or below 0.1. Other epistatic interactions (additive x dominant, dominant x additive and dominant x dominant) were not presented because of the scarcity of available information (the low number of animals and underrepresented haplotypes) to obtain appropriate posterior estimates.

Results

Pgc1α C430S substitution decreases PCK1 mRNA and protein levels and glucose production in HEK293T cells

Cysteine 430 is conserved among species (Fig 1A). We examined in HEK293T cells if pig Pgc1α p.C430S substitution could affect the expression of endogenous human PCK1. Fig 1B shows that this was indeed the case and that Pgc1α WT induced higher levels (2.4 fold) of PCK1 transcription than Pgc1α p.C430S, as measured by mRNA quantification using qPCR (Fig 1B). The amount of PCK1 protein detected by Western blotting in these experiments was also significantly larger when overexpressing Pgc1α WT than Pgc1α p.C430S in HEK293T and HepG2 cells (Fig 1C and 1D), although only by 20%. To determine if these differences in PCK1 protein content had a translation into metabolic performance of HEK293T cells, we analyzed glucose production in cells overexpressing each Pgc1α variant. It has to be remembered here that PCK1 is one of the main regulators of gluconeogenesis. Overexpression of both Pgc1α variants substantially increased glucose production compared to controls transfected with empty vectors, but cells overexpressing Pgc1α WT synthetized 19% more glucose (9.3±0.3 μg glucose/mg total protein) than those overexpressing Pgc1α p.C430S (7.5±0.3 μg/mg) (Fig 1E). This difference in glucose production matches nicely the difference in PCK1 protein amount found when expressing both Pgc1α variants.

Download:

Fig 1. Cys430Ser substitution affects Pgc1α transcriptional activity.

(A) As can be seen in the alignment, pig cysteine 430 (red bold) is well conserved among species and it also lies in a serine/threonine rich region. (B) mRNA levels of endogenous PCK1. HEK293T cells were transfected with either Pgc1α p.C430S or wild type (WT) and after 24h total mRNA was extracted and quantified using qPCR in three individual assays, using o mock transfected cells. Pgc1α WT produces more PCK1 mRNA than Pgc1α p.C430S. Shapiro-Wilk test was performed to check data normality. This experiment was repeated three times with three technical replicates each. ANOVA test showed significant differences between groups (n = 3; p = 1.3x10^-5) and Tukey post-hoc test was performed to find differences among groups. (C) Pgc1α WT overexpression leads to a higher expression of PCK1 in HEK293T cells, as detected by Western blot. The experiment was repeated three times. The blots shown belong to a representative experiment. Bands were analyzed using ImageJ. Cells were treated and data was analyzed as described (n = 3; ANOVA, p = 0.0027). (D) Same experiment as (C) performed in HepG2 cells. (E) Cells overexpressing either Pgc1α p.C430S or WT were incubated in a free-glucose medium. After 6 hours, half of the medium was collected and glucose was determined using a colorimetric assay. Data were normalized to total cell protein content. Pgc1α p.C430S produces more glucose (19%) than Pgc1α WT. Three independent assays were carried out. Data was analyzed as described (n = 3; ANOVA, p = 3.3x10^-8). In all cases, control cells were those transfected with empty vector. Asterisks indicate * p< 0.05, ** p< 0.01 and *** p< 0.001. C430S: Pgc1α p.C430S and WT: Pgc1α wild type.

https://doi.org/10.1371/journal.pone.0179988.g001

Stability and subcellular localization of Pgc1α variants

The metabolic effects of any protein can be modified not only by increasing the amount of the protein or changing its functionality, but also by a change in protein stability. To investigate if the C430S substitution affected Pgc1α stability we analyzed Pgc1α content in HEK293T cells by Western blotting after cycloheximide addition. Cycloheximide blocks protein synthesis and allows following protein decay at different time points. Fig 2A shows that both Pgc1α variants significantly differ in stability: 60 minutes after cycloheximide addition, the amount of Pgc1α WT was reduced down to 26.7±10.8 (%) of the initial protein amount whereas that of Pgc1α p.C430S was only reduced down to 65.4±7.3 (%). This difference was statistically significant (n = 3, p = 0.0095). To rule out the possibility that the substitution could change the subcellular location of Pgc1α and modify in that way its stability, we analyzed the cellular location of GFP-tagged Pgc1α variants. As shown in Fig 2B the subcellular localization was not altered since both variants were localized in the nucleus.

Download:

Fig 2. Cys430Ser substitution in Pgc1α affects protein stability but not subcellular localization.

(A) Protein synthesis was inhibited with cycloheximide to look for differences in protein stability. As can be seen, Pgc1α WT has a shorter half-life. After 60 minutes, this difference is statistically significant (n = 3, p = 0.0095). Image shows a representative western blot of three independent experiments. (B) Subcellular localization of Pgc1α-GFP (green) in HEK293T cells. Hoechst 33342 was used to stain cell nuclei (blue). Overlays of images of green and blue channels are shown. Magnifications of red squares in upper panels are presented in lower panels. White line represents 10 μm. Transfected cells with empty vector show GFP in the cytosol and nucleus. In the case of both Pgc1α p.C430S and WT, they are located only in the nucleus. ** p< 0.01. C430S: Pgc1α p.C430S and WT: Pgc1α wild type.

https://doi.org/10.1371/journal.pone.0179988.g002

The effects of C430S substitution on post-translational modifications of Pgc1α

The larger stability of Pgc1α p.C430S could, at least in theory, counteract its lower capacity to stimulate PCK1 synthesis. We analyzed therefore additional biochemical properties of Pgc1α that could be affected by the substitution and that could help to clarify the interaction between Pgc1α and Pck1. As Pgc1α has been shown to undergo several types of post-translational modifications (PTMs), such as phosphorylation, acetylation or O-GlcNAcylation, we decided to study whether the C430S substitution affected the ability of Pgc1α to undergo these covalent modifications. To do that, we immunoprecipitated both Pgc1α variants with anti-myc antibody and analyzed by Western blot with different specific antibodies their level of phosphorylation, acetylation or O-GlcNAcylation (Fig 3A). No difference was found between Pgc1α variants; indeed, no signal at all of O-GlcNAc was detected. We were somehow surprised because we did not detect any signal of O-GlcNAcylation in our experiments as this modification has been convincingly shown to modify Pgc1α stability [16], the property we were looking for. This led us to follow two different approaches, on one hand cells were co-transfected with both Pgc1α and OGT (Fig 3B), the enzyme that catalyzes O-GlcNAcylation, and on the other hand cells were treated with PUGNAc (Fig 3C), an O-GlcNAcase inhibitor. Fig 3B and 3C show that both cotransfection of Pgc1α and OGT and PUGNAc treatment led to detection of O-GlcNAc in both Pgc1α variants and that the level of O-GlcNAcylation was significantly lower in Pgc1α WT than that in Pgc1α p.C430S.

Download:

Fig 3. Cys430Ser substitution modifies O-GlcNAcylation level of Pgc1α and affects the interaction with PPARγ.

(A) HEK293T cells overexpressing each Pgc1α variant were lysed and Pgc1α was immunoprecipitated using anti-myc. Afterwards, PTMs (glycosylation -O-GlcNAc-, lysine acetylation—AcK- and serine/threonine phosphorylation—P-Ser/Thr-) were studied and no differences were observed between the two Pgc1α isoforms (IP, immunoprecipitation; IB, immunoblotting). In fact, O-GlcNAcylation was too low to be detected under our conditions. (B) HEK293T cells were co-transfected with a DYK-tagged OGT-expressing plasmid and detected with an anti-DYK tag antibody (α-DYK). No differences are observed in the levels of acetylation or phosphorylation but Pgc1α p.C430S shows a higher level of O-GlcNAcylation than Pgc1α WT (n = 3, p = 0.0052). However, detection of PPARγ from immunoprecipitates shows that this protein has stronger interaction with Pgc1α WT (n = 3, p = 0.0079). These differences in O-GlcNAcylation and PPARγ interaction are graphically represented. ** p< 0.01. C430S: Pgc1α p.C430S; WT = Pgc1α wild type. (C) Cells were incubated in the presence of PUGNAc, a O-GlcNAcase inhibitor, for 16h and the levels of each Pgc1α variant were detected by western blotting with anti-myc and total protein O-GlcNAcylation was detected with anti-OGlc-NAc to show the effects of OGT inhibition. Then each Pgc1 variant was immunoprecipitated with anti-myc and detected with anti-myc and with anti-O-Glc-NAc. The lower panel shows a slightly lower level of O-GlcNAcylation of the WT protein with respect to p.C430S. Note that upon detection of immunoprecipitated Pgc1α with anti-myc a background signal is observed also in lanes corresponding to cells that do not express myc-tagged Pgc1α. This assay was carried out once with three technical replicates as an alternative to the experiment shown in Fig 3B.

https://doi.org/10.1371/journal.pone.0179988.g003

The effect of O-GlcNAc level on Pgc1α-PPARγ interaction

PGC1α interacts with PPARγ to stimulate PCK1 expression [17,18]. Since we had found that the C430S variant affected PCK1 transcription (Fig 1), we decided to examine if the different level of O-GlcNAc of 430S and 430C could influence the ability of these variants to interact with PPARγ and affect in that way the expression of PCK1. Overexpressed Pgc1α variants, showing different levels of O-GlcNAcylation (Fig 3B), were immunoprecipitated using an anti-myc antibody and PPARγ was analyzed in the immunoprecipitates using a specific antibody. Fig 3B shows that Pgc1α p.C430S interaction with PPARγ was significantly diminished in comparison to that of Pgc1α WT. Thus, it can be reasonably proposed that this fact will contribute to an increase in PCK1 expression by Pgc1α WT mediated by PPARγ.

Effects of PGC1α C430S substitution on phenotypic traits in pigs

We have previously reported [3], the influence of a PCK1 SNP (c.A2456C, M139L), in some phenotypic traits in pigs, and related it to differences in enzyme kinetics of both variants at the molecular level. This, together with the findings reported above, led us to analyze the influence of the C430S substitution in phenotypic traits in pigs and also a possible epistatic interaction between Pck1 and Pgc1α.

Table 1 shows the Bayesian analysis of the influence of PGC1α c.T1387A (C430S) on conformational and meat quality traits on a previously described pig population [3]. For comparative purposes we present also the published effects of Pck1 SNP on those traits in the same experimental group. S1 Table shows the same analysis for all studied traits. Genotypic frequencies for Pgc1α in this population were 29.8% AA, 56.1% AT and 14.1% TT. Only P values above 0.9 or below 0.1 were considered relevant for the additive and epistatic effects. None, but one, of the Bayesian P values of any traits were close to 0 or to 1 suggesting that the SNP in Pgc1α has no direct influence on those traits (Table 1).

Download:

Table 1. Posterior mean estimate and Bayesian posterior probability of the additive and additive x additive epistatic interaction effect between Pgc1α c.T1378A and Pck1 c.A2456C.

https://doi.org/10.1371/journal.pone.0179988.t001

Epistatic effects between Pgc1α and Pck1

The Bayesian statistical analysis also detects epistatic interactions between Pgc1α and Pck1. Table 1 also shows the results of this analysis on conformational and meat quality traits as in the former section; it shows only epistatic effects with a posterior probability above 0.9 or below 0.1. (S1 Table shows the lack of appreciable epistatic effects on compositional traits). The interaction between Pgc1α and Pck1 (see Table 1) showed a high number (11/19) of these above-threshold interactions, that is almost 60% of total epistatic events. These interactions have some internal coherence: the three exudation traits show the same trend, being increased by the combination of Pck1 c.2456C and Pgc1α c.1378A alleles whereas the same combination decreased the adipose tissue content for four traits (backfat thickness, belly, lard and adipose tissue content in ham). Fig 4 shows graphical representations of these epistatic effects for two traits: variation in the percentage of ham and variation in backfat thickness. In both cases, there is a clear additive effect associated with Pcg1α for genotypes AA and AC at Pck1 locus, but it disappears (backfat thickness) or is diminished (ham content) for individuals with genotype CC.

Download:

Fig 4. Epistasis between Pg1α c.A1378T and Pck1 c.A2456C polymorphisms.

Graphical representation of epistatic effects from data presented in Table 1 for traits (A) backfat thickness and (B) ham content. The additive effect of locus Pgc1α is greater in the Pck1_AA background and disappears (backfat thickness variation) or is diminished (ham content variation) in the Pck1_CC background.

https://doi.org/10.1371/journal.pone.0179988.g004

Discussion

PGC1α regulates the expression of gluconeogenic enzymes [4,5,13] such as PCK1, an enzyme which is in a crossroad of several metabolic pathways [2]. PCK1 expression is regulated, among others, by PPARγ [17] which is activated by PGC1α [1]. We report here how alterations in the interaction between the transcription factors PGC1α and PPARγ caused by different levels of post-translational glycosylation of genetic variants of the former might be responsible for differences in phenotypic traits related to PCK1 activity. Furthermore, these molecular mechanisms can explain, at least in part, the apparent epistatic effects observed between PGC1α and PCK1.

Pig Pgc1α has 95% identity and 98.2% similarity with human PGC1α, which justified using a gluconeogenic human cell line with high transfection efficiency, HEK293T, to study at the molecular level the genetic interaction between PGC1α and PCK1, focusing on a common Pgc1α polymorphism found in pigs, a C430S substitution. PCK1 mRNA levels strongly increased when cells overexpressed each Pgc1α variant but that effect was higher (2.4 fold) when overexpressing Pgc1α WT (Fig 1B). In the same way, when overexpressing Pgc1α WT, cells reached higher PCK1 protein levels (20%) (Fig 1C and 1D) and produced more glucose (19%) (Fig 1E). Therefore, our results indicate that these two Pgc1α variants have different transcriptional activities with respect to PCK1 expression.

In order to identify the molecular mechanism underlying the differential PCK1 expression mediated by Pgc1α variants, we first analyzed their stability in HEK293T cells. To our surprise, the most active variant with respect to PCK1 expression, Pgc1α WT, was less stable than Pgc1α p.C430S. Zhang et al. [19] described the presence of a splicing isoform of PGC1α (NT-PGC1α) that was mainly present in the cytosol and had increased half-life in the cell [20]. Altered subcellular localization could therefore explain the differences we observed. However, our results showed that this was not the case since both variants were clearly located in the nucleus (Fig 2B).

Since the ability of Pgc1α WT to stimulate PCK1 synthesis was higher in spite of its lower stability, we decided to analyze whether differences in post-translational modifications (PTMs) could be responsible for the observed effects. We focused on the most common PTMs in PGC1α: acetylation, phosphorylation and O-GlcNAcylation [21]. Under standard culture conditions, no differences were found for any of these PTMs (Fig 3A) between Pgc1α variants. Furthermore, no signal was detected in the case of O-GlcNAcylation. Taking into account the difficulties in detecting O-GlcNAcylation in PGC1α [16], we co-transfected HEK293T cells with human O-GlcNAc transferase (OGT) with either Pgc1α p.C430S or Pgc1α WT to ensure more OGT activity in the cells. On the other hand, transfected HEK293T cells were treated with PUGNAc, an O-GlcNAcase inhibitor, as previously described [16]. Under these conditions, neither acetylation nor phosphorylation were affected. Nevertheless, Pgc1α p.C430S showed not only higher levels of O-GlcNAcylation (Fig 3B and 3C), but also its interaction with PPARγ was considerably weakened (Fig 3B). Acetylation of some residues has been described to activate [22] or to inactivate [23] PGC1α. The same occurs with phosphorylation [24,25]. In both cases, the detected signals were weak. In the case of O-GlcNAcylation, it has been described to stabilize and to increase PGC1α activity via OGT/HIF1 complex [16]. If this is true, the increased O-GlcNAcylation of p.C430S should contribute to the increased stability of this protein (Fig 2A) and, thus, should facilitate the interaction of this variant of Pgc1α with other proteins. However, we found exactly the opposite: a significantly weaker interaction of Pgc1α p.C430S with PPARγ than that for the WT variant. As this interaction is the one that ultimately leads to increased PCK1 expression, we could reasonably conclude that increased levels of O-GlcNAcylation impair the PPARγ-Pgc1α interaction and explain the higher PCK1 expression upon overexpression of Pgc1α WT (Fig 1B).

The existence of gene variability defines phenotypic differences between individuals. In most cases, the individual effects of these genes are subtle but not negligible. Quantitative traits play a relevant role in animal [26] and plant [27] breeding, tissue development [28] and disease [29]. Part of these effects can be a consequence of the existence of epistatic interactions [30]. The importance of Pck1 for pig breeders has been highlighted by our recent report: the existence of a SNP in the coding region of this gene which resulted in the expression of an enzyme with substantially modified biochemical properties associated with better meat quality traits [3]. Given the importance of quantitative traits and the existence of variability in both pig Pgc1α and Pck1, we studied the epistatic interactions between these two genes and the potential mechanisms underlying this effect, focusing on the Pgc1α c.T1378A (p.C430S) polymorphism described in this work and on the Pck1 c.A2456C one (p.M139L) described in our former work [3].

We first evaluated Pgc1α additive effects on a Duroc x Landrace/Large White population and we found relevant association in just 6 out of 55 traits (noble parts yield, tenderloin content, protein content in loin, protein content in B. femoris, moisture in B. femoris and moisture in ham; S1 Table). This is even less than what could be expected just by chance. Liu et al. [12] described the effects of pig Pgc1α c.T1378A substitution on a Tibetan breed population and found significant effects on IMF content and muscle fiber size. The effects on IMF can be related to the expression of metabolic genes such as Pck1. However, although PGC1α is involved in defining muscle fiber type [6], it has no effects on fiber size in mice [31]. The fact that we did not find differences in similar traits in our pig population may be either a consequence of the low frequency of the Pgc1α TT genotype (14.1%) in our pigs compared to that of the Tibetan breeds studied (63.9%) by Liu et al [12] or to breed-specific effects as has been suggested by other authors [32].

The lack of enough animals representing all the possible haplotypes precluded to extend the study to additive x dominant, dominant x additive and dominant x dominant interactions. Nevertheless, our epistatic analysis of additive x additive interactions under the Bayesian approach showed a total of almost 60% of total epistatic interactions with Bayesian posterior probabilities above 0.9 or below 0.1 (Table 1). This percentage should be compared with the expected 20% of Bayesian posterior probabilities over 0.9 or below 0.1 that can be expected just by chance. Moreover, these relevant posterior probabilities are found in traits mostly related to percentage of adipose tissue content and meat water holding capacity. These traits were also the most affected by the Pck1 polymorphism [3]. The combination of Pgc1α p.C430S and Pck1 p.M139L was associated to diverse effects. On one side, Pgc1α p.C430S reinforced the effects of Pck1 p.M139L decreasing meat quality, that is, yielding higher drip loss and exudation. On the other side, this very same allele combination is associated to lower backfat thickness and lower percentages of cuts rich in adipose tissue (lard and belly). The latter would be more convenient from an economical point of view because these cuts are less valued. Therefore, there is no single combination of Pgc1α and Pck1 alleles with only positive effects on both meat quality and carcass conformation traits as happened with Pck1 [3]. This comes to no surprise because PGC1α could produce its effects not only through PCK1 but through many other protein pathways [4,5]. Thus, the results of the epistatic analysis indicate that the application of these polymorphisms in marker-assisted selection must be taken with caution because, 1) there is not a clear optimum genotypic configuration and 2) it is plausible that the detected epistatic interaction is part of a more complex epistatic network with unpredictable results with the available information.

In conclusion, we report that differential O-GlcNAcylation of Pgc1α regulates PCK1 activity and this molecular mechanism could explain at least in part the epistatic interaction between both genes. This should be taken into account if these genes are used for animal selection in pig breeding.

Supporting information

S1 Table. Posterior mean estimate and Bayesian posterior probability of all analyzed traits.

https://doi.org/10.1371/journal.pone.0179988.s001

(XLSX)

References

1. 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. pmid:9529258
- View Article
- PubMed/NCBI
- Google Scholar
2. Yang J, Kalhan SC, Hanson RW. What is the metabolic role of phosphoenolpyruvate carboxykinase? Journal of Biological Chemistry. 2009. pp. 27025–27029. pmid:19636077
- View Article
- PubMed/NCBI
- Google Scholar
3. Latorre P, Burgos C, Hidalgo J, Varona L, Carrodeguas JA, López-Buesa P. c.A2456C-substitution in Pck1 changes the enzyme kinetic and functional properties modifying fat distribution in pigs. Sci Rep. 2016;6: 19617. pmid:26792594
- View Article
- PubMed/NCBI
- Google Scholar
4. Fernandez-Marcos PJ, Auwerx J. Regulation of PGC-1α, a nodal regulator of mitochondrial biogenesis. American Journal of Clinical Nutrition. 2011. pmid:21289221
- View Article
- PubMed/NCBI
- Google Scholar
5. Finck BN, Kelly DP. PGC-1 coactivators: Inducible regulators of energy metabolism in health and disease. Journal of Clinical Investigation. 2006. pp. 615–622. pmid:16511594
- View Article
- PubMed/NCBI
- Google Scholar
6. Lin J, Wu H, Tarr PT, Zhang C-Y, Wu Z, Boss O, et al. Transcriptional co-activator PGC-1 alpha drives the formation of slow-twitch muscle fibres. Nature. 2002;418: 797–801. pmid:12181572
- View Article
- PubMed/NCBI
- Google Scholar
7. Chen Y, Mu P, He S, Tang X, Guo X, Li H, et al. Gly482Ser mutation impairs the effects of peroxisome proliferator-activated receptorα coactivator-1α on decreasing fat deposition and stimulating phosphoenolpyruvate carboxykinase expression in hepatocytes. Nutr Res. 2013;33: 332–339. pmid:23602251
- View Article
- PubMed/NCBI
- Google Scholar
8. Kunej T, Wu XL, Berlic TM, Michal JJ, Jiang Z, Dovc P. Frequency distribution of a Cys430Ser polymorphism in peroxisome proliferator-activated receptor-gamma coactivator-1 (PPARGC1) gene sequence in Chinese and Western pig breeds. J Anim Breed Genet. 2005;122: 7–11. pmid:16130483
- View Article
- PubMed/NCBI
- Google Scholar
9. Erkens T, Rohrer GA, Van Zeveren A, Peelman LJ. SNP detection in the porcine PPARGC1A promoter region and 3’UTR, and an association study in a Landrace-Duroc-Large White population. Czech J Anim Sci. 2009;54: 408–416.
- View Article
- Google Scholar
10. Erkens T, De Smet S, Van den Maagdenberg K, Stinckens A, Buys N, Van Zeveren A, et al. Association analysis of PPARGC1A mutations with meat quality parameters in a commercial hybrid pig population. Czech J Anim Sci. 2010;55: 200–208.
- View Article
- Google Scholar
11. Gandolfi G, Cinar MU, Ponsuksili S, Wimmers K, Tesfaye D, Looft C, et al. Association of PPARGC1A and CAPNS1 gene polymorphisms and expression with meat quality traits in pigs. Meat Sci. 2011;54: 478–485. pmid:21680104
- View Article
- PubMed/NCBI
- Google Scholar
12. Liu R, Li J, Lv X. Association of PGC-1α Gene with Intramuscular Fat Content and Muscle Fiber Traits and Gene Expression in Tibetan Pigs. J Anim Vet Adv. 2011;10: 2301–2304.
- View Article
- Google Scholar
13. Yoon JC, Puigserver P, Chen G, Donovan J, Wu Z, Rhee J, et al. Control of hepatic gluconeogenesis through the transcriptional coactivator PGC-1. Nature. 2001;413: 131–138. pmid:11557972
- View Article
- PubMed/NCBI
- Google Scholar
14. Alvarez-Castro JM, Carlborg O. A Unified Model for Functional and Statistical Epistasis and Its Application in Quantitative Trait Loci Analysis. Genetics. 2007;176: 1151–1167. pmid:17409082
- View Article
- PubMed/NCBI
- Google Scholar
15. Gelfand A, Smith AFM. Sampling-Based Approaches to Calculating Marginal Densities. J Am Stat Assoc. 1990;85: 398–409.
- View Article
- Google Scholar
16. Bin Ruan H, Han X, Li MD, Singh JP, Qian K, Azarhoush S, et al. O-GlcNAc transferase/host cell factor C1 complex regulates gluconeogenesis by modulating PGC-1α stability. Cell Metab. 2012;16: 226–237. pmid:22883232
- View Article
- PubMed/NCBI
- Google Scholar
17. Tontonoz P, Hu E, Devine J, Beale EG, Spiegelman BM. PPAR gamma 2 regulates adipose expression of the phosphoenolpyruvate carboxykinase gene. Mol Cell Biol. 1995;15: 351–357. pmid:7799943
- View Article
- PubMed/NCBI
- Google Scholar
18. Glorian M, Duplus E, Beale EG, Scott DK, Granner DK, Forest C. A single element in the phosphoenolpyruvate carboxykinase gene mediates thiazolidinedione action specifically in adipocytes. Biochimie. 2001;83: 933–943. pmid:11728630
- View Article
- PubMed/NCBI
- Google Scholar
19. Zhang Y, Huypens P, Adamson AW, Chang JS, Henagan TM, Boudreau A, et al. Alternative mRNA splicing produces a novel biologically active short isoform of PGC-1α. J Biol Chem. 2009;284: 32813–32826. pmid:19773550
- View Article
- PubMed/NCBI
- Google Scholar
20. Trausch-Azar J, Leone TC, Kelly DP, Schwartz AL. Ubiquitin Proteasome-dependent Degradation of the Transcriptional Coactivator PGC-1α via the N-terminal Pathway. J Biol Chem. 2010;285: 40192–40200. pmid:20713359
- View Article
- PubMed/NCBI
- Google Scholar
21. Housley MP, Udeshi ND, Rodgers JT, Shabanowitz J, Puigserver P, Hunt DF, et al. A PGC-1alpha-O-GlcNAc transferase complex regulates FoxO transcription factor activity in response to glucose. J Biol Chem. 2009;284: 5148–5157. pmid:19103600
- View Article
- PubMed/NCBI
- Google Scholar
22. Higashida K, Kim SH, Jung SR, Asaka M, Holloszy JO, Han DH. Effects of Resveratrol and SIRT1 on PGC-1alpha Activity and Mitochondrial Biogenesis: A Reevaluation. PLoS Biol. 2013;11: e1001603. pmid:23874150
- View Article
- PubMed/NCBI
- Google Scholar
23. Lagouge M, Argmann C, Gerhart-Hines Z, Meziane H, Lerin C, Daussin F, et al. Resveratrol Improves Mitochondrial Function and Protects against Metabolic Disease by Activating SIRT1 and PGC-1α. Cell. 2006;127: 1109–1122. pmid:17112576
- View Article
- PubMed/NCBI
- Google Scholar
24. Jäger SS, Handschin CC, St-Pierre JJ, Spiegelman BMBM. AMP-activated protein kinase (AMPK) action in skeletal muscle via direct phosphorylation of PGC-1alpha. Pnas. 2007;104: 12017–12022. pmid:17609368
- View Article
- PubMed/NCBI
- Google Scholar
25. Li X, Monks B, Ge Q, Birnbaum MJ. Akt/PKB regulates hepatic metabolism by directly inhibiting PGC-1alpha transcription coactivator. Nature. 2007;447: 1012–1016. pmid:17554339
- View Article
- PubMed/NCBI
- Google Scholar
26. Dekkers J. Application of Genomics Tools to Animal Breeding. Curr Genomics. 2012;13: 207–212. pmid:23115522
- View Article
- PubMed/NCBI
- Google Scholar
27. St.Clair DA. Quantitative Disease Resistance and Quantitative Resistance Loci in Breeding. Annu Rev Phytopathol. 2010;48: 247–268. pmid:19400646
- View Article
- PubMed/NCBI
- Google Scholar
28. Wood AR, Esko T, Yang J, Vedantam S, Pers TH, Gustafsson S, et al. Defining the role of common variation in the genomic and biological architecture of adult human height. Nat Genet. 2014;46: 1173–86. pmid:25282103
- View Article
- PubMed/NCBI
- Google Scholar
29. Plomin R, Haworth CM a, Davis OSP. Quantitative Traits [Internet]. Genetics. 2009.
- View Article
- Google Scholar
30. Mackay TFC. Epistasis and quantitative traits: using model organisms to study gene-gene interactions. Nat Rev Genet. 2013;15: 22–33. pmid:24296533
- View Article
- PubMed/NCBI
- Google Scholar
31. Takikita S, Schreiner C, Baum R, Xie T, Ralston E, Plotz PH, et al. Fiber type conversion by PGC-1α. activates lysosomal and autophagosomal biogenesis in both unaffected and pompe skeletal muscle. PLoS One. 2010;5. pmid:21179212
- View Article
- PubMed/NCBI
- Google Scholar
32. Wood JD, Nute GR, Richardson RI, Whittington FM, Southwood O, Plastow G, et al. Effects of breed, diet and muscle on fat deposition and eating quality in pigs. Meat Sci. 2004;67: 651–667. pmid:22061815
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. 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. pmid:9529258
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Yang J, Kalhan SC, Hanson RW. What is the metabolic role of phosphoenolpyruvate carboxykinase? Journal of Biological Chemistry. 2009. pp. 27025–27029. pmid:19636077
View Article
PubMed/NCBI
Google Scholar

[6] View Article

[7] PubMed/NCBI

[8] Google Scholar

[ref3] 3. Latorre P, Burgos C, Hidalgo J, Varona L, Carrodeguas JA, López-Buesa P. c.A2456C-substitution in Pck1 changes the enzyme kinetic and functional properties modifying fat distribution in pigs. Sci Rep. 2016;6: 19617. pmid:26792594
View Article
PubMed/NCBI
Google Scholar

[10] View Article

[11] PubMed/NCBI

[12] Google Scholar

[ref4] 4. Fernandez-Marcos PJ, Auwerx J. Regulation of PGC-1α, a nodal regulator of mitochondrial biogenesis. American Journal of Clinical Nutrition. 2011. pmid:21289221
View Article
PubMed/NCBI
Google Scholar

[14] View Article

[15] PubMed/NCBI

[16] Google Scholar

[ref5] 5. Finck BN, Kelly DP. PGC-1 coactivators: Inducible regulators of energy metabolism in health and disease. Journal of Clinical Investigation. 2006. pp. 615–622. pmid:16511594
View Article
PubMed/NCBI
Google Scholar

[18] View Article

[19] PubMed/NCBI

[20] Google Scholar

[ref6] 6. Lin J, Wu H, Tarr PT, Zhang C-Y, Wu Z, Boss O, et al. Transcriptional co-activator PGC-1 alpha drives the formation of slow-twitch muscle fibres. Nature. 2002;418: 797–801. pmid:12181572
View Article
PubMed/NCBI
Google Scholar

[22] View Article

[23] PubMed/NCBI

[24] Google Scholar

[ref7] 7. Chen Y, Mu P, He S, Tang X, Guo X, Li H, et al. Gly482Ser mutation impairs the effects of peroxisome proliferator-activated receptorα coactivator-1α on decreasing fat deposition and stimulating phosphoenolpyruvate carboxykinase expression in hepatocytes. Nutr Res. 2013;33: 332–339. pmid:23602251
View Article
PubMed/NCBI
Google Scholar

[26] View Article

[27] PubMed/NCBI

[28] Google Scholar

[ref8] 8. Kunej T, Wu XL, Berlic TM, Michal JJ, Jiang Z, Dovc P. Frequency distribution of a Cys430Ser polymorphism in peroxisome proliferator-activated receptor-gamma coactivator-1 (PPARGC1) gene sequence in Chinese and Western pig breeds. J Anim Breed Genet. 2005;122: 7–11. pmid:16130483
View Article
PubMed/NCBI
Google Scholar

[30] View Article

[31] PubMed/NCBI

[32] Google Scholar

[ref9] 9. Erkens T, Rohrer GA, Van Zeveren A, Peelman LJ. SNP detection in the porcine PPARGC1A promoter region and 3’UTR, and an association study in a Landrace-Duroc-Large White population. Czech J Anim Sci. 2009;54: 408–416.
View Article
Google Scholar

[34] View Article

[35] Google Scholar

[ref10] 10. Erkens T, De Smet S, Van den Maagdenberg K, Stinckens A, Buys N, Van Zeveren A, et al. Association analysis of PPARGC1A mutations with meat quality parameters in a commercial hybrid pig population. Czech J Anim Sci. 2010;55: 200–208.
View Article
Google Scholar

[37] View Article

[38] Google Scholar

[ref11] 11. Gandolfi G, Cinar MU, Ponsuksili S, Wimmers K, Tesfaye D, Looft C, et al. Association of PPARGC1A and CAPNS1 gene polymorphisms and expression with meat quality traits in pigs. Meat Sci. 2011;54: 478–485. pmid:21680104
View Article
PubMed/NCBI
Google Scholar

[40] View Article

[41] PubMed/NCBI

[42] Google Scholar

[ref12] 12. Liu R, Li J, Lv X. Association of PGC-1α Gene with Intramuscular Fat Content and Muscle Fiber Traits and Gene Expression in Tibetan Pigs. J Anim Vet Adv. 2011;10: 2301–2304.
View Article
Google Scholar

[44] View Article

[45] Google Scholar

[ref13] 13. Yoon JC, Puigserver P, Chen G, Donovan J, Wu Z, Rhee J, et al. Control of hepatic gluconeogenesis through the transcriptional coactivator PGC-1. Nature. 2001;413: 131–138. pmid:11557972
View Article
PubMed/NCBI
Google Scholar

[47] View Article

[48] PubMed/NCBI

[49] Google Scholar

[ref14] 14. Alvarez-Castro JM, Carlborg O. A Unified Model for Functional and Statistical Epistasis and Its Application in Quantitative Trait Loci Analysis. Genetics. 2007;176: 1151–1167. pmid:17409082
View Article
PubMed/NCBI
Google Scholar

[51] View Article

[52] PubMed/NCBI

[53] Google Scholar

[ref15] 15. Gelfand A, Smith AFM. Sampling-Based Approaches to Calculating Marginal Densities. J Am Stat Assoc. 1990;85: 398–409.
View Article
Google Scholar

[55] View Article

[56] Google Scholar

[ref16] 16. Bin Ruan H, Han X, Li MD, Singh JP, Qian K, Azarhoush S, et al. O-GlcNAc transferase/host cell factor C1 complex regulates gluconeogenesis by modulating PGC-1α stability. Cell Metab. 2012;16: 226–237. pmid:22883232
View Article
PubMed/NCBI
Google Scholar

[58] View Article

[59] PubMed/NCBI

[60] Google Scholar

[ref17] 17. Tontonoz P, Hu E, Devine J, Beale EG, Spiegelman BM. PPAR gamma 2 regulates adipose expression of the phosphoenolpyruvate carboxykinase gene. Mol Cell Biol. 1995;15: 351–357. pmid:7799943
View Article
PubMed/NCBI
Google Scholar

[62] View Article

[63] PubMed/NCBI

[64] Google Scholar

[ref18] 18. Glorian M, Duplus E, Beale EG, Scott DK, Granner DK, Forest C. A single element in the phosphoenolpyruvate carboxykinase gene mediates thiazolidinedione action specifically in adipocytes. Biochimie. 2001;83: 933–943. pmid:11728630
View Article
PubMed/NCBI
Google Scholar

[66] View Article

[67] PubMed/NCBI

[68] Google Scholar

[ref19] 19. Zhang Y, Huypens P, Adamson AW, Chang JS, Henagan TM, Boudreau A, et al. Alternative mRNA splicing produces a novel biologically active short isoform of PGC-1α. J Biol Chem. 2009;284: 32813–32826. pmid:19773550
View Article
PubMed/NCBI
Google Scholar

[70] View Article

[71] PubMed/NCBI

[72] Google Scholar

[ref20] 20. Trausch-Azar J, Leone TC, Kelly DP, Schwartz AL. Ubiquitin Proteasome-dependent Degradation of the Transcriptional Coactivator PGC-1α via the N-terminal Pathway. J Biol Chem. 2010;285: 40192–40200. pmid:20713359
View Article
PubMed/NCBI
Google Scholar

[74] View Article

[75] PubMed/NCBI

[76] Google Scholar

[ref21] 21. Housley MP, Udeshi ND, Rodgers JT, Shabanowitz J, Puigserver P, Hunt DF, et al. A PGC-1alpha-O-GlcNAc transferase complex regulates FoxO transcription factor activity in response to glucose. J Biol Chem. 2009;284: 5148–5157. pmid:19103600
View Article
PubMed/NCBI
Google Scholar

[78] View Article

[79] PubMed/NCBI

[80] Google Scholar

[ref22] 22. Higashida K, Kim SH, Jung SR, Asaka M, Holloszy JO, Han DH. Effects of Resveratrol and SIRT1 on PGC-1alpha Activity and Mitochondrial Biogenesis: A Reevaluation. PLoS Biol. 2013;11: e1001603. pmid:23874150
View Article
PubMed/NCBI
Google Scholar

[82] View Article

[83] PubMed/NCBI

[84] Google Scholar

[ref23] 23. Lagouge M, Argmann C, Gerhart-Hines Z, Meziane H, Lerin C, Daussin F, et al. Resveratrol Improves Mitochondrial Function and Protects against Metabolic Disease by Activating SIRT1 and PGC-1α. Cell. 2006;127: 1109–1122. pmid:17112576
View Article
PubMed/NCBI
Google Scholar

[86] View Article

[87] PubMed/NCBI

[88] Google Scholar

[ref24] 24. Jäger SS, Handschin CC, St-Pierre JJ, Spiegelman BMBM. AMP-activated protein kinase (AMPK) action in skeletal muscle via direct phosphorylation of PGC-1alpha. Pnas. 2007;104: 12017–12022. pmid:17609368
View Article
PubMed/NCBI
Google Scholar

[90] View Article

[91] PubMed/NCBI

[92] Google Scholar

[ref25] 25. Li X, Monks B, Ge Q, Birnbaum MJ. Akt/PKB regulates hepatic metabolism by directly inhibiting PGC-1alpha transcription coactivator. Nature. 2007;447: 1012–1016. pmid:17554339
View Article
PubMed/NCBI
Google Scholar

[94] View Article

[95] PubMed/NCBI

[96] Google Scholar

[ref26] 26. Dekkers J. Application of Genomics Tools to Animal Breeding. Curr Genomics. 2012;13: 207–212. pmid:23115522
View Article
PubMed/NCBI
Google Scholar

[98] View Article

[99] PubMed/NCBI

[100] Google Scholar

[ref27] 27. St.Clair DA. Quantitative Disease Resistance and Quantitative Resistance Loci in Breeding. Annu Rev Phytopathol. 2010;48: 247–268. pmid:19400646
View Article
PubMed/NCBI
Google Scholar

[102] View Article

[103] PubMed/NCBI

[104] Google Scholar

[ref28] 28. Wood AR, Esko T, Yang J, Vedantam S, Pers TH, Gustafsson S, et al. Defining the role of common variation in the genomic and biological architecture of adult human height. Nat Genet. 2014;46: 1173–86. pmid:25282103
View Article
PubMed/NCBI
Google Scholar

[106] View Article

[107] PubMed/NCBI

[108] Google Scholar

[ref29] 29. Plomin R, Haworth CM a, Davis OSP. Quantitative Traits [Internet]. Genetics. 2009.
View Article
Google Scholar

[110] View Article

[111] Google Scholar

[ref30] 30. Mackay TFC. Epistasis and quantitative traits: using model organisms to study gene-gene interactions. Nat Rev Genet. 2013;15: 22–33. pmid:24296533
View Article
PubMed/NCBI
Google Scholar

[113] View Article

[114] PubMed/NCBI

[115] Google Scholar

[ref31] 31. Takikita S, Schreiner C, Baum R, Xie T, Ralston E, Plotz PH, et al. Fiber type conversion by PGC-1α. activates lysosomal and autophagosomal biogenesis in both unaffected and pompe skeletal muscle. PLoS One. 2010;5. pmid:21179212
View Article
PubMed/NCBI
Google Scholar

[117] View Article

[118] PubMed/NCBI

[119] Google Scholar

[ref32] 32. Wood JD, Nute GR, Richardson RI, Whittington FM, Southwood O, Plastow G, et al. Effects of breed, diet and muscle on fat deposition and eating quality in pigs. Meat Sci. 2004;67: 651–667. pmid:22061815
View Article
PubMed/NCBI
Google Scholar

[121] View Article

[122] PubMed/NCBI

[123] Google Scholar

O-GlcNAcylation mediates the control of cytosolic phosphoenolpyruvate carboxykinase activity via Pgc1α

O-GlcNAcylation mediates the control of cytosolic phosphoenolpyruvate carboxykinase activity via Pgc1α

Corrections

Figures

Abstract

Introduction

Materials and methods

Cloning

Quantitative PCR assays

Endogenous PCK1 detection

Glucose production assay

Subcellular localization of Pgc1α variants

Protein stability in HEK293T cells

Immunoprecipitation and co-immunoprecipitation

Statistical analysis

Isolation of genomic DNA

Phenotypic recording and analysis of the Pgc1α c.T1378A and Pck1 c.A2456C genotypes

Epistatic analysis

Results

Pgc1α C430S substitution decreases PCK1 mRNA and protein levels and glucose production in HEK293T cells

Stability and subcellular localization of Pgc1α variants

The effects of C430S substitution on post-translational modifications of Pgc1α

The effect of O-GlcNAc level on Pgc1α-PPARγ interaction

Effects of PGC1α C430S substitution on phenotypic traits in pigs

Epistatic effects between Pgc1α and Pck1

Discussion

Supporting information

S1 Table. Posterior mean estimate and Bayesian posterior probability of all analyzed traits.

References