Cory Batenchuk

Model-based rational design of an oncolytic virus with improved therapeutic potential

Oncolytic viruses are complex biological agents that interact at multiple levels with both tumour and normal tissues. Antiviral pathways induced by interferon are known to have a critical role in determining tumour cell sensitivity and normal cell resistance to infection with oncolytic viruses. Here we pursue a synthetic biology approach to identify methods that enhance antitumour activity of oncolytic viruses through suppression of interferon signalling. On the basis of the mathematical analysis of multiple strategies, we hypothesize that a positive feedback loop, established by virus-mediated expression of a soluble interferon-binding decoy receptor, increases tumour cytotoxicity without compromising normal cells. Oncolytic rhabdoviruses engineered to express a secreted interferon antagonist have improved oncolytic potential in cellular cancer models, and display improved therapeutic potential in tumour-bearing mice. Our results demonstrate the potential of this methodology in evaluating potential caveats of viral immune-evasion strategies and improving the design oncolytic viruses.

VEGF-Mediated Induction of PRD1-BF1/Blimp1 Expression Sensitizes Tumor Vasculature to Oncolytic Virus Infection.

Oncolytic viruses designed to attack malignant cells can in addition infect and destroy tumor vascular endothelial cells. We show here that this expanded tropism of oncolytic vaccinia virus to the endothelial compartment is a consequence of VEGF-mediated suppression of the intrinsic antiviral response. VEGF/VEGFR2 signaling through Erk1/2 and Stat3 leads to upregulation, nuclear localization, and activation of the transcription repressor PRD1-BF1/Blimp1. PRD1-BF1 does not contribute to the mitogenic effects of VEGF, but directly represses genes involved in type I interferon (IFN)-mediated antiviral signaling. In vivo suppression of VEGF signaling diminishes PRD1-BF1/Blimp1 expression in tumor vasculature and inhibits intravenously administered oncolytic vaccinia delivery to and consequent spread within the tumor.

Chromosomal Position Effects Are Linked to Sir2-Mediated Variation in Transcriptional Burst Size

Gene expression noise varies with genomic position and is a driving force in the evolution of chromosome organization. Nevertheless, position effects remain poorly characterized. Here, we present a systematic analysis of chromosomal position effects by characterizing single-cell gene expression from euchromatic positions spanning the length of a eukaryotic chromosome. We demonstrate that position affects gene expression by modulating the size of transcriptional bursts, rather than their frequency, and that the histone deacetylase Sir2 plays a role in this process across the chromosome.

Microtubule disruption synergizes with oncolytic virotherapy by inhibiting interferon translation and potentiating bystander killing

In this study, we show that several microtubule-destabilizing agents used for decades for treatment of cancer and other diseases also sensitize cancer cells to oncolytic rhabdoviruses and improve therapeutic outcomes in resistant murine cancer models. Drug-induced microtubule destabilization leads to superior viral spread in cancer cells by disrupting type I IFN mRNA translation, leading to decreased IFN protein expression and secretion. Furthermore, microtubule-destabilizing agents specifically promote cancer cell death following stimulation by a subset of infection-induced cytokines, thereby increasing viral bystander effects. This study reveals a previously unappreciated role for microtubule structures in the regulation of the innate cellular antiviral response and demonstrates that unexpected combinations of approved chemotherapeutics and biological agents can lead to improved therapeutic outcomes.

Non-replicating rhabdovirus-derived particles (NRRPs) eradicate acute leukemia by direct cytolysis and induction of antitumor immunity.

Rhabdoviruses (RVs) are currently being pursued as anticancer therapeutics for various tumor types, notably leukemia. However, modest virion production and limited spread between noncontiguous circulating leukemic cells requires high-dose administration of RVs, which exceeds the maximum tolerable dose of the live virus. Furthermore, in severely immunosuppressed leukemic patients, the potential for uncontrolled live virus spread may compromise the safety of a live virus approach. We hypothesized that the barriers to oncolytic virotherapy in liquid tumors may be overcome by administration of high-dose non-replicating RVs. We have developed a method to produce unique high-titer bioactive yet non-replicating rhabdovirus-derived particles (NRRPs). This novel biopharmaceutical is multimodal possessing direct cytolytic and immunomodulatory activity against acute leukemia. We demonstrate that NRRP resistance in normal cells is mediated by intact antiviral defences including interferon (IFN). This data was substantiated using murine models of blast crisis. The translational promise of NRRPs was demonstrated in clinical samples obtained from patients with high-burden multidrug-resistant acute myeloid leukemia. This is the first successful attempt to eradicate disseminated cancer using a non-replicating virus-derived agent, representing a paradigm shift in our understanding of oncolytic virus-based therapies and their application toward the treatment of acute leukemia.

Quantitative epistasis analysis and pathway inference from genetic interaction data.

Inferring regulatory and metabolic network models from quantitative genetic interaction data remains a major challenge in systems biology. Here, we present a novel quantitative model for interpreting epistasis within pathways responding to an external signal. The model provides the basis of an experimental method to determine the architecture of such pathways, and establishes a new set of rules to infer the order of genes within them. The method also allows the extraction of quantitative parameters enabling a new level of information to be added to genetic network models. It is applicable to any system where the impact of combinatorial loss-of-function mutations can be quantified with sufficient accuracy. We test the method by conducting a systematic analysis of a thoroughly characterized eukaryotic gene network, the galactose utilization pathway in Saccharomyces cerevisiae. For this purpose, we quantify the effects of single and double gene deletions on two phenotypic traits, fitness and reporter gene expression. We show that applying our method to fitness traits reveals the order of metabolic enzymes and the effects of accumulating metabolic intermediates. Conversely, the analysis of expression traits reveals the order of transcriptional regulatory genes, secondary regulatory signals and their relative strength. Strikingly, when the analyses of the two traits are combined, the method correctly infers ∼80% of the known relationships without any false positives.

Identification of response-modulated genetic interactions by sensitivity-based epistatic analysis

BackgroundHigh-throughput genomics has enabled the global mapping of genetic interactions based on the phenotypic impact of combinatorial genetic perturbations. An important next step is to understand how these networks are dynamically remodelled in response to environmental stimuli. Here, we report on the development and testing of a method to identify such interactions. The method was developed from first principles by treating the impact on cellular growth of environmental perturbations equivalently to that of gene deletions. This allowed us to establish a novel neutrality function marking the absence of epistasis in terms of sensitivity phenotypes rather than fitness. We tested the method by identifying fitness- and sensitivity-based interactions involved in the response to drug-induced DNA-damage of budding yeast Saccharomyces cerevisiae using two mutant libraries - one containing transcription factor deletions, and the other containing deletions of DNA repair genes.ResultsWithin the library of transcription factor deletion mutants, we observe significant differences in the sets of genetic interactions identified by the fitness- and sensitivity-based approaches. Notably, among the most likely interactions, only ~50% were identified by both methods. While interactions identified solely by the sensitivity-based approach are modulated in response to drug-induced DNA damage, those identified solely by the fitness-based method remained invariant to the treatment. Comparison of the identified interactions to transcriptional profiles and protein-DNA interaction data indicate that the sensitivity-based method improves the identification of interactions involved in the DNA damage response. Additionally, for the library containing DNA repair mutants, we observe that the sensitivity-based method improves the grouping of functionally related genes, as well as the identification of protein complexes, involved in DNA repair.ConclusionOur results show that the identification of response-modulated genetic interactions can be improved by incorporating the effect of a changing environment directly into the neutrality function marking the absence of epistasis. We expect that this extension of conventional epistatic analysis will facilitate the development of dynamic models of gene networks from quantitative measurements of genetic interactions. While the method was developed for growth phenotype, it should apply equally well for other phenotypes, including the expression of fluorescent reporters.

High-throughput Titration of Luciferase-expressing Recombinant Viruses

Standard plaque assays to determine infectious viral titers can be time consuming, are not amenable to a high volume of samples, and cannot be done with viruses that do not form plaques. As an alternative to plaque assays, we have developed a high-throughput titration method that allows for the simultaneous titration of a high volume of samples in a single day. This approach involves infection of the samples with a Firefly luciferase tagged virus, transfer of the infected samples onto an appropriate permissive cell line, subsequent addition of luciferin, reading of plates in order to obtain luminescence readings, and finally the conversion from luminescence to viral titers. The assessment of cytotoxicity using a metabolic viability dye can be easily incorporated in the workflow in parallel and provide valuable information in the context of a drug screen. This technique provides a reliable, high-throughput method to determine viral titers as an alternative to a standard plaque assay.