Context-dependent host dispersal and habitat fragmentation determine heterogeneity in infected tick burdens: an agent-based modelling study

Borrelia burgdorferi sensu lato prevalence in Ixodes scapularis from Canada: A thirty-year summary and meta-analysis (1990-2020)

Borrelia burgdorferi sensu lato (Bb) are a complex of bacteria genospecies that can cause Lyme disease (LD) in humans after the bite of an infected Ixodes spp. vector tick. In Canada, incidence of LD is increasing in part due to the rapid geographic expansion of Ixodes scapularis across the southcentral and eastern provinces. To better understand temporal and spatial (provincial) prevalence of Bb infection of I. scapularis and how tick surveillance is utilized in Canada to assess LD risk, a literature review was conducted. Tick surveillance studies published between January 1975 to November 2023, that measured the prevalence of Bb in I. scapularis via "passive surveillance" from the public citizenry or "active surveillance" by drag or flag sampling of host-seeking ticks in Canada were included for review. Meta-analyses were conducted via random effects modeling. Forty-seven articles, yielding 26 passive and 28 active surveillance studies, met inclusion criteria. Mean durations of collection for I. scapularis were 2.1 years in active surveillance studies (1999-2020) and 5.5 years by passive surveillance studies (1990-2020). Collectively, data were extracted on 99,528 I. scapularis nymphs and adults collected between 1990-2020 across nine provinces, including Newfoundland & Labrador (33 ticks) and Alberta (208 ticks). More studies were conducted in Ontario (36) than any other province. Across nine provinces, the prevalence of Bb infection in I. scapularis collected by passive surveillance was 14.6% with the highest prevalence in Nova Scotia at 20.5% (minimum studies >1). Among host-seeking I. scapularis collected via active surveillance, Bb infection prevalence was 10.5% in nymphs, 31.9% in adults, and 23.8% across both life stages. Host-seeking I. scapularis nymphs and adults from Ontario had the highest Bb prevalence at 13.6% and 34.8%, respectively. Between 2007-2019, Bb infection prevalence in host-seeking I. scapularis was positively associated over time (p<0.001) which is concurrent with a ∼25-fold increase in the number of annually reported LD cases in Canada over the same period. The prevalence of Bb-infection in I. scapularis has rapidly increased over three decades as reported by tick surveillance studies in Canada which coincides with increasing human incidence for LD. The wide-ranging distribution and variable prevalence of Bb-infected I. scapularis ticks across provinces demonstrates the growing need for long-term standardized tick surveillance to monitor the changing trends in I. scapularis populations and best define LD risk areas in Canada.

IxPopDyMod: an R package to write, run, and analyze tick population and infection dynamics models

Given the increasing prevalence of tick-borne diseases, such as Lyme disease, modeling the population and infection dynamics of tick vectors is an important public health tool. These models have applications for testing the effects of control methods or climate change on tick populations. There is an established history of tick population models, but code for them is rarely shared, especially not in a convenient format for others to modify and use. We present an R package, called IxPopDyMod, intended to function as a flexible and consistent framework for reproducible Ixodidae (hard-bodied ticks) population dynamics models. Here we focus on two key parts of the package: a function to create valid model configurations and a function to run a configured model and return the daily population over time. We provide three examples in appendices: one reproducing an existing Ixodes scapularis population model, one providing a novel Dermacentor albipictus model, and one showing Borrelia burgdorferi infection in ticks. Together these examples show the flexibility of the package to model scenarios of interest to tick researches.

Mechanistic movement models to predict geographic range expansions of ticks and tick-borne pathogens: Case studies with Ixodes scapularis and Amblyomma americanum in eastern North America

The geographic range of the blacklegged tick, Ixodes scapularis, is expanding northward from the United States into southern Canada, and studies suggest that the lone star tick, Amblyomma americanum, will follow suit. These tick species are vectors for many zoonotic pathogens, and their northward range expansion presents a serious threat to public health. Climate change (particularly increasing temperature) has been identified as an important driver permitting northward range expansion of blacklegged ticks, but the impacts of host movement, which is essential to tick dispersal into new climatically suitable regions, have received limited investigation. Here, a mechanistic movement model was applied to landscapes of eastern North America to explore 1) relationships between multiple ecological drivers and the speed of the northward invasion of blacklegged ticks infected with the causative agent of Lyme disease, Borrelia burgdorferi sensu stricto, and 2) its capacity to simulate the northward range expansion of infected blacklegged ticks and uninfected lone star ticks under theoretical scenarios of increasing temperature. Our results suggest that the attraction of migratory birds (long-distance tick dispersal hosts) to resource-rich areas during their spring migration and the mate-finding Allee effect in tick population dynamics are key drivers for the spread of infected blacklegged ticks. The modeled increases in temperature extended the climatically suitable areas of Canada for infected blacklegged ticks and uninfected lone star ticks towards higher latitudes by up to 31% and 1%, respectively, and with an average predicted speed of the range expansion reaching 61 km/year and 23 km/year, respectively. Differences in the projected spatial distribution patterns of these tick species were due to differences in climate envelopes of tick populations, as well as the availability and attractiveness of suitable habitats for migratory birds. Our results indicate that the northward invasion process of lone star ticks is primarily driven by local dispersal of resident terrestrial hosts, whereas that of blacklegged ticks is governed by long-distance migratory bird dispersal. The results also suggest that mechanistic movement models provide a powerful approach for predicting tick-borne disease risk patterns under complex scenarios of climate, socioeconomic and land use/land cover changes.

Context-dependent host dispersal and habitat fragmentation determine heterogeneity in infected tick burdens: an agent-based modelling study

As the incidence of tick-borne diseases has sharply increased over the past decade, with serious consequences for human and animal health, there is a need to identify ecological drivers contributing to heterogeneity in tick-borne disease risk. In particular, the relative importance of animal host dispersal behaviour in its three context-dependent phases of emigration, transfer and settlement is relatively unexplored. We built a spatially explicit agent-based model to investigate how the host dispersal process, in concert with the tick and host demographic processes, habitat fragmentation and the pathogen transmission process, affects infected tick distributions among hosts. A sensitivity analysis explored the impacts of different input parameters on infected tick burdens on hosts and infected tick distributions among hosts. Our simulations indicate that ecological predictors of infected tick burdens differed among the post-egg life stages of ticks, with tick attachment and detachment, tick questing activity and pathogen transmission dynamics identified as key processes, in a coherent way. We also found that the type of host settlement strategy and the proportion of habitat suitable for hosts determined super-spreading of infected ticks. We developed a theoretical mechanistic framework that can serve as a first step towards applied studies of on-the-ground public health intervention strategies.