- Research article
- Open Access
Association and interaction of PPAR-complex gene variants with latent traits of left ventricular diastolic function
BMC Medical Genetics volume 11, Article number: 65 (2010)
Abnormalities in myocardial metabolism and/or regulatory genes have been implicated in left ventricular systolic dysfunction. However, the extent to which these modulate left ventricular diastolic function (LVDF) is uncertain.
Independent component analysis was applied to extract latent LVDF traits from 14 measured echocardiography-derived endophenotypes of LVDF in 403 Caucasians. Genetic association was assessed between measured and latent LVDF traits and 64 single nucleotide polymorphisms (SNPs) in three peroxisome proliferator-activated receptor (PPAR)-complex genes involved in the transcriptional regulation of fatty acid metabolism.
By linear regression analysis, 7 SNPs (4 in PPARA, 2 in PPARGC1A, 1 in PPARG) were significantly associated with the latent LVDF trait, whereas a range of 0-4 SNPs were associated with each of the 14 measured echocardiography-derived endophenotypes. Frequency distribution of P values showed a greater proportion of significant associations with the latent LVDF trait than for the measured endophenotypes, suggesting that analyses of the latent trait improved detection of the genetic underpinnings of LVDF. Ridge regression was applied to investigate within-gene and gene-gene interactions. In the within-gene analysis, there were five significant pair-wise interactions in PPARGC1A and none in PPARA or PPARG. In the gene-gene analysis, significant interactions were found between rs4253655 in PPARA and rs1873532 (p = 0.02) and rs7672915 (p = 0.02), both in PPARGC1A, and between rs1151996 in PPARG and rs4697046 in PPARGC1A (p = 0.01).
Myocardial metabolism PPAR-complex genes, including within and between genes interactions, may play an important role modulating left ventricular diastolic function.
Both animal models of left ventricular (LV) pressure-overload and clinical studies in humans with hypertension implicate abnormal myocardial metabolism in the development of hypertensive heart disease (HHD) which is characterized by phenotypes such as LV hypertrophy (LVH), left ventricular diastolic dysfunction (LVDD), left ventricular systolic dysfunction (LVSD), and/or the development of heart failure (HF). Myocardial fatty acid metabolism is a key modulator of HHD phenotypes and peroxisome proliferator-activated receptor (PPAR)-complex genes play a critical role in regulating these metabolic processes [1–3].
The assessment of LV diastolic function (LVDF) by Doppler echocardiography is challenging due to the age-dependency of measurements, the confounding effect of LVH, and the non-linear distribution of many LVDF parameters such that no single measurement has emerged that defines LVDD [4–6]. As LVDD lies on a continuum between normal LVDF and diastolic heart failure (DHF), it represents a relevant clinical diagnosis for several reasons. LVDD, in the presence of normal LV systolic function, occurs in approximately 25% of adults 45 years of age or older and is a multivariate predictor of all-cause mortality after controlling for age, sex, and LVSD (hazards ratio >8). The most severe form of LVDD, DHF, accounts for 40-50% of the approximately 5.7 million patients with HF in the US.
We have previously shown that independent component analysis (ICA)-derived latent traits of HHD can be used to enhance the detection of genetic association. Herein we apply ICA to derive a single latent LVDF trait from multiple echocardiography-derived indices (or "endophenotypes") of LVDF. We tested the hypothesis that PPAR-complex gene variants modulate LVDF traits by examining individual SNP associations and within-gene and gene-gene interaction effects.
I. Study population
Caucasian adults (n = 403) were analyzed from among a multi-racial cohort of consecutive subjects (n = 543) genotyped as part of a prospective genotype-phenotype association study of hypertensive heart disease at Washington University. Subjects exhibited a wide-range of cardiovascular and/or metabolic phenotypes, and also included healthy volunteers. Exclusion criteria included: 1) Hispanic ethnicity; 2) incomplete echocardiogram; 3) LVH due to conditions other than hypertension (i.e, hypertrophic cardiomyopathy); 4) significant valvular heart disease (regurgitation and/or stenosis > mild); or 5) significant systemic disease (i.e. malignancy, creatinine ≥0.22 mmol/dL [2.5 mg/dL]).
Subject demographics were assessed by interview. Heart rate, blood pressure, and anthropomorphic measurements were obtained according to a standard protocol. All subjects had blood drawn for fasting glucose, insulin, and lipids; DNA was extracted from peripheral blood leukocytes. The study was approved by the Human Research Protection Office at Washington University; all subjects provided written informed consent.
The echocardiographic study was performed with an ultrasound system (Acuson-Siemens Sequoia, Mountain View, CA) using a 3.5-MHz array transducer. LV end-systolic and end-diastolic volumes were calculated according to the "method of discs" to derive LVEF (normal LV systolic function defined as LVEF ≥50%). LV mass was measured by the area-length method and indexed by height2.7 (LVM/Ht2.7) to adjust for body habitus. LVH was defined as LVM/Ht2.7 >51 g/m2.7 for men and >49.5 g/m2.7 for women. Pulse-wave Doppler (PWD)-derived transmitral indices were recorded from the four-chamber view at the mitral valve leaflet tips to determine the early diastolic transmitral (E-wave) and atrial (A-wave) velocities (m/s), the E/A wave velocity ratio, E-wave deceleration time (DT, in ms); the isovolumic relaxation time (IVRT, in ms) was measured by continuous-wave Doppler from the apical 5-chamber view. Tissue Doppler Imaging (TDI)-derived early diastolic myocardial velocity was obtained from a 2.5 mm sample at the septal and lateral mitral annulus in the apical four-chamber view (E'sep and E'lat, A'sep and A'lat, respectively, in cm/s) with the average value represented as a "global" value (E'gl and A'gl, in cm/s) [13–15]. The mitral E-wave/E' ratio was calculated to estimate the LV filling pressure (E/E'sep, E/E'lat, and E/E'gl)[16, 17]. All reported measurements represent the average of three consecutive cardiac cycles obtained by a single observer blinded to clinical status. The intra-class correlation coefficients for echocardiographic indices measured in our lab are: 0.75-0.88 for LV structure (i.e., LV mass/volumes) and 0.82-0.97 for PWD- and TDI-derived indices of LV diastolic function (i.e., E-, A-waves, DT, IVRT, and E').
III. PPAR-complex genes genotyping and quality control
SNPs in three PPAR-complex genes (i.e., PPARA, PPARG, and PPARGC1A) were genotyped using Illumina BeadArray technology; 22 PPARA, 17 PPARG, and 39 PPARGC1A were selected to provide dense coverage and for compatibility in the multiplex reactions. Excluded from analysis were SNPs a) that deviated from Hardy-Weinberg equilibrium (HWE, p < 0.01), b) with minor allele frequency ≤5%, and/or c) with genotype call rates ≤90%. Maximal tolerated individual missing rate was 50%.
IV. Statistical analysis
A. Extract latent LVDF trait by ICA
Initial quality control was performed on each of the echocardiography-derived endophenotypes by univariate analysis; those with an absolute skewness >1.5 or absolute kurtosis >2.0 were log-transformed for the factor analysis. Potential confounding effects were adjusted by regression of the values over a cubic age polynomial within sex group; residuals from the regression were used for ICA and further analyses.
We used a freely available implementation called FastICA (version 1.9) available in R (v. 2.7.0) to analyze the residuals of a matrix of the echocardiography-derived endophenotypes by pre-specified numbers of latent components[8, 18, 19]. For each component, the output from ICA consists of the extracted independent component (IC) represented by column vectors of the matrix of loadings on the echocardiography-derived endophenotypes and the coefficients of the extracted IC for each subject. The coefficients were then used to represent the individual latent LVDF trait values in subsequent analyses. Spearman's correlation between several traditional risk factors (systolic and diastolic blood pressure, body mass index, insulin, total cholesterol, triglycerides, LDL-C, HDL-C) and the latent LVDF trait were performed.
B. Select latent LVDF trait from extracted ICs
Loadings of each IC were examined for overrepresentation of the echocardiography-derived endophenotypes, which may reflect the importance of a particular IC. Individual IC's were treated as a potential latent LVDF trait and its clinical characteristics were characterized by examining risks of LVDF in cohorts defined by the median of corresponding IC coefficients (groups assigned as "High-Risk" and "Low-Risk"). The means of each of the 14 echocardiography-derived phenotypes were compared to gauge the IC's capability of separating risks. Three expert echocardiographers (JMJ, LdlF, VGD-R) independently assessed the distributional characteristics of these measures and selected the IC that yielded the most clinically relevant separation between groups; by consensus the IC was selected to represent the latent LVDF trait for subsequent analysis.
C. Genetic association analysis of the selected latent LVDF trait
HWE was tested by χ2 test with 1 degree of freedom to identify SNPs with potential systematic bias resulting from genotyping errors. Linear regression analysis was performed to test for associations between the selected latent LVDF trait (represented by the IC coefficient) and the PPAR-complex candidate SNPs. In each regression model, covariates for age, sex, hypertension status, and LVM/Ht2.7 were included. All regression analyses of genetic association were performed by SAS (v. 9.1.3, SAS Institute Inc., Cary, NC).
D. Explore potential SNP-SNP (within-gene and gene-gene) interactions
For identification of SNP-SNP interactions of potential importance to LVDD, a two-step procedure was applied: Step 1) detection of interactions within individual candidate genes (within-gene), and Step 2) identification of between-gene interactions (gene-gene). To deal with the large number of variables representing individual SNPs, their interactions, and their potential for collinearity, ridge regression analysis was performed (R package "penalized" version 0.9); ridge regression effectively eliminates many variables that are irrelevant to the trait of interest and effectively eliminates the collinearity problem[20, 21]. In Step 1, for each candidate gene, ridge regression was applied to a model including all SNPs within that gene and their pair-wise interactions. SNPs that had main and/or interaction effects selected by the ridge regression were retained. In Step 2, all retained SNPs and their pairwise interactions (including those between genes) were entered into a general model. Ridge regression was again applied to the general model to identify a final set of important SNPs and interactions that together is significantly associated with the LVDD trait. Finally, conventional linear regression was applied to the selected SNPs and interactions to derive the significance levels of the associations.
E. Statistical power and adjustment for multiple testing
Quanto (v.1.2, http://hydra.usc.edu/gxe) was used to estimate statistical power (set at 80% at a significance level of 0.05 with a 2-sided alternate hypothesis). The power calculation assumed that an associated SNP has minor allele frequency of 0.2 and a locus-specific heritability of at least 2%. To adjust for multiple testing, the method of false discovery rate was applied in SAS.
I. Clinical and echocardiographic characteristics
The clinical and echocardiographic characteristics are shown in Tables 1 and 2, respectively. The majority of subjects had normal LV systolic function; mild LV systolic dysfunction (LVEF 40-49%) was present in only 1.9% of the population; 17% of subjects had LVH. Distributional characteristics of the variables are shown in Additional File 1, Table S1.
II. Extraction and identification of latent LVDF components
ICA analyses were performed by extracting 2, 3, and 6 components. Each extracted latent LVDF component was denoted by a letter "E", followed by a first digit which denoted the predefined number of components, and a secondary digit which denoted the individual component. For example, E21 refers to the first extracted component of a total of 2 latent components (Figure 1; figures for all extracted components in Additional File 1, Figures S1-S3). The collective loading pattern across the echocardiography-derived endophenotypes reflects the capability of a component to summarize a common source of regulation or an underlying interaction observed among genes. The component E61 had the highest loadings on E'sep, E'lat, E'gl, E/E'sep, E/E'lat, and E/E'gl compared to all other components in the loading of plots (Figure 1) and was thus selected as the latent LVDF endophenotype (Additional File 1, Table S2). The Spearman's correlation coefficients between several traditional risk factors (systolic and diastolic blood pressure, body mass index, insulin, total cholesterol, triglycerides, LDL-C, HDL-C) and the latent LVDF trait were determined (Additional File 1, Table S3). With the exception of LDL-C and total cholesterol, the remaining risk factors were significantly correlated with the latent LVDF trait.
III. Ascertainment of quality of selected latent LVDF trait
The entire cohort was grouped as "High Risk" or "Low Risk" based on the medians of the individual latent LVDF trait values (i.e. individual ICs corresponding to each extracted component). Clinical characteristics and echocardiographic measurements in the "High-risk" and "Low-risk" groups are shown in Table 3 and Table 4, respectively. E61 resulted in the highest loadings on the TDI-derived echocardiographic variables, and thus, resulted in the most distinct separation of these LVDF echocardiography-derived measures (Additional File 1, Tables S4-6).
IV. Genetic association analysis of the latent LVDF trait
After applying exclusion criteria, the resultant dataset consisted of 64 SNPs (15 PPARA, 14 PPARG, and 34 PPARGC1A). Individual locus information including allele frequencies and HWE p-values for all 78 SNPs (Additional File 1, Tables S6A-C) and the Haploview linkage disequilibrium displays for each gene (Additional File 1, Figures S4-6) are shown in the supplement. By linear regression analysis of single SNPs (after adjustment for covariates), 7 SNPs (4 in PPARA, 2 in PPARGC1A, 1 in PPARG) were significantly associated with the latent LVDF trait (E61), whereas a range of 0 to 4 SNPs were associated with each of the 14 echocardiography-derived endophenotypes (Table 5). The frequency distribution of the P values derived from the latent LVDF trait exhibited a higher degree of statistical significance than those derived from the echocardiography-derived endophenotypes (Additional File 1, Figure S7).
V. Within-gene and gene-gene interactions
In the within-gene analysis, significant SNPs (8 in PPARA, 12 in PPARGC1A, and 4 in PPARG) for each of the three genes were retained for the final best-fit model (log likelihood: -276.1, -281.2 and -281.6, respectively; Table 6 and Figure 2). PPARGC1A SNPs rs12500214 significantly interacted with both rs2970847 and rs7672915 (p = 0.02 and 0.009, respectively); rs768695 significantly interacted with both rs2970847 and rs2970853 (both p = 0.03); and rs4235308 significantly interacted with rs7672915 (p = 0.007). No significant within-gene interactions were found among PPARA and PPARG SNPs.
In the gene-gene analysis, rs4253655 in PPARA significantly interacted with rs1873532 (p = 0.02) and rs7672915 (p = 0.02), both in PPARGC1A. Significant interactions were also found between rs1151996 in PPARG and rs4697046 in PPARGC1A (p = 0.01). No gene-gene interactions were found between PPARA and PPARG genes.
VI. Power calculation and false discovery rate
The study sample of 403 subjects provided adequate power to detect association of any SNP with a heritability ≥0.02 and an allele frequency ≥0.2 (power ≥81%, significance level α = 0.05). Although there was a significant association for several individual SNPs in single SNP analysis, none remained significant after adjustment for multiple testing by the False Discovery Rate.
In this study, we show genetic association between echocardiography-derived endophenotypes of LVDF and SNPs in three PPAR-complex genes involved in the transcriptional regulation of fatty acid metabolism. ICA was used to extract latent LVDF traits from 14 echocardiography-derived endophenotypes used clinically to evaluate diastolic function. We found that the latent LVDF trait identified relevant profiles of echocardiography-derived endophenotypes that could be clinically recognized as representing individuals of either High- or Low-Risk of LVDD. By use of linear regression analysis, 7 SNPs (4 in PPARA, 2 in PPARGC1A, 1 in PPARG) were significantly associated with the latent LVDF trait. In PPARA, 2 SNPs were clustered around intron 6 and two were in introns 1A and 2A, adjacent to alternate promoters. The significant PPARGC1A SNPs were located in introns 1 and 7. The single significant PPARG SNP (rs3856806) is a synonymous coding SNP in Exon 7 (H447H).
Complex phenotypes such as LVDF arise from multiple gene-gene and gene-environment interactions, each contributing a small effect to the overall expression of the trait. A consideration of epistasis, or gene-gene interactions, may identify a portion of the unexplained risk (i.e., "missing heritability" ) noted in the overwhelming majority of genetic studies of complex diseases. It is therefore of particular importance to investigate potential gene-gene interactions in the complex mechanisms responsible for LVDF. In the present study, an exploratory analysis of SNP-SNP interactions by ridge regression yielded potentially important findings. First, significant associations detected by modeling both main effects and interactions were largely different from those by main effects alone; with the only exception of rs3856806, which had a relatively strong marginal effect, suggesting that the single-SNP scan approach may only be appropriate for variants with strong signals. Second, all gene-gene interactions involved PPARGC1A, probably reflecting the important regulatory role of this gene. Thus, metabolic modulation of left ventricular diastolic function by PPAR-complex genes may be an important mechanism in the development of LVDD and may contribute to the pathogenesis of DHF. This is the first study to show a genetic association between LVDF and PPAR-complex genes in humans.
Regulation of myocardial metabolism and cardiac structure/function by PPAR-complex
Alterations in myocardial fatty acid metabolism have been shown in both animal models and in humans to be an important determinant of the presence and development of hypertensive heart disease-related traits including hypertension, LVH, LVSD, and LVDD[3, 24]. The normal fasting adult mammalian heart derives approximately two-thirds of ATP from the oxidation of fatty acids. Under normal physiologic conditions, the heart is able to switch the energy substrate utilization from fatty acid to carbohydrates during the fed state, however under pathologic conditions the heart loses this metabolic plasticity. For example, our group has previously shown that fatty acid metabolism is decreased in patients with heart failure and is a strong predictor of LV mass in patients with hypertension and in normotensive controls. Metabolic modulation of the heart occurs primarily at the transcriptional level through the coordinated regulation of enzymes and proteins in specific metabolic pathways. PPAR-complex genes, transcription factors and coactivators known to regulate the expression of fatty acid transport and oxidation genes, have been shown in elegant transgenic models to modulate the pathophysiology of LV systolic and diastolic dysfunction[3, 28–36]. The PPARA gene has been identified as a "master switch" for the metabolic remodeling of the heart, particularly when coactivated by PPARGC1A [36–39]. PPARA expression is down-regulated in the adapted, hypertrophied pressure-overloaded heart and reactivation of PPARA by an agonist is associated with contractile dysfunction. Cardiac-restricted overexpression of PPARA results in hypertrophy, the activation of gene markers of pathologic hypertrophic growth, and in systolic dysfunction. PPARG, on the other hand, plays a critical role in modulating the substrate environment of the heart by its actions in adipocytes. Myocardial fatty acid metabolism has also been more specifically associated with indices of LV diastolic function in animal models and in humans[1–3, 42, 43]. In an animal study, PPARA reduced myocardial fibrosis and prevented the development of diastolic dysfunction. Ciglitazole (a PPARG agonist) has also been shown to attenuate LVDD in pressure-overloaded rats. PPARGC1A is a key regulator of cardiac energy metabolism; overexpression of PPARGC1A in the murine heart leads to a modest increase in mitochondrial number, derangements of mitochondrial ultrastructure, and development of cardiomyopathy[32, 46]. Thus, there is ample evidence to support a critical role for PPAR-complex genes in regulating cardiac structure and function. However, the precise molecular mechanisms by which these alterations may modulate LVDF remain uncertain.
Several lines of evidence support the existence of genetic association between LVDF and PPAR-complex genes in humans. Recent reports have shown associations between cardiovascular traits and 4 of the 19 SNPs indentified in the present study. The most widely published among these, a synonymous coding PPARGC1A SNP rs2970847 (T394T), has been associated individually, in combination with other PPARGC1A SNPs in haplotype blocks, and in gene-gene interactions with relevant traits including diabetes risk, glucose uptake, obesity, and DNA damage in a variety of populations [47–54]. Intriguingly, this SNP has also been associated with non-hypertensive LVH in a cohort including 270 hypertrophic cardiomyopathy and 2486 hypertensive patients (with and without LVH), yielding an odds ratio 1.49 (95% confidence interval, 1.15-1.98). The significant PPARGC1A SNP rs4697046 has also been linked with plasma glucose and DNA damage. In our study, this SNP was found to interact with two other PPARGC1A SNPs, although these SNPs were not found to share significant linkage disequilibrium (data not shown), thus implicating cis-acting regulatory elements.
This study also identified two SNPs (rs3856806 in PPARG and rs4253623 in PPARA) as associated with MI risk [57–59]. The rare allele of PPARG SNP rs3856806 has been linked with not only coronary artery disease progression, but also with pro-inflammatory cytokines (MMP9 and TNFα), plasma homocysteine levels, and obesity [60–62]. The minor allele of PPARA alternate promoter SNP, rs4253623, along with another significant PPARA promoter SNP rs4253681, is clustered adjacent to alternatively-spliced untranslated exons 1A and 2A and may play a role in regulating PPARA gene expression by mediating the expression of the 6 alternative spliced variants resulting from 4 different promoters. Thus, it is plausible that variants in PPAR-complex genes may modulate LVDF-related traits in humans. However, all these possible mechanisms deserve to be explored further in animal models.
Utility of ICA for study of genetic association
The present study was performed using state-of-the-art analysis techniques in genetic epidemiology and represents a continuation of our study of the genetics of hypertensive heart disease through analysis of endophenotypes, which lay proximal along the pathway from observed clinical phenotype to genotype. The latent factor analysis by ICA was used to address the multidimensional data in which non-Gaussian structure such as clustering and independence represent important components. Although several echocardiographic measures are considered jointly to clinically ascertain LVDF, the precise diagnosis of LVDD remains controversial and there currently exists no universally accepted diagnostic criteria[4, 64, 65]. Furthermore, many echocardiographic endophenotypes are characterized by age-dependency and non-linear associations with disease severity, presenting additional challenges for genetic association studies of this common, complex disease. The ICA trait is extracted from a panel of measured echocardiography-derived endophenotypes, which themselves represent traits more proximal to the genetic underpinnings, and as such has emerged as an important tool in the study of complex diseases[66, 67]. ICA is well-suited to the study of LVDF because it effectively reduces the number of dimensions by identifying linear representations of the original endophenotypes. The ICA-derived latent LVDF trait showed significant associations with 7 PPAR-complex SNPs, whereas a range of 0 to 4 SNPs were associated with each of the 14 echocardiography-derived endophenotypes. The frequency distribution of P values showed a greater proportion of significant associations with the latent LVDF trait than for the measured echocardiography-derived endophenotypes, suggesting that analyses of the latent trait improved detection of the genetic underpinnings of LVDF. Thus, use of ICA-derived latent LVDF trait improved our ability to detect the genetic underpinning of LVDD, as our group and others have previously shown for other echocardiography-derived endophenotypes[8, 68].
The ridge regression method was used to handle collinearity among SNPs, and thus allowed for testing many SNPs simultaneously. Since these three PPAR-complex genes are all related to fatty acid metabolism, and LVDD is a complex disease, epistasis may contribute to LVDF. Thus, there was a real possibility of collinearity among SNPs, a result of multiple makers exhibiting strong linkage disequilibrium (LD) in a single genomic region. The presence of LD between SNPs at neighboring loci can make it difficult to distinguish functionally relevant variations from nonfunctional variations. However, the application of ridge regression to the evaluation of SNP-SNP interactions provided a robust method to control for the effect of collinearity.
Finally, we investigated within-gene and gene-gene SNP-SNP interactions using a two-stage procedure as outlined in the methods. This led to interesting findings of interaction effects that were not detected by testing for marginal effects alone. This last issue is of particular importance as supported by recent studies in yeast showing that "pure" interactions (with little or negligible marginal effects) are both real and important.
ICA-derived latent LVDF traits improved the ability to detect genetic underpinnings of LVDF. Significant genetic associations were found between PPAR-complex gene variants and LVDF, suggesting that these genes may be involved in the pathogenesis of this complex disease. Although these findings of association are suggestive, these results need to be validated by mechanistic studies in animal models. Pharmacologic modulation of myocardial substrate environment, PPAR-complex genes, and/or both may be an important target to restore myocardial metabolic function and thus normalization of left ventricular diastolic function.
diastolic heart failure
hypertensive heart disease
independent component analysis
isovolumic relaxation time
left ventricular diastolic function
Left Ventricular hypertrophy
left ventricular systolic dysfunction
peroxisome proliferator-activated receptor
single nucleotide polymorphisms
Tissue Doppler Imaging.
Zile MR, Brutsaert DL: New concepts in diastolic dysfunction and diastolic heart failure: Part II: causal mechanisms and treatment. Circulation. 2002, 105 (12): 1503-1508. 10.1161/hc1202.105290.
Tian R, Nascimben L, Ingwall JS, Lorell BH: Failure to maintain a low ADP concentration impairs diastolic function in hypertrophied rat hearts. Circulation. 1997, 96 (4): 1313-1319.
Christoffersen C, Bollano E, Lindegaard ML, Bartels ED, Goetze JP, Andersen CB, Nielsen LB: Cardiac lipid accumulation associated with diastolic dysfunction in obese mice. Endocrinology. 2003, 144 (8): 3483-3490. 10.1210/en.2003-0242.
Nagueh SF, Appleton CP, Gillebert TC, Marino PN, Oh JK, Smiseth OA, Waggoner AD, Flachskampf FA, Pellikka PA, Evangelista A: Recommendations for the evaluation of left ventricular diastolic function by echocardiography. J Am Soc Echocardiogr. 2009, 22 (2): 107-133. 10.1016/j.echo.2008.11.023.
How to diagnose diastolic heart failure. European Study Group on Diastolic Heart Failure. Eur Heart J. 1998, 19 (7): 990-1003. 10.1053/euhj.1998.1057.
Redfield MM, Jacobsen SJ, Burnett JC, Mahoney DW, Bailey KR, Rodeheffer RJ: Burden of systolic and diastolic ventricular dysfunction in the community: appreciating the scope of the heart failure epidemic. JAMA. 2003, 289 (2): 194-202. 10.1001/jama.289.2.194.
Lloyd-Jones D, Adams R, Carnethon M, De Simone G, Ferguson TB, Flegal K, Ford E, Furie K, Go A, Greenlund K, et al: Heart Disease and Stroke Statistics--2009 Update. A Report From the American Heart Association Statistics Committee and Stroke Statistics Subcommittee. Circulation. 2008, 119 (3): e21-181.
Gu CC, Flores HR, de las Fuentes L, Davila-Roman VG: Enhanced detection of genetic association of hypertensive heart disease by analysis of latent phenotypes. Genet Epidemiol. 2008, 32 (6): 528-538. 10.1002/gepi.20326.
de las Fuentes L, Brown AL, Mathews SJ, Waggoner AD, Soto PF, Gropler RJ, Davila-Roman VG: Metabolic syndrome is associated with abnormal left ventricular diastolic function independent of left ventricular mass. Eur Heart J. 2007, 28 (5): 553-559. 10.1093/eurheartj/ehl526.
Schiller NB, Shah PM, Crawford M, DeMaria A, Devereux R, Feigenbaum H, Gutgesell H, Reichek N, Sahn D, Schnittger I, et al: Recommendations for quantitation of the left ventricle by two-dimensional echocardiography. American Society of Echocardiography Committee on Standards, Subcommittee on Quantitation of Two-Dimensional Echocardiograms. J Am Soc Echocardiogr. 1989, 2 (5): 358-367.
de Simone G, Kizer JR, Chinali M, Roman MJ, Bella JN, Best LG, Lee ET, Devereux RB: Normalization for body size and population-attributable risk of left ventricular hypertrophy: the Strong Heart Study. Am J Hypertens. 2005, 18 (2 Pt 1): 191-196. 10.1016/j.amjhyper.2004.08.032.
de Simone G, Daniels SR, Devereux RB, Meyer RA, Roman MJ, de Divitiis O, Alderman MH: Left ventricular mass and body size in normotensive children and adults: assessment of allometric relations and impact of overweight. J Am Coll Cardiol. 1992, 20 (5): 1251-1260. 10.1016/0735-1097(92)90385-Z.
Quinones MA, Otto CM, Stoddard M, Waggoner A, Zoghbi WA: Recommendations for quantification of Doppler echocardiography: a report from the Doppler Quantification Task Force of the Nomenclature and Standards Committee of the American Society of Echocardiography. J Am Soc Echocardiogr. 2002, 15 (2): 167-184. 10.1067/mje.2002.120202.
Wang M, Yip GW, Wang AY, Zhang Y, Ho PY, Tse MK, Lam PK, Sanderson JE: Peak early diastolic mitral annulus velocity by tissue Doppler imaging adds independent and incremental prognostic value. J Am Coll Cardiol. 2003, 41 (5): 820-826. 10.1016/S0735-1097(02)02921-2.
Dumesnil JG, Paulin C, Pibarot P, Coulombe D, Arsenault M: Mitral annulus velocities by Doppler tissue imaging: practical implications with regard to preload alterations, sample position, and normal values. J Am Soc Echocardiogr. 2002, 15 (10 Pt 2): 1226-1231.
Nagueh SF, Middleton KJ, Kopelen HA, Zoghbi WA, Quinones MA: Doppler tissue imaging: a noninvasive technique for evaluation of left ventricular relaxation and estimation of filling pressures. J Am Coll Cardiol. 1997, 30 (6): 1527-1533. 10.1016/S0735-1097(97)00344-6.
Ommen SR, Nishimura RA, Appleton CP, Miller FA, Oh JK, Redfield MM, Tajik AJ: Clinical utility of Doppler echocardiography and tissue Doppler imaging in the estimation of left ventricular filling pressures: A comparative simultaneous Doppler-catheterization study. Circulation. 2000, 102 (15): 1788-1794.
Hyvarinen A, Karhunen J, Oja E: Independent Component Analysis. 2001, J Wiley New York
Hyvarinen A, Oja E: Independent component analysis: algorithms and applications. Neural Netw. 2000, 13 (4-5): 411-430. 10.1016/S0893-6080(00)00026-5.
Team RDC: R: A language and environment for statistical computing. 2005, R Foundation for Statistical Computing, Vienna, Austria, ISBN 3-900051-07-0
Malo N, Libiger O, Schork NJ: Accommodating linkage disequilibrium in genetic-association analyses via ridge regression. Am J Hum Genet. 2008, 82 (2): 375-385. 10.1016/j.ajhg.2007.10.012.
Gauderman W, Morrison J: QUANTO 1.1: A computer program for power and sample size calculations for genetic-epidemiololgystudies. 2006, [http://hydra.usc.edu/gxe/]
Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, Hunter DJ, McCarthy MI, Ramos EM, Cardon LR, Chakravarti A, et al: Finding the missing heritability of complex diseases. Nature. 2009, 461 (7265): 747-753. 10.1038/nature08494.
Chiu HC, Kovacs A, Blanton RM, Han X, Courtois M, Weinheimer CJ, Yamada KA, Brunet S, Xu H, Nerbonne JM, et al: Transgenic expression of fatty acid transport protein 1 in the heart causes lipotoxic cardiomyopathy. Circ Res. 2005, 96 (2): 225-233. 10.1161/01.RES.0000154079.20681.B9.
Bing RJ, Siegel A, Ungar I, Gilbert M: Metabolism of the human heart. II. Studies on fat, ketone and amino acid metabolism. Am J Med. 1954, 16 (4): 504-515. 10.1016/0002-9343(54)90365-4.
Goodwin GW, Taylor CS, Taegtmeyer H: Regulation of energy metabolism of the heart during acute increase in heart work. J Biol Chem. 1998, 273 (45): 29530-29539. 10.1074/jbc.273.45.29530.
de las Fuentes L, Herrero P, Peterson LR, Kelly DP, Gropler RJ, Davila-Roman VG, (eds): Myocardial fatty acid metabolism: independent predictor of left ventricular mass in hypertensive heart disease. 2003
Chiu HC, Kovacs A, Ford DA, Hsu FF, Garcia R, Herrero P, Saffitz JE, Schaffer JE: A novel mouse model of lipotoxic cardiomyopathy. J Clin Invest. 2001, 107 (7): 813-822. 10.1172/JCI10947.
Zhou YT, Grayburn P, Karim A, Shimabukuro M, Higa M, Baetens D, Orci L, Unger RH: Lipotoxic heart disease in obese rats: implications for human obesity. Proc Natl Acad Sci USA. 2000, 97 (4): 1784-1789. 10.1073/pnas.97.4.1784.
Shimabukuro M, Higa M, Zhou YT, Wang MY, Newgard CB, Unger RH: Lipoapoptosis in beta-cells of obese prediabetic fa/fa rats. Role of serine palmitoyltransferase overexpression. J Biol Chem. 1998, 273 (49): 32487-32490. 10.1074/jbc.273.49.32487.
Shimabukuro M, Zhou YT, Levi M, Unger RH: Fatty acid-induced beta cell apoptosis: a link between obesity and diabetes. Proc Natl Acad Sci USA. 1998, 95 (5): 2498-2502. 10.1073/pnas.95.5.2498.
Finck BN, Lehman JJ, Leone TC, Welch MJ, Bennett MJ, Kovacs A, Han X, Gross RW, Kozak R, Lopaschuk GD, et al: The cardiac phenotype induced by PPARA lpha overexpression mimics that caused by diabetes mellitus. J Clin Invest. 2002, 109 (1): 121-130.
Finck BN, Lehman JJ, Leone TC, Welch MJ, Bennett MJ, Kovacs A, Han X, Gross RW, Kozak R, Lopaschuk GD, et al: The cardiac phenotype induced by PPARA lpha overexpression mimics that caused by diabetes mellitus. Journal of Clinical Investigation. 2002, 109 (1): 121-130.
Lehman JJ, Barger PM, Kovacs A, Saffitz JE, Medeiros DM, Kelly DP: Peroxisome proliferator-activated receptor gamma coactivator-1 promotes cardiac mitochondrial biogenesis. J Clin Invest. 2000, 106 (7): 847-856. 10.1172/JCI10268.
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 (12): 1723-1730. 10.1172/JCI9056.
Barger PM, Kelly DP: PPAR signaling in the control of cardiac energy metabolism. Trends Cardiovasc Med. 2000, 10 (6): 238-245. 10.1016/S1050-1738(00)00077-3.
Taegtmeyer H: Switching metabolic genes to build a better heart. Circulation. 2002, 106 (16): 2043-2045. 10.1161/01.CIR.0000036760.42319.3F.
Kelly DP: PPARs of the heart: three is a crowd. Circ Res. 2003, 92 (5): 482-484. 10.1161/01.RES.0000064382.46274.95.
Vega RB, Huss JM, Kelly DP: The coactivator PGC-1 cooperates with peroxisome proliferator-activated receptor alpha in transcriptional control of nuclear genes encoding mitochondrial fatty acid oxidation enzymes. Mol Cell Biol. 2000, 20 (5): 1868-1876. 10.1128/MCB.20.5.1868-1876.2000.
Young ME, Laws FA, Goodwin GW, Taegtmeyer H: Reactivation of peroxisome proliferator-activated receptor alpha is associated with contractile dysfunction in hypertrophied rat heart. J Biol Chem. 2001, 276 (48): 44390-44395. 10.1074/jbc.M103826200.
Semple RK, Chatterjee VK, O'Rahilly S: PPAR gamma and human metabolic disease. J Clin Invest. 2006, 116 (3): 581-589. 10.1172/JCI28003.
Leichman JG, Aguilar D, King TM, Vlada A, Reyes M, Taegtmeyer H: Association of plasma free fatty acids and left ventricular diastolic function in patients with clinically severe obesity. Am J Clin Nutr. 2006, 84 (2): 336-341.
Nakayama H, Morozumi T, Nanto S, Shimonagata T, Ohara T, Takano Y, Kotani J, Watanabe T, Fujita M, Nishio M, et al: Abnormal myocardial free fatty acid utilization deteriorates with morphological changes in the hypertensive heart. Jpn Circ J. 2001, 65 (9): 783-787. 10.1253/jcj.65.783.
Ogata T, Miyauchi T, Sakai S, Takanashi M, Irukayama-Tomobe Y, Yamaguchi I: Myocardial fibrosis and diastolic dysfunction in deoxycorticosterone acetate-salt hypertensive rats is ameliorated by the peroxisome proliferator-activated receptor-alpha activator fenofibrate, partly by suppressing inflammatory responses associated with the nuclear factor-kappa-B pathway. J Am Coll Cardiol. 2004, 43 (8): 1481-1488. 10.1016/j.jacc.2003.11.043.
Henderson BC, Sen U, Reynolds C, Moshal KS, Ovechkin A, Tyagi N, Kartha GK, Rodriguez WE, Tyagi SC: Reversal of systemic hypertension-associated cardiac remodeling in chronic pressure overload myocardium by ciglitazone. Int J Biol Sci. 2007, 3 (6): 385-392.
Russell LK, Mansfield CM, Lehman JJ, Kovacs A, Courtois M, Saffitz JE, Medeiros DM, Valencik ML, McDonald JA, Kelly DP: Cardiac-specific induction of the transcriptional coactivator peroxisome proliferator-activated receptor gamma coactivator-1alpha promotes mitochondrial biogenesis and reversible cardiomyopathy in a developmental stage-dependent manner. Circ Res. 2004, 94 (4): 525-533. 10.1161/01.RES.0000117088.36577.EB.
Bhat A, Koul A, Rai E, Sharma S, Dhar MK, Bamezai RN: PGC-1alpha Thr394Thr and Gly482Ser variants are significantly associated with T2DM in two North Indian populations: a replicate case-control study. Hum Genet. 2007, 121 (5): 609-614. 10.1007/s00439-007-0352-0.
Rai E, Sharma S, Koul A, Bhat AK, Bhanwer AJ, Bamezai RN: Interaction between the UCP2-866G/A, mtDNA 10398G/A and PGC1alpha p.Thr394Thr and p.Gly482Ser polymorphisms in type 2 diabetes susceptibility in North Indian population. Hum Genet. 2007, 122 (5): 535-540. 10.1007/s00439-007-0421-4.
Su Y, Peng SB, Li ZQ, Huang QY: [Association study between PPARGC1A Thr394Thr/Gly482Ser polymorphisms and type 2 diabetes]. Yi Chuan. 2008, 30 (3): 304-308.
Zhang SL, Lu WS, Yan L, Wu MC, Xu MT, Chen LH, Cheng H: Association between peroxisome proliferator-activated receptor-gamma coactivator-1alpha gene polymorphisms and type 2 diabetes in southern Chinese population: role of altered interaction with myocyte enhancer factor 2C. Chin Med J (Engl). 2007, 120 (21): 1878-1885.
Vimaleswaran KS, Radha V, Anjana M, Deepa R, Ghosh S, Majumder PP, Rao MR, Mohan V: Effect of polymorphisms in the PPARGC1A gene on body fat in Asian Indians. Int J Obes (Lond). 2006, 30 (6): 884-891. 10.1038/sj.ijo.0803228.
Vimaleswaran KS, Radha V, Ghosh S, Majumder PP, Deepa R, Babu HN, Rao MR, Mohan V: Peroxisome proliferator-activated receptor-gamma co-activator-1alpha (PGC-1alpha) gene polymorphisms and their relationship to Type 2 diabetes in Asian Indians. Diabet Med. 2005, 22 (11): 1516-1521. 10.1111/j.1464-5491.2005.01709.x.
Pihlajamaki J, Kinnunen M, Ruotsalainen E, Salmenniemi U, Vauhkonen I, Kuulasmaa T, Kainulainen S, Laakso M: Haplotypes of PPARGC1A are associated with glucose tolerance, body mass index and insulin sensitivity in offspring of patients with type 2 diabetes. Diabetologia. 2005, 48 (7): 1331-1334. 10.1007/s00125-005-1800-9.
Hara K, Tobe K, Okada T, Kadowaki H, Akanuma Y, Ito C, Kimura S, Kadowaki T: A genetic variation in the PGC-1 gene could confer insulin resistance and susceptibility to Type II diabetes. Diabetologia. 2002, 45 (5): 740-743. 10.1007/s00125-002-0803-z.
Wang S, Fu C, Wang H, Shi Y, Xu X, Chen J, Song X, Sun K, Wang J, Fan X, et al: Polymorphisms of the peroxisome proliferator-activated receptor-gamma coactivator-1alpha gene are associated with hypertrophic cardiomyopathy and not with hypertension hypertrophy. Clin Chem Lab Med. 2007, 45 (8): 962-967. 10.1515/CCLM.2007.189.
Lai CQ, Tucker KL, Parnell LD, Adiconis X, Garcia-Bailo B, Griffith J, Meydani M, Ordovas JM: PPARGC1A variation associated with DNA damage, diabetes, and cardiovascular diseases: the Boston Puerto Rican Health Study. Diabetes. 2008, 57 (4): 809-816. 10.2337/db07-1238.
Enquobahrie DA, Smith NL, Bis JC, Carty CL, Rice KM, Lumley T, Hindorff LA, Lemaitre RN, Williams MA, Siscovick DS, et al: Cholesterol ester transfer protein, interleukin-8, peroxisome proliferator activator receptor alpha, and Toll-like receptor 4 genetic variations and risk of incident nonfatal myocardial infarction and ischemic stroke. Am J Cardiol. 2008, 101 (12): 1683-1688. 10.1016/j.amjcard.2008.02.052.
Szeto CC, Chow KM, Poon PY, Kwan BC, Li PK: Peroxisome proliferator-activated receptor-gamma gene polymorphism and risk of cardiovascular disease in patients with diabetic nephropathy. Am J Nephrol. 2008, 28 (5): 715-722. 10.1159/000127452.
Wang XL, Oosterhof J, Duarte N: Peroxisome proliferator-activated receptor gamma C161-->T polymorphism and coronary artery disease. Cardiovasc Res. 1999, 44 (3): 588-594. 10.1016/S0008-6363(99)00256-4.
Rhee EJ, Oh KW, Lee WY, Kim SY, Oh ES, Baek KH, Kang MI, Kim SW: Effects of two common polymorphisms of peroxisome proliferator-activated receptor-gamma gene on metabolic syndrome. Arch Med Res. 2006, 37 (1): 86-94. 10.1016/j.arcmed.2005.04.008.
Golledge J, Norman PE: Relationship between two sequence variations in the gene for peroxisome proliferator-activated receptor-gamma and plasma homocysteine concentration. Health in men study. Hum Genet. 2008, 123 (1): 35-40. 10.1007/s00439-007-0446-8.
Liu Y, Yuan Z, Zhang J, Yin P, Wang D, Wang Y, Kishimoto C, Ma A: PPARG amma gene C161T substitution is associated with reduced risk of coronary artery disease and decreased proinflammatory cytokine expression. Am Heart J. 2007, 154 (4): 718-724. 10.1016/j.ahj.2007.06.009.
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 (2): 235-243. 10.1016/S0006-291X(03)00731-9.
Bursi F, Weston SA, Redfield MM, Jacobsen SJ, Pakhomov S, Nkomo VT, Meverden RA, Roger VL: Systolic and diastolic heart failure in the community. JAMA. 2006, 296 (18): 2209-2216. 10.1001/jama.296.18.2209.
Pedersen F, Raymond I, Madsen LH, Mehlsen J, Atar D, Hildebrandt P: Echocardiographic indices of left ventricular diastolic dysfunction in 647 individuals with preserved left ventricular systolic function. Eur J Heart Fail. 2004, 6 (4): 439-447. 10.1016/j.ejheart.2003.12.003.
Dick DM, Jones K, Saccone N, Hinrichs A, Wang JC, Goate A, Bierut L, Almasy L, Schuckit M, Hesselbrock V, et al: Endophenotypes successfully lead to gene identification: results from the collaborative study on the genetics of alcoholism. Behav Genet. 2006, 36 (1): 112-126. 10.1007/s10519-005-9001-3.
Gottesman II, Gould TD: The endophenotype concept in psychiatry: etymology and strategic intentions. Am J Psychiatry. 2003, 160 (4): 636-645. 10.1176/appi.ajp.160.4.636.
Diebold B, Delouche A, Abergel E, Raffoul H, Diebold H, Frouin F: Optimization of factor analysis of the left ventricle in echocardiography for detecting wall motion abnormalities. Ultrasound Med Biol. 2005, 31 (12): 1597-1606. 10.1016/j.ultrasmedbio.2005.07.019.
Gerke J, Lorenz K, Cohen B: Genetic interactions between transcription factors cause natural variation in yeast. Science. 2009, 323 (5913): 498-501. 10.1126/science.1166426.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/11/65/prepub
The authors would like to acknowledge Joann Reagan, RN, BSN, and Sharon Heuerman, RN, for their assistance with subject recruitment and phenotypic characterization. This research was supported in part by NIH grants R01HL71782 (V.G.D.-R.), R21HL094668, KL2RR024994 and K12RR023249 (L.d.l.F.), AHA0855626G (C.C.G), UL1RR024992 (Clinical and Translational Sciences Award to Washington University), and from the Robert Wood Johnson (L.d.l.F.) and Barnes-Jewish Hospital (V.G.D.-R.) Foundations. Genotyping for PPARGC1A, PPARA, and PPARG was provided, in part, by an award from the SeattleSNPs Program for Genomic Applications (L.d.l.F.).
The authors declare that they have no competing interests.
JMJ conceived and designed the study, performed statistical analyses, interpret the data, and drafted the manuscript. LdlF conceived and designed the study, performed statistical analyses, warranted the quality of the genetic data, interpret the data, and provided a substantive review of the manuscript. ADW conceived and designed the study, collected and warranted the quality of the echocardiographic phenotype data and provided a substantive review of the manuscript. CG performed statistical analyses, interpret the data, and provided a substantive review of the manuscript. VGD-R conceived and designed the study, interpret the data, offered the guidance on the whole study and provided a substantive review of the manuscript. All authors have read and approved the final manuscript.
Jyh-Ming Jimmy Juang, Lisa de las Fuentes contributed equally to this work.
Electronic supplementary material
Additional file 1: Supplementary Tables and Figures. Table S1. Eleven LVDD-related echocardiographic parameters in the ICA analysis. Table S2. Absolute echocardiographic parameter loadings of 2-, 3-, 6-component ICA. Table S3. Spearman's correlation coefficients between traditional cardiovascular risk factors and the latent LVDF trait. Table S4. Descriptive statistics for primary echocardiographic endophenotypes (2 components set). Table S5. Descriptive statistics for primary echocardiographic endophenotypes (3 components set). Table S6. Descriptive statistics for primary echocardiographic endophenotypes (6 components set). Table S7A-C. Characteristics of 39 SNPs in the PPARGC1A, PPARA, and PPARG genes. Figure S1. Plots of echocardiographic endophenotype loadings of E21 and E22. Figure S2. Plots of echocardiographic endophenotype loadings of E31, E32 and E33. Figure S3. Plots of echocardiographic endophenotype loadings of E61-E66. Figure S4. Haploview LD display of 34 PPARGC1A SNPs. Figure S5. Haploview LD display of 15 PPARA SNPs. Figure S6. Haploview LD display of the 14 PPARG SNPs. Figure S7. Comparison of the distributions of P values between latent LVDD trait (E61) and primary 14 echocardiographic endophenotypes. (DOC 2 MB)
About this article
Cite this article
Juang, JM.J., de las Fuentes, L., Waggoner, A.D. et al. Association and interaction of PPAR-complex gene variants with latent traits of left ventricular diastolic function. BMC Med Genet 11, 65 (2010). https://doi.org/10.1186/1471-2350-11-65
- Independent Component Analysis
- Latent Trait
- Ridge Regression
- Leave Ventricular Systolic Dysfunction
- Leave Ventricular Diastolic Dysfunction