Predicting distributions of blacklegged ticks (Ixodes scapularis), Lyme disease spirochetes (Borrelia burgdorferi sensu stricto) and human Lyme disease cases in the eastern United States

Lyme disease is the most commonly reported vector-borne disease in the United States (US), with approximately 300,000 -to- 40,000 cases reported annually. The blacklegged tick, Ixodes scapularis, is the primary vector of the Lyme disease-causing spirochete, Borrelia burgdorferi sensu stricto, in high incidence regions in the upper midwestern and northeastern US. Using county-level records of the presence of I. scapularis or presence of B. burgdorferi s.s. infected host-seeking I. scapularis, we generated habitat suitability consensus maps based on an ensemble of statistical models for both acarological risk metrics. Overall accuracy of these suitability models was high (AUC = 0.76 for I. scapularis and 0.86 for B. burgdorferi s.s. infected-I. scapularis). We sought to compare which acarological risk metric best described the distribution of counties reporting high Lyme disease incidence (≥10 confirmed cases/100,000 population) by setting the models to a fixed omission rate (10%). We compared the percent of high incidence counties correctly classified by the two models. The I. scapularis consensus map correctly classified 53% of high and low incidence counties, while the B. burgdorferi s.s. infected-I. scapularis consensus map classified 83% correctly. Counties classified as suitable by the B. burgdorferi s.s. map showed a 91% overlap with high Lyme disease incidence counties with over a 38-fold difference in Lyme disease incidence between high- and low-suitability counties. A total of 288 counties were classified as highly suitable for B. burgdorferi s.s., but lacked records of infected-I. scapularis and were not classified as high incidence. These counties were considered to represent a leading edge for B. burgdorferi s.s. infection in ticks and humans. They clustered in Illinois, Indiana, Michigan, and Ohio. This information can aid in targeting tick surveillance and prevention education efforts in counties where Lyme disease risk may increase in the future.     

Anaplasma phagocytophilum, Babesia microti, and Borrelia burgdorferi in Ixodes scapularis, southern coastal Maine

Ixodes scapularis (deer ticks) from Maine were tested for multiple infections by polymerase chain reaction and immunofluorescence. In 1995, 29.5%, 9.5%, and 1.9% of deer ticks were infected with Borrelia burgdorferi, Anaplasma phagocytophilum, and Babesia microti, respectively. In 1996 and 1997, the number of A. phagocytophilum-infected ticks markedly declined. In 1995 through 1996, 4 (1.3%) of 301 were co-infected.     

Borrelia burgdorferi infection and cutaneous Lyme disease, Mexico

Four patients who had received tick bites while visiting forests in Mexico had skin lesions that met the case definition of erythema migrans, or borrelial lymphocytoma. Clinical diagnosis was supported with histologic, serologic, and molecular tests. This study suggests the Borrelia burgdorferi infection is in Mexico.     

Comparing the relative potential of rodents as reservoirs of the Lyme disease spirochete (Borrelia burgdorferi)

The authors compared the contribution of white-footed mice (Peromyscus leucopus), chipmunks (Tamias striatus), and meadow voles (Microtus pennsylvanicus) to infection of vector ticks with the Lyme disease spirochete, Borrelia burgddorferi. At one Massachusetts location where Lyme disease is endemic, all three species of rodents were found to be infected. Prevalence of infection, however, varied from 90% for mice, and 75% for chipmunks to just 5.5% for meadow voles. Infectivity of these hosts for larval Ixodes dammini also varied, but mice were found to be the most infective, followed by chipmunks and meadow voles. Density estimates of these three hosts, collected between 1981 and 1986 in three coastal Massachusetts locations, revealed that mice were more abundant than the other two rodents in areas in which ticks were abundant. In addition, mice were infested more abundantly with larval I. dammini than the other two host species. Integrating these results, we determined each species' "reservoir potential," a novel term which describes the relative contribution made by a host species to the horizontal infection of a vector population. The authors' findings demonstrate that, at least in coastal Massachusetts, P. leucopus is the most important small mammal reservoir for B. burgdorferi.     

Failure of the Asian longhorned tick,Haemaphysalis longicornis,to serve as an experimental vector of the Lyme disease spirochete,Borrelia burgdorferi sensu stricto

The invasive, human-biting Asian longhorned tick, Haemaphysalis longicornis, was detected in New Jersey in the eastern United States in August of 2017 and by November of 2018 this tick had been recorded from 45 counties across 9 states, primarily along the Eastern Seaboard. The establishment of H. longicornis in the United States has raised the questions of how commonly it will bite humans and which native pathogens may naturally infect this tick. There also is a need for experimental vector competence studies with native pathogens to determine if H. longicornis can acquire a given pathogen while feeding, pass it transstadially, and then transmit the pathogen in the next life stage. In this experimental study, we evaluated the vector competence of a population of H. longicornis originating from the United States (New York) for a native isolate (B31) of the Lyme disease spirochete, Borrelia burgdorferi sensu stricto (s.s.). In agreement with a previous experimental study on the vector competence of H. longicornis for Borrelia garinii, we found that uninfected H. longicornis larvae could acquire B. burgdorferi s.s. while feeding on infected Mus musculus mice (infection prevalence >50% in freshly fed larvae) but that the infection was lost during the molt to the nymphal stage. None of 520 tested molted nymphs were found to be infected, indicating that transstadial passage of B. burgdorferi s.s. is absent or rare in H. longicornis; and based on the potential error associated with the number of nymphs testing negative in this study, we estimate that the upper 95% limit for infection prevalence was 0.73%. An Ixodes scapularis process control showed both effective acquisition of B. burgdorferi s.s. from infected mice by uninfected larvae and transstadial passage to the nymphal stage (infection prevalence of 80–82% for both freshly fed larvae and molted nymphs). We also observed that although H. longicornis larvae could be compelled to feed on mice by placing the ticks within feeding capsules, attachment and feeding success was minimal (<0.5%) when larvae were placed freely on the fur of the mice. We conclude that H. longicornis is unlikely to contribute more than minimally, if at all, to transmission of Lyme disease spirochetes in the United States.     

Population genetic structure of the Lyme disease vector Ixodes scapularis at an apparent spatial expansion front

Modeling and empirical evidence suggests that Lyme disease is undergoing geographic expansion from principal foci in the midwestern and northeastern United States. Virginia is at the southern edge of the current expansion zone and has seen dramatic rise in human Lyme disease cases since 2007, potentially owing to a recent increase in vector abundance. Ixodes scapularis is known throughout the eastern US but behavioral or physiological variation between northern and southern lineages might lead northern-variant ticks to more frequently parasitize humans. We hypothesized that recent spatial and numerical increase in Lyme disease cases is associated with demographic and/or spatial expansion of I. scapularis and that signals of these phenomena would be detectable and discernable in population genetic signals. In summer and fall 2011, we collected nymphal I. scapularis by drag sampling and adult I. scapularis from deer carcasses at hunting check stations at nine sites arranged along an east–west transect through central Virginia. We analyzed 16S mtDNA sequences data from up to 24 I. scapularis individuals collected from each site and detected a total of 24 haplotypes containing 29 segregating sites. We found no evidence for population genetic structure among these sites but we did find strong signals of both demographic and spatial expansion throughout our study system. We found two haplotypes (one individual each) representing a lineage of ticks that is only found in the southeastern United States, with the remaining individuals representing a less genetically diverse clade that is typical of the northern United States, but that has also been detected in the American South. Taken together, these results lead us to conclude that I. scapularis populations in Virginia are expanding and that this expansion may account for recent observed increases in Lyme disease.     

A climate-based model predicts the spatial distribution of the Lyme disease vector Ixodes scapularis in the United States

An understanding of the spatial distribution of the black-legged tick, Ixodes scapularis, is a fundamental component in assessing human risk for Lyme disease in much of the United States. Although a county-level vector distribution map exists for the United States, its accuracy is limited by arbitrary categories of its reported presence. It is unknown whether reported positive areas can support established populations and whether negative areas are suitable for established populations. The steadily increasing range of I. scapularis in the United States suggests that all suitable habitats are not currently occupied. Therefore, we developed a spatially predictive logistic model for I. scapularis in the 48 conterminous states to improve the previous vector distribution map. We used ground-observed environmental data to predict the probability of established I. scapularis populations. The autologistic analysis showed that maximum, minimum, and mean temperatures as well as vapor pressure significantly contribute to population maintenance with an accuracy of 95% (p < 0.0001). A cutoff probability for habitat suitability was assessed by sensitivity analysis and was used to reclassify the previous distribution map. The spatially modeled relationship between I. scapularis presence and large-scale environmental data provides a robust suitability model that reveals essential environmental determinants of habitat suitability, predicts emerging areas of Lyme disease risk, and generates the future pattern of I. scapularis across the United States.     

Risk maps for range expansion of the Lyme disease vector, Ixodes scapularis, in Canada now and with climate change

Background

Lyme disease is the commonest vector-borne zoonosis in the temperate world, and an emerging infectious disease in Canada due to expansion of the geographic range of the tick vector Ixodes scapularis. Studies suggest that climate change will accelerate Lyme disease emergence by enhancing climatic suitability for I. scapularis. Risk maps will help to meet the public health challenge of Lyme disease by allowing targeting of surveillance and intervention activities.

Results

A risk map for possible Lyme endemicity was created using a simple risk algorithm for occurrence of I. scapularis populations. The algorithm was calculated for each census sub-division in central and eastern Canada from interpolated output of a temperature-driven simulation model of I. scapularis populations and an index of tick immigration. The latter was calculated from estimates of tick dispersion distances by migratory birds and recent knowledge of the current geographic range of endemic I. scapularis populations. The index of tick immigration closely predicted passive surveillance data on I. scapularis occurrence, and the risk algorithm was a significant predictor of the occurrence of I. scapularis populations in a prospective field study. Risk maps for I. scapularis occurrence in Canada under future projected climate (in the 2020s, 2050s and 2080s) were produced using temperature output from the Canadian Coupled Global Climate Model 2 with greenhouse gas emission scenario enforcing 'A2' of the Intergovernmental Panel on Climate Change.

Conclusion

We have prepared risk maps for the occurrence of I. scapularis in eastern and central Canada under current and future projected climate. Validation of the risk maps provides some confidence that they provide a useful first step in predicting the occurrence of I. scapularis populations, and directing public health objectives in minimizing risk from Lyme disease. Further field studies are needed, however, to continue validation and refinement of the risk maps.     

全沟硬蜱为内蒙大兴安岭莱姆病螺旋体的主要生物媒介

我们于1990年5～8月在内蒙大兴安岭,对莱姆病螺旋体的生物媒介和人群感染进行了调查,全沟硬蜱(Ixodes persulcatus)是当地蜱类的优势种,用直接荧光抗体法检查全沟硬蜱的中肠涂片,发现其带菌率为40%;从46组(380只)全沟硬蜱中分离出19株莱姆病螺旋体。分离出的菌株在单克隆反应上与美国菌株不同。用间接荧光抗体法检测308名居民血清,阳性率为13.3%。以上研究证明,大兴安岭是莱姆病疫区,全沟硬蜱是将莱姆病螺旋体传播到人和动物的主要生物媒介。     

The Ixodes ricinus tick as a vector of Borrelia burgdorferi in Slovenia

In many parts of Slovenia, the vegetation, climate and animals form an ideal natural environment for the persistence of different tick species. In this study, Ixodes ricinus ticks, the prevalent tick species in our country, are identified as the vectors of Borrelia burgdorferi in Slovenia. Ticks were collected by flagging between May and October 1990. All the ticks collected from different parts of Slovenia were I. ricinus and each one was examined individually. The ticks were tested to determine the presence of B. burgdorferi by indirect immunofluorescent assay. Among the 496 ticks examined for borreliae there were 85 adults and 411 nymphs. Twenty (23.5%) adults and 18 (4.4%) nymphs were found positive for Borreliae. This study discusses the importance of I. ricinus ticks in the role of Lyme borreliosis in Slovenia.     

Synchronous phenology of juvenile Ixodes scapularis,vertebrate host relationships,and associated patterns of Borrelia burgdorferi ribotypes in the midwestern United States

To elucidate features of enzootic maintenance of the Lyme disease bacterium that affect human risk of infection, we conducted a longitudinal study of the phenology of the vector tick, Ixodes scapularis, at a newly invaded site in the north-central United States. Surveys for questing ticks and ticks parasitizing white-footed mice and eastern chipmunks revealed that I. scapularis nymphal and larval activity peaked synchronously in June and exhibited an atypical, unimodal seasonality. Adult seasonal activity was bimodal and distributed evenly in spring and fall. We discuss implications of these phenology data for the duration of the I. scapularis life cycle. Densities of Borrelia burgdorferi-infected, questing nymphs were comparable to those found in endemic areas elsewhere in the midwestern and northeastern U.S. Molecular genetic diversity of B. burgdorferi infecting these ticks and rodents was assessed by analysis of the ribosomal spacer types (RSTs). RST 1, a clade that includes strains with highly pathogenic properties, was relatively uncommon (3.4%) in contrast to the northeastern U.S., whereas less pathogenic ribotypes of the RST 2 and 3 clades were more common. These features of the ecology of this midwestern Lyme disease system likely contribute to the lower incidence of Lyme disease in humans in the Upper Midwest compared with that of the Northeast owing to reduced exposure to pathogenic strains of B. burgdorferi.     

Deer management generally reduces densities of nymphal Ixodes scapularis,but not prevalence of infection with Borrelia burgdorferi sensu stricto

Human Lyme disease–primarily caused by the bacterium Borrelia burgdorferi sensu stricto (s.s.) in North America–is the most common vector-borne disease in the United States. Research on risk mitigation strategies during the last three decades has emphasized methods to reduce densities of the primary vector in eastern North America, the blacklegged tick (Ixodes scapularis). Controlling white-tailed deer populations has been considered a potential method for reducing tick densities, as white-tailed deer are important hosts for blacklegged tick reproduction. However, the feasibility and efficacy of white-tailed deer management to impact acarological risk of encountering infected ticks (namely, density of host-seeking infected nymphs; DIN) is unclear. We investigated the effect of white-tailed deer density and management on the density of host-seeking nymphs and B. burgdorferi s.s. infection prevalence using surveillance data from eight national parks and park regions in the eastern United States from 2014–2022. We found that deer density was significantly positively correlated with the density of nymphs (nymph density increased by 49% with a 1 standard deviation increase in deer density) but was not strongly correlated with the prevalence of B. burgdorferi s.s. infection in nymphal ticks. Further, while white-tailed deer reduction efforts were followed by a decrease in the density of I. scapularis nymphs in parks, deer removal had variable effects on B. burgdorferi s.s. infection prevalence, with some parks experiencing slight declines and others slight increases in prevalence. Our findings suggest that managing white-tailed deer densities alone may not be effective in reducing DIN in all situations but may be a useful tool when implemented in integrated management regimes.     

Lyme disease in Italy: Isolation of Borrelia burgdorferi from a patient

A strain of Borrelia burgdorferi was isolated from a 48-year-old female patient suffering from a chronic form of polyarthritis. Significant titers of specific anti-Borrelia antibodies were not found. This is the first report of isolation of Borrelia burgdorferi from a patient in Italy.     

Increased diversity of zoonotic pathogens and Borrelia burgdorferi strains in established versus incipient Ixodes scapularis populations across the Midwestern United States

The center of origin theory predicts that genetic diversity will be greatest near a specie’s geographic origin because of the length of time for evolution. By corollary, diversity will decrease with distance from the origin; furthermore, invasion and colonization are frequently associated with founder effects that reduce genetic variation in incipient populations. The blacklegged tick, Ixodes scapularis, harbors a suite of zoonotic pathogens, and the geographic range of the tick is expanding in the upper Midwestern United States. Therefore, we posited that diversity of I. scapularis-borne pathogens across its Midwestern range should correlate with the rate of the range expansion of this tick as well as subsequent disease emergence. Analysis of 1565 adult I. scapularis ticks from 13 sites across five Midwestern states revealed that tick infection prevalence with multiple microbial agents (Borrelia burgdorferi, Borrelia miyamotoi, Babesia odocoilei, Babesia microti, and Anaplasma phagocytophilum), coinfections, and molecular genetic diversity of B. burgdorferi all were positively correlated with the duration of establishment of tick populations, and therefore generally support the center of origin – pathogen diversity hypothesis. The observed differences across the gradient of establishment, however, were not strong and were nuanced by the high frequency of coinfections in tick populations at both established and recently-invaded tick populations. These results suggest that the invasion of ticks and their associated pathogens likely involve multiple means of pathogen introduction, rather than the conventionally presented scenario whereby infected, invading ticks are solely responsible for introducing pathogens to naïve host populations.     

Borrelia species in Ixodes affinis and Ixodes scapularis ticks collected from the coastal plain of North Carolina

Ixodes affinis and I. scapularis are tick species that are widely distributed in the coastal plain region of North Carolina. Both tick species are considered enzootic vectors for spirochetal bacteria of the genus Borrelia and specifically for B. burgdorferi s.s., the pathogen most often attributed as the cause of Lyme disease in the USA. Laboratory testing of individual I. affinis and I. scapularis ticks for the presence of Borrelia DNA was accomplished by PCR, targeting 2 regions of the 16S-23S intergenic spacer. In I. affinis, Borrelia DNA was detected in 63.2% of 155 individual ticks. B. burgdorferi s.s. and B. bissettii were identified by DNA sequencing in 33.5% and 27.9% I. affinis, respectively. Statistical differences were found for sex distribution of Borrelia DNA between I. affinis females (76.8%) and I. affinis males (55.6%) where B. burgdorferi s.s. was more prevalent in females (44.6%) than in males (27.3%). In I. scapularis, 298 individually tested ticks yielded no Borrelia PCR-positive results. This study found a higher incidence of Borrelia spp. in I. affinis collected in coastal North Carolina as compared to previous reports for this tick species in other Southern states, highlighting the potential importance of I. affinis in the maintenance of the enzootic transmission cycle of B. burgdorferi s.l. in North Carolina. The lack of Borrelia DNA in I. scapularis highlights the need for additional studies to better define the transmission cycle for B. burgdorferi s.s. in the southeastern USA and specifically in the state of North Carolina.     

Occurrence of Borrelia burgdorferi sensu stricto, Borrelia garinii and Borrelia afzelii in the Ixodes ricinus ticks from Eastern Slovakia

A total of 2816 unfed adults nymphs of Ixodes ricinus ticks were collected from vegetation in Koice (Eastern Slovakia) from 1994 to 1997. Prevalence of Borrelia burgdorferi s. 1. in I. ricinus ticks, detected by dark field microscopy, varied and depended upon the year and the habitat of the collected ticks. The lowest prevalence was observed in 1994 (4.8%). During 1995 it increased to 17.2% and during the next two years decreased to 15.5% and 14.2%. The rate of infection varied from 2.1 to 23.3% within 10 examined habitats of the Koice area. A different value of relative density of ticks was observed in various habitats. It ranged from 9–212 ticks per collecting hour within one flagged area (600m²) which is 1.5–35.5 ticks per 100m². Eight isolates were obtained from the infected ticks. Electrophoresis and immunoblotting with 6 monoclonal antibodies were used for the identification of Borrelia strains. Three tick isolates were identified as B. burgdorferi s. s. and the other three isolates were found to be B. garinii. One strain reacted as a mixed culture of B. burgdorferi s. s., and B. garinii. The strain originated from the Vihorlat Mountains habitat and was detected by PCR-SSCP as B. burgdorferi s. s. with a small amount of B. afzelii. The obtained results emphasize the epidemiological importance not only of B. garinii and B. afzelii but also of B. burgdorferi s. s. in Central Europe.     

Comparison of in vitro activities of tigecycline,doxycycline, and tetracycline against the spirochete Borrelia burgdorferi

Tigecycline is a new glycylcycline that has recently been revealed to be very effective in vitro against a variety of Gram-negative and Gram-positive bacteria including multi-drug resistant microorganisms. Using a standardized microdilution susceptibility testing method, we determined the minimal inhibitory concentrations (MICs) and the minimal bactericidal concentrations (MBCs) of tigecycline, in parallel with doxycycline, tetracycline, and other antibiotic agents relevant for Lyme borreliosis treatment such as ceftriaxone and cefotaxime. The activity of all agents against 16 different Borrelia isolates belonging to all borrelial genospecies known to be pathogenic for humans was investigated and analyzed under standardized conditions.The overall rank order of MIC₉₀s was tigecycline (≤0.016 mg/L) > ceftriaxone (0.03 mg/L) > cefotaxime (≤0.125 mg/L) > doxycycline (0.25 mg/L) > tetracycline (0.25 mg/L). The rank order of MBC₉₀s was tigecycline (0.5 mg/L) > ceftriaxone (2 mg/L) > tetracycline (16 mg/L) > doxycycline (16 mg/L) > cefotaxime (>16 mg/L). High in vitro activity of the new glycylcycline against Borrelia was further substantiated by time-kill experiments performed with B. afzelii isolate EB1. Parallel testing of tigecycline and ceftriaxone demonstrated a bacteriostatic effect for 0.016 mg/L of tigecycline and for 0.03 mg/L for ceftriaxone after 72 h of incubation. Moreover, tigecycline was bactericidal at a concentration of 0.25 mg/L showing a >3 log₁₀ unit reduction of the initial inoculum, whereas for ceftriaxone a concentration of 2 mg/L was needed.     

Predicting the risk of Lyme disease: habitat suitability for Ixodes scapularis in the north central United States

The distribution and abundance of Ixodes scapularis were studied in Wisconsin, northern Illinois, and portions of the Upper Peninsula of Michigan by inspecting small mammals for ticks and by collecting questing ticks at 138 locations in state parks and natural areas. Environmental data were gathered at a local level (i.e., micro and meso levels), and a geographic information system (GIS) was used with several digitized coverages of environmental data to create a habitat profile for each site and a grid map for Wisconsin and Illinois. Results showed that the presence and abundance of I. scapularis varied, even when the host population was adequate. Tick presence was positively associated with deciduous, dry to mesic forests and alfisol-type soils of sandy or loam-sand textures overlying sedimentary rock. Tick absence was associated with grasslands, conifer forests, wet to wet/mesic forests, acidic soils of low fertility and a clay soil texture, and Precambrian bedrock. We performed a discriminant analysis to determine environmental differences between positive and negative tick sites and derived a regression equation to examine the probability of I. scapularis presence per grid. Both analyses indicated that soil order and land cover were the dominant contributors to tick presence. We then constructed a risk map indicating suitable habitats within areas where I. scapularis is already established. The risk map also shows areas of high probability the tick will become established if introduced. Thus, this risk analysis has both explanatory power and predictive capability.