Genetic diversity fuels gene discovery for tobacco and alcohol use.

Tobacco and alcohol use are heritable behaviours associated with 15% and 5.3% of worldwide deaths, respectively, due largely to broad increased risk for disease and injury1-4. These substances are used across the globe, yet genome-wide association studies have focused largely on individuals of European ancestries5. Here we leveraged global genetic diversity across 3.4 million individuals from four major clines of global ancestry (approximately 21% non-European) to power the discovery and fine-mapping of genomic loci associated with tobacco and alcohol use, to inform function of these loci via ancestry-aware transcriptome-wide association studies, and to evaluate the genetic architecture and predictive power of polygenic risk within and across populations. We found that increases in sample size and genetic diversity improved locus identification and fine-mapping resolution, and that a large majority of the 3,823 associated variants (from 2,143 loci) showed consistent effect sizes across ancestry dimensions. However, polygenic risk scores developed in one ancestry performed poorly in others, highlighting the continued need to increase sample sizes of diverse ancestries to realize any potential benefit of polygenic prediction.

Sequencing of 53,831 diverse genomes from the NHLBI TOPMed Program.

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

Whole-genome analysis of plasma fibrinogen reveals population-differentiated genetic regulators with putative liver roles.

Genetic studies have identified numerous regions associated with plasma fibrinogen levels in Europeans, yet missing heritability and limited inclusion of non-Europeans necessitates further studies with improved power and sensitivity. Compared with array-based genotyping, whole genome sequencing (WGS) data provides better coverage of the genome and better representation of non-European variants. To better understand the genetic landscape regulating plasma fibrinogen levels, we meta-analyzed WGS data from the NHLBI's Trans-Omics for Precision Medicine (TOPMed) program (n=32,572), with array-based genotype data from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium (n=131,340) imputed to the TOPMed or Haplotype Reference Consortium panel. We identified 18 loci that have not been identified in prior genetic studies of fibrinogen. Of these, four are driven by common variants of small effect with reported MAF at least 10 percentage points higher in African populations. Three signals (SERPINA1, ZFP36L2, and TLR10) contain predicted deleterious missense variants. Two loci, SOCS3 and HPN, each harbor two conditionally distinct, non-coding variants. The gene region encoding the fibrinogen protein chain subunits (FGG;FGB;FGA), contains 7 distinct signals, including one novel signal driven by rs28577061, a variant common in African ancestry populations but extremely rare in Europeans (MAFAFR=0.180; MAFEUR=0.008). Through phenome-wide association studies in the VA Million Veteran Program, we found associations between fibrinogen polygenic risk scores and thrombotic and inflammatory disease phenotypes, including an association with gout. Our findings demonstrate the utility of WGS to augment genetic discovery in diverse populations and offer new insights for putative mechanisms of fibrinogen regulation.

Genome-wide association meta-analysis identifies 48 risk variants and highlights the role of the stria vascularis in hearing loss.

Hearing loss is one of the top contributors to years lived with disability and is a risk factor for dementia. Molecular evidence on the cellular origins of hearing loss in humans is growing. Here, we performed a genome-wide association meta-analysis of clinically diagnosed and self-reported hearing impairment on 723,266 individuals and identified 48 significant loci, 10 of which are novel. A large proportion of associations comprised missense variants, half of which lie within known familial hearing loss loci. We used single-cell RNA-sequencing data from mouse cochlea and brain and mapped common-variant genomic results to spindle, root, and basal cells from the stria vascularis, a structure in the cochlea necessary for normal hearing. Our findings indicate the importance of the stria vascularis in the mechanism of hearing impairment, providing future paths for developing targets for therapeutic intervention in hearing loss.

Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity.

STUDY QUESTION: Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins? SUMMARY ANSWER: We identified four new loci, GNRH1, FSHR, ZFPM1, and IPO8, in addition to previously identified loci, FSHB and SMAD3. WHAT IS KNOWN ALREADY: The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures. STUDY DESIGN, SIZE, DURATION: We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls). PARTICIPANTS/MATERIALS, SETTING, METHODS: Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation. MAIN RESULTS AND THE ROLE OF CHANCE: This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR, have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle. LARGE SCALE DATA: The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/). LIMITATIONS, REASONS FOR CAUTION: Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility. WIDER IMPLICATIONS OF THE FINDINGS: About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility. STUDY FUNDING/COMPETING INTEREST(S): Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE-European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Sklodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests. TRIAL REGISTRATION NUMBER: N/A.

Broadcasters, receivers, functional groups of metabolites, and the link to heart failure by revealing metabolomic network connectivity.

BACKGROUND AND OBJECTIVE: Blood-based small molecule metabolites offer easy accessibility and hold significant potential for insights into health processes, the impact of lifestyle, and genetic variation on disease, enabling precise risk prevention. In a prospective study with records of heart failure (HF) incidence, we present metabolite profiling data from individuals without HF at baseline. METHODS: We uncovered the interconnectivity of metabolites using data-driven and causal networks augmented with polygenic factors. Exploring the networks, we identified metabolite broadcasters, receivers, mediators, and subnetworks corresponding to functional classes of metabolites, and provided insights into the link between metabolomic architecture and regulation in health. We incorporated the network structure into the identification of metabolites associated with HF to control the effect of confounding metabolites. RESULTS: We identified metabolites associated with higher and lower risk of HF incidence, such as glycine, ureidopropionic and glycocholic acids, and LPC 18:2. These associations were not confounded by the other metabolites due to uncovering the connectivity among metabolites and adjusting each association for the confounding metabolites. Examples of our findings include the direct influence of asparagine on glycine, both of which were inversely associated with HF. These two metabolites were influenced by polygenic factors and only essential amino acids, which are not synthesized in the human body and are obtained directly from the diet. CONCLUSION: Metabolites may play a critical role in linking genetic background and lifestyle factors to HF incidence. Revealing the underlying connectivity of metabolites associated with HF strengthens the findings and facilitates studying complex conditions like HF.

Association of Physical Activity With Bioactive Lipids and Cardiovascular Events.

BACKGROUND: To clarify the mechanisms underlying physical activity (PA)-related cardioprotection, we examined the association of PA with plasma bioactive lipids (BALs) and cardiovascular disease (CVD) events. We additionally performed genome-wide associations. METHODS: PA-bioactive lipid associations were examined in VITAL (VITamin D and OmegA-3 TriaL)-clinical translational science center (REGISTRATION: URL: https://www. CLINICALTRIALS: gov; Unique identifier: NCT01169259; N=1032) and validated in JUPITER (Justification for the Use of statins in Prevention: an Intervention Trial Evaluating Rosuvastatin)-NC (REGISTRATION: URL: https://www. CLINICALTRIALS: gov; Unique identifier: NCT00239681; N=589), using linear models adjusted for age, sex, race, low-density lipoprotein-cholesterol, total-C, and smoking. Significant BALs were carried over to examine associations with incident CVD in 2 nested CVD case-control studies: VITAL-CVD (741 case-control pairs) and JUPITER-CVD (415 case-control pairs; validation). RESULTS: We detected 145 PA-bioactive lipid validated associations (false discovery rate <0.1). Annotations were found for 6 of these BALs: 12,13-diHOME, 9,10-diHOME, lysoPC(15:0), oxymorphone-3b-D-glucuronide, cortisone, and oleoyl-glycerol. Genetic analysis within JUPITER-NC showed associations of 32 PA-related BALs with 22 single-nucleotide polymorphisms. From PA-related BALs, 12 are associated with CVD. CONCLUSIONS: We identified a PA-related bioactive lipidome profile out of which 12 BALs also had opposite associations with incident CVD events.

Cross-Ancestry Investigation of Venous Thromboembolism Genomic Predictors.

BACKGROUND: Venous thromboembolism (VTE) is a life-threatening vascular event with environmental and genetic determinants. Recent VTE genome-wide association studies (GWAS) meta-analyses involved nearly 30 000 VTE cases and identified up to 40 genetic loci associated with VTE risk, including loci not previously suspected to play a role in hemostasis. The aim of our research was to expand discovery of new genetic loci associated with VTE by using cross-ancestry genomic resources. METHODS: We present new cross-ancestry meta-analyzed GWAS results involving up to 81 669 VTE cases from 30 studies, with replication of novel loci in independent populations and loci characterization through in silico genomic interrogations. RESULTS: In our genetic discovery effort that included 55 330 participants with VTE (47 822 European, 6320 African, and 1188 Hispanic ancestry), we identified 48 novel associations, of which 34 were replicated after correction for multiple testing. In our combined discovery-replication analysis (81 669 VTE participants) and ancestry-stratified meta-analyses (European, African, and Hispanic), we identified another 44 novel associations, which are new candidate VTE-associated loci requiring replication. In total, across all GWAS meta-analyses, we identified 135 independent genomic loci significantly associated with VTE risk. A genetic risk score of the significantly associated loci in Europeans identified a 6-fold increase in risk for those in the top 1% of scores compared with those with average scores. We also identified 31 novel transcript associations in transcriptome-wide association studies and 8 novel candidate genes with protein quantitative-trait locus Mendelian randomization analyses. In silico interrogations of hemostasis and hematology traits and a large phenome-wide association analysis of the 135 GWAS loci provided insights to biological pathways contributing to VTE, with some loci contributing to VTE through well-characterized coagulation pathways and others providing new data on the role of hematology traits, particularly platelet function. Many of the replicated loci are outside of known or currently hypothesized pathways to thrombosis. CONCLUSIONS: Our cross-ancestry GWAS meta-analyses identified new loci associated with VTE. These findings highlight new pathways to thrombosis and provide novel molecules that may be useful in the development of improved antithrombosis treatments.

Genome-wide gene-diet interaction analysis in the UK Biobank identifies novel effects on hemoglobin A1c.

Diet is a significant modifiable risk factor for type 2 diabetes (T2D), and its effect on disease risk is under partial genetic control. Identification of specific gene-diet interactions (GDIs) influencing risk biomarkers such as glycated hemoglobin (HbA1c) is a critical step towards precision nutrition for T2D prevention, but progress has been slow due to limitations in sample size and accuracy of dietary exposure measurement. We leveraged the large UK Biobank (UKB) cohort and a diverse group of dietary exposures, including 30 individual dietary traits and 8 empirical dietary patterns, to conduct genome-wide interaction studies in ~340 000 European-ancestry participants to identify novel GDIs influencing HbA1c. We identified five variant-dietary trait pairs reaching genome-wide significance (P < 5 × 10-8): two involved dietary patterns (meat pattern with rs147678157 and a fruit & vegetable-based pattern with rs3010439) and three involved individual dietary traits (bread consumption with rs62218803, dried fruit consumption with rs140270534 and milk type [dairy vs. other] with 4:131148078_TAGAA_T). These were affected minimally by adjustment for geographical and lifestyle-related confounders, and four of the five variants lacked genetic main effects that would have allowed their detection in a traditional genome-wide association study for HbA1c. Notably, multiple loci near transient receptor potential subfamily M genes (TRPM2 and TRPM3) interacted with carbohydrate-containing food groups. These interactions were further characterized using non-European UKB subsets and alternative measures of glycaemia (fasting glucose and follow-up HbA1c measurements). Our results highlight GDIs influencing HbA1c for future investigation, while reinforcing known challenges in detecting and replicating GDIs.

Pleiotropy-Based Decomposition of Genetic Risk Scores: Association and Interaction Analysis for Type 2 Diabetes and CAD.

Genetic risk for a disease in the population may be represented as a genetic risk score (GRS) constructed as the sum of inherited risk alleles, weighted by allelic effects established in an independent population. While this formulation captures overall genetic risk, it typically does not address risk due to specific biological mechanisms or pathways that may nevertheless be important for interpretation or treatment response. Here, a GRS for disease is resolved into independent or nearly independent components pertaining to biological mechanisms inferred from pleiotropic relationships. The component GRSs' weights are derived from the singular value decomposition (SVD) of the matrix of appropriately scaled genetic effects, i.e., beta coefficients, of the disease variants across a panel of the disease-related phenotypes. The SVD-based formalism also associates combinations of disease-related phenotypes with inferred disease pathways. Applied to incident type 2 diabetes (T2D) in the Women's Genome Health Study (N = 23,294), component GRSs discriminate glycemic control and lipid-based genetic risk, while revealing significant interactions between specific components and BMI or physical activity, the latter not observed with a GRS for overall T2D genetic liability. Applied to coronary artery disease (CAD) in both the WGHS and in JUPITER (N = 8,749), a randomized trial of rosuvastatin for primary prevention of CVD, component GRSs discriminate genetic risk associated with LDL-C from risk associated with reciprocal genetic effects on triglycerides and HDL-C. They also inform the pharmacogenetics of statin treatment by demonstrating that benefit from rosuvastatin is as strongly related to genetic risk from triglycerides and HDL-C as from LDL-C.

The Inflammation Biomarker GlycA Reflects Plasma N-Glycan Branching.

BACKGROUND: GlycA is a nuclear magnetic resonance (NMR) signal in plasma that correlates with inflammation and cardiovascular outcomes in large data sets. The signal is thought to originate from N-acetylglucosamine (GlcNAc) residues of branched plasma N-glycans, though direct experimental evidence is limited. Trace element concentrations affect plasma glycosylation patterns and may thereby also influence GlycA. METHODS: NMR GlycA signal was measured in plasma samples from 87 individuals and correlated with MALDI-MS N-glycomics and trace element analysis. We further evaluated the genetic association with GlycA at rs13107325, a single nucleotide polymorphism resulting in a missense variant within SLC39A8, a manganese transporter that influences N-glycan branching, both in our samples and existing genome-wide association studies data from 22 835 participants in the Women's Health Study (WHS). RESULTS: GlycA signal was correlated with both N-glycan branching (r2 ranging from 0.125-0.265; all P < 0.001) and copper concentration (r2 = 0.348, P < 0.0001). In addition, GlycA levels were associated with rs13107325 genotype in the WHS (ß [standard error of the mean] = -4.66 [1.2674], P = 0.0002). CONCLUSIONS: These results provide the first direct experimental evidence linking the GlycA NMR signal to N-glycan branching commonly associated with acute phase reactive proteins involved in inflammation.

Association between ABO haplotypes and the risk of venous thrombosis: impact on disease risk estimation.

Genetic risk score (GRS) analysis is a popular approach to derive individual risk prediction models for complex diseases. In venous thrombosis (VT), such type of analysis shall integrate information at the ABO blood group locus, which is one of the major susceptibility loci. However, there is no consensus about which single nucleotide polymorphisms (SNPs) must be investigated when properly assessing association between ABO locus and VT risk. Using comprehensive haplotype analyses of ABO blood group tagging SNPs in 5425 cases and 8445 controls from 6 studies, we demonstrate that using only rs8176719 (tagging O1) to correctly assess the impact of ABO locus on VT risk is suboptimal, because 5% of rs8176719-delG carriers do not have an increased risk of developing VT. Instead, we recommend the use of 4 SNPs, rs2519093 (tagging A1), rs1053878 (A2), rs8176743 (B), and rs41302905 (O2), when assessing the impact of ABO locus on VT risk to avoid any risk misestimation. Compared with the O1 haplotype, the A2 haplotype is associated with a modest increase in VT risk (odds ratio, âˆ¼1.2), the A1 and B haplotypes are associated with an âˆ¼1.8-fold increased risk, whereas the O2 haplotype tends to be slightly protective (odds ratio, âˆ¼0.80). In addition, although the A1 and B blood groups are associated with increased von Willebrand factor and factor VIII plasma levels, only the A1 blood group is associated with ICAM levels, but in an opposite direction, leaving additional avenues to be explored to fully understand the spectrum of biological effects mediated by ABO locus on cardiovascular traits.

Re-remembering the influence of randomized ß-carotene on cognitive decline.

The wave of individuals impacted by dementia continues to rise rapidly as worldwide lifespan increases. Dietary strategies to slow cognitive decline and prolong time to clinical dementia remain understudied, but with potentially powerful public health consequences. Indeed, previously conducted large, randomized, placebo-controlled trials of micronutrients remain an under-leveraged resource to study changes in cognitive performance. As a motivating example, we highlight an ancillary report from the Physicians' Health Study, where subjects randomized to ß-carotene (a provitamin A carotenoid) had a more attenuated change in longitudinal global cognitive performance and verbal memory, as compared to subjects randomized to placebo. Despite mechanistic evidence from cell and animal studies supporting a vitamin A-mediated role in the biology associated with cognition, limited follow-up work has been conducted. We argue that dietary factors (including provitamin A) deserve a second look, leveraging multi-omic approaches, to elucidate how they may mitigate cognitive decline and dementia risk.

Associations of Mitochondrial and Nuclear Mitochondrial Variants and Genes with Seven Metabolic Traits.

Mitochondria (MT), the major site of cellular energy production, are under dual genetic control by 37 mitochondrial DNA (mtDNA) genes and numerous nuclear genes (MT-nDNA). In the CHARGEmtDNA+ Consortium, we studied genetic associations of mtDNA and MT-nDNA associations with body mass index (BMI), waist-hip-ratio (WHR), glucose, insulin, HOMA-B, HOMA-IR, and HbA1c. This 45-cohort collaboration comprised 70,775 (insulin) to 170,202 (BMI) pan-ancestry individuals. Validation and imputation of mtDNA variants was followed by single-variant and gene-based association testing. We report two significant common variants, one in MT-ATP6 associated (p ≤ 5E-04) with WHR and one in the D-loop with glucose. Five rare variants in MT-ATP6, MT-ND5, and MT-ND6 associated with BMI, WHR, or insulin. Gene-based meta-analysis identified MT-ND3 associated with BMI (p ≤ 1E-03). We considered 2,282 MT-nDNA candidate gene associations compiled from online summary results for our traits (20 unique studies with 31 dataset consortia's genome-wide associations [GWASs]). Of these, 109 genes associated (p ≤ 1E-06) with at least 1 of our 7 traits. We assessed regulatory features of variants in the 109 genes, cis- and trans-gene expression regulation, and performed enrichment and protein-protein interactions analyses. Of the identified mtDNA and MT-nDNA genes, 79 associated with adipose measures, 49 with glucose/insulin, 13 with risk for type 2 diabetes, and 18 with cardiovascular disease, indicating for pleiotropic effects with health implications. Additionally, 21 genes related to cholesterol, suggesting additional important roles for the genes identified. Our results suggest that mtDNA and MT-nDNA genes and variants reported make important contributions to glucose and insulin metabolism, adipocyte regulation, diabetes, and cardiovascular disease.

Exome-Derived Adiponectin-Associated Variants Implicate Obesity and Lipid Biology.

Circulating levels of adiponectin, an adipocyte-secreted protein associated with cardiovascular and metabolic risk, are highly heritable. To gain insights into the biology that regulates adiponectin levels, we performed an exome array meta-analysis of 265,780 genetic variants in 67,739 individuals of European, Hispanic, African American, and East Asian ancestry. We identified 20 loci associated with adiponectin, including 11 that had been reported previously (p < 2 × 10-7). Comparison of exome array variants to regional linkage disequilibrium (LD) patterns and prior genome-wide association study (GWAS) results detected candidate variants (r2 > .60) spanning as much as 900 kb. To identify potential genes and mechanisms through which the previously unreported association signals act to affect adiponectin levels, we assessed cross-trait associations, expression quantitative trait loci in subcutaneous adipose, and biological pathways of nearby genes. Eight of the nine loci were also associated (p < 1 × 10-4) with at least one obesity or lipid trait. Candidate genes include PRKAR2A, PTH1R, and HDAC9, which have been suggested to play roles in adipocyte differentiation or bone marrow adipose tissue. Taken together, these findings provide further insights into the processes that influence circulating adiponectin levels.

Sugar-Sweetened Beverage Consumption and Plasma Lipoprotein Cholesterol, Apolipoprotein, and Lipoprotein Particle Size Concentrations in US Adults.

BACKGROUND: Prospective cohort studies have found a relation between sugar-sweetened beverage (SSB) consumption (sodas and fruit drinks) and dyslipidemia. There is limited evidence linking SSB consumption to emerging features of dyslipidemia, which can be characterized by variation in lipoprotein particle size, remnant-like particle (RLP), and apolipoprotein concentrations. OBJECTIVES: To examine the association between SSB consumption and plasma lipoprotein cholesterol, apolipoprotein, and lipoprotein particle size concentrations among US adults. METHODS: We examined participants from the Framingham Offspring Study (FOS; 1987-1995, n = 3047) and the Women's Health Study (1992, n = 26,218). Concentrations of plasma LDL cholesterol, apolipoprotein B (apoB), HDL cholesterol, apolipoprotein A1 (apoA1), triglyceride (TG), and non-HDL cholesterol, as well as total cholesterol:HDL cholesterol ratio and apoB:apoA1 ratio, were quantified in both cohorts; concentrations of apolipoprotein E, apolipoprotein C3, RLP-TG, and RLP cholesterol (RLP-C) were measured in the FOS only. Lipoprotein particle sizes were calculated from nuclear magnetic resonance signals for lipoprotein particle subclass concentrations (TG-rich lipoprotein particles [TRL-Ps]: very large, large, medium, small, and very small; LDL particles [LDL-Ps]: large, medium, and small; HDL particles [HDL-Ps]: large, medium, and small). SSB consumption was estimated from food frequency questionnaire data. We examined the associations between SSB consumption and all lipoprotein and apoprotein measures in linear regression models, adjusting for confounding factors such as lifestyle, diet, and traditional lipoprotein risk factors. RESULTS: SSB consumption was positively associated with LDL cholesterol, apoB, TG, RLP-TG, RLP-C, and non-HDL cholesterol concentrations and total cholesterol:HDL cholesterol and apoB:apoA1 ratios; and negatively associated with HDL cholesterol and apoA1 concentrations (P-trend range: <0.0001 to 0.008). After adjustment for traditional lipoprotein risk factors, SSB consumers had smaller LDL-P and HDL-P sizes; lower concentrations of large LDL-Ps and medium HDL-Ps; and higher concentrations of small LDL-Ps, small HDL-Ps, and large TRL-Ps (P-trend range: <0.0001 to 0.001). CONCLUSIONS: Higher SSB consumption was associated with multiple emerging features of dyslipidemia that have been linked to higher cardiometabolic risk in US adults.

Gene-educational attainment interactions in a multi-ancestry genome-wide meta-analysis identify novel blood pressure loci.

Educational attainment is widely used as a surrogate for socioeconomic status (SES). Low SES is a risk factor for hypertension and high blood pressure (BP). To identify novel BP loci, we performed multi-ancestry meta-analyses accounting for gene-educational attainment interactions using two variables, "Some College" (yes/no) and "Graduated College" (yes/no). Interactions were evaluated using both a 1 degree of freedom (DF) interaction term and a 2DF joint test of genetic and interaction effects. Analyses were performed for systolic BP, diastolic BP, mean arterial pressure, and pulse pressure. We pursued genome-wide interrogation in Stage 1 studies (N = 117 438) and follow-up on promising variants in Stage 2 studies (N = 293 787) in five ancestry groups. Through combined meta-analyses of Stages 1 and 2, we identified 84 known and 18 novel BP loci at genome-wide significance level (P < 5 × 10-8). Two novel loci were identified based on the 1DF test of interaction with educational attainment, while the remaining 16 loci were identified through the 2DF joint test of genetic and interaction effects. Ten novel loci were identified in individuals of African ancestry. Several novel loci show strong biological plausibility since they involve physiologic systems implicated in BP regulation. They include genes involved in the central nervous system-adrenal signaling axis (ZDHHC17, CADPS, PIK3C2G), vascular structure and function (GNB3, CDON), and renal function (HAS2 and HAS2-AS1, SLIT3). Collectively, these findings suggest a role of educational attainment or SES in further dissection of the genetic architecture of BP.

Identifying blood pressure loci whose effects are modulated by multiple lifestyle exposures.

Although multiple lifestyle exposures simultaneously impact blood pressure (BP) and cardiovascular health, most analysis so far has considered each single lifestyle exposure (e.g., smoking) at a time. Here, we exploit gene-multiple lifestyle exposure interactions to find novel BP loci. For each of 6,254 Framingham Heart Study participants, we computed lifestyle risk score (LRS) value by aggregating the risk of four lifestyle exposures (smoking, alcohol, education, and physical activity) on BP. Using the LRS, we performed genome-wide gene-environment interaction analysis in systolic and diastolic BP using the joint 2 degree of freedom (DF) and 1 DF interaction tests. We identified one genome-wide significant (p < 5 × 10-8 ) and 11 suggestive (p < 1 × 10-6 ) loci. Gene-environment analysis using single lifestyle exposures identified only one of the 12 loci. Nine of the 12 BP loci detected were novel. Loci detected by the LRS were located within or nearby genes with biologically plausible roles in the pathophysiology of hypertension, including KALRN, VIPR2, SNX1, and DAPK2. Our results suggest that simultaneous consideration of multiple lifestyle exposures in gene-environment interaction analysis can identify additional loci missed by single lifestyle approaches.

A genome-wide association study of bitter and sweet beverage consumption.

Except for drinking water, most beverages taste bitter or sweet. Taste perception and preferences are heritable and determinants of beverage choice and consumption. Consumption of several bitter- and sweet-tasting beverages has been implicated in development of major chronic diseases. We performed a genome-wide association study (GWAS) of self-reported bitter and sweet beverage consumption among ~370 000 participants of European ancestry, using a two-staged analysis design. Bitter beverages included coffee, tea, grapefruit juice, red wine, liquor and beer. Sweet beverages included artificially and sugar sweetened beverages (SSBs) and non-grapefruit juices. Five loci associated with total bitter beverage consumption were replicated (in/near GCKR, ABCG2, AHR, POR and CYP1A1/2). No locus was replicated for total sweet beverage consumption. Sub-phenotype analyses targeting the alcohol, caffeine and sweetener components of beverages yielded additional loci: (i) four loci for bitter alcoholic beverages (GCKR, KLB, ADH1B and AGBL2); (ii) five loci for bitter non-alcoholic beverages (ANXA9, AHR, POR, CYP1A1/2 and CSDC2); (iii) 10 loci for coffee; six novel loci (SEC16B, TMEM18, OR8U8, AKAP6, MC4R and SPECC1L-ADORA2A); (iv) FTO for SSBs. Of these 17 replicated loci, 12 have been associated with total alcohol consumption, coffee consumption, plasma caffeine metabolites or BMI in previous GWAS; none was involved in known sweet and bitter taste transduction pathways. Our study suggests that genetic variants related to alcohol consumption, coffee consumption and obesity were primary genetic determinants of bitter and sweet beverage consumption. Whether genetic variants related to taste perception are associated with beverage consumption remains to be determined.

A multi-ancestry genome-wide study incorporating gene-smoking interactions identifies multiple new loci for pulse pressure and mean arterial pressure.

Elevated blood pressure (BP), a leading cause of global morbidity and mortality, is influenced by both genetic and lifestyle factors. Cigarette smoking is one such lifestyle factor. Across five ancestries, we performed a genome-wide gene-smoking interaction study of mean arterial pressure (MAP) and pulse pressure (PP) in 129 913 individuals in stage 1 and follow-up analysis in 480 178 additional individuals in stage 2. We report here 136 loci significantly associated with MAP and/or PP. Of these, 61 were previously published through main-effect analysis of BP traits, 37 were recently reported by us for systolic BP and/or diastolic BP through gene-smoking interaction analysis and 38 were newly identified (P < 5 × 10-8, false discovery rate < 0.05). We also identified nine new signals near known loci. Of the 136 loci, 8 showed significant interaction with smoking status. They include CSMD1 previously reported for insulin resistance and BP in the spontaneously hypertensive rats. Many of the 38 new loci show biologic plausibility for a role in BP regulation. SLC26A7 encodes a chloride/bicarbonate exchanger expressed in the renal outer medullary collecting duct. AVPR1A is widely expressed, including in vascular smooth muscle cells, kidney, myocardium and brain. FHAD1 is a long non-coding RNA overexpressed in heart failure. TMEM51 was associated with contractile function in cardiomyocytes. CASP9 plays a central role in cardiomyocyte apoptosis. Identified only in African ancestry were 30 novel loci. Our findings highlight the value of multi-ancestry investigations, particularly in studies of interaction with lifestyle factors, where genomic and lifestyle differences may contribute to novel findings.