Recreational hazard: Vegetation and host habitat use correlate with changes in tick-borne disease hazard at...

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.