Comparative quantitative monitoring of rabbit haemorrhagic disease viruses in rabbit kittens
Virology Journalvolume 11, Article number: 109 (2014)
Only one strain (the Czech CAPM-v351) of rabbit haemorrhagic disease virus (RHDV) has been released in Australia and New Zealand to control pest populations of the European rabbit O. cuniculus. Antigenic variants of RHDV known as RHDVa strains are reportedly replacing RHDV strains in other parts of the world, and Australia is currently investigating the usefulness of RHDVa to complement rabbit biocontrol efforts in Australia and New Zealand. RHDV efficiently kills adult rabbits but not rabbit kittens, which are more resistant to RHD the younger they are and which may carry the virus without signs of disease for prolonged periods. These different infection patterns in young rabbits may significantly influence RHDV epidemiology in the field and hence attempts to control rabbit numbers.
We quantified RHDV replication and shedding in 4–5 week old rabbits using quantitative real time PCR to assess their potential to shape RHDV epidemiology by shedding and transmitting virus. We further compared RHDV-v351 with an antigenic variant strain of RHDVa in kittens that is currently being considered as a potential RHDV strain for future release to improve rabbit biocontrol in Australia.
Kittens were susceptible to infection with virus doses as low as 10 ID50. Virus growth, shedding and transmission after RHDVa infection was found to be comparable or non-significantly lower compared to RHDV. Virus replication and shedding was observed in all kittens infected, but was low in comparison to adult rabbits. Both viruses were shed and transmitted to bystander rabbits. While blood titres indicated that 4–5 week old kittens mostly clear the infection even in the absence of maternal antibodies, virus titres in liver, spleen and mesenteric lymph node were still high on day 5 post infection.
Rabbit kittens are susceptible to infection with very low doses of RHDV, and can transmit virus before they seroconvert. They may therefore play an important role in RHDV field epidemiology, in particular for virus transmission within social groups during virus outbreaks.
Rabbit haemorrhagic disease (RHD) was first described in 1984 and the infectious agent, Rabbit Haemorrhagic Disease Virus (RHDV) was later identified as a member of the Caliciviridae family of non-enveloped, positive-strand RNA viruses [1–3]. The virus is highly contagious and kills up to 95% of infected wild and domestic adult European rabbits (Oryctolagus cuniculus) causing high losses to wild and farmed rabbits [4–9]. Wild rabbits are an important part of many European ecosystems as well as a popular target for recreational hunting. Domestic rabbits are farmed worldwide for meat and fur . However, in Australia, following the introduction of wild rabbits in 1859, they quickly became and remain one of the worst introduced vertebrate pest species, causing immense losses to the agricultural industries and having a devastating impact on Australia’s biodiversity [9–14].
In 1995, the Czech CAPM-v351 strain of RHDV (herein referred to as RHDV-v351) was unintentionally released in Australia while being evaluated as a biological control agent for rabbits and subsequently deployed to control rabbit populations with enormous initial success . In recent years rabbit populations have recovered  potentially due to a combination of high levels of population immunity to RHDV-v351 [17–19], rabbits developing genetic resistance against the virus [20, 21] and the presence of a non-pathogenic rabbit calicivirus acting as an imperfect natural vaccine [22, 23].
In an effort to complement rabbit control programs in Australia, additional RHDV strains are currently being evaluated to determine their usefulness for release in addition to RHDV-v351. There are indications that members of the RHDVa subgroup of RHDV, which are antigenically distinct to classical RHDV strains like the v351 strain, are replacing the classical RHDV strains in both farmed and wild rabbits in some Asian and European countries [9, 24–27]. RHDVa has so far not been described in wild populations from the Iberian Peninsula [28–31]. It is very rare in French rabbit populations  but common in Italy where it was mainly described in domestic breeds . Since this suggests that RHDVa strains potentially have a competitive advantage over classical RHDV in the field they may be able to complement existing biological rabbit control in Australia. A possible explanation for the increased competitiveness of RHDVa strains is their reported ability to cause fatal infections in rabbits incompletely immunised with a standard RHDV vaccine . Hence, RHDVa strains might be able to partially overcome immunity to classical RHDV strains in field rabbits by avoiding antibodies to some key epitopes . It is further known that classical RHDV and RHDVa strains have a slightly different specificity for histo-blood group antigens, which may enable RHDVa strains to infect rabbits that are naturally less prone to infection with classical RHDV strains [21, 35]. RHDVa strains are hence currently being assessed in adult rabbits to examine their potential to improve rabbit biocontrol in Australia (A.J.R. unpublished).
To most reliably predict the outcome of the release of an additional RHDV strain it is essential to consider as many potential factors influencing RHDV epidemiology as possible and the current knowledge suggests that rabbit kittens could play an important role in RHDV epidemiology. Firstly, kittens are innately more resistant to fatal RHD than adult rabbits [36, 37]. At four weeks of age this protection is almost complete, and gradually wanes until they are fully susceptible at approximately 10 weeks. Maternal antibodies can protect kittens from acquiring the RHDV infection, and delay RHDV exposure in the young, but are not responsible for the innate resistance to lethal infection . The protection can be completely reversed by artificial immunosuppression, and differences in the innate immune mechanisms of young rabbits have been suggested as the underlying cause for this resistance . Secondly, virus transmission from 2 week old kittens to their dams has been observed, even though very limited pathology and virus replication was detected in livers from kittens 2–4 weeks of age via immunofluorescence detection of virus protein [40, 41]. Thirdly, RHDV genomes were detected via polymerase chain reaction (PCR) in tissues of kittens weeks after infection, and systemic and potentially persistent RHDV infections of young rabbits have been suggested , although reactivation of RHDV infection and shedding of infectious virus has not been reported. Lastly, kittens that survive RHDV infection develop antibodies that confer protection from further infection and may become part of an immune breeding population, impacting RHDV epidemiology and rabbit biocontrol [8, 38, 42, 43]. Despite their obvious relevance for field epidemiology, to this date a quantitative assessment of RHDV replication and shedding of RHDV in kittens has not been carried out.
To examine the potential role of rabbit kittens in RHDV epidemiology, we quantitatively monitored infection, seroconversion, replication and excretion of the RHDV-v351 biocontrol strain in kittens over time. To complement the evaluation of RHDVa strains in adult rabbits, we further compared the replication kinetics and the ability to spread to bystander rabbits of RHDV-v351 and a candidate RHDVa strain identified in 2008 in Inchaon, South Korea . While our results indicate that rabbit kittens may play an important role in RHD epidemiology, no significant differences in infection characteristics between the current RHDV-v351 and the RHDVa strain in rabbit kittens were detected.
RHDV infection characteristics in rabbit kittens
To assess the dynamics of viremia and virus excretion of the current Australian biocontrol strain RHDV-v351, we followed the course of infection over time in five rabbit kittens 4–5 weeks of age infected orally with 700 ID50 (as titrated in adult rabbits, A. J. R. unpublished). We monitored virus genome copy numbers in blood and rectal swabs via reverse transcription real-time PCR (RT-qPCR). Median virus genome copy numbers in blood samples were mostly below the detection limit of the assay (16 genomes) and peaked at 1 day post infection (dpi) in three of the five rabbits, and very little viremia was detectable in these kittens 3 dpi and thereafter (Figure 1A). Virus genome copies in rectal swabs peaked between day 1 and 5 and subsequently declined with little virus detectable 10 and 14 dpi (a median of 173 genomes per swab) when kittens were autopsied (Figure 1B). In contrast to the low amounts of virus observed in blood samples at late time points, we detected 350, 800, 14700 18700, 40100 genome copies (average = 15,000) per mg liver at 14 dpi (data not shown). Highest titres of IgM antibodies specific for RHDV were detected at day 7, which slowly decreased until the end of the experiment (day 14, Figure 1C). RHDV-specific IgG was first detected at 5 or 7 dpi and titres increased to day 14 post infection (Figure 1D). No signs of disease were observed during the 14 days the kittens were monitored.
To estimate to what extent we detect inoculum or virus progeny due to a productive infection when measuring genome copy numbers, we also inoculated five kittens orally with heat-inactivated virus. Heat-inactivation has been shown to best preserve genome copy numbers measurable by RT-qPCR compared to chemical inactivation procedures [44, 45]. RT-qPCR analysis showed a drop of 55% in genome copy numbers in heat-inactivated virus stock compared to untreated virus preparation (9×108 genome copies per millilitre heat-inactivated stock, equivalent to 2000 ID50/ml). After inoculation of 1 ml heat-inactivated virus, we measured only trace amounts of virus genomes in daily blood samples (Figure 2A) and in autopsy samples taken at 5 dpi (liver, mesenteric lymph node, spleen, duodenum and bile, data not shown). In contrast, rectal swabs did contain detectable viral genomes with a median of up to 730 genomes per swab between day 1 and 4 (Figure 2B). Therefore, 100-fold higher median viral genome concentrations per swab as measured for untreated virus (compare Figure 1B to Figure 2B) between 1 to 7 dpi, very likely indicate RHDV replication and do not reflect inoculated virus alone. All five kittens inoculated with heat-inactivated virus were serologically negative for RHDV-specific IgM at 5 dpi, indicating that virus replication had not occurred.
In contrast to all previous studies where very high doses (usually >1500 ID50) of virus were used to infect kittens, we next examined if kittens are also susceptible to infection with a dose as low as 10 ID50 (titrated in adult rabbits). Notably, of the five kittens inoculated, one kitten died at 3 dpi presenting typical RHD liver pathology (discoloured liver, enlarged spleen) at autopsy. The remaining four kittens did not show signs of disease and seroconverted, as indicated by IgM titres ≥ 1:1280 at 5 dpi and IgG titres > 1:20000 at 14 dpi, demonstrating that kittens do not require a substantially higher dose than adults to become productively infected.
Comparing the currently used biocontrol strain RHDV-v351 with a potential RHDVa candidate for future release
To compare RHDV-v351 and the potential candidate RHDVa strain, we infected six kittens with 700 ID50 of RHDV and in parallel six kittens with a similar amount of RHDVa (approximately 2000 ID50 as titrated in adult rabbits). Low levels of replication were observed for the RHDV-v351 strain, with peak genome copy numbers in blood after 2 dpi (Figure 3A), except for one kitten that showed a pronounced viremia until the end of the experiment (day 5). RHDVa genome copy numbers of up to 280 genome copies/μl blood were observed, but median genome concentrations were lower than observed for RHDV (Figure 3B). Genome copy numbers in swab samples were comparable at day 4 and 5 between RHDV and RHDVa (Figure 3C, D). Concentrations of genome copy numbers in mesenteric lymph node, spleen, liver, duodenum and bile did not differ significantly in RHDV and RHDVa infected kittens at autopsy, 5 dpi (Figure 3E), indicating a similar tissue specificity. In line with the low level of virus replication, we only saw signs of RHD at autopsy in one RHDV infected kitten (RHDV 3–3) that had a pale liver.
We repeated these experiments using a higher amount of inoculum equalling 4,300 ID50 of RHDV (1.6×109 genomes) and 53,000 ID50 of RHDVa (4.7×1010 genomes). Inoculation of such high amounts of virus resulted in high viral loads in the blood of two of six RHDV infected animals (2.2×105 and 6.8×105 genome copies/μl) at 2 dpi and a third kitten had 2.8×104 genomes/μl blood at day 3. While the amount of viral genomes mostly decreased in blood samples towards later time points, one RHDVa infected kitten showed increasing viremia from 3 dpi onwards, reaching 7×104 genomes/μl on day 5 (Figure 4A, B). Viral genome numbers measured in rectal swabs revealed a biphasic curve with a first peak between day one and two and another one at day 4–5. Median RHDV genome copy numbers in rectal swabs at day 4 to 5 pi were 2–7 times higher compared to RHDVa (Figure 4C, D). Notably, the >10 times higher dose of virus used to inoculate RHDVa was not reflected in the genome copy numbers measured in rectal swab samples. The amount of virus genomes detected in tissues taken 5 dpi varied considerably between individuals, with very high virus loads observed in the kittens that also had high genome copy numbers in blood. Interestingly, median genome copy numbers in duodenum and bile were very similar between strains (Figure 4E). None of the kittens showed clinical signs of RHD during this experiment, but at autopsy clear evidence of RHD (discoloured liver, enlarged spleen) was observed in two RHDV infected kittens (RHDV 4–3 and RHDV 4–4) and in kitten RHDVa 4–6, which still showed increasing viremia at the time of autopsy. The Korean RHDVa strain reached lower blood virus concentrations compared to RHDV, even when inoculated at a much higher dose. No differences in rectal virus excretion and tissue loads were obvious.
Comparative transmission studies of RHDV-v351 and RHDVa
In addition, we examined if transmission of virus to bystander kittens occurs, and if the two strains differ in their transmissibility/infectivity, by group-housing one kitten infected with RHDV and one with RHDVa with four uninfected bystanders. The experiment was carried out in triplicates. Only one of the 12 bystanders became infected with the RHDV-v351 strain, indicated by 13,000 genome copies per mg liver taken at autopsy, 5 days after being housed with the infected kittens (data not shown).In a second transmission study, two kittens infected with the high dose of RHDV and two infected with RHDVa were housed with four uninfected bystanders in one cage, 24 h post infection. Virus replication in the infected animals was monitored daily in blood samples and rectal swabs. The infected animals were euthanized at day 2 (n = 2) and day 4 (n = 2). Bystander rabbits were also monitored daily for virus replication by collecting blood and rectal swabs, and were euthanized at day 5 after their placement with the infected littermates. Livers of all animals were collected at autopsy for virus quantification. This experiment was also carried out in triplicate. Six of the 12 bystander rabbits acquired an infection, shown by detectable RNA in swab, blood and liver samples (Figure 5), demonstrating that kittens can shed infectious virus that is transmitted to littermates. Using PCR assays specific either for the RHDV or RHDVa strain, we determined that four of the six kittens that acquired an infection got infected with the RHDV strain and two with the RHDVa strain. One of the RHDVa infected bystander kittens was also positive for RHDV (with < 10% of detected genome copies being RHDV), indicating a mixed infection. It is worth noting that 3 of the 4 RHDV positive bystanders came from the same cage in which one kitten initially inoculated with RHDV (RHDV 5–5), was found dead at 4 dpi with very high amounts of viral genomes in swab, blood and liver samples (Figure 5C, F, I). One of the three bystanders (BS 5–11) in that cage also was found dead with clear signs of RHD at autopsy and another kitten (BS 5–10) had a high viremia five days after being placed with the infected kittens and a very pale liver at autopsy.
In summary, our results confirm that kittens are susceptible to productive RHDV infection after oral inoculation. RHDV infections in kittens were mostly apathogenic, with little virus detectable in blood and tissues. Virus clearance started at day 3 post infection, at the same time as IgM became detectable. Only a few kittens developed high virus loads and these also displayed signs of disease at autopsy or died. Differences between the RHDV and RHDVa strain were minor and not statistically significant with slightly lower amounts of virus in the blood and less pathogenicity detected for the RHDVa strain.
Our results verify and add to the existing knowledge about RHD in rabbit kittens: It has been shown that kittens can be infected with RHDV if high amounts of virus (at least 1500 LD50) are inoculated intra muscularly or orally [38, 40, 41, 46]. We demonstrate here productive infection and seroconversion in 4 of 5 kittens (and one death) after oral inoculation with 10 ID50, which is at least 150-times less than what had been previously used to inoculate kittens. The infectious dose for domestic seronegative rabbit kittens therefore does not appear to be substantially higher than for adult rabbits.
Virus amounts observed in livers of kittens autopsied at 5 dpi were approximately 1000-times lower compared to adult rabbits that succumbed to infection ( and A.J.R., unpublished). The only exceptions were the eight out of 82 kittens exposed to the virus that died or showed signs of disease at autopsy and may have died if the trial had been continued. In kittens showing signs of disease, virus loads were higher in liver and blood compared to kittens with subclinical infections, which demonstrates directly for the first time that pathogenicity of RHDV in rabbit kittens is linked to high levels of virus replication. Furthermore, our results indicate that in the small proportion of 4–5 week old rabbits that succumb to infection , RHDV does not cause a prolonged disease or a gradual increase of replication levels, but rather shows the same characteristics as observed in adult rabbits. Despite the low levels of replication, a strong adaptive immune response was induced in all kittens tested for seroconversion (5 of 5 kittens infected with 700 ID50 RHDV, shown in Figure 1, and 4 of 4 as mentioned in the text infected with 10 ID50 of RHDV).
While virus RNA can still be detected in tissues of convalescent rabbits months after primary infection, this RNA has never been demonstrated to be infectious [48, 49], and no reactivation of RHDV infections has been described. In our study the amount of virus shed from kittens declined with the onset of the adaptive immune response (Figure 1) and the period of virus excretion was limited to a few days after infection in surviving disease-free kittens. The likelihood that prolonged RHDV infection in kittens may represent a reservoir for RHDV between outbreaks [42, 46] hence appears low. This is also supported by a recent genetic study indicating that subsequent annual outbreaks at one site are caused by genetically different RHDV strains . Nevertheless, our results show that kittens can shed virus in sufficient amounts to potentially transmit RHDV in the field, and may therefore play an important role for transmitting the virus during an outbreak, in particular within warren systems. However, the experimental conditions mean that we cannot completely exclude the possibility that some of the virus that infected the bystanders in the transmission experiment was from inoculum shed by the orally infected kittens and not virus progeny resulting from productive infection.
An RHDVa strain with different infection characteristics compared to RHDV-v351 in rabbit kittens could significantly alter RHD epidemiology and the effectiveness of Australian viral rabbit control efforts, e.g. by causing more kittens to seroconvert or by being spread more effectively due to increased replication in kittens. As a case in point, a recently described RHDV variant termed RHDV2 appears to be replacing other field strains of RHDV and RHDVa in wild rabbit populations and rabbitries in Europe [29, 51–53]. Notably, this new variant has been reported to kill young rabbits more effectively, further pointing to an important role kittens may play in RHDV epidemiology and highlighting the need to study virus isolates considered for release as additional rabbit control agents in Australia in young rabbits. Hence, we also examined the infection characteristics of an RHDVa strain and showed that fewer kittens developed high virus loads and signs of disease after infection and there were no fatalities compared to infection with RHDV. As the amount of RHDV- v351 and RHDVa used for infection in this trial was not identical, a direct comparison is difficult to make. However, lower genome copy numbers in blood and tissues as well as less virus excretion and transmission were observed even when 12-times higher amounts of RHDVa compared to RHDV were inoculated into kittens (Figures 4 and 5), although this difference was not statistically significant. This suggests that the RHDVa strain might be less virulent in New Zealand White rabbit kittens than the current RHDV biocontrol strain, independent of the amount of virus inoculated.
Previous studies have shown that wild rabbits have an increased level of resistance to infection with RHDV-v351 . It is feasible that the replication kinetics of RHDVa also differs in wild rabbits – this will need to be assessed experimentally in future studies. However, replication and transmission properties in kittens alone are not likely to be the sole determining factor for a putative additional biocontrol strain to be considered for release in Australia. More importantly, virulence parameters in adult rabbits (wild and domestic) such as efficient replication, high mortality rates and infectivity/transmissibility, as well as the ability to partially overcome antibodies to RHDV and the benign calicivirus RCV-A1 will be key to defining a suitable additional RHDV strain able to complement rabbit biocontrol in Australia. A number of candidates are currently being investigated.
Materials and methods
Virus strains and preparations
The inoculum of the Czech CAPM-v351 RHDV strain was a commercially available virus preparation used in rabbit biocontrol (RHDV batch 1C, produced at the Elizabeth Macarthur Agricultural Institute (EMAI); the RHDVa inoculum was clarified 10% liver homogenate in PBS of an adult rabbit infected with the RHDVa strain described as 08Q712-1 . Heat-inactivated RHDV stock was generated by incubating the virus preparation in 2 ml screw cap tubes for 180 s at 85°C in a water bath to ensure a drop in infectious virus of more than 99.9% [44, 45]. RT-qPCR analysis of the heat inactivated RHDV suspension indicated a drop in detectable genomes to 45% compared to untreated virus suspension. Virus amounts inoculated into rabbits are given as multiples of ID50 determined in adult rabbits (older than 12 weeks). The RHDV and RHDVa preparations used had comparable ID50 to genome copy ratios.
Animal trials were carried out according to Australian Animal Research Acts and the Australian Code of Practice for the care and use of animals for scientific purposes at the EMAI in Australian Quarantine Inspection Service approved facilities and at the CSIRO Black Mountain site. All procedures were approved by the respective animal ethics committees (M11/09 and M10-09, as well as CESAEC 11–01 and DOMRAB-12). Rabbit kittens were bred from four New Zealand White does, which were not vaccinated against RHDV, and one vaccinated male rabbit in the CSIRO Ecosystem Sciences rabbit breeding colony. The colony is routinely monitored and confirmed free of the benign calicivirus RCV-A1 . Kittens were weaned at 28 days of age and had ad libitum access to food (commercial rabbit pellets and autoclaved lucerne hay) and water. Rabbits from the same litter were group-housed during the experiments whenever possible. At 31 days of age (+/−1 d), kittens were infected with RHDV or RHDVa and monitored daily. Blood samples were taken from the marginal ear vein (50 μl) and pipetted into 50 μl distilled water to lyse blood cells. Swab samples were taken with the copped end of Swisspers® cosmetic cotton tips. Excess faecal material was stripped, before the swabs were eluted in 400 μl PBS. To assess seroconversion of kittens infected with RHDV, additional volumes of blood were collected into Vacuettes® (Greiner Bio-One, with clot activator) and the serological status monitored via ELISA (see below). Tissue samples including mesenteric lymph node, spleen, liver and duodenum, as well as bile, were collected during autopsy after euthanasia via intracardiac barbiturate injection following anaesthesia with Zoletil (50 mg/kg). Samples were frozen at −80°C until further processed.
To compare RHDV and RHDVa strains, equal numbers of kittens from different litters were infected with one of the viruses and littermates infected with the same virus were housed together. Transmission experiments were done by placing infected kittens with uninfected littermates 24 hrs post inoculation, to reduce the chance of infecting bystanders due to spillage of inoculum.
RNA isolation was done using the Qiagen RNeasy (tissues) and Invitrogen’s PureLink Viral RNA/DNA kit (swab, blood and bile) from samples obtained during the experiments with RHDV alone (Figures 1 and 2). RNA was extracted from 50–200 μl of swab eluates, 50 μl supernatant of the 1:1 blood-water mixtures after centrifugation for 10 min at 13000 g−1, 50 μl of bile, and 50 μl of approximately 50 mg tissue lysed in 1 ml PBS via bead beating (glass beads, 50 sec, Mini Bead beater, BioSpec).
For the comparative experiments and transmission studies (Figures 3, 4 and 5), total nucleic acids were extracted from blood, rectal swabs and tissue homogenates using the MagMAX-96 RNA Isolation Kit (Cat AM 1836–5, Ambion, Texas, USA) on a KingFisher™ 96 magnetic particle handling system (Thermo Fisher Scientific, Finland). Lysis-binding buffer (130–200 μl) was placed into a deep well plate (Cat. No. 97040450, Thermo Scientific) followed by sample (50 μl blood, bile and tissue homogenate, up to 200 μl swab eluate). The addition of magnetic beads, washing and elution of nucleic acid followed the manufacturer’s instructions and the manufacturer supplied extraction programs ‘AM1836 DW std’ or ‘AM1836 DW 200v2’ for the KingFisher™ 96 system (for 50 μl or 200 μl sample volume, respectively). Nucleic acids were eluted in 30 or 50 μl H2O. Tissue homogenisation was done as described above but using Zirconia beads and a mini bead-beater (Biospec Products, Bartesville, USA).
Reverse transcription real-time PCR analysis
Standards for reverse transcription-real time PCR were generated from RHDV and RHDVa RNA, using primer 5′ CTC GCC AGT GGT ATT ATA AAT CTT AAC AC-3′ for reverse transcription. The cDNA was amplified using the same primer and 5′ GAC AAC AGG TAA TAC GAC TCA CTA TAG GGA CTG CAA CCA GTA CC-3′, which contains a T7-promoter. The PCR products were purified and used for in vitro transcription (Promega, Riboprobe kit) according to the manufacturer’s instructions. RT-qPCR was done using the One-step rt-PCR kit with Sybr Green (Bio-Rad) on a Bio-Rad CFX96/C1000 Thermal cycler platform, using the following primers: 5′-ACC CAG TAC GGC ACR GGC TCC CAA CCA C-3′ and 5′-CTA TCT CCA TGA AAC CAG ATG CAA AGG T-3′ , which amplified the respective RHDV and RHDVa fragments with comparable efficiency. Strain-specific PCR analysis was carried out as described above using 5′-GTT GTA TTC GGG CAA CCG GAA G-3′/ 5′-CTA TGG AAC ACA AGC AAA CAT CAC CA-3′ (RHDV) and 5′-GGT TGT TTA TGG TGC CTGTAA AGC A-3′/5′-GCA CAA ATT ACC TTG TCC CAA AAG TTT-3′ (RHDVa).
RT-qPCR assays started with 10 min at 50°C for reverse transcription, followed by 5 min denaturation at 95°C. Cycle conditions were 10 s at 95°C denaturation, 40 s at 63°C annealing and elongation and 10 s at 78°C data acquisition repeated 41 times and followed by a melt curve analysis (65-95°C, 0.5°C increments, 5 s per increment). Reactions with a total reaction volume of 10 μl including 1 μl of the extracted RNA were set-up manually in triplicates in 96-well plates (Bio-Rad, Multiplate™) on ice and the plates sealed with Microseal B Adhesive Sealer (Bio-Rad). For cycler control, data analysis and evaluation the CFX manager 2.0 software (Bio-Rad) was used. The assay has a detection limit of 20 genomes/reaction  and the respective threshold for genome copy numbers per mg tissue, swab or μl blood are indicated by a dotted line in the figures or the figure captions. Samples that gave a negative result in the RT-qPCR analysis were set to 1 genome copy for subsequent analysis and illustration.
Titration of RHDV-specific IgG and IgM in rabbit sera
Serum was separated from blood collected in evacuated tubes (Vacuettes®, Greiner Bio-One) from the marginal ear vein (or from the heart of the anaesthetised animals prior to euthanasia) and subsequently analysed via RHDV IgG and IgM specific ELISA  according to the published protocol .
Data analysis tools
Testing for significant differences between the RHDV and the RHDVa strain was done using two-tailed student t-tests or Wilcoxon matched-pairs signed ranked tests (if normal distribution could not be assumed) and GraphPad Prism 6. The most significant difference between the RHDV and RHDVa strain detected was p = 0.065 (exact p-value) for median virus genome copy numbers in blood samples after high dose infection (Figure 4).
Meyers G, Wirblich C, Thiel HJ: Rabbit hemorrhagic-disease virus - molecular-cloning and nucleotide sequencing of a calicivirus genome. Virology 1991, 184: 664-676. 10.1016/0042-6822(91)90436-F
Parra F, Prieto M: Purification and characterization of a calicivirus as the causative agent of a lethal hemorrhagic disease in rabbits. J Virol 1990, 64: 4013-4015.
Ohlinger VF, Haas B, Thiel HJ: Rabbit hemorrhagic disease (RHD): characterization of the causative calicivirus. Vet Res 1993, 24: 103-116.
Chasey D, Trout RC, Edwards S: Susceptibility of wild rabbits (Oryctolagus cuniculus) in the United Kingdom to rabbit haemorrhagic disease (RHD). Vet Res 1997, 28: 271-276.
Lavazza A, Scicluna MT, Capucci L: Susceptibility of hares and rabbits to the European brown hare syndrome virus (EBHSV) and rabbit haemorrhagic disease virus (RHDV) under experimental conditions. Zentralbl Veterinarmed B 1996, 43: 401-410.
Marchandeau S, Chantal J, Portejoie Y, Barraud S, Chaval Y: Impact of viral hemorrhagic disease on a wild population of European rabbits in France. J Wildl Dis 1998, 34: 429-435. 10.7589/0090-3558-34.3.429
Mutze G, Cooke B, Alexander P: The initial impact of rabbit hemorrhagic disease on European rabbit populations in South Australia. J Wildl Dis 1998, 34: 221-227. 10.7589/0090-3558-34.2.221
Xu ZJ, Chen WX: Viral haemorrhagic disease in rabbits: a review. Vet Res Commun 1989, 13: 205-212. 10.1007/BF00142046
McIntosh MT, Behan SC, Mohamed FM, Lu Z, Moran KE, Burrage TG, Neilan JG, Ward GB, Botti G, Capucci L, Metwally SA: A pandemic strain of calicivirus threatens rabbit industries in the Americas. Virol J 2007, 4: 96. 10.1186/1743-422X-4-96
Delibes-Mateos M, Delibes M, Ferreras P, Villafuerte R: Key role of European rabbits in the conservation of the Western Mediterranean basin hotspot. Conserv Biol 2008, 22: 1106-1117. 10.1111/j.1523-1739.2008.00993.x
Sandell PR: Implications of rabbit haemorrhagic disease for the short-term recovery of semi-arid woodland communities in north-west Victoria. Wildlife Res 2002, 29: 591-598. 10.1071/WR00089
Ward VK, Cooke BD, Strive T: Rabbit Haemorrhagic Disease Virus and other Lagoviruses. In Caliciviruses. Edited by: Green K, Hansman G. Norfolk, UK: Caister Academic Press; 2011.
Cooke BD: Rabbits: manageable environmental pests or participants in new Australian ecosystems? Wildlife Res 2012, 39: 279-289.
Rolls EC: They all ran wild. Angus and Robertson: Melbourne; 1969.
Cooke BD, Fenner F: Rabbit haemorrhagic disease and the biological control of wild rabbits, Oryctolagus cuniculus, in Australia and New Zealand. Wildlife Res 2002, 29: 689-706. 10.1071/WR02010
Saunders G, Cooke B, McColl K, Shine R, Peacock T: Modern approaches for the biological control of vertebrate pests: an Australian perspective. Biol Control 2010, 52: 288-295. 10.1016/j.biocontrol.2009.06.014
Cooke BD, Robinson AJ, Merchant JC, Nardin A, Capucci L: Use of ELISAs in field studies of rabbit haemorrhagic disease (RHD) in Australia. Epidemiol Infect 2000, 124: 563-576. 10.1017/S0950268899003994
McPhee SR, Butler KL, Kovaliski J, Mutze G, Capucci L, Cooke BD: Antibody status and survival of Australian wild rabbits challenged with rabbit haemorrhagic disease virus. Wildlife Res 2009, 36: 447-456. 10.1071/WR08137
Mutze G, Kovaliski J, Butler K, Capucci L, McPhee S: The effect of rabbit population control programmes on the impact of rabbit haemorrhagic disease in south-eastern Australia. J Appl Ecol 2010, 47: 1137-1146. 10.1111/j.1365-2664.2010.01844.x
Elsworth PG, Kovaliski J, Cooke BD: Rabbit haemorrhagic disease: are Australian rabbits (Oryctolagus cuniculus) evolving resistance to infection with Czech CAPM 351 RHDV? Epidemiol Infect 2012, 140: 1972-1981. 10.1017/S0950268811002743
Nystrom K, Le Gall-Recule G, Grassi P, Abrantes J, Ruvoen-Clouet N, Le Moullac-Vaidye B, Lopes AM, Esteves PJ, Strive T, Marchandeau S, Dell A, Haslam SM, LePendu J: Histo-blood group antigens act as attachment factors of rabbit hemorrhagic disease virus infection in a virus strain-dependent manner. PloS Pathog 2011, 7: 8.
Strive T, Elsworth P, Liu JN, Wright JD, Kovaliski J, Capucci L: The non-pathogenic Australian rabbit calicivirus RCV-A1 provides temporal and partial cross protection to lethal rabbit haemorrhagic disease virus infection which is not dependent on antibody titres. Vet Res 2013, 44: 51. 10.1186/1297-9716-44-51
Strive T, Wright JD, Robinson AJ: Identification and partial characterisation of a new lagovirus in Australian wild rabbits. Virology 2009, 384: 97-105. 10.1016/j.virol.2008.11.004
Capucci L, Fallacara F, Grazioli S, Lavazza A, Pacciarini ML, Brocchi E: A further step in the evolution of rabbit hemorrhagic disease virus: the appearance of the first consistent antigenic variant. Virus Res 1998, 58: 115-126. 10.1016/S0168-1702(98)00106-3
Kerr PJ, Kitchen A, Holmes EC: Origin and phylodynamics of rabbit hemorrhagic disease virus. J Virol 2009, 83: 12129-12138. 10.1128/JVI.01523-09
Oem JK, Lee KN, Roh IS, Lee KK, Kim SH, Kim HR, Park CK, Joo YS: Identification and characterization of rabbit hemorrhagic disease virus genetic variants isolated in Korea. J Vet Med Sci 2009, 71: 1519-1523. 10.1292/jvms.001519
Wang X, Hao H, Qiu L, Dang R, Du E, Zhang S, Yang Z: Phylogenetic analysis of rabbit hemorrhagic disease virus in China and the antigenic variation of new strains. Arch Virol 2012, 157: 1523-1530. 10.1007/s00705-012-1340-9
Alda F, Gaitero T, Suarez M, Merchan T, Rocha G, Doadrio I: Evolutionary history and molecular epidemiology of rabbit haemorrhagic disease virus in the Iberian Peninsula and Western Europe. BMC Evol Biol 2010, 10: 347. 10.1186/1471-2148-10-347
Abrantes J, Lopes AM, Dalton KP, Melo P, Correia JJ, Ramada M, Alves PC, Parra F, Esteves PJ: New variant of rabbit hemorrhagic disease virus, Portugal, 2012–2013. Emerg Infect Dis 2013, 19: 1900-1902.
Abrantes J, van der Loo W, Le Pendu J, Esteves PJ: Rabbit haemorrhagic disease (RHD) and rabbit haemorrhagic disease virus (RHDV): a review. Vet Res 2012, 43: 12. 10.1186/1297-9716-43-12
Muller A, Freitas J, Silva E, Le Gall-Recule G, Zwingelstein F, Abrantes J, Esteves PJ, Alves PC, Loo WD, Kolodziejek J, Nowotny N, Thompson G: Evolution of rabbit haemorrhagic disease virus (RHDV) in the European rabbit (Oryctolagus cuniculus) from the Iberian Peninsula. Vet Microbiol 2009, 135: 368-373. 10.1016/j.vetmic.2008.09.057
Le Gall-Recule G, Zwingelstein F, Laurent S, de Boisseson C, Portejoie Y, Rasschaert D: Phylogenetic analysis of rabbit haemorrhagic disease virus in France between 1993 and 2000, and the characterisation of RHDV antigenic variants. Arch Virol 2003, 148: 65-81. 10.1007/s00705-002-0908-1
Le Gall-Recule G, Lavazza A, Marchandeau S, Bertagnoli S, Zwingelstein F, Cavadini P, Martinelli N, Lombardi G, Guerin J-L, Lemaitre E, Decors A, Boucher S, Le Normand B, Capucci L: Emergence of a new lagovirus related to rabbit haemorrhagic disease virus. Vet Res 2013, 44: 1-13. 10.1186/1297-9716-44-1
Schirrmeier H, Reimann I, Kollner B, Granzow H: Pathogenic, antigenic and molecular properties of rabbit haemorrhagic disease virus (RHDV) isolated from vaccinated rabbits: detection and characterization of antigenic variants. Arch Virol 1999, 144: 719-735. 10.1007/s007050050538
Nystrom K, Le Moullac-Vaidye B, Ruvoen-Clouet N, Le Pendu J: Shared human/rabbit ligands for rabbit hemorrhagic disease virus. Emerg Infect Dis 2012, 18: 518-519. 10.3201/eid1803.111402
Morisse JP, Le Gall G, Boilletot E: Hepatitis of viral origin in Leporidae: introduction and aetiological hypotheses. Rev Sci Tech 1991, 10: 269-310.
Villafuerte R, Calvete C, Gortazar C, Moreno S: First epizootic of rabbit hemorrhagic disease in free living populations of Oryctolagus cuniculus at Donana National Park, Spain. J Wildl Dis 1994, 30: 176-179. 10.7589/0090-3558-30.2.176
Robinson AJ, So PTM, Muller WJ, Cooke BD, Capucci L: Statistical models for the effect of age and maternal antibodies on the development of rabbit haemorrhagic disease in Australian wild rabbits. Wildlife Res 2002, 29: 663-671. 10.1071/WR00119
Marques RM, Teixeira L, Aguas AP, Ribeiro JC, Costa-e-Silva A, Ferreira PG: Immunosuppression abrogates resistance of young rabbits to Rabbit Haemorrhagic Disease (RHD). Vet Res 2014, 45: 14. 10.1186/1297-9716-45-14
Prieto JM, Fernandez F, Alvarez V, Espi A, Garcia Marin JF, Alvarez M, Martin JM, Parra F: Immunohistochemical localisation of rabbit haemorrhagic disease virus VP-60 antigen in early infection of young and adult rabbits. Res Vet Sci 2000, 68: 181-187. 10.1053/rvsc.1999.0357
Mikami O, Park JH, Kimura T, Ochiai K, Itakura C: Hepatic lesions in young rabbits experimentally infected with rabbit haemorrhagic disease virus. Res Vet Sci 1999, 66: 237-242. 10.1053/rvsc.1998.0266
Shien JH, Shieh HK, Lee LH: Experimental infections of rabbits with rabbit haemorrhagic disease virus monitored by polymerase chain reaction. Res Vet Sci 2000, 68: 255-259. 10.1053/rvsc.1999.0372
Ferreira PG, Dinis M, Costa ESA, Aguas AP: Adult rabbits acquire resistance to lethal calicivirus infection by adoptive transfer of sera from infected young rabbits. Vet Immunol Immunopathol 2008, 121: 364-369. 10.1016/j.vetimm.2007.09.005
Duizer E, Bijkerk P, Rockx B, De Groot A, Twisk F, Koopmans M: Inactivation of caliciviruses. Appl Environ Microbiol 2004, 70: 4538-4543. 10.1128/AEM.70.8.4538-4543.2004
Baert L, Wobus CE, Van Coillie E, Thackray LB, Debevere J, Uyttendaele M: Detection of murine norovirus 1 by using plaque assay, transfection assay, and real-time reverse transcription-PCR before and after heat exposure. Appl Environ Microbiol 2008, 74: 543-546. 10.1128/AEM.01039-07
Ferreira PG, Costa ESA, Oliveira MJ, Monteiro E, Aguas AP: Leukocyte-hepatocyte interaction in calicivirus infection: differences between rabbits that are resistant or susceptible to rabbit haemorrhagic disease (RHD). Vet Immunol Immunopathol 2005, 103: 217-221. 10.1016/j.vetimm.2004.09.028
Strive T, Wright J, Kovaliski J, Botti G, Capucci L: The non-pathogenic Australian lagovirus RCV-A1 causes a prolonged infection and elicits partial cross-protection to rabbit haemorrhagic disease virus. Virology 2010, 398: 125-134. 10.1016/j.virol.2009.11.045
Gall A, Hoffmann B, Teifke JP, Lange B, Schirrmeier H: Persistence of viral RNA in rabbits which overcome an experimental RHDV infection detected by a highly sensitive multiplex real-time RT-PCR. Vet Microbiol 2007, 120: 17-32. 10.1016/j.vetmic.2006.10.006
Gall A, Schirrmeier H: Persistence of rabbit haemorrhagic disease virus genome in vaccinated rabbits after experimental infection. J Vet Med B Infect Dis Vet Public Health 2006, 53: 358-362. 10.1111/j.1439-0450.2006.00986.x
Kovaliski J, Sinclair R, Mutze G, Peacock D, Strive T, Abrantes J, Esteves P, Holmes E: Molecular epidemiology of rabbit haemorrhagic disease virus (RHDV) in Australia: when one became many. Mol Ecol 2014, 23: 408-420. 10.1111/mec.12596
Dalton KP, Nicieza I, Balseiro A, Muguerza MA, Rosell JM, Casais R, Alvarez AL, Parra F: Variant rabbit hemorrhagic disease virus in young rabbits, Spain. Emerg Infect Dis 2012, 18: 2009-2012. 10.3201/eid1812.120341
Le Gall-Recule G, Lavazza A, Marchandeau S, Bertagnoli S, Zwingelstein F, Cavadini P, Martinelli N, Lombardi G, Guerin JL, Lemaitre E, Decors A, Boucher S, Le Normand B, Capucci L: Emergence of a new lagovirus related to rabbit haemorrhagic disease virus. Vet Res 2013, 44: 81. 10.1186/1297-9716-44-81
Dalton KP, Nicieza I, Abrantes J, Esteves PJ, Parra F: Spread of new variant RHDV in domestic rabbits on the Iberian Peninsula. Vet Microbiol 2014, 169: 67-73. 10.1016/j.vetmic.2013.12.015
Liu J, Kerr PJ, Wright JD, Strive T: Serological assays to discriminate rabbit haemorrhagic disease virus from Australian non-pathogenic rabbit calicivirus. Vet Microbiol 2012, 157: 345-354. 10.1016/j.vetmic.2012.01.018
Capucci L, Frigoli G, Ronshold L, Lavazza A, Brocchi E, Rossi C: Antigenicity of the rabbit hemorrhagic disease virus studied by its reactivity with monoclonal antibodies. Virus Res 1995, 37: 221-238. 10.1016/0168-1702(95)00033-M
The authors would like to thank Peter Kirkland for helpful comments and providing facilities, Dr Oem from National Veterinary Research and Quarantine Service, Ministry of Agriculture, Animal Disease Diagnostic Center Division, Anayang, Republic of Korea for supplying the RHDVa strain and Ian McDonald and Mathias Drewer for technical assistance with serological assays.
The authors declare that they have no competing interest.
Designed study: MM, TS, PJK; Carried out experiments: MM, TS, PJK, AJR, SH, JDW, PH; Analysed and interpreted data: MM, TS, PJK; Wrote Manuscript: MM, TS, PJK; contributed reagents/analysis tools: ARJ, PH. All authors read and approved the final manuscript.