Elevated lipoprotein(a) predicts high risk of cardiovascular disease among a modest proportion of healthy individuals, an issue that complicates screening guidelines. To examine spline models, clinica Show more
Elevated lipoprotein(a) predicts high risk of cardiovascular disease among a modest proportion of healthy individuals, an issue that complicates screening guidelines. To examine spline models, clinical thresholds, and percentiles of baseline lipoprotein(a) levels as 30-year determinants of cardiovascular risk. This cohort study was conducted among female health professionals participating in the Women's Health Study, who were followed up prospectively from 1993 to 2023. Women without cardiovascular disease, cancer, and other major chronic illnesses had blood samples taken at baseline. All individuals with lipoprotein(a) measurements and/or of European ancestry with genotype information for the LPA rs3798220 variation were included. Data analyses were performed from January through April 2025. Continuously valued baseline lipoprotein(a), lipoprotein(a) clinical thresholds and percentiles, and LPA rs3798220 genotypes known to predict lipoprotein(a) levels among individuals of European ancestry. The primary outcomes were incident major cardiovascular events, coronary heart disease, ischemic stroke, and cardiovascular death. Age- and multivariable-adjusted cause-specific Cox models were used to calculated hazard ratios for the cardiovascular outcomes. The hypothesis was formulated after collection of the data. A total of 27 748 women with baseline lipoprotein(a) measurements and 23 279 women of European ancestry with rs3798220 genotype information were included (median [IQR] age, 53 [49-60] years), among whom 3707 and 3165 major cardiovascular events, respectively, accrued during a median (IQR) follow-up period of 27.8 (22.8-29.4) years. Among women with lipoprotein(a) measurements, lipoprotein(a) levels above 30 mg/dL or the 75th percentile (31 mg/dL) were associated with increased 30-year risk of major cardiovascular events and coronary heart disease. Levels above 120 mg/dL or the 99th percentile (131 mg/dL) were associated with increased risk of ischemic stroke and cardiovascular death. Multivariable adjusted hazard ratios for levels above 120 mg/dL vs below 10 mg/dL or above the 99th percentile vs below the 50th percentile (11 mg/dL) were 1.54 (95% CI, 1.24-1.92) and 1.74 (95% CI, 1.35-2.25) for major cardiovascular events, 1.80 (95% CI, 1.36-2.37) and 2.06 (95% CI, 1.49-2.84) for coronary heart disease, 1.41 (95% CI, 0.93-2.15) and 1.85 (95% CI, 1.17-2.93) for ischemic stroke, and 1.63 (95% CI, 1.16-2.28) and 1.86 (95% CI, 1.26-2.72) for cardiovascular death, respectively. Among women with genotype information, rs3798220 minor allele carriers had a higher risk of major cardiovascular events. Per the results of this cohort study, very high lipoprotein(a) levels correlated with increased 30-year risk of cardiovascular disease among healthy women. Screening for elevated lipoprotein(a) in the general population may be warranted. Show less
Discovery and translation of gene-environment interactions (GxEs) influencing clinical outcomes is limited by low statistical power and poor mechanistic understanding. Molecular omics data may help ad Show more
Discovery and translation of gene-environment interactions (GxEs) influencing clinical outcomes is limited by low statistical power and poor mechanistic understanding. Molecular omics data may help address these limitations, but their incorporation into GxE testing requires principled analytic approaches. We focused on genetic modification of the established mechanistic link between dietary long-chain omega-3 fatty acid (dN3FA) intake, plasma N3FA (pN3FA), and chronic inflammation as measured by high sensitivity CRP (hsCRP). We considered an approach that decomposes the overall genetic effect modification into components upstream and downstream of a molecular mediator to increase the potential to discover gene-N3FA interactions. Simulations demonstrated improved power of the upstream and downstream tests compared to the standard approach when the molecular mediator for many biologically plausible scenarios. The approach was applied in the UK Biobank (N = 188,700) with regression models that used measures of dN3FA (based on fish and fish oil intake), pN3FA (% of total fatty acids measured by nuclear magnetic resonance), and hsCRP. Mediation analysis showed that pN3FA fully mediated the dN3FA-hsCRP main effect relationship. Next, we separately tested modification of the dN3FA-hsCRP ("standard"), dN3FA-pN3FA ("upstream"), and pN3FA-hsCRP ("downstream") associations. The known FADS1-3 locus variant rs174535 reached p = 1.6 × 10 Show less
Atrial fibrillation (AF) is a prevalent and morbid abnormality of the heart rhythm with a strong genetic component. Here, we meta-analyzed genome and exome sequencing data from 36 studies that include Show more
Atrial fibrillation (AF) is a prevalent and morbid abnormality of the heart rhythm with a strong genetic component. Here, we meta-analyzed genome and exome sequencing data from 36 studies that included 52,416 AF cases and 277,762 controls. In burden tests of rare coding variation, we identified novel associations between AF and the genes MYBPC3, LMNA, PKP2, FAM189A2 and KDM5B. We further identified associations between AF and rare structural variants owing to deletions in CTNNA3 and duplications of GATA4. We broadly replicated our findings in independent samples from MyCode, deCODE and UK Biobank. Finally, we found that CRISPR knockout of KDM5B in stem-cell-derived atrial cardiomyocytes led to a shortening of the action potential duration and widespread transcriptomic dysregulation of genes relevant to atrial homeostasis and conduction. Our results highlight the contribution of rare coding and structural variants to AF, including genetic links between AF and cardiomyopathies, and expand our understanding of the rare variant architecture for this common arrhythmia. Show less
Obesity is a major public health crisis associated with high mortality rates. Previous genome-wide association studies (GWAS) investigating body mass index (BMI) have largely relied on imputed data fr Show more
Obesity is a major public health crisis associated with high mortality rates. Previous genome-wide association studies (GWAS) investigating body mass index (BMI) have largely relied on imputed data from European individuals. This study leveraged whole-genome sequencing (WGS) data from 88,873 participants from the Trans-Omics for Precision Medicine (TOPMed) Program, of which 51% were of non-European population groups. We discovered 18 BMI-associated signals (P < 5 × 10 Show less
Discovery and translation of gene-environment interactions (GxEs) influencing clinical outcomes is limited by low statistical power and poor mechanistic understanding. Molecular omics data may help ad Show more
Discovery and translation of gene-environment interactions (GxEs) influencing clinical outcomes is limited by low statistical power and poor mechanistic understanding. Molecular omics data may help address these limitations, but their incorporation into GxE testing requires principled analytic approaches. We focused on genetic modification of the established mechanistic link between dietary long-chain omega-3 fatty acid (dN3FA) intake, plasma N3FA (pN3FA), and chronic inflammation as measured by high sensitivity CRP (hsCRP). We considered an approach that decomposes the overall genetic effect modification into components upstream and downstream of a molecular mediator to increase the potential to discover gene-N3FA interactions. Simulations demonstrated improved power of the upstream and downstream tests compared to the standard approach when the molecular mediator for many biologically plausible scenarios. The approach was applied in the UK Biobank (N = 188,700) with regression models that used measures of dN3FA (based on fish and fish oil intake), pN3FA (% of total fatty acids measured by nuclear magnetic resonance), and hsCRP. Mediation analysis showed that pN3FA fully mediated the dN3FA-hsCRP main effect relationship. Next, we separately tested modification of the dN3FA-hsCRP ("standard"), dN3FA-pN3FA ("upstream"), and pN3FA-hsCRP ("downstream") associations. The known Show less
Background The risk of atherosclerotic cardiovascular disease (ASCVD) increases sharply with age. Some older individuals, however, remain unaffected despite high predicted risk. These individuals may Show more
Background The risk of atherosclerotic cardiovascular disease (ASCVD) increases sharply with age. Some older individuals, however, remain unaffected despite high predicted risk. These individuals may carry cardioprotective genetic variants that contribute to resilience. Our aim was to assess whether asymptomatic older individuals without prevalent ASCVD carry cardioprotective genetic variants that contribute to ASCVD resilience. Methods and Results We performed a genome-wide association study using a 10-year predicted ASCVD risk score as a quantitative trait, calculated only in asymptomatic older individuals aged ≥70 years without prevalent ASCVD. Our discovery genome-wide association study of N=12 031 ASCVD event-free individuals from the ASPREE (Aspirin in Reducing Events in the Elderly) trial identified 2 independent variants, rs9939224 ( Show less
Identifying genetic determinants of reproductive success may highlight mechanisms underlying fertility and identify alleles under present-day selection. Using data in 785,604 individuals of European a Show more
Identifying genetic determinants of reproductive success may highlight mechanisms underlying fertility and identify alleles under present-day selection. Using data in 785,604 individuals of European ancestry, we identified 43 genomic loci associated with either number of children ever born (NEB) or childlessness. These loci span diverse aspects of reproductive biology, including puberty timing, age at first birth, sex hormone regulation, endometriosis and age at menopause. Missense variants in ARHGAP27 were associated with higher NEB but shorter reproductive lifespan, suggesting a trade-off at this locus between reproductive ageing and intensity. Other genes implicated by coding variants include PIK3IP1, ZFP82 and LRP4, and our results suggest a new role for the melanocortin 1 receptor (MC1R) in reproductive biology. As NEB is one component of evolutionary fitness, our identified associations indicate loci under present-day natural selection. Integration with data from historical selection scans highlighted an allele in the FADS1/2 gene locus that has been under selection for thousands of years and remains so today. Collectively, our findings demonstrate that a broad range of biological mechanisms contribute to reproductive success. Show less
Endometriosis is a common condition associated with debilitating pelvic pain and infertility. A genome-wide association study meta-analysis, including 60,674 cases and 701,926 controls of European and Show more
Endometriosis is a common condition associated with debilitating pelvic pain and infertility. A genome-wide association study meta-analysis, including 60,674 cases and 701,926 controls of European and East Asian descent, identified 42 genome-wide significant loci comprising 49 distinct association signals. Effect sizes were largest for stage 3/4 disease, driven by ovarian endometriosis. Identified signals explained up to 5.01% of disease variance and regulated expression or methylation of genes in endometrium and blood, many of which were associated with pain perception/maintenance (SRP14/BMF, GDAP1, MLLT10, BSN and NGF). We observed significant genetic correlations between endometriosis and 11 pain conditions, including migraine, back and multisite chronic pain (MCP), as well as inflammatory conditions, including asthma and osteoarthritis. Multitrait genetic analyses identified substantial sharing of variants associated with endometriosis and MCP/migraine. Targeted investigations of genetically regulated mechanisms shared between endometriosis and other pain conditions are needed to aid the development of new treatments and facilitate early symptomatic intervention. Show less
Obesity is a major public health crisis associated with high mortality rates. Previous genome-wide association studies (GWAS) investigating body mass index (BMI) have largely relied on imputed data fr Show more
Obesity is a major public health crisis associated with high mortality rates. Previous genome-wide association studies (GWAS) investigating body mass index (BMI) have largely relied on imputed data from European individuals. This study leveraged whole-genome sequencing (WGS) data from 88,873 participants from the Trans-Omics for Precision Medicine (TOPMed) Program, of which 51% were of non-European population groups. We discovered 18 BMI-associated signals ( Show less
While periodontal disease (PD) has been associated with type 2 diabetes (T2D) and osteoporosis, the underlying genetic mechanisms for these associations remain largely unknown. The aim of this study i Show more
While periodontal disease (PD) has been associated with type 2 diabetes (T2D) and osteoporosis, the underlying genetic mechanisms for these associations remain largely unknown. The aim of this study is to apply cross-trait genetic analyses to investigate the potentially shared biology among PD, T2D, and bone mineral density (BMD) by assessing pairwise genetic correlations and searching for shared polymorphisms. We applied cross-trait genetic analyses leveraging genome-wide association study (GWAS) summary statistics for: Periodontitis/loose teeth from the UKBB/GLIDE consortium (PerioLT, N=506594), T2D from the DIAGRAM consortium (N Significant genetic correlations were identified between PerioLT/T2D (Rg=0.23; SE=0.04; p=7.4e This integrative approach identified one colocalized locus and 14 additional candidate loci that are shared between T2D and PD/oral health by comparing effects across PD, T2D and BMD. Future research is needed to independently validate our findings. Show less
Dietary intake is a major contributor to the global obesity epidemic and represents a complex behavioural phenotype that is partially affected by innate biological differences. Here, we present a mult Show more
Dietary intake is a major contributor to the global obesity epidemic and represents a complex behavioural phenotype that is partially affected by innate biological differences. Here, we present a multivariate genome-wide association analysis of overall variation in dietary intake to account for the correlation between dietary carbohydrate, fat and protein in 282,271 participants of European ancestry from the UK Biobank (n = 191,157) and Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium (n = 91,114), and identify 26 distinct genome-wide significant loci. Dietary intake signals map exclusively to specific brain regions and are enriched for genes expressed in specialized subtypes of GABAergic, dopaminergic and glutamatergic neurons. We identified two main clusters of genetic variants for overall variation in dietary intake that were differently associated with obesity and coronary artery disease. These results enhance the biological understanding of interindividual differences in dietary intake by highlighting neural mechanisms, supporting functional follow-up experiments and possibly providing new avenues for the prevention and treatment of prevalent complex metabolic diseases. Show less
Migraine affects over a billion individuals worldwide but its genetic underpinning remains largely unknown. Here, we performed a genome-wide association study of 102,084 migraine cases and 771,257 con Show more
Migraine affects over a billion individuals worldwide but its genetic underpinning remains largely unknown. Here, we performed a genome-wide association study of 102,084 migraine cases and 771,257 controls and identified 123 loci, of which 86 are previously unknown. These loci provide an opportunity to evaluate shared and distinct genetic components in the two main migraine subtypes: migraine with aura and migraine without aura. Stratification of the risk loci using 29,679 cases with subtype information indicated three risk variants that seem specific for migraine with aura (in HMOX2, CACNA1A and MPPED2), two that seem specific for migraine without aura (near SPINK2 and near FECH) and nine that increase susceptibility for migraine regardless of subtype. The new risk loci include genes encoding recent migraine-specific drug targets, namely calcitonin gene-related peptide (CALCA/CALCB) and serotonin 1F receptor (HTR1F). Overall, genomic annotations among migraine-associated variants were enriched in both vascular and central nervous system tissue/cell types, supporting unequivocally that neurovascular mechanisms underlie migraine pathophysiology. Show less
To evaluate phenotypic and genetic relationships between migraine and lipoprotein subfractions. We evaluated phenotypic associations between migraine and 19 lipoprotein subfraction measures in the Wom Show more
To evaluate phenotypic and genetic relationships between migraine and lipoprotein subfractions. We evaluated phenotypic associations between migraine and 19 lipoprotein subfraction measures in the Women's Genome Health Study (n = 22,788). We then investigated genetic relationships between these traits using summary statistics from the International Headache Genetics Consortium for migraine (n There was a significant phenotypic association (odds ratio 1.27 [95% confidence interval 1.12-1.44]) and a significant genetic correlation at 0.18 ( The study supports the association between certain lipoprotein subfractions, especially for TRLP, and migraine in populations of European ancestry. The corresponding shared genetic components may help identify potential targets for future migraine therapeutics. This study provides Class I evidence that migraine is significantly associated with some lipoprotein subfractions. Show less
The HERMES (HEart failure Molecular Epidemiology for Therapeutic targetS) consortium aims to identify the genomic and molecular basis of heart failure. The consortium currently includes 51 studies fro Show more
The HERMES (HEart failure Molecular Epidemiology for Therapeutic targetS) consortium aims to identify the genomic and molecular basis of heart failure. The consortium currently includes 51 studies from 11 countries, including 68 157 heart failure cases and 949 888 controls, with data on heart failure events and prognosis. All studies collected biological samples and performed genome-wide genotyping of common genetic variants. The enrolment of subjects into participating studies ranged from 1948 to the present day, and the median follow-up following heart failure diagnosis ranged from 2 to 116 months. Forty-nine of 51 individual studies enrolled participants of both sexes; in these studies, participants with heart failure were predominantly male (34-90%). The mean age at diagnosis or ascertainment across all studies ranged from 54 to 84 years. Based on the aggregate sample, we estimated 80% power to genetic variant associations with risk of heart failure with an odds ratio of ≥1.10 for common variants (allele frequency ≥ 0.05) and ≥1.20 for low-frequency variants (allele frequency 0.01-0.05) at P < 5 × 10 HERMES is a global collaboration aiming to (i) identify the genetic determinants of heart failure; (ii) generate insights into the causal pathways leading to heart failure and enable genetic approaches to target prioritization; and (iii) develop genomic tools for disease stratification and risk prediction. Show less
ChREBP (carbohydrate responsive element binding protein) is a transcription factor that responds to sugar consumption. Sugar-sweetened beverage (SSB) consumption and genetic variants in the Data from Show more
ChREBP (carbohydrate responsive element binding protein) is a transcription factor that responds to sugar consumption. Sugar-sweetened beverage (SSB) consumption and genetic variants in the Data from 11 cohorts from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium (N=63 599) and the UK Biobank (N=59 220) were used to quantify associations of SSB consumption, genetic variants, and their interaction on HDL-C and triglyceride concentrations using linear regression models. A total of 1606 single nucleotide polymorphisms within or near In a meta-analysis, rs71556729 was significantly associated with higher HDL-C concentrations only among the highest SSB consumers (β, 2.12 [95% CI, 1.16-3.07] mg/dL per allele; Our results identified genetic variants in the Show less
Genome-wide association studies have identified multiple genomic loci associated with coronary artery disease, but most are common variants in non-coding regions that provide limited information on ca Show more
Genome-wide association studies have identified multiple genomic loci associated with coronary artery disease, but most are common variants in non-coding regions that provide limited information on causal genes and etiology of the disease. To overcome the limited scope that common variants provide, we focused our investigation on low-frequency and rare sequence variations primarily residing in coding regions of the genome. Using samples of individuals of European ancestry from ten cohorts within the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, both cross-sectional and prospective analyses were conducted to examine associations between genetic variants and myocardial infarction (MI), coronary heart disease (CHD), and all-cause mortality following these events. For prevalent events, a total of 27,349 participants of European ancestry, including 1831 prevalent MI cases and 2518 prevalent CHD cases were used. For incident cases, a total of 55,736 participants of European ancestry were included (3,031 incident MI cases and 5,425 incident CHD cases). There were 1,860 all-cause deaths among the 3,751 MI and CHD cases from six cohorts that contributed to the analysis of all-cause mortality. Single variant and gene-based analyses were performed separately in each cohort and then meta-analyzed for each outcome. A low-frequency intronic variant (rs988583) in PLCL1 was significantly associated with prevalent MI (OR = 1.80, 95% confidence interval: 1.43, 2.27; P = 7.12 × 10-7). We conducted gene-based burden tests for genes with a cumulative minor allele count (cMAC) ≥ 5 and variants with minor allele frequency (MAF) < 5%. TMPRSS5 and LDLRAD1 were significantly associated with prevalent MI and CHD, respectively, and RC3H2 and ANGPTL4 were significantly associated with incident MI and CHD, respectively. No loci were significantly associated with all-cause mortality following a MI or CHD event. This study identified one known locus (ANGPTL4) and four new loci (PLCL1, RC3H2, TMPRSS5, and LDLRAD1) associated with cardiovascular disease risk that warrant further investigation. Show less
Alcohol intake influences plasma lipid levels, and such effects may be moderated by genetic variants. We aimed to characterize the role of aggregated rare and low-frequency protein-coding variants in Show more
Alcohol intake influences plasma lipid levels, and such effects may be moderated by genetic variants. We aimed to characterize the role of aggregated rare and low-frequency protein-coding variants in gene by alcohol consumption interactions associated with fasting plasma lipid levels. In the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium, fasting plasma triglycerides and high- and low-density lipoprotein cholesterol were measured in 34 153 individuals with European ancestry from 5 discovery studies and 32 277 individuals from 6 replication studies. Rare and low-frequency functional protein-coding variants (minor allele frequency, ≤5%) measured by an exome array were aggregated by genes and evaluated by a gene-environment interaction test and a joint test of genetic main and gene-environment interaction effects. Two dichotomous self-reported alcohol consumption variables, current drinker, defined as any recurrent drinking behavior, and regular drinker, defined as the subset of current drinkers who consume at least 2 drinks per week, were considered. We discovered and replicated 21 gene-lipid associations at 13 known lipid loci through the joint test. Eight loci ( In conclusion, this study applied new gene-based statistical approaches and suggested that rare and low-frequency genetic variants interacted with alcohol consumption on lipid levels. Show less
Blood pressure (BP) was inconsistently associated with migraine and the mechanisms of BP-lowering medications in migraine prophylaxis are unknown. Leveraging large-scale summary statistics for migrain Show more
Blood pressure (BP) was inconsistently associated with migraine and the mechanisms of BP-lowering medications in migraine prophylaxis are unknown. Leveraging large-scale summary statistics for migraine (N Show less
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 Show more
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. Show less
Genome-wide association studies (GWAS) have identified >250 loci for body mass index (BMI), implicating pathways related to neuronal biology. Most GWAS loci represent clusters of common, noncoding var Show more
Genome-wide association studies (GWAS) have identified >250 loci for body mass index (BMI), implicating pathways related to neuronal biology. Most GWAS loci represent clusters of common, noncoding variants from which pinpointing causal genes remains challenging. Here we combined data from 718,734 individuals to discover rare and low-frequency (minor allele frequency (MAF) < 5%) coding variants associated with BMI. We identified 14 coding variants in 13 genes, of which 8 variants were in genes (ZBTB7B, ACHE, RAPGEF3, RAB21, ZFHX3, ENTPD6, ZFR2 and ZNF169) newly implicated in human obesity, 2 variants were in genes (MC4R and KSR2) previously observed to be mutated in extreme obesity and 2 variants were in GIPR. The effect sizes of rare variants are ~10 times larger than those of common variants, with the largest effect observed in carriers of an MC4R mutation introducing a stop codon (p.Tyr35Ter, MAF = 0.01%), who weighed ~7 kg more than non-carriers. Pathway analyses based on the variants associated with BMI confirm enrichment of neuronal genes and provide new evidence for adipocyte and energy expenditure biology, widening the potential of genetically supported therapeutic targets in obesity. Show less
We screened variants on an exome-focused genotyping array in >300,000 participants (replication in >280,000 participants) and identified 444 independent variants in 250 loci significantly associated w Show more
We screened variants on an exome-focused genotyping array in >300,000 participants (replication in >280,000 participants) and identified 444 independent variants in 250 loci significantly associated with total cholesterol (TC), high-density-lipoprotein cholesterol (HDL-C), low-density-lipoprotein cholesterol (LDL-C), and/or triglycerides (TG). At two loci (JAK2 and A1CF), experimental analysis in mice showed lipid changes consistent with the human data. We also found that: (i) beta-thalassemia trait carriers displayed lower TC and were protected from coronary artery disease (CAD); (ii) excluding the CETP locus, there was not a predictable relationship between plasma HDL-C and risk for age-related macular degeneration; (iii) only some mechanisms of lowering LDL-C appeared to increase risk for type 2 diabetes (T2D); and (iv) TG-lowering alleles involved in hepatic production of TG-rich lipoproteins (TM6SF2 and PNPLA3) tracked with higher liver fat, higher risk for T2D, and lower risk for CAD, whereas TG-lowering alleles involved in peripheral lipolysis (LPL and ANGPTL4) had no effect on liver fat but decreased risks for both T2D and CAD. Show less
The ages of puberty, first sexual intercourse and first birth signify the onset of reproductive ability, behavior and success, respectively. In a genome-wide association study of 125,667 UK Biobank pa Show more
The ages of puberty, first sexual intercourse and first birth signify the onset of reproductive ability, behavior and success, respectively. In a genome-wide association study of 125,667 UK Biobank participants, we identify 38 loci associated (P < 5 × 10(-8)) with age at first sexual intercourse. These findings were taken forward in 241,910 men and women from Iceland and 20,187 women from the Women's Genome Health Study. Several of the identified loci also exhibit associations (P < 5 × 10(-8)) with other reproductive and behavioral traits, including age at first birth (variants in or near ESR1 and RBM6-SEMA3F), number of children (CADM2 and ESR1), irritable temperament (MSRA) and risk-taking propensity (CADM2). Mendelian randomization analyses infer causal influences of earlier puberty timing on earlier first sexual intercourse, earlier first birth and lower educational attainment. In turn, likely causal consequences of earlier first sexual intercourse include reproductive, educational, psychiatric and cardiometabolic outcomes. Show less
Circulating trans fatty acids (TFAs), which cannot be synthesized by humans, are linked to adverse health outcomes. Although TFAs are obtained from diet, little is known about subsequent influences (e Show more
Circulating trans fatty acids (TFAs), which cannot be synthesized by humans, are linked to adverse health outcomes. Although TFAs are obtained from diet, little is known about subsequent influences (e.g., relating to incorporation, metabolism, or intercompetition with other fatty acids) that could alter circulating concentrations and possibly modulate or mediate impacts on health. The objective was to elucidate novel biologic pathways that may influence circulating TFAs by evaluating associations between common genetic variation and TFA biomarkers. We performed meta-analyses using 7 cohorts of European-ancestry participants (n = 8013) having measured genome-wide variation in single-nucleotide polymorphisms (SNPs) and circulating TFA biomarkers (erythrocyte or plasma phospholipids), including trans-16:1n-7, total trans-18:1, trans/cis-18:2, cis/trans-18:2, and trans/trans-18:2. We further evaluated SNPs with genome-wide significant associations among African Americans (n = 1082), Chinese Americans (n = 669), and Hispanic Americans (n = 657) from 2 of these cohorts. Among European-ancestry participants, 31 SNPs in or near the fatty acid desaturase (FADS) 1 and 2 cluster were associated with cis/trans-18:2; a top hit was rs174548 (β = 0.0035, P = 4.90 × 10(-15)), an SNP previously associated with circulating n-3 and n-6 polyunsaturated fatty acid concentrations. No significant association was identified for other TFAs. rs174548 in FADS1/2 was also associated with cis/trans-18:2 in Hispanic Americans (β = 0.0053, P = 1.05 × 10(-6)) and Chinese Americans (β = 0.0028, P = 0.002) but not African Americans (β = 0.0009, P = 0.34); however, in African Americans, fine mapping identified a top hit in FADS2 associated with cis/trans-18:2 (rs174579: β = 0.0118, P = 4.05 × 10(-5)). The association between rs174548 and cis/trans-18:2 remained significant after further adjustment for individual circulating n-3 and n-6 fatty acids, except arachidonic acid. After adjustment for arachidonic acid concentrations, the association between rs174548 and cis/trans-18:2 was nearly eliminated in European-ancestry participants (β-coefficient reduced by 86%), with similar reductions in Hispanic Americans and Chinese Americans. Our findings provide novel evidence for genetic regulation of cis/trans-18:2 by the FADS1/2 cluster and suggest that this regulation may be influenced/mediated by concentrations of arachidonic acid, an n-6 polyunsaturated fat. Show less
Tissue concentrations of omega-3 fatty acids may reduce cardiovascular disease risk, and genetic variants are associated with circulating fatty acids concentrations. Whether dietary fatty acids intera Show more
Tissue concentrations of omega-3 fatty acids may reduce cardiovascular disease risk, and genetic variants are associated with circulating fatty acids concentrations. Whether dietary fatty acids interact with genetic variants to modify circulating omega-3 fatty acids is unclear. We evaluated interactions between genetic variants and fatty acid intakes for circulating alpha-linoleic acid, eicosapentaenoic acid, docosahexaenoic acid, and docosapentaenoic acid. We conducted meta-analyses (N = 11 668) evaluating interactions between dietary fatty acids and genetic variants (rs174538 and rs174548 in FADS1 (fatty acid desaturase 1), rs7435 in AGPAT3 (1-acyl-sn-glycerol-3-phosphate), rs4985167 in PDXDC1 (pyridoxal-dependent decarboxylase domain-containing 1), rs780094 in GCKR (glucokinase regulatory protein), and rs3734398 in ELOVL2 (fatty acid elongase 2)). Stratification by measurement compartment (plasma versus erthyrocyte) revealed compartment-specific interactions between FADS1 rs174538 and rs174548 and dietary alpha-linolenic acid and linoleic acid for docosahexaenoic acid and docosapentaenoic acid. Our findings reinforce earlier reports that genetically based differences in circulating fatty acids may be partially due to differences in the conversion of fatty acid precursors. Further, fatty acids measurement compartment may modify gene-diet relationships, and considering compartment may improve the detection of gene-fatty acids interactions for circulating fatty acid outcomes. Show less
Coffee, a major dietary source of caffeine, is among the most widely consumed beverages in the world and has received considerable attention regarding health risks and benefits. We conducted a genome- Show more
Coffee, a major dietary source of caffeine, is among the most widely consumed beverages in the world and has received considerable attention regarding health risks and benefits. We conducted a genome-wide (GW) meta-analysis of predominately regular-type coffee consumption (cups per day) among up to 91,462 coffee consumers of European ancestry with top single-nucleotide polymorphisms (SNPs) followed-up in ~30 062 and 7964 coffee consumers of European and African-American ancestry, respectively. Studies from both stages were combined in a trans-ethnic meta-analysis. Confirmed loci were examined for putative functional and biological relevance. Eight loci, including six novel loci, met GW significance (log10Bayes factor (BF)>5.64) with per-allele effect sizes of 0.03-0.14 cups per day. Six are located in or near genes potentially involved in pharmacokinetics (ABCG2, AHR, POR and CYP1A2) and pharmacodynamics (BDNF and SLC6A4) of caffeine. Two map to GCKR and MLXIPL genes related to metabolic traits but lacking known roles in coffee consumption. Enhancer and promoter histone marks populate the regions of many confirmed loci and several potential regulatory SNPs are highly correlated with the lead SNP of each. SNP alleles near GCKR, MLXIPL, BDNF and CYP1A2 that were associated with higher coffee consumption have previously been associated with smoking initiation, higher adiposity and fasting insulin and glucose but lower blood pressure and favorable lipid, inflammatory and liver enzyme profiles (P<5 × 10(-8)).Our genetic findings among European and African-American adults reinforce the role of caffeine in mediating habitual coffee consumption and may point to molecular mechanisms underlying inter-individual variability in pharmacological and health effects of coffee. Show less
Pleiotropic genetic variants have independent effects on different phenotypes. C-reactive protein (CRP) is associated with several cardiometabolic phenotypes. Shared genetic backgrounds may partially Show more
Pleiotropic genetic variants have independent effects on different phenotypes. C-reactive protein (CRP) is associated with several cardiometabolic phenotypes. Shared genetic backgrounds may partially underlie these associations. We conducted a genome-wide analysis to identify the shared genetic background of inflammation and cardiometabolic phenotypes using published genome-wide association studies (GWAS). We also evaluated whether the pleiotropic effects of such loci were biological or mediated in nature. First, we examined whether 283 common variants identified for 10 cardiometabolic phenotypes in GWAS are associated with CRP level. Second, we tested whether 18 variants identified for serum CRP are associated with 10 cardiometabolic phenotypes. We used a Bonferroni corrected p-value of 1.1×10-04 (0.05/463) as a threshold of significance. We evaluated the independent pleiotropic effect on both phenotypes using individual level data from the Women Genome Health Study. Evaluating the genetic overlap between inflammation and cardiometabolic phenotypes, we found 13 pleiotropic regions. Additional analyses showed that 6 regions (APOC1, HNF1A, IL6R, PPP1R3B, HNF4A and IL1F10) appeared to have a pleiotropic effect on CRP independent of the effects on the cardiometabolic phenotypes. These included loci where individuals carrying the risk allele for CRP encounter higher lipid levels and risk of type 2 diabetes. In addition, 5 regions (GCKR, PABPC4, BCL7B, FTO and TMEM18) had an effect on CRP largely mediated through the cardiometabolic phenotypes. In conclusion, our results show genetic pleiotropy among inflammation and cardiometabolic phenotypes. In addition to reverse causation, our data suggests that pleiotropic genetic variants partially underlie the association between CRP and cardiometabolic phenotypes. Show less
Metabolic syndrome (MetS) has become a health and financial burden worldwide. The MetS definition captures clustering of risk factors that predict higher risk for diabetes mellitus and cardiovascular Show more
Metabolic syndrome (MetS) has become a health and financial burden worldwide. The MetS definition captures clustering of risk factors that predict higher risk for diabetes mellitus and cardiovascular disease. Our study hypothesis is that additional to genes influencing individual MetS risk factors, genetic variants exist that influence MetS and inflammatory markers forming a predisposing MetS genetic network. To test this hypothesis a staged approach was undertaken. (a) We analyzed 17 metabolic and inflammatory traits in more than 85,500 participants from 14 large epidemiological studies within the Cross Consortia Pleiotropy Group. Individuals classified with MetS (NCEP definition), versus those without, showed on average significantly different levels for most inflammatory markers studied. (b) Paired average correlations between 8 metabolic traits and 9 inflammatory markers from the same studies as above, estimated with two methods, and factor analyses on large simulated data, helped in identifying 8 combinations of traits for follow-up in meta-analyses, out of 130,305 possible combinations between metabolic traits and inflammatory markers studied. (c) We performed correlated meta-analyses for 8 metabolic traits and 6 inflammatory markers by using existing GWAS published genetic summary results, with about 2.5 million SNPs from twelve predominantly largest GWAS consortia. These analyses yielded 130 unique SNPs/genes with pleiotropic associations (a SNP/gene associating at least one metabolic trait and one inflammatory marker). Of them twenty-five variants (seven loci newly reported) are proposed as MetS candidates. They map to genes MACF1, KIAA0754, GCKR, GRB14, COBLL1, LOC646736-IRS1, SLC39A8, NELFE, SKIV2L, STK19, TFAP2B, BAZ1B, BCL7B, TBL2, MLXIPL, LPL, TRIB1, ATXN2, HECTD4, PTPN11, ZNF664, PDXDC1, FTO, MC4R and TOMM40. Based on large data evidence, we conclude that inflammation is a feature of MetS and several gene variants show pleiotropic genetic associations across phenotypes and might explain a part of MetS correlated genetic architecture. These findings warrant further functional investigation. Show less
The strong observational association between total homocysteine (tHcy) concentrations and risk of coronary artery disease (CAD) and the null associations in the homocysteine-lowering trials have promp Show more
The strong observational association between total homocysteine (tHcy) concentrations and risk of coronary artery disease (CAD) and the null associations in the homocysteine-lowering trials have prompted the need to identify genetic variants associated with homocysteine concentrations and risk of CAD. We tested whether common genetic polymorphisms associated with variation in tHcy are also associated with CAD. We conducted a meta-analysis of genome-wide association studies (GWAS) on tHcy concentrations in 44,147 individuals of European descent. Polymorphisms associated with tHcy (P < 10(⁻⁸) were tested for association with CAD in 31,400 cases and 92,927 controls. Common variants at 13 loci, explaining 5.9% of the variation in tHcy, were associated with tHcy concentrations, including 6 novel loci in or near MMACHC (2.1 × 10⁻⁹), SLC17A3 (1.0 × 10⁻⁸), GTPB10 (1.7 × 10⁻⁸), CUBN (7.5 × 10⁻¹⁰), HNF1A (1.2 × 10⁻¹²)), and FUT2 (6.6 × 10⁻⁹), and variants previously reported at or near the MTHFR, MTR, CPS1, MUT, NOX4, DPEP1, and CBS genes. Individuals within the highest 10% of the genotype risk score (GRS) had 3-μmol/L higher mean tHcy concentrations than did those within the lowest 10% of the GRS (P = 1 × 10⁻³⁶). The GRS was not associated with risk of CAD (OR: 1.01; 95% CI: 0.98, 1.04; P = 0.49). We identified several novel loci that influence plasma tHcy concentrations. Overall, common genetic variants that influence plasma tHcy concentrations are not associated with risk of CAD in white populations, which further refutes the causal relevance of moderately elevated tHcy concentrations and tHcy-related pathways for CAD. Show less
Numerous obesity loci have been identified using genome-wide association studies. A UK study indicated that physical activity may attenuate the cumulative effect of 12 of these loci, but replication s Show more
Numerous obesity loci have been identified using genome-wide association studies. A UK study indicated that physical activity may attenuate the cumulative effect of 12 of these loci, but replication studies are lacking. Therefore, we tested whether the aggregate effect of these loci is diminished in adults of European ancestry reporting high levels of physical activity. Twelve obesity-susceptibility loci were genotyped or imputed in 111,421 participants. A genetic risk score (GRS) was calculated by summing the BMI-associated alleles of each genetic variant. Physical activity was assessed using self-administered questionnaires. Multiplicative interactions between the GRS and physical activity on BMI were tested in linear and logistic regression models in each cohort, with adjustment for age, age(2), sex, study center (for multicenter studies), and the marginal terms for physical activity and the GRS. These results were combined using meta-analysis weighted by cohort sample size. The meta-analysis yielded a statistically significant GRS × physical activity interaction effect estimate (Pinteraction = 0.015). However, a statistically significant interaction effect was only apparent in North American cohorts (n = 39,810, Pinteraction = 0.014 vs. n = 71,611, Pinteraction = 0.275 for Europeans). In secondary analyses, both the FTO rs1121980 (Pinteraction = 0.003) and the SEC16B rs10913469 (Pinteraction = 0.025) variants showed evidence of SNP × physical activity interactions. This meta-analysis of 111,421 individuals provides further support for an interaction between physical activity and a GRS in obesity disposition, although these findings hinge on the inclusion of cohorts from North America, indicating that these results are either population-specific or non-causal. Show less
Chronic kidney disease (CKD) is an important public health problem with a genetic component. We performed genome-wide association studies in up to 130,600 European ancestry participants overall, and s Show more
Chronic kidney disease (CKD) is an important public health problem with a genetic component. We performed genome-wide association studies in up to 130,600 European ancestry participants overall, and stratified for key CKD risk factors. We uncovered 6 new loci in association with estimated glomerular filtration rate (eGFR), the primary clinical measure of CKD, in or near MPPED2, DDX1, SLC47A1, CDK12, CASP9, and INO80. Morpholino knockdown of mpped2 and casp9 in zebrafish embryos revealed podocyte and tubular abnormalities with altered dextran clearance, suggesting a role for these genes in renal function. By providing new insights into genes that regulate renal function, these results could further our understanding of the pathogenesis of CKD. Show less