Air-quality networks collect environmental DNA with the potential to measure biodiversity at continental scales

One of the biggest planetary challenges is the accelerating loss of biodiversity threatening ecosystem functioning on a global scale. The WWF Living Planet Report (https://livingplanet.panda.org/) estimates a 69% decline in populations since 1970. The Convention on Biological Diversity and related international treaties ask countries to monitor shifts in community composition and assess rates of species decline to quantify extant biodiversity relative to global targets¹. However, quantifying biodiversity is a challenge, and monitoring continual change is impossible at almost any scale due to a lack of standardized data and indicators^2,3. A common problem is that the required infrastructure for such global monitoring does not exist. Here, we challenge this notion by analysing environmental DNA (eDNA) captured along with particulate matter by routine ambient air quality monitoring stations in the UK. In our samples, we identified eDNA from >180 vertebrate, arthropod, plant and fungal taxa representative of local biodiversity. We contend that air monitoring networks are in fact gathering eDNA data reflecting local biodiversity on a continental scale, as a result of their routine function. In some regions, air quality samples are stored for decades, presenting the potential for high resolution biodiversity time series. With minimal modification of current protocols, this material provides the best opportunity to date for detailed monitoring of terrestrial biodiversity using an existing, replicated transnational design and it is already in operation.

Out of thin air: surveying tropical bat roosts through air sampling of eDNA

Understanding roosting behaviour is essential to bat conservation and biomonitoring, often providing the most accurate methods of assessing bat population size and health. However, roosts can be challenging to survey, e.g., physically impossible to access or presenting risks for researchers. Disturbance during monitoring can also disrupt natural bat behaviour and present material risks to the population such as disrupting hibernation cycles. One solution to this is the use of non-invasive monitoring approaches. Environmental (e)DNA has proven especially effective at detecting rare and elusive species particularly in hard-to-reach locations. It has recently been demonstrated that eDNA from vertebrates is carried in air. When collected in semi-confined spaces, this airborne eDNA can provide remarkably accurate profiles of biodiversity, even in complex tropical communities. In this study, we deploy novel airborne eDNA collection for the first time in a natural setting and use this approach to survey difficult to access potential roosts in the neotropics. Using airborne eDNA, we confirmed the presence of bats in nine out of 12 roosts. The identified species matched previous records of roost use obtained from photographic and live capture methods, thus demonstrating the utility of this approach. We also detected the presence of the white-winged vampire bat (Diaemus youngi) which had never been confirmed in the area but was long suspected based on range maps. In addition to the bats, we detected several non-bat vertebrates, including the big-eared climbing rat (Ototylomys phyllotis), which has previously been observed in and around bat roosts in our study area. We also detected eDNA from other local species known to be in the vicinity. Using airborne eDNA to detect new roosts and monitor known populations, particularly when species turnover is rapid, could maximize efficiency for surveyors while minimizing disturbance to the animals. This study presents the first applied use of airborne eDNA collection for ecological analysis moving beyond proof of concept to demonstrate a clear utility for this technology in the wild.

Interpretation and Application of Bat Diversity and Phylogeny

Brock Fenton has devoted his career to the study of the ecological adaptations of bat diversity. In this paper we describe his interest and research on the subject of bat evolution and how he has used phylogenetic hypotheses to revise our understanding of divergences and convergences of specific traits within this mammalian order. While he has always been fascinated by the evolution of echolocation, he has also written about the evolution of reproductive, morphological and behavioural traits. Recent methods of documenting diversity incorporate his lifelong love of photography. We describe taxonomic attempts to honour his contributions to bat species diversity. Brock has inspired hundreds of students and colleagues with his endless enthusiasm and generosity and inspired them with ongoing evolutionary research on bats.

Non-target effects of agri-environmental schemes on solitary bees and fungi in the United Kingdom

Agri-environmental schemes (AES) are used to enhance pollinator diversity on agricultural farms within the UK. Though the impacts of these schemes on archetypal pollinator species such as the bumblebee (Bombus) and honeybee (Apis) are well-studied, the effects on non-target bee species like solitary bees, in the same environment, are generally lacking. One goal of AES is to alter floral provision and taxonomic composition of plant communities to provide better forage for pollinators, however, this may potentially impact other ecological communities such as fungal diversity associated with plant-bee communities. Fungi are integral in these bee communities as they can impact bee species both beneficially and detrimentally. We test the hypothesis that alteration of the environment through provision of novel plant communities has non-target effects on the fungi associated with solitary bee communities. We analyse fungal diversity and ecological networks formed between fungi and solitary bees present on 15 agricultural farms in the UK using samples from brood cells. The farms were allocated to two categories, low and high management, which differ in the number of agri-environmental measures implemented. Using internal transcribed spacer metabarcoding, we identified 456 fungal taxa that interact with solitary bees. Of these, 202 (approximately 44%) could be assigned to functional groups, the majority being pathotrophic and saprotrophic species. A large proportion was Ascosphaeraceae, a family of bee-specialist fungi. We considered the connectance, nestedness, modularity, nestedness overlap and decreasing fill, linkage density and fungal generality of the farms' bee-fungi ecological networks. We found no difference in the structure of bee-fungi ecological networks between low and high management farms, suggesting floral provision by AES has no significant impact on interactions between these two taxonomic groups. However, bee emergence was lower on the low management farms compared to high management, suggesting some limited non-target effects of AES. This study characterizes the fungal community associated with solitary bees and provides evidence that floral provision through AES does not impact fungal interactions.

Cascading effects of habitat loss on ectoparasite-associated bacterial microbiomes

Suitable habitat fragment size, isolation, and distance from a source are important variables influencing community composition of plants and animals, but the role of these environmental factors in determining composition and variation of host-associated microbial communities is poorly known. In parasite-associated microbial communities, it is hypothesized that evolution and ecology of an arthropod parasite will influence its microbiome more than broader environmental factors, but this hypothesis has not been extensively tested. To examine the influence of the broader environment on the parasite microbiome, we applied high-throughput sequencing of the V4 region of 16S rRNA to characterize the microbiome of 222 obligate ectoparasitic bat flies (Streblidae and Nycteribiidae) collected from 155 bats (representing six species) from ten habitat fragments in the Atlantic Forest of Brazil. Parasite species identity is the strongest driver of microbiome composition. To a lesser extent, reduction in habitat fragment area, but not isolation, is associated with an increase in connectance and betweenness centrality of bacterial association networks driven by changes in the diversity of the parasite community. Controlling for the parasite community, bacterial network topology covaries with habitat patch area and exhibits parasite-species specific responses to environmental change. Taken together, habitat loss may have cascading consequences for communities of interacting macro- and microorgansims.

Detection rates of aphid DNA in the guts of larval hoverflies and potential links to the provision of floral resources

Aphidophagous hoverflies (Diptera, Syrphidae, Syrphinae) are common flower visitors and aphid predators in a range of flowering plants, including fruit crops. Here, we investigate whether aphid prey DNA can be detected in the gut contents of hoverfly larvae from a commercial strawberry field as a proof of concept that a molecular approach can be used to measure agricultural biocontrol. We used high-throughput sequencing (HTS) to target insect DNA and compared the resulting data to reference databases containing aphid and hoverfly DNA sequences. We explored what impact incorporating wildflowers within polythene-clad tunnels may have on aphid DNA detection rates in hoverfly larvae. In a randomized block experiment, coriander (Coriandrum sativum), field forget-me-not (Myosotis arvensis) and corn mint (Mentha arvensis) plants were inserted in rows of strawberries. Their effect on aphid DNA detection rates was assessed. Aphid DNA was found in 55 of 149 specimens (37%) validating the method in principle for measuring agricultural services provided by hoverflies. Interestingly, detection rates were higher near plots with forget-me-not than plots with coriander, though detection rates in control plots did not differ significantly from either wildflower species. These findings confirm that hoverflies predate aphids in UK strawberry fields, and that HTS is a viable method of identifying aphid DNA in predatory hoverflies. We comment on the need for further method development to narrow down identifications of both predator and prey. We furthermore provide some evidence that there is an effect of intercropping strawberry crops with wildflowers which may affect aphid consumption in hoverfly larvae.

Shunning the scoop: Sidestepping the race to publish

What happens when a researcher finds out that research very similar to their own is already being conducted? What if they find out that the said research is also very close to being published? First, there is probably anxiety and panic. Maybe, there are frantic calls to collaborators. Perhaps Twitter rants about the phenomenon of scooping that plagues all researchers, especially those early-career researchers who often feel they are in a race to get their best work out to the world.

Measuring biodiversity from DNA in the air

The crisis of declining biodiversity¹ exceeds our current ability to monitor changes in ecosystems. Rapid terrestrial biomonitoring approaches are essential to quantify the causes and consequences of global change. Environmental DNA² has revolutionized aquatic ecology,³ permitting population monitoring⁴ and remote diversity assessments matching or outperforming conventional methods of community sampling.^3-5 Despite this model, similar methods have not been widely adopted in terrestrial ecosystems. Here, we demonstrate that DNA from terrestrial animals can be filtered, amplified, and then sequenced from air samples collected in natural settings representing a powerful tool for terrestrial ecology. We collected air samples at a zoological park, where spatially confined non-native species allowed us to track DNA sources. We show that DNA can be collected from air and used to identify species and their ecological interactions. Air samples contained DNA from 25 species of mammals and birds, including 17 known terrestrial resident zoo species. We also identified food items from air sampled in enclosures and detected taxa native to the local area, including the Eurasian hedgehog, endangered in the United Kingdom. Our data demonstrate that airborne eDNA concentrates around recently inhabited areas but disperses away from sources, suggesting an ecology to airborne eDNA and the potential for sampling at a distance. Our findings demonstrate the profound potential of air as a source of DNA for global terrestrial biomonitoring.

Selective Logging Shows No Impact on the Dietary Breadth of a Generalist Bat Species: The Fawn Leaf-Nosed Bat (Hipposideros cervinus)

Logging activities degrade forest habitats across large areas of the tropics, but the impacts on trophic interactions that underpin forest ecosystems are poorly understood. DNA metabarcoding provides an invaluable tool to investigate such interactions, allowing analysis at a far greater scale and resolution than has previously been possible. We analysed the diet of the insectivorous fawn leaf-nosed bat Hipposideros cervinus across a forest disturbance gradient in Borneo, using a dataset of ecological interactions from an unprecedented number of bat-derived faecal samples. Bats predominantly consumed insects from the orders Lepidoptera, Diptera, Blattodea, and Coleoptera, and the taxonomic composition of their diet remained relatively consistent across sites regardless of logging disturbance. There was little difference in the richness of prey consumed per-bat in each logging treatment, indicating potential resilience of this species to habitat degradation. In fact, bats consumed a high richness of prey items, and intensive sampling is needed to reliably compare feeding ecology over multiple sites. Multiple bioinformatic parameters were used, to assess how they altered our perception of sampling completeness. While parameter choice altered estimates of completeness, a very high sampling effort was always required to detect the entire prey community.

Molecular diet analysis of the marine fish-eating bat (Myotis vivesi) and potential mercury exposure

Mercury is a toxic element acquired by animals through feeding which can accumulate within food chains through biomagnification. This possesses particular risks to higher trophic levels and may unduly impact marine foraging species or individuals. The fish-eating bat (Myotis vivesi Menegaux, 1901) inhabits islands in the Gulf of California and can act as a predator in the marine environment. A predominantly marine diet and a high trophic position increase the risk of mercury exposure owing to increased bioaccumulation. Using molecular techniques to reconstruct diet, we show that M. vivesi regularly feeds on small fishes and crustaceans, particularly on the Californian anchovy (Engraulis mordax Girard, 1854) and a krill species (Nyctiphanes simplex Hansen, 1911). Additionally, we identify significant interannual variation in diet composition within this population, but measured levels of total mercury in faecal samples were not related to dietary diversity or trophic level.

Dung beetles as samplers of mammals in Malaysian Borneo-a test of high throughput metabarcoding of iDNA

Invertebrate-derived DNA (iDNA) sampling in biodiversity surveys is becoming increasingly widespread, with most terrestrial studies relying on DNA derived from the gut contents of blood-feeding invertebrates, such as leeches and mosquitoes. Dung beetles (superfamily Scarabaeoidea) primarily feed on the faecal matter of terrestrial vertebrates and offer several potential benefits over blood-feeding invertebrates as samplers of vertebrate DNA. Importantly, these beetles can be easily captured in large numbers using simple, inexpensive baited traps, are globally distributed, and occur in a wide range of habitats. To build on the few existing studies demonstrating the potential of dung beetles as sources of mammalian DNA, we subjected the large-bodied, Bornean dung beetle (Catharsius renaudpauliani) to a controlled feeding experiment. We analysed DNA from gut contents at different times after feeding using qPCR techniques. Here, we first describe the window of DNA persistence within a dung beetle digestive tract. We found that the ability to successfully amplify cattle DNA decayed over relatively short time periods, with DNA copy number decreasing by two orders of magnitude in just 6 h. In addition, we sampled communities of dung beetles from a lowland tropical rainforest in Sabah, Malaysia, in order to test whether it is possible to identify vertebrate sequences from dung beetle iDNA. We sequenced both the gut contents from large dung beetle species, as well as whole communities of smaller beetles. We successfully identified six mammalian species from our samples, including the bearded pig (Sus barbatus) and the sambar deer (Rusa unicolor)—both vulnerable species on the IUCN red list. Our results represent the first use of dung beetle iDNA to sample Southeast Asian vertebrate fauna, and highlight the potential for dung beetle iDNA to be used in future biodiversity monitoring surveys.

Altered structure of bat-prey interaction networks in logged tropical forests revealed by metabarcoding

Habitat degradation is pervasive across the tropics and is particularly acute in Southeast Asia, with major implications for biodiversity. Much research has addressed the impact of degradation on species diversity; however, little is known about how ecological interactions are altered, including those that constitute important ecosystem functions such as consumption of herbivores. To examine how rainforest degradation alters trophic interaction networks, we applied DNA metabarcoding to construct interaction networks linking forest-dwelling insectivorous bat species and their prey, comparing old-growth forest and forest degraded by logging in Sabah, Borneo. Individual bats in logged rainforest consumed a lower richness of prey than those in old-growth forest. As a result, interaction networks in logged forests had a less nested structure. These network structures were associated with reduced network redundancy and thus increased vulnerability to perturbations in logged forests. Our results show how ecological interactions change between old-growth and logged forests, with potentially negative implications for ecosystem function and network stability.

Exploratory analysis reveals arthropod consumption in 10 lemur species using DNA metabarcoding

Arthropods (insects, spiders, etc.) can fulfill major nutritional requirements for primates, particularly in terms of proteins, fats, vitamins, and minerals. Yet, for many primate species we know very little about the frequency and importance of arthropod consumption. Traditional methods for arthropod prey identification, such as behavioral observations and fecal dissections, offer limited taxonomic resolution and, as a result, underestimate true diversity. Metabarcoding arthropod DNA from primate fecal samples provides a promising but underused alternative. Here, we inventoried arthropod prey diversity in wild lemurs by sequencing two regions of the CO1 gene. Samples were collected opportunistically from 10 species of lemurs inhabiting three national parks in southern Madagascar using a combination of focal animal follows and live trapping. In total, we detected arthropod DNA in 98 of the 170 fecal samples analyzed. Although all lemur species included in these analyses showed evidence of arthropod consumption, those within the family Cheirogaleidae appeared to consume the highest frequency and diversity of arthropods. To our knowledge, this study presents the first evidence of arthropod consumption in Phaner pallescens, Avahi peyrierasi, and Propithecus verreauxi, and identifies 32 families of arthropods as probable food items that have not been published as lemur dietary items to date. Our study emphasizes the importance of arthropods as a nutritional source and the role DNA metabarcoding can play in elucidating an animal's diet.

eDNAir: proof of concept that animal DNA can be collected from air sampling

Environmental DNA (eDNA) is one of the fastest developing tools for species biomonitoring and ecological research. However, despite substantial interest from research, commercial and regulatory sectors, it has remained primarily a tool for aquatic systems with a small amount of work in substances such as soil, snow and rain. Here we demonstrate that eDNA can be collected from air and used to identify mammals. Our proof of concept successfully demonstrated that eDNA sampled from air contained mixed templates which reflect the species known to be present within a confined space and that this material can be accessed using existing sampling methods. We anticipate this demonstration will initiate a much larger research programme in terrestrial airDNA sampling and that this may rapidly advance biomonitoring approaches. Lastly, we outline these and potential related applications we expect to benefit from this development.

Trophic resource partitioning drives fine-scale coexistence in cryptic bat species

Understanding the processes that enable species coexistence has important implications for assessing how ecological systems will respond to global change. Morphology and functional similarity increase the potential for competition, and therefore, co-occurring morphologically similar but genetically unique species are a good model system for testing coexistence mechanisms. We used DNA metabarcoding and high-throughput sequencing to characterize for the first time the trophic ecology of two recently described cryptic bat species with parapatric ranges, Myotis escalerai and Myotis crypticus. We collected fecal samples from allopatric and sympatric regions and from syntopic and allotopic locations within the sympatric region to describe the diets both taxonomically and functionally and compare prey consumption with prey availability. The two bat species had highly similar diets characterized by high arthropod diversity, particularly Lepidoptera, Diptera and Araneae, and a high proportion of prey that is not volant at night, which points to extensive use of gleaning. Diet overlap at the prey item level was lower in syntopic populations, supporting trophic shift under fine-scale co-occurrence. Furthermore, the diet of M. escalerai had a marginally lower proportion of not nocturnally volant prey in syntopic populations, suggesting that the shift in diet may be driven by a change in foraging mode. Our findings suggest that fine-scale coexistence mechanisms can have implications for maintaining broad-scale diversity patterns. This study highlights the importance of including both allopatric and sympatric populations and choosing meaningful spatial scales for detecting ecological patterns. We conclude that a combination of high taxonomic resolution with a functional approach helps identify patterns of niche shift.

Leech blood-meal invertebrate-derived DNA reveals differences in Bornean mammal diversity across habitats

The application of metabarcoding to environmental and invertebrate‐derived DNA (eDNA and iDNA) is a new and increasingly applied method for monitoring biodiversity across a diverse range of habitats. This approach is particularly promising for sampling in the biodiverse humid tropics, where rapid land‐use change for agriculture means there is a growing need to understand the conservation value of the remaining mosaic and degraded landscapes. Here we use iDNA from blood‐feeding leeches (Haemadipsa picta) to assess differences in mammalian diversity across a gradient of forest degradation in Sabah, Malaysian Borneo. We screened 557 individual leeches for mammal DNA by targeting fragments of the 16S rRNA gene and detected 14 mammalian genera. We recorded lower mammal diversity in the most heavily degraded forest compared to higher quality twice logged forest. Although the accumulation curves of diversity estimates were comparable across these habitat types, diversity was higher in twice logged forest, with more taxa of conservation concern. In addition, our analysis revealed differences between the community recorded in the heavily logged forest and that of the twice logged forest. By revealing differences in mammal diversity across a human‐modified tropical landscape, our study demonstrates the value of iDNA as a noninvasive biomonitoring approach in conservation assessments.

Assessing the impact of taxon resolution on network structure

Constructing ecological networks has become an indispensable approach in understanding how different taxa interact. However, the methods used to generate data in network research vary widely among studies, potentially limiting our ability to compare results meaningfully. In particular, methods of classifying nodes vary in their precision, likely altering the architecture of the network studied. For example, rather than being classified as Linnaean species, taxa are regularly assigned to morphospecies in observational studies, or to molecular operational taxonomic units (MOTUs) in molecular studies, with the latter defined based on an arbitrary threshold of sequence similarity. Although the use of MOTUs in ecological networks holds great potential, especially for allowing rapid construction of large data sets of interactions, it is unclear how the choice of clustering threshold can influence the conclusions obtained. To test the impact of taxonomic precision on network architecture, we obtained and analyzed 16 data sets of ecological interactions, inferred from metabarcoding and observations. Our comparisons of networks constructed under a range of sequence thresholds for assigning taxa demonstrate that even small changes in node resolution can cause wide variation in almost all key metric values. Moreover, relative values of commonly used metrics such as robustness were seen to fluctuate continuously with node resolution, thereby potentially causing error in conclusions drawn when comparing multiple networks. In observational networks, we found that changing node resolution could, in some cases, lead to substantial changes to measurements of network topology. Overall, our findings highlight the importance of classifying nodes to the greatest precision possible, and demonstrate the need for caution when comparing networks that differ with respect to node resolution, even where taxonomic groups and interaction types are similar. In such cases, we recommend that comparisons of networks should focus on relative differences rather than absolute values between the networks studied.

Differential effects of fertilisers on pollination and parasitoid interaction networks

Grassland fertilisation drives non-random plant loss resulting in areas dominated by perennial grass species. How these changes cascade through linked trophic levels, however, is not well understood. We studied how grassland fertilisation propagates change through the plant assemblage into the plant-flower-visitor, plant-leaf miner and leaf miner-parasitoid networks using a year's data collection from a long-term grassland fertiliser application experiment. Our experiment had three fertiliser treatments each applied to replicate plots 15 m² in size: mineral fertiliser, farmyard manure, and mineral fertiliser and farmyard manure combined, along with a control of no fertiliser. The combined treatment had the most significant impact, and both plant species richness and floral abundance decreased with the addition of fertiliser. While insect species richness was unaffected by fertiliser treatment, fertilised plots had a significantly higher abundance of leaf miners and parasitoids and a significantly lower abundance of bumblebees. The plant-flower-visitor and plant-herbivore networks showed higher values of vulnerability and lower modularity with fertiliser addition, while leaf miner-parasitoid networks showed a rise in generality. The different groups of insects were impacted by fertilisers to varying degrees: while the effect on abundance was the highest for leaf miners, the vulnerability and modularity of flower-visitor networks was the most affected. The impact on the abundance of leaf miners was positive and three times higher than the impact on parasitoids, and the impact on bumblebee abundance was negative and double the magnitude of impact on flower abundance. Overall, our results show that while insect species richness was unaffected by fertilisers, network structure changed significantly as the replacement of forbs by grasses resulted in changes in relative abundance across trophic levels, with the direction of change depending on the type of network. Synthesis. By studying multiple networks simultaneously, we were able to rank the relative impact of habitat change on the different groups of species within the community. This provided a more holistic picture of the impact of agricultural intensification and provides useful information when deciding on priorities for mitigation.

Wing morphology predicts individual niche specialization in Pteronotus mesoamericanus (Mammalia: Chiroptera)

Morphological variation between individuals can increase niche segregation and decrease intraspecific competition when heterogeneous individuals explore their environment in different ways. Among bat species, wing shape correlates with flight maneuverability and habitat use, with species that possess broader wings typically foraging in more cluttered habitats. However, few studies have investigated the role of morphological variation in bats for niche partitioning at the individual level. To determine the relationship between wing shape and diet, we studied a population of the insectivorous bat species Pteronotus mesoamericanus in the dry forest of Costa Rica. Individual diet was resolved using DNA metabarcoding, and bat wing shape was assessed using geometric morphometric analysis. Inter-individual variation in wing shape showed a significant relationship with both dietary dissimilarity based on Bray-Curtis estimates, and nestedness derived from an ecological network. Individual bats with broader and more rounded wings were found to feed on a greater diversity of arthropods (less nested) in comparison to individuals with triangular and pointed wings (more nested). We conclude that individual variation in bat wing morphology can impact foraging efficiency leading to the observed overall patterns of diet specialization and differentiation within the population.

Chemical characterisation of potential pheromones from the shoulder gland of the Northern yellow-shouldered-bat, Sturnira parvidens (Phyllostomidae: Stenodermatinae)

Bats of the genus Sturnira (Family Phyllostomidae) are characterised by shoulder glands that are more developed in reproductively mature adult males. The glands produce a waxy secretion that accumulates on the fur around the gland, dyeing the fur a dark colour and giving off a pungent odour. These shoulder glands are thought to play a role in their reproductive behaviour. Using gas chromatography–mass spectrometry, we analysed solvent extracts of fur surrounding the shoulder gland in the northern-shouldered bat, Sturnira parvidens to (i) characterise the chemical composition of shoulder gland secretions for the first time, and (ii) look for differences in chemical composition among and between adult males, sub-adult/juvenile males and adult females. Fur solvent extracts were analysed as liquids and also further extracted using headspace solid-phase microextraction to identify volatile components in the odour itself. Odour fingerprint analysis using non-metric multidimensional scaling plots and multivariate analysis revealed clear and significant differences (P < 0.001) between adult males vs both juvenile males and adult females. The chemical components of the shoulder gland secretion included terpenes and phenolics, together with alcohols and esters, most likely derived from the frugivorous diet of the bat. Many of the compounds identified were found exclusively or in elevated quantities among adult (reproductive) males compared with adult females and non-reproductive (juvenile) males. This strongly suggests a specific role in male–female attraction although a function in male–male competition and/or species recognition is also possible.

Molecular diet analysis finds an insectivorous desert bat community dominated by resource sharing despite diverse echolocation and foraging strategies

Interspecific differences in traits can alter the relative niche use of species within the same environment. Bats provide an excellent model to study niche use because they use a wide variety of behavioral, acoustic, and morphological traits that may lead to multi-species, functional groups. Predatory bats have been classified by their foraging location (edge, clutter, open space), ability to use aerial hawking or substrate gleaning and echolocation call design and flexibility, all of which may dictate their prey use. For example, high frequency, broadband calls do not travel far but offer high object resolution while high intensity, low frequency calls travel further but provide lower resolution. Because these behaviors can be flexible, four behavioral categories have been proposed: (a) gleaning, (b) behaviorally flexible (gleaning and hawking), (c) clutter-tolerant hawking, and (d) open space hawking. Many recent studies of diet in bats use molecular tools to identify prey but mainly focus on one or two species in isolation; few studies provide evidence for substantial differences in prey use despite the many behavioral, acoustic, and morphological differences. Here, we analyze the diet of 17 sympatric species in the Chihuahuan desert and test the hypothesis that peak echolocation frequency and behavioral categories are linked to differences in diet. We find no significant correlation between dietary richness and echolocation peak frequency though it spanned close to 100 kHz across species. Our data, however, suggest that bats which use both gleaning and hawking strategies have the broadest diets and are most differentiated from clutter-tolerant aerial hawking species.

Approaches to integrating genetic data into ecological networks

As molecular tools for assessing trophic interactions become common, research is increasingly focused on the construction of interaction networks. Here, we demonstrate three key methods for incorporating DNA data into network ecology and discuss analytical considerations using a model consisting of plants, insects, bats and their parasites from the Costa Rica dry forest. The simplest method involves the use of Sanger sequencing to acquire long sequences to validate or refine field identifications, for example of bats and their parasites, where one specimen yields one sequence and one identification. This method can be fully quantified and resolved and these data resemble traditional ecological networks. For more complex taxonomic identifications, we target multiple DNA loci, for example from a seed or fruit pulp sample in faeces. These networks are also well resolved but gene targets vary in resolution and quantification is difficult. Finally, for mixed templates such as faecal contents of insectivorous bats, we use DNA metabarcoding targeting two sequence lengths (157 and 407 bp) of one gene region and a MOTU, BLAST and BIN association approach to resolve nodes. This network type is complex to generate and analyse, and we discuss the implications of this type of resolution on network analysis. Using these data, we construct the first molecular-based network of networks containing 3,304 interactions between 762 nodes of eight trophic functions and involving parasitic, mutualistic and predatory interactions. We provide a comparison of the relative strengths and weaknesses of these data types in network ecology.

Temporal changes in arthropod activity in tropical anthropogenic forests

Arthropod communities in the tropics are increasingly impacted by rapid changes in land use. Because species showing distinct seasonal patterns of activity are thought to be at higher risk of climate-related extirpation, global warming is generally considered a lower threat to arthropod biodiversity in the tropics than in temperate regions. To examine changes associated with land use and weather variables in tropical arthropod communities, we deployed Malaise traps at three major anthropogenic forests (secondary reserve forest, oil palm forest, and urban ornamental forest (UOF)) in Peninsular Malaysia and collected arthropods continuously for 12 months. We used metabarcoding protocols to characterize the diversity within weekly samples. We found that changes in the composition of arthropod communities were significantly associated with maximum temperature in all the three forests, but shifts were reversed in the UOF compared with the other forests. This suggests arthropods in forests in Peninsular Malaysia face a double threat: community shifts and biodiversity loss due to exploitation and disturbance of forests which consequently put species at further risk related to global warming. We highlight the positive feedback mechanism of land use and temperature, which pose threats to the arthropod communities and further implicates ecosystem functioning and human well-being. Consequently, conservation and mitigation plans are urgently needed.

Spatiotemporal and demographic variation in the diet of New Zealand lesser short-tailed bats (Mystacina tuberculata)

Variation in the diet of generalist insectivores can be affected by site-specific traits including weather, habitat, and season, as well as demographic traits such as reproductive status and age. We used molecular methods to compare diets of three distinct New Zealand populations of lesser short-tailed bats, Mystacina tuberculata. Summer diets were compared between a southern cold-temperate (Eglinton) and a northern population (Puroera). Winter diets were compared between Pureora and a subtropical offshore island population (Hauturu). This also permitted seasonal diet comparisons within the Pureora population. Lepidoptera and Diptera accounted for >80% of MOTUs identified from fecal matter at each site/season. The proportion of orders represented within prey and the Simpson diversity index, differed between sites and seasons within the Pureora population. For the Pureora population, the value of the Simpson diversity index was higher in summer than winter and was higher in Pureora compared to Eglinton. Summer Eglinton samples revealed that juvenile diets appeared to be more diverse than other demographic groups. Lactating females had the lowest dietary diversity during summer in Pureora. In Hauturu, we found a significant negative relationship between mean ambient temperature and prey richness. Our data suggest that M. tuberculata incorporate a narrower diversity of terrestrial insects than previously reported. This provides novel insights into foraging behavior and ecological interactions within different habitats. Our study is the first from the Southern Hemisphere to use molecular techniques to examine spatiotemporal variation in the diet of a generalist insectivore that inhabits a contiguous range with several habitat types and climates.

Counting with DNA in metabarcoding studies: How should we convert sequence reads to dietary data?

Advances in DNA sequencing technology have revolutionized the field of molecular analysis of trophic interactions, and it is now possible to recover counts of food DNA sequences from a wide range of dietary samples. But what do these counts mean? To obtain an accurate estimate of a consumer's diet should we work strictly with data sets summarizing frequency of occurrence of different food taxa, or is it possible to use relative number of sequences? Both approaches are applied to obtain semi-quantitative diet summaries, but occurrence data are often promoted as a more conservative and reliable option due to taxa-specific biases in recovery of sequences. We explore representative dietary metabarcoding data sets and point out that diet summaries based on occurrence data often overestimate the importance of food consumed in small quantities (potentially including low-level contaminants) and are sensitive to the count threshold used to define an occurrence. Our simulations indicate that using relative read abundance (RRA) information often provides a more accurate view of population-level diet even with moderate recovery biases incorporated; however, RRA summaries are sensitive to recovery biases impacting common diet taxa. Both approaches are more accurate when the mean number of food taxa in samples is small. The ideas presented here highlight the need to consider all sources of bias and to justify the methods used to interpret count data in dietary metabarcoding studies. We encourage researchers to continue addressing methodological challenges and acknowledge unanswered questions to help spur future investigations in this rapidly developing area of research.

Mercury bioaccumulation in bats reflects dietary connectivity to aquatic food webs

Mercury (Hg) is a persistent and widespread heavy metal with neurotoxic effects in wildlife. While bioaccumulation of Hg has historically been studied in aquatic food webs, terrestrial consumers can become contaminated with Hg when they feed on aquatic organisms (e.g., emergent aquatic insects, fish, and amphibians). However, the extent to which dietary connectivity to aquatic ecosystems can explain patterns of Hg bioaccumulation in terrestrial consumers has not been well studied. Bats (Order: Chiroptera) can serve as a model system for illuminating the trophic transfer of Hg given their high dietary diversity and foraging links to both aquatic and terrestrial food webs. Here we quantitatively characterize the dietary correlates of long-term exposure to Hg across a diverse local assemblage of bats in Belize and more globally across bat species from around the world with a comparative analysis of hair samples. Our data demonstrate considerable interspecific variation in hair total Hg concentrations in bats that span three orders of magnitude across species, ranging from 0.04 mg/kg in frugivorous bats (Artibeus spp.) to 145.27 mg/kg in the piscivorous Noctilio leporinus. Hg concentrations showed strong phylogenetic signal and were best explained by dietary connectivity of bat species to aquatic food webs. Our results highlight that phylogeny can be predictive of Hg concentrations through similarity in diet and how interspecific variation in feeding strategies influences chronic exposure to Hg and enables movement of contaminants from aquatic to terrestrial ecosystems.

The effects of pastoral intensification on the feeding interactions of generalist predators in streams

Land‐use change can alter trophic interactions with wide‐ranging functional consequences, yet the consequences for aquatic food webs have been little studied. In part, this may reflect the challenges of resolving the diets of aquatic organisms using classical gut contents analysis, especially for soft‐bodied prey. We used next‐generation sequencing to resolve prey use in nearly 400 individuals of two predatory invertebrates (the Caddisfly, Rhyacophila dorsalis, and the Stonefly Dinocras cephalotes) in streams draining land with increasingly intensive livestock farming. Rhyacophila dorsalis occurred in all streams, whereas D. cephalotes was restricted to low intensities, allowing us to test whether: (i) apparent sensitivity to agriculture in the latter species reflects a more specialized diet and (ii) diet in R. dorsalis varied between sites with and without D. cephalotes. DNA was extracted from dissected gut contents, amplified without blocking probes and sequenced using Ion Torrent technology. Both predators were generalists, consuming 30 prey taxa with a preference for taxa that were abundant in all streams or that increased with intensification. Where both predators were present, their diets were nearly identical, and R. dorsalis's diet was virtually unchanged in the absence of D. cephalotes. The loss of D. cephalotes from more intensive sites was probably due to physicochemical stressors, such as sedimentation, rather than to dietary specialization, although wider biotic factors (e.g., competition with other predatory taxa) could not be excluded. This study provides a uniquely detailed description of predator diets along a land‐use intensity gradient, offering new insights into how anthropogenic stressors affect stream communities.

Barcoding the food chain: from Sanger to high-throughput sequencing

Society faces the complex challenge of supporting biodiversity and ecosystem functioning, while ensuring food security by providing safe traceable food through an ever-more-complex global food chain. The increase in human mobility brings the added threat of pests, parasites, and invaders that further complicate our agro-industrial efforts. DNA barcoding technologies allow researchers to identify both individual species, and, when combined with universal primers and high-throughput sequencing techniques, the diversity within mixed samples (metabarcoding). These tools are already being employed to detect market substitutions, trace pests through the forensic evaluation of trace “environmental DNA”, and to track parasitic infections in livestock. The potential of DNA barcoding to contribute to increased security of the food chain is clear, but challenges remain in regulation and the need for validation of experimental analysis. Here, we present an overview of the current uses and challenges of applied DNA barcoding in agriculture, from agro-ecosystems within farmland to the kitchen table.

The effects of parameter choice on defining molecular operational taxonomic units and resulting ecological analyses of metabarcoding data

The combination of DNA barcoding and high-throughput (next-generation) sequencing (metabarcoding) provides many promises but also serious challenges. Generating a reliable comparable estimate of biodiversity remains a central challenge to the application of the technology. Many approaches have been used to turn millions of sequences into distinct taxonomic units. However, the extent to which these methods impact the outcome of simple ecological analyses is not well understood. Here we performed a simple analysis of dietary overlap by skinks and shrews on Ile Aux Aigrettes, Mauritius. We used a combination of filtering thresholds and clustering algorithms on a COI metabarcoding dataset and demonstrate that all bioinformatics parameters will have interacting effects on molecular operational taxonomic unit (MOTU) recovery rates. These effects generated estimates covering two orders of magnitude. However, the effect on a simple ecological analysis was not large and, despite the wide variation in estimates of niche overlap, the same ecological conclusion was drawn in most cases. We advise that a conservative clustering programme coupled with larger sequence divergences to define a cluster, the removal of singletons, rigorous length filtering, and stringent match criteria for Molecular Identifier tags are preferable to avoid MOTU inflation and that the same parameters be used in all comparative analyses.

Acoustic shadows help gleaning bats find prey, but may be defeated by prey acoustic camouflage on rough surfaces

Perceptual abilities of animals, like echolocating bats, are difficult to study because they challenge our understanding of non-visual senses. We used novel acoustic tomography to convert echoes into visual representations and compare these cues to traditional echo measurements. We provide a new hypothesis for the echo-acoustic basis of prey detection on surfaces. We propose that bats perceive a change in depth profile and an 'acoustic shadow' cast by prey. The shadow is more salient than prey echoes and particularly strong on smooth surfaces. This may explain why bats look for prey on flat surfaces like leaves using scanning behaviour. We propose that rather than forming search images for prey, whose characteristics are unpredictable, predators may look for disruptions to the resting surface (acoustic shadows). The fact that the acoustic shadow is much fainter on rougher resting surfaces provides the first empirical evidence for 'acoustic camouflage' as an anti-predator defence mechanism.

Molecular detection of trophic interactions: emerging trends, distinct advantages, significant considerations and conservation applications

The emerging field of ecological genomics contains several broad research areas. Comparative genomic and conservation genetic analyses are providing great insight into adaptive processes, species bottlenecks, population dynamics and areas of conservation priority. Now the same technological advances in high-throughput sequencing, coupled with taxonomically broad sequence repositories, are providing greater resolution and fundamentally new insights into functional ecology. In particular, we now have the capacity in some systems to rapidly identify thousands of species-level interactions using non-invasive methods based on the detection of trace DNA. This represents a powerful tool for conservation biology, for example allowing the identification of species with particularly inflexible niches and the investigation of food-webs or interaction networks with unusual or vulnerable dynamics. As they develop, these analyses will no doubt provide significant advances in the field of restoration ecology and the identification of appropriate locations for species reintroduction, as well as highlighting species at ecological risk. Here, I describe emerging patterns that have come from the various initial model systems, the advantages and limitations of the technique and key areas where these methods may significantly advance our empirical and applied conservation practices.

Diet of the insectivorous bat Pipistrellus nathusii during autumn migration and summer residence

Migration is widespread among vertebrates, yet bat migration has received little attention and only in the recent decades has a better understanding of it been gained. Migration can cause significant changes in behaviour and physiology, due to increasing energy demands and aerodynamic constraints. Dietary shifts, for example, have been shown to occur in birds before onset of migration. For bats, it is not known if a change in diet occurs during migration, although breeding season-related dietary preference has been documented. It is known that a diet rich in fats and the accumulation of fat deposits do increase the flight range of migratory bats. Some bat species can be regarded as long-distance migrants, covering up to 2000 km between summer and winter roosting areas. Pipistrellus nathusii (Vespertilionidae), a European long-distant migrant, travels each year along the Baltic Sea from north-eastern Europe to hibernate in central and southern Europe. This study presents data on the dietary habits of migrating Pipistrellus nathusii compared with those during the breeding season. We analysed faecal samples from bats on fall migration caught at the Ornithological Field Station in Pape, Latvia and from samples collected in North-Latvian summer roosts. We applied both morphological identification and molecular methods, as morphological methods also recognize life stages of prey and can contribute frequency data. The diets of bats on migration and breeding bats were similar, with Diptera and Lepidoptera comprising the major prey categories. However, certain prey groups could be explained by the different hunting habitats exploited during migration vs. summer residence.

Dietary competition between the alien Asian Musk Shrew (Suncus murinus) and a re-introduced population of Telfair's Skink (Leiolopisma telfairii)

Re-introduction of rare species to parts of their historical range is becoming increasingly important as a conservation strategy. Telfair's Skinks (Leiolopisma telfairii), once widespread on Mauritius, were until recently found only on Round Island. There it is vulnerable to stochastic events, including the introduction of alien predators that may either prey upon it or compete for food resources. Consequently, skinks have been introduced to Ile aux Aigrettes, another small Mauritian island that has been cleared of rats. However, the island has been invaded by Asian Musk Shrews (Suncus murinus), a commensal species spread by man well beyond its natural Asian range. Our aim was to use next-generation sequencing to analyse the diets of the shrews and skinks to look for niche competition. DNA was extracted from skink faeces and from the stomach contents of shrews. Application of shrew- and skink-specific primers revealed no mutual predation. The DNA was then amplified using general invertebrate primers with tags to identify individual predators, and then sequenced by 454 pyrosequencing. 119 prey MOTUs (molecular taxonomic units) were isolated, although none could be identified to species. Seeding of cladograms with known sequences allowed higher taxonomic assignments in some cases. Although most MOTUs were not shared by shrews and skinks, Pianka's niche overlap test showed significant prey overlap, suggesting potentially strong competition where food resources are limited. These results suggest that removal of the shrews from the island should remain a priority.

The promise of genomics in the study of plant-pollinator interactions

Flowers exist in exceedingly complex fitness landscapes, in which subtle variation in each trait can affect the pollinators, herbivores and pleiotropically linked traits in other plant tissues. A whole-genome approach to flower evolution will help our understanding of plant-pollinator interactions.

Diversification and reproductive isolation: cryptic species in the only New World high-duty cycle bat, Pteronotus parnellii

BACKGROUND

Molecular techniques are increasingly employed to recognize the presence of cryptic species, even among commonly observed taxa. Previous studies have demonstrated that bats using high-duty cycle echolocation may be more likely to speciate quickly. Pteronotus parnellii is a widespread Neotropical bat and the only New World species to use high-duty cycle echolocation, a trait otherwise restricted to Old World taxa. Here we analyze morphological and acoustic variation and genetic divergence at the mitochondrial COI gene, the 7th intron region of the y-linked Dby gene and the nuclear recombination-activating gene 2, and provide extensive evidence that P. parnellii is actually a cryptic species complex.

RESULTS

Central American populations form a single species while three additional species exist in northern South America: one in Venezuela, Trinidad and western Guyana and two occupying sympatric ranges in Guyana and Suriname. Reproductive isolation appears nearly complete (only one potential hybrid individual found). The complex likely arose within the last ~6 million years with all taxa diverging quickly within the last ~1-2 million years, following a pattern consistent with the geological history of Central and northern South America. Significant variation in cranial measures and forearm length exists between three of the four groups, although no individual morphological character can discriminate these in the field. Acoustic analysis reveals small differences (5-10 kHz) in echolocation calls between allopatric cryptic taxa that are unlikely to provide access to different prey resources but are consistent with divergence by drift in allopatric species or through selection for social recognition.

CONCLUSIONS

This unique approach, considering morphological, acoustic and multi-locus genetic information inherited maternally, paternally and bi-parentally, provides strong support to conclusions about the cessation of gene flow and degree of reproductive isolation of these cryptic species.