ABSTRACT
We performed an unbiased whole-genome CRISPR/Cas9 screen in A549 cells to identify potential regulators involved in cell death triggered by double-stranded RNA (dsRNA). Of several top candidate genes, we identified the RNA-binding gene ELAV like protein 1 (256529), which encodes the protein Hu antigen R (HuR). Depletion of HuR led to less cell death induced by dsRNA. HuR is mainly involved in apoptosis, and all of its RNA recognition motifs are essential for its pro-apoptotic function. We further showed that the HuR depletion had no influence on the mRNA level of the anti-apoptotic gene BCL2, but instead that HuR downregulates BCL2 translation in a cap-independent way. Polysome fractionation studies showed that HuR retarded the BCL2 mRNA in the non-translating pool of polysomes. Moreover, protection from dsRNA-induced apoptosis by HuR depletion required the presence of BCL2, indicating that the pro-apoptotic function of HuR is executed by suppressing BCL2. Consistent with this, HuR regulated apoptosis induced by infection of encephalomyocarditis or Semliki Forest virus. Collectively, our work identified a suite of proteins that regulate dsRNA-induced cell death, and elucidated the mechanism by which HuR acts as a pro-apoptotic factor.
INTRODUCTION
Programmed cell death (PCD) occurs after virus infections and often functions as a host defense strategy (Danthi, 2016; Orzalli and Kagan, 2017). Three major types of programmed cell death induced by virus infection have been documented – apoptosis, necroptosis and pyroptosis (Danthi, 2016; Imre, 2020). Apoptosis is characterized by cell shrinkage, plasma membrane blebbing, DNA fragmentation and chromosomal condensation. Apoptosis occurs via two distinct pathways, an intrinsic mitochondrial pathway regulated by BCL2 family members and an extrinsic pathway mediated by death receptors of the tumor necrosis factor (TNF) family (Orzalli and Kagan, 2017). Necroptosis is a cell death characterized by cell swelling, plasma membrane rupture and release of cytosolic contents. The necroptosis pathway is mediated through receptor-interacting protein (RIP) kinases (Upton et al., 2017). Pyroptosis exhibits characteristics of both apoptosis and necroptosis. Pyroptosis is carried out by caspase-1, which mediates the maturation of proinflammatory cytokines such as interleukin (IL)-1 and IL-18, and gasdermin D (GSDMD). GSDMD causes holes to form in the plasma membrane, eventually leading to cell death (Man et al., 2017).
Cell death can be triggered by virus attachment and entry, nuclei acids or proteins (Danthi, 2016). dsRNA, as a common replication byproduct of both DNA and RNA viruses (Son et al., 2015), is a shared pathogen-associated molecular pattern that triggers cell death (Harashima et al., 2014; Jacobs and Langland, 1996; Kaufman, 1999) and the cellular innate immune response (Takeuchi and Akira, 2009). Indeed, signaling through pattern recognition receptors including TLR3, RIG-I and MDA5, can promote cell death pathway signaling (Chattopadhyay and Sen, 2017; He et al., 2011; Hiscott et al., 2003). Interferons (IFNs), a downstream product of some pathogen detection sensing pathways, can promote the cell death by transcriptionally regulating the expression of key IFN-stimulated genes (ISGs). Several ISG-encoded proteins, including protein kinases regulated by dsRNA (PKR) (Balachandran et al., 1998; Barber, 2005), 2,5-oligoadenylate synthetase (OAS) and RNase L (Castelli et al., 1998; Diaz-Guerra et al., 1997; Zhou et al., 1997), regulate cell death induced by dsRNA or virus infection (Chawla-Sarkar et al., 2003). Although ISG-induced cell death pathways are well understood, the molecular mechanisms underlying the dsRNA-induced cell death require further delineation.
To identify new regulators of virus-triggered cell death pathway, we performed a genome-wide CRISPR/Cas9-based gene editing screen, using the synthetic dsRNA analog poly(I:C) as a stimulus. We focused study on ELAVL1 (also known as HuR), our top ‘hit’ for promoting dsRNA-induced cell death. ELAVL1 encodes for Hu antigen R (HuR), a member of Hu family protein that is ubiquitously expressed in many tissues (Hinman and Lou, 2008). HuR is an RNA-binding protein comprising three RNA recognition motifs (RRMs) – RRM1 and RRM2 are contiguous, whereas the RRM (RRM3) is separated by a nucleocytoplasmic shuttling sequence. By regulating AU-rich element (ARE)-mediated mRNA stability and protein translation (Brennan and Steitz, 2001), HuR regulates multiple cellular processes, including tumorigenesis, stress responses and programmed cell death (Hinman and Lou, 2008; Zhang et al., 2018). Although a context-dependent role of HuR in programmed cell death and survival has been defined in cancer (Hinman and Lou, 2008), the role of HuR in virus-induced apoptosis or survival remains unknown.
As part of our validation of HuR in promoting the cell death induced by dsRNA, we generated gene-edited cell lines lacking expression of HuR protein. Depletion of HuR resulted in less apoptosis being induced by dsRNA. Mechanistically, HuR downregulated the expression of the anti-apoptotic protein BCL2, and the pro-apoptotic effect of HuR requires the presence of BCL2. Overall, our study shows that HuR is also involved in induction of apoptosis caused by encephalomyocarditis or Semliki Forest virus, two unrelated positive strand RNA viruses.
RESULTS
A CRISPR/Cas9 screen for genes involved in the dsRNA-induced cell death
Transfection of poly(I:C) leads to rapid cell death in many human cell types (Harashima et al., 2014; Kibler et al., 1997). To identify cellular proteins involved in the cell death pathway induced by this dsRNA, we performed a genome-wide CRISPR/Cas9 screen. A549 cells stably expressing Cas9 were transduced with the GeCKOv2 library targeting 19,050 human genes (Sanjana et al., 2014), with each gene targeted by six distinct single-guide RNAs (sgRNAs), and enriched by puromycin selection. Stable sgRNA-expressing A549 cells were subjected to two rounds of poly(I:C) transfection, the second occurring 7 days after the first to ensure all of the control non-library transduced cells underwent death. We expanded the few surviving cells containing sgRNAs in culture and recovered their genomic DNA. We determined the targeted genes followed by PCR amplification and deep-sequencing of the sgRNAs from the integrated proviruses (Fig. 1A).
The top-ranking genes with statistically enriched sgRNAs were determined (Fig. 1B). As expected, genes in the type I IFN system, including RNASEL, JAK1 and STAT1 and STAT2, were at the top of the list, consistent with the established role of IFN signaling in modulating cell death (Barber, 2001). We also analyzed the functions of top 100 enriched genes by GOTERM, and classified them into ten distinct categories, including the glycosaminoglycan biosynthetic process, host cell defense, and type I IFN signaling pathway, among others (Fig. 1C). Genes in the glycosaminoglycan biosynthetic process pathway (SLC35B2, NDST1, B4GALT7, XYLT2 and B3GAT3) were enriched in our screen, likely because glycosaminoglycans facilitate liposome attachment to cells, which enables delivery of poly(I:C) (Payne et al., 2007).
We chose representative genes from each clustered biological process and validated their effects on dsRNA-mediated cell death. Lentiviruses carrying individual sgRNAs were used to establish gene-edited A549 bulk cells. Cells expressing non-targeting sgRNAs served as a negative control, and sgRNA targeting IFNAR and RNASEL genes served as positive controls. The efficiency of gene editing was confirmed by quantitative real-time RT-PCR (qRT-PCR) (Fig. 1D). The gene-edited bulk cells were transfected with poly(I:C) and cell viability was determined 24 h later using an MTT assay. The cell viability of poly(I:C)-transfected IFNAR or RNASEL gene-edited cells were higher than that of control cells (Fig. 1E). The disruption of glycosaminoglycan biosynthesis pathway genes, including B3GALT7, B4GLAT7, NDST7, SCL35B2 and XYLT2 (blue columns), also resulted in higher cell viability. Moreover, the cell viability of ELAVL1-, CAB39-, CYFIP1-, GATA6-, OR8G5- and PDCD10-edited cells was enhanced, which supports the data from our CRISPR/Cas9 screen.
HuR is required for the dsRNA-induced apoptosis
Based on the CRISPR/Cas9 screen data, we pursued the ELAVL1 gene for further investigation, because it was a top ‘hit’ with five of six sgRNAs enriched (Fig. 1B). The ELAVL1-edited bulk cells were relatively resistant to poly(I:C)-induced cell death (Fig. 1E), although the protein encoded by ELAVL1, HuR is not known to regulate the virus-induced cell death. Upon dsRNA transfection, HuR protein was partially transported from nuclei to cytoplasm, implying it is actively enrolled in the cellular responses to dsRNA stimuli (Fig. S1). Next, we explored the functions of HuR in the cell death through a loss-of-function strategy. Two HuR-knockout (KO) cell clones (HuR-KO1 and HuR-KO2) were selected from A549 and HeLa cell lines, and disruption of ELAVL1 gene was confirmed by western blotting and genomic DNA sequencing (Fig. 2A; Fig. S2A). The viability of control A549 or HeLa cells was ∼20% at 24 h post transfection with poly(I:C), indicating that most cells did not survive (black columns, Fig. 2B). By contrast, the viability of all four KO cell clones was ∼40%, about 2-fold higher than the control cells (white and grey columns, Fig. 2B). To confirm a role for HuR in dsRNA-induced cell death, we complemented the HuR-KO A549 cells by introducing a ELAVL1-FLAG fused to a fluorescent reporter gene. The Venus-positive HuR-complemented cells (HuR-RES) were sorted by flow cytometry. The HuR level in the HuR-RES cells was largely restored (Fig. 2C). At 24 h post transfection with poly(I:C), the extent of cell death in the HuR-RES cells was similar to control cells, but greater than in the HuR-KO cells (Fig. 2D). Most control and HuR-RES cells (∼90%) showed characteristics of cell death including membrane blebbing, cell shrinkage and cell detachment, similar to control cells (Fig. 2E). In comparison, the morphological changes of HuR-KO cells were less severe than the control cells. These results demonstrate that HuR promotes the cell death induced by dsRNA.
To determine which type of cell death HuR contributed to, we evaluated cell viability in the presence of the inhibitors of apoptosis [z-VAD(OMe)-FMK; zVAD] or necroptosis (Necrostatin-1; Nec-1), or both inhibitors. Cell viability of the control and HuR-RES cells was enhanced by treatment of zVAD or zVAD plus Nec-1 but not Nec-1 alone (Fig. 2F); in comparison, treatment with zVAD only slightly increased the cell viability of poly(I:C)-transfected HuR-KO cells. We also evaluated cleavage of poly ADP-ribose polymerase (PARP; also known as PARP1), a hallmark feature of apoptosis. Cleavage of PARP in poly(I:C)-transfected control cells was detected at 12 h and increased at 24 h. In the poly(I:C)-treated HuR-KO cells, levels of PARP cleavage were lower than in control cells (Fig. 2G). In addition, we tested whether pyroptosis occurs in A549 cells upon dsRNA treatment by using VX765, an inhibitor of caspase-1, which is a principle effector protease of pyroptosis (Wannamaker et al., 2007). A549 cells were transfected with poly(I:C) in the presence of vehicle (DMSO) or VX765. VX765 treatment did not alter the cell morphological changes or cell viability (Fig. S3A,B), indicating pyroptosis barely occurs in the dsRNA-transfected A549 cells. Overall, these results suggest that HuR mainly regulates the apoptosis pathway of cell death induced by dsRNA.
Moreover, we examined whether HuR plays a role in the apoptosis induced by other classic stimuli, including hydrogen peroxide, doxorubicin and staurosporine. Upon hydrogen peroxide stimulation, a large proportion of control cells underwent apoptosis, while the viability of HuR-KO cells was significantly higher (Fig. S4A,B). In contrast, upon treatment with doxorubicin or staurosporine, PARP cleavage levels in HuR-KO cells were slightly increased from those in control cells (Fig. S4C,D), indicating that HuR plays differential roles in the apoptosis induced by different stimuli.
All three domains of HuR are required for its pro-apoptotic function
To evaluate role of individual domain of HuR in mediating its pro-apoptotic activity, we generated bulk cells expressing truncated HuR lacking the RRM1 (ΔRRM1), RRM2 (ΔRRM2) or RRM3 domain (ΔRRM3) in HuR-KO cells (Fig. 3A). Western blotting data confirmed that all three HuR truncates were expressed (Fig. 3B). Upon treatment with poly(I:C), control and HuR-RES cells showed severe characteristics of cell death, including membrane blebbing, cell shrinkage and cell detachment (Fig. 3C). The morphological changes of ΔRRM1, ΔRRM2 and ΔRRM3 cells were less severe than those in control cells, and were comparable to those in HuR-KO cells (Fig. 3C). The cell viabilities were measured using the MTT assay. As expected, the viability of HuR-RES cells was significantly lower than HuR-KO cells. In contrast, the viabilities of ΔRRM1, ΔRRM2 and ΔRRM3 cells were higher than the control cells, and were comparable to HuR-KO cells, (Fig. 3D), suggesting that all three RRM domains of HuR are required for mediating cell death. The cleaved PARP levels in HuR-KO cells, ΔRRM1, ΔRRM2 and ΔRRM3 cells were comparable, and lower than those in control and HuR-RES cells (Fig. 3E). Deletion of each individual RRM domain of HuR was sufficient to impair the pro-apoptotic activity of HuR, suggesting that all three RRM domains are required for the dsRNA-induced apoptosis regulated by HuR.
HuR inhibits the translation of BCL2
Given an important function of HuR is to stabilize target gene expression at the transcriptional level (Brennan and Steitz, 2001), we performed a RNA sequencing (RNA-seq) assay to define genes transcriptionally altered by HuR expression. However, we did not identify any apoptosis-related genes among the top-ranked list, implying that HuR might not promote dsRNA-induced apoptosis by regulating the mRNA levels of apoptosis-related genes (Fig. S5). We then explored whether HuR regulates expression of apoptosis-related genes. As the mRNA stability and translation of BCL2 has been shown to be regulated by HuR (Ishimaru et al., 2009), we compared protein levels of BCL2 in control, HuR-KO, HuR-RES, ΔRRM1, ΔRRM2 and ΔRRM3 cells. Notably, the protein levels of the BCL2 in HuR-KO cells were markedly increased compared to those in the control and HuR-RES cells. Similarly, the BCL2 level in the ΔRRM1, ΔRRM2 and ΔRRM3 cells was higher than in control cells (Fig. 4A).
As HuR binds to the BCL2 mRNA in vitro and in HL60 cells (Ishimaru et al., 2009), we examined whether HuR is also associated with BCL2 mRNA in A549 cells by performing an RNA immunoprecipitation (RIP) assay. The mRNA levels of negative control GAPDH were comparable in the HuR-KO and HuR-RES A549 cells, while positive control β-actin mRNA level in the HuR-RES cells was higher than the HuR-KO cells (Dormoy-Raclet et al., 2007) (Fig. 4B). The BCL2 mRNA levels associated with HuR in HuR-RES cells was significantly higher than HuR-KO cells, suggesting that HuR binds to BCL2 mRNA in A549 cells.
As HuR KO led to higher level of BCL2 protein, we first tested whether HuR regulates the mRNA level of BCL2. To this end, the BCL2 mRNA levels in control, HuR-KO, HuR-RES, ΔRRM1, ΔRRM2 and ΔRRM3 cells were compared. Loss of expression of full-length or truncated HuR did not impact on the overall cellular levels of BCL2 mRNA (Fig. 4C), which is consistent with our RNA-seq data and suggests that HuR does not affect transcription of BCL2 gene in A549 cells. Next, we examined whether HuR regulates the stability of BCL2 mRNA by using actinomycin D to block the cellular transcription. At indicated time points, total RNA was prepared for qRT-PCR to measure mRNA levels of BCL2 or COX-2, the latter a known target gene stabilized by HuR (Sengupta et al., 2003). As expected, the levels of COX-2 mRNA decreased more rapidly in HuR-KO cells than in control cells. In contrast, the BCL2 mRNA levels were similar in control and HuR-KO cells at all time points (Fig. 4D), suggesting that HuR does not affect the BCL2 mRNA stability in A549 cells.
To investigate whether HuR regulated the translation of BCL2, we evaluated the association of BCL2 mRNA with translating ribosomes by performing a polysome profiling assay. Whole-cell extracts were subjected to sucrose density gradient centrifugation, and the BCL2 mRNA levels in each ribosome fraction were measured by qRT-PCR. The BCL2 mRNA in control cells was distributed mainly in the low molecular mass polysomes (gradient top, fractions 1–7, Fig. 4E); however, this was shifted to high molecular mass translating polysomes (gradient bottom, fractions 9–12) in HuR-KO cells (Fig. 4E). The distributions of control GAPDH mRNAs were similar in the two cell lines. This result indicates that HuR depletion leads to a redistribution of BCL2 mRNA from non-translating polysomes to translating polysomes, thereby increasing its rate of translation.
Considering the HuR regulates the cap-independent translation of various genes directed by internal ribosome entry sites (IRESs) (Kullmann et al., 2002; Winkler et al., 2014), HuR might regulate the BCL2 translation in a similar way. To test this possibility, we assessed whether the BCL2 protein level is affected by rapamycin, an inhibitor of cap-dependent translation. Cells were transfected with poly(I:C) in the presence of DMSO or rapamycin, and harvested at 24 h post transfection. The incubation of rapamycin did not affect the protein levels of BCL2 in either control or HuR-KO cells (Fig. 4F), suggesting that HuR regulates the BCL2 translation in a cap-independent way.
We further tested whether HuR affects the post-translational degradation of BCL2 using inhibitors of two major protein degradation pathways, MG132 (which blocks proteasome degradation; Lee and Goldberg, 1998) and NH4Cl (which blocks lysosomal degradation; Seglen and Reith, 1976). Treatment with MG132 or NH4Cl did not alter the protein levels of BCL2 in control and HuR-KO cells (Fig. 4G), indicating that the reduction of BCL2 protein levels mediated by HuR was independent of proteasome- and autophagy-mediated degradation.
BCL2 plays an essential role in mediating the pro-apoptotic effect of HuR
To explore whether BCL2 plays a role in the pro-apoptotic action of HuR, we utilized a highly selective inhibitor of BCL2, Venetoclax (Ashkenazi et al., 2017). Cells were stimulated with poly(I:C) in the presence of Venetoclax ranging from 0 to 4 µM, and cell viability was measured at 24 h post transfection. As shown in Fig. 5A, the cell viability of all tested cells gradually decreased with increasing Venetoclax concentrations, indicating a protective role of BCL2 in the dsRNA-induced apoptosis. When treated with 1 µM Venetoclax, the cell viability of HuR-KO cells was higher than the control and HuR-RES cells, whereas the cell viability of HuR-KO cells was closer to that of the other two cells in the presence of 2 µM Venetoclax. Treatment of 4 µM Venetoclax resulted in similar low percentages of cell viability (under 10%) in all three cells (Fig. 5A). Treatment of Venetoclax alone did not affect the cell viability (Fig. S6A). In addition, we examined whether other anti-apoptotic proteins, such as Bcl-XL and Mcl-1 are involved in the HuR-regulated apoptosis by using their specific inhibitors (WEHI-539 against Bcl-XL and S63845 against Mcl-1; Lessene et al., 2013; Kotschy et al., 2016). Our data showed that treatment with either WEHI-539 or S63845 did not alter the cell viabilities of control and HuR-KO cells, ruling out the involvement of Bcl-XL and Mcl-1 in the dsRNA-induced apoptosis (Fig. S6B).
To validate that BCL2 mediates the pro-apoptotic effect of HuR, we further generated two gene-edited cells: BCL2-KO cells (BKO), as well as HuR and BCL2 double-knockout cells (dKO). Disruption of the BCL2 and HuR gene in BKO and dKO cells was confirmed by western blot and sequencing (Fig. 5B; Fig. S2B). As expected, the viability of BCL2-KO cells was lower than control cells (Fig. 5C). Knockout of BCL2 in addition to knockout of HuR (i.e. dKO cells) resulted in a lower cell viability than in cells with HuR KO alone (Fig. 5C). The BCL2-KO cells showed more severe morphological changes after poly(I:C) transfection than the other three cells, with greater cell shrinkage and loss of adherence (Fig. 5D). Next, whole-cell lysates were collected for western blotting. As shown in Fig. 5B, PARP cleavage dramatically decreased in poly(I:C)-treated HuR-KO cells compared to control cells. In line with this, BCL2 protein was markedly increased in HuR-KO cells compared to the control cells, consistent with the observation in Fig. 4A. In contrast, PARP cleavage was markedly increased in BKO cells over that in control cells, suggesting that BCL2 has anti-apoptotic activity in these cells. As expected, the PARP cleavage in dKO cells was increased over that in HuR-KO cells, suggesting that HuR promotes poly(I:C)-induced apoptosis through BCL2. Collectively, these results suggest that the pro-apoptotic effect of HuR is mediated in part by BCL2.
HuR promotes apoptosis triggered by virus infection
To explore whether HuR confers a pro-apoptotic activity in the context of virus infection, two single-stranded, positive-sense RNA viruses [encephalomyocarditis virus (EMCV) and Semliki Forest virus (SFV)] were employed, as they produce dsRNA during infection and trigger apoptosis (Bauernfried et al., 2020; El Maadidi et al., 2014; Iordanov et al., 2005; Scallan et al., 1997; Urban et al., 2008). The control and HuR-KO cells were inoculated with EMCV or SFV at a multiplicity of infection (MOI) of 1. qRT-PCR analysis at 3, 6, 9 and 12 h post-infection (h.p.i.) showed comparable viral RNA levels between control and HuR-KO cells (Fig. 6A), indicating that HuR depletion does not alter viral replication. Next, we used morphological observation, flow cytometry and western blot analysis to measure apoptosis. Microscopy images showed that after EMCV or SFV infection, the majority of control cells showed typical characteristics of cell death, including cell membrane blebbing and cell detachment, whereas more HuR-KO cells remained intact. In response to either EMCV or SFV infection, a large proportion of BKO cells underwent shrinkage and cytoplasmic blebbing, and cellular morphological change of dKO cells was more severe than for HuR-KO cells (Fig. 6B). The flow cytometry data showed that, after virus infection, ∼30% and 35% of control cells underwent apoptosis at 24 h.p.i., whereas only ∼12% and ∼15% of HuR-KO cells were apoptotic (Fig. 6C). In response to either EMCV or SFV infection, the levels of PARP cleavage in the BKO cells were the highest among the cells. The cleaved PARP level in the dKO cells was lower than BCL2-KO cells, but higher than HuR-KO cells (Fig. 6D), suggesting that the pro-apoptotic effect of HuR requires the presence of BCL2.
DISCUSSION
Cell death elicited by viral infection affects the interplay between host and virus. In this study, we identified that HuR, an ‘old’ regulator of apoptosis in cancer, plays a ‘new’ role in regulating virus-triggered apoptosis. Our work also illustrated the mechanistic basis of this HuR function – by binding the mRNA of BCL2 and suppressing its translation.
Our CRISPR/Cas9 screen showed that the top 100 genes from the screen were clustered in many biological processes including the type I IFN signaling pathway, which is known to play a role in cell death, confirming that our screen was robust and reliable. We validated and then investigated the function of one of the top hits, ELAV1, encoding HuR, in this study. The disruption of the ELAV1 gene led to lower cell death. Conversely, complementation of HuR in HuR-KO cells restored the extent of cell death induced by dsRNA. dsRNA triggers three main processes that lead to cell death, namely apoptosis, necroptosis and pyroptosis in different cells (Danthi, 2016; Imre, 2020). Our study using an inhibitor of caspase-1 indicated that pyroptosis does not occur in dsRNA-transfected A549 cells (Fig. S3); in line with the notion that pyroptosis mainly occurs in immune cells mediating inflammatory responses (Galluzzi et al., 2018). Moreover, neither apoptosis inhibitor nor necroptosis inhibitor could fully recover the cell viability (Fig. 2F), suggesting that A549 cells undergo multiple types of cell death, such as apoptosis, necroptosis and autophagy. Our data gathered from inhibitor assay, western blot and morphological observation demonstrated that HuR is mainly involved in mediating apoptosis induced by dsRNA. As pyroptosis does not occur, the PARP cleavage in poly(I:C)-treated cells was mediated by caspases in the apoptosis pathway (Erener et al., 2012). Interestingly, HuR also confers a pro-apoptotic activity in response to hydrogen peroxide, but not to doxorubicin or staurosporine, in keeping with previous studies (Badawi et al., 2018; Winkler et al., 2014). HuR plays differential roles in the apoptosis under different contexts, implying that multiple apoptosis signaling pathways, which can be HuR dependent or HuR independent, are employed by cells.
The pro-apoptotic effect of HuR does not appear to rely on the IFN signaling pathway, as the HuR depletion barely affected the levels of IFN-β-encoding and ISG genes. Also, HuR barely affected the mRNA levels of apoptosis-related genes, indicating that its regulatory action is not executed through stabilizing mRNA of apoptosis-related genes.
Our work showed that in A549 cells, HuR binds to the mRNA of BCL2, which encodes an antiapoptotic protein, consistent with a previous report (Ishimaru et al., 2009). HuR lacking each of RRM domains could not restore its pro-apoptotic ability, implying that the RNA-binding activity is crucial. Intriguingly, the association of HuR with BCL2 mRNA does not regulate BCL2 transcription, but instead negatively regulates its translation in an IRES-directed cap-independent way, which differs from the observations gathered from other cells. In HL60 cells, HuR promotes both BCL2 mRNA stability and translation (Ishimaru et al., 2009); in 293T cells, HuR knockdown led to a decrease of endogenous BCL2 expression, but an increase of a BCL2-ARE-reporter transcription, suggesting HuR has opposite effects on endogenous and ectopic BCL2 ARE (Ghisolfi et al., 2009). Therefore, we propose that although HuR binding to BCL2 mRNA is commonly found in many cell types, its regulation on BCL2 transcription and translation is cell specific, directing different levels of BCL2 in different cells, and thereby leading to different phenotypes under physiological stresses.
Our work reveals BCL2 is one of major mediators involved in the pro-apoptotic effect of HuR. The inhibition of BCL2 by a specific inhibitor or disruption of the BCL2 gene increased the vulnerability of A549 cells to cell death upon dsRNA transfection, consistent with idea that BCL2 is protective against dsRNA-induced cell death (Diaz-Guerra et al., 1997; Domingo-Gil and Esteban, 2006). In the absence of BCL2, the pro-apoptotic effect of HuR was substantially diminished, indicating that BCL2 likely mediates the pro-apoptotic function of HuR. Notably, HuR not only binds to BCL2 mRNA, but also interacts with BCL2 protein (Ghisolfi et al., 2009), adding another layer of complexity between HuR and BCL2.
We also demonstrated that HuR contributes to the apoptosis triggered by two RNA viruses, EMCV and SFV. Both EMCV and SFV can produce dsRNA during replication, leading to apoptosis of cells in vitro (El Maadidi et al., 2014; Iordanov et al., 2005; Scallan et al., 1997; Urban et al., 2008) and in vivo in mice (Chattopadhyay et al., 2011). Our study revealed that A549 cells also underwent apoptosis upon infection with EMCV or SFV. As BCL2 overexpression allowed survival of SFV-infected cells (Scallan et al., 1997), we propose that the upregulated BCL2 levels in HuR-KO cells partially contributes to the lower level of apoptosis induced by SFV infection. It is important to note that HuR might also regulate the apoptosis in a BCL2-independent way, as the extent of apoptosis was not fully restored in virus-infected dKO cells (Fig. 6D, lane 8). Nonetheless, HuR does not have a pro-apoptotic role in cell death caused by Dengue virus, VSV and Sendai virus (data not shown), which might be explained by these viruses trigging cytolysis via non-dsRNA factors. For example, viral proteins of DENVs including capsid protein, NS2A and NS2B3/NS3, rather than dsRNA, are the main stimuli of apoptosis (Okamoto et al., 2017).
In summary, our CRISPR/Cas9 screen identified proteins that regulate dsRNA-triggered cell death. In response to dsRNA or viral infection, HuR is transported into cytoplasm, binding to BCL2 mRNA and blocking its translation, eventually promoting the apoptosis. These findings underscore the significance of the association between HuR and BCL2 mRNA, and illustrate a new role of HuR in apoptosis induced by viruses.
MATERIALS AND METHODS
Cells and reagents
Human lung carcinoma epithelial cells (A549, ATCC CCL-185), human cervical cancer cells (HeLa, ATCC CCL-2) and human embryonic kidney cells (293T, ATCC CRL-3216) were maintained in Dulbecco modified Eagle's medium (DMEM) (Gibco) supplemented with 10% fetal bovine serum (FBS; Gibco) at 37°C with 5% CO2. Baby hamster kidney cells (BHK21, ATCC CCL-10) and African green monkey kidney cells (Vero, ATCC CCL-81) were maintained in DMEM supplemented with 5% FBS at 37°C in an incubator with 5% CO2. Aedes albopictus cells (C6/36, ATCC CRL-1660) were maintained in DMEM (Gibco) supplemented with 10% FBS (Gibco) and non-essential amino acids (NEAA; Gibco) at 28°C in an incubator with 5% CO2. The media were supplemented with 100 units/ml of streptomycin and penicillin (Invitrogen).
The inhibitor of apoptosis z-VAD(OMe)-FMK (MCE) and inhibitor of necrosis Necrostatin-1 (MCE) were dissolved in DMSO at a stock concentration of 100 mM. Inhibitors of Bcl-XL (WEHI-539; HY-15607, MCE) and of Mcl-1 (S63845; HY-100741, MCE) were dissolved in DMSO at a stock concentration of 10 mM. The inhibitor of cap-dependent translation rapamycin (MCE) was dissolved in DMSO at a stock concentration of 10 mM.
CRISPR/Cas9 genetic screen
The GeCKO v.2 human CRISPR pooled libraries (A and B) encompassing 123,411 different sgRNA targeting 19,050 genes were purchased from GenScript. Lentiviral production was prepared independently for each half-library in HEK293FT cells by co-transfecting sgRNA plasmids with psPAX2 and pCMV-VSV-G at a ratio of 4:3:1 with FuGENE HD (Promega). After 2 days, culture supernatants were passed through a 0.45 µm filter, and used for gene transduction. A549 cells were transduced with each CRISPR-sgRNA lentiviral library at an MOI of 0.3 and a coverage of 500× for sgRNA representation. Cells were selected with puromycin for 8 days and expanded. Then, 6×107 cells from each library were pooled and transfected with poly(I:C) (1 μg/ml; InvivoGen) using Lipofectamine 2000 (Invitrogen). After 7 days, cells were re-transfected with poly(I:C) for a second round. A complete poly(I:C)-induced cell death was observed in control cells at 7 days after the second transfection. Surviving cells were collected at 7 days post-transfection for genomic DNA extraction using a QIAamp DNA column (Qiagen, Hilden, Germany). The inserted sgRNAs were amplified and sequenced using next-generation sequencing on an Illumina MiSeq (Plateforme MGX, Institut Génomique Fonctionelle). The sgRNA sequences were analyzed using the MAGeCK software and the RIGER software (Jeong et al., 2019; Li et al., 2014).
Generation of bulk gene-edited cells
Gene-edited cell clones were generated using the CRISPR/Cas9 system. sgRNAs for the top-ranking genes were individually cloned into the plasmid lentiCRISPR v2 (Addgene #52961). sgRNAs are listed in Table S1. LentiCRISPR v2 containing sgRNA, pSPAX2 (Addgene #12260), and pMD2.G (Addgene #12259) were introduced into 293T cells using FuGENE® HD Reagent (Promega). After 2 days, culture supernatants were passed through a 0.45 µm filter, and used for gene transduction. A549 cells were transduced with lentiviruses encoding sgRNA. At 24 h post-transduction, 1 µg/ml puromycin was added for selection for 7 days to generate bulk gene-edited cells. To isolate single cell clones, the lentivirus was transduced into A549 cells for 24 h. Then, cells were transferred to a 10-cm dish and selected using 1 μg/ml puromycin for 7 days. Puromycin-resistant clones were selected and confirmed by western blotting and genome DNA sequencing.
Genomic DNA of bulk or cloned cells was extracted using a Cell Genomic DNA Extraction kit (Bioteke, Beijing, China). Regions surrounding the sgRNA target sequences were amplified by PCR. PCR products were cloned into pMD-18T (Takara, Kyoto, Japan) for sequencing. Double KO cells were transduced with two lentiviruses carry sgRNA against each gene and selected as above. The sequencing primers are listed in Table S2.
Induction of apoptosis
Cells were transfected with mock or 1 μg/ml poly(I:C) (Sigma) using Lipofectamine 2000 (Invitrogen). Cells were incubated with 3% hydrogen peroxide, 10 μM doxorubicin or 2 μM staurosporine for 24 h.
MTT and CCK8 assay
Cell viability was measured by means of an 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay. MTT (Genview, Beijing, China) was added to cell culture (5 mg/ml dissolved in PBS) and incubated for 4 h at 37°C. After centrifugation at 200 g for 5 min, the supernatant was discarded and DMSO was added to the plate, and gently shaken for 10 min to dissolve the formazan product. The optical density (OD) value was measured at 490 nm using a BioTek Instrument (BioTeke).
Cell Counting Kit-8 (CCK8) was purchased from MCE. The cells were seeded in a 12-well plate; then, 80 µl reagent was added to cell culture at 24 h post-transfection and cell were incubated for 1 h. Then, the OD value was measured at 450 nm using a BioTek Instrument (BioTek).
Generation of HuR-RES-, ΔRRM1-, ΔRRM2- and ΔRRM3-expressing cells
The full length of HuR or a fragment of HuR with a deleted RRM1 (ΔRRM1), RRM2 (ΔRRM2), and RRM3 (ΔRRM3) gene was amplified by PCR. Primer sequences are listed in Table S3. The amplified fragments were cloned into the vector CSII-EF-MCS-IRES2-Venus (Riken, Kobe, Japan). The lentivirus was packaged in 293T cells through transfecting CSII-EFMCS-IRES2-Venus-HuR, CSII-EFMCS-IRES2-Venus-ΔRRM1, CSII-EFMCS-IRES2-Venus-ΔRRM2, or CSII-EFMCS-IRES2-Venus-ΔRRM3, pSPAX2 and pMD2.G into cells using FuGENE®HD reagent. After 2 days, the culture supernatant was collected and passed through a 0.45 μm filter for gene transduction. HuR-KO cells were transduced with lentivirus carrying the full length of HuR, ΔRRM1, ΔRRM2 or ΔRRM3, and Venus-positive cells were sorted by flow cytometry.
Western blotting
Whole cell extracts were prepared in the presence of 1 mM PMSF, 1% (v/v) protease inhibitor cocktail (Sigma) and dithiothreitol (DDT). Proteins were separated by SDS-PAGE and transferred onto nitrocellulose membranes, followed by blocking in 0.1% PBST (PBS with 0.1% Tween 20) with 5% bovine serum albumin (BSA) (New England Biolabs, Massachusetts, USA), and incubating with primary antibodies at 4°C overnight. Primary antibodies included anti-tubulin (1:5000; RM2007, BBI Life Science, Shanghai, China), anti-HuR (1:400; E0418, Santa Cruz Biotechnology), anti-PARP (1:1000; 9542, Cell Signaling Technology, Massachusetts, USA) and anti-BCL2 (1:1000; YM3041, Immunoway, Delaware, USA). Detection was performed with IRDye 800 CW-conjugated anti-rabbit-IgG or IRDye 680 LT-conjugated anti-mouse-IgG secondary antibody according to the manufacturer's protocols (LI-COR, USA). Immunoreactive bands were visualized using an Odyssey infrared imaging system (LI-COR).
qRT-PCR
Total RNA from A549 cells was extracted using TRIzol reagent (Invitrogen) and reverse transcribed using HI Script Q RT SuperMix (Vazyme, Nanjing, China). qRT-PCR was performed using the cDNA templates and the SYBR Select Master Mix for CFX (Applied Bio systems) in Bio-Rad CFX96 machine. Primers used for qRT-PCR are listed in Table S4. The PCR data were analyzed using SDS software (Applied Biosystems). The level of β-actin was measured as an internal control.
RNA-seq
Total RNA of control and HuR-KO cells was extracted with TRIzol. Three biological replicates were performed. The quality of RNA was assessed by using a Agilent bioanalyzer 2100. An RNA-seq transcriptome library were produced using the NEB Next Ultra RNA library Prep Kit and sequenced on an Illumina NextSeq instrument. R analysis was carried out using the systemPiper workflow. Differential expression analysis was performed using DESeq2 and DEXSeq. Functional-enrichment analysis of the Kyoto Encyclopedia of Gene and Genomes (KEGG) pathways were compared with the whole-transcriptome background.
mRNA stability assay
Cells were treated with DMSO or 0.5 μg/ml actinomycin D (Act D; MCE, Shanghai, China) for 2, 4 or 8 h. Total RNA was extracted for cDNA synthesis, and transcript abundances of BCL2 and COX2 mRNAs were measured by qRT-PCR.
Flow cytometry analysis
Apoptosis was assessed with the PE Annexin V Apoptosis Detection Kit I (Roche, Switzerland) using a flow cytometer (LSRFORTESA, Becton Dickinson, USA). Briefly, A549 cells were stimulated with poly(I:C) for 24 h. The cells were harvested and washed twice in PBS, and then resuspended in 1× binding buffer. 5 µl Annexin V conjugated to phycoerythrin (PE) and 5 µl 7-aminoactinomycin D (7-AAD) were added to a 100 µl cell suspension. Next, the mixture was incubated for 15 min at room temperature and resuspended in 1× binding buffer. Labeled cells were analyzed by flow cytometry (BD Biosciences).
Polysome fractionation
Control and HuR-KO cells were treated with 100 mg/ml cycloheximide at 37°C for 5 min. Whole-cell extracts were prepared and layered onto a 10–50% continuous sucrose gradient and centrifuged at 39,000 rpm in a Beckman SW-41Ti rotor for 3 h at 4°C. Samples were collected from the top of the gradient as 12 fractions. Collected fractions were extracted and subjected to qRT-PCR.
Subcellular fractionation assay
Upon indicated treatment, cells were washed twice with PBS and lysed in cytoplasmic extract buffer [10 mm HEPES, 60 mm KCl, 1 mm EDTA, 0.075% (v/v) Nonidet P-40, 1 mm PMSF and 1 mm DTT]. The mixture was spun at 13,800 g for 5 min to get the cytoplasmic extract (supernatant). The precipitate was washed with cytoplasmic extract buffer and then lysed in nuclear extract buffer [20 mm HEPES, 420 mm NaCl, 10 mm KCl, 1 mm PMSF, 1 mm EDTA, and 20% (v/v) glycerol]. The mixture was spun at 13,800 g for 5 min to get the nuclear extract (supernatant).
Virus, virus infection and titration
EMCV and SFV were provided by Prof. Rongbin Zhou at University of Science and Technology of China and Prof. Xi Zhou at Wuhan Institute of Virology. EMCV and SFV were propagated in HeLa and C6/36 cells, respectively. The EMCV-infected cells were collected when the cytopathic effect appeared, and cells were subjected to three cycles of freezing-and-thawing, followed by centrifugation at 200 g for 5 min. The supernatant of SFV-infected cells was collected when cytopathic effect appeared and subjected to centrifugation and filtration. Virus stocks were titered by a standard plaque assay on BHK21 or Vero cells and stored at −80°C. Serial 10-fold dilutions of virus sample were prepared and inoculated on the cells. The incubation medium was removed and cultured in the mixture of 2×DMEM (Invitrogen) and 2% methylcellulose (Sigma) (1:1). Visible plaques were counted at 2–3 days post infection (EMCV) or 24 hpi (SFV).
For viral infection, A549 cells were incubated with EMCV or SFV at an MOI of 1 for 1 h at 37°C. At the indicated time points post infection, cells were collected for qRT-PCR, western blot or flow cytometry detection.
RNA immunoprecipitation assay
Cells were lysed in RIP lysis buffer [100 mM KCl, 5 mM MgCl2, 0.5% (v/v) NP-40, 10 mM HEPES, 1 mM PMSF, 1% protease inhibitor cocktail (C0001, TargetMol) and RNase inhibitor (2313A, TAKARA)]. Immunoprecipitation was carried out using anti-FLAG M2 Affinity Gel (Sigma). Immune complexes were precipitated with protein A/G agarose as described above. Whole RNA extraction was performed using TRIzol reagent (Invitrogen) according to the manufacturer's protocol. BCL2, β-actin and GAPDH mRNA levels were analyzed by qRT-PCR (primer sequences were listed in Table S4).
Statistical analysis
All experiments were independently repeated at least three times. Comparisons between two groups were performed using Kruskal–Wallis one-way ANOVA or one-way ANOVA with Dunnett's multiple comparison test, or Student's t-tests. Graphs were generated using Graph Pad Prism 6.0 software. P-values of 0.05 or lower were considered to be statistically significant.
Footnotes
Author contributions
Conceptualization: H.G., Y.L., C.H., M.S.D., C.L., R.Z., P.Z.; Methodology: H.G., Y.L., C.H., X.L., P.Z.; Software: C.H., X.L., R.Z.; Validation: H.G., Y.L., C.L.; Formal analysis: H.G., Y.L., C.H., X.L.; Investigation: H.G., Y.L., C.H., M.S.D., P.Z.; Resources: R.Z., P.Z.; Data curation: H.G., Y.L., C.H., C.L.; Writing - original draft: H.G., Y.L., C.H., P.Z.; Writing - review & editing: M.S.D., C.L., R.Z., P.Z.; Supervision: M.S.D., C.L., R.Z., P.Z.; Project administration: M.S.D., C.L., R.Z., P.Z.; Funding acquisition: P.Z.
Funding
This research was supported by National Natural Science Foundation of China (grant number 31970887, 92169110), Guangdong Science and Technology Department (grant number 2018A050506029), and Natural Science Foundation of Guangdong province (grant number 2021A1515011491).
Data availability
RNA-Seq data has been deposited in the Gene Expression Omnibus under GSE165625.
References
Competing interests
M.S.D. is a consultant for Inbios, Vir Biotechnology, NGM Biopharmaceuticals, and Carnival Corporation and on the Scientific Advisory Boards of Moderna and Immunome. The Diamond laboratory has received unrelated funding support in sponsored research agreements from Moderna, Vir Biotechnology and Emergent BioSolutions. Other authors declare no conflict of interest.