Wild ungulates as sentinels of flaviviruses and tick-borne zoonotic pathogen circulation: an Italian perspective
BMC Veterinary Research volume 19, Article number: 155 (2023)
Vector-borne zoonotic diseases are a concerning issue in Europe. Lyme disease and tick-borne encephalitis virus (TBEV) have been reported in several countries with a large impact on public health; other emerging pathogens, such as Rickettsiales, and mosquito-borne flaviviruses have been increasingly reported. All these pathogens are linked to wild ungulates playing roles as tick feeders, spreaders, and sentinels for pathogen circulation. This study evaluated the prevalence of TBEV, Borrelia burgdorferi sensu lato, Rickettsia spp., Ehrlichia spp., and Coxiella spp. by biomolecular screening of blood samples and ticks collected from wild ungulates. Ungulates were also screened by ELISA and virus neutralization tests for flaviviral antibody detection.
A total of 274 blood samples were collected from several wild ungulate species, as well as 406 Ixodes ricinus, which were feeding on them. Blood samples tested positive for B. burgdorferi s.l. (1.1%; 0-2.3%) and Rickettsia spp. (1.1%; 0-2.3%) and showed an overall flaviviral seroprevalence of 30.6% (22.1–39.2%): 26.1% (17.9–34.3%) for TBEV, 3.6% (0.1–7.1%) for Usutu virus and 0.9% (0-2.7%) for West Nile virus. Ticks were pooled when possible and yielded 331 tick samples that tested positive for B. burgdorferi s.l. (8.8%; 5.8–11.8%), Rickettsia spp. (26.6%; 21.8–31.2%) and Neoehrlichia mikurensis (1.2%; 0-2.4%). TBEV and Coxiella spp. were not detected in either blood or tick samples.
This research highlighted a high prevalence of several tick-borne zoonotic pathogens and high seroprevalence for flaviviruses in both hilly and alpine areas. For the first time, an alpine chamois tested positive for anti-TBEV antibodies. Ungulate species are of particular interest due to their sentinel role in flavivirus circulation and their indirect role in tick-borne diseases and maintenance as Ixodes feeders and spreaders.
Vector-borne diseases are described as emerging infections due to their massive spread in recent decades. Tick-borne diseases (TBDs) and mosquito-borne diseases represent the most threatening vector-borne infections worldwide . The rise of vector populations, together with the pathogens they may carry, is a deeply studied topic, but there are many aspects yet to be revealed [2,3,4,5].
The diffusion of TBDs is driven by several factors, mainly consequent to human activities. Modern agriculture has led to changes in land use leading to an increase in mountain abandonment, usually followed by natural reforestation, which is a frequently reported phenomenon in the Eastern Alps [6, 7]. In turn, land abandonment stimulates an increase in wild animals, which can boost tick presence .
Wild ungulates are abundant in northern Italy [8, 9] and are known as preferential hosts for tick feeding and reproduction, especially for Ixodes ricinus (I. ricinus) ticks [9, 10]. In Europe, as well as in northern Italy, I. ricinus is one of the most studied tick species due to its high competence in the transmission of zoonotic pathogens, including viruses, bacteria, and parasites, to humans [1, 11,12,13].
The most prevalent zoonotic agents transmitted by I. ricinus in Europe belong to the Borrelia burgdorferi sensu lato complex . Other zoonotic and emerging pathogens transmitted by I. ricinus bite are tick-borne encephalitis virus (TBEV), Rickettsia spp., Anaplasma phagocytophilum (A. phagocytophilum), Ehrlichia spp., and Neoehrlichia spp. [15,16,17,18]. Currently, Borrelia burgdorferi sensu latu (B. burgdorferi s.l.) is the main zoonotic tick-borne pathogen in northern Italian regions , while TBEV , A. phagocytophilum, and Rickettsia spp.  have been only sporadically detected in humans, whereas Ehrlichia spp. and Neoehrlichia spp. . have not yet been identified.
The epidemiological cycles of these infections involve wild ungulates that mainly act as tick feeders and spreaders in the environment, as well as reservoir hosts of A. phagocytophilum . In contrast, I. ricinus acts as both a vector and reservoir for Rickettsia spp. , and with regard to the other infectious agents, different vertebrate species serve as reservoir hosts [21, 23, 24].
Coxiella burnetii (C. burnetii) has been identified in I. ricinus, although its vector capacity is still debated [25, 26]. It may cause disease in humans and livestock , and similar to domestic ruminants, red deer and wild boar are susceptible to C. burnetii and have been suggested to play a role in its maintenance .
Flaviviruses are currently the most widespread vector-borne zoonotic viruses in Europe and are transmitted by hard ticks (TBEV) or mosquitos (Usutu virus-USUV and West Nile virus-WNV). Both USUV and WNV are considered endemic in Italy first being detected in black birds (1996) and horses (1998), respectively [29, 30]. The epidemiological cycle of these flaviviral species includes Culex mosquitos and wild birds. Mammals, including humans, can be infected but are considered dead-end hosts [31,32,33].
Regarding TBEV, Italy is considered a low-risk country, and most cases of human infection are found in the northeastern regions . The main vector is I. ricinus , but the detection of TBEV in ticks is challenging, time-consuming and expensive since the prevalence is restricted to small foci of viral circulation . In areas where TBEV incidence is high in humans, prevalence in ticks rarely exceeds 1% [35, 36]. Recently, wild ungulates have been considered TBEV sentinels, noting that they seroconvert but usually do not develop disease, even though a fatal case has been recently described in a roe deer [34, 37,38,39]. In European countries, TBEV antibodies were found in several wild ungulate species [34, 40,41,42], with seroprevalence ranging between 0.74% in Finland and 63.5% in Poland. Despite the comparable clinical and epidemiological scenarios, fewer studies have investigated the seroprevalence in ungulates for WNV and USUV; seroconversion in asymptomatic animals was reported, suggesting the role of ungulates as sentinels in viral diffusion [36, 43, 44].
The involvement of wild ungulates in the epidemiology of all the abovementioned vector-borne pathogens is noticeable, and although several of these etiological agents are considered endemic to the alpine area of northern Italy, data derived from ungulate populations are currently scarce.
This research aimed to (i) assess the prevalence of a zoonotic virus (TBEV) and zoonotic bacteria (B. burgdorferi s.l., Rickettsia spp., Ehrlichia spp., C. burnetii) in wild ungulate blood and associated tick samples using biomolecular tools; (ii) identify the selected species and genospecies through sequencing; (iii) investigate the sentinel role of ungulates in flaviviral (TBEV, WNV, USUV) circulation by antibody detection; and (iv) evaluate any the statistical significance of the obtained prevalence and seroprevalence related to several variables of interest.
The zoonotic bacterium A. phagocytophilum was purposely excluded from this study since it was previously investigated in the same ungulate and vector population .
Sample collection and molecular biology
A total of 281 blood samples were collected from culled wild ungulates in three different investigated regions: Lombardy (n = 105), Veneto (n = 67), and Friuli Venezia Giulia (n = 109). The species of ungulates included roe deer (Capreolus capreolus, n = 108), red deer (Cervus elaphus, n = 87), wild boar (Sus scrofa, n = 59), mouflon (Ovis orientalis musimon, n = 18), and chamois (Rupicapra rupicapra, n = 9). The DNA/RNA internal control was not detected in seven samples, which were excluded from the study; therefore, 274 blood samples were investigated for tick-borne pathogens (TBPs). Additional details are provided in Table 1.
Three out of 274 (1.1%; 0-2.3%) blood samples collected from two red deer (from Lombardy) and a wild boar (from Veneto) were positive for B. burgdorferi s.l.. Due to the low quantity of bacterial DNA (Cq > 35), none of the samples were confirmed by PCR and sequencing. Three out of 274 (1.1%; 0-2.3%) blood samples collected from a roe deer (from Lombardy), a red deer and a wild boar (from Friuli Venezia Giulia) tested positive for Rickettsia spp. Sequencing analysis confirmed the detection of an unidentified Rickettsia spp. from roe deer blood, with only 96% identity with known rickettsiae. Furthermore, R. helvetica was identified from red deer, while positivity of the wild boar blood was not confirmed by PCR. All blood samples tested negative for Ehrlichia spp., Coxiella spp., and TBEV.
A total of 406 ticks were collected and extracted in pools of two ticks (n = 75) or individually (n = 256). Thus, the total number of screened tick samples was 331 (Table 2). All ticks were identified as I. ricinus and most of the tick samples consisted of adult females (n = 231) or males (n = 90), while few immature stages were found, i.e., nymphs (n = 6) and larvae (n = 4). Most of the ticks were collected from roe deer, red deer, and mouflons. All tick samples tested positive for the DNA/RNA internal control, showing no PCR inhibition.
Twenty-nine out of 331 (8.8%; 5.8–11.8%) tick samples tested positive for B. burgdorferi s.l. Amplicon sequencing identified the following species: Borrelia garinii (B. garinii) (7/29), Borrelia afzelii (B. afzelii) (10/29), Borrelia burgdorferi sensu stricto (B. burgdorferi s.s.) (1/29) and Borrelia valaisiana (B. valaisiana) (2/29); however, some amplicons from tick samples (9/29) did not yield good quality sequences, and thus, genospecies identification was not possible.
A total of 88 out of 331 (26.6%; 21.8–31.4%) of the tick samples were positive for Rickettsia spp. Two zoonotic species were identified, namely, Rickettsia monacensis (R. monacensis) (35/88) and Rickettsia helvetica (R. helvetica) (36/88). Candidatus Rickettsia mendelii (Ca. R. mendelii) was found in two samples from ticks feeding on roe deer.
Four out of 331 tick samples (1.2%; 0-2.4%) collected from roe deer tested positive for Neoehrlichia mikurensis (N. mikurensis). No other Ehrlichia species were identified.
The prevalence of B. burgdorferi s.l. and Rickettsia spp. in tick samples did not show any significant difference when considering the ungulate species from which the ticks were collected. All tick samples tested negative for Coxiella spp. and TBEV. More details and data on the host species from which positive ticks were collected are shown in Table 2.
Ten out of 331 samples (3%; 1.2–4.9%) were positive for more than one pathogen. The presence of concurrent bacterial species was confirmed by Sanger sequencing in 6 of 10 specimens: coinfection between B. afzelii and R. helvetica was found in 2 pooled samples, one pool of male ticks and one of females, and in 3 individually extracted ticks, one female and two males. One of the males also tested positive for a third pathogen, N. mikurensis. A different male tick tested positive for both B. afzelii and R. monacensis.
ELISA and virus neutralization for flaviviruses
A total of 111 serum samples originated from the Veneto (n = 46) and Friuli Venezia Giulia (n = 65) regions; no serum samples were available from Lombardy. All ungulate species were represented, namely, roe deer (n = 45), red deer (n = 31), wild boar (n = 23), chamois (n = 8), and mouflon (n = 4).
In total, 34 out of 111 blood samples tested positive by ELISA showing an overall anti-flavivirus seroprevalence of 30.6% (22.1–39.2%); in detail, 24 samples were positive, and 10 were borderline. All sampled species were positive, except for mouflon. The detailed seroprevalence of the tested species is shown in Table 3.
The virus neutralization test (VNT) confirmed that most of the ELISA-positive sera were positive for TBEV (29/34), resulting in an overall prevalence of 26.1% (17.9–34.3%). Four ungulate species had TBEV antibodies: chamois, wild boar, roe deer and red deer (Table 3). Positive animals were mainly from the alpine and prealpine areas of both investigated regions. Of note, some were shot in a hilly area (altitude is approximately 170–180 m above see level) close to the flatlands. In addition, a roe deer tested positive for both TBEV and WNV.
Of the five TBEV-negative samples, one was negative for all investigated flaviviruses, while the remaining four were USUV positive. These were wild boars culled in a hilly area of the Friuli Venezia Giulia region; three were culled between June and October of the same year in different municipalities but in a range of 10–15 km2. Flavivirus antibody positivity did not show any statistical association with “species”, “sex”, “age”, and “season” variables. Detailed results of ELISA and VNT tests are available in Additional file 3.
The present study evaluated the presence and frequency of several vector-borne pathogens using molecular biology and serology. Except for C. burnetii, all investigated pathogens were identified through direct or indirect testing, highlighting the wide diffusion of both tick-borne and mosquito-borne infections in northern Italy.
Overall, the medical and veterinary importance of I. ricinus was confirmed by its high vector occurrence and remarkably frequent infection with zoonotic pathogens. According to the present results, wild ruminants (except for chamois) appear to be highly infested compared to wild boars, but they seem to play a minor role as reservoirs for tick-borne infectious agents. On the other hand, since most of the feeding I. ricinus were adults, the role of ungulates as tick amplifiers was confirmed. In addition, ungulates appear to be promising sentinels of flavivirus infections.
B. burgdorferi s.l. was detected in both blood (1.1%) and tick (8.8%) samples, although the low amount of Borrelia DNA (high Cq) in blood samples hampered further characterization. The low levels of bacteria are in accordance with previous investigations that outlined how wild ungulates can act as Borrelia dilution hosts, and thus, blood positivity is only sporadically identified without any apparent epidemiological relevance . However, the role of artiodactyls is controversial: on the one hand, they represent an important feeding source for ticks; on the other hand, they have shown a limited contribution to the transmission of the Lyme spirochete, decreasing the spread of infection [45, 46]. Moreover, I. ricinus tick samples showed a noticeable prevalence of Borrelia spp. and several zoonotic species were identified. The detection of four genospecies (B. garinii, B. afzelii, B. burgdorferi s.s., and B. valaisiana) indirectly highlights the presence of suitable reservoirs in the study area, such as wild birds and small mammals, and confirms the endemicity of these pathogens in the study area. A previous study  conducted in northeastern Italy on questing ticks reported the same zoonotic species, although at a lower prevalence, and similar identifications were found in northwestern Italy [47, 48]. The prevalence reported herein is also higher than that described in Spain (2.3%) and Poland (3.3%) [49, 50], likely due to the diversity of the ecological niches in different areas.
Evidence of Ehrlichia, particularly N. mikurensis, was limited to the vectors, with a prevalence of 1.2%, in agreement with a previous study conducted in a nearby area on questing ticks . Despite the recent discovery of N. mikurensis, it seems that rodents could be reservoirs, while wild ungulates may act as tick spreaders rather than amplifiers .
In contrast, Rickettsia spp. was detected in both blood samples and feeding ticks. The prevalence in blood was low (1.1%), and the detection of R. helvetica (in red deer) and an unidentified Rickettsia (in roe deer) is likely to be an occasional finding because wild ruminants are not considered reservoirs of infection . In contrast, the prevalence in tick samples was high (26.6%). Most Rickettsia spp. were R. monacensis and R. helvetica, which are typically found in inland and continental areas . A lower (Slovakia – 6.8%) and similar (Poland – 26.8%) prevalence in I. ricinus collected from wild ungulates was reported by other authors, reflecting the heterogeneity between habitats [12, 50]. Although these rickettsial species are considered pathogenic only in immunocompromised patients, a recent case described the onset of disease in an immunocompetent patient in Portugal .
Two ticks were positive for Ca. R. mendelii, which was first identified in 2016 in I. ricinus in Eastern Europe  and has been reported in I. ricinus questing ticks in Poland and the Czech Republic and in feeding ticks on migratory birds in Italy. The zoonotic potential of Ca. R. mendelii is still unknown [15, 55, 56].
C. burnetii was not found in blood and tick samples, whereas a prevalence of approximately 5% had been reported both in questing and wildlife-collected I. ricinus [25, 26]; other in vivo experimental studies demonstrated the shedding of C. burnetii in I. ricinus faeces, suggesting its potential role as a reservoir in the wild . In the studied area, I. ricinus and wild ungulates do not appear to have a pivotal role in C. burnetii maintenance in the wild. To date, their role in the sylvatic cycle has not yet been clarified: however, several cow and goat farms were found to be positive in the same regions [57, 58].
Similarly, all blood and tick samples were negative for TBEV based on direct identification, in line with the role of wild ungulates as tick hosts instead of reservoirs . Moreover, direct investigation in ticks is a sensitive tool for TBEV identification only when testing high amounts of specimens since the prevalence is usually low. Similar results have been reported by other authors in the same area , where the prevalence in questing ticks was 0.21% after testing more than 2300 samples . Therefore, other more cost-effective methods should be considered to assess the risk of new foci. The present study shows the importance of serological surveys, which yielded a high flavivirus seroprevalence (i.e., 30.6%). Most of the samples (26.1%) were positive for antibodies against TBEV (Table 3) and were collected from chamois, red deer, roe deer and wild boar, revealing a high seroprevalence of TBEV in the investigated area. As reported by other authors, wild ungulates can be considered suitable sentinels for TBEV circulation [41, 60]. TBEV antibodies were found - at a lower percentage - in many artiodactyls all over Europe: in roe deer from 2.1 to 22.9% [40, 43, 60, 61]; in wild boar from 5.6 to 20% [40, 43, 62]; in moose at 0.74%; in white-tailed deer at 0.74% ; and in red deer at 1.4% . Only Krzysiak et al. (2021) found a higher seroprevalence in European bison (62.7%) .
Comparing the different seroprevalence values, a great influence of the investigated area seems in place. Knowing that TBEV foci are localized and not homogeneous, a different pattern between and within nations can be expected. In fact, when sampling a larger area, a lower prevalence is expected, and when sampling smaller areas, high variability is likely [41, 60]. The present study focused on two of the most affected areas in Italy for TBEV; thus, a higher seroprevalence was expected compared to other Italian regions. However, the investigated region is considered a low-risk area in the European context, and the remarkable seroprevalence described is in contrast with the scarce human reports, suggesting probable under reporting of cases .
Imhoff et al. (2015) highlighted some critical aspects of the use of TBEV serological techniques as screening methods, such as haemolysis of sera (related to the shot) and the scarce precision of the geographical data, due to the wide foraging area of wild ungulates . Despite haemolysis in some samples, antibody detection was confirmed by two different diagnostic tests: ELISA and VNT, the gold standard test. As proposed by other authors, precise testing of TBEV antibodies in wild ungulates could be useful to establish risk maps in areas where data based only on human incidence could be biased, considering the high human TBEV vaccination coverage [60, 61]. Indeed, finding positive animals in areas where no human cases were previously reported would be important for the identification of other potential risk areas [24, 60]. Of note, TBEV seroconversion in alpine chamois has never been described in Europe before. This alpine species could be relevant for the early detection of new TBEV foci as an indirect sentinel of positive ticks in areas where other ungulate species are less frequently observed.
The results of the VNTs also highlighted the presence of USUV and WNV antibodies, detected in four wild boars (USUV – 3.6%) and one roe deer (WNV – 0.9%). The preference of Culex mosquitos for wild boars rather than ruminants, when present in the same area, was thus confirmed as already proposed . USUV was known to be present in the studied region: in fact, a human case has been recently described , but it had never been identified in wild boars before . Few studies researched USUV antibodies in wild boars and found similar (3.4% in Serbia) or higher (8% in France) seroprevalences [43, 64].
Regarding WNV antibodies in roe deer, both higher (23.5% in Serbia, 4.8% in Czech Republic) and lower seroprevalence (0% in Spain) were described [44, 64, 65]. The positive roe deer identified in the present study were found in Belluno Province, where until 2020, no human cases were reported . West Nile disease cases were mainly located along the Po Valley and the description of WNV at northern latitudes may reflect the effects of climate change.
When dealing with vector-borne diseases, all epidemiological data may be useful to better understand their prevalence and diffusion. Emerging pathogens are often identified in vectors and/or animals first, and only later the same pathogens are diagnosed in humans [29, 33, 35, 66, 67]. Emerging infectious agents such as N. mikurensis and Ca. R. mendelii were identified, but they have never been diagnosed in humans in the studied area; thus, their presence should be acknowledged going forward. Tick sampling from wild ungulates has proven effective for tick-borne pathogen surveillance, as well as serological surveys on wild ungulates, due to the remarkable susceptibility of these wild species to flaviviral infections. Further epidemiological studies should take these aspects into account, since wild ungulate monitoring could serve as an early warning system for the detection of viral diffusion in areas considered lacking but at risk of vector, both mosquito and tick, expansion.
The study area encompassed the Friuli Venezia Giulia (FVG), Veneto (V), and Lombardy (L) regions, located in northern Italy (Fig. 1). Prealpine and alpine areas located in the provinces of Udine (UD), Belluno (BL) and Varese (VA) were investigated. These areas are approximately located at latitudes between 45°-46° N and are characterized by a wide altimetric excursion, ranging from 180 m a.s.l. in hilly areas up to more than 3000 m a.s.l. of the highest mountain peaks. The climate of the area under study is characterized by the climatic features of the alpine region with relevant temperature excursions between seasons in relation to different altitudes. Higher altitudes are characterized by continental weather: winter months are cold and snowy, while the temperature is mild during warm seasons. In the prealpine area, the climate is overall milder, and during the hunting seasons, the temperatures range from 0 °C to approximately 30 °C, with a maximum humidity reaching 70–98% in several periods (data extrapolated from ARPA Veneto, ARPA FVG and kindly provided by ARPA Lombardia) [68,69,70].
In addition, the alpine region is facing continuous and dramatic changes due to global warming, which highlights not only a change in temperature but also several other modifications regarding snow cover, humidity, precipitation, vegetation and natural hazards [71, 72].
Sampling procedures were carried out between May 2017 and September 2020 during licenced hunting seasons. Blood samples were collected by hunters from wild ungulates in 9 mL Vacumed® tubes with K3EDTA (FL Medical srl, Italy) in hunting check stations or directly in the field. Ungulate carcasses were carefully inspected for ectoparasites by hunters and/or a veterinarian. When present, ticks were collected up to a maximum of 10 specimens per carcass, and each of them was placed in a 1.5 mL sterile microtube. An anamnestic form was filled out for each animal recording species, age, sex, date and place of culling, weight, health condition, and presence/absence of ticks. Both blood and tick samples were stored at 4 °C and sent to the Laboratory of Microbiology and Infectious Diseases of the Animal Medicine, Production and Health Department (Legnaro, Italy). Once in the laboratory, blood samples were divided into 200 µL aliquots while the ectoparasites were first morphologically identified using the identification keys of Manilla and Cringoli [74, 75] and then placed, individually or in pools of a maximum of two specimens, in a new 1.5 mL sterile microtube. Pools included ticks from the same host with the same features (species, sex, stage), and a low engorgement stage. Both blood samples and ectoparasites were stored at -80 °C until processing.
Nucleic acids were extracted from 200 µL of whole blood and tick samples using the All Prep DNA/RNA Mini Kit (QIAGEN GmbH, Germany) following the manufacturer’s instructions. Each tick sample was ground thoroughly with a sterile disposable plastic pestle in a 1.5 mL microtube, resuspended in 350 µL of Buffer RLT Plus, and homogenized by repetitive pipetting.
Before the lysis step during nucleic acid extraction, both DNA and RNA internal controls, supplied by the Quantinova Pathogen + IC kit (QIAGEN GmbH, Germany), were added to all specimens to assess both extraction efficiency and the presence of PCR inhibitors. Extracted nucleic acids were stored at -80 °C (RNA) or -20 °C (DNA) and then screened using previously published PCR and real-time PCR methods to detect the following viral and bacterial pathogens: TBEV , B. burgdorferi s.l. [77, 78], Rickettsia spp. [79, 80], Coxiella spp. , Ehrlichia spp.  and N. mikurensis [83, 84]. Positive and negative controls were included in each run. Details about the methods and procedures are provided in Additional file 1.
Internal control detection and pathogen screening were performed using a Quantinova Pathogen + IC kit (QIAGEN GmbH, Germany) on a LightCycler96 Instrument (Roche, Switzerland) with the Internal Control Assay kit and genus-specific real-time PCR assays (Additional file 1). N. mikurensis HRM real-time PCR assays were performed using 5x HOT FIREPol EvaGreen qPCR Mix Plus (Solis Biodyne) on a MyGo Pro instrument (IT-IS, United Kingdom). Ehrlichia spp. endpoint PCR screening was performed using 1× Phire Hot Start II PCR Master Mix (Thermo Fischer Scientific Baltics, Lithuania) on a Biometra TGradient thermal cycler (Analytic Jena GmbH, Germany).
All samples yielding a positive signal in both the internal control assay and pathogen screening were further investigated by specific end-point PCR assays (Additional file 1), followed by Sanger sequencing of amplicons. PCRs were performed with the same reagents and instrument used for Ehrlichia spp. screening.
PCR products were visualized by electrophoresis on 2% agarose gels stained with SybrSafe DNA Stain (Invitrogen by Thermo Fischer Scientific, USA) and subsequently purified using ExoSap-IT Express PCR Product Cleanup (Thermo Fischer Scientific Baltics, Lithuania) according to the manufacturer’s instructions. Bidirectional Sanger sequencing of all specific PCR products was carried out at the StarSEQ® GmbH facilities (Mainz, Germany) using the same PCR primers. Nucleotide sequences were assembled and edited using ChromasPro v.2.1.8 (Technelysium Pty Ltd, Australia) and were then deposited in GenBank (Accession numbers in Additional file 2) and analysed using the Nucleotide BLAST  search engine (National Center for Biotechnology Information, Bethesda, MD).
ELISA and virus neutralization for flaviviruses
Serum samples were collected in 9 mL Vacumed® tubes without anticoagulants (FL Medical srl, Italy). To assess the seroprevalence of TBEV, serum samples were screened with the two-step ELISA test Immunozym FSME IgG all species (PROGEN, Biotechnik GmbH, Heidelberg, Germany). According to the manufacturer’s recommendations, sera were diluted 1:50. The results of the test were expressed in Vienna International Units and samples were considered positive with > 126 Vienna units/mL, borderline if between 63 and 126 Vienna units/mL, and negative when < 63 Vienna units/mL. Positive and borderline samples were further tested with the gold standard confirmatory test, i.e., the virus neutralization test. VNTs were performed at the National Reference Laboratory for Arboviruses of the Ostrava Public Health Institute (Ostrava, Czech Republic). Anti-TBEV VNTs were performed using sterile 96-well plates. The TBEV strain (Hypr) was cultivated in intracerebrally infected suckling mice, and PS cells (porcine stable kidney cell line) were used as the susceptible cell line; VNTs were performed and the results were expressed as previously described  with minor changes consisting of the use of 25 µL of PS cell suspension (600,000 cells per mL) instead of the CV-1 cell line (African green monkey kidney fibroblasts). For anti-WNV and anti-USUV tests, an identical procedure was used, choosing the CV-1 cell line as the susceptible cell line for both viruses. WNV lineage 2 and USUV lineage Eur3 were used as virus suspensions. In each VNT, the endpoint titre was assessed as the higher serum dilution that inhibited the viral cytopathic effect. Samples showing a titre of anti-Flavivirus antibodies equal to 1:8 or higher were considered positive.
Data regarding sample collection and laboratory analyses were organized in a database on a Microsoft Excel Worksheet and descriptive statistics (counts, percentage and CI95%) were used to summarize results according to pathogen detection or serological results with respect to “species”, “sex”, and “age” of the animals, and “season” and “region” of the collection variables; the statistical analyses to detect significant differences in infection rates were conducted by means of chi-square test or Fisher’s exact test, when appropriate. The level of statistical significance was set for alpha = 0.05.
All data generated or analysed during this study are included in this research article [and its additional files 1, 2, and 3]. Further information may be available from the corresponding author upon request.
- A. phagocytophilum :
- B. burgdorferi s.l.:
Borrelia burgdorferi sensu lato
- B. garinii :
- B. afzelii :
- B. burgdorferi s.s.:
Borrelia burgdorferi sensu stricto
- B. valaisiana :
- Ca. R. mendelii:
Candidatus Rickettsia mendelii
- C. burnetii :
Friuli Venezia Giulia
- I. ricinus :
- m a.s.l.:
metres above sea level
- N. mikurensis :
- R. monacensis :
- R. helvetica :
tick-borne encephalitis virus
virus neutralization test
West Nile virus
Boulanger N, Boyer P, Talagrand-Reboul E, Hansmann Y. Ticks and tick-borne diseases. Med Mal Infect. 2019;49:87–97.
Daszak P, Cunningham AA, Hyatt AD. Emerging infectious Diseases of Wildlife - Threats to biodiversity and human health. Sci Compass. 2000;287:443–9.
Dantas-Torres F. Climate change, biodiversity, ticks and tick-borne diseases: the butterfly effect. Int J Parasitol Parasites Wildl. 2015;4:452–61.
Ostfeld RS, Brunner JL. Climate change and Ixodes tick-borne diseases of humans. Philos Trans R Soc B. 2015;370.
Cunningham AA, Daszak P, Wood JLN. One health, emerging infectious diseases and wildlife: two decades of progress? Philos Trans R Soc B. 2017;372:20160167.
Tasser E, Walde J, Tappeiner U, Teutsch A, Noggler W. Land-use changes and natural reforestation in the Eastern Central Alps. Agric Ecosyst Environ. 2007;118:115–29.
Battaglini L, Bovolenta S, Gusmeroli F, Salvador S, Sturaro E. Environmental sustainability of Alpine livestock farms. Ital J Anim Sci. 2014;13(3155):431–43.
Apollonio M, Andersen R, Putman R. European ungulates and their management in the 21st Century. Cambridge, UK: Cambridge University Press; 2010.
Carnevali L. In: Pedrotti L, Riga F, Toso S, editors. Ungulates in Italy - Status, distribution, abundance, management and hunting of ungulate populations in Italy. Volume Rapporto 2. ISPRA (Istituto Superiore per la Protezione e Ricerca Ambientale); 2009. pp. 1–168.
Gray J, Kahl O, Zintl A. What do we still need to know about Ixodes ricinus? Ticks and Tick-borne Diseases. 2021;12(101682).
Otranto D, Dantas-Torres F, Giannelli A, Latrofa MS, Cascio A, Cazzin S et al. Ticks infesting humans in Italy and associated pathogens. Parasites and Vectors. 2014;7(328).
Kazimírová M, Hamšíková Z, Špitalská E, Minichová L, Mahríková L, Caban R, et al. Diverse tick-borne microorganisms identified in free-living ungulates in Slovakia. Parasit Vectors. 2018;11:495.
Medlock JM, Hansford KM, Bormane A, Derdakova M, Estrada-Peña A, George J-C et al. Driving forces for changes in geographical distribution of Ixodes ricinus ticks in Europe. Parasit Vectors. 2013;6(1).
Wolcott KA, Margos G, Fingerle V, Becker NS. Host association of Borrelia burgdorferi sensu lato: a review. Ticks Tick Borne Dis. 2021;12(101766).
Guccione C, Colomba C, Tolomeo M, Trizzino M, Iaria C, Cascio A. Rickettsiales in Italy. Pathogens. 2021;10(181).
Höper L, Skoog E, Stenson M, Grankvist A, Wass L, Olsen B, et al. Vasculitis due to Candidatus Neoehrlichia mikurensis: a cohort study of 40 swedish patients. Clin Infect Dis. 2021;73(7):e2372–8.
Stuen S, Granquist EG, Silaghi C. Anaplasma phagocytophilum - A widespread multi-host pathogen with highly adaptive strategies. Front Cell Infect Microbiol. 2013;3:31.
Michelitsch A, Wernike K, Klaus C, Dobler G, Beer M. Exploring the Reservoir Hosts of Tick-Borne Encephalitis Virus. Viruses. 2019;11(669).
Trevisan G, Cinco M, Trevisini S, Di Meo N, Chersi K, Ruscio M et al. Borreliae Part 1: Borrelia Lyme Group and Echidna-Reptile Group. Biology. 2021;10(1036).
Riccò M. Epidemiology of tick-borne encephalitis in north-eastern italy (2017–2020): International insights from national notification reports. Acta Biomed. 2021;92(5):e2021229.
Grassi L, Franzo G, Martini M, Mondin A, Cassini R, Drigo M, et al. Ecotyping of Anaplasma phagocytophilum from wild ungulates and ticks shows circulation of zoonotic strains in northeastern Italy. Animals. 2021;11:310.
Tomassone L, Portillo A, Nováková M, De Sousa R, Oteo JA. Neglected aspects of tick-borne rickettsioses. Parasit Vectors. 2018;11(263).
Bertola M, Montarsi F, Obber F, Da Rold G, Carlin S, Toniolo F et al. Occurrence and identification of Ixodes ricinus borne pathogens in northeastern Italy. Pathogens. 2021;10(1181).
Alfano N, Tagliapietra V, Rosso F, Ziegler U, Arnoldi D, Rizzoli A. Tick-borne encephalitis foci in northeast Italy revealed by combined virus detection in ticks, serosurvey on goats and human cases. Emerg Microbes Infect. 2020;9:474–84.
Špitalská E, Sparagano O, Stanko M, Schwarzová K, Zdenko Å, Škultéty Ľ, et al. Diversity of Coxiella-like and Francisella-like endosymbionts, and Rickettsia spp., Coxiella burnetii as pathogens in the tick populations of Slovakia, Central Europe. Ticks Tick Borne Dis. 2018;9:1207–11.
Knap N, Žele2 D, Glinšek Biškup U, Avšič-Županc T, Vengušt G. The prevalence of Coxiella burnetii in ticks and animals in Slovenia. BMC Vet Res. 2019;15(368).
Yessinou ER, Mertens-scholz K, Neubauer H, Farougou S. Prevalence of Coxiella-infections in ticks - review and meta-analysis. Ticks Tick Borne Dis. 2022;13:101926.
González-Barrio D, Ruiz-Fons F. Coxiella burnetii in wild mammals: a systematic review. Transbound Emerg Dis. 2019;66(2):662–71.
Autorino GL, Battisti A, Deubel V, Ferrari G, Forletta R, Giovannini A, et al. West Nile virus Epidemic in Horses, Tuscany Region, Italy. Emerg Infect Dis. 2002;8(12):1372–8.
Weissenböck H, Bakonyi T, Rossi G, Mani P, Nowotny N. Usutu virus, Italy, 1996. Emerg Infect Dis. 2013;19(2):274–7.
Zecchin B, Fusaro A, Milani A, Schivo A, Ravagnan S, Ormelli S, et al. The central role of Italy in the spatial spread of USUTU virus in Europe. Virus Evol. 2021;7(1):veab048.
Vilibic-Cavlek T, Petrovic T, Savic V, Barbic L, Tabain I, Stevanovic V et al. Epidemiology of Usutu Virus: the european scenario. Pathogens. 2020;9(699).
Riccò M, Peruzzi S, Balzarini F. Epidemiology of West Nile Virus Infections in humans, Italy, 2012–2020: a summary of available evidences. Trop Med Infect Dis. 2021;6:61.
Krzysiak MK, Anusz K, Konieczny A, Rola J, Salat J, Strakova P, et al. The european bison (Bison bonasus) as an indicatory species for the circulation of tick-borne encephalitis virus (TBEV) in natural foci in Poland. Ticks Tick Borne Dis. 2021;12:101799.
Holding M, Dowall SD, Medlock JM, Carter DP, Pullan ST, Lewis J, et al. Tick-borne encephalitis virus, United Kingdom. Emerg Infect Dis. 2020;26(1):90–6.
Imhoff M, Hagedorn P, Schulze Y, Hellenbrand W, Pfeffer M, Niedrig M. Review: sentinels of tick-borne encephalitis risk. Ticks Tick Borne Dis. 2015;6:592–600.
Rizzoli A, Hauffe HC, Tagliapietra V, Neteler M, Rosà R. Forest structure and roe deer abundance predict tick-borne encephalitis risk in Italy. PLoS ONE. 2009;4(2):e4336.
Rosà R, Tagliapietra V, Manica M, Arnoldi D, Hauffe HC, Rossi C, et al. Changes in host densities and co-feeding pattern efficiently predict tick-borne encephalitis hazard in an endemic focus in northern Italy. Int J Parasitol. 2019;49:779–87.
Da Rold G, Obber F, Monne I, Milani A, Ravagnan S, Toniolo F et al. Clinical Tick-Borne Encephalitis in a Roe deer (Capreolus capreolus L.). Viruses. 2022;14(300).
Potkonjak A, Petrović T, Ristanović E, Lalić I, Vračar V, Savić S, et al. Molecular Detection and Serological evidence of Tick-Borne Encephalitis Virus in Serbia. Vector-Borne Zoonotic Dis. 2017;17(12):813–20.
Paulsen KM, das Neves CG, Granquist EG, Madslien K, Stuen S, Pedersen BN, et al. Cervids as sentinel-species for tick-borne encephalitis virus in Norway - A serological study. Zoonoses Public Health. 2020;67:342–51.
Tonteri E, Jokelainen P, Matala J, Pusenius J, Vapalahti O. Serological evidence of tick-borne encephalitis virus infection in moose and deer in Finland: sentinels for virus circulation. Parasites and Vectors. 2016;9(54).
Bournez L, Umhang G, Faure E, Boucher J-M, Boué F, Jourdain E et al. Exposure of wild ungulates to the Usutu and Tick-borne encephalitis viruses in France in 2009–2014: evidence of undetected Flavivirus circulation a decade ago. Viruses. 2020;12(10).
García-Bocanegra I, Paniagua J, Gutiérrez-Guzmán AV, Lecollinet S, Boadella M, Arenas-Montes A et al. Spatio-temporal trends and risk factors affecting West Nile virus and related flavivirus exposure in spanish wild ruminants. BMC Vet Res. 2016;12(249).
Fabri ND, Sprong H, Hofmeester TR, Heesterbeek H, Donnars BF, Widemo F et al. Wild ungulate species differ in their contribution to the transmission of Ixodes ricinus-borne pathogens. Parasites and Vectors. 2021;14(360).
Hofmeester T, Coipan E, Evan Wieren S, Prins H, Takken W, Sprong H. Few vertebrate species dominate the Borrelia burgdorferi s.l. life cycle. Environ Res Lett. 2016;11:043001.
Garcia-Vozmediano A, Krawczyk AI, Sprong H, Rossi L, Ramassa E, Tomassone L. Ticks climb the mountains: ixodid tick infestation and infection by tick-borne pathogens in the Western Alps. Ticks Tick Borne Dis. 2020;11(5).
Pistone D, Pajoro M, Novakova E, Vicari N, Gaiardelli C, Viganò R, et al. Ticks and bacterial tick-borne pathogens in Piemonte region, Northwest Italy. Exp Appl Acarol. 2017;73:477–91.
Díaz-Cao JM, Adaszek Ł, Dzięgiel B, Paniagua J, Caballero-Gómez J, Winiarczyk S, et al. Prevalence of selected tick-borne pathogens in wild ungulates and ticks in southern Spain. Transbound Emerg Dis. 2022;69(3):1084–94.
Michalski MM, Kubiak K, Szczotko M, Dmitryjuk M. Tick-borne pathogens in ticks collected from wild ungulates in North-Eastern Poland. Pathogens. 2021;10:587.
Portillo A, Santibáñez P, Palomar AM, Santibáñez S, Oteo JA. Candidatus Neoehrlichia mikurensis ’ in Europe. New Microbes New Infect. 2018;22:30–6.
Stefanidesova K, Kocianova E, Boldis V, Kostanova Z, Kanka P, Nemethova D, et al. Evidence of Anaplasma phagocytophilum and Rickettsia helvetica infection in free-ranging ungulates in central Slovakia. Eur J Wildl Res. 2008;54(3):519–24.
de Sousa R, dos Santos ML, Cruz C, Almeida V, Garrote AR, Ramirez F, et al. Rare case of Rickettsiosis caused by Rickettsia monacensis, Portugal, 2021. Emerg Infect Dis. 2022;28(5):1068–71.
Hajduskova E, Literak I, Papousek I, Costa FB, Novakova M, Labruna MB, et al. Candidatus Rickettsia mendelii, a novel basal group rickettsia detected in Ixodes ricinus ticks in the Czech Republic. Ticks Tick Borne Dis. 2016;7(3):482–6.
Stańczak J, Biernat B, Racewicz M, Zalewska M, Matyjasek A. Prevalence of different Rickettsia spp. in Ixodes ricinus and Dermacentor reticulatus ticks (Acari: Ixodidae) in north-eastern Poland. Ticks Tick Borne Dis. 2018;9(2):427–34.
Pajoro M, Pistone D, Boccazzi IV, Mereghetti V, Bandi C, Fabbi M et al. Molecular screening for bacterial pathogens in ticks (Ixodes ricinus) collected on migratory birds captured in northern Italy. Folia Parasitol (Praha). 2018;65(008).
Ceglie L, Guerrini E, Rampazzo E, Barberio A, Tilburg JJHC, Hagen F, et al. Molecular characterization by MLVA of Coxiella burnetii strains infecting dairy cows and goats of north-eastern Italy. Microbes Infect. 2015;17:776–81.
Barberio A. Coxiella burnetii infection in dairy cows and goats: assessment of diagnostic methods, and evaluation of immune response in shedders. PhD Dissertation Thesis. 2015; Università degli studi di Milano, Facoltà di Medicina Veterinaria, Dipartimento di Scienze veterinarie e Sanità pubblica (Accessed on 29/03/2023:https://air.unimi.it/retrieve/handle/2434/352272/515810/phd_unimi_R10106.pdf).
D’Agaro P, Martinelli E, Burgnich P, Nazzi F, Del Fabbro S, Iob A, et al. Prevalence of Tick-Borne Encephalitis Virus in Ixodes Ricinus from a Novel Endemic Area of North Eastern Italy. J Med Virol. 2009;81:309–16.
Duscher GG, Wetscher M, Baumgartner R, Walder G. Roe deer sera used for TBE surveillance in Austria. Ticks Tick Borne Dis. 2015;6(4):489–93.
Kiffner C, Vor T, Hagedorn P, Niedrig M, Rühe F. Determinants of tick-borne encephalitis virus antibody presence in roe deer (Capreolus capreolus) sera. Med Vet Entomol. 2012;26:18–25.
Balling A, Plessow U, Beer M, Pfeffer M. Prevalence of antibodies against tick-borne encephalitis virus in wild game from Saxony, Germany. Ticks Tick Borne Dis. 2014;5:805–9.
Caracciolo I, Mora-Cardenas E, Aloise C, Carletti T, Segat L, Burali MS, et al. Comprehensive response to usutu virus following first isolation in blood donors in the Friuli Venezia Giulia region of Italy: development of recombinant NS1-based serology and sensitivity to antiviral drugs. PLoS Negl Trop Dis. 2020;14(3):e0008156.
Escribano-Romero E, Lupulović D, Merino-Ramos T, Blázquez AB, Lazić G, Lazić S, et al. West Nile virus serosurveillance in pigs, wild boars, and roe deer in Serbia. Vet Microbiol. 2015;176:365–9.
Hubálek Z, Juricová Z, Straková P, Blazejová H, Betásová L, Rudolf I. Serological survey for West Nile Virus in Wild Artiodactyls, Southern Moravia (Czech Republic). Vector-Borne Zoonotic Dis. 2017;17(9):654–7.
Estrada-Peña A, Palomar AM, Santibáñez P, Sánchez N, Habela MA, Portillo A et al. Crimean-Congo Hemorrhagic Fever Virus in Ticks, Southwestern Europe, 2010. Emerg Infect Dis. 2012;18(1).
Negredo A, de la Calle-Prieto F, Palencia-Herrejón E, Mora-Rillo M, Astray-Mochales J, Sánchez-Seco MP, et al. Autochthonous Crimean–Congo Hemorrhagic Fever in Spain. N Engl J Med. 2017;377(2):154–61.
ARPA FVG - Osservatorio metereologico regionale del Friuli Venezia Giulia. Accessed on 29/03/23:https://www.clima.fvg.it/archivio.php?ln=&p=dati
ARPA Veneto - Agenzia Regionale per la Prevenzione e Protezione Ambientale del Veneto. Accessed on 29/03/23: https://www.arpa.veneto.it/dati-ambientali/dati-storici.
ARPA Lombardia - Agenzia Regionale per la Protezione dell’Ambiente. Accessed on 29/03/23: https://www.arpalombardia.it/Pages/Meteorologia/Richiesta-dati-misurati.aspx#.
Mourey J, Marcuzzi M, Ravanel L, Pallandre F. Effects of climate change on high Alpine mountain environments: evolution of mountaineering routes in the Mont Blanc massif (Western Alps) over half a century. Arct Antarct Alp Res. 2019;51(1):176–89.
Gobiet A, Kotlarski S, Beniston M, Heinrich G, Rajczak J, Stoffel M. 21st century climate change in the european Alps - a review. Sci Total Environ. 2014;493:1138–51.
Map chart. https://www.mapchart.net/italy.html (Accessed 12 June 2023).
Manilla G. Ixodida, Acari (Fauna d’Italia). Bologna - Italy: Edizioni Calderini; 1998. pp. 1–280.
Cringoli G, Iori A, Rinaldi L, Veneziano V, Genchi C. Zecche - Mappe parassitologiche. Rolando editore, Naples, Italy. 2005. 1–263 p.
Schwaiger M, Cassinotti P. Development of a quantitative real-time RT-PCR assay with internal control for the laboratory detection of tick borne encephalitis virus (TBEV) RNA. J Clin Virol. 2003;27(2):136–45.
Courtney JW, Kostelnik LM, Zeidner NS, Massung RF. Multiplex real-time PCR for detection of Anaplasma phagocytophilum and Borrelia burgdorferi. J Clin Microbiol. 2004;42(7):3164–8.
Lee SH, Lee JH, Park HS, Jang WJ, Koh SE, Yang YM, et al. Differentiation of Borrelia burgdorferi sensu lato through groEL gene analysis. FEMS Microbiol Lett. 2003;222(1):51–7.
Rolain JM, Bitam I, Buffet S, Marié JL, Bourry O, Portelli-Clerc C, et al. Presence or absence of plasmid in Rickettsia felis depending on the source of fleas. Eur Soc Clin Microbiol Infect Dis. 2009;15(Suppl 2):296–7.
Choi YJ, Jang WJ, Kim JH, Ryu JS, Lee SH, Park KH, et al. Spotted fever group and typhus group rickettsioses in humans, South Korea. Emerg Infect Dis. 2005;11(2):237–44.
Klee SR, Tyczka J, Ellerbrok H, Franz T, Linke S, Baljer G, et al. Highly sensitive real-time PCR for specific detection and quantification of Coxiella burnetii. BMC Microbiol. 2006;6:1–8.
Guillemi EC, Tomassone L, Farber MD. Tick-borne Rickettsiales: molecular tools for the study of an emergent group of pathogens. J Microbiol Methods. 2015;119:87–97.
Krücken J, Schreiber C, Maaz D, Kohn M, Demeler J, Beck S, et al. A novel high-resolution melt PCR assay discriminates Anaplasma phagocytophilum and Candidatus Neoehrlichia mikurensis. J Clin Microbiol. 2013;51(6):1958–61.
Goodman JL, Nelson C, Vitale B, Madigan JE, Dumler JS, Kurtti TJ, et al. Direct cultivation of the Causative Agent of Human granulocytic ehrlichiosis. N Engl J Med. 1996;334(4):209–15.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic Local Alignment Search Tool. Vol. 215, J Mol Biol. 1990.
Šimánek V, Pecen L, Krátká Z, Fürst T, Řezáčková H, Topolčan O, Fajfrlík K, Sedláček D, Šín R, Pazdiora P, et al. Five commercial Immunoassays for SARS-CoV-2 antibody determination and their comparison and correlation with the Virus neutralization test. Diagnostics. 2021;11:593.
The authors are grateful for the great effort provided by hunters and veterinary colleagues of the provinces of Belluno, Udine and Varese (Italy) that played an essential role in the support of field activities and sample collection.
Open access funding provided by Università degli Studi di Padova. This work was financially supported by the Integrated Budget for Interdepartmental Research BIRD223833/22 of the Department of Animal Medicine, Production and Health (MAPS). “The funding body played no role in the design of the study and collection, analysis, interpretation of data, and in writing the manuscript”.
The authors declare that they have no competing interests.
Ethics approval and consent to participate
Animals included in the study were hunted during licenced hunting seasons and not culled specifically for the purpose of this study. Therefore, ethical committee approval was not required since animals were hunted by authorized hunting associations according to the currently effective law in Italy regarding wildlife protection and hunting procedures (Law n. 157/92 “Norma per la protezione della fauna omeoterma e prelievo venatorio”).
Consent for publication
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Grassi, L., Drigo , M., Zelená, H. et al. Wild ungulates as sentinels of flaviviruses and tick-borne zoonotic pathogen circulation: an Italian perspective. BMC Vet Res 19, 155 (2023). https://doi.org/10.1186/s12917-023-03717-x
- Vector-borne zoonotic pathogens
- Wild ungulates
- Ixodes ricinus
- Molecular biology
- Virus neutralization test