- Research
- Open access
- Published:
Modeling platform to assess the effectiveness of single and integrated Ixodes scapularis tick control methods
Parasites & Vectors volume 17, Article number: 339 (2024)
Abstract
Background
Lyme disease continues to expand in Canada and the USA and no single intervention is likely to curb the epidemic.
Methods
We propose a platform to quantitatively assess the effectiveness of a subset of Ixodes scapularis tick management approaches. The platform allows us to assess the impact of different control treatments, conducted either individually (single interventions) or in combination (combined efforts), with varying timings and durations. Interventions include three low environmental toxicity measures in differing combinations, namely reductions in white-tailed deer (Odocoileus virginianus) populations, broadcast area-application of the entomopathogenic fungus Metarhizium anisopliae, and fipronil-based rodent-targeted bait boxes. To assess the impact of these control efforts, we calibrated a process-based mathematical model to data collected from residential properties in the town of Redding, southwestern Connecticut, where an integrated tick management program to reduce I.xodes scapularis nymphs was conducted from 2013 through 2016. We estimated parameters mechanistically for each of the three treatments, simulated multiple combinations and timings of interventions, and computed the resulting percent reduction of the nymphal peak and of the area under the phenology curve.
Results
Simulation outputs suggest that the three-treatment combination and the bait boxes–deer reduction combination had the overall highest impacts on suppressing I. scapularis nymphs. All (single or combined) interventions were more efficacious when implemented for a higher number of years. When implemented for at least 4 years, most interventions (except the single application of the entomopathogenic fungus) were predicted to strongly reduce the nymphal peak compared with the no intervention scenario. Finally, we determined the optimal period to apply the entomopathogenic fungus in residential yards, depending on the number of applications.
Conclusions
Computer simulation is a powerful tool to identify the optimal deployment of individual and combined tick management approaches, which can synergistically contribute to short-to-long-term, costeffective, and sustainable control of tick-borne diseases in integrated tick management (ITM) interventions.
Graphical Abstract
![](http://media.springernature.com/lw685/springer-static/image/art%3A10.1186%2Fs13071-024-06387-2/MediaObjects/13071_2024_6387_Figa_HTML.png)
Background
Tick-borne diseases are expanding globally [1,2,3]. In Canada and the USA, the number of cases of one of these zoonotic diseases, Lyme disease, is on a steady rise [4,5,6,7,8]. This long-term trend has been attributed to the range expansion of the primary vector to humans, the blacklegged tick, Ixodes scapularis Say [9,10,11]. Complementary approaches have been developed to reduce human exposure to the causative agent of Lyme disease, Borrelia burgdorferi. One group of interventions targets reduction in density of ticks in the environment by environmental and host-targeted acaricides [12,13,14,15,16,17,18], another aims to decrease infection prevalence of hosts and ticks [15, 18, 19], and the other aims to reduce human-tick encounters [17, 20]. Here we focus on methods suppressing I. scapularis abundance in the environment, which include area-wide application of natural and synthetic acaricidal compounds, landscape and vegetative cover modifications and management [14], host-targeted synthetic acaricides such as fipronil-based rodent-targeted bait boxes [15, 16], biological control [17, 18, 20, 21] by entomopathogenic fungal agents [18, 20, 21], and reductions of white-tailed deer, Odocoileus virginianus, populations [22, 23].
Integrated tick management (ITM) has been proposed to be more effective than single measures [24, 25]. However, there is limited empirical support [12, 13, 26,27,28] and knowledge gaps on the effectiveness of such integrated tick control methods remain [20, 29]. A significant barrier in evaluation and optimization of ITM approaches is that intervention outcomes are often assessed as the reduction in density of nymphs/infected nymphs, but there is a lack of mechanistic understanding about what component of the enzootic cycle was impacted. For example, when assessing the effect of white-tailed deer management, the direct effect of deer reduction is often not assessed. Rather, the focus tends to be on deer reduction impacts to the entomological inoculation rate (measured as the density of questing I. scapularis nymphs infected with the target pathogen), which is usually the metric of human health risk. While useful in assessing effectiveness of the intervention in reducing human risk in the particular context assessed, this black-box evaluation limits our ability to expand the application to other settings or optimizing different arms of an ITM approach.
Here we present a robust, empirically calibrated modeling platform that can simulate and compare single and combined interventions, including differences in timing of delivery. Mathematical modeling approaches provide a mechanistic understanding of the tick–host system. They also represent the system’s continuous (i.e., not restricted to specific periods) dynamics in simulations, and thus capture the temporal complexity of intervention implementation. Using an adaptation of an extensively validated climate-driven model [30, 31], we propose a modeling platform to evaluate the impact of individual and multiple tick management approaches on the density of questing I. scapularis nymphs (DON), a key component of Lyme disease risk [32]. Hereafter we refer to I. scapularis as ‘tick’ and specific I. scapularis stages as ‘larvae’, ‘nymphs’ and ‘adults.’
Methods
General approach
The platform integrates a process-based mechanistic model with local weather forcing and a set of three tick-reduction interventions (Fig. 1). We calibrated the mechanistic model using data from previously published studies by our group [13, 28] involving residential yards that underwent single or combined tick-reduction interventions, or were included as reference (no intervention) properties (see Supplementary Table S1). First, we calibrated the model parameters related to tick phenology (in absence of intervention) to the data from reference properties. Second, using the data from properties where interventions were implemented, we calibrated more specifically three parameters quantifying the effects of the three tick-reduction interventions on biological mechanisms. Depending on the duration and intensity of the interventions, we modeled them either as “pulse” (instantaneous perturbation of the system, followed by a gradual return to previous state) or “press” interventions (sustained perturbation of the system at constant intensity) [33]. Third, we ran the fully calibrated model multiple times to represent various tick-reduction intervention scenarios and timings. Finally, we analyzed the outcomes of these data-informed model predictions to quantitatively compare the efficacy of these different ITM approaches for the reduction in the DON.
Sketch diagram of the modeling platform. Blue underlines highlight the endogenous variables of the Ogden et al. [31] process-based model that are affected by weather conditions: time delay for the pre-eclosion period of eggs (q); host-finding probability for questing larvae, questing nymphs, and questing adults (λQL, λQN, and λQA, respectively); temperature-variable factors for questing activity of immature ticks and adult ticks (Θi and Θa, respectively); time delay for engorged larva to nymph development, or premolt period of engorged larvae (s); time delay for engorged nymph to adult development, or premolt period of engorged nymphs (v); and time delay for the pre-oviposition period (x). The orange, red, and green arrows point to the exogenous and endogenous variables that are affected by the three single interventions (see main text): white-tailed deer population reduction (Deer), broadcast area-application of the entomopathogenic fungus Metarhizium anisopliae (Met52), and fipronil-based small rodent bait boxes (FipBox), respectively. A full description of all the parameters of the process-based model is included in the Supplementary material (level variables in Table S2; parameters in Tables S3 and S4; and discrete equations in Table S5)
Process-based model
To represent tick population dynamics, we built an R-version of the Ogden et al. [30] process-based, mechanistic model. The model is based on a system of 12 coupled, nonlinear ordinary differential equations, each representing a given life stage or activity level of the tick population: eggs; hardening, questing, feeding, and engorged larvae; questing, feeding, and engorged nymphs; questing adults; and feeding, engorged, and egg-laying adult females (see Supplementary Tables S2–S5). The mathematical tool was modified to include the weather influence on the mortality (or survival) rate of questing nymphs and the host finding probability, which were both assumed to be constant in the Ogden et al. [30] version of the model.
Integrated vector control data used for calibration
Data used to calibrate the model were derived from an intervention study previously implemented on residential properties in the town of Redding (41.3044°N, 73.3928°W), Fairfield County, in southwestern Connecticut. Specifically, an ITM program to reduce questing nymphs was conducted from January 2013 through September 2016 [13, 28]. During that period, a total of 41 properties were selected to run the program, including 12 reference (no intervention) properties and 29 properties with implemented interventions. Interventions included single treatments and two- and three-way combinations of three measures, namely white-tailed deer population reduction (Deer), broadcast area-application of the entomopathogenic fungus Metarhizium anisopliae (Met52), and fipronil-based small rodent bait boxes (FipBox). Bait boxes were deployed twice at each property during the summer season (May–August) ahead of peak nymphal and larval questing activity, respectively. Sample sizes for each intervention and combination thereof are described in the Supplementary material (Table S1).
Weather data
To characterize weather conditions, we used three data sources: satellite-gauge combined climatic records (specifically, gridded near-surface air temperature datasets), weather observations, and reanalysis data (specifically, gridded near-surface humidity data). A full description of weather datasets used in simulations is provided in the Supplementary Note 1 and Fig. S1. The model was forced with temperature, humidity, and day length data; the first two variables were corrected to reflect environmental conditions under the leaf litter on the basis of linear transformations from near-surface air temperature and relative humidity to leaf litter temperature and moisture. The linear transformations were determined using data collected at 1 m above ground surface (HOBO U23-001 Pro v2 dataloggers) and in the leaf litter (iButtons) in a previous study (Supplementary Note 1). The conversion aimed to account for ticks spending the majority of their life cycle under leaf litter in forested habitats, where temperature and relative humidity are more moderate than for surrounding air [34, 35]. The adapted Ogden et al. [30] mechanistic model implemented here includes: temperature-dependent development rates for eggs and engorged larvae, nymphs, and females; temperature- and humidity-dependent survival rates of free-living ticks ; temperature-, day length-, and host-density-dependent host-finding rates; and density-dependent survival of all tick stages during engorgement.
Host data
The 2013–2016 small mammal, live trapping campaign in Redding [13, 28] yielded a grand total of 1298 hosts for immature tick stages captured from residential properties, including white-footed mouse (Peromyscus leucopus; 86.4% of captures), eastern chipmunk (Tamias striatus; 9.0%), eastern meadow vole (Microtus pennsylvanicus; 2.2%), northern short-tailed shrew (Blarina brevicauda; 1.9%), American red squirrel (Tamiasciurus hudsonicus; 0.2%), and eastern gray squirrel (Sciurus carolinensis; 0.1%). We restricted hosts for immature ticks to white-footed mice (as they were the dominant small mammal) and white-tailed deer as the only host for adults in the simulations; both host populations were assumed to be constant.
Simulating tick population dynamics
The process-based model was run for a time horizon from 1 January 1950 to 31 December 2020 (25,933 days). This long-term period considers an initial 10-year burn-in period for the model to reach equilibrium and additional 10 years to assess day-to-day fluctuations of the tick population after the model had reached stable conditions. The remaining 50 years of the simulation run included the final 2013–2016 experimentation period, for which data are available on the observed densities of questing nymphs expressed per 100 m2, from a total of 1244 tick dragging collections conducted over the research period [13, 28].
Calibration of the model’s tick phenology parameters
Parameters related to tick phenology in the mathematical model were first calibrated to questing nymphal data collected at the reference properties (we averaged the observed questing nymphal data of the 12 properties) by simultaneously minimizing three statistical errors: the percent error (BE) [36], the root mean square error (RMSE) [37], and the general tendency (Pbias) [38] and maximizing the Nash–Sutcliffe efficiency (NSE) [39].
Efficacy of interventions
We then used data collected from the 29 properties where single or combinations of interventions were implemented (see Table S1) to estimate their efficacy. For a given year at a given property, the Deer treatment was assumed to multiply the host-finding probability for questing adults (λQA) between 1 January and 31 December by a factor fD, where fD was calibrated to the data (press function).
We considered that the Met52 treatment affected the per capita mortality rate of questing nymphs (μQN) and questing adults (μQA) of the process-based model as a pulse (instantaneous perturbation) function. Immediately after the treatment was applied, these two parameters were increased by a factor fM (calibrated to the data), followed by a gradual decay in mortality over the effective residual life of Met52, back toward their baseline (without intervention) values at a rate given by:
where Δt depicts the number of days since the treatment application and where ε is a tuning parameter set to 0.1. ε is a tuning parameter that allows the effect of treatment on mortality to decrease over time, as shown in Supplementary Fig. S2 (upper panel). From formula (1), ε affects both the strength of the effect at treatment application (that is, fM (1 − ε2) at Δt = 0) and the shape of the decrease or derivative (that is, fM ε(28− Δt)/14 log(ε)/14). Therefore, ε cannot exactly be set to 0 but needs to be close to 0 to obtain an effect value close to fM at Δt = 0. This is why we arbitrarily set it to 0.1. The effective residual life of Met52 was assumed to be 4 weeks (or 28 days) [40]. The first and second application rounds took place in early June and early July respectively, leading to an impact of Met52 only treatment lasting for two consecutive months (early June to early August), consistent with the broadcast applications during the 2013–2016 control campaign.
Similarly, the FipBox treatment was assumed to impact the per capita mortality rate of feeding larvae (μFL) and feeding nymphs (μFN) of the process-based model as a pulse function (see Fig. S2). These two parameters were increased by a factor fF (calibrated to the data) after boxes were deployed, followed by a gradual decay in mortality over the effective residual life of the compound, back toward their baseline (without intervention) values at a rate given by:
where Δt depicts the number of days since the treatment application and ε is a tuning parameter set to 0.1. The effective residual life of the compound was assumed to be 8 weeks (or 56 days). To account for two deployments of bait boxes during the campaign, we simulated an increase in tick mortality in early May and again after 8 weeks (i.e., approximately 16 weeks in total), the period after which bait boxes were retrieved.
Parameters fD, fM, and fF were simultaneously estimated by minimizing the RMSE for all 29 properties with interventions, which enabled assessment of the effect of each individual intervention despite the simultaneous implementation in some of the properties.
Simulation of the effectiveness of ITM approaches
Using the values of fD, fM, and fF estimated from the data, we simulated the implementation of all combinations of the three treatments: Deer only (D), Met52 only (M), FipBox only (F), Deer + Met52 (D + M), Deer + FipBox (D + F), Met52 + FipBox (M + F), and Deer + Met52 + FipBox (D + M + F).
Simulation of alternative number and timing of interventions
The last set of simulations included the analysis of changes in the Met52 only treatment applications. Starting from the baseline scenario with two applications of treatment in early June and early July, respectively, we assessed the effects of shifting the timing of application and reducing the number of applications. We simulated a change in the dates of both Met52 only applications: the first application taking place between late April and mid-July, and the second application 29 days later. We also simulated these alternatives in a scenario where only the first application of Met52 only was performed.
Metrics for the percent reduction
In these different scenarios, we calculated the percent reduction of the predicted (simulated) nymphal peak, i.e., the percent reduction in the number of feeding nymphs per host at the peak of each nymphal season associated with each single treatment and ITM approaches. We also calculated the percentage change in the area under the phenology curve, which accounts for the seasonal variation in the densities of feeding nymphs compared with reference plots, for simulated single treatments and combined efforts.
Results
Model calibration
For all reference properties (ITM reference), simulation outputs of the process-based model yielded a BE value of 29.1%, a RMSE value of 1.12, a Pbias value of −64.9, and a NSE value of −0.98, when the model was forced with near-surface air temperature weather station data. Predictions of the calibrated model are shown in Figs. 2 and S3. In the absence of interventions, the model predicted yearly feeding nymphal peaks at around three nymphs per mouse (Figs. 2 and S3).
Daily time series of feeding ticks on host simulated by the dynamic tick population model for the historical period 2013–2019, for in situ conditions corrected to reflect temperature and moisture conditions under leaf-litter, for weather-dependent host finding probability, and for the group of properties included in ITM reference. Black and red time series depict the simulation outputs when the model is forced with near-surface air temperature weather station data. The blue dots depict the average number of observed densities of nymphs per 100 m2
Efficacy of interventions on reduction of the peak in feeding nymphs
The point estimate for fM was 1.25, meaning that the Met52 only treatment was estimated to multiply the daily mortality of questing nymphs and adults by 2.25 on the first day of application, an effect then decreasing with time (Eq. 1). The point estimate for fF was 4.11, meaning that the FipBox treatment was estimated to multiply the daily mortality of feeding larvae and nymph by 5.11 on the first day, an effect then decreasing with time (Eq. 2) (Fig. S2). The point estimate for fD was 0.985, meaning that the Deer only intervention was estimated to reduce the host-finding rate for adults in yards by 98.5%. In model simulations utilizing these parameters, the efficacy of all interventions was higher and persisted for a longer time (Figs. 3, 4 and S4) when they were applied for a higher number of consecutive years. Combinations of multiple interventions were always more efficacious than single interventions (Fig. 3).
Percent reduction of the height of the feeding nymphal peak on the last year of intervention, compared with the no intervention scenario. For each intervention (single and combined), different alternatives for the number of years of implementation (from 1 to 4 years) are presented. D: Deer only; M: Met52 only; F: FipBox only; D + M: Deer + Met52; D + F: Deer + FipBox; M + F: Met52 + FipBox; D + M + F: Deer + Met52 + FipBox
Evolution of the percent reduction of the height of the feeding nymphs annual peak for each intervention, compared with the no intervention scenario. In our model simulations, interventions are all implemented in 2013, for 1 to 4 years. For each intervention and number of years, the dot represents the last year of intervention, with the line representing changes in effect size with years after the intervention ended. D: Deer only; M: Met52 only; F: FipBox only; D + M: Deer + Met52; D + F: Deer + FipBox; M + F: Met52 + FipBox; D + M + F: Deer + Met52 + FipBox
Among single interventions implemented for only 1 year, FipBox only was the most efficacious as it decreased the feeding nymphal peak on that year by 53.6% compared with the no intervention scenario (by killing the nymphs feeding on hosts), versus 6.4% and 4.7% for the Deer only and Met52 only, respectively (Fig. 3). However, because the FipBox treatment affects both larvae and nymphs on host, it showed its peak efficacy on the nymphal population the year following the first year of intervention (56.7% reduction of the nymphal peak) (Fig. 4). Similarly, because of its effects on adult ticks, the Deer only intervention reached its peak efficacy 2 years after the first year of intervention, with a 19.6% decrease in nymphal peak (Figs. 4 and S4).
When implemented for several years, Deer only reached a similar efficacy to FipBox only, i.e., a 97.0% versus 94.3% reduction, respectively, of the nymphal peak on the last year of 4 years of intervention, compared with the no intervention scenario. The Met52 only treatment was less efficacious even after 4 years of implementation (16.1% reduction) (Figs. 3, 4).
The combination of the three interventions (Deer + Met52 + FipBox) was always the most efficacious, leading to a decrease of the nymphal peak of 59.6% (after 1 year of implementation) to 99.8% (after 4 years of implementation), compared with the no intervention scenario. However, it was closely followed by the Deer + FipBox intervention, whose efficacy was never more than 1% lower than the Deer + Met52 + FipBox scenario. Assuming four consecutive years of implementation, all interventions (single and combined) reached a reduction in the nymphal peak of more than 94% compared with no intervention, except Met52 only (Figs. 3, 4).
Comparison of percent reduction metrics
The results described above were qualitatively similar when we considered the reduction in the area under the feeding nymphs phenology curve instead of the height of the feeding nymphal peak (Figs. S5 and S6).
Changes in the timing of Met52 only interventions
From our model simulations, the application of the Met52 only treatment on 19 May followed by a second application 29 days later led to the highest reduction of the nymphal peak that year (− 5.85% compared with the no intervention scenario) (Fig. 5). However, applying the Met52 only treatment only once in mid-June led to almost similar results (−5.13% compared with the no intervention scenario) while reducing the amount of input needed (Fig. 5). The advantage of applying a second Met52 only treatment was different when measuring the area under the phenology curve (Fig. S7), with a maximum effect size of −3.61% for two treatments, instead of −1.98% for one treatment, compared with the no intervention scenario.
Percent reduction of the height of the feeding nymphal peak in the only year of intervention, compared with the no intervention scenario, for different variations of the Met52 only treatment. One or two Met52 applications are simulated, starting at different dates between 12 April and 14 July. If done, the second Met52 application always takes place 29 days after the first one. The application of Met52 after the nymphal peak (here, after 7 July) has, as expected, no impact on peak height
Discussion
We found that combined treatments were globally more efficacious than single treatments in reducing DON in residential yards. Among single treatments, our results suggest that deer reduction (Deer) and fipronil-based small rodent bait boxes (FipBox) were more efficacious than the broadcast area-application of the entomopathogenic soil-borne fungus Metarhizium anisopliae (Met52). Because interventions targeted different tick life stages, they impacted the nymphal tick outcome at different time lags, with an immediate effect of Met52, an effect of FipBox reaching its maximum after a 1-year lag, and deer reduction 2 years after the intervention [17]. Globally, higher number of subsequent years of (single or combined) interventions led to better efficacy. Most interventions (except Met52 only) reached a reduction in the feeding nymph peak of more than 94% compared with the no intervention scenario when implemented for at least 4 years. Moreover, we found that the optimal date for a single application of Met52 in Connecticut residential yards would be around mid-June. Preceding this application by another application 1 month earlier (in mid-May) would only slightly reduce the height of the feeding nymph peak (maximum exposure potential) but would reduce by around half the area under the phenology curve across the whole season. The area under the phenology curve is a cumulative measure of human exposure risk because exposure to questing nymphs can occur throughout the period of nymphal activity (results shown in the supplement).
Previous empirical studies investigated the effect of interventions on the DON, tick burden on hosts, or the prevalence of tick infection by various pathogens [13, 17, 18, 22, 28, 41,42,43]. Some published models also theoretically simulated the implementation of interventions (vaccination) in a tick–host system [44, 45]. Our study combined both of these approaches by using data collected from field experiments to estimate intervention effect parameters in a mechanistic model. This method is innovative for tick–host systems, while more common for other vector-borne disease systems [46,47,48], and has several advantages. First, it provides a mechanistic and dynamic understanding of interventions, in contrast with empirical-only studies, which only assess differences in global output (e.g., the density of questing nymphs at a single point in time). For example, we estimated that Met52 only in our study multiplied the daily mortality rate of questing nymphs and adults by 2.25 on the first day of application and could estimate the cumulative impact of the intervention during multiple years. Second, once mechanistic intervention effect parameters are calibrated to the data, it is possible to explore a wide range of simulated treatment scenarios in terms of timing or combinations of interventions. Simulations are hence easier than implementing a large number of treatment options in the field, and can narrow down which combinations, intensity, and timing of interventions should be further tested in experimental studies.
The effect sizes identified in our model were consistent with some previous studies. On the year of intervention, we predicted similar reductions of the number of feeding nymphs per mouse (nymphal burden) than a previously published empirical study [49] for both FipBox (around 50%) and Met52 (less than 15%, not significant in the empirical study). Our results were also consistent with Ostfeld et al. [42], who found a stronger effect of FipBox than Met52 on the nymphal burden on mice.
In our study, despite the fact that the deer population was not entirely removed [28], we predicted a strong effect of this intervention on the nymphal peak when it was implemented for several years. However, our modeling framework did not include the infection process, and the effect of deer population reduction (rather than elimination) on nymph infection prevalence and the density of infected nymphs is still unclear [50, 51]. Immediately following deer reduction, in the absence of numerous large-bodied hosts, sampling efforts for remaining questing ticks often result in a perceived, temporary amplification in both abundance and pathogen infection prevalence as juveniles engaging in questing behavior and are then more likely to obtain a bloodmeal from a reservoir-competent host [22, 50,51,52].
Our study has some limitations. First, Gaff et al. [53] found that assuming constant instead of varying host population densities—as we did—can affect the density of nymphs predicted by a model, although it hardly changes the predicted number of feeding ticks per host, which was the output we considered here. Second, we did not incorporate host movements in our model as in Li et al. [54,55,], Wang et al. [55] and Wang et al. [56]. In the future, adding host behavioral components to the model would more directly describe human exposure to tick-borne disease hazard. Moreover, expanding the model to include infection dynamics and other sets of interventions, including vaccines and reductions in tick habitat suitability through landscape and vegetation management, will allow for estimates of efficacy in reducing the density of infected nymphs. Integrating the costs of interventions would allow for a more thorough cost-effectiveness analysis of tick- and host-targeted interventions in residential yards [57].
Conclusions
Computer simulations allow us to identify optimal control targets to minimize exposure to tick bites and the risk for tick-borne diseases. Our long-term goal is to provide information to stakeholders on the optimal deployment of individual and combined tick management approaches that can synergistically contribute to short-to-long-term, cost-effective, and sustainable control of tick vectors using integrated vector management guidelines.
Availability of data and materials
The data that support the findings of this study are available from the corresponding author upon reasonable request. All images were generated by the authors, so no consent for publication was required.
References
Eisen RJ, Eisen L, Beard CB. County-scale distribution of Ixodes scapularis and Ixodes pacificus (Acari: Ixodidae) in the continental United States. J Med Entomol. 2016;53:349–86.
Ogden NH, Lindsay LR, Leighton PA. Predicting the rate of invasion of the agent of Lyme disease Borrelia burgdorferi. J Appl Ecol. 2013;50:510–8.
Ogden NH, Lindsay LR. Effects of climate and climate change on vectors and vector-borne diseases: ticks are different. Trends Parasitol. 2016;32:646–56.
Rosenberg R, Lindsey NP, Fischer M, Gregory CJ, Hinckley AF, Mead PS, Paz-Bailey G, Waterman SH, Drexler NA, Kersh GJ, Hooks H, Partridge SK, Visser SN, Beard CB, Petersen LR. Vital signs: trends in reported vectorborne disease cases—United States and Territories, 2004–2016. Morb Mortal Wkly Rep. 2018;67:496.
Sonenshine DE. Range expansion of tick disease vectors in North America: implications for spread of tick-borne disease. Int J Environ Res Public Health. 2018;15:478.
Eisen L. Stemming the rising tide of human-biting ticks and tickborne diseases, United States. Emerg Infect Dis. 2020;26:641–7.
Kugeler KJ, Schwartz AM, Delorey MJ, Mead PS, Hinckley AF. Estimating the frequency of lyme disease diagnoses, United States, 2010–2018. Emerg Infect Dis. 2021;27:616–9.
Bisanzio D, Fernández MP, Martello E, Reithinger R, Diuk-Wasser MA. Current and future spatiotemporal patterns of lyme disease reporting in the Northeastern United States. JAMA Netw Open. 2020;3:e200319.
Khatchikian CE, Prusinski MA, Stone M, Backenson PB, Wang IN, Foley E, Seifert SN, Levy MZ, Brisson D. Recent and rapid population growth and range expansion of the Lyme disease tick vector, Ixodes scapularis, in North America. Evolution. 2015;69:1678–89.
Kugeler KJ, Farley GM, Forrester JD, Mead PS. Geographic distribution and expansion of human lyme disease, United States. Emerg Infect Dis. 2015;21:1455–7.
Clow KM, Leighton PA, Ogden NH, Lindsay LR, Michel P, Pearl DL, Jardine CM. Northward range expansion of Ixodes scapularis evident over a short timescale in Ontario, Canada. PLoS ONE. 2017;12:e0189393.
Schulze TL, Jordan RA, Schulze CJ, Healy SP, Jahn MB, Piesman J. Integrated use of 4-Poster passive topical treatment devices for deer, targeted acaricide applications, and Maxforce TMS bait boxes to rapidly suppress populations of Ixodes scapularis (Acari: Ixodidae) in a residential landscape. J Med Entomol. 2007;44:830–9.
Williams SC, Little EAH, Stafford KC, Molaei G, Linske MA. Integrated control of juvenile Ixodes scapularis parasitizing Peromyscus leucopus in residential settings in Connecticut, United States. Ticks Tick-Borne Dis. 2018;9:1310–6.
Stafford KC, Kitron U. Environmental management for Lyme borreliosis control. In: Gray J, Kahl O, Lane RS, Stanek G, editors. Lyme borreliosis: biology, epidemiology and control. CABI; 2002. p. 301–34.
Eisen L. Rodent-targeted approaches to reduce acarological risk of human exposure to pathogen-infected Ixodes ticks. Ticks Tick-borne Dis. 2023;14:102119.
Jordan RA, Schulze TL. Availability and nature of commercial tick control services in three Lyme disease endemic states. J Med Entomol. 2020;57:807–14.
Eisen L, Dolan MC. Evidence for personal protective measures to reduce human contact with blacklegged ticks and for environmentally based control methods to suppress host-seeking blacklegged ticks and reduce infection with Lyme disease spirochetes in tick vectors and rodent reservoirs. J Med Entomol. 2016;53:1063–92.
Ostfeld RS, Adish S, Mowry S, Bremer W, Duerr S, Evans AS, Fischhoff IR, Keating F, Pendleton J, Pfister A, Teator M, Keesing F. Effects of neighborhood-scale acaricidal treatments on infection prevalence of blacklegged ticks (Ixodes scapularis) with three zoonotic pathogens. Pathogens. 2023;12:172.
Gomes-Solecki M. Blocking pathogen transmission at the source: reservoir targeted OspA-based vaccines against Borrelia burgdorferi. Front Cell Infect Microbiol. 2014;4:136.
Eisen L, Stafford KC III. Barriers to effective tick management and tick-bite prevention in the United States (Acari: Ixodidae). J Med Entomol. 2021;58:1588–600.
Paules CI, Marston HD, Bloom ME, Fauci AS. Tickborne diseases—Confronting a growing threat. N Engl J Med. 2018;379:701–3.
Kugeler KJ, Jordan RA, Schulze TL, Griffith KS, Mead PS. Will culling white-tailed deer prevent Lyme disease? Zoonoses Public Health. 2016;63:337–45.
Telford SR III. Deer reduction is a cornerstone of integrated deer tick management. J Integr Pest Manag. 2017;8:25.
Bloemer SR, Mount GA, Morris TA, Zimmerman RH, Barnard DR, Snoddy EL. Management of lone star ticks (Acari: Ixodidae) in recreational areas with acaricide applications, vegetative management, and exclusion of white-tailed deer. J Med Entomol. 1990;27:543–50.
Ginsberg HS, Stafford III KC. Management of ticks and tick-borne diseases. In: Tick-Borne Diseases of Humans [Internet]. John Wiley & Sons, Ltd; 2005 [cited 2023 Oct 20]. p. 65–86. Available from: https://doi.org/10.1128/9781555816490.ch4
Kilpatrick HJ, LaBonte AM, Stafford KC. The relationship between deer density, tick abundance, and human cases of Lyme disease in a residential community. J Med Entomol. 2014;51:777–84.
Stafford KC III. Tick management handbook. Conn Agric Exp Stn Bull. 2007;1010:79.
Williams SC, Stafford KC, Molaei G, Linske MA. Integrated control of nymphal ixodes scapularis: effectiveness of white-tailed deer reduction, the entomopathogenic fungus metarhizium anisopliae, and fipronil-based rodent bait boxes.Vector-Borne Zoonotic Dis. 2018;18:55–64.
Eisen L. Control of ixodid ticks and prevention of tick-borne diseases in the United States: the prospect of a new Lyme disease vaccine and the continuing problem with tick exposure on residential properties. Ticks Tick Borne Dis. 2021;12:101649.
Ogden NH, Bigras-Poulin M, O’Callaghan CJ, Barker IK, Lindsay LR, Maarouf A, Smoyer-Tomic KE, Waltner-Toews D, Charron D. A dynamic population model to investigate effects of climate on geographic range and seasonality of the tick Ixodes scapularis. Int J Parasitol. 2005;35:375–89.
Ogden NH, Bigras-Poulin M, O’callaghan CJ, Barker IK, Kurtenbach K, Lindsay LR, Charron DF. Vector seasonality, host infection dynamics and fitness of pathogens transmitted by the tick Ixodes scapularis. Parasitology. 2007;134:209–27.
Diuk-Wasser MA, VanAcker MC, Fernandez MP. Impact of land use changes and habitat fragmentation on the eco-epidemiology of tick-borne diseases. J Med Entomol. 2021;58:1546–64.
Bender EA, Case TJ, Gilpin ME. Perturbation experiments in community ecology: theory and practice. Ecology. 1984;65:1–13.
Pihlatie MK, Kiese R, Brüggemann N, Butterbach-Bahl K, Kieloaho AJ, Laurila T, Lohila A, Mammarella I, Minkkinen K, Penttilä T, Schönborn J, Vesala T. Greenhouse gas fluxes in a drained peatland forest during spring frost-thaw event. Biogeosciences. 2010;7:1715–27.
Keppel G, Anderson S, Williams C, Kleindorfer S, O’Connell C. Microhabitats and canopy cover moderate high summer temperatures in a fragmented Mediterranean landscape. PLoS ONE. 2017;12:e0183106.
Tornqvist L, Vartia P, Vartia YO. How should relative changes be measured? Am Stat. 1985;39:43–6.
Hyndman RJ, Koehler AB. Another look at measures of forecast accuracy. Int J Forecast. 2006;22:679–88.
Gupta HV, Sorooshian S, Yapo PO. Status of automatic calibration for hydrologic models: comparison with multilevel expert calibration. J Hydrol Eng. 1999;4:135–43.
Nash JE, Sutcliffe JV. River flow forecasting through conceptual models part I—A discussion of principles. J Hydrol. 1970;10:282–90.
Sullivan CF, Parker BL, Skinner M. A review of commercial metarhizium- and beauveria-based biopesticides for the biological control of ticks in the USA. Insects. 2022;13:260.
Martin AM, Buttke D, Raphael J, Taylor K, Maes S, Parise CM, Ginsberg HS, Cross PC. Deer management generally reduces densities of nymphal Ixodes scapularis, but not prevalence of infection with Borrelia burgdorferi sensu stricto. Ticks Tick-Borne Dis. 2023;14:102202.
Ostfeld RS, Mowry S, Bremer W, Duerr S, Evans AS, Fischhoff IR, Hinckley AF, Hook SA, Keating F, Pendleton J, Pfister A, Teator M, Keesing F. Impacts over time of neighborhood-scale interventions to control ticks and tick-borne disease incidence. Vector Borne Zoonotic Dis. 2023;23:89–105.
Tsao JI, Wootton JT, Bunikis J, Luna MG, Fish D, Barbour AG. An ecological approach to preventing human infection: vaccinating wild mouse reservoirs intervenes in the Lyme disease cycle. Proc Natl Acad Sci. 2004;101:18159–64.
Carrera-Pineyro D, Hanes H, Litzler A, McCormack A, Velazquez-Molina J, Mubayi A, Ríos-Soto K, Kribs C. Cost analysis of vaccination in tick-mouse transmission of Lyme disease. J Theor Biol. 2020;494:110245.
Vesga JF, Clark MHA, Ayazi E, Apolloni A, Leslie T, Edmunds WJ, Métras R. Transmission dynamics and vaccination strategies for Crimean-Congo haemorrhagic fever virus in Afghanistan: a modelling study. PLoS Negl Trop Dis. 2022;16:e0010454.
Jourdain F, de Valk H, Noël H, Paty MC, L’Ambert G, Franke F, Mouly D, Desenclos JC, Roche B. Estimating chikungunya virus transmission parameters and vector control effectiveness highlights key factors to mitigate arboviral disease outbreaks. PLoS Negl Trop Dis. 2022;16:e0010244.
Routledge I, Walker M, Cheke RA, Bhatt S, Nkot PB, Matthews GA, Baleguel D, Dobson HM, Wiles TL, Basañez MG. Modelling the impact of larviciding on the population dynamics and biting rates of Simulium damnosum (s.l.): implications for vector control as a complementary strategy for onchocerciasis elimination in Africa. Parasit Vectors. 2018;11:316.
Barbu C, Dumonteil E, Gourbière S. Optimization of control strategies for non-domiciliated triatoma dimidiata, chagas disease vector in the Yucatán Peninsula, Mexico. PLoS Negl Trop Dis. 2009;3:e416.
Keesing F, Mowry S, Bremer W, Duerr S, Evans AS, Fischhoff IR, Hinckley AF, Hook SA, Keating F, Pendleton J, Pfister A, Teator M, Ostfeld RS. Effects of tick-control interventions on tick abundance, human encounters with ticks, and incidence of tickborne diseases in residential neighborhoods, New York, USA. Emerg Infect Dis. 2022;28:957–66.
Gandy S, Kilbride E, Biek R, Millins C, Gilbert L. No net effect of host density on tick-borne disease hazard due to opposing roles of vector amplification and pathogen dilution. Ecol Evol. 2022;12:e9253.
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.
Rand PW, Lubelczyk C, Holman MS, Lacombe EH, Smith RP. Abundance of Ixodes scapularis (Acari: Ixodidae) after the complete removal of deer from an Isolated Offshore Island, endemic for lyme disease. J Med Entomol. 2004;41:779–84.
Gaff H, Eisen RJ, Eisen L, Nadolny R, Bjork J, Monaghan AJ. LYMESIM 2.0: an updated simulation of blacklegged tick (Acari: Ixodidae) population dynamics and enzootic transmission of Borrelia burgdorferi (Spirochaetales: Spirochaetaceae). J Med Entomol. 2020;57:715–27.
Li S, Vanwambeke SO, Licoppe AM, Speybroeck N. Impacts of deer management practices on the spatial dynamics of the tick Ixodes ricinus: a scenario analysis. Ecol Model. 2014;276:1–13.
Wang HH, Grant WE, Teel PD. Simulation of climate–host–parasite–landscape interactions: a spatially explicit model for ticks (Acari: Ixodidae). Ecol Model. 2012;243:42–62.
Wang HH, Teel PD, Grant WE, Schuster G, Pérez de León AA. Simulated interactions of white-tailed deer (Odocoileus virginianus), climate variation and habitat heterogeneity on southern cattle tick (Rhipicephalus (Boophilus) microplus) eradication methods in south Texas, USA. Ecol Model. 2016;342:82–96.
Schulze TL, Eisen L, Russell K, Jordan RA. Community-based integrated tick management programs: cost and feasibility scenarios. J Med Entomol. 2023;60:1048–60.
Acknowledgement
This publication was supported by Department of Defense Army Research Office (ARO) award W911QY1910006. It was also partially supported by the Cooperative Agreement Number U01CK000509-01 between the Centers for Disease Control and Prevention and Northeast Regional Center for Excellence in Vector Borne Diseases, and the National Science Foundation’s Coupled Natural Human Systems 2/Dynamics of Integrated Socio-Environmental Systems (CNH2/DISES) program (award #1924061). Its contents are solely the responsibility of the authors and do not necessarily represent the official views of the Department of Defense, the Centers for Disease Control and Prevention, the Department of Health and Human Services, or the National Science Foundation.
Funding
This publication was supported by Department of Defense Army Research Office (ARO) award W911QY1910006. It was also partially supported by the Cooperative Agreement Number U01CK000509-01 between the Centers for Disease Control and Prevention and Northeast Regional Center for Excellence in Vector Borne Diseases, and the National Science Foundation’s Coupled Natural Human Systems 2/Dynamics of Integrated Socio-Environmental Systems (CNH2/DISES) program (award #1924061).
Author information
Authors and Affiliations
Contributions
M.D.W. conceptualized and supervised the project. S.C.W. provided the data from residential properties. D.R.C. and J.B. developed the code for model analysis and drafted the manuscript. All authors revised and approved the latest version of this article.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Verbal permission was obtained from homeowners to access their properties. Internal Review Board review considered the work exempt because there are no data derived from humans.
Competing interests
The authors declare no competing interests.
Consent for publication
Not Applicable.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
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
Ruiz-Carrascal, D., Bastard, J., Williams, S.C. et al. Modeling platform to assess the effectiveness of single and integrated Ixodes scapularis tick control methods. Parasites Vectors 17, 339 (2024). https://doi.org/10.1186/s13071-024-06387-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13071-024-06387-2