The relevance of the insulin gene (INS) variable number tandem repeat (VNTR) polymorphism to indexes of body size and adult obesity is inconclusive. Given the equivocal reports on the association between the VNTR class genotype at the insulin gene locus and indexes of body size and obesity, we assessed these associations in a series of cohort studies based on 7,999 middle-aged men and women. We found no convincing evidence that INS VNTR class genotype was associated with indexes of body size and adult obesity. These data suggest that INS VNTR class is not an important determinant of size and body weight regulation in middle-aged men and women.
Experimental and observational studies have suggested that insulin signaling is involved in regulation of somatic growth and body weight (1). The variable number tandem repeat (VNTR) polymorphism, lying 396 base pairs 5′ to the insulin gene on chromosome 11p15.5, is thought to influence insulin expression (2,3). This minisatellite polymorphism has also been associated with altered expression of the nearby IGF-II, a key regulator of fetal growth and development (4,5). Collectively, observational studies have reported conflicting associations between the INS VNTR class genotype and indexes of body size in childhood and obesity (4,6–8). Importantly, the relevance of INS VNTR class to indexes of body size and obesity in adults is also inconclusive. Given these equivocal results, we assessed the association between INS VNTR class genotype and indexes of body size and adult obesity in a series of cohort studies based on 7,999 middle-aged men and women.
RESEARCH DESIGN AND METHODS
We used data from the Norfolk cohort of the European Prospective Investigation into Cancer (EPIC) and the Ely Study. EPIC is an ongoing prospective cohort study of ∼25,000 people aged 45–75, resident in Norfolk, and recruited from general practice registers between 1993 and 1997 (9). It is an ethnically homogeneous Caucasian population. These participants completed a health examination. In January 1998, we invited the cohort for a second health examination, and 15,786 people had attended by October 2000. All participants gave informed consent. Details of recruitment, anthropometric measurements, and health examinations following standardized protocols have been published (9). All participants completed a detailed health and lifestyle questionnaire, which included a question on birth weight. A d-loop nonstretch fiberglass tape was used for the head circumference measures. DNA was extracted from the 15,786 whole blood/EDTA samples (Whatman Biosciences, Ely, U.K.) taken at the second health examination. BMI was calculated from the weight and height measures.
For the purposes of this analysis, we examined data from two subcohorts of the 15,786 EPIC study participants who attended the second health examination (9). EPIC1 comprises 3,744 participants with fully arrayed DNA available for genotyping. EPIC2 is a random sample of 5,000 participants who were free of baseline disease (cancer, coronary heart disease, and diabetes) and had fully arrayed DNA available, completed food frequency questionnaire data, and had blood HbA1c (A1C) and BMI measured at the two clinical assessments. We excluded 759 participants from EPIC2 as they were also present in the EPIC1 sample set. However, we used this duplicate set to assess concordance between genotyping methods, as described below.
We also used data from the Medical Research Council Ely Study, a prospective population-based cohort study of the etiology and pathogenesis of type 2 diabetes and related metabolic disorders (10). Again, the study population consists of an ethnically homogeneous Caucasian population in which phenotypic data have been recorded at clinical assessments using standardized protocols for anthropometric measurements. All participants were between 40 and 65 years of age at baseline and completed a health and lifestyle questionnaire, which included a question on birth weight. This cohort was recruited from a population-sampling frame with a high response rate (74%), making it representative of the general population for this area of eastern England.
Study design and genotyping.
The INS VNTR allele length falls into two broad classes in European populations: class I (26–63 repeats) and class III (141–209 repeats). We genotyped two subcohorts of the EPIC-Norfolk study (EPIC1 [n = 3,446] and EPIC2 [n = 3,720]) and 833 participants in the Ely study for the −23Hph1 polymorphism. There is 99.6% concordance between the A to T polymorphism at the −23Hph1 site CCACT at nucleotide 2,401 of the insulin gene (accession #V00565 and rs689) and the adjacent INS VNTR class I and III alleles in European populations (11). Hence −23/Hph1 genotype infers INS VNTR class I/III genotype. Genotype data for the EPIC1 study set were obtained from Invader assays (Third Wave Technologies, Madison, WI) and for the EPIC2 and the Ely study sets with TaqMan chemistry (Applied Biosystems, Warrington, U.K.), as previously reported (12). We assessed concordance between methods by typing a subset of 759 samples from EPIC1 and EPIC2 using both methods.
Statistical analysis.
We tested for Hardy–Weinberg equilibrium using the χ2 test. We calculated the effect of genotype on continuous anthropometric and metabolic characteristics using linear regression. In the primary analysis, we used an additive genetic model, assuming a linear relationship between the number of T (class III) alleles and the trait of interest. In a secondary analysis, we used a recessive model, comparing class III homozygotes with class I carriers, as demonstrated previously (8). Possible interactions between covariates and genotype were assessed with log-likelihood ratio tests. Stata version 8.2 was used for all analyses (Stata Corp, College Station, TX).
RESULTS
In the three cohort studies, the −23Hph1 gene variant was in Hardy–Weinberg equilibrium (EPIC1 P = 0.71, EPIC2 P = 0.24, and Ely P = 0.36). In a subgroup of 759 duplicate samples from EPIC1 and EPIC2, concordance between the two genotyping methods was 99%. There was no material difference in the association between genotype and indexes of body size and obesity in sex-stratified analyses; therefore, all data are presented for men and women combined. Table 1 shows demographic and anthropometric characteristics according to −23Hph1 genotype for the three cohorts. We found that the proportion of women varied by genotype (P = 0.02) in the EPIC1 study set; however, this observation was not evident in the other cohorts and is likely to be a chance finding. There was no association between the −23Hph1 genotype and adult head circumference, BMI, or body fat within the individual cohorts (Table 1). Similarly, in the pooled analysis adjusted for study, we found no substantial association between −23Hph1 genotype and indexes of growth and obesity (Table 2). In additional analyses, we found no consistent association between the −23Hph1 genotype and weight change (data not shown). However, we found marginal evidence for association between the −23Hph1 genotype and percentage body fat in the pooled analysis (P = 0.04). Nevertheless, based on the number of statistical tests, this association may be a false-positive (13). Exclusion of outliers (greater than three SDs for any variables) did not materially alter these findings (data not shown).
DISCUSSION
Our study does not support previous reports suggesting that INS VNTR class is associated with indexes of body size and obesity or related metabolic traits. By contrast, these data suggest that INS VNTR class genotype is not associated with indexes of body size and adult obesity in middle-aged men and women.
Variation in insulin expression influences growth, development, and metabolism. Phenotypes of animal models of insulin gene deletion exhibit reduced viability, growth, and metabolic control (1). Circulating maternal and fetal insulin is also a positive regulator of fetal growth and birth size in humans (14). However, the mechanistic processes underlying any association between INS VNTR class and indexes of growth or body size and metabolic regulation are speculative. In vitro studies have suggested that class I allele carriers may have higher insulin expression (15), but in vivo studies have found inconsistent associations between INS VNTR class and indexes of circulating levels of insulin (4,6,16). Moreover, whereas higher cord blood IGF-II levels, a regulator of fetal growth and development, have been associated with class III homozygote children (4), in vitro data have found no association (17) or have suggested the class III allele is associated with lower IGF-II transcription (5). In the present study, we found no association between INS VNTR class and circulating levels of insulin or IGF-II.
The INS VNTR class I allele may be overrepresented in obese children (18). Investigations of obese subpopulations also suggest that children carrying a paternal class I allele may have a greater rate of weight gain and risk of obesity (18,19). However, this imprinting effect on body weight regulation has not been replicated in studies of children from unselected populations (7) nor in studies of overweight and obese Caucasian and Pima adults (20,21). Nevertheless in unrelated individuals, one study of children and adolescents reported that class III homozygotes may have lower BMI compared with class I carriers (7). Conversely, it has also been shown that class III homozygote children and young adults may have a greater BMI than class I carriers (4,6). By contrast, we found no association between INS VNTR class as inferred by the −23Hph1 genotype and BMI in 7,999 middle-aged men and women, confirming previous reports based on middle-aged study participants (7,22).
One contemporary birth cohort has reported that INS VNTR class is associated with head circumference at birth and in early childhood (23), which has been corroborated in a larger subgroup of this cohort (4). These studies suggest that class III homozygotes may have larger head circumferences at birth and at 7 years. However, a subsequent study of 5,091 participants failed to find any association between the INS VNTR class and head circumference at 1 year (8). Similarly, in our analysis based on 7,149 study participants, we found no evidence for an association between adult head circumference and INS VNTR class.
We did not find a main genotype effect on birth weight in 3,473 study participants, which is consistent with several earlier studies (6,8,23). However, a study of Pima Indians reported that class III carriers may have lower birth weights (21). This observation is at variance with early results from the Avon Longitudinal Study of Parents and Children cohort, which suggested that, compared with class I homozygotes, class III carriers may have relatively higher birth weights in a subgroup of 348 “nonchangers” (23). However, this observation was not replicated in an additional subcohort of the study (4). Subsequent studies of a U.K. population cohort series and a Finnish birth cohort comprising a total of 2,795 participants have also failed to replicate the association between genotype and birth weight among nonchangers (6,8). Collectively, these observations suggest that the INS VNTR class genotype is unlikely to be a principal regulator of birth weight.
One limitation of this investigation was reliance on self-reported birth weight. Random misclassification of birth weight may lead to attenuation of the association among genotype and indexes of early growth. However, self-reports correlate well with birth records and mothers’ reports of birth weight (24). Birth weight in this cohort also showed the expected correlations with adult anthropometry. Moreover, our self-reported birth weight data showed the well-established inverse associations with measures of blood pressure and glucose tolerance (M.S.S., N.J.W., unpublished observations). The high concordance between genotyping methods suggests that it is unlikely that genotyping error could account for the lack of association between INS VNTR class genotype and indexes of body size and adult obesity.
We did not have data relating to other birth characteristics such as length of gestation, birth order, or postnatal growth realignment that might obscure or modify any genotype–growth/obesity index relation. However, to confound the association between INS VNTR and indexes of body size and obesity, these birth characteristics would need to be associated with both the INS VNTR class genotype and the specific outcome, for which there is no evidence (4,8). Associations between INS VNTR and indexes of body size and obesity have been found largely in contemporary birth cohorts. It is possible that these associations characterize the effect of INS VNTR class and changes in intergenerational environmental factors on growth, body size, and weight (4,6). Yet, effect-modification between INS VNTR and indexes of body size and obesity has not been reliably reproduced, and the biological basis for these possible interactions is unclear (6). Moreover, given the size of the present study, any association between INS VNTR class and indexes of body size and adult obesity in specific subgroups is unlikely to be very strong in this population.
In conclusion, in this population cohort series of 7,999 study participants, we found no evidence to suggest that INS VNTR class genotype is an important determinant of body size and adult obesity in middle-aged men and women.
Demographic and anthropometric characteristics of study participants by −23Hph1 genotype and study population
. | −23Hph1 genotype . | . | . | P value . | . | |||
---|---|---|---|---|---|---|---|---|
. | AA (I/I) . | TA (III/I) . | TT (III/III) . | PR (1 df) . | PA (1 df) . | |||
EPIC1 | ||||||||
n (%) | 1,747 (51) | 1,407 (41) | 292 (8) | |||||
Women (%) | 1,108 (63) | 939 (67) | 202 (69) | 0.14 | 0.02 | |||
Age (years) | 64.1 ± 0.2 | 64.0 ± 0.2 | 64.8 ± 0.5 | 0.64 | 0.50 | |||
BMI (kg/m2) | 26.8 ± 0.1 | 26.6 ± 0.1 | 26.4 ± 0.2 | 0.66 | 0.08 | |||
Body fat (%) | 34.3 ± 0.3 | 34.2 ± 0.3 | 33.8 ± 0.7 | 0.79 | 0.61 | |||
Birth weight (kg) | 3.36 ± 0.03 | 3.34 ± 0.03 | 3.33 ± 0.07 | 0.41 | 0.63 | |||
Head circumference (cm) | 56.1 ± 0.1 | 56.0 ± 0.1 | 56.0 ± 0.1 | 0.75 | 0.17 | |||
EPIC2 | ||||||||
n (%) | 1,942 (52) | 1,471 (40) | 307 (8) | |||||
Women (%) | 963 (50) | 709 (48) | 155 (50) | 0.61 | 0.75 | |||
Age (years) | 59.3 ± 0.2 | 59.6 ± 0.3 | 59.2 ± 0.5 | 0.13 | 0.73 | |||
BMI (kg/m2) | 26.5 ± 0.1 | 26.5 ± 0.1 | 26.6 ± 0.2 | 0.16 | 0.50 | |||
Body fat (%) | 31.4 ± 0.3 | 31.1 ± 0.3 | 31.4 ± 0.6 | 0.55 | 0.75 | |||
Birth weight (kg) | 3.37 ± 0.02 | 3.39 ± 0.03 | 3.32 ± 0.06 | 0.86 | 0.87 | |||
Head circumference (cm) | 56.5 ± 0.1 | 56.5 ± 0.1 | 56.5 ± 0.1 | 0.41 | 0.54 | |||
Ely | ||||||||
n (%) | 409 (49) | 358 (43) | 66 (8) | |||||
Women (%) | 240 (59) | 206 (58) | 38 (58) | 0.93 | 0.85 | |||
Age (years) | 54.1 ± 0.5 | 55.1 ± 0.6 | 55.1 ± 1.3 | 0.72 | 0.20 | |||
BMI (kg/m2) | 26.6 ± 0.2 | 26.4 ± 0.2 | 27.3 ± 0.5 | 0.18 | 0.59 | |||
Body fat (%) | 32.3 ± 0.5 | 32.0 ± 0.5 | 32.2 ± 1.1 | 0.93 | 0.78 | |||
Birth weight (kg) | 3.25 ± 0.06 | 3.26 ± 0.06 | 3.44 ± 0.14 | 0.19 | 0.38 |
. | −23Hph1 genotype . | . | . | P value . | . | |||
---|---|---|---|---|---|---|---|---|
. | AA (I/I) . | TA (III/I) . | TT (III/III) . | PR (1 df) . | PA (1 df) . | |||
EPIC1 | ||||||||
n (%) | 1,747 (51) | 1,407 (41) | 292 (8) | |||||
Women (%) | 1,108 (63) | 939 (67) | 202 (69) | 0.14 | 0.02 | |||
Age (years) | 64.1 ± 0.2 | 64.0 ± 0.2 | 64.8 ± 0.5 | 0.64 | 0.50 | |||
BMI (kg/m2) | 26.8 ± 0.1 | 26.6 ± 0.1 | 26.4 ± 0.2 | 0.66 | 0.08 | |||
Body fat (%) | 34.3 ± 0.3 | 34.2 ± 0.3 | 33.8 ± 0.7 | 0.79 | 0.61 | |||
Birth weight (kg) | 3.36 ± 0.03 | 3.34 ± 0.03 | 3.33 ± 0.07 | 0.41 | 0.63 | |||
Head circumference (cm) | 56.1 ± 0.1 | 56.0 ± 0.1 | 56.0 ± 0.1 | 0.75 | 0.17 | |||
EPIC2 | ||||||||
n (%) | 1,942 (52) | 1,471 (40) | 307 (8) | |||||
Women (%) | 963 (50) | 709 (48) | 155 (50) | 0.61 | 0.75 | |||
Age (years) | 59.3 ± 0.2 | 59.6 ± 0.3 | 59.2 ± 0.5 | 0.13 | 0.73 | |||
BMI (kg/m2) | 26.5 ± 0.1 | 26.5 ± 0.1 | 26.6 ± 0.2 | 0.16 | 0.50 | |||
Body fat (%) | 31.4 ± 0.3 | 31.1 ± 0.3 | 31.4 ± 0.6 | 0.55 | 0.75 | |||
Birth weight (kg) | 3.37 ± 0.02 | 3.39 ± 0.03 | 3.32 ± 0.06 | 0.86 | 0.87 | |||
Head circumference (cm) | 56.5 ± 0.1 | 56.5 ± 0.1 | 56.5 ± 0.1 | 0.41 | 0.54 | |||
Ely | ||||||||
n (%) | 409 (49) | 358 (43) | 66 (8) | |||||
Women (%) | 240 (59) | 206 (58) | 38 (58) | 0.93 | 0.85 | |||
Age (years) | 54.1 ± 0.5 | 55.1 ± 0.6 | 55.1 ± 1.3 | 0.72 | 0.20 | |||
BMI (kg/m2) | 26.6 ± 0.2 | 26.4 ± 0.2 | 27.3 ± 0.5 | 0.18 | 0.59 | |||
Body fat (%) | 32.3 ± 0.5 | 32.0 ± 0.5 | 32.2 ± 1.1 | 0.93 | 0.78 | |||
Birth weight (kg) | 3.25 ± 0.06 | 3.26 ± 0.06 | 3.44 ± 0.14 | 0.19 | 0.38 |
Data are means ± SE, unless otherwise indicated. PR, P value for recessive genetic model (III/III versus class I carriers); PA, P value for additive genetic model [per T (III) allele].
Pooled analysis of associations among −23Hph1 genotype, insulin expression, and indices of body size and obesity
. | n . | β* ± SE . | PA . | Between mean difference (III/III versus class I carriers) . | PR . |
---|---|---|---|---|---|
BMI (kg/m2) | 7,999 | −0.03 ± 0.07 | 0.65 | −0.02 ± 0.16 | 0.91 |
Body fat (%) | 7,926 | −0.28 ± 0.14 | 0.04 | −0.46 ± 0.32 | 0.16 |
Head circumference (cm) | 7,149 | 0.01 ± 0.03 | 0.69 | 0.04 ± 0.07 | 0.59 |
Birth weight (kg) | 3,473 | −0.01 ± 0.02 | 0.76 | −0.02 ± 0.05 | 0.71 |
0-h insulin (pmol/l)* | 803 | 0.16 ± 2.16 | 0.94 | 3.07 ± 5.08 | 0.55 |
Incremental insulin Δ30 | 793 | −3.42 ± 13.80 | 0.80 | −11.9 ± 32.94 | 0.72 |
IGF-II (ng/ml) | 554 | 16.85 ± 14.88 | 0.26 | −25.9 ± 34.84 | 0.46 |
. | n . | β* ± SE . | PA . | Between mean difference (III/III versus class I carriers) . | PR . |
---|---|---|---|---|---|
BMI (kg/m2) | 7,999 | −0.03 ± 0.07 | 0.65 | −0.02 ± 0.16 | 0.91 |
Body fat (%) | 7,926 | −0.28 ± 0.14 | 0.04 | −0.46 ± 0.32 | 0.16 |
Head circumference (cm) | 7,149 | 0.01 ± 0.03 | 0.69 | 0.04 ± 0.07 | 0.59 |
Birth weight (kg) | 3,473 | −0.01 ± 0.02 | 0.76 | −0.02 ± 0.05 | 0.71 |
0-h insulin (pmol/l)* | 803 | 0.16 ± 2.16 | 0.94 | 3.07 ± 5.08 | 0.55 |
Incremental insulin Δ30 | 793 | −3.42 ± 13.80 | 0.80 | −11.9 ± 32.94 | 0.72 |
IGF-II (ng/ml) | 554 | 16.85 ± 14.88 | 0.26 | −25.9 ± 34.84 | 0.46 |
Regression coefficient is the mean change (± SE) in the trait of interest for each additional T (class III) allele (adjusted for sex and study, where appropriate). PR, P value for recessive genetic model (III/III versus class I carriers); PA, P value for additive genetic model [per T (III) allele].