1. Introduction
In the complex realm of public health, Lyme disease stands as a formidable challenge, annually impacting approximately 30,000 individuals, with an estimated 476,000 cases in the United States (
http://www.gov.cdc/lyme). The causative agent,
Borrelia burgdorferi sensu stricto, a tick-borne spirochete, is transmitted through the bite of the blacklegged tick (
Ixodes scapularis), colloquially known as the deer tick. This spatiotemporal surge in Lyme disease incidence and geographic spread is particularly concentrated in the Northeast and Upper Midwest, driven by a dynamic interplay of climate, host species populations, and various human-related factors [
1,
2].
Contrary to the prevailing perception that Lyme disease risk is mitigated in highly urbanized settings, recent revelations challenge this notion, unraveling a new narrative, especially in expansive metropolises like New York City. The density of
I. scapularis adults and nymphs, coupled with the infection prevalence of
B. burgdorferi within these ticks, mirrors that of nonurban residential and recreational areas in the highly endemic coastal Northeast [
3,
4]. This paradigm-shifting study seeks to comprehensively evaluate Lyme disease exposure risk in 24 regional parks nestled within the iconic Staten Island of New York City, probing into the influence of inner-city greenspaces in the vicinity of densely populated neighborhoods.
Dispelling the conventional belief that most New York City Lyme disease cases stem from travel outside the urban enclave, a notable upswing in white-tailed deer populations in Staten Island has prompted localized interventions, including innovative sterilization programs [
5]. As the number of ticks escalates in parks and locally acquired Lyme disease cases surge, the impact of neighborhood-level characteristics, such as greenspaces and the diverse fauna inhabiting ticks in Staten Island, becomes increasingly apparent [
6].
Building on the foundation of prior research that underlines the critical importance of understanding tick-borne disease risks at the neighborhood scale, our study pioneers the utilization of the acarological risk index [
7,
8,
9,
10,
11]. This index, rooted in the infection rates of
I. scapularis nymphs by
B. burgdorferi, emerges as a beacon of predictability and correlation for Lyme disease risk. Drawing from meticulously conducted monthly tick surveys spanning from 2019 to 2022 in Staten Island's 24 main parks, our acarological indices transcend temporal and spatial boundaries, offering a nuanced understanding of Lyme disease transmission risks. These indices stand as invaluable tools, empowering us to alert park patrons about the ever-changing levels of Lyme disease transmission risk throughout the year. In the tapestry of public health, this study unfurls a narrative of foresight, awareness, and informed intervention in the face of a mounting health challenge.
2. Material and Methods
The investigation spanned from 2019 to 2022, encompassing monthly assessments conducted across 24 parks in Staten Island. The study focused on two developmental stages of ticks: nymphs and adults. To quantify the acarological risk index (ARI), the collection figures and testing outcomes for both stages were combined.
2.1. Tick Collections
Every year, the sampling regimen targeted twenty-four primary parks in Staten Island (refer to
Table 1). Standard flagging procedures were employed for tick collection. This method involves dragging a 1 square meter section of white corduroy cloth across the leaf litter and low-lying vegetation. Sampling areas were strategically flagged along publicly accessible pathways, within a 5-meter radius on either side of the woodland's edge. Following every 10 steps, the cloth underwent visual inspection, and all identified ticks were promptly removed and stored in labeled vials preserved with 70% ethanol.
Concurrently, pertinent information such as the date, sampling time, current temperature, relative humidity, and other relevant weather conditions were meticulously recorded during the flagging process. Subsequently, ticks were subjected to species, gender, and stage identification using morphological keys [
12]. The identified ticks were then stored in absolute ethanol and dispatched to the laboratory for comprehensive pathogen testing.
2.2. DNA Isolation and PCR Assays
Blacklegged ticks, encompassing nymphs and adults, underwent transportation to the Public Health Laboratory at the Department of Health and Mental Hygiene in New York City. Upon arrival, individual or pooled ticks were meticulously placed in 1.5-2.0 mL microcentrifuge tubes, where they were consistently maintained on dry ice or ice post-removal from the freezer to prevent degeneration of DNA and viral RNA. The extracted total nucleic acid samples derived from these ticks were then utilized in real-time RT-PCR. Both the homogenized ticks and nucleic acid underwent storage at -70°C after examination to prevent degradation during freezing/thawing. The Public Health Laboratory employed the QuantStudio Dx instrument for a multiplex real-time RT-PCR assay capable of detecting five pathogens from individual or pooled ticks [
13]. Executed by competent personnel equipped with appropriate biosafety measures, this nucleic acid amplification assay incorporated five sets of oligonucleotide primers and hydrolysis probes. The target for
B. burgdorferi was a plasmid-borne gene exclusive to Lyme borreliosis-borrelia. Fluorescently labeled probes bound to amplified DNA fragments, and the QuantStudio Dx instrument monitored the fluorescent signal intensity during each PCR cycle. Target amplification was recorded through the observation of an increase in fluorescence over time relative to the background signal.
2.3. Data Analysis
The acarological risk index (ARI) was computed for each park where tick collections took place. Representing the number of nymphs and adults of
I. scapularis infected with
B. burgdorferi collected per minute of flag sampling [
14], the ARI values formed a basis for subsequent statistical analyses. Analyses of variance (ANOVA) were conducted, employing SAS EG 7.1 (SAS Institute, Cary, NC), to compare monthly and yearly average tick densities, average tick infection rates, and average ARI values across the 24 parks. This analysis entailed testing the hypothesis that each sample is drawn from the same underlying probability distribution, against the alternative hypothesis suggesting variations in underlying probability distributions among the samples. Statistical significance was attributed to p-values less than 0.05. Spatial mapping of park locations and spatial ARI interpolation were executed using Arc/GIS (version 10.6.1, Esri, INC).
3. Results
3.1. Tick Collections and Testing for B. burgdorferi
Between 2019 and 2022, a comprehensive survey yielded a total of 1,139
I. scapularis adults and nymphs across the 24 parks in Staten Island (refer to
Table 2). Among the sampled ticks, 662 (58.1%) were nymphs, while 477 (41.8%) were adults. Notably, 330 ticks, comprising 29.0% of all nymphs and adults, were found to be infected with
B. burgdorferi. The yearly average infection rates across these 24 parks exhibited minimal fluctuations, ranging from a low of 26.84% in 2020 to a peak of 34.21% in 2019, resulting in an overall four-year average of 28.97%.
Further analysis revealed a notable disparity between the infection rates of adults and nymphs, with the percentage of infected adults notably higher at 39.4%, compared to 21.5% in nymphs. Interestingly, the infection rates among male and female adults were remarkably close, standing at 40.3% and 38.3%, respectively.
Delving into specific parks, Wolfe's Pond Park emerged with the highest number of tested ticks at 226, followed by 117 in Clay Pit Pond Park, and 109 in Fair View Park. Meanwhile, High Rock Park, Willowbrook Park, La Tourette Golf Course, and Clove Lake Park reported lower numbers, ranging from 50 to 100. The remaining 17 parks displayed the lowest tested numbers, all falling below 50 ticks.
3.2. Tick Densities and Acarological Risk Index (ARI)
The densities of ticks, as detailed in
Table 3 and illustrated in
Figure 1, exhibited significant variations across parks and years. In 2021, a noteworthy increase in ticks was observed, with a yearly average of 0.5 ticks per minute, surpassing the other three years (df=3, F=2.72, p<0.05). Throughout the period from 2019 to 2022, the monthly average infection rates across all 24 parks remained relatively stable, fluctuating modestly between 8% and 10%. However, a more granular examination of individual park infection rates unveiled substantial diversity in infectivity, ranging from 0% to 100% on a monthly basis.
Regarding the number of ticks tested, 71.4% of the 24 parks reported no ticks collected or tested in specific months between 2019 and 2022. Additionally, 24.2% of parks had monthly tick test numbers ranging from 1 to 5, while 3.2% reported numbers from 6 to 10, and only 1.2% recorded figures between 11 and 40 ticks. Although not statistically significant, 2021 stood out with the highest ARI of 0.14, followed by 0.12 in 2022, and comparatively lower values in 2019 and 2020.
Wolfe’s Pond Park, High Rock Park, Clay Pit Pond Park, Clove Lake Park, and Fair View Park (refer to
Figure 1) emerged as focal points with both elevated tick collections (ranging from 0.39 to 1.54) and ARI values (ranging from 0.16 to 0.53). Notably, these parks also exhibited relatively higher infection rates, spanning from 10% to 26%.
In the seasonal analysis (refer to
Figure 2), notably elevated acarological risk indices (ARIs) were identified in June (0.36) and November (0.21) (df=11, F=5.08, p<0.01). This finding aligns with the months of peak tick collection, registering 1.36 and 0.54 ticks per minute in June and November, respectively. Generally, a heightened tick collection trend was observed from March to July and October to November (df=11, F=6.15, p<0.01). However, the ability to detect infectivity persisted throughout the year, given the consistent monthly tick collections (df=11, F=5.77, p<0.01), except for an exception in September, likely attributed to a minimal number of ticks collected during that period.
3.3. Spatial Distribution of Acarological Risk Index (ARI)
Figure 3 illustrates the spatial distribution of Lyme disease transmission risk based on the classification of acarological risk indices (ARIs) across 24 parks. Wolfe’s Pond Park stands out with the highest risk level (ARI = 0.53), succeeded by 13 parks at a medium risk level (High Rock Park, Clay Pit Pond Park, Clove Lake Park, etc., 0.07<ARI<0.18), while the remaining 10 parks (Ocean Breeze Park, Reeds Basket Willow Swamp Park, etc.) exhibit the lowest risk level (ARI<0.07). Notably, four parks (Aesop Park, Lemon Creek Park, Sailor Snag Harbor, and Silver Lake Park) record ARIs as 0, indicating either no ticks were collected or no
B. burgdorferi was detected in these locations.
Leveraging the inverse distance weighting (IDW) model of the Spatial Analyst within Arc/GIS for tick ARI distribution interpolation in Staten Island reveals distinct risk patterns. Southern Staten Island, particularly around Wolfe’s Pond Park, Clay Pit Pond Park, and Fair View Park, emerges as a hotspot with elevated Lyme disease transmission risk. The central part of Staten Island exhibits a medium-level risk, while areas surrounding High Rock Park and Clove Lake Park may experience heightened activity. Conversely, the northeastern and most southwestern regions, characterized by reduced natural areas and higher human population density, demonstrate the lowest risk.
4. Discussion
The results of our investigation
unveil significant variations in blacklegged tick collections, for both nymphs and adults, across various parks and diverse temporal dimensions in Staten Island. The elevated tick collections and infection rates observed in Wolfe’s Pond Park, High Rock Park, Clay Pit Pond Park, Clove Lake Park, and Fair View Park, particularly in recent years, highlight the persistent challenges posed by Lyme disease [
15]. Infections were widespread, affecting nymphs and adults from 20 of the 24 surveyed parks (83.3%), suggesting a heightened risk of Lyme disease transmission throughout the year, except for September due to limited tick collection and testing [
16]. The continuous high infectivity of nymphs, despite a lower infection rate compared to adults, sheds light on the observed surge in locally acquired Lyme disease cases [
17].
Traditionally, studies in the United States primarily focused on low-intensity residential and forested areas to identify factors contributing to human Lyme disease infection [
18,
19,
20]. However, recent attention has shifted towards urban Lyme disease infection risks, challenging the notion that tickborne disease risk is confined to suburban and natural settings [
15]. Our findings advocate for a comprehensive risk assessment that integrates various methods, encompassing all tick stages, and considers human behavior and habits.
Our study highlights the imperative for ongoing assessment and management of Lyme disease risk in Staten Island's greenspaces. Spatial distribution analysis indicates higher transmission risks in southern Staten Island, emphasizing the need for targeted interventions and public awareness campaigns [
21]. However, the limited locally acquired cases and tick surveillance efforts hinder a detailed exploration of the relationship between acarological ARIs and human Lyme disease risk in Staten Island.
ARIs facilitate a meaningful comparison among parks characterized by substantial tick populations, yet continuous monitoring proves indispensable for the 17 parks with lower tick densities to discern population development trends [
17,
22]. Our findings not only function as a critical alert for the implementation of signage in parks and the guidance of tick control measures but also advocate for heightened public awareness. We propose the strategic installation of tick warning signs in City and State parks within Staten Island, particularly in areas designated with medium to high ARI values. Although the reduction of tick populations is generally deemed promising, our study underscores the necessity of evaluating outcomes for individuals, such as disease incidence or tick encounters, to gauge the efficacy of tick reduction methods [
21,
23]. The collaborative endeavors of our research team significantly contribute to the ongoing discourse on Lyme disease, aiming for well-informed public health measures and community safety. Furthermore, our study establishes a crucial baseline for tracking the spread of infection from areas endemic in 2019-2022.
Author Contributions
Waheed Bajwa, Conceptualization and Visualization, Project administration, Writing and reviewing; Liyang Zhou, Writing, reviewing, project administration; Leonid Tsynman, tick surveillance; Kamesan Kanapathipillai, specimen preparation for pathological testing; Zahir Shah, field supervision, specimen preparation for pathological testing. All authors have read and agreed to the published version of the manuscript.
Data Availability Statement
The original contributions presented in the study are included in the article; further inquiries can be directed to the corresponding author.
Conflicts of Interest
The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.
References
- Fouque, F.; Reeder, J.C. Impact of past and on-going changes on climate and weather on vector-borne diseases transmission: a look at the evidence. Infect. Dis. Poverty 2019, 8, 51. [Google Scholar] [CrossRef] [PubMed]
- Tran, T.; Prusinski, M. A.; White, J. L.; Falco, R. C.; Kokas, J.; Vinci, V.; Gall, W. K.; Tober, K. J.; Haight, J.; Oliver, J.; Sporn, L. A.; Meehan, L.; Banker, E.; Backenson, P. B.; Jensen, S. T.; Brisson, D. Predicting spatio-temporal population patterns of Borrelia burgdorferi, the Lyme disease pathogen. J. Appl. Ecol. 2022, 59, 2779–2789. [Google Scholar] [CrossRef] [PubMed]
- Kilpatrick, A.M.; Dobson, A.D.M.; Levi, T.; Salkeld, D.J.; Swei, A.; Ginsberg, H.S.; Kjemtrup, A.; Padgett, K.A.; Jensen, P.; Fish, D.; et al. Lyme disease ecology in a changing world: consensus, uncertainty and critical gaps for improving control. Philos. Trans. R. Soc. B: Biol. Sci. 2017, 372, 20160117. [Google Scholar] [CrossRef] [PubMed]
- Alison, E. S.; Anna, B. M.; Tashi, B.; Shannon, L.T.; Sara, E. M.; Matthew, Q. D.; Joseph, E. D.; Thomas, W. S. Lyme disease risk of exposure to blacklegged ticks (Acari: Ixodidae) infected with Borrelia burgdorferi (Spirochaetales: Spirochaetaceae) in Pittsburgh regional parks. J. Med. Entomol. 2020, 57, 273–280. [Google Scholar]
- Jason, R.B.; Paul, D.C.; Evan, G.C.; Anthony, J.D. Sterilization as an alternative deer control technique: a review. Human-Wildlife interactions 2012, 6, 273–282. [Google Scholar]
- New Yor k City Department of Health and Mental Hygiene. Preventing, diagnosing, and managing tick-borne diseases. City Health Information 2020, 39, 12–18.
- Ilya, R. F.; Sarah, E. B.; Felicia, K.; Richard, S. O. Systematic review and meta-analysis of tick-borne disease risk factors in residential yards, neighborhoods, and beyond. BMC Inf. Dis. 2019, 19, 861. [Google Scholar]
- Mather, T.N.; Nicholson, M.C.; Donnelly, E.F.; Matyas, B.T. Entomologic Index for Human Risk of Lyme Disease. Am. J. Epidemiology 1996, 144, 1066–1069. [Google Scholar] [CrossRef] [PubMed]
- Kirby, C. S.; Matthew, L. C.; Louis, A. M.; Starrhope, E.; Patricia, A. M. Temporal correlations between tick abundance and prevalence of ticks infected with Borrelia burgdorferi and increasing incidence of Lyme disease. J. Clin. Micro. 1998, 36, 1240–1244. [Google Scholar]
- Parasuraman, B.; Pichai, K.; Rajagopal, T.P.; Sivsankaran, S.; Subramaniam, S.; Mahaligam, K. Study on entomological surveillance and its significance during a dengue outbreak in the district of Tirunelveli in Tamil Nadu, India. Osong Pub. Health Res. Perspect. 2013, 4, 152–158. [Google Scholar]
- Milena, K.; Novica, S.; Srdjan, L. Prevalence of Borrelia burgdorferi sensu lato in Ixodes ricinus ticks and assessment of entomological risk index at localities in Belgrade. Vojnosanit Pregl. 2016, 73, 817–824. [Google Scholar]
- Keirans, J.E.; Hutcheson, H.J.; Durden, L.A.; Klompen, J. Ixodes (Ixodes) scapularis (Acari: Ixodidae): Redescription of all Active Stages, Distribution, Hosts, Geographical Variation, and Medical and Veterinary Importance. J. Med Èntomol. 1996, 33, 297–318. [Google Scholar] [CrossRef] [PubMed]
- Tokarz, R.; Tagliafierro, T.; Cucura, D.M.; Rochlin, I.; Sameroff, S.; Lipkin, W.I. Detection of Anaplasma phagocytophilum, Babesia microti, Borrelia burgdorferi, Borrelia miyamotoi, and Powassan Virus in Ticks by a Multiplex Real-Time Reverse Transcription-PCR Assay. mSphere 2017, 2. [Google Scholar] [CrossRef] [PubMed]
- Milena, K.; Novica, S.; Srdjan, L. Prevalence of Borrelia burgdorferi sensu lato in Ixodes ricinus ticks and assessment of entomological risk index at localities in Belgrade. Vojnosanit Pregl. 2016, 73, 817–824. [Google Scholar]
- Adams, D.A.; Thomas, K.R.; Jajosky, R. Summary of notifiable infectious diseases and conditions – United States. MMWR Morb. Mortal. Wkly. Rep. 2016, 63, 1–152. [Google Scholar] [CrossRef] [PubMed]
- Diuk-Wasser, M.A.; Brinkerhoff, R.; Kitron, U.; Fish, D.; Melton, F.; Cislo, P.; Tsao, J.I.; Rowland, M.; Barbour, A.G.; Vourc'H, G.; et al. Human Risk of Infection with Borrelia burgdorferi, the Lyme Disease Agent, in Eastern United States. Am. J. Trop. Med. Hyg. 2012, 86, 320–327. [Google Scholar] [CrossRef] [PubMed]
- Eisen, L.; Eisen, R.J. Critical Evaluation of the Linkage Between Tick-Based Risk Measures and the Occurrence of Lyme Disease Cases: Table 1. J. Med Èntomol. 2016, 53, 1050–1062. [Google Scholar] [CrossRef]
- Meredith, C.V.; Eliza, A.H.L.; Goudarz, M.; Waheed, I.B.; Maria, A.D. Enhancement of risk for Lyme disease by landscape connectivity, New York, New York, USA. Emerg. Inf. Dis. 2019, 25, 1137–1143. [Google Scholar]
- Simon, J. A.; Marrotte, R. R.; Desrosiers, N.; Fiset, J.; Gaitan, J.; Gonzalez, A.; Koffi, J. K.; Lapointe, F.J.; Leighton, P. A.; Lindsay, L. R.; Logan, T.; Milord, F.; Ogden, N. H.; Rogic, A.; Roy-Dufresne, E.; Suter, D.; Tessier, N.; Millien, V. Climate change and habitat fragmentation drive the occurrence of Borrelia burgdorferi, the agent of Lyme disease, at the northeastern limit of its distribution. Evo. App. 2014, 7, 750–764. [Google Scholar] [CrossRef]
- Johnson, T.L.; Bjork, J.K.H.; Neitzel, D.F.; Dorr, F.M.; Schiffman, E.K.; Eisen, R.J. Habitat Suitability Model for the Distribution ofIxodes scapularis(Acari: Ixodidae) in Minnesota. J. Med Èntomol. 2016, 53, 598–606. [Google Scholar] [CrossRef]
- Hinckley, A.F.; Meek, J.I.; Ray, J.A.; Niesobecki, S.A.; Connally, N.P.; Feldman, K.A.; Jones, E.H.; Backenson, P.B.; White, J.L.; Lukacik, G.; Kay, A.B.; Miranda, W.P.; Mead, P.S. Effectiveness of residential acaricides in reducing tick infestations in the United States. J. Med. Entomol. 2014, 51, 880–889. [Google Scholar]
- Donal, B.; Maria, P. F.; Elisa, M.; Richard, R.; Maria, A. D. Current and future spatiotemporal patterns of Lyme disease reporting in the Northeastern United States. JAMA Net. Opn. 2020, 3, 1–12. [Google Scholar]
- Felicia, K.; Stacy, M.; William, B.; Shannon, D.; Andrew, S.E.; Ilya, R.F.; Alison, F.H.; Sarah, A.H.; Fiona, K.; Jennifer, P.; Ashley, P.; Marissa, T.; Richard, S.O. Effects of tick-control interventions on tick abundance, human encounters with ticks, and incidence of tickborne diseases in residential neighborhoods, New York, USA. Emerg. Inf. Dis. 2022, 28, 957–966. [Google Scholar]
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).