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.

Immuno-proteomic profiling reveals aberrant immune cell regulation in the airways of individuals with ongoing post-COVID-19 respiratory disease.

Some patients hospitalized with acute COVID-19 suffer respiratory symptoms that persist for many months. We delineated the immune-proteomic landscape in the airways and peripheral blood of healthy controls and post-COVID-19 patients 3 to 6 months after hospital discharge. Post-COVID-19 patients showed abnormal airway (but not plasma) proteomes, with an elevated concentration of proteins associated with apoptosis, tissue repair, and epithelial injury versus healthy individuals. Increased numbers of cytotoxic lymphocytes were observed in individuals with greater airway dysfunction, while increased B cell numbers and altered monocyte subsets were associated with more widespread lung abnormalities. A one-year follow-up of some post-COVID-19 patients indicated that these abnormalities resolved over time. In summary, COVID-19 causes a prolonged change to the airway immune landscape in those with persistent lung disease, with evidence of cell death and tissue repair linked to the ongoing activation of cytotoxic T cells.

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.

Genomic atlas of the human plasma proteome.

Although plasma proteins have important roles in biological processes and are the direct targets of many drugs, the genetic factors that control inter-individual variation in plasma protein levels are not well understood. Here we characterize the genetic architecture of the human plasma proteome in healthy blood donors from the INTERVAL study. We identify 1,927 genetic associations with 1,478 proteins, a fourfold increase on existing knowledge, including trans associations for 1,104 proteins. To understand the consequences of perturbations in plasma protein levels, we apply an integrated approach that links genetic variation with biological pathway, disease, and drug databases. We show that protein quantitative trait loci overlap with gene expression quantitative trait loci, as well as with disease-associated loci, and find evidence that protein biomarkers have causal roles in disease using Mendelian randomization analysis. By linking genetic factors to diseases via specific proteins, our analyses highlight potential therapeutic targets, opportunities for matching existing drugs with new disease indications, and potential safety concerns for drugs under development.

Genetic Landscape of the ACE2 Coronavirus Receptor.

BACKGROUND: SARS-CoV-2, the causal agent of COVID-19, enters human cells using the ACE2 (angiotensin-converting enzyme 2) protein as a receptor. ACE2 is thus key to the infection and treatment of the coronavirus. ACE2 is highly expressed in the heart and respiratory and gastrointestinal tracts, playing important regulatory roles in the cardiovascular and other biological systems. However, the genetic basis of the ACE2 protein levels is not well understood. METHODS: We have conducted the largest genome-wide association meta-analysis of plasma ACE2 levels in >28 000 individuals of the SCALLOP Consortium (Systematic and Combined Analysis of Olink Proteins). We summarize the cross-sectional epidemiological correlates of circulating ACE2. Using the summary statistics-based high-definition likelihood method, we estimate relevant genetic correlations with cardiometabolic phenotypes, COVID-19, and other human complex traits and diseases. We perform causal inference of soluble ACE2 on vascular disease outcomes and COVID-19 severity using mendelian randomization. We also perform in silico functional analysis by integrating with other types of omics data. RESULTS: We identified 10 loci, including 8 novel, capturing 30% of the heritability of the protein. We detected that plasma ACE2 was genetically correlated with vascular diseases, severe COVID-19, and a wide range of human complex diseases and medications. An X-chromosome cis-protein quantitative trait loci-based mendelian randomization analysis suggested a causal effect of elevated ACE2 levels on COVID-19 severity (odds ratio, 1.63 [95% CI, 1.10-2.42]; P=0.01), hospitalization (odds ratio, 1.52 [95% CI, 1.05-2.21]; P=0.03), and infection (odds ratio, 1.60 [95% CI, 1.08-2.37]; P=0.02). Tissue- and cell type-specific transcriptomic and epigenomic analysis revealed that the ACE2 regulatory variants were enriched for DNA methylation sites in blood immune cells. CONCLUSIONS: Human plasma ACE2 shares a genetic basis with cardiovascular disease, COVID-19, and other related diseases. The genetic architecture of the ACE2 protein is mapped, providing a useful resource for further biological and clinical studies on this coronavirus receptor.

Reduced monocyte and macrophage TNFSF15/TL1A expression is associated with susceptibility to inflammatory bowel disease.

Chronic inflammation in inflammatory bowel disease (IBD) results from a breakdown of intestinal immune homeostasis and compromise of the intestinal barrier. Genome-wide association studies have identified over 200 genetic loci associated with risk for IBD, but the functional mechanisms of most of these genetic variants remain unknown. Polymorphisms at the TNFSF15 locus, which encodes the TNF superfamily cytokine commonly known as TL1A, are associated with susceptibility to IBD in multiple ethnic groups. In a wide variety of murine models of inflammation including models of IBD, TNFSF15 promotes immunopathology by signaling through its receptor DR3. Such evidence has led to the hypothesis that expression of this lymphocyte costimulatory cytokine increases risk for IBD. In contrast, here we show that the IBD-risk haplotype at TNFSF15 is associated with decreased expression of the gene by peripheral blood monocytes in both healthy volunteers and IBD patients. This association persists under various stimulation conditions at both the RNA and protein levels and is maintained after macrophage differentiation. Utilizing a "recall-by-genotype" bioresource for allele-specific expression measurements in a functional fine-mapping assay, we localize the polymorphism controlling TNFSF15 expression to the regulatory region upstream of the gene. Through a T cell costimulation assay, we demonstrate that genetically regulated TNFSF15 has functional relevance. These findings indicate that genetically enhanced expression of TNFSF15 in specific cell types may confer protection against the development of IBD.

Insight into Genotype-Phenotype Associations through eQTL Mapping in Multiple Cell Types in Health and Immune-Mediated Disease.

Genome-wide association studies (GWAS) have transformed our understanding of the genetics of complex traits such as autoimmune diseases, but how risk variants contribute to pathogenesis remains largely unknown. Identifying genetic variants that affect gene expression (expression quantitative trait loci, or eQTLs) is crucial to addressing this. eQTLs vary between tissues and following in vitro cellular activation, but have not been examined in the context of human inflammatory diseases. We performed eQTL mapping in five primary immune cell types from patients with active inflammatory bowel disease (n = 91), anti-neutrophil cytoplasmic antibody-associated vasculitis (n = 46) and healthy controls (n = 43), revealing eQTLs present only in the context of active inflammatory disease. Moreover, we show that following treatment a proportion of these eQTLs disappear. Through joint analysis of expression data from multiple cell types, we reveal that previous estimates of eQTL immune cell-type specificity are likely to have been exaggerated. Finally, by analysing gene expression data from multiple cell types, we find eQTLs not previously identified by database mining at 34 inflammatory bowel disease-associated loci. In summary, this parallel eQTL analysis in multiple leucocyte subsets from patients with active disease provides new insights into the genetic basis of immune-mediated diseases.

Severe localised granulomatosis with polyangiitis (Wegener's granulomatosis) manifesting with extensive cranial nerve palsies and cranial diabetes insipidus: a case report and literature review.

BACKGROUND: Granulomatosis with polyangiitis (GPA, formerly Wegener's granulomatosis) is a multisystem vasculitis of small- to medium-sized blood vessels. Cranial involvement can result in cranial nerve palsies and, rarely, pituitary infiltration. CASE PRESENTATION: We describe the case of a 32 year-old woman with limited but severe GPA manifesting as progressive cranial nerve palsies and pituitary dysfunction. Our patient initially presented with localised ENT involvement, but despite treatment with methotrexate, she deteriorated. Granulomatous inflammatory tissue around the skull base resulted in cavernous sinus syndrome, facial nerve palsy, palsies of cranial nerves IX-XII (Collet-Sicard syndrome), and the rare complication of cranial diabetes insipidus due to pituitary infiltration. The glossopharyngeal, vagus and accessory nerve palsies resulted in severe dysphagia and she required nasogastric tube feeding. Her neurological deficits substantially improved with treatment including high dose corticosteroid, cyclophosphamide and rituximab. CONCLUSIONS: This case emphasises that serious morbidity can arise from localised cranial Wegener's granulomatosis in the absence of systemic disease. In such cases intensive induction immunosuppression is required. Analysis of previously reported cases of pituitary involvement in GPA reveals that this rare complication predominantly affects female patients.

MT-HESS: an efficient Bayesian approach for simultaneous association detection in OMICS datasets, with application to eQTL mapping in multiple tissues.

MOTIVATION: Analysing the joint association between a large set of responses and predictors is a fundamental statistical task in integrative genomics, exemplified by numerous expression Quantitative Trait Loci (eQTL) studies. Of particular interest are the so-called ': hotspots ': , important genetic variants that regulate the expression of many genes. Recently, attention has focussed on whether eQTLs are common to several tissues, cell-types or, more generally, conditions or whether they are specific to a particular condition. RESULTS: We have implemented MT-HESS, a Bayesian hierarchical model that analyses the association between a large set of predictors, e.g. SNPs, and many responses, e.g. gene expression, in multiple tissues, cells or conditions. Our Bayesian sparse regression algorithm goes beyond ': one-at-a-time ': association tests between SNPs and responses and uses a fully multivariate model search across all linear combinations of SNPs, coupled with a model of the correlation between condition/tissue-specific responses. In addition, we use a hierarchical structure to leverage shared information across different genes, thus improving the detection of hotspots. We show the increase of power resulting from our new approach in an extensive simulation study. Our analysis of two case studies highlights new hotspots that would remain undetected by standard approaches and shows how greater prediction power can be achieved when several tissues are jointly considered. AVAILABILITY AND IMPLEMENTATION: C[Formula: see text] source code and documentation including compilation instructions are available under GNU licence at http://www.mrc-bsu.cam.ac.uk/software/.

Comparison of gene expression microarray data with count-based RNA measurements informs microarray interpretation.

BACKGROUND: Although numerous investigations have compared gene expression microarray platforms, preprocessing methods and batch correction algorithms using constructed spike-in or dilution datasets, there remains a paucity of studies examining the properties of microarray data using diverse biological samples. Most microarray experiments seek to identify subtle differences between samples with variable background noise, a scenario poorly represented by constructed datasets. Thus, microarray users lack important information regarding the complexities introduced in real-world experimental settings. The recent development of a multiplexed, digital technology for nucleic acid measurement enables counting of individual RNA molecules without amplification and, for the first time, permits such a study. RESULTS: Using a set of human leukocyte subset RNA samples, we compared previously acquired microarray expression values with RNA molecule counts determined by the nCounter Analysis System (NanoString Technologies) in selected genes. We found that gene measurements across samples correlated well between the two platforms, particularly for high-variance genes, while genes deemed unexpressed by the nCounter generally had both low expression and low variance on the microarray. Confirming previous findings from spike-in and dilution datasets, this "gold-standard" comparison demonstrated signal compression that varied dramatically by expression level and, to a lesser extent, by dataset. Most importantly, examination of three different cell types revealed that noise levels differed across tissues. CONCLUSIONS: Microarray measurements generally correlate with relative RNA molecule counts within optimal ranges but suffer from expression-dependent accuracy bias and precision that varies across datasets. We urge microarray users to consider expression-level effects in signal interpretation and to evaluate noise properties in each dataset independently.

Serial analysis of clinical and imaging indices reveals prolonged efficacy of TNF-α and IL-6 receptor targeted therapies in refractory Takayasu arteritis.

OBJECTIVES: We analysed a large cohort of patients with Takayasu arteritis, seeking robust clinical evidence for prolonged responses to tumour necrosis factor-α (TNF-α) and interleukin-6 receptor (IL-6R) antagonists in severe refractory disease. METHODS: Case notes from ninety-eight patients with Takayasu arteritis were retrospectively reviewed. Drug treatment, laboratory and serial non-invasive imaging data were analysed, and the Indian Takayasu arteritis activity (ITAS) and damage scores (TADs) calculated. RESULTS: Nine patients were treated with biologic therapies. All had previously received high dose prednisolone and ≥1 conventional immunosuppressant. Five patients had failed cyclophosphamide. The patients prescribed biologics had more extensive arterial injury than the remainder of the cohort and persistent active disease (ITAS range 2-9, CRP 12-206 mg/L, TADs 3--1). Eight patients were prescribed anti-TNF-α therapy, three IL-6R blockade. The mean duration of anti-TNF-α treatment was 42 months (maximum 8 years). One patient developed new arterial stenoses while receiving anti-TNF-α and subsequently achieved disease remission with tocilizumab. Two patients have now demonstrated sustained responses to IL-6R inhibition at 19 and 20 months. Following introduction of biologic therapy, serial non-invasive imaging has revealed no significant progression in arterial injury. A significant fall in CRP (p<0.01), prednisolone dose (p<0.01) and ITAS (p<0.01) was observed, with no increase in TADs. CONCLUSIONS: We report for the first time sustained responses to both anti-TNF-α and IL6R antagonists in refractory Takayasu arteritis. As 5/9 patients were cyclophosphamide non-responders, we propose that biologics should now be considered ahead of cyclophosphamide in these young patients.

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 ).

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.

Genetic mechanisms of 184 neuro-related proteins in human plasma.

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

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

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

Somatostatin Receptor PET/MR Imaging of Inflammation in Patients With Large Vessel Vasculitis and Atherosclerosis.

BACKGROUND: Assessing inflammatory disease activity in large vessel vasculitis (LVV) can be challenging by conventional measures. OBJECTIVES: We aimed to investigate somatostatin receptor 2 (SST2) as a novel inflammation-specific molecular imaging target in LVV. METHODS: In a prospective, observational cohort study, in vivo arterial SST2 expression was assessed by positron emission tomography/magnetic resonance imaging (PET/MRI) using 68Ga-DOTATATE and 18F-FET-ßAG-TOCA. Ex vivo mapping of the imaging target was performed using immunofluorescence microscopy; imaging mass cytometry; and bulk, single-cell, and single-nucleus RNA sequencing. RESULTS: Sixty-one participants (LVV: n = 27; recent atherosclerotic myocardial infarction of ≤2 weeks: n = 25; control subjects with an oncologic indication for imaging: n = 9) were included. Index vessel SST2 maximum tissue-to-blood ratio was 61.8% (P < 0.0001) higher in active/grumbling LVV than inactive LVV and 34.6% (P = 0.0002) higher than myocardial infarction, with good diagnostic accuracy (area under the curve: ≥0.86; P < 0.001 for both). Arterial SST2 signal was not elevated in any of the control subjects. SST2 PET/MRI was generally consistent with 18F-fluorodeoxyglucose PET/computed tomography imaging in LVV patients with contemporaneous clinical scans but with very low background signal in the brain and heart, allowing for unimpeded assessment of nearby coronary, myocardial, and intracranial artery involvement. Clinically effective treatment for LVV was associated with a 0.49 ± 0.24 (standard error of the mean [SEM]) (P = 0.04; 22.3%) reduction in the SST2 maximum tissue-to-blood ratio after 9.3 ± 3.2 months. SST2 expression was localized to macrophages, pericytes, and perivascular adipocytes in vasculitis specimens, with specific receptor binding confirmed by autoradiography. SSTR2-expressing macrophages coexpressed proinflammatory markers. CONCLUSIONS: SST2 PET/MRI holds major promise for diagnosis and therapeutic monitoring in LVV. (PET Imaging of Giant Cell and Takayasu Arteritis [PITA], NCT04071691; Residual Inflammation and Plaque Progression Long-Term Evaluation [RIPPLE], NCT04073810).

Multi-omics identify falling LRRC15 as a COVID-19 severity marker and persistent pro-thrombotic signals in convalescence.

Patients with end-stage kidney disease (ESKD) are at high risk of severe COVID-19. Here, we perform longitudinal blood sampling of ESKD haemodialysis patients with COVID-19, collecting samples pre-infection, serially during infection, and after clinical recovery. Using plasma proteomics, and RNA-sequencing and flow cytometry of immune cells, we identify transcriptomic and proteomic signatures of COVID-19 severity, and find distinct temporal molecular profiles in patients with severe disease. Supervised learning reveals that the plasma proteome is a superior indicator of clinical severity than the PBMC transcriptome. We show that a decreasing trajectory of plasma LRRC15, a proposed co-receptor for SARS-CoV-2, is associated with a more severe clinical course. We observe that two months after the acute infection, patients still display dysregulated gene expression related to vascular, platelet and coagulation pathways, including PF4 (platelet factor 4), which may explain the prolonged thrombotic risk following COVID-19.

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

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

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.