Dietary intake of plant- and animal-derived protein and incident cardiovascular diseases: the pan-European EPIC-CVD case-cohort study

BACKGROUND

Epidemiological evidence suggests that a potential association between dietary protein intake and cardiovascular disease (CVD) may depend on the protein source, that is, plant- or animal-derived, but past research was limited and inconclusive.

OBJECTIVES

To evaluate the association of dietary plant- or animal-derived protein consumption with risk of CVD, and its components ischemic heart disease (IHD) and stroke.

METHODS

This analysis in the European Prospective Investigation into Cancer and Nutrition (EPIC)-CVD case-cohort study included 16,244 incident CVD cases (10,784 IHD and 6423 stroke cases) and 15,141 subcohort members from 7 European countries. We investigated the association of estimated dietary protein intake with CVD, IHD, and stroke (total, fatal, and nonfatal) using multivariable-adjusted Prentice-weighted Cox regression. We estimated isocaloric substitutions of replacing fats and carbohydrates with plant- or animal-derived protein and replacing food-specific animal protein with plant protein. Multiplicative interactions between dietary protein and prespecified variables were tested.

RESULTS

Neither plant- nor animal-derived protein intake was associated with incident CVD, IHD, or stroke in adjusted analyses without or with macronutrient-specified substitution analyses. Higher plant-derived protein intake was associated with 22% lower total stroke incidence among never smokers [HR 0.78, 95% confidence intervals (CI): 0.62, 0.99], but not among current smokers (HR 1.08, 95% CI: 0.83, 1.40, P-interaction = 0.004). Moreover, higher plant-derived protein (per 3% total energy) when replacing red meat protein (HR 0.52, 95% CI: 0.31, 0.88), processed meat protein (HR 0.39, 95% CI: 0.17, 0.90), and dairy protein (HR 0.54, 95% CI: 0.30, 0.98) was associated with lower incidence of fatal stroke.

CONCLUSION

Plant- or animal-derived protein intake was not associated with overall CVD. However, the association of plant-derived protein consumption with lower total stroke incidence among nonsmokers, and with lower incidence of fatal stroke highlights the importance of investigating CVD subtypes and potential interactions. These observations warrant further investigation in diverse populations with varying macronutrient intakes and dietary patterns.

Genome-wide analysis in over 1 million individuals of European ancestry yields improved polygenic risk scores for blood pressure traits

Hypertension affects more than one billion people worldwide. Here we identify 113 novel loci, reporting a total of 2,103 independent genetic signals (P < 5 × 10-8) from the largest single-stage blood pressure (BP) genome-wide association study to date (n = 1,028,980 European individuals). These associations explain more than 60% of single nucleotide polymorphism-based BP heritability. Comparing top versus bottom deciles of polygenic risk scores (PRSs) reveals clinically meaningful differences in BP (16.9 mmHg systolic BP, 95% CI, 15.5-18.2 mmHg, P = 2.22 × 10-126) and more than a sevenfold higher odds of hypertension risk (odds ratio, 7.33; 95% CI, 5.54-9.70; P = 4.13 × 10-44) in an independent dataset. Adding PRS into hypertension-prediction models increased the area under the receiver operating characteristic curve (AUROC) from 0.791 (95% CI, 0.781-0.801) to 0.826 (95% CI, 0.817-0.836, ∆AUROC, 0.035, P = 1.98 × 10-34). We compare the 2,103 loci results in non-European ancestries and show significant PRS associations in a large African-American sample. Secondary analyses implicate 500 genes previously unreported for BP. Our study highlights the role of increasingly large genomic studies for precision health research.

Genome-wide characterization of circulating metabolic biomarkers

Genome-wide association analyses using high-throughput metabolomics platforms have led to novel insights into the biology of human metabolism1-7. This detailed knowledge of the genetic determinants of systemic metabolism has been pivotal for uncovering how genetic pathways influence biological mechanisms and complex diseases8-11. Here we present a genome-wide association study for 233 circulating metabolic traits quantified by nuclear magnetic resonance spectroscopy in up to 136,016 participants from 33 cohorts. We identify more than 400 independent loci and assign probable causal genes at two-thirds of these using manual curation of plausible biological candidates. We highlight the importance of sample and participant characteristics that can have significant effects on genetic associations. We use detailed metabolic profiling of lipoprotein- and lipid-associated variants to better characterize how known lipid loci and novel loci affect lipoprotein metabolism at a granular level. We demonstrate the translational utility of comprehensively phenotyped molecular data, characterizing the metabolic associations of intrahepatic cholestasis of pregnancy. Finally, we observe substantial genetic pleiotropy for multiple metabolic pathways and illustrate the importance of careful instrument selection in Mendelian randomization analysis, revealing a putative causal relationship between acetone and hypertension. Our publicly available results provide a foundational resource for the community to examine the role of metabolism across diverse diseases.

Genetic Associations of Circulating Cardiovascular Proteins With Gestational Hypertension and Preeclampsia

Importance

Hypertensive disorders of pregnancy (HDPs), including gestational hypertension and preeclampsia, are important contributors to maternal morbidity and mortality worldwide. In addition, women with HDPs face an elevated long-term risk of cardiovascular disease.

Objective

To identify proteins in the circulation associated with HDPs.

Design, Setting, and Participants

Two-sample mendelian randomization (MR) tested the associations of genetic instruments for cardiovascular disease-related proteins with gestational hypertension and preeclampsia. In downstream analyses, a systematic review of observational data was conducted to evaluate the identified proteins' dynamics across gestation in hypertensive vs normotensive pregnancies, and phenome-wide MR analyses were performed to identify potential non-HDP-related effects associated with the prioritized proteins. Genetic association data for cardiovascular disease-related proteins were obtained from the Systematic and Combined Analysis of Olink Proteins (SCALLOP) consortium. Genetic association data for the HDPs were obtained from recent European-ancestry genome-wide association study meta-analyses for gestational hypertension and preeclampsia. Study data were analyzed October 2022 to October 2023.

Exposures

Genetic instruments for 90 candidate proteins implicated in cardiovascular diseases, constructed using cis-protein quantitative trait loci (cis-pQTLs).

Main Outcomes and Measures

Gestational hypertension and preeclampsia.

Results

Genetic association data for cardiovascular disease-related proteins were obtained from 21 758 participants from the SCALLOP consortium. Genetic association data for the HDPs were obtained from 393 238 female individuals (8636 cases and 384 602 controls) for gestational hypertension and 606 903 female individuals (16 032 cases and 590 871 controls) for preeclampsia. Seventy-five of 90 proteins (83.3%) had at least 1 valid cis-pQTL. Of those, 10 proteins (13.3%) were significantly associated with HDPs. Four were robust to sensitivity analyses for gestational hypertension (cluster of differentiation 40, eosinophil cationic protein [ECP], galectin 3, N-terminal pro-brain natriuretic peptide [NT-proBNP]), and 2 were robust for preeclampsia (cystatin B, heat shock protein 27 [HSP27]). Consistent with the MR findings, observational data revealed that lower NT-proBNP (0.76- to 0.88-fold difference vs no HDPs) and higher HSP27 (2.40-fold difference vs no HDPs) levels during the first trimester of pregnancy were associated with increased risk of HDPs, as were higher levels of ECP (1.60-fold difference vs no HDPs). Phenome-wide MR analyses identified 37 unique non-HDP-related protein-disease associations, suggesting potential on-target effects associated with interventions lowering HDP risk through the identified proteins.

Conclusions and Relevance

Study findings suggest genetic associations of 4 cardiovascular disease-related proteins with gestational hypertension and 2 associated with preeclampsia. Future studies are required to test the efficacy of targeting the corresponding pathways to reduce HDP risk.

Genetic drivers of heterogeneity in type 2 diabetes pathophysiology

Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes1,2 and molecular mechanisms that are often specific to cell type3,4. Here, to characterize the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study data from 2,535,601 individuals (39.7% not of European ancestry), including 428,452 cases of T2D. We identify 1,289 independent association signals at genome-wide significance (P < 5 × 10-8) that map to 611 loci, of which 145 loci are, to our knowledge, previously unreported. We define eight non-overlapping clusters of T2D signals that are characterized by distinct profiles of cardiometabolic trait associations. These clusters are differentially enriched for cell-type-specific regions of open chromatin, including pancreatic islets, adipocytes, endothelial cells and enteroendocrine cells. We build cluster-specific partitioned polygenic scores5 in a further 279,552 individuals of diverse ancestry, including 30,288 cases of T2D, and test their association with T2D-related vascular outcomes. Cluster-specific partitioned polygenic scores are associated with coronary artery disease, peripheral artery disease and end-stage diabetic nephropathy across ancestry groups, highlighting the importance of obesity-related processes in the development of vascular outcomes. Our findings show the value of integrating multi-ancestry genome-wide association study data with single-cell epigenomics to disentangle the aetiological heterogeneity that drives the development and progression of T2D. This might offer a route to optimize global access to genetically informed diabetes care.

Genetically predicted plasma cortisol and common chronic diseases: A Mendelian randomization study

OBJECTIVE

Cushing's syndrome is characterized by hypercortisolaemia and is frequently accompanied by comorbidities such as type 2 diabetes, hypertension, osteoporosis, depression and schizophrenia. It is unclear whether moderate but lifelong hypercortisolaemia is causally associated with these diseases in the general population. We aimed to address this research gap using a Mendelian randomization approach.

METHODS

We used three cortisol-associated genetic variants in the SERPINA6/SERPINA1 region as genetic instruments in a two-sample, inverse-variance-weighted Mendelian randomization analysis. We obtained summary-level statistics for cortisol and disease outcomes from publicly available genetic consortia, and meta-analysed them as appropriate. We conducted a multivariable Mendelian randomization analysis to assess potential mediating effects.

RESULTS

A 1 standard deviation higher genetically predicted plasma cortisol was associated with greater odds of hypertension (odds ratio: 1.12; 95% confidence interval [CI]: 1.05-1.18) as well as higher systolic blood pressure (mean difference [MD]: 0.03 SD change; 95% CI: 0.01-0.05) and diastolic blood pressure (MD: 0.03 SD change; 95% CI: 0.01-0.04). There was no evidence of association with type 2 diabetes, osteoporosis, depression and schizophrenia. The association with hypertension was attenuated upon adjustment for waist circumference, suggesting potential mediation through central obesity.

CONCLUSION

There is strong evidence for a causal association between plasma cortisol and greater risk for hypertension, potentially mediated by obesity.

Dietary amino acids and risk of stroke subtypes: a prospective analysis of 356,000 participants in seven European countries

PURPOSE

Previously reported associations of protein-rich foods with stroke subtypes have prompted interest in the assessment of individual amino acids. We examined the associations of dietary amino acids with risks of ischaemic and haemorrhagic stroke in the EPIC study.

METHODS

We analysed data from 356,142 participants from seven European countries. Dietary intakes of 19 individual amino acids were assessed using validated country-specific dietary questionnaires, calibrated using additional 24-h dietary recalls. Multivariable-adjusted Cox regression models were used to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) of ischaemic and haemorrhagic stroke in relation to the intake of each amino acid. The role of blood pressure as a potential mechanism was assessed in 267,642 (75%) participants.

RESULTS

After a median follow-up of 12.9 years, 4295 participants had an ischaemic stroke and 1375 participants had a haemorrhagic stroke. After correction for multiple testing, a higher intake of proline (as a percent of total protein) was associated with a 12% lower risk of ischaemic stroke (HR per 1 SD higher intake 0.88; 95% CI 0.82, 0.94). The association persisted after mutual adjustment for all other amino acids, systolic and diastolic blood pressure. The inverse associations of isoleucine, leucine, valine, phenylalanine, threonine, tryptophan, glutamic acid, serine and tyrosine with ischaemic stroke were each attenuated with adjustment for proline intake. For haemorrhagic stroke, no statistically significant associations were observed in the continuous analyses after correcting for multiple testing.

CONCLUSION

Higher proline intake may be associated with a lower risk of ischaemic stroke, independent of other dietary amino acids and blood pressure.

Estimating dose-response relationships for vitamin D with coronary heart disease, stroke, and all-cause mortality: observational and Mendelian randomisation analyses

BACKGROUND

Randomised trials of vitamin D supplementation for cardiovascular disease and all-cause mortality have generally reported null findings. However, generalisability of results to individuals with low vitamin D status is unclear. We aimed to characterise dose-response relationships between 25-hydroxyvitamin D (25[OH]D) concentrations and risk of coronary heart disease, stroke, and all-cause mortality in observational and Mendelian randomisation frameworks.

METHODS

Observational analyses were undertaken using data from 33 prospective studies comprising 500 962 individuals with no known history of coronary heart disease or stroke at baseline. Mendelian randomisation analyses were performed in four population-based cohort studies (UK Biobank, EPIC-CVD, and two Copenhagen population-based studies) comprising 386 406 middle-aged individuals of European ancestries, including 33 546 people who developed coronary heart disease, 18 166 people who had a stroke, and 27 885 people who died. Primary outcomes were coronary heart disease, defined as fatal ischaemic heart disease (International Classification of Diseases 10th revision code I20-I25) or non-fatal myocardial infarction (I21-I23); stroke, defined as any cerebrovascular disease (I60-I69); and all-cause mortality.

FINDINGS

Observational analyses suggested inverse associations between incident coronary heart disease, stroke, and all-cause mortality outcomes with 25(OH)D concentration at low 25(OH)D concentrations. In population-wide genetic analyses, there were no associations of genetically predicted 25(OH)D with coronary heart disease (odds ratio [OR] per 10 nmol/L higher genetically-predicted 25(OH)D concentration 0·98, 95% CI 0·95-1·01), stroke (1·01, [0·97-1·05]), or all-cause mortality (0·99, 0·95-1·02). Null findings were also observed in genetic analyses for cause-specific mortality outcomes, and in stratified genetic analyses for all outcomes at all observed levels of 25(OH)D concentrations.

INTERPRETATION

Stratified Mendelian randomisation analyses suggest a lack of causal relationship for 25(OH)D concentrations with both cardiovascular and mortality outcomes for individuals at all levels of 25(OH)D. Our findings suggest that substantial reductions in mortality and cardiovascular morbidity due to long-term low-dose vitamin D supplementation are unlikely even if targeted at individuals with low vitamin D status.

FUNDING

British Heart Foundation, Medical Research Council, National Institute for Health Research, Health Data Research UK, Cancer Research UK, and International Agency for Research on Cancer.

Mendelian randomization for cardiovascular diseases: principles and applications

Large-scale genome-wide association studies conducted over the last decade have uncovered numerous genetic variants associated with cardiometabolic traits and risk factors. These discoveries have enabled the Mendelian randomization (MR) design, which uses genetic variation as a natural experiment to improve causal inferences from observational data. By analogy with the random assignment of treatment in randomized controlled trials, the random segregation of genetic alleles when DNA is transmitted from parents to offspring at gamete formation is expected to reduce confounding in genetic associations. Mendelian randomization analyses make a set of assumptions that must hold for valid results. Provided that the assumptions are well justified for the genetic variants that are employed as instrumental variables, MR studies can inform on whether a putative risk factor likely has a causal effect on the disease or not. Mendelian randomization has been increasingly applied over recent years to predict the efficacy and safety of existing and novel drugs targeting cardiovascular risk factors and to explore the repurposing potential of available drugs. This review article describes the principles of the MR design and some applications in cardiovascular epidemiology.

Life expectancy associated with different ages at diagnosis of type 2 diabetes in high-income countries: 23 million person-years of observation

BACKGROUND

The prevalence of type 2 diabetes is increasing rapidly, particularly among younger age groups. Estimates suggest that people with diabetes die, on average, 6 years earlier than people without diabetes. We aimed to provide reliable estimates of the associations between age at diagnosis of diabetes and all-cause mortality, cause-specific mortality, and reductions in life expectancy.

METHODS

For this observational study, we conducted a combined analysis of individual-participant data from 19 high-income countries using two large-scale data sources: the Emerging Risk Factors Collaboration (96 cohorts, median baseline years 1961-2007, median latest follow-up years 1980-2013) and the UK Biobank (median baseline year 2006, median latest follow-up year 2020). We calculated age-adjusted and sex-adjusted hazard ratios (HRs) for all-cause mortality according to age at diagnosis of diabetes using data from 1 515 718 participants, in whom deaths were recorded during 23·1 million person-years of follow-up. We estimated cumulative survival by applying age-specific HRs to age-specific death rates from 2015 for the USA and the EU.

FINDINGS

For participants with diabetes, we observed a linear dose-response association between earlier age at diagnosis and higher risk of all-cause mortality compared with participants without diabetes. HRs were 2·69 (95% CI 2·43-2·97) when diagnosed at 30-39 years, 2·26 (2·08-2·45) at 40-49 years, 1·84 (1·72-1·97) at 50-59 years, 1·57 (1·47-1·67) at 60-69 years, and 1·39 (1·29-1·51) at 70 years and older. HRs per decade of earlier diagnosis were similar for men and women. Using death rates from the USA, a 50-year-old individual with diabetes died on average 14 years earlier when diagnosed aged 30 years, 10 years earlier when diagnosed aged 40 years, or 6 years earlier when diagnosed aged 50 years than an individual without diabetes. Using EU death rates, the corresponding estimates were 13, 9, or 5 years earlier.

INTERPRETATION

Every decade of earlier diagnosis of diabetes was associated with about 3-4 years of lower life expectancy, highlighting the need to develop and implement interventions that prevent or delay the onset of diabetes and to intensify the treatment of risk factors among young adults diagnosed with diabetes.

FUNDING

British Heart Foundation, Medical Research Council, National Institute for Health and Care Research, and Health Data Research UK.

Genetics of circulating inflammatory proteins identifies drivers of immune-mediated disease risk and therapeutic targets

Circulating proteins have important functions in inflammation and a broad range of diseases. To identify genetic influences on inflammation-related proteins, we conducted a genome-wide protein quantitative trait locus (pQTL) study of 91 plasma proteins measured using the Olink Target platform in 14,824 participants. We identified 180 pQTLs (59 cis, 121 trans). Integration of pQTL data with eQTL and disease genome-wide association studies provided insight into pathogenesis, implicating lymphotoxin-α in multiple sclerosis. Using Mendelian randomization (MR) to assess causality in disease etiology, we identified both shared and distinct effects of specific proteins across immune-mediated diseases, including directionally discordant effects of CD40 on risk of rheumatoid arthritis versus multiple sclerosis and inflammatory bowel disease. MR implicated CXCL5 in the etiology of ulcerative colitis (UC) and we show elevated gut CXCL5 transcript expression in patients with UC. These results identify targets of existing drugs and provide a powerful resource to facilitate future drug target prioritization.

A genome-wide association study of blood cell morphology identifies cellular proteins implicated in disease aetiology

Blood cells contain functionally important intracellular structures, such as granules, critical to immunity and thrombosis. Quantitative variation in these structures has not been subjected previously to large-scale genetic analysis. We perform genome-wide association studies of 63 flow-cytometry derived cellular phenotypes-including cell-type specific measures of granularity, nucleic acid content and reactivity-in 41,515 participants in the INTERVAL study. We identify 2172 distinct variant-trait associations, including associations near genes coding for proteins in organelles implicated in inflammatory and thrombotic diseases. By integrating with epigenetic data we show that many intracellular structures are likely to be determined in immature precursor cells. By integrating with proteomic data we identify the transcription factor FOG2 as an early regulator of platelet formation and α-granularity. Finally, we show that colocalisation of our associations with disease risk signals can suggest aetiological cell-types-variants in IL2RA and ITGA4 respectively mirror the known effects of daclizumab in multiple sclerosis and vedolizumab in inflammatory bowel disease.

An empirical investigation into the impact of winner's curse on estimates from Mendelian randomization

INTRODUCTION

Genetic associations for variants identified through genome-wide association studies (GWASs) tend to be overestimated in the original discovery data set as, if the association was underestimated, the variant may not have been detected. This bias, known as winner's curse, can affect Mendelian randomization estimates, but its severity and potential impact are unclear.

METHODS

We performed an empirical investigation to assess the potential bias from winner's curse in practice. We considered Mendelian randomization estimates for the effect of body mass index (BMI) on coronary artery disease risk. We randomly divided a UK Biobank data set 100 times into three equally sized subsets. The first subset was treated as the 'discovery GWAS'. We compared genetic associations estimated in the discovery GWAS to those estimated in the other subsets for each of the 100 iterations.

RESULTS

For variants associated with BMI at P < 5 × 10-8 in at least one iteration, genetic associations with BMI were up to 5-fold greater in iterations in which the variant was associated with BMI at P < 5 × 10-8 compared with its mean association across all iterations. If the minimum P-value for association with BMI was P = 10-13 or lower, then this inflation was <25%. Mendelian randomization estimates were affected by winner's curse bias. However, bias did not materially affect results; all analyses indicated a deleterious effect of BMI on coronary artery disease risk.

CONCLUSIONS

Winner's curse can bias Mendelian randomization estimates, although its practical impact may not be substantial. If avoiding sample overlap is infeasible, analysts should consider performing a sensitivity analysis based on variants strongly associated with the exposure.

A multi-ancestry polygenic risk score improves risk prediction for coronary artery disease

Identification of individuals at highest risk of coronary artery disease (CAD)-ideally before onset-remains an important public health need. Prior studies have developed genome-wide polygenic scores to enable risk stratification, reflecting the substantial inherited component to CAD risk. Here we develop a new and significantly improved polygenic score for CAD, termed GPSMult, that incorporates genome-wide association data across five ancestries for CAD (>269,000 cases and >1,178,000 controls) and ten CAD risk factors. GPSMult strongly associated with prevalent CAD (odds ratio per standard deviation 2.14, 95% confidence interval 2.10-2.19, P < 0.001) in UK Biobank participants of European ancestry, identifying 20.0% of the population with 3-fold increased risk and conversely 13.9% with 3-fold decreased risk as compared with those in the middle quintile. GPSMult was also associated with incident CAD events (hazard ratio per standard deviation 1.73, 95% confidence interval 1.70-1.76, P < 0.001), identifying 3% of healthy individuals with risk of future CAD events equivalent to those with existing disease and significantly improving risk discrimination and reclassification. Across multiethnic, external validation datasets inclusive of 33,096, 124,467, 16,433 and 16,874 participants of African, European, Hispanic and South Asian ancestry, respectively, GPSMult demonstrated increased strength of associations across all ancestries and outperformed all available previously published CAD polygenic scores. These data contribute a new GPSMult for CAD to the field and provide a generalizable framework for how large-scale integration of genetic association data for CAD and related traits from diverse populations can meaningfully improve polygenic risk prediction.

Risk Factors of Secondary Cardiovascular Events in a Multi-Ethnic Asian Population with Acute Myocardial Infarction: A Retrospective Cohort Study from Malaysia

This retrospective cohort study investigated the incidence and risk factors of major adverse cardiovascular events (MACE) after 1 year of first-documented myocardial infarctions (MIs) in a multi-ethnic Asian population. Secondary MACE were observed in 231 (14.3%) individuals, including 92 (5.7%) cardiovascular-related deaths. Both histories of hypertension and diabetes were associated with secondary MACE after adjustment for age, sex, and ethnicity (HR 1.60 [95%CI 1.22-2.12] and 1.46 [95%CI 1.09-1.97], respectively). With further adjustments for traditional risk factors, individuals with conduction disturbances demonstrated higher risks of MACE: new left-bundle branch block (HR 2.86 [95%CI 1.15-6.55]), right-bundle branch block (HR 2.09 [95%CI 1.02-4.29]), and second-degree heart block (HR 2.45 [95%CI 0.59-10.16]). These associations were broadly similar across different age, sex, and ethnicity groups, although somewhat greater for history of hypertension and BMI among women versus men, for HbA1c control in individuals aged >50 years, and for LVEF ≤ 40% in those with Indian versus Chinese or Bumiputera ethnicities. Several traditional and cardiac risk factors are associated with a higher risk of secondary major adverse cardiovascular events. In addition to hypertension and diabetes, the identification of conduction disturbances in individuals with first-onset MI may be useful for the risk stratification of high-risk individuals.

South Asian medical cohorts reveal strong founder effects and high rates of homozygosity

The benefits of large-scale genetic studies for healthcare of the populations studied are well documented, but these genetic studies have traditionally ignored people from some parts of the world, such as South Asia. Here we describe whole genome sequence (WGS) data from 4806 individuals recruited from the healthcare delivery systems of Pakistan, India and Bangladesh, combined with WGS from 927 individuals from isolated South Asian populations. We characterize population structure in South Asia and describe a genotyping array (SARGAM) and imputation reference panel that are optimized for South Asian genomes. We find evidence for high rates of reproductive isolation, endogamy and consanguinity that vary across the subcontinent and that lead to levels of rare homozygotes that reach 100 times that seen in outbred populations. Founder effects increase the power to associate functional variants with disease processes and make South Asia a uniquely powerful place for population-scale genetic studies.

Genome-wide association meta-analysis of spontaneous coronary artery dissection identifies risk variants and genes related to artery integrity and tissue-mediated coagulation

Spontaneous coronary artery dissection (SCAD) is an understudied cause of myocardial infarction primarily affecting women. It is not known to what extent SCAD is genetically distinct from other cardiovascular diseases, including atherosclerotic coronary artery disease (CAD). Here we present a genome-wide association meta-analysis (1,917 cases and 9,292 controls) identifying 16 risk loci for SCAD. Integrative functional annotations prioritized genes that are likely to be regulated in vascular smooth muscle cells and artery fibroblasts and implicated in extracellular matrix biology. One locus containing the tissue factor gene F3, which is involved in blood coagulation cascade initiation, appears to be specific for SCAD risk. Several associated variants have diametrically opposite associations with CAD, suggesting that shared biological processes contribute to both diseases, but through different mechanisms. We also infer a causal role for high blood pressure in SCAD. Our findings provide novel pathophysiological insights involving arterial integrity and tissue-mediated coagulation in SCAD and set the stage for future specific therapeutics and preventions.

Trans-ethnic polygenic risk scores for body mass index: An international hundred K+ cohorts consortium study

BACKGROUND

While polygenic risk scores hold significant promise in estimating an individual's risk of developing a complex trait such as obesity, their application in the clinic has, to date, been limited by a lack of data from non-European populations. As a collaboration model of the International Hundred K+ Cohorts Consortium (IHCC), we endeavored to develop a globally applicable trans-ethnic PRS for body mass index (BMI) through this relatively new international effort.

METHODS

The polygenic risk score (PRS) model was developed, trained and tested at the Center for Applied Genomics (CAG) of The Children's Hospital of Philadelphia (CHOP) based on a BMI meta-analysis from the GIANT consortium. The validated PRS models were subsequently disseminated to the participating sites. Scores were generated by each site locally on their cohorts and summary statistics returned to CAG for final analysis.

RESULTS

We show that in the absence of a well powered trans-ethnic GWAS from which to derive marker SNPs and effect estimates for PRS, trans-ethnic scores can be generated from European ancestry GWAS using Bayesian approaches such as LDpred, by adjusting the summary statistics using trans-ethnic linkage disequilibrium reference panels. The ported trans-ethnic scores outperform population specific-PRS across all non-European ancestry populations investigated including East Asians and three-way admixed Brazilian cohort.

CONCLUSIONS

Here we show that for a truly polygenic trait such as BMI adjusting the summary statistics of a well powered European ancestry study using trans-ethnic LD reference results in a score that is predictive across a range of ancestries including East Asians and three-way admixed Brazilians.

An atlas of genetic scores to predict multi-omic traits

The use of omic modalities to dissect the molecular underpinnings of common diseases and traits is becoming increasingly common. But multi-omic traits can be genetically predicted, which enables highly cost-effective and powerful analyses for studies that do not have multi-omics1. Here we examine a large cohort (the INTERVAL study2; n = 50,000 participants) with extensive multi-omic data for plasma proteomics (SomaScan, n = 3,175; Olink, n = 4,822), plasma metabolomics (Metabolon HD4, n = 8,153), serum metabolomics (Nightingale, n = 37,359) and whole-blood Illumina RNA sequencing (n = 4,136), and use machine learning to train genetic scores for 17,227 molecular traits, including 10,521 that reach Bonferroni-adjusted significance. We evaluate the performance of genetic scores through external validation across cohorts of individuals of European, Asian and African American ancestries. In addition, we show the utility of these multi-omic genetic scores by quantifying the genetic control of biological pathways and by generating a synthetic multi-omic dataset of the UK Biobank3 to identify disease associations using a phenome-wide scan. We highlight a series of biological insights with regard to genetic mechanisms in metabolism and canonical pathway associations with disease; for example, JAK-STAT signalling and coronary atherosclerosis. Finally, we develop a portal ( https://www.omicspred.org/ ) to facilitate public access to all genetic scores and validation results, as well as to serve as a platform for future extensions and enhancements of multi-omic genetic scores.

A nutritional biomarker score of the Mediterranean diet and incident type 2 diabetes: Integrated analysis of data from the MedLey randomised controlled trial and the EPIC-InterAct case-cohort study

BACKGROUND

Self-reported adherence to the Mediterranean diet has been modestly inversely associated with incidence of type 2 diabetes (T2D) in cohort studies. There is uncertainty about the validity and magnitude of this association due to subjective reporting of diet. The association has not been evaluated using an objectively measured biomarker of the Mediterranean diet.

METHODS AND FINDINGS

We derived a biomarker score based on 5 circulating carotenoids and 24 fatty acids that discriminated between the Mediterranean or habitual diet arms of a parallel design, 6-month partial-feeding randomised controlled trial (RCT) conducted between 2013 and 2014, the MedLey trial (128 participants out of 166 randomised). We applied this biomarker score in an observational study, the European Prospective Investigation into Cancer and Nutrition (EPIC)-InterAct case-cohort study, to assess the association of the score with T2D incidence over an average of 9.7 years of follow-up since the baseline (1991 to 1998). We included 22,202 participants, of whom 9,453 were T2D cases, with relevant biomarkers from an original case-cohort of 27,779 participants sampled from a cohort of 340,234 people. As a secondary measure of the Mediterranean diet, we used a score estimated from dietary-self report. Within the trial, the biomarker score discriminated well between the 2 arms; the cross-validated C-statistic was 0.88 (95% confidence interval (CI) 0.82 to 0.94). The score was inversely associated with incident T2D in EPIC-InterAct: the hazard ratio (HR) per standard deviation of the score was 0.71 (95% CI: 0.65 to 0.77) following adjustment for sociodemographic, lifestyle and medical factors, and adiposity. In comparison, the HR per standard deviation of the self-reported Mediterranean diet was 0.90 (95% CI: 0.86 to 0.95). Assuming the score was causally associated with T2D, higher adherence to the Mediterranean diet in Western European adults by 10 percentiles of the score was estimated to reduce the incidence of T2D by 11% (95% CI: 7% to 14%). The study limitations included potential measurement error in nutritional biomarkers, unclear specificity of the biomarker score to the Mediterranean diet, and possible residual confounding.

CONCLUSIONS

These findings suggest that objectively assessed adherence to the Mediterranean diet is associated with lower risk of T2D and that even modestly higher adherence may have the potential to reduce the population burden of T2D meaningfully.

TRIAL REGISTRATION

Australian New Zealand Clinical Trials Registry (ANZCTR) ACTRN12613000602729 https://www.anzctr.org.au/Trial/Registration/TrialReview.aspx?id=363860.

Multi-ancestry genome-wide study in >2.5 million individuals reveals heterogeneity in mechanistic pathways of type 2 diabetes and complications

Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes. To characterise the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study (GWAS) data from 2,535,601 individuals (39.7% non-European ancestry), including 428,452 T2D cases. We identify 1,289 independent association signals at genome-wide significance (P<5×10-8) that map to 611 loci, of which 145 loci are previously unreported. We define eight non-overlapping clusters of T2D signals characterised by distinct profiles of cardiometabolic trait associations. These clusters are differentially enriched for cell-type specific regions of open chromatin, including pancreatic islets, adipocytes, endothelial, and enteroendocrine cells. We build cluster-specific partitioned genetic risk scores (GRS) in an additional 137,559 individuals of diverse ancestry, including 10,159 T2D cases, and test their association with T2D-related vascular outcomes. Cluster-specific partitioned GRS are more strongly associated with coronary artery disease and end-stage diabetic nephropathy than an overall T2D GRS across ancestry groups, highlighting the importance of obesity-related processes in the development of vascular outcomes. Our findings demonstrate the value of integrating multi-ancestry GWAS with single-cell epigenomics to disentangle the aetiological heterogeneity driving the development and progression of T2D, which may offer a route to optimise global access to genetically-informed diabetes care.

Unraveling Neuro-Proteogenomic Landscape and Therapeutic Implications for Human Behaviors and Psychiatric Disorders

Understanding the genetic basis of neuro-related proteins is essential for dissecting the molecular basis of human behavioral traits and the disease etiology of neuropsychiatric disorders. Here, the SCALLOP Consortium conducted a genome-wide association meta-analysis of over 12,500 individuals for 184 neuro-related proteins in human plasma. The analysis identified 117 cis-regulatory protein quantitative trait loci (cis-pQTL) and 166 trans-pQTL. The mapped pQTL capture on average 50% of each protein's heritability. Mendelian randomization analyses revealed multiple proteins showing potential causal effects on neuro-related traits such as sleeping, smoking, feelings, alcohol intake, mental health, and psychiatric disorders. Integrating with established drug information, we validated 13 out of 13 matched combinations of protein targets and diseases or side effects with available drugs, while suggesting hundreds of re-purposing and new therapeutic targets. This consortium effort provides a large-scale proteogenomic resource for biomedical research on human behaviors and other neuro-related phenotypes.

Genetic regulation of fetal hemoglobin across global populations

Human genetic variation has enabled the identification of several key regulators of fetal-to-adult hemoglobin switching, including BCL11A, resulting in therapeutic advances. However, despite the progress made, limited further insights have been obtained to provide a fuller accounting of how genetic variation contributes to the global mechanisms of fetal hemoglobin (HbF) gene regulation. Here, we have conducted a multi-ancestry genome-wide association study of 28,279 individuals from several cohorts spanning 5 continents to define the architecture of human genetic variation impacting HbF. We have identified a total of 178 conditionally independent genome-wide significant or suggestive variants across 14 genomic windows. Importantly, these new data enable us to better define the mechanisms by which HbF switching occurs in vivo. We conduct targeted perturbations to define BACH2 as a new genetically-nominated regulator of hemoglobin switching. We define putative causal variants and underlying mechanisms at the well-studied BCL11A and HBS1L-MYB loci, illuminating the complex variant-driven regulation present at these loci. We additionally show how rare large-effect deletions in the HBB locus can interact with polygenic variation to influence HbF levels. Our study paves the way for the next generation of therapies to more effectively induce HbF in sickle cell disease and β-thalassemia.

The common IL1A single nucleotide polymorphism rs17561 is a hypomorphic mutation that significantly reduces interleukin-1α release from human blood cells

Interleukin-1 alpha (IL-1α) is a powerful cytokine that drives inflammation and modulates adaptive immunity. Due to these powerful effects, IL-1α is controlled at multiple levels from transcription to cleavage and release from the cell. Genome-wide association studies can identify loci that drive important diseases, although often the functional effect of the variant on phenotype remains unknown or small, with most risk variants in non-coding regions. We find that the common variant rs17561 changes a conserved amino acid in the central region of IL-1α linking the pro piece to the cytokine domain. Using a recall-by-genotype study and whole blood stimulation, we find that minor allele homozygotes release ~50% less IL-1α than the major allele, with IL-1β release equivalent. IL-1α transcript level was identical between groups, implying a post-transcriptional effect, whilst cleavage of recombinant pro-IL-1α by multiple proteases was also equivalent for both forms. Importantly, transfected macrophages also release less minor allele IL-1α upon inflammasome activation, revealing that reduced secretion is directly caused by the missense amino acid substitution and more minor allele IL-1α was retained within the cell. Thus, rs17561 represents a very common hypomorphic mutation in IL-1α. We believe this novel data will be important for determining the potential contribution of IL-1α to disease and/or physiological processes, for example, by Mendelian randomisation, and may aid patient stratification when considering anti-IL-1 therapies.

Genetic mechanisms of 184 neuro-related proteins in human plasma

Understanding the genetic basis of neuro-related proteins is essential for dissecting the disease etiology of neuropsychiatric disorders and other complex traits and diseases. Here, the SCALLOP Consortium conducted a genome-wide association meta-analysis of over 12,500 individuals for 184 neuro-reiated proteins in human plasma. The analysis identified 117 cis-regulatory protein quantitative trait loci (cis-pQTL) and 166 trans-pQTL. The mapped pQTL capture on average 50% of each protein's heritability. Mendelian randomization analyses revealed multiple proteins showing potential causal effects on neuro-reiated traits as well as complex diseases such as hypertension, high cholesterol, immune-related disorders, and psychiatric disorders. Integrating with established drug information, we validated 13 combinations of protein targets and diseases or side effects with available drugs, while suggesting hundreds of re-purposing and new therapeutic targets for diseases and comorbidities. This consortium effort provides a large-scale proteogenomic resource for biomedical research.

Predictive Utility of a Coronary Artery Disease Polygenic Risk Score in Primary Prevention

Importance

The clinical utility of polygenic risk scores (PRS) for coronary artery disease (CAD) has not yet been established.

Objective

To investigate the ability of a CAD PRS to potentially guide statin initiation in primary prevention after accounting for age and clinical risk.

Design, Setting, and Participants

This was a longitudinal cohort study with enrollment starting on January 1, 2006, and ending on December 31, 2010, with data updated to mid-2021, using data from the UK Biobank, a long-term population study of UK citizens. A replication analysis was performed in Biobank Japan. The analysis included all patients without a history of CAD and who were not taking lipid-lowering therapy. Data were analyzed from January 1 to June 30, 2022.

Exposures

Polygenic risk for CAD was defined as low (bottom 20%), intermediate, and high (top 20%) using a CAD PRS including 241 genome-wide significant single-nucleotide variations (SNVs). The pooled cohort equations were used to estimate 10-year atherosclerotic cardiovascular disease (ASCVD) risk and classify individuals as low (<5%), borderline (5-<7.5%), intermediate (7.5-<20%), or high risk (≥20%).

Main Outcomes and Measures

Myocardial infarction (MI) and ASCVD events (defined as incident clinical CAD [including MI], stroke, or CV death).

Results

A total of 330 201 patients (median [IQR] age, 57 [40-74] years; 189 107 female individuals [57%]) were included from the UK Biobank. Over the 10-year follow-up, 4454 individuals had an MI. The CAD PRS was significantly associated with the risk of MI in all age groups but had significantly stronger risk prediction at younger ages (age <50 years: hazard ratio [HR] per 1 SD of PRS, 1.72; 95% CI, 1.56-1.89; age 50-60 years: HR, 1.46; 95% CI, 1.38-1.53; age >60 years: HR, 1.42; 95% CI, 1.37-1.48; P for interaction <.001). In patients younger than 50 years, those with high PRS had a 3- to 4-fold increased associated risk of MI compared with those in the low PRS category. A significant interaction between CAD PRS and age was replicated in Biobank Japan. When CAD PRS testing was added to the clinical ASCVD risk score in individuals younger than 50 years, 591 of 4373 patients (20%) with borderline risk were risk stratified into intermediate risk, warranting initiation of statin therapy and 3198 of 7477 patients (20%) with both borderline or intermediate risk were stratified as low risk, thus not warranting therapy.

Conclusions and Relevance

Results of this cohort study suggest that the predictive ability of a CAD PRS was greater in younger individuals and can be used to better identify patients with borderline and intermediate clinical risk who should initiate statin therapy.

Quality control and removal of technical variation of NMR metabolic biomarker data in ~120,000 UK Biobank participants

Metabolic biomarker data quantified by nuclear magnetic resonance (NMR) spectroscopy in approximately 121,000 UK Biobank participants has recently been released as a community resource, comprising absolute concentrations and ratios of 249 circulating metabolites, lipids, and lipoprotein sub-fractions. Here we identify and characterise additional sources of unwanted technical variation influencing individual biomarkers in the data available to download from UK Biobank. These included sample preparation time, shipping plate well, spectrometer batch effects, drift over time within spectrometer, and outlier shipping plates. We developed a procedure for removing this unwanted technical variation, and demonstrate that it increases signal for genetic and epidemiological studies of the NMR metabolic biomarker data in UK Biobank. We subsequently developed an R package, ukbnmr, which we make available to the wider research community to enhance the utility of the UK Biobank NMR metabolic biomarker data and to facilitate rapid analysis.

Trans-ethnic Polygenic Risk Scores for Body Mass Index: An International Hundred K+ Cohorts Consortium Study

Background

While polygenic risk scores hold significant promise in estimating an individual's risk of developing a complex trait such as obesity, their application in the clinic has, to date, been limited by a lack of data from non-European populations. As a collaboration model of the International Hundred K+ Cohorts Consortium (IHCC), we endeavored to develop a globally applicable trans-ethnic PRS for body mass index (BMI) through this relatively new international effort.

Methods

The PRS model was developed trained and tested at the Center for Applied Genomics (CAG) of The Children's Hospital of Philadelphia (CHOP) based on a BMI meta-analysis from the GIANT consortium. The validated PRS models were subsequently disseminated to the participating sites. Scores were generated by each site locally on their cohorts and summary statistics returned to CAG for final analysis.

Results

We show that in the absence of a well powered trans-ethnic GWAS from which to derive SNPs and effect estimates, trans-ethnic scores can be generated from European ancestry GWAS using Bayesian approaches such as LDpred to adjust the summary statistics using trans-ethnic linkage disequilibrium reference panels. The ported trans-ethnic scores outperform population specific-PRS across all non-European ancestry populations investigated including East Asians and three-way admixed Brazilian cohort.

Conclusions

Widespread use of PRS in the clinic is hampered by a lack of genotyping data in individuals of non-European ancestry for the vast majority of traits. Here we show that for a truly polygenic trait such as BMI adjusting the summary statistics of a well powered European ancestry study using trans-ethnic LD reference results in a score that is predictive across a range of ancestries including East Asians and three-way admixed Brazilians.

Fibroblast Growth Factor-23 and Risk of Cardiovascular Diseases: A Mendelian Randomization Study

BACKGROUND

Fibroblast growth factor-23 (FGF-23) is associated with a range of cardiovascular and noncardiovascular diseases in conventional epidemiological studies, but substantial residual confounding may exist. Mendelian randomization approaches can help control for such confounding.

METHODS

SCALLOP Consortium data of 19,195 participants were used to generate an FGF-23 genetic score. Data from 337,448 UK Biobank participants were used to estimate associations between higher genetically predicted FGF-23 concentration and the odds of any atherosclerotic cardiovascular disease (n=26,266 events), nonatherosclerotic cardiovascular disease (n=12,652), and noncardiovascular diseases previously linked to FGF-23. Measurements of carotid intima-media thickness and left ventricular mass were available in a subset. Associations with cardiovascular outcomes were also tested in three large case-control consortia: CARDIOGRAMplusC4D (coronary artery disease, n=181,249 cases), MEGASTROKE (stroke, n=34,217), and HERMES (heart failure, n=47,309).

RESULTS

We identified 34 independent variants for circulating FGF-23, which formed a validated genetic score. There were no associations between genetically predicted FGF-23 and any of the cardiovascular or noncardiovascular outcomes. In UK Biobank, the odds ratio (OR) for any atherosclerotic cardiovascular disease per 1-SD higher genetically predicted logFGF-23 was 1.03 (95% confidence interval [95% CI], 0.98 to 1.08), and for any nonatherosclerotic cardiovascular disease, it was 1.01 (95% CI, 0.94 to 1.09). The ORs in the case-control consortia were 1.00 (95% CI, 0.97 to 1.03) for coronary artery disease, 1.01 (95% CI, 0.95 to 1.07) for stroke, and 1.00 (95% CI, 0.95 to 1.05) for heart failure. In those with imaging, logFGF-23 was not associated with carotid or cardiac abnormalities.

CONCLUSIONS

Genetically predicted FGF-23 levels are not associated with atherosclerotic and nonatherosclerotic cardiovascular diseases, suggesting no important causal link.

PODCAST

This article contains a podcast at https://dts.podtrac.com/redirect.mp3/www.asn-online.org/media/podcast/CJASN/2023_01_10_CJN05080422.mp3.

Discovery and systematic characterization of risk variants and genes for coronary artery disease in over a million participants

The discovery of genetic loci associated with complex diseases has outpaced the elucidation of mechanisms of disease pathogenesis. Here we conducted a genome-wide association study (GWAS) for coronary artery disease (CAD) comprising 181,522 cases among 1,165,690 participants of predominantly European ancestry. We detected 241 associations, including 30 new loci. Cross-ancestry meta-analysis with a Japanese GWAS yielded 38 additional new loci. We prioritized likely causal variants using functionally informed fine-mapping, yielding 42 associations with less than five variants in the 95% credible set. Similarity-based clustering suggested roles for early developmental processes, cell cycle signaling and vascular cell migration and proliferation in the pathogenesis of CAD. We prioritized 220 candidate causal genes, combining eight complementary approaches, including 123 supported by three or more approaches. Using CRISPR-Cas9, we experimentally validated the effect of an enhancer in MYO9B, which appears to mediate CAD risk by regulating vascular cell motility. Our analysis identifies and systematically characterizes >250 risk loci for CAD to inform experimental interrogation of putative causal mechanisms for CAD.

Gene Sequencing Identifies Perturbation in Nitric Oxide Signaling as a Nonlipid Molecular Subtype of Coronary Artery Disease

BACKGROUND

A key goal of precision medicine is to disaggregate common, complex diseases into discrete molecular subtypes. Rare coding variants in the low-density lipoprotein receptor gene (LDLR) are identified in 1% to 2% of coronary artery disease (CAD) patients, defining a molecular subtype with risk driven by hypercholesterolemia.

METHODS

To search for additional subtypes, we compared the frequency of rare, predicted loss-of-function and damaging missense variants aggregated within a given gene in 41 081 CAD cases versus 217 115 controls.

RESULTS

Rare variants in LDLR were most strongly associated with CAD, present in 1% of cases and associated with 4.4-fold increased CAD risk. A second subtype was characterized by variants in endothelial nitric oxide synthase gene (NOS3), a key enzyme regulating vascular tone, endothelial function, and platelet aggregation. A rare predicted loss-of-function or damaging missense variants in NOS3 was present in 0.6% of cases and associated with 2.42-fold increased risk of CAD (95% CI, 1.80-3.26; P=5.50×10-9). These variants were associated with higher systolic blood pressure (+3.25 mm Hg; [95% CI, 1.86-4.65]; P=5.00×10-6) and increased risk of hypertension (adjusted odds ratio 1.31; [95% CI, 1.14-1.51]; P=2.00×10-4) but not circulating cholesterol concentrations, suggesting that, beyond lipid pathways, nitric oxide synthesis is a key nonlipid driver of CAD risk.

CONCLUSIONS

Beyond LDLR, we identified an additional nonlipid molecular subtype of CAD characterized by rare variants in the NOS3 gene.

Genetically personalised organ-specific metabolic models in health and disease

Understanding how genetic variants influence disease risk and complex traits (variant-to-function) is one of the major challenges in human genetics. Here we present a model-driven framework to leverage human genome-scale metabolic networks to define how genetic variants affect biochemical reaction fluxes across major human tissues, including skeletal muscle, adipose, liver, brain and heart. As proof of concept, we build personalised organ-specific metabolic flux models for 524,615 individuals of the INTERVAL and UK Biobank cohorts and perform a fluxome-wide association study (FWAS) to identify 4312 associations between personalised flux values and the concentration of metabolites in blood. Furthermore, we apply FWAS to identify 92 metabolic fluxes associated with the risk of developing coronary artery disease, many of which are linked to processes previously described to play in role in the disease. Our work demonstrates that genetically personalised metabolic models can elucidate the downstream effects of genetic variants on biochemical reactions involved in common human diseases.

Mild-to-Moderate Kidney Dysfunction and Cardiovascular Disease: Observational and Mendelian Randomization Analyses

BACKGROUND

End-stage renal disease is associated with a high risk of cardiovascular events. It is unknown, however, whether mild-to-moderate kidney dysfunction is causally related to coronary heart disease (CHD) and stroke.

METHODS

Observational analyses were conducted using individual-level data from 4 population data sources (Emerging Risk Factors Collaboration, EPIC-CVD [European Prospective Investigation into Cancer and Nutrition-Cardiovascular Disease Study], Million Veteran Program, and UK Biobank), comprising 648 135 participants with no history of cardiovascular disease or diabetes at baseline, yielding 42 858 and 15 693 incident CHD and stroke events, respectively, during 6.8 million person-years of follow-up. Using a genetic risk score of 218 variants for estimated glomerular filtration rate (eGFR), we conducted Mendelian randomization analyses involving 413 718 participants (25 917 CHD and 8622 strokes) in EPIC-CVD, Million Veteran Program, and UK Biobank.

RESULTS

There were U-shaped observational associations of creatinine-based eGFR with CHD and stroke, with higher risk in participants with eGFR values <60 or >105 mL·min-1·1.73 m-2, compared with those with eGFR between 60 and 105 mL·min-1·1.73 m-2. Mendelian randomization analyses for CHD showed an association among participants with eGFR <60 mL·min-1·1.73 m-2, with a 14% (95% CI, 3%-27%) higher CHD risk per 5 mL·min-1·1.73 m-2 lower genetically predicted eGFR, but not for those with eGFR >105 mL·min-1·1.73 m-2. Results were not materially different after adjustment for factors associated with the eGFR genetic risk score, such as lipoprotein(a), triglycerides, hemoglobin A1c, and blood pressure. Mendelian randomization results for stroke were nonsignificant but broadly similar to those for CHD.

CONCLUSIONS

In people without manifest cardiovascular disease or diabetes, mild-to-moderate kidney dysfunction is causally related to risk of CHD, highlighting the potential value of preventive approaches that preserve and modulate kidney function.

Rare and common genetic determinants of metabolic individuality and their effects on human health

Garrod's concept of 'chemical individuality' has contributed to comprehension of the molecular origins of human diseases. Untargeted high-throughput metabolomic technologies provide an in-depth snapshot of human metabolism at scale. We studied the genetic architecture of the human plasma metabolome using 913 metabolites assayed in 19,994 individuals and identified 2,599 variant-metabolite associations (P < 1.25 × 10^-11) within 330 genomic regions, with rare variants (minor allele frequency ≤ 1%) explaining 9.4% of associations. Jointly modeling metabolites in each region, we identified 423 regional, co-regulated, variant-metabolite clusters called genetically influenced metabotypes. We assigned causal genes for 62.4% of these genetically influenced metabotypes, providing new insights into fundamental metabolite physiology and clinical relevance, including metabolite-guided discovery of potential adverse drug effects (DPYD and SRD5A2). We show strong enrichment of inborn errors of metabolism-causing genes, with examples of metabolite associations and clinical phenotypes of non-pathogenic variant carriers matching characteristics of the inborn errors of metabolism. Systematic, phenotypic follow-up of metabolite-specific genetic scores revealed multiple potential etiological relationships.

Systematic Mendelian randomization using the human plasma proteome to discover potential therapeutic targets for stroke

Stroke is the second leading cause of death with substantial unmet therapeutic needs. To identify potential stroke therapeutic targets, we estimate the causal effects of 308 plasma proteins on stroke outcomes in a two-sample Mendelian randomization framework and assess mediation effects by stroke risk factors. We find associations between genetically predicted plasma levels of six proteins and stroke (P ≤ 1.62 × 10-4). The genetic associations with stroke colocalize (Posterior Probability >0.7) with the genetic associations of four proteins (TFPI, TMPRSS5, CD6, CD40). Mendelian randomization supports atrial fibrillation, body mass index, smoking, blood pressure, white matter hyperintensities and type 2 diabetes as stroke risk factors (P ≤ 0.0071). Body mass index, white matter hyperintensity and atrial fibrillation appear to mediate the TFPI, IL6RA, TMPRSS5 associations with stroke. Furthermore, thirty-six proteins are associated with one or more of these risk factors using Mendelian randomization. Our results highlight causal pathways and potential therapeutic targets for stroke.

Physical activity attenuates but does not eliminate coronary heart disease risk amongst adults with risk factors: EPIC-CVD case-cohort study

AIMS

This study aimed to evaluate the association between physical activity and the incidence of coronary heart disease (CHD) in individuals with and without CHD risk factors.

METHODS AND RESULTS

EPIC-CVD is a case-cohort study of 29 333 participants that included 13 582 incident CHD cases and a randomly selected sub-cohort nested within the European Prospective Investigation into Cancer and Nutrition (EPIC) cohort. Self-reported physical activity was summarized using the Cambridge physical activity index (inactive, moderately inactive, moderately active, and active). Participants were categorized into sub-groups based on the presence or the absence of the following risk factors: obesity (body mass index ≥30 kg/m2), hypercholesterolaemia (total cholesterol ≥6.2 mmol/L), history of diabetes, hypertension (self-reported or ≥140/90 mmHg), and current smoking. Prentice-weighted Cox regression was used to assess the association between physical activity and incident CHD events (non-fatal and fatal).Compared to inactive participants without the respective CHD risk factor (referent), excess CHD risk was highest in physically inactive and lowest in moderately active participants with CHD risk factors. Corresponding excess CHD risk estimates amongst those with obesity were 47% [95% confidence interval (CI) 32-64%] and 21% (95%CI 2-44%), with hypercholesterolaemia were 80% (95%CI 55-108%) and 48% (95%CI 22-81%), with hypertension were 80% (95%CI 65-96%) and 49% (95%CI 28-74%), with diabetes were 142% (95%CI 63-260%), and 100% (95%CI 32-204%), and amongst smokers were 152% (95%CI 122-186%) and 109% (95%CI 74-150%).

CONCLUSIONS

In people with CHD risk factors, moderate physical activity, equivalent to 40 mins of walking per day, attenuates but does not completely offset CHD risk.

Neurocognitive trajectory and proteomic signature of inherited risk for Alzheimer's disease

For Alzheimer's disease-a leading cause of dementia and global morbidity-improved identification of presymptomatic high-risk individuals and identification of new circulating biomarkers are key public health needs. Here, we tested the hypothesis that a polygenic predictor of risk for Alzheimer's disease would identify a subset of the population with increased risk of clinically diagnosed dementia, subclinical neurocognitive dysfunction, and a differing circulating proteomic profile. Using summary association statistics from a recent genome-wide association study, we first developed a polygenic predictor of Alzheimer's disease comprised of 7.1 million common DNA variants. We noted a 7.3-fold (95% CI 4.8 to 11.0; p < 0.001) gradient in risk across deciles of the score among 288,289 middle-aged participants of the UK Biobank study. In cross-sectional analyses stratified by age, minimal differences in risk of Alzheimer's disease and performance on a digit recall test were present according to polygenic score decile at age 50 years, but significant gradients emerged by age 65. Similarly, among 30,541 participants of the Mass General Brigham Biobank, we again noted no significant differences in Alzheimer's disease diagnosis at younger ages across deciles of the score, but for those over 65 years we noted an odds ratio of 2.0 (95% CI 1.3 to 3.2; p = 0.002) in the top versus bottom decile of the polygenic score. To understand the proteomic signature of inherited risk, we performed aptamer-based profiling in 636 blood donors (mean age 43 years) with very high or low polygenic scores. In addition to the well-known apolipoprotein E biomarker, this analysis identified 27 additional proteins, several of which have known roles related to disease pathogenesis. Differences in protein concentrations were consistent even among the youngest subset of blood donors (mean age 33 years). Of these 28 proteins, 7 of the 8 proteins with concentrations available were similarly associated with the polygenic score in participants of the Multi-Ethnic Study of Atherosclerosis. These data highlight the potential for a DNA-based score to identify high-risk individuals during the prolonged presymptomatic phase of Alzheimer's disease and to enable biomarker discovery based on profiling of young individuals in the extremes of the score distribution.

Higher body mass index raises immature platelet count: potential contribution to obesity-related thrombosis

Higher body mass index (BMI) is a risk factor for thrombosis. Platelets are essential for hemostasis but contribute to thrombosis when activated pathologically. We hypothesized that higher BMI leads to changes in platelet characteristics, thereby increasing thrombotic risk. The effect of BMI on platelet traits (measured by Sysmex) was explored in 33 388 UK blood donors (INTERVAL study). Linear regression showed that higher BMI was positively associated with greater plateletcrit (PCT), platelet count (PLT), immature platelet count (IPC), and side fluorescence (SFL, a measure of mRNA content used to derive IPC). Mendelian randomization (MR), applied to estimate a causal effect with BMI proxied by a genetic risk score, provided causal estimates for a positive effect of BMI on both SFL and IPC, but there was little evidence for a causal effect of BMI on PCT or PLT. Follow-up analyses explored the functional relevance of platelet characteristics in a pre-operative cardiac cohort (COPTIC). Linear regression provided observational evidence for a positive association between IPC and agonist-induced whole blood platelet aggregation. Results indicate that higher BMI raises the number of immature platelets, which is associated with greater whole blood platelet aggregation in a cardiac cohort. Higher IPC could therefore contribute to obesity-related thrombosis.

Integrated analyses of growth differentiation factor-15 concentration and cardiometabolic diseases in humans

Growth differentiation factor-15 (GDF15) is a stress response cytokine that is elevated in several cardiometabolic diseases and has attracted interest as a potential therapeutic target. To further explore the association of GDF15 with human disease, we conducted a broad study into the phenotypic and genetic correlates of GDF15 concentration in up to 14,099 individuals. Assessment of 772 traits across 6610 participants in FINRISK identified associations of GDF15 concentration with a range of phenotypes including all-cause mortality, cardiometabolic disease, respiratory diseases and psychiatric disorders, as well as inflammatory markers. A meta-analysis of genome-wide association studies (GWAS) of GDF15 concentration across three different assay platforms (n=14,099) confirmed significant heterogeneity due to a common missense variant (rs1058587; p.H202D) in GDF15, potentially due to epitope-binding artefacts. After conditioning on rs1058587, statistical fine mapping identified four independent putative causal signals at the locus. Mendelian randomisation (MR) analysis found evidence of a causal relationship between GDF15 concentration and high-density lipoprotein (HDL) but not body mass index (BMI). Using reverse MR, we identified a potential causal association of BMI on GDF15 (IVW pFDR = 0.0040). Taken together, our data derived from human population cohorts do not support a role for moderately elevated GDF15 concentrations as a causal factor in human cardiometabolic disease but support its role as a biomarker of metabolic stress.

Flashfm-ivis: interactive visualization for fine-mapping of multiple quantitative traits

SUMMARY

flashfm-ivis provides a suite of interactive visualization plots to view potential causal genetic variants that underlie associations that are shared or distinct between multiple quantitative traits and compares results between single- and multi-trait fine-mapping. Unique features include network diagrams that show joint effects between variants for each trait and regional association plots that integrate fine-mapping results, all with user-controlled zoom features for an interactive exploration of potential causal variants across traits.

AVAILABILITY AND IMPLEMENTATION

flashfm-ivis is an open-source software under the MIT license. It is available as an interactive web-based tool (http://shiny.mrc-bsu.cam.ac.uk/apps/flashfm-ivis/) and as an R package. Code and documentation are available at https://github.com/fz-cambridge/flashfm-ivis and https://zenodo.org/record/6376244#.YjnarC-l2X0. Additional features can be downloaded as standalone R libraries to encourage reuse.

SUPPLEMENTARY INFORMATION

Supplementary information are available at Bioinformatics online.

Genetically Determined Reproductive Aging and Coronary Heart Disease: A Bidirectional 2-sample Mendelian Randomization

BACKGROUND

Accelerated reproductive aging, in women indicated by early natural menopause, is associated with increased coronary heart disease (CHD) risk in observational studies. Conversely, an adverse CHD risk profile has been suggested to accelerate menopause.

OBJECTIVES

To study the direction and evidence for causality of the relationship between reproductive aging and (non-)fatal CHD and CHD risk factors in a bidirectional Mendelian randomization (MR) approach, using age at natural menopause (ANM) genetic variants as a measure for genetically determined reproductive aging in women. We also studied the association of these variants with CHD risk (factors) in men.

DESIGN

Two-sample MR, using both cohort data as well as summary statistics, with 4 methods: simple and weighted median-based, standard inverse-variance weighted (IVW) regression, and MR-Egger regression.

PARTICIPANTS

Data from EPIC-CVD and summary statistics from UK Biobank and publicly available genome-wide association studies were pooled for the different analyses.

MAIN OUTCOME MEASURES

CHD, CHD risk factors, and ANM.

RESULTS

Across different methods of MR, no association was found between genetically determined reproductive aging and CHD risk in women (relative risk estimateIVW = 0.99; 95% confidence interval (CI), 0.97-1.01), or any of the CHD risk factors. Similarly, no associations were found in men. Neither did the reversed analyses show evidence for an association between CHD (risk factors) and reproductive aging.

CONCLUSION

Genetically determined reproductive aging is not causally associated with CHD risk (factors) in women, nor were the genetic variants associated in men. We found no evidence for a reverse association in a combined sample of women and men.

Whole-exome sequencing identifies rare genetic variants associated with human plasma metabolites

Metabolite levels measured in the human population are endophenotypes for biological processes. We combined sequencing data for 3,924 (whole-exome sequencing, WES, discovery) and 2,805 (whole-genome sequencing, WGS, replication) donors from a prospective cohort of blood donors in England. We used multiple approaches to select and aggregate rare genetic variants (minor allele frequency [MAF] < 0.1%) in protein-coding regions and tested their associations with 995 metabolites measured in plasma by using ultra-high-performance liquid chromatography-tandem mass spectrometry. We identified 40 novel associations implicating rare coding variants (27 genes and 38 metabolites), of which 28 (15 genes and 28 metabolites) were replicated. We developed algorithms to prioritize putative driver variants at each locus and used mediation and Mendelian randomization analyses to test directionality at associations of metabolite and protein levels at the ACY1 locus. Overall, 66% of reported associations implicate gene targets of approved drugs or bioactive drug-like compounds, contributing to drug targets' validating efforts.

Association of shorter leucocyte telomere length with risk of frailty

BACKGROUND

Frailty is a multidimensional syndrome of decline that affects multiple systems and predisposes to adverse health outcomes. Although chronological age is the major risk factor, inter-individual variation in risk is not fully understood. Leucocyte telomere length (LTL), a proposed marker of biological age, has been associated with risk of many diseases. We sought to determine whether LTL is associated with risk of frailty.

METHODS

We utilized cross-sectional data from 441 781 UK Biobank participants (aged 40-69 years), with complete data on frailty indicators and LTL. Frailty was defined as the presence of at least three of five indicators: weaker grip strength, slower walking pace, weight loss in the past year, lower physical activity, and exhaustion in the past 2 weeks. LTL was measured using a validated qPCR method and reported as a ratio of the telomere repeat number (T) to a single-copy gene (S) (T/S ratio). Association of LTL with frailty was evaluated using adjusted (chronological age, sex, deprivation, smoking, alcohol intake, body mass index, and multimorbidity) multinomial and ordinal regression models, and results are presented as relative risk (RRR) or odds ratios (OR), respectively, alongside the 95% confidence interval (CI). Mendelian randomization (MR), using 131 genetic variants associated with LTL, was used to assess if the association of LTL with frailty was causal.

RESULTS

Frail participants (4.6%) were older (median age difference (95% CI): 3 (2.5; 3.5) years, P = 2.73 × 10-33 ), more likely to be female (61%, P = 1.97 × 10-129 ), and had shorter LTL (-0.13SD vs. 0.03SD, P = 5.43 × 10-111 ) than non-frail. In adjusted analyses, both age and LTL were associated with frailty (RRR = 1.03 (95% CI: 1.02; 1.04) per year of older chronological age, P = 3.99 × 10-12 ; 1.10 (1.08; 1.11) per SD shorter LTL, P = 1.46 × 10-30 ). Within each age group (40-49, 50-59, 60-69 years), the prevalence of frailty was about 33% higher in participants with shorter (-2SD) versus longer telomeres (+2SD). MR analysis showed an association of LTL with frailty that was directionally consistent with the observational association, but not statistically significant (MR-Median: OR (95% CI): 1.08 (0.98; 1.19) per SD shorter LTL, P = 0.13).

CONCLUSIONS

Inter-individual variation in LTL is associated with the risk of frailty independently of chronological age and other risk factors. Our findings provide evidence for an additional biological determinant of frailty.

Genetic Landscape of the ACE2 Coronavirus Receptor

BACKGROUND

SARS-CoV-2, the causal agent of COVID-19, enters human cells using the ACE2 (angiotensin-converting enzyme 2) protein as a receptor. ACE2 is thus key to the infection and treatment of the coronavirus. ACE2 is highly expressed in the heart and respiratory and gastrointestinal tracts, playing important regulatory roles in the cardiovascular and other biological systems. However, the genetic basis of the ACE2 protein levels is not well understood.

METHODS

We have conducted the largest genome-wide association meta-analysis of plasma ACE2 levels in >28 000 individuals of the SCALLOP Consortium (Systematic and Combined Analysis of Olink Proteins). We summarize the cross-sectional epidemiological correlates of circulating ACE2. Using the summary statistics-based high-definition likelihood method, we estimate relevant genetic correlations with cardiometabolic phenotypes, COVID-19, and other human complex traits and diseases. We perform causal inference of soluble ACE2 on vascular disease outcomes and COVID-19 severity using mendelian randomization. We also perform in silico functional analysis by integrating with other types of omics data.

RESULTS

We identified 10 loci, including 8 novel, capturing 30% of the heritability of the protein. We detected that plasma ACE2 was genetically correlated with vascular diseases, severe COVID-19, and a wide range of human complex diseases and medications. An X-chromosome cis-protein quantitative trait loci-based mendelian randomization analysis suggested a causal effect of elevated ACE2 levels on COVID-19 severity (odds ratio, 1.63 [95% CI, 1.10-2.42]; P=0.01), hospitalization (odds ratio, 1.52 [95% CI, 1.05-2.21]; P=0.03), and infection (odds ratio, 1.60 [95% CI, 1.08-2.37]; P=0.02). Tissue- and cell type-specific transcriptomic and epigenomic analysis revealed that the ACE2 regulatory variants were enriched for DNA methylation sites in blood immune cells.

CONCLUSIONS

Human plasma ACE2 shares a genetic basis with cardiovascular disease, COVID-19, and other related diseases. The genetic architecture of the ACE2 protein is mapped, providing a useful resource for further biological and clinical studies on this coronavirus receptor.

Modifiable traits, healthy behaviours, and leukocyte telomere length: a population-based study in UK Biobank

Background

Telomere length is associated with risk of several age-related diseases and cancers. We aimed to investigate the extent to which telomere length might be modifiable through lifestyle and behaviour, and whether such modification has any clinical consequences.

Methods

In this population-based study, we included participants from UK Biobank who had leukocyte telomere length (LTL) measurement, ethnicity, and white blood cell count data. We investigated associations of LTL with 117 potentially modifiable traits, as well as two indices of healthy behaviours incorporating between them smoking, physical activity, diet, maintenance of a healthy bodyweight, and alcohol intake, using both available and imputed data. To help interpretation, associations were summarised as the number of equivalent years of age-related change in LTL by dividing the trait β coefficients with the age β coefficient. We used mendelian randomisation to test causality of selected associations. We investigated whether the associations of LTL with 22 diseases were modified by the number of healthy behaviours and the extent to which the associations of more healthy behaviours with greater life expectancy and lower risk of coronary artery disease might be mediated through LTL.

Findings

422 797 participants were available for the analysis (227 620 [53·8%] were women and 400 036 [94·6%] were White). 71 traits showed significant (p<4·27 × 10-4) associations with LTL but most were modest, equivalent to less than 1 year of age-related change in LTL. In multivariable analyses of 17 traits with stronger associations (equivalent to ≥2 years of age-related change in LTL), oily fish intake, educational attainment, and general health status retained a significant association of this magnitude, with walking pace and current smoking being additionally significant at this level of association in the imputed models. Mendelian randomisation analysis suggested that educational attainment and smoking behaviour causally affect LTL. Both indices of healthy behaviour were positively and linearly associated with LTL, with those with the most healthy behaviours having longer LTL equivalent to about 3·5 years of age-related change in LTL than those with the least heathy behaviours (p<0·001). However, healthy behaviours explained less than 0·2% of the total variation in LTL and did not significantly modify the association of LTL with risk of any of the diseases studied. Neither the association of more healthy behaviours on greater life expectancy or lower risk of coronary artery disease were substantially mediated through LTL.

Interpretation

Although several potentially modifiable traits and healthy behaviours have a quantifiable association with LTL, at least some of which are likely to be causal, these effects are not of a sufficient magnitude to substantially alter the association between LTL and various diseases or life expectancy. Attempts to change telomere length through lifestyle or behavioural changes might not confer substantial clinical benefit.

Funding

UK Medical Research Council, UK Biotechnology and Biological Sciences Research Council, and British Heart Foundation.