The role of obesity-associated loci identified in genome wide association studies in the determination of pediatric BMI

doi:10.1038/oby.2009.159

Obesity (Silver Spring). Author manuscript; available in PMC 2010 April 28.

Published in final edited form as:

Obesity (Silver Spring). 2009 December; 17(12): 2254–2257.

Published online 2009 May 28. doi: 10.1038/oby.2009.159

PMCID: PMC2860782

NIHMSID: NIHMS195995

The role of obesity-associated loci identified in genome wide association studies in the determination of pediatric BMI

Jianhua Zhao,¹^* Jonathan P. Bradfield,²^* Mingyao Li,³ Kai Wang,² Haitao Zhang,² Cecilia E. Kim,² Kiran Annaiah,² Joseph T. Glessner,² Kelly Thomas,² Maria Garris,² Edward C. Frackelton,² F. George Otieno,² Julie L. Shaner,² Ryan M. Smith,² Rosetta M. Chiavacci,² Robert I. Berkowitz,^4,⁵ Hakon Hakonarson,^1,^2,⁶^** and Struan F.A. Grant^1,^2,⁶^**

¹Division of Human Genetics, The Children's Hospital of Philadelphia, Philadelphia, Pennsylvania 19104, USA

²Center for Applied Genomics, Abramson Research Center, The Children's Hospital of Philadelphia, Philadelphia, Pennsylvania 19104, USA

³Department of Biostatistics and Epidemiology, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA; ³Department of Pediatrics, University of Pennsylvania School of Medicine, Philadelphia, Pennsylvania 19104, USA

⁴Behavioral Health Center and Department of Child and Adolescent Psychiatry, The Children's Hospital of Philadelphia, Philadelphia PA 19104, USA

⁵Center for Weight and Eating Disorders, Department of Psychiatry, University of Pennsylvania, Philadelphia PA 19104, USA

⁶Department of Pediatrics, University of Pennsylvania, Philadelphia PA 19104, USA

^**To whom correspondence should be addressed. Email: grants/at/chop.edu or ; Email: hakonarson/at/chop.edu

^*These authors contributed equally to this work.

See other articles in PMC that cite the published article.

Abstract

The prevalence of obesity in children and adults in the United States has increased dramatically over the past decade. Besides environmental factors, genetic factors are known to play an important role in the pathogenesis of obesity. A number of genetic determinants of adult BMI have already been established through genome wide association studies. In this study, we examined 25 single nucleotide polymorphisms (SNPs) corresponding to thirteen previously reported genomic loci in 6,078 children with measures of BMI. Fifteen of these SNPs yielded at least nominally significant association to BMI, representing nine different loci including INSIG2, FTO, MC4R, TMEM18, GNPDA2, NEGR1, BDNF, KCTD15 and 1q25. Other loci revealed no evidence for association, namely at MTCH2, SH2B1, 12q13 and 3q27. For the 15 associated variants, the genotype score explained 1.12% of the total variation for BMI z-score. We conclude that among thirteen loci that have been reported to associate with adult BMI, at least nine also contribute to the determination of BMI in childhood as demonstrated by their associations in our pediatric cohort.

In the past three years, thirteen genetic loci have been implicated for BMI from the outcomes of genome wide association studies (GWA) studies primarily in adults. Insulin-induced gene 2 (INSIG2) was the first locus to be reported by this method to have a role in obesity¹ but replication attempts have yielded inconsistent outcomes²^–⁶. The second reported locus, the fat mass- and obesity-associated gene (FTO)⁷ has been more robustly observed by others⁸^–¹¹. Subsequent larger studies have uncovered eleven additional genes¹²^–¹⁴, firstly melanocortin 4 receptor (MC4R) from a multi-center meta-analysis¹², then the GIANT consortium revealed six more genes [transmembrane protein 18 (TMEM18), potassium channel tetramerisation domain containing 15 (KCTD15), glucosamine-6-phosphate deaminase 2 (GNPDA2), SH2B adaptor protein 1 (SH2B1), mitochondrial carrier 2 (MTCH2) and neuronal growth regulator 1 (NEGR1)]¹⁴ , five of which were confirmed in the GWA study reported from Iceland (but not GNPDA2 due to an unavailable proxy SNP), who also uncovered and reported loci on 1q25, 3q27 and 12q13¹³ and verified association with the brain-derived neurotrophic factor (BDNF) gene¹⁵.

In this study we aimed at examining these finding in a large pediatric cohort with BMI measures and to determine the relative impact of these variants in childhood. For this purpose, we leveraged genotyping data from our ongoing GWA study of BMI variation in children. The twenty five SNPs corresponding to the thirteen previously reported obesity loci were investigated with respect to their association to normalized pediatric BMI (Table 1; also Supplementary Table 1 for analyses by age categories).

Table 1

Quantitative association results for the candidate loci in the European American BMI cohort (n=6,078), sorted by chromosomal location.

In summary, fifteen of these SNPs yielded at least nominally significant association to BMI (P < 0.05), representing nine different loci with the same direction of effect as previously reported. Of these nine loci, variants at the FTO locus yielded the strongest association with P<10⁻⁴, namely rs8044769 and rs3751812 (P = 7.26×10⁻⁵ and 9.68×10⁻⁵, respectively); in addition, this locus also yielded association with rs8050136 and rs7190492 (P = 1.40×10⁻⁴ and 0.021, respectively) but not with rs6499640.

With a similar magnitude of association to FTO was TMEM18, with rs2867125 yielding a P = 9.72×10⁻⁵, together with almost as strongly associated SNPs, rs7561317 and rs4854344 (P = 1.02×10⁻⁴ and 1.52×10⁻⁴, respectively). Indeed, it is interesting to note that the two adult cohorts that uncovered TMEM18 in obesity also showed it to be second only in significance to FTO¹³^,¹⁴. The third most significantly associated locus was at GNPDA2, (rs13130484; P = 1.32×10⁻⁴).

Overall, in addition to FTO, TMEM18 and GNPDA2, we found evidence for association at the INSIG2, MC4R, NEGR1, 1q25, BDNF and KCTD15 loci. One could argue that we have carried out multiple testing in our BMI cohort for these previously reported SNPs, albeit at a number of magnitudes less than for a full GWA study. If we were to apply the strictest correction, i.e. the Bonferroni correction based on twenty five SNPs, then FTO, TMEM18, GNPDA2 and MC4R would still be considered significant. This is very much in line with the observations made with the pediatric cohort utilized by Willer et al¹⁴; however the one exception is that we do not observe strong association with KCTD15.

It was also observed that SNPs residing at the 12q13, 3q27, MTCH2 and SH2B1 loci did not reveal any evidence of association with BMI in our pediatric cohort.

The positive results for FTO and MC4R come as no surprise as we have previously reported their association with the CDC-defined 95^th percentile of BMI, i.e. obesity, in our pediatric cohort, but limited to ages 2–18 years old⁸^,¹⁶. In this current study, where we utilized BMI z-score on all children 0–18 years old, it is satisfying that they continue to show association at a similar magnitude.

One of the more notable results is the positive association with INSIG2. This association with pediatric BMI, albeit at just the nominal level, will further add to the debate on the relevance of INSIG2 in BMI determination.

For the loci we did not observe any evidence for association for at all may be due to power issues, but could also indicate that they have a less pronounced role in a pediatric setting. Indeed, many of the newly uncovered genes have been implicated in neurological functions, genes which may be more important in BMI determination in adults rather than in children, where other more direct metabolic genes could play a more important role.

Finally, we investigated the fifteen significant SNPs further by testing for association between BMI Z-score and the genotype score by summing the number of BMI increasing alleles across all these SNPs. The resulting P-value for the genotype score was 2.53×10⁻¹⁶ (Figure 1). The genotype score explains 1.12% of the total variation for BMI z-score. We also tested pair-wise interactions between the fifteen significant SNPs, but none of the interaction effects were significant (Supplementary Table 2), suggesting that these fifteen SNPs act additively on the pediatric BMI z-score. As such, we did observe a cumulative effect but not as striking as reported by the GIANT consortium in their adult cohorts¹⁴.

Figure 1

Scatter plot for association between BMI z-score and the genotype score by summing the number of BMI increasing alleles across all fifteen BMI-associated SNPs.

From this analysis, it is clear that a number of loci previously reported from GWA analyses of adult BMI and / or obesity also play a role in our phenotype of interest. While these recently discovered loci unveil several new biomolecular pathways not previously associated with obesity, it is important to note that these well established genetic associations with obesity explain very little of the genetic risk for this pediatric phenotype, suggesting the existence of additional loci whose number and effect size remain unknown. Once our GWA study is complete, we will have the opportunity to look for other variants in the genome that are associated with BMI in childhood.

RESEARCH METHODS AND PROCEDURES

Study Subjects

All subjects were consecutively recruited from the Greater Philadelphia area from 2006 to 2008 at the Children's Hospital of Philadelphia. Our study cohort consisted of 6,078 children of European ancestry with BMI information. All subjects were biologically unrelated and were aged between 0 and 18 years old. This study was approved by the Institutional Review Board of the Children's Hospital of Philadelphia. Parental informed consent was given for each study participant for both the blood collection and subsequent genotyping.

Genotyping

We performed high throughput genome-wide SNP genotyping using either the Illumina Infinium™ II HumanHap550 or Human 610 BeadChip technology in the same manner as our center has reported previously¹⁷. The SNPs analyzed survived the filtering of the genome wide dataset for SNPs with call rates <95%, minor allele frequency <1%, missing rate per person <2% and Hardy-Weinberg equilibrium P < 10⁻⁵.

Loci described from GWA studies published to date have been found using either the Affymetrix or Illumina platform. In the event a locus was reported using both the Illumina and Affymetrix arrays, we used the SNPs present on the Illumina array. In the event of a signal only being described on the Affymetrix array, we either already had that SNP on our Illumina array or we identified and used the best surrogate SNP available. As such, as rs7566605 at INSIG2 and rs10938397 at GNPDA2 were not available on the BeadChip, we employed perfect surrogates for the association analysis i.e. rs17047697 and rs13130484 respectively (r²=1). With rs7498665 at the SH2B1 locus, rs8049439 and rs4788102 are already in strong LD with this SNP so no additional surrogate was selected.

Analysis

From our database of heights and weights for our multi-dimensional scaling (MDS) determined Caucasians, we eliminated BMI outliers using 2% cutoff for each age category in order to remove measurement error and syndromic forms of obesity. As BMI values vary widely across pediatric age groups, each BMI value was adjusted for age and sex then expressed as a z-score.

We queried the data for the indicated SNPs in our pediatric samples. All statistical analyses were carried out using the software package plink¹⁸. By treating BMI as a quantitative trait, association analysis for each SNP was carried out using linear regression with the SNP included as an independent variable (coded as 0, 1, and 2).

Supplementary Material

Supplementary Tables

Click here to view.^{(1.0M, pdf)}

ACKNOWLEDGEMENTS

We would like to thank all participating subjects and families. Hope Thomas, Kisha Harden, Andrew Hill, Kenya Fain, Crystal Johnson-Honesty, Cynthia Drummond, Shanell Harrison and Sarah Wildrick provided expert assistance with genotyping or data collection and management. We would also like to thank Smari Kristinsson, Larus Arni Hermannsson and Asbjörn Krisbjörnsson of Raförninn ehf for their extensive software design and contribution. This research was financially supported by the Children’s Hospital of Philadelphia.

The study is supported in part by a Research Development Award from the Cotswold Foundation (H.H. & S.F.A.G.) and NIH grant 1R01HD056465-01A1.

REFERENCES

1. Herbert A, Gerry NP, McQueen MB, et al. A common genetic variant is associated with adult and childhood obesity. Science. 2006 Apr 14;312(5771):279–283. [PubMed]

2. Loos RJ, Barroso I, O'Rahilly S, Wareham NJ. Comment on "A common genetic variant is associated with adult and childhood obesity". Science. 2007 Jan 12;315(5809):187. author reply 187. [PMC free article] [PubMed]

3. Dina C, Meyre D, Samson C, et al. Comment on "A common genetic variant is associated with adult and childhood obesity". Science. 2007 Jan 12;315(5809):187. author reply 187. [PubMed]

4. Rosskopf D, Bornhorst A, Rimmbach C, et al. Comment on "A common genetic variant is associated with adult and childhood obesity". Science. 2007 Jan 12;315(5809):187. author reply 187. [PubMed]

5. Lyon HN, Emilsson V, Hinney A, et al. The association of a SNP upstream of INSIG2 with body mass index is reproduced in several but not all cohorts. PLoS Genet. 2007 Apr 27;3(4):e61. [PMC free article] [PubMed]

6. Hotta K, Nakamura M, Nakata Y, et al. INSIG2 gene rs7566605 polymorphism is associated with severe obesity in Japanese. J Hum Genet. 2008;53(9):857–862. [PMC free article] [PubMed]

7. Frayling TM, Timpson NJ, Weedon MN, et al. A common variant in the FTO gene is associated with body mass index and predisposes to childhood and adult obesity. Science. 2007 May 11;316(5826):889–894. [PMC free article] [PubMed]

8. Grant SF, Li M, Bradfield JP, et al. Association analysis of the FTO gene with obesity in children of Caucasian and African ancestry reveals a common tagging SNP. PLoS ONE. 2008;3(3):e1746. [PMC free article] [PubMed]

9. Hinney A, Nguyen TT, Scherag A, et al. Genome Wide Association (GWA) Study for Early Onset Extreme Obesity Supports the Role of Fat Mass and Obesity Associated Gene (FTO) Variants. PLoS ONE. 2007;2(12):e1361. [PMC free article] [PubMed]

10. Dina C, Meyre D, Gallina S, et al. Variation in FTO contributes to childhood obesity and severe adult obesity. Nat Genet. 2007 Jun;39(6):724–726. [PubMed]

11. Scuteri A, Sanna S, Chen WM, et al. Genome-Wide Association Scan Shows Genetic Variants in the FTO Gene Are Associated with Obesity-Related Traits. PLoS Genet. 2007 Jul 20;3(7):e115. [PubMed]

12. Loos RJ, Lindgren CM, Li S, et al. Common variants near MC4R are associated with fat mass, weight and risk of obesity. Nat Genet. 2008 Jun;40(6):768–775. [PMC free article] [PubMed]

13. Thorleifsson G, Walters GB, Gudbjartsson DF, et al. Genome-wide association yields new sequence variants at seven loci that associate with measures of obesity. Nat Genet. 2009 Jan;41(1):18–24. [PubMed]

14. Willer CJ, Speliotes EK, Loos RJ, et al. Six new loci associated with body mass index highlight a neuronal influence on body weight regulation. Nat Genet. 2009 Jan;41(1):25–34. [PMC free article] [PubMed]

15. Gunstad J, Schofield P, Paul RH, et al. BDNF Val66Met polymorphism is associated with body mass index in healthy adults. Neuropsychobiology. 2006;53(3):153–156. [PubMed]

16. Grant SF, Bradfield JP, Zhang H, et al. Investigation of the Locus Near MC4R With Childhood Obesity in Americans of European and African Ancestry. Obesity (Silver Spring) 2009 Mar 5; [PMC free article] [PubMed]

17. Hakonarson H, Grant SFA, Bradfield JP, et al. A genome-wide association study identifies KIAA0350 as a type 1 diabetes gene. Nature. 2007;448(7153):591–594. [PubMed]

18. Purcell S, Neale B, Todd-Brown K, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007 Sep;81(3):559–575. [PubMed]