Abstract
Used as a communicative tool for risk management, risk maps provide a service to the public, conveying information that can raise risk awareness and encourage mitigation. Several studies have utilized risk maps to determine risks associated with the distribution of Borrelia burgdorferi, the causal agent of Lyme disease in North America and Europe, as this zoonotic disease can lead to severe symptoms. This literature review focused on the use of risk maps to model distributions of B. burgdorferi and its vector, the blacklegged tick (Ixodes scapularis), in North America to compare variables used to predict these spatial models. Data were compiled from the existing literature to determine which ecological, environmental, and anthropic (i.e., human focused) variables past research has considered influential to the risk level for Lyme disease. The frequency of these variables was examined and analyzed via a non-metric multidimensional scaling analysis to compare different map elements that may categorize the risk models performed. Environmental variables were found to be the most frequently used in risk spatial models, particularly temperature. It was found that there was a significantly dissimilar distribution of variables used within map elements across studies: Map Type, Map Distributions, and Map Scale. Within these map elements, few anthropic variables were considered, particularly in studies that modeled future risk, despite the objective of these models directly or indirectly focusing on public health intervention. Without including human-related factors considering these variables within risk map models, it is difficult to determine how reliable these risk maps truly are. Future researchers may be persuaded to improve disease risk models by taking this into consideration.
Introduction
Tick borne diseases are caused by pathogens transmitted by infected ticks to an uninfected host. As the climate warms, it becomes possible for ticks to have increased abundance, survival, and feeding activity, and to expand their geographic distribution northwards (, ). Human activities and their impacts on natural habitats are further altering the distribution of these disease vectors, and human-wildlife contacts are increasing through means of socio-demographics (globalization, urbanization, etc...) and public health systems [vector control and other health interventions ()]. Due to the constant changes in the environment, it is challenging to assess the risks associated with these infections without predictive modeling. Many current and future scenario predictive risk maps have been developed to monitor zoonotic infections for public health interventions (–). Epidemiologic risk maps are effective visualization tools used to identify geographical areas of high risk for disease transmission and potential future high-risk regions. Used as a communicative tool for risk management, these maps provide a service to the public, conveying information that can raise risk awareness and encourage mitigation strategies (–). Since climate and land use changes are constantly altering the dynamics between vector and host, continuous monitoring of the emergence and expansion of the disease vector is required. This is true of the spread of Lyme disease and other tick-borne diseases (, ).
In eastern North America, Lyme disease is typically caused by an infection of the spirochaete Borrelia burgdorferi via blacklegged ticks [Ixodes scapularis ()]. In humans, the infection can result in a multisystem illness that substantially affects the individual's quality of life, if left untreated (, ). As a common and widespread disease (–) that can lead to severe conditions (), effective communication tools for risk management of this infection are necessary.
Risk maps for Lyme disease in North America tend to focus on the impact climate change has on the distribution of blacklegged ticks (–), since these ectoparasites are vectors for B. burgdorferi () and their geographic range is increasing (, ). The geographic range of blacklegged ticks is heavily dependent on environmental variables such as temperature and precipitation (–). The geographical ranges of their hosts also play a major role (). Blacklegged ticks require a single host for each life stage (, ), and will migrate with these hosts, such as small mammals, birds, and ruminants (, ).
Different variables have been considered to affect Lyme disease distribution. These include ecological variables related to blacklegged ticks [tick density, dispersion (–)] and their small mammal reservoir and migratory avian hosts (, ), as well as environmental variables such as temperature, humidity, and forest fragmentation (, , ). Recent risk map publications have adopted “One Health” approaches, which incorporate sociological, ecological, and biological knowledge into their research (, ). This approach aims to examine and integrate human-related or anthropic variables that may influence human health or risk (, ) including those variables beyond human demographics.
For instance, one may expect that human exposure and the risk of becoming infected by Lyme disease is also dependent on individual human behavior (e.g., knowledge, activity). As such, outdoor workers have been found to be more at risk for zoonotic diseases than those who are outdoors recreationally, due to their degree of exposure to the environment (), and those with immune deficiencies may be more at risk for severe symptoms (). Studies have included surveys of a population to gauge their knowledge on their risks to Lyme disease or tick infections (, , , ). Socio-economic status and ethnicity have also been found to play a role in Lyme disease risk (, , ). Knowing that not all individuals are equally at risk for being infected with diseases, including Lyme disease, it should be expected that studies in which Lyme disease risk maps are developed would include variables associated with human characteristics and behaviors. For these studies to be relevant to public health, variables associated with humans (i.e., social, economic, risk perception) should be taken into consideration, as these factors directly affect the risk posed to the public.
Here, we reviewed the literature to identify the variables past research has considered influential to the distribution of Lyme disease via blacklegged ticks in North America. Variables that researchers routinely included in risk models were examined, and those human variables that were often disregarded but may be informative were highlighted. By calling attention to the lack of human variables found in previous risk maps, future researchers may be persuaded to enhance models by including anthropic factors to improve disease risk prediction.
Methods
Collection of data
We focused on past studies that are comparable due to similarities in geography [same continent, overlapping tick populations (, )], disease vector, and spirochaete strain (B. burgdorferi). For this reason, we focused on one tick vector (black-legged tick) which is endemic to eastern North America (). On March 31, 2023, a descriptive literature review was conducted following methods by Paré and Kitsiou () using Google Scholar, PubMed, and CrossRef with the inclusion criteria terms: “Ixodes scapularis,” “blacklegged ticks,” “risk map” (exact phrase), “Lyme disease,” “risk assessment,” and “B. burgdorferi.” This literature review concentrated on eastern North American (across Canada, the United States, and Mexico) risk assessments of Lyme disease transmitted by blacklegged ticks only. Geographical scale varied across studies, with some focusing on areas at the municipality, provincial/state, regional, or national scale. However, it should be noted that there are several risk assessments for Lyme disease concentrated in Europe and western North America where other tick species and spirochaete vectors can transmit Lyme disease (, , ). Only studies in which analyses included at least one risk map in their results was considered for this review, as we were specifically interested in comparing studies that produced risk map models to evaluate risk. We performed a systematic review following PRISMA (Preferred Reporting Items for Systematic Reviews and Meta-Analyses) guidelines (). The initial search consisted of 145 studies published between 2000 and 2022 with the above criteria. Twenty studies were excluded, as they did not include a risk map in their results. Twelve more studies were excluded since their results were based outside of North America, another two studies were excluded as they focused on different vectors or bacteria and diseases and one article was removed from the literature review as it was retired and no longer considered relevant. Finally, 49 government reports, reviews, and theses were removed, leaving a total of 61 peer-reviewed articles meeting our criteria to be included in analyses (Figure 1; Supplementary Table 1).
Figure 1
Once the literature review was completed, data on the models within the articles were collected. This included collecting information on the number of ecological, environmental, and anthropic variables used in the studies to develop risk map models. Frequency counts for the number of each ecological, environmental, and anthropic variable in a risk map model were recorded. Ecological variables included variables that influence or dictate the relationship between an organism and its environment (i.e., tick occurrence, B. burgdorferi prevalence). Environmental variables included natural resource factors that define an ecosystem or habitat (i.e., temperature, humidity, land cover). Anthropic variables were defined as variables related to human beings (i.e., population density, sex, age). As Lyme disease is primarily transmitted via tick vectors (), these observations do not provide any information on how the disease circulates within a human population. Further, certain predictor variables were simplified to allow comparisons more easily. For example, forest cover and vegetation index were categorized together, as were vapor pressure and humidity, elevation and altitude, human population size and density, and tick abundance and density (Supplementary Table 2).
Additional elements that characterized these risk models were collected and recorded, including the year of publication (Supplementary Table 3), the study's focal location (Country), Map Type (predictive vs. surveillance), the period of the study (year), the focal Scale of study (local, regional, national), the distribution of the study considered– vector (tick) vs. host (human or otherwise) vs. vector and host (both considered)—the Tick Surveillance methodology used for the model (passive vs. active vs. no tick surveillance), and the Tick Life Stage the tick data was based on (immature ticks vs. adult vs. all stages vs. no tick data; Table 1; Supplementary Table 4). Map Type (predictive vs. surveillance) was also included, whereby “predictive” maps referred to future predictive map models, as they predict future scenarios, while current predictive maps will be referred to as “surveillance” models, as they pertain to current risks.
Table 1
| Map element | Anthropic | Environmental | Ecological | |
|---|---|---|---|---|
| Map Type | Surveillance () | 17 | 40 | 52 |
| Predictive () | 3 | 62 | 30 | |
| Map distribution | Host () | 14 | 24 | 19 |
| Vector () | 5 | 62 | 38 | |
| Host and vector () | 2 | 16 | 26 | |
| Map scale | Local () | 14 | 58 | 54 |
| Regional () | 6 | 30 | 18 | |
| National () | 0 | 14 | 10 | |
| Tick life stage | Adult () | 1 | 3 | 4 |
| Immature () | 1 | 30 | 16 | |
| All stages () | 10 | 51 | 49 | |
| No tick data () | 8 | 18 | 13 | |
| Tick surveillance | Passive () | 4 | 23 | 23 |
| Active () | 5 | 27 | 26 | |
| Active and passive () | 0 | 15 | 16 | |
| No surveillance data () | 14 | 22 | 10 | |
| Country | Canada () | 12 | 41 | 34 |
| USA () | 5 | 55 | 45 | |
| USA and Canada () | 2 | 4 | 1 | |
| USA and Mexico () | 1 | 0 | 1 | |
| USA, Canada, and Mexico () | 0 | 2 | 1 |
The cumulative number of anthropic, environmental, and ecological variables used across map elements in the reviewed studies.
The total number of variables for each variable group is the sum of those variables used across studies (n = 61). Numbers in brackets indicate the number of studies for that specific map element.
USA, United States of America.
Most risk map models are predictive as they use modeling to create these maps. However, some models are used to predict risks associated with potential future geographical ranges and distributions of Lyme disease (, , , , ), while others are used to predict the current risk or prevalence (–). In addition, Map Scale was considered and categorized by whether the study focused their spatial scale by country (national), province or state (regional), or a smaller unit (e.g., census division; local).
Synthesis of data
A non-metric multidimensional scaling (NMDS) analysis was performed in R [(); version 1.4.1717] using the “vegan” package (version 2.6-2) to determine any significant differences in the frequency of ecological, environmental, and anthropic variables used across studies (n = 61), depending on their map elements: Map Type, Map Scale, Map Distribution, Year, Tick Life Stage, Map Surveillance method, and Country. Here, the NMDS, which is commonly used as an ordination for community ecology (), was performed where “sites” were the individual studies, and “environmental data” were the map elements that influenced the abundance of “species” (ecological, environmental, and anthropic variables).
Linear models were then used to determine if groups of variables (ecological, environmental, and anthropic) differed in frequency across the different map elements identified by the NMDS. A binomial linear regression was conducted with Map Type as the response variable (surveillance vs. predictive) and the number of ecological, environmental, and anthropic variables used in the studies as predictors. An ordinal linear regression with the “logit” function and “equidistant” threshold was conducted () with Map Distribution (as a factor) as the response variable (host vs. vector. vs. host and vector), using the package “ordinal” [version 2019.10 ()]. Post-hoc tests were conducted, and box plots were used to visualize variation within those map elements that were identified as significant by the NMDS.
Results
Comparatively, certain map elements were more frequently used across the studies assessed in this literature review. For Map Type, there were more predictive maps than surveillance. For Map Distribution, there were more vector (blacklegged tick) distributions considered than human or both human and tick distributions. Maps at the local scale were most common across studies. All stages for Tick Life Stage were most frequently considered vs. specific life stages. Both active and passive data for Tick Surveillance was most frequently included in risk map models and most studies in this assessment were based in the United States (Table 1). The most common variable included in Lyme disease or blacklegged tick risk maps was temperature (n = 28 studies). Of the 10 most frequently used variables in these map models, six were environmental variables, four were ecological, and none were anthropic (Figure 2). The results of the NMDS analysis suggested that there were differences in the frequency of ecological, environmental, and anthropic variables used across Map Type and Map Distribution (Table 2; Figure 3; Supplementary Table 5).
Figure 2
Table 2
| Variables | r2 | Pr (>r) |
|---|---|---|
| Year of publication | 0.0421 | 0.289 |
| Map Type | 0.2268 | 0.001 |
| Distribution | 0.1221 | 0.012 |
| Scale | 0.0350 | 0.372 |
| Tick life stage | 0.0415 | 0.811 |
| Tick surveillance | 0.0352 | 0.653 |
| Country | 0.0468 | 0.730 |
Non-metric multidimensional scaling (NMDS) goodness of fit results of map elements classifying studies based on the frequency of ecological, environmental, and anthropic variables incorporated in each risk map model (n = 61).
Variables include Year of publication, Map Type, Map Distribution, Tick Life Stage, Country. Number of permutations = 999. Statistically significant results are in bold.
Figure 3
When comparing Map Types (surveillance vs. predictive), the difference in use of environmental variables was statistically significant (p < 0.0173), with predictive maps using more of these types of variables. For both types of risk maps, anthropic variables were rarely used (Table 3; Figure 4A). Comparing across Map Distributions, studies that included host distributions in risk models considered more anthropic variables than studies that included vector distributions (p < 0.0152). Meanwhile studies that included both host and vector distributions tended to include more ecological variables (p < 0.0291; Table 4; Figure 4B). There was no significant difference across Year of publication, Surveillance Type used, Map Scale, Country of study origin, or Tick Life Stage focused on in studies (Table 2).
Table 3
| Variable type | Estimate | Standard error | z value | Pr (>|z|) |
|---|---|---|---|---|
| (Intercept) | 0.5878 | 0.8051 | 0.730 | 0.4653 |
| Anthropic | 0.7311 | 0.6579 | 1.111 | 0.2665 |
| Environmental | −0.6633 | 0.2129 | −3.115 | 0.0018 |
| Ecological | 0.6571 | 0.5070 | 1.296 | 0.1950 |
General linear model results for Map Type ~ anthropic + environmental + ecological, where family = “binomial” and the independent variables are counts (the number of variables in each category).
Null deviance = 81.772 on 60 degrees of freedom, residual deviance = 64.676 on 57 degrees of freedom. AIC: 72.676, n = 61. Statistically significant results are in bold.
Figure 4
Table 4
| Variable type | Estimate | Standard error | z value | Pr (>|z|) |
|---|---|---|---|---|
| Anthropic | −1.7432 | 0.6016 | −2.897 | 0.00376 |
| Environmental | −0.3410 | 0.1790 | −1.905 | 0.05681 |
| Ecological | 1.0577 | 0.4504 | 2.348 | 0.01887 |
| Threshold-coefficients | ||||
| Threshold | −1.1201 | 0.7432 | −1.507 | – |
| Spacing | 2.9951 | 0.4817 | 5.657 | – |
Cumulative link model results for map distribution ~ anthropic + environmental + ecological, where link = “logit,” threshold = “equidistant” and the independent variables are counts (the number of variables in each category).
logLik = −51.82, AIC = 113.64, n = 61. Statistically significant results are in bold.
Discussion
Although there is a plethora of literature dedicated to identifying factors that can influence an individual's risk for Lyme disease via tick vectors directly and indirectly (, , ), risk maps that demonstrate the spatial breadth of these risks are less common. The results of this review have shown that within this subset of risk maps, there is no standardized risk score, or variable being used across studies. Temperature was the most common variable used in risk map models, however, it was included in less than half of the risk maps considered. There is extensive research on the relationship between blacklegged ticks and temperature, as it affects a tick's development, survival, and host-seeking behavior (, , ). These are significant factors that influence tick abundance and distribution, and therefore influence the distribution and incidence of Lyme disease (, ). In general, the most common variables included in Lyme disease (or blacklegged tick) risk map models were environmental and ecological, while anthropic were lacking. It should be noted that certain studies, such as Slatculescu et al. () considered many anthropic variables, including population density, walkability scores in an urban setting, median income, and drew conclusions about the contribution of an individual's variability on their Lyme disease risk, however, they did not express these results spatially in a risk map model. Similarly, several other ecological variables often studied and considered influential to Lyme disease and/or blacklegged tick distributions were rarely included in these risk map models. For instance, blacklegged tick distributions are affected by reservoir and migratory host distributions (, ) while Lyme disease distributions can be influenced by the genetic diversity of B. burgdorferi strains (, ). Without considering these variables within risk map models, it is difficult to determine how reliable these risk maps truly are.
Depending on the Map Type and Map Distribution used, studies significantly differed in the number of ecological, environmental, and anthropic variables used to produce risk maps. Studies that utilized different Map Types also significantly varied in their usage of environmental variables, where studies that produced predictive maps used this group of variables more often. Interestingly, predictive map studies also used anthropic variables less often. This suggests that when research focuses on future scenarios, they are reliant on how the environment may change, but do not consider human behavior. Although it is difficult to assess future trends in human behavior, it is still possible to include anthropic variables such as socio-economic status and demographic information () to gain a better understanding of the patterns in human attributes that could influence their risk for Lyme disease.
Across Map Distribution types, studies that included host distributions in their risk maps had more anthropic variables than other distributions (vector only or host and vector distributions). Risk maps that included vector only or both host and vector distributions tended to include more ecological variables. As most host distributions in these risk map models were human distributions, it is logical for human attributed variables to be included in these models. Meanwhile, research that incorporated both host and vector distributions, or only vector distributions may be more likely to include variables from the environment, as they are less directly focused on Lyme disease risk for humans specifically. Rather, they use the result of their risk map model that included disease vector distributions to indirectly make conclusions for public health risks (, , ).
This review was limited to research that included risk maps as results within a study, excluding studies that included maps within introductions and methods for context. In some cases, a single variable (e.g., tick distribution) was used as a proxy for Lyme disease distribution and correlated with environmental or human-related data [e.g., forest cover, urban development; (, )]. Other studies considered several variables, but these distribution models were not applied to develop risk maps (, ). It should also be noted that despite only 20 incidences of human-related variables being included in these studies' models overall (Table 5), more recent research has begun to consider general public influence within their studies by including citizen science or “Google trends” to determine blacklegged tick and Lyme disease distributions (, , ). Other human demographic data can be acquired through public government agencies such as population densities and household incomes (, –). More personal and individualistic information can be acquired via questionnaire surveys (, , , ).
Table 5
| Paper ID | Anthropic variables used | References |
|---|---|---|
| MP9 | - Human population density | Lieske and Lloyd () |
| MP15 | - Human population size per county | Bisanzio et al. () |
| MP18 | - Household income | Little et al. () |
| MS3 | - Age and sex of Lyme disease patients | Tutt-Guerette et al. () |
| MS8 | - Human population density - Hiking behavior and Lyme awareness survey | Tadiri et al. () |
| MS12 | - Human population size per census division | Gasmi et al. () |
| MS17 | - “Lifestyle” categories based on surveys | Ozdenerol et al. () |
| MS18 | - Sociobehaviours; preventive behavior score, knowledge score, and risk perception score - Human population density | Bouchard et al. () |
| MS21 | - Behavioral risk factors - Reported tick exposure (by survey respondents) | Aenishaenslin et al. () |
| MS25 | - Human population density | Self et al. () |
| MS26 | - Human population density | Glavanakov et al. () |
| MS28 | - Human population density - Google trends for Lyme disease focused key words | Kutera et al. () |
| MS35 | - Human population density | Diuk-Wasser et al. () |
| MS36 | - Human population density | Larsen et al. () |
List of peer reviewed articles included in the literature review that incorporated anthropic variables in their risk assessments where a map was produced as a result.
Specific anthropic variables are described below along with their source.
There is a clear understanding that environmental factors heavily influence the distribution of diseases and their vectors (, , , ). However, very few studies consider that human related factors may also influence these distributions (, ). Further, human population growth, urbanization, and travel can affect vector-borne distributions (), particularly when human movement or land development influences animal movement, effectively altering the host dynamics of tick vectors (, ). This is concerning, as the risks related to Lyme disease directly affect human health, and human behavior can affect infection rates (, , , ). Predictive Map Type risk maps especially ignore how human-related factors may influence Lyme disease risks for humans, and this may be due to studies focusing on B. burgdorferi distribution rather than infection rates.
Overall, it was found that risk maps that focused on Lyme disease and blacklegged tick distributions differed across the types of variables used according to the study goal and intended use of the map. Because of these inconsistencies, it is difficult to compare and validate these models to accurately forecast Lyme disease risks geographically or temporally. At the same time, geospatial data related to anthropic factors can be difficult to acquire. Our results bring attention to the fact that there is no consistent “risk” variable or assessment across studies, likely because these studies tend to vary in specific objectives, despite the general intent of public health intervention. For this reason, risk maps should be scrutinized more thoroughly. As our knowledge on blacklegged ticks and B. burgdorferi increases, we must continually re-assess how risk models have predicted their geographic distributions over time. Differences in tick exposure patterns and Lyme disease risk is likely across regions and can depend on the scale and socioeconomic factors included in the assessment (). Future studies should consider improvements for forecasting these risks, as well as exploring risk assessments beyond comparison of blacklegged ticks and Lyme disease. Expanding the scope to other tick-borne diseases or co-infections of other bacteria including Babesia and Anaplasma sp. (, ) may demonstrate further patterns with spatial risk models.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
EF: conceptualization, methodology, formal analysis, investigation, data curation, writing—original draft, visualization, and project administration. MV: resources, writing—review and editing, project administration, and funding acquisition. VM: conceptualization, resources, writing—review and editing, supervision, project administration, and funding acquisition. All authors contributed to the article and approved the submitted version.
Funding
This project was funded through a Research Grant from the Institut Robert Sauvé en Santé du Travail (IRSST) Competitive Research Grant 2019-0001. VM was supported by funds from a NSERC Discovery Grant RGPIN-2017-03839.
Acknowledgments
We thank G. Joanisse and B. Chalghaf as well as the Institut Robert Sauvé en Santé du Travail (IRSST) for their support.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2023.1113024/full#supplementary-material
References
1.
SonenshineDE. 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. 10.3390/ijerph15030478
2.
RocklövJDubrowR. Climate change: an enduring challenge for vector-borne disease prevention and control. Nat Immunol. (2020) 21:479–83. 10.1038/s41590-020-0648-y
3.
ChalaBHamdeF. Emerging and re-emerging vector-borne infectious diseases and the challenges for control: a review. Front Public Health. (2021) 9:715759. 10.3389/fpubh.2021.715759
4.
Fichet-CalvetERogersDJ. Risk maps of lassa fever in West Africa. PLoS Negl Trop Dis. (2009) 3:e388. 10.1371/journal.pntd.0000388
5.
ZhaoXThanapongtharmWLawawirojwongSWeiCTangYZhouYet al. Malaria risk map using spatial multi-criteria decision analysis along yunnan border during the pre-elimination period. Am J Trop Med Hyg. (2020) 103:793–809. 10.4269/ajtmh.19-0854
6.
García-CarrascoJMMuñozAROliveroJSeguraMRealR. Mapping the risk for west nile virus transmission, Africa. Emerg Infect Dis. (2022) 28:777–85. 10.3201/eid2804.211103
7.
MiniganJNHagerHAPeregrineASNewmanJA. Current and potential future distribution of the American dog tick (Dermacentor variabilis, Say) in North America. Ticks Tick Borne Dis. (2018) 9:354–62. 10.1016/j.ttbdis.2017.11.012
8.
DranschDRotzollHPoserK. The contribution of maps to the challenges of risk communication to the public. Int J Digit Earth. (2010) 3:292–311. 10.1080/17538941003774668
9.
WallerLACarlinBP. Disease mapping. Chapman Hall CRC Handb Mod Stat Methods. (2010) 2010:217–43. 10.1201/9781420072884-c14
10.
JordanSJorgensenLMitterhoferH. Performing risk and the project: risk maps as mediating instruments. Manag Account Res. (2013) 24:156–74. 10.1016/j.mar.2013.04.009
11.
ThompsonCSpielmanAKrausePJ. Coinfecting deer-associated zoonoses: Lyme disease, babesiosis, and ehrlichiosis. Clin Infect Dis. (2001) 33:676–85. 10.1086/322681
12.
WormserGPDattwylerRJShapiroEDHalperinJJSteereACKlempnerMSet al. The clinical assessment, treatment, and prevention of Lyme disease, human granulocytic anaplasmosis, and babesiosis: clinical practice guidelines by the infectious diseases society of America. Clin Infect Dis. (2006) 43:1089–134. 10.1086/508667
13.
AucottJNYangTYoonIPowellDGellerSARebmanAWet al. Risk of post-treatment Lyme disease in patients with ideally-treated early Lyme disease: a prospective cohort study. Int J Infect Dis. (2022) 116:230–7. 10.1016/j.ijid.2022.01.033
14.
GrayJSKirsteinFRobertsonJNSteinJKahlO. Borrelia burgdorferi sensu lato in Ixodes ricinus ticks and rodents in a recreational park in south-western Ireland. Exp Appl Acarol. (1999) 23:717–29.
15.
SteereACCoburnJGlicksteinL. The emergence of Lyme disease. J Clin Invest. (2004) 113:1093–101. 10.1172/JCI21681
16.
SchwartzAMHinckleyAFMeadPSHookSAKugelerKJ. Surveillance for Lyme disease-United States, 2008–2015. MMWR Surveill Summ. (2017) 66:1–12. 10.15585/mmwr.ss6622a1
17.
RosenbergRLindseyNPFischerMGregoryCJHinckleyAFMeadPSet al. Vital Signs: trends in reported vectorborne disease case-United States and territories, 2004-2016. NMWR. (2018) 67:496–501. 10.15585/mmwr.mm6717e1
18.
LeightonPAKofflJKPelcatYLindsayLROgdenNH. Predicting the speed of tick invasion: an empirical model of range expansion for the Lyme disease vector Ixodes scapularis in Canada. J App Ecol. (2012) 49:457–64. 10.1111/j.1365-2664.2012.02112.x
19.
EisenRJEisenLOgdenNHBeardCB. Linkages of weather and climate with Ixodes scapularis and Ixodes pacificus (Acari: Ixodidae), enzootic transmission of Borrelia burgdorferi, and Lyme disease in North America. J Med Entomol. (2015) 53:250–61. 10.1093/jme/tjv199
20.
LiSGilbertLHarrisonPARounsevellMDA. Modelling the seasonality of lyme disease risk and the potential impacts of a warming climate within the heterogenous landscapes of Scotland. J R Soc Interface. (2016) 13:2016140. 10.1098/rsif.2016.0140
21.
GinsbergHSAlbertMAcevedoLDyerMCArsnoeIMTsaoJLet al. Environmental factors affecting survival of immature Ixodes scapularis and implications for geographical distribution of Lyme disease: the climate/behavior hypothesis. PLoS ONE. (2017) 12:e0168723. 10.1371/journal.pone.0168723
22.
BouchardCAenishaenslinCReesEEKoffiJKPelcatYRipocheMet al. Integrated social-behavioral and ecological risk maps to prioritize local public health responses to Lyme disease. Environ Health Perspect. (2018) 15:1048. 10.1289/EHP1943
23.
OstfeldRMillerMHazlerK. Causes and consequences of tick (Ixodes scapularis) burdens on white-footed mice (Peromyscus leucopus). J Mamm. (1996) 77:166–273. 10.2307/1382727
24.
OgdenNHLindsayLRMorshedMSockettPNArtsobH. The emergence of Lyme disease in Canada. CMAJ. (2009) 180:1221–4. 10.1503/cmaj.080148
25.
SimonJAMarrotteRRDesrosiersNFisetJGaitanJGonzalezAet al. Climate change and habitat fragmentation drive the occurrence of Borrelia burgdorferi, the agent of Lyme disease, at the northeastern limit of its distribution. Evol Appl. (2014) 7:750–64. 10.1111/eva.12165
26.
BurtisJCSullivanPLeviTOggenfussKFaheyTJOstfeldRSet al. The impact of temperature and precipitation on blacklegged tick activity and Lyme disease incidence in endemic and emerging regions. Parasit Vectors. (2016) 9:606. 10.1186/s13071-016-1894-6
27.
ClowKWeeseMRousseauJSJardineJCM. Microbiota of field-collected Ixodes scapularis and Dermacentor variabilis from eastern and southern Ontario, Canada. Ticks Tick Borne Dis. (2018) 9:235–44. 10.1016/j.ttbdis.2017.09.009
28.
McClureMDiuk-WasserMA. Climate impacts on blacklegged tick host-seeking behavior. Int J Parasit. (2019) 49:37–47. 10.1016/j.ijpara.2018.08.005
29.
KeiransJHutchesonJHDurdenLKlompenJSH. Ixodes (Ixodes) scapularis (Acari: Ixodidae): redescription of all active stages, distribution, hosts, geographical variation, and medical and veterinary importance. J Med Entomol. (1996) 33:297–318. 10.1093/jmedent/33.3.297
30.
HalseySAllanBMillerJ. The role of Ixodes scapularis, Borrelia burgdorferi and wildlife hosts in Lyme disease prevalence: a quantitative review. Ticks Tick Borne Dis. (2018) 9:1103–14. 10.1016/j.ttbdis.2018.04.006
31.
KhatchikianCPrusinskiMStoneMBackensonPBWangIFoleyEet al. Recent and rapid population growth and range expansion of the Lyme disease tick vector, Ixodes scapularis, in North America. Evolution. (2015) 69:1678–89. 10.1111/evo.12690
32.
KadaSMcCoyKBoulinierT. Impact of life stage-dependent dispersal on the colonization dynamics of host patches by ticks and tick-borne infectious agents. Parasit Vectors. (2017) 10:375. 10.1186/s13071-017-2261-y
33.
OgdenNHBigras-PoulinMO'CallaghanCJBarkerIKLindsayLRMaaroufAet al. 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. 10.1016/j.ijpara.2004.12.013
34.
PepinKMEisenRJMeadPSPiesmanJFishDHoenAGet al. Geographic variation in the relationship between human lyme disease incidence and density of infected host-seeking Ixodes scapularis nymphs in the Eastern United States. Am J Trop Med Hyg. (2012) 86:1062–71. 10.4269/ajtmh.2012.11-0630
35.
OgdenNHRadojevicMWuXDuvvuriVRLeightonPAWuJet al. Estimated effects of projected climate change on the basic reproductive number of the lyme disease vector Ixodes scapularis. Environ Health Perspect. (2014) 631–8. 10.1289/ehp.1307799
36.
Diuk-WasserMAVanAckerMCFernandezMP. Impact of land use changes and habitat fragmentation on the eco-epidemiology of tick-borne diseases. J Med Entomol. (2021) 58:1546–64. 10.1093/jme/tjaa209
37.
CullB. Monitoring trends in distribution and seasonality of medically important ticks in North America using online crowdsourced records from iNaturalist. Insects. (2022) 13:404. 10.3390/insects13050404
38.
OgdenNHSt-OngeLBarkerIKBrazeauSBigras-PoulinMCharronDFet al. Risk maps for range expansion of the Lyme disease vector, Ixodes scapularis, in Canada now and with climate change. Int J Health Geogr. (2008) 7:1–5. 10.1186/1476-072X-7-24
39.
KotchiSOBouchardCBrazeauSOgdenNH. Earth observation-informed risk maps of the Lyme disease vector Ixodes scapularis in Central and Eastern Canada. Remote Sens. (2021) 3:524–44. 10.3390/rs13030524
40.
KalaniTJSouthATalmadgeCLeiblerJWhittierCRosenbaumMet al. One map: using geospatial analysis to understand lead exposure across humans, animals, and the environment in an urban US city. One Health. (2021) 13:100341. 10.1016/j.onehlt.2021.100341
41.
BouchardCDumasABaronGBowserNLeightonPALindsayLRet al. Integrated human behavior and tick risk maps to prioritize Lyme disease interventions using a 'One Health' approach. Tick Tick Borne Dis. (2023) 14:102083. 10.1016/j.ttbdis.2022.102083
42.
OzdenerolEBingham-ByrneRMSebolyJD. The effects of lifestyle on the risk of lyme disease in the united states: evaluation of market segmentation systems in prevention and control strategies. Int J Environ Res Public Health. (2021) 18:12883. 10.3390/ijerph182412883
43.
PiacentinoJDSchwartzBS. Occupational risk of Lyme disease: an epidemiological review. Occup Environ Med. (2002) 59:75–84. 10.1136/oem.59.2.75
44.
Gutierrez-HoffmannMGO'MeallyRNColeRNTiniakouEDarrahESoloskiMJet al. Borrelia burgdorferi-induced changes in the class II self-immunopeptidome displayed on HLA-DR molecules expressed by dendritic cells. Front Med. (2020) 7:568. 10.3389/fmed.2020.00568
45.
TadiriCAinsworthNDe BonoNGavinSLiJMilbersKet al. Assessment of human health risk for Lyme disease in a peri urban park in southern Quebec. McGill Sci Undergrad Res J. (2011) 6:56–61. 10.26443/msurj.v6i1.96
46.
AenishaenslinCCharlandKBowserNPerez-TrejoEBaronGMilordFet al. Behavioral risk factors associated with reported tick exposure in a Lyme disease high incidence region in Canada. BMC Public Health. (2022) 22:807. 10.1186/s12889-022-13222-9
47.
FixADPeñaCAStricklandT. Racial differences in reported Lyme disease incidence. Am J Epidemiol. (2000) 152:756–9. 10.1093/aje/152.8.756
48.
DongYZhouGCaoWXuXZhangYJiZet al. Global seroprevalence and sociodemographic characteristics of Borrelia burgdorferi sensu lato in human populations: a systematic review and meta-analysis. BMJ Global Health. (2022) 7:e007744. 10.1136/bmjgh-2021-007744
49.
XuGWielstraBRichSM. Northern and southern blacklegged (deer) ticks are genetically distinct with different histories and Lyme spirochete infection rates. Sci Rep. (2020) 10:10289. 10.1038/s41598-020-67259-0
50.
GinsbergHSHicklingGJBurkeRLOgdenNHBeatiLLeBrunRAet al. Why Lyme disease is common in the northern US, but rare in the south: the roles of host choice, host-seeking behavior, and tick density. PLoS Biol. (2021) 19:e3001396. 10.1371/journal.pbio.3001396
51.
ParéGKitsiouS. Chapter 9: Methods for literature reviews. In:LauFKuziemskyC, editors. Handbook of eHealth Evaluation: An Evidence-based Approach. Victoria, BC: University of Victoria (2017). Available online at: https://www.ncbi.nlm.nih.gov/books/NBK481583/ (accessed August 1, 2023).
52.
KjaerLJSolengAEdgarKSLindstedtHEHPaulsenKMAndreassenAKet al. Predicting and mapping human risk of exposure to Ixodes ricinus nymphs using climatic and environmental data, Denmark, Norway and Sweden, 2016. Euro Surveill. (2019) 24:1800101. 10.2807/1560-7917.ES.2019.24.9.1800101
53.
MoherDShamseerLClarkeMGhersiDLiberatiAPetticrewMet al. Preferred reporting items for systematic review and meta-analysis protocols (PRISMA-P) 2015 statement. Syst Rev. (2015) 4:1. 10.1186/2046-4053-4-1
54.
ClowKMOgdenNHLindsayLRMichelPPearlDLJardineCM. Distribution of ticks and the risk of lyme disease and other tick-borne pathogens of public health significance in Ontario, Canada. Vector Borne Zoonotic Dis. (2016) 16:215–22. 10.1089/vbz.2015.1890
55.
ClowKMLeightonPAOgdenNHLindsayLRMichelPPearlDLet al. Northward range expansion of Ixodes scapularis evident over a short timescale in Ontario, Canada. PLoS ONE. (2017) 12: e0189393. 10.1371/journal.pone.0189393
56.
KoffiJKLeightonPAPelcatYTrudelLLindsayLRMilordFet al. Passive surveillance for I. scapularis ticks: enhanced analysis for early detection of emerging Lyme disease risk. J Med Entomol. (2012) 49:400–9. 10.1603/ME11210
57.
Tokarska-RodakMPlewikDMichalskiAJKolodziejMMelgiesAPanczukAet al. Serological surveillance of vector-borne and zoonotic diseases among hunters in eastern Poland. J Vector Borne Dis. (2016) 53:355–61.
58.
GasmiSOgdenNHLindsayLRBurnsSFlemingSBadcockJet al. Surveillance for Lyme disease in Canada: 2009-2015. Can Commun Dis Rep. (2017) 43:194–9. 10.14745/ccdr.v43i10a01
59.
ChiltonNBCurryPSLindsayLRRochonKLysykTJDergousoffSJet al. Passive and active surveillance for Ixodes scapularis (Acari: Ixodidae) in Saskatchewan, Canada. J Med Entomol. (2020) 57:156–63. 10.1093/jme/tjz155
60.
RStudio Team,. RStudio: Integrated Development Environment for R. Boston, MA: RStudio, PBC (2021). Available online at: http://www.rstudio.com/ (accessed August 12, 2023).
61.
OksanenJBlanchetFGFriendlyMKindtRLegendrePMcGlinnDet al. vegan: Community Ecology Package. R package version 2.5-7. (2020). Available online at: https://CRAN.R-project.org/package=vegan (accessed August 12, 2023).
62.
GuisanAHarrellFE. Ordinal response regression models in ecology. J Veg Sci. (2000) 11:617–26. 10.2307/3236568
63.
ChristensenRHB,. Ordinal - Regression Models for Ordinal Data. R package version 2019.12-10. (2019). Available online at: https://CRAN.R-project.org/package=ordinal (accessed August 12, 2023).
64.
FischhoffIRKeesingFOstfeldRS. Risk factors for bites and diseases associated with black-legged ticks: a meta-analysis. Am J Epidemiol. (2019) 188:1742–50. 10.1093/aje/kwz130
65.
ShararehNBehlerRPRoomeABShepherdJGarrutoRMSabounchiNSet al. Risk factors of Lyme disease: an intersection of environmental ecology and systems science. Healthcare. (2019) 7:66. 10.3390/healthcare7020066
66.
Diuk-WasserMAGatewoodAGCortinasMRYaremych-HamerSTsaoJKitronUet al. Spatiotemporal patterns of host-seeking Ixodes scapularis nymphs (Acari: Ixodidae) in the United States. J Med Entomol. (2006) 43:166–76. 10.1093/jmedent/43.2.166
67.
GrayJSDautelHEstrada-PenaAKahlOLindgrenE. Effects of climate change on ticks and tick-borne diseases in Europe. Interdis Perspect Infect Dis. (2009) 2009:1–12. 10.1155/2009/593232
68.
EndoAAmarasekareP. Predicting the spread of vector-borne diseases in a warming world. Front Ecol Evol. (2022) 10:758277. 10.3389/fevo.2022.758277
69.
SlatculescuAMDuguayCOgdenNHSanderBDesjardinsMWilliam CameronDet al. Spatiotemporal trends and socioecological factors associated with Lyme disease in eastern Ontario, Canada from 2010-2017. BMC Public Health. (2022) 22:736. 10.1186/s12889-022-13167-z
70.
TalbotBSlatculescuAThickstunCRKoffiJKLeightonPAMcKayRet al. Landscape determinants of density of blacklegged ticks, vectors of Lyme disease, at the northern edge of their distribution in Canada. Sci Rep. (2019) 9:16652. 10.1038/s41598-019-50858-x
71.
Feria-ArroyoTPCastro-ArellanoIGordillo-PerezGCavazosALVargas-SandovalMGroverAet al. Implications of climate change on the distribution of the tick vector Ixodes scapularis and risk for Lyme disease in the Texs-Mexico transboundary region. Parasit Vectors. (2014) 7:199. 10.1186/1756-3305-7-199
72.
LarsenAEMacDonaldAJPlantingaAJ. Lyme disease risk influences human settlement in the wildland-urban interface: evidence from a longitudinal analysis of counties in the Northeastern United States. Am J Trop Med Hyg. (2014) 91:747–55. 10.4269/ajtmh.14-0181
73.
BaldwinHLandesmanWJBorgmann-WinterBAllenD. A geographic information system approach to map tick exposure risk at a scale for public health intervention. J Med Entomol. (2022) 59:162–72. 10.1093/jme/tjab169
74.
LieskeDJLloydVK. Combining public participatory surveillance and occupancy modelling to predict the distributional response of Ixodes scapularis to climate change. Ticks Tick Borne Dis. (2018) 9:695–706. 10.1016/j.ttbdis.2018.01.018
75.
LittleEAHAndersonJFStaXord IIIKCEisenLEisenRJMolaeiG. Predicting spatiotemporal patterns of Lyme disease incidence from passively collected surveillance data for Borrelia burgdorferi sensu lato-infected Ixodes scapularis ticks. Ticks Tick Borne Dis. (2019) 10:970–80. 10.1016/j.ttbdis.2019.04.010
76.
BisanzioDFernandezMPMartelloEReithingerRDiuk-WasserMA. Current and future spatiotemporal patterns of Lyme disease reporting in the Northeastern United States. JAMA Network Open. (2020) 3:e200319. 10.1001/jamanetworkopen.2020.0319
77.
FinchCAl-DamlujiMSKrausePJNiccolaiLSteevesTFolsom O'KeefeCet al. Integrated assessment of behavioral and environmental risk factors for Lyme disease infection on Block Island, Rhode Island. PLoS ONE. (2014) 9:e84758. 10.1371/journal.pone.0084758
78.
EisenLEisenRJ. Critical evaluation of the linkage between tick-based risk measures and the occurrence of Lyme Disease cases. J Med Entomol. (2018) 53:1050–62. 10.1093/jme/tjw092
79.
FèvreEMBronsvoortBMHamiltonKACleavelandS. Animal movements and the spread of infectious diseases. Trends Microbiol. (2006) 14:125–31. 10.1016/j.tim.2006.01.004
80.
MuehlenbeinMP. Human-wildlife contact and emerging infectious diseases. Hum.-Environ Interact. (2012) 1:79–94. 10.1007/978-94-007-4780-7_4
81.
MacDonaldAJLarsenAEPlantingaAJ. Missing the people for the trees: identifying coupled natural–human system feedbacks driving the ecology of Lyme disease. J Appl Ecol. (2019) 56:354–64. 10.1111/1365-2664.13289
82.
TelfordS. Human Behaviour Trumps Entomological Risk. CABI Climate Change Series. Wallingford: CABI International. 10.1079/9781789249637.0040
83.
NelderMPRussellCBSheehanNJSanderBMooreSLiYet al. Human pathogens associated with the blacklegged tick Ixodes scapularis: a systematic review. Parasit Vectors. (2016) 9:265. 10.1186/s13071-016-1529-y
84.
Tutt-GueretteMYuanMSzarozDMcKinnonBKestensYGuillotCet al. Modelling spatiotemporal patterns of Lyme disease emergence in Quebec. Int J Environ Res Public Health. (2021) 18:9669. 10.3390/ijerph18189669
85.
GasmiSOgdenNHRipocheMLeightonPALindsayRLNelderMPet al. Detection of municipalities at-risk of Lyme disease using passive surveillance of Ixodes scapularis as an early signal: a province specific indicator in Canada. PLoS ONE. (2019) 14:e0212637. 10.1371/journal.pone.0212637
86.
SelfSCWLiuYNordoneSKYabsleyMJWaldenHSLundRBet al. Canine vector-borne disease: mapping and the accuracy of forecasting using big data from the veterinary community. Anim Health Res Rev. (2019) 20:47–60. 10.1017/S1466252319000045
87.
GlavanakovSWhiteDJCaracoTLapenisARobinsonGRSzymanskiBKet al. Lyme disease in New York State: spatial pattern at a regional scale. AM J Trop Med Hyg. (2001) 65:538–45. 10.4269/ajtmh.2001.65.538
88.
KuteraMBerkeOSobkowichK. Spatial epidemiological analysis of Lyme disease in southern Ontario utilizing Google Trends searches. EHR. (2022) 64:105–10. 10.5864/d2021-025
89.
PorterWTMotykaPJWacharaJBarrandZAHmoodZMcLaughlinMet al. Citizen science informs human-tick exposure in the Northeastern United States. Int J Health Geogr. (2019) 18: 9. 10.1186/s12942-019-0173-0
90.
HartCEBhaskarJRReynoldsEHermanceMEarlMMahoneyMet al. Community engaged tick surveillance and tickMAP as a public health tool to track the emergence of ticks and tick-borne diseases in New York. PLOS Global Public Health. (2022) 2:e0000215. 10.1371/journal.pgph.0000215
Summary
Keywords
blacklegged ticks, data synthesis, human-related, Lyme disease, risk assessment, risk map
Citation
Fellin E, Varin M and Millien V (2023) Risky business: human-related data is lacking from Lyme disease risk models. Front. Public Health 11:1113024. doi: 10.3389/fpubh.2023.1113024
Received
01 February 2023
Accepted
09 October 2023
Published
03 November 2023
Volume
11 - 2023
Edited by
Miguel Angel Garcia Bereguiain, University of the Americas, Ecuador
Reviewed by
Vett Lloyd, Mount Allison University, Canada; Alberto Orlando, Espiritu Santo University, Ecuador
Updates
Copyright
© 2023 Fellin, Varin and Millien.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Erica Fellin erica.fellin@mail.mcgill.ca
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.