Range Expansion of Ixodes scapularis and Borrelia burgdorferi in Ontario, Canada, from 2017 to 2019

Comparative ecological analysis and predictive modeling of tick-borne pathogens

Tick-borne diseases constitute the predominant vector-borne health threat in North America. Recent observations have noted a significant expansion in the range of the black-legged tick (Ixodes scapularis Say, Acari: Ixodidae), alongside a rise in the incidence of diseases caused by its transmitted pathogens: Borrelia burgdorferi Johnson (Spirochaetales: Spirochaetaceae), Babesia microti Starcovici (Piroplasmida: Babesiidae), and Anaplasma phagocytophilium Zhu (Rickettsiales: Anaplasmataceae), the causative agents of Lyme disease, babesiosis, and anaplasmosis, respectively. Prior research identified environmental features that influence the ecological dynamics of I. scapularis and B. burgdorferi that can be used to predict the distribution and abundance of these organisms, and thus Lyme disease risk. In contrast, there is a paucity of research into the environmental determinants of B. microti and A. phagocytophilium. Here, we use over a decade of surveillance data to model the impact of environmental features on the infection prevalence of these increasingly common human pathogens in ticks across New York State (NYS). Our findings reveal a consistent northward and westward expansion of B. microti in NYS from 2009 to 2019, while the range of A. phagocytophilum varied at fine spatial scales. We constructed biogeographic models using data from over 650 site-year visits and encompassing more than 250 environmental variables to accurately forecast infection prevalence for each pathogen to a future year that was not included in model training. Several environmental features were identified to have divergent effects on the pathogens, revealing potential ecological differences governing their distribution and abundance. These validated biogeographic models have applicability for disease prevention efforts.

Broad-scale ecological niches of pathogens vectored by the ticks Ixodes scapularis and Amblyomma americanum in North America

Environmental dimensions, such as temperature, precipitation, humidity, and vegetation type, influence the activity, survival, and geographic distribution of tick species. Ticks are vectors of various pathogens that cause disease in humans, and Ixodes scapularis and Amblyomma americanum are among the tick species that transmit pathogens to humans across the central and eastern United States. Although their potential geographic distributions have been assessed broadly via ecological niche modeling, no comprehensive study has compared ecological niche signals between ticks and tick-borne pathogens. We took advantage of National Ecological Observatory Network (NEON) data for these two tick species and associated bacteria pathogens across North America. We used two novel statistical tests that consider sampling and absence data explicitly to perform these explorations: a univariate analysis based on randomization and resampling, and a permutational multivariate analysis of variance. Based on univariate analyses, in Amblyomma americanum, three pathogens (Borrelia lonestari, Ehrlichia chaffeensis, and E. ewingii) were tested; pathogens showed nonrandom distribution in at least one environmental dimension. Based on the PERMANOVA test, the null hypothesis that the environmental position and variation of pathogen-positive samples are equivalent to those of A. americanum could not be rejected for any of the pathogens, except for the pathogen E. ewingii in maximum and minimum vapor pressure and minimum temperature. For Ixodes scapularis, six pathogens (A. phagocytophilum, Babesia microti, Borrelia burgdorferi sensu lato, B. mayonii, B. miyamotoi, and Ehrlichia muris-like) were tested; only B. miyamotoi was not distinct from null expectations in all environmental dimensions, based on univariate tests. In the PERMANOVA analyses, the pathogens departed from null expectations for B. microti and B. burgdorferi sensu lato, with smaller niches in B. microti, and larger niches in B. burgdorferi sensu lato, than the vector. More generally, this study shows the value of large-scale data resources with consistent sampling methods, and known absences of key pathogens in particular samples, for answering public health questions, such as the relationship of presence and absence of pathogens in their hosts respect to environmental conditions.

Ixodes scapularis density and Borrelia burgdorferi prevalence along a residential-woodland gradient in a region of emerging Lyme disease risk

The environmental risk of Lyme disease, defined by the density of Ixodes scapularis ticks and their prevalence of Borrelia burgdorferi infection, is increasing across the Ottawa, Ontario region, making this a unique location to explore the factors associated with environmental risk along a residential-woodland gradient. In this study, we collected I. scapularis ticks and trapped Peromyscus spp. mice, tested both for tick-borne pathogens, and monitored the intensity of foraging activity by deer in residential, woodland, and residential-woodland interface zones of four neighbourhoods. We constructed mixed-effect models to test for site-specific characteristics associated with densities of questing nymphal and adult ticks and the infection prevalence of nymphal and adult ticks. Compared to residential zones, we found a strong increasing gradient in tick density from interface to woodland zones, with 4 and 15 times as many nymphal ticks, respectively. Infection prevalence of nymphs and adults together was 15 to 24 times greater in non-residential zone habitats. Ecological site characteristics, including soil moisture, leaf litter depth, and understory density, were associated with variations in nymphal density and their infection prevalence. Our results suggest that high environmental risk bordering residential areas poses a concern for human-tick encounters, highlighting the need for targeted disease prevention.

Epidemiology and clinical manifestations of reported Lyme disease cases: Data from the Canadian Lyme disease enhanced surveillance system

Lyme disease cases reported in seven Canadian provinces from 2009 to 2019 through the Lyme Disease Enhanced Surveillance System are described herein by demographic, geography, time and season. The proportion of males was greater than females. Bimodal peaks in incidence were observed in children and older adults (≥60 years of age) for all clinical signs except cardiac manifestations, which were more evenly distributed across age groups. Proportions of disease stages varied between provinces: Atlantic provinces reported mainly early Lyme disease, while Ontario reported equal proportions of early and late-stage Lyme disease. Early Lyme disease cases were mainly reported between May through November, whereas late Lyme disease were reported in December through April. Increased awareness over time may have contributed to a decrease in the proportion of cases reporting late disseminated Lyme disease. These analyses help better describe clinical features of reported Lyme disease cases in Canada.

Isolation and multilocus sequence typing of Borrelia burgdorferi from Ixodes scapularis collected from dogs in Ontario, Canada

OBJECTIVE

To identify the multilocus sequence typing (MLST) sequence types of Borrelia burgdorferi from Ixodes scapularis in Ontario, Canada.

RESULTS

One hundred and eighty-five I. scapularis ticks were submitted from 134 dogs via participating clinics from April 1, 2019, to March 31, 2020. Seventeen MLST sequence types of B. burgdorferi were detected from fifty-eight cultured isolates from 21 ticks. The most common MLST sequence types were 12 and 16. Mixed infections of two MLST sequence types were detected in four ticks. Three sequence types (48, 317, 639) were new detections in Ontario.

Comparing Canadian Lyme disease risk area classification methodologies

Lyme disease risk areas have increased across Canada in recent decades with the ongoing range expansion of Ixodes scapularis and Borrelia burgdorferi. Different methodologies are used by federal and provincial governments to determine local Lyme disease risk, which can make comparisons between regions challenging. In this study, seven Canadian Lyme disease risk classification methodologies were compared with each other to highlight the strengths and limitations of how each definition measured I. scapularis and B. burgdorferi risk. Each methodology was applied to active surveillance data from Ontario, and per cent agreement and kappa statistics were calculated. The methodologies varied in their measurements of the risk of exposure to I. scapularis and B. burgdorferi based on their use of active surveillance techniques, multiple types of collected surveillance data and laboratory confirmation of B. burgdorferi. Most initial Lyme disease risk site classifications were maintained over time. Kappa and per cent agreement statistics highlighted large differences between 8 of the 15 methodology pairings, indicating the presence of inconsistencies between most methodologies. Accurate, consistent surveillance and assessment of the spread of I. scapularis and its pathogens will aid with communicating Lyme disease risk to the public and preventing tick-borne pathogen transmission.

Emerging Tick-Borne Pathogens in Central Canada: Recent Detections of Babesia odocoilei and Rickettsia rickettsii

Background: The spread of emerging tick-borne pathogens has steadily increased in Canada with the widespread establishment of tick vectors and vertebrate hosts. At present, Borrelia burgdorferi, the bacterium causing Lyme disease, is the most common tick-borne pathogen in Canada and primarily transmitted by Ixodes scapularis. A low prevalence of other emerging tick-borne pathogens, such as Anaplasma phagocytophilum, Babesia species, Borrelia miyamotoi, and Francisella tularensis have also been detected through surveillance efforts in Canada. Although Rickettsia rickettsii has been historically detected in Haemaphysalis leporispalustris in Canada, the current prevalence and geographic extent of this pathogen is unknown. Material and Methods: In this study, we assessed the presence and prevalence of several emerging tick-borne pathogens in ticks and hosts collected through tick dragging and small mammal trapping in Central Canada. Results: Nested PCR testing detected three pathogen species in ticks, with Babesia odocoilei and B. burgdorferi in I. scapularis in addition to R. rickettsii in H. leporispalustris. Three pathogen species were detected in small mammals by nested PCR including B. odocoilei in Blarina brevicauda, Babesia microti in Peromyscus leucopus, and a Hepatozoon species in P. leucopus and Peromyscus maniculatus. B. burgdorferi and Babesia species were the pathogens most often detected in our samples, suggesting they are widely distributed across Central Canada. We also detected B. odocoilei and R. rickettsii beyond their known geographic distribution. Conclusions: Our results provide evidence that emerging tick-borne pathogens may be present outside defined risk areas identified by current surveillance efforts in Canada. As a result, emerging tick-borne pathogens introduced by the dispersal of infected ticks by migratory birds or maintained by hosts and vectors through cryptic transmission cycles may go undetected. More comprehensive testing including all tick life stages and additional tick-borne pathogens will help detect the spread and potential risk of emerging or re-emerging tick-borne pathogens for human and wildlife populations throughout Canada.