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.

Inherited polygenic effects on common hematological traits influence clonal selection on JAK2V617F and the development of myeloproliferative neoplasms

Myeloproliferative neoplasms (MPNs) are chronic cancers characterized by overproduction of mature blood cells. Their causative somatic mutations, for example, JAK2V617F, are common in the population, yet only a minority of carriers develop MPN. Here we show that the inherited polygenic loci that underlie common hematological traits influence JAK2V617F clonal expansion. We identify polygenic risk scores (PGSs) for monocyte count and plateletcrit as new risk factors for JAK2V617F positivity. PGSs for several hematological traits influenced the risk of different MPN subtypes, with low PGSs for two platelet traits also showing protective effects in JAK2V617F carriers, making them two to three times less likely to have essential thrombocythemia than carriers with high PGSs. We observed that extreme hematological PGSs may contribute to an MPN diagnosis in the absence of somatic driver mutations. Our study showcases how polygenic backgrounds underlying common hematological traits influence both clonal selection on somatic mutations and the subsequent phenotype of cancer.

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.

The effects of pathogenic and likely pathogenic variants for inherited hemostasis disorders in 140 214 UK Biobank participants

Rare genetic diseases affect millions, and identifying causal DNA variants is essential for patient care. Therefore, it is imperative to estimate the effect of each independent variant and improve their pathogenicity classification. Our study of 140 214 unrelated UK Biobank (UKB) participants found that each of them carries a median of 7 variants previously reported as pathogenic or likely pathogenic. We focused on 967 diagnostic-grade gene (DGG) variants for rare bleeding, thrombotic, and platelet disorders (BTPDs) observed in 12 367 UKB participants. By association analysis, for a subset of these variants, we estimated effect sizes for platelet count and volume, and odds ratios for bleeding and thrombosis. Variants causal of some autosomal recessive platelet disorders revealed phenotypic consequences in carriers. Loss-of-function variants in MPL, which cause chronic amegakaryocytic thrombocytopenia if biallelic, were unexpectedly associated with increased platelet counts in carriers. We also demonstrated that common variants identified by genome-wide association studies (GWAS) for platelet count or thrombosis risk may influence the penetrance of rare variants in BTPD DGGs on their associated hemostasis disorders. Network-propagation analysis applied to an interactome of 18 410 nodes and 571 917 edges showed that GWAS variants with large effect sizes are enriched in DGGs and their first-order interactors. Finally, we illustrate the modifying effect of polygenic scores for platelet count and thrombosis risk on disease severity in participants carrying rare variants in TUBB1 or PROC and PROS1, respectively. Our findings demonstrate the power of association analyses using large population datasets in improving pathogenicity classifications of rare variants.

Genetically predicted androgenic profiles and adverse cardiac markers: a sex-specific Mendelian randomization study

AIMS

Observational evidence suggests associations between sex hormone levels and heart failure (HF). We used sex-specific genetic variants associated with androgenic sex hormone profiles to investigate the causal relevance of androgenic sex hormone profiles on cardiac structure and function and HF using Mendelian randomization (MR).

METHODS AND RESULTS

Sex-specific uncorrelated genome-wide significant (P < 5 × 10-8 ) variants predicting sex hormone-binding globulin (SHBG), total testosterone, and bioavailable testosterone were extracted from summary statistics of genome-wide association study (GWAS) on 425 097 participants in the UK Biobank. Sex-specific gene-outcome association estimates were computed for left ventricular ejection fraction (LVEF), left ventricular end-diastolic and end-systolic volumes (LVEDV and LVESV, respectively), left ventricular stroke volume (LVSV), cardiac index, and cardiac output in 11 528 female and 14 356 male UK Biobank Imaging Study participants and for incident or prevalent HF in an external cohort of 47 309 cases and 930 014 controls. Inverse-variance weighted MR was the primary analysis method. In females, higher genetically predicted bioavailable testosterone was associated with lower LVEDV [β per nmol/L = -0.11 (-0.19 to -0.03), P = 0.006], lower LVESV [β = -0.09 (-0.17 to -0.01), P = 0.022], lower LVSV [β = -0.11 (-0.18 to -0.03), P = 0.005], lower cardiac output [β = -0.08 (-0.16 to 0.00), P = 0.046], and lower cardiac index [β = -0.08 (-0.16 to -0.01), P = 0.034] and a higher risk of HF [odds ratio 1.10 (1.01-1.19), P = 0.026] on external validation analysis in larger scale, sex-adjusted GWAS data. Higher genetically predicted SHBG was associated with higher LVEDV [β per nmol/L = 0.17 (0.08-0.25), P = 2 × 10-4 ], higher LVESV [β = 0.13 (0.05-0.22), P = 0.003], and higher LVSV [β = 0.18 (0.08-0.28), P = 2 × 10-4 ]. In males, higher genetically predicted total and bioavailable testosterone was associated with lower LVESV [β = -0.07 (-0.12 to -0.02), P = 0.007] and LVEF [β = -0.11 (-0.18 to -0.04), P = 0.003], respectively.

CONCLUSIONS

This study supports a causal effect of pro-androgenic sex hormone profiles in females on adverse markers of left ventricular structure and function typically associated with HF with preserved ejection fraction and with HF. There was weaker evidence of association in males.

The relevance of competing risk adjustment in cardiovascular risk prediction models for clinical practice

BACKGROUND

Many models developed for predicting the risk of cardiovascular disease (CVD) are adjusted for the competing risk of non-CVD mortality, which has been suggested to reduce potential overestimation of cumulative incidence in populations where the risk of competing events is high. The objective was to evaluate and illustrate the clinical impact of competing risk adjustment when deriving a CVD prediction model in a high-risk population.

METHODS AND RESULTS

Individuals with established atherosclerotic CVD were included from the Utrecht Cardiovascular Cohort-Secondary Manifestations of ARTerial disease (UCC-SMART). In 8355 individuals, followed for a median of 8.2 years (IQR 4.2-12.5), two similar prediction models for the estimation of 10-year residual CVD risk were derived: with competing risk adjustment using a Fine and Gray model and without competing risk adjustment using a Cox proportional hazards model. On average, predictions were higher from the Cox model. The Cox model predictions overestimated the cumulative incidence [predicted-observed ratio 1.14 (95% CI 1.09-1.20)], which was most apparent in the highest risk quartiles and in older persons. Discrimination of both models was similar. When determining treatment eligibility on thresholds of predicted risks, more individuals would be treated based on the Cox model predictions. If, for example, individuals with a predicted risk > 20% were considered eligible for treatment, 34% of the population would be treated according to the Fine and Gray model predictions and 44% according to the Cox model predictions.

INTERPRETATION

Individual predictions from the model unadjusted for competing risks were higher, reflecting the different interpretations of both models. For models aiming to accurately predict absolute risks, especially in high-risk populations, competing risk adjustment must be considered.

Personalized lifetime prediction of survival and treatment benefit in patients with heart failure with reduced ejection fraction: The LIFE-HF model

AIMS

Although trials have proven the group-level effectiveness of various therapies for heart failure with reduced ejection fraction (HFrEF), important differences in absolute effectiveness exist between individuals. We developed and validated the LIFEtime-perspective for Heart Failure (LIFE-HF) model for the prediction of individual (lifetime) risk and treatment benefit in patients with HFrEF.

METHODS AND RESULTS

Cox proportional hazards functions with age as the time scale were developed in the PARADIGM-HF and ATMOSPHERE trials (n = 15 415). Outcomes were cardiovascular death, heart failure (HF) hospitalization or cardiovascular death, and non-cardiovascular mortality. Predictors were age, sex, New York Heart Association class, prior HF hospitalization, diabetes mellitus, extracardiac vascular disease, systolic blood pressure, left ventricular ejection fraction, N-terminal pro-B-type natriuretic peptide, and glomerular filtration rate. The functions were combined in life-tables to predict individual overall and HF hospitalization-free survival. External validation was performed in the SwedeHF registry, ASIAN-HF registry, and DAPA-HF trial (n = 51 286). Calibration of 2- to 10-year risk was adequate, and c-statistics were 0.65-0.74. An interactive tool was developed combining the model with hazard ratios from trials to allow estimation of an individual's (lifetime) risk and treatment benefit in clinical practice. Applying the tool to the development cohort, combined treatment with a mineralocorticoid receptor antagonist, sodium-glucose cotransporter 2 inhibitor, and angiotensin receptor-neprilysin inhibitor was estimated to afford a median of 2.5 (interquartile range [IQR] 1.7-3.7) and 3.7 (IQR 2.4-5.5) additional years of overall and HF hospitalization-free survival, respectively.

CONCLUSION

The LIFE-HF model enables estimation of lifelong overall and HF hospitalization-free survival, and (lifetime) treatment benefit for individual patients with HFrEF. It could serve as a tool to improve the management of HFrEF by facilitating personalized medicine and shared decision-making.

Improving 10-year cardiovascular risk prediction in apparently healthy people: flexible addition of risk modifiers on top of SCORE2

AIMS

In clinical practice, factors associated with cardiovascular disease (CVD) like albuminuria, education level, or coronary artery calcium (CAC) are often known, but not incorporated in cardiovascular risk prediction models. The aims of the current study were to evaluate a methodology for the flexible addition of risk modifying characteristics on top of SCORE2 and to quantify the added value of several clinically relevant risk modifying characteristics.

METHODS AND RESULTS

Individuals without previous CVD or DM were included from the UK Biobank; Atherosclerosis Risk in Communities (ARIC); Multi-Ethnic Study of Atherosclerosis (MESA); European Prospective Investigation into Cancer, The Netherlands (EPIC-NL); and Heinz Nixdorf Recall (HNR) studies (n = 409 757) in whom 16 166 CVD events and 19 149 non-cardiovascular deaths were observed over exactly 10.0 years of follow-up. The effect of each possible risk modifying characteristic was derived using competing risk-adjusted Fine and Gray models. The risk modifying characteristics were applied to individual predictions with a flexible method using the population prevalence and the subdistribution hazard ratio (SHR) of the relevant predictor. Risk modifying characteristics that increased discrimination most were CAC percentile with 0.0198 [95% confidence interval (CI) 0.0115; 0.0281] and hs-Troponin-T with 0.0100 (95% CI 0.0063; 0.0137). External validation was performed in the Clinical Practice Research Datalink (CPRD) cohort (UK, n = 518 015, 12 675 CVD events). Adjustment of SCORE2-predicted risks with both single and multiple risk modifiers did not negatively affect calibration and led to a modest increase in discrimination [0.740 (95% CI 0.736-0.745) vs. unimproved SCORE2 risk C-index 0.737 (95% CI 0.732-0.741)].

CONCLUSION

The current paper presents a method on how to integrate possible risk modifying characteristics that are not included in existing CVD risk models for the prediction of CVD event risk in apparently healthy people. This flexible methodology improves the accuracy of predicted risks and increases applicability of prediction models for individuals with additional risk known modifiers.

Association of Longer Leukocyte Telomere Length With Cardiac Size, Function, and Heart Failure

Importance

Longer leukocyte telomere length (LTL) is associated with a lower risk of adverse cardiovascular outcomes. The extent to which variation in LTL is associated with intermediary cardiovascular phenotypes is unclear.

Objective

To evaluate the associations between LTL and a diverse set of cardiovascular imaging phenotypes.

Design, Setting, and Participants

This is a population-based cross-sectional study of UK Biobank participants recruited from 2006 to 2010. LTL was measured using a quantitative polymerase chain reaction method. Cardiovascular measurements were derived from cardiovascular magnetic resonance using machine learning. The median (IQR) duration of follow-up was 12.0 (11.3-12.7) years. The associations of LTL with imaging measurements and incident heart failure (HF) were evaluated by multivariable regression models. Genetic associations between LTL and significantly associated traits were investigated by mendelian randomization. Data were analyzed from January to May 2023.

Exposure

LTL.

Main Outcomes and Measures

Cardiovascular imaging traits and HF.

Results

Of 40 459 included participants, 19 529 (48.3%) were men, and the mean (SD) age was 55.1 (7.6) years. Longer LTL was independently associated with a pattern of positive cardiac remodeling (higher left ventricular mass, larger global ventricular size and volume, and higher ventricular and atrial stroke volumes) and a lower risk of incident HF (LTL fourth quartile vs first quartile: hazard ratio, 0.86; 95% CI, 0.81-0.91; P = 1.8 × 10-6). Mendelian randomization analysis suggested a potential causal association between LTL and left ventricular mass, global ventricular volume, and left ventricular stroke volume.

Conclusions and Relevance

In this cross-sectional study, longer LTL was associated with a larger heart with better cardiac function in middle age, which could potentially explain the observed lower risk of incident HF.

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.

Evaluation of interventions to prevent vasovagal reactions among whole blood donors: rationale and design of a large cluster randomised trial

BACKGROUND

Vasovagal reactions (VVRs) are the most common acute complications of blood donation. Responsible for substantial morbidity, they also reduce the likelihood of repeated donations and are disruptive and costly for blood services. Although blood establishments worldwide have adopted different strategies to prevent VVRs (including water loading and applied muscle tension [AMT]), robust evidence is limited. The Strategies to Improve Donor Experiences (STRIDES) trial aims to reliably assess the impact of four different interventions to prevent VVRs among blood donors.

METHODS

STRIDES is a cluster-randomised cross-over/stepped-wedge factorial trial of four interventions to reduce VVRs involving about 1.4 million whole blood donors enrolled from all 73 blood donation sites (mobile teams and donor centres) of National Health Service Blood and Transplant (NHSBT) in England. Each site ("cluster") has been randomly allocated to receive one or more interventions during a 36-month period, using principles of cross-over, stepped-wedge and factorial trial design to assign the sequence of interventions. Each of the four interventions is compared to NHSBT's current practices: (i) 500-ml isotonic drink before donation (vs current 500-ml plain water); (ii) 3-min rest on donation chair after donation (vs current 2 min); (iii) new modified AMT (vs current practice of AMT); and (iv) psychosocial intervention using preparatory materials (vs current practice of nothing). The primary outcome is the number of in-session VVRs with loss of consciousness (i.e. episodes involving loss of consciousness of any duration, with or without additional complications). Secondary outcomes include all in-session VVRs (i.e. with and without loss of consciousness), all delayed VVRs (i.e. those occurring after leaving the venue) and any in-session non-VVR adverse events or reactions.

DISCUSSION

The STRIDES trial should yield novel information about interventions, singly and in combination, for the prevention of VVRs, with the aim of generating policy-shaping evidence to help inform blood services to improve donor health, donor experience, and service efficiency.

TRIAL REGISTRATION

ISRCTN: 10412338. Registration date: October 24, 2019.

Using Polygenic Risk Scores for Prioritizing Individuals at Greatest Need of a Cardiovascular Disease Risk Assessment

Background The aim of this study was to provide quantitative evidence of the use of polygenic risk scores for systematically identifying individuals for invitation for full formal cardiovascular disease (CVD) risk assessment. Methods and Results A total of 108 685 participants aged 40 to 69 years, with measured biomarkers, linked primary care records, and genetic data in UK Biobank were used for model derivation and population health modeling. Prioritization tools using age, polygenic risk scores for coronary artery disease and stroke, and conventional risk factors for CVD available within longitudinal primary care records were derived using sex-specific Cox models. We modeled the implications of initiating guideline-recommended statin therapy after prioritizing individuals for invitation to a formal CVD risk assessment. If primary care records were used to prioritize individuals for formal risk assessment using age- and sex-specific thresholds corresponding to 5% false-negative rates, then the numbers of men and women needed to be screened to prevent 1 CVD event are 149 and 280, respectively. In contrast, adding polygenic risk scores to both prioritization and formal assessments, and selecting thresholds to capture the same number of events, resulted in a number needed to screen of 116 for men and 180 for women. Conclusions Using both polygenic risk scores and primary care records to prioritize individuals at highest risk of a CVD event for a formal CVD risk assessment can efficiently prioritize those who need interventions the most than using primary care records alone. This could lead to better allocation of resources by reducing the number of risk assessments in primary care while still preventing the same number of CVD events.

SCORE2-Diabetes: 10-year cardiovascular risk estimation in type 2 diabetes in Europe

AIMS

To develop and validate a recalibrated prediction model (SCORE2-Diabetes) to estimate the 10-year risk of cardiovascular disease (CVD) in individuals with type 2 diabetes in Europe.

METHODS AND RESULTS

SCORE2-Diabetes was developed by extending SCORE2 algorithms using individual-participant data from four large-scale datasets comprising 229 460 participants (43 706 CVD events) with type 2 diabetes and without previous CVD. Sex-specific competing risk-adjusted models were used including conventional risk factors (i.e. age, smoking, systolic blood pressure, total, and HDL-cholesterol), as well as diabetes-related variables (i.e. age at diabetes diagnosis, glycated haemoglobin [HbA1c] and creatinine-based estimated glomerular filtration rate [eGFR]). Models were recalibrated to CVD incidence in four European risk regions. External validation included 217 036 further individuals (38 602 CVD events), and showed good discrimination, and improvement over SCORE2 (C-index change from 0.009 to 0.031). Regional calibration was satisfactory. SCORE2-Diabetes risk predictions varied several-fold, depending on individuals' levels of diabetes-related factors. For example, in the moderate-risk region, the estimated 10-year CVD risk was 11% for a 60-year-old man, non-smoker, with type 2 diabetes, average conventional risk factors, HbA1c of 50 mmol/mol, eGFR of 90 mL/min/1.73 m2, and age at diabetes diagnosis of 60 years. By contrast, the estimated risk was 17% in a similar man, with HbA1c of 70 mmol/mol, eGFR of 60 mL/min/1.73 m2, and age at diabetes diagnosis of 50 years. For a woman with the same characteristics, the risk was 8% and 13%, respectively.

CONCLUSION

SCORE2-Diabetes, a new algorithm developed, calibrated, and validated to predict 10-year risk of CVD in individuals with type 2 diabetes, enhances identification of individuals at higher risk of developing CVD across Europe.

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.

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.

Impact of a post-donation hemoglobin testing strategy on efficiency and safety of whole blood donation in England: A modeling study

BACKGROUND

Deferrals due to low hemoglobin are time-consuming and costly for blood donors and donation services. Furthermore, accepting donations from those with low hemoglobin could represent a significant safety issue. One approach to reduce them is to use hemoglobin concentration alongside donor characteristics to inform personalized inter-donation intervals.

STUDY DESIGN AND METHODS

We used data from 17,308 donors to inform a discrete event simulation model comparing personalized inter-donation intervals using "post-donation" testing (i.e., estimating current hemoglobin from that measured by a hematology analyzer at last donation) versus the current approach in England (i.e., pre-donation testing with fixed intervals of 12-weeks for men and 16-weeks for women). We reported the impact on total donations, low hemoglobin deferrals, inappropriate bleeds, and blood service costs. Personalized inter-donation intervals were defined using mixed-effects modeling to estimate hemoglobin trajectories and probability of crossing hemoglobin donation thresholds.

RESULTS

The model had generally good internal validation, with predicted events similar to those observed. Over 1 year, a personalized strategy requiring ≥90% probability of being over the hemoglobin threshold, minimized adverse events (low hemoglobin deferrals and inappropriate bleeds) in both sexes and costs in women. Donations per adverse event improved from 3.4 (95% uncertainty interval 2.8, 3.7) under the current strategy to 14.8 (11.6, 19.2) in women, and from 7.1 (6.1, 8.5) to 26.9 (20.8, 42.6) in men. In comparison, a strategy incorporating early returns for those with high certainty of being over the threshold maximized total donations in both men and women, but was less favorable in terms of adverse events, with 8.4 donations per adverse event in women (7.0, 10,1) and 14.8 (12.1, 21.0) in men.

DISCUSSION

Personalized inter-donation intervals using post-donation testing combined with modeling of hemoglobin trajectories can help reduce deferrals, inappropriate bleeds, and costs.

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.

Including measures of chronic kidney disease to improve cardiovascular risk prediction by SCORE2 and SCORE2-OP

AIMS

The 2021 European Society of Cardiology (ESC) guideline on cardiovascular disease (CVD) prevention categorizes moderate and severe chronic kidney disease (CKD) as high and very-high CVD risk status regardless of other factors like age and does not include estimated glomerular filtration rate (eGFR) and albuminuria in its algorithms, systemic coronary risk estimation 2 (SCORE2) and systemic coronary risk estimation 2 in older persons (SCORE2-OP), to predict CVD risk. We developed and validated an 'Add-on' to incorporate CKD measures into these algorithms, using a validated approach.

METHODS

In 3,054 840 participants from 34 datasets, we developed three Add-ons [eGFR only, eGFR + urinary albumin-to-creatinine ratio (ACR) (the primary Add-on), and eGFR + dipstick proteinuria] for SCORE2 and SCORE2-OP. We validated C-statistics and net reclassification improvement (NRI), accounting for competing risk of non-CVD death, in 5,997 719 participants from 34 different datasets.

RESULTS

In the target population of SCORE2 and SCORE2-OP without diabetes, the CKD Add-on (eGFR only) and CKD Add-on (eGFR + ACR) improved C-statistic by 0.006 (95%CI 0.004-0.008) and 0.016 (0.010-0.023), respectively, for SCORE2 and 0.012 (0.009-0.015) and 0.024 (0.014-0.035), respectively, for SCORE2-OP. Similar results were seen when we included individuals with diabetes and tested the CKD Add-on (eGFR + dipstick). In 57 485 European participants with CKD, SCORE2 or SCORE2-OP with a CKD Add-on showed a significant NRI [e.g. 0.100 (0.062-0.138) for SCORE2] compared to the qualitative approach in the ESC guideline.

CONCLUSION

Our Add-ons with CKD measures improved CVD risk prediction beyond SCORE2 and SCORE2-OP. This approach will help clinicians and patients with CKD refine risk prediction and further personalize preventive therapies for CVD.

Estimating individual lifetime risk of incident cardiovascular events in adults with Type 2 diabetes: an update and geographical calibration of the DIAbetes Lifetime perspective model (DIAL2)

AIMS

The 2021 European Society of Cardiology cardiovascular disease (CVD) prevention guidelines recommend the use of (lifetime) risk prediction models to aid decisions regarding intensified preventive treatment options in adults with Type 2 diabetes, e.g. the DIAbetes Lifetime perspective model (DIAL model). The aim of this study was to update the DIAL model using contemporary and representative registry data (DIAL2) and to systematically calibrate the model for use in other European countries.

METHODS AND RESULTS

The DIAL2 model was derived in 467 856 people with Type 2 diabetes without a history of CVD from the Swedish National Diabetes Register, with a median follow-up of 7.3 years (interquartile range: 4.0-10.6 years) and comprising 63 824 CVD (including fatal CVD, non-fatal stroke and non-fatal myocardial infarction) events and 66 048 non-CVD mortality events. The model was systematically recalibrated to Europe's low- and moderate-risk regions using contemporary incidence data and mean risk factor distributions. The recalibrated DIAL2 model was externally validated in 218 267 individuals with Type 2 diabetes from the Scottish Care Information-Diabetes (SCID) and Clinical Practice Research Datalink (CPRD). In these individuals, 43 074 CVD events and 27 115 non-CVD fatal events were observed. The DIAL2 model discriminated well, with C-indices of 0.732 [95% confidence interval (CI) 0.726-0.739] in CPRD and 0.700 (95% CI 0.691-0.709) in SCID.

CONCLUSION

The recalibrated DIAL2 model provides a useful tool for the prediction of CVD-free life expectancy and lifetime CVD risk for people with Type 2 diabetes without previous CVD in the European low- and moderate-risk regions. These long-term individualized measures of CVD risk are well suited for shared decision-making in clinical practice as recommended by the 2021 CVD ESC prevention guidelines.

A Hybrid Architecture (CO-CONNECT) to Facilitate Rapid Discovery and Access to Data Across the United Kingdom in Response to the COVID-19 Pandemic: Development Study

BACKGROUND

COVID-19 data have been generated across the United Kingdom as a by-product of clinical care and public health provision, as well as numerous bespoke and repurposed research endeavors. Analysis of these data has underpinned the United Kingdom's response to the pandemic, and informed public health policies and clinical guidelines. However, these data are held by different organizations, and this fragmented landscape has presented challenges for public health agencies and researchers as they struggle to find relevant data to access and interrogate the data they need to inform the pandemic response at pace.

OBJECTIVE

We aimed to transform UK COVID-19 diagnostic data sets to be findable, accessible, interoperable, and reusable (FAIR).

METHODS

A federated infrastructure model (COVID - Curated and Open Analysis and Research Platform [CO-CONNECT]) was rapidly built to enable the automated and reproducible mapping of health data partners' pseudonymized data to the Observational Medical Outcomes Partnership Common Data Model without the need for any data to leave the data controllers' secure environments, and to support federated cohort discovery queries and meta-analysis.

RESULTS

A total of 56 data sets from 19 organizations are being connected to the federated network. The data include research cohorts and COVID-19 data collected through routine health care provision linked to longitudinal health care records and demographics. The infrastructure is live, supporting aggregate-level querying of data across the United Kingdom.

CONCLUSIONS

CO-CONNECT was developed by a multidisciplinary team. It enables rapid COVID-19 data discovery and instantaneous meta-analysis across data sources, and it is researching streamlined data extraction for use in a Trusted Research Environment for research and public health analysis. CO-CONNECT has the potential to make UK health data more interconnected and better able to answer national-level research questions while maintaining patient confidentiality and local governance procedures.

Obesity and cardiovascular disease: mechanistic insights and management strategies. A joint position paper by the World Heart Federation and World Obesity Federation

The ongoing obesity epidemic represents a global public health crisis that contributes to poor health outcomes, reduced quality of life, and >2.8 million deaths each year. Obesity is relapsing, progressive, and heterogeneous. It is considered a chronic disease by the World Obesity Federation (WOF) and a chronic condition by the World Heart Federation (WHF). People living with overweight/obesity are at greater risk for cardiovascular (CV) morbidity and mortality. Increased adiposity (body fat), particularly visceral/abdominal fat, is linked to CV risk and CV disease (CVD) via multiple direct and indirect pathophysiological mechanisms. The development of CVD is driven, in part, by obesity-related metabolic, endocrinologic, immunologic, structural, humoral, haemodynamic, and functional alterations. The complex multifaceted nature of these mechanisms can be challenging to understand and address in clinical practice. People living with obesity and CVD often have concurrent chronic physical or psychological disorders (multimorbidity) requiring multidisciplinary care pathways and polypharmacy. Evidence indicates that intentional weight loss (particularly when substantial) lowers CVD risk among people with overweight/obesity. Long-term weight loss and maintenance require ongoing commitment from both the individual and those responsible for their care. This position paper, developed by the WOF and the WHF, aims to improve understanding of the direct and indirect links between overweight/obesity and CVD, the key controversies in this area and evidence relating to cardiometabolic outcomes with available weight management options. Finally, an action plan for clinicians provides recommendations to help in identifying and addressing the risks of obesity-related CVD (recognizing resource and support variances between countries).

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.

New Way to "SCORE" Risk: Updates on the ESC Scoring System and Incorporation into ESC Cardiovascular Prevention Guidelines

PURPOSE OF REVIEW

Prior European Society of Cardiology (ESC) guidelines endorsed the SCORE 10-year cardiovascular disease (CVD) risk calculator to inform the total risk approach to CVD prevention, including the use of preventive interventions like lipid lowering therapies. However, SCORE was released in 2003, did not allow for estimation of fatal and non-fatal CVD events, and was limited to adults aged 40 to 70 years. The ESC's Cardiovascular Risk Collaboration (CRC) was tasked with updating SCORE (SCORE2) and with extending the upper age range of adults eligible for risk estimation (SCORE2-OP). This review summarises these two deliverables.

RECENT FINDINGS

Published in 2021, these updated risk scores allow for estimation of 10 year total (fatal + non-fatal) risks of a first atherosclerotic cardiovascular event in adults (SCORE2) and older persons (SCORE2-OP), calibrated for use in four European risk regions. The models account for competing risk of non-CVD death. These were extensively validated with excellent calibration and C-statistics ranging from 0.67 to 0.81. SCORE2 and SCORE2-OP have informed the 2021 ESC Guidelines on cardiovascular disease prevention in clinical practice. In addition to endorsing these two updated risk calculators, these guidelines have, for the first time, recommended the use of age-related risk categories. This change was motivated to prevent overreliance on age when making CVD prevention decisions.

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.

Association of COVID-19 With Major Arterial and Venous Thrombotic Diseases: A Population-Wide Cohort Study of 48 Million Adults in England and Wales

BACKGROUND

Infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) induces a prothrombotic state, but long-term effects of COVID-19 on incidence of vascular diseases are unclear.

METHODS

We studied vascular diseases after COVID-19 diagnosis in population-wide anonymized linked English and Welsh electronic health records from January 1 to December 7, 2020. We estimated adjusted hazard ratios comparing the incidence of arterial thromboses and venous thromboembolic events (VTEs) after diagnosis of COVID-19 with the incidence in people without a COVID-19 diagnosis. We conducted subgroup analyses by COVID-19 severity, demographic characteristics, and previous history.

RESULTS

Among 48 million adults, 125 985 were hospitalized and 1 319 789 were not hospitalized within 28 days of COVID-19 diagnosis. In England, there were 260 279 first arterial thromboses and 59 421 first VTEs during 41.6 million person-years of follow-up. Adjusted hazard ratios for first arterial thrombosis after COVID-19 diagnosis compared with no COVID-19 diagnosis declined from 21.7 (95% CI, 21.0-22.4) in week 1 after COVID-19 diagnosis to 1.34 (95% CI, 1.21-1.48) during weeks 27 to 49. Adjusted hazard ratios for first VTE after COVID-19 diagnosis declined from 33.2 (95% CI, 31.3-35.2) in week 1 to 1.80 (95% CI, 1.50-2.17) during weeks 27 to 49. Adjusted hazard ratios were higher, for longer after diagnosis, after hospitalized versus nonhospitalized COVID-19, among Black or Asian versus White people, and among people without versus with a previous event. The estimated whole-population increases in risk of arterial thromboses and VTEs 49 weeks after COVID-19 diagnosis were 0.5% and 0.25%, respectively, corresponding to 7200 and 3500 additional events, respectively, after 1.4 million COVID-19 diagnoses.

CONCLUSIONS

High relative incidence of vascular events soon after COVID-19 diagnosis declines more rapidly for arterial thromboses than VTEs. However, incidence remains elevated up to 49 weeks after COVID-19 diagnosis. These results support policies to prevent severe COVID-19 by means of COVID-19 vaccines, early review after discharge, risk factor control, and use of secondary preventive agents in high-risk patients.

Incremental value of risk factor variability for cardiovascular risk prediction in individuals with type 2 diabetes: results from UK primary care electronic health records

BACKGROUND

Cardiovascular disease (CVD) risk prediction models for individuals with type 2 diabetes are important tools to guide intensification of interventions for CVD prevention. We aimed to assess the added value of incorporating risk factors variability in CVD risk prediction for people with type 2 diabetes.

METHODS

We used electronic health records (EHRs) data from 83 910 adults with type 2 diabetes but without pre-existing CVD from the UK Clinical Practice Research Datalink for 2004-2017. Using a landmark-modelling approach, we developed and validated sex-specific Cox models, incorporating conventional predictors and trajectories plus variability of systolic blood pressure (SBP), total and high-density lipoprotein (HDL) cholesterol, and glycated haemoglobin (HbA1c). Such models were compared against simpler models using single last observed values or means.

RESULTS

The standard deviations (SDs) of SBP, HDL cholesterol and HbA1c were associated with higher CVD risk (P < 0.05). Models incorporating trajectories and variability of continuous predictors demonstrated improvement in risk discrimination (C-index = 0.659, 95% CI: 0.654-0.663) as compared with using last observed values (C-index = 0.651, 95% CI: 0.646-0.656) or means (C-index = 0.650, 95% CI: 0.645-0.655). Inclusion of SDs of SBP yielded the greatest improvement in discrimination (C-index increase = 0.005, 95% CI: 0.004-0.007) in comparison to incorporating SDs of total cholesterol (C-index increase = 0.002, 95% CI: 0.000-0.003), HbA1c (C-index increase = 0.002, 95% CI: 0.000-0.003) or HDL cholesterol (C-index increase= 0.003, 95% CI: 0.002-0.005).

CONCLUSION

Incorporating variability of predictors from EHRs provides a modest improvement in CVD risk discrimination for individuals with type 2 diabetes. Given that repeat measures are readily available in EHRs especially for regularly monitored patients with diabetes, this improvement could easily be achieved.

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.

Estimation of recurrent atherosclerotic cardiovascular event risk in patients with established cardiovascular disease: the updated SMART2 algorithm

AIMS

The 10-year risk of recurrent atherosclerotic cardiovascular disease (ASCVD) events in patients with established ASCVD can be estimated with the Secondary Manifestations of ARTerial disease (SMART) risk score, and may help refine clinical management. To broaden generalizability across regions, we updated the existing tool (SMART2 risk score) and recalibrated it with regional incidence rates and assessed its performance in external populations.

METHODS AND RESULTS

Individuals with coronary artery disease, cerebrovascular disease, peripheral artery disease, or abdominal aortic aneurysms were included from the Utrecht Cardiovascular Cohort-SMART cohort [n = 8355; 1706 ASCVD events during a median follow-up of 8.2 years (interquartile range 4.2-12.5)] to derive a 10-year risk prediction model for recurrent ASCVD events (non-fatal myocardial infarction, non-fatal stroke, or cardiovascular mortality) using a Fine and Gray competing risk-adjusted model. The model was recalibrated to four regions across Europe, and to Asia (excluding Japan), Japan, Australia, North America, and Latin America using contemporary cohort data from each target region. External validation used data from seven cohorts [Clinical Practice Research Datalink, SWEDEHEART, the international REduction of Atherothrombosis for Continued Health (REACH) Registry, Estonian Biobank, Spanish Biomarkers in Acute Coronary Syndrome and Biomarkers in Acute Myocardial Infarction (BACS/BAMI), the Norwegian COgnitive Impairment After STroke, and Bialystok PLUS/Polaspire] and included 369 044 individuals with established ASCVD of whom 62 807 experienced an ASCVD event. C-statistics ranged from 0.605 [95% confidence interval (CI) 0.547-0.664] in BACS/BAMI to 0.772 (95% CI 0.659-0.886) in REACH Europe high-risk region. The clinical utility of the model was demonstrated across a range of clinically relevant treatment thresholds for intensified treatment options.

CONCLUSION

The SMART2 risk score provides an updated, validated tool for the prediction of recurrent ASCVD events in patients with established ASCVD across European and non-European populations. The use of this tool could allow for a more personalized approach to secondary prevention based upon quantitative rather than qualitative estimates of residual risk.

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.

Investigating Genetic and Other Determinants of First-Onset Myocardial Infarction in Malaysia: Protocol for the Malaysian Acute Vascular Events Risk Study

BACKGROUND

Although the burden of premature myocardial infarction (MI) is high in Malaysia, direct evidence on the determinants of MI in this multi-ethnic population remains sparse.

OBJECTIVE

The Malaysian Acute Vascular Events Risk (MAVERIK) study is a retrospective case-control study established to investigate the genomic, lipid-related, and other determinants of acute MI in Malaysia. In this paper, we report the study protocol and early results.

METHODS

By June 2019, we had enrolled approximately 2500 patients with their first MI and 2500 controls without cardiovascular disease, who were frequency-matched by age, sex, and ethnicity, from 17 hospitals in Malaysia. For each participant, serum and whole blood have been collected and stored. Clinical, demographic, and behavioral information has been obtained using a 200-item questionnaire.

RESULTS

Tobacco consumption, a history of diabetes, hypertension, markers of visceral adiposity, indicators of lower socioeconomic status, and a family history of coronary disease were more prevalent in cases than in controls. Adjusted (age and sex) logistic regression models for traditional risk factors indicated that current smoking (odds ratio [OR] 4.11, 95% CI 3.56-4.75; P<.001), previous smoking (OR 1.34, 95% CI 1.12-1.60; P=.001), a history of high blood pressure (OR 2.13, 95% CI 1.86-2.44; P<.001), a history of diabetes mellitus (OR 2.72, 95% CI 2.34-3.17; P<.001), a family history of coronary heart disease (OR 1.28, 95% CI 1.07-1.55; P=.009), and obesity (BMI >30 kg/m²; OR 1.19, 95% CI 1.05-1.34; P=.009) were associated with MI in age- and sex-adjusted models.

CONCLUSIONS

The MAVERIK study can serve as a useful platform to investigate genetic and other risk factors for MI in an understudied Southeast Asian population. It should help to hasten the discovery of disease-causing pathways and inform regionally appropriate strategies that optimize public health action.

INTERNATIONAL REGISTERED REPORT IDENTIFIER (IRRID)

RR1-10.2196/31885.

An Expanded Genome-Wide Association Study of Fructosamine Levels Identifies RCN3 as a Replicating Locus and Implicates FCGRT as the Effector Transcript

Fructosamine is a measure of short-term glycemic control, which has been suggested as a useful complement to glycated hemoglobin (HbA1c) for the diagnosis and monitoring of diabetes. To date, a single genome-wide association study (GWAS) including 8,951 U.S. White and 2,712 U.S. Black individuals without a diabetes diagnosis has been published. Results in Whites and Blacks yielded different association loci, near RCN3 and CNTN5, respectively. In this study, we performed a GWAS on 20,731 European-ancestry blood donors and meta-analyzed our results with previous data from U.S. White participants from the Atherosclerosis Risk in Communities (ARIC) study (Nmeta = 29,685). We identified a novel association near GCK (rs3757840, βmeta = 0.0062; minor allele frequency [MAF] = 0.49; Pmeta = 3.66 × 10-8) and confirmed the association near RCN3 (rs113886122, βmeta = 0.0134; MAF = 0.17; Pmeta = 5.71 × 10-18). Colocalization analysis with whole-blood expression quantitative trait loci data suggested FCGRT as the effector transcript at the RCN3 locus. We further showed that fructosamine has low heritability (h2 = 7.7%), has no significant genetic correlation with HbA1c and other glycemic traits in individuals without a diabetes diagnosis (P > 0.05), but has evidence of shared genetic etiology with some anthropometric traits (Bonferroni-corrected P < 0.0012). Our results broaden knowledge of the genetic architecture of fructosamine and prioritize FCGRT for downstream functional studies at the established RCN3 locus.

Association of COVID-19 vaccines ChAdOx1 and BNT162b2 with major venous, arterial, or thrombocytopenic events: A population-based cohort study of 46 million adults in England

BACKGROUND

Thromboses in unusual locations after the Coronavirus Disease 2019 (COVID-19) vaccine ChAdOx1-S have been reported, although their frequency with vaccines of different types is uncertain at a population level. The aim of this study was to estimate the population-level risks of hospitalised thrombocytopenia and major arterial and venous thromboses after COVID-19 vaccination.

METHODS AND FINDINGS

In this whole-population cohort study, we analysed linked electronic health records from adults living in England, from 8 December 2020 to 18 March 2021. We estimated incidence rates and hazard ratios (HRs) for major arterial, venous, and thrombocytopenic outcomes 1 to 28 and >28 days after first vaccination dose for ChAdOx1-S and BNT162b2 vaccines. Analyses were performed separately for ages <70 and ≥70 years and adjusted for age, age2, sex, ethnicity, and deprivation. We also prespecified adjustment for anticoagulant medication, combined oral contraceptive medication, hormone replacement therapy medication, history of pulmonary embolism or deep vein thrombosis, and history of coronavirus infection in analyses of venous thrombosis; and diabetes, hypertension, smoking, antiplatelet medication, blood pressure lowering medication, lipid lowering medication, anticoagulant medication, history of stroke, and history of myocardial infarction in analyses of arterial thromboses. We selected further covariates with backward selection. Of 46 million adults, 23 million (51%) were women; 39 million (84%) were <70; and 3.7 million (8.1%) Asian or Asian British, 1.6 million (3.5%) Black or Black British, 36 million (79%) White, 0.7 million (1.5%) mixed ethnicity, and 1.5 million (3.2%) were of another ethnicity. Approximately 21 million (46%) adults had their first vaccination between 8 December 2020 and 18 March 2021. The crude incidence rates (per 100,000 person-years) of all venous events were as follows: prevaccination, 140 [95% confidence interval (CI): 138 to 142]; ≤28 days post-ChAdOx1-S, 294 (281 to 307); >28 days post-ChAdOx1-S, 359 (338 to 382), ≤28 days post-BNT162b2-S, 241 (229 to 253); >28 days post-BNT162b2-S 277 (263 to 291). The crude incidence rates (per 100,000 person-years) of all arterial events were as follows: prevaccination, 546 (95% CI: 541 to 555); ≤28 days post-ChAdOx1-S, 1,211 (1,185 to 1,237); >28 days post-ChAdOx1-S, 1678 (1,630 to 1,726), ≤28 days post-BNT162b2-S, 1,242 (1,214 to 1,269); >28 days post-BNT162b2-S, 1,539 (1,507 to 1,572). Adjusted HRs (aHRs) 1 to 28 days after ChAdOx1-S, compared with unvaccinated rates, at ages <70 and ≥70 years, respectively, were 0.97 (95% CI: 0.90 to 1.05) and 0.58 (0.53 to 0.63) for venous thromboses, and 0.90 (0.86 to 0.95) and 0.76 (0.73 to 0.79) for arterial thromboses. Corresponding aHRs for BNT162b2 were 0.81 (0.74 to 0.88) and 0.57 (0.53 to 0.62) for venous thromboses, and 0.94 (0.90 to 0.99) and 0.72 (0.70 to 0.75) for arterial thromboses. aHRs for thrombotic events were higher at younger ages for venous thromboses after ChAdOx1-S, and for arterial thromboses after both vaccines. Rates of intracranial venous thrombosis (ICVT) and of thrombocytopenia in adults aged <70 years were higher 1 to 28 days after ChAdOx1-S (aHRs 2.27, 95% CI: 1.33 to 3.88 and 1.71, 1.35 to 2.16, respectively), but not after BNT162b2 (0.59, 0.24 to 1.45 and 1.00, 0.75 to 1.34) compared with unvaccinated. The corresponding absolute excess risks of ICVT 1 to 28 days after ChAdOx1-S were 0.9 to 3 per million, varying by age and sex. The main limitations of the study are as follows: (i) it relies on the accuracy of coded healthcare data to identify exposures, covariates, and outcomes; (ii) the use of primary reason for hospital admission to measure outcome, which improves the positive predictive value but may lead to an underestimation of incidence; and (iii) potential unmeasured confounding.

CONCLUSIONS

In this study, we observed increases in rates of ICVT and thrombocytopenia after ChAdOx1-S vaccination in adults aged <70 years that were small compared with its effect in reducing COVID-19 morbidity and mortality, although more precise estimates for adults aged <40 years are needed. For people aged ≥70 years, rates of arterial or venous thrombotic events were generally lower after either vaccine compared with unvaccinated, suggesting that either vaccine is suitable in this age group.

Machine learning optimized polygenic scores for blood cell traits identify sex-specific trajectories and genetic correlations with disease

Genetic association studies for blood cell traits, which are key indicators of health and immune function, have identified several hundred associations and defined a complex polygenic architecture. Polygenic scores (PGSs) for blood cell traits have potential clinical utility in disease risk prediction and prevention, but designing PGS remains challenging and the optimal methods are unclear. To address this, we evaluated the relative performance of 6 methods to develop PGS for 26 blood cell traits, including a standard method of pruning and thresholding (P + T) and 5 learning methods: LDpred2, elastic net (EN), Bayesian ridge (BR), multilayer perceptron (MLP) and convolutional neural network (CNN). We evaluated these optimized PGSs on blood cell trait data from UK Biobank and INTERVAL. We find that PGSs designed using common machine learning methods EN and BR show improved prediction of blood cell traits and consistently outperform other methods. Our analyses suggest EN/BR as the top choices for PGS construction, showing improved performance for 25 blood cell traits in the external validation, with correlations with the directly measured traits increasing by 10%-23%. Ten PGSs showed significant statistical interaction with sex, and sex-specific PGS stratification showed that all of them had substantial variation in the trajectories of blood cell traits with age. Genetic correlations between the PGSs for blood cell traits and common human diseases identified well-known as well as new associations. We develop machine learning-optimized PGS for blood cell traits, demonstrate their relationships with sex, age, and disease, and make these publicly available as a resource.

Cost-effectiveness of alternative minimum recall intervals between whole blood donations

BACKGROUND

The INTERVAL trial showed shorter inter-donation intervals could safely increase the frequency of whole-blood donation. We extended the INTERVAL trial to consider the relative cost-effectiveness of reduced inter-donation intervals.

METHODS

Our within-trial cost-effectiveness analysis (CEA) used data from 44,863 whole-blood donors randomly assigned to 12, 10 or 8 week (males), and 16, 14 or 12 week inter-donation intervals (females). The CEA analysed the number of whole-blood donations, deferrals including low- haemoglobin deferrals, and donors' health-related quality of life (QoL) to report costs and cost-effectiveness over two years.

FINDINGS

The mean number of blood donation visits over two years was higher for the reduced interval strategies, for males (7.76, 6.60 and 5.68 average donations in the 8-, 10- and 12- week arms) and for females (5.10, 4.60 and 4.01 donations in the 12-, 14- and 16- week arms). For males, the average rate of deferral for low haemoglobin per session attended, was 5.71% (8- week arm), 3.73% (10- week), and 2.55% (12- week), and for females the rates were: 7.92% (12-week), 6.63% (14- week), and 5.05% (16- week). Donors' QoL was similar across strategies, although self-reported symptoms were increased with shorter donation intervals. The shorter interval strategies increased average cost, with incremental cost-effectiveness ratios of £9.51 (95% CI 9.33 to 9.69) per additional whole-blood donation for the 8- versus 12- week interval for males, and £10.17 (95% CI 9.80 to 10.54) for the 12- versus 16- week interval arm for females.

CONCLUSIONS

Over two years, reducing the minimum donation interval could provide additional units of whole-blood at a small additional cost, including for those donor subgroups whose blood type is in relatively high demand. However, the significance of self-reported symptoms needs to be investigated further before these policies are expanded.

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, stroke, or all-cause mortality. However, for the participants with vitamin D deficiency (25[OH]D concentration <25 nmol/L), genetic analyses provided strong evidence for an inverse association with all-cause mortality (odds ratio [OR] per 10 nmol/L increase in genetically-predicted 25[OH]D concentration 0·69 [95% CI 0·59-0·80]; p<0·0001) and non-significant inverse associations for stroke (0·85 [0·70-1·02], p=0·09) and coronary heart disease (0·89 [0·76-1·04]; p=0·14). A finer stratification of participants found inverse associations between genetically-predicted 25(OH)D concentrations and all-cause mortality up to around 40 nmol/L.

INTERPRETATION

Stratified Mendelian randomisation analyses suggest a causal relationship between 25(OH)D concentrations and mortality for individuals with low vitamin D status. Our findings have implications for the design of vitamin D supplementation trials, and potential disease prevention strategies.

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.

Integrative analysis of the plasma proteome and polygenic risk of cardiometabolic diseases

Cardiometabolic diseases are frequently polygenic in architecture, comprising a large number of risk alleles with small effects spread across the genome1-3. Polygenic scores (PGS) aggregate these into a metric representing an individual's genetic predisposition to disease. PGS have shown promise for early risk prediction4-7 and there is an open question as to whether PGS can also be used to understand disease biology8. Here, we demonstrate that cardiometabolic disease PGS can be used to elucidate the proteins underlying disease pathogenesis. In 3,087 healthy individuals, we found that PGS for coronary artery disease, type 2 diabetes, chronic kidney disease and ischaemic stroke are associated with the levels of 49 plasma proteins. Associations were polygenic in architecture, largely independent of cis and trans protein quantitative trait loci and present for proteins without quantitative trait loci. Over a follow-up of 7.7 years, 28 of these proteins associated with future myocardial infarction or type 2 diabetes events, 16 of which were mediators between polygenic risk and incident disease. Twelve of these were druggable targets with therapeutic potential. Our results demonstrate the potential for PGS to uncover causal disease biology and targets with therapeutic potential, including those that may be missed by approaches utilizing information at a single locus.

Genetically Predicted Type 2 Diabetes Mellitus Liability, Glycated Hemoglobin and Cardiovascular Diseases: A Wide-Angled Mendelian Randomization Study

(1) Aim: To investigate the causal effects of T2DM liability and glycated haemoglobin (HbA1c) levels on various cardiovascular disease outcomes, both in the general population and in non-diabetic individuals specifically. (2) Methods: We selected 243 variants as genetic instruments for T2DM liability and 536 variants for HbA1c. Linear Mendelian randomization analyses were performed to estimate the associations of genetically-predicted T2DM liability and HbA1c with 12 cardiovascular disease outcomes in 367,703 unrelated UK Biobank participants of European ancestries. We performed secondary analyses in participants without diabetes (HbA1c < 6.5% with no diagnosed diabetes), and in participants without diabetes or pre-diabetes (HbA1c < 5.7% with no diagnosed diabetes). (3) Results: Genetically-predicted T2DM liability was positively associated (p < 0.004, 0.05/12) with peripheral vascular disease, aortic valve stenosis, coronary artery disease, heart failure, ischaemic stroke, and any stroke. Genetically-predicted HbA1c was positively associated with coronary artery disease and any stroke. Mendelian randomization estimates generally shifted towards the null when excluding diabetic and pre-diabetic participants from analyses. (4) Conclusions: This genetic evidence supports causal effects of T2DM liability and HbA1c on a range of cardiovascular diseases, suggesting that improving glycaemic control could reduce cardiovascular risk in a general population, with greatest benefit in individuals with diabetes.

Prediction of Cardiovascular Disease Risk Accounting for Future Initiation of Statin Treatment

Cardiovascular disease (CVD) risk-prediction models are used to identify high-risk individuals and guide statin initiation. However, these models are usually derived from individuals who might initiate statins during follow-up. We present a simple approach to address statin initiation to predict “statin-naive” CVD risk. We analyzed primary care data (2004–2017) from the UK Clinical Practice Research Datalink for 1,678,727 individuals (aged 40–85 years) without CVD or statin treatment history at study entry. We derived age- and sex-specific prediction models including conventional risk factors and a time-dependent effect of statin initiation constrained to 25% risk reduction (from trial results). We compared predictive performance and measures of public-health impact (e.g., number needed to screen to prevent 1 event) against models ignoring statin initiation. During a median follow-up of 8.9 years, 103,163 individuals developed CVD. In models accounting for (versus ignoring) statin initiation, 10-year CVD risk predictions were slightly higher; predictive performance was moderately improved. However, few individuals were reclassified to a high-risk threshold, resulting in negligible improvements in number needed to screen to prevent 1 event. In conclusion, incorporating statin effects from trial results into risk-prediction models enables statin-naive CVD risk estimation and provides moderate gains in predictive ability but had a limited impact on treatment decision-making under current guidelines in this population.

Effects of adiposity on the human plasma proteome: observational and Mendelian randomisation estimates

BACKGROUND

Variation in adiposity is associated with cardiometabolic disease outcomes, but mechanisms leading from this exposure to disease are unclear. This study aimed to estimate effects of body mass index (BMI) on an extensive set of circulating proteins.

METHODS

We used SomaLogic proteomic data from up to 2737 healthy participants from the INTERVAL study. Associations between self-reported BMI and 3622 unique plasma proteins were explored using linear regression. These were complemented by Mendelian randomisation (MR) analyses using a genetic risk score (GRS) comprised of 654 BMI-associated polymorphisms from a recent genome-wide association study (GWAS) of adult BMI. A disease enrichment analysis was performed using DAVID Bioinformatics 6.8 for proteins which were altered by BMI.

RESULTS

Observationally, BMI was associated with 1576 proteins (P < 1.4 × 10-5), with particularly strong evidence for a positive association with leptin and fatty acid-binding protein-4 (FABP4), and a negative association with sex hormone-binding globulin (SHBG). Observational estimates were likely confounded, but the GRS for BMI did not associate with measured confounders. MR analyses provided evidence for a causal relationship between BMI and eight proteins including leptin (0.63 standard deviation (SD) per SD BMI, 95% CI 0.48-0.79, P = 1.6 × 10-15), FABP4 (0.64 SD per SD BMI, 95% CI 0.46-0.83, P = 6.7 × 10-12) and SHBG (-0.45 SD per SD BMI, 95% CI -0.65 to -0.25, P = 1.4 × 10-5). There was agreement in the magnitude of observational and MR estimates (R2 = 0.33) and evidence that proteins most strongly altered by BMI were enriched for genes involved in cardiovascular disease.

CONCLUSIONS

This study provides evidence for a broad impact of adiposity on the human proteome. Proteins strongly altered by BMI include those involved in regulating appetite, sex hormones and inflammation; such proteins are also enriched for cardiovascular disease-related genes. Altogether, results help focus attention onto new proteomic signatures of obesity-related disease.