Exposure to Coxiella burnetii in Wild Lagomorphs in Spanish Mediterranean Ecosystems

Coxiella burnetii is an important zoonotic pathogen of worldwide distribution that can infect a wide range of wild and domestic species. The European wild rabbit (Oryctolagus cuniculus) can play a role as a reservoir for this bacterium in certain epidemiological scenarios, but, to date, a very limited numbers of large-scale serosurveys have been conducted for this species worldwide. Although exposure in hare species has also been described, C. burnetii in Iberian hare (Lepus granatensis) has never been assessed. Here, we aimed to determine the seroprevalence and risk factors associated with C. burnetii exposure in wild lagomorphs in the Mediterranean ecosystems of southern Spain. Between the 2018/2019 and 2021/2022 hunting seasons, blood samples from 638 wild lagomorphs, including 471 wild rabbits and 167 Iberian hares, were collected from 112 hunting grounds distributed across all eight provinces of Andalusia (southern Spain). The overall apparent individual seroprevalence was 8.9% (57/638; 95% CI: 6.8-11.4). Antibodies against C. burnetii were found in 11.3% (53/471; 95% CI: 8.4-14.1) of the wild rabbits and 2.4% (4/167; 95% CI: 0.1-4.7) of the Iberian hares. Seropositive animals were detected for 16 (14.3%; 95% CI: 7.8-20.8) of the 112 hunting grounds tested and in all the hunting seasons sampled. A generalized estimating equations model showed that the geographical area (western Andalusia) and presence of sheep were risk factors potentially associated with C. burnetii exposure in wild lagomorphs. A statistically significant spatial cluster (p < 0.001) was identified in the south-west of Andalusia. Our results provide evidence of moderate, endemic and heterogeneous circulation of C. burnetii in wild lagomorph populations in Spanish Mediterranean ecosystems. Risk-based strategies for integrative surveillance programs should be implemented in these species to reduce the risk of transmission of the bacterium to sympatric species, including humans.

Retrospective serological and molecular survey of myxoma or antigenically related virus in the Iberian hare, Lepus granatensis

The 2018 outbreak of myxomatosis in the Iberian hare (Lepus granatensis) has been hypothesized to originate from a species jump of the rabbit-associated myxoma virus (MYXV), after natural recombination with an unknown poxvirus. Iberian hares were long considered resistant to myxomatosis as no prior outbreaks were reported. To provide insights into the emergence of this recombinant virus (ha-MYXV), we investigated serum samples from 451 Iberian hares collected over two time periods almost two decades apart, 1994-1999 and 2017-2019 for the presence of antibodies and MYXV-DNA. First, we screened all serum samples using a rabbit commercial indirect ELISA (iELISA) and then tested a subset of these samples in parallel using indirect immunofluorescence test (IFT), competitive ELISA (cELISA) and qPCR targeting M000.5L/R gene conserved in MYXV and ha-MYXV. The cut-off of iELISA relative index 10 = 6.1 was selected from a semiparametric finite mixture analysis aiming to minimize the probability of false positive results. Overall, MYXV related-antibodies were detected in 57 hares (12.6%) including 38 apparently healthy hares (n = 10, sampled in 1994-1999, none MYXV-DNA positive, and n = 28 sampled in 2017-2019 of which four were also ha-MYXV-DNA positive) and 19 found-dead and ha-MYXV-DNA-positive sampled in 2018-2019. Interestingly, four seronegative hares sampled in 1997 were MYXV-DNA positive by qPCR, the result being confirmed by sequencing of three of them. For the Iberian hares hunted or live trapped (both apparently health), seroprevalence was significantly higher in 2017-2019 (13.0%, CI95% 9.2-18.2%) than in 1994-1999 (5.4%, CI95% 3.0-9.6%) (p = .009). Within the second period, seroprevalence was significantly higher in 2019 compared to 2017 (24.7 vs 1.7% considering all the sample, p = .007), and lower during the winter than the autumn (p < .001). While our molecular and serological results show that Iberian hares have been in contact with MYXV or an antigenically similar virus at least since 1996, they also show an increase in seroprevalence in 2018-2019. The remote contact with MYXV may have occurred with strains that circulated in rabbits, or with unnoticed strains already circulating in Iberian hare populations. This work strongly suggests the infection of Iberian hares with MYXV or an antigenically related virus, at least 20 years before the severe virus outbreaks were registered in 2018.

Spatiotemporal Overlap between the European Brown Hare and Its Potential Predators and Competitors

Analysis of spatiotemporal partitioning is pivotal to shed light on interspecific coexistence. Most research effort has involved large-sized carnivores and their prey, whereas little attention has been devoted to lagomorphs. We assessed spatiotemporal overlap among the European brown hare Lepus europaeus and its potential competitors and predators through camera-trapping in an area in Central Italy. We estimated the interspecific patterns of the spatiotemporal activity rhythms of brown hares, its potential predators (the red fox Vulpes vulpes, the pine marten Martes martes, the domestic cat Felis catus, and the domestic dog Canis familiaris), and a competitor, the roe deer Capreolus capreolus. Brown hare activity was studied in natural conditions as well as in a fenced area that excluded terrestrial predators and competitors. Free-ranging hares developed a more nocturnal behavior to avoid diurnal predators (i.e., domestic carnivores and martens). Although high temporal overlap was observed between free-ranging brown hares and both red foxes (82%) and roe deer (81%), hares avoided fox by being more active on darkest nights, as well as avoided roe deer through spatial partitioning. We suggest that hares may adapt their spatiotemporal behavior to avoid potential predators and competitors.

Zoonotic Microsporidia in Wild Lagomorphs in Southern Spain

Microsporidia are obligate intracellular protist-like fungal pathogens that infect a broad range of animal species, including humans. This study aimed to assess the presence of zoonotic microsporidia (Enterocytozoon bieneusi, Encephalitozoon intestinalis, Encephalitozoon hellem, and Encephalitozoon cuniculi) in organ meats of European wild rabbit (Oryctolagus cuniculus) and Iberian hare (Lepus granatensis) consumed by humans in Spain. Between July 2015 and December 2018, kidney samples from 383 wild rabbits and kidney and brain tissues from 79 Iberian hares in southern Spain were tested by species-specific PCR for the detection of microsporidia DNA. Enterocytozoon bieneusi infection was confirmed in three wild rabbits (0.8%; 95% CI: 0.0-1.7%) but not in hares (0.0%; 95% CI: 0.0-4.6%), whereas E. intestinalis DNA was found in one wild rabbit (0.3%; 95% CI: 0.0-0.8%) and three Iberian hares (3.8%; 95% CI: 0.0-8.0%). Neither E. hellem nor E. cuniculi infection were detected in the 462 (0.0%; 95% CI: 0.0-0.8%) lagomorphs analyzed. The absence of E. hellem and E. cuniculi infection suggests a low risk of zoonotic foodborne transmission from these wild lagomorph species in southern Spain. To the authors' knowledge, this is the first report of E. intestinalis infection in wild rabbits and Iberian hares. The presence of E. bieneusi and E. intestinalis in organ meats from wild lagomorphs can be of public health concern. Additional studies are required to determine the real prevalence of these parasites in European wild rabbit and Iberian hare.

Monitoring of emerging myxoma virus epidemics in Iberian hares (Lepus granatensis) in Spain, 2018-2020

Myxomatosis is an infectious disease caused by the myxoma virus (MYXV), which has very high mortality rates in European wild rabbits (Oryctolagus cuniculus). While sporadic cases of myxomatosis have also been reported in some hare species, these lagomorphs are considered to have a low susceptibility to MYXV infection. In the present study, we describe the spatiotemporal evolution and main epidemiological findings of novel hare MYXV (ha-MYXV or MYXV-Tol) epidemics in Iberian hares (Lepus granatensis) in Spain. In the period 2018-2020, a total of 487 hares from 372 affected areas were confirmed to be MYXV-infected by PCR. ha-MYXV outbreaks were detected in most of the Spanish regions where the Iberian hare is present. The spatial distribution was not homogeneous, with most outbreaks concentrated in the southern and central parts of Spain. Consecutive outbreaks reported in the last two years suggest endemic circulation in Spain of this emerging virus. A retrospective study carried out just after the first epidemic period (2018-2019) revealed that the virus could have been circulating since June 2018. The number of outbreaks started to rise in July, peaked during the first half of August and October and then decreased sharply until January 2019. The apparent mean mortality rate was 55.4% (median: 70%). The results indicated high susceptibility of the Iberian hare to ha-MYXV infection, but apparent resistance in the sympatric hare species present in Spain and less infectivity in European rabbits. The novel ha-MYXV has had significant consequences on the health status of Iberian hare populations in Spain, which is of animal health and conservation concern. The present study contributes to a better understanding of ha-MYXV emergence and will provide valuable information for the development of control strategies. Further research is warranted to assess the impact of this emerging virus on wild lagomorph populations and to elucidate its ecological implications for Iberian Mediterranean ecosystems.

Leishmania infantum in free-ranging hares, Spain, 2004-2010

Iberian hares (Lepus granatensis) were recently deemed responsible for an outbreak of human leishmaniasis affecting metropolitan Madrid, Spain. However, the reservoir potential of hares in Europe is poorly known. We report a retrospective survey on Leishmania infantum, the causal agent of zoonotic endemic leishmaniasis in the Mediterranean basin, infection status of Iberian, European (Le. europaeus) and Broom (Le. castroviejoi) hares in Spain. Spleen samples from 94 hares were tested by polymerase chain reaction. Sequencing and restriction fragment length polymorphism (RFLP) assays were performed on positive samples and RFLP patterns compared with those of strains reported in the scientific literature. DNA prevalence in hare spleen samples was 43.6% (95% confidence interval: 33.6-53.6). In all six regions studied at least one positive sample was found. RFLP revealed existence of specific hare strains of L. infantum differing from those reported in wild carnivores in Spain. The widespread presence of L. infantum in the most abundant Spanish hare species and the recent evidence of the ability of naturally infected hares to transmit the pathogen to Phlebotomus perniciosus, its main vector in the western Mediterranean, suggest that hares may have an unexpected role in the epidemiology of L. infantum in Spain.

Past, present and future distributions of an Iberian Endemic, Lepus granatensis: ecological and evolutionary clues from species distribution models

The application of species distribution models (SDMs) in ecology and conservation biology is increasing and assuming an important role, mainly because they can be used to hindcast past and predict current and future species distributions. However, the accuracy of SDMs depends on the quality of the data and on appropriate theoretical frameworks. In this study, comprehensive data on the current distribution of the Iberian hare (Lepus granatensis) were used to i) determine the species' ecogeographical constraints, ii) hindcast a climatic model for the last glacial maximum (LGM), relating it to inferences derived from molecular studies, and iii) calibrate a model to assess the species future distribution trends (up to 2080). Our results showed that the climatic factor (in its pure effect and when it is combined with the land-cover factor) is the most important descriptor of the current distribution of the Iberian hare. In addition, the model's output was a reliable index of the local probability of species occurrence, which is a valuable tool to guide species management decisions and conservation planning. Climatic potential obtained for the LGM was combined with molecular data and the results suggest that several glacial refugia may have existed for the species within the major Iberian refugium. Finally, a high probability of occurrence of the Iberian hare in the current species range and a northward expansion were predicted for future. Given its current environmental envelope and evolutionary history, we discuss the macroecology of the Iberian hare and its sensitivity to climate change.

Deviance partitioning of host factors affecting parasitization in the European brown hare (Lepus europaeus)

Deviance partitioning can provide new insights into the ecology of host-parasite interactions. We studied the host-related factors influencing parasite prevalence, abundance, and species richness in European brown hares (Lepus europaeus) from northern Spain. We defined three groups of explanatory variables: host environment, host population, and individual factors. We hypothesised that parasite infection rates and species richness were determined by different host-related factors depending on the nature of the parasite (endo- or ectoparasite, direct or indirect life cycle). To assess the relative importance of these components, we used deviance partitioning, an innovative approach. The explained deviance (ED) was higher for parasite abundance models, followed by those of prevalence and then by species richness, suggesting that parasite abundance models may best describe the host factors influencing parasitization. Models for parasites with a direct life cycle yielded higher ED values than those for indirect life cycle ones. As a general trend, host individual factors explained the largest proportion of the ED, followed by host environmental factors and, finally, the interaction between host environmental and individual factors. Similar hierarchies were found for parasite prevalence, abundance, and species richness. Individual factors comprised the most relevant group of explanatory variables for both types of parasites. However, host environmental factors were also relevant in models for indirect life-cycle parasites. These findings are consistent with the idea of the host as the main habitat of the parasite; whereas, for indirect life-cycle parasites, transmission would be also modulated by environmental conditions. We suggest that parasitization can be used not only as an indicator of individual fitness but also as an indicator of environmental quality for the host. This research underlines the importance of monitoring parasite rates together with environmental, population, and host factors.

Long-term changes in game species over a long period of transformation in the Iberian Mediterranean landscape

Agricultural change has transformed large areas of traditional farming landscapes, leading to important changes in the species community assemblages in most European countries. We suspect that the drastic changes in land-use that have occurred in Andalusia (southern Spain) over recent decades, may have affected the distribution and abundance of game species in this region. This article compares the distribution of the main game species in Andalusia during the 1960s and 1990s, using data from maps available from the Mainland Spanish Fish, Game and National Parks Service and from recent datasets on hunting yield distributions, respectively. Big-game and small-game species were significantly segregated in southern Spain during the 1990s, as two clearly independent chorotypes (groups of species whose abundances are similarly distributed) were obtained from the classification analysis. In contrast, big-game and small-game species were not significantly segregated several decades ago, when there was only one chorotype consisting of small-game species and wild boar. The other three ungulates did not constitute a significant chorotype, as they showed positive correlations with some species in the group mentioned above. These changes seem to be a consequence of the transformations that have occurred in the Iberian Mediterranean landscape over the last few decades. The abandoning of traditional activities, and the consequent formation of dense scrubland and woodland, has led to an expansion of big-game species, and a decrease of small-game species in mountain areas. Moreover, agricultural intensification has apparently depleted small-game species populations in some agricultural areas. On the other hand, the increasingly intensive hunting management could be artificially boosting this segregation between small-game and big-game species. Our results suggest that the conservation and regeneration of traditional agricultural landscapes (like those predominating in the 1960s) should be a priority for the conservation of small-game species.