Genome-Wide Association Study Identifies Novel Genetic Variants Contributing to Variation in Blood Metabolite Levels
Total Page:16
File Type:pdf, Size:1020Kb
ARTICLE Received 15 Aug 2014 | Accepted 20 Apr 2015 | Published 12 Jun 2015 DOI: 10.1038/ncomms8208 Genome-wide association study identifies novel genetic variants contributing to variation in blood metabolite levels Harmen H.M. Draisma1,2,3,*, Rene´ Pool1,2,4,*, Michael Kobl5, Rick Jansen6, Ann-Kristin Petersen5, Anika A.M. Vaarhorst4,7,8, Idil Yet9, Toomas Haller10,Ays¸e Demirkan11,12,To˜nu Esko10,13,14,15, Gu Zhu16, Stefan Bo¨hringer17, Marian Beekman7, Jan Bert van Klinken11, Werner Ro¨misch-Margl18, Cornelia Prehn19, Jerzy Adamski19,20,21, Anton J.M. de Craen22, Elisabeth M. van Leeuwen12, Najaf Amin12, Harish Dharuri11, Harm-Jan Westra23, Lude Franke23, Eco J.C. de Geus1,2, Jouke Jan Hottenga1,2, Gonneke Willemsen1,2, Anjali K. Henders16, Grant W. Montgomery16, Dale R. Nyholt16,24, John B. Whitfield16, Brenda W. Penninx2,25, Tim D. Spector9, Andres Metspalu10, P. Eline Slagboom4,7, Ko Willems van Dijk11,26, Peter A.C. ‘t Hoen11, Konstantin Strauch5,27, Nicholas G. Martin16, Gert-Jan B. van Ommen11, Thomas Illig28,29,30, Jordana T. Bell9, Massimo Mangino9, Karsten Suhre18,31,32, Mark I. McCarthy33,34,35, Christian Gieger5,28,36, Aaron Isaacs12, Cornelia M. van Duijn4,8,12,** & Dorret I. Boomsma1,2,4,** Metabolites are small molecules involved in cellular metabolism, which can be detected in biological samples using metabolomic techniques. Here we present the results of genome-wide association and meta-analyses for variation in the blood serum levels of 129 metabolites as measured by the Biocrates metabolomic platform. In a discovery sample of 7,478 individuals of European descent, we find 4,068 genome- and metabolome-wide significant (Z-test, Po1.09 Â 10 À 9) associations between single-nucleo- tide polymorphisms (SNPs) and metabolites, involving 59 independent SNPs and 85 metabolites. Five of the fifty-nine independent SNPs are new for serum metabolite levels, and were followed-up for replication in an independent sample (N ¼ 1,182). The novel SNPs are located in or near genes encoding metabolite transporter proteins or enzymes (SLC22A16, ARG1, AGPS and ACSL1) that have demonstrated biomedical or pharmaceutical importance. The further characterization of genetic influences on metabolic phenotypes is important for progress in biological and medical research. 1 Department of Biological Psychology, VU University Amsterdam, van der Boechorststraat 1, Amsterdam 1081 BT, The Netherlands. 2 The EMGO þ Institute for Health and Care Research, VU University Medical Center, Van der Boechorststraat 7, Amsterdam 1081 BT, The Netherlands. 3 Neuroscience Campus Amsterdam, De Boelelaan 1085, Amsterdam 1081 HV, The Netherlands. 4 BBMRI- NL: Infrastructure for the Application of Metabolomics Technology in Epidemiology (RP4), S4-P, Postbus 9600, Leiden 2300 RC, The Netherlands. 5 Institute of Genetic Epidemiology, Helmholtz Zentrum Mu¨nchen—German Research Center for Environmental Health, Ingolsta¨dter Landstrae 1, Neuherberg 85764, Germany. 6 Department of Psychiatry, VU University Medical Center, Neuroscience Campus Amsterdam, VUmc, A.J. Ernststraat 1187, Amsterdam 1081 HL, The Netherlands. 7 Department of Molecular Epidemiology, Leiden University Medical Center, PO Box 9600, Leiden 2300 RC, The Netherlands. 8 Netherlands Consortium for Healthy Aging, Leiden University Medical Center, Leiden, The Netherlands. 9 Department of Twin Research and Genetic Epidemiology, King’s College London, Westminster Bridge Road, London SE1 7EH, UK. 10 Estonian Genome Center, University of Tartu, 23b Riia Street, Tartu 51010, Estonia. 11 Department of Human Genetics, Leiden University Medical Center, S4-P, PO Box 9600, Leiden 2300 RC, The Netherlands. 12 Genetic Epidemiology Unit, Department of Epidemiology, Erasmus Medical Center, P.O. Box 2040, Rotterdam 3000 CA, The Netherlands. 13 Divisions of Endocrinology and Genetics and Center for Basic and Translational Obesity Research, Boston Children’s Hospital, 300 Longwood Ave, Boston MA02115, Massachusetts, USA. 14 Medical and Population Genetics Program, Broad Institute of MITand Harvard, Cambridge, Massachusetts 2142, USA. 15 Department of Genetics, Harvard Medical School, 77 Avenue Louis Pasteur, NRB 0330, Boston MA 02115, Massachusetts, USA. 16 Department of Genetics and Computational Biology, QIMR Berghofer Medical Research Institute, 300 Herston Road, Brisbane 4006, Australia. 17 Department of Medical Statistics and Bioinformatics, Leiden University Medical Center, PO Box 9600, Leiden 2300 RC, The Netherlands. 18 Institute of Bioinformatics and Systems Biology, Helmholtz Zentrum Mu¨nchen, Ingolsta¨dter Landstrae 1, Neuherberg 85764, Germany. 19 Institute of Experimental Genetics, Genome Analysis Center, Helmholtz Zentrum Mu¨nchen, Ingolsta¨dter Landstrae1, Neuherberg 85764, Germany. 20 German Center for Diabetes Research at Helmholtz Zentrum Mu¨nchen, Ingolsta¨dter Landstr. 1, Neuherberg 85764, Germany. 21 Lehrstuhl fu¨r Experimentelle Genetik, Technische Universita¨t Mu¨nchen, Freising-Weihenstephan 85350, Germany. 22 Department of Gerontology and Geriatrics, Leiden University Medical Center, PO Box 9600, Leiden 2300 RC, The Netherlands. 23 Department of Genetics, CB50, University Medical Center Groningen, University of Groningen, P.O. Box 30001, Groningen 9700 RB, The Netherlands. 24 Statistical and Genomic Epidemiology, Institute of Health and Biomedical Innovation, Queensland University of Technology, 60 Musk Ave, Kelvin Grove QLD 4059, Queensland, Australia. 25 Department of Psychiatry, VU University Medical Center, A.J. Ernststraat 1187, Amsterdam 1081 HL, The Netherlands. 26 Department of Endocrinology, Leiden University Medical Center, S4-P, PO Box 9600, Leiden 2300 RC, The Netherlands. 27 Institute of Medical Informatics, Biometry and Epidemiology, Chair of Genetic Epidemiology, Ludwig-Maximilians-Universita¨t, Munich 81377, Germany. 28 Research Unit of Molecular Epidemiology, Helmholtz Zentrum Mu¨nchen, Ingolsta¨dter Landstrae 1, Neuherberg 85764, Germany. 29 Hannover Unified Biobank, Hannover Medical School, Hannover 30625, Germany. 30 Institute for Human Genetics, Hannover Medical School, Carl-Neuberg-Strasse 1, Hannover 30625, Germany. 31 Faculty of Biology, Ludwig-Maximilians-Universita¨t, Planegg-Martinsried 82152, Germany. 32 Department of Physiology and Biophysics, Weill Cornell Medical College in Qatar (WCMC-Q), PO Box 24144, Education City—Qatar Foundation, Doha, Qatar. 33 Wellcome Trust Centre for Human Genetics, University of Oxford, Roosevelt Drive, Oxford OX3 7BN, UK. 34 Oxford Centre for Diabetes, Endocrinology and Metabolism, University of Oxford, Churchill Hospital, Headington OX3 7LE, UK. 35 Oxford National Institute for Health Research Biomedical Research Centre, The Joint Research Office, Block 60, Churchill Hospital, Old Road, Headington OX3 7LE, UK. 36 Institute of Epidemiology II, Helmholtz Zentrum Mu¨nchen—German Research Center for Environmental Health, Ingolsta¨dter Landstrae 1, Neuherberg 85764, Germany. * These authors contributed equally to this work. ** These authors jointly supervised this work. Correspondence and requests for materials should be addressed to H.H.M.D. (email: [email protected]). NATURE COMMUNICATIONS | 6:7208 | DOI: 10.1038/ncomms8208 | www.nature.com/naturecommunications 1 & 2015 Macmillan Publishers Limited. All rights reserved. ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8208 etabolite levels in human blood reflect the physiological metabolome-wide significance (Z-test, Po1.09 Â 10 À 9), which state of the body, and may differ between individuals reduced to 123 associations involving 59 independent SNPs and Mbecause of variation in genetic makeup and environ- 85 different metabolites. Of these 123 associations (listed in mental exposure1. The study of the genetic contribution to Supplementary Data 2), 4 represented secondary association variation in metabolite levels is an important basis for improved signals according to approximate conditional analysis. Regional aetiological understanding, prevention, diagnosis and treatment association plots, showing the association signals in the regions of complex disorders1,2. Modern high-throughput metabolomics surrounding the lead metabolomic SNPs, are given for all 123 enables the cost-effective measurement of large metabolite panels associations in Supplementary Fig. 3. SNPs representing in blood samples obtained from many individuals. The data independent association signals were aggregated into 31 genomic generated by such metabolomic experiments have been combined loci, which are listed in Supplementary Data 3. Figure 2 depicts all with genotypic data in several recent genome-wide association associations between loci and metabolites as detected in the (GWA) studies2–12. Indeed, the combined investigation of large discovery phase. numbers of genetic variants and large numbers of metabolic traits Five independent SNPs had not been associated with variation is beginning to draw a systems-wide overview of genetic in serum metabolite levels in previous GWA studies (see Table 1). influences on human metabolism11. However, the heritability To further interpret the association of the remaining 54 SNPs estimates from twin and family studies9–11,13 suggest that with serum metabolite concentrations, we compared our findings additional genetic variants influencing variation in serum with those from 11 published GWA studies2–12 for which at metabolite levels remain to be found in GWA studies. least one of the included metabolites overlapped with the In the current study, we set out to further characterize