VODKA2: a fast and accurate method to detect non-standard viral genomes from large RNA-seq data sets.

During viral replication, viruses carrying an RNA genome produce non-standard viral genomes (nsVGs), including copy-back viral genomes (cbVGs) and deletion viral genomes (delVGs), that play a crucial role in regulating viral replication and pathogenesis. Because of their critical roles in determining the outcome of RNA virus infections, the study of nsVGs has flourished in recent years, exposing a need for bioinformatic tools that can accurately identify them within next-generation sequencing data obtained from infected samples. Here, we present our data analysis pipeline, Viral Opensource DVG Key Algorithm 2 (VODKA2), that is optimized to run on a parallel computing environment for fast and accurate detection of nsVGs from large data sets.

VODKA2: A fast and accurate method to detect non-standard viral genomes from large RNA-seq datasets.

During viral replication, viruses carrying an RNA genome produce non-standard viral genomes (nsVGs), including copy-back viral genomes (cbVGs) and deletion viral genomes (delVGs), that play a crucial role in regulating viral replication and pathogenesis. Because of their critical roles in determining the outcome of RNA virus infections, the study of nsVGs has flourished in recent years exposing a need for bioinformatic tools that can accurately identify them within Next-Generation Sequencing data obtained from infected samples. Here, we present our data analysis pipeline, Viral Opensource DVG Key Algorithm2 (VODKA2), that is optimized to run on a High Performance Computing (HPC) environment for fast and accurate detection of nsVGs from large data sets.

Accumulation of copy-back viral genomes during respiratory syncytial virus infection is preceded by diversification of the copy-back viral genome population followed by selection.

RNA viruses generate nonstandard viral genomes during their replication, including viral genomes of the copy-back (cbVGs) type that cannot replicate in the absence of a standard virus. cbVGs play a crucial role in shaping virus infection outcomes due to their ability to interfere with virus replication and induce strong immune responses. However, despite their critical role during infection, the principles that drive the selection and evolution of cbVGs within a virus population are poorly understood. As cbVGs are dependent on the virus replication machinery to be generated and replicated, we hypothesized that host factors that affect virus replication exert selective pressure on cbVGs and drive their evolution within a virus population. To test this hypothesis, we used respiratory syncytial virus (RSV) as a model and took an experimental evolution approach by serially passaging RSV in immune-competent human lung adenocarcinoma A549 control and immune-deficient A549 Signal transducer and activator of transcription 1 (STAT1) KO cells, which allow higher levels of virus replication. As predicted, we observed that virus populations accumulated higher amounts of cbVGs in the more permissive A549 STAT1 KO cells over time; however, unexpectedly, the predominant cbVG species after passages in the two conditions were different. While A549 STAT1 KO cells accumulated relatively short cbVGs, A549 control cells mainly contained cbVGs of much longer predicted size, which have not been described previously. These long cbVGs were predominant at first in both cell lines in vitro and the predominant ones observed in samples from RSV-infected patients. Although sustained high replication levels are associated with cbVG generation and accumulation, our data show that sustained high levels of virus replication are critical for cbVG population diversification, a process that precedes the generation of shorter cbVGs that selectively accumulate over time. Taken together, we show that selection and evolution of cbVGs within a virus population are shaped by how resistant or permissive a host is to RSV.

A Virus Is a Community: Diversity within Negative-Sense RNA Virus Populations.

Negative-sense RNA virus populations are composed of diverse viral components that interact to form a community and shape the outcome of virus infections. At the genomic level, RNA virus populations consist not only of a homogeneous population of standard viral genomes but also of an extremely large number of genome variants, termed viral quasispecies, and nonstandard viral genomes, which include copy-back viral genomes, deletion viral genomes, mini viral RNAs, and hypermutated RNAs. At the particle level, RNA virus populations are composed of pleomorphic particles, particles missing or having additional genomes, and single particles or particle aggregates. As we continue discovering more about the components of negative-sense RNA virus populations and their crucial functions during virus infection, it will become more important to study RNA virus populations as a whole rather than their individual parts. In this review, we will discuss what is known about the components of negative-sense RNA virus communities, speculate how the components of the virus community interact, and summarize what vaccines and antiviral therapies are being currently developed to target or harness these components.

Harnessing defective viruses to fight infections.

One of the most pressing challenges in medicine is to develop fast and effective strategies to combat infections. Xiao et al.1 report proof-of-concept experiments for the use of defective viral genomes as broad-spectrum antivirals by harnessing their ability to stimulate the innate immune response.

Detection of respiratory syncytial virus defective genomes in nasal secretions is associated with distinct clinical outcomes.

Respiratory syncytial virus (RSV) causes respiratory illness in children, immunosuppressed individuals and the elderly. However, the viral factors influencing the clinical outcome of RSV infections remain poorly defined. Defective viral genomes (DVGs) can suppress virus replication by competing for viral proteins and by stimulating antiviral immunity. We studied the association between detection of DVGs of the copy-back type and disease severity in three RSV A-confirmed cohorts. In hospitalized children, detection of DVGs in respiratory samples at or around the time of admission associated strongly with more severe disease, higher viral load and a stronger pro-inflammatory response. Interestingly, in experimentally infected adults, the presence of DVGs in respiratory secretions differentially associated with RSV disease severity depending on when DVGs were detected. Detection of DVGs early after infection associated with low viral loads and mild disease, whereas detection of DVGs late after infection, especially if DVGs were present for prolonged periods, associated with high viral loads and severe disease. Taken together, we demonstrate that the kinetics of DVG accumulation and duration could predict clinical outcome of RSV A infection in humans, and thus could be used as a prognostic tool to identify patients at risk of worse clinical disease.

Correction: A specific sequence in the genome of respiratory syncytial virus regulates the generation of copy-back defective viral genomes.

[This corrects the article DOI: 10.1371/journal.ppat.1007707.].

A specific sequence in the genome of respiratory syncytial virus regulates the generation of copy-back defective viral genomes.

Defective viral genomes of the copy-back type (cbDVGs) are the primary initiators of the antiviral immune response during infection with respiratory syncytial virus (RSV) both in vitro and in vivo. However, the mechanism governing cbDVG generation remains unknown, thereby limiting our ability to manipulate cbDVG content in order to modulate the host response to infection. Here we report a specific genomic signal that mediates the generation of a subset of RSV cbDVG species. Using a customized bioinformatics tool, we identified regions in the RSV genome frequently used to generate cbDVGs during infection. We then created a minigenome system to validate the function of one of these sequences and to determine if specific nucleotides were essential for cbDVG generation at that position. Further, we created a recombinant virus unable to produce a subset of cbDVGs due to mutations introduced in this sequence. The identified sequence was also found as a site for cbDVG generation during natural RSV infections, and common cbDVGs originated at this sequence were found among samples from various infected patients. These data demonstrate that sequences encoded in the viral genome determine the location of cbDVG formation and, therefore, the generation of cbDVGs is not a stochastic process. These findings open the possibility of genetically manipulating cbDVG formation to modulate infection outcome.

Recent advances in vesicular stomatitis virus-based oncolytic virotherapy: a 5-year update.

Oncolytic virus (OV) therapy is an anti-cancer approach that uses viruses that preferentially infect, replicate in and kill cancer cells. Vesicular stomatitis virus (VSV, a rhabdovirus) is an OV that is currently being tested in the USA in several phase I clinical trials against different malignancies. Several factors make VSV a promising OV: lack of pre-existing human immunity against VSV, a small and easy to manipulate genome, cytoplasmic replication without risk of host cell transformation, independence of cell cycle and rapid growth to high titres in a broad range of cell lines facilitating large-scale virus production. While significant advances have been made in VSV-based OV therapy, room for improvement remains. Here we review recent studies (published in the last 5 years) that address 'old' and 'new' challenges of VSV-based OV therapy. These studies focused on improving VSV safety, oncoselectivity and oncotoxicity; breaking resistance of some cancers to VSV; preventing premature clearance of VSV; and stimulating tumour-specific immunity. Many of these approaches were based on combining VSV with other therapeutics. This review also discusses another rhabdovirus closely related to VSV, Maraba virus, which is currently being tested in Canada in phase I/II clinical trials.

Ruxolitinib and Polycation Combination Treatment Overcomes Multiple Mechanisms of Resistance of Pancreatic Cancer Cells to Oncolytic Vesicular Stomatitis Virus.

Vesicular stomatitis virus (VSV) is a promising oncolytic virus (OV). Although VSV is effective against a majority of pancreatic ductal adenocarcinoma cell (PDAC) cell lines, some PDAC cell lines are highly resistant to VSV, and the mechanisms of resistance are still unclear. JAK1/2 inhibitors (such as ruxolitinib and JAK inhibitor I) strongly stimulate VSV replication and oncolysis in all resistant cell lines but only partially improve the susceptibility of resistant PDACs to VSV. VSV tumor tropism is generally dependent on the permissiveness of malignant cells to viral replication rather than on receptor specificity, with several ubiquitously expressed cell surface molecules playing a role in VSV attachment to host cells. However, as VSV attachment to PDAC cells has never been tested before, here we examined if it was possibly inhibited in resistant PDAC cells. Our data show a dramatically weaker attachment of VSV to HPAF-II cells, the most resistant human PDAC cell line. Although sequence analysis of low-density lipoprotein (LDL) receptor (LDLR) mRNA did not reveal any amino acid substitutions in this cell line, HPAF-II cells displayed the lowest level of LDLR expression and dramatically lower LDL uptake. Treatment of cells with various statins strongly increased LDLR expression levels but did not improve VSV attachment or LDL uptake in HPAF-II cells. However, LDLR-independent attachment of VSV to HPAF-II cells was dramatically improved by treating cells with Polybrene or DEAE-dextran. Moreover, combining VSV with ruxolitinib and Polybrene or DEAE-dextran successfully broke the resistance of HPAF-II cells to VSV by simultaneously improving VSV attachment and replication.IMPORTANCE Oncolytic virus (OV) therapy is an anticancer approach that uses viruses that selectively infect and kill cancer cells. This study focuses on oncolytic vesicular stomatitis virus (VSV) against pancreatic ductal adenocarcinoma (PDAC) cells. Although VSV is effective against most PDAC cells, some are highly resistant to VSV, and the mechanisms are still unclear. Here we examined if VSV attachment to cells was inhibited in resistant PDAC cells. Our data show very inefficient attachment of VSV to the most resistant human PDAC cell line, HPAF-II. However, VSV attachment to HPAF-II cells was dramatically improved by treating cells with polycations. Moreover, combining VSV with polycations and ruxolitinib (which inhibits antiviral signaling) successfully broke the resistance of HPAF-II cells to VSV by simultaneously improving VSV attachment and replication. We envision that this novel triple-combination approach could be used in the future to treat PDAC tumors that are highly resistant to OV therapy.

Novel biomarkers of resistance of pancreatic cancer cells to oncolytic vesicular stomatitis virus.

Vesicular stomatitis virus (VSV) based recombinant viruses (such as VSV-ΔM51) are effective oncolytic viruses (OVs) against a majority of pancreatic ductal adenocarcinoma (PDAC) cell lines. However, some PDAC cell lines are highly resistant to VSV-ΔM51. We recently showed that treatment of VSV-resistant PDAC cells with ruxolitinib (JAK1/2 inhibitor) or TPCA-1 (IKK-ß inhibitor) breaks their resistance to VSV-ΔM51. Here we compared the global effect of ruxolitinib or TPCA-1 treatment on cellular gene expression in PDAC cell lines highly resistant to VSV-ΔM51. Our study identified a distinct subset of 22 interferon-stimulated genes (ISGs) downregulated by both ruxolitinib and TPCA-1. Further RNA and protein analyses demonstrated that 4 of these genes (MX1, EPSTI1, XAF1, and GBP1) are constitutively co-expressed in VSV-resistant, but not in VSV-permissive PDACs, thus serving as potential biomarkers to predict OV therapy success. Moreover, shRNA-mediated knockdown of one of such ISG, MX1, showed a positive effect on VSV-ΔM51 replication in resistant PDAC cells, suggesting that at least some of the identified ISGs contribute to resistance of PDACs to VSV-ΔM51. As certain oncogene and tumor suppressor gene variants are often associated with increased tropism of OVs to cancer cells, we also analyzed genomic DNA in a set of PDAC cell lines for frequently occurring cancer associated mutations. While no clear correlation was found between such mutations and resistance of PDACs to VSV-ΔM51, the analysis generated valuable genotypic data for future studies.

Breaking resistance of pancreatic cancer cells to an attenuated vesicular stomatitis virus through a novel activity of IKK inhibitor TPCA-1.

Vesicular stomatitis virus (VSV) is an effective oncolytic virus against most human pancreatic ductal adenocarcinoma (PDAC) cell lines. However, some PDAC cell lines are highly resistant to oncolytic VSV-ΔM51 infection. To better understand the mechanism of resistance, we tested a panel of 16 small molecule inhibitors of different cellular signaling pathways, and identified TPCA-1 (IKK-ß inhibitor) and ruxolitinib (JAK1/2 inhibitor), as strong enhancers of VSV-ΔM51 replication and virus-mediated oncolysis in all VSV-resistant PDAC cell lines. Both TPCA-1 and ruxolitinib similarly inhibited STAT1 and STAT2 phosphorylation and decreased expression of antiviral genes MxA and OAS. Moreover, an in situ kinase assay provided biochemical evidence that TPCA-1 directly inhibits JAK1 kinase activity. Together, our data demonstrate that TPCA-1 is a unique dual inhibitor of IKK-ß and JAK1 kinase, and provide a new evidence that upregulated type I interferon signaling plays a major role in resistance of pancreatic cancer cells to oncolytic viruses.

Induction of apoptosis in pancreatic cancer cells by vesicular stomatitis virus.

Effective oncolytic virus (OV) therapy is dependent on the ability of replication-competent viruses to kill infected cancer cells. We previously showed that human pancreatic ductal adenocarcinoma (PDAC) cell lines are highly heterogeneous in their permissiveness to vesicular stomatitis virus (VSV), in part due to differences in type I interferon (IFN) signaling. Here, using 10 human PDAC cell lines and three different VSV recombinants (expressing ΔM51 or wild type matrix protein), we examined cellular and viral factors affecting VSV-mediated apoptosis activation in PDACs. In most cell lines, VSVs activated both extrinsic and intrinsic apoptosis pathways, and VSV-ΔM51 primarily activated the type II extrinsic pathway. In cells with defective IFN signaling, all VSV recombinants induced robust apoptosis, whereas VSV-ΔM51 was a more effective apoptosis activator in PDACs with virus-inducible IFN signaling. Three cell lines constitutively expressing high levels of IFN-stimulated genes (ISGs) were resistant to apoptosis under most experimental conditions, even when VSV replication levels were dramatically increased by Jak inhibitor I treatment. Two of these cell lines also poorly activated apoptosis when treated with Fas activating antibody, suggesting a general defect in apoptosis.