European Brown Hare Syndrome in Poland: Current Epidemiological Situation.

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.

Phylogenetic Analysis of European Brown Hare Syndrome Virus Strains from Poland (1992-2004).

European brown hare syndrome (EBHS) is lethal to several species of free-living hares worldwide. The genetic characterization of its virus (EBHSV) strains in European circulation and epidemiological knowledge of EBHSV infections is not yet complete. The study determined the nucleotide sequences of the genomes of EBHSV strains from Poland and analyzed their genetic and phylogenetic relationships to a group of hare lagoviruses. The genome of five virus strains detected in Poland between 1992 and 2004 was obtained by RT-PCR and sequencing of the obtained amplicons. The genetic relationships of the EBHSV strains were analyzed using the full genome and VP60 gene sequences. Additionally, the amino acid sequence of the VP60 gene was analyzed to identify mutations specific to recognized EBHSV subgroups. Partial amplification of the virus open reading frame (ORF)1 and ORF2 regions obtained nearly complete nucleotide genome sequences of the EBHSV strains. Phylogenetic analysis placed them in a GII.1 cluster with other European strains related to nonpathogenic hare caliciviruses. VP60 gene analysis allocated these EBHSV strains to the G1.2, G2.2-2.3 or G3 virus genetic groups. The amino acid sequence differences in the entire genome ranged from 1.1 to 2.6%. Compared to a reference French EBHSV-GD strain, 22 variable amino acid sites were identified in the VP60 region of the Polish strains, but only six were in VP10. Single amino acid changes appeared in different sequence positions among Polish and other European virus strains from different genetic groups, as well as in VP10 sequences of nonpathogenic hare caliciviruses. The results of the study showed a high genetic homogeneity of EBHSV strains from Poland despite their different location occurrence and initial detection times. These strains are also phylogenetically closely related to other EBHSV strains circulating in Europe, likely confirming the slow evolutionary dynamics of this lagovirus species.

Evidence of independent introductions of RHDV2 strains in Poland based on the genome analysis of viral isolates from 2016-2018.

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.

Phylogenetic analysis of rabbit haemorrhagic disease virus (RHDV) strains isolated in Poland.

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.