- Open Access
A 117-year retrospective analysis of Pennsylvania tick community dynamics
© The Author(s) 2019
- Received: 12 November 2018
- Accepted: 20 April 2019
- Published: 29 April 2019
Tick-borne diseases have been increasing at the local, national, and global levels. Researchers studying ticks and tick-borne diseases need a thorough knowledge of the pathogens, vectors, and epidemiology of disease spread. Both active and passive surveillance approaches are typically used to estimate tick population size and risk of tick encounter. Our data consists of a composite of active and long-term passive surveillance, which has provided insight into spatial variability and temporal dynamics of ectoparasite communities and identified rarer tick species. We present a retrospective analysis on compiled data of ticks from Pennsylvania over the last 117 years.
We compiled data from ticks collected during tick surveillance research, and from citizen-based submissions. The majority of the specimens were submitted by citizens. However, a subset of the data was collected through active methods (flagging or dragging, or removal of ticks from wildlife). We analyzed all data from 1900–2017 for tick community composition, host associations, and spatio-temporal dynamics.
In total there were 4491 submission lots consisting of 7132 tick specimens. Twenty-four different species were identified, with the large proportion of submissions represented by five tick species. We observed a shift in tick community composition in which the dominant species of tick (Ixodes cookei) was overtaken in abundance by Dermacentor variabilis in the early 1990s and then replaced in abundance by I. scapularis. We analyzed host data and identified overlaps in host range amongst tick species.
We highlight the importance of long-term passive tick surveillance in investigating the ecology of both common and rare tick species. Information on the geographical distribution, host-association, and seasonality of the tick community can help researchers and health-officials to identify high-risk areas.
- Passive surveillance
- Museum collections
- Community composition
Surveillance data collected over multiple decades may reveal spatio-temporal changes in ectoparasite communities . Data such as spatial distribution and occurrence of both abundant and rare species of ticks can be correlated with land use (e.g. habitat loss, fragmentation, management), fluctuating environmental conditions, or changes in human or animal behavior (e.g. encroachment may bring reservoir hosts in closer proximity) [3, 6]. Long-term surveillance data can also reveal shifts in the temporal dynamics of tick populations and communities . Although the seasonality of known tick species has been described, year-to-year distribution of tick species may be influenced by inter-annual variability (e.g. local climate) and biotic factors (e.g. local reservoir species abundance). These data can be useful for developing predictive models that accurately measure the risks of tick-borne zoonotic agents.
We present a retrospective analysis of tick collection data in Pennsylvania from the early 1900s to June of 2017. Some of the data prior to 1968 had been published in a progress report on ticks from Pennsylvania but were presented in a format that included anecdotes and overall percentages rather than raw number breakdowns by species . We revisit these specimens and utilize the raw data from both these time periods (1900–1967 and 1968–June 2017) to identify shifts in tick community composition, phenological patterns, and host associations. We used our database to map the distribution of major tick species at the county level, investigate tick community spatiotemporal dynamics, and explore host associations by tick species.
The state of Pennsylvania (PA) is located in the Northeast/Mid-Atlantic region of the USA (State Center: 40.9699889, −77.7278831 . Climate in Pennsylvania is variable by location, but broadly classified as a continental type with warm, humid summers (mean temperature ranges of 17.7–23.33°C, . The majority of Pennsylvania’s land-use is dedicated to agriculture (both croplands and pastures), forestland, with some dense urban areas [9, 10]. There have been significant changes in the human populations of PA from 1960 to 2010, but a large proportion of the PA population has remained heavily clustered around Philadelphia and Pittsburgh, which are located in south-eastern and in south-western PA respectively (Additional file 1: Figure S1).
The PSU Frost Entomological Museum (‘Frost Museum’) houses arthropod samples collected by researchers, teaching collections, and samples submitted by the public for identification. We present our analysis of the tick specimens from 1900 to June 2017, although some Frost Museum collections date as far back as the late 1800s. Because tick samples were submitted to the Department of Entomology or the Frost Museum over a period of 117 years, they represent multiple collection/submission periods (early 1900–1959; 1960–1969; 1970–1988; Tick Research Lab (TRL) submissions from 1990 to 1993 and 1995–present). Two public campaigns account for the majority of the citizen-submitted specimens. The first campaign (between 1963–1967) was conducted by Dr Robert Snetsinger. He enlisted the help of the public through advertisements in radio, television, and newspapers to obtain 700 specimens . Additionally, he utilized active surveillance methods to collect approximately 500 ticks using a combination of dragging, sweeping, live animal trapping, and roadkill examinations of mammals and birds to assess tick abundance in localized areas . A second funded campaign dedicated to estimating tick abundance and species diversity was launched by Steven B. Jacobs (second author) from 1990 to 1993 (TRL), in which he cataloged, identified, and labeled each specimen. Specimens were accompanied by additional data: date of tick discovery, location and vegetation type of tick encounter, and host species.
Data from both campaigns and subsequent submissions were combined into a single dataset for our analyses. For analysis on the distribution of tick species over time, we used total tick counts. For quantifying host association, however, we used “submission” number which we defined as a vial or lot containing one or more ticks. We chose to use this more conservative measure rather than total specimen count to avoid skews in abundance by hosts. For example, a submission lot of 1 tick versus 50 ticks from a host were both classified as “one submission”.
While most specimens were collected within state boundaries, a few were declared from people either visiting or returning from visiting other states. Tick specimens identified as species that are not commonly found in Pennsylvania were later discovered to have been imported from other states/countries or found on exotic animals. Non-PA data were excluded from state-wide analyses, but were included in Additional file 2: Tables S1 and Additional file 3: Table S2).
Ticks were morphologically identified to species and life stage using taxonomic keys for Argasidae, Ixodidae east of the Mississippi, Dermacentor, nymphal Ixodes, and nymphs of Amblyomma [11–17]. Species-level identification is crucial since at least 3 Dermacentor species, 3 species of Amblyomma, and 9 different Ixodes species have been reported in Pennsylvania. If diagnostic characters were missing due to damage to the specimen, the next level of taxonomic identification was used (e.g. samples with missing mouthparts that were clearly Prostriata were identified as “Ixodes spp.”). In a few cases, samples were not identified beyond “tick” and were designated “Ixodidae” for hard ticks or “Argasidae” for soft ticks. Unusual specimens or those that were difficult to identify were sent to the National Tick Collection, Georgia Southern University for confirmation (by Dr James Oliver at the time of confirmation).
We focused on the geographical distribution of the five most abundant species of significant public health and veterinary importance: Amblyomma americanum (Linnaeus), Dermacentor variabilis (Say), Ixodes cookei (Packard), Ixodes scapularis (Say), and Rhipicephalus sanguineus (Latreille). Working on the assumption that counties with higher populations would submit more specimens than less populated counties, we estimated the prevalence rate (the total numbers of individual ticks per 100,000 people). This was done by adjusting the total tick count of each species by the countyʼs total population. We looked at relevant time periods during the surveillance programme: 1960–1969; 1990–1999; 2000–2009; and 2010–2018. For each time period, we used the United States Census data for 1960, 1990, 2000 and 2010, respectively, to calculate the tick prevalence rate (Additional file 1: Figure S1) [18–20]
We used the annual sum of all individual tick specimens to investigate how annual submission rates changed over time. We then analyzed the temporal dynamics of the five most abundant taxa (A. americanum, D. variabilis, I. cookei, I. scapularis and R. sanguineus). We did not evaluate total counts by year as these varied drastically due to the active campaigning for citizen submissions or the introduction of identification fees. Therefore, we looked at the proportional contribution of each species to the annual summed counts of all the five major species. To detect if there have been any monotonic trends (i.e. gradual shifts in abundance), we ran a non-parametric, two-sided Mann-Kendall trend test on the yearly proportion of each of the species between 1900–2017.
We investigated the seasonal distribution of the tick community by analyzing the monthly frequency of submissions. Citizen submissions were to include the date of discovery, but for specimens that lacked these, we used the date that a given submission was received. We analyzed the seasonal distribution of motile life stages (larvae, nymphs and adults) for the five most abundant taxa. The proportions were calculated by comparing the monthly abundance of each life stage (larvae, nymphs and adults) to the cumulative sum of all stages.
Host information (combined by family, except for dog, cat, human and groundhog) was available for many of the tick submissions. Host data were classified as either domestic or wildlife. We analyzed the host-tick data by summing the total submissions by both the tick species and the host group. We constructed a circular network map to visualize the relationships between tick species and hosts. All host association analyses were done with R (Version 3.4.13) with the packages Kendall for the Mann-Kendall test and the circlize package for chord diagrams of host association mapping [21, 22].
The total submissions to the PSU Department of Entomology/Frost Entomological Museum from 1900 to 2017. Generic names that have been changed since the submission date are shown in parentheses
No. of submission lots
No of specimens
Amblyomma (Aponomma) transversale
Carios (Ornithodoros) kelleyi
Dermacentor variabilis distribution was largely localized to the southern parts of the state. From 1990 to 2000, the highest proportion of D. variabilis submissions came from Greene County, the most southeastern county of Pennsylvania (865.45 submissions per 100,000). Other southern counties with high tick loads per capita included Fulton County (350.60 per 100,000) and Franklin County (117.26 per 100,000). Ixodes cookei was more evenly distributed throughout Pennsylvania, although like I. scapularis, it was more highly abundant in the northern counties. In 1990–2000, Forest County had the highest prevalence rates of Ixodes cookei with 80.87 per 100,000. R. sanguineus and A. americanum had very few submissions and their distribution was mostly scattered across Pennsylvania (Fig. 2).
For the tick species with less than 150 submissions across 1900 to 2017, we aggregated the submissions by genus. Multiple species within the genera Ixodes and Dermacentor were widely distributed across Pennsylvania (I. scapularis, I. cookei, D. andersoni and D. albipictus) (Additional file 4: Figure S2). Other species in the genera Amblyomma, Argas, Carios (Ornithodoros) and Haemaphysalis were not as widely distributed, possibly because these species are not commonly encountered or because the specimens were introduced from their native geographical ranges (Additional file 4: Figure S2, Additional file 5: Figure S3).
Temporal shifts in species abundance
Prior to 1990 I. scapularis were very rare so we have only shown their lifestage-specific seasonal abundance since 1990. Overall, the nymphal and larval submissions showed a unimodal pattern with the highest proportion of submissions received in June. For the adult submissions, there were prominent bimodal peaks in May and October with similar proportions of submissions received in both seasons.
The majority of I. cookei submissions were nymphs (n = 521), followed by adults (n = 182) and larvae (n = 88). The submission patterns indicate that I. cookei specimens can be encountered year-round, but that nymphs were the most commonly encountered lifestage. Across all lifestages, we see that the distributions are unimodal with peaks in early summer between May and June.
There were too few lifestage submissions for A. americanum and R. sanguineus to make sufficient comparisons with the seasonality of these species.
Our data is unique in that it contains details about tick community composition and spatio-temporal dynamics from Pennsylvania over a 117-year period. Subsets of our data had been reported as percentages or combined with data from other museums and literature reviews to estimate the distribution of one or more tick species across the state of Pennsylvania [7, 23]. To our knowledge, this is the first time that these data have been compiled in their entirety and analyzed in this format. We were able to detect seasonality, shifts in tick community composition, and host associations that have not been well-documented in a quantitative manner. The seasonality data for the five most abundant tick species inferred by our passive surveillance data is consistent with previous records of seasonality described by other researchers [24–27], demonstrating that these types of passive data contain biologically meaningful signal.
Shifts in tick community composition
In the 1960s, PA tick communities consisted predominantly of three species: I. cookei, D. variabilis, and R. sanguineus (Fig. 2). The most abundant species at that time, I. cookei, often referred to as a groundhog tick, is actually a broad- host tick feeding principally on medium mammals, although humans and dogs will also be parasitized . The second most abundant species, D. variabilis, was widely distributed and eventually became the dominant species submitted over I. cookei in the 1990s. By 1991 D. variabilis had been identified from all but 4 counties (Fig. 2) . After 1995, D. variabilis annual submission rates declined as I. scapularis submission rates increased (Fig. 2). Because there were gaps in submission numbers for certain years, we cannot say for certain why there were shifts in the abundance of these species.
Although we cannot directly infer a causal negative relationship between D. variabilis and I. cookei with I. scapularis, this pattern was also observed in neighboring Ohio. The passive surveillance programme run by Ohio Department of Health (started in 1978) did not detect I. scapularis (=formerly I. dammini) until 1989 . At that time, the dominant species were D. variabilis (~97% of submissions) and I. cookei (1.2%) . Between 1989 until 2008, I. scapularis accounted for less than 1% of the total submissions, but after 2009, the abundance began to increase. By 2012 they accounted for 24.8% of ticks submitted to the Ohio Department of Health .
Ixodes cookei abundance was highest prior to 1990 but has since become very rare in our dataset. Yet, we know that the abundance of I. cookei in the Maine passive surveillance programme has been constant, even as I. scapularis submissions have increased . It is possible that I. cookei abundance in PA has also remained relatively stable, but that we lack sufficient power to detect I. cookei. This inconsistency in submission rates may be due to the shift from free to per-submission charges for tick identifications that occurred in the mid-late 1990s.
The fourth and fifth most abundant tick species in our database were R. sanguineus, the brown dog tick (287 submission lots consisting of 332 specimens) and A. americanum (183 submission lots consisting of 196 specimens) (Table 1). Although R. sanguineus originated in Africa, it is a cosmopolitan urban pest species found worldwide in association with humans and their canine companions . Snetsinger  suggested in 1968 that R. sanguineus had established breeding populations in Pennsylvania, but the abundance tapered off after 1968 and none exist in our database since 2002. In contrast, the submission rates of A. americanum increased and then leveled off from 1960 to 2000. Although not commonly encountered, we have specimens from as recent as 2016.
Vector-host associations (including host specificity) are important for predicting the risk of pathogen transmission and identifying key players in a sylvatic disease cycle. The specificity to host varies with the tick species. In our dataset generalist tick species (e.g. D. variabilis, I. scapularis, I. cookei) were found parasitizing a wide range of vertebrate hosts, while specialist species (e.g. I. dentatus, I. marxi, I. muris) were associated with a single host or limited to host size (e.g. small mammals or birds). This was consistent with other data in the literature.
Less commonly encountered tick species can sometimes lead to incorrect assumptions about host preferences and subsequent risk of pathogen transmission. While some tick species may be presumed to hold strict host preferences, they may bite humans if given the opportunity. For instance, I. dentatus bit humans in cabins that had been inhabited by their squirrel hosts in Maine and Vermont, and I. cookei was found on humans in West Virginia [32, 33]. Ixodes texanus was only found from raccoons in our dataset, but this species is known to feed on several mammalian hosts [24, 34, 35]. Commonly held misconceptions about host associations of certain tick species (e.g. I. cookei as “groundhog ticks”) based on lack of encounter may result in ignoring a potentially epidemiologically important vector (e.g. I. cookei is a vector of Powassan encephalitis virus and may potentially be another vector of Ehrlichia muris ).
Multi-faceted approach to tick surveillance
Surveillance can be a powerful tool for the detection of introduced species (transient or established), emergent arthropod-borne pathogens, and disease risks due to increases or changes in vector community composition. Both passive and active surveillance strategies have their strengths and weaknesses but combined, they provide a more complete picture of tick community dynamics. Active tick surveillance approaches such as dragging, flagging, CO2-trapping, or live animal capture, can be very effective for assessing tick load by habitat [7, 24, 34, 37]. It can, however, be labor-intensive, costly, and difficult to implement over a wide geographical area. Passive surveillance is more cost-effective and less labor-intensive and can provide insight into ectoparasite abundance, host associations, or habitat associations across a wider geographical area . Passive surveillance (particularly based on submissions by citizens) may run the risk of under-representing certain taxa or reflect a bias toward certain host associations. However, citizen-submitted tick collections can provide valuable baseline data on prevalence and likelihood of tick encounters and may be more strongly correlated with reported human cases of tick-borne diseases than active surveillance alone [39–42]. A community engagement programme that actively recruits ticks submitted by citizens should be coupled with support for a rigorously curated database of tick submissions.
Utilizing complementary strategies can help fill in knowledge gaps about tick prevalence. In a study using a combination of retrospective literature review, data compilation of specimens from archival collections, and active collection (dry ice, dragging and flagging) in counties presumed to be free of A. americanum, 68 of 77 counties of Oklahoma were determined to be colonized . The metadata associated with a multi-pronged approach to tick surveillance (assuming proper data curation and management) can provide insight into tick-host associations, vegetation, seasonality, and shifts in population structure that can be used for modeling disease risk. Archival tick samples (or their DNA) can be useful for retrospective mining for research on the population genetics of ticks to detect gene flow, host shifts, or on their microbial inhabitants.
Implicit in any tick surveillance strategy is having trained tick biologists who can readily distinguish species by morphological and molecular characteristics. In the last 20 years, I. scapularis has become the most abundant tick species in Pennsylvania. While distinguishing Ixodes from other genera of ticks is fairly simple, species-level identification requires more detailed morphological examination, since there are six endemic species of Ixodes and three exotic species that could potentially be misidentified as I. scapularis. More generally, although many tick species are incompetent vectors of B. burgdorferi, they may be vectors and/or reservoirs of other pathogens/parasites, or acquire pathogens during co-feeding [24, 43–46]. It is therefore important to correctly identify tick species, not only for the determination of disease risk but also because the treatments for the pathogens they transmit may differ significantly.
An ideal tick surveillance programme would not only utilize multiple approaches and have a dedicated tick biologist proficient at species identification on staff, but it would also take a proactive stance that is not limited strictly to immediate threats. Since 1993 (~25 years) there have been 28 publications on ticks from Pennsylvania, and 22 of them were focused on I. scapularis and/or the microbiota (mostly on pathogens) [23, 25, 46–71]. While I. scapularis is an important vector that warrants this attention, other tick species are being ignored. The consequences of neglecting other potential ticks of epidemiological significance include missing shifts in tick biodiversity, not identifying the potential causes of said shifts, not monitoring changes in range expansions of vectors, and not detecting the presence of introduced or established species. Recently, the presence of H. longicornis was reported in Pennsylvania, but we have no data on whether it has been introduced previously, or whether it has established populations. Given the potential risks that this parthenogenic tick poses (wide host range, possible vector of multiple pathogens, may induce meat allergies), this example highlights the potential dangers of focusing on only one vector-pathogen system [70–73]. Hybrid datasets from tick collections derived from multiple sources represent a powerful tool for mining past ecological and epidemiological events. Many states maintain county records on passive tick submissions to veterinary or medical health officials, but there may be other cryptic collections (and associated data) housed in museums, universities, government institutions, or with private individuals. Combining these data with other ectoparasite databases and currently unexplored collections will provide ectoparasite researchers the robust dataset needed for a massive meta-analysis. These cryptic ectoparasite collections will provide the basis for exploring hypotheses such as: (i) are shifts in tick populations correlated with increasing human encroachment on natural habitats; (ii) what are some phenological reasons for the increase in I. scapularis abundance; or (iii) if displacement of a dominant tick community species occurs, what are the implications for tick-borne disease risk? We anticipate that participating in such a study will fill in the gaps of knowledge about less-studied tick species as well as highlight the intrinsic value of museum collections of ectoparasites.
The authors would like to thank Dr István Mikó for critical manuscript review, citizens and researchers who have contributed ticks to PSU, The Frost Entomological Museum for access to the collection, and critical suggestions from reviewers.
This work was supported by the USDA National Institute of Food and Agriculture and Hatch Appropriations under Project #PEN04691 and Accession #1018545. Additional support came from NSF GRFP DGE1255832 (to DP), and start-up funds from the Huck Institutes of Life Sciences, and the Penn State College of Agriculture to JMS.
Availability of data and materials
Data supporting the conclusions of this article are provided within the article and its additional files. The datasets supporting the conclusions of this article are available at GitHub in the passive_tick_surveillance_2018repository, https://github.com/pakdamie/passive_surveillance_tick_2018/tree/master/MAIN_DAT, and at Zenodo (https://zenodo.org/record/1476091).
SBJ compiled that data into a digital database, collected metadata, identified tick specimens from 1989 until 2017 and approved manuscript revisions. DP generated data visualizations, conducted statistical analysis, and contributed to manuscript writing and revisions. JMS curated the database, directed analyses, and contributed literature review, manuscript writing, and revisions. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Not applicable as these specimens were in a museum collection without any identifying information beyond host, date of collection, and collection location.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Rosenberg R, Lindsey NP, Fischer M, Gregory CJ, Hinckley AF, Mead PS, et al. Vital Signs: Trends in reported vectorborne disease cases—United States and Territories, 2004–2016. Morb Mortal Wkly Rep. 2018;67:496–501.View ArticleGoogle Scholar
- Sonenshine DE. Range expansion of tick disease vectors in North America: implications for spread of tick-borne disease. Int J Environ Res Public Health. 2018;15(3):478.PubMed CentralView ArticleGoogle Scholar
- Simon JA, Marrotte RR, Desrosiers N, Fiset J, Gaitan J, Gonzalez A, et 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Eisen RJ, Kugeler KJ, Eisen L, Beard CB, Paddock CD. Tick-borne zoonoses in the United States: persistent and emerging threats to human health. ILAR J. 2017;58:319–35.PubMedPubMed CentralView ArticleGoogle Scholar
- Rand PW, Lacombe EH, Dearborn R, Cahill B, Elias S, Lubelczyk CB, et al. Passive surveillance in Maine, an area emergent for tick-borne diseases. J Med Entomol. 2007;44:1118–29.PubMedView ArticleGoogle Scholar
- Brownstein JS, Skelly DK, Holford TR, Fish D. Forest fragmentation predicts local scale heterogeneity of Lyme disease risk. Oecologia. 2005;146:469–75.PubMedView ArticleGoogle Scholar
- Snetsinger RJ. Distribution of ticks and tick-borne diseases in Pennsylvania. University Park, PA: Pennsylvania State University; 1968 June p. 1–8 Report No.: 288.Google Scholar
- Node: Pennsylvania (316987717). OpenStreetMap. https://www.openstreetmap.org/node/316987717
- Pennsylvania State Climatologist. http://climate.met.psu.edu/data/state/
- Land Use | NRCS Pennsylvania. https://www.nrcs.usda.gov/wps/portal/nrcs/main/pa/technical/landuse/
- Keirans JE, Litwak TR. Pictorial key to the adults of hard ticks, family Ixodidae (Ixodida: Ixodoidea), east of the Mississippi River. J Med Entomol. 1989;26:435–48.PubMedView ArticleGoogle Scholar
- Cooley RA, Kohls GM. The genus Amblyomma (Ixodidae) in the United States. J Parasitol. 1944;30:77.View ArticleGoogle Scholar
- Keirans JE, Clifford CM. The genus Ixodes in the United States: a scanning electron microscope study and key to the adults. J Med Entomol. 1978;15:1–38.View ArticleGoogle Scholar
- Yunker CE, Keirans JE, Clifford CM, Easton EA. Dermacentor ticks (Acari: Ixodoidea: Ixodidae) of the New World: a scanning electron microscope atlas. Proc Entomol Soc Wash. 1986;88:609–27.Google Scholar
- Durden LA, Keirans JE. Nymphs of the genus Ixodes of the United States. Lanham, Maryland, USA: Entomological Society of America; 1996.Google Scholar
- Keirans JE, Durden LA. Illustrated key to nymphs of the tick genus Amblyomma (Acari: Ixodidae) found in the United States. J Med Entomol. 1998;35:489–95.PubMedView ArticleGoogle Scholar
- Cooley RA, Kohls GM. The Argasidae of North America, Central America and Cuba. Notre Dame, USA: The University Press; 1944. http://www.biodiversitylibrary.org/bibliography/4511
- U.S. Census data from 1960 to 1990. https://www.census.gov/population/cencounts/pa190090.txt
- U.S. Census data: U.S. Census Bureau; Table DP-1, Census 2000, Profile of General Demographic Characteristics Summary File 1, https://factfinder.census.gov/bkmk/table/1.0/en/DEC/10_DP/DPDP1/0400000US42
- US Census data: U.S. Census Bureau; Table DP-1, Census 2010, Profile of General Demographic Characteristics Summary File 1, https://factfinder.census.gov/bkmk/table/1.0/en/DEC/00_SF1/DP1/0400000US42
- McLeod AI. Kendall: Kendall rank correlation and Mann-Kendall trend test. 2011. https://CRAN.R-project.org/package=Kendall
- Gu Z, Gu L, Eils R, Schlesner M, Brors B. circlize implements and enhances circular visualization in R. Bioinforma Oxf Engl. 2014;30:2811–2.View ArticleGoogle Scholar
- Snetsinger R, Jacobs SB, Kim KC, Tavris D. Extension of the range of Dermacentor variabilis (Acari: Ixodidae) in Pennsylvania. J Med Entomol. 1993;30:795–8.PubMedView ArticleGoogle Scholar
- Kollars TM, Oliver JH. Host associations and seasonal occurrence of Haemaphysalis leporispalustris, Ixodes brunneus, I. cookei, I. dentatus, and I. texanus (Acari: Ixodidae) in southeastern Missouri. J Med Entomol. 2003;40:103–7.PubMedView ArticleGoogle Scholar
- Simmons TW, Shea J, Myers-Claypole MA, Kruise R, Hutchinson ML. Seasonal activity, density, and collection efficiency of the blacklegged tick (Ixodes scapularis) (Acari: Ixodidae) in mid-western Pennsylvania. J Med Entomol. 2015;52:1260–9.PubMedView ArticleGoogle Scholar
- Burg JG. Seasonal activity and spatial distribution of host-seeking adults of the tick Dermacentor variabilis. Med Vet Entomol. 2001;15:413–21.PubMedView ArticleGoogle Scholar
- Kollars TM, Oliver JH, Durden LA, Kollars PG. Host associations and seasonal activity of Amblyomma americanum (Acari: Ixodidae) in Missouri. J Parasitol. 2000;86:1156–9.PubMedView ArticleGoogle Scholar
- Bishopp FC, Trembley HL. Distribution and hosts of certain North American ticks. J Parasitol. 1945;31:1–54.View ArticleGoogle Scholar
- Pretzman C, Daugherty N, Poetter K, Ralph D. The distribution and dynamics of Rickettsia in the tick population of Ohio. Ann N Y Acad Sci. 1990;590:227–36.PubMedView ArticleGoogle Scholar
- Wang P, Glowacki MN, Hoet AE, Needham GR, Smith KA, Gary RE, et al. Emergence of Ixodes scapularis and Borrelia burgdorferi, the Lyme disease vector and agent, in Ohio. Front Cell Infect Microbiol. 2014;4:70.PubMedPubMed CentralView ArticleGoogle Scholar
- Brites-Neto J, Duarte KMR, Martins TF. Tick-borne infections in human and animal population worldwide. Vet World. 2015;8:301–15.PubMedPubMed CentralView ArticleGoogle Scholar
- Hall JE, Amrine JW, Gais RD, Kolanko VP, Hagenbuch BE, Gerencser VF, et al. Parasitization of humans in West Virginia by Ixodes cookei (Acari: Ixodidae), a potential vector of Lyme borreliosis. J Med Entomol. 1991;28:186–9.PubMedView ArticleGoogle Scholar
- Lubelczyk C, Cahill BK, Hanson T, Turmel J, Lacombe E, Rand PW, et al. Tick (Acari: Ixodidae) infestation at two rural, seasonal camps in Maine and Vermont. J Parasitol. 2010;96:442–3.PubMedView ArticleGoogle Scholar
- Brillhart DB, Fox LB, Upton SJ. Ticks (Acari: Ixodidae) collected from small and medium-sized Kansas mammals. J Med Entomol. 1994;31:500–4.PubMedView ArticleGoogle Scholar
- Cohen SB, Freye JD, Dunlap BG, Dunn JR, Jones TF, Moncayo AC. Host associations of Dermacentor, Amblyomma, and Ixodes (Acari: Ixodidae) ticks in Tennessee. J Med Entomol. 2010;47:415–20.PubMedGoogle Scholar
- Xu G, Pearson P, Rich SM. Ehrlichia muris in Ixodes cookei ticks, northeastern United States, 2016–2017. Emerg Infect Dis. 2018;24:1143–4.PubMedPubMed CentralView ArticleGoogle Scholar
- Bouchard C, Leighton PA, Beauchamp G, Nguon S, Trudel L, Milord F, et al. Harvested white-tailed deer as sentinel hosts for early establishing Ixodes scapularis populations and risk from vector-borne zoonoses in southeastern Canada. J Med Entomol. 2013;50:384–93.PubMedView ArticleGoogle Scholar
- Oliver JD, Bennett SW, Beati L, Bartholomay LC. Range expansion and increasing Borrelia burgdorferi infection of the tick Ixodes scapularis (Acari: Ixodidae) in Iowa, 1990–2013. J Med Entomol. 2017;54:1727–34.PubMedView ArticleGoogle Scholar
- Ripoche M, Gasmi S, Adam-Poupart A, Koffi JK, Lindsay LR, Ludwig A, et al. Passive tick surveillance provides an accurate early signal of emerging Lyme disease risk and human cases in southern Canada. J Med Entomol. 2018;55:1016–26.PubMedView ArticleGoogle Scholar
- Cortinas R, Spomer SM. Occurrence and county-level distribution of ticks (Acari: Ixodoidea) in Nebraska using passive surveillance. J Med Entomol. 2014;51:352–9.PubMedView ArticleGoogle Scholar
- Barrett AW, Noden BH, Gruntmeir JM, Holland T, Mitcham JR, Martin JE, et al. County scale distribution of Amblyomma americanum (Ixodida: Ixodidae) in Oklahoma: addressing local deficits in tick maps based on passive reporting. J Med Entomol. 2015;52:269–73.PubMedView ArticleGoogle Scholar
- Xu G, Mather TN, Hollingsworth CS, Rich SM. Passive surveillance of Ixodes scapularis (Say), their biting activity, and associated pathogens in Massachusetts. Vector Borne Zoonotic Dis. 2016;16:520–7.PubMedPubMed CentralView ArticleGoogle Scholar
- Baer-Lehman ML, Light T, Fuller NW, Barry-Landis KD, Kindlin CM, Stewart RL. Evidence for competition between Ixodes scapularis and Dermacentor albipictus feeding concurrently on white-tailed deer. Exp Appl Acarol. 2012;58:301–14.PubMedView ArticleGoogle Scholar
- Brown SM, Lehman PM, Kern RA, Henning JD. Detection of Borrelia burgdorferi and Anaplasma phagocytophilum in the black-legged tick, Ixodes scapularis, within southwestern Pennsylvania. J Vector Ecol J Soc Vector Ecol. 2015;40:180–3.View ArticleGoogle Scholar
- Campagnolo ER, Tewari D, Farone TS, Livengood JL, Mason KL. Evidence of Powassan/deer tick virus in adult black-legged ticks (Ixodes scapularis) recovered from hunter-harvested white-tailed deer (Odocoileus virginianus) in Pennsylvania: a public health perspective. Zoonoses Public Health. 2018;65:589–94.PubMedView ArticleGoogle Scholar
- Courtney JW, Dryden RL, Montgomery J, Schneider BS, Smith G, Massung RF. Molecular characterization of Anaplasma phagocytophilum and Borrelia burgdorferi in Ixodes scapularis ticks from Pennsylvania. J Clin Microbiol. 2003;41:1569–73.PubMedPubMed CentralView ArticleGoogle Scholar
- Anderson JF, Duray PH, Magnarelli LA. Borrelia burgdorferi and Ixodes dammini prevalent in the greater Philadelphia area. J Infect Dis. 1990;161:811–2.PubMedView ArticleGoogle Scholar
- Crowder CD, Carolan HE, Rounds MA, Honig V, Mothes B, Haag H, et al. Prevalence of Borrelia miyamotoi in Ixodes ticks in Europe and the United States. Emerg Infect Dis. 2014;20:1678–82.PubMedPubMed CentralView ArticleGoogle Scholar
- Daniels TJ, Fish D, Levine JF, Greco MA, Eaton AT, Padgett PJ, et al. Canine exposure to Borrelia burgdorferi and prevalence of Ixodes dammini (Acari: Ixodidae) on deer as a measure of Lyme disease risk in the northeastern United States. J Med Entomol. 1993;30:171–8.PubMedView ArticleGoogle Scholar
- Devevey G, Brisson D. The effect of spatial heterogeneity on the aggregation of ticks on white-footed mice. Parasitology. 2012;139:915–25.PubMedView ArticleGoogle Scholar
- Dick CW, Gannon MR, Little WE, Patrick MJ. Ectoparasite associations of bats from central Pennsylvania. J Med Entomol. 2003;40:813–9.PubMedView ArticleGoogle Scholar
- Edwards MJ, Barbalato LA, Makkapati A, Pham KD, Bugbee LM. Relatively low prevalence of Babesia microti and Anaplasma phagocytophilum in Ixodes scapularis ticks collected in the Lehigh Valley region of eastern Pennsylvania. Ticks Tick-Borne Dis. 2015;6:812–9.PubMedView ArticleGoogle Scholar
- Han GS, Stromdahl EY, Wong D, Weltman AC. Exposure to Borrelia burgdorferi and other tick-borne pathogens in Gettysburg National Military Park, South-Central Pennsylvania, 2009. Vector Borne Zoonotic Dis. 2014;14:227–33.PubMedView ArticleGoogle Scholar
- Hutchinson ML, Strohecker MD, Simmons TW, Kyle AD, Helwig MW. Prevalence rates of Borrelia burgdorferi (Spirochaetales: Spirochaetaceae), Anaplasma phagocytophilum (Rickettsiales: Anaplasmataceae), and Babesia microti (Piroplasmida: Babesiidae) in host-seeking Ixodes scapularis (Acari: Ixodidae) from Pennsylvania. J Med Entomol. 2015;52:693–8.PubMedView ArticleGoogle Scholar
- Lo Re V, Occi JL, MacGregor RR. Identifying the vector of Lyme disease. Am Fam Phys. 2004;69:1935–7.Google Scholar
- Lord RD, Lord VR, Humphreys JG, McLean RG. Distribution of Borrelia burgdorferi in host mice in Pennsylvania. J Clin Microbiol. 1994;32:2501–4.PubMedPubMed CentralGoogle Scholar
- Magnarelli LA, Stafford KC, Mather TN, Yeh MT, Horn KD, Dumler JS. Hemocytic rickettsia-like organisms in ticks: serologic reactivity with antisera to ehrlichiae and detection of DNA of agent of human granulocytic ehrlichiosis by PCR. J Clin Microbiol. 1995;33:2710–4.PubMedPubMed CentralGoogle Scholar
- Miller MJ, Esser HJ, Loaiza JR, Herre EA, Aguilar C, Quintero D, et al. Molecular ecological insights into Neotropical bird-tick interactions. PloS One. 2016;11:e0155989.PubMedPubMed CentralView ArticleGoogle Scholar
- Rogers MB, Cui L, Fitch A, Popov V, Travassos APA, Vasilakis N, et al. Whole genome analysis of Sierra Nevada virus, a novel mononegavirus in the family nyamiviridae. Am J Trop Med Hyg. 2014;91:159–64.PubMedPubMed CentralView ArticleGoogle Scholar
- Sakamoto JM, Ng TFF, Suzuki Y, Tsujimoto H, Deng X, Delwart E, et al. Bunyaviruses are common in male and female Ixodes scapularis ticks in central Pennsylvania. PeerJ. 2016;4:e2324.PubMedPubMed CentralView ArticleGoogle Scholar
- Sakamoto JM, Goddard J, Rasgon JL. Population and demographic structure of Ixodes scapularis Say in the eastern United States. PloS One. 2014;9:e101389.PubMedPubMed CentralView ArticleGoogle Scholar
- Schoelkopf L, Hutchinson CE, Bendele KG, Goff WL, Willette M, Rasmussen JM, et al. New ruminant hosts and wider geographic range identified for Babesia odocoilei (Emerson & Wright 1970). J Wildl Dis. 2005;41:683–90.PubMedView ArticleGoogle Scholar
- Serfass TL, Rymon LM, Brooks RP. Ectoparasites from river otters in Pennsylvania. J Wildl Dis. 1992;28:138–40.PubMedView ArticleGoogle Scholar
- Shock BC, Moncayo A, Cohen S, Mitchell EA, Williamson PC, Lopez G, et al. Diversity of piroplasms detected in blood-fed and questing ticks from several states in the United States. Ticks Tick Borne Dis. 2014;5:373–80.PubMedView ArticleGoogle Scholar
- Springer YP, Jarnevich CS, Barnett DT, Monaghan AJ, Eisen RJ. Modeling the present and future geographic distribution of the lone star tick, Amblyomma americanum (Ixodida: Ixodidae), in the continental United States. Am J Trop Med Hyg. 2015;93:875–90.PubMedPubMed CentralView ArticleGoogle Scholar
- Steiner FE, Pinger RR, Vann CN, Grindle N, Civitello D, Clay K, et al. Infection and co-infection rates of Anaplasma phagocytophilum variants, Babesia spp., Borrelia burgdorferi, and the rickettsial endosymbiont in Ixodes scapularis (Acari: Ixodidae) from sites in Indiana, Maine, Pennsylvania, and Wisconsin. J Med Entomol. 2008;45:289–97.PubMedView ArticleGoogle Scholar
- Stromdahl E, Hamer S, Jenkins S, Sloan L, Williamson P, Foster E, et al. Comparison of phenology and pathogen prevalence, including infection with the Ehrlichia muris-like (EML) agent, of Ixodes scapularis removed from soldiers in the midwestern and the northeastern United States over a 15-year period (1997–2012). Parasit Vectors. 2014;7:553.PubMedPubMed CentralView ArticleGoogle Scholar
- Waits K, Edwards MJ, Cobb IN, Fontenele RS, Varsani A. Identification of an anellovirus and genomoviruses in ixodid ticks. Virus Genes. 2018;54:155–9.PubMedView ArticleGoogle Scholar
- Yeh MT, Bak JM, Hu R, Nicholson MC, Kelly C, Mather TN. Determining the duration of Ixodes scapularis (Acari: Ixodidae) attachment to tick-bite victims. J Med Entomol. 1995;32:853–8.PubMedView ArticleGoogle Scholar
- Chen X, Xu S, Yu Z, Guo L, Yang S, Liu L, et al. Multiple lines of evidence on the genetic relatedness of the parthenogenetic and bisexual Haemaphysalis longicornis (Acari: Ixodidae). Infect Genet Evol J Mol Epidemiol Evol Genet Infect Dis. 2014;21:308–14.Google Scholar
- Chinuki Y, Ishiwata K, Yamaji K, Takahashi H, Morita E. Haemaphysalis longicornis tick bites are a possible cause of red meat allergy in Japan. Allergy. 2016;71:421–5.PubMedView ArticleGoogle Scholar
- Rainey T, Occi JL, Robbins RG, Egizi A. Discovery of Haemaphysalis longicornis (Ixodida: Ixodidae) parasitizing a sheep in New Jersey, United States. J Med Entomol. 2018;55:757–9.PubMedView ArticleGoogle Scholar
- Tufts DM, VanAcker MC, Fernandez MP, DeNicola A, Egizi A, Diuk-Wasser MA. Distribution, host-seeking phenology, and host and habitat associations of Haemaphysalis longicornis ticks, Staten Island, New York, USA. Emerg Infect Dis. 2019;25:792–6.PubMedPubMed CentralView ArticleGoogle Scholar