Journal List > J Korean Med Sci > v.32(2) > 1108640

Gong, Li, Ren, Wang, Wang, Yang, and Zhang: The FOXO1 Gene-Obesity Interaction Increases the Risk of Type 2 Diabetes Mellitus in a Chinese Han Population

Abstract

Here, we aimed to study the effect of the forkhead box O1-insulin receptor substrate 2 (FOXO1-IRS2) gene interaction and the FOXO1 and IRS2 genes-environment interaction for the risk of type 2 diabetes mellitus (T2DM) in a Chinese Han population. We genotyped 7 polymorphism sites of FOXO1 gene and IRS2 gene in 780 unrelated Chinese Han people (474 cases of T2DM, 306 cases of healthy control). The risk of T2DM in individuals with AA genotype for rs7986407 and CC genotype for rs4581585 in FOXO1 gene was 2.092 and 2.57 times higher than that with GG genotype (odds ratio [OR] = 2.092; 95% confidence interval [CI] = 1.178–3.731; P = 0.011) and TT genotype (OR = 2.571; 95% CI = 1.404–4.695; P = 0.002), respectively. The risk of T2DM in individuals with GG genotype for Gly1057Asp in IRS2 gene was 1.42 times higher than that with AA genotype (OR = 1.422; 95% CI = 1.037–1.949; P = 0.029). The other 4 single nucleotide polymorphisms (SNPs) had no significant association with T2DM (P > 0.05). Multifactor dimensionality reduction (MDR) analysis showed that the interaction between SNPs rs7986407 and rs4325426 in FOXO1 gene and waist was the best model confirmed by interaction analysis, closely associating with T2DM. There was an increased risk for T2DM in the case of non-obesity with genotype combined AA/CC, AA/AC or AG/AA for rs7986407 and rs4325426, and obesity with genotype AA for rs7986407 or AA for rs4325426 (OR = 3.976; 95% CI = 1.156–13.675; P value from sign test [Psign] = 0.025; P value from permutation test [Pperm] = 0.000–0.001). Together, this study indicates an association of FOXO1 and IRS2 gene polymorphisms with T2DM in Chinese Han population, supporting FOXO1-obesity interaction as a key factor for the risk of T2DM.

Graphical Abstract

jkms-32-264-ab001

INTRODUCTION

Type 2 diabetes mellitus (T2DM) is a common and frequent disorder, which is seriously harmful to human health with a clear genetic susceptibility (123). By the year 2030, an estimated 350 million individuals worldwide will suffer from diabetes (4), with approximately 100 million people with T2DM in China alone (5). However, this disease is usually caused by the combination of gene and environment. Therefore, it's helpful to identify vulnerable populations with high risk of T2DM by studying the either gene-gene interaction or gene-environment interaction (67), thereby managing early prevention.
As known, insulin resistance (IR) and β-cell dysfunction contribute to T2DM (8). It has been suggested that forkhead box O1 (FOXO1) is closely related to the 2 pathogenesis of T2DM. Activated FOXO1 leads to IR by regulating transcription of target gene, resulting in an increase of the gluconeogenesis and hepatic glucose production in liver; it also leads to β-cell dysfunction by decreasing β-cell proliferation and compensatory ability in islet (910). Furthermore, activity of FOXO1 is mainly regulated through insulin/insulin-like growth factor signaling pathway insulin receptor substrate (IRS)-AKT-FOXO1 (11). Furthermore, IRS2 knockout mice have a progressive development of T2DM, with IR and progressive deterioration of β-cell (12). The IRS +/− mice with heterozygous knockout FOXO1 can significantly improve sensitivity of hepatic insulin and prevent diabetes (13), as well as prevent function failure of β-cell (14). Thus, FOXO1 and IRS2 may be co-involved in IR and β-cell dysfunction, and their interaction may promote the occurrence of T2DM.
Previously, different results have been shown in association of single nucleotide polymorphisms (SNPs) at IRS2 gene and FOXO1 gene with T2DM in multiple populations (1516171819), which may be due to the various ethnic groups. However, it is still unclear whether the FOXO1-IRS2 interaction or the 2 genes-environment interaction could increase T2DM susceptibility, especially in East Asian population. This study aimed to investigate the impact of FOXO1-IRS2 interaction or the 2 genes-environment interactions on T2DM susceptibility in Chinese Han population in Chongqing.

MATERIALS AND METHODS

Subjects

We collected a total of 780 cases from a Chinese Han population in Chongqing, including 474 cases of T2DM and 306 cases of healthy controls. There was no blood relationship among these individuals. Cases of healthy control group were subjected to oral glucose tolerance testing (OGTT) to exclude diabetes and impaired glucose regulation (IGR). Diabetes and IGR were diagnosed according to 1999 World Health Organization (WHO) standard. Criteria for the diagnosis of diabetes was fasting plasma glucose (FPG) ≥ 126 mg/dL (7.0 mmol/L) or the 2-hour plasma glucose (2-hr PG) ≥ 200 mg/dL (11.1 mmol/L) during an OGTT. Criteria for the diagnosis of IGR was FPG ≥ 110 mg/dL (6.1 mmol/L) or the 2-hr PG ≥ 140 mg/dL (7.8 mmol/L) during an OGTT. Exclusion criteria included: 1) patients with a history of ketoacidosis; 2) patients requiring continuous insulin treatment; 3) patients having exocrine pancreatic disease; or 4) patients with exceptionally early age of disease onset (< 30 years). Moreover, the obesity was defined by waist circumference of 90 and 80 cm for men and women respectively according to the definition of International Diabetes Federation (IDF). This study was approved by the ethics committee of our institute, and all the subjects signed the informed consent.

Clinical measurements

All participants underwent a demographic, anthropometric, and biochemical investigation as described previously (20). Briefly, anthropometric parameters such as height, weight, blood pressure, waist and hip circumference were measured. Body mass index (BMI) and the waist-to-hip ratio (WHR) were calculated. Subjects were conducted with 75 g oral glucose tolerance test. Fasting blood samples were taken at 2 hours after serving sugar, and then DNA isolation and biochemical indices detection were performed. Glucose oxidase method was used to detect blood glucose, radioimmunoassay was used to detect insulin, and automatic biochemical analyzer was used to detect blood lipids. Homeostasis model assessment (HOMA) calculated by FPG and insulin was used for estimating IR index and β-cell function (21).

Selection of SNPs loci

According to public SNP database of National Center for Biotechnology Information (NCBI), we selected 15 SNPs from missense mutation of exon region, promoter, introns and a 3'-none-coding region with high degree of heterozygosity (heterozygosity > 0.1) in FOXO1 gene. In the 105 cases of Chinese Han population in Chongqing (48 healthy cases, 57 T2DM cases), we analyzed the 15 SNPs by polymerase chain reaction-based restriction fragment length polymorphism (PCR-RFLP) and gene sequencing. Polymorphism was not found with the reported SNPs from the 3 databases (including 2 missense mutations). Finally, we conducted association analysis of the 5 loci from the screened genes of the 12 SNPs that were coincided with Hardy-Weinberg equilibrium (48 cases of healthy individuals and 57 cases of T2DM patients, P > 0.05) (Table 1), and they were rs17592236, rs9577066, rs7986407, rs4581585, and rs4325426, respectively (standardized LD coefficient [D'] > 0.7). Then, linkage disequilibrium (LD) analysis was performed again with much more cases (204 cases of healthy individuals and 414 cases of T2DM patients) to confirm the 5 selected loci (Fig. 1). Furthermore, according to literatures (222324), Gly1057Asp and −769C/T loci in IRS2 gene were selected to conduct the further detection.
Table 1

LD coefficients (D' and r2) among 12 SNPs in FOXO1

jkms-32-264-i001
D
r2 rs17592236 rs9577066 rs2755209 rs2721069 rs6563840 rs7986407 rs7335520 rs4581585 rs4325426 rs9549244 rs2297627 rs9577097
rs17592236 - 1.0000 1.0000 0.7908 1.0000 1.0000 0.8811 1.0000 0.8399 1.0000 1.0000 1.0000
rs9577066 0.0149 - 0.6705 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
rs2755209 0.1968 0.0341 - 1.0000 1.0000 1.0000 1.0000 1.0000 0.9657 1.0000 1.0000 0.6705
rs2721069 0.1613 0.0578 0.7628 - 0.9692 1.0000 1.0000 0.9692 0.8554 1.0000 1.0000 1.0000
rs6563840 0.2296 0.0650 0.8568 0.8363 - 1.0000 0.9377 0.9702 0.9389 1.0000 1.0000 1.0000
rs7986407 0.2229 0.0669 0.8829 0.8640 0.9704 - 0.9382 1.0000 0.9688 1.0000 1.0000 1.0000
rs7335520 0.1678 0.0690 0.9101 0.8381 0.8277 0.8538 - 0.9377 0.9042 0.9382 0.9382 1.0000
rs4581585 0.2296 0.0650 0.8568 0.8363 0.9413 0.9704 0.8277 - 0.9389 1.0000 1.0000 1.0000
rs4325426 0.1768 0.0595 0.7320 0.7112 0.8076 0.8343 0.7050 0.8076 - 0.9688 0.9688 1.0000
rs9549244 0.2229 0.0669 0.8829 0.8640 0.9704 1.0000 0.8538 0.9704 0.8343 - 1.0000 1.0000
rs2297627 0.2229 0.0669 0.8829 0.8640 0.9704 1.0000 0.8538 0.9704 0.8343 1.0000 - 1.0000
rs9577097 0.0149 1.0000 0.0341 0.0578 0.0650 0.0669 0.0690 0.0650 0.0595 0.0669 0.0669 -
LD = linkage disequilibrium, D' = standardized LD coefficient, SNP = single nucleotide polymorphism, FOXO1 = forkhead box O1.
Fig. 1
LD maps for SNPs genotyped in FOXO 1 gene.
Shades of gray show the strength of the pairwise LD based on D', and numbers indicate the value of r 2 expressed as a percentage.
LD = linkage disequilibrium, SNP = single nucleotide polymorphism, FOXO1 = forkhead box O1, D' = standardized LD coefficient.
jkms-32-264-g001

Genotype identification of SNPs loci

Phenol-chloroform method was used to extract genomic DNA. The primers were designed and synthesized based on public SNP databases and relevant literatures. SNPs genotypes were identified by PCR-RFLP. The PCR products with different enzyme genotype were randomly selected and confirmed by sequencing. PCR primer, reaction condition, and restrictive enzyme of each SNP locus are listed in Table 2.
Table 2

PCR primers, reaction condition, and restriction enzymes

jkms-32-264-i002
Gene names SNPs PCR product length, bp Primers PCR Tm, ℃ Restriction enzymes Restriction enzyme digested fragments, bp
FOXO1 rs4325426 170 5'TGCTAGGATATGTAGCAGATGCAACCAGGC 3' 61.0 BsuRI 170, 140, 30
5'TGGTGACAGGCATGAGAG ATACCTTTTTGG 3'
rs4581585 151 5'TTTGAGCATCATGTTGCCACTCAAGAAGTT 3' 62.5 DraI 151, 125, 26
5'CGAAGCCCACAACCCACT GAGCATTT 3'
rs7986407 325 5'CTACTCGATAGCTCCCTACTCT 3' 54.0 Hin6I 325, 223, 102
5'ACGTACTGCTGGCAACTGACT 3'
rs9577066 229 5'GTACGCATAGTTCAGTAGAAG 3' 56.5 MvaI 229, 121, 108
5'TATACATGCTTAACTGGTTTG 3'
rs17592236 238 5'GGGAGAATGAGATGAAGTAT 3' 49.0 Eco47I 238, 164, 74
5'GACCTCTGTAGTCCTGGGAG 3'
IRS2 Gly1057Asp 257 5'CAAAAGCCATCTCGGTGTAGT 3' 56.5 HhaI 257, 195, 62, 54
5'GCTCTCCGACTACATGAACCTC 3'
−769C/T 308 5'CTCTTCCGCGCCCCTTTTCCC 3' 59.0 Btg I 308, 198, 110
PCR = polymerase chain reaction, SNP = single nucleotide polymorphism, Tm = melting temperature, FOXO1 = forkhead box O1, IRS2 = insulin receptor substrate 2.

Statistical analysis

Statistical analysis was performed by SPSS 11.5 software (SPSS Inc, Chicago, IL, USA): categorical data was represented by rate (%); continuous data was represented by mean plus or minus standard deviation (± s); comparison of clinical data between the 2 groups was conducted by t-test; comparison between allele and genotypic distribution was conducted by χ2 test; relative T2DM risk of with SNP genotypes was analyzed by non-conditional logistic regression analysis, which was represented by odds ratio (OR) and 95% confidence interval (CI), and the relative risk OR was used to adjust age and gender. Linkage Disequilibrium Analyzer (LDA) 1.0 software (Chinese National Human Genome Center, Beijing, China) was used to perform Hardy-Weinberg equilibrium and LD analysis. Multifactor dimensionality reduction (MDR) was used for analysis of the interaction between genes or gene and environment by MDR_3.0.2 software (Computational Genetics Laboratory of the University of Pennsylvania, Philadelphia, PA, USA). This method effectively reduced high-dimensional data into 2-dimensional data and established multi-factor interaction model with different combinations. The model with highest cross-validation consistency (CVC) and the maximum testing balanced accuracy (TBA) was considered as a best model. If interactive model had statistical significance (P < 0.05), then MDRpt-1.0_beta_2 software was used to conduct permutation test, which further confirmed the reliability of the result.

Ethics statement

This study was approved by the Institutional Review Board of the First Affiliated Hospital of Chongqing Medical University (2013044). Informed consent was signed and submitted by every subject.

RESULTS

Clinical characteristics

Clinical features of the subjects are shown in Table 3. There was no significant difference in age, sex, diastolic blood pressure (DBP), total cholesterol, and low-density lipoprotein cholesterol between the 2 groups. However, T2DM patients had significantly higher BMI, WHR, waist, systolic blood pressure (SBP), triglycerides (TG), fasting glucose, 2-hour blood glucose after OGTT, fasting insulin, and 2-hour insulin after oral glucose compared with control subjects. In contrast, high-density lipoprotein cholesterol and homa-β were lower in the T2DM group compared to controls. As expected, T2DM patients had IR compared to controls, as evaluated by homeostasis model assessment of insulin resistance index (HOMA-IR).
Table 3

Clinical characteristics of the subjects

jkms-32-264-i003
Characteristics Control group T2DM group
Cases (male/female) 306 (179/127) 474 (258/216)
Age, yr 48.71 ± 10.89 50.54 ± 12.97
BMI, kg/m2 22.85 ± 2.81 24.67 ± 3.37*
WHR 0.84 ± 0.07 0.90 ± 0.07*
Waist, cm 77.50 ± 8.84 84.85 ± 9.40*
SBP, mmHg 111.55 ± 16.41 127.04 ± 20.05*
DBP, mmHg 73.01 ± 9.11 76.90 ± 10.65
TG, mmol/L 1.23 ± 0.59 2.03 ± 1.57*
Total cholesterol, mmol/L 5.05 ± 1.02 5.09 ± 1.24
High-density lipoprotein cholesterol, mmol/L 1.54 ± 0.39 1.38 ± 0.48*
Low-density lipoprotein cholesterol, mmol/L 2.75 ± 0.84 2.73 ± 0.97
FPG, mmol/L 5.06 ± 0.48 8.44 ± 3.33*
2-hour blood glucose after oral glucose tolerance test, mmol/L 6.01 ± 0.97 17.82 ± 6.66*
Fasting insulin, mU/L 11.54 ± 5.88 15.83 ± 11.19*
2-hour insulin after oral glucose tolerance test, mU/L 48.90 ± 28.58 52.71 ± 41.09*
HOMR-IR 2.60 ± 1.35 5.83 ± 4.98*
HOMR-β 163.00 ± 112.14 102.25 ± 139.08*
HOMA-IR = (FPG × Fasting insulin)/22.5; HOMA-β = 20 × Fasting insulin/(FPG − 3.5). After adjusting of age, BMI, WHR, and waist were compared; after adjusting age, BMI, and WHR, the other indicators were compared.
T2DM = type 2 diabetes mellitus, BMI = body mass index, WHR = waist-to-hip ratio, SBP = systolic blood pressure, DBP = diastolic blood pressure, TG = triglycerides, FPG = fasting plasma glucose, HOMA-IR = homeostasis model assessment of insulin resistance index, HOMA-β = homeostasis model assessment of β-cell function.
Compared with normal control: *P<0.01.

Rs7986407 and rs4581585 in FOXO1 gene and IRS2 gene are associated with T2DM

A allele frequency (31.2% vs. 24.2%; χ2 = 7.798; P = 0.005) of rs7986407 locus and AA genotype (12.0% vs. 7.4%; χ2 = 7.070; P = 0.029) in T2DM group (414 cases), as well as CC genotype frequency of rs4581585 locus were significantly higher than those in control group (290 cases). However, there was no significant difference in other loci between the 2 groups (P > 0.05). Association analysis showed that the T2MD risk of individuals with rs7986407 AA genotypes was 2.092 times higher than that with GG genotype, and the T2MD risk of individuals with rs4581585 CC genotypes was 2.57 times higher than those with TT genotype. Taken together, it indicates that both rs7986407 locus and rs4581585 locus are closely related to the risk of T2DM (Table 4).
Table 4

Association between polymorphisms in FOXO1 and the risk of T2DM

jkms-32-264-i004
SNPs Genotypes Control (n = 290) T2DM (n = 414) Additive Dominant Recessive
OR (95% CI) P OR (95% CI) P OR (95% CI) P
rs17592236 TT 39 (14.9) 59 (15.1) - - - - - -
CT 121 (46.2) 170 (43.5) 1.088 (0.676–1.759) 0.729 1.085 (0.654–1.600) 0.922 0.896 (0.647–1.241) 0.508
CC 102 (38.9) 162 (41.4) - - - - - -
rs9577066 AA 0 0 - - - - - -
AC 8 (3.1) 10 (2.7) - - - - 0.811 (0.312–2.110) 0.668
CC 252 (96.9) 364 (97.3) - - - - - -
rs7986407 AA 20 (7.4) 47 (12.0) - - - - - -
AG 91 (33.6) 150 (38.4) 2.092 (1.178–3.731) 0.012 1.789 (1.024–3.125) 0.041 1.585 (1.149–2.188) 0.005
GG 160 (59.0) 194 (49.6) - - - - - -
rs4581585 CC 17 (6.5) 51 (14.1) - - - - - -
CT 132 (50.6) 157 (43.5) 2.571 (1.404–4.695) 0.002 2.457 (1.374–4.405) 0.002 1.071 (0.770–1.488) 0.686
TT 112 (42.9) 153 (42.4) - - - - - -
rs4325426 CC 16 (6.2) 25 (6.7) - - - - - -
AC 121 (46.5) 147 (39.2) 1.272 (0.642–2.519) 0.491 1.075 (0.556–2.083) 0.828 0.758 (0.549–1.046) 0.092
AA 123 (47.3) 203 (54.1) - - - - - -
Genotype distributions were shown with number (%). ORs, 95% CI, and P values were from logistic regression analyses with additive, dominant, and recessive models controlling for age and sex as covariates. In additive models, ORs was expressed per difference in number of rare allele. Genotype was given codes of (0, 1, and 2), (0, 1, and 1), (0, 0, and 1) in additive, dominant and recessive models, respectively.
FOXO1 = forkhead box O1, T2DM = type 2 diabetes mellitus, SNP = single nucleotide polymorphism, OR = odds ratio, CI = confidence interval.

Gly1057Asp locus in IRS2 gene is associated with T2DM

A allele and AA genotypic frequency of Gly1057Asp locus in T2DM group (474 cases) were significantly lower than those in control group (306 cases), 32.8% vs. 38.5% (χ2 = 5.383; P = 0.02) and 10.5% vs. 13.4% (χ2 = 6.743; P = 0.034), respectively (Table 5). The T2DM risk of individuals with GG genotype of Gly1057Asp was 1.42 times higher compared to AA genotype (additive model: OR = 1.422; 95% CI = 1.037–1.949; P = 0.029), indicating a significant association of GG genotype of Gly1057Asp with T2DM. Furthermore, allele and genotypic frequency distribution of −769C/T locus had neither a significant difference between the 2 groups, nor an apparent association with T2DM (P > 0.05) (Table 5).
Table 5

Associations between polymorphisms in IRS2 and the risk of T2DM

jkms-32-264-i005
SNPs Genotypes T2DM Additive Dominant Recessive
OR (95% CI) P OR (95% CI) P OR (95% CI) P
Gly1057Asp n = 305 n = 467
 AA 41 (13.4) 49 (10.5) 1.422 (1.037–1.949) 0.029 1.461 (1.081–1.973) 0.014 1.339 (0.851–2.106) 0.207
 GA 153 (50.2) 208 (44.5) - - - - - -
 GG 111 (36.4) 210 (45.0) - - - - - -
−769C/T n = 293 n = 455
 CC 54 (18.4) 72 (15.8) 1.216 (0.870–1.700) 0.252 1.243 (0.906–1.705) 0.177 1.182 (0.793–1.761) 0.412
 CT 143 (48.8) 217 (47.7) - - - - - -
 TT 96 (32.8) 166 (36.5) - - - - - -
Genotype distributions are shown as number (%). ORs, 95% CI, and P values were from logistic regression analyses with additive, dominant, and recessive models controlling for age and sex as covariates. In additive models, ORs are expressed per difference in number of rare allele. Genotype were given codes of (0, 1, and 2), (0, 1, and 1), (0, 0, and 1) in additive, dominant and recessive models, respectively.
IRS2 = insulin receptor substrate 2, T2DM = type 2 diabetes mellitus, SNP = single nucleotide polymorphism, OR = odds ratio, CI = confidence interval.

T2DM susceptibility is increased due to the interaction between rs7986407 or rs4325426 in FOXO1 gene and obesity

Analysis of inter-gene interaction was conducted by MDR method for the 7 SNPs at FOXO1 and IRS2, showing no interaction model for T2DM (P > 0.05) (data not shown). It indicates that the FOXO1-IRS2 gene interaction is irrelevant to the prevalence of T2DM
Next, we analyzed the gene-environment interaction between the 7 SNPs and environment factors (aging, obesity, and hypertriglyceridemia). MDR analysis suggested that the interaction model composed of the 3 factors, rs7986407, rs4325426, and waist, was the best model with statistical significance (P value from sign test [Psign] = 0.0251; P value from permutation test [Pperm] = 0.000–0.001) (Table 6 and Fig. 2). Moreover, case/control ratio of more than 1.23 was defined as high-risk, and the ratio of less than 1.23 was defined as low-risk. T2DM incidence in high-risk group was increased significantly compared to low-risk group, and its OR was 3.976 (95% CI = 1.156–13.675), indicating that rs7986407, rs4325426, and waist are associated with T2DM. In case of non-obesity, the risk of T2DM was increased in individuals with AA of rs7986407 combined with CC or AC of rs4325426, as well as AG of rs7986407 combined AA of rs4325426; while in obesity, individuals with AA of rs7986407 and any genotype of rs4325426, or individuals with AA of rs4325426 and any genotype of rs7986407 were prone to T2DM.
Table 6

MDR models of gene-environmental interactions on T2DM

jkms-32-264-i006
Models Testing Bal. Acc. CVC OR (95% CI) Psign Pperm
Waist 0.6572 10/10 4.007 (1.131–14.204) 0.0274 0.000–0.001
rs7986407, waist 0.6584 9/10 3.819 (1.123–12.994) 0.0285 0.000–0.001
rs7986407, rs4325426, waist 0.6616 8/10 3.976 (1.156–13.675) 0.0251 0.000–0.001
MDR = multifactor dimensionality reduction, T2DM = type 2 diabetes mellitus, CVC = cross-validation consistency, OR = odds ratio, CI = confidence interval, Psign = P value from sign test, Pperm = P value from permutation test.
Fig. 2
Distribution of high-risk and low-risk genotypes in the best 3-factor model.
High-risk genotypes are in dark gray and low-risk genotypes are in light gray. The number of T2DM patients (left bar in boxes) and control subjects (right bar in boxes) are shown for each 3-factor combination. rs7986407: 0 for AA genotype, 1 for AG genotype, 2 for GG genotype; rs4325426: 0 for AA genotype, 1 for AC genotype, 2 for CC genotype; Age: 0 for ≤ 40 years, 1 for > 40 years and ≤ 60 years, 2 for > 60 years; BMI: 0 for < 25 kg/m2, 1 for ≥ 25 kg/m2; Waist: 0 for < 90 cm (male) or 80 cm (female), 1 for ≥ 90 cm (male) or 80 cm (female); TG: 0 for < 1.7 mmol/L, 1 for ≥ 1.7 mmol/L.
T2DM = type 2 diabetes mellitus, BMI = body mass index, TG = triglycerides.
jkms-32-264-g002

DISCUSSION

In this study, we showed that there were significant differences in genotype distributions of rs7986407 and rs4581585 in FOXO1 gene between control group and T2DM group, suggesting that FOXO1 gene mutation is associated with T2DM. AA genotype for rs7986407 and CC genotype for rs4581585 increase the risk of T2DM in Chinese Han population in Chongqing. The association of FOXO1 gene mutation with T2DM was not found in Caucasians and African-Americans in the past (17), but the study of a population in Germany and Finland suggested that FOXO1 gene mutation was associated with the increasing risk of T2DM. Thus, one possibility could be that the inconsistence of these findings may be associated with ethnic differences. However, the study on the Chinese Han population by Li and colleagues (19) suggested that SNPs at FOXO1 were not associated with T2DM. The selected SNPs in our study almost covered the whole FOXO1 gene. Among the SNPs, the rs17592236 was the only one detected by Li et al. (19). Therefore, the inconsistent results between our study and the Li and colleagues' study may be related to the different selection of SNPs. Previously, it has been reported in Asian population (including other population besides Chinese population) that IRS2 gene was associated with T2DM (1516). Although there was obvious ethnic difference, in line with the Asian population study, we found that IRS2 gene mutation was also associated with T2DM, and GG genotype for Gly1057Asp locus increased the risk of T2DM in Chinese Han population in Chongqing.
Based on the above findings, polymorphism of FOXO1 and IRS2 may be both associated with the risk of T2DM, we asked ourselves whether there is an interaction between them, which would increase susceptibility to T2DM. To answer the question, we analyzed 7 SNPs at FOXO1 and IRS2 by MDR, suggesting that there is no association of T2DM prevalence with FOXO1-IRS2 genes interaction. Although FOXO1 and IRS2 were associated and interacted with each other by pathophysiology and signaling transduction (111314), they may not be related to each other through molecular biology directly; and the interaction between them had no effect on prevalence of T2DM.
Genetic susceptibility is an important factor to T2DM, but the result of genome-wide association studies (GWAS) can only explain about 10% of the risk of T2DM susceptibility. This is because the prevalence of T2DM is also related to environmental factors, and the interaction between genetic factors and environment also has an important influence on individual illness (62526). Therefore, it is more beneficial to identify the T2DM susceptible population by studying the gene-environment interaction. Aging, high-fat diet, lack of exercise, and obesity are environmental factors that affect the prevalence of T2DM. In this study, we analyzed the interaction between SNPs of either FOXO1 or IRS2 and aging, obesity, and high TG hyperlipidemia. It suggested that the risk of T2DM prevalence can be increased by FOXO1-obesity interaction.
In this study, we analyzed the interactions between either rs7986407 locus or rs4325426 locus in FOXO1 and obesity. In the case of non-obesity, T2DM susceptibility was increased in individuals with AA of rs7986407 locus in combination with heterozygous rs4325426 or individuals with AA of rs4325426 locus in combination with heterozygous rs7986407; while in obesity, if obesity-related environment factors existed, T2DM susceptibility would be increased in individuals with AA genotype for rs7986407 or with AA genotype for rs4325426. These data suggest that the FOXO1-obesity interaction increases susceptibility to T2DM. Interestingly, the result of correlation analysis was inconsistent with the interactive loci of FOXO1 gene and obesity to the prevalence of T2DM. Correlation analysis suggested that there was no association between obesity and rs4581585 locus that may increase risk of T2DM; while rs4325426 locus was not related to the risk of T2DM, but associated with obesity in T2DM susceptibility. Thus, it can be seen that the gene-gene interaction and gene-environment interaction are complex, and genetic impact on incidence of individual would change dependent on environmental factors (27). The risk of individual susceptibility may not be truly and fully reflected by only studying the genetic influence on the prevalence of T2DM (67). Therefore, the study in the influence of gene-gene interaction and gene-environment interaction on T2DM susceptibility would be possible to find real high-risk susceptible populations.
Our study also has limitations. One of the limitations in this study was that the numbers of SNPs were not sufficient to capture most genetic information of the candidate genes. Another limitation was that we used MDR method to analyze gene-gene interaction and gene-environment interaction. However, the MDR method can only be used in association studies for qualitative traits. Future studies should elaborate on our findings by including more haplotype tagging or functional SNPs in the studied candidate genes by multiple analysis methods.
In summary, this study suggests that FOXO1 and IRS2 genes are T2DM susceptibility genes in Chinese Han population in Chongqing; obesity with genotypes of rs7986407 locus AA or rs4325426 locus AA in FOXO1 gene might be a high-risk factor of T2DM. Different interactions between T2DM susceptibility genes and environment factors would lead to different rate of T2DM susceptibility between individuals. Therefore, by studying the influence of gene-environment interaction on T2DM susceptibility, finding high-risk T2DM population, and focusing prevention and intervention on them, it may be an effective strategy to reduce the prevalence of T2DM.

Notes

Funding This study was funded by National Key Clinical Specialist Construction Program of China ([2011]170).

DISCLOSURE The authors have no potential conflicts of interest to disclose.

Author Contributions

  • Conceptualization: Gong L.

  • Data curation: Gong L, Li R, Ren W, Zhang S.

  • Investigation: Gong L, Wang Z, Yang M.

  • Writing - review & editing: Gong L.

References

1. Edelman SV. Type II diabetes mellitus. Adv Intern Med. 1998; 43:449–500.
2. Voight BF, Scott LJ, Steinthorsdottir V, Morris AP, Dina C, Welch RP, Zeggini E, Huth C, Aulchenko YS, Thorleifsson G, et al. Twelve type 2 diabetes susceptibility loci identified through large-scale association analysis. Nat Genet. 2010; 42:579–589.
3. O’Rahilly S, Barroso I, Wareham NJ. Genetic factors in type 2 diabetes: the end of the beginning? Science. 2005; 307:370–373.
4. Wild S, Roglic G, Green A, Sicree R, King H. Global prevalence of diabetes: estimates for the year 2000 and projections for 2030. Diabetes Care. 2004; 27:1047–1053.
5. Chan JC, Malik V, Jia W, Kadowaki T, Yajnik CS, Yoon KH, Hu FB. Diabetes in Asia: epidemiology, risk factors, and pathophysiology. JAMA. 2009; 301:2129–2140.
6. Bouret S, Levin BE, Ozanne SE. Gene-environment interactions controlling energy and glucose homeostasis and the developmental origins of obesity. Physiol Rev. 2015; 95:47–82.
7. Cornelis MC, Hu FB. Gene-environment interactions in the development of type 2 diabetes: recent progress and continuing challenges. Annu Rev Nutr. 2012; 32:245–259.
8. Stumvoll M, Goldstein BJ, van Haeften TW. Type 2 diabetes: principles of pathogenesis and therapy. Lancet. 2005; 365:1333–1346.
9. Haeusler RA, Hartil K, Vaitheesvaran B, Arrieta-Cruz I, Knight CM, Cook JR, Kammoun HL, Febbraio MA, Gutierrez-Juarez R, Kurland IJ, et al. Integrated control of hepatic lipogenesis versus glucose production requires FoxO transcription factors. Nat Commun. 2014; 5:5190.
10. Kitamura T. The role of FOXO1 in beta-cell failure and type 2 diabetes mellitus. Nat Rev Endocrinol. 2013; 9:615–623.
11. Guo S. Insulin signaling, resistance, and the metabolic syndrome: insights from mouse models into disease mechanisms. J Endocrinol. 2014; 220:T1–23.
12. Oliveira JM, Rebuffat SA, Gasa R, Gomis R. Targeting type 2 diabetes: lessons from a knockout model of insulin receptor substrate 2. Can J Physiol Pharmacol. 2014; 92:613–620.
13. Nakae J, Biggs WH 3rd, Kitamura T, Cavenee WK, Wright CV, Arden KC, Accili D. Regulation of insulin action and pancreatic beta-cell function by mutated alleles of the gene encoding forkhead transcription factor Foxo1. Nat Genet. 2002; 32:245–253.
14. Kitamura T, Nakae J, Kitamura Y, Kido Y, Biggs WH 3rd, Wright CV, White MF, Arden KC, Accili D. The forkhead transcription factor Foxo1 links insulin signaling to Pdx1 regulation of pancreatic beta cell growth. J Clin Invest. 2002; 110:1839–1847.
15. Haghani K, Bakhtiyari S. The study on the relationship between IRS-1 Gly972Arg and IRS-2 Gly1057Asp polymorphisms and type 2 diabetes in the Kurdish ethnic group in West Iran. Genet Test Mol Biomarkers. 2012; 16:1270–1276.
16. Okazawa K, Yoshimasa Y, Miyamoto Y, Takahashi-Yasuno A, Miyawaki T, Masuzaki H, Hayashi T, Hosoda K, Inoue G, Nakao K. The haplotypes of the IRS-2 gene affect insulin sensitivity in Japanese patients with type 2 diabetes. Diabetes Res Clin Pract. 2005; 68:39–48.
17. Karim MA, Craig RL, Wang X, Hale TC, Elbein SC. Analysis of FOXO1A as a candidate gene for type 2 diabetes. Mol Genet Metab. 2006; 88:171–177.
18. Müssig K, Staiger H, Machicao F, Stancáková A, Kuusisto J, Laakso M, Thamer C, Machann J, Schick F, Claussen CD, et al. Association of common genetic variation in the FOXO1 gene with beta-cell dysfunction, impaired glucose tolerance, and type 2 diabetes. J Clin Endocrinol Metab. 2009; 94:1353–1360.
19. Li T, Wu X, Zhu X, Li J, Pan L, Li P, Xin Z, Liu Y. Association analyses between the genetic polymorphisms of HNF4A and FOXO1 genes and Chinese Han patients with type 2 diabetes. Mol Cell Biochem. 2011; 353:259–265.
20. Jia WP, Pang C, Chen L, Bao YQ, Lu JX, Lu HJ, Tang JL, Wu YM, Zuo YH, Jiang SY, et al. Epidemiological characteristics of diabetes mellitus and impaired glucose regulation in a Chinese adult population: the Shanghai diabetes studies, a cross-sectional 3-year follow-up study in Shanghai urban communities. Diabetologia. 2007; 50:286–292.
21. Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, Turner RC. Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia. 1985; 28:412–419.
22. Lautier C, El Mkadem SA, Renard E, Brun JF, Gris JC, Bringer J, Grigorescu F. Complex haplotypes of IRS2 gene are associated with severe obesity and reveal heterogeneity in the effect of Gly1057Asp mutation. Hum Genet. 2003; 113:34–43.
23. Bodhini D, Radha V, Deepa R, Ghosh S, Majumder PP, Rao MR, Mohan V. The G1057D polymorphism of IRS-2 gene and its relationship with obesity in conferring susceptibility to type 2 diabetes in Asian Indians. Int J Obes. 2007; 31:97–102.
24. Villuendas G, Botella-Carretero JI, Roldán B, Sancho J, Escobar-Morreale HF, San Millán JL. Polymorphisms in the insulin receptor substrate-1 (IRS-1) gene and the insulin receptor substrate-2 (IRS-2) gene influence glucose homeostasis and body mass index in women with polycystic ovary syndrome and non-hyperandrogenic controls. Hum Reprod. 2005; 20:3184–3191.
25. Murea M, Ma L, Freedman BI. Genetic and environmental factors associated with type 2 diabetes and diabetic vascular complications. Rev Diabet Stud. 2012; 9:6–22.
26. Burguera B, Agusti A, Arner P, Baltasar A, Barbe F, Barcelo A, Breton I, Cabanes T, Casanueva FF, Couce ME, et al. Critical assessment of the current guidelines for the management and treatment of morbidly obese patients. J Endocrinol Invest. 2007; 30:844–852.
27. Romao I, Roth J. Genetic and environmental interactions in obesity and type 2 diabetes. J Am Diet Assoc. 2008; 108:S24–8.
TOOLS
ORCID iDs

Lilin Gong
https://orcid.org/http://orcid.org/0000-0003-4110-1935

Rong Li
https://orcid.org/http://orcid.org/0000-0001-6365-5782

Wei Ren
https://orcid.org/http://orcid.org/0000-0003-4222-8217

Zengchan Wang
https://orcid.org/http://orcid.org/0000-0002-0525-5315

Zhihong Wang
https://orcid.org/http://orcid.org/0000-0002-8259-9947

Maosheng Yang
https://orcid.org/http://orcid.org/0000-0003-3059-4690

Suhua Zhang
https://orcid.org/http://orcid.org/0000-0003-2033-473X

Similar articles