- Short report
- Open Access
Enhanced susceptibility of cancer cells to oncolytic rhabdo-virotherapy by expression of Nodamura virus protein B2 as a suppressor of RNA interference
Journal for ImmunoTherapy of Cancervolume 6, Article number: 62 (2018)
Antiviral responses are barriers that must be overcome for efficacy of oncolytic virotherapy. In mammalian cells, antiviral responses involve the interferon pathway, a protein-signaling cascade that alerts the immune system and limits virus propagation. Tumour-specific defects in interferon signaling enhance viral infection and responses to oncolytic virotherapy, but many human cancers are still refractory to oncolytic viruses. Given that invertebrates, fungi and plants rely on RNA interference pathways for antiviral protection, we investigated the potential involvement of this alternative antiviral mechanism in cancer cells. Here, we detected viral genome-derived small RNAs, indicative of RNAi-mediated antiviral responses, in human cancer cells. As viruses may encode suppressors of the RNA interference pathways, we engineered an oncolytic vesicular stomatitis virus variant to encode the Nodamura virus protein B2, a known inhibitor of RNAi-mediated immune responses. B2-expressing oncolytic virus showed enhanced viral replication and cytotoxicity, impaired viral genome cleavage and altered microRNA processing in cancer cells. Our data establish the improved therapeutic potential of our novel virus which targets the RNAi-mediated antiviral defense of cancer cells.
Oncolytic viruses (OVs) possess an intrinsic or engineered ability to selectively target, replicate in and kill cancer cells . These promising cancer therapeutics exploit cellular defects that promote tumour growth , destroy tumour-associated vasculature , induce anti-tumour immunity  and synergize with other treatments . OVs must overcome barriers triggered by viral infection, including those mounted by cancer cells and components of the tumour microenvironment . The type I interferon (IFN) pathway is a well-characterized mammalian signaling cascade triggered upon viral attack to protect the surrounding cells and alert the immune system to contain infection . The production of type I IFNs promotes an antiviral and anti-proliferative state in addition to inducing innate and adaptive immunity . This antiviral response represents a major barrier to virus replication and spread in healthy tissues and is necessary to the safety of OV therapy . Interestingly, the genetic alterations promoting tumourigenesis are associated with enhanced cancer cell susceptibility to viral infection . Many pathways activated in response to infection that inhibit cell growth, activate apoptosis, and alert the immune system, are incompatible with malignant growth and are often defective in cancer cells . Since these defects are common in cancer cells, they facilitate the targeted killing of cancer cells by certain OVs. Despite having IFN pathway defects, many cancers are still quite resistant to OV therapy . For example, the vesicular stomatitis virus (VSV) is an OV platform with promising potential for clinical translation . A VSV variant with an improved therapeutic index (VSV∆51) is impaired in its capacity to block the IFN response and infect normal tissues . The degree of susceptibility to VSV∆51 killing varies among human cancers , due to the IFN status of the cancer cells and the potential involvement of other antiviral mechanisms within resistant tumours.
An alternative antiviral strategy relies on RNA interference (RNAi) , in order to combat infection in plants, fungi and invertebrates. This system is similar to the microRNA (miRNA) processing pathway used for post-transcriptional regulation in most eukaryotes. Viral double-stranded RNA generated during replication and transcription is bound and cleaved by the host cytoplasmic enzyme Dicer to form 22–23 nucleotides long RNA fragments . These short RNA fragments are loaded into the RNA-induced silencing complex (RISC) where a single strand is selected to target homologous viral RNA and therefore prevent viral genome replication and translation . To counteract this RNAi-mediated antiviral response, many plant and insect viruses have evolved viral suppressors of RNAi (VSRs) . One such virus is the Nodamura virus, which primarily infects insects but is also highly virulent to certain mammals like suckling mice and hamsters [11,12,13]. Nodamura virus encodes a VSR known as B2, which binds double-stranded RNA and inhibits processing by Dicer which prevents the production of antiviral siRNAs [14,15,16].
RNAi- and protein-mediated immunity were long thought to be non-overlapping mechanisms, with insects and invertebrates using one strategy and mammals using the other. Interestingly, recent discoveries suggest that these mechanisms might not be mutually exclusive. In fact, antiviral RNAi has been shown to function in mammalian embryonic or undifferentiated cells . Considering the discovery of mammalian antiviral RNAi in embryonic stem cells and the genetic similarities between cancer cells and embryonic stem cells [18,19,20,21,22]; we hypothesized a role for antiviral RNAi in cancer cells. To investigate antiviral RNAi and its effects on OV therapy, we engineered a recombinant VSV∆51 to express the Nodamura virus B2 protein. Herein, we characterize this novel OV and demonstrate the interplay between the B2-expressing virus and RNA processing pathways in cancer. Our results show enhanced cancer-specific killing by our virus as well as improved viral replication in vivo. Together, our data strongly suggest the involvement of the RNAi pathway in the antiviral defense of cancer cells.
Evidence of a functional antiviral RNAi mechanism detected in VSV∆51 infected cancer cells
To investigate the potential involvement of antiviral RNAi mechanisms in cancer cells, we infected human cancer cell lines with a VSV variant impaired in its ability to block IFN response (VSVΔ51) and performed small-RNA deep sequencing. We showed that virus-derived small RNAs (vsRNAs) have a length bias towards 22-mers in several cell lines (Fig. 1a and Additional file 1: Figure S1A), consistent with the size of Dicer cleavage products. Importantly, this enrichment for 22-mers is present in positive strand vsRNAs, suggesting the occurrence of double stranded RNA cleavage. This is also characteristic of Dicer products, and likely acts during the synthesis of positive strands in VSV genome replication.
B2 protein enhances VSV∆51 replication in virus-resistant cancer cells
We reasoned that if antiviral RNAi is triggered in mammalian cancer cells upon VSV∆51 virus infection, expression of an RNAi viral suppressor should significantly enhance virus growth and cytotoxicity. To this end, we investigated the effects of the VSR protein B2 on VSV∆51 replication in cancer cells, and characterized two human cancer cell lines (melanoma M14 cells and renal carcinoma 786-O cells) transfected with fluorescently-tagged B2 or empty vector (mock control). M14 and 786-O cell lines were selected as models for further characterization as they are both resistant to VSV infection and have functional type I IFN pathways. Upon drug-selection and sorting of positive cells, we confirmed ectopic expression of fluorescently-tagged constructs, as shown in fluorescence microscopy images (Additional file 1: Figure S1B). Viral outputs from B2-expressing M14 and 786-O cells infected with VSV∆51 were significantly higher relative to mock controls (Fig. 1b), suggesting that B2 may enhance viral production.
VSV∆51-mediated expression of B2 enhances cytotoxicity in cancer cell lines
As ectopic expression of the B2 protein enhanced VSV∆51 titers in both M14 and 786-O cell lines, we engineered a VSV∆51 virus variant encoding His-tagged B2 to assess effects of virus-mediated B2 expression. B2, or GFP as a control, was cloned between the G and L genes of the VSV∆51 backbone (Fig. 1c), using a strategy previously shown to support expression of transgenes without impairing virus replication [7, 23]. We infected Vero cells with VSVΔ51-B2 and confirmed transgene expression by immunoblotting for His-tagged B2. As predicted, expression of B2 enhanced expression of VSV viral proteins (Fig. 1d).
While the mechanism of action of VSVΔ51-encoded B2 on mammalian cancer cells remains to be elucidated, previous studies have shown that B2 blocks processing of small RNAs by Dicer [24, 25]. Given that B2 enhanced VSV∆51 production by mammalian cancer cells, we investigated whether miRNA levels were affected by VSV∆51-B2 using quantitative PCR (qPCR) for various miRNAs from infected 786-O cells. For the majority of miRNAs tested, including miR-1, miR-16, miR-24, and miR-93, the miRNA expression levels measured in VSV∆51-B2 infected samples were significantly lower compared to VSV∆51-GFP samples (Fig. 1e), suggesting inhibition of small RNA processing by B2.
To determine if VSV∆51-B2 could kill cancer cells more efficiently than VSV∆51-GFP, we screened a panel of 38 different human cancer cell lines. The cells were infected at a multiplicity of infection (MOI) of 1 and cell viability was assessed. B2-expressing virus showed enhanced killing in the majority of cancer cell lines tested, including M14 and 786-O cells (Fig. 2a). In our screen, we also included an additional virus variant encoding vaccinia Copenhagen virus VP55, a different VSR , which similarly enhanced virus-mediated cell killing in our studies (Fig. 2a).
To determine whether VSV∆51-B2 could affect virus production, we assessed viral replication at multiple time points in M14 and 786-O cells, and found that VSV∆51-B2 significantly enhanced replication over time, relative to VSV∆51-GFP (Fig. 2b). Further, we infected GM38 fibroblasts with VSVΔ51-B2 or VSVΔ51-GFP to investigate if B2 expression affects virus replication in healthy cells, and found that VSVΔ51-B2 infection did not significantly increase viral cytotoxicity at an MOI of 1 (Additional file 1: Figure S2A).
VSV∆51-B2 prevents VSV genome cleavage in cancer cells
To determine whether B2 protects VSV from genome cleavage, we performed small-RNA deep-sequencing on VSVΔ51-B2 infected cancer cells similar to the experiment in Fig. 1. We showed that B2 expression in VSVΔ51 virus abrogates genomic cleavage as 22-mer vsRNAs are no longer prominent in different cell lines (Fig. 2c and Additional file 1: Figure. S2B). Interestingly, VSVΔ51-B2 derived vsRNAs display a bias towards positive strand reads compared to VSVΔ51 vsRNAs in cancer cells (Fig. 2c and Additional file 1: Figure S2B). Since the viral positive strand consists of viral mRNAs and positive sense genome copies, a larger bias for positive sense vRNAs suggests more efficient mRNA transcription in our VSVΔ51-B2 virus. In addition, we infected M14 and 786-O cells with VSVΔ51-GFP or VSVΔ51-B2 and showed that the expression level of virus genomes for each sample was enhanced in response to VSVΔ51-B2, relative to VSVΔ51-GFP (Fig. 2d). Taken together, these data suggest that VSV∆51-B2 inhibits direct cleavage of the viral genome and host RNA processing pathways.
VSV∆51-B2 and the type I IFN response
To characterize the effect of B2 expression on mammalian cancer cells at a transcriptome level, we performed a microarray analysis on samples from M14 cells infected with either VSV∆51-GFP or VSV∆51-B2. Our results show that at a low MOI, VSV∆51-B2 virus induced the expression of a variety of immune related genes, which were not affected by VSV∆51 infection (Fig. 3a). Through GO-term analysis, we detected an upregulation of genes by at least four-fold in response to VSV∆51-B2, but not VSV∆51-GFP. We also showed enrichment of cytokine and cytokine activity predominantly associated with an IFN response (Fig. 3b). Interestingly, at high MOI, most immune genes upregulated by VSV∆51-B2 at low MOI remain unchanged with few visible differences between the viruses (Fig. 3a).
Given that the IFN response is an important antiviral mechanism in mammalian cells, we investigated the potential impact of B2 on IFN responsiveness by qPCR after infection of 786-O cells. Consistent with our microarray results, we observed a significant increase in IFN-ß levels 24 h post-infection with VSV∆51-B2 relative to the control virus (Fig. 3c). We also showed that VSV∆51-B2 enhances IFN-ß secretion in 786-O cells, relative to VSV∆51-GFP, by ELISA (Fig. 3d). Lastly, we investigated if the production of VSV∆51-B2 could be further enhanced by blocking the IFN pathway. We have previously shown that B19R, a soluble type I IFN scavenger expressed by vaccinia Copenhagen virus, enhances VSV∆51 production . In order to block the antiviral effect of the IFN that would be produced in response to VSV∆51 infection, we generated conditioned-media from vaccinia virus-infected HeLa cells and pre-treated 786-O cells with B19R-containing media. We found that viral titers were significantly higher for VSV∆51-B2 with both control media and vaccinia Copenhagen virus conditioned-media compared to VSV∆51-GFP (Fig. 3e); however, the absolute increase in virus titers after exposing the cells to vaccinia Copenhagen virus-conditioned media was similar for both VSVΔ51-GFP and VSVΔ51-B2. Although induction of antiviral RNAi has been shown in mature mammalian cells [27, 28], our data suggest that the RNAi-based pathway is an IFN-independent antiviral mechanism in cancer cells as B2 expression alone enhances viral titers, which can be further enhanced by blocking the IFN response.
Virus-mediated B2 expression enhances replication and cytokine production in vivo
To establish an in vivo mouse model, we first screened RENCA mouse renal carcinoma cells in vitro to determine whether VSV∆51-B2 enhanced cytotoxicity. We showed that, as observed using human cell lines, RENCA cells were more efficiently killed by the B2-expressing virus (Additional file 1: Figure S2C). In vivo, we subsequently tested the RENCA cell line which is syngeneic to Balb/c mice. In addition, we used the human M14 melanoma cell line as a xenograft model in nude mice. For both models, VSV∆51-B2 titers from subcutaneous tumours harvested 24 h post-intratumoural viral injection were significantly higher compared to VSV∆51-GFP (Fig. 4a). Consistent with our microarray analysis, the IFN-γ, TNF-α and MCP-1 concentrations from the serum of the RENCA tumour-bearing Balb/c mice were significantly increased for the VSV∆51-B2 treated mice compared to VSV∆51-GFP treated mice (Fig. 4b). In contrast, Il-6 concentrations did not significantly increase (Fig. 4b). Biodistribution analyses following intravenous administration revealed unchanged amounts of virus obtained from the different organs for VSV∆51 and VSV∆51-B2 (Additional file 1: Figure S2D-E).
In this study, we demonstrate that B2 expression is sufficient to enhance VSV∆51 replication and cytotoxicity (Fig. 1-2) in mammalian cancer cells. Our screen of human cancer cells demonstrated enhanced efficacy of VSV∆51-B2 in killing the majority of cancer cell lines tested. The enhanced cytotoxic ability of VSV∆51-B2 suggests that RNAi may be an important factor hindering virus replication in resistant cancer cells. It is important to note that the cell lines in which there was no difference in cytotoxicity were the most sensitive to viral infection (Fig. 2), suggesting that the lack of improvement could potentially be attributed to the already maximal virus production by these cells. An antiviral RNAi system may still function in these cells but may only be apparent in certain conditions, such as within the tumour microenvironment where many factors come together and create additional barriers to infection.
Interestingly, virus variants encoding B2 and VP55, two VSRs that impair the RNAi response by different mechanisms, show the same improvement in killing ability for all cell lines tested (Fig. 2). The mechanism of B2 involves binding of small RNA fragments which could either prevent their processing by Dicer or loading into RISC. On the other hand, VP55 polyadenylates miRNAs which targets them for degradation . Given that both VSRs improve VSV∆51-mediated killing to the same extent, this suggests that inhibition of the RNAi pathway improves virus replication regardless of the mechanism through which inhibition is achieved. VP55 does not polyadenylate all small RNAs and key features such as the presence of a 2’O methyl group protect a subset of small RNAs from degradation . Notably, vsRNAs have been previously shown to be 2’O methylated, which protects them from degradation . This may be advantageous for cell lines in which direct viral genome cleavage occurs as ideally, cleaved genome fragments can provide additional protection through targeting homologous viral genomes and transcripts.
To begin exploring the mechanism of B2 on VSV∆51 replication, we investigated the potential impact of B2 on the IFN response. A number of cellular proteins such as Toll-like receptor 3, retinoic acid inducible gene I, 2′-5′ oligoadenylate synthetase and protein kinase R recognize dsRNA and trigger a potent antiviral immune response . Therefore, B2 may sequester the dsRNA substrates of these antiviral factors or interact with these proteins to prevent the sensing of dsRNA. As such, we investigated the effects of B2 on the IFN response, which is downstream of these pathways. We demonstrate that B2 does not suppress but actually significantly increases IFN-ß production (Fig. 3) compared to control virus, which is likely a result of enhanced replication. More specifically, at low MOI we detect an upregulation of immune genes with VSV∆51-B2 due to the ability of the VSV∆51-B2 virus to replicate in M14 cells and establish a successful infection which triggers a more robust IFN response (Fig. 3). However, at high MOI we do not see upregulation of immune genes by either VSV∆51-B2 or VSV∆51-GFP virus (Fig. 3). Higher MOI is often used to overcome resistance to infectivity and often leads to faster cell death, suggesting there was insufficient time to mount a type I IFN response as most, if not all, cells were infected in the first round of replication. Of note, low MOI levels are more comparable to in vivo systems, in which VSV∆51-B2 virus may be more immunogenic. This prediction is in line with higher cytokine (IFN-γ, TNF-α and MCP-1) levels found in the serum of tumour bearing VSV∆51-B2 treated immune-competent mice (Fig. 4).
Additionally, our vaccinia Copenhagen virus-conditioned media transfer experiments demonstrated that the blocking of type I IFN did not further enhance VSV∆51-B2 replication compared to VSV∆51-GFP. This is an indirect way to neutralize IFN-1  and further suggested that IFN and B2 have different mechanisms of action (Fig. 3). Consistent with this notion, a similar induction of IFN-stimulated genes was observed in wildtype and RNAi-defective mouse embryonic fibroblasts . In addition, suppression of RNAi by Nodamura virus protein B2 protein does not change the expression levels of IFN-stimulated genes in infected mice . Importantly, our data do not eliminate the possibility of IFN-stimulated miRNAs limiting VSV∆51 efficacy. Despite stimulating a type I IFN response, VSV∆51-B2 has enhanced replication, suggesting that B2 expression is sufficient to overcome the antiviral response it stimulates.
One study has shown production of abundant influenza viral siRNAs in IFN-competent A549 cancer cells , but the existence of antiviral RNAi in cancer cells remains largely unexplored. We were able to detect viral genome fragments (vsRNAs) following infection with VSV∆51 in a number of cell lines (Fig. 1 and Additional file 1: Figure S1), suggesting that viral genome cleavage is occurring. This may be facilitated by RNAi machinery being repurposed during viral infection. Interestingly, VSV∆51-B2 infection resulted in a decrease in the percentage of vsRNAs (Fig. 2 and Additional file 1: Figure S2), suggesting the efficient prevention of this antiviral mechanism. A number of recent reports support the notion of antiviral genome cleavage in mammalian cells. It has been demonstrated that infection with RNA viruses can trigger the production of viral siRNAs, presumably as a result of direct virus genome cleavage [17, 32]. The fact that many mammalian viruses encode VSRs, further supports the concept of a mammalian RNAi system. For example, Influenza A encodes the NS1 protein [33, 34], Ebola encodes VP35 [35, 36], HIV-1 encodes Tat [37, 38] and vaccinia Copenhagen virus encodes VP55 . These proteins all have VSR-like functions, which suggests an evolutionary advantage of blocking antiviral RNAi.
It is possible that both direct inhibition of viral genome cleavage and inhibition of cellular antiviral miRNA production are co-occurring within certain cancer cell lines, as we observe changes in both VSV-specific and total RNA read length distributions (Additional file 1: Figure S1 and S2 and data not shown). The potential interplay between viruses and infected host cell miRNAs is a concept that is supported by several studies. For example, miR-29 has been reported to bind the 3’ UTR of HIV mRNA which inhibits its translation and results in sequestering the mRNA into processing bodies . IFN-ß itself induces transcription of a number of miRNAs in hepatocytes that are complimentary to the hepatitis C virus genomic RNA and inhibit viral replication . In fact, activation of the IFN pathway has been shown to lead to the upregulation of a number of miRNAs including miR-1 [40,41,42], miR-129 , miR-146  and miR-155 [42, 44, 45], some of which likely function to control infection. Perhaps most relevant to our study, 2 miRNAs (miR-24 and miR-93) have been previously shown to directly target the VSV genome and limit VSV replication . Our results show that VSV∆51-B2 infection leads to the downregulation of both of these miRNAs, providing a potential explanation for the increased virus production using the B2 virus.
Overall, we demonstrate a novel role of the RNAi pathway as an intrinsic antiviral mechanism in cancer cells and how RNAi inhibition may be used to improve OV replication. Mechanistically, inhibition of direct viral genome cleavage and/or modulation of miRNA processing contribute to enhancing VSV∆51 infection in a cell line specific manner. This work provides insight into the basic biology of viral defense mechanisms in cancer and promises to improve current OV therapies by tailoring viruses to overcome alternative antiviral mechanisms.
Cell lines and culture
All cell lines were purchased from the American Type Culture Collections (Manassas, VA). Mammalian cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM) (Corning cellgro, Manassas, VA) or RPMI-1640 (K562, OVCAR3, OVCAR4, OVCAR8) (Corning cellgro, Manassas, VA) supplemented with 10% fetal bovine serum (FBS) (Sigma life science, St-Louis, MO) and maintained at 37°C with 5% CO2. Drosophila melanogaster Schneider line 2 (S2) cells were cultured in SF900II serum-free medium (Invitrogen) at 25 °C under atmospheric pressures.
DNA constructs and viral constructs
The pcDNA3.1-puro B2 (Nodamura gene) plasmid used for generating the B2-expressing stable cell lines was made available by the Christopher Sullivan lab (Addgene plasmid # 17228). The pEGFP-N1 plasmid (Catalog # 6085–1) was purchased from Clontech (Moutain View, CA).
The B2 gene was amplified by PCR of the Nodamura virus genome. The primers were designed to include the XhoI and NheI restriction sites as well as to insert a 6× histidine tag at the 5′ end of the B2 sequence. The digested PCR fragment was cloned into the XhoI and NheI digested VSVΔ51 backbone, as previously described. The primer pairs for inserting B2 into the VSV backbone are listed in Additional file 2: Table S1.
VP55 was PCR amplified from the Copenhagen strain of vaccinia virus and subcloned into pcDNA3.1 with an N-terminal Flag epitope. Flag-VP55 was subsequently PCR amplified and cloned into VSVΔ51M using the same strategy.
Transfection and selection of cell lines
M14 and 786-O cells were transfected using Lipofectamine 2000 (Invitrogen, Carlsbad, CA) according to manufacturer instructions. Briefly, cells were plated in 6 well format 1 day prior to transfection. Plasmid and Lipofectamine reagent were incubated for 20 min and then added to plated cells in OptiMEM (Thermoscientific, Waltham, MA). 24 h post-transfection, medium was replaced by DMEM with 10% FBS and cultured for 48 h. Cells were then subjected to drug selection by the addition of Geneticin (800 μg/mL) (Thermo Fisher, Carlsbad, CA). Cells were expanded and GFP- or YFP-positive cells were sorted twice by FACS (MoFlo Astrios).
Viral titers were obtained by plaque assays. Serial dilutions of the samples were prepared in serum-free DMEM. The dilutions were then transferred to monolayers of Vero cells and incubated at 37 °C for 1 h. After the incubation, cells were overlaid with 0.5% agarose in DMEM supplemented with 10% FBS. Plates were incubated for 24 h at 37 °C with 5% CO2 and plaques were counted.
Virus rescue and purification
Virus rescues were performed as previously described. Vero cells were infected with T7 polymerase-expressing vaccinia Copenhagen virus at an MOI of 3. Following a 2 h incubation, media was removed and cells were transfected with T7-driven plasmids encoding VSV N, P, and L genes as well as the VSV∆51-B2 backbone. Supernatants collected 48 h post-transfection were passed through a 0.22 μm filter (MillexGP, Carrigtwohill, Ireland) to remove vaccinia Copenhagen virus.
For expansion and purification of the viral preparations, Vero cells were infected at an MOI of 0.001 and culture supernatants were collected 24 h post-infection. Supernatants were then filtered through a 0.2 μm bottle top filter (Millipore, Etobicoke, Canada) and centrifuged at 30,100 g for 90 min. The supernatant was discarded and pelleted virus was resuspended in Dulbecco’s phosphate-buffered saline (Corning cellgro, Manassas, VA). Purified virus was kept at − 80°C.
Deep sequencing of vsRNAs
Total RNA was extracted with TRIzol reagent (Invitrogen) according to manufacturer instructions. Library preparation for Illumina sequencing was performed (TCAG DNA Sequencing Facility, Toronto, ON). Briefly, RNA was enriched for 15–25 nt sizes before strand-specific, small-RNA library preparation and 50 bp single end read sequencing. Adapter trimming was done with Trimmomatic  following default parameters. Before read mapping, VSV∆51 genome was constructed from the VSV reference genome (NC_001560), manually edited to delete the 51st methionine amino acid in the M gene. Reads were mapped to VSV∆51 genome using bbmap.sh script from the BBMap toolkit (http://sourceforge.net/projects/bbmap) with a minimum alignment identity of 100%. SAMtools was used to separate positive sense mapping from bbmap produced sam files . Lastly, positive sense reads were analyzed for size distribution using the reformat.sh script from the BBMap toolkit.
Cell pellets were lysed on ice for 30 min using complete protease inhibitor cocktail (Roche, Mississauga, Ontario, Canada) supplemented lysis buffer (1% NP40, 150 mM NaCl, 5 mM EDTA, 50 mM Tris pH 7.4). Lysates were centrifuged for 10 min at 16,000 g and cleared supernatants were mixed with dithiothreitol-supplemented loading buffer (250 mM Tris-HCl pH 6.8, 10% SDS, 30% glycerol, 5% β-Mercaptoethanol, 0.02% bromophenol blue). The samples were migrated on Bio-Rad Mini Protean 4–15% TGX Protein Gels (Bio-Rad, Mississauga, ON) and transferred onto PVDF membranes (GE Healthcare, Buckinghanshire, UK) prior to blocking with 5% skim milk powder (Oxoid Ltd., Basingstoke, UK) in Tris-buffered saline (TBS)with 0.1% Tween-20. The membranes were proofed using specific rabbit antibodies for 6× His tag (Abcam, Cambridge, UK), VSV (polyclonal anti-VSV serum for hyperimmune rabbits) . Rabbit anti-GAPDH (Abcam, Cambridge, UK) and rat anti-tubulin (Novus Biologicals, Oakville, ON) antibodies were used as loading controls. Membranes were then probed with a horse radish peroxidase-coupled goat anti-rabbit secondary antibody (Millipore, Etobicoke, Canada) or goat anti-rat secondary antibody (Life Technologies, Carlsbad, CA) and the signal was revealed using Amersham ECL Western Blotting Detection Reagent (GE Healthcare, Buckinghamshire, UK). The gels were analyzed using FluorChem FC2 (Alpha Innotech, San Leandro, CA).
Cell viability assay
Relative metabolic activity of cells was used as a readout of cell viability and was determined using alamarBlue reagent (Bio-Rad, Mississauga, Ontario, Canada). The assays were performed according to manufacturer instructions. Briefly, cells were plated in 96-well plates and infected with the different viruses 24 h later. 48 h after virus infection, alamarBlue was added to each well to a final concentration of 1:10. The samples were incubated 1 to 5 h and the fluorescence readings (excitation and emission wavelengths of 530 nm and 590 nm, respectively) were taken using a Fluoroskan Ascent FL (Thermo Labsystems, Beverly, MA).
For miRNA qPCRs, RNA was extracted from infected cell pellets using TRIzol reagent (Life Technologies, Carlsbad, CA) according to the manufacturer instructions. The RNA concentration and purity was assessed using a NanoDrop ND-1000 spectrophotometer (Thermoscientific, Waltham, MA) prior to reverse transcription using Quanta miRNA cDNA synthesis kit (Gaithersburg, MD).
For all other qPCRs, RNA was extracted using RNAeasy RNA extraction kit (QIAGEN, Toronto, ON, Canada) according to the manufacturer’s instructions. The RNA concentration and purity was assessed using a NanoDrop ND-1000 spectrophotometer (Thermoscientific, Waltham, MA) prior to reverse transcription using RevertAid H Minus First Strand cDNA Synthesis kit (Thermoscientific, Waltham, MA).
qRT -PCR was performed on the non-pooled triplicate samples. Following conversion to cDNA by Superscript RT II (Invitrogen, Carlsbad, CA), qRT-PCR was carried out using Sybergreen (Invitrogen) according to the manufacturer’s instructions. Analyses were performed on a Rotor-Gene RG-3000A machine (Corbett Research, Mortlake, AU) according to the manufacture instructions. The primer pairs specific for various gene products used in our experiments are listed in Additional file 2: Table S1. qRT-PCR measurements were normalized to the U6 or GAPDH house-keeping genes for miRNA or RNA transcripts, respectively, using the Pfaffl method .
Monolayers of M14 cells were treated at an MOI of 0.1 or 3 for 24 h with either VSV∆51 or VSV∆51 encoding B2 gene virus Total RNA was extracted with TRIzol reagent (Invitrogen) according to the manufacturer instructions. Experimental triplicate total RNA samples were processed by The Center for Applied Genomics at The Hospital for Sick Children for microarray analysis on a Human Prime View chip. Raw files were analyzed using the Transcriptome Analysis Console v3.0 (Affymetrix) software. Normalized transcript expression values further processed with R. Heatmaps were produced using the R package “pheatmap” v1.0.8. GO Term Enrichment analysis (http://CRAN.R-project.org/package=pheatmap) was performed using the online EnrichR tool . Genes selected for enrichment analysis are the subset of genes upregulated by the expression of the B2 gene in VSV∆51 by at least 4-fold.
The concentration of IFN-β was determined using the human IFN- β ELISA kit (R&D Systems, Minneapolis, MN) according to the manufacturer’s instructions.
Supernatant transfer experiments
Vaccinia stocks were propagated in U-2 OS cells and cell-associated virus was collected by repeat  freeze–thaw cycles. Purification of viral stocks was done by centrifugation at 20,700 g through a 36% sucrose cushion (in 1 mM Tris) before resuspension in 1 mM Tris, pH 9.
To generate infected cell conditioned media, U-2 OS cells were either mock treated or infected with VVdd-mCherry at a multiplicity of 10 PFU/cell for 24 h, harvested and then pelleted by centrifugation. Supernatants were collected and passed through a 0.22 μm filter to eliminate cell-free vaccinia virions. To test for factors enhancing VSV infectivity, tumour cell monolayers were pre-treated for 2 h with conditioned U-2 OS supernatant. Tumour cells were then infected with VSV in the presence of conditioned medium.
In vivo experiments and tumour models
6–8 weeks old female Balb/c or nude mice (Charles River Laboratories, Wilmington, MA) were used. For the Balb/c mice 5 × 105 RENCA tumour cells were implanted subcutaneously 21 days prior to treatment. For the nude mice 1 × 108 M14 tumour cells were implanted subcutaneously 14 days prior to treatment. A single intratumoural injection of 1E8 PFU of VSVΔ51-GFP or VSVΔ51-B2 was performed. Tumours were harvested 24 h post-treatment, weighed and homogenized in PBS using a Powergen 125 tissue homogenizer (Fischer Scientific, Hampton, New Hampshire). Serial dilutions of tumour homogenates were tittered to obtain intratumoural titer data. For biodistribution experiments, naïve Balb/c mice were treated with a single intravenous injection of 1E8 PFU of VSVΔ51-GFP or VSVΔ51-B2. Mice were sacrificed 24 h or 48 h post-treatment and organs were harvested, frozen and homogenized in PBS as above. All experiments were approved by the University of Ottawa animal care and veterinary services (ME-2258).
Serum was obtained upon centrifugation of blood samples collected using lithium-heparin coated capillary tubes (Sarstedt, Newton) at 16000 g for 5 min. The concentrations of IFNγ, TNFα, MCP-1 and IL-6 were measured using a mouse inflammation cytometric bead array kit (BD Biosciences, San Jose, CA) according to manufacturer instructions. The results were acquired on a LSR Fortessa flow cytometer and analyzed using the FCAP array software (BD Biosciences).
RNA-induced silencing complex
Viral suppressors of RNAi
Virus-derived small RNAs
Vesicular stomatitis virus
Lichty BD, Breitbach CJ, Stojdl DF, Bell JC. Going viral with cancer immunotherapy. Nat Rev Cancer. 2014;14(8):559–67.
Ilkow CS, Swift SL, Bell JC, Diallo JS. From scourge to cure: tumour-selective viral pathogenesis as a new strategy against cancer. PLoS Pathog. 2014;10(1):e1003836.
Arulanandam R, Batenchuk C, Angarita FA, Ottolino-Perry K, Cousineau S, Mottashed A, Burgess E, Falls TJ, De Silva N, Tsang J, Howe GA, Bourgeois-Daigneault MC, Conrad DP, Daneshmand M, Breitbach CJ, et al. VEGF-mediated induction of PRD1-BF1/Blimp1 expression sensitizes tumor vasculature to oncolytic virus infection. Cancer Cell. 2015;28(2):210–24.
Bourgeois-Daigneault MC, Roy DG, Aitken AS, El Sayes N, Martin NT, Varette O, Falls T, St-Germain LE, Pelin A, Lichty BD, Stojdl DF, Ungerechts G, Diallo JS, Bell JC. Neoadjuvant oncolytic virotherapy before surgery sensitizes triple-negative breast cancer to immune checkpoint therapy. Sci Transl Med. 2018;10(422)
Ilkow CS, Marguerie M, Batenchuk C, Mayer J, Ben Neriah D, Cousineau S, Falls T, Jennings VA, Boileau M, Bellamy D, Bastin D, de Souza CT, Alkayyal A, Zhang J, Le Boeuf F, et al. Reciprocal cellular cross-talk within the tumor microenvironment promotes oncolytic virus activity. Nat Med. 2015;21(5):530–6.
Ivashkiv LB, Donlin LT. Regulation of type I interferon responses. Nat Rev Immunol. 2014;14(1):36–49.
Stojdl DF, Lichty BD, tenOever BR, Paterson JM, Power AT, Knowles S, Marius R, Reynard J, Poliquin L, Atkins H, Brown EG, Durbin RK, Durbin JE, Hiscott J, Bell JC. VSV strains with defects in their ability to shutdown innate immunity are potent systemic anti-cancer agents. Cancer Cell. 2003;4(4):263–75.
Stojdl DF, Lichty B, Knowles S, Marius R, Atkins H, Sonenberg N, Bell JC. Exploiting tumor-specific defects in the interferon pathway with a previously unknown oncolytic virus. Nat Med. 2000;6(7):821–5.
Wilson RC, Doudna JA. Molecular mechanisms of RNA interference. Annu Rev Biophys. 2013;42:217–39.
Yan N, Chen ZJ. Intrinsic antiviral immunity. Nat Immunol. 2012;13(3):214–22.
Jiang L, Wei C, Li Y. Viral suppression of RNA silencing. Sci China Life Sci. 2012;55(2):109–18.
Bailey L, Newman JF, Porterfield JS. The multiplication of Nodamura virus in insect and mammalian cell cultures. J Gen Virol. 1975;26(1):15–20.
Bailey L, Scott HA. The pathogenicity of Nodamura virus for insects. Nature. 1973;241(5391):545.
Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, Sweet-Cordero A, Ebert BL, Mak RH, Ferrando AA, Downing JR, Jacks T, Horvitz HR, Golub TR. MicroRNA expression profiles classify human cancers. Nature. 2005;435(7043):834–8.
Sullivan CS, Ganem D. A virus-encoded inhibitor that blocks RNA interference in mammalian cells. J Virol. 2005;79(12):7371–9.
Aliyari R, Wu Q, Li HW, Wang XH, Li F, Green LD, Han CS, Li WX, Ding SW. Mechanism of induction and suppression of antiviral immunity directed by virus-derived small RNAs in Drosophila. Cell Host Microbe. 2008;4(4):387–97.
Li Y, Lu J, Han Y, Fan X, Ding SW. RNA interference functions as an antiviral immunity mechanism in mammals. Science. 2013;342(6155):231–4.
Boumahdi S, Driessens G, Lapouge G, Rorive S, Nassar D, Le Mercier M, Delatte B, Caauwe A, Lenglez S, Nkusi E, Brohee S, Salmon I, Dubois C, del Marmol V, Fuks F, et al. SOX2 controls tumour initiation and cancer stem-cell functions in squamous-cell carcinoma. Nature. 2014;511(7508):246–50.
Ben-Porath I, Thomson MW, Carey VJ, Ge R, Bell GW, Regev A, Weinberg RA. An embryonic stem cell-like gene expression signature in poorly differentiated aggressive human tumors. Nat Genet. 2008;40(5):499–507.
Somervaille TC, Matheny CJ, Spencer GJ, Iwasaki M, Rinn JL, Witten DM, Chang HY, Shurtleff SA, Downing JR, Cleary ML. Hierarchical maintenance of MLL myeloid leukemia stem cells employs a transcriptional program shared with embryonic rather than adult stem cells. Cell Stem Cell. 2009;4(2):129–40.
Schoenhals M, Kassambara A, De Vos J, Hose D, Moreaux J, Klein B. Embryonic stem cell markers expression in cancers. Biochem Biophys Res Commun. 2009;383(2):157–62.
Mizuno H, Spike BT, Wahl GM, Levine AJ. Inactivation of p53 in breast cancers correlates with stem cell transcriptional signatures. Proc Natl Acad Sci U S A. 2010;107(52):22745–50.
Bourgeois-Daigneault MC, Roy DG, Falls T, Twumasi-Boateng K, St-Germain LE, Marguerie M, Garcia V, Selman M, Jennings VA, Pettigrew J, Amos S, Diallo JS, Nelson B, Bell JC. Oncolytic vesicular stomatitis virus expressing interferon-gamma has enhanced therapeutic activity. Mol Ther Oncolytics. 2016;3:16001.
Chao JA, Lee JH, Chapados BR, Debler EW, Schneemann A, Williamson JR. Dual modes of RNA-silencing suppression by flock house virus protein B2. Nat Struct Mol Biol. 2005;12(11):952–7.
Johnson KL, Price BD, Eckerle LD, Ball LA. Nodamura virus nonstructural protein B2 can enhance viral RNA accumulation in both mammalian and insect cells. J Virol. 2004;78(12):6698–704.
Le Boeuf F, Diallo JS, McCart JA, Thorne S, Falls T, Stanford M, Kanji F, Auer R, Brown CW, Lichty BD, Parato K, Atkins H, Kirn D, Bell JC. Synergistic interaction between oncolytic viruses augments tumor killing. Mol Ther. 2010;18(5):888–95.
Li Y, Basavappa M, Lu J, Dong S, Cronkite DA, Prior JT, Reinecker HC, Hertzog P, Han Y, Li WX, Cheloufi S, Karginov FV, Ding SW, Jeffrey KL. Induction and suppression of antiviral RNA interference by influenza a virus in mammalian cells. Nat Microbiol. 2016;2:16250.
Qiu Y, Xu Y, Zhang Y, Zhou H, Deng YQ, Li XF, Miao M, Zhang Q, Zhong B, Hu Y, Zhang FC, Wu L, Qin CF, Zhou X. Human virus-derived small RNAs can confer antiviral immunity in mammals. Immunity. 2017;46(6):992–1004. e1005
Backes S, Shapiro JS, Sabin LR, Pham AM, Reyes I, Moss B, Cherry S, tenOever BR. Degradation of host microRNAs by poxvirus poly(a) polymerase reveals terminal RNA methylation as a protective antiviral mechanism. Cell Host Microbe. 2012;12(2):200–10.
Daffis S, Szretter KJ, Schriewer J, Li J, Youn S, Errett J, Lin TY, Schneller S, Zust R, Dong H, Thiel V, Sen GC, Fensterl V, Klimstra WB, Pierson TC, et al. 2'-O methylation of the viral mRNA cap evades host restriction by IFIT family members. Nature. 2010;468(7322):452–6.
Gantier MP, Williams BR. The response of mammalian cells to double-stranded RNA. Cytokine Growth Factor Rev. 2007;18(5–6):363–71.
Maillard PV, Ciaudo C, Marchais A, Li Y, Jay F, Ding SW, Voinnet O. Antiviral RNA interference in mammalian cells. Science. 2013;342(6155):235–8.
Lu Y, Wambach M, Katze MG, Krug RM. Binding of the influenza virus NS1 protein to double-stranded RNA inhibits the activation of the protein kinase that phosphorylates the elF-2 translation initiation factor. Virology. 1995;214(1):222–8.
Bucher E, Hemmes H, de Haan P, Goldbach R, Prins M. The influenza a virus NS1 protein binds small interfering RNAs and suppresses RNA silencing in plants. J Gen Virol. 2004;85(Pt 4):983–91.
Cardenas WB, Loo YM, Gale M Jr, Hartman AL, Kimberlin CR, Martinez-Sobrido L, Saphire EO, Basler CF. Ebola virus VP35 protein binds double-stranded RNA and inhibits alpha/beta interferon production induced by RIG-I signaling. J Virol. 2006;80(11):5168–78.
Haasnoot J, de Vries W, Geutjes EJ, Prins M, de Haan P, Berkhout B. The Ebola virus VP35 protein is a suppressor of RNA silencing. PLoS Pathog. 2007;3(6):e86.
Bennasser Y, Le SY, Benkirane M, Jeang KT. Evidence that HIV-1 encodes an siRNA and a suppressor of RNA silencing. Immunity. 2005;22(5):607–19.
Bennasser Y, Jeang KT. HIV-1 tat interaction with dicer: requirement for RNA. Retrovirology. 2006;3:95.
Nathans R, Chu CY, Serquina AK, Lu CC, Cao H, Rana TM. Cellular microRNA and P bodies modulate host-HIV-1 interactions. Mol Cell. 2009;34(6):696–709.
Pedersen IM, Cheng G, Wieland S, Volinia S, Croce CM, Chisari FV, David M. Interferon modulation of cellular microRNAs as an antiviral mechanism. Nature. 2007;449(7164):919–22.
Scagnolari C, Zingariello P, Vecchiet J, Selvaggi C, Racciatti D, Taliani G, Riva E, Pizzigallo E, Antonelli G. Differential expression of interferon-induced microRNAs in patients with chronic hepatitis C virus infection treated with pegylated interferon alpha. Virol J. 2010;7:311.
Bruni R, Marcantonio C, Tritarelli E, Tataseo P, Stellacci E, Costantino A, Villano U, Battistini A, Ciccaglione AR. An integrated approach identifies IFN-regulated microRNAs and targeted mRNAs modulated by different HCV replicon clones. BMC Genomics. 2011;12:485.
Zhang J, Li S, Yan Q, Chen X, Yang Y, Liu X, Wan X. Interferon-beta induced microRNA-129-5p down-regulates HPV-18 E6 and E7 viral gene expression by targeting SP1 in cervical cancer cells. PLoS One. 2013;8(12):e81366.
Zhang J, Zhao H, Chen J, Xia B, Jin Y, Wei W, Shen J, Huang Y. Interferon-beta-induced miR-155 inhibits osteoclast differentiation by targeting SOCS1 and MITF. FEBS Lett. 2012;586(19):3255–62.
O'Connell RM, Taganov KD, Boldin MP, Cheng G, Baltimore D. MicroRNA-155 is induced during the macrophage inflammatory response. Proc Natl Acad Sci U S A. 2007;104(5):1604–9.
Otsuka M, Jing Q, Georgel P, New L, Chen J, Mols J, Kang YJ, Jiang Z, Du X, Cook R, Das SC, Pattnaik AK, Beutler B, Han J. Hypersusceptibility to vesicular stomatitis virus infection in Dicer1-deficient mice is due to impaired miR24 and miR93 expression. Immunity. 2007;27(1):123–34.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. Genome project data processing S. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
Roy DG, Power AT, Bourgeois-Daigneault MC, Falls T, Ferreira L, Stern A, Tanese de Souza C, McCart JA, Stojdl DF, Lichty BD, Atkins H, Auer RC, Bell JC, Le Boeuf F. Programmable insect cell carriers for systemic delivery of integrated cancer biotherapy. J Control Release. 2015;220(Pt A):210–21.
Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29(9):e45.
Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, Koplev S, Jenkins SL, Jagodnik KM, Lachmann A, McDermott MG, Monteiro CD, Gundersen GW, Ma'ayan A. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016;44(W1):W90–7.
We would like to thank the exceptional technical support of Catia Cemeus, Christiano Tanese de Souza and Julia Petryk, as well as members of the Bell, Auer, Atkins, Ilkow and Diallo laboratories for feedback on this project.
This work was funded by grants from the Terry Fox Research Foundation and the Canadian Institutes of Health Research to J.C.B. and C.S.I. J.C.B. is also supported by the Ontario Institute for Cancer Research and the Ottawa Regional Cancer Foundation. A.A. is funded by the Canadian Institutes of Health Research Frederick Banting and Charles Best Master’s Award. A.P. is funded by the Joseph and Wolf Lebovic Cancer Genomics and Immunity Fellowship and Ontario Graduate Scholarship. L.P. and M.C.B.D. are funded by the Canadian Institutes of Health Research post-doctoral Fellowship.
Availability of data and materials
Deep-sequencing datasets of VSV genome reads generated during the current study are available in the Sequence Read Archive (SRA) repository (Accession number: SRP132124). Microarray data are deposited in the NCBI Gene Expression Omnibus (Accession number: GSE115138).
All animal studies were approved by the institutional animal care committee of the University of Ottawa (Protocol ID: ME-2258) and carried out in accordance with guidelines of the National Institutes of Health and the Canadian Council on animal care.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. B2 selectively enhances VSV∆51 replication in other cancer cells. (A) MCF7, HT-29, or SF-295 cells were infected with VSVΔ51 virus and small-RNA deep sequencing was performed. Virus-derived small RNAs have a length bias towards 22-mers. The enrichment for 22-mers is indicated for positive strand vsRNAs. (B) Fluorescence microscopy images of M14 or 786-O cells stably expressing EGFP-B2 or fluorescently-tagged empty vector (mock control). Figure S2. VSVΔ51-B2 does not enhance viral replication in non-cancer healthy cells or alter biodistribution in various organs. (A) Relative metabolic activity of GM38 fibroblasts infected with VSVΔ51-GFP or VSVΔ51-B2 for 48 h at an MOI of 1. The results are expressed as a percentage of the signal obtained compared to mock treatment. NS: P > 0.1, *P < 0.1, **P < 0.01, ***P < 0.001, using Student’s t-test. Only significantly different pairs are indicated on the fig. (B) We performed small-RNA deep-sequencing using MCF7, HT-29, or SF-295 cells infected with VSVΔ51-B2 at an MOI of 0.1 for 18 h. B2 expression in VSVΔ51 virus abrogates genomic cleavage as 22-mer vsRNAs are no longer prominent. VSVΔ51-B2 derived vsRNAs display a bias towards positive strand reads in M14 and 786-O cells. (C) Relative metabolic activity of RENCA cells infected with VSVΔ51-GFP or VSVΔ51-B2 for 48 h at an MOI of 1. The results are expressed as a percentage of the signal obtained compared to mock treatment. (D&E) Biodistribution of VSVΔ51-B2 in tumour-bearing C57BL/6 mice. Viral titers obtained from organs of tumour-bearing C57BL/6 mice, D] 24 or E] 48 hpi. Virus was administered intravenously at a dose of 1E9 pfu of VSVΔ51-GFP or VSVΔ51-B2. For organs where virus was undetectable, the titer was considered to be the value of the limit of detection of titering for this assay (5E1 pfu/organ). NS: P > 0.1, *P < 0.1, **P < 0.01, ***P < 0.001, using Student’s t-test. (PDF 5336 kb)
Table S1. List of primers used in qRT-PCR assays. (DOCX 14 kb)