The Foraging Ecology of the Mountain Long-Eared Bat Plecotus macrobullaris Revealed with DNA Mini-Barcodes

Abstract

Molecular analysis of diet overcomes the considerable limitations of traditional techniques for identifying prey remains in bat faeces. We collected faeces from individual Mountain Long-eared Bats Plecotus macrobullaris trapped using mist nets during the summers of 2009 and 2010 in the Pyrenees. We analysed their diet using DNA mini-barcodes to identify prey species. In addition, we inferred some basic features of the bat's foraging ecology that had not yet been addressed. P. macrobullaris fed almost exclusively on moths (97.8%). As prey we detected one dipteran genus (Tipulidae) and 29 moth taxa: 28 were identified at species level (23 Noctuidae, 1 Crambidae, 1 Geometridae, 1 Pyralidae, 1 Sphingidae, 1 Tortricidae), and one at genus level (Rhyacia sp., Noctuidae). Known ecological information about the prey species allowed us to determine that bats had foraged at elevations between 1,500 and 2,500 m amsl (above mean sea level), mostly in subalpine meadows, followed by other open habitats such as orophilous grasslands and alpine meadows. No forest prey species were identified in the diet. As 96.4% of identified prey species were tympanate moths and no evidence of gleaning behaviour was revealed, we suggest P. macrobullaris probably forages by aerial hawking using faint echolocation pulses to avoid detection by hearing moths. As we could identify 87.8% of the analysed sequences (64.1% of the MOTUs, Molecular Operational Taxonomic Units) at species level, we conclude that DNA mini-barcodes are a very useful tool to analyse the diet of moth-specialist bats.

About the Authors:

Antton Alberdi

Affiliation: Department of Zoology and Animal Cell Biology, Faculty of Science and Technology, University of the Basque Country UPV/EHU, Sarriena z.g., Leioa, The Basque Country

Inazio Garin

Affiliation: Department of Zoology and Animal Cell Biology, Faculty of Science and Technology, University of the Basque Country UPV/EHU, Sarriena z.g., Leioa, The Basque Country

Ostaizka Aizpurua

Affiliation: Department of Zoology and Animal Cell Biology, Faculty of Science and Technology, University of the Basque Country UPV/EHU, Sarriena z.g., Leioa, The Basque Country

Joxerra Aihartza

* E-mail: [email protected]

Affiliation: Department of Zoology and Animal Cell Biology, Faculty of Science and Technology, University of the Basque Country UPV/EHU, Sarriena z.g., Leioa, The Basque Country

Introduction

The study of the trophic resources used by a species and the habitats where they are consumed are key aspects of addressing foraging ecology, which will provide a basic understanding of the relationships among consumers, resources, and environment [1]. Elucidating diet and habitat preferences is, therefore, paramount for any management or conservation purposes. Nevertheless, such studies are often beset by problems such as difficulty in obtaining information on elusive animals, capture/handling restrictions imposed by conservation status, or other methodological constraints. Thorough investigations of diet may provide sufficient information on most animals' foraging requirements, but such detailed estimates have hitherto been hard to achieve due to limitations in methods used.

Bat diet studies began with analysis of prey remains collected under feeding perches (e.g. [2], [3]) and inspection of stomach contents from sacrificed bats [4], [5]. The former technique digs out mostly the biggest prey [6], i.e. those dismembered prior to ingestion, and is limited to bat species that use perches and to individuals with known perching sites. The latter method entails sacrifice of animals that usually belong to protected species, and it is no longer used.

In recent decades, most dietary studies on bats have been carried out through morphological identification of prey fragments enduring in faeces. This approach depicts the diet realistically [7], [8], and it has allowed researchers to analyse allocation of food resources among separate bat species [9]–[11] or to study prey selection through comparison of diet and prey abundance in foraging areas (e.g. [12], [13]). In fact, a comprehensive feeding ecology of nearly all bat species has been described using this technique (e.g. [14]–[16]). Underrepresentation of soft-bodied prey is usually assumed, however, because their identifiable parts are less likely to persist [17]. Moreover, the key morphological features used to identify lower taxa are fatally damaged through digestion. Consequently, prey remains are seldom identified below the ordinal or family level [7], [18]. To increase resolving power, some authors have combined morphological identification with stable isotopic analysis of faeces [19], though they were only able to suggest the families of ingested items.

The development of molecular techniques has now taken the analysis of diet a qualitative step forward [20]. Along with the standardisation of a single molecular marker and the development of a large reference database (BOLD Systems, www.boldsystems.org; [21]), diet studies increasingly employ DNA barcodes [22]–[24]. The animal DNA barcode is a small fragment of the mitochondrial genome (COI gene) that is being sequenced in many species as an identity label [25], [26], which allows species-level description of an animal's diet.

Beyond a mere listing of prey, species-level identification affords the opportunity to test hunting-strategy hypotheses that could not previously be addressed with analysis techniques lacking the required resolution. For instance, the predator-prey relationship between bats and eared moths is one of the best exponents of coevolutionary arms race [27]. Although sensorial adaptations and related evasive and attacking behaviours have been extensively studied [28]–[35], strikingly, the actual contribution of eared moths to any bat species' diet has seldom been revealed to the family level [36]–[38]. Traditional techniques are generally unable to identify the species of consumed moths, so insight on precise predator-prey relationship could only be occasionally ascertained, mostly when culled parts were available [39]–[41]. Conversely, molecular identification of a tympanate moth as one of the main prey of Barbastella barbastellus [24] enabled an investigation showing that some aerial hawking bats hunt tympanate moths by casting faint echolocation calls that overcome the moths' hearing ability [42]. As many biological and environmental factors contribute to the diet of bats, molecular analysis of prey can also be a handy tool to spot a wide spectrum of habitat parameters. For example, Clare et al. [43] inferred ecosystem-level features by assessing the quality of water environments where bats fed through a quality ranking of the source environments of prey identified by molecular means.

The Mountain Long-eared Bat Plecotus macrobullaris, Kuzjakin 1965, was accepted as a species in 2002 [44]–[46] (Fig. S1). Most of the published research on the species has referred to its morphology or phylogenetics [47]–[51]. Its ecology is still poorly described, and its foraging preferences are controversial. Breeding colonies have been found in man-made buildings at 400–1,300 m in the Alps [51], [52] and as low as at sea level in Croatia [48]. P. macrobullaris has been captured by mist net in pastures above 1,800 m and as high as 2,800 m in the Pyrenees [53], [54] and in montane open areas with sparse shrubbery in the Middle East [50]. Accordingly, a recent study [51] documented the resemblance of the species' trophic niche to that of the open-land forager P. austriacus [55]. On the other hand, the echolocation and morphological characteristics of P. macrobullaris have prompted some authors to believe it is suited to forage in narrow rather than open spaces [51], [56]. This latter hypothesis seems to be supported by mist-net captures in Croatian forests and roost placements in Switzerland [48], [57]. Nevertheless, in a recently published radio-tracking study carried out in a nursery colony in Northern Italy bats avoided woods when foraging, selecting open rural areas and ecotones [58].

Recent studies on the diet of P. macrobullaris have revealed the importance of moths [16], [59], though the exact species composing its diet remain unidentified. Fortunately, Lepidoptera is one of the insect taxa with an extensive database of DNA barcodes available online (more than 70,000 species). Therefore, in this study we aimed to unveil the foraging behaviour of the Mountain Long-eared Bat P. macrobullaris. Particularly, we broadly tested the bat's alpine nature by investigating the elevational distribution and habitat of prey. Meanwhile, through detailed inspection of the prey species, we tested whether the Mountain Long-eared Bat feeds on hearing and non-flying prey. We used molecular tools to identify consumed prey at species level, and relied on available information about the habitats of moths and associated host plants. We have also assessed the usefulness and reliability of molecular tools in diet analysis.

Materials and Methods

Ethics Statement

Capture and handling protocols met the guidelines for treatment of animals in research and teaching [60], were approved by the Regional Council, and met Spanish legal requirements. The Government of Aragón (Spain) gave the necessary permits to carry out the fieldwork (LC/mp 24/2009/2958 and LC/ehv 25/2010/3234). Bats were released after being measured and faeces collected. To minimise stress, retention time never exceeded 30 minutes.

Sample collection

Sampling was carried out during the summers of 2009 and 2010 in 25 different locations of the Pyrenees mountain range. We captured bats by mist nets placed at commuting paths, water points, and putative foraging areas 1,000–2,100 m amsl. Habitat types sampled were coniferous and holm oak forests (1,000–1,400 m, n = 6 sampling sites), subalpine meadows (1,500–1,900 m, n = 7), alpine meadows (1,700–2,100 m, n = 7), and rocky areas with poor vegetation cover (1,900–2,600 m, n = 5). In the Pyrenees, as in the Alps and many other European mountain ranges, the tree line has been lowered for pastoral activities by clearing subalpine shrubs and forests [61]–[64]. Consequently, the lower boundary of some alpine meadows in the Pyrenees has been lowered as much as 500 m.

Each trapped bat was held individually in a clean cloth bag for 10–25 min to collect faeces. Droppings were stored dry and frozen within 6 h of collection.

DNA isolation and amplification

Each individual bat was considered as a sampling unit [65]. For each bat specimen, 20–40 mg of faeces (2–6 pieces of guano) were used for DNA extraction with the DNA Stool Mini Kit (Qiagen), following the manufacturer's instructions with some modifications [24].

Although a 650-bp sequence has been standardised as the molecular barcode [66], DNA degradation in droppings commonly forces the use of shorter sequences, even less than 200 bp [67]. Those short sequences perform at 90% species resolution in taxa with high species-level specificity [68], [69]. Moreover, their usefulness for prey identification in droppings has successfully been tested [24]. Consequently, a 157-bp fragment of COI gene was amplified using ZBJ-ArtF1c and ZBJ-ArtR2c primers developed by Zeale et al. [24]. The PCRs were performed using the Biotaq PCR kit (Bioline, www.bioline.com) in 10-µl-volume reactions. Each reaction contained 1.05 µl deionised water, 1 µl NH4×10 buffer, 0.4 µl 50 mM MgCl₂, 1 µl 2 mM dNTPs, 0.75 µl 10 µM forward primer, 0.75 µl 10 µM reverse primer, 0.05 µl BIOTAQ DNA polymerase, and 5 µl of DNA sample extracted from the faeces. We used the PCR protocol for amplification developed by Zeale et al. [24]. PCR products were visualised on a 1.5% agarose gel to determine the amplicon∶primer-dimers ratio and choose the correct reaction products for cloning. PCR reactions were optimised to avoid the purification step of PCR-products, preventing DNA loss, and raising the ligation efficiency.

DNA cloning and sequencing

The pGEM-T Easy Vector System and high-efficiency competent cells (≥1×10⁸ cfu/µg DNA) (Promega www.promega.com) were used for cloning the PCR products. A 3∶1 insert-to-vector molar ratio was used for ligation and reactions were carried out with overnight incubation at 4°C. Transformation was performed following the manufacturer's instructions. Twenty colony-forming units (cfu) were selected from each plate, and cloned inserts were liberated from cells and amplified using the universal M13 primers. Each PCR contained 20.8 µl deionised water, 2.5 µl NH4×10 buffer, 0.75 µl 50 mM MgCl₂, 0.5 µl 10 mM dNTPs, 0.1 µl 0.2 mM M13 forward primer (5′-CCCAGTCACGACGTTGTAAAACG-3′), 0.1 µl 0.2 mM M13 reverse primer (5′-AGCGGATAACAATTTCACACAGG-3′), and 0.25 µl Taq DNA polymerase to complete a reaction volume of 25 µl. PCR thermal cycling conditions were as follows: 5 min at 94°C (longer than usual, to lyse bacteria), followed by 35 cycles of 45 s at 94°C, 45 s at 50°C, and 45 s at 72°C. The cycles ended with one final extension of 8 min at 72°C. Resulting amplicons were unidirectionally sequenced using M13 forward primer.

Data analysis

The obtained sequences were aligned and edited with BioEdit [70], and deposited in the Dryad database (http://dx.doi.org/10.5061/dryad.611310kt). We used the software jMOTU [71] to collapse sequences into Molecular Operational Taxonomic Units (MOTUs) setting the threshold-value to 1.3% (2 nucleotides). This is less conservative than the threshold-values used by Bohmann et al. [23], Clare et al. [43] and Razgour et al. [72] –2.5%, 2.0% and 2.0% respectively–, but given the small sampling range we anticipated less diversity among available insects. The most common sequences from each MOTU were used for species identification. BLAST was used for sequence similarity analysis, comparing query sequences to reference sequences stored in GenBank derived from BOLD [21]. Taxonomic level thresholds were chosen based on our own data and bibliography [22], [43]. Species-level identification was assigned when similarity between query and reference sequences was greater than 98.7%, i.e. equal or less than 2 mutations. Similarities from 97.4% to 98.7% were identified at genus level. Similarities under 97.4% were not used to assign any sequence to a higher taxonomic level because they might be error prone [43]. As some MOTUs matched more than one species above 98.7% level, we gave different confidence levels to the identifications, based on the following criteria:

1. Single match: the sequence matches a single species with the selected similarity threshold (2 mutations),

2. Several global/single local matches: the sequence matches more than one species with the selected similarity threshold, but just one of them is present in Western Europe.

3. Several global/several local matches: the sequence matches more than one species present in Western Europe with the selected similarity threshold, but just one of them matches at 100%.

Rarefaction curves with 1,000 iterations were carried out for each sample with EstimateS 8.2.0 [73]. Dietary diversity was determined at the MOTU level using the Shannon's diversity index, where p_i was defined as the number of samples in which MOTU i was found divided by sum of occurrences of all MOTUs [72].

Prey item analysis

We used the relationship between body mass and wing-length in noctuid moths published by Wood et al. [74], to assess average body mass of consumed prey. Moths consumed by P. macrobullaris were grouped according to their elevational preference as follows: (1) widespread species occurring in a wide elevational range from sea level up to alpine meadows (0–2,800 m), (2) orophilous-alpine species appearing in both montane and alpine areas (600–2,500 m), or (3) strictly alpine species (>1,500 m).

Identified moth species were also classified according to their known foraging habitat preferences ([75]–[81] and www.nmnhs.com/butterfly_areas_bg/, accessed in April 2011). We distinguished 6 habitat types: (1) alpine meadows, grasslands extending up from the tree line, i.e. the imaginary line delimiting the highest elevations where trees at least 3 m tall are found in distinct patches [82]; (2) subalpine meadows, grasslands extending from the timberline to the tree line [82]; (3) orophilous grasslands, meadows below the timberline; (4) rockies, barren lands with vegetation cover less than 15% [83]; (5) shrublands, areas with shrub canopy (heath, juniper, etc.) greater than 20% of total vegetation; and (6) forests, either coniferous or deciduous woods.

Results

In 12 out of the 25 sampled locations we captured 40 Mountain Long-eared Bats. DNA with sufficient quality for analysis was extracted from 31 of 40 individual faecal samples (77%), and at least 20 colonies were grown in 29 of those 31 (93%); thus, 102 faecal pellets produced by 29 individuals from 12 different locations were analysed. Rarefaction curves for individuals' samples showed that 22 out of 29 clearly reached the asymptote, 4 almost did, and 3 clearly did not (Fig. 1).

[Figure omitted. See PDF.]

Figure 1. Diversity of prey sequences detected in samples in relation to number of clones sequenced.

The rarefaction curves for individuals' samples showed that 22 out of 29 clearly reached the asymptote (green), four almost did (yellow), and three clearly did not (red), which means that the amount of clones sequenced per sample was high enough to detect the great majority of prey sequences in it.

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

We analysed 20 colonies/sequences from each individual faecal sample, examining a total of 580 DNA sequences, among which we found 90 different sequences that were collapsed into 54 MOTUs. We identified 34 of the MOTUs to the likely species and 11 to genus level. All except one came from insects (Table 1): 44 in the Order Lepidoptera and one in the Order Diptera (Genus Tipula). Eight MOTUs could not be identified to species or genus level and the remaining one may correspond to bacteria, fungi, or another microbial element.

[Figure omitted. See PDF.]

Table 1. Prey taxa detected in the diet of Plecotus macrobullaris.

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

We successfully identified 28 moth species (539 sequences, 31 MOTUs) and one moth genus (Rhyacia) (Table 1). In some cases more than one different MOTU was assigned to the same species: three MOTUs were identified as Agrotis simplonia, and two as Apamea zeta; similarly, 11 MOTUs were linked to five genera, and some of them were subsequently identified at species level. Among the moths identified at species level, 27 species (23 Noctuidae, 1 Geometridae, 1 Crambidae, 1 Sphingidae, 1 Pyralidae) and one genus (Noctuidae) were tympanate lepidopterans, while one species (Tortricidae) was not; 98.1% of species-level-identified MOTUs belonged to eared moths, while 1.9% belonged to non-tympanate lepidopterans.

The average wingspan of prey moth species in the diet of the Mountain Long-eared Bat was 37.60 mm (SD ± 7.87, range 21–55). We calculated that P. macrobullaris feeds on moth species with an average body mass of 177 mg (SD ± 82, range 35–281).

The average number of prey species per bat was 2.75 (SD ± 1.6, range 1–8), more than 50% of the individuals had only one or two species in their faeces, and 28 of the 29 individuals showed five species or fewer, while the faeces of one exceptional bat contained 8 species. The average number of prey species per sampling site was 5.58 (SD ± 3.5, range 1–13). Shannon's index of dietary diversity was 3.76.

The identified prey species occupy very diverse elevational ranges (Fig. 2). Strictly alpine species were found in 15 samples and orophilous-alpine species in 13 samples. Species with wide elevational ranges appeared together with strictly alpine species in five individuals, with orophilous-alpine species in four individuals, and alone in six. Most of the sequences belonged to species with broad elevation range (42.8%), followed by orophilous-alpine species (32.3%) and strictly alpine species (25.0%). The leading habitat types used by moths consumed by P. macrobullaris are subalpine meadows, followed by orophilous meadows and alpine meadows (Fig. 3). Thus, almost all the habitat types used by consumed moth species are open habitats. In addition, 76% of the detected moths use host plants from grasslands.

[Figure omitted. See PDF.]

Figure 2. Elevational distribution of the moth prey species identified in faeces.

a) Number of identified prey species present in each elevational range; b) elevational ranges of each prey identified at species level in South Europe. Species with wide elevational ranges are drawn in light green, orophilous-alpine species in yellow, and strictly alpine species in magenta. Bibliographic sources are given in the main text.

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

[Figure omitted. See PDF.]

Figure 3. Foraging habitats of the prey species identified in faeces.

a) Number of identified prey species foraging in each habitat type; b) habitat types used by each identified prey species in South Europe. Used habitats are drawn in black. Bibliographic sources are given in the main text.

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

Discussion

Diet composition

As all but one of the MOTUs identified at species or genus level from the faeces of the Mountain Long-eared Bat belonged to the Order Lepidoptera, P. macrobullaris can be referred to as a moth specialist, in accordance with previous data reported across distant areas such as Turkey and the Alps [16], [59]. In this work, however, we were able to identify the species of moths composing its diet, following the path of recent studies that have also described the prey of bats by molecular techniques, achieving parallel taxonomic resolution [22], [24], [43], [72]. The most frequent prey species in our sample was Agrotis clavis, appearing in grasslands and open habitats of mountain areas [78], [81], which widely overlap the distribution range of P. macrobullaris. It was followed by Apamea maillardi, defined as arctic-alpine [84], and the ubiquitous species Autographa gamma and A. monoglypha. Other minor prey previously reported for P. macrobullaris –e.g. coleopterans and dipterans [59]–, were either absent or infrequent in the obtained sequences.

Although Plecotus species base their diet mainly on lepidopterans, and share common and ubiquitous prey species such us Autographa gamma or Noctua pronuba [72], there are several differences among them. The number of MOTUs (54), orders detected (2), the number of prey taxa per dropping (2.75) and the dietary diversity (3.76) detected in P. macrobullaris is smaller than in the case of P. auritus and P. austriacus [72]. This fact may be the result of two factors. On the one hand, high-mountain environments are generally less diverse than lowlands [85], so the results may reflect the available food diversity. On the other hand, the method (high throughput pyrosequencing) used by Razgour et al. [72] detects almost all the sequences present in the sample, while the cloning approach used by us may fail to detect the less abundant ones as it is constrained by the number of clones used for analysis of each sample.

Alpine nature

Plecotus macrobullaris fed readily on prey found in subalpine and alpine elevational belts, which suggests that in summer they indeed forage between 1,500 and 2,500 m in the Pyrenees. In fact, all the prey species identified in this study have been reported to fly up to at least 1,800 m. Strictly alpine or subalpine prey species were detected in half of the samples. Moreover, the concurrent appearance in some faeces of moths with wide elevational ranges and those strictly alpine suggests that the former prey were also captured at high elevation. In fact, despite being more common in the lowland, some altitudinally ubiquitous moth species can often appear at higher elevations [86], [87]. Among the prey of P. macrobullaris, the moth Autographa gamma shows high migratory activity in Europe [73], [88], Apamea monoglypha has the migratory potential to cover large distances at high elevation [73], and the abundance of so-called “lowland moths” such as Noctua pronuba can be very high in alpine habitats as well [81], [86], [87]. Similarly, the altitudinally widely distributed Hadena compta is the main pollinator of the rock pink (Dianthus sylvestris), commonly found on rocky slopes at 1,000–2,800 m in the Alps [89], [90]. Additionally, most of the moth species identified as prey only appear in open areas, and more than 75% of their host plants are linked to meadows. We did not find any moth with strict forest habits, and only a handful of the generalist moths that appeared in faeces (Apamea monoglypha, Noctua pronuba, Nomophila noctuella) can be found in forests as well as in open areas. These findings are in accordance with the results obtained by Preatoni et al. [58], who observed that female Mountain Long-eared Bats avoid woods and make use of open areas and ecotones, showing a habitat use pattern more similar to P. austriacus than to P. auritus.

Although no bat species has been so far defined as purely alpine, a number of bats exploit high mountain environments [64], [91]–[93]. Alpine habitats are rich environments [94], [95] holding a great diversity of vertebrates and insects [87], [98]. Thus, moths linked to high-elevation meadows constitute a suitable resource for P. macrobullaris, which is likely most profitable when the highest meadows bloom as summer progresses. At the same time, exploiting this resource may help P. macrobullaris avoid and/or reduce interspecific competition with other moth specialists such as P. austriacus [59], [72] in lower open habitats.

However, the alpine habitat is known to exhibit extreme abiotic conditions regarding temperature, wind, solar radiation, and hypoxia [91], [96]–[98]. With low temperatures and a short vegetative growth period, some putative foraging areas of P. macrobullaris are covered with snow until June [82]. Hence, sprouting or blossoming and the matching phenology of insects are also delayed at high elevation, which, altogether, shortens the lifespan of adult phases [99]. In fact, most of the moths identified in this study only fly after June. Consequently, alpine habitats are suitable for bats for a short period that includes summer months only, as no moth species are found above 1,600 m until June at higher elevations [85]. Therefore, we should expect some differences in diet composition in spring, likely reflecting the bat's foraging at lower elevations, as do many other alpine species such as the Chamois Rupycapra sp. or the Alpine Accentor Prunella collaris [100], [101].

Hunting

All the prey species identified in this study were nocturnal flying insects, and thus gleaning behaviour cannot be inferred for P. macrobullaris, in accordance with results recently published by Ashrafi et al. [59]. Further, almost all were tympanate moths, mostly of the family Noctuidae (23 out of 28). Rather than indicating any selective behaviour, these results likely reflect the actual prey availability at high elevation, because noctuids and geometrids are almost the only lepidopteran taxa present in alpine environments, the former eight times more abundant than the latter [85]. Moths of the families Noctuidae, Geometridae, Crambidae and Pyralidae, as well as Notodontidae, are able to hear ultrasounds [28], [102], [103]. That sensitivity allows the moths to detect and escape bats, and therefore reduces the risk of predation [37], [104]. To overcome this ability, the echolocation frequency of some moth-eating bats is displayed outside the sensitive frequency range of prey [31], [105]; but unlike them, long-eared bats of the genus Plecotus studied so far use low-intensity echolocation calls to avoid detection when approaching moths [106], [107]. Plausibly, P. macrobullaris may also emit faint echolocation signals when hunting by aerial hawking.

Methodology

Mini-barcodes are a useful tool for studies with degraded DNA, such as diet analyses of bats. We could identify many prey sequences at species level, so the achieved taxonomic resolution considerably exceeds that attained by morphological analysis of faecal content, both in P. macrobullaris [16], [59] and in other moth specialists (e.g. [13], [108]).

Lepidopterans—all the sequences but one identified from faeces of P. macrobullaris—as well as coleopterans and orthopterans have yielded good identification success using DNA barcodes [109], [110]; in contrast, other taxa such as dipterans or hemipterans do not yield such effective results [111], [112]. Moths, beetles, grasshoppers, and crickets represent the bulk of food resources for many bat species, and hence, DNA barcodes will surely provide a great amount of information on these bats' ecology.

Nevertheless, there are some methodological limits worth mentioning. Although we lowered the threshold value for defining MOTUs comparing to other similar studies [23], [43], [72], anticipating less diversity among insects in high elevations, in 70% of the cases the achieved relationship among MOTUs and species was not 1∶1. This strongly contrasts the success or MOTU:species agreement observed by Razgour et al. [72]. In fact, several MOTUs (9.4%) were identified as belonging to the same species suggesting this threshold over-split actual diversity, though some MOTUs (23%) matched more than one species. It is not clear what (or if) an appropriate threshold would be. An intensive study of regional insect diversity is required to establish whether a MOTU approach can correctly approximate insect diversity in this region. In most of the cases (18 out of 23%) it was possible to identify the consumed species, because the rest of the species were not present in Western Europe. We found two MOTUs matching Agrotis simplonia, but if the first one was exclusive of this species, the second one matched other moth species as well. The only MOTU we assigned to Hadena compta matched several other species as well, some of them also present in Western Europe. As Hadena compta was the only species that matched at 100%, species-level identification was accepted. Razgour et al. [72] reported similar phenomena, showing the high similarity among sequences belonging to different species and a potential false positive using a short amplify sequence region. Consequently, it is advisable to use very conservative identification-criteria when assigning MOTUs to species.

Eight MOTUs we obtained from faeces did not match any sequence in the available reference databases. Although we can say they belong to 8 different species, we could not identify them. Nevertheless, these had likely very weak presence in the diet (4.3% of sequences), and are probably rare species limited to some high-mountain areas, which have not yet been sampled for barcoding purposes.

Correct identification of species may be conditioned by the presence of nuclear pseudogenes of the target sequence of DNA barcodes [113], [114]. Although several authors have demonstrated the presence of numts in insect taxa [114], we did not find any clear sign of nuclear pseudogenes. We took the measures proposed by Zhang & Hewitt [115] to detect numts, but no evidences such as suspicious mutations (indels, frame-shift mutations, or in-frame stop codons) or phylogenetic incongruences were detected. Even so, it is sometimes difficult to differentiate them from their mitochondrial paralogues [113], so their presence cannot be completely ruled out.

The individualised accumulation curves analysing 20 sequences per individual show that our results covered a great part of the specific diversity in diet for the sampling season. According to those curves, 100% of the estimated number of species was detected in 22 individuals. In the remaining 7 individuals, we likely did not uncover all prey content in the faeces. Lack of reliability of nonparametric species-richness estimators did not make it advisable to calculate the total species, especially when sample size was as small as in our case [116]. Therefore, the number of species that remained hidden in each dropping sample was not estimated. Nevertheless, the number of clones chosen in this study is balanced in terms of cost-efficiency, as all the estimated species were presumably detected in 72% of samples, and many more clones would have been required to produce a small increase in the number of species detected. High throughput pyrosequencing overcomes this constraint produced by clone numbers [72], but its cost/effectiveness depends on the objectives of the study, the economical constraints and the technical facilities available. Both methods are useful to detect the main bulk of species, but pyrosequencing may allow detecting also species with very low occurrence [72].

Conclusions

The molecular analysis of diet, and particularly DNA mini-barcodes, has proven a highly useful tool to obtain an accurate description of the trophic resources used by animals, and also to unveil many aspects of their foraging behaviour based on the ecology of the prey species. This may be extremely helpful when dealing with elusive and hard-to-track species such as bats, or when implementation of non-invasive techniques is necessary or advisable due to conservation or ethical constraints. In our study, results show that the Mountain Long-eared Bat P. macrobullaris is a moth specialist that forages in high mountain meadows and rocky areas during summer. However, our data are restricted to summer and we should not dismiss the possibility that P. macrobullaris forages at lower elevations and in different habitats during other seasons, especially spring, when alpine habitats are moth-impoverished. Further research will be needed to unravel such questions. Moreover, detailed inspection of the diet content revealed that P. macrobullaris feeds mainly on tympanate moths, without any data supporting it hunts by gleaning prey. Thus, we conclude that the Mountain Long-eared Bat forages by aerial hawking, emitting faint echolocation pulses to avoid early detection by moths.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

The Mountain Long-eared Bat, Plecotus macrobullaris (photo A. Alberdi).

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

(TIF)

Acknowledgments

We wish to thank all the people who assisted in sampling, and David Guixé who kindly gave bat faecal samples. We are very grateful to Matt Zeale for his help with protocols, to Elizabeth Clare for her valuable comments and assistance with the analysis of MOTUs, and to Orly Razgour for her comments to a last version of this paper. We thank very much the research group “Cell Biology in Environmental Toxicology” at the University of The Basque Country, for the logistic support, and especially Eider Bilbao, for her invaluable assistance in the laboratory. We also thank the Tiroler Landesmuseen of Austria for kindly providing the sequences of moth species. Kate Johnson proofread the English of the manuscript before submission.

Author Contributions

Conceived and designed the experiments: IG JA AA. Performed the experiments: AA OA. Analyzed the data: AA OA IG JA. Contributed reagents/materials/analysis tools: IG JA. Wrote the paper: AA IG JA. Fieldwork: AA OA IG JA.

Citation: Alberdi A, Garin I, Aizpurua O, Aihartza J (2012) The Foraging Ecology of the Mountain Long-Eared Bat Plecotus macrobullaris Revealed with DNA Mini-Barcodes. PLoS ONE 7(4): e35692. https://doi.org/10.1371/journal.pone.0035692

References

1. Stephens DW, Krebs JR (1986) Foraging Theory. Princeton, New Jersey: Princeton University Press. 247 p.DW StephensJR Krebs1986Foraging TheoryPrinceton, New JerseyPrinceton University Press247

2. Poulton EB (1929) British insectivorous bats and their prey. Proc Zool Soc London 19: 277–303.EB Poulton1929British insectivorous bats and their prey.Proc Zool Soc London19277303

3. Eisentraut M (1951) Die Ernährung der Fledermäuse (Microchiroptera). Zool Jahrb Abt Syst Geograph Biol Tiere 79: 114–177.M. Eisentraut1951Die Ernährung der Fledermäuse (Microchiroptera).Zool Jahrb Abt Syst Geograph Biol Tiere79114177

4. Ross A (1967) Ecological aspects of the food habits of insectivorous bats. Proc West Found Vert Zool 1(4): 205–264.A. Ross1967Ecological aspects of the food habits of insectivorous bats.Proc West Found Vert Zool14205264

5. Easterla DA, Whitaker JO Jr (1972) Food habits of some bats from Big Bend National Park. J Mammal 53(4): 887–890.DA EasterlaJO Whitaker Jr1972Food habits of some bats from Big Bend National Park.J Mammal534887890

6. Vaughan N (1997) The diets of British bats (Chiroptera). Mammal Rev 27(2): 77–94.N. Vaughan1997The diets of British bats (Chiroptera).Mammal Rev2727794

7. Kunz TH, Whittaker JO (1983) An evaluation of fecal analysis for determining food habits of insectivorous bats. Can J Zool 61: 1317–1321.TH KunzJO Whittaker1983An evaluation of fecal analysis for determining food habits of insectivorous bats.Can J Zool6113171321

8. Dickman CR, Huang C (1988) The reliability of fecal analysis as a method for determining the diet of insectivorous mammals. J Mammal 69: 108–113.CR DickmanC. Huang1988The reliability of fecal analysis as a method for determining the diet of insectivorous mammals.J Mammal69108113

9. Gajdosik M, Gaisler J (2004) Diet of two Eptesicus bat species in Moravia (Czech Republic). Folia Zool 53: 7–16.M. GajdosikJ. Gaisler2004Diet of two Eptesicus bat species in Moravia (Czech Republic).Folia Zool53716

10. Burles DW, Brigham RM, Ring RA, Reimchen TE (2008) Diet of two insectivorous bats, Myotis lucifugus and Myotis keenii, in relation to arthropod abundance in a temperate Pacific Northwest rainforest environment. Can J Zool 86: 1367–1375.DW BurlesRM BrighamRA RingTE Reimchen2008Diet of two insectivorous bats, Myotis lucifugus and Myotis keenii, in relation to arthropod abundance in a temperate Pacific Northwest rainforest environment.Can J Zool8613671375

11. Fukui D, Okazaki K, Maeda K (2009) Diet of three sympatric insectivorous bat species on Ishigaki Island, Japan. ESR 8: 117–128.D. FukuiK. OkazakiK. Maeda2009Diet of three sympatric insectivorous bat species on Ishigaki Island, Japan.ESR8117128

12. Almenar D, Aihartza J, Goiti U, Salsamendi E, Garin I (2008) Diet and prey selection in the trawling long-fingered bat. J Zool 274: 340–348.D. AlmenarJ. AihartzaU. GoitiE. SalsamendiI. Garin2008Diet and prey selection in the trawling long-fingered bat.J Zool274340348

13. Salsamendi E, Garin I, Almenar D, Goiti U, Napal M, et al. (2008) Diet and prey selection in Mehelyi's horseshoe bat Rhinolophus mehelyi (Chiroptera, Rhinolophidae) in the south-western Iberian Peninsula. Acta Chiropterol 10: 279–286.E. SalsamendiI. GarinD. AlmenarU. GoitiM. Napal2008Diet and prey selection in Mehelyi's horseshoe bat Rhinolophus mehelyi (Chiroptera, Rhinolophidae) in the south-western Iberian Peninsula.Acta Chiropterol10279286

14. Beck A (1995) Fecal analyses of European bat species. Myotis 32–33: 109–119.A. Beck1995Fecal analyses of European bat species.Myotis32–33109119

15. Rydell J, Natuschke G, Theiler A, Zingg PE (1996) Food habits of the barbastelle bat Barbastella barbastellus. Ecography 19: 62–66.J. RydellG. NatuschkeA. TheilerPE Zingg1996Food habits of the barbastelle bat Barbastella barbastellus.Ecography196266

16. Whitaker JO, Karatas A (2009) Food and feeding habitats of some bats from Turkey. Acta Chiropterol 11: 393–403.JO WhitakerA. Karatas2009Food and feeding habitats of some bats from Turkey.Acta Chiropterol11393403

17. Rabinowitz AR, Tuttle MD (1982) A test of the validity of two currently used methods of determining bat prey preferences. Acta Theriol 27: 283–293.AR RabinowitzMD Tuttle1982A test of the validity of two currently used methods of determining bat prey preferences.Acta Theriol27283293

18. Lee YF, McCracken GF (2005) Dietary variation of brazilian free-tailed bats links to migratory populations of pest insect. J Mammal 86: 67–76.YF LeeGF McCracken2005Dietary variation of brazilian free-tailed bats links to migratory populations of pest insect.J Mammal866776

19. Painter ML, Chambers CL, Siders M, Doucett RR, Whitaker JO Jr, et al. (2009) Diet of spotted bats (Euderma maculatum) in Arizona as indicated by faecal analysis and stable isotopes. Can J Zool 87: 865–875.ML PainterCL ChambersM. SidersRR DoucettJO Whitaker Jr2009Diet of spotted bats (Euderma maculatum) in Arizona as indicated by faecal analysis and stable isotopes.Can J Zool87865875

20. King RA, Read DS, Traugott M, Symondson WOC (2008) Molecular analysis of predation: a review of best practice for DNA-based approaches. Mol Ecol 17: 947–963.RA KingDS ReadM. TraugottWOC Symondson2008Molecular analysis of predation: a review of best practice for DNA-based approaches.Mol Ecol17947963

21. Ratnasingham S, Hebert P (2007) BOLD: The Barcode of Life Data System (www.barcodinglife.org). Molecular Ecol Notes 7: 355–364.S. RatnasinghamP. Hebert2007BOLD: The Barcode of Life Data System (www.barcodinglife.org).Molecular Ecol Notes7355364Available: http://www.barcodinglife.org. Accessed 2012 Mar 24. Available: http://www.barcodinglife.org. Accessed 2012 Mar 24.

22. Clare EL, Fraser EE, Braid HE, Fenton MB, Hebert PDN (2009) Species on the menu of a generalist predator, the eastern red bat (Lasiurus borealis): using a molecular approach to detect arthropod prey. Mol Ecol 18: 2532–2542.EL ClareEE FraserHE BraidMB FentonPDN Hebert2009Species on the menu of a generalist predator, the eastern red bat (Lasiurus borealis): using a molecular approach to detect arthropod prey.Mol Ecol1825322542

23. Bohmann K, Monadjem A, Noer CL, Rasmussen M, Zeale MRK, et al. (2011) Molecular diet analysis of two African free-tailed bats (Molossidae) using high throughput sequencing. PLoS ONE 6: e21441.K. BohmannA. MonadjemCL NoerM. RasmussenMRK Zeale2011Molecular diet analysis of two African free-tailed bats (Molossidae) using high throughput sequencing.PLoS ONE6e21441

24. Zeale MRK, Butlin RK, Barker GLA, Lees DC, Jones G (2011) Taxon-specific PCR for DNA barcoding arthropod prey in bat faeces. Mol Ecol Resour 11: 236–244.MRK ZealeRK ButlinGLA BarkerDC LeesG. Jones2011Taxon-specific PCR for DNA barcoding arthropod prey in bat faeces.Mol Ecol Resour11236244

25. Hebert PDN, Cywinska A, Ball SL, deWaard JR (2003) Biological identifications through DNA barcodes. Proc R Soc -B Biol Sci 270: 313–321.PDN HebertA. CywinskaSL BallJR deWaard2003Biological identifications through DNA barcodes.Proc R Soc -B Biol Sci270313321

26. Hebert PDN, Ratnasingham S, deWaard JR (2003) Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proc R Soc -B Biol Sci 270: S96–S99.PDN HebertS. RatnasinghamJR deWaard2003Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species.Proc R Soc -B Biol Sci270S96S99

27. Miller LA, Surlykke A (2001) How Some Insects Detect and Avoid Being Eaten by Bats: Tactics and Countertactics of Prey and Predator. Biosci 51: 570–581.LA MillerA. Surlykke2001How Some Insects Detect and Avoid Being Eaten by Bats: Tactics and Countertactics of Prey and Predator.Biosci51570581

28. Roeder KD (1967) Nerve Cells and Insect Behavior. Cambridge (MA): Harvard University Press. 238 p.KD Roeder1967Nerve Cells and Insect BehaviorCambridge (MA)Harvard University Press238

29. Werner TK (1981) Responses of nonflying moths to ultrasound: the threat of gleaning bats. Can J Zool 59: 525–529.TK Werner1981Responses of nonflying moths to ultrasound: the threat of gleaning bats.Can J Zool59525529

30. Miller LA (1984) Hearing in green lacewings and their responses to the cries of bats. In: Canard M, Séméria Y, New TR, editors. Biology of Chrysopidae. The Hague: DRW Junk Publishers. pp. 134–149.LA Miller1984Hearing in green lacewings and their responses to the cries of bats.M. CanardY. SémériaTR NewBiology of ChrysopidaeThe HagueDRW Junk Publishers134149

31. Fullard JH (1988) The tuning of moth ears. Experientia 44: 423–428.JH Fullard1988The tuning of moth ears.Experientia44423428

32. Faure PA, Fullard JH, Barclay MR (1990) The response of tympanate moths to the echolocation calls of a substrate gleaning bat, Myotis evotis. J Comp Physiol A 166: 843–849.PA FaureJH FullardMR Barclay1990The response of tympanate moths to the echolocation calls of a substrate gleaning bat, Myotis evotis.J Comp Physiol A166843849

33. Waters DA, Jones G (1995) Echolocation call structure and intensity in five species of insectivorous bats. J Exp Biol 198: 475–489.DA WatersG. Jones1995Echolocation call structure and intensity in five species of insectivorous bats.J Exp Biol198475489

34. Hartbauer M, Ofner E, Grossauer V, Siemers BM (2010) The Cercal Organ May Provide Singing Tettigoniids a Backup Sensory System for the Detection of Eavesdropping Bats. PLoS ONE 5: e12698.M. HartbauerE. OfnerV. GrossauerBM Siemers2010The Cercal Organ May Provide Singing Tettigoniids a Backup Sensory System for the Detection of Eavesdropping Bats.PLoS ONE5e12698

35. ter Hofstede HM, Kalko EKV, Fullard JH (2010) Auditory-based defence against gleaning bats in neotropical katydids (Orthoptera: Tettigoniidae). J Comp Physiol A 196: 249–358.HM ter HofstedeEKV KalkoJH Fullard2010Auditory-based defence against gleaning bats in neotropical katydids (Orthoptera: Tettigoniidae).J Comp Physiol A196249358

36. Rydell J, Arlettaz R (1994) Low-frequency echolocation enables the bat Tadarida teniotis to feed on tympanate insects. Proc R Soc London B 257: 175–178.J. RydellR. Arlettaz1994Low-frequency echolocation enables the bat Tadarida teniotis to feed on tympanate insects.Proc R Soc London B257175178

37. Pavey CR, Burwell CJ (1998) Bat predation on eared moths: a test of the allotonic frequency hypothesis. Oikos 81: 143–151.CR PaveyCJ Burwell1998Bat predation on eared moths: a test of the allotonic frequency hypothesis.Oikos81143151

38. Valdez EW, Cryan PM (2009) Food habits of the hoary bat (Lasiurus cinereus) during spring migration through New Mexico. Southwestern Nat 54: 195–200.EW ValdezPM Cryan2009Food habits of the hoary bat (Lasiurus cinereus) during spring migration through New Mexico.Southwestern Nat54195200

39. Rostovskaya MS, Zhukova DV, Illarionova AE, Ustyugova SV, Borissenko AV, et al. (2000) Insect prey of the long-eared bat Plecotus auritus (L.) (Chiroptera: Vespertilionidae) in Central Russia. Russian Entomol J 9: 185–189.MS RostovskayaDV ZhukovaAE IllarionovaSV UstyugovaAV Borissenko2000Insect prey of the long-eared bat Plecotus auritus (L.) (Chiroptera: Vespertilionidae) in Central Russia.Russian Entomol J9185189

40. Lacki M, Ladeur K (2001) Seasonal use of lepidopteran prey by Rafinesque's big-eared bats (Corynorhinus rafinesquii). Am Midland Nat 145: 213–217.M. LackiK. Ladeur2001Seasonal use of lepidopteran prey by Rafinesque's big-eared bats (Corynorhinus rafinesquii).Am Midland Nat145213217

41. Dodd LE, Lacki MJ (2007) Prey consumed by Corynorhinus townsendii ingens in the Ozark Mountain region. Acta Chiropterol 9: 451–461.LE DoddMJ Lacki2007Prey consumed by Corynorhinus townsendii ingens in the Ozark Mountain region.Acta Chiropterol9451461

42. Goerlitz HR, ter Hofstede HM, Zeale MRK, Jones G, Holderied MW (2010) An aerial-hawking bat uses stealth echolocation to counter moth hearing. Curr Biol 20: 1568–1572.HR GoerlitzHM ter HofstedeMRK ZealeG. JonesMW Holderied2010An aerial-hawking bat uses stealth echolocation to counter moth hearing.Curr Biol2015681572

43. Clare EL, Barber B, Sweeney B, Hebert PDN, Fenton MB (2011) Eating local: influences of habitat on the diet of little brown bats (Myotis lucifugus). Mol Ecol 20: 1772–1780.EL ClareB. BarberB. SweeneyPDN HebertMB Fenton2011Eating local: influences of habitat on the diet of little brown bats (Myotis lucifugus).Mol Ecol2017721780

44. Kiefer A, Veith M (2001) A new species of long-eared bat from Europe (Chiroptera: Vespertilionidae). Myotis 39: 5–16.A. KieferM. Veith2001A new species of long-eared bat from Europe (Chiroptera: Vespertilionidae).Myotis39516

45. Spitzenberger F, Haring E, Tvrtokovic N (2002) Plecotus microdontus (Mammalia, Vespertilionidae), a new bat species from Austria. Nat Croatica 11: 1–18.F. SpitzenbergerE. HaringN. Tvrtokovic2002Plecotus microdontus (Mammalia, Vespertilionidae), a new bat species from Austria.Nat Croatica11118

46. Spitzenberger F, Strelkov P, Haring E (2003) Morphology and mitochondrial DNA sequences show that Plecotus alpinus Kiefer & Veith, 2002 and Plecotus microdontus Spitzenberger, 2002 are synonyms of Plecotus macrobullaris Kuzjakin, 1965. Nat Croatica 12: 39–53.F. SpitzenbergerP. StrelkovE. Haring2003Morphology and mitochondrial DNA sequences show that Plecotus alpinus Kiefer & Veith, 2002 and Plecotus microdontus Spitzenberger, 2002 are synonyms of Plecotus macrobullaris Kuzjakin, 1965.Nat Croatica123953

47. Juste J, Ibáñez C, Muñoz D, Trujillo D, Benda P, et al. (2004) Mitochondrial phylogeography of the long-eared bats (Plecotus) in the Mediterranean Palaearctic and Atlantic Islands. Mol Phyl Evol 31: 1114–1126.J. JusteC. IbáñezD. MuñozD. TrujilloP. Benda2004Mitochondrial phylogeography of the long-eared bats (Plecotus) in the Mediterranean Palaearctic and Atlantic Islands.Mol Phyl Evol3111141126

48. Pavlinic I, Tvrtkovic N (2004) Altitudinal distribution of four Plecotus species (Mammalia, Vespertilionidae) occurring in Croatia. Nat Croatica 13: 395–401.I. PavlinicN. Tvrtkovic2004Altitudinal distribution of four Plecotus species (Mammalia, Vespertilionidae) occurring in Croatia.Nat Croatica13395401

49. Spitzenberger F, Strelkov P, Winkler H, Haring E (2006) A preliminary revision of the genus Plecotus (Chiroptera, Vespertilionidae) based on genetic and morphological results. Zool Scr 35: 187–230.F. SpitzenbergerP. StrelkovH. WinklerE. Haring2006A preliminary revision of the genus Plecotus (Chiroptera, Vespertilionidae) based on genetic and morphological results.Zool Scr35187230

50. Benda P, Georgiakakis P, Dietz C, Hanák V, Galanaki K, et al. (2008) Bats (Mammalia: Chiroptera) of the Eastern Mediterranean and Middle East. Part 7. The bat fauna of Crete, Greece. Acta Soc Zool Bohemoslovacae Bohemicae 72: 105–190.P. BendaP. GeorgiakakisC. DietzV. HanákK. Galanaki2008Bats (Mammalia: Chiroptera) of the Eastern Mediterranean and Middle East. Part 7. The bat fauna of Crete, Greece.Acta Soc Zool Bohemoslovacae Bohemicae72105190

51. Ashrafi S, Bontadina F, Kiefer A, Pavlinic I, Arlettaz R (2010) Multiple morphological characters needed for field identification of cryptic long-eared bat species around the Swiss Alps. J Zool 281: 241–248.S. AshrafiF. BontadinaA. KieferI. PavlinicR. Arlettaz2010Multiple morphological characters needed for field identification of cryptic long-eared bat species around the Swiss Alps.J Zool281241248

52. Mattei-Roesli M (2010) Situazione del genere Plecotus (Chiroptera) nel Cantone Ticino (Svizzera). Boll Soc ticinese sci nat 98: 31–34.M. Mattei-Roesli2010Situazione del genere Plecotus (Chiroptera) nel Cantone Ticino (Svizzera).Boll Soc ticinese sci nat983134

53. Garin I, García-Mudarra JL, Aihartza J, Goiti U, Juste J (2003) Presence of Plecotus macrobullaris (Chiroptera: Vespertilionidae) in the Pyrenees. Acta Chiropterol 5: 243–250.I. GarinJL García-MudarraJ. AihartzaU. GoitiJ. Juste2003Presence of Plecotus macrobullaris (Chiroptera: Vespertilionidae) in the Pyrenees.Acta Chiropterol5243250

54. Alcalde JT, Trujillo D, Artázco A, Agirre-Mendi PT (2008) Distribución y estado de conservación de los quirópteros en Aragón. Graellsia 64: 3–16.JT AlcaldeD. TrujilloA. ArtázcoPT Agirre-Mendi2008Distribución y estado de conservación de los quirópteros en Aragón.Graellsia64316

55. Boye P, Dietz M (2005) Research Report No 661: Development of good practice guidelines for woodland management for bats. Peterborough: English Nature. 90 p.P. BoyeM. Dietz2005Research Report No 661: Development of good practice guidelines for woodland management for batsPeterboroughEnglish Nature90

56. Dietrich S, Szameitat DP, Kiefer A, Schnitzler HU, Denzinger A (2006) Echolocation signals of the plecotine bat, Plecotus macrobullaris Kuzyakin, 1965. Acta Chiropterol 8: 465–475.S. DietrichDP SzameitatA. KieferHU SchnitzlerA. Denzinger2006Echolocation signals of the plecotine bat, Plecotus macrobullaris Kuzyakin, 1965.Acta Chiropterol8465475

57. Ashrafi S (2010) Resource partitioning in three cryptic, sympatric bat species (Plecotus spp.) with contrasting conservation status. S. Ashrafi2010Resource partitioning in three cryptic, sympatric bat species (Plecotus spp.) with contrasting conservation status.PhD Thesis, Universitaät Bern, Switzerland. PhD Thesis, Universitaät Bern, Switzerland.

58. Preatoni DG, Spada M, Wauters LA, Tosi G, Martinoli A (2011) Habitat use in the female Alpine long-eared bat (Plecotus macrobullaris): does breeding make the difference? Acta Chiropterologica 13(2): 355–365.DG PreatoniM. SpadaLA WautersG. TosiA. Martinoli2011Habitat use in the female Alpine long-eared bat (Plecotus macrobullaris): does breeding make the difference?Acta Chiropterologica132355365

59. Ashrafi S, Beck A, Rutishauser M, Arletta R, Bontadina F (2011) Trophic niche partitioning of cryptic species of long-eared bats in Switzerland: implications for conservation. Eur J Wildl Res 57: 843–849.S. AshrafiA. BeckM. RutishauserR. ArlettaF. Bontadina2011Trophic niche partitioning of cryptic species of long-eared bats in Switzerland: implications for conservation.Eur J Wildl Res57843849

60. Animal Behaviour Society (2006) Guidelines for the treatment of animals in behavioural research and teaching. An Behav 71: 245–253.Animal Behaviour Society2006Guidelines for the treatment of animals in behavioural research and teaching.An Behav71245253

61. Miehe G, Miehe S (2000) Comparative high mountain research on the treeline ecotone under human impact. Erdkunde 54: 34–50.G. MieheS. Miehe2000Comparative high mountain research on the treeline ecotone under human impact.Erdkunde543450

62. Olsson EGA, Austrheim G, Grenne SvN (2000) Landscape change patterns in mountains, land use and environmental diversity, Mid-Norway 1960–1993. Landscape Ecol 15: 155–170.EGA OlssonG. AustrheimSvN Grenne2000Landscape change patterns in mountains, land use and environmental diversity, Mid-Norway 1960–1993.Landscape Ecol15155170

63. Coldea G (2003) The alpine flora and vegetation of the South-Eastern Carpathians. In: Nagy L, Grabherr G, Körner C, Thompson DBA, editors. Alpine Biodiversity in Europe. Berlin: Springer. pp. 65–72.G. Coldea2003The alpine flora and vegetation of the South-Eastern Carpathians.L. NagyG. GrabherrC. KörnerDBA ThompsonAlpine Biodiversity in EuropeBerlinSpringer6572

64. García-González R (2008) Management of Natura 2000 habitats. 6170 Alpine and subalpine calcareous grasslands. European Commission. Environment DG. R. García-González2008Management of Natura 2000 habitats. 6170 Alpine and subalpine calcareous grasslands. European Commission. Environment DG.Available: http://ec.europa.eu/environment/nature/natura2000/management/habitats/models_en.htm. Accessed 2012 Apr 4. Available: http://ec.europa.eu/environment/nature/natura2000/management/habitats/models_en.htm. Accessed 2012 Apr 4.

65. Whitaker JO, Neefus C Jr, Kunz TH (1996) Dietary variation in the Mexican free-tailed bat (Tadarida brasiliensis mexicana). J Mammal 77: 716–724.JO WhitakerC. Neefus JrTH Kunz1996Dietary variation in the Mexican free-tailed bat (Tadarida brasiliensis mexicana).J Mammal77716724

66. Hebert PDN, Ratnasingham S, Dooh R (2004) Barcodes of Life. PDN HebertS. RatnasinghamR. Dooh2004Barcodes of Life.Available: http://www.barcodinglife.com/. Accessed 2012 Mar 24. Available: http://www.barcodinglife.com/. Accessed 2012 Mar 24.

67. Deagle BE, Eveson JP, Jarman SN (2006) Quantification of damage in DNA recovered from highly degraded samples – a case study on DNA in faces. Front Zool 3: 11.BE DeagleJP EvesonSN Jarman2006Quantification of damage in DNA recovered from highly degraded samples – a case study on DNA in faces.Front Zool311

68. Hajibabaei M, Smith M, Janzen D, Rodriguez J, Whitfield J, et al. (2006) A minimalist barcode can identify a specimen whose DNA is degraded. Mol Ecol Notes 6: 959–964.M. HajibabaeiM. SmithD. JanzenJ. RodriguezJ. Whitfield2006A minimalist barcode can identify a specimen whose DNA is degraded.Mol Ecol Notes6959964

69. Meusnier I, Singer G, Landry JF, Hickey D, Hebert P, et al. (2008) A universal DNA mini-barcode for biodiversity analysis. BMC Genomics 9: 214.I. MeusnierG. SingerJF LandryD. HickeyP. Hebert2008A universal DNA mini-barcode for biodiversity analysis.BMC Genomics9214

70. Hall TA (1999) BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser 41: 95–98.TA Hall1999BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT.Nucleic Acids Symp Ser419598

71. Jones M, Ghoorah A, Blaxter M (2011) jMOTU and taxonerator: turning DNA barcode sequences into annotated operational taxonomic units. PLos ONE 6: e19259.M. JonesA. GhoorahM. Blaxter2011jMOTU and taxonerator: turning DNA barcode sequences into annotated operational taxonomic units.PLos ONE6e19259

72. Razgour O, Clare EL, Zeale MRK, Hanmer J, Schnell IB, et al. (2011) High-throughput sequencing offers insight into mechanisms of resource partitioning in cryptic bat species. Ecol Evol 1: 556–570.O. RazgourEL ClareMRK ZealeJ. HanmerIB Schnell2011High-throughput sequencing offers insight into mechanisms of resource partitioning in cryptic bat species.Ecol Evol1556570

73. Colwell RK (2006) EstimateS: Statistical estimation of species richness and shared species from samples. Version 8. RK Colwell2006EstimateS: Statistical estimation of species richness and shared species from samples. Version 8.Available: http://viceroy.eeb.uconn.edu/EstimateS. Accessed 2012 Mar 24. Available: http://viceroy.eeb.uconn.edu/EstimateS. Accessed 2012 Mar 24.

74. Wood C, Reynolds D, Wells P, Barlow J, Woiwod I, et al. (2009) Flight periodicity and the vertical distribution of high-altitude moth migration over southern Britain. Bull Entomol Res 99: 525–535.C. WoodD. ReynoldsP. WellsJ. BarlowI. Woiwod2009Flight periodicity and the vertical distribution of high-altitude moth migration over southern Britain.Bull Entomol Res99525535

75. Pérez De-Gregorio JJ, Muñoz J, Rondós M (2001) Atlas fotográfico de los lepidópteros macroheteróceros íbero-baleares 2. Barcelona: Argania Editio. 210 p.JJ Pérez De-GregorioJ. MuñozM. Rondós2001Atlas fotográfico de los lepidópteros macroheteróceros íbero-baleares 2BarcelonaArgania Editio210

76. Dodok I (2003) Noctuidae (Lepidoptera) of the Uzice region (Western Serbia). Acta Entomol Serbica 8: 1–13.I. Dodok2003Noctuidae (Lepidoptera) of the Uzice region (Western Serbia).Acta Entomol Serbica8113

77. Waring P, Townsend M (2003) Field Guide to the Moths of Great Britain and Ireland. Rotherwick, Hampshire: British Wildlife Publishing. 432 p.P. WaringM. Townsend2003Field Guide to the Moths of Great Britain and IrelandRotherwick, HampshireBritish Wildlife Publishing432

78. Robineau R (2007) Guide des papillons nocturnes de France: plus de 1620 espèces décrites et illustrées. Paris: Delachaux et Niestlé. 287 p.R. Robineau2007Guide des papillons nocturnes de France: plus de 1620 espèces décrites et illustréesParisDelachaux et Niestlé287

79. Dinca V, Székely L, Kovács S, Kovács Z, Vila R (2008) Pyrgus andromedae (Wallengren, 1853) (Hesperiidae) in the Romanian Carpathians. Nota Lepidopterol 31: 263–272.V. DincaL. SzékelyS. KovácsZ. KovácsR. Vila2008Pyrgus andromedae (Wallengren, 1853) (Hesperiidae) in the Romanian Carpathians.Nota Lepidopterol31263272

80. Ortiz A, Garre M, Guerrero JJ, Rubio RM, Calle JA (2009) Nuevas citas de la Familia Noctuidae en Cuenca (España) (Lepidoptera: Noctuidae). SHILAP Rev Lepidopterol 37: 105–114.A. OrtizM. GarreJJ GuerreroRM RubioJA Calle2009Nuevas citas de la Familia Noctuidae en Cuenca (España) (Lepidoptera: Noctuidae).SHILAP Rev Lepidopterol37105114

81. Scalercio S (2009) On top of a Mediterranean massif: climate change and conservation of orophilous moths at the southern boundary of their range (Lepidoptera: Macroheterocera). Eur J Entomol 106: 231–239.S. Scalercio2009On top of a Mediterranean massif: climate change and conservation of orophilous moths at the southern boundary of their range (Lepidoptera: Macroheterocera).Eur J Entomol106231239

82. Körner C (1999) Alpine Plant Life. Functional Plant Ecology of High Mountain Ecosystems. Berlin: Springer-Verlag. 338 p.C. Körner1999Alpine Plant Life. Functional Plant Ecology of High Mountain EcosystemsBerlinSpringer-Verlag338

83. EPA (2001) National Land Cover Data: Land Cover Class Definitions. EPA2001National Land Cover Data: Land Cover Class Definitions.Available: http://www.epa.gov/mrlc/nlcd-2001.html. Accessed 2012 March 24. Available: http://www.epa.gov/mrlc/nlcd-2001.html. Accessed 2012 March 24.

84. Varga ZS, Schmitt T (2008) Types of oreal and oreotundral disjunctions in the western Palearctic. Biol J Linn Soc 93: 415–430.ZS VargaT. Schmitt2008Types of oreal and oreotundral disjunctions in the western Palearctic.Biol J Linn Soc93415430

85. Beck J, Altermatt F, Hagmann R, Lang S (2010) Seasonality in the altitude–diversity pattern of Alpine moths. Basic Appl Ecol 11: 714–722.J. BeckF. AltermattR. HagmannS. Lang2010Seasonality in the altitude–diversity pattern of Alpine moths.Basic Appl Ecol11714722

86. Spitzer K, Jaros J, Vanek J (2003) Migrations of Moths in the Krkonoše Mts. (The Giant Mts.). Opera Corcontica 40: 275–286.K. SpitzerJ. JarosJ. Vanek2003Migrations of Moths in the Krkonoše Mts. (The Giant Mts.).Opera Corcontica40275286

87. Schnepf J (2010) Diversity of diurnal Lepidoptera across an elevational gradient above the timberline: a case study from the Austrian Alps. J. Schnepf2010Diversity of diurnal Lepidoptera across an elevational gradient above the timberline: a case study from the Austrian Alps.PhD Thesis, Universität Wien, Austria. PhD Thesis, Universität Wien, Austria.

88. Chapman JW, Reynolds DR, Mouritsen H, Hill JK, Riley JR, et al. (2008) Wind selection and drift compensation optimise migratory pathways in a high-flying moth. Curr Biol 18: 514–518.JW ChapmanDR ReynoldsH. MouritsenJK HillJR Riley2008Wind selection and drift compensation optimise migratory pathways in a high-flying moth.Curr Biol18514518

89. Collin C, Shykoff J (2003) Outcrossing rates in the gynomonoecious-gynodioecious species Dianthus sylvestris (Caryophyllaceae). Am J Bot 90: 579–585.C. CollinJ. Shykoff2003Outcrossing rates in the gynomonoecious-gynodioecious species Dianthus sylvestris (Caryophyllaceae).Am J Bot90579585

90. Collin C, Shykoff J (2010) Flowering phenology and female fitness: impact of a pre-dispersal seed predator on a sexually polymorphic species. Plant Ecol 206: 1–13.C. CollinJ. Shykoff2010Flowering phenology and female fitness: impact of a pre-dispersal seed predator on a sexually polymorphic species.Plant Ecol206113

91. Martin K (2001) Wildlife in alpine and sub-alpine habitats. In: Johnson DH, O'Neil TA, editors. Wildlife-Habitat Relationships in Oregon and Washington. Corvallis (OR): Oregon State University Press. pp. 285–310.K. Martin2001Wildlife in alpine and sub-alpine habitatsDH JohnsonTA O'NeilWildlife-Habitat Relationships in Oregon and WashingtonCorvallis (OR)Oregon State University Press285310

92. Nappi A (2002) Vertical distribution of the snow vole Chionomys nivalis (Martins, 1842) (Rodentia, Arvicolidae) in Italy. Hystrix 13: 45–52.A. Nappi2002Vertical distribution of the snow vole Chionomys nivalis (Martins, 1842) (Rodentia, Arvicolidae) in Italy.Hystrix134552

93. Huemer P (2005) Die Kanisfluh im Bregenzerwald (Vorarlberg), ein «Hot Spot» der Biodiversität für Schmetterlinge (Lepidoptera). Vorarlberger Naturschau 16: 9–92.P. Huemer2005Die Kanisfluh im Bregenzerwald (Vorarlberg), ein «Hot Spot» der Biodiversität für Schmetterlinge (Lepidoptera).Vorarlberger Naturschau16992

94. Grabherr G, Gottfried M, Pauli H (2000) GLORIA: A Global Observation Research Initiative in Alpine Environments. Mountain Res Develop 20: 190–192.G. GrabherrM. GottfriedH. Pauli2000GLORIA: A Global Observation Research Initiative in Alpine Environments.Mountain Res Develop20190192

95. Väre H, Lampinen R, Humphries C, Williams P (2003) Taxonomic diversity of vascular plants in the European alpine areas. In: Nagy L, Grabherr G, Körner C, Thompson DBA, editors. Alpine Biodiversity in Europe. Berlin: Springer. pp. 133–148.H. VäreR. LampinenC. HumphriesP. Williams2003Taxonomic diversity of vascular plants in the European alpine areas.L. NagyG. GrabherrC. KörnerDBA ThompsonAlpine Biodiversity in EuropeBerlinSpringer133148

96. Davidson P, Wilson S, Bears H, Camfield A, Martin K (2007) Birds in High Places. Birdwatch Can 41: 4–7.P. DavidsonS. WilsonH. BearsA. CamfieldK. Martin2007Birds in High Places.Birdwatch Can4147

97. Wiebe KL, Martin K (1998) Age specific patterns of reproduction in white-tailed and willow ptarmigan Lagopus leucurus and L. lagopus. Ibis 140: 14–24.KL WiebeK. Martin1998Age specific patterns of reproduction in white-tailed and willow ptarmigan Lagopus leucurus and L. lagopus.Ibis1401424

98. Hodkinson I (2005) Terrestrial insects along elevation gradients: species and community responses to altitude. Biol Rev 80: 489–513.I. Hodkinson2005Terrestrial insects along elevation gradients: species and community responses to altitude.Biol Rev80489513

99. Shapiro AM (1975) The temporal component of butterfly species diversity. 181–195. In: Cody MJ, Diamond JM, editors. Ecology and evolution of communities. Cambridge: Harvard University Press. pp. 181–195.AM Shapiro1975The temporal component of butterfly species diversity. 181–195.MJ CodyJM DiamondEcology and evolution of communitiesCambridgeHarvard University Press181195

100. Pepin D, Menaut P, Desneux L, Cargnelutti B (1992) Seasonal changes in the use of space by isards (Rupicapra pyrenaica) in a protected area. In: Spitz F, Janeau G, Gonzalez G, Aulagnier S, editors. pp. 327–330.D. PepinP. MenautL. DesneuxB. Cargnelutti1992Seasonal changes in the use of space by isards (Rupicapra pyrenaica) in a protected area.F. SpitzG. JaneauG. GonzalezS. Aulagnier327330Ongules/Ungulates 91, Toulouse, France. Ongules/Ungulates 91, Toulouse, France.

101. Henry PY (2011) Differential migration in the polygynandrous Alpine accentor Prunella collaris. Bird Study. PY Henry2011Differential migration in the polygynandrous Alpine accentor Prunella collaris.Bird StudyDOI:10.1080/00063657.2010.543948. DOI:10.1080/00063657.2010.543948.

102. Roeder KD, Treat AE (1957) Ultrasonic reception by the tympanic organ of noctuid moths. J Exp Zool 134: 127–157.KD RoederAE Treat1957Ultrasonic reception by the tympanic organ of noctuid moths.J Exp Zool134127157

103. Scoble MJ (1992) The Lepidoptera. Form, Function and Diversity. Oxford: Oxford University Press. 416 p.MJ Scoble1992The Lepidoptera. Form, Function and DiversityOxfordOxford University Press416

104. Acharya L, Fenton MB (1992) Echolocation behaviour of vespertilionid bats (Lasiurus cinereus and L. borealis) attacking airborne targets including arctiid moths. Can J Zool 70: 1292–1298.L. AcharyaMB Fenton1992Echolocation behaviour of vespertilionid bats (Lasiurus cinereus and L. borealis) attacking airborne targets including arctiid moths.Can J Zool7012921298

105. Bogdanowicz W, Fenton MB, Daleszczyk K (1999) The relationships between echolocation calls, morphology and diet in insectivorous bats. Journal of Zoology, London 247: 381–393.W. BogdanowiczMB FentonK. Daleszczyk1999The relationships between echolocation calls, morphology and diet in insectivorous bats.Journal of Zoology, London247381393

106. Anderson ME, Racey PA (1991) Feeding behavior of captive brown long-eared bats, Plecotus auritus. An Behav 42: 489–493.ME AndersonPA Racey1991Feeding behavior of captive brown long-eared bats, Plecotus auritus.An Behav42489493

107. Waters DA, Jones G (1996) The peripheral auditory characteristics of noctuid moths: responses to the search-phase echolocation calls of bats. J Exp Biol 199: 847–856.DA WatersG. Jones1996The peripheral auditory characteristics of noctuid moths: responses to the search-phase echolocation calls of bats.J Exp Biol199847856

108. Goiti U, Aihartza J, Garin I (2004) Diet and prey selection in the Mediterranean horseshoe bat Rhinolophus euryale (Chiroptera, Rhinolophidae) during the pre-breeding season. Mammalia 68: 397–402.U. GoitiJ. AihartzaI. Garin2004Diet and prey selection in the Mediterranean horseshoe bat Rhinolophus euryale (Chiroptera, Rhinolophidae) during the pre-breeding season.Mammalia68397402

109. Hajibabaei M, Janzen DH, Burns JM, Hallwachs W, Hebert PDN (2006) DNA barcodes distinguish species of tropical Lepidoptera. Proc Natl Acad Sci U S A 103: 968–971.M. HajibabaeiDH JanzenJM BurnsW. HallwachsPDN Hebert2006DNA barcodes distinguish species of tropical Lepidoptera.Proc Natl Acad Sci U S A103968971

110. Raupach MJ, Astrin JJ, Hannig K, Peters MK, Stoeckle MY, et al. (2010) Molecular species identification of Central European ground beetles (Coleoptera: Carabidae) using nuclear rDNA expansion segments and DNA barcodes. Front Zool 7: 26.MJ RaupachJJ AstrinK. HannigMK PetersMY Stoeckle2010Molecular species identification of Central European ground beetles (Coleoptera: Carabidae) using nuclear rDNA expansion segments and DNA barcodes.Front Zool726

111. Meier R, Shiyang K, Vaidya G, Ng PKL (2006) DNA barcoding and taxonomy in Diptera: a tale of high intraspecific variability and low identification success. Syst Biol 55: 715–728.R. MeierK. ShiyangG. VaidyaPKL Ng2006DNA barcoding and taxonomy in Diptera: a tale of high intraspecific variability and low identification success.Syst Biol55715728

112. Virgilio M, Backeljau T, Nevado B, De Meyer M (2010) Comparative performances of DNA barcoding across insect orders. BMC Bioinformatics 11: 206.M. VirgilioT. BackeljauB. NevadoM. De Meyer2010Comparative performances of DNA barcoding across insect orders.BMC Bioinformatics11206

113. Song H, Buhay JE, Whiting MF, Crandall KA (2008) Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified. Proc Natl Acad Sci U S A 105: 13486–13491.H. SongJE BuhayMF WhitingKA Crandall2008Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified.Proc Natl Acad Sci U S A1051348613491

114. Moulton MJ, Song H, Whiting M (2010) Assessing the effects of primer specificity on eliminating numt coamplification in DNA barcoding: a case study from Orthoptera (Arthropoda: Insecta). Mol Ecol Resour 10: 615–627.MJ MoultonH. SongM. Whiting2010Assessing the effects of primer specificity on eliminating numt coamplification in DNA barcoding: a case study from Orthoptera (Arthropoda: Insecta).Mol Ecol Resour10615627

115. Zhang DX, Hewitt GM (1996) Highly conserved nuclear copies of the mitochondrial control region in the desert locust Schistocerca gregaria: some implications for population studies. Mol Ecol 5: 295–300.DX ZhangGM Hewitt1996Highly conserved nuclear copies of the mitochondrial control region in the desert locust Schistocerca gregaria: some implications for population studies.Mol Ecol5295300

116. Chiarucci A, Enright NJ, Perry GLW, Miller BP, Lamont BB (2003) Performance of nonparametric species richness estimators in a high diversity plant community. Divers Distrib 9: 283–295.A. ChiarucciNJ EnrightGLW PerryBP MillerBB Lamont2003Performance of nonparametric species richness estimators in a high diversity plant community.Divers Distrib9283295