- Research
- Open access
- Published:
Predicting habitat suitability for Ixodes ricinus and Ixodes persulcatus ticks in Finland
Parasites & Vectors volume 15, Article number: 310 (2022)
Abstract
Background
Ticks are responsible for transmitting several notable pathogens worldwide. Finland lies in a zone where two human-biting tick species co-occur: Ixodes ricinus and Ixodes persulcatus. Tick densities have increased in boreal regions worldwide during past decades, and tick-borne pathogens have been identified as one of the major threats to public health in the face of climate change.
Methods
We used species distribution modelling techniques to predict the distributions of I. ricinus and I. persulcatus, using aggregated historical data from 2014 to 2020 and new tick occurrence data from 2021. By aiming to fill the gaps in tick occurrence data, we created a new sampling strategy across Finland. We also screened for tick-borne encephalitis virus (TBEV) and Borrelia from the newly collected ticks. Climate, land use and vegetation data, and population densities of the tick hosts were used in various combinations on four data sets to estimate tick species’ distributions across mainland Finland with a 1-km resolution.
Results
In the 2021 survey, 89 new locations were sampled of which 25 new presences and 63 absences were found for I. ricinus and one new presence and 88 absences for I. persulcatus. A total of 502 ticks were collected and analysed; no ticks were positive for TBEV, while 56 (47%) of the 120 pools, including adult, nymph, and larva pools, were positive for Borrelia (minimum infection rate 11.2%, respectively). Our prediction results demonstrate that two combined predictor data sets based on ensemble mean models yielded the highest predictive accuracy for both I. ricinus (AUC = 0.91, 0.94) and I. persulcatus (AUC = 0.93, 0.96). The suitable habitats for I. ricinus were determined by higher relative humidity, air temperature, precipitation sum, and middle-infrared reflectance levels and higher densities of white-tailed deer, European hare, and red fox. For I. persulcatus, locations with greater precipitation and air temperature and higher white-tailed deer, roe deer, and mountain hare densities were associated with higher occurrence probabilities. Suitable habitats for I. ricinus ranged from southern Finland up to Central Ostrobothnia and North Karelia, excluding areas in Ostrobothnia and Pirkanmaa. For I. persulcatus, suitable areas were located along the western coast from Ostrobothnia to southern Lapland, in North Karelia, North Savo, Kainuu, and areas in Pirkanmaa and Päijät-Häme.
Conclusions
This is the first study conducted in Finland that estimates potential tick species distributions using environmental and host data. Our results can be utilized in vector control strategies, as supporting material in recommendations issued by public health authorities, and as predictor data for modelling the risk for tick-borne diseases.
Background
In temperate and boreal forests, ticks are among the most important arthropods transmitting viruses and other pathogens causing diseases in humans. In Northern Europe, among the most prevalent tick-borne pathogens are Borrelia burgdorferi (s.l.) bacteria causing Lyme borreliosis (LB) and tick-borne encephalitis virus (TBEV) causing tick-borne encephalitis (TBE). The abundances and distributions of vectors and their hosts are mainly restricted by their habitat suitability, determined largely by climatic and environmental conditions, similarly to other arthropods. Finland lies in a zone with overlapping geographical distributions of two tick species, the castor bean tick Ixodes ricinus Linnaeus 1758 and the taiga tick Ixodes persulcatus Schulze 1930. Tick abundances have increased in Finland, which has been mirrored by a rising incidence of tick-borne diseases (TBDs) in the country during the past decades [1,2,3,4]. The distribution of I. ricinus extends throughout Europe from Sweden to North Africa and from Ireland to the Urals [5, 6], and the distribution of I. persulcatus extends from Fennoscandia to Japan. During the last 2 decades, the taiga tick, I. persulcatus, has spread to new areas in Northern Europe: from eastern Finland to the north-western coast of Finland and all the way to eastern Sweden [7,8,9]. Currently, ticks occur throughout Finland, excluding most of Lapland [9]. Ixodes ricinus is predominant in southern Finland while I. persulcatus prevails in the eastern regions of Finland and on the north-western coast up to southern Lapland, although I. persulcatus may be spreading southwards, with recent records from the capital region along the south coast [10]. Central Finland is considered a sympatric area where both species occur.
Knowledge of tick seasonalities, i.e. the time periods when ticks are active, is important because of its public health relevance, but also for understanding how environmental and other factors influence their distributions. Ixodes ricinus is active from May to September in Finland, with two activity peaks. Adult I. ricinus activity peaks in July–August and nymph activity peaks in September [9, 11], while I. ricinus nymphs and adult females have been found to peak in May–June [12] and in August–September [12, 13]. The seasonal activity of I. persulcatus adults in Finland is unimodal, as they are active in May–July, with the highest activity peak in May [9, 11, 14]. Tick reproduction and abundance and pathogen transmission are dependent on the presence of a suitable vertebrate host. Adult and nymphal I. ricinus and I. persulcatus often feed on small to medium-sized animals, such as rodents, birds, foxes, and hares, but also on large-sized animals including deer, moose, and humans [15,16,17]. In contrast to I. ricinus nymphs, I. persulcatus nymphs rarely feed on larger hosts, including humans [14]. Larval I. ricinus and I. persulcatus feed mainly on terrestrial vertebrates, such as voles or birds, as they are numerous and move close to the ground in the vicinity of tick larvae [18, 19]. Tick densities vary largely depending on geographical location. Generally, 100–1000/100 m2 is considered high density for I. ricinus nymphs in Central Europe, while 10–13/100 m2 and 10/100 m2 already indicate high densities for I. ricinus [2, 20] and I. persulcatus nymphs or adults in Finland [14], respectively. High density numbers in Finland coincide with those observed for I. persulcatus in northern Sweden [21].
Weather, climate change, and the environment influence habitat suitability, vector activity, and the rate of vector development [22, 23]. Ticks need sufficiently high temperatures to complete their development [24]. Ixodes ricinus and I. persulcatus differ in terms of their temperature-related activity, with I. persulcatus being active at lower temperatures [25,26,27]. Generally, warm temperatures and increased rainfall positively affect tick densities [22], but extremely high temperatures combined with decreased rainfall may reduce tick populations [28]. Water vapour is the main source of moisture for active unfed ticks, and it is only absorbed at sufficiently high relative air humidity [24]. Ixodes ricinus requires high relative humidity (> 80%) to survive during off-host periods [29]. In the northern distribution limit of the species, relative humidity was found to positively correlate with the abundances of I. ricinus nymphs and adults [2]. Saturation deficit was also found to positively influence the abundance of questing I. ricinus larvae and to positively influence adult abundance until the optimal value (3.16 mmHg) [12]. Despite less research conducted on I. persulcatus than on I. ricinus in their northern distribution limits, findings suggest that the climatic requirements of I. persulcatus differ from those of I. ricinus. Ixodes persulcatus occurs in areas where growing season length varies between 140 and 150 days, while I. ricinus prefer areas with a growing season of > 180 days [8]. Furthermore, I. persulcatus occurs in areas with less precipitation and a lower humidity index and temperature sum than the northern distribution limits of I. ricinus [8, 24]. Snow cover influences tick development, acting as an insulator that protects the nymphal ticks and larvae from freezing [30].
Species distribution modelling (SDM) methods can be used to determine the habitat suitability areas of the two species and to indicate the most influential factors affecting their distributions [31]. SDM techniques have widely been used to study I. ricinus occurrences [32,33,34,35,36,37] in Europe. In addition to studies focused on modelling I. ricinus occurrences, multiple studies have modelled I. ricinus densities [38,39,40,41], which requires more comprehensive and continuous tick sampling than needed for presence-absence modelling. No earlier SDM studies have estimated I. persulcatus distributions in Europe or Asia, although a few studies have estimated their probable occurrences [42, 43]. Earlier research in Finland on tick distributions shows occurrence point patterns [9] and I. ricinus [11, 12, 20] and I. persulcatus [11, 14] abundances at smaller scales, but SDM methods have not been conducted previously. Host data are difficult to obtain and are thus missing from the majority of SDM studies examining ticks. However, several statistical analyses on I. ricinus distributions have included tick host data [44,45,46,47].
In this study, we aim to (1) analyse the ticks collected in 2021 for species and possible pathogens: TBEV and Borrelia, (2) estimate the potential distributions of I. ricinus and I. persulcatus in mainland Finland and assess the model performance using various variable compositions, and (3) identify the most influential factors driving the spatial patterns. The study aim is to predict tick species occurrences instead of tick densities, and thus habitat suitability in the Results refers to the probability of species presence instead of the estimated abundance. Our study is based on data obtained from historical tick data sets from 2014 to 2020 and on new tick collections from 2021.
Methods
Study area
Finland (59°50′N, 20°38′E, 70°09′N, 31°30′E) is located in Northern Europe between Sweden and Russia (Fig. 1). This study covers Finland, excluding the Åland Islands and most parts of Lapland. Only the south-western municipalities of Lapland were included, Simo, Keminmaa, Tornio, and Kemi, because of the presence of both tick species in the region [9].
Tick data
Historical tick data
Historical tick data were collected by the Universities of Turku, Jyväskylä, and Helsinki. The majority of the data consisted of national crowdsourcing data, where citizens collected ticks across Finland from April to November, 2015, and sent them to the University of Turku for identification and further analysis [9]. The coordinates for the crowdsourcing data were manually gathered and recorded by biologists from the University of Turku, based on information provided in the letters used to send the ticks, along with their own assessment of suitable areas for ticks (based on topographic and satellite imagery). As for pinpointing the location where each tick got on its vessel, most ticks were sent in by dog owners after being removed from the dogs. These dog owners reported consistently walking approximately the same routes with their pets, and thus we would expect the ticks to have gotten on the dogs within a 500-m diameter of the mentioned walking areas. These factors naturally lead to uncertainty in the correctness of the coordinate points and, more notably, to uncertainty regarding the accuracy to which the coordinates depict the exact location for acquiring the tick. However, we would expect the coordinate points and surrounding environments to represent the tick acquisition area at a 1-km resolution in most cases. In addition, tick collections at smaller scales were conducted by the Universities of Helsinki, Turku, and Jyväskylä from 2014 to 2020. Aggregated historical data consisted of 4152 presences for I. ricinus and 986 presences for I. persulcatus.
Sampling strategy for tick collections in 2021
Historical I. ricinus data covered the entire study area well, but I. persulcatus occurrence data were spatially clustered. We therefore targeted the activity season of I. persulcatus from May to June in 2021 to fill the gaps in environmental space covered by the data. We created a GIS-based sampling strategy using both ESRI ArcGIS (version 10.3.1) (ESRI, Redlands, CA, USA) and VECMAP® software [48] with the following criteria (Additional file 1: Fig. S1a). We created "regulated areas" as buffer zones around existing I. persulcatus occurrences within a 5-km radius to exclude them from the sampling strategy. Then, we created 500-m buffer zones around the roads to ensure accessibility. We used the CORINE Land cover CLC2018 data set available for Finland [49] to extract land cover classes known to be suitable habitats for ticks: a forest category including deciduous, coniferous, and mixed forests and a meadow category including sparsely wooded areas. We used subdivisions of landscape provinces to distinguish areas already well covered with I. persulcatus occurrences from four larger areas where further collections were needed. For each of the four areas, a random sample of 25 collections points (a total of 100 locations) was created in VECMAP® based on the relative shares of the two land cover categories in the area (approximately 86% forests and 14% meadows). As ecotones are locations with most tick abundance [50], each created collection point was moved to the closest ecotone using a Google Maps Satellite map.
Tick collections and processing
Additional sampling of 100 locations was carried out from the beginning of May to the end of June 2021, which is considered the activity season of I. persulcatus (peak in May) in Finland [9, 11, 14]. Ticks were collected by slowly dragging a 1.0 × 1.5-m or a 1.0 × 1.0-m cotton cloth over 10-m sections, with a total of 400-m dragging sessions in each locality. This design was defined together with Finnish tick specialists. If ticks were not detected after 400 m of dragging, the site was considered an absence location. Sampling was conducted during the day or night but not on rainy days. Sampling date or time may affect the sampling results. All ticks attached to the cloth were collected using tweezers. The larvae, nymphs, and adults were separated, placed into 15-ml Greiner tubes with a grass stalk inserted in each tube, and transported to the Department of Virology at the University of Helsinki. Ticks were stored alive at 4 °C until homogenization. Ticks were pooled according to species, development stage, and collection site (adults 1; nymphs 1–5; larvae 1–50). RNA was extracted using a QIAamp Viral RNA Mini Kit (Qiagen) and DNA using a DNeasy Blood and Tissue DNA extraction kit (Qiagen, Hilden, Germany) according to the manufacturers’ instructions. Tick species was confirmed by molecular identification, using a species-specific duplex TaqMan real-time PCR assay, as previously described [1]. TBEV was screened using a quantitative real-time PCR assay [51], with the following modifications; assays were carried out in 20 μl of reaction volume, including 5 μl Taqman fast virus 1-step mastermix (Thermo Fischer Scientific), 0.4 μl forward primer, 0.6 μl reverse primer, 0.16 μl probe, 8.84 μl H2O, and 5 μl RNA. The thermal cycling profile was 50 °C for 5 min, followed by 20 s at 95 °C, prior to amplification (95 °C 15 s, 60 °C 30 s, 40 cycles). For Borrelia, we used a quantitative real-time PCR assay, as previously described [9]. All PCR reactions were performed on a T100™ Thermal Cycler (Bio-Rad, Germany).
Data pre-processing
Spatial autocorrelation (SA) refers to data or residuals correlated with themselves rather than being independent [52], and it may inflate the effective sample size and bias parameter estimates. To reduce SA, tick occurrence data were thinned using R package Wallace [53], which uses the spThin approach [54]. For I. ricinus presence-absence data, we used 10 km of spatial thinning for I. ricinus and 3 km of spatial thinning for I. persulcatus, which were considered I. ricinus absences with true absences (N = 88) obtained from the summer sampling in 2021. For I. persulcatus presence-absence data, 3 km of spatial thinning for I. persulcatus was used to obtain the highest possible number of presences and 5 km of spatial thinning for I. ricinus. These data were aggregated with data of true absences described above. After data thinning, I. ricinus data consisted of 622 presences and 637 absences, and I. persulcatus data consisted of 509 presences and 1289 absences (Fig. 2a, b).
Geospatial data
Environmental, climate, and other predictors were selected based on factors known to influence the I. ricinus and I. persulcatus distributions. Environmental data for Finland were obtained from various sources and included data derived from satellite imagery, GIS layers, or interpolated data. Details of the predictor data are provided in Table 1. The initial data sets for both species included 25 predictors before running a multicollinearity analysis. Final predictors in each variable combination are seen in Fig. 3 and Additional file 5: Table S1. For monthly climate data (relative air humidity, air temperature, and precipitation), we calculated the mean for the activity seasons for both species: May–September for I. ricinus and April–June for I. persulcatus.
Host data
Host abundance estimates were based on snow track counts from the wildlife triangle census, which is a national monitoring scheme coordinated by the Natural Resources Institute Finland (LUKE). Abundance was approximated using a snow track index (I = tracks/10 km/day) for each species and triangle in an area, i.e., the number of snow tracks observed per 10 km of counted transect, further divided by the number of days during which snow tracks have accumulated (days since the last snowfall or since a round of “pre-mapping”). For each species occurrence, we calculated the average annual snow track index within a 50-km radius, and the annual indexes for that area were further averaged over the study period (2014–2021). The layer of species-specific distributions was formed by repeating this procedure over a 1 × 1-km grid (Table 1). The data consist of calculated average abundance indexes for European hare (Lepus europaeus), mountain hare (Lepus timidus), moose (Alces alces), red fox (Vulpes vulpes), roe deer (Capreolus capreolus), and white-tailed deer (Odocoileus virginianus).
Data preparation and data analysis
We used four variable compositions for the SDM predictions: environment only, host only, combined data on environmental and host variables, and combined data on environmental and host variables with suitability data of either tick species. In the fourth variable composition, the suitability data of either tick species were produced based on environmental and host data to reveal any potential effects of I. ricinus occurrence to I. persulcatus presence and vice versa. We decided to include this predictor because of field observations in which the one species has not been found in the vicinity of the other species. We used the biomod2 platform in R (version 3.4.6) [61, 62] to create SDMs to identify areas with suitable habitat conditions for I. ricinus and I. persulcatus. All geospatial data sets, including environmental and other data, were processed in ESRI ArcGIS (version 10.3.1) (ESRI, Redlands, CA, USA) and were set to the same spatial extent, geographic coordinate system (EUREF FIN TM35FIN, epsg:3067), and resolution (1 × 1 km). Multicollinearity of the variables was investigated using Variance Inflation Factors (VIFs), as implemented in R package usdm [63, 64]. Predictor VIFs were calculated and correlated variables were excluded in a stepwise procedure using a commonly applied threshold value of 10 [65, 66]. The following eight predictive modelling techniques were employed: generalized linear models (GLM) [67], generalized additive models (GAM) [68], classification tree analysis (CTA) [69], artificial neural networks (ANN) [70], multivariate adaptive regression splines (MARS) [71], generalized boosting models (GBM) [72], random forest (RF) [73], and maximum entropy (MAXENT) [74]. Flexible discriminant analysis (FDA) and surface range envelope (SRE) were excluded due to generally poor predictive performance [75, 76]. Models were mostly run using the default settings of biomod2, with the following exception: we used the GAM function in the mgcv package, with k = 3 as the basis dimension for the thin plate smoothing terms [77]. We used a cross-validation technique where we split the thinned data set into two subsets, one to calibrate the models (70%) and another to evaluate the models (30%) [31]. We repeated the calibration and evaluation sets 50 times for each model (400 model evaluation runs in total) [31]. To reduce the uncertainty related to the choice of a single modelling technique, we built ensemble predictions using the ensemble mean method that averages predictions across the best-performing individual models with the selected threshold (0.7 < AUC < 1.0) [62].
The accuracy assessment
Sensitivity (the proportion of correctly predicted presences) and specificity (the proportion of correctly predicted absences) were calculated to quantify omission errors [78]. Sensitivity is calculated by dividing true presences (TP) by the sum of TP and false absences (FA). Similarly, specificity is calculated by dividing true absences (TA) by the sum of TA and false positives (FP). The area under the curve (AUC), used here to assess model performance, is the measure of a model’s ability to distinguish between these presence and absence classes. The AUC scores range from 0 to 1, with 0.5 being the threshold for predictions better than random [78] and > 0.7 being an acceptable threshold for predictions [79]. The variable importance of predictors, based on decreasing accuracy, was extracted from the biomod2 output [56]. Each importance value was normalized by dividing the value by the total sum of importance values to compare the most powerful variables. Partial dependency plots were generated to show the predictors’ estimated effects on tick distributions. Prediction uncertainty was assessed by the coefficient of variation (CV) of predictions, with a high CV value indicating high uncertainty in the predicted distributions [56].
Results
Additional tick collections in 2021
We identified the tick species from the additional data collections in 2021, and the ticks were screened for TBEV and Borrelia pathogens. We sampled 89 new locations, which included 25 new presences and 63 absences for I. ricinus, and only one presence and 88 absences for I. persulcatus. Collections were not conducted in 11 locations because of the tight collection schedule. A total of 502 ticks (35 adults, 375 nymphs, and 92 larvae) were collected from 26 locations and were pooled into 120 pools according to locality, species and development stage. Five hundred one ticks were identified as I. ricinus and only one nymph as I. persulcatus. We found no tick pools positive for TBEV. In contrast, 56 (47%) of the 120 tick pools were positive for B. burgdorferi (s.l.), consisting of 220 ticks, all identified as I. ricinus. The minimum infection rate (MIR) calculated for pooled ticks was 11.2%. Results from the additional collections in 2021 are shown in Additional file 1: Fig. S1b.
Predictive model accuracies
We used four variable compositions for the predictions: environment only, host only, combined data on environmental and host variables, and combined data on environmental and host variables with suitability data of either tick species. The suitability data of I. ricinus and I. persulcatus for the previous variable composition were first produced using combined data of environmental and host variables.
The means and ranges of the predictive performances over 400 model runs by eight individual modelling techniques and four variable compositions are shown in Additional files 2 and 3: Figs. S2–S3. In all variable compositions, the GBM and RF models gained the highest model performances. The mean ensemble models were built over several modelling techniques, and the number of times they contributed to the final ensemble are seen in Additional file 4: Table S1. The mean ensembles estimating the habitat suitabilities yielded AUC values in the range of 0.90–0.94 (Table 2) for I. ricinus and 0.91–0.96 for I. persulcatus. The highest predictive performances from the four variable compositions were obtained with the combined data sets for both species in both the individual models and mean ensembles, and the lowest predictive performances were obtained with environmental data only (Additional files 2 and 3: Figs. S2–S3; Table 2).
Sensitivity and specificity (by AUC) for estimating the distributions of I. ricinus and I. persulcatus varied based on the variable compositions. To estimate I. ricinus distributions, the mean ensemble model better identified unsuitable environments in all variable compositions (sensitivity = 75.2–78.3%, specificity = 84.2–89.5%), excluding the variable composition with environmental, host, and habitat suitability data of I. persulcatus (sensitivity = 92.2%, specificity = 79.2%). In contrast, when estimating I. persulcatus distributions, we found a better ability to identify suitable environments (sensitivity 89.3–92.7%) than unsuitable environments (specificity 78.1–83.8%). The highest sensitivity and specificity rates were obtained from the variable composition that contained all the variables for I. persulcatus (sensitivity = 91.9%, specificity = 83.8%).
Predictor contributions to tick species distributions
Variable importance, which indicates the influence of the variable to the mean ensemble model, is referred to here as the relative contribution of the predictor (%). The higher the value, the more influence the variable has on the model. The relative contributions varied between the species and the variable compositions. Here, we present the relative contributions of the predictors used in the combined data sets (Fig. 3a, b), which gained the highest model performances. The relative contributions of the environmental only and host only data sets are shown in Additional file 5: Fig. S4a–b.
The highest relative contributions for I. ricinus predictions based on combined host and environmental data were obtained from the mean temperature of the I. ricinus activity season (42%), precipitation sum (16%), and red fox (15%), white-tailed deer (6%), and roe deer (4%) densities. For I. persulcatus, red fox density (24%), precipitation sum (23%), the densities of mountain hare (13%), white-tailed deer (10%), and European hare (9%), the mean precipitation for the I. persulcatus (5%) activity season, and roe deer density (4%) were the variables with the highest relative contributions. Predictions based on host, environmental, and habitat suitability data of the other species (Fig. 3b) indicated slightly different relative contributions. When predicting I. ricinus distributions, the highest relative contributions were obtained from the habitat suitability for I. persulcatus (74%), the mean temperature of the I. ricinus activity season (6%), and the red fox (5%), mountain hare (4%), and white-tailed deer (3%) densities. In contrast, the habitat suitability for I. ricinus (80%), night land surface temperature (6%), elevation (3%), and roe deer density (3%) were the most important predictors for I. persulcatus.
We analysed partial dependency plots for I. ricinus and I. persulcatus based on the combined data set of environmental and host data (Fig. 4) because of the high influence of other species in the combined data set of environmental, host, and suitability data for the other species (see Fig. 3b). Partial dependency plots based on environmental only data, host only data, and a combined data set of environmental, host, and suitability data for the other species are shown in Additional file 6: Figs. S5–S7. Locations with higher relative humidity (> 75%), mean air temperature during the activity season (> 13 °C), precipitation sum (> 300 mm), middle-infrared reflectance (MIR, > 0.25), and higher white-tailed deer (> 30 individuals/km2), European hare (10–40/km2), and red fox (> 10/km2) densities were associated with higher habitat suitability for I. ricinus (Fig. 4a). However, higher mountain hare density (> 10/km2) was negatively associated with I. ricinus suitability. For I. persulcatus, higher mean precipitation (> 40 mm) and white-tailed deer (> 5/km2), roe deer (> 5/km2), and mountain hare (> 10/km2) densities were associated with high suitabilities for I. persulcatus (Fig. 4b). However, an excessively high precipitation sum (> 320 mm), higher mean air temperature during the activity season (> 9 °C), and day land surface temperature (DLST, > 10 °C) at the locations began negatively influencing the suitability for I. persulcatus. Furthermore, higher moose (> 8/km2) and red fox (> 10/km2) densities indicated lower suitabilities for I. persulcatus in the locations.
Habitat suitabilities for I. ricinus and I. persulcatus based on variable compositions
The habitat suitability maps for I. ricinus and I. persulcatus in four variable compositions are shown in Figs. 5 and 6. We focus on analysing the results from the combined data set of environmental and host data for I. ricinus (Fig. 5c) and I. persulcatus (Fig. 6c) based on the high predictive performance and the influence of the predictors contributing to the final ensemble model. In this study, low suitability for species presence is interpreted as 0–30%, moderate suitability as 31–60%, and high suitability as 61–100%. Figure 5c shows that the areas with moderate to high suitability for I. ricinus were located southwards from Central Ostrobothnia, with the following exceptions: narrow areas located in southern Pirkanmaa and southern coast of Ostrobothnia. Northern regions of North Savo, North Karelia, North Ostrobothnia, and Kainuu were estimated to have a low to moderate suitability for I. ricinus. Areas with moderate to high suitability for I. persulcatus were located mainly northwards from Ostrobothnia up to southern Lapland, including areas along the western coast, and eastern Finland (Fig. 6c). The moderate to high suitability areas for I. persulcatus in southern Finland were located across Pirkanmaa and in narrow areas of Kanta-Häme, Päijät-Häme, South Karelia, South Savo, and Uusimaa. Other areas in southern Finland, the northern part of Kainuu, and the eastern part of North Ostrobothnia were estimated to have a low to moderate probability for I. persulcatus.
Discussion
Study validity
No earlier SDM studies have estimated tick species distributions in Finland. However, there is a greater need for studying tick distributions when tick abundances and TBD incidences are increasing [80]. We used the mean ensemble model approach over several modelling methods, which are generally found to yield more accurate estimates than single-model estimates [81, 82] and are widely used to estimate the potential distributions of vectors and TBDs [83,84,85]. We assessed model performance using variable compositions and predictions uncertainties by using the coefficient of variation approach. Prediction uncertainties are often not considered in SDM studies, although the uncertainty assessment is highly recommended in SDM literature [86, 87]. This study obtained results at a 1-km resolution, which provides higher accuracy than generally in the SDM studies. The analyses utilized host data, which are often unavailable and thus are not included in a majority of tick SDM studies. Furthermore, species data were not aggregated from big data repositories, such as the Global Biodiversity Information Facility (GBIF), which are widely used but may affect results, e.g. due to spatial bias [88]. The results can be utilized as predictor data in future studies estimating the TBD risk in Finland.
Predictive performance
Predictive performance of the mean ensembles over several modelling methods yielded model performances with minimum AUC values of 0.90 for both I. ricinus (Fig. 5) and I. persulcatus (Fig. 6). The mean ensemble models showed improvements of the AUC values compared with individual models (Additional files 2 and 3: Figs. S2–S3). Model performances can be classified into different categories: AUC values of 0.7–0.8 are considered fair model performance, 0.8–0.9 are good model performance, and > 0.9 are excellent model performance [89]. Mean ensembles yielded good to excellent model performances, while individual models resulted in fair to good model performances. Of the four variable compositions, ensemble models with combined data sets yielded the highest predictive accuracy for both I. ricinus (Fig. 5) and I. persulcatus (Fig. 6). Ecological and vegetation data, together with climatic data, are recommended for inclusion in the predictor data set when modelling species and disease distributions [86]. The ensemble model better identified unsuitable environments (i.e. specificity) in I. ricinus predictions in contrast to I. persulcatus predictions, in which suitable environments were better identified (i.e. sensitivity, Table 2). The high specificity rate for I. ricinus ensured that the proportion of true absences predicted as presences was minimized, and in contrast, the high sensitivity rate for I. persulcatus indicated that the proportion of true presences predicted as absences was minimized.
The uncertainty assessment
Other influential factors are not included in this study such as microclimate and vertebrate density data. Vertebrate data, including e.g. rodent, bird, and raccoon dog densities, were not available in similar accuracy or spatial scale. Bank vole abundance in particular has previously been found to closely relate to questing I. ricinus abundances in Finland [12]. Along with rodents, birds are important hosts for Ixodes ticks [90]. However, bird density data covering the whole country are not available. Furthermore, some uncertainty arose from tick absences. The collections in 2021 resulted in 88 true absences for I. persulcatus. However, to obtain enough absences for I. persulcatus, we needed to use randomly selected points from the I. ricinus presences, obtained with the Wallace package. Only one thinned data set was obtained with Wallace, so we were unable to assess the occurrence-related uncertainty. True absences are difficult to obtain, and we are aware that the absences used in this study may not be true absences but may rather depend on e.g. collection time or weather conditions. As the ensemble model was better at identifying suitable than unsuitable environments for I. persulcatus, further studies are needed to bring new insights to whether this is an indication that the species is still spreading or whether this is due to modelling-specific reasons. Using presence-absence data instead of tick abundance data loses information on the relative suitability of habitats when all presences are treated as equal, regardless of the abundance of the individuals that the habitat supports. We excluded the Åland Islands and most of Lapland from the study due to the following reasons. The Åland Islands have abundant tick populations, and we did not have enough tick samples to reliably estimate their distributions. Only occasional tick observations have been made from Lapland, excluding the south-eastern municipalities of Kemi, Tornio, Keminmaa, and Simo.
Hereafter, we focus on analysing the results based on the combined data set of environmental and host data due to high predictive performance (Figs. 5, 6), variable influences in the ensembles (Fig. 3a, b), and the recommendation to use ecological and vegetation factors with climatic factors in SDM [86]. The model-driven uncertainty of the predictions was presented using the coefficient of variation to assess the uncertainty of the suitability maps (Figs. 7, 8). In this study, low uncertainty is interpreted as 0–30%, moderate uncertainty as 31–60%, and high uncertainty as 61–100%. The areas with highest uncertainty for I. ricinus suitability were mainly located in North Ostrobothnia, southern Lapland, and Kainuu, and in narrow areas of Ostrobothnia, Pirkanmaa, and southern Uusimaa (Fig. 7c). High uncertainty in the northern and eastern parts of the study area may be associated with climatic factors of the areas characterized by shorter growing seasons, lower relative air humidity rates, and lower mean monthly air temperatures than in southern Finland. The lowest uncertainties were obtained from southern, eastern, central, and south-western Finland, indicating more reliable prediction results. For I. persulcatus, the areas with high uncertainty were somewhat larger in geographical extent than in the I. ricinus predictions (Fig. 8c). The highest prediction uncertainties were located across southern Finland, South Ostrobothnia, Central Finland, South Savo, Northern Kainuu, and north-eastern parts of North Ostrobothnia. In contrast, the areas with highest prediction confidence for I. persulcatus were located in coastal Finland northwards from Ostrobothnia up to southern Lapland, in Pirkanmaa, southern Kainuu, North Savo, and North Karelia.
Influential factors
Consistent with previous research, the environmental, host, and climatic variables were important determinants for I. ricinus and I. persulcatus occurrence. Our study suggests that climatic factors, such as higher relative humidity, mean air temperature, and precipitation sum during the growing season, were associated with higher I. ricinus occurrence. Higher air temperatures [36, 37] and precipitation, especially in spring [91], have previously been found to positively influence I. ricinus presence. Based on our study results, higher red fox, white-tailed deer, and European hare densities were associated with higher habitat suitabilities for I. ricinus, which is mainly consistent with earlier findings (Fig. 3a) [92,93,94,95]. Although previous studies have found red foxes to be suitable hosts for ticks in Europe [96, 97], they were not considered suitable hosts in a recent study from Norway [93]. Mountain hare density had a negative effect on I. ricinus presence. This is obvious, as mountain hare densities have significantly declined in southern and western Finland during the past 30 years. However, as host densities may be directly and indirectly affected by climate, it is difficult to separate causal and confounding factors from one another [98, 99]. Consistent with previous research [36], our study suggests that high middle-infrared reflectance (MIR) levels positively influenced I. ricinus occurrence.
Similar to I. ricinus, I. persulcatus suitabilities were higher in locations with higher mean precipitation and air temperature during the activity season. However, when precipitation, mean air temperature, and daily land surface temperature (DLST) increased to a particular point, the suitability for I. persulcatus began decreasing. This finding may demonstrate the characteristics of I. persulcatus to prosper in slightly drier and colder habitats than I. ricinus [24]. However, the expansion of I. persulcatus' range has been found to correlate with the increase in mean annual air temperatures, which determine compatible temperature conditions for I. persulcatus establishment in new territories [100, 101]. As warmer winters and hotter summers are estimated to change the dynamics and pattern of seasonal tick activity [102], it remains to be seen whether I. persulcatus adapts to warmer and wetter habitat conditions. Some adaptation has already occurred with recent range expansion to southern parts of Finland [10]. In contrast to I. ricinus, higher red fox density negatively affected I. persulcatus. This is probably related to spatial variation in red fox density across the country; highest densities occur in southern parts decreasing steadily northward, to I. persulcatus-dominated area. Our study also suggests that higher densities of white-tailed deer, roe deer, and mountain hare were associated with higher habitat suitabilities for I. persulcatus, which is in line with previous findings [14, 99, 103]. Notably, white-tailed deer and roe deer populations have rapidly increased in southern Finland during the past few years, especially in the south-west [104], which may potentially affect not only the rise in I. ricinus abundances but also the southward spread of I. persulcatus.
Habitat suitabilities for I. ricinus and I. persulcatus in Finland
Ixodes ricinus and I. persulcatus have their own environmental and other limits for surviving and reproducing, which restrict their geographical distributions. Based on our study results, moderate to high suitability areas for I. ricinus occurred throughout southern and Central Finland up to Central Ostrobothnia (64°N), excluding narrow areas in Ostrobothnia and Pirkanmaa (Fig. 6c). In neighbouring Sweden, only areas southwards from the capital region (60°N) were predicted to have abundant I. ricinus populations [36], despite I. ricinus having been found up to 66°N [99]. Based on recent I. ricinus studies from Russian Karelia, the species was absent already north of 63°N [105]. The narrow areas in Ostrobothnia and Pirkanmaa, considered sympatric areas, were estimated to have low suitability for I. ricinus, which may, partly, be explained by the moderate uncertainty in the prediction (Fig. 8c). The dominance area of I. persulcatus is known to be more northerly than that of I. ricinus. We note that the suitability for I. persulcatus was highest northwards from Ostrobothnia along the northern coast up to southern Lapland, Kainuu, North Savo, and North Karelia (Fig. 7c). Southern Finland, excluding areas in Pirkanmaa, western Päijät-Häme, northern Kanta-Häme, and southern Uusimaa, were estimated to have low habitat suitability for I. persulcatus. Also, there is a moderate to high uncertainty in the predictions for I. persulcatus across southern Finland (Fig. 8c), which may indicate that more areas with high suitabilities for the taiga tick may exist. Notably, I. persulcatus have only been found up to 63°N in Russian Karelia [43, 106]. In Sweden, the species was first observed in 2015 close to the Finnish border at ≈ 66°N [8].
Conclusions
An increased risk of vector-borne pathogens and the spread of invasive and naturally spreading species due to changing weather patterns adds to the needs and requirements for increased research and concrete actions. In our study, habitat suitability areas for I. ricinus and I. persulcatus were identified for the first time in Finland. From additional tick collections in 2021, 25 new presences and 63 absences were found for I. ricinus, and 1 presence and 88 absences for I. persulcatus. A total of 502 ticks were analysed for pathogens, with no ticks positive for TBEV and ≈ 47% of tick pools positive for Borrelia burgdorferi (s.l.). High suitability areas for I. ricinus occurred throughout southern and Central Finland up to Central Ostrobothnia, excluding narrow areas in Ostrobothnia and Pirkanmaa. For I. persulcatus, the regions northwards from Ostrobothnia along the northern coast up to southern Lapland, Kainuu, North Savo, North Karelia, and areas in Pirkanmaa and Päijät-Häme were estimated to be suitable areas. Based on the predictions, locations with higher air temperature, relative humidity, precipitation sum, and MIR and higher densities of white-tailed deer, European hare, and red fox were suitable for I. ricinus. For I. persulcatus, higher mean precipitation and higher densities of white-tailed deer, roe deer, and mountain hare indicated a higher occurrence probability. The data produced in this study have implications for improving knowledge on disease prevention and for assisting authorities in decision-making concerning vector control strategies. Our results can be used as predictor data to estimate the risk for TBDs in Finland. In future studies, our aim is to focus on studying tick distributions at various spatial scales: in microhabitats and at larger scales covering all of Scandinavia. Our aims are to use in situ measurements to achieve more accurate microclimate data, to conduct more extensive tick sampling to enable tick abundance modelling, to model the distribution of TBPs in ticks, and to forecast tick distributions in the future climate.
Availability of data and materials
Data supporting the conclusions of this article are included within the article and its additional files. Raw data are available from the corresponding author upon request.
Abbreviations
- AUC:
-
Area under the ROC curve
- LB:
-
Lyme borreliosis
- MIR:
-
Middle-infrared reflectance
- PCR:
-
Polymerase chain reaction
- SDM:
-
Species distribution modelling
- TBD:
-
Tick-borne disease
- TBE:
-
Tick-borne encephalitis
- TBEV:
-
Tick-borne encephalitis virus
- TBP:
-
Tick-borne pathogen
- VIF:
-
Variance inflation factor
References
Sormunen JJ, Penttinen R, Klemola T, Hänninen J, Vuorinen I, Laaksonen M, et al. Tick-borne bacterial pathogens in southwestern Finland. Parasit Vectors. 2016;1:168. https://doi.org/10.1186/s13071-016-1449-x.
Sormunen JJ, Klemola T, Vesterinen EJ, Vuorinen I, Hytönen J, Hänninen J, et al. Assessing the abundance, seasonal questing activity, and Borrelia and tick-borne encephalitis virus (TBEV) prevalence of Ixodes ricinus ticks in a Lyme borreliosis endemic area in Southwest Finland. Ticks Tick Borne Dis. 2016;7:208–15. https://doi.org/10.1016/j.ttbdis.2015.10.011.
Sormunen JJ, Klemola T, Hänninen J, Mäkelä S, Vuorinen I, Penttinen R, et al. The importance of study duration and spatial scale in pathogen detection—evidence from a tick-infested island. Emerg Microb Infect. 2018;7:1–11. https://doi.org/10.1038/s41426-018-0188-9.
NIDR. Finnish National Infectious Diseases Register. Finnish National Institute of Health and Welfare. 2021. https://thl.fi/fi/web/infektiotaudit-ja-rokotukset/seurantajarjestelmat-ja-rekisterit/tartuntatautirekisteri. Accessed 15 Jan 2022.
Estrada-Pena A. Forecasting habitat suitability for ticks and prevention of tick-borne diseases. Vet Parasitol. 2001;98:111–32. https://doi.org/10.1016/s0304-4017(01)00426-5.
Estrada-Pena A, Farkas R, Jaenson TG, Koenen F, Madder M, Pascucci I, et al. Association of environmental traits with the geographic ranges of ticks (Acari: Ixodidae) of medical and veterinary importance in the western Palearctic. A digital data set. Exp Appl Acarol. 2013;59:351–66. https://doi.org/10.1007/s10493-012-9600-7.
Jääskeläinen AE, Tonteri E, Sironen T, Pakarinen L, Vaheri A, Vapalahti O. European subtype tick-borne encephalitis virus in Ixodes persulcatus ticks. Emerg Infect Dis J. 2011;17:323. https://doi.org/10.3201/eid1702.101487.
Jaenson TG, Varv K, Frojdman I, Jaaskelainen A, Rundgren K, Versteirt V, et al. First evidence of established populations of the taiga tick Ixodes persulcatus (Acari: Ixodidae) in Sweden. Parasit Vectors. 2016;9:377. https://doi.org/10.1186/s13071-016-1658-3.
Laaksonen M, Sajanti E, Sormunen JJ, Penttinen R, Hanninen J, Ruohomaki K, et al. Crowdsourcing-based nationwide tick collection reveals the distribution of Ixodes ricinus and I. persulcatus and associated pathogens in Finland. Emerg Microbes Infect. 2017;6:e31. https://doi.org/10.1038/emi.2017.17.
Zakham F, Jaaskelainen AJ, Castren J, Sormunen JJ, Uusitalo R, Smura T, et al. Molecular detection and phylogenetic analysis of Borrelia miyamotoi strains from ticks collected in the capital region of Finland. Ticks Tick Borne Dis. 2021;12:101608. https://doi.org/10.1016/j.ttbdis.2020.101608.
Sormunen JJ, Andersson T, Aspi J, Bäck J, Cederberg T, Haavisto N, et al. Monitoring of ticks and tick-borne pathogens through a nationwide research station network in Finland. Ticks Tick Borne Dis. 2020;11:101449. https://doi.org/10.1016/j.ttbdis.2020.101449.
Cayol C, Koskela E, Mappes T, Siukkola A, Kallio ER. Temporal dynamics of the tick Ixodes ricinus in northern Europe: epidemiological implications. Parasit Vectors. 2017;10:166. https://doi.org/10.1186/s13071-017-2112-x.
Sormunen J, Klemola T, Vesterinen E, Vuorinen I, Hytönen J, Hänninen J, et al. Assessing the abundance, seasonal questing activity, and Borrelia and tick-borne encephalitis virus (TBEV) prevalence of Ixodes ricinus ticks in a Lyme borreliosis endemic area in Southwest Finland. Ticks Tick Borne Dis. 2015;7:208–15. https://doi.org/10.1016/j.ttbdis.2015.10.011.
Pakanen VM, Sormunen JJ, Sippola E, Blomqvist D, Kallio ER. Questing abundance of adult taiga ticks Ixodes persulcatus and their Borrelia prevalence at the north-western part of their distribution. Parasit Vectors. 2020;13:384. https://doi.org/10.1186/s13071-020-04259-z.
Gray JS. The ecology of Lyme borreliosis vector ticks. In: Bruin J, van der Geest LPS, Sabelis MW, editors. Ecology and evolution of the Acari: Proceedings of the 3rd Symposium of the European Association of Acarologists 1–5 July 1996, Amsterdam, The Netherlands. Dordrecht: Springer Netherlands; 1999. p. 269–76.
Uspensky I. The taiga tick Ixodes persulcatus (Acari: Ixodidae), the main vector of Borrelia burgdorferi sensu lato in Eurasia. 2016. p. 16.
Brugger K, Walter M, Chitimia-Dobler L, Dobler G, Rubel F. Seasonal cycles of the TBE and Lyme borreliosis vector Ixodes ricinus modelled with time-lagged and interval-averaged predictors. Exp Appl Acarol. 2017;73:439. https://doi.org/10.1007/s10493-017-0197-8.
Matuschka FR, Fischer P, Musgrave K, Richter D, Spielman A. Hosts on which nymphal Ixodes ricinus most abundantly feed. Am J Trop Med Hyg. 1991;44:100–7. https://doi.org/10.4269/ajtmh.1991.44.100.
Randolph SE, Storey K. Impact of microclimate on immature tick-rodent host interactions (Acari: Ixodidae): implications for parasite transmission. J Med Entomol. 1999;36:741–8. https://doi.org/10.1093/jmedent/36.6.741.
Klemola T, Sormunen JJ, Mojzer J, Mäkelä S, Vesterinen EJ. High tick abundance and diversity of tick-borne pathogens in a Finnish city. Urban Ecosystems. 2019;22:817–26. https://doi.org/10.1007/s11252-019-00854-w.
Jaenson TGT, Wilhelmsson P. First records of tick-borne pathogens in populations of the taiga tick Ixodes persulcatus in Sweden. Parasit Vectors. 2019;12:559. https://doi.org/10.1186/s13071-019-3813-0.
Gray JS, Dautel H, Estrada-Peña A, Kahl O, Lindgren E. Effects of climate change on ticks and tick-borne diseases in Europe. Interdiscip Perspect Infect Dis. 2009;2009:593232. https://doi.org/10.1155/2009/593232.
ECDC. European Centre for Disease Prevention and Control. Vector-borne diseases. https://www.ecdc.europa.eu/en/climate-change/climate-change-europe/vector-borne-diseases. 2021. Accessed 15 Jan 2022.
Sirotkin MB, Korenberg EI. Influence of abiotic factors on different developmental stages of the taiga tick Ixodes persulcatus and the Sheep Tick Ixodes ricinus. Entomological Review. 2018;98:496–513. https://doi.org/10.1134/S0013873818040115.
Daniel M, Malý M, Danielová V, Kříž B, Nuttall P. Abiotic predictors and annual seasonal dynamics of Ixodes ricinus, the major disease vector of Central Europe. Parasit Vectors. 2015;8:478. https://doi.org/10.1186/s13071-015-1092-y.
Daniel M, Danielová V, Fialová A, Malý M, Kříž B, Nuttall PA. Increased relative risk of tick-borne encephalitis in warmer weather. Front Cell Infect Microbiol. 2018;8:90. https://doi.org/10.3389/fcimb.2018.00090.
Lutta AS, Heysin EM, Shulman RE. Ticks in KarSSR and methods of control—Karel. Fil. Akad. Nauk SSSR: Petrozavodsk; 1959. pp. 68 [In Russian].
Brown L, Medlock J, Murray V. Impact of drought on vector-borne diseases—how does one manage the risk? Public Health. 2014;128:29–37. https://doi.org/10.1016/j.puhe.2013.09.006.
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. https://doi.org/10.1186/1756-3305-6-1.
Vollack K, Sodoudi S, Névir P, Müller K, Richter D. Influence of meteorological parameters during the preceding fall and winter on the questing activity of nymphal Ixodes ricinus ticks. Int J Biometeorol. 2017;61:1787–95. https://doi.org/10.1007/s00484-017-1362-9.
Guisan A, Zimmermann NE. Predictive habitat distribution models in ecology. Ecol Model. 2000;135:147–86. https://doi.org/10.1016/S0304-3800(00)00354-9.
Alfredsson M, Olafsson E, Eydal M, Unnsteinsdottir ER, Hansford K, Wint W, et al. Surveillance of Ixodes ricinus ticks (Acari: Ixodidae) in Iceland. Parasit Vectors. 2017;10:466. https://doi.org/10.1186/s13071-017-2375-2.
Williams HW, Cross DE, Crump HL, Drost CJ, Thomas CJ. Climate suitability for European ticks: assessing species distribution models against null models and projection under AR5 climate. Parasit Vectors. 2015;8:440. https://doi.org/10.1186/s13071-015-1046-4.
Alkishe AA, Peterson AT, Samy AM. Climate change influences on the potential geographic distribution of the disease vector tick Ixodes ricinus. PLoS ONE. 2017;12:e0189092. https://doi.org/10.1371/journal.pone.0189092.
Signorini M, Stensgaard A-S, Drigo M, Simonato G, Marcer F, Montarsi F, et al. Towards improved, cost-effective surveillance of Ixodes ricinus ticks and associated pathogens using species distribution modelling. Geospatial Health. 2019;14. https://doi.org/10.4081/gh.2019.745.
Kjær LJ, Soleng A, Edgar KS, Lindstedt HEH, Paulsen KM, Andreassen ÅK, et al. Predicting and mapping human risk of exposure to Ixodes ricinus nymphs using climatic and environmental data, Denmark, Norway and Sweden, 2016. Eurosurveillance. 2019;24:1800101. https://doi.org/10.2807/1560-7917.ES.2019.24.9.1800101.
Rochat E, Vuilleumier S, Aeby S, Greub G, Joost S. Nested species distribution models of chlamydiales in Ixodes ricinus (Tick) hosts in Switzerland. Appl Environ Microbiol. 2020;87. https://doi.org/10.1128/aem.01237-20.
Boehnke D, Brugger K, Pfäffle M, Sebastian P, Norra S, Petney T, et al. Estimating Ixodes ricinus densities on the landscape scale. Int J Health Geogr. 2015;14:23. https://doi.org/10.1186/s12942-015-0015-7.
Brugger K, Boehnke D, Petney T, Dobler G, Pfeffer M, Silaghi C, et al. A density map of the tick-borne encephalitis and lyme borreliosis vector Ixodes ricinus (Acari: Ixodidae) for Germany. J Med Entomol. 2016;53:1292–302. https://doi.org/10.1093/jme/tjw116.
Jung Kjær L, Soleng A, Edgar KS, Lindstedt HEH, Paulsen KM, Andreassen ÅK, et al. Predicting the spatial abundance of Ixodes ricinus ticks in southern Scandinavia using environmental and climatic data. Sci Rep. 2019;9:18144. https://doi.org/10.1038/s41598-019-54496-1.
Vogelgesang JR, Walter M, Kahl O, Rubel F, Brugger K. Long-term monitoring of the seasonal density of questing ixodid ticks in Vienna (Austria): setup and first results. Exp Appl Acarol. 2020;81:409–20. https://doi.org/10.1007/s10493-020-00511-4.
Yasyukevich VV, Kazakova EV, Popov IO, Semenov SM. Distribution of Ixodes ricinus L., 1758 and Ixodes persulcatus Shulze, 1930 (Parasitoformes, Ixodidae) in Russia and adjacent countries in view of observable climate changes. Dokl Earth Sci. 2009;427:1030–4. https://doi.org/10.1134/S1028334X09060312.
Popov IO, Popova EN. Identification of climate conditions restricting the distribution of the Taiga Tick Ixodes Persulcatus on the territory of Russia and in neighboring countries. Dokl Earth Sci. 2020;493:562–4. https://doi.org/10.1134/S1028334X20070144.
Hoch T, Monnet Y, Agoulon A. Influence of host migration between woodland and pasture on the population dynamics of the tick Ixodes ricinus: a modelling approach. Ecol Model. 2010;221:1798–806. https://doi.org/10.1016/j.ecolmodel.2010.04.008.
Dobson ADM, Randolph SE. Modelling the effects of recent changes in climate, host density and acaricide treatments on population dynamics of Ixodes ricinus in the UK. J Appl Ecol. 2011;48:1029–37. https://doi.org/10.1111/j.1365-2664.2011.02004.x.
Gilbert L, Maffey GL, Ramsay SL, Hester AJ. The effect of deer management on the abundance of Ixodes ricinus in Scotland. Ecol Appl. 2012;22:658–67.
Mysterud A, Easterday WR, Qviller L, Viljugrein H, Ytrehus B. Spatial and seasonal variation in the prevalence of Anaplasma phagocytophilum and Borrelia burgdorferi sensu lato in questing Ixodes ricinus ticks in Norway. Parasit Vectors. 2013;6:187. https://doi.org/10.1186/1756-3305-6-187.
Kruijff M, Hendrickx G, Wint W, Ginati A. Mapping habitats for vectors of infectious disease: VECMAP. 2011.
SYKE. Finnish Environment Institute. Corine maanpeite 2018. 2018. https://ckan.ymparisto.fi/dataset/corine-maanpeite-2018. Accessed 15 Sep 2021.
Semtner PJ, Howell DE, Hair JA. The ecology and behavior of the lone star tick (Acarina: Ixodidae) I. The relationship between vegetative habitat type and tick abundance and distribution in Cherokee Co., Oklahoma1. J Med Entomol. 1971;8:329–35. https://doi.org/10.1093/jmedent/8.3.329.
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:136–45. https://doi.org/10.1016/S1386-6532(02)00168-3.
Drew CA, Wiersma YF, Huettmann F. Predictive species and habitat modeling in landscape ecology: concepts and applications. New York: Springer; 2011.
Kass JM, Vilela B, Aiello-Lammens ME, Muscarella R, Merow C, Anderson RP. Wallace: a flexible platform for reproducible modeling of species niches and distributions built for community expansion. Methods Ecol Evol. 2018;9:1151–6. https://doi.org/10.1111/2041-210X.12945.
Aiello-Lammens ME, Boria RA, Radosavljevic A, Vilela B, Anderson RP. spThin: an R package for spatial thinning of species occurrence records for use in ecological niche models. Ecography. 2015;38:541–5. https://doi.org/10.1111/ecog.01132.
Aalto J, Pirinen P, Jylhä K. New gridded daily climatology of Finland: permutation-based uncertainty estimates and temporal trends in climate. J Geophys Res Atmos. 2016;121:3807–23. https://doi.org/10.1002/2015JD024651.
Pirinen P, Simola H, Aalto J, Kaukoranta JP, Karlsson P, Ruuhela R. Tilastoja Suomen ilmastosta 1981–2010 (Climatological statistics of Finland 1981–2010). Finnish Meteorological Institute Reports. 2012;2012.
Statistics Finland. https://tilastokeskus-kartta.swgis.fi/#2021. Accessed 10 Sep 2021.
NLS of Finland. National Land Survey of Finland. Digital elevation model 2019. 10 m x 10 m. CSC–IT Center for Science Ltd. 2019. http://urn.fi/urn:nbn:fi:csc-kata00001000000000000622. Accessed 15 Sep 2021.
Scharlemann JP, Benz D, Hay SI, Purse BV, Tatem AJ, Wint GR, et al. Global data for ecology and epidemiology: a novel algorithm for temporal Fourier processing MODIS data. PLoS ONE. 2008;3:e1408. https://doi.org/10.1371/journal.pone.0001408.
LUKE. National Resources Institute Finland. Snow-track data. 2021. https://www.riistakolmiot.fi/riistakolmio/. Accessed 15 Nov 2021.
Thuiller W, Lafourcade B, Engler R, Araújo M. BIOMOD—a platform for ensemble forecasting of species distributions. Ecography. 2009;32:369–73. https://doi.org/10.1111/j.1600-0587.2008.05742.x.
Thuiller W, Georges D, Gueguen M, Engler R, Breiner F. Biomod2: Ensemble platform for species distribution modeling. R package version 3.5.1. 2021. https://cran.r-project.org/package=biomod2. Accessed 15 Sep 2021.
Belsley DA, Kuh E, Welsch RE. Regression diagnostics: identifying influential data and sources of collinearity. New York: Wiley; 1980. https://doi.org/10.1002/0471725153.
Naimi B, Hamm NAS, Groen TA, Skidmore AK, Toxopeus AG. Where is positional uncertainty a problem for species distribution modelling? Ecography. 2014;37:191–203. https://doi.org/10.1111/j.1600-0587.2013.00205.x.
Chatterjee S, Hadi AS. Regression analysis by example. New Jersey: Wiley; 2013.
Sulaiman MS, Abood MM, Sinnakaudan SK, Shukor MR, You GQ, Chung XZ. Assessing and solving multicollinearity in sediment transport prediction models using principal component analysis. ISH J Hydraulic Eng. 2021;27:343–53. https://doi.org/10.1080/09715010.2019.1653799.
McCullagh P, Nelder JA. Generalized linear models. 2nd ed. London: Chapman & Hall; 1989.
Hastie TJ, Tibshirani R. Generalized additive models. 1st ed. London: Chapman & Hall; 1990.
Breiman L, Friedman JH, Olshen RA, Stone CJ. Classification and regression trees. New York: Chapman & Hall; 1984.
Ripley BD. Pattern recognition and neural networks. Cambridge University Press: Cambridge; 1996. https://doi.org/10.1017/CBO9780511812651.
Friedman JH. Multivariate adaptive regression splines. Ann Stat. 1991;19:1–67. https://doi.org/10.1214/aos/1176347963.
Ridgeway G. The state of boosting. Comput Sci Stat. 1999;31:172–81.
Breiman L. Random forests. Mach Learn. 2001;45:5–32. https://doi.org/10.1023/A:1010933404324.
Phillips S, Dudík M, Schapire R. A maximum entropy approach to species distribution modelling. Chemom Intell Lab Syst. 2004;146:406.
Zhao C, Gao F. A nested-loop Fisher discriminant analysis algorithm. Chemom Intell Lab Syst. 2015;146:396–406. https://doi.org/10.1016/j.chemolab.2015.06.008.
Elith J, Graham CH, Anderson PR, Dudík M, Ferrier S, Guisan A, et al. Novel methods improve prediction of species’ distributions from occurrence data. Ecography. 2006;29:129–51. https://doi.org/10.1111/j.2006.0906-7590.04596.x.
Wood SN. Fast stable restricted maximum likelihood and marginal likelihood estimation of semiparametric generalized linear models. J R Statis Soc Ser B Statis Methodol. 2011;73:3–36.
Fielding AH, Bell JF. A review of methods for the assessment of prediction errors in conservation presence/absence models. Environ Conserv. 1997;24:38–49. https://doi.org/10.1017/S0376892997000088.
Morán-Ordóñez A, Lahoz-Monfort JJ, Elith J, Wintle BA. Evaluating 318 continental-scale species distribution models over a 60-year prediction horizon: what factors influence the reliability of predictions? Glob Ecol Biogeogr. 2017;26:371–84. https://doi.org/10.1111/geb.12545.
Ministry of Social Affairs and Health. Raportteja ja muistioita 2021(16). Uusien vektorivälitteisten tautien mahdollinen saapuminen Suomeen ilmastonmuutoksen ja ihmisten liikkuvuuden kylkiäisinä. 2021. http://urn.fi/URN:ISBN:978-952-00-5402-1. Accessed 15 Sep 2021.
Araújo MB, New M. Ensemble forecasting of species distributions. Trends Ecol Evol. 2007;22:42–7. https://doi.org/10.1016/j.tree.2006.09.010.
Willcock S, Hooftman DAP, Blanchard R, Dawson TP, Hickler T, Lindeskog M, et al. Ensembles of ecosystem service models can improve accuracy and indicate uncertainty. Sci Total Environ. 2020;747:141006. https://doi.org/10.1016/j.scitotenv.2020.141006.
Leta S, Fetene E, Mulatu T, Amenu K, Jaleta MB, Beyene TJ, et al. Modeling the global distribution of Culicoides imicola: an Ensemble approach. Sci Rep. 2019;9:14187. https://doi.org/10.1038/s41598-019-50765-1.
Messina JP, Brady OJ, Golding N, Kraemer MUG, Wint GRW, Ray SE, et al. The current and future global distribution and population at risk of dengue. Nat Microbiol. 2019;4:1508–15. https://doi.org/10.1038/s41564-019-0476-8.
Sinka ME, Pironon S, Massey NC, Longbottom J, Hemingway J, Moyes CL, et al. A new malaria vector in Africa: predicting the expansion range of Anopheles stephensi and identifying the urban populations at risk. Proc Natl Acad Sci USA. 2020;117:24900–8. https://doi.org/10.1073/pnas.2003976117.
Franklin J. Mapping species distributions: spatial inference and prediction. Cambridge: Cambridge University Press; 2010.
Araújo MB, Anderson RP, Barbosa AM, Beale CM, Dormann CF, Early R, et al. Standards for distribution models in biodiversity assessments. Sci Adv. 2019;5:eaat4858. https://doi.org/10.1126/sciadv.aat4858.
Beck J, Böller M, Erhardt A, Schwanghart W. Spatial bias in the GBIF database and its effect on modeling species’ geographic distributions. Ecol Informat. 2013;19:10–5. https://doi.org/10.1016/j.ecoinf.2013.11.002.
Swets JA. Measuring the accuracy of diagnostic systems. Science. 1988;240:1285–93. https://doi.org/10.1126/science.3287615.
Marsot M, Henry PY, Vourc’h G, Gasqui P, Ferquel E, Laignel J, et al. Which forest bird species are the main hosts of the tick, Ixodes ricinus, the vector of Borrelia burgdorferi sensu lato, during the breeding season? Int J Parasitol. 2012;42:781–8. https://doi.org/10.1016/j.ijpara.2012.05.010.
Jore S, Vanwambeke SO, Viljugrein H, Isaksen K, Kristoffersen AB, Woldehiwet Z, et al. Climate and environmental change drives Ixodes ricinus geographical expansion at the northern range margin. Parasit Vectors. 2014;7:11. https://doi.org/10.1186/1756-3305-7-11.
Hofmeester TR, Sprong H, Jansen PA, Prins HHT, van Wieren SE. Deer presence rather than abundance determines the population density of the sheep tick, Ixodes ricinus, in Dutch forests. Parasit Vectors. 2017;10:433. https://doi.org/10.1186/s13071-017-2370-7.
Mysterud A, Hügli C, Viljugrein H. Tick infestation on medium–large-sized mammalian hosts: are all equally suitable to Ixodes ricinus adults? Parasit Vectors. 2021;14:254. https://doi.org/10.1186/s13071-021-04775-6.
Tälleklint L, Jaenson TG. Infestation of mammals by Ixodes ricinus ticks (Acari: Ixodidae) in south-central Sweden. Exp Appl Acarol. 1997;21:755–71. https://doi.org/10.1023/a:1018473122070.
Jaenson TGT, Petersson EH, Jaenson DGE, Kindberg J, Pettersson JHO, Hjertqvist M, et al. The importance of wildlife in the ecology and epidemiology of the TBE virus in Sweden: incidence of human TBE correlates with abundance of deer and hares. Parasit Vectors. 2018;11:477. https://doi.org/10.1186/s13071-018-3057-4.
Wodecka B, Skotarczak B. Identification of host blood-meal sources and Borrelia in field-collected Ixodes ricinus ticks in north-western Poland. Ann Agric Environ Med. 2016;23:59–63. https://doi.org/10.5604/12321966.1196853.
Cadenas FM, Rais O, Humair P-F, Douet V, Moret J, Gern L. Identification of host bloodmeal source and Borrelia burgdorferi Sensu Lato in Field-Collected Ixodes ricinus Ticks in Chaumont (Switzerland). J Med Entomol. 2007;44:1109–17. https://doi.org/10.1093/jmedent/44.6.1109.
Scharlemann JPW, Johnson PJ, Smith AA, Macdonald DW, Randolph SE. Trends in ixodid tick abundance and distribution in Great Britain. Med Vet Entomol. 2008;22:238–47. https://doi.org/10.1111/j.1365-2915.2008.00734.x.
Jaenson TGT, Jaenson DGE, Eisen L, Petersson E, Lindgren E. Changes in the geographical distribution and abundance of the tick Ixodes ricinus during the past 30 years in Sweden. Parasit Vectors. 2012;5:8. https://doi.org/10.1186/1756-3305-5-8.
Tokarevich NK, Tronin AA, Blinova OV, Buzinov RV, Boltenkov VP, Yurasova ED, et al. The impact of climate change on the expansion of Ixodes persulcatus habitat and the incidence of tick-borne encephalitis in the north of European Russia. Glob Health Action. 2011;4:8448. https://doi.org/10.3402/gha.v4i0.8448.
Bugmyrin SV, Bespyatova LA, Korotkov YS. Long-term dynamics of Ixodes persulcatus (Acari: Ixodidae) abundance in the north–west of its range (Karelia, Russia). Exp Appl Acarol. 2019;77:229–40. https://doi.org/10.1007/s10493-019-00342-y.
Lindgren E, Gustafson R. Tick-borne encephalitis in Sweden and climate change. Lancet. 2001;358:16–8. https://doi.org/10.1016/s0140-6736(00)05250-8.
Kim HC, Han SH, Chong ST, Klein TA, Choi C-Y, Nam H-Y, et al. Ticks collected from selected mammalian hosts surveyed in the Republic of Korea during 2008–2009. Korean J Parasitol. 2011;49:331–5. https://doi.org/10.3347/kjp.2011.49.3.331.
LUKE. National Resources Institute Finland. Valkohäntäpeurakanta talvella 2020 – 2021. 2021. http://wordpress1.luke.fi/riistahavainnot-hirvielaimet/wp-content/uploads/sites/5/2021/03/Valkoh%C3%A4nt%C3%A4peurakanta_talvella_2020_2021.pdf. Accessed 15 Dec 2021.
Bespyatova LA. On the Distribution of the Castor Bean Tick Ixodes ricinus (Acarina, Ixodidae) in the Republic of Karelia, Russia. Entomol Revue. 2021;101:581–91. https://doi.org/10.1134/S0013873821040084.
Bugmyrin SV, Bespyatova LA, Korotkov YS, Burenkova LA, Belova OA, Romanova L, et al. Distribution of Ixodes ricinus and I. persulcatus ticks in southern Karelia (Russia). Ticks Tick Borne Dis. 2013;4:57–62. https://doi.org/10.1016/j.ttbdis.2012.07.004.
Acknowledgements
We would like to thank the researchers and institutions who provided and collected data for this study. We thank Saana Sipari, Lauri Kirjalainen, Petteri Puonti, Roope Nykänen, Robert Castren, Kirsi Aaltonen, Ravi Kant, Ella Lintunen, and Amalia Kokkonen who participated in data collection efforts during early summer 2021.
Funding
This study was funded by the Doctoral Programme in Interdisciplinary Environmental Sciences (DENVI) of the University of Helsinki and by the Academy of Finland through the VECLIMIT project (decision No #329323).
Author information
Authors and Affiliations
Contributions
RU, MS, GH, CM and JA designed and supervised the study; RU, GH and MS designed the sampling; RU, JS, EK and EMK conducted the sampling; RU performed the species identifications and pathogen screening; RU analysed the data; RU drafted the manuscript; all authors contributed to the final manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors have declared that no competing interests exist.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Figure S1.
(a) The sampling strategy for new collections in 2021 was created based on the following criteria. Subdivisions of landscape areas (Area1–Area4), CORINE land cover 2018, a 5-km buffer around existing I. persulcatus occurrences (grey circles), and a 500-m buffer around roads were used to delimit the four sampling areas (light grey lines). For each sampling area, a random sample of 25 collection locations was created depending on the relative shares of forest and meadow categories in each area. (b) The map showing the 2021 results indicates the locations where I. ricinus was found with B. burgdorferi (s.l.)-positive locations.
Additional file 2: Figure S2.
The range (lines) and mean (dots) of model performances over 50 model runs in each model algorithm estimating habitat suitabilities for I. ricinus in different variable compositions: (a) environmental only, (b) host only, (c) environmental and host, and (d) environmental, host, and suitability for I. ricinus.
Additional file 3: Figure S3.
The range (lines) and mean (dots) of model performances over 50 model runs in each model algorithm estimating habitat suitabilities for I. persulcatus in different variable compositions: (a) environmental only, (b) host only, (c) environmental and host, and (d) environmental, host, and suitability for I. ricinus.
Additional file 4: Figure S4.
The relative contributions of the explanatory variables in the data set of (a) host only, (b) environment only based on the mean ensemble model.
Additional file 5: Table S1.
The number of times each model contributed to the final ensemble in different data sets.
Additional file 6: Figure S5.
Partial dependency plots for (a) I. ricinus and (b) I. persulcatus solely based on environmental data.
Additional file 7: Figure S6.
Partial dependency plots for (a) I. ricinus and (b) I. persulcatus solely based on host data.
Additional file 8: Figure S7.
Partial dependency plots for (a) I. ricinus and (b) I. persulcatus based on combined host and environmental data, and habitat suitability data for the other tick species.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Uusitalo, R., Siljander, M., Lindén, A. et al. Predicting habitat suitability for Ixodes ricinus and Ixodes persulcatus ticks in Finland. Parasites Vectors 15, 310 (2022). https://doi.org/10.1186/s13071-022-05410-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13071-022-05410-8