Abstract

Hepatitis E virus (HEV) is an emerging zoonotic pathogen in Europe. Even though swine species are considered the main host of the zoonotic HEV-3 genotype, rabbits are recognized as the main reservoir of the divergent HEV-3ra subtype. However, the role of wild lagomorphs in the epidemiology of this virus in Mediterranean ecosystems is under debate. The aims of this study were to assess exposure of HEV in wild rabbit (Oryctolagus cuniculus) and Iberian hare (Lepus granatensis) populations in southern Spain and to determine potential risk factors associated with HEV exposure in these species. Between 2018 and 2021, blood and fecal samples from 370 wild rabbits and 60 Iberian hares were collected. A total of 29 (6.7%; 95%CI: 4.4–9.1) out of 430 sampled animals showed anti-HEV antibodies. By species, the seroprevalences in wild rabbit and Iberian hare were 6.8% (29/370; 95%CI: 4.2–9.3) and 6.7% (4/60; 95%CI: 0.4–13.0), respectively. Seropositive animals were detected on 17 (26.2%; 95%CI: 15.4–36.8) of the 65 sampled hunting estates. The generalized estimating equations model showed that geographical area was a risk factor potentially associated with HEV exposure in wild lagomorphs in the study region. HEV RNA was not detected in any of the 242 (0.0%; 95%CI: 0.0–1.5) fecal samples tested. This is the first large-scale serosurvey performed in wild rabbits in the Iberian Peninsula and in Iberian hares worldwide. Our results provide evidence of low, widespread, and heterogeneous distribution of HEV among wild rabbit and Iberian hare populations in Spanish Mediterranean ecosystems, which indicates a limited role of wild lagomorphs in the maintenance of the virus and a low risk of transmission of HEV to other species, including humans.

1. Introduction

Paslahepevirus balayani (also known as hepatitis E virus (HEV); family Hepeviridae) is the main cause of human acute viral hepatitis worldwide [1]. Eight different genotypes of HEV have been recognized so far (HEV-1 to 8), of which HEV-3, HEV-4 and HEV-7 are zoonotic. Among them, HEV-3 is the genotype with the greatest geographical distribution, being the most prevalent in many high-income regions, including Europe [2]. In this continent, HEV-3 human cases have considerably increased during the last few years, and contact with infected animals and the consumption of raw or undercooked animal products have been shown to be the major transmission routes of HEV infection [3]. Although wild boar (Sus scrofa) and domestic pig (Sus scrofa domesticus) are recognized as the main hosts of HEV-3 [4], rabbits, which are also susceptible to swine HEV-3 related strains [5, 6], represent the major reservoir of the divergent but also emergent and zoonotic HEV-3ra subtype [7].

Since the first description of HEV-3ra in domestic and wild rabbits in France in 2012 [8], HEV RNA and anti-HEV antibodies have been detected in lagomorphs in several European countries, including France [9], United Kingdom, Italy [6, 10], The Netherlands [11], Poland [12], and Germany [5, 1315] with prevalence and seroprevalence values ranged from 5.0% to 60.0% and from 2.6% to 37.3%, respectively. In the Iberian Peninsula, the European wild rabbit (Oryctolagus cuniculus) and the Iberian hare (Lepus granatensis) are two endemic species, being the most important small game lagomorph species in terms of abundance and hunting interest. These species are considered an important source of food for humans, especially in rural areas and usually for self-consumption, as rabbit meat is one of the most nutritional white meats [16]. Over the last decade, circulation of HEV has been confirmed in sympatric wildlife from this European region, such as wild boar (Sus scrofa) [17], red deer (Cervus elaphus) [18], and also the Iberian lynx (Lynx pardinus) [19], whose staple prey is the European wild rabbit. Even though active HEV infection was not found in wild lagomorphs in southern Spain [20], only a single local study has assessed HEV exposure in a limited number of wild lagomorphs in the Iberian Peninsula (southern Portugal), detecting seropositivity [21]. Therefore, the involvement of these species in the epidemiology of HEV in the Iberian Peninsula is under debate. The aims of the present study were to assess the role of wild lagomorphs in the epidemiology of HEV in Mediterranean ecosystems of southern Spain and to determine potential risk factors associated with HEV exposure in these species.

2. Material and Methods

2.1. Study Design and Sampling

A cross-sectional study was carried out in Andalusia (southern Spain; 36°N-38°60′N, 1°75′W-7°25′W), one of the most important regions in Spain in terms of small game population [22]. About 1.4 million wild rabbits and 251 thousand hares are harvested annually in Andalusia, being the second Spanish community with the highest number of lagomorphs hunted annually [22, 23]. The Andalusian Mediterranean ecosystems are characterized by the “dehesa” agroforestry system interspersed with Mediterranean forest, where different land uses such as agriculture, farming, and/or hunting are simultaneously exploited. The climate is a continental thermos Mediterranean climate with hot, dry summers and mild winters. The western region presents higher mean humidity and less extreme mean temperatures than the central and eastern regions.

Whenever possible, 60 wild rabbits were sampled per province in order to ensure a 95% probability of detecting at least one positive animal, assuming a minimum seroprevalence of 5% in each province [24]. In addition, for each province, sampling sites (hunting states) were randomly selected, and between 1 and 15 (mean: 8.6), wild rabbits were kindly provided by hunters in each of the hunting estates to collect samples (Figure 1). Between 2020 and 2021, serum and fecal samples were collected from a total of 370 wild rabbits sampled in 43 hunting estates. In addition, serum and feces from 60 Iberian hares were taken between 2018 and 2021 in 28 hunting estates in the same study region using a convenience sampling.

Blood samples from all animals were taken from the heart or thoracic cavity, and sera were then obtained after blood centrifugation at 400 for 10 min. Feces were collected directly from the rectum. Samples were stored at −80°C until laboratorial analysis.

Information about each animal, including location, sampling date, age, and sex, was recorded whenever possible (Table 1). Bodyweight and length were used as indicators of age [25]. Three groups of age were considered: young (<40 days old), subadult (from 40 days to 8 months), or adult (over 8 months). Epidemiological data related to hunting estates were also gathered through personal interviews with the gamekeepers using a standardized questionnaire (Table 1).

2.2. Laboratory Analysis

The presence of anti-HEV antibodies was assessed using a commercial double-antigen multispecies ELISA (HEV ELISA 4.0v; MP Diagnostics, Illkirch, France), in accordance with the manufacturer’s instructions. Samples were considered positive when the optical density (OD) at 450 nm of the sample was superior to the cut-off value (cut-off value = mean OD of the negative control + 0.2). This assay is based on the highly conserved recombinant protein ET2.1 [26] and detects the presence of total antibodies (IgM, IgG, and IgA) against the virus in sera or plasma in a wide range of animal species.

RNA from feces of wild lagomorphs from seropositive hunting states was extracted using the IndiSpin Pathogen Kit (Indical Biosciences, Germany) following the manufacturer’s instructions. For HEV RNA detection, two RT-PCR assays were carried out in parallel. A broad-spectrum real-time RT-PCR (CFX connect real-time PCR system) capable of detecting all P. balayani, genotypes, including HEV-3ra subtype, with the QIAGEN one-step RT-PCR kit (QIAGEN, Hilden, Germany) was used as previously described [27]. In addition, a nested broad-spectrum RT-PCR that is capable to amplify the four genera of hepevirus was also carried out according with Johne et al. [28]. For the first round, the QIAGEN one-step RT-PCR kit was used whereas the nested PCR was carried out with the premixed 2X solution Taq DNA polymerase, dNTPs, and reaction buffer kit (Promega). The second PCR products were examined on 1.5% agarose gel stained with RedSafe™ Nucleic Acid Staining Solution.

2.3. Statistical Analyses

The prevalence and seroprevalence against HEV was determined from the proportion of positive and seropositive animals, respectively, to the total number of examined, using the two-sided exact binomial test, 95% CI. Firstly, associations between results and independent variables were screened using the Pearson’s Chi-square or Fisher’s exact test, as appropriate. All variables with a in the bivariate analysis were selected for further analyses. Collinearity between pairs of variables was tested by Cramer’s V coefficient. Finally, a generalized estimating equation (GEE) analysis was carried out to study the effect of the variables selected based on the bivariate analysis. The number of positive animals was assumed to follow a binomial distribution, and “hunting estate” was included as the subject variable. The model was rerun until all remaining variables showed statistically significant values (). SPSS 25.0 software (Statistical Package for Social Sciences, Inc., Chicago, IL, USA) was used for all statistical analyses.

3. Results and Discussion

In recent decade, new viral strains of HEV have emerged in different European countries, which may have important clinical and epidemiological implications [2932]. Since 2012, an increasing number of HEV-3ra human cases has been confirmed in France, Belgium, Switzerland, and, recently, also in Spain [8, 29, 3335], and the high homology between HEV-3ra isolates from humans and rabbits supports the hypothesis of zoonotic transmission linked to consumption of these animals [3]. In addition, HEV-3ra seems to produce chronic hepatitis E cases in immunosuppressed patients more often than other subtypes [33]. In this context, understanding the role of lagomorphs in the epidemiology of HEV-3ra is an important key issue to control this emerging genotype.

The present study provides new epidemiological data on HEV in wild lagomorph populations in Mediterranean ecosystems of southern Spain. The major strengths of the survey include a large sample size, which is representative of the whole study area. A total of 29 (6.7%; 95%CI: 4.4–9.1) out of 430 sampled animals showed anti-HEV antibodies, which confirm that wild lagomorphs are naturally exposed to HEV in Iberian Mediterranean ecosystems. The distribution of seropositivity according to the explanatory variables is shown in Table 1. The seroprevalence in wild rabbit was 6.8% (29/370; 95% CI: 4.2–9.3). Our result is consistent with that previously found in this species in Australia (9.1%; 33/362) [36] and slightly higher to those reported in the United Kingdom (3.3%; 1/30) [6] and also in a local study conducted in southern Portugal, where 3 out of 74 showed anti-HEV antibodies (4.1%) [21]. By contrast, higher seropositivity values were obtained in wild rabbits in Japan (33.3%; 20/60) [37], Germany (30.8%; 4/13–37.3%; 47/126) [5, 13], Italy (42.9%; 15/35) [6], and Burkina Faso (60.0%; 60/100) [38].

In the present study, seropositivity was detected in four (6.7%; 95%CI: 0.4–13.0) of the 60 Iberian hares sampled. It should be noted that only a very few numbers of serosurveys have been conducted to date in hares worldwide. Two of them were carried out in Germany in European brown hare (Lepus europaeus) and found lower seroprevalences (2.2%; 14/624−2.6%; 25/944) [5, 15], whereas higher seropositivity values were found in Lepus africana in Burkina Faso (52.6%; 10/19) [38]. Comparisons between studies should be made with caution, given the differences in the serological methods used, study design and species, and number of animals sampled. Nevertheless, we would like to state that the seroprevalence of HEV in wild rabbits and Iberian hares in the study area should be considered low.

Statistical differences among ages were not found () (Table 1). In line, seropositivity was detected in 6.0% (7/117) of subadults and 12.8% (5/39) of young lagomorphs in all the years of the study period. Although the presence of maternal antibodies in yearling individuals cannot be ruled out, our results suggest endemic HEV circulation in wild lagomorph populations in Iberian Mediterranean ecosystems during the study period. Seropositive animals were found in all the provinces of the study area and in 17 (26.2%; 95%CI: 15.4–36.8) of the 65 sampled hunting estates (Figure 1).

Although variables related with the hunting estate (“geographical area,” “presence of wildcat,” “presence of goat,” “presence of swamps,” and “lagomorph density”) showed association () with the dependent variable in the bivariate analyses (Table 1), only the variable “geographical area” was retained in the GEE model. Significantly higher seropositivity was detected in eastern (11.2%; 95%CI: 6.9–15.5; ; OR = 6.5 [1.6–26.4]) compared to central Andalusia (2.2%; 95%CI: 0.0–4.7) but not with western region (3.4%; 95%CI: 0.0–7.1; ; OR = 1.9 [0.3–10.1]). These findings denote a widespread but not homogeneous distribution of HEV in wild rabbit and Iberian hare populations from southern Spain.

Experimental studies reported that HEV can be detected in lagomorphs’ feces between two- and ten-weeks postinfection [39]. Previous studies have confirmed natural circulation of HEV in wild lagomorphs with prevalence values that ranged between 1.7% and 16.0% [8, 11, 37]. In the present study, none of the 242 (0.0%; 95%CI: 0.0–1.5) feces analyzed were positive to HEV RNA, which indicate a limited active HEV circulation in wild lagomorphs in southern Spain. Our results are also consistent with the absence of infection reported previously in liver from these species in the same study region [20] and confirm a limited risk of transmission of HEV from wild lagomorphs to other sympatric species, including humans.

In conclusion, the serological results obtained in the present study indicate natural exposure to HEV but low, widespread, and heterogeneous viral circulation in wild rabbits and Iberian hare populations in southern Spain. The low seroprevalence as well as the absence of active HEV infection point that these wild lagomorphs play a limited role in the epidemiology of this virus in Mediterranean ecosystems of southern Spain. Further studies in different regions of the Iberian Peninsula are warranted to get a deeper and broader understanding of the epidemiological situation of HEV in these lagomorph species.

Data Availability

The data that support the findings of this study are available from the authors upon reasonable request.

Ethical Approval

No ethical approval was necessary.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors gratefully acknowledge the help of the personal of the Epidemiological Surveillance Program in Wildlife (Regional Government of Andalusia) in the collection of samples and epidemiological information. This work has benefited from the financial aid of research grants funded by the Spanish Ministry of Science and Innovation (project LagoHealth; REF: PID2019-111080RB-C21 and LAGOMED project REF: PRIMAS2-11-PCI2019-103698) and by the University of Córdoba (REF: UCO-FEDER-1264967). The Ministerio de Sanidad (RD12/0017/0012) integrated into the Plan Nacional de I+D+I and co-financed by the ISCIII-Subdirección General de Evaluación and the Fondo Europeo de Desarrollo Regional (FEDER); the Fundación para la Investigación en Salud (FIS) del Instituto Carlos III (PI19/00864). In addition, this research was supported by the CIBER-Consorcio Centro de Investigación Biomédica en Red (CB 2021/13/00083), Instituto de Salud Carlos III, Ministerio de Ciencia e Innovación and Unión Europea-NextGenerationEU. S. Castro-Scholten is supported by an FPU grant from the Spanish Ministry of Universities (FPU19/06026). A. Rivero-Juarez is the recipient of a Miguel Servet Research Contract by the Ministerio de Ciencia, Promoción y Universidades of Spain (CP18/00111). D. Jiménez-Martín holds a Ph.D. contract granted by Own Research Plan of the University of Cordoba. (Open Access funding enabled and organized by CRUE-CBUA Gold).