Diabetes 56:1454-1459, 2007 DOI: 10.2337/db06-0640 © 2007 by the American Diabetes Association
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
1 Diabetes Molecular Genetics Section, Phoenix Epidemiology and Clinical Research Branch, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Department of Health and Human Services, Phoenix, Arizona Address correspondence and reprint requests to Leslie J. Baier, PhD, Diabetes Molecular Genetics Section, Phoenix Epidemiology and Clinical Research Branch, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, 445 N. Fifth St., Suite 210, Phoenix, AZ 85004. E-mail: lbaier{at}phx.niddk.nih.gov
Abbreviations:
IRS, insulin receptor substrate; LD, linkage disequilibrium; SNP, single nucleotide polymorphism
A prior genome-wide linkage scan in Pima Indians indicated a young-onset (aged <45 years) type 2 diabetes susceptibility locus on chromosome 1q21-q23. ARHGEF11, which encodes the Rho guanine nucleotide exchange factor 11, was analyzed as a positional candidate gene for this linkage because this protein may stimulate Rho-dependent signals, such as the insulin signaling cascade. The ARHGEF11 gene, and two adjacent genes NTRK1 and INSRR, were sequenced in 24 Pima Indians who were not first-degree relatives. Sequencing of the coding regions, 5' and 3' untranslated regions and putative promoter regions of these genes, identified 28 variants in ARHGEF11, 11 variants in NTRK1, and 8 variants in INSSR. These 47 variants, as well as 84 additional public database variants within/between these genes, were genotyped for association analysis in the same group of Pima Indians who had participated in the linkage study (n = 1,228). An R1467H in ARHGEF11, and several additional noncoding variants that were in high linkage disequilibrium with this variant, were nominally associated with young-onset type 2 diabetes (P = 0.01; odds ratio 3.39) after adjusting for sex, family membership, and Pima heritage. The risk allele H had a frequency of 0.10. In a subgroup of 262 nondiabetic, full-heritage Pima Indians who had undergone detailed metabolic testing, the risk allele H also was associated with a lower mean insulin-mediated glucose disposal rate and a lower mean nonoxidative glucose storage rate after adjusting for age, sex, nuclear family membership, and percentage of body fat (P 0.01). These findings suggest that variation within ARHGEF11 nominally increases risk of type 2 diabetes, possibly as a result of increased insulin resistance. The Pima Indians of Arizona have an extremely high prevalence of type 2 diabetes (1). Their diabetes is characterized by obesity, dysfunction of insulin secretion, insulin resistance (decreased insulin-mediated glucose disposal), and increased rates of postabsorptive endogenous glucose output (2). Studies (3–5) have shown that type 2 diabetes, insulin resistance, the acute insulin response, and obesity are highly heritable in this population. A prior genome-wide linkage scan in Pima Indians indicated a susceptibility locus for young-onset (onset age <45 years) type 2 diabetes on chromosome 1q21-q23 at D1S1677 (6). Linkage to type 2 diabetes on chromosome 1q21–23 has subsequently been observed in seven other populations who have formed the International Type 2 Diabetes 1q Consortium (7–14). Within this region of linkage, there is a cluster of three genes (ARHGEF11, NTRK1, and INSRR) located between 153 and 154 Mb, which encode proteins that have putative roles in the insulin signaling system. ARHGEF11 encodes the Rho guanine nucleotide exchange factor ARHGEF11 (also called PDZ-RhoGEF and KIAA0380). ARHGEF11 interacts with small GTPases (G-protein, guanine nucleotide-binding proteins), such as Rho, that function as molecular switches in signaling pathways that include the insulin signaling cascade (15,16). NTRK1 encodes the neurotrophic tyrosine kinase receptor type 1, which recruits insulin receptor substrate (IRS)-1 and IRS-2 (17), and INSRR encodes the insulin receptor–related receptor, which potentially phosphorylates IRS-1 and IRS-2 (available at http://genecards.bcgsc.bc.ca/cgi-bin/carddisp?insrr). In this study, ARHGEF11 and the two adjacent genes, NTRK1 and INSRR, were analyzed as positional candidate genes for type 2 diabetes in the Pima Indians.
The subjects who were sequenced (n = 24) and genotyped (n = 1,228) were participants in our prior genome-wide linkage study for diabetes susceptibility loci in Pima Indians (6) and are part of our ongoing longitudinal study of the etiology of type 2 diabetes among the Gila River Indian Community in central Arizona (1). All individuals in the longitudinal study are invited to participate in a standardized health examination biannually. To determine diabetes status, a 75-g orally administered glucose tolerance test is given and the results are interpreted according to the criteria of the World Health Organization (18). All studies were approved by the tribal council of the Gila River Indian Council and the institutional review board of the National Institutes of Diabetes and Digestive and Kidney Diseases. Among 1,228 subjects that were genotyped, 262 of the nondiabetic subjects had additionally been studied for measures of pre-diabetic phenotypes in our clinical research center. Only individuals (aged 18–45 years) who are confirmed to be healthy by medical history, physical examination, and routine laboratory tests and are not taking medications are studied. Oral glucose tolerance is measured after 2–3 days on a weight-maintaining diet of mixed composition. Blood for plasma glucose and insulin measurements is drawn before ingesting 75 g of glucose and at 30, 60, 120, and 180 min thereafter. Subjects also receive a 25-g intravenous injection of glucose over 3 min to measure the acute insulin response. Blood samples were collected before infusion and at 3, 4, 5, 6, 8, and 10 min after infusion for determination of plasma glucose and insulin concentrations. The acute insulin response was calculated as the mean increment in plasma insulin concentrations from 3 to 5 min (19). The hyperinsulinemic-euglycemic clamp technique was used to determine basal glucose appearance and insulin-stimulated glucose disappearance (uptake) rates (19). Briefly, insulin was infused to achieve physiologic plasma insulin concentrations (137 ± 3 µU/ml) for 100 min. Plasma glucose concentrations were held constant at 100 mg/dl by a variable 20% glucose infusion. Tritiated glucose was infused for 2 h before the insulin infusion to calculate glucose disappearance rates during the insulin infusion. During the last 40 min of the insulin infusion, the rate of insulin-stimulated glucose disposal was calculated, adjusted for steady-state plasma glucose and insulin concentration, and normalized to estimated metabolic body size (fat-free mass plus 17.7 kg) as described (19,20). Ventilated-hood indirect calorimetry was used to estimate rates of glucose and lipid oxidation before and during the insulin infusions (20). Glucose and lipid oxidation rates were normalized to estimated metabolic body size (fat-free mass plus 17.7) as described (19,20). Body composition was estimated by underwater weighing until January 1996 and currently is measured by dual-energy X-ray absorptiometry (DPX-1; Lunar Radiation, Madison, WI). A conversion equation derived from comparative analyses is used to make estimates of body composition comparable between methods (21).
Single nucleotide polymorphism identification and genotyping.
Statistical analysis.
Plasma insulin concentrations and glucose disposal rate during the physiological plasma concentration were log transformed before analyses to approximate a normal distribution. For analysis under an additive model, homozygotes for the major allele (1/1) and heterozygotes (1/2) and homozygotes for the minor allele (2/2) were coded to a continuous numeric variable for genotype (as 0, 1, and 2). The recessive model was defined as contrasting genotypic groups 1/1 vs. 1/2 + 2/2, where allele 1 is defined as the major allele. Although the general estimating equation approach accounts for the correlation among family members, it does not provide a specific within-family test of association and, thus, still is potentially confounded by population stratification. Therefore, data also were analyzed using a modification of the method of Abecasis et al. (22) to test for within-family association. In brief, this method partitions the association into between- and within-family components, represented, respectively, by the sibship mean of the continuous numeric variable for genotype and each individual's deviation from this mean. The test of the significance of the within-family components is a test of cotransmission among siblings, which is robust to population stratification (although it is less powerful than the more general test). Analyses of all traits were adjusted for potentially confounding variables (e.g., age, sex, percentage of body fat); all adjustments were specified a priori based on previous studies of the determinants of these traits. To examine pairwise linkage disequilibrium (LD), haplotype frequencies were estimated with the estimating haplotype (EH) program (Xie and Ott) (available at http://linkage.rockefeller.edu/ott/), and these haplotype frequencies were used to calculate D' and
Sequencing and genotyping of ARHGEF11. Sequencing of the 41 exons, exon-intron boundary regions, and 2 kb of the 5' (putative promoter) region of ARHGEF11 in 24 non–first-degree–related Pima Indians identified 28 variants, 3 of which were nonsynonymous amino acid substitutions (R293Q, S1456G, and R1467H). These 28 variants, and 66 additional database SNPs positioned outside of the regions that were sequenced, were genotyped for association analysis in the same Pima Indian subjects who were used for our prior type 2 diabetes linkage study (online appendix Table 1 [available at http://dx.doi.org/10.2337/db06-0640]). Three coding variants in ARHGEF11 were nominally associated with young-onset type 2 diabetes. The H allele (frequency = 0.10) of R1467H was associated with young-onset type 2 diabetes (P = 0.01; odds ratio 3.4 [95% CI 1.29–8.93], recessive model) after adjusting for sex, family membership, and Pima heritage (Fig. 1; Table 1). The G allele (frequency = 0.26) of S1456G also was associated with young-onset type 2 diabetes (P = 0.005, 1.85 [1.20–2.85], additive model; P = 0.012, 1.93 [1.16–3.22], recessive model) after adjusting for sex, family membership, and Pima heritage (Fig. 1; Table 1). The Q allele of R293Q was extremely rare (frequency = 0.002), and although the association data for this variant are included in Table 1, these data may not be reliable due to the limited sample power. Thirteen noncoding SNPs also were nominally associated with young-onset type 2 diabetes (listed in Table 1 legend). All of these noncoding SNPs, with the exception of rs3753213 (shown in Table 1), were in high LD (D' = 0.99, r2 = 0.99) with either R1467H or S1456G (online appendix Fig. 3).
We further assessed whether the R1467H and S1456G were associated with metabolic traits predictive of type 2 diabetes in Pima Indians who had not yet developed diabetes. In a subgroup of 262 nondiabetic, full-heritage Pima Indians who have undergone detailed metabolic testing, including measurements of body composition, oral glucose tolerance, insulin secretory function, insulin-stimulated glucose uptake, and indirect calorimetry, the risk allele (H) for R1467H was nominally associated with a lower mean glucose disposal rate (P = 0.005) and a lower mean nonoxidative glucose storage rate (P = 0.01) during a euglycemic-hyperinsulinemic clamp under physiologic concentrations of plasma insulin infusion, after adjusting for age, sex, nuclear family membership, and percentage of body fat (Table 2; Fig. 2). The R1467H was not associated with insulin secretion as assessed by the acute insulin response to a 25-g intravenous bolus injection of glucose or the insulin level at 30 min during an oral glucose tolerance test (Table 2). No association was found between S1456G and metabolic predictors of type 2 diabetes in nondiabetic, full-heritage Pima subjects (data not shown).
Sequencing and genotyping of NTRK1 and INSRR. Direct sequencing of all exons, exon-intron boundary regions of NTRK1 and INSRR, and 2.5 kb of the common 5' flanking regulatory region between these two genes in 24 non–first-degree–related Pima Indians identified 19 variants (11 in NTRK1 and 8 in INSRR; online appendix Table 1). Four of these variants predicted nonsynonymous amino acid changes (Phe316Leu, Arg869Leu, Pro928Leu, and Arg999STOP) located in INSRR. None of these coding variants were detected when INSRR previously was screened in Pima Indians using the technique of denaturing high-performance liquid chromatography on pooled DNA samples (23). The 19 variants, and 18 additional database SNPs positioned outside of the regions that were sequenced, were genotyped for association analysis in the same Pima Indian subjects who were used for the type 2 diabetes linkage study. None of the amino acid changes were associated with young-onset type 2 diabetes or metabolic predictors of type 2 diabetes in Pima Indians. Only SNP rs3753213, located in the shared common 5' flanking regulatory region (2.5 kb) between INSRR and NTRK1, showed a nominal association with young-onset type 2 diabetes (P = 0.04; Table 1; Fig. 1), but this variant was not associated with any metabolic predictors of type 2 diabetes in nondiabetic, full-heritage Pima Indians (data not shown). The variants detected in INSRR and NTRK1 are in a distinct LD block compared with variants in ARHGEF11 in Pima Indians (online appendix Fig. 3).
Haplotype analyses for ARHGEF11 and INSRR.
Our data suggest that the H allele at R1467H in ARHGEF11 (or a variant that is carried on the H allele in Pima Indians) increases risk of type 2 diabetes by reducing insulin-stimulated glucose uptake. Among the nondiabetic, full-heritage Pima Indians, subjects with an H allele (homozygotes and heterozygotes) had a lower mean glucose disposal rate during a hyperinsulinemic-euglycemic clamp after adjusting for age, sex, nuclear family membership, and percentage of body fat. This lower glucose disposal rate mainly was attributed to reduced nonoxidative glucose storage because the other component of glucose disposal rate, glucose oxidation, was similar between subjects with an H allele and subjects who were homozygous for the R allele. These findings indicate that ARHGEF11 may alter insulin-mediated glucose storage pathways, such as glycogen synthesis in the muscle and liver. ARHGEF11 is ubiquitously expressed and is present in both liver and muscle, which are important organs for insulin-mediated glucose metabolism (24). ARHGEF11 is a member of a subfamily of RhoGEFs that contain regulator of G-protein signaling domains. LARG (Leukemia-associated Rho guanine nucleotide exchange factor; also called ARHGEF12) also is a member of this regulator of G-protein signaling subfamily. We previously studied LARG as a candidate gene for type 2 diabetes and identified a functional Tyr1306Cys variant that also was associated with reduced insulin-stimulated glucose uptake among nondiabetic Pima Indians (16).
Although nominally significant associations were observed in our analysis of ARHGEF11, it should be cautioned that multiple comparisons inherent to association studies can lead to the discovery of false-positives, and for any given variant the prior probability of a true association is probably low. For these two reasons, particularly stringent P value thresholds for declaring a significant association have been proposed (25,26). The P values given in this study are unadjusted, and, if corrected for the 131 variants studied (which could be argued to be an overcorrection due to high LD among many of the variants), the P values of
Replication of associations in other populations also can provide evidence that associations are not spurious. The region spanning ARHGEF11, NTRK1, and INSRR also has been sequenced and densely genotyped in the Amish population, where variants within ARHGEF11 also were found to be associated with type 2 diabetes and oral glucose levels in response to an oral glucose tolerance test. However, the specific variants in ARHGEF11 demonstrating the best associations differed between the Amish and the Pima populations (28). The chromosomal region spanning ARHGEF11 also has been genotyped at a In summary, we propose that ARHGEF11 may have a minor, or modifying, role in influencing susceptibility to type 2 diabetes via its effect upon insulin action in Pima Indians. However, due to the relatively low frequency of risk allele H versus the high prevalence of type 2 diabetes in Pima Indians and the fact that this gene shows little evidence for association in other populations with linkage to type 2 diabetes on chromosome 1q21–23, except for the Amish, we speculate that there are additional type 2 diabetes susceptibility genes, with larger and more universal effects, on chromosome 1q21-q23.
This study was supported by the intramural research program of the National Institute of Diabetes and Digestive and Kidney Diseaes, National Institutes of Health (NIH). Funding was also provided by NIH Grant R01 DK073490. L.M. was supported by a mentor grant from the American Diabetes Association.
Published ahead of print at http://diabetes.diabetesjournals.org on 29 January 2007. DOI: 10.2337/db06-0640. Additional information for this article can be found in an online appendix at http://dx.doi.org/10.2337/db06-0640. The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. Received for publication May 9, 2006 and accepted in revised form January 22, 2007
This article has been cited by other articles:
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||