PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2010; 5(9): e12584.
Published online 2010 September 3. doi:  10.1371/journal.pone.0012584
PMCID: PMC2933242

A PPARα Promoter Variant Impairs ERR-Dependent Transactivation and Decreases Mortality after Acute Coronary Ischemia in Patients with Diabetes

Annarosa Leri, Editor

Abstract

Activation of peroxisome proliferator-activated receptor alpha (PPARα) occurs in animal models of diabetes (DM) and is implicated in pathological responses to myocardial ischemia. Using bioinformatics, we identified a single nucleotide polymorphism (SNP) in the PPARα gene promoter (PPARA −54,642 G>A; rs135561) that altered the consensus sequence for a nuclear receptor binding site. Electrophoretic mobility shift assays showed that the domain bound two known PPARA transcriptional activators, estrogen-related receptor (ERR)-α and -γ and that PPARA G bound with greater affinity than PPARA A (>2-fold; P<0.05). Likewise, promoter-reporter analyses showed enhanced transcriptional activity for PPARA G vs. PPARA A for both ERR-α and -γ (3.1 vs.1.9-fold; P<0.05). Since PPARα activation impairs post-ischemic cardiac function in experimental models of DM, we tested whether decreased PPARA transcription in PPARA A carriers favorably impacted outcome after acute coronary ischemia in 705 patients hospitalized with acute coronary syndromes (ACS; 552 Caucasian, 106 African American). PPARA A allele frequencies were similar to non-diseased subjects. However, PPARA genotype correlated with 5-year mortality in diabetic (22.2% AA vs. 18.8% AG vs. 39.5% GG; P = 0.008), but not non-diabetic (P = 0.96) subjects (genotype by diabetes interaction P = 0.008). In the diabetic ACS subjects, PPARA A carriers had strikingly reduced all-cause mortality compared to PPARA G homozygotes, (unadjusted HR 0.44, 95% CI 0.26–0.75; P = 0.003; adjusted HR 0.48, 95% CI 0.27–0.83; P = 0.009). Consistent with previous descriptions of PPARα in experimental models and human disease, we describe a novel PPARA promoter SNP that decreases transcriptional activation of PPARA and protects against mortality in diabetic patients after ACS.

Introduction

Peroxisome proliferator-activated receptor alpha (PPARα) is a transcription factor involved in the regulated expression of genes that control fatty acid uptake and metabolism. PPARα is highly expressed in the heart and other tissues that rely on fatty acid oxidation as their primary energy substrate [1]. Under normal conditions, regulation of myocardial PPARα expression and activity contributes to maintaining homeostatic balance between cellular fatty acid and glucose utilization via PPARα-mediated activation of target genes. This balance is perturbed by acute coronary ischemia [2][11], where myocardial PPARα expression can be a powerful determinant of functional recovery. For example, studies in which PPARα gene expression was experimentally manipulated in mouse models have demonstrated that forced myocardial PPARα overexpression decreases cardiac recovery after ischemia-reperfusion injury [12], whereas lowering PPARα expression protects against ischemic damage [12], [13]. Accordingly, genetic and epigenetic factors that influence PPARα gene (PPARA) expression have the potential to modify outcomes following acute coronary ischemia in patients.

Relevant to the human condition, we and others have suggested that the pathophysiological impact of PPARα may be greater in diabetic heart disease [7], [8], [14]. Insulin resistance and diabetes mellitus cause a shift of myocardial energy metabolism away from glucose and toward increased reliance on fatty acids that is driven in part by activation of PPARα [4], [15]. Animal models of diabetes and obesity, including STZ-induced diabetic rats, db/db mice, and ob/ob mice, show increased cardiac PPARα activation [15][17]. Given that cardiac PPARα is activated in DM, such a cardiac disease-modifier role might have special importance in diabetic patients. Indeed, a previously described single nucleotide polymorphism (SNP) within intron 1 of PPARA has been associated with the age of onset, and progression to insulin monotherapy, in patients with type 2 DM [18], although the biological mechanism for this effect has not been elucidated. Here, we identify and characterize a pair of linked SNPs in the 5′-flanking promoter region of human PPARA that modify a consensus nuclear hormone receptor response element to decrease PPARA transcription. We found that this genotype confers protection against early post-ischemic mortality that is specific to diabetic patients. Since these SNPs are also in linkage disequilibrium (LD) with the previously described intronic polymorphism, these results also suggest a mechanism for the clinical effects of that SNP.

Methods

Ethics Statement

The study was approved by the University of Missouri Kansas City Adult Health Sciences Institutional Review Board, the Saint Luke's Hospital of Kansas City Institutional Review Board and the Washington University Human Research Protection Office Institutional Review Board. Written informed consent was provided by each participant. A separate consent form for the acquisition of blood for genetic analysis was signed by each participant.

Subjects

Patients were prospectively enrolled into the INFORM ACS registry at two Kansas City hospitals as previously described [19], [20]. Between 3/1/01 and 10/31/02, 1199 patients met the criteria for ACS using standard, accepted definitions of MI (n = 680) and unstable angina (n = 519) [21], [22]. MI patients were defined by an elevated troponin blood test [21]. Patients with DM were defined by having been given a diagnosis of DM by the referring physician and/or being treated with oral hypoglycemic agents or insulin. Three physicians reviewed the charts of all patients with diagnostic uncertainty and attained consensus on the final diagnosis.

Although there were no differences in gender (93.2% of men vs. 92.1% of women), Caucasians (91.5% vs. 98.3%, p<0.001) and older patients (mean age for those consenting = 61±13 vs. 65±13 for those not consenting, p = 0.004) were less likely to consent to DNA testing. 726 patients were enrolled in the genetic portion of this registry, constituting the cohort for the present analysis.

Mortality Assessment

The Social Security Administration Death Master File was queried to determine patients' vital status as of 03/12/2008 (http://www.ntis.gov/products/ssa-dmf.asp).

Bioinformatic Analysis

As an initial approach to discovering functionally significant polymorphisms of the PPARA promoters, we used the web-based program PromoLign (http://polly.wustl.edu/promolign/main.html) [23], to identify nucleotide variants within 10 kb of the transcription start site. This program identifies promoter polymorphisms that are within human-mouse homologous blocks and/or occur within, and alter, the consensus sequences for putative trans-activating factor binding sites (output shown in Supplemental Figure S1). We considered that four essential criteria were necessary to assign functional importance to PPARA promoter polymorphisms: 1. The bioinformatic SNP had to be validated in a multi-ethnic human cohort; 2. The putative nuclear receptor binding domain in which it is located had to bind a transcription factor; 3.The SNP had to alter binding of that transcription factor; and 4. The SNP had to change gene transcriptional activity.

Targeted Resequencing

Targeted resequencing of the putative nuclear receptor domain was performed using pyrosequencing in 380 DNA samples from unrelated, random, healthy blood donors, 95 each from Caucasians, African-Americans, Han Chinese, and Mexicans (Coriell Institute for Medical Research, Camden NJ).

Genotyping

DNA was isolated and extracted using the Puregene genomic DNA purification kit (Gentra, Minneapolis, MN). The DNA segments containing the region of interest were amplified with the polymerase chain reaction (PCR). PCR primers were designed using Primer3 online software (http://fokker.wi.mit.edu/cgi-bin/primer3/primer3_www.cgi) [24], and pyrosequencing primers were designed using the Pyrosequencing SNP Primer Design Version 1.01 software (http://www.pyrosequencing.com). Before use, PCR primer sequences were screened across the human genome using the NCBI Blast program to ensure their specificity for the gene of interest. PCR was carried out using Amplitaq Gold PCR master mix (ABI, Foster City, CA), 1 pmole of each primer (IDT, Coralville, IA), and 1ng DNA. The PCR primers and conditions are listed in Supplemental Table S1. Pyrosequencing was performed using the PSQ HS 96A system with MA v2.0 software as previously described [25]. Data were automatically transferred from the PSQ HS 96A to a Microsoft Access database for permanent storage and merging with the clinical datasets through SAS v9.1. Pair wise linkage (D′) and haplotype analysis was carried out using the Polymorphism and Haplotype Analysis Suite [26].

In vitro Studies

Electrophoretic mobility shift assays used standard methodologies that have been previously described [27], [28]. Complementary oligonucleotides corresponding to the PPARA promoter region encompassing the SNP of interest were annealed to generate double-stranded fragments for radiolabeling or cloning. Oligonucleotide sequences were as follows: PPARA −54,642 G forward: gatcTAGCTCCTGCAGGTTCTCAAGGTTGTAGCCCGCACCCTGCT

PPARA −54,642 A forward: gatcTAGCTCCTGCAGGTTCTCAAGGTTGTAGTCCACACCCTGCT. Note that these sequences differed at two sites, the −54,642 SNP site and at a site 3 nucleotides upstream of this SNP site, due to the upstream site being in complete LD with the −54,642 SNP site; see results. The positive control element consisted of the previously characterized estrogen-related receptor (ERR) responsive element [28]. Probes were synthesized by Klenow fill-in reaction with [γ-32P]dCTP using the double-stranded fragments and used at 15,000 cpm per reaction. Binding reactions were performed as described [28] using recombinant ERRα and ERRγ generated from TNT Quick Coupled T7 reticulocyte lysate (Promega). In competition assays cold competitor was added at the indicated concentrations simultaneously with probe and protein. Gels were imaged on a Storm phosphorimager and band intensities were quantified using ImageQuant software (Molecular Dynamics).

To clone heterologous PPARA variant promoter-reporter constructs, (−54,642A)2c.TK.Luc and (−54,642G)2c.TK.Luc, double-stranded fragments corresponding to the probes used in EMSA were 5′ phosphorylated then ligated into the BamHI site of a luciferase expression vector immediately upstream of the thymidine kinase minimal promoter of the pGL2.TK.Luc reporter plasmid. Clones were screened by PCR and those carrying 2 element copies were verified by sequence analysis. Transient transfection studies were performed as previously described [28], [29]. PPARA promoter-reporter constructs were cotransfected with empty expression vector or vector expressing ERRα or ERRγ. Transient transfections in CV1 cells were performed using Lipofectamine 2000 (Invitrogen) with 2.7 µg/ml reporter, 0.3 µg/ml each expression vector, and 0.3 µg/ml pRL-CMV to control for transfection efficiency. C2C12 myoblasts were transfected with 4 µg/ml reporter, 0.5 µg/ml each expression vector, and 0.5 µg/ml pRL-CMV using the calcium phosphate precipitation method. Firefly luciferase activity normalized to that of Renilla luciferase was measured 48h post-transfection on a Clarity luminescence microplate reader (BioTek) with triplicate determinations in four trials.

Statistical Analyses

Baseline patient characteristics were summarized and compared stratified by genotype. Hardy-Weinberg equilibrium was assessed by chi-square in Caucasians and African Americans separately. Continuous variables were reported as mean ± standard deviation and compared using t-tests. Lipid values had skewed distributions, were log-transformed prior to analysis, and were summarized by median and interquartile range. Categorical variables were summarized by frequency and percent and compared using chi-square or Fisher's exact tests, when the expected cell size was <5.

The primary outcome was time to all-cause mortality through 60 months. Event rates by PPARA −54,642 genotype were calculated using Kaplan-Meier analysis and compared using log-rank tests. The relative hazard associated with each genotype was estimated using Cox proportional hazard models, both in crude analysis and adjusting for demographics (age, gender, race) and characteristics that were significantly different across genotypes (systolic BP and treatment strategy).

Proportional hazard assumptions were verified using Schoenfeld residuals. P-values<0.05 were considered statistically significant. Analyses were performed using SAS version 9.1 (SAS Institute, Inc., Cary, NC) and R version 2.2.0.(http://www.R-project.org).

Luciferase activity and binding level values were summarized by mean and standard error and compared using t-tests. P values<0.05 were considered statistically significant.

Results

Bioinformatics analysis of PPARA promoters A–D

The PPARα gene (PPARA) has a complex 5′ structure, consisting of four distinct promoter regions designated promoters A–D [30]. 13 nucleotide variants were identified within 10 kb of the transcription start site as described in Methods (Table 1 and Supplemental Figure S1). Of the 13 PPARA promoter polymorphisms in the database within human-mouse homologous blocks, only 2 single nucleotide polymorphisms (SNPs) located 3 nucleotides apart from one another, occurred within, and altered the consensus sequences for, putative trans-activating factor binding sites. Of particular interest because of the previously established important effects of nuclear hormone receptors on PPARA expression [28], [31], was a G to A substitution located 6,862 nucleotides upstream of the major PPARA transcription start site of promoter A that disrupts a nuclear receptor half-site, designated PPARA −54,642 G>A (based on its position relative to the translation initiation site; rs135561). This domain is highly conserved across vertebrate species, including mouse, rat, rabbit, dog, chicken, elephant and human (Supplemental Figure S2) and was therefore selected for further genotyping and analysis.

Table 1
13 nucleotide variants identified within 10 kb of the PPARA transcription start site located within human-mouse homologous blocks and/or within consensus sequences for putative trans-activating factor binding sites.

The PPARA −54,642 A variant is in complete LD with a PPARA −54,645 T variant and is the major allele in African Americans

Targeted resequencing of the putative nuclear receptor domain was performed using pyrosequencing in 380 DNA samples from random, healthy blood donors, 95 each from Caucasians, African-Americans, Han Chinese, and Mexicans (Coriell Institute for Medical Research, Camden NJ). The PPARA −54,642 G>A substitution was confirmed, as was a C to T substitution three nucleotides 5′ (PPARA −54,645 C>T; rs135562) that was in complete LD (D′ = 1, r2 = 93, P<0.001) in all ethnic cohorts. Frequencies of the variant PPARA −54,642 A/PPARA −54,645 T allele were 0.27 in the Caucasian cohort, 0.68 in the African-American cohort, 0.15 in the Chinese cohort, and 0.12 in the Mexican cohort, and were each consistent with predictions of Hardy-Weinberg equilibrium. Thus, the PPARA −54,642 A variant (and its linked upstream partner) is relatively common in Caucasians, and is the most common allele in African Americans, but is relatively uncommon among Mexicans and Asians (≤15%).

The PPARA −54,642 G>A SNP alters ERRα and ERRγ binding and PPARA transcriptional activation

Having confirmed the bioinformatics in a multi-racial human cohort, the next step was to determine if there were any effects of the PPARA promoter SNP on transcription factor binding and transcriptional activity. We considered that the most likely factors to bind to the PPARA nuclear receptor half-site were estrogen related receptor (ERR)α and ERRγ, based on the characteristic nucleotide sequence [32], [33], on known similarities in tissue expression between ERR and PPARα [34][36], and on ERR regulatory activities for genes downstream of PPARα [28], [35], [36]. To test ERR binding, a 49 bp double-stranded (ds) oligonucleotide was generated corresponding to the ‘wild-type’ form of PPARA −54,642 (G) and −54,645 (C). Electrophoretic mobility shift assays (EMSA) were performed to compare binding of recombinant ERRα and ERRγ proteins between 32P-labeled PPARA and authentic, previously described ERR binding sequences [28]. Both ERRα and ERRγ bound the wild-type PPARA −54,642 G sequence with activities equal to, or greater than, the control ERR responsive element (Probe G vs. +control; Figure 1A).

Figure 1
PPARA −54,642 A variant binds ERRα and ERRγ with less affinity.

We next examined relative ERR binding in EMSA studies performed with the same wild-type PPARA −54,642 G probe described above and a ‘variant’ PPARA −54,642 A probe, identical in sequence except for the nucleotides encoding the two linked SNPs (A at nucleotide −54,642 and T at nucleotide −54,645, respectively). The apparent binding of both ERRα and ERRγ for the variant PPARA −54,642 A probe was qualitatively less than for the wild-type G probe (Probe A vs. Probe G; Figure 1A). The relative affinity of ERRα for the wild-type and variant sequences was assessed in competition binding studies using increasing amounts of unlabeled G or A dsDNA probe to displace wild-type G probe in the EMSA binding reactions. The variant A probe competed less effectively for binding to the 32P-labelled wild-type sequence (IC50 of A probe 2.3-fold greater than G probe, P<0.05), demonstrating a lower affinity for ERRα (Figure 1B and 1C) and ERRγ (data not shown). Thus, the G>A substitution decreases PPARA promoter affinity for the transcription factors ERRα and ERRγ.

To determine whether the decrease in ERR binding afforded by the PPARA −54,642 G>A polymorphism altered transcriptional activity, we performed reporter assays. The reporter constructs contained two copies of either the wild-type G or the variant A 49 bp promoter element inserted upstream of a minimal thymidine kinase (TK) promoter-luciferase reporter [29]. When co-transfected with an ERRα expression plasmid into C2C12 myoblasts, which are a useful in vitro model of skeletal and cardiac muscle [37], the variant A promoter-reporter was ~25% less responsive to ERRα-mediated activation than the wild-type G promoter-reporter (p<0.05; n = 4; Figure 2A). Similar decreased activation of variant PPARA −54,642 A was observed with ERRγ co-transfection (data not shown).

Figure 2
PPARA −54,642 A variant is less responsive to ERR-mediated co-activation.

To control for the possibility that endogenous ERRα in C2C12 myoblasts was masking PPARA/TK A promoter-reporter activation by transfected ERRα, we repeated the experiments in CV1 cells that are null for endogenous ERRα and ERRγ [28]. In the CV1 system the wild-type G PPARA promoter-reporter had ~60% greater responsiveness to ERRα-mediated activation than the A variant promoter-reporter (3.1-fold vs.1.9-fold; P<0.05; n = 4, Figure 2B). Again, the results were similar with co-transfection of ERRγ (data not shown). These studies demonstrate that the variant PPARA promoter sequence exhibits diminished ERR binding and decreased ERR-mediated transcriptional activation, compared to the wild-type PPARA −54,642 G allele.

In diabetic patients, PPARA -54,642 A allele carriers have deceased mortality after acute coronary ischemia

The above studies demonstrate clear biological effects of the PPARA −54,642 G>A promoter polymorphism: The more common (in Caucasians) PPARA G promoter is more responsive to ERR, while the variant A promoter is less effective at binding ERR and is therefore resistant to ERR-mediated transcription. Since we have previously shown that cardiac PPARα is activated in experimental DM [15], and that such PPARα activation is deleterious during experimental myocardial ischemia [38][49], we hypothesized that diminished PPARα gene activity conferred by the −54,642 A variant might have the opposite effect, i.e. be protective, in human subjects with acute coronary ischemic syndromes (ACS). We tested this idea in the 726 subject genetic sub-study of the INFORM trial of outcomes after ACS [19]. Baseline clinical and demographic data for the 705 patients (96%) for whom genotypes were obtained are listed in Table 2, and for patients with DM stratified by PPARA promoter genotype in Table 3. The PPARA −54,642 A allele frequency in ACS patients was 0.30 in Caucasians and 0.60 in African-Americans, which is not significantly different from the frequencies in the healthy volunteers used for SNP validation (P = 0.53 and  = 0.26, respectively; see above), and suggests that the PPARA SNP is not an independent risk factor for cardiac ischemic disease. As in the non-diseased cohort, the adjacent PPARA −54,645 C>T SNP (rs135562) was found to be in complete LD (D′ = 1, r2 = 93, P<0.001) in both Caucasian and African-American subgroups, and did not deviate significantly from Hardy-Weinberg equilibrium.

Table 2
Clinical Characteristics of the INFORM ACS Cohort.
Table 3
Clinical Characteristics of the Diabetic INFORM ACS Cohort according to PPARA −54,642 genotype.

To determine if the PPARA promoter variant modified outcome in patients with acute coronary ischemic syndromes we examined the relationship between PPARA −54,642 (rs135561) genotype and five year mortality in the INFORM study. When DM status was not a co-variate, there was no apparent association. However, when the cohort was stratified by DM status, PPARA −54,642 genotype was associated with increased 5-year mortality in diabetic (p = 0.008) but not in non-diabetic (p = 0.96) ACS patients (Figure 3; genotype by DM interaction p = 0.008). In non-diabetic ACS patients, 5-year mortality was 13.5% for wild-type PPARA GG homozygous subjects, and 12.7% for variant PPARA A allele carriers (Figure 3B). Consistent with the known adverse effects on ACS outcome of having diabetes [50][52], homozygous −54,642 GG (wild-type) PPARA ACS patients had a 5-year mortality rate of 39.5% (Figure 3A). By contrast, PPARA variant A allele carriers (AG or AA) had a 5-year mortality rate following admission for ACS of 19.9% (Figure 3A and Table 4; HR 0.44, 95% CI 0.26–0.75; P = 0.003), demonstrating a striking protective effect of the PPARA promoter polymorphism. When adjusted for age, race, and gender, genotype remained an independent predictor of 5-year mortality, with PPARA −54,642 A allele carriers having less than half the relative risk of death within 5 years of presenting with ACS, compared to wild-type (homozygous PPARA −54,642 GG genotype) diabetics (HR 0.48, 95% CI 0.27–0.83; p = 0.009). This relationship was not changed by further adjusting for any other variable that differed significantly between genotype groups (HR 0.48, 95% CI 0.26–0.89; p = 0.019). Notwithstanding the differences in PPARA −54,642 allele frequency between Caucasians and African-Americans, in Caucasian DM subjects alone (the largest racial group) PPARA −54,642 A carriers were significantly protected against 5-year mortality compared to wild-type GG allele patients (unadjusted HR 0.47, 95% CI 0.24–0.90; P = 0.023; adjusted HR 0.46, 95% CI 0.24–0.89; P = 0.021; Table 4).

Figure 3
In diabetic patients, PPARA −54,642 A allele carriers have deceased mortality after acute coronary ischemia.
Table 4
Unadjusted and adjusted 5-year mortality for PPARA −54,642 A allele carriers.

As mentioned above, an A>C SNP (rs135539) located in intron 1 (promoter B) of PPARA, 35,014 base pairs upstream of the translation start site, has been reported to influence the age of onset and progression to insulin monotherapy in patients with type 2 DM [18]. Although the authors provided no mechanism to explain these findings, they speculated that C allele carriers had reduced PPARα expression [18]. To evaluate whether PPARA −35,014 genotype impacted our findings, we genotyped the 726 subject genetic substudy of the INFORM trial for this polymorphism and found that the two SNPs (PPARA −54,642 and PPARA −35,014) were in moderate LD (D′ = 0.58, r2 = 17, P<0.001 Caucasians; D′ = 0.78, r2 = 39, P<0.001 African Americans). In addition, we found a borderline significant association with mortality in diabetic subjects within our cohort (unadjusted P = 0.051; AA vs. AC vs CC). Although this association only achieved borderline significance, there was a highly significant genotype by diabetes interaction (p = 0.009).

Discussion

The current results define a new genetic disease modifier of outcome after cardiac ischemia that is specific for diabetes and suggests a biologic mechanism for the previously described PPARA intronic SNP. We used bioinformatics to identify candidate SNPs with the potential to modify transcription factor binding and high homology to mouse promoter regions in the four PPARA promoters. One candidate, the PPARA −54,642 (rs135561) G>A polymorphism, interrupted a putative nuclear factor biding domain and altered binding of two critical regulatory factors for PPARA expression, ERRα and ERRγ. In vitro promoter-reporter assays revealed that this effect changed transcriptional activity of the PPARA promoter, and our clinical association studies revealed a significant impact on 5-year mortality in patients with DM and acute coronary ischemia.

We took advantage of an experiment of nature, a common pair of SNPs that alter PPARA promoter activity, to demonstrate that PPARA expression level helps determine outcome after myocardial infarction and unstable angina in human subjects with DM. This result is consistent with previous findings in genetic mouse models that PPARα expression level inversely correlates with post-ischemic cardiac function [12] and extends the previous report of an association between PPARA −35,014 (rs135539) and onset and progression of DM [18]. Taken together, the human, mouse, and cell culture studies represent a comprehensive evaluation of PPARα effects in ischemic hearts.

A growing body of evidence describes associations between the PPARα gene and cardiovascular disease, risk or outcomes in subjects with type 2 DM [18], [53][58]. Of particular relevance to the current findings is the A>C SNP (rs135539) located in intron 1 of PPARA reported to influence the age of onset and progression to insulin monotherapy in patients with type 2 DM [18]. No mechanism has been forthcoming to explain these findings but the authors speculated that C allele carriers had reduced PPARα expression [18]. The current results suggest an alternate explanation for the previously described associations with the PPARA intronic SNP and provide evidence to support the authors' hypothesis that C allele carriers have reduced PPARα expression.

Increased risk of mortality after ACS in patients with DM has been recognized for many years [38][41], [43], [45], but the specific reasons for this adverse outcome are unknown. It is notable that the recent therapeutic advances that have improved general outcomes in ACS, have not impacted the disparity in outcomes in diabetic individuals [42], [44], [46][49]. Data from clinical trials also suggest that PPAR activation with pharmacological agonists increases the risk of myocardial infarction and adverse cardiovascular events in patients with DM [59], [60]. Since our results identify a gene-disease interaction that alters mortality specifically in the diabetic population after cardiac ischemia, they suggest that a therapeutic approach targeting PPAR-responsive genes and metabolic pathways might favorably affect outcomes.

Our studies suggest that modified ERR binding is a mechanism for altered PPARA expression in the polymorphic promoter, but there is scant available data on factors that affect ERRα and ERRγ mRNA or protein expression levels in vivo. We previously found that ERRα and ERRγ activate PPARA expression via direct binding to PPARA promoter B (−49,076 to −49,064 upstream of the translation start site) [28] and observed that myocardial ERRα transcript expression is developmentally regulated in parallel with PPARα [61]. Interestingly, myocardial ERRα transcript expression is also up-regulated in animals fed high fat diets and in animal models of insulin-deficient DM (Huss, unpublished observation), providing further circumstantial evidence for a link between PPARα, ERR, and DM.

Our findings should be considered in the context of some potential limitations. Although genotype remained an independent predictor of 5-year mortality in both unadjusted and adjusted analyses, as well as in subgroup analyses performed in our largest racial group (Caucasians), replication in other studies is still warranted. Furthermore, our findings cannot be extrapolated to racial groups other than those represented in INFORM (Caucasians and African Americans) and should, therefore, be examined independently in other racial groups.

In conclusion, we have identified a PPARA promoter polymorphism that effects binding and activation of PPARA co-activators and is strongly associated with 5-year mortality in diabetic ACS patients. These findings suggest a genetic mechanism for the unfavorable cardiovascular outcomes in diabetic patients after ACS.

Supporting Information

Table S1

PCR and Pyrosequencing primers and conditions.

(0.04 MB DOC)

Figure S1

Output of the web-based program PromoLign (http://polly.wustl.edu/promolign/main.html) [23], showing the 13 nucleotide variants identified within 10 kb of the transcription start site of PPARA.

(2.85 MB TIF)

Figure S2

Sequence comparison of PPARA −54,642 SNP and PPARA −54,645 SNP sites showing conservation of sequence between species (dbSNP build 126; http://genome.ucsc.edu/).

(2.04 MB TIF)

Acknowledgments

We would like to thank Derek Van Booven for invaluable informatics assistance.

Footnotes

Competing Interests: Patent pending. Daniel P. Kelly is on the scientific advisory board for Lilly and Johnson and Johnson. Howard L. McLeod is a consultant to the Food and Drug Administration. This does not alter the authors' adherence to all the PLoS ONE policies on sharing data and materials.

Funding: This work has been funded by the National Institutes of Health (NIH) Specialized Center for Clinically-Oriented Research (SCCOR) in Cardiac Dysfunction and Disease # 1 P50 HL077113 and NIH Pharmacogenetics Research Network grant U01 GM63340. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. Braissant O, Foufelle F, Scotto C, Dauca M, Wahli W. Differential expression of peroxisome proliferator-activated receptors (PPARs): tissue distribution of PPAR-alpha, -beta, and -gamma in the adult rat. Endocrinology. 1996;137:354–366. [PubMed]
2. Barger PM, Brandt JM, Leone TC, Weinheimer CJ, Kelly DP. Deactivation of peroxisome proliferator-activated receptor-alpha during cardiac hypertrophic growth. J Clin Invest. 2000;105:1723–1730. [PMC free article] [PubMed]
3. Djouadi F, Weinheimer CJ, Saffitz JE, Pitchford C, Bastin J, et al. A gender-related defect in lipid metabolism and glucose homeostasis in peroxisome proliferator- activated receptor alpha- deficient mice. J Clin Invest. 1998;102:1083–1091. [PMC free article] [PubMed]
4. Finck BN, Han X, Courtois M, Aimond F, Nerbonne JM, et al. A critical role for PPARalpha-mediated lipotoxicity in the pathogenesis of diabetic cardiomyopathy: modulation by dietary fat content. Proc Natl Acad Sci U S A. 2003;100:1226–1231. [PubMed]
5. Guerre-Millo M, Rouault C, Poulain P, Andre J, Poitout V, et al. PPAR-alpha-null mice are protected from high-fat diet-induced insulin resistance. Diabetes. 2001;50:2809–2814. [PubMed]
6. Huss JM, Levy FH, Kelly DP. Hypoxia inhibits the peroxisome proliferator-activated receptor alpha/retinoid X receptor gene regulatory pathway in cardiac myocytes: a mechanism for O2-dependent modulation of mitochondrial fatty acid oxidation. J Biol Chem. 2001;276:27605–27612. [PubMed]
7. Huss JM, Kelly DP. Nuclear receptor signaling and cardiac energetics. Circ Res. 2004;95:568–578. [PubMed]
8. Huss JM, Kelly DP. Mitochondrial energy metabolism in heart failure: a question of balance. J Clin Invest. 2005;115:547–555. [PMC free article] [PubMed]
9. Lemberger T, Desvergne B, Wahli W. Peroxisome proliferator-activated receptors: a nuclear receptor signaling pathway in lipid physiology. Annu Rev Cell Dev Biol. 1996;12:335–363. [PubMed]
10. Lemberger T, Saladin R, Vazquez M, Assimacopoulos F, Staels B, et al. Expression of the peroxisome proliferator-activated receptor alpha gene is stimulated by stress and follows a diurnal rhythm. J Biol Chem. 1996;271:1764–1769. [PubMed]
11. Leone TC, Weinheimer CJ, Kelly DP. A critical role for the peroxisome proliferator-activated receptor alpha (PPARalpha) in the cellular fasting response: the PPARalpha-null mouse as a model of fatty acid oxidation disorders. Proc Natl Acad Sci U S A. 1999;96:7473–7478. [PubMed]
12. Sambandam N, Morabito D, Wagg C, Finck BN, Kelly DP, et al. Chronic activation of PPARalpha is detrimental to cardiac recovery after ischemia. Am J Physiol Heart Circ Physiol. 2006;290:H87–H95. [PubMed]
13. Panagia M, Gibbons GF, Radda GK, Clarke K. PPAR-alpha activation required for decreased glucose uptake and increased susceptibility to injury during ischemia. Am J Physiol Heart Circ Physiol. 2005;288:H2677–H2683. [PubMed]
14. Evans RM, Barish GD, Wang YX. PPARs and the complex journey to obesity. Nat Med. 2004;10:355–361. [PubMed]
15. Finck BN, Lehman JJ, Leone TC, Welch MJ, Bennett MJ, et al. The cardiac phenotype induced by PPARalpha overexpression mimics that caused by diabetes mellitus. J Clin Invest. 2002;109:121–130. [PMC free article] [PubMed]
16. Buchanan J, Mazumder PK, Hu P, Chakrabarti G, Roberts MW, et al. Reduced cardiac efficiency and altered substrate metabolism precedes the onset of hyperglycemia and contractile dysfunction in two mouse models of insulin resistance and obesity. Endocrinology. 2005;146:5341–5349. [PubMed]
17. Sharma S, Adrogue JV, Golfman L, Uray I, Lemm J, et al. Intramyocardial lipid accumulation in the failing human heart resembles the lipotoxic rat heart. FASEB J. 2004;18:1692–1700. [PubMed]
18. Flavell DM, Ireland H, Stephens JW, Hawe E, Acharya J, et al. Peroxisome proliferator-activated receptor alpha gene variation influences age of onset and progression of type 2 diabetes. Diabetes. 2005;54:582–586. [PubMed]
19. Lanfear DE, Jones PG, Marsh S, Cresci S, McLeod HL, et al. Beta2-adrenergic receptor genotype and survival among patients receiving beta-blocker therapy after an acute coronary syndrome. JAMA. 2005;294:1526–1533. [PubMed]
20. Spertus J, Safley D, Garg M, Jones P, Peterson ED. The influence of race on health status outcomes one year after an acute coronary syndrome. J Am Coll Cardiol. 2005;46:1838–1844. [PubMed]
21. Alpert JS, Thygesen K, Antman E, Bassand JP. Myocardial infarction redefined–a consensus document of The Joint European Society of Cardiology/American College of Cardiology Committee for the redefinition of myocardial infarction. J Am Coll Cardiol. 2000;36:959–969. [PubMed]
22. Braunwald E. Unstable angina. A classification. Circulation. 1989;80:410–414. [PubMed]
23. Zhao T, Chang LW, McLeod HL, Stormo GD. PromoLign: a database for upstream region analysis and SNPs. Hum Mutat. 2004;23:534–539. [PubMed]
24. Rozen S, Skaletsky H. Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000;132:365–386. [PubMed]
25. Marsh S, King CR, Garsa AA, McLeod HL. Pyrosequencing of clinically relevant polymorphisms. Methods Mol Biol. 2005;311:97–114. [PubMed]
26. Schaid DJ, Rowland CM, Tines DE, Jacobson RM, Poland GA. Score tests for association between traits and haplotypes when linkage phase is ambiguous. Am J Hum Genet. 2002;70:425–434. [PubMed]
27. Cresci S, Clabby ML, Kelly DP. Evidence for a novel cardiac-enriched retinoid X receptor partner. J Biol Chem. 1999;274:25668–25674. [PubMed]
28. Huss JM, Torra IP, Staels B, Giguere V, Kelly DP. Estrogen-related receptor alpha directs peroxisome proliferator-activated receptor alpha signaling in the transcriptional control of energy metabolism in cardiac and skeletal muscle. Mol Cell Biol. 2004;24:9079–9091. [PMC free article] [PubMed]
29. Gulick T, Cresci S, Caira T, Moore DD, Kelly DP. The peroxisome proliferator-activated receptor regulates mitochondrial fatty acid oxidative enzyme gene expression. Proc Natl Acad Sci U S A. 1994;91:11012–11016. [PubMed]
30. Chew CH, Samian MR, Najimudin N, Tengku-Muhammad TS. Molecular characterisation of six alternatively spliced variants and a novel promoter in human peroxisome proliferator-activated receptor alpha. Biochem Biophys Res Commun. 2003;305:235–243. [PubMed]
31. Pineda TI, Jamshidi Y, Flavell DM, Fruchart JC, Staels B. Characterization of the human PPARalpha promoter: identification of a functional nuclear receptor response element. Mol Endocrinol. 2002;16:1013–1028. [PubMed]
32. Wende AR, Huss JM, Schaeffer PJ, Giguere V, Kelly DP. PGC-1alpha coactivates PDK4 gene expression via the orphan nuclear receptor ERRalpha: a mechanism for transcriptional control of muscle glucose metabolism. Mol Cell Biol. 2005;25:10684–10694. [PMC free article] [PubMed]
33. Mootha VK, Handschin C, Arlow D, Xie X, St PJ, et al. Erralpha and Gabpa/b specify PGC-1alpha-dependent oxidative phosphorylation gene expression that is altered in diabetic muscle. Proc Natl Acad Sci U S A. 2004;101:6570–6575. [PubMed]
34. Heard DJ, Norby PL, Holloway J, Vissing H. Human ERRgamma, a third member of the estrogen receptor-related receptor (ERR) subfamily of orphan nuclear receptors: tissue-specific isoforms are expressed during development and in the adult. Mol Endocrinol. 2000;14:382–392. [PubMed]
35. Vega RB, Kelly DP. A role for estrogen-related receptor alpha in the control of mitochondrial fatty acid beta-oxidation during brown adipocyte differentiation. J Biol Chem. 1997;272:31693–31699. [PubMed]
36. Sladek R, Bader JA, Giguere V. The orphan nuclear receptor estrogen-related receptor alpha is a transcriptional regulator of the human medium-chain acyl coenzyme A dehydrogenase gene. Mol Cell Biol. 1997;17:5400–5409. [PMC free article] [PubMed]
37. Yaffe D, Saxel O. Serial passaging and differentiation of myogenic cells isolated from dystrophic mouse muscle. Nature. 1977;270:725–727. [PubMed]
38. Behar S, Boyko V, Reicher-Reiss H, Goldbourt U. Ten-year survival after acute myocardial infarction: comparison of patients with and without diabetes. SPRINT Study Group. Secondary Prevention Reinfarction Israeli Nifedipine Trial. Am Heart J. 1997;133:290–296. [PubMed]
39. Haffner SM, Lehto S, Ronnemaa T, Pyorala K, Laakso M. Mortality from coronary heart disease in subjects with type 2 diabetes and in nondiabetic subjects with and without prior myocardial infarction. N Engl J Med. 1998;339:229–234. [PubMed]
40. Malmberg K, Ryden L. Myocardial infarction in patients with diabetes mellitus. Eur Heart J. 1988;9:259–264. [PubMed]
41. Miettinen H, Lehto S, Salomaa V, Mahonen M, Niemela M, et al. Impact of diabetes on mortality after the first myocardial infarction. The FINMONICA Myocardial Infarction Register Study Group. Diabetes Care. 1998;21:69–75. [PubMed]
42. Gu K, Cowie CC, Harris MI. Diabetes and decline in heart disease mortality in US adults. JAMA. 1999;281:1291–1297. [PubMed]
43. Malmberg K, Yusuf S, Gerstein HC, Brown J, Zhao F, et al. Impact of diabetes on long-term prognosis in patients with unstable angina and non-Q-wave myocardial infarction: results of the OASIS (Organization to Assess Strategies for Ischemic Syndromes) Registry. Circulation. 2000;102:1014–1019. [PubMed]
44. Fujiwara K, Hiasa Y, Takahashi T, Yamaguchi K, Ogura R, et al. Influence of diabetes mellitus on outcome in the era of primary stenting for acute myocardial infarction. Circ J. 2002;66:800–804. [PubMed]
45. Svensson AM, Abrahamsson P, McGuire DK, Dellborg M. Influence of diabetes on long-term outcome among unselected patients with acute coronary events. Scand Cardiovasc J. 2004;38:229–234. [PubMed]
46. Madsen MM, Busk M, Sondergaard HM, Bottcher M, Mortensen LS, et al. Does diabetes mellitus abolish the beneficial effect of primary coronary angioplasty on long-term risk of reinfarction after acute ST-segment elevation myocardial infarction compared with fibrinolysis? (A DANAMI-2 substudy). Am J Cardiol. 2005;96:1469–1475. [PubMed]
47. Bouhanick B, Cambou JP, Ferrieres J, Amelineau E, Guize L. Characteristics and six-month outcomes in a cohort of 8288 diabetic and non-diabetic patients with previous history of acute coronary syndrome or stroke: the French PREVENIR 3 survey. Diabetes Metab. 2006;32:460–466. [PubMed]
48. Hansen HH, Joensen AM, Riahi S, Malczynski J, Molenberg D, et al. Short and long-term outcome in diabetic patients with acute myocardial infarction in the invasive era. Scand Cardiovasc J. 2007;41:19–24. [PubMed]
49. Cubbon RM, Wheatcroft SB, Grant PJ, Gale CP, Barth JH, et al. Temporal trends in mortality of patients with diabetes mellitus suffering acute myocardial infarction: a comparison of over 3000 patients between 1995 and 2003. Eur Heart J. 2007;28:540–545. [PubMed]
50. Franklin K, Goldberg RJ, Spencer F, Klein W, Budaj A, et al. Implications of diabetes in patients with acute coronary syndromes. The Global Registry of Acute Coronary Events. Arch Intern Med. 2004;164:1457–1463. [PubMed]
51. Mak KH, Moliterno DJ, Granger CB, Miller DP, White HD, et al. Influence of diabetes mellitus on clinical outcome in the thrombolytic era of acute myocardial infarction. GUSTO-I Investigators. Global Utilization of Streptokinase and Tissue Plasminogen Activator for Occluded Coronary Arteries. J Am Coll Cardiol. 1997;30:171–179. [PubMed]
52. Donahoe SM, Stewart GC, McCabe CH, Mohanavelu S, Murphy SA, et al. Diabetes and mortality following acute coronary syndromes. JAMA. 2007;298:765–775. [PubMed]
53. Flavell DM, Pineda TI, Jamshidi Y, Evans D, Diamond JR, et al. Variation in the PPARalpha gene is associated with altered function in vitro and plasma lipid concentrations in Type II diabetic subjects. Diabetologia. 2000;43:673–680. [PubMed]
54. Tai ES, Corella D, Demissie S, Cupples LA, Coltell O, et al. Polyunsaturated fatty acids interact with the PPARA-L162V polymorphism to affect plasma triglyceride and apolipoprotein C-III concentrations in the Framingham Heart Study. J Nutr. 2005;135:397–403. [PubMed]
55. Tai ES, Demissie S, Cupples LA, Corella D, Wilson PW, et al. Association between the PPARA L162V polymorphism and plasma lipid levels: the Framingham Offspring Study. Arterioscler Thromb Vasc Biol. 2002;22:805–810. [PubMed]
56. Chen S, Tsybouleva N, Ballantyne CM, Gotto AM, Jr, Marian AJ. Effects of PPARalpha, gamma and delta haplotypes on plasma levels of lipids, severity and progression of coronary atherosclerosis and response to statin therapy in the lipoprotein coronary atherosclerosis study. Pharmacogenetics. 2004;14:61–71. [PubMed]
57. Jamshidi Y, Montgomery HE, Hense HW, Myerson SG, Torra IP, et al. Peroxisome proliferator–activated receptor alpha gene regulates left ventricular growth in response to exercise and hypertension. Circulation. 2002;105:950–955. [PubMed]
58. Doney AS, Fischer B, Lee SP, Morris AD, Leese G, et al. Association of common variation in the PPARA gene with incident myocardial infarction in individuals with type 2 diabetes: A Go-DARTS study. Nucl Recept. 2005;3:4. [PMC free article] [PubMed]
59. Nissen SE, Wolski K, Topol EJ. Effect of muraglitazar on death and major adverse cardiovascular events in patients with type 2 diabetes mellitus. JAMA. 2005;294:2581–2586. [PubMed]
60. Nissen SE, Wolski K. Effect of rosiglitazone on the risk of myocardial infarction and death from cardiovascular causes. N Engl J Med. 2007;356:2457–2471. [PubMed]
61. Huss JM, Kopp RP, Kelly DP. Peroxisome proliferator-activated receptor coactivator-1alpha (PGC-1alpha) coactivates the cardiac-enriched nuclear receptors estrogen-related receptor-alpha and -gamma. Identification of novel leucine-rich interaction motif within PGC-1alpha. J Biol Chem. 2002;277:40265–40274. [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science