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 obesity1
but replication attempts have yielded inconsistent outcomes2–6
. The second reported locus, the fat mass- and obesity-associated gene (FTO
has been more robustly observed by others8–11
. Subsequent larger studies have uncovered eleven additional genes12–14
, firstly melanocortin 4 receptor (MC4R
) from a multi-center meta-analysis12
, 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 12q1313
and verified association with the brain-derived neurotrophic factor (BDNF
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 (; also Supplementary Table 1
for analyses by age categories).
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−4, namely rs8044769 and rs3751812 (P = 7.26×10−5 and 9.68×10−5, respectively); in addition, this locus also yielded association with rs8050136 and rs7190492 (P = 1.40×10−4 and 0.021, respectively) but not with rs6499640.
With a similar magnitude of association to FTO
, with rs2867125 yielding a P
, together with almost as strongly associated SNPs, rs7561317 and rs4854344 (P
, 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 FTO13, 14
. The third most significantly associated locus was at GNPDA2
, (rs13130484; P
Overall, in addition to FTO
, we found evidence for association at the INSIG2
, 1q25, BDNF
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
would still be considered significant. This is very much in line with the observations made with the pediatric cohort utilized by Willer et al14
; 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
come as no surprise as we have previously reported their association with the CDC-defined 95th
percentile of BMI, i.e. obesity, in our pediatric cohort, but limited to ages 2–18 years old8, 16
. 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−16
(). 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 cohorts14
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.