The effect of spatial heterogenity on the aggregation of ticks on white-footed mice

Parasite infestation patterns differ between ticks and chigger mites on two rodent host species in Taiwan

Parasites are typically concentrated on a few host individuals, and identifying the mechanisms underlying aggregated distribution can facilitate a more targeted control of parasites. We investigated the infestation patterns of hard ticks and chigger mites on two rodent species, the striped field mouse, Apodemus agrarius, and the lesser ricefield rat, Rattus losea, in Taiwan. We also explored abiotic and biotic factors that were important in explaining variation in the abundance of ticks and chiggers on rodent hosts. Ticks were more aggregated than chiggers on both rodent species. Factors important for the variation in parasitic loads, especially biotic factors, largely differed between ticks and chiggers. Variation partitioning analyses revealed that a larger proportion of variation in chiggers than in ticks can be explained, especially by abiotic factors. If, as proposed, the higher number of parasites in males is due to a larger range area or immunity being suppressed by testosterone, when A. agrarius males host more ticks, they are expected to also host more chiggers, given that chiggers adopt a similar host finding approach to that of ticks. Instead, the similar abundance of chiggers in male and female A. agrarius implies that a large home range or suppressed immunity does not predispose males to inevitably host more parasites. More variations were explained by abiotic than biotic factors, suggesting that controlling practices are more likely to be successful by focusing on factors related to the environment instead of host traits. Our study indicated that the extent of parasitism is rarely determined by a sole factor, but is an outcome of complex interactions among animal physiology, animal behavior, characteristics of parasites, and the environments.

Individual heterogeneity in ixodid tick infestation and prevalence of Borrelia burgdorferi sensu lato in a northern community of small mammalian hosts

Heterogeneous aggregation of parasites between individual hosts is common and regarded as an important factor in understanding transmission dynamics of vector-borne diseases. Lyme disease is vectored by generalist tick species, yet we have a limited understanding of how individual heterogeneities within small mammal host populations affect the aggregation of ticks and likelihood of infection. Male hosts often have higher parasite and infection levels than females, but whether this is linked to sexual body size dimorphism remains uncertain. Here, we analysed how host species, sex, and body mass influenced Ixodes ricinus tick infestations and the infection prevalence of Borrelia burgdorferi sensu lato (s.l.) in three species of small mammals involved in the enzootic transmission cycle of Lyme disease in Norway from 2018 to 2022. Larval and nymphal ticks were found on 98% and 34% of all individual hosts, respectively. In bank voles and wood mice, both larval and nymphal tick infestation and infection probability increased with body mass, and it increased more with mass for males than for females. Tick infestation in the common shrew increased with body mass and was higher in males, while pathogen infection was higher in females. Sex-biases in infestation did not correspond with level of sexual body mass dimorphism across species. This study contributes to our understanding of how individual heterogeneity among small mammalian hosts influences I. ricinus tick aggregation and prevalence of B. burgdorferi s.l. at northern latitudes.

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.

Variation among strains of Borrelia burgdorferi in host tissue abundance and lifetime transmission determine the population strain structure in nature

Pathogen life history theory assumes a positive relationship between pathogen load in host tissues and pathogen transmission. Empirical evidence for this relationship is surprisingly rare due to the difficulty of measuring transmission for many pathogens. The comparative method, where a common host is experimentally infected with a set of pathogen strains, is a powerful approach for investigating the relationships between pathogen load and transmission. The validity of such experimental estimates of strain-specific transmission is greatly enhanced if they can predict the pathogen population strain structure in nature. Borrelia burgdorferi is a multi-strain, tick-borne spirochete that causes Lyme disease in North America. This study used 11 field-collected strains of B. burgdorferi, a rodent host (Mus musculus, C3H/HeJ) and its tick vector (Ixodes scapularis) to determine the relationship between pathogen load in host tissues and lifetime host-to-tick transmission (HTT). Mice were experimentally infected via tick bite with 1 of 11 strains. Lifetime HTT was measured by infesting mice with I. scapularis larval ticks on 3 separate occasions. The prevalence and abundance of the strains in the mouse tissues and the ticks were determined by qPCR. We used published databases to obtain estimates of the frequencies of these strains in wild I. scapularis tick populations. Spirochete loads in ticks and lifetime HTT varied significantly among the 11 strains of B. burgdorferi. Strains with higher spirochete loads in the host tissues were more likely to infect feeding larval ticks, which molted into nymphal ticks that had a higher probability of B. burgdorferi infection (i.e., higher HTT). Our laboratory-based estimates of lifetime HTT were predictive of the frequencies of these strains in wild I. scapularis populations. For B. burgdorferi, the strains that establish high abundance in host tissues and that have high lifetime transmission are the strains that are most common in nature.

Borrelia burgdorferi strain and host sex influence pathogen prevalence and abundance in the tissues of a laboratory rodent host

Experimental infections with different pathogen strains give insight into pathogen life history traits. The purpose of the present study was to compare variation in tissue infection prevalence and spirochete abundance among strains of Borrelia burgdorferi in a rodent host (Mus musculus, C3H/HeJ). Male and female mice were experimentally infected via tick bite with one of 12 strains. Ear tissue biopsies were taken at days 29, 59 and 89 postinfection, and seven tissues were collected at necropsy. The presence and abundance of spirochetes in the mouse tissues were measured by quantitative polymerase chain reaction. To determine the frequencies of our strains in nature, their multilocus sequence types were matched to published data sets. For the infected mice, 56.6% of the tissues were infected with B. burgdorferi. The mean spirochete load in the mouse necropsy tissues varied 4.8-fold between the strains. The mean spirochete load in the ear tissue biopsies decreased rapidly over time for some strains. The percentage of infected tissues in male mice (65.4%) was significantly higher compared to female mice (50.5%). The mean spirochete load in the seven tissues was 1.5× higher in male mice compared to female mice; this male bias was 15.3× higher in the ventral skin. Across the 11 strains, the mean spirochete loads in the infected mouse tissues were positively correlated with the strain-specific frequencies in their tick vector populations. The study suggests that laboratory-based estimates of pathogen abundance in host tissues can predict the strain composition of this important tick-borne pathogen in nature.

Clobetasol increases the abundance of Borrelia burgdorferi in the skin 70 times more in male mice compared to female mice

Lyme borreliosis is caused by the spirochete Borrelia burgdorferi and is transmitted among vertebrate hosts by Ixodes scapularis ticks in eastern North America. Treatment with topical corticosteroids increases the abundance of B. burgdorferi in the skin of lab mice that have been experimentally infected via needle inoculation. In the present study, female and male C3H/HeJ mice were infected with B. burgdorferi via nymphal tick bite. Infected mice were treated with clobetasol on the skin of the right hindleg on days 35 and 36 post-infection and euthanized at days -2, 1, 3, 5, and 7 post-treatment; a group of control mice was infected but not treated with clobetasol. The spirochete abundance was quantified in 8 mouse tissues including bladder, heart, left hindleg skin, right hindleg skin, dorsal skin, ventral skin, left ear and right ear. Averaged across the 8 mouse tissues, the abundance of B. burgdorferi on days 3 and 5 were 21.4x and 14.4x higher in mice treated with clobetasol compared to the untreated control mice, but there were large differences among tissues. There was a dramatic sex-specific effect of the clobetasol treatment; the peak abundance of B. burgdorferi in the skin (left hindleg, right hindleg, dorsal, ventral) was 72.6x higher in male mice compared to female mice. In contrast, there was little difference between the sexes in the tissue spirochete load in the ears, bladder, and heart. Topical application of clobetasol could increase the sensitivity of direct diagnostic methods (e.g., culture, PCR) to detect B. burgdorferi in host skin biopsies.

Seasonal dynamics of tick burden and associated Borrelia burgdorferi s.l. and Borrelia miyamotoi infections in rodents in a Dutch forest ecosystem

Ixodes ricinus ticks transmit Borrelia burgdorferi sensu lato (s.l.) as well as Borrelia miyamotoi. Larvae become infected when feeding on infected rodents, with horizontal transmission of B. burgdorferi and horizontal and vertical transmission of B. miyamotoi. We studied seasonal dynamics of infection rates of I. ricinus and their rodent hosts, and hence transmission risk of these two distinctly different Borrelia species. Rodents were live-trapped and inspected for ticks from May to November in 2013 and 2014 in a forest in The Netherlands. Trapped rodents were temporarily housed in the laboratory and detached ticks were collected. Borrelia infections were determined from the trapped rodents and collected ticks. Borrelia burgdorferi s.l. and B. miyamotoi were found in ticks as well as in rodents. Rodent density was higher in 2014, whereas tick burden as well as the Borrelia infection rates in rodents were higher in 2013. The density of B. miyamotoi-infected nymphs did not differ between the years. Tick burdens were higher on Apodemus sylvaticus than on Myodes glareolus, and higher on males than on females. Borrelia-infection rate of rodents varied strongly seasonally, peaking in summer. As the larval tick burden also peaked in summer, the generation of infected nymphs was highest in summer. We conclude that the heterogeneity of environmental and host-specific factors affects the seasonal transmission of Borrelia spp., and that these effects act more strongly on horizontally transmitted B. burgdorferi spp. than on the vertically transmitted B. miyamotoi.

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.

Borrelia burgdorferi (Spirochaetales: Spirochaetaceae) Infection Prevalence and Host Associations of Ticks Found on Peromyscus spp. in Maryland

Lyme disease, caused by Borrelia burgdorferi sensu stricto and most commonly transmitted by Ixodes scapularis Say (Ixodida: Ixodidae), is the most common tick-borne disease in Maryland. Because B. burgdorferi s.s. is maintained in enzootic cycles among wild mice (Peromyscus spp) and Ixodes spp ticks, differing patterns of parasitism of ticks on mice could impact the infection prevalence with B. burgdorferi. We determined the infection prevalence of Peromyscus spp as well as questing and partially engorged nymphal ticks collected at six sites on private land in five counties in Maryland from May to August 2020. Questing nymph infection prevalence (NIP) was 14%. We trapped 1258 mice and collected 554 ticks and 413 ear tissue samples. The prevalence of infested Peromyscus spp varied based on host age and sex, with older and male mice more likely to be infested. We detected a significant difference amongst the proportion of attached Ixodes and the location of trapping. Similarly, the prevalence of B. burgdorferi infected Peromyscus spp mice varied between locations (average mouse infection prevalence was 40%), with the highest prevalence in locations where Ixodes were the most commonly found ticks. The B. burgdorferi infection prevalence in partially engorged I. scapularis nymphs retrieved from Peromyscus spp was ~36% which lends further support to the host infection prevalence. Local differences in distribution of infected vectors and reservoirs are important factors to consider when planning interventions to reduce Lyme disease risk.

Context-dependent host dispersal and habitat fragmentation determine heterogeneity in infected tick burdens: an agent-based modelling study

As the incidence of tick-borne diseases has sharply increased over the past decade, with serious consequences for human and animal health, there is a need to identify ecological drivers contributing to heterogeneity in tick-borne disease risk. In particular, the relative importance of animal host dispersal behaviour in its three context-dependent phases of emigration, transfer and settlement is relatively unexplored. We built a spatially explicit agent-based model to investigate how the host dispersal process, in concert with the tick and host demographic processes, habitat fragmentation and the pathogen transmission process, affects infected tick distributions among hosts. A sensitivity analysis explored the impacts of different input parameters on infected tick burdens on hosts and infected tick distributions among hosts. Our simulations indicate that ecological predictors of infected tick burdens differed among the post-egg life stages of ticks, with tick attachment and detachment, tick questing activity and pathogen transmission dynamics identified as key processes, in a coherent way. We also found that the type of host settlement strategy and the proportion of habitat suitable for hosts determined super-spreading of infected ticks. We developed a theoretical mechanistic framework that can serve as a first step towards applied studies of on-the-ground public health intervention strategies.

Context-dependent host dispersal and habitat fragmentation determine heterogeneity in infected tick burdens: an agent-based modelling study

As the incidence of tick-borne diseases has sharply increased over the past decade, with serious consequences for human and animal health, there is a need to identify ecological drivers contributing to heterogeneity in tick-borne disease risk. In particular, the relative importance of animal host dispersal behaviour in its three context-dependent phases of emigration, transfer and settlement is relatively unexplored. We built a spatially explicit agent-based model to investigate how the host dispersal process, in concert with the tick and host demographic processes, habitat fragmentation and the pathogen transmission process, affects infected tick distributions among hosts. A sensitivity analysis explored the impacts of different input parameters on infected tick burdens on hosts and infected tick distributions among hosts. Our simulations indicate that ecological predictors of infected tick burdens differed among the post-egg life stages of ticks, with tick attachment and detachment, tick questing activity and pathogen transmission dynamics identified as key processes, in a coherent way. We also found that the type of host settlement strategy and the proportion of habitat suitable for hosts determined super-spreading of infected ticks. We developed a theoretical mechanistic framework that can serve as a first step towards applied studies of on-the-ground public health intervention strategies.

The role of host phenology for parasite transmission

Phenology is a fundamental determinant of species distributions, abundances, and interactions. In host–parasite interactions, host phenology can affect parasite fitness due to the temporal constraints it imposes on host contact rates. However, it remains unclear how parasite transmission is shaped by the wide range of phenological patterns observed in nature. We develop a mathematical model of the Lyme disease system to study the consequences of differential tick developmental-stage phenology for the transmission of B. burgdorferi. Incorporating seasonal tick activity can increase B. burgdorferi fitness compared to continuous tick activity but can also prevent transmission completely. B. burgdorferi fitness is greatest when the activity period of the infectious nymphal stage slightly precedes the larval activity period. Surprisingly, B. burgdorferi is eradicated if the larval activity period begins long after the end of nymphal activity due to a feedback with mouse population dynamics. These results highlight the importance of phenology, a common driver of species interactions, for the fitness of a parasite.

Host infection and community composition predict vector burden

Lyme disease is the most prevalent vector-borne disease in the United States, yet critical gaps remain in our understanding of tick and host interactions that shape disease dynamics. Rodents such as deer mice (Peromyscus spp.) and dusky-footed woodrats (Neotoma fuscipes) are key reservoirs for Borrelia burgdorferi, the etiological bacterium of Lyme disease, and can vary greatly in abundance between habitats. The aggregation of Ixodes pacificus, the western black-legged tick, on rodent hosts is often assumed to be constant across various habitats and not dependent on the rodent or predator communities; however, this is rarely tested. The factors that determine tick burdens on key reservoir hosts are important in estimating Lyme disease risk because larger tick burdens can amplify pathogen transmission. This study is the first to empirically measure I. pacificus larval burdens on competent reservoir hosts as a function of community factors such as rodent diversity, predator diversity, and questing tick abundance. Rodents were live trapped at oak woodland sites to collect tick burdens and tissue samples to test for infection with Borrelia burgdorferi sensu lato. We found that N. fuscipes tick burdens were negatively correlated with predator diversity, but positively correlated with questing I. pacificus larvae. In addition, rodent hosts that were infected with B. burgdorferi sensu lato tend to have higher burdens of larval ticks. These results demonstrate that tick burdens can be shaped by variability between individuals, species, and the broader host community with consequences for transmission and prevalence of tick-borne pathogens.

Dynamic rodent behavioral response to predation risk: implications for disease ecology

Prey modify their behavior in response to variation in predation risk, and such modifications can affect trophic processes such as disease transmission. However, variation in predation risk is complex, arising from direct risk from the predator itself and indirect risk due to the environment. Moreover, direct risk typically stems from multiple predators and varies over timescales (e.g., a predator nearby vs. its seasonal activities). We implemented a field-based experiment to disentangle these sources of risk and relate them to antipredator behavior in rodents. We modeled rodent occurrence and activity as a function of short- and long-term risk from a primary predator, red foxes (Vulpes vulpes), long-term risk from a second predator, coyotes (Canis latrans), and environmental variables. We found that long-term red fox activity strongly reduced rodent occurrence and that cues of nearby red fox presence decreased rodent activity by > 50%. In addition, this activity reduction was dynamic in that varied according to the background level of long-term red fox activity. Importantly, rodents did not respond to environmental variables (moonlight, temperature, and habitat) or long-term coyote activity. These results bear upon recent work that suggests predators can alter tick-borne disease dynamics via induced antipredator behavior of rodents, which are hosts for pathogens and ticks. Specifically, our study corroborates the hypothesis that red foxes act as important proximal agents in regulating tick-borne diseases by reducing rodent activity. More generally, this study highlights the need to consider the dynamic nature of prey antipredator response across landscapes with variable long-term predation risk.

Landscapes within landscapes: A parasite utilizes different ecological niches on the host landscapes of two host species

Parasites are distributed across populations of hosts, but also across microhabitats on or inside hosts: together the host population distribution and "host landscape" distribution comprise a part of the ecological niche of a parasite. In this paper we examine how a generalist parasite, the tick Ixodes holocyclus, is distributed at both the host population and host landscape scales in two species of host (Perameles nasuta and Rattus rattus) that co-occur. We anaesthetized wildlife to locate ticks from five generalized host body regions; we then analysed the distribution of ticks among the populations of hosts (aggregation) and the distribution of ticks among the available host body regions as niches. Ixodes holocyclus is more aggregated in the R. rattus population. At the host landscape scale, I. holocyclus's utilized niche includes the entire surface of P. nasuta equally, while the niche on R. rattus is focused on the head. Differences in tick aggregation between host species may reflect tick habitat suitability at the host landscape scale, as well as differences in ecological and evolutionary histories. By exploring the distribution of parasites at multiple scales, including host landscapes, we can better understand the complex ecology of parasites.

A 117-year retrospective analysis of Pennsylvania tick community dynamics

Background

Tick-borne diseases have been increasing at the local, national, and global levels. Researchers studying ticks and tick-borne diseases need a thorough knowledge of the pathogens, vectors, and epidemiology of disease spread. Both active and passive surveillance approaches are typically used to estimate tick population size and risk of tick encounter. Our data consists of a composite of active and long-term passive surveillance, which has provided insight into spatial variability and temporal dynamics of ectoparasite communities and identified rarer tick species. We present a retrospective analysis on compiled data of ticks from Pennsylvania over the last 117 years.

Methods

We compiled data from ticks collected during tick surveillance research, and from citizen-based submissions. The majority of the specimens were submitted by citizens. However, a subset of the data was collected through active methods (flagging or dragging, or removal of ticks from wildlife). We analyzed all data from 1900–2017 for tick community composition, host associations, and spatio-temporal dynamics.

Results

In total there were 4491 submission lots consisting of 7132 tick specimens. Twenty-four different species were identified, with the large proportion of submissions represented by five tick species. We observed a shift in tick community composition in which the dominant species of tick (Ixodes cookei) was overtaken in abundance by Dermacentor variabilis in the early 1990s and then replaced in abundance by I. scapularis. We analyzed host data and identified overlaps in host range amongst tick species.

Conclusions

We highlight the importance of long-term passive tick surveillance in investigating the ecology of both common and rare tick species. Information on the geographical distribution, host-association, and seasonality of the tick community can help researchers and health-officials to identify high-risk areas.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3451-6) contains supplementary material, which is available to authorized users.

Multiflora rose invasion amplifies prevalence of Lyme disease pathogen, but not necessarily Lyme disease risk

BACKGROUND

Forests in urban landscapes differ from their rural counterparts in ways that may alter vector-borne disease dynamics. In urban forest fragments, tick-borne pathogen prevalence is not well characterized; mitigating disease risk in densely-populated urban landscapes requires understanding ecological factors that affect pathogen prevalence. We trapped blacklegged tick (Ixodes scapularis) nymphs in urban forest fragments on the East Coast of the United States and used multiplex real-time PCR assays to quantify the prevalence of four zoonotic, tick-borne pathogens. We used Bayesian logistic regression and WAIC model selection to understand how vegetation, habitat, and landscape features of urban forests relate to the prevalence of B. burgdorferi (the causative agent of Lyme disease) among blacklegged ticks.

RESULTS

In the 258 nymphs tested, we detected Borrelia burgdorferi (11.2% of ticks), Borrelia miyamotoi (0.8%) and Anaplasma phagocytophilum (1.9%), but we did not find Babesia microti (0%). Ticks collected from forests invaded by non-native multiflora rose (Rosa multiflora) had greater B. burgdorferi infection rates (mean = 15.9%) than ticks collected from uninvaded forests (mean = 7.9%). Overall, B. burgdorferi prevalence among ticks was positively related to habitat features (e.g. coarse woody debris and total understory cover) favorable for competent reservoir host species.

CONCLUSIONS

Understory structure provided by non-native, invasive shrubs appears to aggregate ticks and reservoir hosts, increasing opportunities for pathogen transmission. However, when we consider pathogen prevalence among nymphs in context with relative abundance of questing nymphs, invasive plants do not necessarily increase disease risk. Although pathogen prevalence is greater among ticks in invaded forests, the probability of encountering an infected tick remains greater in uninvaded forests characterized by thick litter layers, sparse understories, and relatively greater questing tick abundance in urban landscapes.

The impact of strain-specific immunity on Lyme disease incidence is spatially heterogeneous

Lyme disease, caused by the bacterium Borrelia burgdorferi, is the most common tick-borne infection in the US. Recent studies have demonstrated that the incidence of human Lyme disease would have been even greater were it not for the presence of strain-specific immunity, which protects previously infected patients against subsequent infections by the same B. burgdorferi strain. Here, spatial heterogeneity is incorporated into epidemiological models to accurately estimate the impact of strain-specific immunity on human Lyme disease incidence. The estimated reduction in the number of Lyme disease cases is greater in epidemiologic models that explicitly include the spatial distribution of Lyme disease cases reported at the county level than those that utilize nationwide data. strain-specific immunity has the greatest epidemiologic impact in geographic areas with the highest Lyme disease incidence due to the greater proportion of people that have been previously infected and have developed strain-specific immunity.

A Borrelia afzelii Infection Increases Larval Tick Burden on Myodes glareolus (Rodentia: Cricetidae) and Nymphal Body Weight of Ixodes ricinus (Acari: Ixodidae)

Several microorganisms have been shown to manipulate their host or vector to enhance their own transmission. Here we examined whether an infection with Borrelia afzelii affects its transmission between its bank vole (Myodes glareolus, Schreber, 1780) host and tick vector. Captive-bred bank voles were inoculated with B. afzelii or control medium, after which host preference of Ixodes ricinus L. nymphs was determined in a Y-tube olfactometer. Thereafter, infected and uninfected bank voles were placed in a semifield arena containing questing larvae to measure larval tick attachment. Engorged larvae were collected from these bank voles, molted into nymphs, weighed, and analyzed for infection by PCR.Nymphs were attracted to the odors of a bank vole compared to ambient air and preferred the odors of an infected bank vole over that of an uninfected bank vole. In the semifield arena, infected male bank voles had greater larval tick burdens then uninfected males, while similar larval tick burdens were observed on females regardless of infection status. Nymphal ticks that acquired a B. afzelii infection had higher body weight than nymphs that did not acquire an infection regardless of the infection status of the vole. These results show that a B. afzelii infection in bank voles increases larval tick burden and that a B. afzelii infection in larvae increases nymphal body weight. This finding provides novel ecological insights into the enzootic cycle of B. afzelii.

Multistrain Infections with Lyme Borreliosis Pathogens in the Tick Vector

Mixed or multiple-strain infections are common in vector-borne diseases and have important implications for the epidemiology of these pathogens. Previous studies have mainly focused on interactions between pathogen strains in the vertebrate host, but little is known about what happens in the arthropod vector. Borrelia afzelii and Borrelia garinii are two species of spirochete bacteria that cause Lyme borreliosis in Europe and that share a tick vector, Ixodes ricinus Each of these two tick-borne pathogens consists of multiple strains that are often differentiated using the highly polymorphic ospC gene. For each Borrelia species, we studied the frequencies and abundances of the ospC strains in a wild population of I. ricinus ticks that had been sampled from the same field site over a period of 3 years. We used quantitative PCR (qPCR) and 454 sequencing to estimate the spirochete load and the strain diversity within each tick. For B. afzelii, there was a negative relationship between the two most common ospC strains, suggesting the presence of competitive interactions in the vertebrate host and possibly the tick vector. The flat relationship between total spirochete abundance and strain richness in the nymphal tick indicates that the mean abundance per strain decreases as the number of strains in the tick increases. Strains with the highest spirochete load in the nymphal tick were the most common strains in the tick population. The spirochete abundance in the nymphal tick appears to be an important life history trait that explains why some strains are more common than others in nature.

IMPORTANCE

Lyme borreliosis is the most common vector-borne disease in the Northern Hemisphere and is caused by spirochete bacteria that belong to the Borrelia burgdorferi sensu lato species complex. These tick-borne pathogens are transmitted among vertebrate hosts by hard ticks of the genus Ixodes Each Borrelia species can be further subdivided into genetically distinct strains. Multiple-strain infections are common in both the vertebrate host and the tick vector and can result in competitive interactions. To date, few studies on multiple-strain vector-borne pathogens have investigated patterns of cooccurrence and abundance in the arthropod vector. We demonstrate that the abundance of a given strain in the tick vector is negatively affected by the presence of coinfecting strains. In addition, our study suggests that the spirochete abundance in the tick is an important life history trait that can explain why some strains are more common in nature than others.

Multi-trophic interactions driving the transmission cycle of Borrelia afzelii between Ixodes ricinus and rodents: a review

The tick Ixodes ricinus is the main vector of the spirochaete Borrelia burgdorferi sensu lato, the causal agent of Lyme borreliosis, in the western Palearctic. Rodents are the reservoir host of B. afzelii, which can be transmitted to I. ricinus larvae during a blood meal. The infected engorged larvae moult into infected nymphs, which can transmit the spirochaetes to rodents and humans. Interestingly, even though only about 1 % of the larvae develop into a borreliae-infected nymph, the enzootic borreliae lifecycle can persist. The development from larva to infected nymph is a key aspect in this lifecycle, influencing the density of infected nymphs and thereby Lyme borreliosis risk. The density of infected nymphs varies temporally and geographically and is influenced by multi-trophic (tick-host-borreliae) interactions. For example, blood feeding success of ticks and spirochaete transmission success differ between rodent species and host-finding success appears to be affected by a B. afzelii infection in both the rodent and the tick. In this paper, we review the major interactions between I. ricinus, rodents and B. afzelii that influence this development, with the aim to elucidate the critical factors that determine the epidemiological risk of Lyme borreliosis. The effects of the tick, rodent and B. afzelii on larval host finding, larval blood feeding, spirochaete transmission from rodent to larva and development from larva to nymph are discussed. Nymphal host finding, nymphal blood feeding and spirochaete transmission from nymph to rodent are the final steps to complete the enzootic B. afzelii lifecycle and are included in the review. It is concluded that rodent density, rodent infection prevalence, and tick burden are the major factors affecting the development from larva to infected nymph and that these interact with each other. We suggest that the B. afzelii lifecycle is dependent on the aggregation of ticks among rodents, which is manipulated by the pathogen itself. Better understanding of the processes involved in the development and aggregation of ticks results in more precise estimates of the density of infected nymphs, and hence predictions of Lyme borreliosis risk.

Factors affecting larval tick feeding success: host, density and time

BACKGROUND

Ectoparasites rely on blood-feeding to sustain activity, support development and produce offspring. Blood-feeding is also a route for transmission of diverse vector-borne pathogens. The likelihood of successfully feeding is thus an important aspect of ectoparasite population dynamics and pathogen transmission. Factors that affect blood-feeding include ectoparasite density, host defenses, and ages of the host and ectoparasite. How these factors interact to affect feeding success is not well understood.

METHODS

We monitored blood-feeding success of larval Rocky Mountain wood ticks (RMWTs; Dermacentor andersoni) on deer mice (Peromyscus maniculatus) in several experiments to determine how tick density, host defense, and ages of mice and ticks interact to influence feeding success. In the first experiment, tick-naive deer mice were infested with one of several densities of RMWT larvae, while a second cohort of mice were infested with 50 larvae each. Two weeks after ticks dropped off, mice in the first cohort were re-exposed to 50 larvae each and mice in the second cohort were re-exposed to varying densities of larvae. In the second experiment mice of different ages (45-374 days old) were exposed to 50 larvae each. Two weeks later mice were re-exposed to 50 larvae each. We combined data from these and several similar experiments to test the generality of the patterns we observed. Lastly, we tested whether tick feeding success was consistent on individual mice that were challenged on four occasions.

RESULTS

Mice acquired resistance such that feeding success declined dramatically from the first to the second infestation. Feeding success also declined with tick density and tick age. Mice, however, became more permissive with age. The sizes of these effects were similar and additive. Surprisingly, over successive infestations the relative resistance among mice changed among hosts within a cohort.

CONCLUSIONS

We predict that larval blood-feeding success, and thus development to the nymph stage, will change due to variation in tick age and density, as well as the age and history of the host. Incorporating these biotic factors into modeling of tick population dynamics may improve predictions of tick-borne pathogen transmission.

First arrived takes all: inhibitory priority effects dominate competition between co-infecting Borrelia burgdorferi strains

Background

Within-host microbial communities and interactions among microbes are increasingly recognized as important factors influencing host health and pathogen transmission. The microbial community associated with a host is indeed influenced by a complex network of direct and indirect interactions between the host and the lineages of microbes it harbors, but the mechanisms are rarely established. We investigated the within-host interactions among strains of Borrelia burgdorferi, the causative agent of Lyme disease, using experimental infections in mice. We used a fully crossed-design with three distinct strains, each group of hosts receiving two sequential inoculations. We used data from these experimental infections to assess the effect of coinfection on bacterial dissemination and fitness (by measuring the transmission of bacteria to xenodiagnostic ticks) as well as the effect of coinfection on host immune response compared to single infection.

Results

The infection and transmission data strongly indicate a competitive interaction among B. burgdorferi strains within a host in which the order of appearance of the strain is the main determinant of the competitive outcome. This pattern is well described by the classic priority effect in the ecological literature. In all cases, the primary strain a mouse was infected with had an absolute fitness advantage primarily since it was transmitted an order of magnitude more than the secondary strain. The mechanism of exclusion of the secondary strain is an inhibition of the colonization of mouse tissues, even though 29% of mice showed some evidence of infection by secondary strain. Contrary to expectation, the strong and specific adaptive immune response evoked against the primary strain was not followed by production of immunoglobulins after the inoculation of the secondary strain, neither against strain-specific antigen nor against antigens common to all strains. Hence, the data do not support a major role of the immune response in the observed priority effect.

Conclusion

The strong inhibitory priority effect is a dominant mechanism underlying competition for transmission between coinfecting B. burgdorferi strains, most likely through resource exploitation. The observed priority effect could shape bacterial diversity in nature, with consequences in epidemiology and evolution of the disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-015-0381-0) contains supplementary material, which is available to authorized users.

The lyme disease pathogen has no effect on the survival of its rodent reservoir host

Zoonotic pathogens that cause devastating morbidity and mortality in humans may be relatively harmless in their natural reservoir hosts. The tick-borne bacterium Borrelia burgdorferi causes Lyme disease in humans but few studies have investigated whether this pathogen reduces the fitness of its reservoir hosts under natural conditions. We analyzed four years of capture-mark-recapture (CMR) data on a population of white-footed mice, Peromyscus leucopus, to test whether B. burgdorferi and its tick vector affect the survival of this important reservoir host. We used a multi-state CMR approach to model mouse survival and mouse infection rates as a function of a variety of ecologically relevant explanatory factors. We found no effect of B. burgdorferi infection or tick burden on the survival of P. leucopus. Our estimates of the probability of infection varied by an order of magnitude (0.051 to 0.535) and were consistent with our understanding of Lyme disease in the Northeastern United States. B. burgdorferi establishes a chronic avirulent infection in their rodent reservoir hosts because this pathogen depends on rodent mobility to achieve transmission to its sedentary tick vector. The estimates of B. burgdorferi infection risk will facilitate future theoretical studies on the epidemiology of Lyme disease.

Borrelia burgdorferi promotes the establishment of Babesia microti in the northeastern United States

Babesia microti and Borrelia burgdorferi, the respective causative agents of human babesiosis and Lyme disease, are maintained in their enzootic cycles by the blacklegged tick (Ixodes scapularis) and use the white-footed mouse (Peromyscus leucopus) as primary reservoir host. The geographic range of both pathogens has expanded in the United States, but the spread of babesiosis has lagged behind that of Lyme disease. Several studies have estimated the basic reproduction number (R₀) for B. microti to be below the threshold for persistence (<1), a finding that is inconsistent with the persistence and geographic expansion of this pathogen. We tested the hypothesis that host coinfection with B. burgdorferi increases the likelihood of B. microti transmission and establishment in new areas. We fed I. scapularis larva on P. leucopus mice that had been infected in the laboratory with B. microti and/or B. burgdorferi. We observed that coinfection in mice increases the frequency of B. microti infected ticks. To identify the ecological variables that would increase the probability of B. microti establishment in the field, we integrated our laboratory data with field data on tick burden and feeding activity in an R₀ model. Our model predicts that high prevalence of B. burgdorferi infected mice lowers the ecological threshold for B. microti establishment, especially at sites where larval burden on P. leucopus is lower and where larvae feed simultaneously or soon after nymphs infect mice, when most of the transmission enhancement due to coinfection occurs. Our studies suggest that B. burgdorferi contributes to the emergence and expansion of B. microti and provides a model to predict the ecological factors that are sufficient for emergence of B. microti in the wild.

Co-feeding transmission in Lyme disease pathogens

This review examines the phenomenon of co-feeding transmission in tick-borne pathogens. This mode of transmission is critical for the epidemiology of several tick-borne viruses but its importance for Borrelia burgdorferi sensu lato, the causative agents of Lyme borreliosis, is still controversial. The molecular mechanisms and ecological factors that facilitate co-feeding transmission are therefore examined with particular emphasis on Borrelia pathogens. Comparison of climate, tick ecology and experimental infection work suggests that co-feeding transmission is more important in European than North American systems of Lyme borreliosis, which potentially explains why this topic has gained more traction in the former continent than the latter. While new theory shows that co-feeding transmission makes a modest contribution to Borrelia fitness, recent experimental work has revealed new ecological contexts where natural selection might favour co-feeding transmission. In particular, co-feeding transmission might confer a fitness advantage in the Darwinian competition among strains in mixed infections. Future studies should investigate the ecological conditions that favour the evolution of this fascinating mode of transmission in tick-borne pathogens.

When is a parasite not a parasite? Effects of larval tick burdens on white-footed mouse survival

Many animal species can carry considerable burdens of ectoparasites: parasites living on the outside of a host's body. Ectoparasite infestation can decrease host survival, but the magnitude and even direction of survival effects can vary depending on the type of ectoparasite and the nature and duration of the association. When ectoparasites also serve as vectors of pathogens, the effects of ectoparasite infestation on host survival have the potential to alter disease dynamics by regulating host populations and stabilizing transmission. We quantified the impact of larval Ixodes scapularis tick burdens on both within-season and overwinter survival of white-footed mice (Peromyscus leucopus) using a hierarchical Bayesian capture-mark-recapture model. I. scapularis and P. leucopus are, respectively, vectors and competent reservoirs for the causative agents of Lyme disease, anaplasmosis, and babesiosis. Using a data set of 5587 individual mouse capture histories over sixteen years, we found little evidence for any effect of tick burdens on either within-season or overwinter mouse survival probabilities. In male mice, tick burdens were positively correlated with within-season survival probabilities. Mean maximum tick burdens were also positively correlated with population rates of change during the concurrent breeding season. The apparent indifference of mice to high tick burdens may contribute to their effectiveness as reservoir hosts for several human zoonotic pathogens.

Natural infestation of Hydrochoerus hydrochaeris by Amblyomma dubitatum ticks

Natural infestation of Amblyomma dubitatum in relation to individual specific attributes of Hydrochoerus hydrochaeris such as sex, body mass and body condition was analyzed. The anatomical distribution of A. dubitatum on H. hyrochaeris was also evaluated. Prevalence of adults and nymphs were significantly higher than prevalence of larvae. Non-significant differences in the infestation levels were found among host sex. Multiple regression analysis did not show any statistically significant association among the level of infestation with ticks and body mass and body condition of the host. All parasitic tick stages were collected in all five anatomical areas of the host, but they exhibited significant differences in feeding site preference. Factors associated to the host which determine the high levels of infestation with A. dubitatum could be assigned to a combination of population-level properties of the host as abundance, ubiquity and aggregation, rather than individual specific attributes related to body condition, body mass or sex.

A simple model for the establishment of tick-borne pathogens of Ixodes scapularis: a global sensitivity analysis of R0

The basic reproduction number of a pathogen, R0, determines whether a pathogen will spread (R0>1), when introduced into a fully susceptible population or fade out (R0<1), because infected hosts do not, on average, replace themselves. In this paper we develop a simple mechanistic model for the basic reproduction number for a group of tick-borne pathogens that wholly, or almost wholly, depend on horizontal transmission to and from vertebrate hosts. This group includes the causative agent of Lyme disease, Borrelia burgdorferi, and the causative agent of human babesiosis, Babesia microti, for which transmission between co-feeding ticks and vertical transmission from adult female ticks are both negligible. The model has only 19 parameters, all of which have a clear biological interpretation and can be estimated from laboratory or field data. The model takes into account the transmission efficiency from the vertebrate host as a function of the days since infection, in part because of the potential for this dynamic to interact with tick phenology, which is also included in the model. This sets the model apart from previous, similar models for R0 for tick-borne pathogens. We then define parameter ranges for the 19 parameters using estimates from the literature, as well as laboratory and field data, and perform a global sensitivity analysis of the model. This enables us to rank the importance of the parameters in terms of their contribution to the observed variation in R0. We conclude that the transmission efficiency from the vertebrate host to Ixodes scapularis ticks, the survival rate of Ixodes scapularis from fed larva to feeding nymph, and the fraction of nymphs finding a competent host, are the most influential factors for R0. This contrasts with other vector borne pathogens where it is usually the abundance of the vector or host, or the vector-to-host ratio, that determine conditions for emergence. These results are a step towards a better understanding of the geographical expansion of currently emerging horizontally transmitted tick-borne pathogens such as Babesia microti, as well as providing a firmer scientific basis for targeted use of acaricide or the application of wildlife vaccines that are currently in development.

Reductions in human Lyme disease risk due to the effects of oral vaccination on tick-to-mouse and mouse-to-tick transmission

Vaccinating wildlife is becoming an increasingly popular method to reduce human disease risks from pathogens such as Borrelia burgdorferi, the causative agent of Lyme disease. To successfully limit human disease risk, vaccines targeting the wildlife reservoirs of B. burgdorferi must be easily distributable and must effectively reduce pathogen transmission from infected animals, given that many animals in nature will be infected prior to vaccination. We assessed the efficacy of an easily distributable oral bait vaccine based on the immunogenic outer surface protein A (OspA) to protect uninfected mice from infection and to reduce transmission from previously infected white-footed mice, an important reservoir host of B. burgdorferi. Oral vaccination of white-footed mice effectively reduces transmission of B. burgdorferi at both critical stages of the Lyme disease transmission cycle. First, oral vaccination of uninfected white-footed mice elicits an immune response that protects mice from B. burgdorferi infection. Second, oral vaccination of previously infected mice significantly reduces the transmission of B. burgdorferi to feeding ticks despite a statistically nonsignificant immune response. We used the estimates of pathogen transmission to and from vaccinated and unvaccinated mice to model the efficacy of an oral vaccination campaign targeting wild white-footed mice. Projection models suggest that the effects of the vaccine on both critical stages of the transmission cycle of B. burgdorferi act synergistically in a positive feedback loop to reduce the nymphal infection prevalence, and thus human Lyme disease risk, well below what would be expected from either effect alone. This study suggests that oral immunization of wildlife with an OspA-based vaccine can be a promising long-term strategy to reduce human Lyme disease risk.

Geographical and environmental factors driving the increase in the Lyme disease vector Ixodes scapularis

The population densities of many organisms have changed dramatically in recent history. Increases in the population density of medically relevant organisms are of particular importance to public health as they are often correlated with the emergence of infectious diseases in human populations. Our aim is to delineate increases in density of a common disease vector in North America, the blacklegged tick, and to identify the environmental factors correlated with these population dynamics. Empirical data that capture the growth of a population are often necessary to identify environmental factors associated with these dynamics. We analyzed temporally- and spatially-structured field collected data in a geographical information systems framework to describe the population growth of blacklegged ticks (Ixodes scapularis) and to identify environmental and climatic factors correlated with these dynamics. The density of the ticks increased throughout the study's temporal and spatial ranges. Tick density increases were positively correlated with mild temperatures, low precipitation, low forest cover, and high urbanization. Importantly, models that accounted for these environmental factors accurately forecast future tick densities across the region. Tick density increased annually along the south-to-north gradient. These trends parallel the increases in human incidences of diseases commonly vectored by I. scapularis. For example, I. scapularis densities are correlated with human Lyme disease incidence, albeit in a non-linear manner that disappears at low tick densities, potentially indicating that a threshold tick density is needed to support epidemiologically-relevant levels of the Lyme disease bacterium. Our results demonstrate a connection between the biogeography of this species and public health.