The aim of this study was to characterise the nucleotide and amino acid sequence of complete genomes (7.5 kb) from RHDV strains isolated in Poland and estimate the genetic variability in different elements of the viral RNA. In addition, the sequence of Polish RHDV isolates isolated from 1988-2015 was compared with the sequences of other European RHDV, including the RHDVa and RHDV2/RHDVb subtypes. The complete sequence was developed by the compilation of partial nucleotide sequences. This sequence consisted of approximately 7428 nucleotides. For comparison of nucleotide sequences and the development of phylogenetic trees of Polish RHDV isolates and reference RHDV strains representing the main phylogenetic groups of classical RHDV, RHDVa and RHDV2 as well as the non-pathogenic rabbit lagovirus RCV, the BLAST software with blastn and MEGA6 with neighbour-joining method was applied. The complete nucleotide sequence of Polish isolates of RHDV has also been entered into GenBank. For comparative analysis, nineteen complete sequences representing the main RHDV genetic types available in GenBank were used. The results of phylogenetic analysis of Polish RHDV strains reveals the presence of three classical RHDV genogroups (G2, G4 and G5) and an RHDVa variant (G6). The oldest RHDV isolates (KGM 1988, PD 1989 and MAL 1994) belong to genogroup G2. It can be assumed that the elimination of these strains from the environment probably occurred at the turn of 1994 and 1995. Genogroup G2 was replaced by the phylogenetically younger BLA 1994 and OPO 2004 strains from genogroup G4, which probably originated from the G3 lineage, represented by the Italian strains BS89. The last representatives of classical RHDV in Poland are isolates GSK 1988 and ZD0 2000 from genogroup G5. A single clade contains the Polish RHDV strains from 2004-2015 (GRZ 2004, KRY 2004, L145 2004, W147 2005, SKO 2013, GLE 2013, RED1 2013, STR 2012, STR2 2013, STR 2014, BIE 2015) identified as RHDVa, which clustered into genogroup G6, as represented by the RHDV strain Triptis 1996. All recent isolated RHDV isolates belong exclusively to RHDVa and no RHDV2 was diagnosed in Poland.

European brown hare syndrome (EBHS) is one of the main causes of mortality in brown hares (Lepus europaeus) and mountain hares (Lepus timidus) in Europe. Since the mid-1990s, this highly lethal and contagious plague has been widespread in many European countries, contributing to a drastic decline in the number of free-living and farmed hares. A second lagovirus, able to infect some species of hares is rabbit haemorrhagic disease virus 2 (RHDV2; GI.2) recognised in 2010, a new viral emergence of RHDV (GI.1) which is known to be responsible for haemorrhagic disease in rabbits—RHD. The aim of this study was to evaluate the current EBHS epidemiological situation on the basis of the presence of antibodies to European brown hare syndrome virus (EBHSV) and anti-RHDV2 antibodies in sera collected from free-ranging hares in Central and Southeastern Poland in 2020–2021. Additionally, studies on the presence of EBHSV and RHDV2 antigens or their genetic material in the blood and internal organs taken from brown hares between 2014–2021 have been carried out. The results of the serological examination showed nearly 88% of tested blood samples were positive for EBHSV antibodies. No EBHSV was identified in the examined hares using virological and molecular tests. The positive results of EBHS serological studies confirmed the circulation and maintenance of EBHSV in free-living brown hares in Poland. However, no serological, virological or molecular evidence was obtained indicating that the brown hares tested had been in contact with RHDV2.

The aim of this study was the molecular epidemiology of independently introduced RHDV2 strains in Poland. The nucleotide sequences of RHDV2 diagnosed in domestic rabbits in 2018 in the voivodeships of Swietokrzyskie (strain PIN), Malopolskie (strain LIB) and Mazowieckie (strain WAK), and RHDVa from 2015 (strain F77-3) recognized in wild rabbits in Kujawsko-Pomorskie voivodeship were compared to the genome sequences of the first native RHDV2 strains from 2016–2017. The reference sequences available in public databases, the representative for a classical RHDV (G1-G5 genogroups), RHDVa (G6), non-pathogenic caliciviruses (RCV, GI.3 and GI.4) as well as original and recombinant RHDV2 isolates were included for this analysis. Nucleotide sequence similarity among the most distanced RHDV2 strains isolated in Poland in 2018 was from 92.3% to 98.2% in the genome sequence encoding ORF1, ORF2 and 3’UTR, between 94.8–98.7% in the VP60 gene and between 91.3-98.1% in non-structural proteins (NSP) region. The diversity between three RHDV2 and RHDVa from 2015 was up to 16.3% in the VP60 region. Similarities are shown for the VP60 tree within the RHDV2 group, however, the nucleotide analysis of NSP region revealed the differences between older and new native RHDV2 strains. The Polish RHDV2 isolates from 2016-2017 clustered together with RHDV G1/RHDV2 recombinants, first identified in the Iberian Peninsula in 2012, while all strains from 2018 are close to the original RHDV2. The F77-3 strain clustered to well supported RHDVa (G6) genetic group, together with other Polish and European RHDVa isolates. Based on the results of phylogenetic characterization of RHDV2 strains detected in Poland between 2016–2018 and the chronology of their emergence it can be concluded that RHDV2 strains of 2018 and RHDV2 strains of 2016–2017 were introduced independently thus confirming their different origin and simultaneous pathway of spreading.

Seroprevalence studies of RHDV antibodies in domestic rabbits were conducted between 2008-2014. A total of 12,169 sera from the provinces of central, southern and south-east Poland, including 7,570 samples collected from mixed-breed rabbits reared in smallholder farms and nearly 4,600 sera taken mainly from unvaccinated rabbits kept in industrial farms, were examined using ELISA tests. Additionally, cross-reactivity of selected tested and control archival sera using both classic RHDV and RHDVa antigens was determined by HI assay. The overall seroprevalence was 13.3%. In rabbits with unkown history of immunisation or RHD infection which came from small farms, RHDV antibodies were detected in 6.1% ranging between 1.0% to 17.2% of animals. In rabbits of the same group, but with a declared vaccination status, or confirmed exposure to an infectious virus, or coming from exposed females, the seroprevalence ranged from 83% to 100%. Among unvaccinated meat rabbits aged 71 to 90 days from industrial farms, low (1.85%, 4.17%, 11%), medium (34%, 54%) or high rates (98.7%) of seropositivity were detected. The seroconversion recorded in adult vaccinated females from industrial farms was 70% and 95%. Generally, the antibody levels examined by ELISAs and HI were comparable. However, a number of sera from the rabbits from small farms, as well as archival sera, showed clear differences. Several-fold differences in antibody titers, evidenced mainly in the postoutbreak sera, indictaed the contact of animals with RHDVa antigen. The overall results of the survey revealed a great proportion of seronegative rabbits potentially highly susceptible to RHD infection. In combination with the emergence of a novel pathogenic RHD virus type (RHDV2), it poses a severe risk of a next wave of fatal disease cases spreading in the native population of domestic rabbits, especially in farms with a traditional system of husbandry.

In order to determine the genetic variability of Polish RHD virus strains and to confirm the presence of genetic variant (RHDVa) subtype the partial nucleotide sequences of capsid protein gene, including two highly variable regions C and E, were examined. Phylogenetic analyses of 15 viral strains obtained over 18 years revealed the presence of three genetic groups. The oldest RHDV strains exhibit very close amino acid sequence similarity (98-99%) to the German FRG89 reference strain and most of European strains of the same period, as well as Chinese isolate from 1984. The HA-negative strains and isolates with variable reactivity in the HA test belong to the second subgroup and exhibit an intermediate level of variability (about 3%) in the analysed VP60 gene fragment. The most genetically variable strains (6-7%) clustered to RHDVa subtype. The analysis of nucleotides and amino acid sequences demonstrated three pairs of well conserved RHDV strains, isolated over 3, 6 and 10-year period.

The highly fatal rabbit hemorrhagic disease (RHD) that first emerged in 1984 in China has spread worldwide and affects both domestic and wild rabbits. The disease was originally caused by RHD virus (Lagovirus europaeus, L.europaeus) of GI.1 genotype, but over the years, two further pathogenic forms, known as the antigenic and genetic variant RHDVa (GI.1a) and RHDV2 (genotype GI.2), have been identified. RHD was first reported in Poland in 1988, when two RHDV strains were isolated, currently classified as GI.1c, while RHDVa and RHDV2 emerged in 2003 and 2016, respectively. In this study, using virological and molecular methods, we characterized five new RHDV strains belonging to GI.1 (RHDV)/GI.1a (RHDVa) and GI.2 (RHDV2) genotypes isolated in Poland in 2020–2022, in domestic rabbits from backyard farm and companion animals. We showed that two strains of L. europaeus (NRU 2020 and LIB 2020) from 2020 in the phylogenies of nonstructural proteins (NSP) and structural capsid protein (SP-VP60) clustered in a homogeneous GI.1a variant group. We stated that three strains of L. europaeus from 2020 to 2022 (KOB 2020, ZWO 2021, WAE 2022) in the VP60 phylogeny were positioned in the GI.2 (RHDV2) genotype, while in the NSP phylogeny, they are genetically related to recombinants with the GI.3/GI.2 genotype. Unexpectedly, in two RHD cases identified in the same small geographical area of south-eastern Poland (Libusza and Kobylanka), the close coexistence of RHDVa (LIB2020) and RHDV2 (KOB2020) strains capable of causing independent infections at the same time was found. This leads to the conclusion that the close natural coexistence of RHDV strains belonging to different genotypes does not necessarily have to directly lead to the emergence of new genetic or antigenic variants, which confirms the distinctness of both genetic forms and indicates different evolutionary paths leading to the best possible adaptation to the host.