The effect of rabbit population control programmes on the impact of rabbit haemorrhagic disease in south-eastern Australia

Utilizing Molecular Epidemiology and Citizen Science for the Surveillance of Lagoviruses in Australia

Australia has multiple lagoviruses with differing pathogenicity. The circulation of these viruses was traditionally determined through opportunistic sampling events. In the lead up to the nationwide release of RHDVa-K5 (GI.1aP-GI.1a) in 2017, an existing citizen science program, RabbitScan, was augmented to allow members of the public to submit samples collected from dead leporids for lagovirus testing. This study describes the information obtained from the increased number of leporid samples received between 2015 and 2022 and focuses on the recent epidemiological interactions and evolutionary trajectory of circulating lagoviruses in Australia between October 2020 and December 2022. A total of 2771 samples were tested from January 2015 to December 2022, of which 1643 were lagovirus-positive. Notable changes in the distribution of lagovirus variants were observed, predominantly in Western Australia, where RHDV2-4c (GI.4cP-GI.2) was detected again in 2021 after initially being reported to be present in 2018. Interestingly, we found evidence that the deliberately released RHDVa-K5 was able to establish and circulate in wild rabbit populations in WA. Overall, the incorporation of citizen science approaches proved to be a cost-efficient method to increase the sampling area and enable an in-depth analysis of lagovirus distribution, genetic diversity, and interactions. The maintenance of such programs is essential to enable continued investigations of the critical parameters affecting the biocontrol of feral rabbit populations in Australia, as well as to enable the detection of any potential future incursions.

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.

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.

Climate as a control on foredune mode in Southern Australia

Foredunes are formed by aeolian sand deposition in vegetation on the backshore of beaches. In this paper, the foredune mode (nebkha, discontinuous foredune, and continuous foredune), and transgressive dunefield development is studied along the Great Australian Bight (GAB), 2668 km of coastline. Orthophotos are used to classify the foredune mode, coastal landforms and the vegetation, through geographic information systems (GIS), with fieldwork support. The results show that the foredune mode is strongly controlled by rainfall and temperature with respect to latitude, and to drift potential with respect to longitude across the GAB. Between 200 and 300 mm annual rainfall, nebkha predominate. When the annual rainfall is between 300 and 400, at latitude 32°, a clear pattern is not observed in foredune mode and this is identified as a transition zone. Discontinuous foredunes and continuous foredunes are strongly represented in regions experiencing above 400 mm annual rainfall. The main contribution of this study is the identification of foredune modes which are not only related to a climatic gradient and latitude, but also related to variations in longitude, vegetation cover and diversity, and dune mobility indices. Finally, there are other environmental relationships between the wind and longitude, where the geomorphology of the bay could be playing an important role.

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.

Effects of myxoma virus and rabbit hemorrhagic disease virus on the physiological condition of wild European rabbits: Is blood biochemistry a useful monitoring tool?

Myxomatosis and rabbit hemorrhagic disease (RHD) are the major viral diseases that affect the wild European rabbit (Oryctolagus cuniculus). These diseases arrived in Europe within the last decades and have caused wild rabbit populations to decline dramatically. Both viruses are currently considered to be endemic in the Iberian Peninsula; periodic outbreaks that strongly impact wild populations regularly occur. Myxoma virus (MV) and rabbit hemorrhagic disease virus (RHDV) alter the physiology of infected rabbits, resulting in physical deterioration. Consequently, the persistence and viability of natural populations are affected. The main goal of our study was to determine if blood biochemistry is correlated with serostatus in wild European rabbits. We carried out seven live-trapping sessions in three wild rabbit populations over a two-year period. Blood samples were collected to measure anti-MV and anti-RHDV antibody concentrations and to measure biochemical parameters related to organ function, protein metabolism, and nutritional status. Overall, we found no significant relationships between rabbit serostatus and biochemistry. Our main result was that rabbits that were seropositive for both MV and RHDV had low gamma glutamyltransferase concentrations. Given the robustness of our analyses, the lack of significant relationships may indicate that the biochemical parameters measured are poor proxies for serostatus. Another explanation is that wild rabbits might be producing attenuated physiological responses to these viruses because the latter are now enzootic in the study area.

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

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

Comparative quantitative monitoring of rabbit haemorrhagic disease viruses in rabbit kittens

Background

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.

Methods

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.

Results

Kittens were susceptible to infection with virus doses as low as 10 ID₅₀. 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.

Conclusions

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.

Multi-event capture-recapture modeling of host-pathogen dynamics among European rabbit populations exposed to myxoma and Rabbit Hemorrhagic Disease Viruses: common and heterogeneous patterns

Host-pathogen epidemiological processes are often unclear due both to their complexity and over-simplistic approaches used to quantify them. We applied a multi-event capture-recapture procedure on two years of data from three rabbit populations to test hypotheses about the effects on survival of, and the dynamics of host immunity to, both myxoma virus and Rabbit Hemorrhagic Disease Virus (MV and RHDV). Although the populations shared the same climatic and management conditions, MV and RHDV dynamics varied greatly among them; MV and RHDV seroprevalences were positively related to density in one population, but RHDV seroprevalence was negatively related to density in another. In addition, (i) juvenile survival was most often negatively related to seropositivity, (ii) RHDV seropositives never had considerably higher survival, and (iii) seroconversion to seropositivity was more likely than the reverse. We suggest seropositivity affects survival depending on trade-offs among antibody protection, immunosuppression and virus lethality. Negative effects of seropositivity might be greater on juveniles due to their immature immune system. Also, while RHDV directly affects survival through the hemorrhagic syndrome, MV lack of direct lethal effects means that interactions influencing survival are likely to be more complex. Multi-event modeling allowed us to quantify patterns of host-pathogen dynamics otherwise difficult to discern. Such an approach offers a promising tool to shed light on causative mechanisms.

Recovery of South Australian rabbit populations from the impact of rabbit haemorrhagic disease

Context Recovery of Australian rabbit populations from the impact of rabbit haemorrhagic disease virus (RHDV) contrasts with more prolonged suppression of wild rabbits in Europe, and has been widely discussed in the scientific community, but not yet documented in formal scientific literature. The underlying causes of recovery remain unclear, but resistance to RHDV infection has been reported in laboratory studies of wild-caught rabbits. Aims We document numerical changes in two South Australian wild rabbit populations that were initially suppressed by RHDV, and examine serological data to evaluate several alternative hypotheses for the cause of recovery. Methods Rabbit numbers were assessed from spotlight transect counts and dung mass transects between 1991 and 2011, and age and RHDV antibody sero-prevalence were estimated from rabbits shot in late summer. Key results Rabbit numbers were heavily suppressed by RHDV between 1995 and 2002, then increased 5- to 10-fold between 2003 and 2010. During the period of increase, annual RHDV infection rates remained stable or increased slightly, average age of rabbits remained stable and annual rainfall was below average. Conclusions Rabbit populations recovered but neither avoidance of RHDV infection, gradual accumulation of long-lived RHD-immune rabbits, nor high pasture productivity were contributing factors. This leaves increased annual survival from RHDV infection as the most likely cause of recovery. Implications Previously documented evidence of resistance to RHDV infection may be of little consequence to post-RHD recovery in rabbit numbers, unless the factors that influence the probability of infection also shape the course of infection and affect survival of infected rabbits.

Recolonisation of rabbit warrens following coordinated ripping programs in Victoria, south-eastern Australia

Context Warren ripping has been demonstrated to be an effective tool for controlling rabbit populations. However, few studies have examined factors influencing the rate at which ripped warrens are likely to be recolonised (i.e. be re-opened). Aims To examine factors influencing the recolonisation of ripped warrens by rabbits by using data collected on 555 warrens for up to 15 years following coordinated ripping programs at 12 sites in Victoria, south-eastern Australia. Methods Warren-monitoring data (number of active and inactive warren entrances) were analysed using discrete-time survival analysis to determine the effects of warren-level and site-level covariates on the recolonisation of ripped warrens. Key results Warren recolonisation was related to the distance between the ripped warren and the nearest active warren, the number of active entrances in the nearest warren, the initial number of active entrances in the ripped warren and the rabbit spotlight abundance index at the site. The probability of warren recolonisation was highest for ripped warrens within 1 km of an active warren and negligible beyond 3 km. The probability of warren recolonisation also increased by 22% for every increase in the rabbit spotlight count at the site by 10 rabbits km–1. Conclusions The recolonisation of ripped warrens was highly influenced by both the distance to, and size of, neighbouring active warrens. Larger warrens also appear to be preferentially recolonised compared with smaller warrens, suggesting that recolonisation of ripped areas may be related to habitat quality. The present results are consistent with ideas from classical metapopulation theory predicting that the rates of colonisation of vacant patches are dependent on both the proximity and size of the source population as well as the quality of habitat patches. Implications Although coordinated warren ripping programs are effective at achieving long-term control of rabbits, their efficiency at maintaining low rabbit populations can be increased by adopting an adaptive monitoring program that incorporates warren size and the spatial relationships among warrens, and using this information to better target maintenance-control activities.

A rapid survey method for estimating population density of European rabbits living in native vegetation

European rabbits are severe environmental pests in Australia but reporting of density-damage relationships has been hindered by a lack of simple methods to estimate the density of rabbit populations in native vegetation. A methodology for quantifying rabbit densities suitable for use in sparse populations of rabbits in conjunction with surveys of the condition of native vegetation is proposed. Dung pellets were counted in 11 629 regularly-spaced quadrats of 0.1 m2 in semiarid, coastal and cool-temperate areas of southern Australia. Mean pellet counts in latrines and the relationship between dung counts outside of latrines and the proportion of quadrats falling on latrines were quantified. This allowed density of dung pellets to be estimated by using a correction factor for latrines instead of attempting to count all pellets within quadrats that included parts of latrines. Rabbit density was calculated from pellet density based on mean pellet size, pellet breakdown rate and estimates of rabbits’ dry matter intake and digestive efficiency. Results were validated against estimates of rabbit density from long-term studies using a combination of spotlight transect counts and burrow entrance counts. The proposed methodology allows estimates of rabbit density in native vegetation to be obtained from just a few hours work and can be used in conjunction with surveys of the condition of native vegetation to quantify rabbit impacts. This methodology is seen as particularly useful in providing a tool to allow rabbit densities to be estimated and then compared with the thresholds, determined separately, at which damage occurs for given ecosystems.

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

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

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

The endemic non-pathogenic Australian rabbit calicivirus RCV-A1 is known to provide some cross protection to lethal infection with the closely related Rabbit Haemorrhagic Disease Virus (RHDV). Despite its obvious negative impacts on viral biocontrol of introduced European rabbits in Australia, little is known about the extent and mechanisms of this cross protection. In this study 46 rabbits from a colony naturally infected with RCV-A1 were exposed to RHDV. Survival rates and survival times did not correlate with titres of serum antibodies specific to RCV-A1 or cross reacting to RHDV, but were instead influenced by the time between infection with the two viruses, demonstrating for the first time that the cross protection to lethal RHDV infection is transient. These findings are an important step towards a better understanding of the complex interactions of co-occurring pathogenic and non-pathogenic lagoviruses.

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

BACKGROUND

Antibodies to non-pathogenic rabbit caliciviruses (RCVs) cross-react in serological tests for rabbit hemorrhagic disease virus (RHDV) and vice versa, making epidemiological studies very difficult where both viruses occur. It is important to understand the distribution and interaction of the two viruses because the highly pathogenic RHDV has been used as a biocontrol agent for wild rabbits in Australia and New Zealand for the past 17 years. The presence of the benign RCV Australia 1 (RCV-A1) is considered a key factor for the failure of RHDV mediated rabbit control in cooler, wetter areas of Australia.

RESULTS

A highly sensitive and specific blocking ELISA was developed for the detection of RCV-A1 antibodies. When sera from rabbits with a known infection history for either RCV-A1 or RHDV were tested, this assay showed 100% sensitivity and no cross-reactivity with RHDV sera (100% specificity).

CONCLUSIONS

This new ELISA not only allows the detection of RCV-A1 at a population level, but also permits the serological status of individual rabbits to be determined more reliably than previously described methods. This robust and simple to perform assay is therefore the tool of choice for studying RCV-A1 epidemiology in Australian wild rabbit populations.

Rabbit haemorrhagic disease: advantages of cELISA in assessing immunity in wild rabbits (Oryctolagus cuniculus)

Rabbit haemorrhagic disease (RHD) is an acute fatal disease of domestic and wild European rabbits (Oryctolagus cuniculus) caused by RHD virus (RHDV). Accurate assessment of immunity is of great importance for the conservation and control of wild rabbits. We evaluated a competitive ELISA (cELISA) against isotype ELISAs for assessing the protective immunity against the disease by challenging 50 wild-caught rabbits with a lethal dose of RHDV. Death or survival to the challenge was used as a criterion to determine the performance characteristics of the assay for the assessment of immunity in rabbits. At 1:10 dilution, a serum exhibiting ≥ 25% inhibition (1:10(25)) was regarded as the presence of RHDV-specific antibodies. Eleven of 16 (68.8%) rabbits with antibodies at 1:10(25) (<1:40) died of RHD. When the cut-off was moved from 25% to 50% inhibition (1:10(50)) at 1:10 serum dilution, the assay sensitivity, specificity and accuracy for the protective immunity were improved from 84%, 54.2% and 69.4% to 84%, 100% and 91.8%, respectively. We also demonstrated at the epitope amino acid sequence level why the presence of the RHDV-cross reactive benign rabbit calicivirus, which interfered with isotype ELISAs, had little impact on the specificity of the cELISA for the diagnosis of RHDV infection. The presence of RHDV-specific antibody at 1:10(50) by the cELISA is a reliable indicator for the protective immunity. In contrast to isotype ELISAs, the cELISA is a valuable specific tool for monitoring the herd immunity to RHD for the conservation and management of wild rabbits in the field.