Factors influencing detection of eDNA from a stream‐dwelling amphibian

Environmental DNA (eDNA) methods for detecting and estimating abundance of aquatic species are emerging rapidly, but little is known about how processes such as secretion rate, environmental degradation, and time since colonization or extirpation from a given site affect eDNA measurements. Using stream‐dwelling salamanders and quantitative PCR (qPCR) analysis, we conducted three experiments to assess eDNA: (i) production rate; (ii) persistence time under different temperature and light conditions; and (iii) detectability and concentration through time following experimental introduction and removal of salamanders into previously unoccupied streams. We found that 44–50 g individuals held in aquaria produced 77 ng eDNA/h for 2 h, after which production either slowed considerably or began to equilibrate with degradation. eDNA in both full‐sun and shaded treatments degraded exponentially to <1% of the original concentration after 3 days. eDNA was no longer detectable in full‐sun samples after 8 days, whereas eDNA was detected in 20% of shaded samples after 11 days and 100% of refrigerated control samples after 18 days. When translocated into unoccupied streams, salamanders were detectable after 6 h, but only when densities were relatively high (0.2481 individuals/m²) and when samples were collected within 5 m of the animals. Concentrations of eDNA detected were very low and increased steadily from 6–24 h after introduction, reaching 0.0022 ng/L. Within 1 h of removing salamanders from the stream, eDNA was no longer detectable. These results suggest that eDNA detectability and concentration depend on production rates of individuals, environmental conditions, density of animals, and their residence time.     

The importance of including imperfect detection models in eDNA experimental design

Environmental DNA (eDNA) is DNA that has been isolated from field samples, and it is increasingly used to infer the presence or absence of particular species in an ecosystem. However, the combination of sampling procedures and subsequent molecular amplification of eDNA can lead to spurious results. As such, it is imperative that eDNA studies include a statistical framework for interpreting eDNA presence/absence data. We reviewed published literature for studies that utilized eDNA where the species density was known and compared the probability of detecting the focal species to the sampling and analysis protocols. Although biomass of the target species and the volume per sample did not impact detectability, the number of field replicates and number of samples from each replicate were positively related to detection. Additionally, increased number of PCR replicates and increased primer specificity significantly increased detectability. Accordingly, we advocate for increased use of occupancy modelling as a method to incorporate effects of sampling effort and PCR sensitivity in eDNA study design. Based on simulation results and the hierarchical nature of occupancy models, we suggest that field replicates, as opposed to molecular replicates, result in better detection probabilities of target species.     

Patterns of population genetic variation in sympatric chiltoniid amphipods within a calcrete aquifer reveal a dynamic subterranean environment

Calcrete aquifers from the Yilgarn region of arid central Western Australia contain an assemblage of obligate groundwater invertebrate species that are each endemic to single aquifers. Fine-scale phylogeographic and population genetic analyses of three sympatric and independently derived species of amphipod (Chiltoniidae) were carried out to determine whether there were common patterns of population genetic structure or evidence for past geographic isolation of populations within a single calcrete aquifer. Genetic diversity in amphipod mitochondrial DNA (cytochrome c oxidase subunit I gene) and allozymes were examined across a 3.5 km² region of the Sturt Meadows calcrete, which contains a grid of 115 bore holes (=wells). Stygobiont amphipods were found to have high levels of mitochondrial haplotype diversity coupled with low nucleotide diversity. Mitochondrial phylogeographic structuring was found between haplogroups for one of the chiltoniid species, which also showed population structuring for nuclear markers. Signatures of population expansion in two of the three species, match previous findings for diving beetles at the same site, indicating that the system is dynamic. We propose isolation of populations in refugia within the calcrete, followed by expansion events, as the most likely source of intraspecific genetic diversity, due to changes in water level influencing gene flow across the calcrete.     

Residual eDNA detection sensitivity assessed by quantitative real‐time PCR in a river ecosystem

Several studies have demonstrated that environmental DNA (eDNA) can be used to detect the presence of aquatic species, days to weeks after the target species has been removed. However, most studies used eDNA analysis in lentic systems (ponds or lakes), or in controlled laboratory experiments. While eDNA degrades rapidly in all aquatic systems, it also undergoes dilution effects and physical destruction in flowing systems, complicating detection in rivers. However, some eDNA (i.e. residual eDNA) can be retained in aquatic systems, even those subject to high flow regimes. Our goal was to determine residual eDNA detection sensitivity using quantitative real‐time polymerase chain reaction (qRT–PCR), in a flowing, uncontrolled river after the eDNA source was removed from the system; we repeated the experiment over 2 years. Residual eDNA had the strongest signal strength at the original source site and was detectable there up to 11.5 h after eDNA source removal. Residual eDNA signal strength decreased as sampling distance downstream from the eDNA source site increased, and was no longer detectable at the source site 48 h after the eDNA source water was exhausted in both experiments. This experiment shows that residual eDNA sampled in surface water can be mapped quantitatively using qRT–PCR, which allows a more accurate spatial identification of the target species location in lotic systems, and relative residual eDNA signal strength may allow the determination of the timing of the presence of target species.     

Aliens in caves: the global dimension of biological invasions in subterranean ecosystems

Alien species are a significant threat to natural ecosystems and human economies. Despite global efforts to address this challenge, the documented number of alien species is rapidly increasing worldwide. However, the magnitude of the impact of alien species may vary significantly across habitats. For example, some habitats are naturally less prone to biological invasions due to stringent abiotic and biotic characteristics, selecting for a limited number of introduced species possessing traits closely related to the native organisms. Subterranean ecosystems are quintessential examples of habitats with strong environmental filters (e.g. lack of light and scarcity of food), driving convergent adaptations in species that have successfully adapted to life in darkness. Despite these stringent environmental constraints, the number of records of alien species in subterranean ecosystems has increased in recent decades, but the relevant literature remains largely fragmented and mostly anecdotal. Therefore, even though caves are generally considered very fragile ecosystems, their susceptibility to impacts by alien species remains untested other than for some very specific cases. We provide the first systematic literature survey to synthesise available knowledge on alien species in subterranean ecosystems globally. This review is supported by a database summarising the available literature, aiming to identify gaps in the distribution and spread of alien invertebrate species in subterranean habitats, and laying the foundations for future management practices and interventions. First, we quantitatively assessed the current knowledge of alien species in subterranean ecosystems to shed light on broader questions about taxonomic biases, geographical patterns, modes of dispersal, pathways for introductions and potential impacts. Secondly, we collected species-specific traits for each recorded alien species and tested whether subterranean habitats act as ecological filters for their establishment, favouring organisms with pre-adaptive traits suitable for subterranean life. We found information on the presence of 246 subterranean alien species belonging to 18 different classes. The dominant alien species were invertebrates, especially insects and arachnids. Most species were reported in terrestrial subterranean habitats from all continents except Antarctica. Palaearctic and Nearctic biogeographic regions represented the main source of alien species. The main routes of introductions into the recipient country are linked to commercial activities (84.3% of cases for which there was information available). Negative impacts have been documented for a small number of case studies (22.7%), mostly related to increased competition with native species. For a limited number of case studies (6.1%), management strategies were reported but the effectiveness of these interventions has rarely been quantified. Accordingly, information on costs is very limited. Approximately half of the species in our database can be considered established in subterranean habitats. According to our results, the presence of suitable traits grants access to the stringent environmental filter posed by subterranean environments, facilitating establishment in the new habitat. We recommend that future studies deepen the understanding of invasiveness into subterranean habitats, raising public and scientific community awareness of preserving these fragile ecosystems.     

The ancient Britons: groundwater fauna survived extreme climate change over tens of millions of years across NW Europe

Global climate changes during the Cenozoic (65.5–0 Ma) caused major biological range shifts and extinctions. In northern Europe, for example, a pattern of few endemics and the dominance of wide‐ranging species is thought to have been determined by the Pleistocene (2.59–0.01 Ma) glaciations. This study, in contrast, reveals an ancient subsurface fauna endemic to Britain and Ireland. Using a Bayesian phylogenetic approach, we found that two species of stygobitic invertebrates (genus Niphargus) have not only survived the entire Pleistocene in refugia but have persisted for at least 19.5 million years. Other Niphargus species form distinct cryptic taxa that diverged from their nearest continental relative between 5.6 and 1.0 Ma. The study also reveals an unusual biogeographical pattern in the Niphargus genus. It originated in north‐west Europe approximately 87 Ma and underwent a gradual range expansion. Phylogenetic diversity and species age are highest in north‐west Europe, suggesting resilience to extreme climate change and strongly contrasting the patterns seen in surface fauna. However, species diversity is highest in south‐east Europe, indicating that once the genus spread to these areas (approximately 25 Ma), geomorphological and climatic conditions enabled much higher diversification. Our study highlights that groundwater ecosystems provide an important contribution to biodiversity and offers insight into the interactions between biological and climatic processes.     

Genome‐wide adaptive evolution to underground stresses in subterranean mammals: Hypoxia adaption,immunity promotion,and sensory specialization

Life underground has provided remarkable examples of adaptive evolution in subterranean mammals; however, genome‐wide adaptive evolution to underground stresses still needs further research. There are approximately 250 species of subterranean mammals across three suborders and six families. These species not only inhabit hypoxic and dark burrows but also exhibit evolved adaptation to hypoxia, cancer resistance, and specialized sensory systems, making them an excellent model of evolution. The adaptive evolution of subterranean mammals has attracted great attention and needs further study. In the present study, phylogenetic analysis of 5,853 single‐copy orthologous gene families of five subterranean mammals (Nannospalax galili, Heterocephalus glaber, Fukomys damarensis, Condylura cristata, and Chrysochloris asiatica) showed that they formed fou distinct clusters. This result is consistent with the traditional systematics of these species. Furthermore, comparison of the high‐quality genomes of these five subterranean mammalian species led to the identification of the genomic signatures of adaptive evolution. Our results show that the five subterranean mammalian did not share positively selected genes but had similar functional enrichment categories, including hypoxia tolerance, immunity promotion, and sensory specialization, which adapted to the environment of underground stresses. Moreover, variations in soil hardness, climate, and lifestyles have resulted in different molecular mechanisms of adaptation to the hypoxic environment and different degrees of visual degradation. These results provide insights into the genome‐wide adaptive evolution to underground stresses in subterranean mammals, with special focus on the characteristics of hypoxia adaption, immunity promotion, and sensory specialization response to the life underground.     

The distribution of coastal fish eDNA sequences in the Anthropocene

Aim

Coastal fishes have a fundamental role in marine ecosystem functioning and contributions to people, but face increasing threats due to climate change, habitat degradation and overexploitation. The extent to which human pressures are impacting coastal fish biodiversity in comparison with geographic and environmental factors at large spatial scale is still under scrutiny. Here, we took advantage of environmental DNA (eDNA) metabarcoding to investigate the relationship between fish biodiversity, including taxonomic and genetic components, and environmental but also socio-economic factors.

Location

Tropical, temperate and polar coastal areas.

Time period

Present day.

Major taxa studied

Marine fishes.

Methods

We analysed fish eDNA in 263 stations (samples) in 68 sites distributed across polar, temperate and tropical regions. We modelled the effect of environmental, geographic and socio-economic factors on α- and β-diversity. We then computed the partial effect of each factor on several fish biodiversity components using taxonomic molecular units (MOTU) and genetic sequences. We also investigated the relationship between fish genetic α- and β-diversity measured from our barcodes, and phylogenetic but also functional diversity.

Results

We show that fish eDNA MOTU and sequence α- and β-diversity have the strongest correlation with environmental factors on coastal ecosystems worldwide. However, our models also reveal a negative correlation between biodiversity and human dependence on marine ecosystems. In areas with high dependence, diversity of all fish, cryptobenthic fish and large fish MOTUs declined steeply. Finally, we show that a sequence diversity index, accounting for genetic distance between pairs of MOTUs, within and between communities, is a reliable proxy of phylogenetic and functional diversity.

Main conclusions

Together, our results demonstrate that short eDNA sequences can be used to assess climate and direct human impacts on marine biodiversity at large scale in the Anthropocene and can further be extended to investigate biodiversity in its phylogenetic and functional dimensions.     

Evaluating the effects of laboratory protocols on eDNA detection probability for an endangered freshwater fish

The effectiveness and accuracy of detection using environmental DNA (eDNA) is dependent on understanding the influence laboratory methods such as DNA extraction and PCR strategies have on detection probability. Ideally choice of sampling and extraction method will maximize eDNA yield and detection probability. Determining the survey effort required to reach a satisfactory detection probability (via increased PCR replicates or more sampling) could compensate for a lower eDNA yield if the sampling and extraction method has other advantages for a study, species or system. I analysed the effect of three different sampling and extraction methods on eDNA yield, detection probability and PCR replication for detecting the endangered freshwater fish Macquaria australasica from water samples. The impact of eDNA concentration, PCR strategy, target amplicon size and two marker regions: 12S (a mitochondrial gene) and 18S (a nuclear gene) was also assessed. The choice of sampling and extraction method and PCR strategy, rather than amplicon size and marker region, had the biggest effect on detection probability and PCR replication. The PCR replication effort required to achieve a detection probability of 0.95, ranged from 2 to 6 PCR replicates depending on the laboratory method used. As all methods yielded eDNA from which M. australasica was detected using the three target amplicons, differences in eDNA yield and detection probability between the three methods could be mitigated by determining the appropriate PCR replication effort. Evaluating the effect sampling and extraction methods will have on the detection probability and determining the laboratory protocols and PCR replication required to maximize detection and minimize false positives and negatives is a useful first step for eDNA occupancy studies.     

Fine‐tuning for the tropics: application of eDNA technology for invasive fish detection in tropical freshwater ecosystems

Invasive species pose a major threat to aquatic ecosystems. Their impact can be particularly severe in tropical regions, like those in northern Australia, where >20 invasive fish species are recorded. In temperate regions, environmental DNA (eDNA) technology is gaining momentum as a tool to detect aquatic pests, but the technology's effectiveness has not been fully explored in tropical systems with their unique climatic challenges (i.e. high turbidity, temperatures and ultraviolet light). In this study, we modified conventional eDNA protocols for use in tropical environments using the invasive fish, Mozambique tilapia (Oreochromis mossambicus) as a detection model. We evaluated the effects of high water temperatures and fish density on the detection of tilapia eDNA, using filters with larger pores to facilitate filtration. Large‐pore filters (20 μm) were effective in filtering turbid waters and retaining sufficient eDNA, whilst achieving filtration times of 2–3 min per 2‐L sample. High water temperatures, often experienced in the tropics (23, 29, 35 °C), did not affect eDNA degradation rates, although high temperatures (35 °C) did significantly increase fish eDNA shedding rates. We established a minimum detection limit for tilapia (1 fish/0.4 megalitres/after 4 days) and found that low water flow (3.17 L/s) into ponds with high fish density (>16 fish/0.4 megalitres) did not affect eDNA detection. These results demonstrate that eDNA technology can be effectively used in tropical ecosystems to detect invasive fish species.     

Ribosomal DNA non-transcribed spacer polymorphism in subterranean mole rats: genetic differentiation,environmental correlates and phylogenetic relationships

Summary An analysis is presented of genetic differentiation in the non-transcribed spacers of ribosomal DNA (NTS rDNA). Diversity, environmental correlates and the phylogenetic relationships are examined within and between species of the actively speciating subterranean mole rat, superspeciesSpalax ehrenbergi (2n=52, 54, 58, 60) in Israel. This analysis is based on a previous study of the geographic distribution of restriction fragment length polymorphisms of NTS rDNA. Here we present results indicating that NTS rDNA diversity exists mostly (66%) within populations, while 20% is between populations within species, and 14% between species. Multivariate discriminant analysis succeeded in separating 10 of the 13 populations (77%) into their correct chromosomal species, on the basis of the combination of three NTS rDNA repetypes. The phylogenetic relationships suggest that the complex involves two pairs of closely related species (2n=52–54 and 2n=58–60). NTS rDNA diversity, as well as the decrease southward in frequency of repetype C, are correlated with climatic factors of humidity and temperature. These data are discussed in terms of the evolutionary forces of migration and selection which may cause NTS rDNA differentiation. Climatic selection appears to be the major differentiating factor of NTS rDNA.     

Variation of Japanese eel eDNA in sequentially changing conditions and in different sample volumes

Environmental DNA (eDNA) surveys have been conducted to evaluate the distribution and abundance of Japanese eels. However, various environmental and biological factors may influence eDNA concentrations. An experiment was conducted using three water sample replicates (50, 100 and 200 ml) of the same group of eels in a tank that were exposed to sequential nonfeeding/feeding and low/high temperature conditions. Slightly higher concentrations occurred at higher temperature (22–23°C) with nonfeeding, and the highest concentrations occurred when feeding started even though it was in the lower temperature (16–17°C) condition, but sample volume had no effect.     

Methodological considerations for detection of terrestrial small‐body salamander eDNA and implications for biodiversity conservation

Environmental DNA (eDNA) can be used as an assessment tool to detect populations of threatened species and provide fine‐scale data required to make management decisions. The objectives of this project were to use quantitative PCR (qPCR) to: (i) detect spiked salamander DNA in soil, (ii) quantify eDNA degradation over time, (iii) determine detectability of salamander eDNA in a terrestrial environment using soil, faeces, and skin swabs, (iv) detect salamander eDNA in a mesocosm experiment. Salamander eDNA was positively detected in 100% of skin swabs and 66% of faecal samples and concentrations did not differ between the two sources. However, eDNA was not detected in soil samples collected from directly underneath wild‐caught living salamanders. Salamander genomic DNA (gDNA) was detected in all qPCR reactions when spiked into soil at 10.0, 5.0, and 1.0 ng/g soil and spike concentration had a significant effect on detected concentrations. Only 33% of samples showed recoverable eDNA when spiked with 0.25 ng/g soil, which was the low end of eDNA detection. To determine the rate of eDNA degradation, gDNA (1 ng/g soil) was spiked into soil and quantified over seven days. Salamander eDNA concentrations decreased across days, but eDNA was still amplifiable at day 7. Salamander eDNA was detected in two of 182 mesocosm soil samples over 12 weeks (n = 52 control samples; n = 65 presence samples; n = 65 eviction samples). The discrepancy in detection success between experiments indicates the potential challenges for this method to be used as a monitoring technique for small‐bodied wild terrestrial salamander populations.     

Distance,flow and PCR inhibition: eDNA dynamics in two headwater streams

Environmental DNA (eDNA) detection has emerged as a powerful tool for monitoring aquatic organisms, but much remains unknown about the dynamics of aquatic eDNA over a range of environmental conditions. DNA concentrations in streams and rivers will depend not only on the equilibrium between DNA entering the water and DNA leaving the system through degradation, but also on downstream transport. To improve understanding of the dynamics of eDNA concentration in lotic systems, we introduced caged trout into two fishless headwater streams and took eDNA samples at evenly spaced downstream intervals. This was repeated 18 times from mid‐summer through autumn, over flows ranging from approximately 1–96 L/s. We used quantitative PCR to relate DNA copy number to distance from source. We found that regardless of flow, there were detectable levels of DNA at 239.5 m. The main effect of flow on eDNA counts was in opposite directions in the two streams. At the lowest flows, eDNA counts were highest close to the source and quickly trailed off over distance. At the highest flows, DNA counts were relatively low both near and far from the source. Biomass was positively related to eDNA copy number in both streams. A combination of cell settling, turbulence and dilution effects is probably responsible for our observations. Additionally, during high leaf deposition periods, the presence of inhibitors resulted in no amplification for high copy number samples in the absence of an inhibition‐releasing strategy, demonstrating the necessity to carefully consider inhibition in eDNA analysis.     

The application of eDNA for monitoring of the Great Crested Newt in the UK

Current ecological surveys for great crested newts are time‐consuming and expensive and can only be carried out within a short survey window. Additional survey methods which would facilitate the detection of rare or protected species such as the great crested newt (Triturus cristatus) would be extremely advantageous. Environmental DNA (eDNA) analysis has been utilized for the detection of great crested newts in Denmark. Here, the same methodology has been applied to water samples taken from UK ponds concurrently with conventional field surveying techniques. Our eDNA analysis exhibited an 84% success rate with a kappa coefficient of agreement between field and eDNA surveys of 0.86. One pond determined to be negative for great crested newt by field survey was positive by eDNA analysis, revealing the potential for improved detection rates using this methodology. Analysis of water samples collected in late summer indicates that eDNA analysis could be used to detect great crested newt after the optimal survey window for current field techniques had passed. Consequently, eDNA analysis could augment currently stipulated techniques for great crested newt surveying as a relatively quick and inexpensive tool for collecting great crested newt presence and distribution data within the UK instead of or prior to full field surveys.     

Retrospective eDNA assessment of potentially harmful algae in historical ship ballast tank and marine port sediments

Microalgal bloom events can cause major ecosystem disturbances, devastate local marine economies, and endanger public health. Therefore, detecting and monitoring harmful microalgal taxa is essential to ensure effective risk management in waterways used for fisheries, aquaculture, recreational activity, and shipping. To fully understand the current status and future direction of algal bloom distributions, we need to know how populations and ecosystems have changed over time. This baseline knowledge is critical for predicting ecosystem responses to future anthropogenic change and will assist in the future management of coastal ecosystems. We explore a NGS metabarcoding approach to rapidly identify potentially harmful microalgal taxa in 63 historic and modern Australian marine port and ballast tank sediment samples. The results provide a record of past microalgal distribution and important baseline data that can be used to assess the efficacy of shipping guidelines, nutrient pollution mitigation, and predict the impact of climate change. Critically, eDNA surveys of archived sediments were able to detect harmful algal taxa that do not produce microscopic fossils, such as Chattonella, Heterosigma, Karlodinium, and Noctiluca. Our data suggest a potential increase in Australian harmful microalgal taxa over the past 30 years, and confirm ship ballast tanks as key dispersal vectors. These molecular mapping tools will assist in the creation of policies aimed at reducing the global increase and spread of harmful algal taxa and help prevent economic and public‐health problems caused by harmful algal blooms.     

Phylogeographical structure in the subterranean tuco-tuco Ctenomys talarum (Rodentia: Ctenomyidae): contrasting the demographic consequences of regional and habitat-specific histories

In this work we examined the phylogeography of the South American subterranean herbivorous rodent Ctenomys talarum (Talas tuco-tuco) using mitochondrial DNA (mtDNA) control region (D-loop) sequences, and we assessed the geographical genetic structure of this species in comparison with that of subterranean Ctenomys australis, which we have shown previously to be parapatric to C. talarum and to also live in a coastal sand dune habitat. A significant apportionment of the genetic variance among regional groups indicated that putative geographical barriers, such as rivers, substantially affected the pattern of genetic structure in C. talarum. Furthermore, genetic differentiation is consistent with a simple model of isolation by distance, possibly evidencing equilibrium between gene flow and local genetic drift. In contrast, C. australis showed limited hierarchical partitioning of genetic variation and departed from an isolation-by-distance pattern. Mismatch distributions and tests of neutrality suggest contrasting histories of these two species: C. talarum appears to be characterized by demographic stability and no significant departures from neutrality, whereas C. australis has undergone a recent demographic expansion and/or departures from strict neutrality in its mtDNA.     

Unlocking biodiversity and conservation studies in high‐diversity environments using environmental DNA (eDNA): A test with Guianese freshwater fishes

Determining the species compositions of local assemblages is a prerequisite to understanding how anthropogenic disturbances affect biodiversity. However, biodiversity measurements often remain incomplete due to the limited efficiency of sampling methods. This is particularly true in freshwater tropical environments that host rich fish assemblages, for which assessments are uncertain and often rely on destructive methods. Developing an efficient and nondestructive method to assess biodiversity in tropical freshwaters is highly important. In this study, we tested the efficiency of environmental DNA (eDNA) metabarcoding to assess the fish diversity of 39 Guianese sites. We compared the diversity and composition of assemblages obtained using traditional and metabarcoding methods. More than 7,000 individual fish belonging to 203 Guianese fish species were collected by traditional sampling methods, and ~17 million reads were produced by metabarcoding, among which ~8 million reads were assigned to 148 fish taxonomic units, including 132 fish species. The two methods detected a similar number of species at each site, but the species identities partially matched. The assemblage compositions from the different drainage basins were better discriminated using metabarcoding, revealing that while traditional methods provide a more complete but spatially limited inventory of fish assemblages, metabarcoding provides a more partial but spatially extensive inventory. eDNA metabarcoding can therefore be used for rapid and large‐scale biodiversity assessments, while at a local scale, the two approaches are complementary and enable an understanding of realistic fish biodiversity.