Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes COVID-19. SARS-CoV-2 relies on cellular RNA-binding proteins (RBPs) to replicate and spread, although which RBPs control its life cycle remains largely unknown. Here, we employ a multi-omic approach to identify systematically and comprehensively the cellular and viral RBPs that are involved in SARS-CoV-2 infection. We reveal that SARS-CoV-2 infection profoundly remodels the cellular RNA-bound proteome, which includes wide-ranging effects on RNA metabolic pathways, non-canonical RBPs and antiviral factors. Moreover, we apply a new method to identify the proteins that directly interact with viral RNA, uncovering dozens of cellular RBPs and six viral proteins. Amongst them, several components of the tRNA ligase complex, which we show regulate SARS-CoV-2 infection. Furthermore, we discover that available drugs targeting host RBPs that interact with SARS-CoV-2 RNA inhibit infection. Collectively, our results uncover a new universe of host-virus interactions with potential for new antiviral therapies against COVID-19.

SignificanceThere is a strong need for the development of new antiviral therapies, and this study sheds light on a host–virus interaction that is significant for a number of medically important human viruses. The study also suggests that the RNA-binding zinc finger CCCH-type containing 11A (ZC3H11A) protein takes part in a mechanism that facilitates nuclear export of mRNA, particularly under cellular stress, a mechanism that has been “hijacked” by several nuclear-replicating viruses to promote their replication. ZC3H11A is therefore a potential target for development of an antiviral therapy.

Human adenoviruses (HAds) encode for one or two highly abundant virus-associated RNAs, designated VA RNAI and VA RNAII, which fold into stable hairpin structures resembling miRNA precursors. Here we show that the terminal stem of the VA RNAs originating from Ad4, Ad5, Ad11 and Ad37, all undergo Dicer dependent processing into virus-specific miRNAs (so-called mivaRNAs). We further show that the mivaRNA duplex is subjected to a highly asymmetric RISC loading with the 3′-strand from all VA RNAs being the favored strand, except for the Ad37 VA RNAII, where the 5′-mivaRNAII strand was preferentially assembled into RISC. Although the mivaRNA seed sequences are not fully conserved between the HAds a bioinformatics prediction approach suggests that a large fraction of the VA RNAII-, but not the VA RNAI-derived mivaRNAs still are able to target the same cellular genes. Using small RNA deep sequencing we demonstrate that the Dicer processing event in the terminal stem of the VA RNAs is not unique and generates 3′-mivaRNAs with a slight variation of the position of the 5′ terminal nucleotide in the RISC loaded guide strand. Also, we show that all analyzed VA RNAs, except Ad37 VA RNAI and Ad5 VA RNAII, utilize an alternative upstream A start site in addition to the classical +1 G start site. Further, the 5′-mivaRNAs with an A start appears to be preferentially incorporated into RISC. Although the majority of mivaRNA research has been done using Ad5 as the model system our analysis demonstrates that the mivaRNAs expressed in Ad11- and Ad37-infected cells are the most abundant mivaRNAs associated with Ago2-containing RISC. Collectively, our results show an unexpected variability in Dicer processing of the VA RNAs and a serotype-specific loading of mivaRNAs into Ago2-based RISC.