SIRT1 rs10823108 and FOXO1 rs17446614 responsible for genetic susceptibility to diabetic nephropathy

SIRT1 and FOXO1 play an important role in the pathogenesis of diabetic nephropathy (DN). However, the association between genetic polymorphisms and susceptibility to type 2 DN (T2DN) has not been explored. In this study, a total of 1066 patients with type 2 diabetes mellitus (T2DM) (413 without and 653 with DN) were enrolled. The genotypes of three htSNPs (rs3818292, rs4746720, rs10823108) within SIRT1 and two htSNPs (rs2721068, rs17446614) in FOXO1 were determined by PCR-RFLP. HbA1C, LDL, HDL, TC, and TG levels were also examined. SIRT1 rs10823108 AA genotype was significantly associated with a decreased risk of DN (OR = 0.60, 95%CI: 0.38–0.97), while GA genotype (OR = 1.77, 95%CI: 1.33–2.35) and AA genotype (OR = 2.32, 95%CI: 1.25–4.34) of FOXO1 rs17446614 was associated with an increased T2DN risk. The interactions among rs1744 6614, BMI and duration of diabetes (OR: 2.63, 95%CI: 1.23–4.31) were also observed. Subsequent haplotype analysis revealed that two haplotype defined by AC (OR: 1.50, 95%CI: 1.15–1.94) and AT (OR: 1.79, 95%CI: 1.06–2.80) within FOXO1 gene may increase the risk of T2DN. In conclusion, genetic variant rs10823108 in SIRT1 and variant rs17446614 in FoxO1 may contribute to the risk of DN in T2DM patients.

suggested that the FoxO1 activity was significantly reduced and with a concomitant decrease in the expression of FoxO1 target gene, catalase in diabetic kidney 16 . The FoxO1 downregulation correlated with an increase in the generation of malondialdehyde (MDA), a decrease in the activity of SOD and an increase in the expression of collagen IV and fibronectin proteins in renal cortex 10 . These results indicate that SIRT1/FOXO1 axis-mediated signaling pathway may be considered as a potential therapeutic target for DN. Therefore, SIRT1/FOXO1 are a very plausible candidate gene potentially contributing to the genetic polymorphism of DN. The limited finding identified 4 SNPs and haplotype consisting of the 11 SNPs within SIRT1 locus were nominally or significantly associated with DN 17 . Meanwhile,FOXO1 gene has been found to be involved in the susceptibility to type 2 diabetes in German 18 . However, it has not been reported that whether SIRT1 or FOXO1 genetic polymorphism was associated with DN in Chinese population.
In order to explore the genetic role of SIRT1 and FOXO1, and environmental interactions in diabetic nephropathy susceptibility in Chinese Han subjects with T2DN, we investigated the association between single nucleotide polymorphisms (SNPs) within the SIRT1 and FOXO1 gene and DN in Chinese Han subjects.

Results
Demographic and clinical characteristics of study samples. Baseline characteristics (including sociodemographic factors, diabetes related variables, lifestyle behaviors, biochemical variables, and history of diseases) of study subjects (653 T2DN cases and 413 T2DM without DN controls) were illustrated in Table 1. The mean age among cases (62.46 ± 11.04) was significantly older than that of controls (58.35 ± 12.07). The mean durations of diabetes were significantly different between cases (11.54 ± 7.47) and controls (14.29 ± 4.95). T2DN cases had significant higher BMI, HbA1C, total cholesterol level, and lower HDL-C level than controls. In addition, increased risks of T2DN were associated with smoking (OR = 1.40, 95% CI:1.04, 1.88) and hypertension (OR = 1.95, 95% CI:1.51, 2.50).
Effects of htSNPs in SIRT1, FOXO1 and T2DN risk. No significant deviations from Hardy-Weinberg equilibrium for the five polymorphisms was found among the 413 controls (P > 0.05). The results of single locus analysis revealed that AA genotype carriers of SIRT1 gene locus rs10823108 had decreased risk of suffering from DN (OR = 0.60, 95%CI: 0.38-0.97) than GG genotype carriers. And the decreased risk were still significant in the recessive genetic model (OR = 0.58, 95%CI: 0.37-0.91). GA genotype (OR = 1.77, 95%CI: 1.33-2.35) and AA genotype (OR = 2.32, 95%CI: 1. 25-4.34) of FOXO1 gene locusrs17446614 had increased risk of suffering from T2DN than GG genotype carriers. The significant association between T2DN risk and rs17446614 were significant in both dominant genetic model (OR = 1.86, 95%CI: 1.42-2.43) and recessive genetic model (OR = 1.96, 95%CI: 1.06-3.63). After being adjusted for age, gender, duration of diabetes, BMI, HbA1C, total cholesterol, triglycerides, HDL-C, LDL-C, hypertension, DM family history, smoking and drinking status, all the associations above were still significant. However, no significant difference was observed in the frequency distribution of rs4746720, rs3818292 and rs2721068 polymorphisms between T2DN patients and controls (Table 2).  Table 3. Stratification analysis of the five SNPs polymorphisms and DN susceptibility. a Adjusted for age, gender,duration of diabetes, BMI, HbA1C, total cholesterol, triglycerides, HDL-C, LDL-C,history of Hypertension, DM family history, smoking and drinking (the stratified factor in each stratum excluded). b P value for heterogeneity test.
to individuals with 0-2 mutation allele. And there was no increased dose-dependent manner on the combined effect of the five htSNPs.

The interaction of SIRT1 and FOXO1 gene on the pathogenesis of T2DN. Multifactor
Dimensionality Reduction (MDR) analysis was performed to analyze the gene-environment factors interaction with SNPs of SIRT1 and FOXO1(rs3818292, rs4746720, rs10823108, rs2721068, and rs17446614), age, gender, duration of diabetes, BMI, HbA1C, total cholesterol, triglycerides, HDL-C, LDL-C, Hypertension, DM family history, smoking and drinking status (

Discussion
It has been reported that genetic predisposition is one of the main risk factors for the development of DN 19 .
Evidence of familial aggregation of diabetic nephropathy in T2DM and the heritability confirm that genetic factors play a significant role in the development of DN 20,21 . There is now an increasing body of data to suggest that genetic factors may contribute to the initiation and development of diabetic nephropathy. Nevertheless, there were only a very limited number of studies that have characterized SIRT1 genes for DN in Japanese Caucasians, while no study focused on the association between FOXO1 gene and diabetic nephropathy.
In the present study, we investigated the role of SIRT1 and FOXO1 variants and their haplotypes in T2DN among unrelated Chines Caucasians. Our results suggested that rs10823108 AA genotype within SIRT1 has significant associations with about a 0.60-fold decreased risk of DN (OR = 0.60,95%CI:0.38-0.97) than GG genotype    Table 6. Interaction results between the five SNPs and risk factors by MDR.
carriers, even after adjustments for age, gender, duration of diabetes, BMI, HbA1C,total cholesterol, triglycerides, HDL-C, LDL-C, history of Hypertension, DM family history, smoking and drinking status (OR = 0.57, 95%CI: 0.34-0.94), which is similar with the result in Japanese population 17 . We can't see any significant difference between rs4746720 and rs3818292 SNPs and diabetic nephropathy in univariate analysis and multivariate analysis, which is different with the results in Japanese, which may be due to ethnic differences and genetic backgrounds. In the next stratified analysis, rs3818292 AG + GG genotypes had a significant increased risk of T2DN in subjects with TC level > 5 (OR: 2.11; 95%CI: 1.06-4.21).
Haplotype was accepted as a group of correlated SNPs that are located in the same homologous chromosome, and passed on to descendants together as a whole 23 . Haplotype contained multiple SNPs information, and the statistical power of haplotype analysis was stronger than analysis of single SNP 24 . Furthermore, we detected that G rs17446614 T rs2721068 (OR:0.77, 95%CI:0.63-0.93) and G rs17446614 C rs2721068 (OR:0.78, 95%CI:0.65-0.93) within FOXO1 could reduce the risk of DN in T2DM patients. Therefore, our results suggest that the SNPs within SIRT1 and FOXO1 might be involved in the process of diabetic nephropathy in Henan Han Chinese population. In this study, the gene-reproductive factors interaction on T2DN susceptibility was examined by using a MDR method. Anemically significant interaction was found for rs1746614, BMI and duration of diabetes (OR:4.01, 95%CI:3.02-5.31). To our knowledge, this is the first study to examine the effect of FOXO1 polymorphisms on T2DN risk with an emphasis on the gene-environment factors interaction.
Although, the mechanism by which the SIRT1 polymorphism contributes to conferring susceptibility to diabetic nephropathy remains to be elucidated, combining the present finding with a previous report, SIRT1 and FOXO1 may be considered a good new candidate gene for diabetic nephropathy. The effects of SIRT1 polymorphisms on susceptibility to DN might be mediated by differences in the metabolic state among individuals, including glycemic control, obesity, blood pressure, respectively [25][26][27] . The evidence in Wild-type (C57BL/6 J) and SIRT1 transgenic mice (C57BL/6-Actb tm3.1(Sirt1)Npa /J) showed SIRT1 mediated protection against renal and retinal injury in DN 28 . FOXO1, as an important target for insulin and growth factor signaling, has been reported to play key roles in regulating metabolism, cell proliferation, cell cycle, apoptosis and oxidative stress response by regulating the expression of functionally important target genes 29,30 . It is shown that GWAS in Pima Indians identified variation within FOXO1A that modestly associated with early-onset T2DN. Common variation in FOXO1A may modestly affect risk for T2D and obesity in American Indians 31 . Interestingly, our previous studies have confirmed that the transcriptional activity of FOXO1 was decreased in both renal cortex of DN rats and MCs cultured under high-glucose (HG) conditions, which indicated that FOXO1 may be considered a good new candidate gene for diabetic nephropathy 14,15 .
Our study has several notable strengths. First, although the role of SIRT1 and FOXO1 in multiple diseases has been reported, this is the first study investigating genetic variation of FOXO1 and SIRT1 in relation to T2DN susceptibility. Second, in order to draw a convincing conclusion, only patients who had suffered T2DM for more than 8 years and no diabetic retinopathy were eligible to participate in the control group, because the duration of DM plays a very important role in the development of DN, and diabetic retinopathy has a complicated relationship with DN 32 . To increase comparability between DN cases and controls, the DN cases and the controls were much more homogeneous in terms of gender, age, and durations of diabetes, and family history, after excluding subjects with extreme characteristics. Third, our genotyping results were accurate and reliable, since more than 30% of the samples were genotyped in duplicate. 10% of the samples were randomly selected from both cases and controls to confirm the genotyping results, and the result of confirmation were found to be 100% concordant.
Some limitations of this study should also be noted. One limitation was the limited size sample in this study, which might be the reason why we didn't see any statistical significance in gene-gene interaction analysis. Secondly, although we identified two novel SNPs that are susceptibility variants of DN in a Han Chinese population in Henan, screening more loci within SIRT1 gene and FOXO1 gene to further clarify the association between them and DN is clearly warranted. In order to validate our findings, further research needs to be conducted, including further functional evaluation with other ethnic populations and animal experiments.
BMI might be a factor in affecting SIRT-1 activation and reversing epigenetic metabolic imprinting or silencing. In the recent report, the clinical evidence also support the relationship of BMI (Epicardial fat thickness) and circulating SIRT1 33 . In our study, the data of BMI was 26.19 ± 2.55 kg/m 2 in case group while 25.25 ± 2.74 kg/m 2 in control group. Interestingly, their findings are consistent with our results. Of couse, more evidence will be confrimed in animal experiments. Additionally, the activation of SIRT1 might be associated with Zinc levels. In the study of Sheline C.T., he/she found that SIRT1-mediated NAD + loss in Zn 2+ , STZ, or cytokine toxicities of MIN6, and in NOD or DM animal models 34 . These findings indicated that β-cell Zn 2+ levels play a regulator in DM and its vascular complications, including DN. The interesting evidence need to be provided in our following experiments.
Despite the several limitations mentioned above, to the best of our knowledge, our data first investigated the association between SIRT1 and FOXO1 gene and T2DN in Chinese Han population. This study has made a valuable contribution to draw a plausible link between SIRT1 and FOXO1 gene polymorphisms and DN. However, replication in additional DM patient cohorts, case control studies, and further functional work to dissect out the role of SIRT1 and FOXO1 polymorphisms in vitro conditions, are clearly needed. Furthermore, on the basis of the current genotyping and clinical staging, prospective studies in the disease outcome in different stages, especially in patients with early diabetic nephropathy group, further verify the role of each gene locus in diabetic nephropathy, which will be more interesting and meaningful.
We identified two novel SNPs (rs10823108 in SIRT1, rs17446614 in FOXO1) and two haplotypes within FOXO1 gene (defined by GT and GC), that are susceptibility variants of DN in a Han Chinese population in Henan. The findings reveal that SIRT1/FOXO1 may involve in the initiation and development of DN in Han Chinese patients with type 2 diabetes. Our results can provide new insights into the etiology of DN.

Materials and Methods
Study subjects. A case-control study was performed in Han Chinese population located in Henan province. Individuals with type 2 diabetes and aged over 20 years were recruited from June 2012 to Sep 2013 using the American Diabetes Association (ICD-9-CM, Diagnosis code 250) criteria for diagnosis of type 2 diabetes. A total of 1066 type 2 diabetes subjects, recruited from the first affiliated hospital of Zhengzhou University, located in the middle part of China, were involved in this study. All patients with type 2 diabetes were of Han Chinese origin. Subjects who tested positive for anti-glutamic acid decarboxylase (GAD) antibodies and those diagnosed with mitochondrial disease (mitochondrial myopathy, encephalopathy, lactic acidosis, and stroke-like episodes [MELAS]) or maturity-onset diabetes of the young (MODY), gestational diabetes and secondary diabetes were not included.
Diabetic nephropathy was determined by estimated 24 h albumin excretion rates (AERs) or urinary albumin-to-creatinine ratios (ACR) 35,36 . AERs or ACR were measured twice, and the mean value was recorded for each patient. Patients having microalbuminuria, indicated 30 mg/24 h ≤ AERs < 300 mg/24 h or 30 mg/g ≤ ACR < 300 mg/g, and with background or greater diabetic retinopathy, or proteinuria indicated AERs ≥ 300 mg/24 h or ACR ≥ 300 mg/g, or ESRD diagnosed at least 5 years before initiating renal replacement therapy, were diagnosed with diabetic nephropathy (n = 653). Patients with no diabetic retinopathy, and normal albuminuria, indicated AERs < 30 mg/24 h or ACR < 30 mg/g, and a long T2DM duration over 8 years were selected as a control group (n = 413). All the patients participating in this study provided written informed consent, and the study protocol was approved by the ethics committees of the first affiliated hospital of Zhengzhou University. We confirmed that all experiments were performed in accordance with relevant guidelines and regulations.

Measurements.
Family history was defined as the presence of T2DM in either a parent or a sibling; smoking was defined as having smoked more than 1 cigarette every day and at least one year in one's lifetime; similarly, drinking was defined as alcohol intake more than once per week. Anthropometric measurements included weight (kg); height (cm); BMI was calculated as weight divided by squared height. The fasting serum biochemical profile includes total cholesterol (mmol/L); LDL (mmol/L); HDL (mmol/L), and triglycerides (mmol/L), were determined by using cobas c8000 clinical Chemistry System (Roche, German). After an early morning spot urine sample collection, urinary albumin and creatinine concentrations were measured by Afinion AS100 Analyzer and the urine albumin-to-creatinine ratio (UACR) was defined as the urinary albumin value divided by the creatinine concentration (mg/g). SNP selection and genotyping. Reference sequences of SIRT1 and FOXO1 were identified with National Center for Biotechnology Information (NCBI) database (http://www.ncbi.nlm.nih.gov/). Haplotype-tagging SNPs (htSNPs) were used to analyze the SIRT1 and FOXO1 polymorphisms globally. htSNPs selection was conducted with the pairwise option of the Haploview software 4.2, minor allele frequency (MAF) > 0.1 and the pairwise correlation was set as r 2 > 0.8, based on the Chinese Han population data of HapMap (HapMap Data Rel 28 Phase II + III, August 10, on NCBI B36 assembly, dbSNP b126). In summary, five htSNPs (rs3818292 A > G, rs4746720 T > C, rs10823108 G > A in SIRT1 and rs2721068 T > C, rs17446614 G > A in FOXO1) that captured most of the known common SNPs located in the chromosome locus transcribed into SIRT1 and FOXO1 and its flanking region (2000 bp upstream and 2000 bp downstream, respectively) were selected.
DNA was extracted from a peripheral blood sample using a QIAamp (Qiagen, Germany) kit, and analyzed by electrophoresis in 0.8% agarose gels stained with ethidium bromide and visualized in a Gel Doc 2000 (BIORAD CA, USA). DNA concentration was determined using a VICTOR3 1420 spectrophotometer (Perkin-Elmer, Germany). The five htSNPs were genotyped with restriction fragment length polymorphism method (PCR-RFLP). Primers and restriction endonucleases were shown in Table 7. The volume of PCR was 15 μl, which contained 7.5 μl 2 × Tap PCR Master Mix, 0.8 μM each primer, 50 ng DNA, and 5.7 μl deionized water. The PCR profile included an initial melting step of 95 °C for 5 min, followed by 35 cycles of denaturation at 95 °C for 30 s, annealing at 59 °C, 59 °C, 59 °C, 57 °C, and 57 °C for rs3818292, rs4746720, rs10823108, rs2721068, and rs17446614, respectively for 45 s and extension at 72 °C for 1 min and a final extension step of 72 °C for 10 min. PCR products were digested for 3-16 h at 37 °C by VspI, Hpy1881, Hpy1881, TspRI, and CviJI (Fermentas) for rs3818292, rs4746720, rs10823108, rs2721068, and rs17446614, respectively. The digested PCR products were separated on 3% agarose gel electrophoresis and stained with ethidium bromide. For each SNP in the groups of study, a concordance of 100% was observed when 30-50% of the samples were genotyped in duplicate. 10% samples were randomly selected from both cases and controls to confirm the genotyping results. The direct sequencing (BGI Sequencing, Beijing) was also utilized for the genotype confirmation, and the result of confirmation were found to be 100% concordant.
Statistical analysis. Difference in the distribution of demographic characteristics and genotype frequencies between cases and controls were evaluated using the χ 2 -test (Chi-square test). The associations between rs3818292, rs4746720, rs10823108, rs2721068, rs17446614 genotypes and risk of T2DN were estimated by computing the odds ratio (ORs) and their 95% confidence interval(CI). Hardy-Weinberg equilibrium (HWE) for the genotype distribution of each SNP was tested by a goodness-of-fit Chi squared test to compare the observed genotype frequencies with the expected ones among the control subjects. The data were further stratified by demographic characteristics and biochemical factors to evaluate the stratum variable-related ORs among the five SNPs. MDR methods were used to detect the gene-environment factors interaction. The online SHEsis (http:// analysis.bio-x.cn/myAnalysis.php) was taken for Haplotype Analysis. Statistical analysis was performed by using SPSS 17.0 and SAS 9.1 software package. All tests were two sided and a P value < 0.05 was considered to be statistically significant.  Table 7. Detail of PCR Primers and restriction endonucleases primers. a based on the Chinese Han population data of HapMap.