Using faecal metabarcoding to examine consumption of crop pests and beneficial arthropods in communities of generalist avian insectivores

Generalist insectivorous birds can provide ecosystem services in agricultural landscapes by consuming arthropod pests, or they can provide disservices when they consume beneficial arthropods. To examine bird impacts on arthropod communities, including pest control services, we need to know which arthropods birds commonly consume. Faecal metabarcoding is an emerging technique that can be used to identify prey from faecal samples, often to the species level. We used faecal metabarcoding to study diets of birds inhabiting the ecotone between soybean fields and adjacent grasslands in a largely agricultural landscape in Illinois, USA, during the summer of 2017. Whereas previous studies have used faecal metabarcoding to compare bird diets among species or among capture sites, we analysed samples from multiple species within a community at replicate sites. We collected and sequenced DNA from 132 faecal samples from 25 bird species captured at six sites. We found that birds consumed an extremely large and varied diet that differed among both species and sites, suggesting that birds were consuming prey opportunistically as available at each site. Of the nine most commonly detected prey species, three are known pests of soybeans. Bird diets also contained significantly more species of herbivorous prey than natural enemies. Finally, we discovered that American Goldfinches Spinus tristis, a highly granivorous species, may consume arthropods more frequently than expected and thus may provide ecosystem services in agricultural landscapes. Our study demonstrates that birds within this system consume a large variety of prey, suggesting that they may be able to respond quickly to pest outbreaks and contribute to agricultural resiliency.     

DNA metabarcoding multiplexing and validation of data accuracy for diet assessment: application to omnivorous diet

Ecological understanding of the role of consumer–resource interactions in natural food webs is limited by the difficulty of accurately and efficiently determining the complex variety of food types animals have eaten in the field. We developed a method based on DNA metabarcoding multiplexing and next‐generation sequencing to uncover different taxonomic groups of organisms from complex diet samples. We validated this approach on 91 faeces of a large omnivorous mammal, the brown bear, using DNA metabarcoding markers targeting the plant, vertebrate and invertebrate components of the diet. We included internal controls in the experiments and performed PCR replication for accuracy validation in postsequencing data analysis. Using our multiplexing strategy, we significantly simplified the experimental procedure and accurately and concurrently identified different prey DNA corresponding to the targeted taxonomic groups, with ≥60% of taxa of all diet components identified to genus/species level. The systematic application of internal controls and replication was a useful and simple way to evaluate the performance of our experimental procedure, standardize the selection of sequence filtering parameters for each marker data and validate the accuracy of the results. Our general approach can be adapted to the analysis of dietary samples of various predator species in different ecosystems, for a number of conservation and ecological applications entailing large‐scale population level diet assessment through cost‐effective screening of multiple DNA metabarcodes, and the detection of fine dietary variation among samples or individuals and of rare food items.     

A from‐benchtop‐to‐desktop workflow for validating HTS data and for taxonomic identification in diet metabarcoding studies

The main objective of this work was to develop and validate a robust and reliable “from‐benchtop‐to‐desktop” metabarcoding workflow to investigate the diet of invertebrate‐eaters. We applied our workflow to faecal DNA samples of an invertebrate‐eating fish species. A fragment of the cytochrome c oxidase I (COI) gene was amplified by combining two minibarcoding primer sets to maximize the taxonomic coverage. Amplicons were sequenced by an Illumina MiSeq platform. We developed a filtering approach based on a series of nonarbitrary thresholds established from control samples and from molecular replicates to address the elimination of cross‐contamination, PCR/sequencing errors and mistagging artefacts. This resulted in a conservative and informative metabarcoding data set. We developed a taxonomic assignment procedure that combines different approaches and that allowed the identification of ~75% of invertebrate COI variants to the species level. Moreover, based on the diversity of the variants, we introduced a semiquantitative statistic in our diet study, the minimum number of individuals, which is based on the number of distinct variants in each sample. The metabarcoding approach described in this article may guide future diet studies that aim to produce robust data sets associated with a fine and accurate identification of prey items.     

DNA metabarcoding for diet analysis and biodiversity: A case study using the endangered Australian sea lion (Neophoca cinerea)

The analysis of apex predator diet has the ability to deliver valuable insights into ecosystem health, and the potential impacts a predator might have on commercially relevant species. The Australian sea lion (Neophoca cinerea) is an endemic apex predator and one of the world's most endangered pinnipeds. Given that prey availability is vital to the survival of top predators, this study set out to understand what dietary information DNA metabarcoding could yield from 36 sea lion scats collected across 1,500 km of its distribution in southwest Western Australia. A combination of PCR assays were designed to target a variety of potential sea lion prey, including mammals, fish, crustaceans, cephalopods, and birds. Over 1.2 million metabarcodes identified six classes from three phyla, together representing over 80 taxa. The results confirm that the Australian sea lion is a wide‐ranging opportunistic predator that consumes an array of mainly demersal fauna. Further, the important commercial species Sepioteuthis australis (southern calamari squid) and Panulirus cygnus (western rock lobster) were detected, but were present in <25% of samples. Some of the taxa identified, such as fish, sharks and rays, clarify previous knowledge of sea lion prey, and some, such as eel taxa and two gastropod species, represent new dietary insights. Even with modest sample sizes, a spatial analysis of taxa and operational taxonomic units found within the scat shows significant differences in diet between many of the sample locations and identifies the primary taxa that are driving this variance. This study provides new insights into the diet of this endangered predator and confirms the efficacy of DNA metabarcoding of scat as a noninvasive tool to more broadly define regional biodiversity.     

What you need is what you eat? Prey selection by the bat Myotis daubentonii

Optimal foraging theory predicts that predators are selective when faced with abundant prey, but become less picky when prey gets sparse. Insectivorous bats in temperate regions are faced with the challenge of building up fat reserves vital for hibernation during a period of decreasing arthropod abundances. According to optimal foraging theory, prehibernating bats should adopt a less selective feeding behaviour – yet empirical studies have revealed many apparently generalized species to be composed of specialist individuals. Targeting the diet of the bat Myotis daubentonii, we used a combination of molecular techniques to test for seasonal changes in prey selectivity and individual‐level variation in prey preferences. DNA metabarcoding was used to characterize both the prey contents of bat droppings and the insect community available as prey. To test for dietary differences among M. daubentonii individuals, we used ten microsatellite loci to assign droppings to individual bats. The comparison between consumed and available prey revealed a preference for certain prey items regardless of availability. Nonbiting midges (Chironomidae) remained the most highly consumed prey at all times, despite a significant increase in the availability of black flies (Simuliidae) towards the end of the season. The bats sampled showed no evidence of individual specialization in dietary preferences. Overall, our approach offers little support for optimal foraging theory. Thus, it shows how novel combinations of genetic markers can be used to test general theory, targeting patterns at both the level of prey communities and individual predators.     

Effect of Brood Age on Nestling Diet and Prey Composition in a Hedgerow Specialist Bird,the Barred Warbler Sylvia nisoria

The composition and quality of food provided to nestling birds influence their growth and development and offers key insight into the ecological requirements of birds. One bird species whose feeding ecology is poorly understood is the Barred Warbler (Sylvia nisoria), which utilizes semi-natural shrubby vegetation in agroecosystems. Because Barred Warbler nestlings vary greatly in body mass we hypothesised that diet and prey properties (size, diversity, taxonomic composition, and chitin content and resulting body hardness and digestibility) would differ as the nestlings aged. We quantified the diet based on faecal analysis, sampling faecal sacs from the nestlings pooled into three age classes: 2-3 days old, 4-6 d old, and 7-9 d old. Nestlings were provided a wide diversity of food and a strong relationship existed between food characteristics and nestling age. The youngest nestlings (2-3 d old) had the lowest values of each dietary characteristic (diversity, number and total biomass of prey, and individual prey weight), that were significantly lower than the oldest nestlings (7-9 d old). Nestlings aged 4-6 d exhibited intermediate dietary characteristics. Differences in dietary composition of the six major food types showed marked differences between the individual broods and age categories. Percentages of the number and biomass of soft-bodied prey were highest in the diet of 2-3 d and 4-6 d old nestlings, and decreased with increasing age, whereas the opposite trend was observed in the percentage of intermediately and heavily chitinised prey. Parent Barred Warblers probably preferentially select soft-bodied prey for the youngest nestlings, and satisfy the greater energy demands of the older ones by providing them with a greater variety of prey containing more chitin, as well as plant food. The provisioning of less-readily digestible prey to older nestlings suggests that as the quality of food decreases the quantity increases, implying that the youngest nestlings may be physiologically limited as regards their ability to digest more heavily chitinised prey.     

DNA metabarcoding quantifies the relative biomass of arthropod taxa in songbird diets: Validation with camera‐recorded diets

Ecological research is often hampered by the inability to quantify animal diets. Diet composition can be tracked through DNA metabarcoding of fecal samples, but whether (complex) diets can be quantitatively determined with metabarcoding is still debated and needs validation using free‐living animals. This study validates that DNA metabarcoding of feces can retrieve actual ingested taxa, and most importantly, that read numbers retrieved from sequencing can also be used to quantify the relative biomass of dietary taxa. Validation was done with the hole‐nesting insectivorous Pied Flycatcher whose diet was quantified using camera footage. Size‐adjusted counts of food items delivered to nestlings were used as a proxy for provided biomass of prey orders and families, and subsequently, nestling feces were assessed through DNA metabarcoding. To explore potential effects of digestion, gizzard and lower intestine samples of freshly collected birds were subjected to DNA metabarcoding. For metabarcoding with Cytochrome Oxidase subunit I (COI), we modified published invertebrate COI primers LCO1490 and HCO1777, which reduced host reads to 0.03%, and amplified Arachnida DNA without significant changing the recovery of other arthropod taxa. DNA metabarcoding retrieved all commonly camera‐recorded taxa. Overall, and in each replicate year (N = 3), the relative scaled biomass of prey taxa and COI read numbers correlated at R = .85 (95CI:0.68–0.94) at order level and at R = .75 (CI:0.67–0.82) at family level. Similarity in arthropod community composition between gizzard and intestines suggested limited digestive bias. This DNA metabarcoding validation demonstrates that quantitative analyses of arthropod diet is possible. We discuss the ecological applications for insectivorous birds.     

DNA metabarcoding reveals changes in the contents of carnivorous plants along an elevation gradient

Resource variation along abiotic gradients influences subsequent trophic interactions and these effects can be transmitted through entire food webs. Interactions along abiotic gradients can provide clues as to how organisms will face changing environmental conditions, such as future range shifts. However, it is challenging to find replicated systems to study these effects. Phytotelmata, such as those found in carnivorous plants, are isolated aquatic communities and thus form a good model for the study of replicated food webs. Due to the degraded nature of the prey, molecular techniques provide a useful tool to study these communities. We studied the pitcher plant Sarracenia purpurea L. in allochthonous populations along an elevational gradient in the Alps and Jura. We predicted that invertebrate richness in the contents of the pitcher plants would decrease with increasing elevation, reflecting harsher environmental conditions. Using metabarcoding of the COI gene, we sequenced the invertebrate contents of these pitcher plants. We assigned Molecular Operational Taxonomic Units at ordinal level as well as recovering species‐level data. We found small but significant changes in community composition with elevation. These recovered sequences could belong to invertebrate prey, rotifer inquilines, pollinators and other animals possibly living inside the pitchers. However, we found no directional trend or site‐based differences in MOTU richness with elevational gradient. Use of molecular techniques for dietary or contents analysis is a powerful way to examine numerous degraded samples, although factors such as DNA persistence and the relationship with species presence still have to be completely determined.     

Metabarcoding for the parallel identification of several hundred predators and their prey: Application to bat species diet analysis

Assessing diet variability is of main importance to better understand the biology of bats and design conservation strategies. Although the advent of metabarcoding has facilitated such analyses, this approach does not come without challenges. Biases may occur throughout the whole experiment, from fieldwork to biostatistics, resulting in the detection of false negatives, false positives or low taxonomic resolution. We detail a rigorous metabarcoding approach based on a short COI minibarcode and two‐step PCR protocol enabling the “all at once” taxonomic identification of bats and their arthropod prey for several hundreds of samples. Our study includes faecal pellets collected in France from 357 bats representing 16 species, as well as insect mock communities that mimic bat meals of known composition, negative and positive controls. All samples were analysed using three replicates. We compare the efficiency of DNA extraction methods, and we evaluate the effectiveness of our protocol using identification success, taxonomic resolution, sensitivity and amplification biases. Our parallel identification strategy of predators and prey reduces the risk of mis‐assigning prey to wrong predators and decreases the number of molecular steps. Controls and replicates enable to filter the data and limit the risk of false positives, hence guaranteeing high confidence results for both prey occurrence and bat species identification. We validate 551 COI variants from arthropod including 18 orders, 117 family, 282 genus and 290 species. Our method therefore provides a rapid, resolutive and cost‐effective screening tool for addressing evolutionary ecological issues or developing “chirosurveillance” and conservation strategies.     

Do the evolutionary interactions between moths and bats promote niche partitioning between bats and birds?

Ecological theory suggests that the coexistence of species is promoted by the partitioning of available resources, as in dietary niche partitioning where predators partition prey. Yet, the mechanisms underlying dietary niche partitioning are not always clear. We used fecal DNA metabarcoding to investigate the diets of seven nocturnal insectivorous bird and bat species. Low diet overlap (2%–22%) supported resource partitioning among all species. Differences in diet corresponded with species identity, prey detection method, and foraging behavior of predators. Insects with ultrasonic hearing capabilities were consumed significantly more often by birds than bats, consistent with an evolved avoidance of echolocating strategies. In turn, bats consumed a greater proportion of noneared insects such as spruce budworms. Overall, our results suggest that evolutionary interactions among bats and moths translate to dietary niche partitioning and coexistence among bats and nocturnal birds.     

Temperature affects both the Grinnellian and Eltonian dimensions of ecological niches – A tale of two Arctic wolf spiders

To better understand the consequences of global warming for species and their distribution, we need studies quantifying how environmental change affects communities and interaction networks. Where studies to date have mainly focused on climatic effects on species distribution (the Grinnellian dimension of the niche), recent research has emphasised how the environment shapes ecological interactions among species (the Eltonian dimension).Here, we explore both dimensions in a system consisting of two wolf spider species – Pardosa palustris and Pirata piraticus – and their prey. Drawing on a natural experiment consisting of differential geothermal heating of soil, we describe the effects of temperature on the abundance of each species and on its interactions with its prey (using metabarcoding of gut contents). The two spider species differed substantially in their Grinnellian niche, with a peak in the abundance of P. palustris around 10 °C and in P. piraticus around 22 °C. While P. piraticus consumed more prey taxa on average than did P. palustris, both predators maintained their diet breadth and taxon richness of consumed prey across the temperature gradient. This indicates that effects of temperature on metabolic demands did not alter the dietary specialisation of the two predators. Nevertheless, we did also detect effects of temperature on the Eltonian niche, with significant changes in the prey community consumed by the two spider species across the temperature gradient, and a greater turnover of prey taxa in their diet with increasing soil temperature. Importantly, this suggests that the Eltonian niche of species may be conditional on the environment, and that prey use by generalist predators may thus be modified by climate change.     

From feces to data: A metabarcoding method for analyzing consumed and available prey in a bird‐insect food web

Diets play a key role in understanding trophic interactions. Knowing the actual structure of food webs contributes greatly to our understanding of biodiversity and ecosystem functioning. The research of prey preferences of different predators requires knowledge not only of the prey consumed, but also of what is available. In this study, we applied DNA metabarcoding to analyze the diet of 4 bird species (willow tits Poecile montanus, Siberian tits Poecile cinctus, great tits Parus major and blue tits Cyanistes caeruleus) by using the feces of nestlings. The availability of their assumed prey (Lepidoptera) was determined from feces of larvae (frass) collected from the main foraging habitat, birch (Betula spp.) canopy. We identified 53 prey species from the nestling feces, of which 11 (21%) were also detected from the frass samples (eight lepidopterans). Approximately 80% of identified prey species in the nestling feces represented lepidopterans, which is in line with the earlier studies on the parids' diet. A subsequent laboratory experiment showed a threshold for fecal sample size and the barcoding success, suggesting that the smallest frass samples do not contain enough larval DNA to be detected by high‐throughput sequencing. To summarize, we apply metabarcoding for the first time in a combined approach to identify available prey (through frass) and consumed prey (via nestling feces), expanding the scope and precision for future dietary studies on insectivorous birds.     

Deciphering the diet of a wandering spider (Phoneutria boliviensis; Araneae: Ctenidae) by DNA metabarcoding of gut contents

Arachnids are the most abundant land predators. Despite the importance of their functional roles as predators and the necessity to understand their diet for conservation, the trophic ecology of many arachnid species has not been sufficiently studied. In the case of the wandering spider, Phoneutria boliviensis F. O. Pickard‐Cambridge, 1897, only field and laboratory observational studies on their diet exist. By using a DNA metabarcoding approach, we compared the prey found in the gut content of males and females from three distant Colombian populations of P. boliviensis. By DNA metabarcoding of the cytochrome c oxidase subunit I (COI), we detected and identified 234 prey items (individual captured by the spider) belonging to 96 operational taxonomic units (OTUs), as prey for this wandering predator. Our results broaden the known diet of P. boliviensis with at least 75 prey taxa not previously registered in fieldwork or laboratory experimental trials. These results suggest that P. boliviensis feeds predominantly on invertebrates (Diptera, Lepidoptera, Coleoptera, and Orthoptera) and opportunistically on small squamates. Intersex and interpopulation differences were also observed. Assuming that prey preference does not vary between populations, these differences are likely associated with a higher local prey availability. Finally, we suggest that DNA metabarcoding can be used for evaluating subtle differences in the diet of distinct populations of P. boliviensis, particularly when predation records in the field cannot be established or quantified using direct observation.     

18S rRNA metabarcoding diet analysis of a predatory fish community across seasonal changes in prey availability

Predator–prey relationships are important ecological interactions, affecting biotic community composition and energy flow through a system, and are of interest to ecologists and managers. Morphological diet analysis has been the primary method used to quantify the diets of predators, but emerging molecular techniques using genetic data can provide more accurate estimates of relative diet composition. This study used sequences from the 18S V9 rRNA barcoding region to identify prey items in the gastrointestinal (GI) tracts of predatory fishes. Predator GI samples were taken from the Black River, Cheboygan Co., MI, USA (n = 367 samples, 12 predator species) during periods of high prey availability, including the larval stage of regionally threatened lake sturgeon (Acipenser fulvescens Rafinesque 1817) in late May/early June of 2015 and of relatively lower prey availability in early July of 2015. DNA was extracted and sequenced from 355 samples (96.7%), and prey DNA was identified in 286 of the 355 samples (80.6%). Prey were grouped into 33 ecologically significant taxonomic groups based on the lowest taxonomic level sequences that could be identified using sequences available on GenBank. Changes in the makeup of diet composition, dietary overlap, and predator preference were analyzed comparing the periods of high and low prey abundance. Some predator species exhibited significant seasonal changes in diet composition. Dietary overlap was slightly but significantly higher during the period of high prey abundance; however, there was little change in predator preference. This suggests that change in prey availability was the driving factor in changing predator diet composition and dietary overlap. This study demonstrates the utility of molecular diet analysis and how temporal variability in community composition adds complexity to predator–prey interactions.     

Disparities in second‐generation DNA metabarcoding results exposed with accessible and repeatable workflows

Different second‐generation sequencing technologies may have taxon‐specific biases when DNA metabarcoding prey in predator faeces. Our major objective was to examine differences in prey recovery from bat guano across two different sequencing workflows using the same faecal DNA extracts. We compared results between the Ion Torrent PGM and the Illumina MiSeq with similar library preparations and the same analysis pipeline. We focus on repeatability and provide an R Notebook in an effort towards transparency for future methodological improvements. Full documentation of each step enhances the accessibility of our analysis pipeline. We tagged DNA from insectivorous bat faecal samples, targeted the arthropod cytochrome c oxidase I minibarcode region and sequenced the product on both second‐generation sequencing platforms. We developed an analysis pipeline with a high operational taxonomic unit (OTU) clustering threshold (i.e., ≥98.5%) followed by copy number filtering to avoid merging rare but genetically similar prey into the same OTUs. With this workflow, we detected 297 unique prey taxa, of which 74% were identified at the species level. Of these, 104 (35%) prey OTUs were detected by both platforms, 176 (59%) OTUs were detected by the Illumina MiSeq system only, and 17 (6%) OTUs were detected using the Ion Torrent system only. Costs were similar between platforms but the Illumina MiSeq recovered six times more reads and four additional insect orders than did Ion Torrent. The considerations we outline are particularly important for long‐term ecological monitoring; a more standardized approach will facilitate comparisons between studies and allow faster recognition of changes within ecological communities.     

Temperature shapes opposing latitudinal gradients of plant taxonomic and phylogenetic β diversity

Latitudinal and elevational richness gradients have received much attention from ecologists but there is little consensus on underlying causes. One possible proximate cause is increased levels of species turnover, or β diversity, in the tropics compared to temperate regions. Here, we leverage a large botanical dataset to map taxonomic and phylogenetic β diversity, as mean turnover between neighboring 100 × 100 km cells, across the Americas and determine key climatic drivers. We find taxonomic and tip‐weighted phylogenetic β diversity is higher in the tropics, but that basal‐weighted phylogenetic β diversity is highest in temperate regions. Supporting Janzen's ‘mountain passes’ hypothesis, tropical mountainous regions had higher β diversity than temperate regions for taxonomic and tip‐weighted metrics. The strongest climatic predictors of turnover were average temperature and temperature seasonality. Taken together, these results suggest β diversity is coupled to latitudinal richness gradients and that temperature is a major driver of plant community composition and change.     

Stable isotope analysis as a tool to monitor dietary trends in little penguins Eudyptula minor

Long‐term dietary monitoring of seabirds can be used to relate population fluctuations to at‐sea events. Stomach flushing is a conventional dietary monitoring technique, but has a number of disadvantages. Stable isotope analysis (SIA) is a less invasive method that provides unbiased dietary information over a longer period. We evaluated stable isotope analysis as a potential tool for monitoring long‐term little penguin Eudyptula minor diet. We determined diet composition during the chick feeding stage using stomach flushing and SIA at three separate colonies, using spatial variation in diet as a surrogate for potential temporal variation. Bayesian isotopic mixing models were generated for blood and feathers to evaluate their ability to discriminate broad‐scale (fish, squid, crustaceans) and fine‐scale (individual prey species) diet composition. Differences in stable carbon and nitrogen isotope ratios were found between colonies: broad‐scale isotopic mixing models predicted different proportional contributions of broad taxa (fish, cephalopod, crustacean) to diet than was indicated by stomach samples, reflecting the bias incurred by one‐off stomach contents analysis. Fine‐scale isotopic mixing models predicted proportional contributions of prey items with less certainty. Blood isotopic mixing models had narrower confidence intervals than models for feathers, but trends in δ¹⁵N for feathers mirrored those for blood. Our results suggest that relying on stomach contents analysis to detect shifts in prey consumption in little penguins could be very misleading, resulting in a less‐than‐complete idea of total prey consumption. SIA of little penguin tissues could be used to monitor dietary shifts across dissimilar taxa that may affect population numbers, but would fail to detect shifts between fish species.     

Long‐term relationship between diet breadth and breeding success in a declining population of Egyptian Vultures Neophron percnopterus

Between 2000 and 2009 we studied the diet and breeding success of Egyptian Vultures Neophron percnopterus in southern Spain. Wild species accounted for 74.9% of prey items (n = 1071) with a predominance of mammals (62.3%), followed by birds (20.8%) and reptiles (13.1%). Spatially, the diet was highly varied and not restricted to carcasses of livestock; wild Rabbits Oryctolagus cuniculus accounted for 54% of the overall remains. The spatial variability may reflect regional and local disparity in the availability of main prey. The temporal relationship between variation in trophic diversity and Vulture nesting productivity (both values showing a long‐term decrease) might suggest a causal link between variation in diet and reproductive output. We hypothesize that high turnover rates could explain productivity variation as a consequence of the recruitment of less experienced individuals to the breeding population. This could in turn generate covariation between diet and reproductive output.     

Seasonal and ontological variation in diet and age‐related differences in prey choice,by an insectivorous songbird

The diet of an individual animal is subject to change over time, both in response to short‐term food fluctuations and over longer time scales as an individual ages and meets different challenges over its life cycle. A metabarcoding approach was used to elucidate the diet of different life stages of a migratory songbird, the Eurasian reed warbler (Acrocephalus scirpaceus) over the 2017 summer breeding season in Somerset, the United Kingdom. The feces of adult, juvenile, and nestling warblers were screened for invertebrate DNA, enabling the identification of prey species. Dietary analysis was coupled with monitoring of Diptera in the field using yellow sticky traps. Seasonal changes in warbler diet were subtle, whereas age class had a greater influence on overall diet composition. Age classes showed high dietary overlap, but significant dietary differences were mediated through the selection of prey; (i) from different taxonomic groups, (ii) with different habitat origins (aquatic vs. terrestrial), and (iii) of different average approximate sizes. Our results highlight the value of metabarcoding data for enhancing ecological studies of insectivores in dynamic environments.