Diabetes 56:2643-2648, 2007 DOI: 10.2337/db06-1023 © 2007 by the American Diabetes Association
A Genome-Wide Linkage Scan of Insulin Level–Derived TraitsThe Amish Family Diabetes Study
1 Department of Medicine, School of Medicine, University of California, San Francisco, California Address correspondence and reprint requests to Wen-Chi Hsueh, MPH, PhD, 2200 Post St., C433, San Francisco, CA 94143-1640. E-mail: wen-chi.hsueh{at}ucsf.edu
Abbreviations:
AUC, area under the curve; HOMA-IR, homeostasis model assessment of insulin resistance; IFG, impaired fasting glucose; IGT, impaired glucose tolerance; OGTT, oral glucose tolerance test; STR, short tandem repeat
OBJECTIVE—Serum insulin levels are altered in insulin resistance and insulin deficiency, states that are associated with the development of type 2 diabetes. The goal of our study was to identify chromosomal regions that are likely to harbor genetic determinants of these traits. RESEARCH DESIGN AND METHODS—We conducted a series of genetic analyses, including genome-wide and fine-mapping linkage studies, based on insulin levels measured during an oral glucose tolerance test (OGTT) in 552 nondiabetic participants in the Amish Family Diabetes Study. Indices of insulin secretion included the insulinogenic index and insulin at 30 min postglucose load (insulin 30), while indices of insulin resistance included homeostasis model assessment of insulin resistance (HOMA-IR) and fasting insulin. Insulin area under the curve, a measure of both insulin secretion and insulin resistance, was also examined.
RESULTS—All traits were modestly heritable, with heritability estimates ranging from 0.1 to 0.4 (all P < 0.05). There was significant genetic correlation between fasting insulin and HOMA-IR ( CONCLUSIONS—These results suggest that there may be different genes influencing variation in OGTT measures of insulin secretion and insulin resistance. Type 2 diabetes is a classic example of a complex disease that results from the interaction of multiple genetic and environmental factors. One strategy to identify the genes regulating the type 2 diabetes phenotype is to focus on related subclinical (intermediate) phenotypes, which are likely to be less genetically complex and involve fewer alleles. The pathophysiology of type 2 diabetes involves defects in insulin sensitivity and/or insulin secretion thus making these measures excellent intermediate traits to study in order to dissect the genetic underpinnings of type 2 diabetes. Circulating levels of plasma insulin, measured in either the fasting state or in response to a glucose load, vary considerably in nondiabetic individuals, with higher levels predicting future development of diabetes. Several indices of insulin secretion and insulin resistance can be derived from insulin and glucose levels measured at different time points during an oral glucose tolerance test (OGTT). OGTT-derived measures of insulin secretion include the insulinogenic index and insulin at 30 min following a 75-g oral glucose load (insulin 30), while measures of insulin resistance include fasting insulin and the homeostasis model assessment of insulin resistance (HOMA-IR). Insulin area under the curve (insulin AUC) during the OGTT reflects elements of both insulin secretion and insulin resistance. Several studies have investigated the genetic epidemiology of fasting insulin levels and OGTT-derived indices (1,2). Genome-wide scans have identified several chromosomal regions linked to these traits (3–13), but few have been replicated. The goal of this study was to characterize the genetic epidemiology of five insulin-related traits, including fasting insulin, insulin 30, insulinogenic index, HOMA-IR, and insulin AUC. We first estimated their heritability, then assessed whether and to what degree these traits may share common genetic influence, and performed genome-wide and fine-mapping linkage analyses of these traits. We found that these insulin traits were significantly heritable and that there may be different genetic influences underlying these OGTT-derived measures of insulin secretion and resistance. Furthermore, our genome scan provides evidence of linkage on chromosomes 2p and 15q to measures of insulin secretion. Clinical characteristics of 552 nondiabetic subjects are shown in Table 1. The overall prevalence of impaired glucose tolerance (IGT) and/or impaired fasting glucose (IFG) was 22% in this sample. Insulin resistance–related measures in those with IGT/IFG were significantly higher (P < 0.001) than in euglycemic individuals, whereas insulin secretion–related measures were similar between the two groups (data not shown).
Heritability (h2) estimates for insulin-related traits ranged from 0.12 to 0.36 (all P < 0.05) (Table 1). We next examined the extent to which these traits share common genetic influences in order to assist our interpretation of subsequent linkage analyses. As shown in Table 2, high genetic correlations were observed between fasting insulin and HOMA-IR, suggesting that these two measures share very substantial genetic components (estimated G = 0.86, P < 0.05). These findings are not unexpected since HOMA-IR is derived from the fasting insulin. Insulin 30 had a strong genetic correlation with insulinogenic index ( G = 0.81, P < 0.001). Insulin AUC was more genetically correlated with insulin 30 than with any other trait ( G = 0.82, P < 0.001). There was no significant genetic correlation between fasting insulin and insulin 30 or the insulinogenic index.
Our genome-wide linkage analysis identified two chromosomal regions with suggestive evidence for linkage (defined as P < 0.001 or logarithm of odds [LOD] >2.07) to insulin secretion traits. The first region was for insulin 30 (LOD 2.53, P = 0.00032), occurring at 73 cM on chromosome 15q23 (nearest marker D15S131) (Fig. 1A), which was also supported by the two-point analysis of D15S131 (LOD 2.82, P = 0.00016). The second region was for insulinogenic index at 121 cM on chromosome 2p13 (LOD 2.51, P = 0.00034; nearest marker D2S139) (Fig. 1B), as supported by two-point analysis of D2S139 (LOD 1.45, P = 0.0049). Similar results were obtained when the analysis was performed without adjustment for BMI. Nine other chromosomal regions (on chromosomes 1q21, 2q21-q22, 6p24-p23, 7q11-q21, 7q31, 10q11, 11q21-q23, 12p13-p12, and 19p13) showed linkage signals to one or more insulin traits with LOD 1.18 (P < 0.01) (Table 3). Bivariate linkage analysis for the two pairs of highly correlated insulin secretion traits was conducted. The maximum LOD score for (insulin 30 + the insulinogenic index) was 2.21 at 69 cM from the p-ter of chromosome 15. Insulin 30 + insulin AUC showed a lower LOD score (1.58) at 73 cM, also on chromosome 15. These results suggest that the genes in this region may contribute more to the insulin 30 phenotype than to the other two traits. The complete genome scan results can be viewed in the online appendix (available at http://dx.doi.org/10.2337/db06-1023).
To follow up on the linkage signal to insulin 30 on chromosome 15q, 14 additional short tandem repeat (STR) markers were genotyped to increase the information content in the region. These markers reduced the marker density from 11.1 to 3.2 cM and further refined the linkage peak to near D15S153 (at 68 cM, LOD 2.38), 9 cM closer to the centromere compared with the original genome scan results, which was well within the 1-LOD interval.
Insulin secretion and insulin sensitivity are important determinants of glucose homeostasis and diabetes. The current gold standard for quantifying insulin secretion is the acute insulin response to intravenous glucose (AIRg) test and, for insulin sensitivity, the hyperinsulinemic-euglycemic clamp. Since it is not practical to obtain these labor-intensive and costly phenotypes in large numbers of subjects for genetic studies, many have used surrogate measures of insulin secretion and insulin sensitivity. Previous studies have shown that the correlation between insulin-related indices from OGTT investigated in this study and those measured in the clamp were moderate (r = 0.3–0.7) but significant (14,15). In contrast, Bergman et al. (16) showed that in African and Hispanic Americans, the Our systematic genetic analyses of plasma insulin levels during a 3-h OGTT and indices derived from these measures provide insights that are likely to be useful in the search for genetic influences on insulin secretion and sensitivity. We demonstrate modest levels of heritability for a number of insulin secretion–related traits that have not been previously reported, such as insulin AUC, insulinogenic index, and insulin 30. The level of heritability for fasting insulin levels in the Amish was relatively low (h2 = 0.13) compared with values reported in Caucasians (0.37–0.47) (17,18), African Americans (0.28) (18), Mexican Americans (0.38–0.53) (4,19), Pima Indians (0.26) (15), and Asians (0.43) (20), although not much different from that estimated in African Americans and Hispanics from the Insulin Resistance and Atherosclerosis Study (h2 = 0.08–0.17) (16,21). As expected, based on the fact that HOMA-IR is derived from fasting insulin, both of these traits are phenotypically and genetically highly correlated. These two traits may reflect relatively more of the insulin-resistance phenotype. HOMA-IR is commonly used in epidemiological studies as a proxy measurement of insulin resistance, yet it did not appear to be more informative compared with fasting insulin in our genetic analysis. This is likely due to the fact that fasting insulin is a primary component of HOMA-IR calculation, particularly in nondiabetic subjects for whom there is less variation in glucose levels. Our results suggest that a simpler measurement (fasting insulin) may serve as well as a more complicated composite measurement (HOMA-IR) for genetic studies of insulin resistance. On the other hand, insulin 30 shared a significant genetic component with the insulinogenic index—both measures are thought to be more related to insulin secretion. Importantly, there was no significant genetic correlation between fasting insulin and insulin 30, suggesting that these two traits are likely to have different genetic influences. These findings are consistent with the biological concept that defects in insulin sensitivity and insulin secretion are genetically distinct. From our genome-wide linkage analysis, we identified two chromosomal regions, 15q23 for insulin 30 and 2p13 for insulinogenic index, with suggestive evidence of linkage. For the linkage signal for insulin 30 on chromosome 15q, the 1-LOD support interval is a 31.4-cM region defined by markers D15S117 and D15S158. Within this interval, there are 302 genes (263 named, 59 predicted [NCBI Build 35.1]). Several genes within our region of linkage on 15q23 may potentially be associated with ß-cell development and function. Some of the candidate genes in this region are involved with hormone secretion (SCAMP2 and SCAMP5) (22,23), others are growth factors (NRG4) (24), while others such as ISL2 are homologous to proteins known to be involved in ß-cell development (ISL1) (25). For the linkage signal for insulinogenic index on chromosome 2p, the 1-LOD support interval is a 26-cM region defined by markers D2S139 and D2S347. Within this interval, there are 194 genes (167 named, 27 predicted). A number of candidate genes for ß-cell development or function lie within the 1-LOD region, including two genes encoding secretory vesicle associated membrane proteins VAMP5 and VAMP8. It is somewhat surprising that there was little overlap of linkage signals for insulin 30 and insulinogenic index on chromosomes 2 and 15 given that genetic correlation was high. However, these traits were not perfectly correlated and thus may not have an identical genetic basis. Furthermore, gene-by-gene or gene-by-environment interactions, which were not accounted for in estimates of genetic correlation, may contribute to the lack of significant overlapping linkage signals in these regions.
Thus far, genome-wide linkage studies of diabetes-related traits, including insulin traits, have been reported in several ethnic groups (3–9). Studies conducted in Japanese and Mexican Americans have also revealed linkage to type 2 diabetes on chromosome 15q. In a Japanese study (8), a suggestive linkage signal to type 2 diabetes was observed on 15q13-q21 (maximum LOD score = 2.19), whose 1-LOD region overlaps with the 1-LOD region for the linkage to insulin 30 in the Amish. A small linkage signal to type 2 diabetes (two-point maximum LOD score for D15S119 = 1.50) was also observed in the same region in a study of Mexican Americans (9). On the other hand, evidence for replication of our linkage on chromosome 2p13 to insulinogenic index (LOD 2.51 or empirical P = 0.00034, at 121 cM) is more limited. A meta-linkage analysis utilizing information from four ethnic groups of the National Heart, Lung, and Blood Institute Family Blood Pressure Program (3) reported suggestive linkage signals to both fasting insulin and HOMA-IR (LOD 2.3–2.6 or empirical P = 0.03–0.06, at
We also observed LOD scores In summary, we observed evidence for linkage of insulin traits to regions of chromosome 2p13 and 15q23—both regions with a number of ß-cell candidate genes. Further examination of these regions through linkage disequilibrium mapping and positional candidate gene analysis will be necessary to identify the genes and their functional variants, which should be relevant not only to their influences on insulin levels but also to susceptibility to type 2 diabetes.
Our study was based on 691 members of the Amish Family Diabetes Study. Details of subject recruitment have been previously reported (26). The study protocol was approved by the institutional review board at the University of Maryland, Baltimore, and informed consent was obtained from each participant.
Phenotypes.
Genotypes.
Statistical analysis.
All analyses were conducted using a pedigree-based variance components method. We first estimated heritability then genetic correlation. Bivariate modeling was used to partition the phenotypic correlation (
This study was supported in part by a research grant from GlaxoWellcome, Inc. and Axys Pharmaceuticals; National Institutes of Health grants R01 DK54361, K24 DK02673, U01 DK58026, R01 AG023692, R01 DK068495, and K01 AG022782; and the American Diabetes Association. Funding and support was also provided by the University of Maryland General Clinical Research Center (Grant M01 RR 16500), General Clinical Research Centers Program, National Center for Research Resources, and the Baltimore Veterans Administration Geriatric Research and Education Clinical Center. We thank the Amish Research Clinic Staff for their energetic efforts in study subject recruitment and characterization, Drs. Alejandro Schäffer and Richa Agarwala for assistance in pedigree construction, and Dr. Pamela St. Jean for helpful comments on the manuscript. This study would not have been possible without the outstanding cooperation of the Amish community.
Published ahead of print at http://diabetes.diabetesjournals.org on 23 July 2007. DOI: 10.2337/db06-1023. C.J.B. is currently affiliated with EmerGen, Inc., Salt Lake City, Utah. Additional information for this article can be found in an online appendix at http://dx.doi.org/10.2337/db06-1023. 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 July 24, 2006 and accepted in revised form July 17, 2007
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||