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.

During the most recent decade, environmental DNA metabarcoding approaches have been both developed and improved to minimize the biological and technical biases in these protocols. However, challenges remain, notably those relating to primer design. In the current study, we comprehensively assessed the performance of ten COI and two 16S primer pairs for eDNA metabarcoding, including novel and previously published primers. We used a combined approach of in silico, in vivo‐mock community (33 arthropod taxa from 16 orders), and guano‐based analyses to identify primer sets that would maximize arthropod detection and taxonomic identification, successfully identify the predator (bat) species, and minimize the time and financial costs of the experiment. We focused on two insectivorous bat species that live together in mixed colonies: the greater horseshoe bat (Rhinolophus ferrumequinum) and Geoffroy's bat (Myotis emarginatus). We found that primer degeneracy is the main factor that influences arthropod detection in silico and mock community analyses, while amplicon length is critical for the detection of arthropods from degraded DNA samples. Our guano‐based results highlight the importance of detecting and identifying both predator and prey, as guano samples can be contaminated by other insectivorous species. Moreover, we demonstrate that amplifying bat DNA does not reduce the primers' capacity to detect arthropods. We therefore recommend the simultaneous identification of predator and prey. Finally, our results suggest that up to one‐third of prey occurrences may be unreliable and are probably not of primary interest in diet studies, which may decrease the relevance of combining several primer sets instead of using a single efficient one. In conclusion, this study provides a pragmatic framework for eDNA primer selection with respect to scientific and methodological constraints.

Dietary plasticity can be a determining factor allowing species to cope with environmental changes. Consequently, it is an important issue to consider in conservation biology. Despite this, it remains rarely addressed in the literature, potentially due to methodologies which were until recently rather limited. The advent of molecular approaches now makes it possible to get a precise picture of diet and its plasticity, even for endangered and elusive species. Here, we focused on the greater horseshoe bat (Rhinolophus ferrumequinum) in Western France, where this insectivorous species has been classified as “Vulnerable” on the Regional Red List in 2016. We applied an eDNA metabarcoding approach to 1986 fecal samples collected in six maternity colonies on three sampling dates. We described the diet and investigated whether the landscape surrounding colonies and the different phases of the maternity cycle influenced the diversity and the composition of this diet. We showed that R. ferrumequinum feed on a much more diverse prey spectrum than expected from previous studies, highlighting how eDNA metabarcoding can improve our knowledge on the dietary habits of elusive species. Our approach also revealed that the diet of R. ferrumequinum seems to be composed of two distinct features: the core diet consisting of a few preferred taxa shared by all the colonies (25% of the occurrences) and the secondary diet consisting of numerous rare prey taxa that were highly different between colonies and sampling dates (75% of the occurrences). Constraints associated with the greater horseshoe bat life cycle, as well as insect phenology and landscape features, strongly influenced the diversity and composition of both the core diet and the diet as a whole. Further research should now explore the relationships between R. ferrumequinum dietary plasticity and fitness, to better assess the impact of core prey decline on R. ferrumequinum population viability.

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 preys 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 preys reduces the risk of mis-assigning preys 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.

The greater horseshoe bat (Rhinolophus ferrumequinum) is among the most widespread bat species in Europe but it has experienced severe declines, especially in Northern Europe. This species is listed Near Threatened in the European IUCN Red List of Threatened Animals, and it is considered to be highly sensitive to human activities and particularly to habitat fragmentation. Therefore, understanding the population boundaries and demographic history of populations of this species is of primary importance to assess relevant conservation strategies. In this study, we used 17 microsatellite markers to assess the genetic diversity, the genetic structure, and the demographic history of R. ferrumequinum colonies in the western part of its distribution. We identified one large population showing high levels of genetic diversity and large population size. Lower estimates were found in England and northern France. Analyses of clustering and isolation by distance suggested that the Channel and the Mediterranean seas could impede R. ferrumequinum gene flow. These results provide important information to improve the delineation of R. ferrumequinum management units. We suggest that a large management unit corresponding to the population ranging from Spanish Basque Country to northern France must be considered. Particular attention should be given to mating territories as they seem to play a key role in maintaining high levels of genetic mixing between colonies. Smaller management units corresponding to English and northern France colonies must also be implemented. These insular or peripheral colonies could be at higher risk of extinction in the near future.

This last decade, environmental DNA metabarcoding approaches have been developed and improved to minimize biological and technical biases; some challenges, however, remain, as the design of primers. Here we have performed a comprehensive assessment of ten COI and two 16S primer sets. We have combined in silico, in vivo-mock community of 33 arthropod taxa from 16 orders and guano analyses to identify primer sets that should maximize arthropod detection and taxonomic identification, whilst identifying bat species and minimizing labour time and cost. We have focused on two insectivorous bat species living in mixed-colonies, the greater horseshoe bat (Rhinolophus ferrumequinum) and Geoffroy’s bat (Myotis emarginatus). We have found that the level of primer degeneracy is the main factor influencing arthropod detection for in silico and mock community analyses, while the amplicon length is critical for the detection of arthropods from degraded DNA samples. Our results confirm the importance of performing predator detection and taxonomic identification, simultaneously with arthropod sequencing, as faeces samples can be contaminated by different insectivorous species. Moreover, amplifying bat DNA does not affect the primers’ capacity to detect arthropods. We therefore recommend the systematic simultaneous identification of predator and prey. Finally, we evidenced that one third of the prey occurrences are unreliable and probably not of primary interest in diet studies, which might decrease the relevance of combining several primer sets instead of using one efficient primer set. In conclusion, this study provides general criteria enabling the selection of primers whilst considering different scientific and methodological constraints.

Background Climate change has driven shifts in breeding phenology of many amphibians, causing phenological mismatches (e.g., predator-prey interactions), and potentially population declines. Collecting data with high spatiotemporal sensitivity on hibernation emergence and breeding times can inform conservation best practices. However, monitoring the phenology of amphibians can be challenging because of their cryptic nature over much of their life cycle. Moreover, most salamanders and caecilians do not produce conspicuous breeding calls like frogs and toads do, presenting additional monitoring challenges. Methods In this study, we designed and evaluated the performance of an environmental DNA (eDNA) droplet digital PCR (ddPCR) assay as a non-invasive tool to assess the breeding phenology of a Western Chorus Frog population (Pseudacris maculata mitotype) in Eastern Ontario and compared eDNA detection patterns to hourly automatic acoustic monitoring. For two eDNA samples with strong PCR inhibition, we tested three methods to diminish the effect of inhibitors: diluting eDNA samples, adding bovine serum albumin to PCR reactions, and purifying eDNA using a commercial clean-up kit. Results We recorded the first male calling when the focal marsh was still largely frozen. Chorus frog eDNA was detected on April 6th, 6 days after acoustic monitoring revealed this first calling male, but only 2 days after males attained higher chorus activity. eDNA signals were detected at more sampling locales within the marsh and eDNA concentrations increased as more males participated in the chorus, suggesting that eDNA may be a reasonable proxy for calling assemblage size. Internal positive control revealed strong inhibition in some samples, limiting detection probability and quantification accuracy in ddPCR. We found diluting samples was the most effective in reducing inhibition and improving eDNA quantification. Conclusions Altogether, our results showed that eDNA ddPCR signals lagged behind male chorusing by a few days; thus, acoustic monitoring is preferable if the desire is to document the onset of male chorusing. However, eDNA may be an effective, non-invasive monitoring tool for amphibians that do not call and may provide a useful complement to automated acoustic recording. We found inhibition patterns were heterogeneous across time and space and we demonstrate that an internal positive control should always be included to assess inhibition for eDNA ddPCR signal interpretations.

Dietary plasticity is an important issue for conservation biology as it may be essential for species to cope with environmental changes. However, this process still remains scarcely addressed in the literature, potentially because diet studies have long been constrained by methodological limits. The advent of molecular approaches now makes it possible to get a precise picture of diet and its plasticity, even for endangered and elusive species. Here we focused on the greater horseshoe bat (Rhinolophus ferrumequinum) in Western France, where this insectivorous species has been classified as ‘Vulnerable’ on the Regional Red List in 2016. We applied an eDNA metabarcoding approach on 1986 fecal samples collected in six maternity colonies at three sampling dates. We described its diet and investigated whether the landscape surrounding colonies and the different phases of the maternity cycle influenced the diversity and the composition of this diet. We showed that R. ferrumequinum feed on a highly more diverse spectrum of prey than expected from previous studies, therefore highlighting how eDNA metabarcoding can help improving diet knowledge of a flying elusive endangered species. Our approach also revealed that R. ferrumequinum diet is composed of two distinct features: the core diet consisting in a few preferred taxa shared by all the colonies (25% of the occurrences) and the secondary diet consisting in numerous rare prey that were highly different between colonies and sampling dates (75% of the occurrences). Energetic needs and constraints associated with the greater horseshoe bat life-cycle, as well as insect phenology and landscape features, strongly influenced the diversity and composition of both the core and whole diets. Further research should now explore the relationships between R. ferrumequinum dietary plasticity and fitness, to better assess the impact of core prey decline on R. ferrumequinum populations viability.