Oncolytic Viruses and the Immune System: The Dynamic Duo.

Oncolytic viruses (OVs) constitute a new and promising immunotherapeutic approach toward cancer treatment. This therapy takes advantage of the natural propensity of most tumor cells to be infected by specific OVs. Besides the direct killing potential (oncolysis), what makes OV administration attractive for the present cancer immunotherapeutic scenario is the capacity to induce two new overlapping, but distinct, immunities: anti-tumoral and anti-viral. OV infection and oncolysis naturally elicit both innate and adaptive immune responses (required for long-term anti-tumoral immunity); at the same time, the viral infection prompts an anti-viral response. In this review, we discuss the dynamic interaction between OVs and the triggered responses of the immune system. The anti-OV immunological events that lead to viral clearance and the strategies to deal with such potential loss of the therapeutic virus are discussed. Additionally, we review the immune stimulatory actions induced by OVs through different inherent strategies, such as modulation of the tumor microenvironment, the role of immunogenic cell death, and the consequences of genetically modifying OVs by arming them with therapeutic transgenes. An understanding of the balance between the OV-induced anti-tumoral versus anti-viral immunities will provide insight when choosing the appropriate virotherapy for any specific cancer.

Not so unique to Primates: The independent adaptive evolution of TRIM5 in Lagomorpha lineage.

The plethora of restriction factors with the ability to inhibit the replication of retroviruses have been widely studied and genetic hallmarks of evolutionary selective pressures in Primates have been well documented. One example is the tripartite motif-containing protein 5 alpha (TRIM5α), a cytoplasmic factor that restricts retroviral infection in a species-specific fashion. In Lagomorphs, similarly to what has been observed in Primates, the specificity of TRIM5 restriction has been assigned to the PRYSPRY domain. In this study, we present the first insight of an intra-genus variability within the Lagomorpha TRIM5 PRYSPRY domain. Remarkably, and considering just the 32 residue-long v1 region of this domain, the deduced amino acid sequences of Daurian pika (Ochotona dauurica) and steppe pika (O. pusilla) evidenced a high divergence when compared to the remaining Ochotona species, presenting values of 44% and 66% of amino acid differences, respectively. The same evolutionary pattern was also observed when comparing the v1 region of two Sylvilagus species members (47% divergence). However, and unexpectedly, the PRYSPRY domain of Lepus species exhibited a great conservation. Our results show a high level of variation in the PRYSPRY domain of Lagomorpha species that belong to the same genus. This suggests that, throughout evolution, the Lagomorpha TRIM5 should have been influenced by constant selective pressures, likely as a result of multiple different retroviral infections.

Genetic Characterization of a Recombinant Myxoma Virus in the Iberian Hare (Lepus granatensis).

Myxomatosis is a lethal disease in wild European and domestic rabbits (Oryctolagus cuniculus), which is caused by a Myxoma virus (MYXV) infection-a leporipoxvirus that is found naturally in some Sylvilagus rabbit species in South America and California. The introduction of MYXV into feral European rabbit populations of Australia and Europe, in the early 1950s, demonstrated the best-documented field example of host-virus coevolution, following a cross-species transmission. Recently, a new cross-species jump of MYXV has been suggested in both Great Britain and Spain, where European brown hares (Lepus europaeus) and Iberian hares (Lepus granatensis) were found dead with lesions consistent with those observed in myxomatosis. To investigate the possibility of a new cross-species transmission event by MYXV, tissue samples collected from a wild Iberian hare found dead in Spain (Toledo region) were analyzed and deep sequenced. Our results reported a new MYXV isolate (MYXV Toledo) in the tissues of this species. The genome of this new virus was found to encode three disruptive genes (M009L, M036L, and M152R) and a novel ~2.8 kb recombinant region, which resulted from an insertion of four novel poxviral genes towards the 3' end of the negative strand of its genome. From the open reading frames inserted into the MYXV Toledo virus, a new orthologue of a poxvirus host range gene family member was identified, which was related to the MYXV gene M064R. Overall, we confirmed the identity of a new MYXV isolate in Iberian hares, which, we hypothesized, was able to more effectively counteract the host defenses in hares and start an infectious process in this new host.

Parallel adaptation of rabbit populations to myxoma virus.

In the 1950s the myxoma virus was released into European rabbit populations in Australia and Europe, decimating populations and resulting in the rapid evolution of resistance. We investigated the genetic basis of resistance by comparing the exomes of rabbits collected before and after the pandemic. We found a strong pattern of parallel evolution, with selection on standing genetic variation favoring the same alleles in Australia, France, and the United Kingdom. Many of these changes occurred in immunity-related genes, supporting a polygenic basis of resistance. We experimentally validated the role of several genes in viral replication and showed that selection acting on an interferon protein has increased the protein's antiviral effect.

Ex Vivo Oncolytic Virotherapy with Myxoma Virus Arms Multiple Allogeneic Bone Marrow Transplant Leukocytes to Enhance Graft versus Tumor.

Allogeneic stem cell transplant-derived T cells have the potential to seek and eliminate sites of residual cancer that escaped primary therapy. Oncolytic myxoma virus (MYXV) exhibits potent anti-cancer efficacy against human cancers like multiple myeloma (MM) and can arm transplant-derived T cells to become more effective cancer killers in vitro and in an immunodeficient xenotransplant murine model. Here, we tested ex vivo MYXV virotherapy against residual murine MM in immunocompetent mice using an allogeneic mouse-mouse model. In contrast to all human MM cell lines previously tested, the murine MM cell line tested here was highly resistant to direct MYXV infection and oncolysis in vitro. Despite this in vitro resistance, we found that ex vivo MYXV-armed allogeneic bone marrow (BM) transplantation dramatically ablated pre-seeded residual MM in vivo. Unexpectedly, we show that both neutrophils and activated T cells from the donor function as virus-armed carrier cells, and MYXV-preloaded cells enhanced MM killing. Our results demonstrate a novel therapeutic paradigm for residual cancer, in which multiple classes of allotransplant leukocytes can be armed by MYXV ex vivo to enhance the graft-versus-tumor effects.

An overview of the lagomorph immune system and its genetic diversity.

Our knowledge of the lagomorph immune system remains largely based upon studies of the European rabbit (Oryctolagus cuniculus), a major model for studies of immunology. Two important and devastating viral diseases, rabbit hemorrhagic disease and myxomatosis, are affecting European rabbit populations. In this context, we discuss the genetic diversity of the European rabbit immune system and extend to available information about other lagomorphs. Regarding innate immunity, we review the most recent advances in identifying interleukins, chemokines and chemokine receptors, Toll-like receptors, antiviral proteins (RIG-I and Trim5), and the genes encoding fucosyltransferases that are utilized by rabbit hemorrhagic disease virus as a portal for invading host respiratory and gut epithelial cells. Evolutionary studies showed that several genes of innate immunity are evolving by strong natural selection. Studies of the leporid CCR5 gene revealed a very dramatic change unique in mammals at the second extracellular loop of CCR5 resulting from a gene conversion event with the paralogous CCR2. For the adaptive immune system, we review genetic diversity at the loci encoding antibody variable and constant regions, the major histocompatibility complex (RLA) and T cells. Studies of IGHV and IGKC genes expressed in leporids are two of the few examples of trans-species polymorphism observed outside of the major histocompatibility complex. In addition, we review some endogenous viruses of lagomorph genomes, the importance of the European rabbit as a model for human disease studies, and the anticipated role of next-generation sequencing in extending knowledge of lagomorph immune systems and their evolution.

Endogenization of mouse mammary tumor virus (MMTV)-like elements in genomes of pikas (Ochotona sp.).

Despite the finding in European rabbit and other leporid genomes of the first ever described endogenous lentivirus and of a European rabbit exclusive endogenous gammaretrovirus, until now no exogenous retroviruses have been isolated in Lagomorpha species. Nevertheless, looking for the presence of endogenous retroviruses (ERVs) in the species genomes could lead to the discovery of retroviral lineages yet to be found in Lagomorpha. Different mammalian genomes harbor endogenous viral sequences phylogenetically close to the betaretrovirus mouse mammary tumor virus (MMTV), propelling us to look for such retroviral "fossil" in American pika (Ochotona princeps) and European rabbit (Oryctolagus cuniculus) genomes. By performing genomic mining using MMTV gag and LTR as query sequences, we found that such viral elements were absent from the European rabbit genome. Oppositely, significant matches were found in American pika, and more importantly, a nearly complete MMTV-like virus (Pika-BERV) was identified. Using Pika-BERV gag and LTR as templates, we found similar sequences endogenized in different pika (Ochotona sp.) species. The orthology of the LTR flanking region between some pika species supported shared ancestry of specific endogenous betaretroviruses, while in other pika species similar sequences, but not orthologous, should have resulted from independent insertions. Our study supports the possible existence of infecting exogenous betaretroviruses for a long term, after the divergence of Ochotonidae from Leporidae, but yet to be identified.

The phylogeny of pikas (Ochotona) inferred from a multilocus coalescent approach.

The clarification of the systematics of pikas (genus Ochotona) has been hindered by largely overlapping morphological characters among species and the lack of a comprehensive molecular phylogeny. Here we estimate the first multilocus phylogeny of the genus to date, by analysing 12 nuclear DNA markers (total of 7.5Kb) in 11 species of pikas from the four classified subgenera (Pika, Ochotona, Lagotona and Conothoa) using a multispecies coalescent-based framework. The species-tree confirmed the subgeneric classification by retrieving as monophyletic the subgenera represented here by more than one species. Contrary to previous phylogenies based on mtDNA alone, Lagotona was found to be sister to Pika. Also, support for the monophyly of the alpina group was not strong, thus caution should be used in future analyses of this group. A relaxed molecular clock calibrated using the Ochotonidae-Leporidae divergence resulted in more recent estimates of divergence times relative to previous studies. Strong concordance with inferences based on fossil records was found, suggesting that the initial diversification of the genus took place by the end of late Miocene. Finally, this work sets up methodologies and gathers molecular markers that can be used to extend the understanding of the evolutionary history of the genus.

Evolution of viral sensing RIG-I-like receptor genes in Leporidae genera Oryctolagus, Sylvilagus, and Lepus.

One of the most severe European rabbit (Oryctolagus cuniculus) pathogens is myxoma virus (MYXV), a rabbit-specific leporipoxvirus that causes the highly lethal disease myxomatosis. Other leporid genera, Sylvilagus and Lepus, encompass species with variable susceptibilities to MYXV, but these do not develop the lethal form of the disease. The protective role of the retinoic acid-inducible gene-I (RIG-I/DDX58) in sensing MYXV in nonpermissive human myeloid cells prompted the study of the RIG-I-like receptor (RLR) family evolution in the three leporid genera. This viral-sensor family also includes the melanoma differentiation-associated factor 5 (MDA5/IFIH1), and the laboratory of genetics and physiology 2 (LGP2/DHX58). Considering specifically the MYXV susceptible host (European rabbit) and one of the virus natural long-term hosts (Sylvilagus bachmani, brush rabbit), the amino acid differences of positively selected sites in RIG-I between the two species were located in the protein region responsible for viral RNA recognition and binding, the repressor domain. Such differences might play a determinant role in how MYXV is sensed. When looking for episodic selection on MDA5 and LGP2 of the eastern cottontail (Sylvilagus floridanus), we also uncovered evidence of selective pressures that might be exerted by a species-specific leporipoxvirus, the Shope fibroma virus. Finally, a putative alternative splicing case was identified in Oryctolagus and Lepus MDA5 isoforms, corresponding to the deletion of one specific exon. This study provided the first insights into the evolution of the leporid RLR gene family that helps illuminate the origins of the species-specific innate responses to pathogens and more specifically to MYXV.

Positive evolutionary selection on the RIG-I-like receptor genes in mammals.

The mammalian RIG-I-like receptors, RIG-I, MDA5 and LGP2, are a family of DExD/H box RNA helicases responsible for the cytoplasmic detection of viral RNA. These receptors detect a variety of RNA viruses, or DNA viruses that express unusual RNA species, many of which are responsible for a great number of severe and lethal diseases. Host innate sentinel proteins involved in pathogen recognition must rapidly evolve in a dynamic arms race with pathogens, and thus are subjected to long-term positive selection pressures to avoid potential infections. Using six codon-based Maximum Likelihood methods, we were able to identify specific codons under positive selection in each of these three genes. The highest number of positively selected codons was detected in MDA5, but a great percentage of these codons were located outside of the currently defined protein domains for MDA5, which likely reflects the imposition of both functional and structural constraints. Additionally, our results support LGP2 as being the least prone to evolutionary change, since the lowest number of codons under selection was observed for this gene. On the other hand, the preponderance of positively selected codons for RIG-I were detected in known protein functional domains, suggesting that pressure has been imposed by the vast number of viruses that are recognized by this RNA helicase. Furthermore, the RIG-I repressor domain, the region responsible for recognizing and binding to its RNA substrates, exhibited the strongest evidence of selective pressures. Branch-site analyses were performed and several species branches on the three receptor gene trees showed evidence of episodic positive selection. In conclusion, by looking for evidence of positive evolutionary selection on mammalian RIG-I-like receptor genes, we propose that a multitude of viruses have crafted the receptors biological function in host defense, specifically for the RIG-I gene, contributing to the innate species-specific resistance/susceptibility to diverse viral pathogens.

An intriguing shift occurs in the novel protein phosphatase 1 binding partner, TCTEX1D4: evidence of positive selection in a pika model.

T-complex testis expressed protein 1 domain containing 4 (TCTEX1D4) contains the canonical phosphoprotein phosphatase 1 (PPP1) binding motif, composed by the amino acid sequence RVSF. We identified and validated the binding of TCTEX1D4 to PPP1 and demonstrated that indeed this protein is a novel PPP1 interacting protein. Analyses of twenty-one mammalian species available in public databases and seven Lagomorpha sequences obtained in this work showed that the PPP1 binding motif 90RVSF93 is present in all of them and is flanked by a palindromic sequence, PLGS, except in three species of pikas (Ochotona princeps, O. dauurica and O. pusilla). Furthermore, for the Ochotona species an extra glycosylation site, motif 96NLS98, and the loss of the palindromic sequence were observed. Comparison with other lagomorphs suggests that this event happened before the Ochotona radiation. The dN/dS for the sequence region comprising the PPP1 binding motif and the flanking palindrome highly supports the hypothesis that for Ochotona species this region has been evolving under positive selection. In addition, mutational screening shows that the ability of pikas TCTEX1D4 to bind to PPP1 is maintained, although the PPP1 binding motif is disrupted, and the N- and C-terminal surrounding residues are also abrogated. These observations suggest pika as an ideal model to study novel PPP1 complexes regulatory mechanisms.

Evolution and divergence of the mammalian SAMD9/SAMD9L gene family.

BACKGROUND: The physiological functions of the human Sterile Alpha Motif Domain-containing 9 (SAMD9) gene and its chromosomally adjacent paralogue, SAMD9-like (SAMD9L), currently remain unknown. However, the direct links between the deleterious mutations or deletions in these two genes and several human disorders, such as inherited inflammatory calcified tumors and acute myeloid leukemia, suggest their biological importance. SAMD9 and SAMD9L have also recently been shown to play key roles in the innate immune responses to stimuli such as viral infection. We were particularly interested in understanding the mammalian evolutionary history of these two genes. The phylogeny of SAMD9 and SAMD9L genes was reconstructed using the Maximum Likelihood method. Furthermore, six different methods were applied to detect SAMD9 and SAMD9L codons under selective pressure: the site-specific model M8 implemented in the codeml program in PAML software and five methods available on the Datamonkey web server, including the Single Likelihood Ancestor Counting method, the Fixed Effect Likelihood method, the Random Effect Likelihood method, the Mixed Effects Model of Evolution method and the Fast Unbiased Bayesian AppRoximation method. Additionally, the house mouse (Mus musculus) genome has lost the SAMD9 gene, while keeping SAMD9L intact, prompting us to investigate whether this loss is a unique event during evolution. RESULTS: Our evolutionary analyses suggest that SAMD9 and SAMD9L arose through an ancestral gene duplication event after the divergence of Marsupialia from Placentalia. Additionally, selection analyses demonstrated that both genes have been subjected to positive evolutionary selection. The absence of either SAMD9 or SAMD9L genes from some mammalian species supports a partial functional redundancy between the two genes. CONCLUSIONS: To the best of our knowledge, this work is the first study on the evolutionary history of mammalian SAMD9 and SAMD9L genes. We conclude that evolutionary selective pressure has acted on both of these two genes since their divergence, suggesting their importance in multiple cellular processes, such as the immune responses to viral pathogens.