Association of Nuclear Factor Erythroid-2-Related Actor 2 Gene Polymorphisms with Diabetic Nephropathy in Chinese Patients.

AIM: Nuclear factor erythroid-2-related factor 2 (NRF2) has emerged as a therapeutic target in many diseases. To explore this further, we evaluated the relationships between the -617C/A (rs6721961) polymorphisms within the NRF2 promoter and diabetic nephropathy (DN) in Chinese Han patients with type 2 diabetes mellitus (T2DM).
METHODS: A total of 883 subjects with T2DM (500 without and 383 with DN) were enrolled in this study. Multivariable linear regression models were carried out to assess the association of DN with the -617C/A (rs6721961) polymorphisms.
RESULTS: The AA genotype frequencies in patients with DN were significantly lower than those in patients without DN (χ2 = 8.04, p = 0.018). Multivariate logistic regression analyses showed that individuals with the AA genotype had a significantly lower risk for DN (OR 0.52; 95% CI 0.28, 0.94; p = 0.029) than those with the CC genotype. Moreover, AA carriers had a significantly lower risk of DN (OR 0.46; 95% CI 0.26, 0.82; p = 0.009) relative to those with the CC + CA genotype, even after adjusting for known DN risk factors.
CONCLUSION: Our study indicated that the -617C/A polymorphism within the NRF2 promoter was significantly associated with DN in Chinese Han patients with T2DM.

Introduction

Diabetic nephropathy (DN) is one of the microvascular complications of diabetes and is the most common cause of end-stage renal disease (ESRD), and the cost for the treatment of diabetic patients with ESRD is increasing.1 Statistically, approximately 35% of people with diabetes will develop evidence of DN, and a substantial number of them would require dialysis.2 Therefore, an improved understanding of the factors involved in the onset and progression of DN is urgently required. To date, the major risk factors leading to DN include hyperglycemia, disorder of the rennin-angiotensin system, oxygen stress, inflammatory reaction and gene polymorphism.3 Recent evidence also suggests that nuclear factor erythroid-2-related factor 2 (NRF2) plays an important role in the development of type 2 diabetes mellitus (T2DM) as well as its complications.4

The master antioxidant regulator NRF2 is a member of the cap-n-collar basic leucine zipper family that regulates the expression of many antioxidant genes in response to oxidative stress.5,6 Under physiological conditions, Kelch-like ECH associating protein 1 (Keap1) combines with NRF2, triggering its proteasomal degradation.7 In the presence of oxidative stress, NRF2 escapes such degradation and undergoes translocation to the nucleus, where it binds to antioxidant responsive elements (AREs) and upregulates the expression of various downstream genes encoding antioxidant enzymes, including heme oxygenase-1, NAD(P)H quinine oxidoreductase, superoxide dismutase 3, and peroxiredoxin 1, to avoid oxidative damage.8–10 To date, NRF2 has emerged as a potential therapeutic strategy for various diseases because it can induce antioxidant enzymes and other cytoprotective enzymes.4,11–13 A few single nucleotide polymorphisms (SNPs) have been identified in the NRF2 gene, such as rs264723, rs13001694, rs10497511, rs1806649.4 In particular, SNP −617C/A (rs6721961), a variant of the NRF2 gene in the upstream promoter region, is considered to play a pivotal role in the positive feedback loop of transcriptional activation of the NRF2 gene to regulate the NRF2 protein level.14

The −617C/A polymorphism has been gaining attention in the field of pharmacology and especially in the chemoprevention of diseases such as cancer and metabolic and cardiovascular diseases.15–18 A polymorphism in the NRF2 promoter (−617 C/A) diminished the promoter activity compared with the wild type (−617 CC) and was associated with a higher risk for developing acute lung injury, a disease related to reduced antioxidative cell defense.19 Furthermore, the rs6721961 polymorphism has been associated with oxidative stress and the risk of newly diagnosed T2DM.16 However, the role of the −617C/A (rs6721961) polymorphism within NRF2 in the regulation of DN in humans has never been investigated.

In this study, therefore, we evaluated the potential associations of the −617C/A (rs6721961) polymorphism with DN in a relatively large Chinese population.

Patients and Methods Study Population

This cross-sectional genotyping study enrolled consecutive Han Chinese patients with T2DM in the Department of Endocrinology and Metabolism, First Affiliated Hospital of Nanchang University, Nanchang, Jiangxi Province, China, between June 2011 and March 2015. Family trees were used to ensure their Han Chinese ethnic background. T2DM was diagnosed according to the 2003 American Diabetes Association diagnostic criteria for diabetes, based on a fasting blood glucose (FBG) >7 mmol/l, a causal blood glucose >11.1 mmol/l, or a postprandial 2-h blood glucose (2hBG) >11.1 mmol/l following a 75-g oral glucose tolerance test, or a history of therapy for diabetes. Subjects were divided into two groups: without diabetic nephropathy (n = 500) and with diabetic nephropathy (n = 383) according to their 24 hours microalbuminuria (MAU) and estimated glomerular filtration rate (eGFR). The patients without diabetic nephropathy who had at least 5 years of known duration of diabetes had no albuminuria (MAU < 30 mg/24 h) and an eGFR > 60 mL min−1 1.73 m−2 and were not receiving anti-hypertension treatment. The patients with diabetic nephropathy had overt albuminuria (MAU > 300 mg/24 h) or eGFR < 60 mL min−1 1.73 m−2 (no end-stage renal disease or kidney transplantation), without any clinical or laboratory evidence of other kidney diseases. All patients underwent a complete eye examination that included dilated retinal examination and fundus photography or fundus fluorescein angiography. Diabetic retinopathy was evaluated by an experienced ophthalmologist. This study was approved by the ethics committee of the First Affiliated Hospital of Nanchang University (no. CDYFY2014009) and written informed consent was obtained from all individuals.

Biochemical Analyses

All of the patients underwent a standardized clinical and laboratory evaluation. Fasting blood samples were taken for measurement of FBG, hemoglobin A1c (HbA1c), total cholesterol (TC), low density lipoprotein cholesterol (LDL-C), high density lipoprotein (HDL-C), total triglyceride (TG), blood urea nitrogen (BUN), and serum creatinine. FBG, TC, HDL-C, LDL-C, TG, and serum creatinine were analyzed using the automated Olympus AU5421 chemistry analyzer (Olympus, Shizuoka, Japan). HbA1c content was measured using a Bio-Rad D-10 glycated hemoglobin analyzer (Bio-Rad, Hercules, USA). Consider eGFR (mL min−1 1.73 m−2) = 175 × [serum creatinine (mg/dL)−1.234 × age (years)−0.179] × (0.79, if female).20

Genotyping

For genotype analysis, genomic DNA was extracted from peripheral blood leukocytes of each individual using a DNA isolation kit (Bioteke, Beijing, China). The genomic DNA was subjected to polymerase chain reaction (PCR) with the following primers: Primer-F: 5ʹ-GGTTCCCGTTTTTCTCCCAGCTCT-3ʹ, Primer-R: 5ʹ-GGGTGGGGGGGGCTAAAGATTT-3ʹ. PCR was performed in a 20 μL reaction volume containing 150 ng of template DNA, 0.2 μmol/L of forward and reverse primers, 0.225 U Taq DNA polymerase, 200 μmol/L of each deoxynucleotide phosphate, and 2.0 mmol/L Mg2+ on the ABI 7300 PCR system (Applied Biosystems, Foster City, CA, USA). The amplification protocol consisted of 94°C for 5 min, followed by 35 cycles of denaturation at 94°C for 20 s, annealing at 58°C for 30 s, and extension at 72°C for 30 s, with a final extension at 72°C for 5 min. The PCR products were randomly selected for DNA sequencing analysis by Nanjing GenScript Biotech Corporation.

Statistical Analyses

Continuous variables were analyzed using a t-test and are presented as the mean ± standard deviation (SD). For categorical variables, a chi-square test or Fischer exact test was used, and data are presented as percentages. Categorical variables and Hardy-Weinberg equilibrium were analyzed using a chi-square test. Intergroup comparisons of genotype and allele distribution were analyzed using Fisher’s exact test and chi-square test, respectively. A two-sided p value < 0.05 was considered statistically significant.

Multivariable linear regression models were carried out to assess the association of diabetes events with the specific polymorphism. Odds ratios (OR) and 95% confidence intervals (95% CI) were adjusted for known risk factors for age, sex, body mass index (BMI), SBP, DBP, FBG, HbA1c, TC, LDL-C, HDL-C, TG. Statistical analyses were performed using R software ().

Results

Clinical and Laboratory Characteristics of Patients with T2DM

A total of 883 subjects with T2DM were enrolled in the study, including 383 patients in the DN group and 500 in the without DN group. The demographic and clinical characteristics for the DN and without DN groups are shown in Table 1. Participants with DN had a higher age, male gender, SBP, proportion of retinopathy, and levels of BUN, microalbuminuria, and serum creatinine than the without DN subjects. The frequencies of the AA genotype were significantly lower in patients with DN than in the without DN subjects. There were no significant differences in HbA1c percentages, FBG, BMI, DBP, TC, triglyceride, LDL-C, or HDL-C between the two groups.

Table 1

Clinical and Laboratory Characteristics of Patients with Type 2 Diabetes with and without Diabetic Nephropathy

	Total (n=883)	Without Diabetic Nephropathy (n=500)	With Diabetic Nephropathy (n=383)	P value*
Male, No.(%)	462 (52.3)	242 (48.4)	220 (57.4)	0.008
Age, year	60.6 (11.0)	59.4 (10.6)	62.1 (11.4)	<0.001
Body mass index, kg/m2	23.9 (3.5)	23.9 (3.7)	23.9 (3.3)	0.868
Nrf2 polymorphisms, No.(%)				0.018
CC	383 (43.4)	222 (44.4)	161 (42.0)
CA	434 (49.2)	231 (46.2)	203 (53.0)
AA	66 (7.5)	47 (9.4)	19 (5.0)
Retinopathy, No.(%)	245 (32.2)	110 (25.9)	135 (40.1)	<0.001
Blood pressure, mmHg
Systolic BP	135.4 (18.6)	132.8 (18.1)	138.8 (18.7)	<0.001
Diastolic BP	82.0 (11.5)	81.6 (11.6)	82.4 (11.3)	0.303
Laboratory results
FBG, mmol/L	8.5 (3.9)	8.4 (3.7)	8.7 (4.3)	0.173
Hemoglobin A1c (%)	8.7 (2.4)	8.6 (2.4)	8.8 (2.4)	0.123
Total cholesterol, mmol/L	4.8 (11.7)	5.1 (15.4)	4.5 (1.3)	0.454
Triglyceride, mmol/L	1.6 (1.7)	1.5 (1.8)	1.7 (1.5)	0.106
HDL-C, mmol/L	1.3 (3.7)	1.4 (4.8)	1.2 (1.2)	0.481
LDL-C, mmol/L	3.0 (7.8)	2.6 (0.9)	3.6 (11.8)	0.054
BUN, mmol/L	6.9 (5.2)	5.4 (1.9)	8.9 (7.1)	<0.001
Microalbuminuria, mg/24h	42.6 (54.5)	6.7 (6.0)	87.5 (54.9)	<0.001
Serum creatinine, μmol/L	87.4 (76.6)	60.1 (16.0)	122.5 (104.4)	<0.001

Notes: Data are the mean ± SD, or n (%). *P values were obtained by an unpaired Student’s t-test analysis or χ2 test, as appropriate.

Abbreviations: FBG, fasting blood glucose; HDL-C, high density lipoprotein; LDL-C, low density lipoprotein cholesterol; BUN, blood urea nitrogen.

Clinical and Laboratory Characteristics of Patients with Type-2 Diabetic Nephropathy Classified According to −617C/A Genotype

As shown in Table 2, the levels of triglyceride in the AA genotype was significantly higher than that in the patients with the CC or CA genotype. There were no significant differences in other clinical and laboratory characteristics between the genotype groups.

Table 2

Clinical and Laboratory Characteristics of Patients with Type-2 Diabetic Nephropathy Classified According to Their −617C/A Genotypes

	CC (n=161)	CA (n=203)	AA (n=19)	P value*
Male, No.(%)	103 (64.0)	221 (52.2)	11 (57.9)	0.079
Age, year	61.3 (11.1)	62.6 (11.8)	63.5 (9.0)	0.451
Body mass index, kg/m2	23.8 (3.3)	24.0 (3.4)	23.9 (3.3)	0.920
Retinopathy, No.(%)	64 (42.4)	65 (38.9)	6 (31.6)	0.607
Blood pressure, mmHg
Systolic BP	138.0 (19.5)	139.5 (18.6)	138.4 (13.6)	0.747
Diastolic BP	83.0 (12.0)	82.0 (11.1)	83.2 (8.3)	0.676
Laboratory results
FBG, mmol/L	9.1 (4.7)	8.6 (3.9)	6.9 (3.4)	0.092
Hemoglobin A1c (%)	8.8 (2.4)	8.7 (2.4)	9.3 (2.2)	0.649
Total cholesterol, mmol/L	4.6 (1.5)	4.5 (1.2)	4.1 (1.1)	0.294
Triglyceride, mmol/L	1.8 (1.5)	1.5 (1.0)	2.5 (3.5)	0.008
HDL-C, mmol/L	1.1 (0.4)	1.3 (1.6)	1.1 (0.4)	0.400
LDL-C, mmol/L	4.8 (18.1)	2.8 (1.1)	2.5 (1.0)	0.239
BUN, mmol/L	9.3 (7.6)	8.8 (7.0)	6.8 (2.6)	0.343
Microalbuminuria, mg/24h	86.1 (54.7)	88.2 (54.6)	91.5 (62.9)	0.905
Serum creatinine, μmol/L	131.4 (119.0)	115.0 (94.2)	127.2 (67.4)	0.328

Notes: Data are the mean ± SD, or n (%). *P values were obtained by an unpaired Student’s t-test analysis or χ2 test, as appropriate.

Abbreviations: FBG, fasting blood glucose; HDL-C, high density lipoprotein; LDL-C, low density lipoprotein cholesterol; BUN, blood urea nitrogen.

Association of DN with −617C/A Polymorphism in Patients with T2DM

As shown in Table 3, the genotype frequencies of the −617C/A polymorphism in the control group were in Hardy-Weinberg equilibrium (χ2 = 1.40, p = 0.236). There were significant differences in the allelic frequency of the −617C/A polymorphism between DN and without DN. The frequencies of the CC, CA, and AA genotypes were 44.4, 46.2, and 9.4%, respectively, in the group without DN and 42.0, 53.0, and 5.0%, respectively, in the group with DN (χ2 = 8.04, p = 0.018).

Table 3

The Genotypic Distribution of the −617C/A Polymorphism in Type 2 Diabetic Patients with and without Diabetic Nephropathy

	Without Diabetic Nephropathy n (%)	With Diabetic Nephropathy n (%)	χ2	P value*
Genotype			8.04	0.018
CC	222 (44.4)	161 (42.0)
CA	231 (46.2)	203 (53.0)
AA	47 (9.4)	19 (5.0)

Note: *P values were obtained by an χ2 test.

Individuals carrying with the AA genotype had a significantly lower risk for DN (OR 0.52; 95% CI 0.28, 0.94; p = 0.029) than those with the CC genotype; Compared to those with the CC + CA genotype, AA carriers had a significantly lower risk of DN (OR 0.46; 95% CI 0.26, 0.82; p = 0.009) after adjusting for age, sex, BMI, SBP, DBP, FBG, HbA1c, TC, LDL-C, HDL-C, TG. (Table 4).

Table 4

Univariate and Multivariate Odds Ratios of Diabetic Nephropathy According to the −617C/A Polymorphism

	No. of Cases/Control Subjects	OR (95% CI)		P value*
		Unadjusted	Adjusted
Genotype
CC	161/222	1.00	1.00
CA	203/231	1.21 (0.92, 1.60)	1.23 (0.92, 1.64)	0.162
AA	19/47	0.56 (0.32, 0.99)	0.52 (0.28, 0.94)	0.029
CC+CA	364/453	1.00	1.00
AA	19/47	0.50 (0.29, 0.87)	0.46 (0.26, 0.82)	0.009

Notes: *Adjusted for age, sex, body mass index, systolic blood pressure, diastolic blood pressure, fasting blood glucose, hemoglobin A1c, total cholesterol, low density lipoprotein cholesterol, high density lipoprotein, total triglyceride.

Stratified Analyses by Important Covariables

Stratified analyses were performed to further assess the association between −617C/A polymorphism and DN in age and sex subgroups. None of the following variables, including age (<60 vs ≥60 years; P for interaction=0.558), sex (P for interaction=0.653) significantly modified the association between −617C/A polymorphism and DN (Table 5).

Table 5

Subgroup Analysis of Diabetic Nephropathy According to the −617C/A Polymorphism

Source	OR (95% CI)		P value*	P for Interaction*
	CC+CA	AA
Age, years				0.558
<60	1.00	0.58 (0.25, 1.37)	0.215
≥60	1.00	0.36 (0.16, 0.83)	0.016
Sex				0.653
Male	1.00	0.51 (0.23, 1.14)	0.107
Female	1.00	0.38 (0.16, 0.93)	0.034

Notes: *Adjusted for age, sex, body mass index, systolic blood pressure, diastolic blood pressure, fasting blood glucose, hemoglobin A1c, total cholesterol, low density lipoprotein cholesterol, high density lipoprotein, total triglyceride.

Discussion

NRF2 is a key regulator of DN through the suppression of oxidative stress assault,21 which is a major contributor to the development and progression of DN. To our knowledge, this is the first report providing clinical evidence that the 617C/A SNP, which is located in the promoter region of the NRF2 gene, is associated with DN in Chinese patients. Individuals with AA homozygotes had a significantly lower risk of DN than those with CC homozygotes, even after adjusting for known DN risk factors.

Previous research has paid little attention to the association between NRF2 single nucleotide polymorphisms and DN, and there are only limited data and disputable results linking the NRF2 gene with diabetes. In line with our findings, a study conducted in Mexico showed that the rs6721961 (−617C/A) polymorphism of the NRF2 gene was associated with diabetes in male subjects, AA carriers had lower glucose concentrations, and the OR revealed that A carriers had a lower risk of developing diabetes.15 However, Wang et al observed that AA carriers in the NRF2 rs6721961 polymorphism had lower total antioxidative capacity, superoxide dismutase, catalase, glutathione, and glutathione peroxidase activity and had a significantly higher risk of developing T2DM (OR 1.56; 95% CI 1.11, 2.20; p = 0.011) than those with the CC genotype.16 In the previous two studies,15,16 they evaluated the rs6721961 (−617C/A) polymorphism with T2DM versus the −617C/A polymorphism with microvascular complication of T2DM, which was assessed in the present study. Additionally, the included population and other characteristics, such as the patients who were selected and the sample size, might explain the inconsistent results, but we believe that the NRF2 rs6721961 polymorphism is worthy of further examination.

There are some reports linking NRF2 genetic variations and metabolic disease in different racial populations. Polymorphisms in the promoter region of NRF2 have been previously reported in Japanese hemodialysis patients. It was found that the mutant allele of rs6721961 (−617C/A) is significantly associated with systolic and diastolic blood pressure.17 In a study including white and African American populations, −617A variant allele carriers had significantly higher forearm vascular resistance at baseline compared with wild-type individuals in whites.22 One study observed that T2DM patients with complications had a higher frequency of four mutant NRF2 single nucleotide polymorphisms (rs2364723 C allele, rs10497511 G allele, rs1962142 A allele and rs6726395 G allele) than T2DM patients without complications in a cohort of Chinese participants.4

Oxidative stress is emerging as an important contributing factor of DN pathology and progression.3,23 NRF2 is a pivotal transcription factor and functions to upregulate the expression of antioxidant genes in response to oxidative stress. In a mouse model of diabetes induced by methylglyoxal, treatment with resveratrol, which has been proposed as an effective treatment that helps lower the risk of developing complications of diabetes, markedly improved blood glucose levels based on the oral glucose tolerance test and promoted NRF2 phosphorylation of the pancreas.24 Accordingly, the treatment of normal human renal epithelial HK-2 cells with NRF2 inducers effectively elevated the expression of aldo-keto reductase levels.25 Targeted disruption of NRF2 significantly reduced antioxidant capacity in mice and thus increased high-level susceptibility to a variety of toxic chemicals and carcinogens.26,27 Thus, lower antioxidant capacity due to lower expression of the NRF2 gene is a risk factor for complications of diabetes. The findings of the present study revealed that high activity of NRF2, as shown in AA allele carriers of the – 617C/A SNP, seems to be associated with a low risk of DN. This indicates that the AA genotype plays a prominent role in the renal protection of humans against oxidative stress. However, a more complete understanding of the mechanisms by which NRF2 protects against DN is needed.

There are some limitations in this study. First, we were not able to measure plasma reactive oxygen species concentrations in all of the involved participants, particularly in DN subjects. Second, all of the participants in our study were Han Chinese. Whether the results can be extrapolated to other populations requires further verification. Third, we analyzed only one SNP, rs6721961, and haplotypes of the promoter region only, which does not exclude any association of other regions around or within the gene. Finally, the present findings obtained in a cross-sectional study are explorative in nature, and replication in independent prospective population-based studies with different ethnicities is needed to determine whether the −617C/A polymorphism influences patients with DN.

In conclusion, the present study indicated that the −617C/A polymorphism in the NRF2 gene was significantly associated with DN in Chinese Han patients with T2DM. Further studies are necessary to validate the impact of the polymorphism in the NRF2 gene on oxidative stress and the risk of T2DM and its complications.