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

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.

Comment on Shah et al. Genetic Characteristics and Phylogeographic Dynamics of Lagoviruses, 1988-2021. Viruses 2023, 15, 815

Shah and colleagues [...].

Epidemiological characterization and risk assessment of rabbit haemorrhagic disease virus 2 (RHDV2/b/GI.2) in the world

A novel variant of rabbit haemorrhagic disease virus, designated RHDV2/b/GI.2, was first discovered in France in 2010. Subsequently, RHDV2 rapidly spread to Africa, North America, Australia, and Asia. RHDV2 outbreaks have resulted in significant economic losses in the global rabbit industry and disrupted the balance of natural ecosystems. Our study investigated the seasonal characteristics of RHDV2 outbreaks using seasonal indices. RHDV2 is prone to causing significant outbreaks within domestic and wild rabbit populations during the spring season and is more likely to induce outbreaks within wild rabbit populations during late autumn in the Southern Hemisphere. Furthermore, based on outbreak data for domestic and wild rabbits and environmental variables, our study established two MaxEnt models to explore the relationship between RHDV2 outbreaks and the environmental factors and conducted outbreak risk predictions for RHDV2 in global domestic and wild rabbit populations. Both models demonstrated good predictive performance, with AUC values of 0.960 and 0.974, respectively. Road density, isothermality, and population density were identified as important variables in the outbreak of RHDV2 in domestic rabbits, while road density, normalized difference vegetation index, and mean annual solar radiation were considered key variables in the outbreak of RHDV2 in wild rabbits. The environmental factors associated with RHDV2 outbreaks identified in our study and the outbreak risk prediction maps generated in our study will aid in the formulation of appropriate RHDV2 control measures to reduce the risk of morbidity in domestic and wild rabbits.

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.

A new HaCV-EBHSV recombinant lagovirus circulating in European brown hares (Lepus europaeus) from Catalonia, Spain

In 2020/2021, several European brown hare syndrome virus (EBHSV) outbreaks were recorded in European hares (Lepus europaeus) from Catalonia, Spain. Recombination analysis combined with phylogenetic reconstruction and estimation of genetic distances of the complete coding sequences revealed that 5 strains were recombinants. The recombination breakpoint is located within the non-structural protein 2C-like RNA helicase (nucleotide position ~ 1889). For the genomic fragment upstream of the breakpoint, a non-pathogenic EBHSV-related strain (hare calicivirus, HaCV; GII.2) was the most closely related sequence; for the rest of the genome, the most similar strains were the European brown hare syndrome virus (EBHSV) strains recovered from the same 2020/2021 outbreaks, suggesting a recent origin. While the functional impact of the atypical recombination breakpoint remains undetermined, the novel recombinant strain was detected in different European brown hare populations from Catalonia, located 20-100 km apart, and seems to have caused a fatal disease both in juvenile and adult animals, confirming its viability and ability to spread and establish infection. This is the first report of a recombination event involving HaCV and EBHSV and, despite the recombination with a non-pathogenic strain, it appears to be associated with mortality in European brown hares, which warrants close monitoring.

Recombination between non-structural and structural genes as a mechanism of selection in lagoviruses: The evolutionary dead-end of an RHDV2 isolated from European hare

The genus Lagovirus, belonging to the family Caliciviridae, emerged around the 1980s. It includes highly pathogenic species, rabbit hemorrhagic disease virus (RHDV/GI.1) and European brown hare syndrome virus (EBHSV/GII.1), which cause fatal hepatitis, and nonpathogenic viruses with enteric tropism, rabbit calicivirus (RCV/GI.3,4) and hare calicivirus (HaCV/GII.2). Lagoviruses have evolved along two independent genetic lineages: GI (RHDV and RCV) in rabbits and GII (EBHSV and HaCV) in hares. To be emphasized is that genomes of lagoviruses, like other caliciviruses, are highly conserved at RdRp-VP60 junctions, favoring intergenotypic recombination events at this point. The recombination between an RCV (genotype GI.3), donor of non-structural (NS) genes, and an unknown virus, donor of structural (S) genes, likely led to the emergence of a new lagovirus in the European rabbit, called RHDV type 2 (GI.2), identified in Europe in 2010. New RHDV2 intergenotypic recombinants isolated in rabbits in Europe and Australia originated from similar events between RHDV2 (GI.2) and RHDV (GI.1) or RCV (GI.3,4). RHDV2 (GI.2) rapidly spread worldwide, replacing RHDV and showing several lagomorph species as secondary hosts. The recombination events in RHDV2 viruses have led to a number of viruses with very different combinations of NS and S genes. Recombinant RHDV2 with NS genes from hare lineage (GII) was recently identified in the European hare. This study investigated the first RHDV2 (GI.2) identified in Italy in European hare (RHDV2_Bg12), demonstrating that it was a new virus that originated from the recombination between RHDV2, as an S-gene donor and a hare lagovirus, not yet identified but presumably nonpathogenic, as an NS gene donor. When rabbits were inoculated with RHDV2_Bg12, neither deaths nor seroconversions were recorded, demonstrating that RHDV2_Bg12 cannot infect the rabbit. Furthermore, despite intensive and continuous field surveillance, RHDV2_Bg12 has never again been identified in either hares or rabbits in Italy or elsewhere. This result showed that the host specificity of lagoviruses can depend not only on S genes, as expected until today, but potentially also on some species-specific NS gene sequences. Therefore, because RHDV2 (GI.2) infects several lagomorphs, which in turn probably harbor several specific nonpathogenic lagoviruses, the possibility of new speciation, especially in those other than rabbits, is real. RHDV2 Bg_12 demonstrated this, although the attempt apparently failed.

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.

Diagnosis of a Rabbit Hemorrhagic Disease Virus 2 (RHDV2) and the Humoral Immune Protection Effect of VP60 Vaccine

Rabbit hemorrhagic disease (RHD) is known as rabbit plague and hemorrhagic pneumonia. It is an acute, septic, and highly fatal infectious disease caused by the Lagovirus rabbit hemorrhagic disease virus (RHDV) in the family Caliciviridae that infects wild and domestic rabbits and hares (lagomorphs). At present, RHDV2 has caused huge economic losses to the commercial rabbit trade and led to a decline in the number of wild lagomorphs worldwide. We performed a necropsy and pathological observations on five dead rabbits on a rabbit farm in Tai'an, China. The results were highly similar to the clinical and pathological changes of typical RHD. RHDV2 strain was isolated and identified by RT-PCR, and partial gene sequencing and genetic evolution analysis were carried out. There were significant differences in genetic characteristics and antigenicity between RHDV2 and classical RHDV strain, and the vaccine prepared with the RHDV strain cannot effectively prevent rabbit infection with RHDV2. Therefore, we evaluated the protective efficacy of a novel rabbit hemorrhagic virus baculovirus vector inactivated vaccine (VP60) in clinical application by animal regression experiment. The result showed that VP60 could effectively induce humoral immunity in rabbits. The vaccine itself had no significant effect on the health status of rabbits. This study suggested that the clinical application of VP60 may provide new ideas for preventing the spread of RHD2.

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.

Confirmation of the Rabbit Hemorrhagic Disease Virus Type 2 (GI.2) Circulation in North Africa

Rabbit hemorrhagic disease (RHD) is a highly contagious viral disease that causes fatal acute hepatitis in domestic and wild lagomorphs. It has taken on major economic importance in countries like Morocco. In addition to the classical virus (RHDV), a novel emerged genotype (RHDV2) is circulating, especially in the north shore of the Mediterranean basin since 2010. Many small animal farmers reported clinical cases from several rabbitries in Agadir (Morocco) despite systematic vaccination against the RHDV. The main objective was to characterize the current RHDV strains circulating in the studied area to help to choose an adequate vaccine. For that, we extracted viral RNA from rabbit livers, carried out the PCR analyses, and we sequenced the viral structural capsid protein (VP60) of the RHDV. The phylogenetic analysis results allowed us to state that the novel genotype (RHDV2) is circulating in the studied geographical area, and to characterize the isolated sequences. As a conclusion, we recommend updating RHD epidemiological relating data and reviewing the vaccine protocols by both targeting RHDV (GI.1) and RHDV2 (GI.2) in any future preventive program.

Nucleolin interacts with the rabbit hemorrhagic disease virus replicase RdRp, nonstructural proteins p16 and p23, playing a role in virus replication

Rabbit hemorrhagic disease virus (RHDV) is a member of the Caliciviridae family and cannot be propagated in vitro, which has impeded the progress of investigating its replication mechanism. Construction of an RHDV replicon system has recently provided a platform for exploring RHDV replication in host cells. Here, aided by this replicon system and using two-step affinity purification, we purified the RHDV replicase and identified its associated host factors. We identified rabbit nucleolin (NCL) as a physical link, which mediating the interaction between other RNA-dependent RNA polymerase (RdRp)-related host proteins and the viral replicase RdRp. We found that the overexpression or knockdown of NCL significantly increased or severely impaired RHDV replication in RK-13 ​cells, respectively. NCL was identified to directly interact with RHDV RdRp, p16, and p23. Furthermore, NCL knockdown severely impaired the binding of RdRp to RdRp-related host factors. Collectively, these results indicate that the host protein NCL is essential for RHDV replication and acts as a physical link between viral replicase and host proteins.

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.

Comparative susceptibility of eastern cottontails and New Zealand white rabbits to classical rabbit haemorrhagic disease virus (RHDV) and RHDV2

Rabbit haemorrhagic disease virus (RHDV) is associated with high morbidity and mortality in the European rabbit (Oryctolagus cuniculus). In 2010, a genetically distinct RHDV named RHDV2 emerged in Europe and spread to many other regions, including North America in 2016. Prior to this study it was unknown if eastern cottontails (ECT(s); Sylvilagus floridanus), one of the most common wild lagomorphs in the United States, were susceptible to RHDV2. In this study, 10 wild-caught ECTs and 10 New Zealand white rabbits (NZWR(s); O. cuniculus) were each inoculated orally with either RHDV (RHDVa/GI.1a; n = 5 per species) or RHDV2 (a recombinant GI.1bP-GI.2; n = 5 per species) and monitored for the development of disease. Three of the five ECTs that were infected with RHDV2 developed disease consistent with RHD and died at 4 and 6 days post-inoculation (DPI). The RHDV major capsid protein/antigen (VP60) was detected in the livers of three ECTs infected with RHDV2, but none was detected in the ECTs infected with RHDV. Additionally, RHD viral RNA was detected in the liver, spleen, intestine and blood of ECTs infected with RHDV2, but not in the ECTs infected with RHDV. RHD viral RNA was detected in urine, oral swabs and rectal swabs in at least two of five ECTs infected with RHDV2. One ECT inoculated with RHDV2 seroconverted and developed a high antibody titre by the end of the experimental period (21 DPI). ECTs inoculated with the classic RHDV did not seroconvert. In comparison, NZWRs inoculated with RHDV2 exhibited high mortality (five of five) at 2 DPI and four of five NZWRs inoculated with RHDV either died or were euthanized at 2 DPI indicating both of these viruses were highly pathogenic to this species. This experiment indicates that ECTs are susceptible to RHDV2 and can shed viral RNA, thereby suggesting this species could be involved in the epidemiology of this virus.

Hemoglobin subunit beta interacts with the capsid, RdRp and VPg proteins, and antagonizes the replication of rabbit hemorrhagic disease virus

Rabbit hemorrhagic disease virus (RHDV) causes a highly contagious disease in rabbits that is associated with high mortality. Because of the lack of a suitable cell culture system for RHDV, its pathogenic mechanism and replication remain unclear. This study found that the expression level of host protein rabbit hemoglobin subunit beta (HBB) was significantly downregulated in RHDV-infected cells. To investigate the role of HBB in RHDV replication, small interfering RNAs for HBB and HBB eukaryotic expression plasmids were used to change the expression level of HBB in RK-13 cells and the results showed that the RHDV replication level was negatively correlated with the expression level of HBB. It was also verified that HBB inhibited RHDV replication using constructed HBB stable overexpression cell lines and HBB knockout cell lines. The interaction of HBB with viral capsid protein VP60, replicase RdRp, and VPg protein was confirmed, as was the activation of the expression of interferon γ by HBB. The results of this study indicated that HBB may be an important host protein in host resistance to RHDV infection.

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.

Full-genome sequencing of German rabbit haemorrhagic disease virus uncovers recombination between RHDV (GI.2) and EBHSV (GII.1)

Rabbit haemorrhagic disease virus (RHDV; genotypes GI.1 and GI.2) and European brown hare syndrome virus (EBHSV; genotype GII.1) are caliciviruses belonging to the genus Lagovirus. These viruses pose a serious threat to wild and domestic rabbit and hare populations around the world. In recent years, an expanding genetic diversity has been described within the genus, with recombination events occurring between the different genotypes. Here, we generated and analysed 56 full-genome sequences of RHDV and EBHSV from rabbit and hare livers, collected in Germany between the years 2013 and 2020. We could show that genotype Gl.2 (RHDV-2) almost entirely replaced Gl.1 (classical RHDV) in the German rabbit population. However, GI.1 is still present in Germany and has to be included into disease control and vaccination strategies. Three recombinant strains were identified from rabbit samples that contain the structural genes of genotype Gl.2 and the non-structural genes of genotype Gl.1b. Of special interest is the finding that sequences from two hare samples showed recombination events between structural genes of RHDV Gl.2 and non-structural genes of EBHSV GII.1, a recombination between different genogroups that has not been described before. These findings lead to the assumption that also a recombination of the non-structural genes of RHDV Gl.2 with the structural genes of EBHSV Gll.1 might be possible and therefore increase the potential genetic variability of lagoviruses immensely. Our findings underline the importance of whole genome analysis with next-generation sequencing technology as one of new tools now available for in-depth studies that allow in depth molecular epidemiology with continuous monitoring of the genetic variability of viruses that would otherwise likely stay undetected if only routine diagnostic assays are used.

RPS5 interacts with the rabbit hemorrhagic disease virus 3' extremities region and plays a role in virus replication

Rabbit hemorrhagic disease virus (RHDV), a member of Caliciviridae family, causes a highly contagious disease in rabbits. The RHDV replication mechanism is poorly understood due to the lack of a suitable culture system in vitro. This study identified RHDV 5' and 3' extremities (Ex) RNA binding proteins from the rabbit kidney cell line RK-13 based on a pull-down assay by applying a tRNA scaffold streptavidin aptamer. Using mass spectrometry (MS), several host proteins were discovered which interact with RHDV 5' and 3' Ex RNA. The ribosomal protein S5 (RPS5) was shown to interact with RHDV 3' Ex RNA directly by RNA-pulldown and confocal microscopy. To further investigate the role of RPS5 in RHDV replication, small interfering RNAs for RPS5 and RPS5 eukaryotic expression plasmids were used to change the expression level of RPS5 in RK-13 cells and the results showed that the RHDV replication and translation levels were positively correlated with the expression level of RPS5. It was also verified that RPS5 promoted RHDV replication by constructing RPS5 stable overexpression cell lines and RPS5 knockdown cell lines. In summary, it has been identified that RPS5 interacted with the RHDV 3' Ex RNA region and played a role in virus replication. These results will help to understand the mechanism of RHDV replication.

Novel hepaci- and pegi-like viruses in native Australian wildlife and non-human primates

The Flaviviridae family of positive-sense RNA viruses contains important pathogens of humans and other animals, including Zika virus, dengue virus, and hepatitis C virus. The Flaviviridae are currently divided into four genera-Hepacivirus, Pegivirus, Pestivirus, and Flavivirus-each with a diverse host range. Members of the genus Hepacivirus are associated with an array of animal species, including humans, non-human primates, other mammalian species, as well as birds and fish, while the closely related pegiviruses have been identified in a variety of mammalian taxa, also including humans. Using a combination of total RNA and whole-genome sequencing we identified four novel hepaci-like viruses and one novel variant of a known hepacivirus in five species of Australian wildlife. The hosts infected comprised native Australian marsupials and birds, as well as a native gecko (Gehyra lauta). From these data we identified a distinct marsupial clade of hepaci-like viruses that also included an engorged Ixodes holocyclus tick collected while feeding on Australian long-nosed bandicoots (Perameles nasuta). Distinct lineages of hepaci-like viruses associated with geckos and birds were also identified. By mining the SRA database we similarly identified three new hepaci-like viruses from avian and primate hosts, as well as two novel pegi-like viruses associated with primates. The phylogenetic history of the hepaci- and pegi-like viruses as a whole, combined with co-phylogenetic analysis, provided support for virus-host co-divergence over the course of vertebrate evolution, although with frequent cross-species virus transmission. Overall, our work highlights the diversity of the Hepacivirus and Pegivirus genera as well as the uncertain phylogenetic distinction between.

Rapid detection of H146-like goose calicivirus using real-time RT-PCR with a Taqman minor groove binder probe

H146-like goose-origin calicivirus (H146-like GCV) is a novel Caliciviridae family member in the Sanovirus genus that was associated with gosling growth retardation syndrome growth retardation syndrome complicated by visceral urate deposition. However, there is no accurate and high throughput real-time quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR) available for the rapid and highly sensitive identification of H146-like GCV. In this study, a pair of specific primers and a TaqMan minor groove binder (MGB) probe were designed based on a conserved region in the nonstructural (NS) gene sequence. The TaqMan-MGB probe-based one-step qRT-PCR assay was capable of detecting quite low number of targeting nucleic acid as low as 5.07 copies/μL and had excellent intra-assay and inter-assay repeatability with the coefficient of variation (CV) value from 0.558% to 1.293%. The assay was highly specific for H146-like GCV, without cross-reactions with other non-targeted goose-origin viruses, and 62 suspicious tissue samples infected with H146-like GCV from different regions of Fujian Province were used in this study to verify the feasibility and effectiveness of this assay in clinical diagnosis. The results indicated that our assay for the diagnosis and quantification of H146-like GCV was highly sensitive and specific, and should provide a reliable real-time tool for epidemiological and pathogenetic study of H146-like GCV infection, enabling researchers to better understand the epidemiology and clinical presentation of this disease.

Retrospective Analysis Shows That Most RHDV GI.1 Strains Circulating Since the Late 1990s in France and Sweden Were Recombinant GI.3P-GI.1d Strains

Recombination is one of the major sources of genetic variation in viruses. RNA viruses, such as rabbit hemorrhagic disease virus (RHDV), are among the viruses with the highest recombination rates. Several recombination events have been described for RHDV, mostly as a consequence of their genomic architecture. Here, we undertook phylogenetic and recombination analyses of French and Swedish RHDV strains from 1994 to 2016 and uncovered a new intergenotypic recombination event. This event occurred in the late 1990s/early 2000s and involved nonpathogenic GI.3 strains as donors for the nonstructural part of the genome of these recombinants, while pathogenic GI.1d strains contributed to the structural part. These GI.3P-GI.1d recombinant strains did not entirely replace GI.1d (nonrecombinant) strains, but became the dominant strains in France and Sweden, likely due to a fitness advantage associated with this genomic architecture. GI.3P-GI.1d (P stands for polymerase) strains persisted until 2013 and 2016 in Sweden and France, respectively, and cocirculated with the new genotype GI.2 in France. Since strains from the first GI.2 outbreaks were GI.3P-GI.2, we hypothesize that GI.3P-GI.1d could be the parental strain. Our results confirm the outstanding recombination ability of RHDV and its importance in the evolution of lagoviruses, which was only revealed by studying complete genomic sequences.

Comparative Analysis of RNA Virome Composition in Rabbits and Associated Ectoparasites

Ectoparasites play an important role in the transmission of many vertebrate-infecting viruses, including Zika and dengue viruses. Although it is becoming increasingly clear that invertebrate species harbor substantial virus diversity, it is unclear how many of the viruses carried by invertebrates have the potential to infect vertebrate species. We used the European rabbit (Oryctolagus cuniculus) as a model species to compare virome compositions in a vertebrate host and known associated ectoparasite mechanical vectors, in this case, fleas and blowflies. In particular, we aimed to infer the extent of viral transfer between these distinct types of host. Our analysis revealed that despite extensive viral diversity in both rabbits and associated ectoparasites, and the close interaction of these vertebrate and invertebrate species, biological viral transmission from ectoparasites to vertebrate species is rare. We did, however, find evidence to support the idea of a role of blowflies in transmitting viruses without active replication in the insect.

Ectoparasites play an important role in virus transmission among vertebrates. Little, however, is known about the nature of those viruses that pass between invertebrates and vertebrates. In Australia, flies and fleas support the mechanical transmission of two viral biological controls against wild rabbits—rabbit hemorrhagic disease virus (RHDV) and myxoma virus. We compared virome compositions in rabbits and these ectoparasites, sequencing total RNA from multiple tissues and gut contents of wild rabbits, fleas collected from these rabbits, and flies trapped sympatrically. Meta-transcriptomic analyses identified 50 novel viruses from multiple RNA virus families. Rabbits and their ectoparasites were characterized by markedly different viromes, with virus abundance greatest in flies. Although viral contigs from six virus families/groups were found in both rabbits and ectoparasites, they clustered in distinct host-dependent lineages. A novel calicivirus and a picornavirus detected in rabbit cecal content were vertebrate specific; the newly detected calicivirus was distinct from known rabbit caliciviruses, while the picornavirus clustered with sapeloviruses. Several picobirnaviruses were also identified that fell in diverse phylogenetic positions, compatible with the idea that they are associated with bacteria. Further comparative analysis revealed that the remaining viruses found in rabbits, and all those from ectoparasites, were likely associated with invertebrates, plants, and coinfecting endosymbionts. While no full genomes of vertebrate-associated viruses were detected in ectoparasites, small numbers of reads from rabbit astrovirus, RHDV, and other lagoviruses were present in flies. This supports a role for flies in the mechanical transmission of RHDV, while their involvement in astrovirus transmission merits additional exploration.

IMPORTANCE Ectoparasites play an important role in the transmission of many vertebrate-infecting viruses, including Zika and dengue viruses. Although it is becoming increasingly clear that invertebrate species harbor substantial virus diversity, it is unclear how many of the viruses carried by invertebrates have the potential to infect vertebrate species. We used the European rabbit (Oryctolagus cuniculus) as a model species to compare virome compositions in a vertebrate host and known associated ectoparasite mechanical vectors, in this case, fleas and blowflies. In particular, we aimed to infer the extent of viral transfer between these distinct types of host. Our analysis revealed that despite extensive viral diversity in both rabbits and associated ectoparasites, and the close interaction of these vertebrate and invertebrate species, biological viral transmission from ectoparasites to vertebrate species is rare. We did, however, find evidence to support the idea of a role of blowflies in transmitting viruses without active replication in the insect.

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.

Carrion flies as sentinels for monitoring lagovirus activity in Australia

Lagoviruses are an essential tool for managing wild rabbit populations in Australia. Our understanding of lagovirus epidemiology in Australia currently depends on members of the public submitting liver samples from dead lagomorphs (i.e. rabbits and hares) through a monitoring program called Rabbitscan. However, many wild lagomorphs die in inaccessible locations or are scavenged before sampling can occur, leading to considerable sampling bias. In this study, we screened field-caught carrion flies for the presence of lagoviruses to monitor virus circulation patterns in the landscape, with an aim to establish a less biased epidemiological surveillance tool. Carrion flies were collected from two study sites over a 22-month period and these samples were used to optimize and validate molecular testing methods in this sample type for the currently circulating lagovirus variants. Virus was clearly detectable in field-caught carrion flies using optimized SYBR-green RT-qPCR and RT-PCR assays. However, variant identification was frequently hindered by the low virus loads present in carrion fly samples and spurious RT-PCR amplification. This was overcome by frequent sampling, which effectively acts as replicate sampling to verify inconclusive results. There was generally good correlation between virus detections and variant identification in carrion flies and in samples recovered from wild lagomorphs. The methods reported here provide an additional surveillance tool to monitor lagovirus spread and circulation at a landscape scale, which in turn can help to guide more effective rabbit management programs.

Molecular methods in detection and epidemiologic studies of rabbit and hare viruses: a review

Various PCR-based assays for rabbit viruses have gradually replaced traditional virologic assays, such as virus isolation, because they offer high-throughput analysis, better test sensitivity and specificity, and allow vaccine and wild-type virus strains to be fully typed and differentiated. In addition, PCR is irreplaceable in the detection of uncultivable or fastidious rabbit pathogens or those occurring in low quantity in a tested sample. We provide herein an overview of the current state of the art in the molecular detection of lagomorph viral pathogens along with details of their targeted gene or nucleic acid sequence and recommendations for their application. Apart from the nucleic acids-based methods used for identification and comprehensive typing of rabbit viruses, novel methods such as microarray, next-generation sequencing, and mass spectrometry (MALDI-TOF MS) could also be employed given that they offer greater throughput in sample screening for viral pathogens. Molecular methods should be provided with an appropriate set of controls, including an internal amplification control, to confirm the validity of the results obtained.

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.

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.