The Movement Ecology of the Straw-Colored Fruit Bat, Eidolon helvum, in Sub-Saharan Africa Assessed by Stable Isotope Ratios

Abstract

Flying foxes (Pteropodidae) are key seed dispersers on the African continent, yet their migratory behavior is largely unknown. Here, we studied the movement ecology of the straw-colored fruit bat, Eidolon helvum, and other fruit bats by analyzing stable isotope ratios in fur collected from museum specimens. In a triple-isotope approach based on samples of two ecologically similar non-migratory pteropodids, we first confirmed that a stable isotope approach is capable of delineating between geographically distinct locations in Sub-Saharan Africa. A discriminant function analysis assigned 84% of individuals correctly to their capture site. Further, we assessed how well hydrogen stable isotope ratios (δ²H) of fur keratin collected from non-migratory species (n = 191 individuals) records variation in δ²H of precipitation water in sub-Saharan Africa. Overall, we found positive, negative and no correlations within the six studied species. We then developed a reduced major axis regression equation based on individual data of non-migratory species to predict where potentially migratory E. helvum (n = 88) would come from based on their keratin δ²H. Across non-migratory species, δ²H of keratin and local water correlated positively. Based on the isoscape origin model, 22% of E. helvum were migratory, i.e. individuals had migrated over at least 250 km prior to their capture. Migratory individuals came from locations at a median distance of about 860 km from the collection site, four even from distances of at least 2,000 km. Ground-truthing of our isoscape origin model based on keratin δ²H of extant E. helvum (n = 76) supported a high predictive power of assigning the provenance of African flying foxes. Our study highlights that stable isotope ratios can be used to explain the migratory behavior of flying foxes, even on the isotopically relatively homogenous African continent, and with material collected by museums many decades or more than a century ago.

About the Authors:

Gonzalo Ossa

Affiliation: Leibniz Institute for Zoo and Wildlife Research, Berlin, Germany

Stephanie Kramer-Schadt

Affiliation: Leibniz Institute for Zoo and Wildlife Research, Berlin, Germany

Alison J. Peel

Affiliations Department of Veterinary Medicine, University of Cambridge, Cambridge, United Kingdom, Institute of Zoology, Zoological Society of London, Regent’s Park, London, United Kingdom

Anne K. Scharf

Affiliation: Leibniz Institute for Zoo and Wildlife Research, Berlin, Germany

Christian C. Voigt

* E-mail: [email protected]

Affiliations Leibniz Institute for Zoo and Wildlife Research, Berlin, Germany, Behavioral Biology, Freie Universität Berlin, Germany

Introduction

Flying foxes (family Pteropodidae) are central ecosystem service providers in the Old World [1], [2]. More than 280 plant species benefit from flying foxes either by pollination and/or seed dispersal. Ecosystem services offered by pteropodids result in about 450 economically valuable plant products like fruits or timber [1], [3], [4]. Flying foxes are essential for maintaining the genetic connectivity of plants among fragmented patches of rainforests and distant habitats due to their capability to move over long distances [5]–[7]. Pteropodids play also a prominent role as a reservoir for zoonotic diseases [8]–[11]. Thus, understanding the movement ecology of pteropodids is important, not only from a conservation point of view but also from an epidemiological and human health standpoint. In this study, we looked at the movement ecology of Eidolon helvum to better understand their potential role in the large-scale distribution of seeds across the African continent and to shed light on the spatial use of flying foxes in context to zoonotic diseases. Secondly, we tested for the first time whether variation in hydrogen stable isotope ratios of fur keratin can be used to assess the origin of flying foxes in Africa.

The primary habitat of E. helvum encompasses the tropical forests and savannahs in central Africa, from where it migrates following seasonal fluctuations in precipitation and associated changes in food abundance [12]. This species provides important services of pollination and promotes out-crossing in flowering plants [13]. For example, Daïnou and colleagues [14] suggested that migratory movements of E. helvum largely affect the geographic distribution of Milicia excelsa (Moraceae), one of the most important sources for timber in Africa [15]. Baker highlighted the ecological importance of E. helvum for West African ecosystems by cross-pollinating Ceiba pentandra (Malvaceae) [16]. Due to the large population sizes, E. helvum is potentially a key species for regenerating African forests [17]. Populations of E. helvum seem to follow the seasonal variations in fruit abundance [18], [19], i.e. migrating along a south – north axis in West Africa, from Ivory Coast to Niger. Other studies describe specific sedentary periods during a year when individuals aggregate in large colonies [20]–[22]. Most recent studies using satellite telemetry in three individuals at the Kasanka National Park in Zambia showed that E. helvum is capable of migrating several thousands of kilometers in a period of five months [5]. Here, we used a novel approach to the study of migration in E. helvum of sub-Saharan Africa by analyzing stable isotope ratios in fur collected from museum specimens and from live animals in their natural habitat.

Hydrogen stable isotope ratios have been increasingly used as endogenous markers to study the movement ecology of birds [23], [24] and bats [25]–[27]. The major advantages of this approach are that the movement ecology of animals can be assessed from a single capture event [28], [29], and that even historical specimens from museum collections can be studied. The stable isotope approach in the study of animal migration is based on the following two premises: (1) Stable hydrogen isotope ratios of rainfall vary geographically according to temperature [30]–[32]. This is because when gaseous water travels in the atmosphere along temperature gradients, water molecules labeled with heavy isotopes condensate faster than those labeled with light isotopes [33]. (2) Biologically inert body products, such as feather and fur, retain the stable isotope signature of the location where they grew [25], [26], [28], [29]. Recently, other elements such as carbon and nitrogen have been included in studies to improve the geographical resolution for predicting the provenance of animals [26], [34]. Environmental variation in carbon stable isotopes is usually related to geographic variation in plant species compositions and their specific photosynthetic metabolism of CO₂ fixation (C3, C4 and CAM). Nitrogen isotope ratios of ecosystems are known to vary with water availability [35], and, therefore, recent studies included nitrogen as well in triple-isotope approaches for geolocating animals, e.g. European bats [26] or Afrotropical birds [34].

In this study, we wanted to elucidate the migratory behavior of the straw-colored fruit bat, E. helvum, in sub-Saharan Africa based on fur material from museum specimens collected over the course of the last century. We organized our study in three parts. Firstly, we used a triple-isotope approach (carbon, nitrogen and hydrogen) to assess whether fur keratin varies geographically in non-migratory flying foxes in general, and whether geographically distinct populations can be separated based on stable isotope of these three elements. We then focused on hydrogen stable isotope ratios, because currently we lack a geographically complete isoscape for nitrogen and carbon on the African continent; a prerequisite for following a triple-isotope approach. Secondly, we tested within each and across non-migratory species whether hydrogen stable isotope ratios of fur keratin records variation of hydrogen stable isotope ratios extrapolated for precipitation water of the sites where individuals were collected for museums. Based on data of these non-migratory species, we developed a reduced major axis (RMA) regression model. Thirdly, using the RMA regression model, we assessed the likely migratory movements of E. helvum specimens in the period prior to their collection. This assessment was based on information of a) site of capture, b) hydrogen stable isotope ratios of fur keratin in museum specimens, and c) the extrapolated hydrogen isoscape of the African continent [31]. Lastly, we tested the quality of our isoscape origin model by applying it to data obtained from contemporary E. helvum individuals captured in two Tanzanian cities. We preferred to validate our approach with live animals over museum specimens, because we obtained better information about the time and exact location of capture for contemporary bats than for specimens of museum collections. Here, we show that stable isotope ratios can be used in the study of migratory vertebrates even in the relatively homogenous isoscape of the African continent. Further, we show that large portions of E. helvum populations are non-migratory in Sub-Saharan Africa, but that some individuals can travel as far as several thousand kilometers. This supports the potential of E. helvum for providing long-distance seed dispersal across the African continent.

Materials and Methods

Selection of Species and Collection of Fur Material from Museums

In our study, we focused on E. helvum (n = 88 individuals) which is assumed to migrate over large parts of the sub-Saharan African continent, and as a reference we selected six likely non-migratory pteropodid species: Rousettus aegyptiacus (n = 34 individuals), Lissonycteris angolensis (n = 17), Epomophorus wahlbergi (n = 23), Hypsignathus monstrosus (n = 30), Epomops franqueti (n = 30), and Epomophorus crypturus (n = 57). Non-migratory species were selected so to maximize sample size for each species, based on those available in the collections of the Natural History Museum in Paris (MNHN Paris), France, and of Berlin (MfN Berlin), Germany. Since past studies on pteropodids provided only very limited insights into the movement ecology of these species [36], we were not able to choose a priori bat species with documented non-migratory habits, except R. aegyptiacus and L. angolensis. These two species are known to roost in cave shelters during daytime and are considered to be strictly sedentary [37]. Daily foraging movements of R. aegyptiacus are usually limited to less than 50 km [7], yet one study suggests long-distance migration in a single individual of R. aegyptiacus [38]. Both species have a broad but patchy distribution in sub-Saharan Africa, probably related to the patchy distribution of cave roosts across the African continent [36]. E. wahlbergi and H. monstrosus perform only short distance movements of a few kilometers [39], [40]. We lack any information about the movement ecology of E. franqueti, a pteropodid with a wide distributional range from Senegal to Ethiopia [41]. E. crypturus is considered to travel “considerable distances” in search for food, yet more accurate descriptions of the exact distances and whether movements are performed seasonally are lacking [41]. The target species of this study, E. helvum, is classified as Near Threatened by the International Union for Conservation of Nature (IUCN) in response to the considerable decline of local populations due to hunting for food, use in traditional medicine and owing to the large-scale deforestation in many areas of its distributional range [42]. African populations of E. helvum are classified in the appendix II of the Convention on Migratory Species [43], according to their broad distribution in sub-Saharan Africa, and their ecological service function by pollination and seed dispersal [43]. Since there is currently no indication for a gender-specific migratory behavior in E. helvum, we did not differentiate between males and females in our analysis.

We obtained fur samples from alcohol preserved specimens from the two museums (Table 1) that came from collections conducted over a period of more than 100 years (1881–2011; Table S1). Allowance to collect fur samples from museum specimens was granted by the responsible mammal curator; Specimens remained in the corresponding collection during and after sample collection. We chose specimens for which information about the location and preferentially the year of collection was available. The locations of bats included in this study are given in the electronic supplementary material (Table S1, Figure S1). From each of the selected specimens we collected a small fur samples (about 2 mg) from the inter-scapular region of the bats’ back. Samples were placed in plastic tubes (0.5 ml), marked with the serial number of the respective museum, stored at ambient temperature and shipped to the stable isotope laboratory of the Leibniz Institute for Zoo and Wildlife Research (IZW). We determined the geographic coordinates and altitude for the collection localities using Google Earth [44].

[Figure omitted. See PDF.]

Table 1. Number of fur samples collected from museum specimens of seven species of Pteropodidae according to museum origin.

https://doi.org/10.1371/journal.pone.0045729.t001

Collection of Fur Material from Live Specimens for Model Validation

In addition to the samples from museum specimens, we also collected 76 fur samples of extant E. helvum from two colonies in Tanzania: 44 samples from a colony in Dar es Salaam (S 06.80134, E 039.28249) and 32 from Morogoro (S 06.80134, E 039.28249). Bats were captured using mistnets between 28^th of August and 29^th of September 2009 and fur samples were collected from the interscapular region of the dorsal fur, transferred to plastic vials and shipped to the IZW. Local residents reported that the E. helvum colonies were present year-round; yet it is unknown whether migratory individuals may join the non-migratory population at certain times. Migration has been reported to occur between October and December in East Africa [18]. Since we collected samples in August and September, we categorized the captured animals as non-migratory. Collection of samples from extant E. helvum was granted by the Tanzanian authorities (COSTECH permit: 2009-225-NA-2009-160). Prior to this, sample collection was approved by the Animal and Care and Ethics Committee of the IZW. All bats were released at the site of capture after fur collection.

Hydrogen Stable Isotopic Analysis

All analyses of hydrogen, carbon and nitrogen stable isotope ratios were carried out at the Stable Isotope Laboratory of the IZW (Berlin, Germany). Fur samples were cleaned with 2∶1 chloroform methanol solution for 24 hours to remove dirt and surface oils. The solvent was discarded, and the tubes with fur samples deposited in a drying oven (Heraeus T6, Thermo Electronic Corporation, Bremen, Germany) at 50°C for at least one day to eliminate the remainder of the solvent. We analyzed the hydrogen stable isotope ratio of samples in sequence with keratin reference material of known isotope ratio as described by [45]. This method controls for the exchange of hydrogen in keratin with ambient air. Thus, we refer to the non-exchangeable portion of hydrogen in fur keratin in the remainder of the text. The references used were powdered sheep fur from Sweden (SWE-SHE; 0.250±0.01 mg; mean ± one standard deviation) and from Spain (ESP-SHE; 0.249±0.01 mg) and goat fur from Tanzania (AFR-GOA; 0.255±0.01 mg). Each of these references was obtained from a single domestic animal (sheep or goat). The stable isotope ratios of the non–exchangeable hydrogen of the reference samples were −167.9±1‰ (SWE-SHE; (mean ± one standard deviation), −108.3±1‰ (ESP-SHE) and −66.3±0.9‰ (AFR-GOA) [26], [46].

From each of the fur and reference samples, we transferred a subsample of on average 0.274±0.01 mg into 3.5×5.0 mm silver foil capsules (OEA Laboratories, Callington, Cornwell, UK). Loaded capsules were formed as small cubes, transferred into a 96 port microtiter tray (69 bat fur subsamples +27 references of keratin hydrogen isotope) and stored in a drying oven at 50°C for at least 24 hours to allow subsamples and references to equilibrate with ambient air. Loaded capsules were then transferred to an autosampler carousel (Zero Blank autosampler – Costech Analytical Technologies Inc. Italy; [26], [47]).

Samples were pyrolized in a high-temperature elemental analyzer at 1350°C (HTO Elementaranalysator, HEKAtech GmbH, Wegberg, Germany). The resulting gases were carried by chemically pure helium (Linde, Luna, Germany) of H₂, N₂ and CO (N₂ and CO are sub products) were separated in a chromatographic column at 80°C. Hydrogen was carried to an isotope ratio mass spectrometer (Delta V Advantage IRMS, ThermoFisher Scientific, Bremen, Germany). All hydrogen stable isotopic ratios were expressed using the δ notation as units per mill (‰) and normalized on the Vienna Standard Mean Ocean Water (V-SMOW) standard scale. The precision of repeated measurements of laboratory standards was better than 1.2 ‰ (one Standard deviation) for hydrogen. In the remainder of the text, we will refer to the hydrogen stable isotope ratio of the non-exchangeable portion as δ²H_K and to the hydrogen stable isotope ratio of extrapolated precipitation water as δ²H_P. To obtain the mean annual δ²H_P from precipitation for each locality [31], we used the raster dataset in ArcGrid format provided by Bowen [46].

Carbon and Nitrogen Stable Isotope Analysis

For carbon and nitrogen stable isotope ratios, we used protein laboratory standards as references (Tyrosin: Tyr-b1; 0.35–0.45 mg; Leucine: Leu-b1; 0.25–0.30 mg) to monitor possible instrumental drift. From each laboratory standard and fur sample, we placed subsamples of 0.40±0.05 mg into 3.5×5.0 mm in tin foil capsules (OEA Laboratories, Callington, Cornwall, UK), and then transferred these packages to the autosampler carousal (MAS 200 autosampler, Thermo Electron Corporation, Bremen, Germany). Samples were oxidized in the reactor of the elemental analyzer (Flash EA 1112, Thermo Electron Corporation, Bremen, Germany). The resulting gaseous products were separated into CO₂ and N₂. The ratios of the stable isotopes ¹³C/¹²C and ¹⁵N/¹⁴N of the resulting organic combustion were measured using a Delta V Advantage isotope ratio mass spectrometer (Thermo Finnigan, Bremen, Germany) and were reported in δ notation in units per mill (‰) relative to the international Vienna-PeeDeeBelemnite (V-PDB) or ambient nitrogen (AIR-N₂) standard respectively. Laboratory standards were previously calibrated with standards (NBS 19, NBS 22, USGS 24 and L-SVEC for carbon; IAEA-N-1, IAEA-N-2 and IAEA-NO-3 for nitrogen). The δ¹³C and δ¹⁵N values of the laboratory standards were −24.0‰ and 4.4‰ respectively for Tyr-b1 and −30.3‰ and 11.0‰ respectively for Leu-b1. The precision of repeated measurements of laboratory standards was better than 0.18 ‰ (one standard deviation) for nitrogen and 0.14 ‰ (one SD) for carbon.

Statistical Analyses

All statistical analyses were performed with R 2.14.2 [48], and the significance values of all tests were set at 0.05. Values are presented as means ± one standard deviation unless stated otherwise. We tested for a normal distribution of data using the Shapiro-Wilk normality test and checked for heterogeneity of variances by plotting the data. In order to meet the assumptions of linear models, we excluded outliers with extreme values.

In our triple-isotope approach, we studied the explanatory power of δ²H, δ¹³C and δ¹⁵N to delineate between populations of non-migratory species from various geographical locations [26]. We focused on two species, E. crypturus and E. wahlbergi, that were similar in their ecology and were most likely to be non-migratory [49]. We only included data in our analysis that originated from locations where at least three individuals had been collected. To test for differences in δ²H_K, δ¹³C_, and δ¹⁵N among localities, we used a multivariate analysis of variance (MANOVA). Then, we used a discriminant function analysis (package MASS) to determine if site assignment was improved using the triple-isotope approach and by determining the proportion of correct assignments.

To assess the predictive power of hydrogen stable isotope ratios in fur keratin of non-migratory species for localizing species’ origins, we compared δ²H_K with predicted δ²H_P (obtained from http://www.waterisotopes.org). Since we usually lacked information about the Julian date of specimen collection, we used the annual average of δ²H_P for the collecting sites. We calculated Pearson’s product moment correlation (r) for the relationship between δ²H_K of non-migratory species and δ²H_P.

We calculated a reduced major axis (RMA) regression model based on the data of the non-migratory species to assess the movement ecology of E. helvum. For this, we linked δ²H_K of non-migratory specimens (n = 140) and annual mean δ²H_P values in Africa [45] (package lmodel2, [50]) to account for errors in both measures. Once the δ²H_k values were translated into δ²H_P, we modeled the uncertainty related to the analytical error using the δ²H_K standard deviations (σ_location) of the bat populations of known origin (non-migratory species) included in the regression model as error estimate and fitted a gamma distribution Γ to σ_location using maximum likelihood. We fitted Γ to 20 values of σ_location, having a minimum of 3 samples for each location. The shape of this function was k = 4.13 and the scale (θ) was 1.64.

To assess the migratory behavior of flying foxes based on stable isotope ratios measured in museum specimens, we modeled the origin of E. helvum using their individual δ²H_K values, as outlined in [27] and modified from [29] where each value of δ²H_K was associated to a location according to the African isoscape of hydrogen [45]. We reproduced the error of the distribution by a simulation of 5,000 possible values of δ²H_P for each individual’s δ²H_K drawing randomly from a normal distribution N(μ, σ), where μ corresponds to the exact value of δ²H_P obtained from the RMA regression, and σ was randomly generated from the gamma distribution Γ. From this, we created a probability density map using a probability density function (PDFs, package MASS) for each migratory individual (i.e. 88 maps for E. helvum). PDFs for each individual were standardized to range between 0 and 1. Using an ASCII remap file for each individual, we created the final map of origin for each individual. The reclassified raster surface was mapped using ArcGIS v. 9.3.1 with a 10′×10′ (arc minutes) resolution; distribution range of E. helvum was downloaded from the IUCN Red List of Threatened Species (www.iucnredlist.org). For each individual, we visually selected the maps showing the same origin probability (i.e. maps corresponding to individuals with similar δ²H_K from fur samples) and then calculated a mean probability map using the raster calculator from ArcGIS v. 9.3.1. For each mean probability map, we measured the distance of collection site of individuals to the nearest pixel with a minimum probability of origin of 0.7 inside the IUCN distribution range of the species, i.e. we considered probability values higher than 70% to be of sufficient discriminatory power.

We validated the isoscape origin model using stable hydrogen isotope ratios measured in fur keratin of contemporary E. helvum. The data of these live non-migratory individuals were not used in the RMA regression. We performed the same spatial analysis as described above based on annual δ²H_P values. We preferred the data set of annual δ²H_P over monthly δ²H_P, since we planned to validate the isoscape origin model that was based on annual δ²H_P. Lastly, we calculated the probability that colony members originated from the two collection sites based on the isoscape origin model.

Results

A Triple-isotope Approach to Geolocate Flying Foxes in Africa

The sampling localities of the 279 fur samples analyzed from African pteropodid bats ranged between 18.08° and −34.36° in latitude, between −40.11° and −16.92° in longitude and between 0–2,355 m above sea level in altitude. Predicted δ²H_P from the sampling locations ranged between −32.8‰ and −3.4‰ according to [31]. Stable isotope ratios in fur keratin of pteropodids ranged from −160.0‰ to −47.5‰ for δ²H_K, from −27.4‰ to −18.3‰ for δ¹³C, and from 3.9‰ to 18.3‰ for δ¹⁵N (Table 1). A MANOVA revealed significant differences in stable isotope ratios (δ²H_k, δ¹³C, and δ¹⁵N) among locations (F_6,57 = 23.4; p<0.001) in the two Epomophorus species. The discriminant function analysis showed that by using a triple-isotope approach, 84% of the E. wahlbergi and E. crypturus samples were assigned correctly to their location of capture (Fig. 1). Some locations formed isolated groups, e.g. those from city of Dar es Salaam and the Usambara and Lubumbashi region, while others presented a less clear separation in the plot, mostly those from the Democratic Republic of Congo (Fig. 1).

[Figure omitted. See PDF.]

Figure 1. Relationships between stable isotope ratios of carbon, nitrogen and hydrogen measured in fur keratin of Epomophorus wahlbergi and E. crypturus.

(A) Relationship between hydrogen and nitrogen stable isotope ratios, (B) relationship between hydrogen and carbon stable isotope ratios, (C) relationship between carbon and nitrogen stable isotope ratios. Each symbol represents one individual. Localities are indicated according to the legend. Abbreviations: DS = Dar es Salaam (Tanzania), KIN (RC) = Kinshasa (Republic of Congo), TAN/USAM (TAN) = Tanga, Usambara (Tanzania), FO/KA = Forzt/Kaloba (Democratic Republic of Congo), KIM/MUK (DRC) = Kimbongo/Mukulakula (Democratic Republic of Congo), KUD (ZIM) = Kudzwe (Zimbabwe), LUB = Lubumbashi (Democratic Republic of Congo).

https://doi.org/10.1371/journal.pone.0045729.g001

The Predictive Power of Hydrogen Stable Isotope Ratios in Fur Keratin

δ²H_K and δ²H_P were significantly correlated in three out of the six non-migratory species. E. wahlbergi and R. aegyptiacus exhibited a positive correlation between δ²H_K and δ²H_P (Pearson’s r = 0.43 and 0.61; p = 0.038 and <0.001 respectively), while the four other species showed either no (E. crypturus: r = −0.23, p = 0.081; L. angolensis: r = −0.4, p = 0.107) or a negative correlation (H. monstrosus: r = −0.47; p = 0.007; and E. franqueti: r = −0.41; p = 0.023). The RMA regression calculated for data of all non-migratory species read: δ²H_P = 35.7+0.66 * δ²H_K. The slope of the regression was significant (P<0.001; Fig. 2).

[Figure omitted. See PDF.]

Figure 2. Relationship between δ²H_P and δ²H_K for non-migratory species (n = 193). δ²H_P was obtained from [44].

The solid line indicates the RMA regression.

https://doi.org/10.1371/journal.pone.0045729.g002

Assessing the Migratory Behavior of Eidolon helvum Based on Hydrogen Stable Isotope Ratios

Based on the developed RMA regression equation and by taking error propagation into account, we transformed values of δ²H_K of E. helvum (n = 88) to create standardized probability density plots. We observed a large variance of δ²H_K in individuals of the same location, suggesting a variable migratory behavior in the study populations. We mapped nine mean probability maps for E. helvum (Fig. 3). In three cases, we obtained probability maps with a low probability of origin for the whole distribution of the species (map G for E. helvum). We were not able to assign a migratory distance to these three individuals. We categorized individuals with a minimum distance of ≤250 km between the capture location and the nearest most probable point of origin (threshold probability: 70%) as non-migratory, and individuals with a minimum distance of >250 km as migratory. In museum specimens of E. helvum, 77.6% (69 out of 85) of the individuals were non-migratory individuals with a median estimated minimum travel distance of 15 km (range: 0–250 km). We considered the remaining proportion (22.4%; 16 out of 85) as migratory, exhibiting a median minimum travel distance of about 860 km (range: 270–3,000 km). Some E. helvum individuals captured at the same locality and in the same campaign (month and/or year labeled at the museum) showed very different origin probabilities. For example, individuals captured in Nouakchott (Mauritania), Sao Tome and Principe, and Eidolon Land (Democratic Republic of Congo), which were classified in different mean probability groups (Fig. 3).

[Figure omitted. See PDF.]

Figure 3. Mean probability maps for 88 individuals of Eidolon helvum (maps A - I).

The IUCN species distribution range shown as an outline with diagonal lines. Based on our isoscape origin model, land areas are shaded on a scale from high (red) to low (blue) probability of origin, with intermediate probabilities indicated by pale yellow shading.

https://doi.org/10.1371/journal.pone.0045729.g003

The stable isotope analysis of contemporary samples revealed that δ²H_k averaged −80.2±5.1‰ (minimum: −93.5‰; maximum: −62.4‰) for individuals from Dar es Salaam and −78.5±6.3 ‰ (minimum: −92.5‰; maximum: −70.4‰) for individuals from Morogoro. The subtle difference was not statistically significant (Student t-test: t₇₄ = 1.19, P = 0.239). Nonetheless, we performed isotope origin model validation for the two sites separately, because δ²H_P differed between the two sites. These analyses showed a high discriminative power for both locations. The likelihood of E. helvum coming from the colony in Morogoro equaled 0.99 and for the Dar es Salaam colony 0.56. For the latter colony, the distance to the nearest spatial cluster with probability values higher than 70% was about 70 km.

Discussion

We studied the migratory behavior of the straw-colored flying fox (Eidolon helvum) an important ecosystem service provider in sub-Saharan Africa, by means of stable isotope analysis and isoscape origin modeling. In our study, we used three steps that we will discuss under separate sub-headings.

A Triple-isotope Approach to Geolocate Flying Foxes in Africa

We used a triple-isotope approach to evaluate whether the Sub-Saharan African isoscape was sufficiently heterogeneous to warrant the delineation of geographically distinct populations based on stable isotope ratios. We were able to distinguish among geographically distinct populations of E. crypturus and E. wahlbergi using stable isotopes of three elements (δ²H_K, δ¹³C_, and δ¹⁵N). The probability of assigning the geographic origin of a fur sample correctly averaged 84% for the seven localities analyzed, which we consider to be satisfactory given the overall low isotopic variation in large areas of the African continent; at least for stable isotopes of hydrogen [31], [34]. The discriminant function analysis showed a clear separation of most populations, even though the majority of samples were from a relatively small area in central Africa (Republic of Congo and Democratic Republic of Congo). Populations with a clear separation from others were those from the town of Dar es Salaam and the Usambara region (both Tanzania), and those collected from Lubumbashi (Democratic Republic of Congo). Samples coming from Lubumbashi were labeled as coming from the zoological garden. Fur keratin of these animals was ¹⁵N-enriched compared with fur keratin of most other pteropodids for which nitrogen isotope ratios are known [51], suggesting that they may have fed on agricultural plants that grew on fertilized soil. Possibly, they were also forming a captive population in the zoological garden and individuals were fed with fruits from some place distant from Lubumbashi or were facing nutritional stress. Nutritional stress is known to increase the variation or to cause systematic shifts in nitrogen isotope ratios [52]–[54]. In general, we assume that the spatial resolution of the multi-isotope approach for African flying foxes could be improved for African flying foxes by including more locations and also more elements in the analysis, as suggested for Afrotropical birds [34].

The Predictive Power of Hydrogen Stable Isotope Ratios in Fur Keratin

Our general aim was to predict δ²H_P of the place of origin for the 88 E. helvum of the two museum collections based on δ²H_k. To this end, we assessed whether δ²H_k is related to annual average δ²H_P of the location where six non-migratory species were captured. Overall, we found positive (E. wahlbergi, R. aegyptiacus), negative (H. monstrosus and E. franqueti) and lack of correlations (E. crypturus, L. angolensis). Likely explanations for the various direction of the relationships are (1) a biased concentration of samples, for example, from only one location or various neighboring locations, resulting in insufficient resolution to differentiate isotopic values between two localities (E. crypturus, E. franqueti and H. monstrosus), (2) low sample sizes per locality, with possible outlier effects (L. angolensis), (3) a high variance of isotopic data obtained from fur of a given capture site, possibly reflecting large variation of δ²H_P for a small area (L. angolensis), (4) using annual δ²H_P instead of monthly δ²H_P as the reference, and/or (5) simply because some species are not truly non-migratory but rather regional migratory species. For example, E. helvum did not show a correlation between δ²H_K and δ²H_P, most likely because sampled E. helvum fur did not grow in the area where animals were captured. In addition, it is important to remember that large areas of the African continent are not surveyed for δ²H of precipitation water [30]. Therefore, predicted δ²H_P may not match the true values of precipitation water when extrapolated based on temperature, rainfall and altitude. Using annual δ²H_P instead of δ²H_P of the month when bats molt may be another important factor that may have obscured a positive correlation between δ²H_K and δ²H_P. Seasonal variation in δ²H_P is large on the African continent [31], and therefore it is possible that we referred to a false value for δ²H_P. Nonetheless, when analyzed together, we found a statistically significant correlation between δ²H_k and δ²H_P across all non-migratory species, even though intra-specific correlations of δ²H_k and δ²H_P did not show a consistent pattern, and even though we were not able to take the strong seasonal variation of δ²H_P in Sub-Saharan Africa into account.

Assessing the Migratory Behavior of Eidolon helvum Based on Hydrogen Stable Isotope Ratios

For each E. helvum, we created predictive maps of origin using the established RMA and gamma distributions. This allowed the inaccuracy of measuring hydrogen isotopic ratios from fur samples and in predicting the hydrogen isotope ratios of local precipitation water to be taken into account. Our results showed that most E. helvum (78%; 69 out of 85) were captured from non-migratory populations. We estimated a median minimum travel distance of 15 km for the non-migratory animals, which is similar to what has been observed as daily commuting distances in R. aegyptiacus [7]. The ground-truthing of our isoscape origin model, based on stable isotope analysis in extant individuals of E. helvum, showed a high accuracy in predicting the origin of the inland colony in the city of Morogoro and a medium probability of origin for the individuals of the coastal site located in the city of Dar es Salaam. Possibly, the predicted δ²H_P of the coastal site is more error-prone than values from inland sites because of the high humidity and various influxes of water at coastal sites. Also, according to our estimate non-migratory E. helvum from the Dar es Salaam colony could have travelled for 70 km and reached δ²H_P that were similar to their δ²H_k. The successful validation of our isoscape origin model for E. helvum is noteworthy because it was based on two large assumptions. Firstly, we referred to annual δ²H_P instead of monthly δ²H_K values and this should have increased the variance or added a systematic bias to our extrapolation. Secondly, we lacked information about the exact molting period of E. helvum. We assumed that molting would occur before a possible migration period. The successful validation of our isoscape origin model therefore supports the notion that these two assumptions were not largely violated in our stable isotope study of E. helvum.

We categorized 16 out of 85 E. helvum (22%) as migratory. The median minimum travel distance of migratory animals was 860 km, and four individuals seem to have traveled over a distance of more than 2,000 km prior to capture. This observation corroborates earlier findings that E. helvum is capable of performing long-distance migrations when following the seasonal fluctuation in fruit abundance [5]. Based on our literature review and the capture information from some museum specimens included in this study (only available for 7 localities), we assume that E. helvum has a core distribution in equatorial Africa, with migrations in northern direction, e.g. Mauritania, and Niger, from May to September and towards the south, e.g. Tanzania, Zimbabwe and Zambia during the months of October and December. Individuals captured at the same time and at the same locality in Central Africa exhibited contrasting isotopic signatures in their fur, and in consequence are likely to have originated from different locations. For example, this pattern was apparent in E. helvum individuals captured in Nouakchott (Mauritania; n = 4), where two individuals presented an isotopic origin from the same locality, yet two individuals originated from localities at around 250 km and 850 km away from the collection site. Individuals captured at Eidolon land, near the National Park of Kundelungu (Democratic Republic of Congo; n = 10) also showed this pattern, with eight non-migratory individuals, and two having migrated over a distance of about 700 km.

For three individuals, we obtained a mean probability map with a low probability of origin for the whole continent (E. helvum group G). This was caused by low δ²H_k values that were not typical for the average annual δ²H_P of the African continent, i.e. once transformed with the RMA and simulated with the gamma function, these values remained outside the data range of δ²H_P of Africa. One possible explanation is that our isoscape origin model was based on annual mean δ²H_P of Africa [45], and these individuals came from regions where stable hydrogen isotope ratios vary largely during a short period of time, as it is the case for example in West Africa in August (www.waterisotopes.org). Consequently, these bats may have carried a stable hydrogen isotope signature that was representative for a small area or a brief period of time for the specific capture location. Alternatively, these specimens may have been maintained in captivity and fed with non-local fruits between collection and shipment to the museum collections. This raises an important point about using museum specimens in stable isotope studies: circumstances of capture and the exact Julian date of capture are usually not documented in museum specimens. Also, it is not usually documented how specimens were treated after sacrificing and shipping them to the museums, e.g. by using different preservatives at the site of collection, during shipment and for permanent storage in the museums.

Museum Collections and Stable Isotope Studies

Specimens from museum collections are increasingly used to study the isotopic ecology of vertebrates. In our study, the use of specimens from museum collection was constrained, because of the lack of information about the Julian date of capture or the change in the naming of locations over the past century. Also, it was unknown for all studied museum specimens whether bats were captured in their natural habitat or, for example, purchased from local markets. Our isotopic origin model would have suffered if bats collected for museums originated from other locations than indicated by the collector. However, since it is unlikely that wrong location assignments caused a systematic bias, we are confident that the general picture of our study remains unaffected by this possibility. In addition, ethanol preserving could have increased the variance in stable isotope ratios in museum specimens, if ethanol interacts with the keratinous protein. Yet, we found a positive correlation between δ²H_k and δ²H_P in non-migratory species, suggesting that effects of preservatives were probably absent or negligible. In summary, museum collections allowed us to obtain a high number of samples at a low cost and to provide new data about the movement ecology of an ecological and economical key species on the African continent.

Conclusions

Our study is one of the first to show that stable isotope ratios can be used to predict the provenance of animals on the African continent. One recent study in a long-distant bird migrant carrying geolocators confirmed that feather δ²H_K of northern wheatear, Oenanthe oenanthe, recorded the variation of the predicted δ²H_P of the overwintering sites [55]. Other studies have also supported a large potential of using stable isotope ratios for investigating the movement ecology of Afro-European birds [34], [56]. Yet except for these and some other papers, stable isotope studies on migratory African birds or bats have been rare or absent respectively. The major reason for this is that the isoscape of the African continent has been considered to be less informative than those of, for example, Europe and North America [31]. Our study shows that a triple-isotope approach (δ²H_k, δ¹³C and δ¹⁵N) is useful to discriminate between fur samples of flying foxes that originated from various locations in Central Africa. However, since this approach requires a complete baseline map for the stable isotope ratios in keratin of non-migratory species, it may not be feasible on a larger geographic scale. Therefore, we have validated that δ²H_K is correlated positively with δ²H_P across individuals of non-migratory flying foxes. Based on the developed RMA regression and the consideration of error propagation, we used an isoscape origin model to estimate the movement ecology of E. helvum. Using historical specimens of this species, we confirmed that a large portion of the population was non-migratory and that some individuals performed long-distance migrations of more than 2,000 km; a finding that corroborates earlier studies using satellite tracking [18]. Overall, we observed travel distances of only about 870 km in migratory individuals. This could be attributed to the observation that E. helvum may perform several small-scale migrations of a few 100 km from one stepping stone to the other, following the seasonal abundance of local food resources on the Sub-Saharan African continent.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Location of collection sites for each of the seven flying fox species. The African continent is coloured in relation the annual average stable hydrogen isotope ratio of precipitation water according to [31]. The distribution range of flying foxes is marked as a shaded area.

https://doi.org/10.1371/journal.pone.0045729.s001

(TIF)

Table S1.

Reference numbers of the sampled specimens at the National Museum for Natural History in Paris (MNHN) and the Natural History Museum in Berlin (MFN) and the corresponding stable isotope ratios for carbon (δ¹³C), nitrogen (δ¹⁵N) and hydrogen (δ²H_K). Abbreviations for countries are: ZIM = Zimbabwe. DRC = Democratic Republic of Congo. GHA = Ghana. CAM = Cameroon. NIG = Nigeria. ANG = Angola. MAU = Mauritania. TAN = Tanzania. SAF = South Africa. STP = Sao Tome and Principe. IVO = Ivory coast. CAR = Central African Republic. RC = Republic of Congo (na = not available, Alt. = altitude, Lat. = latitude, Long. = longitude).

https://doi.org/10.1371/journal.pone.0045729.s002

(DOCX)

Acknowledgments

We thank Jean-Marc Pons (MNHN-Paris), Frieder Mayer (MFN-Berlin) and Nora Lange (MFN-Berlin) for the generous authorizations to use specimens from the museum collections, Jean-Francois Julien and Julie Marmet for their help to obtain samples in Paris, Anja Luckner and Doris Fichte for their guidance and help in the preparation of samples and for the isotopic analyses, and Ana Popa-Lisseanu, Karin Schneeberger and Daniel Lewanzik for comments on earlier versions of the manuscript. For permission and assistance with contemporary sampling in Tanzania, we thank the staff at the Tanzanian Wildlife Research Institute (TAWIRI), the Ministry of Livestock Development and Fisheries and Sokoine University of Agriculture. We also thank two anonymous reviewers for constructive comments on the manuscript.

Author Contributions

Conceived and designed the experiments: CCV SKS. Performed the experiments: GO AJP AKS. Analyzed the data: CCV SKS GO AKS AJP. Contributed reagents/materials/analysis tools: CCV AJP. Wrote the paper: CCV GO SKS.

Citation: Ossa G, Kramer-Schadt S, Peel AJ, Scharf AK, Voigt CC (2012) The Movement Ecology of the Straw-Colored Fruit Bat, Eidolon helvum, in Sub-Saharan Africa Assessed by Stable Isotope Ratios. PLoS ONE 7(9): e45729. https://doi.org/10.1371/journal.pone.0045729

References

1. Kunz TH, Braun de Torrez E, Bauer DM, Lobova TA, Fleming TH (2011) Ecosystem services provided by bats In: Ostfeld RA, Schlesinger WH, editors. The Year in Ecology and Conservation 2011. New York, USA, Annals of the New York Academy of Sciences, Wiley-Blackwell, 38 p.

2. Ghanem SJ, Voigt CC (2012) Increasing awareness of ecosystem services provided by bats. Advances in the Study of Behavior 44: 279–302.

3. Asenyu ES (1974) Plant and bat interactions in West Africa. Annals Miss Bot Gard 61: 702–727.

4. Fujita MS, Tuttle MD (1991) Flying foxes (Chiroptera: Pteropodidae): Threatened animals of key ecological and economic importance. Cons Biol 5: 455–463.

5. Richter HV, Cumming GS (2008) First application of satellite telemetry to track African straw-coloured fruit bat migration. J Zool 275: 172–176.

6. Smith CS, Epstein JH, Breed AC, Plowright RK, Olival KJ, et al. (2011) Satellite telemetry and long-range bat movements PLoS ONE. 6: e14969.

7. Tsoar A, Nathan R, Bartan Y, Vyssotski A, Dell’Omo G, et al. (2011) Large-scale navigational map in a mammal. Proc Nat Acad Sci U S A 108: E718–E724.

8. Leroy EM, Kumulungui B, Pourrut X, Rouquet P, Hassanin A, et al. (2005) Fruit bats as reservoirs of Ebola virus. Nature 438: 535–576.

9. Calisher CH, Childs JE, Field HE, Holes KV, Schountz T (2006) Bats: important reservoir hosts of emerging viruses. Clin Microbiol Rev 19: 531–545.

10. Wibbelt G, Moore MS, Schountz T, Voigt CC (2010) Emerging diseases in Chiroptera: Why bats? Biol Lett 6: 438–440.

11. Drexler JF, Corman VM., Müller MA, Maganga GD, Vallo P, et al. (2012) Bats host major mammalian paramyxoviruses. Nat Comm doi 10.1038/ncomms1796.

12. DeFrees SL, Wilson DE (1988) Eidolon helvum. Mammal Species. 312: 1–5.

13. Harris BJ, Baker HG (1959) Pollination of flowers by bats in Ghana. Nigerian Field 24: 151–159.

14. Daïnou K, Bizoux JP, Doucet JL, Mahy G, Hardy OJ, et al. (2010) Forest refugia revisited: nSSRs and cpDNA sequences support historical isolation in a wide-spread African tree with high colonization capacity, Milicia excelsa (Moraceae). Mol Ecol 19: 4462–4477.

15. World Conservation Monitoring Center (1998) Milicia excelsa In: IUCN 2011 IUCN Red List of Threatened Species Version 20112. Available: http://wwwiucnredlistorg Accessed 2012 February 10.

16. Baker HG (1963) Evolutionary mechanisms in pollination ecology. Science 139: 877–883.

17. Racey P (2004) 8,000,000 fruit bats: Africa’s best-kept wildlife secret. Bats 22: 1–5.

18. Richter HV, Cumming GS (2005) Food availability and annual migration of the straw-colored fruit bat (Eidolon helvum). J Zool 268: 35–44.

19. Thomas D (1983) The annual migration of three species of West African fruit bats. Can J Zool 61: 2266–2272.

20. Cosson JF, Tranier M, Colas F (1996) On the occurrence and possible migratory behaviour of the fruit bat Eidolon helvum in Mauritania, African J Afr Zool. 110: 369–371.

21. Henry M, Barrière P, Gautier-Hion A, Colyn M (2004) Species composition, abundance and vertical stratification of a bat community (Megachiroptera: Pteropodidae) in a West African rain forest. J Trop Ecol 20: 21–29.

22. Niamien CJM, Yaokokoré-Béibro HK, Koné I, N’Goran KE (2010) Données préliminaires sur l’écologie des chauves-souris frugivores de la commune du Plateau (Abidjan, Côte d’Ivoire). Sciences Nature 7: 21–30.

23. Kelly JF, Atudorei V, Sharp ZD, Finch DM (2002) Insights into Wilson’s Warbler migration from analyses of hydrogen stable-isotope ratios. Oecologia 130: 216–221.

24. Hobson KA (2011) Isotopic ornithology: A perspective. J Ornith 152: 49–66.

25. Cryan PM, Bogan MA, Rye RO, Landis GP, Kester CL (2004) Stable hydrogen isotope analysis of bat hair as evidence for seasonal molt and long-distance migration. J Mammal 85: 995–1001.

26. Popa-Lisseanu AG, Sörgel K, Luckner A, Wassenaar L, Ibañez C, et al. (2012) A triple approach to predict the breeding origins of European bats. PLos ONE 7: e30388.

27. Voigt CC, Popa-Lisseanu AG, Niermann I, Kramer-Schadt S (2012) The catchment area of wind farms for European bats: A plea for international regulations. Biol Conserv 153: 80–86.

28. Hobson KA, Norris DR (2008) Animal migration: a context for using new techniques and approaches In: Hobson KA, LI Wassenaar, editors. Tracking animal migration using stable isotopes. London, Academic Press.

29. Hobson KA, Barnett-Johnson R, Cerling T (2009) Using isoscapes to track animal migration In: West JB, Bowen GJ, Dawson TE, Tu KP, editors. Isoscapes: Understanding movement, pattern, and process on earth through isotope mapping. New York, Springer. 299–318.

30. Bowen GJ, Revenaugh J (2003) Interpolating the isotopic composition of modern meteoric precipitation. Water Res 39: 1299.

31. Bowen GJ, Wassenaar LI, Hobson KA (2005) Global application of stable hydrogen and oxygen isotopes to wildlife forensics. Oecologia 143: 337–348.

32. Bowen GJ (2009) Statistical and geostatistical mapping of precipitation water isotope ratios. In: West JB, Bowen GJ, Dawson TE, Tu KP, editors. Isoscapes: Understanding movement, pattern, and process on earth through isotope mapping. New York, Springer.

33. Craig H (1961) Isotopic variations in meteoric waters. Science 133: 1702–1703.

34. Hobson KA, Van Wilgenburg SK, Wassenaar LI, Powell RK, Still CJ, et al. (2012) A multi-isotope (δ13C, δ15N, δ2H) feather isoscape to assign Afrotopical migrant birds to origins. Ecosphere 3: 44.

35. Hartmann G, Danin A (2010) Isotopic values of plants in relation to water availability in the Eastern Mediterranean region. Oecologia 162: 837–852.

36. African Chiroptera Report (2011) African Chiroptera proyect. Available: http://www.Africanbats.org accessed on 10 February 2012.

37. Kwiecinski GG, Griffiths TA (1999) Rousettus aegyptiacus Mammal Species. 611: 1–9.

38. Jacobsen NHG, Du Plessis E (1976) Observations on the ecology and biology of the Cape fruit bat Roussetus aegyptiacus leachi in the Eastern Transvaal. South Afr J Sci 72: 270–273.

39. Acharya L (1992) Epomophorus wahlbergi. Mammalian Species. 394: 1–4.

40. Langevin P, Barclay RMR (1990) Hypsignathus monstrosus. Mammal Species 357: 1–4.

41. Mickleburgh S, Hutson AM, Bergmans W, Fahr J, Juste J (2008) Epomops franqueti In: UICN 2011 UICN Red List of Threatened Species Version 20112 Available: www.iucnredlist.org.Accessed 2012 March 5.

42. Mickleburgh S, Hutson AM, Bergmans B, Fahr J, Racey PA (2008) Eidolon helvum In: IUCN 2011 IUCN Red List of Threatened Species Version 20112. Available: www.iucnredlist.org. Accessed 2012 February 10.

43. CMS (2012) Convention on Migratory Species. Available: http://www.cms.int. Accessed 2012 March 10.

44. Google Inc. (2012) Google Earth (version 6.1.0.5001) [Software]. Available: www.google.com/earth/index.html. Accessed on 2012 February 10.

45. Wassenaar LI, Hobson KA (2003) Comparative equilibration and online technique for determination of non-exchangeable hydrogen of keratins for use in animal migration studies. Isot Env Health Stud 39: 211–217.

46. Bowen GJ (2012) Gridded maps of the isotopic composition of meteoric waters. Available: http://www.waterisotopes.org accessed 10 May 2012.

47. Erzberger A, Popa-Lisseanu A, Lehmann GUC, Voigt CC (2011) Potential and limits in detecting altitudinal movements of bats using stable hydrogen isotope ratios of fur keratin. Acta Chiropterol 13: 431–438.

48. R development Core Team (2012) R: A language and environment for statistical computing R foundation for statistical computing Vienna, Austria ISBN 3-900051-07-0. Available: www.R-project.org. Accessed 2012 May 3.

49. Taylor PJ, Monadjem A (2008) Maxillary shape as a diagnostic tool for identifying fruit bats, Epomophorus crypturus and E wahlbergi from museum specimens and in the field. South Afr J Wildl Res 38: 22–27.

50. Legendre P (2008) Lmodel2: Model II Regression R package version 16–3.

51. Voigt CC, Akbar Z, Kunz TH, Kingston T (2011) The origin of assimilated proteins in Old and New-World phytophagous bats. Biotropica 43: 108–113.

52. Voigt CC, Matt F (2004) Nitrogen stress causes unpredictable enrichments of 15N in two nectar-feeding bat species. J Exp Biol 207: 1741–1748.

53. Fuller BT, Fuller JL, Sage NE, Harris DA, O'Connell TC, et al. (2005) Nitrogen balance and 15N: why you're not what you eat during nutritional stress. Rapid Communications in Mass Spectrometry 19: 2497–2506.

54. Fuller BT, Fuller JL, Sage NE, Harris DA, O'Connell TC, et al. (2004) Nitrogen balance and 15N: why you're not what you eat during pregnancy. Rapid Communications in Mass Spectrometry 18: 2889–2896.

55. Bairlein F, Norris DR, Nagel R, Bulte M, Voigt CC, et al. (2012) Cross-hemisphere migration of a 25-gram songbird. Biol Lett doi: 10.1098/rsbl.2011.1223 1744-957X.

56. Chamberlain CP, Bensch S, Feng X, Åkesson S, Anderson T (2000) Stable isotopes examinated aross a migratory divide in Scandinavian willow warblers (Phylloscopus trochilus and Phylloscopus trochilus acredula) reflect their African winter quarters. Proc Roy Soc Lond B 267: 43–48.