Abstract1. Metabarcoding of environmental samples has many challenges and limitations that require carefully considered laboratory and analysis workflows to ensure reliable results. We explore how decisions regarding study design, laboratory set-up, and bioinformatic processing affect the final results, and provide guidelines for reliable study of environmental samples.2. We evaluate the performance of four primer sets targeting COI and 16S regions characterizing arthropod diversity in bat faecal samples, and investigate how metabarcoding results are affected by parameters including: (1) number of PCR replicates per sample, (2) sequencing depth, (3) PCR replicate processing strategy (i.e. either additively, by combining the sequences obtained from the PCR replicates, or restrictively, by only retaining sequences that occur in multiple PCR replicates for each sample), (4) minimum copy number for sequences to be retained, (5) chimera removal, and (6) similarity thresholds for Operational Taxonomic Unit (OTU) clustering. Lastly, we measure within-and between-taxa dissimilarities when using sequences from public databases to determine the most appropriate thresholds for OTU clustering and taxonomy assignment.3. Our results show that the use of multiple primer sets reduces taxonomic biases and increases taxonomic coverage. Taxonomic profiles resulting from each primer set are principally affected by how many PCR replicates are carried out per sample and how sequences are filtered across them, the sequence copy number threshold and the OTU clustering threshold. We also report considerable diversity differences between PCR replicates from each sample. Sequencing depth increases the dissimilarity between PCR replicates unless the bioinformatic strategies to remove allegedly artefactual sequences are adjusted according to the number of analysed sequences. Finally, we show that the appropriate identity thresholds for OTU clustering and taxonomy assignment differ between markers.4. Metabarcoding of complex environmental samples ideally requires (1) investigation of whether more than one primer sets targeting the same taxonomic group is needed to offset primer biases, (2) more than one PCR replicate per sample, (3) bioinformatic processing of sequences that balance diversity detection with removal of artefactual sequences, and (4) empirical selection of OTU clustering and taxonomy assignment thresholds tailored to each marker and the obtained taxa.

No abstract

The application of high‐throughput sequencing‐based approaches to DNA extracted from environmental samples such as gut contents and faeces has become a popular tool for studying dietary habits of animals. Due to the high resolution and prey detection capacity they provide, both metabarcoding and shotgun sequencing are increasingly used to address ecological questions grounded in dietary relationships. Despite their great promise in this context, recent research has unveiled how a wealth of biological (related to the study system) and technical (related to the methodology) factors can distort the signal of taxonomic composition and diversity. Here, we review these studies in the light of high‐throughput sequencing‐based assessment of trophic interactions. We address how the study design can account for distortion factors, and how acknowledging limitations and biases inherent to sequencing‐based diet analyses are essential for obtaining reliable results, thus drawing appropriate conclusions. Furthermore, we suggest strategies to minimize the effect of distortion factors, measures to increase reproducibility, replicability and comparability of studies, and options to scale up DNA sequencing‐based diet analyses. In doing so, we aim to aid end‐users in designing reliable diet studies by informing them about the complexity and limitations of DNA sequencing‐based diet analyses, and encourage researchers to create and improve tools that will eventually drive this field to its maturity.

The interaction between agricultural production and wildlife can shape, and even condition, the functioning of both systems. In this study, we i) explored the degree to which a widespread European bat, namely the common bent-wing bat Miniopterus schreibersii, consumes crop-damaging insects at a continental scale, and ii) tested whether its dietary niche is shaped by the extension and type of agricultural fields. We employed a dual-primer DNA metabarcoding approach to characterize arthropod 16S and COI DNA sequences within bat faecal pellets collected across 16 Southern European localities, to first characterize the bat species' dietary niche, second measure the incidence of agricultural pests across their ranges and third assess whether geographical dietary variation responds to climatic, landscape diversity, agriculture type and vegetation productivity factors. We detected 12 arthropod orders, among which lepidopterans were predominant. We identified >200 species, 44 of which are known to cause agricultural damage. Pest species were detected at all but one sampling site and in 94% of the analysed samples. Furthermore, the dietary diversity of M. schreibersii exhibited a negative linear relation with the area of intensive agricultural fields, thus suggesting crops restrict the dietary niche of bats to prey taxa associated with agricultural production within their foraging range. Overall, our results imply that M. schreibersii might be a valuable asset for biological pest suppression in a variety of agricultural productions and highlight the dynamic interplay between wildlife and agricultural systems.

Summary From ontogenesis to homeostasis, the phenotypes of complex organisms are shaped by the bidirectional interactions between the host organisms and their associated microbiota. Current technology can reveal many such interactions by combining multi-omic data from both hosts and microbes. However, exploring the full extent of these interactions requires careful consideration of study design for the efficient generation and optimal integration of data derived from (meta)genomics, (meta)transcriptomics, (meta)proteomics, and (meta)metabolomics. In this perspective, we introduce the holo-omic approach that incorporates multi-omic data from both host and microbiota domains to untangle the interplay between the two. We revisit the recent literature on biomolecular host-microbe interactions and discuss the implementation and current limitations of the holo-omic approach. We anticipate that the application of this approach can contribute to opening new research avenues and discoveries in biomedicine, biotechnology, agricultural and aquacultural sciences, nature conservation, as well as basic ecological and evolutionary research.

Molecular analysis of diet overcomes the considerable limitations of traditional techniques for identifying prey remains in bat faeces. We collected faeces from individual Mountain Long-eared Bats Plecotus macrobullaris trapped using mist nets during the summers of 2009 and 2010 in the Pyrenees. We analysed their diet using DNA mini-barcodes to identify prey species. In addition, we inferred some basic features of the bat's foraging ecology that had not yet been addressed. P. macrobullaris fed almost exclusively on moths (97.8%). As prey we detected one dipteran genus (Tipulidae) and 29 moth taxa: 28 were identified at species level (23 Noctuidae, 1 Crambidae, 1 Geometridae, 1 Pyralidae, 1 Sphingidae, 1 Tortricidae), and one at genus level (Rhyacia sp., Noctuidae). Known ecological information about the prey species allowed us to determine that bats had foraged at elevations between 1,500 and 2,500 m amsl (above mean sea level), mostly in subalpine meadows, followed by other open habitats such as orophilous grasslands and alpine meadows. No forest prey species were identified in the diet. As 96.4% of identified prey species were tympanate moths and no evidence of gleaning behaviour was revealed, we suggest P. macrobullaris probably forages by aerial hawking using faint echolocation pulses to avoid detection by hearing moths. As we could identify 87.8% of the analysed sequences (64.1% of the MOTUs, Molecular Operational Taxonomic Units) at species level, we conclude that DNA mini-barcodes are a very useful tool to analyse the diet of moth-specialist bats.

Inferences of the interactions between species' ecological niches and spatial distribution have been historically based on simple metrics such as low-resolution dietary breadth and range size, which might have impeded the identification of meaningful links between niche features and spatial patterns. We analysed the relationship between dietary niche breadth and spatial distribution features of European bats, by combining continent-wide DNA metabarcoding of faecal samples with species distribution modelling. Our results show that while range size is not correlated with dietary features of bats, the homogeneity of the spatial distribution of species exhibits a strong correlation with dietary breadth. We also found that dietary breadth is correlated with bats' hunting flexibility. However, these two patterns only stand when the phylogenetic relations between prey are accounted for when measuring dietary breadth. Our results suggest that the capacity to exploit different prey types enables species to thrive in more distinct environments and therefore exhibit more homogeneous distributions within their ranges.

Formerly thought to be a strictly insectivorous trawling bat, recent studies have shown that Myotis capaccinii also preys on fish. To determine whether differences exist in bat flight behaviour, prey handling and echolocation characteristics when catching fish and insects of different size, we conducted a field experiment focused on the last stage of prey capture. We used synchronized video and ultrasound recordings to measure several flight and dip features as well as echolocation characteristics, focusing on terminal buzz phase I, characterized by a call rate exceeding 100 Hz, and buzz phase II, characterized by a drop in the fundamental well below 20 kHz and a repetition rate exceeding 150 Hz. When capturing insects, bats used both parts of the terminal phase to the same extent, and performed short and superficial drags on the water surface. In contrast, when preying on fish, buzz I was longer and buzz II shorter, and the bats made longer and deeper dips. These variations suggest that lengthening buzz I and shortening buzz II when fishing is beneficial, probably because buzz I gives better discrimination ability and the broader sonar beam provided by buzz II is useless when no evasive flight of the prey is expected. Additionally, bats continued emitting calls beyond the theoretical signal-overlap zone, suggesting that they might obtain information even when they have surpassed that threshold, at least initially. This study shows that M. capaccinii can regulate the temporal components of its feeding buzzes and modify prey capture technique according to the target.