Maladaptation after a virus host switch leads to increased activation of the pro-inflammatory NF-κB pathway.

Myxoma virus (MYXV) causes localized cutaneous fibromas in its natural hosts, tapeti and brush rabbits; however, in the European rabbit, MYXV causes the lethal disease myxomatosis. Currently, the molecular mechanisms underlying this increased virulence after cross-species transmission are poorly understood. In this study, we investigated the interaction between MYXV M156 and the host protein kinase R (PKR) to determine their crosstalk with the proinflammatory nuclear factor kappa B (NF-κB) pathway. Our results demonstrated that MYXV M156 inhibits brush rabbit PKR (bPKR) more strongly than European rabbit PKR (ePKR). This moderate ePKR inhibition could be improved by hyperactive M156 mutants. We hypothesized that the moderate inhibition of ePKR by M156 might incompletely suppress the signal transduction pathways modulated by PKR, such as the NF-κB pathway. Therefore, we analyzed NF-κB pathway activation with a luciferase-based promoter assay. The moderate inhibition of ePKR resulted in significantly higher NF-κB­dependent reporter activity than complete inhibition of bPKR. We also found a stronger induction of the NF-κB target genes TNFα and IL-6 in ePKR-expressing cells than in bPKR-expressing cells in response to M156 in both transfection and infections assays. Furthermore, a hyperactive M156 mutant did not cause ePKR-dependent NF-κB activation. These observations indicate that M156 is maladapted for ePKR inhibition, only incompletely blocking translation in these hosts, resulting in preferential depletion of short­half-life proteins, such as the NF-κB inhibitor IκBα. We speculate that this functional activation of NF-κB induced by the intermediate inhibition of ePKR by M156 may contribute to the increased virulence of MYXV in European rabbits.

Effect of Myxoma Virus Species Jump on Iberian Hare Populations.

The myxoma virus species jump from European rabbits (Oryctolagus cuniculus) to Iberian hares (Lepus granatensis) has raised concerns. We assess the decline suffered by Iberian hare populations on the Iberian Peninsula and discuss the association between the effect of myxomatosis and the average abundance index, which we estimated by using hunting bags.

High Levels of Extracellular Potassium Can Delay Myxoma Virus Replication by Preventing Release of Virions from the Endosomes.

Potassium (K+) is one of the most abundant cations in the human body. Under normal conditions, the vast majority of K+ is found within cells, and the extracellular [K+] is tightly regulated to within 3.0 to 5.0 mM. However, it has recently been shown that high levels of localized necrosis can increase the extracellular concentration of K+ to above 50 mM. This raises the possibility that elevated extracellular K+ might influence a variety of biological processes that occur within regions of necrotic tissue. For example, K+ has been shown to play a central role in the replication cycles of numerous viral families, and in cases of lytic infection, localized regions containing large numbers of necrotic cells can be formed. Here, we show that the replication of the model poxvirus myxoma virus (MYXV) is delayed by elevated levels of extracellular K+. These increased K+ concentrations alter the cellular endocytic pathway, leading to increased phagocytosis but a loss of endosomal/lysosomal segregation. This slows the release of myxoma virus particles from the endosomes, resulting in delays in genome synthesis and infectious particle formation as well as reduced viral spread. Additionally, mathematical modeling predicts that the extracellular K+ concentrations required to impact myxoma virus replication can be reached in viral lesions under a variety of conditions. Taken together, these data suggest that the extracellular [K+] plays a role in determining the outcomes of myxoma infection and that this effect could be physiologically relevant during pathogenic infection. IMPORTANCE Intracellular K+ homeostasis has been shown to play a major role in the replication of numerous viral families. However, the potential impact of altered extracellular K+ concentrations is less well understood. Our work demonstrates that increased concentrations of extracellular K+ can delay the replication cycle of the model poxvirus MYXV by inhibiting virion release from the endosomes. Additionally, mathematical modeling predicts that the levels of extracellular K+ required to impact MYXV replication can likely be reached during pathogenic infection. These results suggest that localized viral infection can alter K+ homeostasis and that these alterations might directly affect viral pathogenesis.

Comparative analysis of two novel complete genomes of myxoma virus vaccine strains.

Myxoma virus (MYXV) is a double-stranded DNA-containing virus of the family Poxviridae, genus Leporipoxvirus. MYXV is an important model virus for evolutionary and immunological research and a promising oncolytic. In this study, we sequenced and analyzed two complete genomes of MYXV virus vaccine strains B-82 and Rabbivac-B, which are widely used for vaccine production in Russia. Here, we first show that MYXV vaccine strains B-82 and Rabbivac-B share a common origin with the American recombinant MYXV MAV vaccine strain. In addition, our data suggest that the MYXV B-82 and Rabbivac-B strains contain a number of genes at the 5' and 3' ends that are identical to the virulent MYXV Lausanne strain. Several unique genetic signatures were identified in the M013L, M017L, M023, and M121R genes, helping to achieve high genetic resolution between vaccine strains. Overall, these findings highlight the evolutionary flexibility of certain genes in the MYXV genome and provide insights into the molecular epidemiology of the virus and subsequent vaccine development.

Divergent Evolutionary Pathways of Myxoma Virus in Australia: Virulence Phenotypes in Susceptible and Partially Resistant Rabbits Indicate Possible Selection for Transmissibility.

To characterize the ongoing evolution of myxoma virus in Australian rabbits, we used experimental infections of laboratory rabbits to determine the virulence and disease phenotypes of recent virus isolates. The viruses, collected between 2012 and 2015, fell into three lineages, one of which, lineage c, experienced a punctuated increase in evolutionary rate. All viruses were capable of causing acute death with aspects of neutropenic septicemia, characterized by minimal signs of myxomatosis, the occurrence of pulmonary edema and bacteria invasions throughout internal organs, but with no inflammatory response. For the viruses of highest virulence all rabbits usually died at this point. In more attenuated viruses, some rabbits died acutely, while others developed an amyxomatous phenotype. Rabbits that survived for longer periods developed greatly swollen cutaneous tissues with very high virus titers. This was particularly true of lineage c viruses. Unexpectedly, we identified a line of laboratory rabbits with some innate resistance to myxomatosis and used these in direct comparisons with the fully susceptible rabbit line. Importantly, the same disease phenotype occurred in both susceptible and resistant rabbits, although virulence was shifted toward more attenuated grades in resistant animals. We propose that selection against inflammation at cutaneous sites prolongs virus replication and enhances transmission, leading to the amyxomatous phenotype. In some virus backgrounds this creates an immunosuppressive state that predisposes to high virulence and acute death. The alterations in disease pathogenesis, particularly the overwhelming bacterial invasions that characterize the modern viruses, suggest that their virulence grades are not directly comparable with earlier studies. IMPORTANCE The evolution of the myxoma virus (MYXV) following its release as a biological control for European rabbits in Australia is the textbook example of the coevolution of virus virulence and host resistance. However, most of our knowledge of MYXV evolution only covers the first few decades of its spread in Australia and often with little direct connection between how changes in virus phenotype relate to those in the underlying virus genotype. By conducting detailed experimental infections of recent isolates of MYXV in different lines of laboratory rabbits, we examined the ongoing evolution of MYXV disease phenotypes. Our results reveal a wide range of phenotypes, including an amyxomatous type, as well as the impact of invasive bacteria, that in part depended on the level of rabbit host resistance. These results provide a unique insight into the complex virus and host factors that combine to shape disease phenotype and viral evolution.

Quantifying resistance to myxomatosis in wild rabbits produces novel evolutionary insights.

Wild rabbits in Australia developed genetic resistance to the myxoma virus, which was introduced as a biological control agent. However, little is known about the rate at which this evolutionary change occurred. We collated data from challenge trials that estimated rabbit resistance to myxomatosis in Australia and expressed resistance on a continuous scale, enabling trends in its development to be assessed over 45 years up to 1995. Resistance initially increased rapidly, followed by a plateau lasting ten years, before a second rapid increase occurred associated with the introduction of European rabbit fleas as myxoma virus vectors. By contrast, in the United Kingdom, where rabbit flea vectors were already present when the myxoma virus initially spread, resistance developed more slowly. No estimates of rabbit resistance to myxomatosis have been made for almost 30 years, despite other highly lethal rabbit pathogens becoming established worldwide. Continued testing of wild-caught rabbits in Australia to determine current levels of resistance to myxomatosis is recommended to assess its current effectiveness for managing pest rabbits. Given the economic and environmental significance of invasive rabbits, it would be remiss to manage such biological resources and ecosystem services poorly.

RNA Helicase A/DHX9 Forms Unique Cytoplasmic Antiviral Granules That Restrict Oncolytic Myxoma Virus Replication in Human Cancer Cells.

RNA helicase A/DHX9 is required for diverse RNA-related essential cellular functions and antiviral responses and is hijacked by RNA viruses to support their replication. Here, we show that during the late replication stage in human cancer cells of myxoma virus (MYXV), a member of the double-stranded DNA (dsDNA) poxvirus family that is being developed as an oncolytic virus, DHX9, forms unique granular cytoplasmic structures, which we named "DHX9 antiviral granules." These DHX9 antiviral granules are not formed if MYXV DNA replication and/or late protein synthesis is blocked. When formed, DHX9 antiviral granules significantly reduced nascent protein synthesis in the MYXV-infected cancer cells. MYXV late gene transcription and translation were also significantly compromised, particularly in nonpermissive or semipermissive human cancer cells where MYXV replication is partly or completely restricted. Directed knockdown of DHX9 significantly enhanced viral late protein synthesis and progeny virus formation in normally restrictive cancer cells. We further demonstrate that DHX9 is not a component of the canonical cellular stress granules. DHX9 antiviral granules are induced by MYXV, and other poxviruses, in human cells and are associated with other known cellular components of stress granules, dsRNA and virus encoded dsRNA-binding protein M029, a known interactor with DHX9. Thus, DHX9 antiviral granules function by hijacking poxviral elements needed for the cytoplasmic viral replication factories. These results demonstrate a novel antiviral function for DHX9 that is recruited from the nucleus into the cytoplasm, and this step can be exploited to enhance oncolytic virotherapy against the subset of human cancer cells that normally restrict MYXV. IMPORTANCE The cellular DHX9 has both proviral and antiviral roles against diverse RNA and DNA viruses. In this article, we demonstrate that DHX9 can form unique antiviral granules in the cytoplasm during myxoma virus (MYXV) replication in human cancer cells. These antiviral granules sequester viral proteins and reduce viral late protein synthesis and thus regulate MYXV, and other poxviruses, that replicate in the cytoplasm. In addition, we show that in the absence of DHX9, the formation of DHX9 antiviral granules can be inhibited, which significantly enhanced oncolytic MYXV replication in human cancer cell lines where the virus is normally restricted. Our results also show that DHX9 antiviral granules are formed after viral infection but not by common nonviral cellular stress inducers. Thus, our study suggests that DHX9 has antiviral activity in human cancer cells, and this pathway can be targeted for enhanced activity of oncolytic poxviruses against even restrictive cancer cells.

Oncolytic Myxoma virus infects and damages the tegument of the human parasitic flatworm Schistosoma mansoni.

Schistosomiasis is a devastating disease caused by parasitic flatworms of the genus Schistosoma. Praziquantel (PZQ), the current treatment of choice, is ineffective against immature worms and cannot prevent reinfection. The continued reliance on a single drug for treatment increases the risk of the development of PZQ-resistant parasites. Reports of PZQ insusceptibility lends urgency to the need for new therapeutics. Here, we report that Myxoma virus (MYXV), an oncolytic pox virus which is non-pathogenic in all mammals except leporids, infects and replicates in S. mansoni schistosomula, juveniles, and adult male and female worms. MYXV infection results in the shredding of the tegument and reduced egg production in vitro, identifying MYXV as the first viral pathogen of schistosomes. MYXV is currently in preclinical studies to manage multiple human cancers, supporting its use in human therapeutics. Our findings raise the exciting possibility that MYXV virus represents a novel and safe class of potential anthelmintic therapeutics.

A Quadruplex qPCR for Detection and Differentiation of Classic and Natural Recombinant Myxoma Virus Strains of Leporids.

A natural recombinant myxoma virus (referred to as ha-MYXV or MYXV-Tol08/18) emerged in the Iberian hare (Lepus granatensis) and the European rabbit (Oryctolagus cuniculus) in late 2018 and mid-2020, respectively. This new virus is genetically distinct from classic myxoma virus (MYXV) strains that caused myxomatosis in rabbits until then, by acquiring an additional 2.8 Kbp insert within the m009L gene that disrupted it into ORFs m009L-a and m009L-b. To distinguish ha-MYXV from classic MYXV strains, we developed a robust qPCR multiplex technique that combines the amplification of the m000.5L/R duplicated gene, conserved in all myxoma virus strains including ha-MYXV, with the amplification of two other genes targeted by the real-time PCR systems designed during this study, specific either for classic MYXV or ha-MYXV strains. The first system targets the boundaries between ORFs m009L-a and m009L-b, only contiguous in classic strains, while the second amplifies a fragment within gene m060L, only present in recombinant MYXV strains. All amplification reactions were validated and normalized by a fourth PCR system directed to a housekeeping gene (18S rRNA) conserved in eukaryotic organisms, including hares and rabbits. The multiplex PCR (mPCR) technique described here was optimized for Taqman® and Evagreen® systems allowing the detection of as few as nine copies of viral DNA in the sample with an efficiency > 93%. This real-time multiplex is the first fast method available for the differential diagnosis between classic and recombinant MYXV strains, also allowing the detection of co-infections. The system proves to be an essential and effective tool for monitoring the geographical spread of ha-MYXV in the hare and wild rabbit populations, supporting the management of both species in the field.

Myxoma virus M013 protein antagonizes NF-κB and inflammasome pathways via distinct structural motifs.

Among the repertoire of immunoregulatory proteins encoded by myxoma virus, M013 is a viral homologue of the viral pyrin domain-only protein (vPOP) family. In myeloid cells, M013 protein has been shown to inhibit both the inflammasome and NF-κB signaling pathways by direct binding to ASC1 and NF-κB1, respectively. In this study, a three-dimensional homology model of the M013 pyrin domain (PYD) was built based on similarities to known PYD structures. A distinctive feature of the deduced surface electrostatic map of the M013 PYD is the presence of a negatively region consisting of numerous aspartate and glutamate residues in close proximity. Single-site mutations of aspartate and glutamate residues reveal their role in interactions with ASC-1. The biological significance of charge complementarity in the M013-ASC-1 interaction was further confirmed by functional assays of caspase-1 activation and subsequent secretion of cytokines. M013 also has a unique 33-residue C-terminal tail that follows the N-terminal PYD, and it is enriched in positively charged residues. Deletion of the tail of M013 significantly inhibited the interactions between M013 and NF-κB1, thus compromising the ability of the viral protein to suppress the secretion of pro-inflammatory cytokines. These results demonstrate that vPOP M013 exploits distinct structural motifs to regulate both the inflammasome and NF-κB pathways.

Punctuated Evolution of Myxoma Virus: Rapid and Disjunct Evolution of a Recent Viral Lineage in Australia.

Myxoma virus (MYXV) has been evolving in a novel host species-European rabbits-in Australia since 1950. Previous studies of viruses sampled from 1950 to 1999 revealed a remarkably clock-like evolutionary process across all Australian lineages of MYXV. Through an analysis of 49 newly generated MYXV genome sequences isolated in Australia between 2008 and 2017, we show that MYXV evolution in Australia can be characterized by three lineages, one of which exhibited a greatly elevated rate of evolutionary change and a dramatic breakdown of temporal structure. Phylogenetic analysis revealed that this apparently punctuated evolutionary event occurred between 1996 and 2012. The branch leading to the rapidly evolving lineage contained a relatively high number of nonsynonymous substitutions, and viruses in this lineage reversed a mutation found in the progenitor standard laboratory strain (SLS) and all previous sequences that disrupts the reading frame of the M005L/R gene. Analysis of genes encoding proteins involved in DNA synthesis or RNA transcription did not reveal any mutations likely to cause rapid evolution. Although there was some evidence for recombination across the MYXV phylogeny, this was not associated with the increase in the evolutionary rate. The period from 1996 to 2012 saw significant declines in wild rabbit numbers, due to the introduction of rabbit hemorrhagic disease and prolonged drought in southeastern Australia, followed by the partial recovery of populations. It is therefore possible that a rapidly changing environment for virus transmission changed the selection pressures faced by MYXV, altering the course and pace of virus evolution.IMPORTANCE The coevolution of myxoma virus (MYXV) and European rabbits in Australia is one of the most important natural experiments in evolutionary biology, providing insights into virus adaptation to new hosts and the evolution of virulence. Previous studies of MYXV evolution have also shown that the virus evolves both relatively rapidly and in a strongly clock-like manner. Using newly acquired MYXV genome sequences from Australia, we show that the virus has experienced a dramatic change in evolutionary behavior over the last 20 years, with a breakdown in clock-like structure, the appearance of a rapidly evolving virus lineage, and the accumulation of multiple nonsynonymous and indel mutations. We suggest that this punctuated evolutionary event may reflect a change in selection pressures as rabbit numbers declined following the introduction of rabbit hemorrhagic disease virus and drought in the geographic regions inhabited by rabbits.

Myxoma Virus M083 Is a Virulence Factor Which Mediates Systemic Dissemination.

Poxviruses are large, DNA viruses whose protein capsid is surrounded by one or more lipid envelopes. Embedded into these lipid envelopes are three conserved viral proteins which are thought to mediate binding of virions to target cells. While the function of these proteins has been studied in vitro, their specific roles during the pathogenesis of poxviral disease remain largely unclear. Here we present data demonstrating that the putative chondroitin binding protein M083 from the leporipoxvirus myxoma virus is a significant virulence factor during infection of susceptible Oryctolagus rabbits. Removal of M083 results in a reduced capacity of virus to spread beyond the regional lymph nodes and completely eliminates infection-mediated mortality. In vitro, removal of M083 results in only minor intracellular replication defects but causes a significant reduction in the ability of myxoma virus to spread from infected epithelial cells onto primary lymphocytes. We hypothesize that the physiological role of M083 is therefore to mediate the spread of myxoma virus onto rabbit lymphocytes, allowing these cells to disseminate virus throughout infected rabbits.IMPORTANCE Poxviruses represent both a class of human pathogens and potential therapeutic agents for the treatment of human malignancy. Understanding the basic biology of these agents is therefore significant to human health in a variety of ways. While the mechanisms mediating poxviral binding have been well studied in vitro, how these mechanisms impact poxviral pathogenesis in vivo remains unclear. The current study advances our understanding of how poxviral binding impacts viral pathogenesis by demonstrating that the putative chondroitin binding protein M083 plays a critical role during the pathogenesis of myxoma virus in susceptible Oryctolagus rabbits by impacting viral dissemination through changes in the transfer of virions onto primary splenocytes.

RNA Polymerase Mutations Selected during Experimental Evolution Enhance Replication of a Hybrid Vaccinia Virus with an Intermediate Transcription Factor Subunit Replaced by the Myxoma Virus Ortholog.

High-throughput DNA sequencing enables the study of experimental evolution in near real time. Until now, mutants with deletions of nonessential host range genes were used in experimental evolution of vaccinia virus (VACV). Here, we guided the selection of adaptive mutations that enhanced the fitness of a hybrid virus in which an essential gene had been replaced with an ortholog from another poxvirus genus. Poxviruses encode a complete system for transcription, including RNA polymerase and stage-specific transcription factors. The abilities of orthologous intermediate transcription factors from other poxviruses to substitute for those of VACV, as determined by transfection assays, corresponded with the degree of amino acid identity. VACV in which the A8 or A23 intermediate transcription factor subunit gene was replaced by the myxoma (MYX) virus ortholog exhibited decreased replication. During three parallel serial passages of the hybrid virus with the MYXA8 gene, plaque sizes and virus yields increased. DNA sequencing of virus populations at passage 10 revealed high frequencies of five different single nucleotide mutations in the two largest RNA polymerase subunits, RPO147 and RPO132, and two different Kozak consensus sequence mutations predicted to increase translation of the MYXA8 mRNA. Surprisingly, there were no mutations within either intermediate transcription factor subunit. Based on homology with Saccharomyces cerevisiae RNA polymerase, the VACV mutations were predicted to be buried within the internal structure of the enzyme. By directly introducing single nucleotide substitutions into the genome of the original hybrid virus, we demonstrated that both RNA polymerase and translation-enhancing mutations increased virus replication independently.IMPORTANCE Previous studies demonstrated the experimental evolution of vaccinia virus (VACV) following deletion of a host range gene important for evasion of host immune defenses. We have extended experimental evolution to essential genes that cannot be deleted but could be replaced by a divergent orthologous gene from another poxvirus. Replacement of a VACV transcription factor gene with one from a distantly related poxvirus led to decreased fitness as evidenced by diminished replication. Serially passaging the hybrid virus at a low multiplicity of infection provided conditions for selection of adaptive mutations that improved replication. Notably, these included five independent mutations of the largest and second largest RNA polymerase subunits. This approach should be generally applicable for investigating adaptation to swapping of orthologous genes encoding additional essential proteins of poxviruses as well as other viruses.

Myxoma virus M156 is a specific inhibitor of rabbit PKR but contains a loss-of-function mutation in Australian virus isolates.

Myxoma virus (MYXV) is a rabbit-specific poxvirus, which is highly virulent in European rabbits. The attenuation of MYXV and the increased resistance of rabbits following the release of MYXV in Australia is one of the best-documented examples of host-pathogen coevolution. To elucidate the molecular mechanisms that contribute to the restriction of MYXV infection to rabbits and MYXV attenuation in the field, we have studied the interaction of the MYXV protein M156 with the host antiviral protein kinase R (PKR). In yeast and cell-culture transfection assays, M156 only inhibited rabbit PKR but not PKR from other tested mammalian species. Infection assays with human HeLa PKR knock-down cells, which were stably transfected with human or rabbit PKR, revealed that only human but not rabbit PKR was able to restrict MYXV infection, whereas both PKRs were able to restrict replication of a vaccinia virus (VACV) strain that lacks the PKR inhibitors E3 and K3. Inactivation of M156R led to MYXV virus attenuation in rabbit cells, which was rescued by the ectopic expression of VACV E3 and K3. We further show that a mutation in the M156 encoding gene that was identified in more than 50% of MYXV field isolates from Australia resulted in an M156 variant that lost its ability to inhibit rabbit PKR and led to virus attenuation. The species-specific inhibition of rabbit PKR by M156 and the M156 loss-of-function in Australian MYXV field isolates might thus contribute to the species specificity of MYXV and to the attenuation in the field, respectively.

Genetic Variability of Myxoma Virus Genomes.

Myxomatosis is a recurrent problem on rabbit farms throughout Europe despite the success of vaccines. To identify gene variations of field and vaccine strains that may be responsible for changes in virulence, immunomodulation, and immunoprotection, the genomes of 6 myxoma virus (MYXV) strains were sequenced: German field isolates Munich-1, FLI-H, 2604, and 3207; vaccine strain MAV; and challenge strain ZA. The analyzed genomes ranged from 147.6 kb (strain MAV) to 161.8 kb (strain 3207). All sequences were affected by several mutations, covering 24 to 93 open reading frames (ORFs) and resulted in amino acid substitutions, insertions, or deletions. Only strains Munich-1 and MAV revealed the deletion of 10 ORFs (M007L to M015L) and 11 ORFs (M007L to M008.1L and M149R to M008.1R), respectively. Major differences were observed in the 27 immunomodulatory proteins encoded by MYXV. Compared to the reference strain Lausanne, strains FLI-H, 2604, 3207, and ZA showed the highest amino acid identity (>98.4%). In strains Munich-1 and MAV, deletion of 5 and 10 ORFs, respectively, was observed, encoding immunomodulatory proteins with ankyrin repeats or members of the family of serine protease inhibitors. Furthermore, putative immunodominant surface proteins with homology to vaccinia virus (VACV) were investigated in the sequenced strains. Only strain MAV revealed above-average frequencies of amino acid substitutions and frameshift mutations. Finally, we performed recombination analysis and found signs of recombination in vaccine strain MAV. Phylogenetic analysis showed a close relationship of strain MAV and the MSW strain of Californian MYXV. However, in a challenge model, strain MAV provided full protection against lethal challenges with strain ZA. IMPORTANCE: Myxoma virus (MYXV) is pathogenic for European rabbits and two North American species. Due to sophisticated strategies in immune evasion and oncolysis, MYXV is an important model virus for immunological and pathological research. In its natural hosts, MYXV causes a benign infection, whereas in European rabbits, it causes the lethal disease myxomatosis. Since the introduction of MYXV into Australia and Europe for the biological control of European rabbits in the 1950s, a coevolution of host and pathogen has started, selecting for attenuated virus strains and increased resistance in rabbits. Evolution of viruses is a continuous process and influences the protective potential of vaccines. In our analyses, we sequenced 6 MYXV field, challenge, and vaccine strains. We focused on genes encoding proteins involved in virulence, host range, immunomodulation, and envelope composition. Genes affected most by mutations play a role in immunomodulation. However, attenuation cannot be linked to individual mutations or gene disruptions.

Reverse Engineering Field Isolates of Myxoma Virus Demonstrates that Some Gene Disruptions or Losses of Function Do Not Explain Virulence Changes Observed in the Field.

The coevolution of myxoma virus (MYXV) and wild European rabbits in Australia and Europe is a paradigm for the evolution of a pathogen in a new host species. Genomic analyses have identified the mutations that have characterized this evolutionary process, but defining causal mutations in the pathways from virulence to attenuation and back to virulence has not been possible. Using reverse genetics, we examined the roles of six selected mutations found in Australian field isolates of MYXV that fall in known or potential virulence genes. Several of these mutations occurred in genes previously identified as virulence genes in whole-gene knockout studies. Strikingly, no single or double mutation among the mutations tested had an appreciable impact on virulence. This suggests either that virulence evolution was defined by amino acid changes other than those analyzed here or that combinations of multiple mutations, possibly involving epistatic interactions or noncoding sequences, have been critical in the ongoing evolution of MYXV virulence. In sum, our results show that single-gene knockout studies of a progenitor virus can have little power to predict the impact of individual mutations seen in the field. The genetic determinants responsible for this canonical case of virulence evolution remain to be determined.IMPORTANCE The species jump of myxoma virus (MYXV) from the South American tapeti to the European rabbit populations of Australia and Europe is a canonical example of host-pathogen coevolution. Detailed molecular studies have identified multiple genes in MYXV that are critical for virulence, and genome sequencing has revealed the evolutionary history of MYXV in Australia and Europe. However, it has not been possible to categorically identify the key mutations responsible for the attenuation of or reversion to virulence during this evolutionary process. Here we use reverse genetics to examine the role of mutations in viruses isolated early and late in the Australian radiation of MYXV. Surprisingly, none of the candidate mutations that we identified as likely having roles in attenuation proved to be important for virulence. This indicates that considerable caution is warranted when interpreting the possible role of individual mutations during virulence evolution.

Disentangling synergistic disease dynamics: Implications for the viral biocontrol of rabbits.

European rabbits (Oryctolagus cuniculus) have been exposed to rabbit haemorrhagic disease virus (RHDV) and myxoma virus (MYXV) in their native and invasive ranges for decades. Yet, the long-term effects of these viruses on rabbit population dynamics remain poorly understood. In this context, we analysed 17 years of detailed capture-mark-recapture data (2000-2016) from Turretfield, South Australia, using a probabilistic state-space hierarchical modelling framework to estimate rabbit survival and epidemiological dynamics. While RHDV infection and disease-induced death were most prominent during annual epidemics in winter and spring, we found evidence for continuous infection of susceptible individuals with RHDV throughout the year. RHDV-susceptible rabbits had, on average, 25% lower monthly survival rates compared to immune individuals, while the average monthly force of infection in winter and spring was ~38%. These combined to result in an average infection-induced mortality rate of 69% in winter and spring. Individuals susceptible to MYXV and immune to RHDV had similar survival probabilities to those having survived infections from both viruses, whereas individuals susceptible to both RHDV and MYXV had higher survival probabilities than those susceptible to RHDV and immune to MYXV. This suggests that MYXV may reduce the future survival rates of individuals that endure initial MYXV infection. There was no evidence for long-term changes in disease-induced mortality and infection rates for either RHDV or MYXV. We conclude that continuous, year-round virus perpetuation (and perhaps heterogeneity in modes of transmission and infectious doses during and after epidemics) acts to reduce the efficiency of RHDV and MYXV as biocontrol agents of rabbits in their invasive range. However, if virulence can be maintained as relatively constant through time, RHDV and MYXV will likely continue realizing strong benefits as biocontrol agents.

Ex vivo virotherapy with myxoma virus does not impair hematopoietic stem and progenitor cells.

BACKGROUND: Relapsing disease is a major challenge after hematopoietic cell transplantation for hematological malignancies. Myxoma virus (MYXV) is an oncolytic virus that can target and eliminate contaminating cancer cells from auto-transplant grafts. The aims of this study were to examine the impact of MYXV on normal hematopoietic stem and progenitor cells and define the optimal treatment conditions for ex vivo virotherapy. METHODS: Bone marrow (BM) and mobilized peripheral blood stem cells (mPBSCs) from patients with hematologic malignancies were treated with MYXV at various time, temperature and incubation media conditions. Treated BM cells from healthy normal donors were evaluated using flow cytometry for MYXV infection, long-term culture-initiating cell (LTC-IC) assay and colony-forming cell (CFC) assay. RESULTS: MYXV initiated infection in up to 45% of antigen-presenting monocytes, B cells and natural killer cells; however, these infections were uniformly aborted in >95% of all cells. Fresh graft sources showed higher levels of MYXV infection initiation than cryopreserved specimens, but in all cases less than 10% of CD34(+) cells could be infected after ex vivo MYXV treatment. MYXV did not impair LTC-IC colony numbers compared with mock treatment. CFC colony types and numbers were also not impaired by MYXV treatment. MYXV incubation time, temperature or culture media did not significantly change the percentage of infected cells, LTC-IC colony formation or CFC colony formation. CONCLUSIONS: Human hematopoietic cells are non-permissive for MYXV. Human hematopoietic stem and progenitor cells were not infected and thus unaffected by MYXV ex vivo treatment.

Prevention of EBV lymphoma development by oncolytic myxoma virus in a murine xenograft model of post-transplant lymphoproliferative disease.

Epstein-Barr virus (EBV) has been associated with a variety of epithelial and hematologic malignancies, including B-, T- and NK cell-lymphomas, Hodgkin's disease (HD), post-transplant lymphoproliferative diseases (LPDs), nasopharyngeal and gastric carcinomas, smooth muscle tumors, and HIV-associated lymphomas. Currently, treatment options for EBV-associated malignancies are limited. We have previously shown that myxoma virus specifically targets various human solid tumors and leukemia cells in a variety of animal models, while sparing normal human or murine tissues. Since transplant recipients of bone marrow or solid organs often develop EBV-associated post-transplant LPDs and lymphoma, myxoma virus may be of utility to prevent EBV-associated malignancies in immunocompromised transplant patients where treatment options are frequently limited. In this report, we demonstrate the safety and efficacy of myxoma virus purging as a prophylactic strategy for preventing post-transplant EBV-transformed human lymphomas, using a highly immunosuppressed mouse xenotransplantation model. This provides support for developing myxoma virus as a potential oncolytic therapy for preventing EBV-associated LPDs following transplantation of bone marrow or solid organ allografts.

Epidemiologic, clinicopathologic, and diagnostic findings in pet rabbits with myxomatosis caused by the California MSW strain of myxoma virus: 11 cases (2022-2023).

OBJECTIVE: To determine epidemiologic features of naturally occurring myxomatosis in domestic rabbits in California and to characterize clinicopathologic and diagnostic findings. ANIMALS: 11 client-owned rabbits, Oryctolagus cuniculus subsp domesticus. CLINICAL PRESENTATION: A prospective study of pet rabbits with myxomatosis seen at an exotic animal specialty clinic in Santa Cruz county, California, was conducted between January 1, 2022, and December 31, 2023. Rabbits were included in the study if they had bilateral blepharedema and were PCR positive for myxoma virus. RESULTS: All infected rabbits had spent time outdoors. Common clinical signs included bilateral blepharedema (11/11), anogenital edema (10/11), rectal temperature ≥ 39.7 °C (5/9), and sudden death (4/11). Eyelid biopsies from all rabbits (11/11) were positive for myxoma virus by qualitative PCR followed by Sanger sequencing (100% nucleotide identity to strain MSW, also known as California/San Francisco 1950 [Genbank accession KF148065]). Most rabbits had keratinocytes containing eosinophilic intracytoplasmic viral inclusions in biopsies of edematous skin (8/11) and lymphocyte necrosis in the spleen (10/11). Immunohistochemistry identified myxoma virus in samples of skin, heart, lung, ileum, spleen, and lymph node. CLINICAL RELEVANCE: Clinical signs of myxomatosis caused by the MSW strain of myxoma virus are distinctive but subtle. Cases occur regularly in the Santa Cruz and San Jose regions of California. As infection with this virus is almost 100% fatal and no vaccine is available in the US, owners of domestic rabbits in endemic areas should keep their pets indoors or behind mosquito screens. Myxomatosis is a reportable disease in the US, and the appropriate state or federal agencies should be contacted when outbreaks occur.