Oxidative Medicine and Cellular Longevity
Volume 2016 (2016), Article ID 8516748, 7 pages
http://dx.doi.org/10.1155/2016/8516748

Research Article

Interplay between Superoxide Dismutase, Glutathione Peroxidase, and Peroxisome Proliferator Activated Receptor Gamma Polymorphisms on the Risk of End-Stage Renal Disease among Han Chinese Patients

Chia-Ter Chao,1,2 Yen-Ching Chen,³ Chih-Kang Chiang,2,4 Jenq-Wen Huang,⁵ Cheng-Chung Fang,⁶ Chen-Chih Chang,⁵ and Chung-Jen Yen5,7

¹Department of Medicine, National Taiwan University Hospital Jinshan Branch, New Taipei City, Taiwan
²Graduate Institute of Toxicology, College of Medicine, National Taiwan University, No. 1 Jen-Ai Road, Section 1, 10051 Taipei, Taiwan
³Institute of Epidemiology and Preventive Medicine, College of Public Health, National Taiwan University, No. 17 Xu-Zhou Road, 10055 Taipei, Taiwan
⁴Department of Integrated Diagnostics & Therapeutics, National Taiwan University Hospital, No. 7 Chung-Shan South Road, 10002 Taipei, Taiwan
⁵Department of Internal Medicine, National Taiwan University Hospital, No. 7 Chung-Shan South Road, 10002 Taipei, Taiwan
⁶Department of Emergency Medicine, National Taiwan University Hospital, No. 7 Chung-Shan South Road, 10002 Taipei, Taiwan
⁷Department of Geriatric Medicine and Gerontology, National Taiwan University Hospital, No. 7 Chung-Shan South Road, 10002 Taipei, Taiwan

Received 5 September 2015; Accepted 25 November 2015

Academic Editor: Ada Popolo

Copyright © 2016 Chia-Ter Chao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Background. Single nucleotide polymorphisms (SNPs) of antioxidants, including superoxide dismutase 2 (SOD2) and glutathione peroxidase 1 (GPX1), play an important role in the risk for cancer and metabolic disorders. However, little is known regarding the effect of antioxidant SNPs on renal events. Methods. We prospectively enrolled multicenter patients with end-stage renal disease (ESRD) and those without chronic kidney disease (CKD) of Han Chinese origin, with SOD2 (Val16Ala), GPX1 (Pro197Leu), and PPAR-γ (Pro12Ala, C161T) genotyped. Multiple regression analyses were conducted to evaluate the significant risk determinants for ESRD. Results. Compared to ESRD patients, non-CKD subjects were more likely to have T allele at SOD2 Val16Ala () and CC genotype at PPAR-γ Pro12Ala (). Regression analysis showed that TT genotype of SOD2 Val16Ala conferred significantly lower ESRD risk among patients without diabetes (odds ratio 0.699; ). GPX1 SNP alone did not alter the risk. We detected significant interactions between SNPs including PPAR-γ Pro12Ala, C161T, and GPX1 regarding the risk of ESRD. Conclusion. This is the first and largest study on the association between adverse renal outcomes and antioxidant SNPs among Han Chinese population. Determination of SOD2 and PPAR-γ SNPs status might assist in ESRD risk estimation.

1. Introduction

Patients with persistent renal failure are increasing worldwide, and the number is expected to rise even higher as the population age with multimorbidity. This creates tremendous burden on the existent healthcare resources. In USA, the annual expenditure for the care of end-stage renal disease (ESRD) patients reaches nearly 50 billion dollars in 2011, with an incremental range of 2.8 to 3.4% per year [1]. Given this unavoidable trend, researches focusing on the risk factors and mortality determinants for ESRD assume increasing importance.

Oxidative stress is proposed to play an important role in the pathogenesis of diseases including cardiometabolic diseases, neurodegenerative diseases, neoplasm formation, and, most important of all, chronic kidney disease (CKD) [2–4]. Surge of oxidative stress is frequently linked to reactive oxygen species (ROS) formation and/or lowering antioxidant capacity, incurring damage to the surrounding tissues and organ injury. Defending against oxidative stress, antioxidants could be categorized as exogenous (carotene, tocopherol, ascorbate, etc.) and endogenous, especially thiol group-containing macromolecules and enzymes [5]. There are reports suggesting that prooxidant levels increase and antioxidant enzyme activities decrease stepwise as renal function deteriorates [6]. Moreover, ROS accumulation could not only perpetuate immunogenic glomerular/interstitial injury, but also contribute to endothelial dysfunction over the long run, further deranging the intraglomerular and systemic hemodynamics [7].

Among the endogenous antioxidants, superoxide dismutase (SOD) and glutathione peroxidase (GPX) are particularly important players. Several studies have investigated the risk of CKD conferred by the genetic heterogeneity of these two genes. Both SOD1 and SOD2 single nucleotide polymorphisms (SNPs) are found to be associated with 2- to 3-fold higher risk of developing diabetic nephropathy in type 1 and 2 diabetes mellitus (DM) patients [8, 9]. CKD patients carrying SOD2 Ala16Val genotype also had greater decline of estimated glomerular filtration rate (eGFR) during follow-up [10]. Similarly, plasma GPX activity also correlated with annual eGFR change in CKD patients [10]. However, very few studies focused on Asian population, and none of the past reports addressed the influence of SNPs of these antioxidants, especially SOD and GPX, on the risk of dialysis-dependent ESRD patients. Previously, we have identified that peroxisome proliferator activated receptor- (PPAR-) γ SNPs, Pro12Ala, and C161T significantly modified the risk of developing ESRD and prognosis of ESRD patients of Han origin [11]. In the current study, we aimed to determine the relationship between SOD2, GPX1, PPAR-γ genotypes, and the risk of developing ESRD among patients of Han origin, using a large multicenter cohort.

2. Material and Methods

2.1. Study Population

In this multicenter study, we identified ESRD patients () who initiated long-term hemodialysis or peritoneal dialysis from 3 hospitals and 9 clinics in the northern Taiwan between 2002 and 2003 [12–15]. Another 780 patients without chronic kidney disease (CKD) were recruited from the health checkup services during the same time period. Blood samples were collected from each subject at the time of enrollment, and individuals were prospectively followed up. The current study was approved by the institutional review boards of all participating institutions.

All subjects completed a questionnaire about their sociodemographic status (e.g., birthday, sex, and education), comorbidity history (e.g., hypertension and diabetes), and lifestyle factors (e.g., smoking and alcohol consumption) with the assistance of general practitioners. Diabetes was determined based on a history of physician-diagnosed diabetes or 2 separate episodes of fasting blood glucose levels ≥126 mg/dL (7.0 mmole/L). The causes of ESRD and the vintage were recorded for all ESRD cases at enrollment. Diabetic nephropathy as the cause of ESRD was recognized if the patient had a long-standing diabetes history (at least 5 years) before dialysis start, presence of diabetic retinopathy, and/or proteinuria (≥100 mg/dL by urinalysis) without other causes.

2.2. Laboratory Assays (Genotyping)

Genomic DNA was extracted from peripheral blood leukocytes using the Chemagic DNA Blood kits (Chemagen AG, Baesweiler, Germany). Genotypes of the SNPs were determined by the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP). Fragments encompassing the SNPs (V16A of SOD2, rs1050450 of GPX1, P12A, and C161T of PPAR-γ) were generated from genomic DNA by PCR using the primer pairs shown in Table 1. The PCR products were then digested with restriction enzymes (New England BioLabs, Beverly, MA, USA), electrophoresed on a 2.5% agarose gel, and stained with ethidium bromide. The size of the digested PCR products for each SNP variant was also summarized in Table 1.

Table 1: Primers and restriction enzymes for PCR-RFLP.

Gene/SNP
Primers for PCR
Restriction enzyme
Size of digestion products (base pairs)

SOD2/V16A
5′-CTGACCGGGCTGTGCTTTCT-3′ 
BsaWI
TT: ~50, 175 

5′-CAACGCCTCCTGGTACTTCT-3′
CT: ~50, 175, 225 

CC: ~225

GPX1/P197L
5′-GCTTCCAGACCATTGACATC-3′ 
ApaI
CC: ~53, 261 

5′-TCCCAAATGACAATGACACAG-3′
CT: ~53, 261, 314 

TT: ~314

PPAR-γ/P12A
5′-GCCAATTCAAGCCCAGTC-3′ 
BstUI
CC: ~270 

5′-GATATGTTTGCAGACAGTGTATC 
CG: ~43, 227, 270 

AGTGAAGGAATCGCTTTCCG-3′
GG: ~43, 227

PPAR-γ/C161T
5′-TTTGACTGAACCCCCTGTTG-3′ 
NlaIII
CC: ~101, 245 

5′-CAGAATAGTGCAACTGGAAGA-3′
CT: ~41, 101, 204, 245 

TT: ~41, 101, 204

2.3. Statistical Analysis

Data analysis was performed using the SAS, 9.1.3 software (SAS Institute Inc., Cary, NC, USA) and the R, 2.11.1 software (R Foundation for Statistical Computing, Vienna, Austria). In statistical testing, two-sided value ≤ 0.05 was considered statistically significant. The distributional properties of continuous variables were expressed as mean ± standard deviation (SD), whereas categorical variables were presented as frequency and percentage.

The potential risk factors for ESRD were examined in univariate analysis using chi-square test, Fisher’s exact test, two-sample -test, Wilcoxon rank-sum test, and log-rank test, as appropriate. Multivariate analysis was subsequently conducted using weighted logistic regression model to estimate the effects of risk factors on the probability of long-term dialysis.

To ensure the quality of analysis results, basic model-fitting techniques for (1) variable selection, (2) goodness-of-fit (GOF) assessment, and (3) regression diagnostics were used in our regression analyses. Specifically, the stepwise variable selection procedure (with iterations between the forward and backward steps) was applied to obtain the candidate final regression model. All the univariate significant and nonsignificant relevant covariates listed in Table 3 and their interaction variables were selected and the significance levels for entry (SLE) and for stay (SLS) were set to 0.15 or larger (up to 0.35). The best final regression model was identified by reducing the significance levels to 0.05. Both the GOF measures, including the percentage of concordant pairs, estimated area under the receiver operating characteristic (ROC) curve, and adjusted generalized , and the GOF tests, including the deviance GOF test, Pearson chi-squared GOF test, and Hosmer-Lemeshow GOF test (for logistic regression model), and, if feasible, Grønnesby-Borgan GOF test (for Cox’s proportional hazards model), were examined to assess the GOF of the fitted regression models. The generalized additive models (GAM) were applied to detect nonlinear effects of continuous covariates.

3. Results

3.1. Clinical Features of Enrollees

The recruited numbers of ESRD patients and their non-CKD counterparts are comparable. Comparisons of the minor allele frequencies of the four SNPs and examinations of the Hardy-Weinberg equilibrium (HWE) among ESRD patients and non-CKD ones were shown in Table 2. Both SOD2 and GPX1 SNPs satisfied the HWE in non-CKD and ESRD patients, while PPAR-γ Pro12Ala among ESRD cases showed deviation from HWE (). This suggests a potential association between PPAR-γ Pro12Ala and the presence of ESRD, as demonstrated in our previous study [11].

Table 2: Comparisons of the minor allele frequencies of the four SNPs and examinations of the Hardy-Weinberg equilibrium among the entire cohort.

SNPs
rs #
Nucleotide Change (amino acid change)
Location
Controls ( = 780)
Cases ( = 671)

MAF
HWE value
MAF
HWE value

SOD2
rs4880
T → C 
(Val → Ala)
exon 2
13.97%
0.768
17.36%
0.423

GPX1
rs1050450
C → T 
(Pro → Leu)
exon 2
5.38%
0.720
5.29%
0.249

P12A
rs1801282
C → G 
(Pro → Ala)
exon B
3.08%
0.526
4.99%
0.020

C161T
rs3856806
C → T 
(His → His)
exon 6
24.29%
0.697
26.53%
0.488

HWE = Hardy-Weinberg equilibrium.

MAF = minor allele frequency.

Table 3: Characteristics of patients with ESRD and those without CKD.

Variables
Non-CKD
ESRD patients

Number of subjects
780
671

Age (years)
60.0 ± 19.2
58.9 ± 14.6
<0.001

Male
410 (52.6%)
327 (48.7%)
0.155

DM
83 (10.6%)
253 (37.7%)
<0.001

Serum creatinine (mg/dL)
1.0 ± 0.3
11.2 ± 2.5
<0.001

SOD2 exon 2 genotype

0.036

    TT
576 (73.85%)
455 (67.8%)

    CT
190 (24.36%)
199 (29.66%)

    CC
14 (1.79%)
17 (2.53%)

GPX1 exon 2 genotype

1.000

    CC
697 (89.36%)
600 (89.42%)

    CT
82 (10.51%)
71 (10.58%)

    TT
1 (0.13%)
0 (0.00%)

PPAR- exon B genotype

0.028

    CC
733 (93.90%)
609 (90.99%)

    CG
46 (5.97%)
57 (8.30%)

    GG
1 (0.13%)
5 (0.72%)

PPAR- exon 6 genotype

0.382

    CC
449 (57.69%)
366 (55.08%)

    CT
283 (36.21%)
254 (37.63%)

    TT
48 (6.10%)
51 (7.30%)

Continuous variables were tested by Kruskal-Wallis test, whereas categorical variables were tested by Fisher’s exact test.

The genotype distributions and the other clinical characteristics were compared between the enrolled ESRD and non-CKD patients in Table 3. The mean serum creatinine (mg/dL) was significantly higher in ESRD patients (non-CKD versus ESRD, versus ; ). In addition, the proportion of patients with DM was significantly higher in ESRD patients (non-CKD versus ESRD, 10.6% versus 37.7%; ). The distributions of the SOD2 exon 2 and PPAR-γ exon B genotypes were significantly different between the two groups. TT genotype of SOD2 exon 2 was predominant in non-CKD patients (), while CC genotype of PPAR-γ exon B occurred more frequently in non-CKD patients (). There was no significant heterogeneity regarding GPX1 and PPAR-γ exon 6 genotypes between ESRD and non-CKD patients.

We further stratified the ESRD patients according to their ESRD etiology as DM or non-DM related, since PPAR-γ polymorphism is reportedly associated with development of DM nephropathy (DMN) (Table 4) [16]. About 27% patients had DMN-related ESRD. No significant difference in genotype distributions of SOD2, GPX1, and PPAR-γ was noted between those with or without DMN.

Table 4: Characteristics of ESRD patients due to DMN and non-DMN related ESRD.

Variables
Non-DMN-related ESRD patients
DMN-related ESRD patients

Number of subjects
492
179

SOD2 exon 2 genotype

0.716

    TT
330 (67.1%)
125 (69.8%)

    CT
150 (30.5%)
49 (27.4%)

    CC
12 (2.4%)
5 (2.8%)

GPX1 exon 2 genotype

0.258

    CC
444 (90.2%)
156 (87.2%)

    TC
48 (9.8%)
23 (12.8%)

    TT
0 (0.00%)
0 (0.00%)

PPAR- exon B genotype

0.284

    CC
451 (91.7%)
158 (88.3%)

    CG
37 (7.5%)
20 (11.2%)

    GG
4 (0.8%)
1 (0.6%)

PPAR- exon 6 genotype

0.454

    CC
263 (53.5%)
103 (57.5%)

    CT
193 (39.2%)
61 (34.1%)

    TT
36 (7.3%)
51 (8.4%)

Continuous variables were tested by Kruskal-Wallis test, whereas categorical variables were tested by Fisher’s exact test.

3.2. Risk of ESRD Associated with SNPs of SOD2 and GPX1

We first investigated the influence of SOD2 and GPX1 genetic polymorphism on developing ESRD. As shown in Table 5, multiple logistic regression analysis incorporating factors in Table 3 identified five significant predictors of ESRD developing among the entire cohort. Patients with DM and SOD2 exon 2 CC genotype had significantly higher risk of developing ESRD. Comorbid DM played an important role in determining the influence on ESRD risk brought by SOD2 and GPX1 genotypes. It is interesting to observe that patients without DM but carrying SOD2 exon 2 TT genotypes had lower chances of developing ESRD (OR 0.699, 95% confidence interval [CI] 0.52–0.94; ).

Table 5: Multivariate analysis of the predictors of dialysis by fitting multiple logistic regression model with the stepwise variable selection method, including SOD2 and GPX1 polymorphisms#.

Covariates
Regression coefficient
Standard error
value
value
Odds ratio
95% confidence interval

Intercept

<0.001
—
—
—

Age (between 28 and 65)

<0.001
4.847
3.801
6.180

Male

0.040
0.777
0.611
0.988

DM

<0.001
5.792
4.035
8.316

Non-DM × SOD2 exon 2 TT

0.018
0.699
0.520
0.939

Multiple logistic regression model: , adjusted generalized , estimated area under the receiver operating characteristic (ROC) curve = 0.762, and the Hosmer and Lemeshow goodness-of-fit chi-square test (df = 8).

#Variables with odds ratio (OR) of extreme values were listed below: DM × SOD2 exon 2 CC, OR 3.6 × 10⁵ ().

3.3. Risk of ESRD Associated with SNPs of SOD2, GPX1, PPAR-γ Pro12Ala, and C161T

We further combined the previously reported PPAR-γ SNPs (i.e., Pro12Ala and C161T) and SOD2 as well as GPX1 SNPs in statistical models for a more in-depth analysis. Regression analysis identified multiple significant predictors of developing ESRD (Table 6). Individuals without DM but with SOD2 exon 2 TT genotype were still at lower risk of ESRD (OR 0.69, 95% CI 0.51–0.93; ). Additionally, we discovered that PPAR-γ SNPs and GPX1 SNPs interacted with regard to their effect on ESRD risk. Patients with PPAR-γ exon 6 TT genotype and exon B GG genotype displayed a high risk of ESRD (), while exon 6 CC genotype and exon B CC genotype were protective against ESRD (OR 0.78, 95% CI 0.61–0.99; ). Patients with simultaneous GPX1 exon 2 CC genotype and PPAR-γ exon B GG genotype had significantly lower ESRD risk ().

Table 6: Multivariate analysis of the predictors of dialysis by fitting multiple logistic regression model with the stepwise variable selection method, including SOD2, GPX1, and PPAR- polymorphisms#.

Covariates
Regression coefficient
Standard error
value
value
Odds ratio
95% confidence interval

Intercept

0.162

<0.001
—
—
—

Age (between 28 and 65)

0.125

<0.001
4.879
3.820
6.232

Male

0.123

0.051
0.786
0.618
1.001

DM

0.186

<0.001
5.957
4.141
8.569

Non-DM × SOD2 exon 2 TT

0.151

0.014
0.691
0.514
0.928

PPAR-γ exon 6 CC × PPAR-γ exon B CC

0.123

0.040
0.778
0.609
0.988

Multiple logistic regression model: = 1451, adjusted generalized = 0.304, estimated area under the receiver operating characteristic (ROC) curve = 0.769, and the Hosmer and Lemeshow goodness-of-fit chi-square test = 0.0047 (df = 8).

#Variables with odds ratio (OR) of extreme values were listed below: DM × SOD2 exon 2 CC, OR 2.7 × 10⁶ ( < 0.001); non-DM × PPAR-γ exon B GG, OR 3.8 × 10¹³ ( < 0.001); PPAR-γ exon 6 TT × PPAR-γ exon B GG, OR 8.2 × 10¹¹ ( < 0.001); GPX1 exon 2 CC × PPAR-γ exon B GG, OR 1.9 × 10⁻⁷ ( < 0.001).

4. Discussion

In the current study, we addressed the issue whether antioxidative enzymes (SOD2 and GPX1) SNPs affect the risk of developing ESRD, in a large cohort of patients of Han origin, and their potential interactions with other known ESRD risk-modifying SNPs (PPAR-γ). We discovered that risk conferred by SOD2 and PPAR-γ SNPs in Han Chinese population is modified in the presence of DM, and interactions exist among different SNPs. Particularly, PPAR-γ exon B GG genotypes increased risk of ESRD in patients without DM and those with exon 6 TT genotype. In addition, the influence of SOD2 SNPs on the risk of ESRD differed between patients with and without DM. These findings suggested a renewed explanation for the associations between both antioxidants, PPAR-γ SNPs and DM, non-DM related ESRD in the literature.

The SOD2 SNPs, rs4880, substitute a C to T at position 2734 and change amino acid coding from alanine to valine at position 16 (Ala16Val). This alteration leads to mRNA instability and reduced SOD2 expression. SOD2 Ala16Val SNPs have been found to correlate variably with different cancer susceptibility. A meta-analysis concluded that Ala allele conferred a borderline elevated risk of breast cancer [17]. Other researchers identified a higher risk of prostate cancer (odds ratio, 1.16) and increased aggressiveness among Ala alleles of SOD2 exon 2 SNP carriers [18]. Similarly, carriers of SOD2 Ala16Val exhibited a higher risk of bladder cancer, particularly among concurrent smokers [19]. On the other hand, SOD2 Ala16Val SNPs might also play a role in the susceptibility of multiple metabolic disorders. Ala allele has been linked with elevated risk of diabetic macular edema and diabetic retinopathy in Asian population [20, 21]. Moreover, SOD2 Ala16Val SNPs also serve as an important risk factor for atherosclerotic burden, cardiac events, and nonalcoholic steatohepatitis [22–24]. These vascular and metabolic abnormalities among SOD2 Ala16Val carriers, coupling with the susceptibility to urologic malignancies, might pave the way for incident CKD development, progression, and ESRD [25]. Indeed, subsequent study revealed that Val allele of SOD2 exon 2 carriers had a 3-4-fold higher risk of rapid renal function decline compared to the others, but this was established only in Caucasians [10]. Based upon a large cohort of Han Chinese patients, we discovered that CC genotype of SOD2 exon 2 portended significantly higher risk of ESRD among those with DM, while TT genotype decreased the risk among non-DM ones (Tables 5 and 6). Although our identified at-risk allele (Ala) in Han Chinese population differs from that of Caucasians (Val) in most studies, both findings underline the importance of SOD2 exon 2 SNPs in modifying the risk of CKD progression as well as ESRD.

GPX1, a selenoenzyme that catalyzes the breakdown of various peroxides into O₂ and H₂O, is produced by nearly all tissues and serves as a cellular barrier against oxidative stress. Past studies focusing on the role of GPX1 SNPs mostly included GPX1 SNP rs1050450, which causes a C to T mutation [18]. GPX1 SNP rs1050450 is found to influence the risk of lung cancer and bladder cancer [26]. Risk of vascular calcification and atherosclerosis is also affected by Leu allele of GPX1 SNPs [24, 27]. However, available evidence did not agree that GPX1 SNP rs1050450 plays a significant role in CKD progression or renal allograft dysfunction [10, 28]. We similarly did not detect significant influences of GPX1 Pro197Leu SNPs on the risk of ESRD in Han Chinese population (Tables 5 and 6). This is compatible with the literature and further supports the finding that GPX1 SNPs alone did not affect the risk of ESRD in most reports. However, a subgroup of patients with GPX1 CC genotype did demonstrate lower risk of ESRD, that is, those with concurrent PPAR-γ exon B GG (Table 6). A previously undetected interaction between antioxidants SNPs and PPAR-γ SNPs might exist and warrants further investigation.

The interactions between DM/antioxidant SNPs and within different antioxidant SNPs on the risk of ESRD are interesting. From our findings, we propose that DM and CC genotype might act synergistically to elevate the risk of subsequent ESRD or that CC genotype might promote the development of diabetic nephropathy (Table 5). However, past studies failed to disclose any significant effect of SOD2 Ala16Val SNPs on the presence of DM, but, on the contrary, TT genotype carrier of SOD2 exon 2 was more likely to be found in Caucasian and Japanese patients with DM nephropathy than those without [8, 20, 29]. It is likely that SOD2 exon 2 SNPs, specifically CC genotypes, contribute to the risk of DM nephropathy among Han Chinese population. Furthermore, we also discovered that patients without DM were at significantly lower risk of ESRD if they carried TT genotypes of SOD2 exon 2 (Tables 5 and 6). This could be a circumstantial evidence supporting the detrimental effect of CC genotypes on ESRD, while it alternatively could suggest a protective effect brought by TT genotypes against non-DM ESRD. Finally, our findings that PPAR-γ Pro12Ala interacted with non-DM status on the effect of risk for ESRD could also be a token of gene-environmental interaction, as PPAR-γ Pro12Ala has been found to act together with body mass index or dietary fatty acid intake to increase cardiometabolic risk [30, 31].

Several lines of evidence suggested that interactions between different SNPs could modify the risk of CKD and its progression. FIND (Family Investigation of Nephropathy and Diabetes) study group found that α-actinin-4 polymorphism could interact with APOL1 SNPs with regard to the risk of non-DM related ESRD in African Americans [32]. The same group subsequently utilized a large cohort study to identify podocin (NPHS2) and APOL1 SNPs might jointly promote nondiabetic nephropathy in African Americans, although the mechanisms of synergism still remained undetermined [33]. Importantly, PPAR-γ activation has been reported to regulate renal mesangial inflammation and modulate renal fibrosis through p65 and NF-κB pathways [34]. Furthermore, experimental evidence also suggested that PPAR-γ participates in regulating oxidative stress and acute kidney injury [35]. Consequently, we could expect that PPAR-γ SNPs and even antioxidant SNPs might interact with regard to their influence on the risk of ESRD. In this study, we discovered that PPAR-γ SNPs, Pro12Ala, C161T, and GPX1 SNPs might interact with each other, lending support to our hypothesis that there could be interplay between different SNPs concerning their clinical effects.

Our study has its strength and limitations. The case number was large, and our investigation on the reciprocal relationship between different SNPs among a population receiving less attention, Han Chinese, made this study more appealing. However, the current study was also limited by its clinical setting and the limited inclusion of antioxidant SNPs. Further research utilizing more comprehensive methods for SNPs survey might be needed in the future.

5. Conclusion

Using a large Han Chinese-based cohort, we studied the modifying effect of SOD2, GPX1, and PPAR-γ SNPs on the risk of ESRD. We found that CC genotype of SOD2 elevated the risk of ESRD in patients with DM, while TT genotype of SOD2 decreased the risk in patients without DM. In addition, interactions between PPAR-γ exon B/exon 6 and PPAR-γ exon B/GPX1 exon 2 as ESRD risk determinants were also detected. These findings pave the way for a better understanding of the genetic influences on the susceptibility to ESRD among the Han Chinese, a population with high ESRD prevalence.

Conflict of Interests

The authors declare that they have no competing interests.

Authors’ Contribution

Chia-Ter Chao, Yen-Ching Chen, Chih-Kang Chiang, Jenq-Wen Huang, Cheng-Chung Fang, Chen-Chih Chang, and Chung-Jen Yen conceived the idea and collected the required data. Chia-Ter Chao, Yen-Ching Chen, Chih-Kang Chiang, Jenq-Wen Huang, Fu-Chang Hu, Cheng-Chung Fang, Chen-Chih Chang, and Chung-Jen Yen carried out the statistical analyses, interpreted the data, and drafted the paper. All authors approved the submission of this paper.

Acknowledgments

This work was supported in part by Department of Medical Research in National Taiwan University Hospital. The authors thank the staff in charge of health checkup and dialysis services for collecting clinical data and blood specimens from the volunteers. This study was funded by National Taiwan University Hospital (NTUH-92S027), the Ta-Tung Kidney Foundation, and the Mrs. Hsiu-Chin Lee Kidney Research Fund. The authors are extremely grateful to Dr. Fu-Chang Hu of International Harvard Statistical Consultant Agency for his excellent contribution to their analysis in this paper.

References

“Précis: an introduction to end-stage renal disease in the U.S.,” American Journal of Kidney Diseases, vol. 63, no. 1, pp. e169–e192, 2014.

B. P. Oberg, E. McMenamin, F. L. Lucas et al., “Increased prevalence of oxidant stress and inflammation in patients with moderate to severe chronic kidney disease,” Kidney International, vol. 65, no. 3, pp. 1009–1016, 2004.

T. Heitzer, T. Schlinzig, K. Krohn, T. Meinertz, and T. Münzel, “Endothelial dysfunction, oxidative stress, and risk of cardiovascular events in patients with coronary artery disease,” Circulation, vol. 104, no. 22, pp. 2673–2678, 2001.

R. Thanan, S. Oikawa, Y. Hiraku et al., “Oxidative stress and its significant roles in neurodegenerative diseases and cancer,” International Journal of Molecular Sciences, vol. 16, no. 1, pp. 193–217, 2014.

A. Glasauer and N. S. Chandel, “Targeting antioxidants for cancer therapy,” Biochemical Pharmacology, vol. 92, no. 1, pp. 90–101, 2014.

H. F. Tbahriti, A. Kaddous, M. Bouchenak, and K. Mekki, “Effect of different stages of chronic kidney disease and renal replacement therapies on oxidant-antioxidant balance in uremic patients,” Biochemistry Research International, vol. 2013, Article ID 358985, 6 pages, 2013.

R. P. Brandes, “Endothelial dysfunction and hypertension,” Hypertension, vol. 64, no. 5, pp. 924–928, 2014.

T. Nomiyama, Y. Tanaka, L. Piao et al., “The polymorphism of manganese superoxide dismutase is associated with diabetic nephropathy in Japanese type 2 diabetic patients,” Journal of Human Genetics, vol. 48, no. 3, pp. 138–141, 2003.

H. Al-Kateb, A. P. Boright, L. Mirea et al., “Multiple superoxide dismutase 1/splicing factor serine alanine 15 variants are associated with the development and progression of diabetic nephropathy: the diabetes control and complications trial/epidemiology of diabetes interventions and complications genetics study,” Diabetes, vol. 57, no. 1, pp. 218–228, 2008.

A. Crawford, R. G. Fassett, J. S. Coombes et al., “Glutathione peroxidase, superoxide dismutase and catalase genotypes and activities and the progression of chronic kidney disease,” Nephrology Dialysis Transplantation, vol. 26, no. 9, pp. 2806–2813, 2011.

C. T. Chao, Y. C. Chen, C. K. Chiang et al., “Sequence variants of peroxisome proliferator-activated receptor-gamma gene and the clinical courses of patients with end-stage renal disease,” Disease Markers, vol. 2015, Article ID 763459, 7 pages, 2015.

C. T. Chao, S. Y. Lee, W. S. Yang et al., “Serum vitamin D levels are positively associated with varicella zoster immunity in chronic dialysis patients,” Scientific Reports, vol. 4, article 7371, 2014.

C.-T. Chao, V.-C. Wu, C.-F. Lai et al., “Advanced age affects the outcome-predictive power of RIFLE classification in geriatric patients with acute kidney injury,” Kidney International, vol. 82, no. 8, pp. 920–927, 2012.

C. T. Chao, Y. H. Hsu, P. Y. Chang et al., “Simple self-report FRAIL scale might be more closely associated with dialysis complications than other frailty screening instruments in rural chronic dialysis patients,” Nephrology, vol. 20, no. 5, pp. 321–328, 2015.

C. T. Chao, H. B. Tsai, C. Y. Wu, Y. F. Lin, N. C. Hsu, J. S. Chen et al., “The severity of initial acute kidney injury at admission of geriatric patients significantly correlates with subsequent in-hospital complications,” Scientific Reports, vol. 5, article 13925, 2015.

H. Zhang, S. Zhu, J. Chen et al., “Peroxisome proliferator-activated receptor γ polymorphism Pro12Ala is associated with nephropathy in type 2 diabetes: evidence from meta-analysis of 18 studies,” Diabetes Care, vol. 35, no. 6, pp. 1388–1393, 2012.

Consortium TBCA, “Commonly studied single-nucleotide polymorphisms and breast cancer: results from the breast cancer association consortium,” Journal of the National Cancer Institute, vol. 98, no. 19, pp. 1382–1396, 2006.

A. Crawford, R. G. Fassett, D. P. Geraghty et al., “Relationships between single nucleotide polymorphisms of antioxidant enzymes and disease,” Gene, vol. 501, no. 2, pp. 89–103, 2012.

R. J. Hung, P. Boffetta, P. Brennan et al., “Genetic polymorphisms of MPO, COMT, MnSOD, NQO1, interactions with environmental exposures and bladder cancer risk,” Carcinogenesis, vol. 25, no. 6, pp. 973–978, 2004.

S. J. Lee and M. G. Choi, “Association of manganese superoxide dismutase gene polymorphism (V16A) with diabetic macular edema in Korean type 2 diabetic patients,” Metabolism: Clinical and Experimental, vol. 55, no. 12, pp. 1681–1688, 2006.

T. Kangas-Kontio, S. Vavuli, S. J. Kakko et al., “Polymorphism of the manganese superoxide dismutase gene but not of vascular endothelial growth factor gene is a risk factor for diabetic retinopathy,” British Journal of Ophthalmology, vol. 93, no. 10, pp. 1401–1406, 2009.

S. Kakko, M. Päivänsalo, P. Koistinen, Y. A. Kesäniemi, V. L. Kinnula, and M. J. Savolainen, “The signal sequence polymorphism of the MnSOD gene is associated with the degree of carotid atherosclerosis,” Atherosclerosis, vol. 168, no. 1, pp. 147–152, 2003.

J. C. Charniot, A. Sutton, D. Bonnefont-Rousselot et al., “Manganese superoxide dismutase dimorphism relationship with severity and prognosis in cardiogenic shock due to dilated cardiomyopathy,” Free Radical Research, vol. 45, no. 4, pp. 379–388, 2011.

C. Namikawa, Z. Shu-Ping, J. R. Vyselaar et al., “Polymorphisms of microsomal triglyceride transfer protein gene and manganese superoxide dismutase gene in non-alcoholic steatohepatitis,” Journal of Hepatology, vol. 40, no. 5, pp. 781–786, 2004.

G. Musso, R. Gambino, J. H. Tabibian et al., “Association of non-alcoholic fatty liver disease with chronic kidney disease: a systematic review and meta-analysis,” PLoS Medicine, vol. 11, no. 7, Article ID e1001680, 2014.

A. Rosenberger, T. Illig, K. Korb et al., “Do genetic factors protect for early onset lung cancer? A case control study before the age of 50 years,” BMC Cancer, vol. 8, no. 1, article 60, 2008.

M. Nemoto, R. Nishimura, T. Sasaki et al., “Genetic association of glutathione peroxidase-1 with coronary artery calcification in type 2 diabetes: a case control study with multi-slice computed tomography,” Cardiovascular Diabetology, vol. 6, no. 1, article 23, 2007.

G. Dutkiewicz, A. Binczak-Kuleta, A. Pawlik et al., “Lack of association of C599T polymorphism in the glutathione peroxidase (GPX1) gene with delayed graft function, acute kidney graft rejection and chronic allograft nephropathy,” Annals of Transplantation, vol. 15, no. 3, pp. 30–34, 2010.

A. Möllsten, A. Jorsal, M. Lajer, N. Vionnet, and L. Tarnow, “The V16A polymorphism in SOD2 is associated with increased risk of diabetic nephropathy and cardiovascular disease in type 1 diabetes,” Diabetologia, vol. 52, no. 12, pp. 2590–2593, 2009.

A. K. Manning, M. LaValley, C.-T. Liu et al., “Meta-analysis of gene-environment interaction: joint estimation of SNP and SNP × environment regression coefficients,” Genetic Epidemiology, vol. 35, no. 1, pp. 11–18, 2011.

A. AlSaleh, S. D. O'Dell, G. S. Frost et al., “Interaction of PPARG Pro12Ala with dietary fat influences plasma lipids in subjects at cardiometabolic risk,” Journal of Lipid Research, vol. 52, no. 12, pp. 2298–2303, 2011.

M. A. Bostrom, P. Perlegas, L. Lu et al., “Relevance of the ACTN4 gene in African-Americans with non-diabetic end-stage renal disease,” American Journal of Nephrology, vol. 36, no. 3, pp. 252–260, 2012.

M. A. Bostrom, W. H. L. Kao, M. Li et al., “Genetic association and gene-gene interaction analyses in African American dialysis patients with nondiabetic nephropathy,” American Journal of Kidney Diseases, vol. 59, no. 2, pp. 210–221, 2012.

X. Wen, Y. Li, and Y. Liu, “Opposite action of peroxisome proliferator-activated receptor-γ in regulating renal inflammation: functional switch by its ligand,” The Journal of Biological Chemistry, vol. 285, no. 39, pp. 29981–29988, 2010.

Z. Yousefipour, A. Oyekan, and M. Newaz, “Interaction of oxidative stress, nitric oxide and peroxisome proliferator activated receptor γ in acute renal failure,” Pharmacology & Therapeutics, vol. 3, pp. 436–445, 2010.