History and complexity in tick-host dynamics: discrepancies between ‘real’ and ‘visible’ tick populations
© Dobson; licensee BioMed Central Ltd. 2014
Received: 1 April 2014
Accepted: 6 May 2014
Published: 19 May 2014
Numerical responses of ticks to changes in densities of their hosts can be complex and apparently unpredictable. Manipulations even of deterministic models can produce counter-intuitive results, including tick populations that either rise or fall under increasing host densities, depending on initial conditions.
In this paper I use an established simulation model to demonstrate a wide range of numerical responses to different scenarios of host changes, and to examine the basic mechanisms that drive them.
The rate and direction of change of host densities affects the extent to which questing tick numbers reflect those of their hosts. Numerical responses differ profoundly between dynamic tick-host systems and those allowed to reach equilibrium.
The key to understanding tick-host dynamics is to understand the difference between ‘real’ and ‘visible’ tick populations. An appreciation of the implications of this difference – and of the conditions that influence it - will benefit the effective interpretation of field data.
Among the global suite of arthropod disease vectors that includes mosquitoes, tsetse flies and triatomine bugs, Ixodid ticks present a distinct challenge to epidemiologists. The substantial alterations necessary to adapt the R0 (basic reproductive number) equation for tick-borne pathogens from that for pathogens vectored by insects [1, 2] reflects the many ways in which ticks differ from most insects, including a relatively insignificant degree of mobility, a characteristic of parasitising only one (hard ticks) or very few (soft ticks) hosts per instar, and a long interval between the acquisition and transmission of infective agents.
In short (and leaving aside all considerations of temperature and humidity), where hosts are abundant, ticks are more likely to complete their three blood meals and breed successfully - thereby augmenting the tick population - yet host abundance also determines the rate at which ticks are removed from the vegetation (and hence become unavailable for counting by biologists) [10–12]. The fact that host abundance drives both factors explains the otherwise potentially counter-intuitive observation that the wholesale removal of major host populations initially results in an increase in the questing tick population . The ‘momentum’ of the population, derived from egg-laying females that fed successfully before the hosts were removed, causes the usual emergence of questing ticks, which are then removed from the questing population at a much lower rate than before [10, 14]. Regardless of the presence and/or prevalence of alternative hosts, the questing population of larvae in the proceeding generation will therefore be even higher than normal. If one or more of the principal hosts for larvae remain (for example, if deer were the only hosts removed), then some of this population momentum is likely to be maintained through to the nymph stage. The life-history of Ixodes ticks, protracted as it is by diapause, means that significant nymphal activity could occur three or even four years after the removal of a major host population. Elsewhere in this paper I shall refer to this surge of questing ticks following host removal as the ‘latent questing population’, and to the extent to which the questing population exceeds levels expected, given current host density, as the ‘degree of latency’.
The contradictory forces of population augmentation (via the provision of meals to reproductive females) and questing tick removal – both imposed by hosts – will result in a variety of outcomes, in terms of questing ticks, under natural conditions. The relative influence of each pressure is likely to be determined by the density-dependent feeding success (perhaps more conveniently thought of as density-dependent mortality) of ticks. Here, density refers to the burden of ticks per host, rather than absolute density of ticks (or hosts) per se, since this mortality is mediated by the individual host’s acquired immune response [15, 16] and/or grooming/avoidance behaviour .
The fact that density dependence operates at the level of ticks per host means that if a host population steadily increases, the tick population may increase in a similar manner without any great modulating effect of density dependence. Indeed, density-dependent mortality will not rise unless the ratio of ticks to hosts increases. Superficially, as long as the tick population increase is directly caused by a host population increase, the ratio should not be altered. However, in practice, the ratio does increase, because the proportion of ticks that successfully feed before dying from moisture stress, predation or disease also goes up, since opportunities to feed become more abundant. Each egg-laying female therefore gives rise to more successfully feeding immature ticks. Whether or not this per-capita increase in survival – and thus greater ratio of ticks to hosts – leads to increased density-dependent mortality will depend upon individual host responses.
Readers should also note the apparent existence of the opposite effect under certain circumstances – density-dependent facilitation, as reported for domestic sheep  – but I do not consider this further here.
Density-dependence causes there to exist – theoretically at least - a host-species-specific threshold of host density, above which the real tick population ceases to grow (and may thereafter even decrease, if density-dependence is more than compensatory). Meanwhile, however, the per-tick host-contact rate will continue to increase with increasing host density (unless hosts adopt avoidance behaviour of tick-infested areas, as observed in domestic cattle ). At this threshold, recruitment to the questing population is overtaken by removal from it, and the questing population – though not necessarily the real population - will decline [14, 20]. Ultimately, density dependence ensures that, for many host species, the relationship between host density and tick density will not be linear across all values of the former [14, 20, 21].
There is a difference between ‘real’ and ‘visible’ tick populations, wherein the former is a measure of the reproductive potential (e.g. the number of egg-laying females), and the latter is the number of ticks that are typically available for counting – i.e. the unfed, host-seeking individuals. Both are relevant to considerations of tick-borne disease, but only ‘visible’ populations are measured in the field.
Density-dependent host responses to tick infestations may occur.
The often-protracted life cycle of Ixodid ticks causes the existence of a temporal lag between changes in the biotic environment and the numerical response of tick populations. This lag ensures that the relationship between host density and tick density is very strongly dependent on the state of the system – i.e. whether the host density is static or changing.
The aim of this paper is to demonstrate the behaviour of tick populations under different scenarios of host density change, and to explain and discuss them with reference to the above three phenomena (which are interconnected to a large extent). The material herein pertains mainly to Ixodes ricinus, but many of the principles outlined will be more widely applicable, particularly to the free-living (as opposed to nidiculous) Ixodes species, especially I. persulcatus in Eurasia and I. scapularis and I. pacificus in North America.
In order to demonstrate a range of tick responses to host population change, outputs were produced from the model of Dobson et al. . This deterministic population model is based around a modified stage-classified Leslie matrix , and simulates relative numbers of each tick instar (larva, nymph and adult) in each physiological state according to inputs of biotic and abiotic variables. Simple scenarios of host population change were applied, each of which comprised increasing, constant, then decreasing host densities. Initial host contact rates were taken from the model’s validation site that had the lowest natural values (see ) in order to demonstrate the behaviour of the system across the widest possible range of this parameter, and to ensure that the maximum rates reached were still safely within biologically realistic limits. Rates were scaled for graphical purposes such that the initial rate was 1. Climate data from the same site were used.
Changes in large host density were simulated under the following conditions:
Scenario 1. Large host contact rate increasing daily over 20 years, reaching a maximum of twice the starting contact rate; holding constant for 20 years; decreasing daily back to the starting contact rate over 20 years; holding constant for 20 years.
Scenario 2. Large host contact rate increasing daily over 20 years, reaching a maximum of five times the starting contact rate; holding constant for 20 years; decreasing daily back to the starting contact rate over 20 years; holding constant for 20 years.
Scenario 3. Large host contact rate increasing daily over 5 years, reaching a maximum of five times the starting contact rate; holding constant for 20 years; decreasing daily back to the starting contact rate over five years; holding constant for 20 years.
Scenario 4. Large host contact rate increasing daily over 5 years, reaching a maximum of twice the starting contact rate; holding constant for 20 years; decreasing daily back to the starting contact rate over five years; holding constant for 20 years.
In these plots, the questing nymph index is calculated by summing the daily questing abundance across the year (described elsewhere as the area under the graph of questing ticks, or AUTG ). (Though I do not directly discuss the implications, the reader should note that individual ticks will be counted more than once in this metric – which is equivalent to the seasonal questing profile derived from field studies - and that the lower the host contact rate, the more times an individual tick will be counted, since there is a lower likelihood that it will have been removed from the questing population since the previous day).
Results and discussion
The plots in Figure 2 illustrate various points of commonality and difference between the effects of the host-change scenarios, which are discussed in turn.
Variable relationships between questing nymphs and host contact rate
The undoubted truisms that (i) lots of hosts mean lots of ticks, and (ii) a lack of hosts means no ticks, nonetheless lie at very simple extremes that are linked by anything but a simple continuum of intermediates. A cursory glance at all plots in the fourth row of Figure 2 would suggest that there is no clear relationship between host density and tick density. However, Figure 3a demonstrates very clearly that the relationship is simple and linear over the same range of host contact rates if one looks at final densities of tick populations that have been allowed to reach equilibrium. (Here, densities are not yet high enough to trigger a density-dependent response).
This difference is critical – populations at equilibrium behave very differently from those that have not yet reached it, and the rate and direction of change in host density are both influential. The scatterplot points of the fourth row of plots in Figure 2 have been separated into those that derive from years in the simulations when host contact rates are increasing or immediately post-increase (closed circles), and those where the rates are decreasing or post-decrease (open circles). Identifying the derivation of these points – the ‘population history’, if such a term is helpful - is the first step in understanding the underlying mechanisms. The patterns and processes described below should provide a context in which these plots become less surprising than they might otherwise be.
A rising nymph population with a declining host contact rate
The sudden loss of a host of adult ticks is often followed by a surge in questing in the next few years as the momentum of the tick population carries through before fading (e.g. ). Results from Scenario 2 show that if the host population is steadily reduced, rather than suddenly removed, this surge in questing ticks may be relatively long-lasting (approximately ten years, in this case). Indeed, the rate of host removal has a direct impact on the magnitude of the surge in questing ticks, as can be seen by comparing the second row of plots in Figure 2 between Scenario pairs 2 & 3, and 1 & 4. In both cases, a reduction in host contact rate of the same extent, but different speed, results in questing nymph peaks of differing magnitudes. The explanation is largely quite simple and mechanical; the faster the reduction in host contact rate, the smaller the proportion of the latent questing population gets removed by hosts from the vegetation, leaving more to be counted. The population momentum is thus more apparent when the host decline is steeper.
Special variants of this general phenomenon appear to operate when major hosts are only removed/excluded from part of a larger area infested with ticks. As one should predict, experimental exclosures that prevent entry of deer to field sites are subject to ‘edge effects’ that mask, to a greater or lesser extent, the effect of having removed a major host of the reproductive stage of ticks [24, 25]. Some studies report increases in tick activity in exclosures relative to control plots outside, even after sufficient time has elapsed to nullify the effects of the population momentum described above. This putative phenomenon has significant implications for epidemiology, and merits further investigation. In the only meta-analysis to address questing tick amplification in exclosures, Perkins et al.  argued that there was a threshold area for deer exclosures (about 2.5 ha), below which amplification occurred. However, there are at least two good reasons to be wary of the empirical data. First, the four data points underlying the above estimation of the threshold area represent varying combinations of three instars of two species from different genera (Amblyomma americanum and I. scapularis), when it might be expected that different instars and different species would display distinct responses to deer removal. Secondly, three of the four data points are derived from studies that are either partly (N = 2) or entirely (N = 1) based on exclosures established within 1 or 2 years of tick data collection. This leaves the data extremely vulnerable to the effects of population momentum, and hence not indicative of the long- or even medium-term effect of exclosure. Nonetheless, the literature does offer more than one example of such apparent amplification (e.g. ), so its occurrence needs to be explained.
An amplification in exclosures should be counterintuitive if one starts with the basic assumption that tick populations require a host to feed the reproductive stage (i.e. adult females), and that deer and other large herbivores typically fulfil this role; relatively few adult ticks feed successfully on rodents under natural conditions [27, 28]. Of course, it should be quite quickly apparent that the movement of rodents across exclosure fences will allow for the ‘importation’ of engorged larvae that will later quest as nymphs within the exclosure (even if the vast majority that end up as questing adults in the exclosure will ultimately fail to reproduce); however, a further explanation is required to account for actual increases in – not merely persistence of - questing nymphs relative to the outside in small exclosures.
The most parsimonious mechanism for the amplification of ticks in exclosures relative to the unfenced outside is a form of constantly maintained population momentum (which is separate from the initial surge of questing ticks that will have followed immediately upon the erection of exclosure fencing). Pugliese & Rosà  provided a mathematical demonstration of the phenomenon (though they treated ticks per host rather than questing ticks as the variable of interest), but a brief – and non-mathematical - explanation here is appropriate.
There are fewer hosts inside an exclosure than outside, so there are more engorged ticks deposited in a given unit area outside exclosures than inside. Hosts that are able to cross the fence – e.g. rodents and small birds - import and export engorged ticks. After development to the next stage, ticks dropped in exclosures are less likely to be picked up by hosts than those that were dropped outside, meaning that they are more likely to be available for counting at any given time. There are therefore two pressures acting in opposite directions – low host numbers inside exclosures mean fewer opportunities to feed, but also a greater likelihood of being counted – a higher degree of ‘latency’. With no movement of hosts across exclosure fences, the tick population either side would be broadly proportional to host density. However, this situation begins to change if hosts move across fences, since the import of ticks should exceed export (given that hosts moving out of an exclosure will likely have picked up fewer ticks than those moving in, assuming initially higher densities of ticks outside). Exported ticks are more likely to find hosts than imported ticks, because of higher host densities outside, so the number of questing ticks available to be counted in exclosures will increase relative to that outside. Tick densities in exclosures can continue to rise relative to the outside even if import/export becomes balanced, as long as overall host densities remain lower in exclosures, since this will result in more unfed ticks in exclosures.
Amplification in exclosures is therefore theoretically feasible, but Perkins et al.’s  suggestion of a general threshold area below which it occurs is unrealistic, on the same, purely theoretical grounds. Indeed, the extent of amplification of questing tick numbers in exclosures ought to be dependent upon four factors, the first two of which – degree of fence-crossing by small hosts and the exclosure size – are inversely related to one another and can be considered as alternative, reciprocal measures of the same fundamental process . The smaller the exclosure and the greater the degree of fence-crossing, the greater the potential for import of engorged ticks which will then produce questing ticks of the next stage that are less likely to be removed by hosts (higher latency) than those that were dropped outside.
The third factor is density dependence. As mentioned in the Introduction, the strength of this effect determines the point of balance between the two host-mediated forces - population augmentation (via the provision of meals to reproductive females) and questing tick reduction (via successful feeding and therefore removal from the questing population). To re-cap: Inside the exclosure, removal by hosts is lower (because there are fewer hosts), but augmentation via successfully fed adults is also lower (because adults have fewer feeding opportunities). When density-dependent mortality is operating (i.e. when lower host densities mean higher per-tick feeding success inside than outside), the population augmentation force becomes relatively less dominant than the questing tick removal force. The discrepancy in host density between inside and outside therefore has less influence on tick reproduction, resulting in a further relative increase in questing numbers within the exclosure.
In their two-host system, Pugliese & Rosà  also demonstrated an effect of the fourth factor, the contact rate between ticks and the larger (i.e. non-fence-crossing) hosts; the higher the large host contact rate, the higher the likelihood of amplification inside exclosures. This parameter is better thought of as the ratio of tick-host contact rate inside and outside exclosures; since large hosts do not enter exclosures, raising the large host contact rate means that the host contact rate overall becomes relatively higher outside, but is unchanged inside. This causes greater tick population augmentation outside, and hence a greater degree of engorged tick import into exclosures, whilst simultaneously increasing unfed tick removal outside (masking the population augmentation) without altering it inside (thereby allowing imported ticks to result in greater questing numbers).
It should be obvious, therefore, that whilst there ought to be a threshold point for any scenario, above which exclosures show amplified tick densities, this threshold will vary according to the exclosure area, degree of fence-crossing, nature of density dependence and ratio of host contact rates inside and outside. The shape of the exclosure (and more specifically, the edge-to-area ratio) should also be influential, since it will affect patterns of fence-crossing and the spatial distribution of ticks across the exclosure.
‘Real’ versus ‘visible’ populations
Comparison of the laying adult and questing nymph plots in the second row of Figure 2 indicates that the number of laying adults follows the host contact rate (top row of plots) in a far more consistent manner than does the number of questing nymphs, particularly in the time period marked ‘c’. All other things being equal (as they are in these simulations), host contact rate should drive the tick population in a dominant manner – with some discrepancy due to the effects of density dependence – so the broad agreement between laying adults and host contact rate confirms that the former represents a good metric of the ‘real’ population (as opposed to the ‘visible’, questing population). This also means that the ratio of laying adults to questing nymphs (third row of plots) is a function of the discrepancy between real and questing populations. The ratio provides another way of interpreting and rationalising the peak in questing nymphs that continues as host contact rates decline (sections ‘c’, second row of plots). Across all scenarios, the size of this nymphal peak is proportional to the magnitude of the change in the ratio of laying adults to questing nymphs during the earlier phase of host density increase (sections ‘a’, third row of plots). In other words, the disparity between real and questing populations gives a measure of the population’s ‘capital’, which is released as population momentum when host densities decrease. Thinking in these terms helps to explain why the questing nymph peak is higher when the host contact rate decline is faster; the capital is released over a much shorter time period in Scenarios 3 & 4 than in 1 & 2, hence the consequent surge of questing ticks is more acute.
Synthesis and applications
“All models are wrong, but some are useful” ; here, population models are useful because they give us the opportunity to investigate demographic processes in a way that is effectively impossible in an ordinary empirical study. The model of Dobson et al.  is among a small number of models that attempt to capture the convoluted life-history of Ixodes spp. and formally characterise population processes in simple mathematical terms. Providing that they are structured in a manner that sufficiently reflects the realities of the system they aim to simulate, these models may be thought of as strictly controlled experimental environments in which the effects of individual factors may be assessed in isolation . Here, this rarefied experimental environment has illuminated conditions that we must take into account when interpreting tick field data - namely the influence of ‘population history’, whereby the direction of response to host density change may be influenced by previous, unmeasured states.
The dominant issue outlined in this paper – the fact that reliable measurement of questing ticks will not necessarily lead to a reliable understanding of the population under study – is not one that may be practically resolved in the field. Questing ticks can be accurately surveyed by blanket dragging, assuming that an appropriate sampling regime is employed , and there are various techniques for measuring the abundance of different host species, but there is no realistic way of recording the density of laying adult female ticks without going through the extremely laborious and inaccurate process of wholesale extraction of large samples of vegetation and topsoil and picking through it to look for them.
It is possible that some insight into the problems caused by population history may be gained from looking at precise patterns of questing ticks and those on hosts through the year, without reverting to direct examination of ‘real’ population size. Any mechanistically realistic population model would allow detailed examination of these variables and thereby potentially aid the correct interpretation of confusing field data; whether or not enough empirical data yet exist for model parameterisation remains to be seen. Models that explicitly incorporate host individuals (as opposed to host contact rates for tick individuals) would be capable of revealing the epidemiological implications of different values of the ratio of laying adults to questing nymphs (i.e. the ratio of real to visible population size), and should be a focus for future study.
One facet of this particular population model may not be universally applicable, and would benefit from further empirical research: Density dependence is here depicted as a linear response, after a certain threshold has been reached. This is not necessarily realistic. If tick feeding were highly clustered on a few hosts, a threshold model of host response could result in greatly over-compensating density dependence, and possibly fluctuating tick population dynamics, at almost any value of tick-host ratio. It is hoped that future incarnations of the population model will incorporate such mechanisms, following additional empirical input.
In many cases, of course, a given tick population under natural conditions will not experience great fluctuations in numbers of its hosts (the obvious exception being the cyclic dynamics of microtine rodents), meaning that most differences in population size over time will be determined more by a subtle shifting of the host community structure - or some abiotic factor – than by the sorts of processes outlined here. However, the majority of tick-borne disease research interest is currently focussed upon areas where ecological change is, or is predicted to be, extensive – especially regarding the question of whether biodiversity loss will amplify the density of infected ticks [30–32]. Data from field sites where host densities have not been stable for at least ten years may lead researchers to erroneous conclusions if they do not have an appreciation of the considerable complexity of tick population biology. Specifically, tick biologists must acknowledge that numerical responses of ticks to changes in host density are not necessarily consistent between populations. As a consequence, it will be almost impossible to infer much information from ‘snapshot’ sampling of a tick population and those of its hosts, especially if any recent change is suspected in the latter. It is also fundamentally important to recognise that data points drawn from the same population in different years are unlikely to be statistically independent unless the temporal separation is large. Likewise, a set of tick abundance samples taken simultaneously from several separate sites may not be meaningfully related to the corresponding set of host density (or meteorological) data.
Tick population biology is complex. Changes in host density may have counterintuitive effects on the abundance of questing ticks; the key to understanding these numerical responses is to appreciate the difference between real and visible tick populations. An understanding of the suite of mechanisms that drive tick population dynamics – and the range of outcomes they may produce – will aid interpretation of data from field studies, and should warn against over-extrapolation from any particular set of results. Further research on epidemiological implications is required.
I am very grateful to Sarah Randolph, who read several early drafts and made numerous essential suggestions. The comments of two anonymous referees led to a much improved final manuscript. The work was supported by an Impact research fellowship from the University of Stirling; development of the tick population model was funded by the UK Department of the Environment and Rural Affairs (Defra) (project SE4105), the Research Councils’ Rural Economy and Land Use (RELU) Programme (funded by the Economic and Social Research Council, Biotechnology and Biological Sciences Research Council and Natural Environment Research Council, with additional funding from Defra and the Scottish Executive Environment and Rural Affairs Department) (project RES-229-25-0007) and EU grant GOCE-2003-010284 EDEN.
- Randolph SE: Ticks are not insects: consequences of contrasting vector biology for transmission potential. Parasitol Today. 1998, 14: 186-192. 10.1016/S0169-4758(98)01224-1.View ArticlePubMedGoogle Scholar
- Hartemink NA, Randolph SE, Davis SA, Heesterbeek JAP: The basic reproduction number for complex disease systems: defining R 0 for tick-borne infections. Am Nat. 2008, 171: 743-754. 10.1086/587530.View ArticlePubMedGoogle Scholar
- Dantas-Torres F, Chomel BB, Otranto D: Ticks and tick-borne diseases: a one health perspective. Trends Parasitol. 2012, 28: 437-446. 10.1016/j.pt.2012.07.003.View ArticlePubMedGoogle Scholar
- Biology of Ticks Volume 2. Edited by: Sonenshine DE, Roe RM. 2014, New York: Oxford University Press, 2Google Scholar
- Eisen L, Eisen RJ, Lane RS: Seasonal activity patterns of Ixodes pacificus nymphs in relation to climate. Med Vet Ent. 2002, 16: 235-244. 10.1046/j.1365-2915.2002.00372.x.View ArticleGoogle Scholar
- Estrada-Peña A, Venzal JM, Sanchez Acedo C: The tick Ixodes ricinus: distribution and climate preferences in the western Palaearctic. Med Vet Ent. 2006, 20: 189-197. 10.1111/j.1365-2915.2006.00622.x.View ArticleGoogle Scholar
- Knap N, Durmiši E, Saksida A, Korva M, Petrovec M, Avšič-Županc T: Influence of climatic factors on dynamics of questing Ixodes ricinus ticks in Slovenia. Vet Parasitol. 2009, 164: 275-281. 10.1016/j.vetpar.2009.06.001.View ArticlePubMedGoogle Scholar
- Medlock JM, Hansford KM, Bormane A, Derdakova M, Estrada-Peña A, George J-C, Golovljova I, Jaenson TGT, Jensen JK, Jensen PM, Kazimirova M, Oteo JA, Papa A, Pfister K, Plantard O, Randolph SE, Rizzoli A, Santos-Silva MM, Sprong H, Vial L, Hendrickx G, Zeller H, Van Bortel W: Driving forces for changes in geographical distribution of Ixodes ricinus ticks in Europe. Parasit Vectors. 2013, 6: 1-10.1186/1756-3305-6-1.PubMed CentralView ArticlePubMedGoogle Scholar
- Dobson ADM: Ticks in the wrong boxes: assessing error in blanket-drag studies due to occasional sampling. Parasit Vectors. 2013, 6: 344-10.1186/1756-3305-6-344.PubMed CentralView ArticlePubMedGoogle Scholar
- Ginsberg HS, Zhioua E: Influence of deer abundance on the abundance of questing adult Ixodes scapularis (Acari: Ixodidae). J Med Ent. 1999, 36: 376-381.View ArticleGoogle Scholar
- Randolph SE: Tick ecology: processes and patterns behind the epidemiological risk posed by ixodid ticks as vectors. Parasitol. 2004, 129: S37-S66. 10.1017/S0031182004004925.View ArticleGoogle Scholar
- Bolzoni L, Rosà R, Cagnacci F, Rizzoli A: Effect of deer density on tick infestation of rodents and the hazard of tick-borne encephalitis. II: population and infection models. Int J Parasitol. 2012, 42: 373-381. 10.1016/j.ijpara.2012.02.006.View ArticlePubMedGoogle Scholar
- Rand PW, Lubelczyk C, Holman MS, Lacombe EH, Smith RP: Abundance of ixodes scapularis (acari: ixodidae) after the complete removal of deer from an isolated offshore island, endemic for lyme disease. J Med Ent. 2004, 41: 779-784. 10.1603/0022-2585-41.4.779.View ArticleGoogle Scholar
- Dobson ADM, Randolph SE: Modelling the effects of recent changes in climate, host density and acaricide treatments on population dynamics of Ixodes ricinus in the UK. J Appl Ecol. 2011, 48: 1029-1037. 10.1111/j.1365-2664.2011.02004.x.View ArticleGoogle Scholar
- Sutherst RW, Wagland BM, Roberts JA: Effect of density on survival of Boophilus microplus on previously unexposed cattle. Int J Parasitol. 1978, 8: 321-324. 10.1016/0020-7519(78)90097-8.View ArticlePubMedGoogle Scholar
- Randolph SE: Population dynamics and density-dependent seasonal mortality indices of the tick Rhipicephalus appendiculatus in eastern and southern Africa. Med Vet Ent. 1994, 8: 351-368. 10.1111/j.1365-2915.1994.tb00099.x.View ArticleGoogle Scholar
- Levin ML, Fish D: Density-dependent factors regulating feeding success of Ixodes scapularis larvae (Acari: Ixodidae). J Parasitol. 1998, 84: 36-43. 10.2307/3284526.View ArticlePubMedGoogle Scholar
- Ogden NH, Casey AN, French NP, Adams JD, Woldehiwet Z: Field evidence for density-dependent facilitation amongst Ixodes ricinus ticks feeding on sheep. Parasitol. 2002, 124: 117-125.Google Scholar
- Sutherst RW, Floyd RB, Bourne AS, Dallwitz MJ: Cattle grazing behavior regulates tick populations. Experientia. 1986, 42: 194-196. 10.1007/BF01952465.View ArticleGoogle Scholar
- Rosà R, Pugliese A: Effects of tick population dynamics and host densities on the persistence of tick-borne infections. Math Biosci. 2007, 208: 216-240. 10.1016/j.mbs.2006.10.002.View ArticlePubMedGoogle Scholar
- Pugliese A, Rosà R: Effect of host populations on the intensity of ticks and the prevalence of tick-borne pathogens: how to interpret the results of deer exclosure experiments. Parasitol. 2008, 135: 1531-1544. 10.1017/S003118200800036X.View ArticleGoogle Scholar
- Dobson ADM, Finnie TJR, Randolph SE: A modified matrix model to describe the seasonal population ecology of the European tick Ixodes ricinus. J Appl Ecol. 2011, 48: 1017-1028. 10.1111/j.1365-2664.2011.02003.x.View ArticleGoogle Scholar
- Leslie PH: On the use of matrices in certain population mathematics. Biometrika. 1945, 33: 183-212. 10.1093/biomet/33.3.183.View ArticlePubMedGoogle Scholar
- Daniels TJ, Fish D, Schwartz I: Reduced abundance of Ixodes scapularis (Acari: Ixodidae) and Lyme disease risk by deer exclusion. J Med Ent. 1993, 30: 1043-1049.View ArticleGoogle Scholar
- Daniels TJ, Fish D: Effect of deer exclusion on the abundance of immature Ixodes scapularis (Acari: Ixodidae) parasitizing small and medium-sized mammals. J Med Ent. 1995, 32: 5-11.View ArticleGoogle Scholar
- Perkins SE, Cattadori IM, Tagliapietra V, Rizzoli AP, Hudson PJ: Localised deer absence leads to tick amplification. Ecology. 2006, 87 (8): 1981-1986. 10.1890/0012-9658(2006)87[1981:LDALTT]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Schulze TL, Bowen GS, Lakat MF, Parkin WE, Shisler JK: Seasonal abundance and hosts of Ixodes dammini (acari: ixodidae) and other ixodid ticks from an endemic lyme disease focus in New jersey, USA. J Med Ent. 1986, 23 (1): 105-109.View ArticleGoogle Scholar
- Talleklint L, Jaenson TGT: Transmission of Borrelia burgdorferi s.l. from mammal reservoirs to the primary vector of Lyme borreliosis, Ixodes ricinus (Acari: Ixodidae), in Sweden. J Med Ent. 1994, 31: 880-886.View ArticleGoogle Scholar
- Box GEP, Draper NR: Empirical Model Building and Response Surfaces. 1987, New York: John Wiley & SonsGoogle Scholar
- Ogden NH, Tsao JI: Biodiversity and Lyme disease: dilution or amplification?. Epidemics. 2009, 1: 196-206. 10.1016/j.epidem.2009.06.002.View ArticlePubMedGoogle Scholar
- Randolph SE, Dobson ADM: Pangloss revisited: a critique of the dilution effect and the biodiversity-buffers-disease paradigm. Parasitol. 2012, 139: 847-863. 10.1017/S0031182012000200.View ArticleGoogle Scholar
- Wood CL, Lafferty K: Biodiversity and disease: a synthesis of ecological perspectives on Lyme disease transmission. Trends Ecol Evol. 2013, 28: 239-247. 10.1016/j.tree.2012.10.011.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.