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

Two decades of occurrence of non-pathogenic rabbit lagoviruses in Italy and their genomic characterization

Lagoviruses are viruses of the Caliciviridae family affecting lagomorphs. Both pathogenic and non-pathogenic lagoviruses affect the European rabbit (Oryctolagus cuniculus), and they are phylogenetically distinguished. Rabbit Hemorrhagic Disease Virus (RHDV/GI.1) and Rabbit Hemorrhagic Disease Virus-2 (RHDV-2/GI.2) belong to the first group, while in the second group, several genotypes of Rabbit Calicivirus (RCV/GI.3-GI.4) are present. The first RCV strain was described in Italy in 1996, and since then, several RCV strains have been characterised in Europe and Australia. RCVs, different from the pathogenic hepatotropic RHDVs, have an enteric tropism and could be identified from the duodenum/intestine and faeces. This study aimed firstly to indirectly show through a seroepidemiological survey from 1998 to 2008 the circulation of RCVs strains in rabbit farms and then to genetically characterise RCV strains diagnosed in Italy in faecal and intestinal samples of wild and farmed rabbits collected in various regions in the following years (2000-2022). Of 262 analysed samples, 69 resulted in RT-PCR positive for lagovirus but negative for RHDV. Eleven RCV strains were characterised by complete vp60 sequencing. Phylogenetic analysis showed that the Italian RCV strains are grouped in European (RCV_E1/GI.3) and Australian (RCV_E2/GI.4) RCV clusters, with an estimated country prevalence of 26%. Based on the proposed genotype classification, considering the nucleotide differences of vp60 higher than 15%, we can hypothesise that two other genotypes, GI.5 and GI.6, might exist within the cluster of non-pathogenic viruses.

Surveillance of Wildlife Viruses: Insights from South Australia's Monitoring of Rabbit Haemorrhagic Disease Virus (RHDV GI.1 and GI.2)

Surveillance of wildlife virus impacts can be passive or active. Both approaches have their strengths and weaknesses, especially regarding cost and knowledge that can be gained. Monitoring of rabbit haemorrhagic disease virus (GI.1 and GI.2) in South Australia has utilised both strategies and their methods and gained insights are discussed. Active strategies to monitor the continuing impact of rabbit haemorrhagic disease virus 2 (GI.2) on susceptible lagomorphs in countries such as the USA, Mexico, South Africa, Spain, France and Portugal are encouraged to gain critical insights into the evolution, spread and impact of this virus. Furthermore, there are lessons here for the international monitoring of diseases in wildlife, particularly where there is a risk of them becoming zoonotic.

Characterisation of Lagovirus europaeus GI-RHDVs (Rabbit Haemorrhagic Disease Viruses) in Terms of Their Pathogenicity and Immunogenicity

Rabbit haemorrhagic disease viruses (RHDV) belong to the family Caliciviridae, genus Lagovirus europaeus, genogroup GI, comprising four genotypes GI.1-GI.4, of which the genotypes GI.1 and GI.2 are pathogenic RHD viruses, while the genotypes GI.3 and GI.4 are non-pathogenic RCV (Rabbit calicivirus) viruses. Among the pathogenic genotypes GI.1 and GI.2 of RHD viruses, an antigenic variant of RHDV, named RHDVa-now GI.1a-RHDVa, was distinguished in 1996; and in 2010, a variant of RHDV-named RHDVb, later RHDV2 and now GI.2-RHDV2/b-was described; and recombinants of these viruses were registered. Pathogenic viruses of the genotype GI.1 were the cause of a disease described in 1984 in China in domestic (Oryctolagus (O.) cuniculus domesticus) and wild (O. cuniculus) rabbits, characterised by a very rapid course and a mortality rate of 90-100%, which spread in countries all over the world and which has been defined since 1989 as rabbit haemorrhagic disease. It is now accepted that GI.1-RHDV, including GI.1a-RHDVa, cause the predetermined primary haemorrhagic disease in domestic and wild rabbits, while GI.2-RHDV2/b cause it not only in rabbits, including domestic rabbits' young up to 4 weeks and rabbits immunised with rabbit haemorrhagic disease vaccine, but also in five various species of wild rabbits and seven different species of hares, as well as wild ruminants: mountain muskoxen and European badger. Among these viruses, haemagglutination-positive, doubtful and harmful viruses have been recorded and described and have been shown to form phylogenogroups, immunotypes, haematotypes and pathotypes, which, together with traits that alter and expand their infectious spectrum (rabbit, hare, wild ruminant, badger and various rabbit and hare species), are the determinants of their pathogenicity (infectivity) and immunogenicity and thus shape their virulence. These relationships are the aim of our consideration in this article.

First Detection and Circulation of RHDV2 in New Zealand

Rabbit haemorrhage disease virus 2 (RHDV2) is a highly pathogenic lagovirus that causes lethal disease in rabbits and hares (lagomorphs). Since its first detection in Europe in 2010, RHDV2 has spread worldwide and has been detected in over 35 countries so far. Here, we provide the first detailed report of the detection and subsequent circulation of RHDV2 in New Zealand. RHDV2 was first detected in New Zealand in 2018, with positive samples retrospectively identified in December 2017. Subsequent time-resolved phylogenetic analysis suggested a single introduction into the North Island between March and November 2016. Genetic analysis identified a GI.3P-GI.2 variant supporting a non-Australian origin for the incursion; however, more accurate identification of the source of the incursion remains challenging due to the wide global distribution of the GI.3P-GI.2 variant. Furthermore, our analysis suggests the spread of the virus between the North and South Islands of New Zealand at least twice, dated to mid-2017 and around 2018. Further phylogenetic analysis also revealed a strong phylogeographic pattern. So far, no recombination events with endemic benign New Zealand rabbit caliciviruses have been identified. This study highlights the need for further research and surveillance to monitor the distribution and diversity of lagoviruses in New Zealand and to detect incursions of novel variants.

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.

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.

First report of GI.1aP-GI.2 recombinants of rabbit hemorrhagic disease virus in domestic rabbits in China

The rabbit hemorrhagic disease virus 2 (RHDV2 or GI.2) is a highly contagious agent leading to lethal disease in rabbits. It frequently recombines with other Lagovirus genus, generating epidemical variants with high pathogenicity. In this study, twenty-two liver samples tested positive for GI.2 VP60 gene, were collected in rabbit farms from several geographical regions in China. All GI.2 positive specimens were submitted for RT-PCR detection, nucleotide sequencing and phylogenetic analysis. In addition, suspected GI.2 recombinants were evaluated for virus virulence. The results showed that nine presumptive recombinants were identified by testing for RdRp-VP60 recombination. In these recombinants, four were selected to fully characterize the genome of novel GI.2 recombinant variants, which were described as GI.1aP-GI.2. The nucleotide sequence of these novel variants showed unique recombination pattern and phylogenetic features compared to currently prevalent GI.2 variants. Furthermore, this distinctive recombination of new variant SCNJ-2021 moderately enhanced the virulence of GI.2, even for rabbits vaccinated against parental GI.2. In conclusion, the novel GI.1aP-GI.2 recombinants were identified in rabbit industry in China for the first time, which expanded the knowledge on the phylodynamics and genomic diversity of GI.2 genotype. The rapid molecular evolution and varied pathogenicity of these virus recombinants highlight the urgent need for epidemiological surveillance and for future prevention of these neglected GI.2 variants.

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.

Structural Basis for Rabbit Hemorrhagic Disease Virus Antibody Specificity

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.

Frequent intergenotypic recombination between the non-structural and structural genes is a major driver of epidemiological fitness in caliciviruses

The diversity of lagoviruses (Caliciviridae) in Australia has increased considerably in recent years. By the end of 2017, five variants from three viral genotypes were present in populations of Australian rabbits, while prior to 2014 only two variants were known. To understand the evolutionary interactions among these lagovirus variants, we monitored their geographical distribution and relative incidence over time in a continental-scale competition study. Within 3 years of the incursion of rabbit haemorrhagic disease virus 2 (RHDV2, denoted genotype GI.1bP-GI.2 [polymerase genotype]P-[capsid genotype]) into Australia, two novel recombinant lagovirus variants emerged: RHDV2-4e (genotype GI.4eP-GI.2) in New South Wales and RHDV2-4c (genotype GI.4cP-GI.2) in Victoria. Although both novel recombinants contain non-structural genes related to those from benign, rabbit-specific, enterotropic viruses, these variants were recovered from the livers of both rabbits and hares that had died acutely. This suggests that the determinants of host and tissue tropism for lagoviruses are associated with the structural genes, and that tropism is intricately connected with pathogenicity. Phylogenetic analyses demonstrated that the RHDV2-4c recombinant emerged independently on multiple occasions, with five distinct lineages observed. Both the new RHDV2-4e and -4c recombinant variants replaced the previous dominant parental RHDV2 (genotype GI.1bP-GI.2) in their respective geographical areas, despite sharing an identical or near-identical (i.e. single amino acid change) VP60 major capsid protein with the parental virus. This suggests that the observed replacement by these recombinants was not driven by antigenic variation in VP60, implicating the non-structural genes as key drivers of epidemiological fitness. Molecular clock estimates place the RHDV2-4e recombination event in early to mid-2015, while the five RHDV2-4c recombination events occurred from late 2015 through to early 2017. The emergence of at least six viable recombinant variants within a 2-year period highlights the high frequency of these events, detectable only through intensive surveillance, and demonstrates the importance of recombination in lagovirus evolution.

First detection of rabbit haemorrhagic disease virus (RHDV2) in Singapore

Rabbit haemorrhagic disease (RHD) is a significant viral disease caused by infection with Rabbit haemorrhagic disease virus (RHDV). The first documented cases of RHDV in Singapore occurred in adult pet European rabbits (Oryctolagus cuniculus) in September 2020. Rabbits presented with acute hyporexia, lethargy, huddled posture, and varying degrees of pyrexia and tachypnoea. Clinical pathology consistently reflected markedly elevated alanine aminotransferase (ALT), aspartate aminotransferase (AST) and alkaline phosphatase (ALKP). Hepatic lobe torsion was ruled out using ultrasonography and colour Doppler studies in all patients. A total of 11 rabbits owned by 3 families were presented to the clinics; 8/11 rabbits died within 48 hr of presentation, while the remaining two rabbits had recovered after prolonged hospitalization and one rabbit was aclinical. Histopathology revealed acute, marked diffuse hepatocellular necrosis and degeneration, findings which were suggestive for RHDV infection and prompted the undertaking of further molecular diagnostics. Subsequent polymerase chain reaction of the liver samples detected RHDV RNA. Molecular characterization of viral genomes by whole genome sequencing revealed that the outbreak strain was of the genotype GI.2 (RHDV2/RHDVb). Nucleotide sequences of the VP60 gene were compared with various RHDV variants using phylogenetic analysis. The sample genome shared highest sequence identity with a GI.2-genotyped virus from GenBank (RHDV isolate Algarve 1 polyprotein and minor structural protein (VP10) genes, GenBank accession KF442961). The combination of clinical, histopathological, molecular and sequencing technologies enabled rapid detection and detailed genetic characterization of the RHDV virus causing the present outbreak for prompt implementation of disease control measures in Singapore. Further epidemiological investigations of potential virus introduction into Singapore are ongoing.

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.

Genomic insights into a population of introduced European rabbits Oryctolagus cuniculus in Australia and the development of genetic resistance to rabbit hemorrhagic disease virus

The European rabbit (Oryctolagus cuniculus) is one of the most devastating invasive species in Australia. Since the 1950s, myxoma virus (MYXV) and rabbit haemorrhagic disease virus (RHDV) have been used to manage overabundant rabbit populations. Resistance to MYXV was observed within a few years of the release. More recently, resistance to lethal RHDV infection has also been reported, undermining the efficiency of landscape-scale rabbit control. Previous studies suggest that genetic resistance to lethal RHDV infection may differ locally between populations, yet the mechanisms of genetic resistance remain poorly understood. Here, we used genotyping by sequencing (GBS) data representing a reduced representation of the genome, to investigate Australian rabbit populations. Our aims were to understand the relationship between populations and identify possible genomic signatures of selection for RHDV resistance. One population we investigated had previously been reported to show levels of resistance to lethal RHDV infection. This population was compared to three other populations with lower or no previously reported RHDV resistance. We identified a set of novel candidate genes that could be involved in host-pathogen interactions such as virus binding and infection processes. These genes did not overlap with previous studies on RHDV resistance carried out in different rabbit populations, suggesting that multiple mechanisms are feasible. These findings provide useful insights into the different potential mechanisms of genetic resistance to RHDV virus which will inform future functional studies in this area.

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.

Rabbit Hemorrhagic Disease Virus Isolated from Diseased Alpine Musk Deer (Moschus sifanicus)

Rabbit hemorrhagic disease virus (RHDV) is the causative agent of rabbit hemorrhagic disease (RHD), and its infection results in mortality of 70-90% in farmed and wild rabbits. RHDV is thought to replicate strictly in rabbits. However, there are also reports showing that gene segments from the RHDV genome or antibodies against RHDV have been detected in other animals. Here, we report the detection and isolation of a RHDV from diseased Alpine musk deer (Moschussifanicus). The clinical manifestations in those deer were sudden death without clinical signs and hemorrhage in the internal organs. To identify the potential causative agents of the disease, we used sequence independent single primer amplification (SISPA) to detect gene segments from viruses in the tissue samples collected from the dead deer. From the obtained sequences, we identified some gene fragments showing very high nucleotide sequence similarity with RHDV genome. Furthermore, we identified caliciviral particles using an electron microscope in the samples. The new virus was designated as RHDV GS/YZ. We then designed primers based on the genome sequence of an RHDV strain CD/China to amplify and sequence the whole genome of the virus. The genome of the virus was determined to be 7437 nucleotides in length, sharing the highest genome sequence identity of 98.7% with a Chinese rabbit strain HB. The virus was assigned to the G2 genotype of RHDVs according to the phylogenetic analyses based on both the full-length genome and VP60 gene sequences. Animal experiments showed that GS/YZ infection in rabbits resulted in the macroscopic and microscopic lesions similar to that caused by the other RHDVs. This is the first report of RHDV isolated from Alpine musk deer, and our findings extended the epidemiology and host range of RHDV.

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.

Biocontrol of the Common Carp (Cyprinus carpio) in Australia: A Review and Future Directions

Invasive pest species are recognized as one of the important drivers of reduced global biodiversity. In Australia, the 267 invasive plant, animal and microbial species, established since European colonization in the 1770s, have been unequivocally declared the most important threat to species diversity in this country. One invasive pest, the common carp (Cyprinus carpio), has been targeted in an integrated pest management plan that might include cyprinid herpesvirus 3 (CyHV-3) as a potential biocontrol agent. The species-specificity of the released virus (and of field variants that will inevitably arise) has been assessed, and the virus judged to be safe. It has also been hypothesised that, because the virulence of the CyHV-3 will likely decline following release, the virus should be used strategically: initially, the aim would be to markedly reduce numbers of carp in naive populations, and then some other, as yet uncertain, complementary broad-scale control measure would knock-down carp numbers even further. Brief results are included from recent studies on the modelling of release and spread of the virus, the ecological and social concerns associated with virus release, and the restoration benefits that might be expected following carp control. We conclude that, while further work is required (on the virus, the target species, environmental issues, and especially the identification of a suitable broad-scale complementary control measure), optimism must prevail in order to ensure an eventual solution to this important environmental problem.

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.

Bioinformatics analysis of capsid protein of different subtypes rabbit hemorrhagic disease virus

Background

Rabbit Hemorrhagic Disease Virus (RHDV) belongs to the Caliciviridae family, is a highly lethal pathogen to rabbits. Increasing numbers of studies have demonstrated the existence of antigenic variation in RHDV, leading to the emergence of a new RHDV isolate (RHDVb). However, the underlying factors determining the emergence of the new RHDV and its unpredictable epidemiology remain unclear. To investigate these issues, we selected more than 184 partial and/or complete genome sequences of RHDV from GenBank and analyzed their phylogenetic relationships, divergence, and predicted protein modification sites.

Results

Phylogenetic analysis showed that classic RHDV isolates, RHDVa, and RHDVb formed different clades. It’s interesting to note that RHDVa being more closely related to classic RHDV than RHDVb, while RHDVb had a closer genetic relationship to Rabbit Calicivirus (RCV) than to classic RHDV isolates. Moreover, divergence analysis suggested that the accumulation of amino acid (aa) changes might be a consequence of adaptive diversification of capsid protein (VP60) during the division between classical RHDV, RHDVa, RHDVb, and RCV. Notably, the prediction of N-glycosylation sites suggested that RHDVb subtypes had two unique N-glycosylation sites (aa 301, 362) but lacked three other N-glycosylation sites (aa 45, 308, 474) displayed in classic RHDV and RHDVa VP60 implying this divergence of N-glycosylation sites in RHDV might affect viral virulence. Analysis of phosphorylation sites also indicated that some phosphorylation sites in RHDVa and RHDVb differed from those in classic RHDV, potentially related to antigenic variation in RHDV.

Conclusion

The genetic relationship between RHDVb and RCV was closer than classic RHDV isolates. Moreover, compared to RHDV and RHDVa, RHDVb had two unique N-glycosylation sites but lacked three sites, which might affect the virulence of RHDV. These results may provide new clues for further investigations of the origin of new types of RHDV and the mechanisms of genetic variation in RHDV.

Calicivirus RNA-Dependent RNA Polymerases: Evolution, Structure, Protein Dynamics, and Function

The Caliciviridae are viruses with a positive-sense, single-stranded RNA genome that is packaged into an icosahedral, environmentally stable protein capsid. The family contains five genera (Norovirus, Nebovirus, Sapovirus, Lagovirus, and Vesivirus) that infect vertebrates including amphibians, reptiles, birds, and mammals. The RNA-dependent RNA polymerase (RdRp) replicates the genome of RNA viruses and can speed up evolution due to its error-prone nature. Studying calicivirus RdRps in the context of genuine virus replication is often hampered by a lack of suitable model systems. Enteric caliciviruses and RHDV in particular are notoriously difficult to propagate in cell culture; therefore, molecular studies of replication mechanisms are challenging. Nevertheless, research on recombinant proteins has revealed several unexpected characteristics of calicivirus RdRps. For example, the RdRps of RHDV and related lagoviruses possess the ability to expose a hydrophobic motif, to rearrange Golgi membranes, and to copy RNA at unusually high temperatures. This review is focused on the structural dynamics, biochemical properties, kinetics, and putative interaction partners of these RdRps. In addition, we discuss the possible existence of a conserved but as yet undescribed structural element that is shared amongst the RdRps of all caliciviruses.

The discovery of three new hare lagoviruses reveals unexplored viral diversity in this genus

Our knowledge of mammalian viruses has been strongly skewed toward those that cause disease in humans and animals. However, recent metagenomic studies indicate that most apparently healthy organisms carry viruses, and that these seemingly benign viruses may comprise the bulk of virus diversity. The bias toward studying viruses associated with overt disease is apparent in the lagoviruses (family Caliciviridae) that infect rabbits and hares: although most attention has been directed toward the highly pathogenic members of this genus—rabbit haemorrhagic disease virus and European brown hare syndrome virus—a number of benign lagoviruses have also been identified. To determine whether wild European brown hares in Australia might also carry undetected benign viruses, we used a meta-transcriptomics approach to explore the gut and liver RNA viromes of these invasive animals. This led to the discovery of three new lagoviruses. While one was only detected in a single hare, the other two viruses were detected in 20 per cent of all animals tested. All three viruses were most closely related to other hare lagoviruses, but were phylogenetically distinct from both known viruses and from each other, indicating that lagoviruses have circulated for longer than previously assumed. Their evolution was also characterised by complex recombination events. Mapping mutations onto the lagovirus phylogeny revealed no amino acid changes that were consistently associated with virulence phenotype. Overall, our study points to extensive unsampled diversity in this genus, such that additional metagenomic studies are needed to fill gaps in the lagovirus phylogeny and better understand the evolutionary history of this important group of mammalian viruses.

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.

Detection and Circulation of a Novel Rabbit Hemorrhagic Disease Virus in Australia

The highly virulent rabbit hemorrhagic disease virus (RHDV) has been widely used in Australia and New Zealand since the mid-1990s to control wild rabbits, an invasive vertebrate pest in these countries. In January 2014, an exotic RHDV was detected in Australia, and 8 additional outbreaks were reported in both domestic and wild rabbits in the 15 months following its detection. Full-length genomic analysis revealed that this virus is a recombinant containing an RHDVa capsid gene and nonstructural genes most closely related to nonpathogenic rabbit caliciviruses. Nationwide monitoring efforts need to be expanded to assess if the increasing number of different RHDV variants circulating in the Australian environment will affect biological control of rabbits. At the same time, updated vaccines and vaccination protocols are urgently needed to protect pet and farmed rabbits from these novel rabbit caliciviruses.

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.

Rabbit Hemorrhagic Disease Virus 2 (RHDV2; GI.2) Is Replacing Endemic Strains of RHDV in the Australian Landscape within 18 Months of Its Arrival

Rabbit hemorrhagic disease virus 2 (RHDV2; Lagovirus GI.2) is a pathogenic calicivirus that affects European rabbits (Oryctolagus cuniculus) and various hare (Lepus) species. GI.2 was first detected in France in 2010 and subsequently caused epidemics in wild and domestic lagomorph populations throughout Europe. In May 2015, GI.2 was detected in Australia. Within 18 months of its initial detection, GI.2 had spread to all Australian states and territories and rapidly became the dominant circulating strain, replacing Rabbit hemorrhagic disease virus (RHDV/GI.1) in mainland Australia. Reconstruction of the evolutionary history of 127 Australian GI.2 isolates revealed that the virus arrived in Australia at least several months before its initial description and likely circulated unnoticed in wild rabbit populations in the east of the continent prior to its detection. GI.2 sequences isolated from five hares clustered with sequences from sympatric rabbit populations sampled contemporaneously, indicating multiple spillover events into hares rather than an adaptation of the Australian GI.2 to a new host. Since the presence of GI.2 in Australia may have wide-ranging consequences for rabbit biocontrol, particularly with the release of the novel biocontrol agent GI.1a/RHDVa-K5 in March 2017, ongoing surveillance is critical to understanding the interactions of the various lagoviruses in Australia and their impact on host populations.IMPORTANCE This study describes the spread and distribution of Rabbit hemorrhagic disease virus 2 (GI.2) in Australia since its first detection in May 2015. Within the first 18 months following its detection, RHDV2 spread from east to west across the continent and became the dominant strain in all mainland states of Australia. This has important implications for pest animal management and for owners of pet and farmed rabbits, as there currently is no effective vaccine available in Australia for GI.2. The closely related RHDV (GI.1) is used to control overabundant wild rabbits, a serious environmental and agricultural pest in this country, and it is currently unclear how the widespread circulation of GI.2 will impact ongoing targeted wild rabbit management operations.

A strain-specific multiplex RT-PCR for Australian rabbit haemorrhagic disease viruses uncovers a new recombinant virus variant in rabbits and hares

Rabbit haemorrhagic disease virus (RHDV, or GI.1) is a calicivirus in the genus Lagovirus that has been widely utilized in Australia as a biological control agent for the management of overabundant wild European rabbit (Oryctolagus cuniculus) populations since 1996. Recently, two exotic incursions of pathogenic lagoviruses have been reported in Australia; GI.1a-Aus, previously called RHDVa-Aus, is a GI.1a virus detected in January 2014, and the novel lagovirus GI.2 (previously known as RHDV2). Furthermore, an additional GI.1a strain, GI.1a-K5 (also known as 08Q712), was released nationwide in March 2017 as a supplementary tool for wild rabbit management. To discriminate between these lagoviruses, a highly sensitive strain-specific multiplex RT-PCR assay was developed, which allows fast, cost-effective and sensitive detection of the four pathogenic lagoviruses currently known to be circulating in Australia. In addition, we developed a universal RT-qPCR assay to be used in conjunction with the multiplex assay that broadly detects all four viruses and facilitates quantification of viral RNA load in samples. These assays enable rapid detection, identification and quantification of pathogenic lagoviruses in the Australian context. Using these assays, a novel recombinant lagovirus was detected in rabbit tissue samples, which contained the non-structural genes of GI.1a-Aus and the structural genes of GI.2. This variant was also recovered from the liver of a European brown hare (Lepus europaeus). The impact of this novel recombinant on Australian wild lagomorph populations and its competitiveness in relation to circulating field strains, particularly GI.2, requires further studies.

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.

Benign Rabbit Calicivirus in New Zealand

The Czech v351 strain of rabbit hemorrhagic disease virus (RHDV1) is used in Australia and New Zealand as a biological control agent for rabbits, which are important and damaging introduced vertebrate pests in these countries. However, nonpathogenic rabbit caliciviruses (RCVs) can provide partial immunological cross-protection against lethal RHDV infection and thus interfere with effective rabbit biocontrol. Antibodies that cross-reacted against RHDV antigens were found in wild rabbits before the release of RHDV1 in New Zealand in 1997, suggesting that nonpathogenic RCVs were already present in New Zealand. The aim of this study was to confirm the presence of nonpathogenic RCV in New Zealand and describe its geographical distribution. RCV and RHDV antibody assays were used to screen serum samples from 350 wild rabbits from 14 locations in New Zealand. The serological survey indicated that both RCV and RHDV are widespread in New Zealand wild rabbits, with antibodies detected in 10 out of 14 and 12 out of 14 populations, respectively. Two closely related RCV strains were identified in the duodenal tissue from a New Zealand wild rabbit (RCV Gore-425A and RCV Gore-425B). Both variants are most closely related to Australian RCV strains, but with 88% nucleotide identity, they are genetically distinct. Phylogenetic analysis revealed that the New Zealand RCV strains fall within the genetic diversity of the Australian RCV isolates, indicating a relatively recent movement of RCVs between Australia and New Zealand.IMPORTANCE Wild rabbits are important and damaging introduced vertebrate pests in Australia and New Zealand. Although RHDV1 is used as a biological control agent, some nonpathogenic RCVs can provide partial immunological cross-protection against lethal RHDV infection and thus interfere with its effectiveness for rabbit control. The presence of nonpathogenic RCVs in New Zealand wild rabbits has been long hypothesized, but earlier attempts to isolate a New Zealand RCV strain have been unsuccessful. Therefore, it is important to determine if such nonpathogenic viruses exist in New Zealand rabbits, especially considering the proposed introduction of new RHDV strains into New Zealand as biocontrols.

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.

RNA-Dependent RNA Polymerases of Both Virulent and Benign Rabbit Caliciviruses Induce Striking Rearrangement of Golgi Membranes

The extremely pathogenic Rabbit haemorrhagic disease virus (RHDV) and the completely benign Rabbit calicivirus (RCV) are closely related members of the genus Lagovirus (family Caliciviridae). The molecular mechanisms that determine the dramatic difference in virulence are unknown, but indirect evidence suggests that different properties of their RNA-dependent RNA polymerases (RdRps) may at least partially be responsible for the contrasting phenotypes. Here we report that the unusual ability of the RHDV RdRp to induce a striking rearrangement of the Golgi network is not specific to RHDV, but a common feature of virulent and benign rabbit caliciviruses alike. Expression of rabbit calicivirus RdRps induced a redistribution of both cis/medial and medial/trans Golgi membrane markers, but not that of an endoplasmic reticulum membrane marker. Inactivating mutations in the conserved GDD motif did not abolish the ability of RHDV RdRp to rearrange the Golgi network, suggesting that polymerase activity and metal co-factors are not required for this function. Finally, we discuss possible implications of RdRp-induced membrane rearrangements on virus replication and host immune responses.

Rabbit haemorrhagic disease: Macquarie Island rabbit eradication adds to knowledge on both pest control and epidemiology

Rabbit haemorrhagic disease virus (RHDV), introduced into in Australia and New Zealand as a biological-control agent for wild rabbits, is least efficacious in cool humid areas where a non-pathogenic calicivirus (RCV-A1) also circulates. Heavy rabbit mortality following release of RHDV on cold sub-Antarctic Macquarie Island, where RCV-A1 was apparently absent, not only complemented the planned rabbit eradication operations, especially by reducing secondary poisoning of sea-birds from aerial baiting, but also ruled out cool or humid climate as a major limiting factor of disease spread. In turn, this has advanced the idea that RCV-A1 antibodies inhibit RHDV spread as well as reducing disease severity and mortality.