Surveying Europe's Only Cave-Dwelling Chordate Species (Proteus anguinus) Using Environmental DNA

Enhancing amphibian biomonitoring through eDNA metabarcoding

Surveying biodiversity has taken a quantum leap with environmental DNA (eDNA) metabarcoding, an immensely powerful approach lauded for its efficiency, sensitivity, and non-invasiveness. This approach emerges as a game-changer for the elusive realm of endangered and rare species-think nocturnal, environmentally elusive amphibians. Here, we have established a framework for constructing a reliable metabarcoding pipeline for amphibians, covering primer design, performance evaluation, laboratory validation, and field validation processes. The Am250 primer, located on the mitochondrial 16S gene, was optimal for the eDNA monitoring of amphibians, which demonstrated higher taxonomic resolution, smaller species amplification bias, and more extraordinary detection ability compared to the other primers tested. Am250 primer exhibit an 83.8% species amplification rate and 75.4% accurate species identification rate for Chinese amphibians in the in silico PCR and successfully amplified all tested species of the standard samples in the in vitro assay. Furthermore, the field-based mesocosm experiment showed that DNA can still be detected by metabarcoding even days to weeks after organisms have been removed from the mesocosm. Moreover, field mesocosm findings indicate that eDNA metabarcoding primers exhibit different read abundances, which can affect the relative biomass of species. Thus, appropriate primers should be screened and evaluated by three experimental approaches: in silico PCR simulation, target DNA amplification, and mesocosm eDNA validation. The selection of a single primer set or multiple primers' combination should be based on the monitoring groups to improve the species detection rate and the credibility of results.

Using eDNA to survey amphibians: Methods, applications, and challenges

In recent years, environmental DNA (eDNA) has received attention from biologists due to its sensitivity, convenience, labor and material efficiency, and lack of damage to organisms. The extensive application of eDNA has opened avenues for the monitoring and biodiversity assessment of amphibians, which are frequently small and difficult to observe in the field, in areas such as biodiversity survey assessment and detection of specific, rare and threatened, or alien invasive species. However, the accuracy of eDNA can be influenced by factors such as ambient temperature, pH, and false positives or false negatives, which makes eDNA an adjunctive tool rather than a replacement for traditional surveys. This review provides a concise overview of the eDNA method and its workflow, summarizes the differences between applying eDNA for detecting amphibians and other organisms, reviews the research progress in eDNA technology for amphibian monitoring, identifies factors influencing detection efficiency, and discusses the challenges and prospects of eDNA. It aims to serve as a reference for future research on the application of eDNA in amphibian detection.

Current status and topical issues on the use of eDNA-based targeted detection of rare animal species

Animal detection through DNA present in environmental samples (eDNA) is a valuable tool for detecting rare species, that are difficult to observe and monitor. eDNA-based tools are underpinned by molecular evolutionary principles, key to devising tools to efficiently single out a targeted species from an environmental sample. Here, we present a comprehensive review of the use of eDNA-based methods for the detection of targeted animal species, such as rare, endangered, or invasive species, through the analysis of 549 publications (2008-2022). Aquatic ecosystems have been the most surveyed, in particular, freshwaters (74 %), and to a less extent marine (14 %) and terrestrial systems (10 %). Vertebrates, in particular, fish (38 %), and endangered species, have been the focus of most of these studies, and Cytb and COI are the most employed markers. Among invertebrates, assays have been mainly designed for Mollusca and Crustacea species (21 %), in particular, to target invasive species, and COI the most employed marker. Targeted molecular approaches, in particular qPCR, have been the most adopted (75 %), while eDNA metabarcoding has been rarely used to target single or few species (approx. 6 %). However, less attention has been given in these studies to the effects of environmental factors on the amount of shed DNA, the differential amount of shed DNA among species, or the sensitivity of the markers developed, which may impact the design of the assays, particularly to warrant the required detection level and avoid false negatives and positives. The accuracy of the assays will also depend on the availability of genetic data and vouchered tissue or DNA samples from closely related species to assess both marker and primers' specificity. In addition, eDNA-based assays developed for a particular species may have to be refined for use in a new geographic area taking into account site-specific populations, as well as any intraspecific variation.

Development and validation of an environmental DNA assay to detect federally threatened groundwater salamanders in central Texas

The molecular detection of DNA fragments that are shed into the environment (eDNA) has become an increasingly applied tool used to inventory biological communities and to perform targeted species surveys. This method is particularly useful in habitats where it is difficult or not practical to visually detect or trap the target organisms. Central Texas Eurycea salamanders inhabit both surface and subterranean aquatic environments. Subterranean surveys are challenging or infeasible, and the detection of salamander eDNA in water samples is an appealing survey technique for these situations. Here, we develop and validate an eDNA assay using quantitative PCR for E. chisholmensis, E. naufragia, and E. tonkawae. These three species are federally threatened and constitute the Septentriomolge clade that occurs in the northern segment of the Edwards Aquifer. First, we validated the specificity of the assay in silico and with DNA extracted from tissue samples of both target Septentriomolge and non-target amphibians that overlap in distribution. Then, we evaluated the sensitivity of the assay in two controls, one with salamander-positive water and one at field sites known to be occupied by Septentriomolge. For the salamander-positive control, the estimated probability of eDNA occurrence (ψ) was 0.981 (SE = 0.019), and the estimated probability of detecting eDNA in a qPCR replicate (p) was 0.981 (SE = 0.011). For the field control, the estimated probability of eDNA occurring at a site (ψ) was 0.938 (95% CRI: 0.714-0.998). The estimated probability of collecting eDNA in a water sample (θ) was positively correlated with salamander relative density and ranged from 0.371 (95% CRI: 0.201-0.561) to 0.999 (95% CRI: 0.850- > 0.999) among sampled sites. Therefore, sites with low salamander density require more water samples for eDNA evaluation, and we determined that our site with the lowest estimated θ would require seven water samples for the cumulative collection probability to exceed 0.95. The estimated probability of detecting eDNA in a qPCR replicate (p) was 0.882 (95% CRI: 0.807-0.936), and our assay required two qPCR replicates for the cumulative detection probability to exceed 0.95. In complementary visual encounter surveys, the estimated probability of salamanders occurring at a known-occupied site was 0.905 (SE = 0.096), and the estimated probability of detecting salamanders in a visual encounter survey was 0.925 (SE = 0.052). We additionally discuss future research needed to refine this method and understand its limitations before practical application and incorporation into formal survey protocols for these taxa.

A validated protocol for eDNA-based monitoring of within-species genetic diversity in a pond-breeding amphibian

In light of the dramatic decline in amphibian biodiversity, new cost-efficient tools to rapidly monitor species abundance and population genetic diversity in space and time are urgently needed. It has been amply demonstrated that the use of environmental DNA (eDNA) for single-species detection and characterization of community composition can increase the precision of amphibian monitoring compared to traditional (observational) approaches. However, it has been suggested that the efficiency and accuracy of the eDNA approach could be further improved by more timely sampling; in addition, the quality of genetic diversity data derived from the same DNA has been confirmed in other vertebrate taxa, but not amphibians. Given the availability of previous tissue-based genetic data, here we use the common frog Rana temporaria Linnaeus, 1758 as our target species and an improved eDNA protocol to: (i) investigate differences in species detection between three developmental stages in various freshwater environments; and (ii) study the diversity of mitochondrial DNA (mtDNA) haplotypes detected in eDNA (water) samples, by amplifying a specific fragment of the COI gene (331 base pairs, bp) commonly used as a barcode. Our protocol proved to be a reliable tool for monitoring population genetic diversity of this species, and could be a valuable addition to amphibian conservation and wetland management.

Molecular analyses of pseudoscorpions in a subterranean biodiversity hotspot reveal cryptic diversity and microendemism

Nested within the Mediterranean biodiversity hotspot, the Dinaric Karst of the western Balkans is one of the world's most heterogeneous subterranean ecosystems and renowned for its highly diverse and mostly endemic fauna. The evolutionary processes leading to both endemism and diversity remain insufficiently understood, and large-scale analyses on taxa that are abundant in both subterranean and surface habitats remain infrequent. Here, we provide the first comprehensive molecular study on Croatian pseudoscorpions, a lineage of arachnids that is common and diverse in both habitats. Phylogenetic reconstructions using 499 COI sequences derived from 128 morphospecies collected across the Dinaric Karst show that: (i) occurrence in karstic microhabitats boosters speciation and endemism in the most diverse genera Chthonius C.L. Koch, 1843 (37 morphospecies) and Neobisium Chamberlin, 1930 (34 morphospecies), (ii) evidence for ongoing diversification is found in many species and species complexes through low optimal thresholds (OTs) and species delineation analyses, and (iii) landscape features, such as mountain ranges, correlate with patterns of genetic diversity in the diverse genus Neobisium. We present two synonymies: Protoneobisium Ćurčić, 1988 = Neobisium, syn. nov., and Archaeoroncus Ćurčić and Rađa, 2012 = Roncus L. Koch, 1873, syn. nov. Overall, our study suggests that karstic microhabitats promote diversification in soil- and cave-dwelling arthropods at all taxonomic levels, but also provide important refugia for invertebrates in past and present periods of environmental change.

Out of thin air: surveying tropical bat roosts through air sampling of eDNA

Understanding roosting behaviour is essential to bat conservation and biomonitoring, often providing the most accurate methods of assessing bat population size and health. However, roosts can be challenging to survey, e.g., physically impossible to access or presenting risks for researchers. Disturbance during monitoring can also disrupt natural bat behaviour and present material risks to the population such as disrupting hibernation cycles. One solution to this is the use of non-invasive monitoring approaches. Environmental (e)DNA has proven especially effective at detecting rare and elusive species particularly in hard-to-reach locations. It has recently been demonstrated that eDNA from vertebrates is carried in air. When collected in semi-confined spaces, this airborne eDNA can provide remarkably accurate profiles of biodiversity, even in complex tropical communities. In this study, we deploy novel airborne eDNA collection for the first time in a natural setting and use this approach to survey difficult to access potential roosts in the neotropics. Using airborne eDNA, we confirmed the presence of bats in nine out of 12 roosts. The identified species matched previous records of roost use obtained from photographic and live capture methods, thus demonstrating the utility of this approach. We also detected the presence of the white-winged vampire bat (Diaemus youngi) which had never been confirmed in the area but was long suspected based on range maps. In addition to the bats, we detected several non-bat vertebrates, including the big-eared climbing rat (Ototylomys phyllotis), which has previously been observed in and around bat roosts in our study area. We also detected eDNA from other local species known to be in the vicinity. Using airborne eDNA to detect new roosts and monitor known populations, particularly when species turnover is rapid, could maximize efficiency for surveyors while minimizing disturbance to the animals. This study presents the first applied use of airborne eDNA collection for ecological analysis moving beyond proof of concept to demonstrate a clear utility for this technology in the wild.

Behavioral observations of the olm (Proteus anguinus) in a karst spring via direct observations and camera trapping

The olm (Proteus anguinus), an endemic amphibian of the Dinarides’ underground waters (Europe), is one of the world’s most widely known subterranean species. Although various aspects of olm biology have been extensively studied, the data on their behavior in the wild remain scarce mostly due to inaccessibility of their natural habitat. Yet, olms also occur in several karstic springs during nighttime. These are easier to access and present an exciting opportunity to study olm behavior in nature. Here, we report on systematic observations of olms in one such spring in Slovenia, where we observed them for nine consecutive summer nights, coupling direct on-site observations with IR camera trap recordings. We used IR camera trap recordings to construct simple ethograms, as well as to quantify olm movement activity by video-tracking. Olms regularly occurred on the surface during the night, and dawn appeared to be a key stimulus for their retreat underground. They were constantly active, but rarely swam far from the spring. Despite the short-term nature of the study, we collected new occurrence and movement data, and at the same time tested the usability of IR cameras for surveying olm presence and behavior in nature. Experience gained through the study may prompt long-term and more complex behavioral studies using similar approaches.

Towards evidence-based conservation of subterranean ecosystems

Subterranean ecosystems are among the most widespread environments on Earth, yet we still have poor knowledge of their biodiversity. To raise awareness of subterranean ecosystems, the essential services they provide, and their unique conservation challenges, 2021 and 2022 were designated International Years of Caves and Karst. As these ecosystems have traditionally been overlooked in global conservation agendas and multilateral agreements, a quantitative assessment of solution-based approaches to safeguard subterranean biota and associated habitats is timely. This assessment allows researchers and practitioners to understand the progress made and research needs in subterranean ecology and management. We conducted a systematic review of peer-reviewed and grey literature focused on subterranean ecosystems globally (terrestrial, freshwater, and saltwater systems), to quantify the available evidence-base for the effectiveness of conservation interventions. We selected 708 publications from the years 1964 to 2021 that discussed, recommended, or implemented 1,954 conservation interventions in subterranean ecosystems. We noted a steep increase in the number of studies from the 2000s while, surprisingly, the proportion of studies quantifying the impact of conservation interventions has steadily and significantly decreased in recent years. The effectiveness of 31% of conservation interventions has been tested statistically. We further highlight that 64% of the reported research occurred in the Palearctic and Nearctic biogeographic regions. Assessments of the effectiveness of conservation interventions were heavily biased towards indirect measures (monitoring and risk assessment), a limited sample of organisms (mostly arthropods and bats), and more accessible systems (terrestrial caves). Our results indicate that most conservation science in the field of subterranean biology does not apply a rigorous quantitative approach, resulting in sparse evidence for the effectiveness of interventions. This raises the important question of how to make conservation efforts more feasible to implement, cost-effective, and long-lasting. Although there is no single remedy, we propose a suite of potential solutions to focus our efforts better towards increasing statistical testing and stress the importance of standardising study reporting to facilitate meta-analytical exercises. We also provide a database summarising the available literature, which will help to build quantitative knowledge about interventions likely to yield the greatest impacts depending upon the subterranean species and habitats of interest. We view this as a starting point to shift away from the widespread tendency of recommending conservation interventions based on anecdotal and expert-based information rather than scientific evidence, without quantitatively testing their effectiveness.

Cultivable Skin Mycobiota of Healthy and Diseased Blind Cave Salamander (Proteus anguinus)

Proteus anguinus is a neotenic cave salamander, endemic to the Dinaric Karst and a symbol of world natural heritage. It is classified as “vulnerable” by the International Union for Conservation of Nature (IUCN) and is one of the EU priority species in need of strict protection. Due to inaccessibility of their natural underground habitat, scientific studies of the olm have been conducted mainly in captivity, where the amphibians are particularly susceptible to opportunistic microbial infections. In this report, we focused on the diversity of cultivable commensal fungi isolated from the skin of asymptomatic and symptomatic animals obtained from nature (20 specimens) and captivity (22 specimens), as well as from underground water of two karstic caves by direct water filtration and by exposure of keratin-based microbial baits and subsequent isolation from them. In total 244 fungal isolates were recovered from the animals and additional 153 isolates were obtained from water samples. Together, these isolates represented 87 genera and 166 species. Symptomatic animals were colonized by a variety of fungal species, most of them represented by a single isolate, including genera known for their involvement in chromomycosis, phaeohyphomycosis and zygomycosis in amphibians: Acremonium, Aspergillus, Cladosporium, Exophiala, Fusarium, Mucor, Ochroconis, Phialophora and Penicillium. One symptomatic specimen sampled from nature was infected by the oomycete Saprolegnia parasitica, the known causative agent of saprolegniosis. This is the first comprehensive report on cultivable skin mycobiome of this unique amphibian in nature and in captivity, with an emphasis on potentially pathogenic fungi and oomycetes.

Environmental DNA of insects and springtails from caves reveals complex processes of eDNA transfer in soils

Subterranean environments host a substantial amount of biodiversity, however assessing the distribution of species living underground is still extremely challenging. Environmental DNA (eDNA) metabarcoding is a powerful tool to estimate biodiversity in poorly known environments and has excellent performance for soil organisms. Here, we tested 1) whether eDNA metabarcoding from cave soils/sediments allows to successfully detect springtails (Hexapoda: Collembola) and insects (Hexapoda: Insecta); 2) whether eDNA mostly represents autochthonous (cave-dwelling) organisms or it also incorporates information from species living in surface environments; 3) whether eDNA detection probability changes across taxa with different ecology. Environmental DNA metabarcoding analyses detected a large number of Molecular Operational Taxonomic Units (MOTUs) for both insects and springtails. For springtails, detection probability was high, with a substantial proportion of hypogean species, suggesting that eDNA provides good information on the distribution of these organisms in caves. Conversely, for insects most of MOTUs represented taxa living outside caves, and the majority of them represented taxa/organisms living in freshwater environments (Ephemeroptera, Plecoptera and Trichoptera). The eDNA of freshwater insects was particularly abundant in deep sectors of caves, far from the entrance. Furthermore, average detection probability of insects was significantly lower than the one of springtails. This suggests that cave soils/sediments act as "conveyer belts of biodiversity information", possibly because percolating water lead to the accumulation of eDNA of organisms living in nearby areas. Cave soils hold a complex mix of autochthonous and allochthonous eDNA. eDNA provided unprecedented information on the understudied subterranean cave organisms; analyses of detection probability and occupancy can help teasing apart local eDNA from the eDNA representing spatially-integrated biodiversity for whole landscape.

Rediscovery and phylogenetic analysis of the Shelta Cave Crayfish (Orconectes sheltae Cooper & Cooper, 1997), a decapod (Decapoda, Cambaridae) endemic to Shelta Cave in northern Alabama, USA

The Shelta Cave Crayfish (Orconectes sheltae) is a small, cave-obligate member of the genus Orconectes (family Cambaridae) endemic to a single cave system—Shelta Cave—in northwest Huntsville, Madison Co., Alabama, USA. Although never abundant, this stygobiont was regularly observed in the 1960s and early 1970s before the population and aquatic community in general at Shelta Cave collapsed likely in response to groundwater contamination and the loss of energetic inputs from a Grey Bat (Myotis grisescens) maternity colony that abandoned the cave after installation of a poorly designed cave gate. We conducted 20 visual surveys of aquatic habitats at Shelta Cave between October 2018 and July 2021. Although the aquatic community has not recovered, we did confirm the continued existence of O. sheltae, which had not been observed in 31 years, with observations of an adult female on 31 May 2019 and an adult male on 28 August 2020. We conducted the first phylogenetic analyses of O. sheltae and discovered that the species is most closely related to other geographically proximate stygobiotic crayfishes in the genus Cambarus in northern Alabama than members of the genus Orconectes. We advocate for recognition of this species as Cambarus sheltae to more accurately reflect evolutionary relationships of this single-cave endemic and offer recommendations for its management, conservation, and future research, as this species remains at high risk of extinction.

eDNA in subterranean ecosystems: Applications, technical aspects, and future prospects

Monitoring of biota is pivotal for the assessment and conservation of ecosystems. Environments worldwide are being continuously and increasingly exposed to multiple adverse impacts, and the accuracy and reliability of the biomonitoring tools that can be employed shape not only the present, but more importantly, the future of entire habitats. The analysis of environmental DNA (eDNA) metabarcoding data provides a quick, affordable, and reliable molecular approach for biodiversity assessments. However, while extensively employed in aquatic and terrestrial surface environments, eDNA-based studies targeting subterranean ecosystems are still uncommon due to the lack of accessibility and the cryptic nature of these environments and their species. Recent advances in genetic and genomic analyses have established a promising framework for shedding new light on subterranean biodiversity and ecology. To address current knowledge and the future use of eDNA methods in groundwaters and caves, this review explores conceptual and technical aspects of the application and its potential in subterranean systems. We briefly introduce subterranean biota and describe the most used traditional sampling techniques. Next, eDNA characteristics, application, and limitations in the subsurface environment are outlined. Last, we provide suggestions on how to overcome caveats and delineate some of the research avenues that will likely shape this field in the near future. We advocate that eDNA analyses, when carefully conducted and ideally combined with conventional sampling techniques, will substantially increase understanding and enable crucial expansion of subterranean community characterisation. Given the importance of groundwater and cave ecosystems for nature and humans, eDNA can bring to the surface essential insights, such as study of ecosystem assemblages and rare species detection, which are critical for the preservation of life below, as well as above, the ground.

Current cave monitoring practices, their variation and recommendations for future improvement in Europe: A synopsis from the 6th EuroSpeleo Protection Symposium

This manuscript summarizes the outcomes of the 6th EuroSpeleo Protection Symposium. Special emphasis was laid on presenting and discussing monitoring activities under the umbrella of the Habitats Directive (EU Council Directive 92/43/EEC) for habitat type 8310 "Caves not open to the public" and the Emerald Network. The discussions revealed a high level of variation in the currently conducted underground monitoring activities: there is no uniform definition of what kind of underground environments the "cave" habitat should cover, how often a specific cave has to be monitored, and what parameters should be measured to evaluate the conservation status. The variation in spatial dimensions in national definitions of caves further affects the number of catalogued caves in a country and the number of caves to be monitored. Not always participants are aware of the complete national monitoring process and that data sets should be freely available or easily accessible. The discussions further showed an inherent dilemma between an anticipated uniform monitoring approach with a coherent assessment methodology and, on the contrary, the uniqueness of caves and subterranean biota to be assessed – combined with profound knowledge gaps and a lack of resources. Nevertheless, some good practices for future cave monitoring activities have been identified by the participants: (1) Cave monitoring should focus on bio- and geodiversity elements alike; (2) Local communities should be involved, and formal agreements envisaged; (3) Caves must be understood as windows into the subterranean realm; (4) Touristic caves should not be excluded ad-hoc from regular monitoring; (5) New digital tools and open FAIR data infrastructures should be implemented; (6) Cave biomonitoring should focus on a large(r) biological diversity; and (7) DNA-based tools should be integrated. Finally, the importance of the 'forgotten' Recommendation No. 36 from the Bern Convention as a guiding legal European document was highlighted.

From Cave Dragons to Genomics: Advancements in the Study of Subterranean Tetrapods

Throughout most of the kingdom Animalia, evolutionary transitions from surface life to a life permanently bound to caves and other subterranean habitats have occurred innumerous times. Not so in tetrapods, where a mere 14 cave-obligate species—all plethodontid and proteid salamanders—are known. We discuss why cave tetrapods are so exceptional and why only salamanders have made the transition. Their evolution follows predictable and convergent, albeit independent pathways. Among the many known changes associated with transitions to subterranean life, eye degeneration, starvation resistance, and longevity are especially relevant to human biomedical research. Recently, sequences of salamander genomes have become available opening up genomic research for cave tetrapods. We discuss new genomic methods that can spur our understanding of the evolutionary mechanisms behind convergent phenotypic change, the relative roles of selective and neutral evolution, cryptic species diversity, and data relevant for conservation such as effective population size and demography.

Inferring predator-prey interaction in the subterranean environment: a case study from Dinaric caves

Predator–prey interactions are among the most important biotic interactions shaping ecological communities and driving the evolution of defensive traits. These interactions and their effects on species received little attention in extreme and remote environments, where possibilities for direct observations and experimental manipulation of the animals are limited. In this paper, we study such type of environment, namely caves of the Dinarides (Europe), combining spatial and phylogenetic methods. We focused on several species of Niphargus amphipods living in phreatic lakes, as some of them use the dorsal spines as putative morphological defensive traits. We predicted that these spines represent a defense strategy against the olm (Proteus anguinus), a top predator species in the subterranean waters. We tested for spatial overlap of the olm and Niphargus species and showed that spined species live in closer proximity to and co-occur more frequently with the olm than non-spined species. Modeling of the evolution of the spines onto Niphargus phylogeny implies coevolution of this trait in the presence of olm. We conclude that these spines likely evolved as defensive traits in a predator–prey arms race. Combining multiple analyses, we provide an example for a methodological framework to assess predator–prey interactions when in-situ or laboratory observations are not possible.

Refining sampling protocols for cavefishes and cave crayfishes to account for environmental variation

Subterranean habitats represent focal habitats in many conservation strategies; however, these environments are some of the most difficult to sample. New sampling methods, such as environmental DNA (eDNA), show promise to improve stygobiont detection, but sources of sampling bias are poorly understood. Therefore, we determined the factors affecting detection probability using traditional visual surveys and eDNA surveys for both cavefishes and cave crayfishes and demonstrated how detection affects survey efforts for these taxa. We sampled 40 sites (179 visual and 183 eDNA surveys) across the Ozark Highlands ecoregion. We estimated the detection probability of cave crayfishes and cavefishes using both survey methods under varying environmental conditions. The effectiveness of eDNA or visual surveys varied by environmental conditions (i.e., water volume, prevailing substrate, and water velocity) and the target taxa. When sampling in areas with average water velocity, no flow, and coarse substrate, eDNA surveys had a higher detection probability (0.49) than visual surveys (0.35) for cavefishes and visual surveys (0.67) had a higher detection probability than eDNA surveys (0.40) for cave crayfishes. Under the same sampling conditions, 5 visual surveys compared to 10 eDNA surveys would be needed to confidently detect cave crayfishes and 9 visual surveys compared to 4 eDNA surveys for cavefishes. Environmental DNA is a complementary tool to traditional visual surveys; however, the limitations we identified indicate eDNA currently cannot replace visual surveys in subterranean environments. Although sampling designs that account for imperfect sampling are particularly useful, they may not be practical; thus, increasing sampling efforts to offset known detection bias would benefit conservation strategies.

Optimising sampling and analysis protocols in environmental DNA studies

Ecological surveys risk incurring false negative and false positive detections of the target species. With indirect survey methods, such as environmental DNA, such error can occur at two stages: sample collection and laboratory analysis. Here we analyse a large qPCR based eDNA data set using two occupancy models, one of which accounts for false positive error by Griffin et al. (J R Stat Soc Ser C Appl Stat 69: 377-392, 2020), and a second that assumes no false positive error by Stratton et al. (Methods Ecol Evol 11: 1113-1120, 2020). Additionally, we apply the Griffin et al. (2020) model to simulated data to determine optimal levels of replication at both sampling stages. The Stratton et al. (2020) model, which assumes no false positive results, consistently overestimated both overall and individual site occupancy compared to both the Griffin et al. (2020) model and to previous estimates of pond occupancy for the target species. The inclusion of replication at both stages of eDNA analysis (sample collection and in the laboratory) reduces both bias and credible interval width in estimates of both occupancy and detectability. Even the collection of > 1 sample from a site can improve parameter estimates more than having a high number of replicates only within the laboratory analysis.

Development and validation of a quantitative qPCR assay for detecting Natterjack toad (Epidalea calamita) eDNA samples

The Natterjack toad (Epidalea calamita) is the rarest amphibian species in Ireland, regionally Red-Listed as Endangered. We applied an eDNA approach to detect species presence in breeding pond water samples. We developed a species-specific qPCR assay targeting the cytochrome c oxidase subunit I (COI). The assay was tested in silico, in vitro (DNA extracted from tissue) and in vivo (DNA extracted from water samples). Water samples were collected from five ponds with known Natterjack toad presence or absence to validate the sensitivity and specificity of the assay. The assay was shown to be highly specific to the Natterjack toad and tested positive only against toad tissue samples and eDNA samples from ponds with known species presence. We believe this method can be used for rapid assessment of species occurrence.

Development and application of eDNA-based tools for the conservation of white-clawed crayfish

eDNA-based methods represent non-invasive and cost-effective approaches for species monitoring and their application as a conservation tool has rapidly increased within the last decade. Currently, they are primarily used to determine the presence/absence of invasive, endangered or commercially important species, but they also hold potential to contribute to an improved understanding of the ecological interactions that drive species distributions. However, this next step of eDNA-based applications requires a thorough method development. We developed an eDNA assay for the white-clawed crayfish (Austropotamobius pallipes), a flagship species of conservation in the UK and Western Europe. Multiple subsequent in-situ and ex-situ validation tests aimed at improving method performance allowed us to apply eDNA-based surveys to evaluate interactions between white-clawed crayfish, crayfish plague and invasive signal crayfish. The assay performed well in terms of specificity (no detection of non-target DNA) and sensitivity, which was higher compared to traditional methods (in this case torching). The eDNA-based quantification of species biomass was, however, less reliable. Comparison of eDNA sampling methods (precipitation vs. various filtration approaches) revealed that optimal sampling method differed across environments and might depend on inhibitor concentrations. Finally, we applied our methodology together with established assays for crayfish plague and the invasive signal crayfish, demonstrating their significant interactions in a UK river system. Our analysis highlights the importance of thorough methodological development of eDNA-based assays. Only a critical evaluation of methodological strengths and weaknesses will allow us to capitalise on the full potential of eDNA-based methods and use them as decision support tools in environmental monitoring and conservation practice.

Using environmental DNA methods to survey for rare groundwater fauna: Detection of an endangered endemic cave crayfish in northern Alabama

The conservation and management of subterranean biodiversity is hindered by a lack of knowledge on the true distributions for many species, e.g., the Wallacean shortfall. In recent years, several studies have demonstrated the potential of environmental DNA (eDNA) as an effective approach to detect and monitor biodiversity, including rare, threatened, and endangered taxa. However, there are few eDNA studies of groundwater fauna. Here we report the results of the development and implementation of an eDNA assay targeting a short fragment of the mitochondrial CO1 locus of a critically imperiled cave crayfish, the Sweet Home Alabama Cave Crayfish (Cambarus speleocoopi), known from just four cave systems in the Interior Plateau karst region of northern Alabama. We detected C. speleocoopi DNA from water samples collected at 5 of 16 sites sampled (caves and springs), including two historical sites as well as three additional and potentially new sites in Marshall County, Alabama. All three of these sites were within 2 km of historical sites. Our study is the first to detect a groundwater crustacean in the Interior Plateau karst region. Additionally, our study contributes to the growing literature that eDNA is a viable complementary tool for detection and monitoring of a fauna that is difficult to survey and study using traditional approaches.

Seasonal monitoring of Hida salamander Hynobius kimurae using environmental DNA with a genus-specific primer set

The diversity and the abundance of amphibians have dramatically declined globally over the past 30 years, and the monitoring and conservation of their habitats is essential. However, traditional methods such as bait trapping and mark-recapture are costly, and morphological identification usually requires a high level of taxonomic expertise. Here, seasonal surveillances of Hida salamanderHynobius kimuraewere performed by means of environmental DNA (eDNA) analysis withHynobius-specific primers and a species-specific TaqMan probe. Water sampling and visual surveys were conducted seasonally in a stream in Kyoto Prefecture, Japan. Detection rates of eDNA were then calculated by real-time PCR, and eDNA site occupancy probability was estimated by multi-scale occupancy modeling. The eDNA-based detection rate of Hida salamander was 76.7%, whereas the visual survey-based detection rate was 23.3%, and target eDNA was detected at almost all sites where the presence of target species was visually confirmed. Moreover, factors relating to the site- and sample-level occurrence probabilities of the target eDNA differed depending on the developmental stage of the target species. Our findings support previous studies showing that eDNA analysis enables an effective assessment of amphibian distributions without damaging the organisms or their habitat, and we compare for the first time the site occupancy probability of amphibian eDNA throughout the life cycle of an amphibian species. The present study contributes to the development of eDNA analysis as a tool for understanding the distribution and seasonal activity of amphibian species and will thus aid in the planning of conservation measures and habitat restoration for these species.

Dampened virulence and limited proliferation of Batrachochytrium salamandrivorans during subclinical infection of the troglobiont olm (Proteus anguinus)

Emerging infections add to existing threats to the survival of amphibians worldwide. The olm (Proteus anguinus) is a vulnerable, troglobiont urodele species with a small European range and restricted to underground karstic systems. Population declines to emerging threats like the chytrid fungus Batrachochytrium salamandrivorans, are likely to go unnoticed due to inaccessibility of the species’ habitat. We here studied the interaction between olms and B. salamandrivorans. Experimental inoculation of olms resulted in low-level, asymptomatic but persistent infections, with limbs as predilection sites. The lack of exponential fungal growth in the olms’ epidermis correlated with limited fungal proliferation and dampened virulence gene expression after exposure to olm skin compounds. The olm is one of few western Palearctic urodeles that is tolerant to B. salamandrivorans infection and may act as a subterranean disease reservoir, yet costs of subclinical infection may compromise olm fitness on the long term.

A new non-invasive in situ underwater DNA sampling method for estimating genetic diversity

DNA-based methods form the cornerstone of contemporary evolutionary biology and they are highly valued tools in conservation biology. The development of non-invasive sampling methods can be crucial for both gathering sample sizes needed for robust ecological inference and to avoid a negative impact on small and/or endangered populations. Such sampling is particularly challenging in working with aquatic organisms, if the goal is to minimize disturbance and to avoid even temporary removal of individuals from their home range. We developed an in situ underwater method of DNA sampling and preservation that can be applied during diving in less than a minute of animal handling. We applied the method on a Herzegovinian population of olm (Proteus anguinus, Caudata), an endangered aquatic cave-dwelling vertebrate, which makes it an excellent model to test the method under the harshest conditions. We sampled 22 adults during cave-diving and extracted sufficient quantity and quality of DNA from all individuals. We amplified 10 species-specific microsatellite loci, with PCR success varying between 6 and 10 loci (median: 7 loci). Fragment length analyses on 9 loci revealed a single allele at all loci across all individuals. This is in stark contrast to four Croatian populations studied with the same 10 loci previously that showed high within-population genetic variation. Our population and the four Croatian populations were genetically highly divergent. We propose that our method can be widely used to sample endangered aquatic populations, or in projects where the disturbance of individuals must be kept minimal for conservation and scientific purposes.

Cave Communities: From the Surface Border to the Deep Darkness

The discipline of subterranean biology has provided us incredible information on the diversity, ecology and evolution of species living in different typologies of subterranean habitats. However, a general lack of information on the relationships between cave species still exists, leaving uncertainty regarding the dynamics that hold together cave communities and the roles of specific organisms (from the least to the most adapted species) for the community, as well as the entire ecosystem. This Special Issue aims to stimulate and gather studies which are focusing on cave communities belonging to all different typologies of subterranean habitats, with the overarching goal to corroborate the key role of the subterranean biology in ecological and evolutionary studies.

eDNA metabarcoding survey reveals fine-scale coral reef community variation across a remote, tropical island ecosystem

Environmental DNA (eDNA) metabarcoding, a technique for retrieving multispecies DNA from environmental samples, can detect a diverse array of marine species from filtered seawater samples. There is a growing potential to integrate eDNA alongside existing monitoring methods in order to establish or improve the assessment of species diversity. Remote island reefs are increasingly vulnerable to climate-related threats and as such there is a pressing need for cost-effective whole-ecosystem surveying to baseline biodiversity, study assemblage changes and ultimately develop sustainable management plans. We investigated the utility of eDNA metabarcoding as a high-resolution, multitrophic biomonitoring tool at the Cocos (Keeling) Islands, Australia (CKI)-a remote tropical coral reef atoll situated within the eastern Indian Ocean. Metabarcoding assays targeting the mitochondrial 16S rRNA and CO1 genes, as well as the 18S rRNA nuclear gene, were applied to 252 surface seawater samples collected from 42 sites within a 140 km² area. Our assays successfully detected a wide range of bony fish and elasmobranchs (244 taxa), crustaceans (88), molluscs (37) and echinoderms (7). Assemblage composition varied significantly between sites, reflecting habitat partitioning across the island ecosystem and demonstrating the localisation of eDNA signals, despite extensive tidal and oceanic movements. In addition, we document putative new occurrence records for 46 taxa and compare the efficiency of our eDNA approach to visual survey techniques at CKI. Our study demonstrates the utility of a multimarker metabarcoding approach in capturing multitrophic biodiversity across an entire coral reef atoll and sets an important baseline for ongoing monitoring and management.

Combining ddPCR and environmental DNA to improve detection capabilities of a critically endangered freshwater invertebrate

Isogenus nubecula is a critically endangered Plecoptera species. Considered extinct in the UK, I. nubecula was recently rediscovered (in one location of the River Dee, Wales), after 22 years of absence. In a similar way to many other species of Perlodidae, I. nubecula could be utilised as a bio-indicator, for assessing water quality and health status of a given freshwater system. However, conventional monitoring of invertebrates via kick-sampling, is invasive and expensive (time consuming). Further, such methods require a high level of taxonomic expertise. Here, we compared the traditional kick-sampling method with the use of eDNA detection using qPCR and ddPCR-analyses. In spring 2018, we sampled eDNA from twelve locations on the River Dee. I. nubecula was detected using kick-sampling in five of these locations, three locations using both eDNA detection and kick-sampling and one location using eDNA detection alone – resulting in a total of six known and distinct populations of this critically endangered species. Interestingly, despite the eDNA assay being validated in vitro and in silico, and results indicating high sensitivity, qPCR analysis of the eDNA samples proved to be ineffective. In contrast, ddPCR analyses resulted in a clear detection of I. nubecula at four locations suggesting that inhibition most likely explains the large discrepancy between the obtained qPCR and ddPCR results. It is therefore important to explore inhibition effects on any new eDNA assay. We also highlight that ddPCR may well be the best option for the detection of aquatic organisms which are either rare or likely to shed low levels of eDNA into their environment.

Concurrent visual encounter sampling validates eDNA selectivity and sensitivity for the endangered wood turtle (Glyptemys insculpta)

Environmental DNA (eDNA) has been used to record the presence of many different organisms in several different aquatic and terrestrial environments. Although eDNA has been demonstrated as a useful tool for the detection of invasive and/or cryptic and declining species, this approach is subject to the same considerations that limit the interpretation of results from traditional survey techniques (e.g. imperfect detection). The wood turtle is a cryptic semi-aquatic species that is declining across its range and, like so many chelonian species, is in-need of a rapid and effective method for monitoring distribution and abundance. To meet this need, we used an eDNA approach to sample for wood turtle presence in northern Virginia streams. At the same time, we used repeat visual encounter surveys in an occupancy-modelling framework to validate our eDNA results and reveal the relationship of detection and occupancy for both methods. We sampled 37 stream reaches of varying size within and beyond the known distribution of the wood turtle across northern Virginia. Wood turtle occupancy probability was 0.54 (0.31, 0.76) and while detection probability for wood turtle occupancy was high (0.88; 0.58, 0.98), our detection of turtle abundance was markedly lower (0.28; 0.21, 0.37). We detected eDNA at 76% of sites confirmed occupied by VES and at an additional three sites where turtles were not detected but were known to occur. Environmental DNA occupancy probability was 0.55 (0.29, 0.78); directly comparable to the VES occupancy estimate. Higher probabilities of detecting wood turtle eDNA were associated with higher turtle densities, an increasing number of days since the last rainfall, lower water temperatures, and lower relative discharges. Our results suggest that eDNA technology holds promise for sampling aquatic chelonians in some systems, even when discharge is high and biomass is relatively low, when the approach is validated and sampling error is quantified.

Detection of elusive fire salamander larvae (Salamandra salamandra) in streams via environmental DNA

In the face of the global biodiversity crisis, the monitoring of species richness and diversity is experiencing an increased demand entailing a raise in cost and time investment. The analysis of species-specific DNA fragments in environmental samples (eDNA) such as from water or soil, facilitate the molecular detection of species without the specific sampling of individuals. The invasive chytrid fungus Batrachochytrium salamandrivorans (Bsal) is infecting natural fire salamander populations (Salamandra salamandra) and causes chytridiomycosis resulting in infrequent regional extinctions of populations across Central Europe. With regard to the expanding distribution of Bsal over the last years, cost-effective monitoring of fire salamanders is important for the conservation of this species. Based on a real-time quantitative PCR (qPCR) assay, we developed a new protocol to detect S. salamandra larvae in streams via eDNA, using species-specific primers of the mitochondrial control region (D-loop). We tested the efficiency of qPCR primer sets for six combinations of DNA extraction kits coupled with subsequent PCR inhibitor removal kits for obtaining qPCR-detectable S. salamandra eDNA from water filters, that were taken both from natural streams and artificial water tanks in the laboratory as positive controls. We found that the DNeasy Blood & Tissue Kit in combination with the DNeasy PowerClean CleanUp Kit performed best for detecting salamander larvae from natural streams. Our experimental protocol paves the way for resource-saving approaches to monitor S. salamandra larvae, but also confirms the limits to this eDNA approach in that it requires optimized laboratory protocols.

Conventional versus real-time quantitative PCR for rare species detection

Detection of species in nature at very low abundance requires innovative methods. Conventional PCR (cPCR) and real-time quantitative PCR (qPCR) are two widely used approaches employed in environmental DNA (eDNA) detection, though lack of a comprehensive comparison of them impedes method selection. Here we test detection capacity and false negative rate of both approaches using samples with different expected complexities. We compared cPCR and qPCR to detect invasive, biofouling golden mussels (Limnoperna fortunei), in samples from laboratory aquaria and irrigation channels where this mussel was known to occur in central China. Where applicable, the limit of detection (LoD), limit of quantification (LoQ), detection rate, and false negative rate of each PCR method were tested. Quantitative PCR achieved a lower LoD than cPCR (1 × 10-7 vs. 10-6 ng/μl) and had a higher detection rate for both laboratory (100% vs. 87.9%) and field (68.6% vs. 47.1%) samples. Field water samples could only be quantified at a higher concentration than laboratory aquaria and total genomic DNA, indicating inhibition with environmental samples. The false negative rate was inversely related to the number of sample replicates. Target eDNA concentration was negatively related to distance from sampling sites to the water (and animal) source. Detection capacity difference between cPCR and qPCR for genomic DNA and laboratory aquaria can be translated to field water samples, and the latter should be prioritized in rare species detection. Field environmental samples may involve more complexities-such as inhibitors-than laboratory aquaria samples, requiring more target DNA. Extensive sampling is critical in field applications using either approach to reduce false negatives.

Ecological patterns in anchialine caves

Anchialine caves are characterized by high levels of endemism and extreme conditions. However, few ecological studies have been conducted in these ecosystems. This study integrates biotic and abiotic parameters of two sets of cave systems with contrasting high and low species richness. Seven ecological patterns are used to explain the expected species richness and density in an anchialine cave. In addition, the population size for conspicuous macrofauna was estimated. The high impact that single-events have on anchialine fauna are also reported. These findings reinforce the conclusions of previous studies of the high extinction risk of anchialine cave fauna, and substantiate the necessity of ad hoc conservation strategies for anchialine caves.

Is the Italian stream frog (Rana italica Dubois, 1987) an opportunistic exploiter of cave twilight zone?

Studies on frogs exploiting subterranean environments are extremely scarce, as these Amphibians are usually considered accidental in these environments. However, according to recent studies, some anurans actively select subterranean environments on the basis of specific environmental features, and thus are able to inhabit these environments throughout the year. We present the first study on the abundance and spatial use of the Italian stream frog, Ranaitalica, in subterranean environments. We monthly collected data from 66 cave sectors during a whole year (2013), recording > 120 detections of R.italica. Frogs were more frequently found close to the cave entrance, without significant differences between age classes or sexes. Adults generally were observed being higher up along cave walls compared to juveniles. Frogs abundance was higher in areas showing specific environmental features, such as warm temperature, low incident light and the presence of potential prey. Ranaitalica likely occupies subterranean areas characterized by a combination of microclimatic suitability and prey availability.

Seasonal variation in environmental DNA detection in sediment and water samples

The use of aquatic environmental DNA (eDNA) to detect the presence of species depends on the seasonal activity of the species in the sampled habitat. eDNA may persist in sediments for longer than it does in water, and analysing sediment could potentially extend the seasonal window for species assessment. Using the great crested newt as a model, we compare how detection probability changes across the seasons in eDNA samples collected from both pond water and pond sediments. Detection of both aquatic and sedimentary eDNA varied through the year, peaking in the summer (July), with its lowest point in the winter (January): in all seasons, detection probability of eDNA from water exceeded that from sediment. Detection probability of eDNA also varied between study areas, and according to great crested newt habitat suitability and sediment type. As aquatic and sedimentary eDNA show the same seasonal fluctuations, the patterns observed in both sample types likely reflect current or recent presence of the target species. However, given the low detection probabilities found in the autumn and winter we would not recommend using either aquatic or sedimentary eDNA for year-round sampling without further refinement and testing of the methods.

Detection of the European pond turtle (Emys orbicularis) by environmental DNA: is eDNA adequate for reptiles?

Recent studies have demonstrated the potential of combining molecular technologies with environmental sampling to detect various vertebrate species in aquatic ecosystems. The European pond turtle (Emys orbicularis) is a threatened and elusive aquatic reptile with shy behaviour. We aimed to develop and evaluate a methodology to detect the presence of this secretive aquatic reptile in ponds from environmental water samples. First, we determined that reptilian DNA can be isolated and amplified from water samples in artificial and natural ponds with known turtle density. Then we compared the potential of two water sampling methods (through filtration or precipitation) and found no significant differences between these approaches. Finally, we demonstrated that the eDNA concentration detected is not correlated with the number of E. orbicularis individuals or biomass. Detection of eDNA was higher in artificial ponds with small volumes of water or in the shallow waters of natural ponds. The eDNA-based methodology aims to detect the presence of specific species, even at low density, with better accuracy than visual observation. However, our study indicates that this method of population monitoring should be applied with caution to aquatic reptiles.

ednaoccupancy: An r package for multiscale occupancy modelling of environmental DNA data

In this article, we describe ednaoccupancy, an r package for fitting Bayesian, multiscale occupancy models. These models are appropriate for occupancy surveys that include three nested levels of sampling: primary sample units within a study area, secondary sample units collected from each primary unit and replicates of each secondary sample unit. This design is commonly used in occupancy surveys of environmental DNA (eDNA). ednaoccupancy allows users to specify and fit multiscale occupancy models with or without covariates, to estimate posterior summaries of occurrence and detection probabilities, and to compare different models using Bayesian model-selection criteria. We illustrate these features by analysing two published data sets: eDNA surveys of a fungal pathogen of amphibians and eDNA surveys of an endangered fish species.

Environmental DNA metabarcoding: Transforming how we survey animal and plant communities

The genomic revolution has fundamentally changed how we survey biodiversity on earth. High-throughput sequencing ("HTS") platforms now enable the rapid sequencing of DNA from diverse kinds of environmental samples (termed "environmental DNA" or "eDNA"). Coupling HTS with our ability to associate sequences from eDNA with a taxonomic name is called "eDNA metabarcoding" and offers a powerful molecular tool capable of noninvasively surveying species richness from many ecosystems. Here, we review the use of eDNA metabarcoding for surveying animal and plant richness, and the challenges in using eDNA approaches to estimate relative abundance. We highlight eDNA applications in freshwater, marine and terrestrial environments, and in this broad context, we distill what is known about the ability of different eDNA sample types to approximate richness in space and across time. We provide guiding questions for study design and discuss the eDNA metabarcoding workflow with a focus on primers and library preparation methods. We additionally discuss important criteria for consideration of bioinformatic filtering of data sets, with recommendations for increasing transparency. Finally, looking to the future, we discuss emerging applications of eDNA metabarcoding in ecology, conservation, invasion biology, biomonitoring, and how eDNA metabarcoding can empower citizen science and biodiversity education.

Is the detection of aquatic environmental DNA influenced by substrate type?

The use of environmental DNA (eDNA) to assess the presence-absence of rare, cryptic or invasive species is hindered by a poor understanding of the factors that can remove DNA from the system. In aquatic systems, eDNA can be transported out either horizontally in water flows or vertically by incorporation into the sediment. Equally, eDNA may be broken down by various biotic and abiotic processes if the target organism leaves the system. We use occupancy modelling and a replicated mesocosm experiment to examine how detection probability of eDNA changes once the target species is no longer present. We hypothesise that detection probability falls faster with a sediment which has a large number of DNA binding sites such as topsoil or clay, over lower DNA binding capacity substrates such as sand. Water removed from ponds containing the target species (the great crested newt) initially showed high detection probabilities, but these fell to between 40% and 60% over the first 10 days and to between 10% and 22% by day 15: eDNA remained detectable at very low levels until day 22. Very little difference in detection was observed between the control group (no substrate) and the sand substrate. A small reduction in detection probability was observed between the control and clay substrates, but this was not significant. However, a highly significant reduction in detection probability was observed with a topsoil substrate. This result is likely to have stemmed from increased levels of PCR inhibition, suggesting that incorporation of DNA into the sentiment is of only limited importance. Surveys of aquatic species using eDNA clearly need to take account of substrate type as well as other environmental factors when collecting samples, analysing data and interpreting the results.

Environmental DNA in subterranean biology: range extension and taxonomic implications for Proteus

Europe’s obligate cave-dwelling amphibian Proteus anguinus inhabits subterranean waters of the north-western Balkan Peninsula. Because only fragments of its habitat are accessible to humans, this endangered salamander’s exact distribution has been difficult to establish. Here we introduce a quantitative real time polymerase chain reaction-based environmental DNA (eDNA) approach to detect the presence of Proteus using water samples collected from karst springs, wells or caves. In a survey conducted along the southern limit of its known range, we established a likely presence of Proteus at seven new sites, extending its range to Montenegro. Next, using specific molecular probes to discriminate the rare black morph of Proteus from the closely related white morph, we detected its eDNA at five new sites, thus more than doubling the known number of sites. In one of these we found both black and white Proteus eDNA together. This finding suggests that the two morphs may live in contact with each other in the same body of groundwater and that they may be reproductively isolated species. Our results show that the eDNA approach is suitable and efficient in addressing questions in biogeography, evolution, taxonomy and conservation of the cryptic subterranean fauna.