The influence of life history characteristics on flea (Siphonaptera) species distribution models
© van der Mescht et al. 2016
Received: 13 February 2016
Accepted: 21 March 2016
Published: 29 March 2016
Ectoparasites exhibit pronounced variation in life history characteristics such as time spent on the host and host range. Since contemporary species distribution (SD) modelling does not account for differences in life history, the accuracy of predictions of current and future species’ ranges could differ significantly between life history groups.
SD model performance was compared between 21 flea species that differ in microhabitat preferences and level of host specificity. Distribution models generally performed well, with no significant differences in model performance based on either microhabitat preferences or host specificity. However, the relative importance of predictor variables was significantly related to host specificity, with the distribution of host-opportunistic fleas strongly limited by thermal conditions and host-specific fleas more associated with conditions that restrict their hosts’ distribution. The importance of temperature was even more pronounced when considering microhabitat preference, with the distribution of fur fleas being strongly limited by thermal conditions and nest fleas more associated with variables that affect microclimatic conditions in the host nest.
Contemporary SD modelling, that includes climate and landscape variables, is a valuable tool to study the biogeography and future distributions of fleas and other parasites taxa. However, consideration of life history characteristics is cautioned as species may be differentially sensitive to environmental conditions.
Ectoparasites exhibit pronounced variation in life history strategies with parasite-host associations ranging from one-to-one symbiosis (host-specific) to multi-partner symbiosis (host generalist). Furthermore, life history also differs between taxa with some parasite species being only temporarily associated with the body of a host (e.g. ticks and fleas), while others are more permanently linked with a host (e.g. lice) [1–3]. Consequently, it can be argued that the distribution of parasite species with multiple free-living stages (i.e. temporary parasites) are likely more strongly affected by the off-host environment (e.g. climatic and landscape features), whereas the distribution of permanent parasites may be indirectly driven by factors affecting host assemblages (e.g. shelter and food) (see [1, 2]). Temporary parasite taxa, however, are also characterized by varying levels of host association [1, 2] with these subtle differences in life history characteristics also potentially adding further complexity to patterns of parasite species distributions [2, 4]. These complex relationships between parasite, host and environment may make it difficult to achieve accurate range predictions for ectoparasites.
The increasing threat posed by emerging infectious diseases [5, 6] coupled with an increase in the availability of species occurrence records, has stimulated renewed interest in predicting the current and future distributions of arthropod vectors. Species distribution (SD) modelling has proved useful for this purpose, e.g. [7–9], with particular success for several medically and veterinary important arthropod vectors, e.g. [8, 10–12], including a single study on fleas . In the latter study the authors used standard climatic variables (temperature, relative humidity and precipitation) in a GARP modelling approach to estimate the regional distribution of 18 flea species that act as vectors of Yersinia pestis, the bacterial agent of plague in California, USA. Given the recent re-emergence of certain flea-borne diseases  any improvements in our understanding of changes in flea vector distributions is valuable for the field of epidemiology.
Fleas are obligate ectoparasites of terrestrial vertebrates, and are regarded as "permanent satellites" of their hosts, due to the intimate association between fleas and hosts [13, 14]. In general, fleas spend part of their life-cycle in the host’s nest (egg, larvae and pupae) while adults occur on the body of the host. The length of time that adults spend on the host varies between flea taxa [1, 14, 15] and this difference in microhabitat preference allows fleas to be categorized either as “fur” (adults spend more time on the host), “nest” (adults spend more time in the nest of the host) or “fur/nest” species (adults spend more or less equal amounts on the host and in the nest of the host) [1, 14, 15]. Flea species also differ in terms of host specificity, which ranges from host-specific (recorded from ≤ 2 host species) to host-opportunistic (recorded from > 2 host species) [1, 14, 16]. There appears to be no relationship between microhabitat preference and level of host specificity exhibited by flea species (i.e. nest fleas are not generally regarded as having a higher level of host specificity). However, differences in level of host association between a flea and a host may have profound implications for the level of exposure to environmental features (e.g. climate and landscape), see [17–20], highlighting the need to assess the importance of these life history characteristics in SD modelling studies.
The aim of this study was therefore to compare SD model performance and the relative importance of predictor variables between flea species with different microhabitat preferences (fur vs nest) and level of host specificity (opportunistic vs specific). We predicted that fur fleas will be more accurately modelled due to being more strongly affected by variables associated with regional environmental conditions (e.g. climate), while nest fleas will be less accurately modelled due to being affected by conditions within the host nest (e.g. soil conditions and microclimate; since all life stages spend the majority of their life-cycle off the host and have limited dispersal capabilities). Furthermore, we predict that host-specific fleas will be more strongly associated with the abiotic variables constraining their host’s distribution and therefore will be more accurately modelled than host-opportunistic fleas. This is expected because host-specific fleas are assumed to be adapted to the immediate environment of their specific host and thus are expected to tolerate a narrower range of physical conditions compared to host-opportunistic fleas (see [14, 21, 22]).
Microhabitat preference, host specificity, number of unique occurrence records and two measures of model performance (AUC and TSS) for each flea species
Chiastopsylla mulleri simplex
Listropsylla chelura chelura
Preliminary climate and landscape variables were selected based on our knowledge of flea ecology, limiting candidate variables to only include predictors that are considered ecologically relevant to flea species (following [30–32]). All predictor variables and flea occurrence data were converted to Quarter Degree Grid Cell (QDGC) scale and cropped to the borders of South Africa. Five remotely-sensed climate-based variables (daytime land surface temperature (hereafter referred to as day temperature), Leaf Area Index (LAI), Normalised Difference Vegetation Index (NDVI), rainfall, water vapour, and soil characteristics) and one landscape feature variable (Topography) were extracted from the NASA-NEO website (http://neo.sci.gsfc.nasa.gov/about/) as potential predictor variables (missing values were estimated as the average of contiguous cells). Climate is known to generally influence flea populations to a greater extent than host species, especially at regional and local scales [20, 33], with air temperature, rainfall and relative humidity being important for flea survival, see [1, 10, 33–36]. NDVI is widely used in arthropod vector distribution modelling and is a measure of primary productivity (plant photosynthetic activity), and therefore can be considered as a proxy for general arthropod habitat conditions [37, 38]. Furthermore, NDVI has also been successfully used in small mammal resource and population dynamics studies [38–40] and therefore may also be a surrogate for host availability. LAI is a measure of plant canopy structure and can influence incident radiation and evapotranspiration at the soil surface . Additionally soil data, including soil organic carbon content, pH, cation exchange capacity, percentage sand and bulk density, were extracted from the SoilGrids database  (http://www.soilgrids.org/) at a depth of 60–100 cm. These soil characteristics may be expected to have direct (via microhabitat) and indirect (small mammal burrowing conditions) effects on flea species distributions [35, 43–46].
To remove collinearity between climate-based predictors and to summarize seasonality, we performed harmonic regressions for all of the climate variables using monthly data from January 2001 until June 2014 (following the methods of ) (see Additional file 1 for R script). These Fourier-transformed variables (“harmonic variables” hereafter) represent key temporal climate trends, reflecting different measures of seasonality . Maximum, minimum, mean, range, and coefficient of variation (CV) values were also calculated for each climate variable. The correlation between all predictor variables was then calculated to identify collinear predictors, with the most strongly correlated variable excluded and the process repeated until the strongest correlation was weaker than |0.7|. When choosing between two strongly correlated variables, the preference was to drop harmonic variables rather than the other more biologically-interpretable variables (i.e. maximum, minimum, mean, range, CV, topography, and soil variables). Through this process an initial set of 68 predictor variables was reduced to 19 predictors with minimal collinearity (see Additional files 2 and 3) and with clear ecological relevance, avoiding the inclusion of variables that are irrelevant and/or will inflate models [7, 9, 47, 48]. The raw data associated with predictor variables are available upon request from the authors. All analyses were conducted in R v3.1.3  and ArcGIS 10.1 .
Species distribution modelling relates species presence (or presence and absence) data to environmental variables to predict the distribution of species over a specified geographic range . In this study MaxEnt models, based on the maximum entropy algorithm, were used with presence-only data (MaxEnt v3.3.3; [52, 53]). To account for potentially spatially-biased sampling of fleas across our study region (e.g. due to more studies being conducted in protected areas), MaxEnt models were adjusted for uneven sampling by incorporating background data reflecting patterns in sampling effort [31, 32] (see Fig. 1). All 1064 localities from which flea species occurrences have been published were therefore included as background points to distinguish false absences (lack of survey data) from true absences (species was not recorded) (Fig. 1).
For each flea species variable importance (i.e. relative contribution of each predictor variable) was calculated using the full dataset. The mean importance of each variable type (the average of all related individual predictor variables) was determined for each group of flea species (microhabitat preference and host specificity). Model performance was evaluated using 10-fold cross-validation to calculate the area under the curve (AUC) of the receiver operating characteristic (ROC; ) in MaxEnt and the true skill statistic calculated from MaxEnt output in R (TSS; , see Additional file 4). Analysis of variation (ANOVA) was used to test if AUC and TSS values differed according to species’ microhabitat preference and host specificity. Non-metric multi-dimensional scaling (NMDS) and analysis of similarity (ANOSIM; implemented from the vegan package in R) were used to test if differences in variable importance were related to microhabitat preference or host specificity.
A total of 21 flea species from small mammals were selected with known microhabitat preference, host specificity, and sufficient occurrence records (Table 1). Overall, model performance was highly variable between flea species (AUC from 0.545 to 0.978; TSS from 0.067 to 0.804; Table 1), but was good to excellent on average (mean ± SE, AUC = 0.799 ± 0.026; TSS = 0.464 ± 0.045).
Final list of variable type and individual predictor variables used for modelling flea species with different life histories
6th harmonic component of daytime land surface temperature
7th harmonic component of daytime land surface temperature
Minimum daytime land surface temperature
3rd harmonic component of daily LAI
4th harmonic component of daily LAI
4th harmonic component of daily NDVI
7th harmonic component of daily NDVI
4th harmonic component of daily rainfall
6th harmonic component of daily rainfall
7th harmonic component of daily rainfall
4th harmonic component of daily water vapour
5th harmonic component of daily water vapour
6th harmonic component of daily water vapour
Minimum water vapour
Soil percentage sand
Soil organic carbon
The study confirms that contemporary SD modelling, that includes climatic and landscape variables, has potential for improving predictions of changes in the distribution of flea species. The model performance for fleas was good overall, with life history not having a significant effect on model performance. However, the importance of predictor variables differed considerably between species with different life history strategies, suggesting differential sensitivity to climate (temperature and rainfall) and landscape feature (NDVI) variables among groups of fleas.
Abiotic conditions such as air temperature, relative humidity and precipitation are important for fleas [1, 10, 33–36] especially at the regional scale [20, 56]. In particular, air temperature and relative humidity can have a direct (influence on egg production and rate of oviposition, development and survival of pre-imaginal fleas, and survival of imagoes) [34–36] and/or indirect affect (influence hosts abundance and distribution) on flea survival [1, 14, 20]. However, flea species also vary in terms of preferred temperature range (reviewed by ) and although there is evidence that phylogenetic relatedness may explain this pattern for certain flea taxa it appears not to be consistent across all taxa [14, 20, 57]. In the present study we found that temperature is an important predictor of flea species ranges, irrespective of life history, but that the relative importance of temperature differed between flea life histories. In particular, the higher importance of temperature for host-opportunistic compared to host-specific fleas may be related to their physiological limits [34–36]. Host-opportunistic fleas generally have a wider tolerance range (climate and host composition), are geographically more widespread and are likely to experience more environmental heterogeneity throughout their geographical range [14, 20–22]. In contrast, host-specific fleas generally have narrower tolerance ranges (climate and host composition), are less widespread and are likely to experience more environmental homogeneity throughout their geographical range [14, 20–22]. Thus, although temperature is important, especially for the development of immature nest stages in general, it is evidently less important for host-specific fleas in this study.
Host-specific fleas seems to be more influenced by variables related to host availability as evident from the higher NDVI importance for host-specific compared to host-opportunistic fleas in the current study. NDVI has been used as a measure of the amount of suitable habitat for arthropod vectors [37, 38] and may also be a good proxy for several aspects of habitat quality that are of relevance to fleas. Abiotic factors influence vegetation and thus food supply [58–62] and actual or perceived predation risk [63–65] which will affect small mammal host abundance and distribution [38–40], this in turn may directly influence flea abundance and distribution . In the current study, environmental stability could facilitate small mammal host population stability and subsequent flea population stability and specialization [66, 67], because host-specific fleas tolerate restricted abiotic and biotic conditions . In contrast, host-opportunistic fleas can tolerate variable environments (e.g. seasonal environments that have greater variability in NDVI) that could facilitate less stable small mammal host populations and are thus less affected by NDVI .
Predictor importance for individual variables also differed between species with different microhabitat preferences. Temperature was again overall the most important predictor for both fur and nest fleas. However, the difference in the importance of temperature was even more pronounced when compared with level of host-specificity (as discussed above), with temperature being more important for fur fleas as opposed to nest fleas. Fur fleas spend more time on the body of the host and as a result are potentially exposed to higher fluctuations in ambient temperature and relative humidity compared to nest fleas (see ). It has been suggested that nest fleas may have evolved to spend more time in nests due to more constant and buffered microclimatic conditions brought about by the physical properties of nests [1, 69, 70]. Although discrepancies in the fur versus nest dichotomy are mainly attributed to ambient temperature in the literature (see  and references within), it has also been suggested that within-host among-flea difference can be explained by this dichotomy, whereas between-host within-flea differences are better explained by between-host difference in nest construction [71, 72].
Studies on different nest types (burrows, above-ground nests and nest within rock crevices) all demonstrate that microclimatic conditions are more stable in nests compared to the external environment [43–46]. Other factors that can facilitate higher and more stable humidity levels in the nest include, the presence of nest material [46, 73] and higher soil water capacity of the mineral and organic enriched nest soil , due to the activities of hosts in nests [46, 73]. In support of this, rainfall contributed significantly more towards explaining the distribution of nest compared to fur fleas in the current study. Rainwater naturally filters down into soil layers which can contribute to maintaining higher and more stable humidity levels in the nest of hosts [43, 45, 46]. The type and complexity of host nests are influenced by the soil texture [46, 75]. It is therefore not surprising that soil was more important for nest compared to fur fleas. It is evident from this study that the interplay between temperature, rainfall and potentially relative humidity with nest construction can facilitate the separation of fleas into different microhabitat types.
Our study suggests that SD modelling can be a useful tool for studying the drivers of flea species distributions and also the underlying ecology of these species, but caution needs to be taken when deciding which predictor variables to include. While our results highlight how contemporary models can perform well, it is unclear to what extent the inclusion of biotic interactions (e.g. host availability and competition) could further improve model accuracy and transferability [7, 56]. Specifically, modelling and comparing flea species with different levels of host-specificity could possibly benefit from including accurate host species occurrence data. Furthermore, it is important to remember that ecological patterns are affected by processes that act at different scales . For example, the assembly of flea compound (all species infesting a host community) communities is strongly affected by host filters (e.g. evolutionary, biogeographic and historical forces) at the continental scale, while at the regional and local scale it is more strongly affected by the abiotic filters (e.g. topography, NDVI, and climate) (see [17–20]. Therefore, the value of the inclusion of host species data (and the predominant importance of temperature in our models) may be contingent on the spatial scale of analysis.
Five of the flea species in our study (Chiastopsylla rossi, Dinopsyllus lypusus, Listropsylla dorripae, Xenopsylla pirei, and Xenopsylla versuta) have been implicated as possible vectors of diseases in South Africa (e.g. plague ). As a consequence, accurate forecasts of the future distributions of these species are valuable for the field of epidemiology. Our results suggest that despite differences in their degree of host specificity, SD models should perform well for all of these species. However, due to differential sensitivity to different groups of climatic and landscape variables, host-specific and generalist flea species are likely to respond very differently to changes in abiotic conditions. As a result, our results suggest the importance of explicitly considering species life history as a potential mediating variable when predicting flea species distributions.
Private landowners and nature conservation agencies are thanked for granting permission to trap on their property or in their reserves with the following permit numbers (Western Cape, 0035-AAA007-00423; Northern Cape, FAUNA 1076/2011; Eastern Cape, CRO37/11CR; KZN wildlife, OP4990/2010 and Gauteng, CPF 6–0153). Adriaan Engelbrecht, Karlien Malan and Jannie Groenewald are thanked for assistance in the field. The project was made possible by funding from the National Research Foundation (NRF) through research grants (GUN 80764 and 85718) and a scholarship to LVDM, and Stellenbosch University. The Grant holder acknowledges that opinions, findings and conclusions or recommendations expressed in any publication generated by the NRF-supported research are those of the authors, and that the NRF accepts no liability whatsoever in this regard. We would like to thank two anonymous reviewers for their helpful comments.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), 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 (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Marshall AG. The ecology of ectoparasite insects. London: Academic; 1981.Google Scholar
- Morand S, Krasnov BR, Poulin R, editors. Micromammals and macroparasites: from evolutionary ecology to management. Berlin: Springer; 2006.Google Scholar
- Alvarez N, McKey D, Kjelberg F, Hossaert-McKey M. Phylogeography and historical biogeography of obligate specific mutualism. In: Morand S, Krasnov BR, editors. The biogeography of host-parasite interactions. New York: Oxford University Press; 2010. p. 31–9.Google Scholar
- Froeschke G, van der Mescht L, McGeoch M, Matthee S. Life history strategy influences parasite responses to habitat fragmentation. Int J Parasitol. 2013;45:1109–18.View ArticleGoogle Scholar
- Daszak P, Cunningham AA. Emerging infectious diseases of wildlife - threats to biodiversity and human health. Science. 2000;287:443–9.View ArticlePubMedGoogle Scholar
- Bitam I, de la Cruz K-D, Parola P, Whiting MF, Raoult D. Flea and flea-borne diseases. Int J Infect Dis. 2010;14:667–76.View ArticleGoogle Scholar
- Pérez-Rodríguez A, Fernández-González S, De le Hera I. 2013. Finding the appropriate variables to model the distribution of vector-borne parasites with different environmental preferences: climate is not enough. Glob Change Biol. 2013;19:3245–53.Google Scholar
- Poretta D, Mastrantonio V, Amendolia S, Gaiarsa S, Epis S, Genchi C, Bandi C, Otranto C, Urbanelli S. Effects of global changes on the climatic niche of the tick Ixodes ricinus inferred by species distribution modelling. Parasit Vectors. 2013;6:271.Google Scholar
- Estrada-Peña A, Estrada-Sánchez A, Estrada-Sánchez D, de la Fuente J. Assessing the effects of variables and background selection on the capture of the tick climate niche. Int J Health Geogr. 2013;12:43.View ArticlePubMedPubMed CentralGoogle Scholar
- Adjemian JCZ, Girvetz EH, Beckett L, Foley JE. Analysis of genetic algorithm for rule-Set production (GARP) modeling approach for predicting distributions of fleas implicated as vectors of plague, Yersinia pestis, in California. J Med Entomol. 2006;43:93–103.PubMedGoogle Scholar
- Rose H, Wall R. Modelling the impact of climate change on spatial patterns of disease risk: sheep blowfly strike by Lucilia sericata in great Britain. Int J Parasitol. 2011;41:739–46.View ArticlePubMedGoogle Scholar
- Wardhana AH, Cecchi G, Muharsini S, Cameron MM, Ready PD, Hall MJR. Environmental and phylogeographical determinants of the distribution of the Old world screwworm fly in Indonesia. Acta Trop. 2014;138:62–8.View ArticleGoogle Scholar
- Medvedev SG, Krasnov BR. Fleas: permanent satellites of small mammals. In: Morand S, Krasnov BR, Poulin R, editors. Micromammals and macroparasites: from evolutionary ecology to management. Berlin: Springer; 2006. p. 161–77.View ArticleGoogle Scholar
- Krasnov BR. Functional and evolutionary ecology of fleas: a model for ecological parasitology. New York: Cambridge University Press; 2008.View ArticleGoogle Scholar
- Ioff IG. 1941. Ecology of fleas in relevance to their medical importance. Pyatygorsk Publishers, Pyatygorsk, USSR (in Russian). In: Krasnov BR, Burdelova N, Shenbrot GI, Khokhlova IS. Annual cycles of four flea species in the central Negev desert. Med Vet Entomol. 2002;16:266–76.Google Scholar
- Combes C. Parasitism. The ecology and evolution of intimate interactions. Chicago: University of Chicago Press; 2001.Google Scholar
- Krasnov BR, Mouillot D, Shenbrot GI, Khokhlova IS, Vinarski MV, Korallo-Vinarskaya NP, Poulin R. Similarity in ectoparasite faunas of Palaearctic rodents as a function of host phylogenetic, geographic or environmental distances: which matters the most? Int J Parasitol. 2010;40:807–17.View ArticlePubMedGoogle Scholar
- Krasnov BR, Pilosof S, Shenbrot GI, Khokhlova IS. 2013. Spatial variation in the phylogenetic structure of flea assemblages across geographic ranges of small mammalian hosts in the palearctic. Int J Parasitol. 2013;43:763–70.View ArticlePubMedGoogle Scholar
- Krasnov BR, Pilosof S, Shenbrot GI, Khokhlova IS, Degen AA. Phylogenetic structure of host spectra in Palaearctic fleas: stability versus spatial variation in widespread, generalist species. Parasitology. 2014;141:181–91.View ArticlePubMedGoogle Scholar
- Krasnov BR, Shenbrot GI, Khokhlova IS, Stanko M, Morand S, Mouillot D. Assembly rules of ectoparasite communities across scales: combining patterns of abiotic factors, host composition, geographic space, phylogeny and traits. Ecography. 2015;38:184–97.View ArticleGoogle Scholar
- Krasnov BR, Sarfati M, Arakelyan MS, Khokhlova IS, Burdelova NV, Degen AA. Host-specificity and foraging efficiency in blood-sucking parasite: feeding patterns of the flea Parapulex chephrenis on two species of rodents. Parasitol Res. 2003;90:393–9.View ArticlePubMedGoogle Scholar
- Shenbrot G, Krasnov B, Lu L. Geographical range size and host specificity in ectoparasites: a case study with Amphipsylla fleas and rodent hosts. J Biogeogr. 2007;34:1679–90.View ArticleGoogle Scholar
- Segerman J. Siphonaptera of southern Africa. Handbook for the identification of fleas. Publications of the south African institute for medical research No. 57. Johannesburg, South Africa: South African Institute for Medical Research; 1995.Google Scholar
- van der Mescht L, Matthee S, Matthee CA. Comparative phylogeography between two generalist flea species reveal a complex interaction between parasite life history and host vicariance: Parasite-host intimacy matters. BMC Evol Biol. 2015;15:105.View ArticlePubMedPubMed CentralGoogle Scholar
- van der Mescht L, Matthee S, Matthee CA. A genetic perspective on the taxonomy and evolution of the medically important flea, Dinopsyllus ellobius (Siphonaptera: Dinopsyllinae), and the resurrection of Dinopsyllus abaris. Biol J Linn Soc. 2015;116:541–57.View ArticleGoogle Scholar
- Stockwell DRB, Peterson AT. Effects of sample size on accuracy of species distribution models. Ecol Model. 2002;148:1–13.View ArticleGoogle Scholar
- Wisz MS, Hijmans RJ, Li J, Peterson AT, Graham CH, Guisan A. NCEAS Predicting Species Distributions Working Group. Effects of sample size on the performance of species distribution models. Divers Distrib. 2008;14:763–73.View ArticleGoogle Scholar
- Robe JL, Machado S, Bolzan AR, dos Santos JPJ, Valer FB, dos Santos AP, Blauth ML, Gottschalk MS. Comparative ecological niche modeling and evolutionary ecology of Neotropical mycophagous Drosophilidae (Diptera) species. Stud Neotropical Fauna Environ. 2014;49:79–94.Google Scholar
- de Meillon B, Davis D, Hardy F. Plaque in southern Africa, vol.1. The Siphonaptera (excluding Ischnopsyllidae). Pretoria, South Africa: Government Printer; 1961.Google Scholar
- Austin M. 2007. Species distribution models and ecological theory: A critical assessment and some possible new approaches. Ecol Model. 2007;200:1–19.View ArticleGoogle Scholar
- Elith J, Phillips SJ, Hastie T, Dudik M, Chee YE, Yates CJ. A statistical explanation MaxEnt for ecologists. Divers Distrib. 2011;17:43–57.View ArticleGoogle Scholar
- Merow C, Smith MJ, Silander JA. A practical guide to MaxEnt for modelling species’ distributions: what it does, and why inputs and settings matter. Ecography. 2013;36:1058–69.View ArticleGoogle Scholar
- Stark HE. Population dynamics of adult fleas (Siphonaptera) on hosts and in nests of the California vole. J Med Entomol. 2002;39:818–24.View ArticlePubMedGoogle Scholar
- Krasnov BR, Khokhlova IS, Fielden LJ, Burdelova NV. The effect of temperature and humidity on the survival of pre-imaginal stages of two flea species (Siphonaptera: Pulicicdae). J Med Entomol. 2001;38:629–37.View ArticlePubMedGoogle Scholar
- Krasnov BR, Khokhlova IS, Fielden LJ, Burdelova NV. The effect of substrate on survival and development of two species of desert fleas (Siphonaptera: Pulicidae). Parasite. 2002;9:135–42.View ArticlePubMedGoogle Scholar
- Krasnov BR, Khokhlova IS, Fielden LJ, Burdelova NV. Time of survival under starvation in two flea species (Siphonaptera: Pulicidae) at different air temperatures and relative humidities. J Vector Ecol. 2002;27:70–81.PubMedGoogle Scholar
- Reisen WK. Landscape epidemiology of vector-borne diseases. Annu Rev Entomol. 2010;55:461–83.View ArticlePubMedGoogle Scholar
- Pettorelli N, Ryan S, Mueller T, Bunnefeld N, Jedrzejewska B, Lima M, Kausrud K. The Normalized Difference Vegetation Index (NDVI): unforeseen successes in animal ecology. Clim Res. 2011;46:15–27.View ArticleGoogle Scholar
- Andreo V, Lima M, Provensal C, Priotto J, Polop J. Population dynamics of two rodent species in agro-ecosystems of central Argentina: intra-specific competition, land-use, and climate effects. Pop Ecol. 2009;51:297–306.View ArticleGoogle Scholar
- Andreo V, Provensal C, Scavuzzo M, Lamfri M, Polop J. Environmental factors and population fluctuations of Akodon azarae (Muridae: Sigmontinae) in central Argentina. Austral Ecol. 2009;34:132–42.View ArticleGoogle Scholar
- Myneni RB, Hoffman S, Knyazikhin Y, Privette JL, Glassy J, Tian Y, Running SW. Global products of vegetation leaf area and fraction absorbed PAR from year one of MODIS data. Remote Sens Environ. 2002;83:214–31.View ArticleGoogle Scholar
- Hengl T, de Jesus JM, MacMillan RA, Batjes NH, Heuvelink GBM, Ribeiro E, Samuel-Rosa A, Kempen B, Leenaars JGB, Walsh MG, Gonzalez MR. SoilGrids1km—Global soil information based on automated mapping. PLoS One. 2014;9:105992.Google Scholar
- Cook BD. Rabbit burrows as environments for European rabbit fleas, Spilopsyllus cuniculi (Dale), in arid South Australia. Aust J Zool. 1990;38:317–25.View ArticleGoogle Scholar
- du Plessis A, Kerley GIH, deo Winter PE. Refuge microclimates of rodents: a surface nesting Otomys unisulacatus and a burrowing Parotomys brantsii. Acta Theriol. 1992;37:351–8.View ArticleGoogle Scholar
- Osacar-Jimenez JJ, Lucientes-Curdi J, Calvete-Margolles C. Abiotic factors influencing the ecology of wild rabbit fleas in north-eastern Spain. Med Vet Entomol. 2001;15:157–66.View ArticlePubMedGoogle Scholar
- Shenbrot G, Kransov B, Khokhlova I, Demidova T, Fielden L. Habitat-dependence difference in architecture and microclimate of the burrows of Sundevall’s jird (Merione crassus) (Rodentia: Gerbillinae) in the Negev Desert Israel. J Arid Environ. 2002;51:265–79.View ArticleGoogle Scholar
- Estrada-Peña A, Estrada-Sánchez A, de la Fuente J. A global set of Fourier-transformed remotely sensed covariates for the description of abiotic niche in epidemiological studies of tick vector species. Parasit Vectors. 2014;7:302.View ArticlePubMedPubMed CentralGoogle Scholar
- Cumming GS, van Vuuren DP. Will climate change affect ectoparasite species ranges? Glob Ecol Biogeogr. 2006;15:486–97.View ArticleGoogle Scholar
- R Development Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. ISBN 3-900051-07-0, URL http://www.R-project.org/ (2014).
- ArcGIS Desktop: Release 10. Redlands, CA: Environmental Systems Research Institute; 2012.Google Scholar
- Franklin J. Mapping species distributions: spatial inference and prediction. New York: Cambridge University Press; 2009.Google Scholar
- Phillips SJ, Anderson RP, Schapire RE. Maximum entropy modelling of species geographic distributions. Ecol Model. 2006;190:231–59.View ArticleGoogle Scholar
- Phillips SJ, Dudik M. Modeling of species distributions with MaxEnt: new extensions and comprehensive evaluation. Ecography. 2008;31:161–75.View ArticleGoogle Scholar
- Peterson AT, Papeş M, Eaton M. Transferability and model evaluation in ecological niche modelling: a comparison of GARP and MaxEnt. Ecography. 2007;30:550–60.View ArticleGoogle Scholar
- Allouche O, Tsoar A, Kadmon R. Assessing the accuracy of species distribution models: prevalence, kappa, and the true skill statistic (TSS). J Appl Ecol. 2006;43:1223–32.View ArticleGoogle Scholar
- Wisz MS, Pottier J, Kissling WD, Pellissier L, Lenoir J, Damgaard CF, Dormann CF, Forchhammer MC, Grytnes J-A, Guisan A, Heikkinen RK, Høye TT, Kühn I, Luoto M, Maiorano L, Nilsson M-C, Normand S, Öckinger E, Schmidt NM, Termansen M, Timmerman A, Wardle DA, Aastrup P, Svenning J-C. The role of biotic interactions in shaping distributions and realised assemblages of species: implications for species distribution modelling. Biol Rev. 2013;88:15–30.Google Scholar
- Vashchenok VS. 1988. Fleas—vectors of pathogens causing diseases in humans and animals. Nuaka, in Russian. In: Krasnov BR, Shenbrot GI, Khokhlova IS, Stanko M, Morand S, Mouillot D. 2015. Assembly rules of ectoparasite communities across scales: combining patterns of abiotic factors, host composition, geographic space, phylogeny and traits. Ecography. 2015;38:184–97.View ArticleGoogle Scholar
- Singleton G, Krebs CJ, Davis S, Chambers L, Brown P. Reproductive changes in fluctuating house Mouse populations in south eastern Australia. Proc R Soc Lond B. 2001;268:1741–8.View ArticleGoogle Scholar
- Brown JH, Ernest SKM. Rain and rodents: complex dynamics of desert consumers. Bioscience. 2002;52:979–87.View ArticleGoogle Scholar
- Jaksic FM, Lima M. Myths and facts on ratadas: bamboo blooms, rainfall peaks and rodent outbreaks in South America. Austral Ecol. 2003;28:237–51.View ArticleGoogle Scholar
- Meserve PL, Kelt DA, Milstead WB, Gutiérrez JR. Thirteen years of shifting top-down and bottom-up control. Bioscience. 2003;53:633–46.View ArticleGoogle Scholar
- Zhang Z, Pech R, Davis S, Shi D, Wan X, Zhong W. Extrinsic and intrinsic factors determine the eruptive dynamics of Brandt's voles Microtus brandti in Inner Mongolia, China. Oikos. 2003;100:299–310.View ArticleGoogle Scholar
- Rosenzweig ML. Habitat selection experiments with a pair of coexisting heteromyid rodents species. Ecology. 1973;54:111–7.View ArticleGoogle Scholar
- Norrdahl K, Korpimäki E. Does predation risk constrain maturation in cyclic vole populations? Oikos. 1995;72:263–72.View ArticleGoogle Scholar
- Spencer R-J, Cavanough VC, Baxter GS, Kennedy MS. Adult free zones in small mammal populations: response of Australian native rodents to reduced cover. Austral Ecol. 2005;30:876–84.View ArticleGoogle Scholar
- MacArthur RH. Fluctuations of animal populations and a measure of community stability. Ecology. 1955;36:533–6.View ArticleGoogle Scholar
- MacArthur RH. Geographical ecology. Princeton: Princeton University Press; 1972.Google Scholar
- Brown JH. Macroecology. Chicago: University of Chicago Press; 1995.Google Scholar
- Rothschild M, Clay T. Fleas, Flukes and Cuckoos: A Study of Bird Parasites. 3rd edition London: Collins; 1952. In: Krasnov BR, editor. Functional and Evolutionary Ecology of Fleas: A Model for Ecological Parasitology. 3rd ed. New York: Cambridge University Press; 2008.Google Scholar
- Hopkins GHE. 1957. Host associations of Siphonaptera. In: Krasnov BR. Functional and Evolutionary Ecology of Fleas: A Model for Ecological Parasitology. New York: Cambridge University Press; 2008.Google Scholar
- Krasnov BR, Shenbrot GI, Medvedev SG, Vatschenok VS, Khokhlova IS. Host-habitat relations as an important determinant of spatial distribution of flea assemblages (Siphonaptera) on rodents in the Negev Desert. Parasitology. 1997;114:159–73.View ArticlePubMedGoogle Scholar
- Krasnov BR, Shenbrot GI, Medvedev SG, Khokhlova IS, Vatschenok VS. Habitat dependence of a parasite-host relationship: Flea (Siphonaptera) assemblages in two gerbil species of the Negev Desert. J Med Entomol. 1998;35:303–13.View ArticlePubMedGoogle Scholar
- Burda H, Šumbera R, Begall S. Microclimate in burrows of subterranean rodents—Revisited. In: Begall S, Burda H, Schleich CE, editors. Subterranean rodents: News from underground. Berlin: Springer; 2007.Google Scholar
- Hudson BD. Soil organic matter and available water capacity. J Soil Water Conserv. 1994;49:189.Google Scholar
- Laundre JW, Reynolds TD. Effects of soil structure on burrow characteristics of five small mammal species. Great Basin Nat. 1993;53:358–66.Google Scholar
- Coreau A, Martin J-L. Multi-scale study of bird species distribution and of their response to vegetation change: a Mediterranean example. Landsc Ecol. 2007;22:747–64.View ArticleGoogle Scholar