A single nucleotide polymorphism in the coding region of PGC-1α is a male-specific modifier of Huntington disease age-at-onset in a large European cohort
© Weydt et al.; licensee BioMed Central Ltd. 2014
Received: 8 September 2013
Accepted: 10 December 2013
Published: 2 January 2014
Genetic modifiers are important clues for the identification of therapeutic targets in neurodegenerative diseases. Huntington disease (HD) is one of the most common autosomal dominant inherited neurodegenerative diseases. The clinical symptoms include motor abnormalities, cognitive decline and behavioral disturbances. Symptom onset is typically between 40 and 50 years of age, but can vary by several decades in extreme cases and this is in part determined by modifying genetic factors. The metabolic master regulator PGC-1α, coded by the PPARGC1A gene, coordinates cellular respiration and was shown to play a role in neurodegenerative diseases, including HD.
Using a candidate gene approach we analyzed a large European cohort (n = 1706) from the REGISTRY study for associations between PPARGC1A genotype and age at onset (AO) in HD.
We report that a coding variant (rs3736265) in PPARGC1A is associated with an earlier motor AO in men but not women carrying the HD mutation.
These results further strengthen the evidence for a role of PGC-1α in HD and unexpectedly suggest a gender effect.
KeywordsGenetic modifier Gender effect Neurodegeneration Huntington disease
Huntington disease (HD) is an autosomal-dominantly inherited, neurodegenerative disorder caused by the abnormal expansion of a CAG repeat tract in the huntingtin gene. The clinical picture of HD is characterized by a combination of a progressive movement disorder (chorea), cognitive decline and behavioral disturbances typically presenting in mid-life. However, disease onset and symptoms vary widely between individual patients. The length of the CAG repeat in the huntingtin gene and the resulting polyglutamine tract in the huntingtin protein is the major determinant of age at motor symptom onset, accounting for about 70% of the variability . The remaining variability is shaped by genetic and environmental factors that remain incompletely understood . Since genetic and environmental modifiers may provide clues for the identification of therapeutic targets in HD and possibly other neurodegenerative diseases, substantial research efforts are being invested into their discovery and characterization . Genetic modifiers can be identified with two alternative general methodologies: the hypothesis-driven candidate gene approach and unbiased genomewide screens.
The transcriptional co-activator PGC-1α is a metabolic master regulator that orchestrates the adaptation of the cell to changing energy demands . Originally described as regulator of mitochondrial respiration in brown adipose tissue , PGC-1α is now known to participate in many cell types in a wide range of essential and ancilliary metabolic processes, such as angiogenesis, antioxidative defense and autophagy. The human PPARGC1A gene was cloned in 1999 and subsequently found to be a genetic modifier in several metabolic diseases [10, 11]. In addition to HD, PGC-1α was also implicated as an important modifier of other neurodegenerative diseases, including Parkinson disease and amyotrophic lateral sclerosis (ALS) [12, 13]. The role of PGC-1α in the pathophysiology of neurodegenerative diseases is likely to be complex and differs between disease entities . A route towards resolving this complexity is the detailed characterization of the PPARGC1A locus, its variant transcripts and SNPs that are associated with disease onset. This strategy has resulted in the identification of a novel promoter, new brain specific PGC-1α isoforms and new variants associated with HD age at symptom onset .
In HD, independent analysis of three European HD cohorts suggested a protective effect of SNP 7665116 on age of onset in HD [16–18]. Further SNPs, for which an effect on AO was reported, include rs2970870, rs2970848 and rs6821591 [15, 19, 20]. These SNPs all map to non-coding portions of the gene.
As non-synonymous SNPs may be associated with functional changes of the protein, we have previously studied associations of rs8192678, the most frequent SNP in the coding region (minor allele frequency ~0.35) with HD AO, but found no associations in a smaller Italian population . As the current population is much larger and therefore increases the power to detect effects of less common SNPs, we studied the next common coding SNP, rs 3736265 (minor allele frequency ~0.05). Here we report, that this SNP is significantly associated with AO and that this association is male-specific.
We analyzed DNA samples and clinical data provided by the European HD network REGISTRY study, one of the largest cohorts available for this kind of study ; (http://www.euro-hd.net). In DNA from 1706 unrelated symptomatic HD patients (869 men/837 women), rs3736265 was typed using a TaqMan Genotyping Assay (C_25935762_20, Applied Biosystems). Correct genotyping was verified in 10% of the samples. Allele frequencies were estimated by gene counting. Agreement with Hardy-Weinberg equilibrium was ascertained using a χ2 goodness-of-fit test. Association of rs3736265 with motor AO was ascertained in a linear model. As we observed a skewed relationship between the expanded CAG repeat size and age of onset as expected, we used logarithmically transformed AO as the dependent variable and rs3736265 (coded as 0, 1 and 2) and normal and expanded CAG repeat sizes as well as their interactions as independent variables (Duyao M et al. 1993, Soyal SM et al. 2012). We observed a linear relationship between HD CAG repeat size and log-transformed AO. Furthermore, the residuals appeared to follow a normal distribution as judged by the Kolmogorov-Smirnov test. The access to the databases used in our work is publicly available. Participants gave informed written consent according to the International Conference on Harmonisation-Good Clinical Practice (ICH-GCP) guidelines (http://www.ich.org/LOB/media/MEDIA482.pdf). For participants, who lacked capacity to consent, study sites adhered to country-specific guidelines for obtaining consent. Minors assented with both parents consenting for them. Ethical approval was obtained from the local ethics committee for each study site contributing to REGISTRY .
Results and discussion
Effect of PPARGC1A rs3736265 on AO
HD CAG repeat size
non-HD CAG repeat size
HD CAG repeat size
non-HD CAG repeat size
HD CAG repeat size
non-HD CAG repeat size
The fact that rs3736265 is a low frequency SNP implies that small differences in analysis might affect the significance of the association. We therefore scrutinized the robustness of the association. First, three men were homozygous, possessing HD alleles with 54 and 37, 50 and 37 or 44 and 36 CAG repeats. Exclusion of these subjects from the analysis showed p-values of 0.0040 and 0.0036 for rs3736265 in the additive and dominant model, respectively. Second, differences in the relationship between CAG repeat length and AO were noted between the sixteen study centers. Adjustment for these differences did not affect the significance of the association of rs3736265 with logarithmically transformed AO (p-values 0.0047 and 0.0044 for additive and dominant model, respectively). Third, according to a previous study, the normal CAG repeat has no effect on AO . We therefore omitted the normal allele and its product with the HD allele from the analysis and found similar p-levels (0.0037 and 0.0034). Fourth, few samples at the end of the distribution with large repeats may disproportionally drive the results obtained. To address this possibility, we omitted the 21 subjects with HD-CAG size > 55 and observed similar p-values as in the entire male cohort (0.0044 and 0.0040). Additional exclusion of 12 patients with HD-CAG size <40 resulted in p-values of 0.0075 and 0.0070. We also addressed the issue of modifying effects in women with <55 CAG repeats and found no modifying effect of rs3736265 (P-values 0.882 and 0.843 for the additive and dominant models). Finally, the collection of subjects from different European regions and in different centers is likely to result in a reduced signal-to-noise ratio compared to smaller single-region or single-center studies. Hence, the association observed is robust and fulfills the expectations of linear models as linear relationships between variable and a normal distribution of residuals were observed.
The PPARGC1A genomic locus is of potential importance for understanding the mechanisms of PGC-1α action in HD pathogenesis. We thus examined this locus in more detail and demonstrate that a SNP in the coding region of PPARGC1A has a significant effect on age of onset in HD.
SNP rs3736265 is located in exon 9 of PPARGC1A, hence contained in many different tissue specific PGC-1α isoforms reported to date, i.e. PGC-1α1-4 and PGC-1αB1-5 [15, 23]. It results in a Thr612Met amino acid exchange. Hence, a polar, aliphatic amino acid (threonine) is substituted by a sulphur-containing hydrophobic, aliphatic amino acid (methionine). The amino acid stretch comprising the more common threonine is predicted to have a random coil or extended strand structure. The substitution by methionine is likely to result in an α-helical structure of five amino acids . Co-activation of PPARγ on the acyl-CoA-binding promoter did not differ between PGC-1α proteins harboring threonine or methionine at amino acid residue 612 . However, a recent study reported significant associations of rs3736265 with neovascularization in age-related macular degeneration (AMD) as well as interactions with the AMD-associated genes encoding complement factor B and VEGFA .
Surprisingly, the described SNP effect is clearly male-specific, while HD - unlike other neurodegenerative diseases - does not show a gender bias. This is reflected in our study population with a balanced male to female ratio (1.04). Some gender differences in HD are emerging: In the REGISTRY population women tend to have poorer motor and functional scores and a faster rate of disease progression than matched male controls . This however did not translate into differences in age at symptom onset. Mochel et al. found a lower body mass index (BMI) and lean body-mass despite a higher caloric intake in men but not in women . Similarly, a transgenic rat model of HD displayed a marked male-specific gender-bias of symptoms . This study includes a gene-expression profiling study on human HD striata showing only 42 genes to be differentially regulated between men vs. women with HD, compared to 343 genes differentially regulated between male vs. female controls. A possible explanation for the weakness of this signal is that the 1.5 fold-change cut-off in the analysis is too strict to allow for PGC-1α connected changes to become evident . Thus, a net-work based analysis may have been superior to a bioinformatics approach that looks at each transcript in isolation. Finally, sex-specific influences of genetic modifiers in HD have been reported for apolipoprotein E  and for the NR2A and NR2B receptor genes , but these results were not confirmed in other studies [32, 33]. An important difference between these studies and our present work is that our HD cohort is larger by an order of magnitude (n = 138 and n = 250 vs. n = 1706, respectively). Nevertheless, the association we report here awaits independent confirmation.
In summary our result further strengthens the accumulating evidence for a role of the PGC-1α-pathway in HD in particular and neurodegeneration in general and underlines that gender differences in neurodegenerative diseases might be due to differences at the molecular level .
We are indebted to the patients and their families for their participation in this project. We thank the clinicians who referred the patients to us. This study was supported by grants from the Fonds zur Förderung der wissenschaftlichen Forschung (FWF Project P19893-B05) to W.P., the European Huntington’s Disease Network (EHDN) REGISTRY to G.B.L., W.P. and P.W., and Paracelsus Medical University and the Verein für Medizinische Forschung Salzburg, Austria, to W.P. P.W. is an investigator of the Virtual Institute of the Deutsches Zentrum für Neurodegenerative Erkrankungen (DZNE).
- Arning L, Epplen JT: Genetic modifiers of Huntington’s disease: beyond CAG. Future Neurol. 2012, 7: 93-109. 10.2217/fnl.11.65.View Article
- Gusella JF, MacDonald ME: Huntington's disease: the case for genetic modifiers. Genome Med. 2009, 1: 80-10.1186/gm80.PubMedPubMed CentralView Article
- Lin J, Wu P-H, Tarr PT, Lindenberg KS, St-Pierre J, Zhang C-Y, Mootha VK, Jäger S, Vianna CR, Reznick RM, Cui L, Manieri M, Donovan MX, Wu Z, Cooper MP, Fan MC, Rohas LM, Zavacki AM, Cinti S, Shulman GI, Lowell BB, Krainc D, Spiegelman BM: Defects in Adaptive Energy Metabolism with CNS-Linked Hyperactivity in PGC-1α Null Mice. Cell. 2004, 119: 121-135. 10.1016/j.cell.2004.09.013.PubMedView Article
- Leone TC, Lehman JJ, Finck BN, Schaeffer PJ, Wende AR, Boudina S, Courtois M, Wozniak DF, Sambandam N, Bernal-Mizrachi C, Chen Z, Holloszy JO, Medeiros DM, Schmidt RE, Saffitz JE, Abel ED, Semenkovich CF, Kelly DP: PGC-1α Deficiency Causes Multi-System Energy Metabolic Derangements: Muscle Dysfunction Abnormal Weight Control and Hepatic Steatosis. Plos Biol. 2005, 3: e101-10.1371/journal.pbio.0030101.PubMedPubMed CentralView Article
- Weydt P, Pineda VV, Torrence AE, Libby RT, Satterfield TF, Lazarowski ER, Gilbert ML, Morton GJ, Bammler TK, Strand AD, Cui L, Beyer RP, Easley CN, Smith AC, Krainc D, Luquet S, Sweet IR, Schwartz MW, La Spada AR: Thermoregulatory and metabolic defects in Huntington's disease transgenic mice implicate PGC-1α in Huntington's disease neurodegeneration. Cell Metab. 2006, 4: 349-362. 10.1016/j.cmet.2006.10.004.PubMedView Article
- Cui L, Jeong H, Borovecki F, Parkhurst CN, Tanese N, Krainc D: Transcriptional Repression of PGC-1α by Mutant Huntingtin Leads to Mitochondrial Dysfunction and Neurodegeneration. Cell. 2006, 127: 59-69. 10.1016/j.cell.2006.09.015.PubMedView Article
- Szalardy L, Zadori D, Plangar I, Vecsei L, Weydt P, Ludolph A, Klivényi P, Kovacs G: Neuropathology of partial PGC-1α deficiency recapitulates features of mitochondrial encephalopathies but not of neurodegenerative diseases. Neurodegener Dis. 2012, 1-22.
- Handschin C: The biology of PGC-1α and its therapeutic potential. Trends Pharmacol Sci. 2009, 30: 322-329. 10.1016/j.tips.2009.03.006.PubMedView Article
- 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.PubMedView Article
- Esterbauer H, Oberkofler H, Krempler F, Patsch W: Human Peroxisome Proliferator Activated Receptor Gamma Coactivator 1 (PPARGC1) Gene: cDNA Sequence, Genomic Organization, Chromosomal Localization, and Tissue Expression. Genomics. 1999, 62: 98-102. 10.1006/geno.1999.5977.PubMedView Article
- Soyal S, Krempler F, Oberkofler H, Patsch W: PGC-1α: a potent transcriptional cofactor involved in the pathogenesis of type 2 diabetes. Diabetologia. 2006, 49: 1477-1488. 10.1007/s00125-006-0268-6.PubMedView Article
- Zheng B, Liao Z, Locascio JJ, Lesniak KA, Roderick SS, Watt ML, Eklund AC, Zhang-James Y, Kim PD, Hauser MA, Grunblatt E, Moran LB, Mandel SA, Riederer P, Miller RM, Federoff HJ, Wullner U, Papapetropoulos S, Youdim MB, Cantuti-Castelvetri I, Young AB, Vance JM, Davis RL, Hedreen JC, Adler CH, Beach TG, Graeber MB, Middleton FA, Rochet J-C, Scherzer CR: PGC-1α, A Potential Therapeutic Target for Early Intervention in Parkinson's Disease. Sci Transl Med. 2010, 2: 52ra73-PubMedPubMed CentralView Article
- Eschbach J, Schwalenstocker B, Soyal SM, Bayer H, Wiesner D, Akimoto C, Nilsson AC, Birve A, Meyer T, Dupuis L, Danzer KM, Andersen PM, Witting A, Ludolph AC, Patsch W, Weydt P: PGC-1α is a male-specific disease modifier of human and experimental amyotrophic lateral sclerosis. Hum Mol Genet. 2013, 22: 3477-3484. 10.1093/hmg/ddt202.PubMedView Article
- Róna-Vörös K, Weydt P: The role of PGC-1α in the pathogenesis of neurodegenerative disorders. Curr Drug Targets. 2010, 11: 1262-1269. 10.2174/1389450111007011262.PubMedView Article
- Soyal SM, Felder TK, Auer S, Hahne P, Oberkofler H, Witting A, Paulmichl M, Landwehrmeyer GB, Weydt P, Patsch W, For the European Huntington Disease Network: A greatly extended PPARGC1A genomic locus encodes several new brain-specific isoforms and influences Huntington disease age of onset. Hum Mol Genet. 2012, 21: 3461-3473. 10.1093/hmg/dds177.PubMedView Article
- Weydt P, Soyal S, Gellera C, DiDonato S, Weidinger C, Oberkofler H, Landwehrmeyer G, Patsch W: The gene coding for PGC-1α modifies age at onset in Huntington's Disease. Mol Neurodegener. 2009, 4: 3-10.1186/1750-1326-4-3.PubMedPubMed CentralView Article
- Taherzadeh-Fard E, Saft C, Andrich J, Wieczorek S, Arning L: PGC-1α as modifier of onset age in Huntington disease. Mol Neurodegener. 2009, 4: 10-10.1186/1750-1326-4-10.PubMedPubMed CentralView Article
- Che HVB, Metzger S, Portal E, Deyle C, Riess O, Nguyen HP: Localization of sequence variations in PGC-1α influence their modifying effect in Huntington disease. Mol Neurodegener. 2011, 6: 1-10.1186/1750-1326-6-1.PubMedPubMed CentralView Article
- Chen Y, Brew B: Characterization of the kynurenine pathway in nsc 34 cell line: implications for amyotrophic lateral sclerosis. J Neurochem.
- Taherzadeh-Fard E, Saft C, Akkad DA, Wieczorek S, Haghikia A, Chan A, Epplen JT, Arning L: PGC-1alpha downstream transcription factors NRF-1 and TFAM are genetic modifiers of Huntington disease. Mol Neurodegener. 2011, 6: 32-10.1186/1750-1326-6-32.PubMedPubMed CentralView Article
- Orth M, The European Huntington 's Disease Network: Observing Huntington's disease: the European Huntington‘s Disease Network’s REGISTRY. J Neurol Neurosurg Psychiatry. 2011, 82: 1409-1412. 10.1136/jnnp.2010.209668.PubMedView Article
- Lee J-M, Ramos EM, Lee J-H, Gillis T, Mysore JS, Hayden MR, Warby SC, Morrison P, Nance M, Ross CA, Margolis RL, Squitieri F, Orobello S, Di Donato S, Gomez-Tortosa E, Ayuso C, Suchowersky O, Trent RJA, Mccusker E, Novelletto A, Frontali M, Jones R, Ashizawa T, Frank S, Saint-Hilaire MH, Hersch SM, Rosas HD, Lucente D, Harrison MB, Zanko A, et al: CAG repeat expansion in Huntington disease determines age at onset in a fully dominant fashion. Neurology. 2012, 78: 690-695. 10.1212/WNL.0b013e318249f683.PubMedPubMed CentralView Article
- Ruas JL, White JP, Rao RR, Kleiner S, Brannan KT, Harrison BC, Greene NP, Wu J, Estall JL, Irving BA, Lanza IR, Rasbach KA, Okutsu M, Nair KS, Yan Z, Leinwand LA, Spiegelman BM: A PGC-1α Isoform Induced by Resistance Training Regulates Skeletal Muscle Hypertrophy. Cell. 2012, 151: 1319-1331. 10.1016/j.cell.2012.10.050.PubMedPubMed CentralView Article
- Combet C, Blanchet C, Geourjon C, Deléage G: NPS@: network protein sequence analysis. Trends Biochem Sci. 2000, 25: 147-150. 10.1016/S0968-0004(99)01540-6.PubMedView Article
- Nitz I, Ewert A, Klapper M, Döring F: Analysis of PGC-1α variants Gly482Ser and Thr612Met concerning their PPARγ2-coactivation function. Biochem Biophys Res Commun. 2007, 353: 481-486. 10.1016/j.bbrc.2006.12.042.PubMedView Article
- SanGiovanni JP, Chen J, Sapieha P, Aderman CM, Stahl A, Clemons TE, Chew EY, Smith LEH: DNA Sequence Variants in PPARGC1A, a Gene Encoding a Coactivator of the ω-3 LCPUFA Sensing PPAR-RXR Transcription Complex Are Associated with NV AMD and AMD-Associated Loci in Genes of Complement and VEGF Signaling Pathways. PLoS ONE. 2013, 8: e53155-10.1371/journal.pone.0053155.PubMedPubMed CentralView Article
- Zielonka D, Marinus J, Roos RAC, De Michele G, Di Donato S, Putter H, Marcinkowski J, Squitieri F, Bentivoglio AR, Landwehrmeyer GB: The influence of gender on phenotype and disease progression in patients with Huntington's disease. Parkinsonism Relat Disord. 2012, 1-6.
- Mochel F, Charles P, Seguin F, Barritault J, Coussieu C, Perin L, Le Bouc Y, Gervais C, Carcelain G, Vassault A: Early energy deficit in Huntington disease: identification of a plasma biomarker traceable during disease progression. PLoS ONE. 2007, 2: 647-10.1371/journal.pone.0000647.View Article
- Bode FJ, Stephan M, Suhling H, Pabst R, Straub RH, Raber KA, Bonin M, Nguyen H-P, Riess O, Bauer A, Sjoberg C, Petersen A, Horsten Von S: Sex differences in a transgenic rat model of Huntington's disease: decreased 17 -estradiol levels correlate with reduced numbers of DARPP32+ neurons in males. Hum Mol Genet. 2008, 17: 2595-2609. 10.1093/hmg/ddn159.PubMedView Article
- Kehoe P, Krawczak M, Harper PS, Owen MJ, Jones AL: Age of onset in Huntington disease: sex specific influence of apolipoprotein E genotype and normal CAG repeat length. J Med Genet. 1999, 36: 108-111.PubMedPubMed Central
- Arning L, Saft C, Wieczorek S, Andrich J, Kraus PH, Epplen JT: NR2A and NR2B receptor gene variations modify age at onset in Huntington disease in a sex-specific manner. Hum Genet. 2007, 122: 175-182. 10.1007/s00439-007-0393-4.PubMedView Article
- Saft C: Apolipoprotein E genotypes do not influence the age of onset in Huntington's disease. J Neurol Neurosurg Psychiatry. 2004, 75: 1692-1696. 10.1136/jnnp.2003.022756.PubMedPubMed CentralView Article
- Andresen JM, Gayan J, Cherny SS, Brocklebank D, Alkorta-Aranburu G, Addis EA, Cardon LR, Housman DE, Wexler NS, The US-Venezuela Collaborative Research Group: Replication of twelve association studies for Huntington's disease residual age of onset in large Venezuelan kindreds. J Med Genet. 2006, 44: 44-50. 10.1136/jmg.2006.045153.PubMedPubMed CentralView Article
- Jazin E, Cahill L: Sex differences in molecular neuroscience: from fruit flies to humans. Nat Rev Neurosci. 2010, 11: 9-17.PubMedView Article
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2377/14/1/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.