Skip directly to site content Skip directly to page options Skip directly to A-Z link Skip directly to A-Z link Skip directly to A-Z link
Volume 19, Number 11—November 2013

New Variant of Rabbit Hemorrhagic Disease Virus, Portugal, 2012–2013

On This Page
Article Metrics
citations of this article
EID Journal Metrics on Scopus

Cite This Article

To the Editor: During November 2012–February 2013, rabbit hemorrhagic disease virus (RHDV) strains belonging to the new variant RHDV were isolated in Portugal from wild European rabbits (Oryctolagus cuniculus subsp. algirus). The major capsid protein, VP60, of these strains was partially characterized. RHDV had previously been detected in Portugal in 1989 (1). Before 2011, RHDV outbreaks in wild European rabbit (O. cuniculus) populations in the Iberian Peninsula were exclusively caused by strains belonging to genogroup 1 (2,3).

In the Iberian Peninsula, 2 subspecies of European rabbit are found, O. cuniculus subsp. algirus and O. cuniculus subsp. cuniculus. These subspecies are equally susceptible to RHDV (3). In 2011, a new variant was isolated in young rabbits belonging to O. cuniculus subsp. cuniculus from a rabbitry in the province of Navarra, Spain (4). The topology of the phylogenetic tree that included this variant and the susceptibility of kits <2 months old suggest that this strain is similar to that described in France in 2010 (5).

Before the new variant of RHDV emerged and, on the basis of phylogenetic relationships, RHDV strains had been divided into 6 genogroups (G1–G6) (1), with strains of G6, or RHDVa, having a distinct antigenic profile (6). All of these strains replicate in the liver and are responsible for causing death in rabbits >2 months of age. Nonpathogenic and weakly pathogenic RHDV-related strains have also been described. The nonpathogenic and weakly pathogenic strains are phylogenetically distinct from the G1–G6 strains with ≈20% of nucleotide divergence (7); they typically replicate in the intestines (8,9). New variant RHDV causes death in kits as young as 30 days old and affects vaccinated and unvaccinated animals (4). Phylogenetically, this new variant falls between the nonpathogenic groups (4,5).

During November 2012–February 2013, our laboratory, CIBIO, Universidade do Porto, Portugal, received liver samples from wild adult rabbits and kits, belonging to O. cuniculus subsp. algirus, from 3 areas of Portugal, Valpaços, Barrancos, and Algarve. The rabbits had appeared dead and had clinical signs suggesting rabbit hemorrhagic disease (RHD). We analyzed the samples for RHDV by reverse transcription PCR. For this process, total RNA was extracted by using the RNeasy Mini Kit (QIAGEN, Hilden, Germany), according to the manufacturer’s instructions. Reverse transcription was performed by using oligo(dT) as primer (Invitrogen, Carlsbad, CA, USA) and SuperScript III reverse transcription (Invitrogen) as recommended by the manufacturer. Screening of the samples consisted of PCR with a pair of primers as described by Dalton et al. (4). This pair amplifies a 738-bp fragment of the gene encoding the capsid protein, VP60 (PCR conditions are available on 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 same pair of primers.

The virus was detected in 15 samples, 5 from each locality. The obtained sequences were aligned with those available from public databases. Retrieved sequences represent the RHDV groups G1–G6, the nonpathogenic groups, and the new variant (GenBank accession nos. KF442960–KF442964). A phylogenetic tree was inferred in MEGA5 (10) by using a maximum-likelihood (ML) approach. Reliability of the nodes was assessed with a bootstrap resampling procedure consisting of 500 replicates of the ML trees. The best-fit nucleotide substitution model was determined by using MEGA5.


Thumbnail of Maximum-likelihood phylogenetic tree of 95 partial sequences of the rabbit hemorrhagic disease virus (RHDV) capsid gene. Bootstrap values appear next to the nodes and are shown only for the major groups: G1–G6 (GenBank accession nos. AB300693, AF231353, AF258618, AF453761, AJ006019, AJ302016, AJ303106, AJ495856, AJ535092, AJ535094, AJ969628, AM085133, AY269825, AY523410, AY926883, AY928268, AY928269, DQ069280, DQ069281, DQ069282, DQ189077, DQ189078, DQ205345, DQ280493, DQ530363, DQ8

Figure. . . Maximum-likelihood phylogenetic tree of 95 partial sequences of the rabbit hemorrhagic disease virus (RHDV) capsid gene. Bootstrap values appear next to the nodes and are shown only for the...

Our sequences exhibit the highest nucleotide sequence identity with the RHDV N11 strain from Spain (99%; GenBank accession no. JX133161.1), which corresponds to the new RHDV variant. Thirteen nucleotide substitutions were detected in comparison to the Spanish sequence, 3 of which were nonsynonymous. The inferred ML phylogenetic tree is in agreement with those published (1,3,9). G1–G6 (pathogenic) RHDV strains and nonpathogenic and weakly pathogenic RHDV-related strains (generally referred to as RCV) form 2 groups (Figure). The nonpathogenic strain from Australia (RCV-A1_Australia_MIC-07) does not cluster with other nonpathogenic groups and European brown hare syndrome virus (EBHSV_France) appears in a basal position in the tree. As described, the new variant (N11_Spain) appears between RCV and the nonpathogenic Australian strain (4,5). The strains isolated from rabbits in Portugal cluster with the new variant and form a highly supported group (bootstrap value 1.00). These results support the conclusion that the virus recovered in Portugal belongs to the new variant RHDV described in Spain and France.

This confirms the presence of the virus in wild rabbits on the Iberian Peninsula. We also confirm that both European rabbit subspecies are susceptible to the new variant. The appearance and rapid spread of the new variant RHDV into the Iberian wild rabbit populations raise concern for the survival of these populations in this region. These conservation concerns are particular highlighted for the O. cuniculus subsp. algirus, because it only occurs in the southwestern part of the Iberian Peninsula, and it is a key prey species for several carnivores, namely, for the most endangered feline, the Iberian Lynx (Lynx pardinus). Therefore, monitoring the spread and evolution of this new variant is crucial in determining the most appropriate conservation measures.



The samples from Valpaços and Algarve were provided by Jorge Pires and Vitor Palmilha, respectively.

The Portuguese Foundation for Science and Technology supported the doctoral fellowship of A.M.L. (SFRH/BD/78738/2011) and the postdoctoral fellowship of J.A. (SFRH/BPD/73512/2010). The Portuguese Foundation for Science and Technology Projects PTDC/CVT/108490/2008 and FCT-ANR/BIA-BIC/0043/2012 supported this work. Also, project “Genomics Applied to Genetic Resources,” cofinanced by North Portugal Regional Operational Programme 2007/2013 (ON.2–O Novo Norte), under the National Strategic Reference Framework, through the European Regional Development Fund, supported this work. F.P. and K.P.D. gratefully acknowledge the financial support of Organización Interprofesional Cunicola.


Joana Abrantes, Ana M. Lopes, Kevin P. Dalton, Pedro Melo, Jorge J. Correia, Margarida Ramada, Paulo C. Alves, Francisco Parra, and Pedro J. EstevesComments to Author 
Author affiliations: Centro de Investigação em Biodiversidade e Recursos Genéticos/InBIO Laboratório Associado (CIBIO), Vairão, Portugal (J. Abrantes, A.M. Lopes, P.C. Alves, P.J. Esteves); ;Université de Nantes, Nantes, France (J. Abrantes, A.M. Lopes); ; Faculdade de Ciências da Universidade do Porto, Porto, Portugal (A.M. Lopes, P.C. Alves); Universidad de Oviedo, Asturias, Spain (K.P. Dalton, F. Parra); VetNatura, Lisbon, Portugal (P.M., M. Ramada); ; Faculdade de Medicina Veterinária–Universidade Técnica de Lisboa, Lisbon (J.J. Correia); University of Montana, Missoula Montana, USA (P.C. Alves); Instituto de Investigação e Formação Avançada em Ciências e Tecnologias da Saúde (CESPU), Gandra, Portugal (P.J. Esteves)



  1. Abrantes  J, van der Loo  W, Le Pendu  J, Esteves  PJ. Rabbit haemorrhagic disease (RHD) and rabbit haemorrhagic disease virus (RHDV): a review. Vet Res. 2012;43:12. DOIPubMedGoogle Scholar
  2. Alda  F, Gaitero  T, Suarez  M, Merchan  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. DOIPubMedGoogle Scholar
  3. Muller  A, Freitas  J, Silva  E, Le Gall-Recule  G, Zwingelstein  F, Abrantes  J, Evolution of rabbit haemorrhagic disease virus (RHDV) in the European rabbit (Oryctolagus cuniculus) from the Iberian Peninsula. Vet Microbiol. 2009;135:36873. DOIPubMedGoogle Scholar
  4. Dalton  KP, Nicieza  I, Balseiro  A, Muguerza  MA, Rosell  JM, Casais  R, Variant rabbit hemorrhagic disease virus in young rabbits, Spain. Emerg Infect Dis. 2012;18:200912 . DOIPubMedGoogle Scholar
  5. Le Gall-Reculé  G, Zwingelstein  F, Boucher  S, Le Normand  B, Plassiart  G, Portejoie  Y, Detection of a new variant of rabbit haemorrhagic disease virus in France. Vet Rec. 2011;168:1378. DOIPubMedGoogle Scholar
  6. Capucci  L, Fallacara  F, Grazioli  S, Lavazza  A, Pacciarini  ML, Brocchi  E. A further step in the evolution of rabbit hemorrhagic disease virus: the appearance of the first consistent antigenic variant. Virus Res. 1998;58:11526. DOIPubMedGoogle Scholar
  7. Abrantes  J, Esteves  PJ. Not-so-novel Michigan rabbit calicivirus. Emerg Infect Dis. 2010;16:13312. DOIPubMedGoogle Scholar
  8. Capucci  L, Fusi  P, Lavazza  A, Pacciarini  ML, Rossi  C. Detection and preliminary characterization of a new rabbit calicivirus related to rabbit hemorrhagic disease virus but nonpathogenic. J Virol. 1996;70:861423 .PubMedGoogle Scholar
  9. Strive  T, Wright  JD, Robinson  AJ. Identification and partial characterisation of a new Lagovirus in Australian wild rabbits. Virology. 2009;384:97105. DOIPubMedGoogle Scholar
  10. Tamura  K, Peterson  D, Peterson  N, Stecher  G, Nei  M, Kumar  S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:27319. DOIPubMedGoogle Scholar




Cite This Article

DOI: 10.3201/eid1911.130908

Reference has only first page number. Please provide the last page number if article is longer than one page. (in reference 1 "Abrantes, van der Loo, Le Pendu, Esteves, 2012").

Reference has only first page number. Please provide the last page number if article is longer than one page. (in reference 2 "Alda, Gaitero, Suarez, Merchan, Rocha, Doadrio, 2010").

Related Links


Table of Contents – Volume 19, Number 11—November 2013

EID Search Options
presentation_01 Advanced Article Search – Search articles by author and/or keyword.
presentation_01 Articles by Country Search – Search articles by the topic country.
presentation_01 Article Type Search – Search articles by article type and issue.



Please use the form below to submit correspondence to the authors or contact them at the following address:

Pedro J. Esteves, CIBIO-UP, Centro de Investigação em Biodiversidade e Recursos Genéticos, Universidade do Porto, InBIO, Laboratório Associado, Campus Agrário de Vairão, Rua Padre Armando Quintas, Vairão, Portugal

Send To

10000 character(s) remaining.


Page created: October 31, 2013
Page updated: October 31, 2013
Page reviewed: October 31, 2013
The conclusions, findings, and opinions expressed by authors contributing to this journal do not necessarily reflect the official position of the U.S. Department of Health and Human Services, the Public Health Service, the Centers for Disease Control and Prevention, or the authors' affiliated institutions. Use of trade names is for identification only and does not imply endorsement by any of the groups named above.