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,68). 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.

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).

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).

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.

TABLE 1

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/m226.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/m226.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/m226.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/m226.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/m226.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/m226.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].

TABLE 2

Pooled analysis of associations among −23Hph1 genotype, insulin expression, and indices of body size and obesity

nβ* ± SEPABetween mean difference (III/III versus class I carriers)PR
BMI (kg/m27,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β* ± SEPABetween mean difference (III/III versus class I carriers)PR
BMI (kg/m27,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].

1.
Kitamura T, Kahn CR, Accili D: Insulin receptor knockout mice.
Annu Rev Physiol
65
:
313
–332,
2003
2.
Vafiadis P, Bennett ST, Todd JA, Nadeau J, Grabs R, Goodyer CG, Wickramasinghe S, Colle E, Polychronakos C: Insulin expression in human thymus is modulated by INS VNTR alleles at the IDDM2 locus.
Nat Genet
15
:
289
–292,
1997
3.
Kennedy GC, German MS, Rutter WJ: The minisatellite in the diabetes susceptibility locus IDDM2 regulates insulin transcription.
Nat Genet
9
:
293
–298,
1995
4.
Ong KK, Petry CJ, Barratt BJ, Ring S, Cordell HJ, Wingate DL, Pembrey ME, Todd JA, Dunger DB: Maternal-fetal interactions and birth order influence insulin variable number of tandem repeats allele class associations with head size at birth and childhood weight gain.
Diabetes
53
:
1128
–1133,
2004
5.
Paquette J, Giannoukakis N, Polychronakos C, Vafiadis P, Deal C: The INS 5′ variable number of tandem repeats is associated with IGF2 expression in humans.
J Biol Chem
273
:
14158
–14164,
1998
6.
Mitchell SM, Hattersley AT, Knight B, Turner T, Metcalf BS, Voss LD, Davies D, McCarthy A, Wilkin TJ, Smith GD, Ben Shlomo Y, Frayling TM: Lack of support for a role of the insulin gene variable number of tandem repeats minisatellite (INS-VNTR) locus in fetal growth or type 2 diabetes-related intermediate traits in United Kingdom populations.
J Clin Endocrinol Metab
89
:
310
–317,
2004
7.
Heude B, Dubois S, Charles MA, Deweirder M, Dina C, Borys JM, Ducimetiere P, Froguel P: VNTR polymorphism of the insulin gene and childhood overweight in a general population.
Obes Res
12
:
499
–504,
2004
8.
Bennett AJ, Sovio U, Ruokonen A, Martikainen H, Pouta A, Taponen S, Hartikainen AL, King VJ, Elliott P, Jarvelin MR, McCarthy MI: Variation at the insulin gene VNTR (variable number tandem repeat) polymorphism and early growth: studies in a large Finnish birth cohort.
Diabetes
53
:
2126
–2131,
2004
9.
Day N, Oakes S, Luben R, Khaw KT, Bingham S, Welch A, Wareham N: EPIC-Norfolk: study design and characteristics of the cohort: European Prospective Investigation of Cancer.
Br J Cancer
80 (Suppl. 1)
:
95
–103,
1999
10.
Wareham NJ, Byrne CD, Williams R, Day NE, Hales CN: Fasting proinsulin concentrations predict the development of type 2 diabetes.
Diabetes Care
22
:
262
–270,
1999
11.
Bennett ST, Todd JA: Human type 1 diabetes and the insulin gene: principles of mapping polygenes.
Annu Rev Genet
30
:
343
–370,
1996
12.
Barratt BJ, Payne F, Lowe CE, Hermann R, Healy BC, Harold D, Concannon P, Gharani N, McCarthy MI, Olavesen MG, McCormack R, Guja C, Ionescu-Tirgoviste C, Undlien DE, Ronningen KS, Gillespie KM, Tuomilehto-Wolf E, Tuomilehto J, Bennett ST, Clayton DG, Cordell HJ, Todd JA: Remapping the insulin gene/IDDM2 locus in type 1 diabetes.
Diabetes
53
:
1884
–1889,
2004
13.
Dahlman I, Eaves IA, Kosoy R, Morrison VA, Heward J, Gough SC, Allahabadia A, Franklyn JA, Tuomilehto J, Tuomilehto-Wolf E, Cucca F, Guja C, Ionescu-Tirgoviste C, Stevens H, Carr P, Nutland S, McKinney P, Shield JP, Wang W, Cordell HJ, Walker N, Todd JA, Concannon P: Parameters for reliable results in genetic association studies in common disease.
Nat Genet
30
:
149
–150,
2002
14.
Hattersley AT, Tooke JE: The fetal insulin hypothesis: an alternative explanation of the association of low birthweight with diabetes and vascular disease.
Lancet
353
:
1789
–1792,
1999
15.
Vafiadis P, Bennett ST, Colle E, Grabs R, Goodyer CG, Polychronakos C: Imprinted and genotype-specific expression of genes at the IDDM2 locus in pancreas and leucocytes.
J Autoimmun.
9
:
397
–403,
1996
16.
Ahmed S, Bennett ST, Huxtable SJ, Todd JA, Matthews DR, Gough SC: INS VNTR allelic variation and dynamic insulin secretion in healthy adult non-diabetic Caucasian subjects.
Diabet Med
16
:
910
–917,
1999
17.
Vafiadis P, Bennett ST, Todd JA, Grabs R, Polychronakos C: Divergence between genetic determinants of IGF2 transcription levels in leukocytes and of IDDM2-encoded susceptibility to type 1 diabetes.
J Clin Endocrinol Metab
83
:
2933
–2939,
1998
18.
Le Stunff C, Fallin D, Schork NJ, Bougneres P: The insulin gene VNTR is associated with fasting insulin levels and development of juvenile obesity.
Nat Genet
26
:
444
–446,
2000
19.
Le Stunff C, Fallin D, Bougneres P: Paternal transmission of the very common class I INS VNTR alleles predisposes to childhood obesity.
Nat Genet
29
:
96
–99,
2001
20.
Huxtable SJ, Saker PJ, Haddad L, Walker M, Frayling TM, Levy JC, Hitman GA, O’Rahilly S, Hattersley AT, McCarthy MI: Analysis of parent-offspring trios provides evidence for linkage and association between the insulin gene and type 2 diabetes mediated exclusively through paternally transmitted class III variable number tandem repeat alleles.
Diabetes
49
:
126
–130,
2000
21.
Lindsay RS, Hanson RL, Wiedrich C, Knowler WC, Bennett PH, Baier LJ: The insulin gene variable number tandem repeat class I/III polymorphism is in linkage disequilibrium with birth weight but not type 2 diabetes in the Pima population.
Diabetes
52
:
187
–193,
2003
22.
O’Dell SD, Bujac SR, Miller GJ, Day IN: Associations of IGF2 ApaI RFLP and INS VNTR class I allele size with obesity.
Eur J Hum Genet
7
:
821
–827,
1999
23.
Dunger DB, Ong KK, Huxtable SJ, Sherriff A, Woods KA, Ahmed ML, Golding J, Pembrey ME, Ring S, Bennett ST, Todd JA: Association of the INS VNTR with size at birth: ALSPAC Study Team: Avon Longitudinal Study of Pregnancy and Childhood.
Nat Genet
19
:
98
–100,
1998
24.
Sanderson M, Williams MA, White E, Daling JR, Holt VL, Malone KE, Self SG, Moore DE: Validity and reliability of subject and mother reporting of perinatal factors.
Am J Epidemiol
147
:
136
–140,
1998