Population dynamics of the Lyme disease bacterium, Borrelia burgdorferi, during rapid range expansion in New York State.

Recent changes in climate and human land-use have resulted in alterations of the geographic range of many species, including human pathogens. Geographic range expansion and population growth of human pathogens increase human disease risk. Relatively little empirical work has investigated the impact of range changes on within-population variability, a contributor to both colonization success and adaptive potential, during the precise time in which populations are colonized. This is likely due to the difficulties of collecting appropriate natural samples during the dynamic phase of migration and colonization. We systematically collected blacklegged ticks (Ixodes scapularis) across New York State (NY), USA, between 2006 and 2019, a time period coinciding with a rapid range expansion of ticks and their associated pathogens including Borrelia burgdorferi, the etiological agent of Lyme disease. These samples provide a unique opportunity to investigate the population dynamics of human pathogens as they expand into novel territory. We observed that founder effects were short-lived, as gene flow from long-established populations brought almost all B. burgdorferi lineages to newly colonized populations within just a few years of colonization. By 7 years post-colonization, B. burgdorferi lineage frequency distributions were indistinguishable from long-established sites, indicating that local B. burgdorferi populations experience similar selective pressures despite geographic separation. The B. burgdorferi lineage dynamics elucidate the processes underlying the range expansion and demonstrate that migration into, and selection within, newly colonized sites operate on different time scales.

Assessing the impact of areal unit selection and the modifiable areal unit problem on associative statistics between cases of tick-borne disease and entomological indices.

The modifiable areal unit problem (MAUP) is a cause of statistical and visual bias when aggregating data according to spatial units, particularly when spatial units may be changed arbitrarily. The MAUP is a concern in vector-borne disease research when entomological metrics gathered from point-level sampling data are related to epidemiological data aggregated to administrative units like counties or ZIP Codes. Here, we assess the statistical impact of the MAUP when calculating correlations between randomly aggregated cases of anaplasmosis in New York State during 2017 and a geostatistical layer of an entomological risk index for Anaplasma phagocytophilum in blacklegged ticks (Ixodes scapularis Say, Acari: Ixodidae) collected during the fall of 2017. Correlations were also calculated using various administrative boundaries for comparison. We also demonstrate the impact of the MAUP on data visualization using choropleth maps and offer pycnophylactic interpolation as an alternative. Polygon simulations indicate that increasing the number of polygons decreases correlation coefficients and their variability. Correlation coefficients calculated using ZIP Code tabulation area and Census tract polygons were beyond 4 standard deviations from the mean of the simulated correlation coefficients. These results indicate that using smaller polygons may not best incorporate the geographical context of the tick-borne disease system, despite the tendency of researchers to strive for more granular spatial data and associations.

Spatiotemporal distribution of Borrelia miyamotoi (Spirochaetales: Spirochaetaceae) and coinfection with other tick-borne pathogens in host-seeking Ixodes scapularis (Acari: Ixodidae) from New York State, USA.

Blacklegged ticks (Ixodes scapularis Say, Acari: Ixodidae) were collected from 432 locations across New York State (NYS) during the summer and autumn of 2015-2020 to determine the prevalence and geographic distribution of Borrelia miyamotoi (Spirochaetales: Spirochaetaceae) and coinfections with other tick-borne pathogens. A total of 48,386 I. scapularis were individually analyzed using a multiplex real-time polymerase chain reaction assay to simultaneously detect the presence of Bo. miyamotoi, Borrelia burgdorferi (Spirochaetales: Spirochaetaceae), Anaplasma phagocytophilum (Rickettsiales: Anaplasmataceae), and Babesia microti (Piroplasmida: Babesiidae). Overall prevalence of Bo. miyamotoi in host-seeking nymphs and adults varied geographically and temporally at the regional level. The rate of polymicrobial infection in Bo. miyamotoi-infected ticks varied by developmental stage, with certain co-infections occurring more frequently than expected by chance. Entomological risk of exposure to Bo. miyamotoi-infected nymphal and adult ticks (entomological risk index [ERI]) across NYS regions in relation to human cases of Bo. miyamotoi disease identified during the study period demonstrated spatial and temporal variation. The relationship between select environmental factors and Bo. miyamotoi ERI was explored using generalized linear mixed effects models, resulting in different factors significantly impacting ERI for nymphs and adult ticks. These results can inform estimates of Bo. miyamotoi disease risk and further our understanding of Bo. miyamotoi ecological dynamics in regions where this pathogen is known to occur.

Associations of Anaplasma phagocytophilum Bacteria Variants in Ixodes scapularis Ticks and Humans, New York, USA.

Anaplasmosis, caused by the tickborne bacterium Anaplasma phagocytophilum, is an emerging public health threat in the United States. In the northeastern United States, the blacklegged tick (Ixodes scapularis) transmits the human pathogenic genetic variant of A. phagocytophilum (Ap-ha) and a nonpathogenic variant (Ap-V1). New York has recently experienced a rapid and geographically focused increase in cases of anaplasmosis. We analyzed A. phagocytophilum-infected I. scapularis ticks collected across New York during 2008-2020 to differentiate between variants and calculate an entomological risk index (ERI) for each. Ap-ha ERI varied between regions and increased in all regions during the final years of the study. Space-time scan analyses detected expanding clusters of Ap-ha located within documented anaplasmosis hotspots. Ap-ha ERI was more positively correlated with anaplasmosis incidence than non-genotyped A. phagocytophilum ERI. Our findings help elucidate the relationship between the spatial ecology of A. phagocytophilum variants and anaplasmosis.

First Detection of the Invasive Asian Longhorned Tick (Acari: Ixodidae) on Migratory Passerines in the Americas.

The Asian longhorned tick (Haemaphysalis longicornis Neumann), native to East Asia, was first reported in the United States in 2017 and is now established in at least 17 states. Haemaphysalis longicornis feeds on birds in its range outside of the United States, and migratory birds disperse this tick and tick-borne pathogens. However, early studies in the United States did not find H. longicornis on migrating passerine birds. The transport of the parthenogenetic H. longicornis on birds has the potential to greatly expand its range. We report the first discovery of H. longicornis on migratory passerine birds in the Americas.

Phylogeographic dynamics of the arthropod vector, the blacklegged tick (Ixodes scapularis).

BACKGROUND: The emergence of vector-borne pathogens in novel geographic areas is regulated by the migration of their arthropod vectors. Blacklegged ticks (Ixodes scapularis) and the pathogens they vector, including the causative agents of Lyme disease, babesiosis and anaplasmosis, continue to grow in their population sizes and to expand in geographic range. Migration of this vector over the previous decades has been implicated as the cause of the re-emergence of the most prevalent infectious diseases in North America. METHODS: We systematically collected ticks from across New York State (hereafter referred to as New York) from 2004 to 2017 as part of routine tick-borne pathogen surveillance in the state. This time frame corresponds with an increase in range and incidence of tick-borne diseases within New York. We randomly sampled ticks from this collection to explore the evolutionary history and population dynamics of I. scapularis. We sequenced the mitochondrial genomes of each tick to characterize their current and historical spatial genetic structure and population growth using phylogeographic methods. RESULTS: We sequenced whole mitochondrial genomes from 277 ticks collected across New York between 2004 and 2017. We found evidence of population genetic structure at a broad geographic scale due to differences in the relative abundance, but not the composition, of haplotypes among sampled ticks. Ticks were often most closely related to ticks from the same and nearby collection sites. The data indicate that both short- and long-range migration events shape the population dynamics of blacklegged ticks in New York. CONCLUSIONS: We detailed the population dynamics of the blacklegged tick (Ixodes scapularis) in New York during a time frame in which tick-borne diseases were increasing in range and incidence. Migration of ticks occurred at both coarse and fine scales in the recent past despite evidence of limits to gene flow. Past and current tick population dynamics have implications for further range expansion as habitat suitability for ticks changes due to global climate change. Analyses of mitochondrial genome sequencing data will expound upon previously identified drivers of tick presence and abundance as well as identify additional drivers. These data provide a foundation on which to generate testable hypotheses on the drivers of tick population dynamics occurring at finer scales.

Predicting spatio-temporal population patterns of Borrelia burgdorferi, the Lyme disease pathogen.

The causative bacterium of Lyme disease, Borrelia burgdorferi, expanded from an undetected human pathogen into the etiologic agent of the most common vector-borne disease in the United States over the last several decades. Systematic field collections of the tick vector reveal increases in the geographic range and prevalence of B. burgdorferi-infected ticks that coincided with increases in human Lyme disease incidence across New York State.We investigate the impact of environmental features on the population dynamics of B. burgdorferi. Analytical models developed using field collections of nearly 19,000 nymphal Ixodes scapularis and spatially and temporally explicit environmental features accurately explained the variation in the nymphal infection prevalence of B. burgdorferi across space and time.Importantly, the model identified environmental features reflecting landscape ecology, vertebrate hosts, climatic metrics, climate anomalies and surveillance efforts that can be used to predict the biogeographical patterns of B. burgdorferi-infected ticks into future years and in previously unsampled areas.Forecasting the distribution and prevalence of a pathogen at fine geographic scales offers a powerful strategy to mitigate a serious public health threat. Synthesis and applications. A decade of environmental and tick data was collected to create a model that accurately predicts the infection prevalence of Borrelia burgdorferi over space and time. This predictive model can be extrapolated to create a high-resolution risk map of the Lyme disease pathogen for future years that offers an inexpensive approach to improve both ecological management and public health strategies to mitigate disease risk.

Epidemiology and Spatial Emergence of Anaplasmosis, New York, USA, 2010‒2018.

Human granulocytic anaplasmosis, a tickborne disease caused by the bacterium Anaplasma phagocytophilum, was first identified during 1994 and is now an emerging public health threat in the United States. New York state (NYS) has experienced a recent increase in the incidence of anaplasmosis. We analyzed human case surveillance and tick surveillance data collected by the NYS Department of Health for spatiotemporal patterns of disease emergence. We describe the epidemiology and growing incidence of anaplasmosis cases reported during 2010-2018. Spatial analysis showed an expanding hot spot of anaplasmosis in the Capital Region, where incidence increased >8-fold. The prevalence of A. phagocytophilum increased greatly within tick populations in the Capital Region over the same period, and entomologic risk factors were correlated with disease incidence at a local level. These results indicate that anaplasmosis is rapidly emerging in a geographically focused area of NYS, likely driven by localized changes in exposure risk.

A Comparative Spatial and Climate Analysis of Human Granulocytic Anaplasmosis and Human Babesiosis in New York State (2013-2018).

Human granulocytic anaplasmosis (HGA) and human babesiosis are tick-borne diseases spread by the blacklegged tick (Ixodes scapularis Say, Acari: Ixodidae) and are the result of infection with Anaplasma phagocytophilum and Babesia microti, respectively. In New York State (NYS), incidence rates of these diseases increased concordantly until around 2013, when rates of HGA began to increase more rapidly than human babesiosis, and the spatial extent of the diseases diverged. Surveillance data of tick-borne pathogens (2007 to 2018) and reported human cases of HGA (n = 4,297) and human babesiosis (n = 2,986) (2013-2018) from the New York State Department of Health (NYSDOH) showed a positive association between the presence/temporal emergence of each pathogen and rates of disease in surrounding areas. Incidence rates of HGA were higher than human babesiosis among White and non-Hispanic/non-Latino individuals, as well as all age and sex groups. Human babesiosis exhibited higher rates among non-White individuals. Climate, weather, and landscape data were used to build a spatially weighted zero-inflated negative binomial (ZINB) model to examine and compare associations between the environment and rates of HGA and human babesiosis. HGA and human babesiosis ZINB models indicated similar associations with forest cover, forest land cover change, and winter minimum temperature; and differing associations with elevation, urban land cover change, and winter precipitation. These results indicate that tick-borne disease ecology varies between pathogens spread by I. scapularis.

Optimal Collection Methods for Asian Longhorned Ticks (Ixodida: Ixodidae) in the Northeast United States.

The Asian longhorned tick, Haemaphysalis longicornis Neumann, is an invasive species in the United States. Since its earliest recorded presence in West Virginia in 2010, H. longicornis has been reported from 15 states. While its public health significance in the United States is unclear, globally it transmits pathogens that infect livestock and humans, causing economic losses and substantial morbidity. Management and control of H. longicornis requires knowledge of its biology, ecology, and distribution. Here, we address the need for effective collection methods for host-seeking H. longicornis as an important step for accurately assessing tick abundance and potential disease risk. The number of H. longicornis collected were compared across three collection methods (dragging, sweeping, CO2 traps) and three tick check distances (5 m, 10 m, and 20 m) were compared for dragging and sweeping. Field collections were conducted from June through August 2019 in Westchester County, New York, and ticks were grouped by life stage to assess collection method efficiency. Results indicated that implementing shorter (5 m) tick check distance was ideal for adult and nymphal collections. The dragging method proved better than sweeping for adult collections; however, there was no significant difference between the methods for nymphal collections, at any tick check distance evaluated. CO2 traps attracted H. longicornis, but additional research is necessary to devise an effective tick retaining method before the traps can be implemented in the field. The results are presented to inform and support H. longicornis surveillance and control programs across the nation.

Spatio-temporal variation in environmental features predicts the distribution and abundance of Ixodes scapularis.

Many species have experienced dramatic changes in both geographic range and population sizes in recent history. Increases in the geographic range or population size of disease vectors have public health relevance as these increases often precipitate the emergence of infectious diseases in human populations. Accurately identifying environmental factors affecting the biogeographic patterns of vector species is a long-standing analytical challenge, stemming from a paucity of data capturing periods of rapid changes in vector demographics. We systematically investigated the occurrence and abundance of nymphal Ixodes scapularis ticks at 532 sampling locations throughout New York State (NY), USA, between 2008 and 2018, a time frame that encompasses the emergence of diseases vectored by these ticks. Analyses of these field-collected data demonstrated a range expansion into northern and western NY during the last decade. Nymphal abundances increased in newly colonised areas, while remaining stable in areas with long-standing populations over the last decade. These trends in the geographic range and abundance of nymphs correspond to both the geographic expansion of human Lyme disease cases and increases in incidence rates. Analytic models fitted to these data incorporating time, space, and environmental factors, accurately identified drivers of the observed changes in nymphal occurrence and abundance. These models accounted for the spatial and temporal variation in the occurrence and abundance of nymphs and can accurately predict nymphal population patterns in future years. Forecasting disease risk at fine spatial scales prior to the transmission season can influence both public health mitigation strategies and individual behaviours, potentially impacting tick-borne disease risk and subsequently human disease incidence.

Seasonal Activity of Haemaphysalis longicornis (Acari: Ixodidae) in Southern New York State.

The Asian longhorned tick, Haemaphysalis longicornis Neumann, is a species native to eastern Asia that has recently been discovered in the United States. In its native range, H. longicornis transmits pathogens that cause disease in humans and livestock. It is currently unknown whether H. longicornis will act as a vector in the United States. Understanding its seasonal activity patterns will be important in identifying which times of the year represent greatest potential risk to humans and livestock should this species become a threat to animal or public health. A study site was established in Yonkers, NY near the residence associated with the first reported human bite from H. longicornis in the United States. Ticks were collected once each week from July 2018 to November 2019. Haemaphysalis longicornis larvae were most active from August to November, nymphs from April to July, and adult females from June to September. This pattern of activity suggests that H. longicornis is capable of completing a generation within a single year and matches the patterns observed in its other ranges in the northern hemisphere. The data presented here contribute to a growing database for H. longicornis phenology in the northeastern United States. Potential implications of the short life cycle for the tick's vectorial capacity are discussed.

First glimpse into the origin and spread of the Asian longhorned tick, Haemaphysalis longicornis, in the United States.

Established populations of Asian longhorned ticks (ALT), Haemaphysalis longicornis, were first identified in the United States (US) in 2017 by sequencing the mitochondrial cytochrome c oxidase subunit I (cox1) 'barcoding' locus followed by morphological confirmation. Subsequent investigations detected ALT infestations in 12, mostly eastern, US states. To gain information on the origin and spread of US ALT, we (1) sequenced cox1 from ALT populations across 9 US states and (2) obtained cox1 sequences from potential source populations [China, Japan and Republic of Korea (ROK) as well as Australia, New Zealand and the Kingdom of Tonga (KOT)] both by sequencing and by downloading publicly available sequences in NCBI GenBank. Additionally, we conducted epidemiological investigations of properties near its initial detection locale in Hunterdon County, NJ, as well as a broader risk analysis for importation of ectoparasites into the area. In eastern Asian populations (China/Japan/ROK), we detected 35 cox1 haplotypes that neatly clustered into two clades with known bisexual versus parthenogenetic phenotypes. In Australia/New Zealand/KOT, we detected 10 cox1 haplotypes all falling within the parthenogenetic cluster. In the United States, we detected three differentially distributed cox1 haplotypes from the parthenogenetic cluster, supporting phenotypic evidence that US ALT are parthenogenetic. While none of the source populations examined had all three US cox1 haplotypes, a phylogeographic network analysis supports a northeast Asian source for the US populations. Within the United States, epidemiological investigations indicate ALT can be moved long distances by human transport of animals, such as horses and dogs, with smaller scale movements on wildlife. These results have relevant implications for efforts aimed at minimizing the spread of ALT in the United States and preventing additional exotic tick introductions.

Molecular Characterization of Haemaphysalis Species and a Molecular Genetic Key for the Identification of Haemaphysalis of North America.

Haemaphysalis longicornis (Acari: Ixodidae), the Asian longhorned tick, is native to East Asia, but has become established in Australia and New Zealand, and more recently in the United States. In North America, there are other native Haemaphysalis species that share similar morphological characteristics and can be difficult to identify if the specimen is damaged. The goal of this study was to develop a cost-effective and rapid molecular diagnostic assay to differentiate between exotic and native Haemaphysalis species to aid in ongoing surveillance of H. longicornis within the United States and help prevent misidentification. We demonstrated that restriction fragment length polymorphisms (RFLPs) targeting the 16S ribosomal RNA and the cytochrome c oxidase subunit I (COI) can be used to differentiate H. longicornis from the other Haemaphysalis species found in North America. Furthermore, we show that this RFLP assay can be applied to Haemaphysalis species endemic to other regions of the world for the rapid identification of damaged specimens. The work presented in this study can serve as the foundation for region specific PCR-RFLP keys for Haemaphysalis and other tick species and can be further applied to other morphometrically challenging taxa.

First Recognized Human Bite in the United States by the Asian Longhorned Tick, Haemaphysalis longicornis.

We present the case summary of the first human recognized to have been bitten by the Haemaphysalis longicornis tick in the United States, which occurred in New York State. Subsequent field studies confirmed that this tick was present in multiple geographic locations near the patient's residence, including on manicured lawns.

Environmental Determinants of Aedes albopictus Abundance at a Northern Limit of Its Range in the United States.

Aedes albopictus is a vector of arboviruses with high rates of morbidity and mortality. The northern limit of Ae. albopictus in the northeastern United States runs through New York state (NYS) and Connecticut. We present a landscape-level analysis of mosquito abundance measured by daily counts of Ae. albopictus from 338 trap sites in 12 counties during May-September 2017. During the study period, the mean number of Ae. albopictus caught per day of trapping across all sites was 3.21. We constructed four sets of negative binomial generalized linear models to evaluate how trapping methodology, land cover, as well as temperature and precipitation at multiple time intervals influenced Ae. albopictus abundance. Biogents-Sentinel (BGS) traps were 2.78 times as efficient as gravid traps and 1.49 times as efficient as CO2-baited CDC light traps. Greater proportions of low- and medium-intensity development and low proportions of deciduous cover around the trap site were positively associated with increased abundance, as were minimum winter temperature and March precipitation. The cumulative precipitation within a 28-day time window before the date of collection had a nonlinear relationship with abundance, such that greater cumulative precipitation was associated with increased abundance until approximately 70 mm, above which there was a decrease in abundance. We concluded that populations are established in Nassau, Suffolk, and New York City counties in NYS; north of these counties, the species is undergoing population invasion and establishment. We recommend that mosquito surveillance programs monitoring the northward invasion of Ae. albopictus place BGS traps at sites chosen with respect to land cover.

Multistate Infestation with the Exotic Disease-Vector Tick Haemaphysalis longicornis - United States, August 2017-September 2018.

Haemaphysalis longicornis is a tick indigenous to eastern Asia and an important vector of human and animal disease agents, resulting in such outcomes as human hemorrhagic fever and reduction of production in dairy cattle by 25%. H. longicornis was discovered on a sheep in New Jersey in August 2017 (1). This was the first detection in the United States outside of quarantine. In the spring of 2018, the tick was again detected at the index site, and later, in other counties in New Jersey, in seven other states in the eastern United States, and in Arkansas. The hosts included six species of domestic animals, six species of wildlife, and humans. To forestall adverse consequences in humans, pets, livestock, and wildlife, several critical actions are indicated, including expanded surveillance to determine the evolving distribution of H. longicornis, detection of pathogens that H. longicornis currently harbors, determination of the capacity of H. longicornis to serve as a vector for a range of potential pathogens, and evaluation of effective agents and methods for the control of H. longicornis.

Assessment of Duration of Tick Feeding by the Scutal Index Reduces Need for Antibiotic Prophylaxis After Ixodes scapularis Tick Bites.

Doxycycline is recommended for persons with Ixodes scapularis tick bites in certain geographic areas, if the tick had fed for at least 36 hours. Based on the scutal index, over 40% of I. scapularis tick bites from patients seen at the Lyme Disease Diagnostic Center did not warrant antibiotic prophylaxis.

Transient influence of blood meal and natural environment on blacklegged tick bacterial communities.

Blacklegged ticks (Ixodes scapularis) spend the majority of their life cycle off host, typically in woodland habitat, but require a blood meal at each of three life stages (larva, nymph, adult) to reach maturity and reproduce. Blood feeding usually lasts for several days each time and as blood is imbibed, a range of known pathogens from the host may also be acquired. Using next generation sequencing of 16S rRNA gene amplicons, we examined the influence of host blood meal on the internal bacterial community within nymphal blacklegged ticks across host-seeking, feeding, blood meal digestion, and after molting into the adult stage. Results demonstrate bacterial community structuring across host and ticks with 287 taxa found exclusively in ticks, suggesting the field environment plays a significant role in shaping the internal tick microbiome. A decrease in bacterial diversity was noted from unfed nymphs through feeding/digestion and after molting into adults, suggesting that bacterial species are lost during the corresponding physiological changes. The similarity in biochemical pathways across the different tick categories suggests that the loss of bacterial taxa does not mirror a large change in microbial function. Ticks likely lose bacterial taxa after feeding, but continual exposure to bacteria from the field environment counters this loss.

Energy Usage of Known-Age Blacklegged Ticks (Acari: Ixodidae): What Is the Best Method for Determining Physiological Age?

Ticks expend energy while host-seeking and must consume blood to advance to the next life stage. The energy required for activity is derived from the tick's lipid reserves, a valuable resource that sustains the tick until it finds the next host and can take another bloodmeal. The amount of lipid reserves in an unfed tick has been proposed as an index of tick biological age. Two different methods for aging nymphal blacklegged ticks, Ixodes scapularis Say, were analyzed in this study. To study lipid usage, colony-raised nymphs were held in lab-controlled chambers at 21.0 °C and ≥95% relative humidity, with a photoperiod of 12:12 (L:D) h. Samples of ticks were frozen at -80 °C every 2-3 wk, starting at 12-wk postmolt and continuing until 38-wk postmolt. Lipid reserves were determined indirectly through measurements of "physiological age" that estimate the energy a tick has based on the evaluation of morphometric size ratios of the tick scutum and alloscutum, and quantified directly through chloroform extractions of lipid from individual ticks. Morphometric age ratios and lipid amounts were compared to determine if morphometric measurements accurately estimated a tick's physiological state. Although the morphometric age ratio did correlate significantly with total tick lipid content, the predictive value of the ratio was not reliable; chloroform extraction results showed that lipid amounts declined steadily through the study and more accurately characterized the physiological condition of nymphal I. scapularis. The study of physiological aging of blacklegged ticks may lead to a better understanding of how changing environmental conditions affect tick longevity.