Abstract
Chronic kidney disease (CKD) is well known as a strong risk factor for both end stage renal disease and cardiovascular disease. To clarify the association of polymorphisms in the PPAR genes (PPARD, PPARG, and PPARGC1A) with the risk of CKD in Japanese, we examined this association among the Japanese subjects using the cross-sectional data of J-MICC (Japan Multi-Institutional Collaborative Cohort) Study. The subjects for this analysis were 3,285 men and women, aged 35–69 years, selected from J-MICC Study participants; genotyping was conducted by multiplex polymerase chain reaction-based Invader assay. The prevalence of CKD was determined for CKD stages 3–5 (defined as eGFR < 60 ml/min/1.73 m2). Participants with CKD accounted for 17.3% of the study population. When those with PPARD T-842C T/T were defined as reference, those with PPARD T-842C T/C and C/C demonstrated the OR for CKD of 1.26 (95%CI 1.04–1.53) and 1.31 (95%CI 0.83–2.06), respectively. There were no significant associations between the polymorphisms in other PPAR genes and the risk of CKD. The present study found a significantly increased risk of CKD in those with the C allele of PPARD T-842C, which may suggest the possibility of personalized risk estimation of this life-limiting disease in the near future.
1. Introduction
Chronic kidney disease (CKD) is recently attracting attention as a leading cause of end-stage renal disease (ESRD) and a potential risk factor for cardiovascular disease (CVD). The prevalence of CKD is increasing over time in Japan, reaching about 22% of the adult population in 2002 [1–4]. Although the prevalence of CKD is shown to be specifically higher in Japan compared to those in other countries [5], it is also reported that about 10–15% of adults are affected by this disease in the developed Western countries [2], suggesting that CKD is becoming a general major public health problem worldwide, making its prevention a pressing universal issue.
Meanwhile, metabolic disorders are shown to play an important role in the genesis of CKD mainly through the mechanisms of insulin resistance, resultant hyperinsulinemia, and subsequent increase in plasma renin activity and plasma level of the renal vasoconstrictor angiotensin II [6]. While recent reports suggest that the MetS is an independent predictor of CKD development and progression [7], there are some other possible pathogenic factors causing CKD, such as smoking, hyperuricemia, and homocysteinemia [8]. Peroxisome proliferator-activated receptors (PPARs) are the firstly identified genetic sensor responsive to fatty acid ligands and ligand-activated transcription factors belonging to nuclear hormone receptor superfamily, named after its ability to bind chemicals known to induce peroxisome proliferation [9–11]. There are three PPAR families, PPAR-α, PPAR-γ, and PPAR-δ, each of which has distinct tissue distribution and they are also functionally at odds with each other [9]. These PPARs have been linked to many systemic and cellular functions ranging far beyond the process after which they were initially named; they function as the master regulators of glucose, fatty acid, and lipoprotein metabolism, energy balance, cell proliferation and differentiation, inflammation, and atherosclerosis [11, 12].
Compounds with agonistic activity on PPAR-α, such as fibrates, and PPAR-γ, such as thiazolidinediones, are widely used for the respective treatment of hyperlipidemia and insulin resistance. Although the activators of PPAR-δ are not in clinical use, emerging roles of PPAR-δ in regulating metabolism, inflammation, and antioxidant responses in various cell types and tissues argue in favor of the use of specific agonists to treat some aspects of metabolic dysfunctions [13]. Recently, PPARs have also been increasingly recognized as a key player in the pathogenesis of renal complications associated with metabolic disorders [14]. Meanwhile, there are a number of functional genetic variants in genes encoding PPARs, among which PPARD A65G (Asn163Asn: rs2076167) in exon 7, T-48444C (rs6902123) and T-842C (rs2267668) in exon 3, PPARG C161T (His477His: rs3856806) and Pro12Ala (rs1801282), and PPARGC1A Thr394Thr (rs2970847) and Gly482Ser (rs8192678) are the representative polymorphisms that are shown to modulate the function of PPAR pathways and studied well in association with the risks of various chronic diseases or health conditions [15, 16].
We conducted the Japan Multi-Institutional Collaborative Cohort (J-MICC) Study, a large genome cohort study to confirm and detect gene-environment interactions in lifestyle-related diseases, mainly cancer, launched in 2005 supported by a research grant for Scientific Research on Special Priority Areas of Cancer from the Japanese Ministry of Education, Culture, Sports, Science and Technology [17, 18].
Considering the potential roles of PPARs in various etiologies including metabolic disorders, inflammation, cell proliferation, and atherosclerosis, it would be plausible to hypothesize that genetic polymorphisms modulating the biological functions of PPARs will also affect CKD risk in humans. Accordingly, to clarify the association of polymorphisms in genes encoding PPARs (PPARD, PPARG, and PPARGC1A) with the risk of CKD in Japanese, we examined this association among the Japanese subjects using the cross-sectional data of this J-MICC Study.
2. Methods
2.1. Study Subjects
Subjects were the participants of the J-MICC Study, initially conducted in 10 areas of Japan, in which about 75,000 voluntarily enrolled participants aged 35–69 years provided their blood, health check-up data, and their lifestyle data based on the questionnaire after writting informed consent [17].
In this analysis, 4,519 randomly selected participants (about 500 subjects from each of the 10 areas) were analyzed, for whom 108 selected polymorphisms were genotyped [18]. Firstly, six subjects were excluded due to withdrawal from the study. Serum creatinine (SCr) was measured in 3,327 respondents from 8 areas out of 10. Twelve subjects were excluded because of genotyping failure, and the remaining 3,315 subjects were included in the analyses.
2.2. Evaluation of Lifestyle Exposure
Lifestyle exposures were evaluated with a self-administered questionnaire checked by trained staffs. The questionnaire included items on smoking status, alcohol consumption, and medical history. Smoking status was classified as current, former, or never, and the level of exposure was evaluated in pack-years. Former smokers were defined as people who had quitted smoking for at least 1 year. Alcohol consumption of each type of beverage was determined by average number of drinks per day and then converted into the Japanese sake unit: “gou” (180 mL), which is equivalent to 23 g of ethanol.
2.3. Estimated Glomerular Filtration Rate (eGFR) and Definitions of CKD
Serum creatinine (SCr) was measured in all participants using an enzymatic method. The eGFR of each participant was calculated based on SCr, age, and sex using the following Japanese eGFR equation proposed by the Japanese Society of Nephrology [19]: eGFR (mL/min/1.73 m2) = 194 × SCr (mg/dL)−1.094 × age−0.287 (×0.739 if female). The prevalence of CKD was determined for CKD stages 3–5 (defined as eGFR <60 mL/min/1.73 m2).
2.4. Genotyping of Polymorphisms
DNA was extracted from buffy coat with a BioRobot M48 Workstation (QIAGEN Group, Tokyo). The genotyping of PPARD A65G in exon 7 (Asn163Asn) (rs2076167), PPARD T-48444C in exon 3 (rs6902123), PPARD T-842C in exon 3 (rs2267668), PPARG C161T (His477His) (rs3856806), PPARG Pro12Ala (rs1801282), PPARGC1A Thr394Thr (rs2970847), and PPARGC1A Gly482Ser (rs8192678) polymorphisms was conducted by the RIKEN institute using multiplex polymerase chain reaction-based Invader assay (Third Wave Technologies, Madison, WI, USA) as described previously [20]. The genotype distributions of all the 108 polymorphisms examined in this cross-sectional study are shown in the recently published data [18].
2.5. Statistical Analysis
Logistic regression analysis was performed for estimating age- and sex-adjusted odds ratios (aORs) and 95% confidence intervals (CIs) for CKD by genotype. All the other potential confounders of BMI, systolic blood pressure, diastolic blood pressure, use of anti-hypertensive drugs, fasting plasma glucose, HbA1c, use of glucose-lowering drugs, total cholesterol, HDL cholesterol, triglyceride, use of lipid-lowering drugs, uric acid, past history of cardiovascular diseases, past history of cerebrovascular diseases, smoking status, and drinking status were tested for change in estimate (CIE) [21] to see if any of these covariates produces significant change in estimates. We decided not to include any of these variables because none of them fulfilled the criteria of CIE 0.1. As it is known that metabolic factors such as HbA1c and presence of DM will affect CKD risk, they are considered to be causal intermediate of (PPARD) SNP and CKD risk. Adjusting for even a partially causal intermediate phenotype will incorrectly remove a true association and could potentially bias the true association [22]; thus we decided not to adjust for these variables. Gene-environment interactions were assessed by the logistic model, which included a multiplicative interaction term as well as variables for each genotype, age, sex, and smoking and drinking habits. Age adjustments in the analyses were done with ages regarded as continuous variables. Trend analyses by genotypes were done with genotypes for each polymorphism coded as ordinal-categorical variables. Differences in the distribution of the values of each characteristic variable between those with CKD and those without were examined by Student’s -test or by the test. Accordance with the Hardy-Weinberg’s equilibrium, which indicates an absence of discrepancy between genotype and allele frequencies, was checked using the test. Haplotype analysis using genotypes in two loci was conducted by the “haplologit” command of STATA adjusted for age and sex based on the EM algorithm [23]. The linkage disequilibrium (LD) between the polymorphisms in two loci ( and ) was estimated by the “pwld” command of STATA. All the calculations were done using the STATA version 10 (Stata Corp, College Station, TX, USA).
3. Results
3.1. Characteristics of the Subjects and Allele Frequency of the PPARD, PPARG, and PPARGC1A Polymorphisms
The characteristics of the subjects are summarized by CKD status in Table 1. The mean age ± standard deviation was 56.7 ± 8.6 years, and the males were 48.7% of the whole number of subjects. Subjects with CKD accounted for 17.3% (575/3,315) of the entire study population.
The genotype frequencies among the genotyped subjects included in the analyses were in accordance with Hardy-Weinberg’s equilibrium for all of the PPARD A65G ( allele = 0.221 (minor allele frequency), , and ), PPARD T-48444C ( allele = 0.020, , ), PPARD T-842C ( allele = 0.201, , ), PPARG C161T ( allele = 0.152, , ), PPARG Pro12Ala (=C/G; allele = 0.031, , ), PPARGC1A Thr394Thr (=C/T; allele = 0.218, , ), and PPARGC1A Gly482Ser (=G/A; allele = 0.462, , ). The allele frequencies were similar to those among all the genotyped 4,519 subjects: 0.215 for PPARD 65G, 0.020 for PPARD , 0.195 for PPARD , 0.152 for PPARG 161, 0.031 for PPARG 12Ala (), 0.219 for PPARGC1A Thr394Thr , and 0.460 for PPARGC1A 482Ser () [18]. Genotype call rate was more than 99.8% for each genotype among those with SCr data ().
3.2. PPARD, PPARG, and PPARGC1A Polymorphisms and Risk of CKD
When those with PPARD T-842C (rs2267668) were defined as reference, those with PPARD T-842C and demonstrated the OR for CKD of 1.26 (95%CI 1.04–1.53) and 1.31 (95%CI 0.83–2.06), respectively, with the significant trend for increased OR with the increasing number of allele (). There were no significant associations between the polymorphisms in other polymorphisms in genes encoding PPARs, PPARD A65G in exon 7 (Asn163Asn) (rs2076167), PPARD T-48444C in exon 3 (rs6902123), PPARG C161T (His477His) (rs3856806), PPARG Pro12Ala (rs1801282), PPARGC1A Thr394Thr (rs2970847), and PPARGC1A Gly482Ser (rs8192678), with the risk of CKD. When those with at least one minor allele of each PPAR SNP were combined together and compared with the references (subjects with homozygous major allele) (dominant model), those with PPARD T-842C demonstrated the significantly increased risk of CKD with the aORs of 1.27 (1.05–1.53), while no other SNP turned out to be significant in this analysis (Table 2). The eGFRs by GCK and GCKR genotypes were substantially in the same trend to this analysis, with the marginally significant decreasing effect of PPARD T-842C on eGFR, suggesting the influence of GCK and GCKR polymorphisms on human renal functions (Table 3).
Because one single nucleotide polymorphism (SNP) in PPARD gene (rs2267668) was significantly correlated with the CKD risk, we also estimated the LD within PPARD polymorphisms, which revealed that the 3 polymorphisms in PPARD gene investigated, PPARD A65G in exon 7 (Asn163Asn) (rs2076167), PPARD T-48444C in exon 3 (rs6902123), and PPARD T-842C in exon 3 (rs2267668), were in linkage disequilibrium to each other (Figure 1). Haplotype analysis of the PPARD polymorphisms on CKD risk did not reveal any significant association of PPARD haplotypes with the risk of CKD.
We also tested possible interaction between PPAR genotypes and lifestyle/etiologic factors including smoking, alcohol drinking, hypertension, dysglycemia, dyslipidemia, and high uric acid, none of which resulted in statistical significance (data not shown). As this study was held in 10 institutions, of which 8 had data for eGFR, we also conducted the analyses adjusted for institutions, the results of which were not substantially different from the unadjusted ones.
4. Discussion
To date, multiple factors have been identified to play key roles in the genesis of CKD induced by metabolic disorders, among which it is well recognized that insulin resistance and hyperinsulinemia are important factors [23]. In the pathogenesis of diabetic nephropathy patients, it is reported that renal structural features precede the appearance of overt diabetes, suggesting that renal injury is initiated in the phase of prediabetes [6]. In this process, the effect of insulin to stimulate both the renin-angiotensin system and nitric oxide in the renal vasculature is considered to be important [6]. PPARD plays an important role in energy homeostasis and overexpression of PPARD in skeletal muscle, and adipocytes are shown to increase fat catabolism in both of these tissues [24]. Given the important roles of FFA in the genesis of type 2 diabetes, the enhanced fatty acid oxidation as well as the upregulation of the oxidative phosphorylation pathway by PPARD may be beneficial in alleviation of insulin resistance and adiposity [24]. The fact that the PPARD agonists decrease insulin and glucose levels by increasing glucose transport suggests that genetic variations in PPARD gene will modify the glucose metabolism and possibly affect subsequent CKD risks. In the recent study, the rare single nucleotide polymorphisms (SNPs) of PPARD in LD with rs2267668 SNP and the PPARGC1A SNP were shown to have additive effects on the risk of type 2 DM, suggesting the etiological roles of this SNP also in the genesis of CKD [25]. The present study suggested that the polymorphism in the PPARD gene, a polymorphism in the gene encoding the key molecule controlling the expression of genes involved in fatty acid oxidation and energy uncoupling in skeletal muscles [26], is significantly associated with the risk of CKD in Japanese. As far as we know, this is the second report that investigated whether the genetic variations in PPAR genes may influence the renal functions, that is, the risk of CKD in Japanese [27], and the first one that indicated that the polymorphism in the PPARD gene may influence the risk of CKD in humans. It is also reported that common polymorphisms in PPARD gene were associated with fasting glucose, insulin resistance, and risk of conversion from impaired glucose tolerance to type 2 diabetes in Europeans [28], and some functional variants in PPARD gene were reportedly associated with fasting glucose and BMI in Koreans [29] and fasting glucose levels and insulin sensitivity in Chinese [30, 31]. Our study results suggested the trend that the subjects with the allele of PPARD T-842C polymorphism were at an increased risk of CKD, which was in line with our hypothesis considering that the PPARD T-842C polymorphism found to be significantly associated with the risk of CKD in the present study is in tight linkage disequilibrium with other functional polymorphisms in PPARD gene, PPARD T294C (rs2016520) ( and in Japanese [JPT]) [32]. This functional polymorphism is located in the untranslated exon 4 of the PPARD gene influencing the expression of PPARD through altering the sequence of the DNA-binding site for Sp-1 [33, 34], and that demonstrated the possibility that the compromised function of PPARD due to this polymorphism might increase the risk of CKD as well as other vascular diseases. The allele of the PPARD T-842C polymorphism is tightly linked to the allele of the PPARD T294C polymorphism, carriers of which are shown to have a higher plasma LDL-cholesterol levels and it reportedly has a slightly higher promoter activity than the allele in one in vitro study, which may promote lipid accumulation in macrophages as well as lipid uptake and storage [35]. The present study demonstrated that the subjects with the allele of PPARD T-842C polymorphism were at an increased risk of CKD, which is considered to be biologically relevant, given the upregulation of lipid uptake/storage in the allele carrier of the PPARD T294C polymorphism. However, considering that no evident associations between these PPARD SNPs and metabolic risk factors were observed in the present study, there seems to be a possibility that other factors such as inflammation are also involved in this CKD risk modification by the PPARD SNPs [36, 37].
Associations between polymorphisms in PPARG and risk of CKD including diabetic nephropathy (DN) are also reported; in a meta-analysis of 18 studies, it is demonstrated that the Ala allele carriers of the PPARG Pro12Ala polymorphism [38] had the reduced risk of DN, although this association was not significant when restricted only to Asians. Another study in Japanese showed that PPARG C-681G (rs10865710) polymorphism was significantly associated with the risk of CKD among subjects with hypertension, suggesting the biological roles of PPARs in the genesis of CKD also in Asian populations [27]. Regarding PPARGC1A polymorphisms, there is only one study that examined its association with DN in Asian Indians, which found that the Gly482Ser polymorphism of the PPARGC1A was significantly associated with DN [39]. Meanwhile, there are several studies that examined the association of these polymorphisms with the risk of type 2 diabetes (T2DM); the meta-analysis of the reported study results revealed that the PPARGC1A Thr394Thr polymorphism was not significantly associated with the T2DM risk in East Asian populations [40]. Our study results revealed no significant association of CKD risk with the polymorphisms in these genes encoding PPARs, PPARG C161T (His477His), PPARG Pro12Ala, PPARGC1A Thr394Thr, and PPARGC1A Gly482Ser, suggesting that the roles of these genetic variations are limited in Japanese, or the sample sizes of the present study were relatively underpowered to detect these associations.
A recent report from the STOP-NIDDM trial showed that the rare polymorphism of PPARD gene in linkage disequilibrium with the rs2267668 polymorphism was studied here and the 482Ser allele of the PPARGC1A rs8192678 polymorphism had the additive effects on the risk of type 2 diabetes [22], which makes the combined effects of these polymorphisms on the risk of CKD also an issue of interest. Although the present study results were considered to be in a similar trend for the effect of each of these polymorphisms on CKD risk, we could not detect any statistical interaction between these polymorphisms on the risk of CKD, which necessitates the investigation of this association with larger sample sizes in the near future.
There are some potential limitations in this study. All of the CKD cases are diagnosed based on the SCr data, which might potentially be different from the actual GFR based on the renal measurement, and CKD diagnosis itself may also be different from that defined by the recent KDIGO 2012 Guideline, including the presence of albuminuria [41], which might have diluted the effect of each genotype on CKD risk. Different etiologies causing CKD should also have been considered in the analyses, but it could not be examined in the present study because of the limited availability of the clinical information including the pathologic one. Adjustments for multiple comparisons may be another issue. There are also a number of criticisms suggesting that correction of multiple comparisons by Bonferroni procedures is sometimes too conservative and aggravates the researchers’ tendency not to conduct or present more tests [42], and considering that the present study was conducted under rather exploratory context, especially for the gene-environment interaction part, we decided not to adopt these adjustments in this study. Further investigations with better study designs will also be required.
In summary, the present study found a significantly increased risk of CKD in those with the C allele of PPARD T-842C, which may suggest the future possibility of personalized risk estimation of this life-limiting disease in the near future.
Conflict of Interests
The authors have no conflict of interests to declare.
Acknowledgments
The authors wish to thank Mr. Kyota Ashikawa and Ms. Tomomi Aoi at the Laboratory of Genotyping Development, Center for Genomic Medicine, RIKEN, for genotyping. The authors also thank Ms. Yoko Mitsuda and Ms. Keiko Shibata at Nagoya University for their technical assistance. This study was supported in part by a Grant-in-Aid for Scientific Research on Priority Areas of Cancer (no. 17015018) and Scientific Support Programs for Cancer Research, Grant-in-Aid for Scientific Research on Innovative Areas (no. 221S0001) from the Japanese Ministry of Education, Culture, Sports, Science and Technology. The authors were with the J-MICC Study Group.