Surfaceome CRISPR screen identifies OLFML3 as a rhinovirus-inducible IFN antagonist
Genome Biology volume 22, Article number: 297 (2021)
Rhinoviruses (RVs) cause more than half of common colds and, in some cases, more severe diseases. Functional genomics analyses of RVs using siRNA or genome-wide CRISPR screen uncovered a limited set of host factors, few of which have proven clinical relevance.
Herein, we systematically compare genome-wide CRISPR screen and surface protein-focused CRISPR screen, referred to as surfaceome CRISPR screen, for their efficiencies in identifying RV host factors. We find that surfaceome screen outperforms the genome-wide screen in the success rate of hit identification. Importantly, using the surfaceome screen, we identify olfactomedin-like 3 (OLFML3) as a novel host factor of RV serotypes A and B, including a clinical isolate. We find that OLFML3 is a RV-inducible suppressor of the innate immune response and that OLFML3 antagonizes type I interferon (IFN) signaling in a SOCS3-dependent manner.
Our study suggests that RV-induced OLFML3 expression is an important mechanism for RV to hijack the immune system and underscores surfaceome CRISPR screen in identifying viral host factors.
The emerging genome engineering technologies, including clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein 9 (Cas9), have transformed basic and translational biomedical research . In particular, functional genomics using CRISPR screen provides unprecedented approaches to establishing and understanding phenotype-genotype relationships . For example, CRISPR screen has been widely used to identify and dissect the cellular host factors for a variety of viruses  including noroviruses , human immunodeficiency virus (HIV) , flaviviruses [6, 7], influenza viruses [8, 9], picornaviruses [10,11,12], alphaviruses [13, 14], and others.
RVs are known as the prevalent pathogen causing common cold  and have also been found to be associated with other severe respiratory symptoms including asthma exacerbations  and chronic obstructive pulmonary disease . Despite the increasing number of RV-associated severe respiratory diseases, the causal link between RV infection and clinical outcome remains poorly understood. Particularly, the diverse categories of RVs make it extremely sophisticated to dissect host-pathogen interactions. Functional genomics has been employed to understand RV infections, including RNAi or haploid cell-based genetic screen [18, 19] and emerging CRISPR screens . However, conventional genome-wide genetic screen appeared to have limited efficiency and only a few novel host factors of RVs such as EXOC4 and SETD3 have been identified and validated to have clinical relevance [12, 18].
It has been known that cell proliferation and cell cycle-related proteins may complicate the screening process and analyses of CRISPR genetic screen [20,21,22]. One strategy to overcome this problem is to employ focused screening. To date, a number of focused CRISPR libraries have been constructed to realize genetic screen on kinome [23, 24], epigenome [25, 26], and cancer-related [27, 28] genes. Importantly, cell surface protein-focused CRISPR libraries have been explored in previous studies [11, 29, 30]. However, these studies did not systematically analyze the efficiencies of genome-wide CRISPR screen and surfaceome CRISPR screen. A thorough comparative study can elucidate the difference between genome-wide and surfaceome screens for their efficiencies in identifying viral host factors.
In the present study, we constructed genome-wide and surfaceome CRISPR libraries using identical algorithms and performed genome-wide and surfaceome CRISPR screens in parallel to screen for RV host factors. In contrast to the low success rate of genome-wide CRISPR screen, surfaceome screen identified a set of cell surface host factors that were important for RV infection. Notably, OLFML3 was found to be an RV-inducible dependency factor that promoted the infection. We showed that OLFML3 suppressed the innate immune response of host cells via SOCS-mediated negative regulatory pathway of type I IFN signaling.
Design of sgRNA for genome-wide and surfaceome CRISPR-Cas9 libraries
We constructed CRISPR libraries using previously established algorithms for evaluating on-target  and off-target  activities. The genome-wide library contained 18,421 genes with 12 sgRNA for each gene and was divided into three sub-libraries A, B, and C according to the scores of sgRNA (Additional file 1: Fig. S1a and Additional file 3). The collection of cell surface proteins was defined using a mass spectrometry (MS) database that investigated cell surface proteins in 41 individual human cells . The union set contains 1344 genes and the 12 sgRNA targeting to these genes were extracted from the genome-wide library to construct the surfaceome library (Additional file 1: Fig. S1a and Additional file 4). Cross reference of the 1344 cell surface proteins by gene ontology dataset for cellular component showed that the majority of these proteins are located on plasma membrane and/or extracellular matrix (Fig. 1a).
Construction of genome-wide and surfaceome CRISPR-Cas9 libraries in H1-Hela cells
To normalize for the quality of sgRNAs in the library, we performed side-by-side screening with RVs using genome-wide sub-library A and surfaceome library, which exhibited similar on-target and off-target scores (Additional file 1: Fig. S1b). This sub-library A contains 73,527 sgRNA targeting to 18,421 genes with an average of four sgRNA per gene. In comparison, surfaceome CRISPR screen involved 16,975 sgRNA targeting to 1344 genes. These libraries were packaged into lentiviruses (LVs) and transduced with a multiplicity of infection (MOI) of 0.3 into H1-Hela cells that are known to support RV infection  (Fig. 1b). Next-generation sequencing (NGS) analyses of the PCR amplicons of genome-integrated sgRNA revealed that plasmid and cell libraries of genome-wide and surfaceome screens had full coverage of sgRNA with optimum distribution, as evidenced by Gini coefficients of less than 0.1  (Additional file 1: Fig. S2a). The cumulative frequency of sgRNA reads in these four libraries were consistent with the predicted distribution (Additional file 1: Fig. S2b-c).
In order to have further quality control of genome-wide and surfaceome screens, we compared the shift of sgRNAs targeting to the essential and nonessential genes . We evaluated the enriched and depleted sgRNAs at 10 days after LV transduction in comparison with the baseline sgRNAs at 3 days after LV transduction. MAGeCK analysis showed that non-targeting sgRNAs or sgRNAs targeting to nonessential genes have significantly different fold change from sgRNAs targeting to essential genes at the population level (Additional file 1: Fig. S3a and Additional files 5, 6). In addition, the receiver operating characteristic (ROC) analysis [37, 38] showed that both genome-wide and surfaceome screens had good performance in differentiating essential genes from nonessential genes, with the area under the curve (AUC) values of 0.964 and 0.963 respectively (Additional file 1: Fig. S3b). Collectively, these results validated the quality of the constructed CRISPR libraries for subsequent RV challenge.
Identification of RV host factors on H1-Hela using genome-wide and surfaceome CRISPR-Cas9 libraries
To identify RV host factors, we performed the screen using a previously reported RV-induced cell death model . H1-Hela cells harboring genome-wide or surfaceome CRISPR library were challenged with RV serotype B14 (RV-B14) at an MOI of 1.0 for 48 h. Upon completion of RV-B14 challenge, survived cells were collected and sgRNAs were PCR-amplified and then analyzed by NGS. MAGeCK analyses  of the RV infection groups in comparison with mock groups (Fig. 1b) identified a set of enriched sgRNA (Additional file 1: Fig. S3c). Using a false discovery rate (FDR) of less than 0.01 as a cutoff for sgRNA enrichment, surfaceome screen identified 3.2% of the total 1344 surface proteins bearing enriched sgRNA, while genome-wide screen uncovered sgRNA enrichment in 0.54% of the 18,421 genes in the library (Fig. 1c). Of the enriched sgRNA in the genome-wide screen, less than 20% were targeted to the surface proteins (Fig. 1d). It was noted that surfaceome screen identified more absolute number of enriched sgRNA for surface proteins than genome-wide screen (Fig. 1d).
Next we determined the genes with enriched sgRNA using modified robust rank aggregation (a-RRA) analyses in the MAGeCK pipeline (MAGeCK score) and displayed top 10 candidate hits in each library (Fig. 1e and Additional files 7, 8). In order to justify MAGeCK and the choice of FDR, we also analyzed the sgRNA data with CB2, which is a beta-binomial distribution-based algorithm for analyzing sgRNA shift during CRISPR screen . Using an FDR cutoff of 0.001, MAGeCK and CB2 identified 45 and 19 genes respectively from the genome-wide library, and 17 and 52 genes respectively from the surfaceome library (Additional file 1: Fig. S4a-b and Additional files 9, 10). We found that a considerable fraction of the top 10 hits from MAGeCK analysis could be also detected by CB2 method. These results justified the use of MAGeCK analysis and the choice of top 10 hits as representatives for validation. In addition, we found that genome-wide and surfaceome screens identified different sets of hit genes for surface proteins (Additional file 1: Fig. S4c), thus supporting the use of surfaceome screen as an alternative screening strategy.
Importantly, both genome-wide and surfaceome screens identified intercellular adhesion molecule 1 (ICAM-1), a known receptor for RV serotypes A and B , in the top 10 hits (Fig. 1e). Except for ICAM-1, these screens resulted in no overlap in the top 10 hits. It has to be noted that the surfaceome and genome-wide CRISPR libraries contain 12 and 4 sgRNA respectively for each gene. To exclude the possibility that the differential hits in genome-wide and surfaceome screens were due to the effects of different number of sgRNAs, we performed a pseudo screen with the surfaceome library where sgRNA 1 to 4 in the sub-library A (Additional file 1: Fig. S1), rather than all the 12 sgRNA, were included for each surface protein. MAGeCK analyses of the pseudo screen revealed a generally consistent rank of the top 10 hits with the experimental screen (Additional file 1: Fig. S4d). In addition, we found there is minor difference between genome-wide and surfaceome libraries for the distribution of on-target and off-target scores of sgRNAs of the top 10 hits (Additional file 1: Fig. S4e). Collectively, our results suggested that both genome-wide and surfaceome screens were robust enough to identify strong candidates such as RV receptor and that surfaceome screen was more efficient in identifying candidate hits for surface proteins.
Validation of the top 10 candidate hits in genome-wide and surfaceome CRISPR screens
Next we sought to validate the top 10 hits from genome-wide and surfaceome screens by constructing knockout cell lines for individual genes. Two sgRNAs were designed for each gene and the gene disruption efficiency was determined by T7E1 assay and Sanger sequencing (Additional file 1: Fig. S5 and Additional file 11). Non-targeting sgRNA was used as mock for subsequent phenotypical analyses. In similar experimental settings with the virus challenge screen, we examined RV-induced cell death on individual cell lines of the top 10 hits identified from surfaceome and genome-wide screens. It was found that in the genome-wide screen only ICAM-1 knockout exhibited consistently improved cell viability with both sgRNAs compared to non-targeting sgRNA, whereas surfaceome screen identified 6 gene knockout out that displayed protective effects against RV-14 infection (Fig. 2a). Next we analyzed the effects of knockout of the identified genes on viral replication. Individual knockout cells were infected with RV-14 at an MOI of 2 and viral RNA was extracted from medium supernatant or cell lysates at 24 h post infection. Among the top 10 hits from genome-wide screen, only ICAM-1 knockout resulted in significantly reduced viral loads in medium supernatant in comparison with non-targeting sgRNA (Fig. 2b). Investigation of viral loads in cell lysates uncovered generally consistent results (Additional file 1: Fig. S6a). These results suggested that except for ICAM-1 the top 10 candidate genes identified by genome-wide screen were most likely false positive hits as RV host factors. By contrast, the 6 gene knockout from surfaceome screen that exhibited protective effects (Fig. 2a) showed consistently reduced viral loads in medium supernatant (Fig. 2b). The function of these 6 genes was further examined by immunofluorescence staining of RV-B14 envelope protein for evaluation of the effects of gene knockout on viral infection rate. It was found that in addition to ICAM-1, the knockout of IGSF3, RAB5C, and OLFML3 showed consistently reduced RV-B14 infection rates with both sgRNA at an MOI of 1 or 2 (Fig. 2c and Additional file 1: Fig. S6b-c).
To further validate the function of identified genes from surfaceome CRISPR screen, we constructed single clones for knocking out ICAM-1, RAB5C, OLFML3, IGSF3, SLC4A7, and ATP6AP1 respectively. Except for ATP6AP1, we obtained single clones with knockout at both alleles for all identified surface proteins (Additional file 1: Fig. S7 and Additional file 12). ATP6AP1 is an important ATP-driven proton pump , the knockout of which might be lethal to cells . The following experiments were thus performed with ATP6AP1−/+ cells. It was found that ICAM-1−/−, RAB5C−/−, OLFML3−/−, IGSF3−/−, SLC4A7−/−, and ATP6AP1−/+ cells all exhibited markedly reduced viral loads in medium supernatant (Fig. 2d) and cell lysates (Fig. 2e) upon RV-B14 infection, as compared with non-targeting sgRNA-treated H1-Hela cells. Collectively, these results confirmed that the surface proteins ICAM-1, RAB5C, OLFML3, IGSF3, SLC4A7, and ATP6AP1 identified from surfaceome CRISPR screen played important roles during RV infection.
Validation of RAB5C and OLFML3 as RV dependency factors
Although a series of candidate dependency factors were identified for RV-B14 from surfaceome screen, we chose to focus subsequent analyses on RAB5C and OLFML3. RAB5C was reported to involve in the entry process of flaviviruses . OLFML3 was found to function as an immunosuppressive molecule during the carcinogenesis of glioblastoma . These two proteins may represent novel dependency factors of RV-B14, none of which have been reported to involve in RV infection. As a known receptor for RV serotypes A and B [45, 46], ICAM-1 was included as positive control for the following experiments.
It was found that ICAM-1−/−, RAB5C−/−, and OLFML3−/− cells all exhibited significantly increased cell viability upon RV-B14 challenge compared with non-targeting sgRNA-treated H1-Hela cells (Additional file 1: Fig. S8a-b). These knockout cells showed markedly reduced viral loads in medium supernatant (Fig. 3a) and cell lysates (Fig. 3b) over monitored time course, where ICAM-1 or RAB5C knockout nearly abolished viral replication (Fig. 3a-b). To further validate the importance of RAB5C and OLFML3 for RV infection, we performed rescue experiments by overexpressing these genes in corresponding knockout cells. It was found that RAB5C and OLFML3 overexpression rescue could restore RV-B14-induced cell death in H1-Hela knockout cells (Fig. 3c). Importantly, consistent results were observed with a different RV serotype RV-A16 (Fig. 3d). RAB5C and OLFML3 overexpression also rescued the replication of RV-B14 or RV-A16 in corresponding knockout cells, as determined by viral loads in medium supernatant or cell lysates (Fig. 3e–h). Similarly, ICAM-1 overexpression could rescue the susceptibility of ICAM-1−/− cells to RV infection (Additional file 1: Fig. S8c-f). Collectively, these results demonstrate that RAB5C and OLFML3 are host dependency factors of RV.
Next we sought to examine the effects of RAB5C or OLFML3 knockout on the infection of clinical RV strain. A RV strain was isolated from the oropharyngeal swab of a patient with upper respiratory tract infection. This RV isolate induced notable cytopathic effects (CPEs) in H1-Hela cells (Fig. 3i) as laboratory strains did. ICAM-1−/−, RAB5C−/−, and OLFML3−/− H1-Hela cells all exhibited significant resistance to the clinical RV strain (Fig. 3i, j). Phylogenetic analyses using the VP4 gene suggested that the clinical RV isolate was closely related with the serotype A RV with less similarity with RV-B and RV-C strains (Fig. 3k). Further investigation of viral loads in medium supernatant (Fig. 3l) and cell lysates (Fig. 3m) suggested that ICAM-1, RAB5C, or OLFML3 knockout abolished or inhibited the replication of the clinical RV strain. These results thus verified the consistent functions of RAB5C and OLFML3 during the infection of laboratory and clinical RV strains.
Dissection of the functions of RAB5C and OLFML3 in RV infection
To analyze the functions of RAB5C and OLFML3, we first sought to determine whether these proteins were involved in virus attachment and entry. For virus attachment assay, RV-B14 was incubated with mock and knockout cell lines at 4 °C and cells with attached virus were harvested and lysed for RT-qPCR quantification of viral loads. Consistent with the known function of ICAM-1 as the receptor of RV serotypes A and B, ICAM-1−/− cells exhibited significantly reduced virus attachment. In contrast, knockout of RAB5C or OLFML3 did not lead to reduced virus attachment (Additional file 1: Fig. S9a). For virus entry assay, RV-B14 was incubated with mock and knockout cell lines first at 4 °C and then the internalization of RV-B14 was initiated by incubation at 37 °C. Surface-bound RV-B14 was removed by extensive washing and internalized virus was quantified by RT-qPCR. It was found that RAB5C−/− or OLFML3−/− cells did not show reduced virus entry (Additional file 1: Fig. S9b). Consistent results were observed with the functions of RAB5C and OLFML3 in RV-A16 attachment and entry (Additional file 1: Fig. S9c-d). Fluorescence in situ hybridization (FISH) detection of internalized RV-B14 showed similar findings for the roles of RAB5C and OLFML3 in virus entry (Additional file 1: Fig. S9e-f). These results together suggested that RAB5C or OLFML3 did not affect RV attachment or entry.
Next we sought to analyze whether RAB5C and OLFML3 were involved in the life cycle of RV after uncoating of viral genome. The RNA encoding RV-A16 genome was transfected into cells to bypass virus entry and uncoating processes. It was found that mock and knockout cells exhibited similar cell viability (Additional file 1: Fig. S9g) and viral loads (Fig. 4a, b). This suggested that RAB5C or OLFML3 was not involve in the processes after genome uncoating. This experiment along with the above virus attachment and entry assays excluded the functions of RAB5C or OLFML3 in many processes during RV replication. In order to investigate whether RAB5C and OLFML3 participated in RV genome RNA uncoating, we used a previously established guanidine hydrochloride (GuHCl)-mediated mRNA synthesis suppression assay  to analyze the dynamics of intracellular viral RNA. Mock and knockout H1-Hela cells displayed comparable cell viability in the presence of 2 mM GuHCl (Additional file 1: Fig. S9h). During the monitored time course of 24 h, RV-B14 RNA increased in a time-dependent manner in wild-type H1-Hela cells. By contrast, treatment of cells with 2 mM GuHCl led to notable time-dependent RNA decay (Additional file 1: Fig. S9i), suggesting of inhibited synthesis of viral RNA and active cellular RNA degradation machinery. In the presence of 2 mM GuHCl, RV-B14 and RV-A16 infection in ICAM-1−/− cells resulted in little viral loads in cell lysates at all time points (Fig. 4c, d), which was in consistency with the above results of ICAM-1−/− cells (Fig. 3a, b). With 2 mM GuHCl treatment, the dynamics of viral RNA in OLFML3−/− cells displayed similar pattern to that of non-targeting sgRNA-treated cells, while RAB5C−/− cells had significantly higher viral RNA level at 3 and 6 h post RV infection (Fig. 4c, d). Given inhibited RNA synthesis, the slower decay rate of viral RNA in RAB5C−/− cells could be attributed to the protection of viral RNA from cellular RNA degradation machinery as a result of perturbed uncoating process of viral genome. Collectively, these results suggested that RAB5C, but not OLFML3, was involved in RV genome uncoating.
To further validate the role of RAB5C during RV infection, we performed transcriptome-wide RNA sequencing (RNA-Seq) analyses with wild-type and RAB5C−/− cells. The RNA-Seq data were collected from three independent experimental replicates and were found to have high degree of correlation (Additional file 1: Fig. S10a). Analyses of mapped read counts of RV-B14 genome showed significant inhibition of viral replication in RAB5C−/− cells (Additional file 1: Fig. S10b), consistent with above results of RAB5C knockout-mediated inhibition of viral replication. Analyses of differentially expressed genes (DEG) showed a large number of downregulated genes upon RV-B14 challenge in RAB5C−/− cells (Additional file 1: Fig S10c-d). Gene ontology (GO) analyses of the biological processes of DEGs revealed significantly altered cellular response to virus infection in RAB5C−/− cells (Additional file 1: Fig. S10e). Collectively, RNA-Seq results confirmed that RAB5C played an important role during RV infection.
Identification of OLFML3 as a negative regulator of innate immune response
The above results showed that OLFML3 was a critical RV dependency factor but did not participate in RV attachment, entry, or genome uncoating. This prompted us to elucidate the mechanism of action of OLFML3 in RV infection using RNA-Seq. RNA-Seq analyses were performed with wild-type and OLFML3−/− cells before and after RV-B14 infection. The RNA-Seq data from three independent experimental replicates were collected and the replicates within each condition were found to have high degree of correlation (Additional file 1: Fig. S11a). Analyses of mapped read counts of RV-B14 genome revealed significant, transcriptome-wide inhibition of viral gene expression in OLFML3−/− cells (Additional file 1: Fig. S11b), consistent with the above results of OLFML3 knockout-mediated inhibition of viral replication (Fig. 3). Analyses of DEGs showed that removal of OLFML3 induced the upregulation of a broad range of genes upon RV-B14 challenge (Fig. 4e and Additional file 1: Fig. S11c). GO analyses of the biological processes and molecular function of DEGs revealed markedly changed innate immune response (Fig. 4f and Additional file 1: Fig. S11d). To validate the RNA-Seq results, we infected non-targeting sgRNA-transduced and OLFML3−/− cells with RV-B14 and RV-A16 at an MOI of 2 for 24 h. RT-qPCR quantification uncovered a series of upregulated IFN-stimulating genes (ISGs) in OLFML3−/− cells, but not in mock cells, upon RV infection (Fig. 4g, h and Additional file 1: Fig S12). These results were consistent with the RNA-Seq analyses and strongly indicated that OLFML3 promoted RV infection in H1-Hela cells by antagonizing innate immune response.
OLFML3 is a RV-inducible IFN suppressor
Although our results have suggested a possible role of OLFML3 in the innate immune response, the signaling pathway OLFML3 is involved in remains elusive. OLFML3 is secreted glycoprotein consisting of approximately 400 amino acids. OLFML3 belongs to the olfactomedin (OLF) superfamily and bears a C-terminal olfactomedin-like (OLFML) domain (Fig. 5a). In human, there are five OLFML members and the functions of OLFML1, OLFML2A, OLFML2B, and OLFML4 have been illustrated . However, the role of OLFML3 is poorly understood and there have been no studies reporting its functions with viral infection.
To understand the function of OLFML3, we quantified its expression in H1-Hela. Interestingly, OLFML3 had very low mRNA expression under uninfected conditions, and RV-B14 and RV-A16 infection upregulated the expression of OLFML3 by more than 500 and 400 folds respectively (Fig. 5b). Importantly, in the presence of OLFML3, ISGs could not be efficiently activated by RV infection over a course of 12 h. By contrast, removal of OLFML3 allowed time-dependent upregulation of ISGs in response to RV infection (Fig. 5c). These observations are consistent with the above results (Fig. 4) and demonstrated that OLFML3 promoted RV infection (Fig. 3) by suppressing innate immune response. Moreover, IFN and STAT1/2 expression underwent minor changes in response to RV infection in the presence of OLFML3 while OLFML3 knockout sensitized IFN and STAT1/2 upregulation in response to RV infection (Fig. 5d), which was consistent with the RNA-Seq results (Fig. 4). Collectively, our results suggested that RV-induced OLFML3 activation and OLFML3-mediated inhibition of type I IFN signaling might be a critical mechanism for RV to escape the innate immune system.
OLFML3 antagonizes type I IFN signaling in an SOCS3-dependent mechanism
Consistent with the RNA-Seq results (Fig. 4), OLFML3 knockout significantly reduced the expression of suppressor of cytokine signaling-3 (SOCS3) (Fig. 6a, b), a well characterized suppressor of IFN signaling , at 24 h post RV infection. SOCS3 knockdown using siRNA (Fig. 6b) prevented RV replication in H1-Hela cells (Fig. 6c). Interestingly, SOCS3 knockdown seemed more efficient than OLFML3 knockout in inhibiting viral replication (Fig. 6c). Meanwhile, SOCS3 knockdown and OLFML3 knockout did not have additive effects (Fig. 6c). These results suggested that OLFML3 and SOCS3 were in the same signaling pathway and that SOCS3 functioned as the predominating molecule or was downstream of OLFML3.
Next we analyzed the effects of SOCS3 knockdown on the expression of IFN signaling molecules. It was found that SOCS3 knockdown had more prominent impacts than OLFML3 knockout on activating IFN signaling molecules including STAT1/2, IFNB, and ISGs. In consistency with the above viral infection, SOCS3 knockdown in OLFML3−/− cells did not further enhance type I IFN signaling (Fig. 6d). These results consistently suggested that SOCS3 was a downstream molecule of OLFML3 that predominated the negative regulatory pathway of IFN signaling. Therefore, our results collectively suggested that OLFML3 exerted its IFN-inhibiting activity through OLFML3-SOCS3-STAT1/2 axis (Fig. 6e).
The widespread RV infections and its remarkable phenotypic diversity have precluded the development of effective vaccines and antiviral therapeutics. This obstacle has rendered host-directed therapy (HDT) , an attractive option for treating RV infections, where host-virus interactions are interfered. HDT of RVs requires systemic dissection of cellular host factors supporting viral infection. Most importantly, as potential drug targets, the identified host factors should be readily accessible to antiviral therapeutics including macromolecular drugs. This prompted us to investigated surface protein-focused functional genomics approaches for uncovering druggable cell surface host factors of RVs. It is surprising to find that genome-wide and surfaceome CRISPR screens identified different sets of surface proteins. Among the top 10 identified hits, 6 out of 10 from the surfaceome screen were validated, in comparison with 1 out of 10 from the genome-wide screen. These results suggest that surfaceome screen could be used as an efficient alternative screening approach for functional genomics. However, it must be noted that the surfaceome library reported in this study contains only approximately 1400 surface proteins and is thus not guaranteed to cover all the genes of interest.
The identified RV host factor RAB5C is an endosome-localized small GTPase and is deemed to be involved in cellular trafficking . It has been reported that RAB5C plays critical roles during the infection of flaviviruses including Zika virus and Dengue virus  as well as other viruses . However, the roles of RAB5C during RV infection remain elusive. We found that RAB5C did not affect RV attachment or entry nor did it affect the replication of RV genome. Further analyses suggested that RAB5C might involve in the uncoating process of RV genome. Considering the reported cellular localization and function of RAB5C, our results indicate that RAB5C may be a key regulator for the endosomal release of RV genome.
Unlike RAB5C, OLFML3 does not participate in RV genome uncoating nor virus attachment or entry. OLFML3, also known as hOLF44, is one of the five OLFML proteins in human. OLFML1 is found to be associated with cell proliferation and autonomy in human cancer cells. OLFML2A and OLFML2B function as photomedin proteins. OLFML4 is known to have anti-inflammatory and anti-apoptotic activities. However, existing studies do not clearly define the biological functions of OLFML3. Some studies have suggested that OLFML3 may be involved in development and tumorigenesis  but none has reported its role in viral infection. In the present study, we found that OLFML3 promotes RV infection by acting as a type I IFN antagonist. The expression of OLFML3 is induced by RV infection, which may be an important mechanism for RV to escape from the innate immunity of host cells. We have shown that OLFML3-mediated inhibition of type I IFN signaling is dependent on SOCS3, a known suppressor of IFN signaling. Our study is the first to define OLFML3 as an IFN signaling inhibitor and to connect the functions of OLFML proteins to viral infections. Nevertheless, it is not completely unreasonable to expect the activities of OLFML3 in viral infections given the general roles of other OLFML proteins in inflammation and apoptosis and the signaling molecules these pathways share.
In future studies, it would be interesting to explore how OLFML3 initiates the negative regulatory pathway of IFN signaling. Particularly, it is important to investigate whether OLFML3 triggers the signaling through certain cellular receptors or whether OLFML3 directly interfere with the interactions between IFN and IFN receptors. In addition, it remains unclear how OLFML3 affects the expression of SOCS3. The fact the OLFML3 affected RV-activated SOCS3 expression at late time point (24 h after infection) suggested an indirect or delayed interaction between OLFML3 and SOCS3. Furthermore, it will be interesting to investigate whether OLFML3 functions as an antagonist of IFN signaling during the infection of other enteroviruses.
Finally, since OLFML3 is a secreted glycoprotein that can be readily accessed by drug molecules, it is worth developing small molecule or antibody therapeutics targeting OLFML3. On the one hand, targeting OLFML3 with drug molecules can help elucidate its signaling pathway during RV infection. On the other hand, blocking OLFML3 with drug molecules can help assess the feasibility of OLFML3 as a drug target for treating RV infection. Importantly, because we have shown that OLFML3 does not function as attachment or entry receptor, simply blocking the solvent-exposed surface of OLFML3 may not be sufficient for therapeutics to have anti-RV effects. Thus, a screening platform may be needed to identify functional epitopes on OLFML3 for drug design and development.
In the present study, we have shown that surfaceome CRISPR screen outperforms genome-wide screen in identifying RV host factors. Surfaceome screen has identified OLFML3 as a RV-inducible IFN suppressor that exerts the IFN-inhibiting activity through the OLFML3-SOCS3-STAT1/2 axis. Our study has thus underscored surfaceome CRISPR screen for rapid dissection of host-pathogen interactions.
Design and construction of genome-wide and surfaceome CRISPR libraries
sgRNAs were designed to target to protein-coding regions (NCBI CCDS data, released on 8-Sep-2016)  and optimized by two steps. First, off-target scores were calculated according to an established algorithm . Second, on-target scores were calculated using Rule Set 2 . sgRNAs were ranked by on-target scores and the top 12 sgRNAs with off-target scores of less than 20 were selected for each gene. If less than 12 sgRNAs were obtained, the cutoff of off-target scores was increased sequentially to 40, 60, 80, and 100 until 12 sgRNAs were obtained. The genome-wide CRISPR library was divided into three sub-libraries A, B, and C according to sgRNA rank. Non-targeting sgRNAs were included in each library.
Pooled sgRNA oligonucleotides were synthesized as 76-mers by Custom Array (Bothell, WA, USA) and were amplified by PCR with NEBNext High-Fidelity PCR Master Mix (New England BioLabs, NEB, Ipswich, MA, USA) using customized primers (Additional file 2: Table S1). The PCR products were purified using MinElute PCR purification kit (Qiagen, Dusseldorf, Germany). Lentiviral vector LentiCRISPR-v2 was digested with Esp3I (Thermo Fisher Scientific, Waltham, MA, USA) at 37 °C for 3 h and gel-purified using Gel Extraction kit (Omega, Norcross, GA, USA). Purified digestion products were ligated to Esp3I-treated LentiCRISPR-v2 using Gibson assembly kit (NEB) following the manufacturer’s instructions. The ligation product was purified by isopropanol precipitation and then transformed into electrocompetent Escherichia coli (Lucigen, Middleton, WI, USA). Transformed cells were plated on to 15-cm Luria-Bertani (LB) agar plates supplemented with 50 μg/mL ampicillin (Sangon Biotech, Shanghai, China). Approximately 1–3 × 107 colonies were collected for each library to ensure 500-fold coverage. Plasmid DNA was extracted as pooled libraries using NucleoBond Xtra Maxi EF kit (Macherey-Nagel, Duere, Germany) and stored at − 20 °C.
H1-Hela cells were obtained from the American Type Culture Collection (ATCC). HEK293T cells were obtained from the Cell Bank of Shanghai Institutes for Biological Science (SIBS) and were validated by VivaCell Biosciences (Shanghai, China). H1-Hela and HEK293T were grown in Dulbecco’s modified Eagle’s medium (DMEM, Thermo) supplemented with 10% fetal bovine serum (FBS, Thermo) and 1% penicillin-streptomycin (Thermo) and maintained at 37 °C in a fully humidified incubator containing 5% CO2. All cells were confirmed by PCR to be free of mycoplasma contamination.
LV production and transduction
To produce LVs with high titers, HEK293T cells were seeded on to 6-well plates with 2 × 106 cells per well for single sgRNA, or 10-cm petri dishes with 107 cells for library construction. At 24 h after seeding, HEK293T cells at a confluence of 70–90% were transfected with LV packaging plasmid pMD2.G, envelope plasmid psPAX, and transfer plasmid pLentiCRISPR-v2 that carries single or pooled sgRNAs using Lipofectamine 3000 (Thermo). At 6 h after transfection, the medium was replaced with fresh medium. The medium supernatant containing LVs was harvested at 48 h post transfection by centrifugation at 2000 rpm for 10 min, filtrated through a 0.45-μm filter and stored at − 80 °C.
H1-Hela cells were transduced with LVs at an MOI of 0.3 using spinfection. Briefly, H1-Hela cells were washed with phosphate-buffered saline (PBS) and incubated with LVs in serum-free DMEM under centrifugation at 2000 rpm for 2 h. Upon completion of spinfection, LV-containing medium was removed and cells were incubated in DMEM supplemented with 10% FBS and 2 μg/mL puromycin (Thermo) for 2 to 3 days to purge empty cells containing no LVs or sgRNAs. For each library, cells of more than 500-fold coverage of the library size were collected.
RV production and infection
The full-length cDNA clones of RV-A16 (pR16.11, Cat. No. VRMC-8) and RV-B14 (pWR3.26, Cat. No. VRMC-7) were obtained from ATCC. To produce infectious viral RNA, RV-A16 and RV-B14 plasmids were linearized by SacI (NEB) digestion and then in vitro transcribed using HiScribe T7 Transcription Kit (NEB). The RNA transcripts were extracted using Trizol (Thermo) and chloroform (Titan, Shanghai, China), followed by isopropanol precipitation. Viral RNA was transfected into H1-HeLa cells using Lipofectamine 3000 (Thermo) to generate infectious RV-A16 or RV-B14 particles. At 48 h post transfection, the supernatant containing RVs were collected for further infection on H1-Hela cells to produce RVs with higher titers. The aliquots of purified RVs were stored at − 80 °C. Virus titers were determined by the 50% tissue culture infectious dose (TCID50) assay.
For RV infection, H1-Hela cells were seeded on to 96- or 12-well plates with a density of 2 × 104 or 1.5 × 104 cells per well, respectively. Unless noted otherwise, at 24 h after seeding, cells were infected with RVs at an MOI of 2 for 1.5 h, washed with PBS for three times and then cultured in DMEM (Thermo) supplemented with 10% fetal bovine serum (FBS, Thermo) and 1% penicillin-streptomycin (Thermo) for 24 h.
CRISPR screen using RV-B14
H1-Hela cells carrying genome-wide sub-library A or surfaceome library were seeded on to 15-cm petri dishes. Approximately 1.5 × 107 cells were seeded to ensure more than 200-fold coverage of sgRNA. These cells were infected with RV-14 at an MOI of 1 for 48 h. At the end point of RV-14 challenge, the cells were washed with PBS and attached cells were collected from the plates. Two biological replicates were performed for mock and test groups. Genomic DNA of the harvested cells was extracted using phenol: chloroform: isoamyl alcohol (v/v/v, 25:24:1) and then purified using ethanol precipitation.
Next-generation sequencing (NGS) analyses of sgRNA enrichment
Genome-integrated sgRNAs were amplified from extracted genomic DNA by PCR using the primers containing Illumina adaptor (Additional file 2: Table S1). The PCR product was gel-purified and then analyzed on Illumina HiSeq 3000 platform by Genewiz (Suzhou, Jiangsu, China). After removing the adaptors, the 20 bp sgRNA was mapped to the reference sgRNA libraries with one nucleotide mismatch allowed for each sgRNA. Gini index was calculated to analyze the distribution of sgRNAs. The raw read counts were subjected to MAGeCK analyses  and CB2 analyses  to determine the enriched sgRNA and genes. Enrichment of sgRNAs and genes was analyzed using MAGeCK (v0.5.7) and CB2 (v1.3.4) by comparing the read counts from the cells infected with RV-B14 with those from uninfected cells. A false discovery rate (FDR) of less than 0.01 was applied to identify significantly enriched sgRNAs or gene knockout.
Generation of CRISPR-Cas9 knockout cells
The transfer plasmid pLentiCRISPR-v2 carrying single sgRNA for targeted gene knockout was constructed as described above. Briefly, forward and reverse oligonucleotides encoding the 20-bp sgRNA (Additional file 2: Table S2) were annealed to generate double-stranded DNA with overhang that matched the sticky ends of Esp3I (Thermo)-treated pLeniCRISPR-v2 vector. Annealed sgRNA sequence was ligated into digested pLeniCRISPR-v2 and then transformed into DH5α E. coli (Tsingke, Beijing, China).
The LVs carrying single sgRNA were packaged and transduced on to cells as described above. To evaluate the knockout efficiency, the genomic DNA of edited cells were extracted using Quick Extraction kit (Lucigen). Target sites carrying gene-edited sequences were PCR amplified using gene-specific primers (Additional file 2: Table S3). The knockout efficiency of each sgRNA was determined using T7E1 analysis and Sanger sequencing. Sanger sequencing results were analyzed by Synthego website (https://tools.synthego.com/). Single clones were obtained by limited dilution and genotyped by Sanger sequencing to determine the mutations at each allele. Knockout of target proteins was verified by western blotting (WB).
Cell viability assay
H1-Hela cells were seeded into 96- or 24-well plates with a density of 5000 or 50,000 cells per well respectively. At 24 h after seeding, cells were infected with RVs at an MOI of 1 for 1.5 h, washed with PBS for three times, and cultured in fresh medium for 24 h unless noted otherwise. Cell counting Kit-8 (CCK-8, Dojindo, Kumamoto, Japan) was applied to determine cell viability according to the manufacturer’s instructions. The absorbance at 450 nm was determined by Enspire multimode plate reader (PerkinElmer, Waltham, MA, USA).
Immunofluorescence (IF) and fluorescence in situ hybridization (FISH)
IF staining of RV-B14 was performed using mouse anti-RV VP3 antibody (1:50, clone G47A, Thermo) and Alexa Fluor Plus 488 goat anti-mouse IgG (H + L) (1:1000, A32723, Thermo). IF images were acquired and analyzed using Operatta high-content analysis system (PerkinElmer). At least 2000 fluorescent cells were imaged and quantified for each replicate.
FISH was performed using an RNAscope Multiplex Fluorescent V2 Assay kit (Advanced Cell Diagnostics, Newark, USA) according to the manufacturer’s instructions. After fixation and pretreatment, RV-B14 RNA was detected using an RVB RNA probe (Advanced Cell Diagnostics, Cat. No. 447141) and TSA Plus Fluorescein (PerkinElmer, Cat. No. NEL741E001KT). FISH images were acquired and analyzed using a TissueFAXS 200 flow-type tissue quantitative analyzer (TissueGnostics GmbH, Vienna, Austria). At least 5000 cells in each replicate were included in analyses.
Real-time quantitative PCR (RT-qPCR)
To determine viral loads, the medium supernatant or cell lysate-containing RVs were harvested at 24 h after infection with RVs at an MOI of 2 unless noted otherwise. The total RNA from supernatant or cell lysate was purified using Trizol (Thermo) and chloroform (Titan), followed by purification using isopropanol precipitation. Purified viral RNA (vRNA) was reverse transcribed into cDNA using PrimeScript RT reagent Kit with gDNA Eraser (Takara Bio Inc., Shiga, Japan). The number of RV genome copy in the medium supernatant was determined using RT-qPCR with general or serotype-specific Taqman probe and primers (Additional file 2: Table S4) on Applied Biosystems Q6 Real-Time PCR cycler. The absolute viral titers were calculated based on a standard curve of RV genome with known TCID50, and the R square of curve-fitting was guaranteed to be more than 0.99. The mRNA levels of RVs and IFN-stimulating genes (ISGs) in cell lysate were determined using RT-qPCR with SYBR green dye (Thermo) and specific primers (Additional file 2: Table S5) on Applied Biosystems Q6 Real-Time PCR cycler. All SYBR Green primers were validated with dissociation curves. The expression of vRNA and host genes in cell lysate is normalized to ribosomal gene RPLP0 (36b4).
Gene knockdown using siRNA
H1-Hela cells were seeded on to 6-well plates with a density of 5 × 105 cells per well. At 24 h after seeding, cells were transfected with 100 pmol SOCS3 siRNA (Genepharma, Shanghai, China) (Additional file 2: Table S6) using 7.5 μL Lipofectamine 2000 (Thermo) for 6 h, washed with PBS, and then cultured in fresh DMEM (Thermo) supplemented with 10% fetal bovine serum (FBS, Thermo). At 48 h post transfection, cells were infected with RV-B14 at an MOI of 2 for 1.5 h, washed with PBS for three times and then cultured in fresh medium for 24 h. The medium supernatant or cell lysate-containing RVs were harvested and lysed for total RNA extraction, and then the mRNA levels of RVs, SOCS3, STAT1/2, and IFN-stimulating genes (ISGs) in cell lysate were determined using RT-qPCR as described above.
Virus attachment and entry assays
H1-Hela cells were seeded in a 12-well plate with a density of 200,000 cells per well and incubated overnight. For virus attachment assay, cells were incubated with RV-B14 or RV-A16 at an MOI of 20 in cold medium on ice for 60 min, washed by PBS for three times and then the total RNA was extracted using Trizol (Thermo). For virus entry assay, cells were incubated with RV-B14 or RV-A16 at an MOI of 20 in cold medium on ice for 60 min, washed by PBS for three times, and then treated with pre-warmed medium for 40 min at 37 °C. Then cells were washed with PBS for three times and treated with 0.25% trypsin for 2 min (Thermo, Cat. No. 25200072) to remove surface-bound viral particles. The internalized viral RNA was extracted using Trizol (Thermo) and viral loads were determined using RT-qPCR and FISH.
Rescue experiments by overexpression
ICAM1, RAB5C, and OLFML3 genes were codon-optimized for expression in human cells and synthesized by Genewiz. The 20 bp sgRNA-targeting sites and PAM sequences were mutated with silent mutations. Myc and FLAG tags were added to the C-terminus of these genes for WB detection. These genes were cloned into the EcoRI and XhoI sites of pCAGG plasmid that carries a separate mScarlet fluorescent protein as a transfection reporter.
For WB analysis, cells were lysed with RIPA buffer (Beyotime Biotechnology, Beijing, China) on ice for 10 min. The total protein concentration in cell lysate was determined using the BCA Protein Assay Kit (Thermo). Cell lysate was mixed with SDS-PAGE loading buffer (Takara) containing 200 mM DTT, incubated at 95 °C for 10 min, and resolved on NuPAGE 4-12% Bis-Tris gels (Thermo). Protein samples were transferred onto nitrocellulose membranes using an iBlot gel transfer system (Thermo). The following primary and secondary antibodies were used in WB: anti-ICAM1 rabbit antibody (Cell Signaling Technology, Cat. No. 4915S, Danvers, USA), anti-RAB5C rabbit antibody (Thermo, Cat. No. PA551932), anti-OLFML3 rabbit antibody (Thermo, Cat. No. PA531581), HRP-conjugated anti-rabbit IgG (CST, Cat. No. 7074S). Anti-β actin antibody conjugated with HRP (Abcam, Cat. No. ab49900, Cambridge, UK) was used as an internal control.
RNA-Seq analysis of mock and knockout cells
Non-targeting sgRNA-treated mock cells, OLFML3−/− cells and RAB5C−/− were seeded on to 10-cm plates at a density of 1 × 106 cells per plate. At 24 h after seeding, the cells were infected with RV-B14 at an MOI of 1 or with PBS as mock infection. After 24 h of RV challenge, the remaining cells were washed with PBS and collected by Trizol (Thermo) treatment. Three biological triplicates were prepared for each group. The whole-transcriptome RNA sequencing was performed by Genewiz.
RNA-Seq short reads were aligned to the human genome (GRCh37) using Hisat2 (v2.0.1) . Gene expression was counted as the number of short reads fully or partially aligned to the annotated gene model using HTSeq (v0.6.1) . Genotype (non-targeting sgRNA-transduced mock cells and OLFML3−/− cells) and treatment conditions (mock and RV-B14 infection) were the factor variables in our RNA-Seq data. DESeq2 (v1.26.0)  was used to examine the difference between the response of mock and knockout cells to RV infection, which was captured by the interaction term. RV infection-induced gene upregulation and downregulation were calculated, and the differentially upregulated or downregulated genes in mock and knockout cells were defined as DEGs. Significant DEGs were filtrated with an adjusted P value of less than 0.05 and a fold change value of more than 2.
GO enrichment analysis was performed using cluster Profiler (v3.14.3) by comparing DEGs to a list of all human genes . Adjusted P value of less than 0.001 or 0.05 was set as the filter for biological process and molecular function terms respectively.
Isolation of clinical RV strain
Nasopharyngeal swab samples were collected from the hospitalized patients bearing respiratory infection symptoms such as fever, cough, pharyngalgia, and others. Nasopharyngeal swab samples were maintained in viral-transport medium and all samples were conserved on − 80 °C until analyses. The swab samples were diagnosed for respiratory viruses and the remaining samples were used for isolation of RVs. All procedures were complied with the Measures for the Ethical Review of Biomedical Research Involving Human Subjects issued by the National Health and Family Planning Commission of the People’s Republic of China. The Ruijin Hospital Ethics Committee, Shanghai Jiaotong University School of Medicine, approved the sample collection protocol with a permit number of Ruijin Hospital Ethics Committee 2018-48.
For isolation of RVs, the specimens were centrifuged at 2000 rpm for 10 min and the supernatant was collected for RT-qPCR analysis using universal RV primers (Supplementary Table 4). The samples with high RV loads were selected for subsequent RV isolation in H1-Hela cells. H1-Hela cells were seeded on to 12-well plates at a cell density of 100,000 cells per well. Specimen supernatant (700 μL) was mixed with 300 μL of fresh medium and incubated with cells under centrifugation at 2000 rpm for 2 h. Upon completion of spinfection, the supernatant was replaced with fresh medium and cells were incubated for another 2 days. Thereafter, medium supernatant and attached cells were collected and subjected to quick freeze-thaw cycles for three times to release viral particles. The medium supernatant or cell lysate from above was added to cells on 12-well plates. At 24 h after incubation, the supernatant was replaced with fresh medium and cells were incubated for 2 to 3 days before reaching 100% confluence. Isolation of RV strain was confirmed by CPEs and sequencing results.
All data are the results from at least three biological replicates and are shown as mean ± SD unless noted otherwise. Statistical analyses and graphing were performed with GraphPad Prism 7.0. The P values were determined using two-tailed unpaired Student’s t test unless otherwise noted.
Availability of data and materials
The raw NGS data, including sgRNA library and transcriptome sequencing data, have been deposited into NCBI SRA (Sequence Read Archive) database with the accession codes SRP26965 (https://www.ncbi.nlm.nih.gov/sra/?term=SRP269650) . Additional information or materials supporting the findings in this study can be provided upon reasonable request to the corresponding author J. Liu.
Doudna JA. The promise and challenge of therapeutic genome editing. Nature. 2020;578(7794):229–36. https://doi.org/10.1038/s41586-020-1978-5.
Doench JG. Am I ready for CRISPR? A user’s guide to genetic screens. Nat Rev Genet. 2018;19(2):67–80. https://doi.org/10.1038/nrg.2017.97.
Puschnik AS, Majzoub K, Ooi YS, Carette JE. A CRISPR toolbox to study virus-host interactions. Nat Rev Microbiol. 2017;15(6):351–64. https://doi.org/10.1038/nrmicro.2017.29.
Orchard RC, Wilen CB, Doench JG, Baldridge MT, McCune BT, Lee YC, et al. Discovery of a proteinaceous cellular receptor for a norovirus. Science. 2016;353(6302):933–6. https://doi.org/10.1126/science.aaf1220.
Park RJ, Wang T, Koundakjian D, Hultquist JF, Lamothe-Molina P, Monel B, et al. A genome-wide CRISPR screen identifies a restricted set of HIV host dependency factors. Nat Genet. 2017;49(2):193–203. https://doi.org/10.1038/ng.3741.
Zhang R, Miner JJ, Gorman MJ, Rausch K, Ramage H, White JP, et al. A CRISPR screen defines a signal peptide processing pathway required by flaviviruses. Nature. 2016;535(7610):164–8. https://doi.org/10.1038/nature18625.
Marceau CD, Puschnik AS, Majzoub K, Ooi YS, Brewer SM, Fuchs G, et al. Genetic dissection of Flaviviridae host factors through genome-scale CRISPR screens. Nature. 2016;535(7610):159–63. https://doi.org/10.1038/nature18631.
Li B, Clohisey SM, Chia BS, Wang B, Cui A, Eisenhaure T, et al. Genome-wide CRISPR screen identifies host dependency factors for influenza A virus infection. Nat Commun. 2020;11(1):164. https://doi.org/10.1038/s41467-019-13965-x.
Han J, Perez JT, Chen C, Li Y, Benitez A, Kandasamy M, et al. Genome-wide CRISPR/Cas9 screen identifies host factors essential for influenza virus replication. Cell Rep. 2018;23(2):596–607. https://doi.org/10.1016/j.celrep.2018.03.045.
Staring J, von Castelmur E, Blomen VA, van den Hengel LG, Brockmann M, Baggen J, et al. PLA2G16 represents a switch between entry and clearance of Picornaviridae. Nature. 2017;541(7637):412–6. https://doi.org/10.1038/nature21032.
Zhao X, Zhang G, Liu S, Chen X, Peng R, Dai L, et al. Human neonatal Fc receptor is the cellular uncoating receptor for enterovirus B. Cell. 2019;177(6):1553–1565 e1516. https://doi.org/10.1016/j.cell.2019.04.035.
Diep J, Ooi YS, Wilkinson AW, Peters CE, Foy E, Johnson JR, et al. Enterovirus pathogenesis requires the host methyltransferase SETD3. Nat Microbiol. 2019;4(12):2523–37. https://doi.org/10.1038/s41564-019-0551-1.
Zhang R, Kim AS, Fox JM, Nair S, Basore K, Klimstra WB, et al. Mxra8 is a receptor for multiple arthritogenic alphaviruses. Nature. 2018;557(7706):570–4. https://doi.org/10.1038/s41586-018-0121-3.
Meertens L, Hafirassou ML, Couderc T, Bonnet-Madin L, Kril V, Kummerer BM, et al. FHL1 is a major host factor for chikungunya virus infection. Nature. 2019;574(7777):259–63. https://doi.org/10.1038/s41586-019-1578-4.
Jacobs SE, Lamson DM, St George K, Walsh TJ. Human rhinoviruses. Clin Microbiol Rev. 2013;26(1):135–62. https://doi.org/10.1128/CMR.00077-12.
Castillo JR, Peters SP, Busse WW. Asthma exacerbations: pathogenesis, prevention, and treatment. J Allergy Clin Immunol Pract. 2017;5(4):918–27. https://doi.org/10.1016/j.jaip.2017.05.001.
Seemungal T, Harper-Owen R, Bhowmik A, Moric I, Sanderson G, Message S, et al. Respiratory viruses, symptoms, and inflammatory markers in acute exacerbations and stable chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2001;164(9):1618–23. https://doi.org/10.1164/ajrccm.164.9.2105011.
Perreira JM, Meraner P, Brass AL. Functional genomic strategies for elucidating human-virus interactions: will CRISPR knockout RNAi and haploid cells? Adv Virus Res. 2016;94:1–51. https://doi.org/10.1016/bs.aivir.2015.11.001.
Perreira JM, Aker AM, Savidis G, Chin CR, McDougall WM, Portmann JM, et al. RNASEK is a V-ATPase-associated factor required for endocytosis and the replication of rhinovirus, influenza A virus, and Dengue virus. Cell Rep. 2015;12(5):850–63. https://doi.org/10.1016/j.celrep.2015.06.076.
Haapaniemi E, Botla S, Persson J, Schmierer B, Taipale J. CRISPR-Cas9 genome editing induces a p53-mediated DNA damage response. Nat Med. 2018;24(7):927–30. https://doi.org/10.1038/s41591-018-0049-z.
Brown KR, Mair B, Soste M, Moffat J. CRISPR screens are feasible in TP53 wild-type cells. Mol Syst Biol. 2019;15(8):e8679. https://doi.org/10.15252/msb.20188679.
Bowden AR, Morales-Juarez DA, Sczaniecka-Clift M, Agudo MM, Lukashchuk N, Thomas JC, et al. Parallel CRISPR-Cas9 screens clarify impacts of p53 on screen performance. Elife. 2020;9. https://doi.org/10.7554/eLife.55325.
Wang C, Wang H, Lieftink C, du Chatinier A, Gao D, Jin G, et al. CDK12 inhibition mediates DNA damage and is synergistic with sorafenib treatment in hepatocellular carcinoma. Gut. 2020;69(4):727–36. https://doi.org/10.1136/gutjnl-2019-318506.
Tian R, Gachechiladze MA, Ludwig CH, Laurie MT, Hong JY, Nathaniel D, et al. CRISPR interference-based platform for multimodal genetic screens in human iPSC-derived neurons. Neuron. 2019;104(2):239–255 e212. https://doi.org/10.1016/j.neuron.2019.07.014.
Li F, Huang Q, Luster TA, Hu H, Zhang H, Ng WL, et al. In vivo epigenetic CRISPR screen identifies Asf1a as an immunotherapeutic target in Kras-mutant lung adenocarcinoma. Cancer Discov. 2020;10(2):270–87. https://doi.org/10.1158/2159-8290.CD-19-0780.
Xu G, Chhangawala S, Cocco E, Razavi P, Cai Y, Otto JE, et al. ARID1A determines luminal identity and therapeutic response in estrogen-receptor-positive breast cancer. Nat Genet. 2020;52(2):198–207. https://doi.org/10.1038/s41588-019-0554-0.
Inoue D, Chew GL, Liu B, Michel BC, Pangallo J, D'Avino AR, et al. Spliceosomal disruption of the non-canonical BAF complex in cancer. Nature. 2019;574(7778):432–6. https://doi.org/10.1038/s41586-019-1646-9.
Bersuker K, Hendricks JM, Li Z, Magtanong L, Ford B, Tang PH, et al. The CoQ oxidoreductase FSP1 acts parallel to GPX4 to inhibit ferroptosis. Nature. 2019;575(7784):688–92. https://doi.org/10.1038/s41586-019-1705-2.
Chong ZS, Ohnishi S, Yusa K, Wright GJ. Pooled extracellular receptor-ligand interaction screening using CRISPR activation. Genome Biol. 2018;19(1):205. https://doi.org/10.1186/s13059-018-1581-3.
Ye L, Park JJ, Dong MB, Yang Q, Chow RD, Peng L, et al. In vivo CRISPR screening in CD8 T cells with AAV-Sleeping Beauty hybrid vectors identifies membrane targets for improving immunotherapy for glioblastoma. Nat Biotechnol. 2019;37(11):1302–13. https://doi.org/10.1038/s41587-019-0246-4.
Doench JG, Fusi N, Sullender M, Hegde M, Vaimberg EW, Donovan KF, et al. Optimized sgRNA design to maximize activity and minimize off-target effects of CRISPR-Cas9. Nat Biotechnol. 2016;34(2):184–91. https://doi.org/10.1038/nbt.3437.
Hsu PD, Scott DA, Weinstein JA, Ran FA, Konermann S, Agarwala V, et al. DNA targeting specificity of RNA-guided Cas9 nucleases. Nat Biotechnol. 2013;31(9):827–32. https://doi.org/10.1038/nbt.2647.
Bausch-Fluck D, Hofmann A, Bock T, Frei AP, Cerciello F, Jacobs A, et al. A mass spectrometric-derived cell surface protein atlas. PLoS One. 2015;10(4):e0121314. https://doi.org/10.1371/journal.pone.0121314.
Xie J, Yea K, Zhang H, Moldt B, He L, Zhu J, et al. Prevention of cell death by antibodies selected from intracellular combinatorial libraries. Chem Biol. 2014;21(2):274–83. https://doi.org/10.1016/j.chembiol.2013.12.006.
Li W, Koster J, Xu H, Chen CH, Xiao T, Liu JS, et al. Quality control, modeling, and visualization of CRISPR screens with MAGeCK-VISPR. Genome Biol. 2015;16(1):281. https://doi.org/10.1186/s13059-015-0843-6.
Hart T, Tong AHY, Chan K, Van Leeuwen J, Seetharaman A, Aregger M, et al. Evaluation and design of genome-wide CRISPR/SpCas9 knockout screens. G3 (Bethesda). 2017;7:2719–27.
Grau J, Grosse I, Keilwagen J. PRROC: computing and visualizing precision-recall and receiver operating characteristic curves in R. Bioinformatics. 2015;31(15):2595–7. https://doi.org/10.1093/bioinformatics/btv153.
Ong SH, Li Y, Koike-Yusa H, Yusa K. Optimised metrics for CRISPR-KO screens with second-generation gRNA libraries. Sci Rep. 2017;7(1):7384. https://doi.org/10.1038/s41598-017-07827-z.
Jeong HH, Kim SY, Rousseaux MWC, Zoghbi HY, Liu ZD. Beta-binomial modeling of CRISPR pooled screen data identifies target genes with greater sensitivity and fewer false negatives. Genome Res. 2019;29(6):999–1008. https://doi.org/10.1101/gr.245571.118.
Bella J, Kolatkar PR, Marlor CW, Greve JM, Rossmann MG. The structure of the two amino-terminal domains of human ICAM-1 suggests how it functions as a rhinovirus receptor and as an LFA-1 integrin ligand. Proc Natl Acad Sci U S A. 1998;95(8):4140–5. https://doi.org/10.1073/pnas.95.8.4140.
Abbas YM, Wu D, Bueler SA, Robinson CV, Rubinstein JL. Structure of V-ATPase from the mammalian brain. Science. 2020;367:1240- +.
Hart T, Chandrashekhar M, Aregger M, Steinhart Z, Brown KR, MacLeod G, et al. High-resolution CRISPR screens reveal fitness genes and genotype-specific cancer liabilities. Cell. 2015;163(6):1515–26. https://doi.org/10.1016/j.cell.2015.11.015.
Savidis G, McDougall WM, Meraner P, Perreira JM, Portmann JM, Trincucci G, et al. Identification of Zika virus and Dengue virus dependency factors using functional genomics. Cell Rep. 2016;16(1):232–46. https://doi.org/10.1016/j.celrep.2016.06.028.
Chen P, Hsu WH, Chang A, Tan Z, Lan Z, Zhou A, et al. Circadian regulator CLOCK recruits immune-suppressive microglia into the GBM tumor microenvironment. Cancer Discov. 2020;10(3):371–81. https://doi.org/10.1158/2159-8290.CD-19-0400.
Staunton DE, Merluzzi VJ, Rothlein R, Barton R, Marlin SD, Springer TA. A cell adhesion molecule, ICAM-1, is the major surface receptor for rhinoviruses. Cell. 1989;56(5):849–53. https://doi.org/10.1016/0092-8674(89)90689-2.
Greve JM, Davis G, Meyer AM, Forte CP, Yost SC, Marlor CW, et al. The major human rhinovirus receptor is ICAM-1. Cell. 1989;56(5):839–47. https://doi.org/10.1016/0092-8674(89)90688-0.
Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35(6):1547–9. https://doi.org/10.1093/molbev/msy096.
Jin Y, Li JL. Olfactomedin-like 3: possible functions in embryonic development and tumorigenesis. Chin Med J. 2019;132(14):1733–8. https://doi.org/10.1097/CM9.0000000000000309.
Rothlin CV, Ghosh S, Zuniga EI, Oldstone MB, Lemke G. TAM receptors are pleiotropic inhibitors of the innate immune response. Cell. 2007;131(6):1124–36. https://doi.org/10.1016/j.cell.2007.10.034.
Schwegmann A, Brombacher F. Host-directed drug targeting of factors hijacked by pathogens. Sci Signal. 2008;1:re8.
Barbera S, Nardi F, Elia I, Realini G, Lugano R, Santucci A, et al. The small GTPase Rab5c is a key regulator of trafficking of the CD93/Multimerin-2/beta1 integrin complex in endothelial cell adhesion and migration. Cell Commun Signal. 2019;17(1):55. https://doi.org/10.1186/s12964-019-0375-x.
Ziegler CM, Bruce EA, Kelly JA, King BR, Botten JW. The use of novel epitope-tagged arenaviruses reveals that Rab5c-positive endosomal membranes are targeted by the LCMV matrix protein. J Gen Virol. 2018;99(2):187–93. https://doi.org/10.1099/jgv.0.001004.
Pruitt KD, Harrow J, Harte RA, Wallin C, Diekhans M, Maglott DR, et al. The consensus coding sequence (CCDS) project: Identifying a common protein-coding gene set for the human and mouse genomes. Genome Res. 2009;19(7):1316–23. https://doi.org/10.1101/gr.080531.108.
Li W, Xu H, Xiao T, Cong L, Love MI, Zhang F, et al. MAGeCK enables robust identification of essential genes from genome-scale CRISPR/Cas9 knockout screens. Genome Biol. 2014;15(12):554. https://doi.org/10.1186/s13059-014-0554-4.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60. https://doi.org/10.1038/nmeth.3317.
Anders S, Pyl PT, Huber W. HTSeq--a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9. https://doi.org/10.1093/bioinformatics/btu638.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550. https://doi.org/10.1186/s13059-014-0550-8.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7. https://doi.org/10.1089/omi.2011.0118.
Mei H, Zha Z, Wang W, Xie YS, Huang YG, Li WP, et al. Surfaceome CRISPR screen identifies OLFML3 as a rhinovirus-inducible IFN antagonist: NCBI SRA (Sequence Read Archive); January 1, 2021. https://www.ncbi.nlm.nih.gov/sra/?term=SRP269650
We thank the High-throughput Screening (HTS) Platform at Shanghai Institute for Advanced Immunochemical Studies (SIAIS), ShanghaiTech University for the support of cell viability, IF and FISH experiments and Biomedical Big Data Platform for the design of sgRNA libraries, MAGeCK analyses, and analyses of RNA-Seq data. We thank Professor Jincun Zhao (Guangzhou Medical University) and Professor Jia Xie (University of Miami) for providing RV-B14 and RV-A16 genome plasmids.
The review history is available as Additional file 13.
Peer review information
Kevin Pang was the primary editor of this article and managed its editorial process and peer review in collaboration with the rest of the editorial team.
This work is supported by National Key R&D Program of China International Collaboration Project (2018YFE0200402 to J.L.), Zhejiang University special scientific research fund for COVID-19 prevention and control (2020XGZX011 to J.Q. and J.L.), China Postdoctoral Science Foundation (2017 M621551 to H.M), Shanghai Municipal Key Clinical Specialty (shslczdzk02202 to J.Q.), Shanghai Top-Priority Clinical Key Disciplines Construction Project (2017ZZ02014 to J.Q.), Shanghai Shenkang Hospital Development Center Clinical Science and Technology Innovation Project (SHDC12018102 to J.Q.), National Innovative Research Team of High-level Local Universities in Shanghai (to J.Q.), and ShanghaiTech University Startup Fund (2019F0301-000-01 to J.L.).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The original version of this article was revised: Additional file 1 has been updated to correct the author affiliations.
Supplementary Figures. Fig. S1. Construction of CRISPR genome-wide and surfaceome libraries. Fig. S2. Quality analyses of constructed genome-wide and surfaceome CRISPR libraries. Fig. S3. Evaluation of surfaceome and genome-wide CRISPR libraries. Fig. S4. Validation of the screening results. Fig. S5. Determination of gene modification efficiency. Fig. S6. Validation of the top 10 hits from surfaceome and genome-wide screens. Fig. S7. Construction and validation of single clones of ICAM-1-/-, RAB5C-/-, OLFML3-/-, SLC4A7-/- and ATP6AP1-/+ H1-Hela cells. Fig. S8. Validation of the effects of ICAM-1, RAB5C and OLFML3 on RV infection, related to Fig. 3. Fig. S9. Dissection of the functions of RAB5C and OLFML3 in RV infection. Fig. S10. RNA-seq analyses of the effects of RAB5C knockout on RV infection. Fig. S11. RNA-Seq analyses of the effects of OLFML3 on RV infection (related to Fig. 4). Fig. S12. Bar plots showing RT-qPCR quantification of ISG expression in mock and OLML3-/- cells at 24 h post infection of RV-B14 (a) and RV-A16 (b) at an MOI of 2.
Supplementary Tables. Table S1. Primers for construction of CRISPR library and RNA-Seq analyses. Table S2. Primers for construction of sgRNA plasmid. Table S3. Primers for PCR amplification of sgRNA targeted sites for T7E1 analyses. Table S4. RV-B14-specific or RV universal primers for Taqman RT-qPCR. Table S5. Primers for RT-qPCR. Table S6. siRNA sequences. Table S7. Codon-optimized DNA sequences of OLFML3 gene. Table S8. Codon-optimized DNA sequences of RAB5C gene.
Design of genome-wide CRISPR library A.
Design of surfaceome CRISPR library.
MAGeCK analysis of essential and nonessential genes of genome-wide screen.
MAGeCK analysis of essential and nonessential genes of surfaceome screen.
MAGeCK analysis of genome-wide screen after RV challenge.
MAGeCK analysis of surfaceome screen after RV challenge.
CB2 analysis of genome-wide screen after RV challenge.
CB2 analysis of surfaceome screen after RV challenge.
Raw data of Sanger sequencing chromatography.
Uncropped WB images.
About this article
Cite this article
Mei, H., Zha, Z., Wang, W. et al. Surfaceome CRISPR screen identifies OLFML3 as a rhinovirus-inducible IFN antagonist. Genome Biol 22, 297 (2021). https://doi.org/10.1186/s13059-021-02513-w