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

Rabbit hemorrhagic disease virus 2, 2010-2023: a review of global detections and affected species

Rabbit hemorrhagic disease virus 2/genotype GI.2 (RHDV2/GI.2; Caliciviridae, Lagovirus) causes a highly contagious disease with hepatic necrosis and disseminated intravascular coagulation in several Leporidae species. RHDV2 was first detected in European rabbits (Oryctolagus cuniculus) in France in 2010 and has since spread widely. We gather here data on viral detections reported in various countries and affected species, and discuss pathology, genetic differences, and novel diagnostic aspects. RHDV2 has been detected almost globally, with cases reported in Europe, Africa, Oceania, Asia, and North America as of 2023. Since 2020, large scale outbreaks have occurred in the United States and Mexico and, at the same time, cases have been reported for the first time in previously unaffected countries, such as China, Japan, Singapore, and South Africa, among others. Detections have been notified in domestic and wild European rabbits, hares and jackrabbits (Lepus spp.), several species of cottontail and brush rabbits (Sylvilagus spp.), pygmy rabbits (Brachylagus idahoensis), and red rock rabbits (Pronolagus spp.). RHDV2 has also been detected in a few non-lagomorph species. Detection of RHDV2 causing RHD in Sylvilagus spp. and Leporidae species other than those in the genera Oryctolagus and Lepus is very novel. The global spread of this fast-evolving RNA virus into previously unexploited geographic areas increases the likelihood of host range expansion as new species are exposed; animals may also be infected by nonpathogenic caliciviruses that are disseminated by almost all species, and with which genetic recombination may occur.

First Detection and Circulation of RHDV2 in New Zealand

Rabbit haemorrhage disease virus 2 (RHDV2) is a highly pathogenic lagovirus that causes lethal disease in rabbits and hares (lagomorphs). Since its first detection in Europe in 2010, RHDV2 has spread worldwide and has been detected in over 35 countries so far. Here, we provide the first detailed report of the detection and subsequent circulation of RHDV2 in New Zealand. RHDV2 was first detected in New Zealand in 2018, with positive samples retrospectively identified in December 2017. Subsequent time-resolved phylogenetic analysis suggested a single introduction into the North Island between March and November 2016. Genetic analysis identified a GI.3P-GI.2 variant supporting a non-Australian origin for the incursion; however, more accurate identification of the source of the incursion remains challenging due to the wide global distribution of the GI.3P-GI.2 variant. Furthermore, our analysis suggests the spread of the virus between the North and South Islands of New Zealand at least twice, dated to mid-2017 and around 2018. Further phylogenetic analysis also revealed a strong phylogeographic pattern. So far, no recombination events with endemic benign New Zealand rabbit caliciviruses have been identified. This study highlights the need for further research and surveillance to monitor the distribution and diversity of lagoviruses in New Zealand and to detect incursions of novel variants.

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.

Distribution and Genetic Diversity of Hepatitis E Virus in Wild and Domestic Rabbits in Australia

In 2020, Hepatitis E virus (HEV) was detected for the first time in Australian rabbits. To improve our understanding of the genetic diversity and distribution of the virus, 1635 rabbit liver samples from locations across Australia were screened via RT-qPCR for HEV. HEV genomes were amplified and sequenced from 48 positive samples. Furthermore, we tested 380 serum samples from 11 locations across Australia for antibodies against HEV. HEV was detected in rabbits from all states and territories, except the Northern Territory. Seroprevalence varied between locations (from 0% to 22%), demonstrating that HEV is widely distributed in rabbit populations across Australia. Phylogenetic analyses showed that Australian HEV sequences are genetically diverse and that HEV was likely introduced into Australia independently on several occasions. In summary, this study broadens our understanding of the genetic diversity of rabbit HEV globally and shows that the virus is endemic in both domestic and wild rabbit populations in Australia.

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

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

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.

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.

Rabbit Hemorrhagic Disease Virus Isolated from Diseased Alpine Musk Deer (Moschus sifanicus)

Rabbit hemorrhagic disease virus (RHDV) is the causative agent of rabbit hemorrhagic disease (RHD), and its infection results in mortality of 70-90% in farmed and wild rabbits. RHDV is thought to replicate strictly in rabbits. However, there are also reports showing that gene segments from the RHDV genome or antibodies against RHDV have been detected in other animals. Here, we report the detection and isolation of a RHDV from diseased Alpine musk deer (Moschussifanicus). The clinical manifestations in those deer were sudden death without clinical signs and hemorrhage in the internal organs. To identify the potential causative agents of the disease, we used sequence independent single primer amplification (SISPA) to detect gene segments from viruses in the tissue samples collected from the dead deer. From the obtained sequences, we identified some gene fragments showing very high nucleotide sequence similarity with RHDV genome. Furthermore, we identified caliciviral particles using an electron microscope in the samples. The new virus was designated as RHDV GS/YZ. We then designed primers based on the genome sequence of an RHDV strain CD/China to amplify and sequence the whole genome of the virus. The genome of the virus was determined to be 7437 nucleotides in length, sharing the highest genome sequence identity of 98.7% with a Chinese rabbit strain HB. The virus was assigned to the G2 genotype of RHDVs according to the phylogenetic analyses based on both the full-length genome and VP60 gene sequences. Animal experiments showed that GS/YZ infection in rabbits resulted in the macroscopic and microscopic lesions similar to that caused by the other RHDVs. This is the first report of RHDV isolated from Alpine musk deer, and our findings extended the epidemiology and host range of RHDV.

Adaptive changes in the genomes of wild rabbits after 16 years of viral epidemics

Since its introduction to control overabundant invasive European rabbits (Oryctolagus cuniculus), the highly virulent rabbit haemorrhagic disease virus (RHDV) has caused regular annual disease outbreaks in Australian rabbit populations. Although initially reducing rabbit abundance by 60%, continent-wide, experimental evidence has since indicated increased genetic resistance in wild rabbits that have experienced RHDV-driven selection. To identify genetic adaptations, which explain the increased resistance to this biocontrol virus, we investigated genome-wide SNP (single nucleotide polymorphism) allele frequency changes in a South Australian rabbit population that was sampled in 1996 (pre-RHD genomes) and after 16 years of RHDV outbreaks. We identified several SNPs with changed allele frequencies within or close to genes potentially important for increased RHD resistance. The identified genes are known to be involved in virus infections and immune reactions or had previously been identified as being differentially expressed in healthy versus acutely RHDV-infected rabbits. Furthermore, we show in a simulation study that the allele/genotype frequency changes cannot be explained by drift alone and that several candidate genes had also been identified as being associated with surviving RHD in a different Australian rabbit population. Our unique data set allowed us to identify candidate genes for RHDV resistance that have evolved under natural conditions, and over a time span that would not have been feasible in an experimental setting. Moreover, it provides a rare example of host genetic adaptations to virus-driven selection in response to a suddenly emerging infectious disease.

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.

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

The highly virulent rabbit hemorrhagic disease virus (RHDV) has been widely used in Australia and New Zealand since the mid-1990s to control wild rabbits, an invasive vertebrate pest in these countries. In January 2014, an exotic RHDV was detected in Australia, and 8 additional outbreaks were reported in both domestic and wild rabbits in the 15 months following its detection. Full-length genomic analysis revealed that this virus is a recombinant containing an RHDVa capsid gene and nonstructural genes most closely related to nonpathogenic rabbit caliciviruses. Nationwide monitoring efforts need to be expanded to assess if the increasing number of different RHDV variants circulating in the Australian environment will affect biological control of rabbits. At the same time, updated vaccines and vaccination protocols are urgently needed to protect pet and farmed rabbits from these novel rabbit caliciviruses.

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

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

•

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

•

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

•

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

•

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

•

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

Navigare necesse est! Research to understand our body and soul, to heal the patients, to find ourselves

The image of a researcher is a distinguished, enthusiastic doctor in a fresh, washed and ironed white lab coat working in a clinic or research-centre. The research is well-financed, supported by different scientific and/or economic firms, and the aim is to understand the human body and its physiologic processes in atomic level for getting the best, mostly very expensive, and sometimes uneasy medical treatment for the patient. Nowadays on top of the most modern sciences there is the specialist, who lives in an ivory tower and knows almost everything about diseases and sciences. Try to get off to this land!

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.

A strain-specific multiplex RT-PCR for Australian rabbit haemorrhagic disease viruses uncovers a new recombinant virus variant in rabbits and hares

Rabbit haemorrhagic disease virus (RHDV, or GI.1) is a calicivirus in the genus Lagovirus that has been widely utilized in Australia as a biological control agent for the management of overabundant wild European rabbit (Oryctolagus cuniculus) populations since 1996. Recently, two exotic incursions of pathogenic lagoviruses have been reported in Australia; GI.1a-Aus, previously called RHDVa-Aus, is a GI.1a virus detected in January 2014, and the novel lagovirus GI.2 (previously known as RHDV2). Furthermore, an additional GI.1a strain, GI.1a-K5 (also known as 08Q712), was released nationwide in March 2017 as a supplementary tool for wild rabbit management. To discriminate between these lagoviruses, a highly sensitive strain-specific multiplex RT-PCR assay was developed, which allows fast, cost-effective and sensitive detection of the four pathogenic lagoviruses currently known to be circulating in Australia. In addition, we developed a universal RT-qPCR assay to be used in conjunction with the multiplex assay that broadly detects all four viruses and facilitates quantification of viral RNA load in samples. These assays enable rapid detection, identification and quantification of pathogenic lagoviruses in the Australian context. Using these assays, a novel recombinant lagovirus was detected in rabbit tissue samples, which contained the non-structural genes of GI.1a-Aus and the structural genes of GI.2. This variant was also recovered from the liver of a European brown hare (Lepus europaeus). The impact of this novel recombinant on Australian wild lagomorph populations and its competitiveness in relation to circulating field strains, particularly GI.2, requires further studies.

Characterization of old RHDV strains by complete genome sequencing identifies a novel genetic group

Rabbit hemorrhagic disease (RHD) is a veterinary disease that affects the European rabbit and has a significant economic and ecological negative impact. In Portugal, rabbit hemorrhagic disease virus (RHDV) was reported in 1989 and still causes enzootic outbreaks. Several recombination events have been detected in RHDV strains, including in the first reported outbreak. Here we describe the occurrence of recombination in RHDV strains recovered from rabbit and Iberian hare samples collected in the mid-1990s in Portugal. Characterization of full genomic sequences revealed the existence of a single recombination breakpoint at the boundary of the non-structural and the structural encoding regions, further supporting the importance of this region as a recombination hotspot in lagoviruses. Phylogenetic analysis showed that in the structural region, the recombinant strains were similar to pathogenic G1 strains, but in the non-structural region they formed a new group that diverged ~13% from known strains. No further reports of such group exist, but this recombination event was also detected in an Iberian hare that was associated with the earliest species jump in RHDV. Our results highlight the importance of the characterization of full genomes to disclose RHDV evolution and show that lagoviruses’ diversity has been significantly undersampled.

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.

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.

An update on the rabbit hemorrhagic disease virus (RHDV) strains circulating in Portugal in the 1990s: earliest detection of G3-G5 and G6

Rabbit hemorrhagic disease virus (RHDV) causes devastating effects on European rabbit (Oryctolagus cuniculus) populations in the Iberian Peninsula. According to the information available, only genogroup 1 strains were circulating in Iberian wild rabbits until 2011; the antigenic variant G6 has been sporadically detected in rabbitries since 2007. Here, we show for the first time that G3-G5 strains were already present in mainland Portugal in 1998 and that G6 has been circulating since at least 1999. Moreover, we report a G3-G5 strain from the Azores collected in 1998, which is the likely ancestor of Azorean G3-G5like strains. These observations improve the current knowledge on RHDV epidemiology in the Iberian Peninsula and the Azores.

An in vivo system for directed experimental evolution of rabbit haemorrhagic disease virus

The calicivirus Rabbit haemorrhagic disease virus (RHDV) is widely used in Australia as a biocontrol agent to manage wild European rabbit (Oryctolagus cuniculus) populations. However, widespread herd immunity limits the effectiveness of the currently used strain, CAPM V-351. To overcome this, we developed an experimental platform for the selection and characterisation of novel RHDV strains. As RHDV does not replicate in cell culture, variant viruses were selected by serially passaging a highly virulent RHDV field isolate in immunologically naïve laboratory rabbits that were passively immunised 18–24 hours post-challenge with a neutralising monoclonal antibody. After seven passages, two amino acid substitutions in the P2 domain of the capsid protein became fixed within the virus population. Furthermore, a synonymous substitution within the coding sequence of the viral polymerase appeared and was also maintained in all subsequent passages. These findings demonstrate proof-of-concept that RHDV evolution can be experimentally manipulated to select for virus variants with altered phenotypes, in this case partial immune escape.

Experimental transmission of rabbit haemorrhagic disease virus (RHDV) from rabbit to wild mice (Mus spretus and Apodemus sylvaticus) under laboratory conditions

Rabbit haemorrhagic disease (RHD) is a highly lethal and contagious viral disease that produces haemorrhagic lesions in liver and lungs of domestic and wild rabbits (Oryctolagus cuniculus). This study investigates the transmission of RHDV from infected rabbits to mice, based on the detection of viral RNA. Sixteen wild mice (Mus spretus, n=12 and Apodemus sylvaticus, n=4) were put in contact with nine rabbits inoculated with RHDV. No mice died following exposure to RHDV-infected rabbits or developed macroscopic haemorrhagic lesions. On the fourth day of contact, RHDV was detected by RT-PCR in the faeces of three of the four mice killed and in the livers of two of them. Three days after contact period with the inoculated rabbits (7th day of the experiment), RHDV was detected by RT-PCR in 100% (n=4) of the faeces and 50% (n=2) of the livers of euthanized animals. Ten days after contact period (14th day of the experiment), RHDV was not detected in the faeces or liver from any of the mice euthanized. However, 64days after contact period, RHDV was detected in the faeces of one mouse (1 of 4). We demonstrate cross-species transmission of RHDV-RNA from rabbit to rodent and the capability of RHDV-RNA to persist in mice for at least 10days after contact, and potentially up to two months, although viral replication within the rodent and/or infectivity was not evaluated in the present study.

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.

Cross-validation to select Bayesian hierarchical models in phylogenetics

Background

Recent developments in Bayesian phylogenetic models have increased the range of inferences that can be drawn from molecular sequence data. Accordingly, model selection has become an important component of phylogenetic analysis. Methods of model selection generally consider the likelihood of the data under the model in question. In the context of Bayesian phylogenetics, the most common approach involves estimating the marginal likelihood, which is typically done by integrating the likelihood across model parameters, weighted by the prior. Although this method is accurate, it is sensitive to the presence of improper priors. We explored an alternative approach based on cross-validation that is widely used in evolutionary analysis. This involves comparing models according to their predictive performance.

Results

We analysed simulated data and a range of viral and bacterial data sets using a cross-validation approach to compare a variety of molecular clock and demographic models. Our results show that cross-validation can be effective in distinguishing between strict- and relaxed-clock models and in identifying demographic models that allow growth in population size over time. In most of our empirical data analyses, the model selected using cross-validation was able to match that selected using marginal-likelihood estimation. The accuracy of cross-validation appears to improve with longer sequence data, particularly when distinguishing between relaxed-clock models.

Conclusions

Cross-validation is a useful method for Bayesian phylogenetic model selection. This method can be readily implemented even when considering complex models where selecting an appropriate prior for all parameters may be difficult.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0688-y) contains supplementary material, which is available to authorized users.

Resolving the Origin of Rabbit Hemorrhagic Disease Virus: Insights from an Investigation of the Viral Stocks Released in Australia

To resolve the evolutionary history of rabbit hemorrhagic disease virus (RHDV), we performed a genomic analysis of the viral stocks imported and released as a biocontrol measure in Australia, as well as a global phylogenetic analysis. Importantly, conflicts were identified between the sequences determined here and those previously published that may have affected evolutionary rate estimates. By removing likely erroneous sequences, we show that RHDV emerged only shortly before its initial description in China.