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.
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.
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.
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).
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 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.
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).
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.
Figures and Tables
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.
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.
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.
3. Rich SS. Mapping genes in diabetes. Genetic epidemiological perspective. Diabetes. 1990. 39:1315–1319.
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.
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.
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.
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.
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.
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.
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.
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.
12. Etienne-Manneville S, Hall A. Rho GTPases in cell biology. Nature. 2002. 420:629–635.
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.
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.
15. Jaffe AB, Hall A. Rho GTPases: biochemistry and biology. Annu Rev Cell Dev Biol. 2005. 21:247–269.
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.
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.
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.
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.
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.
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.
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.
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.
24. Report of the Expert Committee on the Diagnosis and Classification of Diabetes Mellitus. Diabetes Care. 1997. 20:1183–1197.
25. Faber JJ, Green TJ. Foetal placental blood flow in the lamb. J Physiol. 1972. 223:375–393.
26. Levy JC, Matthews DR, Hermans MP. Correct homeostasis model assessment (HOMA) evaluation uses the computer program. Diabetes Care. 1998. 21:2191–2192.
27. International HapMap Consortium. The International HapMap Project. Nature. 2003. 426:789–796.
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.