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

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

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.

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 and the Leporipoxviruses: an evolutionary paradigm

Myxoma virus (MYXV) is the type species of the Leporipoxviruses, a genus of Chordopoxvirinae, double stranded DNA viruses, whose members infect leporids and squirrels, inducing cutaneous fibromas from which virus is mechanically transmitted by biting arthropods. However, in the European rabbit (Oryctolagus cuniculus), MYXV causes the lethal disease myxomatosis. The release of MYXV as a biological control for the wild European rabbit population in Australia, initiated one of the great experiments in evolution. The subsequent coevolution of MYXV and rabbits is a classic example of natural selection acting on virulence as a pathogen adapts to a novel host species. Slightly attenuated mutants of the progenitor virus were more readily transmitted by the mosquito vector because the infected rabbit survived longer, while highly attenuated viruses could be controlled by the rabbit immune response. As a consequence, moderately attenuated viruses came to dominate. This evolution of the virus was accompanied by selection for genetic resistance in the wild rabbit population, which may have created an ongoing co-evolutionary dynamic between resistance and virulence for efficient transmission. This natural experiment was repeated on a continental scale with the release of a separate strain of MYXV in France and its subsequent spread throughout Europe. The selection of attenuated strains of virus and resistant rabbits mirrored the experience in Australia in a very different environment, albeit with somewhat different rates. Genome sequencing of the progenitor virus and the early radiation, as well as those from the 1990s in Australia and Europe, has shown that although MYXV evolved at high rates there was no conserved route to attenuation or back to virulence. In contrast, it seems that these relatively large viral genomes have the flexibility for multiple pathways that converge on a similar phenotype.

Myxomatosis in Australia and Europe: a model for emerging infectious diseases

Myxoma virus is a poxvirus naturally found in two American leporid (rabbit) species (Sylvilagus brasiliensis and Sylvilagus bachmani) in which it causes an innocuous localised cutaneous fibroma. However, in European rabbits (Oryctolagus cuniculus) the same virus causes the lethal disseminated disease myxomatosis. The introduction of myxoma virus into the European rabbit population in Australia in 1950 initiated the best known example of what happens when a novel pathogen jumps into a completely naïve new mammalian host species. The short generation time of the rabbit and their vast numbers in Australia meant evolution could be studied in real time. The carefully documented emergence of attenuated strains of virus that were more effectively transmitted by the mosquito vector and the subsequent selection of rabbits with genetic resistance to myxomatosis is the paradigm for pathogen virulence and host-pathogen coevolution. This natural experiment was repeated with the release of a separate strain of myxoma virus in France in 1952. The subsequent spread of the virus throughout Europe and its coevolution with the rabbit essentially paralleled what occurred in Australia. Detailed molecular studies on myxoma virus have dissected the role of virulence genes in the pathogenesis of myxomatosis and when combined with genomic data and reverse genetics should in future enable the understanding of the molecular evolution of the virus as it adapted to its new host. This review describes the natural history and evolution of myxoma virus together with the molecular biology and experimental pathogenesis studies that are informing our understanding of evolution of emerging diseases.

Sequence mapping of the Californian MSW strain of Myxoma virus

Partial sequence mapping of the MSW Californian strain of Myxoma virus was performed by cloning EcoRI and SalI restriction fragments of viral DNA and sequencing the ends of these. In this way, regions of 74 MSW open reading frames were sequenced and mapped onto the complete genome sequences of the related leporipoxviruses South American Myxoma virus and Rabbit fibroma virus to form a partial map of the MSW strain. In general, gene locations and sequences were conserved between the three viruses. However the Californian Myxoma virus was more closely related to South American myxoma virus than to Rabbit fibroma virus based on sequence comparisons and the presence of three genes that have been lost from the Rabbit fibroma virus genome. Compared to the other two viruses, the main difference found in the MSW genome was that the terminal inverted repeats were extended with the duplication of 5 complete open reading frames (M151R, M152R, M153R, M154L, M156R) and partial duplication of one open reading frame (M150R). This rearrangement was associated with the loss of the majority of the M009L open reading frame. Three known virulence genes, including the serine proteinase inhibitor (SERPIN) genes M151R and M152R and leukemia associated protein (LAP) gene M153R, and the potential virulence gene M156R are now present in two copies.

Monitoring the spread of myxoma virus in rabbit Oryctolagus cuniculus populations on the southern tablelands of New South Wales, Australia. III. Release, persistence and rate of spread of an identifiable strain of myxoma virus

An identifiable strain of myxoma virus was introduced into four local populations of wild rabbits Oryctolagus cuniculus on the southern tablelands of New South Wales (NSW) and its spread in the presence of other field strains was monitored for 6 months. The main vector in this region was considered to be the European rabbit flea Spilopsyllis cuniculi. Each population of rabbits was of a high density and living in groups of warrens covering areas from 59 to 87 hectares. Rabbits occupying centrally located warrens were inoculated with the virus in late September or early October (spring) and the subsequent appearance of myxomatosis across the sites monitored by trapping, shooting and visual observations. Samples, taken from rabbits with myxomatosis, were examined by polymerase chain reaction (PCR) that allowed identification of the introduced strain. On all four sites the introduced virus spread from the inoculated rabbits in the centrally located warrens to rabbits in surrounding warrens. On Sites 1 and 3, this spread continued across the entire site persisting for at least 118 and 174 days respectively. On Sites 2 and 4, the virus was detected for 78 and 62 days respectively and the subsequent inability to detect the introduced virus correlated with the appearance of an unrelated field strain. Using three different methods of calculation, rates of spread ranged from 3.7 to 17.8 m d(-1).

Studies on the poxvirus Cotia

The poxvirus Cotia was studied by electron microscopy and by serological and biochemical analyses. Thin-sectioned preparations of infected Vero cells indicated that Cotia virus morphogenesis was similar to other mammalian poxviruses; unique filamentous structures and inclusion matrices were apparent in the cytoplasm. Complement fixation tests that included purified Cotia virions showed a reciprocal cross-reaction with rabbit myxoma virus and no cross-reaction with vaccinia virus. Serological results coupled with gradient polyacrylamide gel electropherograms of the structural proteins of purified Cotia, vaccinia, myxoma and fibroma viruses suggested that Cotia virus was similar to the latter two viruses. Agarose gel electropherograms of cleavage fragments of each of these virus DNAs digested with three separate restriction endonucleases showed that each of these viruses had a unique DNA gel profile.