# Large-scale lagovirus disease outbreaks in European brown hares (Lepus europaeus) in France caused by RHDV2 strains spatially shared with rabbits (Oryctolagus cuniculus)

Ghislaine Le Gall-Reculé, Evelyne Lemaitre, Stéphane Bertagnoli, Céline Hubert, Sokunthea Top, Stéphane Marchandeau, Jean-Sébastien Guitton

► To cite this version:

Ghislaine Le Gall-Reculé, Evelyne Lemaitre, Stéphane Bertagnoli, Céline Hubert, Sokunthea Top, et al.. Large-scale lagovirus disease outbreaks in European brown hares (Lepus europaeus) in France caused by RHDV2 strains spatially shared with rabbits (Oryctolagus cuniculus). Veterinary Research, 2017, 48 (1), pp.70. 10.1186/s13567-017-0473-y. hal-01625773

HAL Id: hal-01625773

https://hal.science/hal-01625773

Submitted on 28 Oct 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

Open Access

<figure><img src="image_3.png" /></figure>

CrossMark

# Large-scale lagovirus disease outbreaks in European brown hares (*Lepus europaeus*) in France caused by RHDV2 strains spatially shared with rabbits (*Oryctolagus cuniculus*)

Ghislaine Le Gall-Reculé<fnref n="1" />,²*, Évelyne Lemaitre<fnref n="1" />,², Stéphane Bertagnoli³, Céline Hubert<fnref n="1" />,², Sokunthea Top³, Anouk Decors⁴, Stéphane Marchandeau⁵ and Jean-Sébastien Guitton⁵

## Abstract

Rabbit haemorrhagic disease virus (RHDV) is a lagovirus that causes rabbit haemorrhagic disease (RHD) in European rabbits (*Oryctolagus cuniculus*). In 2010, a new genotype called RHDV2 emerged in France. It exhibits a larger host range than classical RHDV strains by sporadically infecting different hare species, including the European hare (*Lepus europaeus*). Phylogenetic analyses revealed that closely related RHDV2 strains circulate locally in both hares and rabbits, and therefore that RHDV2 strains infecting hares do not belong to a lineage that has evolved only in this species. We showed that RHDV2 is widely distributed in France and that it was responsible for more than a third of cases of lagovirus disease in European hare populations in 2015. The oldest RHDV2 positive hare was sampled in November 2013 and we reported two hares co-infected by EBHSV and RHDV2. All together, our results raise important epide miological and evolutionary issues. In particular, along with the potential emergence of recombinant EBHSV/RHDV2 strains in hares, the enlargement of the host range changes the host population structure of RHDV2 and may alter the impact of the virus on rabbit and hare populations.

## Introduction

Rabbit haemorrhagic disease (RHD) is a highly infectious and often fatal viral disease of European rabbits (*Oryctolagus cuniculus*). The disease was first reported in China in 1984 and subsequently in many countries throughout the world [1]. The aetiological agent is a non-enveloped single-stranded positive-sense RNA virus belonging to the genus *Lagovirus* of the family *Caliciviridae* and called rabbit haemorrhagic disease virus (RHDV). RHD induces high mortality in both domestic and wild rabbit popula tions. Clinical signs and histopathological characteristics include severe necrotizing hepatitis and massive dissemi nated intravascular coagulation [2].

Until 2010, phylogenetic analyses revealed the progres sive emergence of distinct but closely related variants of RHDV [3–5]. In 2010, a new genotype of RHDV, named RHDV2, was identified in France with distinct patho genic, genetic and antigenic profiles [6, 7]. RHDV2 rap idly spread throughout Europe, causing significant losses in both domestic and wild rabbits [6, 8–11], and was then detected in different islands, in Australia, in North ern America and in Africa [12–18]. In Western Europe, RHDV2 is responsible for almost all cases of RHD in both domestic and wild rabbit, cases of RHD caused by former RHDV strains being rare [6, 9, 13, 19, 20]. RHDV2 quickly evolved and recombination with rabbit lagovi ruses that co-circulate in Europe has been evidenced in the Iberian Peninsula [21, 22]. These recombinations include the region of the genome encoding the structural VP60 and VP10 proteins of RHDV2 with the upstream region of the genome encoding the non-structural

<footnote>*Correspondence: ghislaine.legall-recule@anses.fr</footnote>

<footnote n="1">¹ Avian and Rabbit Virology, Immunology and Parasitology Unit, Anses, French Agency for Food, Environmental and Occupational Health & Safety, Ploufragan-Plouzané Laboratory, BP 53, 22440 Ploufragan, France
Full list of author information is available at the end of the article</footnote>
proteins, either from the RHDV genogroup G1 or from non-pathogenic lagoviruses.

RHDV2 exhibits a broader host range than classical RHDV by infecting not only rabbits but also different hare species (*Lepus capensis mediterraneus*, *Lepus corsicanus*, *Lepus europaeus*, *Lepus timidus*) [23–27]. RHDV2-infected hares show very similar clinical signs to those induced by European brown hare syndrome (EBHS): hyperaemic trachea sometimes containing uncoagulated blood, hepatitis necrosis, splenomegaly and congestion of the other organs and tissues [23, 25, 26]. EBHS was first described in 1980 [28] and is caused by another lagovirus called EBHSV, which is known to infect and cause disease in several hare species (*L. europaeus*, *L. timidus* and *L. corsicanus*). In the European hare (*L. europaeus*), some cases of hares showing lesions similar to EBHS (EBHS-like disease) caused by RHDV2 were reported in Italy in 2012, in Spain in 2014 and in Australia in 2016 [24, 26]. To date it is unclear whether RHDV2 may cause outbreaks in hare populations. Indeed, Velarde et al. [26] described sporadic cases in Italy (one hare) and Spain (two hares) and the data recorded in Australia do not indicate whether the five positive hares correspond to outbreaks or to spillover events [24].

Assessing the spreading level of RHDV2 into hare species, and especially *L. europaeus*, is an important issue because it may have significant epidemiological and evolutionary implications in areas where European hares are widespread like France for instance. For example, a large-scale spread into hare populations would increase the probability of recombination events between RHDV2 and EBHSV. It would also change the spatial structure of RHDV2 host population, which is known to influence the evolution of pathogens [29, 30]. Yet rabbit populations are patchily distributed [31, 32], while hare populations are more continuously distributed [33].

Therefore, the present study aims at detecting RHDV2 in European hares which died of lagovirus disease in France and at assessing the relative proportion of cases due to EBHSV and RHDV2. We searched for RHDV2 in a sample of 50 dead hares collected between 2013 and

early 2015 and compared these strains to those involved in RHD outbreaks in rabbits during the same period. We then carried out a large scale molecular epidemiological study to estimate the proportion of lagovirus disease caused by RHDV2 in 2015 in France. This latter analysis was performed on 208 liver samples taken on hares that died from an EBHS-like disease in 2015.

## Materials and methods

### Detection and molecular characterization of RHDV2 in hares

#### European hare samples

We randomly sampled 50 European hares collected by the SAGIR network between January 2013 and February 2015 displaying lesions compatible with EBHS at necropsy (Table 1). The livers were stored at -20 °C before being sent to the Anses' laboratory for further molecular analyses.

SAGIR is a generalist incident-based surveillance network for epidemiological vigilance towards wildlife diseases dealing with early detection and early warning [34]. The vigilance relies on a diagnostic process, based on a transdisciplinary approach (epidemiology, ecology, toxicology and pathology). A gross pathologic examination is performed on each submitted animal. According to field and epidemiological clues as well as a gross pathologic and clinical picture, additional tests (parasitology, bacteriology, virology, histology) helpful to investigate the aetiology of death are prioritized. Calicivirus investigation relies therefore on epidemiological or/and ongoing process type (acute) or/and evocative lesion.

#### Search for lagoviruses by RT-PCR-sequencing

Each frozen liver sample was thawed and total RNA was extracted from 100 µL of the obtained exudate using the NucleoSpin® RNA kit (Macherey-Nagel) according to the manufacturer's instructions. RNAs were analyzed to obtain molecular data using the RT-PCR U38/EBHS9 [35] that can amplify a part of the capsid protein gene of both EBHSV and RHDV/RHDV2. Reverse transcriptions were performed using oligo(dT) as primers (Invitrogen)

<table><caption>**Table 1** Summary of samples, methods and results on hares</caption><thead><tr><td>Objective</td><td>Sample size</td><td>Collection period</td><td>Methods used</td><td>Virus strains</td></tr></thead><tbody><tr><td>Detection and molecular characterization of RHDV2 in hares</td><td>50<sup>a</sup></td><td>01/2013 to 02/2015</td><td>RT-PCRs U38/EBHS9 and VP60, search for recombinant strains (RT-PCR 1U/1L), sequencing</td><td>7 RHDV2 (of which 3 recombinant strains and 1 coinfection RHDV2/EBHSV)</td></tr><tr><td>Relative proportion of RHDV2 and EBHSV in lagovirus infections</td><td>208<sup>a</sup></td><td>2015</td><td>rRT-PCR RHDV2 and rRT-PCR EBHSV</td><td>32 RHDV2, 54 EBHSV and 1 coinfection RHDV2/EBHSV</td></tr></tbody></table>

Complete information on the lagovirus positive samples is reported in the Additional files 1 and 2.

<footnote>a Ten samples have been used for both objectives.</footnote>
and SuperScript™ II Reverse Transcriptase (Invitrogen). The cDNA were amplified using the AmpliTaq Gold DNA polymerase (Applera Applied Biosystems) (10 min at 94 °C, followed by 30 cycles at 94 °C for 30 s, 57.6 °C for 30 s and 72 °C for 45 s, then 10 min at 72 °C) and visualized by electrophoresis on agarose gel.

When they existed, PCR products were purified using MinElute™ PCR Purification Kit (QIAGEN). DNA sequences were determined in both directions using the PCR primers and Big Dye Terminator v3.1 (Life Technologies) as recommended by the manufacturer, and analyzed with an “ABI Prism 3130 Genetic Analyzer” (Applied Biosystems). Nucleotide sequences were compiled using Vector NTI Advance 11.5 (Invitrogen). In order to genotype the viruses, the sequences were aligned against the nucleotide sequences present in databases using the standard nucleotide BLAST (blastn) program available in the National Center for Biotechnology Information (NCBI) web BLAST service.

## Amplification and sequencing of the capsid protein genes of RHDV2

The complete sequence of the gene encoding the capsid protein of seven RHDV2 were obtained following specific full-length gene PCR amplification using the Expand™ High Fidelity enzyme (Roche Applied Science) and the primers 12U and 15L [6]. For one RHDV2 (E14-73), the capsid protein gene amplification was obtained using the primers 12U2 (5'-GGATCGTCTCGGTAGTACC, 5263–5282) and Calici7255-Rev (5'-CTACACTAGCAT CATTATGCAT, 7268–7247). The cycle conditions were those indicated by the manufacturer, with an annealing temperature of 56.8 and 54 °C, respectively, during 30 s.

The amplified DNA were purified and sequenced as mentioned above using the PCR primers as well as several primers designed in the inner region of the templates (primers sequences available upon request). The sequences were compiled using Vector NTI Advance 11.5 (Invitrogen).

## Phylogenetic analysis

Phylogenetic relationships were inferred using the VP60 gene sequences of seven RHDV2 characterized in brown hares and rabbit lagovirus VP60 sequences available in databases. We also included the VP60 gene sequences of nine French RHDV2 circulating in wild rabbits and farms between January 2013 and November 2014 (Additional file 1). The sequence of the French reference EBHSV strain EBHSV-GD (Z69620) was used as an outgroup to root the trees. Phylogenetic analyses were conducted in MEGA software version 5 [36]. The neighbor-joining method was implemented with the pairwise deletion option and based on the Kimura 2-parameter

model including transition and transversion substitutions. Codon positions included were 1st + 2nd + 3rd. Search for the trees was implemented with the close-neighbor-interchange algorithm. Reliability of the trees was assessed by bootstrap with 1000 replicates. Pairwise nucleotide distance comparisons based on the *p*-distance model were conducted in MEGA5.

## Nucleotide VP60 sequence accession numbers

Nucleotide VP60 sequences of six RHDV2 characterized in French hares (E13-37, E14-03, E14-40, E14-45, E14-115 and E15-12) and nine RHDV2 characterized in French rabbits (13–22, 13–69, 13–71, 13–100, 13–122, 13–165, 14–02, 14–59 and 14–114) are available in databases under the EMBL/GenBank Accession Numbers LT168833 to LT168847. The RHDV2 E14-73 VP60 sequence is available in databases under the Accession Number LT549473.

## Search for RHDV2 recombinant strains

To determine whether the characterized RHDV2 strains were recombinant or not with RHDV G1 or RCV-A1-like viruses as are some RHDV2 from the Iberian Peninsula [21], we genotyped a 518 bp-long region of the genome upstream of the VP60 gene and close to the 5' end of the genome (nt 13 to nt 530). This region presents a low sequence similarity between the different rabbit lagovirus genotypes (RHDV, RHDV2, RCV-A1 and RCV-like non-pathogenic lagoviruses) which enables them to be distinguished. Even if the sequenced region does not cover the main recombination breakpoints identified in a region close to the cleavage site of the capsid gene [21, 37], this PCR is sufficient to detect recombinant RHDV2 strains in hares. Indeed, since there is no evidence of European hare sensibility to RHDV G1 or RCV-A1-like viruses, positive results cannot be interpreted as the detection of a co-infection. Amplification was carried out using AmpliTaq Gold DNA polymerase (Applera Applied Biosystems) (10 min at 94 °C, followed by 30 cycles at 94 °C for 30 s, 54 °C for 30 s and 72 °C for 45 s, then 10 min at 72 °C) and the primers 1U (5'-GATTAG GCCGTGAAARTTATG, 1–12) and 1L (5'-CAACGT CAACAAACTTGTCC, 550–531). The amplified DNA were purified and sequenced as mentioned above using the PCR primers. The sequences were compiled using Vector NTI Advance 11.5 (Invitrogen) and genotyping as previously described.

## Relative proportion of RHDV2 and EBHSV in lagovirus infections in 2015

### European hare samples

Among the European hares collected throughout France in 2015 by the SAGIR network displaying lesions
compatible with EBHS at necropsy, we sampled the 208 hares for which a virological diagnostic has been performed. Ten of those 208 samples also belong to the 50 samples that have been previously analysed for detection and molecular characterization of RHDV2 in hares (Table 1). The livers were stored at -20 °C before being sent to the Anses' laboratory for further molecular analyses.

### EBHSV and RHDV2 screening by real-time RT-PCRs

Each frozen liver sample was thawed and total RNA was extracted from 100 µL of the obtained exudate using the NucleoSpin® RNA (individual column) or NucleoSpin® 8 RNA (8-well strip) kits (Macherey-Nagel) according to the manufacturer's instructions.

The rRT-PCR EBHSV developed to amplify the target gene and the internal control (Intype IC-RNA, QIAGEN) in a single tube was performed using the thermal cycler "7500 Real Time PRC System" (Applied Biosystems). Amplification was carried out using the QuantiTect Probe RT-PCR Kit (QIAGEN) according to the manufacturer's instructions. Each reaction contained 5 µL of extracted RNA template and a final concentration of 100 nM of probe EBHS-Pro2 (6-FAM 5'-CTGACGCCCCTGGCACCGCTA BHQ-1®, 5302-5322) and IC-Probe (Cy 5® 5'-AGCACCCAGTCCGC CCTGAGCA BHQ-2®), of 400 nM of the target gene primers EBHS-6Fwd (5'-AATGTTATGGAGGGTAA GCC, 5277-5296) and EBHS + 115Rev (5'-GATGC TATAACGTTGTCAGC, 5416-5397), and of 400 nM of IC primers IC-1F (5'-GACCACTACCAGCAGAACAC) and IC-1R (5'-GAACTCCAGCAGGACCATG). Thermal cycling conditions used were 30 min at 48 °C, 10 min at 95 °C, followed by 40 cycles of amplification (95 °C for 15 s, 60 °C for 1 min). PCR results were analyzed using the 7500 Software (Applied Biosystems).

The rRT-PCR RHDV2 was performed on the same machine as described above by using the Power SYBR® Green RNA-to-CT™ 1-step kit (Applied Biosystems) according to the manufacturer's instructions. Each reaction contained 5 µL of extracted RNA template and a final concentration of 200 nM of primers RHDV2S-F (5'-CGACAACAGATGGAATGGTG, 6063-6082) and RHDV2S-R (5'-AACCATCTGGAG CAATTTGG, 6310-6291). Thermal cycling conditions used were 30 min at 48 °C, 10 min at 95 °C, followed by 40 cycles of amplification (95 °C for 15 s, 60 °C for 1 min) finished with a default melting curve. Analyses of cycle threshold (Ct) and melting temperature (Tm) were calculated using 7500 Software (Applied Biosystems).

Doubtful results (Ct values close to the threshold) or results suggesting an EBHSV-RHDV2 co-infection were

analysed again by RT-PCR-sequencing in order to obtain genotyping data.

## Results

### RHDV2 isolated in hares collected between 2013 and early 2015

The genotyping results obtained by RT-PCR-sequencing U38-EBHS9 performed on 50 hare samples collected between January 2013 and February 2015 revealed the presence of RHDV2 in seven samples (E13-37, E14-03, E14-40, E14-45, E14-73, E14-115, and E15-12) (Table 1). The oldest sample displaying RHDV2 infection has been collected in November 2013. These seven samples have been collected in different hare populations throughout France (Figure 1). In the E14-40 sample, our analysis showed a co-infection by both RHDV2 and EBHSV.

### Genetic and spatial relationships between RHDV2 strains characterized in hares and rabbits showed a strong epidemiological link

We obtained the complete VP60 sequences of the seven RHDV2 strains isolated in hares. Phylogenetic analyses showed that these sequences clustered into the RHDV2 phylogenetic group and thus confirmed that the viruses isolated in hares are RHDV2 (Figure 2). These RHDV2 sequences are closely related to each other but distinct. They significantly clustered into two groups, the first one (E14-40, E14-03 and E13-37) showing an average nucleotide identity of 96.3% with the second one (E14-45, E14-115, E15-12 and E14-73). The two groups contained both hare and rabbit viruses, including several RHDV2 characterized in rabbits during the same period of time and closely related to RHDV2 characterized in hares (99-99.2% nucleotide identity). They also revealed a spatial structure, with a first group that gathered viruses circulating in north-western France and a second group that gathered RHDV2 from two nearby departments in southwestern France (Figure 1).

### Search for recombinant strains

We identified recombinant strains in the south-western group only (Figure 1). The 1U/1L sequences of five RHDV2 strains, three isolated in hares (E14-45, E14-115 and E15-12) and two isolated in rabbits (13-165 and 14-59) showed ~ 98% nucleotide identity with Iberian RHDV-G1/RHDV2 recombinant strain sequences, suggesting a similar genetic origin.

### The molecular epidemiological survey revealed a high proportion and a large geographical distribution of lagovirus disease due to RHDV2 in 2015

Prevalence of EBHSV and RHDV2 in hares that died of an EBHS-like disease in 2015 was estimated by analyzing
## Hare samples

1: E14-40 (27/10/2014)

2: E14-03 (06/12/2013)

3: E13-37 (28/11/2013)

4: E14-45* (28/11/2014)

5: E14-115* (25/07/2014)

6: E15-12* (07/01/2015)

7: E14-73† (09/11/2014)

## Rabbit samples

A: 13-69 (07/07/2013)

B: 13-100 (07/10/2013)

C: 13-71 (19/06/2013)

D: 14-02 (08/12/2013)

E: 14-114 (August 2014)

F: 13-165* (15/09/2013)

G: 14-59* (25/11/2014)

H: 13-122† (October 2013)

I: 13-22 (31/01/2013)

<figure><img src="image_4.png" /><figcaption>**Figure 1** Spatial distribution of analyzed RHDV2 outbreaks in French hare and rabbit populations (January 2013–February 2015). Data are pooled by department (administrative area) and the color indicates the species involved in the analyzed outbreaks. Black: *Lepus europaeus*; light grey: *Oryctolagus cuniculus*; dark grey: the two leporid species. White: no data. The numbers refer to RHDV2 strains characterized in hares and the letters to RHDV2 strains characterized in rabbits. The order of the numbers and the letters corresponds to the order in which the VP60 sequences are positioned in the phylogenetic tree from top to bottom (Figure 2). *Recombinant RHDV-G1/RHDV2 strains, †unsuccessful 1U/1L PCR amplification.</figcaption></figure>

208 hare samples both by rRT-PCR EBHSV and rRT-PCR RHDV2. We showed that 42% (87/208) were infected by a lagovirus. Among them, 62% (54/87) were infected by EBHSV, 37% (32/87) by RHDV2, and 1% (1/87) co-infected by both viruses (Table 1). This survey also showed that EBHSV and RHDV2 were present in hare populations distributed throughout France (Figure 3). In addition, SAGIR records showed that at least three dead hares (E14-03, E14-45 and E15-12) were collected during severe outbreaks with several dead hares found in a restricted area. However, as this information was not systematically recorded, the outbreak severity of the four other cases remains uncertain.

## Discussion

Our study confirms the sensitivity of the European hare to RHDV2 and shows that closely related strains circulate locally in both rabbit and hare populations. It also demonstrates that RHDV2 is responsible for a large proportion of lagovirus disease cases in hare populations in France. These findings have several implications in terms of lagoviruses evolution and diagnosis strategy.

Our study confirms (i) the susceptibility of European hares to RHDV2 and (ii) the fact that RHDV2 and EBHSV cause similar diseases in this species, as

previously observed [24, 26]. In addition, our data show that RHDV2 may cause severe outbreaks in European hares, which had not been previously reported [24, 26]. RHDV2 is therefore highly virulent in at least three lagomorph species, the European rabbit (*O. cuniculus*), the European hare (*L. europaeus*) and the Cape hare (*L. capensis mediterraneus*), and is probably less virulent in the Italian hare (*L. corsicanus*) [7, 23, 25]. Even though one case of RHDV2 infection has recently been reported in the Mountain hare (*L. timidus*), its virulence in this species remains undocumented [27].

Our results also provide new insights into RHDV2 epidemiology in hares since they demonstrate that RHDV2 strains infecting hares do not belong to a lineage that has evolved only in hares from a single species barrier crossing. Indeed, phylogenetic analyses show that RHDV2 strains characterized in hares do not cluster into specific branches but belong to several lineages including also strains isolated in rabbits: their genetic diversity shows that RHDV2 outbreaks in hares correspond to independent infections. On the other hand, their phylogenetic position within RHDV2 strains isolated in rabbits at the same period and in the same geographical areas suggests a strong epidemiological link between cases in rabbits and hares. This finding is strengthened by the
<figure><img src="image_5.png" /><figcaption>**Figure 2** Phylogenetic tree derived for VP60 nucleotide sequences of 381 rabbit lagoviruses including 15 French RHDV2. The tree was obtained using the Neighbor-Joining method and was drawn to a scale of nucleotide substitutions per site. The percentages greater than 70% of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) are given in italics before each major branch node. The European brown hare syndrome virus (EBHSV) strain GD (Z69620) was used as an outgroup to root the tree. The RHDV2 strains characterized in hares are underlined. The 15 French RHDV2 strains characterized in this study (in rabbits and hares) are in bold and their names are followed by the number (for hares) or the letter (for rabbits) mentioned in the Figure 1 for each strain. The RHDV, the non-pathogenic and weakly pathogenic lagoviruses, most of the Iberian and Azorean RHDV2 and the Australian non-pathogenic viruses were collapsed to lighten the figure. *Recombinant RHDV-G1/RHDV2 strains, †unsuccessful 1U/1L PCR amplification.*</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 3** Spatial distributions of EBHSV and RHDV2 outbreaks in hare populations in France in 2015. Data are pooled by department (administrative area) and the color indicates the lagoviruses involved in the outbreaks. Department in white: negative samples for lagovirus or absence of sampling.</figcaption></figure>

characteristics of the strains isolated in south-western France in both species. Indeed, in this area, recombinant RHDV-G1/RHDV2 strains in rabbits and hares could be isolated. Altogether, these two features support the hypothesis of a local sharing of strains between the two species.

Another major finding is that hare infection by RHDV2 is common. Indeed, RHDV2 was involved in nearly 40% of cases of lagovirus disease in hares recorded in France in 2015 and was spatially widespread throughout France. This high proportion of RHDV2 in lagovirus disease cases may be due to intra-specific transmission among hares and also to inter-specific transmission from rabbits to hares, made possible by both the wide distribution of rabbits in France and the large spatial distribution of RHDV2 among rabbit populations.

The first case described in our study occurred in November 2013, which is later than the first case recorded in Italy in November 2012 [26]. However, since crossing the species barrier appears to be frequent, one cannot exclude that some hare infection by RHDV2 in France occurred earlier than 2013, shortly after its emergence in rabbits in 2010 [6]. Such a short time laps between RHDV2 detection in rabbits and hares was observed in Sardinia where RHDV2 was detected in Cape hares (*L. capensis mediterraneus*) only 2 months after the first case recorded in rabbits [25]. A large-scale survey on older samples is currently in progress to assess more precisely the date of the first crossing of species barrier in France.

Recombination is a common feature of rabbit lagoviruses, increasing their genetic diversity and driving their evolution [21, 37–39]. The sharing of the same virus by rabbits and hares may have major evolutionary consequences since it increases the number of viruses possibly involved and the number of hosts in which these recombination events may occur. This hypothesis is supported by the detection of two hares co-infected by EBHSV and RHDV2, which makes possible the recombination between these viruses previously considered to be specific to hares and rabbits, respectively.

The global host population size, but also its spatial structure, strongly influence the evolution of pathogens virulence and their demographic impact [29, 30]. The jump of RHDV2 to the European hare increases the connectivity between host populations because rabbit populations are patchily distributed while hares display a more continuous spatial distribution. Thus, the enlargement of host range may have contrasted effects on RHDV2 virulence and impact on rabbit and hare populations. First, a less fragmented host population could enable a steady circulation of the virus in both hosts and maintain higher immunity levels [40]. By increasing the herd immunity, the change of host range may therefore lead to a decrease of RHDV2 impact in rabbits. Conversely, the change in the spatial structure of the global host population may contribute to increasing the mean virulence of RHDV2. Indeed, Fouchet et al. [41], by modelling, suggest that “the mean virulence of RHDV increases with connectivity”. One may therefore expect a higher virulence of RHDV2, which is supported by recent data [42].

This also raises the question of the competition between three viruses (RHDV, RHDV2 and EBHSV) two of which are genus specific and one is able to infect both rabbit and hare species. Yet, RHDV2 is replacing RHDV in rabbit populations in south-western Europe [6, 9, 13, 19, 20]. Thus one may wonder about a possible future replacement of EBHSV by RHDV2 in hare populations. In this context, the role, evolution and spreading of recombinant strains, such as RHDV-G1/RHDV2, may also be of particular interest to assess whether recombinant strains benefit from a selective advantage.

Another consequence of RHDV2 infection in hares concerns the diagnostic strategy which has to be redefined in order to look for both EBHSV and RHDV2 in case of EBHSV suspicion. A new rRT-PCR has been recently developed to amplify EBHSV and RHDV2 genes along with the exogenous internal control in a single tube (Lemaitre et al., unpublished data). This new method will allow the monitoring of the circulation of the two viruses in hares during epidemiological surveys and to detect the possible replacement of EBHSV by RHDV2. Moreover, in
the case of surveys aiming at identifying the viral strains involved in lagovirus disease outbreaks, the detection of recombinant strains would also be useful for a better understanding of RHDV2 evolution in this species. Regardless of lagovirus strains identification, our data show that a high proportion (121/208 in 2015) of hares displaying lesions compatible with EBHS have not died of a lagovirus disease. Other causes may induce similar or resembling lesions and a better diagnosis of those cases will be another important issue for the future.

Lastly, RHDV2 ability to infect a wide range of species hosts and to cause in European hare a disease similar to EBHS should lead one to reconsider both viruses and diseases taxonomy. Until now, EBHSV and RHDV have been considered as two different virus species mainly on the basis of their distinct host range, hares and rabbits, respectively. The jump of RHDV2 to different hare species makes this classification obsolete. A new nomenclature based on genetic criteria and therefore more powerful has recently been proposed [43]. Similarly, the names of the diseases, EBHS and RHD, are related to the affected host. The passage of RHDV2 from the rabbit to the hare, in which it causes a disease in all respects similar to EBHS, shows the limits of this nomenclature. We are now dealing with a disease—EBHS—whose name calls back to a virus that is not always responsible for this disease. We therefore propose to rename this disease hare haemorrhagic disease (HHD) or hare lagovirus disease (HLV), independently of the viral strains that cause it. This nomenclature is similar to that adopted in the rabbit where RHD designates the disease, irrespective of the virus responsible for it, RHDV or RHDV2.

## Additional files

**Additional file 1. Lagovirus positive samples used for the detection and molecular characterization of RHDV2 in hares.** The table gives the sample number, the year of collection, the French department number, the genotype and the EMBL/Genbank accession number in nucleotide databases of the RHDV2 VP60 sequence of the lagoviruses characterized in hares (n = 7) and of nine RHDV2 characterized in dead rabbits during the same period of time. \*Recombinant RHDV-G1/RHDV2 strains.

**Additional file 2. Lagovirus positive samples used for estimating the proportion of RHDV2 in lagovirus infection in hares.** The table gives the sample number, the year of collection, the French department number and the genotype of the lagoviruses characterized during the screening performed to estimate the relative proportion of RHDV2 and EBHSV in 208 hares displaying lesions compatible with EBHSV at necropsy.

## Competing interests

The authors declare that they have no competing interests.

## Authors' contributions

Study development and design: GLG, SM, SB, JSG. Molecular genetic studies and phylogenetic analyses: EL, CH, ST, GLG, SB. Hare sampling: AD, SM, JSG.

Drafting of the manuscript: GLG, SM, JSG. Writing of the manuscript: GLG, SM, JSG, EL, SB. All authors read and approved the final manuscript.

## Acknowledgements

We thank the departmental public veterinary laboratories which performed EBHS diagnostic and sent us liver samples to characterize the viruses responsible for the observed mortalities: "Inovalyss-Angers" (Angers, France) and "Laboratoire de l'Environnement et de l'Alimentation de la Vendée" (La Roche-sur-Yon, France). We thank the departmental services of the ONCFS, the French departmental, regional and national hunter federations, the departmental public veterinary laboratories, the Anses Laboratory for Rabies and Wildlife of Nancy (LRFSN) and all the people who reported wildlife mortality events. We also thank Julie Tucolet (ONCFS) for her assistance with SAGIR data. We are grateful to Jérémy Dufay and Marina Chatel who contributed to developing the real-time RT-PCR EBHSV and to screening the lagoviruses for the molecular epidemiological survey, respectively.

## Funding

This work was performed within the framework of the ECALEP project supported by the ERA-Net Anihwa (Animal Health and Welfare), a coordination action funded under the European Commission's ERA-Net scheme within the Seventh Framework Programme (Contract No. 291815). The ECALEP project is funded by the ANR (France, contracts ANR-14-ANWA-0004-01, ANR-14-ANWA-0004-02, ANR-14-ANWA-0004-03, ANR-14-ANWA-0004-04), the Ministry of Health, Dep. for Veterinary Public Health, Nutrition & Food Safety (Italy) and the Research council FORMAS (Sweden, contract FORMAS 221-2014-1841). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

## Author details

¹ Avian and Rabbit Virology, Immunology and Parasitology Unit, Anses, French Agency for Food, Environmental and Occupational Health & Safety, Ploufragan-Plouzané Laboratory, BP 53, 22440 Ploufragan, France. ² University Bretagne Loire, CS 54417, 35044 Rennes Cedex, France. ³ IHAP, University of Toulouse, INRA, ENVT, 31076 Toulouse, France. ⁴ Research Department, ONCFS, National Hunting and Wildlife Agency, USF, BP 20, 78610 Le Perray-en-Yvelines, France. ⁵ Research Department, ONCFS, National Hunting and Wildlife Agency, UFP, CS 42355, 44323 Nantes Cedex 3, France.

## Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Received: 10 April 2017 Accepted: 11 October 2017

Published online: 28 October 2017

## References

1. Abrantes J, van der Loo W, Le Pendu J, Esteves P (2012) Rabbit haemorrhagic disease (RHD) and rabbit haemorrhagic disease virus (RHDV): a review. Vet Res 43:12

2. Marcato PS, Benazzi C, Vecchi G, Galeotti M, Della Salda L, Sarli G, Lucidi P (1991) Clinical and pathological features of viral haemorrhagic disease of rabbits and the European brown hare syndrom. Rev Sci Tech 10:371-392

3. Le Gall-Reculé G, Zwingelstein F, Laurent S, De Boisseson C, Portejoie Y, Rasschaert D (2003) Phylogenetic analysis of Rabbit haemorrhagic disease virus in France between 1993 and 2000, and the characterisation of RHDV antigenic variants. Arch Virol 148:65-81

4. Kerr PJ, Kitchen A, Holmes EC (2009) The origin and phylodynamics of rabbit haemorrhagic disease virus. J Virol 83:12129-12138

5. Kinnear M, Linde CC (2010) Capsid gene divergence in rabbit hemorrhagic disease virus. J Gen Virol 91:174-181

6. Le Gall-Recule G, Lavazza A, Marchandeau S, Bertagnoli S, Zwingelstein F, Cavadini P, Martinelli N, Lombardi G, Guerin J-L, Lemaitre E, Decors A, Boucher S, Le Normand B, Capucci L (2013) Emergence of a new lagovirus related to rabbit haemorrhagic disease virus. Vet Res 44:81

7. Le Gall-Reculé G, Zwingelstein F, Boucher S, Le Normand B, Plassiart G, Portejoie Y, Decors A, Bertagnoli S, Guérin JL, Marchandeau S (2011)
Detection of a new variant of rabbit haemorrhagic disease virus in France. Vet Rec 168:137–138

8. Abrantes J, Lopes AM, Dalton KP, Melo P, Correia JJ, Ramada M, Alves PC, Parra F, Esteves PJ (2013) New variant of rabbit hemorrhagic disease virus, Portugal, 2012–2013. Emer Infect Dis 19:1900–1902

9. Dalton KP, Nicieza I, Abrantes J, Esteves PJ, Parra F (2014) Spread of new variant RHDV in domestic rabbits on the Iberian Peninsula. Vet Microbiol 169:67–73

10. Westcott DG, Choudhury B (2015) Rabbit haemorrhagic disease virus 2-like variant in Great Britain. Vet Rec 176:74

11. Neimanis AS, Ahola H, Zohari S, Larsson Pettersson U, Bröjer C, Capucci L, Gavier-Widén D. Arrival of rabbit haemorrhagic disease virus 2 to northern Europe: emergence and outbreaks in wild and domestic rabbits (*Oryctolagus cuniculus*) in Sweden. Transbound Emerg Dis, in press. doi:10.1111/tbed.12650

12. Martin-Alonso A, Martin-Carrillo N, Garcia-Livia K, Valladares B, Foronda P (2016) Emerging rabbit haemorrhagic disease virus 2 (RHDV2) at the gates of the African continent. Infect Genet Evol 44:46–50

13. Duarte M, Carvalho C, Bernardo S, Barros SV, Benevides S, Flor L, Monteiro M, Marques I, Henriques M, Barros SC, Fagulha T, Ramos F, Luís T, Fevereiro M (2015) Rabbit haemorrhagic disease virus 2 (RHDV2) outbreak in Azores: disclosure of common genetic markers and phylogenetic segregation within the European strains. Infect Genet Evol 35:163–171

14. Hall RN, Mahar JE, Haboury S, Stevens V, Holmes EC, Strive T (2015) Emerging rabbit hemorrhagic disease virus 2 (RHDVb), Australia. Emerg Infect Dis 21:2276–2278

15. World Animal Health Information Database (WAHIS interface). http://www.oie.int/wahis_2/public/wahid.php/Reviewreport/Review?page_ref er=MapEventSummary&reportid=18202. Accessed 22 Jul 2015

16. World Animal Health Information Database (WAHIS interface). http://www.oie.int/wahis_2/public/wahid.php/Reviewreport/Review?page_ref er=MapFullEventReport&reportid=20799. Accessed 25 Aug 2016

17. World Animal Health Information Database (WAHIS interface). http://www.oie.int/wahis_2/public/wahid.php/Reviewreport/Review?page_ref er=MapEventSummary&reportid=19903. Accessed 15 Mar 2016

18. World Animal Health Information Database (WAHIS interface). http://www.oie.int/wahis_2/public/wahid.php/Reviewreport/Review?page_ref er=MapFullEventReport&reportid=20852. Accessed 30 Aug 2016

19. Lopes A, Correia J, Abrantes J, Melo P, Ramada M, Magalhães M, Alves P, Esteves P (2015) Is the new variant RHDV replacing genogroup 1 in Portuguese wild rabbit populations? Viruses 7:27–36

20. Calvete C, Sarto P, Calvo AJ, Monroy F, Calvo JH (2014) Could the new rabbit haemorrhagic disease virus variant (RHDVb) be fully replacing classical RHD strains in the Iberian Peninsula? World Rabbit Sci 22:91

21. Lopes AM, Dalton KP, Magalhães MJ, Parra F, Esteves PJ, Holmes EC, Abrantes J (2015) Full genomic analysis of new variant rabbit hemorrhagic disease virus revealed multiple recombination events. J Gen Virol 96:1309–1319

22. Almeida T, Lopes AM, Magalhães MJ, Neves F, Pinheiro A, Gonçalves D, Leitão M, Esteves PJ, Abrantes J (2015) Tracking the evolution of the G1/RHDVb recombinant strains introduced from the Iberian Peninsula to the Azores islands, Portugal. Infect Genet Evol 34:307–313

23. Camarda A, Pugliese N, Cavadini P, Circella E, Capucci L, Caroli A, Legretto M, Mallia E, Lavazza A (2014) Detection of the new emerging rabbit haemorrhagic disease type 2 virus (RHDV2) in Sicily from rabbit (*Oryctolagus cuniculus*) and Italian hare (*Lepus corsicanus*). Res Vet Sci 97:642–645

24. Hall RN, Peacock DE, Kovaliski J, Mahar JE, Mourant R, Piper M, Strive T (2017) Detection of RHDV2 in European brown hares (*Lepus europaeus*) in Australia. Vet Rec 180:121

25. Puggioni G, Cavadini P, Maestrale C, Scivoli R, Botti G, Ligios C, Le Gall-Reculé G, Lavazza A, Capucci L (2013) The new French 2010 rabbit hemorrhagic disease virus causes an RHD-like disease in the Sardinian Cape hare (*Lepus capensis mediterraneus*). Vet Rec 44:96

26. Velarde R, Cavadini P, Neimanis A, Cabezón O, Chiari M, Gaffuri A, Lavín S, Grilli G, Gavier-Widén D, Lavazza A, Capucci L. Spillover events of infection of brown hares (*Lepus europaeus*) with rabbit haemorrhagic disease type 2 virus (RHDV2) caused sporadic cases of an European brown hare syndrome-like disease in Italy and Spain. Transbound Emerg Dis, in press. doi:10.1111/tbed.12562

27. World Animal Health Information Database (WAHIS interface). http://www.oie.int/wahis_2/public/wahid.php/Reviewreport/Review?page_ref er=MapFullEventReport&reportid=22452. Accessed 20 Jan 2017

28. Gavier-Widen D, Mörner T (1991) Epidemiology and diagnosis of the European brown hare syndrome in Scandinavian countries: a review. Rev Sci Tech 10:453–458

29. Read JM, Keeling MJ (2003) Disease evolution on networks: the role of contact structure. Proc Biol Sci 270:699–708

30. Boots M, Hudson PJ, Sasaki A (2004) Large shifts in pathogen virulence relate to host population structure. Science 303:842–844

31. Rogers PM (1981) Ecology of the European wild rabbit *Oryctolagus cuniculus* (L.) in Mediterranean habitats. II. Distribution in the landscape of the Camargue, S. France. J Appl Ecol 18:355–371

32. Rogers PM, Myers K (1979) Ecology of the European wild rabbit, *Oryctolagus cuniculus* (L.), in Mediterranean habitats. I. Distribution in the landscape of the Coto Donana, S. Spain. J Appl Ecol 16:691–703

33. Lewandowski K, Nowakowski JJ (1993) Studies on the European hare. 49. Spatial distribution of brown hare *Lepus europaeus* populations in habitats of various types of agriculture. Acta Theriol 38:435–442

34. Millot F, Decors A, Mastain O, Quintainte T, Berny P, Vey D, Lasseur R, Bro E (2017) Field evidence of bird poisonings by imidacloprid-treated seeds: a review of incidents reported by the French SAGIR network from 1995 to 2014. Environ Sci Pollut Res Int 24:5469–5485

35. Lopes AM, Gavier-Widén D, Le Gall-Reculé G, Esteves PJ, Abrantes J (2013) Complete coding sequences of European brown hare syndrome virus (EBHSV) strains isolated in 1982 in Sweden. Arch Virol 158:2193–2196

36. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S (2011) MEGAS: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol 28:2731–2739

37. Forrester NL, Moss SR, Turner SL, Schirrmeier H, Gould EA (2008) Recombination in rabbit haemorrhagic disease virus: possible impact on evolution and epidemiology. Virology 376:390–396

38. Abrantes J, Esteves PJ, Van der Loo W (2008) Evidence for recombination in the major capsid gene VP60 of the rabbit haemorrhagic disease virus (RHDV). Arch Virol 153:329–335

39. Hu B, Wang F, Fan Z, Song Y, Abrantes J, Zuo Y, Esteves PJ (2017) Recombination between G2 and G6 strains of rabbit hemorrhagic disease virus (RHDV) in China. Arch Virol 162:269–272

40. Fouchet D, Marchandeau S, Bahi-Jaber N, Pontier D (2007) The role of maternal antibodies in the emergence of severe disease as a result of fragmentation. J R Soc Interface 4:479–489

41. Fouchet D, Le Pendu J, Guitton JS, Guiserix M, Marchandeau S, Pontier D (2009) Evolution of microparasites in spatially and genetically structured host populations: the example of RHDV infecting rabbits. J Theor Biol 257:212–227

42. Capucci L, Cavadini P, Schiavitto M, Lombardi G, Lavazza A (2017) Increased pathogenicity in rabbit haemorrhagic disease virus type 2 (RHDV2). Vet Rec 180:426

43. Le Pendu J, Abrantes J, Bertagnoli S, Guitton J-S, Le Gall-Reculé G, Lopes AM, Marchandeau S, Alda F, Almeida T, Célio AP, Bárcena J, Burmakina G, Blanco E, Calvete C, Cavadini P, Cooke B, Dalton K, Delibes Mateos M, Deptula W, Eden JS, Wang F, Ferreira CC, Ferreira P, Foronda P, Gonçalves D, Gavier-Widén D, Hall R, Hukowska-Szematowicz B, Kerr P, Kovaliski J (2017) Proposal for a unified classification system and nomenclature of lagoviruses. J Gen Virol 98:1658–1666