Bayesian evaluation of temporal changes in sensitivity and specificity of three serological tests for multiple circulating strains of rabbit haemorrhagic disease virus.

Competition and indirect ELISAs are currently being used to monitor rabbit haemorrhagic disease viruses (RHDV1 and RHDV2) in rabbits worldwide. Temporal changes in the sensitivity (Se) and specificity (Sp) of RHDV1 competition-ELISA (cELISA1), RHDV2 competition-ELISA (cELISA2), and RHDV1 Immunoglobulin G (IgG1) ELISA, were investigated using Bayesian Latent Class models (BCLM) in the Australian wild rabbit population where both viruses circulate simultaneously and a long-term serological dataset exists. When cELISA1 was compared to IgG1 ELISA, the Se of cELISA1 improved while the Sp of IgG1 ELISA declined over the 2011-21. This corresponded with a decline in the true RHDV1 prevalence in 2018-21, suggesting that a large proportion of RHDV1 exposed rabbits survived the introduction and dominance of RHDV2 up to approximately 2017/2018, after which they died and were not replaced. The Se and Sp estimates for 2014-15 for both cELISA1 and IgG1 ELISA, and the true prevalence when analysing all three tests together were similar to those obtained from the analysis of cELISA1/IgG1 ELISA. The same was also true for the Se and Sp of cELISA2 and IgG1 ELISA estimates from 2018 onwards. This suggests that RHDV1 was the dominant infection type in 2014-15, but RHDV2 was the dominant infection type in 2018-21. Further, the increase in Se of cELISA2 and the low Sp of IgG1 ELISA in the cELISA2/IgG1 ELISA analysis, compared to the Se of cELISA2 and Sp of IgG1 ELISA when analysing all three tests together suggests that the underlying infection status was more influenced by RHDV2 and that the higher Se of IgG1 ELISA is due to cross-reaction of RHDV2 antibodies on IgG1 ELISA. The true prevalence data suggest that RHDV2 exposure peaked in 2017. Our findings show that test characteristics changed in response to the changing virus prevalences over time. IgG1 ELISA, currently having a high Se, should be used to monitor both viruses and will perform better than both cELISAs.

Using genetic analysis to determine the distribution, prevalence and diversity of Eimeria species in pest rabbits (Oryctolagus cuniculus) in Australia.

To genetically assess the Australian distribution and frequency of Eimeria species in wild rabbits, with a primary focus on Eimeria intestinalis and Eimeria flavescens as possible additional agents of rabbit biocontrol, the distal colon and faecal samples from wild rabbits sourced from 26 Australian locations with mean annual rainfalls of between 252 mm and 925 mm were analysed using amplicon sequencing of the ITS1 region. Contrary to previous microscopy studies which had only detected E. flavescens on mainland Australia at Wellstead in south-west Western Australia, we detected this species at all 23 effectively sampled sites. The more pathogenic E. intestinalis was only found at 52.2% of sites. Three unique Eimeria genotypes were detected that did not align to the 11 published sequences using a pairwise-match threshold of 90%, and may represent unsequenced known species or novel species. One genotype we termed E. Au19SH and was detected at 20 sites, E. Au19CO was detected at eight sites, and E. Au19CN was detected in one rabbit at Crows Nest (Qld). Site diversity ranged from only five Eimeria species at Boboyan (ACT) to 13 unique sequences at Cargo (NSW). Eimeria diversity in individual rabbits ranged from 11 unique sequences in a rabbit at Wellstead (WA) and a rabbit at Cargo (NSW), to one in 17 rabbits and zero in six rabbits. The three rabbit age classes averaged 4.3 Eimeria species per rabbit. No relationship was found between the number of Eimeria species detected and mean annual rainfall. As Eimeria species were found to be fairly ubiquitous at most sites they appear to be an unlikely additional candidate to assist the control of pest rabbits in Australia.

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 pre-screening 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 out of 22), followed by K5 (53.8%, seven out of 13) and CZECH (40.0%, two out of five), 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 out of 52) followed by CZECH (46.0%, 23 out of 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 with 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.

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.

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 Australian National Rabbit Database: 50 yr of population monitoring of an invasive species.

With ongoing introductions into Australia since the 1700s, the European rabbit (Oryctolagus cuniculus) has become one of the most widely distributed and abundant vertebrate pests, adversely impacting Australia's biodiversity and agroeconomy. To understand the population and range dynamics of the species and its impacts better, occurrence and abundance data have been collected by researchers and citizens from sites covering a broad spectrum of climatic and environmental conditions in Australia. The lack of a common and accessible repository for these data has, however, limited their use in determining important spatiotemporal drivers of the structure and dynamics of the geographical range of rabbits in Australia. To meet this need, we created the Australian National Rabbit Database, which combines more than 50 yr of historical and contemporary survey data collected from throughout the range of the species in Australia. The survey data, obtained from a suite of complementary monitoring methods, were combined with high-resolution weather, climate, and environmental information, and an assessment of data quality. The database provides records of rabbit occurrence (689,265 records) and abundance (51,241 records, >120 distinct sites) suitable for identifying the spatiotemporal drivers of the rabbit's distribution and for determining spatial patterns of variation in its key life-history traits, including maximum rates of population growth. Because all data are georeferenced and date stamped, they can be coupled with information from other databases and spatial layers to explore the potential effects of rabbit occurrence and abundance on Australia's native wildlife and agricultural production. The Australian National Rabbit Database is an important tool for understanding and managing the European rabbit in its invasive range and its effects on native biodiversity and agricultural production. It also provides a valuable resource for addressing questions related to the biology, success, and impacts of invasive species more generally. No copyright or proprietary restrictions are associated with the use of this data set other than citation of this Data Paper.

Rabbit Hemorrhagic Disease Virus 2 (RHDV2; GI.2) Is Replacing Endemic Strains of RHDV in the Australian Landscape within 18 Months of Its Arrival.

Rabbit hemorrhagic disease virus 2 (RHDV2; Lagovirus GI.2) is a pathogenic calicivirus that affects European rabbits (Oryctolagus cuniculus) and various hare (Lepus) species. GI.2 was first detected in France in 2010 and subsequently caused epidemics in wild and domestic lagomorph populations throughout Europe. In May 2015, GI.2 was detected in Australia. Within 18 months of its initial detection, GI.2 had spread to all Australian states and territories and rapidly became the dominant circulating strain, replacing Rabbit hemorrhagic disease virus (RHDV/GI.1) in mainland Australia. Reconstruction of the evolutionary history of 127 Australian GI.2 isolates revealed that the virus arrived in Australia at least several months before its initial description and likely circulated unnoticed in wild rabbit populations in the east of the continent prior to its detection. GI.2 sequences isolated from five hares clustered with sequences from sympatric rabbit populations sampled contemporaneously, indicating multiple spillover events into hares rather than an adaptation of the Australian GI.2 to a new host. Since the presence of GI.2 in Australia may have wide-ranging consequences for rabbit biocontrol, particularly with the release of the novel biocontrol agent GI.1a/RHDVa-K5 in March 2017, ongoing surveillance is critical to understanding the interactions of the various lagoviruses in Australia and their impact on host populations.IMPORTANCE This study describes the spread and distribution of Rabbit hemorrhagic disease virus 2 (GI.2) in Australia since its first detection in May 2015. Within the first 18 months following its detection, RHDV2 spread from east to west across the continent and became the dominant strain in all mainland states of Australia. This has important implications for pest animal management and for owners of pet and farmed rabbits, as there currently is no effective vaccine available in Australia for GI.2. The closely related RHDV (GI.1) is used to control overabundant wild rabbits, a serious environmental and agricultural pest in this country, and it is currently unclear how the widespread circulation of GI.2 will impact ongoing targeted wild rabbit management operations.

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

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

Proposal for a unified classification system and nomenclature of lagoviruses.

Lagoviruses belong to the Caliciviridae family. They were first recognized as highly pathogenic viruses of the European rabbit (Oryctolagus cuniculus) and European brown hare (Lepus europaeus) that emerged in the 1970-1980s, namely, rabbit haemorrhagic disease virus (RHDV) and European brown hare syndrome virus (EBHSV), according to the host species from which they had been first detected. However, the diversity of lagoviruses has recently expanded to include new related viruses with varying pathogenicity, geographic distribution and host ranges. Together with the frequent recombination observed amongst circulating viruses, there is a clear need to establish precise guidelines for classifying and naming lagovirus strains. Therefore, here we propose a new nomenclature based on phylogenetic relationships. In this new nomenclature, a single species of lagovirus would be recognized and called Lagovirus europaeus. The species would be divided into two genogroups that correspond to RHDV- and EBHSV-related viruses, respectively. Genogroups could be subdivided into genotypes, which could themselves be subdivided into phylogenetically well-supported variants. Based on available sequences, pairwise distance cutoffs have been defined, but with the accumulation of new sequences these cutoffs may need to be revised. We propose that an international working group could coordinate the nomenclature of lagoviruses and any proposals for revision.

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

UNLABELLED: 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.

Increased virulence of rabbit haemorrhagic disease virus associated with genetic resistance in wild Australian rabbits (Oryctolagus cuniculus).

The release of myxoma virus (MYXV) and Rabbit Haemorrhagic Disease Virus (RHDV) in Australia with the aim of controlling overabundant rabbits has provided a unique opportunity to study the initial spread and establishment of emerging pathogens, as well as their co-evolution with their mammalian hosts. In contrast to MYXV, which attenuated shortly after its introduction, rapid attenuation of RHDV has not been observed. By studying the change in virulence of recent field isolates at a single field site we show, for the first time, that RHDV virulence has increased through time, likely because of selection to overcome developing genetic resistance in Australian wild rabbits. High virulence also appears to be favoured as rabbit carcasses, rather than diseased animals, are the likely source of mechanical insect transmission. These findings not only help elucidate the co-evolutionary interaction between rabbits and RHDV, but reveal some of the key factors shaping virulence evolution.

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.

Rabbit haemorrhagic disease: virus persistence and adaptation in Australia.

In Australia, the rabbit haemorrhagic disease virus (RHDV) has been used since 1996 to reduce numbers of introduced European rabbits (Oryctolagus cuniculus) which have a devastating impact on the native Australian environment. RHDV causes regular, short disease outbreaks, but little is known about how the virus persists and survives between epidemics. We examined the initial spread of RHDV to show that even upon its initial spread, the virus circulated continuously on a regional scale rather than persisting at a local population level and that Australian rabbit populations are highly interconnected by virus-carrying flying vectors. Sequencing data obtained from a single rabbit population showed that the viruses that caused an epidemic each year seldom bore close genetic resemblance to those present in previous years. Together, these data suggest that RHDV survives in the Australian environment through its ability to spread amongst rabbit subpopulations. This is consistent with modelling results that indicated that in a large interconnected rabbit meta-population, RHDV should maintain high virulence, cause short, strong disease outbreaks but show low persistence in any given subpopulation. This new epidemiological framework is important for understanding virus-host co-evolution and future disease management options of pest species to secure Australia's remaining natural biodiversity.

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.

Rabbit haemorrhagic disease: applying Occam's razor to competing hypotheses.

Rabbit haemorrhagic disease virus (RHDV) is a highly virulent lagovirus endemic in Europe and Australasian populations of the European rabbit, Oryctolagus cuniculus. It has also caused several unexplained disease outbreaks in domestic European rabbits in North America. Non-pathogenic spread of RHDV leading to persistent infection which later reactivated has recently been proposed as the cause of overt disease and death of a pet rabbit in Canada, the first confirmed case of Rabbit haemorrhagic disease in that country. We suggest that there is little evidence to support non-pathogenic spread of virulent RHDV, some evidence that is contradictory, and evidence to support a simpler alternative hypothesis. RHDV can be spread over long distances between sparse rabbit populations by fomites or flying insects. Although highly pathogenic, RHDV can be limited in its spread within rabbit populations, or its presence masked by closely related but non-pathogenic lagoviruses which can provide protection against acute disease. In the absence of any evidence from clinical studies to support reactivation of persistent RHDV infection, the simpler explanation seems more likely to be correct.