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

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.

First Detection of Benign Rabbit Caliciviruses in Chile

Pathogenic lagoviruses (Rabbit hemorrhagic disease virus, RHDV) are widely spread across the world and are used in Australia and New Zealand to control populations of feral European rabbits. The spread of the non-pathogenic lagoviruses, e.g., rabbit calicivirus (RCV), is less well studied as the infection results in no clinical signs. Nonetheless, RCV has important implications for the spread of RHDV and rabbit biocontrol as it can provide varying levels of cross-protection against fatal infection with pathogenic lagoviruses. In Chile, where European rabbits are also an introduced species, myxoma virus was used for localised biocontrol of rabbits in the 1950s. To date, there have been no studies investigating the presence of lagoviruses in the Chilean feral rabbit population. In this study, liver and duodenum rabbit samples from central Chile were tested for the presence of lagoviruses and positive samples were subject to whole RNA sequencing and subsequent data analysis. Phylogenetic analysis revealed a novel RCV variant in duodenal samples that likely originated from European RCVs. Sequencing analysis also detected the presence of a rabbit astrovirus in one of the lagovirus-positive samples.

Effect of epidemic diseases on wild animal conservation

Under the background of global species extinction, the impact of epidemic diseases on wild animal protection is increasingly prominent. Here, we review and synthesize the literature on this topic, and discuss the relationship between diseases and biodiversity. Diseases usually reduce species diversity by decreasing or extinction of species populations, but also accelerate species evolution and promote species diversity. At the same time, species diversity can regulate disease outbreaks through dilution or amplification effects. The synergistic effect of human activities and global change is emphasized, which further aggravates the complex relationship between biodiversity and diseases. Finally, we emphasize the importance of active surveillance of wild animal diseases, which can protect wild animals from potential diseases, maintain population size and genetic variation, and reduce the damage of diseases to the balance of the whole ecosystem and human health. Therefore, we suggest that a background survey of wild animal populations and their pathogens should be carried out to assess the impact of potential outbreaks on the population or species level. The mechanism of dilution and amplification effect between species diversity and diseases of wild animals should be further studied to provide a theoretical basis and technical support for human intervention measures to change biodiversity. Most importantly, we should closely combine the protection of wild animals with the establishment of an active surveillance, prevention, and control system for wild animal epidemics, in an effort to achieve a win-win situation between wild animal protection and disease control.

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.

Evaluation of Commercial Myxomatosis Vaccines against Recombinant Myxoma Virus (ha-MYXV) in Iberian Hare and Wild Rabbit

The recent emergence of a new myxoma virus capable of causing disease in the Iberian hare (Lepus granatensis) has resulted in numerous outbreaks with high mortality leading to the reduction, or even the disappearance, of many local populations of this wild species in the Iberian Peninsula. Currently, the available vaccines that prevent myxomatosis in domestic rabbits caused by classic strains of myxoma virus have not been assessed for use in Iberian hares. The main objective of this study was to evaluate the efficacy of commercial rabbit vaccines in Iberian hares and wild rabbits against the natural recombinant myxoma virus (ha-MYXV), bearing in mind its application in specific scenarios where capture is possible, such as genetic reserves. The study used a limited number of animals (pilot study), 15 Iberian hares and 10 wild rabbits. Hares were vaccinated with Mixohipra-FSA vaccine (Hipra) and Mixohipra-H vaccine (Hipra) using two different doses, and rabbits were vaccinated with the Mixohipra-H vaccine or the Nobivac Myxo-RHD PLUS (MSD Animal Health) using the recommended doses for domestic rabbits. After the vaccination trials, the animals were challenged with a wild type strain of ha-MYXV. The results showed that no protection to ha-MYXV challenge was afforded when a commercial dose of Mixohipra-FSA or Mixohipra-H vaccine was used in hares. However, the application of a higher dose of Mixohipra-FSA vaccine may induce protection and could possibly be used to counteract the accelerated decrease of wild hare populations due to ha-MYXV emergence. The two commercial vaccines (Mixohipra-H and Nobivac Myxo-RHD PLUS) tested in wild rabbits were fully protective against ha-MYXV infection. This knowledge gives more insights into ha-MYXV management in hares and rabbits and emphasises the importance of developing a vaccine capable of protecting wild populations of Iberian hare and wild rabbit towards MYXV and ha-MYXV strains.

HUMAN ADENOVIRUS TYPE 4 COMPRISES TWO MAJOR PHYLOGROUPS WITH DISTINCT REPLICATIVE FITNESS AND VIRULENCE PHENOTYPES

Human adenovirus type 4 (HAdV-E4) is the only type (and serotype) classified within species Human mastadenovirus E that has been isolated from a human host to the present. Recent phylogenetic analysis of whole genome sequences of strains representing the spectrum of intratypic genetic diversity described to date identified two major evolutionary lineages designated phylogroups (PG) I, and II, and validated the early clustering of HAdV-E4 genomic variants into two major groups by low resolution restriction fragment length polymorphism analysis. In this study we expanded our original analysis of intra- and inter-PG genetic variability, and used a panel of viruses representative of the spectrum of genetic diversity described for HAdV-E4 to examine the magnitude of inter- and intra-PG phenotypic diversity using an array of cell-based assays and a cotton rat model of HAdV respiratory infection. Our proteotyping of HAdV-E strains using concatenated protein sequences in selected coding regions including E1A, E1B-19K and -55K, DNA polymerase, L4-100K, various E3 proteins, and E4-34K confirmed that the two clades encode distinct variants/proteotypes at most of these loci. Our in vitro and in vivo studies demonstrated that PG I and PG II differ in their growth, spread, and cell killing phenotypes in cell culture and in their pulmonary pathogenic phenotypes. Surprisingly, the differences in replicative fitness documented in vitro between PGs did not correlate with the differences in virulence observed in the cotton rat model. This body of work is the first reporting phenotypic correlates of naturally occurring intratypic genetic variability for HAdV-E4. IMPORTANCE Human adenovirus type 4 (HAdV-E4) is a prevalent causative agent of acute respiratory illness of variable severity and of conjunctivitis and comprises two major phylogroups that carry distinct coding variations in proteins involved in viral replication and modulation of host responses to infection. Our data show that PG I and PG II are intrinsically different regarding their ability to grow and spread in culture and to cause pulmonary disease in cotton rats. This is the first report of phenotypic divergence among naturally occurring known genetic variants of a HAdV type of medical importance. This research reveals readily detectable phenotypic differences between strains representing phylogroups I and II, and it introduces a unique experimental system for the elucidation of the genetic basis of adenovirus fitness and virulence and thus for increasing our understanding of the implications of intratypic genetic diversity in the presentation and course of HAdV-E4-associated disease.

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.

Monkeypox virus emergence in wild chimpanzees reveals distinct clinical outcomes and viral diversity

Monkeypox is a viral zoonotic disease on the rise across endemic habitats. Despite the growing importance of monkeypox virus, our knowledge on its host spectrum and sylvatic maintenance is limited. Here, we describe the recent repeated emergence of monkeypox virus in a wild, human-habituated western chimpanzee (Pan troglodytes verus, hereafter chimpanzee) population from Taï National Park, Ivory Coast. Through daily monitoring, we show that further to causing its typical exanthematous syndrome, monkeypox can present itself as a severe respiratory disease without a diffuse rash. By analysing 949 non-invasively collected samples, we identify the circulation of at least two distinct monkeypox virus lineages and document the shedding of infectious particles in faeces and flies, suggesting that they could mediate indirect transmission. We also show that the carnivorous component of the Taï chimpanzees' diet, mainly consisting of the sympatric monkeys they regularly hunt, did not change nor shift towards rodent consumption (the presumed reservoir) before the outbreaks, suggesting that the sudden emergence of monkeypox virus in this population is probably due to changes in the ecology of the virus itself. Using long-term mortality surveillance data from Taï National Park, we provide evidence of little to no prior viral activity over at least two decades. We conclude that great ape sentinel systems devoted to the longitudinal collection of behavioural and health data can help clarify the epidemiology and clinical presentation of zoonotic pathogens.

Adaptive changes in the genomes of wild rabbits after 16 years of viral epidemics

Since its introduction to control overabundant invasive European rabbits (Oryctolagus cuniculus), the highly virulent rabbit haemorrhagic disease virus (RHDV) has caused regular annual disease outbreaks in Australian rabbit populations. Although initially reducing rabbit abundance by 60%, continent-wide, experimental evidence has since indicated increased genetic resistance in wild rabbits that have experienced RHDV-driven selection. To identify genetic adaptations, which explain the increased resistance to this biocontrol virus, we investigated genome-wide SNP (single nucleotide polymorphism) allele frequency changes in a South Australian rabbit population that was sampled in 1996 (pre-RHD genomes) and after 16 years of RHDV outbreaks. We identified several SNPs with changed allele frequencies within or close to genes potentially important for increased RHD resistance. The identified genes are known to be involved in virus infections and immune reactions or had previously been identified as being differentially expressed in healthy versus acutely RHDV-infected rabbits. Furthermore, we show in a simulation study that the allele/genotype frequency changes cannot be explained by drift alone and that several candidate genes had also been identified as being associated with surviving RHD in a different Australian rabbit population. Our unique data set allowed us to identify candidate genes for RHDV resistance that have evolved under natural conditions, and over a time span that would not have been feasible in an experimental setting. Moreover, it provides a rare example of host genetic adaptations to virus-driven selection in response to a suddenly emerging infectious disease.

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.

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.

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.

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.

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.