Evolutionary history and attenuation of myxoma virus on two continents

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.

Identification and Characterisation of a Myxoma Virus Detected in the Italian Hare (Lepus corsicanus)

Myxoma virus (MYXV) is a Leporipoxvirus (genus) belonging to the family Poxviridae; it is characterised by a genome of approximately 161 kb dsDNA encoding for several proteins that play an essential role in both host spectrum determination and immunomodulation. The healthy reservoir of the virus is Sylvilagus spp. At the same time, in wild and domestic European rabbits (Oryctolagus cuniculus), MYXV is the etiologic agent of myxomatosis, a disease with an extremely high mortality rate. In 2014, an interspecies jump of MYXV was reported in Lepus europaeus in the UK. In 2018, myxomatosis induced by a new recombinant strain called MYXV-To was identified during a large outbreak in Iberian hares (Lepus granatensis) in Spain. Here, we describe the case of myxomatosis in another hare species: an adult male Italian hare (Lepus corsicanus) found dead in 2018 in Sicily with lesions suggestive of myxomatosis and treponema infection. Laboratory tests, e.g., end-point PCR and negative staining electron microscopy, confirmed the presence of both pathogens. MYXV was then isolated from tissue samples in permissive cells and sequenced using NGS technology. Main genomic differences concerning known MYXV strains are discussed.

Darwin review: the evolution of virulence in human pathogens

By definition, all pathogens cause some level of harm to their hosts. If this harm occurs while the pathogen is transmitting, it can negatively affect the pathogen's fitness by shortening the duration over which transmission can occur. However, many of the factors that increase virulence (i.e. harm to host) also promote transmission, driving the pathogen population towards an optimal state of intermediate virulence. A wider spectrum of virulence may be maintained among pathogen populations which are structured into multiple discrete strains though direct resource and immune-mediated competition. These various evolutionary outcomes, and the effects of medical and public health interventions, are best understood within a framework that recognizes the complex relationship between transmission and virulence in the context of the antigenic diversity of the pathogen population.

Evolution of pathogen virulence: Studying the complex interplay of pathogen interactions, virulence and transmission helps us understand how they evolve and spread

Navigating the intricate interplay of competitive and co-operative interactions and the complex relationship between virulence and transmission pose challenges for our understanding of how pathogens evolve and spread.

Lumpy Skin Disease Virus Genome Sequence Analysis: Putative Spatio-Temporal Epidemiology, Single Gene versus Whole Genome Phylogeny and Genomic Evolution

Lumpy Skin Disease virus is a poxvirus from the genus Capripox that mainly affects bovines and it causes severe economic losses to livestock holders. The Lumpy Skin Disease virus is currently dispersing in Asia, but little is known about detailed phylogenetic relations between the strains and genome evolution. We reconstructed a whole-genome-sequence (WGS)-based phylogeny and compared it with single-gene-based phylogenies. To study population and spatiotemporal patterns in greater detail, we reconstructed networks. We determined that there are strains from multiple clades within the previously defined cluster 1.2 that correspond with recorded outbreaks across Eurasia and South Asia (Indian subcontinent), while strains from cluster 2.5 spread in Southeast Asia. We concluded that using only a single gene (cheap, fast and easy to routinely use) for sequencing lacks phylogenetic and spatiotemporal resolution and we recommend to create at least one WGS whenever possible. We also found that there are three gene regions, highly variable, across the genome of LSDV. These gene regions are located in the 5' and 3' flanking regions of the LSDV genome and they encode genes that are involved in immune evasion strategies of the virus. These may provide a starting point to further investigate the evolution of the virus.

Sequencing of monkeypox virus from infected patients reveals viral genomes with APOBEC3-like editing, gene inactivation, and bacterial agents of skin superinfection

A large outbreak of Monkeypox virus (MPXV) infections has arisen in May 2022 in nonendemic countries. Here, we performed DNA metagenomics using next-generation sequencing with Illumina or Nanopore technologies for clinical samples from MPXV-infected patients diagnosed between June and July 2022. Classification of the MPXV genomes and determination of their mutational patterns were performed using Nextclade. Twenty-five samples from 25 patients were studied. A MPXV genome was obtained for 18 patients, essentially from skin lesions and rectal swabbing. All 18 genomes were classified in clade IIb, lineage B.1, and we identified four B.1 sublineages (B.1.1, B.1.10, B.1.12, B.1.14). We detected a high number of mutations (range, 64-73) relatively to a 2018 Nigerian genome (genome GenBank Accession no. NC_063383.1), which were harbored by a large part of a set of 3184 MPXV genomes of lineage B.1 recovered from GenBank and Nextstrain; and we detected 35 mutations relatively to genome ON563414.3 (a B.1 lineage reference genome). Nonsynonymous mutations occurred in genes encoding central proteins, among which transcription factors and core and envelope proteins, and included two mutations that would truncate a RNA polymerase subunit and a phospholipase d-like protein, suggesting an alternative start codon and gene inactivation, respectively. A large majority (94%) of nucleotide substitutions were G > A or C > U, suggesting the action of human APOBEC3 enzymes. Finally, >1000 reads were identified as from Staphylococcus aureus and Streptococcus pyogenes for 3 and 6 samples, respectively. These findings warrant a close genomic monitoring of MPXV to get a better picture of the genetic micro-evolution and mutational patterns of this virus, and a close clinical monitoring of skin bacterial superinfection in monkeypox patients.

Molecular Mechanisms of Poxvirus Evolution

Poxviruses are often thought to evolve relatively slowly because they are double-stranded DNA pathogens with proofreading polymerases. However, poxviruses have highly adaptable genomes and can undergo relatively rapid genotypic and phenotypic change, as illustrated by the recent increase in human-to-human transmission of monkeypox virus. Advances in deep sequencing technologies have demonstrated standing nucleotide variation in poxvirus populations, which has been underappreciated. There is also an emerging understanding of the role genomic architectural changes play in shaping poxvirus evolution. These mechanisms include homologous and nonhomologous recombination, gene duplications, gene loss, and the acquisition of new genes through horizontal gene transfer. In this review, we discuss these evolutionary mechanisms and their potential roles for adaption to novel host species and modulating virulence.

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.

The Ecology of Viral Emergence

The coronavirus disease 2019 (COVID-19) pandemic has had a profound impact on human health, economic well-being, and societal function. It is essential that we use this generational experience to better understand the processes that underpin the emergence of COVID-19 and other zoonotic diseases. Herein, I review the mechanisms that determine why and how viruses emerge in new hosts, as well as the barriers to this process. I show that traditional studies of virus emergence have an inherent anthropocentric bias, with disease in humans considered the inevitable outcome of virus emergence, when in reality viruses are integral components of a global ecosystem characterized by continual host jumping with humans also transmitting their viruses to other animals. I illustrate these points using coronaviruses, including severe acute respiratory syndrome coronavirus 2, as a case study. I also outline the potential steps that can be followed to help mitigate and prevent future pandemics, with combating climate change a central component. Expected final online publication date for the Annual Review of Virology, Volume 9 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Identification of a Novel Myxoma Virus C7-Like Host Range Factor That Enabled a Species Leap from Rabbits to Hares

Myxoma virus (MYXV) is naturally found in rabbit Sylvilagus species and is known to cause lethal myxomatosis in European rabbits (Oryctolagus cuniculus). In 2019, an MYXV strain (MYXV strain Toledo [MYXV-Tol]) causing myxomatosis-like disease in Iberian hares (Lepus granatensis) was identified. MYXV-Tol acquired a recombinant region of ∼2.8 kb harboring several new genes, including a novel host range gene (M159) that we show to be an orthologous member of the vaccinia virus C7 host range family. Here, to test whether M159 alone has enabled MYXV to alter its host range to Iberian hares, several recombinant viruses were generated, including an MYXV-Tol ΔM159 (knockout) strain. While MYXV-Tol underwent fully productive infection in hare HN-R cells, neither the wild-type MYXV-Lau strain (lacking M159) nor vMyxTol-ΔM159 (deleted for M159) was able to infect and replicate, showing that the ability of MYXV-Tol to infect these cells and replicate depends on the presence of M159. Similar to other C7L family members, M159 was shown to be expressed as an early/late gene but was translocated into the nucleus at later time points, indicating that further studies are needed to elucidate its role in the nucleus. Finally, in rabbit cells, the M159 protein did not contribute to increased replication but was able to upregulate the replication levels of MYXV in nonpermissive and semipermissive human cancer cells, suggesting that the M159-targeted pathway is conserved across mammalian species. Altogether, these observations demonstrate that the M159 protein plays a critical role in determining the host specificity of MYXV-Tol in hare and human cells by imparting new host range functions. IMPORTANCE The coevolution of European rabbit populations and MYXV is a textbook example of an arms race between a pathogen and a host. Recently, a recombinant MYXV (MYXV-Tol) crossed the species barrier by jumping from leporid species to another species, causing lethal myxomatosis-like disease. Given the highly pathogenic nature of this new virus in hares and the incidences of other poxvirus cross-species spillovers into other animals, including humans, it is important to understand how and why MYXV-Tol was able to become virulent in a new host species. The results presented clearly demonstrate that M159 is the key factor allowing MYXV-Tol replication in hare cells by imparting new host range functions. These results have the potential to improve current knowledge about the virulence of poxviruses and provide a platform to better understand the new MYXV-Tol, rendering the virus capable of leaping into a new host species.

Viral haemorrhagic disease: RHDV type 2 ten years later

Until the early 1980s, it was totally unknown that lagomorphs were the hosts of several caliciviruses, which were included in the genus Lagovirus by the International Committee on Taxonomy of Viruses (ICTV) in 2000. In those years, two new diseases appeared, with very similar clinical and pathological profiles and associated high mortality rates: rabbit haemorrhagic disease (RHD) in rabbits and European Brown Hare Syndrome (EBHS) in European brown hares. It took a few years to ascertain that both diseases, actually acute and fatal forms of hepatitis, were caused by two genetically related caliciviruses, but they were finally classified by ICTV into two distinct viral species on the basis of their molecular characterisation and epidemiological data: RHDV in rabbit and EBHSV in brown hare. RHD has had a devastating effect on rabbit farms, causing great economic damage, especially in China, where RHD was first noticed around 1982, and in Europe. RHD has also severely affected wild rabbit populations, whose drastic decline has caused serious ecological imbalances in territories such as Spain, where rabbits are a central link in the wildlife food chain. Since the early 1990s, with the increased availability on the market of RHDV vaccines effective in protecting rabbits from RHD, the impact of the disease on rabbit farms has been significantly reduced. In the following years, also considering that RHDV is an endemic virus that cannot be eradicated, farmers learned how to manage the continuous use of RHDV vaccine in relation to the epidemiological situation, the type of breeding farm and the costs of vaccination prophylaxis. Although precarious, management of the RHD risk for rabbit farmers reached an acceptable equilibrium, which was, however, completely upset starting from 2010 by the emergence of another lagovirus also causing RHD. The genome of the newly emerged virus shows limited differences from that of RHDV, but the phenotypic traits of the two viruses are distinctive in at least three main respects: 1) The antigenic profile of the virus (the “face” of the virus recognised by the antibodies) is largely different from that of RHDV. 2) Newborn rabbits only a couple of weeks old die of RHD when infected with the new virus, while RHDV infections run asymptomatic until 7-8 wk of age. 3) The new virus, which started in Europe, has spread over the years to several continents, affecting wild and/or domestic rabbit populations. During this worldwide distribution, the new virus infected several lagomorph species and was shown to cause RHD in most of them. Considering these marked differences and the fact that the new virus is not a variant of RHDV, we proposed the name RHDV type 2 (RHDV2). All these main distinctive traits that differentiate RHDV from RHDV2 have the following consequences in practice: 1) The antigenic difference between RHDV and RHDV2 (their ‘faces’) is so great that we need “new” specific vaccines to control RHDV2 (i.e. RHDV2 is a new serotype). 2) In the event of an RHDV2 infection in suckling rabbits, the presence of maternal antibodies to RHDV2 in the blood is the only way to prevent RHD. In contrast, newborns are naturally resistant to RHD if infected with RHDV and therefore, in terms of protection, the presence of maternal antibodies is useless. 3) When RHD outbreaks occur in territories where rabbits live in sympatry with populations of other lagomorphs, viral contamination in the environment reaches sufficiently high levels to facilitate the transmission of RHDV2 to other lagomorphs, including those with a lower susceptibility to infection than the rabbit. Taken together, these phenotypic traits characteristic of RHDV2 are the reason for its rapid spread across the territory and the concomitant disappearance of RHDV. Probably the most striking example of the epidemiological consequences related to the peculiar features of RHDV2 is its rapid spread in the USA and Mexico, where it is now practically endemic. There, despite repeated isolated outbreaks of RHD caused by RHDV from 2000 onwards in small rabbit farms, RHDV has never been able to become endemic.

Genome Plasticity of African Swine Fever Virus: Implications for Diagnostics and Live-Attenuated Vaccines

African swine fever (ASF) is a highly contagious transboundary viral hemorrhagic disease of domestic and wild pigs presenting a significant threat to the global swine industry. Following its introduction in Caucasus, Georgia, in 2007, the genome of the genotype II of African swine fever virus (ASFV) strain Georgia-07 and its derivatives accumulated significant mutations, resulting in the emergence of genetic variants within short epidemiological timescales as it spreads and infects different hosts in diverse ecosystems, causing outbreaks in Europe, South Asia, South East Asia and the Caribbean. This suggests that ASFV, with a comparatively large and complex DNA genome, is susceptible to genetic mutations including deletions and that although the virus is environmentally stable, it is genetically unstable. This has implications for the development of vaccines and diagnostic tests for disease detection and surveillance. Analysis of the ASFV genome revealed recombination hotspots, which in double-stranded DNA (dsDNA) viruses represent key drivers of genetic diversity. The ability of pox virus, a dsDNA virus with a genome complexity similar to ASFV, regaining virulence following the deletion of a virulence gene via gene amplification, coupled with the recent emergence and spread of live-attenuated ASFV vaccine strains causing disease and death in pigs in China, raise legitimate concerns around the use of live-attenuated ASFV vaccines in non-endemic regions to control the potential introduction. Further research into the risk of using live-attenuated ASFV in non-endemic regions is highly needed.

The three Ts of virulence evolution during zoonotic emergence

There is increasing interest in the role that evolution may play in current and future pandemics, but there is often also considerable confusion about the actual evolutionary predictions. This may be, in part, due to a historical separation of evolutionary and medical fields, but there is a large, somewhat nuanced body of evidence-supported theory on the evolution of infectious disease. In this review, we synthesize this evolutionary theory in order to provide a framework for clearer understanding of the key principles. Specifically, we discuss the selection acting on zoonotic pathogens' transmission rates and virulence at spillover and during emergence. We explain how the direction and strength of selection during epidemics of emerging zoonotic disease can be understood by a three Ts framework: trade-offs, transmission, and time scales. Virulence and transmission rate may trade-off, but transmission rate is likely to be favoured by selection early in emergence, particularly if maladapted zoonotic pathogens have 'no-cost' transmission rate improving mutations available to them. Additionally, the optimal virulence and transmission rates can shift with the time scale of the epidemic. Predicting pathogen evolution, therefore, depends on understanding both the trade-offs of transmission-improving mutations and the time scales of selection.

Virus shedding kinetics and unconventional virulence tradeoffs

Tradeoff theory, which postulates that virulence provides both transmission costs and benefits for pathogens, has become widely adopted by the scientific community. Although theoretical literature exploring virulence-tradeoffs is vast, empirical studies validating various assumptions still remain sparse. In particular, truncation of transmission duration as a cost of virulence has been difficult to quantify with robust controlled in vivo studies. We sought to fill this knowledge gap by investigating how transmission rate and duration were associated with virulence for infectious hematopoietic necrosis virus (IHNV) in rainbow trout (Oncorhynchus mykiss). Using host mortality to quantify virulence and viral shedding to quantify transmission, we found that IHNV did not conform to classical tradeoff theory. More virulent genotypes of the virus were found to have longer transmission durations due to lower recovery rates of infected hosts, but the relationship was not saturating as assumed by tradeoff theory. Furthermore, the impact of host mortality on limiting transmission duration was minimal and greatly outweighed by recovery. Transmission rate differences between high and low virulence genotypes were also small and inconsistent. Ultimately, more virulent genotypes were found to have the overall fitness advantage, and there was no apparent constraint on the evolution of increased virulence for IHNV. However, using a mathematical model parameterized with experimental data, it was found that host culling resurrected the virulence tradeoff and provided low virulence genotypes with the advantage. Human-induced or natural culling, as well as host population fragmentation, may be some of the mechanisms by which virulence diversity is maintained in nature. This work highlights the importance of considering non-classical virulence tradeoffs.

Evaluating alternative hypotheses to explain the downward trend in the indices of the COVID-19 pandemic death rate

BACKGROUND

In the ongoing Covid-19 pandemic, in the global data on the case fatality ratio (CFR) and other indices reflecting death rate, there is a consistent downward trend from mid-April to mid-November. The downward trend can be an illusion caused by biases and limitations of data or it could faithfully reflect a declining death rate. A variety of explanations for this trend are possible, but a systematic analysis of the testable predictions of the alternative hypotheses has not yet been attempted.

METHODOLOGY

We state six testable alternative hypotheses, analyze their testable predictions using public domain data and evaluate their relative contributions to the downward trend.

RESULTS

We show that a decline in the death rate is real; changing age structure of the infected population and evolution of the virus towards reduced virulence are the most supported hypotheses and together contribute to major part of the trend. The testable predictions from other explanations including altered testing efficiency, time lag, improved treatment protocols and herd immunity are not consistently supported, or do not appear to make a major contribution to this trend although they may influence some other patterns of the epidemic.

CONCLUSION

The fatality of the infection showed a robust declining time trend between mid April to mid November. Changing age class of the infected and decreasing virulence of the pathogen were found to be the strongest contributors to the trend.

Viruses for Landscape-Scale Therapy: Biological Control of Rabbits in Australia

Viral diseases, whether of animals or humans, are normally considered as problems to be managed. However, in Australia, two viruses have been used as landscape-scale therapeutics to control European rabbits (Oryctolagus cuniculus), the preeminent invasive vertebrate pest species. Rabbits have caused major environmental and agricultural losses and contributed to extinction of native species. It was not until the introduction of Myxoma virus that effective control of this pest was obtained at a continental scale. Subsequent coevolution of rabbit and virus saw a gradual reduction in the effectiveness of biological control that was partially ameliorated by the introduction of the European rabbit flea to act as an additional vector for the virus. In 1995, a completely different virus, Rabbit hemorrhagic disease virus (RHDV), escaped from testing and spread through the Australian rabbit population and again significantly reduced rabbit numbers and environmental impacts. The evolutionary pressures on this virus appear to be producing quite different outcomes to those that occurred with myxoma virus and the emergence and invasion of a novel genotype of RHDV in 2014 have further augmented control. Molecular studies on myxoma virus have demonstrated multiple proteins that manipulate the host innate and adaptive immune response; however the molecular basis of virus attenuation and reversion to virulence are not yet understood.

Analysis of the evolution and transmission dynamics of the field MDV in China during the years 1995-2020, indicating the emergence of a unique cluster with the molecular characteristics of vv+ MDV that has become endemic in southern China

Marek's disease (MD) continues to threaten the sustainability of the world poultry industry. In this study, the sequences of the meq gene of 220 MDV strains isolated during the years 1964-2020 were analysed, including 50 from our group plus 170 isolates from the GenBank. Analyses, using phylogenetic trees, amino acid (aa)-mutation screening, evolutionary studies and transmission dynamics were all performed. All strains were divided into two clusters (Clusters 1 and 2), and Cluster 1 includes the mild strains, the vaccine strains and the foreign virulent strains, while Cluster 2 was dominated by the Chinese field strains. Our study identified that the Chinese field strains in Cluster 2 during the years 1995-2020 likely originated in the 1980s from abroad, and the estimated genetic diversity of these strains experienced two growth phases in the years 2005-2007.5 and 2015-2017. Viral phylogeography identified 3 major geographic provincial regions for the Chinese field strains of Cluster 2: the Northeastern Region (Jilin, Liaoning and Heilongjiang), the East-central Region (Henan, Shandong and Jiangsu) and the Southern Region (Guangxi, Guangdong and Yunnan). The spread of Northeastern strains to East-central chicken flocks and the further spread from Guangxi to Guangdong are strongly indicated. The emergence of the mutations A88T and Q93R together in the Southern strains during the years 2017-2020 with molecular characteristics of vv+ MDV were also found later than those in the Northern strains. Overall, the Chinese field strains in Cluster 2 in southern China in recent years have been rapidly evolving. Guangxi Province has become an epicentre for these viruses and the chicken flocks in the Southern region have been facing the adverse effects of the emerging vv+ MDV.

RNA helicase A as co-factor for DNA viruses during replication

RNA helicase A (RHA) is a ubiquitously expressed DExH-box helicase enzyme that is involved in a wide range of biological processes including transcription, translation, and RNA processing. A number of RNA viruses recruit RHA to the viral RNA to facilitate virus replication. DNA viruses contain a DNA genome and replicate using a DNA-dependent DNA polymerase. RHA has also been reported to associate with some DNA viruses during replication, in which the enzyme acts on the viral RNA or protein products. As shown for Epstein-Barr virus and Kaposi's sarcoma-associated herpesvirus, RHA has potential to allow the virus to control a switch in cellular gene expression to modulate the antiviral response. While the study of the interaction of RHA with DNA viruses is still at an early stage, preliminary evidence indicates that the underlying molecular mechanisms are diverse. We now review the current status of this emerging field.

Myxoma Virus-Encoded Host Range Protein M029: A Multifunctional Antagonist Targeting Multiple Host Antiviral and Innate Immune Pathways

Myxoma virus (MYXV) is the prototypic member of the Leporipoxvirus genus of the Poxviridae family of viruses. In nature, MYXV is highly restricted to leporids and causes a lethal disease called myxomatosis only in European rabbits (Oryctologous cuniculus). However, MYXV has been shown to also productively infect various types of nonrabbit transformed and cancer cells in vitro and in vivo, whereas their normal somatic cell counterparts undergo abortive infections. This selective tropism of MYXV for cancer cells outside the rabbit host has facilitated its development as an oncolytic virus for the treatment of different types of cancers. Like other poxviruses, MYXV possesses a large dsDNA genome which encodes an array of dozens of immunomodulatory proteins that are important for host and cellular tropism and modulation of host antiviral innate immune responses, some of which are rabbit-specific and others can function in nonrabbit cells as well. This review summarizes the functions of one such MYXV host range protein, M029, an ortholog of the larger superfamily of poxvirus encoded E3-like dsRNA binding proteins. M029 has been identified as a multifunctional protein involved in MYXV cellular and host tropism, antiviral responses, and pathogenicity in rabbits.

Coinfections of Novel Polyomavirus, Anelloviruses and a Recombinant Strain of Myxoma Virus-MYXV-Tol Identified in Iberian Hares

Viruses are ubiquitous in nature; however, very few have been identified in the Leporid species. In the fall of 2018, an outbreak of myxomatosis in Iberian hares (Lepus granatensis) was reported in Spain and a novel recombinant myxoma virus strain (MYXV-Tol) was identified. To investigate variability within the recombinant region of the MYXV-Tol and identify any potential viral coinfections, samples (ear, eyelid or vaginal) of Iberian hares were collected from Spain and analyzed. The presence of the recombinant region of the MYXV-Tol was confirmed in six out of eleven samples analyzed. Additionally, a polyomavirus (family Polyomaviridae), representing a putative new species, and anelloviruses (family Anelloviridae) belonging to two putative species were identified, some as coinfection with the recombinant MYXV-Tol. The two polyomavirus genomes were identified in two hares and share >99% genome-wide identity. Based on the analysis of their large T-antigen, the new polyomavirus clusters in a distant clade from other mammals sharing <64% amino acid identity. A total of 14 anelloviruses were identified, which share 63-99% genome-wide identity. Overall, our results show a coinfection of different DNA viruses in the studied samples and raise awareness regarding the extensive unsampled diversity of viruses in hares.

Species-Specific Host-Virus Interactions: Implications for Viral Host Range and Virulence

A growing number of studies indicate that host species-specific and virus strain-specific interactions of viral molecules with the host innate immune system play a pivotal role in determining virus host range and virulence. Because interacting proteins are likely constrained in their evolution, mutations that are selected to improve virus replication in one species may, by chance, alter the ability of a viral antagonist to inhibit immune responses in hosts the virus has not yet encountered. Based on recent findings of host-species interactions of poxvirus, herpesvirus, and influenza virus proteins, we propose a model for viral fitness and host range which considers the full interactome between a specific host species and a virus, resulting from the combination of all interactions, positive and negative, that influence whether a virus can productively infect a cell and cause disease in different hosts.

Myxoma virus jumps species to the Iberian hare

The study of myxoma virus (MYXV) infections in the European rabbit (Oryctolagus cuniculus) has produced one of the most accepted host-pathogen evolutionary models. To date, myxomatosis has been limited to the European rabbit with sporadic reports in hares. However, reports of widespread mortalities in the Iberian hare (Lepus granatensis) with myxomatosis-like clinical signs indicate a potential species jump has occurred. The presence of MYXV DNA was confirmed by PCR in 244 samples received from regional veterinary services, animal health laboratories, hunters or rangers over a 5-month period. PCR analysis of 4 MYXV positive hare samples revealed a 2.8 kb insertion located within the M009 gene with respect to MYXV. The presence of this insertion was subsequently confirmed in 20 samples from 18 Spanish provinces. Sanger sequencing and subsequent analysis show that the insert contained 4 ORFs which are phylogenetically related to MYXV genes M060, M061, M064 and M065. The complete MYXV genome from hare tissue was sequenced using Ion torrent next-generation technology and a summary of the data presented here. With the exception of the inserted region, the virus genome had no large scale modifications and 110 mutations with respect to the MYXV reference strain Lausanne were observed. The next phase in the evolution of MYXV has taken place as a host species jump from the European rabbit to the Iberian hare an occurrence which could have important effects on this naïve population.

Detection of viral DNA of myxoma virus using a validated PCR method with an internal amplification control

Recognition of myxomatosis is usually based on clinical symptoms, but amyxomatous cases of the disease require the use of laboratory methods. Nowadays PCR assays are routinely employed for detection of MYXV DNA, but none of them have had their diagnostic usefulness conclusively confirmed through validation. The aim of the study was the development and validation of a PCR with an internal amplification control (IAC) for intravital and postmortem detection of viral DNA of myxoma virus. To avoid false negative results a chimeric internal amplification control (IAC) was prepared and incorporated into the PCR and amplified by the same primer set as the target DNA (M071L). The optimal concentration of particular ingredients in the PCR mixture (including IAC concentration and volume of DNA sample) was determined. To minimize the risk of amplicon carry-over contamination, uracil N-glycosylase was added to the reaction. Before proper validation the robustness of the IAC-PCR was verified. Validation of the method encompassed the following parameters: the analytical and diagnostic specificity (ASp, DSp) and sensitivity (ASe, DSe) of the assay, repeatability, and intra-laboratory reproducibility. The assay LOD was established at 2 TCIU of the virus particles/0.2 ml tissue homogenate with a 100% capacity to detect different MYXV strains (ASp). The method was characterized by good DSp of 0.955 (0.839-0.999 CI) and DSe of 0.976 (0.914-1.00 CI). In addition, it was repeatable and reproducible and confirmed its suitability for the detection of MYXV in clinical material. The IAC-PCR developed meets OIE validation requirements for virological methods and can be used in diagnostic or epidemiological studies of rabbit myxomatosis.

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.

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.

Fogarty International Center collaborative networks in infectious disease modeling: Lessons learnt in research and capacity building

Due to a combination of ecological, political, and demographic factors, the emergence of novel pathogens has been increasingly observed in animals and humans in recent decades. Enhancing global capacity to study and interpret infectious disease surveillance data, and to develop data-driven computational models to guide policy, represents one of the most cost-effective, and yet overlooked, ways to prepare for the next pandemic. Epidemiological and behavioral data from recent pandemics and historic scourges have provided rich opportunities for validation of computational models, while new sequencing technologies and the 'big data' revolution present new tools for studying the epidemiology of outbreaks in real time. For the past two decades, the Division of International Epidemiology and Population Studies (DIEPS) of the NIH Fogarty International Center has spearheaded two synergistic programs to better understand and devise control strategies for global infectious disease threats. The Multinational Influenza Seasonal Mortality Study (MISMS) has strengthened global capacity to study the epidemiology and evolutionary dynamics of influenza viruses in 80 countries by organizing international research activities and training workshops. The Research and Policy in Infectious Disease Dynamics (RAPIDD) program and its precursor activities has established a network of global experts in infectious disease modeling operating at the research-policy interface, with collaborators in 78 countries. These activities have provided evidence-based recommendations for disease control, including during large-scale outbreaks of pandemic influenza, Ebola and Zika virus. Together, these programs have coordinated international collaborative networks to advance the study of emerging disease threats and the field of computational epidemic modeling. A global community of researchers and policy-makers have used the tools and trainings developed by these programs to interpret infectious disease patterns in their countries, understand modeling concepts, and inform control policies. Here we reflect on the scientific achievements and lessons learnt from these programs (h-index = 106 for RAPIDD and 79 for MISMS), including the identification of outstanding researchers and fellows; funding flexibility for timely research workshops and working groups (particularly relative to more traditional investigator-based grant programs); emphasis on group activities such as large-scale modeling reviews, model comparisons, forecasting challenges and special journal issues; strong quality control with a light touch on outputs; and prominence of training, data-sharing, and joint publications.

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.

Parallel molecular evolution and adaptation in viruses

Parallel molecular evolution is the independent evolution of the same genotype or phenotype from distinct ancestors. The simple genomes and rapid evolution of many viruses mean they are useful model systems for studying parallel evolution by natural selection. Parallel adaptation occurs in the context of several viral behaviours, including cross-species transmission, drug resistance, and host immune escape, and its existence suggests that at least some aspects of virus evolution and emergence are repeatable and predictable. We introduce examples of virus parallel evolution and summarise key concepts. We outline the difficulties in detecting parallel adaptation using virus genomes, with a particular focus on phylogenetic and structural approaches, and we discuss future approaches that may improve our understanding of the phenomenon.

Evolutionary Virology at 40

RNA viruses are diverse, abundant, and rapidly evolving. Genetic data have been generated from virus populations since the late 1970s and used to understand their evolution, emergence, and spread, culminating in the generation and analysis of many thousands of viral genome sequences. Despite this wealth of data, evolutionary genetics has played a surprisingly small role in our understanding of virus evolution. Instead, studies of RNA virus evolution have been dominated by two very different perspectives, the experimental and the comparative, that have largely been conducted independently and sometimes antagonistically. Here, we review the insights that these two approaches have provided over the last 40 years. We show that experimental approaches using in vitro and in vivo laboratory models are largely focused on short-term intrahost evolutionary mechanisms, and may not always be relevant to natural systems. In contrast, the comparative approach relies on the phylogenetic analysis of natural virus populations, usually considering data collected over multiple cycles of virus-host transmission, but is divorced from the causative evolutionary processes. To truly understand RNA virus evolution it is necessary to meld experimental and comparative approaches within a single evolutionary genetic framework, and to link viral evolution at the intrahost scale with that which occurs over both epidemiological and geological timescales. We suggest that the impetus for this new synthesis may come from methodological advances in next-generation sequencing and metagenomics.

The phylogenomics of evolving virus virulence

How virulence evolves after a virus jumps to a new host species is central to disease emergence. Our current understanding of virulence evolution is based on insights drawn from two perspectives that have developed largely independently: long-standing evolutionary theory based on limited real data examples that often lack a genomic basis, and experimental studies of virulence-determining mutations using cell culture or animal models. A more comprehensive understanding of virulence mutations and their evolution can be achieved by bridging the gap between these two research pathways through the phylogenomic analysis of virus genome sequence data as a guide to experimental study.

The life history strategy of a fur seal hookworm in relation to pathogenicity and host health status

The strategies that parasites use to exploit their hosts can lead to adverse effects on human and animal populations. Here, we describe the life cycle, epidemiology, and consequences of hookworm (Uncinaria sp.) disease in South American fur seals (Arctocephalus australis), and propose that hookworm adaptation to fur seal life history traits has led to maximizing transmission at high levels of parasite-induced anemia and mortality. Fur seal pups acquire hookworms during their first days of life through their mothers' colostrum and most adult hookworms are expelled from the pups' intestine 30-65 days later. This gives hookworms little time to feed and reproduce. However, despite reaching high within-host densities, female hookworms do not decrease egg output, therefore pups with high hookworm burden contribute disproportionately to parasite egg shedding. These heavily infected pups also suffer severe anemia and high levels of hookworm-induced mortality. Alternative strategies to maximize total egg shedding and/or transmission, such as increased environmental survival of larval stages or avoidance of clearance, have not been developed by this hookworm. We propose that fur seal hookworms exploit a live fast-die young life history strategy, which translates to the highest levels of host anemia and mortality recorded among hookworms.

Oncolytic Viruses for Canine Cancer Treatment

Oncolytic virotherapy has been investigated for several decades and is emerging as a plausible biological therapy with several ongoing clinical trials and two viruses are now approved for cancer treatment in humans. The direct cytotoxicity and immune-stimulatory effects make oncolytic viruses an interesting strategy for cancer treatment. In this review, we summarize the results of in vitro and in vivo published studies of oncolytic viruses in different phases of evaluation in dogs, using PubMed and Google scholar as search platforms, without time restrictions (to date). Natural and genetically modified oncolytic viruses were evaluated with some encouraging results. The most studied viruses to date are the reovirus, myxoma virus, and vaccinia, tested mostly in solid tumors such as osteosarcomas, mammary gland tumors, soft tissue sarcomas, and mastocytomas. Although the results are promising, there are issues that need addressing such as ensuring tumor specificity, developing optimal dosing, circumventing preexisting antibodies from previous exposure or the development of antibodies during treatment, and assuring a reasonable safety profile, all of which are required in order to make this approach a successful therapy in dogs.

An Agent-Based Metapopulation Model Simulating Virus-Based Biocontrol of Heterodera Glycines

With recently discovered soybean cyst nematode (SCN) viruses, biological control of the nematodes is a theoretical possibility. This study explores the question of what kinds of viruses would make useful biocontrol agents, taking into account evolutionary and population dynamics. An agent-based model, Soybean Cyst Nematode Simulation (SCNSim), was developed to simulate within-host virulence evolution in a virus-nematode-soybean ecosystem. SCNSim was used to predict nematode suppression under a range of viral mutation rates, initial virulences, and release strategies. The simulation model suggested that virus-based biocontrol worked best when the nematodes were inundated with the viruses. Under lower infection prevalence, the viral burden thinned out rapidly due to the limited mobility and high reproductive rate of the SCN. In accordance with the generally accepted trade-off theory, SCNSim predicted the optimal initial virulence for the maximum nematode suppression. Higher initial virulence resulted in shorter lifetime transmission, whereas viruses with lower initial virulence values evolved toward avirulence. SCNSim also indicated that a greater viral mutation rate reinforced the virulence pathotype, suggesting the presence of a virulence threshold necessary to achieve biocontrol against SCN.

Variable changes in nematode infection prevalence and intensity after Rabbit Haemorrhagic Disease Virus emerged in wild rabbits in Scotland and New Zealand

The myxoma virus (a microparasite) reduced wild rabbit numbers worldwide when introduced in the 1950s, and is known to interact with co-infecting helminths (macroparasites) causing both increases and decreases in macroparasite population size. In the 1990s Rabbit Haemorrhagic Disease Virus (RHDV) infected rabbits and also significantly reduced rabbit numbers in several countries. However, not much is known about RHDV interactions with macroparasites. In this study, we compare prevalence and intensity of infection for three gastrointestinal nematode species (Trichostrongylus retortaeformis, Graphidium strigosum and Passalurus ambiguus) before and after RHDV spread across host populations in Scotland and New Zealand. During one common season, autumn, prevalence of T. retortaeformis was higher after RHDV spread in both locations, whereas it was lower for G. strigosum and P. ambiguus after RHDV arrived in New Zealand, but higher in Scotland. Meanwhile, intensity of infection for all species decreased after RHDV arrived in New Zealand, but increased in Scotland. The impact of RHDV on worm infections was generally similar across seasons in Scotland, and also similarities in seasonality between locations suggested effects on infection patterns in one season are likely similar year-round. The variable response by macroparasites to the arrival of a microparasite into Scottish and New Zealand rabbits may be due to differences in the environment they inhabit, in existing parasite community structure, and to some extent, in the relative magnitude of indirect effects. Specifically, our data suggest that bottom-up processes after the introduction of a more virulent strain of RHDV to New Zealand may affect macroparasite co-infections by reducing the availability of their shared common resource, the rabbits. Clearly, interactions between co-infecting micro- and macroparasites vary in host populations with different ecologies, and significantly impact parasite community structure in wildlife.

•

Nematode communities in Scotland and New Zealand were compared pre and post Rabbit Haemorrhagic Disease Virus introduction.

•

Similar species occur in both rabbit populations, but prevalence and intensity changed in opposing directions after RHDV.

•

RHDV had a major impact on rabbit populations, and our data show differing impacts on macroparasites in the two countries.

•

Variability in rabbit environment, parasite community structure, and indirect interaction processes may explain differences.

•

Results can help understand interactions between co-infecting parasites and their epidemiology in wild and domestic animals.

Bacterial Pathogen Emergence Requires More than Direct Contact with a Novel Passerine Host

While direct contact may sometimes be sufficient to allow a pathogen to jump into a new host species, in other cases, fortuitously adaptive mutations that arise in the original donor host are also necessary. Viruses have been the focus of most host shift studies, so less is known about the importance of ecological versus evolutionary processes to successful bacterial host shifts. Here we tested whether direct contact with the novel host was sufficient to enable the mid-1990s jump of the bacterium Mycoplasma gallisepticum from domestic poultry to house finches (Haemorhous mexicanus). We experimentally inoculated house finches with two genetically distinct M. gallisepticum strains obtained either from poultry (Rlow) or from house finches (HF1995) during an epizootic outbreak. All 15 house finches inoculated with HF1995 became infected, whereas Rlow successfully infected 12 of 15 (80%) inoculated house finches. Comparisons among infected birds showed that, relative to HF1995, Rlow achieved substantially lower bacterial loads in the host respiratory mucosa and was cleared faster. Furthermore, Rlow-infected finches were less likely to develop clinical symptoms than HF1995-infected birds and, when they did, displayed milder conjunctivitis. The lower infection success of Rlow relative to HF1995 was not, however, due to a heightened host antibody response to Rlow. Taken together, our results indicate that contact between infected poultry and house finches was not, by itself, sufficient to explain the jump of M. gallisepticum to house finches. Instead, mutations arising in the original poultry host would have been necessary for successful pathogen emergence in the novel finch host.

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.

Understanding Immunity through the Lens of Disease Ecology

As we describe the immune system in ever more exquisite detail, we might find that no matter how successful, this approach will not be sufficient to understand the spread of infectious agents, their susceptibility to vaccine therapy, and human disease resistance. Compared with the strict reductionism practiced as a means of characterizing most biological processes, I propose that the progression and outcome of disease-causing host-parasite interactions will be more clearly understood through a focus on disease ecology.

A phylogenomic analysis of Marek's disease virus reveals independent paths to virulence in Eurasia and North America

Virulence determines the impact a pathogen has on the fitness of its host, yet current understanding of the evolutionary origins and causes of virulence of many pathogens is surprisingly incomplete. Here, we explore the evolution of Marek's disease virus (MDV), a herpesvirus commonly afflicting chickens and rarely other avian species. The history of MDV in the 20th century represents an important case study in the evolution of virulence. The severity of MDV infection in chickens has been rising steadily since the adoption of intensive farming techniques and vaccination programs in the 1950s and 1970s, respectively. It has remained uncertain, however, which of these factors is causally more responsible for the observed increase in virulence of circulating viruses. We conducted a phylogenomic study to understand the evolution of MDV in the context of dramatic changes to poultry farming and disease control. Our analysis reveals evidence of geographical structuring of MDV strains, with reconstructions supporting the emergence of virulent viruses independently in North America and Eurasia. Of note, the emergence of virulent viruses appears to coincide approximately with the introduction of comprehensive vaccination on both continents. The time‐dated phylogeny also indicated that MDV has a mean evolutionary rate of ~1.6 × 10⁻⁵ substitutions per site per year. An examination of gene‐linked mutations did not identify a strong association between mutational variation and virulence phenotypes, indicating that MDV may evolve readily and rapidly under strong selective pressures and that multiple genotypic pathways may underlie virulence adaptation in MDV.

Next step in the ongoing arms race between myxoma virus and wild rabbits in Australia is a novel disease phenotype

In host-pathogen arms races, increases in host resistance prompt counteradaptation by pathogens, but the nature of that counteradaptation is seldom directly observed outside of laboratory models. The best-documented field example is the coevolution of myxoma virus (MYXV) in European rabbits. To understand how MYXV in Australia has continued to evolve in wild rabbits under intense selection for genetic resistance to myxomatosis, we compared the phenotypes of the progenitor MYXV and viral isolates from the 1950s and the 1990s in laboratory rabbits with no resistance. Strikingly, and unlike their 1950s counterparts, most virus isolates from the 1990s induced a highly lethal immune collapse syndrome similar to septic shock. Thus, the next step in this canonical case of coevolution after a species jump has been further escalation by the virus in the face of widespread host resistance.

Genomic and phenotypic characterization of myxoma virus from Great Britain reveals multiple evolutionary pathways distinct from those in Australia

The co-evolution of myxoma virus (MYXV) and the European rabbit occurred independently in Australia and Europe from different progenitor viruses. Although this is the canonical study of the evolution of virulence, whether the genomic and phenotypic outcomes of MYXV evolution in Europe mirror those observed in Australia is unknown. We addressed this question using viruses isolated in the United Kingdom early in the MYXV epizootic (1954-1955) and between 2008-2013. The later UK viruses fell into three distinct lineages indicative of a long period of separation and independent evolution. Although rates of evolutionary change were almost identical to those previously described for MYXV in Australia and strongly clock-like, genome evolution in the UK and Australia showed little convergence. The phenotypes of eight UK viruses from three lineages were characterized in laboratory rabbits and compared to the progenitor (release) Lausanne strain. Inferred virulence ranged from highly virulent (grade 1) to highly attenuated (grade 5). Two broad disease types were seen: cutaneous nodular myxomatosis characterized by multiple raised secondary cutaneous lesions, or an amyxomatous phenotype with few or no secondary lesions. A novel clinical outcome was acute death with pulmonary oedema and haemorrhage, often associated with bacteria in many tissues but an absence of inflammatory cells. Notably, reading frame disruptions in genes defined as essential for virulence in the progenitor Lausanne strain were compatible with the acquisition of high virulence. Combined, these data support a model of ongoing host-pathogen co-evolution in which multiple genetic pathways can produce successful outcomes in the field that involve both different virulence grades and disease phenotypes, with alterations in tissue tropism and disease mechanisms.

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.

Myxoma Virus dsRNA Binding Protein M029 Inhibits the Type I IFN-Induced Antiviral State in a Highly Species-Specific Fashion

Myxoma virus (MYXV) is a Leporipoxvirus that possesses a specific rabbit-restricted host tropism but exhibits a much broader cellular host range in cultured cells. MYXV is able to efficiently block all aspects of the type I interferon (IFN)-induced antiviral state in rabbit cells, partially in human cells and very poorly in mouse cells. The mechanism(s) of this species-specific inhibition of type I IFN-induced antiviral state is not well understood. Here we demonstrate that MYXV encoded protein M029, a truncated relative of the vaccinia virus (VACV) E3 double-stranded RNA (dsRNA) binding protein that inhibits protein kinase R (PKR), can also antagonize the type I IFN-induced antiviral state in a highly species-specific manner. In cells pre-treated with type I IFN prior to infection, MYXV exploits M029 to overcome the induced antiviral state completely in rabbit cells, partially in human cells, but not at all in mouse cells. However, in cells pre-infected with MYXV, IFN-induced signaling is fully inhibited even in the absence of M029 in cells from all three species, suggesting that other MYXV protein(s) apart from M029 block IFN signaling in a species-independent manner. We also show that the antiviral state induced in rabbit, human or mouse cells by type I IFN can inhibit M029-knockout MYXV even when PKR is genetically knocked-out, suggesting that M029 targets other host proteins for this antiviral state inhibition. Thus, the MYXV dsRNA binding protein M029 not only antagonizes PKR from multiple species but also blocks the type I IFN antiviral state independently of PKR in a highly species-specific fashion.

Emergence of a Viral RNA Polymerase Variant during Gene Copy Number Amplification Promotes Rapid Evolution of Vaccinia Virus

Viruses are under relentless selective pressure from host immune defenses. To study how poxviruses adapt to innate immune detection pathways, we performed serial vaccinia virus infections in primary human cells. Independent courses of experimental evolution with a recombinant strain lacking E3L revealed several high-frequency point mutations in conserved poxvirus genes, suggesting important roles for essential poxvirus proteins in innate immune subversion. Two distinct mutations were identified in the viral RNA polymerase gene A24R, which seem to act through different mechanisms to increase virus replication. Specifically, a Leu18Phe substitution encoded within A24R conferred fitness trade-offs, including increased activation of the antiviral factor protein kinase R (PKR). Intriguingly, this A24R variant underwent a drastic selective sweep during passaging, despite enhanced PKR activity. We showed that the sweep of this variant could be accelerated by the presence of copy number variation (CNV) at the K3L locus, which in multiple copies strongly reduced PKR activation. Therefore, adaptive cases of CNV can facilitate the accumulation of point mutations separate from the expanded locus. This study reveals how rapid bouts of gene copy number amplification during accrual of distant point mutations can potently facilitate poxvirus adaptation to host defenses.

IMPORTANCE

Viruses can evolve quickly to defeat host immune functions. For poxviruses, little is known about how multiple adaptive mutations emerge in populations at the same time. In this study, we uncovered a means of vaccinia virus adaptation involving the accumulation of distinct genetic variants within a single population. We identified adaptive point mutations in the viral RNA polymerase gene A24R and, surprisingly, found that one of these mutations activates the nucleic acid sensing factor PKR. We also found that gene copy number variation (CNV) can provide dual benefits to evolving virus populations, including evidence that CNV facilitates the accumulation of a point mutation distant from the expanded locus. Our data suggest that transient CNV can accelerate the fixation of mutations conferring modest benefits, or even fitness trade-offs, and highlight how structural variation might aid poxvirus adaptation through both direct and indirect actions.

17th Century Variola Virus Reveals the Recent History of Smallpox

Smallpox holds a unique position in the history of medicine. It was the first disease for which a vaccine was developed and remains the only human disease eradicated by vaccination. Although there have been claims of smallpox in Egypt, India, and China dating back millennia [1-4], the timescale of emergence of the causative agent, variola virus (VARV), and how it evolved in the context of increasingly widespread immunization, have proven controversial [4-9]. In particular, some molecular-clock-based studies have suggested that key events in VARV evolution only occurred during the last two centuries [4-6] and hence in apparent conflict with anecdotal historical reports, although it is difficult to distinguish smallpox from other pustular rashes by description alone. To address these issues, we captured, sequenced, and reconstructed a draft genome of an ancient strain of VARV, sampled from a Lithuanian child mummy dating between 1643 and 1665 and close to the time of several documented European epidemics [1, 2, 10]. When compared to vaccinia virus, this archival strain contained the same pattern of gene degradation as 20^th century VARVs, indicating that such loss of gene function had occurred before ca. 1650. Strikingly, the mummy sequence fell basal to all currently sequenced strains of VARV on phylogenetic trees. Molecular-clock analyses revealed a strong clock-like structure and that the timescale of smallpox evolution is more recent than often supposed, with the diversification of major viral lineages only occurring within the 18^th and 19^th centuries, concomitant with the development of modern vaccination.

Benign Rabbit Caliciviruses Exhibit Evolutionary Dynamics Similar to Those of Their Virulent Relatives

Two closely related caliciviruses cocirculate in Australia: rabbit hemorrhagic disease virus (RHDV) and rabbit calicivirus Australia 1 (RCV-A1). RCV-A1 causes benign enteric infections in the European rabbit (Oryctolagus cuniculus) in Australia and New Zealand, while its close relative RHDV causes a highly pathogenic infection of the liver in the same host. The comparison of these viruses provides important information on the nature and trajectory of virulence evolution, particularly as highly virulent strains of RHDV may have evolved from nonpathogenic ancestors such as RCV-A1. To determine the evolution of RCV-A1 we sequenced the full-length genomes of 44 RCV-A1 samples isolated from healthy rabbits and compared key evolutionary parameters to those of its virulent relative, RHDV. Despite their marked differences in pathogenicity and tissue tropism, RCV-A1 and RHDV have evolved in a very similar manner. Both viruses have evolved at broadly similar rates, suggesting that their dynamics are largely shaped by high background mutation rates, and both exhibit occasional recombination and an evolutionary environment dominated by purifying selection. In addition, our comparative analysis revealed that there have been multiple changes in both virulence and tissue tropism in the evolutionary history of these and related viruses. Finally, these new genomic data suggest that either RCV-A1 was introduced into Australia after the introduction of myxoma virus as a biocontrol agent in 1950 or there was drastic reduction of the rabbit population, and hence of RCV-A1 genetic diversity, perhaps coincident with the emergence of myxoma virus.

IMPORTANCE

The comparison of closely related viruses that differ profoundly in propensity to cause disease in their hosts offers a powerful opportunity to reveal the causes of changes in virulence and to study how such changes alter the evolutionary dynamics of these pathogens. Here we describe such a novel comparison involving two closely related RNA viruses that cocirculate in Australia, the highly virulent rabbit hemorrhagic disease virus (RHDV) and the nonpathogenic rabbit calicivirus Australia 1 (RCV-A1). Both viruses infect the European rabbit, but they differ in virulence, tissue tropism, and mechanisms of transmission. Surprisingly, and despite these fundamental differences, RCV-A1 and RHDV have evolved at very similar (high) rates and with strong purifying selection. Furthermore, candidate key mutations were identified that may play a role in virulence and/or tissue tropism and therefore warrant further investigation.

Host and Parasite Evolution in a Tangled Bank

Most hosts and parasites exist in diverse communities wherein they interact with other species, spanning the parasite-mutualist continuum. These additional interactions have the potential to impose selection on hosts and parasites and influence the patterns and processes of their evolution. Yet, host-parasite interactions are almost exclusively studied in species pairs. A wave of new research has incorporated a multispecies community context, showing that additional ecological interactions can alter components of host and parasite fitness, as well as interaction specificity and virulence. Here, we synthesize these findings to assess the effects of increased species diversity on the patterns and processes of host and parasite evolution. We argue that our understanding of host-parasite interactions would benefit from a richer biotic perspective.

Link of a ubiquitous human coronavirus to dromedary camels

The four human coronaviruses (HCoVs) are globally endemic respiratory pathogens. The Middle East respiratory syndrome (MERS) coronavirus (CoV) is an emerging CoV with a known zoonotic source in dromedary camels. Little is known about the origins of endemic HCoVs. Studying these viruses' evolutionary history could provide important insight into CoV emergence. In tests of MERS-CoV-infected dromedaries, we found viruses related to an HCoV, known as HCoV-229E, in 5.6% of 1,033 animals. Human- and dromedary-derived viruses are each monophyletic, suggesting ecological isolation. One gene of dromedary viruses exists in two versions in camels, full length and deleted, whereas only the deleted version exists in humans. The deletion increased in size over a succession starting from camelid viruses via old human viruses to contemporary human viruses. Live isolates of dromedary 229E viruses were obtained and studied to assess human infection risks. The viruses used the human entry receptor aminopeptidase N and replicated in human hepatoma cells, suggesting a principal ability to cause human infections. However, inefficient replication in several mucosa-derived cell lines and airway epithelial cultures suggested lack of adaptation to the human host. Dromedary viruses were as sensitive to the human type I interferon response as HCoV-229E. Antibodies in human sera neutralized dromedary-derived viruses, suggesting population immunity against dromedary viruses. Although no current epidemic risk seems to emanate from these viruses, evolutionary inference suggests that the endemic human virus HCoV-229E may constitute a descendant of camelid-associated viruses. HCoV-229E evolution provides a scenario for MERS-CoV emergence.