Proteome-wide Mendelian randomization identifies causal links between blood proteins and severe COVID-19

Alish B. Palmos, Vincent Millischer, David K. Menon, Timothy R. Nicholson, Leonie Taams, Benedict Michael, COVID Clinical Neuroscience Study Consortium, Christopher Hübel, Gerome Breen
doi: https://doi.org/10.1101/2021.03.09.21253206

Abstract

The COVID-19 pandemic death toll now surpasses two million individuals and there is a need for early identification of individuals at increased risk of mortality. Host genetic variation partially drives the immune and biochemical responses to COVID-19 that lead to risk of mortality. We identify and prioritise blood proteins and biomarkers that may indicate increased risk for severe COVID-19, via a proteome Mendelian randomization approach by collecting genome-wide association study (GWAS) summary statistics for >4,000 blood proteins. After multiple testing correction, troponin I3, cardiac type (TNNI3) had the strongest effect (odds ratio (O.R.) of 6.86 per standard deviation increase in protein level), with proteinase 3 (PRTN3) (O.R.=2.48), major histocompatibility complex, class II, DQ alpha 2 (HLA-DQA2) (O.R.=2.29), the C4A-C4B heterodimer (O.R.=1.76) and low-density lipoprotein receptor-related protein associated protein 1 (LRPAP1) (O.R.=1.73) also being associated with higher odds of severe COVID-19. Conversely, major histocompatibility complex class I polypeptide-related sequence A (MHC1A) (O.R.=0.6) and natural cytotoxicity triggering receptor 3 (NCR3) (O.R.=0.46) were associated with lower odds. These proteins are involved in heart muscle contraction, natural killer and antigen presenting cells, and the major histocompatibility complex. Based on these findings, it may be possible to better predict which patients may develop severe COVID-19 and to design better treatments targeting the implicated mechanisms.

Introduction

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was identified in late 2019 in Wuhan, China, as the cause of what is now termed the coronavirus disease 2019 (COVID-19) that rapidly evolved into a global pandemic (1,2). As of February 2021, more than 100 million cases have been confirmed worldwide with total deaths exceeding 2 million (3). COVID-19 pathology encompasses a wide spectrum of clinical manifestations from asymptomatic, over mild, moderate, to severe infections (around 15%) (1,2). Severe COVID-19 commonly requires hospitalization and intensive care with assisted respiratory support, and subsequent respiratory failure is the most common reason for COVID-19 associated mortality (4).

A dysregulated pro- and anti-inflammatory immunomodulatory response drives much of the severe pathophysiology of COVID-19 and comprises alveolar damage, lung inflammation and pathology of an acute respiratory distress syndrome (1,2),(5,6). While initial descriptions and a response to corticosteroids led to a label of a “cytokine storm” (7) in severe disease cases, recent reports suggest that the abnormalities in host response may be more complex and nuanced (8). Given that the innate immune response has an individual-level genetic basis, genetic variants carried by an individual could play an important role in the individual-level immune response and, therefore, may influence progression and severity of COVID-19. Most biological systems are based on the interaction of biological and environmental factors. The influence of environmental factors can confound measured associations between biological markers and diseases outcomes. Therefore, studying immunomodulatory blood proteins in COVID-19 patients is difficult due to potential confounding effects of factors, such as initial viral exposure/inoculum, smoking behaviour and high body mass index (BMI). These factors themselves are associated with increased pro-inflammatory cytokine levels and may represent independent risk factors for severe COVID-19 (9–11).

Numerous genome-wide association studies (GWASs) in healthy populations associated genetic variants with immunomodulatory blood proteins (12–14). In addition, the COVID-19 Host Genetics Initiative carried out GWASs on COVID-19 subtypes to understand the role of host genetic factors in susceptibility and severity of COVID-19 (15). These data represent a powerful source of information to find new biomarkers and therapeutic leads. The method of Mendelian randomization can be applied to investigate the relationship between immunomodulatory blood proteins and severe COVID-19 infection. This exploits the fact that alleles are randomly inherited from parent to offspring in a manner analogous to a randomised-controlled trial, and allows estimation of putative causal effects of an exposure on a disease while avoiding confounding environmental effects, thus overcoming some of the limitations of observational studies. Recent advancements in Mendelian randomization methods allow use of GWAS summary statistics data, to identify genetic proxies (i.e., instrumental variables) of modifiable risk factors and testing their association with disease outcomes (16). We, therefore, conducted Mendelian randomization analyses between the levels of a large number of blood proteins and severe COVID-19 to identify putative causal associations that could be prioritised in clinical studies.

Methods

Blood protein GWAS data

In total, we amassed 5,305 sets of GWAS summary statistics for blood biomarkers (12,13,84–89). A systematic search was performed based on the ontology lookup service (OLS; www.ebi.ac.uk/ols/index) using R and the packages ‘rols’ and ‘gwasrapid’ between July 7th and July 27th, 2020. OLS is a repository for biomedical ontologies, such as gene ontology (GO) or the experimental factor ontology (EFO), including a systematic description of many experimental variables. First, all subnodes of the EFO ‘protein measurements’ (EFO:0004747) were determined using an iterative process based on the package ‘rols’. Overall, 628 unique EFO IDs were determined. Subsequently, all genetic associations reported in the GWAS Catalog (90) (www.ebi.ac.uk/gwas/) were determined and linked to the corresponding study using the package ‘gwasrapid’. One hundred and seventy-eight unique GWAS catalog accession IDs with available summary statistics were curated manually before inclusion. In order to expand the dataset, studies published at a later date were included manually and the first and the last author of studies without publicly available summary statistics were contacted.

In total, we included eight publications for which summary data was readily available and processed those using standard GWAS summary statistics quality control metrics including removal of incomplete genetic variants, variants with information metrics of lower than 0.6 and allele frequencies more extreme than 0.005 or 0.995. Allele frequencies were taken from the 1,000 Genomes Project data where needed (91). See Supplementary Table 1 for a full list of the studies included in these analyses.

COVID-19 GWAS data

We downloaded the COVID-19 Host Genome Initiative GWAS meta-analysis (15) of “very severe respiratory confirmed covid vs not hospitalized covid” (release 4, October 2020; www.covid19hg.org/results/). Cases and controls were defined as SARS-CoV-2 infected individuals within a clinical window of 15-90 days from diagnosis; with cases defined as mortality or the need for respiratory support (intubation, continuous positive airway pressure or bilevel positive airway pressure) and controls defined as surviving individuals who did not require respiratory support at any point since diagnosis.(15) The European sample consisted of 269 cases and 688 controls.

This GWAS was selected due to its extreme phenotyping, focusing specifically on mortality and respiratory failure (cases) compared to discharge from hospital without respiratory complications (controls). Subsequent COVID-19 GWAS are based on different, much less severe outcomes and broader phenotypes such as hospitalisation versus non hospitalisation that do not directly address the encompassing genetic basis for COVID-19 mortality or respiratory failure.

Mendelian randomization

To examine the influence of blood proteins on the risk of developing severe COVID-19, we selected genetic variants, predominantly single nucleotide polymorphisms (SNPs), that were strongly associated with actual blood protein levels in 3,609 genome-wide analyses of single proteins using robust methodologies (see supplement for more details on how the proteins were measured). Using these genetic loci as proxies for protein levels we performed an analysis using Mendelian randomization, a method that enables tests of putative causal associations of these blood proteins with the development of severe COVID-19. We used the Generalised Summary-based Mendelian randomization (GSMR) method as the base method (92).

For all GWASs, SNPs were used as instrumental variables and were selected by applying a suggestive p-value threshold (p < 5 × 10−6), in order to identify enough SNPs in common between the exposure (i.e., blood marker) and outcome (i.e. severe COVID-19). Where possible (due to SNPs in common between the exposure and the outcome), bidirectional analyses were performed between 3,609 blood markers and severe COVID-19, resulting in a total of 7,406 tests. To account for multiple testing, we calculated false discovery rate (FDR) corrected Q values using the p.adjust function in R (pFDR = 0.1).(93)

Sensitivity analyses

To test for robustness, we performed sensitivity analyses on our significant results from the GSMR analyses using additional Mendelian randomization methods, including the weighted median (WM), inverse variance-weighted (IVW), robust adjusted profile score (MR-RAPS), and MR-Egger methods (94–96). In order to pass our sensitivity analyses, at least three of these four methods must agree with the GSMR results.

Clinical application

The protein with the highest effect size was investigated for its clinical application. For this we used the Open Targets Platform (20). Specifically we focus on gene ontology and drug targets.

Results

Out of the possible 7,406 associations tested, we tested 3,609 associations with severe COVID-19 as the exposure and 3,797 associations with severe COVID-19 as the outcome. No tests with severe COVID-19 as the exposure reached a level of statistical significance and are therefore not reported. See Table 1 for significantly associated blood markers alongside test statistics.

View this table:
Table 1. Mendelian randomization results

This table details significant, false discovery rate-(FDR)-corrected Mendelian randomization results, using the Generalised Summary Data-based Mendelian randomization (GSMR) method. Seven blood markers were associated with a significantly increased risk for severe COVID-19 while two blood markers were associated with a significantly decreased risk for severe COVID-19. CREB3L4 and DEFB119 did not survive our sensitivity analyses. The genetic instruments for C4a-C4b and HLA-DQA2 were in low to moderate linkage disequilibrium but this did not influence results (see Supplementary Table 4). The table presents the log odds statistics and corresponding standard error as well as odds ratios, 95% confidence intervals, and the FDR-adjusted Q values (pFDR = 0.1).

Elevated risk for severe COVID-19

After multiple testing correction (pFDR = 0.1), seven blood markers were associated with a statistically significantly elevated risk of severe COVID-19 (see Figure 1). One standard deviation (SD) increase in troponin I, cardiac muscle (TNNI3) was associated with 589% (95% CI: 5.91-7.80, q = 0.06) higher odds of severe COVID-19. One SD increase in myeloblastin (PRTN3) was associated with 148% (95% CI: 2.04-2.91, q = 0.04) higher odds of severe COVID-19. One SD increase in major histocompatibility complex (MHC), class II, DQ alpha 2 (HLA-DQA2) was associated with 129% (95% CI: 1.94-2.64, q = 0.01) higher odds of severe COVID-19. One SD increase in complement components 4A and 4B (C4a-C4b) was associated with 76% (95% CI: 1.51-2.01, q = 0.02) higher odds of severe COVID-19. One SD increase in low-density lipoprotein (LDL) receptor-related protein associated protein 1 (LRPAP1) was associated with 73% (95% CI: 1.49-1.97, q = 0.02) higher odds of severe COVID-19 infection.

Figure 1.
Figure 1. Blood markers putatively causally associated with severe COVID-19

Summary figure of the false discovery rate-corrected (pFDR = 0.1) Mendelian randomization results using the Generalised Summary data-based Mendelian randomization (GSMR) method. Odds ratios (ORs) of the blood markers associated with severe COVID-19 are displayed on the x-axis (with 95% confidence intervals). The blood markers are displayed on the y-axis. The dashed red line represents an odds ratio of one (i.e., no effect). Seven blood markers were associated with a significantly increased risk for severe COVID-19 and two blood markers were associated with a significantly decreased risk for severe COVID-19 (qFDR ≤ 0.1). Markers which additionally survived our sensitivity analyses are marked with an asterisk (*). C4A = complement C4A (Rodgers blood group), C4B = complement C4B (Chido blood group), CREB3L4 = cAMP responsive element binding protein 3 like 4, DEFB119 = beta-defensin 119, HLA-DQA2 = major histocompatibility complex (MHC), class II, DQ alpha 2, LRPAP1 = low-density lipoprotein (LDL) receptor-related protein associated protein 1, MICA = MHC class I polypeptide-related sequence A, NCR3 = natural cytotoxicity triggering receptor 3, PRTN3 = proteinase 3, TNNI3 = troponin I3, cardiac type.

Decreased risk for severe COVID-19

A one SD increase in MHC class I polypeptide-related sequence A (MICA) was associated with 40% (95% CI: 0.41-0.78, q = 0.001) lower odds and one SD increase in natural cytotoxicity triggering receptor 3 (NCR3) was associated with 54% (95% CI: 0.06-0.86, q = 0.09) lower odds of severe COVID-19 (see Figure 1).

Sensitivity analyses

To further increase confidence in the findings, we additionally filtered our results for those concurrent with results calculated by additional Mendelian randomization methods (WM, IVW and MR-RAPS). Our sensitivity analyses revealed that the TNNI3, PRTN3, HLA-DQA2, C4A-C4B, LRPAP1, MICA and NCR3 associations with severe COVID-19 were confirmed. The two additional markers, cAMP responsive element binding protein 3 like 4 (CREB3L4) and beta-defensin 119 (DEFB119), associated with an increased risk for severe COVID-19 with the GSMR method, did not survive our sensitivity analyses. See Figure 2 for Mendelian randomization scatter plots of the sensitivity analyses and Supplementary Table 2 for a full results table of the sensitivity analyses.

Figure 2.
Figure 2. Sensitivity analyses for blood markers putatively causally associated with severe COVID-19

Mendelian randomization scatter plots using the inverse variance weighted, robust adjusted profile score, MR-Egger and weighted median methods associating each blood marker with severe COVID-19. Single nucleotide polymorphism (SNP) effects of the blood markers (i.e., exposures) are displayed on the x-axis and SNP effects of severe COVID-19 (outcome) are displayed on the y-axis. The horizontal line of each point represents the standard error of the SNP-exposure association and the vertical line the standard error of the SNP-outcome association. The effect of (A) TNNI3; (B) PRTN3; (C) HLA-DQ2; (D) DEFB119; (E) CREB3L4; (F) C4a-C4b; (G) LRPAP1; (H) MICA; (I) NCR3. SNP = single nucleotide polymorphism; C4A = complement C4A (Rodgers blood group), C4B = complement C4B (Chido blood group), CREB3L4 = cAMP responsive element binding protein 3 like 4, DEFB119 = beta-defensin 119, HLA-DQA2 = major histocompatibility complex (MHC), class II, DQ alpha 2, LRPAP1 = low-density lipoprotein (LDL) receptor-related protein associated protein 1, MICA = MHC class I polypeptide-related sequence A, NCR3 = natural cytotoxicity triggering receptor 3, PRTN3 = proteinase 3, TNNI3 = troponin I3, cardiac type.

Several of the proteins (MHC1A, NC3R, C4a-C4b, and HLA-DQA2) are encoded by non-overlapping genes in separate locations within the MHC (major histocompatibility complex) region, located on chromosome 6 (17–19). We verified that the associations observed are independent and not driven by linkage disequilibrium between the genetic instrumental variables used (see Supplementary Table 3). HLA proteins are vital in evoking immune responses to antigenic stimuli thereby orchestrating cellular and humoral immunity. These molecules facilitate the recognition of self or altered-self vs. non-self proteins by the body’s immune system (17–19).

Clinical application

Due to its large effect, TNNI3 was further investigated for potential clinical application using the Open Targets Platform (20). We highlight the upstream effect that TNNI3 may be having on heart muscle contraction, which may be having a negative downstream effect on a wide range of biological processes in the nervous, respiratory and digestive systems (see Supplementary Figure 1 for a bubble plot of genetically associated processes with TNNI3). In addition, we identify levosimendan, a positive inotropic agent having ATP-dependent potassium-channel-opening and calcium-sensitizing effects by binding to a calcium sensing myocardial complex consisting of cardiac troponin C and troponin I(20),(21,22). This is a clinically approved drug used for heart failure, kidney failure and SARS (https://www.targetvalidation.org/summary?drug=CHEMBL2051955). We highlight the potential role levosimendan for severe COVID-19 (see Supplementary Figure 2 for a bubble plot of druggable targets associated with TNNI3 and levosimendan), but emphasise the need for further exploration of this relationship.

Discussion

Our proteome-wide analysis of the influence of the blood levels of 3,609 proteins on risk for severe COVID-19 using Mendelian randomization, yielded evidence consistent with higher blood levels of TNNI3, PRTN3, HLADQA2, C4a-C4b and LRPAP1 being putative causal risk factors for severe COVID-19, and higher blood levels of MHC1A and NCR3 being protective against severe COVID-19. All of these proteins have detectable blood plasma or serum levels; we summarised functions of these blood markers in Supplementary Table 4. It is important to note that we did not identify the typical canonical immune proteins (such as interleukin-6 or C-reactive protein)(23,24), suggesting that with a larger database we may be able to pinpoint proteins that are upstream of major inflammatory mechanisms.

TNNI3, which codes for cardiac troponin I, showed the largest effect on COVID-19 severity in our Mendelian randomization analysis (OR per SD increase in TNNI3 level = 6.86, 95% C.I. 5.91-7.80). Cardiac troponin I is highly specific for heart muscle and higher concentrations in the blood are associated with higher general in-hospital mortality for a number of illnesses (25). It serves as a sensor of intracellular calcium, regulating the interaction between thick and thin sarcomere filaments during the heart muscle contraction (26). We hypothesise that this protein may be relevant to COVID-19 as ∼40% of deaths in a cohort from Wuhan (27) were attributed to myocardial damage and heart failure, as well as increased mortality in individuals with both COVID-19 and acute cardiac injury (28). Between 12-28% of patients with COVID-19 present with elevated troponin (29–32) are more likely to require intensive care (29,33), have higher in-hospital mortality(28,30,32,34–36) and show myocardial injury after recovery (37). Different cardiovascular complications, including myocarditis,(38) microangiopathy, myocardial infarction, or a dysregulated host immune response have been proposed to explain the high mortality in COVID-19(39) and troponin may be a predictor for mortality independent of inflammation.(40) In addition, TNNI3 has been associated with angiotensin-converting enzyme 2 (ACE2) upregulation in obstructive hypertrophic cardiomyopathy(41) and is shown to be a related target of ACE2(20). Given the importance of ACE2 in COVID-19,(42) this further implicates TNNI3 as an important target of severe disease pathophysiology (see Supplementary Table 5 for TNNI3 associated proteins). Levosimendan has been shown to exert action on the myocardial troponin complex,(20) (21,22)which involves complex interactions between troponin I and troponin C (43,44). Levosimendan is a calcium sensitiser, which canonically also exerts action on troponin C that is implicated in several diseases, including hypertrophic cardiomyopathy (45). Both troponin C and I (as well troponin T) are part of the troponin complex, with troponin I interacting with troponin C to regulate calcium sensitivity and contractility.(45) The translation of these findings to physiology is not straightforward since the effects of these target TNNI3 polymorphisms on the contribution of troponin I to myocardial Ca sensing are not established. Consequently, this finding may indicate either impaired calcium sensing in the presence of these SNPs, in which case a calcium sensitiser (such as levosimendan) may be helpful. Alternatively, these SNPs may drive greater Ca myocardial Ca sensitivity, and the increased blood levels of troponin I and poor outcomes seen in their presence may be due to myocardial necrosis in the presence of excessive myocardial stress – in which case Ca channel blockers might be indicated. In any event, our findings make a strong case for further exploration of this physiology and, potentially, clinical trials of promising interventions.(46)

MHC1A, also known as MICA, encodes a 62 kDa cell surface glycoprotein which is expressed on endothelial, dendritic, epithelial cells and fibroblasts. It is a target for cellular and humoral immune responses and becomes upregulated in response to stimuli of cellular stress like DNA damage and viral infection.(47) MHC1A binds to the natural killer group 2D (NKG2D) transmembrane protein activating natural killer cells, γδ T cells and CD8+ αβ T cells implicated in the antiviral response (48). Severe COVID-19 patients show high counts of NKG2C cells and low counts of NKG2D cells (49,50). Our findings suggest that a genetic propensity for higher blood MHC1A may protect against severe COVID-19, potentially through activating cytolytic cells (OR = 0.60, 95% CI: 0.41-0.78).

The natural cytotoxicity triggering receptor 3 (NCR3, also known as NKp30) is expressed by NK cells as well as by some peripheral blood CD8+ and γδ T cells and mediates NK cytotoxicity by facilitating the crosstalk between NK and dendritic cells.(51) NCR3 becomes upregulated in response to interferon γ and bacterial infection, suggesting a role in inflammatory and infectious diseases.(52) A separate transcriptomic analysis showed that reduced NCR3 transcription is associated with more severe influenza(53) while NCR3 plays important roles in both immunity against viral infections and in the development of immune evasion by viruses.(54) Our findings suggest that a genetic predisposition for higher blood NCR3 may protect against severe COVID-19 (OR = 0.46, 95% CI: 0.06-0.86).(55–57) However, we note that this association is with the soluble form of NCR3 found in plasma, complicating the interpretation of the results.

A one SD increase in levels of C4a-C4b was associated with COVID-19 severity (OR = 1.76, 95% CI 1.51-2.01). The C4a-C4b heterodimer is part of the complement cascade, which plays an important role in the host defense response against viral infections but has also been robustly associated with schizophrenia(58,59), with the associated C4 structural variants leading to possible altered neural C4A expression, excessive or inappropriate synaptic pruning(60,61). However, while a diagnosis of schizophrenia has been reported to be associated with a >7-fold increased risk for infection and almost a 3-fold increase in risk of mortality for COVID-19(62,63), additional sensitivity analyses did not show any significant Mendelian randomization between the general genetics of schizophrenia(64) and severe COVID-19 (see Supplementary Table 6). Activation of the classical, lectin and alternative complement pathways is associated with severity of COVID-19,(65–67) and microvascular injury and thrombosis in fatal cases.(67,68) Small non-randomised studies(69,70) have suggest the C5 inhibitor, eculizumab is a potential therapeutic for COVID-19. Finally in the MHC region, we also found evidence that HLA-DQA2 is causally associated with severe COVID-19 (OR = 2.29, 95% CI 1.94-2.64). In contrast to other HLA class II molecules, HLA-DQA2 is less polymorphic with a more discrete expression profile.(71),(72).

Turning to non-MHC proteins, we find proteinase 3 (PRTN3; OR = 2.48, 95% CI: 2.04-2.91) is causally associated with severe COVID-19. PRTN3 is a serine proteinase expressed in neutrophil granulocytes and a constituent of azurophil granules, that contain proteases and antibacterial proteins. PRTN3 is involved in neutrophil degranulation, highly overexpressed (29-fold) in naso-oropharyngeal samples of COVID-19 patients,(73) and associated with greater COVID-19 severity.(74) Neutrophil degranulation itself is dysregulated in COVID-19(75) while the myelopoiesis produces immature and dysfunctional neutrophils in severe COVID-19.(76) Given that PRTN3 is also implicated in immune dysregulation, sepsis and acute kidney injury,(77) this molecule may represent an interesting target to modulate immune responses against SARS-CoV-2.

The LDL receptor-related protein associated protein 1 (LRPAP1; OR = 1.73, 95% CI: 1.49-1.97) is involved as an escort protein in trafficking members of the LDL receptor family through the endoplasmic reticulum and Golgi apparatus preventing premature binding of ligands with these receptors.(78) In COVID-19, LDL cholesterol below or equal to 69 mg/dl has been associated with poor clinical outcomes.(79) Regarding a potential involvement in the development of neurological and cardiovascular symptoms, genomic variants in LRPAP1 have been associated with late-onset Alzheimer’s and Parkinson’s disease while serum anti-LRPAP1 is commonly used as a biomarker for atherosclerotic diseases.(80–82)

Our study has several limitations. Firstly, although we required confirmation of our findings by several Mendelian randomization methods, the p-value threshold for selecting genetic variants as instruments for our analyses was set at suggestive significance (p < 5 × 10−6) to allow enough to be identified for each protein, meaning that some of the genetic variants may not have true associations with protein levels. Secondly, some blood marker GWASs were excluded from our analyses due to unavailability, therefore, we may have missed associations with these markers. Thirdly, the severe COVID-19 GWAS used is of small sample size, reducing power to detect associations. However, it was a GWAS of the extreme phenotype of COVID-19 mortality or ventilator support, which may compensate for this.(83) Finally, our findings are limited to ancestral European populations due to data availability. Future endeavours should include participants from diverse ancestry groups.

Our results highlight the utility of applying large scale Mendelian randomization analyses to identify blood markers that may be causal for severe COVID-19. We demonstrate that TNNI3, PRTN3, HLADQA2, C4a-C4b and LRPAP1 may increase risk for severe COVID-19 while MHC1A and NCR3 may protect against severe COVID-19. These blood markers may play a role in the severity of COVID-19, indicating possible avenues to develop prognostic biomarkers and therapeutics for COVID-19.

Data Availability

All data and software is publicly available via original sources.

Funding

COVID-19 Clinical Neuroscience Study funded by the Medical Research Council (Grant Code: MR/V03605X/1).

Author contributions

ABP, VM and CH collected, processed and analysed the data. ABP, VM, CH and GB wrote the manuscript. DM, LT and BM helped write the manuscript and added expert knowledge. This work was carried out as part of the COVID Clinical Neuroscience Study consortium.

Competing interests

GB has received consultancy fees from Compass Pathways Ltd and Otsuka Ltd. There are no other conflicts to disclose.

Data availability

All data and software is publicly available via original sources.

Acknowledgements

This study is supported by the National Institute for Health Research (NIHR) Mental Health Translational Research Collaboration (MH-TRC). The views expressed are those of the authors and not necessarily those of the NIHR or the Department of Health and Social Care.

Footnotes

  • ↵† Indicates joint first authorship

  • ↵‡ Indicates joint last authorship

Abbreviations

C4A
complement C4A (Rodgers blood group)
C4B
complement C4B (Chido blood group)
CREB3L4
cAMP responsive element binding protein 3 like 4
COVID-19
coronavirus disease 2019
DEFB119
beta-defensin 119
GWAS
genome-wide association study
HLA-DQA2
major histocompatibility complex (MHC), class II, DQ alpha 2
LRPAP1
low-density lipoprotein (LDL) receptor-related protein associated protein 1
MICA
MHC class I polypeptide-related sequence A
NCR3
natural cytotoxicity triggering receptor 3
OR
odds ratio
PRTN3
proteinase 3
TNNI3
troponin I3, cardiac type
SARS-CoV-2
severe acute respiratory syndrome coronavirus 2
SE
standard error
SNP
single nucleotide polymorphism

References

  1. 1.↵
    Blanco-Melo, D., Nilsson-Payant, B. E., Liu, W.-C., et al. (2020) Imbalanced Host Response to SARS-CoV-2 Drives Development of COVID-19. Cell, 181, 1036–1045.e9.
    OpenUrlCrossRefPubMed
  2. 2.↵
    García, L. F. (2020) Immune Response, Inflammation, and the Clinical Spectrum of COVID-19. Front. Immunol., 11, 1441.
    OpenUrlCrossRefPubMed
  3. 3.↵
    Dong, E., Du, H. and Gardner, L. (2020) An interactive web-based dashboard to track COVID-19 in real time. Lancet Infect. Dis., 20, 533–534.
    OpenUrlCrossRefPubMed
  4. 4.↵
    Berlin, D. A., Gulick, R. M. and Martinez, F. J. (2020) Severe Covid-19. N. Engl. J. Med., 383, 2451–2460.
  5. 5.↵
    Kadkhoda, K. (2020) COVID-19: an Immunopathological View. mSphere, 5.
  6. 6.↵
    Ovsyannikova, I. G., Haralambieva, I. H., Crooke, S. N., et al. (2020) The role of host genetics in the immune response to SARS-CoV-2 and COVID-19 susceptibility and severity. Immunol. Rev., 296, 205–219.
    OpenUrlCrossRef
  7. 7.↵
    Fajgenbaum, D. C. and June, C. H. (2020) Cytokine Storm. N. Engl. J. Med., 383, 2255–2273.
    OpenUrlCrossRefPubMed
  8. 8.↵
    Sarma, A., Christenson, S., Mick, E., et al. (2021) COVID-19 ARDS is characterized by a dysregulated host response that differs from cytokine storm and is modified by dexamethasone. Res Sq.
  9. 9.↵
    Palmos, A. B., Chung, R., Frissa, S., et al. (2020) Reconsidering the reasons for heightened inflammation in major depressive disorder. J. Affect. Disord.
  10. 10.
    Palmos, A. B., Watson, S., Hughes, T., et al. (2019) Associations between childhood maltreatment and inflammatory markers. BJPsych Open, 5.
  11. 11.↵
    Monteiro, A. C., Suri, R., Emeruwa, I. O., et al. (2020) Obesity and smoking as risk factors for invasive mechanical ventilation in COVID-19: A retrospective, observational cohort study. PLoS One, 15, e0238552.
    OpenUrl
  12. 12.↵
    Sliz, E., Kalaoja, M., Ahola-Olli, A., et al. (2019) Genome-wide association study identifies seven novel loci associating with circulating cytokines and cell adhesion molecules in Finns. J. Med. Genet., 56, 607–616.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    Ahola-Olli, A. V., Würtz, P., Havulinna, A. S., et al. (2017) Genome-wide Association Study Identifies 27 Loci Influencing Concentrations of Circulating Cytokines and Growth Factors. Am. J. Hum. Genet., 100, 40–50.
    OpenUrl
  14. 14.↵
    Höglund, J., Rafati, N., Rask-Andersen, M., et al. (2019) Improved power and precision with whole genome sequencing data in genome-wide association studies of inflammatory biomarkers. Sci. Rep., 9, 16844.
    OpenUrlCrossRefPubMed
  15. 15.↵
    COVID-19 Host Genetics Initiative (2020) The COVID-19 Host Genetics Initiative, a global initiative to elucidate the role of host genetic factors in susceptibility and severity of the SARS-CoV-2 virus pandemic. Eur. J. Hum. Genet., 28, 715–718.
    OpenUrlPubMed
  16. 16.↵
    Burgess, S., Foley, C. N., Allara, E., et al. (2020) A robust and efficient method for Mendelian randomization with hundreds of genetic variants. Nat. Commun., 11, 376.
    OpenUrlPubMed
  17. 17.↵
    Sumitran-Holgersson, S. (2008) Beyond ABO and human histocompatibility antigen: other histocompatibility antigens with a role in transplantation. Curr. Opin. Organ Transplant., 13, 425–429.
    OpenUrlPubMed
  18. 18.
    Howell, W. M., Carter, V. and Clark, B. (2010) The HLA system: immunobiology, HLA typing, antibody screening and crossmatching techniques. J. Clin. Pathol., 63, 387–390.
    OpenUrlAbstract/FREE Full Text
  19. 19.↵
    Klein, J. and Sato, A. (2000) The HLA system. First of two parts. N. Engl. J. Med., 343, 702–709.
    OpenUrlCrossRefPubMedWeb of Science
  20. 20.↵
    Ghoussaini, M., Mountjoy, E., Carmona, M., et al. (2021) Open Targets Genetics: systematic identification of trait-associated genes using large-scale genetics and functional genomics. Nucleic Acids Res., 49, D1311–D1320.
    OpenUrl
  21. 21.↵
    Kachooei, E., Cordina, N. M., Potluri, P. R., et al. (2021) Phosphorylation of Troponin I finely controls the positioning of Troponin for the optimal regulation of cardiac muscle contraction. J. Mol. Cell. Cardiol., 150, 44–53.
    OpenUrl
  22. 22.↵
    Marston, S. and Zamora, J. E. (2020) Troponin structure and function: a view of recent progress. J. Muscle Res. Cell Motil., 41, 71–89.
    OpenUrl
  23. 23.↵
    McGowan, L. M., Smith, G. D., Gaunt, T. R., et al. Integrating Mendelian randomization and multiple-trait colocalization to uncover cell-specific inflammatory drivers of autoimmune and atopic disease. Integrating Mendelian randomization and multiple-trait colocalization to uncover cell-specific inflammatory drivers of autoimmune and atopic disease.
  24. 24.↵
    Georgakis, M. K., Gill, D., Rannikmäe, K., et al. (2019) Genetically Determined Levels of Circulating Cytokines and Risk of Stroke. Genetically Determined Levels of Circulating Cytokines and Risk of Stroke. Circulation (2019), 139, 256–268.
  25. 25.↵
    Reynolds, T., Cecconi, M., Collinson, P., et al. (2012) Raised serum cardiac troponin I concentrations predict hospital mortality in intensive care unit patients. Br. J. Anaesth., 109, 219–224.
    OpenUrlCrossRefPubMed
  26. 26.↵
    Pantou, M. P., Gourzi, P., Gkouziouta, A., et al. (2019) A case report of recessive restrictive cardiomyopathy caused by a novel mutation in cardiac troponin I (TNNI3). BMC Med. Genet., 20, 61.
    OpenUrl
  27. 27.↵
    Ruan, Q., Yang, K., Wang, W., et al. (2020) Clinical predictors of mortality due to COVID-19 based on an analysis of data of 150 patients from Wuhan, China. Intensive Care Med., 46, 846–848.
    OpenUrlCrossRefPubMed
  28. 28.↵
    Shi, S., Qin, M., Shen, B., et al. (2020) Association of Cardiac Injury With Mortality in Hospitalized Patients With COVID-19 in Wuhan, China. JAMA Cardiol, 5, 802–810.
    OpenUrl
  29. 29.↵
    Huang, C., Wang, Y., Li, X., et al. (2020) Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet, 395, 497–506.
    OpenUrlCrossRefPubMed
  30. 30.↵
    Zhou, F., Yu, T., Du, R., et al. (2020) Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. Lancet, 395, 1054–1062.
    OpenUrlCrossRefPubMed
  31. 31.
    Yang, X., Yu, Y., Xu, J., et al. (2020) Clinical course and outcomes of critically ill patients with SARS-CoV-2 pneumonia in Wuhan, China: a single-centered, retrospective, observational study. Lancet Respir Med, 8, 475–481.
    OpenUrl
  32. 32.↵
    Guo, T., Fan, Y., Chen, M., et al. (2020) Cardiovascular Implications of Fatal Outcomes of Patients With Coronavirus Disease 2019 (COVID-19). JAMA Cardiol, 5, 811–818.
    OpenUrl
  33. 33.↵
    Wang, D., Hu, B., Hu, C., et al. (2020) Clinical Characteristics of 138 Hospitalized Patients With 2019 Novel Coronavirus–Infected Pneumonia in Wuhan, China. JAMA, 323, 1061–1069.
    OpenUrlCrossRefPubMed
  34. 34.↵
    Chen, T., Wu, D., Chen, H., et al. (2020) Clinical characteristics of 113 deceased patients with coronavirus disease 2019: retrospective study. BMJ, 368, m1091.
    OpenUrlAbstract/FREE Full Text
  35. 35.
    Nie, S.-F., Yu, M., Xie, T., et al. (2020) Cardiac Troponin I Is an Independent Predictor for Mortality in Hospitalized Patients With COVID-19. Circulation, 142, 608–610.
    OpenUrl
  36. 36.↵
    Lombardi, C. M., Carubelli, V., Iorio, A., et al. (2020) Association of Troponin Levels With Mortality in Italian Patients Hospitalized With Coronavirus Disease 2019: Results of a Multicenter Study. JAMA Cardiol, 5, 1274–1280.
    OpenUrl
  37. 37.↵
    Knight, D. S., Kotecha, T., Razvi, Y., et al. (2020) COVID-19: Myocardial Injury in Survivors. Circulation, 142, 1120–1122.
    OpenUrlPubMed
  38. 38.↵
    Doyen, D., Moceri, P., Ducreux, D., et al. (2020) Myocarditis in a patient with COVID-19: a cause of raised troponin and ECG changes. Lancet, 395, 1516.
    OpenUrlPubMed
  39. 39.↵
    Tersalvi, G., Vicenzi, M., Calabretta, D., et al. (2020) Elevated Troponin in Patients With Coronavirus Disease 2019: Possible Mechanisms. J. Card. Fail., 26, 470–475.
    OpenUrlCrossRefPubMed
  40. 40.↵
    Majure, D. T., Gruberg, L., Saba, S. G., et al. (2021) Usefulness of Elevated Troponin to Predict Death in Patients With COVID-19 and Myocardial Injury. Am. J. Cardiol., 138, 100–106.
    OpenUrl
  41. 41.↵
    Bos, J. M., Hebl Virginia Bs, Oberg, A.L., et al. (2020) Marked Up-Regulation of ACE2 in Hearts of Patients With Obstructive Hypertrophic Cardiomyopathy: Implications for SARS-CoV-2eMediated COVID-19. Mayo Clinic Proceedings; Rochester, 95, 1354–1368.
    OpenUrl
  42. 42.↵
    Ni, W., Yang, X., Yang, D., et al. (2020) Role of angiotensin-converting enzyme 2 (ACE2) in COVID-19. Crit. Care, 24, 422.
    OpenUrl
  43. 43.↵
    Genchev, G. Z., Kobayashi, M., Kobayashi, T., et al. (2020) Molecular dynamics provides new insights into the mechanism of calcium signal transduction and inter-domain interactions in cardiac troponin. FEBS Open Bio.
  44. 44.↵
    Marston, S. B. and Redwood, C. S. (2003) Modulation of thin filament activation by breakdown or isoform switching of thin filament proteins: physiological and pathological implications. Circ. Res., 93, 1170–1178.
    OpenUrlAbstract/FREE Full Text
  45. 45.↵
    Robertson, I. M., Baryshnikova, O. K., Li, M. X., et al. (2008) Defining the binding site of levosimendan and its analogues in a regulatory cardiac troponin C-troponin I complex. Biochemistry, 47, 7485–7495.
    OpenUrlCrossRefPubMedWeb of Science
  46. 46.↵
    Akhmerov, A. and Marbán, E. (2020) COVID-19 and the Heart. Circ. Res., 126, 1443–1455.
    OpenUrlCrossRefPubMed
  47. 47.↵
    Ahmed, Z. and Askar, M. (2017) MICA (MHC class I polypeptide-related sequence A). Atlas Genet. Cytogenet. Oncol. Haematol.
  48. 48.↵
    Chen, J. L.-Y., Chang, C.-C., Huang, Y.-S., et al. (2018) Persistently elevated soluble MHC class I polypeptide-related sequence A and transforming growth factor-β1 levels are poor prognostic factors in head and neck squamous cell carcinoma after definitive chemoradiotherapy. PLoS One, 13, e0202224.
    OpenUrl
  49. 49.↵
    Maucourant, C., Filipovic, I., Ponzetta, A., et al. (2020) Natural killer cell immunotypes related to COVID-19 disease severity. Sci Immunol, 5.
  50. 50.↵
    Maucourant, C., Filipovic, I., Ponzetta, A., et al. (2020) Natural killer cell activation related to clinical outcome of COVID-19. medRxiv.
  51. 51.↵
    Mantovani, S., Oliviero, B., Lombardi, A., et al. (2019) Deficient natural killer cell NKp30-mediated function and altered NCR3 splice variants in hepatocellular carcinoma: Hepatology. Hepatology, 69, 1165–1179.
    OpenUrl
  52. 52.↵
    Mulcahy, H., O’Rourke, K. P., Adams, C., et al. (2006) LST1 and NCR3 expression in autoimmune inflammation and in response to IFN-gamma, LPS and microbial infection. Immunogenetics, 57, 893–903.
    OpenUrl
  53. 53.↵
    Hoang, L. T., Tolfvenstam, T., Ooi, E. E., et al. (2014) Patient-based transcriptome-wide analysis identify interferon and ubiquination pathways as potential predictors of influenza A disease severity. PLoS One, 9, e111640.
    OpenUrl
  54. 54.↵
    Hammer, Q., Rückert, T., Borst, E. M., et al. (2018) Peptide-specific recognition of human cytomegalovirus strains controls adaptive natural killer cells. Nat. Immunol., 19, 453–463.
    OpenUrlCrossRefPubMed
  55. 55.↵
    Jarahian, M., Fiedler, M., Cohnen, A., et al. (2011) Modulation of NKp30- and NKp46-Mediated Natural Killer Cell Responses by Poxviral Hemagglutinin. Modulation of NKp30- and NKp46-Mediated Natural Killer Cell Responses by Poxviral Hemagglutinin. PLoS Pathogens (2011), 7, e1002195.
    OpenUrlCrossRefPubMed
  56. 56.
    Semeraro, M., Rusakiewicz, S., Minard-Colin, V., et al. (2015) Clinical impact of the NKp30/B7-H6 axis in high-risk neuroblastoma patients. Sci. Transl. Med., 7, 283ra55.
    OpenUrlAbstract/FREE Full Text
  57. 57.↵
    Han, B., Mao, F.-Y., Zhao, Y.-L., et al. (2018) Altered NKp30, NKp46, NKG2D, and DNAM-1 Expression on Circulating NK Cells Is Associated with Tumor Progression in Human Gastric Cancer. J Immunol Res, 2018, 6248590.
    OpenUrl
  58. 58.↵
    Sekar, A., Bialas, A. R., de Rivera, H., et al. (2016) Schizophrenia risk from complex variation of complement component 4. Nature, 530, 177–183.
    OpenUrlCrossRefPubMed
  59. 59.↵
    Li, N., Zhang, J., Liao, D., et al. (2017) Corrigendum: Association between C4, C4A, and C4B copy number variations and susceptibility to autoimmune diseases: a meta-analysis. Sci. Rep., 7, 46785.
    OpenUrl
  60. 60.↵
    Woo, J. J., Pouget, J. G., Zai, C. C., et al. (2020) The complement system in schizophrenia: where are we now and what’s next? Mol. Psychiatry, 25, 114–130.
    OpenUrl
  61. 61.↵
    Nimgaonkar, V. L., Prasad, K. M., Chowdari, K. V., et al. (2017) The complement system: a gateway to gene-environment interactions in schizophrenia pathogenesis. Mol. Psychiatry, 22, 1554–1561.
    OpenUrlCrossRefPubMed
  62. 62.↵
    Wang, Q., Xu, R. and Volkow, N. D. (2021) Increased risk of COVID-19 infection and mortality in people with mental disorders: analysis from electronic health records in the United States. World Psychiatry, 20, 124–130.
    OpenUrl
  63. 63.↵
    Nemani, K., Li, C., Olfson, M., et al. (2021) Association of Psychiatric Disorders With Mortality Among Patients With COVID-19. JAMA Psychiatry.
  64. 64.↵
    Pardiñas, A. F., Holmans, P., Pocklington, A. J., et al. (2018) Common schizophrenia alleles are enriched in mutation-intolerant genes and in regions under strong background selection. Nat. Genet., 50, 381–389.
    OpenUrlCrossRefPubMed
  65. 65.↵
    Holter, J. C., Pischke, S. E., de Boer, E., et al. Systemic Complement Activation Is Associated with Respiratory Failure in COVID-19 Hospitalized Patients: A Prospective Cohort Study. Systemic Complement Activation Is Associated with Respiratory Failure in COVID-19 Hospitalized Patients: A Prospective Cohort Study. SSRN Electronic Journal.
  66. 66.
    Malaquias, M. A. S., Gadotti, A. C., Motta-Junior, J. da S., et al. (2020) The role of the lectin pathway of the complement system in SARS-CoV-2 lung injury. Transl. Res.
  67. 67.↵
    Macor, P., Durigutto, P., Mangogna, A., et al. (2021) Multi-organ complement deposition in COVID-19 patients. medRxiv.
  68. 68.↵
    Magro, C., Mulvey, J. J., Berlin, D., et al. (2020) Complement associated microvascular injury and thrombosis in the pathogenesis of severe COVID-19 infection: A report of five cases. Transl. Res., 220, 1–13.
    OpenUrlCrossRefPubMed
  69. 69.↵
    Annane, D., Heming, N., Grimaldi-Bensouda, L., et al. (2020) Eculizumab as an emergency treatment for adult patients with severe COVID-19 in the intensive care unit: A proof-of-concept study. EClinicalMedicine, 28, 100590.
    OpenUrl
  70. 70.↵
    Laurence, J., Mulvey, J. J., Seshadri, M., et al. (2020) Anti-complement C5 therapy with eculizumab in three cases of critical COVID-19. Clin. Immunol., 219, 108555.
    OpenUrl
  71. 71.↵
    Yu, L. P. and Sheehy, M. J. (1991) The cryptic HLA-DQA2 (“DX alpha”) gene is expressed in human B cell lines. J. Immunol., 147, 4393–4397.
    OpenUrlAbstract
  72. 72.↵
    Lenormand, C., Bausinger, H., Gross, F., et al. (2012) HLA-DQA2 and HLA-DQB2 genes are specifically expressed in human Langerhans cells and encode a new HLA class II molecule. J. Immunol., 188, 3903–3911.
    OpenUrlAbstract/FREE Full Text
  73. 73.↵
    Akgun, E., Tuzuner, M. B., Sahin, B., et al. (2020) Proteins associated with neutrophil degranulation are upregulated in nasopharyngeal swabs from SARS-CoV-2 patients. PLoS One, 15, e0240012.
    OpenUrl
  74. 74.↵
    Gisby, J., Clarke, C. L., Medjeral-Thomas, N., et al. (2020) Longitudinal proteomic profiling of high-risk patients with COVID-19 reveals markers of severity and predictors of fatal disease. Longitudinal proteomic profiling of high-risk patients with COVID-19 reveals markers of severity and predictors of fatal disease. bioRxiv (2020).
  75. 75.↵
    Overmyer, K. A., Shishkova, E., Miller, I. J., et al. (2021) Large-Scale Multi-omic Analysis of COVID-19 Severity. Cell Syst, 12, 23–40.e7.
    OpenUrl
  76. 76.↵
    Schulte-Schrepping, J., Reusch, N., Paclik, D., et al. (2020) Severe COVID-19 Is Marked by a Dysregulated Myeloid Cell Compartment. Cell, 182, 1419–1440.e23.
    OpenUrlPubMed
  77. 77.↵
    Wang, C., Wang, Q., Li, R., et al. (2019) LTF, PRTN3, and MNDA in Synovial Fluid as Promising Biomarkers for Periprosthetic Joint Infection: Identification by Quadrupole Orbital-Trap Mass Spectrometry. J. Bone Joint Surg. Am., 101, 2226–2234.
    OpenUrl
  78. 78.↵
    Van Leuven, F., Thiry, E., Stas, L., et al. (1998) Analysis of the human LRPAP1 gene coding for the lipoprotein receptor-associated protein: identification of 22 polymorphisms and one mutation. Genomics, 52, 145–151.
    OpenUrlPubMed
  79. 79.↵
    Aparisi, A., Iglesias-Echeverria, C., Ybarra-Falcon, C., et al. (2020) Low-density lipoprotein cholesterol levels are associated with poor clinical outcomes in COVID-19. medRxiv.
  80. 80.↵
    Sumazaki, M., Shimada, H., Ito, M., et al. (2020) Serum anti-LRPAP1 is a common biomarker for digestive organ cancers and atherosclerotic diseases. Cancer Sci., 111, 4453–4464.
    OpenUrl
  81. 81.
    Singh, N. K., Banerjee, B. D., Bala, K., et al. (2014) APOE and LRPAP1 gene polymorphism and risk of Parkinson’s disease. Neurol. Sci., 35, 1075–1081.
    OpenUrl
  82. 82.↵
    Sánchez, L., Alvarez, V. and González, P. (2001) Variation in the LRP-associated protein gene (LRPAP1) is associated with late-onset Alzheimer disease. American journal of.
  83. 83.↵
    Zabaneh, D., Krapohl, E., Gaspar, H. A., et al. (2018) A genome-wide association study for extremely high intelligence. Mol. Psychiatry, 23, 1226–1232.
    OpenUrl
  84. 84.↵
    Sun, B. B., Maranville, J. C., Peters, J. E., et al. (2018) Genomic atlas of the human plasma proteome. Nature, 558, 73–79.
    OpenUrlCrossRefPubMed
  85. 85.
    Suhre, K., Arnold, M., Bhagwat, A. M., et al. (2017) Connecting genetic risk to disease end points through the human blood plasma proteome. Nat. Commun., 8, 14357.
    OpenUrlCrossRefPubMed
  86. 86.
    Wood, A. R., Perry, J. R. B., Tanaka, T., et al. (2013) Imputation of variants from the 1000 Genomes Project modestly improves known associations and can identify low-frequency variant-phenotype associations undetected by HapMap based imputation. PLoS One, 8, e64343.
    OpenUrlCrossRefPubMed
  87. 87.
    Folkersen, L., Fauman, E., Sabater-Lleal, M., et al. (2017) Mapping of 79 loci for 83 plasma protein biomarkers in cardiovascular disease. PLoS Genet., 13, e1006706.
    OpenUrlCrossRefPubMed
  88. 88.
    Folkersen, L., Gustafsson, S., Wang, Q., et al. (2020) Genomic and drug target evaluation of 90 cardiovascular proteins in 30,931 individuals. Nat Metab, 2, 1135–1148.
    OpenUrl
  89. 89.↵
    Bretherick, A. D., Canela-Xandri, O., Joshi, P. K., et al. (2020) Linking protein to phenotype with Mendelian Randomization detects 38 proteins with causal roles in human diseases and traits. PLoS Genet., 16, e1008785.
    OpenUrl
  90. 90.↵
    MacArthur, J., Bowler, E., Cerezo, M., et al. (2017) The new NHGRI-EBI Catalog of published genome-wide association studies (GWAS Catalog). Nucleic Acids Res., 45, D896–D901.
    OpenUrlCrossRefPubMed
  91. 91.↵
    1000 Genomes Project Consortium, The 1000 Genomes Project, Auton, A., Brooks, L. D., et al. (2015) A global reference for human genetic variation. Nature, 526, 68–74.
    OpenUrlCrossRefPubMed
  92. 92.↵
    Zhu, Z., Zheng, Z., Zhang, F., et al. (2018) Causal associations between risk factors and common diseases inferred from GWAS summary data. Nat. Commun., 9, 224.
    OpenUrl
  93. 93.↵
    Benjamini, Y. and Hochberg, Y. (1995) Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Series B Stat. Methodol., 57, 289–300.
    OpenUrlCrossRefPubMed
  94. 94.↵
    Hemani, G., Bowden, J. and Davey Smith, G. (2018) Evaluating the potential role of pleiotropy in Mendelian randomization studies. Hum. Mol. Genet., 27, R195–R208.
    OpenUrlPubMed
  95. 95.
    Bowden, J., Del Greco MF., Minelli, C., et al. (2017) A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization. Stat. Med., 36, 1783–1802.
    OpenUrlCrossRefPubMed
  96. 96.↵
    Teumer, A. (2018) Common Methods for Performing Mendelian Randomization. Front Cardiovasc Med, 5, 51.
    OpenUrl
Back to top
Posted March 11, 2021.
Download PDF

Supplementary Material

Data/Code
Email
Proteome-wide Mendelian randomization identifies causal links between blood proteins and severe COVID-19
Alish B. Palmos, Vincent Millischer, David K. Menon, Timothy R. Nicholson, Leonie Taams, Benedict Michael, COVID Clinical Neuroscience Study Consortium, Christopher Hübel, Gerome Breen
medRxiv 2021.03.09.21253206; doi: https://doi.org/10.1101/2021.03.09.21253206
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area