Search

Background Circulating white blood cell and platelet traits are clinically linked to various disease outcomes and differ across individuals and ancestry groups. Genetic factors play an important role in determining these traits and many loci have been identified. However, most of these findings were identified in populations of European ancestry (EA), with African Americans (AA), Hispanics/Latinos (HL), and other races/ethnicities being severely underrepresented. Results We performed ancestry-combined and ancestry-specific genome-wide association studies (GWAS) for white blood cell and platelet traits in the ancestrally diverse Population Architecture using Genomics and Epidemiology (PAGE) Study, including 16,201 AA, 21,347 HL, and 27,236 EA participants. We identified six novel findings at suggestive significance (P < 5E-8), which need confirmation, and independent signals at six previously established regions at genome-wide significance (P < 2E-9). We confirmed multiple previously reported genome-wide significant variants in the single variant association analysis and multiple genes using PrediXcan. Evaluation of loci reported from a Euro-centric GWAS indicated attenuation of effect estimates in AA and HL compared to EA populations. Conclusions Our results highlighted the potential to identify ancestry-specific and ancestry-agnostic variants in participants with diverse backgrounds and advocate for continued efforts in improving inclusion of racially/ethnically diverse populations in genetic association studies for complex traits.

Background Quantitative red blood cell (RBC) traits are highly polygenic clinically relevant traits, with approximately 500 reported GWAS loci. The majority of RBC trait GWAS have been performed in European- or East Asian-ancestry populations, despite evidence that rare or ancestry-specific variation contributes substantially to RBC trait heritability. Recently developed combined-phenotype methods which leverage genetic trait correlation to improve statistical power have not yet been applied to these traits. Here we leveraged correlation of seven quantitative RBC traits in performing a combined-phenotype analysis in a multi-ethnic study population. Results We used the adaptive sum of powered scores (aSPU) test to assess combined-phenotype associations between ~ 21 million SNPs and seven RBC traits in a multi-ethnic population (maximum n = 67,885 participants; 24% African American, 30% Hispanic/Latino, and 43% European American; 76% female). Thirty-nine loci in our multi-ethnic population contained at least one significant association signal (p 5%) across all ancestral populations. Nineteen additional independent association signals were identified at seven known loci (HFE, KIT, HBS1L/MYB, CITED2/FILNC1, ABO, HBA1/2, and PLIN4/5). For example, the HBA1/2 locus contained 14 conditionally independent association signals, 11 of which were previously unreported and are specific to African and Amerindian ancestries. One variant in this region was common in all ancestries, but exhibited a narrower LD block in African Americans than European Americans or Hispanics/Latinos. GTEx eQTL analysis of all independent lead SNPs yielded 31 significant associations in relevant tissues, over half of which were not at the gene immediately proximal to the lead SNP. Conclusion This work identified seven loci containing multiple independent association signals for RBC traits using a combined-phenotype approach, which may improve discovery in genetically correlated traits. Highly complex genetic architecture at the HBA1/2 locus was only revealed by the inclusion of African Americans and Hispanics/Latinos, underscoring the continued importance of expanding large GWAS to include ancestrally diverse populations.

Background: Thousands of genetic variants have been associated with hematological traits, though target genes remain unknown at most loci. Also, limited analyses have been conducted in African ancestry and Hispanic/Latino populations; hematological trait associated variants more common in these populations have likely been missed. Methods: To derive gene expression prediction models, we used ancestry-stratified datasets from the Multi-Ethnic Study of Atherosclerosis (MESA, including N=229 African American and N=381 Hispanic/Latino participants, monocytes) and the Depression Genes and Networks study (DGN, N = 922 European ancestry participants, whole blood). We then performed a transcriptome-wide association study (TWAS) for platelet count, hemoglobin, hematocrit, and white blood cell count in African (N = 27,955) and Hispanic/Latino (N = 28,324) ancestry participants. Results: Our results revealed 24 suggestive signals (p < 1×10^(-4)) that were conditionally distinct from known GWAS identified variants and successfully replicated these signals in European ancestry subjects from UK Biobank. We found modestly improved correlation of predicted and measured gene expression in an independent African American cohort (the Genetic Epidemiology Network of Arteriopathy (GENOA) study (N=802), lymphoblastoid cell lines) using the larger DGN reference panel; however, some genes were well predicted using MESA but not DGN. Conclusions: These analyses demonstrate the importance of performing TWAS and other genetic analyses across diverse populations and of balancing sample size and ancestry background matching when selecting a TWAS reference panel.

Hematological measures are important intermediate clinical phenotypes for many acute and chronic diseases. Hematological measures are highly heritable, and although genome-wide association studies (GWAS) have identified thousands of loci containing trait-associated variants, the causal genes underlying these associations are often uncertain. To better understand the underlying genetic regulatory mechanisms, we performed a transcriptome-wide association study (TWAS) using PrediXcan to systematically investigate the association between genetically-predicted gene expression and hematological measures in 54,542 individuals of European ancestry from the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort. We found 239 significant gene-trait associations with hematological measures. Among this set of 239 associations, we replicated 71 at p < 0.05 with same direction of effect for the blood cell trait in a meta-analysis of TWAS results consisting of up to 35,900 European ancestry individuals from the Women’s Health Initiative (WHI), the Atherosclerosis Risk in Communities Study (ARIC), and BioMe Biobank. We further attempted to refine this list of candidate genes by performing conditional analyses, adjusting for individual variants previously associated with these hematological measures, and performed further fine-mapping of TWAS loci. To assist with the interpretation of TWAS findings, we designed an R Shiny application to interactively visualize TWAS results, one genomic locus at a time, by integrating our TWAS results with additional genetic data sources (GWAS, TWAS from other gene expression reference panels, conditional analyses, known GWAS variants, etc.). Our results and R Shiny application highlight frequently overlooked challenges with TWAS and illustrate the complexity of TWAS fine-mapping efforts.Author SummaryTranscriptome-wide association studies (TWAS) have shown great promise in furthering our understanding of the genetic regulatory mechanisms underlying complex trait variation. However, interpreting TWAS results can be incredibly complex, especially in large-scale analyses where hundreds of signals appear throughout the genome, with multiple genes often identified in a single chromosomal region. Our research demonstrates this complexity through real data examples from our analysis of hematological traits, and we provide a useful web application to visualize TWAS results in a broadly approachable format. Together, our results and web application illustrate the importance of interpreting TWAS studies in context and highlight the need to carefully examine results in a region-wide context to draw reasonable conclusions and formulate mechanistic hypotheses.

Fear extinction is an adaptive behavioral process critical for organism’s survival, but deficiency in extinction may lead to PTSD. While the amygdala and its neural circuits are critical for fear extinction, the molecular identity and organizational logic of cell types that lie at the core of these circuits remain unclear. Here we report that mice deficient for amygdala-enrichedgastrin-releasing peptidegene (Grp-/-) exhibit enhanced neuronal activity in the basolateral amygdala (BLA) and stronger fear conditioning, as well as deficient extinction in stress-enhanced fear learning (SEFL). rAAV2-retro-based tracing combined with visualization of the GFP knocked in theGrpgene showed that BLA receives GRPergic or conditioned stimulus projections from the indirect auditory thalamus-to-auditory cortex pathway, ventral hippocampus and ventral tegmental area. Transcription of dopamine-related genes was decreased in BLA ofGrp-/-mice following SEFL extinction recall, suggesting that the GRP may mediate fear extinction regulation by dopamine.Impact statementMice deficient for the amygdala-enrichedgastrin-releasing peptidegene are susceptible to stress-enhanced fear, a behavioral protocol with relevance to PTSD, and show a decrease in dopamine-related gene transcription.

Dubowitz syndrome (DubS) is considered a recognizable syndrome characterized by a distinctive facial appearance and deficits in growth and development. There have been over 200 individuals reported with Dubowitz or a "Dubowitz-like" condition, although no single gene has been implicated as responsible for its cause. We have performed exome (ES) or genome sequencing (GS) for 31 individuals clinically diagnosed with DubS. After genome-wide sequencing, rare variant filtering and computational and Mendelian genomic analyses, a presumptive molecular diagnosis was made in 13/27 (48%) families. The molecular diagnoses included biallelic variants in SKIV2L, SLC35C1, BRCA1, NSUN2; de novo variants in ARID1B, ARID1A, CREBBP, POGZ, TAF1, HDAC8, and copy-number variation at1p36.11(ARID1A), 8q22.2(VPS13B), Xp22, and Xq13(HDAC8). Variants of unknown significance in known disease genes, and also in genes of uncertain significance, were observed in 7/27 (26%) additional families. Only one gene, HDAC8, could explain the phenotype in more than one family (N = 2). All but two of the genomic diagnoses were for genes discovered, or for conditions recognized, since the introduction of next-generation sequencing. Overall, the DubS-like clinical phenotype is associated with extensive locus heterogeneity and the molecular diagnoses made are for emerging clinical conditions sharing characteristic features that overlap the DubS phenotype.

Most genome-wide association and fine-mapping studies to date have been conducted in individuals of European descent, and genetic studies of populations of Hispanic/Latino and African ancestry are limited. In addition, these populations have more complex linkage disequilibrium structure. In order to better define the genetic architecture of these understudied populations, we leveraged >100,000 phased sequences available from deep-coverage whole genome sequencing through the multi-ethnic NHLBI Trans-Omics for Precision Medicine (TOPMed) program to impute genotypes into admixed African and Hispanic/Latino samples with genome-wide genotyping array data. We demonstrated that using TOPMed sequencing data as the imputation reference panel improves genotype imputation quality in these populations, which subsequently enhanced gene-mapping power for complex traits. For rare variants with minor allele frequency (MAF) < 0.5%, we observed a 2.3- to 6.1-fold increase in the number of well-imputed variants, with 11–34% improvement in average imputation quality, compared to the state-of-the-art 1000 Genomes Project Phase 3 and Haplotype Reference Consortium reference panels. Impressively, even for extremely rare variants with minor allele count 86%. Subsequent association analyses of TOPMed reference panel-imputed genotype data with hematological traits (hemoglobin (HGB), hematocrit (HCT), and white blood cell count (WBC)) in ~21,600 African-ancestry and ~21,700 Hispanic/Latino individuals identified associations with two rare variants in the HBB gene (rs33930165 with higher WBC [p = 8.8x10-15] in African populations, rs11549407 with lower HGB [p = 1.5x10-12] and HCT [p = 8.8x10-10] in Hispanics/Latinos). By comparison, neither variant would have been genome-wide significant if either 1000 Genomes Project Phase 3 or Haplotype Reference Consortium reference panels had been used for imputation. Our findings highlight the utility of the TOPMed imputation reference panel for identification of novel rare variant associations not previously detected in similarly sized genome-wide studies of under-represented African and Hispanic/Latino populations., Author summary Admixed African and Hispanic/Latino populations remain understudied in genetic studies of complex diseases. These populations have more complex linkage disequilibrium (LD) structure that can impair mapping of variants. Genotype imputation represents an approach to improve genome coverage, especially for rare or ancestry-specific variation; however, these understudied populations also have smaller relevant imputation reference panels. In this study, we leveraged >100,000 phased sequences generated from the multi-ethnic NHLBI TOPMed project for imputation in ~21,600 individuals of African ancestry (AAs) and ~21,700 Hispanics/Latinos. We demonstrated substantially higher imputation quality for low frequency and rare variants in comparison to the 1000 Genomes Project and Haplotype Reference Consortium reference panels. Analysis of quantitative hematological traits led to the discovery of associations with two rare variants in the HBB gene; one of these variants was replicated in an independent sample, and the other is known to cause anemia in the homozygous state. By comparison, the same HBB variants would not have been genome-wide significant using current reference panels due to lower imputation quality. Our findings demonstrate the power of TOPMed whole genome sequencing data for imputation and subsequent association analysis in admixed African and Hispanic/Latino populations.

Background: Thousands of genetic variants have been associated with hematological traits, though target genes remain unknown at most loci. Moreover, limited analyses have been conducted in African ancestry and Hispanic/Latino populations; hematological trait associated variants more common in these populations have likely been missed. Methods: To derive gene expression prediction models, we used ancestry-stratified datasets from the Multi-Ethnic Study of Atherosclerosis (MESA, including n = 229 African American and n = 381 Hispanic/Latino participants, monocytes) and the Depression Genes and Networks study (DGN, n = 922 European ancestry participants, whole blood). We then performed a transcriptome-wide association study (TWAS) for platelet count, hemoglobin, hematocrit, and white blood cell count in African (n = 27,955) and Hispanic/Latino (n = 28,324) ancestry participants. Results: Our results revealed 24 suggestive signals (p < 1 × 10−4) that were conditionally distinct from known GWAS identified variants and successfully replicated these signals in European ancestry subjects from UK Biobank. We found modestly improved correlation of predicted and measured gene expression in an independent African American cohort (the Genetic Epidemiology Network of Arteriopathy (GENOA) study (n = 802), lymphoblastoid cell lines) using the larger DGN reference panel; however, some genes were well predicted using MESA but not DGN. Conclusions: These analyses demonstrate the importance of performing TWAS and other genetic analyses across diverse populations and of balancing sample size and ancestry background matching when selecting a TWAS reference panel.

Summary/AbstractGenome-wide association studies (GWAS) have laid the foundation for investigations into the biology of complex traits, drug development, and clinical guidelines. However, the dominance of European-ancestry populations in GWAS creates a biased view of the role of human variation in disease, and hinders the equitable translation of genetic associations into clinical and public health applications. The Population Architecture using Genomics and Epidemiology (PAGE) study conducted a GWAS of 26 clinical and behavioral phenotypes in 49,839 non-European individuals. Using strategies designed for analysis of multi-ethnic and admixed populations, we confirm 574 GWAS catalog variants across these traits, and find 38 secondary signals in known loci and 27 novel loci. Our data shows strong evidence of effect-size heterogeneity across ancestries for published GWAS associations, substantial benefits for fine-mapping using diverse cohorts, and insights into clinical implications. We strongly advocate for continued, large genome-wide efforts in diverse populations to reduce health disparities.

Physical activity (PA) may modify the genetic effects that give rise to increased risk of obesity. To identify adiposity loci whose effects are modified by PA, we performed genome-wide interaction meta-analyses of BMI and BMI-adjusted waist circumference and waist-hip ratio from up to 200,452 adults of European (n = 180,423) or other ancestry (n = 20,029). We standardized PA by categorizing it into a dichotomous variable where, on average, 23% of participants were categorized as inactive and 77% as physically active. While we replicate the interaction with PA for the strongest known obesity-risk locus in the FTO gene, of which the effect is attenuated by ~30% in physically active individuals compared to inactive individuals, we do not identify additional loci that are sensitive to PA. In additional genome-wide meta-analyses adjusting for PA and interaction with PA, we identify 11 novel adiposity loci, suggesting that accounting for PA or other environmental factors that contribute to variation in adiposity may facilitate gene discovery.

Physical activity (PA) may modify the genetic effects that give rise to increased risk of obesity. To identify adiposity loci whose effects are modified by PA, we performed genome-wide interaction meta-analyses of BMI and BMI-adjusted waist circumference and waist-hip ratio from up to 200,452 adults of European (n = 180,423) or other ancestry (n = 20,029). We standardized PA by categorizing it into a dichotomous variable where, on average, 23% of participants were categorized as inactive and 77% as physically active. While we replicate the interaction with PA for the strongest known obesity-risk locus in the FTO gene, of which the effect is attenuated by ~30% in physically active individuals compared to inactive individuals, we do not identify additional loci that are sensitive to PA. In additional genome-wide meta-analyses adjusting for PA and interaction with PA, we identify 11 novel adiposity loci, suggesting that accounting for PA or other environmental factors that contribute to variation in adiposity may facilitate gene discovery., The views expressed in this manuscript are those of the authors and do not necessarily represent the views of the National Heart, Lung, and Blood Institute; the National Institutes of Health; or the U.S. Department of Health and Human Services. Funding for this study was provided by the Aase and Ejner Danielsens Foundation; Academy of Finland (102318; 104781, 120315, 123885, 129619, 286284, 134309, 126925, 121584, 124282, 129378, 117787, 250207, 258753, 41071, 77299, 124243, 1114194, 24300796); Accare Center for Child and Adolescent Psychiatry; Action on Hearing Loss (G51); Agence Nationale de la Recherche; Agency for Health Care Policy Research (HS06516); Age UK Research into Ageing Fund; Åke Wiberg Foundation; ALF/LUA Research Grant in Gothenburg; ALFEDIAM; ALK-Abello´ A/S (Hørsholm, Denmark); American Heart Association (13POST16500011, 10SDG269004); Ardix Medical; Arthritis Research UK; Association Diabète Risque Vasculaire; AstraZeneca; Australian Associated Brewers; Australian National Health and Medical Research Council (241944, 339462, 389927, 389875, 389891, 389892, 389938, 442915, 442981, 496739, 552485, 552498); Avera Research Institute; Bayer Diagnostics; Becton Dickinson; Biobanking and Biomolecular Resources Research Infrastructure (BBMRI –NL, 184.021.007); Biocentrum Helsinki; Boston Obesity Nutrition Research Center (DK46200); British Heart Foundation (RG/10/12/28456, SP/04/002); Canada Foundation for Innovation; Canadian Institutes of Health Research (FRN-CCT-83028); Cancer Research UK; Cardionics; Center for Medical Systems Biology; Center of Excellence in Complex Disease Genetics and SALVECenter of Excellence in Genomics (EXCEGEN); Chief Scientist Office of the Scottish Government; City of Kuopio; Cohortes Santé TGIR; Contrat de Projets État-Région; Croatian Science Foundation (8875); Danish Agency for Science, Technology and Innovation; Danish Council for Independent Research (DFF–1333-00124, DFF–1331-007308); Danish Diabetes Academy; Danish Medical Research Council; Department of Psychology and Education of the VU University Amsterdam; Diabetes Hilfs- und Forschungsfonds Deutschland; Dutch Brain Foundation; Dutch Ministry of Justice; Emil Aaltonen Foundation; Erasmus Medical Center; Erasmus University; Estonian Government (IUT20-60, IUT24-6); Estonian Ministry of Education and Research (3.2.0304.11-0312); European Commission (230374, 284167, 323195, 692145, FP7 EurHEALTHAgeing-277849, FP7 BBMRI-LPC 313010, nr 602633, HEALTH-F2-2008-201865-GEFOS, HEALTH-F4-2007-201413, FP6 LSHM-CT-2004-005272, FP5 QLG2-CT-2002-01254, FP6 LSHG-CT-2006-01947, FP7 HEALTH-F4-2007-201413, FP7 279143, FP7 201668, FP7 305739, FP6 LSHG-CT-2006-018947, HEALTH-F4-2007-201413, QLG1-CT-2001-01252); European Regional Development Fund; European Science Foundation (EuroSTRESS project FP-006, ESF, EU/QLRT-2001-01254); Faculty of Biology and Medicine of Lausanne; Federal Ministry of Education and Research (01ZZ9603, 01ZZ0103, 01ZZ0403, 03ZIK012, 03IS2061A); Federal State of Mecklenburg - West Pomerania; Fédération Française de Cardiologie; Finnish Cultural Foundation; Finnish Diabetes Association; Finnish Foundation of Cardiovascular Research; Finnish Heart Association; Food Standards Agency; Fondation de France; Fonds Santé; Genetic Association Information Network of the Foundation for the National Institutes of Health; German Diabetes Association; German Federal Ministry of Education and Research (BMBF, 01ER1206, 01ER1507); German Research Council (SFB-1052, SPP 1629 TO 718/2-1); GlaxoSmithKline; Göran Gustafssons Foundation; Göteborg Medical Society; Health and Safety Executive; Heart Foundation of Northern Sweden; Icelandic Heart Association; Icelandic Parliament; Imperial College Healthcare NHS Trust; INSERM, Réseaux en Santé Publique, Interactions entre les déterminants de la santé; Interreg IV Oberrhein Program (A28); Italian Ministry of Economy and Finance; Italian Ministry of Health (ICS110.1/RF97.71); John D and Catherine T MacArthur Foundation; Juho Vainio Foundation; King's College London; Kjell och Märta Beijers Foundation; Kuopio University Hospital; Kuopio, Tampere and Turku University Hospital Medical Funds (X51001); Leiden University Medical Center; Lilly; LMUinnovativ; Lundbeck Foundation; Lundberg Foundation; Medical Research Council of Canada; MEKOS Laboratories (Denmark); Merck Santé; Mid-Atlantic Nutrition Obesity Research Center (P30 DK72488); Ministère de l’Économie, de l’Innovation et des Exportations; Ministry for Health, Welfare and Sports of the Netherlands; Ministry of Cultural Affairs of the Federal State of Mecklenburg-West Pomerania; Ministry of Education and Culture of Finland (627;2004-2011); Ministry of Education, Culture and Science of the Netherlands; MRC Human Genetics Unit; MRC-GlaxoSmithKline Pilot Programme Grant (G0701863); Municipality of Rotterdam; Netherlands Bioinformatics Centre (2008.024); Netherlands Consortium for Healthy Aging (050-060-810); Netherlands Genomics Initiative; Netherlands Organisation for Health Research and Development (904-61-090, 985-10-002, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, Middelgroot-911-09-032, Spinozapremie 56-464-14192); Netherlands Organisation for Health Research and Development (2010/31471/ZONMW); Netherlands Organisation for Scientific Research (10-000-1002, GB-MW 940-38-011, 100-001-004, 60-60600-97-118, 261-98-710, GB-MaGW 480-01-006, GB-MaGW 480-07-001, GB-MaGW 452-04-314, GB-MaGW 452-06-004, 175.010.2003.005, 175.010.2005.011, 481-08-013, 480-05-003, 911-03-012); Neuroscience Campus Amsterdam; NHS Foundation Trust; Novartis Pharmaceuticals; Novo Nordisk; Office National Interprofessionel des Vins; Paavo Nurmi Foundation; Påhlssons Foundation; Päivikki and Sakari Sohlberg Foundation; Pierre Fabre; Republic of Croatia Ministry of Science, Education and Sport (108-1080315-0302); Research Centre for Prevention and Health, the Capital Region for Denmark; Research Institute for Diseases in the Elderly (014-93-015, RIDE2); Roche; Russian Foundation for Basic Research (NWO-RFBR 047.017.043); Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06); Sanofi-Aventis; Scottish Executive Health Department (CZD/16/6); Siemens Healthcare; Social Insurance Institution of Finland (4/26/2010); Social Ministry of the Federal State of Mecklenburg-West Pomerania; Société Francophone du Diabète; State of Bavaria; Stroke Association; Swedish Diabetes Association; Swedish Foundation for Strategic Research; Swedish Heart-Lung Foundation (20140543); Swedish Research Council (2015-03657); Swedish Medical Research Council (K2007-66X-20270-01-3, 2011-2354); Swedish Society for Medical Research; Swiss National Science Foundation (33CSCO-122661, 33CS30-139468, 33CS30-148401); Tampere Tuberculosis Foundation; The Marcus Borgström Foundation; The Royal Society; The Wellcome Trust (084723/Z/08/Z, 088869/B/09/Z); Timber Merchant Vilhelm Bangs Foundation; Topcon; Torsten and Ragnar Söderberg's Foundation; UK Department of Health; UK Diabetes Association; UK Medical Research Council (MC_U106179471, G0500539, G0600705, G0601966, G0700931, G1002319, K013351, MC_UU_12019/1); UK National Institute for Health Research BioResource Clinical Research Facility and Biomedical Research Centre; UK National Institute for Health Research (NIHR) Comprehensive Biomedical Research Centre; UK National Institute for Health Research (RP-PG-0407-10371); Umeå University Career Development Award; United States – Israel Binational Science Foundation Grant (2011036); University Hospital Oulu (75617); University Medical Center Groningen; University of Tartu (SP1GVARENG); National Institutes of Health (AG13196, CA047988, HHSN268201100046C, HHSN268201100001C, HHSN268201100002C, HHSN268201100003C, HHSN268201100004C, HHSC271201100004C, HHSN268200900041C, HHSN268201300025C, HHSN268201300026C, HHSN268201300027C, HHSN268201300028C, HHSN268201300029C, HHSN268201500001I, HL36310, HG002651, HL034594, HL054457, HL054481, HL071981, HL084729, HL119443, HL126024, N01-AG12100, N01-AG12109, N01-HC25195, N01-HC55015, N01-HC55016, N01-HC55018, N01-HC55019, N01-HC55020, N01-HC55021, N01-HC55022, N01-HD95159, N01-HD95160, N01-HD95161, N01-HD95162, N01-HD95163, N01-HD95164, N01-HD95165, N01-HD95166, N01-HD95167, N01-HD95168, N01-HD95169, N01-HG65403, N02-HL64278, R01-HD057194, R01-HL087641, R01-HL59367, R01HL-086694, R01-HL088451, R24-HD050924, U01-HG-004402, HHSN268200625226C, UL1-RR025005, UL1-RR025005, UL1-TR-001079, UL1-TR-00040, AA07535, AA10248, AA11998, AA13320, AA13321, AA13326, AA14041, AA17688, DA12854, MH081802, MH66206, R01-D004215701A, R01-DK075787, R01-DK089256, R01-DK8925601, R01-HL088451, R01-HL117078, R01-DK062370, R01-DK072193, DK091718, DK100383, DK078616, 1Z01-HG000024, HL087660, HL100245, R01DK089256, 2T32HL007055-36, U01-HL072515-06, U01-HL84756, NIA-U01AG009740, RC2-AG036495, RC4-AG039029, R03 AG046389, 263-MA-410953, 263-MD-9164, 263-MD-821336, U01-HG004802, R37CA54281, R01CA63, P01CA33619, U01-CA136792, U01-CA98758, RC2-MH089951, MH085520, R01-D0042157-01A, MH081802, 1RC2-MH089951, 1RC2-MH089995, 1RL1MH08326801, U01-HG007376, 5R01-HL08767902, 5R01MH63706:02, HG004790, N01-WH22110, U01-HG007033, UM1CA182913, 24152, 32100-2, 32105-6, 32108-9, 32111-13, 32115, 32118-32119, 32122, 42107-26, 42129-32, 44221); USDA National Institute of Food and Agriculture (2007-35205-17883); Västra Götaland Foundation; Velux Foundation; Veterans Affairs (1 IK2 BX001823); Vleugels Foundation; VU University’s Institute for Health and Care Research (EMGO+, HEALTH-F4-2007-201413) and Neuroscience Campus Amsterdam; Wellcome Trust (090532, 091551, 098051, 098381); Wissenschaftsoffensive TMO; and Yrjö Jahnsson Foundation. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Few genome-wide association studies (GWAS) account for environmental exposures, like smoking, potentially impacting the overall trait variance when investigating the genetic contribution to obesity-related traits. Here, we use GWAS data from 51,080 current smokers and 190,178 nonsmokers (87% European descent) to identify loci influencing BMI and central adiposity, measured as waist circumference and waist-to-hip ratio both adjusted for BMI. We identify 23 novel genetic loci, and 9 loci with convincing evidence of gene-smoking interaction (GxSMK) on obesity-related traits. We show consistent direction of effect for all identified loci and significance for 18 novel and for 5 interaction loci in an independent study sample. These loci highlight novel biological functions, including response to oxidative stress, addictive behaviour, and regulatory functions emphasizing the importance of accounting for environment in genetic analyses. Our results suggest that tobacco smoking may alter the genetic susceptibility to overall adiposity and body fat distribution., A full list of acknowledgments appears in the Supplementary Note 4. Co-author A.J.M.d.C. recently passed away while this work was in process. This work was performed under the auspices of the Genetic Investigation of ANthropometric Traits (GIANT) consortium. We acknowledge the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium for encouraging CHARGE studies to participate in this effort and for the contributions of CHARGE members to the analyses conducted for this research. Funding for this study was provided by the Aase and Ejner Danielsens Foundation; Academy of Finland (41071, 77299, 102318, 110413, 117787, 121584, 123885, 124243, 124282, 126925, 129378, 134309, 286284); Accare Center for Child and Adolescent Psychiatry; Action on Hearing Loss (G51); Agence Nationale de la 359 Recherche; Agency for Health Care Policy Research (HS06516); ALF/LUA research grant in Gothenburg; ALFEDIAM; ALK-Abelló A/S; Althingi; American Heart Association (13POST16500011); Amgen; Andrea and Charles Bronfman Philanthropies; Ardix Medical; Arthritis Research UK; Association Diabète Risque Vasculaire; Australian National Health and Medical Research Council (241944, 339462, 389875, 389891, 389892, 389927, 389938, 442915, 442981, 496739, 552485, 552498); Avera Institute; Bayer Diagnostics; Becton Dickinson; BHF (RG/14/5/30893); Boston Obesity Nutrition Research Center (DK46200), Bristol-Myers Squibb; British Heart Foundation (RG/10/12/28456, RG2008/08, RG2008/014, SP/04/002); Medical Research Council of Canada; Canadian Institutes for Health Research (FRCN-CCT-83028); Cancer Research UK; Cardionics; Cavadis B.V., Center for Medical Systems Biology; Center of Excellence in Genomics; CFI; CIHR; City of Kuopio; CNAMTS; Cohortes Santé TGIR; Contrat de Projets État-Région; Croatian Science Foundation (8875); Danish Agency for Science, Technology and Innovation; Danish Council for Independent Research (DFF-1333-00124, DFF-1331-00730B); County Council of Dalarna; Dalarna University; Danish Council for Strategic Research; Danish Diabetes Academy; Danish Medical Research Council; Department of Health, UK; Development Fund from the University of Tartu (SP1GVARENG); Diabetes Hilfs- und Forschungsfonds Deutschland; Diabetes UK; Diabetes Research and Wellness Foundation Fellowship; Donald W. Reynolds Foundation; Dr Robert Pfleger-Stiftung; Dutch Brain Foundation; Dutch Diabetes Research Foundation; Dutch Inter University Cardiology Institute; Dutch Kidney Foundation (E033); Dutch Ministry of Justice; the DynaHEALTH action No. 633595, Economic Structure Enhancing Fund of the Dutch Government; Else Kröner-Fresenius-Stiftung (2012_A147, P48/08//A11/08); Emil Aaltonen Foundation; Erasmus University Medical Center Rotterdam; Erasmus MC and Erasmus University Rotterdam; the Municipality of Rotterdam; Estonian Government (IUT20-60, IUT24-6); Estonian Research Roadmap through the Estonian Ministry of Education and Research (3.2.0304.11-0312); European Research Council (ERC Starting Grant and 323195:SZ-245 50371-GLUCOSEGENES-FP7-IDEAS-ERC); European Regional Development Fund; European Science Foundation (EU/QLRT-2001-01254); European Commission (018947, 018996, 201668, 223004, 230374, 279143, 284167, 305739, BBMRI-LPC-313010, HEALTH-2011.2.4.2-2-EU-MASCARA, HEALTH-2011-278913, HEALTH-2011-294713-EPLORE, HEALTH-F2-2008-201865-GEFOS, HEALTH-F2-2013-601456, HEALTH-F4-2007-201413, HEALTH-F4-2007-201550-HYPERGENES, HEALTH-F7-305507 HOMAGE, IMI/115006, LSHG-CT-2006-018947, LSHG-CT-2006-01947, LSHM-CT-2004-005272, LSHM-CT-2006-037697, LSHM-CT-2007-037273, QLG1-CT-2002-00896, QLG2-CT-2002-01254); Faculty of Biology and Medicine of Lausanne; Federal Ministry of Education and Research (01ZZ0103, 01ZZ0403, 01ZZ9603, 03IS2061A, 03ZIK012); Federal State of Mecklenburg-West Pomerania; Fédération Française de Cardiologie; Finnish Cultural Foundation; Finnish Diabetes Association; Finnish Foundation of Cardiovascular Research; Finnish Heart Association; Fondation Leducq; Food Standards Agency; Foundation for Strategic Research; French Ministry of Research; FRSQ; Genetic Association Information Network (GAIN) of the Foundation for the NIH; German Federal Ministry of Education and Research (BMBF, 01ER1206, 01ER1507); GlaxoSmithKline; Greek General Secretary of Research and Technology; Göteborg Medical Society; Health and Safety Executive; Healthcare NHS Trust; Healthway; Western Australia; Heart Foundation of Northern Sweden; Helmholtz Zentrum München—German Research Center for Environmental Health; Hjartavernd; Ingrid Thurings Foundation; INSERM; InterOmics (PB05 MIUR-CNR); INTERREG IV Oberrhein Program (A28); Interuniversity Cardiology Institute of the Netherlands (ICIN, 09.001); Italian Ministry of Health (ICS110.1/RF97.71); Italian Ministry of Economy and Finance (FaReBio di Qualità); Marianne and Marcus Wallenberg Foundation; the Ministry of Health, Welfare and Sports, the Netherlands; J.D.E. and Catherine T, MacArthur Foundation Research Networks on Successful Midlife Development and Socioeconomic Status and Health; Juho Vainio Foundation; Juvenile Diabetes Research Foundation International; KfH Stiftung Präventivmedizin e.V.; King's College London; Knut and Alice Wallenberg Foundation; Kuopio University Hospital; Kuopio, Tampere and Turku University Hospital Medical Funds (X51001); La Fondation de France; Leenaards Foundation; Lilly; LMUinnovativ; Lundberg Foundation; Magnus Bergvall Foundation; MDEIE; Medical Research Council UK (G0000934, G0601966, G0700931, MC_U106179471, MC_UU_12019/1); MEKOS Laboratories; Merck Santé; Ministry for Health, Welfare and Sports, The Netherlands; Ministry of Cultural Affairs of Mecklenburg-West Pomerania; Ministry of Economic Affairs, The Netherlands; Ministry of Education and Culture of Finland (627;2004-2011); Ministry of Education, Culture and Science, The Netherlands; Ministry of Science, Education and Sport in the Republic of Croatia (108-1080315-0302); MRC centre for Causal Analyses in Translational Epidemiology; MRC Human Genetics Unit; MRC-GlaxoSmithKline pilot programme (G0701863); MSD Stipend Diabetes; National Institute for Health Research; Netherlands Brain Foundation (F2013(1)-28); Netherlands CardioVascular Research Initiative (CVON2011-19); Netherlands Genomics Initiative (050-060-810); Netherlands Heart Foundation (2001 D 032, NHS2010B280); Netherlands Organization for Scientific Research (NWO) and Netherlands Organisation for Health Research and Development (ZonMW) (56-464-14192, 60-60600-97-118, 100-001-004, 261-98-710, 400-05-717, 480-04-004, 480-05-003, 481-08-013, 904-61-090, 904-61-193, 911-11-025, 985-10-002, Addiction-31160008, BBMRI–NL 184.021.007, GB-MaGW 452-04-314, GB-MaGW 452-06-004, GB-MaGW 480-01-006, GB-MaGW 480-07-001, GB-MW 940-38-011, Middelgroot-911-09-032, NBIC/BioAssist/RK 2008.024, Spinozapremie 175.010.2003.005, 175.010.2007.006); Neuroscience Campus Amsterdam; NHS Foundation Trust; National Institutes of Health (1RC2MH089951, 1Z01HG000024, 24152, 263MD9164, 263MD821336, 2R01LM010098, 32100-2, 32122, 32108, 5K99HL130580-02, AA07535, AA10248, AA11998, AA13320, AA13321, AA13326, AA14041, AA17688, AG13196, CA047988, DA12854, DK56350, DK063491, DK078150, DK091718, DK100383, DK078616, ES10126, HG004790, HHSN268200625226C, HHSN268200800007C, HHSN268201200036C, HHSN268201500001I, HHSN268201100046C, HHSN268201100001C, HHSN268201100002C, HHSN268201100003C, HHSN268201100004C, HHSN271201100004C, HL043851, HL45670, HL080467, HL085144, HL087660, HL054457, HL119443, HL118305, HL071981, HL034594, HL126024, HL130114, KL2TR001109, MH66206, MH081802, N01AG12100, N01HC55015, N01HC55016, N01C55018, N01HC55019, N01HC55020, N01HC55021, N01HC55022, N01HC85079, N01HC85080, N01HC85081, N01HC85082, N01HC85083, N01HC85086, N01HC95159, N01HC95160, N01HC95161, N01HC95162, N01HC95163, N01HC95164, N01HC95165, N01HC95166, N01HC95167, N01HC95168, N01HC95169, N01HG65403, N01WH22110, N02HL6‐4278, N01-HC-25195, P01CA33619, R01HD057194, R01HD057194, R01AG023629, R01CA63, R01D004215701A, R01DK075787, R01DK062370, R01DK072193, R01DK075787, R01DK089256, R01HL53353, R01HL59367, R01HL086694, R01HL087641, R01HL087652, R01HL103612, R01HL105756, R01HL117078, R01HL120393, R03 AG046389, R37CA54281, RC2AG036495, RC4AG039029, RPPG040710371, RR20649, TW008288, TW05596, U01AG009740, U01CA98758, U01CA136792, U01DK062418, U01HG004402, U01HG004802, U01HG007376, U01HL080295, UL1RR025005, UL1TR000040, UL1TR000124, UL1TR001079, 2T32HL007055-36, T32GM074905, HG002651, HL084729, N01-HC-25195, UM1CA182913); NIH, National Institute on Aging (Intramural funding, NO1-AG-1-2109); Northern Netherlands Collaboration of Provinces; Novartis Pharma; Novo Nordisk; Novo Nordisk Foundation; Nutricia Research Foundation (2016-T1); ONIVINS; Parnassia Bavo group; Pierre Fabre; Province of Groningen; Päivikki and Sakari Sohlberg Foundation; Påhlssons Foundation; Paavo Nurmi Foundation; Radboud Medical Center Nijmegen; Research Centre for Prevention and Health, the Capital Region of Denmark; the Research Institute for Diseases in the Elderly; Research into Ageing; Robert Dawson Evans Endowment of the Department of Medicine at Boston University School of Medicine and Boston Medical Center; Roche; Royal Society; Russian Foundation for Basic Research (NWO-RFBR 047.017.043); Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06); Sanofi-Aventis; Scottish Government Health Directorates, Chief Scientist Office (CZD/16/6); Siemens Healthcare; Social Insurance Institution of Finland (4/26/2010); Social Ministry of the Federal State of Mecklenburg-West Pomerania; Société Francophone du 358 Diabète; State of Bavaria; Stiftelsen för Gamla Tjänarinnor; Stockholm County Council (560183, 592229); Strategic Cardiovascular and Diabetes Programmes of Karolinska Institutet and Stockholm County Council; Stroke Association; Swedish Diabetes Association; Swedish Diabetes Foundation (2013-024); Swedish Foundation for Strategic Research; Swedish Heart-Lung Foundation (20120197, 20150711); Swedish Research Council (0593, 8691, 2012-1397, 2012-1727, and 2012-2215); Swedish Society for Medical Research; Swiss Institute of Bioinformatics; Swiss National Science Foundation (3100AO-116323/1, 31003A-143914, 33CSCO-122661, 33CS30-139468, 33CS30-148401, 51RTP0_151019); Tampere Tuberculosis Foundation; Technology Foundation STW (11679); The Fonds voor Wetenschappelijk Onderzoek Vlaanderen, Ministry of the Flemish Community (G.0880.13, G.0881.13); The Great Wine Estates of the Margaret River Region of Western Australia; Timber Merchant Vilhelm Bangs Foundation; Topcon; Tore Nilsson Foundation; Torsten and Ragnar Söderberg's Foundation; United States – Israel Binational Science Foundation (Grant 2011036), Umeå University; University Hospital of Regensburg; University of Groningen; University Medical Center Groningen; University of Michigan; University of Utrecht; Uppsala Multidisciplinary Center for Advanced Computational Science (UPPMAX) (b2011036); Velux Foundation; VU University’s Institute for Health and Care Research; Västra Götaland Foundation; Wellcome Trust (068545, 076113, 079895, 084723, 088869, WT064890, WT086596, WT098017, WT090532, WT098051, 098381); Wissenschaftsoffensive TMO; Yrjö Jahnsson Foundation; and Åke Wiberg Foundation. The views expressed in this manuscript are those of the authors and do not necessarily represent the views of the National Heart, Lung, and Blood Institute (NHLBI); the National Institutes of Health (NIH); or the U.S. Department of Health and Human Services.

Using a phenome-wide association study (PheWAS) approach, we comprehensively tested genetic variants for association with phenotypes available for 70,061 study participants in the Population Architecture using Genomics and Epidemiology (PAGE) network. Our aim was to better characterize the genetic architecture of complex traits and identify novel pleiotropic relationships. This PheWAS drew on five population-based studies representing four major racial/ethnic groups (European Americans (EA), African Americans (AA), Hispanics/Mexican-Americans, and Asian/Pacific Islanders) in PAGE, each site with measurements for multiple traits, associated laboratory measures, and intermediate biomarkers. A total of 83 single nucleotide polymorphisms (SNPs) identified by genome-wide association studies (GWAS) were genotyped across two or more PAGE study sites. Comprehensive tests of association, stratified by race/ethnicity, were performed, encompassing 4,706 phenotypes mapped to 105 phenotype-classes, and association results were compared across study sites. A total of 111 PheWAS results had significant associations for two or more PAGE study sites with consistent direction of effect with a significance threshold of p, Author Summary In phenome-wide association studies (PheWAS) all potential genetic variants in a dataset are systematically tested for association with all available phenotypes and traits that have been measured in study participants. By investigating the relationship between genetic variation and a diversity of phenotypes, there is the potential for uncovering novel relationships between single nucleotide polymorphisms (SNPs), phenotypes, and networks of interrelated phenotypes. PheWAS also can expose pleiotropy, provide novel mechanistic insights, and foster hypothesis generation. This approach is complementary to genome-wide association studies (GWAS) that test the association between hundreds of thousands, to over a million, single nucleotide polymorphisms and a single phenotype or limited phenotypic domain. The Population Architecture using Genomics and Epidemiology (PAGE) network has measures for a wide array of phenotypes and traits, including prevalent and incident status for clinical conditions and risk factors, as well as clinical parameters and intermediate biomarkers. We performed tests of association between a series of genome-wide association study (GWAS)–identified SNPs and a comprehensive range of phenotypes from the PAGE network in a high-throughput manner. We replicated a number of previously reported associations, validating the PheWAS approach. We also identified novel genotype–phenotype associations possibly representing pleiotropic effects.