Skip to main content

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



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.


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.


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.


Urbanization affects many aspects of vector-borne disease ecology [1]. In the case of tick-borne disease systems such as Lyme disease (caused by Borrelia burgdorferi) in forested ecosystems, urbanization alters habitat suitability for vectors (i.e. ticks), vertebrate hosts, and as a result, pathogens. Human development in the Lyme disease endemic, mid-Atlantic region of the United States reduces overall forest cover and average patch size while increasing the area of edge and impervious surface. Reduced forest patch size, in particular, results in predictable changes to host community composition that increase acarological risk in terms of nymphal infection prevalence and density of infected nymphs [2,3,4,5]. Yet in human-dominated landscapes, patch size may have a smaller or perhaps unpredictable influence on host community relative to other effects of urbanization on forested ecosystems. How ecological characteristics of urban forest fragments affect acarological risk has not been well explored.

Complex land use histories in human-dominated landscapes form networks of diverse, heterogeneous forest fragments. In the urban mid-Atlantic region, clear-cutting, intensive agriculture, and urban sprawl have created a variety of forest fragment types on a spectrum between remnants of mature (> 100 yr. old) forests and forest fragments that have regrown from fallow agricultural land set aside while surrounding areas were developed [6]. In the latter case, native tree species have competed with and grown alongside non-native species that were part of the agricultural landscape or subsequent development. As a result, regrown urban forest patches have closed canopies of mostly native trees with thick understories composed of non-native, invasive species [7]. These two extremes of urban forest fragment types both face serious ecological problems (e.g. loss of native understory or reduced regeneration), with implications for tick-borne disease risk.

Due to changes to below-ground processes and browsing pressure from high density white-tailed deer (Odocoileus virginianus) populations (in Delaware, recent county surveys estimate between 18 and 52 deer/km2 [8]), mature forests may have sparse or no woody understories and cannot replace many species of dead or dying trees [9]. Although they maintain a thick litter layer and low soil pH, which may help buffer mature forests from invasion by non-native plants [10], many native woody plants cannot regenerate. The thick litter layer maintained in these forests provides suitable habitat for blacklegged ticks (Ixodes scapularis), which are found in greater abundance in mature forests relative to other urbanized forest fragment types [11]. In contrast, forest fragments with significant non-native plant invasion in the understory have high densities of invasive earthworms and very little leaf litter [12,13,14,15], which constrains tick abundance [11, 16, 17]. However, the dense understory structure provided by invasive plants may aggregate immature ticks and infective hosts, potentially amplifying acarological risk in invaded forest fragments [18,19,20,21].

Recent studies have identified greater pathogen prevalence in ticks and reservoir hosts associated with invasive shrubs [18,19,20, 22]. However, because leaf litter loss, which constrains tick abundance, is also associated with non-native plant invasion, it is unclear how tick-borne disease risk differs in regrown, invaded forest fragments compared to mature, uninvaded fragments. To contribute to our understanding of tick-borne disease ecology in urbanized landscapes, we designed a study in urban forest fragments with three objectives: (i) characterize B. burgdorferi and emerging tick-borne pathogen prevalence among questing ticks; (ii) test for differences in pathogen prevalence between forests invaded by non-native understory plants and uninvaded forests; and (iii) determine which habitat and landscape features influence pathogen prevalence.


Study area and tick collection

We collected nymphal Ixodes scapularis ticks from April to July, 2013 and 2014, using CO2-baited traps in forest fragments around New Castle County, Delaware. Although drag-sampling or flagging is more commonly used to capture I. scapularis, we used CO2-baited traps to avoid confounding results from Rosa multiflora’s dense structure (Fig. 1). Even thick, canvas cloth becomes snagged on R. multiflora thorns, preventing effective sampling of tick habitat. CO2-baited traps are unbiased by habitat structure [23]. We built traps following the design of Kensinger & Allan [23], by drilling four holes in 6-quart Coleman® coolers and bolting the coolers to plywood squares. We baited traps for 24 h with 1.4 kg of pelleted dry ice and lined the plywood base with doubled-sided carpet tape (3 M, Maplewood, USA).

Fig. 1
figure 1

Photographs of representative forest sites with R. multiflora invasion (a) and without R. multiflora (b)

Forest fragments (6–16 ha) consisted of mixed deciduous hardwood stands and varied in understory woody species composition, particularly in the extent of non-native R. multiflora invasion (Fig. 1, Additional file 1: Figure S1). Each year we trapped ticks in eight forest fragments, four of which had understories with 10–58% of total area covered by R. multiflora invasion (hereafter: invaded), and four fragments lacked R. multiflora invasion (< 1%), (hereafter: uninvaded). Within invaded sites, we captured ticks at four sets of paired traps: one trap within R. multiflora cover and its pair 25 m away, not in R. multiflora. Paired traps were separated by 25 m to eliminate the possibility that both traps could be attracting the same ticks [24]. In uninvaded sites, we deployed traps at four random points, separated by at least 25 m. We used a total of 64 trap locations over the 2 yr. study, and half of the traps were active on any given trap night. To avoid weather-related impacts on tick questing behavior, we always deployed paired traps together, and baited in equal numbers of invaded and uninvaded fragments on the same nights. We transported all captured ticks to the laboratory live, in individual microcentrifuge tubes, froze them at -80 °C, and later identified them to species and life stage with dichotomous keys [25,26,27].

Covariate data collection

We identified a set of 25 variables that we expected would influence B. burgdorferi infection rates by increasing or decreasing interactions between larval ticks and competent reservoir hosts (Tables 1 and 2). Further detail on field and computational methods used to collect these data were published in [11, 28]. We surveyed understory vegetation characteristics within a 12.5 m radius surrounding each trap location, which included: estimating the percent of ground covered by R. multiflora, the percent of ground covered by coarse woody debris, and the density of understory vegetation using a 2.0 m high Nudds board for which observers estimated the percentage of each of four 0.5 m panels obscured by vegetation from a distance of 12.5 m [29]. We chose a 12.5 m radius to correspond with the approximate home range size of Peromyscus leucopus (an important reservoir for B. burgdorferi), while avoiding overlap with paired traps [30]. We also estimated the percent of ground covered by R. multiflora within a 2.5 m radius of the trap to more directly represent the effective trapping radius for I. scapularis [24]. We quantified leaf litter volume for all litter collected within a 0.5 m2 quadrat next to each trap.

Table 1 Summary of vegetation and landscape covariates measured at the trap scale by forest type and location, modified from [11]. Covariates are summarized as mean ± standard error. Different superscript letters A, B, C denote significant differences among groups (P < 0.05) detected using analysis of variance (ANOVA), blocking on site, followed up with Tukey’s post-hoc comparisons when there were more than two groups
Table 2 Summary of vegetation and landscape covariates measured at the patch scale by forest type (invaded or uninvaded), modified from [11]. Covariates are summarized as mean ± standard error. Superscript letters A, B denote significant differences among groups (P < 0.05) detected using analysis of variance (ANOVA), blocking on site

We measured landscape variables at each trap location in ArcGIS using a 2007 Delaware land use land cover layer [31], focusing on variables that could influence habitat suitability for ticks and/or hosts and that reflected the human-dominated landscape context of the study area [32, 33]: distance to nearest road, stream, agriculture, forest edge and residential development. We also used data from prior [11, 28] and concurrent studies (Adalsteinsson et al., unpublished data) to quantify abundance of ticks, potential hosts, and host-tick interactions in the study area. Tick abundance at the trap-level was the number of I. scapularis nymphs captured at a given trap, standardized by effort (number of trap nights). The densities of ground-foraging bird territories in forest fragments were estimated from spot-mapping surveys conducted during two breeding seasons [28]. Concurrent studies of P. leucopus abundance and parasitism by immature ticks (Adalsteinsson et al., unpublished data) provided estimates of mouse abundance and parasitism rates at the trap and forest fragment scale. To study P. lecuopus abundance, we checked 15 nest boxes per forest fragment once per month; for trap-level estimates, mouse abundance was the mean of the number of mice caught during fall (larval tick season) at the two nest boxes nearest to the trapping location. For patch-level estimates, the number of mice caught at nest boxes in fall was averaged across all 15 nest boxes in a given forest fragment. Larval tick burdens were the mean number of larvae per mouse at either the two closest nest boxes (trap-level) or across all 15 nest boxes (patch-level).

We also included data collected previously to characterize vegetation at the patch level: proportions of Fagus grandifolia, Acer spp., Quercus spp., Liriodendron tulipifera, or Liquidambar styraciflua as dominant canopy trees; percent of total area covered by R. multiflora; mean leaf litter volume measured at 15 locations in the patch; percent of ground covered by understory plants (all spp.); percentage of understory woody stems that were non-native; and year of canopy closure [28].

Pathogen testing

We used a modified version of the DNeasy Blood & Tissue Kit (Qiagen, Venlo, Netherlands) protocol to extract DNA from ticks. Here, we explain the steps in which we deviated from the manufacturer’s protocol. First, we used sterile pipette tips to manually crush each I. scapularis nymph individually in 20 μl of Hyclone Dulbecco’s phosphate buffer saline solution (Thermo Fisher Scientific, Waltham, USA). Next, we incubated samples with lysis buffer ATL and proteinase K in a 56 °C hot water bath for 3 h. We performed an extra spin step at 13,000× rpm to remove trace ethanol after the Buffer AW2 wash. Finally, we modified the last step by eluting our samples twice (50 μl each time), for a final product of 100 μl. We checked concentrations of a subset of our samples using a NanoDrop UV-Vis spectrophotometer (Thermo Fisher Scientific, Waltham, USA) to confirm successful DNA extractions.

We tested ticks for the presence of Borrelia burgdorferi (sensu lato), Anaplasma phagocytophilum, and Babesia microti using a previously described multiplex PCR assay [34]. In addition, ticks were also tested for the presence of Borrelia miyamotoi in a TaqMan PCR assay using the following primers and probe: F770-5′-ACC TGC AAC CTT CGG ATT C-3′; R771-5′-TGG TTG TAG CTC AGT TGG TAG-3′; P1277-CalRd610-5′-CTT GTA TCG AAC TAC ACC CAT AGC TC-3′-BHQ2.

Data analysis

A sufficient number of ticks tested positive for Borrelia burgdorferi to allow statistical analyses; however, infection prevalence was too low for the remaining pathogens to determine patterns related to invasion and other habitat and landscape features. We tested for spatial autocorrelation in B. burgdorferi prevalence across forest fragments using a spline correlogram in package ncf [35] in R [36].

Sample sizes were uneven because of differences in tick abundance in invaded and uninvaded sites [11], so we used logistic regression models in a Bayesian framework to compare pathogen (i.e. B. burgdorferi) prevalence among ticks collected from invaded and uninvaded forest fragments. This approach allowed us to include uncertainty in pathogen prevalence due to varying sample sizes rather than comparing raw proportions. We first tested whether pathogen prevalence differed among treatments (in rose vs not in rose) within invaded forest fragments. Next, we tested whether pathogen prevalence differed between invaded and uninvaded forests. For both of these questions, we used the following model structure:

$$ {\displaystyle \begin{array}{l}{Y}_i\sim \mathrm{Binomial}\left(1,{p}_i\right),\\ {}\mathrm{logit}\left({p}_i\right)=\alpha +{\beta}^{\ast }{X}_i,\end{array}} $$

where Y is the infection status of each individual (i) tick (0 or 1), which has a binomial distribution with probability p. With the logit link, the probability of infection (p) was estimated as a linear function of baseline infection probability (α) plus the effect of being from one of two treatment groups (β*X), where the variable X indicates whether the sampling site was either in rose or not in rose, or in an invaded forest or not (0 or 1), depending on the question. We specified vague priors for the parameters α and β including a normal distribution, and we used a Hamiltonian Monte Carlo sampler to run the analysis with STAN [37] using packages rstan [38] and rethinking [39] in R [36]. Using the difference between posterior distributions (89% highest posterior density intervals) of the estimated pathogen prevalence in two groups (in rose or not in rose, and invaded or uninvaded forest), we calculated the probability of infection prevalence in one treatment group being larger than the other, given our data and the specified model [40].

To understand the influence of habitat and landscape factors on B. burgdorferi prevalence, we used aggregated logistic regression models in a Bayesian framework [40]:

$$ {\displaystyle \begin{array}{l}{Y}_i\sim \mathrm{Binomial}\left(n,{p}_i\right),\\ {}\mathrm{logit}\left({p}_i\right)=\alpha +{\beta_1}^{\ast }{X}_{1i}+{\beta_2}^{\ast }{X}_{2i}+\dots +{\beta_z}^{\ast }{X}_{zi}.\end{array}} $$

Counts of B. burgdorferi-positive ticks (Y i ) were aggregated by trap (i), and the number of trials (n) was the number of ticks tested from each trap. We developed multivariate model sets according to variable types (vegetation, habitat, or landscape) and spatial scale (trap or forest fragment). To select between models, we used Watanabe-Akaike Information Criterion (WAIC) [40,41,42] and converged on several “best models” that included similar combinations of variables (Table 3). We compared the relative influence of individual variables on the performance of the model by systematically dropping variables and comparing changes in WAIC scores and their standard errors [43].

Table 3 WAIC table of best models. Variables in the model set are coarse woody debris (CWD), leaf litter volume (litter), distance to the nearest road (dist. Road), mouse abundance in fall (mice), total understory cover (total cover), and nymphal tick capture rate (tick abundance). Field headings refer to the effective number of parameters (pWAIC), the difference between WAIC estimates for each model and the top-ranked model (ΔWAIC), the Akaike weight (Weight), the standard error of the WAIC estimate (SE), and the standard error of ΔWAIC value (ΔSE)


We tested 258 I. scapularis nymphs, from which we successfully extracted DNA (determined through NanoDrop and PCR methods). Twenty-nine ticks (11.2%) were positive for B. burgdorferi, five (1.9%) were positive for Anaplasma phagocytophilum, and two (0.8%) were positive for Borrelia miyamotoi. Only one tick was co-infected with B. burgdorferi and A. phagocytophilum. We did not find any Babesia microti-positive ticks. Spatial autocorrelation in B. burgdorferi prevalence across forest fragments was not significant at any distance; at 0 m, the correlation coefficient was -0.61 (95% CI: -1.79–0.72). Within invaded forests, 11 of 75 ticks (14.6%) captured within R. multiflora and 6 of 32 ticks (18.7%) captured outside of R. multiflora were B. burgdorferi-positive. Of 151 ticks tested from uninvaded sites, 12 (7.9%) were positive for B. burgdorferi. Based on differences between posterior distributions, we estimated that there was a 64% probability that infection prevalence was greater outside of R. multiflora patches within invaded sites. However, there was a 97% probability that infection prevalence was higher among nymphs from invaded sites compared to nymphs from uninvaded sites (Fig. 2).

Fig. 2
figure 2

Distributions of differences between estimated Borrelia burgdorferi infection prevalence (posterior distributions) in nymphs between invaded and uninvaded forest fragments (a) and within invaded forests, between ticks captured within Rosa multiflora and outside of it (b). Y-axes display the density of samples from posterior distributions. Between invaded and uninvaded forests (a), there is 97% probability that B. burgdorferi prevalence is greater in invaded forests. Within invaded forests (b), there is only 64% probability that B. burgdorferi prevalence is lower within R. multiflora stands. Thus, we find support for a difference in B. burgdorferi prevalence at the forest-fragment scale (a), but not within invaded fragments (b)

Among the 50 models tested, the best models for B. burgdorferi infection prevalence (Table 3) included woody debris, leaf litter, distance to road, mouse abundance, tick abundance (all at trap-scale), and total cover (fragment-scale). Woody debris, distance to road, mouse abundance, and total cover were positively related to infection prevalence. Leaf litter and tick abundance were negatively related to infection prevalence (Fig. 3).

Fig. 3
figure 3

Mean model-averaged partial predicted responses (with 89% posterior probability intervals) of Borrelia burgdorferi prevalence among ticks (proportion of infected ticks) to six different variables: woody debris (%) (a); leaf litter volume (L/m2) (b); tick abundance (nymphs/24 h) (c); understory cover (%) (d); mouse abundance (mice per nest box per check) (e); distance to nearest road (m) (f). Overall, B. burgdorferi prevalence is predicted to increase with increasing woody debris, understory cover, mouse abundance, and distance to nearest road, while increasing leaf litter volume and tick abundance should decrease B. burgdorferi prevalence among ticks


In our comparison of invaded and uninvaded forest fragments, we found that B. burgdorferi prevalence among questing ticks did not differ within invaded forests, but that the infection prevalence in ticks from invaded forests was almost double that in ticks from uninvaded forests. Borrelia burgdorferi was the most common pathogen detected in nymphal I. scapularis from our study sites, followed by A. phagocytophilum and B. miyamotoi. Only one I. scapularis nymph was co-infected with B. burgdorferi and A. phagocytophilum, and we did not detect Ba. microti in any of the ticks tested. At finer scales within both invaded and uninvaded sites, infection prevalence was positively related to coarse woody debris, distance to the nearest road, mouse abundance, and extent of understory cover within the forest fragment. We found a negative relationship between infection prevalence and both leaf litter and tick abundance. Rosa multiflora invasion and the additional factors positively influencing pathogen prevalence point to suitable habitat characteristics for small mammal and bird hosts that are competent pathogen reservoirs.

Invaded and uninvaded fragments represent two extremes of different, degraded habitat fragment types that can be separated by the presence/absence of R. multiflora invasion in our landscape. Uninvaded sites have deep litter layers, sparse understory, high densities of questing nymphs, and relatively low infection prevalence (mean = 0.079). Invaded sites have very little leaf litter, dense understory structure, fewer questing nymphs, and roughly double the infection prevalence (mean = 0.159). Our modeling results showed that the total understory cover in a forest fragment positively influences pathogen prevalence. Understory structure, which is provided almost exclusively by invasive plants, may aggregate immature ticks and infective hosts, resulting in increased pathogen prevalence among ticks in invaded forest fragments [19,20,21]. Because B. burgdorferi is not transmitted transovarially [44], infected free-living nymphs acquire the bacteria by feeding on an infected host during their larval stage. Similarly, potential pathogen hosts must acquire B. burgdorferi by being fed upon by an infected nymph. Therefore, both immature stages of ticks must interact with infected hosts to elevate pathogen prevalence among nymphs [45].

Understory structure facilitates interactions between immature ticks and competent B. burgdorferi reservoir hosts [22, 46, 47], but see [48]. White-footed mouse (Peromyscus leucopus) and breeding bird densities are positively correlated with understory structure [47, 49,50,51] [i.e. invasive plants, in our landscape (unpublished data)]. Within invaded forests, immature ticks are aggregated in stands of invasive shrubs [11, 20, 21]. We hypothesize that larval ticks in uninvaded sites derive a greater proportion of blood meals from larger-bodied hosts that are less-competent B. burgdorferi reservoirs [52, 53]. We expect that this is in contrast to larval tick blood meals in invaded sites, which we predict are composed of a greater proportion of small-bodied hosts that are positively affected by understory structure [46] and are competent B. burgdorferi reservoirs [53,54,55]. Future work should use blood meal analysis or identification of ospC types in B. burgdorferi-positive ticks to understand how non-native plant invasion affects the interaction between specific hosts and ticks, and the resulting implications for transmission of human-invasive B. burgdorferi strains [56,57,58].

An additional hypothesis to explain greater nymphal infection prevalence in invaded sites concerns tick overwinter survival. Invaded habitats lack the litter layer that comprises suitable off-host tick habitat [11, 16, 17]. Ticks depend on the high humidity microclimate within the litter to conserve moisture and to buffer themselves from environmental fluctuations [59]. However, saturated soils coupled with extremely low temperatures may also lead to decreased overwinter survival [60]. Recent studies show that ixodid ticks infected with B. burgdorferi have greater energy reserves and are more robust to desiccation [61,62,63,64]. Therefore, the harsh litter-free environment of invaded forests may exert stronger pressure against over-winter survival of uninfected ticks, thus increasing overall infection prevalence.

The negative relationships of nymphal infection prevalence with leaf litter and tick abundance raise questions about our understanding of Lyme disease ecology in over-browsed, mature forest fragments. Uninvaded, mature forest fragments that lack understory structure have greater litter volumes and questing tick abundance than invaded forests. We hypothesize that the lack of understory structure in uninvaded fragments shifts the composition of blood meal hosts toward reservoir-incompetent species such as white-tailed deer or other large-bodied hosts [53, 65]. Talleklint & Jaenson [66] also detected a negative relationship between tick density and infection prevalence at high tick densities (> 20 nymphs/m2), which they attributed to greater roe deer (Capreolus capreolus) densities. Elevated deer densities could account for both greater tick density and lesser infection prevalence if deer act as both reproductive hosts and the dominant blood meal source [66, 67]. The close proximity among our study sites suggests that deer do not account for differences in tick abundance; most sites are close enough to be within a single deer’s home range [68,69,70] (Additional file 1: Figure S1). However, deer may reduce infection prevalence by shifting blood meals away from reservoir competent hosts that do not find suitable understory cover in over-browsed, uninvaded fragments.

The importance of invasion, habitat, and landscape variables from our models suggest that understory structure and woody debris aggregate infectious hosts and larval ticks, increasing pathogen transmission. Coarse woody debris, total understory cover, distance to road, and white-footed mouse abundance, variables that directly or indirectly represent the distribution of reservoir hosts, were positively related to infection prevalence. Coarse woody debris provides cover, nest sites, movement corridors, and foraging opportunities for immature tick hosts such as white-footed mice, Sorex and Blarina shrews, and ground-foraging birds [71,72,73,74,75]. Shrews, in particular, are often overlooked in terms of their importance in the Lyme disease system, despite evidence that they can feed and infect more ticks than white-footed mice [76]. Outside of the Pacific Northwest and southern Appalachian regions of the USA, there is a dearth of studies on habitat associations of shrews [74]; in regions where shrews have been well studied, coarse woody debris appears to be an important habitat component [77,78,79,80]. Similarly, total understory cover represents the structure available to white-footed mice and shrub-nesting birds [47, 49, 50]. The importance of distance to road suggests that perhaps small mammals and birds avoid hard edges near roads in our landscape, or at least that larval ticks encounter infectious hosts farther from roads.


Although nymphal infection prevalence was greater in invaded forests, acarological risk in terms of density of infected nymphs may be higher in uninvaded sites; the uninvaded sites examined in this study supported ~3 times as many questing nymphs compared to invaded sites [11]. Although uninvaded sites lack understory structure and therefore support lower densities of immature tick hosts, their relatively intact litter layers may allow nymphal ticks to survive longer [81] and quest more often [82], creating more opportunities to attach to humans than in invaded forests. Perhaps in uninvaded fragments, restoration of native understory structure [83] that promotes greater host diversity could reduce densities of questing infected nymphs.


  1. 1.

    Bradley CA, Altizer S. Urbanization and the ecology of wildlife diseases. Trends Ecol Evol. 2007;22:95–102.

    Article  PubMed  Google Scholar 

  2. 2.

    Allan BF, Keesing F, Ostfeld RS. Effect of forest fragmentation on Lyme disease risk. Conserv Biol. 2003;17:267–72.

    Article  Google Scholar 

  3. 3.

    LoGiudice K, Ostfeld RS, Schmidt KA, Keesing F. The ecology of infectious disease: effects of host diversity and community composition on Lyme disease risk. Proc Natl Acad Sci USA. 2003;100:567–71.

  4. 4.

    LoGiudice K, Duerr ST, Newhouse MJ, Schmidt KA, Killilea ME, Ostfeld RS. Impact of host community composition on Lyme disease risk. Ecology. 2008;89:2841–9.

    Article  PubMed  Google Scholar 

  5. 5.

    Zolnik CP, Falco RC, Kolokotronis S-O, Daniels TJ. No observed effect of landscape fragmentation on pathogen infection prevalence in blacklegged ticks (Ixodes scapularis) in the northeastern United States. PLoS One. 2015;10:e0139473.

    Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Vellend M, Verheyen K, Flinn KM, Jacquemyn H, Kolb A, Van Calster H, et al. Homogenization of forest plant communities and weakening of species-environment relationships via agricultural land use. J Ecol. 2007;95:565–73.

    Article  Google Scholar 

  7. 7.

    Huebner CD, Steinman J, Hutchinson TF, Ristau TE, Royo AA. The distribution of a non-native (Rosa multiflora) and native (Kalmia latifolia) shrub in mature closed-canopy forests across soil fertility gradients. Plant Soil. 2014;377:259–76.

    CAS  Article  Google Scholar 

  8. 8.

    Rogerson J, Globetti M, Hossler R, Moore E, Reynolds K, Hotton D, et al. Delaware Deer Management Plan. Delaware Department of Natural Resources and Environmental Control. 2010. Accessed 20 Oct 2015.

  9. 9.

    Rossell CR, Patch S, Salmons S. Effects of deer browsing on native and non-native vegetation in a mixed oak-beech forest on the Atlantic coastal plain. Northeast Nat. 2007;14:61–72.

    Article  Google Scholar 

  10. 10.

    Bernard MJ, Neatrour MA, McCay TS. Influence of soil buffering capacity on earthworm growth, survival, and community composition in the western Adirondacks and Central New York. Northeast Nat. 2009;16:269–84.

    Article  Google Scholar 

  11. 11.

    Adalsteinsson SA, D’Amico V, Shriver WG, Brisson D, Buler JJ. Scale-dependent effects of nonnative plant invasion on host-seeking tick abundance. Ecosphere. 2016;7(3):e01317.

    Article  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Lawrence B, Fisk MC, Fahey TJ, Suárez ER. Influence of nonnative earthworms on mycorrhizal colonization of sugar maple (Acer saccharum). New Phytol. 2003;157:145–53.

    Article  Google Scholar 

  13. 13.

    Suárez ER, Fahey TJ, Yavitt JB, Groffman PM, Bohlen PJ. Patterns of litter disappearance in a northern hardwood forest invaded by exotic earthworms. Ecol Appl. 2006;16:154–65.

    Article  PubMed  Google Scholar 

  14. 14.

    Hale CM, Frelich LE, Reich PB. Changes in hardwood forest understory plant communities in response to European earthworm invasions. Ecology. 2006;87:1637–49.

    Article  PubMed  Google Scholar 

  15. 15.

    Nuzzo VA, Maerz JC, Blossey B. Earthworm invasion as the driving force behind plant invasion and community change in northeastern north American forests. Conserv Biol. 2009;23:966–74.

    Article  PubMed  Google Scholar 

  16. 16.

    Schulze TL, Jordan RA, Hung RW. Suppression of subadult Ixodes scapularis (Acari: Ixodidae) following removal of leaf litter. J Med Entomol. 1995;32:730–3.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Burtis JC, Fahey TJ, Yavitt JB. Impact of invasive earthworms on Ixodes scapularis and other litter-dwelling arthropods in hardwood forests, central New York state, USA. Appl Soil Ecol. 2014;84:148–57.

    Article  Google Scholar 

  18. 18.

    Lubelczyk CB, Elias SP, Rand PW, Holman MS, Lacombe EH, Smith RP. Habitat associations of Ixodes scapularis (Acari: Ixodidae) in Maine. Environ Entomol. 2004;33:900–6.

    Article  Google Scholar 

  19. 19.

    Elias SP, Lubelczyk CB, Rand PW, Lacombe EH, Holman MS, Smith RP. Deer browse resistant exotic-invasive understory: an indicator of elevated human risk of exposure to Ixodes scapularis (Acari: Ixodidae) in southern coastal Maine woodlands. J Med Entomol. 2006;43:1142–52.

    Article  PubMed  Google Scholar 

  20. 20.

    Williams SC, Ward JS, Worthley TE, Stafford KC. Managing Japanese barberry (Ranunculales: Berberidaceae) infestations reduces blacklegged tick (Acari: Ixodidae) abundance and infection prevalence with Borrelia burgdorferi (Spirochaetales: Spirochaetaceae). Environ Entomol. 2009;38:977–84.

    Article  PubMed  Google Scholar 

  21. 21.

    Allan BF, Dutra HP, Goessling LS, Barnett K, Chase JM, Marquis RJ, et al. Invasive honeysuckle eradication reduces tick-borne disease risk by altering host dynamics. Proc Natl Acad Sci USA. 2010;107:18523–7.

  22. 22.

    Prusinski MA, Chen H, Drobnack JM, Kogut SJ, Means RG, Howard JJ, et al. Habitat structure associated with Borrelia burgdorferi prevalence in small mammals in New York state. Environ Entomol. 2006;35:308–19.

    Article  Google Scholar 

  23. 23.

    Kensinger BJ, Allan BF. Efficacy of dry ice-baited traps for sampling Amblyomma americanum (Acari: Ixodidae) varies with life stage but not habitat. J Med Entomol. 2011;48:708–11.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Falco RC, Fish D. Horizontal movement of adult Ixodes dammini (Acari: Ixodidae) attracted to CO2-baited traps. J Med Entomol. 1991;28:726–9.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Keirans JE, Durden LA. Illustrated key to nymphs of the tick genus Amblyomma (Acari: Ixodidae) found in the United States. J Med Entomol. 1998;35:489–95.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Keirans JE, Litwak TR. Pictorial key to the adults of hard ticks, family Ixodidae (Ixodida: Ixodoidea), east of the Mississippi River. J Med Entomol. 1989;26:435–48.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Durden LA, Keirans JE. Nymphs of the genus Ixodes (Acari: Ixodidae) of the United States: taxonomy, identification key, distribution, hosts, and medical/veterinary importance. Entomol Monogr. 1996;vol:1–50.

  28. 28.

    Rega C. Impacts of soil calcium availability and non-native plant invasions on an urban forest bird community. University of Delaware 2012. Accessed 2 Jun 2015.

  29. 29.

    Nudds TD. Quantifying the vegetative structure of wildlife cover. Wildl Soc Bull. 1977;5:113–7.

    Google Scholar 

  30. 30.

    Wolff JO. The effects of density, food, and interspecific inference on home range size in Peromyscus leucopus and Peromyscus maniculatus. Can J Zool. 1985;63:2657–62.

    Article  Google Scholar 

  31. 31.

    State of Delaware, Office of Management and Budget, Delaware Geographic Committee. 2007 Delaware land use and land cover. 1st Edition. Dover, Delaware: State of Delaware, Office of Management and Budget, Delaware Geographic Data Committee. 2007. Accessed 23 Jan 2014.

  32. 32.

    Nicholson MC, Mather TN. Methods for evaluating Lyme disease risks using geographic information systems and geospatial analysis. J Med Entomol. 1996;33:711–20.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Bunnell JE, Price SD, Das A, Shields TM, Glass GE. Geographic information systems and spatial analysis of adult Ixodes scapularis (Acari: Ixodidae) in the middle Atlantic region of the USA. J Med Entomol. 2003;40:570–6.

    Article  PubMed  Google Scholar 

  34. 34.

    Hojgaard A, Lukacik G, Piesman J. Detection of Borrelia burgdorferi, Anaplasma phagocytophilum and Babesia microti, with two different multiplex PCR assays. Ticks Tick-Borne Dis. 2014;5:349–51.

    Article  PubMed  Google Scholar 

  35. 35.

    Bjornstad, ON. ncf: Spatial nonparametric covariance functions. 2016.

  36. 36.

    Development Core R, Team R, Language A. Environment for statistical computing. Vienna, Austria: R Foundation for Statistical. Computing. 2014;

  37. 37.

    Stan Development Team. Stan: A C++ Library for Probability and Sampling. 2015.

  38. 38.

    Guo J, Lee D, Sakrejda K, Gabry J, Goodrich B, de Guzman J, et al. rstan: R Interface to Stan. 2016.

  39. 39.

    McElreath R. Rethinking: an R package for fitting and manipulating Bayesian models. 2016.

  40. 40.

    McElreath R. Statistical rethinking: a Bayesian course with examples in R and Stan. 1st ed. Boca Raton, FL, USA: CRC Press; 2016.

    Google Scholar 

  41. 41.

    Watanabe S. Asymptotic equivalence of Bayes cross-validation and widely applicable information criterion in singular learning theory. J Mach Learn Res. 2010;11:3571–94.

    Google Scholar 

  42. 42.

    Gelman A, Hwang J, Vehtari A. Understanding predictive information criteria for Bayesian models. Stat Comput. 2014;24:997–1016.

    Article  Google Scholar 

  43. 43.

    Yamashita T, Yamashita K, Kamimura RA. Stepwise AIC method for variable selection in linear regression. Commun Stat Theory Methods. 2007;36:2395–403.

    Article  Google Scholar 

  44. 44.

    Rollend L, Fish D, Childs JE. Transovarial transmission of Borrelia spirochetes by Ixodes scapularis: a summary of the literature and recent observations. Ticks Tick-Borne Dis. 2013;4:46–51.

    Article  PubMed  Google Scholar 

  45. 45.

    Buskirk JV, Ostfeld RS. Controlling Lyme disease by modifying the density and species composition of tick hosts. Ecol Appl. 1995;5:1133.

    Article  Google Scholar 

  46. 46.

    Adler GH, Telford SR, Wilson ML, Spielman A. Vegetation structure influences the burden of immature Ixodes dammini on its main host, Peromyscus leucopus. Parasitology. 1992;105:105–10.

    Article  PubMed  Google Scholar 

  47. 47.

    Willson MF, Comet TA. Bird communities of northern forests: ecological correlates of diversity and abundance in the understory. Condor. 1996;98:350–62.

    Article  Google Scholar 

  48. 48.

    Devevey G, Brisson D. The effect of spatial heterogenity on the aggregation of ticks on white-footed mice. Parasitology. 2012;139:915–25.

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    Adler GH, Wilson ML. Demography of a habitat generalist, the white-footed mouse, in a heterogeneous environment. Ecology. 1987;68:1785.

    Article  Google Scholar 

  50. 50.

    Leston LFV, Rodewald AD. Are urban forests ecological traps for understory birds? An examination using northern cardinals. Biol Conserv. 2006;131:566–74.

    Article  Google Scholar 

  51. 51.

    Croci S, Butet A, Georges A, Aguejdad R, Clergeau P. Small urban woodlands as biodiversity conservation hot-spot: a multi-taxon approach. Landsc Ecol. 2008;23:1171–86.

    Article  Google Scholar 

  52. 52.

    Cagnacci F, Bolzoni L, Rosà R, Carpi G, Hauffe HC, Valent M, et al. Effects of deer density on tick infestation of rodents and the hazard of tick-borne encephalitis. I: empirical assessment. Int J Parasitol. 2012;42:365–72.

    CAS  Article  PubMed  Google Scholar 

  53. 53.

    Barbour AG, Bunikis J, Fish D, Hanincová K. Association between body size and reservoir competence of mammals bearing Borrelia burgdorferi at an endemic site in the northeastern United States. Parasit Vectors. 2015;8:299.

    Article  PubMed  PubMed Central  Google Scholar 

  54. 54.

    Donahue JG, Piesman J, Spielman A. Reservoir competence of white-footed mice for Lyme disease spirochetes. Am J Trop Med Hygeine. 1987;36:92–6.

    CAS  Article  Google Scholar 

  55. 55.

    Ginsberg HS, Buckley PA, Balmforth MG, Zhioua E, Mitra S, Buckley FG. Reservoir competence of native north American birds for the Lyme disease spirochete Borrelia burgdorferi. J Med Entomol. 2005;42:445–9.

    Article  PubMed  Google Scholar 

  56. 56.

    Brisson D, Dykhuizen DE. ospC diversity in Borrelia burgdorferi: different hosts are different niches. Genetics. 2004;168:713–22.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  57. 57.

    Önder Ö, Shao W, Kemps BD, Lam H, Brisson D. Identifying sources of tick blood meals using unidentified tandem mass spectral libraries. Nat Commun. 2013;4:1746.

    Article  PubMed  PubMed Central  Google Scholar 

  58. 58.

    Vuong H, Canham CD, Fonseca DM, Brisson D, Morin PJ, Smouse PE, et al. Occurrence and transmission efficiencies of Borrelia burgdorferi ospC types in avian and mammalian wildlife. Infect Genet Evol. 2014;27:594–600.

    Article  PubMed  Google Scholar 

  59. 59.

    Stafford KC. Survival of immature Ixodes scapularis (Acari: Ixodidae) at different relative humidities. J Med Entomol. 1994;31:310–4.

    Article  PubMed  Google Scholar 

  60. 60.

    Brunner JL, Killilea M, Ostfeld RS. Overwintering survival of nymphal Ixodes scapularis (Acari: Ixodidae) under natural conditions. J Med Entomol. 2012;49:981–7.

    Article  PubMed  Google Scholar 

  61. 61.

    Naumov RL. Longevity of forest and taiga ticks (Ixodidae) infected and non-infected with Borrelia burgdorferi groups. Parazitologiya. 2003;37:527–32.

    CAS  Google Scholar 

  62. 62.

    Herrmann C, Gern L. Survival of Ixodes ricinus (Acari: Ixodidae) under challenging conditions of temperature and humidity is influenced by Borrelia burgdorferi sensu lato infection. J Med Entomol. 2010;47:1196–204.

    Article  PubMed  Google Scholar 

  63. 63.

    Herrmann C, Gern L. Search for blood or water is influenced by Borrelia burgdorferi in Ixodes ricinus. Parasit Vectors. 2015;8:6.

    Article  PubMed  PubMed Central  Google Scholar 

  64. 64.

    Herrmann C, Voordouw MJ, Gern L. Ixodes ricinus ticks infected with the causative agent of Lyme disease, Borrelia burgdorferi sensu lato, have higher energy reserves. Int J Parasitol. 2013;43:477–83.

    CAS  Article  PubMed  Google Scholar 

  65. 65.

    Telford SR III, Mather TN, Moore SI, Wilson ML, Spielman A. Incompetence of deer as reservoirs of the Lyme disease spirochete. Am J Trop Med Hygeine. 1988;39:105–9.

    Article  Google Scholar 

  66. 66.

    Talleklint L, Jaenson TGT. Relationship between Ixodes ricinus density and prevalence of infection with Borrelia-like spirochetes and density of infected ticks. J Med Entomol. 1996;33:805–11.

    CAS  Article  PubMed  Google Scholar 

  67. 67.

    Gray JS, Kahl O, Janetzki C, Stein J, Guy E. The spatial distribution of Borrelia burgdorferi-infected Ixodes ricinus in the Connemara region of county Galway, Ireland. Exp Appl Acarol. 1995;19:163–72.

    CAS  Article  PubMed  Google Scholar 

  68. 68.

    Rhoads CL, Bowman JL, Eyler B. Home range and movement rates of female exurban white-tailed deer. J Wildl Manag. 2010;74:987–94.

    Article  Google Scholar 

  69. 69.

    Etter DR, Hollis KM, Deelen TRV, Ludwig DR, Chelsvig JE, Anchor CL, et al. Survival and movements of white-tailed deer in suburban Chicago, Illinois. J Wildl Manag. 2002;66:500.

    Article  Google Scholar 

  70. 70.

    Grund MD, McAninch JB, Wiggers EP. Seasonal movements and habitat use of female white-tailed deer associated with an urban park. J Wildl Manag. 2002;66:123.

    Article  Google Scholar 

  71. 71.

    Roche BE, Schulte-Hostedde AI, Brooks RJ. Route choice by deer mice (Peromyscus maniculatus): reducing the risk of auditory detection by predators. Am Midl Nat. 1999;142:194–7.

    Article  Google Scholar 

  72. 72.

    Greenberg CH. Response of white-footed mice (Peromyscus leucopus) to coarse woody debris and microsite use in southern Appalachian treefall gaps. For Ecol Manag. 2002;164:57–66.

    Article  Google Scholar 

  73. 73.

    Lohr SM, Gauthreaux SA, Kilgo JC. Importance of coarse woody debris to avian communities in loblolly pine forests. Conserv Biol. 2002;16:767–77.

    Article  Google Scholar 

  74. 74.

    McCay TS, Komoroski MJ. Demographic responses of shrews to removal of coarse woody debris in a managed pine forest. For Ecol Manag. 2004;189:387–95.

    Article  Google Scholar 

  75. 75.

    Jones CG, Lindquist ES. Utilization of woody debris by Peromyscus leucopus in a fragmented urban forest. Southeast Nat. 2012;11:689–98.

    Article  Google Scholar 

  76. 76.

    Brisson D, Dykhuizen DE, Ostfeld RS. Conspicuous impacts of inconspicuous hosts on the Lyme disease epidemic. Proc R Soc B Biol Sci. 2008;275:227–35.

    Article  Google Scholar 

  77. 77.

    Carey AB, Johnson ML. Small mammals in managed, naturally young, and old-growth forests. Ecol Appl. 1995;5:336–52.

    Article  Google Scholar 

  78. 78.

    McCay TS, Laerm J, Menzel MA, Ford WM. Methods used to survey shrews (Insectivora: Soricidae) and the importance of forest-floor structure. Brimleyana. 1998;25:110–9.

    Google Scholar 

  79. 79.

    Brannon MP. Niche relationships of two syntopic species of shrews, Sorex fumeus and S. cinereus, in the southern Appalachian Mountains. J Mammal. 2000;81:1053–61.

    Article  Google Scholar 

  80. 80.

    Butts SR, McComb WC. Associations of forest-floor vertebrates with coarse woody debris in managed forests of western Oregon. J Wildl Manag. 2000;64:95.

    Article  Google Scholar 

  81. 81.

    Yuval B, Spielman A. Duration and regulation of the developmental cycle of Ixodes dammini (Acari: Ixodidae). J Med Entomol. 1990;27:196–201.

    CAS  Article  PubMed  Google Scholar 

  82. 82.

    Randolph SE, Storey K. Impact of microclimate on immature tick-rodent host interactions (Acari: Ixodidae): implications for parasite transmission. J Med Entomol. 1999;36:741–8.

    CAS  Article  PubMed  Google Scholar 

  83. 83.

    Morlando S, Schmidt SJ, LoGiudice K. Reduction in Lyme disease risk as an economic benefit of habitat restoration. Restor Ecol. 2012;20:498–504.

    Article  Google Scholar 

Download references


We are grateful for the cooperation of land agencies involved in this study: Newark City Parks, New Castle County Parks, Delaware State Parks, and Mt. Cuba Center. We thank the following individuals who helped collect field data: Z. Ladin, K. Handley, J. Nimmerichter, A. Lutto, K. Serno, J. Bondi, C. Piazza, L. Newton, and J. Curry. We are grateful for advice on DNA extraction protocols from S. Seifert, E. Stromdahl, and R. Nadolny, and for laboratory assistance from M. Brown. We also wish to thank C. Graham for assisting in the development of the Borrelia miyamotoi PCR assay.


The University of Delaware, USDA McIntire Stennis, and the Aeroecology Program at University of Delaware supported this work. DB was supported by a grant from the National Science Foundation (DEB 1354184).

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author upon reasonable request.

Author information




SAA, WGS, JJB, VD and DB designed the study. SAA trapped ticks, collected field data, and extracted DNA from ticks. JLB obtained permits to collect vertebrate data. WGS, JJB and VD contributed habitat and landscape data. AH performed PCR assays on DNA samples from the ticks. SAA analyzed the data. SAA, WGS and JJB wrote the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Solny A. Adalsteinsson.

Ethics declarations

Ethics approval and consent to participate

Data on white-footed mice was collected under Delaware State Scientific Collecting Permit #2013-007 W and with approval from University of Delaware’s Institutional Animal Care and Use Committee under protocol #1249, both of which were issued to JLB.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1: Figure S1.

Reproduced without modifications from Adalsteinsson et al. Ecosphere. 2016;7(3):e01317 [11] under a Creative Commons license (CC BY 3.0). Map of study area in New Castle County, Delaware. Forest cover is green; agriculture is pale yellow; blue is water; and human development is white. Fragments designated as “rose-invaded” and “uninvaded” refer to the presence or absence of Rosa multiflora invasion (TIFF 7506 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Adalsteinsson, S.A., Shriver, W.G., Hojgaard, A. et al. Multiflora rose invasion amplifies prevalence of Lyme disease pathogen, but not necessarily Lyme disease risk. Parasites Vectors 11, 54 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Lyme disease
  • Borrelia burgdorferi
  • Borrelia miyamotoi
  • Anaplasma phagocytophilum
  • Babesia microti
  • Invasive species
  • Urbanization
  • Forest fragment