# Characterization and functional interrogation of the SARS-CoV-2 RNA interactome

Athéna Labeau, Luc Fery-Simonian, Alain Lefevre-Utile, Marie Pourcelot, Lucie Bonnet-Madin, Vassili Soumelis, Vincent Lotteau, Pierre-Olivier Vidalain, Ali Amara, Laurent Meertens

<figure><img src="image_1.png" /></figure>

To cite this version:

Athéna Labeau, Luc Fery-Simonian, Alain Lefevre-Utile, Marie Pourcelot, Lucie Bonnet-Madin, et al. Characterization and functional interrogation of the SARS-CoV-2 RNA interactome. Cell Reports, 2022, 39 (4), pp.110744. 10.1016/j CELREP.2022.110744. hal-03834949

HAL Id: hal-03834949

https://u-paris.hal.science/hal-03834949

Submitted on 4 Nov 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Cell Reports

# Characterization and functional interrogation of the SARS-CoV-2 RNA interactome

## Graphical abstract

<figure><img src="image_3.png" /></figure>

## Authors

Athéna Labeau, Luc Fery-Simonian,
Alain Lefevre-Utile, ...,
Pierre-Olivier Vidalain, Ali Amara,
Laurent Meertens

## Correspondence

ali.amara@inserm.fr (A.A.),
laurent.meertens@inserm.fr (L.M.)

## In brief

Labeau et al. combine a ChIRP-MS approach with a functional interrogation screen to provide an unbiased and comprehensive view of the host cell factors that interact with the viral genome and regulate SARS-CoV-2 infection. They show that these host dependency RNA-binding proteins (RBPs) represent pharmacological targets against SARS-CoV-2.

## Highlights

* ChIRP-MS of SARS-CoV-2 RNA identifies cellular and viral RNA interactors

* siRNA perturbation screen identifies both SARS-CoV-2 restriction and dependency factors

* RBPs HNRNPA2B1, ILF3, QKI, and SFPQ bind SARS-CoV-2 RNA to promote its replication

* Pharmacological inhibition of some of these RBPs can abolish SARS-CoV-2 infection

<figure><img src="image_4.png" /></figure>
Report

Characterization and functional interrogation
of the SARS-CoV-2 RNA interactome

Athéna Labeau,<fnref n="1" /> Luc Fery-Simonian,<fnref n="1" /> Alain Lefevre-Utile,<fnref n="2" /> Marie Pourcelot,<fnref n="1" /> Lucie Bonnet-Madin,<fnref n="1" /> Vassili Soumelis,<fnref n="2" />
Vincent Lotteau,<fnref n="3" /> Pierre-Olivier Vidalain,<fnref n="3" /> Ali Amara,<fnref n="1" />,\* and Laurent Meertens<fnref n="1" />,<fnref n="4" />,\*

SUMMARY

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the COVID-19 pandemic, which has led to a devastating global health crisis. The emergence of variants that escape neutralizing responses emphasizes the urgent need to deepen our understanding of SARS-CoV-2 biology. Using a comprehensive identification of RNA-binding proteins (RBPs) by mass spectrometry (ChIRP-MS) approach, we identify 107 high-confidence cellular factors that interact with the SARS-CoV-2 genome during infection. By systematically knocking down their expression in human lung epithelial cells, we find that the majority of the identified RBPs are SARS-CoV-2 proviral factors. In particular, we show that HNRNPA2B1, ILF3, QKI, and SFPQ interact with the SARS-CoV-2 genome and promote viral RNA amplification. Our study provides valuable resources for future investigations into the mechanisms of SARS-CoV-2 replication and the identification of host-centered antiviral therapies.

INTRODUCTION

Severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2), which belongs to the *Coronaviridae* family, has been identified as the causative agent of the ongoing pandemic of coronavirus disease 2019 (COVID-19) (Zhu et al., 2020). One year after the emergence of SARS-CoV-2, several vaccines in Phase III clinical trials (Poland et al., 2020) received emergency approval. Most of these vaccines are primarily designed to elicit immune responses to the viral spike protein. However, several SARS-CoV-2 variants of concern (VOCs) have begun to emerge, which contain mutations in the spike protein that increase transmissibility and reduce susceptibility to natural or vaccine-induced neutralizing immunity (i.e., breakthrough infection) (Cele et al., 2021; Garcia-Beltran et al., 2021; Liu et al., 2021a, 2021b; Zhou et al., 2021). Thus, a deeper understanding of SARS-CoV-2 replication mechanisms is required to develop additional therapeutic strategies aimed at preventing the emergence of new SARS-CoV-2 variants with an increased rate of reinfection, transmissibility, and disease severity.

The genome of SARS-CoV-2 is composed of a single-stranded positive RNA with a length of approximately 30 kb (Zhou et al., 2020). In the cytoplasm of infected cells, the viral RNA is translated by the host machinery in 2 overlapping replicase polyproteins that undergo subsequent viral protease-mediated processing to generate 16 non-structural proteins (NSPs) (de Wilde et al., 2018). These NSPs assemble and recruit poorly

characterized host cell factors to form the replication and trans-
scription complex (RTC), which localizes in a virus-induced
network of double-membrane vesicles (DMVs) (Cortese et al.,
2020; Eymieux et al., 2021). Within the RTC, the RNA-
dependent-RNA polymerase (RdRp) complex uses the genomic
RNA as a template to generate negative-strand intermediates
that synthesize both progeny genomic and subgenomic viral
RNAs (sgvRNAs) (Snijder et al., 2016). SgvRNAs are then trans-
lated in 4 structural (S, E, M, and N) and 7 accessory proteins
(ORF3a-9b) (Kim et al., 2020). The viral genomic RNA associated
with the nucleocapsid (N) forms the ribonucleocapsid, which
together with the other structural proteins, drives the assembly
of new viral particles that bud in the lumen of the endoplasmic re-
ticulum-Golgi intermediate compartment (ERGIC).

Like all viruses, SARS-CoV-2 heavily relies on cellular proteins
to accomplish its infectious life cycle. This host dependency rep-
resents a potential Achilles' heel that could be exploited to
develop new host-centered approaches to treat SARS-CoV-2
infection. In this context, affinity purification combined with
mass spectrometry (MS) has allowed the identification of
SARS-CoV-2 virus-host protein-protein interactions (PPIs) (Gor-
don et al., 2020a, 2020b; Li et al., 2020), while CRISPR-Cas9
phenotypic screens have provided valuable information about
the host genes and cellular pathways important for the SARS-
CoV-2 life cycle (Daniloski et al., 2020; Hoffmann et al., 2021;
Schneider et al., 2020; Wang et al., 2021b; Wei et al., 2020).
These works have led to the identification of potential cellular

<figure><img src="image_5.png" /></figure>

<footnote n="1">¹Université Paris Cité, INSERM U944 CNRS 7212, Institut de Recherche Saint-Louis, Hôpital Saint-Louis, 75010 Paris, France</footnote>

<footnote n="2">²Université Paris Cité, INSERM U976, Institut de Recherche Saint-Louis, Hôpital Saint-Louis, 75010 Paris, France</footnote>

<footnote n="3">³Centre International de Recherche en Infectiologie (CIRI), Univ Lyon, Inserm U1111, Université Claude Bernard Lyon 1, CNRS UMR5308, ENS de Lyon, 69007 Lyon, France</footnote>

<footnote n="4">⁴Lead contact</footnote>

<footnote>*Correspondence: ali.amara@inserm.fr (A.A.), laurent.meertens@inserm.fr (L.M.)</footnote>

<footnote>https://doi.org/10.1016/j.celrep.2022.110744</footnote>
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>B</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>C</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>D</figcaption></figure>

<figure><img src="image_10.png" /></figure>

<table><caption>E</caption><thead><tr><th colspan="11">Number of unique peptide/viral protein</th></tr><tr><th colspan="5">Non-structural proteins</th><th colspan="2">Structural proteins</th><th colspan="4">Auxiliary proteins</th></tr><tr><th>NSP1</th><th>NSP2</th><th>NSP3</th><th>NSP5</th><th>NSP8</th><th>NSP9</th><th>NSP13</th><th>NSP14</th><th>S</th><th>M</th><th>N</th><th>ORF3a</th><th>ORF6</th><th>ORF9b</th></tr></thead><tbody><tr><td>2</td><td>1</td><td>11</td><td>2</td><td>1</td><td>2</td><td>1</td><td>1</td><td>31</td><td>5</td><td>49</td><td>1</td><td>1</td><td>4</td></tr></tbody></table>

<figure><img src="image_11.png" /><figcaption>F</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 1. Purification of the SARS-CoV-2 RNA interactome**</figcaption></figure>

(A) Schematic illustrating the ChIRP-MS approach used to identify host factors bound to SARS-CoV-2 RNA.

(B) Proteins enriched by ChIRP from uninfected and SARS-CoV-2-infected cells were resolved by SDS-PAGE and visualized by silver staining.

(C) Quantification of the peptides ratio as the average of peptides count for a prey in SARS-CoV-2 divided by the average peptides count for the same prey in the control (zero counts are replaced by 1). Viral proteins are marked by red circles (n = 5 biological replicates).

(D) Top: Schematic representing unique peptide distribution along the enriched viral proteins. Bottom: Table showing the number of unique peptides for each viral protein. Data from 1 representative biological replicate are shown.

(legend continued on next page)
targets for drug repurposing, such as cholesterol homeostasis, the phosphatidylinositol 3-kinase (PI3K) complex or the sigma-1 and -2 receptors (Daniloski et al., 2020; Gordon et al., 2020a; Hoffmann et al., 2021; Wang et al., 2021b).

In this study, we undertook a comprehensive identification of RNA-binding protein by MS (ChIRP-MS) approach to establish a global map of the SARS-CoV-2 genome-associated host proteins. By combining these proteomics data with gene-silencing experiments, we identified a set of SARS-CoV-2 host dependency factors that affect infection and could be targetable for antiviral therapies.

## RESULTS AND DISCUSSION

### Characterization of the SARS-CoV-2 RNA interactome in human cells

To characterize the host factors that interact with the SARS-CoV-2 genome during infection, we performed a ChIRP-MS (Figure 1A). Briefly, we challenged 293T cells stably expressing the human angiotensin-converting enzyme 2 (ACE2) receptor (293T-ACE2) with a primary clinical SARS-CoV-2 220/95 strain (EPI_ISL_469284) isolated by our laboratory in March 2020 (Onodi et al., 2021). Forty-eight hours later, SARS-CoV-2-infected and -uninfected cells were treated with formaldehyde to crosslink RNA with RNA-binding proteins (RBPs) to preserve ribonucleoprotein (RNP) complexes. To specifically capture SARS-CoV-2 RNA, we designed 129 biotinylated antisense probes (Table S1) complementary to the full-length viral RNA, and subsequently used streptavidin beads to pull down the RNPs. Co-immunoprecipitated proteins were separated by SDS-PAGE, visualized by silver staining (Figure 1B), and analyzed by MS (Table S2). Correlation experiments showed that our experimental replicates (n = 5) are highly consistent (Pearson's correlation coefficient R > 0.97 between replicates). As expected, the SARS-CoV-2 N protein, which is the main viral RBP, was the most abundant viral protein (Figures 1B and 1C; Table S2). The replicase polyprotein pp1ab, as well as the structural M and S and the accessory ORF9b proteins were also significantly and reproducibly enriched in our experiments (Figures 1C and 1D; Table S3). We also identified 9 of the 16 pp1ab-derived NSPs (Figure 1D), which are part of the coronavirus replication/transcription complex (RTC) (V'kovski et al., 2019). Among them, NSP2-NSP16 play essential functions for RNA replication (V'kovski et al., 2021). NSP1 is known to interact both with the 40S ribosomal subunit to inhibit cellular translation, and the 5' untranslated region (UTR) of the vRNA to promote translation (Schubert et al., 2020). The M protein associates with the N protein, within the ribonucleocapsid, to favor RNA genome packaging (Narayanan et al., 2000). M also orchestrates the assembly of viral particles at the ERGIC by interacting with the other structural proteins and ensures the retention of S within this compartment, and thereby its incorporation within the virion

(Schoeman and Fielding, 2019). Thus, we speculate that enrichment of M and S may be a consequence of their association with the ribonucleocapsid during viral particle assembly. Finally, ORF9b has been shown to localize at the mitochondrial membrane and to antagonize the retinoic acid-inducible gene-I-mitochondrial antiviral-signaling protein (RIG-I-MAVS) antiviral signaling through its association with TOM70 (Gordon et al., 2020b; Jiang et al., 2020). Together, these data provide strong evidence that our ChIRP-MS approach reproducibly and specifically pulled down viral RNP complexes involved in different steps of the SARS-CoV-2 life cycle, from replication/transcription to particle assembly.

### A view of the SARS-CoV-2 RNA-host protein interactome

We used the SAINTEXpress scoring algorithm to compare the MS data from the 5 replicates (SARS-CoV-2 over Mock-infected cells) and defined a high-confidence interactome (Bayesian false discovery rate [BFDR] ≤0.05 and fold change [FC] ≥5; Table S3). We identified 107 human proteins that reproducibly associated with the SARS-CoV-2 RNA (Figure 1E). We also defined an expanded SARS-CoV-2 interactome (BFDR <0.1 and FC ≥5; Table S3) and identified 39 additional proteins. Next, we compared our high-confidence interactome with a SARS-CoV-2 interactome in which we have used SARS-CoV-2-infected cells without cross-link as negative control (Figure S1A; Tables S2 and S3). We found 78 common enriched factors and defined a high-confidence core interactome (Figures S1B and S1C). Because the majority of our high-confidence hits (78/107; 73%) were confirmed by a second approach, we could likely exclude obvious shortcomings in our enrichment calculation approach (SARS-CoV-2 over Mock-infected cells). To rule out biases resulting from infection-induced host proteome changes, we monitored the expression levels of selected hits before and after infection. We selected a panel of 9 proteins chosen among the top-ranked (e.g., G3BP1, HNRNPK, HNRNPA2B1, KHSRP, ILF3, NONO, and SFPQ) and lower-ranked (QKI, EWSR1) proteins in our MS analysis (Table S2) for which antibodies were available in our laboratory. As shown in Figure S1D, none of the studied proteins were upregulated upon SARS-CoV-2 infection of either 293T or A549 cells stably expressing ACE2.

Among the proteins included in the expanded interactome, 88% (128/146) were annotated by the Gene Ontology (GO) term "RNA binding" (GO: 0003723), suggesting a direct interaction with the viral genome. To test this hypothesis, we performed native viral RNA immunoprecipitation (vRIP) experiments on a subset of 5 selected factors (HNRNPA2B1, ILF3, QKI, SFPQ, and EWSR1). As negative controls, we included HSPB1, a heat shock protein not enriched in our ChIRP-MS, and the heterogeneous GFP protein, which unlikely binds RNA. We also added TRIM28, a protein for which no RNA-binding activity was reported, which likely was co-purified with vRNA in our

<footnote>(E) Fold change (FC; SAINTEXpress FC score) of the high-confidence and the expanded interactomes identified in our SARS-CoV-2 ChIRP-MS. (F) vRIP assays in 293T-ACE2 cells transfected with plasmid encoding the indicated proteins bearing a hemagglutinin (HA) or FLAG tag and infected with SARS-CoV-2. FLAG-HSPB1 and FLAG-GFP were used as negative controls. Enrichment of immunoprecipitated vRNA was calculated as $2^{(-\Delta\Delta\text{Ct} [\text{normalized RIP/normalized HSPB1}])}$. The data shown are means ± SDs of 2 biological replicates (with technical duplicates). Adjusted p values were calculated by the Kruskal-Wallis test with Benjamini, Krieger, and Yekutieli correction (n.s, non-significant; \*p < 0.05, \*\*p < 0.001, and \*\*\*p < 0.0001).</footnote>
ChIRP-MS through protein-mediated interaction. We confirmed that for the 5 proteins tested, there was a significant enrichment of vRNA relative to the negative controls (Figures 1F and S2A), whereas there was no vRNA enrichment in the TRIM28 sample, strongly pointing to the existence of specific interaction. These results were also confirmed in A549-ACE2 stably expressing the different proteins (Figures S2B and S2C). Altogether, these observations demonstrate that our ChIRP approach can identify proteins that associate with vRNA.

Next, we compared our RNA interactome with recently published datasets obtained by either RNA antisense purification coupled with MS (RAP-MS) or ChIRP-MS approaches in different SARS-CoV-2-infected cell lines (Flynn et al., 2021; Lee et al., 2021; Schmidt et al., 2020). We found that 36 (Schmidt et al.: 34.2%), 79 (Flynn et al.: 34.6%), and 74 (Lee et al.: 38.3%) SARS-CoV-2 RNA host factors identified in these studies were also enriched in our ChIRP-MS interactome (Table S4). Interestingly, we highlighted 58 common hits that defined a reproducible set of cellular factors binding SARS-CoV-2 RNA across 3 cell lines (293T, VeroE6, and Huh7.5; Figure S3). These encompass multiple RBPs involved in nucleic acid metabolism, including splicing, mRNA stability, and transcriptional regulation such as several heterogeneous nuclear RNPs (hnRNPs) (Geuens et al., 2016), the serine/arginine-rich splicing factor (SRSF) (Shepard and Hertel, 2009), poly(A) binding protein (PABP) (Mangus et al., 2003), and the insulin-like growth factor 2 mRNA binding protein (IGF2BP) (Bell et al., 2013). In our expanded interactome, we identified several components of the small 40S ribosome subunit (RPS3, RPS3A, RPS16, RPS19, and RPS20), the translation initiation factor 4F complex (elf4G1, elf3F, and elf4H), and translation elongation factors (EEF1D and EEF2). In agreement with other published RNA-centric SARS-CoV-2 interactomes (Flynn et al., 2021; Kamel et al., 2021; Lee et al., 2021; Schmidt et al., 2020), the eukaryotic initiation factor 4E (elf4E), the mRNA CAP binding protein, was not found. Although we cannot exclude the inability of RNA-centric methods to capture the vRNA-EIF4E complex, it is tempting to speculate that SARS-CoV-2 may either use an alternate CAP-dependent EIF4E-independent translation mechanism (de la Parra et al., 2018) to bypass host translation shutoff during the course of infection or switch from a CAP-dependent to a CAP-independent translation initiation, as proposed for dengue virus (Edgil et al., 2006). Chaperone proteins (HSP90AB1, HSPA2, CCT2, and CCT8) were also significantly enriched in our study. It is likely that they are indirectly associated with the vRNA either by contributing to nascent viral protein co-translational folding (Cabrita et al., 2016) or by interacting with viral NS proteins, as shown for dengue virus NS1 (Hafirassou et al., 2017) or hepatitis C virus NS5B or influenza PB2 viral polymerases (Fislová et al., 2010; Inoue et al., 2011). Characteristic RBPs involved in RNP granules such as paraspeckles (RBM14, NONO, MATR3, SFPQ) and stress granules (SGs) (G3BP1, G3BP2, TIA1, CAPRIN1, ATXN2L) (Nakagawa et al., 2018; Yang et al., 2020) were among the highest scoring candidates. Many paraspeckles-interacting proteins were also highly represented in our ChIRP-MS dataset (DAZAP1, EWSR1, FAM98A, FUS, RBMX, UBAP2L, as well as hnRNPA1, A1L2, H1, H3, K, R, UL1) (Naganuma et al., 2012). We also identified several components of the tRNA-splicing

ligase complex (DDX1, FAM98A, RTCB, C14ORF166/hCLE), which is involved in XBP1 cytoplasmic splicing (Filipowicz, 2014), a key regulator of the unfolded protein response that is induced by coronavirus replication (Fung and Liu, 2014). Finally, we identified several proteins with antiviral functions among the hits (ZC3HAV1, MOV10, ILF3, ADAR). For instance, ZC3HAV1 is an interferon-stimulated gene (ISG) that directly binds the viral RNA and promotes viral RNA decay (Guo et al., 2007). MOV10 was recently shown to bind bunyavirus N protein to interfere with RNP assembly and RNA synthesis (Mo et al., 2020). It was also found to interact with coronaviruses N protein and to restrict Middle East respiratory syndrome (MERS)-CoV infection (Wang et al., 2021a). These proteins contribute to the host cell response to viral infection, and their enrichment strongly suggests that SARS-CoV-2 infection triggers an innate immunity response in 293T-ACE2 cells.

Functional analyses of the vRNA-interacting proteins revealed more than 150 statistically enriched GO terms, which clustered into 23 functional groups (Figure 2A) that are shared by the 3 other published SARS-CoV-2 RNA interactome studies (Figure 2B). Among them, we found host molecules that participate in RNA metabolic process, translation, RNA stabilization, RNP complex assembly, viral process, and innate immune response (Figures 2A and 2C). Based on its predominant nuclear localization and its role in mRNA maturation, the significant enrichment of the RNA splicing GO term (p = 3.55 × 10⁻⁴³) was unexpected. Pre-mRNA splicing is critical for mRNA nucleo-cytoplasmic export, thereby contributing to the regulation of gene expression (Reed and Hurt, 2002). We infer that SARS-CoV-2 RNA may sequester host RBPs involved in splicing to hamper cellular mRNA nuclear export, thus favoring the cytoplasmic translation of viral RNAs. This hypothesis is consistent with recent findings showing that SARS-CoV-2 NPS16 disrupts mRNA splicing to antagonize innate immunity at a post-transcriptional level (Banerjee et al., 2020).

Finally, we asked whether the vRNA-binding cellular factors identified in our study also interact with viral proteins. By examining a reference coronavirus interactome (Perrin-Cocon et al., 2020), we found that 81 (55%) of the 146 hits identified by our ChIRP-MS approach are known interactors of 1 or several coronaviruses viral proteins (Figure S4). This suggests that many of the proteins identified in our SARS-CoV-2 interactome could interact with viral proteins within RNP complexes.

## Functional interrogation of the SARS-CoV-2 RBP identified in our ChIRP-MS approach

To further pinpoint the function of the identified host factors in SARS-CoV-2 infection, we analyzed the consequences of silencing their expression by RNA interference (RNAi) on virus replication. We designed a library of small interfering RNA (siRNA) pools targeting 138 RBPs identified in our study and performed a loss-of-function screen in the human lung epithelial cell line A549 stably expressing ACE2. siRNA targeting the cathepsin L protease (CTSL) and the vacuolar ATPase subunit ATP6V1B2, 2 host molecules important for SARS-CoV-2 viral entry (Daniloski et al., 2020), were used as positive controls, while a non-targeting siRNA pool (siNT) was used as negative control. siRNA transfected cells were challenged with SARS-CoV-2 for 24 h and viral
<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

Immune response

<figure><img src="image_16.png" /></figure>

Virus biology

<figure><img src="image_17.png" /></figure>

Biology of stress granule

<figure><img src="image_18.png" /></figure>

<figure><img src="image_19.png" /></figure>

<figure><img src="image_20.png" /><figcaption>**Figure 2. Biological analysis of the SARS-CoV-2 RNA interactome**</figcaption></figure>

(A) GO enrichment analysis of SARS-CoV-2 RNA interactome proteins. Circles represent the enriched function for an annotated ontology term and circle size corresponds to the number of enriched proteins within that term.

(B) Heatmap of shared GO:biological processes enriched across all SARS-CoV-2 RNA interactomes. Horizontally, biological functions have been clustered according to their enrichment significance among cohorts. Vertically, cohorts have been clustered according to their functional enrichment similarities.

(C) Selected biological process networks of SARS-CoV-2 interactome proteins.
A

<figure><img src="image_21.png" /></figure>

B

Viral titer (RT-qPCR)

<figure><img src="image_22.png" /></figure>

C

Replication (Flow cytometry)

<figure><img src="image_23.png" /></figure>

D

(Spearman r=0.742 p< 0.0001)

<figure><img src="image_24.png" /></figure>

E

<figure><img src="image_25.png" /><figcaption>Figure 3. Functional interrogation of the SARS-CoV-2 RNA interactome and compounds screening</figcaption></figure>

(A) Schematic illustrating the loss-of-function screen procedure.

(B and C) A549-ACE2 cells were transfected with an arrayed siRNA library and challenged with SARS-CoV-2 (MOI 0.05) for 24 h.

(B) Yield of viral particles released in the supernatant of infected cells was quantified by qRT-PCR and normalized to the siNT-transfected cells. The data shown are the means of 2 biological replicates (with technical duplicates).

(C) Viral replication was assessed by flow cytometry using anti-N protein monoclonal antibody (mAb) and normalized to the siNT-transfected cells. The data shown are the means of 2 biological replicates (with technical duplicates). Adjusted p values were calculated by 1-way ANOVA with Benjamini and Hochberg

(legend continued on next page)
infection was assessed either by quantifying the intracellular accumulation of the N protein by flow cytometry (viral replication) or the vRNA released in the supernatant of infected cells by qRT-PCR (release of viral particles) (Figure 3A). All readouts were normalized to the siNT-transfected cells. Transfected cells were also monitored for cell viability, and the siRNA with cytotoxic effects (below 75% of cell viability in control) were removed from the analysis (7%, 9/138). As expected, siRNA-mediated depletion of either CTSL or ATP6V1B2 resulted in a strong decrease in infection in both assays (Figures 3B and 3C). qRT-PCR experiments revealed that 45 siRNAs reduced viral particle release by at least 50%, whereas 4 siRNAs increased it by 200% (Figure 3B). Flow cytometry studies showed that 21 siRNAs decreased the percentage of infected cells by more than 50%, while 12 enhanced it by 150% (Figure 3C). By intersecting the data from both assays, we identified 25 proteins that affected viral infection by acting either as host dependency factors (HDFs) (19/25) or host restriction factors (HRFs) (6/25) (Figure 3D). This contrasts with the data recently published by Flynn et al. (2021) showing that the majority of SARS-CoV-2 RBP5 identified in their study has an antiviral function. It is likely that RNP complexes are dynamically regulated and that some RBP5 could function as pro- or antiviral host factors at different stages of the viral life cycle. We observed a high correlation coefficient ($r = 0.742$) between our qRT-PCR and flow cytometry datasets indicating that most of the identified proteins are required for genome replication (Figure 3D), and only very few act at a late step (assembly/particle egress) of the SARS-CoV-2 life cycle (Figure 3D, green circles). Interestingly, several HDFs such as ILF3, HNRNPU, HNRNPK, FAM120A, and QKI were predicted to bind the 5'- or 3'-UTR of SARS-CoV-2 RNA (Sun et al., 2021). Others, including SFPQ, DDX5, EWSR1, HNRNPK, ADAR, and PA2G4 were previously reported to promote the replication of other RNA viruses (Garcia-Moreno et al., 2019; Gélinas et al., 2011; Landeras-Bueno et al., 2011; Li et al., 2013; Oakland et al., 2013; Thompson et al., 2020; Zhou et al., 2019a), supporting their role in SARS-CoV-2 replication. Interestingly, ILF3 was shown to participate in the antiviral type I interferon response by promoting the expression of IFNB1 and ISGs during host translational shutoff (Watson et al., 2020). QKI was shown to repress host interferon response by downregulating MAVS expression at a post-transcriptional level (Liao et al., 2020). This suggests that SARS-CoV-2 has evolved several mechanisms to dampen the type I interferon response by targeting specific RBP5.

We then looked at the DrugBank, ChEMBL, DGIdb (Drug-Gene Interaction Database), and PanDrugs databases and identified 274 small molecules targeting 30 vRNA-binding proteins.

Among them, we selected camptothecin (XRCC5), cerulenin (FASN), dabigatran (HNRNPC), geldanamycin (HSP90AB1), phenethyl isothiocyanate (HNRNPK), sunitinib malate (EWSR1), and artemisinin (an antimalarial interacting with multiple RBPs), which were predicted to interact with HDF having a log2 (titer ratio) $\le 0.5$ in our siRNA screen. To test their effect on infection, A549-ACE2 cells were incubated with each compound at a concentration of 10 and 1 µM and then challenged with SARS-CoV-2. We used equivalent amounts of DMSO as negative control and similar concentrations of remdesivir, an inhibitor SARS-CoV-2 replication (Wang et al., 2020), as positive control. We monitored both the release of progeny viruses in the supernatant and cell viability 24 h post-infection. Of all of the compounds tested, we observed a significant inhibitory effect for phenethyl isothiocyanate (40% inhibition) and sunitinib malate (80% inhibition) at a concentration of 1 µM without toxicity (Figure 3E). These data suggest that the host factors identified in our ChIRP-MS study could represent promising antiviral drug targets.

## HNRPNA2B1, ILF3, QKI, and SFPQ mediate SARS-CoV-2 replication

Among the HDF identified, we were particularly interested in ILF3 and QKI (for which limited information on viral infection are available), as well as HNRNPA2B1 and SFPQ, which are likely able to associate with SARS-CoV-2 vRNA (Figures 1F and S2C). Consistent with our screen, siRNA-mediated depletion of HNRPNA2B1, ILF3, QKI, or SFPQ (Figure 4A, left panel) significantly reduced the percentage of infected cells and progeny viruses release (Figure 4A, center and right panels). To reinforce these data, we generated QKI knockout (QKI$^{KO}$) cells using the CRISPR-Cas9 technology and performed infection assays (Figure 4B). We found that QKI$^{KO}$ cells were approximately 2-fold less permissive than parental cells, while trans-complementation with the human QKI cDNA rescued the infection (Figure 4C). Of note, our efforts to generate KO cells for the other selected HDFs (HNRPNA2B1, ILF3, SFPQ) were unsuccessful, suggesting that these genes are essential for the survival of A549 cells. To dissect the viral steps where these host factors are required, we infected A549-ACE2 cells knocked down for the expression of HNRPNA2B1, ILF3, QKI, and SFPQ, and quantified vRNA at different time points (Figure 4D). We did not observe any significant differences in vRNA levels in all of the conditions at 2 hpi, suggesting that these proteins are not involved in viral particle uptake. In contrast, we observed a large reduction in vRNA levels at 6 hpi, which increased at 24 hpi. Since in A549-ACE2, the beginning of vRNA amplification occurred between 4 and 6 hpi (data not shown), it is likely that HNRPNA2B1, ILF3, QKI, and

<footnote>correction. Host dependency factors are marked in blue and host restriction factors in are marked in red. Positive controls (CTSL1 and ATP6V1B2) are highlighted in yellow.</footnote>

<footnote>(D) Intersection of the data obtained from N protein quantification by flow cytometry and viruses released in the supernatant of infected cells by qRT-PCR. The data shown are the means of 2 biological replicates (with technical duplicates). Host dependency factors are marked in blue and host restriction factors are marked in red.</footnote>

<footnote>(E) A549-ACE2 were infected with SARS-CoV-2 (MOI 0.05) in the continuous presence of compounds (10 and 1 µM). Viruses released in supernatant were quantified 24 hpi by qRT-PCR (top panel). Cell viability was assessed in parallel (bottom panel). The data shown are the means ± SDs of 3 biological replicates (with technical duplicates). Significance was calculated using a 2-way ANOVA statistical test with Dunnett's multiple comparisons test (n.s, not significant; **p < 0.01, ***p < 0.001, ****p < 0.0001).</footnote>
<figure><img src="image_26.png" /><figcaption>B</figcaption></figure>

<figure><img src="image_27.png" /><figcaption>C</figcaption></figure>

<figure><img src="image_28.png" /><figcaption>D</figcaption></figure>

<figure><img src="image_29.png" /><figcaption>Figure 4. HNRNPA2B1, ILF3, QKI, and SFPQ mediate SARS-CoV-2 replication</figcaption></figure>

(A) A549-ACE2 cells were reverse transfected with the indicated siRNA pool. In 48 h, protein knockdown was ascertained by western blot (left panel) and cells were challenged with SARS-CoV-2 (MOI 0.01) for 24 h. Viral replication was assessed by flow cytometry using anti-N protein mAb and normalized to the siNT-transfected cells (center panel); the yield of viral particles released in the supernatant of infected cells was quantified by qRT-PCR (right panel). The data shown are the means ± SDs of 3 biological replicates (with technical triplicates). Significance was calculated using 1-way ANOVA statistical test with Dunnett's multiple comparisons test.

(legend continued on next page)
SFPQ are involved in vRNA replication. Nonetheless, because SARS-CoV-2 vRNA translation and replication are tight, intricate steps, we could not rule out a role of some of these proteins in vRNA translation.

In conclusion, our data provide a landscape of functional interactions that the SARS-CoV-2 RNA genome establishes with the host cell during infection. We provide robust evidence showing that SARS-CoV-2 interacts with several cellular RBPs to facilitate viral replication, some of which could represent promising targets for antiviral intervention. All of this information constitutes a valuable resource for future investigations into the mechanism of SARS-CoV-2 replication and the identification of antiviral drugs to manage the emergence of future variants of concern or new human coronavirus(es).

## Limitations of the study

This study provides a valuable resource to improve our understanding of the molecular interactions that SARS-CoV-2 establishes with the infected cells. However, our ChIRP-MS was performed in a cell line that is not derived from the lung; thus, the vRNA-host protein interactions shown in this study could not fully cover the whole interactions that take place during natural infection. It will be interesting to extend this study to more physiological cellular models, such as primary cells or cells derived from nasopharyngeal or lung epithelium, to derive an overall view of the SARS-CoV-2 interactome.

## STAR★METHODS

Detailed methods are provided in the online version of this paper and include the following:

* **KEY RESOURCES TABLE**

* **RESOURCE AVAILABILITY**

* Lead contact

* Materials availability

* Data and code availability

* **EXPERIMENTAL MODEL AND SUBJECT DETAILS**

* Cell lines

* **METHOD DETAILS**

* Virus preparation and titration

* Generation of cell lines stably expressing ACE2

* Flow cytometry analysis of ACE2 cell surface expression

* Comprehensive identification of RNA binding proteins by mass spectrometry

* High-confidence RNA binding protein scoring

* Identification of known interactors of coronaviruses in the ChIRP dataset

* siRNA screen assay

* Drug treatment and SARS-CoV-2 infection

* Viral infection quantification assays

* Kinetic of infection by qPCR assay

* Viral RNA immunoprecipitation (vRIP)

* Gene editing and trans-complementation experiments

* Cell viability assay

* **QUANTIFICATION AND STATISTICAL ANALYSIS**

* Statistical analyses

## SUPPLEMENTAL INFORMATION

Supplemental information can be found online at https://doi.org/10.1016/j.celrep.2022.110744.

## ACKNOWLEDGMENTS

This work was supported by the French government's Investissement d'Avenir program, Laboratoire d'Excellence "Integrative Biology of Emerging Infectious Diseases" (grant no. ANR-10-LABX-62-IBEID), the Fondation pour la Recherche Médicale (grant no. FRM-EQU202003010193), the Agence Nationale de la Recherche (ANR/FRM Flash COVID project IDISCOVR and ANR-20-CO11-0004 project FISHBP), and University Paris Cité (Plan de Soutien Covid-19: RACPL20FIR01-COVID-SOUL). The authors thank Alessia Zamborlini for critical reading of the manuscript. L.M. and A.A. dedicate this work to the memory of Professor Jean-Louis Virelizier (Unité d'Immunologie Virale, Institut Pasteur, Paris) and Professor Renaud Mahieux (Ecole Normale Supérieure, Lyon, France), who left us during the SARS-CoV-2 epidemic.

## AUTHOR CONTRIBUTIONS

L.M. and A.A. conceived of the project; L.M. supervised the project. L.M. performed the ChIRP experiments and analyzed the data with the help of M.P. and P.-O.V. A.L.-U., P.-O.V., and V.S. performed the bioinformatic analysis. L.B.-M. and L.M. produced the SARS-CoV-2 virus stocks used in this study. P.-O.V. and V.L. analyzed the viral protein-RBP protein interactions. A.L. and L.M. performed the siRNA screen and the vRIP assays. A.L.-U. screened the compounds. A.L.-U., L.F.-S., and L.M. performed the knockdown and KO experiments. L.M., P.-O.V., A.L.-U., and A.A. wrote the manuscript, with input of all of the authors. A.A. and L.M. secured the funding for the study.

## DECLARATION OF INTERESTS

The authors declare no competing interests.

Received: March 19, 2021

Revised: July 28, 2021

Accepted: April 7, 2022

Published: April 26, 2022

## REFERENCES

Banerjee, A.K., Blanco, M.R., Bruce, E.A., Honson, D.D., Chen, L.M., Chow, A., Bhat, P., Ollikainen, N., Quinodoz, S.A., Loney, C., et al. (2020). SARS-CoV-2 disrupts splicing, translation, and protein trafficking to suppress host defenses. Cell 183, 1325-1339.e21.

(B) Immunoblot with anti-FLAG or anti-QKI mAb in parental, QKI<sup>KO</sup> A549-ACE2 cells, and QKI<sup>KO</sup> cells trans-complemented with GFP or QKI cDNA.

(C) Indicated cells were inoculated with SARS-CoV-2 (MOI 0.01). Infection was quantified 24 h later by flow cytometry using anti-N protein mAb. The data shown are the means ± SDs of 3 biological replicates (with technical duplicates). Significance was calculated by the Kruskal-Wallis test with Dunnett's multiple comparisons test.

(D) A549-ACE2 cells were treated as in (A). At the indicated time points, cells were treated with trypsin to remove cell surface-bound viral particles and viral RNA was quantified by qRT-PCR. The data shown are the means ± SDs of 3 biological replicates (with technical triplicates). The p values were calculated by 1-way ANOVA with Dunnett's multiple comparisons test (n.s, non-significant; \*p < 0.05, \*\*p < 0.001, \*\*\*p < 0.0001).
Bell, J.L., Wächter, K., Mühleck, B., Pazaitis, N., Köhn, M., Lederer, M., and Hüttelmaier, S. (2013). Insulin-like growth factor 2 mRNA-binding proteins (IGF2BPs): post-transcriptional drivers of cancer progression? *Cell Mol. Life Sci.* 70, 2657–2675.

Bindea, G., Mlecnik, B., Hackl, H., Charoentong, P., Tosolini, M., Kirilovsky, A., Fridman, W.-H., Pagès, F., Trajanoski, Z., and Galon, J. (2009). ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. *Bioinformatics* 25, 1091–1093.

Cabrita, L.D., Cassaignau, A.M.E., Launay, H.M.M., Waudby, C.A., Wlodarski, T., Camilloni, C., Karyadi, M.-E., Robertson, A.L., Wang, X., Wentink, A.S., et al. (2016). A structural ensemble of a ribosome–nascent chain complex during cotranslational protein folding. *Nat. Struct. Mol. Biol.* 23, 278–285.

Cele, S., Gazy, I., Jackson, L., Hwa, S.-H., Tegally, H., Lustig, G., Giandhari, J., Pillay, S., Wilkinson, E., Naidoo, Y., et al. (2021). Escape of SARS-CoV-2 501Y.V2 from neutralization by convalescent plasma. *Nature* 593, 142–146.

Corman, V.M., Landt, O., Kaiser, M., Molenkamp, R., Meijer, A., Chu, D.K., Bleicker, T., Brünink, S., Schneider, J., Schmidt, M.L., et al. (2020). Detection of 2019 novel coronavirus (2019-nCoV) by real-time RT-PCR. *Eurosurveillance* 25, 2000045.

Cortese, M., Lee, J.-Y., Cerikan, B., Neufeldt, C.J., Oorschot, V.M.J., Köhrer, S., Hennies, J., Schieber, N.L., Ronchi, P., Mizzon, G., et al. (2020). Integrative imaging reveals SARS-CoV-2-induced reshaping of subcellular morphologies. *Cell Host Microbe* 28, 853–866.e5.

Daniloski, Z., Jordan, T.X., Wessels, H.-H., Hoagland, D.A., Kasela, S., Legut, M., Maniatis, S., Mimitou, E.P., Lu, L., Geller, E., et al. (2020). Identification of required host factors for SARS-CoV-2 infection in human cells. *Cell* 184, 92–105.e16.

Edgil, D., Polacek, C., and Harris, E. (2006). Dengue virus utilizes a novel strategy for translation initiation when cap-dependent translation is inhibited. *J. Virol.* 80, 2976–2986.

Eymieux, S., Rouillé, Y., Terrier, O., Seron, K., Blanchard, E., Rosa-Calatrava, M., Duboisson, J., Belouzard, S., and Roingeard, P. (2021). Ultrastructural modifications induced by SARS-CoV-2 in Vero cells: a kinetic analysis of viral factory formation, viral particle morphogenesis and virion release. *Cell. Mol. Life Sci.* 78, 3565–3576.

Filipowicz, W. (2014). Making ends meet: a role of RNA ligase RTCB in unfolded protein response. *EMBO J.* 33, 2887–2889.

Fislová, T., Thomas, B., Graef, K.M., and Fodor, E. (2010). Association of the influenza virus RNA polymerase subunit PB2 with the host chaperonin CCT. *J. Virol.* 84, 8691–8699.

Flynn, R.A., Belk, J.A., Qi, Y., Yasumoto, Y., Wei, J., Alfajaro, M.M., Shi, Q., Mumbach, M.R., Limaye, A., DeWeirdt, P.C., et al. (2021). Discovery and functional interrogation of SARS-CoV-2 RNA-host protein interactions. *Cell* 184, 2394–2411.e16.

Fung, T.S., and Liu, D.X. (2014). Coronavirus infection, ER stress, apoptosis and innate immunity. *Front. Microbiol.* 5, 296.

Garcia-Beltran, W.F., Lam, E.C., St. Denis, K., Nitido, A.D., Garcia, Z.H., Hauser, B.M., Feldman, J., Pavlovic, M.N., Gregory, D.J., Poznansky, M.C., et al. (2021). Multiple SARS-CoV-2 variants escape neutralization by vaccine-induced humoral immunity. *Cell* 184, 2372–2383.e9.

Garcia-Moreno, M., Noerenberg, M., Ni, S., Järvelin, A.I., González-Almela, E., Lenz, C.E., Bach-Pages, M., Cox, V., Avolio, R., Davis, T., et al. (2019). System-wide profiling of RNA-binding proteins uncovers key regulators of virus infection. *Mol. Cell* 74, 196–211.e11.

Gélinas, J.-F., Clerzius, G., Shaw, E., and Gatignol, A. (2011). Enhancement of replication of RNA viruses by ADAR1 via RNA editing and inhibition of RNA-activated protein kinase. *J. Virol.* 85, 8460–8466.

Geuens, T., Bouhy, D., and Timmerman, V. (2016). The hnRNP family: insights into their role in health and disease. *Hum. Genet.* 135, 851–867.

Gordon, D.E., Jang, G.M., Bouhaddou, M., Xu, J., Obernier, K., White, K.M., O'Meara, M.J., Rezelj, V.V., Guo, J.Z., Swaney, D.L., et al. (2020a). A SARS-CoV-2 protein interaction map reveals targets for drug repurposing. *Nature* 583, 459–468.

Gordon, D.E., Hiatt, J., Bouhaddou, M., Rezelj, V.V., Ulferts, S., Braberg, H., Jureka, A.S., Obernier, K., Guo, J.Z., Batra, J., et al. (2020b). Comparative host-coronavirus protein interaction networks reveal pan-viral disease mechanisms. *Science* 370, eabe9403.

Guo, X., Ma, J., Sun, J., and Gao, G. (2007). The zinc-finger antiviral protein recruits the RNA processing exosome to degrade the target mRNA. *Proc. Natl. Acad. Sci. U S A* 104, 151–156.

Hafirassou, M.L., Meertens, L., Umaña-Diaz, C., Labeau, A., Dejarnac, O., Bonnet-Madin, L., Kümmerer, B.M., Delaugerre, C., Roingeard, P., Vidalain, P.-O., et al. (2017). A global interactome map of the dengue virus NS1 identifies virus restriction and dependency host factors. *Cell Rep.* 21, 3900–3913.

Hoffmann, H.-H., Sánchez-Rivera, F.J., Schneider, W.M., Luna, J.M., Soto-Feliciano, Y.M., Ashbrook, A.W., Le Pen, J., Leal, A.A., Ricardo-Lax, I., Michailidis, E., et al. (2021). Functional interrogation of a SARS-CoV-2 host protein interactome identifies unique and shared coronavirus host factors. *Cell Host Microbe* 29, 267–280.e5.

Inoue, Y., Aizaki, H., Hara, H., Matsuda, M., Ando, T., Shimoji, T., Murakami, K., Masaki, T., Shoji, I., Homma, S., et al. (2011). Chaperonin TRiC/CCT participates in replication of hepatitis C virus genome via interaction with the viral NS5B protein. *Virology* 410, 38–47.

Jiang, H., Zhang, H., Meng, Q., Xie, J., Li, Y., Chen, H., Zheng, Y., Wang, X., Qi, H., Zhang, J., et al. (2020). SARS-CoV-2 Orf9b suppresses type I interferon responses by targeting TOM70. *Cell Mol. Immunol.* 17, 998–1000.

Kamel, W., Noerenberg, M., Cerikan, B., Chen, H., Järvelin, A.I., Kammoun, M., Lee, J.Y., Shuai, N., Garcia-Moreno, M., Andrejeva, A., et al. (2021). Global analysis of protein-RNA interactions in SARS-CoV-2-infected cells reveals key regulators of infection. *Mol. Cell* 81, 2851–2867.e7.

Kim, D., Lee, J.-Y., Yang, J.-S., Kim, J.W., Kim, V.N., and Chang, H. (2020). The architecture of SARS-CoV-2 transcriptome. *Cell* 181, 914–921.e10.

Landeras-Bueno, S., Jorba, N., Pérez-Cidoncha, M., and Ortín, J. (2011). The splicing factor proline-glutamine rich (SFPQ/PSF) is involved in influenza virus transcription. *PLoS Pathog.* 7, e1002397.

Lee, S., Lee, Y., Choi, Y., Son, A., Park, Y., Lee, K.-M., Kim, J., Kim, J.-S., and Kim, V.N. (2021). The SARS-CoV-2 RNA interactome. *Mol. Cell* 81, 2838–2850.e6.

Li, C., Ge, L., Li, P., Wang, Y., Sun, M., Huang, L., Ishag, H., Di, D., Shen, Z., Fan, W., et al. (2013). The DEAD-box RNA helicase DDX5 acts as a positive regulator of Japanese encephalitis virus replication by binding to viral 3' UTR. *Antivir. Res.* 100, 487–499.

Li, J., Guo, M., Tian, X., Wang, X., Yang, X., Wu, P., Liu, C., Xiao, Z., Qu, Y., Yin, Y., et al. (2020). Virus-host interactome and proteomic survey reveal potential virulence factors influencing SARS-CoV-2 pathogenesis. *Med* 12, 99–112.e7.

Liao, K.-C., Chuo, V., Fagg, W.S., Bradrick, S.S., Pompon, J., and Garcia-Blanco, M.A. (2020). The RNA binding protein Quaking represses host interferon response by downregulating MAVS. *RNA Biol.* 17, 366–380.

Liu, C., Ginn, H.M., Dejnirattisai, W., Supasa, P., Wang, B., Tuekprakhon, A., Nutalai, R., Zhou, D., Mentzer, A.J., Zhao, Y., et al. (2021a). Reduced neutralization of SARS-CoV-2 B.1.617 by vaccine and convalescent serum. *Cell* 184, 4220–4236.e13.

Liu, L., Iketani, S., Guo, Y., Chan, J.F.-W., Wang, M., Liu, L., Luo, Y., Chu, H., Huang, Y., Nair, M.S., et al. (2021b). Striking antibody evasion manifested by the omicron variant of SARS-CoV-2. *Nature* 602, 678–681.

Mangus, D.A., Evans, M.C., and Jacobson, A. (2003). Poly(A)-binding proteins: multifunctional scaffolds for the post-transcriptional control of gene expression. *Genome Biol.* 4, 223.

Marmisolle, F.E., García, M.L., and Reyes, C.A. (2018). RNA-binding protein immunoprecipitation as a tool to investigate plant miRNA processing interference by regulatory proteins of diverse origin. *Plant Methods* 14, 9.

Meertens, L., Hafirassou, M.L., Couderc, T., Bonnet-Madin, L., Kril, V., Kümmerer, B.M., Labeau, A., Brugier, A., Simon-Loriere, E., Burlaud-Gaillard, J., et al. (2019). FHL1 is a major host factor for chikungunya virus infection. *Nature* 574, 259–263.
Mo, Q., Xu, Z., Deng, F., Wang, H., and Ning, Y.-J. (2020). Host restriction of emerging high-pathogenic bunyaviruses via MOV10 by targeting viral nucleoprotein and blocking ribonucleoprotein assembly. PLoS Pathog. 16, e1009129.

Naganuma, T., Nakagawa, S., Tanigawa, A., Sasaki, Y.F., Goshima, N., and Hirose, T. (2012). Alternative 3'-end processing of long noncoding RNA initiates construction of nuclear paraspeckles. EMBO J. 31, 4020-4034.

Nakagawa, S., Yamazaki, T., and Hirose, T. (2018). Molecular dissection of nuclear paraspeckles: towards understanding the emerging world of the RNP milieu. Open Biol. 8, 180150.

Narayanan, K., Maeda, A., Maeda, J., and Makino, S. (2000). Characterization of the coronavirus M protein and nucleocapsid interaction in infected cells. J. Virol. 74, 8127-8134.

Oakland, T.E., Haselton, K.J., and Randall, G. (2013). EWSR1 binds the hepatitis C virus cis-acting replication element and is required for efficient viral replication. J. Virol. 87, 6625-6634.

Onodi, F., Bonnet-Madin, L., Meertens, L., Karpf, L., Poirot, J., Zhang, S.-Y., Picard, C., Puel, A., Jouanguy, E., Zhang, Q., et al. (2021). SARS-CoV-2 induces human plasmacytoid predendritic cell diversification via UNC93B and IRAK4. Preprint at bioRxiv. https://doi.org/10.1101/2020.07.10.197343.

Ooi, Y.S., Majzoub, K., Flynn, R.A., Mata, M.A., Diep, J., Li, J.K., van Buuren, N., Rumachik, N., Johnson, A.G., Puschnik, A.S., et al. (2019). An RNA-centric dissection of host complexes controlling flavivirus infection. Nat. Microbiol. 4, 2369-2382.

de la Parra, C., Ernlund, A., Alard, A., Ruggles, K., Ueberheide, B., and Schneider, R.J. (2018). A widespread alternate form of cap-dependent mRNA translation initiation. Nat. Commun. 9, 3068.

Perrin-Cocon, L., Diaz, O., Jacquemin, C., Barthel, V., Ogire, E., Ramière, C., André, P., Lotteau, V., and Vidalain, P.-O. (2020). The current landscape of coronavirus-host protein-protein interactions. J. Transl. Med. 18, 319.

Poland, G.A., Ovsyannikova, I.G., and Kennedy, R.B. (2020). SARS-CoV-2 immunity: review and applications to phase 3 vaccine candidates. Lancet 396, 1595-1606.

Reed, R., and Hurt, E. (2002). A conserved mRNA export machinery coupled to pre-mRNA splicing. Cell 108, 523-531.

Schmidt, N., Lareau, C.A., Keshishian, H., Ganskih, S., Schneider, C., Hennig, T., Melanson, R., Werner, S., Wei, Y., Zimmer, M., et al. (2020). The SARS-CoV-2 RNA-protein interactome in infected human cells. Nat. Microbiol. 6, 339-353.

Schneider, W.M., Luna, J.M., Hoffmann, H.-H., Sánchez-Rivera, F.J., Leal, A.A., Ashbrook, A.W., Le Pen, J., Ricardo-Lax, I., Michailidis, E., Peace, A., et al. (2020). Genome-scale identification of SARS-CoV-2 and pan-coronavirus host factor networks. Cell 184, 120-132.e14.

Schoeman, D., and Fielding, B.C. (2019). Coronavirus envelope protein: current knowledge. Virol. J. 16, 69.

Schubert, K., Karousis, E.D., Jomaa, A., Scaiola, A., Echeverria, B., Gurzeler, L.-A., Leibundgut, M., Thiel, V., Mühlemann, O., and Ban, N. (2020). SARS-CoV-2 Nsp1 binds the ribosomal mRNA channel to inhibit translation. Nat. Struct. Mol. Biol. 27, 959-966.

Shannon, P., Markiel, A., Ozier, O., Baliga, N.S., Wang, J.T., Ramage, D., Amin, N., Schwikowski, B., and Ideker, T. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498-2504.

Shepard, P.J., and Hertel, K.J. (2009). The SR protein family. Genome Biol. 10, 242.

Snijder, E.J., Decroly, E., and Ziebuhr, J. (2016). Chapter three - the nonstructural proteins directing coronavirus RNA synthesis and processing. In Advances in Virus Research, J. Ziebuhr, ed. (Academic Press), pp. 59-126.

Sun, L., Li, P., Ju, X., Rao, J., Huang, W., Ren, L., Zhang, S., Xiong, T., Xu, K., Zhou, X., et al. (2021). In vivo structural characterization of the SARS-CoV-2 RNA genome identifies host proteins vulnerable to repurposed drugs. Cell 184, 1865-1883.e20.

Teo, G., Liu, G., Zhang, J., Nesvizhskii, A.I., Gingras, A.-C., and Choi, H. (2014). SAINTexpress: improvements and additional features in Significance Analysis of INTERactome software. J. Proteomics 100, 37-43.

Thompson, M.G., Dittmar, M., Mallory, M.J., Bhat, P., Ferretti, M.B., Fontoura, B.M., Cherry, S., and Lynch, K.W. (2020). Viral-induced alternative splicing of host genes promotes influenza replication. Elife 9, e55500.

V'kovski, P., Gerber, M., Kelly, J., Pfaender, S., Ebert, N., Braga Lagache, S., Simillion, C., Portmann, J., Stalder, H., Gaschen, V., et al. (2019). Determination of host proteins composing the microenvironment of coronavirus replicase complexes by proximity-labeling. Elife 8, e42037.

V'kovski, P., Kratzel, A., Steiner, S., Stalder, H., and Thiel, V. (2021). Coronavirus biology and replication: implications for SARS-CoV-2. Nat. Rev. Microbiol. 19, 155-170.

Wang, L., Sola, I., Enjuanes, L., and Zuñiga, S. (2021a). MOV10 helicase interacts with coronavirus nucleocapsid protein and has antiviral activity. mBio 12, e0131621.

Wang, M., Cao, R., Zhang, L., Yang, X., Liu, J., Xu, M., Shi, Z., Hu, Z., Zhong, W., and Xiao, G. (2020). Remdesivir and chloroquine effectively inhibit the recently emerged novel coronavirus (2019-nCoV) in vitro. Cell Res. 30, 269-271.

Wang, R., Simoneau, C.R., Kulsuptrakul, J., Bouhaddou, M., Travisano, K.A., Hayashi, J.M., Carlson-Stevermer, J., Zengel, J.R., Richards, C.M., Fozouni, P., et al. (2021b). Genetic screens identify host factors for SARS-CoV-2 and common cold coronaviruses. Cell 184, 106-119.e14.

Watson, S.F., Bellora, N., and Macias, S. (2020). ILF3 contributes to the establishment of the antiviral type I interferon program. Nucleic Acids Res. 48, 116-129.

Wei, J., Alfajaro, M.M., DeWeirdt, P.C., Hanna, R.E., Lu-Culligan, W.J., Cai, W.L., Strine, M.S., Zhang, S.-M., Graziano, V.R., Schmitz, C.O., et al. (2020). Genome-wide CRISPR screens reveal host factors critical for SARS-CoV-2 infection. Cell 184, 76-91.e13.

de Wilde, A.H., Snijder, E.J., Kikkert, M., and van Hemert, M.J. (2018). Host factors in coronavirus replication. Curr. Top Microbiol. Immunol. 419, 1-42.

Yang, P., Mathieu, C., Kolaitis, R.-M., Zhang, P., Messing, J., Yurtsever, U., Yang, Z., Wu, J., Li, Y., Pan, Q., et al. (2020). G3BP1 is a tunable switch that triggers phase separation to assemble stress granules. Cell 181, 325-345.e28.

Zhou, B., Wu, F., Han, J., Qi, F., Ni, T., and Qian, F. (2019a). Exploitation of nuclear protein SFPQ by the encephalomyocarditis virus to facilitate its replication. Biochem. Biophys. Res. Commun. 510, 65-71.

Zhou, D., Dejnirattisai, W., Supasa, P., Liu, C., Mentzer, A.J., Ginn, H.M., Zhao, Y., Duyvesteyn, H.M.E., Tuekprakhon, A., Nutlai, R., et al. (2021). Evidence of escape of SARS-CoV-2 variant B.1.351 from natural and vaccine-induced sera. Cell 184, 2348-2361.e6.

Zhou, P., Yang, X.-L., Wang, X.-G., Hu, B., Zhang, L., Zhang, W., Si, H.-R., Zhu, Y., Li, B., Huang, C.-L., et al. (2020). A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 579, 270-273.

Zhou, Y., Zhou, B., Pache, L., Chang, M., Khodabakhshi, A.H., Tanaseichuk, O., Benner, C., and Chanda, S.K. (2019b). Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat. Commun. 10, 1523.

Zhu, N., Zhang, D., Wang, W., Li, X., Yang, B., Song, J., Zhao, X., Huang, B., Shi, W., Lu, R., et al. (2020). A novel coronavirus from patients with pneumonia in China, 2019. N. Engl. J. Med. 382, 727-733.
## STAR★METHODS

### KEY RESOURCES TABLE

<table><thead><tr><th>REAGENT or RESOURCE</th><th>SOURCE</th><th>IDENTIFIER</th></tr></thead><tbody><tr><td colspan="3"><strong>Antibodies</strong></td></tr><tr><td>Mouse monoclonal anti-Human/Rat/Hamster ACE2</td><td>Biotechne</td><td>Cat #MAB9332</td></tr><tr><td>Mouse monoclonal SARS-CoV/SARS-CoV-2 Nucleocapsid</td><td>Sino Biological</td><td>Cat #40143-MM05; RRID:AB_2827977</td></tr><tr><td>Alexa Fluor® 647 AffiniPure Donkey Anti-Mouse IgG (H + L)</td><td>Jackson ImmunoResearch</td><td>Cat #715-605-150; RRID:AB_2340862</td></tr><tr><td>Monoclonal ANTI-FLAG M2 antibody</td><td>Sigma Aldrich</td><td>Cat #F1804; RRID:AB_262044</td></tr><tr><td>Rabbit monoclonal HA-Tag (C29F4)</td><td>Cell Signaling Technology</td><td>Cat #3724S; RRID:AB_1549585</td></tr><tr><td>Rabbit polyclonal GAPDH (0411) Antibody</td><td>SantaCruz Biotechnology</td><td>Cat #sc-47724; RRID:AB_627678</td></tr><tr><td>Peroxidase AffiniPure Donkey anti-Rabbit IgG</td><td>Jackson ImmunoResearch</td><td>Cat #711-035-152; RRID:AB_10015282</td></tr><tr><td>Rabbit anti-mouse Immunoglobulins/HRP</td><td>Dako</td><td>Cat #P0260; AB_2636929</td></tr><tr><td>Rabbit monoclonal anti-ILF3</td><td>Abcam</td><td>Cat #Ab 92355; RRID:AB_2049804</td></tr><tr><td>Rabbit monoclonal anti-QKI</td><td>Abcam</td><td>Cat #Ab 126742; RRID:AB_11129508</td></tr><tr><td>Rabbit monoclonal anti-HNRNPA2B1</td><td>Abcam</td><td>Cat #Ab 259894</td></tr><tr><td>Rabbit monoclonal anti-SFPQ</td><td>Abcam</td><td>Cat #Ab 177149</td></tr><tr><td colspan="3"><strong>Bacterial and virus strains</strong></td></tr><tr><td>SARS-CoV-2 isolate 220_95</td><td>Isolated from patient<br>(Onodi et al., 2021)</td><td>EPI_ISL_469284</td></tr><tr><td colspan="3"><strong>Chemicals, peptides, and recombinant proteins</strong></td></tr><tr><td>32% Paraformaldehyde (formaldehyde) aqueous solution</td><td>Electron Microscopy Sciences</td><td>Cat #15714</td></tr><tr><td>DMSO</td><td>Sigma-Aldrich</td><td>Cat #D5879</td></tr><tr><td>Triton X-100</td><td>Sigma-Aldrich</td><td>Cat #T9284</td></tr><tr><td>Artemisinin</td><td>Sigma-Aldrich</td><td>Cat #361593</td></tr><tr><td>(S)-(+)-Camphotecin</td><td>Sigma-Aldrich</td><td>Cat #C9911</td></tr><tr><td>Cerulenin</td><td>Sigma-Aldrich</td><td>Cat #C2389</td></tr><tr><td>Dabigatran</td><td>Sigma-Aldrich</td><td>Cat #SML2370</td></tr><tr><td>Geldanamycin from Streptomyces hygroscopicus</td><td>Sigma-Aldrich</td><td>Cat #G3381</td></tr><tr><td>Phenethyl isothiocyanate</td><td>Sigma-Aldrich</td><td>Cat #253731</td></tr><tr><td>Remdesivir</td><td>Selleckchem</td><td>Cat #SE-S8932</td></tr><tr><td>Sunitinib malate</td><td>Sigma-Aldrich</td><td>Cat #PZ0012</td></tr><tr><td>MyOne C1 beads</td><td>Thermo Fisher Scientific</td><td>Cat #65001</td></tr><tr><td>50 mM D-biotin</td><td>Thermo Fisher Scientific</td><td>Cat #B20656</td></tr><tr><td>SUPERase-in</td><td>Thermo Fisher Scientific</td><td>Cat #AM2696</td></tr><tr><td>NA-Deoxycholate</td><td>Sigma Merck</td><td>Cat #30970</td></tr><tr><td>N-Laurosylsarcosine sodium salt</td><td>Sigma Merck</td><td>Cat #L7414</td></tr><tr><td>EDTA ULTROL grade</td><td>Sigma Merck</td><td>Cat #324504</td></tr><tr><td>Trichloroacetic acid</td><td>Sigma Merck</td><td>Cat #T0699</td></tr><tr><td>Proteinase K</td><td>Qiagen</td><td>Cat #19133</td></tr><tr><td>TRizol LS Reagent</td><td>Thermo Fisher Scientific</td><td>Cat #10296028</td></tr><tr><td>Dynabeads Protein G</td><td>Thermo Fisher Scientific</td><td>Cat #10004D</td></tr><tr><td>3XFLAG peptide</td><td>Sigma Merck</td><td>Cat #F4799</td></tr><tr><td>HA peptide</td><td>Sigma Merck</td><td>Cat #11666975001</td></tr><tr><td>GlycoBlue Coprecipitant</td><td>Thermo Fisher Scientific</td><td>Cat #AM9515</td></tr><tr><td>Lipofectamine RNAiMax Transfection Reagent</td><td>Thermo Fisher Scientific</td><td>Cat #13778150</td></tr><tr><td>Lipofectamine 3000 Transfection Reagent</td><td>Thermo Fisher Scientific</td><td>Cat #L3000015</td></tr><tr><td>Bolt 4-12% Bis-Tris plus</td><td>Thermo Fisher Scientific</td><td>Cat #NW04122BOX</td></tr></tbody></table>

(Continued on next page)
Continued

<table><thead><tr><td>REAGENT or RESOURCE</td><td>SOURCE</td><td>IDENTIFIER</td></tr></thead><tbody><tr><td>Bolt 10% Bis-Tris plus</td><td>Thermo Fisher Scientific</td><td>Cat #NW00100BOX</td></tr><tr><td>2X Laemmli Sample Buffer</td><td>Biorad</td><td>Cat #1610737</td></tr></tbody></table>

Critical commercial assays

<table><tr><td>CellTiter-Glo® 2.0 Assay</td><td>Promega</td><td>Cat #G9242</td></tr><tr><td>Profection® Mammalian Transfection System</td><td>Promega</td><td>Cat #E1200</td></tr><tr><td>Luna® Universal One-Step RT-qPCR Kit</td><td>New England Biolabs</td><td>Cat #E3005E</td></tr><tr><td>Maxima First Strand cDNA Synthesis Kit</td><td>Thermo Fisher Scientific</td><td>Cat #K1671</td></tr></table>

Experimental models: Cell lines

<table><tr><td>Lenti-X 293T</td><td>Takara Bio</td><td>Cat #632180</td></tr><tr><td>Lenti-X 293T-ACE2</td><td>This manuscript</td><td>N/A</td></tr><tr><td>A549</td><td>ATCC</td><td>Cat #CCL-185</td></tr><tr><td>A549-ACE2</td><td>This manuscript</td><td>N/A</td></tr><tr><td>Vero E6</td><td>ATCC</td><td>Cat #CRL-1586</td></tr></table>

Oligonucleotides

<table><tr><td>ChIRP tilling oligos, see Table S1</td><td>This paper</td><td>N/A</td></tr><tr><td>E_Sarbeco Forward: ACA GGT ACG TTA</td><td>Corman et al., (2020)</td><td>N/A</td></tr><tr><td>ATA GTT AAT AGC GT</td><td></td><td></td></tr><tr><td>E_Sarbeco Reverse: ATA TTG CAG CAG</td><td>Corman et al., (2020)</td><td>N/A</td></tr><tr><td>TAC GCA CAC A</td><td></td><td></td></tr><tr><td>EWSR1 HA Forward: TAA ACT ATG AAT TCA</td><td>This paper</td><td>N/A</td></tr><tr><td>GGA TGG CGT CCA CGG ATT ACA GTA CCT</td><td></td><td></td></tr><tr><td>ATA GCC AAG CTG CAG CGC AGC AGG GC</td><td></td><td></td></tr><tr><td>EWSR1 HA Reverse: ATA GTT TAG CGG CCG</td><td>This paper</td><td>N/A</td></tr><tr><td>CTC AAG CGT AAT CTG GAA CAT CGT ATG</td><td></td><td></td></tr><tr><td>GGT ATC CTC CAG CGG CCG AGT AGG GCC</td><td></td><td></td></tr><tr><td>GAT CTC TGC GC</td><td></td><td></td></tr><tr><td>ACE2 Forward: CCG CTC GAG CGG GCC</td><td>This paper</td><td>N/A</td></tr><tr><td>ATG TCA AGC TCT TCC TGG CTC</td><td></td><td></td></tr><tr><td>ACE2 Reverse: ATA AGA ATG CGG CCG CTT</td><td>This paper</td><td>N/A</td></tr><tr><td>AAA AGG AGG TCT GAA CAT CAT CAG TG</td><td></td><td></td></tr><tr><td>ON-TARGETplus Human ILF3 siRNA - SMARTpool</td><td>Horizon Discovery</td><td>Cat #L-012442-00</td></tr><tr><td>ON-TARGETplus Human HNRNPA2B1 siRNA -SMARTpool</td><td>Horizon Discovery</td><td>Cat #L-011690-01</td></tr><tr><td>ON-TARGETplus Human QKI siRNA - SMARTpool</td><td>Horizon Discovery</td><td>Cat #L-024905-01</td></tr><tr><td>ON-TARGETplus Human SFPQ siRNA - SMARTpool</td><td>Horizon Discovery</td><td>Cat #L-006455-00</td></tr><tr><td>ON-TARGETplus Human CTSL siRNA SMART Pool</td><td>Horizon Discovery</td><td>Cat #L-005841-00</td></tr><tr><td>ON-TARGETplus Non-targeting Pool</td><td>Horizon Discovery</td><td>Cat #D-001810-10-05</td></tr><tr><td>QKI sgRNA targeting sequence:</td><td>This paper</td><td>N/A</td></tr><tr><td>GGATCTTCAACCACCTCGAG</td><td></td><td></td></tr></table>

Recombinant DNA

<table><tr><td>pLVX-EF1alpha-IRES-Puro</td><td>Takara Bio</td><td>Cat #631988</td></tr><tr><td>pLVX-ACE2</td><td>This manuscript</td><td>N/A</td></tr><tr><td>MGC Human EWSR1 Sequence-Verified cDNA</td><td>Horizon Discovery</td><td>Cat #MHS6278-202759357</td></tr><tr><td>3xFlag SFPQ_pUC57-mini Plasmid</td><td>Genscript</td><td>N/A</td></tr><tr><td>3xFlag QKI_pUC57 Plasmid</td><td>Genscript</td><td>N/A</td></tr><tr><td>3xFlag ILF3_pUC57-Mini Plasmid</td><td>Genscript</td><td>N/A</td></tr><tr><td>3xFlag HNRNPA2B1_pUC57 Plasmid</td><td>Genscript</td><td>N/A</td></tr><tr><td>3xFlag HSPB1_pUC57 Plasmid</td><td>Genscript</td><td>N/A</td></tr><tr><td>3xFlag GFP_pUC57 Plasmid</td><td>Genscript</td><td>N/A</td></tr><tr><td>3xFlag TRIM28_pUC57-mini Plasmid</td><td>Genscript</td><td>N/A</td></tr></table>

(Continued on next page)
Continued

<table><thead><tr><th>REAGENT or RESOURCE</th><th>SOURCE</th><th>IDENTIFIER</th></tr></thead><tbody><tr><td colspan="3"><strong>Software and algorithms</strong></td></tr><tr><td>Attune NxT</td><td>Thermo Fisher Scientific</td><td>Version 2.6</td></tr><tr><td>FlowJo</td><td>Tree Star</td><td>Version 10.0.8</td></tr><tr><td>Prism</td><td>GraphPad Software</td><td>Version 7.0d</td></tr><tr><td>Serial Cloner</td><td>SerialBasics</td><td>Version 2.6.1</td></tr><tr><td>ClueGo</td><td>Cytoscape</td><td>Version 2.5.7</td></tr><tr><td>Genamania</td><td>Cytoscape</td><td>Version 3.5.2</td></tr><tr><td colspan="3"><strong>Other</strong></td></tr><tr><td>Attune NxT</td><td>Thermo Fisher Scientific</td><td>N/A</td></tr><tr><td>TriStar<sup>2</sup> LB 942</td><td>BERTHOLD Technologies</td><td>N/A</td></tr><tr><td>LightCycler 480</td><td>Roche</td><td>N/A</td></tr></tbody></table>

RESOURCE AVAILABILITY

Lead contact

Further information and requests or resources and reagents should be directed to and will be fulfilled by the lead contact, Laurent Meertens (laurent.meertens@inserm.fr).

Materials availability

All the plasmids generated in this study, the 293T-ACE2 and A549-ACE2 are available on request.

Data and code availability

* The datasets of the mass spectrometry analyses are provided in supplemental files (Table S2). Additional Supplemental Items are available from Mendeley Data at https://data.mendeley.com/datasets/9pmrxdxdpyg/1.

* This paper does not report original code.

* Any additional information required to reanalyze the data reported in this paper is available from the Lead contact upon request.

EXPERIMENTAL MODEL AND SUBJECT DETAILS

Cell lines

Lenti-X 293T cells (human embryonic kidney cells), A549 cells (Human alveolar basal epithelial carcinoma), and Vero E6 cells (African green monkey kidney cells) were maintained in Dulbecco Modified Eagle Medium (DMEM; Invitrogen Life Technologies) supplemented with 10% heat-inactivated fetal bovine serum (FBS) and 1% penicillin/streptomycin (P/S) and 1% GlutaMAX (Life Technologies). All cell lines were periodically tested negative for mycoplasma contamination prior to use in experiments.

METHOD DETAILS

Virus preparation and titration

SARS-CoV-2 strain 220_95 (EPI_ISL_469284) was isolated from nasopharyngeal swab specimens collected at Service de Virologie (Hospital Saint Louis, Paris) (Onodi et al., 2021). Virus was propagated on Vero E6 in DMEM-2% (DMEM supplemented with 2% FBS, 1% P/S, 1% GlutaMAX, and 25 mM Hepes). Viruses were passaged three times before being used for experiments. For the last passage, viruses were purified through a 20% sucrose cushion by ultracentrifugation at 80,000 × g for 2 h at 4°C. Pellets were resuspended in TNE1X pH7.4 (Tris 50 mM, NaCl 100mM, EDTA 0.5 mM), aliquoted and stored at -80°C.

Viruses titer was determined by plaque assays in Vero E6 cells and expressed as PFU/mL. Cells were incubated for one hour with 10-fold serial dilution of viral stocks. The inoculum was then replaced with Avicel 2.4% mixed at equal volume with DMEM supplemented with 4% FBS, 2% Glutamax, 50mM MgCl₂, 0.225% of NaHCO₃, and incubated 3 days at 37°C. Then, Vero E6 cells were washed twice with phosphate-buffered saline (PBS) and stained with an aqueous solution containing 1% formaldehyde, 0.5% crystal violet, and 0.89% sodium chloride, for 15 min at room temperature. Vero E6 cells were then washed with water to visualize lysis plaques.

Generation of cell lines stably expressing ACE2

The human ACE2 coding DNA sequence (CDS) was amplified by RT-PCR from pLenti-ACE2 plasmid (a gift from Olivier Schwartz, Virus and Immunity, Institut Pasteur) using the ACE2 forward and reverse primers. Amplification product was cloned into
pLVX-IRES-ZsGreen1 vector using Xhol/NotI restriction enzymes. The ACE2 protein-expressing lentiviral vectors pseudotyped with vesicular stomatitis virus glycoprotein G (VSV-G) were generated by transfecting lenti-X 293T (5 × 10⁶) cells with pLVX-ACE2-ZsGreen, psPAX2, and pVSV-G (4:3:1 ratio) using calcium phosphate transfection (Promega) following the manufacturer's instructions. Supernatants were harvested 48 h after transfection, cleared by centrifugation, and filtered. Lentiviral vectors were purified as described above (see 'virus preparation and titration' section). 293T and A549 cells (2 × 10⁵) were seeded in 6 well plate and transduced by spinofection at 1,000 × g for 2 h at 33°C. Cell populations were sorted based on GFP expression with a BD FACSAria II (Becton Dickinson) with FACSDiva 6.1.2 software (Becton Dickinson).

## Flow cytometry analysis of ACE2 cell surface expression

Cells were detached with 2.5 mM EDTA in PBS and incubated with anti-human ACE2 mAb (5 mg/mL) in 100 mL of PBS with 0.02% NaN₃ and 5% FBS for 1 h at 4°C. Following the primary staining, cells were washed and stained with incubated with the Alexa 647-conjugated goat anti-mouse IgG (Jackson ImmunoResearch) for 30 min at 4°C. Acquisition was performed on an Attune NxT Flow Cytometer (Thermo Fisher Scientific) and analysis was done by using FlowJo software (Tree Star).

## Comprehensive identification of RNA binding proteins by mass spectrometry

ChIRP-M/S was performed as previously described by Ooi and colleagues (Ooi et al., 2019). SARS-CoV-2 probes were designed using the online tool (https://www.biosearchtech.com/stellaris-designer), with a repeat masking setting of 5 and even coverage of the whole viral genome (1 oligos every 200 bp). Oligos were synthesized with 3' biotin-TEG modification at Eurofins Genomics. The 129 oligos were resuspended at 100 µM and mixed at equal volume to prepare the probes pool.

293T cells stably expressing ACE2 were seeded at 9 × 10⁶ cells per T150 flask (3 flasks per condition) and inoculated with SARS-CoV-2 at a multiplicity of infection (MOI) of 0.001 or mock-treated, aiming for 50–60% of infected cells as monitored by flow cytometry (see "viral infection quantification assays" section) using a fraction of cells collected at 48 hpi. The media was removed 48 h later and the cells were rinsed once with PBS, trypsinized and pooled for each condition. Cells were centrifugated at 290 × g for 5 min and pellets were washed twice with PBS. The cells were resuspended with 35 mL of PBS containing 3% formaldehyde and incubated under agitation for 30 min at room temperature (RT). Of note, all the buffers used throughout the ChIRP experiments were freshly prepared. Cross-linking was terminated by adding glycine to a final concentration of 0.125 mM for 5 min at RT. Cross-linked cells were centrifugated at 1,000 × g for 5 min, washed once with PBS, weighted, then snap frozen and kept at –80°C until used. Cells were resuspended in 1 mL of lysis buffer (50 mM Tris-HCl, pH7.0, 10 mM EDTA, 1% SDS supplemented with SUPERase-IN and protease inhibitor) per 100 mg of cell pellet. The lysates were sonicated using a Diagenode Bioruptor. At this stage, 10 µL of lysates were kept as sample inputs. Lysates were precleared by adding 30 µL of washed MyOne C1 beads per ml of lysate and rocked for 30 min at 37°C. Beads were removed from the lysates using a magnetic stand and 2 mL of ChIRP hybridization buffer (750 mM NaCl, 1% SDS, 50 mM Tris-HCl pH 7.0, 1 mM EDTA, 15% formamide supplemented with 1 mM PMSF, protease inhibitor and SUPERase-in) and 1.5 µL of probes pool were added for every ml of lysate. After mixing, hybridization took place for 16h at 37°C under rotation. Then, 100 µL of washed MyOne C1 beads per microliter of probes pool were added to each sample and incubated for 45 min at 37°C. Magnetic beads were collected on a magnetic stand and washed 5 times with 1 mL of washing buffer (2x SSC, 0.5% SDS, add 1 mM PMSF) for 5 min at 37°C. To eluate the enriched ribonucleoprotein complexes, the beads were incubated with 300 µL of biotin elution buffer (12.5 mM D-biotin, 7.5mM Hepes pH 7.5, 75 mM NaCl, 1.5 mM EDTA, 0.15% SDS, 0.075% sarkosyl and 0.02% Na-Doexycholate) for 20 min at RT under rotation and then 15 min at 65°C with shaking. Eluate was collected in a fresh tube and the beads subject to elution again. The two eluents were pooled (600 µL) and 150 µL of trichloroacetic acid (25% of the total eluate volume) was added, mixed vigorously, and let overnight at 4°C for precipitation. The next day, the samples were centrifugated at 20,000 × g for 45 min at 4°C. Pellets were washed once with ice-cold acetone and spun at 20,000 × g for 5 min at 4°C. The acetone supernatant was removed, the tubes centrifuged briefly to carefully remove the residual acetone and the pellets were left to air-dry on the bench top. The pellets were solubilized in 30 µL of 2X Laemmli buffer with 20 mM DTT and boiled at 95°C for 30 min with a periodic mixing (every 5 min) to reverse the cross-linking.

One-tenth of the protein samples and their corresponding inputs (pre-diluted 1/10) were size-separated on a bolt 4–12% Bis-Tris plus (Thermo Fisher Scientific) and stained with the silver stain kit (Thermo Fisher Scientific) as per manufacturer's instructions. The remaining samples were size-separated on a bolt 10% Bis-Tris plus and stained with the Imperial protein stain (Thermo Fisher Scientific) according to manufacturer's instructions. For each of the 5 ChIRP experiments one gel slice from the SARS-CoV-2 and mock-infected samples were cut from the SDS-page. Gel slices from 5 independent replicates were sent for mass spectrometry analysis at Taplin Biological Mass Spectrometry Facility (Harvard Medical School) as previously described (Hafirassou et al., 2017).

## High-confidence RNA binding protein scoring

ChIRP-M/S data were analyzed using the online server (https://reprint-apms.org/) to run SAINTexpress (Teo et al., 2014). Total peptides count from each of the five technical replicates, including both mock- and SARS-CoV-2-infected cells, were used as quantitative values to run the analysis. To determine a core high-confidence interactome, we applied stringent filtering criteria with BFDR ≤ 0.05 and a fold change ≥ 5. We also defined an expanded interactome by loosen the BFDR requirement (BFDR < 0.1). Data shown in Figure 1E correspond to the fold change scores (FC) calculated by SAINTexpress from the average mean of spectral counts.
## Identification of known interactors of coronaviruses in the ChIRP dataset

A recently published meta-analysis by Perrin-Cocon et al. was used as a reference for known interactions between coronavirus and host proteins (Perrin-Cocon et al., 2020). In this report, the literature was manually curated to assemble a dataset of 1,311 coronavirus-host interactions encompassing 1,140 distinct cellular proteins. Host proteins in the ChIRP dataset overlapping with those present in this reference interactome were considered as previously reported or « known » interactors of coronaviruses. They were displayed together with the interacting-coronavirus proteins as an interaction network using Cytoscape (Shannon et al., 2003).

### siRNA screen assay

An arrayed ON-TARGETplus SMARTpool siRNA library targeting 138 of the 146 RNA-interacting proteins significantly enriched in our ChIRP-M/S was purchased from Horizon Discovery. SiRNA targeting 8 genes were not assessed because either no predesigned siRNA pools were available or proteins were not considered as candidates in our first analysis. These proteins are: hnRNP C-like 2, histone H2B type1-H, histone H3.2, putative 40S ribosomal protein S10-like, HSPA2, EEF1D, TKT and ACTN4. A549 cells stably expressing ACE2 were seeded in 24-well plate format (viral infection) or a 96-well plate (viability) and then were reverse transfected in duplicate with 30 nM of siRNA using the Lipofectamine RNAiMAX reagent (Invitrogen) according to manufacturer's instruction. 48 h post-transfection, cells seeded in a 24-well format were challenged with SARS-CoV-2 at a MOI of 0.05 in DMEM containing 5% FBS. After three hours at 37°C, virus inoculum was removed and replace with fresh medium. 24 h post-infection, cells and 200 µL of supernatants were collected for N protein staining and viral RNA quantification as described in "Viral infection quantification assays". Throughout the screening assay, siRNA-transfected cells seeded in 96-well format were incubated for 72 h at 37°C and used to assess cell viability.

### Drug treatment and SARS-CoV-2 infection

Using DGIdb, DrugBank, ChEMBL, and PanDrugs databases, we identified 274 compounds targeting 30 RBPs. Compounds were first filtered to select molecules described either as an inhibitor or with the highest interaction score. Finally, we selected the compounds targeting RBPs with a log2 (titer ratio) ≤ 0.5 in our siRNA screen.

Cells grown in 24-well plates were treated with the indicated compound concentrations for two hours prior infection and compounds were maintained throughout the course of infection. Cells treated with equivalent concentration of DMSO were used as control. The cells were challenged with SARS-CoV-2 at a MOI of 0.05 for 3 h, then washed once with PBS and incubated with fresh media. 24 h post-infection, cells and 200 µL of supernatants were collected for N protein staining and viral RNA quantification, as described in "Viral infection quantification assays". In parallel cells seeded in 96-well format were incubated 24 h with the drugs to assess cell viability.

### Viral infection quantification assays

For infection quantification by flow cytometry analysis, 24 h post-infection cells were trypsinized and fixed with 2% (v/v) paraformaldehyde (PFA) diluted in PBS for 15 min at room temperature. Cells were incubated for 45 min at 4°C with 1 µg/mL of the anti-N mAb diluted in permeabilization flow cytometry buffer (PBS supplemented with 5% FBS, 0.5% (w/v) saponin, 0.1% sodium azide). After washing, cells were incubated with 1 µg/mL of Alexa Fluor 647-conjugated goat anti-mouse IgG diluted in permeabilization flow cytometry buffer for 30 min at 4°C. Acquisition was performed on an Attune Nxt Flow Cytometer (Thermo Fisher Scientific) and analysis was done by using FlowJo software (Tree Star). To assess infectious viral particle release during infection by RT-qPCR, viruses were first inactivated by incubating the supernatants v/v with 1% Triton X-100 (Sigma) in PBS for 30 min under agitation at RT. Yields of viral RNA were quantified by real-time qPCR by using SARS-CoV-2 specific primers targeting the E gene with the Luna® Universal One-Step RT-qPCR Kit (New England Biolabs) in a LightCycler 480 thermocycler (Roche) according to the manufacturer's protocol. The number of viral genomes is expressed as PFU equivalent/mL and was calculated by performing a standard curve with a similarly treated supernatant from a viral stock with a known titer as described by Gordon et al. (Gordon et al., 2020b).

### Kinetic of infection by qPCR assay

A549-ACE2 cells, plated on 6 well plates and reverse transfected with the indicated siRNA, were inoculated with SARS-CoV-2 (MOI of 0.01). At indicated time point cells were washed twice with PBS, incubated with trypsin 0.25% for 5 min at 37°C to remove cells surface bound particles, and total RNA was extracted using the RNeasy plus mini kit (Qiagen) according to manufacturer's instruction. cDNAs were generated from 400 ng total RNA by using the Maxima First Strand Synthesis Kit following manufacturer's instruction (Thermo Fisher Scientific). Amplification products were incubated with 1 Unit of RNAse H for 20 min at 37°C, followed by 10 min at 72°C for enzyme inactivation, and diluted 10-fold in DNAse/RNAse free water. Real time quantitative PCR was performed using a Power Syber green PCR master Mix (Fisher Thermo Scientific) on a Light Cycler 480 (Roche). The primers used for qPCR were the E_Sarbeco primers for viral RNA quantification (Corman et al., 2020), and Quantitect primers for GAPDH were purchased from Qiagen. The relative expression quantification was performed based on the comparative threshold cycle (C~T~) method, using GAPDH as endogenous reference control.

### Viral RNA immunoprecipitation (v RIP)

The human EWSR1 coding DNA sequence (CDS) was amplified by PCR from MGC Human EWSR1 cDNA plasmid (Horizon Discovery) using the EWSR1-HA forward and reverse primers. Amplification product was cloned into pLVX-EF1alpha-IRES-Puro vector
using EcoRI/Notl restriction enzymes. Plasmids expressing FLAG-tagged version of HNRNPA2B1, ILF3, QKI, SFPQ, HSPB1, GFP and TRIM28 cloned in pUC57 vector were purchased from GenScript, and subcloned into pLVX-EF1alpha-IRES-Puro vector using EcoRI/Notl restriction enzymes, with the exception of SFPQ (EcoRI/BamHI). For A549-ACE2 transductions, lentiviral particles production and spinfection, were performed as described above for ACE2-expressing cells generation.

293T cells stably expressing ACE2 were seeded in 10 cm² dishes and transfected with 20 µg of plasmid using the Lipofectamine 3000 reagent (Thermo Fisher Scientific) according to manufacturer's instruction. 24 h post-transfection, cells were challenged with SARS-CoV-2 at a MOI of 0.1. Thirty hours post-infection, cells were lysed for 30 min at 4°C in IP lysis buffer (Thermo Fisher Scientific) supplemented with 1% Halt Protease and Phosphatase Inhibitor Cocktail (Thermo Fisher Scientific) and 0.5% SUPERase-in, then cleared by centrifugation for 30 min at 13,000 × g. For each condition, 1 mg of total proteins was transferred in a new tube in a 900 µL final volume for IP, from which 2 × 20 µL (2%) were kept as input for RNA extraction and immunoblot. Protein G magnetic beads (30 µL per conditions) were washed twice with 1 mL of IP Lysis buffer supplemented with 1% Triton and 0.01% SDS. Beads were resuspended in IP lysis buffer (100 µL per conditions) and incubated with anti-HA or anti-FLAG antibodies (3 µg per conditions) for one hour at room temperature under rotation. After 5 washes with IP lysis buffer, the beads were resuspended in IP lysis buffer (50 µL per conditions) supplemented with 1% Halt™ Protease and Phosphatase Inhibitor Cocktail and 0.5% SUPERase-in. One mg of total proteins was incubated 3 h at 4°C with 50 µL with the magnetic beads. Beads were washed five times with RIP buffer (KCl 150 mM, Tris-HCl 25 mM, EDTA 5mM, IPGAL 0.5%, 1% protease and phosphatase inhibitor and 0.5% SUPERase-in). Proteins were eluted twice (2 × 50 µL) with either HA (400 mg/mL) or 3XFLAG (200 mg/mL) peptides for 30 min at room temperature. Twenty microliters of the eluate were kept for the immunoblot, whilst the remaining eluate and 20 µL of the input were treated with proteinase K for 30 min at 50°C. Immunoprecipitated RNA was purified using Trizol (Qiagen) following the manufacturers recommendations.

For immunoblot, 20 µL of the inputs and eluates were resuspended in LDS Sample Buffer 4X (Thermo Fisher Scientific) containing 25 mM dithiothreitol (DTT) and were separated by SDS PAGE. Immunoblotting was performed as previously described (Hafirassou et al., 2017) with an anti-HA rabbit monoclonal antibody (Cell Signaling) or an anti-FLAG M2 mouse monoclonal antibody (Sigma).

For vRNA quantification, reverse transcription was performed on 2 µL of input and immunoprecipitated RNA as described above. Real time quantitative PCR was performed as described above with the E_Sarbeco primers. All our analyses were performed as previously described (Marmisolle et al., 2018) using the ΔΔCT method. Briefly, we normalized each RIP fractions' Ct to the corresponding input fraction Ct average for the same qPCR assay (ΔCt [normalized RIP]) to account for sample preparation differences. Then we adjusted the normalized RIP fraction Ct value for the normalized WT Ct value (ΔΔCt = ΔCt [normalized RIP]- ΔCt [normalized HSPB1 negative control]). Finally, we performed a linear conversion of the ΔΔCt (2^[-ΔΔCt]) to calculate the fold change over the HSPB1 protein binding to vRNA.

## Gene editing and trans-complementation experiments

sgRNA targeting QKI were designed using the CRISPOR software (http://crispor.org). Sequences for all the sgRNAs are listed in the key resources table. The sgRNAs were cloned into the plasmid lentiCRISPR v2 (Addgene) according to the recommendations of members of the Zhang laboratory. A549-ACE2 gene-edited cells were generated as previously described (Meertens et al., 2019). Gene-edited cells were transduced with the FLAG-GFP and FLAG-QKI cDNA cloned in the lentiviral vector pLVX-EF1alpha-IRES-Puro.

## Cell viability assay

Cell viability and proliferation were assessed using the CellTiter-Glo 2.0 Assay (Promega). Briefly, cells were plated in 96-well plates (5 × 10³). At specific times, 100 µL of CellTiter-Glo reagent mixed v/v with PBS were added to each well. After 10 min incubation, 90 µL from each well were transferred to an opaque 96-well plate (Cellstar, Greiner bio-one) and luminescence was measured on a TriStar2 LB 942 (Berthold) with 0.1 s integration time.

## QUANTIFICATION AND STATISTICAL ANALYSIS

### Statistical analyses

Graphical representation and statistical analyses were performed using Prism7 software (GraphPad Software). Unless otherwise stated, results are shown as means +/- standard deviation (SD) from at least 2 independent experiments in duplicates. Differences were tested for statistical significance using One-way or two-way ANOVA with multiple comparisons post-test.

Protein functional enrichment was performed in Cytoscape (Version: 3.8.2). Gene Ontology (GO) Biological process were enriched with the ClueGO (Bindea et al., 2009) extension (Version: 2.5.7) using a two-sided hypergeometric test with Bonferroni correction. Adjusted p-values < 0.5 were considered as significant. p-values and percentage of represented genes per GO terms were computed to select the most enriched functional families. This was illustrated with ggplot2 package (Version: 3.3.2) in R Studio interface (Version: 1.2.1335). Most representative terms were represented in networks using the Cytoscape's Genemania extension (Version: 3.5.2). Comparison of functional enrichment between our datasets and previous studies was performed in Metascape (Zhou et al., 2019b). GO Biological process with p-value < 0.01, enrichment score >1.5 and overlap = 4 were selected.