jueves, 21 de diciembre de 2017

Detection and Circulation of a Novel Rabbit Hemorrhagic Disease Virus in Australia - Volume 24, Number 1—January 2018 - Emerging Infectious Disease journal - CDC

Detection and Circulation of a Novel Rabbit Hemorrhagic Disease Virus in Australia - Volume 24, Number 1—January 2018 - Emerging Infectious Disease journal - CDC

Volume 24, Number 1—January 2018


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

Jackie E. Mahar, Andrew J. Read, Xingnian Gu, Nadya Urakova1, Roslyn Mourant, Melissa Piper, Stéphanie Haboury, Edward C. Holmes, Tanja Strive, and Robyn N. HallComments to Author 
Author affiliations: Commonwealth Scientific and Industrial Research Organisation (CSIRO), Acton, Australian Capital Territory, Australia (J.E. Mahar, N. Urakova, R. Mourant, M. Piper, S. Haboury, T. Strive, R.N. Hall)The University of Sydney, Sydney, New South Wales, Australia (J.E. Mahar, E.C. Holmes)Elizabeth Macarthur Agricultural Institute, Menangle, New South Wales, Australia (A.J. Read, X. Gu)Invasive Animals Cooperative Research Centre, Bruce, Australian Capital Territory, Australia (N. Urakova, T. Strive, R.N. Hall)


The highly virulent rabbit hemorrhagic disease virus (RHDV) has been widely used in Australia and New Zealand since the mid-1990s to control wild rabbits, an invasive vertebrate pest in these countries. In January 2014, an exotic RHDV was detected in Australia, and 8 additional outbreaks were reported in both domestic and wild rabbits in the 15 months following its detection. Full-length genomic analysis revealed that this virus is a recombinant containing an RHDVa capsid gene and nonstructural genes most closely related to nonpathogenic rabbit caliciviruses. Nationwide monitoring efforts need to be expanded to assess if the increasing number of different RHDV variants circulating in the Australian environment will affect biological control of rabbits. At the same time, updated vaccines and vaccination protocols are urgently needed to protect pet and farmed rabbits from these novel rabbit caliciviruses.
Rabbit hemorrhagic disease virus (RHDV) is a calicivirus in the genus Lagovirus causing acute hepatic necrosis in European rabbits (Oryctolagus cuniculus), leading to disseminated intravascular coagulation and death in >90% of susceptible animals (1). RHDV was first reported in China in 1984 and subsequently spread to South Korea by 1985, various European countries by 1986, and Mexico in 1988 (1). RHDV is now enzootic in domestic and wild rabbits in parts of Asia and Europe, and sporadic outbreaks continue to occur in the Americas, the Middle East, and Africa (1). RHDV also has had substantial ecological impacts in certain regions such as the Iberian Peninsula, where the rabbit is a keystone species, serving as a major prey species and actively altering the environment through grazing and burrowing behaviors (2).
Conversely, in Australia and New Zealand, the European rabbit is a major vertebrate pest, threatening the survival of native plants and animals, facilitating erosion through burrowing and grazing, and causing massive economic losses to agricultural sectors. In 1991, a strain of RHDV from Czechoslovakia (CAPM V-351) was imported into Australia for assessment of its suitability as a rabbit biocontrol agent; it has been used for this purpose since 1995, resulting in substantial benefits for both the economy and the environment (35). The same CAPM V-351 strain appeared in New Zealand in 1997, after it was presumably deliberately released by farmers (6). Circulating field strains (descendants of CAPM V-351) now cause regular natural outbreaks in both countries, and the original CAPM V-351 strain is continually rereleased through rabbit control programs (7).
Thumbnail of Genomic organization of RHDV detected in rabbits, Australia. A) RHDV has a polyadenylated single-stranded positive-sense gRNA of ≈7.4 kb consisting of 2 ORFs (open boxes), with the VP (open circle) covalently attached to the 5′ end. ORF1 encodes a polyprotein that is proteolytically cleaved to produce the major capsid protein, VP60, and the nonstructural proteins p16, p23, 2C-like protein (2C-like), p29, VPg, the viral protease, and the RdRp. ORF2 encodes the minor structural protei
Figure 1. Genomic organization of RHDV. Top panel: RHDV has a polyadenylated single-stranded positive-sense gRNA of ≈7.4 kb consisting of 2 ORFs (open boxes), with the VP (open circle) covalently attached to the...
RHDV has a 7.4-kb single-stranded positive-sense RNA genome (gRNA) comprising 2 open reading frames (ORFs) (1). ORF1 encodes a single polyprotein that is cleaved by the viral protease into 7 nonstructural proteins and the major capsid protein, viral protein (VP) 60, whereas ORF2 encodes a minor structural protein, (VP10) (Figure 1) (1). Ninety dimers of VP60 self-assemble to form the RHDV capsid, whereas VP10 is present in much smaller amounts, and its function is not well established (1,8). In addition to the gRNA, a 2.1-kb subgenomic RNA (sgRNA) can also be detected in virus preparations and encodes both structural proteins VP60 and VP10 (Figure 1) (1). The first 16 nt at the 5′ ends of the RHDV gRNA and sgRNA are identical, with the exception of 2 sites; this structure may facilitate template switching at the junction of RNA-dependent RNA polymerase (RdRp) and VP60 (9).
RHDV strains are currently all classified as a single serotype, although phylogenetic analyses have led to the viruses being subclassified into 6 genogroups (10,11). The CAPM V-351 strain used for biocontrol in Australia and its descendants are genogroup G2 viruses. In the late 1990s, new genetically and antigenically distinct RHDV variants were reported and were designated as subtype RHDVa or genogroup G6 on the basis of differences in reactivity profiles with monoclonal antibodies and on sequence divergence (10,12,13). Despite this difference, vaccines developed for classic RHDV strains were still found to be protective against RHDVa when administered at the correct dose and within recommended intervals (10,12,13). RHDVa strains were first detected in Italy and were subsequently detected in multiple countries, including Portugal, the United States, China, Cuba, and South Korea (12; reviewed in 1). In some geographic regions, it was suggested that RHDVa had replaced previous enzootic viruses (1417).
Subsequently, a novel lagovirus, designated RHDV2, was reported from France in 2010 (18). This virus showed greater genetic and antigenic diversity than RHDVa compared with classic RHDV, with ≈82% nucleotide sequence similarity in VP60 between wild-type RHDV2 and RHDV (11,18,19). Early RHDV2 pathogenesis studies demonstrated that this virus had a longer disease duration, lower virulence, and lower mortality rates than RHDV and that it could cause disease in vaccinated rabbits (18). RHDV2 rapidly spread to Italy, Spain, Portugal, Great Britain, Scotland, the Azores, and Australia, and ongoing surveillance suggests that in many areas RHDV2 replaced existing circulating viruses, including RHDVa (1826). RHDV2 has also been shown to cause disease in other lagomorph species, in contrast to RHDV and RHDVa, which are considered species-specific for the European rabbit (2730).
In addition to the virulent RHDV, RHDVa, and RHDV2, several nonpathogenic lagoviruses have been described that are genetically distinct and have historically been called rabbit caliciviruses (RCVs). RCVs have been reported from various geographic locations; in contrast to the pathogenic hepatotropic lagoviruses, most RCVs mainly cause a nonclinical infection of the small intestine (1). Until 2014, only the benign rabbit calicivirus Australia-1 (RCV-A1) and the RHDV biocontrol strain CAPM V-351 (and its descendants) were circulating in Australia (7,31). In May 2015, RHDV2 was detected in Australia for the first time (23). Before the detection of RHDV2, another exotic RHDV incursion was reported to the World Organisation for Animal Health (OIE) in January 2014 (32). We describe the detection, spread, and evolution of this first exotic virus, RHDVa, in Australia during December 2013–May 2015, before the arrival of RHDV2.

Materials and Methods

Liver samples collected from rabbits suspected to have died from RHDV were submitted to the Elizabeth Macarthur Agricultural Institute (EMAI) (Menangle, New South Wales [NSW], Australia) or to CSIRO (Black Mountain, Australian Capital Territory [ACT], Australia) by veterinarians or NSW Local Land Service rangers. Samples were collected during December 1, 2013–April 30, 2015. No animal ethics permit is required in Australia for sampling of rabbits that are found dead.
For serologic analysis of healthy wild rabbits, we collected 20 wild rabbits from the ACT; these animals were shot from a vehicle using a 0.22-caliber rifle targeting the head or chest. We collected samples (blood, liver, duodenum, and bile) postmortem; all sampling was conducted according to the Australian Code for the Care and Use of Animals for Scientific Purposes as approved by the CSIRO Ecosystem Sciences animal ethics committee (approval #12–15).
RHDV Antigen Capture ELISA
Liver samples submitted to EMAI were prepared by homogenizing ≈1 g of tissue in 5 mL of phosphate-buffered saline and centrifuging at 3,000 rpm for 10 minutes. Supernatant was tested in the RHDV antigen-capture ELISA as previously reported (33).
Real-time Reverse Transcription PCR Assay
We screened samples for RHDV and RHDVa using real-time reverse transcription PCR (rRT-PCR). We applied sterile cotton swabs to the freshly cut surface of liver samples and transferred them to 5-mL vials containing 3 mL of sterile phosphate-buffered gelatin saline, then extracted viral RNA from 50-µL samples from these vials using the MagMAX-96 Viral Isolation System (Ambion, Austin, TX, USA) on a Kingfisher 96 magnetic particle handling system (ThermoFisher Scientific, Waltham, MA, USA).
We tested purified viral RNA (5 µL) in an RHDV-specific Taqman rRT-PCR assay using the AgPath-ID One-Step RT-PCR kit (ThermoFisher Scientific), based on detection of a conserved sequence of the VP60 gene with modifications (Table 1) (34). We developed an in-house Taqman rRT-PCR to specifically detect the exotic RHDVa virus (Table 1). Cycling conditions are available from the authors upon request.
We extracted RNA from 20–30 mg of liver tissue or bile using the RNeasy Mini Kit (QIAGEN, Hilden, Germany); the Maxwell 16 LEV simplyRNA Tissue Kit (Promega, Madison, WI, USA); or the Invitrogen Purelink viral RNA/DNA Mini Kit (ThermoFisher Scientific) following the manufacturers’ instructions. We performed first-strand cDNA synthesis from 5 μL of RNA with 500 ng of OligodT (18mer) (Geneworks, Thebarton, South Australia, Australia) using Invitrogen Superscript III (Life Technologies, Carlsbad, CA, USA) following the manufacturer’s directions. We conducted initial detection and typing of RHDV samples submitted to CSIRO (>1 isolate from each outbreak) using the universal lagovirus RT-PCR, as described previously (35), followed by Sanger sequencing of amplified products. RHDV G2 genomes were amplified as described previously (7). We amplified overlapping fragments of the RHDVa genome using Platinum Taq DNA polymerase high fidelity (Life Technologies) using primer sets detailed in Table 2 (cycling conditions available from the authors upon request). We purified PCR products using the QIAquick PCR purification kit (QIAGEN), quantified with the Qubit dsDNA BR assay (Life Technologies) and pooled in equimolar ratios.
Library Preparation and Illumina Sequencing
We prepared DNA libraries of PCR amplicons using the Nextera XT DNA Sample Prep Kit (Illumina, San Diego, CA, USA). Sequencing was performed on an Illumina MiSeq using the 300 cycle paired-end MiSeq Reagent Kit v2 (Illumina), as described previously (7).
NGS Genome Assembly and Data Analyses
We performed read quality assessment and trimming and merged overlapping paired-end reads, as described previously (7). We then mapped individual cleaned reads to the RHDV reference genome (Genbank accession no. M67473.1) using the Geneious in-built mapper tool as available in Geneious version 8.1.6 (39) and generated a majority consensus sequence. Primer sequences were trimmed from the 5′ and 3′ genomic termini. We named sequences using the syntax country/state/isolate identifier/date, with date in the format yyyy/mm (e.g., AUS/NSW/BER-1/2013/12). Consensus sequences were aligned with representative RHDV, RHDV2, and RCV-A1 sequences for recombination and phylogenetic analyses. We submitted all sequences generated in this study to GenBank (accession nos. KY628306–21.
We inferred maximum likelihood phylogenies for both the nonstructural genes (44 sequences, 5,238 nt) and the VP60 gene (47 sequences, 1,743 nt) using PhyML (40). We estimated trees using the general time reversible model with rate heterogeneity among sites (4 discrete rate categories) and a proportion of invariant sites (determined as most appropriate using jModel test [41]), with a combination of nearest neighbor interchange and subtree pruning and regrafting branch swapping methods employed to search for the optimal topology. We rooted trees using a European brown hare syndrome virus sequence (EBHSV; GenBank accession no. Z69620) and estimated branch support using 1,000 bootstrap replicates.
We screened a full genome alignment containing the 15 newly acquired sequences and 12 reference genomes (27 sequences, 7,303 nt) for recombination using the RDP, GENECONV, MaxChi, and BootScan methods available in the RDP4 package (42), with p<0.05 representing significant evidence of recombination. We confirmed recombination events by phylogenetic analyses on either side of the proposed breakpoint.
Serologic Analysis
We tested serum specimens from the shot wild rabbits for the presence of RCV-A1 antibodies using a specific blocking ELISA (43). We also tested for antibodies against RHDV using a competition ELISA and used ELISAs for RHDV IgA and IgM subclass antibodies as described previously (44).


In late December 2013 and January 2014, a property in northern Sydney, NSW, Australia, experienced the sudden deaths of 30 domestic rabbits from a population of 80. The treating veterinarian submitted liver samples to EMAI for routine RHDV testing in January 2014. The samples showed little or no reactivity in the RHDV antigen capture ELISA but tested strongly positive in the RHDV-specific rRT-PCR. This result was in contrast to the typical reactivity of field strains of RHDV from Australia, which react strongly in both the RHDV antigen capture ELISA and the RHDV-specific rRT-PCR (Table 3). To further characterize this calicivirus, a 326 bp fragment of the VP60 capsid gene was amplified using a universal lagovirus PCR (35). Sanger sequencing of this fragment identified the virus as a member of the RHDVa group (RHDV G6), a virus not previously reported in Australia.
Thumbnail of RHDV and RHDVa detections in Australia, January 2014–March 2015. Sites where RHDVa and Australian RHDV field strains were detected are indicated on the map and numbered according to the order in which the outbreaks occurred. Inset shows location of NSW and ACT in Australia. ACT, Australian Capital Territory; NSW, New South Wales.
Figure 2. RHDV and RHDVa detections in eastern Australia, January 2014–March 2015. Sites where RHDVa and Australian RHDV field strains were detected are indicated on the map and numbered according to the order...
From January 2014 through March 2015, the exotic RHDVa strain was detected by rRT-PCR in 7 reported outbreaks in domestic rabbit breeding facilities (outbreaks 1, 2, 3, 5, 8, 9, and 10; Table 4Figure 2), as well as in wild rabbits from 2 locations (outbreaks 7 and 11; Table 4Figure 2). Overall, >70 rabbit deaths were reported in domestic rabbits from properties where RHDVa was confirmed.
Outbreak 11, from Mulligan’s Flat Nature Reserve near Canberra, ACT, Australia, was suspected when rangers detected a sudden reduction in rabbit numbers during their rabbit eradication program, although no rabbit carcasses were recovered. Subsequent serologic analysis of 20 healthy shot rabbits revealed a high proportion of rabbits positive for RHDV-specific IgM (7/20) and IgA (15/20), indicating a recent virus outbreak (45). However, despite moderate to high IgA and IgM titers, the serum specimens tested low or negative in the RHDV competition ELISA and RCV-A1 blocking ELISA, which are known to be more strain-specific than the IgM and IgA ELISAs (data not shown). This finding indicated possible exposure of the Mulligan’s Flat rabbit population to a lagovirus that was antigenically distinct from both RHDV and RCV-A1. Because RHDV RNA can be detected in the bile of recovering animals for as long as 15 weeks postinfection (34), we extracted RNA from the bile of IgM-positive individuals and analyzed it using the universal lagovirus PCR (35). We amplified a lagovirus-specific fragment from 1 animal (AUS/ACT/MF-109/2015/03; Table 4), and Sanger sequencing confirmed the presence of the exotic RHDVa in this rabbit.
We also detected Australian field strains (RHDV G2) during the same sampling period on 3 occasions, 2 from wild rabbits (outbreaks 4 and 6; Table 4), and 1 from a domestic rabbit facility (outbreak 1; Table 4). The strain from the domestic rabbit facility was identified in samples from the first reported outbreak of the exotic RHDVa, indicating that both viruses were active simultaneously on this property.
Thumbnail of RDP plot (42) demonstrating recombination in a representative RHDVa strain from Australia. The pairwise identity of the recombinant, KYO-1, with the putative parental strains, RHDVa/AB300693.2/JPN/Hokkaido/2002 (black) and RCV-A1/EU871528.1/AUS/MIC-07(1–4)/2007 (dark gray), is plotted according to genome position (nt). A clear crossover event can be observed at the junction of RNA-dependent RNA polymerase and viral protein 60. The window size was set to 30. A schematic representatio
Figure 3. Recombination detection program plot (42) demonstrating recombination in a representative rabbit hemorraghic disease virus type a (RHDVa) strain from Australia. The pairwise identity of the recombinant, KYO-1, with the...
Thumbnail of Phylogenetic analysis of nonstructural (n = 44) and VP60 capsid (n = 47) genes of RHDV strains from Australia and reference sequences. Maximum likelihood phylogenies of the A) VP60 capsid genes and B) nonstructural genes were prepared from an alignment of the newly sequenced RHDV samples (bold) along with published sequences (accession numbers of published sequences indicated in the taxa name). The JN165235/FP/China/2009 and JN165234/XA/China/2010 sequences were restricted to the ca
Figure 4. Phylogenetic analysis of viral protein 60 (VP60) capsid (n = 47) and nonstructural (n = 44) genes of RHDV strains from Australia and reference sequences. Maximum likelihood phylogenies of the A)...
We conducted full-genome sequencing for 15 RHDV isolates we sampled during December 2013–March 2015. We could not recover a full-length genome from AUS/ACT/MF-109/2015/03. Recombination and phylogenetic analyses of the nonstructural genes indicated that all RHDVa viruses were recombinants between RHDVa and an RCV-A1-like virus, with strong bootstrap support (Figures 34). Specifically, we found strong evidence for recombination with a putative breakpoint located between the RdRp and VP60 genes (genome position 5,304 in Genbank accession no. M67473.1 sequence numbering), detected by all recombination detection methods in RDP4 (Figure 3). Phylogenetic analysis of the complete VP60 capsid gene revealed that the RHDV G2 viruses clustered with other recent RHDV G2 samples from Australia and that the RHDVa samples clustered (99% nucleotide sequence identity) with exotic RHDVa variants, particularly a variant from China, XA/China/2010 (Genbank accession no. JN165234) (Figure 4, panel A) (14). In contrast, the nonstructural genes of the RHDVa viruses found in Australia clustered with the RCV-A1–like clade that contains several other lagoviruses, including RCV-A1/RHDV2 recombinant viruses (46) and Michigan RCV (47), neither of which have been detected in Australia (Figure 4, panel B). Although the nonstructural genes of these viruses cluster most closely with RCV-A1, they are quite distinct, sharing only 83.3%–84.4% nucleotide identity with RCV-A1 sequences in this region, and could therefore represent a novel lagovirus genotype.


We describe the detection and characterization of RHDVa in Australia, the first of 2 recent incursions of exotic lagoviruses (23). The route of entry of these viruses into Australia is unclear. The spread of RHDVa from Europe to distant locations such as Asia, the Americas, Reunion Island, and Australia (1) highlights the ease with which this virus can be disseminated, despite strict quarantine regulations. Caliciviruses are known to be highly environmentally stable, they replicate to very high titers in infected rabbits, and they are efficiently transmitted (1), suggesting that they would be relatively easy to disseminate inadvertently, analogous to the situation observed with the emergence of canine parvovirus (48).
Unexpectedly, full-genome sequencing revealed that all RHDVa viruses in Australia had recombinant genomes, with the capsid gene closely related to a 2010 virus from China and nonstructural genes related to RCV-A1–like viruses (Figure 4). This finding highlights the importance of recombination in the generation of genetic diversity in this virus family. The emergence of this RCV-A1–like/RHDVa recombinant is of particular interest given the differences in tissue tropisms of RHDVa, which replicates in the liver, and nonpathogenic RCVs, which replicate in the small intestine (1,35). Recombinant viruses containing RCV-A1–like nonstructural genes and RHDV2 structural genes have also been reported (46). Although RHDV is widely disseminated in an infected animal, there is no evidence to suggest it can actively replicate in the small intestine. Conversely, RCV-A1 has only been shown to infect the columnar epithelial cells of the small intestine (35). Therefore, we hypothesize that the recombination event likely occurred during early viral replication in macrophages, which has been suggested for both RHDV (1) and RCV-A1 (T. Strive, unpub. data). Although the exact parental virus and source of the RCV-A1–like nonstructural genes were not identified, the distance between the nonstructural genes of RHDVa in Australia and currently sampled RCV-A1 nonstructural genes (83.3%–84.4% nt identity, Figure 4B) suggests either that the recombination event occurred before the virus arrived in Australia with an as-yet uncharacterized calicivirus or that recombination did occur in Australia but the parental virus belongs to a divergent RCV-A1 lineage that has not yet been identified.
The exotic RHDVa variant was first detected in Australia in samples from an RHD outbreak in December 2013 in domestic rabbits in a suburb of northern Sydney (32). During January 2014–March 2015, this virus was detected in a further 8 outbreaks. Detections were localized to southeastern Australia, predominantly in the Sydney basin. In contrast, the initial spread of the CAPM V-351 virus in 1995 was estimated to be 50 km/wk (1). Seven of the 9 documented RHDVa outbreaks were in domestic rabbits, although this probably reflects sampling bias, because very limited surveillance occurred in NSW and ACT in wild rabbits during the sampling period.
The seroprevalence of RHDV and RCV-A1 antibodies was shown to be very high in wild rabbit populations in southeastern Australia (49), which may limit RHDVa infection and transmission rates in wild rabbits compared with unvaccinated domestic rabbits, due to partial or complete immunological cross-protection from previous infection with RHDV, RCV-A1, or both. In addition, RHDV outbreaks are much more readily detected in domestic rabbit facilities where animals are checked daily, whereas outbreaks in wild rabbits often go unnoticed, particularly as carcasses are rapidly removed by scavengers. Movement of domestic rabbits for breeding and showing purposes may also have facilitated spread of the virus in this population within the Sydney basin. Nevertheless, RHDVa was confirmed in 2 wild rabbit populations, in Girvan, NSW (outbreak 7), and Mulligan’s Flat (outbreak 11; Table 4), suggesting infection and transmission of RHDVa was occurring in wild populations.
In Mulligan’s Flat (outbreak 11), we recovered RHDVa RNA from a healthy shot rabbit following an observed decline in rabbit numbers. Subsequent serologic analyses from a sample of healthy shot rabbits suggests that high levels of cross-reacting isotype antibodies (IgA, IgM, and IgG), in combination with low or negative test results for specific competition ELISAs for RHDV and RCV-A1, may be used as a tool to infer previous exposure of a rabbit population to RHDVa. This particular serologic profile may prove useful for large-scale field epidemiology studies attempting to track the spread and impact of RHDVa in the Australian landscape, even in the absence of a specific serologic test. Furthermore, our study demonstrates that PCR analysis of bile in recovering, recently RHDV-infected (IgM-positive) rabbits can be an additional tool to determine to which virus a wild rabbit population had been exposed.
Several of the RHDVa cases in domestic rabbits were reported to have occurred in previously vaccinated rabbits, although details about the vaccination regimens and intervals were difficult to retrieve. The RHDVa variant JN165234/XA/China/2010, which closely resembles the VP60 sequence of the RHDVa from Australia we describe, was reported to be able to partially overcome protection induced by vaccination against classic RHDV (14). However, the vaccine used in that study was not a commercial preparation, and dose rates were not provided (14). The same study demonstrated that vaccination against classic RHDV induced 100% protection against another RHDVa virus from China (14). RHDV2 is considered to be a second lagovirus serotype, but the OIE Terrestrial Manual does not discriminate between RHDV and RHDVa with respect to vaccination and conservatively recommends that breeder rabbits be vaccinated every 6 months to ensure protection in commercial rabbitries (11). Australian authorities have issued updated vaccination recommendations subsequent to the detection of RHDV2 to reflect those of the OIE (11). Ultimately, a polyvalent vaccine effective against all viruses circulating in Australia (RHDV and RHDV2) would be preferred for domestic rabbits.
Before the detection of RHDVa, the deliberately released CAPM V-351 and its descendants were the only RHDV variants circulating in Australia, despite both RHDVa and RHDV2 circulating overseas since 1997 and 2010, respectively (7,12,18). As such, the RHDV antigen capture ELISA was routinely used for diagnostic purposes. Given that the RHDVa virus we report here shows little to no reactivity in this ELISA, it is recommended that the monoclonal antibodies used in diagnostic ELISAs be updated, as described in the OIE Terrestrial Manual, or that alternative detection methods such as RT-PCR be used for RHD diagnosis in Australia (11). At this stage, it is unclear what effects the presence of this new virus will have on Australia’s wild rabbit populations and on biocontrol efforts using RHDV. With the recent (March 2017) release of another RHDVa virus (an RHDVa from South Korea), it is critical that ongoing monitoring and surveillance systems are in place to explore the interactions of these new viruses as they establish under Australian field conditions (50). In this context, it will be particularly interesting to determine if the recombinant RHDVa described here persists in the face of the exotic incursion of RHDV2, which has reportedly replaced RHDV and RHDVa viruses in many parts of Europe (2426). Once again, the epidemiology of these viruses demonstrates the utility of Australia’s rabbit populations serving as a model system for disease emergence, viral competition, and evolution.
Dr. Mahar is a postdoctoral associate at the University of Sydney and collaborates closely with CSIRO on a project that aims to understand the evolution of virulence in rabbit hemorrhagic disease virus. She has a general interest in molecular microbiology and is particularly interested in virus evolution and virus discovery.


We thank Barry Richardson, Oliver Orgill, and the ACT parks and conservation service rangers of Mulligan’s Flat Woodland Sanctuary for their assistance with sampling.
J.E.M. is supported by the Australian Research Council (grant DP140103362). E.C.H. is supported by an NHMRC Australia Fellowship (GNT1037231).


  1. Abrantes Jvan der Loo WLe Pendu JEsteves PJRabbit haemorrhagic disease (RHD) and rabbit haemorrhagic disease virus (RHDV): a review. Vet Res (Faisalabad)2012;43:12DOIPubMed
  2. Delibes-Mateos MFerreira CCarro FEscudero MAGortázar CEcosystem effects of variant rabbit hemorrhagic disease virus, Iberian Peninsula.Emerg Infect Dis2014;20:21668DOIPubMed
  3. Cooke BChudleigh PSimpson SSaunders GThe economic benefits of the biological control of rabbits in Australia, 1950–2011. Aust Econ Hist Rev2013;53:91107DOI
  4. Pedler RDBrandle RRead JLSouthgate RBird PMoseby KERabbit biocontrol and landscape-scale recovery of threatened desert mammals.Conserv Biol2016;30:77482DOIPubMed
  5. Cooke BDFenner FRabbit haemorrhagic disease and the biological control of wild rabbits, Oryctolagus cuniculus, in Australia and New Zealand.Wildl Res2002;29:689706DOI
  6. O’Hara PThe illegal introduction of rabbit haemorrhagic disease virus in New Zealand. Rev Sci Tech2006;25:11923DOIPubMed
  7. Eden J-SKovaliski JDuckworth JASwain GMahar JEStrive Tet al. Comparative phylodynamics of rabbit hemorrhagic disease virus in Australia and New Zealand. J Virol2015;89:954858DOIPubMed
  8. Wirblich CThiel HJMeyers GGenetic map of the calicivirus rabbit hemorrhagic disease virus as deduced from in vitro translation studies. J Virol1996;70:797483.PubMed
  9. Meyers GWirblich CThiel HJGenomic and subgenomic RNAs of rabbit hemorrhagic disease virus are both protein-linked and packaged into particles. Virology1991;184:67786DOIPubMed
  10. Le Gall-Reculé GZwingelstein FLaurent Sde Boisséson CPortejoie YRasschaert DPhylogenetic analysis of rabbit haemorrhagic disease virus in France between 1993 and 2000, and the characterisation of RHDV antigenic variants. Arch Virol2003;148:6581DOIPubMed
  11. Lavazza ACapucci L. Chapter 2.6.2. Rabbit haemorrhagic disease. In: OIE manual of diagnostic tests and vaccines for terrestrial animals. 2016 [cited 2017 Jan 13]. http://www.oie.int/fileadmin/Home/eng/Health_standards/tahm/2.06.02_RHD.pdf
  12. Capucci LFallacara FGrazioli SLavazza APacciarini MLBrocchi EA further step in the evolution of rabbit hemorrhagic disease virus: the appearance of the first consistent antigenic variant. Virus Res1998;58:11526DOIPubMed
  13. Schirrmeier HReimann IKöllner BGranzow HPathogenic, antigenic and molecular properties of rabbit haemorrhagic disease virus (RHDV) isolated from vaccinated rabbits: detection and characterization of antigenic variants. Arch Virol1999;144:71935DOIPubMed
  14. Wang XHao HQiu LDang RDu EZhang Set al. Phylogenetic analysis of rabbit hemorrhagic disease virus in China and the antigenic variation of new strains. Arch Virol2012;157:152330DOIPubMed
  15. Oem JKLee KNRoh ISLee KKKim SHKim HRet al. Identification and characterization of rabbit hemorrhagic disease virus genetic variants isolated in Korea. J Vet Med Sci2009;71:151923DOIPubMed
  16. Burmakina GMalogolovkina NLunitsin ATitov ITsybanov SMalogolovkin AComparative analysis of rabbit hemorrhagic disease virus strains originating from outbreaks in the Russian Federation. Arch Virol2016;161:19739DOIPubMed
  17. Abrantes JLopes AMDalton KPParra FEsteves PJDetection of RHDVa on the Iberian Peninsula: isolation of an RHDVa strain from a Spanish rabbitry. Arch Virol2014;159:3216DOIPubMed
  18. Le Gall-Reculé GLavazza AMarchandeau SBertagnoli SZwingelstein FCavadini Pet al. Emergence of a new lagovirus related to Rabbit Haemorrhagic Disease Virus. Vet Res (Faisalabad)2013;44:81DOIPubMed
  19. Dalton KPNicieza IBalseiro AMuguerza MARosell JMCasais Ret al. Variant rabbit hemorrhagic disease virus in young rabbits, Spain. Emerg Infect Dis2012;18:200912DOIPubMed
  20. Abrantes JLopes AMDalton KPMelo PCorreia JJRamada Met al. New variant of rabbit hemorrhagic disease virus, Portugal, 2012-2013. Emerg Infect Dis2013;19:19002DOIPubMed
  21. Baily JLDagleish MPGraham MMaley MRocchi MSRHDV variant 2 presence detected in Scotland. Vet Rec2014;174:411DOIPubMed
  22. Duarte MHenriques MBarros SCFagulha TRamos FLuís Tet al. Detection of RHDV variant 2 in the Azores. Vet Rec2015;176:130DOIPubMed
  23. Hall RNMahar JEHaboury SStevens VHolmes ECStrive TEmerging rabbit hemorrhagic disease virus 2 (RHDVb), Australia. Emerg Infect Dis2015;21:22768DOIPubMed
  24. Dalton KPNicieza IAbrantes JEsteves PJParra FSpread of new variant RHDV in domestic rabbits on the Iberian Peninsula. Vet Microbiol2014;169:6773DOIPubMed
  25. Lopes AMCorreia JAbrantes JMelo PRamada MMagalhães MJet al. Is the new variant RHDV replacing genogroup 1 in Portuguese wild rabbit populations? Viruses2014;7:2736DOIPubMed
  26. Westcott DGFrossard JPEverest DDastjerdi ADuff JPSteinbach Fet al. Incursion of RHDV2-like variant in Great Britain. Vet Rec2014;174:333DOIPubMed
  27. Puggioni GCavadini PMaestrale CScivoli RBotti GLigios Cet al. The new French 2010 Rabbit Hemorrhagic Disease Virus causes an RHD-like disease in the Sardinian Cape hare (Lepus capensis mediterraneus). Vet Res (Faisalabad)2013;44:96DOIPubMed
  28. Camarda APugliese NCavadini PCircella ECapucci LCaroli Aet al. Detection of the new emerging rabbit haemorrhagic disease type 2 virus (RHDV2) in Sicily from rabbit (Oryctolagus cuniculus) and Italian hare (Lepus corsicanus). Res Vet Sci2014;97:6425DOIPubMed
  29. Lopes AMMarques SSilva EMagalhães MJPinheiro AAlves PCet al. Detection of RHDV strains in the Iberian hare (Lepus granatensis): earliest evidence of rabbit lagovirus cross-species infection. Vet Res2014;45:94.PubMed
  30. Hall RNPeacock DEKovaliski JMahar JEMourant RPiper Met al. Detection of RHDV2 in European brown hares (Lepus europaeus) in Australia.Vet Rec2017;180:121DOIPubMed
  31. Mahar JENicholson LEden JSDuchêne SKerr PJDuckworth Jet al. Benign rabbit caliciviruses exhibit evolutionary dynamics similar to those of their virulent relatives. J Virol2016;90:931729DOIPubMed
  32. Office International des Epizooties. Rabbit haemorrhagic disease, Australia—immediate notification report. Ref OIE = 14719. 2014 [cited 2017 Jan 16]. http://www.oie.int/wahis_2/public/wahid.php/Reviewreport/Review?page_refer=MapFullEventReport&reportid=14719
  33. Collins BJWhite JRLenghaus CMorrissy CJWestbury HAPresence of rabbit haemorrhagic disease virus antigen in rabbit tissues as revealed by a monoclonal antibody dependent capture ELISA. J Virol Methods1996;58:14554DOIPubMed
  34. Gall AHoffmann BTeifke JPLange BSchirrmeier HPersistence of viral RNA in rabbits which overcome an experimental RHDV infection detected by a highly sensitive multiplex real-time RT-PCR. Vet Microbiol2007;120:1732DOIPubMed
  35. Strive TWright JDRobinson AJIdentification and partial characterisation of a new Lagovirus in Australian wild rabbits. Virology2009;384:97105DOIPubMed
  36. Jahnke MHolmes ECKerr PJWright JDStrive TEvolution and phylogeography of the nonpathogenic calicivirus RCV-A1 in wild rabbits in Australia. J Virol2010;84:12397404DOIPubMed
  37. Elsworth PCooke BDKovaliski JSinclair RHolmes ECStrive TIncreased virulence of rabbit haemorrhagic disease virus associated with genetic resistance in wild Australian rabbits (Oryctolagus cuniculus). Virology2014;464-465:41523DOIPubMed
  38. Hall RNCapucci LMatthaei MEsposito SKerr PJFrese Met al. An in vivo system for directed experimental evolution of rabbit haemorrhagic disease virus. PLoS One2017;12:e0173727DOIPubMed
  39. Kearse MMoir RWilson AStones-Havas SCheung MSturrock Set al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics2012;28:16479DOIPubMed
  40. Guindon SGascuel ORannala BA simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol2003;52:696704DOIPubMed
  41. Posada DjModelTest: phylogenetic model averaging. Mol Biol Evol2008;25:12536DOIPubMed
  42. Martin DPMurrell BGolden MKhoosal AMuhire B. RDP4: Detection and analysis of recombination patterns in virus genomes. Virus Evol. 2015;1:vev003 .
  43. Liu JKerr PJStrive TA sensitive and specific blocking ELISA for the detection of rabbit calicivirus RCV-A1 antibodies. Virol J2012;9:182DOIPubMed
  44. Liu JKerr PJWright JDStrive TSerological assays to discriminate rabbit haemorrhagic disease virus from Australian non-pathogenic rabbit calicivirus. Vet Microbiol2012;157:34554DOIPubMed
  45. Cooke BDRobinson AJMerchant JCNardin ACapucci LUse of ELISAs in field studies of rabbit haemorrhagic disease (RHD) in Australia.Epidemiol Infect2000;124:56376DOIPubMed
  46. Lopes AMDalton KPMagalhães MJParra FEsteves PJHolmes ECet al. Full genomic analysis of new variant rabbit hemorrhagic disease virus revealed multiple recombination events. J Gen Virol2015;96:130919DOIPubMed
  47. Bergin ILWise AGBolin SRMullaney TPKiupel MMaes RKNovel calicivirus identified in rabbits, Michigan, USA. Emerg Infect Dis2009;15:195562DOIPubMed
  48. MacLachlan NJ. Parvoviridae A2. In: Dubovi EJ, editor. Fenner’s veterinary virology. 5th edition. Boston: Academic Press; 2017. p. 245–57.
  49. Cox TELiu Jde Ven RVStrive TDifferent serological profiles to co-occurring pathogenic and non-pathogenic caliciviruses in wild European rabbits (Oryctolagus cuniculus) across Australia. J Wildl Dis2017;53:47281DOIPubMed
  50. Matthaei MKerr PJRead AJHick PHaboury SWright JDet al. Comparative quantitative monitoring of rabbit haemorrhagic disease viruses in rabbit kittens. Virol J2014;11:109DOIPubMed



Cite This Article

DOI: 10.3201/eid2401.170412
1Current affiliation: University of Alabama at Birmingham School of Medicine, Birmingham, Alabama, USA.

No hay comentarios:

Publicar un comentario