Using occupancy modelling to compare environmental DNA to traditional field methods for regional‐scale monitoring of an endangered aquatic species

Environmental DNA (eDNA) monitoring approaches promise to greatly improve detection of rare, endangered and invasive species in comparison with traditional field approaches. Herein, eDNA approaches and traditional seining methods were applied at 29 research locations to compare method‐specific estimates of detection and occupancy probabilities for endangered tidewater goby (Eucyclogobius newberryi). At each location, multiple paired seine hauls and water samples for eDNA analysis were taken, ranging from two to 23 samples per site, depending upon habitat size. Analysis using a multimethod occupancy modelling framework indicated that the probability of detection using eDNA was nearly double (0.74) the rate of detection for seining (0.39). The higher detection rates afforded by eDNA allowed determination of tidewater goby occupancy at two locations where they have not been previously detected and at one location considered to be locally extirpated. Additionally, eDNA concentration was positively related to tidewater goby catch per unit effort, suggesting eDNA could potentially be used as a proxy for local tidewater goby abundance. Compared to traditional field sampling, eDNA provided improved occupancy parameter estimates and can be applied to increase management efficiency across a broad spatial range and within a diversity of habitats.     

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.     

Environmental DNA analysis for estimating the abundance and biomass of stream fish

相似文献   

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.     

Seasonal changes in fish catch and environmental variables in a large Tropical Lake,Volta, Ghana

Fish species assemblage and selected environmental variables were monitored in Lake Volta from September 2014 to August 2016 to determine seasonal variability in species composition, catch and environmental variables that determine the structure of the fish community. A total of 1,557 individual fish belonging to 41 species, 25 genera and 13 families were recorded. The important fish species with respect to frequency of occurrence, abundance and weight, respectively, were as follows: Chrysichthys (100%; 43.03%; 17.93%), tilapias (100%; 28.97%; 17.86%), Alestes (100%; 14.13%; 32.10%) and Bagrus (91.67%; 5.65%; 12.80%) in that order. The composition of fish species and their temporal variation in experimental catches were similar to that of commercial catches. There were no significant differences (p > 0.05) in species abundance, weight and diversity indices between the dry and wet seasons. The modal class of length frequency distribution of the dominant species in the catch, Chrysichthys, reduced from 125 mm in 2006 to 95 mm in the current study indicating overfishing. Environmental variables considered showed little variation and within optimal ranges for fish survival and did not differ significantly between seasons. Canonical Correspondence Analysis showed that environmental variables explained 43.30% of the variation in species abundance with Lake water level, nitrite‐nitrogen and total dissolved solids being the main environmental factors influencing the structure of the fish community.     

Combining genetic and demographic data for prioritizing conservation actions: insights from a threatened fish species

Prioritizing and making efficient conservation plans for threatened populations requires information at both evolutionary and ecological timescales. Nevertheless, few studies integrate multidisciplinary approaches, mainly because of the difficulty for conservationists to assess simultaneously the evolutionary and ecological status of populations. Here, we sought to demonstrate how combining genetic and demographic analyses allows prioritizing and initiating conservation plans. To do so, we combined snapshot microsatellite data and a 30‐year‐long demographic survey on a threatened freshwater fish species (Parachondrostoma toxostoma) at the river basin scale. Our results revealed low levels of genetic diversity and weak effective population sizes (<63 individuals) in all populations. We further detected severe bottlenecks dating back to the last centuries (200–800 years ago), which may explain the differentiation of certain populations. The demographic survey revealed a general decrease in the spatial distribution and abundance of P. toxostoma over the last three decades. We conclude that demo‐genetic approaches are essential for (1) identifying populations for which both evolutionary and ecological extinction risks are high; and (2) proposing conservation plans targeted toward these at risk populations, and accounting for the evolutionary history of populations. We suggest that demo‐genetic approaches should be the norm in conservation practices.     

Reliable eDNA detection and quantification of the European weather loach (Misgurnus fossilis)

The European weather loach (Misgurnus fossilis) is a cryptic and poorly known fish species of high conservation concern. The species is experiencing dramatic population collapses across its native range to the point of regional extinction. Although environmental DNA (eDNA)-based approaches offer clear advantages over conventional field methods for monitoring rare and endangered species, accurate detection and quantification remain difficult and quality assessment is often poorly incorporated. In this study, we developed and validated a novel digital droplet PCR (ddPCR) eDNA-based method for reliable detection and quantification, which allows accurate monitoring of M. fossilis across a number of habitat types. A dilution experiment under laboratory conditions allowed the definition of the limit of detection (LOD) and the limit of quantification (LOQ), which were set at concentrations of 0.07 and 0.14 copies μl^–1, respectively. A series of aquarium experiments revealed a significant and positive relationship between the number of individuals and the eDNA concentration measured. During a 3 year survey (2017–2019), we assessed 96 locations for the presence of M. fossilis in Flanders (Belgium). eDNA analyses on these samples highlighted 45% positive detections of the species. On the basis of the eDNA concentration per litre of water, only 12 sites appeared to harbour relatively dense populations. The other 31 sites gave a relatively weak positive signal that was typically situated below the LOQ. Combining sample-specific estimates of effective DNA quantity (Q_e) and conventional field sampling, we concluded that each of these weak positive sites still likely harboured the species and therefore they do not represent false positives. Further, only seven of the classified negative samples warrant additional sampling as our analyses identified a substantial risk of false-negative detections (i.e., type II errors) at these locations. Finally, we illustrated that ddPCR outcompetes conventional qPCR analyses, especially when target DNA concentrations are critically low, which could be attributed to a reduced sensitivity of ddPCR to inhibition effects, higher sample concentrations being accommodated and higher sensitivity obtained.     

Microsatellites from the vairone Leuciscus souffia (Pisces: Cyprinidae) and their application to closely related species

Eleven microsatellites were isolated from the vairone Leuciscus souffia (Risso 1826), an endangered fish that inhabits river systems in and around the Alps in Europe. The level of genetic diversity was assessed in 29 individuals of the subspecies L. s. souffia, and their variability was further estimated in seven individuals of a different subspecies, L. s. muticellus. Eight of these microsatellite loci were also applied to seven closely related cyprinid species. Availability of the reported microsatellite loci will facilitate the investigation of population genetic structure of these species with applications for the development of conservation strategies and phylogeographical approaches.     

Sequence capture and next‐generation sequencing of ultraconserved elements in a large‐genome salamander

Amidst the rapid advancement in next‐generation sequencing (NGS) technology over the last few years, salamanders have been left behind. Salamanders have enormous genomes—up to 40 times the size of the human genome—and this poses challenges to generating NGS data sets of quality and quantity similar to those of other vertebrates. However, optimization of laboratory protocols is time‐consuming and often cost prohibitive, and continued omission of salamanders from novel phylogeographic research is detrimental to species facing decline. Here, we use a salamander endemic to the southeastern United States, Plethodon serratus, to test the utility of an established protocol for sequence capture of ultraconserved elements (UCEs) in resolving intraspecific phylogeographic relationships and delimiting cryptic species. Without modifying the standard laboratory protocol, we generated a data set consisting of over 600 million reads for 85 P. serratus samples. Species delimitation analyses support recognition of seven species within P. serratus sensu lato, and all phylogenetic relationships among the seven species are fully resolved under a coalescent model. Results also corroborate previous data suggesting nonmonophyly of the Ouachita and Louisiana regions. Our results demonstrate that established UCE protocols can successfully be used in phylogeographic studies of salamander species, providing a powerful tool for future research on evolutionary history of amphibians and other organisms with large genomes.     

Genetic variation in blue whales in the eastern pacific: implication for taxonomy and use of common wintering grounds

Many aspects of blue whale biology are poorly understood. Some of the gaps in our knowledge, such as those regarding their basic taxonomy and seasonal movements, directly affect our ability to monitor and manage blue whale populations. As a step towards filling in some of these gaps, microsatellite and mtDNA sequence analyses were conducted on blue whale samples from the Southern Hemisphere, the eastern tropical Pacific (ETP) and the northeast Pacific. The results indicate that the ETP is differentially used by blue whales from the northern and southern eastern Pacific, with the former showing stronger affinity to the region off Central America known as the Costa Rican Dome, and the latter favouring the waters of Peru and Ecuador. Although the pattern of genetic variation throughout the Southern Hemisphere is compatible with the recently proposed subspecies status of Chilean blue whales, some discrepancies remain between catch lengths and lengths from aerial photography, and not all blue whales in Chilean waters can be assumed to be of this type. Also, the range of the proposed Chilean subspecies, which extends to the Galapagos region of the ETP, at least seasonally, perhaps should include the Costa Rican Dome and the eastern North Pacific as well.     

Morphometric and genetic differentiation among populations of flat‐headed cusimanse (Crossarchus platycephalus) in Nigeria

Geographic barriers can partition genetic diversity among populations and drive evolutionary divergence between populations, promoting the speciation process and affecting conservation goals. We integrated morphological and genomic data to assess the distribution of variation in the flat‐headed cusimanse (Crossarchus platycephalus), a species of least conservation concern, on either side of the River Niger in Nigeria. Ecological disturbances affect the conservation status of many other animals in this region. The two populations were differentiated in the snout and fore limbs, with greater morphological diversity in the western population. We used Restriction site Associated DNA sequencing (RAD‐seq) and identified two genotypic clusters in a STRUCTURE analysis. Individuals from the eastern population are almost entirely assigned to one cluster, whereas genotypes from the western population are a mixture of the two clusters. The population from west of the River Niger also had higher heterozygosity. The morphological and population genetic data are therefore in agreement that the population from west of the River Niger is more diverse than the eastern population, and the eastern population contains a subset of the genetic variation found in the western population. Our results demonstrate that combining morphological and genotypic measures of diversity can provide a congruent picture of the distribution of intraspecific variation. The results also suggest that future work should explore the role of the River Niger as a natural barrier to migration in Nigeria.     

Detecting selection along environmental gradients: analysis of eight methods and their effectiveness for outbreeding and selfing populations

Thanks to genome‐scale diversity data, present‐day studies can provide a detailed view of how natural and cultivated species adapt to their environment and particularly to environmental gradients. However, due to their sensitivity, up‐to‐date studies might be more sensitive to undocumented demographic effects such as the pattern of migration and the reproduction regime. In this study, we provide guidelines for the use of popular or recently developed statistical methods to detect footprints of selection. We simulated 100 populations along a selective gradient and explored different migration models, sampling schemes and rates of self‐fertilization. We investigated the power and robustness of eight methods to detect loci potentially under selection: three designed to detect genotype–environment correlations and five designed to detect adaptive differentiation (based on F_ST or similar measures). We show that genotype–environment correlation methods have substantially more power to detect selection than differentiation‐based methods but that they generally suffer from high rates of false positives. This effect is exacerbated whenever allele frequencies are correlated, either between populations or within populations. Our results suggest that, when the underlying genetic structure of the data is unknown, a number of robust methods are preferable. Moreover, in the simulated scenario we used, sampling many populations led to better results than sampling many individuals per population. Finally, care should be taken when using methods to identify genotype–environment correlations without correcting for allele frequency autocorrelation because of the risk of spurious signals due to allele frequency correlations between populations.     

Blue whale population structure along the eastern South Pacific Ocean: evidence of more than one population

Blue whales (Balaenoptera musculus) were among the most intensively exploited species of whales in the world. As a consequence of this intense exploitation, blue whale sightings off the coast of Chile were uncommon by the end of the 20th century. In 2004, a feeding and nursing ground was reported in southern Chile (SCh). With the aim to investigate the genetic identity and relationship of these Chilean blue whales to those in other Southern Hemisphere areas, 60 biopsy samples were collected from blue whales in SCh between 2003 and 2009. These samples were genotyped at seven microsatellite loci and the mitochondrial control region was sequenced, allowing us to identify 52 individuals. To investigate the genetic identity of this suspected remnant population, we compared these 52 individuals to blue whales from Antarctica (ANT, n = 96), Northern Chile (NCh, n = 19) and the eastern tropical Pacific (ETP, n = 31). No significant differentiation in haplotype frequencies (mtDNA) or among genotypes (nDNA) was found between SCh, NCh and ETP, while significant differences were found between those three areas and Antarctica for both the mitochondrial and microsatellite analyses. Our results suggest at least two breeding population units or subspecies exist, which is also supported by other lines of evidence such as morphometrics and acoustics. The lack of differences detected between SCh/NCh/ETP areas supports the hypothesis that eastern South Pacific blue whales are using the ETP area as a possible breeding area. Considering the small population sizes previously reported for the SCh area, additional conservation measures and monitoring of this population should be developed and prioritized.     

Consequences of temperature and temperature variability on swimming activity,group structure,and predation of endangered delta smelt

相似文献   

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.     

Molecular monitoring of Escherichia coli O157: H7 sterilization rate using qPCR and propidium monoazide treatment

Propidium monoazide is a DNA‐intercalating dye. PMA‐qPCR has been reported as a novel method to detect live bacteria in complex samples. In this study, this method was used to monitor the sterilization effects of UHP, ultrasound and high PEF on Escherichia coli O157:H7. Our results showed that all three sterilization techniques are successful to kill viable E. coli O157:H7 cells under their appropriate conditions. PMA‐qPCR can effectively monitor the amount of DNA released from viable E. coli O157:H7 cells, and the results from PMA‐qPCR were highly consistent with those from plate counting after treatment with UHP, ultrasound and high PEF. The maximal ΔC_t between PMA‐qPCR and qPCR obtained in this study was 10·39 for UHP, 5·76 for ultrasound and 2·30 for high PEF. The maximal sterilization rates monitored by PMA‐qPCR were 99·92% for UHP, 99·99% for ultrasound and 100% for high PEF. Thus, PMA‐qPCR can be used to detect the sterilization effect on food and water supplies after treatment with UHP, ultrasound and high PEF.

Significance and Impact of the Study

The reliable detection of viable foodborne pathogenic bacteria in water and food is of great importance in our daily life. However, the traditional bacteria cultivation‐based methods are time‐consuming and difficult to monitor all viable bacteria because of the limitation of cultivation conditions. This study demonstrated that PMA‐qPCR technique is very effective to monitor viable E. coli O157:H7 after sterilization and will help to monitor the viable bacteria in food and water.     

Metazoan parasites of African annual killifish (Nothobranchiidae): abundance,diversity, and their environmental correlates

Estimates of biodiversity and its global patterns are affected by parasite richness and specificity. Despite this, parasite communities are largely neglected in biodiversity estimates, especially in the tropics. We studied the parasites of annual killifish of the genus Nothobranchius that inhabit annually desiccating pools across the African savannah and survive the dry period as developmentally arrested embryos. Their discontinuous, non‐overlapping generations make them a unique organism in which to study natural parasite fauna. We investigated the relationship between global (climate and altitude) and local (pool size, vegetation, host density and diversity, and diversity of potential intermediate hosts) environmental factors and the community structure of killifish parasites. We examined metazoan parasites from 21 populations of four host species (Nothobranchius orthonotus, N. furzeri, N. kadleci, and N. pienaari) across a gradient of aridity in Mozambique. Seventeen parasite taxa were recorded, with trematode larval stages (metacercariae) being the most abundant taxa. The parasites recorded were both allogenic (life cycle includes non‐aquatic host; predominantly trematodes) and autogenic (cycling only in aquatic hosts; nematodes). The parasite abundance was highest in climatic regions with intermediate aridity, while parasite diversity was associated with local environmental characteristics and positively correlated with fish species diversity and the amount of aquatic vegetation. Our results suggest that parasite communities of sympatric Nothobranchius species are similar and dominated by the larval stages of generalist parasites. Therefore, Nothobranchius serve as important intermediate or paratenic hosts of parasites, with piscivorous birds and predatory fish being their most likely definitive hosts.     

Endemic species may have complex histories: within‐refugium phylogeography of an endangered Iberian vole

Glacial refugia protected and promoted biodiversity during the Pleistocene, not only at a broader scale, but also for many endemics that contracted and expanded their ranges within refugial areas. Understanding the evolutionary history of refugial endemics is especially important in the case of endangered species to recognize the origins of their genetic structure and thus produce better informed conservation practices. The Iberian Peninsula is an important European glacial refugium, rich in endemics of conservation concern, including small mammals, such as the Cabrera vole (Microtus cabrerae). This near‐threatened rodent is characterized by an unusual suite of genetic, life history and ecological traits, being restricted to isolated geographic nuclei in fast‐disappearing Mediterranean subhumid herbaceous habitats. To reconstruct the evolutionary history of the Cabrera vole, we studied sequence variation at mitochondrial, autosomal and sex‐linked loci, using invasive and noninvasive samples. Despite low overall mitochondrial and nuclear nucleotide diversities, we observed two main well‐supported mitochondrial lineages, west and east. Phylogeographic modelling in the context of the Cabrera vole's detailed fossil record supports a demographic scenario of isolation of two populations during the Last Glacial Maximum from a single focus in the southern part of the Iberian Peninsula. In addition, our data suggest subsequent divergence within the east, and secondary contact and introgression of the expanding western population, during the late Holocene. This work emphasizes that refugial endemics may have a phylogeographic history as rich as that of more widespread species, and conservation of such endemics includes the preservation of that genetic legacy.     

Invasion genetics from eDNA and thousands of larvae: A targeted metabarcoding assay that distinguishes species and population variation of zebra and quagga mussels

Identifying species and population genetic compositions of biological invasions at early life stages and/or from environmental (e)DNA using targeted high‐throughput sequencing (HTS) metabarcode assays offers powerful and cost‐effective means for early detection, analysis of spread patterns, and evaluating population changes. The present study develops, tests, and applies this method with a targeted sequence assay designed to simultaneously identify and distinguish between the closely related invasive Eurasian zebra and quagga mussels (Dreissena polymorpha and D. rostriformis) and their relatives and discern their respective population genetic patterns. Invasions of these dreissenid mussel species have markedly changed freshwater ecosystems throughout North America and Europe, exerting severe ecological and economic damage. Their planktonic early life stages (eggs and larvae) are morphologically indistinguishable, yet each species exerts differential ecological effects, with the quagga often outcompeting the zebra mussel as adults. Our targeted assay analyzes genetic variation from a diagnostic sequence region of the mitochondrial (mt)DNA cytochrome oxidase I (COI) gene, to assess temporal and spatial inter‐ and intra‐specific genetic variability. The assay facilitates analysis of environmental (e)DNA from water, early life stages from thousands of individuals, and simultaneous analysis of 50–100 tagged field‐collected samples. Experiments evaluated its accuracy and performance using: (a) mock laboratory communities containing known DNA quantities per taxon, (b) aquaria with mixed‐species/haplotype compositions of adults, and (c) field‐collected water and plankton versus traditional sampling of adult communities. Results delineated species compositions, relative abundances, and population‐level diversity differences among ecosystems, habitats, time series, and life stages from two allopatric concurrent invasions in the Great Lakes (Lake Erie) and the Hudson River, which had separate founding histories. Findings demonstrate application of this targeted assay and our approach to accurately and simultaneously discern species‐ and population‐level differences across spatial and temporal scales, facilitating early detection and ecological understanding of biological invasions.