High-Resolution Ecological Niche Modeling of Ixodes scapularis Ticks Based on Passive Surveillance Data at the Northern Frontier of Lyme Disease Emergence in North America

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.

Refining Ixodes scapularis (Acari: Ixodidae) distribution models: a comparison of current methods to an established protocol

Blacklegged ticks (Ixodes scapularis Say) pose an enormous public health risk in eastern North America as the vector responsible for transmitting 7 human pathogens, including those causing the most common vector-borne disease in the United States, Lyme disease. Species distribution modeling is an increasingly popular method for predicting the potential distribution and subsequent risk of blacklegged ticks, however, the development of such models thus far is highly variable and would benefit from the use of standardized protocols. To identify where standardized protocols would most benefit current distribution models, we completed the "Overview, Data, Model, Assessment, and Prediction" (ODMAP) distribution modeling protocol for 21 publications reporting 22 blacklegged tick distribution models. We calculated an average adherence of 73.4% (SD ± 29%). Most prominently, we found that authors could better justify and connect their selection of variables and associated spatial scales to blacklegged tick ecology. In addition, the authors could provide clearer descriptions of model development, including checks for multicollinearity, spatial autocorrelation, and plausibility. Finally, authors could improve their reporting of variable effects to avoid undermining the models' utility in informing species-environment relationships. To enhance future model rigor and reproducibility, we recommend utilizing several resources including the ODMAP protocol, and suggest that journals make protocol compliance a publication prerequisite.

Models and data used to predict the abundance and distribution of Ixodes scapularis (blacklegged tick) in North America: a scoping review

Tick-borne diseases (TBD) remain prevalent worldwide, and risk assessment of tick habitat suitability is crucial to prevent or reduce their burden. This scoping review provides a comprehensive survey of models and data used to predict I. scapularis distribution and abundance in North America. We identified 4661 relevant primary research articles published in English between January 1st, 2012, and July 18th, 2022, and selected 41 articles following full-text review. Models used data-driven and mechanistic modelling frameworks informed by diverse tick, hydroclimatic, and ecological variables. Predictions captured tick abundance (n = 14, 34.1%), distribution (n = 22, 53.6%) and both (n = 5, 12.1%). All studies used tick data, and many incorporated both hydroclimatic and ecological variables. Minimal host- and human-specific data were utilized. Biases related to data collection, protocols, and tick data quality affect completeness and representativeness of prediction models. Further research and collaboration are needed to improve prediction accuracy and develop effective strategies to reduce TBD.

Resource selection by New York City deer reveals the effective interface between wildlife, zoonotic hazards and humans

Although the role of host movement in shaping infectious disease dynamics is widely acknowledged, methodological separation between animal movement and disease ecology has prevented researchers from leveraging empirical insights from movement data to advance landscape scale understanding of infectious disease risk. To address this knowledge gap, we examine how movement behaviour and resource utilization by white-tailed deer (Odocoileus virginianus) determines blacklegged tick (Ixodes scapularis) distribution, which depend on deer for dispersal in a highly fragmented New York City borough. Multi-scale hierarchical resource selection analysis and movement modelling provide insight into how deer's movements contribute to the risk landscape for human exposure to the Lyme disease vector-I. scapularis. We find deer select highly vegetated and accessible residential properties which support blacklegged tick survival. We conclude the distribution of tick-borne disease risk results from the individual resource selection by deer across spatial scales in response to habitat fragmentation and anthropogenic disturbances.

Comparison of acarological risk metrics derived from active and passive surveillance and their concordance with tick-borne disease incidence

Tick-borne diseases continue to threaten human health across the United States. Both active and passive tick surveillance can complement human case surveillance, providing spatio-temporal information on when and where humans are at risk for encounters with ticks and tick-borne pathogens. However, little work has been done to assess the concordance of the acarological risk metrics from each surveillance method. We used data on Ixodes scapularis and its associated human pathogens from Connecticut (2019-2021) collected through active collections (drag sampling) or passive submissions from the public to compare county estimates of tick and pathogen presence, infection prevalence, and tick abundance by life stage. Between the surveillance strategies, we found complete agreement in estimates of tick and pathogen presence, high concordance in infection prevalence estimates for Anaplasma phagocytophilum, Borrelia burgdorferi sensu stricto, and Babesia microti, but no consistent relationships between actively and passively derived estimates of tick abundance or abundance of infected ticks by life stage. We also compared nymphal metrics (i.e., pathogen prevalence in nymphs, nymphal abundance, and abundance of infected nymphs) with reported incidence of Lyme disease, anaplasmosis, and babesiosis, but did not find any consistent relationships with any of these metrics. The small spatial and temporal scale for which we had consistently collected active and passive data limited our ability to find significant relationships. Findings are likely to differ if examined across a broader spatial or temporal coverage with greater variation in acarological and epidemiological outcomes. Our results indicate similar outcomes between some actively and passively derived tick surveillance metrics (tick and pathogen presence, pathogen prevalence), but comparisons were variable for abundance estimates.

Passive Surveillance of Human-Biting Ixodes scapularis Ticks in Massachusetts from 2015-2019

This study aimed to analyze human-biting Ixodes scapularis ticks submitted to TickReport tick testing service from 2015-2019 in Massachusetts to (1) examine possible patterns of pathogen-positive adult and nymphal ticks over time and (2) explore how socioeconomic factors can influence tick submissions. A passive surveillance data set of ticks and tick-borne pathogens was conducted over 5 years (2015-2019) in Massachusetts. The percentages of four tick-borne pathogens: Borrelia burgdorferi, Anaplasma phagocytophilum, Babesia microti, and Borrelia miyamotoi were determined by Massachusetts county and by month and year. Regression models were used to examine the association between zip-code-level socioeconomic factors and submissions. A total of 13,598 I. scapularis ticks were submitted to TickReport from Massachusetts residents. The infection rate of B. burgdorferi, A. phagocytophilum, and B. microti was 39%, 8%, and 7% in adult ticks; 23%, 6%, and 5% in nymphal ticks, respectively. A relatively higher level of education was associated with high tick submission. Passive surveillance of human-biting ticks and associated pathogens is important for monitoring tick-borne diseases, detecting areas with potentially high risks, and providing public information. Socioeconomic factors should be considered to produce more generalizable passive surveillance data and to target potentially underserved areas.

Rurality, Socioeconomic Status, and Residence in Environmental Risk Areas Associated with Increased Lyme Disease Incidence in Ontario, Canada: A Case-Control Study

Background: Lyme disease (LD) is the most common tick-borne illness in North America. LD is acquired through exposure to the tick vector, Ixodes scapularis, known as the blacklegged tick. In Canada, LD is rapidly emerging, with the establishment of I. scapularis in many newly endemic regions posing a growing risk to local communities. In the Canadian context, many environmental and socioeconomic risk factors for human LD infection are yet to be ascertained and the degree of risk associated with residential and community exposure to ticks is not well known. Methods: We conducted a matched case-control study in southeastern Ontario, using LD patient data from provincial laboratory databases and uninfected population controls from 2014 to 2018. We aimed to identify area-level risk factors for LD and associations with residence in environmental risk areas, defined as areas with high model-predicted probability of I. scapularis occurrence, using the neighborhood dissemination area as the unit of analysis. Results: Using multivariable conditional logistic regression analysis, we identified that patients with LD had higher odds (odds ratio, OR; 95% confidence interval, CI) of living in neighborhoods with high probability of tick occurrence in the environment (OR = 2.2; 95% CI: 2.0-2.5), low walkability (OR = 1.6; 95% CI: 1.2-2.1), low material deprivation (OR = 1.4; 95% CI: 1.2-1.7), and low ethnic concentration (OR = 8.1; 95% CI: 6.7-9.9). We also found that the odds of LD infection for individuals residing in environmental risk areas was highest for those living in public health units (PHUs) with <250,000 population (OR = 3.0; 95% CI: 2.4-3.9) compared to those living in PHUs with >1,000,000 population (OR = 1.5; 95% CI: 1.1-2.1). Conclusion: This study shows that odds of human LD infection in Ontario, Canada is higher in less urbanized areas with higher socioeconomic status and indicates that exposure to ticks around the home residence or neighborhood is linked to increased odds of LD.

Surveillance of Rocky Mountain wood ticks (Dermacentor andersoni) and American dog ticks (Dermacentor variabilis) in Colorado

Ticks pose an emerging threat of infectious pathogen transmission in the United States in part due to expanding suitable habitat ranges in the wake of climate change. Active and passive tick surveillance can inform maps of tick distributions to warn the public of their risk of exposure to ticks. In Colorado, widespread active surveillance programs have difficulty due to the state's diverse terrain. However, combining multiple citizen science techniques can create a more accurate representation of tick distribution than any passive surveillance dataset alone. Our study uses county-level tick distribution data from Northern Arizona University, the Colorado Department of Public Health and the Environment, and veterinary surveillance in addition to literature data to assess the distribution of the Rocky Mountain wood tick, Dermacentor andersoni, and the American dog tick, Dermacentor variabilis. We found that D. andersoni for the most part inhabits counties at higher elevations than D. variabilis in Colorado.

Environmental Drivers of Gulf Coast Tick (Acari: Ixodidae) Range Expansion in the United States

In the United States, the Gulf Coast tick (Amblyomma maculatum Koch) is a species of growing medical and veterinary significance, serving as the primary vector of the pathogenic bacterium, Rickettsia parkeri (Rickettsiales: Rickettsiaceae), in humans and the apicomplexan parasite, Hepatozoon americanum, in canines. Ongoing reports of A. maculatum from locations outside its historically reported distribution in the southeastern United States suggest the possibility of current and continuing range expansion. Using an ecological niche modeling approach, we combined new occurrence records with high-resolution climate and land cover data to investigate environmental drivers of the current distribution of A. maculatum in the United States. We found that environmental suitability for A. maculatum varied regionally and was primarily driven by climatic factors such as annual temperature variation and seasonality of precipitation. We also found that presence of A. maculatum was associated with open habitat with minimal canopy cover. Our model predicts large areas beyond the current distribution of A. maculatum to be environmentally suitable, suggesting the possibility of future northward and westward range expansion. These predictions of environmental suitability may be used to identify areas at potential risk for establishment and to guide future surveillance of A. maculatum in the United States.

A Scoping Review of Species Distribution Modeling Methods for Tick Vectors

Background

Globally, tick-borne disease is a pervasive and worsening problem that impacts human and domestic animal health, livelihoods, and numerous economies. Species distribution models are useful tools to help address these issues, but many different modeling approaches and environmental data sources exist.

Objective

We conducted a scoping review that examined all available research employing species distribution models to predict occurrence and map tick species to understand the diversity of model strategies, environmental predictors, tick data sources, frequency of climate projects of tick ranges, and types of model validation methods.

Design

Following the PRISMA-ScR checklist, we searched scientific databases for eligible articles, their references, and explored related publications through a graphical tool (www.connectedpapers.com). Two independent reviewers performed article selection and characterization using a priori criteria.

Results

We describe data collected from 107 peer-reviewed articles that met our inclusion criteria. The literature reflects that tick species distributions have been modeled predominantly in North America and Europe and have mostly modeled the habitat suitability for Ixodes ricinus (n = 23; 21.5%). A wide range of bioclimatic databases and other environmental correlates were utilized among models, but the WorldClim database and its bioclimatic variables 1–19 appeared in 60 (56%) papers. The most frequently chosen modeling approach was MaxEnt, which also appeared in 60 (56%) of papers. Despite the importance of ensemble modeling to reduce bias, only 23 papers (21.5%) employed more than one algorithm, and just six (5.6%) used an ensemble approach that incorporated at least five different modeling methods for comparison. Area under the curve/receiver operating characteristic was the most frequently reported model validation method, utilized in nearly all (98.9%) included studies. Only 21% of papers used future climate scenarios to predict tick range expansion or contraction. Regardless of the representative concentration pathway, six of seven genera were expected to both expand and retract depending on location, while Ornithodoros was predicted to only expand beyond its current range.

Conclusion

Species distribution modeling techniques are useful and widely employed tools for predicting tick habitat suitability and range movement. However, the vast array of methods, data sources, and validation strategies within the SDM literature support the need for standardized protocols for species distribution and ecological niche modeling for tick vectors.

Effects of Climate on the Variation in Abundance of Three Tick Species in Illinois

The range of ticks in North America has been steadily increasing likely, in part, due to climate change. Along with it, there has been a rise in cases of tick-borne disease. Among those medically important tick species of particular concern are Ixodes scapularis Say (Acari: Ixodidae), Dermacentor variabilis Say (Acari: Ixodidae), and Amblyomma americanum Linneaus (Acari: Ixodidae). The aim of this study was to determine if climate factors explain existing differences in abundance of the three aforementioned tick species between two climatically different regions of Illinois (Central and Southern), and if climate variables impact each species differently. We used both zero-inflated regression approaches and Bayesian network analyses to assess relationships among environmental variables and tick abundance. Results suggested that the maximum average temperature and total precipitation are associated with differential impact on species abundance and that this difference varied by region. Results also reinforced a differential level of resistance to desiccation among these tick species. Our findings help to further define risk periods of tick exposure for the general public, and reinforce the importance of responding to each tick species differently.

Occupational Lyme Disease: A Systematic Review and Meta-Analysis

Lyme disease (LD) can have significant consequences for the health of workers. The frequency of infection can be estimated by using prevalence and incidence data on antibodies against Borrelia Burgdoferi (BB). A systematic search of studies published in English between 2002 and 2021 and a meta-analysis were conducted in PubMed/Medline, Web of Science, Scopus, and Google Scholar databases. Out of a total of 1125 studies retrieved, 35 articles were included in the systematic review. Overall, in these studies, outdoor workers showed a 20.5% BB seroprevalence rate. Meta-analysis, performed on 15 studies (3932 subjects), revealed a significantly increased risk in outdoor activities (OR 1.93 95%CI 1.15–3.23), with medium-level heterogeneity (I² = 69.2%), and non-significant publication bias. The estimated OR in forestry and agricultural workers was 2.36 (CI95% 1.28; 4.34) in comparison with the controls, while a non-significant increase in risk (OR = 1.05, CI95% 0.28; 3.88) was found in the remaining categories of workers (veterinarians, animal breeders, soldiers). The estimated pooled risk was significantly higher in the studies published until 2010 (OR 3.03 95%CI 1.39–6.61), while in more recent studies the odds became non-significant (OR 1.08 95% CI 0.63–1.85). The promotion of awareness campaigns targeting outdoor workers in endemic areas, and the implementation of local programs aimed at controlling range expansion of vectors, are key strategies for protecting workers.

A Geographic Information System Approach to Map Tick Exposure Risk at a Scale for Public Health Intervention

Tick-borne disease control and prevention have been largely ineffective compared to the control of other vector-borne diseases. Although control strategies exist, they are costly or ineffective at large spatial scales. We need tools to target these strategies to places of highest tick exposure risk. Here we present a geographic information system (GIS) method for mapping predicted tick exposure risk at a 200 m by 200 m resolution, appropriate for public health intervention. We followed the approach used to map tick habitat suitability over large areas. We used drag-cloth sampling to measure the density of nymphal blacklegged ticks (Ixodes scapularis, Say (Acari: Ixodidae)) at 24 sites in Addison and Rutland Counties, VT, United States. We used a GIS to average habitat, climatological, land-use/land-cover, and abiotic characteristics over 100 m, 400 m, 1,000 m, and 2,000 m buffers around each site to evaluate which characteristic at which buffer size best predicted density of nymphal ticks (DON). The relationships between predictor variables and DON were determined with random forest models. The 100 m buffer model performed best and explained 37.7% of the variation in DON, although was highly accurate at classifying sites as having below or above average DON. This model was applied to Addison County, VT, to predict tick exposure risk at a 200 m resolution. This GIS approach to map predicted DON over a small area with fine resolution, could be used to target public health campaigns and land management practices to reduce human exposure to ticks.

Projecting the Potential Distribution Areas of Ixodes scapularis (Acari: Ixodidae) Driven by Climate Change

Ixodes scapularis is a vector of tick-borne diseases. Climate change is frequently invoked as an important cause of geographic expansions of tick-borne diseases. Environmental variables such as temperature and precipitation have an important impact on the geographical distribution of disease vectors. We used the maximum entropy model to project the potential geographic distribution and future trends of I. scapularis. The main climatic variables affecting the distribution of potential suitable areas were screened by the jackknife method. Arc Map 10.5 was used to visualize the projection results to better present the distribution of potential suitable areas. Under climate change scenarios, the potential suitable area of I. scapularis is dynamically changing. The largest suitable area of I. scapularis is under SSP3-7.0 from 2081 to 2100, while the smallest is under SSP5-8.5 from 2081 to 2100, even smaller than the current suitable area. Precipitation in May and September are the main contributing factors affecting the potential suitable areas of I. scapularis. With the opportunity to spread to more potential suitable areas, it is critical to strengthen surveillance to prevent the possible invasion of I. scapularis.

Estimating disease vector population size from citizen science data

Citizen science projects have the potential to address hypotheses requiring extremely large datasets that cannot be collected with the financial and labour constraints of most scientific projects. Data collection by the general public could expand the scope of scientific enquiry if these data accurately capture the system under study. However, data collection inconsistencies by the untrained public may result in biased datasets that do not accurately represent the natural world. In this paper, we harness the availability of scientific and public datasets of the Lyme disease tick vector to identify and account for biases in citizen science tick collections. Estimates of tick abundance from the citizen science dataset correspond moderately with estimates from direct surveillance but exhibit consistent biases. These biases can be mitigated by including factors that may impact collector participation or effort in statistical models, which, in turn, result in more accurate estimates of tick population sizes. Accounting for collection biases within large-scale, public participation datasets could update species abundance maps and facilitate using the wealth of citizen science data to answer scientific questions at scales that are not feasible with traditional datasets.

Monitoring Risk: Tick and Borrelia burgdorferi Public Participatory Surveillance in the Canadian Maritimes, 2012-2020

Ticks are vectors of many diseases, including Lyme disease (Ld). Lyme disease is an emerging disease in Canada caused by infection with the Lyme borreliosis (Lb) members of the Borrelia genus of spirochaete bacteria, of which Borrelia&nbsp;burgdorferi is regionally the most prevalent. The primary tick vector in central and eastern Canada, Ixodes scapularis, is increasing in numbers and in the geographical extent of established populations. This study documents the distribution of ticks recovered by passive surveillance, and their B.&nbsp;burgdorferi infection prevalence, in three Canadian Maritime provinces from 2012-2020. These regions represent areas in which tick populations are widely established, establishing, and considered non-established. Using a community science approach by partnering with veterinarians and members of the public, we collected over 7000 ticks from the 3 provinces. The three species found most often on companion animals and humans were I. scapularis (76.9%), Ixodes cookei (10.4%) and Dermacentor variabilis (8.9%). The most common hosts were dogs (60.5%), cats (16.8%) and humans (17.6%). As is typical of passive surveillance tick collections, the majority of ticks recovered were adult females; for I. scapularis 90.2%, 5.3%, 3.9% and 0.6% of the total of 5630 ticks recovered for this species were adult females, adult males, nymphs and larvae, respectively. The majority of B. burgdorferi-infected ticks were I. scapularis, as expected. Borrelia infection prevalence in I scapularis was higher in Nova Scotia (20.9%), the province with the most endemic regions, than New Brunswick (14.1%) and Prince Edward Island (9.1%), provinces thought to have established and non-established tick populations, respectively. The province-wide Borrelia infection prevalence generally increased in these latter tow provinces over the course of the study. The host did not have a significant effect on B. burgdorferi infection prevalence; I. scapularis ticks from dogs, cats, humans was, 13.3% (n = 3622), 15.6% (n = 817), 17.9% (n = 730), respectively. No I. scapularis larvae were found infected (n = 33) but B. burgdorferi was detected in 14.8% of both adults (n = 5140) and nymphs (n = 215). The incidence of B. burgdorferi infection also did not differ by engorgement status 15.0% (n = 367), 15.1% (n = 3101) and 14.4% (n = 1958) of non-engorged, engorged and highly engorged ticks, respectively, were infected. In New Brunswick, at the advancing front of tick population establishment, the province-wide infection percentages generally increased over the nine-year study period and all health district regions showed increased tick recoveries and a trend of increased percentages of Borrelia-infected ticks over the course of the study. Within New Brunswick, tick recoveries but not Borrelia infection prevalence were significantly different from endemic and non-endemic regions, suggesting cryptic endemic regions existed prior to their designation as a risk area. Over the 9 years of the study, tick recoveries increased in New Brunswick, the primary study region, and I. scapularis recoveries spread northwards and along the coast, most but not all new sites of recoveries were predicted by climate-based models, indicating that ongoing tick surveillance is necessary to accurately detect all areas of risk. Comparison of tick recoveries and public health risk areas indicates a lag in identification of risk areas. Accurate and timely information on tick distribution and the incidence of Borrelia and other infections are essential for keeping the public informed of risk and to support disease prevention behaviors.

Addressing the global snakebite crisis with geo-spatial analyses - Recent advances and future direction

Venomous snakebite is a neglected tropical disease that annually leads to hundreds of thousands of deaths or long-term physical and mental ailments across the developing world. Insufficient data on spatial variation in snakebite risk, incidence, human vulnerability, and accessibility of medical treatment contribute substantially to ineffective on-ground management. There is an urgent need to collect data, fill knowledge gaps and address on-ground management problems. The use of novel, and transdisciplinary approaches that take advantage of recent advances in spatio-temporal models, 'big data', high performance computing, and fine-scale spatial information can add value to snakebite management by strategically improving our understanding and mitigation capacity of snakebite. We review the background and recent advances on the topic of snakebite related geospatial analyses and suggest avenues for priority research that will have practical on-ground applications for snakebite management and mitigation. These include streamlined, targeted data collection on snake distributions, snakebites, envenomings, venom composition, health infrastructure, and antivenom accessibility along with fine-scale models of spatio-temporal variation in snakebite risk and incidence, intraspecific venom variation, and environmental change modifying human exposure. These measures could improve and 'future-proof' antivenom production methods, antivenom distribution and stockpiling systems, and human-wildlife conflict management practices, while simultaneously feeding into research on venom evolution, snake taxonomy, ecology, biogeography, and conservation.

Genetic Melting Pot in Blacklegged Ticks at the Northern Edge of their Expansion Front

Blacklegged ticks (Ixodes scapularis) are considered to be the main vector of Lyme disease in eastern North America. They may parasitize a wide range of bird and mammal hosts. Northward dispersal of blacklegged ticks has been attributed largely to movement of hosts to areas outside of the current range of the tick, in conjunction with climate change. To better understand the drivers of range expansion in the blacklegged tick, we need investigations of the genetic connectivity and differentiation of tick populations at a fine spatial scale using appropriate markers. In this study, we investigated genetic connectivity and differentiation in blacklegged ticks, in an area of putatively recent advance in Ontario and Quebec, Canada, using microsatellite markers. Our findings suggest patchy differentiation of alleles, no spatial pattern of genetic structure, and genetic subdivision within sites, which are consistent with the very limited evidence available near the leading edge of range expansion of blacklegged ticks into Canada. These findings are consistent with the prevailing hypothesis, drawn from a variety of fields of study, suggesting that migratory birds from a variety of regions may be bringing hitchhiking ticks northward into Canada.

Statewide Passive Surveillance of Ixodes scapularis and Associated Pathogens in Maine

The blacklegged tick, Ixodes scapularis, is the primary vector of multiple human pathogens, including the causative agents of Lyme disease, anaplasmosis, and babesiosis. Both I. scapularis and its associated pathogens have expanded their geographic range throughout the northeastern Unites States and into northern New England. Through this study, we present an updated distribution of I. scapularis in Maine and report the first statewide passive surveillance infection and coinfection prevalence of Borrelia burgdorferi, Anaplasma phagocytophilum, and Babesia microti within the state's I. scapularis population. In 2019, we collected 2016 ticks through a passive surveillance program, in which Maine residents submitted tick samples for identification and/or pathogen testing. We used a single multiplex quantitative PCR assay to detect tickborne pathogens in 1901 tick samples. At the state level, we found that Bo. burgdorferi and A. phagocytophilum infection rates of adults (42.4%, 11.1%) were nearly double that of nymphs (26.9%, 6.7%), whereas B. microti prevalence was similar for both adults (6.5%) and nymphs (5.2%). Spatially, we found an uneven distribution of both tick activity and pathogen prevalence, with both increasing on a north to south gradient. We also noted a potential association between the ratio of adult to nymphal ticks and the incidence of tickborne disease in human populations, with counties that exhibit high rates of human disease also maintaining low adult to nymph ratios.

The Potential Distribution of Pythium insidiosum in the Chincoteague National Wildlife Refuge, Virginia

Pythium insidiosum is a widespread pathogen that causes pythiosis in mammals. Recent increase in cases reported in North America indicates a need to better understand the distribution and persistence of the pathogen in the environment. In this study, we reconstructed the distribution of P. insidiosum in the Chincoteague National Wildlife Refuge, located on Assateague Island, Virginia, and based on 136 environmental water samples collected between June and September of 2019. The Refuge hosts two grazing areas for horses, also known as the Chincoteague Ponies. In the past 3 years, 12 horses have succumbed to infection by P. insidiosum. Using an ecological niche model framework, we estimated and mapped suitable areas for P. insidiosum throughout the Refuge. We found P. insidiosum throughout much of the study area. Our results showed significant monthly variation in the predicted suitability, where the most influential environmental predictors were land-surface water and temperature. We found that June, July, and August were the months with the highest suitability for P. insidiosum across the Refuge, while December through March were less favorable months. Likewise, significant differences in suitability were observed between the two grazing areas. The suitability map provided here could also be used to make management decisions, such as monitoring horses for lesions during high risk months.

A multi-year assessment of blacklegged tick (Ixodes scapularis) population establishment and Lyme disease risk areas in Ottawa, Canada, 2017-2019

Canadians face an emerging threat of Lyme disease due to the northward expansion of the tick vector, Ixodes scapularis. We evaluated the degree of I. scapularis population establishment and Borrelia burgdorferi occurrence in the city of Ottawa, Ontario, Canada from 2017–2019 using active surveillance at 28 sites. We used a field indicator tool developed by Clow et al. to determine the risk of I. scapularis establishment for each tick cohort at each site using the results of drag sampling. Based on results obtained with the field indicator tool, we assigned each site an ecological classification describing the pattern of tick colonization over two successive cohorts (cohort 1 was comprised of ticks collected in fall 2017 and spring 2018, and cohort 2 was collected in fall 2018 and spring 2019). Total annual site-specific I. scapularis density ranged from 0 to 16.3 ticks per person-hour. Sites with the highest density were located within the Greenbelt zone, in the suburban/rural areas in the western portion of the city of Ottawa, and along the Ottawa River; the lowest densities occurred at sites in the suburban/urban core. B. burgdorferi infection rates exhibited a similar spatial distribution pattern. Of the 23 sites for which data for two tick cohorts were available, 11 sites were classified as “high-stable”, 4 were classified as “emerging”, 2 were classified as “low-stable”, and 6 were classified as “non-zero”. B. burgdorferi-infected ticks were found at all high-stable sites, and at one emerging site. These findings suggest that high-stable sites pose a risk of Lyme disease exposure to the community as they have reproducing tick populations with consistent levels of B. burgdorferi infection. Continued surveillance for I. scapularis, B. burgdorferi, and range expansion of other tick species and emerging tick-borne pathogens is important to identify areas posing a high risk for human exposure to tick-borne pathogens in the face of ongoing climate change and urban expansion.

A Generalized Additive Model Correlating Blacklegged Ticks With White-Tailed Deer Density, Temperature, and Humidity in Maine, USA, 1990-2013

Geographical range expansions of blacklegged tick [Ixodes scapularis Say (Acari: Ixodidae)] populations over time in the United States have been attributed to a mosaic of factors including 20th century reforestation followed by suburbanization, burgeoning populations of the white-tailed deer [Odocoileus virginianus Zimmerman (Artiodactyla: Cervidae)], and, at the northern edge of I. scapularis' range, climate change. Maine, a high Lyme disease incidence state, has been experiencing warmer and shorter winter seasons, and relatively more so in its northern tier. Maine served as a case study to investigate the interacting impacts of deer and seasonal climatology on the spatial and temporal distribution of I. scapularis. A passive tick surveillance dataset indexed abundance of I. scapularis nymphs for the state, 1990-2013. With Maine's wildlife management districts as the spatial unit, we used a generalized additive model to assess linear and nonlinear relationships between I. scapularis nymph abundance and predictors. Nymph submission rate increased with increasing deer densities up to ~5 deer/km2 (13 deer/mi2), but beyond this threshold did not vary with deer density. This corroborated the idea of a saturating relationship between I. scapularis and deer density. Nymphs also were associated with warmer minimum winter temperatures, earlier degree-day accumulation, and higher relative humidity. However, nymph abundance only increased with warmer winters and degree-day accumulation where deer density exceeded ~2 deer/km2 (~6/mi2). Anticipated increases in I. scapularis in the northern tier could be partially mitigated through deer herd management.

Predicting the current and future distribution of the western black-legged tick, Ixodes pacificus, across the Western US using citizen science collections

In the twenty-first century, ticks and tick-borne diseases have expanded their ranges and impact across the US. With this spread, it has become vital to monitor vector and disease distributions, as these shifts have public health implications. Typically, tick-borne disease surveillance (e.g., Lyme disease) is passive and relies on case reports, while disease risk is calculated using active surveillance, where researchers collect ticks from the environment. Case reports provide the basis for estimating the number of cases; however, they provide minimal information on vector population or pathogen dynamics. Active surveillance monitors ticks and sylvatic pathogens at local scales, but it is resource-intensive. As a result, data are often sparse and aggregated across time and space to increase statistical power to model or identify range changes. Engaging public participation in surveillance efforts allows spatially and temporally diverse samples to be collected with minimal effort. These citizen-driven tick collections have the potential to provide a powerful tool for tracking vector and pathogen changes. We used MaxEnt species distribution models to predict the current and future distribution of Ixodes pacificus across the Western US through the use of a nationwide citizen science tick collection program. Here, we present niche models produced through citizen science tick collections over two years. Despite obvious limitations with citizen science collections, the models are consistent with previously-predicted species ranges in California that utilized more than thirty years of traditional surveillance data. Additionally, citizen science allows for an expanded understanding of I. pacificus distribution in Oregon and Washington. With the potential for rapid environmental changes instigated by a burgeoning human population and rapid climate change, the development of tools, concepts, and methodologies that provide rapid, current, and accurate assessment of important ecological qualities will be invaluable for monitoring and predicting disease across time and space.

Nested Species Distribution Models of Chlamydiales in Ixodes ricinus (Tick) Hosts in Switzerland

The tick Ixodes ricinus is the vector of various pathogens, including Chlamydiales bacteria, which potentially cause respiratory infections. In this study, we modeled the spatial distribution of I. ricinus and associated Chlamydiales over Switzerland from 2009 to 2019. We used a total of 2,293 ticks and 186 Chlamydiales occurrences provided by a Swiss Army field campaign, a collaborative smartphone application, and a prospective campaign. For each tick location, we retrieved from Swiss federal data sets the environmental factors reflecting the topography, climate, and land cover. We then used the Maxent modeling technique to estimate the suitability of particular areas for I. ricinus and to subsequently build the nested niche of Chlamydiales bacteria. Results indicate that I. ricinus habitat suitability is determined by higher temperature and normalized difference vegetation index (NDVI) values, lower temperature during the driest months, and a higher percentage of artificial and forest areas. The performance of the model was improved when extracting the environmental variables for a 100-m radius buffer around the sampling points and when considering the climatic conditions of the 2 years previous to the sampling date. Chlamydiales bacteria were favored by a lower percentage of artificial surfaces, drier conditions, high precipitation during the coldest months, and short distances to wetlands. From 2009 to 2018, we observed an extension of areas suitable to ticks and Chlamydiales, associated with a shift toward higher altitude. The importance of considering spatiotemporal variations in the environmental conditions for obtaining better prediction was also demonstrated.IMPORTANCEIxodes ricinus is the vector of pathogens including the agent of Lyme disease, the tick-borne encephalitis virus, and the less well-known Chlamydiales bacteria, which are responsible for certain respiratory infections. In this study, we identified the environmental factors influencing the presence of I. ricinus and Chlamydiales in Switzerland and generated maps of their distribution from 2009 to 2018. We found an important expansion of suitable areas for both the tick and the bacteria during the last decade. Results also provided the environmental factors that determine the presence of Chlamydiales within ticks. Distribution maps as generated here are expected to bring valuable information for decision makers in controlling tick-borne diseases in Switzerland and establishing prevention campaigns. The methodological framework presented could be used to predict the distribution and spread of other host-pathogen pairs to identify environmental factors driving their distribution and to develop control or prevention strategies accordingly.

Applying species distribution models in public health research by predicting snakebite risk using venomous snakes' habitat suitability as an indicating factor

Snakebite envenoming is an important public health problem in Iran, despite its risk not being quantified. This study aims to use venomous snakes' habitat suitability as an indicator of snakebite risk, to identify high-priority areas for snakebite management across the country. Thus, an ensemble approach using five distribution modelling methods: Generalized Boosted Models, Generalized Additive Models, Maximum Entropy Modelling, Generalized Linear Models, and Random Forest was applied to produce a spatial snakebite risk model for Iran. To achieve this, four venomous snakes' habitat suitability (Macrovipera lebetinus, Echis carinatus, Pseudocerastes persicus and Naja oxiana) were modelled and then multiplied. These medically important snakes are responsible for the most snakebite incidents in Iran. Multiplying habitat suitability models of the four snakes showed that the northeast of Iran (west of Khorasan-e-Razavi province) has the highest snakebite risk in the country. In addition, villages that were at risk of envenoming from the four snakes were identified. Results revealed that 51,112 villages are at risk of envenoming from M. lebetinus, 30,339 from E. carinatus, 51,657 from P. persicus and 12,124 from N. oxiana. Precipitation seasonality was identified as the most important variable influencing distribution of the P. persicus, E. carinatus and M. lebetinus in Iran. Precipitation of the driest quarter was the most important predictor of suitable habitats of the N. oxiana. Since climatic variables play an important role in shaping the distribution of the four venomous snakes in Iran, thus their distribution may alter with changing climate. This paper demonstrates application of species distribution modelling in public health research and identified potential snakebite risk areas in Iran by using venomous snakes' habitat suitability models as an indicating factor. Results of this study can be used in snakebite and human-snake conflict management in Iran. We recommend increasing public awareness of snakebite envenoming and education of local people in areas which identified with the highest snakebite risk.

Habitat suitability model of endangered Latidens salimalii and the probable consequences of global warming

Salim Ali's fruit bat, Latidens salimalii, is a monotypic endangered fruit bat endemic to Western Ghats (WG) with an ambiguous distribution. The distribution range, habitat suitability, and biology of this species are still uncertain. Endemic species inhabiting the high elevation of  WG like L. salimalii are threatened due to climatic change and seeks urgent management interventions. Hence, we developed a habitat suitability model for L. salimalii using MaxEnt in the current climate condition and projected their distribution for three Representation Concentration Pathway (RCP 4.5, 6.0, and 8.5) climate scenarios of the 2070 time frame. The results show that 9531 km²of habitat in WG is suitable for L. salimalii at present, while all the future scenarios estimates propose complete loss of highly suitable habitat. The significant factors influencing the distribution of L. salimalii are the precipitation of the driest month, tree density, rain in the coldest quarter, canopy height, and altitude. The study pioneers in predicting the suitable habitat and emphasis the need to develop strategies for the long-term conservation of endangered L. salimalii in WG under global warming scenarios.

Species distribution models for the eastern blacklegged tick, Ixodes scapularis, and the Lyme disease pathogen, Borrelia burgdorferi, in Ontario, Canada

The blacklegged tick, Ixodes scapularis, is established in several regions of Ontario, Canada, and continues to spread into new geographic areas across the province at a rapid rate. This poses a significant public health risk since I. scapularis transmits the Lyme disease-causing bacterium, Borrelia burgdorferi, and other pathogens of potential public health concern. The objective of this study was to develop species distribution models for I. scapularis and B. burgdorferi to predict and compare the potential distributions of the tick vector and the Lyme disease pathogen as well as the ecological factors most important for species establishment. Ticks were collected via tick dragging at 120 sites across southern, central, and eastern Ontario between 2015 and 2018 and tested for tick-borne pathogens. A maximum entropy (Maxent) approach was used to model the potential distributions of I. scapularis and B. burgdorferi. Two independent datasets derived from tick dragging at 25 new sites in 2019 and ticks submitted by the public to local health units between 2015 and 2017 were used to validate the predictive accuracy of the models. The model for I. scapularis showed high suitability for blacklegged ticks in eastern Ontario and some regions along the shorelines of the Great Lakes, and moderate suitability near Algonquin Provincial Park and the Georgian Bay with good predictive accuracy (tick dragging 2019: AUC = 0.898; ticks from public: AUC = 0.727). The model for B. burgdorferi showed a similar predicted distribution but was more constrained to eastern Ontario, particularly between Ottawa and Kingston, and along Lake Ontario, with similarly good predictive accuracy (tick dragging 2019: AUC = 0.958; ticks from public: AUC = 0.863. The ecological variables most important for predicting the distributions of I. scapularis and B. burgdorferi included elevation, distance to deciduous and coniferous forest, proportions of agricultural land, water, and infrastructure, mean summer/spring temperature, and cumulative annual degree days above 0°C. Our study presents a novel application of species distribution modelling for I. scapularis and B. burgdorferi in Ontario, Canada, and provides an up to date projection of their potential distributions for public health knowledge users.

MaxEnt Modeling of Dermacentor marginatus (Acari: Ixodidae) Distribution in Xinjiang, China

Dermacentor marginatus Sulkzer is a common tick species found in the Xinjiang Uygur Autonomous Region (XUAR) of China, and is a vector for a variety of pathogens. To determine the potential distribution of this tick species in Xinjiang, a metadata containing 84 D. marginatus presence records combined with four localities from field collection were used for MaxEnt modeling to predict potential distribution of this tick species. Identification of tick samples showed 756 of 988 (76%) were D. marginatus. MaxEnt modeling results indicated that the potential distribution of this tick species was mainly confined to northern XUAR. Highly suitable areas included west side of Altay mountain, west rim of Junggar basin, and Yili River valley in the study area. The model showed an AUC value of 0.838 ± 0.063 (SD), based on 10-fold cross-validation. Although tick presence records used for modeling were limited, this is the first regional tick distribution model for D. marginatus in Xinjiang. The model will be helpful in assessing the risk of tick-borne diseases to human and animals in the region.

Modelling habitat suitability for occurrence of human tick-borne encephalitis (TBE) cases in Finland

The numbers of reported human tick-borne encephalitis (TBE) cases in Europe have increased in several endemic regions (including Finland) in recent decades, indicative of an increasing threat to public health. As such, it is important to identify the regions at risk and the most influential factors associated with TBE distributions, particularly in understudied regions. This study aimed to identify the risk areas of TBE transmission in two different datasets based on human TBE disease cases from 2007 to 2011 (n = 86) and 2012-2017 (n  = 244). We also examined which factors best explain the presence of human TBE cases. We used ensemble modelling to determine the relationship of TBE occurrence with environmental, ecological, and anthropogenic factors in Finland. Geospatial data including these variables were acquired from several open data sources and satellite and aerial imagery and, were processed in GIS software. Biomod2, an ensemble platform designed for species distribution modelling, was used to generate ensemble models in R. The proportion of built-up areas, field, forest, and snow-covered land in November, people working in the primary sector, human population density, mean precipitation in April and July, and densities of European hares, white-tailed deer, and raccoon dogs best estimated distribution of human TBE disease cases in the two datasets. Random forest and generalized boosted regression models performed with a very good to excellent predictive power (ROC = 0.89-0.96) in both time periods. Based on the predictive maps, high-risk areas for TBE transmission were located in the coastal regions in Southern and Western Finland (including the Åland Islands), several municipalities in Central and Eastern Finland, and coastal municipalities in Southern Lapland. To explore potential changes in TBE distributions in future climate, we used bioclimatic factors with current and future climate forecast data to reveal possible future hotspot areas. Based on the future forecasts, a slightly wider geographical extent of TBE risk was introduced in the Åland Islands and Southern, Western and Northern Finland, even though the risk itself was not increased. Our results are the first steps towards TBE-risk area mapping in current and future climate in Finland.

Landscape determinants of density of blacklegged ticks, vectors of Lyme disease, at the northern edge of their distribution in Canada

In eastern North America, including Canada, Lyme disease is caused by Borrelia burgdorferi sensu stricto and transmitted to humans by the blacklegged tick, Ixodes scapularis. The last decade has seen a growing incidence of Lyme disease in Canada, following the northward range expansion of I. scapularis tick populations from endemic areas in eastern United States. This may be attributable to movement of the many hosts that they parasitize, including songbirds, deer and small mammals. In this study, we wanted to test the effect of spatial, temporal and ecological variables, on blacklegged tick density and infection rates, near the northern limit of their distribution in Ontario and Quebec, Canada. We found an effect of both proportion of forested areas and distance to roads, on density of I. scapularis ticks and prevalence of infection by B. burgdorferi. We also found an effect of both sampling year and ordinal sampling data on prevalence of infection by B. burgdorferi. In six adjacent sites showing evidence of reproducing I. scapularis populations, we found that forest composition and structure influenced density of I. scapularis ticks. Our results suggest that blacklegged tick density and infection rate in Canada may be influenced by a variety of factors.

High-resolution remote sensing-based spatial modeling for the prediction of potential risk areas of schistosomiasis in the Dongting Lake area, China

The geographical distribution of snail (i.e., the intermediate host of schistosomiasis) is consistent with that of endemic areas. The suitable snail habitus requires necessary environmental conditions for snail population. The high-resolution remote sensing provides an important tool for the spatio-temporal analysis of disease monitoring and prediction. This study conducted a typical schistosomiasis epidemic area in the marshland and lake regions along the Yangtze River, Yueyang City, Hunan Province of China. And three types of environmental factors, i.e., NDVI, soil moisture, and shortest distance to water body, associated with the geographical distribution of snail population, were extracted from the high-resolution remoting sensing data. The predicted distribution of snail habitus from the high-resolution environmental factors were compared with the data of annual program of snail survey. The results have shown that the application of high-resolution remote sensing can improve the accuracy of the modeled and predicted the potential risk areas of schistosomiasis, and may become an important tool for the ongoing national schistosomiasis control program.

High-resolution remote sensing-based spatial modeling for the prediction of potential risk areas of schistosomiasis in the Dongting Lake area, China

The geographical distribution of snail (i.e., the intermediate host of schistosomiasis) is consistent with that of endemic areas. The suitable snail habitus requires necessary environmental conditions for snail population. The high-resolution remote sensing provides an important tool for the spatio-temporal analysis of disease monitoring and prediction. This study conducted a typical schistosomiasis epidemic area in the marshland and lake regions along the Yangtze River, Yueyang City, Hunan Province of China. And three types of environmental factors, i.e., NDVI, soil moisture, and shortest distance to water body, associated with the geographical distribution of snail population, were extracted from the high-resolution remoting sensing data. The predicted distribution of snail habitus from the high-resolution environmental factors were compared with the data of annual program of snail survey. The results have shown that the application of high-resolution remote sensing can improve the accuracy of the modeled and predicted the potential risk areas of schistosomiasis, and may become an important tool for the ongoing national schistosomiasis control program.

Ecological and Socioeconomic Factors Associated with Bartonella henselae Exposure in Dogs Tested for Vector-Borne Diseases in North Carolina

Bartonella henselae is a zoonotic vector-borne pathogen affecting both humans and dogs. Little is known about the epidemiology of B. henselae in dogs, including risk factors associated with exposure. The objectives of this study were to map the current distribution of B. henselae in dogs in North Carolina (NC) and to identify ecological and socioeconomic factors influencing B. henselae seroreactivity. Results from 4446 B. henselae serology samples from dogs in NC submitted by veterinarians for clinical diagnostic testing to the North Carolina State University College of Veterinary Medicine Vector Borne Disease Diagnostic Laboratory between January 1, 2004 and December 31, 2015 were retrospectively reviewed. These results were used to generate a map of B. henselae seroreactivity. To account for sparsely sampled areas, statistical smoothing using head banging and areal interpolation kriging was performed. Using previously described risk factors for exposure to canine tick-borne diseases, eight multivariable logistic regression models based on biologically plausible hypotheses were tested, and a final model was selected using an Akaike's Information Criterion weighted-average approach. Seroreactivity among dogs tested for vector-borne disease was variable across the state: higher along the southern/eastern coastal plains and eastern Piedmont, and lower in the western mountains. Of 25 explanatory factors considered, the model combining demographic, socioeconomic, climatic, and land use variables fits best. Based on this model, female intact sex and increasing percentage of the county with low-intensity development and evergreen forest were associated with higher seroreactivity. Conversely, moderate development, increasing median household income, and higher temperature range and relative humidity were associated with lower seroreactivity. This model could be improved, however, by including local and host-scale factors that may play a significant role in dogs' exposure.

What is the real number of Lyme disease cases in Canada?

BACKGROUND

Lyme disease is emerging in Canada due to expansion of the range of the tick vector Ixodes scapularis from the United States. National surveillance for human Lyme disease cases began in Canada in 2009. Reported numbers of cases increased from 144 cases in 2009 to 2025 in 2017. It has been claimed that few (< 10%) Lyme disease cases are reported associated with i) supposed under-diagnosis resulting from perceived inadequacies of serological testing for Lyme disease, ii) expectation that incidence in Canadian provinces and neighbouring US states should be similar, and iii) analysis of serological responses of dogs to the agent of Lyme disease, Borrelia burgdorferi. We argue that performance of serological testing for Lyme disease is well studied, and variations in test performance at different disease stages are accounted for in clinical diagnosis of Lyme disease, and in surveillance case definitions. Extensive surveillance for tick vectors has taken place in Canada providing a clear picture of the emergence of risk in the Canadian environment. This surveillance shows that the geographic scope of I. scapularis populations and Lyme disease risk is limited but increasing in Canada. The reported incidence of Lyme disease in Canada is consistent with this pattern of environmental risk, and the differences in Lyme disease incidence between US states and neighbouring Canadian provinces are consistent with geographic differences in environmental risk. Data on serological responses in dogs from Canada and the US are consistent with known differences in environmental risk, and in numbers of reported Lyme disease cases, between the US and Canada.

CONCLUSION

The high level of consistency in data from human case and tick surveillance, and data on serological responses in dogs, suggests that a high degree of under-reporting in Canada is unlikely. We speculate that approximately one third of cases are reported in regions of emergence of Lyme disease, although prospective studies are needed to fully quantify under-reporting. In the meantime, surveillance continues to identify and track the ongoing emergence of Lyme disease, and the risk to the public, in Canada.

An Ecological Framework for Modeling the Geography of Disease Transmission

Ecological niche modeling (ENM) is widely employed in ecology to predict species’ potential geographic distributions in relation to their environmental constraints and is rapidly becoming the gold-standard method for disease risk mapping. However, given the biological complexity of disease systems, the traditional ENM framework requires reevaluation. We provide an overview of the application of ENM to disease systems and propose a theoretical framework based on the biological properties of both hosts and parasites to produce reliable outputs resembling disease system distributions. Additionally, we discuss the differences between biological considerations when implementing ENM for distributional ecology and epidemiology. This new framework will help the field of disease ecology and applications of biogeography in the epidemiology of infectious diseases.

Infectious diseases greatly impact human health, biodiversity, and global economies, highlighting the need to understand and predict their distributions.

Ecological niche modeling (ENM) was not originally designed to explicitly reconstruct complex biological phenomena such as diseases or parasitism, requiring a reevaluation of the traditional framework.

We provide an integrative ENM framework for disease systems that considers suitable host availability, parasite ecologies, and different scales of modeling.

Disease transmission is driven by factors related to parasite availability and host exposure and susceptibility, which can be incorporated in ENM frameworks.

Citizen science informs human-tick exposure in the Northeastern United States

BACKGROUND

Tick-borne disease is the result of spillover of pathogens into the human population. Traditionally, literature has focused on characterization of tick-borne disease pathogens and ticks in their sylvatic cycles. A limited amount of research has focused on human-tick exposure in this system, especially in the Northeastern United States. Human-tick interactions are crucial to consider when assessing the risk of tick-borne disease since a tick bite is required for spillover to occur.

METHODS

Citizen scientists collected ticks from the Northeastern US through a free nationwide program. Submitted ticks were identified to species, stage, and sex. Blacklegged ticks, Ixodes scapularis, were tested for the presence of Borrelia burgdorferi sensu lato (s.l.) and hard-tick relapsing fever Borrelia. Seasonality of exposure and the citizen science activity during tick exposure was recorded by the citizen scientist. A negative binomial model was fit to predict county level CDC Lyme disease cases in 2016 using citizen science Ixodes scapularis submissions, state, and county population as predictor variables.

RESULTS

A total of 3740 submissions, comprising 4261 ticks, were submitted from the Northeastern US and were reported to be parasitizing humans. Of the three species submitted, blacklegged ticks were the most prevalent followed by American dog ticks and lone star ticks. Submissions peaked in May with the majority of exposure occurring during every-day activities. The most common pathogen in blacklegged ticks was B. burgdorferi s.l. followed by hard-tick relapsing fever Borrelia. Negative binomial model performance was best in New England states followed by Middle Atlantic states.

CONCLUSIONS

Citizen science provides a low-cost and effective methodology for describing the seasonality and characteristics of human-tick exposure. In the Northeastern US, everyday activities were identified as a major mechanism for tick exposure, supporting the role of peri-domestic exposure in tick-borne disease. Citizen science provides a method for broad pathogen and tick surveillance, which is highly related to human disease, allowing for inferences to be made about the epidemiology of tick-borne disease.

Predicting spatiotemporal patterns of Lyme disease incidence from passively collected surveillance data for Borrelia burgdorferi sensu lato-infected Ixodes scapularis ticks

Lyme disease is the most prevalent vector-borne disease in the United States. Ixodes scapularis, commonly referred to as the blacklegged tick, is the primary vector of Lyme disease spirochetes, Borrelia burgdorferi sensu lato (s.l.), in the eastern United States. Connecticut has pervasive populations of I. scapularis and remains a hotspot for Lyme disease. A primary aim of this study was to determine if passively collected data on human-biting I. scapularis ticks in Connecticut could serve as a useful proxy for Lyme disease incidence based on the cases reported by the Connecticut Department of Public Health (CDPH). Data for human-biting I. scapularis ticks submitted to the Tick Testing Laboratory at the Connecticut Agricultural Experiment Station (CAES-TTL), and tested for infection with B. burgdorferi s.l., were used to estimate the rate of submitted nymphs, nymphal infection prevalence, and the rate of submitted infected nymphs. We assessed spatiotemporal patterns in tick-based measures and Lyme disease incidence with generalized linear and spatial models. In conjunction with land cover and household income data, we used generalized linear mixed effects models to examine the association between tick-based risk estimates and Lyme disease incidence. Between 2007 and 2017, the CAES-TTL received 26,116 I. scapularis tick submissions and the CDPH reported 23,423 Lyme disease cases. The rate of submitted nymphs, nymphal infection prevalence, the rate of submitted infected nymphs, and Lyme disease incidence all decreased over time during this eleven-year period. The rate of submitted nymphs, the rate of submitted infected nymphs, and Lyme disease incidence were spatially correlated, but nymphal infection prevalence was not. Using a mixed modeling approach to predict Lyme disease incidence and account for spatiotemporal structuring of the data, we found the best fitting tested model included a strong, positive association with the rate of submitted infected nymphs and a negative association with the percent of developed land for each county. We show that within counties, submissions of B. burgdorferi s.l. infected nymphs were strongly and positively associated with inter-annual variation in reported Lyme disease cases. Tick-based passive surveillance programs may be useful in providing independent measures of entomological risk, particularly in settings where Lyme disease case reporting practices change substantially over time.

Ixodes scapularis tick distribution and infection rates in Ottawa, Ontario, 2017

Background

The incidence of Lyme disease has increased in many regions of Canada in recent years, including in Ottawa, Ontario. To date there has been limited active tick surveillance in the region.

Objectives

To estimate both the distribution and density of Ixodes scapularis ticks in the city of Ottawa, and the infection rates of ticks with Borrelia burgdorferi (that causes Lyme disease) and other tick-borne pathogens.

Methods

Between June and October 2017, tick surveillance was conducted by drag sampling at 23 sites in Ottawa municipal parks, recreational trails and forests. Blacklegged ticks were tested for B. burgdorferi, Borrelia miyamotoi and Anaplasma phagocytophilum using quantitative polymerase chain reaction protocols.

Results

I. scapularis ticks were found in 16 of the 23 sites (70%). Recreational trails, conservation areas/forests and the provincial park within the city of Ottawa had significantly higher tick densities than municipal parks (p<0.01). Of the 194 adult and 26 nymphal I. scapularis tested, prevalence of infection was 29.5% for B. burgdorferi, 0.45% for B. miyamotoi and 0.91% for A. phagocytophilum.

Conclusion

Almost 30% of I. scapularis ticks tested in suburban and rural areas of the city of Ottawa were infected with B. burgdorferi, known to cause Lyme disease. Other types of infection, known to cause anaplasmosis and tick-borne relapsing fever, were also detected, although were very rare. Conducting active tick surveillance at the local level may help to inform risk assessment and public health actions.