# Is the New Variant RHDV Replacing Genogroup 1 in Portuguese Wild Rabbit Populations?

Ana M. Lopes<sup>1,2,3,†</sup>, Jorge Correia<sup>4,†</sup>, Joana Abrantes<sup>1</sup>, Pedro Melo<sup>5</sup>, Margarida Ramada<sup>5</sup>, Maria J. Magalhães<sup>1</sup>, Paulo C. Alves<sup>1,2,6</sup> and Pedro J. Esteves<sup>1,2,7,*</sup>

**Abstract:** The *Lagovirus* rabbit hemorrhagic disease virus (RHDV), a member of the family *Caliciviridae*, severely affects European rabbit (*Oryctolagus cuniculus*) populations by causing rabbit hemorrhagic disease (RHD). RHDV is subdivided in six genogroups but, more recently, a new RHDV variant with a unique genetic and antigenic profile emerged. We performed a study in rabbits found dead in the field during 2013 and 2014 in Portugal to determine the prevalence of this new variant *versus* the classical RHDV. Fifty-seven liver

<footnote n="1">¹ CIBIO, InBIO—Research Network in Biodiversity and Evolutionary Biology, Universidade do Porto, Campus de Vairão, Rua Padre Armando Quintas, 4485-661 Vairão, Portugal; E-Mails: analopes@cibio.up.pt (A.M.L.); jabrantes@cibio.up.pt (J.A.); magalhaes.mjtr@hotmail.com (M.J.M.); pcalves@fc.up.pt (P.C.A.)</footnote>

<footnote n="2">² Departamento de Biologia, Faculdade de Ciências da Universidade do Porto, Rua do Campo Alegre, s/n, 4169-007, Porto, Portugal</footnote>

<footnote n="3">³ INSERM—Institut national de la santé et de la recherche médicale, U892, Université de Nantes, 44007 Nantes, France</footnote>

<footnote n="4">⁴ CIISA—Centro de Investigação Interdisciplinar em Sanidade Animal, Faculdade de Medicina Veterinária, Universidade de Lisboa, 1300-477 Lisboa, Portugal; E-Mail: jcorreia@fmv.ulisboa.pt</footnote>

<footnote n="5">⁵ VETNATURA, 1600-176 Lisboa, Portugal; E-Mails: vetnatura@gmail.com (P.M.); margarida.ramada@dcv.eu.com (M.R.)</footnote>

<footnote n="6">⁶ Wildlife Biology Program, University of Montana, 32 Campus Drive, Missoula, MT 59812, USA</footnote>

<footnote n="7">⁷ CITS—Centro de Investigação em Tecnologias da Saúde, IPSN, CESPU, 4585-116 Gandra, Portugal</footnote>

<footnote>† These authors contributed equally to this work.</footnote>

<footnote>* Author to whom correspondence should be addressed; E-Mail: pjesteves@cibio.up.pt; Tel.: +351-252-660-411.</footnote>
samples were screened for the presence of RHDV and positive samples were genotyped.
All cases of RHDV infection were caused by the new variant. The only former genogroup
circulating in Portugal, G1, was not detected. We hence conclude that the new RHDV variant
is replacing G1 in Portugal, probably due to a selective advantage. This sudden and rapid
replacement emphasizes the necessity of continued monitoring of wild rabbit populations.

**Keywords:** rabbit hemorrhagic disease virus (RHDV); rabbit hemorrhagic disease (RHD); new variant RHDV; genogroup 1; replacement; Portugal

**1. Introduction**

Rabbit hemorrhagic disease virus (RHDV), a member of the genus *Lagovirus*, family *Caliciviridae*, is responsible for causing rabbit hemorrhagic disease (RHD). RHD was first reported in China in 1984 [1] and is a highly contagious and fatal disease for European rabbit (*Oryctolagus cuniculus*). The acute form of the disease can cause up to 80%–100% of mortality within 1–3 days [2,3]. RHD is associated with hepatic necrosis, hemorrhages and congestion in several organs, including the upper respiratory tract and lungs, and splenomegaly [2]. The disease became endemic in several European countries, where it was first reported in 1986 [4], and in Australia and New Zealand, where it was introduced as biological control agent [5]. In Portugal, RHD was described in 1989 (reviewed in [3]).

RHDV is a non-enveloped, positive-sense, single-stranded RNA virus, with a genome of ~7.4 kb organized in two ORFs, with ORF1 encoding the major structural capsid protein gene, VP60. Historically, six RHDV genogroups were recognized based on the VP60: genogroups 1–5 (G1–G5) and the antigenic variant RHDVa or G6 [6,7]. However, more recent studies showed that G3, G4 and G5 belong to the same genetic cluster [8,9]. Furthermore, weak pathogenic and non-pathogenic forms of the virus exist, which are genetically different from the pathogenic RHDV, exhibiting ~20% divergence [10–14]. The other member of the genus *Lagovirus*, the closely related European brown hare syndrome virus (EBHSV), has a similar organization of the genome, albeit presenting ~30% divergence from RHDV [15].

In 2010, a new RHDV variant emerged in France [16], tentatively named RHDV2 or RHDVb, that rapidly spread to several European countries including Portugal, where it was described in 2012 [17–21]. With a unique antigenic profile, it causes mortality in young rabbits (less than two months) and vaccinated rabbits, that are typically not susceptible to G1–G6 strains [17]. Vaccinated rabbits appear to be less susceptible to this new variant, albeit protection conferred by the vaccine is not complete [19]. Before the emergence of the new variant, only G1 strains were known to circulate in Portugal [22–24]. RHDVa was also described in the Iberian Peninsula [25], but this description refers only to domestic rabbits and no reported cases exist in wild animals.

The Iberian Peninsula is inhabited by two subspecies of *O. cuniculus*: one occupying the northeastern part (*O. c. cuniculus*) and the other restricted to the southwest and the islands of Madeira, Azores and Canaries (*O. c. algirus*) [reviewed in 26]. Both subspecies are equally affected by RHDV [22], but the susceptibility of the subspecies *O. c. algirus* acquires special importance since its distribution is
restricted to the Iberian Peninsula. Thus, this subspecies is more vulnerable to severe outbreaks, which
might have implications on its conservation status.

In this work, we performed a genetic study to determine the prevalence of the new variant comparing
to G1 in Portuguese rabbit populations and a histopathological study to characterize the lesions found in
young rabbits. Our results confirm that the new variant has replaced G1 in Portugal in the last two years,
possibly reflecting a selective advantage.

**2. Materials and Methods**

Fifty-seven liver samples were collected in 12 localities covering from north to south of Portugal, between 2013 and 2014 (Figure 1), and frozen at -20 °C. All rabbits were found dead in the field with lesions consistent with RHD and, based on the geographical distribution, belong to the subspecies *O. c. algirus* [26]. In addition, for seven young wild rabbits, carcasses were submitted to necropsy and samples of lung, liver, thymus, trachea, kidney and spleen were collected for histopathological exam. Tissue samples were fixed in 10% buffered formalin and embedded in paraffin. Sections of 3 µm thick were stained with haematoxylin and eosin for routine microscopical examination.

Total RNA was extracted from the liver of all samples with the RNeasy Mini Kit (Qiagen, Hilden,
Germany). Up to 30 mg of liver were collected and tissue was homogenized in a rotor-stator
homogenizer (Mixer Mill MM400 from Retsch, Haan, Germany) at 30 Hz for 7 min. RNA extraction
followed the manufacturer's instructions. Viral cDNA was synthesized using oligo(dT) primers and
SuperScriptTM III Reverse Transcriptase (Invitrogen, Carlsbad, CA, USA) following the manufacturer's
instructions. A screening to test for the presence of RHDV was performed with two pairs of primers: the
first pair is specific for the new variant: RHNavF 5'-GTTCGTCAAATGTACTTGAGC-3' and
RHNavR 5'-GTGTACGTAATGGCACTACTG-3' [27]. Samples that tested negative were screened
with a second more general pair of primers (RHDV4831F 5'-GTGTATGCCATGACTCCGAT-3' and
RHDV_VP60_0467R 5'-GCGTCGATGACAACATGAG-3'). These primers were used in previous
studies and successfully amplified several RHDV strains including G1, RHDVa and G3-G5
like [24,25,28,29]. For positive samples, the complete capsid protein gene, VP60, with 1740 base pairs,
was amplified (PCR primers and conditions are available upon request). After purification, PCR
products were sequenced on an automatic sequencer ABI PRISM 310 Genetic Analyzer (PE Applied
Biosystems, Foster City, CA, USA) with the amplification primers. These sequences were deposited in
GenBank with the accession numbers KM115667-KM115716 and aligned in BioEdit version
7.0.9.0 [30]. Each sequence was compared and blasted against publicly available sequences in NBCI
(http://blast.ncbi.nlm.nih.gov/Blast.cgi). In order to facilitate visualization and the subsequent analyses,
a maximum of ten sequences of the main RHDV groups were used (G1-G6, the new variant,
non-pathogenic RCV-A1 and weakly pathogenic MRCV; EBHSV was also included). Genetic distances
between our samples and each RHDV group were calculated in MEGA6 (p-distance method,
1000 replicates) [31]. Additionally, a Maximum Likelihood (ML) tree, using a representative of each
unique haplotype found, was estimated using the same software. The best-fit nucleotide substitution
model as determined in MEGA6 (GTR+G+I) and a bootstrap resampling analysis of 1000 replicates,
were used to reconstruct the phylogeny.
<figure><img src="image_1.png" /></figure>

**Figure 1.** Map of Portugal with sampling localities indicated. Localities sampled in 2013 appear underlined; in parenthesis, the number of positive samples per total number of samples is shown. The approximate locations of the first descriptions of the new variant's presence in Portugal according to [18] are marked by the black stars.

**3. Results and Discussion**

Between 2013 and 2014, Portuguese wild rabbit populations suffered severe RHD outbreaks. Indeed, during that period, we received in our laboratory rabbit samples from several Portuguese localities (Figure 1) that were shown to be positive for the new RHDV variant. This new variant was described in Portugal in 2012 [18]. Interestingly, before that, G1 was the only genogroup known to be circulating in Portugal [22–24].

While the macroscopic and histopathological alterations caused by RHDV are well-known, little is known on how the new RHDV variant affects wild animals. As this is a new and emergent virus, monitoring these alterations might give insights on eventual changes in lesions. At necropsy, the
alterations and histopathological lesions observed were similar to those observed for other animals infected with the new RHDV variant and that are typical of a RHDV infection [16,17,19]. Necropsy of the seven young rabbits revealed epistaxis, hemorrhagic tracheitis with sero-hemorrhagic exudate in the lumen of the trachea, sero-hemorrhagic pleurisy, pulmonary congestion (predominantly in the dorsal region), and multiple hemorrhages scattered throughout the organs. Livers were predominantly pale in appearance and in some cases were congested. The thymus, spleen and kidneys showed also congestion. Histopathology revealed hemorrhagic pneumonia and tracheitis, congestion of the liver and diffuse necrotizing hepatitis (Figure 2A,B). This first establishment of a diagnosis and consequent cause of death was then confirmed by the molecular techniques.

<figure><img src="image_2.png" /><figcaption>**Figure 2.** Histopathological characterization of the new rabbit hemorrhagic disease virus (RHDV) variant in Portuguese samples. (A) Liver section of a healthy rabbit. Staining H&E, scale bar = 20 µm; (B) Liver with lesions of microvacuolar degeneration and necrotizing hepatitis with karyolysis. Staining H&E, scale bar = 20 µm.</figcaption></figure>

Of the 57 rabbit liver samples, 50 were positive for the new variant (87.7%; Figure 1); all but one were successfully amplified for the entire capsid gene. The 143 nucleotide differences defined 33 unique haplotypes, with an overall identity of 98.7% (±0.1%). None of the samples was positive for G1 or for any other genogroup. The blast analysis revealed 96%–98% of identity with publicly available sequences of the new variant from Italy and France (GenBank accession numbers KC741409, HE800529-32, FR819781, HE819400, JQ929052, KC345611-13, JX106022-23, KC907712), and an overall identity of 82% with G1 strains. The analysis of the nucleotide identities between our dataset and genogroups showed an identity of 81.6% (± 0.8%) with G2, 82.1% (± 0.8%) with RHDVa and 82.2% (± 0.8%) with G3-G5. Furthermore, our sequences presented a nucleotide identity of 81% (± 1.0%) with the weakly pathogenic MRCV, 79.2% (± 0.8%) with the non-pathogenic RCV-A1 (strain MIC-07) and 70.0% (± 0.9%) with EBHSV (overall identity with strains O4021-9, BS89 and EBHSV-GD; GenBank accession numbers of these sequences are indicated in Figure 3). The ML tree also confirmed that all samples belong to the new RHDV variant group (Figure 3), since all RHDV samples recovered from Portuguese rabbits cluster with publicly available sequences described as belonging to the new RHDV variant (bootstrap value of 99; Figure 3). Comparison with the partial sequences of the capsid obtained in the first reported cases of the new RHDV variant in Portugal (GenBank accession numbers KF442960-64) [18], revealed that our sequences present between 0.0% (± 0.2%) and 2.5% (± 0.6%)
divergence. However, the small length of the analyzed fragment does not allow us to make any conclusions about the epidemiology and dispersal of the virus since the initial outbreaks.

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

**Figure 3.** Maximum Likelihood tree of the capsid gene. Bootstrap values of the main nodes (≥0.80) are shown next to them. To facilitate visualization, the new variant group is collapsed (previously identified new variant strains have their accession numbers indicated). Samples were named according to the sampling locality, the year of collection and the number of the isolate (cf. Figure 1); when known, month and year of collection are indicated. Scale bar indicates the nucleotide substitutions per site. GenBank accession numbers of the publicly available sequences used in this study are indicated.
Indeed, the lack of epidemiological data on RHDV prevalence in Portuguese wild rabbit populations impairs a precise comparison between G1 and the new variant virulence/pathogenicity. However, the number of dead animals found in the field in 2013 and 2014 seems to have substantially increased, as has been described in Spain by Delibes-Mateos et al [32]. Our results show that these RHD outbreaks were caused by the new variant and no apparent cases of G1 (or other genogroup) were detected. This indicates that the new variant is the main variant currently circulating in Portugal and has potentially replaced G1; the factors contributing to this replacement are still unclear, but are currently under study. A similar situation has been described in France and Spain, where the new RHDV variant appears to be replacing G5 and G1, respectively [19,33,34].

After the first description of the new variant in Portugal in three Portuguese localities in 2012 and 2013 [18], the new RHDV variant was able to disperse to the center and south of Portugal (Figure 1), suggesting a high capacity of dispersal in a very short period of time. The weak dispersal to the north of Portugal is possibly due to low rabbit densities. By contrast, in the south, the high densities of rabbits facilitate virus dispersal. The high capacity of dispersal, associated with the replacement of G1, suggests a selective advantage of the new variant over other genogroups, possibly by overcoming existing immunity to older strains. Le Gall-Reculé and colleagues (2013) reported that the new RHDV variant is less virulent than classical RHDV in experimental conditions and that this could be advantageous for the new variant, by giving more time to the virus to spread. The different mortality rates observed between studies (e.g., [17,19]) make possible another scenario, where the new variant is still evolving to acquire its optimal virulence. Another possibility is that the new RHDV variant might be evolving towards higher levels of virulence, as observed in a recent study in Australia, where RHDV virulence has been increasing due to development of genetic resistance in rabbits in Australia [35].

The threat presented by the new RHDV variant to the European rabbit might be particularly high for *O. c. algirus*, since distribution of this subspecies is limited to the southwest of Iberian Peninsula. In addition, the presence and circulation of the new RHDV variant compromises the reintroduction of the most endangered felid in the world (Iberian lynx, *Lynx pardinus*) in Portugal, as rabbits represent 90% of its diet. The European rabbit is also economically important as game species, being a traditionally coveted prey, with high social value. Since the epidemiological outcome of G1 replacement by the new variant is unpredictable, and the new variant is causing important ecological and economical losses in Portugal, rabbit translocations should also be highly restricted, and continued monitoring of natural rabbit populations is strongly recommended.

## Acknowledgments

The authors are thankful to Prof. Francisco Parra and Dr. Kevin P. Dalton for sharing the sequence of the primers specific for the new variant. We are also grateful to Fençaça (Jacinto Amaro), Federação de Caçadores do Algarve (Vitor Palmilha), Clube dos Caçadores da Freguesia da Madalena (Julio Dias Pouseiro), José Manuel Maia and José Franco for providing wild rabbit samples for this study.

This work is funded by National Funds through FCT—Foundation for Science and Technology under the Project FCT-ANR/BIA-BIC/0043/2012. FCT also supported the doctoral fellowship of AML (ref.: SFRH/BD/78738/2011) and the Investigator grant of JA (ref.: IF/01396/2013). "Genomics Applied To Genetic Resources" co-financed by North Portugal Regional Operational Programme 2007/2013
(ON.2-O Novo Norte), under the National Strategic Reference Framework (NSRF), through the European Regional Development Fund (ERDF) also supported the work.

## Author Contributions

Ana M. Lopes performed the genetic analysis, analyzed the data and drafted the manuscript; Jorge Correia, Pedro Melo and Margarida Ramada carried out the necropsies and the histopathological analysis; Jorge Correia also drafted the manuscript; Maria J. Magalhães performed the genetic analysis; Paulo C. Alves participated in the sampling and helped writing the manuscript; Joana Abrantes and Pedro J. Esteves discussed the data and helped writing the manuscript.

## Conflicts of Interest

The authors declare no conflict of interest.

## References and Notes

1. Liu, S.; Xue, H.; Pu, B.; Qian, N. A new viral disease in rabbit. *Anim. Husb. Vet. Med.* **1984**, *16*, 253-255.

2. Marcato, P.S.; Benazzi, C.; Vecchi, G.; Galeotti, M.; Salda, L.D.; Sarli, G.; Lucidi, P. Clinical and pathological features of viral haemorrhagic disease of rabbits and the European brown hare syndrome. *Rev. Sci. Tech.* **1991**, *10*, 371-392.

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

4. Cancellotti, F.M.; Renzi, M. Epidemiology and current situation of viral haemorrhagic disease of rabbits and the European brown hare syndrome in Italy. *Rev. Sci. Tech. Off. Int. Epiz.* **1991**, *10*, 409-422.

5. Mutze, G.; Kovaliski, J.; Butler, K.; Capucci, L.; McPhee, S. The effect of rabbit population control programmes on the impact of rabbit haemorrhagic disease in south-eastern Australia. *J. Appl. Ecol.* **2010**, *47*, 1137-1146.

6. Le Gall, G.; Arnauld, C.; Boilletot, E.; Morisse, J.P.; Rasschaert, D. Molecular epidemiology of rabbit haemorrhagic disease virus outbreaks in France during 1988 to 1995. *J. Gen. Virol.* **1998**, *79*, 11-16.

7. Le Gall-Reculé, G.; Zwingelstein, F.; Laurent, S.; de Boisséson, C.; Portejoie, Y.; Rasschaert, D. Phylogenetic analysis of rabbit haemorrhagic disease virus in France between 1993 and 2000, and the characterisation of RHDV antigenic variants. *Arch. Virol.* **2003**, *148*, 65-81.

8. Kerr, P.J.; Kitchen, A.; Holmes, E.C. Origin and phylodynamics of rabbit hemorrhagic disease virus. *J. Virol.* **2009**, *83*, 12129-12138.

9. Kinnear, M.; Linde, C.C. Capsid gene divergence in rabbit hemorrhagic disease virus. *J. Gen. Virol.* **2010**, *91*, 174-181.

10. Capucci, L.; Fusi, P.; Lavazza, A.; Pacciarini, M.L.; Rossi, C. Detection and preliminary characterization of a new rabbit calicivirus related to rabbit hemorrhagic disease virus but nonpathogenic. *J. Virol.* **1996**, *70*, 8614-8623.
11. Strive, T.; Wright, J.D.; Robinson, A.J. Identification and partial characterisation of a new lagovirus in Australian wild rabbits. *Virology* **2009**, *384*, 97–105.

12. Bergin, I.L.; Wise, A.G.; Bolin, S.R.; Mullaney, T.P.; Kiupel, M.; Maes, R.K. Novel calicivirus identified in rabbits, Michigan, USA. *Emerg. Infect. Dis.* **2009**, *15*, 1955–1962.

13. Abrantes, J.; Esteves, P.J. Not-so-novel Michigan rabbit calicivirus. *Emerg. Infect. Dis.* **2010**, *16*, 1331–1332.

14. Le Gall-Reculé, G.; Zwingelstein, F.; Fages, M.-P.; Bertagnoli, S.; Gelfi, J.; Aubineau, J.; Roobrouck, A.; Botti, G.; Lavazza, A.; Marchandeau, S. Characterisation of a non-pathogenic and non-protective infectious rabbit lagovirus related to RHDV. *Virology* **2011**, *410*, 395–402.

15. Le Gall, G.; Huguet, S.; Vende, P.; Vautherot, J.-F.; Rasschaert, D. European brown hare syndrome virus: molecular cloning and sequencing of the genome. *J. Gen. Virol.* **1996**, *77*, 1693–1697.

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

17. Dalton, K.P.; Nicieza, I.; Balseiro, A.; Muguerza, M.A.; Rosell, J.M.; Casais, R.; Álvarez, Á.L.; Parra, F. Variant rabbit hemorrhagic disease virus in young rabbits, Spain. *Emerg. Infect. Dis.* **2012**, *18*, 2009–2012.

18. Abrantes, J.; Lopes, A.M.; Dalton, K.P.; Melo, P.; Correia, J.J.; Ramada, M.; Alves, P.C.; Parra, F.; Esteves, P.J. New variant of rabbit hemorrhagic disease virus, Portugal, 2012–2013. *Emerg. Infect. Dis.* **2013**, *19*, 1900–1902.

19. Le Gall-Reculé, G.; Lavazza, A.; Marchandeau, S.; Bertagnoli, S.; Zwingelstein, F.; Cavadini, P.; Martinelli, N.; Lombardi, G.; Guerin, J.-L.; Lemaitre, E.; et al. Emergence of a new lagovirus related to Rabbit Haemorrhagic Disease Virus. *Vet. Res.* **2013**, *44*, 81.

20. Westcott, D.; Frossard, J.-P.; Everest, D.; Dastjerdi, A.; Duff, J.P.; Steinbach, F.; Choudhury, B. Incursion of RHDV2-like variant in Great Britain. *Vet. Rec.* **2014**, *174*, 333.

21. Baily, J.L.; Dagleish, M.P.; Graham, M.; Maley, M.; Rocchi, M.S. RHDV variant 2 presence detected in Scotland. *Vet. Rec.* **2014**, *174*, 411.

22. Müller, A.; Freitas, J.; Silva, E.; Le Gall-Reculé, G.; Zwingelstein, F.; Abrantes, J.; Esteves, P.J.; Alves, P.C.; van der Loo, W.; Kolodziejek, J.; et al. Evolution of rabbit haemorrhagic disease virus (RHDV) in the European rabbit (*Oryctolagus cuniculus*) from the Iberian Peninsula. *Vet. Microbiol.* **2009**, *135*, 368–373.

23. Alda, F.; Gaitero, T.; Suárez, M.; Merchán, T.; Rocha, G.; Doadrio, I. Evolutionary history and molecular epidemiology of rabbit haemorrhagic disease virus in the Iberian Peninsula and Western Europe. *BMC Evol. Biol.* **2010**, *10*, 347–356.

24. Abrantes, J.; Lopes, A.M.; Esteves, P.J. Complete genomic sequences of rabbit hemorrhagic disease virus G1 strains isolated in the European rabbit original range. *J. Virol.* **2012**, *86*, 13886.

25. Abrantes, J.; Lopes, A.M.; Dalton, K.P.; Parra, F.; Esteves, P.J. Detection of RHDVa on the Iberian Peninsula: Isolation of an RHDVa strain from a Spanish rabbitry. *Arch. Virol.* **2014**, *159*, 321–326.

26. Ferrand, N.; Branco, M. The evolutionary history of the European rabbit (*Oryctolagus cuniculus*): Major patterns of population differentiation and geographic expansion inferred from protein polymorphism. In *Phylogeography of Southern European Refugia*; Weiss, S., Ferrand, N., Eds. Springer: Dordrecht, the Netherlands, 2007; p. 207.
27. Dalton, K.P.; Abrantes, J.; Lopes, A.M.; Nicieza, I.; Álvarez, Á.L.; Esteves, P.J.; Parra, F. Complete genome sequence of two variant RHDVb isolates detected on the Iberian Peninsula. *Arch. Virol.*, accepted.

28. Lopes, A.M.; Marques, S.; Silva, E.; Magalhães, M.J.; Pinheiro, A.; Alves, P.C.; Le Pendu, J.; Esteves, P.J.; Thompson, G.; Abrantes, J. Detection of RHDV strains in the Iberian hare (*Lepus granatensis*): Earliest evidence of rabbit lagovirus cross-species infection. *Vet. Res.* **2014**, *45*, 94.

29. Esteves, P.J.; Lopes, A.M.; Magalhães, M.J.; Pinheiro, A.; Gonçalves, D.; Abrantes, J. Rabbit hemorrhagic disease virus detected in Pico, Azores, Portugal, revealed a unique endemic strain with more than 17 years of independent evolution. *Viruses* **2014**, *6*, 2698–2707.

30. Hall, T.A. BioEdit: A user-friendly biological sequence alignment editor analysis program for Windows 95/98/NT. *Nucleic Acids Symp. Ser.* **1999**, *41*, 95–98.

31. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. *Mol. Biol. Evol.* **2013**, *30*, 2725–2729.

32. Delibes-Mateos, M.; Ferreira, C.; Carro, F.; Escudero, M.A.; Gortázar, C. Ecosystem effects of variant rabbit hemorrhagic disease virus, Iberian Peninsula. *Emerg. Infect. Dis.* **2014**, *20*, 2166–2168.

33. Calvete, C.; Sarto, P.; Calvo, A.J.; Monroy, F.; Calvo, J.H. Could the new rabbit haemorrhagic disease virus variant (RHDVb) be fully replacing classical RHD strains in the Iberian Peninsula? *World Rabbit Sci.* **2014**, *22*, 91.

34. Dalton, K.P.; Nicieza, I.; Abrantes, J.; Esteves, P.J.; Parra, F. Spread of new variant RHDV in domestic rabbits on the Iberian Peninsula. *Vet. Microbiol.* **2014**, *169*, 67–73.

35. Elsworth, P.; Cooke, B.D.; Kovaliski, J.; Sinclair, R.; Holmes, E.C.; Strive, T. Increased virulence of rabbit haemorrhagic disease virus associated with genetic resistance in wild Australian rabbits (*Oryctolagus cuniculus*). *Virology* **2014**, *464–465*, 415–423.

© 2014 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/4.0/).