Mendelian randomization analysis identified genes pleiotropically associated with the risk and prognosis of COVID-19 ==================================================================================================================== * Di Liu * Jingyun Yang * Bowen Feng * Wenjin Lu * Chuntao Zhao * Lizhuo Li ## Abstract **Objectives** COVID-19 has caused a large global pandemic. Patients with COVID-19 exhibited considerable variation in disease behavior. Pervious genome-wide association studies have identified potential genetic variants involved in the risk and prognosis of COVID-19, but the underlying biological interpretation remains largely unclear. **Methods** We applied the summary data-based Mendelian randomization (SMR) method to identify genes that were pleiotropically associated with the risk and various outcomes of COVID-19, including severe respiratory confirmed COVID-19 and hospitalized COVID-19. **Results** In blood, we identified 2 probes, ILMN\_1765146 and ILMN\_1791057 tagging *IFNAR2*, that showed pleiotropic association with hospitalized COVID-19 (β [SE] = 0.42 [0.09], *P* = 4.75×10−06 and β [SE] = –0.48 [0.11], *P* = 6.76×10−06, respectively). Although no other probes were significant after correction for multiple testing in both blood and lung, multiple genes as tagged by the top 5 probes were involved in inflammation or antiviral immunity, and several other tagged genes, such as *PON2* and *HPS5*, were involved in blood coagulation. **Conclusions** We identified *IFNAR2* and other potential genes that could be involved in the susceptibility or prognosis of COVID-19. These findings provide important leads to a better understanding of the mechanisms of cytokine storm and venous thromboembolism in COVID-19 and potential therapeutic targets for the effective treatment of COVID-19. Key words: * coronavirus disease 2019 * *IFNAR2* * gene expression quantitative trait loci * summary Mendelian randomization ## Introduction Coronavirus disease 2019 (COVID-19), which is caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has created a large global pandemic and poses a serious threat to public health (1, 2). As of August 7, 2020, there were more than 19.3 million confirmed cases worldwide, with the total deaths exceeding 719,830 (3). SARS-CoV-2 is a highly pathogenic and transmissible coronavirus that primarily spreads through respiratory droplets and close contact (4). Seeking solutions to control the spread of COVID-19 and exploring effective treatments are of utmost importance to address the global challenge posed by COVID-19. Therefore, there is pressing urgency to further identify the pathological mechanisms underlying COVID-19. Patients with COVID-19 exhibited considerable variation in disease behavior. Recently, genome-wide association studies (GWAS) have been performed to identify genetic variants associated with diagnosis of and prognosis of COVID-19 (5, 6), but biological interpretation of their findings remains largely unclear. Previous research found that approximately 88% of trait-associated genetic variants detected by GWAS resided in non-coding regions of the genome and might have regulatory functions on gene expression (7). In the context of COVID-19 research, it is therefore important to explore genes whose expressions were pleiotropically/potentially causally associated with susceptibility or the development of COVID-19. Different from conventional randomized controlled trials (RCTs), Mendelian randomization (MR), which uses genetic variants as the proxy to randomization (8), is a promising tool to search for pleotropic/potentially causal effect of an exposure (e.g., gene expression) on the outcome (e.g., COVID-19 susceptibility). MR minimizes confounding and reverse causation that are commonly encountered in traditional association studies (8, 9), and has been successful in identifying gene expression sites or DNA methylation loci that are pleiotropically/potentially causally associated with various phenotypes, such as cardiovascular diseases, BMI, and rheumatoid arthritis (10-13). In this paper, we applied the summary data-based MR (SMR) method integrating summarized GWAS data for COVID-19 and cis- eQTL (expression quantitative trait loci) data to prioritize genes that are pleiotropically/potentially causally associated with the risk and prognosis of COVID-19. ## Methods ### Data sources #### eQTL data In the SMR analysis, cis-eQTL genetic variants were used as the instrumental variables (IVs) for gene expression. We performed SMR analysis for gene expression in blood and lung separately. For blood, we used the CAGE eQTL summarized data(14), which included 2765 participants. For lung, we used the V7 release of the GTEx eQTL summarized data (15), which included 278 participants. The eQTL data for blood and lung can be downloaded at [https://cnsgenomics.com/data/SMR/#eQTLsummarydata](https://cnsgenomics.com/data/SMR/#eQTLsummarydata). #### GWAS data for COVID-19 The GWAS summarized data were provided by *the COVID-19 host genetics initiative* (6) and can be downloaded at [https://www.covid19hg.org/results/](https://www.covid19hg.org/results/). Three phenotypes were examined, including severe respiratory confirmed COVID-19, COVID-19 and hospitalized COVID-19. Details on definition of the phenotypes can be found in **Table S1**. The control groups were subjects from the general population without the specific phenotype, subjects who were COVID-19 negative based on prediction or self-report, or subjects who had COVID-19 without hospitalization, making a total of five comparisons: severe respiratory confirmed COVID-19 (n = 536) vs. population (n = 329,391; hereafter severe COVID-19); COVID-19 (n = 6,696) vs. population (n = 1073072; hereafter COVID-19), COVID-19 (n = 3,523) vs. COVID-19 negative (n = 36,634; hereafter COVID-19 negative); hospitalized COVID-19 (n = 3,199) vs. population (n = 897,488; hereafter hospitalized COVID-19), and hospitalized COVID-19 (n = 928) vs. COVID-19 without hospitalization (n = 2,028; hereafter COVID-19 without hospitalization). In addition, we repeated the SMR analysis for severe respiratory confirmed COVID-19 using GWAS summarized data from the Severe COVID-19 GWAS Group (5) (hereafter severe COVID-19 NEJM). The study included 1,160 patients who had severe respiratory confirmed COVID-19, and 2,205 participants from the general population without COVID-19 as the control (**Table 1**). The definition of severe respiratory confirmed COVID-19 in this study was different from the one as defined by *the COVID-19 host genetics initiative* (**Table S1**). The GWAS summarized data can be downloaded at [https://ikmb.shinyapps.io/COVID-19\_GWAS_Browser/](https://ikmb.shinyapps.io/COVID-19_GWAS_Browser/). View this table: [Table 1.](http://medrxiv.org/content/early/2020/09/03/2020.09.02.20187179/T1) Table 1. Basic information of the GWAS data and eQTL data in the SMR analyses*. ## SMR analysis MR was undertaken with cis-eQTL as the IV, gene expression as the exposure, and each of the three phenotypes as the outcome (in comparison with different control groups). MR analysis was done using the method as implemented in the software SMR. Detailed information regarding the SMR method was reported in a previous publication (10). In brief, SMR applies the principles of MR to jointly analyze GWAS and eQTL summary statistics in order to test for pleotropic association between gene expression and a trait due to a shared and potentially causal variant at a locus. The heterogeneity in dependent instruments (HEIDI) test was performed to evaluate the existence of linkage in the observed association. Rejection of the null hypothesis (i.e., *P*HEIDI< 0.05) indicates that the observed association might be due to two distinct genetic variants in high linkage disequilibrium with each other. We adopted the default settings in SMR (e.g., *P*eQTL < 5 × 10−8, minor allele frequency [MAF] > 0.01, removing SNPs in very strong linkage disequilibrium [LD, r2 > 0.9] with the top associated eQTL, and removing SNPs in low LD or not in LD [r2 < 0.05] with the top associated eQTL), and used false discovery rate (FDR) to adjust for multiple testing. Annotations of transcripts were based on the Affymetrix exon array S1.0 platforms. To functionally annotate putative transcripts, we conducted functional enrichment analysis using the functional annotation tool “Metascape” (16) for the top tagged genes in blood and lung, separately. Gene symbols corresponding to putative genes *(P*< 0.05) were used as the input of the gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Data cleaning and statistical/bioinformatical analysis was performed using R version 4.0.0 ([https://www.r-project.org/](https://www.r-project.org/)), PLINK 1.9 ([https://www.coggenomics.org/plink/1.9/](https://www.coggenomics.org/plink/1.9/)) and SMR ([https://cnsgenomics.com/software/smr/](https://cnsgenomics.com/software/smr/)). ## Results ### Basic information of the summarized data The number of cases and controls varied dramatically among different analyses. The number of probes were approximately 8,500 and 5,500 in the analysis for blood and lung, respectively. The detailed information was shown in **Table 1**. ### SMR analysis in blood Information of the top 5 probes for each phenotype was presented in **Table 2**. Using GWAS summarized data from *the COVID-19 host genetics initiative*, we identified 2 probes, ILMN\_1765146 and ILMN\_1791057 tagging *IFNAR2*, that showed pleiotropic association with hospitalized COVID-19 (β [SE] = 0.42 [0.09], *P* = 4.75×10−06 and β [SE] = –0.48 [0.11], *P* = 6.76×10−06, respectively), but no significant pleiotropic association for other phenotypes after correction for multiple testing **(Figure 1, Figure S1–5)**. The same two probes were also among the top 5 probes in the analysis of severe COVID-19 NEJM and COVID-19, although they did not reach statistical significance after correction for multiple testing (**Table 2**). In addition, we found that multiple probes tagging *TRIM5* were among the top 5 probes in the analysis of severe COVID-19 NEJM and hospitalized COVID-19 in blood. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/09/03/2020.09.02.20187179/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/09/03/2020.09.02.20187179/F1) Figure 1. Prioritizing genes around *IFNAR2* in association with hospitalized COVID-19 in blood. Top plot, grey dots represent the –log10*(P* values) for SNPs from the GWAS of hospitalized COVID-19, and rhombuses represent the –log10*(P* values) for probes from the SMR test with solid rhombuses indicating that the probes pass HEIDI test and hollow rhombuses indicating that the probes do not pass the HEIDI test. Middle plot, eQTL results in blood for two probes, ILMN\_1765146 and ILMN\_1791057 probes, tagging *IFNAR2*. Bottom plot, location of genes tagged by the probes. Highlighted in maroon indicates probes that pass SMR threshold. COVID-19, coronavirus disease 2019; *IFNAR2*, interferon alpha and beta receptor subunit 2; GWAS, genome-wide association studies; SMR, summary data-based Mendelian randomization; HEIDI, heterogeneity in dependent instruments; eQTL, expression quantitative trait loci View this table: [Table 2.](http://medrxiv.org/content/early/2020/09/03/2020.09.02.20187179/T2) Table 2. The top five probes identified in the SMR analysis for blood. GO enrichment analysis of biological process and molecular function showed that the genes tagged by the top 5 probes were involved in four GO terms, including regulation of type I interferon-mediated signaling pathway (GO: 0060338), regulation of vesicle-mediated transport (GO: 0060627), response to endoplasmic reticulum stress (GO: 0034976); and response to oxidative stress (GO:0006979). Multiple genes as tagged by the top 5 probes, including *IFNAR2, TRIM5, NLRC5, MCTP1, PTPN1, FCER1G, ATF4* and *PON2*, were involved in inflammation or antiviral immunity(**Table S2**). ### SMR analysis in lung We didn’t identify any significant pleiotropic association after correction for multiple testing **(Figure S6–11)**. Information of the top 5 probes for each phenotype was presented in **Table 3**. We found that 3 probes tagging *AP006621* were among the top probes in the analysis of severe COVID-19 NEJM (**Table 3**). In addition, multiple genes, including *C7orf25, ITGAD, FCER1G* and *ZNF589*, were tagged by the top 5 probes in both blood and lung. View this table: [Table 3.](http://medrxiv.org/content/early/2020/09/03/2020.09.02.20187179/T3) Table 3. The top five probes identified in the SMR analysis for lung. The genes tagged by the top 5 probes were involved in two GO terms, including stimulatory C-type lectin receptor signaling pathway (GO: 0002223) and leukocyte activation involved in immune response (GO: 0002366). Multiple genes as tagged by the top 5 probes, including *ARSA, FCER1G, XOSC6* and *PSMD13*, were involved inflammation or antiviral immunity (**Table S3**). ## Discussion In the present study, we integrated GWAS and eQTL data in the MR analysis to explore putative genes that showed pleiotropic/potentially causal association the susceptibility/prognosis of COVID-19. We identified 2 probes tagging *IFNAR2* showing pleiotropic association with hospitalized COVID-19 in blood. Multiple genes as tagged by the top 5 probes were involved in inflammation and antiviral immunity in both blood and lung. Several genes tagged by the top probes were in blood coagulation. Our findings provided important leads to a better understanding of the mechanism of cytokine storm and venous thromboembolism in COVID-19 and revealed potential therapeutic targets for the effective treatment of COVID-19. Interferons (IFNs) refer to a group of signaling proteins made and released by the host cells in response to viral invasion (17). There are three types of IFNs: type I IFNs (IFN-α/β), type II IFNs (IFN-γ) and type III IFNs (IFN-λ) (18-20). *IFNAR2* (interferon alpha and beta receptor subunit 2), located at 21q22.11, encodes one of the two type I IFNs (21). In the cascade of host’s response to coronavirus, IFNs play an essential role in the establishment of antiviral state and in intensifying the antiviral response (22). It was found that IFNs could have both beneficial and detrimental effect on SARS-Cov-2 replication (23). A recent retrospective study of 77 adults found that IFN-α2b treatment with or without arbidol significantly reduced the duration of detectable virus in the upper respiratory tract of COVID-19 patients (24). However, the COVID-19 Treatment Guidelines Panel recommend against the use of interferons for the treatment of patients with severe and critical COVID-19, except in a clinical trial, due to insufficient data to support the beneficial or detrimental effects of interferons (25). Because *IFNAR2* was tagged by the top 5 probes for multiple phenotypes, we think that it is likely involved in determining COVID-19 severity and could be a potential therapeutic target for the treatment of COVID-19. More studies are needed to elucidate the mechanisms underlying the dual role of IFNs during SARS-Cov-2 and whether/how *IFNAR2* is involved in this process. We found that multiple genes as tagged by the top 5 probes were involved in inflammation. *ATF4* (activating transcription factor 4) is an endoplasmic reticulum stress sensor that defends lungs via induction of heme oxygenase 1 (26). ATF4 were decreased in inflamed intestinal mucosa from patients with active Crohn’s disease or active ulcerative colitis, and ATF4 deficiency promotes intestinal inflammation in mice (27). *ATF4* was downregulated in the alveolar type II cells of the elderly, compared with the young (28). *ATF4* has more than 20 candidate downstream factors, with the majority of them being significantly downregulated in the elderly, who demonstrated compromised ATF4-dependent ability to respond to endoplasmic reticulum stress. Lung-specific delivery of ATF4-related antioxidants has the potential to work in synergy with promising antiviral drugs to further improve COVID-19 outcomes in the elderly (28). Some of the tagged genes, such as *TRIM5, NLRC5, MCTP1, PTPN1, ARSA* and *FCER1G*, are also involved antiviral immunity (29-35). These genes have not been reported to be associated with COVID-19 in previous studies. It was estimated that approximately 40% of COVID-19 patients were considered as having a high risk of venous thromboembolism risk, and 11% of them developed venous thromboembolism without prophylaxis (36). Abnormal coagulation in patients with COVID-19 was associated with an increased risk of death (37). Recent research indicated that dysregulated platelets and neutrophils cooperated to drive a systemic prothrombotic state in SARS-CoV-2 infection (38). However, the underlying mechanisms in general, and genetic contributions in particular, remain to be further explored. In our study, we found two genes, as tagged by the top five probes, that were involved in blood coagulation: *PON2* and *HPS5* in the analysis of COVID-19 without hospitalization in blood (**Table 2**) and severe COVID-19 in lung (**Table 3**), respectively. *PON2* (paraoxonase 2) encodes a member of the paraoxonase gene family, and acts as a cellular antioxidant, protecting cells from oxidative stress (39). It was shown that deregulated redox regulation in PON2 deficiency caused vascular inflammation and abnormalities in blood coagulation (40). *HPS5* encodes a protein that is involved in organelle biogenesis associated with melanosomes and platelet dense granules (41). Mutations in *HPS5* are associated with Hermansky-Pudlak syndrome type 5 (HPS5), a subtype of a series of disorders characterized by oculocutaneous albinism and prolonged bleeding (42). Although the pleotropic association was not statistically significant in the SMR analysis after correction for multiple testing (P = 7.56×10−4 for *PON2* and 4.21×10−4 for *HPS5)*, our findings suggested potential involvement of these genes in the pathogenesis of venous thromboembolism in COVID-19 patients. More studies are needed to explore the functions of these two genes in response to SARS-CoV-2 infection. Our study has some limitations. The GWAS analyses did not control confounding factors which might affect the outcome. It is also unclear whether selection of the subjects in the GWAS studies was a representative of the exposure-outcome distributions in the overall population, and therefore, the possibility of selection bias, which can affect estimation accuracy, could not be ruled out. The GWAS studies only examined the short-term effect of COVID-19 due to the limited duration of the COVID-19 pandemic, and we were unable to assess the long-term outcomes/lingering effects of COVID-19. Similarly, we could not analyze the genetic contribution of other interesting phenotypes, such as different disease behaviors among children/teens, adults and the elderly patients, and asymptomatic COVID-19, due to a lack of the corresponding GWAS summarized data. We only performed analyses using blood and lung eQTL data, more studies are needed to explore tissue-and cell-type-specific genes involved in the host responses to COVID-19 infection. Due to a lack of individual eQTL data, we could not quantify the changes in gene expression in patients with COVID-19 in comparison with the control. ## Conclusion We identified *IFNAR2* and other potential genes that could be involved in the susceptibility/prognosis of COVID-19. These findings provide important leads to a better understanding of the mechanism of cytokine storm in COVID-19 and reveals potential therapeutic targets for the effective treatment of COVID-19. ## Data Availability All data generated or analyzed during this study are included in this published article and its supplementary information files. ## Author contribution DL, JY and LL designed and registered the study. DL, BF and WL analyzed data and performed data interpretation. DL and JY wrote the initial draft and BF, WL, CZ and LL contributed writing to subsequent versions of the manuscript. All authors reviewed the study findings and read and approved the final version before submission. ## Availability of data and materials All data generated or analyzed during this study are included in this published article and its supplementary information files. ## Disclosure of Potential Conflicts of Interest No potential conflicts of interest were disclosed by the authors. ## Acknowledgements Dr. Lizhuo Li’s research was supported by the National Natural Science Foundation of China (Grant No. 81660329). The study was supported by NIH/NIA grants P30AG10161, R01AG15819, R01AG17917, R01AG36042, U01AG61356 and 1RF1AG064312–01. Di Liu was supported by China Scholarship Council (CSC 201908110339). The authors confirmed that all authors have reviewed the contents of the article being submitted, approved its contents, and validated the accuracy of the data. * Received September 2, 2020. * Revision received September 2, 2020. * Accepted September 3, 2020. * © 2020, 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. ## References 1. 1.Li Q, Guan X, Wu P, Wang X, Zhou L, Tong Y, et al. Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia. N Engl J Med. 2020 Mar 26;382(13):1199–207. PubMed PMID: 31995857. Pubmed Central PMCID: PMC7121484. Epub 2020/01/30. eng. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2001316&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31995857&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 2. 2.Thomson B. The COVID-19 Pandemic: A Global Natural Experiment. Circulation. 2020 Jul 7;142(1):14–6. PubMed PMID: 32324429. Epub 2020/04/24. eng. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32324429&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 3. 3.Dong E, Du H, Gardner L. An interactive web-based dashboard to track COVID-19 in real time. The Lancet Infectious diseases. 2020 May;20(5):533–4. PubMed PMID: 32087114. Pubmed Central PMCID: PMC7159018. Epub 2020/02/23. eng. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(20)30120-1&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32087114&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 4. 4.Coronaviridae Study Group of the International Committee on Taxonomy of Viruses. The species Severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARSCoV-2. Nature microbiology. 2020 Apr;5(4):536–44. PubMed PMID: 32123347. Pubmed Central PMCID: PMC7095448. Epub 2020/03/04. eng. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41564-020-0695-z&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32123347&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 5. 5.Ellinghaus D, Degenhardt F, Bujanda L, Buti M, Albillos A, Invernizzi P, et al. Genomewide Association Study of Severe Covid-19 with Respiratory Failure. N Engl J Med. 2020 Jun 17. PubMed PMID: 32558485. Pubmed Central PMCID: PMC7315890. Epub 2020/06/20. eng. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2020283&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32558485&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 6. 6.COVID-19 Host Genetics Initiative. 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. European journal of human genetics: EJHG. 2020 Jun;28(6):715–8. PubMed PMID: 32404885. Pubmed Central PMCID: PMC7220587. Epub 2020/05/15. eng. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41431-020-0636-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32404885&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 7. 7.Edwards SL, Beesley J, French JD, Dunning AM. Beyond GWASs: illuminating the dark road from association to function. Am J Hum Genet. 2013 Nov 7;93(5):779–97. PubMed PMID: 24210251. Pubmed Central PMCID: PMC3824120. Epub 2013/11/12. eng. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2013.10.012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24210251&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 8. 8.Davey Smith G, Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet. 2014 Sep 15;23(R1):R89–98. PubMed PMID: 25064373. Pubmed Central PMCID: PMC4170722. Epub 2014/07/30. eng. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddu328&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25064373&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000349825700013&link_type=ISI) 9. 9.Emdin CA, Khera AV, Kathiresan S. Mendelian Randomization. Jama. 2017 Nov 21;318(19):1925–6. PubMed PMID: 29164242. Epub 2017/11/23. eng. [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%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 10. 10.Zhu Z, Zhang F, Hu H, Bakshi A, Robinson MR, Powell JE, et al. Integration of summary data from GWAS and eQTL studies predicts complex trait gene targets. Nat Genet. 2016 May;48(5):481–7. PubMed PMID: 27019110. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3538&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27019110&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 11. 11.Huan T, Joehanes R, Song C, Peng F, Guo Y, Mendelson M, et al. Genome-wide identification of DNA methylation QTLs in whole blood highlights pathways for cardiovascular disease. Nat Commun. 2019 Sep 19;10(1):4267. PubMed PMID: 31537805. Pubmed Central PMCID: PMC6753136. Epub 2019/09/21. eng. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-12228-z&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31537805&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 12. 12.Porcu E, Rüeger S, Lepik K, Santoni FA, Reymond A, Kutalik Z. Mendelian randomization integrating GWAS and eQTL data reveals genetic determinants of complex and clinical traits. Nat Commun. 2019 Jul 24;10(1):3300. PubMed PMID: 31341166. Pubmed Central PMCID: PMC6656778. Epub 2019/07/26. eng. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-10936-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31341166&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 13. 13.Pavlides JM, Zhu Z, Gratten J, McRae AF, Wray NR, Yang J. Predicting gene targets from integrative analyses of summary data from GWAS and eQTL studies for 28 human complex traits. Genome Med. 2016 Aug 9;8(1):84. PubMed PMID: 27506385. Pubmed Central PMCID: PMC4979185. Epub 2016/08/11. eng. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27506385&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 14. 14.Lloyd-Jones LR, Holloway A, McRae A, Yang J, Small K, Zhao J, et al. The Genetic Architecture of Gene Expression in Peripheral Blood. Am J Hum Genet. 2017 Feb 2;100(2):228–37. PubMed PMID: 28065468. Pubmed Central PMCID: PMC5294670. Epub 2017/01/10. eng. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28065468&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 15. 15.Battle A, Brown CD, Engelhardt BE, Montgomery SB. Genetic effects on gene expression across human tissues. Nature. 2017 Oct 11;550(7675):204–13. PubMed PMID: 29022597. Pubmed Central PMCID: PMC5776756. Epub 2017/10/13. eng. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature24277&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29022597&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000412829500039&link_type=ISI) 16. 16.Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019 Apr 3;10(1):1523. PubMed PMID: 30944313. Pubmed Central PMCID: PMC6447622. Epub 2019/04/05. eng. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-09234-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30944313&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 17. 17.De Andrea M, Ravera R, Gioia D, Gariglio M, Landolfo S. The interferon system: an overview. European journal of paediatric neurology: EJPN: official journal of the European Paediatric Neurology Society. 2002;6 Suppl A:A41–6; discussion A55–8. PubMed PMID: 12365360. Epub 2002/10/09. eng. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1053/ejpn.2002.0573&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12365360&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 18. 18.Müller U, Steinhoff U, Reis LF, Hemmi S, Pavlovic J, Zinkernagel RM, et al. Functional role of type I and type II interferons in antiviral defense. Science. 1994 Jun 24;264(5167):1918–21. PubMed PMID: 8009221. Epub 1994/06/24. eng. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIyNjQvNTE2Ny8xOTE4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDkvMDMvMjAyMC4wOS4wMi4yMDE4NzE3OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 19. 19.Stetson DB, Medzhitov R. Type I interferons in host defense. Immunity. 2006 Sep;25(3):373–81. PubMed PMID: 16979569. Epub 2006/09/19. eng. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.immuni.2006.08.007&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16979569&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000240862200004&link_type=ISI) 20. 20.Broggi A, Granucci F, Zanoni I. Type III interferons: Balancing tissue tolerance and resistance to pathogen invasion. The Journal of experimental medicine. 2020 Jan 6;217(1). PubMed PMID: 31821443. Pubmed Central PMCID: PMC7037241. Epub 2019/12/11. eng. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1084/jem.20190295&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31821443&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 21. 21.Duncan CJ, Mohamad SM, Young DF, Skelton AJ, Leahy TR, Munday DC, et al. Human IFNAR2 deficiency: Lessons for antiviral immunity. Science translational medicine. 2015 Sep 30;7(307):307ra154. PubMed PMID: 26424569. Pubmed Central PMCID: PMC4926955. Epub 2015/10/02. eng. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6InNjaXRyYW5zbWVkIjtzOjU6InJlc2lkIjtzOjE0OiI3LzMwNy8zMDdyYTE1NCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA5LzAzLzIwMjAuMDkuMDIuMjAxODcxNzkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 22. 22.Schindler C, Levy DE, Decker T. JAK-STAT signaling: from interferons to cytokines. The Journal of biological chemistry. 2007 Jul 13;282(28):20059–63. PubMed PMID: 17502367. Epub 2007/05/16. eng. [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEyOiIyODIvMjgvMjAwNTkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wOS8wMy8yMDIwLjA5LjAyLjIwMTg3MTc5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 23. 23.Sa Ribero M, Jouvenet N, Dreux M, Nisole S. Interplay between SARS-CoV-2 and the type I interferon response. PLoS pathogens. 2020 Jul;16(7):e1008737. PubMed PMID: 32726355. Pubmed Central PMCID: PMC7390284. Epub 2020/07/30. eng. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.ppat.1008737&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32726355&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 24. 24.Zhou Q, Chen V, Shannon CP, Wei XS, Xiang X, Wang X, et al. Interferon-α2b Treatment for COVID-19. Front Immunol. 2020;11:1061. PubMed PMID: 32574262. Pubmed Central PMCID: PMC7242746. Epub 2020/06/24. eng. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fimmu.2020.01061&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32574262&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 25. 25.Coronavirus Disease 2019 (COVID-19) Treatment Guidelines. Coronavirus Disease 2019 (COVID-19) Treatment Guidelines 2020. Available from: [https://www.covid19treatmentguidelines.nih.gov/immune-basedtherapy/immunomodulators/interferons/](https://www.covid19treatmentguidelines.nih.gov/immune-basedtherapy/immunomodulators/interferons/). 26. 26.Duan XH, Chang JR, Zhang J, Zhang BH, Li YL, Teng X, et al. Activating transcription factor 4 is involved in endoplasmic reticulum stress-mediated apoptosis contributing to vascular calcification. Apoptosis: an international journal on programmed cell death. 2013 Sep;18(9):1132–44. PubMed PMID: 23686245. Epub 2013/05/21. eng. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10495-013-0861-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23686245&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 27. 27.Hu X, Deng J, Yu T, Chen S, Ge Y, Zhou Z, et al. ATF4 Deficiency Promotes Intestinal Inflammation in Mice by Reducing Uptake of Glutamine and Expression of Antimicrobial Peptides. Gastroenterology. 2019 Mar;156(4):1098–111. PubMed PMID: 30452920. Epub 2018/11/20. eng. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 28. 28.Abouhashem AS, Singh K, Azzazy HME, Sen CK. Is Low Alveolar Type II Cell SOD3 in the Lungs of Elderly Linked to the Observed Severity of COVID-19? Antioxidants & redox signaling. 2020 Jul 10;33(2):59–65. PubMed PMID: 32323565. Pubmed Central PMCID: PMC7307702. Epub 2020/04/24. eng. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32323565&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 29. 29.Neerincx A, Lautz K, Menning M, Kremmer E, Zigrino P, Hösel M, et al. A role for the human nucleotide-binding domain, leucine-rich repeat-containing family member NLRC5 in antiviral responses. The Journal of biological chemistry. 2010 Aug 20;285(34):26223–32. PubMed PMID: 20538593. Pubmed Central PMCID: PMC2924034. Epub 2010/06/12. eng. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEyOiIyODUvMzQvMjYyMjMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wOS8wMy8yMDIwLjA5LjAyLjIwMTg3MTc5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 30. 30.de Silva S, Wu L. TRIM5 acts as more than a retroviral restriction factor. Viruses. 2011 Jul;3(7):1204–9. PubMed PMID: 21866272. Pubmed Central PMCID: PMC3159411. Epub 2011/08/26. eng. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21866272&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 31. 31.Duhan V, Hamdan TA, Xu HC, Shinde P, Bhat H, Li F, et al. NK cell-intrinsic FcεRIγ limits CD8+ T-cell expansion and thereby turns an acute into a chronic viral infection. PLoS pathogens. 2019 Jun;15(6):e1007797. PubMed PMID: 31220194. Pubmed Central PMCID: PMC6605677. Epub 2019/06/21. eng. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31220194&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 32. 32.Sweeney NP, Regan C, Liu J, Galleu A, Dazzi F, Lindemann D, et al. Rapid and Efficient Stable Gene Transfer to Mesenchymal Stromal Cells Using a Modified Foamy Virus Vector. Molecular therapy: the journal of the American Society of Gene Therapy. 2016 Aug;24(7):1227–36. PubMed PMID: 27133965. Pubmed Central PMCID: PMC4982542. Epub 2016/05/03. eng. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27133965&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 33. 33.Wang J, Lin J, Chang Y, Li P, Yang Y. MCM3AP, a novel HBV integration site in hepatocellular carcinoma and its implication in hepatocarcinogenesis. Journal of Huazhong University of Science and Technology Medical sciences = Hua zhong ke ji da xue xue bao Yi xue Ying De wen ban = Huazhong keji daxue xuebao Yixue Yingdewen ban. 2010 Aug;30(4):425–9. PubMed PMID: 20714864. Epub 2010/08/18. eng. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20714864&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 34. 34.Ovsyannikova IG, White SJ, Larrabee BR, Grill DE, Jacobson RM, Poland GA. Leptin and leptinrelated gene polymorphisms, obesity, and influenza A/H1N1 vaccine-induced immune responses in older individuals. Vaccine. 2014 Feb 7;32(7):881–7. PubMed PMID: 24360890. Pubmed Central PMCID: PMC3922536. Epub 2013/12/24. eng. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.vaccine.2013.12.009&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24360890&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 35. 35.Farid AS, Horii Y. Modulation of paraoxonases during infectious diseases and its potential impact on atherosclerosis. Lipids Health Dis. 2012 Jul 23;11:92. PubMed PMID: 22824324. Pubmed Central PMCID: PMC3457911. Epub 2012/07/25. eng. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22824324&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 36. 36.Wang T, Chen R, Liu C, Liang W, Guan W, Tang R, et al. Attention should be paid to venous thromboembolism prophylaxis in the management of COVID-19. The Lancet Haematology. 2020 May;7(5):e362–e3. PubMed PMID: 32278361. Pubmed Central PMCID: PMC7158946. Epub 2020/04/13. eng. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32278361&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 37. 37.Tang N, Li D, Wang X, Sun Z. Abnormal coagulation parameters are associated with poor prognosis in patients with novel coronavirus pneumonia. Journal of thrombosis and haemostasis: JTH. 2020 Apr;18(4):844–7. PubMed PMID: 32073213. Pubmed Central PMCID: PMC7166509. Epub 2020/02/20. eng. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32073213&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 38. 38.Nicolai L, Leunig A, Brambs S, Kaiser R, Weinberger T, Weigand M, et al. Immunothrombotic Dysregulation in COVID-19 Pneumonia is Associated with Respiratory Failure and Coagulopathy. Circulation. 2020 Jul 28. PubMed PMID: 32755393. Epub 2020/08/07. eng. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1161/CIRCULATIONAHA.120.048488&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32755393&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 39. 39.Altenhöfer S, Witte I, Teiber JF, Wilgenbus P, Pautz A, Li H, et al. One enzyme, two functions: PON2 prevents mitochondrial superoxide formation and apoptosis independent from its lactonase activity. The Journal of biological chemistry. 2010 Aug 6;285(32):24398–403. PubMed PMID: 20530481. Pubmed Central PMCID: PMC2915675. Epub 2010/06/10. eng. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEyOiIyODUvMzIvMjQzOTgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wOS8wMy8yMDIwLjA5LjAyLjIwMTg3MTc5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 40. 40.Ebert J, Wilgenbus P, Teiber JF, Jurk K, Schwierczek K, Döhrmann M, et al. Paraoxonase-2 regulates coagulation activation through endothelial tissue factor. Blood. 2018 May 10;131(19):2161–72. PubMed PMID: 29439952. Epub 2018/02/15. eng. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImJsb29kam91cm5hbCI7czo1OiJyZXNpZCI7czoxMToiMTMxLzE5LzIxNjEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wOS8wMy8yMDIwLjA5LjAyLjIwMTg3MTc5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 41. 41.Huizing M, Hess R, Dorward H, Claassen DA, Helip-Wooley A, Kleta R, et al. Cellular, molecular and clinical characterization of patients with Hermansky-Pudlak syndrome type 5. Traffic (Copenhagen, Denmark). 2004 Sep;5(9):711–22. PubMed PMID: 15296495. Epub 2004/08/07. eng. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1600-0854.2004.00208.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15296495&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) 42. 42.Huizing M, Anikster Y, Gahl WA. Hermansky-Pudlak syndrome and related disorders of organelle formation. Traffic (Copenhagen, Denmark). 2000 Nov;1(11):823–35. PubMed PMID: 11208073. Epub 2001/02/24. eng. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1034/j.1600-0854.2000.011103.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11208073&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F03%2F2020.09.02.20187179.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000165627400003&link_type=ISI)