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.     

Adaptive changes in the genomes of wild rabbits after 16 years of viral epidemics

Since its introduction to control overabundant invasive European rabbits (Oryctolagus cuniculus), the highly virulent rabbit haemorrhagic disease virus (RHDV) has caused regular annual disease outbreaks in Australian rabbit populations. Although initially reducing rabbit abundance by 60%, continent‐wide, experimental evidence has since indicated increased genetic resistance in wild rabbits that have experienced RHDV‐driven selection. To identify genetic adaptations, which explain the increased resistance to this biocontrol virus, we investigated genome‐wide SNP (single nucleotide polymorphism) allele frequency changes in a South Australian rabbit population that was sampled in 1996 (pre‐RHD genomes) and after 16 years of RHDV outbreaks. We identified several SNPs with changed allele frequencies within or close to genes potentially important for increased RHD resistance. The identified genes are known to be involved in virus infections and immune reactions or had previously been identified as being differentially expressed in healthy versus acutely RHDV‐infected rabbits. Furthermore, we show in a simulation study that the allele/genotype frequency changes cannot be explained by drift alone and that several candidate genes had also been identified as being associated with surviving RHD in a different Australian rabbit population. Our unique data set allowed us to identify candidate genes for RHDV resistance that have evolved under natural conditions, and over a time span that would not have been feasible in an experimental setting. Moreover, it provides a rare example of host genetic adaptations to virus‐driven selection in response to a suddenly emerging infectious disease.     

Distribution and Prevalence of the Australian Non-Pathogenic Rabbit Calicivirus Is Correlated with Rainfall and Temperature

Background

Australia relies heavily on rabbit haemorrhagic disease virus (RHDV) for the biological control of introduced European wild rabbits Oryctolagus cuniculus, which are significant economic and environmental pests. An endemic non-pathogenic rabbit calicivirus termed RCV–A1 also occurs in wild rabbits in Australian and provides partial protection against lethal RHDV infection, thus interfering with effective rabbit control. Despite its obvious importance for rabbit population management, little is known about the epidemiology of this benign rabbit calicivirus.

Methods

We determined the continent-wide distribution and prevalence of RCV-A1 by analysing 1,805 serum samples from wild rabbit populations at 78 sites across Australia for the presence of antibodies to RCV-A1 using a serological test that specifically detects RCV-A1 antibodies and does not cross-react with co-occurring RHDV antibodies. We also investigated possible correlation between climate variables and prevalence of RCV-A1 by using generalised linear mixed effect models.

Results

Antibodies to RCV-A1 were predominantly detected in rabbit populations in cool, high rainfall areas of the south-east and south-west of the continent. There was strong support for modelling RCV-A1 prevalence as a function of average annual rainfall and minimum temperature. The best ranked model explained 26% of the model structural deviance. According to this model, distribution and prevalence of RCV-A1 is positively correlated with periods of above average rainfall and negatively correlated with periods of drought.

Implications

Our statistical model of RCV-A1 prevalence will greatly increase our understanding of RCV-A1 epidemiology and its interaction with RHDV in Australia. By defining the environmental conditions associated with the prevalence of RCV-A1, it also contributes towards understanding the distribution of similar viruses in New Zealand and Europe.     

表达兔出血症病毒VP60蛋白的重组兔粘液瘤病毒的构建与初步评价

兔出血症病毒 (Rabbit hemorrhagic disease virus,RHDV) 及兔粘液瘤病毒 (Myxoma virus,MYXV) 分别引起兔出血症 (兔瘟) 和兔粘液瘤病,是两种严重危害家兔养殖业以及导致原产地欧洲野兔-穴兔 (Oryctolagus cuniculus) 种群近濒危的重要病原。VP60为构成RHDV衣壳的主要抗原蛋白。为研制能同时免疫预防该两种疫病的重组二联疫苗,本研究分别以MYXV和其复制非必需基因——胸腺激酶 (Thymidine kinase,TK) 基因为重组载体和同源重组靶基因,构建穿梭载体p7.5-VP60-GFP。将p7.5-VP60-GFP载体转染被MYXV感染的兔肾细胞株RK13,经同源重组后,在荧光显微镜下筛选出表达GFP的重组病毒,并将其命名为rMV-VP60-GFP。通过PCR和Western blotting进行重组病毒vp60基因特异性插入和表达验证结果显示,vp60和gfp基因成功插入MYXV基因组中并且可成功表达,表明成功构建了表达RHDV衣壳蛋白基因vp60的重组MYXV。动物攻毒保护试验表明,制备的重组病毒能保护家兔抵抗MYXV的致死性攻击,这为后续疫苗的研发奠定了基础。     

First case of rabbit haemorrhagic disease in Canada: contaminated flying insect, vs. long-term infection hypothesis

Following the announcement of the first case of rabbit haemorrhagic disease (RHD) in a pet rabbit, housed indoors in Canada for more than 1 year, I submitted an evidence‐based explanation to ProMed explaining how RHD might have caused the death of ‘one’ of the three pet rabbits. I suggested with supporting evidence, that it may have been persistently infected with rabbit haemorrhagic disease virus (RHDV) which may have reactivated to cause the fatal disease. However, in this issue, Peacock et al. have proposed an alternative ‘hypothesis’ for the appearance of RHD in the pet rabbit. They hypothesise that a non‐identified insect or fomite might have become contaminated by a Chinese strain of RHDV somewhere in the US. This insect/fomite then flew or was windborne, from the US to Canada where it entered the house containing three pet rabbits and infected one of them. RHD is non‐endemic and is rarely reported in the US, where it has only been observed in domestic European rabbits, held in rabbitries. My proposal was based on the details provided by ProMed, the veterinary report from Canada, where RHDV has never previously been identified and the epidemiological, ecological and evolutionary history of RHDV which includes serological and phylogenetic evidence that ancestral RHDV lineages circulated before 1984. The flying insect hypothesis of Peacock et al. is based on circumstantial evidence and, I believe, has a lower probability of being correct than my evidence‐based long‐term infection proposal.     

A sensitive and specific blocking ELISA for the detection of rabbit calicivirus RCV-A1 antibodies

ABSTRACT: BACKGROUND: Antibodies to non-pathogenic rabbit caliciviruses (RCVs) cross-react in serological tests for rabbit hemorrhagic disease virus (RHDV) and vice versa, making epidemiological studies very difficult where both viruses occur. It is important to understand the distribution and interaction of the two viruses because the highly pathogenic RHDV has been used as a biocontrol agent for wild rabbits in Australia and New Zealand for the past 17 years. The presence of the benign RCV Australia 1 (RCV-A1) is considered a key factor for the failure of RHDV mediated rabbit control in cooler, wetter areas of Australia. RESULTS: A highly sensitive and specific blocking ELISA was developed for the detection of RCV-A1 antibodies. When sera from rabbits with a known infection history for either RCV-A1 or RHDV were tested, this assay showed 100% sensitivity and no cross-reactivity with RHDV sera (100% specificity). CONCLUSIONS: This new ELISA not only allows the detection of RCV-A1 at a population level, but also permits the serological status of individual rabbits to be determined more reliably than previously described methods. This robust and simple to perform assay is therefore the tool of choice for studying RCV-A1 epidemiology in Australian wild rabbit populations.     

Oxidative stress in wild European rabbits naturally infected with myxoma virus and rabbit haemorrhagic disease virus

The European rabbit (Oryctolagus cuniculus) is one of the most important vertebrate species in the Mediterranean Basin ecosystem. Over the last 60 years, the arrival of two viral diseases, myxomatosis and rabbit haemorrhagic disease, have led to dramatic declines in wild rabbit populations across the Iberian Peninsula. These diseases are currently endemic. Periodic outbreaks occur and have significant impacts on wild populations. Both infection types have diverse physiological effects on their hosts that are rooted in aerobic metabolic processes. To fight off these viruses, rabbits activate their immune systems. However, the production of immune defences generates reactive oxygen species that may consequently damage host tissues. Hypothesising that immune responses increase oxidative stress, we examined whether wild rabbits naturally infected with myxoma virus (MV) and rabbit haemorrhagic disease virus (RHDV) had high oxidative stress. Using blood samples, we measured anti-MV and anti-RHDV antibody concentrations and different oxidative stress markers (i.e., glutathione peroxidase, glutathione reductase, superoxide dismutase, catalase, and malondialdehyde). Our results show that rabbits that were seropositive for both MV and RHDV had high concentrations of malondialdehyde. Age and body condition were also positively related to dual seropositivity. No significant relationships were observed between serostatus and the concentrations of the other oxidative stress markers. Although we expected infection with MV and RHDV to be correlated with oxidative stress, the influence of external sources of oxidative stress (e.g., climatic conditions) likely made it more difficult to detect such relationships in wild rabbits.     

Resistance to RHD virus in wild Australian rabbits: Comparison of susceptible and resistant individuals using a genomewide approach

Deciphering the genes involved in disease resistance is essential if we are to understand host–pathogen coevolutionary processes. The rabbit haemorrhagic disease virus (RHDV) was imported into Australia in 1995 as a biocontrol agent to manage one of the most successful and devastating invasive species, the European rabbit (Oryctolagus cuniculus). During the first outbreaks of the disease, RHDV caused mortality rates of up to 97%. Recently, however, increased genetic resistance to RHDV has been reported. Here, we have aimed to identify genomic differences between rabbits that survived a natural infection with RHDV and those that died in the field using a genomewide next‐generation sequencing (NGS) approach. We detected 72 SNPs corresponding to 133 genes associated with survival of a RHD infection. Most of the identified genes have known functions in virus infections and replication, immune responses or apoptosis, or have previously been found to be regulated during RHD. Some of the genes identified in experimental studies, however, did not seem to play a role under natural selection regimes, highlighting the importance of field studies to complement the genomic background of wildlife diseases. Our study provides a set of candidate markers as a tool for the future scanning of wild rabbits for their resistance to RHDV. This is important both for wild rabbit populations in southern Europe where RHD is regarded as a serious problem decimating the prey of endangered predator species and for assessing the success of currently planned RHDV variant biocontrol releases in Australia.     

Serological evidence for the presence of non-pathogenic rabbit haemorrhagic disease virus-like strains in rabbits (Oryctolagus cuniculus) of the Kerguelen archipelago

Antibodies raised against a Lagovirus were found in healthy rabbits Oryctolagus cuniculus sampled in 2003 and 2004 in the Kerguelen archipelago. The serological test we used enabled the detection of antibodies due to both pathogenic and non-pathogenic viruses related to the rabbit haemorrhagic disease virus (RHDV). The overall proportion of seropositive rabbits was 35% and differed between sites. Since previous studies have failed to detect antibodies raised against pathogenic RHDV strains, the antibodies detected in the present study were likely due to non-pathogenic strains. The interest of these strains circulating in the Kerguelen archipelago is that they have evolved independently of those detected in the other parts of the world and should derive from an ancestral RHDV precursor. Their characterization may help understanding the evolution of the virus and the emergence of the disease.     

Histo-blood group antigens act as attachment factors of rabbit hemorrhagic disease virus infection in a virus strain-dependent manner

Rabbit Hemorrhagic disease virus (RHDV), a calicivirus of the Lagovirus genus, and responsible for rabbit hemorrhagic disease (RHD), kills rabbits between 48 to 72 hours post infection with mortality rates as high as 50-90%. Caliciviruses, including noroviruses and RHDV, have been shown to bind histo-blood group antigens (HBGA) and human non-secretor individuals lacking ABH antigens in epithelia have been found to be resistant to norovirus infection. RHDV virus-like particles have previously been shown to bind the H type 2 and A antigens. In this study we present a comprehensive assessment of the strain-specific binding patterns of different RHDV isolates to HBGAs. We characterized the HBGA expression in the duodenum of wild and domestic rabbits by mass spectrometry and relative quantification of A, B and H type 2 expression. A detailed binding analysis of a range of RHDV strains, to synthetic sugars and human red blood cells, as well as to rabbit duodenum, a likely gastrointestinal site for viral entrance was performed. Enzymatic cleavage of HBGA epitopes confirmed binding specificity. Binding was observed to blood group B, A and H type 2 epitopes in a strain-dependent manner with slight differences in specificity for A, B or H epitopes allowing RHDV strains to preferentially recognize different subgroups of animals. Strains related to the earliest described RHDV outbreak were not able to bind A, whereas all other genotypes have acquired A binding. In an experimental infection study, rabbits lacking the correct HBGA ligands were resistant to lethal RHDV infection at low challenge doses. Similarly, survivors of outbreaks in wild populations showed increased frequency of weak binding phenotypes, indicating selection for host resistance depending on the strain circulating in the population. HBGAs thus act as attachment factors facilitating infection, while their polymorphism of expression could contribute to generate genetic resistance to RHDV at the population level.     

Evolution of microparasites in spatially and genetically structured host populations: The example of RHDV infecting rabbits

Several studies have shown that classical results of microparasite evolution could not extend to the case where the host species shows an important spatial structure. Rabbit haemorrhagic disease virus (RHDV), responsible for rabbit haemorrhagic disease (RHD), which recently emerged in rabbits, has strains within a wide range of virulence, thus providing an interesting example of competition between strains infecting a host species with a metapopulation structure. In addition, rabbits may show a genetic diversity regarding RHDV susceptibility. In the present paper we use the example of the rabbit-RHDV interaction to study the competition between strains of a same microparasite in a host population that is both spatially and genetically structured. Using metapopulation models we show that the evolution of the microparasite is guided by a trade-off between its capacity to invade subpopulations potentially infected by other strains and its capacity to persist within the subpopulation. In such a context, host genetic diversity acts by reducing the number of hosts susceptible to each strain, often favouring more persistent—and generally less virulent—strains. We also show that even in a stochastic context where host genes regularly go locally extinct, the microparasite pressure helps maintain the genetic diversity in the long term while reinforcing gene loss risk in the short term. Finally, we study how different demographic and epidemiologic parameters affect the coevolution between the rabbit and RHDV.     

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.     

High adaptive variability and virus-driven selection on major histocompatibility complex (MHC) genes in invasive wild rabbits in Australia

The rabbit haemorrhagic disease virus (RHDV) was imported into Australia in 1995 as a biocontrol agent to manage one of the most successful and devastating invasive species, the European rabbit (Oryctolagus cuniculus cuniculus). During the first disease outbreaks, RHDV caused mortality rates of up to 97% and reduced Australian rabbit numbers to very low levels. However, recently increased genetic resistance to RHDV and strong population growth has been reported. Major histocompatibility complex (MHC) class I immune genes are important for immune responses against viruses, and a high MHC variability is thought to be crucial in adaptive processes under pathogen-driven selection. We asked whether strong population bottlenecks and presumed genetic drift would have led to low MHC variability in wild Australian rabbits, and if the retained MHC variability was enough to explain the increased resistance against RHD. Despite the past bottlenecks we found a relatively high number of MHC class I sequences distributed over 2–4 loci. We identified positive selection on putative antigen-binding sites of the MHC. We detected evidence for RHDV-driven selection as one MHC supertype was negatively associated with RHD survival, fitting expectations of frequency-dependent selection. Gene duplication and pathogen-driven selection are possible (and likely) mechanisms that maintained the adaptive potential of MHC genes in Australian rabbits. Our findings not only contribute to a better understanding of the evolution of invasive species, they are also important in the light of planned future rabbit biocontrol in Australia.     

Atomic Model of Rabbit Hemorrhagic Disease Virus by Cryo-Electron Microscopy and Crystallography

Rabbit hemorrhagic disease, first described in China in 1984, causes hemorrhagic necrosis of the liver. Its etiological agent, rabbit hemorrhagic disease virus (RHDV), belongs to the Lagovirus genus in the family Caliciviridae. The detailed molecular structure of any lagovirus capsid has yet to be determined. Here, we report a cryo-electron microscopic (cryoEM) reconstruction of wild-type RHDV at 6.5 Å resolution and the crystal structures of the shell (S) and protruding (P) domains of its major capsid protein, VP60, each at 2.0 Å resolution. From these data we built a complete atomic model of the RHDV capsid. VP60 has a conserved S domain and a specific P2 sub-domain that differs from those found in other caliciviruses. As seen in the shell portion of the RHDV cryoEM map, which was resolved to ∼5.5 Å, the N-terminal arm domain of VP60 folds back onto its cognate S domain. Sequence alignments of VP60 from six groups of RHDV isolates revealed seven regions of high variation that could be mapped onto the surface of the P2 sub-domain and suggested three putative pockets might be responsible for binding to histo-blood group antigens. A flexible loop in one of these regions was shown to interact with rabbit tissue cells and contains an important epitope for anti-RHDV antibody production. Our study provides a reliable, pseudo-atomic model of a Lagovirus and suggests a new candidate for an efficient vaccine that can be used to protect rabbits from RHDV infection.     

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.     

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.     

Seroepidemiology of rabbit haemorrhagic disease (RHD) in wild rabbits (Oryctolagus cuniculus) in the United Kingdom

An extensive collection of blood samples from adult wild rabbits, from areas where Rabbit Haemorrhagic Disease (RHD) had not been recorded, was obtained from sites across the U.K. and parts of Eire over the winter of 1994/95. Sera from 946 animals were examined for antibodies to RHD. Antibody was found in all populations, varying between 20 and 100% depending on locality, and overall 64% of rabbits were seropositive in midwinter, supporting the view that non-pathogenic RHD or RHD-like caliciviruses, which are not producing clinical disease, circulate within rabbit populations. The serological response to these agents appears to confer significant immunity on rabbits subsequently exposed to the virulent RHD virus, so the disease will not be as devastating as myxomatosis was in the early 1950s (99% mortality) but at least one-third of the c. 40 million rabbits in the U.K. are susceptible.     

Protection against myxomatosis and rabbit viral hemorrhagic disease with recombinant myxoma viruses expressing rabbit hemorrhagic disease virus capsid protein.

Two myxoma virus-rabbit hemorrhagic disease virus (RHDV) recombinant viruses were constructed with the SG33 strain of myxoma virus to protect rabbits against myxomatosis and rabbit viral hemorrhagic disease. These recombinant viruses expressed the RHDV capsid protein (VP60). The recombinant protein, which is 60 kDa in size, was antigenic, as revealed by its reaction in immunoprecipitation with antibodies raised against RHDV. Both recombinant viruses induced high levels of RHDV- and myxoma virus-specific antibodies in rabbits after immunization. Inoculations by the intradermal route protected animals against virulent RHDV and myxoma virus challenges.