The feeding habits of insectivorous bats are of great interest to people, because they are considered to be important in the control of insect pests. Here we present a study showing the relationship of bat morphology to differences in prey selection by various bat species. We compared dietary data from 10,884 faecal pellets, bodyweight, cranial length, forearm length and echolocation calls from published peer-reviewed studies for 92 bat species. We demonstrated that insectivorous bats tend to prefer certain insect orders which we have grouped as soft bodied insects, hard bodied insects and Lepidoptera. Wing characteristics which we measured by bodyweight-forearm ratio showed the strongest relationship with hard insects followed by longest cranial length. The content of soft insects in bat diets was negatively related to bodyweight, forearm length and longest cranial length. Lepidoptera content was positively related to the echolocation frequency with the maximum intensity (FMAXE), bats with high FMAXE fed on more Lepidoptera than those with low frequencies. We propose that a combination of dietary analysis and morphological analysis is needed to make strong inference about prey preference rather than comparing the dietary analysis with the insect abundance at the location were the bat or faecal pellets were collected
Chiroptera; diet; cranial morphology; bodyweight-forearm ratio; Echolocation.
The feeding habits of insectivorous bats are of great interest to people, because they are considered to be important in the control of insects ranging from agricultural pests such as moths (Lepidoptera) [1,2] to disease vectors such as mosquitoes (Culicidae) [3,4]. A thorough understanding of feeding habits and prey preference is still lacking regardless of the several hundred dietary analyses that have been conducted and published over the last few decades. Prey preference is difficult to determine in the field due to the large distances that bats fly to their foraging habitats . Analysing faecal pellets or stomach content does give an understanding of what the species is eating, but there are no studies known in which diet has been compared to the actual food supply. Many studies compare the dietary analysis with the insect abundance at the location were the bat or faecal pellets were collected, but this is in most cases not the location were the bat is foraging .
A few studies have combined the data from dietary analysis of many species to look for patterns related to biological characteristics such as skull morphology, wing span, bodyweight and/or echolocation which can help us predict the feeding habits and prey preference [6-8]. Food hardness, for example, is related to the bite force [8,10,11], which depends on skull morphology [12,13] and relates to body size or weight . Wing morphology influences feeding habits, it having been shown that different wing types cause different flight styles affecting agility and thus hunting behaviour . Bogdanowicz and colleagues  showed that the peak frequency (frequency of maximum intensity) of echolocation is also strongly related to different prey types. When peak frequencies become higher the amount of Lepidoptera in the diet increases; however, when the frequencies get lower, the amount of Coleoptera increases. Bats that use higher call frequencies are capable of detecting much smaller insects in comparison to bats that use low frequencies following the general physics of acoustics [16,17]. The wavelength has to be shorter than the size of an insect for a bat to be able to detect the prey . The small differences in skull morphology, flight style and echolocation have made bats very successful in many parts of the world because they can exploit a wide range of niches .
Here we present a study of the relative importance of different morphological characteristics in determining the dietary content of different insect types. We first conducted a comprehensive analysis of the co-occurrence of food items in the diets of insectivorous bat species. The content of ten insect orders in the diets of 92 bat species were compared. These data were acquired from existing peer-reviewed studies. We expected to find two clusters of insect orders; one with hard bodied insects and one with softer insects [8,10,11]. Finally we used these co-occurring insect orders to assess how different morphological characteristics are affecting prey selection in bats. We expected that body size, echolocation, skull size and forearm length are important factors that determine the feeding behaviour [7-9]. Diets for single bat species can differ tremendously between seasons and locality [2,19]. Therefore, we expected to find large variation in our dataset that could not be explained by these morphological characteristics. In comparison to previous studies, this study does not focus on a narrow range of food items, but includes an analysis of the complete diet of bats.
We analysed dietary content, bodyweight, longest cranial length, forearm length and echolocation calls based on data from published peer-reviewed studies. Data regarding diets were only considered when food items were presented as volume percentage and were based on faecal pellets. We excluded studies based on stomach content to avoid biases caused by using two different methods. Insectivorous bats digest their food very rapidly and do not thoroughly chew their prey. Therefore, faecal pellets often contain many parts of prey items that can still be used to identify the prey to the taxonomic order and sometimes to higher taxonomic levels [6,20]. When dietary data for multiple species were available, we calculated a weighted mean in which the number of studied pellets was used as weight. The final dataset consisted of 31 food categories: 25 insect orders, spiders (Araneae), two insect families, fish, plants and unidentified. These categories were later reduced to the ten most common invertebrate orders. For echolocation data we only considered the frequency with the maximum energy (FMAXE). Data regarding bodyweight, longest cranial length and forearm length were also acquired from published data. We used forearm length and longest cranial length as variables representing wing and skull morphology, because these are most often documented in existing literature When multiple values for the same species were encountered, we calculated the mean.
Nonmetric multidimensional scaling
We used Non-metric Multidimensional Scaling (NMS) based on Euclidean distances as a tool for creating ordination plots. The goodness of fit of a NMS can be assessed through the stress-values. A stress-values of 0.0 to 0.1 is considered a good fit, a stress value of 0.1 to 0.2 is considered a moderate fit and a higher stress values are considered to be a poor fit . We used dietary content of different insects as our input variables to assess the similarity of bat species based on their diets. Only those insects orders were used that were most common. This included all orders that occurred in the diet of 20 or more bat species. First these food items were plotted for all 92 species. Then we added grouping polygons based on the standard error of the weighted average. These polygons displayed the dominant food item based on three categories: hard bodied insects, soft bodied insects and Lepidoptera. Hard bodied insects included Coleoptera, Hemiptera, Homoptera and Orthoptera, while soft bodied insects included Diptera, Trichoptera, Neuroptera and Araneae. A category was considered dominant when the percentage was highest for this group in comparison to the other two categories. We also added polygons that represented the four most common bat genera in our dataset: Pippistrellus (n=8), Myotis (n=24), Hipposideros (n=6) and Rhinolophus (n=12). We compared the differences of bodyweight, longest cranial length, forearm length and FMAXE between these four genera using an ANOVA procedure with a multiple comparison test based on the Tukey-Kramer test proposed by Herberich et al. (2010). This procedure allows data to be heteroscedastic and unbalanced . Where necessary, data were log-transformed to obtain normality. All analyses were conducted with RStudio version 0.96.331  built on R version 3.0.1 . The 'vegan' package was used for NMS , while the 'sandwich' package was used for the multi-comparison procedure .
Generalized Linear Models
Food habits were modelled using weighted Generalized Linear Model (GLM) assuming a Tweedie compound Poisson-gamma distribution and a log-link function . The Tweedie distribution does not assume a constant variance, but relates the variance to the mean following the power law [27,28]. We used this distribution, because the data were following a Poisson shaped distribution but are continuous rather than count data. The Tweedie compound Poisson-gamma distribution allowed us to avoid any transformations of the response variables, which made it easier to compare the different models. All variables were separately modelled because of collinearity among explanatory variables. We did not want to exclude any variable or use a factor or dummy variable because we were interested in the relationship of each separate variable with the diet. Collinearity among explanatory variables was visualized in scatterplots to which we added fitted lines and confidence intervals based on generalized additive models using local regression smoothers. The weights that were included in the models were based on the number of pellets investigated for each bat species. The total dataset was based on a number of 10,884 faecal pellets. Prior to the analysis, all variables were standardized to a mean of zero and a standard deviation of one . The 'tweedie' package  was used for fitting the generalized linear models.
The NMS ordination of all main insect orders clearly showed that certain food items in bat faeces often co-occur (Figure 1A). Coleoptera, Hemiptera and Homoptera (hard bodied) were often found in faeces together. Pipistrellus and Hipposideros species tended to feed more often on these harder insects (Figure 1B). Soft bodied prey e.g. Diptera, Araneae and Trichoptera, were more often found in the faeces of other bat species. Lepidoptera deviated strongly from all the other insect orders and were not commonly found with other food items. Myotis species fed on considerably softer insects. Rhinolophus species fed most on Lepidoptera in comparison to the other three genera, while Pipistrellus species fed the least on Lepidoptera. Pipistrellus species were significantly smaller than species from the other three genera (Table 1) which was also reflected in forearm length and the longest cranial length. Differences in bodyweight between the other genera were not significant. Rhinolophus and Hipposideros species had significant higher FMAXE compared to Pipistrellus and Myotis species. Rhinolophus species had generally a lower but non-significant bodyweight than Hipposideros species, however, this was not reflected in the longest cranial length. Bat species with low bodyweight-forearm ratios (e.g. Pipistrellus) generally fed on softer insects in comparison to species with high bodyweight-forearm ratio (e.g. Hipposideros, Figure 1C).
Figure 1: Ordination plots based on nonmetric multidimensional scaling. The stress value for all plots is 0.161. The size of ellipses displays the standard error of the weighted averages. A) Dominant food types (30% or more) is in the diets of bats; hard insects, soft insects and Lepidoptera. B) The main bat generafor which diets were compared. C) Differences in diets of bats with low bodyweight-forearm ratios (0.0 – 0.4) and high bodyweight-forearm ratios (greater than 0.4). Ara = Araneae, Col = Coleoptera, Di = Diptera, Hem = Hemiptera, Hom = Homoptera, Hym = Hymenoptera, Lep = Lepidoptera, Neu = Neuroptera, Ort = Orthoptera, Tri = Trichoptera.
The GLMs how that the content of hard insects in diets is related to the bodyweight, cranial length, forearm length, bodyweight-forearm ratio and call frequency. Bodyweight-forearm ratio showed the strongest relationship with hard insects followed by longest cranial length (Table 2). Larger bats with relatively large skulls where most likely to have the highest content of hard insects in their diet. The FMAXE was negatively related to the content of hard insects in the diet. Nevertheless, the model based on FMAXE scored only slightly better than the null model. The content of soft insects in bat diets was related to bodyweight, forearm length and longest cranial length but not to FMAXE. Cranial length showed the strongest negative relationship to soft insect content. Smaller bats with relatively small skulls were most likely to have high amounts of soft insects in their diet. Lepidoptera content was positively related to FMAXE, bats with high FMAXE fed on more Lepidoptera than those with low frequencies. Lepidoptera content was not related to any of the other variables.
Bodyweight showed a non-linear relationship with longest cranial length (F = 19.16, d.f. = 88.3, P < 0.001, Figure 2) and forearm length (F = 33.4, d.f. = 88.3, P < 0.001). Cranial length and forearm length were showing a strong linear relationship (F = 262.7, d.f. = 90, P < 0.001). Larger bats generally had lower peak frequencies although this non linear relationship was not significant (F = 2.64, d.f. = 89.5, P = 0.118). There was a weak non linear relationship between FMAXE and cranial length (F = 9.75, d.f. = 89.5, P = 0.012) and between FMAXE and forearm length (F =5.62, d.f. = 89.5, P = 0.041).
Figure 2: Scatterplots of all explanatory variables. The solid lines display the fit of a generalized additive model with the confidence interval (dashed lines). Local regression smoothers (LOESS) were used to fit the models. A) Bodyweight vs. longest cranial length with a LOESS interval of 80%. B) Bodyweight vs. forearm length with a LOESS interval of 80%. C) Longest cranial length vs. forearm length in which the solid line represent a linear model. D) Bodyweight vs. FMAXE with a LOESS interval of 100%. E) Longest cranial length vs. FMAXE with a LOESS interval of 100%. F) Forearm length vs. FMAXE with a LOESS interval of 100%.
We demonstrated that insectivorous bats tend to prefer certain invertebrate orders which we have grouped as soft bodied invertebrates, hard bodied insects and Lepidoptera. Although bats do specialize and prefer certain prey, they do not exclude other prey from their diets. Preference of certain prey items can be explained by certain characteristics of the bats. The most important characteristics are bodyweight-forearm ratio, cranial length and FMAXE. Bodyweight-forearm ratio was most important in predicting the content of hard insects, the longest cranial length was most important for soft insects and FMAXE was most important for Lepidoptera. The differentiation of these characteristics in bats is likely to have evolved from trophic specialization [31,32]; as is hypothesized for other species [33,34]. Where several species of bats co-occur, it is in the interest of the individual species to specialize on prey that are most abundant or are not available for competitors [34,35]. Differences in flight style, cranial morphology and echolocation could therefore have found their origin within resource partitioning .
Wing morphology is really important in shaping the flight style e.g. aerial-hawking, trawling, gleaning or perch hunting . Most insectivorous bats catch their prey in flight, therefore flight style and wing morphology are very important in prey selection [38-40]. Agility and manoeuvrability of flight are for a great part defined by the wing loading, bodyweight-forearm ratio is an alternative parameter that closely represents the wing loading [40,41]. Species with a lower bodyweight-forearm ratio can be considered to be more agile, because they have less mass in comparison to the size of the wing . A higher bodyweight-forearm ratio relates to harder and larger insects in diets. These insects are easier detected from larger distances and therefore agility is less important. Pipistrellus and Hipposideros species differed significantly in their bodyweight-forearm ratio, from which we can infer that Pipistrellus species are generally more agile allowing them to catch smaller prey such as Diptera. This is reflected in the diets of both genera; where Hipposideros species mainly feed on harder insects and Lepidoptera, the diet of Pipistrellus species often consists of softer species.
The large differences in skull morphology among bat species is an important factor in predicting and understanding the diets of bats [43,44]. Harder insects are more often consumed by bats that have robuster jaws and thus have a stronger bite force [11,12,45]. Bite force is also related to the thickness of the teeth enamels. Therefore bats with thicker enamels often feed on harder insects . Species that lack these traits are therefore more likely to select softer prey. We used the longest cranial length to represent bite force. Our results show that cranial size and thus bite force is a better predictor than bodyweight for the diets containing soft and hard insects. This indicates that not just the size of the bat determines the diet but that certain morphological characteristics of the skull are much more important [11-13,44].
Lepidoptera are relatively large, and bats don't need to be very agile to catch these insect in comparison to smaller insects. To detect larger insects, a low frequency should normally suffice  should normally suffice to detect larger insects, however, high frequency calls are actually needed to feed efficiently on Lepidoptera. Several moth species are known to easily detect lower call frequencies and subsequently avoid the bats [16,47,48]. Other Lepidoptera species can even jam the echolocation call which makes it very difficult for the bat to detect the prey [16,49]. Call frequency was relative strongly related to Lepidoptera content, while cranial length and bodyweight did not show any relationship with the dietary content of Lepidoptera. With higher peak frequencies the amount of Lepidoptera increased and the amount of other insect reduced. Especially Rhinolophus and Hipposideros species fed most on Lepidoptera.
We have shown that bodyweight-forearm ratio and cranial morphology rather than echolocation are good predictors for the content of hard or soft insects in the diets of bats. Whereas echolocation is the one parameter that determines Lepidoptera content. These traits relate strongly to prey selectivity. While some morphological characteristics can used to quantified diets to some degree, other external variables are also important in determining the diets of insectivorous bats e.g. season, local insect community composition or geographic range [2,19,50,51]. These external factors make it difficult to only use dietary analysis for assessing food preferences in bats. Additional difficulties arise from large biases in assessing the food supply, and comparing this to the consumed food items . We propose that a combination of dietary analysis and morphological analysis is needed to make stronger inference about prey preference.
The authors would like to thank the faculty of Agriculture, Natural Resources and Environment of Naresuan University and the Cat Drop Foundation for financial support of this study. We kindly thank Mr. William Newcomb for proofreading an early draft of this manuscript.