High burdens of Ixodes scapularis larval ticks on white-tailed deer may limit Lyme disease risk in a low biodiversity setting

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.

Host adaptation drives genetic diversity in a vector-borne disease system

The range of hosts a pathogen can infect is a key trait, influencing human disease risk and reservoir host infection dynamics. Borrelia burgdorferi sensu stricto (Bb), an emerging zoonotic pathogen, causes Lyme disease and is widely considered a host generalist, commonly infecting mammals and birds. Yet the extent of intraspecific variation in Bb host breadth, its role in determining host competence, and potential implications for human infection remain unclear. We conducted a long-term study of Bb diversity, defined by the polymorphic ospC locus, across white-footed mice, passerine birds, and tick vectors, leveraging long-read amplicon sequencing. Our results reveal strong variation in host breadth across Bb genotypes, exposing a spectrum of genotype-specific host-adapted phenotypes. We found support for multiple niche polymorphism, maintaining Bb diversity in nature and little evidence of temporal shifts in genotype dominance, as would be expected under negative frequency-dependent selection. Passerine birds support the circulation of several human-invasive strains (HISs) in the local tick population and harbor greater Bb genotypic diversity compared with white-footed mice. Mouse-adapted Bb genotypes exhibited longer persistence in individual mice compared with nonadapted genotypes. Genotype communities infecting individual mice preferentially became dominated by mouse-adapted genotypes over time. We posit that intraspecific variation in Bb host breadth and adaptation helps maintain overall species fitness in response to transmission by a generalist vector.

Tick Densities and Infection Prevalence on Coastal Islands in Massachusetts, USA: Establishing a Baseline

Tick-borne diseases and a tick-induced red meat allergy have become increasingly common in the northeastern USA and elsewhere. At the scale of local communities, few studies have documented tick densities or infection levels to characterize current conditions and provide a baseline for further monitoring. Using the town of Nantucket, MA, as a case study, we recorded tick densities by drag sampling along hiking trails in nature preserves on two islands. Nymphal blacklegged ticks (Ixodes scapularis Say) were most abundant at shadier sites and least common in grasslands and scrub oak thickets (Quercus ilicifolia). Lone star ticks (Amblyomma americanum L.) were common on Tuckernuck Island and rare on Nantucket Island, while both tick species were more numerous in 2021 compared to 2020 and 2022. We tested for pathogens in blacklegged nymphs at five sites over two years. In 2020 and 2021, infection levels among the four Nantucket Island sites averaged 10% vs. 19% for Borrelia burgdorferi, 11% vs. 15% for Babesia microti, and 17% (both years) for Anaplasma phagocytophilum, while corresponding levels were significantly greater on Tuckernuck in 2021. Our site-specific, quantitative approach represents a practical example of how potential exposure to tick-borne diseases can be monitored on a local scale.

Horses as Sentinels for the Circulation of Flaviviruses in Eastern-Central Germany

Since 2018, autochthonous West Nile virus (WNV) infections have been regularly reported in eastern-central Germany. While clinically apparent infections in humans and horses are not frequent, seroprevalence studies in horses may allow the tracing of WNV and related flaviviruses transmission, such as tick-borne encephalitis virus (TBEV) and Usutu virus (USUV), and consequently help to estimate the risk of human infections. Hence, the aim of our study was to follow the seropositive ratio against these three viruses in horses in Saxony, Saxony Anhalt, and Brandenburg and to describe their geographic distribution for the year 2021. In early 2022, i.e., before the virus transmission season, sera from 1232 unvaccinated horses were tested using a competitive pan-flavivirus ELISA (cELISA). In order to estimate the true seropositive ratio of infection with WNV, TBEV, and USUV for 2021, positive and equivocal results were confirmed by a virus neutralization test (VNT). In addition, possible risk factors for seropositivity using questionnaires were analyzed using logistic regression based on questionnaires similar to our previous study from 2020. In total, 125 horse sera reacted positive in the cELISA. Based on the VNT, 40 sera showed neutralizing antibodies against WNV, 69 against TBEV, and 5 against USUV. Three sera showed antibodies against more than one virus, and eight were negative based on the VNT. The overall seropositive ratio was 3.3% (95% CI: 2.38-4.40) for WNV, 5.6% (95% CI: 4.44-7.04) for TBEV, and 0.4% (95% CI: 0.14-0.98) for USUV infections. While age and number of horses on the holding were factors predicting TBEV seropositivity, no risk factors were discovered for WNV seropositivity. We conclude that horses are useful sentinels to determine the flavivirus circulation in eastern-central Germany, as long as they are not vaccinated against WNV.

Understanding Habitats and Environmental Conditions of White-Tailed Deer Population Density and Public Health Data to Aid in Assessing Human Tick-Borne Disease Risk

The extent of tick-borne diseases (TBDs) in the United States is largely unknown and underreported. Equitable diagnostic and treatment options may vary by geographic location. Triangulating multi-modal data sources informed by a One Health approach provides robust proxies for human TBD risk. Using data from the Indiana Department of Natural Resources collected from hunters during the white-tailed deer (Odocoileus virginianus) hunting season and other sources, we employ a mixed-methods approach based on thematic mapping and mixed effects modelling to determine if deer population density aligns with official disease data at the county level from (1) positive canine serological reports for, anaplasmosis, and Lyme Disease (LD); (2) positive human cases of ehrlichiosis, anaplasmosis, LD, and Spotted Fever rickettsioses; and (3) tick infectivity. We propose the need for multimodal data analysis using a variety of potential proxies to better estimate disease risk and inform public health policy and practice. We find similar spatial distributions between deer population density and human and canine TBDs in northeastern and southern Indiana, which are rural and mixed geographic areas. Overall, LD is more prevalent in the northwest, central-western, and southeastern counties, while ehrlichiosis is more common in the southern counties. These findings hold true across humans, canines, and deer.

Risk factors for Lyme disease: A scale-dependent effect of host species diversity and a consistent negative effect of host phylogenetic diversity

Biodiversity can influence disease risk. One example of a diversity-disease relationship is the dilution effect, which suggests higher host species diversity (often indexed by species richness) reduces disease risk. While numerous studies support the dilution effect, its generality remains controversial. Most studies of diversity-disease relationships have overlooked the potential importance of phylogenetic diversity. Furthermore, most studies have tested diversity-disease relationships at one spatial scale, even though such relationships are likely scale dependent. Using Lyme disease as a model system, we investigated the effects of host species richness and phylogenetic relatedness on the number of reported Lyme disease cases in humans in the U.S.A. at two spatial scales (the county level and the state level) using piecewise structural equation modelling. We also accounted for relevant climatic and habitat-related factors and tested their correlations with the number of Lyme disease cases. We found that species assemblages with more related species (i.e., host species in the order Rodentia) were associated with more Lyme disease cases in humans. Host species richness correlated negatively with the number of Lyme disease cases at the state level (i.e., a dilution effect), a pattern that might be explained by the higher number of reservoir-incompetent species at high levels of species richness at this larger spatial scale. In contrast, a positive correlation was found between species richness and the number of Lyme disease cases at the county level, where a higher proportion of rodent species was associated with higher levels of species richness, potentially amplifying the disease risk. Our results highlight that analyse at a single spatial scale can miss some impacts of biodiversity on human health. Thus, multi-scale analyses with consideration of host phylogenetic diversity are critical for improving our understanding of diversity-disease relationships.

An ecological approach to predict areas with established populations of Ixodes scapularis in Quebec, Canada

Public health management of Lyme disease (LD) is a dynamic challenge in Canada. Climate warming is driving the northward expansion of suitable habitat for the tick vector, Ixodes scapularis. Information about tick population establishment is used to inform the risk of LD but is challenged by sampling biases from surveillance data. Misclassifying areas as having no established tick population underestimates the LD risk classification. We used a logistic regression model at the municipal level to predict the probability of I. scapularis population establishment based on passive tick surveillance data during the period of 2010-2017 in southern Quebec. We tested for the effect of abiotic and biotic factors hypothesized to influence tick biology and ecology. Additional variables controlled for sampling biases in the passive surveillance data. In our final selected model, tick population establishment was positively associated with annual cumulative degree-days > 0°C, precipitation and deer density, and negatively associated with coniferous and mixed forest types. Sampling biases from passive tick surveillance were controlled for using municipal population size and public health instructions on tick submissions. The model performed well as indicated by an area under the curve (AUC) of 0.92, sensitivity of 86% and specificity of 81%. Our model enables prediction of I. scapularis population establishment in areas which lack data from passive tick surveillance and may improve the sensitivity of LD risk categorization in these areas. A more sensitive system of LD risk classification is important for increasing awareness and use of protective measures employed against ticks, and decreasing the morbidity associated with LD.

Patterns of white-tailed deer movements in suburban Maryland: implications for zoonotic disease mitigation

Understanding the ecology of the often dense white-tailed deer populations in urban and suburban landscapes is important for mitigating a variety of conflicts that arise with dense human populations, such as issues surrounding zoonotic disease mitigation. We collared white-tailed deer in highly suburban areas of Howard County, Maryland with high-resolution GPS collars. Then, we created autocorrelated kernel density home ranges for broader land use analyses and concurrently used general additive models to characterize fine-scale hourly measures of speed, activity, and proximity to residential buildings. Suburban deer home ranges encompassed approximately 35% direct residential land, and an average of 71 and 129 residential properties were found within female and male core ranges, respectively. Sex, time of day, and day of year all influenced fine-scale speeds, activity levels, and proximity to residential property buildings. Deer moved into residential areas nightly, especially in winter, and exhibited bouts of increased speed and activity shortly after sunrise and sunset, with distinctive seasonal changes. We discuss how variation in home ranges and movements may influence population management success and explore year-round periods of increased risk of deer transporting ticks to residential areas. These findings focus our broad understanding of deer movements in suburban and urban landscapes to improve deer population management and to mitigate the spread of ticks into residential areas.

Recreation reduces tick density through fine-scale risk effects on deer space-use

Altered interactions between pathogens, their hosts and vectors have potential consequences for human disease risk. Notably, tick-borne pathogens, many of which are associated with growing deer abundance, show global increasing prevalence and pose increasing challenges for disease prevention. Human activities can largely affect the patterns of deer space-use and can therefore be potential management tools to alleviate human-wildlife conflicts. Here, we tested how deer space-use patterns are influenced by human recreational activities, and how this in turn affects the spatial distribution of the sheep tick (Ixodes ricinus), a relevant disease vector of zoonoses such as Lyme borrelioses. We compared deer dropping and questing tick density on transects near (20 m) and further away from (100 m) forest trails that were either frequently used (open for recreation) or infrequently used (closed for recreation, but used by park managers). In contrast to infrequently used trails, deer dropping density was 31% lower near (20 m) than further away from (100 m) frequently used trails. Similarly, ticks were 62% less abundant near (20 m) frequently used trails compared to further away from (100 m) these trails, while this decline in tick numbers was only 14% near infrequently used trails. The avoidance by deer of areas close to human-used trails was thus associated with a similar reduction in questing tick density near these trails. As tick abundance generally correlates to pathogen prevalence, the use of trails for recreation may reduce tick-borne disease risk for humans on and near these trails. Our study reveals an unexplored effect of human activities on ecosystems and how this knowledge could be potentially used to mitigate zoonotic disease risk.

No net effect of host density on tick-borne disease hazard due to opposing roles of vector amplification and pathogen dilution

To better understand vector‐borne disease dynamics, knowledge of the ecological interactions between animal hosts, vectors, and pathogens is needed. The effects of hosts on disease hazard depends on their role in driving vector abundance and their ability to transmit pathogens. Theoretically, a host that cannot transmit a pathogen could dilute pathogen prevalence but increase disease hazard if it increases vector population size. In the case of Lyme disease, caused by Borrelia burgdorferi s.l. and vectored by Ixodid ticks, deer may have dual opposing effects on vectors and pathogen: deer drive tick population densities but do not transmit B. burgdorferi s.l. and could thus decrease or increase disease hazard. We aimed to test for the role of deer in shaping Lyme disease hazard by using a wide range of deer densities while taking transmission host abundance into account. We predicted that deer increase nymphal tick abundance while reducing pathogen prevalence. The resulting impact of deer on disease hazard will depend on the relative strengths of these opposing effects. We conducted a cross‐sectional survey across 24 woodlands in Scotland between 2017 and 2019, estimating host (deer, rodents) abundance, questing Ixodes ricinus nymph density, and B. burgdorferi s.l. prevalence at each site. As predicted, deer density was positively associated with nymph density and negatively with nymphal infection prevalence. Overall, these two opposite effects canceled each other out: Lyme disease hazard did not vary with increasing deer density. This demonstrates that, across a wide range of deer and rodent densities, the role of deer in amplifying tick densities cancels their effect of reducing pathogen prevalence. We demonstrate how noncompetent host density has little effect on disease hazard even though they reduce pathogen prevalence, because of their role in increasing vector populations. These results have implications for informing disease mitigation strategies, especially through host management.

Relevance of Spatial and Temporal Trends in Nymphal Tick Density and Infection Prevalence for Public Health and Surveillance Practice in Long-Term Endemic Areas: A Case Study in Monmouth County, NJ

Tick-borne diseases are a growing public health problem in the United States, and the US northeast has reported consistently high case rates for decades. Monmouth County, New Jersey, was one of the earliest jurisdictions to report Lyme disease cases in 1979 and reports several hundred cases per year nearly 40 yr later. In the time since, however, tick-borne health risks have expanded far beyond Lyme disease to include a variety of other bacterial pathogens and viruses, and additional vectors, necessitating a continually evolving approach to tick surveillance. In 2017, Monmouth County initiated an active surveillance program targeting sites across three ecological regions for collection of Ixodes scapularis Say (Acari: Ixodidae) and Amblyomma americanum L. (Acari: Ixodidae) as well as testing via qPCR for associated bacterial pathogens. During the first five years of this program (2017-2021), we report high levels of spatiotemporal variability in nymphal density and infection prevalence in both species, limiting the granularity with which human risk can be predicted from acarological data. Nonetheless, broader patterns emerged, including an ongoing trend of A. americanum dominance, risks posed by Borrelia miyamotoi, and the frequency of coinfected ticks. We present some of the first county-level, systematic surveillance of nymphal A. americanum density and infection prevalence in the northeastern US. We also documented a temporary decline in Borrelia burgdorferi that could relate to unmeasured trends in reservoir host populations. We discuss the implications of our findings for tick-borne disease ecology, public health communication, and tick surveillance strategies in endemic areas.

Effect of Urbanization on Presence, Abundance, and Coinfection of Bacteria and Protozoa in Ticks in the US Great Plains

Urbanization alters components of natural ecosystems which can affect tick abundance and tick-borne disease prevalence. Likely due to these changes, tick-borne pathogen prevalence has increased in many U.S. urban areas. Despite the growing public health importance of tick-borne diseases, little is known about how they are influenced by urbanization in North America, especially in the central U.S. where several pathogens occur at or near their highest levels of incidence nationally. To determine whether urban development influences tick infection with bacteria and protozoa, we collected ticks at 16 parks across a gradient of urbanization intensity in Oklahoma City, Oklahoma, USA and tested them using a variety of PCR assays. Adult ticks tested positive for Rickettsia parkeri, R. amblyommatis, R. rhiphicephali, 'Candidatus R. andeanae', Ehrlichia chaffeensis, E. ewingii, Panola Mountain Ehrlichia, 'Borrelia lonestari', Theileria cervi, Babesia spp. Coco, and Cytauxzoon felis. These results indicate the presence of a high diversity of tick-borne bacteria and protozoa across an expanding urban area in the U.S. Great Plains. Although there appeared to be some risk of encountering tick-borne microorganisms across the entire urbanization gradient, E. chaffeensis, E. ewingii, and T. cervi-infected ticks and microbe diversity decreased with increasing urbanization intensity. We identified a low rate of coinfection between different microorganisms, with coinfected ticks mainly collected from sites in the least-urbanized areas. This study suggests the need for awareness of tick-borne disease risk throughout urban areas in the central U.S., and highlights a need for studies of tick host habitat use and movement in cities.

A single tick screening for infectious pathogens using targeted mass spectrometry

The black-legged tick, Ixodes scapularis, is a well-known vector for the Lyme disease-causing pathogen (Borrelia burgdorferi) but can also carry other disease-causing pathogens such as Rickettsia, Anaplasma, Bartonella, Ehrlichia, and Theileria. Hence, tick screening using highly specific protein signatures for specific pathogens will help assess the prevalence of infected ticks and understand the pathogen-tick interactions in a specific geographic area. In this study, we used data-dependent acquisition to key pathogen protein signatures in black-legged ticks collected from the Southern Tier New York. Bottom-up proteomic analysis of extract from five combined ticks identified 2,052 tick proteins and 41 pathogen proteins with high confidence (≥ 99% C.I.). Results show high peptide spectral match counts for Rickettsia species and Borrelia species and lower counts for other rarer pathogens such as Anaplasma phagocytophilum. Parallel reaction monitoring performed on protein extracts from individual ticks (n = 10) revealed that 8 out of the 10 screened ticks carried Rickettsia species, 5 carried Borrelia species, 3 carried both pathogens, and only 1 tick carried no detectable bacteria. Mass spectrometry-based proteomics is a highly specific way to define the expression of different types of pathogen proteins in infected ticks. This might bring insights into the tick-pathogen interactions at the molecular level and especially expression pathogen surface proteins in ticks.

An eco-epidemiological modeling approach to investigate dilution effect in two different tick-borne pathosystems

Disease (re)emergence appears to be driven by biodiversity decline and environmental change. As a result, it is increasingly important to study host-pathogen interactions within the context of their ecology and evolution. The dilution effect is the concept that higher biodiversity decreases pathogen transmission. It has been observed especially in zoonotic vector-borne pathosystems, yet evidence against it has been found. In particular, it is still debated how the community (dis)assembly assumptions and the degree of generalism of vectors and pathogens affect the direction of the biodiversity-pathogen transmission relationship. The aim of this study was to use empirical data and mechanistic models to investigate dilution mechanisms in two rodent-tick-pathogen systems differing in their vector degree of generalism. A community was assembled to include ecological interactions that expand from purely additive to purely substitutive. Such systems are excellent candidates to analyze the link between vector ecology, community (dis)assembly dynamics, and pathogen transmission. To base our mechanistic models on empirical data, rodent live-trapping, including tick sampling, was conducted in Wales across two seasons for three consecutive years. We have developed a deterministic single-vector, multi-host compartmental model that includes ecological relationships with non-host species, uniquely integrating theoretical and observational approaches. To describe pathogen transmission across a gradient of community diversity, the model was populated with parameters describing five different scenarios differing in ecological complexity; each based around one of the pathosystems: Ixodes ricinus (generalist tick)-Borrelia burgdorferi and I. trianguliceps (small mammals specialist tick)-Babesia microti. The results suggested that community composition and interspecific dynamics affected pathogen transmission with different dilution outcomes depending on the vector degree of generalism. The model provides evidence that dilution and amplification effects are not mutually exclusive in the same community but depend on vector ecology and the epidemiological output considered (i.e., the "risk" of interest). In our scenarios, more functionally diverse communities resulted in fewer infectious rodents, supporting the dilution effect. In the pathosystem with generalist vector we identified a hump shaped relationship between diversity and infections in hosts, while for that characterized by specialist tick, this relationship was more complex and more dependent upon specific parameter values.

Limited capacity of deer to serve as zooprophylactic hosts for Borrelia burgdorferi in northeastern United States

Because deer are considered to be incompetent reservoirs of the agent of Lyme disease (Borrelia burgdorferi sensu stricto) in the northeastern U.S., they may serve as zooprophylactic or "dilution" hosts if larvae of the deer tick vector (Ixodes dammini, "northern" clade of Ixodes scapularis) frequently feed on them. To determine whether host-seeking nymphal deer ticks commonly feed on deer as larvae, we used a real time PCR host bloodmeal remnant identification assay to identify the host on which these ticks had fed. Nymphal Lone star ticks (Amblyomma americanum) were collected simultaneously in our sites and provided an index of the availability of deer in these sites. At 3 of the 4 sites, Ixodes nymphs had fed as larvae on a variety of hosts, including mice, birds and shrews, but rarely on deer (<6% for all sites); in contrast, Lone star tick nymphs had commonly fed on deer (31-78%). Deer were common larval hosts for Ixodes ticks (39% of bloodmeals) in only one site. The prevalence of B. burgdorferi in host seeking nymphal deer ticks was associated with mouse-fed ticks (p=0.007) but there was no association with deer-fed ticks (p=0.5). The diversity and prevalence of hosts that were identified differed between deer ticks and Lone star ticks that were collected simultaneously, demonstrating that there is no confounding of host bloodmeal identification by contaminating environmental DNA (eDNA). We conclude that deer were not common hosts for larval deer ticks, thus limiting their zooprophylactic role in our sites. Importance: Because deer are incompetent reservoirs for B. burgdorferi, their presence may modulate the force of enzootic transmission by serving as zoophrophylatic or "dilution" hosts. Such an effect would depend on the extent to which subadult deer ticks feed on other hosts. We used bloodmeal analysis on nymphal deer ticks to identify the host upon which larvae had fed. We found that Lone star ticks collected at the same time as deer ticks commonly fed on deer, but deer ticks did not. We conclude that deer are not a preferred host for larval deer ticks and thus are not necessarily zooprophylactic.

Tick infestation effects on haemoglobin levels of deer mice (Peromyscus maniculatus)

Deer mice (Peromyscus maniculatus) are hosts to ixodid ticks as well as the associated tick-borne pathogens they can spread. As the ranges of black-legged ticks (Ixodes scapularis) and American dog ticks (Dermacentor variabilis) expand northwards, naïve host populations of deer mice are likely to become infested by ticks and experience the physiological effects that ticks can have on them via blood-feeding. The prevalence of these haematophagous ticks can affect the haemoglobin levels of the mice they infest. Haemoglobin levels were compared and analysed in deer mice populations at three different sites with varying tick exposure. These results suggested that without confounding effects, the abundance of black-legged and American dog ticks on individual mice had a significant negative effect on the hosts' haemoglobin levels, but only in an area with high tick infestation. This was seen across the average haemoglobin levels between populations, where there was a significant difference between the source population with the longest established tick populations and the source population where neither black-legged nor American dog ticks were prevalent. As the ticks' ranges expand and they become more abundant, it is important to understand how their prevalence and intensity can alter host physiology, potentially affecting their own range expansion and the spread of the diseases they may carry.

Is composition of vertebrates an indicator of the prevalence of tick-borne pathogens?

Communities of vertebrates tend to appear together under similar ranges of environmental features. This study explores whether an explicit combination of vertebrates and their contact rates with a tick vector might constitute an indicator of the prevalence of a pathogen in the quest for ticks at the western Palearctic scale. We asked how ‘indicator’ communities could be ‘markers’ of the actual infection rates of the tick in the field of two species of Borrelia (a bacterium transmitted by the tick Ixodes ricinus). We approached an unsupervised classification of the territory to obtain clusters on the grounds of abundance of each vertebrate and contact rates with the tick. Statistical models based on Neural Networks, Random Forest, Gradient Boosting, and AdaBoost were detect the best correlation between communities’ composition and the prevalence of Borrelia afzelii and Borrelia gariniii in questing ticks. Both Gradient Boosting and AdaBoost produced the best results, predicting tick infection rates from the indicator communities. A ranking algorithm demonstrated that the prevalence of these bacteria in the tick is correlated with indicator communities of vertebrates on sites selected as a proof-of-concept. We acknowledge that our findings are supported by statistical outcomes, but they provide consistency for a framework that should be deeper explored at the large scale.

The bacterial biome of ticks and their wildlife hosts at the urban-wildland interface

Advances in sequencing technologies have revealed the complex and diverse microbial communities present in ticks (Ixodida). As obligate blood-feeding arthropods, ticks are responsible for a number of infectious diseases that can affect humans, livestock, domestic animals and wildlife. While cases of human tick-borne diseases continue to increase in the northern hemisphere, there has been relatively little recognition of zoonotic tick-borne pathogens in Australia. Over the past 5 years, studies using high-throughput sequencing technologies have shown that Australian ticks harbour unique and diverse bacterial communities. In the present study, free-ranging wildlife (n=203), representing ten mammal species, were sampled from urban and peri-urban areas in New South Wales (NSW), Queensland (QLD) and Western Australia (WA). Bacterial metabarcoding targeting the 16S rRNA locus was used to characterize the microbiomes of three sample types collected from wildlife: blood, ticks and tissue samples. Further sequence information was obtained for selected taxa of interest. Six tick species were identified from wildlife: Amblyomma triguttatum, Ixodes antechini, Ixodes australiensis, Ixodes holocyclus, Ixodes tasmani and Ixodes trichosuri. Bacterial 16S rRNA metabarcoding was performed on 536 samples and 65 controls, generating over 100 million sequences. Alpha diversity was significantly different between the three sample types, with tissue samples displaying the highest alpha diversity (P<0.001). Proteobacteria was the most abundant taxon identified across all sample types (37.3 %). Beta diversity analysis and ordination revealed little overlap between the three sample types (P<0.001). Taxa of interest included Anaplasmataceae, Bartonella, Borrelia, Coxiellaceae, Francisella, Midichloria, Mycoplasma and Rickettsia. Anaplasmataceae bacteria were detected in 17.7% (95/536) of samples and included Anaplasma, Ehrlichia and Neoehrlichia species. In samples from NSW, 'Ca. Neoehrlichia australis', 'Ca. Neoehrlichia arcana', Neoehrlichia sp. and Ehrlichia sp. were identified. A putative novel Ehrlichia sp. was identified from WA and Anaplasma platys was identified from QLD. Nine rodent tissue samples were positive for a novel Borrelia sp. that formed a phylogenetically distinct clade separate from the Lyme Borrelia and relapsing fever groups. This novel clade included recently identified rodent-associated Borrelia genotypes, which were described from Spain and North America. Bartonella was identified in 12.9% (69/536) of samples. Over half of these positive samples were obtained from black rats (Rattus rattus), and the dominant bacterial species identified were Bartonella coopersplainsensis and Bartonella queenslandensis. The results from the present study show the value of using unbiased high-throughput sequencing applied to samples collected from wildlife. In addition to understanding the sylvatic cycle of known vector-associated pathogens, surveillance work is important to ensure preparedness for potential zoonotic spillover events.

Evidence for Vertical Transmission of Babesia odocoilei (Piroplasmida: Babesiidae) in Ixodes scapularis (Acari: Ixodidae)

Limited evidence suggests that the cervid parasite, Babesia odocoilei, is transovarially transmitted from adult female Ixodes scapularis Say to offspring. The prevalence of B. odocoilei in unfed larval I. scapularis and whether vertical transmission is crucial to pathogen maintenance are currently unknown. To investigate these questions, 275 unfed larvae from two Wisconsin counties were tested for B. odocoilei genetic material. Sixteen of 29 pools were positive for the parasite. The maximum likelihood estimation for overall larval infection prevalence was 7.8% (95% confidence interval: 4.7-12). This vertically acquired infection appears to be sustained transstadially in nymphal ticks the following year; however, our relatively small sample and replicate size warrants additional evaluation. Our study revealed further evidence of vertical transmission, a low and consistent infection prevalence in larvae, and the potential importance of vertical transmission in B. odocoilei maintenance.

Research effort on birds' reservoir host potential for Lyme borreliosis: A systematic review and perspectives

Zoonotic tick-borne diseases threat human and animal health. Understanding the role of hosts in the production of infected ticks in an epidemiological system is essential to be able to design effective measures to reduce the exposure of humans and animals to infectious tick bites. The reservoir host potential, that is, number of infected ticks produced by a host species, depends on three components: tick production, realized reservoir competence and host density. The parameters and factors that determine the reservoir host potential need to be characterized to achieve a robust understanding of the dynamics of pathogen-tick-host systems, and thus to mitigate the acarological risk of emerging infections. Few studies have investigated the role of birds in the local spread of Lyme borreliosis Borrelia. Knowledge of the research effort on the reservoir host potential of birds in Lyme borreliosis Borrelia circulation is necessary to prioritize future research on this topic. We provide a systematic review of the research effort on components of the reservoir host potential of wild birds for Lyme borreliosis Borrelia circulation, and factors that modulate these components in the European epidemiological system. Our review of 242 selected publications showed that tick production has been 1.4 and 21 times more studied than realized reservoir competence and bird density respectively. Only one study achieved to characterize the global host reservoir potential of birds in a given epidemiological system. Investigated factors were mostly related to bird species identity, individual characteristics of birds and tick characteristics, whereas the influence of bird life-history traits have been largely under-investigated. Because simultaneous characterization of all parameters is notoriously complex, interdisciplinary research is needed to combine and accumulate independent field and laboratory investigations targeting each parameter on specific epidemiological system or host species. This can help gain an integrated appraisal of the functioning of the studied system at a local scale.

Coinfection of Ixodes scapularis (Acari: Ixodidae) Nymphs With Babesia spp. (Piroplasmida: Babesiidae) and Borrelia burgdorferi (Spirochaetales: Spirochaetaceae) in Wisconsin

Borrelia burgdorferi, the spirochete that causes Lyme disease, is endemic and widespread in Wisconsin. Research in the northeastern United States has revealed a positive association between Babesia microti, the main pathogen that causes babesiosis in humans, and Bo. burgdorferi in humans and in ticks. This study was conducted to examine associations between the disease agents in the Upper midwestern United States. Ixodes scapularis Say nymphs (N = 2,858) collected between 2015 and 2017 from nine locations in Wisconsin were tested for Babesia spp. and Borrelia spp. using real-time PCR. Two species of Babesia were detected; Ba. microti and Babesia odocoilei (a parasite of members of the family Cervidae). Prevalence of infection at the nine locations ranged from 0 to 13% for Ba. microti, 11 to 31% for Bo. burgdorferi sensu stricto, and 5.7 to 26% for Ba. odocoilei. Coinfection of nymphs with Bo. burgdorferi and Ba. odocoilei was detected in eight of the nine locations and significant positive associations were observed in two of the eight locations. The prevalence of nymphal coinfection with both and Bo. burgdorferi and Ba. microti ranged from 0.81 to 6.5%. These two pathogens were significantly positively associated in one of the five locations where both pathogens were detected. In the other four locations, the observed prevalence of coinfection was higher than expected in all but one site-year. Clinics and healthcare providers should be aware of the association between Ba. microti and Bo. burgdorferi pathogens when treating patients who report tick bites.

Impact of Land Use Changes and Habitat Fragmentation on the Eco-epidemiology of Tick-Borne Diseases

The incidence of tick-borne diseases has increased in recent decades and accounts for the majority of vector-borne disease cases in temperate areas of Europe, North America, and Asia. This emergence has been attributed to multiple and interactive drivers including changes in climate, land use, abundance of key hosts, and people's behaviors affecting the probability of human exposure to infected ticks. In this forum paper, we focus on how land use changes have shaped the eco-epidemiology of Ixodes scapularis-borne pathogens, in particular the Lyme disease spirochete Borrelia burgdorferi sensu stricto in the eastern United States. We use this as a model system, addressing other tick-borne disease systems as needed to illustrate patterns or processes. We first examine how land use interacts with abiotic conditions (microclimate) and biotic factors (e.g., host community composition) to influence the enzootic hazard, measured as the density of host-seeking I. scapularis nymphs infected with B. burgdorferi s.s. We then review the evidence of how specific landscape configuration, in particular forest fragmentation, influences the enzootic hazard and disease risk across spatial scales and urbanization levels. We emphasize the need for a dynamic understanding of landscapes based on tick and pathogen host movement and habitat use in relation to human resource provisioning. We propose a coupled natural-human systems framework for tick-borne diseases that accounts for the multiple interactions, nonlinearities and feedbacks in the system and conclude with a call for standardization of methodology and terminology to help integrate studies conducted at multiple scales.

Dilution and amplification effects in Lyme disease: Modeling the effects of reservoir-incompetent hosts on Borrelia burgdorferi sensu stricto transmission

The literature on Lyme disease includes a lively debate about the paradoxical role of changing deer populations. A decrease in the number of deer will both (1) reduce the incidence of Lyme disease by decreasing the host populations for ticks and therefore tick populations, and (2) enhance the incidence of Lyme disease by offering fewer reservoir-incompetent hosts for ticks, forcing the vector to choose reservoir-competent, and therefore possibly diseased, hosts to feed on. A review of field studies exploring the net impact of changing deer populations shows mixed results. In this manuscript, we investigate the hypothesis that the balance of these two responses to changing deer populations depends on the relative population sizes of reservoir-competent vs. reservoir-incompetent hosts and the presence of host preference in larval and adult stages. A temperature driven seasonal model of Borrelia burgdorferi sensu stricto (cause of Lyme disease) transmission among three host types (reservoir-competent infected and uninfected hosts, and reservoir-incompetent hosts) is constructed as a system of nonlinear ordinary differential equations. The model, which produces biologically reasonable results for both the tick vector Ixodes scapularis Say 1921 and the hosts, is used to investigate the effects of reservoir-incompetent host removal on both tick populations and disease prevalence for various relative population sizes of reservoir-competent hosts vs. reservoir-incompetent hosts. In summary, the simulation results show that the model with host preference appears to be more accurate than the one with no host preference. Given these results, we found that removal of adult I. scapularis(Say) hosts is likely to reduce questing nymph populations. At very low levels questing adult abundance may rise with lack of adult hosts. There is a dilution effect at low reservoir-competent host populations and there is an amplification effect at high reservoir-competent host populations.

Masting by beech trees predicts the risk of Lyme disease

BACKGROUND

The incidence of Lyme borreliosis and other tick-borne diseases is increasing in Europe and North America. There is currently much interest in identifying the ecological factors that determine the density of infected ticks as this variable determines the risk of Lyme borreliosis to vertebrate hosts, including humans. Lyme borreliosis is caused by the bacterium Borrelia burgdorferi sensu lato (s.l.) and in western Europe, the hard tick Ixodes ricinus is the most important vector.

METHODS

Over a 15-year period (2004-2018), we monitored the monthly abundance of I. ricinus ticks (nymphs and adults) and their B. burgdorferi s.l. infection status at four different elevations on a mountain in western Switzerland. We collected climate variables in the field and from nearby weather stations. We obtained data on beech tree seed production (masting) from the literature, as the abundance of Ixodes nymphs can increase dramatically 2 years after a masting event. We used generalized linear mixed effects models and AIC-based model selection to identify the ecological factors that influence inter-annual variation in the nymphal infection prevalence (NIP) and the density of infected nymphs (DIN).

RESULTS

We found that the NIP decreased by 78% over the study period. Inter-annual variation in the NIP was explained by the mean precipitation in the present year, and the duration that the DNA extraction was stored in the freezer prior to pathogen detection. The DIN decreased over the study period at all four elevation sites, and the decrease was significant at the top elevation. Inter-annual variation in the DIN was best explained by elevation site, year, beech tree masting index 2 years prior and the mean relative humidity in the present year. This is the first study in Europe to demonstrate that seed production by deciduous trees influences the density of nymphs infected with B. burgdorferi s.l. and hence the risk of Lyme borreliosis.

CONCLUSIONS

Public health officials in Europe should be aware that masting by deciduous trees is an important predictor of the risk of Lyme borreliosis.

Ixodes scapularis (Acari: Ixodidae) Nymphal Survival and Host-Finding Success in the Eastern United States

The blacklegged tick (Ixodes scapularis Say) is the primary vector of Borrelia burgdorferi sensu stricto (Spirochaetales: Spirochaetaceae), the Lyme disease agent in North America. The basic reproduction number (R0) for B. burgdorferi in I. scapularis in the Northeast is highly sensitive to the probability that engorged larvae survive the winter, molt into nymphs, and find a host. These processes are dependent on local environmental variables, including climate, host population size and movement, and tick behavior. A simple model is presented for estimating host-finding success from the ratio of tick abundance in two subsequent years, accounting for overwinter survival and possible differences in host associations between nymphs and larvae. This model was parameterized using data from two sites in mainland Connecticut and two on Block Island, RI. Host abundance and tick burdens were estimated via mark-recapture trapping of the primary host, Peromyscus leucopus Rafinesque. Overwintering survival was estimated using engorged larvae placed in field enclosures at each site. Only nymphs were recovered alive, and no significant differences in model parameters were observed between Connecticut and Block Island. Host-finding success was predicted to be high across a wide range of host association patterns at three of four sites. Assuming equivalent host association between larvae and nymphs, R0 was also estimated to be greater than one at three of four sites, suggesting these conditions allow for the persistence of B. burgdorferi. The model output was highly sensitive to differences between nymphal and larval host associations.

Game Animal Density, Climate, and Tick-Borne Encephalitis in Finland, 2007-2017

Tick-borne encephalitis (TBE) is an endemic infection of public health importance in Finland. We investigated the effect of ecologic factors on 2007-2017 TBE trends. We obtained domestic TBE case data from the National Infectious Diseases Register, weather data from the US National Oceanic and Atmospheric Administration, and data from the Natural Resources Institute in Finland on mammals killed by hunters yearly in game management areas. We performed a mixed-effects time-series analysis with time lags on weather and animal parameters, adding a random effect to game management areas. During 2007-2017, a total of 395/460 (86%) domestic TBE cases were reported with known place of exposure and date of sampling. Overall, TBE incidence increased yearly by 15%. After adjusting for the density of other animals and minimum temperatures, we found thatTBE incidence was positively associated with white-tailed deer density. Variation in host animal density should be considered when assessing TBE risks and designing interventions.

Lactobacilli and other gastrointestinal microbiota of Peromyscus leucopus, reservoir host for agents of Lyme disease and other zoonoses in North America

The cricetine rodent Peromyscus leucopus is an important reservoir for several human zoonoses, including Lyme disease, in North America. Akin to hamsters, the white-footed deermouse has been unevenly characterized in comparison to the murid Mus musculus. To further understanding of P. leucopus' total genomic content, we investigated gut microbiomes of an outbred colony of P. leucopus, inbred M. musculus, and a natural population of P. leucopus. Metagenome and whole genome sequencing were combined with microbiology and microscopy approaches. A focus was the genus Lactobacillus, four diverse species of which were isolated from forestomach and feces of colony P. leucopus. Three of the species-L. animalis, L. reuteri, and provisionally-named species "L. peromysci"-were identified in fecal metagenomes of wild P. leucopus but not discernibly in samples from M. musculus. L. johnsonii, the fourth species, was common in M. musculus but absent or sparse in wild P. leucopus. Also identified in both colony and natural populations were a Helicobacter sp. in feces but not stomach, and a Tritrichomonas sp. protozoan in cecum or feces. The gut metagenomes of colony P. leucopus were similar to those of colony M. musculus at the family or higher level and for major subsystems. But there were multiple differences between species and sexes within each species in their gut metagenomes at orthologous gene level. These findings provide a foundation for hypothesis-testing of functions of individual microbial species and for interventions, such as bait vaccines based on an autochthonous bacterium and targeting P. leucopus for transmission-blocking.

LYMESIM 2.0: An Updated Simulation of Blacklegged Tick (Acari: Ixodidae) Population Dynamics and Enzootic Transmission of Borrelia burgdorferi (Spirochaetales: Spirochaetaceae)

Lyme disease is the most commonly reported vector-borne disease in the United States, and the number of cases reported each year continues to rise. The complex nature of the relationships between the pathogen (Borrelia burgdorferi sensu stricto), the tick vector (Ixodes scapularis Say), multiple vertebrate hosts, and numerous environmental factors creates challenges for understanding and predicting tick population and pathogen transmission dynamics. LYMESIM is a mechanistic model developed in the late 1990s to simulate the life-history of I. scapularis and transmission dynamics of B. burgdorferi s.s. Here we present LYMESIM 2.0, a modernized version of LYMESIM, that includes several modifications to enhance the biological realism of the model and to generate outcomes that are more readily measured under field conditions. The model is tested for three geographically distinct locations in New York, Minnesota, and Virginia. Model-simulated timing and densities of questing nymphs, infected nymphs, and abundances of nymphs feeding on hosts are consistent with field observations and reports for these locations. Sensitivity analysis highlighted the importance of temperature in host finding for the density of nymphs, the importance of transmission from small mammals to ticks on the density of infected nymphs, and temperature-related tick survival for both density of nymphs and infected nymphs. A key challenge for accurate modeling of these metrics is the need for regionally representative inputs for host populations and their fluctuations. LYMESIM 2.0 is a useful public health tool that downstream can be used to evaluate tick control interventions and can be adapted for other ticks and pathogens.

Local abundance of Ixodes scapularis in forests: Effects of environmental moisture, vegetation characteristics, and host abundance

Ixodes scapularis is the primary vector of Lyme disease spirochetes in eastern and central North America, and local densities of this tick can affect human disease risk. We sampled larvae and nymphs from sites in Massachusetts and Wisconsin, USA, using flag/drag devices and by collecting ticks from hosts, and measured environmental variables to evaluate the environmental factors that affect local distribution and abundance of I. scapularis. Our sites were all forested areas with known I. scapularis populations. Environmental variables included those associated with weather (e.g., temperature and relative humidity), vegetation characteristics (at canopy, shrub, and ground levels), and host abundance (small and medium-sized mammals and reptiles). The numbers of larvae on animals at a given site and season showed a logarithmic relationship to the numbers in flag/drag samples, suggesting limitation in the numbers on host animals. The numbers of nymphs on animals showed no relationship to the numbers in flag/drag samples. These results suggest that only a small proportion of larvae and nymphs found hosts because in neither stage did the numbers of host-seeking ticks decline with increased numbers on hosts. Canopy cover was predictive of larval and nymphal numbers in flag/drag samples, but not of numbers on hosts. Numbers of small and medium-sized mammal hosts the previous year were generally not predictive of the current year's tick numbers, except that mouse abundance predicted log numbers of nymphs on all hosts the following year. Some measures of larval abundance were predictive of nymphal numbers the following year. The mean number of larvae per mouse was well predicted by measures of overall larval abundance (based on flag/drag samples and samples from all hosts), and some environmental factors contributed significantly to the model. In contrast, the mean numbers of nymphs per mouse were not well predicted by environmental variables, only by overall nymphal abundance on hosts. Therefore, larvae respond differently than nymphs to environmental factors. Furthermore, flag/drag samples provide different information about nymphal numbers than do samples from hosts. Flag/drag samples can provide information about human risk of acquiring nymph-borne pathogens because they provide information on the densities of ticks that might encounter humans, but to understand the epizootiology of tick-borne agents both flag/drag and host infestation data are needed.

Genetically Engineering Wild Mice to Combat Lyme Disease: An Ecological Perspective

Genetic engineering of wild populations has been proposed for reducing human diseases by altering pathogens’ hosts. For example, CRISPR-based genome editing may be used to create white-footed mice (Peromyscus leucopus) that are resistant to the Lyme disease spirochete vectored by blacklegged ticks (Ixodes scapularis). Toward this goal, academic researchers are developing Lyme-resistant and tick-resistant white-footed mice, which are a primary pathogen reservoir for Lyme disease in the United States. If field trials on small, experimental islands are successful, the project would scale up to the larger islands of Nantucket and Martha's Vineyard, Massachusetts, and possibly to the mainland, most likely with a local gene drive to speed the traits’ proliferation, pending approvals from relevant constituents. Despite considerable publicity, this project has yet to be evaluated by independent professional ecologists. In the present article, I discuss key ecological and evolutionary questions that should be considered before such genetically engineered mice are released into natural habitats.