Skip to main content

Second generation sequencing and morphological faecal analysis reveal unexpected foraging behaviour by Myotis nattereri (Chiroptera, Vespertilionidae) in winter

Abstract

Background

Temperate winters produce extreme energetic challenges for small insectivorous mammals. Some bat species inhabiting locations with mild temperate winters forage during brief inter-torpor normothermic periods of activity. However, the winter diet of bats in mild temperate locations is studied infrequently. Although microscopic analyses of faeces have traditionally been used to characterise bat diet, recently the coupling of PCR with second generation sequencing has offered the potential to further advance our understanding of animal dietary composition and foraging behaviour by allowing identification of a much greater proportion of prey items often with increased taxonomic resolution. We used morphological analysis and Illumina-based second generation sequencing to study the winter diet of Natterer’s bat (Myotis nattereri) and compared the results obtained from these two approaches. For the first time, we demonstrate the applicability of the Illumina MiSeq platform as a data generation source for bat dietary analyses.

Results

Faecal pellets collected from a hibernation site in southern England during two winters (December-March 2009–10 and 2010–11), indicated that M. nattereri forages throughout winter at least in a location with a mild winter climate. Through morphological analysis, arthropod fragments from seven taxonomic orders were identified. A high proportion of these was non-volant (67.9% of faecal pellets) and unexpectedly included many lepidopteran larvae. Molecular analysis identified 43 prey species from six taxonomic orders and confirmed the frequent presence of lepidopteran species that overwinter as larvae.

Conclusions

The winter diet of M. nattereri is substantially different from other times of the year confirming that this species has a wide and adaptable dietary niche. Comparison of DNA derived from the prey to an extensive reference dataset of potential prey barcode sequences permitted fine scale taxonomic resolution of prey species. The high occurrence of non-volant prey suggests that gleaning allows prey capture at low ambient temperatures when the abundance of flying insects may be substantially reduced. Interesting questions arise as to how M. nattereri might successfully locate and capture some of the non-volant prey species encountered in its faeces. The consumption of lepidopteran larvae such as cutworms suggests that M. nattereri eats agricultural pest species.

Background

The onset of winter with associated decreases in ambient temperature leads to increased energetic challenges for mammals in many temperate regions. Small insectivorous mammals face two key challenges. First, their high surface area-to-volume ratio renders them vulnerable to high levels of heat loss and energy expenditure at low temperatures [1, 2]. Second, they face a reduction in available prey because reduced temperatures reduce arthropod abundance and activity [3, 4]. Heterothermic mammals such as bats cope with these challenges using hibernation, that is, by lowering their body temperature, metabolic rate and therefore energetic requirements over extended periods of time [2].

All hibernators arouse periodically, probably to meet behavioural or physiological requirements that are brought about by their reduced metabolism [5]. Studies on hibernating bats in cold continental locations (where winter temperatures rarely rise above freezing) show that although arousals do occur, they are relatively infrequent [68]. In contrast, bats arouse from hibernation more frequently in locations where winters are mild [913]. The relatively high energetic costs of winter arousal [2, 14] will dictate that periods of euthermia are kept to the minimum necessary.

To obtain a net energetic gain from foraging, insect prey need to be available in sufficient quantities to offset the costs of activity. Insect flight activity declines dramatically as ambient temperature drops below 6-10°C [3, 15], but some arthropods remain active on the ground or amongst foliage at temperatures below those suitable for flight [3, 16, 17]. The sensory ecology of a predator, and therefore the range of foraging techniques that it may utilise, will play an important role in its ability to exploit potential opportunities for feeding in winter [17, 18]. For example, foraging may remain profitable for bats that glean prey from vegetation (gleaning bats) even during nights when ambient temperatures drop below the threshold for insect flight [13].

Evidence for winter foraging

In temperate locations where winter temperatures are relatively mild, bat activity outside of the hibernacula occurs frequently during winter. For instance, in New South Wales, Australia, tree-roosting bats, Nyctophilus geoffroyi and N. gouldi, were active for longer on warmer winter nights suggesting occasional winter foraging [11]. Similarly, in the south of England, feeding buzzes (increase in echolocation pulse rate as a bat attempts to catch prey) were recorded during winter from foraging pipistrelles (Pipistrellus spp.) [19] and noctules (Nyctalus noctula) [20], and a single Natterer’s bat (Myotis nattereri) was radio tracked foraging outside of the hibernaculum over ten winter nights in the south of England [13]. Furthermore, studies on the winter diet of rhinolophid bats confirmed that these species feed throughout the winter in the United Kingdom [21, 22] and continental Europe [23]. The majority of arousals [10, 13] and peaks in subsequent activity [24, 25] are synchronised around sunset, when insects will be most abundant [11, 26]. Activity is also more likely on milder nights when insect abundance and activity is greater [27]. In contrast, arousals [8] and subsequent activity [28] in cold continental climates do not appear to be synchronised with zeitgebers such as photoperiod or temperature.

Until now, the winter diet of vespertilionid bats in locations with mild temperate climate has been largely overlooked, but studies of the activity and diet of winter flying northern bats (Myotis septentrionalis) and little brown bats (Myotis lucifugus) in Indiana, USA, demonstrated that these bats did not feed during the winter [29, 30]. Dehydration or the selection of more suitable roosting location may be a driving force behind the winter flights of bats in such cold continental regions [31].

Studies of trophic ecology

Traditionally, microscopic analysis has been used to identify prey insects remains in faecal pellets produced by insectivorous bats [32]. This technique is complicated by the bats’ tendency to rapidly and thoroughly chew their prey and spit out hard parts, resulting in limited taxonomic resolution of the diet [33]. Consequently, during the last five years, molecular diet analyses, which can yield great taxonomic resolution of prey if the reference database is comprehensive, have increased in popularity. Initially, such molecular studies of bat diet were typically conducted by PCR amplification of prey DNA followed by molecular cloning so as to enable subsequent Sanger sequencing of the composite prey sequences (e.g. [34, 35]). Subsequently, however, the introduction of second generation sequencing approaches such as the GS FLX [36, 37] and Ion-torrent PGM [3840] increased sequencing power while decreasing the cost of metabarcoding, and were therefore rapidly adopted for bat diet analyses. To our knowledge, at the time of writing, the current study is the first time the Illumina MiSeq platform is used as a data generation source for bat dietary analyses.

In this study, we used both Illumina MiSeq-based second generation sequencing and microscope-based morphological analyses to identify the constituent prey in faecal pellets from Myotis nattereri sensu stricto [41] collected during two winters. The overall aim was to investigate the winter diet of M. nattereri at this hibernation site, while comparing the scope of the two analytical techniques. We find that quantitative representation of diet scored by morphological and molecular analyses are significantly positively correlated at the order level, and that M. nattereri feed on a high proportion of non-flying invertebrates in winter, indicating a gleaning foraging strategy. We discuss the potential challenges that M. nattereri may face in detecting non-flying invertebrate prey items.

Results

Amplified Cyt b sequences were successfully obtained from 73 (81%) of the 90 extracted faecal samples. A match to bat species (using the criteria set out in Methods) was obtained for 71 (79%) of these, with 62 identified as Myotis nattereri sensu stricto and the remainder as Myotis daubentonii (GenBank accession numbers: KJ719267- KJ719283). Samples that yielded no positive match to a single bat species were discarded. All subsequent DNA and morphological analyses were conducted using only the M. nattereri samples. Of these 62 analysed droppings, 25 (40%) were collected in December, 21 (34%) in January and February and 16 (26%) during March. Each period included samples collected during both winters (2009–10 and 2010–11) of the study.

Morphological analysis

Microscopic analysis revealed arthropod fragments from seven taxonomic orders (Table 1). When all samples were pooled (n = 62) representing both winters, Lepidoptera was the most frequently occurring order measured by % occurrence (%O) (Table 1), though a high proportion of these remains were attributed to lepidopteran larvae. Araneae was the next most frequently encountered order in pellets, followed by Isopoda and Diptera. Coleoptera, Hymenoptera and Hemiptera were less frequently encountered over the entire winter period. By proportion of pooled data (%F), 67.9% of taxa recorded were either non-volant invertebrates or identified as being in a non-volant stage of their life cycle (i.e. lepidopteran larvae).

Table 1 Percentage of prey orders encountered within M.nattereri faecal pellets

The percentage frequency (%F) of adult Lepidoptera was highest during January and February (26%), while the proportion of lepidopteran larvae increased as winter progressed from 18% in December, to 26% in January/February up to 34% in March. The percentage frequency (%F) of Diptera was highest during December (20%), and then dropped to its lowest value during January and February (2%). The proportion (%F) of non-volant invertebrates was at its highest during January/February (70%).

Molecular analysis

A total of 99 Operational Taxonomical Units (OTUs) were identified in the 62 Myotis nattereri faecal droppings sequenced on the Illumina MiSeq platform. The assigned OTUs were divided into three main groups: Arachnida, Insecta, and Malacostraca; which contained 297, 968, and 897 collapsed haplotypes respectively. Furthermore, 1,834 collapsed haplotypes could be assigned to species level. The number of prey taxa (OTUs) per faecal sample ranged between one and 21 with a mean of 8.1 (SD ± 3.9). We successfully identified arthropod DNA from six taxonomic orders (Table 1), and were able to assign OTUs to 43 arthropod species (and an additional six OTUs to one of two species) from 24 families (Table 2). DNA analysis failed to identify only one order (Hemiptera) that had been identified through morphological analysis.

Table 2 List of prey identified within M.natterei faecal pellets through high-throughput sequencing

The percentage frequency of arthropod orders shows a similar seasonal trend to the morphological results. The proportion (%F) of Lepidoptera remained relatively stable over the winter (December: 30%; January/February: 32%; March: 33%), while the proportion (%F) of arachnids decreased in early spring (December: 23%; January/February: 20%; March: 6%). The proportion (%F) of Diptera decreased from December (18%) to 13% in January/February then increased again during March (18%).

Isopods comprised a substantial part of the diet throughout the sampling period and were recorded in all but three (n = 59) of the faecal samples. The large numbers of isopods found conflicts with the results from the morphological analysis where this order was recorded in less than one-third (n = 20) of all faecal samples. Coleoptera, Hymenoptera and Hemiptera formed a very small component part of the winter diet of M. nattereri making up <4% of morphological and <6% of molecular dietary composition.

The molecular results show that the percent frequency (%F) of lepidopteran larvae identified was greater than that of adults in all winter months (Figure 1a). Morphological findings (Figure 1b) suggest a higher proportion of larvae in December and March compared to adults.

Figure 1
figure 1

Percent frequency of lepidopteran adults and larvae in the diet of M. nattereri during three winter periods using data from 2009–10 and 2010–11 combined. a) Results from molecular analysis, whereby adults or larvae are determined by the proportion of OTUs allotted to either adult or larval stages according to the known phenology of species encountered. b) Results from morphological analysis where adult or larvae are determined from lepidopteran body fragments recovered from faecal pellets.

The proportion of Lepidoptera (%F) by species (Figure 2) shows that the putative larval sequences identified were dominated by the square spot rustic, Xanthia xanthographa, and the large yellow underwing, Noctua pronuba.

Figure 2
figure 2

The proportion of Lepidoptera (%F) by species in the winter diet of M. nattereri from DNA analysis of faecal pellets collected from within Greywell Tunnel in December-March 2009–10 and 2010–11 combined. Open bars indicate species predicted to overwinter as larvae and black bars indicate those species that overwinter as adults [44, 45].

We compared the percentage frequency (%F) of prey orders in faecal pellets between the morphological and molecular data (Figure 3). Although Isopoda were proportionately more frequent in the molecular analysis, there was nevertheless a significant positive correlation in the occurrence of prey orders described by the two methods (r = 0.816, n = 6, P <0.05).

Figure 3
figure 3

Comparison between percentage frequency (%F) of prey orders in faecal pellets found by using morphological and molecular analyses (OTUs identified to species using BOLD then assigned to individual pellets). Labels identify the taxonomic orders as follows, a) Hymenoptera, b) Coleoptera, c) Diptera, d) Araneae, e) Isopoda and f) Lepidoptera.

Discussion

Using second generation sequencing approaches to assess diet can be a very powerful tool [36, 4649], also see review by Pompanon et al. [50]. These approaches can detect prey items that are not visible under a microscope, and if the reference database is comprehensive, they can be used to identify prey species at a much finer taxonomic level than morphological analysis [34]. If the reference database is not comprehensive, however, the difference between morphological and molecular analyses in the taxonomic level of prey identification becomes less pronounced and as with morphological analysis, the results from the molecular analysis will have to be presented at order or family level while species numbers will have to be inferred by OTUs [36]. Despite the data generation power of second generation sequencing and the ability to present results at a much higher taxonomic resolution, it does have its shortcomings compared to microscopic analysis. For example, molecular approaches are unable to identify different life stages of prey, as the DNA of e.g. a moth will be the same no matter if the moth is a larva or an adult. Also, identification of the noctuid parasite, the tachinid fly (Ramonda spathula) (Table 3) in 11.3% of analysed pellets shows that second generation sequencing can even detect internal prey parasites. For both parasites and secondary predation, molecular (and to a lesser extent morphological) analyses of diet can make it difficult to discern between prey, which were eaten directly and prey that were consumed because they were in turn eaten by prey [51] or because they parasitised the prey.

Table 3 Phenology of lepidopteran prey species found in the winter diet of M. nattereri

Skews in amplification of multi-template extracts can arise due to PCR stochasticity and primer preference [50, 53]. In this study, generic arthropod primers (ZBJ-ArtF1c and ZBJ-ArtR2c) amplifying a 157 bp mini-barcode fragment were used [34]. When designed and tested, feeding trials and comparisons between molecular and morphological analyses suggested that the primers did not suffer from significant amplification biases and showed that they are able to co-amplify insect DNA in mixed templates [34]. Subsequently, the primers have been used in molecular bat diet studies where they have amplified insect DNA across orders e.g. [34, 36, 37, 54, 55]. In the current study, each extract was amplified twice with a unique barcode combination and during subsequent analyses only sequences appearing in both PCRs were computationally retained. As PCR and sequencing errors arise independently, this approach was chosen in order to minimise the occurrence of artefacts in the final results. This conservative approach has the drawback that sequences originating from rare templates which might not be amplified in both PCRs due to PCR stochasticity [53, 55] or templates which the primers are less likely to amplify due to primer preferences, might be discarded alongside sequences with PCR and/or sequencing errors. Therefore, many metabarcoding studies take the opposite approach when dealing with mixed templates. That is, after amplifying each extract a given number of times, all sequences originating from these amplifications are retained e.g. [56, 57]. Future metabarcoding bat diet studies might benefit from a combination of these two approaches, i.e. by performing a study-dependent number (>2) of PCR replicates per extract, based on the number of replicates in which a sequence should appear, in order to be retained can be determined. This might help balance the trade-offs between detecting true biological sequences from low-templates and eliminating PCR artefacts.

Amplified Cyt b sequences were successfully obtained from 81% of extracted faecal samples. This is lower than that reported by Boston et al. [58] who successfully amplified a larger fragment size (ca. 1200) with a greater degree of success 89%. Our lower amplification success of a shorter fragment may be due to differences in extraction protocol or higher DNA degradation due to non-optimal storage conditions.

Although a correlation was apparent between the percentage frequency of prey orders in molecular and morphological analyses, some major differences in the representation of prey as revealed by these methods occurred (Table 1). For example, Isopoda were more highly represented in the molecular analysis, and arachnids and lepidopterans were relatively more represented in the morphological analysis. Both molecular and morphological methods for dietary analysis have biases [32]. Controlled feeding experiments in which bats are fed measured amounts of known prey types and their faecal droppings are analysed with both methods would be informative for better understanding these biases.

We found 99 OTUs in the bat faecal pellets of which we assigned 72 to an insect species (Table 2). Since the microscopic analysis was not able to go beyond family level, we cannot compare the molecular and morphological results at either genus or species level and without having performed controlled feeding experiments, we can only assume that the obtained number of OTUs gives an accurate representation of the insect species present in the analysed faecal pellets.

Previous studies on the diet of M .nattereri

Previous dietary studies (See [59] for a review of studies prior to 1995; [18, 60, 61]) revealed that M. nattereri has a wide dietary niche breadth, taking prey from a diverse range of arthropod orders [59]. M. nattereri eats predominantly diurnal dipteran species [59, 62], which presumably are gleaned from their nightly resting places [59, 62]. However, these studies were based on examination of faecal pellets usually gathered from below nursery colonies between May and September [60, 61, 63], and none assessed the diet of M. nattereri during winter.

Due to the seasonal and geographic variation in the abundance of prey [60] and changes in bats’ seasonal requirements, considerable differences in diet have been observed both among study sites and at different times of the year in rhinolophid (horseshoe) bats [64, 65]. In a study of the winter diet of the greater horseshoe bat (Rhinolophus ferrumequinum), Ransome [21] found that insects in the winter diet were broadly similar in type (but not importance) to those found in their summer diet. As during summer, there was variation in prey species taken between study sites [64]. However, prey species taken at sites (but not proportions) remained consistent across years.

Predation of Lepidoptera

The tympanic organs of many moths appear to have evolved mainly to detect echolocating bats [66]. In free flight, M. nattereri emits broadband calls that pass through the frequency range where moth hearing is at its most sensitive [67]. These factors probably render M. nattereri poorly equipped for capturing adult tympanate moths. This is largely supported by dietary studies of M. nattereri during the summer where Lepidoptera make up a relatively small component (5.4% [63] and 1.2% [60]) of the diet. However, in contrast with previous studies, our findings showed that at least at this particular study site, adult Lepidoptera make a relatively high proportion (>16%) of the diet of M. nattereri during winter. The temperature threshold for flight in many moth species is 6-10°C [15] (though the threshold may be lower in species that fly during winter). During our study periods, ambient temperatures ranged from −13.0–18.3°C. The ability of an insect to fly occurs abruptly at a critical temperature [3]. Therefore, moths grounded by a drop in temperature to below their flight threshold may become vulnerable to gleaning bats. Furthermore, many moth species flutter their wings when warming up, and the sound of wing-fluttering is exploited by gleaning bats [68]. Together, these features may explain our finding that adult Lepidoptera make up a larger part of the diet of M. nattereri during winter.

Adult lepidopteran prey remains were often encountered during the morphological analysis and the molecular study identified six lepidopteran species that would have been adult when eaten during winter (Table 3). These comprised two noctuid, two geometrid, one tortricid and one ypsolophid species (Table 2; 3). Both noctuid and geometrid moths are known to possess tympanic organs [42], and the geometrid winter moth Operophtera brumata initiates evasive flight in response to ultrasound, indicating bat avoidance behaviour [69]. However, tympanic organs have not been found in the Tortricidae [70] and there is no published work on the presence of tympanic organs in the Ypsolophidae.

Unexpectedly, we detected large quantities of Lepidopteran larvae in the winter diet of M. nattereri. Lepidopteran larvae have previously been recorded within the diets of bats [7173] including M. nattereri[59, 74]. The proportion of lepidopteran larvae recorded in the morphological analysis should be considered a conservative estimate due to a reduced likelihood of identifiable body parts of soft-bodied prey surviving passage through the bats’ digestive systems [75, 76]. Many Lepidoptera overwinter in their larval stages and are active on mild winter nights [44, 77], and some species are able to withstand temperatures below freezing [7880]. Once located, caterpillars should be easily captured by a gleaning bat and provide considerable nutrition [81]. Gleaning bats may eat larvae hanging from silk threads [73] when they would be an obvious target readily detectable by echolocation [73]. It seems unlikely, though, that larvae would be hanging from threads during winter nights. We have not identified where the M. nattereri in our study caught the larvae identified in the diet nor the foraging strategy used to gather them. However, the larvae we encountered in the diet of M. nattereri do feed on mild winter nights, and they are most commonly encountered in open grassy lowland habitats [44].

Globally, bats are known to provide a range of valuable ecosystem services [82], including reducing the numbers of insect pests [83, 84]. We found that Noctuidae was the most frequently encountered lepidopteran family represented mainly by the square spot rustic, X. xanthographa, which were identified in 80.6% of the pellets in the molecular analysis, and the large yellow underwing, Noctua pronuba, which were found in 63% of the pellets in the molecular analysis. N. pronuba larvae are often referred to as ‘cutworms’ and are considered serious agricultural and horticultural pests in Europe [85] and North America [86]. Although the total abundance of UK moths decreased by 28% between 1968 and 2007, the UK population of N. pronuba increased by 186% over the same period [87]. Further studies are needed to determine whether M. nattereri provides valuable ecosystem services by consuming N. pronuba and related agricultural pest species.

The consumption of other prey species

Some arachnid species are active during winter [88, 89] even at very low temperatures [90, 91]. We identified 14 arachnid species in the diet of M. nattereri. Of these, approximately half are web-building while the others are non-web-building or hunting spiders. The species identified are found on the leaves of herbs, shrubs and trees both within the canopy and close to the ground (M. Nyffeler, personal communication). The lycosid species (Pardosa sp) can usually be found on the ground or on low herbs in woodlands, grassland and agricultural fields (M. Nyffeler, personal communication). Arachnids have been recorded frequently in the diet of M. nattereri during summer [6062, 76]. They are presumably gleaned from the ground or other surfaces, or taken directly from webs [92] where they may be obvious targets for echolocating bats. Linyphiid spiders have also been observed ‘ballooning’ in favourable weather conditions during winter [93] so the capture of some species by aerial hawking cannot be ruled out. The possibility of some arachnid species appearing in the diet as a result of secondary predation must be considered [50]. However, analysis of the contents from individual faecal samples shows that although some arachnids were found within the same sample, many were not; also morphological results confirm prey remains from different arachnid families. If secondary predation was involved, it might be expected that the prey were digested completely by the predatory spider and therefore not visible under the microscope.

Isopod fragments were recorded in both the morphological and molecular analysis. Molecular analysis identified fragments as belonging to the common striped woodlouse (Philoscia muscorum). Woodlouse body fragments are not commonly encountered within bat faecal remains, leading some to conclude that bats do not eat woodlice [94]. Woodlouse body parts have, however, been recorded previously in the diets of long-eared bats (Plecotus auritus and Plecotus austriacus) [37, 95], M. nattereri[95] and the lesser horseshoe bat (Rhinolophus hipposideros) [22]. Woodlice are known to be active at night and remain active throughout winter [96], and they were commonly observed to be active on winter nights within our study area (P. Hope, unpublished observations).

Woodlice appeared in a higher proportion in the DNA results than in morphological findings, which may be caused by a number of factors. It is possible that faecal pellets may have been contaminated with woodlice DNA prior to collection. However, the fact that woodlouse body parts were identified in faecal pellets during the microscopic analysis, verifies that they are consumed during winter. The lower proportion of woodlouse remains identified in the microscopic analysis may be because certain body parts are discarded (e.g. the pereon) by bats due to difficulty in digestion or have low nutritional value or energy content. This is supported by the fact that we only found woodlouse leg parts, antennae and mouthparts during the microscopic analysis of the faecal pellets.

In previous studies, Diptera was the most frequently recorded order in the diet of M. nattereri. We identified eight species from six families and over the whole winter period, Diptera made up 11.5% (morphological) and 16.4% (molecular) of the bats’ diet. There was a tenfold decrease in Diptera between December and January-February (2%) in the morphological findings, with a less pronounced decrease in the molecular findings (18% and 13% respectively). In the particularly mild maritime climate of Cornwall, UK, Williams et al. [22] recorded Diptera as the most commonly encountered prey item over winter in the diet of the lesser horseshoe bat (R. hipposideros), as is the case during summer [97]. Diptera were either less abundant during winter in our study area or M. nattereri may have selected alternative prey items that may be easier to capture or have greater nutritional value (e.g. Lepidoptera adults and larvae).

The majority of dipteran species identified within our study are known to be active during winter. Dung flies (Scathophaga stercoraria) were recorded as a dominant prey species in the winter diet of greater horseshoe bat (R. ferrumequinum) [21] and cluster flies (Pollenia rudis and Pollenia pediculate), which overwinter as adults, become active during spells of mild weather and are commonly found in bat roosts in the UK [98]. Also, dipteran species from the family Muscidae (of which three species were identified in our study) are active as adults during winter [99].

We identified the tachinid fly (Ramonda spathulata) among the prey. It is a known parasite of noctuid moth caterpillars [100]. Ramonda spathulata was identified molecularly in seven of the bat faecal samples and was always found in conjunction with noctuid larvae or an unidentified lepidopteran species, which suggests that it was consumed indirectly.

Sensory ecology and its influence on foraging

The ability to detect stationary or slow moving prey on or close to surfaces using echolocation would be highly advantageous for foraging in winter when low temperatures reduce insect activity. Our finding that M. nattereri eats a high proportion of non-volant prey concurs with studies conducted between the months of May and September. Shiel et al. [63] concluded that 68% of the prey items eaten by M. nattereri were non-volant. Differences in sensory ecology among gleaning bat species contributes to niche differentiation [18]. Some gleaning bat species listen for prey-generated sound cues to detect and locate prey [101, 102]. Mouse-eared bats (Myotis myotis) and (Myotis blythii) enter a ‘whispering mode’ of echolocation during the final stages of prey approach [103], while the brown long-eared bat (Plecotus auritus) ceases to use echolocation when attempting to capture prey in cluttered environments [61]. M. nattereri in contrast emits a feeding buzz [61] and uses very broadband echolocation pulses to locate immobile prey close to vegetation [104].

Some gleaning bat species appear unable to detect non-moving prey that remain on the surface of vegetation [105] while others can [106]. The detection of stationary prey amongst dense vegetation (e.g. leaf litter), however, would be more difficult due to the effects of increased acoustic clutter [107]. Some bats may be able to discern texture when using echolocation because spectral interference of echoes reflected from surfaces allows discrimination of targets [108]. It would be interesting to establish whether M. nattereri is able to detect sufficient textural detail with its echolocation calls to distinguish slow moving or stationary prey from background clutter. M. nattereri may use hairs fringing the outer margin of the uropatagium to detect prey at very close range [109], and bats of this species have been observed trawling through vegetation using their tail membrane to gather prey [61]. Myotis nattereri has also been recorded using quadrupedal movement on the ground to chase and capture prey [61, 74]. This behaviour may also be used to locate and capture non-volant or immobile prey during winter.

Conclusions

Until now, the winter diet of vespertilionid bats in locations with mild temperate climate has been largely overlooked. This study indicates that, at least in a location with a mild winter climate, M. nattereri forages throughout winter. Morphological data show that a high proportion of non-volant arthropods are taken as prey throughout the winter, while the molecular data provide a very detailed picture of the winter diet of M. nattereri, identifying 43 prey species from six taxonomic orders. Previous studies on the diet of M. nattereri in summer revealed high incidence of Diptera in the diet, and studies in a flight tent revealed that the species catches prey close to surfaces by using broadband echolocation calls. Unexpectedly, we detected large quantities of lepidopteran larvae, which must have been captured by gleaning. The high occurrence of other non-volant prey also suggests that a gleaning foraging strategy is employed during winter, which may be effective at low ambient temperatures when the abundance of flying insects may be reduced substantially.

Materials and methods

Study site

Our study site was in Greywell Tunnel, which is a 1,125 m-long canal tunnel on the Basingstoke Canal, Hampshire, Southern England, United Kingdom (UK) (51.266380, −0.988030 (WGS1984)). Greywell Tunnel has not been used for navigation since the 1930s when a roof fall blocked a central section [110]. Both east and west portals of the tunnel have been grilled to prevent unauthorised human access, though grilles are designed to allow bats to freely use the site. The tunnel was designated a Site of Special Scientific Interest (SSSI) in 1985 in recognition of its importance as a bat hibernation site [111].

Average, minimum and maximum daily temperature readings from December 1 to March 31 2009–10 and 2010–11 were measured at a weather station at a Royal Air Force (RAF) base (RAF Odiham) 2 km from Greywell Tunnel. Mean daily temperatures (with SD and range in parentheses) between December 1 and March 31 were: in 2009–10, 3.4 ±3.6°C (−8.6–16.0°C) and in 2010–11, 4.3 ±4.1°C (−13.0–18.3°C).

Sample collection and preparation

Bat faecal samples were collected from within the western end of Greywell Tunnel during the winters of 2009–10 and 2010–11. Over the collection period, the number of visible M. nattereri ranged from ten (February 2011) to 63 (January 2011). A radio telemetry study conducted during 2008–09 and 2009–10 [13] showed that there were changes in the hibernating population with new individuals coming to and leaving the site over the winter. Bats are also known to roost in cavities behind the tunnel’s brick lining [25], so counts of visual bats may underestimate the total number of bats present.

Heavy-gauge white polythene sheeting (3.9 m × 4.9 m) was laid down on a clay surface close to the roof fall at approximately 130 m from the western portal. Faecal samples were collected between December-March, to minimise disturbance to hibernating bats samples were collected at the end of each month. During collection, each faecal sample was placed in a sterile 1.5 ml tube. After sample collection, the sheeting was cleared of bat faeces to avoid any mixing of faecal material between survey months. Fewer samples were collected during January and February of both study years so the samples from these months were pooled to produce a mid-winter sample size sufficiently large for comparison with December and March samples. Samples were stored at room temperature.

DNA extraction

Fifteen droppings were selected randomly from each of the three survey periods for both study years (n = 90). DNA was extracted from each faecal pellet using QIAamp DNA Stool Mini Kits (Qiagen Ltd., Crawley, West Sussex, UK) following the methodology as set out by Zeale et al. [34], all extractions were conducted during a single, two week period. After DNA extraction, faecal pellet remains were placed in 2 ml tubes and frozen for later morphological analysis.

Assigning pellets to bat species

To identify bat species to pellets, faecal extracts were amplified using forward (R3.1-F 5′- TGA AAA ACC ATC GTT GTA TTT CAA CTA CAA-3′ [112]) and reverse (MVZ04 5′-GCA GCC CCT CAG AAT GAT ATT TGT CCT C -3′ [113]) primers amplifying a 405 bp fragment of the mitochondrial DNA (mtDNA) cytochrome b gene. A HotStar Taq Plus Master Mix Kit (Qiagen Ltd., Crawley, West Sussex, UK) was used. PCR mix for each sample (20 μl) was as follows: 1 μl of extracted faecal DNA, 1 μl each of forward and reverse primers (final concentration in reaction 0.1 μM), 2 μl of loading buffer (CoraLoad), 10 μl HotStarTaq Master Mix (HotStar Taq plus DNA polymerase, PCR buffer (with 1.5 mM MgCl2) and 200 μM each dNTP) and 5 μl H2O. PCRs were carried out on a PTC-200 Thermo Cycler (MJ Research, Reno, Nevada) with the following conditions: 95°C for 5 minutes followed by 34 cycles of 94°C for 30 seconds, 50°C for 30 seconds, 72°C for 1 minute, followed by a final extension at 72°C for 10 minutes. Six blank negative controls were included in the PCR. PCR products were visualised on a 1.5% agarose gel stained with Web Green DNA stain (Web Scientific). Sanger sequencing of PCR products was performed in one direction (using the primer MVZ04) and carried out by a commercial laboratory (LGC Genomics, Germany).

Sequences were compared with reference sequences deposited in GenBank using the Basic Local Alignment Tool (BLAST) to obtain ‘closest match’ identifications. Matches were classed as positively identified to bat species if maximum identity was at ≥98% and query coverage ≥70%.

Morphological analysis

Morphological analysis of faecal material was conducted using traditional microscopic methods [32, 63, 76]. Faecal material was re-suspended in a 90 mm diameter Petri dish; a 2x2 mm piece of graph paper was fixed to the bottom of the dish to aid judgement of fragment size. One drop of glycerol, then one drop of alcohol (70%) were added to help spread prey fragments evenly over the Petri dish. Fragments were examined using a binocular microscope (MX3, 20x and 40x magnifications, Brunel, Chippenham, Wiltshire, UK). All fragments considered suitable for identification were mounted on a microscope slide in glycerinated gelatine under an 18x18 mm cover slip.

Arthropod fragments were compared with examples in published keys or identification guides [76, 114] and attempts were made to identify fragments to insect order, or where possible family. The relative importance of different prey orders in the diet was quantified in two ways. First, we calculated the percentage occurrence (%O) of prey orders, which is the number of pellets that an order was found in / total number of pellets in samples x 100 [76]. Second, we calculated the percentage frequency (%F), which is defined as the number of occurrences of an order divided by total occurrences of all orders x 100 [76]. Whereas %O produces values that exceed 100% when all taxa are combined (as pellets can contain more than one taxon), %F gives a proportional representation of the relative importance of prey taxa rounded to 100%.

Molecular analysis

In order to identify insect-prey in the bat droppings, amplifications of ca. 157 bp prey-insect-COI fragments were performed on the faecal extracts following a modification of the methods described in Bohmann et al. [36]. The assay described here was customised for Illumina MiSeq sequencing as opposed to GS FLX sequencing used by Bohmann et al. [36]. Specifically, insect-generic COI mini-barcode primers (ZBJ-ArtF1c and ZBJ-ArtR2c [34] were 5′ nucleotide barcoded [115], yielding a set of 20 forward and 20 reverse primers, all varying by the 8 bp of the 5′ barcode (Additional file 1), producing a final amplicon size of 227 bp (157 bp insert, 70 bp primers). Subsequently, all PCRs performed on the DNA extracts were performed using different combinations of the 5′ barcoded forward and reverse primers, ensuring that all resulting positive amplicons from each PCR reaction were uniquely labelled. To aid downstream bioinformatics processing of the data, and in particular to enable discrimination of PCR and/or sequencing artefacts from true biological sequences, each extract was subjected to PCR twice, in each case incorporating a unique barcode combination to each reaction. PCR reactions were performed in 25 μl reactions using the Amplitaq Gold enzyme system (Roche, Basel, Switzerland). For every eight reactions, a PCR blank was included as negative control. Furthermore, extraction blanks were amplified. Each 25 ul reaction contained: 1 μl DNA, 1x PCR Gold buffer, 2.5 mM MgCl2 solution, 0.2 mM dNTPs, 0.2 μl AmpliTaq Gold, and 0.4 μM of each primer. Amplifications were performed on a 2720 Thermal Cycler (Applied Biosystems) with the following conditions: 95°C for 5 minutes, then 40 cycles of 95°C for 15 seconds, 52°C for 30 seconds and 72°C for 30 seconds, followed by a final extension at 72°C for 7 minutes. 5 μl of each of the PCR products were transferred to a tube under a flow hood before being visualised with GelRed Nucleic Acid Stain (Biotium) on 2% agarose gels against a 50 bp ladder. Negative controls were negative and were not included in further analyses.

All PCR products were pooled into a single pool before library building, at approximately equimolar ratios (as determined by gel band strength). Specifically, bands were visually discriminated into three categories: bright, semi-bright and faint, from which PCR products were added to the pool with 1 μl, 5 μl and 15 μl, respectively. The pools of PCR products were purified using the MinElute PCR Purification Kit (Qiagen) following the MinElute Handbook 03/2008 MinElute PCR Purification Kit Protocol (using a microcentrifuge) with the following modifications: in step 4, the sample was centrifuged at 6000 g. In step 6, 700 μl Buffer PE was added, and the column was centrifuged at 10,000 g. In step 7, the empty column was centrifuged at 13,000 g for 2 minutes. In step 9, after Buffer EB was added to the centre of the membrane, the column was incubated for 5 minutes at 37°C. The column was centrifuged for 1 minute at 13,000 g, and the tube was turned and the spin repeated before DNA was collected.

The purified PCR pool was size-selected (227 bp +/−10%) on a LabChip XT (Caliper), before being purified using MinElute columns following the above-mentioned protocol and eluted in 42.5 μl Qiagen Buffer EB. The concentration of the purified, size-selected PCR pool was measured using a Qubit Flourometer (Invitrogen). The pool was subsequently converted into an Illumina sequencing library, using the NEBNext DNA Library Prep Master Mix Set for 454 (#E6070L) although using blunt end Illumina adapters [116] in place of Roche/454 FLX adaptors. The NEBNext End Repair Module Protocol was followed with the following modifications: end repair was performed in a 50 ul reaction using 42.5 μl size-selected, purified PCR pool. A control library blank was constructed in which 42.5 μl H2O was added in place of DNA. They were incubated in a 2720 Thermal Cycler (Applied Biosystems) for 20 minutes at 12°C and 15 minutes at 37°C, before being purified on MinElute columns following the above-mentioned protocol. In the Quick Ligation Module procedure, 0.5 μM Illumina sequencing adaptors [116] were included in the 50 μl reaction, and the mix was incubated on a 2720 Thermal Cycler (Applied Biosystems) for 20 minutes at 20°C before being purified on MinElute column following the above-mentioned protocol and eluted in 42 μl Buffer EB. In the fill-in reaction procedure, only step 8 was performed in which 5 μl Adapter Fill-in Reaction Buffer and 3 μl Bst DNA polymerase was added to the 42 μl DNA eluted above. This was incubated in a 2720 Thermal Cycler (Applied Biosystems) for 20 minutes at 65°C, followed by 20 minutes at 80°C to inactivate the enzyme.

The library was subsequently subjected to index PCR in 50 μl reactions using the Amplitaq Gold enzyme system (Roche, Basel, Switzerland): the library was amplified in five separate reactions each containing 5 μl library, alongside a reaction containing 5 μl of the library blank and one PCR blank in which 5 μl H2O was added. Furthermore, each reaction contained 1x PCR Gold Buffer, 2.5 mM MgCl2 solution, 0.25 mM dNTP’s, 5x purified BSA (BioLabs), 0.4 μl AmpliTaq Gold, 29.6 H2O, 0.2 μM forward index primer paired end (InPe 1.0) and 0.2 μM reverse index primer paired end. The five reactions containing the library were amplified with the same reverse index, while the library blank and the PCR blank were amplified with different reverse index primers. The index PCR was run on a 2720 Thermal Cycler (Applied Biosystems) at the following conditions: 95°C for 5 minutes, then 15 cycles of 95°C for 30 seconds, 60°C for 30 seconds and 72°C for 30 seconds, followed by a final extension at 72°C for 7 minutes. The five 50 μl indexed libraries were combined to a total of 250 μl before the library, library blank and PCR blank were purified on Qiagen QIAquick columns following the above-mentioned MinElute purification protocol and eluted in 30 μl Buffer EB. The library, the library blank and the PCR blank were visualised on a 2100 Bioanalyzer (Agilent Technologies). As no target fragments were detectable in either the library blank or the index PCR blank, only the library was sequenced. Sequencing was performed on the Illumina MiSeq platform where tagged amplicons were sequenced 150 bp paired-end, taking up ca. 15% of a run.

Initial bioinformatics analysis was conducted using the Quantitative Insights Into Microbial Ecology software package (QIIME) [117]. First, the quality of the paired-end reads was assessed with FastQC v 0.10.0 [118] and then low quality bases were trimmed at the 3′ end (q<20) with PRINSEQ v0.15 [119]. Although q<30 is commonly used in studies with modern DNA (good quality material), the DNA in our study showed more damage so we used q<20 which is normally used for ancient DNA which is more degraded. Afterwards, using customized Perl scripts (program available upon request) the corresponding paired reads (363,844) were merged if the overlap was 100% identical producing 140,009 sequences. In order to discriminate PCR and sequencing artefacts from true biological variation, only sequences present with 100% identity in both PCRs from the same extract were kept; this step yielded 114,345 sequences.

Also, the sequences were collapsed so that each sequence was only present once and reads shorter than 157 bp were filtered out. This gave a final filtered total of 1130 unique sequences that were used for the rest of the pipeline.

COI barcode sequences from different taxonomic levels (Arachnida, Chilopoda, Insecta, and Isopoda) were downloaded from the BOLD v3 [120] Public Data Portal to create a database and its corresponding taxonomy map to be used for QIIME v1.6.0 [117].

Species identification was made with QIIME with the uclust_ref OTU (Operational Taxonomic Unit) picking method (uclust v1.2.22q), using 98.5% identity (two mutations), with max_accepts and max_rejects set to 0 to force a search of the entire database with every query, thus guaranteeing that the best hit will be found if one exists. The biom table obtained from the OTU picking was summarised through plots with QIIME. 267 sequences failed to be assigned to a cluster (11%), and 2162 (89%) were assigned an OTU.

After the sequence data were obtained, we made the taxonomic assignments. OTUs were matched to most likely arthropod species using the identification tool within the BOLD database, using four confidence levels based on sequence similarity largely following Razgour et al. [37], whereby confidence level 1 = solid match to one species or genus (>98.5% sequence similarity); level = 2 match to more than one species (98.5%), only one of which was a UK species; level 3 = matched to two UK species of the same family (>98.5%) and level 4 = match to several species of different genera, or to reference sequences only identified to family (>98%). See Zeale et al. [34] for justification of the use of these levels of sequence similarity to taxonomic levels. As Noctua pronuba and Noctua janthe only differ by one base pair yielding a 99.4% similarity, these were clustered together at level 3. Since Noctua pronuba is of special interest due to its status as a serious agricultural and horticultural pest [85, 86], these sequences were separated at 100% similarity to identify the extent to which M. nattereri had preyed on Noctua pronuba at this site.

Field guides [44, 45, 77] and information from local records [52] were used to ascertain the phenology and life cycle stage of lepidopteran species identified genetically within the diet (Table 3).

References

  1. Aschoff J: Thermal conductance in mammals and birds: its dependence on body size and circadian phase. Comp Biochem Physiol A Physiol. 1981, 69a: 611-619.

    Google Scholar 

  2. Geiser F: Metabolic rate and body temperature reduction during hibernation and daily topor. Annu Rev Physiol. 2004, 66: 239-274. 10.1146/annurev.physiol.66.032102.115105.

    PubMed  CAS  Google Scholar 

  3. Taylor LR: Analysis of the effects of temperature on insects in flight. J Anim Ecol. 1963, 32: 99-117. 10.2307/2520.

    Google Scholar 

  4. Bale JS, Masters GJ, Hodkinson ID, Awmack C, Bezemer TM, Brown VK, Butterfield J, Buse A, Coulson JC, Farrar J, Good JEG, Harrington R, Hartley S, Jones TH, Lindroth RL, Press MC, Symrnioudis I, Watt AD, Whittaker JB: Herbivory in global climate change research: direct effects of rising temperatures on insect herbivores. Glob Change Biol. 2002, 8: 1-16. 10.1046/j.1365-2486.2002.00451.x.

    Google Scholar 

  5. Daan S: Activity during natural hibernation in three species of vespertilionid bats. Neth J Zool. 1973, 23: 1-71.

    Google Scholar 

  6. Thomas DW: The physiological ecology of hibernation in vespertilionid bats. Symp Zool Soc Lond. 1995, 67: 233-244.

    Google Scholar 

  7. Jonasson KA, Willis CK: Hibernation energetics of free-ranging little brown bats. J Exp Biol. 2012, 215: 2141-2149. 10.1242/jeb.066514.

    PubMed  Google Scholar 

  8. Czenze ZJ, Park AD, Willis CKR: Staying cold through dinner: cold-climate bats re-warm with conspecifics but not sunset during hibernation. J Comp Physiol B. 2013, 183: 859-866. 10.1007/s00360-013-0753-4.

    PubMed  Google Scholar 

  9. Ransome RD: Effect of ambient temperature on arousal frequency of hibernating greater horseshoe bat, Rhinolophus ferrumequinum, in relation to site and hibernation state. J Zool. 1971, 164: 353-371.

    Google Scholar 

  10. Park KJ, Jones G, Ransome RD: Torpor, arousal and activity of hibernating greater horseshoe bats (Rhinolophus ferrumequinum). Funct Ecol. 2000, 14: 580-588. 10.1046/j.1365-2435.2000.t01-1-00460.x.

    Google Scholar 

  11. Turbill C, Geiser F: Hibernation by tree roosting bats. J Comp Physiol. 2008, 178: 597-605. 10.1007/s00360-007-0249-1.

    Google Scholar 

  12. Liu J-N, Karasov WH: Hibernation in warm hibernacula by free-ranging Formosan leaf-nosed bats, Hipposideros terasensis, in subtropical Taiwan. J Comp Physiol. 2011, 181: 125-135. 10.1007/s00360-010-0509-3.

    Google Scholar 

  13. Hope PR, Jones G: Warming up for dinner: torpor and arousal in hibernating Natterer’s bats (Myotis nattereri) studied by radio telemetry. J Comp Physiol. 2012, 182: 569-578. 10.1007/s00360-011-0631-x.

    Google Scholar 

  14. Thomas DW, Dorais M, Bergeron J: Winter energy budgets and cost of arousals for hibernating little brown bats Myotis lucifugus. J Mammal. 1990, 71: 475-479. 10.2307/1381967.

    Google Scholar 

  15. Jones G, Duvergé PL, Ransome RD: Conservation biology of an endangered species: field studies of greater horseshoe bats. Symp Zool Soc Lond. 1995, 67: 309-324.

    Google Scholar 

  16. Mellanby K: Low temperature and insect activity. Proc R Soc Lond B Biol Sci. 1939, 127: 473-487. 10.1098/rspb.1939.0035.

    Google Scholar 

  17. Waters DA, Jones G: Echolocation call structure and intensity in five species of insectivorous bats. J Exp Biol. 1995, 198: 475-489.

    PubMed  CAS  Google Scholar 

  18. Siemers BM, Swift SM: Differences in sensory ecology contribute to resource partitioning in the bats Myotis bechsteinii and Myotis nattereri (Chiroptera: Vespertilionidae). Behav Ecol Sociobiol. 2006, 59: 373-380. 10.1007/s00265-005-0060-5.

    Google Scholar 

  19. Avery MI: Winter activity of pipistrelle bats. J Anim Ecol. 1985, 54: 721-738. 10.2307/4374.

    Google Scholar 

  20. Avery MI: The winter activity of noctule bats (Nyctalus noctula). J Zool. 1986, 209: 296-299. 10.1111/j.1469-7998.1986.tb03590.x.

    Google Scholar 

  21. Ransome RD: Winter feeding studies on greater horseshoe bats. English Nature Research Report 449. 2002, English Nature: Peterborough

    Google Scholar 

  22. Williams C, Salter L, Jones G: The winter diet of the lesser horseshoe bat (Rhinolophus hipposideros) in Britain and Ireland. Ital J Mammal. 2011, 22: 159-166.

    Google Scholar 

  23. Miková E, Varcholová K, Boldogh S, Uhrin M: Winter diet analysis in Rhinolophus euryale (chiroptera). Cent Eur J Biol. 2013, 8: 848-853. 10.2478/s11535-013-0199-9.

    Google Scholar 

  24. Park KJ, Jones G, Ransome RD: Winter activity of a population of greater horseshoe bats (Rhinolophus ferrumequinum). J Zool. 1999, 248: 419-427.

    Google Scholar 

  25. Hope PR, Jones G: An entrained circadian cycle of peak activity in a population of hibernating bats. J Mammal. 2013, 94: 497-505. 10.1644/12-MAMM-A-095.1.

    Google Scholar 

  26. Jones G, Rydell J: Foraging strategy and predation risk as factors influencing emergence time in echolocating bats. Philos T Roy Soc B. 1994, 346: 445-455. 10.1098/rstb.1994.0161.

    Google Scholar 

  27. Williams CB: The times of activity of certain nocturnal insects, chiefly Lepidoptera, as indicated by a light-trap. T Roy Entomol Soc. 1935, 83: 523-555.

    Google Scholar 

  28. Thomas DW: Lack of evidence for a biological alarm clock in bats (Myotis spp.) hibernating under natural conditions. Can J Zool. 1993, 71: 1-3. 10.1139/z93-001.

    Google Scholar 

  29. Whitaker JO, Rissler LJ: Winter activity of bats at a mine entrance in Vermillion County, Indiana. Am Midl Nat. 1992, 127: 52-59. 10.2307/2426321.

    Google Scholar 

  30. Whitaker JO, Rissler LJ: Do bats feed in winter?. Am Midl Nat. 1993, 129: 200-203. 10.2307/2426448.

    Google Scholar 

  31. Lausen CL, Barclay RMR: Winter bat activity in the Canadian prairies. Can J Zool. 2006, 84: 1079-1086. 10.1139/z06-093.

    Google Scholar 

  32. Whitaker JO, McCracken GF, Siemers BM: Food habits analysis of insectivorous bats. Ecological and behavioral methods for the study of bats 2nd edition. Edited by: Kunz TH, Parsons S. 2009, Baltimore: Johns Hopkins University Press, 567-592.

    Google Scholar 

  33. Clare EL, Fraser EE, Braid HE, Fenton MB, Hebert PDN: Species on the menu of a generalist predator, the eastern red bat (Lasiurus borealis): using a molecular approach to detect arthropod prey. Mol Ecol. 2009, 18: 2532-2542. 10.1111/j.1365-294X.2009.04184.x.

    PubMed  Google Scholar 

  34. Zeale MRK, Butlin RK, Barker GLA, Lees DC, Jones G: Taxon-specific PCR for DNA barcoding arthropod prey in bat species. Mol Ecol Resour. 2011, 11: 236-244. 10.1111/j.1755-0998.2010.02920.x.

    PubMed  CAS  Google Scholar 

  35. Alberdi A, Garin I, Aizpurua O, Aihartza J: The foraging ecology of the mountain long-eared bat Plecotus macrobullaris revealed with DNA mini-barcodes. PLoS ONE. 2012, 7: e35692-10.1371/journal.pone.0035692.

    PubMed  CAS  PubMed Central  Google Scholar 

  36. Bohmann K, Monadjem A, Noer CL, Rasmussen M, Zeale MRK, Clare E, Jones G, Willerslev E, Gilbert MTP: Molecular diet analysis of two African free-tailed bats (Molossidae) using high throughput sequencing. PLoS ONE. 2011, 6: e21441-10.1371/journal.pone.0021441.

    PubMed  CAS  PubMed Central  Google Scholar 

  37. Razgour O, Clare EL, Zeale MRK, Hanmer J, Schnell IB, Rasmussen M, Gilbert MTP, Jones G: High-throughput sequencing offers insight into mechanisms of resource partitioning in cryptic bat species. Ecol Evol. 2011, 1: 556-570. 10.1002/ece3.49.

    PubMed  PubMed Central  Google Scholar 

  38. Clare EL, Symondson WOC, Fenton MB: An inordinate fondness for beetles? Variation in seasonal dietary preferences of night roosting big brown bats (Eptesicus fuscus). Mol Ecol. 2013, doi:10.1111/mec.12519

    Google Scholar 

  39. Krüger F, Clare EL, Symondson WOC, Keišs O, Pētersons G: Diet of the insectivorous bat Pipistrellus nathusii during autumn migration and summer residence. Mol Ecol. 2013, doi:10.1111/mec.12547

    Google Scholar 

  40. Vesterinen EJ, Lilley T, Laine VN, Wahlberg N: Next Generation Sequencing of fecal DNA reveals the dietary diversity of the widespread insectivorous predator Daubenton’s Bat (Myotis daubentonii) in southwest Finland. PLoS ONE. 2013, 8 (11): e82168-10.1371/journal.pone.0082168.

    PubMed  PubMed Central  Google Scholar 

  41. Puechmaille SJ, Allegrini B, Boston ESM, Dubourg-Savage M-J, Evin A, Knochel A, Le Bris Y, Lecoq V, Lemaire M, Rist D, Teeling EC: Genetic analysis reveals further cryptic lineages within the Myotis nattereri species complex. Mamm Biol. 2012, 77: 224-228.

    Google Scholar 

  42. Roeder KD: Responses of the less sensitive acoustic sense cells in the tympanic organs of some noctuid and geometrid moths. J Insect Physiol. 1974, 20: 55-66. 10.1016/0022-1910(74)90123-1.

    PubMed  CAS  Google Scholar 

  43. Robert D, Edgecomb RS, Read MP, Hoy RR: Tympanal hearing in tachinid flies (Diptera, Tachinidae, Ormiini): the comparative morphology of an innovation. Cell Tissue Res. 1996, 284: 435-448. 10.1007/s004410050604.

    PubMed  CAS  Google Scholar 

  44. Waring P, Townsend M: Field guide to the moths of Great Britain and Ireland. 2009, Gillingham: British Wildlife Publishing Ltd

    Google Scholar 

  45. Sterling P, Parsons M: Field guide to the micro moths of Great Britain and Ireland. 2012, Gillingham: British Wildlife Publishing Ltd

    Google Scholar 

  46. Deagle BE, Kirkwood R, Jarman SN: Analysis of Australian fur seal diet by pyrosequencing prey DNA in faeces. Mol Ecol. 2009, 18: 2022-2038. 10.1111/j.1365-294X.2009.04158.x.

    PubMed  CAS  Google Scholar 

  47. Pegard A, Miquel C, Valentini A, Coissac E, Bouvier F, François D, Taberlet P, Engel E, Pompanon F: Universal DNA-based methods for assessing the diet of grazing livestock and wildlife from feces. J Agric Food Chem. 2009, 57: 5700-5706. 10.1021/jf803680c.

    PubMed  CAS  Google Scholar 

  48. Rasmussen M, Cummings LS, Gilbert MTP, Bryant V, Smith C: Response to comment by Goldberg et al. on “DNA from pre-clovis human coprolites in Oregon, North America.”. Science. 2009, 325: 148-148. 10.1126/science.1167531.

    CAS  Google Scholar 

  49. Soininen EM, Valentini A, Coissac E, Miquel C, Gielly L, Brockmann C, Brystring AK, Sønstebø JH, Ims RA, Yoccoz NG, Taberlet P: Analysing diet of small herbivores: the efficiency of DNA barcoding coupled with high-throughput pyrosequencing for deciphering the composition of complex plant mixtures. Front Zool. 2009, 6: 16-10.1186/1742-9994-6-16.

    PubMed  PubMed Central  Google Scholar 

  50. Pompanon F, Deagle BE, Symondson WOC, Brown DS, Jarman SN, Taberlet P: Who is eating what: diet assessment using next generation sequencing. Mol Ecol. 2012, 21: 1931-1950. 10.1111/j.1365-294X.2011.05403.x.

    PubMed  CAS  Google Scholar 

  51. Sheppard SK, Bell J, Sunderland KD, Fenlon J, Skervin D, Symondson WOC: Detection of secondary predation by PCR analysis of the gut contents of invertebrate generalist predators. Mol Ecol. 2005, 14: 4461-4468. 10.1111/j.1365-294X.2005.02742.x.

    PubMed  CAS  Google Scholar 

  52. Hantsmoths. Mike Wall [http//www.hantsmoths.org.uk]

  53. Anderson-Carpenter LL, McLachlan JS, Jackson ST, Kuch M, Lumibao CY, Poinar HN: Ancient DNA from lake sediments: bridging the gap between paleoecology and genetics. BMC Evol Biol. 2011, 11: 30-10.1186/1471-2148-11-30.

    PubMed  CAS  PubMed Central  Google Scholar 

  54. Emrich MA, Clare EL, Symondson WOC, Koenig SE, Fenton MB: Resource partitioning by insectivorous bats in Jamaica. Mol Ecol. 2014, doi:10.1111/mec.12504

    Google Scholar 

  55. Boessenkool S, Epp LS, Haile J, Bellemain E, Edwards M, Coissac E, Willerslev E, Brochmann C: Blocking human contaminant DNA during PCR allows amplification of rare mammal species from sedimentary ancient DNA. Mol Ecol. 2012, 21: 1806-1815. 10.1111/j.1365-294X.2011.05306.x.

    PubMed  CAS  Google Scholar 

  56. Coghlan ML, Haile J, Houston J, Murray DC, White NE, Moolhuijzen P, Bellgard MI, Bunce M: Deep sequencing of plant and animal DNA contained within traditional Chinese medicines reveals legality issues and health safety concerns. PLoS Genet. 2012, 8: e1002657-10.1371/journal.pgen.1002657.

    PubMed  CAS  PubMed Central  Google Scholar 

  57. Murray DC, Pearson SG, Fullagar R, Chase BM, Houston J, Atchison J, White NE, Bellgard MI, Clarke E, Macphail M, Gilbert MTP, Haile J, Bunce M: High-throughput sequencing of ancient plant and mammal DNA preserved in herbivore middens. Quat Sci Rev. 2012, 58 (c): 135-145.

    Google Scholar 

  58. Boston ESM, Puechmaille SJ, Scott DD, Buckley DJ, Lundy MG, Montgomery IW, Prodöhl PA, Teeling EC: Empirical assessment of non-invasive population genetics in bats: comparison of DNA quality from faecal and tissue samples. Acta Chiropterol. 2012, 14: 45-52. 10.3161/150811012X654259.

    Google Scholar 

  59. Vaughan N: The diets of British bats (Chiroptera). Mammal Rev. 1997, 27: 77-94. 10.1111/j.1365-2907.1997.tb00373.x.

    Google Scholar 

  60. Swift SM: Roosting and foraging behaviour of Natterer’s bats (Myotis nattereri, Kuhl, 1818) close to the northern border of their distribution. J Zool. 1997, 242: 375-384.

    Google Scholar 

  61. Swift SM, Racey PA: Gleaning as a foraging strategy in Natterer’s bat Myotis nattereri. Behav Ecol Sociobiol. 2002, 52: 408-416. 10.1007/s00265-002-0531-x.

    Google Scholar 

  62. Gregor F, Bauerovà Z: The role of Diptera in the diet of Natterer’s bat, Myotis nattereri. Folia Zool. 1987, 36: 13-19.

    Google Scholar 

  63. Shiel CB, McAney CM, Fairley JS: Analysis of the diet of Natterer’s bat Myotis nattereri and the common long-eared bat Plecotus auritus in the west of Ireland. J Zool. 1991, 233: 299-305.

    Google Scholar 

  64. Jones G: Prey selection by the greater horseshoe bat (Rhinolophus ferrumequinum): optimal foraging by echolocation?. J Anim Ecol. 1990, 59: 587-602. 10.2307/4882.

    Google Scholar 

  65. Knight T: The use of landscape features and habitats by the lesser horseshoe bat (Rhinolophus hipposideros). 2006, PhD thesis. University of Bristol, Biological Sciences

    Google Scholar 

  66. Rydell J, Jones G, Waters D: Bats and hearing moths: who are the winners?. Oikos. 1995, 73: 419-424. 10.2307/3545970.

    Google Scholar 

  67. Fenton BM, Fullard JH: The influence of moth hearing on bat echolocation strategies. J Comp Physiol. 1979, 132: 77-86. 10.1007/BF00617734.

    Google Scholar 

  68. Werner TK: Responses of nonflying moths to ultrasound: the threat of gleaning bats. Can J Zool. 1981, 59: 525-529. 10.1139/z81-076.

    Google Scholar 

  69. Svensson GE, Rydell J, Brown R: Bat predation and flight timing of winter moths, Epirrita and Operophtera species (Lepidoptera, Geometridae). Oikos. 1999, 84: 193-198. 10.2307/3546713.

    Google Scholar 

  70. Scoble MJ: The Lepidoptera. Form, function and diversity. 1992, Oxford: Natural History Museum and Oxford University Press

    Google Scholar 

  71. Jones G, Morton M, Hughes PM, Budden RM: Echolocation, flight morphology and foraging strategies of some West African hipposiderid bats. J Zool. 1993, 230: 385-400. 10.1111/j.1469-7998.1993.tb02691.x.

    Google Scholar 

  72. Shiel CB, Duvergé PL, Smiddy P, Fairley JS: Analysis of the diet of Leisler’s bat (Nyctalus leisleri) in Ireland with some comparative analysis from England and Germany. J Zool. 1998, 246: 417-425. 10.1111/j.1469-7998.1998.tb00173.x.

    Google Scholar 

  73. Wilson JM, Barclay MR: Consumption of caterpillars by bats during an outbreak of western spruce budworm. Am Midl Nat. 2006, 155: 244-249. 10.1674/0003-0031(2006)155[0244:COCBBD]2.0.CO;2.

    Google Scholar 

  74. Arlettaz R: Foraging behaviour of the gleaning bat Myotis nattereri (Chrioptera, Vespertilionidae) in the Swiss Alps. Mammalia. 1996, 60: 181-186.

    Google Scholar 

  75. Kunz TH: An evaluation of fecal analysis for determining food habits of insectivorous bats. Can J Zool. 1983, 61: 1317-1321. 10.1139/z83-177.

    Google Scholar 

  76. Shiel CB, McAney CM, Sullivan C, Fairley JS: Identification of arthropod fragments in bat droppings. No.17. 1997, London: Mammal Society

    Google Scholar 

  77. Carter DJ, Hargreaves B: Collins field guide: caterpillars of Britain and Europe. 1986, London: HarperCollins

    Google Scholar 

  78. Layne JR, Blakeley DL: Effect of freeze temperature on ice formation and long-term survival of the woolly bear caterpillar (Pyrrharctia isabella). J Insect Physiol. 2002, 48: 1133-1137. 10.1016/S0022-1910(02)00206-8.

    PubMed  CAS  Google Scholar 

  79. Sinclair BJ, Chown SL: Rapid responses to high temperature and desiccation but not to low temperature in the freeze tolerant sub-Antarctic caterpillar Pringleophaga marioni (Lepidoptera, Tineidae). J Insect Physiol. 2003, 49: 45-52. 10.1016/S0022-1910(02)00225-1.

    PubMed  CAS  Google Scholar 

  80. Layne JR: Freeze tolerance and cryoprotection in caterpillars of the giant leopard moth (Ecpantheria scribonia Lepidoptera: Arctiidae). J Therm Biol. 2005, 30: 267-271. 10.1016/j.jtherbio.2004.12.006.

    CAS  Google Scholar 

  81. Ashiru MO: The food value of larvae of Anaphe venata Butler (Lepidoptera: Notodontidae). Ecol Food Nutr. 1989, 22: 313-320. 10.1080/03670244.1989.9991080.

    Google Scholar 

  82. Kunz TH, Braun de Torrez E, Bauer D, Lobova T, Fleming TH: Ecosystem services provided by bats. Ann NY Acad Sci. 2011, 1233: 1-38. 10.1111/j.1749-6632.2011.06169.x.

    Google Scholar 

  83. Boyles JG, Cryan PM, McCracken GF, Kunz TH: Economic importance of bats in agriculture. Science. 2011, 332: 41-42. 10.1126/science.1201366.

    PubMed  Google Scholar 

  84. McCracken GF, Westbrook JK, Brown VA, Eldridge M, Federico P, Kunz TH: Bats track and exploit changes in insect pest populations. PLoS ONE. 2012, 7: e43839-10.1371/journal.pone.0043839.

    PubMed  CAS  PubMed Central  Google Scholar 

  85. Difonzo C, Russell H: Noctua pronuba (Lepidoptera: Noctuidae): an outbreak in emails. J Integr Pest Manag. 2010, 1: B1-B6.

    Google Scholar 

  86. Brou V: A Gulf Coast record of the European cutworm, Noctua pronuba (L.). News Lepid Soc. 1997, 19: 3-

    Google Scholar 

  87. Fox R, Parsons MS, Chapman JW, Woiwod IP, Warren MS, Brooks DR: The state of Britain’s larger moths 2013. 2013, Butterfly Conservation and Rothamsted Research: Wareham

    Google Scholar 

  88. Vanin S, Turchetto M: Winter activity of spiders and pseudoscorpions in South-Eastern alps (Italy). Ital J Zool. 2007, 74: 31-38. 10.1080/11250000601017233.

    Google Scholar 

  89. Aitchison CW: The phenology of winter-active spiders. J Arachnol. 1984, 12: 249-271.

    Google Scholar 

  90. Husby JA, Zachariassen KE: Antifreeze agents in the body fluid of winter active insects and spiders. Experientia. 1980, 36: 963-964. 10.1007/BF01953821.

    CAS  Google Scholar 

  91. Aitchison CW: Low temperature feeding by winter-active spiders. J Arachnol. 1984, 12: 297-305.

    Google Scholar 

  92. Kervyn T, Godin M-C, Jocqué R, Grootaert P, Libois R: Web-building spiders and blood-feeding flies as prey of the notch-eared bat (Myotis emarginatus). Belg J Zool. 2012, 142: 59-67.

    Google Scholar 

  93. Foelix R: Biology of spiders. 2010, New York: Oxford University Press

    Google Scholar 

  94. Siemers BM, Güttinger R: Prey conspicuousness can explain apparent prey selectivity. Curr Biol. 2006, 16: R157-R159. 10.1016/j.cub.2006.02.056.

    PubMed  CAS  Google Scholar 

  95. Gregory S: What eats woodlice (continued!). Myriapod and Isopod Group Newsletter 2004. Edited by: Lee P. 2004, 8-1–2

    Google Scholar 

  96. Tuf IH, Jeřábková E: Diurnal epigeic activity of terrestrial isopods (Isopoda: Oniscidea). Proceedings of the International Symposium on Terrestrial Isopod Biology. 2008, Aachen: Shaker Verlag, 167-172.

    Google Scholar 

  97. McAney CM, Fairley JS: Analysis of the diet of the lesser horseshoe bat Rhinolophus hipposideros, in the west of Ireland. J Zool. 1989, 217: 491-498. 10.1111/j.1469-7998.1989.tb02504.x.

    Google Scholar 

  98. Mitchell-Jones AJ: Timber treatment, pest control and building work. Bat workers’ manual. Edited by: Mitchell-Jones AJ, McLeish AP. 2004, Peterborough: Joint Nature Conservation Committee, 95-110. 3

    Google Scholar 

  99. Hwang C, Turner BD: Spatial and temporal variability of necrophagous Diptera from urban to rural areas. Med Vet Entomol. 2005, 19: 379-391. 10.1111/j.1365-2915.2005.00583.x.

    PubMed  CAS  Google Scholar 

  100. Ford TH, Shaw MR: Host records of some Western Palaearctic Tachinidae (Diptera). Entomol Rec J Var. 1991, 103: 23-38.

    Google Scholar 

  101. Fiedler J: Prey catching with and without echolocation in the Indian false vampire (Megaderma lyra). Behav Ecol Sociobiol. 1979, 6: 155-160. 10.1007/BF00292562.

    Google Scholar 

  102. Anderson EM, Racey PA: Feeding behaviour of captive brown long-eared bats (Plecotus auritus). Anim Behav. 1991, 42: 489-493. 10.1016/S0003-3472(05)80048-X.

    Google Scholar 

  103. Russo D, Jones G, Arlettaz R: Echolocation and passive listening by foraging mouse-eared bats Myotis myotis and Myotis blythii. J Exp Biol. 2007, 210: 166-176. 10.1242/jeb.02644.

    PubMed  Google Scholar 

  104. Siemers BM, Schnitzler H-U: Natterer’s bat (Myotis nattereri Kuhl, 1818) hawks prey close to vegetation using echolocation signals of very broad bandwidth. Behav Ecol Sociobiol. 2000, 47: 400-412. 10.1007/s002650050683.

    Google Scholar 

  105. Marimuthu G: Stationary prey ensures life and moving prey ensures death during the hunting flight of gleaning bats. Curr Sci. 1997, 72: 928-931.

    Google Scholar 

  106. Geipel I, Jung K, Kalko EKV: Perception of silent and motionless prey on vegetation by echolocation in a gleaning bat Micronycteris microtis. P Roy Soc Lond B Bio. 2013, 280: 20122830-10.1098/rspb.2012.2830.

    Google Scholar 

  107. Jones G: Sensory biology: listening in the dark for echoes from silent and stationary prey. Curr Biol. 2013, 23: R249-R251. 10.1016/j.cub.2013.02.016.

    PubMed  CAS  Google Scholar 

  108. Schmidt S: Evidence for a spectral basis of texture perception in bat sonar. Nature. 1988, 331: 617-619. 10.1038/331617a0.

    PubMed  CAS  Google Scholar 

  109. Siemers BM: Finding prey by associative learning in gleaning bats: experiments with a Natterer’s bat Myotis nattereri. Acta Chiropterol. 2001, 3: 211-215.

    Google Scholar 

  110. Vine PAL: London’s lost route to Basingstoke. The story of the Basingstoke Canal. 1994, David & Charles: Newton Abbot

    Google Scholar 

  111. Stebbings RE: Greywell Tunnel: an internationally important haven for bats. 1993, English Nature: Peterborough

    Google Scholar 

  112. Puechmaille SJ, Gouilh MA, Piyapan P, Yokubol M, Mie KM, Bates PJ, Satasook C, Nwe T, Bu SSH, Mackie IJ, Petit EJ, Teeling EC: The evolution of sensory divergence in the context of limited gene flow in the bumblebee bat. Nat Commun. 2011, 2: 573-

    PubMed  PubMed Central  Google Scholar 

  113. Smith MF, Patton JL: Variation in mitochondrial cytochrome b sequences in natural populations of South American akodontine rodents (Muridae: Sigmodontinae). Mol Biol Evol. 1992, 8: 85-103.

    Google Scholar 

  114. Chinery M: Collins field guide to the insects of Britain and western Europe. 2000, London: HarperCollins

    Google Scholar 

  115. Binladen J, Gilbert MTP, Bollback JP, Panitz F, Bendixen C, Nielsen R, Willerslev E: The use of coded PCR primers enables high-throughput sequencing of multiple homolog amplification products by 454 parallel sequencing. PLoS ONE. 2011, 2: e197-1–9

    Google Scholar 

  116. Meyer M, Kircher M: Illumina sequencing library preparation for highly multiplexed target capture and sequencing. Cold Spring Harb Protoc. 2010, 6: pdb.prot5448. doi:10.1101/pdb.prot5448

    Google Scholar 

  117. Quantitative Insights into Microbial Ecology website. [http://qiime.org/]

  118. Babraham Bioinformatics website. [http://www.bioinformatics.babraham.ac.uk/projects/fastqc/]

  119. Schmieder R, Edwards R: Quality control and preprocessing of metagenomic datasets. Bioinformatics. 2011, 27: 863-864. 10.1093/bioinformatics/btr026.

    PubMed  CAS  PubMed Central  Google Scholar 

  120. Ratnasingham S, Hebert PDN: BOLD: The barcode of life data system (http://www.barcodinglife.org). Mol Ecol Resour. 2007, 7: 355-364.

    CAS  Google Scholar 

Download references

Acknowledgements

We would like to thank the following people for their help with fieldwork, Roy Champion, Colleen Hope and volunteers from the Hampshire Bat Group. Thanks to Hampshire County Council and the team on the Basingstoke Canal for assistance and providing access to Greywell Tunnel. We thank Matt Zeale and Beth Clare for assistance with some lab based procedures, also Carol Williams and Marian James for assistance in identification of some insect fragments extracted from faecal pellets. The research in this paper was carried out under license from Natural England. It was supported by funding from Hampshire County Council (thanks to Andrew Davidson), the Basingstoke Canal Authority and Natural England via a Higher Level Stewardship (HLS) grant, a NERC grant to use the NERC Biomolecular Analysis Facilities – Sheffield and Edinburgh (MGF 242 to GJ ) and Danish National Research Foundation grant DNRF94. We also thank the Danish National High-Throughput DNA Sequencing Centre for their invaluable advice and technical assistance with regards to generation of the MiSeq data. We thank three anonymous reviewers for their valuable comments on the manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Paul R Hope.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

PRH conceived the study, conducted the morphological analysis, DNA extraction, amplification of bat DNA and drafted the manuscript. KB conducted the amplification and second generation sequencing of arthropod DNA and drafted part of the manuscript. MTPG supervised the amplification and second generation sequencing of arthropod DNA and commented on the manuscript. MLZM conducted the initial bioinformatics analysis of the second generation sequencing results and drafted part of the manuscript. OR assisted with extraction and PCR of bat DNA and commented on the manuscript. GJ conceived and supervised the study and assisted in drafting the manuscript. All authors read and approved the final manuscript.

Paul R Hope, Kristine Bohmann contributed equally to this work.

Electronic supplementary material

12983_2013_324_MOESM1_ESM.xlsx

Additional file 1:Zbjart_short primers. xls. Primers. 20 forward and 20 reverse primers, Illumina sequencing. (XLSX 15 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hope, P.R., Bohmann, K., Gilbert, M.T.P. et al. Second generation sequencing and morphological faecal analysis reveal unexpected foraging behaviour by Myotis nattereri (Chiroptera, Vespertilionidae) in winter. Front Zool 11, 39 (2014). https://doi.org/10.1186/1742-9994-11-39

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1742-9994-11-39

Keywords