Molecular epidemiology of rabbit haemorrhagic disease virus outbreaks in France during 1988 to 1995

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.

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.

A Review on the Methods Used for the Detection and Diagnosis of Rabbit Hemorrhagic Disease Virus (RHDV)

Since the early 1980s, the European rabbit (Oryctolagus cuniculus) has been threatened by the rabbit hemorrhagic disease (RHD). The disease is caused by a lagovirus of the family Caliciviridae, the rabbit hemorrhagic disease virus (RHDV). The need for detection, identification and further characterization of RHDV led to the development of several diagnostic tests. Owing to the lack of an appropriate cell culture system for in vitro propagation of the virus, much of the methods involved in these tests contributed to our current knowledge on RHD and RHDV and to the development of vaccines to contain the disease. Here, we provide a comprehensive review of the RHDV diagnostic tests used since the first RHD outbreak and that include molecular, histological and serological techniques, ranging from simpler tests initially used, such as the hemagglutination test, to the more recent and sophisticated high-throughput sequencing, along with an overview of their potential and their limitations.

An outbreak of rabbit hemorrhagic disease (RHD) caused by Lagovirus europaeus GI.2/rabbit hemorrhagic disease virus 2 (RHDV2) in Ehime, Japan

A total of ten 1–2-year-old rabbits died within 2 weeks at a facility in Ehime prefecture in May 2019. Necropsy revealed liver discoloration and fragility, hemorrhage of some organs and blood coagulation failure. On histopathologic examination, necrotizing hepatitis was a common finding, together with fibrin thrombi in the small vessels and hemorrhage in some organs. Rabbit hemorrhagic disease (RHD) virus gene was detected in liver samples, and viral particles of approximately 32 nm in diameter were found in the cytoplasm of degenerated hepatocytes by electron microscopy. Phylogenetic analysis based on the partial VP60 gene sequence classified it as Lagovirus europaeus GI.2/RHDV2. This is the first confirmed outbreak of RHD caused by globally emerging GI.2/RHDV2 in Japan.

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.

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.

Detection of rabbit Haemorrhagic disease virus 2 during the wild rabbit (Oryctolagus cuniculus) eradication from the Berlengas archipelago, Portugal

BACKGROUND

In the regular wildlife monitoring action carried out in the summer of the past few years at the Berlenga Island, wild rabbits (Oryctolagus cuniculus) have been repeatedly found dead. However, the origin of those deaths was never investigated. Our aim was to investigate the cause of death of 11 rabbits collected between April and May 2016.

RESULTS

While screening samples from rabbit carcasses for the major viral rabbit pathogens, five tested positive to RHDV2 but all were negative for RHDV and myxoma virus (MYXV). For six RHDV2-negative specimens, emaciation and parasitism were considered the most probable cause of death. Lesions identified in the RHDV2-positive rabbits included non-suppurative diffuse hepatic necrosis and pulmonary lesions varying from congestion and oedema of the lungs to interstitial pneumonia. Sequencing analysis of the vp60 gene obtained from two specimens showed identical vp60 sequences. Comparison with other known RHDV2 strains from public databases through BLAST analysis revealed a closer similarity with strains from Alentejo collected during 2013. Maximum Likelihood and Bayesian phylogenetic analysis showed that the 2016 strains from the archipelago have a higher resemblance with a group of strains mostly collected in the South of Portugal between 2013 and 2014.

CONCLUSION

The results suggest that RHDV2 may have been introduced on the Berlenga Island a few years ago, having evolved separately from mainland strains due to insularity.

Adaptive diversification between the classic rabbit hemorrhagic disease virus (RHDV) and the RHDVa isolates: A genome-wide perspective

Rabbit haemorrhagic disease virus (RHDV) is a highly infectious pathogen that causes high mortality in wild and domestic rabbits. RHDV could be divided into two subtypes, classic RHDV and RHDVa, which present clear genetic, antigenic, and epidemiological differences. To further understand the nature of the diversity, we performed a genome-wide evolutionary study on the classic RHDV and RHDVa isolates. The results show that RHDV had experienced adaptive diversification with the dividing process of these subtypes. Furthermore, amino acid changes relevant to the adaptive diversification mainly cluster in viral capsid protein VP60. These results might be beneficial for a further understanding the function of VP60 and provide helpful hints for the genetic basis of RHDV emergence and re-emergence.

Emergence of rabbit haemorrhagic disease virus 2 in the archipelago of Madeira, Portugal (2016-2017)

We report the detection of rabbit haemorrhagic disease virus 2 (RHDV2) in the Madeira archipelago, Portugal. Viral circulation was confirmed by RT-qPCR and vp60 sequencing. Epidemiological data revealed the outbreak initiated in October 2016 in Porto Santo affecting wild and domestic rabbits. It was then detected three months later on the island of Madeira. Five haplotypes were identified and a genetic overall similarity of 99.54 to 99.89% was observed between the two viral populations. Unique single nucleotide polymorphisms were recognised in the Madeira archipelago strains, two of which resulting in amino acid substitutions at positions 480 and 570 in the VP60 protein. Phylogenetic investigation by Maximum Likelihood showed all the vp60 sequences from the Madeira archipelago group together with high bootstraps. The analysis also showed that the Madeira archipelago strains are closely related to the strains detected in the south of mainland Portugal in 2016, suggesting a possible introduction from the mainland. The epidemiological data and high genetic similarity indicate a common source for the Porto Santo and Madeira RHDV2 outbreaks. Human activity related to hunting was most probably at the origin of the Madeira outbreak.

Progression of rabbit haemorrhagic disease virus 2 upon vaccination in an industrial rabbitry: a laboratorial approach

<p>Rabbit haemorrhagic disease virus 2 (RHDV2) emerged recently in several European countries, leading to extensive economic losses in the industry. In response to this new infection, specific inactivated vaccines were developed in Europe and full and rapid setup of protective immunity induced by vaccination was reported. However, data on the efficacy of these vaccines in an ongoing-infection scenario is unavailable. In this study we investigated an infected RHDV2 indoor industrial meat rabbitry, where fatalities continued to occur after the implementation of the RHDV2 vaccination, introduced to control the disease. The aim of this study was to understand if these mortalities were RHDV2-related, to discover if the dead animals showed any common features such as age or time distance from vaccination, and to identify the source of the outbreak. Anatomo-pathological analysis of vaccinated animals with the virus showed lesions compatible with systemic haemorrhagic disease and RHDV2-RNA was detected in 85.7% of the animals tested. Sequencing of the <em>vp60</em> gene amplified from liver samples led to the recognition of RHDV2 field strains demonstrating that after the implementation of vaccination, RHDV2 continued to circulate in the premises and to cause sporadic deaths. A nearby, semi-intensive, RHDV2 infected farm belonging to the same owner was identified as the most probable source of the virus. The main risk factors for virus introduction in these two industries were identified. Despite the virus being able to infect a few of the vaccinated rabbits, the significant decrease in mortality rate observed in vaccinated adult rabbits clearly reflects the efficacy of the vaccination. Nonetheless, the time taken to control the infection also highlights the importance of RHDV2 vaccination prior to the first contact with the virus, highly recommendable in endemic areas, to mitigate the infection’s impact on the industry.</p>

Emergence of new virulent rabbit hemorrhagic disease virus strains in Saudi Arabia

Rabbit hemorrhagic disease is an acute fatal highly contagious viral infectious disease that causes high losses among rabbitries. The disease was first reported in China in 1984 and later on in Saudi Arabia in 1996. The aim of this study was to investigate the emergence and pathogenicity of new rabbit hemorrhagic disease virus (RHDV) strains in Saudi Arabia. The pathogenicity was confirmed by inoculation in susceptible rabbits. Three RHDV strains were detected by reverse transcriptase polymerase chain reaction (RT-PCR) using primers targeting VP60 capsid protein gene in infected rabbitries during 2012 and 2013. These strains clustered into two genetically distinct genogroups related to year of isolation (G2 and G3). All new Saudi Arabia viruses clustered with the European strains, while the old strains clustered with strains from China and America. Based on amino acids and nucleotide sequences, the Saudi Arabia strains (RHD/1/SA/2012, RHD/2/SA/2012, and RHD/3/SA /2013) had high identity with Mexico89, Ca11-ITA, and 00-13,FRA virus; on the other hand, there was a relatively high identity with Bahrain strain. The evolutionary relationship of Saudi RHDVs strains revealed significant nucleotides and amino acid substitutions in hypervariable region E, suggesting the emergence of new RHDVs circulating in Saudi Arabia rabbitries. These antigenic changes represented by the antigenic index might be a potential cause of vaccination failure and raises the need to review the vaccination strategies against RHD.

Identification of novel rabbit hemorrhagic disease virus B-cell epitopes and their interaction with host histo-blood group antigens

Rabbit haemorrhagic disease, caused by rabbit hemorrhagic disease virus (RHDV), results in the death of millions of adult rabbits worldwide, with a mortality rate that exceeds 90%. The sole capsid protein, VP60, is divided into shell (S) and protruding (P) domains, and the more exposed P domain likely contains determinants for cell attachment and antigenic diversity. Nine mAbs against VP60 were screened and identified. To map antigenic epitopes, a set of partially overlapping and consecutive truncated proteins spanning VP60 were expressed. The minimal determinants of the linear B-cell epitopes of VP60 in the P domain, N(326)PISQV(331), D(338)MSFV(342) and K(562)STLVFNL(569), were recognized by one (5H3), four (1B8, 3D11, 4C2 and 4G2) and four mAbs (1D4, 3F7, 5G2 and 6B2), respectively. Sequence alignment showed epitope D(338)MSFV(342) was conserved among all RHDV isolates. Epitopes N(326)PISQV(331) and K(562)STLVFNL(569) were highly conserved among RHDV G1-G6 and variable in RHDV2 strains. Previous studies demonstrated that native viral particles and virus-like particles (VLPs) of RHDV specifically bound to synthetic blood group H type 2 oligosaccharides. We established an oligosaccharide-based assay to analyse the binding of VP60 and epitopes to histo-blood group antigens (HBGAs). Results showed VP60 and its epitopes (aa 326-331 and 338-342) in the P2 subdomain could significantly bind to blood group H type 2. Furthermore, mAbs 1B8 and 5H3 could block RHDV VLP binding to synthetic H type 2. Collectively, these two epitopes might play a key role in the antigenic structure of VP60 and interaction of RHDV and HBGA.

Is the new variant RHDV replacing genogroup 1 in Portuguese wild rabbit populations?

The Lagovirus rabbit hemorrhagic disease virus (RHDV), a member of the family Caliciviridae, severely affects European rabbit (Oryctolagus cuniculus) populations by causing rabbit hemorrhagic disease (RHD). RHDV is subdivided in six genogroups but, more recently, a new RHDV variant with a unique genetic and antigenic profile emerged. We performed a study in rabbits found dead in the field during 2013 and 2014 in Portugal to determine the prevalence of this new variant versus the classical RHDV. Fifty-seven liver samples were screened for the presence of RHDV and positive samples were genotyped. All cases of RHDV infection were caused by the new variant. The only former genogroup circulating in Portugal, G1, was not detected. We hence conclude that the new RHDV variant is replacing G1 in Portugal, probably due to a selective advantage. This sudden and rapid replacement emphasizes the necessity of continued monitoring of wild rabbit populations.

Structural analysis of a rabbit hemorrhagic disease virus binding to histo-blood group antigens

Rabbit hemorrhagic disease virus (RHDV) is a member of the Caliciviridae family (Lagovirus genus). RHDV is highly contagious and attaches to epithelial cells in the digestive or respiratory tract, leading to massive lesions with high mortality rates. A new variant of RHDV (termed RHDVb) recently has emerged, and previously vaccinated rabbits appear to have little protection against this new strain. Similar to human norovirus (Caliciviridae, Norovirus genus), RHDV binds histo-blood group antigens (HBGAs), and this is thought to be important for infection. Here, we report the HBGA binding site on the RHDVb capsid-protruding domain (P domain) using X-ray crystallography. The HBGA binding pocket was located in a negatively charged patch on the side of the P domain and at a dimeric interface. Residues from both monomers contributed to the HBGA binding and involved a network of direct hydrogen bonds and water-mediated interactions. An amino acid sequence alignment of different RHDV strains indicated that the residues directly interacting with the ABH-fucose of the HBGAs (Asp472, Asn474, and Ser479) were highly conserved. This result suggested that different RHDV strains also could bind HBGAs at the equivalent pocket. Moreover, several HBGA binding characteristics between RHDVb and human genogroup II norovirus were similar, which indicated a possible convergent evolution of HBGA binding interactions. Further structural studies with other RHDV strains are needed in order to better understand the HBGA binding mechanisms among the diverse RHDV strains.

IMPORTANCE

We identified, for the first time, the HBGA binding site on an RHDVb P domain using X-ray crystallography. Our results showed that RHDVb and human genogroup II noroviruses had similar HBGA binding interactions. Recently, it was discovered that synthetic HBGAs or HBGA-expressing enteric bacteria could enhance human genogroup II norovirus infection in B cells. Considering that RHDVb and genogroup II norovirus similarly interacted with HBGAs, it may be possible that a comparable cell culture system also could work with RHDVb. Taken together, these new findings will extend our understanding of calicivirus HBGA interactions and may help to elucidate the specific roles of HBGAs in calicivirus infections.

Molecular evolution and antigenic variation of European brown hare syndrome virus (EBHSV)

European brown hare syndrome virus (EBHSV) is the aetiological agent of European brown hare syndrome (EBHS), a disease affecting Lepus europaeus and Lepus timidus first diagnosed in Sweden in 1980. To characterize EBHSV evolution we studied hare samples collected in Sweden between 1982 and 2008. Our molecular clock dating is compatible with EBHSV emergence in the 1970s. Phylogenetic analysis revealed two lineages: Group A persisted until 1989 when it apparently suffered extinction; Group B emerged in the mid-1980s and contains the most recent strains. Antigenic differences exist between groups, with loss of reactivity of some MAbs over time, which are associated with amino acid substitutions in recognized epitopes. A role for immune selection is also supported by the presence of positively selected codons in exposed regions of the capsid. Hence, EBHSV evolution is characterized by replacement of Group A by Group B viruses, suggesting that the latter possess a selective advantage.

Molecular epidemiology of Rabbit Haemorrhagic Disease Virus in Australia: when one became many

Rabbit Haemorrhagic Disease Virus (RHDV) was introduced into Australia in 1995 as a biological control agent against the wild European rabbit (Oryctolagus cuniculus). We evaluated its evolution over a 16-year period (1995-2011) by examining 50 isolates collected throughout Australia, as well as the original inoculum strains. Phylogenetic analysis of capsid protein VP60 sequences of the Australian isolates, compared with those sampled globally, revealed that they form a monophyletic group with the inoculum strains (CAPM V-351 and RHDV351INOC). Strikingly, despite more than 3000 rereleases of RHDV351INOC since 1995, only a single viral lineage has sustained its transmission in the long-term, indicative of a major competitive advantage. In addition, we find evidence for widespread viral gene flow, in which multiple lineages entered individual geographic locations, resulting in a marked turnover of viral lineages with time, as well as a continual increase in viral genetic diversity. The rate of RHDV evolution recorded in Australia -4.0 (3.3-4.7) × 10(-3) nucleotide substitutions per site per year - was higher than previously observed in RHDV, and evidence for adaptive evolution was obtained at two VP60 residues. Finally, more intensive study of a single rabbit population (Turretfield) in South Australia provided no evidence for viral persistence between outbreaks, with genetic diversity instead generated by continual strain importation.

Emergence of a new lagovirus related to Rabbit Haemorrhagic Disease Virus

Since summer 2010, numerous cases of Rabbit Haemorrhagic Disease (RHD) have been reported in north-western France both in rabbitries, affecting RHD-vaccinated rabbits, and in wild populations. We demonstrate that the aetiological agent was a lagovirus phylogenetically distinct from other lagoviruses and which presents a unique antigenic profile. Experimental results show that the disease differs from RHD in terms of disease duration, mortality rates, higher occurrence of subacute/chronic forms and that partial cross-protection occurs between RHDV and the new RHDV variant, designated RHDV2. These data support the hypothesis that RHDV2 is a new member of the Lagovirus genus. A molecular epidemiology study detected RHDV2 in France a few months before the first recorded cases and revealed that one year after its discovery it had spread throughout the country and had almost replaced RHDV strains. RHDV2 was detected in continental Italy in June 2011, then four months later in Sardinia.

Detection of RHDVa on the Iberian Peninsula: isolation of an RHDVa strain from a Spanish rabbitry

Rabbit haemorrhagic disease virus (RHDV), genus Lagovirus, family Caliciviridae, causes a large number of deaths in wild and domestic adult European rabbits (Oryctolagus cuniculus). The first documented outbreak dates from 1984 in China, but the virus rapidly dispersed worldwide. In 1997, an antigenic variant was detected in Italy and designated RHDVa. Despite causing symptoms similar to those caused by classic RHDV strains, marked antigenic and genetic differences exist. In some parts of Europe, RHDVa is replacing classic strains. Here, we report the presence of RHDVa on the Iberian Peninsula, where this variant was thought not to contribute to viral diversity.

Analysis of genetic variability and phylogenetic analysis of selected Czech and French strains of rabbit haemorrhagic disease virus (RHDV)

The objective of this study was to analyse the genetic variability and phylogenetic analysis of six strains of rabbit haemorrhagic disease virus (RHDV), including four Czech strains (CAMPV-351, CAMPV-561, CAMPV-562, CAMPV-558) and two French strains (Fr-1, Fr-2), on the basis of a fragment of the VP60 capsid structural protein-coding gene N-terminal region. The results of our own studies were compared to 26 RHDV strains obtained from GenBank. The analysis of the genetic variability of six RHDV strains indicated that the CAMPV-561 strain is the most genetically variable. Less variable were the Fr-1 and Fr-2 strains, while the least variable was CAMPV-351. In turn, the genetic distance among the six analysed strains and 26 strains obtained from GenBank was the greatest for CAMPV-351 and Whn/China [11.3 % according to the observed divergence (OD) method and 12.2 % according to the maximum likelihood (ML) method], while it was the lowest for CAMPV-351 and FRG (0.8 % in both the OD and ML methods). In turn, the scale of the genetic distances among the six analysed strains and five RHDVa strains (99-05, NY-01, Whn/China, Triptis, Iowa2000) ranged from 9.3–10.3 % in the OD method to 10.3–13.7 % in the ML method. The image of phylogenetic dependencies generated for the strains analysed and those obtained from GenBank revealed their distribution to be in five genetic groups (G1–G5), whereas the analysed strains were included in genetic groups 2 and 3.

Phylogenetic analysis of rabbit hemorrhagic disease virus in China and the antigenic variation of new strains

This study aimed to investigate rabbit hemorrhagic disease virus (RHDV) in China. VP60 sequences of five RHDVs collected by our team, as well as those of 16 other published Chinese RHDV strains, were analyzed. Polygenic analysis using MEGA 4 software showed that 20 of the 21 Chinese strains could be clustered in the RHDVa subgroup, and WX/China/1984 was different from them. The Chinese RHDV strains were further classified into four subgroups, CH1 to CH4. Subgroup CH1, represented by the WX/China/1984 strain, was not prevalent in China after the first RHDV epidemic strain was reported. The CH2, CH3, and CH4 subgroups were far different from the CH1 subgroup, formed three separate clusters, and were distributed according to the time the strains were collected. Recently collected strains formed a new subgroup (CH4), represented by new RHDV varieties identified by challenging immunized rabbits and by comparison of genomic sequences. The present work is the first comprehensive analysis of Chinese RHDV and reveals a new RHDV variation that should be carefully monitored.

Rabbit haemorrhagic disease (RHD) and rabbit haemorrhagic disease virus (RHDV): a review

Rabbit haemorrhagic disease virus (RHDV) is a calicivirus of the genus Lagovirus that causes rabbit haemorrhagic disease (RHD) in adult European rabbits (Oryctolagus cuniculus). First described in China in 1984, the virus rapidly spread worldwide and is nowadays considered as endemic in several countries. In Australia and New Zealand where rabbits are pests, RHDV was purposely introduced for rabbit biocontrol. Factors that may have precipitated RHD emergence remain unclear, but non-pathogenic strains seem to pre-date the appearance of the pathogenic strains suggesting a key role for the comprehension of the virus origins. All pathogenic strains are classified within one single serotype, but two subtypes are recognised, RHDV and RHDVa. RHD causes high mortality in both domestic and wild adult animals, with individuals succumbing between 48-72 h post-infection. No other species has been reported to be fatally susceptible to RHD. The disease is characterised by acute necrotising hepatitis, but haemorrhages may also be found in other organs, in particular the lungs, heart, and kidneys due to disseminated intravascular coagulation. Resistance to the disease might be explained in part by genetically determined absence or weak expression of attachment factors, but humoral immunity is also important. Disease control in rabbitries relies mainly on vaccination and biosecurity measures. Such measures are difficult to be implemented in wild populations. More recent research has indicated that RHDV might be used as a molecular tool for therapeutic applications. Although the study of RHDV and RHD has been hampered by the lack of an appropriate cell culture system for the virus, several aspects of the replication, epizootology, epidemiology and evolution have been disclosed. This review provides a broad coverage and description of the current knowledge on the disease and the virus.

Evolutionary history and molecular epidemiology of rabbit haemorrhagic disease virus in the Iberian Peninsula and Western Europe

BACKGROUND

Rabbit haemorrhagic disease virus (RHDV) is a highly virulent calicivirus, first described in domestic rabbits in China in 1984. RHDV appears to be a mutant form of a benign virus that existed in Europe long before the first outbreak. In the Iberian Peninsula, the first epidemic in 1988 severely reduced the populations of autochthonous European wild rabbit. To examine the evolutionary history of RHDV in the Iberian Peninsula, we collected virus samples from wild rabbits and sequenced a fragment of the capsid protein gene VP60. These data together with available sequences from other Western European countries, were analyzed following Bayesian Markov chain Monte Carlo methods to infer their phylogenetic relationships, evolutionary rates and demographic history.

RESULTS

Evolutionary relationships of RHDV revealed three main lineages with significant phylogeographic structure. All lineages seem to have emerged at a common period of time, between ~1875 and ~1976. The Iberian Peninsula showed evidences of genetic isolation, probably due to geographic barriers to gene flow, and was also the region with the youngest MRCA.Overall, demographic analyses showed an initial increase and stabilization of the relative genetic diversity of RHDV, and a subsequent reduction in genetic diversity after the first epidemic breakout in 1984, which is compatible with a decline in effective population size.

CONCLUSIONS

Results were consistent with the hypothesis that the current Iberian RHDV arose from a single infection between 1869 and 1955 (95% HPD), and rendered a temporal pattern of appearance and extinction of lineages. We propose that the rising positive selection pressure observed throughout the history of RHDV is likely mediated by the host immune system as a consequence of the genetic changes that rendered the virus virulent. Consequently, this relationship is suggested to condition RHDV demographic history.

Identification and partial characterisation of a new Lagovirus in Australian wild rabbits

Rabbit Haemorrhagic Disease Virus (RHDV) is widely used in Australia to control feral rabbit populations. Before RHDV was released on the Australian continent in 1996, antibodies cross-reacting in RHDV specific ELISAs were found in Australian wild rabbits, leading to the hypothesis that a non-pathogenic calicivirus had been circulating in rabbit populations in Australia, potentially providing some level of cross-immunoprotection to RHDV infection. For the detection of this putative virus, a universal lagovirus PCR test was developed to screen a variety of different tissues of wild caught rabbits. We identified a new lagovirus in the intestinal tissues of three apparently healthy young wild rabbits. Quantitative Real Time PCR analysis revealed high concentrations of viral RNA in intestinal tissues and suggests a faecal-oral mode of transmission. Genome organisation and phylogenetic analysis following the sequencing of the entire viral genome revealed a new member of the genus Lagovirus within the family Caliciviridae.

Review of companion animal viral diseases and immunoprophylaxis

In this article we review important established, newly emergent and potential viral diseases of cats, dogs and rabbits. Topics covered include virus epidemiology, disease pathogenesis, existing and prospective immunoprophylaxis against the viruses. For some feline viruses, notably the immunodeficiency virus, leukaemia virus and peritonitis virus, available vaccines are poorly efficacious but there are good prospects for this. A further challenge for the industry is likely to be due to viruses jumping species and the emergence of more virulent variants of established viruses resulting from mutations as has been the case for the canine parvovirus, coronaviruses and feline calicivirus.

Detection of positive selection in the major capsid protein VP60 of the rabbit haemorrhagic disease virus (RHDV)

Mutations were analysed in the major capsid protein VP60 of the rabbit haemorrhagic disease virus (RHDV), a calicivirus responsible for high mortality rates in both wild and domestic European rabbits (Oryctolagus cuniculus). Likelihood of positive selection was estimated using the PAML software applied to 43 non-identical complete sequences of the major capsid protein. Three codons showed signs of positive selection (with posterior probabilities over 95%), one of them is located in the region containing the major antigenic determinants (region E). The presence of positively selected codons (PSCs) in other regions may suggest the existence of other antigenic regions on the major capsid protein that stimulate protective immune responses. At all the 3 PSCs, variation contributes to putative N-glycosylation sites of the protein. An N-glycosylation site is deleted in the non-pathogenic strain RCV. Some of the substitutions at PSCs may alter the polarity and the charge of the protein with possible implications in the protein structure and host interaction. The detection of PSCs should allow a better understanding of the interaction between RHDV and the rabbit immune system.

Molecular and antigenic characterization of rabbit hemorrhagic disease virus isolated in Cuba indicates a distinct antigenic subtype

Phylogenetic analyses conducted on isolates of rabbit hemorrhagic disease virus (RHDV) from throughout the world have shown well-defined genogroups comprising representative strains of the virus and antigenic variants. In this work, we have isolated and characterized RHDV from the major epizootic that occurred in Cuba in 2004-2005. Sequence analysis of the capsid protein gene and antigenic characterization of this strain has allowed its inclusion as a member of the distinct RHDVa subtype. We also found that specific antibodies directed against RHDV reference strains bound to the Cuban isolate in a competition ELISA and inhibited virus hemagglutination in vitro. This is the second report on the molecular characterization of RHDVa circulating in the American region.

Towards a unique and transmissible vaccine against myxomatosis and rabbit haemorrhagic disease for rabbit populations

Currently available vaccines against myxomatosis and rabbit hemorrhagic disease virus (RHDV) are not suited to immunise wild rabbit populations, as vaccines need to be delivered individually by conventional veterinary practices. As an alternative approach, research in Spain has focused on the development of a transmissible vaccine. A recombinant virus has been constructed based on a naturally attenuated myxoma virus (MV) field strain, expressing the RHDV capsid protein (VP60). Following inoculation of rabbits, the recombinant virus (MV-VP60) induced specific antibody responses against MV and RHDV, conferring protection against lethal challenges with both viruses. Furthermore, the recombinant MV-VP60 virus showed a limited horizontal transmission capacity, either by direct contact or in a flea-mediated process, promoting immunisation of contact uninoculated animals. Efficacy and safety of the vaccine have been extensively evaluated under laboratory conditions and in a limited field trial. The development of the transmissible vaccine strategy and the steps being taken to obtain the marketing authorisation for the vaccine in the European Union are presented in this review.

Molecular epidemiology of European brown hare syndrome virus in France between 1989 and 2003

Genetic diversity between French European brown hare syndrome (EBHS) viruses since the disease appeared has been evaluated. Nucleotide sequencing of the partial capsid protein genes of 169 EBHS viruses collected from various parts of France between 1989 and 2003, three reference strains, and a Greek EBHSV collected in 2002 revealed a maximum nucleotide divergence of 11.7%, indicating a high level of conservation between viruses. Two major groups were identified. The first group contained EBHS viruses collected since 1989 from different parts of France, the reference strains, and all of the viruses located in the far north of France. In this group, three genogroups were clearly identified as mainly related to their geographic origin. The distribution of the viruses suggests that the early viruses have not disappeared and have slowly evolved in their area of origin. The second group, supported by a significant bootstrap value, contained the Greek EBHSV with the French EBHS viruses collected between 1999 and 2003 from regions of southern France. It constitutes a newly identified genogroup. Our results demonstrate strong differences in genetic evolution between EBHSV and rabbit haemorrhagic disease virus, with persistence of the earlier EBHS viruses and interaction between the geographical and temporal distributions.

Phylogenetic analysis of rabbit haemorrhagic disease virus (RHDV) strains isolated between 1988 and 2003 in eastern Hungary

To define the genetic variability of RHDV strains collected in eastern Hungary, liver samples from rabbits that had died of RHD were collected between 1988 and 2003. The phylogenetic analysis of a 528-nucleotide-long portion of the gene encoding the VP60 capsid protein assigned the strains into three genogroups. The first group contained viruses from 1988-1993, and a second group comprised isolates from 1994-2002. A third group comprised all of the tested representatives of the RHDVa subtype and a Hungarian isolate from 2003. These findings were supported by the alignments of the deduced amino acid sequences of the VP60 gene and strongly suggest the presence of the RHDVa subtype in Hungary.

Phylogenetic analysis of rabbit haemorrhagic disease virus strains from the Arabian Peninsula: did RHDV emerge simultaneously in Europe and Asia?

Rabbit haemorrhagic disease virus (RHDV) emerged in 1984 in China and subsequently a single strain apparently dispersed worldwide killing millions of rabbits. Two isolates that caused outbreaks in Saudi Arabia and Bahrain have been sequenced and analysed phylogenetically. The Saudi Arabian lineage is directly descended from the Chinese strain, but the Bahrain isolate occupies a distinct and more divergent lineage than the Chinese virus implying that epidemic RHDV strains have emerged at least twice during the past 20 years and are co-circulating in both domestic and wild rabbits.

Immunocapture-RT-PCR assay for detection and molecular epidemiology studies of Rabbit Haemorrhagic Disease and European Brown Hare Syndrome viruses

Rabbit Haemorrhagic Disease Virus and European Brown Hare Syndrome Virus are two members of the genus Lagovirus in the family Caliciviridae. They are the causative agents of highly contagious and fatal diseases of rabbits and hares respectively. We adjusted one assay for the detection and the genomic characterisation of each virus, based on viral purification by immunocapture and genomic amplification by reverse transcription-polymerase chain reaction (IC-RT-PCR). It is carried out directly with the liver exudate obtained after thawing and suppresses the viral nucleic acid preparation step. This assay combines the advantages of an ELISA test (rapidity) because immunocapture and the RT reaction were carried out in the same microtitre plate, and the advantages of PCR (sensitivity). The procedure described allows the processing of large numbers of samples and is suitable for phylogenetic studies of lagomorphs caliciviruses. In addition, it was compared with sandwich-ELISA used for Rabbit Haemorrhagic Disease or European Brown Hare Syndrome diagnosis. A good correlation was found between ELISA and IC-RT-PCR results for Rabbit Haemorrhagic Disease diagnosis, whereas for European Brown Hare Syndrome diagnosis, the results confirmed the higher sensitivity of the molecular method.

The emergence of rabbit haemorrhagic disease virus: will a non-pathogenic strain protect the UK?

Rabbit haemorrhagic disease virus emerged in China in 1984, and has killed hundreds of millions of wild rabbits in Australia and Europe. In the UK there appears to be an endemic non-pathogenic strain, with high levels of seroprevalence being recorded, in the absence of associated mortality. Using a seasonal, age-structured model we examine the hypothesis that differences in rabbit population demography differentially affect the basic reproductive rates (R(0)) of the pathogenic and non-pathogenic strains, leading to each dominating in some populations and not others. The strain with the higher R(0) excluded the other, with the dynamics depending upon the ratio of the two R(0) values. When the non-pathogenic strain dominated, the pathogenic strain caused only transient mortality, although this could be significant when the two R(0) values were similar. When the pathogenic strain dominated, repeated epidemics led to host eradication. Seroprevalence data suggest that the non-pathogenic strain may be protecting some, but not all UK populations, with half being 'at risk' from invasion by the pathogenic strain and a fifth prone to significant transient mortality. We identify key questions for empirical research to test this prediction.

Horizontal transmissible protection against myxomatosis and rabbit hemorrhagic disease by using a recombinant myxoma virus

We have developed a new strategy for immunization of wild rabbit populations against myxomatosis and rabbit hemorrhagic disease (RHD) that uses recombinant viruses based on a naturally attenuated field strain of myxoma virus (MV). The recombinant viruses expressed the RHDV major capsid protein (VP60) including a linear epitope tag from the transmissible gastroenteritis virus (TGEV) nucleoprotein. Following inoculation, the recombinant viruses induced specific antibody responses against MV, RHDV, and the TGEV tag. Immunization of wild rabbits by the subcutaneous and oral routes conferred protection against virulent RHDV and MV challenges. The recombinant viruses showed a limited horizontal transmission capacity, either by direct contact or in a flea-mediated process, promoting immunization of contact uninoculated animals.

Nucleotide sequence of UK and Australian isolates of feline calicivirus (FCV) and phylogenetic analysis of FCVs

We have determined the first complete genome sequence and capsid gene sequences of feline calicivirus (FCV) isolates from the UK and Australia. These were compared with other previously published sequences. The viruses used in the comparisons were isolated between 1957 and 1995 from various geographical locations and obtained from cats showing a range of clinical signs. Despite these diverse origins, comparisons between all strains showed a similar degree of sequence variation within both ORF1 (non-structural polyprotein) and ORF2 (major capsid protein) (amino acid distances of 7.7-13.0% and 8.8-18.6%, respectively). In contrast, ORF3 (putative minor structural protein) sequences indicated a more heterogenous distribution of FCV relatedness (amino acid distances of 1.9-17.9%). Phylogenetic analysis suggested that, unlike some other caliciviruses, FCV isolates within the current data set fall into one diverse genogroup. Within this group, there was an overall lack of geographic or temporal clustering which may be related to the epidemiology of FCV infection in cats. Analysis of regions of variability in the genome has shown that, as well as the previously identified variable regions in ORF2, similar domains exist within ORFs 1 and 3 also, although to a lesser extent. In ORF1, these variable domains largely fall between the putative non-structural protein functional domains.

A further step in the evolution of rabbit hemorrhagic disease virus: the appearance of the first consistent antigenic variant

Rabbit hemorrhagic disease virus (RHDV) is a noncultivable calicivirus that infects rabbits (Oryctolagus cuniculus) and causes epidemics of an acute fatal hepatitis. In 1997 we identified two RHDV isolates from geographically distant Italian regions, which differed antigenically from the reference strain RHDV.Bs89. In fact, they were not reactive with mAb 1H8, that is able to protect rabbits from RHD and showed a low reactivity with the rabbit convalescent serum raised against RHDV.Bs89. Experimental infection of rabbits with either RHDV isolates confirmed their high pathogenicity and their peculiar antigenic profile; nevertheless, rabbits vaccinated with the current vaccine were protected against challenge infection with these isolates. Sequence comparison definitely demonstrated that the two isolates originated from the same RHDV variant and that the similarity of their structural protein (VP60) sequences with the RHDV.Bs89 is equal to 93%. This variant was named RHDVa since shows consistent genetic and antigenic differences from the wild-type RHDV. In particular, 44% of amino acid substitutions in RHDVa VP60 were located between amino acids 344 and 370, where the similarity with RHDV.Bs89 drops to 70%, suggesting that this region probably contains the epitope recognized by mAb 1H8. In addition, this paper presents preliminary data concerning the amino acids of VP60 involved in the hemagglutination site of the virus.