Shared genetic etiology and causality between COVID-19 and venous thromboembolism: evidence from genome-wide cross trait analysis and bi-directional Mendelian randomization study
==================================================================================================================================================================================

* Xin Huang
* Minhao Yao
* Peixin Tian
* Jason Y.Y. Wong
* Zilin Li
* Zhonghua Liu
* Jie V. Zhao

## Abstract

Venous thromboembolism (VTE) occurs in up to one third patients with COVID-19. VTE and COVID-19 may share a common genetic architecture, which has not been clarified yet. To fill this gap, we leveraged summary-level genetic data from the latest COVID-19 host genetics consortium and UK Biobank and examined the shared genetic etiology and causal relationship between COVID-19 and VTE. The cross-trait analysis identified 8, 11, and 7 shared loci between VTE and severe COVID-19, COVID-19 hospitalization, SARS-CoV-2 infection respectively, in 13 genes involved in coagulation and immune function and enriched in the lung. Co-localization analysis identified eight shared loci in *ABO, ADAMTS13* and *FUT2* genes. Bi-direction Mendelian randomization suggested that VTE was associated with higher risks of all COVID-19 related traits, and SARS-CoV-2 infection was associated with higher risk of VTE. Our study provided timely evidence and novel insights into the genetic etiology between COVID-19 and VTE.

## Introduction

Coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, has led to a worldwide pandemic since March 2020, and caused repeated waves of outbreaks across the globe. Venous thromboembolism (VTE) is one common and serious comorbidity. Notably, thrombotic events occur in up to one-third of patients with COVID-19.1 People with more pronounced thrombotic symptoms are more likely to develop severe COVID-19,2 and thrombotic complications are a well-established predictor of mortality in people with COVID-19.1 Consistently, several genome-wide association studies (GWASs) of COVID-19 related traits (severe COVID-19, COVID-19 hospitalization and SARS-CoV-2 infection) have identified genetic loci in *ABO*,3-7 an established gene related to thrombosis. These evidence suggests that VTE and COVID-19 may share common genetic architecture. Identifying shared genetic factors contributing to both COVID-19 and VTE can provide novel insights into disease pathogenesis and pinpoint targets for therapeutic development or drug repurposing. However, to our knowledge, this problem has not been comprehensively investigated.

Meanwhile, the causal relationship between VTE and COVID-19 has not been clarified. On the one hand, VTE predicts COVID-19 severity and mortality.1,2 On the other hand, the cytokine storm and excessive inflammation caused by SARS-CoV-2 infection is hypothesized to lead to the systemic coagulation dysfunction;8-11 therefore, SARS-CoV-2 infection may increase the risk of VTE. Many institutional evidence-based guidelines supported the use of prophylactic treatments such as anticoagulants for thromboprophylaxis in COVID-19 patients.12 However, the recommendation is mainly based on observational studies, which might be subject to unmeasured confounding bias. To resolve this issue, Mendelian randomization (MR) study uses genetic variants as instrument for causal inference and can provide unbiased causal effect estimate even in the presence of unmeasured confounding.13 One previous MR study in a relatively smaller GWAS for COVID-19 suggested that genetically predicted VTE is associated with higher risk of COVID-19 hospitalization and SARS-CoV-2 infection,14 however, the association with severe COVID-19 is uncertain possibly due to the limited sample size. The reverse association, i.e., the association of COVID-19 with VTE has not been examined. The recently updated GWAS for COVID-19 with doubled sample size provided a well-powered study dataset to investigate the relationship between VTE and COVID-19.

Taken together, to fill this knowledge gap, we performed a genome-wide cross-trait analysis for VTE and three COVID-19 related traits to estimate their genetic correlations and shared genetic components using publicly available summary level data from GWASs of VTE and COVID-19. We also applied bi-directional Mendelian randomization methods to study the causal relationship between VTE and three COVID-19 related traits.

## Results

### Genetic correlation of VTE with COVID-19 related traits

The heritability (h2) estimated by Linkage disequilibrium score regression (LDSC) analysis suggested severe COVID-19, COVID-19 hospitalization, SARS-CoV-2 infection and VTE are heritable (P<0.05 shown in Supplementary Table 1). We found positive overall genetic correlation of VTE with COVID-19 hospitalization (rg = 0.2320, P-value= 0.0092) (shown in Figure 1 and Supplementary Table 2). The genetic correlation of VTE and severe COVID-19 and SARS-CoV-2 were positive but not significant at p value of 0.05 level. Further partitioned LDSC analysis found that severe COVID-19, COVID-19 hospitalization and SARS-CoV-2 infection are genetically correlated with VTE in 5 (Fetal DNase I hypersensitivity sites (DHS), H3K4me1, H3K9ac H3K27ac and transcription factor-binding site (TFBS)) (Figure 2), 11 (except for super enhancers regions) and 10 (except for conserved and super enhancers regions) of the 12 functional categories, respectively.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/06/01/2022.05.21.22275413/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2022/06/01/2022.05.21.22275413/F1)

Figure 1. Genetic correlation of VTE with COVID-19 related traits

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/06/01/2022.05.21.22275413/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2022/06/01/2022.05.21.22275413/F2)

Figure 2. Partitioned genetic correlation between VTE and COVID-19 related traits.
The x-axis represents the 12 functional categories and the y-axis represents the estimated partitioned genetic correlation. The significant functional categories (P<0.05) are starred. DGF: DNaseI digital genomic footprinting; DHS: DNase I hypersensitivity site; TFBS: transcription factor binding sites.

### Multi-Trait Analysis of GWAS (MTAG)

Based on the COVID-19 host genetics initiative (HGI) updated GWAS meta-analysis (release 7), there were 45 genome-wide significant (P<5×10−8) and uncorrelated (r2<0.01) loci for severe COVID-19, 46 for COVID-19 hospitalization and 26 for SARS-CoV-2 infection (Supplementary Table 3), and 12 for VTE in the GWAS from UK Biobank (Supplementary Table 4). Compared with these previously identified SNPs, using MTAG incorporating information from GWAS of COVID-19 and VTE we did not identify additional genome-wide significant loci for COVID-19 traits or VTE.

As shown in Table 1, we identified 8 shared novel genetic loci associated with both VTE and severe COVID-19, 11 with both VTE and COVID-19 hospitalization and 7 with both VTE and SARS-CoV-2 infection (Pmeta < 5×10−8; single trait P <0.05). Figure 3 displays the Manhattan plots of these results. In line with previous studies,3-6,15,16 we identified *ABO* gene and *FUT2* gene, which contributed to both VTE and COVID-19. Notably, we identified seven novel genes which have not been reported yet, including *LINC00970*, and six protein-coding genes (*TSPAN15, ADAMTS13, F5, DNAJB4, SLC39A8* and *OBSCN)*.

View this table:
[Table 1.](http://medrxiv.org/content/early/2022/06/01/2022.05.21.22275413/T1)

Table 1. Genome-wide significant loci associated with the COVID-19 traits and VTE in cross trait meta-analysis

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/06/01/2022.05.21.22275413/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2022/06/01/2022.05.21.22275413/F3)

Figure 3. Circular Manhattan plot of cross-trait analysis for (a) VTE and severe COVID-19, (b) VTE and COVID-19 hospitalization, (c) VTE and SARS-CoV-2 infection.
Each point represents a SNP, and shared significant loci with meta-analysis P-value < 5×10−8 and single-trait P-value < 0.05 are colored in red. SNPs are arranged according to the chromosome position.

For severe COVID-19 with VTE, the strongest association signals were localized to the *LINC00970* gene (index SNP: rs114101204 for severe COVID-19 and VTE, Pmeta=1.44×10−37) at locus 1q24.2. For COVID-19 hospitalization and SARS-CoV-2 infection with VTE, the strongest association signals were localized on or near the *ABO* gene (index SNP: rs11244061 for COVID-19 hospitalization, Pmeta=3.64×10−28; rs550057 for SARS-CoV-2 infection, Pmeta=1.56×10−51) at locus 9q34.2.

### Fine-mapping and colocalization analysis identify shared causal variants

For each of the shared loci of VTE with COVID-19, Supplementary Tables 5-7 listed all SNPs within 500 kb of these variants in the 99% credible sets. Co-localization analysis shows that those colocalized genetic loci for VTE and severe COVID-19 were located in *ABO* (index SNPs: rs11244061 and rs138683771), *ADAMTS13* (index SNP: rs149181677), the colocalized loci for VTE and COVID-19 hospitalization were located in *ABO* (index SNPs: rs11244061and rs8176686) and *ADAMTS13* (index SNP: rs149181677), the colocalized loci for VTE and SARS-CoV-2 infection were located in *ABO* (index SNPs: rs550057, rs9411367 and rs71503180) and *FUT2* (index SNP: rs492602) (Table 1).

### GTEx tissue-specific expression analysis (TSEA) and over-representation enrichment analysis of shared genes

The identified shared genes for severe COVID-19 and COVID-19 hospitalization with VTE were significantly enriched for expression in the lung tissue; however, no highly enriched tissues were found for SARS-CoV-2 infection and VTE (Figure 4). Gene ontology (GO) analysis highlighted several significant shared biological processes between three COVID-19 traits and VTE, such as “calcium-mediated signaling”, “second-messenger-mediated signaling”, “chemokine-mediated signaling pathway”, “response to chemokine”, “response to type I interferon” (Supplementary Tables 8-10).

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/06/01/2022.05.21.22275413/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2022/06/01/2022.05.21.22275413/F4)

Figure 4. GTEx tissue enrichment analysis for expression of shared significant genes (meta-analysis P-value < 5×10−6) for (a) VTE and severe COVID-19, (b) VTE and COVID-19 hospitalization, (c) VTE and SARS-CoV-2 infection.
P-values of Fisher’s exact test after Benjamin-Hochberg correction are presented in – log10 scale. Red represents significant tissue enrichment (Lung, P-value= 7.16×10−4 for severe COVID-19 and VTE, P-value= 8.20×10−3 for COVID-19 hospitalization and VTE).

### Bi-directional MR analysis

We found that genetically predicted VTE was positively associated with higher risk of all three COVID-19 related traits (OR= 1.08 for severe COVID-19, 95% CI: 1.04-1.12, P-value= 3.22×10−5; OR=1.08 for COVID-19 hospitalization, 95% CI: 1.05-1.11, P-value =9.09×10−10; OR=1.05 for SARS-CoV-2 infection, 95% CI: 1.04-1.07, P-value= 1.02×10−9) (Figure 5). In the reverse direction, we found that a positive association of genetically predicted SARS-CoV-2 infection with higher risk of VTE (OR=1.80, 95% CI: 1.41-2.29, P-value=1.98×10−6), and null associations of genetically predicted severe COVID-19 or COVID-19 hospitalization with VTE (Figure 6). These results are robust to different MR methods (Supplementary Figures 1 and 2).

![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/06/01/2022.05.21.22275413/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2022/06/01/2022.05.21.22275413/F5)

Figure 5. Mendelian randomization analysis on the association of genetically predicted VTE with the risk of COVID-19 related traits using inverse variance weighting (IVW).
The estimates are presented as odds ratios (OR) with 95% confidence intervals (CI).

![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/06/01/2022.05.21.22275413/F6.medium.gif)

[Figure 6.](http://medrxiv.org/content/early/2022/06/01/2022.05.21.22275413/F6)

Figure 6. Mendelian randomization analysis on the association of genetically predicted COVID-19 related traits with the risk of VTE using IVW.
The estimates are presented as odds ratios (OR) with 95% confidence intervals (CI). An arrow represents the estimate is out of boundary in this direction.

## Discussion

In this study, we comprehensively investigated the shared genetic etiology between three COVID-19 related traits and VTE based on the latest data from COVID-19 HGI and UK Biobank. For the first time, we identified 8 shared loci for severe COVID-19, 11 loci for COVID-19 hospitalization and 7 loci for SARS-CoV-2 infection with VTE, with mapped genes involved in coagulation and immune function. Co-localization analysis identified eight shared genetic loci in *ABO, ADAMTS13* and *FUT2* genes. Enrichment analysis suggested that those gene expressions were enriched mainly in the lung tissue, and supported pathways related to coagulation and immunity. Moreover, we examined their causal relationship using bi-directional MR, which suggested that VTE may increase the risk of severe COVID-19, COVID-19 hospitalization and SARS-CoV-2 infection, and interestingly, SARS-CoV-2 infection may increase the risk of VTE.

Our study found genetic correlations between VTE and COVID-19. LDSC analysis showed a significant positive genetic correlation between VTE and COVID-19 hospitalization. Although LDSC analysis did not find significant overall genetic correlation between severe COVID-19 and VTE, however partitioned LDSC analysis found severe COVID-19 and VTE were positively correlated in some functional categories, including TFBS, Fetal DHS, H3K4me1, H3K9ac and H3K27ac, which are associated with the control of transcription and the status of *cis*-regulatory elements such as promoters and enhancers within gene regulatory regions.17-20

We identified shared genetic loci between VTE and COVID-19 using cross-trait meta-analysis and colocalization analysis. We identified shared loci in 13 genes (*LINC00970, ABO, ADAMTS13, TSPAN15, NAPSA, IFNAR2, F5, WNT3, DNAJB4, SLC39A8, OBSCN, SLC6A20, FUT2*) and co-localized genetic loci in three genes (*ABO, ADAMTS13* and *FUT2*), involved in coagulation and immune function. Consistent with previous report, we confirmed the role of *ABO, NAPSA, IFNAR2* in severe COVID-19 and COVID-19 hospitalization,3,4,7,21,22 *WNT3* in COVID-19 hospitalization23, and *SLC6A20* in SARS-CoV-2 infection4. Notably, we found several novel genes, including *ADAMTS13, TSPAN15, WNT3, DNAJB4, SLC39A8, OBSCN*, which have not been reported yet. These novel findings provide new insights into genetic etiology of both VTE and COVID-19 related traits, and may help pinpoint therapeutic targets.

Interestingly, we found several genes related to coagulation. In cross-trait meta-analysis, we found a strong shared genetic signal for severe COVID-19 and VTE located at locus 1q24.2 in *LINC00970* close to *NME7* gene. Locus 1q24.2 has been found to be associated with VTE previously,24 but has not been linked to severe COVID-19. This locus (1q24.2) also includes shared genetic variant in *F5* gene. *F5* gene encodes plasma procoagulant factor (F)V, a coagulation factor involved in thrombin activation, and its impaired downregulation is key to thrombosis.25 Similarly, *ABO, FUT2*, and *ADAMST13* are well-established VTE related genes,24,26 and were found in the colocalization analysis. Our findings are consistent with previous genetic studies showing *ABO* gene is associated with severe COVID-19 and SARS-CoV-2 infection,3-5,15 possibly by regulating thrombosis.3,27 In line with the role of *ABO* gene, we found that *FUT2* gene, a fucosyltransferase gene involved in ABO blood group antigen synthesis, was associated with both VTE and SARS-CoV-2 infection. We also found a novel shared gene, *ADAMST13*, which encodes protein ADAMST13. The function of ADAMTS13 might be affected by ABO blood group,28 and a previous MR study showed genetically predicted ADAMST13 is associated with severe COVID-19.29

Moreover, several shared genes are involved in immune function. The signals shared by severe COVID-19 and COVID-19 hospitalization with VTE were also mapped to the *IFNAR2* gene at 21q22.11 and *NAPSA* gene at 4q24. *IFNAR2* gene encodes interferon receptor subunit that mediates the early host immune response to viral infection.7 *NAPSA* gene is associated with damage-associated transient progenitors promoted by inflammation.21,30 Other shared genes have been known to be associated with innate antiviral defenses (*SLC6A20*),3,31 immune cell signaling regulation (*WNT3, TSPAN15)*32,33, inflammatory lung injury *(SLC39A8)*34 or tumor activation (*TSPAN15, DNAJB4, OBSCN*).35-37 These results imply that immune function may be related to both COVID-19 and VTE.

The findings from enrichment analysis also suggest immune function might be involved in the shared etiology. In the enrichment analysis, we found several shared pathways related to immune function, such as calcium-mediated signaling and chemokine/interferon related response. Previous studies suggested that calcium signaling may play a role in hemostasis and thrombosis.38,39 Calcium signaling is also of paramount importance in immune cells,40,41 and chemokine and interferon play a major role in activating host immune and inflammatory responses.42-44 Our TSEA reported that shared genes for severe COVID-19 and COVID-19 hospitalization with VTE were mainly enriched for gene expressions in the lung tissue. Consistently, pulmonary vascular endothelial injury and immunothrombosis (most of which occur within the lung microvessels) are key drivers of severe events after SARS-CoV-2 infection, such as acute respiratory failure and ARDS.45,46 So, these evidence suggests that the common pathways between COVID-19 and VTE may relate to immunity, endothelial cell function, and coagulation.

Our study also provided timely evidence regarding the causal relationship between VTE and COVID-19 traits. Interestingly, our MR analysis showed a positive association between genetic susceptibility for VTE and the risk of severe COVID-19, COVID-19 hospitalization and SARS-CoV-2 infection, as well as a positive association of SARS-CoV-2 infection with VTE. Despite the statistical significance, our finding still needs to be interpreted with caution as we cannot exclude that this positive association may have occurred by chance. However, the positive association was also shown in a previous MR study using COVID-19 HGI release 5 data showing that genetically predicted VTE was associated with higher risk of COVID-19 hospitalization and SARS-CoV-2 infection.14 Using COVID-19 HGI release 7 with doubled sample size of release 5, we added by showing a positive association of genetically predicted VTE with higher risk of severe COVID-19.

Our study delivers an important information that VTE possibly increase the risk of COVID-19, which matters in the consideration of COVID-19 management strategies, especially for people with VTE. For example, taking vaccination to lower the risk for severe COVID-19 in people with VTE. Meanwhile, our MR analysis suggested that genetically predicted SARS-CoV-2 infection was associated with a higher risk of VTE. Partly consistent with our results, in a cohort of 153,760 individuals with COVID-19, as well as two sets of control cohorts with 5,637,647 (contemporary controls) and 5,859,411 (historical controls) individuals, people with SARS-CoV-2 infection had higher risk of cardiovascular events, including VTE.47 Nevertheless, a replication study of MR analysis in a larger GWAS of VTE will be worthwhile.

Our study is the first to use large-scale genetic data to explore the shared genetic architecture between COVID-19 related traits and VTE, providing timely evidence and novel insights into the genetic etiology between them. The GWAS summary-level data for COVID-19 was obtained from the COVID19 HGI release 7 summary statistics, the largest and latest GWAS of COVID-19 available, which doubles the sample size from the previous version, and improves study power.

We also acknowledge several limitations of our study. First, the GWAS data for COVID-19 and VTE used in this study were derived from the European population, so the associations may not be generalizable to other ancestries. Second, although GWAS summary statistics conducted study-specific quality control, misclassification of COVID-19 might exist. Third, the summary statistics limit us to assess sex and age-specific genetic effects. Fourth, the GWAS of COVID-19 related traits might be conducted at different time periods, so there might be differences in the SARS-CoV-2 infection strain. However, our study does not aim to assess the association with specific strain of SARS-CoV-2 infection. Finally, although our study provides evidence of genetic correlation and genetic overlap between COVID-19 and VTE, the underlying biological mechanisms are still unclear, and further studies are still needed for validation.

In conclusion, our findings provided novel evidence of genetic correlations between severe COVID-19, COVID-19 hospitalization, SARS-CoV-2 infection and VTE, and highlighted their common genetic architecture, with shared genes closely related to coagulation and immunity. Our work contributes to the understanding of COVID-19 and VTE etiology and provides new insights into the prevention and comorbidity management of COVID-19.

## Methods

### Study population

The GWAS summary statistics for COVID-19 of European ancestry were provided by the COVID-19 host genetics initiative round 7 ([https://www.covid19hg.org/results/](https://www.covid19hg.org/results/), release date: April 08, 2022). We included three COVID-19 related traits: (1) Severe COVID-19, defined as COVID-19-confirmed individuals with very severe respiratory symptoms or those who died from the disease (up to 13,769 cases and 1,072,442 controls); (2) COVID-19 hospitalization defined as individuals who were hospitalized for related infection symptoms, with laboratory-confirmed SARS-CoV-2 infection (up to 32,519 cases and 2,062,805 controls); (3) SARS-CoV-2 infection defined as all individuals who reported positive (laboratory diagnosis, physician diagnosis or self-report) for SARS-CoV-2 infection (up to 122,616 cases and 2,475,240 controls).

For VTE, we used summary statistics on the GWAS of VTE (3,900 cases and 369,592 controls of European ancestry) from the UK Biobank. Participating cohorts of COVID-19 and VTE were described in detail in the Supplementary Note.

### Linkage disequilibrium score regression (LDSC) analysis

LDSC analysis was conducted to assess the heritability for a single trait and genetic correlations between two traits. The analysis was conducted using the LDSC software based on the GWAS summary statistics. LD scores of 1000G European ancestry was used as reference.48 An estimate of the heritability or genetic correlation can be obtained by regressing the χ2 statistics or the products of z-scores on LD scores, respectively. LDSC can also correct for the inflation of test statistics caused by polygeneicity.49

### Partitioned LDSC analysis

We performed partitioned LDSC analysis to estimate the genetic correlation between two traits within each of the following 12 functional categories18: conserved region, DNaseI digital genomic footprinting region (DGF), DNase I hypersensitivity sites (DHS), fetal DHS, H3K4me1, H3K4me3, H3K9ac and H3K27ac, intron region, super enhancers, transcription factor-binding site (TFBS) and transcribed region. The re-calculated LD scores of the SNPs classified in each specific annotation category allowed us to find out which functional categories accounted for the majority of the overall genetic correlation.

### Multi-Trait Analysis of GWAS (MTAG)

We applied MTAG to identify novel loci with strong signals for COVID-19, and to detect shared genetic variants for VTE and three COVID-19 traits.50-52 MTAG can improve the effect estimates for each COVID-19 trait by incorporating the weighted sum of GWAS estimates for VTE. We also used MTAG to conduct cross-trait meta-analysis, which utilizes sample size-weighted, fixed-effect model together with genetic covariance modeling from all sources to combine evidence of the association between individual variants for VTE and COVID-19.

### Fine-mapping and co-localization analysis

For the significant shared loci between two traits, we used the Bayesian fine-mapping algorithm to identify a 99%-credible set of causal variants for each of the shared loci within 500kb of the index SNP.53 Then, we conducted the co-localization analysis to check whether COVID-19 and VTE association signals co-localized at shared loci, by calculating the probability that two traits shared the same causal variant (P(H4)).54 If P(H4) is greater than 0.5 for any SNP, we labelled it as a co-localized genetic variant.

### Tissue-specific expression analysis (TSEA) and over-representation enrichment analysis

In order to test whether the shared genes are overly expressed in a specific tissue, we performed TSEA with the HUGO Gene Nomenclature Committee (HGNC) name of genes that correspond with the shared loci in cross-trait meta-analysis. We used the R package ‘deTS’, which uses GTEx RNA-seq data and the ENCODE panel as a reference panel, and calculates the corresponding z-score for each tissue.55 We also used the WEB-based GEne SeT AnaLysis Toolkit to further assess the overrepresented enrichment of the same shared genes in Gene ontology (GO) biological process.56

### Mendelian randomization (MR) analysis

We conducted a bidirectional MR analysis to assess the causal association between each COVID-19 related trait and VTE. Genetic instruments for VTE were obtained from a published genetic instruments published in a large GWAS of VTE ((23,151 cases, 553,439 controls).57 The GWAS provided 297 SNPs (P-value<1×10−5, r2<0.2) as instrument for VTE (Supplementary Table 11).57 The genetic associations of these SNPs with COVID-19 related traits were obtained from COVID-19 HGI release 7 GWAS meta-analysis, as shown in study population. Using the same process of identifying instruments, based on COVID-19 HGI release 7 GWAS meta-analysis we obtained 221, 282 and 177 SNPs as instruments for severe COVID-19, COVID-19 hospitalization and SARS-CoV-2 infection respectively (Supplementary Tables 12-14). The genetic associations of these SNPs with VTE were obtained from the GWAS of VTE (3,900 cases and 369,592 controls of European ancestry) in the UK Biobank. We used inverse variance weighted (IVW) analysis in the main analysis,58 and used several other MR methods that are robust to pleiotropy, including the weighted median,59 MR-Egger,60 MR-PRESSO,61 and MR-RAPS,62 in sensitivity analysis. The details of methods were provided in the Supplemental Note.

## Data Availability

Genetic associations with VTE were from the UK biobank GWAS results provided by Lee Lab (https://www.leelabsg.org/resources). Genetic associations with severe COVID-19, COVID-19 hospitalization and SARS-CoV-2 infection was obtained from COVID-19 host genetics consortium GWAS meta-analyses round 7, downloaded from https://www.covid19hg.org/results/r7/.

[https://www.covid19hg.org/results/](https://www.covid19hg.org/results/) 

[https://www.leelabsg.org/resources](https://www.leelabsg.org/resources) 

## Statements and Declarations

### Funding

The authors reported no funding received for this study.

### Competing Interests

The authors have no competing interests.

### Author Contributions

The manuscript was drafted by XH, with great help of JVZ and ZL (edited by JYYW and ZLL). Data analysis was conducted by MY and PT. XH and MY interpreted the results, with the help of JYYW, ZLL, ZL and JVZ. ZL and JVZ developed the study conception, designed the analysis, and critically revised the manuscript for important intellectual content. All authors read and approved the final manuscript.

### Ethics approval

The study is an analysis using publicly available summary data that does not require ethical approval.

## Data availability

Genetic associations with VTE were from the UK biobank GWAS results provided by Lee Lab ([https://www.leelabsg.org/resources](https://www.leelabsg.org/resources)). Genetic associations with severe COVID-19, COVID-19 hospitalization and SARS-CoV-2 infection were obtained from COVID-19 host genetics consortium GWAS meta-analyses round 7, downloaded from [https://www.covid19hg.org/results/r7/](https://www.covid19hg.org/results/r7/).

## Acknowledgments

This research was conducted using the summary statistics from COVID-19 host genetics consortium and UK Biobank. The authors would like to thank all participants in the study and investigators for sharing the valuable data.

## Footnotes

*   We fixed some typos and polished writing.

*   Received May 21, 2022.
*   Revision received May 31, 2022.
*   Accepted June 1, 2022.


*   © 2022, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## Reference

1.  1.Malas, M.B. et al. Thromboembolism risk of COVID-19 is high and associated with a higher risk of mortality: A systematic review and meta-analysis. eClinicalMedicine 29(2020).
    
    

2.  2.Zhou, F. et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. Lancet 395, 1054–1062 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30566-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

3.  3.Severe Covid, G.G. et al. Genomewide Association Study of Severe Covid-19 with Respiratory Failure. N Engl J Med 383, 1522–1534 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/nejmoa2020283&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

4.  4.Niemi, M.E.K. et al. Mapping the human genetic architecture of COVID-19. Nature 600, 472–477 (2021).
    
    

5.  5.Shelton, J.F. et al. Trans-ancestry analysis reveals genetic and nongenetic associations with COVID-19 susceptibility and severity. Nat Genet 53, 801–808 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-021-00854-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

6.  6.Kousathanas, A. et al. Whole genome sequencing reveals host factors underlying critical Covid-19. Nature (2022).
    
    

7.  7.Pairo-Castineira, E. et al. Genetic mechanisms of critical illness in COVID-19. Nature 591, 92–98 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-03065-y&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33307546&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

8.  8.Middeldorp, S. et al. Incidence of venous thromboembolism in hospitalized patients with COVID-19. J Thromb Haemost 18, 1995–2002 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10;1111/jth.14888&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

9.  9.Tang, N., Li, D., Wang, X. & Sun, Z. Abnormal coagulation parameters are associated with poor prognosis in patients with novel coronavirus pneumonia. J Thromb Haemost 18, 844–847 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/jth.14768&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

10. 10.Driggin, E. et al. Cardiovascular Considerations for Patients, Health Care Workers, and Health Systems During the COVID-19 Pandemic. J Am Coll Cardiol 75, 2352–2371 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0735-1097(20)32979-X&link_type=DOI) 

11. 11.Mehta, P. et al. COVID-19: consider cytokine storm syndromes and immunosuppression. Lancet 395, 1033–1034 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30628-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

12. 12.Cuker, A. et al. American Society of Hematology living guidelines on the use of anticoagulation for thromboprophylaxis in patients with COVID-19: July 2021 update on post-discharge thromboprophylaxis. Blood Adv, 2021005945 (2021).
    
    

13. 13.Lawlor, D.A., Harbord, R.M., Sterne, J.A., Timpson, N. & Davey Smith, G. Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat Med 27, 1133–63 (2008).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/sim.3034&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17886233&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

14. 14.Xiang, Y., Chau, C.K.-L., Qiu, J., Rao, S. & So, H.-C. Exploring causal relationships between COVID-19 and cardiometabolic disorders: A bi-directional Mendelian randomization study. medRxiv, 2021.03.20.21254008 (2021).
    
    

15. 15.Roberts, G.H.L. et al. AncestryDNA COVID-19 Host Genetic Study Identifies Three Novel Loci. medRxiv, 2020.10.06.20205864 (2020).
    
    

16. 16.Thibord, F. et al. Cross-Ancestry Investigation of Venous Thromboembolism Genomic Predictors. medRxiv, 2022.03.04.22271003 (2022).
    
    

17. 17.Dunham, I. et al. An integrated encyclopedia of DNA elements in the human genome. Nature 489, 57–74 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature11247&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22955616&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000308347000039&link_type=ISI) 

18. 18.Finucane, H.K. et al. Partitioning heritability by functional annotation using genome-wide association summary statistics. Nat Genet 47, 1228–35 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3404&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26414678&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

19. 19.Trynka, G. et al. Chromatin marks identify critical cell types for fine mapping complex trait variants. Nat Genet 45, 124–130 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2504&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23263488&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

20. 20.Vinson, C., Chatterjee, R. & Fitzgerald, P. Transcription factor binding sites and other features in human and Drosophila proximal promoters. Subcell Biochem 52, 205–222 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/978-90-481-9069-0_10&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21557085&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

21. 21.Degenhardt, F. et al. New susceptibility loci for severe COVID-19 by detailed GWAS analysis in European populations. medRxiv, 2021.07.21.21260624 (2021).
    
    

22. 22.COVID-19 HGI & Ganna, A. Mapping the human genetic architecture of COVID-19: an update. medRxiv, 2021.11.08.21265944 (2021).
    
    

23. 23.Wu, L., Zhu, J., Liu, D., Sun, Y. & Wu, C. An integrative multiomics analysis identifies putative causal genes for COVID-19 severity. Genet Med 23, 2076–2086 (2021).
    
    

24. 24.Heit, J.A. et al. A genome-wide association study of venous thromboembolism identifies risk variants in chromosomes 1q24.2 and 9q. J Thromb Haemost 10, 1521–1531 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1538-7836.2012.04810.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22672568&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

25. 25.Bos, M.H.A. & Camire, R.M. Blood coagulation factors V and VIII: Molecular Mechanisms of Procofactor Activation. J Coagul Disord 2, 19–27 (2010).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21165149&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

26. 26.Lotta, L.A. et al. Next-generation sequencing study finds an excess of rare, coding single-nucleotide variants of ADAMTS13 in patients with deep vein thrombosis. J Thromb Haemost 11, 1228–39 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/jth.12291&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23648131&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

27. 27.Wang, L. et al. An atlas connecting shared genetic architecture of human diseases and molecular phenotypes provides insight into COVID-19 susceptibility. Genome Med 13, 83–83 (2021).
    
    

28. 28.McGrath, R.T. et al. Expression of terminal alpha2-6-linked sialic acid on von Willebrand factor specifically enhances proteolysis by ADAMTS13. Blood 115, 2666–73 (2010).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImJsb29kam91cm5hbCI7czo1OiJyZXNpZCI7czoxMToiMTE1LzEzLzI2NjYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wNi8wMS8yMDIyLjA1LjIxLjIyMjc1NDEzLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

29. 29.Zhou, Y. et al. Coagulation factors and the incidence of COVID-19 severity: Mendelian randomization analyses and supporting evidence. Signal Transduct Target Ther 6, 222–222 (2021).
    
    

30. 30.Melms, J.C. et al. A molecular single-cell lung atlas of lethal COVID-19. Nature 595, 114–119 (2021).
    
    

31. 31.Vuille-dit-Bille, R.N. et al. Human intestine luminal ACE2 and amino acid transporter expression increased by ACE-inhibitors. Amino Acids 47, 693–705 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00726-014-1889-6&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25534429&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

32. 32.Haseeb, M., Pirzada, R.H., Ain, Q.U. & Choi, S. Wnt Signaling in the Regulation of Immune Cell and Cancer Therapeutics. Cells 8, 1380 (2019).
    
    

33. 33.Cabañas, C., Yáñez-Mó, M. & van Spriel, A.B. Editorial: Functional Relevance of Tetraspanins in the Immune System. Front Immunol 10(2019).
    
    

34. 34.Besecker, B. et al. The human zinc transporter SLC39A8 (Zip8) is critical in zinc-mediated cytoprotection in lung epithelia. Am J Physiol Lung Cell Mol Physiol 294, L1127–36 (2008).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/ajplung.00057.2008&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18390834&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000256381000015&link_type=ISI) 

35. 35.Acun, T. et al. HLJ1 (DNAJB4) Gene Is a Novel Biomarker Candidate in Breast Cancer. OMICS J Integr Biol. 21, 257–265 (2017).
    
    

36. 36.Rajendran, B.K. & Deng, C.-X. A comprehensive genomic meta-analysis identifies confirmatory role of OBSCN gene in breast tumorigenesis. Oncotarget 8, 102263–102276 (2017).
    
    

37. 37.Sidahmed-Adrar, N. et al. Tspan15 Is a New Stemness-Related Marker in Hepatocellular Carcinoma. Proteomics 19, e1900025 (2019).
    
    

38. 38.Nesbitt, W.S. et al. Intercellular calcium communication regulates platelet aggregation and thrombus growth. J Cell Biol 160, 1151–1161 (2003).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNiIjtzOjU6InJlc2lkIjtzOjEwOiIxNjAvNy8xMTUxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDYvMDEvMjAyMi4wNS4yMS4yMjI3NTQxMy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

39. 39.Varga-Szabo, D., Braun, A. & Nieswandt, B. Calcium signaling in platelets. J Thromb Haemost 7, 1057–66 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1538-7836.2009.03455.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19422456&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

40. 40.Vig, M. & Kinet, J.-P. Calcium signaling in immune cells. Nat Immunol 10, 21–27 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ni.f.220&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19088738&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000261788800008&link_type=ISI) 

41. 41.Trebak, M. & Kinet, J.-P. Calcium signalling in T cells. Nat Rev Immunol 19, 154–169 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41577-018-0110-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30622345&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

42. 42.Moser, B. & Willimann, K. Chemokines: role in inflammation and immune surveillance. Ann Rheum Dis 63, ii84 (2004).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6ImFubnJoZXVtZGlzIjtzOjU6InJlc2lkIjtzOjE1OiI2My9zdXBwbF8yL2lpODQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wNi8wMS8yMDIyLjA1LjIxLjIyMjc1NDEzLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

43. 43.McNab, F., Mayer-Barber, K., Sher, A., Wack, A. & O’Garra, A. Type I interferons in infectious disease. Nat Rev Immunol 15, 87–103 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nri3787&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25614319&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

44. 44.Charo, I.F. & Ransohoff, R.M. The Many Roles of Chemokines and Chemokine Receptors in Inflammation. N Engl J Med 354, 610–621 (2006).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMra052723&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16467548&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000235177500010&link_type=ISI) 

45. 45.McFadyen, J.D., Stevens, H. & Peter, K. The Emerging Threat of (Micro)Thrombosis in COVID-19 and Its Therapeutic Implications. Circ Res 127, 571–587 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1161/CIRCRESAHA.120.317447&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

46. 46.Bonaventura, A. et al. Endothelial dysfunction and immunothrombosis as key pathogenic mechanisms in COVID-19. Nat Rev Immunol 21, 319–329 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41577-021-00536-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

47. 47.Xie, Y., Xu, E., Bowe, B. & Al-Aly, Z. Long-term cardiovascular outcomes of COVID-19. Nat Med 28, 583–590 (2022).
    
    

48. 48.Bulik-Sullivan, B.K. et al. LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat Genet 47, 291–5 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3211&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25642630&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

49. 49.Yang, J. et al. Genomic inflation factors under polygenic inheritance. Eur J Hum Genet 19, 807–812 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ejhg.2011.39&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21407268&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

50. 50.Turley, P. et al. Multi-trait analysis of genome-wide association summary statistics using MTAG. Nat Genet 50, 229–237 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-017-0009-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29292387&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

51. 51.Chung, J. et al. Genome-wide association study of cerebral small vessel disease reveals established and novel loci. Brain 142, 3176–3189 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/brain/awz233&link_type=DOI) 

52. 52.Guo, Y. et al. A genome-wide cross-phenotype meta-analysis of the association of blood pressure with migraine. Nat Commun 11, 3368 (2020).
    
    

53. 53.Farh, K.K. et al. Genetic and epigenetic fine mapping of causal autoimmune disease variants. Nature 518, 337–43 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature13835&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25363779&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

54. 54.Giambartolomei, C. et al. Bayesian test for colocalisation between pairs of genetic association studies using summary statistics. PLoS Genet 10, e1004383 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1004383&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24830394&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

55. 55.Pei, G., Dai, Y., Zhao, Z. & Jia, P. deTS: tissue-specific enrichment analysis to decode tissue specificity. Bioinformatics 35, 3842–3845 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btz138&link_type=DOI) 

56. 56.Wang, J., Duncan, D., Shi, Z. & Zhang, B. WEB-based GEne SeT AnaLysis Toolkit (WebGestalt): update 2013. Nucleic Acids Res 41, W77–83 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkt439&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23703215&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000323603200013&link_type=ISI) 

57. 57.Klarin, D. et al. Genome-wide association analysis of venous thromboembolism identifies new risk loci and genetic overlap with arterial vascular disease. Nat Genet 51, 1574–1579 (2019).
    
    

58. 58.Burgess, S., Butterworth, A. & Thompson, S.G. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol 37, 658–65 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/gepi.21758&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24114802&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

59. 59.Emdin, C.A., Khera, A.V. & Kathiresan, S. Mendelian Randomization. JAMA 318, 1925–1926 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2017.17219&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29164242&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

60. 60.Bowden, J., Davey Smith, G. & Burgess, S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol 44, 512–25 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyv080&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26050253&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

61. 61.Verbanck, M., Chen, C.Y., Neale, B. & Do, R. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet 50, 693–698 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0099-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29686387&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F01%2F2022.05.21.22275413.atom) 

62. 62.Zhao, Q., Wang, J., Hemani, G., Bowden, J. & Small, D.S. Statistical inference in two-sample summary-data Mendelian randomization using robust adjusted profile score. Ann Stat 48, 1742–1769, 28 (2020).