- Research article
- Open Access
No evidence for association between BMI and 10 candidate genes at ages 4, 7 and 10 in a large UK sample of twins
BMC Medical Genetics volume 9, Article number: 12 (2008)
Over the last decade, associations between Body Mass Index (BMI) and a variety of candidate genes have been reported, but samples have almost all been limited to adults. The purpose of the present study was to test the developmental origins of some of these associations in a large longitudinal sample of children.
For 10 single-nucleotide polymorphisms (SNPs) in candidate genes reported to be associated with BMI in adults, we examined associations with BMI in a sample of 5000 children (2500 twin pairs) with BMI data at 4, 7 and 10 years. Association analyses were performed using the Quantitative Transmission Disequilibrium Test and we corrected for multiple testing using the False Discovery Rate.
Despite having 80% power to detect associations that account for as little as 0.2% of the variance of BMI, none of the 10 SNPs were significantly associated with BMI at any age, although two SNPs showed trends in the expected direction.
The lack of association for these ten previously reported associations, despite our large sample size, is typical of associations between candidate genes and complex traits. However, some of the reported SNP associations with BMI might emerge as we continue to follow the sample into adolescence and adulthood. This report highlights the importance of developmentally appropriate candidate genes.
The rise in obesity is not only seen in the adult population – obesity is also becoming more prevalent in childhood and adolescence [see e.g. ]. It is possible that obesity that starts in childhood results in additional health and psychological consequences [2, 3]. Twin research, across the lifespan, consistently points to substantial genetic influence on individual differences in BMI and obesity [4–6]. Due to the rising levels of overweight and obesity in childhood, and the obvious long-term health problems and costs that early-onset and persistent obesity will cause, there has been an increased interest into the developmental patterns of BMI in childhood and adolescence. Very little is known about the developmental etiology of BMI and obesity in early and middle childhood. Recent research into the prevalence of obesity and overweight in an adolescent sample has demonstrated that weight is a stable phenotype by this stage of development, with little movement between groups of children classified as normal weight or overweight at age 11 and then followed up for five years . Since persistent obesity is established before the age of 11, research into BMI and obesity should be targeted at younger children to assess the developmental patterns and etiology of BMI and obesity .
Over the last decade associations between Body Mass Index (BMI) and a variety of candidate genes have been reported , but samples have almost all been limited to adults. When genes are identified that are associated with complex traits such as obesity, one direction for top-down 'behavioral genomic' research  is to explore the developmental origins of the associations such as how early in the lifespan the associations appear. The purpose of the present study was to test some of these associations in a large longitudinal sample of children. We investigated associations with BMI in a longitudinal sample of 5000 children assessed at 4, 7 and 10 years for ten polymorphisms previously associated with obesity ; (see Table 1). Because genes largely contribute to stability for complex traits , we predicted that genes associated with BMI in adulthood would show associations with BMI in childhood, even as early as 4 years of age. The sample size provides 80% power to detect associations accounting for as little as 0.2% of the variance (p = .01; one-tailed because we only accept results in the same direction as the original reports) .
The sample was the Twins Early Development Study (TEDS), a study of twins born in the UK between 1994 – 1996 . The TEDS sample has been shown to be reasonably representative of the UK population [12, 13]. All twins and parents involved in TEDS provide informed consent for each stage of the study. Ethical approval for the Twins Early Development Study has been provided by the King's College London ethics committee, reference number: 05/Q0706/228. Response rates based on active families were 65%, 63%, and 65%, respectively at the three ages. For the purposes of the current study, we excluded from the analyses families in which at least one member of the twin pair had a specific medical syndrome or was an extreme outlier for perinatal problems such as extreme low birth weight. Data was included from monozygotic (MZ) male and female twins; dizygotic (DZ) male and female twins, and opposite-sex DZ twins. Across all genotypes the sample was 42% MZ (MZ female: 23%; MZ male: 19%), 32% DZ (DZ female: 17%; DZ male: 15%), and 26% opposite-sex DZ.
The mean ages of the twins when questionnaires were returned were 4.05 (s.d. = .16), 7.05 (s.d. = .25), and 9.93 (s.d. = .87) respectively. Zygosity was assessed through a parent questionnaire of physical similarity, which has been shown to be over 95% accurate when compared to DNA testing . For cases where zygosity was unclear from this questionnaire, DNA testing was conducted.
Height and weight data were collected by postal questionnaire from the parents when their children were 4, 7 and 10 years of age. Correlations with measured heights and weights in a sub-sample of these children at age 11 were 0.83 and 0.90. The sample included 5000 children at each age. BMI was calculated from the height and weight data (BMI = weight (kg)/(height (m)2)) and converted to BMI z-scores. BMI z-scores take into consideration the child's age and sex. They were based on 1990 UK growth reference curves  and were calculated using the program lmsGrowth (available from http://homepage.mac.com/tjcole). Children were categorized as normal weight, overweight or obese based on International Obesity Task Force (IOTF) criteria ; which are based on age-specific and sex-specific cut-off points linked to growth curves, and correspond to BMI criteria for normal weight, overweight and obesity in adults . Although there are disagreements about the correct cut-off for childhood obesity [see for example [17–19]] we opted for the IOTF criteria because it was important that the cut-offs used in this study were internationally meaningful, in addition, we are not using the cut-offs for clinical purposes. Finally, the focus of our paper is not on the categorical analyses, but on the more powerful continuous analyses.
The percentages of the sample in each category (normal weight, overweight and obese) were 83%, 12% and 5% at 4 years; 88%, 9% and 3% at 7 years; and 86%, 11% and 3% at 10 years.
The mean raw BMI scores (i.e. unstandardized for age and sex) were 15.83 (sd = 1.96) at 4 years; 15.73 (sd = 1.99) at 7 years; and 17.24 (sd = 2.93) at 10 years. The raw BMIs of the 3 groups (normal weight, overweight and obese) at each age were: 4 years: 15.19 (sd = 1.29); 18.06 (sd = 0.53); and 20.95 (sd = 1.81); 7 years: 15.22 (sd = 1.29); 18.92 (sd = 0.76); and 22.53 (sd = 1.96); and 10 years: 16.36 (sd = 1.68); 21.60 (sd = 1.41); and 27.41 (sd = 5.07).
Candidate SNP selection and genotyping
Candidate SNPs were selected based on the 2005 update of the human obesity gene map . Nine SNPs were selected on the following criteria: i) the association must have been shown in a sample with at least 500 participants, ii) the measured phenotype had to be BMI, waist circumference or obesity, and iii) the SNP had to have a minor allele frequency greater than 0.1 because statistical power to detect associations is too low for alleles with lower frequencies. The SNPs that were selected are shown in Table 1, which also includes information about the genes in which they reside, the risk allele (that is, which of the two alleles is associated with higher BMI), and the phenotype of interest. Genotyping of these nine SNPs was performed by Kbiosciences, UK  using the competitive allele specific PCR (KASPar) method. The error rate of this method is low (< 0.5% in our analyses of blind duplicates).
In addition, SNP rs7566605 was genotyped in an attempt to replicate a recent report of an association between rs7566605 and obesity in four out of five populations , although it subsequently failed to replicate in several studies . Genotyping for rs7566605 was performed in-house using Taqman on a subsample of 3000 children. The overall genotyping call rate was 93%, and we only called the genotypes if the TaqMan quality value exceeded 95. The allele frequencies for all 10 candidate SNPs are included in Table 1; all the frequencies in our sample are similar to those for Europeans in HapMap.
Results were analyzed using the QTDT total association model , which takes into account the twin pair structure by comparing variance between and within pairs. We did not have parental genotypes, so were unable to include IBD status in the model. The total association model was chosen after we had verified that there was no significant population stratification. In the QTDT analyses we specifically modeled a polygenic variance component, as well as a non-shared environment and twin-specific environmental variance components. Variance components were tested for significance before being included in the model. The heritability of BMI estimated by the QTDT analyses was 41% at 4 years, 60% at 7 years and 75% at 10 years.
Our sample of 5000 (2500 twin pairs) provides 80% power to detect an effect size (r2) of 0.2% (p = .01, one-tailed, as calculated using ). We corrected for multiple testing using the False Discovery Rate . This multiple testing correction was applied across all ages (i.e. 30 tests) rather than at each age separately, because the measures of BMI at each age are not independent. We analyzed the entire sample of 5000 in order to maximize power to detect associations of small effect size. In addition, because the original reports compared obese and overweight cases versus controls, we also divided the sample into groups (obese, overweight, and normal) and compared genotypic frequency differences between the groups using chi-square analysis. We compared the obese group to the normal weight group and also the obese and overweight groups (combined) were compared to the normal-weight group. These analyses were done on one randomly selected member of each twin pair.
A specific limitation of the present study is the use of parental reports of children's heights and weights, although previous research has shown that they are reasonably reliable . This was confirmed in our sample on a sub-sample at 11 years: Home measurements of height and weight, taken by trained researchers, correlated .83 and .90, respectively, with parental report. A further limitation is that we were only able to study BMI. It would have been beneficial to include other measures of adiposity, but this was not feasible in this large sample. In addition BMI was measured only at ages 4, 7 and 10 years. It is possible that we could have got different results if we had more frequent measurements of BMI during childhood, so that we could more accurately estimate when in development the association with the candidate gene emerged. However, as we found no significant associations, this limitation is unlikely to impact on the interpretation of these results.
Finally, this analysis does not consider genetic influences that impact on change during development. This could have been done using growth curve analyses, similar to Podolsky et al. . However, growth curve analyses are best performed with at least four measurements, so we were unable to include such analyses in this study. All of the results should be interpreted in light of these limitations.
Results and Discussion
We verified that none of the SNPs deviated significantly from Hardy-Weinberg Equilibrium using the exact test in Pedstats . Table 2 shows the mean (and SD) of BMI by genotype for each of the 10 SNPs. Table 3 presents results from the Quantitative Transmission Disequilibrium Test (QTDT)  total association test. Of 30 comparisons, only one yielded a p-value less than .05. After correcting the data for multiple testing using a False Discovery Rate (FDR) of 0.05 , none of the associations were significant.
The only notable genotypic differences in Table 2 involved SNPs rs1801282 (at ages 7 and 10) and rs1800206 (at ages 4, 7 and 10), which demonstrated trends in the expected direction for homozygotes for the risk allele. SNP rs1801282 showed 0.3SD greater BMI for children homozygous for the risk allele at ages 7 and 10, although not at age 4 and SNP rs1800206 showed a minimum 0.5SD greater BMI for children homozygous for the risk allele. These differences are nominally significant in comparisons between children homozygous for the risk allele and the rest of the sample (p values = 0.001 – 0.028). However, although they are in the direction expected on the basis of the previous reports, these associations could not survive a false discovery rate test. Also, because the risk alleles for these two SNPs have low frequencies (average 12% and 7%, respectively), these two genotypes are among those with the smallest sample size and thus have the largest standard errors.
We also ran the analyses separately for males and females and found no gender-specific associations, although with only half the sample we had reduced power to detect such associations. Results from these analyses can be found in Table 4. We present the p-value from the gender-specific QTDT analyses (performed in the same way as the analyses on the full sample). Before correction for multiple testing, there were two nominally significant associations for the female analyses (rs5082 at 4 years and rs6659176 at 10 years). These significant values did not survive the correction for multiple testing. There were no significant associations in the male analyses.
Because the original reports of association involved comparisons between cases and controls rather than continuous variation in BMI throughout the distribution, we also compared genotypic frequencies for obese versus normal-weight and for obese + overweight (combined) versus normal-weight, even though these are less powerful analyses. Chi-square analyses indicated no significant genotypic frequency differences after correcting for multiple testing using FDR. There were two significant p-values (< .05) in each set of 30 chi-square analyses. The p-values (uncorrected for multiple testing) ranged from .019–.932 for the analysis of obese+overweight versus controls (average p-value = .46), and from .019–1.00 for the comparison between obese and normal weight individuals (average p-value = .53). However, the two SNPs mentioned above again showed trends in the expected direction with the risk alleles showing greater genotypic frequencies in the obese and overweight groups.
These negative results reflect the larger issue that candidate gene associations have proven difficult to replicate for complex traits  including BMI and obesity . It is likely that some of the non-replication of candidate gene associations is due to studies being underpowered to detect associations of small effect . But studies with large sample sizes have also failed to replicate obesity SNPs [for example ], in particular SNP rs7566605 . Replication in multiple samples is essential (e.g., the recent FTO gene study ) and may soon be a requirement by journals before publication.
Although the lack of an association between these ten candidate genes and BMI in our sample is typical of association studies for complex traits, it also highlights the importance of developmental candidate genes. We cannot assume that the same genetic variants are influencing complex traits, such as BMI, throughout the developmental process. Childhood BMI is marked by change, and thus may not correspond to adult BMI, which is relatively stable over time. There is the possibility that positive associations for these candidate genes may emerge later in development as we follow these children through adolescence and into adulthood.
Baker JL, Olsen LW, Sorensen TIA: Childhood Body-Mass Index and the Risk of Coronary Heart Disease in Adulthood. New England Journal of Medicine. 2007, 357: 2329-2337. 10.1056/NEJMoa072515.
Wardle J: Understanding the aetiology of childhood obesity: implications for treatment. Proc Nutr Soc. 2005, 64: 73-79. 10.1079/PNS2004411.
Whitaker RC, Wright JA, Pepe MS, Seidel KD, Dietz WH: Predicting Obesity in Young Adulthood from Childhood and Parental Obesity. New England Journal of Medicine. 1997, 337: 869-10.1056/NEJM199709253371301.
Haworth CMA, Plomin R, Carnell S, Wardle J: Childhood Obesity: Genetic and Environmental Overlap with Normal-Range BMI. Obesity.
Koeppen-Schomerus G, Wardle J, Plomin R: A genetic analysis of weight and overweight in 4-year-old twin pairs. Int J Obes Relat Metab Disord. 2001, 25: 838-844. 10.1038/sj.ijo.0801589.
Maes HH, Neale MC, Eaves LJ: Genetic and environmental factors in relative body weight and human adiposity. Behav Genet. 1997, 27: 325-351. 10.1023/A:1025635913927.
Wardle J, Brodersen NH, Cole TJ, Jarvis MJ, Boniface DR: Development of adiposity in adolescence: five year longitudinal study of an ethnically and socioeconomically diverse sample of young people in Britain. BMJ. 2006, 332: 1130-1135. 10.1136/bmj.38807.594792.AE.
Rankinen T, Zuberi A, Chagnon YC, Weisnagel SJ, Argyropoulos G, Walts B, Perusse L, Bouchard C: The human obesity gene map: the 2005 update. Obesity (Silver Spring). 2006, 14: 529-644.
Plomin R, DeFries JC, Craig IW, McGuffin P: Behavioral genetics in the postgenomic era. 2003, Washington, DC: American Psychological Association
Cardon LR: Genetic influence on body mass index in early childhood. Behavior genetic approaches in behavioral medicine. Edited by: Turner JR, Cardon LR, Hewitt JK. 1995, New York: Plenum, 133-143.
Cohen J: Statistical power analysis for the behavioral sciences. 1988, Hillsdale, New Jersey: Lawrence Erlbaum Associates, 2
Oliver BR, Plomin R: Twins Early Development Study (TEDS): A multivariate, longitudinal genetic investigation of language, cognition and behavior problems from childhood through adolescence. Twin Research and Human Genetics. 2007, 10: 96-105. 10.1375/twin.10.1.96.
Kovas Y, Haworth CMA, Dale PS, Plomin R: The genetic and environmental origins of learning abilities and disabilities in the early school years. Monogr Soc Res Child Dev. 2007, 72 (3): vii-160. 10.1111/j.1540-5834.2007.00453.x.
Price TS, Freeman B, Craig IW, Petrill SA, Ebersole L, Plomin R: Infant zygosity can be assigned by parental report questionnaire data. Twin Research. 2000, 3: 129-133. 10.1375/136905200320565391.
Cole TJ, Freeman JV, Preece MA: Body mass index reference curves for the UK, 1990. Arch Dis Child. 1995, 73: 25-29.
Cole TJ, Bellizzi MC, Flegal KM, Dietz WH: Establishing a standard definition for child overweight and obesity worldwide: international survey. BMJ. 2000, 320: 1240-1245. 10.1136/bmj.320.7244.1240.
Chinn S: Definitions of childhood obesity: current practice. European Journal of Clinical Nutrition. 2006, 60: 1189-1194. 10.1038/sj.ejcn.1602436.
Reilly JJ: Assessment of Childhood Obesity: National Reference Data or International Approach?. Obesity Research. 2002, 10: 838-840.
Reilly JJ: Diagnostic accuracy of the BMI for age in paediatrics. Int J Obes (Lond). 2006, 30: 595-597. 10.1038/sj.ijo.0803301.
Kbiosciences, UK. [http://www.kbioscience.co.uk]
Herbert A, Gerry NP, McQueen MB, Heid IM, Pfeufer A, Illig T, Wichmann HE, Meitinger T, Hunter D, Hu FB, Colditz G, Hinney A, Hebebrand J, Koberwitz K, Zhu X, Cooper R, Ardlie K, Lyon H, Hirschhorn JN, Laird NM, Lenburg ME, Lange C, Christman MF: A common genetic variant is associated with adult and childhood obesity. Science. 2006, 312: 279-283. 10.1126/science.1124779.
Herbert A, Gerry NP, McQueen MB, Heid IM, Pfeufer A, Illig T, Wichmann HE, Meitinger T, Hunter D, Hu FB, Colditz G, Hinney A, Hebebrand J, Koberwitz K, Zhu X, Cooper R, Ardlie K, Lyon H, Hirschhorn JN, Laird NM, Lenburg ME, Lange C, Christman MF: Response to comments on "A common genetic variant is associated with adult and childhood obesity". Science. 2007, 315: 187e-10.1126/science.1129763.
Abecasis GR, Cardon LR, Cookson WO: A general test of association for quantitative traits in nuclear families. Am J Hum Genet. 2000, 66: 279-292. 10.1086/302698.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society B. 1995, 57: 289-300.
Reed DR, Price RA: Estimates of the heights and weights of family members: accuracy of informant reports. Int J Obes Relat Metab Disord. 1998, 22: 827-835. 10.1038/sj.ijo.0800666.
Podolsky RH, Barbeau P, Kang HS, Zhu H, Treiber FA, Snieder H: Candidate genes and growth curves for adiposity in African-and European-American youth. International Journal of Obesity. 2007, 31: 1491-1499. 10.1038/sj.ijo.0803673.
Wigginton JE, Cutler DJ, Abecasis GR: A note on exact tests of Hardy-Weinberg equilibrium. Am J Hum Genet. 2005, 76: 887-893. 10.1086/429864.
Lohmueller KE, Pearce CL, Pike M, Lander ES, Hirschhorn JN: Meta-analysis of genetic association studies supports a contribution of common variants to susceptibility to common disease. Nat Genet. 2003, 33: 177-182. 10.1038/ng1071.
Lyon HN, Hirschhorn JN: Genetics of common forms of obesity: a brief overview. Am J Clin Nutr. 2005, 82: 215S-217S.
Cardon LR, Bell JI: Association study designs for complex diseases. Nat Rev Genet. 2001, 2: 91-99. 10.1038/35052543.
Lyon HN, Florez JC, Bersaglieri T, Saxena R, Winckler W, Almgren P, Lindblad U, Tuomi T, Gaudet D, Zhu X, Cooper R, Ardlie KG, Daly MJ, Altshuler D, Groop L, Hirschhorn JN: Common variants in the ENPP1 gene are not reproducibly associated with diabetes or obesity. Diabetes. 2006, 55: 3180-3184. 10.2337/db06-0407.
Frayling TM, Timpson NJ, Weedon MN, Zeggini E, Freathy RM, Lindgren CM, Perry JRB, Elliott KS, Lango H, Rayner NW, Shields B, Harries LW, Barrett JC, Ellard S, Groves CJ, Knight B, Patch A, Ness AR, Ebrahim S, Lawlor DA, Ring SM, Ben-Shlomo Y, Jarvelin MR, Sovio U, Bennett AJ, Melzer D, Ferrucci L, Loos RJF, Barroso I, Wareham NJ, Karpe F, Owen KR, Cardon LR, Walker M, Hitman GA, Palmer CNA, Doney ASF, Morris AD, Smith GD, The Wellcome Trust Case Control Consortium, Hattersley AT, McCarthy MI: A Common Variant in the FTO Gene Is Associated with Body Mass Index and Predisposes to Childhood and Adult Obesity. Science. 2007, 316: 889-894. 10.1126/science.1141634.
Robitaille J, Despres JP, Perusse L, Vohl MC: The PPAR-gamma P12A polymorphism modulates the relationship between dietary fat intake and components of the metabolic syndrome: results from the Quebec Family Study. Clin Genet. 2003, 63: 109-116. 10.1034/j.1399-0004.2003.00026.x.
van 't Hooft FM, Ruotolo G, Boquist S, de FU, Eggertsen G, Hamsten A: Human evidence that the apolipoprotein a-II gene is implicated in visceral fat accumulation and metabolism of triglyceride-rich lipoproteins. Circulation. 2001, 104: 1223-1228. 10.1161/hc3601.095709.
Meirhaeghe A, Helbecque N, Cottel D, Amouyel P: Beta2-adrenoceptor gene polymorphism, body weight, and physical activity. Lancet. 1999, 353: 896-10.1016/S0140-6736(99)00251-2.
Bosse Y, Despres JP, Bouchard C, Perusse L, Vohl MC: The peroxisome proliferator-activated receptor alpha L162V mutation is associated with reduced adiposity. Obes Res. 2003, 11: 809-816.
Hung CC, Farooqi IS, Ong K, Luan J, Keogh JM, Pembrey M, Yeo GS, Dunger D, Wareham NJ, O'Rahilly S: Contribution of variants in the small heterodimer partner gene to birthweight, adiposity, and insulin levels: mutational analysis and association studies in multiple populations. Diabetes. 2003, 52: 1288-1291. 10.2337/diabetes.52.5.1288.
Brand E, Wang JG, Herrmann SM, Staessen JA: An epidemiological study of blood pressure and metabolic phenotypes in relation to the Gbeta3 C825T polymorphism. J Hypertens. 2003, 21: 729-737. 10.1097/00004872-200304000-00016.
Meyre D, Bouatia-Naji N, Tounian A, Samson C, Lecoeur C, Vatin V, Ghoussaini M, Wachter C, Hercberg S, Charpentier G, Patsch W, Pattou F, Charles MA, Tounian P, Clement K, Jouret B, Weill J, Maddux BA, Goldfine ID, Walley A, Boutin P, Dina C, Froguel P: Variants of ENPP1 are associated with childhood and adult obesity and increase the risk of glucose intolerance and type 2 diabetes. Nat Genet. 2005, 37: 863-867. 10.1038/ng1604.
Corella D, Guillen M, Portoles O, Sorli JV, Alonso V, Folch J, Saiz C: Gender specific associations of the Trp64Arg mutation in the beta3-adrenergic receptor gene with obesity-related phenotypes in a Mediterranean population: interaction with a common lipoprotein lipase gene variation. J Intern Med. 2001, 250: 348-360. 10.1046/j.1365-2796.2001.00894.x.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/9/12/prepub
We gratefully acknowledge the ongoing contribution of the parents and children in the Twins' Early Development Study (TEDS). TEDS is supported by a program grant (G0500079) from the U.K. Medical Research Council and the work on obesity in TEDS is supported in part by a grant from the Biotechnology and Biological Sciences Research Council (31/D19086).
The author(s) declare that they have no competing interests.
CMAH performed the statistical analyses and drafted the manuscript. LMB selected the candidate SNPs, and organized the genotyping by Kbiosciences. SJD participated in the genotyping for the study. JW and RP conceived and designed the study. All authors contributed to the final critical revision of the manuscript.
About this article
Cite this article
Haworth, C.M., Butcher, L.M., Docherty, S.J. et al. No evidence for association between BMI and 10 candidate genes at ages 4, 7 and 10 in a large UK sample of twins. BMC Med Genet 9, 12 (2008). https://doi.org/10.1186/1471-2350-9-12
- Body Mass Index
- Risk Allele
- Twin Pair
- International Obesity Task Force
- Quantitative Transmission Disequilibrium Test