Integration of risk variants from GWAS with SARS-CoV-2 RNA interactome prioritizes FUBP1 and RAB2A as risk genes for COVID-19

The role of host genetic factors in COVID-19 outcomes remains unclear despite various genome-wide association studies (GWAS). We annotate all significant variants and those variants in high LD (R(2) > 0.8) from the COVID-19 host genetics initiative (HGI) and identify risk genes by recognizing genes intolerant nonsynonymous mutations in coding regions and genes associated with cis-expression quantitative trait loci (cis-eQTL) in non-coding regions. These genes are enriched in the immune response pathway and viral life cycle. It has been found that host RNA binding proteins (RBPs) participate in different phases of the SARS-CoV-2 life cycle. We collect 503 RBPs that interact with SARS-CoV-2 RNA concluded from in vitro studies. Combining risk genes from the HGI with RBPs, we identify two COVID-19 risk loci that regulate the expression levels of FUBP1 and RAB2A in the lung. Due to the risk allele, COVID-19 patients show downregulation of FUBP1 and upregulation of RAB2A. Using single-cell RNA sequencing data, we show that FUBP1 and RAB2A are expressed in SARS-CoV-2-infected upper respiratory tract epithelial cells. We further identify NC_000001.11:g.77984833C>A and NC_000008.11:g.60559280T>C as functional variants by surveying allele-specific transcription factor sites and cis-regulatory elements and performing motif analysis. To sum up, our research, which associates human genetics with expression levels of RBPs, identifies FUBP1 and RAB2A as two risk genes for COVID-19 and reveals the anti-viral role of FUBP1 and the pro-viral role of RAB2A in the infection of SARS-CoV-2.