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

Characterisation of Lagovirus europaeus GI-RHDVs (Rabbit Haemorrhagic Disease Viruses) in Terms of Their Pathogenicity and Immunogenicity

Rabbit haemorrhagic disease viruses (RHDV) belong to the family Caliciviridae, genus Lagovirus europaeus, genogroup GI, comprising four genotypes GI.1-GI.4, of which the genotypes GI.1 and GI.2 are pathogenic RHD viruses, while the genotypes GI.3 and GI.4 are non-pathogenic RCV (Rabbit calicivirus) viruses. Among the pathogenic genotypes GI.1 and GI.2 of RHD viruses, an antigenic variant of RHDV, named RHDVa-now GI.1a-RHDVa, was distinguished in 1996; and in 2010, a variant of RHDV-named RHDVb, later RHDV2 and now GI.2-RHDV2/b-was described; and recombinants of these viruses were registered. Pathogenic viruses of the genotype GI.1 were the cause of a disease described in 1984 in China in domestic (Oryctolagus (O.) cuniculus domesticus) and wild (O. cuniculus) rabbits, characterised by a very rapid course and a mortality rate of 90-100%, which spread in countries all over the world and which has been defined since 1989 as rabbit haemorrhagic disease. It is now accepted that GI.1-RHDV, including GI.1a-RHDVa, cause the predetermined primary haemorrhagic disease in domestic and wild rabbits, while GI.2-RHDV2/b cause it not only in rabbits, including domestic rabbits' young up to 4 weeks and rabbits immunised with rabbit haemorrhagic disease vaccine, but also in five various species of wild rabbits and seven different species of hares, as well as wild ruminants: mountain muskoxen and European badger. Among these viruses, haemagglutination-positive, doubtful and harmful viruses have been recorded and described and have been shown to form phylogenogroups, immunotypes, haematotypes and pathotypes, which, together with traits that alter and expand their infectious spectrum (rabbit, hare, wild ruminant, badger and various rabbit and hare species), are the determinants of their pathogenicity (infectivity) and immunogenicity and thus shape their virulence. These relationships are the aim of our consideration in this article.

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.

Comparison of the Impact between Classical and Novel Strains of Rabbit Haemorrhagic Disease on Wild Rabbit Populations in Spain

The outbreaks of two strains of rabbit haemorrhagic disease (RHD) (GI.1 and GI.2) in the Iberian Peninsula have caused substantial economic losses in commercial rabbitries and have affected the conservation of rabbit-sensitive predators due to the dramatic decline of their natural populations. However, the assessment of the impact of both RHD strains on wild rabbit populations has been limited to a few small-scale studies. Little is known about the overall impact within its native range. In this study, we described and compared the effects of GI.1 and GI.2 countrywide by using time series of hunting bag data widely available across the country and compared their trend during the first eight years after the first outbreak of GI.1 (i.e., 1998) and GI.2 (i.e., 2011), respectively. We used Gaussian generalised additive models (GAM) with the number of hunted rabbits as the response variable and year as the predictor to evaluate the non-linear temporal dynamics of the population at the national and regional community levels. The first GI.1 caused a population decline of around 53%, affecting most Spanish regional communities where the disease occurred. The positive trend observed after GI.1 in Spain ended with the initial outbreak of GI.2, which did not appear to cause a national population decline. In contrast, we found significant variability in the rabbit population trend among regional communities, where some increased, and others decreased. Such a disparity is unlikely to be explained by a single factor; rather, it appears to result from several factors, such as climatic conditions, host resistance improvement, virulence attenuation, or population density. Our study suggests that a national comprehensive hunting bag series could aid in elucidating the differences in the impact of emerging diseases on a large scale. Future research should focus on national longitudinal serological studies to shed light on the immunological status of rabbit populations in different regions to better understand the evolution of RHD strains and the resistance gained by the wild populations.

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.

Comparative Epidemiology of Rabbit Haemorrhagic Disease Virus Strains from Viral Sequence Data

Since their introduction in 1859, European rabbits (Oryctolagus cuniculus) have had a devastating impact on agricultural production and biodiversity in Australia, with competition and land degradation by rabbits being one of the key threats to agricultural and biodiversity values in Australia. Biocontrol agents, with the most important being the rabbit haemorrhagic disease virus 1 (RHDV1), constitute the most important landscape-scale control strategies for rabbits in Australia. Monitoring field strain dynamics is complex and labour-intensive. Here, using phylodynamic models to analyse the available RHDV molecular data, we aimed to: investigate the epidemiology of various strains, use molecular data to date the emergence of new variants and evaluate whether different strains are outcompeting one another. We determined that the two main pathogenic lagoviruses variants in Australia (RHDV1 and RHDV2) have had similar dynamics since their release, although over different timeframes (substantially shorter for RHDV2). We also found a strong geographic difference in their activities and evidence of overall competition between the two viruses.

Immunological Cross-Protection between Different Rabbit Hemorrhagic Disease Viruses—Implications for Rabbit Biocontrol and Vaccine Development

The use of rabbit hemorrhagic disease virus (RHDV) as a biocontrol agent to control feral rabbit populations in Australia, in combination with circulating endemic strains, provides a unique environment to observe the interactions between different lagoviruses competing for the same host. Following the arrival of RHDV2 (GI.2) in Australia, it became necessary to investigate the potential for immunological cross-protection between different variants, and the implications of this for biocontrol programs and vaccine development. Laboratory rabbits of various immune status—(1) rabbits with no detectable immunity against RHDV; (2) rabbits with experimentally acquired immunity after laboratory challenge; (3) rabbits immunised with a GI.2-specific or a multivalent RHDV inactivated virus prototype vaccine; or (4) rabbits with naturally acquired immunity—were challenged with one of three different RHDV variants (GI.1c, GI.1a or GI.2). The degree of cross-protection observed in immune rabbits was associated with the variant used for challenge, infectious dose of the virus and age, or time since acquisition of the immunity, at challenge. The immune status of feral rabbit populations should be determined prior to intentional RHDV release because of the high survival proportions in rabbits with pre-existing immunity. In addition, to protect domestic rabbits in Australia, a multivalent RHDV vaccine should be considered because of the limited cross-protection observed in rabbits given monovalent vaccines.

Passive Immunisation against RHDV2 Induces Protection against Disease but Not Infection

Rabbit haemorrhagic disease virus 2 (RHDV2) is a lagovirus in the family Caliciviridae. The closely related Rabbit haemorrhagic disease virus (RHDV, termed RHDV1 throughout this manuscript for clarity) has been used extensively as a biocontrol agent in Australia since the mid-1990s to manage wild rabbit populations, a major economic and environmental pest species. Releasing RHDV1 into populations with a high proportion of rabbits less than 8–10 weeks of age leads to non-lethal infection in many of these young animals, with subsequent seroconversion and long-term immunity against reinfection. In contrast, RHDV2 causes lethal disease even in young rabbits, potentially offering substantial benefits for rabbit management programs over RHDV1. However, it is not clear how acquired resistance from maternal antibodies may influence immunity after RHDV2 infection. In this study, we assessed serological responses after RHDV2 challenge in young rabbits of three different ages (5-, 7-, or 9-weeks-old) that were passively immunised with either high- (titre of 2560 by RHDV IgG ELISA; 2.41 mg/mL total protein) or low- (titre of 160–640 by RHDV IgG ELISA; 1.41 mg/mL total protein) dose RHDV2 IgG to simulate maternal antibodies. All rabbits treated with a high dose and 75% of those treated with a low dose of RHDV2 IgG survived virus challenge. Surviving animals developed robust lagovirus-specific IgA, IgM, and IgG responses within 10 days post infection. These findings demonstrate that the protection against RHDV2 conferred by passive immunisation is not sterilising. Correspondingly, this suggests that the presence of maternal antibodies in wild rabbit populations may impede the effectiveness of RHDV2 as a biocontrol.

A Review on the Methods Used for the Detection and Diagnosis of Rabbit Hemorrhagic Disease Virus (RHDV)

Since the early 1980s, the European rabbit (Oryctolagus cuniculus) has been threatened by the rabbit hemorrhagic disease (RHD). The disease is caused by a lagovirus of the family Caliciviridae, the rabbit hemorrhagic disease virus (RHDV). The need for detection, identification and further characterization of RHDV led to the development of several diagnostic tests. Owing to the lack of an appropriate cell culture system for in vitro propagation of the virus, much of the methods involved in these tests contributed to our current knowledge on RHD and RHDV and to the development of vaccines to contain the disease. Here, we provide a comprehensive review of the RHDV diagnostic tests used since the first RHD outbreak and that include molecular, histological and serological techniques, ranging from simpler tests initially used, such as the hemagglutination test, to the more recent and sophisticated high-throughput sequencing, along with an overview of their potential and their limitations.

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.

Immunity against Lagovirus europaeus and the Impact of the Immunological Studies on Vaccination

In the early 1980s, a highly contagious viral hemorrhagic fever in rabbits (Oryctolagus cuniculus) emerged, causing a very high rate of mortality in these animals. Since the initial occurrence of the rabbit hemorrhagic disease virus (RHDV), several hundred million rabbits have died after infection. The emergence of genetically-different virus variants (RHDV GI.1 and GI.2) indicated the very high variability of RHDV. Moreover, with these variants, the host range broadened to hare species (Lepus). The circulation of RHDV genotypes displays different virulences and a limited induction of cross-protective immunity. Interestingly, juvenile rabbits (<9 weeks of age) with an immature immune system display a general resistance to RHDV GI.1, and a limited resistance to RHDV GI.2 strains, whereas less than 3% of adult rabbits survive an infection by either RHDV GI.1. or GI.2. Several not-yet fully understood phenomena characterize the RHD. A very low infection dose followed by an extremely rapid viral replication could be simplified to the induction of a disseminated intravascular coagulopathy (DIC), a severe loss of lymphocytes—especially T-cells—and death within 36 to 72 h post infection. On the other hand, in animals surviving the infection or after vaccination, very high titers of RHDV-neutralizing antibodies were induced. Several studies have been conducted in order to deepen the knowledge about the virus’ genetics, epidemiology, RHDV-induced pathology, and the anti-RHDV immune responses of rabbits in order to understand the phenomenon of the juvenile resistance to this virus. Moreover, several approaches have been used to produce efficient vaccines in order to prevent an infection with RHDV. In this review, we discuss the current knowledge about anti-RHDV resistance and immunity, RHDV vaccination, and the further need to establish rationally-based RHDV vaccines.

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.

Culture and differentiation of rabbit intestinal organoids and organoid-derived cell monolayers

Organoids emulate many aspects of their parental tissue and are therefore used to study pathogen-host interactions and other complex biological processes. Here, we report a robust protocol for the isolation, maintenance and differentiation of rabbit small intestinal organoids and organoid-derived cell monolayers. Our rabbit intestinal spheroid and monolayer cultures grew most efficiently in L-WRN-conditioned medium that contained Wnt, R-spondin and Noggin, and that had been supplemented with ROCK and TGF-β inhibitors. Organoid and monolayer differentiation was initiated by reducing the concentration of the L-WRN-conditioned medium and by adding ROCK and Notch signalling inhibitors. Immunofluorescence staining and RT-qPCR demonstrated that our organoids contained enterocytes, enteroendocrine cells, goblet cells and Paneth cells. Finally, we infected rabbit organoids with Rabbit calicivirus Australia-1, an enterotropic lagovirus that—like many other caliciviruses—does not grow in conventional cell culture. Despite testing various conditions for inoculation, we did not detect any evidence of virus replication, suggesting either that our organoids do not contain suitable host cell types or that additional co-factors are required for a productive infection of rabbit organoids with Rabbit calicivirus Australia-1.

Genomic insights into a population of introduced European rabbits Oryctolagus cuniculus in Australia and the development of genetic resistance to rabbit hemorrhagic disease virus

The European rabbit (Oryctolagus cuniculus) is one of the most devastating invasive species in Australia. Since the 1950s, myxoma virus (MYXV) and rabbit haemorrhagic disease virus (RHDV) have been used to manage overabundant rabbit populations. Resistance to MYXV was observed within a few years of the release. More recently, resistance to lethal RHDV infection has also been reported, undermining the efficiency of landscape-scale rabbit control. Previous studies suggest that genetic resistance to lethal RHDV infection may differ locally between populations, yet the mechanisms of genetic resistance remain poorly understood. Here, we used genotyping by sequencing (GBS) data representing a reduced representation of the genome, to investigate Australian rabbit populations. Our aims were to understand the relationship between populations and identify possible genomic signatures of selection for RHDV resistance. One population we investigated had previously been reported to show levels of resistance to lethal RHDV infection. This population was compared to three other populations with lower or no previously reported RHDV resistance. We identified a set of novel candidate genes that could be involved in host-pathogen interactions such as virus binding and infection processes. These genes did not overlap with previous studies on RHDV resistance carried out in different rabbit populations, suggesting that multiple mechanisms are feasible. These findings provide useful insights into the different potential mechanisms of genetic resistance to RHDV virus which will inform future functional studies in this area.

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.

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.

Identification of a new European rabbit IgA with a serine-rich hinge region

In mammals, the most striking IgA system belongs to Lagomorpha. Indeed, 14 IgA subclasses have been identified in European rabbits, 11 of which are expressed. In contrast, most other mammals have only one IgA, or in the case of hominoids, two IgA subclasses. Characteristic features of the mammalian IgA subclasses are the length and amino acid sequence of their hinge regions, which are often rich in Pro, Ser and Thr residues and may also carry Cys residues. Here, we describe a new IgA that was expressed in New Zealand White domestic rabbits of IGHVa1 allotype. This IgA has an extended hinge region containing an intriguing stretch of nine consecutive Ser residues and no Pro or Thr residues, a motif exclusive to this new rabbit IgA. Considering the amino acid properties, this hinge motif may present some advantage over the common IgA hinge by affording novel functional capabilities. We also sequenced for the first time the IgA14 CH2 and CH3 domains and showed that IgA14 and IgA3 are expressed.

Elucidation of the pathology and tissue distribution of Lagovirus europaeus GI.2/RHDV2 (rabbit haemorrhagic disease virus 2) in young and adult rabbits (Oryctolagus cuniculus)

Lagovirus europaeus GI.2, also known as RHDV2 or RHDVb, is an emerging virus that causes rabbit haemorrhagic disease (RHD) in European rabbits (Oryctolagus cuniculus). In contrast to L. europaeus GI.1 (or RHDV/RHDVa) viruses that are only pathogenic for adults, GI.2 causes clinical disease in both adults and kittens. However, detailed descriptions of the pathology of this virus that may provide insight into its pathogenicity and emergence are lacking. Using an Australian GI.2 field strain isolated in 2015, we provide the first detailed description of pathology, viral antigen distribution and tissue load of GI.2 in adult and 5-week old New Zealand white rabbits using histology, immunohistochemistry and RT-qPCR. Liver was the target organ, but in contrast to GI.1 viruses, lesions and inflammatory responses did not differ between adults and kittens. Lymphocytic inflammation, proposed to be protective in kittens infected with GI.1, was notably absent. We also present the first descriptions of bone marrow changes in RHD, including decreased myeloid-to-erythroid ratio. Consistent with other pathogenic lagoviruses, intracellular viral antigen was demonstrated in hepatocytes and cells of the mononuclear phagocytic system. In terminal stages of disease, viral loads were highest in liver, serum and spleen. Despite the small sample size, our data suggest that unlike early European GI.2 strains, the pathogenicity of the Australian GI.2 virus is similar to GI.1 viruses. Additionally, GI.2 was fatal for all (n = 5) inoculated kittens in this study. This may significantly alter RHD epidemiology in the field, and may impact biocontrol programs for invasive rabbits in Australia where GI.1 viruses are intentionally released.

Rabbit haemorrhagic disease: Cross-protection and comparative pathogenicity of GI.2/RHDV2/b and GI.1b/RHDV lagoviruses in a challenge trial

European rabbits (Oryctolagus cuniculus) are severely affected by rabbit haemorrhagic disease (RHD). Caused by a lagovirus, the disease leads to losses in the rabbit industry and has implications for wildlife conservation. Past RHD outbreaks have been caused by GI.1/RHDV genotype viruses. A new virus belonging to the GI.2/RHDV2/b genotype emerged in 2010, quickly spreading and replacing the former in several countries; however, limited data are available on its pathogenicity and epidemiological factors. The present work extends these issues and evaluates cross-protection between both genotypes. Ninety-four and 88 domestic rabbits were challenged with GI.2/RHDV2/b and GI.1b/RHDV variant isolates, respectively. Cross-protection was determined by a second challenge on survivors with the corresponding strain. Mortality by GI.2/RHDV2/b was highly variable due to unknown individual factors, whereas mortality by GI.1b/RHDV was associated with age. Mortality in rabbits < 4 weeks old was 84%, higher than previously reported. Cross-protection was not identical between the two viruses because the ratio of mortality rate ratios for the first and second challenges was 3.80 ± 2.68 times higher for GI.2/RHDV2/b than it was for GI.1b/RHDV. Rabbit susceptibility to GI.2/RHDV2/b varied greatly and appeared to be modulated by the innate functionality of the immune response and/or its prompt activation by other pathogens. GI.1b/RHDV pathogenicity appeared to be associated with undetermined age-related factors. These results suggest that GI.2/RHDV2/b may interact with other pathogens at the population level but does not satisfactorily explain the GI.1b/RHDV virus's quick replacement.

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.

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.

Timing and severity of immunizing diseases in rabbits is controlled by seasonal matching of host and pathogen dynamics

Infectious diseases can exert a strong influence on the dynamics of host populations, but it remains unclear why such disease-mediated control only occurs under particular environmental conditions. We used 16 years of detailed field data on invasive European rabbits (Oryctolagus cuniculus) in Australia, linked to individual-based stochastic models and Bayesian approximations, to test whether (i) mortality associated with rabbit haemorrhagic disease (RHD) is driven primarily by seasonal matches/mismatches between demographic rates and epidemiological dynamics and (ii) delayed infection (arising from insusceptibility and maternal antibodies in juveniles) are important factors in determining disease severity and local population persistence of rabbits. We found that both the timing of reproduction and exposure to viruses drove recurrent seasonal epidemics of RHD. Protection conferred by insusceptibility and maternal antibodies controlled seasonal disease outbreaks by delaying infection; this could have also allowed escape from disease. The persistence of local populations was a stochastic outcome of recovery rates from both RHD and myxomatosis. If susceptibility to RHD is delayed, myxomatosis will have a pronounced effect on population extirpation when the two viruses coexist. This has important implications for wildlife management, because it is likely that such seasonal interplay and disease dynamics has a strong effect on long-term population viability for many species.

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.

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

BACKGROUND

Australia relies heavily on rabbit haemorrhagic disease virus (RHDV) for the biological control of introduced European wild rabbits Oryctolagus cuniculus, which are significant economic and environmental pests. An endemic non-pathogenic rabbit calicivirus termed RCV-A1 also occurs in wild rabbits in Australian and provides partial protection against lethal RHDV infection, thus interfering with effective rabbit control. Despite its obvious importance for rabbit population management, little is known about the epidemiology of this benign rabbit calicivirus.

METHODS

We determined the continent-wide distribution and prevalence of RCV-A1 by analysing 1,805 serum samples from wild rabbit populations at 78 sites across Australia for the presence of antibodies to RCV-A1 using a serological test that specifically detects RCV-A1 antibodies and does not cross-react with co-occurring RHDV antibodies. We also investigated possible correlation between climate variables and prevalence of RCV-A1 by using generalised linear mixed effect models.

RESULTS

Antibodies to RCV-A1 were predominantly detected in rabbit populations in cool, high rainfall areas of the south-east and south-west of the continent. There was strong support for modelling RCV-A1 prevalence as a function of average annual rainfall and minimum temperature. The best ranked model explained 26% of the model structural deviance. According to this model, distribution and prevalence of RCV-A1 is positively correlated with periods of above average rainfall and negatively correlated with periods of drought.

IMPLICATIONS

Our statistical model of RCV-A1 prevalence will greatly increase our understanding of RCV-A1 epidemiology and its interaction with RHDV in Australia. By defining the environmental conditions associated with the prevalence of RCV-A1, it also contributes towards understanding the distribution of similar viruses in New Zealand and Europe.

Viral biocontrol: grand experiments in disease emergence and evolution

Although viral emergence is commonly associated with cross-species transmission, the processes and determinants of viral evolution in a novel host environment are poorly understood. We address key questions in virus emergence and evolution using data generated from two unique natural experiments: the deliberate release of myxoma virus (MYXV) and rabbit hemorrhagic disease virus (RHDV) as biological control (biocontrol) agents against the European rabbit in Australia, and which have been of enormous benefit to Australia's ecosystem and agricultural industries. Notably, although virulence evolution in MYXV and RHDV followed different trajectories, a strongly parallel evolutionary process was observed in Australia and Europe. These biocontrol agents were also characterized by a lack of transmission to nontarget host species, suggesting that there are major barriers to successful emergence.

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.

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.

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.

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.