| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Diabetes and Genetic Epidemiology Unit (S.K., V.H., J.E., J.T.), Department of Epidemiology and Health Promotion, and Department of Mental Health and Alcohol Research (J.K.), National Public Health Institute, FIN-00300 Helsinki, Finland; Division of Diabetes and Endocrinology (S.K., N.T.), Department of Internal Medicine, Jikei University School of Medicine, 105-8461 Tokyo, Japan; Department of Medicine (M.L.), Helsinki University Hospital, FIN-00029 HUS Helsinki, Finland; Department of Public Health (J.K., J.T., S.K., J.E.), University of Helsinki, Helsinki FI-00014, Finland; Department of Public Health (M.K.), University of Turku, FI-20520 Turku, Finland; and South Ostrobotnia Central Hospital (J.T.), FIN-60101 Seinäjoki, Finland
Address all correspondence and requests for reprints to: Jaakko Kaprio, M.D., Ph.D., Department of Public Health, P.O. Box 41 (Mannerheimintie 172), University of Helsinki, FIN-00014 Helsinki, Finland. E-mail: jaakko.kaprio{at}helsinki.fi.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
There is no clear consensus on the degree to which genes account for the variance of insulin sensitivity and secretion. Both methodologies to measure traits and the selection of study subjects influence the results. A study with 209 middle-aged Dutch twin pairs estimated that heritability (h2), the proportion of the total variance accounted for by genes, was 2025% for fasting insulin, which is a proxy for insulin resistance (13). Danish researchers studied the insulin and glucose responses for oral glucose tolerance test (OGTT) in 303 twin pairs and concluded that only 19% of the variance in fasting insulin could be attributed to genetic factors (20). The former of these two studies included a bivariate variance component analysis of fasting glucose and insulin values; no postchallenge data, however, were available. The latter of the studies involved both fasting and postchallenge concentrations of glucose and insulin, but no bivariate analysis based on these values was performed.
We, on the other hand, have shown, by applying metabolic techniques more precise than OGTT to a limited number of twin pairs, that polygenic factors accounted for more than 55% of the variance in insulin secretion (iv glucose tolerance test) and 37% of the variance in glucose uptake (euglycemic hyperinsulinemic clamp) (19).
The objective for the current study was to test whether our observations (19) on the heritability of insulin sensitivity and insulin secretion could be replicated when the number of twin pairs was doubled and a more clinical study approach (OGTT) applied. The size of the study population also allowed for the evaluation of sex-associated differences of h2 estimates as well as the study of possible genetic covariance between any of the traits measured both in the fasting and postchallenge state of OGTT.
| Subjects and Methods |
|---|
|
|
|---|
|
Height and weight were recorded with light indoor clothes and body mass index (BMI) was calculated as weight/height2 (kilograms per square meter). Waist circumference (to the nearest centimeter without clothes) was measured with a nonelastic soft tape midway between the lowest rib and the iliac crest while subjects were standing. Hip circumference (to the nearest cm without clothes) was measured over the widest part of the gluteal region. The waist to hip ratio (WHR) was also calculated.
A standard 75-g OGTT was carried out. Blood samples for the analysis of serum insulin (Ins) and plasma glucose (Gluc) were drawn in the fasting state (Ins_0, Gluc_0), and 2 h after the oral glucose load (Ins_120, Gluc_120), respectively.
Analytical methods
Serum insulin was measured using RIA (Pharmacia, Uppsala, Sweden). Plasma glucose was measured by the glucose oxidase method (Glucose Analyzer II, Beckman Instruments, Fullerton, CA, and Cobas Mira Plus analyzer, Roche Instrument Center, Rotkreuz, Switzerland).
Statistical methods
The homeostasis model insulin resistance index (HOMA-R) [fasting glucose (millimoles per liter) x fasting insulin (milliunits per liter)/22.5] and ß-cell function (HOMA-BETA) [fasting insulin (milliunits per liter) x 20/(fasting glucose (millimoles per liter) 3.5)] were calculated as suggested by Matthews et al. (28). For the calculation of phenotypic correlations between variables, all except waist (square root) and HOMA-BETA values (reciprocal of the value) were log transformed. The descriptive statistics and phenotypic correlations were calculated using STATA 7-software (Stata Corp. LP, College Station, TX), which also enables the computation of correct SE and P values in a clustered population (i.e. by taking into account that the study sample consisted of twin pairs).
For the estimation of variance components, log-transformed values were used in all variables except glucose, which fitted best as such. Intraclass correlations, i.e. within twin pair-correlation coefficients (intraclass r), were calculated to first test for the existence of a genetic variance component for each trait; higher MZ than DZ correlations indicate that genetic effects may exist. Intraclass correlations were computed using the Twinan software (29).
After that, standard univariate analyses were carried out using Mx (30), a custom-built program for the modeling of genetically informative data (http://www.vcu.edu/mx/). Finally, bivariate analyses were performed to study the degree of correlation of the latent genetic and environmental factors influencing both Ins_0 and Ins_120. Variance components were estimated by applying the maximum likelihood method to fit the models to the raw data, using age and sex as covariates. Models were fitted to the data, assuming that variation in the phenotype was attributable to a combination of a nonshared environmental component (E), a shared environmental component (C), an additive genetic component (A), and genetic effects due to dominance (D). The estimates of the variance components were calculated as the proportion of variance divided by the total variance. The Akaikes information criterion value was used to assess the goodness of fit of the model (31). The model with the lowest Akaikes information criterion value reflects both the goodness of fit of a model and its parsimony, i.e. the ability to account for the observed data with few parameters. The superiority of alternative, hierarchically nested models was tested by hierarchic
2 tests because the difference between the value for the reduced model and that of the full model (
2) is
2 distributed, with degrees of freedom (df) equal to the differences in df of the models to be compared. In the case in which covariance structure was different between men and women, analyses were done separately for men and women (32). The proportion of A in relation to total variance (percent) was called heritability, and expressed as h2. Heritability is a population-specific characteristic that has no interpretation on the level of the individual or family. To interpret the results of model fitting, a comparison of the trait intraclass correlation (intraclass r, between twin pair members) coefficient values between MZ and DZ twins was applied (33).
| Results |
|---|
|
|
|---|
Table 1
shows the number of twin pairs and the summary of descriptive statistics. The proportion of female subjects was 58% among MZ and 61% among DZ pairs, respectively. There were no significant differences in means or variances between individuals from MZ and DZ pairs in study variables, except for Ins_120, which was lower in female (P = 0.04) MZ pairs, compared with DZ pairs. Also, both systolic (P = 0.01) and diastolic (P = 0.01) blood pressure mean values were higher in individuals from MZ pairs, compared with DZ pairs.
The phenotypic correlations among BMI, waist, Gluc_0, Gluc_120, Ins_0, Ins_120, HOMA-R, and HOMA-BETA are shown in Table 2
. All correlations were statistically significant at the level of P < 0.01, with the exception of that between HOMA-BETA and Gluc_120 (P = 0.91).
|
The correlation coefficient between waist and WHR was 0.82, and the paired correlation coefficients of WHR with other traits were almost identical with the corresponding correlations between waist and other traits (data not shown). Therefore, only waist was used in the intraclass r and variance component analyses described below.
Intrapair correlations
All intrapair correlation coefficients (r) were significant, and higher, in MZ than DZ twins, with the exception of Ins_120 in male twins (MZ intraclass r = 0.24, N.S.; DZ 0.48, P = 0.002) and male and female twins together (MZ intraclass r = 0.32, P = 0.006; DZ 0.40, P < 0.001), respectively. For BMI, female waist, Ins_0 (male and female twins together), and HOMA-R, the MZ intraclass r was more than twice the DZ value. The MZ intraclass r was less than twice the DZ intraclass r in male waist and Ins_120 in both sexes as well as waist, Gluc_0, Gluc_120, and HOMA-BETA in men and women together (Table 3
).
|
Univariate analysis.
Table 4
shows the best-fitting models for each trait, based on the univariate analysis of the variance components. The model with component E (unique and error variance) alone could be rejected in all traits (data not shown). The AE model was the best fitting for BMI, waist female, waist overall, Ins_0, female Ins_120, Gluc_0, Gluc_120, HOMA-R, and HOMA-BETA. The CE model fitted best for male waist and Ins_120 in men and overall. The h2 estimates were 68% for both BMI and for waist when male and female twins were analyzed together. In female twins, the h2 for waist was 70%. Whereas the corresponding estimate could not be estimated for male twins, in which the best-fitting model was CE. The h2 for Ins_0 was 43% in male and 42% in female twins, respectively. When the most parsimonious model was used for Ins_120, only female h2 (51%) could be estimated. The h2 were 45% for Gluc_0, 35% for Gluc_120, 42% for HOMA-R, and 38% for HOMA-BETA, respectively.
|
2 (CE)
2 (ACE) = 828.034 827.200 = 0.834, difference in df = 1; P = 0.361) or common environmental influence alone (for
2 (AE)
2 (ACE) = 828.932 827.200 = 1.731, difference in df = 1; P = 0.188) could be discarded but not both of them together (for
2 (E)
2 (ACE) = 844.163 827.200 = 16.963 with difference in df = 2; P < 0.0001). However, testing E alone against the AE model (for
2 (E)
2 (AE) = 844.163 829.304 = 14.859 with difference in df = 1; P = 0.00012) disclosed that Ins_0 and Ins_120 were influenced by a shared genetic factor. Using this AE model as the final model, the genetic correlation coefficient between the latent polygenetic factors, influencing the variances of both fasting and 2-h (postchallenge) insulin values, was estimated to be 0.812 [95% confidence interval (CI) 0.5910.990] and the correlation between the unique environmental influences 0.540 (95% CI 0.3640.680), respectively. | Discussion |
|---|
|
|
|---|
We observed a difference in heritability estimates between men and women for visceral fat accumulation, measured as waist circumference, and postchallenge (120-min) insulin levels, both of which are important risk factors for coronary heart disease (34, 35). Although the overall phenotypic correlation coefficient (Pearsons r) between waist and 120-min insulin was only 0.25 (P < 0.0001) (Table 2
), it was possible to detect a difference in the respective sex-specific values: 0.09 (P = 0.32) for male and 0.45 (P < 0.0001) for female twins, respectively (data not shown in the results). This observation is in line with that of Van Gaal et al. (36), who reported that the level of insulin after 120 min could be explained up to 4080% by sex hormones and WHR in women. Fasting serum insulin is mainly cleared by the liver, whereas postchallenge insulin is cleared by several of the hormones target organs (i.e. liver, muscle, and adipocytes) (37, 38, 39, 40). Is it possible that visceral fat accumulation affects liver function and postchallenge insulin clearance (36, 41, 42, 43) in women more than in men? A sex-specific genetic difference of abdominal obesity followed by decreased plasma adiponectin levels and increased oxidative stress might correlate with peripheral insulin resistance and postprandial insulin levels (44, 45). But the sex-specific genes have not been clearly identified; therefore, to determine whether there is a sex-specific genetic difference, we would need information on opposite-sex (one male, one female) twin pairs or other sibling pairs of different sexes. So-called sex-limitation models could then be used to test whether such sex-specific effects exist in a future study. For such an assessment, a much larger number of twin pairs would, however, be needed, with opposite-sex twin pairs.
The current h2 for insulin resistance, measured here both by the use of fasting insulin and of HOMA-R, was in line with results obtained in several previous studies (8, 46, 47, 48) and somewhat higher than that estimated by us when a clamp technique was applied (19). This difference may in part be explained by the fact that fasting insulin and HOMA-R represent the fasting state, whereas clamp glucose uptake reflects insulin sensitivity during more or less constant hyperinsulinemia.
Decomposition of variance components of insulin and glucose concentrations during the OGTT has been performed by several research groups (15, 21, 36, 47). Genetic covariance between trait values measured at two separate moments during the OGTT has, however, not been evaluated. The observation of a strong correlation between genetic variance components of fasting and 120-min insulin is in line with a recent study of Finnish families in which the common genetic architecture of type 2 diabetes was identified to include both fasting and 120-min insulin (47). The observations are, furthermore, compatible with findings by Schumacher et al. (49), suggesting a major gene locus contributing to both fasting and 60-min insulin concentrations.
Conclusions
The strong observed input from environmental factors to peripheral insulin resistance justifies strategies that address the importance of lifestyle changes in the prevention of type 2 diabetes.
Whether abdominal obesity and postprandial insulin levels are unequally associated in men and women, and whether this difference is due to a sex-specific genetic covariance between the traits, could be resolved in studies specifically designed for such a purpose.
| Footnotes |
|---|
First Published Online February 22, 2005
1 S.K. and M.L. share the first authorship ![]()
Abbreviations: A, Additive genetic component; BMI, body mass index; C, environmental component; CI, confidence interval; df, degrees of freedom; DZ, dizygotic; Gluc, plasma glucose; E, environmental component; h2, heritability; HOMA-BETA, homeostasis model of ß-cell function; HOMA-R, homeostasis model of insulin resistance index; Ins, serum insulin; MZ, monozygotic; OGTT, oral glucose tolerance test; WHR, waist to hip ratio.
Received December 16, 2004.
Accepted February 11, 2005.
| References |
|---|
|
|
|---|
gene on glucose tolerance and plasma insulin profiles in monozygotic and dizygotic twins: thrifty genotype, thrifty phenotype, or both? Diabetes 52:194198This article has been cited by other articles:
![]() |
K. Beardsall, K. K. Ong, N. Murphy, M. L. Ahmed, J. H. Zhao, M. W. Peeters, and D. B. Dunger Heritability of Childhood Weight Gain from Birth and Risk Markers for Adult Metabolic Disease in Prepubertal Twins J. Clin. Endocrinol. Metab., October 1, 2009; 94(10): 3708 - 3713. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Uribarri, W. Cai, M. Peppa, S. Goodman, L. Ferrucci, G. Striker, and H. Vlassara Circulating Glycotoxins and Dietary Advanced Glycation Endproducts: Two Links to Inflammatory Response, Oxidative Stress, and Aging J. Gerontol. A Biol. Sci. Med. Sci., April 1, 2007; 62(4): 427 - 433. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. F. Lawler, K. Chase, R. Teckenbrock, and K. G. Lark Heritable Components of Feline Hematology, Clinical Chemistry, and Acid-Base Profiles J. Hered., December 8, 2006; (2006) esl041v1. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Endocrinology | Endocrine Reviews | J. Clin. End. & Metab. |
| Molecular Endocrinology | Recent Prog. Horm. Res. | All Endocrine Journals |