# Novel loci for childhood body mass index and shared heritability with adult cardiometabolic traits

Suzanne Vogelezang, Jonathan Bradfield, Tarunveer Ahluwalia, John Curtin, Timo Lakka, Niels Grarup, Markus Scholz, Peter van Der Most, Claire Monnereau, Evie Stergiakouli, et al.

► To cite this version:

Suzanne Vogelezang, Jonathan Bradfield, Tarunveer Ahluwalia, John Curtin, Timo Lakka, et al.,

Novel loci for childhood body mass index and shared heritability with adult cardiometabolic traits.

PLoS Genetics, 2020, 16 (10), pp.e1008718. 10.1371/journal.pgen.1008718. hal-03269687

HAL Id: hal-03269687

https://hal.inrae.fr/hal-03269687

Submitted on 28 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

Novel loci for childhood body mass index and
shared heritability with adult cardiometabolic
traits

Suzanne Vogelezang<sup>1,2,3‡</sup>, Jonathan P. Bradfield<sup>4,5‡</sup>, Tarunveer S. Ahluwalia<sup>6,7,8,9</sup>, John A. Curtin<sup>10</sup>, Timo A. Lakka<sup>11,12,13</sup>, Niels Grarup<sup>8</sup>, Markus Scholz<sup>14,15</sup>, Peter J. van der Most<sup>16</sup>, Claire Monnereau<sup>1,2,3</sup>, Evie Stergiakouli<sup>17,18,19</sup>, Anni Heiskala<sup>20</sup>, Momoko Horikoshi<sup>21,22,23</sup>, Iryna O. Fedko<sup>24</sup>, Natalia Vilor-Tejedor<sup>25,26,27,28</sup>, Diana L. Cousminer<sup>29,30</sup>, Marie Standl<sup>31</sup>, Carol A. Wang<sup>32</sup>, Jorma Viikari<sup>33,34</sup>, Frank Geller<sup>35</sup>, Carmen Íñiguez<sup>28,36,37</sup>, Niina Pitkänen<sup>38,39</sup>, Alessandra Chesi<sup>29</sup>, Jonas Bacelis<sup>40,41</sup>, Loic Yengo<sup>42,43</sup>, Maties Torrent<sup>44,45</sup>, Ioanna Ntalla<sup>46</sup>, Øyvind Helgeland<sup>47,48</sup>, Saskia Selzam<sup>49</sup>, Judith M. Vonk<sup>50</sup>, Mohammed H. Zafarmand<sup>51,52,53</sup>, Barbara Heude<sup>54</sup>, Ismaa Sadaf Farooqi<sup>55</sup>, Akram Alyass<sup>56</sup>, Robin N. Beaumont<sup>57</sup>, Christian T. Have<sup>8</sup>, Peter Rzehak<sup>58</sup>, Jose Ramon Bilbao<sup>59,60,61</sup>, Theresia M. Schnurr<sup>8</sup>, Inês Barroso<sup>62,63</sup>, Klaus Bønnelykke<sup>6</sup>, Lawrence J. Beilin<sup>64</sup>, Lisbeth Carstensen<sup>35</sup>, Marie-Aline Charles<sup>54</sup>, Bo Chawes<sup>6</sup>, Karine Clément<sup>65</sup>, Ricardo Closa-Monasterolo<sup>66</sup>, Adnan Custovic<sup>67</sup>, Johan G. Eriksson<sup>68,69</sup>, Joaquin Escribano<sup>66</sup>, Maria Groen-Blokhuis<sup>24</sup>, Veit Grote<sup>58</sup>, Dariusz Gruszfeld<sup>70</sup>, Hakon Hakonarson<sup>5,71</sup>, Torben Hansen<sup>8</sup>, Andrew T. Hattersley<sup>57,72</sup>, Mette Hollensted<sup>8,73</sup>, Jouke-Jan Hottenga<sup>24</sup>, Elina Hyppönen<sup>74,75,76</sup>, Stefan Johansson<sup>47,77</sup>, Raimo Joro<sup>11</sup>, Mika Kähönen<sup>78,79</sup>, Ville Karhunen<sup>80,81</sup>, Wieland Kiess<sup>82</sup>, Bridget A. Knight<sup>72</sup>, Berthold Koletzko<sup>58</sup>, Andreas Kühnapfel<sup>14,15</sup>, Kathrin Landgraf<sup>82,83</sup>, Jean-Paul Langhendries<sup>84</sup>, Terho Lehtimäki<sup>85,86</sup>, Jaakko T. Leinonen<sup>87</sup>, Aihuali Li<sup>56</sup>, Virpi Lindi<sup>88†</sup>, Estelle Lowry<sup>20,89</sup>, Mariona Bustamante<sup>25,28,90</sup>, Carolina Medina-Gomez<sup>1,3,91</sup>, Mads Melbye<sup>35,92</sup>, Kim F. Michaelsen<sup>93</sup>, Camilla S. Morgen<sup>94,95</sup>, Trevor A. Mori<sup>64</sup>, Tenna R. H. Nielsen<sup>96,97</sup>, Harri Niilikoski<sup>98,99</sup>, Albertine J. Oldehinkel<sup>100</sup>, Katja Pahkala<sup>38,39,101</sup>, Kalliope Panoutsopoulou<sup>102</sup>, Oluf Pedersen<sup>8</sup>, Craig E. Pennell<sup>32</sup>, Christine Power<sup>103</sup>, Sijmen A. Reijneveld<sup>104</sup>, Fernando Rivadeneira<sup>3,91</sup>, Angela Simpson<sup>10</sup>, Peter D. Sly<sup>105,106</sup>, Jakob Stokholm<sup>6</sup>, Kook K. Teo<sup>107</sup>, Elisabeth Thiering<sup>31,58</sup>, Nicholas J. Timpson<sup>17,18</sup>, André G. Uitterlinden<sup>3,91,108</sup>, Catharina E. M. van Beijsterveldt<sup>24</sup>, Barbera D. C. van Schaik<sup>109</sup>, Marc Vaudel<sup>47,77</sup>, Elvira Verduci<sup>110</sup>, Rebecca K. Vinding<sup>6</sup>, Mandy Vogel<sup>82,83</sup>, Eleftheria Zeggini<sup>102,111,112</sup>, Sylvain Sebert<sup>20,89,113</sup>, Mads V. Lind<sup>93</sup>, Christopher D. Brown<sup>114</sup>, Loreto Santa-Marina<sup>115,116,117</sup>, Eva Reischl<sup>118</sup>, Christine Frithioff-Bojsøe<sup>8,97,119</sup>, David Meyre<sup>56</sup>, Eleanor Wheeler<sup>62,63</sup>, Ken Ong<sup>120</sup>, Ellen A. Nohr<sup>121</sup>, Tanja G. M. Vrijkotte<sup>51</sup>, Gerard H. Koppelman<sup>122</sup>, Robert Plomin<sup>49</sup>, Pål R. Njølstad<sup>47,123,124</sup>, George D. Dedoussis<sup>125</sup>, Philippe Froguel<sup>42,113</sup>, Thorkild I. A. Sørensen<sup>8,17,94</sup>, Bo Jacobsson<sup>40,41,48</sup>, Rachel M. Freathy<sup>57,126</sup>, Babette S. Zemel<sup>71,127</sup>, Olli Raitakari<sup>38,39,128</sup>, Martine Vrijheid<sup>25,28,90</sup>, Bjarke Feenstra<sup>35</sup>, Leo-Pekka Lyytikäinen<sup>85,86,129</sup>, Harold Snieder<sup>16</sup>, Holger Kirsten<sup>14,15</sup>, Patrick G. Holt<sup>130</sup>, Joachim Heinrich<sup>31,131,132</sup>, Elisabeth Widén<sup>87</sup>, Jordi Sunyer<sup>25,28,90,133</sup>, Dorret I. Boomsma<sup>24,134</sup>, Marjo-Riitta Järvelin<sup>20,80,81</sup>, Antje Körner<sup>82,83</sup>, George Davey Smith<sup>17,18</sup>, Jens-Christian Holm<sup>8,97,119</sup>, Mustafa Atalay<sup>11</sup>, Clare Murray<sup>10</sup>, Hans Bisgaard<sup>6</sup>, Mark I. McCarthy<sup>21,22,135e</sup>, Early Growth Genetics Consortium, Vincent W. V. Jaddoe<sup>1,2‡</sup>, Struan F. A. Grant<sup>5,29,30,70,136‡</sup>, Janine F. Felix<sup>1,2‡*</sup>

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Vogelezang S, Bradfield JP, Ahluwalia TS, Curtin JA, Lakka TA, Grarup N, et al. (2020) Novel loci for childhood body mass index and shared heritability with adult cardiometabolic traits. PLoS Genet 16(10): e1008718. https://doi.org/10.1371/journal.pgen.1008718

**Editor:** Gregory P. Copenhaver, The University of North Carolina at Chapel Hill, UNITED STATES

**Received:** September 6, 2019

**Accepted:** March 16, 2020

**Published:** October 12, 2020

**Copyright:** © 2020 Vogelezang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** GWAS summary data are available at ftp://ftp.ebi.ac.uk/pub/databases/gwas/summary_statistics/GCST90002409 and has been deposited at the EGG Consortium website (https://egg-consortium.org/). Access to individual study-level data may be subject to local rules and regulations.

**Funding:** CP received funding from the National Institute for Health Research Biomedical Research Centre at Great Ormond Street Hospital for Children NHS Foundation Trust and University

1. The Generation R Study Group, Erasmus MC, University Medical Center, Rotterdam, the Netherlands,

2. Department of Pediatrics, Erasmus MC, University Medical Center, Rotterdam, the Netherlands,

3. Department of Epidemiology, Erasmus University Medical Center, Rotterdam, the Netherlands,

4. Quantinuum Research LLC, San Diego, California, United States of America, 5. Center for Applied Genomics, Division of Human Genetics, Children's Hospital of Philadelphia, Philadelphia, Pennsylvania, United States of America, 6. Copenhagen Prospective Studies on Asthma in Childhood, Herlev and Gentofte Hospital, University of Copenhagen, Copenhagen, Denmark, 7. Steno Diabetes Center Copenhagen, Gentofte, Denmark, 8. Novo Nordisk Foundation Center for Basic Metabolic Research, Faculty of Health and
College London. MMcC is a Wellcome Senior Investigator and an NIHR Senior Investigator. MMcC received funding from Wellcome (090532, 106130, 098381, 203141, 212259), NIDDK (U01-DK105535), and NIHR (NF-SI-0617-10090). TGMV was supported by ZonMW (TOP 40-00812-98-11010). DLC was supported by the American Diabetes Association Grant 1-17-PDF-077. SFAG is supported by the Daniel B. Burke Chair for Diabetes Research and NIH Grant R01 HD058886. NVT is funded by a pre-doctoral grant from the Agência de Gestió d'Ajuts Universitaris i de Recerca (2017 FI_B 00636), Generalitat de Catalunya - Fons Social Europeu. BK received personal funding from the European Research Council Advanced Grant META-GROWTH (ERC-2012-AdG - no. 322605). BF was supported by an Oak Foundation Fellowship. RMF and RNB are supported by Sir Henry Dale Fellowship (Wellcome Trust and Royal Society grant: WT104150). ATH is supported by a Wellcome Trust Senior Investigator award (grant number 098395/Z/12/Z). DM is supported by a Canada Research Chair. DLC was supported by the American Diabetes Association Grant 1-17-PDF-077. JTL was supported by the Finnish Cultural Foundation. DIB received a KNAW Academy Professor Award (PAH/6635). MH received PhD scholarship funding from TARGET (http://target.ku.dk), The Danish Diabetes Academy (http://danshdiabetesacademy.dk) and the Copenhagen Graduate School of Health and Medical Sciences. VWJ received funding from the Netherlands Organization for Health Research and Development (VIDI 016.136.361) and the European Research Council (ERC-2014-CoG-648916). ISF was supported by the European Research Council, Wellcome Trust (098497/Z/12/Z), Medical Research Council (MRC MC_UU_12012/5), the NIHR Cambridge Biomedical Research Centre, the Botnar Foundation, the Bernard Wolfe Health Neuroscience Endowment and the European Community's Seventh Framework Programme (FP7/2007-2013) project Beta-JUDO n°279153. IB acknowledges funding from Wellcome (WT206194). ES works in a unit that receives funding from the University of Bristol and the UK Medical Research Council (MC_UU_00011/1, MC_UU_00011/3). GDS works in the Medical Research Council Integrative Epidemiology Unit at the University of Bristol, which is supported by the Medical Research Council (MC_UU_00011/1). KC received funds from the French National Agency of Research, F-CRIN/FORCE. All funding for the studies that provided the underlying data for this manuscript can be found in S1 Text. The funders had no role in study design, data collection and

Medical Sciences, University of Copenhagen, Copenhagen, Denmark, **9** The Bioinformatics Center, Department of Biology, University of Copenhagen, Copenhagen, Denmark, **10** Division of Infection Immunity and Respiratory Medicine, School of Biological Sciences, The University of Manchester, Manchester Academic Health Science Centre, and Manchester University NHS Foundation Trust, Manchester, United Kingdom, **11** Institute of Biomedicine, Physiology, University of Eastern Finland, Kuopio, Finland, **12** Foundation for Research in Health Exercise and Nutrition, Kuopio Research Institute of Exercise Medicine, Kuopio, Finland, **13** Department of Clinical Physiology and Nuclear Medicine, Kuopio University Hospital, Kuopio, Finland, **14** Institute for Medical Informatics, Statistics and Epidemiology, University of Leipzig, Leipzig, Germany, **15** LIFE Research Center for Civilization Diseases, University of Leipzig, Leipzig, Germany, **16** Department of Epidemiology, University of Groningen, University Medical Center Groningen, Groningen, the Netherlands, **17** MRC Integrative Epidemiology Unit at the University of Bristol, Bristol, United Kingdom, **18** Population Health Sciences, Bristol Medical School, University of Bristol, Bristol, United Kingdom, **19** School of Oral and Dental Sciences, University of Bristol, Bristol, United Kingdom, **20** Center for Life Course Health Research, University of Oulu, Oulu, Finland, **21** Wellcome Centre for Human Genetics, University of Oxford, Oxford, United Kingdom, **22** Oxford Centre for Diabetes, Endocrinology and Metabolism, University of Oxford, Oxford, United Kingdom, **23** RIKEN Center for Integrative Medical Sciences, Yokohama, Kanagawa, Japan, **24** Department of Biological Psychology, Vrije Universiteit Amsterdam, Amsterdam, the Netherlands, **25** ISGlobal, Barcelona, Spain, **26** Centre for Genomic Regulation (CRG), The Barcelona Institute for Science and Technology, Barcelona, Spain, **27** BarcelonaBeta Brain Research Center (BBRC), Pasqual Maragall Foundation, Barcelona, Spain, **28** CIBER Epidemiología y Salud Pública (CIBERESP), Madrid, Spain, **29** Division of Human Genetics, Children's Hospital of Philadelphia, Philadelphia, Pennsylvania, United States of America, **30** Department of Genetics, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America, **31** Institute of Epidemiology, Helmholtz Zentrum München- German Research Center for Environmental Health, Neuherberg, Germany, **32** School of Medicine and Public Health, Faculty of Medicine and Health, The University of Newcastle, Newcastle, Australia, **33** Department of Medicine, University of Turku, Turku, Finland, **34** Division of Medicine, Turku University Hospital, Turku, Finland, **35** Department of Epidemiology Research, Statens Serum Institut, Copenhagen, Denmark, **36** Department of Statistics and Computational Research—Universitat de València, València, Spain, **37** Epidemiology and Environmental Health Joint Research Unit, FISABIO-Universitat Jaume I-Universitat de València, València, Spain, **38** Research Centre of Applied and Preventive Cardiovascular Medicine, University of Turku, Turku, Finland, **39** Centre for Population Health Research, University of Turku and Turku University Hospital, Turku, Finland, **40** Department of Obstetrics and Gynecology, Institute of Clinical Science, Sahlgrenska Academy, University of Gothenburg, Gothenburg Sweden, **41** Region Västra Götaland, Sahlgrenska University Hospital, Department of Obstetrics and Gynecology, Gothenburg Sweden, **42** University Lille, Centre National de la Recherche Scientifique, Institut Pasteur de Lille, UMR 8199—European Genomic Institute for Diabetes, Lille, France, **43** Institute for Molecular Bioscience, The University of Queensland, Brisbane, Australia, **44** Área de Salut de Menorca ib-salut, Menorca, Spain, **45** Institut d'Investigació Sanità Illes Balears (IdISBa), Palma de Mallorca, Spain, **46** William Harvey Research Institute, Barts and the London School of Medicine and Dentistry, Queen Mary University of London, London, United Kingdom, **47** KG Jebsen Center for Diabetes Research, Department of Clinical Science, University of Bergen, Bergen, Norway, **48** Department of Genetics and Bioinformatics, Health Data and Digitalization, Norwegian Institute of Public Health, Oslo, Norway, **49** Social, Genetic and Developmental Psychiatry Centre, Institute of Psychiatry, Psychology and Neuroscience, King's College London, London, United Kingdom, **50** Department of Epidemiology, GRIAC (Groningen Research Institute for Asthma and COPD), University of Groningen, University Medical Center Groningen, Groningen, the Netherlands, **51** Department of Public Health, Amsterdam Public Health Research Institute, Amsterdam UMC, University of Amsterdam, Amsterdam, the Netherlands, **52** Department of Obstetrics & Gynecology, Amsterdam UMC, University of Amsterdam, Amsterdam, the Netherlands, **53** Department of Clinical Epidemiology, Biostatistics and Bioinformatics, Amsterdam Public Health Research Institute, Amsterdam UMC, University of Amsterdam, Amsterdam, the Netherlands, **54** Université de Paris, CRESS, INSERM, INRA, Paris, France, **55** University of Cambridge Metabolic Research Laboratories and NIHR Cambridge Biomedical Research Centre, Wellcome Trust-MRC Institute of Metabolic Science, Addenbrooke's Hospital, Cambridge, United Kingdom, **56** Department of Health Research Methods, Evidence, and Impact, McMaster University, Hamilton, Canada, **57** Institute of Biomedical and Clinical Science, College of Medicine and Health, University of Exeter, Exeter, United Kingdom, **58** Division of Metabolic and Nutritional Medicine, Dr. von Hauner Children's Hospital, Ludwig-Maximilians Universität München (LMU), Munich, Germany, **59** University of the Basque Country (UPV/EHU), Leioa, Spain, **60** Biocrues-Bizkaia Health Research Institute, Barakaldo, Spain, **61** CIBER Diabetes y Enfermedades Metabólicas (CIBERDEM), Spain, **62** Wellcome Sanger Institute, Cambridge, United Kingdom, **63** MRC Epidemiology Unit, University of Cambridge, Cambridge, United Kingdom, **64** Medical School, The University of Western Australia, Perth, Western Australia, Australia, **65** Nutrition and Obesities; systemic approaches research unit, Sorbonne University, INSERM, Pitie-Salpêtrière Hospital, Assistance Publique hôpital de Paris, Paris, France, **66** Pediatrics, Nutrition and Development Research Unit, Universitat Rovira i Virgili, IISPV, Reus, Spain,
analysis, data interpretation, decision to publish, or
preparation of the manuscript.

**Competing interests:** I have read the journal's policy and the authors of this manuscript have the following competing interests: MMcC: The views expressed in this article are those of the author(s) and not necessarily those of the NHS, the NIHR, or the Department of Health. He serves on advisory panels for Pfizer, NovoNordisk, Zoe Global; has received honoraria from Merck, Pfizer, NovoNordisk and Eli Lilly; has stock options in Zoe Global; has received research funding from Abbvie, Astra Zeneca, Boehringer Ingelheim, Eli Lilly, Janssen, Merck, NovoNordisk, Pfizer, Roche, Sanofi Aventis, Servier & Takeda. MS receives funding from Pfizer Inc. for a project not related to this research. IB and spouse own stock in GlaxoSmithKline and Incyte Corp. ZE and CDB currently serve on the editorial board of PLOS Genetics. AC reports personal fees from Novartis, personal fees from Thermo Fisher Scientific, personal fees from Philips, personal fees from Sanofi, personal fees from Stallergenes Greer, outside the submitted work. KC in involved in consultancy for Danone Research, LNC-therapeutic and Confo-therapeutic but no personal funding is received and this activity not linked to the present research.

67 National Heart and Lung Institute, Imperial College London, London, United Kingdom, 68 Department of General Practice and Primary Health Care, University of Helsinki and Helsinki University Hospital, Helsinki, Finland, 69 Department of Obstetrics & Gynecology, Yong Loo Lin School of Medicine, National University of Singapore, Singapore, 70 Neonatal Department, Children's Memorial Health Institute, Warsaw, Poland, 71 Department of Pediatrics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, United States of America, 72 NIHR Exeter Clinical Research Facility, College of Medicine and Health, University of Exeter, and Royal Devon and Exeter NHS Foundation Trust, Exeter, United Kingdom, 73 The Danish Diabetes Academy, Odense, Denmark, 74 Australian Centre for Precision Health, University of South Australia Cancer Research Institute, Adelaide, Australia, 75 Population, Policy and Practice, UCL Institute of Child Health, University College London, London, United Kingdom, 76 South Australian Health and Medical Research Institute, Adelaide, Australia, 77 Department of Medical Genetics, Haukeland University Hospital, Bergen, Norway, 78 Department of Clinical Physiology, Tampere University Hospital, Tampere, Finland, 79 Department of Clinical Physiology, Finnish Cardiovascular Research Center—Tampere, Faculty of Medicine and Health Technology, Tampere University, Tampere, Finland, 80 Department of Epidemiology and Biostatistics, School of Public Health, Imperial College, London, United Kingdom, 81 MRC-PHE Centre for Environment and Health, School of Public Health, Imperial College, London, United Kingdom, 82 Center for Pediatric Research, University Hospital for Children and Adolescents, University of Leipzig, Leipzig, Germany, 83 Integrated Research and Treatment Center (IFB) Adiposity Diseases, University of Leipzig, Leipzig, Germany, 84 CHC-Health Group, Liège, Belgium, 85 Department of Clinical Chemistry, Fimlab Laboratories, Tampere, Finland, 86 Department of Clinical Chemistry, Finnish Cardiovascular Research Center—Tampere, Faculty of Medicine and Health Technology, Tampere University, Tampere, Finland, 87 Institute For Molecular Medicine Finland, FIMM, University of Helsinki, Helsinki, Finland, 88 University of Eastern Finland Library Kuopio, Kuopio, Finland, 89 Biocenter Oulu, Oulu University Hospital, Oulu, Finland, 90 Universitat Pompeu Fabra (UPF), Barcelona, Spain, 91 Department of Internal Medicine, Erasmus MC, University Medical Center, Rotterdam, the Netherlands, 92 Department of Medicine, Stanford School of Medicine, Stanford, California, United States of America, 93 Department of Nutrition, Exercise and Sports, Faculty of Science, University of Copenhagen, Copenhagen, Denmark, 94 Department of Public Health, Section of Epidemiology, University of Copenhagen, Copenhagen, Denmark, 95 National Institute of Public Health, University of Southern Denmark, Copenhagen, Denmark, 96 Department of Pediatrics, Hvidovre Hospital, Hvidovre, Denmark, 97 The Children's Obesity Clinic, Department of Pediatrics, Copenhagen University Hospital Holbæk, Holbæk, Denmark, 98 Department of Physiology, University of Turku, Turku, Finland, 99 Department of Pediatrics, University of Turku, Turku, Finland, 100 Interdisciplinary Center Psychopathology and Emotion Regulation, University of Groningen, University Medical Center, Groningen, the Netherlands, 101 Paavo Nurmi Centre, Sports and Exercise Medicine Unit, Department of Physical Activity and Health, University of Turku, Turku, Finland, 102 Wellcome Trust Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridgeshire, United Kingdom, 103 Population, Policy and Practice, UCL Great Ormond Street Institute of Child Health, University College London, London, United Kingdom, 104 Department of Health Sciences, University of Groningen, University Medical Center Groningen, Groningen, the Netherlands, 105 Child Health Research Centre, University of Queensland, Brisbane, Australia, 106 World Health Organization, WHO Collaborating Centre for Children's Health and Environment, Brisbane, Queensland, Australia, 107 Department of Medicine, McMaster University, Hamilton, Canada; Department of Health Research Methods, Evidence, and Impact, McMaster University, Hamilton, Canada, 108 Netherlands Genomics Initiative (NGI)-sponsored Netherlands Consortium for Healthy Aging NCHA), Leiden, the Netherlands, 109 Bioinformatics Laboratory, Department of Clinical Epidemiology, Biostatistics and Bioinformatics, Amsterdam Public Health Research Institute, Amsterdam UMC, University of Amsterdam, Amsterdam, the Netherlands, 110 Department of Pediatrics, San Paolo Hospital, University of Milan, Milan, Italy, 111 Institute of Translational Genomics, Helmholtz Zentrum München-German Research Center for Environmental Health, Neuherberg, Germany, 112 TUM School of Medicine, Technical University of Munich and Klinikum Rechts der Isar, Munich, Germany, 113 Section of Genomics of Common Disease, Department of Medicine, Imperial College London, London, United Kingdom, 114 Department of Genetics, University of Pennsylvania, Perelman School of Medicine, Philadelphia, Pennsylvania, United States of America, 115 Consortium for Biomedical Research in Epidemiology and Public Health (CIBER en Epidemiología y Salud Publica-CIBERESP), Barcelona, Spain, 116 Biodonostia Health Research Institute, San Sebastian, Spain, 117 Subdirección Salud Pública de Gipuzkoa, San Sebastian, Spain, 118 Research Unit of Molecular Epidemiology, Institute of Epidemiology, Helmholtz Zentrum Muenchen, Munich, Germany, 119 University of Copenhagen, Faculty of Health and Medical Sciences, Copenhagen N, Denmark, 120 Medical Research Council Epidemiology Unit & Department of Paediatrics, University of Cambridge, Addenbrooke's Hospital, Cambridge, England, 121 Research Unit for Gynaecology and Obstetrics, Institute of Clinical Research, University of Southern Denmark, Odense, Denmark, 122 University Medical Center Groningen, University of Groningen, Department of Pediatric Pulmonology and Pediatric Allergology, Beatrix Children's Hospital, GRIAC (Groningen Research Institute for Asthma and COPD), Groningen, the Netherlands, 123 Department of Pediatrics and Adolescents, Haukeland University Hospital, Bergen, Norway, 124 Broad Institute of Harvard and MIT, Cambridge, Massachusetts, United States of America, 125 Department of Nutrition and Dietetics,
School of Health Science and Education, Harokopio University, Athens, Greece, **126** Medical Research Council Integrative Epidemiology Unit, University of Bristol, Bristol, United Kingdom, **127** Division of Gastroenterology, Hepatology and Nutrition, The Children's Hospital of Philadelphia, Philadelphia, Pennsylvania, United States of America, **128** Department of Clinical Physiology and Nuclear Medicine, Turku University Hospital, Turku, Finland, **129** Department of Cardiology, Heart Center, Tampere University Hospital, Tampere, Finland, **130** Telethon Kids Institute, The University of Western Australia, Perth, Western Australia, Australia, **131** Institute and Outpatient Clinic for Occupational, Social and Environmental Medicine, Inner City Clinic, University Hospital Munich, Ludwig-Maximilians-Universität of Munich, Munich, Germany, **132** Allergy and Lung Health Unit, Melbourne School of Population and Global Health, The University of Melbourne, Melbourne, Australia, **133** Hospital del Mar Medical Research Institute (IMIM), Barcelona, Spain, **134** Amsterdam Public Health research institute and Amsterdam Reproduction & Development research Institute, Amsterdam, the Netherlands, **135** Oxford National Institute for Health Research (NIHR) Biomedical Research Centre, Churchill Hospital, Oxford, United Kingdom, **136** Center for Spatial and Functional Genomics, Children's Hospital of Philadelphia, Philadelphia, Pennsylvania, United States of America

© These authors contributed equally to this work.

† Deceased.

‡ These authors shared position on this work.

* j.felix@erasmusmc.nl

## Abstract

The genetic background of childhood body mass index (BMI), and the extent to which the well-known associations of childhood BMI with adult diseases are explained by shared genetic factors, are largely unknown. We performed a genome-wide association study meta-analysis of BMI in 61,111 children aged between 2 and 10 years. Twenty-five independent loci reached genome-wide significance in the combined discovery and replication analyses. Two of these, located near *NEDD4L* and *SLC45A3*, have not previously been reported in relation to either childhood or adult BMI. Positive genetic correlations of childhood BMI with birth weight and adult BMI, waist-to-hip ratio, diastolic blood pressure and type 2 diabetes were detected ($R_g$ ranging from 0.11 to 0.76, P-values <0.002). A negative genetic correlation of childhood BMI with age at menarche was observed. Our results suggest that the biological processes underlying childhood BMI largely, but not completely, overlap with those underlying adult BMI. The well-known observational associations of BMI in childhood with cardio-metabolic diseases in adulthood may reflect partial genetic overlap, but in light of previous evidence, it is also likely that they are explained through phenotypic continuity of BMI from childhood into adulthood.

## Author summary

Although twin studies have shown that body mass index (BMI) is highly heritable, many common genetic variants involved in the development of BMI have not yet been identified, especially in children. We studied associations of more than 40 million genetic variants with childhood BMI in 61,111 children aged between 2 and 10 years. We identified 25 genetic variants that were associated with childhood BMI. Two of these have not been implicated for BMI previously, located close to the genes *NEDD4L* and *SLC45A3*. We also show that the genetic background of childhood BMI overlaps with that of birth weight, adult BMI, waist-to-hip-ratio, diastolic blood pressure, type 2 diabetes, and age at menarche. Our results suggest that the biological processes underlying childhood BMI largely overlap with those underlying adult BMI. However, the overlap is not complete.
Additionally, the genetic backgrounds of childhood BMI and other cardio-metabolic phenotypes are overlapping. This may mean that the associations of childhood BMI and later cardio-metabolic outcomes are partially explained by shared genetics, but it could also be explained by the strong association of childhood BMI with adult BMI.

## Introduction

Childhood obesity is a major public health problem with impact on health in both the short and the long term [1]. Besides the well-established lifestyle and behavioral factors, genetics influence the risk of obesity, with reported heritability estimates from twin studies for body mass index (BMI) ranging from 40 to 70% [2,3]. An estimated 17 to 27% seems to be explained by common variants [4–6]. Large genome-wide association studies (GWAS) have identified 941 loci associated with adult BMI, accounting for 5% of the phenotypic variation [7]. Less is known about the genetic background of childhood BMI. A previous GWAS of BMI among 35,668 children identified 15 associated loci, accounting for 2% of the phenotypic variance [8]. Of these loci, 12 were also associated with adult BMI [9,10]. The remaining 3 identified genetic loci, specifically associated with childhood BMI, suggest possible age-specific differences between the two stages of life or could indicate stronger effects for these genetic loci in childhood BMI than in adult BMI [11–13]. Thus far, most common variants explaining the genetic variability of childhood BMI remain undetected. It is well known that obesity in early-life tends to track into later life [14]. Furthermore, childhood obesity has been associated with a lower age at menarche and with non-communicable diseases in later life, including hypertension, dyslipidemia, type 2 diabetes, neurodegenerative disease and asthma [15–19]. Findings from recent studies suggest a shared genetic background for BMI in childhood and adulthood [8,20,21]. To which extent the associations of childhood BMI with common adult diseases are genetically explained, has not been explored in detail.

We aimed to study the genetic background of childhood BMI by performing a two-stage GWAS meta-analysis consisting of 41 studies with a total sample size of 61,111 children of European ancestry. We also examined the genetic correlations of childhood BMI with anthropometric, cardio-metabolic, respiratory, neurocognitive and endocrinological traits in adults, using GWAS summary statistics from various consortia.

## Results

### Identification of genome-wide significant loci for childhood BMI

Sex- and age-adjusted Standard Deviation Scores (SDS) were created for BMI at the latest time point (oldest age, if multiple measurements were available) between 2 and 10 years using the same software and external reference across all studies (LMS growth; Pan H, Cole TJ, 2012; http://www.healthforallchildren.co.uk). Individual study characteristics are shown in S1 Table. The discovery meta-analysis included data from 26 studies ($N_{discovery} = 39,620$) with data imputed to the 1000 Genomes Project or The Haplotype Reference Consortium (HRC). We performed a fixed-effects inverse variance-weighted meta-analysis and performed conditional analyses based on summary-level statistics and Linkage Disequilibrium (LD) estimation between SNPs in Genome-wide Complex Trait Analysis (GCTA) to select independently associated SNPs at each locus on the basis of conditional P-values [22]. Seventeen independent SNPs reached genome-wide significance (P-values <5 × 10⁻⁸) and thirty SNPs showed suggestive association with childhood BMI (P-values >5 × 10⁻⁸ and <5 × 10⁻⁶). A Manhattan plot of
<figure><img src="image_4.png" /><figcaption>**Fig 1.** Manhattan plot of results of the discovery meta-analysis of 26 single study GWAS. On the x-axis the chromosomes are shown. On the y-axis the-log 10 of the P-value is shown. Novel SNPs are shown in green. Independent SNPs are shown in blue. Known SNPs are shown in black. The genome wide significance cutoff of $5 \times 10^{-8}$ is represented by the grey dotted line.</figcaption></figure>

https://doi.org/10.1371/journal.pgen.1008718.g001

the discovery meta-analysis is shown in Fig 1. No evidence of inflation due to population stratification or cryptic relatedness or other confounders was observed (genomic inflation factor ($\lambda$да) = 1.05; LD-score regression intercept = 1.0) (S1 Fig) [23]. All 47 independent SNPs identified in the discovery meta-analysis were taken forward for analysis in 15 replication cohorts ($N_{replication}$ = 21,491) and results of the two stages were then combined. Results of the discovery, replication and combined meta-analyses are shown in Table 1 and S2 Table and S3 Table. Results of the discovery analysis for SNPs with P-values $<5 \times 10^{-6}$ are shown in S4 Table. As the replication stage might lack power to replicate SNPs from the discovery analysis, we consider the joint analysis as the primary analysis.

In total, 25 loci achieved genome-wide significance in the combined meta-analysis. We defined a SNP as representing a known BMI-locus if it was within 500 kb of and in LD ($r^2 \ge 0.2$) with a previously reported BMI-associated signal. Of the 25 SNPs, two were novel and had not been previously associated with BMI in either adults or children: rs1094647 near *SLC45A3* and rs184566112 near *NEDD4L*. Per additional risk allele (G, allele frequency = 0.55) of rs1094647 (*SLC45A3*), childhood BMI increased by 0.04 SDS (Standard Error (SE) = 0.01, P-value = $7.20 \times 10^{-10}$), equal to 0.09 kg/m$^2$. Per additional risk allele (A, allele frequency = 0.84) of rs184566112 (*NEDD4L*), childhood BMI increased by 0.06 SDS (SE = 0.01; P-value = $4.24 \times 10^{-8}$), equal to 0.11 kg/m$^2$. Regional plots of the 2 novel SNPs are shown in Fig 2.

Despite the fact that these novel SNPs were not associated with either childhood or adult BMI previously, they have been reported to be associated with other anthropometric phenotypes. Rs1094647 (*SLC45A3*) has been associated with both height and whole-body fat-free mass in adulthood [24–26]. Additionally, rs708724, which is in high LD with rs1094647 ($r^2 = 0.70$) was associated with adult weight [24–26]. Rs184566112 (*NEDD4L*) is located in the same
<table><caption>Table 1. Results of the discovery, replication and combined analyses for the 47 loci with P-values <5 x 10⁻⁶ in the discovery phase.</caption><thead><tr><th>SNP</th><th>CHR</th><th>Position</th><th>Nearest gene</th><th>EA/ non_EA</th><th>P-value discovery</th><th>P-value replication</th><th>EAF<sup>a</sup></th><th>Beta combined</th><th>SE combined</th><th>P-value combined</th></tr></thead><tbody><tr><td>rs11676272<sup>b,c</sup></td><td>2</td><td>25141538</td><td>ADCY3</td><td>G/A</td><td>2.37 x 10<sup>-21</sup></td><td>1.88 x 10<sup>-10</sup></td><td>0.46</td><td>0.071</td><td>0.006</td><td><strong>3.79 x 10<sup>-30</sup></strong></td></tr><tr><td>rs7138803<sup>b,c</sup></td><td>12</td><td>50247468</td><td>BCDIN3D</td><td>A/G</td><td>7.12 x 10<sup>-20</sup></td><td>7.39 x 10<sup>-12</sup></td><td>0.37</td><td>0.072</td><td>0.006</td><td><strong>4.23 x 10<sup>-30</sup></strong></td></tr><tr><td>rs939584<sup>b,c,d</sup></td><td>2</td><td>621558</td><td>TMEM18</td><td>T/C</td><td>8.85 x 10<sup>-26</sup></td><td>2.52 x 10<sup>-6</sup></td><td>0.83</td><td>0.092</td><td>0.008</td><td><strong>3.73 x 10<sup>-29</sup></strong></td></tr><tr><td>rs17817449<sup>b,c</sup></td><td>16</td><td>53813367</td><td>FTO</td><td>G/T</td><td>1.69 x 10<sup>-17</sup></td><td>3.21 x 10<sup>-11</sup></td><td>0.40</td><td>0.069</td><td>0.006</td><td><strong>2.98 x 10<sup>-27</sup></strong></td></tr><tr><td>rs12042908<sup>b,c</sup></td><td>1</td><td>74997762</td><td>FPGT-TNNI3K, TNNI3K</td><td>A/G</td><td>2.77 x 10<sup>-14</sup></td><td>2.01 x 10<sup>-12</sup></td><td>0.46</td><td>0.064</td><td>0.006</td><td><strong>6.37 x 10<sup>-25</sup></strong></td></tr><tr><td>rs543874<sup>b,c</sup></td><td>1</td><td>177889480</td><td>SEC16B</td><td>G/A</td><td>1.61 x 10<sup>-15</sup></td><td>6.96 x 10<sup>-8</sup></td><td>0.19</td><td>0.075</td><td>0.008</td><td><strong>6.02 x 10<sup>-22</sup></strong></td></tr><tr><td>rs56133711<sup>b</sup></td><td>11</td><td>27723334</td><td>BDNF</td><td>A/G</td><td>2.00 x 10<sup>-10</sup></td><td>3.69 x 10<sup>-6</sup></td><td>0.24</td><td>0.056</td><td>0.007</td><td><strong>3.75 x 10<sup>-15</sup></strong></td></tr><tr><td>rs2076308<sup>b,c</sup></td><td>6</td><td>50791640</td><td>TFAP2B</td><td>C/G</td><td>1.59 x 10<sup>-12</sup></td><td>0.01</td><td>0.19</td><td>0.058</td><td>0.008</td><td><strong>3.07 x 10<sup>-13</sup></strong></td></tr><tr><td>rs4477562<sup>b,c,e</sup></td><td>13</td><td>54104968</td><td>LINC00558</td><td>T/C</td><td>8.29 x 10<sup>-13</sup></td><td>0.01</td><td>0.13</td><td>0.065</td><td>0.009</td><td><strong>5.81 x 10<sup>-13</sup></strong></td></tr><tr><td>rs571312<sup>b,c</sup></td><td>18</td><td>57839769</td><td>MC4R</td><td>A/C</td><td>2.00 x 10<sup>-10</sup></td><td>4.84 x 10<sup>-4</sup></td><td>0.23</td><td>0.052</td><td>0.007</td><td><strong>8.80 x 10<sup>-13</sup></strong></td></tr><tr><td>rs12641981<sup>b,c</sup></td><td>4</td><td>45179883</td><td>GNPDA2</td><td>T/C</td><td>4.19 x 10<sup>-8</sup></td><td>7.12 x 10<sup>-6</sup></td><td>0.44</td><td>0.045</td><td>0.006</td><td><strong>1.29 x 10<sup>-12</sup></strong></td></tr><tr><td>rs62107261<sup>f</sup></td><td>2</td><td>422144</td><td>FAM150B</td><td>T/C</td><td>3.17 x 10<sup>-7</sup></td><td>6.64 x 10<sup>-6</sup></td><td>0.95</td><td>0.121</td><td>0.018</td><td><strong>9.93 x 10<sup>-12</sup></strong></td></tr><tr><td>rs114285994<sup>b</sup></td><td>16</td><td>19935763</td><td>GPRC5B</td><td>G/A</td><td>1.41 x 10<sup>-10</sup></td><td>5.98 x 10<sup>-3</sup></td><td>0.87</td><td>0.063</td><td>0.009</td><td><strong>1.11 x 10<sup>-11</sup></strong></td></tr><tr><td>rs144376234<sup>b,c</sup></td><td>1</td><td>110114504</td><td>GNAI3</td><td>T/C</td><td>1.35 x 10<sup>-8</sup></td><td>2.21 x 10<sup>-3</sup></td><td>0.04</td><td>0.111</td><td>0.017</td><td><strong>1.38 x 10<sup>-10</sup></strong></td></tr><tr><td>rs1094647</td><td>1</td><td>205655378</td><td>SLC45A3</td><td>G/A</td><td>2.46 x 10<sup>-6</sup></td><td>7.45 x 10<sup>-5</sup></td><td>0.55</td><td>0.038</td><td>0.006</td><td><strong>7.20 x 10<sup>-10</sup></strong></td></tr><tr><td>rs76227980<sup>f</sup></td><td>18</td><td>58036384</td><td>MC4R</td><td>C/T</td><td>1.71 x 10<sup>-6</sup></td><td>1.19 x 10<sup>-4</sup></td><td>0.98</td><td>0.140</td><td>0.023</td><td><strong>8.68 x 10<sup>-10</sup></strong></td></tr><tr><td>rs13107325<sup>b</sup></td><td>4</td><td>103188709</td><td>SLC39A8</td><td>T/C</td><td>3.51 x 10<sup>-8</sup></td><td>4.84 x 10<sup>-3</sup></td><td>0.07</td><td>0.082</td><td>0.014</td><td><strong>1.38 x 10<sup>-9</sup></strong></td></tr><tr><td>rs62500888<sup>c</sup></td><td>8</td><td>28061823</td><td>ELP3</td><td>A/G</td><td>6.91 x 10<sup>-10</sup></td><td>0.09</td><td>0.57</td><td>0.037</td><td>0.006</td><td><strong>1.81 x 10<sup>-9</sup></strong></td></tr><tr><td>rs114670539<sup>c</sup></td><td>2</td><td>207064335</td><td>GPR1</td><td>T/C</td><td>3.16 x 10<sup>-8</sup></td><td>0.01</td><td>0.05</td><td>0.088</td><td>0.015</td><td><strong>1.92 x 10<sup>-9</sup></strong></td></tr><tr><td>rs61765651<sup>b</sup></td><td>1</td><td>72754314</td><td>NEGR1</td><td>C/T</td><td>9.50 x 10<sup>-9</sup></td><td>0.03</td><td>0.83</td><td>0.047</td><td>0.008</td><td><strong>4.99 x 10<sup>-9</sup></strong></td></tr><tr><td>rs7719067<sup>b</sup></td><td>5</td><td>153538241</td><td>GALNT10</td><td>A/G</td><td>3.96 x 10<sup>-7</sup></td><td>4.36 x 10<sup>-3</sup></td><td>0.43</td><td>0.036</td><td>0.006</td><td><strong>6.54 x 10<sup>-9</sup></strong></td></tr><tr><td>rs11030391<sup>f</sup></td><td>11</td><td>28644626</td><td>METTL15</td><td>A/G</td><td>4.73 x 10<sup>-7</sup></td><td>6.49 x 10<sup>-3</sup></td><td>0.63</td><td>0.036</td><td>0.006</td><td><strong>1.51 x 10<sup>-8</sup></strong></td></tr><tr><td>rs184566112</td><td>18</td><td>55943926</td><td>NEDD4L</td><td>A/T</td><td>4.40 x 10<sup>-6</sup></td><td>1.24 x 10<sup>-3</sup></td><td>0.84</td><td>0.057</td><td>0.011</td><td><strong>4.24 x 10<sup>-8</sup></strong></td></tr><tr><td>rs116664060</td><td>6</td><td>31592524</td><td>PRRC2A</td><td>C/G</td><td>3.03 x 10<sup>-6</sup></td><td>3.28 x 10<sup>-3</sup></td><td>0.18</td><td>0.049</td><td>0.009</td><td><strong>4.63 x 10<sup>-8</sup></strong></td></tr><tr><td>rs11215427<sup>c</sup></td><td>11</td><td>115093438</td><td>CADM1</td><td>G/C</td><td>1.25 x 10<sup>-7</sup></td><td>0.05</td><td>0.74</td><td>0.039</td><td>0.007</td><td><strong>4.64 x 10<sup>-8</sup></strong></td></tr><tr><td>rs1336980<sup>c</sup></td><td>9</td><td>129377855</td><td>LMX1B</td><td>C/G</td><td>6.61 x 10<sup>-7</sup></td><td>0.03</td><td>0.36</td><td>0.033</td><td>0.006</td><td>1.17 x 10<sup>-7</sup></td></tr><tr><td>rs146823532<sup>f</sup></td><td>1</td><td>74979126</td><td>FPGT-TNNI3K, TNNI3K</td><td>A/G</td><td>4.09 x 10<sup>-7</sup></td><td>0.03</td><td>0.97</td><td>0.114</td><td>0.022</td><td>1.33 x 10<sup>-7</sup></td></tr><tr><td>rs79386556</td><td>13</td><td>71229046</td><td>LINC00348</td><td>A/G</td><td>1.84 x 10<sup>-6</sup></td><td>0.04</td><td>0.04</td><td>0.095</td><td>0.019</td><td>4.83 x 10<sup>-7</sup></td></tr><tr><td>rs9942489</td><td>6</td><td>35323709</td><td>PPARD</td><td>A/T</td><td>3.62 x 10<sup>-6</sup></td><td>0.02</td><td>0.04</td><td>0.081</td><td>0.016</td><td>5.56 x 10<sup>-7</sup></td></tr><tr><td>rs17086809</td><td>9</td><td>86708695</td><td>RMI1</td><td>C/T</td><td>3.06 x 10<sup>-6</sup></td><td>0.05</td><td>0.33</td><td>0.038</td><td>0.007</td><td>9.30 x 10<sup>-7</sup></td></tr><tr><td>rs80332495</td><td>5</td><td>19191677</td><td>CDH18</td><td>A/G</td><td>3.78 x 10<sup>-7</sup></td><td>0.26</td><td>0.94</td><td>0.090</td><td>0.019</td><td>1.14 x 10<sup>-6</sup></td></tr><tr><td>rs4594227</td><td>15</td><td>84497207</td><td>ADAMTSL3</td><td>A/G</td><td>4.46 x 10<sup>-6</sup></td><td>0.07</td><td>0.56</td><td>0.030</td><td>0.006</td><td>1.75 x 10<sup>-6</sup></td></tr><tr><td>rs2457463</td><td>10</td><td>70315687</td><td>TET1</td><td>G/T</td><td>1.79x 10<sup>-7</sup></td><td>0.43</td><td>0.04</td><td>0.159</td><td>0.034</td><td>2.69 x 10<sup>-6</sup></td></tr><tr><td>rs11865086<sup>b</sup></td><td>16</td><td>30130493</td><td>MAPK3</td><td>C/A</td><td>2.23 x 10<sup>-7</sup></td><td>0.47</td><td>0.53</td><td>0.029</td><td>0.006</td><td>3.31 x 10<sup>-6</sup></td></tr><tr><td>rs1565356</td><td>6</td><td>34046065</td><td>GRM4</td><td>C/A</td><td>1.76 x 10<sup>-6</sup></td><td>0.25</td><td>0.92</td><td>0.057</td><td>0.012</td><td>4.01 x 10<sup>-6</sup></td></tr><tr><td>rs2952863<sup>b</sup></td><td>4</td><td>130759647</td><td>C4orf33</td><td>T/G</td><td>3.10 x 10<sup>-6</sup></td><td>0.15</td><td>0.30</td><td>0.031</td><td>0.007</td><td>4.27 x 10<sup>-6</sup></td></tr><tr><td>rs2358954</td><td>12</td><td>66379504</td><td>HMGA2</td><td>T/G</td><td>3.07 x 10<sup>-6</sup></td><td>0.17</td><td>0.68</td><td>0.031</td><td>0.007</td><td>4.77 x 10<sup>-6</sup></td></tr><tr><td>rs7652876</td><td>3</td><td>179831733</td><td>PEX5L</td><td>A/C</td><td>3.22 x 10<sup>-6</sup></td><td>0.19</td><td>0.25</td><td>0.033</td><td>0.007</td><td>6.93 x 10<sup>-6</sup></td></tr><tr><td>rs4923207</td><td>11</td><td>24757325</td><td>LUZP2</td><td>T/G</td><td>1.07 x 10<sup>-6</sup></td><td>0.86</td><td>0.81</td><td>0.040</td><td>0.009</td><td>7.05 x 10<sup>-6</sup></td></tr><tr><td>rs7757288</td><td>6</td><td>55205502</td><td>GFRAL</td><td>G/A</td><td>6.64 x 10<sup>-7</sup></td><td>0.47</td><td>0.39</td><td>0.028</td><td>0.006</td><td>1.03 x 10<sup>-5</sup></td></tr><tr><td>rs6876477</td><td>5</td><td>50878621</td><td>ISL1</td><td>A/T</td><td>3.03 x 10<sup>-6</sup></td><td>0.35</td><td>0.76</td><td>0.032</td><td>0.007</td><td>1.31 x 10<sup>-5</sup></td></tr><tr><td>rs28599560</td><td>5</td><td>91791853</td><td>FLJ42709</td><td>A/G</td><td>3.99 x 10<sup>-7</sup></td><td>0.85</td><td>0.61</td><td>0.027</td><td>0.006</td><td>3.48 x 10<sup>-5</sup></td></tr><tr><td>rs117281273</td><td>8</td><td>42981400</td><td>SGK196</td><td>C/G</td><td>2.17 x 10<sup>-7</sup></td><td>0.94</td><td>0.97</td><td>0.081</td><td>0.020</td><td>3.48 x 10<sup>-5</sup></td></tr><tr><td>rs9695734</td><td>9</td><td>96407983</td><td>PHF2</td><td>C/T</td><td>1.12 x 10<sup>-6</sup></td><td>0.81</td><td>0.84</td><td>0.035</td><td>0.009</td><td>4.99 x 10<sup>-5</sup></td></tr><tr><td>rs72833479</td><td>17</td><td>45960449</td><td>SP2</td><td>A/G</td><td>7.96 x 10<sup>-7</sup></td><td>0.99</td><td>0.25</td><td>0.029</td><td>0.007</td><td>6.54 x 10<sup>-5</sup></td></tr></tbody></table>

(Continued)
<table><caption>**Table 1.** (Continued)</caption><thead><tr><th>SNP</th><th>CHR</th><th>Position</th><th>Nearest gene</th><th>EA/ non_EA</th><th>P-value discovery</th><th>P-value replication</th><th>EAF<sup>a</sup></th><th>Beta combined</th><th>SE combined</th><th>P-value combined</th></tr></thead><tbody><tr><td>rs142367753</td><td>2</td><td>128938956</td><td>UGGT1</td><td>C/G</td><td>4.03 x 10<sup>-6</sup></td><td>0.24</td><td>0.98</td><td>0.088</td><td>0.027</td><td>1.34 x 10<sup>-3</sup></td></tr><tr><td>rs6896578</td><td>5</td><td>76423090</td><td>ZBED3-AS1</td><td>C/T</td><td>3.52 x 10<sup>-6</sup></td><td>0.17</td><td>0.84</td><td>0.025</td><td>0.009</td><td>4.1 x 10<sup>-3</sup></td></tr></tbody></table>

CHR, chromosome; EA, effect allele; EAF, effect allele frequency; SE, standard error.

Bolded P-values indicate genome-wide significance in the combined analysis.

Detailed information on beta and SE of the discovery and replication stage separately can be found in S2 Table

a From combined analysis

b Locus previously reported for adult BMI

c Locus previously reported for childhood BMI

d Locus previously reported for adult body fat

e Locus previously reported for childhood obesity

f Independent SNP at the same locus selected by conditional analysis

https://doi.org/10.1371/journal.pgen.1008718.t001

region as rs6567160 (distance = 448 kb, r² <0.2), previously associated with adult body fat [27]. In the current study, we did not observe evidence for association between rs184566112 (NEDD4L, effect allele = G, allele frequency = 0.84) and body fat percentage measured by Dual energy X-ray Absorptiometry (age range 24 to 120 months) in 2,698 children from 4 cohorts (0.03 SDS (SE = 0.04, P-value = 0.51)). Individual study characteristics of studies with data on body fat percentage are shown in S5 Table. No evidence of association with childhood obesity was found for the two novel SNPs (P-values >0.11) [28].

We additionally identified 2 independent SNPs (METTL15 and PRRC2A) within 500 kb of previously reported SNPs associated with adult BMI, but only in weak LD with prior reported signals (r² <0.2). Similarly, we found 2 independent SNPs in regions that are known for both childhood and adult BMI (FAM150B and MC4R) [7,8,10]. Regional plots of the 4 independent SNPs at known loci are shown in S2 Fig. Of the remaining 19 SNPs, 6 mapped to loci previously associated with adult BMI (BDNF, GPRC5B, SLC39A8, NEGR1, GALNT10, and CADM1), 2 mapped to loci previously associated with childhood BMI only (ELP3 and GPR1) and 11 SNPs mapped to loci known to be associated with both adult and childhood BMI (ADCY3, BCDIN3D, TMEM18, FTO, FPGT-TNNI3K/TNNI3K, SEC16B, TFAP2B, LINC00558, MC4R, GNPDA2, and GNAI3) [7,8,10].

Overall, there was low heterogeneity between studies for the 25 SNPs, except for FTO (S2 Table) [29]. The broad age range included in the discovery meta-analysis of this study may conceal age-specific effects. Therefore, we performed a sensitivity analysis excluding studies of children aged <6 years (remaining Nsensitivity analysis = 55,354), which showed similar results (S6 Table) [30]. Additionally, we ran a sensitivity analysis excluding case-control studies and one excluding studies with a sample size <n = 500, showing similar results (S7 Table).

## Functional characterization

We used several strategies to gain insight into the functional characterization of the 25 SNPs leading the association signals with childhood BMI. A summary of relevant information from all strategies can be found in S8 Table.

First, we examined gene expression profiles of the nearest genes to the 25 SNPs from the combined meta-analysis with GTEx v7 in 53 tissues, using the tool FUMA [3,31]. We found differential expression of the 25 nearest genes in brain and salivary gland. In a second analysis of gene expression profiles in GTEx, we considered all genes in a region of 500 kb to either
A

<figure><img src="image_5.png" /></figure>

B

<figure><img src="image_6.png" /></figure>

**Fig 2.** Locus zoom plots of the 2 novel SNPs Regional association plot of the 2 novel SNPs in the 26 childhood BMI discovery studies. SNPs are plotted with their P-values (as-log10; left y-axis) as a function of genomic position (x-axis). Estimated recombination rates (right y-axis) taken from 1000 Genomes, March 2012 release are plotted to
reflect the local LD-structure around the top associated SNP (indicated with purple color) and the correlated proxies (indicated in colors). A. rs1094647 B. rs184566112.

https://doi.org/10.1371/journal.pgen.1008718.g002

side of the 25 SNPs. Using this strategy, we additionally found differential expression in liver, heart, kidney, pancreas, muscle, skin and adipose tissue [31].

Second, we assessed whether the 25 SNPs were associated with gene expression in whole adipose tissue, isolated adipocytes, and isolated stroma-vascular cells from the Leipzig Adipose Tissue Childhood Cohort [32]. Full results can be found in S9 Table. We observed differential gene expression associated with multiple SNPs. Rs1094647 (nearest gene: *SLC45A3*) was associated with gene expression of *PM20D1* (*P*<sub>FDR</sub> <0.05) in whole adipose tissue. We additionally found associations of rs114285994 (nearest gene: *GPRC5B*) with expression of *C16orf88* in isolated adipocytes. Rs115181845, which is in moderate LD (*r*<sup>2</sup> = 0.47) with rs144376234 (nearest gene: *GNAI3*), was associated with expression of *GSTM1* and *GSTM2* in whole adipose tissue, isolated adipocytes and isolated stroma-vascular cells (S8 Table and S9 Table). No associations with gene expression were observed for any of the other 22 SNPs.

Third, we used Bayesian colocalization analysis to examine evidence for colocalization between GWAs and eQTL signals and to identify additional candidate genes for the 25 SNPs (GTEx v7). Briefly, GWAS summary statistics were extracted for each eQTL for all SNPs that were present in the meta-analysis and that were in common to both GWAS and eQTL studies. In most pairs, no evidence for association was found with either trait. To define colocalization we used restriction to pairs of childhood BMI and eQTL signals with a high posterior probability for colocalization (See Methods and Materials for details) [33]. We found significant colocalizations at 6 loci (*ADCY3*, *DNAJC27-AS1*, *CENPO*, *ADAM23*, *LIN7C*, *TFAP2B*) across a range of tissues (S8 Table and S10A and S10B Table) [8].

Fourth, to explore biological processes, we used DAVID, with the 25 nearest genes as input, using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database [34,35]. Pathway analysis revealed one enriched biological process, cAMP signaling (P-value = 0.03).

Fifth, we performed a look-up in mouse-knockout data of the 25 nearest genes and, additionally, the genes that were indicated by colocalization and gene expression analysis. Mice in which *NEDD4L* was knocked out displayed neuronal abnormalities [36]. No related phenotypes were shown for *SLC45A3* or any of the 4 independent loci (*METTL1*, *PRRC2A*, *FAM150B*, and *MC4R*). Of the 19 known loci, *ADCY3* showed an association with increased total body fat in female heterozygous knockout mice, whereas *NEGR1* was associated with decreased lean body mass in male and female homozygous knockout mice (S8 Table). Full results can be found in S8 Table.

Sixth, among the 25 top SNPs, combined annotation-dependent depletion (CADD) scores >12.37, indicating potential pathogenicity of a SNP, were observed for rs13107325 (*SLC39A8*), rs56133711 (*BDNF*) and rs17817449 (*FTO*) (CADD scores of 34, 15.3 and 15.3, respectively) (S8 Table) [3,37].

## Genetic correlations of childhood BMI with adult phenotypes

First, to estimate the SNP heritability and the genetic correlations between childhood BMI and other traits from external GWAS meta-analysis data, we used LD-score regression [20]. SNP heritability was 0.23. There were positive genetic correlations of childhood BMI with several anthropometric and cardio-metabolic traits, including adult BMI (*R*<sub>g</sub> = 0.76, P-value = 1.45 × 10<sup>-112</sup>), waist-to-hip ratio (*R*<sub>g</sub> = 0.39, P-value = 1.57 × 10<sup>-20</sup>), body fat percentage (*R*<sub>g</sub> = 0.46, P-value = 7.99 × 10<sup>-44</sup>), diastolic blood pressure (*R*<sub>g</sub> = 0.11, P-value = 0.002), type 2 diabetes (*R*<sub>g</sub> = 0.19, P-value = 0.002), and coronary artery disease (*R*<sub>g</sub> = 0.14, P-value = 0.002).
value = 0.001) (Fig 3 and S11 Table). For birth weight, there were positive genetic correlations with childhood BMI, both when using fetal genetic effects on birth weight and when using maternal genetic effects on birth weight ($R_g$ = 0.20, P-value = 3.19 x 10⁻⁵ and $R_g$ = 0.12, P-value = 0.002 for fetal and maternal effects, respectively). Negative genetic correlations were observed between childhood BMI and total cholesterol ($R_g$ = -0.15, P-value = 0.001), high-density lipoprotein (HDL) ($R_g$ = -0.22, P-value = 8.65 x 10⁻⁶), and age at menarche ($R_g$ = -0.42, P-value = 1.03 x 10⁻³²). We did not find genetic correlations with any of the respiratory and neurocognitive phenotypes. Genetic correlations of childhood BMI with a selection of phenotypes that show evidence of association in observational studies are shown in Fig 3. Full results can be found in S11 Table.

Second, we did a look-up of the 25 SNPs in the adult BMI GWAS [7]. In total, 12 SNPs and 8 proxy SNPs ($r^2 \ge 0.87$) were available in the adult BMI study comprising ~700,000 individuals. No information was available on five loci, *FAM150B*, *GPR1*, *NEGR1*, *NEDD4L*, and *PRRC2A*. The directions of effect of all 20 SNPs were the same in adults as in children. Of these, 18 were genome-wide significantly associated with adult BMI (P-value < 5 x 10⁻⁸) and the other 2 SNPs, *SLC45A3* and *METTL15*, showed suggestive evidence of association (P-values < 2.1 x 10⁻⁶) (S12 Table). Effect sizes of these 20 SNPs for adult BMI were highly correlated with those for childhood BMI ($r^2$ = 0.86).

Third, we calculated a combined childhood BMI genetic risk score (GRS) of the 25 genome-wide significant SNPs, summing the number of BMI-increasing alleles weighted by their effect sizes from the combined meta-analysis. The GRS was associated with childhood BMI (P-value = 2.84 x 10⁻¹¹) in 1,169 children from the Tracking Adolescents' Individual Lives Survey (TRAILS) Cohort, aged 7 years, one of the largest replication cohorts (Fig 4). For each additional average risk allele in the GRS, childhood BMI increased by 0.06 SDS (SE = 0.009). This GRS explained 3.6% of the variance in childhood BMI. When calculating the risk score for the TRAILS cohort, effect estimates from the combined meta-analysis were used after excluding TRAILS from the meta-analysis. We additionally tested the GRS for

<figure><img src="image_7.png" /><figcaption>**Fig 3.** Genome-wide genetic correlations between childhood BMI and adult traits and diseases. On the x-axis the traits and diseases are shown. On the y-axis the genetic correlations ($R_g$) and corresponding standard errors, indicated by error bars, between childhood BMI and each trait were shown, estimated by LD score regression. The genetic correlation estimates ($R_g$) are colored according to their intensity and direction. Red indicates positive correlation, blue indicates negative correlation. References can be found in S11 Table.</figcaption></figure>

https://doi.org/10.1371/journal.pgen.1008718.g003
<figure><img src="image_8.png" /><figcaption>**Fig 4.** Associations of the weighted risk score with childhood BMI. Along the x-axis, categories of the risk score are shown together with the mean SDS BMI on the y-axis on the right and a line representing the regression of the mean SDS childhood BMI values for each category of the risk score. Along the y-axis on the left a histogram represents the number of individuals in each risk-score category. P-value is based on the continuous risk score. Analysis was performed in the TRAILS cohort (N = 1,169).</figcaption></figure>

https://doi.org/10.1371/journal.pgen.1008718.g004

association with adult BMI in the three sub-cohorts of the Rotterdam Study [38] (RS-I-1; n = 5,957, RS-II-1; n = 2,147 and RS-III-1; n = 2,998). We found the GRS to be associated with adult BMI in all study samples (P-values = 5.09 × 10⁻⁹, 0.02, and 1.49 × 10⁻¹⁰, respectively). Per additional average risk allele, adult BMI increased by 0.03 SDS (SE = 0.005), 0.02 SDS (SE = 0.009) and 0.04 SDS (SE = 0.007), explaining 0.6%, 0.2%, and 1.3% of the variance in adult BMI, respectively. No association was found of the GRS with birth weight and cardio-metabolic phenotypes, including insulin, triglycerides, low-density lipoprotein, HDL, total cholesterol, diastolic blood pressure and systolic blood pressure in 2,831 children aged 6 years from the Generation R Study if considering a Bonferroni corrected P-value of 0.00625 (S13 Table).

Discussion

In this large GWAS meta-analysis of childhood BMI among >60,000 children aged 2–10 years, we identified 25 genome-wide significant loci. Two of these loci, rs1094647 near *SLC45A3* and rs184566112 near *NEDD4L* had not been associated with BMI before. We observed moderate to strong genetic correlations of childhood BMI with several anthropometric, cardio-metabolic, and endocrinological traits in adulthood, suggesting a shared genetic background.
The closest genes to the two novel loci, *SLC45A3* and *NEDDL4*, have not been strongly linked to obesity in previous studies and databases, indicating that functional studies are needed to identify possible biological pathways. *SLC45A3*, encoding the solute carrier family 45, member 3 protein, also known as prostate cancer-associated protein 6, has been related to prostate-specific antigen serum concentrations and prostate cancer [39–42]. *NEDD4L*, ubiquitin protein ligase Nedd4-like, known for its role in the regulation of ion channel internalization and turnover, is suggested to play a role in the regulation of respiratory, cardiovascular, renal, and neuronal functions [36,43–45]. The independent SNPs identified at loci known from previous studies on adult or childhood BMI may represent fully independent signals, although due to the low LD, these SNPs might still tag the same causal variant as the previously identified SNPs.

Since there is no strong previous evidence supporting the closest genes to the 25 SNPs as the causal genes, we took multiple approaches for further functional characterization. As many different tissues have been implicated to play a role in body composition we chose to include all available tissues in the gene expression analysis. Using GTEx, we found differential expression of the 25 nearest genes in brain. This may be of interest as appetite regulation might play a role in the development of obesity [46–48]. Gene expression data revealed an association between one of the novel SNPs, rs1094647 (nearest gene: *SLC45A3*), and expression of *PM20D1* in whole adipose tissue. *PM20D1*, Peptidase M20 domain-containing 1, previously identified as a factor secreted by thermogenic adipose cells, is known for its association with insulin resistance, glucose intolerance and enhanced defense of body temperature in cold when knocked out in mice. Furthermore, increased circulating PM20D1, together with adiponectin, a well-known adipocyte-derived factor, leads to a higher energy expenditure and reduced adiposity in mice [49,50]. We used colocalization analysis to further identify candidate causal genes. This did not identify specific potential causal genes for rs1094647 (*SLC45A3*) and rs184566112 (*NEDD4L*). However, we identified ADCY3, DNAJC27-AS1, CENPO, ADAM23, LIN7C, TFAP2B as candidate genes for known loci across different tissues, including tibial nerve tissue, tibial artery tissue and the skin. No candidate genes were detected in biologically more relevant tissues, including subcutaneous or visceral adipose tissue.

Information on rs184566112 near *NEDD4L* was available in 24 out of 26 discovery cohorts that primarily used 1000 Genomes phase 1 imputed data (N = 37,104), thus clearly surviving our pre-set filter of having information in at least 50% of the number of studies and at least 50% of the total sample size in the discovery analysis. However, it was available in only less than half of the replication studies, mainly using 1000 Genomes phase 3 or HRC imputed data (N = 5,518) as this SNP was not included in these more recent reference panels. No other SNPs in high LD were available as proxy for this SNP in the replication analysis. Therefore, this signal needs to be interpreted with caution. However, no heterogeneity of this SNP between the discovery stage studies ($I^2 = 0$; P-value for heterogeneity = 0.98), a high imputation quality (weighted mean $R^2 = 0.89$) and the known association of another locus in the same region with adult body fat percentage might lends credibility to this signal, although further work is needed to unravel the details [27]. Previous studies have shown that variants might have strong age-dependent effects across childhood [15,51]. We performed a sensitivity analyses excluding children aged <6 years, as the approximate age of the adiposity rebound [30]. However, no difference in main results with the full meta-analysis were observed.

Observational studies suggest that childhood obesity is not only related to several anthropometric and cardio-metabolic phenotypes in later life, such as type 2 diabetes, but also to respiratory and neurocognitive traits, including asthma and Parkinson's disease and to a lower age at menarche [14–18,52–58]. In observational studies, effect estimates may be influenced by confounding factors and reverse causation, potentially evoking spurious associations [59,60].
Genetic studies can provide more insight into the etiology of complex diseases. We observed a strong positive genetic correlation of childhood BMI with adult BMI. This is in line with previous studies [8,20,21]. We additionally observed positive genetic correlations between childhood BMI and several cardio-metabolic phenotypes in later life, including waist-to-hip ratio, diastolic blood pressure, type 2 diabetes, and coronary artery disease. Negative genetic correlations were found between childhood BMI and HDL-C and age at menarche. These results may suggest that the associations reported in observational studies are partly explained by genetic factors [15-17,58,61]. However, there is also evidence from previous work to support that the associations of childhood BMI with cardiometabolic phenotypes in adulthood are explained by the continuity of a high BMI from childhood until later ages, rather than by an independent effect of childhood BMI on adult cardiometabolic phenotypes [15,62]. From our data, we are not able to distinguish this. Childhood BMI was not genetically correlated with asthma and Parkinson's disease. This may indicate that the observational associations between childhood BMI and these phenotypes are not strongly explained by shared genetics [18,56,57,63].

The GRS combining the 25 top SNPs was not associated with cardiometabolic phenotypes in children aged 6 years. This may indicate that there is no shared genetic basis between childhood BMI and these phenotypes in childhood. However, the GRS analyses in children had a much lower sample size than the LD score regression analyses in adults and phenotypic variation in these phenotypes is more limited in children, leading to a much lower power to detect associations in these analyses. Additionally, the GRS was composed of the top-associated SNPs, whereas the genetic correlation estimated from the LDSR examined variation genome-wide.

We observed a SNP heritability of 0.23 which is consistent with previous findings [4-6].

Secular trends in obesity across populations and age groups can influence the heritability estimates across distinct population settings, requiring careful interpretation. This contention is also relevant for the interpretation of the genetic correlations estimated between traits. Environmental influences like those giving rise to the increase obesity in the last decades, can influence heritability estimates and hence, the power to identify significant genetic correlations. Before concluding unequivocal absence of some degree of “shared heritability” between childhood BMI and some of the adult traits, genetic correlations should be interpreted in the context of power limitations. Increasingly larger environmental influences along the life-course can result in lower heritability, but recent work has also shown that the increase in phenotypic variance accompanying increasing prevalence of obesity occurs alongside an increase in genetic variance [64-67]. This results in relatively stable (broad sense) heritability estimates across measurement years, as recently shown by a large-scale meta-analysis of adult twin data.

The 25-SNP GRS was positively associated with both childhood and adult BMI, showing slightly larger effect estimates in children suggesting that these specific genetic variants affect BMI in both childhood and adulthood, but with stronger effects at younger ages. A recent study, using genome-wide polygenic scores of 2.1 million common variants, found that the overall effect of those variants on weight starts in early childhood and increases over time [4]. Two previous studies also describe specific genetic variants associated with BMI in infancy only, and overlapping patterns of genetic variants with those in adults emerging from childhood onwards. Three SNPs associated with infant BMI from these studies were not genome-wide significantly associated with childhood BMI in our data (P-values >0.02), which supports their infancy-specific effects [68,69].

Although many of the associated variants from the current study overlap between children and adults, the relative order of the signals differs. Additionally, SLC45A3, one of the novel loci did not show genome-wide association in adult data [7]. However, suggestive association of this locus with adult BMI was observed (P-value = 2.7 x 10⁻⁵). Overall, the effect estimates of the 25 SNPs in childhood were highly correlated with those in adulthood (r² = 0.86). Taken
together, evidence from the current and previous studies suggests that biological processes
underlying BMI are similar from childhood onwards, but their relative influence may differ
depending on the life stage.

Conclusions

In conclusion, we identified 25 loci for childhood BMI, together explaining 3.6% of the variance in childhood BMI. Two of these are novel and four represent independent SNPs at loci known to be associated with adult or childhood BMI. A strong positive genetic correlation of childhood BMI with adult BMI and related cardio-metabolic phenotypes was observed. Our results suggest that the biological processes underlying childhood BMI largely, but not completely, overlap with those underlying adult BMI. The well-known observational associations of BMI in childhood with cardio-metabolic diseases in adulthood may reflect partial genetic overlap, but in light of previous evidence, it is also likely that they are explained through phenotypic continuity of BMI from childhood into adulthood.

**Materials and methods**

Ethics statement

All individual studies got approval by their medical ethics review committees. All participants gave written informed consent. Study-specific ethics statements are given in S1 Text.

Study design

We conducted a two-stage meta-analysis in children of European ancestry to identify genetic loci associated with childhood BMI. Sex- and age-adjusted standard deviation scores were cre- ated for BMI at the latest time point (oldest age, if multiple measurements were available) between 2 and 10 years using the same software and external reference across all studies (LMS growth; Pan H, Cole TJ, 2012; http://www.healthforallchildren.co.uk). In the case of twin pairs and siblings, only one of each twin or sibling pair was included, either randomly or based on genotyping or imputation quality.

In the discovery stage, we performed a meta-analysis of 26 studies (N = 39,620), including the Avon Longitudinal Study of Parents and Children (ALSPAC, N = 6790), the Bone Mineral Density in Childhood Study (BMDCS, N = 543), BRain dEvelopment and Air polluTion ultra-fine particles in scHool childrEn (BREATHE, N = 1633), the Children's Hospital of Philadelphia (CHOP, N = 5488), the Copenhagen Prospective Studies on Asthma in Childhood 2000 (COPSAC2000, N = 327), the Copenhagen Prospective Studies on Asthma in Childhood 2010 (COPSAC2010, N = 571), the Danish National Birth Cohort- preterm birth study (DNBC-PTB, N = 1007), the French Young Study (French young, N = 304 cases and N = 144 controls), the Generation R Study (GenerationR, N = 2071), the Genetics of Overweight Young Adults (GOYA Male, N = 319), the Helsinki Birth Cohort Study (HBCS, N = 1575), the INfancia y Medio Ambiente [Environment and Childhood] Project, with two subcohorts that were entered into the meta-analysis seperately (INMA-Sabadell and Valencia subcohort, N = 650, and INMA-Menorca subcohort, N = 300), German Infant Study on the influence of Nutrition Intervention PLUS environmental and genetic influences on allergy development & Influence of life-style factors on the development of the immune system and allergies in East and West Germany (GINIplus&LISA, N = 1471), the Manchester Asthma and Allergy Study (MAAS, N = 784), the Norwegian Mother, Father and Child Cohort (MoBa, N = 522), the Northern Finland Birth Cohort 1966 (NFBC 1966, N = 3949), the Northern Finland Birth Cohort 1986 (NFBC 1986, N = 1056), the Netherlands Twin Register (NTR, N = 1767), the
Physical Activity and Nutrition in Children Study (PANIC, N = 374), The Danish Childhood Obesity Data and Biobank (TDCOB, N = 158 controls), the Raine Study (the Raine Study (Generation 2), N = 1458), the Special Turku Coronary Risk factor Intervention Project (STRIP, N = 551), the Young Finns Study (YFS, N = 1134), the TEENs of Attica: Genes and Environment (TEENAGE, N = 252), and the British 1958 Birth Cohort Study (1958BC-T1DGC, N = 2081 and 1958BC-WTCCC, N = 2341).

In the replication stage, we included 14 studies (n = 21,491), which did not have genome-wide association data available at the time of the discovery analysis: 888 additional children from the Danish National Birth Cohort- Goya offspring (N = 407 offspring from obese mothers, N = 481 offspring from randomly selected mothers), 294 additional children from the INfancia y Medio Ambiente [Environment and Childhood] (INMA) Project (INMA- Gipuzkoa subcohort, N = 314), 6,828 additional children from the Norwegian Mother, Father and Child Cohort (MoBa, N = 6828), 753 additional children from TDCOB (N = 344 controls and N = 409 cases), the Amsterdam Born Children and their Development- Genetic Enrichment (ABCD, N = 1154), The European Childhood Obesity Project (CHOP Study, N = 369), the Family Atherosclerosis Monitoring In earLY life (FAMILY) study (the FAMILY study, N = 543), the Etude des Déterminants pré- et postnatals précoces du développement et de la santé de l'Enfant (EDEN) mother-child cohort (EDEN, N = 821), the Exeter Family Study of Childhood Health (EFSOCH, N = 542), the Leipzig Research Center for Civilization Diseases —Child study (LIFE-Child, N = 1318), the Prevention and Incidence of Asthma and Mite Allergy study (PIAMA, N = 1958), the Screening of older women for prevention of fracture Study (SCOOP, N = 685), the Småbørns Kost Og Trivsel study (SKOT1, N = 236), the Twin Early Development Study (TEDS, N = 3933), the Tracking Adolescents' Individual Lives Survey cohort (TRAILS-population cohort, N = 1169). In the EDEN mother-child cohort, information was available about three SNPs only (rs7138803, rs13107325, and rs987237).

Characteristics of discovery and replication studies can be found in S1 Table and S1 Text.

## Study-level analyses

Genome-wide association analyses were first run in all discovery cohorts separately. Studies used high-density Illumina or Affymetrix SNP arrays, followed by imputation to the 1000 Genomes Project or HRC. Before imputation, studies applied study specific quality filters on sample and SNP call rate, minor allele frequency and Hardy-Weinberg disequilibrium (see S1 Table for details). Linear regression models assuming an additive genetic model were run in each study to assess the association of each SNP with BMI SDS, adjusting for principal components if this was deemed needed in the individual studies. As BMI SDS is age and sex specific, no further adjustments were made. Before the meta-analysis, we applied quality filters to each study, filtering out SNPs with a minor allele frequency (MAF) below 1% and SNPs with poor imputation quality (MACH r2Hat ≤0.3, IMPUTE proper_info ≤0.4 or info ≤0.4).

## Meta-analysis

We performed fixed-effects inverse-variance weighted meta-analysis of all discovery samples using Metal [70]. Genomic control was applied to every study before the meta-analysis. Individual study lambdas before genomic control ranged from 0.993 to 1.036 (S1 Table). The lambda of the discovery meta-analysis is shown in S1 Fig. After the meta-analysis, we excluded SNPs for which information was available in less than 50% of the studies and less than 50% of the total sample size. We report I2 and p-value for heterogeneity for all findings.

The final dataset consisted of 8,228,795 autosomal SNPs. Genome-wide Complex Trait Analysis (GCTA) was used to select the independent SNPs for each locus [22]. We performed
conditional analyses based on summary-level statistics and LD estimation between SNPs from the Generation R Study as a reference sample to select independently associated SNPs based on conditional P-values [22]. Forty-seven genome-wide significant or suggestive loci (P-values <5 × 10⁻⁸ and <5 × 10⁻⁶, respectively) were taken forward for replication in 14 replication cohorts. Fixed-effects inverse variance meta-analysis was performed for these 47 SNPs combining the discovery samples and all replication samples, giving a combined analysis beta, standard error and P-value (Table 1). SNPs that reached genome-wide significance in the combined analysis were considered to be genome-wide significant.

## Functional mapping and annotation of genetic associations (FUMA)

To obtain predicted functional consequences for our 25 SNPs, we used SNP2FUNC in FUMA, a web-based platform to facilitate and visualize functional annotation of GWAS results (http://fuma.ctglab.nl) [3]. By matching chromosome, position, and reference and alternative alleles, combined annotation-dependent depletion (CADD) scores were annotated, indicating the del-teriousness of a SNP [37].

To annotate the nearest genes of the 25 SNPs in biological context, we used the GENE2-FUNC option in FUMA, which provides hypergeometric tests of enrichment of a list of genes in 53 GTEx tissue-specific gene expression sets (GTEx v 7) [3,31]. We used GENE2FUNC for two sets of genes: 1. Nearest genes of 25 SNPs; 2. Genes located in a region of 500 kb to either side of the 25 SNPs.

## Look-up of the 25 SNPs in expression data

We studied the associations of the 25 SNPs associated with childhood BMI with gene expression levels in adipose tissue samples from the Leipzig Adipose Tissue Childhood Cohort [32]. These associations were examined in the following tissues: whole adipose tissue, isolated adipocytes and isolated stroma-vascular cells using genome-wide expression analysis (Illumina HumanHT-12 v4 arrays). Gene expression raw data of all 47,231 probes was extracted by Illumina GenomeStudio without additional background correction. Data was further processed within R / Bioconductor. Expression values were log2-transformed and quantile-normalised [71,72]. Batch effects of expression BeadChips were corrected using an empirical Bayes method [73]. Within pre-processing, gene-expression probes detected by Illumina GenomeStudio as expressed in less than 5% of the samples were excluded as well as probes still found to be significantly associated with batch effects after Bonferroni-correction. Furthermore, gene-expression probes with poor mapping on the human transcriptome [74] were also excluded. In summary, these filters resulted in 23354, 21258, and 22637 valid gene-expression probes from which 20672, 18956, and 20230 probes corresponded to 14455, 13518, and 14256 genes mapping to a unique position in the human genome (hg19) for whole adipose tissue, adipocytes, and stroma/vascular cells, respectively. Three criteria were used to remove samples of low quality: First, the number of detected gene-expression probes of a sample was required to be within ± 3 interquartile ranges (IQR) from the median. Second, the Mahalanobis distance of several quality characteristics of each sample had to be lower than median + 4 x IQR. Third, Euclidean distances of expression values as described [71] had to be lower than median + 4 x IQR. Overall, of the assayed samples, 2, 4, and 2 samples were excluded for quality reasons leaving 203, 63 and 69 unique individuals having also valid data for eQTL analysis for whole adipose tissue, adipocytes, and stroma/vascular cells, respectively. Associations between the genotype and gene expression of genes in *cis* (respective gene area +/- 1 Mb regarding transcription start and transcription end) were analyzed using a gene-dose based linear regression model adjusted for age and sex as implemented in MatrixEQTL [75]. Analysis of variants
within one haplotype were done through analyses of linkage-disequilibrium using 1000 Genomes Phase 1 Version 3 and HapMap r28 hg19 CEU as references.

## Colocalization analysis

We used Bayesian colocalization analysis to examine evidence for colocalization between childhood BMI and eQTL signals (GTEx v7). Colocalization analyses were conducted using the R package coloc, hht://cran.r-project.org/web/packages/coloc, as described previously [33]. Briefly, in each of the GTEx v7 tissues, all cis-eQTLs at FDR <5% were identified. For each eQTL, GWAS summary statistics were extracted for all SNPs that were present in >50% of the studies and >50% of the total sample size and that were in common to both GWAS and eQTL studies, within 1 MB of the transcription start site of the gene. For each such locus, colocalization analyses were done with default parameters, testing the following hypotheses [33]:

$H_0$: No association with either trait;

$H_1$: Association with childhood BMI only;

$H_2$: Association with gene expression only;

$H_3$: Association with childhood BMI and gene expression, two distinct causal variants;

$H_4$: Association with childhood BMI and gene expression, one shared causal variant.

Support for each hypothesis was quantified in terms of posterior probabilities, defined at SNP level and indicated by $PP_0, PP_1, PP_2, PP_3$ or $PP_4$, corresponding to the five hypotheses and measuring how likely these hypotheses were. S10B Table shows the above-mentioned posterior probabilities for all pairs. In most pairs, no evidence for association was found with either trait. In case association was observed, it was mostly with a single trait. To define colocalization we used restriction to pairs of childhood BMI and eQTL signals with a high posterior probability for colocalization, indicated by a $PP_4/(PP_3+PP_4) >0.9$ (S10A Table).

## DAVID

To explore biological processes, we used DAVID, with the 25 nearest genes as input, using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database [34,35].

## Linkage-disequilibrium score regression

The use of LD score regression to estimate genetic correlations between two phenotypes has been described in detail previously [20]. Briefly, LD score is a measure of how much a genetic variation is tagged by each variant. A high LD score indicates that a variant is in high LD with many nearby polymorphisms. Variants with high LD scores are more likely to contain true signals and have a higher chance of overlap with genuine signals between GWAS. To estimate LD scores, summary statistics from GWAS meta-analysis are used to calculate the cross-product of test statistics per SNP, which is regressed on the LD score. The slope of the regression is a function of the genetic covariance between traits [20]:

$$E(z1jz2j) = \frac{\sqrt{N1N2\rho_g}}{M} l_j + \frac{\rho N_s}{\sqrt{N1N2}}$$

where $N_i$ is the sample size of study i, $\rho_g$ is the genetic covariance, $M$ is the number of SNPs in the reference panel with a MAF between 5% and 50%, $l_j$ is the LD score for SNP j, $N_s$ quantifies the number of individuals that overlap both studies, and $\rho$ is the phenotypic correlation amongst the $N_s$ of overlapping samples. A sample overlap or cryptic relatedness between
samples will only affect the intercept from the regression but not the slope. Estimates of genetic covariance will therefore not be biased by overlapping samples. Similarly, in case of population stratification, the intercept will be affected but it will have only minimal impact on the slope since population stratification does not correlate with LD between variants.

Because of the correlation between the imputation quality and LD score, imputation quality is a confounder for LD score regression. Therefore, SNPs were excluded for the following reasons: MAF <0.01 and INFO ≤0.9. The filtered GWAS results were uploaded on http://ldsc.broadinstitute.org/ldhub/, a website with many GWAS meta-analyses available on which LD score regression has been implemented by the developers of the LD score regression method. In case multiple GWAS meta-analyses were available for the same phenotype, the genetic correlation with childhood BMI was estimated using the most recent meta-analysis. Genetic correlations are shown in Fig 2 and S11 Table.

## Genetic risk score and percentage of variance explained

We combined the 25 genome-wide significant SNPs from the combined meta-analysis into a GRS by summing up the number of BMI SDS-increasing alleles, weighted by the effect sizes from the combined meta-analysis. The GRS was rescaled to a range from 0 to 50, which is the maximum number of BMI SDS increasing alleles and rounded to the nearest integer. Linear regression analysis was used to examine the associations of the risk score with childhood and adult BMI. For these analyses data from the TRAILS cohort (N = 1169), one of the largest replication cohorts, and data from the Rotterdam Study (RS-I-1; n = 5,957, RS-II-1; n = 2,147 and RS-III-1; n = 2,998) were used. Additionally, linear regression analysis was used to examine the associations of the GRS with birth weight and childhood metabolic phenotypes in Generation R in which detailed information on these phenotypes was available. When calculating the risk score for the TRAILS cohort and Generation R, effect estimates from the combined meta-analysis were used after excluding TRAILS and Generation R, respectively, from the meta-analysis. The variance explained was estimated by the adjusted R² of the models.

## Supporting information

<caption>S1 Table. Descriptive characteristics of studies—Separate Excel sheets.</caption>

(XLSX)

<caption>S2 Table. Results of the discovery, replication and combined analyses for all 47 independent loci with *P*-values <5 x 10<sup>-6</sup> in the discovery analysis—Separate Excel sheet.</caption>

(XLSX)

<caption>S3 Table. Directions of effect for the individual discovery and replication studies for all 47 loci with *P*-values <5 x 10<sup>-6</sup> in the discovery analysis.</caption>

(DOCX)

<caption>S4 Table. Results of the discovery analysis for all SNPs with *P*-values <5 x 10<sup>-6</sup>—Separate Excel sheet.</caption>

(XLSX)

<caption>S5 Table. Descriptive characteristics of studies with data on fat mass percentage available—Separate Excel sheets.</caption>

(XLSX)

<caption>S6 Table. Results for the 25 SNPs with *P*-values <5 x 10<sup>-8</sup> in the combined analysis, restricted to children aged ≥ 6 years.</caption>

(DOCX)
S7 Table. Results of the combined analyses for the 25 SNPs excluding case-control studies and excluding studies with a sample size smaller than n = 500.

(DOCX)

S8 Table. Look up of the 25 nearest genes and candidate genes indicated by colocalization analysis- *Separate Excel sheet*.

(XLSX)

S9 Table. Expression levels of the 25 childhood BMI SNPs in the Adipose Tissue Childhood Cohort- *Separate Excel sheet*.

(XLSX)

S10 Table. Colocalization analyses with GTEx v7—*Separate Excel sheet*.

(XLSX)

S11 Table. Summary results of LD score regression analyses between childhood BMI and various metabolic, cardiovascular and lifestyle-related phenotypes.

(DOCX)

S12 Table. Results of adult BMI GWAS meta-analysis for 20 loci with P-values <5 x 10⁻⁸ in the combined analysis.

(DOCX)

S13 Table. Results of the GRS for birth weight and childhood cardio-metabolic phenotypes.

(DOCX)

S1 Fig. Quantile-Quantile plot of the SNPs in the discovery meta-analysis of 26 studies.

(TIF)

S2 Fig. Locus zoom plots of the 4 independent SNPs at known loci.

(TIF)

S1 Text. Cohort information.

(DOCX)

## Author Contributions

**Conceptualization:** Suzanne Vogelezang, Niels Grarup, Momoko Horikoshi, Jorma Viikari, Ioanna Ntalla, Peter Rzehak, Klaus Bønnelykke, Lawrence J. Beilin, Ricardo Closa-Monas-terolo, Joaquin Escribano, Veit Grote, Dariusz Gruszfeld, Hakon Hakonarson, Torben Hansen, Elina Hyppönen, Wieland Kiess, Berthold Koletzko, Kathrin Landgraf, Jean-Paul Langhendries, Mads Melbye, Kim F. Michaelsen, Trevor A. Mori, Harri Niinikoski, Oluf Pedersen, Craig E. Pennell, Peter D. Sly, Elvira Verduci, Sylvain Sebert, Ellen A. Nohr, Tanja G. M. Vrijkotte, George D. Dedoussis, Philippe Froguel, Thorkild I. A. Sørensen, Babette S. Zemel, Olli Raitakari, Patrick G. Holt, Joachim Heinrich, Dorret I. Boomsma, Marjo-Riitta Järvelin, Antje Körner, Jens-Christian Holm, Hans Bisgaard, Vincent W. V. Jaddoe, Struan F. A. Grant, Janine F. Felix.

**Data curation:** Suzanne Vogelezang, Tarunveer S. Ahluwalia, John A. Curtin, Timo A. Lakka, Niels Grarup, Markus Scholz, Peter J. van der Most, Claire Monnereau, Evie Stergiakouli, Anni Heiskala, Marie Standl, Jorma Viikari, Frank Geller, Carmen Íñiguez, Ioanna Ntalla, Øyvind Helgeland, Mohammed H. Zafarmand, Barbara Heude, Christian T. Have, Theresia
M. Schnurr, Klaus Bønnelykke, Lawrence J. Beilin, Marie-Aline Charles, Bo Chawes, Karine Clément, Ricardo Closa-Monasterolo, Adnan Custovic, Johan G. Eriksson, Joaquin Escribano, Veit Grote, Dariusz Gruszfeld, Hakon Hakonarson, Torben Hansen, Andrew T. Hattersley, Mette Hollensted, Jouke-Jan Hottenga, Elina Hyppönen, Stefan Johansson, Mika Kähönen, Wieland Kiess, Bridget A. Knight, Berthold Koletzko, Kathrin Landgraf, Jean-Paul Langhendries, Terho Lehtimäki, Virpi Lindi, Mariona Bustamante, Carolina Medina-Gomez, Mads Melbye, Kim F. Michaelsen, Camilla S. Morgen, Trevor A. Mori, Tenna R. H. Nielsen, Harri Niinikoski, Albertine J. Oldehinkel, Katja Pahkala, Oluf Pedersen, Craig E. Pennell, Christine Power, Sijmen A. Reijneveld, Fernando Rivadeneira, Angela Simpson, Peter D. Sly, Jakob Stokholm, Kook K. Teo, Elisabeth Thiering, Nicholas J. Timpson, André G. Uitterlinden, Catharina E. M. van Beijsterveldt, Marc Vaudel, Elvira Verduci, Rebecca K. Vinding, Mandy Vogel, Eleftheria Zeggini, Sylvain Sebert, Mads V. Lind, Loreto Santa-Marina, Eva Reischl, Christine Frithioff-Bøjsøe, David Meyre, Eleanor Wheeler, Ken Ong, Ellen A. Nohr, Tanja G. M. Vrijkotte, Robert Plomin, Pål R. Njølstad, George D. Dedoussis, Thorkild I. A. Sørensen, Rachel M. Freathy, Babette S. Zemel, Olli Raitakari, Martine Vrijheid, Bjarke Feenstra, Harold Snieder, Patrick G. Holt, Joachim Heinrich, Elisabeth Widén, Jordi Sunyer, Marjo-Riitta Järvelin, Antje Körner, George Davey Smith, Jens-Christian Holm, Mustafa Atalay, Clare Murray, Hans Bisgaard, Mark I. McCarthy, Vincent W. V. Jaddoe, Struan F. A. Grant, Janine F. Felix.

**Formal analysis:** Suzanne Vogelezang, Jonathan P. Bradfield, Tarunveer S. Ahluwalia, John A. Curtin, Niels Grarup, Markus Scholz, Peter J. van der Most, Claire Monnereau, Evie Stergiakouli, Anni Heiskala, Momoko Horikoshi, Iryna O. Fedko, Natalia Vilor-Tejedor, Diana L. Cousminer, Carol A. Wang, Frank Geller, Niina Pitkänen, Alessandra Chesi, Loic Yengo, Ioanna Ntalla, Øyvind Helgeland, Saskia Selzam, Judith M. Vonk, Mohammed H. Zafarmand, Barbara Heude, Akram Alyass, Robin N. Beaumont, Christian T. Have, Peter Rzehak, Jose Ramon Bilbao, Theresia M. Schnurr, Lisbeth Carstensen, Maria Groen-Blokhuis, Raimo Joro, Ville Karhunen, Andreas Kühnapfel, Jaakko T. Leinonen, Aihuali Li, Estelle Lowry, Carolina Medina-Gomez, Kalliope Panoutsopoulou, Elisabeth Thiering, Barbera D. C. van Schaik, Marc Vaudel, Eleftheria Zeggini, Eleanor Wheeler, Bjarke Feenstra, Leo-Pekka Lyytikäinen, Holger Kirsten, Janine F. Felix.

**Funding acquisition:** Natalia Vilor-Tejedor, Diana L. Cousminer, Maties Torrent, Ismaa Sadaf Farooqi, Inês Barroso, Andrew T. Hattersley, Mette Hollensted, Berthold Koletzko, Jaakko T. Leinonen, Christine Power, Loreto Santa-Marina, Eva Reischl, Christine Frithioff-Bøjsøe, David Meyre, Ken Ong, Ellen A. Nohr, Tanja G. M. Vrijkotte, Gerard H. Kopelman, Robert Plomin, Pål R. Njølstad, George D. Dedoussis, Philippe Froguel, Thorkild I. A. Sørensen, Bo Jacobsson, Rachel M. Freathy, Babette S. Zemel, Olli Raitakari, Martine Vrijheid, Bjarke Feenstra, Leo-Pekka Lyytikäinen, Harold Snieder, Holger Kirsten, Patrick G. Holt, Joachim Heinrich, Elisabeth Widén, Jordi Sunyer, Dorret I. Boomsma, Marjo-Riitta Järvelin, Antje Körner, George Davey Smith, Jens-Christian Holm, Mustafa Atalay, Clare Murray, Hans Bisgaard, Mark I. McCarthy, Vincent W. V. Jaddoe, Struan F. A. Grant.

**Investigation:** Suzanne Vogelezang, Janine F. Felix.

**Methodology:** Suzanne Vogelezang, Jonathan P. Bradfield, Peter J. van der Most, Claire Monnereau, Fernando Rivadeneira, Christopher D. Brown, Janine F. Felix.

**Project administration:** Suzanne Vogelezang, Janine F. Felix.

**Resources:** Janine F. Felix.
**Supervision:** Maties Torrent, Loreto Santa-Marina, Eva Reischl, Christine Frithioff-Bojsøe, David Meyre, Ken Ong, Ellen A. Nohr, Tanja G. M. Vrijkotte, Gerard H. Koppelman, Robert Plomin, Pål R. Njølstad, George D. Dedoussis, Philippe Froguel, Thorkild I. A. Sørensen, Bo Jacobsson, Rachel M. Freathy, Babette S. Zemel, Olli Raitakari, Martine Vrijheid, Bjarke Feenstra, Leo-Pekka Lyytikäinen, Harold Snieder, Holger Kirsten, Patrick G. Holt, Joachim Heinrich, Elisabeth Widén, Jordi Sunyer, Dorret I. Boomsma, Marjo-Riitta Järvelin, Antje Körner, George Davey Smith, Jens-Christian Holm, Mustafa Atalay, Clare Murray, Hans Bisgaard, Mark I. McCarthy, Vincent W. V. Jaddoe, Struan F. A. Grant, Janine F. Felix.

**Visualization:** Jonathan P. Bradfield, Carolina Medina-Gomez.

**Writing – original draft:** Suzanne Vogelezang, Vincent W. V. Jaddoe, Janine F. Felix.

**Writing – review & editing:** Jonathan P. Bradfield, Tarunveer S. Ahluwalia, John A. Curtin, Timo A. Lakka, Niels Grarup, Markus Scholz, Peter J. van der Most, Claire Monnereau, Evie Stergiakouli, Anni Heiskala, Momoko Horikoshi, Iryna O. Fedko, Natalia Vilor-Tejedor, Diana L. Cousminer, Marie Standl, Carol A. Wang, Jorma Viikari, Frank Geller, Carmen Íñiguez, Niina Pitkänen, Alessandra Chesi, Jonas Bacelis, Loic Yengo, Maties Torrent, Ioanna Ntalla, Øyyvind Helgeland, Saskia Selzam, Judith M. Vonk, Mohammed H. Zafarmand, Barbara Heude, Ismaa Sadaf Farooqi, Akram Alyass, Robin N. Beaumont, Christian T. Have, Peter Rzehak, Jose Ramon Bilbao, Theresia M. Schnurr, Inês Barroso, Klaus Bønnelykke, Lawrence J. Beilin, Lisbeth Carstensen, Marie-Aline Charles, Bo Chawes, Karine Clément, Ricardo Closa-Monasterolo, Adnan Custovic, Johan G. Eriksson, Joaquin Escribano, Maria Groen-Blokhuis, Veit Grote, Dariusz Gruszfeld, Hakon Hakonarson, Torben Hansen, Andrew T. Hattersley, Mette Hollensted, Jouke-Jan Hottenga, Elina Hyppönen, Stefan Johansson, Raimo Joro, Mika Kähönen, Ville Karhunen, Wieland Kiess, Bridget A. Knight, Berthold Koletzko, Andreas Kühnapfel, Kathrin Landgraf, Jean-Paul Langhendries, Terho Lehtimäki, Jaakko T. Leinonen, Aihuali Li, Virpi Lindi, Estelle Lowry, Mariona Bustamante, Carolina Medina-Gomez, Mads Melbye, Kim F. Michaelsen, Camilla S. Morgen, Trevor A. Mori, Tenna R. H. Nielsen, Harri Niinikoski, Albertine J. Oldehinkel, Katja Pahkala, Kalliope Panoutsopoulou, Oluf Pedersen, Craig E. Pennell, Christine Power, Sijmen A. Reijneveld, Fernando Rivadeneira, Angela Simpson, Peter D. Sly, Jakob Stokholm, Kook K. Teo, Elisabeth Thiering, Nicholas J. Timpson, André G. Uitterlinden, Catharina E. M. van Beijsterveldt, Barbera D. C. van Schaik, Marc Vaudel, Elvira Verduci, Rebecca K. Vinding, Mandy Vogel, Eleftheria Zeggini, Sylvain Sebert, Mads V. Lind, Christopher D. Brown, Loreto Santa-Marina, Eva Reischl, Christine Frithioff-Bojsøe, David Meyre, Eleanor Wheeler, Ken Ong, Ellen A. Nohr, Tanja G. M. Vrijkotte, Gerard H. Koppelman, Robert Plomin, Pål R. Njølstad, George D. Dedoussis, Philippe Froguel, Thorkild I. A. Sørensen, Bo Jacobsson, Rachel M. Freathy, Babette S. Zemel, Olli Raitakari, Martine Vrijheid, Bjarke Feenstra, Leo-Pekka Lyytikäinen, Harold Snieder, Holger Kirsten, Patrick G. Holt, Joachim Heinrich, Elisabeth Widén, Jordi Sunyer, Dorret I. Boomsma, Marjo-Riitta Järvelin, Antje Körner, George Davey Smith, Jens-Christian Holm, Mustafa Atalay, Clare Murray, Hans Bisgaard, Mark I. McCarthy, Vincent W. V. Jaddoe, Struan F. A. Grant, Janine F. Felix.

**References**

1. NCD Risk Factor Collaboration. Worldwide trends in body-mass index, underweight, overweight, and obesity from 1975 to 2016: a pooled analysis of 2416 population-based measurement studies in 128.9 million children, adolescents, and adults. Lancet. 2017; 390(10113):2627–42. https://doi.org/10.1016/S0140-6736(17)32129-3 PMID: 29029897

2. Silventoinen K, Jelenkovic A, Sund R, Hur YM, Yokoyama Y, Honda C, et al. Genetic and environmental effects on body mass index from infancy to the onset of adulthood: an individual-based pooled analysis
of 45 twin cohorts participating in the COllaborative project of Development of Anthropometrical measures in Twins (CODATwins) study. Am J Clin Nutr. 2016; 104(2):371-9. https://doi.org/10.3945/ajcn.116.130252 PMID: 27413137

3. Silventoinen K, Jelenkovic A, Sund R, Yokoyama Y, Hur YM, Cozen W, et al. Differences in genetic and environmental variation in adult BMI by sex, age, time period, and region: an individual-based pooled analysis of 40 twin cohorts. Am J Clin Nutr. 2017; 106(2):457-66. https://doi.org/10.3945/ajcn.117.153643 PMID: 28679550

4. Khera AV, Chaffin M, Wade KH, Zahid S, Brancale J, Xia R, et al. Polygenic Prediction of Weight and Obesity Trajectories from Birth to Adulthood. Cell. 2019; 177(3):587-96 e9. https://doi.org/10.1016/j.cell.2019.03.028 PMID: 31002795

5. Yang J, Bakshi A, Zhu Z, Hemani G, Vinkhuyzen AA, Lee SH, et al. Genetic variance estimation with imputed variants finds negligible missing heritability for human height and body mass index. Nat Genet. 2015; 47(10):1114-20. https://doi.org/10.1038/ng.3390 PMID: 26323059

6. Yang J, Manolio TA, Pasquale LR, Boerwinkle E, Caporaso N, Cunningham JM, et al. Genome partitioning of genetic variation for complex traits using common SNPs. Nat Genet. 2011; 43(6):519-25. https://doi.org/10.1038/ng.823 PMID: 21552263

7. Yengo L, Sidorenko J, Kemper KE, Zheng Z, Wood AR, Weedon MN, et al. Meta-analysis of genome-wide association studies for height and body mass index in approximately 700000 individuals of European ancestry. Hum Mol Genet. 2018.

8. Felix JF, Bradfield JP, Monnereau C, van der Valk RJ, Stergiakouli E, Chesi A, et al. Genome-wide association analysis identifies three new susceptibility loci for childhood body mass index. Hum Mol Genet. 2016; 25(2):389-403. https://doi.org/10.1093/hmg/ddv472 PMID: 26604143

9. Berndt SI, Gustafsson S, Magi R, Ganna A, Wheeler E, Feitosa MF, et al. Genome-wide meta-analysis identifies 11 new loci for anthropometric traits and provides insights into genetic architecture. Nat Genet. 2013; 45(5):501-12. https://doi.org/10.1038/ng.2606 PMID: 23563607

10. Locke AE, Kahali B, Berndt SI, Justice AE, Pers TH, Day FR, et al. Genetic studies of body mass index yield new insights for obesity biology. Nature. 2015; 518(7538):197-206. https://doi.org/10.1038/nature14177 PMID: 25673413

11. Winkler TW, Justice AE, Graff M, Barata L, Feitosa MF, Chu S, et al. Correction: The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study. PLoS Genet. 2016; 12(6):e1006166. https://doi.org/10.1371/journal.pgen.1006166 PMID: 27355579

12. Haworth S, Shapland CY, Hayward C, Prins BP, Felix JF, Medina-Gomez C, et al. Low-frequency variation in TP53 has large effects on head circumference and intracranial volume. Nat Commun. 2019; 10(1):357. https://doi.org/10.1038/s41467-018-07863-x PMID: 30664637

13. Graff M, Ngwa JS, Workalemahu T, Homuth G, Schipf S, Teumer A, et al. Genome-wide analysis of BMI in adolescents and young adults reveals additional insight into the effects of genetic loci over the life course. Hum Mol Genet. 2013; 22(17):3597-607. https://doi.org/10.1093/hmg/ddt205 PMID: 23669352

14. Singh AS, Mulder C, Twisk JW, van Mechelen W, Chinapaw MJ. Tracking of childhood overweight into adulthood: a systematic review of the literature. Obes Rev. 2008; 9(5):474-88. https://doi.org/10.1111/j.1467-789X.2008.00475.x PMID: 18331423

15. Bjerregaard LG, Jensen BW, Angquist L, Osler M, Sorensen TIA, Baker JL. Change in Overweight from Childhood to Early Adulthood and Risk of Type 2 Diabetes. N Engl J Med. 2018; 378(14):1302-12. https://doi.org/10.1056/NEJMoa1713231 PMID: 29617589

16. Zhang T, Zhang H, Li Y, Sun D, Li S, Fernandez C, et al. Temporal Relationship Between Childhood Body Mass Index and Insulin and Its Impact on Adult Hypertension: The Bogalusa Heart Study. Hypertension. 2016; 68(3):818-23. https://doi.org/10.1161/HYPERTENSIONAHA.116.07991 PMID: 27432860

17. Sinaiko AR, Donahue RP, Jacobs DR Jr., Prineas RJ. Relation of weight and rate of increase in weight during childhood and adolescence to body size, blood pressure, fasting insulin, and lipids in young adults. The Minneapolis Children's Blood Pressure Study. Circulation. 1999; 99(11):1471-6. https://doi.org/10.1161/01.cir.99.11.1471 PMID: 10086972

18. Mamun AA, Lawlor DA, Alati R, O'Callaghan MJ, Williams GM, Najman JM. Increasing body mass index from age 5 to 14 years predicts asthma among adolescents: evidence from a birth cohort study. Int J Obes (Lond). 2007; 31(4):578-83.

19. Zimmermann E, Bjerregaard LG, Gamborg M, Vaag AA, Sorensen TIA, Baker JL. Childhood body mass index and development of type 2 diabetes throughout adult life-A large-scale danish cohort study. Obesity (Silver Spring). 2017; 25(5):965-71.
20. Bulik-Sullivan B, Finucane HK, Anttila V, Gusev A, Day FR, Loh PR, et al. An atlas of genetic correlations across human diseases and traits. Nat Genet. 2015; 47(11):1236-41. https://doi.org/10.1038/ng.3406 PMID: 26414676

21. Pickrell JK, Berisa T, Liu JZ, Segurel L, Tung JY, Hinds DA. Detection and interpretation of shared genetic influences on 42 human traits. Nat Genet. 2016; 48(7):709-17. https://doi.org/10.1038/ng.3570 PMID: 27182965

22. Yang J, Ferreira T, Morris AP, Medland SE, Genetic Investigation of ATC, Replication DIG, et al. Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits. Nat Genet. 2012; 44(4):369-75, S1-3. https://doi.org/10.1038/ng.2213 PMID: 22426310

23. Bulik-Sullivan BK, Loh PR, Finucane HK, Ripke S, Yang J, Schizophrenia Working Group of the Psychiatric Genomics C, et al. LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat Genet. 2015; 47(3):291-5. https://doi.org/10.1038/ng.3211 PMID: 25642630

24. Staley JR, Blackshaw J, Kamat MA, Ellis S, Surendran P, Sun BB, et al. PhenoScanner: a database of human genotype-phenotype associations. Bioinformatics. 2016; 32(20):3207-9. https://doi.org/10.1093/bioinformatics/btw373 PMID: 27318201

25. Sudlow C, Gallacher J, Allen N, Beral V, Burton P, Danesh J, et al. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med. 2015; 12(3):e1001779. https://doi.org/10.1371/journal.pmed.1001779 PMID: 25826379

26. Bycroft C, Freeman C, Petkova D, Band G, Elliott LT, Sharp K, et al. The UK Biobank resource with deep phenotyping and genomic data. Nature. 2018; 562(7726):203-9. https://doi.org/10.1038/s41586-018-0579-z PMID: 30305743

27. Lu Y, Day FR, Gustafsson S, Buchkovich ML, Na J, Bataille V, et al. New loci for body fat percentage reveal link between adiposity and cardiometabolic disease risk. Nat Commun. 2016; 7:10495. https://doi.org/10.1038/ncomms10495 PMID: 26833246

28. Bradfield JP. Personal communication, June 20, 2019.

29. Higgins JP, Thompson SG. Quantifying heterogeneity in a meta-analysis. Stat Med. 2002; 21(11):1539-58. https://doi.org/10.1002/sim.1186 PMID: 12111919

30. Rolland-Cachera MF, Deheeger M, Bellisle F, Sempe M, Guilloud-Bataille M, Patois E. Adiposity rebound in children: a simple indicator for predicting obesity. Am J Clin Nutr. 1984; 39(1):129-35. https://doi.org/10.1093/ajcn/39.1.129 PMID: 6691287

31. Consortium GTEx. Human genomics. The Genotype-Tissue Expression (GTEx) pilot analysis: multisite gene regulation in humans. Science. 2015; 348(6235):648-60. https://doi.org/10.1126/science.1262110 PMID: 25954001

32. Landgraf K, Rockstroh D, Wagner IV, Weise S, Tauscher R, Schwartze JT, et al. Evidence of early alterations in adipose tissue biology and function and its association with obesity-related inflammation and insulin resistance in children. Diabetes. 2015; 64(4):1249-61. https://doi.org/10.2337/db14-0744 PMID: 25392242

33. Giambartolomei C, Vukcevic D, Schadt EE, Franke L, Hingorani AD, Wallace C, et al. Bayesian test for colocalisation between pairs of genetic association studies using summary statistics. PLoS Genet. 2014; 10(5):e1004383. https://doi.org/10.1371/journal.pgen.1004383 PMID: 24830394

34. UniProt C. UniProt: a hub for protein information. Nucleic Acids Res. 2015; 43 part 2:2001-2012. https://doi.org/10.1093/nar/gku989 PMID: 25348405

35. Huang da W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009; 37(1):1-13. https://doi.org/10.1093/nar/gkn923 PMID: 19033363

36. Dickinson ME, Flenniken AM, Ji X, Teboul L, Wong MD, White JK, et al. High-throughput discovery of novel developmental phenotypes. Nature. 2016; 537(7621):508-14. https://doi.org/10.1038/nature19356 PMID: 27626380

37. Kircher M, Witten DM, Jain P, O'Reak BJ, Cooper GM, Shendure J. A general framework for estimating the relative pathogenicity of human genetic variants. Nat Genet. 2014; 46(3):310-5. https://doi.org/10.1038/ng.2892 PMID: 24487276

38. Ikram MA, Brusselle GGO, Murad SD, van Duijn CM, Franco OH, Goedegebure A, et al. The Rotterdam Study: 2018 update on objectives, design and main results. Eur J Epidemiol. 2017; 32(9):807-50. https://doi.org/10.1007/s10654-017-0321-4 PMID: 29064009

39. Walker MG, Volkmuth W, Sprinzak E, Hodgson D, Klingler T. Prediction of gene function by genome-scale expression analysis: prostate cancer-associated genes. Genome Res. 1999; 9(12):1198-203. https://doi.org/10.1101/gr.9.12.1198 PMID: 10613842
40. Sun J, Tao S, Gao Y, Peng T, Tan A, Zhang H, et al. Genome-wide association study identified novel genetic variant on SLC45A3 gene associated with serum levels prostate-specific antigen (PSA) in a Chinese population. Hum Genet. 2013; 132(4):423-9. https://doi.org/10.1007/s00439-012-1254-3 PMID: 23269536

41. Pin E, Henjes F, Hong MG, Wiklund F, Magnusson P, Bjartell A, et al. Identification of a Novel Autoimmune Peptide Epitope of Prostein in Prostate Cancer. J Proteome Res. 2017; 16(1):204-16. https://doi.org/10.1021/acs.jproteome.6b00620 PMID: 27700103

42. Kim S, Shin C, Jee SH. Genetic variants at 1q32.1, 10q11.2 and 19q13.41 are associated with prostate-specific antigen for prostate cancer screening in two Korean population-based cohort studies. Gene. 2015; 556(2):199-205. https://doi.org/10.1016/j.gene.2014.11.059 PMID: 25434496

43. Yanpallewar S, Wang T, Koh DC, Quarta E, Fulgenzi G, Tessarollo L. Nedd4-2 haploinsufficiency causes hyperactivity and increased sensitivity to inflammatory stimuli. Sci Rep. 2016; 6:32957. https://doi.org/10.1038/srep32957 PMID: 27604420

44. Goel P, Manning JA, Kumar S. NEDD4-2 (NEDD4L): the ubiquitin ligase for multiple membrane proteins. Gene. 2015; 557(1):1-10. https://doi.org/10.1016/j.gene.2014.11.051 PMID: 25433090

45. Donovan P, Poronnik P. Nedd4 and Nedd4-2: ubiquitin ligases at work in the neuron. Int J Biochem Cell Biol. 2013; 45(3):706-10. https://doi.org/10.1016/j.biocel.2012.12.006 PMID: 23262292

46. Llewellyn CH, Trzaskowski M, van Jaarsveld CHM, Plomin R, Wardle J. Satiety mechanisms in genetic risk of obesity. JAMA Pediatr. 2014; 168(4):338-44. https://doi.org/10.1001/jamapediatrics.2013.4944 PMID: 24535189

47. Wardle J, Carnell S, Haworth CM, Farooqi IS, O'Rahilly S, Plomin R. Obesity associated genetic variation in FTO is associated with diminished satiety. J Clin Endocrinol Metab. 2008; 93(9):3640-3. https://doi.org/10.1210/jc.2008-0472 PMID: 18583465

48. Magno F, Guarana HC, Fonseca ACP, Cabello GMK, Carneiro JRI, Pedrosa AP, et al. Influence of FTO rs9939609 polymorphism on appetite, ghrelin, leptin, IL6, TNFalpha levels, and food intake of women with morbid obesity. Diabetes Metab Syndr Obes. 2018; 11:199-207. https://doi.org/10.2147/DMSO.S154978 PMID: 29785132

49. Long JZ, Roche AM, Berdan CA, Louie SM, Roberts AJ, Svensson KJ, et al. Ablation of PM20D1 reveals N-acyl amino acid control of metabolism and nociception. Proc Natl Acad Sci U S A. 2018; 115 (29):E6937-E45. https://doi.org/10.1073/pnas.1803389115 PMID: 29967167

50. Long JZ, Svensson KJ, Bateman LA, Lin H, Kamenecka T, Lokurkar IA, et al. The Secreted Enzyme PM20D1 Regulates Lipidated Amino Acid Uncouplers of Mitochondria. Cell. 2016; 166(2):424-35. https://doi.org/10.1016/j.cell.2016.05.071 PMID: 27374330

51. Sovio U, Mook-Kanamori DO, Warrington NM, Lawrence R, Briollais L, Palmer CN, et al. Association between common variation at the FTO locus and changes in body mass index from infancy to late childhood: the complex nature of genetic association through growth and development. PLoS Genet. 2011; 7(2):e1001307. https://doi.org/10.1371/journal.pgen.1001307 PMID: 21379325

52. Baker JL, Olsen LW, Sorensen TI. Childhood body-mass index and the risk of coronary heart disease in adulthood. N Engl J Med. 2007; 357(23):2329-37. https://doi.org/10.1056/NEJMoa072515 PMID: 18057335

53. Tirosh A, Shai I, Afek A, Dubnov-Raz G, Ayalon N, Gordon B, et al. Adolescent BMI trajectory and risk of diabetes versus coronary disease. N Engl J Med. 2011; 364(14):1315-25. https://doi.org/10.1056/NEJMoa1006992 PMID: 21470009

54. Raitakari OT, Juonala M, Kahonen M, Taittonen L, Laitinen T, Maki-Torkko N, et al. Cardiovascular risk factors in childhood and carotid artery intima-media thickness in adulthood: the Cardiovascular Risk in Young Finns Study. Jama. 2003; 290(17):2277-83. https://doi.org/10.1001/jama.290.17.2277 PMID: 14600186

55. Juonala M, Magnussen CG, Berenson GS, Venn A, Burns TL, Sabin MA, et al. Childhood adiposity, adult adiposity, and cardiovascular risk factors. N Engl J Med. 2011; 365(20):1876-85. https://doi.org/10.1056/NEJMoa1010112 PMID: 22087679

56. Burgess JA, Walters EH, Byrnes GB, Giles GG, Jenkins MA, Abramson MJ, et al. Childhood adiposity predicts adult-onset current asthma in females: a 25-yr prospective study. Eur Respir J. 2007; 29 (4):668-75. https://doi.org/10.1183/09031936.00080906 PMID: 17251231

57. Wang YL, Wang YT, Li JF, Zhang YZ, Yin HL, Han B. Body Mass Index and Risk of Parkinson's Disease: A Dose-Response Meta-Analysis of Prospective Studies. PLoS One. 2015; 10(6):e0131778. https://doi.org/10.1371/journal.pone.0131778 PMID: 26121579

58. Umer A, Kelley GA, Cottrell LE, Giacobbi P Jr, Innes KE, Lilly CL. Childhood obesity and adult cardiovascular disease risk factors: a systematic review with meta-analysis. BMC Public Health. 2017; 17 (1):683. https://doi.org/10.1186/s12889-017-4691-z PMID: 28851330
59. Smith GD, Ebrahim S. 'Mendelian randomization': can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003; 32(1):1-22. https://doi.org/10.1093/ije/dyg070 PMID: 12689998

60. Davey Smith G, Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet. 2014; 23(R1):R89-98. https://doi.org/10.1093/hmg/ddu328 PMID: 25064373

61. Currie C, Ahluwalia N, Godeau E, Nic Gabhainn S, Due P, Currie DB. Is obesity at individual and national level associated with lower age at menarche? Evidence from 34 countries in the Health Behaviour in School-aged Children Study. J Adolesc Health. 2012; 50(6):621-6. https://doi.org/10.1016/j.jadohealth.2011.10.254 PMID: 22626490

62. Richardson TG, Sanderson E, Elsworth B, Tilling K, Davey Smith G. Can the impact of childhood adiposity on disease risk be reversed? A Mendelian randomization study. medRxiv. 2019. https://doi.org/10.1101/19008011

63. Roos E, Grotta A, Yang F, Bellocco R, Ye W, Adami HO, et al. Body mass index, sitting time, and risk of Parkinson disease. Neurology. 2018; 90(16):e1413-e7. https://doi.org/10.1212/WNL. 0000000000005328 PMID: 29661903

64. Rokholm B, Silventoinen K, Angquist L, Skytthe A, Kyvik KO, Sorensen TI. Increased genetic variance of BMI with a higher prevalence of obesity. PLoS One. 2011; 6(6):e20816. https://doi.org/10.1371/journal.pone.0020816 PMID: 21738588

65. Rokholm B, Silventoinen K, Tynelius P, Gamborg M, Sorensen TI, Rasmussen F. Increasing genetic variance of body mass index during the Swedish obesity epidemic. PLoS One. 2011; 6(11):e27135. https://doi.org/10.1371/journal.pone.0027135 PMID: 22087252

66. Ajslev TA, Angquist L, Silventoinen K, Baker JL, Sorensen TI. Trends in parent-child correlations of childhood body mass index during the development of the obesity epidemic. PLoS One. 2014; 9(10): e109932. https://doi.org/10.1371/journal.pone.0109932 PMID: 25329656

67. Ajslev TA, Angquist L, Silventoinen K, Baker JL, Sorensen TI. Stable intergenerational associations of childhood overweight during the development of the obesity epidemic. Obesity (Silver Spring). 2015; 23 (6):1279-87.

68. Alves AC, De Silva NMG, Karhunen V, Sovio U, Das S, Taal HR, et al. GWAS on longitudinal growth traits reveal different genetic factors influencing infant, child and adult BMI. Science Advances, 2019 Sep 4;5(9).

69. Helgeland Ø, Vaudel M, Juliusson PB, Holmen OL, Juodakis J, Johansson S, et al. Genome-wide association study reveals a dynamic role of common genetic variation in infant and early childhood growth. Nature Communications, 10, 4448 (2019). https://doi.org/10.1038/s41467-019-12308-0 PMID: 31575865

70. Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010; 26(17):2190-1. https://doi.org/10.1093/bioinformatics/btq340 PMID: 20616382

71. Du P, Kibbe WA, Lin SM. Iumi: a pipeline for processing Illumina microarray. Bioinformatics. 2008; 24 (13):1547-8. https://doi.org/10.1093/bioinformatics/btn224 PMID: 18467348

72. Schmid R, Baum P, Ittrich C, Fundel-Clemens K, Huber W, Brors B, et al. Comparison of normalization methods for Illumina BeadChip HumanHT-12 v3. BMC Genomics. 2010; 11:349. https://doi.org/10.1186/1471-2164-11-349 PMID: 20525181

73. Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007; 8(1):118-27. https://doi.org/10.1093/biostatistics/kxj037 PMID: 16632515

74. Barbosa-Morais NL, Dunning MJ, Samarajiwa SA, Darot JF, Ritchie ME, Lynch AG, et al. A re-annotation pipeline for Illumina BeadArrays: improving the interpretation of gene expression data. Nucleic Acids Res. 2010; 38(3):e17. https://doi.org/10.1093/nar/gkp942 PMID: 19923232

75. Shabalin AA. Matrix eQTL: ultra fast eQTL analysis via large matrix operations. Bioinformatics. 2012; 28(10):1353-8. https://doi.org/10.1093/bioinformatics/bts163 PMID: 22492648