Persistence at distributional edges: Columbia spotted frog habitat in the arid Great Basin,USA

A common challenge in the conservation of broadly distributed, yet imperiled species is understanding which factors facilitate persistence at distributional edges, locations where populations are often vulnerable to extirpation due to changes in climate, land use, or distributions of other species. For Columbia spotted frogs (Rana luteiventris) in the Great Basin (USA), a genetically distinct population segment of conservation concern, we approached this problem by examining (1) landscape‐scale habitat availability and distribution, (2) water body‐scale habitat associations, and (3) resource management‐identified threats to persistence. We found that areas with perennial aquatic habitat and suitable climate are extremely limited in the southern portion of the species’ range. Within these suitable areas, native and non‐native predators (trout and American bullfrogs [Lithobates catesbeianus]) are widespread and may further limit habitat availability in upper‐ and lower‐elevation areas, respectively. At the water body scale, spotted frog occupancy was associated with deeper sites containing abundant emergent vegetation and nontrout fish species. Streams with American beaver (Castor canadensis) frequently had these structural characteristics and were significantly more likely to be occupied than ponds, lakes, streams without beaver, or streams with inactive beaver ponds, highlighting the importance of active manipulation of stream environments by beaver. Native and non‐native trout reduced the likelihood of spotted frog occupancy, especially where emergent vegetation cover was sparse. Intensive livestock grazing, low aquatic connectivity, and ephemeral hydroperiods were also negatively associated with spotted frog occupancy. We conclude that persistence of this species at the arid end of its range has been largely facilitated by habitat stability (i.e., permanent hydroperiod), connectivity, predator‐free refugia, and a commensalistic interaction with an ecosystem engineer. Beaver‐induced changes to habitat quality, stability, and connectivity may increase spotted frog population resistance and resilience to seasonal drought, grazing, non‐native predators, and climate change, factors which threaten local or regional persistence.     

Hybridization between two parapatric ranid frog species in the northern Sierra Nevada,California, USA

Contact zones between species provide a unique opportunity to test whether taxa can hybridize or not. Cross‐breeding or hybridization between closely related taxa can promote gene flow (introgression) between species, adaptation, or even speciation. Though hybridization events may be short‐lived and difficult to detect in the field, genetic data can provide information about the level of introgression between closely related taxa. Hybridization can promote introgression between species, which may be an important evolutionary mechanism for either homogenization (reversing initial divergence between species) or reproductive isolation (potentially leading to speciation). Here, we used thousands of genetic markers from nuclear DNA to detect hybridization between two parapatric frog species (Rana boylii and Rana sierrae) in the Sierra Nevada of California. Based on principal components analysis, admixture, and analysis of heterozygosity at species diagnostic SNPs, we detected two F1 hybrid individuals in the Feather River basin, as well as a weak signal of introgression and gene flow between the frog species compared with frog populations from two other adjacent watersheds. This study provides the first documentation of hybridization and introgression between these two species, which are of conservation concern.     

Evaluating Sierra Nevada Yellow-Legged Frog Distribution Using Environmental DNA

Genetic tools that identify species from trace DNA samples could supplement traditional survey methods to clarify distributional limits of rare species. For species with legal habitat protection, elevational limits of distributions are used to determine where management actions may affect endangered species. The endangered Sierra Nevada yellow-legged frog (Rana sierrae) generally is found down to 1,370 m, but in the Plumas National Forest, California, USA, there are a number of historical records below this elevation, resulting in protections extending to 1,067 m. This species is phenotypically similar to the foothill yellow-legged frog (R. boylii), with which it occasionally hybridizes. We used a combination of genetic methods to investigate the fine-scale distribution of the Sierra Nevada yellow-legged frog in the Plumas National Forest. We collected and analyzed environmental DNA (eDNA) samples from all accessible lower elevation sites with records of Sierra Nevada yellow-legged frog (n = 17) and swabbed 220 individuals for genetic identification from 2016–2018 to clarify the distribution of this endangered species. We created a climatic suitability model using the validated Sierra Nevada yellow-legged frog records and current (1970–2000) climate models to assess additional highly suitable localities for Sierra Nevada yellow-legged frog presence using eDNA capture. We did not confirm detection of Sierra Nevada yellow-legged frog eDNA at any historical sites and identified all swabbed individuals from below 1,370 m (n = 144) as foothill yellow-legged frogs. We located a new Sierra Nevada yellow-legged frog site (at 1,919 m) during surveys guided by the climatic suitability model. It does not appear after extensive eDNA and genetic sampling that the Sierra Nevada yellow-legged frog occurs below 1,370 m in this portion of their range at present. Our results show that eDNA sampling can be used as an effective management tool to evaluate historical locations and previously unknown suitable localities for current presence of a species of interest. © 2021 The Authors. The Journal of Wildlife Management published by Wiley Periodicals LLC on behalf of The Wildlife Society.     

Environmental DNA metabarcoding of lake fish communities reflects long‐term data from established survey methods

Organisms continuously release DNA into their environments via shed cells, excreta, gametes and decaying material. Analysis of this ‘environmental DNA’ (eDNA) is revolutionizing biodiversity monitoring. eDNA outperforms many established survey methods for targeted detection of single species, but few studies have investigated how well eDNA reflects whole communities of organisms in natural environments. We investigated whether eDNA can recover accurate qualitative and quantitative information about fish communities in large lakes, by comparison to the most comprehensive long‐term gill‐net data set available in the UK. Seventy‐eight 2L water samples were collected along depth profile transects, gill‐net sites and from the shoreline in three large, deep lakes (Windermere, Bassenthwaite Lake and Derwent Water) in the English Lake District. Water samples were assayed by eDNA metabarcoding of the mitochondrial 12S and cytochrome b regions. Fourteen of the 16 species historically recorded in Windermere were detected using eDNA, compared to four species in the most recent gill‐net survey, demonstrating eDNA is extremely sensitive for detecting species. A key question for biodiversity monitoring is whether eDNA can accurately estimate abundance. To test this, we used the number of sequence reads per species and the proportion of sampling sites in which a species was detected with eDNA (i.e. site occupancy) as proxies for abundance. eDNA abundance data consistently correlated with rank abundance estimates from established surveys. These results demonstrate that eDNA metabarcoding can describe fish communities in large lakes, both qualitatively and quantitatively, and has great potential as a complementary tool to established monitoring methods.     

Environmental DNA detects a rare large river crayfish but with little relation to local abundance

相似文献   

Environmental DNA for the enumeration and management of Pacific salmon

Pacific salmon are a keystone resource in Alaska, generating annual revenues of well over ~US$500 million/year. Due to their anadromous life history, adult spawners distribute amongst thousands of streams, posing a huge management challenge. Currently, spawners are enumerated at just a few streams because of reliance on human counters and, rarely, sonar. The ability to detect organisms by shed tissue (environmental DNA, eDNA) promises a more efficient counting method. However, although eDNA correlates generally with local fish abundances, we do not know if eDNA can accurately enumerate salmon. Here we show that daily, and near‐daily, flow‐corrected eDNA rate closely tracks daily numbers of returning sockeye and coho spawners and outmigrating sockeye smolts. eDNA thus promises accurate and efficient enumeration, but to deliver the most robust numbers will need higher‐resolution stream‐flow data, at‐least‐daily sampling, and a focus on species with simple life histories, since shedding rate varies amongst jacks, juveniles, and adults.     

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.     

Predictive modeling for allopatric Strix (Strigiformes: Strigidae) owls in South America: determinants of their distributions and ecological niche‐based processes

Strix (Strigidae) is a worldwide genus of 17 owl species typical of forested habitats, including Rusty‐barred Owls (S. hylophila), Chaco Owls (S. chacoensis), and Rufous‐legged Owls (S. rufipes) in South America. These species are distributed allopatrically, but the ecological traits that determine their distributions remain largely unknown and their phylogenetic relationships are unclear. We used species distribution models (SDMs) to identify variables explaining their distribution patterns and test hypotheses about ecological divergence and conservatism based on niche overlap analysis. For Rusty‐barred Owls and Chaco Owls, climatic factors related to temperature played a major role, whereas a rainfall variable was more important for Rufous‐legged Owls. When niche overlaps were compared, accounting for regional similarities in the habitat available to each species, an ecological niche divergence process was supported for Chaco Owl‐Rusty‐barred Owl and Chaco Owl‐Rufous‐legged Owl, whereas a niche conservatism process was supported for Rusty‐barred Owl‐Rufous‐legged Owl. Different ecological requirements support current species delimitation, but they are in disagreement with the two main hypotheses currently envisaged about their phylogenetic relationships (Chaco Owls as the sister taxa of either Rufous‐legged Owls or Rusty‐barred Owls) and support a new phylogenetic hypothesis (Rufous‐legged Owls as sister taxa of Rusty‐barred Owls). Our findings suggest that speciation of Rusty‐barred Owls and Rufous‐legged Owls was a vicariant event resulting from Atlantic marine transgressions in southern South America in the Miocene, but their niche was conserved because habitat changed little in their respective ranges. In contrast, Chaco Owls diverged ecologically from the other two species as a result of their adaptations to the habitat they currently occupy. Ecological and historical approaches in biogeography can be embedded to explain distribution patterns, and results provided by SDMs can be used to infer historical and ecological processes in an integrative way.     

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.     

Microgeographic morphological variation across larval wood frog populations associated with environment despite gene flow

Gene flow has historically been thought to constrain local adaptation; yet, recent research suggests that populations can diverge despite exchanging genes. Here I use a common garden experiment to assess the combined effects of gene flow and natural selection on morphological variation of 16 wood frog (Rana sylvatica) populations, a species known to experience divergent selection pressures in open‐ and closed‐canopy ponds across relatively small geographic scales. Wood frog tadpoles from different ponds showed significant morphological variation associated with canopy type with a trade‐off between tail length and body depth consistent with previous research. In contrast, neutral genetic differentiation of nine microsatellite loci as measured by Jost's D was not associated with canopy type, indicating no pattern of isolation by environment. Genetic structure analyses indicated some substructure across the 16 ponds (K = 4); however, three out of four assigned clusters included both open‐ and closed‐canopy ponds. Together, these results suggest that morphological divergence among these wood frog populations is occurring despite gene flow and that selection within these environments is strong. Furthermore, morphological variation among ponds differed across two sampling periods during larval development, demonstrating the importance of evaluating phenotypic divergence over multiple time periods and at a time relevant to the processes being studied.     

Stepping into the past to conserve the future: Archived skin swabs from extant and extirpated populations inform genetic management of an endangered amphibian

Moving animals on a landscape through translocations and reintroductions is an important management tool used in the recovery of endangered species, particularly for the maintenance of population genetic diversity and structure. Management of imperiled amphibian species rely heavily on translocations and reintroductions, especially for species that have been brought to the brink of extinction by habitat loss, introduced species, and disease. One striking example of amphibian declines and associated management efforts is in California's Sequoia and Kings Canyon National Parks with the mountain yellow‐legged frog species complex (Rana sierrae/muscosa). Mountain yellow‐legged frogs have been extirpated from more than 93% of their historic range, and limited knowledge of their population genetics has made long‐term conservation planning difficult. To address this, we used 598 archived skin swabs from both extant and extirpated populations across 48 lake basins to generate a robust Illumina‐based nuclear amplicon data set. We found that samples grouped into three main genetic clusters, concordant with watershed boundaries. We also found evidence for historical gene flow across watershed boundaries with a north‐to‐south axis of migration. Finally, our results indicate that genetic diversity is not significantly different between populations with different disease histories. Our study offers specific management recommendations for imperiled mountain yellow‐legged frogs and, more broadly, provides a population genetic framework for leveraging minimally invasive samples for the conservation of threatened species.     

Moving eDNA surveys onto land: Strategies for active eDNA aggregation to detect invasive forest insects

The use of environmental DNA (eDNA) surveys to monitor terrestrial species has been relatively limited, with successful implementations still confined to sampling DNA from natural or artificial water bodies and soil. Sampling water for eDNA depends on proximity to or availability of water, whereas eDNA from soil is limited in its spatial scale due to the large quantities necessary for processing and difficulty in doing so. These challenges limit the widespread use of eDNA in several systems, such as surveying forests for invasive insects. We developed two new eDNA aggregation approaches that overcome the challenges of above‐ground terrestrial sampling and eliminate the dependency on creating or utilizing pre‐existing water bodies to conduct eDNA sampling. The first, “spray aggregation,” uses spray action to remove eDNA from surface substrates and was developed for shrubs and other understorey vegetation, while the second, “tree rolling,” uses physical transfer via a roller to remove eDNA from the surface of tree trunks and large branches. We tested these approaches by surveying for spotted lanternfly, Lycorma delicatula, a recent invasive pest of northeastern USA that is considered a significant ecological and economic threat to forests and agriculture. We found that our terrestrial eDNA surveys matched visual surveys, but also detected L. delicatula presence ahead of visual surveys, indicating increased sensitivity of terrestrial eDNA surveys over currently used methodology. The terrestrial eDNA approaches we describe can be adapted for use in surveying a variety of forest insects and represent a novel strategy for surveying terrestrial biodiversity.     

Environmental DNA facilitates accurate,inexpensive, and multiyear population estimates of millions of anadromous fish

Although environmental DNA shed from an organism is now widely used for species detection in a wide variety of contexts, mobilizing environmental DNA for management requires estimation of population size and trends in addition to assessing presence or absence. However, the efficacy of environmental‐DNA‐based indices of abundance for long‐term population monitoring have not yet been assessed. Here we report on the relationship between six years of mark‐recapture population estimates for eulachon (Thaleichthys pacificus) and “eDNA rates” which are calculated from the product of stream flow and DNA concentration. Eulachon are a culturally and biologically important anadromous fish that have significantly declined in the southern part of their range but were historically rendered into oil and traded. Both the peak eDNA rate and the area under the curve of the daily eDNA rate were highly predictive of the mark‐recapture population estimate, explaining 84.96% and 92.53% of the deviance, respectively. Even in the absence of flow correction, the peak of the daily eDNA concentration explained an astonishing 89.53% while the area under the curve explained 90.74% of the deviance. These results support the use of eDNA to monitor eulachon population trends and represent a >80% cost savings over mark‐recapture, which could be further increased with automated water sampling, reduced replication, and focused temporal sampling. Due to its logistical ease and affordability, eDNA sampling can facilitate monitoring a larger number of rivers and in remote locations where mark‐recapture is infeasible.     

Microsatellite DNA markers for the Chinese wood frog (Rana chensinensis) and tests for their cross-utility in 15 ranid frog species

We developed 22 microsatellite markers for the Chinese wood frog (Rana chensinensis) to study the impact of landscape features on its population structure. Thirty‐four individuals from one breeding site were examined and 14 loci were polymorphic. The number of alleles, expected heterozygosity and observed heterozygosity varied from two to 14, from 0.0833 to 0.9118, and from 0.1376 to 0.8667, respectively. Cross‐species amplification was tested for 15 ranid frog species. The Plateau brown frog, Rana kukunoris (n = 23), was successfully amplified at 18 loci, and 15 were polymorphic with number of alleles varying from two to 18. Ten other species were also amplified at a limited number of loci.     

Effects of changing climate on aquatic habitat and connectivity for remnant populations of a wide‐ranging frog species in an arid landscape

Amphibian species persisting in isolated streams and wetlands in desert environments can be susceptible to low connectivity, genetic isolation, and climate changes. We evaluated the past (1900–1930), recent (1981–2010), and future (2071–2100) climate suitability of the arid Great Basin (USA) for the Columbia spotted frog (Rana luteiventris) and assessed whether changes in surface water may affect connectivity for remaining populations. We developed a predictive model of current climate suitability and used it to predict the historic and future distribution of suitable climates. We then modeled changes in surface water availability at each time period. Finally, we quantified connectivity among existing populations on the basis of hydrology and correlated it with interpopulation genetic distance. We found that the area of the Great Basin with suitable climate conditions has declined by approximately 49% over the last century and will likely continue to decline under future climate scenarios. Climate conditions at currently occupied locations have been relatively stable over the last century, which may explain persistence at these sites. However, future climates at these currently occupied locations are predicted to become warmer throughout the year and drier during the frog's activity period (May – September). Fall and winter precipitation may increase, but as rain instead of snow. Earlier runoff and lower summer base flows may reduce connectivity between neighboring populations, which is already limited. Many of these changes could have negative effects on remaining populations over the next 50–80 years, but milder winters, longer growing seasons, and wetter falls might positively affect survival and dispersal. Collectively, however, seasonal shifts in temperature, precipitation, and stream flow patterns could reduce habitat suitability and connectivity for frogs and possibly other aquatic species inhabiting streams in this arid region.     

Hybridization and mitochondrial genome introgression between Rana chensinensis and R. kukunoris

Mitochondrial genome (mito‐genome) introgression among metazoans is commonplace, and several biological processes may promote such introgression. We examined two proposed processes for the mito‐genome introgression between Rana chensinensis and R. kukunoris: natural hybridization and sex‐biased dispersal. We sampled 477 individuals from 28 sites in the potential hybrid zone in the western Tsinling Mountains. Mitochondrial gene (cyt‐b) trees were used to examine the introgression events. Microsatellite DNA loci, cyt‐b and morphological data were used to identify hybrids and to examine the extent of natural hybridization. We detected rampant bidirectional introgressions, both ancient and recent, between the two species. Furthermore, we found a wide hybrid zone, and frequent and asymmetric hybridization. The hybrid zone cline analysis revealed a clear mitochondrial–nuclear discordance; while most nuclear markers displayed similar and steep clines, cyt‐b had a displaced cline centre and a more gradual and wider cline. We also detected strong and asymmetric historical maternal gene flow across the hybrid zone. This widespread hybridization and detected low mito‐nuclear conflicts may, at least partially, explain the high frequency of introgression. Lastly, microsatellite data and population genetic methods were used to assess sex‐biased dispersal. A weak pattern of female‐biased dispersal was detected in both species, suggesting it may not play an important role in the observed introgression. Our data are consistent with the hybridization hypothesis, but support for the sex‐biased dispersal hypothesis is weak. We further suggest that selective advantages of the R. kukunoris‐type mito‐genome in thermal adaptation may also contribute to the introgression between the two species.     

An exception to the matched filter hypothesis: A mismatch of male call frequency and female best hearing frequency in a torrent frog

The matched filter hypothesis proposes that the tuning of auditory sensitivity and the spectral character of calls will match in order to maximize auditory processing efficiency during courtship. In this study, we analyzed the acoustic structure of male calls and both male and female hearing sensitivities in the little torrent frog (Amolops torrentis), an anuran species who transmits acoustic signals across streams. The results were in striking contradiction to the matched filter hypothesis. Auditory brainstem response results showed that the best hearing range was 1.6–2 kHz consistent with the best sensitive frequency of most terrestrial lentic taxa, yet completely mismatched with the dominant frequency of conspecific calls (4.3 kHz). Moreover, phonotaxis tests show that females strongly prefer high‐frequency (4.3 kHz) over low‐frequency calls (1.6 kHz) regardless of ambient noise levels, although peripheral auditory sensitivity is highest in the 1.6–2 kHz range. These results are consistent with the idea that A. torrentis evolved from nonstreamside species and that high‐frequency calls evolved under the pressure of stream noise. Our results also suggest that female preferences based on central auditory system characteristics may evolve independently of peripheral auditory system sensitivity in order to maximize communication effectiveness in noisy environments.     

Optimizing techniques to capture and extract environmental DNA for detection and quantification of fish

Few studies have examined capture and extraction methods for environmental DNA (eDNA) to identify techniques optimal for detection and quantification. In this study, precipitation, centrifugation and filtration eDNA capture methods and six commercially available DNA extraction kits were evaluated for their ability to detect and quantify common carp (Cyprinus carpio) mitochondrial DNA using quantitative PCR in a series of laboratory experiments. Filtration methods yielded the most carp eDNA, and a glass fibre (GF) filter performed better than a similar pore size polycarbonate (PC) filter. Smaller pore sized filters had higher regression slopes of biomass to eDNA, indicating that they were potentially more sensitive to changes in biomass. Comparison of DNA extraction kits showed that the MP Biomedicals FastDNA SPIN Kit yielded the most carp eDNA and was the most sensitive for detection purposes, despite minor inhibition. The MoBio PowerSoil DNA Isolation Kit had the lowest coefficient of variation in extraction efficiency between lake and well water and had no detectable inhibition, making it most suitable for comparisons across aquatic environments. Of the methods tested, we recommend using a 1.5 μm GF filter, followed by extraction with the MP Biomedicals FastDNA SPIN Kit for detection. For quantification of eDNA, filtration through a 0.2–0.6 μm pore size PC filter, followed by extraction with MoBio PowerSoil DNA Isolation Kit was optimal. These results are broadly applicable for laboratory studies on carps and potentially other cyprinids. The recommendations can also be used to inform choice of methodology for field studies.     

A probabilistic model for designing and assessing the performance of eDNA sampling protocols

Environmental DNA (eDNA) sampling, the detection of species‐specific genetic material in water samples, is an emerging tool for monitoring aquatic invasive species. Optimizing eDNA sampling protocols can be challenging because there is imperfect understanding of how each step of the protocol influences its sensitivity. This paper develops a probabilistic model that characterizes each step of an eDNA sampling protocol to evaluate the protocol's overall detection sensitivity for one sample. The model is then applied to analyse how changes over time made to the eDNA sampling protocol to detect bighead (BH) and silver carp (SC) eDNA have influenced its sensitivity, and hence interpretation of the results. The model shows that changes to the protocol have caused the sensitivity of the protocol to fluctuate. A more efficient extraction method in 2013, new species‐specific markers with a qPCR assay in 2014, and a more efficient capture method in 2015 have improved the sensitivity, while switching to a larger elution volume in 2013 and a smaller sample volume in 2015 have reduced the sensitivity. Overall, the sensitivity of the current protocol is higher for BH eDNA detection and SC eDNA detection compared to the original protocol used from 2009 to 2012. The paper shows how this model of eDNA sampling can be used to evaluate the effect of proposed changes in an eDNA sampling and analysis protocol on the sensitivity of that protocol to help researchers optimize their design.