A species-specific digital PCR assay for the endangered blue racer (Coluber constrictor foxii) in Canada.

The only population of the endangered blue racer (Coluber constrictor foxii) in Canada occurs on Pelee Island, Ontario. The species is threatened by multiple factors, including habitat degradation and loss, road mortality, persecution, and potentially predation. We designed and evaluated the performance of an environmental DNA droplet digital PCR assay that can be used for multiple facets of conservation of this species. We tested the assay in silico and in vitro using DNA of blue racers and co-occurring snake species and estimated the LOD and LOQ using synthetic DNA. As wild turkey predation has been suggested to negatively affect racers, we tested the assay on eight wild turkey faecal samples. Our assay is specific, can detect the target species at very low levels of concentration (0.002 copies/µL), and can accurately quantify copy numbers ≥ 0.26 copies/µL. We detected no racer DNA in any wild turkey faecal sample. More faecal samples collected at strategic locations during snake peak activity on Pelee Island would enable a more thorough assessment of the possibility of turkey predation. Our assay should be effective for other environmental samples and can be used for investigating other factors negatively affecting blue racers, for example, helping to quantify blue racer habitat suitability and site occupancy.

Assessing the breeding phenology of a threatened frog species using eDNA and automatic acoustic monitoring.

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.

In silico and empirical evaluation of twelve metabarcoding primer sets for insectivorous diet analyses.

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.

Integrating population genetics to define conservation units from the core to the edge of Rhinolophus ferrumequinum western range.

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.

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.