Abstract
There is much aggregation of parasites in nature but the factors influencing ectoparasite variation in mammals and how this affects populations are unclear due to costly compensatory behaviours performed, such as grooming. Bat research is lacking due to their sensitivity and protection. This study uses body mass, forearm length and ectoparasite count data collected over decades at low-frequency during the hibernation period of the Greater Horseshoe Bat, to elucidate the different factors affecting average relative mite density (RMD) in males and females of high and low reproductive condition. On average, females were larger than males, carrying more mites, with females of high reproductive condition carrying significantly more than males and non-reproductive females, suggesting lowered immunocompetence or a grooming/reproduction trade-off. Using the Akaike’s Information Criterion and significance of Confidence Intervals, Ordinal Logistic Regression models, show that differing factors are likely to determine RMD for male and female Greater Horseshoe bats with possible implications for their conservation.
Key words: Ectoparasite, variation, reproductive, condition, trade-off
Introduction
Parasites are known to be heterogenetic in mammals (Wilson et al. 2001) but information about the variation in magnitude and probability of infestation within bat populations is lacking (Lucan 2006; Frank et al. 2015; Orlova et al. 2015).Ectoparasite infestation is likely to depend on exposure level and immunocompetence (Hart 1992; Luguterah and Lawer 2015) and females are often more heavily infected than males (Zahn and Rupp 2004; Patterson, Dick and Dittma 2008; Presley and Willig 2008) due to fundamental life-history differences such as sexual size-dimorphism or age differences (Stebbings 1966), and behaviour such as aggregation by reproductive females not exhibited by males (Williams and Findley 1979; Ransome 1991). Reproductive condition and related hormones can also have an effect on overall bat parasite load (Christe, Arlettaz and Vogel 2000).
Bat are rare and highly protected in the UK, (Joint Nature Conservation Committee 2014; The Bat Conservation Trust 2016) carrying a variety of ectoparasites including mites (Sachanowicz, Kristofik and Ciechanowski 2014; Frank et al. 2015) which need contact between bats for transmission (Reckhart and Kerth 2009) The understanding of bats is of ecological importance as they contribute to essential ecosystem services such as insect control and pollination (Kunz et al. 2011; Ghanem and Voigt 2012; Wanger 2014), and this study aims to investigate variation in average parasitic mite load per individual (relative mite density) in the Greater Horseshoe bat, (Rhinolophus ferrumequinum) using data collected over 34 years from 39 locations, via low-disturbance methods, during autumn and winter months.
In the UK, populations of R. ferrumequinum are confined to southwest England (The Bat Conservation Trust 2016). This species is particularly sensitive to site disturbance (Hutson, Mickleburgh and Racey 2001) but can be handled during torpor (e.g. Park, Jones and Ransome 2000) in hibernacula over autumn and winter Similarly to other bat species females are larger and likely have higher relative mite density (RMD) than males.
The Body Condition Index (BCI), a ratio calculated by dividing body mass by forearm length, has been used to investigate ectoparasite load (Speakman and Racey 1986; Lucan 2006; Sharifi et al. 2013; Postawa and Szubert-Kruszyriska 2014) and preliminary analyses, using scatterplots with linear regression lines and the Spearman’s rank correlation coefficients for the data of this study, demonstrated that there is no significant correlation between RMD and body condition for R. ferrumequinum males or females (Appendix 1: Table 1, Fig. 1), supporting similar findings in other bat species (Zahn and Rupp 2004; Lucan 2006; Sharifi et al. 2013, Postawa and Furman 2014;Postawa and Szubert-Kruszyriska 2014)
Body mass, forearm length, and RMD were compared for high and low reproductive condition males and females (indicated by visibility of false nipples in females and testes in males) and regression models were created to explain and predict RMD variation for the sample population as a whole, as well as for male and female data separately, to reveal potential differences in factors determining RMD between males and females.
First, greater body mass, forearm length and RMD was predicted for high reproductive condition (RC) females than for both females of low RC and all males (Sharifi et al. 2013). This is hypothesised to be due to the dynamics between parasite and host that take place within summer colonial roosts of adult females and some juveniles, but which adult males attend rarely and only until females give birth (Piksa et al. 2014). The higher density individuals within a roosting colony can have increases probabilities of horizontal and vertical ectoparasite transmission (Patterson, Dick and Dittmar 2008; Presley and Willig 2008).
RC in males and females is associated with age and body mass (Ward et al., 2014) and parasitic mites may have adapted their reproductive cycle to that of their host (Christe, Arlettaz and Vogel 2000; Lourenco and Palmeirim 2008). High RC (possibly pregnant) females may have a greater RMD than low RC, possibly juvenile, females during the autumn and winter, until the time of birth (summer) when parasites are thought to move onto male and female juveniles in the roost (Christe, Arlettaz and Vogel 2000). Therefore the second prediction of this study is that males of lower RC would have a higher RMD, but lower body mass and forearm length than both higher RC males and low RC females.
Finally, due to the possibility of different underlying causal mechanisms (Ward et al. 2014) influencing their parasite load (Lucan 2006), female RMD variation is predicted to be more dependent on RC than body mass or forearm length, while male RMD variation is predicted to be more dependent on body mass than on RC or forearm length.
Materials and Methods
a. Data Collection
R. ferrumequinum were captured while in torpor, weighed (g), sexed, aged by rating the ossification of joints (Kunz 1988) from 0-4, forearm length measured (mm) and number of wing mites and tail mites counted. The visibility of false nipples was recorded for females, the absence of which indicating non-reproductive condition (Dietz and Dietz 2007). For males, degree of testicular visibility was recorded, and larger or smaller testicle size can be an indicator of higher or lower spermatogenesis, respectively, during autumn and winter (Haarsma 2008). Bats were released by placement back on the cave wall in a suitable position. Hibernacula temperature, outside temperature and humidity were recorded at least once for most of the sites, except four which contributed very little data in total.
Data was collected between 1976 and 2010 from approximately 39 locations in Southwest England at varying frequencies (Wiltshire Bat Group 2015). Hibernacula at each site were usually visited, between November and April but in total, the majority of data was collected from 7 main locations with most data being collected between 2002 and 2010. Thus, spatial and population dependence is expected of the data. This low-frequency sampling by volunteering members of the Wiltshire bat group, and handling of bats during torpor (e.g. Park, Jones and Ransome 2000) under the supervision of a licenced bat handler, likely limited disruption.
b. Data Selection and Manipulation
Due to the large number of data entries for which two or more variables were unrecorded, gaps in the data were filled in using the mean of the recorded values (assuming no trend in the data in order to prevent bias), resulting in 1293 males and 1236 female data points for each recorded variable. This enabled statistical analyses without eliminating large parts of useful data.
New variables were systematically created from the raw recorded data. Male testes visibility, originally ranked from 0-4, was categorised into males with more regressed testes (visibility 0, 1 and half of 2) and males with more pronounced testes (visibility 3, 4 and the other half of 2). These categories, along with visibility of false nipples in females, were then given a rank of either 0, representing low RC (males with more regressed testes and females without visible false nipples), or 1, representing high RC (males with more pronounced testes and females with visible false nipples).
RMD was calculated, as the number of tail mites added to wing mites per bat, and since these data contained a high frequency of zero values, they were ranked into 8 ordinal groups (rank 1 = 0-2 mites/bat, rank 2 = 3-5 mites/bat, rank 3 = 6-8 mites/bat, rank 4 = 9-11 mites/bat, rank 5 = 12-14 mites/bat, rank 6 = 15-17 mites/bat, rank 7 = 18-20 mites/bat and rank 8 = 21 or more mites/bat) in order to avoid problems caused by high frequencies of zero values during statistical analyses.
c. Data Analysis
Effects of body mass, forearm length and RC on RMD were investigated using the statistical analysis software R (R Core Team, 2015), with all tests carried out at a 95% confidence interval. Before doing analyses, the data was explored, assessing outliers, distribution normality and homoscedasticity. Six extreme high and low outliers for body mass and forearm length data were removed as these were possibly caused by measurement error. Extreme outliers for RMD were all larger than the average and were not removed as individuals in natural animal populations can experience extreme infestation (e.g. Boulinier, Ives and Danchin 1996; Pierce et al. 2014; Frank et al. 2015).
Since the data for RMD was greatly right skewed, Kruskall-Wallis and Mann-Whitney U tests were done to investigate the differences in ranked RMD between males and females, and between RC groups (Kowalski, 1972). Body mass and forearm length data were found to be nearly normally distributed with similar homogeneity of variances, leading to the use of a one-way ANOVA and two-sample unpaired Student’s t-Tests to examine differences in the two variables between males and females, and between the RC groups.
Ordinal Logistic Regression (OLR) models were fitted using the Proportional Odds Regression function and MASS package in R (Venables and Ripley 2002), to explore the causal relationship between the ordinal dependent variable, RMD and the covariates (predictor variables) body mass, forearm length and RC on variation in RMD. Nine OLR models using the combined male and female data, male-only and female-only data were run and the Aike Information Criteria (AIC), Confidence Intervals (CI) and Odds Ratios were used to select an appropriate model that could explain the variation in RMD for each set of data.
To validate the models, evaluation of the Proportional Odds Assumption was carried out by running a series of binary logistic regressions and assessing the equality of the difference in coefficients of the logits for every outcome of the predictor variable, at each level of the dependent variable. Finally, the predicted probabilities for the new models, of a given dependent variable level being the case, at every given predictor variable level was plotted using the packages reshape2 (Wickham 2007) and ggplot2 (Wickham 2009) and compared to graphs of actual data in order to assess the models’ effectiveness.
Results
(For all statistical tests except OLR, results were significant unless otherwise stated)
a. Abiotic Factors and Bat Age
Outside temperature ranged from 0-18°C (8.5°C average). Hibernacula temperature ranges 0-18°C (9.1°C average), and humidity ranges 9.6-100% (86.1% average). Bat joint ossification ratings ranged 0-4 but the average was 3.
b. RMD Variation
Females had 1.14 times (Fig. 2) greater average RMD than males (Appendix 2, Table Set 2.a). Females of high RC had 1.063 times significantly greater RMD than females of low RC. Males of lower RC had 1.071 times greater average RMD than males of higher RC but this difference was not significant. High RC females had greater average RMD than all males while there was no significant difference between low RC females and males of either RC (Appendix 3,Table 3.b) , supporting the trend in Fig. 2.
c. Body Mass Variation
Female bats had 1.050 times (Fig.3) greater average body mass than males (Appendix 2, Table 2.b). High RC females had 1.063 times greater average body mass than low RC females, and males of both higher and lower RC. Higher RC males had 1.062 times significantly greater average body mass than lower RC males. There was no significant difference in body mass between low RC females and higher RC males but low RC females had significantly greater average body mass than lower RC males (Appendix 3, Table 4.b). These results support the trend in Fig. 3.
d. Forearm Length Variation
On average, female bats had 1.017 times (Fig. 4) longer forearms than males (Appendix2, Table 2.b). High RC females had 1.003 times greater average forearm length than low RC females and males of both higher and lower RC (as did females of lower RC). Higher RC males had 1.006 times greater forearm length than lower RC males (Appendix 3, Table 4.c). These results support the trend in Fig.4.
e. Ordinal Logistic Regression modelling
Nine OLR models were created with each predictor variable in turn, using combined, as well as male-only and female-only data (Appendix 4, Table 5). The CI was significant for models using 1. body mass and RC for combined data, 2. body mass for male-only data, and 3) RC and forearm length for female-only data. The AIC and Residual Deviance for the models that did not have significant CI was also relatively high. In the female-only data, forearm length Odds Ratio showed an inverse relationship with RMD compared to the Odds Ratios for the other significant variables. The predictor variables of these significant models were used to fit three OLR models, M1 (combined data), M2 (male) and M3 (female), that may predict the variation in RMD of R. ferrumequinum (Appendix 5, Table 6) for the sample population.
For model validation, the Proportional Odds Assumption was found to be violated for RC in M1 but not for Body Mass (Appendix 6, Table 7). It was not violated in M2 (Appendix 6, Table 8), and was violated for Forearm Length in M3 but not for RC (Appendix 6, Table 9). The predicted probabilities of a given level of the dependent variable being the case, at any given level of each predictor variable, were then plotted for each model and compared to graphs of actual data (Fig.5, 6, and 7). The models were found to effectively predict the change in RMD in relation to 1. Body Mass and RC for the combined data, 2. mass in the male data and 3.forearm length and reproductive condition in the female data .
Model 1 used Body Mass (g) and RC (0=lower, 1=higher) to explain ranked RMD (levels 1-8). This model had the highest AIC, with the effect of RC being found not significant. The Body Mass variable had an odds ratio of 1.06: the odds for a higher outcome in RMD is 1.06 times greater at any higher value of Body Mass compared to the next value down, given that all other variables stay constant (Appendix 5, Table 6). Model 2 used only Body Mass to explain ranked RMD in male bats. This was the model with the lowest AIC, with the same general outcome as for Model 1(Fig. 5.), but with a greater Odds Ratio of 1.09 (Appendix 5, Table 6). Model 3 used RC and Forearm Length (mm) to explain ranked RMD in female bats and also had a relatively low AIC, with the highest Odds Ratio for RC of 1.57 and the lowest odds ratio for forearm length at 0.84, meaning the odds for a higher outcome in RMD is 1.57 times greater at higher RC than at low RC, and approximately 1.2 times greater at a lower forearm length compared to the next value up, given that all other variables stay constant (Appendix 5, Table 6).
Discussion
Ectoparasite load can carry implications for reproduction and survivorship (Giorgi et al. 2001; Simon et al. 2004; Hawlena, Abramsky and Krasnow 2006). Direct impact may occur through resource consumption (Hawlena, Abramsky, Krasnov 2006) and imposing immunocompetence costs (Navarro-Gonzalez et al. 2011), as well as indirectly through pathogen transmission, and changes in habitat selection (Reckardt and Kerth 2007).
Grooming behaviour functions as defence against ectoparasite infestation but carries an energetic and time cost (Giorgi et al. 2001; Godinho et al. 2013) that may also impact the ability of an animal to rest, forage, defend its territory or mate, reproduce and care for young (Giorgi et al. 2001). This energetic cost may also vary according to differences in the cost that different parasite species draw from hosts (Godinho et al. 2013).
For this study of R. ferrumequinum during autumn and winter in southwest England, average RMD was greater in high RC females than any other group, while there was no significant RMD difference between low RC females and any males, supported by other studies (Christe et al. 2007; Encarnacao, Baulechner and Bekcer 2012; Sundari et al. 2012).
Grooming behaviour may be energetically and time restricted in reproductive females and inefficient in low RC individuals that are likely to be juveniles (McLean and Speakman 1997; Viljoen et al. 2011; Encarnacao, Baulechner and Bekcer 2012). Reproductive females may also have reduced immunocompetence due to hormonal changes (Christe et al. 2000) or to a resource allocation trade-off for reproduction (Dzal and Bringham 2013). In addition, reproductive females tend to live in close colonies in maternity roosts for part of the year, often with some juveniles, while non-reproductive females and males roost solitarily (Ransome 1991) leading to a greater probability of parasite transmission within the reproductive female population (Christe, Arlettaz & Vogel 2000; Patterson, Dick and Dittmar 2008).
Modelling found that high RC females are more likely to have a higher level of mite infestation than low RC females (and all males) but longer female forearms mean lower probability of high-level mite infestation and vice versa (Fig. 7). Greater body mass may mean lower probability of low-level mite infestation and higher probability of higher-level infestation for all bats (Fig. 5) as is the case for males separately (Fig. 6) but not for females separately. Possible reasons include higher mass providing more surface area and more resources for mites to inhabit (Viljoen et al. 2011) and less energy spent on costly grooming when larger individuals are able to carry the cost of infestation better than small individuals (Giorgi et al. 2001).
Differences in body mass, forearm length and RMD between the study groups were extremely small, as were the predicted increases or decreases in RMD, and there is a general need for future study of factors affecting on mite infection intensity (average mites per infected bat) and prevalence (percentage of infested bats).
In general, ectoparasite load does not affect overall bat health ( this study’s Appendix 1: Table 1, Fig. 1; Zahn and Rupp 2004; Lucan 2006; Sharifi et al. 2008; Postawa and Szubert-Kruszyriska 2014) but at times of more intense demand on resources, even a moderate parasite load may affect survival. Fluctuations in climate, lack of suitable roosts and decreased food availability such as those faced by bats in the UK (Ransome 1968; Hutson, Mickleburgh, Racey 2001; The Bat Conservation Trust, 2014) can affect ectoparasite prevalence due to temperature (Chen and Mullins 2008; Lourenco and Palmeirim 2008) and humidity changes (Moyer et al. 2002; Zhang et al. 2010), as well as body mass and immunocompetence (Navarro-Gonzalez et al. 2011; Dlugosz et al. 2014). Heterogeneity of infestation impact due to varying energetic and time demands may carry implications for bat population dynamics and conservation that need further investigation over the course of all seasons (Lourenco and Palmeirim 2008; Sharifi et al. 2008).
References
Boulinier T, Ivers AR and Danchin E (1996) Measuring aggregation of parasites at different host population levels. Parasitology 112:581-587
Chen BL and Mullens BA (2008) Temperature and humidity effects on off-host survival of the Northern fowl mite (Acari:Macronyssidae) and the chicken body louse (Phthiraptera:Menoponidae). J Econ Entomol 101:637-346
Christe P, Arlettaz R and Vogel P (2000) Variation in intensity of a parasitic mite (Spinturnix myoti) in relation to the reproductive cycle and immunocompetence of its bat host (Myotis myotis). Ecol Lett 3:207-212. DOI: 10.1046/j.1461-0248.2000.00142.x
Christe P, Glaizot O, Evanno G, Bruyndoncks N, Devevey G, Yannic G, Patthey P, Maeder A, Vogel P and Arlettaz R (2007) Host sex and ectoparasites choice: preferences for, and higher survival on female hosts. J Anim Ecol 76:703-710
Dietz C and Dietz I (2007) Age classification and assessment of reproductive condition in female Greater Horseshoe bats (Rhinolophus ferrumequinum, Chiroptera:Rhinolophidae). In: Dietz C (2007), Aspects of ecomorphology in the five European Horseshoe Bats (Chiroptera:Rhinolophidae), Dissertation, University of Tübingen
Dlugosz EM, de Bellocq JG, Khokhlova IS, Degen AA, Pinshow B and Krasnov BR (2014) Age at weaning, immunocompetence and ectoparasite performance in precocial desert rodent. J Exp Biol 217:3078-3038. doi: 10.1242/jeb.106005
Dzal YA and Bringham RM (2013) The tradeoff between torpor use and reproduction in little brown bats (Myotis licifugus). J Comp Physiol B 183:279-288
Frank R, Kuhn T, Werblow A, Liston A, Kochmann J, and Klimpel S (2015) Parasite diversity of European Myotis species with special emphasis on Myotis myotis (Micropchiroptera, Vespertilionidae) from a typical nursery roost. Parasite Vector 8:1-13. doi: 10.1186/s13071-015-0707-7
Giorgi MS, Arlettaz R, Christe P and Vogel P (2001) The energetic grooming costs imposed by a parasitic mite (Spinturnix myoti) upon its bat host (Myotis myotis). Proc. R. Soc. Lond. B 268:2071-2075. dio: 10.1098/rspb.2001.1686
Ghanem SJ, and Voigt CC (2012). Increasing Awareness of Ecosystem Services Provided by Bats. In:Brockman HJ, Roper TJ, Naguib M, Mitani JC and Simmons LW (eds) Advances in the Study of Behavior, vol. 44. Elsevier, USA, pp 279-302
Godhino LN, Cripps JK, Coulson G and Lumsden LF (2013) The effect of ectoparasites on the grooming behaviour of Gould’s Wattled bat (Chalinolobus Gouldii): An experimental study. Acta Chiropterol 15:463-472. doi:10.3163/150811013×679080
Haarsma AJ (2008) Manual for assessment of reproductive status, age and helwth in European Vespertilionid bats. Haarsma AJ, Holland http://vleermuizenvangen.nl/index.php/achtergrond-info/voorbeeld-examen Accessed 18 January 2016
Hart BL (1992) Behavioural adaptations to parasites: an ethological approach. J Parasitol 78:256-265
Hawlena H, Abramsky Z, Krasnov BR (2006) Ectoparasites and age-dependent survival in a desert rodent. Oecologia 148:30-39
Hutson AM, Mickleburgh SP and Racey PA (2001) Microchiropteran bats: global status survey and conservation action plan. IUCN / SSC Chiroptera Specialist Group. IUCN, Gland, Switzerland and Cambridge, U.K.
Joint Nature Conservation Committee (2014) EC Habitats Directive. JNCC. http://jncc.defra.gov.uk/page-1374. Accessed 10 January 2016
Kowalski CJ (1972) On the effects of non-normality on the distribution of the Sample Product-moment Correlation Coefficient. J R Stat Soc 21:1-12. doi:10.2307/2346598
Kunz TH (1988) Ecological and behavioural methods for the study of bats. Smithsonian Institution Press, Washington D.C.
Kunz TH, de Torrez EB, Bauer D, Lobova T and Fleming TH (2011) Ecosystem services provided by bats. Ann. N.Y. Acad. Sci. 1223:1-38 . doi: 10.1111/j.1749-6632.2011.06004.x
Lucan RK (2006) Relationships between the parasitic mite Spinturnix andegavinus (Acari: Spinturnicidae) and its bat host, Myotis daubentonii (Chiroptera: Vespertilionidae): seasonal, sex- and age-related variation in infestation and possible impact of the parasite on the host condition and roosting behaviour. Folia Parasitol 53:147-153
Lourenco S and Palmeirim JM (2008) Which factors regulate the reproduction of ectoparasites of temperate-zone cave-dwelling bats? Parasitol Res 104:127-134
Luguterah A and Lawer EA (2015) Effect of dietary guild (frugivory and insectivory) and other host characteristics on ectoparasite abundance (mite and nycteribiid) of chiropterans. Folia Parasit 62:1-7. doi: 10.14411/fp.2015.021
McLean JA and Speakman JR (1997) Non-nutritional maternal support in the brown long-eared bat. Anim Behav 54:1193-1204
Moyer BR; Drown DM and Clayton DH (2002) Low humidity reduced ectoparasite pressure: implications for host life history evolution. Oikos 97:223-228. doi: 10.1034/j.1600-0706.2002.970208.x
Navarro-Gonzalez N, Verheyden H, Hoste H, Cargnelutti B, Lourtet B, Merlet J, Daufresne T, Lavin S and Hewison AJM (2011) Diet quality and immunocompetence influence parasite load of roe dee in a fragmented landscape. Eur J Wildlife Res 57:639-645
Orlova MV, Zhigalina AV, Orlov OL, Kruskop SV and Bogdanove II (2015) Contribution to the ectoparasite fauna of rare and poor studied bat species of Southern Siberia. Biology Bull 42:254-259
Park KJ, Jones G and Ransome RD (2000) Torpor, arousal and activity of hibernating Greater Horseshoe Bats (Rhinolophus ferrumequinum) Funct Ecol 4:580-588
Patterson BD, Dick CW and Dittmar K (2008) Sex biases in parasitism of neotropical bats by bat flies (diptera:Streblidae). J Tropic Ecol 24:387-396. doi:10.1017/S0266467408005117
Pierce AA, de Roode JC, Altizer S and Bartel RA (2014) Extreme Heterogeneity in parasitism despite low population genetic structure among Monarch butterflies inhabiting the Hawaiian Islands. PLOS. doi: 10.1371/journal.pone.0100061
Piksa K, Gorz A, Nowak-Chmura M and Siuda K (2014) The patterns of seasonal activity of Ixodes vespertilionis (Acari:Ixodidae) on Rhinolophus hipposideros ins nursery colonies. Ticks and Tick-borne Diseases 5:69-74
Postawa T and Furman A (2014) Abundance patterns of ectoparasites infesting different populations of Miniopterus species in their contact zone in Asia Minor. Acta Chiropterl 16:387-39
Postawa T and Szubert-Kruszyriska A (2014) Is parasite load dependent on host aggregation size? The case of the greater mouse-eared bat Myotis Myotis (Mammalia: Chiroptera) and its parasitic mite Spinturnix myoti (Acari:Gamasida). Parisitol Res 113:1803-1811. doi:10.1007/s00436-014-3826-8
Presley SJ and Willig MR (2008) Intraspecific patterns of ectoparasite abundances on Paraguayan bats: effects of host sex and body size. J Tropic Ecol 24:75-83
R Core Team (2015). R: A language and environment for statistical computing. R Foundation for statistical computing, Vienna, Austria. https://www.R-project.org/
Ransome R D (1968) The distribution of the Greater horse-shoe bat, Rhinolophus ferrum-equinum, during hibernation, in relation to environmental factors. J Zool 154:77–112. doi:10.1111/j.1469-7998.1968.tb05040.x
Reckardt K and Kerth G (2007) Roost selection and roost switching of female Bechstein’s bat (Myotis bechsteinii) as a strategy of parasite avoidance. Oecologia 154:581-588 doi: 10.1007/s00442-007-0843-7
Ransome RD (1991) Greater horseshoe bat. In: Corbet GB and Harris S (eds), The handbook of British mammals, pp. 88-94, Blackwell, Oxford, UK.
Reckardt K and Kerth G (2009) Does the mode of transmission between hosts affect the host choice strategies of parasites? Implications from a field study on bat fly and wing mite infestation of Bechstein’s bats. Oikos 118: 183–190
Sachanowicz K, Kristofik J and Ciechanowski M (2014) Spinturnicid mites of bats in Albania – host spectrum and morphometrics as a tool of species separation. J Nat Hist 48:2661-2674. doi: 10.1080/00222933.2014.939729
Sharifi M, Fatema M, Najmeh T and Hossain J (2008) Variation in the ectoparasite load reflects life history traits in the lesser mouse-eared bat Myotis Blythii (Chiroptra: Vespertilionidae) in Western Iran. J Parasitol 94:622-625
Sharifi M, Taghinezhad N, Mozafari F and Vaissi S (2013) Variation in ectoparasite load in the Mehely’s horseshoe bat, Rhinolophus mehelui (Chiroptera:Rhinolophidae) in a nursery colony in western Iran. Acta Parasitol 58:180-184. doi: 10.2478/s11686-013-0122-1
Simon A, Thomas D, Blondel J, Perret P and Lambrechts MM (2004) Physiological ecology of mediterranean Blue Tits (Parus caeruleus L.): effects of ectoparasites (Protocalliphora spp.) and food abundance on metabolic capacity of nestlings. Physiol Biochem Zool 77:492-501
Speakman JR and Racey PA (1986) The influence of body condition on sexual development of male brown long-eared bats (Plectus auritus) in the wild. J Zool (Lond) 210:515-525
Stebbings RE (1966) A population study of bats of the genus Plecotus. J Zool 150: 53–75. doi:10.1111/j.1469-7998.1966.tb02998.x
Sundari AA, Bogdanowicz W, Varman DR, Marimuthu G and Rajan KE (2012 Ectoparasite Raymondia lobulata infestation in relation to the reproductive cycle of its host – the Greater False Vampire Bat Megaderma lyra. J Parasitol 98:60-62. doi:10.1645/GE-2833.1
The Bat Conservation Trust (2016) Greater Horseshoe Bat. BCT. http://www.bats.org.uk/pages/-greater_horseshoe_bat-800.html#legal status. Accessed 10 January 2016
The Bat Conservation Trust (2014) Conservation Research: BCT Science Strategy. http://www.bats.org.uk/pages/research.html Accessed 23 November 2015
Thompson RCA, Lymbery AJ, Smith A (2010) Parasites, emerging disease and wildlife conservation. Int J Parasitol 40:1163-1170
Venables WN and Ripley BD (2002) Modern Applied Statistics with S. 4th Edition. Springer, New York
Viljoen H, Bennett NC, Ueckermann EA, Lutermann H (2011) The Role of Host Traits, Season and Group Size on Parasite Burdens in a Cooperative Mammal. PLoS ONE 6: e27003. doi:10.1371/journal.pone.0027003
Wanger TC, Darras K, Bumrungsri S, Tscharntke T and Klein AM (2014) Bat pest control contributes to food security in Thailand. Biol Conserv 171:220-223
Ward HL, Ransome RD, Jones G, Rossiter SJ (2014) Determinants and Patterns of Reproductive Success in the Greater Horseshoe Bat during a Population Recovery. PLoS ONE 9: e87199. doi:10.1371/journal.pone.0087199
Wickham H (2007) Reshaping data with the reshape package. J Stat Softw, 21:1-20 http://www.jstatsoft.org/v21/i12/
Wickham H (2009) ggplot2: elegant graphics for data analysis. Spring, New York
Williams DF and Findley JS (1979) Sexual size dimorphisms in Vespertilionid bats. The American Midland Naturalist 102:113-126. doi:10.2307/2425072
Wilson K, Bjornstad ON, Dobson AP, Merler S, Poglaen G, et al. (2001) Heterogeneities in macroparasite infections: patterns and processes. In: Hudson PJ, Rizzoli A, Grenfell BT, Heesterbeek H, Dobson AP, eds. The ecology of wildlife diseases, pp6-44. Oxford: Oxford University Press
Wiltshire Bat Group (2015). Wiltshire Bat Group Web Page. Wiltshire and Swindon Biological Records Centre. http://www.wsbrc.org.uk/Links/LocalGroups/WBG/wideTemplate.aspx. Accessed 23 November 2015
Zahn A and Rupp D (2004) Ectoparasite load in European vespertilionid bats. J Zool 262:383-391
Zhang L, Parsons S, Daszak P, Wei L, Zhu G and Zhang S (2010) Variation in the abundance of ectoparasitic mites of flat-headed bats. J Mammal 91:136-143
Appendix 1.
Preliminary finding for the effects of an immunity indicator on R. ferrumequinum Relative Mite Density:
Appendix 2.
Results of statistical tests for differences between male and female R. ferrumequinum.
Appendix 3.
Results of statistical tests for differences within male and female R. ferrumequinum: between individuals with false nipples, without false nipples, and between individuals with larger testes and without larger testes.
Appendix 4.
Results of Ordinal Logistic Regression modelling used to select the models which best describe variation in RMD.
Appendix 5.
Results of Ordinal Logistic Regression modelling for the three selected models which best describe variation in RMD.
Appendix 6.
Predicted linear logits.