R1467H Variants of Rho Guanine Nucleotide Exchange Factor 11 (ARHGEF11) are Associated with Type 2 Diabetes Mellitus in Koreans

Article information

Diabetes Metab J. 2010;34(6):368-373
Publication date (electronic) : 2010 December 31
doi : https://doi.org/10.4093/kdj.2010.34.6.368
1Diabetes Clinical Research Center, Inha University Hospital, Incheon, Korea.
2Department of Internal Medicine, Inha University College of Medicine, Incheon, Korea.
3Department of Preventive and Social Medicine, Inha University College of Medicine, Incheon, Korea.
Corresponding author: Moonsuk Nam. Department of Internal Medicine, Inha University Hospital, Inha University College of Medicine, 7-206 Sinheung-dong 3-ga, Jung-gu, Incheon 400-711, Korea. namms@inha.ac.kr
Received 2010 May 20; Accepted 2010 September 28.

Abstract

Background

The human Rho guanine nucleotide exchange factor 11 (ARHGEF11) functions as an activator of Rho GTPases and is thought to influence insulin signaling. The R1467H variant of ARHGEF11 has been reported to be associated with susceptibility to type 2 diabetes mellitus (T2DM) in Western populations.

Methods

We investigated the effects of the R1467H variant on susceptibility to T2DM as well as related traits in a Korean population. We genotyped the R1467H (rs945508) of ARHGEF11 in 689 unrelated T2DM patients and 249 non-diabetic individuals and compared the clinical and biochemical characteristics according to different alleles.

Results

The H allele was significantly more frequent in T2DM cases than in controls (P = 0.037, 17.1% and 13.1%; respectively). H homozygocity was associated with a higher risk of T2DM compared to those with R/R or R/H genotype (odds ratio, 5.24; 95% confidence interval, 1.06 to 25.83; P = 0.042). The fasting plasma glucose, HbA1c, fasting insulin, HOMA2-IR and HOMA2-%β levels did not differ significantly between different genotypes.

Conclusion

Our study replicated associations of the ARHGEF11 polymorphism with increased risk of T2DM in a Korean population and thus supports previous data implicating a potential role of ARHGEF11 in the etiology of T2DM. Further studies revealing the underlying mechanism for this association are needed.

INTRODUCTION

Type 2 diabetes mellitus (T2DM) is a multi-factorial trait, in which individual risk is defined by the complex interplay of genetic and environmental factors [1-3]. To search for genetic loci that contribute to T2DM, several genomic linkage scans have been completed worldwide [4,5]. Several previous studies have reported evidence for linkage between T2DM and the chromosome 1q21-q24 region [6-8] and have also reported that this region contributes to impaired glucose tolerance (IGT) [9]. In the Pima Indians and Amish populations, investigators genotyped multiple single nucleotide polymorphisms (SNPs) within this region, and results pointed to a region including the Rho guanine nucleotide exchange factor 11 (ARHGEF11) gene [10,11]. ARHGEF11 is an activator of Rho GTPases and may play a role in the G protein signaling involved in β-cell apoptosis, insulin secretion and insulin signaling. ARHGEF11 is expressed in the pancreas, liver, muscle, and adipose tissue [12-14]. Previous studies have implicated Rho subfamily G proteins in physiological insulin secretion [13,14], and Rho GT-Pases have been suggested to play a role in insulin signaling through the activation of p38 mitogen activated protein kinase and Jun NH2-terminal kinase pathways [15-17]. Moreover, Rho family small GTP-binding protein TC10 activation and phosphatidylinositol 3-kinase activation regulates the dynamic actin rearrangement required for insulin-stimulated translocation of GLUT4 [18,19]. Together, these data suggest that proteins involved in G protein signaling, such as ARHGEF11, may play an important role in glucose homeostasis.

The full-length ARHGEF11 protein contains 1,522 amino acids, has a calculated molecular mass of 168.6 kD and includes 41 exons [20]. The C-terminal (amino acids 1,181-1,522) region of ARHGEF11 protein has been shown to regulate ARHGEF11 protein activity [21]. Recently, a non-synonymous amino acid substitution variant within the C-terminal region of ARHGEF11 R1467H (rs945508, G/A), was reported to be associated with T2DM in Western sample populations [11,22,23]. In the Amish, the 1467R allele was the risk allele for T2DM, while the 1467H allele was the risk allele for insulin resistance and T2DM in Pima Indians [11,22]. Moreover, replicated association of the 1467H allele as the risk allele for T2DM and/or IGT in a German Caucasian population have been reported [23]. This inconsistency in the effects of genetic variants may be due to different genetic backgrounds between the different ethnic groups. Therefore, in the present study, we evaluated the effects of this variant on the susceptibility to T2DM as well as its relation to metabolic parameters in a Korean population.

METHODS

Study subjects

A total of 689 unrelated T2DM patients and 249 non-diabetic individuals were recruited at the Inha University Hospital of Korea from April 2004 to December 2007. Diabetes was defined according to the 1997 American Diabetes Association diagnostic criteria [24]. Exclusion criteria were type 1 diabetes mellitus, current malignancy, history of ketoacidosis, positive for glutamic acid decarboxylase antibody, severe kidney or liver disease, and corticosteroid use. All subjects gave written informed consent before taking part in the study, which was approved by the Inha University Hospital Institutional Review Board.

Single nucleotide polymorphism genotyping

Genomic DNA was extracted from leukocytes in the whole-blood samples using a Gentra Puregene Blood Kit Plus (Qiagen, Valencia, CA, USA). Genotyping of the R1467H was achieved using the TaqMan allelic discrimination assay (Applied Biosystems, Foster City, CA, USA). To assess genotyping reproducibility, a random -10% selection of the samples was re-genotyped in all SNPs with 100% concordance. The SNPs were in accordance with Hardy-Weinberg equilibrium in the T2DM subjects (P > 0.10) and in the control subjects (P > 0.20).

Clinical and biochemical measurements

Baseline laboratory data and anthropometric measurements were collected for each patient. Low density lipoprotein cholesterol levels were calculated using the Friedewald formula [25], and plasma glucose was measured using a glucose oxidase method. Hemoglobin A1c (HbA1c) values were determined using high-performance liquid chromatography, and insulin concentrations were measured using a radioimmunoassay kit (IRMA kit; DAINABOT, Tokyo, Japan). The insulin resistance and insulin secretion indexes were calculated using the HOMA2 model (HOMA2-IR and HOMA2-%β, respectively), i.e., the updated computer model for the pairing of fasting plasma glucose and insulin (The HOMA2 model is available from www.ocdem.ox.ac.uk) [26].

Statistical analysis

Statistical analyses were performed using SPSS for Windows version 18.0 (SPSS Inc., Chicago, IL, USA), and data are shown as means ± standard deviation, unless otherwise indicated. Associations between categorical variables were analyzed using χ2 tests and continuous variables were assessed using Student's t-test. The differences in allelic and genotypic frequencies between the diabetic and control individuals were analyzed using Pearson's χ2 test. Fasting insulin, HOMA2-IR and HOMA2-%β variables were log-transformed to correct for skew. The odds ratio (OR) was assessed by counting the number of risk alleles for each individual, and association of genotypes with T2DM was assessed via a logistic regression analysis adjusted for age, gender, and body mass index (BMI). All P < 0.05 were considered to be significant.

RESULTS

We performed a case-control study including 249 non-diabetic individuals and 689 unrelated T2DM patients in Korea. The group of patients with T2DM included 383 men and 306 women, and the control group included 128 men and 121 women (P = 0.256). Non-diabetic individuals were older than T2DM patients (P < 0.001, 63.25 ± 7.12, 52.23 ± 9.70 years, respectively) and had a lower BMI 23.59 ± 2.89 kg/m2 (P < 0.001, 25.21 ± 3.87 kg/m2). The clinical features of the participants are summarized in Table 1. Allelic and genotypic frequencies of the SNP in ARHGEF11 are shown in Table 2. The H allele was significantly more frequent in T2DM cases than in controls (17.1% in T2DM vs. 13.1% in control group, P = 0.037). Genotypic frequencies were significantly different in T2DM patients compared to those of the controls (P = 0.040, R/R 479 [69.5%] vs. 186 [74.7%], R/H 184 [26.7%] vs. 61 [24.5%], H/H 26 [3.8%] vs. 2 [0.8%], respectively). Moreover, in patients with T2DM, the probability of having the H/H genotype vs. R/R + H/H was significantly higher compared to that of the control subjects (P = 0.018, 3.8% vs. 0.8%), and this association was significant after adjustment for age, gender and BMI (OR, 5.24; 95% confidence intervals [CI], 1.06-25.83; P = 0.042). The odds ratio estimates along with 95% CI values for the ARHGEF11 R1467H allele in various studies conducted in T2DM patients are presented in Table 3. In order to further examine the role of the H allele in T2DM subjects, we assessed whether the H homozygote was associated with insulin secretion or insulin resistance in T2DM subjects. There was no significant difference between R/R or R/H subjects and H/H subjects with regard to fasting plasma glucose, HbA1c, fasting insulin, HOMA2-IR or HOMA2-%β (P > 0.05, Table 4).

Table 1

Clinical characteristics of the study population

Table 2

Allelic and genotypic distributions of the R1467H variant in T2DM and control subjects

Table 3

Odds ratio estimates with confidence intervals of the R1467H variant in multiple ethnic populations

Table 4

Clinical and biochemical characteristics between R/R+R/H and H/H genotypes of ARHGEF11

DISCUSSION

In this study, we found a significant association between the variation at R1467H in ARHGEF11 and T2DM in a Korean population, similar to the results previously reported in Pima Indians and German Caucasians [10,23]. It is noteworthy that the allele frequency and odds ratio in our study are closer to those in Pima Indians (Table 3) and that the association between the variation at R1467H in ARHGEF11 and T2DM which has been observed in other ethnic groups has been replicated, supporting the reliabilities of both previous studies as well as that of our present study. However, the association with type 2 diabetes may have a stronger impact in the German Caucasians than in Koreans due to the higher frequency of the R1467H risk allele (H) in German Caucasians (46% vs. 13%). The lower allelic frequency in Koreans is similar to the data reported by the HapMap project [27]. A different association of the common R allele with T2DM was shown in a study of the Amish [11]. This may be due to population differences; however, considering the Swiss-German ancestry of the Amish, it was surprising that data from German Caucasians more closely adhered to the Pima findings rather than those reported in the Amish. Although, there is no evidence that the R1467H variant may exert an effect on ARHGEF11 expression, an Arg to His amino acid change in the C-terminal region of ARHGEF11 may be brought about by the variant. Thus, we believe that our present findings are consistent with those found in Pima Indians and German Caucasians, strongly suggesting the role of the variant in the pathophysiology of T2DM.

The H allele at R1467H in ARHGEF11 increases the risk for T2DM in Pima Indians by reducing insulin-stimulated glucose uptake. Among the non-diabetic, full-heritage Pima Indians, subjects with an H allele had a lower mean glucose disposal rate during a hyperinsulinemic-euglycemic clamp. To further illuminate the role of the H allele in a population of Korean T2DM patients, related traits were compared in subjects with H/H homozygocity and in those carrying the major R allele (R/H or R/R). Our study did not reveal any significant differences in fasting plasma glucose, HbA1c, fasting insulin, HOMA2-IR and HOMA2-%β levels in T2DM subjects between the two groups. The H allele also did not increase fasting plasma glucose in the control patients, possibly due to the various drug effects experienced by T2DM patients or that the HOMA parameters are a relatively crude analysis method compared to the clamp method [28]. It may also be due to different ethnic and environmental backgrounds. Therefore, additional investigations in independent populations, as well as functional studies, will be necessary to clarify the influence of the ARHGEF11 variant on T2DM. There is also a possibility that the R1467H variant may not be the true functional variant but rather a variant in linkage disequilibrium with the true functional variant.

This study had several limitations. First, our study was similar to others in Pima Indians and had similar associations between the variations at R1467H in ARHGEF11 and T2DM but could not determine a significant association between the R/R+R/H and H/H groups of T2DM with regard to fasting plasma glucose, HbA1c, fasting insulin, HOMA2-IR or HOMA2-%β level. In Pima Indians, the investigator demonstrated that the H allele at R1467H in ARHGEF11 increases the risk of T2DM by reducing insulin-stimulated glucose uptake in non-diabetic controls. However, our study did not measure insulin levels in non-diabetic control subjects. Second, we used fasting insulin concentration and the HOMA2 model as markers of insulin resistance. The HOMA2 model has been significantly correlated with the hyperinsulinemic-euglycemic clamp, which is accepted as the standard in defining insulin resistance; however, the HOMA2 model is described for approximate estimation [28]. Finally, because our study was conducted in a Korean population, we cannot generalize these findings to other ethnic groups.

In conclusion, our study replicated associations of the ARHGEF11 polymorphism with increased risk of T2DM in a Korean population and thus supports previous data implicating a potential role of ARHGEF11 in the etiology of T2DM. Further replication studies in different ethnic populations are necessary to confirm the effect of the R1467H variant on T2DM.

ACKNOWLEDGEMENT

This study was supported by a grant of the Korean Health 21 Research and Development Project, Ministry of Health & Welfare, Republic of Korea (A050463, A070001).

References

1. Stumvoll M, Goldstein BJ, van Haeften TW. Type 2 diabetes: principles of pathogenesis and therapy. Lancet 2005;365:1333–1346. 15823385.
2. Pezzolesi MG, Nam M, Nagase T, Klupa T, Dunn JS, Mlynarski WM, Rich SS, Warram JH, Krolewski AS. Examination of candidate chromosomal regions for type 2 diabetes reveals a susceptibility locus on human chromosome 8p23.1. Diabetes 2004;53:486–491. 14747302.
3. Rich SS. Mapping genes in diabetes. Genetic epidemiological perspective. Diabetes 1990;39:1315–1319. 2227105.
4. Sladek R, Rocheleau G, Rung J, Dina C, Shen L, Serre D, Boutin P, Vincent D, Belisle A, Hadjadj S, Balkau B, Heude B, Charpentier G, Hudson TJ, Montpetit A, Pshezhetsky AV, Prentki M, Posner BI, Balding DJ, Meyre D, Polychronakos C, Froguel P. A genome-wide association study identifies novel risk loci for type 2 diabetes. Nature 2007;445:881–885. 17293876.
5. Saxena R, Voight BF, Lyssenko V, Burtt NP, de Bakker PI, Chen H, Roix JJ, Kathiresan S, Hirschhorn JN, Daly MJ, Hughes TE, Groop L, Altshuler D, Almgren P, Florez JC, Meyer J, Ardlie K, Bengtsson Bostrom K, Isomaa B, Lettre G, Lindblad U, Lyon HN, Melander O, Newton-Cheh C, Nilsson P, Orho-Melander M, Rastam L, Speliotes EK, Taskinen MR, Tuomi T, Guiducci C, Berglund A, Carlson J, Gianniny L, Hackett R, Hall L, Holmkvist J, Laurila E, Sjogren M, Sterner M, Surti A, Svensson M, Tewhey R, Blumenstiel B, Parkin M, Defelice M, Barry R, Brodeur W, Camarata J, Chia N, Fava M, Gibbons J, Handsaker B, Healy C, Nguyen K, Gates C, Sougnez C, Gage D, Nizzari M, Gabriel SB, Chirn GW, Ma Q, Parikh H, Richardson D, Ricke D, Purcell S. Genome-wide association analysis identifies loci for type 2 diabetes and triglyceride levels. Science 2007;316:1331–1336. 17463246.
6. Hanson RL, Ehm MG, Pettitt DJ, Prochazka M, Thompson DB, Timberlake D, Foroud T, Kobes S, Baier L, Burns DK, Almasy L, Blangero J, Garvey WT, Bennett PH, Knowler WC. An autosomal genomic scan for loci linked to type II diabetes mellitus and body-mass index in Pima Indians. Am J Hum Genet 1998;63:1130–1138. 9758619.
7. Elbein SC, Hoffman MD, Teng K, Leppert MF, Hasstedt SJ. A genome-wide search for type 2 diabetes susceptibility genes in Utah Caucasians. Diabetes 1999;48:1175–1182. 10331426.
8. Vionnet N, Hani EH, Dupont S, Gallina S, Francke S, Dotte S, De Matos F, Durand E, Lepretre F, Lecoeur C, Gallina P, Zekiri L, Dina C, Froguel P. Genomewide search for type 2 diabetes-susceptibility genes in French whites: evidence for a novel susceptibility locus for early-onset diabetes on chromosome 3q27-qter and independent replication of a type 2-diabetes locus on chromosome 1q21-q24. Am J Hum Genet 2000;67:1470–1480. 11067779.
9. Xiang K, Wang Y, Zheng T, Jia W, Li J, Chen L, Shen K, Wu S, Lin X, Zhang G, Wang C, Wang S, Lu H, Fang Q, Shi Y, Zhang R, Xu J, Weng Q. Genome-wide search for type 2 diabetes/impaired glucose homeostasis susceptibility genes in the Chinese: significant linkage to chromosome 6q21-q23 and chromosome 1q21-q24. Diabetes 2004;53:228–234. 14693720.
10. Ma L, Hanson RL, Que LN, Cali AM, Fu M, Mack JL, Infante AM, Kobes S, Bogardus C, Shuldiner AR, Baier LJ. Variants in ARHGEF11, a candidate gene for the linkage to type 2 diabetes on chromosome 1q, are nominally associated with insulin resistance and type 2 diabetes in Pima Indians. Diabetes 2007;56:1454–1459. 17287471.
11. Fu M, Sabra MM, Damcott C, Pollin TI, Ma L, Ott S, Shelton JC, Shi X, Reinhart L, O'Connell J, Mitchell BD, Baier LJ, Shuldiner AR. Evidence that Rho guanine nucleotide exchange factor 11 (ARHGEF11) on 1q21 is a type 2 diabetes susceptibility gene in the Old Order Amish. Diabetes 2007;56:1363–1368. 17369523.
12. Etienne-Manneville S, Hall A. Rho GTPases in cell biology. Nature 2002;420:629–635. 12478284.
13. Hirosumi J, Tuncman G, Chang L, Gorgun CZ, Uysal KT, Maeda K, Karin M, Hotamisligil GS. A central role for JNK in obesity and insulin resistance. Nature 2002;420:333–336. 12447443.
14. Geiger PC, Wright DC, Han DH, Holloszy JO. Activation of p38 MAP kinase enhances sensitivity of muscle glucose transport to insulin. Am J Physiol Endocrinol Metab 2005;288:E782–E788. 15585585.
15. Jaffe AB, Hall A. Rho GTPases: biochemistry and biology. Annu Rev Cell Dev Biol 2005;21:247–269. 16212495.
16. Nevins AK, Thurmond DC. A direct interaction between Cdc42 and vesicle-associated membrane protein 2 regulates SNARE-dependent insulin exocytosis. J Biol Chem 2005;280:1944–1952. 15537656.
17. Kowluru A, Veluthakal R. Rho guanosine diphosphate-dissociation inhibitor plays a negative modulatory role in glucose-stimulated insulin secretion. Diabetes 2005;54:3523–3529. 16306371.
18. Khayat ZA, Tong P, Yaworsky K, Bloch RJ, Klip A. Insulin-induced actin filament remodeling colocalizes actin with phosphatidylinositol 3-kinase and GLUT4 in L6 myotubes. J Cell Sci 2000;113(Pt 2):279–290. 10633079.
19. Chiang SH, Baumann CA, Kanzaki M, Thurmond DC, Watson RT, Neudauer CL, Macara IG, Pessin JE, Saltiel AR. Insulin-stimulated GLUT4 translocation requires the CAP-dependent activation of TC10. Nature 2001;410:944–948. 11309621.
20. Rumenapp U, Blomquist A, Schworer G, Schablowski H, Psoma A, Jakobs KH. Rho-specific binding and guanine nucleotide exchange catalysis by KIAA0380, a dbl family member. FEBS Lett 1999;459:313–318. 10526156.
21. Chikumi H, Barac A, Behbahani B, Gao Y, Teramoto H, Zheng Y, Gutkind JS. Homo- and hetero-oligomerization of PDZ-Rho-GEF, LARG and p115RhoGEF by their C-terminal region regulates their in vivo Rho GEF activity and transforming potential. Oncogene 2004;23:233–240. 14712228.
22. Hsueh WC, St Jean PL, Mitchell BD, Pollin TI, Knowler WC, Ehm MG, Bell CJ, Sakul H, Wagner MJ, Burns DK, Shuldiner AR. Genome-wide and fine-mapping linkage studies of type 2 diabetes and glucose traits in the Old Order Amish: evidence for a new diabetes locus on chromosome 14q11 and confirmation of a locus on chromosome 1q21-q24. Diabetes 2003;52:550–557. 12540634.
23. Bottcher Y, Schleinitz D, Tonjes A, Bluher M, Stumvoll M, Kovacs P. R1467H variant in the rho guanine nucleotide exchange factor 11 (ARHGEF11) is associated with impaired glucose tolerance and type 2 diabetes in German Caucasians. J Hum Genet 2008;53:365–367. 18231709.
24. Report of the Expert Committee on the Diagnosis and Classification of Diabetes Mellitus. Diabetes Care 1997;20:1183–1197. 9203460.
25. Faber JJ, Green TJ. Foetal placental blood flow in the lamb. J Physiol 1972;223:375–393. 5039279.
26. Levy JC, Matthews DR, Hermans MP. Correct homeostasis model assessment (HOMA) evaluation uses the computer program. Diabetes Care 1998;21:2191–2192. 9839117.
27. International HapMap Consortium. The International HapMap Project. Nature 2003;426:789–796. 14685227.
28. Kang ES, Yun YS, Park SW, Kim HJ, Ahn CW, Song YD, Cha BS, Lim SK, Kim KR, Lee HC. Limitation of the validity of the homeostasis model assessment as an index of insulin resistance in Korea. Metabolism 2005;54:206–211. 15690315.

Article information Continued

Table 1

Clinical characteristics of the study population

Table 1

Values are presented as means ± standard deviation or number (%) unless otherwise indicated.

T2DM, type 2 diabetes mellitus; BMI, body mass index; LDL-C, low density lipoprotein cholesterol; NA, not applicable.

Table 2

Allelic and genotypic distributions of the R1467H variant in T2DM and control subjects

Table 2

Values are presented as number (%).

Recessive model: wild homozygote genotype vs. (risk homozygote genotype + heterozygote genotype).

T2DM, type 2 diabetes mellitus.

Table 3

Odds ratio estimates with confidence intervals of the R1467H variant in multiple ethnic populations

Table 3

T2DM, type 2 diabetes mellitus; IGT, impaired glucose tolerance; OR, odds ratio; CI, confidence interval.

aAnalysis of the association of the R1467H variant with patients was calculated using an additive model under a logistic regression model after adjusting for age, gender, and body mass index.

Table 4

Clinical and biochemical characteristics between R/R+R/H and H/H genotypes of ARHGEF11

Table 4

Data are means ± standard deviation unless otherwise indicated. Associations between categorical variables were analyzed using χ2 tests and continuous variables were assessed using Student's t-test.

T2DM, type 2 diabetes mellitus; BMI, body mass index.

aFasting insulin, HOMA2-IR and HOMA2-%β variables were log-transformed to correct for skew. In the R/R+R/H and H/H groups, fasting insulin, HOMA2-IR and HOMA2-%β measurements were available for only 534 and 21 patients, respectively.