Perturbative diffraction methods resolve a conformational switch that facilitates a two-step enzymatic mechanism.

Enzymes catalyze biochemical reactions through precise positioning of substrates, cofactors, and amino acids to modulate the transition-state free energy. However, the role of conformational dynamics remains poorly understood due to poor experimental access. This shortcoming is evident with Escherichia coli dihydrofolate reductase (DHFR), a model system for the role of protein dynamics in catalysis, for which it is unknown how the enzyme regulates the different active site environments required to facilitate proton and hydride transfer. Here, we describe ligand-, temperature-, and electric-field-based perturbations during X-ray diffraction experiments to map the conformational dynamics of the Michaelis complex of DHFR. We resolve coupled global and local motions and find that these motions are engaged by the protonated substrate to promote efficient catalysis. This result suggests a fundamental design principle for multistep enzymes in which pre-existing dynamics enable intermediates to drive rapid electrostatic reorganization to facilitate subsequent chemical steps.

Correcting systematic errors in diffraction data with modern scaling algorithms.

X-ray diffraction enables the routine determination of the atomic structure of materials. Key to its success are data-processing algorithms that allow experimenters to determine the electron density of a sample from its diffraction pattern. Scaling, the estimation and correction of systematic errors in diffraction intensities, is an essential step in this process. These errors arise from sample heterogeneity, radiation damage, instrument limitations and other aspects of the experiment. New X-ray sources and sample-delivery methods, along with new experiments focused on changes in structure as a function of perturbations, have led to new demands on scaling algorithms. Classically, scaling algorithms use least-squares optimization to fit a model of common error sources to the observed diffraction intensities to force these intensities onto the same empirical scale. Recently, an alternative approach has been demonstrated which uses a Bayesian optimization method, variational inference, to simultaneously infer merged data along with corrections, or scale factors, for the systematic errors. Owing to its flexibility, this approach proves to be advantageous in certain scenarios. This perspective briefly reviews the history of scaling algorithms and contrasts them with variational inference. Finally, appropriate use cases are identified for the first such algorithm, Careless, guidance is offered on its use and some speculations are made about future variational scaling methods.

Resolving conformational changes that mediate a two-step catalytic mechanism in a model enzyme.

Enzymes catalyze biochemical reactions through precise positioning of substrates, cofactors, and amino acids to modulate the transition-state free energy. However, the role of conformational dynamics remains poorly understood due to lack of experimental access. This shortcoming is evident with E. coli dihydrofolate reductase (DHFR), a model system for the role of protein dynamics in catalysis, for which it is unknown how the enzyme regulates the different active site environments required to facilitate proton and hydride transfer. Here, we present ligand-, temperature-, and electric-field-based perturbations during X-ray diffraction experiments that enable identification of coupled conformational changes in DHFR. We identify a global hinge motion and local networks of structural rearrangements that are engaged by substrate protonation to regulate solvent access and promote efficient catalysis. The resulting mechanism shows that DHFR's two-step catalytic mechanism is guided by a dynamic free energy landscape responsive to the state of the substrate.

A unifying Bayesian framework for merging X-ray diffraction data.

Novel X-ray methods are transforming the study of the functional dynamics of biomolecules. Key to this revolution is detection of often subtle conformational changes from diffraction data. Diffraction data contain patterns of bright spots known as reflections. To compute the electron density of a molecule, the intensity of each reflection must be estimated, and redundant observations reduced to consensus intensities. Systematic effects, however, lead to the measurement of equivalent reflections on different scales, corrupting observation of changes in electron density. Here, we present a modern Bayesian solution to this problem, which uses deep learning and variational inference to simultaneously rescale and merge reflection observations. We successfully apply this method to monochromatic and polychromatic single-crystal diffraction data, as well as serial femtosecond crystallography data. We find that this approach is applicable to the analysis of many types of diffraction experiments, while accurately and sensitively detecting subtle dynamics and anomalous scattering.

Structural Basis for Rabbit Hemorrhagic Disease Virus Antibody Specificity.

Rabbit hemorrhagic disease virus (RHDV) typically causes a fatal disease in rabbits. In Australia, RHDV was imported to control the feral rabbit population, while it poses a severe threat to native rabbits in other countries. RHDV variants are genetically diverse and serological studies using antibodies isolated from infected rabbits or raised against RHDV virus-like particles (VLPs) have found RHDV variants antigenically distinct. In this study, we determined the X-ray crystal structure of an RHDV GI.2 (N11 strain) protruding (P) domain in complex with a diagnostic monoclonal antibody (2D9) Fab. We showed that 2D9 interacted with conserved and variable residues on top of the P domain with nanomolar affinity. To better illustrate 2D9 specificity, we determined the X-ray crystal structure of an RHDV GI.1b (Ast89 strain) that was a 2D9 non-binder. Structural analysis indicated that amino acid substitutions on the GI.1b P domain likely restricted 2D9 binding. Interestingly, a model of the GI.2 P domain-Fab complex superimposed onto a cryo-EM structure of an RHDV VLP revealed that 2D9 Fab molecules clashed with neighboring Fabs and indicated that there was a reduced antibody binding occupancy. Moreover, the RHDV GI.2 histo-blood group antigen (HBGA) co-factor binding site appeared obstructed when 2D9 was modeled on the VLP and suggested that 2D9 might also function by blocking HBGA attachment. Overall, this new data provides the first structural basis of RHDV antibody specificity and explains how amino acid variation at the binding site likely restricts 2D9 cross-reactivity. IMPORTANCE Isolated RHDV antibodies have been used for decades to distinguish between antigenic variants, monitor temporal capsid evolution, and examine neutralizing capacities. In this study, we provided the structural basis for an RHDV GI.2 specific diagnostic antibody (2D9) binding and reveal that a small number of amino acid substitutions at the binding site could differentiate between RHDV GI.2 and GI.1b. This novel structural information provides a framework for understanding how RHDV displays a specific antigenic epitope and engages an antibody at the atomic level. Importantly, part of the 2D9 binding region was earlier reported to contain a neutralizing epitope and our structural modeling as well as recent human norovirus antibody-mediated neutralization studies, suggest that the 2D9 antibody has the potential to block HBGA attachment. These new findings should aid in characterizing antigenic variants and advance the development of novel monoclonal antibodies for diagnostics and therapeutics.

Native SAD phasing at room temperature.

Single-wavelength anomalous diffraction (SAD) is a routine method for overcoming the phase problem when solving macromolecular structures. This technique requires the accurate measurement of intensities to determine differences between Bijvoet pairs. Although SAD experiments are commonly conducted at cryogenic temperatures to mitigate the effects of radiation damage, such temperatures can alter the conformational ensemble of the protein and may impede the merging of data from multiple crystals due to non-uniform freezing. Here, a strategy is presented to obtain high-quality data from room-temperature, single-crystal experiments. To illustrate the strengths of this approach, native SAD phasing at 6.55 keV was used to solve four structures of three model systems at 295 K. The resulting data sets allow automatic phasing and model building, and reveal alternate conformations that reflect the structure of proteins at room temperature.

Virucidal Properties of Photocatalytic Coating on Glass against a Model Human Coronavirus.

The antimicrobial properties of photocatalysts have long been studied. However, most of the available literature describes their antibacterial properties, while knowledge of their antiviral activity is rather scarce. Since the outset of the coronavirus disease 2019 (COVID-19) pandemic, an increasing body of research has suggested their antiviral potential and highlighted the need for further research in this area. In this study, we investigated the virucidal properties of a commercial TiO2-coated photocatalytic glass against a model human coronavirus. Our findings demonstrate that the TiO2-coated glass consistently inactivates coronaviruses upon contact under daylight illumination, in a time-dependent manner. A 99% drop in virus titer was achieved after 3.9 h. The electron micrographs of virus-covered TiO2-glass showed a reduced number of virions compared to control glass. Morphological alterations of TiO2-exposed viruses included deformation, disruption of the viral envelope, and virion ghosts, endorsing the application of this material in the construction of protective elements to mitigate the transmission of viruses. To the best of our knowledge, this is the first report showing direct visual evidence of human coronaviruses being damaged and morphologically altered following exposure to this photocatalyst. IMPORTANCE Surface contamination is an important contributor to SARS-CoV-2 spread. The use of personal protective elements and physical barriers (i.e., masks, gloves, and indoor glass separators) increases safety and has proven invaluable in preventing contagion. Redesigning these barriers so that the virus cannot remain infectious on them could make a difference in COVID-19 epidemiology. The introduction of additives with virucidal activity could potentiate the protective effects of these barriers to serve not only as physical containment but also as virus killers, reducing surface contamination after hand touch or aerosol deposition. We performed in-depth analysis of the kinetics of photocatalysis-triggered coronavirus inactivation on building glass coated with TiO2. This is the first report showing direct visual evidence (electron microscopy) of coronaviruses being morphologically damaged following exposure to this photocatalyst, demonstrating the high potential of this material to be incorporated into daily-life high-touch surfaces, giving them an added value in decelerating the virus spread.

Co-infection by classic MYXV and ha-MYXV in Iberian hare (Lepus granatensis) and European wild rabbit (Oryctolagus cuniculus algirus).

Myxomatosis is an emergent disease in the Iberian hare (Lepus granatensis). In this species, the disease is caused by a natural recombinant virus (ha-myxoma virus [MYXV]) identified for the first time in 2018 and has since been responsible for a large number of outbreaks in Spain and Portugal. The ha-MYXV, which harbours a 2.8 Kb insert-disrupting gene M009L, can also infect and cause disease in wild and domestic rabbits, despite being less frequently identified in rabbits. During the laboratory investigations of wild leporids found dead in Portugal carried out within the scope of a Nacional Surveillance Plan (Dispatch 4757/17, MAFDR), co-infection events by classic (MYXV) and naturally recombinant (ha-MYXV) strains were detected in both one Iberian hare and one European wild rabbit (Oryctolagus cuniculus algirus). These two cases were initially detected by a multiplex qPCR detection of MYXV and ha-MYXV and subsequently confirmed by conventional PCR and sequencing of the M009L gene, which contains an ha-MYXV-specific insertion. To our knowledge, this is the first documented report of co-infection by classic MYXV and ha-MYXV strains either in Iberian hare or in European wild rabbit. It is also the first report of infection of an Iberian hare by a classic MYXV strain. These findings highlight the continuous evolution of the MYXV and the frequent host range changes that justify the nonstop monitoring of the sanitary condition of wild Leporidae populations in the Iberian Peninsula.

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.

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.

Spillover events of rabbit haemorrhagic disease virus 2 (recombinant GI.4P-GI.2) from Lagomorpha to Eurasian badger.

Rabbit haemorrhagic disease (RHD) is a major threat to domestic and wild European rabbits. Presently, in Europe, the disease is caused mainly by Rabbit haemorrhagic disease virus 2 (RHDV2/b or Lagovirus europaeus GI.2), the origin of which is still unclear, as no RHDV2 reservoir hosts were identified. After the RHDV2 emergence in 2010, viral RNA was detected in a few rodent species. Furthermore, RHDV2 was found to cause disease in some hare species resembling the disease in rabbits, evidencing the ability of the virus to cross the species barrier. In this study, through molecular, histopathologic, antigenic and morphological evidences, we demonstrate the presence and replication of RHDV2 in Eurasian badgers (Meles meles) found dead in the district of Santarém, Portugal, between March 2017 and January 2020. In these animals, we further classify the RHDV2 as a Lagovirus europaeus recombinant GI.4P-GI.2. Our results indicate that Meles meles is susceptible to RHDV2, developing systemic infection, and excreting the virus in the faeces. Given the high viral loads seen in several organs and matrices, we believe that transmission to the wild rabbit is likely. Furthermore, transmission electron microscopy data show the presence of calicivirus compatible virions in the nucleus of hepatocytes, which constitutes a paradigm shift for caliciviruses' replication cycle.

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.

reciprocalspaceship: a Python library for crystallographic data analysis.

Crystallography uses the diffraction of X-rays, electrons or neutrons by crystals to provide invaluable data on the atomic structure of matter, from single atoms to ribosomes. Much of crystallography's success is due to the software packages developed to enable automated processing of diffraction data. However, the analysis of unconventional diffraction experiments can still pose significant challenges - many existing programs are closed source, sparsely documented, or challenging to integrate with modern libraries for scientific computing and machine learning. Described here is reciprocalspaceship, a Python library for exploring reciprocal space. It provides a tabular representation for reflection data from diffraction experiments that extends the widely used pandas library with built-in methods for handling space groups, unit cells and symmetry-based operations. As is illustrated, this library facilitates new modes of exploratory data analysis while supporting the prototyping, development and release of new methods.

Chimeric VLPs Bearing VP60 from Two Serotypes of Rabbit Haemorrhagic Disease Virus Are Protective against Both Viruses.

The VP60 capsid protein from rabbit haemorrhagic disease virus (RHDV), the causative agent of one of the most economically important disease in rabbits worldwide, forms virus-like particles (VLPs) when expressed using heterologous protein expression systems such as recombinant baculovirus, yeasts, plants or mammalian cell cultures. To prevent RHDV dissemination, it would be beneficial to develop a bivalent vaccine including both RHDV GI.1- and RHDV GI.2-derived VLPs to achieve robust immunisation against both serotypes. In the present work, we developed a strategy of production of a dual-serving RHDV vaccine co-expressing the VP60 proteins from the two RHDV predominant serotypes using CrisBio technology, which uses Tricholusia ni insect pupae as natural bioreactors, which are programmed by recombinant baculovirus vectors. Co-infecting the insect pupae with two baculovirus vectors expressing the RHDV GI.1- and RHDV GI.2-derived VP60 proteins, we obtained chimeric VLPs incorporating both proteins as determined by using serotype-specific monoclonal antibodies. The resulting VLPs showed the typical size and shape of this calicivirus as determined by electron microscopy. Rabbits immunised with the chimeric VLPs were fully protected against a lethal challenge infection with the two RHDV serotypes. This study demonstrates that it is possible to generate a dual cost-effective vaccine against this virus using a single production and purification process, greatly simplifying vaccine manufacturing.

Molecular detection of myxoma virus in the environment of vaccinated rabbitries.

Myxoma virus (MYXV) is the aetiological agent of myxomatosis, a systemic, mostly lethal disease that affects European rabbits. Vaccination against it, although widespread, has not been completely effective and disease outbreaks still take place on farms which carry out vaccination programmes. Since some of these cases have been attributed to airborne transmission or the spread of the virus via inanimate vectors, the aims of this study were to determine MYXV contamination levels and distribution in the environment of vaccinated farms and to ascertain whether the detected virus corresponded to field strains. For that, environmental samples from several areas, tools and employees from four (three infected and one uninfected) rabbitries were taken and analysed by qPCR. MYXV was detected in the environment of all the infected farms, whereas all the samples from the non-infected farm were negative. Furthermore, all the positive samples contained viral DNA compatible with field strains of the virus. These results lead us to believe that the administration of currently available commercial vaccines does not prevent infected animals from shedding the field virus. Moreover, viral DNA was also found in items that are not in direct contact with the animals, which could play a role in the transmission of the infection throughout the farm and to other farms. Therefore, this study proves that current vaccination schemes on their own are not sufficient to prevent this disease and should be accompanied by adequate biosecurity measures.

Reverse Genetics System for Rabbit vesivirus.

Most caliciviruses are refractory to replication in cell culture and only a few members of the family propagate in vitro. Rabbit vesivirus (RaV) is unique due to its ability to grow to high titers in several animal and human cell lines. This outstanding feature makes RaV an ideal candidate for reverse genetics studies, an invaluable tool to understand the molecular basis of virus replication, the biological functions of viral genes and their roles in pathogenesis. The recovery of viruses from a cDNA clone is a prerequisite for reverse genetics studies. In this work, we constructed a RaV infectious cDNA clone using a plasmid expression vector, under the control of bacteriophage T7 RNA-polymerase promoter. The transfection of permissive cells with this plasmid DNA in the presence of T7 RNA-polymerase, provided in trans by a helper recombinant poxvirus, led to de novo synthesis of RNA transcripts that emulated the viral genome. The RaV progeny virus produced the typical virus-induced cytopathic effect after several passages of cell culture supernatants. Similarly, infectious RaV was recovered when the transcription step was performed in vitro, prior to transfection, provided that a 5'-cap structure was added to the 5' end of synthetic genome-length RNAs. In this work, we report an efficient and consistent RaV rescue system based on a cDNA transcription vector, as a tool to investigate calicivirus biology through reverse genetics.

Author Correction: Water sports could contribute to the translocation of ranaviruses.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

REP-X: An Evolution-guided Strategy for the Rational Design of Cysteine-less Protein Variants.

Site-specific labeling of proteins is often a prerequisite for biophysical and biochemical characterization. Chemical modification of a unique cysteine residue is among the most facile methods for site-specific labeling of proteins. However, many proteins have multiple reactive cysteines, which must be mutated to other residues to enable labeling of unique positions. This trial-and-error process often results in cysteine-free proteins with reduced activity or stability. Herein we describe a general methodology to rationally engineer cysteine-less proteins. Briefly, natural variation across orthologues is exploited to identify suitable cysteine replacements compatible with protein activity and stability. As a proof-of-concept, we recount the successful engineering of a cysteine-less mutant of the group II chaperonin from methanogenic archaeon Methanococcus maripaludis. A webapp, REP-X (Replacement at Endogenous Positions from eXtant sequences), which enables users to design their own cysteine-less protein variants, will make this rational approach widely available.

Myxomatosis and Rabbit Haemorrhagic Disease: A 30-Year Study of the Occurrence on Commercial Farms in Spain.

In this retrospective study, we describe the relative occurrence of clinical myxomatosis, and rabbit haemorrhagic disease (RHD), on 1714 commercial farms visited in Spain, between 1988 and 2018. We determined the annual prevalence based on 817 visits to 394 farms affected by myxomatosis. Myxomatosis was more prevalent from August to March, being lowest in June (3%) and highest in September (8.9%). With regard to RHD, we assessed 253 visits to 156 affected farms. We analyzed mean annual and monthly incidence. Two important RHD epidemics occurred; the first in 1988-1989 due to RHDV GI.1 (also known as RHDV), and the second from 2011 to 2013 due to RHDV GI.2 (RHDV2 or RHDVb). These epidemics occurred at times when effective vaccination had not been carried out. Relative monthly incidence in 2011-2018 was higher from April to August (p < 0.001). The results we obtained from 1404 necropsies on 102 farms did not clearly relate serosanguinous nasal discharge in rabbits with disease caused by GI.2 infection. We also assessed vaccination schedules used on 200 doe farms visited from the end of 2014 to 2018; 95.5% vaccinated against myxomatosis and 97.5% against RHD. Both diseases remain prevalent; however, effective vaccination has produced a steady decline in myxomatosis and RHDV GI.1 and GI.2 on-farm detection. The maintenance of high hygienic standards will be needed to continue and improve this control. However, further studies are required to investigate the causes of sustained virus presence and vaccine breaks.

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.