Quantitative ecology of ticks as a basis for transmission models of tick-borne pathogens

Effect of rodent density on tick and tick-borne pathogen populations: consequences for infectious disease risk

BACKGROUND

Rodents are considered to contribute strongly to the risk of tick-borne diseases by feeding Ixodes ricinus larvae and by acting as amplifying hosts for pathogens. Here, we tested to what extent these two processes depend on rodent density, and for which pathogen species rodents synergistically contribute to the local disease risk, i.e. the density of infected nymphs (DIN).

METHODS

In a natural woodland, we manipulated rodent densities in plots of 2500 m² by either supplementing a critical food source (acorns) or by removing rodents during two years. Untreated plots were used as controls. Collected nymphs and rodent ear biopsies were tested for the presence of seven tick-borne microorganisms. Linear models were used to capture associations between rodents, nymphs, and pathogens.

RESULTS

Investigation of data from all plots, irrespective of the treatment, revealed a strong positive association between rodent density and nymphal density, nymphal infection prevalence (NIP) with Borrelia afzelii and Neoehrlichia mikurensis, and hence DIN's of these pathogens in the following year. The NIP, but not the DIN, of the bird-associated Borrelia garinii, decreased with increasing rodent density. The NIPs of Borrelia miyamotoi and Rickettsia helvetica were independent of rodent density, and increasing rodent density moderately increased the DINs. In addition, NIPs of Babesia microti and Spiroplasma ixodetis decreased with increasing rodent density, which had a non-linear association with DINs of these microorganisms.

CONCLUSIONS

A positive density dependence for all rodent- and tick-associated tick-borne pathogens was found, despite the observation that some of them decreased in prevalence. The effects on the DINs were variable among microorganisms, more than likely due to contrasts in their biology (including transmission modes, host specificity and transmission efficiency). The strongest associations were found in rodent-associated pathogens that most heavily rely on horizontal transmission. Our results draw attention to the importance of considering transmission mode of a pathogen while developing preventative measures to successfully reduce the burden of disease.

Impact of vertebrate communities on Ixodes ricinus-borne disease risk in forest areas

Background

The density of questing ticks infected with tick-borne pathogens is an important parameter that determines tick-borne disease risk. An important factor determining this density is the availability of different wildlife species as hosts for ticks and their pathogens. Here, we investigated how wildlife communities contribute to tick-borne disease risk. The density of Ixodes ricinus nymphs infected with Borrelia burgdorferi (sensu lato), Borrelia miyamotoi, Neoehrlichia mikurensis and Anaplasma phagocytophilum among 19 forest sites were correlated to the encounter probability of different vertebrate hosts, determined by encounter rates as measured by (camera) trapping and mathematical modeling.

Result

We found that the density of any tick life stage was proportional to the encounter probability of ungulates. Moreover, the density of nymphs decreased with the encounter probability of hare, rabbit and red fox. The density of nymphs infected with the transovarially-transmitted B. miyamotoi increased with the density of questing nymphs and the encounter probability of bank vole. The density of nymphs infected with all other pathogens increased with the encounter probability of competent hosts: bank vole for Borrelia afzelii and N. mikurensis, ungulates for A. phagocytophilum and blackbird for Borrelia garinii and Borrelia valaisiana. The negative relationship we found was a decrease in the density of nymphs infected with B. garinii and B. valaisiana with the encounter probability of wood mouse.

Conclusions

Only a few animal species drive the densities of infected nymphs in forested areas. There, foxes and leporids have negative effects on tick abundance, and consequently on the density of infected nymphs. The abundance of competent hosts generally drives the abundances of their tick-borne pathogen. A dilution effect was only observed for bird-associated Lyme spirochetes.

Research on the ecology of ticks and tick-borne pathogens--methodological principles and caveats

Interest in tick-transmitted pathogens has experienced an upsurge in the past few decades. Routine application of tools for the detection of fragments of foreign DNA in ticks, together with a high degree of interest in the quantification of disease risk for humans, has led to a marked increase in the number of reports on the eco-epidemiology of tick-borne diseases. However, procedural errors continue to accumulate in the scientific literature, resulting in misleading information. For example, unreliable identification of ticks and pathogens, erroneous interpretations of short-term field studies, and the hasty acceptance of some tick species as vectors have led to ambiguities regarding the vector role of these arthropods. In this review, we focus on the ecological features driving the life cycle of ticks and the resulting effects on the eco-epidemiology of tick-transmitted pathogens. We review the factors affecting field collections of ticks, and we describe the biologically and ecologically appropriate procedures for describing tick host-seeking activity and its correlation with environmental traits. We detail the climatic variables that have biological importance on ticks and explain how they should be properly measured and analyzed. We also provide evidence to critically reject the use of some environmental traits that are being increasingly reported as the drivers of the behavior of ticks. With the aim of standardization, we propose unambiguous definitions of the status of hosts and ticks regarding their ability to maintain and spread a given pathogen. We also describe laboratory procedures and standards for evaluating the vectorial capacity of a tick or the reservoir role of a host. This approach should provide a coherent framework for the reporting of research findings concerning ticks and tick-borne diseases.

The paradox of the binomial Ixodes ricinus activity and the observed unimodal Lyme borreliosis season in Hungary

The change of ambient temperature plays a key role in determining the run of the annual Lyme season. Our aim was to explain the apparent contradiction between the annual unimodal Lyme borreliosis incidence and the bimodal Ixodes ricinus tick activity run--both observed in Hungary--by distinguishing the temperature-dependent seasonal human and tick activity, the temperature-independent factors, and the multiplicative effect of human outdoor activity in summer holiday, using data from Hungary in the period of 1998-2012. This separation was verified by modeling the Lyme incidence based on the separated factors, and comparing the run of the observed and modeled incidence. We demonstrated the bimodality of tick season by using the originally unimodal Lyme incidence data. To model the outdoor human activity, the amount of camping guest nights was used, which showed an irregular run from mid-June to September. The human outdoor activity showed a similar exponential correlation with ambient temperature to that what the relative incidence did. It was proved that summer holiday has great influence on Lyme incidence.

Seasonal activity and tick-borne pathogen infection rates of Ixodes ricinus ticks in Hungary

Ixodes ricinus is the most important tick species in Europe as it is most widely distributed and transmits the majority of tick-borne zoonotic pathogens. As limited data are available for Hungary, the aim of the present study was to investigate the seasonal timing of questing by I. ricinus and the infection rate of this tick species with all major tick-borne zoonotic pathogens. Monthly collections of I. ricinus were carried out over 3 consecutive years by dragging a blanket in 6 biotopes representing different areas of Hungary. Altogether, 1800 nymphs (300 per collection point) were screened as pooled samples (each of 5 specimens) by PCR-based methods for tick-borne pathogens. I. ricinus larvae, nymphs, and adults had bimodal activity patterns with a major peak in the spring. As newly moulted ticks of all stages are thought to emerge in the autumn of each year, it appears that most newly emerged ticks delayed their questing until the following spring. The minimum prevalence of Borrelia burgdorferi sensu lato was 2.5%. Borr. afzelii, Borr. burgdorferi sensu stricto, Borr. garinii, Borr. lusitaniae, and Borr. valaisiana were identified by hybridization. The minimum infection rate with spotted fever group rickettsiae was 1.9%. Rickettsia helvetica was identified in all biotopes. The minimum prevalence of Anaplasma phagocytophilum, Babesia divergens and Bab. microti was low (0.3-0.5%). Bartonella spp.-, Francisella tularensis-, and TBE virus-specific amplification products were not detected. Relative to the results of comparable studies carried out in the Carpathian Basin, the prevalence of tick-borne pathogens was low in Hungary. This might be attributed to the climatic difference between the lowland areas of Hungary and submountain areas of the surrounding countries involved in the studies.

Transmission dynamics of Borrelia lusitaniae and Borrelia afzelii among Ixodes ricinus, lizards, and mice in Tuscany, central Italy

To estimate the basic reproduction number (R(0)) of Borrelia lusitaniae and Borrelia afzelii, we formulated a mathematical model considering the interactions among the tick vector, vertebrate hosts, and pathogens in a 500-ha enclosed natural reserve on Le Cerbaie hills, Tuscany, central Italy. In the study area, Ixodes ricinus were abundant and were found infected by B. lusitaniae and B. afzelii. Lizards (Podarcis spp.) and mice (Apodemus spp.), respectively, are the reservoir hosts of these two Borrelia burgdorferi sensu lato (s.l.) genospecies and compete for immature ticks. B. lusitaniae R(0) estimation is in agreement with field observations, indicating the maintenance and diffusion of this genospecies in the study area, where lizards are abundant and highly infested by I. ricinus immature stages. In fact, B. lusitaniae shows a focal distribution in areas where the tick vector and the vertebrate reservoir coexist. Mouse population dynamics and their relatively low suitability as hosts for nymphs seem to determine, on the other hand, a less efficient transmission of B. afzelii, whose R(0) differs between scenarios in the study area. Considering host population dynamics, the proposed model suggests that, given a certain combination of the two host population sizes, both spirochete genospecies can coexist in our study area. Additional incompetent hosts for B. burgdorferi s.l. have a negative effect on B. afzelii maintenance, whose R(0) results > 1 only with high mouse population densities and/or low lizards abundance, but they do not seem to influence B. lusitaniae transmission cycle on Le Cerbaie. Secondly, our model confirms the importance of nymphs' infestation, of host population density and diversity, and spirochetes host association for the maintenance of the transmission cycle of B. burgdorferi s.l.

Influence of climatic factors on dynamics of questing Ixodes ricinus ticks in Slovenia

Ixodes ricinus is a vector of pathogens that cause many diseases in Europe and Slovenia: tick-borne encephalitis (TBE), anaplasmosis, borreliosis, babesiosis and others. The risk for contracting these diseases depends strongly on the density of the infected questing ticks and many studies have investigated tick population dynamics and the parameters affecting them. They have shown a clear influence of climatic and landscape arrangements in the microhabitat on tick abundance and dynamics and therefore on transmission of pathogens important in human and veterinary medicine. In our study we assessed the influence of climatic factors on questing activity of ticks over a three-year period at 7 locations in Slovenia. Locations were selected in endemic foci of TBE with different intensity, which were identified according to the presence of human disease. Sites differ according to various abiotic and biotic factors, such as climate, amount of rain, height above sea level, vegetation and wildlife. All three stages of ticks were collected monthly over a three-year period (2005-2007). Temperature, humidity and precipitation data were collected for these years. The purpose of our study was to relate observed differences in I. ricinus ticks questing activity to local climate. We found a correlation between the decrease of questing ticks in the summer and the combination of air temperatures and humidity in the form of saturation deficit.

Temporal variations in the usefulness of normalized difference vegetation index as a predictor for Ixodes ricinus (Acari: Ixodidae) in a Borrelia lusitaniae focus in Tuscany, central Italy

Host-seeking ticks were collected during monthly dragging sessions from November 2004 through October 2006 in Tuscany, central Italy. Normalized difference vegetation index (NDVI), which was calculated from Landsat ETM+ 7 remote sensing data recorded in August 2001, was significantly correlated with numbers of host-seeking immature Ixodes ricinus L. (Acari: Ixodidae) during periods of relatively low rainfall such as summer 2005 (Spearman's p = 0.78, P < 0.001 for nymphs in July) and to a lower extent in spring-summer 2006. In spring 2005, when rainfall was relatively high, the correlation was weak and not statistically significant. Generalized estimating equations (GEEs), taking into account repeated sampling of the same dragging sites, were used to model the effects of NDVI and season on counts of host-seeking I. ricinus nymphs. Seasonal variations of the effect of NDVI yielded a significant NDVI-by-season interaction in the first year of the study (November 2004-October 2005), but not in the second year (November 2005-October 2006) when there was a 2.5-fold increase of the number of nymphs per 100-m dragging for every 0.1 unit increase in NDVI (95% confidence interval = 1.6, 3.0). Risk maps that were obtained based on GEE results confirmed that the predicted number of I. ricinus nymphs per 100 m was relatively homogeneous through the study area during the 2005 spring peak of activity. Conversely, in 2006, the predicted abundance of nymphs was greater in moist bottomland habitat (characterized by high NDVI) than in dry, typically Mediterranean, upland habitat.

Temporal variations in the usefulness of normalized difference vegetation index as a predictor for Ixodes ricinus (Acari: Ixodidae) in a Borrelia lusitaniae focus in Tuscany, central Italy

Host-seeking ticks were collected during monthly dragging sessions from November 2004 through October 2006 in Tuscany, central Italy. Normalized difference vegetation index (NDVI), which was calculated from Landsat ETM+ 7 remote sensing data recorded in August 2001, was significantly correlated with numbers of host-seeking immature Ixodes ricinus L. (Acari: Ixodidae) during periods of relatively low rainfall such as summer 2005 (Spearman's p = 0.78, P < 0.001 for nymphs in July) and to a lower extent in spring-summer 2006. In spring 2005, when rainfall was relatively high, the correlation was weak and not statistically significant. Generalized estimating equations (GEEs), taking into account repeated sampling of the same dragging sites, were used to model the effects of NDVI and season on counts of host-seeking I. ricinus nymphs. Seasonal variations of the effect of NDVI yielded a significant NDVI-by-season interaction in the first year of the study (November 2004-October 2005), but not in the second year (November 2005-October 2006) when there was a 2.5-fold increase of the number of nymphs per 100-m dragging for every 0.1 unit increase in NDVI (95% confidence interval = 1.6, 3.0). Risk maps that were obtained based on GEE results confirmed that the predicted number of I. ricinus nymphs per 100 m was relatively homogeneous through the study area during the 2005 spring peak of activity. Conversely, in 2006, the predicted abundance of nymphs was greater in moist bottomland habitat (characterized by high NDVI) than in dry, typically Mediterranean, upland habitat.

The tick Ixodes ricinus: distribution and climate preferences in the western Palaearctic

In this study, multivariate spatial clustering on monthly normalized difference vegetation index (NDVI) maps is used to classify ecological regions over the western Palaearctic. This classification is then used to delineate the distribution and climate preferences of populations (clades) of the tick Ixodes ricinus L. (Acari: Ixodidae) from a geographically extensive dataset of tick records and a gridded 2.5-km resolution climate dataset. Using monthly layers of the NDVI, regions of similar ecological attributes were defined and nine populations with significant differences in critical climate parameters (P< 0.005) were detected. Grouping of tick records according to other categories, such as political divisions, a 4 degrees x 4 degrees grid overlying the study area, or the CORINE) and USGS) vegetation classification schemes did not provided significantly separated populations (P = 0.094-0.304). Factor analysis and hierarchical tree clustering provided an ecological overview of these tick clades: two Mediterranean and one Scandinavian (western) clades are clearly separated from a node that includes clades of different parts of central Europe and the British Isles, with contrasting affinities between the different clades. The capture records of these ecologically separated clades produce a clear bias when bioclimate envelope modelling is applied to the mapping of habitat suitability for the tick in the western Palaearctic. The best-performing methods (Cohen's kappa = 0.834-0.912) use partial models developed with data from each ecoregion, which are then overlapped over the region of study. It is concluded that the use of ecologically derived ecoregions is an objective step in assessing the presence of ecologically different clades, and provides a guide in the development of data partitioning for habitat suitability modelling.