Synergistic insights into human health from aptamer- and antibody-based proteomic profiling

A review on trends in development and translation of omics signatures in cancer

The field of cancer genomics and transcriptomics has evolved from targeted profiling to swift sequencing of individual tumor genome and transcriptome. The steady growth in genome, epigenome, and transcriptome datasets on a genome-wide scale has significantly increased our capability in capturing signatures that represent both the intrinsic and extrinsic biological features of tumors. These biological differences can help in precise molecular subtyping of cancer, predicting tumor progression, metastatic potential, and resistance to therapeutic agents. In this review, we summarized the current development of genomic, methylomic, transcriptomic, proteomic and metabolic signatures in the field of cancer research and highlighted their potentials in clinical applications to improve diagnosis, prognosis, and treatment decision in cancer patients.

Proteomic aging clock predicts mortality and risk of common age-related diseases in diverse populations

Circulating plasma proteins play key roles in human health and can potentially be used to measure biological age, allowing risk prediction for age-related diseases, multimorbidity and mortality. Here we developed a proteomic age clock in the UK Biobank (n = 45,441) using a proteomic platform comprising 2,897 plasma proteins and explored its utility to predict major disease morbidity and mortality in diverse populations. We identified 204 proteins that accurately predict chronological age (Pearson r = 0.94) and found that proteomic aging was associated with the incidence of 18 major chronic diseases (including diseases of the heart, liver, kidney and lung, diabetes, neurodegeneration and cancer), as well as with multimorbidity and all-cause mortality risk. Proteomic aging was also associated with age-related measures of biological, physical and cognitive function, including telomere length, frailty index and reaction time. Proteins contributing most substantially to the proteomic age clock are involved in numerous biological functions, including extracellular matrix interactions, immune response and inflammation, hormone regulation and reproduction, neuronal structure and function and development and differentiation. In a validation study involving biobanks in China (n = 3,977) and Finland (n = 1,990), the proteomic age clock showed similar age prediction accuracy (Pearson r = 0.92 and r = 0.94, respectively) compared to its performance in the UK Biobank. Our results demonstrate that proteomic aging involves proteins spanning multiple functional categories and can be used to predict age-related functional status, multimorbidity and mortality risk across geographically and genetically diverse populations.

Whole genome sequencing based analysis of inflammation biomarkers in the Trans-Omics for Precision Medicine (TOPMed) consortium

Inflammation biomarkers can provide valuable insight into the role of inflammatory processes in many diseases and conditions. Sequencing based analyses of such biomarkers can also serve as an exemplar of the genetic architecture of quantitative traits. To evaluate the biological insight, which can be provided by a multi-ancestry, whole-genome based association study, we performed a comprehensive analysis of 21 inflammation biomarkers from up to 38 465 individuals with whole-genome sequencing from the Trans-Omics for Precision Medicine (TOPMed) program (with varying sample size by trait, where the minimum sample size was n = 737 for MMP-1). We identified 22 distinct single-variant associations across 6 traits-E-selectin, intercellular adhesion molecule 1, interleukin-6, lipoprotein-associated phospholipase A2 activity and mass, and P-selectin-that remained significant after conditioning on previously identified associations for these inflammatory biomarkers. We further expanded upon known biomarker associations by pairing the single-variant analysis with a rare variant set-based analysis that further identified 19 significant rare variant set-based associations with 5 traits. These signals were distinct from both significant single variant association signals within TOPMed and genetic signals observed in prior studies, demonstrating the complementary value of performing both single and rare variant analyses when analyzing quantitative traits. We also confirm several previously reported signals from semi-quantitative proteomics platforms. Many of these signals demonstrate the extensive allelic heterogeneity and ancestry-differentiated variant-trait associations common for inflammation biomarkers, a characteristic we hypothesize will be increasingly observed with well-powered, large-scale analyses of complex traits.

Comparative Analysis of Protein Quantification by the SomaScan Assay versus Orthogonal Methods in Urine from People with Diabetic Kidney Disease

To our knowledge, calibration curves or other validations for thousands of SomaScan aptamers are not publicly available. Moreover, the abundance of urine proteins obtained from these assays is not routinely validated with orthogonal methods (OMs). We report an in-depth comparison of SomaScan readout for 23 proteins in urine samples from patients with diabetic kidney disease (n = 118) vs OMs, including liquid chromatography-targeted mass spectrometry (LC-MS), ELISA, and nephelometry. Pearson correlation between urine abundance of the 23 proteins from SomaScan 3.2 vs OMs ranged from -0.58 to 0.86, with a median (interquartile ratio, [IQR]) of 0.49 (0.18, 0.53). In multivariable linear regression, the SomaScan readout for 6 of the 23 examined proteins (26%) was most strongly associated with the OM-derived abundance of the same (target) protein. For 3 of 23 (13%), the SomaScan and OM-derived abundance of each protein were significantly associated, but the SomaScan readout was more strongly associated with OM-derived abundance of one or more "off-target" proteins. For the remaining 14 proteins (61%), the SomaScan readouts were not significantly associated with the OM-derived abundance of the targeted proteins. In 6 of the latest group, the SomaScan readout was not associated with urine abundance of any of the 23 quantified proteins. To sum, over half of the SomaScan results could not be confirmed by independent orthogonal methods.

Blood protein assessment of leading incident diseases and mortality in the UK Biobank

The circulating proteome offers insights into the biological pathways that underlie disease. Here, we test relationships between 1,468 Olink protein levels and the incidence of 23 age-related diseases and mortality in the UK Biobank (n = 47,600). We report 3,209 associations between 963 protein levels and 21 incident outcomes. Next, protein-based scores (ProteinScores) are developed using penalized Cox regression. When applied to test sets, six ProteinScores improve the area under the curve estimates for the 10-year onset of incident outcomes beyond age, sex and a comprehensive set of 24 lifestyle factors, clinically relevant biomarkers and physical measures. Furthermore, the ProteinScore for type 2 diabetes outperforms a polygenic risk score and HbA1c-a clinical marker used to monitor and diagnose type 2 diabetes. The performance of scores using metabolomic and proteomic features is also compared. These data characterize early proteomic contributions to major age-related diseases, demonstrating the value of the plasma proteome for risk stratification.

Targeted proteomics involved in cardiovascular health and heart rate variability in children with overweight/obesity

BACKGROUND

Children with overweight/obesity often exhibit alterations in their plasma protein profiles and reduced heart rate variability (HRV). Plasma proteomics is at the forefront of identifying biomarkers for various clinical conditions. We aimed to examine the association between plasma-targeted proteomics involved in cardiovascular health and resting vagal-related HRV parameters in children with overweight/obesity.

METHODS

Forty-four children with overweight/obesity (10.2 ± 1.1 years old; 52% boys) participated in the study. Olink's technology was used to quantify 92 proteins involved in cardiovascular health. HRV was measured using a heart rate monitor (Polar RS800CX). Four resting vagal-related HRV parameters were derived in time- and frequency-domain.

RESULTS

Eight proteins (KIM1, IgG Fc receptor II-b, IDUA, BOC, IL1RL2, TNFRSF11A, VSIG2, and TF) were associated with at least one out of the four vagal-related HRV parameters (β values ranging from -0.188 to 0.288; all p < .05), while KIM1, IDUA, and BOC associated with ≥ three vagal-related HRV parameters. Multiple hypothesis testing corrections did not reach statistical significance (false discovery rate [FDR >0.05]).

CONCLUSION

Plasma-targeted proteomics suggested novel biomarkers for resting vagal-related HRV parameters in children with overweight/obesity. Future studies using larger cohorts and longitudinal designs should confirm our findings and their potential clinical implications.

Dysregulated STAT3 signaling and T cell immunometabolic dysfunction define a targetable, high mortality subphenotype of critically ill children

Sepsis is the leading cause of death of hospitalized children worldwide. Despite the established link between immune dysregulation and mortality in pediatric sepsis, it remains unclear which host immune factors contribute causally to adverse sepsis outcomes. Identifying modifiable pathobiology is an essential first step to successful translation of biologic insights into precision therapeutics. We designed a prospective, longitudinal cohort study of 88 critically ill pediatric patients with multiple organ dysfunction syndrome (MODS), including patients with and without sepsis, to define subphenotypes associated with targetable mechanisms of immune dysregulation. We first assessed plasma proteomic profiles and identified shared features of immune dysregulation in MODS patients with and without sepsis. We then employed consensus clustering to define three subphenotypes based on protein expression at disease onset and identified a strong association between subphenotype and clinical outcome. We next identified differences in immune cell frequency and activation state by MODS subphenotype and determined the association between hyperinflammatory pathway activation and cellular immunophenotype. Using single cell transcriptomics, we demonstrated STAT3 hyperactivation in lymphocytes from the sickest MODS subgroup and then identified an association between STAT3 hyperactivation and T cell immunometabolic dysregulation. Finally, we compared proteomics findings between patients with MODS and patients with inborn errors of immunity that amplify cytokine signaling pathways to further assess the impact of STAT3 hyperactivation in the most severe patients with MODS. Overall, these results identify a potentially pathologic and targetable role for STAT3 hyperactivation in a subset of pediatric patients with MODS who have high severity of illness and poor prognosis.

Proteomic analysis of cardiorespiratory fitness for prediction of mortality and multisystem disease risks

Despite the wide effects of cardiorespiratory fitness (CRF) on metabolic, cardiovascular, pulmonary and neurological health, challenges in the feasibility and reproducibility of CRF measurements have impeded its use for clinical decision-making. Here we link proteomic profiles to CRF in 14,145 individuals across four international cohorts with diverse CRF ascertainment methods to establish, validate and characterize a proteomic CRF score. In a cohort of around 22,000 individuals in the UK Biobank, a proteomic CRF score was associated with a reduced risk of all-cause mortality (unadjusted hazard ratio 0.50 (95% confidence interval 0.48-0.52) per 1 s.d. increase). The proteomic CRF score was also associated with multisystem disease risk and provided risk reclassification and discrimination beyond clinical risk factors, as well as modulating high polygenic risk of certain diseases. Finally, we observed dynamicity of the proteomic CRF score in individuals who undertook a 20-week exercise training program and an association of the score with the degree of the effect of training on CRF, suggesting potential use of the score for personalization of exercise recommendations. These results indicate that population-based proteomics provides biologically relevant molecular readouts of CRF that are additive to genetic risk, potentially modifiable and clinically translatable.

Identifying novel biomarkers using proteomics to predict cancer-associated thrombosis

Comprehensive protein analyses of plasma are made possible by high-throughput proteomic screens, which may help find new therapeutic targets and diagnostic biomarkers. Patients with cancer are frequently affected by venous thromboembolism (VTE). The limited predictive accuracy of current VTE risk assessment tools highlights the need for new, more targeted biomarkers. Although coagulation biomarkers for the diagnosis, prognosis, and treatment of VTE have been investigated, none of them have the necessary clinical validation or diagnostic accuracy. Proteomics holds the potential to uncover new biomarkers and thrombotic pathways that impact the risk of thrombosis. This review explores the fundamental methods used in proteomics and focuses on particular biomarkers found in VTE and cancer-associated thrombosis.

Ratio of plasma IL-13/TNF- ∝ and CXCL10/CCL17 predicts mepolizumab and omalizumab response in asthma better than eosinophil count or immunoglobulin E level

To date, most studies to identify biomarkers associated with response to the anti-interleukin 5 agent, mepolizumab, and to the anti-immunoglobulin E agent, omalizumab have focused on clinically available biomarkers, such as the peripheral blood eosinophil counts (BEC) and total immunoglobulin E (IgE). However, these biomarkers often have low predictive accuracy, with many patients with eosinophilic or allergic asthma failing to demonstrate clinical response to mepolizumab or omalizumab respectively. In this study, we evaluated the association of baseline pre-biologic plasma levels of 26 cytokines and chemokines, including T-helper 1 (Th1)-, Th2-, Th17-related cytokines, and their ratios with subsequent clinical response to mepolizumab or omalizumab. We defined clinical response as a reduction in the baseline annual exacerbation rate by half or more over the one-year period following initiation of the biologic. Baseline levels of plasma IL-13 were differentially elevated in responders versus non-responders to mepolizumab and plasma CXCL10 levels were differentially elevated in responders to omalizumab. The ratio of IL-13/TNF-α had the best sensitivity and specificity in predicting response to mepolizumab and CXCL10/CCL17 to omalizumab, and these performed better as predictive biomarkers of response than BEC and IgE. Cytokines and chemokines associated with airway eosinophilia, allergic inflammation, or Th2 inflammation, such as IL-13 and CXCL10, may be better predictors of clinical response to mepolizumab and omalizumab, than IL-5 or IgE, the targets of mepolizumab and omalizumab.

Harnessing the power of proteomics in precision diabetes medicine

Precision diabetes medicine (PDM) aims to reduce errors in prevention programmes, diagnosis thresholds, prognosis prediction and treatment strategies. However, its advancement and implementation are difficult due to the heterogeneity of complex molecular processes and environmental exposures that influence an individual's disease trajectory. To address this challenge, it is imperative to develop robust screening methods for all areas of PDM. Innovative proteomic technologies, alongside genomics, have proven effective in precision cancer medicine and are showing promise in diabetes research for potential translation. This narrative review highlights how proteomics is well-positioned to help improve PDM. Specifically, a critical assessment of widely adopted affinity-based proteomic technologies in large-scale clinical studies and evidence of the benefits and feasibility of using MS-based plasma proteomics is presented. We also present a case for the use of proteomics to identify predictive protein panels for type 2 diabetes subtyping and the development of clinical prediction models for prevention, diagnosis, prognosis and treatment strategies. Lastly, we discuss the importance of plasma and tissue proteomics and its integration with genomics (proteogenomics) for identifying unique type 2 diabetes intra- and inter-subtype aetiology. We conclude with a call for action formed on advancing proteomics technologies, benchmarking their performance and standardisation across sites, with an emphasis on data sharing and the inclusion of diverse ancestries in large cohort studies. These efforts should foster collaboration with key stakeholders and align with ongoing academic programmes such as the Precision Medicine in Diabetes Initiative consortium.

Proteomics applications in next generation induced pluripotent stem cell models

INTRODUCTION

Induced pluripotent stem (iPS) cell technology has transformed biomedical research. New opportunities now exist to create new organoids, microtissues, and body-on-a-chip systems for basic biology investigations and clinical translations.

AREAS COVERED

We discuss the utility of proteomics for attaining an unbiased view into protein expression changes during iPS cell differentiation, cell maturation, and tissue generation. The ability to discover cell-type specific protein markers during the differentiation and maturation of iPS-derived cells has led to new strategies to improve cell production yield and fidelity. In parallel, proteomic characterization of iPS-derived organoids is helping to realize the goal of bridging in vitro and in vivo systems.

EXPERT OPINIONS

We discuss some current challenges of proteomics in iPS cell research and future directions, including the integration of proteomic and transcriptomic data for systems-level analysis.

Systemic proteome adaptions to 7-day complete caloric restriction in humans

Surviving long periods without food has shaped human evolution. In ancient and modern societies, prolonged fasting was/is practiced by billions of people globally for religious purposes, used to treat diseases such as epilepsy, and recently gained popularity as weight loss intervention, but we still have a very limited understanding of the systemic adaptions in humans to extreme caloric restriction of different durations. Here we show that a 7-day water-only fast leads to an average weight loss of 5.7 kg (±0.8 kg) among 12 volunteers (5 women, 7 men). We demonstrate nine distinct proteomic response profiles, with systemic changes evident only after 3 days of complete calorie restriction based on in-depth characterization of the temporal trajectories of ~3,000 plasma proteins measured before, daily during, and after fasting. The multi-organ response to complete caloric restriction shows distinct effects of fasting duration and weight loss and is remarkably conserved across volunteers with >1,000 significantly responding proteins. The fasting signature is strongly enriched for extracellular matrix proteins from various body sites, demonstrating profound non-metabolic adaptions, including extreme changes in the brain-specific extracellular matrix protein tenascin-R. Using proteogenomic approaches, we estimate the health consequences for 212 proteins that change during fasting across ~500 outcomes and identified putative beneficial (SWAP70 and rheumatoid arthritis or HYOU1 and heart disease), as well as adverse effects. Our results advance our understanding of prolonged fasting in humans beyond a merely energy-centric adaptions towards a systemic response that can inform targeted therapeutic modulation.

Cross-Platform Comparison of Highly Sensitive Immunoassays for Inflammatory Markers in a COVID-19 Cohort

A variety of commercial platforms are available for the simultaneous detection of multiple cytokines and associated proteins, often employing Ab pairs to capture and detect target proteins. In this study, we comprehensively evaluated the performance of three distinct platforms: the fluorescent bead-based Luminex assay, the proximity extension-based Olink assay, and a novel proximity ligation assay platform known as Alamar NULISAseq. These assessments were conducted on human serum samples from the National Institutes of Health IMPACC study, with a focus on three essential performance metrics: detectability, correlation, and differential expression. Our results reveal several key findings. First, the Alamar platform demonstrated the highest overall detectability, followed by Olink and then Luminex. Second, the correlation of protein measurements between the Alamar and Olink platforms tended to be stronger than the correlation of either of these platforms with Luminex. Third, we observed that detectability differences across the platforms often translated to differences in differential expression findings, although high detectability did not guarantee the ability to identify meaningful biological differences. Our study provides valuable insights into the comparative performance of these assays, enhancing our understanding of their strengths and limitations when assessing complex biological samples, as exemplified by the sera from this COVID-19 cohort.

Decoding the glycoproteome: a new frontier for biomarker discovery in cancer

Cancer early detection and treatment response prediction continue to pose significant challenges. Cancer liquid biopsies focusing on detecting circulating tumor cells (CTCs) and DNA (ctDNA) have shown enormous potential due to their non-invasive nature and the implications in precision cancer management. Recently, liquid biopsy has been further expanded to profile glycoproteins, which are the products of post-translational modifications of proteins and play key roles in both normal and pathological processes, including cancers. The advancements in chemical and mass spectrometry-based technologies and artificial intelligence-based platforms have enabled extensive studies of cancer and organ-specific changes in glycans and glycoproteins through glycomics and glycoproteomics. Glycoproteomic analysis has emerged as a promising tool for biomarker discovery and development in early detection of cancers and prediction of treatment efficacy including response to immunotherapies. These biomarkers could play a crucial role in aiding in early intervention and personalized therapy decisions. In this review, we summarize the significant advance in cancer glycoproteomic biomarker studies and the promise and challenges in integration into clinical practice to improve cancer patient care.

Large-scale circulating proteome association study (CPAS) meta-analysis identifies circulating proteins and pathways predicting incident hip fractures

Hip fractures are associated with significant disability, high cost, and mortality. However, the exact biological mechanisms underlying susceptibility to hip fractures remain incompletely understood. In an exploratory search of the underlying biology as reflected through the circulating proteome, we performed a comprehensive Circulating Proteome Association Study (CPAS) meta-analysis for incident hip fractures. Analyses included 6430 subjects from two prospective cohort studies (Cardiovascular Health Study and Trøndelag Health Study) with circulating proteomics data (aptamer-based 5 K SomaScan version 4.0 assay; 4979 aptamers). Associations between circulating protein levels and incident hip fractures were estimated for each cohort using age and sex-adjusted Cox regression models. Participants experienced 643 incident hip fractures. Compared with the individual studies, inverse-variance weighted meta-analyses yielded more statistically significant associations, identifying 23 aptamers associated with incident hip fractures (conservative Bonferroni correction 0.05/4979, P < 1.0 × 10-5). The aptamers most strongly associated with hip fracture risk corresponded to two proteins of the growth hormone/insulin growth factor system (GHR and IGFBP2), as well as GDF15 and EGFR. High levels of several inflammation-related proteins (CD14, CXCL12, MMP12, ITIH3) were also associated with increased hip fracture risk. Ingenuity pathway analysis identified reduced LXR/RXR activation and increased acute phase response signaling to be overrepresented among those proteins associated with increased hip fracture risk. These analyses identified several circulating proteins and pathways consistently associated with incident hip fractures. These findings underscore the usefulness of the meta-analytic approach for comprehensive CPAS in a similar manner as has previously been observed for large-scale human genetic studies. Future studies should investigate the underlying biology of these potential novel drug targets.

Exercise-Dependent Modulation of Immunological Response Pathways in Endurance Athletes With and Without Atrial Fibrillation

BACKGROUND

Atrial fibrillation (AF) is a common arrhythmia characterized by uncoordinated atrial electrical activity. Lone AF occurs in the absence of traditional risk factors and is frequently observed in male endurance athletes, who face a 2- to 5-fold higher risk of AF compared with healthy, moderately active males. Our understanding of how endurance exercise contributes to the pathophysiology of lone AF remains limited. This study aimed to characterize the circulating protein fluctuations during high-intensity exercise as well as explore potential biomarkers of exercise-associated AF.

METHODS AND RESULTS

A prospective cohort of 12 male endurance cyclists between the ages of 40 and 65 years, 6 of whom had a history of exercise-associated AF, were recruited to participate using a convenience sampling method. The circulating proteome was subsequently analyzed using multiplex immunoassays and aptamer-based proteomics before, during, and after an acute high-intensity endurance exercise bout to assess temporality and identify potential markers of AF. The endurance exercise bout resulted in significant alterations to proteins involved in immune modulation (eg, growth/differentiation factor 15), skeletal muscle metabolism (eg, α-actinin-2), cell death (eg, histones), and inflammation (eg, interleukin-6). Subjects with AF differed from those without, displaying modulation of proteins previously known to have associations with incident AF (eg, C-reactive protein, insulin-like growth factor-1, and angiopoietin-2), and also with proteins having no previous association (eg, tapasin-related protein and α2-Heremans-Schmid glycoprotein).

CONCLUSIONS

These findings provide insights into the proteomic response to acute intense exercise, provide mechanistic insights into the pathophysiology behind AF in athletes, and identify targets for future study and validation.

Proteomic profiles of cytokines and chemokines in moderate to severe depression: Implications for comorbidities and biomarker discovery

Objective

This study assessed the proteomic profiles of cytokines and chemokines in individuals with moderate to severe depression, with or without comorbid medical disorders, compared to healthy controls. Two proteomic multiplex platforms were employed for this purpose.

Metods

An immunofluorescent multiplex platform and an aptamer-based method were used to evaluate 32 protein analytes from 153 individuals with moderate to severe major depressive disorder (MDD) and healthy controls (HCs). The study focused on determining the level of agreement between the two platforms and evaluating the ability of individual analytes and principal components (PCs) to differentiate between the MDD and HC groups. Additionally, the study investigated the relationship between PCs consisting of chemokines and cytokines and comorbid inflammatory and cardiometabolic diseases.

Findings

Analysis revealed a small or moderate correlation between 47% of the analytes measured by the two platforms. Two proteomic profiles were identified that differentiated individuals with moderate to severe MDD from HCs. High eotaxin, age, BMI, IP-10, or IL-10 characterized profile 1. This profile was associated with several cardiometabolic risk factors, including hypertension, hyperlipidemia, and type 2 diabetes. Profile 2 is characterized by higher age, BMI, interleukins, and a strong negative loading for eotaxin. This profile was associated with inflammation but not cardiometabolic risk factors.

Conclusion

This study provides further evidence that proteomic profiles can be used to identify potential biomarkers and pathways associated with MDD and comorbidities. Our findings suggest that MDD is associated with distinct profiles of proteins that are also associated with cardiometabolic risk factors, inflammation, and obesity. In particular, the chemokines eotaxin and IP-10 appear to play a role in the relationship between MDD and cardiometabolic risk factors. These findings suggest that a focus on the interplay between MDD and comorbidities may be useful in identifying potential targets for intervention and improving overall health outcomes.

Proteome-wide Mendelian randomization identifies causal plasma proteins in lung cancer

Plasma proteins are promising biomarkers and potential drug targets in lung cancer. To evaluate the causal association between plasma proteins and lung cancer, we performed proteome-wide Mendelian randomization meta-analysis (PW-MR-meta) based on lung cancer genome-wide association studies (GWASs), protein quantitative trait loci (pQTLs) of 4,719 plasma proteins in deCODE and 4,775 in Fenland. Further, causal-protein risk score (CPRS) was developed based on causal proteins and validated in the UK Biobank. 270 plasma proteins were identified using PW-MR meta-analysis, including 39 robust causal proteins (both FDR-q < 0.05) and 78 moderate causal proteins (FDR-q < 0.05 in one and p < 0.05 in another). The CPRS had satisfactory performance in risk stratification for lung cancer (top 10% CPRS:Hazard ratio (HR) (95%CI):4.33(2.65-7.06)). The CPRS [AUC (95%CI): 65.93 (62.91-68.78)] outperformed the traditional polygenic risk score (PRS) [AUC (95%CI): 55.71(52.67-58.59)]. Our findings offer further insight into the genetic architecture of plasma proteins for lung cancer susceptibility.

Nanoparticle enrichment mass-spectrometry proteomics identifies protein-altering variants for precise pQTL mapping

Proteogenomics studies generate hypotheses on protein function and provide genetic evidence for drug target prioritization. Most previous work has been conducted using affinity-based proteomics approaches. These technologies face challenges, such as uncertainty regarding target identity, non-specific binding, and handling of variants that affect epitope affinity binding. Mass spectrometry-based proteomics can overcome some of these challenges. Here we report a pQTL study using the Proteograph™ Product Suite workflow (Seer, Inc.) where we quantify over 18,000 unique peptides from nearly 3000 proteins in more than 320 blood samples from a multi-ethnic cohort in a bottom-up, peptide-centric, mass spectrometry-based proteomics approach. We identify 184 protein-altering variants in 137 genes that are significantly associated with their corresponding variant peptides, confirming target specificity of co-associated affinity binders, identifying putatively causal cis-encoded proteins and providing experimental evidence for their presence in blood, including proteins that may be inaccessible to affinity-based proteomics.

Circadian protein expression patterns in healthy young adults

OBJECTIVES

To explore how the blood plasma proteome fluctuates across the 24-hour day and identify a subset of proteins that show endogenous circadian rhythmicity.

METHODS

Plasma samples from 17 healthy adults were collected hourly under controlled conditions designed to unmask endogenous circadian rhythmicity; in a subset of 8 participants, we also collected samples across a day on a typical sleep-wake schedule. A total of 6916 proteins were analyzed from each sample using the SomaScan aptamer-based multiplexed platform. We used differential rhythmicity analysis based on a cosinor model with mixed effects to identify a subset of proteins that showed circadian rhythmicity in their abundance.

RESULTS

One thousand and sixty-three (15%) proteins exhibited significant daily rhythmicity. Of those, 431 (6.2%) proteins displayed consistent endogenous circadian rhythms on both a sleep-wake schedule and under controlled conditions: it included both known and novel proteins. When models were fitted with two harmonics, an additional 259 (3.7%) proteins exhibited significant endogenous circadian rhythmicity, indicating that some rhythmic proteins cannot be solely captured by a simple sinusoidal model. Overall, we found that the largest number of proteins had their peak levels in the late afternoon/evening, with another smaller group peaking in the early morning.

CONCLUSIONS

This study reveals that hundreds of plasma proteins exhibit endogenous circadian rhythmicity in humans. Future analyses will likely reveal novel physiological pathways regulated by circadian clocks and pave the way for improved diagnosis and treatment for patients with circadian disorders and other pathologies. It will also advance efforts to include knowledge about time-of-day, thereby incorporating circadian medicine into personalized medicine.

Association between circulating inflammatory markers and adult cancer risk: a Mendelian randomization analysis

BACKGROUND

Tumour-promoting inflammation is a "hallmark" of cancer and conventional epidemiological studies have reported links between various inflammatory markers and cancer risk. The causal nature of these relationships and, thus, the suitability of these markers as intervention targets for cancer prevention is unclear.

METHODS

We meta-analysed 6 genome-wide association studies of circulating inflammatory markers comprising 59,969 participants of European ancestry. We then used combined cis-Mendelian randomization and colocalisation analysis to evaluate the causal role of 66 circulating inflammatory markers in risk of 30 adult cancers in 338,294 cancer cases and up to 1,238,345 controls. Genetic instruments for inflammatory markers were constructed using genome-wide significant (P < 5.0 × 10-8) cis-acting SNPs (i.e., in or ±250 kb from the gene encoding the relevant protein) in weak linkage disequilibrium (LD, r2 < 0.10). Effect estimates were generated using inverse-variance weighted random-effects models and standard errors were inflated to account for weak LD between variants with reference to the 1000 Genomes Phase 3 CEU panel. A false discovery rate (FDR)-corrected P-value ("q-value") <0.05 was used as a threshold to define "strong evidence" to support associations and 0.05 ≤ q-value < 0.20 to define "suggestive evidence". A colocalisation posterior probability (PPH4) >70% was employed to indicate support for shared causal variants across inflammatory markers and cancer outcomes. Findings were replicated in the FinnGen study and then pooled using meta-analysis.

FINDINGS

We found strong evidence to support an association of genetically-proxied circulating pro-adrenomedullin concentrations with increased breast cancer risk (OR: 1.19, 95% CI: 1.10-1.29, q-value = 0.033, PPH4 = 84.3%) and suggestive evidence to support associations of interleukin-23 receptor concentrations with increased pancreatic cancer risk (OR: 1.42, 95% CI: 1.20-1.69, q-value = 0.055, PPH4 = 73.9%), prothrombin concentrations with decreased basal cell carcinoma risk (OR: 0.66, 95% CI: 0.53-0.81, q-value = 0.067, PPH4 = 81.8%), and interleukin-1 receptor-like 1 concentrations with decreased triple-negative breast cancer risk (OR: 0.92, 95% CI: 0.88-0.97, q-value = 0.15, PPH4 = 85.6%). These findings were replicated in pooled analyses with the FinnGen study. Though suggestive evidence was found to support an association of macrophage migration inhibitory factor concentrations with increased bladder cancer risk (OR: 2.46, 95% CI: 1.48-4.10, q-value = 0.072, PPH4 = 76.1%), this finding was not replicated when pooled with the FinnGen study. For 22 of 30 cancer outcomes examined, there was little evidence (q-value ≥0.20) that any of the 66 circulating inflammatory markers examined were associated with cancer risk.

INTERPRETATION

Our comprehensive joint Mendelian randomization and colocalisation analysis of the role of circulating inflammatory markers in cancer risk identified potential roles for 4 circulating inflammatory markers in risk of 4 site-specific cancers. Contrary to reports from some prior conventional epidemiological studies, we found little evidence of association of circulating inflammatory markers with the majority of site-specific cancers evaluated.

FUNDING

Cancer Research UK (C68933/A28534, C18281/A29019, PPRCPJT∖100005), World Cancer Research Fund (IIG_FULL_2020_022), National Institute for Health Research (NIHR202411, BRC-1215-20011), Medical Research Council (MC_UU_00011/1, MC_UU_00011/3, MC_UU_00011/6, and MC_UU_00011/4), Academy of Finland Project 326291, European Union's Horizon 2020 grant agreement no. 848158 (EarlyCause), French National Cancer Institute (INCa SHSESP20, 2020-076), Versus Arthritis (21173, 21754, 21755), National Institutes of Health (U19 CA203654), National Cancer Institute (U19CA203654).

Proteomics for heart failure risk stratification: a systematic review

BACKGROUND

Heart failure (HF) is a complex clinical syndrome with persistently high mortality. High-throughput proteomic technologies offer new opportunities to improve HF risk stratification, but their contribution remains to be clearly defined. We aimed to systematically review prognostic studies using high-throughput proteomics to identify protein signatures associated with HF mortality.

METHODS

We searched four databases and two clinical trial registries for articles published from 2012 to 2023. HF proteomics studies measuring high numbers of proteins using aptamer or antibody-based affinity platforms on human plasma or serum with outcomes of all-cause or cardiovascular death were included. Two reviewers independently screened articles, extracted data, and assessed the risk of bias. A third reviewer resolved conflicts. We assessed the risk of bias using the Risk Of Bias In Non-randomized Studies-of Exposure tool.

RESULTS

Out of 5131 unique articles identified, nine articles were included in the review. The nine studies were observational; three used the aptamer platform, and six used the antibody platform. We found considerable heterogeneity across studies in measurement panels, HF definitions, ejection fraction categorization, follow-up duration, and outcome definitions, and a lack of risk estimates for most protein associations. Hence, we proceeded with a systematic review rather than a meta-analysis. In two comparable aptamer studies in patients with HF with reduced ejection fraction, 21 proteins were identified in common for the association with all-cause death. Among these, one protein, WAP four-disulfide core domain protein 2 was also reported in an antibody study on HFrEF and for the association with CV death. We proposed standardized reporting criteria to facilitate the interpretation of future studies.

CONCLUSIONS

In this systematic review of nine studies evaluating the association of proteomics with mortality in HF, we identified a limited number of proteins common across several studies. Heterogeneity across studies compromised drawing broad inferences, underscoring the importance of standardized approaches to reporting.

Longitudinal plasma proteomics reveals biomarkers of alveolar-capillary barrier disruption in critically ill COVID-19 patients

The pathobiology of respiratory failure in COVID-19 consists of a complex interplay between viral cytopathic effects and a dysregulated host immune response. In critically ill patients, imatinib treatment demonstrated potential for reducing invasive ventilation duration and mortality. Here, we perform longitudinal profiling of 6385 plasma proteins in 318 hospitalised patients to investigate the biological processes involved in critical COVID-19, and assess the effects of imatinib treatment. Nine proteins measured at hospital admission accurately predict critical illness development. Next to dysregulation of inflammation, critical illness is characterised by pathways involving cellular adhesion, extracellular matrix turnover and tissue remodelling. Imatinib treatment attenuates protein perturbations associated with inflammation and extracellular matrix turnover. These proteomic alterations are contextualised using external pulmonary RNA-sequencing data of deceased COVID-19 patients and imatinib-treated Syrian hamsters. Together, we show that alveolar capillary barrier disruption in critical COVID-19 is reflected in the plasma proteome, and is attenuated with imatinib treatment. This study comprises a secondary analysis of both clinical data and plasma samples derived from a clinical trial that was registered with the EU Clinical Trials Register (EudraCT 2020-001236-10, https://www.clinicaltrialsregister.eu/ctr-search/trial/2020-001236-10/NL ) and Netherlands Trial Register (NL8491, https://www.trialregister.nl/trial/8491 ).

Proteome Analysis of Serum Purified Using Solanum tuberosum and Lycopersicon esculentum Lectins

Serum and plasma exhibit a broad dynamic range of protein concentrations, posing challenges for proteome analysis. Various technologies have been developed to reduce this complexity, including high-abundance depletion methods utilizing antibody columns, extracellular vesicle enrichment techniques, and trace protein enrichment using nanobead cocktails. Here, we employed lectins to address this, thereby extending the scope of biomarker discovery in serum or plasma using a novel approach. We enriched serum proteins using 37 different lectins and subjected them to LC-MS/MS analysis with data-independent acquisition. Solanum tuberosum lectin (STL) and Lycopersicon esculentum lectin (LEL) enabled the detection of more serum proteins than the other lectins. STL and LEL bind to N-acetylglucosamine oligomers, emphasizing the significance of capturing these oligomer-binding proteins when analyzing serum trace proteins. Combining STL and LEL proved more effective than using them separately, allowing us to identify over 3000 proteins from serum through single-shot proteome analysis. We applied the STL/LEL trace-protein enrichment method to the sera of systemic lupus erythematosus model mice. This revealed differences in >1300 proteins between the systemic lupus erythematosus model and control mouse sera, underscoring the utility of this method for biomarker discovery.

Analytical validation of the PROphet test for treatment decision-making guidance in metastatic non-small cell lung cancer

The blood proteome, consisting of thousands of proteins engaged in various biological processes, acts as a valuable source of potential biomarkers for various medical applications. PROphet is a plasma proteomics-based test that serves as a decision-support tool for non-small cell lung cancer (NSCLC) patients, combining proteomic profiling using SomaScan technology and subsequent computational algorithm. PROphet was implemented as a laboratory developed test (LDT). Under the Clinical Laboratory Improvement Amendments (CLIA) and Commission on Office Laboratory Accreditation (COLA) regulations, prior to releasing patient test results, a clinical laboratory located in the United States employing an LDT must examine its performance characteristics with regard to analytical validity. This study describes the experimental and computational analytical validity of the PROphet test, as required by CLIA/COLA regulations. Experimental precision analysis displayed a median coefficient of variation (CV) of 3.9 % and 4.7 % for intra-plate and inter-plate examination, respectively, and the median accuracy rate between sites was 88 %. Computational precision exhibited a high accuracy rate, with 93 % of samples displaying complete concordance in results. A cross-platform comparison between SomaScan and other proteomics platforms yielded a median Spearman's rank correlation coefficient of 0.51, affirming the consistency and reliability of the SomaScan platform as used under the PROphet test. Our study presents a robust framework for evaluating the analytical validity of a platform that combines an experimental assay with subsequent computational algorithms. When applied to the PROphet test, strong analytical performance of the test was demonstrated.

Associations of Baseline and Longitudinal Serum Uromodulin With Kidney Failure and Mortality: Results From the African American Study of Kidney Disease and Hypertension (AASK) Trial

RATIONALE & OBJECTIVE

Uromodulin (UMOD) is the most abundant protein found in urine and has emerged as a promising biomarker of tubule health. Circulating UMOD is also detectable, but at lower levels. We evaluated whether serum UMOD levels were associated with the risks of incident kidney failure with replacement therapy (KFRT) and mortality.

STUDY DESIGN

Prospective cohort.

SETTING & PARTICIPANTS

Participants in AASK (the African American Study of Kidney Disease and Hypertension) with available stored serum samples from the 0-, 12-, and 24-month visits for biomarker measurement.

PREDICTORS

Baseline log-transformed UMOD and change in UMOD over 2 years.

OUTCOMES

KFRT and mortality.

ANALYTICAL APPROACH

Cox proportional hazards and mixed-effects models.

RESULTS

Among 500 participants with baseline serum UMOD levels (mean age, 54y; 37% female), 161 KFRT events occurred during a median of 8.5 years. After adjusting for baseline demographic factors, clinical factors, glomerular filtration rate, log-transformed urine protein-creatinine ratio, and randomized treatment groups, a 50% lower baseline UMOD level was independently associated with a 35% higher risk of KFRT (adjusted HR, 1.35; 95% CI, 1.07-1.70). For annual UMOD change, each 1-standard deviation lower change was associated with a 67% higher risk of KFRT (adjusted HR, 1.67; 95% CI, 1.41-1.99). Baseline UMOD and UMOD change were not associated with mortality. UMOD levels declined more steeply for metoprolol versus ramipril (P<0.001) as well as for intensive versus standard blood pressure goals (P = 0.002).

LIMITATIONS

Small sample size and limited generalizability.

CONCLUSIONS

Lower UMOD levels at baseline and steeper declines in UMOD over time were associated with a higher risk of subsequent KFRT in a cohort of African American adults with chronic kidney disease and hypertension.

PLAIN-LANGUAGE SUMMARY

Prior studies of uromodulin (UMOD), the most abundant protein in urine, and kidney disease have focused primarily on urinary UMOD levels. The present study evaluated associations of serum UMOD levels with the risks of kidney failure with replacement therapy (KFRT) and mortality in a cohort of African American adults with hypertension and chronic kidney disease. It found that participants with lower levels of UMOD at baseline were more likely to experience KFRT even after accounting for baseline kidney measures. Similarly, participants who experienced steeper annual declines in UMOD also had a heightened risk of kidney failure. Neither baseline nor annual change in UMOD was associated with mortality. Serum UMOD is a promising biomarker of kidney health.

Lower myostatin and higher MUC1 levels are associated with better response to mepolizumab and omalizumab in asthma: a protein-protein interaction analyses

INTRODUCTION

Biomarkers are needed to inform the choice of biologic therapy in patients with asthma given the increasing number of biologics. We aimed to identify proteins associated with response to omalizumab and mepolizumab.

METHODS

Aptamer-based proteomic profiling (SomaScan) was used to assess 1437 proteins from 51 patients with moderate to severe asthma who received omalizumab (n = 29) or mepolizumab (n = 22). Response was defined as the change in asthma-related exacerbations in the 12 months following therapy initiation. All models were adjusted for age, sex, and pre-treatment exacerbation rate. Additionally, body mass index was included in the omalizumab model and eosinophil count in the mepolizumab model. We evaluated the association between molecular signatures and response using negative binomial regression correcting for the false discovery rate (FDR) and gene set enrichment analyses (GSEA) to identify associated pathways.

RESULTS

Over two-thirds of patients were female. The average age for omalizumab patients was 42 years and 57 years for mepolizumab. At baseline, the average exacerbation rate was 1.5/year for omalizumab and 2.4/year for mepolizumab. Lower levels of LOXL2 (unadjusted p: 1.93 × 10E-05, FDR-corrected: 0.028) and myostatin (unadjusted: 3.87 × 10E-05, FDR-corrected: 0.028) were associated with better response to mepolizumab. Higher levels of CD9 antigen (unadjusted: 5.30 × 10E-07, FDR-corrected: 0.0006) and MUC1 (unadjusted: 1.15 × 10E-06, FDR-corrected: 0.0006) were associated with better response to omalizumab, and LTB4R (unadjusted: 1.12 × 10E-06, FDR-corrected: 0.0006) with worse response. Protein-protein interaction network modeling showed an enrichment of the TNF- and NF-kB signaling pathways for patients treated with mepolizumab and multiple pathways involving MAPK, including the FcER1 pathway, for patients treated with omalizumab.

CONCLUSIONS

This study provides novel fundamental data on proteins associated with response to mepolizumab or omalizumab in severe asthma and warrants further validation as potential biomarkers for therapy selection.

Evaluation of circulating plasma proteins in breast cancer using Mendelian randomisation

Biomarkers for early detection of breast cancer may complement population screening approaches to enable earlier and more precise treatment. The blood proteome is an important source for biomarker discovery but so far, few proteins have been identified with breast cancer risk. Here, we measure 2929 unique proteins in plasma from 598 women selected from the Karolinska Mammography Project to explore the association between protein levels, clinical characteristics, and gene variants, and to identify proteins with a causal role in breast cancer. We present 812 cis-acting protein quantitative trait loci for 737 proteins which are used as instruments in Mendelian randomisation analyses of breast cancer risk. Of those, we present five proteins (CD160, DNPH1, LAYN, LRRC37A2 and TLR1) that show a potential causal role in breast cancer risk with confirmatory results in independent cohorts. Our study suggests that these proteins should be further explored as biomarkers and potential drug targets in breast cancer.

Identification of circulating proteins associated with general cognitive function among middle-aged and older adults

Identifying circulating proteins associated with cognitive function may point to biomarkers and molecular process of cognitive impairment. Few studies have investigated the association between circulating proteins and cognitive function. We identify 246 protein measures quantified by the SomaScan assay as associated with cognitive function (p < 4.9E-5, n up to 7289). Of these, 45 were replicated using SomaScan data, and three were replicated using Olink data at Bonferroni-corrected significance. Enrichment analysis linked the proteins associated with general cognitive function to cell signaling pathways and synapse architecture. Mendelian randomization analysis implicated higher levels of NECTIN2, a protein mediating viral entry into neuronal cells, with higher Alzheimer's disease (AD) risk (p = 2.5E-26). Levels of 14 other protein measures were implicated as consequences of AD susceptibility (p < 2.0E-4). Proteins implicated as causes or consequences of AD susceptibility may provide new insight into the potential relationship between immunity and AD susceptibility as well as potential therapeutic targets.

Differences in the Circulating Proteome in Individuals with versus without Sickle Cell Trait

BACKGROUND

Sickle cell trait affects approximately 8% of Black individuals in the United States, along with many other individuals with ancestry from malaria-endemic regions worldwide. While traditionally considered a benign condition, recent evidence suggests that sickle cell trait is associated with lower eGFR and higher risk of kidney diseases, including kidney failure. The mechanisms underlying these associations remain poorly understood. We used proteomic profiling to gain insight into the pathobiology of sickle cell trait.

METHODS

We measured proteomics ( N =1285 proteins assayed by Olink Explore) using baseline plasma samples from 592 Black participants with sickle cell trait and 1:1 age-matched Black participants without sickle cell trait from the prospective Women's Health Initiative cohort. Age-adjusted linear regression was used to assess the association between protein levels and sickle cell trait.

RESULTS

In age-adjusted models, 35 proteins were significantly associated with sickle cell trait after correction for multiple testing. Several of the sickle cell trait-protein associations were replicated in Black participants from two independent cohorts (Atherosclerosis Risk in Communities study and Jackson Heart Study) assayed using an orthogonal aptamer-based proteomic platform (SomaScan). Many of the validated sickle cell trait-associated proteins are known biomarkers of kidney function or injury ( e.g. , hepatitis A virus cellular receptor 1 [HAVCR1]/kidney injury molecule-1 [KIM-1], uromodulin [UMOD], ephrins), related to red cell physiology or hemolysis (erythropoietin [EPO], heme oxygenase 1 [HMOX1], and α -hemoglobin stabilizing protein) and/or inflammation (fractalkine, C-C motif chemokine ligand 2/monocyte chemoattractant protein-1 [MCP-1], and urokinase plasminogen activator surface receptor [PLAUR]). A protein risk score constructed from the top sickle cell trait-associated biomarkers was associated with incident kidney failure among those with sickle cell trait during Women's Health Initiative follow-up (odds ratio, 1.32; 95% confidence interval, 1.10 to 1.58).

CONCLUSIONS

We identified and replicated the association of sickle cell trait with a number of plasma proteins related to hemolysis, kidney injury, and inflammation.

Cross-platform comparison of highly-sensitive immunoassays for inflammatory markers in a COVID-19 cohort 1

A variety of commercial platforms are available for the simultaneous detection of multiple cytokines and associated proteins, often employing antibody pairs to capture and detect target proteins. In this study, we comprehensively evaluated the performance of three distinct platforms: the fluorescent bead-based Luminex assay, the proximity extension-based Olink assay, and a novel proximity ligation assay platform known as Alamar NULISAseq. These assessments were conducted on serum samples from the NIH IMPACC study, with a focus on three essential performance metrics: detectability, correlation, and differential expression. Our results reveal several key findings. Firstly, the Alamar platform demonstrated the highest overall detectability, followed by Olink and then Luminex. Secondly, the correlation of protein measurements between the Alamar and Olink platforms tended to be stronger than the correlation of either of these platforms with Luminex. Thirdly, we observed that detectability differences across the platforms often translated to differences in differential expression findings, although high detectability did not guarantee the ability to identify meaningful biological differences. Our study provides valuable insights into the comparative performance of these assays, enhancing our understanding of their strengths and limitations when assessing complex biological samples, as exemplified by the sera from this COVID-19 cohort.

Liquid-biopsy proteomics combined with AI identifies cellular drivers of eye aging and disease in vivo

Single-cell analysis in living humans is essential for understanding disease mechanisms, but it is impractical in non-regenerative organs, such as the eye and brain, because tissue biopsies would cause serious damage. We resolve this problem by integrating proteomics of liquid biopsies with single-cell transcriptomics from all known ocular cell types to trace the cellular origin of 5,953 proteins detected in the aqueous humor. We identified hundreds of cell-specific protein markers, including for individual retinal cell types. Surprisingly, our results reveal that retinal degeneration occurs in Parkinson's disease, and the cells driving diabetic retinopathy switch with disease stage. Finally, we developed artificial intelligence (AI) models to assess individual cellular aging and found that many eye diseases not associated with chronological age undergo accelerated molecular aging of disease-specific cell types. Our approach, which can be applied to other organ systems, has the potential to transform molecular diagnostics and prognostics while uncovering new cellular disease and aging mechanisms.

Time-varying and tissue-dependent effects of adiposity on leptin levels: A Mendelian randomization study

Background

Findings from Mendelian randomization (MR) studies are conventionally interpreted as lifelong effects, which typically do not provide insight into the molecular mechanisms underlying the effect of an exposure on an outcome. In this study, we apply two recently developed MR approaches (known as 'lifecourse' and 'tissue-partitioned' MR) to investigate lifestage-specific effects and tissues of action in the relationship between adiposity and circulating leptin levels.

Methods

Genetic instruments for childhood and adult adiposity were incorporated into a multivariable MR (MVMR) framework to estimate lifestage-specific effects on leptin levels measured during early life (mean age: 10 y) in the Avon Longitudinal Study of Parents and Children and in adulthood (mean age: 55 y) using summary-level data from the deCODE Health study. This was followed by partitioning body mass index (BMI) instruments into those whose effects are putatively mediated by gene expression in either subcutaneous adipose or brain tissues, followed by using MVMR to simultaneously estimate their separate effects on childhood and adult leptin levels.

Results

There was strong evidence that childhood adiposity has a direct effect on leptin levels at age 10 y in the lifecourse (β = 1.10 SD change in leptin levels, 95% CI = 0.90-1.30, p=6 × 10-28), whereas evidence of an indirect effect was found on adulthood leptin along the causal pathway involving adulthood body size (β = 0.74, 95% CI = 0.62-0.86, p=1 × 10-33). Tissue-partitioned MR analyses provided evidence to suggest that BMI exerts its effect on leptin levels during both childhood and adulthood via brain tissue-mediated pathways (β = 0.79, 95% CI = 0.22-1.36, p=6 × 10-3 and β = 0.51, 95% CI = 0.32-0.69, p=1 × 10-7, respectively).

Conclusions

Our findings demonstrate the use of lifecourse MR to disentangle direct and indirect effects of early-life exposures on time-varying complex outcomes. Furthermore, by integrating tissue-specific data, we highlight the etiological importance of appetite regulation in the effect of adiposity on leptin levels.

Funding

This work was supported by the Integrative Epidemiology Unit, which receives funding from the UK Medical Research Council and the University of Bristol (MC_UU_00011/1).

High-Throughput CSF Proteomics and Machine Learning to Identify Proteomic Signatures for Parkinson Disease Development and Progression

BACKGROUND AND OBJECTIVES

This study aimed to identify CSF proteomic signatures characteristic of Parkinson disease (PD) and evaluate their clinical utility.

METHODS

This observational study used data from the Parkinson's Progression Markers Initiative (PPMI), which enrolled patients with PD, healthy controls (HCs), and non-PD participants carrying GBA1, LRRK2, and/or SNCA pathogenic variants (genetic prodromals) at international sites. Study participants were chosen from PPMI enrollees based on the availability of aptamer-based CSF proteomic data, quantifying 4,071 proteins, and classified as patients with PD without GBA1, LRRK2, and/or SNCA pathogenic variants (nongenetic PD), HCs, patients with PD carrying the aforementioned pathogenic variants (genetic PD), or genetic prodromals. Differentially expressed protein (DEP) analysis and the least absolute shrinkage and selection operator (LASSO) were applied to the data from nongenetic PD and HCs. Signatures characteristics of nongenetic PD were quantified as a PD proteomic score (PD-ProS), validated internally and then externally using data of 1,556 CSF proteins from the LRRK2 Cohort Consortium (LCC). We further tested the PD-ProS in genetic PD and genetic prodromals and examined associations with clinical progression.

RESULTS

Data from 279 patients with nongenetic PD (mean ± SD, age 62.0 ± 9.6 years; male 67.7%) and 141 HCs (age 60.5 ± 11.9 years; male 64.5%) were used for PD-ProS derivation. From 23 DEPs, LASSO determined weights of 14 DEPs for the PD-ProS (area under the curve [AUC] 0.83, 95% CI 0.78-0.87), validated in an independent internal validation cohort of 71 patients with nongenetic PD and 35 HCs (AUC 0.81, 95% CI 0.73-0.90). In the LCC, only 5 of the 14 DEPs were also measured. Notably, these 5 DEPs still distinguished 34 patients with nongenetic PD from 31 HCs with the same weights (AUC 0.75, 95% CI 0.63-0.87). Furthermore, the PD-ProS distinguished 258 patients with genetic PD from 365 genetic prodromals. Finally, regardless of genetic status, the PD-ProS independently predicted both cognitive and motor decline in PD (dementia, adjusted hazard ratio in the highest quintile [aHR-Q5] 2.8 [95% CI 1.6-5.0]; Hoehn and Yahr stage IV, aHR-Q5 2.1 [95% CI 1.1-4.0]).

DISCUSSION

By integrating high-throughput proteomics with machine learning, we identified PD-associated CSF proteomic signatures crucial for PD development and progression.

TRIAL REGISTRATION INFORMATION

ClinicalTrials.gov (NCT01176565). A link to the trial registry page is clinicaltrials.gov/ct2/show/NCT01141023.

CLASSIFICATION OF EVIDENCE

This study provides Class II evidence that the CSF proteome contains clinically important information regarding the development and progression of Parkinson disease that can be deciphered by a combination of high-throughput proteomics and machine learning.

Large-scale plasma proteomics comparisons through genetics and disease associations

High-throughput proteomics platforms measuring thousands of proteins in plasma combined with genomic and phenotypic information have the power to bridge the gap between the genome and diseases. Here we performed association studies of Olink Explore 3072 data generated by the UK Biobank Pharma Proteomics Project1 on plasma samples from more than 50,000 UK Biobank participants with phenotypic and genotypic data, stratifying on British or Irish, African and South Asian ancestries. We compared the results with those of a SomaScan v4 study on plasma from 36,000 Icelandic people2, for 1,514 of whom Olink data were also available. We found modest correlation between the two platforms. Although cis protein quantitative trait loci were detected for a similar absolute number of assays on the two platforms (2,101 on Olink versus 2,120 on SomaScan), the proportion of assays with such supporting evidence for assay performance was higher on the Olink platform (72% versus 43%). A considerable number of proteins had genomic associations that differed between the platforms. We provide examples where differences between platforms may influence conclusions drawn from the integration of protein levels with the study of diseases. We demonstrate how leveraging the diverse ancestries of participants in the UK Biobank helps to detect novel associations and refine genomic location. Our results show the value of the information provided by the two most commonly used high-throughput proteomics platforms and demonstrate the differences between them that at times provides useful complementarity.

Plasma proteomic associations with genetics and health in the UK Biobank

The Pharma Proteomics Project is a precompetitive biopharmaceutical consortium characterizing the plasma proteomic profiles of 54,219 UK Biobank participants. Here we provide a detailed summary of this initiative, including technical and biological validations, insights into proteomic disease signatures, and prediction modelling for various demographic and health indicators. We present comprehensive protein quantitative trait locus (pQTL) mapping of 2,923 proteins that identifies 14,287 primary genetic associations, of which 81% are previously undescribed, alongside ancestry-specific pQTL mapping in non-European individuals. The study provides an updated characterization of the genetic architecture of the plasma proteome, contextualized with projected pQTL discovery rates as sample sizes and proteomic assay coverages increase over time. We offer extensive insights into trans pQTLs across multiple biological domains, highlight genetic influences on ligand-receptor interactions and pathway perturbations across a diverse collection of cytokines and complement networks, and illustrate long-range epistatic effects of ABO blood group and FUT2 secretor status on proteins with gastrointestinal tissue-enriched expression. We demonstrate the utility of these data for drug discovery by extending the genetic proxied effects of protein targets, such as PCSK9, on additional endpoints, and disentangle specific genes and proteins perturbed at loci associated with COVID-19 susceptibility. This public-private partnership provides the scientific community with an open-access proteomics resource of considerable breadth and depth to help to elucidate the biological mechanisms underlying proteo-genomic discoveries and accelerate the development of biomarkers, predictive models and therapeutics1.

oFVSD: a Python package of optimized forward variable selection decoder for high-dimensional neuroimaging data

The complexity and high dimensionality of neuroimaging data pose problems for decoding information with machine learning (ML) models because the number of features is often much larger than the number of observations. Feature selection is one of the crucial steps for determining meaningful target features in decoding; however, optimizing the feature selection from such high-dimensional neuroimaging data has been challenging using conventional ML models. Here, we introduce an efficient and high-performance decoding package incorporating a forward variable selection (FVS) algorithm and hyper-parameter optimization that automatically identifies the best feature pairs for both classification and regression models, where a total of 18 ML models are implemented by default. First, the FVS algorithm evaluates the goodness-of-fit across different models using the k-fold cross-validation step that identifies the best subset of features based on a predefined criterion for each model. Next, the hyperparameters of each ML model are optimized at each forward iteration. Final outputs highlight an optimized number of selected features (brain regions of interest) for each model with its accuracy. Furthermore, the toolbox can be executed in a parallel environment for efficient computation on a typical personal computer. With the optimized forward variable selection decoder (oFVSD) pipeline, we verified the effectiveness of decoding sex classification and age range regression on 1,113 structural magnetic resonance imaging (MRI) datasets. Compared to ML models without the FVS algorithm and with the Boruta algorithm as a variable selection counterpart, we demonstrate that the oFVSD significantly outperformed across all of the ML models over the counterpart models without FVS (approximately 0.20 increase in correlation coefficient, r, with regression models and 8% increase in classification models on average) and with Boruta variable selection algorithm (approximately 0.07 improvement in regression and 4% in classification models). Furthermore, we confirmed the use of parallel computation considerably reduced the computational burden for the high-dimensional MRI data. Altogether, the oFVSD toolbox efficiently and effectively improves the performance of both classification and regression ML models, providing a use case example on MRI datasets. With its flexibility, oFVSD has the potential for many other modalities in neuroimaging. This open-source and freely available Python package makes it a valuable toolbox for research communities seeking improved decoding accuracy.

Mesenchymal Stem Cell-Derived Exosomes: A Novel Approach to Diabetes-Associated Cognitive Impairment

The progression of diabetes frequently results in a myriad of neurological disorders, including ischemic stroke, depression, blood-brain barrier impairment, and cognitive dysfunction. Notably, diabetes-associated cognitive impairment, a prevalent comorbidity during the course of diabetes, progressively affects patients' cognitive abilities and may reciprocally influence diabetes management, thereby severely impacting patients' quality of life. Extracellular vesicles, particularly nanoscale exosomes, have garnered considerable attention in recent years. These exosomes carry and transfer various functional molecules, such as proteins, lipids, and diverse non-coding RNAs, serving as novel regulators and communicators in intercellular interactions. Of particular interest, mesenchymal stem cell-derived exosomes (MSC-Exos) have been reported to traverse the blood-brain barrier and ameliorate intracerebral pathologies. This review elucidates the role of MSC-Exos in diabetes-related cognitive impairment, with a focus on their applications as biomarkers, modulation of neuronal regeneration and synaptic plasticity, anti-inflammatory properties, antioxidative effects, and their involvement in regulating the functionality of β-amyloid proteins during the course of cognitive impairment. The immense therapeutic potential of MSC-Exos in the treatment of diabetes-induced cognitive dysfunction is emphasized.

Plasma proteins and onset of type 2 diabetes and diabetic complications: Proteome-wide Mendelian randomization and colocalization analyses

We conduct proteome-wide Mendelian randomization and colocalization analyses to decipher the associations of blood proteins with the risk of type 2 diabetes and diabetic complications. Genetic data on plasma proteome are obtained from 54,306 UK Biobank participants and 35,559 Icelanders. Summary-level data on type 2 diabetes are obtained from the DIAGRAM (DIAbetes Genetics Replication And Meta-analysis consortium) consortium (74,124 cases) and FinnGen study (33,043 cases). Data on 10 diabetic complications are obtained from FinnGen and corresponding studies. Among 1,886 proteins, genetically predicted levels of 47 plasma proteins are associated with type 2 diabetes. Eleven of these proteins have strong support of colocalization. Seventeen proteins are associated with at least one diabetic complication, although a few have colocalization support. HLA-DRA, AGER, HSPA1A, and HSPA1B are associated with most microvascular complications. This study reveals causal proteins for the onset of type 2 diabetes and diabetic complications, which enhances the understanding of molecular etiology and development of therapeutics.

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.

Pan-cancer proteogenomics connects oncogenic drivers to functional states

Cancer driver events refer to key genetic aberrations that drive oncogenesis; however, their exact molecular mechanisms remain insufficiently understood. Here, our multi-omics pan-cancer analysis uncovers insights into the impacts of cancer drivers by identifying their significant cis-effects and distal trans-effects quantified at the RNA, protein, and phosphoprotein levels. Salient observations include the association of point mutations and copy-number alterations with the rewiring of protein interaction networks, and notably, most cancer genes converge toward similar molecular states denoted by sequence-based kinase activity profiles. A correlation between predicted neoantigen burden and measured T cell infiltration suggests potential vulnerabilities for immunotherapies. Patterns of cancer hallmarks vary by polygenic protein abundance ranging from uniform to heterogeneous. Overall, our work demonstrates the value of comprehensive proteogenomics in understanding the functional states of oncogenic drivers and their links to cancer development, surpassing the limitations of studying individual cancer types.

Plasma Proteomic Kinetics in Response to Acute Exercise

Regular exercise has many favorable effects on human health, which may be mediated in part by the release of circulating bioactive factors during each bout of exercise. Limited data exist regarding the kinetic responses of plasma proteins during and after acute exercise. Proteomic profiling of 4163 proteins was performed using a large-scale, affinity-based platform in 75 middle-aged adults who were referred for treadmill exercise stress testing. Plasma proteins were quantified at baseline, peak exercise, and 1-h postexercise, and those with significant changes at both exercise timepoints were further examined for their associations with cardiometabolic traits and change with aerobic exercise training in the Health, Risk Factors, Exercise Training and Genetics Family Study, a 20-week exercise intervention study. A total of 765 proteins changed (false discovery rate < 0.05) at peak exercise compared to baseline, and 128 proteins changed (false discovery rate < 0.05) at 1-h postexercise. The 56 proteins that changed at both timepoints included midkine, brain-derived neurotrophic factor, metalloproteinase inhibitor 4, and coiled-coil domain-containing protein 126 and were enriched for secreted proteins. The majority had concordant direction of change at both timepoints. Across all proteins assayed, gene set enrichment analysis showed increased abundance of coagulation-related proteins at 1-h postexercise. Forty-five proteins were associated with at least one measure of adiposity, lipids, glucose homeostasis, or cardiorespiratory fitness in Health, Risk Factors, Exercise Training and Genetics Family Study, and 20 proteins changed with aerobic exercise training. We identified hundreds of novel proteins that change during acute exercise, most of which resolved by 1 h into recovery. Proteins with sustained changes during exercise and recovery may be of particular interest as circulating biomarkers and pathways for further investigation in cardiometabolic diseases. These data will contribute to a biochemical roadmap of acute exercise that will be publicly available for the entire scientific community.

Integrated multiple-microarray analysis and mendelian randomization to identify novel targets involved in diabetic nephropathy

Background

Diabetic nephropathy (DN), which is the main cause of renal failure in end-stage renal disease, is becoming a common chronic renal disease worldwide. Mendelian randomization (MR) is a genetic tool that is widely used to minimize confounding and reverse causation when identifying the causal effects of complex traits. In this study, we conducted an integrated multiple microarray analysis and large-scale plasma proteome MR analysis to identify candidate biomarkers and evaluate the causal effects of prospective therapeutic targets in DN.

Methods

Five DN gene expression datasets were selected from the Gene Expression Omnibus. The robust rank aggregation (RRA) method was used to integrate differentially expressed genes (DEGs) of glomerular samples between patients with DN and controls, followed by functional enrichment analysis. Protein quantitative trait loci were incorporated from seven different proteomic genome-wide association studies, and genetic association data on DN were obtained from FinnGen (3676 cases and 283,456 controls) for two-sample MR analysis. External validation and clinical correlation were also conducted.

Results

A total of 82 DEGs (53 upregulated and 29 downregulated) were identified through RRA integrated analysis. The enriched Gene Ontology annotations and Kyoto Encyclopedia of Genes and Genomes pathways of the DEGs were significantly enriched in neutrophil degranulation, neutrophil activation, proteoglycan binding, collagen binding, secretory granule lumen, gluconeogenesis, tricarboxylic acid cycle, and pentose phosphate pathways. MR analysis revealed that the genetically predicted levels of MHC class I polypeptide-related sequence B (MICB), granzyme A (GZMA), cathepsin S (CTSS), chloride intracellular channel protein 5, and ficolin-1 (FCN1) were causally associated with DN risk. Expression validation and clinical correlation analysis showed that MICB, GZMA, FCN1, and insulin-like growth factor 1 may participate in the development of DN, and carbonic anhydrase 2 and lipoprotein lipase may play protective roles in patients with DN.

Conclusion

Our integrated analysis identified novel biomarkers, including MICB and GZMA, which may help further understand the complicated mechanisms of DN and identify new target pathways for intervention.

Renoprotective effects of genetically proxied fibroblast growth factor 21: Mendelian randomization, proteome-wide and metabolome-wide association study

BACKGROUND

Fibroblast growth factor 21 (FGF21) has demonstrated efficacy for reducing liver fat and reversing non-alcoholic steatohepatitis in phase 2 clinical trials. It is also postulated to have anti-fibrotic effects and therefore may be amenable to repurposing for the prevention and treatment of chronic kidney disease (CKD).

METHODS

We leverage a missense genetic variant, rs739320 in the FGF21 gene, that associates with magnetic resonance imaging-derived liver fat as a clinically validated and biologically plausible instrumental variable for studying the effects of FGF21 analogs. Performing Mendelian randomization, we ascertain associations between instrumented FGF21 and kidney phenotypes, cardiometabolic disease risk factors, as well as the circulating proteome (Somalogic, 4907 aptamers) and metabolome (Nightingale platform, 249 metabolites).

RESULTS

We report consistent renoprotective associations of genetically proxied FGF21 effect, including higher glomerular filtration rates (p = 1.9 × 10^-4), higher urinary sodium excretion (p = 5.1 × 10^-11), and lower urine albumin-creatinine ratio (p = 3.6 × 10^-5). These favorable effects translated to lower CKD risk (odds ratio per rs739320 C-allele, 0.96; 95%CI, 0.94-0.98; p = 3.2 × 10^-4). Genetically proxied FGF21 effect was also associated with lower fasting insulin, waist-to-hip ratio, blood pressure (systolic and diastolic BP, p < 1.0 × 10^-07) and blood lipid (low-density lipoprotein cholesterol, triglycerides and apolipoprotein B, p < 6.5 × 10^-24) profiles. The latter associations are replicated in our metabolome-wide association study. Proteomic perturbations associated with genetically predicted FGF21 effect were consistent with fibrosis reduction.

CONCLUSION

This study highlights the pleiotropic effects of genetically proxied FGF21 and supports a re-purposing opportunity for the treatment and prevention of kidney disease specifically. Further work is required to triangulate these findings, towards possible clinical development of FGF21 towards the treatment and prevention of kidney disease.

Elevated plasma complement factor H related 5 protein is associated with venous thromboembolism

Venous thromboembolism (VTE) is a common, multi-causal disease with potentially serious short- and long-term complications. In clinical practice, there is a need for improved plasma biomarker-based tools for VTE diagnosis and risk prediction. Here we show, using proteomics profiling to screen plasma from patients with suspected acute VTE, and several case-control studies for VTE, how Complement Factor H Related 5 protein (CFHR5), a regulator of the alternative pathway of complement activation, is a VTE-associated plasma biomarker. In plasma, higher CFHR5 levels are associated with increased thrombin generation potential and recombinant CFHR5 enhanced platelet activation in vitro. GWAS analysis of ~52,000 participants identifies six loci associated with CFHR5 plasma levels, but Mendelian randomization do not demonstrate causality between CFHR5 and VTE. Our results indicate an important role for the regulation of the alternative pathway of complement activation in VTE and that CFHR5 represents a potential diagnostic and/or risk predictive plasma biomarker.

The Human Gastric Juice: A Promising Source for Gastric Cancer Biomarkers

Gastric cancer (GC) is a major public health problem worldwide, with high mortality rates due to late diagnosis and limited treatment options. Biomarker research is essential to improve the early detection of GC. Technological advances and research methodologies have improved diagnostic tools, identifying several potential biomarkers for GC, including microRNA, DNA methylation markers, and protein-based biomarkers. Although most studies have focused on identifying biomarkers in biofluids, the low specificity of these markers has limited their use in clinical practice. This is because many cancers share similar alterations and biomarkers, so obtaining them from the site of disease origin could yield more specific results. As a result, recent research efforts have shifted towards exploring gastric juice (GJ) as an alternative source for biomarker identification. Since GJ is a waste product during a gastroscopic examination, it could provide a "liquid biopsy" enriched with disease-specific biomarkers generated directly at the damaged site. Furthermore, as it contains secretions from the stomach lining, it could reflect changes associated with the developmental stage of GC. This narrative review describes some potential biomarkers for gastric cancer screening identified in gastric juice.

Association between circulating inflammatory markers and adult cancer risk: a Mendelian randomization analysis

Background

Tumour-promoting inflammation is a "hallmark" of cancer and conventional epidemiological studies have reported links between various inflammatory markers and cancer risk. The causal nature of these relationships and, thus, the suitability of these markers as intervention targets for cancer prevention is unclear.

Methods

We meta-analysed 6 genome-wide association studies of circulating inflammatory markers comprising 59,969 participants of European ancestry. We then used combined cis-Mendelian randomization and colocalisation analysis to evaluate the causal role of 66 circulating inflammatory markers in risk of 30 adult cancers in 338,162 cancer cases and up to 824,556 controls. Genetic instruments for inflammatory markers were constructed using genome-wide significant (P < 5.0 x 10-8) cis-acting SNPs (i.e. in or ±250 kb from the gene encoding the relevant protein) in weak linkage disequilibrium (LD, r2 < 0.10). Effect estimates were generated using inverse-variance weighted random-effects models and standard errors were inflated to account for weak LD between variants with reference to the 1000 Genomes Phase 3 CEU panel. A false discovery rate (FDR)-corrected P-value ("q-value") < 0.05 was used as a threshold to define "strong evidence" to support associations and 0.05 ≤ q-value < 0.20 to define "suggestive evidence". A colocalisation posterior probability (PPH4) > 70% was employed to indicate support for shared causal variants across inflammatory markers and cancer outcomes.

Results

We found strong evidence to support an association of genetically-proxied circulating pro-adrenomedullin concentrations with increased breast cancer risk (OR 1.19, 95% CI 1.10-1.29, q-value=0.033, PPH4=84.3%) and suggestive evidence to support associations of interleukin-23 receptor concentrations with increased pancreatic cancer risk (OR 1.42, 95% CI 1.20-1.69, q-value=0.055, PPH4=73.9%), prothrombin concentrations with decreased basal cell carcinoma risk (OR 0.66, 95% CI 0.53-0.81, q-value=0.067, PPH4=81.8%), macrophage migration inhibitory factor concentrations with increased bladder cancer risk (OR 1.14, 95% CI 1.05-1.23, q-value=0.072, PPH4=76.1%), and interleukin-1 receptor-like 1 concentrations with decreased triple-negative breast cancer risk (OR 0.92, 95% CI 0.88-0.97, q-value=0.15), PPH4=85.6%). For 22 of 30 cancer outcomes examined, there was little evidence (q-value ≥ 0.20) that any of the 66 circulating inflammatory markers examined were associated with cancer risk.

Conclusion

Our comprehensive joint Mendelian randomization and colocalisation analysis of the role of circulating inflammatory markers in cancer risk identified potential roles for 5 circulating inflammatory markers in risk of 5 site-specific cancers. Contrary to reports from some prior conventional epidemiological studies, we found little evidence of association of circulating inflammatory markers with the majority of site-specific cancers evaluated.

Causal factors underlying diabetes risk informed by Mendelian randomisation analysis: evidence, opportunities and challenges

Diabetes and its complications cause a heavy disease burden globally. Identifying exposures, risk factors and molecular processes causally associated with the development of diabetes can provide important evidence bases for disease prevention and spur novel therapeutic strategies. Mendelian randomisation (MR), an epidemiological approach that uses genetic instruments to infer causal associations between an exposure and an outcome, can be leveraged to complement evidence from observational and clinical studies. This narrative review aims to summarise the evidence on potential causal risk factors for diabetes by integrating published MR studies on type 1 and 2 diabetes, and to reflect on future perspectives of MR studies on diabetes. Despite the genetic influence on type 1 diabetes, few MR studies have been conducted to identify causal exposures or molecular processes leading to increased disease risk. In type 2 diabetes, MR analyses support causal associations of somatic, mental and lifestyle factors with development of the disease. These studies have also identified biomarkers, some of them derived from the gut microbiota, and molecular processes leading to increased disease risk. These studies provide valuable data to better understand disease pathophysiology and explore potential therapeutic targets. Because genetic association studies have mostly been restricted to participants of European descent, multi-ancestry cohorts are needed to examine the role of different types of physical activity, dietary components, metabolites, protein biomarkers and gut microbiome in diabetes development.

Genetic evidence implicating natriuretic peptide receptor-3 in cardiovascular disease risk: a Mendelian randomization study

BACKGROUND

C-type natriuretic peptide (CNP) is a known target for promoting growth and has been implicated as a therapeutic opportunity for the prevention and treatment of cardiovascular disease (CVD). This study aimed to explore the effect of CNP on CVD risk using the Mendelian randomization (MR) framework.

METHODS

Instrumental variables mimicking the effects of pharmacological intervention on CNP were identified as uncorrelated genetic variants located in the genes coding for its primary receptors, natriuretic peptide receptors-2 and 3 (NPR2 and NPR3), that associated with height. We performed MR and colocalization analyses to investigate the effects of NPR2 signalling and NPR3 function on CVD outcomes and risk factors. MR estimates were compared to those obtained when considering height variants from throughout the genome.

RESULTS

Genetically-proxied reduced NPR3 function was associated with a lower risk of CVD, with odds ratio (OR) 0.74 per standard deviation (SD) higher NPR3-predicted height, and 95% confidence interval (95% CI) 0.64-0.86. This effect was greater in magnitude than observed when considering height variants from throughout the genome. For CVD subtypes, similar MR associations for NPR3-predicted height were observed when considering the outcomes of coronary artery disease (0.75, 95% CI 0.60-0.92), stroke (0.69, 95% CI 0.50-0.95) and heart failure (0.77, 95% CI 0.58-1.02). Consideration of CVD risk factors identified systolic blood pressure (SBP) as a potential mediator of the NPR3-related CVD risk lowering. For stroke, we found that the MR estimate for NPR3 was greater in magnitude than could be explained by a genetically predicted SBP effect alone. Colocalization results largely supported the MR findings, with no evidence of results being driven by effects due to variants in linkage disequilibrium. There was no MR evidence supporting effects of NPR2 on CVD risk, although this null finding could be attributable to fewer genetic variants being identified to instrument this target.

CONCLUSIONS

This genetic analysis supports the cardioprotective effects of pharmacologically inhibiting NPR3 receptor function, which is only partly mediated by an effect on blood pressure. There was unlikely sufficient statistical power to investigate the cardioprotective effects of NPR2 signalling.