- Open Access
Sarcoptes scabiei infestation does not alter the stability of ectoparasite communities
© The Author(s). 2016
- Received: 15 March 2016
- Accepted: 21 June 2016
- Published: 1 July 2016
The host represents a heterogeneous ecosystem where multiple parasite species co-occur and interact with each other for space and resources. Although these interactions may rule the features of an infracommunity and may shape the infracommunity response to external perturbations, the resilience of ectoparasite communities to new infestations remains poorly explored.
We analysed the composition of the ectoparasite communities found on 214 individual Iberian ibexes (Capra pyrenaica) inhabiting the Sierra Nevada Natural Space, southern Spain. Using classification and regression trees, we explored how the presence of Sarcoptes scabiei (a highly contagious mite), the off-host environment and the host sex govern the prevalence and abundance of lice and ticks. Null model analysis was applied to assess the impact of S. scabiei on the structure of the ectoparasite communities.
Our results suggest that S. scabiei infestation acts in tandem with off-host environment and host sex to define the prevalence and abundance of lice and ticks. We also provided evidence for differences in species co-occurrence only at the early stages of S. scabiei infestation. Regarding species diversity, we recorded that ectoparasite communities in scabietic ibexes reached a high richness faster than those in healthy individuals.
Even though we show that ectoparasite burden is correlated with S. scabiei infestation, off-host environment and host sex, the species response to S. scabiei infestation and climate seem to be highly variable and influenced by ectoparasite life-history traits. Ectoparasite communities also appear resilient to perturbations which is in agreement with what was previously reported for endoparasites. Future refinement of sample collection and the incorporation of ecological and epidemiological-related variables may allow us to establish causal effects and deepen the knowledge about the mechanisms and consequences of ectoparasite interactions.
- Capra pyrenaica
- Community ecology
- Ecosystem engineer
- Sarcoptes scabiei
The host represents a heterogeneous ecosystem encompassing a wide range of linked microhabitats (Table 1). From the host’s skin surface to the host’s inner organs, each microhabitat has particular characteristics that influence the prevalence and abundance of parasite populations . However, the tissue tropism and the fundamental niche of a parasite species are not only influenced by microhabitat suitability. The host immune and defensive responses (top-down regulation) as well as the within-host competition (exploitation, apparent and interference) for space and resources (bottom-up regulation) are believed to act, synergistically or not, in the maintenance of parasite specificity and parasite burdens . Gaining a broader understanding of parasite interactions within host is of paramount importance not only because such interactions drive several features of an infracommunity but also because may shape disease dynamics .
Although the vast majority of studies reporting interspecific interactions has been achieved at laboratory level and focused on helminths and microparasites , recent experiments have started to explore the interactions between endo and ectoparasites [5, 6]. However, thus far the mechanisms of parasite co-infection and the resilience of infracommunities to external perturbations remain poorly understood [7, 8]. For instance, infracommunity interactions between different ectoparasite taxa have only recently begun to receive attention. By inducing external perturbations through the application of acaricides in the ectoparasite community of a wild mammal, recent research showed that ectoparasites interactions affect host susceptibility to other ectoparasites as well as the distribution of ectoparasites among hosts [9, 10]. It also suggested that species with differing life-history traits are affected by perturbations in different ways [9, 10]. Yet how ectoparasite communities respond to infestations caused by highly contagious mites, namely those infestations able to alter hosts’ phenotype, remains unclear. Several parasite species may be involved, directly or indirectly, in ecosystem engineering processes either causing immune and structural changes in their hosts or altering host’s traits . Parasite-mediated changes can range from antagonistic to facilitating the establishment and growth of existing or subsequent parasites . Consequently, the infestation by one virulent parasite may have important implications in the dynamics of within-host ectoparasite communities.
This study complies with the Spanish and the Andalusian laws regarding bioethics and animal welfare. The Sierra Nevada National Park approved this study.
The Sierra Nevada Natural Space (SNNS; Latitude 36°00'–37°10'N, Longitude 2°34'–3°40'W) holds the largest and best known population of Iberian ibex in Andalusia . This area is characterised by a Mediterranean subarctic climate experiencing seasonal and altitudinal gradients of temperature and precipitation . The average monthly temperature ranges from -5 °C in February to 17 °C in July and may vary between 12 and 16 °C below 1,500 m and 0 °C above 3,000 m. Annual average precipitation is approximately 600 mm . The snow is generally present between December and May; vegetation growth mainly occurs during summer months (June–August).
The data were collected from 214 ibexes (131 males and 83 females) shot-harvested between 2002 and 2008 on a monthly basis as part of a sarcoptic mange control program devoted to manage ibex density and the spread of mange in the SNNS. Iberian ibexes of both sexes and all ages were shot-harvested whenever mange lesions were detected. Additionally, some ibexes were selectively harvested in the context of a population management program devoted to ensure the equilibrium of the Iberian ibex population. Sex was assigned by visual inspection while age determination was performed through horn-segment counts . Ibexes were grouped into three age classes: kids (0 ≤ age ≤ 1 year, n = 20), yearlings (1 < age ≤ 2 years, n = 42) and adults (> 2 years, n = 152); kids were excluded from the analysis. The severity of S. scabiei infestation was assessed by measuring the surface of scabietic lesions and by digesting five skin fragments (6.25 cm2 each ). The digestion procedure was carried out overnight using a 5 % potassium hydroxide (KOH) at 40 °C. Each fragment was digested separately. The resulting products were re-suspended and analysed using a stereomicroscope for mite counts. The mean of the five counts was used as a proxy of mite load. The animals were categorised as healthy (ibexes without skin lesions), mildly infested (skin surface affected ≤ 50 %) and severely infested (skin surface affected > 50 %) . Ectoparasites (lice and ticks) were systematically removed and counted by three observers during 15 min. Therefore, the number of lice and ticks represents an abundance index. Lice were treated with lactic acid and mounted with DePex medium. Ticks were fixed in 70 % ethanol. Lice and ticks were identified to the species level using the available descriptions [32–34].
Definition of terms
Concept coined by  that refers to organisms that directly or indirectly alter the environment in which they occur. Their action may modify or create new habitats through their own physical structures (autogenic engineers) or through their activities (allogenic engineers).
Differences in tissue preferences or resources used by a species in response to interspecific competition . It is an indicator that different parasite species are interacting.
Includes all the infrapopulations that colonise a single host or an organ and interact with each other for space and resources .
Small-scale environment where an organism or an assemblage of organisms naturally occur and interact, both directly and indirectly, with the biotic and abiotic elements.
Positive or negative variation of a particular infrapopulation size when facing the presence of another species . It is an indicator that different parasite species are interacting.
Effects of an infection in host fitness, i.e. severity of the disease caused by a particular organism (damage-related concept ). It could also express the transmissibility of an organism, i.e. its capacity to grow and proliferate within a host .
Variables selected by the models and its respective measure units. In this study, we defined four calendar-based seasons
Group of ectoparasites
Average temperature for the month of sampling (°C)
Average NDVI for the month of sampling (mm)
Average summer temperature for the year of sampling (°C)
Average winter temperature for the year of sampling (°C)
Average autumn temperature for the year of sampling (°C)
Average spring NDVI for the year of sampling (mm)
Average autumn NDVI for the year of sampling (mm)
Lice and Ticks
Degree of S. scabiei infestation
Prevalence and mean abundance of ectoparasite species
The prevalence and mean abundance  of ectoparasite (lice and ticks) species were estimated using the software Quantitative Parasitology 3.0 . The confidence interval (CI) for prevalence was calculated through Sterne method , while bootstrap (BCa, 2,000 replications) were used to estimate the CI for the mean abundance.
Classification and regression trees
The contribution of S. scabiei infestation, off-host environment and host sex to the prevalence and mean abundance of lice and ticks (numerical response) were explored through classification (prevalence) and regression (abundance) trees (CART ). CART’s main advantage relies on its flexibility to handle interactions and nonlinearities among variables, predictive power and easy interpretation. The two main issues in constructing a reliable and informative tree model are to find good data splits and to avoid data over-fitting. The first can be addressed through the determination of information gain or node impurity measures (e.g. entropy, Gini index of diversity or misclassification error) whereas model over-fitting is reduced by pruning the tree . In our analysis, the information gain was applied to select the best split and the complexity parameter was used in order to prune the tree and represent the data as simple and interpretable as possible. The prediction error rate in cross-validation procedure was used to assess the model reliability. CART models were fitted using the ‘rpart’ library  and plotted using the ‘rpart.plot’ library , R statistical software version 3.1.3 .
Null model analysis was used to explore whether pairwise associations among ectoparasite (lice and ticks) species occurred by chance or were influenced by the severity of S. scabiei infestation (functional response); these were also used to study the patterns of species diversity (richness) among healthy and scabietic hosts.
To assess patterns in species co-occurrence, three datasets, each representing a degree of S. scabiei infestation, were organised as r × c matrices where each row represents a species of ectoparasite and each column an individual ibex. The cells in the matrix denote the presence (1) or the absence (0) of a particular species of ectoparasite in a particular ibex. We assume individual ibexes as replicates to unveil repeated patterns of species co-occurrence. Co-occurrence was assessed by computing the number of checkerboard species pairs, the default co-occurrence index (C-score) and the number of species combinations. The number of checkerboard species was originally defined as an indicator of species competition . This statistic identifies the number of pairs that do not occur together in any site, i.e. two or more ectoparasite species (lice and ticks) have mutually exclusive distributions among the sampled ibexes. The C-score is a measure of “checkerboardedness” that expands the concept of “checkerboard distributions”. C-score quantifies the average number of checkerboard "units" for each species pair . In our study, C-score determines the randomness of the distribution of n ectoparasite species (lice and ticks) through the sampled ibexes. Higher values of C-score are associated to a more segregated matrix/distribution, whereas a lower C-score is associated to a more aggregated matrix/distribution. The number of species combinations calculates the number of unique combinations represented in each ibex. In a competitively structured community, the following assumptions must be fulfilled: (i) the C-score should be significantly greater than expected by chance (O > E), and/or (ii) the number of checkerboard pairs of species should be larger than expected by chance (O > E), and/or (iii) the number of species combinations should be lower than expected by chance (O < E). In our analysis, we used a fixed-equiprobable algorithm in which we kept the observed rows fixed, i.e. the number of occurrences of each species in the null communities and in the original dataset is the same, and the columns were treated as equally likely, i.e. all ibexes share the same suitability to be invaded by an ectoparasite. The indices were calculated for each presence/absence matrix and compared with the expected indices computed for 5,000 randomised communities through Monte Carlo procedures.
Species richness is the most straightforward and easy-to-interpret measure of species diversity. In this study, one single dataset was created to develop species accumulation curves for each stage of S. scabiei infestation. The dataset was organised as a matrix where each row represents an ectoparasite species and each column the abundance of ectoparasites by stage of S. scabiei infestation. Species richness was assessed through a species accumulation curve in which the number of parasite species is measured as a function of sampling effort, e.g. the number of parasite individuals identified . The shape and slope of the curve towards its asymptote are influenced by the prevalence and abundance of each parasite species and may vary considerably between infracommunities. The analyses were performed using the EcoSim 7.72 software .
Prevalence and mean abundance of ectoparasite (lice and ticks) species collected from 194 Iberian ibexes (yearling and adults) of both sexes. The confidence interval (CI 95 %) for prevalence and mean abundance are also presented
Total no. examined
Prevalence (95 % CI)
Mean abundance (95 % CI)
Our classification and regression tree analyses indicated that host sex may interact with S. scabiei infestation and environmental drivers in shaping the prevalence and abundance of one louse species, L. stenopsis, which mainly occurred in infested female hosts. We did not, however, detect a clear influence of the host sex on tick burdens (Fig. 2).
Co-occurrence analysis considering healthy ibexes and two stages of S. scabiei infestation/mange severity (mildly and severely infested hosts). The significance between observed (O) and expected by chance (E) values of the C-score, number of checkerboards and number of species combinations is presented for presence/absence matrices of ectoparasite communities of 214 Iberian ibexes
Mean of simulated indices
Variance of simulated indices
P-value (O ≤ E)
P-value (O ≥ E)
Number of checkerboards
Number of species combinations
Co-occurring ectoparasite species interact with each other and with on-host and off-host environment. Despite such interactions may govern all features of an infracommunity, their synergistic or antagonistic effects have seldom been assessed. We suggest that the combination of on-host and off-host features is necessary to understand the dynamics of ectoparasite communities in the wild.
Beyond the effects of S. scabiei infestation, our results suggest that the prevalence and abundance of lice and ticks follow a pronounced seasonal pattern and vary with off-host environment and host sex. We recorded slightly higher abundances of lice at colder and dryer periods. This trend was previously reported in lice populations occurring in sheep , feral ponies  and donkeys . For ticks, we reported that D. marginatus and Haemaphysalis sp. were more prevalent and abundant under average low temperatures, which coincided with the winter/early spring periods. The genus Rhipicephalus recorded the highest prevalence and abundance among the species identified. The increasing activity of this genus seems to be correlated with humid and warm winter/spring seasons which contrasts with the results reported in other studies (northern Iran ; northern Greece ). We believe that host senility and immunocompromised state caused by S. scabiei infestation coupled with diet impoverishment noticeable in the winter contribute to higher burdens of parasitism during this period. Collectively, our results corroborate empirical evidence about the role of climate in the regulation of ticks’ prevalence and abundance. However, we realise that each tick species has its own phenology which poses complications for comparative studies and for disentangling the effects of life-history traits, S. scabiei infestation and environmental conditions on the features of ticks’ communities.
We reported a weak but female-biased ectoparasitism in the prevalence L. stenopsis. For the others ectoparasite (lice and ticks) species, no patterns were observed. The most likely reasons for this result may be related to gender differences in the social structure, phenology and prevalence of sarcoptic mange. Female ibexes are generally more gregarious which should increase the number and duration of body-to-body contacts and therefore the parasite prevalence, abundance and spread . Almost one-third of females were sampled during the periparturient period. At this time, females are particularly susceptible to infestations due to immunosuppression. Further, the prevalence and severity of mange is significantly higher in males than in females  which may influence the burdens of ectoparasites (lice and ticks) among host sexes. The immune response of Iberian ibexes to S. scabiei is not influenced by host age  therefore, this variable was not included in the analysis. Additionally, our data also does not support that differences in vulnerability to parasitism are related to the age classes considered (yearlings and adults).
Regarding the expected changes that S. scabiei infestation may promote in ectoparasite community structure and composition, we reported differences in species co-occurrence only at early stages of S. scabiei infestation. This fact may indicate that infracommunities are resistant and resilient to external perturbation because after an early change of ectoparasite communities in response to S. scabiei infestation we detected a quick return to the initial condition, i.e. we recorded no significant differences in parasite co-occurrence patterns between healthy and severely infested ibexes. The stability of within-host communities in wild mammals was already experimentally demonstrated for endoparasites . Knowles et al. , by reducing the burdens of nematode infestation through anthelminthic administration, reported a surprisingly stability of non-target parasite communities to perturbation. Here, we observed a similar response of ectoparasite communities when exposed to a natural and external perturbation such as the infestation by a highly contagious species. In addition, we recorded that scabietic ibexes reached higher parasite richness faster than healthy ibexes, i.e. the asymptote of species accumulation curves was reached first in scabietic ibexes. We conclude that S. scabiei may influence the diversity of infracommunities through alteration of host energy needs, thermoregulation and spatial behaviour. Sarcoptes scabiei infestation causes an increase of temperature in the affected skin due to hyperemia which makes infested hosts more prone to detection by ectoparasites. Note that one way of tick host-finding is by sensing heat loss from hosts’ body. Further, the disruption of homeostasis by S. scabiei made the movement of ibexes increasingly difficult, i.e. infested ibexes tend to exhibit abnormal movements and to move less than healthy ones . Consequently, the home range of infested ibexes tends to decrease and this may create more chances for body-to-body contact and, therefore, for parasite transmission and spread among hosts . Additionally, scabietic ibexes with lower home ranges use specific areas more intensively; this poses an increased risk of new infestations and several opportunities for some ectoparasites to attach. Using radio-collars it may be possible to confirm how the parasite diversity varies in relation to host home range and to which extent the severity of S. scabiei infestation affects the spatial behaviour of infested hosts.
Our study broadens our understanding of the dynamics of infracommunities and constitutes the first insight into the mechanisms underpinning co-infestation in a wild Iberian host experiencing a highly contagious mite infestation. More precisely, we explored how ecosystem engineering and allogenic processes carried out by a highly contagious parasite may impact the features of ectoparasite community. In our study, the numerical response of ectoparasite communities to external perturbations was highly variable and our results suggested that the presence of one highly contagious parasite may not negatively influence the presence and abundance of lice and ticks. Future refinement of sample collection aiming to control for uneven sampling effort and the potential bias caused by host and parasite phenology would serve to increase the reliability of our results and to assess the role of climatic fluctuations on the critical stages of tick’s development. The incorporation of further drivers, either ecologically related variables such as the habitat structure and altitude or epidemiologically related variables such as host population density and individual home range, will allow to make better predictions about the influence of external factors on both parasite counts and community structure. The quantification of attachment sites should be included in future studies to test the spatial competition between ectoparasite species within host. We advocate that manipulative experiments need to be run in order to clarify cause-effect relationships. Such detailed information will allow us to use advanced analytical tools (e.g. mechanistic models) to analyse, in greater detail, the mechanisms of parasite interactions and competitive exclusion at host level.
CART, Classification and Regression Trees; CI, Confidence Interval; NDVI, Normalized Difference Vegetation Index; MODIS, Moderate-Resolution Imaging Spectroradiometer; SNNS, Sierra Nevada Natural Space
The research was supported by the project CGL2012-40043-C02-01 (MEC, Spanish Goverment). The authors’ research activities are partially supported by the Andalusian Goverment (RNM-118 group). João Carvalho was supported by a PhD Grant (SFRH/BD/98387/2013) co-financed by Fundação para a Ciência e a Tecnologia (FCT), European Social Fund (ESF) and Ministério da Educação e Ciência (MEC) National Funds. Emmanuel Serrano was supported by the postdoctoral program (SFRH/BPD/96637/2013) of the FCT, Portugal. The authors acknowledge the support of the Sierra Nevada Natural Space staff for providing biological samples used in this study. We are grateful to João Gonçalo who kindly designed our Fig. 4. We would like to thank University of Aveiro (Department of Biology) and FCT/MEC for the financial support to CESAM RU (UID/AMB/50017) through national funds and co-financed by the FEDER, within the PT2020 Partnership Agreement.
João Carvalho was supported by a PhD Grant (SFRH/BD/98387/2013) co-financed by the Fundação para a Ciência e a Tecnologia (FCT), European Social Fund (ESF) and Ministério da Educação e Ciência (MEC) national funds. Emmanuel Serrano was supported by the postdoctoral program (SFRH/BPD/96637/2013) of the FCT, Portugal.
Availability of data and material
The datasets are currently being used in other studies and will be shared as soon as possible.
Designed the study: JC, ES, NP, CF and JMP. Performed the ibex sampling: JEG, MAH, SO and JMP. Analysed the data: JC, ES, NP, CF and JMP. Wrote the paper: JC, ES, NP, CF and JMP. All authors read and approved the final version of the manuscript.
The authors declare that they have no competing interests.
Ethics approval and consent to participate
This study complies with the Spanish and the Andalusian laws regarding bioethics and animal welfare. The Sierra Nevada National Park approved this study.
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.
- Rynkiewicz EC, Pedersen AB, Fenton A. An ecosystem approach to understanding and managing within-host parasite community dynamics. Trends Parasitol. 2015;31:212–21.View ArticlePubMedGoogle Scholar
- Mideo N. Parasite adaptations to within-host competition. Trends Parasitol. 2009;25:261–8.View ArticlePubMedGoogle Scholar
- Johnson PTJ, de Roode JC, Fenton A. Why infectious disease research needs community ecology. Science 2015;349:1259504. doi:10.1126/science.1259504.
- Telfer S, Lambin X, Birtles R, Beldomenico P, Burthe S, Paterson S, Begon M. Species interactions in a parasite community drive infection risk in a wildlife population. Science. 2010;330:243–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Bandilla M, Valtonen ET, Suomalainen L-R, Aphalo PJ, Hakalahti T. A link between ectoparasite infection and susceptibility to bacterial disease in rainbow trout. Int J Parasitol. 2006;36:987–91.View ArticlePubMedGoogle Scholar
- Knowles SCL, Fenton A, Petchey OL, Jones TR, Barber R, Pedersen AB. Stability of within-host parasite communities in a wild mammal system. Proc R Soc B. 2013;280:20130598.View ArticlePubMedPubMed CentralGoogle Scholar
- Pedersen AB, Fenton A. Emphasizing the ecology in parasite community ecology. Trends Ecol Evol. 2007;22:133–9.View ArticlePubMedGoogle Scholar
- Pedersen AB, Antonovics J. Anthelmintic treatment alters the parasite community in a wild mouse host. Biol Lett. 2013;9:20130205.View ArticlePubMedPubMed CentralGoogle Scholar
- Lutermann H, Fagir DM, Bennett NC. Complex interactions within the ectoparasite community of the eastern rock sengi (Elephantulus myurus). Int J Parasitol Parasites Wildl. 2015;4:148–58.View ArticlePubMedPubMed CentralGoogle Scholar
- Hoffmann S, Horak IG, Bennett NC, Lutermann H. Evidence for interspecific interactions in the ectoparasite infracommunity of a wild mammal. Parasit Vectors. 2016;9:58.View ArticlePubMedPubMed CentralGoogle Scholar
- Thomas F, Poulin R, Meeüs T, Guégan J-F, Renaud F. Parasites and ecosystem engineering: what roles could they play? Oikos. 1999;84:167–71.View ArticleGoogle Scholar
- Labaude S, Rigaud T, Cézilly F. Host manipulation in the face of environmental changes: Ecological consequences. Int J Parasitol Parasites Wildl. 2015;4:442–51.View ArticlePubMedPubMed CentralGoogle Scholar
- Turchetto S, Obber F, Permunian R, Vendrami S, Lorenzetto M, Ferré N, et al. Spatial and temporal explorative analysis of sarcoptic mange in Alpine chamois (Rupicapra r. rupicapra). Hystrix. 2014;25:25–30.Google Scholar
- Pérez JM, Ruiz-Martinez I, Granados JE, Soriguer RC, Fandos P. The dynamics of sarcoptic mange in the ibex population of Sierra Nevada in Spain - Influence of climatic factors. J Wild Res. 1997;2:86–9.Google Scholar
- León-Vizcaíno L, de Ybáñez MR R, Cubero MJ, Ortíz JM, Espinosa J, Pérez L, et al. Sarcoptic mange in Spanish ibex from Spain. J Wildl Dis. 1999;35:647–59.View ArticlePubMedGoogle Scholar
- Pence DB, Ueckerman E. Sarcoptic mange in wildlife. Rev Sci Tech Off Int Epiz. 2002;21:385–98.View ArticleGoogle Scholar
- Wall R, Shearer D. Veterinary Ectoparasites. Biology, Pathology & Control. 2nd ed. Oxford: Blackwell Science; 2001.View ArticleGoogle Scholar
- Ellse L, Burden FA, Wall R. Seasonal infestation of donkeys by lice: Phenology, risk factors and management. Vet Parasitol. 2014;203:303–9.View ArticlePubMedGoogle Scholar
- Estrada-Peña A, de la Fuente J. The ecology of ticks and epidemiology of tick-borne viral diseases. Antiviral Res. 2014;108:104–28.View ArticlePubMedGoogle Scholar
- Krasnov BR, Bordes F, Khokhlova IS, Morand S. Gender-biased parasitism in small mammals: patterns, mechanisms, consequences. Mammalia. 2012;76:1–13.View ArticleGoogle Scholar
- Kiffner C, Stanko M, Morand S, Khokhlova IS, Shenbrot GI, Laudisoit A, et al. Variable effects of host characteristics on species richness of flea infracommunities in rodents from three continents. Parasitol Res. 2014;113:2777–88.Google Scholar
- Navarro-Gonzalez N, Verheyden H, Hoste H, Cargnelutti B, Lourtet B, Merlet J, et al. Diet quality and immunocompetence influence parasite load of roe deer in a fragmented landscape. Eur J Wildlife Res. 2011;57:639–45.View ArticleGoogle Scholar
- Boyer N, Réale D, Marmet J, Pisanu B, Chapuis J-L. Personality, space use and tick load in an introduced population of Siberian chipmunks Tamias sibiricus. J Anim Ecol. 2010;79:538–47.View ArticlePubMedGoogle Scholar
- Moore SL, Wilson K. Parasites as a viability cost of sexual selection in natural populations of mammals. Science. 2002;297:2015–8.View ArticlePubMedGoogle Scholar
- Granados JE, Pérez JE, Márquez FJ, Serrano E, Soriguer RC, Fandos P. La cabra montés (Capra pyrenaica, Schinz 1838). Galemys. 2001;13:3–37.Google Scholar
- Poulin R. Interactions between species and the structure of helminth communities. Parasitology. 2001;122:3–11.View ArticleGoogle Scholar
- Pérez JM, Granados JE, Soriguer RC, Fandos P, Márquez FJ, Crampe JP. Distribution, status and conservation problems of the Spanish Ibex, Capra pyrenaica (Mammalia: Artiodactyla). Mammal Rev. 2002;32:26–39.View ArticleGoogle Scholar
- Kottek M, Grieser J, Beck C, Rudolf B, Rubel F. World Map of the Köppen-Geiger climate classification updated. Meteorol Z. 2006;15:259–63.View ArticleGoogle Scholar
- Pérez-Luque AJ, Bonet FJ, Pérez-Pérez R, Aspizua R, Lorite J, Zamora R. Sinfonevada: Dataset of floristic diversity in Sierra Nevada forests (SE Spain). PhytoKeys. 2014;35:1–15.View ArticlePubMedGoogle Scholar
- Fandos P. La cabra montés Capra pyrenaica en el Parque Natural de las Sierras de Cazorla, Segura y Las Villas. Madrid: ICONA-CSIC; 1991.Google Scholar
- Pérez JM, Granados JE, Sarasa M, Serrano E. Usefulness of estimated surface area of damaged skin as a proxy of mite load in the monitoring of sarcoptic mange in free-ranging populations of Iberian wild goat, Capra pyrenaica. Vet Parasitol. 2011;176:258–64.View ArticlePubMedGoogle Scholar
- Rodríguez F, Jiménez A, Martín-Mateo MP. Primeras citas de malófagos parásitos de Capra pyrenaica hispanica. Nouv Rev Entomol. 1980;10:363–71.Google Scholar
- Pajot FX. Les poux (Insecta, Anoplura) de la région afrotropicale. Paris: IRD Éditions; 2000.Google Scholar
- Habela M, Peña J, Corchero E, Sevilla RG. Garrapatas y hemoparásitos transmitidos de interes veterinario en España. Manual práctico para su identificación. 1st ed. Madrid: Schering Plough Animal Health; 2000.Google Scholar
- Calvete C, Estrada R, Lucientes J, Estrada A. Ectoparasite ticks and chewing lice of red-legged partridge, Alectoris rufa, in Spain. Med Vet Entomol. 2003;17:33–7.View ArticlePubMedGoogle Scholar
- Pettorelli N. The Normalized Difference Vegetation Index. 1st ed. Oxford: Oxford University Press; 2013.View ArticleGoogle Scholar
- 37.Gu Y, Hunt E, Wardlow B, Basara JB, Brown JF, Verdin JP. Evaluation of MODIS NDVI and NDWI for vegetation drought monitoring using Oklahoma Mesonet soil moisture data. Geophys Res Lett. 2008;35. doi:10.1029/2008GL035772.
- Carvalho J, Granados JE, López Olvera JR, Pérez JM, Fandos P, Soriguer RC, Velarde R, Fonseca C, Pettorelli N, Serrano E. Sarcoptic mange breaks up bottom-up regulation of body condition in a large herbivore population. Parasit Vectors. 2015. doi:10.1186/s13071-015-1188-4.Google Scholar
- Zar JH. Biostatistical analysis. 4th ed. New Jersey: Prentice Hall; 1999.Google Scholar
- Bush AO, Lafferty KD, Lotz JM, Shostak AW. Parasitology meets ecology on its own terms: Margolis et al. revisited. J Parasitol. 1997;83:575–83.View ArticlePubMedGoogle Scholar
- Reiczigel J, Lang Z, Rózsa L, Tóthmérész B. Properties of crowding indices and statistical tools to analyze crowding data. J Parasitol. 2005;91:245–52.View ArticlePubMedGoogle Scholar
- Reiczigel J. Confidence intervals for the binomial parameter: some new considerations. Stat Med. 2003;22:611–21.View ArticlePubMedGoogle Scholar
- Breiman L, Friedman JH, Olshen RA, Stone CJ. Classification and regression trees. 1st ed. Wadsworth Statistics/Probability: Belmont; 1984.Google Scholar
- Williams G. Data Mining with Rattle and R. The Art of Excavating Data for Knowledge Discovery. 1st ed. New York: Springer Verlag; 2011.Google Scholar
- Therneau T, Atkinson B, Ripley B. Package “rpart”: Recursive partitioning and regression trees. 2013. http://cran.r-project.org/web/packages/rpart/index.html. Accessed 10 Nov 2015.
- Milborrow S. Package “rpart.plot”: Plot rpart Models. An Enhanced Version of plot.rpart. 2012. http://cran.r-project.org/web/packages/rpart.plot/index.html. Accessed 12 Nov 2015.
- R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing. 2014. http://www.R-project.org/. Accessed 30 Oct 2015.
- Diamond JM. Assembly of species communities. In: Cody ML, Diamond JM, editors. Ecology and evolution of communities. Harvard University: Belknap; 1995. p. 342–444.Google Scholar
- Stone L, Roberts A. The checkerboard score and species distributions. Oecologia. 1990;85:74–9.View ArticleGoogle Scholar
- Dove ADM, Cribb TH. Species accumulation curves and their applications in parasite ecology. Trends Parasitol. 2006;22:568–74.View ArticlePubMedGoogle Scholar
- Gotelli NJ, Entsminger GL. Ecosim: Null Models Software for Ecology, Version 7.72. 2001; Acquired Intelligence Inc, & Kesey-Bear.Google Scholar
- MacArthur R, Levins R. The limiting similarity, convergence, and divergence of coexisting species. Am Nat. 1967;101:377–85.View ArticleGoogle Scholar
- Samuel WM, Pybus MJ, Kocan AA. Parasitic Diseases of Wild Mammals. 2nd ed. Iowa: Iowa State University Press; 2001.View ArticleGoogle Scholar
- James PJ, Moon RD, Brown DR. Seasonal dynamics and variation among sheep in densities of the sheep biting louse. Bovicola Ovis Int J Parasitol. 1998;28:283–92.View ArticlePubMedGoogle Scholar
- Gawler R, Coles GC, Stafford KA. Prevalence and distribution of the horse louse, Werneckiella equi equi, on hides collected at a horse abattoir in south-west England. Vet Rec. 2005;157:419–20.View ArticlePubMedGoogle Scholar
- Noori NV, Rahbari S, Bokaei S. The seasonal activity of Rhipicephalus bursa in Cattle in Amol (Northern Iran). World Appl Sci J. 2012;18:590–3.Google Scholar
- Papazahaiadou MG, Papadopoulos EG, Himonas CA. Seasonal activity of ixodid ticks on goats in northern Greece. Vet Rec. 1995;136:585–8.View ArticlePubMedGoogle Scholar
- Altizer S, Nunn CL, Thrall PH, Gittleman JL, Antonovics J, Cunningham AA, et al. Social organization and parasite risk in mammals: Integrating theory and empirical studies. Annu Rev Ecol Evol Syst. 2003;34:517–47.View ArticleGoogle Scholar
- López-Olvera JR, Serrano E, Armenteros A, Pérez JM, Fandos P, Carvalho J, et al. Sex-biased severity of sarcoptic mange at the same biological cost in a sexually dimorphic ungulate. Parasit Vectors. 2015;8:583.View ArticlePubMedPubMed CentralGoogle Scholar
- Sarasa M, Rambozzi L, Rossi L, Meneguz PG, Serrano E, Granados J-E, et al. Sarcoptes scabiei: Specific immune response to sarcoptic mange in the Iberian ibex Capra pyrenaica depends on previous exposure and sex. Exp Parasitol. 2010;124:265–71.View ArticlePubMedGoogle Scholar
- Bordes F, Morand S, Kelt DA, Van Vuren DH. Home range and parasite diversity in mammals. Am Nat. 2009;173:467–74.View ArticlePubMedGoogle Scholar
- Jones CG, Lawton JH, Shachak M. Organisms as ecosystem engineers. Oikos. 1994;69:373–86.View ArticleGoogle Scholar
- Alizon S, Hurford A, Mideo N, Baalen V. Virulence evolution and the trade-off hypothesis: history, current state of affairs and the future. J Evol Biol. 2009;22:245–59.View ArticlePubMedGoogle Scholar
- Casadevall A, Pirofski L. Host-pathogen interactions: the attributes of virulence. J Infect Dis. 2001;184:337–44.View ArticlePubMedGoogle Scholar