Year-to-year variation in the density of Ixodes ricinus ticks and the prevalence of the rodent-associated human pathogens Borrelia afzelii and B. miyamotoi in different forest types

Tick abundance and infection with three zoonotic bacteria are heterogeneous in a Belgian peri-urban forest

Ixodes ricinus is a vector of several pathogens of public health interest. While forests are the primary habitat for I. ricinus, its abundance and infection prevalence are expected to vary within forest stands. This study assesses the spatio-temporal variations in tick abundance and infection prevalence with three pathogens in and around a peri-urban forest where human exposure is high. Ticks were sampled multiple times in 2016 and 2018 in multiple locations with a diversity of undergrowth, using the consecutive drags method. Three zoonotic pathogens were screened for, Borrelia burgdorferi s.l., Coxiella burnetii, and Francisella tularensis. The influence of season, type of site and micro-environmental factors on tick abundance were assessed with negative binomial generalized linear mixed-effects models. We collected 1642 nymphs and 181 adult ticks. Ticks were most abundant in the spring, in warmer temperatures, and where undergrowth was higher. Sites with vegetation unaffected by human presence had higher abundance of ticks. Forest undergrowth type and height were significant predictors of the level of tick abundance in a forest. The consecutive drags method is expected to provide more precise estimates of tick abundance, presumably through more varied contacts with foliage. Borrelia burgdorferi s.l. prevalence was estimated from pooled ticks at 5.33%, C. burnetii was detected in six pools and F. tularensis was not detected. Borrelia afzelii was the dominant B. burgdorferi genospecies. Tick abundance and B. burgdorferi s.l. infection prevalence were lower than other estimates in Belgian forests.

Recreational hazard: Vegetation and host habitat use correlate with changes in tick-borne disease hazard at infrastructure within forest stands

Studies on density and pathogen prevalence of Ixodes ricinus indicate that vegetation and local host community drive much of their variation between green spaces. Contrarily, micro-geographic variation is understudied, although its understanding could reduce disease risk. We studied the density of infectious nymphal Ixodes sp. ("DIN", proxy for disease hazard), density of questing nymphs ("DON") and nymphal infection prevalence ("NIP") near recreational forest infrastructure. Drag sampling within forest stands and at adjacent benches and trails was combined with vegetation surveys, camera trapping hosts and pathogen screening of ticks. We analysed Borrelia burgdorferi s.l. and its genospecies, with complementary analyses on Rickettsia sp., Anaplasma phagocytophilum, Neoehrlichia mikurensis and Borrelia miyamotoi. DIN was highest in forest interior and at trails enclosed by forest. Lower disease hazard was observed at benches and trails at forest edges. This infrastructure effect can be attributed to variation in vegetation characteristics and the habitat use of tick hosts, specifically roe deer, rodents and songbirds. DON is the main driver of DIN at micro-geographic scale and negatively affected by infrastructure and forest edges. A positive association with vegetation cover in understorey and canopy was observed, as were positive trends for local rodent and songbird abundance. NIP of different pathogens was affected by different drivers. Lower B. burgdorferi s.l. prevalence in the interior of forest stands, driven by its most prevalent genospecies B. afzelii, points towards higher density of uninfected hosts there. B. afzelii was positively associated with understorey containing tall species and with high canopy cover, whereas local bird community composition predicts B. garinii prevalence. A positive effect of songbird abundance and a negative effect of pigeons were observed. Our findings support amplification and inhibition mechanisms within forest stands and highlight that the effect of established drivers of DIN may differ based on the considered spatial scale.

Tick-borne diseases under the radar in the North Sea Region

The impact of tick-borne diseases caused by pathogens such as Anaplasma phagocytophilum, Neoehrlichia mikurensis, Borrelia miyamotoi, Rickettsia helvetica and Babesia species on public health is largely unknown. Data on the prevalence of these pathogens in Ixodes ricinus ticks from seven countries within the North Sea Region in Europe as well as the types and availability of diagnostic tests and the main clinical features of their corresponding diseases is reported and discussed. Raised awareness is needed to discover cases of these under-recognized types of tick-borne disease, which should provide valuable insights into these diseases and their clinical significance.

Biodiversity in the Lyme-light: ecological restoration and tick-borne diseases in Europe

Biodiversity loss and the emergence of zoonotic diseases are two major global challenges. An urgent question is how ecosystems and wildlife communities can be restored whilst minimizing the risk of zoonotic diseases carried by wildlife. Here, we evaluate how current ambitions to restore Europe's natural ecosystems may affect the hazard of diseases vectored by the tick Ixodes ricinus at different scales. We find that effects of restoration efforts on tick abundance are relatively straightforward but that the interacting effects of vertebrate diversity and abundance on pathogen transmission are insufficiently known. Long-term integrated surveillance of wildlife communities, ticks, and their pathogens is needed to understand their interactions and to prevent nature restoration from increasing tick-borne disease (TBD) hazard.

The influence of local environmental factors in southwestern Poland on the abundance of Ixodes ricinus and prevalence of infection with Borrelia burgdorferi s.l. and B. miyamotoi

Ticks are important ectoparasites and vectors of pathogens that cause disease in humans and animals. The natural habitat of Ixodes ricinus ticks is forests, which are convenient habitats to search for hosts, including reservoir hosts, and therefore can be an important habitat source of tick-borne pathogens. The aim of the study was to assess the usefulness of detailed forest habitat-type maps to estimate the tick-borne risk at a local scale (Lower Silesia, SW Poland). For the purposes of estimating tick abundance, we used the land cover maps available from the Forest Data Bank. For I. ricinus collection, nine sites located in three forest habitat types were chosen: broadleaf forest, mixed broadleaf and coniferous forest and coniferous forest. Ticks were collected once a month from April to June 2018 and 2019 using the standard flagging method. At each of the nine sites, ticks were collected in four plots, of 100 m² each. Tick abundance was analysed using general linear mixed models (GLMM). A total of 2196 (10.1/100 m²) ticks were collected, including 2093 Ixodes ricinus (95.3%; 9.6/100 m²), 46 Dermacentor reticulatus (2.1%; 0.2/100 m²) and 57 Haemaphysalis concinna (2.6%; 0.3/100 m²). Among the collected I. ricinus were 589 larvae (28.1%; 2.7/100 m²), 1261 nymphs (60.3%; 5.8/100 m²), 128 females (6.1%; 0.6/100 m²) and 115 males (5.5%; 0.5/100 m²). We found a highly significant effect of forest habitat type on the density of ticks for broadleaf forest (coefficient = 1.87267, p-value = 2.79e - 07). Additionally, a significant influence of air temperature and relative humidity on the abundance of ticks was observed. During spring, the peak activity of I. ricinus was recorded in May and June. For DNA amplification of Borrelia burgdorferi s.l., a nested PCR method was used. Out of 494 I. ricinus, 83 (16.8%) were positive for Borrelia spp. The RFLP method showed the occurrence of five species including four belonging to the B. burgdorferi s.l. complex: B. afzelii (30.1%), B. garinii (38.6%), B. valaisiana (2.4%) and B. lusitaniae (18.1%). Furthermore, B. miyamotoi (9.6%), a species belonging to bacteria that cause relapsing fever as well as co-infection of B. miyamotoi/B.lusitaniae (1.2%) were found. The differences in the infection level of Borrelia spp. between broadleaf forest and mixed broadleaf and coniferous forest were statistically significant.

Seasonal patterns and spatial variation of Borrelia burgdorferi (sensu lato) infections in Ixodes ricinus in the Netherlands

BACKGROUND

The incidence of Lyme borreliosis varies over time and space through as yet incompletely understood mechanisms. In Europe, Lyme borreliosis is caused by infection with a Borrelia burgdorferi (s.l.) genospecies, which is primarily transmitted by a bite of Ixodes ricinus nymphs. The aim of this study was to investigate the spatial and temporal variation in nymphal infection prevalence of B. burgdorferi (s.l.) (NIP), density of questing nymphs (DON) and the resulting density of infected nymphs (DIN).

METHODS

We investigated the infection rates in I. ricinus nymphs that were collected monthly between 2009 and 2016 in 12 locations in the Netherlands. Using generalized linear mixed models, we explored how the NIP, DON and DIN varied during the seasons, between years and between locations. We also determined the genospecies of the Borrelia infections and investigated whether the genospecies composition differed between locations.

RESULTS

The overall NIP was 14.7%. A seasonal pattern in infection prevalence was observed, with higher estimated prevalences in the summer than in the spring and autumn. This, combined with higher nymphal densities in summer, resulted in a pronounced summer peak in the estimated DIN. Over the 7.5-year study period, a significant decrease in infection prevalence was found, as well as a significant increase in nymphal density. These two effects appear to cancel each other out; the density of infected nymphs, which is the product of NIP × DON, showed no significant trend over years. Mean infection prevalence (NIP, averaged over all years and all months) varied considerably between locations, ranging from 5 to 26%. Borrelia genospecies composition differed between locations: in some locations almost all infections consisted of B. afzelii, whereas other locations had more diverse genospecies compositions.

CONCLUSION

In the Netherlands, the summer peak in DIN is a result of peaks in both NIP and DON. No significant trend in DIN was observed over the years of the study, and variations in DIN between locations were mostly a result of the variation in DON. There were considerable differences in acarological risk between areas in terms of infection prevalence and densities of ticks as well as in Borrelia genospecies composition.

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.

A Suspected Case of Lyme Borreliosis in a Dog from Belgium

A 5-year-old Jack Russell Spaniel was presented in December 2017 to his veterinarian in Belgium for sudden weakness, reluctance to move, and pain. Blood analysis showed no deviations and serum increased levels of B. burgdorferi s.l. antibodies were detected. The dog recovered a few days after the onset of doxycycline treatment. This case illustrates the possible relationship between tick-borne diseases and orthopedic problems.

Vertical transmission rates of Borrelia miyamotoi in Ixodes scapularis collected from white-tailed deer

Borrelia miyamotoi is a relapsing fever spirochete transmitted by ticks in the Ixodes ricinus complex. In the eastern United States, B. miyamotoi is transmitted by I. scapularis, which also vectors several other pathogens including B. burgdorferi sensu stricto. In contrast to Lyme borreliae, B. miyamotoi can be transmitted vertically from infected female ticks to their progeny. Therefore, in addition to nymphs and adults, larvae can vector B. miyamotoi to wildlife and human hosts. Two widely varying filial infection prevalence (FIP) estimates - 6% and 73% - have been reported previously from two vertically infected larval clutches; to our knowledge, no other estimates of FIP or transovarial transmission (TOT) rates for B. miyamotoi have been described in the literature. Thus, we investigated TOT and FIP of larval clutches derived from engorged females collected from hunter-harvested white-tailed deer in 2015 (n = 664) and 2016 (n = 599) from Maine, New Hampshire, Tennessee, and Wisconsin. After engorged females oviposited in the lab, they (n = 492) were tested for B. miyamotoi infection by PCR. Subsequently, from each clutch produced by an infected female, larval pools, as well as 100 individual eggs or larvae, were tested. The TOT rate of the 11 infected females was 90.9% (95% CI; 57.1-99.5%) and the mean FIP of the resulting larval clutches was 84.4% (95% CI; 68.1-100%). Even though the overall observed vertical transmission rate (the product of TOT and FIP; 76.7%, 95% CI; 44.6-93.3%) was high, additional horizontal transmission may be required for enzootic maintenance of B. miyamotoi based on the results of a previously published deterministic model. Further investigation of TOT and FIP variability and the underlying mechanisms, both in nature and the laboratory, will be needed to resolve this question. Meanwhile, studies quantifying the acarological risk of Borrelia miyamotoi disease need to consider not only nymphs and adults, but larval I. scapularis as well.

Role of mustelids in the life-cycle of ixodid ticks and transmission cycles of four tick-borne pathogens

BACKGROUND

Elucidating which wildlife species significantly contribute to the maintenance of Ixodes ricinus populations and the enzootic cycles of the pathogens they transmit is imperative in understanding the driving forces behind the emergence of tick-borne diseases. Here, we aimed to quantify the relative contribution of four mustelid species in the life-cycles of I. ricinus and Borrelia burgdorferi (sensu lato) in forested areas and to investigate their role in the transmission of other tick-borne pathogens. Road-killed badgers, pine martens, stone martens and polecats were collected in Belgium and the Netherlands. Their organs and feeding ticks were tested for the presence of tick-borne pathogens.

RESULTS

Ixodes hexagonus and I. ricinus were found on half of the screened animals (n = 637). Pine martens had the highest I. ricinus burden, whereas polecats had the highest I. hexagonus burden. We detected DNA from B. burgdorferi (s.l.) and Anaplasma phagocytophilum in organs of all four mustelid species (n = 789), and Neoehrlichia mikurensis DNA was detected in all species, except badgers. DNA from B. miyamotoi was not detected in any of the investigated mustelids. From the 15 larvae of I. ricinus feeding on pine martens (n = 44), only one was positive for B. miyamotoi DNA, and all tested negative for B. burgdorferi (s.l.), N. mikurensis and A. phagocytophilum. The two feeding larvae from the investigated polecats (n = 364) and stone martens (n = 39) were negative for all four pathogens. The infection rate of N. mikurensis was higher in feeding nymphs collected from mustelids compared to questing nymphs, but not for B. burgdorferi (s.l.), B. miyamotoi or A. phagocytophilum.

CONCLUSIONS

Although all stages of I. ricinus can be found on badgers, polecats, pine and stone martens, their relative contribution to the life-cycle of I. ricinus in forested areas is less than 1%. Consequently, the relative contribution of mustelids to the enzootic cycles of I. ricinus-borne pathogens is negligible, despite the presence of these pathogens in organs and feeding ticks. Interestingly, all four mustelid species carried all stages of I. hexagonus, potentially maintaining enzootic cycles of this tick species apart from the cycle involving hedgehogs as main host species.

Seroprevalence of Borrelia burgdorferi in Belgian forestry workers and associated risk factors

Background

As forest is the preferred environment for ticks, forestry workers are exposed to tick bites and tick-borne diseases. We assessed the seroprevalence of anti-Borrelia burgdorferi (Bb) antibodies and investigated, using an integrated landscape approach, the individual and environmental factors associated with the seroprevalence of Bb in Belgian forestry workers, a high-risk group in Belgium.

Methods

A group of 310 Belgian forest workers was examined to assess the seroprevalence of anti-Borrelia IgG antibodies. Using principal component analysis and binary logistic regression, the joint effects of individual characteristics and environmental characteristics were examined.

Results

Sixty-seven of the 310 workers were seropositive for Lyme disease (LD), leading to a seroprevalence of 21.6%. The seroprevalence was higher among forest workers visiting forests more frequently (P = 0.003) or who reported over 100 tick bites (P-value < 0.001). The intensity of tick bites and the use of protection measures against tick bites have a positive impact on LD seroprevalence while the quantity of shadow from trees at ground level had a negative one.

Conclusions

This study showed that forest workers are a population at risk for LD and, by extension, at risk for various tick-borne diseases. In addition to the role of the environment, our results also showed the importance of considering exposure when predicting the risk of infection by Bb.