The application of DNA metabarcoding to dietary analysis of trophic generalists requires using multiple markers in order to overcome problems of primer specificity and bias. However, limited attention has been given to the integration of information from multiple markers, particularly when they partly overlap in the taxa amplified, and vary in taxonomic resolution and biases. Here, we test the use of a mix of universal and specific markers, provide criteria to integrate multi‐marker metabarcoding data and a python script to implement such criteria and produce a single list of taxa ingested per sample. We then compare the results of dietary analysis based on morphological methods, single markers, and the proposed combination of multiple markers. The study was based on the analysis of 115 faeces from a small passerine, the Black Wheatears (Oenanthe leucura). Morphological analysis detected far fewer plant taxa (12) than either a universal 18S marker (57) or the plant trnL marker (124). This may partly reflect the detection of secondary ingestion by molecular methods. Morphological identification also detected far fewer taxa (23) than when using 18S (91) or the arthropod markers IN16STK (244) and ZBJ (231), though each method missed or underestimated some prey items. Integration of multi‐marker data provided far more detailed dietary information than any single marker and estimated higher frequencies of occurrence of all taxa. Overall, our results show the value of integrating data from multiple, taxonomically overlapping markers in an example dietary data set.

DNA metabarcoding is increasingly used in dietary studies to estimate diversity, composition and frequency of occurrence of prey items. However, few studies have assessed how technical and biological replication affect the accuracy of diet estimates. This study addresses these issues using the European free‐tailed bat Tadarida teniotis, involving high‐throughput sequencing of a small fragment of the COI gene in 15 separate faecal pellets and a 15‐pellet pool per each of 20 bats. We investigated how diet descriptors were affected by variability among (a) individuals, (b) pellets of each individual and (c) PCRs of each pellet. In addition, we investigated the impact of (d) analysing separate pellets vs. pellet pools. We found that diet diversity estimates increased steadily with the number of pellets analysed per individual, with seven pellets required to detect ~80% of prey species. Most variation in diet composition was associated with differences among individual bats, followed by pellets per individual and PCRs per pellet. The accuracy of frequency of occurrence estimates increased with the number of pellets analysed per bat, with the highest error rates recorded for prey consumed infrequently by many individuals. Pools provided poor estimates of diet diversity and frequency of occurrence, which were comparable to analysing a single pellet per individual, and consistently missed the less common prey items. Overall, our results stress that maximizing biological replication is critical in dietary metabarcoding studies and emphasize that analysing several samples per individual rather than pooled samples produce more accurate results.

In multifunctional landscapes, diverse communities of flying vertebrate predators provide vital services of insect pest control. In such landscapes, conservation biocontrol should benefit service‐providing species to enhance the flow, stability and resilience of pest control services supporting the production of food and fiber. However, this would require identifying key service providers, which may be challenging when multiple predators interact with multiple pests. Here we provide a framework to identify the functional role of individual species to pest control in multifunctional landscapes. First, we used DNA metabarcoding to provide detailed data on pest species predation by diverse predator communities. Then, these data were fed into an extensive network analysis, in which information relevant for conservation biocontrol is gained from parameters describing network structure (e.g., modularity) and species roles in such network (e.g., centrality, specialization). We applied our framework to a Mediterranean landscape, where 19 bat species were found to feed on 132 insect pest species. Metabarcoding data revealed potentially important bats that consumed insect pest species in high frequency and/or diversity. Network analysis showed a modular structure, indicating sets of bat species that are required to regulate specific sets of insect pests. A few generalist bats had particularly important roles, either at network or module levels. Extinction simulations highlighted six bats, including species of conservation concern, which were sufficient to ensure that over three‐quarters of the pest species had at least one bat predator. Combining DNA metabarcoding and ecological network analysis provides a valuable framework to identify individual species within diverse predator communities that might have a disproportionate contribution to pest control services in multifunctional landscapes. These species can be regarded as candidate targets for conservation biocontrol, although additional information is needed to evaluate their actual effectiveness in pest regulation.

In bats, sexual segregation has been described in relation to differential use of roosting and foraging habitats. It is possible that variation may also exist between genders in the use of different prey types. However, until recently this idea was difficult to test owing to poorly resolved taxonomy of dietary studies. Here, we use high-throughput sequencing to describe gender-related variation in diet composition of the European free-tailed bat ( Tadarida teniotis ), while controlling for effects of age and season. We analysed guano pellets collected from 143 individuals mist-netted from April to October 2012 and 2013, in northeast Portugal. Moths (Lepidoptera; mainly Noctuidae and Geometridae) were by far the most frequently recorded prey, occurring in nearly all samples and accounting for 96 out of 115 prey taxa. There were significant dietary differences between males and females, irrespective of age and season. Compared to males, females tended to consume larger moths and more moths of migratory behaviour (e.g. Autographa gamma ). Our study provides the first example of gender-related dietary variation in bats, illustrating the value of novel molecular tools for revealing intraspecific variation in food resource use in bats and other insectivores.