Distribution and prevalence of the Australian non-pathogenic rabbit calicivirus is correlated with rainfall and temperature

First Detection of Benign Rabbit Caliciviruses in Chile

Pathogenic lagoviruses (Rabbit hemorrhagic disease virus, RHDV) are widely spread across the world and are used in Australia and New Zealand to control populations of feral European rabbits. The spread of the non-pathogenic lagoviruses, e.g., rabbit calicivirus (RCV), is less well studied as the infection results in no clinical signs. Nonetheless, RCV has important implications for the spread of RHDV and rabbit biocontrol as it can provide varying levels of cross-protection against fatal infection with pathogenic lagoviruses. In Chile, where European rabbits are also an introduced species, myxoma virus was used for localised biocontrol of rabbits in the 1950s. To date, there have been no studies investigating the presence of lagoviruses in the Chilean feral rabbit population. In this study, liver and duodenum rabbit samples from central Chile were tested for the presence of lagoviruses and positive samples were subject to whole RNA sequencing and subsequent data analysis. Phylogenetic analysis revealed a novel RCV variant in duodenal samples that likely originated from European RCVs. Sequencing analysis also detected the presence of a rabbit astrovirus in one of the lagovirus-positive samples.

Sustained Impact of RHDV2 on Wild Rabbit Populations across Australia Eight Years after Its Initial Detection

Following the arrival of rabbit haemorrhagic disease virus 2 (RHDV2) in Australia, average rabbit population abundances were reduced by 60% between 2014 and 2018 based on monitoring data acquired from 18 sites across Australia. During this period, as the seropositivity to RHDV2 increased, concurrent decreases were observed in the seroprevalence of both the previously circulating RHDV1 and RCVA, a benign endemic rabbit calicivirus. However, the detection of substantial RHDV1 seropositivity in juvenile rabbits suggested that infections were continuing to occur, ruling out the rapid extinction of this variant. Here we investigate whether the co-circulation of two pathogenic RHDV variants was sustained after 2018 and whether the initially observed impact on rabbit abundance was still maintained. We monitored rabbit abundance and seropositivity to RHDV2, RHDV1 and RCVA at six of the initial eighteen sites until the summer of 2022. We observed sustained suppression of rabbit abundance at five of the six sites, with the average population reduction across all six sites being 64%. Across all sites, average RHDV2 seroprevalence remained high, reaching 60-70% in adult rabbits and 30-40% in juvenile rabbits. In contrast, average RHDV1 seroprevalence declined to <3% in adult rabbits and 5-6% in juvenile rabbits. Although seropositivity continued to be detected in a low number of juvenile rabbits, it is unlikely that RHDV1 strains now play a major role in the regulation of rabbit abundance. In contrast, RCVA seropositivity appears to be reaching an equilibrium with that of RHDV2, with RCVA seroprevalence in the preceding quarter having a strong negative effect on RHDV2 seroprevalence and vice versa, suggesting ongoing co-circulation of these variants. These findings highlight the complex interactions between different calicivirus variants in free-living rabbit populations and demonstrate the changes in interactions over the course of the RHDV2 epizootic as it has moved towards endemicity. While it is encouraging from an Australian perspective to see sustained suppression of rabbit populations in the eight years following the arrival of RHDV2, it is likely that rabbit populations will eventually recover, as has been observed with previous rabbit pathogens.

Structural Basis for Rabbit Hemorrhagic Disease Virus Antibody Specificity

Isolated RHDV antibodies have been used for decades to distinguish between antigenic variants, monitor temporal capsid evolution, and examine neutralizing capacities. In this study, we provided the structural basis for an RHDV GI.2 specific diagnostic antibody (2D9) binding and reveal that a small number of amino acid substitutions at the binding site could differentiate between RHDV GI.2 and GI.1b.

Cross-protection, infection, and case fatality rates in wild European rabbits experimentally challenged with different rabbit haemorrhagic disease viruses

Rabbit haemorrhagic disease virus 2 (RHDV2) is now the dominant calicivirus circulating in wild rabbit populations in Australia. This study compared the infection and case fatality rates of RHDV2 and two RHDVs in wild rabbits, as well as their ability to overcome immunity to the respective other strains. Wild rabbits were allocated to groups either blindly or based on prescreening for RHDV/RHDV2 antibodies at capture. Rabbits were monitored regularly until their death or humane killing at 7 days post infection. Liver and eyeball samples were collected for lagovirus testing and aging rabbits, respectively. At capture, rabbits showed high seroprevalence to RHDV2 but not to RHDV. In RHDV/RHDV2 seronegative rabbits at capture, infection rates were highest in those inoculated with RHDV2 (81.8%, 18/22), followed by K5 (53.8%, 7/13) and CZECH (40.0%, 2/5), but these differences were not statistically significant. In rabbits with previous exposure to RHDV2 at capture, infection rates were highest when inoculated with K5 (59.6%, 31/52) followed by CZECH (46.0%, 23/50), with infection rates higher in younger rabbits for both viruses. In RHDV/RHDV2 seronegative rabbits at capture, case fatality rates were highest for those inoculated with K5 (71.4%), followed by RHDV2 (50.0%) and CZECH (50.0%). In rabbits with previous exposure to RHDV2 at capture, case fatality rates were highest in rabbits inoculated with K5 (12.9%) followed by CZECH (8.7%), with no case fatalities following RHDV2 inoculation. Case fatality rates did not differ significantly between inoculums in either serostatus group at capture. Based on multivariable modelling, time to death post RHDV inoculation increased in rabbits with recent RHDV2 exposure compared to seronegative rabbits and with age. The results suggest that RHDV2 may cause higher mortalities than other variants in seronegative rabbit populations but that K5 may be more effective in reducing rabbit populations in an RHDV2-dominant landscape. This article is protected by copyright. All rights reserved.

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

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

Viruses for Landscape-Scale Therapy: Biological Control of Rabbits in Australia

Viral diseases, whether of animals or humans, are normally considered as problems to be managed. However, in Australia, two viruses have been used as landscape-scale therapeutics to control European rabbits (Oryctolagus cuniculus), the preeminent invasive vertebrate pest species. Rabbits have caused major environmental and agricultural losses and contributed to extinction of native species. It was not until the introduction of Myxoma virus that effective control of this pest was obtained at a continental scale. Subsequent coevolution of rabbit and virus saw a gradual reduction in the effectiveness of biological control that was partially ameliorated by the introduction of the European rabbit flea to act as an additional vector for the virus. In 1995, a completely different virus, Rabbit hemorrhagic disease virus (RHDV), escaped from testing and spread through the Australian rabbit population and again significantly reduced rabbit numbers and environmental impacts. The evolutionary pressures on this virus appear to be producing quite different outcomes to those that occurred with myxoma virus and the emergence and invasion of a novel genotype of RHDV in 2014 have further augmented control. Molecular studies on myxoma virus have demonstrated multiple proteins that manipulate the host innate and adaptive immune response; however the molecular basis of virus attenuation and reversion to virulence are not yet understood.

Changes in virus transmission dynamics following the emergence of RHDV2 shed light on its competitive advantage over previously circulating variants

Rabbit haemorrhagic disease virus (RHDV) is highly pathogenic to European rabbits. Until recently, only one serotype of RHDV was known, GI.1/RHDV. RHDV2/GI.2 is a novel virus that has rapidly spread and become the dominant pathogenic calicivirus in wild rabbits worldwide. It is speculated that RHDV2 has three competitive advantages over RHDV: (a) the ability to partially overcome immunity to other variants; (b) the ability to clinically infect young rabbits; and (c) a wider host range. These differences would be expected to influence virus transmission dynamics. We used markers of recent infection (IgM/IgA antibodies) to investigate virus transmission dynamics pre and post the arrival of RHDV2. Our data set contained over 3,900 rabbits sampled across a 7-year period at 12 Australian sites. Following the arrival of RHDV2, seasonal peaks in IgM and IgA seropositivity shifted forward one season, from winter to autumn and spring to winter, respectively. Contrary to predictions, we found only weak effects of rabbit age, seropositivity to non-pathogenic calicivirus RCV-A1 and population abundance on IgM/IgA seropositivity. Our results demonstrate that RHDV2 enters rabbit populations shortly after the commencement of annual breeding cycles. Upon entering, the population RHDV2 undergoes extensive replication in young rabbits, causing clinical disease, high virus shedding, mortality and the creation of virus-laden carcasses. This results in high virus contamination in the environment, furthering the transmission of RHDV2 and initiating outbreaks, whilst simultaneously removing the susceptible cohort required for the effective transmission of RHDV. Although RHDV may enter the population at the same time point, it is sub-clinical in young rabbits, causing minimal virus shedding and low environmental contamination. Our results demonstrate a major shift in epidemiological patterns in virus transmission, providing the first evidence that RHDV2's ability to clinically infect young rabbits is a key competitive advantage in the field.

Field validation of phylodynamic analytical methods for inference on epidemiological processes in wildlife

Amongst newly developed approaches to analyse molecular data, phylodynamic models are receiving much attention because of their potential to reveal changes to viral populations over short periods. This knowledge can be very important for understanding disease impacts. However, their accuracy needs to be fully understood, especially in relation to wildlife disease epidemiology, where sampling and prior knowledge may be limited. The release of the rabbit haemorrhagic disease virus (RHDV) as biological control in naïve rabbit populations in Australia in 1996 provides a unique data set with which to validate phylodynamic models. By comparing results obtained from RHDV sequence data with our current understanding of RHDV epidemiology in Australia, we evaluated the performances of these recently developed models. In line with our expectations, coalescent analyses detected a sharp increase in the virus population size in the first few months after release, followed by a more gradual increase. Phylodynamic analyses using a birth-death model generated effective reproductive number estimates (the average number of secondary infections per each infectious case, R_e ) larger than one for most of the epochs considered. However, the possible range of the initial R_e included estimates lower than one despite the known rapid spread of RHDV in Australia. Furthermore, the analyses that accounted for geographical structuring failed to converge. We argue that the difficulties that we encountered most likely stem from the fact that the samples available from 1996 to 2014 were too sparse with respect to both geographic and within outbreak coverage to adequately infer some of the model parameters. In general, while these phylodynamic analyses proved to be greatly informative in some regards, we caution that their interpretation may not be straightforward. We recommend further research to evaluate the robustness of these models to assumption violations and sensitivity to sampling regimes.

Models of spatiotemporal variation in rabbit abundance reveal management hot spots for an invasive species

The European rabbit (Oryctolagus cuniculus) is a notorious economic and environmental pest species in its invasive range. To better understand the population and range dynamics of this species, 41 yr of abundance data have been collected from 116 unique sites across a broad range of climatic and environmental conditions in Australia. We analyzed this time series of abundance data to determine whether interannual variation in climatic conditions can be used to map historic, contemporary, and potential future fluctuations in rabbit abundance from regional to continental scales. We constructed a hierarchical Bayesian regression model of relative abundance that corrected for observation error and seasonal biases. The corrected abundances were regressed against environmental and disease variables in order to project high spatiotemporal resolution, continent-wide rabbit abundances. We show that rabbit abundance in Australia is highly variable in space and time, being driven primarily by internnual variation in temperature and precipitation in concert with the prevalence of a non-pathogenic virus. Moreover, we show that internnual variation in local spatial abundances can be mapped effectively at a continental scale using highly resolved spatiotemporal predictors, allowing "hot spots" of persistently high rabbit abundance to be identified. Importantly, cross-validated model performance was fair to excellent within and across distinct climate zones. Long-term monitoring data for invasive species can be used to map fine-scale spatiotemporal fluctuations in abundance patterns when accurately accounting for inherent sampling biases. Our analysis provides ecologists and pest managers with a clearer understanding of the determinants of rabbit abundance in Australia, offering an important new approach for predicting spatial abundance patterns of invasive species at the near-term temporal scales that are directly relevant to resource management.

Retrospective serological analysis reveals presence of the emerging lagovirus RHDV2 in Australia in wild rabbits at least five months prior to its first detection

The lagovirus rabbit haemorrhagic disease virus (RHDV) has been circulating in Australia since the mid-1990s when it was released to control overabundant rabbit populations. In recent years, the viral diversity of different RHDVs in Australia has increased, and currently four different types of RHDV are known to be circulating. To allow for ongoing epidemiological studies and impact assessments of these viruses on Australian wild rabbit populations, it is essential that serological tools are updated. To this end, reference sera were produced against all four virulent RHDVs (RHDV, RHDV2 and two different strains of RHDVa) known to be present in Australia and tested in a series of available immunological assays originally developed for the prototype RHDV, to assess patterns of cross-reactivity and the usefulness of these assays to detect lagovirus antibodies, either in a generic or specific manner. Enzyme-linked immunosorbent assays (ELISAs) developed to detect antibody isotypes IgM, IgA and IgG were sufficiently cross-reactive to detect antibodies raised against all four virulent lagoviruses. For the more specific detection of antibodies to the antigenically more different RHDV2, a competition ELISA was adapted using RHDV2-specific monoclonal antibodies in combination with Australian viral antigen. Archival serum banks from a long-term rabbit monitoring site where rabbits were sampled quarterly over a period of 6 years were re-screened using this assay and revealed serological evidence for the arrival of RHDV2 in this population at least 5 months prior to its initial detection in Australia in a dead rabbit in May 2015. The serological methods and reference reagents described here will provide valuable tools to study presence, prevalence and impact of RHDV2 on Australian rabbit populations; however, the discrimination of different antigenic variants of RHDVs as well as mixed infections at the serological level remains challenging.

The impact of RHDV-K5 on rabbit populations in Australia: an evaluation of citizen science surveys to monitor rabbit abundance

The increasing popularity of citizen science in ecological research has created opportunities for data collection from large teams of observers that are widely dispersed. We established a citizen science program to complement the release of a new variant of the rabbit biological control agent, rabbit haemorrhagic disease virus (RHDV), known colloquially as K5, across Australia. We evaluated the impact of K5 on the national rabbit population and compared citizen science and professionally-collected spotlight count data. Of the citizen science sites (n = 219), 93% indicated a decrease in rabbit abundance following the release of K5. The overall finite monthly growth rate in rabbit abundance was estimated as 0.66 (95%CI, 0.26, 1.03), averaging a monthly reduction of 34% at the citizen science sites one month after the release. No such declines were observed at the professionally monitored sites (n = 22). The citizen science data submissions may have been unconsciously biased or the number of professional sites may have been insufficient to detect a change. Citizen science participation also declined by 56% over the post-release period. Future programs should ensure the use of blinded trials to check for unconscious bias and consider how incentives and/or the good will of the participants can be maintained throughout the program.

Myxomatosis and Rabbit Haemorrhagic Disease: A 30-Year Study of the Occurrence on Commercial Farms in Spain

In this retrospective study, we describe the relative occurrence of clinical myxomatosis, and rabbit haemorrhagic disease (RHD), on 1714 commercial farms visited in Spain, between 1988 and 2018. We determined the annual prevalence based on 817 visits to 394 farms affected by myxomatosis. Myxomatosis was more prevalent from August to March, being lowest in June (3%) and highest in September (8.9%). With regard to RHD, we assessed 253 visits to 156 affected farms. We analyzed mean annual and monthly incidence. Two important RHD epidemics occurred; the first in 1988-1989 due to RHDV GI.1 (also known as RHDV), and the second from 2011 to 2013 due to RHDV GI.2 (RHDV2 or RHDVb). These epidemics occurred at times when effective vaccination had not been carried out. Relative monthly incidence in 2011-2018 was higher from April to August (p < 0.001). The results we obtained from 1404 necropsies on 102 farms did not clearly relate serosanguinous nasal discharge in rabbits with disease caused by GI.2 infection. We also assessed vaccination schedules used on 200 doe farms visited from the end of 2014 to 2018; 95.5% vaccinated against myxomatosis and 97.5% against RHD. Both diseases remain prevalent; however, effective vaccination has produced a steady decline in myxomatosis and RHDV GI.1 and GI.2 on-farm detection. The maintenance of high hygienic standards will be needed to continue and improve this control. However, further studies are required to investigate the causes of sustained virus presence and vaccine breaks.

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

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

Variable changes in nematode infection prevalence and intensity after Rabbit Haemorrhagic Disease Virus emerged in wild rabbits in Scotland and New Zealand

The myxoma virus (a microparasite) reduced wild rabbit numbers worldwide when introduced in the 1950s, and is known to interact with co-infecting helminths (macroparasites) causing both increases and decreases in macroparasite population size. In the 1990s Rabbit Haemorrhagic Disease Virus (RHDV) infected rabbits and also significantly reduced rabbit numbers in several countries. However, not much is known about RHDV interactions with macroparasites. In this study, we compare prevalence and intensity of infection for three gastrointestinal nematode species (Trichostrongylus retortaeformis, Graphidium strigosum and Passalurus ambiguus) before and after RHDV spread across host populations in Scotland and New Zealand. During one common season, autumn, prevalence of T. retortaeformis was higher after RHDV spread in both locations, whereas it was lower for G. strigosum and P. ambiguus after RHDV arrived in New Zealand, but higher in Scotland. Meanwhile, intensity of infection for all species decreased after RHDV arrived in New Zealand, but increased in Scotland. The impact of RHDV on worm infections was generally similar across seasons in Scotland, and also similarities in seasonality between locations suggested effects on infection patterns in one season are likely similar year-round. The variable response by macroparasites to the arrival of a microparasite into Scottish and New Zealand rabbits may be due to differences in the environment they inhabit, in existing parasite community structure, and to some extent, in the relative magnitude of indirect effects. Specifically, our data suggest that bottom-up processes after the introduction of a more virulent strain of RHDV to New Zealand may affect macroparasite co-infections by reducing the availability of their shared common resource, the rabbits. Clearly, interactions between co-infecting micro- and macroparasites vary in host populations with different ecologies, and significantly impact parasite community structure in wildlife.

•

Nematode communities in Scotland and New Zealand were compared pre and post Rabbit Haemorrhagic Disease Virus introduction.

•

Similar species occur in both rabbit populations, but prevalence and intensity changed in opposing directions after RHDV.

•

RHDV had a major impact on rabbit populations, and our data show differing impacts on macroparasites in the two countries.

•

Variability in rabbit environment, parasite community structure, and indirect interaction processes may explain differences.

•

Results can help understand interactions between co-infecting parasites and their epidemiology in wild and domestic animals.

Benign Rabbit Calicivirus in New Zealand

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

DIFFERENT SEROLOGICAL PROFILES TO CO-OCCURRING PATHOGENIC AND NONPATHOGENIC CALICIVIRUSES IN WILD EUROPEAN RABBITS (ORYCTOLAGUS CUNICULUS) ACROSS AUSTRALIA

Rabbit hemorrhagic disease virus (RHDV) was released in Australia as a biocontrol agent for wild European rabbits ( Oryctolagus cuniculus ) in 1995-96; however, its effects were variable across Australia with the greatest population reductions seen in lower annual rainfall areas (<400 mm). There is speculation that the reduced effectiveness observed at higher annual rainfall sites is at least partially due to the presence of a nonpathogenic calicivirus (RCV-A1). The RCV-A1 is related to RHDV and confers partial and transient protection against lethal RHDV infection in laboratory tests. What is not well understood is where, how, and to what degree RCV-A1 impedes the effect of RHDV-mediated rabbit control under field conditions. We investigated seven wild rabbit populations across six states and territories representing different seasonal rainfall zones across Australia, four times during 2011-12, to investigate if the presence and prevalence of RCV-A1 coincided with a change in RHDV immunity status within these populations. Besides serology, tissue samples from both trapped and shot rabbits were collected for virus detection by reverse transcription PCR. Overall, 52% (n=258) of the total samples (n=496) tested positive for RHDV antibodies and 42% (n=208) positive for RCV-A1 antibodies; 30% (n=150) of the sera contained antibodies to both viruses. The proportion of rabbits with RHDV antibodies increased significantly at sites where RCV-A1 antibodies were present (χ²₁, α=0.1, P<0.001). Evidence that preinfection of RCV-A1 may lead to a higher proportion of sampled rabbits with antibodies to both viruses was found at only one site.

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

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

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

Two closely related caliciviruses cocirculate in Australia: rabbit hemorrhagic disease virus (RHDV) and rabbit calicivirus Australia 1 (RCV-A1). RCV-A1 causes benign enteric infections in the European rabbit (Oryctolagus cuniculus) in Australia and New Zealand, while its close relative RHDV causes a highly pathogenic infection of the liver in the same host. The comparison of these viruses provides important information on the nature and trajectory of virulence evolution, particularly as highly virulent strains of RHDV may have evolved from nonpathogenic ancestors such as RCV-A1. To determine the evolution of RCV-A1 we sequenced the full-length genomes of 44 RCV-A1 samples isolated from healthy rabbits and compared key evolutionary parameters to those of its virulent relative, RHDV. Despite their marked differences in pathogenicity and tissue tropism, RCV-A1 and RHDV have evolved in a very similar manner. Both viruses have evolved at broadly similar rates, suggesting that their dynamics are largely shaped by high background mutation rates, and both exhibit occasional recombination and an evolutionary environment dominated by purifying selection. In addition, our comparative analysis revealed that there have been multiple changes in both virulence and tissue tropism in the evolutionary history of these and related viruses. Finally, these new genomic data suggest that either RCV-A1 was introduced into Australia after the introduction of myxoma virus as a biocontrol agent in 1950 or there was drastic reduction of the rabbit population, and hence of RCV-A1 genetic diversity, perhaps coincident with the emergence of myxoma virus.

IMPORTANCE

The comparison of closely related viruses that differ profoundly in propensity to cause disease in their hosts offers a powerful opportunity to reveal the causes of changes in virulence and to study how such changes alter the evolutionary dynamics of these pathogens. Here we describe such a novel comparison involving two closely related RNA viruses that cocirculate in Australia, the highly virulent rabbit hemorrhagic disease virus (RHDV) and the nonpathogenic rabbit calicivirus Australia 1 (RCV-A1). Both viruses infect the European rabbit, but they differ in virulence, tissue tropism, and mechanisms of transmission. Surprisingly, and despite these fundamental differences, RCV-A1 and RHDV have evolved at very similar (high) rates and with strong purifying selection. Furthermore, candidate key mutations were identified that may play a role in virulence and/or tissue tropism and therefore warrant further investigation.

Comparative Phylodynamics of Rabbit Hemorrhagic Disease Virus in Australia and New Zealand

The introduction of rabbit hemorrhagic disease virus (RHDV) into Australia and New Zealand during the 1990s as a means of controlling feral rabbits is an important case study in viral emergence. Both epidemics are exceptional in that the founder viruses share an origin and the timing of their release is known, providing a unique opportunity to compare the evolution of a single virus in distinct naive populations. We examined the evolution and spread of RHDV in Australia and New Zealand through a genome-wide evolutionary analysis, including data from 28 newly sequenced RHDV field isolates. Following the release of the Australian inoculum strain into New Zealand, no subsequent mixing of the populations occurred, with viruses from both countries forming distinct groups. Strikingly, the rate of evolution in the capsid gene was higher in the Australian viruses than in those from New Zealand, most likely due to the presence of transient deleterious mutations in the former. However, estimates of both substitution rates and population dynamics were strongly sample dependent, such that small changes in sample composition had an important impact on evolutionary parameters. Phylogeographic analysis revealed a clear spatial structure in the Australian RHDV strains, with a major division between those viruses from western and eastern states. Importantly, RHDV sequences from the state where the virus was first released, South Australia, had the greatest diversity and were diffuse throughout both geographic lineages, such that this region was likely a source population for the subsequent spread of the virus across the country.

IMPORTANCE

Most studies of viral emergence lack detailed knowledge about which strains were founders for the outbreak or when these events occurred. Hence, the human-mediated introduction of rabbit hemorrhagic disease virus (RHDV) into Australia and New Zealand from known starting stocks provides a unique opportunity to understand viral evolution and emergence. Within Australia, we revealed a major phylogenetic division between viruses sampled from the east and west of the country, with both regions likely seeded by viruses from South Australia. Despite their common origins, marked differences in evolutionary rates were observed between the Australian and New Zealand RHDV, which led to conflicting conclusions about population growth rates. An analysis of mutational patterns suggested that evolutionary rates have been elevated in the Australian viruses, at least in part due to the presence of low-fitness (deleterious) variants that have yet to be selectively purged.