Emergence of Pathogenicity in Lagoviruses: Evolution from Pre-existing Nonpathogenic Strains or through a Species Jump?

Detection of a New Recombinant Rabbit Hemorrhagic Disease Virus 2 in China and Development of Virus-like Particle-Based Vaccine

Rabbit hemorrhagic disease virus (RHDV) is a very virulent virus of the genus Lagovirus causing severe and fatal hepatitis in the European rabbit (Oryctolagus cuniculus). RHDV has two distinct genotypes: GI.1 (RHDV) and GI.2 (RHDV2). The first RHDV2/GI.2 strain was identified as a recombinant virus between a non-pathogenic (GI.3P) and a pathogenic (GI.2) lagovirus, and the recombination is thought to have been a key mechanism in the emergence and evolution of RHDV2. Here, a new variant of RHDV2 was identified affecting domestic rabbits on Chinese farms, with a mortality rate of 70-80%. Phylogenetic analysis indicated that the nonstructural portion of this newly identified strain's genome clustered with the GI.1a variants. In contrast, the capsid gene shared the highest nucleotide identity of 97.9% with the North American GI.2 strains, suggesting a possible introduction in China of North American strains and recombination with the GI.1a strains circulating in China. We have produced a recombinant vaccine using the first Chinese GI.2 strain, SC2020/0401, by cloning the vp60 gene into a baculovirus expression vector. Virus-like particles (VLPs) were then produced in Sf9 insect cells, and a challenge study was performed. Rabbits immunized with the VLP vaccine survived 7 d after being challenged with the new virus. The results indicate that commercial vaccines are urgently required in China to control the circulation of RHDV2 variants.

New Insights into Rabbit Viral Diseases

Viruses are responsible for many devastating rabbit diseases that impact their health and welfare and put their conservation and economic revenue at risk [...].

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.

Serological characterisation of Lagovirus virus-like particles originating from native and mutated VP60 of rabbit haemorrhagic disease virus 2 and European brown hare syndrome virus

Introduction

Since lagoviruses cannot be cultivated in vitro, using expression systems is an alternative and promising way of producing diagnostic viral antigens. It opens up their use as active immunogens for vaccine production.

Material and Methods

Virus-like particles (VLPs) were produced in a baculovirus expression system in Spodoptera frugiperda 9 (Sf9) insect cells based on wild-type and mutated variants of the virus capsid VP60 protein from a Polish strain of European brown hare syndrome virus (EBHSV) and wild-type and mutated versions of this protein from a Polish strain of rabbit haemorrhagic disease virus 2 (RHDV2). The mutations were the substitution of an arginylglycylaspartic acid (Arg-Gly-Asp/RGD) motif in the P2 subdomain and, in the S or P2 domain, the substitution of three lysines. The VLPs were purified with sucrose gradient ultracentrifugation.

Results

Protein production was confirmed by Western blot analysis using rabbit or hare sera and ELISA tests with different types of monoclonal antibody. The haemagglutination properties of some VLPs were also evaluated. Electron microscopy of wild-type EBHSV, wild-type RHDV2 and the four VP60 variants produced in this experiment revealed the formation of characteristic VLP structures.

Conclusion

For the first time, mutated VLPs of RHDV2 with an RGD motif in the VP60 sequence were obtained, which could potentially be used to deliver cargo to eukaryotic cells. Virus-like particles based on the VP60 proteins of EBHSV and RHDV with a three-lysine substitution in the S or P2 domains were also obtained. Potential exists for VLPs of EBHSV and RHDV2 as vaccine candidates.

A new HaCV-EBHSV recombinant lagovirus circulating in European brown hares (Lepus europaeus) from Catalonia, Spain

In 2020/2021, several European brown hare syndrome virus (EBHSV) outbreaks were recorded in European hares (Lepus europaeus) from Catalonia, Spain. Recombination analysis combined with phylogenetic reconstruction and estimation of genetic distances of the complete coding sequences revealed that 5 strains were recombinants. The recombination breakpoint is located within the non-structural protein 2C-like RNA helicase (nucleotide position ~ 1889). For the genomic fragment upstream of the breakpoint, a non-pathogenic EBHSV-related strain (hare calicivirus, HaCV; GII.2) was the most closely related sequence; for the rest of the genome, the most similar strains were the European brown hare syndrome virus (EBHSV) strains recovered from the same 2020/2021 outbreaks, suggesting a recent origin. While the functional impact of the atypical recombination breakpoint remains undetermined, the novel recombinant strain was detected in different European brown hare populations from Catalonia, located 20-100 km apart, and seems to have caused a fatal disease both in juvenile and adult animals, confirming its viability and ability to spread and establish infection. This is the first report of a recombination event involving HaCV and EBHSV and, despite the recombination with a non-pathogenic strain, it appears to be associated with mortality in European brown hares, which warrants close monitoring.

On the virome's hidden diversity: lessons from RHDV

Emerging infectious diseases are a major challenge to human and animal health. While predicting the emergence of pathogens is complex, the advent of high-throughput sequencing technologies has allowed the rapid identification of unknown microbiology diversity within organisms. Here, we discuss an example of a metatranscriptomics output to decipher viral evolution.

Novel Recombinant Rabbit Hemorrhagic Disease Virus 2 (RHDV2) is Circulating in China within 12 Months after Original RHDV2 Arrival

Rabbit hemorrhagic disease (RHD) causes lethal fulminant hepatitis in rabbits. Two different genotypes (GI.1 and GI.2) responsible for RHD are reported. GI.2 was first detected in France in 2010 and subsequently spread to other countries in Europe. In April 2020, GI.2 was detected in China. In this study, we report a novel recombinant strain of fatal rabbit hemorrhagic disease virus 2 (RHDV2 or GI.2) detected from domestic rabbits in three provinces in China in 2020-2021. Full-length genomic analysis has revealed that the recombinant virus contained an RHDV2 capsid gene and nonstructural genes from an unclassified lagovirus genotype. This type of virus emerged and circulated throughout China within a year after the initial detection of the original RHDV2. Compared with the original strain, the new virus showed a longer infected time and lower mortality rate but almost the same viral load at the moribund stage of infection. This might have resulted in high virus contamination in the environment, facilitating virus transmission. As the consequences of the presence of novel recombinant strains are unpredictable, the circulation of the novel variant in the population should be carefully monitored in China.

Genetic Characteristics and Phylogeographic Dynamics of Lagoviruses, 1988-2021

Rabbit haemorrhagic disease virus (RHDV), European brown hare syndrome virus (EBHSV), rabbit calicivirus (RCV), and hare calicivirus (HaCV) belong to the genus Lagovirus of the Caliciviridae family that causes severe diseases in rabbits and several hare (Lepus) species. Previously, Lagoviruses were classified into two genogroups, e.g., GI (RHDVs and RCVs) and GII (EBHSV and HaCV) based on partial genomes, e.g., VP60 coding sequences. Herein, we provide a robust phylogenetic classification of all the Lagovirus strains based on full-length genomes, grouping all the available 240 strains identified between 1988 and 2021 into four distinct clades, e.g., GI.1 (classical RHDV), GI.2 (RHDV2), HaCV/EBHSV, and RCV, where the GI.1 clade is further classified into four (GI.1a-d) and GI.2 into six sub-clades (GI.2a-f). Moreover, the phylogeographic analysis revealed that the EBHSV and HaCV strains share their ancestor with the GI.1, while the RCV shares with the GI.2. In addition, all 2020-2021 RHDV2 outbreak strains in the USA are connected to the strains from Canada and Germany, while RHDV strains isolated in Australia are connected with the USA-Germany haplotype RHDV strain. Furthermore, we identified six recombination events in the VP60, VP10, and RNA-dependent RNA polymerase (RdRp) coding regions using the full-length genomes. The amino acid variability analysis showed that the variability index exceeded the threshold of 1.00 in the ORF1-encoded polyprotein and ORF2-encoded VP10 protein, respectively, indicating significant amino acid drift with the emergence of new strains. The current study is an update of the phylogenetic and phylogeographic information of Lagoviruses that may be used to map the evolutionary history and provide hints for the genetic basis of their emergence and re-emergence.

Rabbit haemorrhagic disease virus 2 from Singapore 2020 outbreak revealed an Australian recombinant variant

Rabbit haemorrhagic disease (RHD) is a significant and debilitating viral disease affecting lagomorphs. In September 2020, Singapore reported its first cases of RHD virus (RHDV) infection in domesticated rabbits. The initial findings reported that the outbreak strain belonged to genotype GI.2 (RHDV2/RHDVb), and epidemiological investigations could not identify the definitive source of the virus origin. Further recombination detection and phylogenetic analyses of the Singapore outbreak strain revealed that the RHDV was a GI.2 structural (S)/GI.4 non-structural (NS) recombinant variant. Sequence analyses on the National Centre for Biotechnology Information (NCBI) database showed high homology to recently emerged Australian variants, which were prevalent in local Australian lagomorph populations since 2017. Time-structured and phylogeographic analyses for the S and NS genes revealed a close genetic relationship between the Singapore RHDV strain and the Australian RHDV variants. More thorough epidemiological inquiries are necessary to ascertain how an Australian RHDV was introduced into the Singapore rabbit population, and opportune development of RHDV diagnostics and vaccines will be important to safeguard lagomorphs from future RHDV infection and disease management.

Pathogenicity of the newly emerged Lagovirus europaeus GI.2 strain in China in experimentally infected rabbits

In April 2020, rabbit hemorrhagic virus type 2 (Lagovirus europaeus GI.2), which causes highly infectious fatal rabbit hemorrhagic disease, was emerged in China. The phylogenetic analyses of the complete genome sequence of GI.2 showed that it belonged to the non-recombinant GI.3/GI.2 genotype. However, the pathogenicity of this GI.2 strain differed from that of early typical GI.2 strains in Europe. To prevent the spread of the new strain in China, its pathogenicity urgently needs to be studied. Thus, viral shedding and distribution as well as clinical symptoms, histopathological changes, and serum cytokines were studied in experimentally GI.2/SC2020-infected rabbit adults and kits. The kit group showed a shorter survival time after the challenge than the adult group did. The mortality rate was higher in the kits (80 %) than in the adults (30 %). Viral RNA could be detected in both nasal and fecal swabs, and the main dissemination route appeared to be the fecal route. Viral RNA rapidly increased in the blood of the adults and kits at 6 h post-infection, indicating that blood viral load testing can be used for early diagnosis. The most affected organs were the liver and spleen, and the lesions were more severe in the kits than in the adults. The liver contained the highest viral RNA levels. Moreover, serum interleukin (IL)-6, IL-8, IL-10, and tumor necrosis factor-alpha levels were increased in the infected rabbits. In conclusion, our findings will help to understand the evolutionary trends and pathogenic characteristics of GI.2 strains in China.

Immunity against Lagovirus europaeus and the Impact of the Immunological Studies on Vaccination

In the early 1980s, a highly contagious viral hemorrhagic fever in rabbits (Oryctolagus cuniculus) emerged, causing a very high rate of mortality in these animals. Since the initial occurrence of the rabbit hemorrhagic disease virus (RHDV), several hundred million rabbits have died after infection. The emergence of genetically-different virus variants (RHDV GI.1 and GI.2) indicated the very high variability of RHDV. Moreover, with these variants, the host range broadened to hare species (Lepus). The circulation of RHDV genotypes displays different virulences and a limited induction of cross-protective immunity. Interestingly, juvenile rabbits (<9 weeks of age) with an immature immune system display a general resistance to RHDV GI.1, and a limited resistance to RHDV GI.2 strains, whereas less than 3% of adult rabbits survive an infection by either RHDV GI.1. or GI.2. Several not-yet fully understood phenomena characterize the RHD. A very low infection dose followed by an extremely rapid viral replication could be simplified to the induction of a disseminated intravascular coagulopathy (DIC), a severe loss of lymphocytes—especially T-cells—and death within 36 to 72 h post infection. On the other hand, in animals surviving the infection or after vaccination, very high titers of RHDV-neutralizing antibodies were induced. Several studies have been conducted in order to deepen the knowledge about the virus’ genetics, epidemiology, RHDV-induced pathology, and the anti-RHDV immune responses of rabbits in order to understand the phenomenon of the juvenile resistance to this virus. Moreover, several approaches have been used to produce efficient vaccines in order to prevent an infection with RHDV. In this review, we discuss the current knowledge about anti-RHDV resistance and immunity, RHDV vaccination, and the further need to establish rationally-based RHDV vaccines.

Recombination at the emergence of the pathogenic rabbit haemorrhagic disease virus Lagovirus europaeus/GI.2

Rabbit haemorrhagic disease is a viral disease that emerged in the 1980s and causes high mortality and morbidity in the European rabbit (Oryctolagus cuniculus). In 2010, a new genotype of the rabbit haemorrhagic disease virus emerged and replaced the former circulating Lagovirus europaeus/GI.1 strains. Several recombination events have been reported for the new genotype Lagovirus europaeus/GI.2, with pathogenic (variants GI.1a and GI.1b) and benign (genotype GI.4) strains that served as donors for the non-structural part while GI.2 composed the structural part; another recombination event has also been described at the p16/p23 junction involving GI.4 strains. In this study, we analysed new complete coding sequences of four benign GI.3 strains and four GI.2 strains. Phylogenetic and recombination detection analyses revealed that the first GI.2 strains, considered as non-recombinant, resulted from a recombination event between GI.3 and GI.2, with GI.3 as the major donor for the non-structural part and GI.2 for the structural part. Our results indicate that recombination contributed to the emergence, persistence and dissemination of GI.2 as a pathogenic form and that all described GI.2 strains so far are the product of recombination. This highlights the need to study full-genomic sequences of lagoviruses to understand their emergence and evolution.

Retrospective Analysis Shows That Most RHDV GI.1 Strains Circulating Since the Late 1990s in France and Sweden Were Recombinant GI.3P-GI.1d Strains

Recombination is one of the major sources of genetic variation in viruses. RNA viruses, such as rabbit hemorrhagic disease virus (RHDV), are among the viruses with the highest recombination rates. Several recombination events have been described for RHDV, mostly as a consequence of their genomic architecture. Here, we undertook phylogenetic and recombination analyses of French and Swedish RHDV strains from 1994 to 2016 and uncovered a new intergenotypic recombination event. This event occurred in the late 1990s/early 2000s and involved nonpathogenic GI.3 strains as donors for the nonstructural part of the genome of these recombinants, while pathogenic GI.1d strains contributed to the structural part. These GI.3P-GI.1d recombinant strains did not entirely replace GI.1d (nonrecombinant) strains, but became the dominant strains in France and Sweden, likely due to a fitness advantage associated with this genomic architecture. GI.3P-GI.1d (P stands for polymerase) strains persisted until 2013 and 2016 in Sweden and France, respectively, and cocirculated with the new genotype GI.2 in France. Since strains from the first GI.2 outbreaks were GI.3P-GI.2, we hypothesize that GI.3P-GI.1d could be the parental strain. Our results confirm the outstanding recombination ability of RHDV and its importance in the evolution of lagoviruses, which was only revealed by studying complete genomic sequences.

Genetic diversity and evolution of Hare Calicivirus (HaCV), a recently identified lagovirus from Lepus europaeus

First recognized as highly pathogenic viruses, hare lagoviruses belonging to genotype GII.1 (EBHSV) infect various Lepus species. Genetically distinct benign lagoviruses (Hare Calicivirus, HaCV) have recently been identified but few data have been available so far on these strains. The analysis of 199 samples from hunted hares collected throughout France allowed the detection of 20 HaCV and showed that they were widely distributed in this country. Ten HaCV capsid protein gene sequences were characterized. A first HaCV capsid protein structural model was proposed, revealing a global structure similar to that of a pathogenic GII.1 strain. The HaCV sequences showed an even higher genetic diversity than previously appreciated, with the characterization of two genotypes (GII.2, GII.3) and several additional putative genotypes. The most recent common ancestor for HaCV VP60 gene was estimated to be much older than that for GII.1 pathogenic strains. These results give new insights into the phylogenetic relationships of HaCV within the Lagovirus genus.

Sequencing of VDJ genes in Lepus americanus confirms a correlation between VHn expression and the leporid species continent of origin

Leporid VH genes used in the generation of their primary antibody repertoire exhibit highly divergent lineages. For the European rabbit (Oryctolagus cuniculus) four VHa lineages have been described, the a1, a2, a3 and a4. Hares (Lepus spp.) and cottontail (Sylvilagus floridanus) express one VHa lineage each, the a2L and the a5, respectively, along with a more ancient lineage, the Lepus spp. sL and S. floridanus sS. Both the European rabbit and the Lepus europaeus use a third lineage, VHn, in a low proportion of their VDJ rearrangements. The VHn genes are a conserved ancestral polymorphism that is being maintained in the leporid genome.Their usage in a low proportion of VDJ rearrangements by both European rabbit and L. europaeus but not S. floridanus has been argued to be a remnant of an ancient European leporid immunologic response to pathogens. To address this hypothesis, in this study we sequenced VDJ rearranged genes for another North American leporid, L. americanus. Our results show that L. americanus expressed these genes less frequently and in a highly modified fashion compared to the European Lepus species. Our results suggest that the American leporid species use a different VH repertoire than the European species which may be related with an immune adaptation to different environmental conditions, such as different pathogenic agents.

First Complete Genome Sequence of a Hare Calicivirus Strain Isolated from Lepus europaeus

The first full-genome sequence of a hare calicivirus (HaCV), recently characterized as a novel member of the Caliciviridae, is described. This presumed nonpathogenic lagovirus is 7,433 nucleotides long, shows the same genomic organization as that of other lagoviruses, and has the highest nucleotide identity (79%) with pathogenic European brown hare syndrome viruses.

The first full-genome sequence of a hare calicivirus (HaCV), recently characterized as a novel member of the Caliciviridae, is described. This presumed nonpathogenic lagovirus is 7,433 nucleotides long, shows the same genomic organization as that of other lagoviruses, and has the highest nucleotide identity (79%) with pathogenic European brown hare syndrome viruses.

Identification of a new European rabbit IgA with a serine-rich hinge region

In mammals, the most striking IgA system belongs to Lagomorpha. Indeed, 14 IgA subclasses have been identified in European rabbits, 11 of which are expressed. In contrast, most other mammals have only one IgA, or in the case of hominoids, two IgA subclasses. Characteristic features of the mammalian IgA subclasses are the length and amino acid sequence of their hinge regions, which are often rich in Pro, Ser and Thr residues and may also carry Cys residues. Here, we describe a new IgA that was expressed in New Zealand White domestic rabbits of IGHVa1 allotype. This IgA has an extended hinge region containing an intriguing stretch of nine consecutive Ser residues and no Pro or Thr residues, a motif exclusive to this new rabbit IgA. Considering the amino acid properties, this hinge motif may present some advantage over the common IgA hinge by affording novel functional capabilities. We also sequenced for the first time the IgA14 CH2 and CH3 domains and showed that IgA14 and IgA3 are expressed.

First complete genome sequence of a European non-pathogenic rabbit calicivirus (lagovirus GI.3)

We report the full genome sequence of the non-pathogenic rabbit lagovirus Lagovirus europaeus/GI.3/O cun/FR/2006/06-11 (GI.3/06-11), collected from a healthy French domestic rabbit in 2006, and initially described as 06-11 strain. The sequence reveals a genomic organization similar to lagoviruses. It was 7,436 bases long and contained two open reading frames (ORF). A dipeptide variation at the potential p23/2C-like helicase cleavage site (EE instead of ED) was observed, a feature only shared with non-recombinant pathogenic lagoviruses in GI.2 and with two European brown hare syndrome viruses (EBHSV) collected in 1982 in Sweden. GI.3/06-11 has only one initiation codon at the beginning of the ORF2 like the avirulent Italian rabbit calicivirus (RCV) and EBHSV. Previous genetic analyses based on the capsid gene sequences showed that GI.3/06-11 was closer to the RCV and pathogenic lagoviruses GI.1 strains than other lagoviruses. This study, by revealing that GI.3/06-11 genome sequence significantly clustered with pathogenic GI.2 strains, gives prominence of new genetic relationship among lagoviruses and should contribute to understand the emergence of pathogenic strains.

Detecting European Rabbit ( Oryctolagus cuniculus) Disease Outbreaks by Monitoring Digital Media

Digital media and digital search tools offer simple and effective means to monitor for pathogens and disease outbreaks in target organisms. Using tools such as Rich Site Summary feeds, and Google News and Google Scholar specific key word searches, international digital media were actively monitored from 2012 to 2016 for pathogens and disease outbreaks in the taxonomic order Lagomorpha, with a specific focus on the European rabbit ( Oryctolagus cuniculus). The primary objective was identifying pathogens for assessment as potential new biocontrol agents for Australia's pest populations of the European rabbit. A number of pathogens were detected in digital media reports. Additional benefits arose in the regular provision of case reports and research on myxomatosis and rabbit haemorrhagic disease virus that assisted with current research.

Host-Specific Glycans Are Correlated with Susceptibility to Infection by Lagoviruses, but Not with Their Virulence

Rabbit hemorrhagic disease virus (RHDV) and European brown hare syndrome virus (EBHSV) are two lagoviruses from the family Caliciviridae that cause fatal diseases in two leporid genera, Oryctolagus and Lepus, respectively. In the last few years, several examples of host jumps of lagoviruses among leporids were recorded. In addition, a new pathogenic genotype of RHDV emerged, and many nonpathogenic strains of lagoviruses have been described. The molecular mechanisms behind host shifts and the emergence of virulence are unknown. Since RHDV uses glycans of the histo-blood group antigen type as attachment factors to initiate infection, we studied if glycan specificities of the new pathogenic RHDV genotype, nonpathogenic lagoviruses, and EBHSV potentially play a role in determining the host range and virulence of lagoviruses. We observed binding to A, B, or H antigens of the histo-blood group family for all strains known to primarily infect European rabbits (Oryctolagus cuniculus), which have recently been classified as GI strains. However, we could not explain the emergence of virulence, since similar glycan specificities were found in several pathogenic and nonpathogenic strains. In contrast, EBHSV, recently classified as GII.1, bound to terminal β-linked N-acetylglucosamine residues of O-glycans. Expression of these attachment factors in the upper respiratory and digestive tracts in three lagomorph species (Oryctolagus cuniculus, Lepus europaeus, and Sylvilagus floridanus) showed species-specific patterns regarding susceptibility to infection by these viruses, indicating that species-specific glycan expression is likely a major contributor to lagovirus host specificity and range.IMPORTANCE Lagoviruses constitute a genus of the family Caliciviridae comprising highly pathogenic viruses, RHDV and EBHSV, that infect rabbits and hares, respectively. Recently, nonpathogenic strains were discovered and new pathogenic strains have emerged. In addition, host jumps between lagomorphs have been observed. The mechanisms responsible for the emergence of pathogenicity and host species range are unknown. Previous studies showed that RHDV strains attach to glycans expressed in the upper respiratory and digestive tracts of rabbits, the likely portals of virus entry. Here, we studied the glycan-binding properties of novel pathogenic and nonpathogenic strains looking for a link between glycan binding and virulence or between glycan specificity and host range. We found that glycan binding did not correlate with virulence. However, expression of glycan motifs in the upper respiratory and digestive tracts of lagomorphs revealed species-specific patterns associated with the host ranges of the virus strains, suggesting that glycan diversity contributes to lagovirus host ranges.

Characterization of old RHDV strains by complete genome sequencing identifies a novel genetic group

Rabbit hemorrhagic disease (RHD) is a veterinary disease that affects the European rabbit and has a significant economic and ecological negative impact. In Portugal, rabbit hemorrhagic disease virus (RHDV) was reported in 1989 and still causes enzootic outbreaks. Several recombination events have been detected in RHDV strains, including in the first reported outbreak. Here we describe the occurrence of recombination in RHDV strains recovered from rabbit and Iberian hare samples collected in the mid-1990s in Portugal. Characterization of full genomic sequences revealed the existence of a single recombination breakpoint at the boundary of the non-structural and the structural encoding regions, further supporting the importance of this region as a recombination hotspot in lagoviruses. Phylogenetic analysis showed that in the structural region, the recombinant strains were similar to pathogenic G1 strains, but in the non-structural region they formed a new group that diverged ~13% from known strains. No further reports of such group exist, but this recombination event was also detected in an Iberian hare that was associated with the earliest species jump in RHDV. Our results highlight the importance of the characterization of full genomes to disclose RHDV evolution and show that lagoviruses’ diversity has been significantly undersampled.

Proposal for a unified classification system and nomenclature of lagoviruses

Lagoviruses belong to the Caliciviridae family. They were first recognized as highly pathogenic viruses of the European rabbit (Oryctolagus cuniculus) and European brown hare (Lepus europaeus) that emerged in the 1970-1980s, namely, rabbit haemorrhagic disease virus (RHDV) and European brown hare syndrome virus (EBHSV), according to the host species from which they had been first detected. However, the diversity of lagoviruses has recently expanded to include new related viruses with varying pathogenicity, geographic distribution and host ranges. Together with the frequent recombination observed amongst circulating viruses, there is a clear need to establish precise guidelines for classifying and naming lagovirus strains. Therefore, here we propose a new nomenclature based on phylogenetic relationships. In this new nomenclature, a single species of lagovirus would be recognized and called Lagovirus europaeus. The species would be divided into two genogroups that correspond to RHDV- and EBHSV-related viruses, respectively. Genogroups could be subdivided into genotypes, which could themselves be subdivided into phylogenetically well-supported variants. Based on available sequences, pairwise distance cutoffs have been defined, but with the accumulation of new sequences these cutoffs may need to be revised. We propose that an international working group could coordinate the nomenclature of lagoviruses and any proposals for revision.

Phylogenetic analysis of rabbit haemorrhagic disease virus (RHDV) strains isolated in Poland

The aim of this study was to characterise the nucleotide and amino acid sequence of complete genomes (7.5 kb) from RHDV strains isolated in Poland and estimate the genetic variability in different elements of the viral RNA. In addition, the sequence of Polish RHDV isolates isolated from 1988-2015 was compared with the sequences of other European RHDV, including the RHDVa and RHDV2/RHDVb subtypes. The complete sequence was developed by the compilation of partial nucleotide sequences. This sequence consisted of approximately 7428 nucleotides. For comparison of nucleotide sequences and the development of phylogenetic trees of Polish RHDV isolates and reference RHDV strains representing the main phylogenetic groups of classical RHDV, RHDVa and RHDV2 as well as the non-pathogenic rabbit lagovirus RCV, the BLAST software with blastn and MEGA6 with neighbour-joining method was applied. The complete nucleotide sequence of Polish isolates of RHDV has also been entered into GenBank. For comparative analysis, nineteen complete sequences representing the main RHDV genetic types available in GenBank were used. The results of phylogenetic analysis of Polish RHDV strains reveals the presence of three classical RHDV genogroups (G2, G4 and G5) and an RHDVa variant (G6). The oldest RHDV isolates (KGM 1988, PD 1989 and MAL 1994) belong to genogroup G2. It can be assumed that the elimination of these strains from the environment probably occurred at the turn of 1994 and 1995. Genogroup G2 was replaced by the phylogenetically younger BLA 1994 and OPO 2004 strains from genogroup G4, which probably originated from the G3 lineage, represented by the Italian strains BS89. The last representatives of classical RHDV in Poland are isolates GSK 1988 and ZD0 2000 from genogroup G5. A single clade contains the Polish RHDV strains from 2004-2015 (GRZ 2004, KRY 2004, L145 2004, W147 2005, SKO 2013, GLE 2013, RED1 2013, STR 2012, STR2 2013, STR 2014, BIE 2015) identified as RHDVa, which clustered into genogroup G6, as represented by the RHDV strain Triptis 1996. All recent isolated RHDV isolates belong exclusively to RHDVa and no RHDV2 was diagnosed in Poland.

An update on the rabbit hemorrhagic disease virus (RHDV) strains circulating in Portugal in the 1990s: earliest detection of G3-G5 and G6

Rabbit hemorrhagic disease virus (RHDV) causes devastating effects on European rabbit (Oryctolagus cuniculus) populations in the Iberian Peninsula. According to the information available, only genogroup 1 strains were circulating in Iberian wild rabbits until 2011; the antigenic variant G6 has been sporadically detected in rabbitries since 2007. Here, we show for the first time that G3-G5 strains were already present in mainland Portugal in 1998 and that G6 has been circulating since at least 1999. Moreover, we report a G3-G5 strain from the Azores collected in 1998, which is the likely ancestor of Azorean G3-G5like strains. These observations improve the current knowledge on RHDV epidemiology in the Iberian Peninsula and the Azores.

Experimental transmission of rabbit haemorrhagic disease virus (RHDV) from rabbit to wild mice (Mus spretus and Apodemus sylvaticus) under laboratory conditions

Rabbit haemorrhagic disease (RHD) is a highly lethal and contagious viral disease that produces haemorrhagic lesions in liver and lungs of domestic and wild rabbits (Oryctolagus cuniculus). This study investigates the transmission of RHDV from infected rabbits to mice, based on the detection of viral RNA. Sixteen wild mice (Mus spretus, n=12 and Apodemus sylvaticus, n=4) were put in contact with nine rabbits inoculated with RHDV. No mice died following exposure to RHDV-infected rabbits or developed macroscopic haemorrhagic lesions. On the fourth day of contact, RHDV was detected by RT-PCR in the faeces of three of the four mice killed and in the livers of two of them. Three days after contact period with the inoculated rabbits (7th day of the experiment), RHDV was detected by RT-PCR in 100% (n=4) of the faeces and 50% (n=2) of the livers of euthanized animals. Ten days after contact period (14th day of the experiment), RHDV was not detected in the faeces or liver from any of the mice euthanized. However, 64days after contact period, RHDV was detected in the faeces of one mouse (1 of 4). We demonstrate cross-species transmission of RHDV-RNA from rabbit to rodent and the capability of RHDV-RNA to persist in mice for at least 10days after contact, and potentially up to two months, although viral replication within the rodent and/or infectivity was not evaluated in the present study.

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.

Comparative analysis of rabbit hemorrhagic disease virus strains originating from outbreaks in the Russian Federation

Since the first introduction of rabbit hemorrhagic disease (RHD) in 1986, disease outbreaks have been continuously reported in different regions of Russia. Despite extensive vaccination, sporadic RHD cases are still reported. Here, we examine eleven RHDV strains originating from disease outbreaks occurring between 2003 and 2012 and one widely used vaccine strain. Notable phenotypic and genetic heterogeneity among RHDV strains was observed. The RHDV strains Tambov-2010, Perm-2010, Manihino-09 showed different hemagglutinating activity (HA) at 4 °C and room temperature. While all RHDV field strains were identified as hemagglutinating virulent viruses of the RHDVa variant, the vaccine strain was assigned as a "classical" RHDV. These data indicate that since 2003, RHDVa has become the predominant variant circulating in Russia.

An overview of the lagomorph immune system and its genetic diversity

Our knowledge of the lagomorph immune system remains largely based upon studies of the European rabbit (Oryctolagus cuniculus), a major model for studies of immunology. Two important and devastating viral diseases, rabbit hemorrhagic disease and myxomatosis, are affecting European rabbit populations. In this context, we discuss the genetic diversity of the European rabbit immune system and extend to available information about other lagomorphs. Regarding innate immunity, we review the most recent advances in identifying interleukins, chemokines and chemokine receptors, Toll-like receptors, antiviral proteins (RIG-I and Trim5), and the genes encoding fucosyltransferases that are utilized by rabbit hemorrhagic disease virus as a portal for invading host respiratory and gut epithelial cells. Evolutionary studies showed that several genes of innate immunity are evolving by strong natural selection. Studies of the leporid CCR5 gene revealed a very dramatic change unique in mammals at the second extracellular loop of CCR5 resulting from a gene conversion event with the paralogous CCR2. For the adaptive immune system, we review genetic diversity at the loci encoding antibody variable and constant regions, the major histocompatibility complex (RLA) and T cells. Studies of IGHV and IGKC genes expressed in leporids are two of the few examples of trans-species polymorphism observed outside of the major histocompatibility complex. In addition, we review some endogenous viruses of lagomorph genomes, the importance of the European rabbit as a model for human disease studies, and the anticipated role of next-generation sequencing in extending knowledge of lagomorph immune systems and their evolution.

Genetic diversity comparison of the DQA gene in European rabbit (Oryctolagus cuniculus) populations

The European rabbit (Oryctolagus cuniculus) natural populations within the species native region, the Iberian Peninsula, are considered a reservoir of genetic diversity. Indeed, the Iberia was a Pleistocene refuge to the species and currently two subspecies are found in the peninsula (Oryctolagus cuniculus cuniculus and Oryctolagus cuniculus algirus). The genes of the major histocompatibility complex (MHC) have been substantially studied in wild populations due to their exceptional variability, believed to be pathogen driven. They play an important function as part of the adaptive immune system affecting the individual fitness and population viability. In this study, the MHC variability was assessed by analysing the exon 2 of the DQA gene in several European rabbit populations from Portugal, Spain and France and in domestic breeds. Twenty-eight DQA alleles were detected, among which 18 are described for the first time. The Iberian rabbit populations are well differentiated from the French population and domestic breeds. The Iberian populations retained the higher allelic diversity with the domestic breeds harbouring the lowest; in contrast, the DQA nucleotide diversity was higher in the French population. Signatures of positive selection were detected in four codons which are putative peptide-binding sites and have been previously detected in other mammals. The evolutionary relationships showed instances of trans-species polymorphism. Overall, our results suggest that the DQA in European rabbits is evolving under selection and genetic drift.