Simultaneous species detection and discovery with environmental DNA metabarcoding: A freshwater mollusk case study.

Environmental DNA (eDNA) sampling is a powerful tool for rapidly characterizing biodiversity patterns for specious, cryptic taxa with incomplete taxonomies. One such group that are also of high conservation concern are North American freshwater gastropods. In particular, springsnails of the genus Pyrgulopsis (Family: Hydrobiidae) are prevalent throughout the western United States where >140 species have been described. Many of the described species are narrow endemics known from a single spring or locality, and it is believed that there are likely many additional species which have yet to be described. The distribution of these species across the landscape is of interest because habitat loss and degradation, climate change, groundwater mining, and pollution have resulted in springsnail imperilment rates as high as 92%. Determining distributions with conventional sampling methods is limited by the fact that these snails are often <5 mm in length with few distinguishing morphological characters, making them both difficult to detect and to identify. We developed an eDNA metabarcoding protocol that is both inexpensive and capable of rapid, accurate detection of all known Pyrgulopsis species. When compared with conventional collection techniques, our pipeline consistently resulted in detection at sites previously known to contain Pyrgulopsis springsnails and at a cost per site that is likely to be substantially less than the conventional sampling and individual barcoding that has been done historically. Additionally, because our method uses eDNA extracted from filtered water, it is non-destructive and suitable for the detection of endangered species where "no take" restrictions may be in effect. This effort represents both a tool which is immediately applicable to taxa of high conservation concern across western North America and a case study in the broader application of eDNA sampling for landscape assessments of cryptic taxa of conservation concern.

The unknown unknown: A framework for assessing environmental DNA assay specificity against unsampled taxa.

Taxon-specific quantitative PCR (qPCR) assays are commonly used for environmental DNA sampling-based inference of animal presence. These assays require thorough validation to ensure that amplification truly indicates detection of the target taxon, but a thorough validation is difficult when there are potentially many non-target taxa, some of which may have incomplete taxonomies. Here, we use a previously published, quantitative model of cross-amplification risk to describe a framework for assessing qPCR assay specificity when there is missing information and it is not possible to assess assay specificity for each individual non-target confamilial. In this framework, we predict assay specificity against unsampled taxa (non-target taxa without sequence data available) using the sequence information that is available for other confamilials. We demonstrate this framework using four case study assays for: (1) An endemic, freshwater arthropod (meltwater stonefly; Lednia tumana), (2) a globally distributed, marine ascidian (Didemnum perlucidum), (3) a continentally distributed freshwater crustacean (virile crayfish; Faxonius virilis, deanae and nais species complex) and (4) a globally distributed freshwater teleost (common carp; Cyprinus carpio and its close relative C. rubrofuscus). We tested the robustness of our approach to missing information by simulating application of our framework for all possible subsamples of 20-all non-target taxa. Our results suggest that the modelling framework results in estimates which are largely concordant with observed levels of cross-amplification risk using all available sequence data, even when there are high levels of data missingness. We explore potential limitations and extensions of this approach for assessing assay specificity and provide users with an R Markdown template for generating reproducible reports to support their own assay validation efforts.

Broad-scale eDNA sampling for describing aquatic species distributions in running waters: Pacific lamprey Entosphenus tridentatus in the upper Snake River, USA.

One of the most fundamental yet challenging tasks for aquatic ecologists is to precisely delineate the range of species, particularly those that are broadly distributed, require specialized sampling methods, and may be simultaneously declining and increasing in different portions of their range. An exemplar is the Pacific lamprey Entosphenus tridentatus, a jawless anadromous fish of conservation concern that is actively managed in many coastal basins in western North America. To efficiently determine its distribution across the accessible 56,168 km of the upper Snake River basin in the north-western United States, we first delimited potential habitat by using predictions from a species distribution model based on conventionally collected historical data and from the distribution of a potential surrogate, Chinook salmon Oncorhynchus tshawytscha, which yielded a potential habitat network of 10,615 km. Within this area, we conducted a two-stage environmental DNA survey involving 394 new samples and 187 archived samples collected by professional biologists and citizen scientists using a single, standardized method from 2015 to 2021. We estimated that Pacific lamprey occupied 1875 km of lotic habitat in this basin, of which 1444 km may have been influenced by recent translocation efforts. Pacific lamprey DNA was consistently present throughout most river main stems, although detections became weaker or less frequent in the largest and warmest downstream channels and near their headwater extent. Pacific lamprey were detected in nearly all stocked tributaries, but there was no evidence of indigenous populations in such habitats. There was evidence of post-stocking movement because detections were 1.8-36.0 km upstream from release sites. By crafting a model-driven spatial sampling template and executing an eDNA-based sampling campaign led by professionals and volunteers, supplemented by previously collected samples, we established a benchmark for understanding the current range of Pacific lamprey across a large portion of its range in the interior Columbia River basin. This approach could be tailored to refine range estimates for other wide-ranging aquatic species of conservation concern.

eDNAssay: A machine learning tool that accurately predicts qPCR cross-amplification.

Environmental DNA (eDNA) sampling is a highly sensitive and cost-effective technique for wildlife monitoring, notably through the use of qPCR assays. However, it can be difficult to ensure assay specificity when many closely related species co-occur. In theory, specificity may be assessed in silico by determining whether assay oligonucleotides have enough base-pair mismatches with nontarget sequences to preclude amplification. However, the mismatch qualities required are poorly understood, making in silico assessments difficult and often necessitating extensive in vitro testing-typically the greatest bottleneck in assay development. Increasing the accuracy of in silico assessments would therefore streamline the assay development process. In this study, we paired 10 qPCR assays with 82 synthetic gene fragments for 530 specificity tests using SYBR Green intercalating dye (n = 262) and TaqMan hydrolysis probes (n = 268). Test results were used to train random forest classifiers to predict amplification. The primer-only model (SYBR Green results) and full-assay model (TaqMan probe-based results) were 99.6% and 100% accurate, respectively, in cross-validation. We further assessed model performance using six independent assays not used in model training. In these tests the primer-only model was 92.4% accurate (n = 119) and the full-assay model was 96.5% accurate (n = 144). The high performance achieved by these models makes it possible for eDNA practitioners to more quickly and confidently develop assays specific to the intended target. Practitioners can access the full-assay model online via eDNAssay (https://NationalGenomicsCenter.shinyapps.io/eDNAssay), a user-friendly tool for predicting qPCR cross-amplification.

Do metapopulations and management matter for relict headwater bull trout populations in a warming climate?

Mountain headwater streams have emerged as important climate refuges for native cold-water species due to their slow climate velocities and extreme physical conditions that inhibit non-native invasions. Species persisting in refuges often do so as fragmented, relict populations from broader historical distributions that are subject to ongoing habitat reductions and increasing isolation as climate change progresses. Key for conservation planning is determining where remaining populations will persist and how habitat restoration strategies can improve biological resilience to enhance the long-term prospects for species of concern. Studying bull trout, a headwater species in the northwestern USA, we developed habitat occupancy models using a data set of population occurrence in 991 natal habitat patches with a suite of novel geospatial covariates derived from high-resolution hydroclimatic scenarios and other sources representing watershed and instream habitat conditions, patch geometry, disturbance, and biological interactions. The best model correctly predicted bull trout occupancy status in 82.6% of the patches and included effects for: patch size estimated as habitat volume, extent of within-patch reaches <9°C mean August temperature, distance to nearest occupied patch, road density, invasive brook trout prevalence, patch slope, and frequency of high winter flows. The model was used to assess 16 scenarios of bull trout occurrence within the study streams that represented a range of restoration strategies under three climatic conditions (baseline, moderate change, and extreme change). Results suggested that regional improvements in bull trout status were difficult to achieve in realistic restoration strategies due to the pervasive nature of climate change and the limited extent of restoration actions given their high costs. However, occurrence probabilities in a subset of patches were highly responsive to restoration actions, suggesting that targeted investments to improve the resilience of some populations may be contextually beneficial. A possible strategy, therefore, is focusing effort on responsive populations near more robust population strongholds, thereby contributing to local enclaves where dispersal among populations further enhances resilience. Equally important, strongholds constituted a small numerical percentage of patches (5%-21%), yet encompassed the large majority of occupied habitat by volume (72%-89%) and their protection could have significant conservation benefits for bull trout.

Molecular species delimitation refines the taxonomy of native and nonnative physinine snails in North America.

Being able to associate an organism with a scientific name is fundamental to our understanding of its conservation status, ecology, and evolutionary history. Gastropods in the subfamily Physinae have been especially troublesome to identify because morphological variation can be unrelated to interspecific differences and there have been widespread introductions of an unknown number of species, which has led to a speculative taxonomy. To resolve uncertainty about species diversity in North America, we targeted an array of single-locus species delimitation methods at publically available specimens and new specimens collected from the Snake River basin, USA to generate species hypotheses, corroborated using nuclear analyses of the newly collected specimens. A total-evidence approach delineated 18 candidate species, revealing cryptic diversity within recognized taxa and a lack of support for other named taxa. Hypotheses regarding certain local endemics were confirmed, as were widespread introductions, including of an undescribed taxon likely belonging to a separate genus in southeastern Idaho for which the closest relatives are in southeast Asia. Overall, single-locus species delimitation was an effective first step toward understanding the diversity and distribution of species in Physinae and to guiding future investigation sampling and analyses of species hypotheses.

Pliocene-Early Pleistocene Geological Events Structure Pacific Martens (Martes caurina).

The complex topography, climate, and geological history of Western North America have shaped contemporary patterns of biodiversity and species distributions in the region. Pacific martens (Martes caurina) are distributed along the northern Pacific Coast of North America with disjunct populations found throughout the Northwestern Forested Mountains and Marine West Coast Forest ecoregions of the West Coast. Martes in this region have been classified into subspecies; however, the subspecific designation has been extensively debated. In this study, we use genomic data to delineate conservation units of Pacific marten in the Sierra-Cascade-Coastal montane belt in the western United States. We analyzed the mitochondrial genome for 94 individuals to evaluate the spatial distribution and divergence times of major lineages. We further genotyped 401 individuals at 13 microsatellite loci to investigate major patterns of population structure. Both nuclear and mitochondrial DNA suggest substantial genetic substructure concordant with historical subspecies designations. Our results revealed that the region contains 2 distinct mitochondrial lineages: a Cascades/Sierra lineage that diverged from the Cascades/coastal lineage 2.23 (1.48-3.14 mya), consistent with orogeny of the Cascade Mountain chain. Interestingly, Pacific Martes share phylogeographic patterns similar with other sympatric taxa, suggesting that the complex geological history has shaped the biota of this region. The information is critical for conservation and management efforts, and further investigation of adaptive diversity is warranted following appropriate revision of conservation management designations.

Decreased total psoas muscle area after neoadjuvant therapy is a predictor of increased mortality in patients undergoing oesophageal cancer resection.

BACKGROUND: Oesophagectomy for locally advanced cancer carries high rates of morbidity and mortality. Patients require a thorough risk assessment alongside preoperative counselling. Total psoas area (TPA) measurements have been used as a surrogate marker of sarcopenia to predict post-operative complications in oesophageal cancer patients. No studies to date have determined whether there is an association between the proportion of TPA lost during neoadjuvant therapy and post-operative outcomes. METHODS: Clinical data and imaging of patients who underwent neoadjuvant therapy followed by open two-stage oesophagectomy between January 2008 and April 2018 were analysed retrospectively. Patients who did not undergo restaging computed tomography scan prior to surgery were excluded from the study. The TPA was measured on two cross-sectional slices at L4 on computed tomography scans pre- and post-neoadjuvant therapy. RESULTS: A total of 53 patients who met inclusion criteria were identified. The mean loss of TPA was 7.3%. Patients who had a decrease of TPA of more than 4% had significantly increased 30-day mortality compared to those who lost 4% or less (24% versus 0%, P = 0.02). Patients aged over 65 years who also had a loss of TPA >4% had significantly increased 30-day mortality (37% versus 2.9%, odds ratio 19, P = 0.008). CONCLUSION: A decrease in TPA of >4% is associated with a significantly higher risk of post-operative mortality in patients undergoing neoadjuvant therapy followed by oesophagectomy. Measuring the loss of TPA during neoadjuvant treatment could be a novel aid to preoperative risk assessment.

Integrative taxonomy refutes a species hypothesis: The asymmetric hybrid origin of Arsapnia arapahoe (Plecoptera, Capniidae).

Molecular tools are commonly directed at refining taxonomies and the species that constitute their fundamental units. This has been especially insightful for groups for which species hypotheses are ambiguous and have largely been based on morphological differences between certain life stages or sexes, and has added importance when taxa are a focus of conservation efforts. Here, we examine the taxonomic status of Arsapnia arapahoe, a winter stonefly in the family Capniidae that is a species of conservation concern because of its limited abundance and restricted range in northern Colorado, USA. Phylogenetic analyses of sequences of mitochondrial and nuclear genes of this and other capniid stoneflies from this region and elsewhere in western North America indicated extensive haplotype sharing, limited genetic differences, and a lack of reciprocal monophyly between A. arapahoe and the sympatric A. decepta, despite distinctive and consistent morphological differences in the sexual apparatus of males of both species. Analyses of autosomal and sex-linked single nucleotide polymorphisms detected using genotyping by sequencing indicated that all individuals of A. arapahoe consisted of F1 hybrids between female A. decepta and males of another sympatric stonefly, Capnia gracilaria. Rather than constitute a self-sustaining evolutionary lineage, A. arapahoe appears to represent the product of nonintrogressive hybridization in the limited area of syntopy between two widely distributed taxa. This offers a cautionary tale for taxonomists and conservation biologists working on the less-studied components of the global fauna.

An improved environmental DNA assay for bull trout (Salvelinus confluentus) based on the ribosomal internal transcribed spacer I.

The majority of environmental DNA (eDNA) assays for vertebrate species are based on commonly analyzed regions of the mitochondrial genome. However, the high degree of mitochondrial similarity between two species of charr (Salvelinus spp.), southern Dolly Varden (S. malma lordii) and bull trout (Salvelinus confluentus), precludes the development of a mitochondrial eDNA assay to distinguish them. Presented here is an eDNA assay to detect bull trout based on the first ribosomal internal transcribed spacer (ITSI), a nuclear marker. This assay successfully detects bull trout and avoids detecting Dolly Varden as well as brook trout (S. fontinalis), Arctic char (S. alpinus), and lake trout (S. namaycush). In addition, this assay was compared with an extensively used mitochondrial bull trout assay and it was found that the ITSI-based assay produced higher detectability. Our results suggest this assay should out-perform the published mtDNA assay across the range of bull trout, while the added specificity allows reliable bull trout detection in areas where bull trout co-occur with other charr such as Dolly Varden. While clearly a superior assay in this instance, basing assays on ITSI is not without problems. For vertebrates, there are fewer ITSI sequences available than commonly sequenced regions of the mitochondrial genome. Thus, the initial in silico screening of candidate assays must be preceded by much more extensive sampling and sequencing of sympatric or closely related taxa. Further, all copies of the internal transcribed spacers within an individual may not be identical, which can lead to complications. Lastly, the copy number for ITSI varies widely across taxa; the greater detectability associated with this assay cannot be assumed for other species.

Capture enrichment of aquatic environmental DNA: A first proof of concept.

Environmental DNA (eDNA) sampling-the detection of genetic material in the environment to infer species presence-has rapidly grown as a tool for sampling aquatic animal communities. A potentially powerful feature of environmental sampling is that all taxa within the habitat shed DNA and so may be detectable, creating opportunity for whole-community assessments. However, animal DNA in the environment tends to be comparatively rare, making it necessary to enrich for genetic targets from focal taxa prior to sequencing. Current metabarcoding approaches for enrichment rely on bulk amplification using conserved primer annealing sites, which can result in skewed relative sequence abundance and failure to detect some taxa because of PCR bias. Here, we test capture enrichment via hybridization as an alternative strategy for target enrichment using a series of experiments on environmental samples and laboratory-generated, known-composition DNA mixtures. Capture enrichment resulted in detecting multiple species in both kinds of samples, and postcapture relative sequence abundance accurately reflected initial relative template abundance. However, further optimization is needed to permit reliable species detection at the very low-DNA quantities typical of environmental samples (<0.1 ng DNA). We estimate that our capture protocols are comparable to, but less sensitive than, current PCR-based eDNA analyses.

Repurposing environmental DNA samples-detecting the western pearlshell (Margaritifera falcata) as a proof of concept.

Information on the distribution of multiple species in a common landscape is fundamental to effective conservation and management. However, distribution data are expensive to obtain and often limited to high-profile species in a system. A recently developed technique, environmental DNA (eDNA) sampling, has been shown to be more sensitive than traditional detection methods for many aquatic species. A second and perhaps underappreciated benefit of eDNA sampling is that a sample originally collected to determine the presence of one species can be re-analyzed to detect additional taxa without additional field effort. We developed an eDNA assay for the western pearlshell mussel (Margaritifera falcata) and evaluated its effectiveness by analyzing previously collected eDNA samples that were annotated with information including sample location and deposited in a central repository. The eDNA samples were initially collected to determine habitat occupancy by nonbenthic fish species at sites that were in the vicinity of locations recently occupied by western pearlshell. These repurposed eDNA samples produced results congruent with historical western pearlshell surveys and permitted a more precise delineation of the extent of local populations. That a sampling protocol designed to detect fish was also successful for detecting a freshwater mussel suggests that rapidly accumulating collections of eDNA samples can be repurposed to enhance the efficiency and cost-effectiveness of aquatic biodiversity monitoring.

qPCR detection of Sturgeon chub (Macrhybopsis gelida) DNA in environmental samples.

The Sturgeon chub (Macrhybopsis gelida) is a cyprinid fish native to the Missouri and Mississippi River basins of the U.S. Suspected long-term declines in the size of its distribution have prompted a review of its conservation status by the U.S. Fish and Wildlife Service, a process which depends on reliable methods to delineate the distribution and status of extant populations. To facilitate monitoring of Sturgeon chub populations, we developed a quantitative PCR assay to detect Sturgeon chub DNA in environmental samples. The assay consistently detected Sturgeon chub DNA in concentrations as low as 2 copies per reaction, and did not amplify DNA from non-target fish species that are sympatric in the upper Missouri River basin. Field tests of this assay with environmental samples successfully detected Sturgeon chub from sites known to be occupied. This assay offers an extremely sensitive methodology that can be applied to determine the range of Sturgeon chub, regardless of variation in habitat characteristics.

Ecological segregation moderates a climactic conclusion to trout hybridization.

For decades, it has been assumed that introgressive hybridization between introduced rainbow trout and native cutthroat trout in western North America will lead to genomic extinction of the latter. A broad-scale re-examination of their interaction indicates that ecological differences between these species and demographic processes are dictating the location and extent of their hybrid zones, and that runaway introgression between these taxa is unlikely.

Environmental DNA assays for the sister taxa sauger (Sander canadensis) and walleye (Sander vitreus).

Sauger (Sander canadensis) and walleye (S. vitreus) are percid fishes that naturally co-occur throughout much of the eastern United States. The native range of sauger extends into the upper Missouri River drainage where walleye did not historically occur, but have been stocked as a sport fish. Sauger populations have been declining due to habitat loss, fragmentation, and competition with non-native species, such as walleye. To effectively manage sauger populations, it is necessary to identify areas where sauger occur, and particularly where they co-occur with walleye. We developed quantitative PCR assays that can detect sauger and walleye DNA in filtered water samples. Each assay efficiently detected low quantities of target DNA and failed to detect DNA of non-target species with which they commonly co-occur.

A Noninvasive Tool to Assess the Distribution of Pacific Lamprey (Entosphenus tridentatus) in the Columbia River Basin.

The Pacific lamprey (Entosphenus tridentatus) is an anadromous fish once abundant throughout coastal basins of western North America that has suffered dramatic declines in the last century due primarily to human activities. Here, we describe the development of an environmental DNA (eDNA) assay to detect Pacific lamprey in the Columbia River basin. The eDNA assay successfully amplified tissue derived DNA of Pacific lamprey collected from 12 locations throughout the Columbia River basin. The assay amplifies DNA from other Entosphenus species found outside of the Columbia River basin, but is species-specific within this basin. As a result, the assay presented here may be useful for detecting Entosphenus spp. in geographic range beyond the Columbia River Basin. The assay did not amplify tissue or synthetically derived DNA of 14 commonly sympatric non-target species, including lampreys of the genus Lampetra, which are morphologically similar to Pacific lamprey in the freshwater larval stage.

Big biology meets microclimatology: defining thermal niches of ectotherms at landscape scales for conservation planning.

Temperature profoundly affects ecology, a fact ever more evident as the ability to measure thermal environments increases and global changes alter these environments. The spatial structure of thermalscapes is especially relevant to the distribution and abundance of ectothermic organisms, but the ability to describe biothermal relationships at extents and grains relevant to conservation planning has been limited by small or sparse data sets. Here, we combine a large occurrence database of >23 000 aquatic species surveys with stream microclimate scenarios supported by an equally large temperature database for a 149 000-km mountain stream network to describe thermal relationships for 14 fish and amphibian species. Species occurrence probabilities peaked across a wide range of temperatures (7.0-18.8°C) but distinct warm- or cold-edge distribution boundaries were apparent for all species and represented environments where populations may be most sensitive to thermal changes. Warm-edge boundary temperatures for a native species of conservation concern were used with geospatial data sets and a habitat occupancy model to highlight subsets of the network where conservation measures could benefit local populations by maintaining cool temperatures. Linking that strategic approach to local estimates of habitat impairment remains a key challenge but is also an opportunity to build relationships and develop synergies between the research, management, and regulatory communities. As with any data mining or species distribution modeling exercise, care is required in analysis and interpretation of results, but the use of large biological data sets with accurate microclimate scenarios can provide valuable information about the thermal ecology of many ectotherms and a spatially explicit way of guiding conservation investments.

Climate, Demography, and Zoogeography Predict Introgression Thresholds in Salmonid Hybrid Zones in Rocky Mountain Streams.

Among the many threats posed by invasions of nonnative species is introgressive hybridization, which can lead to the genomic extinction of native taxa. This phenomenon is regarded as common and perhaps inevitable among native cutthroat trout and introduced rainbow trout in western North America, despite that these taxa naturally co-occur in some locations. We conducted a synthetic analysis of 13,315 genotyped fish from 558 sites by building logistic regression models using data from geospatial stream databases and from 12 published studies of hybridization to assess whether environmental covariates could explain levels of introgression between westslope cutthroat trout and rainbow trout in the U.S. northern Rocky Mountains. A consensus model performed well (AUC, 0.78-0.86; classification success, 72-82%; 10-fold cross validation, 70-82%) and predicted that rainbow trout introgression was significantly associated with warmer water temperatures, larger streams, proximity to warmer habitats and to recent sources of rainbow trout propagules, presence within the historical range of rainbow trout, and locations further east. Assuming that water temperatures will continue to rise in response to climate change and that levels of introgression outside the historical range of rainbow trout will equilibrate with those inside that range, we applied six scenarios across a 55,234-km stream network that forecast 9.5-74.7% declines in the amount of habitat occupied by westslope cutthroat trout populations of conservation value, but not the wholesale loss of such populations. We conclude that introgression between these taxa is predictably related to environmental conditions, many of which can be manipulated to foster largely genetically intact populations of westslope cutthroat trout and help managers prioritize conservation activities.

Correction: Environmental DNA Marker Development with Sparse Biological Information: A Case Study on Opossum Shrimp (Mysis diluviana).

[This corrects the article DOI: 10.1371/journal.pone.0161664.].

Quantitative PCR Assays for Detecting Loach Minnow (Rhinichthys cobitis) and Spikedace (Meda fulgida) in the Southwestern United States.

Loach minnow (Rhinichthys cobitis) and spikedace (Meda fulgida) are legally protected with the status of Endangered under the U.S. Endangered Species Act and are endemic to the Gila River basin of Arizona and New Mexico. Efficient and sensitive methods for monitoring these species' distributions are critical for prioritizing conservation efforts. We developed quantitative PCR assays for detecting loach minnow and spikedace DNA in environmental samples. Each assay reliably detected low concentrations of target DNA without detection of non-target species, including other cyprinid fishes with which they co-occur.