The future of biotic indices in the ecogenomic era: Integrating (e)DNA metabarcoding in biological assessment of aquatic ecosystems

Using DNA metabarcoding to characterize national scale diatom-environment relationships and to develop indicators in streams and rivers of the United States

Recent advancements in DNA techniques, metabarcoding, and bioinformatics could help expand the use of benthic diatoms in monitoring and assessment programs by providing relatively quick and increasingly cost-effective ways to quantify diatom diversity in environmental samples. However, such applications of DNA-based approaches are relatively new, and in the United States, unknowns regarding their applications at large scales exist because only a few small-scale studies have been done. Here, we present results from the first nationwide survey to use DNA metabarcoding (rbcL) of benthic diatoms, which were collected from 1788 streams and rivers across nine ecoregions spanning the conterminous USA. At the national scale, we found that diatom assemblage structure (1) was strongly associated with total phosphorus and total nitrogen concentrations, conductivity, and pH and (2) had clear patterns that corresponded with differences in these variables among the nine ecoregions. These four variables were strong predictors of diatom assemblage structure in ecoregion-specific analyses, but our results also showed that diatom-environment relationships, the importance of environmental variables, and the ranges of these variables within which assemblage changes occurred differed among ecoregions. To further examine how assemblage data could be used for biomonitoring purposes, we used indicator species analysis to identify ecoregion-specific taxa that decreased or increased along each environmental gradient, and we used their relative abundances of gene reads in samples as metrics. These metrics were strongly correlated with their corresponding variable of interest (e.g., low phosphorus diatoms with total phosphorus concentrations), and generalized additive models showed how their relationships compared among ecoregions. These large-scale national patterns and nine sets of ecoregional results demonstrated that diatom DNA metabarcoding is a robust approach that could be useful to monitoring and assessment programs spanning the variety of conditions that exist throughout the conterminous United States.

Taxonomic accuracy and complementarity between bulk and eDNA metabarcoding provides an alternative to morphology for biological assessment of freshwater macroinvertebrates

Determining biological status of freshwater ecosystems is critical for ensuring ecosystem health and maintaining associated services to such ecosystems. Freshwater macroinvertebrates respond predictably to environmental disturbances and are widely used in biomonitoring programs. However, many freshwater species are difficult to capture and sort from debris or substrate and morphological identification is challenging, especially larval stages, damaged specimens, or hyperdiverse groups such as Diptera. The advent of high throughput sequencing technologies has enhanced DNA barcoding tools to automatise species identification for whole communities, as metabarcoding is increasingly used to monitor biodiversity. However, recent comparisons have revealed little congruence between morphological and molecular-based identifications. Using broad range universal primers for DNA barcode marker cox1, we compare community composition captured between morphological and molecular-based approaches from different sources - tissue-based (bulk benthic and bulk drift samples) and environmental DNA (eDNA, filtered water) metabarcoding - for samples collected along a gradient of anthropogenic disturbances. For comparability, metabarcoding taxonomic assignments were filtered by taxa included in the standardised national biological metric IBMWP. At the family level, bulk benthic metabarcoding showed the highest congruence with morphology, and the most abundant taxa were captured by all techniques. Richness captured by morphology and bulk benthic metabarcoding decreased along the gradient, whereas richness recorded by eDNA remained constant and increased downstream when sequencing bulk drift. Estimates of biological metrics were higher using molecular than morphological identification. At species level, diversity captured by bulk benthic samples were higher than the other techniques. Importantly, bulk benthic and eDNA metabarcoding captured different and complementary portions of the community - benthic versus water column, respectively - and their combined use is recommended. While bulk benthic metabarcoding can likely replace morphology using similar benthic biological indices, water eDNA will require new metrics because this technique sequences a different portion of the community.

Assessing the potential of nematode metabarcoding for benthic monitoring of offshore oil platforms

Environmental DNA metabarcoding is gaining momentum as a time and cost-effective tool for biomonitoring and environmental impact assessment. Yet, its use as a replacement for the conventional marine benthic monitoring based on morphological analysis of macrofauna is still challenging. Here we propose to study the meiofauna, which is much better represented in sediment DNA samples. We focus on nematodes, which are the most numerous and diverse group of meiofauna. Our aim is to assess the potential of nematode metabarcoding to monitor impacts associated with offshore oil platform activities. To achieve this goal, we used nematode-optimized marker (18S V1V2-Nema) and universal eukaryotic marker (18S V9) region to analyse 252 sediment DNA samples collected near three offshore oil platforms in the North Sea. For both markers, we analysed changes in alpha and beta diversity in relation to distance from the platforms and environmental variables. We also defined three impact classes based on selected environmental variables that are associated with oil extraction activities and used random forest classifiers to compare the predictive performance of both datasets. Our results show that alpha- and beta-diversity of nematodes varies with the increasing distance from the platforms. The variables directly related to platform activity, such as Ba and THC, strongly influence the nematode community. The nematode metabarcoding data provide more robust predictive models than eukaryotic data. Furthermore, the nematode community appears more stable in time and space, as illustrated by the overlap of nematode datasets obtained from the same platform three years apart. A significative negative correlation between distance and Shannon diversity also advocates for higher performance of the V1V2-Nema over the V9. Overall, these results suggest that the sensitivity of nematodes is higher compared to the eukaryotic community. Hence, nematode metabarcoding has the potential to become an effective tool for benthic monitoring in marine environment.

Environmental heterogeneity caused by large-scale cultivation of Pyropia haitanensis shapes multi-group biodiversity distribution in coastal areas

The response of marine biodiversity to mariculture has long been a research focus in marine ecology. However, the effects of seaweed cultivation on biological community assembly are poorly understood, especially in diverse communities with distinct ecological characteristics. In this study, we used environmental DNA metabarcoding to investigate the spatial distribution patterns of bacterial, protistan, and metazoan diversity, aiming to reveal the mechanisms of community assembly in the Pyropia haitanensis cultivation zone along the Fujian coast, China. We found that, compared with the biological communities in control zones, those in P. haitanensis cultivation zones exhibited stronger geographic distance-decay patterns and displayed more complex and stable network structures. Deterministic processes (environmental selection) played a more important role in the assembly of bacterial, protistan, and metazoan communities in P. haitanensis cultivation zones, especially metazoan communities. Variance partitioning analysis showed that environmental variables made greater contributions to the diversity of the three types of communities within the P. haitanensis cultivation zones than in the control zones. Partial least squares path modeling analysis identified nitrate‑nitrogen (NO3-N), pH, particulate organic carbon (POC), and dissolved organic carbon (DOC) as the key environmental variables affecting biodiversity. Overall, the environmental heterogeneity caused by the large-scale cultivation of P. haitanensis could be the crucial factor influencing the composition and structure of various biological communities. Our results highlight the importance of the responses of multi-group organisms to the cultivation of seaweed, and provide insights into the coexistence patterns of biodiversity at the spatial scale.

Ciliate diversity in rodrigo de freitas lagoon (Rio de Janeiro, Brazil) from an integrative standpoint

The Rodrigo de Freitas Lagoon is a highly eutrophic lacustrine system and has one of the longest histories of exploration and anthropic alteration in Brazil. Despite its relevance, limited studies explored the diversity of micro-eukaryotes in the lagoon. Ciliates (Alveolata, Ciliophora) are overlooked in environmental microbiology, especially in tropical and subtropical ecosystems, resulting in limited knowledge about their diversity and functional relevance in South American habitats, particularly in coastal lagoons. To fill this gap, here we investigated the diversity of ciliates in a brackish coastal lagoon in an urban area of Rio de Janeiro, Brazil, applying and comparing the performance of morphological and metabarcoding approaches. The metabarcoding analysis, based on high-throughput sequencing of the hipervariable region V4 of the 18S rRNA genes detected 37 molecular operational taxonomic units (MOTUs) assigned to Ciliophora, representing only about a half (56.9%) of the diversity detected by microscopy, which counted 65 ciliate morphotypes. The most representative classes in both approaches were Spirotrichea and Oligohymenophorea. The metabarcoding analysis revealed that 35.3% of the ciliate MOTUs had less than 97% similarity to available sequences in the NCBI database, indicating that more than one-third of these MOTUs potentially represents still not represented or undescribed ciliate species in current databases. Our findings indicate that metabarcoding techniques can significantly enhance the comprehension of ciliate diversity in tropical environments, but the scarcity of reference sequences of brackish ciliates in molecular databases represents a challenge to the taxonomic assignment of the MOTUs. This study provides new insights into the diversity of ciliates in a threatened coastal lagoon, revealing a vast array of still unknown and rare ciliate taxonomic units in tropical environments.

Relating target fish DNA concentration to community composition analysis in freshwater fish via metabarcoding

Metabarcoding has been widely accepted as a useful tool for biodiversity assessment based on eDNA. The method allows for the detection of entire groups of organisms in a single sample, making it particularly applicable in aquatic habitats. The high sensitivity of the molecular approaches is especially beneficial in detecting elusive and rare fish species, improving biodiversity assessments. Numerous biotic and abiotic factors that affect the persistence and availability of fish DNA in surface waters and therefore affecting species detectability, have been identified. However, little is known about the relationship between the total fish DNA concentration and the detectability of differential abundant species. In this study three controlled mock-community DNA samples (56 individual samples) were analyzed by (i) metabarcoding (MiSeq) of 12S rDNA (175 bp) and by (ii) total freshwater fish DNA quantification (via qPCR of 12S rDNA). We show that the fish DNA quantity affects the relative abundance of species-specific sequences and the detectability of rare species. In particular we found that samples with a concentration between 1000 pg/μL down to 10 pg/μL of total fish DNA revealed a stable relative frequency of DNA sequences obtained for a specific fish species, as well as a low variability between replicates. Additionally, we observed that even in complex mock-community DNA samples, a total fish DNA concentration of 23 pg/μL was sufficient to reliably detect all species in every replicate, including three rare species with proportions of ≤0.5 %. We also found that the DNA barcode similarity between species can affect detectability, if evenness is low. Our data suggest that the total DNA concentration of fish is an important factor to consider when analyzing and interpreting relative sequence abundance data. Therefore, the workflow proposed here will contribute to an ecologically and economically efficient application of metabarcoding in fish biodiversity assessment.

Water quality assessment methods of the highland Andean rivers: A scoping systematic review

Water is a resource that influences sustainable development in different ways in social, economic, and environmental aspects, being the Andes the major provider of this resource. However, they have been affected mainly by anthropogenic activities due to the proximity of settlements in the watersheds, so they tend to have more significant contamination, and their evaluation is essential to mitigate problems for those who consume them. However, despite being a fundamental resource and one of the main contributors of water, it is not so studied, so the present study aims to determine the studies based on the water quality of the high mountain rivers of the Andes by using a PRISMA methodology with the scoping review extension, based on search techniques, inclusion and exclusion criteria, and monitoring tables, in order to maintain a line of research attached to the objective of the study. After using the methodology, ten articles were obtained, which were analyzed after a bibliometric analysis to determine features of interest, such as countries in which the studies were carried out, years of publication, methodologies used, and authors' consensus. High Andean rivers' importance, the need for more studies within these areas, and the lack of suitable indexes for these unique ecosystems are highlighted.

From DNA to diagnostics: A case study using macroinvertebrate metabarcoding to assess the effectiveness of restoration measures in a Dutch stream

Stream ecosystems are under pressure due to multiple stressors. Restoration measures can halt further degradation and improve their ecological status. However, assessment of the effectiveness of the implemented measures is often insufficient because of logistic and financial constraints. DNA-metabarcoding has been proposed to scale up sample processing, although its application as a diagnostic tool has received less attention. The aim of our study was to evaluate if DNA-metabarcoding of stream macroinvertebrates can be used to compute a stressor-specific index to assess the effectiveness of a stream restoration project. For this purpose, we sampled the upstream, restored, and downstream section of a recently restored lowland stream in the Netherlands. At each site, we applied three different methods of macroinvertebrate identification: morphological identification of bulk samples (morphology), DNA-metabarcoding of the same bulk samples (DNA) and metabarcoding of eDNA extracted from the water (eDNA). First, we compared the community composition identified by each method. The communities identified by morphology and DNA were highly similar, whereas the communities generated by the eDNA differed. Second, we analysed whether the identification methods could be used to assess the effectiveness of the restoration project, focussing on a stressor-specific index for flow as the restoration measures aimed at improving flow conditions. Both the morphology and bulk DNA samples indicated improved flow conditions in the restored section of the stream (i.e., less stress from the reduction or absence of flow than in the unrestored sections). Contrary, the eDNA-water samples did not differentiate the amount of stress throughout the catchment, although applying recent developments in eDNA sampling could lead to more robust results. In conclusion, this study forms proof of concept that DNA from bulk samples can be utilized to assess the effectiveness of restoration measures, showing the added value of this approach for water managers.

Development and validation of a DNA-based multi-species biomonitoring toolkit using a high-throughput qPCR platform: A case study of Irish shellfish species

Biomonitoring of marine life has been enhanced in recent years by the integration of innovative DNA-based approaches, which offer advantages over more laborious techniques (e.g. microscopy). However, trade-offs between throughput, sensitivity and quantitative measurements must be made when choosing between the prevailing molecular methodologies (i.e. metabarcoding or qPCR/dPCR). Thus, the aim of the present study was to demonstrate the utility of a microfluidic-enabled high-throughput quantitative PCR platform (HTqPCR) for the rapid and cost-effective development and validation of a DNA-based multi-species biomonitoring toolkit, using larvae of 23 commercially targeted bivalve and crustacean species as a case study. The workflow was divided into three main phases: definition of (off-) target taxa and establishment of reference databases (PHASE 1); selection/development and assessment of molecular assays (PHASE 2); and protocol optimization and field validation (PHASE 3). 42 assays were eventually chosen and validated. Genetic signal not only showed good correlation with direct visual counts by microscopy but also showed the ability to provide quantitative data at the highest taxonomic resolution (species level) in a time- and cost-effective fashion. This study developed a biomonitoring toolkit, demonstrating the considerable advantages of this state-of-the-art technology in boosting the developmental testing and application of panels of molecular assays for the monitoring and management of natural resources. Once developed, this approach provides a cost and time-effective alternative compared to other multi-species approaches (e.g. metabarcoding). In addition, it is transferable to a wide range of species and will aid future monitoring programmes.

Improving the efficiency of adaptive management methods in multiple fishways using environmental DNA

Dams and weirs impede the continuity of rivers and transit of migratory fish. To overcome this obstacle, fishways are installed worldwide; however, management after installation is important. The Miyanaka Intake Dam has three fish ladders with different flow velocities and discharges and has been under adaptive management since 2012. Fish catch surveys, conducted as an adaptive management strategy, place a heavy burden on fish. Furthermore, a large number of investigators must be mobilized during the 30-day investigation period. Thus, a monitoring method using environmental DNA that exerts no burden on fish and requires only a few surveyors (to obtain water samples) and an in-house analyst was devised; however, its implementation in a fishway away from the point of analysis and with limited flow space and its effective water sampling frequency have not been reported. Therefore, in 2019, we started a trial aiming to evaluate the methods and application conditions of environmental DNA surveys for the continuous and long-term monitoring of various fish fauna upstream and downstream of the Miyanaka Intake Dam. To evaluate the fish fauna, the results of an environmental DNA survey (metabarcoding method) for 2019 to 2022 were compared to those of a catch survey in the fishway from 2012 to 2022. The results confirmed the use of environmental DNA surveys in evaluating the contribution of fishways to biodiversity under certain conditions and introduced a novel method for sample collection.

TICI: a taxon-independent community index for eDNA-based ecological health assessment

Global biodiversity is declining at an ever-increasing rate. Yet effective policies to mitigate or reverse these declines require ecosystem condition data that are rarely available. Morphology-based bioassessment methods are difficult to scale, limited in scope, suffer prohibitive costs, require skilled taxonomists, and can be applied inconsistently between practitioners. Environmental DNA (eDNA) metabarcoding offers a powerful, reproducible and scalable solution that can survey across the tree-of-life with relatively low cost and minimal expertise for sample collection. However, there remains a need to condense the complex, multidimensional community information into simple, interpretable metrics of ecological health for environmental management purposes. We developed a riverine taxon-independent community index (TICI) that objectively assigns indicator values to amplicon sequence variants (ASVs), and significantly improves the statistical power and utility of eDNA-based bioassessments. The TICI model training step uses the Chessman iterative learning algorithm to assign health indicator scores to a large number of ASVs that are commonly encountered across a wide geographic range. New sites can then be evaluated for ecological health by averaging the indicator value of the ASVs present at the site. We trained a TICI model on an eDNA dataset from 53 well-studied riverine monitoring sites across New Zealand, each sampled with a high level of biological replication (n = 16). Eight short-amplicon metabarcoding assays were used to generate data from a broad taxonomic range, including bacteria, microeukaryotes, fungi, plants, and animals. Site-specific TICI scores were strongly correlated with historical stream condition scores from macroinvertebrate assessments (macroinvertebrate community index or MCI; R2 = 0.82), and TICI variation between sample replicates was minimal (CV = 0.013). Taken together, this demonstrates the potential for taxon-independent eDNA analysis to provide a reliable, robust and low-cost assessment of ecological health that is accessible to environmental managers, decision makers, and the wider community.

Insights into freshwater ciliate diversity through high throughput DNA metabarcoding

The freshwater bodies of India are highly biodiverse but still understudied, especially concerning ciliates. Ciliates constitute a significant portion of eukaryotic diversity and play crucial roles in microbial loops, nutrient recycling, and ecosystem maintenance. The present study aimed to elucidate ciliate diversity in three freshwater sites in the Delhi region of India: Okhla Bird Sanctuary (OBS), Sanjay Lake (SL), and Raj Ghat pond (RJ). This study represents the first investigation into the taxonomic diversity and richness of freshwater ciliates in India using a high-throughput DNA metabarcoding approach. For the analysis, total environmental DNA was extracted from the three freshwater samples, followed by sequencing of the 18S V4 barcode region and subsequent phylogenetic analyses. Operational taxonomic units (OTU) analyses revealed maximum species diversity in OBS (106), followed by SL (104) and RJ (99) sites. Ciliates from the classes Oligohymenophorea, Prostomatea, and Spirotrichea were dominant in the three sites. The study discusses the ability of the metabarcoding approach to uncover unknown and rare species. The study highlights the need for refined reference databases and cautious interpretation of the high-throughput sequencing-generated data while emphasizing the complementary nature of molecular and morphological approaches in studying ciliate diversity.

From biomarkers to community composition: Negative effects of UV/chlorine-treated reclaimed urban wastewater on freshwater biota

The use of urban wastewater reclaimed water has recently increased across the globe to restore stream environmental flows and mitigate the effects of water scarcity. Reclaimed water is disinfected using different treatments, but their effects into the receiving rivers are little studied. Physiological bioassays and biomarkers can detect sub-lethal effects on target species, but do not provide information on changes in community structure. In contrast, official monitoring programs use community structure information but often at coarse taxonomic resolution level that may fail to detect species level impacts. Here, we combined commonly used biomonitoring approaches from organism physiology to community species composition to scan a broad range of effects of disinfection of reclaimed water by UV-light only and both UV/chlorine on the biota. We (1) performed bioassays in one laboratory species (water flea Daphnia magna) and measured biomarkers in two wild species (caddisfly Hydropsyche exocellata and the barbel Luciobarbus graellsii), (2) calculated standard indices of biotic quality (IBQ) for diatoms, benthic macroinvertebrates, and fishes, and (3) analysed community species composition of eukaryotes determined by Cytochrome Oxidase C subunit I (cox1) metabarcoding. Only the UV/chlorine treatment caused significant changes in feeding rates of D. magna and reduced antioxidant defenses, increased anaerobic metabolism and altered the levels of lipid peroxidiation in H. exocellata. However, inputs of reclaimed water were significantly associated with a greater proportion of circulating neutrophils and LG-PAS cells in L. graellsii. Despite IBQ did not discriminate between the two water treatments, metabarcoding data detected community composition changes upon exposure to UV/chlorine reclaimed water. Overall, despite the effects of UV/chlorine-treated water were transient, our study suggests that UV-light treated is less harmful for freshwater biota than UV/chlorine-treated reclaimed water, but those effects depend of the organizational level.

Environmental DNA metabarcoding reveals the effect of environmental selection on phytoplankton community structure along a subtropical river

The Dongjiang River, a major tributary of the Pearl River system that supplies water to more than 40 million people in Guangdong Province and neighboring regions of China, harbors rich biodiversity, including many endemic and endangered species. However, human activities such as urbanization, agriculture, and industrialization have posed serious threats to its water quality and biodiversity. To assess the status and drivers of phytoplankton diversity, which is a key indicator of aquatic ecosystem health, this study used Environmental DNA (eDNA) metabarcoding combined with machine learning methods to explore spatial variations in the composition and structure of phytoplankton communities along the Dongjiang River, including its estuary. The results showed that phytoplankton diversity exhibited spatial distribution patterns, with higher community structure similarity and lower network complexity in the upstream than in the downstream regions. Environmental selection was the main mechanism shaping phytoplankton community composition, with natural factors driving the dominance of Pyrrophyta, Ochrophyta, and Cryptophyta in the upstream regions and estuaries. In contrast, the downstream regions was influenced by high concentrations of pollutants, resulting in increased abundance of Cryptophyta. The random forest model identified temperature, dissolved oxygen, chlorophyll a, NO2-, and NH4+ as the main factors influencing the primary phytoplankton communities and could be used to predict changes during wet periods. This study provides valuable insights into the factors influencing phytoplankton diversity and community composition in the Dongjiang River, and demonstrates the application value of eDNA metabarcoding technique in large-scale, long-distance river biodiversity monitoring.

A broad-taxa approach as an important concept in ecotoxicological studies and pollution monitoring

Aquatic invertebrates play a pivotal role in (eco)toxicological assessments because they offer ethical, cost-effective and repeatable testing options. Additionally, their significance in the food chain and their ability to represent diverse aquatic ecosystems make them valuable subjects for (eco)toxicological studies. To ensure consistency and comparability across studies, international (eco)toxicology guidelines have been used to establish standardised methods and protocols for data collection, analysis and interpretation. However, the current standardised protocols primarily focus on a limited number of aquatic invertebrate species, mainly from Arthropoda, Mollusca and Annelida. These protocols are suitable for basic toxicity screening, effectively assessing the immediate and severe effects of toxic substances on organisms. For more comprehensive and ecologically relevant assessments, particularly those addressing long-term effects and ecosystem-wide impacts, we recommended the use of a broader diversity of species, since the present choice of taxa exacerbates the limited scope of basic ecotoxicological studies. This review provides a comprehensive overview of (eco)toxicological studies, focusing on major aquatic invertebrate taxa and how they are used to assess the impact of chemicals in diverse aquatic environments. The present work supports the use of a broad-taxa approach in basic environmental assessments, as it better represents the natural populations inhabiting various ecosystems. Advances in omics and other biochemical and computational techniques make the broad-taxa approach more feasible, enabling mechanistic studies on non-model organisms. By combining these approaches with in vitro techniques together with the broad-taxa approach, researchers can gain insights into less-explored impacts of pollution, such as changes in population diversity, the development of tolerance and transgenerational inheritance of pollution responses, the impact on organism phenotypic plasticity, biological invasion outcomes, social behaviour changes, metabolome changes, regeneration phenomena, disease susceptibility and tissue pathologies. This review also emphasises the need for harmonised data-reporting standards and minimum annotation checklists to ensure that research results are findable, accessible, interoperable and reusable (FAIR), maximising the use and reusability of data. The ultimate goal is to encourage integrated and holistic problem-focused collaboration between diverse scientific disciplines, international standardisation organisations and decision-making bodies, with a focus on transdisciplinary knowledge co-production for the One-Health approach.

Comparison of different macroinvertebrates bioassessment indices in a large near-natural watershed under the context of metacommunity theory

The metacommunity theory proposes that community structure and biodiversity are influenced by both local processes (such as environmental filtering) and regional processes (such as dispersal). Despite the extensive use of traditional bioassessments based on species-environment relationships, the impact of dispersal processes on these assessments has been largely overlooked. This study aims to compare correlations between various bioassessment indices, including Shannon Weiner (H'), Biological Monitoring Working Party (BMWP), average score per taxon (ASPT), biotic index (BI), and EPT taxa index (EPT), based on macroinvertebrates collected from 147 sampling sites in a subtropical Chinese near-natural catchment. Modified indices were calculated by removing species strongly influenced by dispersal processes to address the influence of dispersal processes. Their relationship with environmental factors was then compared to the original indices. The study employed random forest regression (RFR) to compare the explanatory power of environmental factors using the two sets of indices. The spearman rank correlation analysis was conducted to examine the correlation between indices and environmental factors. The river health assessment was performed based on both modified and original indices. The results reveal significant differences between original and modified indices (especially H' and BI) providing a more accurate reflection of environmental conditions. Furthermore, the sensitivity of the different indices to various environmental factors varied, leading to differences in the bioassessment results between the modified and the original indices. Notably, original H', BMWP, and ASPT overestimated the bioassessment results, whereas the original BI underestimated them. These findings offer valuable insights into bioassessment and river health assessment evaluation within the catchment and other interconnected freshwater ecosystems, such as lakes, reservoirs, and wetlands. Our study underscores the importance of assessing and mitigating the impact of dispersal processes on bioassessment to obtain a more precise representation of the status of freshwater ecosystems.

The Evolution, Assembly, and Dynamics of Marine Holobionts

The holobiont concept (i.e., multiple living beings in close symbiosis with one another and functioning as a unit) is revolutionizing our understanding of biology, especially in marine systems. The earliest marine holobiont was likely a syntrophic partnership of at least two prokaryotic members. Since then, symbiosis has enabled marine organisms to conquer all ocean habitats through the formation of holobionts with a wide spectrum of complexities. However, most scientific inquiries have focused on isolated organisms and their adaptations to specific environments. In this review, we attempt to illustrate why a holobiont perspective-specifically, the study of how numerous organisms form a discrete ecological unit through symbiosis-will be a more impactful strategy to advance our understanding of the ecology and evolution of marine life. We argue that this approach is instrumental in addressing the threats to marine biodiversity posed by the current global environmental crisis.

DNA Sequencing Technologies and DNA Barcoding

DNA barcodes are short, standardized DNA segments that geneticists can use to identify all living taxa. On the other hand, DNA barcoding identifies species by analyzing these specific regions against a DNA barcode reference library. In its initial years, DNA barcodes sequenced by Sanger's method were extensively used by taxonomists for the characterization and identification of species. But in recent years, DNA barcoding by next-generation sequencing (NGS) has found broader applications, such as quality control, biomonitoring of protected species, and biodiversity assessment. Technological advancements have also paved the way to metabarcoding, which has enabled massive parallel sequ.encing of complex bulk samples using high-throughput sequencing techniques. In future, DNA barcoding along with high-throughput techniques will show stupendous progress in taxonomic classification with reference to available sequence data.

Sedimentary DNA reveals the link between microbial community dynamics and climate during the late last glaciation in the offshore region of the Zambezi River, Southwest Indian Ocean

Reconstructing the relationship between microbial communities and past abrupt climate change is of great importance for understanding current biodiversity patterns and predicting changes under future climate scenarios. However, little is currently known about how microbial communities respond to changes in key environmental stages due to a lack of research in this area. Here, we examine the variability in the communities of bacteria, archaea, and fungi from sediments deposited offshore region of the Zambezi River between 21.7 and 9.6 thousand years ago (ka) (covering the last glacial maximum, or LGM, and the early Holocene) using DNA metabarcoding approach via high-throughput sequencing. The results showed that (1) microbial assemblages differed across three key time intervals, with the last deglaciation having the most homogeneous prokaryotic assemblages, while for fungal communities in the LGM, and the early Holocene and LGM differing the most; (2) the warm early Holocene showed the highest diversity, whereas the lowest diversity was found in the LGM; and (3) the selected indicator species better reflected the climatic characteristics of different environmental stages. These results highlight the power of ancient sedimentary DNA to refine our understanding of microbial dynamics in marine sedimentary systems near large rivers, thus providing a basis for better modeling ecological processes in further research.

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

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

Simultaneous analysis of seven 16S rRNA hypervariable gene regions increases efficiency in marine bacterial diversity detection

Environmental DNA sequencing is the gold standard to reveal microbial community structures. In most applications, a one-fragment PCR approach is applied to amplify a taxonomic marker gene, usually a hypervariable region of the 16S rRNA gene. We used a new reverse complement (RC)-PCR-based assay that amplifies seven out of the nine hypervariable regions of the 16S rRNA gene, to interrogate bacterial communities in sediment samples collected from different coastal marine sites with an impact gradient. In parallel, we employed a traditional one-fragment analysis of the hypervariable V3-V4 region to investigate whether the RC-PCR reveals more of the 'unseen' diversity obtained by the one-fragment approach. As a benchmark for the full deck of diversity, we subjected the samples to PCR-free metagenomic sequencing. None of the two PCR-based approaches recorded the full taxonomic repertoire obtained from the metagenomics datasets. However, the RC-PCR approach detected 2.8 times more bacterial genera compared to the near-saturation sequenced V3-V4 samples. RC-PCR is an ideal compromise between the standard one-fragment approach and metagenomics sequencing and may guide future environmental sequencing studies, in which bacterial diversity is a central subject.

Insights in Pharmaceutical Pollution: The Prospective Role of eDNA Metabarcoding

Environmental pollution is a growing threat to natural ecosystems and one of the world's most pressing concerns. The increasing worldwide use of pharmaceuticals has elevated their status as significant emerging contaminants. Pharmaceuticals enter aquatic environments through multiple pathways related to anthropogenic activity. Their high consumption, insufficient waste treatment, and the incapacity of organisms to completely metabolize them contribute to their accumulation in aquatic environments, posing a threat to all life forms. Various analytical methods have been used to quantify pharmaceuticals. Biotechnology advancements based on next-generation sequencing (NGS) techniques, like eDNA metabarcoding, have enabled the development of new methods for assessing and monitoring the ecotoxicological effects of pharmaceuticals. eDNA metabarcoding is a valuable biomonitoring tool for pharmaceutical pollution because it (a) provides an efficient method to assess and predict pollution status, (b) identifies pollution sources, (c) tracks changes in pharmaceutical pollution levels over time, (d) assesses the ecological impact of pharmaceutical pollution, (e) helps prioritize cleanup and mitigation efforts, and (f) offers insights into the diversity and composition of microbial and other bioindicator communities. This review highlights the issue of aquatic pharmaceutical pollution while emphasizing the importance of using modern NGS-based biomonitoring actions to assess its environmental effects more consistently and effectively.

Supervised machine learning improves general applicability of eDNA metabarcoding for reservoir health monitoring

Effective and standardized monitoring methodologies are vital for successful reservoir restoration and management. Environmental DNA (eDNA) metabarcoding sequencing offers a promising alternative for biomonitoring and can overcome many limitations of traditional morphological bioassessment. Recent attempts have even shown that supervised machine learning (SML) can directly infer biotic indices (BI) from eDNA metabarcoding data, bypassing the cumbersome calculation process of BI regardless of the taxonomic assignment of eDNA sequences. However, questions surrounding the general applicability of this taxonomy-free approach to monitoring reservoir health remain unclear, including model stability, feature selection, algorithm choice, and multi-season biomonitoring. Here, we firstly developed a novel biological integrity index (Me-IBI) that integrates multitrophic interactions and environmental information, based on taxonomy-assigned eDNA metabarcoding data. The Me-IBI can better distinguish the actual health status of the Three Gorges Reservoir (TGR) than physicochemical assessments and have a clear response to human activity. Then, taking this reliable Me-IBI as a supervised label, we compared the impact of selecting different numbers of features and SML algorithms on the stability and predictive performance of the model for predicting ecological conditions in multiple seasons using taxonomy-free eDNA metabarcoding data. We discovered that even with a small number of features, different SML algorithms can establish a stable model and obtain excellent predictive performance. Finally, we proposed a four-step strategy for standardized routine biomonitoring using SML tools. Our study firstly explores the general applicability problem of the taxonomy-free eDNA-SML approach and establishes a solid foundation for the large-scale and standardized biomonitoring application.

Environmental DNA detects biodiversity and ecological features of phytoplankton communities in Mediterranean transitional waters

Climate changes and anthropogenic pressures are causing a biodiversity decline in terms of species number and genetic diversity, reducing the adaptability and evolvability of natural communities. Transitional water ecosystems are more sensitive to habitat reduction and degradation and, thus, are more exposed to biodiversity declines requiring biodiversity monitoring programs for their conservation. Environmental DNA (eDNA) metabarcoding represents a high-throughput tool for biodiversity assessment that is facilitating data collection for biodiversity monitoring. In this study, we applied, for the first time, eDNA metabarcoding in a Mediterranean coastal lagoon to assess the ecological features of eukaryotic phytoplankton communities. We sampled water in seven different lagoon sites and amplified the extracted DNA with primers targeting the variable region 4 (V4) of the 18S rRNA gene marker. The results demonstrated the validity of eDNA studies to provide insights into lagoon phytoplankton composition, establish the structure and spatial variation of phytoplankton communities, and evaluate its correlation to abiotic factors. Finally, the genetic distances analysis suggests that the different spatial distribution of OTUs, at least for the Tetraselmis genus, reflects the genetic background.

Environmental DNA metabarcoding serves as a promising method for aquatic species monitoring and management: A review focused on its workflow, applications, challenges and prospects

Marine and freshwater biodiversity is under threat from both natural and manmade causes. Biological monitoring is currently a top priority for biodiversity protection. Given present limitations, traditional biological monitoring methods may not achieve the proposed monitoring aims. Environmental DNA metabarcoding technology reflects species information by capturing and extracting DNA from environmental samples, using molecular biology techniques to sequence and analyze the DNA, and comparing the obtained information with existing reference libraries to obtain species identification. However, its practical application has highlighted several limitations. This paper summarizes the main steps in the environmental application of eDNA metabarcoding technology in aquatic ecosystems, including the discovery of unknown species, the detection of invasive species, and evaluations of biodiversity. At present, with the rapid development of big data and artificial intelligence, certain advanced technologies and devices can be combined with environmental DNA metabarcoding technology to promote further development of aquatic species monitoring and management.

A combination of machine-learning and eDNA reveals the genetic signature of environmental change at the landscape levels

The current advances of environmental DNA (eDNA) bring profound changes to ecological monitoring and provide unique insights on the biological diversity of ecosystems. The very nature of eDNA data is challenging yet also revolutionizing how biological monitoring information is analysed. In particular, new metrics and approaches should take full advantage of the extent and detail of molecular data produced by genetic methods. In this perspective, machine learning algorithms are particularly promising as they can capture complex relationships between the multiple environmental pressures and the diversity of biological communities. We investigated the potential of a new generation of biomonitoring tools that implement machine-learning techniques to fully exploit eDNA datasets. We trained a machine learning model to discriminate between reference and impacted communities of freshwater macroinvertebrates and assessed its performances using a large eDNA dataset collected at 64 standard federal monitoring sites across Switzerland. We show that a model trained on eDNA is significantly better than a naive model and performs similarly to a model trained on traditional data. Our proof-of-concept shows that such a combination of eDNA and machine learning approaches has the potential to complement or even replace traditional environmental monitoring, and could be scaled along temporal or spatial dimensions.

Bacterial bioindicators enable biological status classification along the continental Danube river

Despite the importance of bacteria in aquatic ecosystems and their predictable diversity patterns across space and time, biomonitoring tools for status assessment relying on these organisms are widely lacking. This is partly due to insufficient data and models to identify reliable microbial predictors. Here, we show metabarcoding in combination with multivariate statistics and machine learning allows to identify bacterial bioindicators for existing biological status classification systems. Bacterial beta-diversity dynamics follow environmental gradients and the observed associations highlight potential bioindicators for ecological outcomes. Spatio-temporal links spanning the microbial communities along the river allow accurate prediction of downstream biological status from upstream information. Network analysis on amplicon sequence veariants identify as good indicators genera Fluviicola, Acinetobacter, Flavobacterium, and Rhodoluna, and reveal informational redundancy among taxa, which coincides with taxonomic relatedness. The redundancy among bacterial bioindicators reveals mutually exclusive taxa, which allow accurate biological status modeling using as few as 2-3 amplicon sequence variants. As such our models show that using a few bacterial amplicon sequence variants from globally distributed genera allows for biological status assessment along river systems.

It Is Not Only Data-Freshwater Invertebrates Misused in Biological Monitoring

The article presents and discusses the issues of the use of free-living invertebrates to assess the ecological status of freshwater environments with different methods of biological monitoring. Invertebrates are excluded from ethical consideration in the procedures of environmental protection, which results in the killing of many more individuals during sampling than necessary. Biomonitoring is used as a routine method for environmental protection that results in the cruel death of even millions of aquatic animals annually. In many cases, the mortality of animals used in such types of activities has been shown as excessive, e.g., because the vast majority die due to unnecessary subsampling procedures. Improperly planned and conducted procedures which result in excessive mortality have or may have a negative impact on the environment and biodiversity. Their existence as sensitive beings is reduced to an information function; they become only data useful for biomonitoring purposes. The main problem when trying to determine the mortality of invertebrates due to biomonitoring activities and its impact on natural populations seems to be the lack of access to raw data presenting how many animals were killed during sampling.

Using eDNA to understand predator-prey interactions influenced by invasive species

Invasive predatory species may alter population dynamic processes of their prey and impact biological communities and ecosystem processes. Revealing biotic interactions, however, including the relationship between predator and prey, is a difficult task, in particular for species that are hard to monitor. Here, we present a case study that documents the utility of environmental DNA analysis (eDNA) to assess predator-prey interactions between two invasive fishes (Lepomis gibbosus, Pseudorasbora parva) and two potential amphibian prey species, (Triturus cristatus, Pelobates fuscus). We used species-specific TaqMan assays for quantitative assessment of eDNA concentrations from water samples collected from 89 sites across 31 ponds during three consecutive months from a local amphibian hotspot in Germany. We found a negative relationship between eDNA concentrations of the predators (fishes) and prey (amphibians) using Monte-Carlo tests. Our study highlights the potential of eDNA application to reveal predator-prey interactions and confirms the hypothesis that the observed local declines of amphibian species may be at least partly caused by recently introduced invasive fishes. Our findings have important consequences for local conservation management and highlight the usefulness of eDNA approaches to assess ecological interactions and guide targeted conservation action.

Haplotype-level metabarcoding of freshwater macroinvertebrate species: A prospective tool for population genetic analysis

Metabarcoding is a molecular-based tool capable of large quantity high-throughput species identification from bulk samples that is a faster and more cost-effective alternative to conventional DNA-sequencing approaches. Still, further exploration and assessment of the laboratory and bioinformatics strategies are required to unlock the potential of metabarcoding-based inference of haplotype information. In this study, we assessed the inference of freshwater macroinvertebrate haplotypes from metabarcoding data in a mock sample. We also examined the influence of DNA template concentration and PCR cycle on detecting true and spurious haplotypes. We tested this strategy on a mock sample containing twenty individuals from four species with known haplotypes based on the 658-bp Folmer region of the mitochondrial cytochrome c oxidase gene. We recovered fourteen zero-radius operational taxonomic units (zOTUs) of 421-bp length, with twelve zOTUs having a 100% match with the Sanger haplotype sequences. High-quality reads relatively increased with increasing PCR cycles, and the relative abundance of each zOTU was consistent for each cycle. This suggests that increasing the PCR cycles from 24 to 64 did not affect the relative abundance of each zOTU. As metabarcoding becomes more established and laboratory protocols and bioinformatic pipelines are continuously being developed, our study demonstrated the method's ability to infer intraspecific variability while highlighting the challenges that must be addressed before its eventual application for population genetic studies.

Genetic Markers for Metabarcoding of Freshwater Microalgae: Review

The metabarcoding methods for studying the diversity of freshwater microalgae and routine biomonitoring are actively used in modern research. A lot of experience has been accumulated already, and many methodological questions have been solved (such as the influence of the methods and time of sample conservation, DNA extraction and bioinformatical processing). The reproducibility of the method has been tested and confirmed. However, one of the main problems-choosing a genetic marker for the study-still lacks a clear answer. We analyzed 70 publications and found out that studies on eukaryotic freshwater microalgae use 12 markers (different nuclear regions 18S and ITS and plastids rbcL, 23S and 16S). Each marker has its peculiarities; they amplify differently and have various levels of efficiency (variability) in different groups of algae. The V4 and V9 18S and rbcL regions are used most often. We concentrated especially on the studies that compare the results of using different markers and microscopy. We summarize the data on the primers for each region and on how the choice of a marker affects the taxonomic composition of a community.

Pest Alert Tool-a web-based application for flagging species of concern in metabarcoding datasets

Advances in high-throughput sequencing (HTS) technologies and their increasing affordability have fueled environmental DNA (eDNA) metabarcoding data generation from freshwater, marine and terrestrial ecosystems. Research institutions worldwide progressively employ HTS for biodiversity assessments, new species discovery and ecological trend monitoring. Moreover, even non-scientists can now collect an eDNA sample, send it to a specialized laboratory for analysis and receive in-depth biodiversity record from a sampling site. This offers unprecedented opportunities for biodiversity assessments across wide temporal and spatial scales. The large volume of data produced by metabarcoding also enables incidental detection of species of concern, including non-indigenous and pathogenic organisms. We introduce an online app-Pest Alert Tool-for screening nuclear small subunit 18S ribosomal RNA and mitochondrial cytochrome oxidase subunit I datasets for marine non-indigenous species as well as unwanted and notifiable marine organisms in New Zealand. The output can be filtered by minimum length of the query sequence and identity match. For putative matches, a phylogenetic tree can be generated through the National Center for Biotechnology Information's BLAST Tree View tool, allowing for additional verification of the species of concern detection. The Pest Alert Tool is publicly available at https://pest-alert-tool-prod.azurewebsites.net/.

Zeta diversity patterns in metabarcoded lotic algal assemblages as a tool for bioassessment

Assessments of the ecological health of algal assemblages in streams typically focus on measures of their local diversity and classify individuals by morphotaxonomy. Such assemblages are often connected through various ecological processes, such as dispersal, and may be more accurately assessed as components of regional-, rather than local-scale assemblages. With recent declines in the costs of sequencing and computation, it has also become increasingly feasible to use metabarcoding to more accurately classify algal species and perform regional-scale bioassessments. Recently, zeta diversity has been explored as a novel method of constructing regional bioassessments for groups of streams. Here, we model the use of zeta diversity to investigate whether stream health can be determined by the landscape diversity of algal assemblages. We also compare the use of DNA metabarcoding and morphotaxonomy classifications in these zeta diversity-based bioassessments of regional stream health. From 96 stream samples in California, we used various orders of zeta diversity to construct models of biotic integrity for multiple assemblages of diatoms, as well as hybrid assemblages of diatoms in combination with soft-bodied algae, using taxonomy data generated with both DNA sequencing as well as traditional morphotaxonomic approaches. We compared our ability to evaluate the ecological health of streams with the performance of multiple algal indices of biological condition. Our zeta diversity-based models of regional biotic integrity were more strongly correlated with existing indices for algal assemblages classified using metabarcoding compared to morphotaxonomy. Metabarcoding for diatoms and hybrid algal assemblages involved rbcL and 18S V9 primers, respectively. Importantly, we also found that these algal assemblages, independent of the classification method, are more likely to be assembled under a process of niche differentiation rather than stochastically. Taken together, these results suggest the potential for zeta diversity patterns of algal assemblages classified using metabarcoding to inform stream bioassessments.

Progress in understanding the vulnerability of freshwater ecosystems

The ability to collect and synthesize long-term environmental monitoring data is essential for the effective management of freshwater ecosystems. Progress has been made in assessment and monitoring approaches that have integrated routine monitoring programs into more holistic watershed-scale vulnerability assessments. While the concept of vulnerability assessment is well-defined for ecosystems, complementary and sometimes competing concepts of adaptive management, ecological integrity, and ecological condition complicate the communication of results to a broader audience. Here, we identify progress in freshwater assessments that can contribute to the identification and communication of freshwater vulnerability. We review novel methods that address common challenges associated with: 1) a lack of baseline information, 2) variability associated with a spatial context, and 3) the taxonomic sufficiency of biological indicators used to make inferences about ecological conditions. Innovation in methods and communication are discussed as a means to highlight meaningful cost-effective results that target policy towards heuristic ecosystem-management.

Assessing arthropod diversity metrics derived from stream environmental DNA: spatiotemporal variation and paired comparisons with manual sampling

Background

Benthic invertebrate (BI) surveys have been widely used to characterize freshwater environmental quality but can be challenging to implement at desired spatial scales and frequency. Environmental DNA (eDNA) allows an alternative BI survey approach, one that can potentially be implemented more rapidly and cheaply than traditional methods.

Methods

We evaluated eDNA analogs of BI metrics in the Potomac River watershed of the eastern United States. We first compared arthropod diversity detected with primers targeting mitochondrial 16S (mt16S) and cytochrome c oxidase 1 (cox1 or COI) loci to that detected by manual surveys conducted in parallel. We then evaluated spatial and temporal variation in arthropod diversity metrics with repeated sampling in three focal parks. We also investigated technical factors such as filter type used to capture eDNA and PCR inhibition treatment.

Results

Our results indicate that genus-level assessment of eDNA compositions is achievable at both loci with modest technical noise, although database gaps remain substantial at mt16S for regional taxa. While the specific taxa identified by eDNA did not strongly overlap with paired manual surveys, some metrics derived from eDNA compositions were rank-correlated with previously derived biological indices of environmental quality. Repeated sampling revealed statistical differences between high- and low-quality sites based on taxonomic diversity, functional diversity, and tolerance scores weighted by taxon proportions in transformed counts. We conclude that eDNA compositions are efficient and informative of stream condition. Further development and validation of scoring schemes analogous to commonly used biological indices should allow increased application of the approach to management needs.

Phytoplankton and Microzooplankton Community Structure and Assembly Mechanisms in Northwestern Pacific Ocean Estuaries with Environmental Heterogeneity and Geographic Segregation

Phytoplankton and microzooplankton are crucial players in marine ecosystems and first responders to environmental changes, but their community structures and how they are shaped by environmental conditions have rarely been studied simultaneously. In this study, we conducted an eDNA metabarcoding sequencing combined with multiple statistical methods to simultaneously analyze the phytoplankton and microzooplankton in Liaohe (LH) and Yalujiang (YLJ) estuaries. The major objective was to examine how plankton community structure and assembly mechanism may differ between two estuaries with similar latitudinal position and climate but geographical segregation and differential level of urbanization (more in LH). Clear differences in diversity and composition of phytoplankton and microzooplankton communities between LH and YLJ estuaries were observed. Richness of phytoplankton was significantly higher in LH than YLJ, while richness of microzooplankton was higher in YLJ. The magnitude of intrahabitat variations in phytoplankton communities was significantly stronger than that of microzooplankton. Some phytoplankton and microzooplankton taxa also showed interhabitat differences in their relative abundances. Phytoplankton showed a stronger geographic distance-decay of similarity than microzooplankton, while significant environmental distance-decay of similarity in microzooplankton was found in the less urbanized YLJ estuary. Community assembly of phytoplankton was, based on the neutral community models, driven primarily by stochastic processes, while deterministic processes contributed more for microzooplankton. Furthermore, we detected wider habitat niche breadths and stronger dispersal abilities in phytoplankton than in microzooplankton. These results suggest that passive dispersal shapes the phytoplankton community whereas environmental selection shapes the microzooplankton community. IMPORTANCE Understanding the underlying mechanisms shaping a metacommunity is useful to management for improving the ecosystem function. The research presented in the manuscript mainly tried to address the effects of habitat geography and environmental conditions on the phytoplankton and microzooplankton communities, and the underlying mechanisms of community assembly in temperate estuaries. In order to achieve this purpose, we developed a metabarcoding sequencing method based on 18S rRNA gene. The phytoplankton and microzooplankton communities from two estuaries with similar latitude and climatic conditions but obvious geographical segregation and significant environmental heterogeneity were investigated. The results of our study could lay a solid foundation for ascertaining phytoplankton and microzooplankton communities in estuaries with obvious environmental heterogeneity and geographic segregation and mechanisms underlying community assembly.

Future of DNA-based insect monitoring

Insects are crucial for ecosystem health but climate change and pesticide use are driving massive insect decline. To mitigate this loss, we need new and effective monitoring techniques. Over the past decade there has been a shift to DNA-based techniques. We describe key emerging techniques for sample collection. We suggest that the selection of tools should be broadened, and that DNA-based insect monitoring data need to be integrated more rapidly into policymaking. We argue that there are four key areas for advancement, including the generation of more complete DNA barcode databases to interpret molecular data, standardisation of molecular methods, scaling up of monitoring efforts, and integrating molecular tools with other technologies that allow continuous, passive monitoring based on images and/or laser imaging, detection, and ranging (LIDAR).

Maximising signal-to-noise ratios in environmental DNA-based monitoring

Man's impacts on global ecosystems are increasing and there is a growing demand that these activities be appropriately monitored. Monitoring requires measurement of a response metric ('signal') that changes maximally and consistently in response to the monitored activity irrespective of other factors ('noise'), thus maximising the signal-to-noise ratio. Indices derived from time-consuming morphology-based taxonomic identification of organisms are a core part of many monitoring programmes. Metabarcoding is an alternative to morphology-based identification and involves the sequencing of short fragments of DNA ('markers') from multiple taxa simultaneously. DNA suitable for metabarcoding includes that extracted from environmental samples (eDNA). Metabarcoding outputs DNA sequences that can be identified (annotated) by matching them against archived annotated sequences. However, sequences from most organisms are not archived - preventing annotation and potentially limiting metabarcoding in monitoring applications. Consequently, there is growing interest in using unannotated sequences as response metrics in monitoring programmes. We compared the sequences from three commonly used markers (16S (V3/V4 regions), 18S (V1/V2 regions) and COI) and, sampling along steep impact gradients, showed that the 16S and COI sequences were associated with the largest and smallest signal-to-noise ratio respectively. We trialled four separate, intuitive, noise-reduction approaches and demonstrated that removing less frequent sequences improved the signal-to-noise ratio, partitioning an additional 25 % from noise to explanatory factors in non-parametric ANOVA (NPA) and reducing dispersion in the data. For the 16S marker, retaining only the most frequently observed sequence, per sample, resulting in nine sequences across 150 samples, generated a near-maximal signal-to-noise ratio (95 % of the variance explained in NPA). We recommend that NPA, combined with rigorous elimination of less frequent sequences, be used to pre-filter sequences/taxa being used in monitoring applications. Our approach will simplify downstream analysis, for example the identification of key taxa and functional associations.

The Neolithic site "La Marmotta": DNA metabarcoding to identify the microbial deterioration of waterlogged archeological wood

Introduction

The evaluation of biological degradation of waterlogged archeological wood is crucial to choose the conservative and protective treatments to be applied to the wooden material. The waterlogged environmental conditions are characterized by oxygen scarcity, only allowing the growth of adapted microbes capable to degrade the organic wooden material, mainly erosion bacteria and soft-rot fungi. In this work, we characterized and evaluated the biodegradation state and the microbial communities of wooden fragments preserved in storage tanks. These were preserved by waterlogging within the Neolithic village "La Marmotta," currently found under the Bracciano Lake (Lazio, Italy).

Methods

The waterlogged wood samples were first identified taxonomically with an optical microscope, also allowing an evaluation of their preservation state. The microbial community was then evaluated through the sequencing of Internal Transcribed Spacer sequences for fungi and 16S for bacteria with the Oxford Nanopore Technologies (ONT) MinION platform.

Results

The identified microbial community appears to be consistent with the waterlogged samples, as many bacteria attributable to the erosion of wood and ligninolytic fungi have been sequenced.

Discussion

The reported results highlight the first use of targeted metabarcoding by ONT applied to study the biodeterioration of waterlogged archeological wood.

Metabarcoding reveals high genetic diversity of harmful algae in the coastal waters of Texas, Gulf of Mexico

Environmental-DNA (eDNA) for metabarcoding is a rapid and effective means to investigate microplankton community composition and species diversity. The objective of this study was to examine the genetic diversity of the phytoplankton community in the Gulf of Mexico, with particular emphasis on harmful algal bloom species. Samples were collected at stations along the coast of Texas in September-October 2017 that were inundated by low salinity waters in the aftermath of Hurricane Harvey. Metabarcodes were generated from the eDNA targeting both the V4 and V8-V9 regions of the 18S rDNA gene. Evaluation of the metabarcodes revealed an unexpectedly high number of harmful algal species during this short period, including five that had not been documented in this region previously. A total of 36 harmful algal species could be differentiated based on V4 and V8-V9 metabarcode markers. Using a phylogenetic approach, the taxonomic resolution of each marker differed and not all species could be differentiated using solely one marker. The V4 region resolved species within some genera (e.g., Heterocapsa), while the V8-V9 marker was necessary to resolve species within other genera (e.g., Chattonella). In other cases, species differentiation within a genus required a combination of both markers (e.g., Prorocentrum, Karenia), or another marker will be needed to resolve all species (e.g., Alexandrium, Dinophysis). We conclude that no single marker can delineate all species, so it is recommended HAB monitoring programs use more than one marker. Overall, the observed diversity of HAB species along the Texas coast using metabarcoding exceeded reports from other parts of the world.

Progress in 'taxonomic sufficiency' in aquatic biological investigations

The 'taxonomic sufficiency' (TS) approach has been applied to algae, protists, invertebrates, and vertebrates, generally by aggregating species-level abundance data to a higher taxonomic level, where genus-level data are often highly correlated with species-level data and are a valid proxy level. The TS approach offers the possibility of a comparison of data from different geographical areas and highlights the effects of contaminants. The TS approach is stable in the face of different researchers and in the comparison of long-term biological survey data. The effectiveness of the TS approach may increase with increasing environmental gradients or spatial area. The TS approach should be avoided when the spatial area is small and small differences in species-level data are considered important, so as not to cancel out the distribution patterns specific to the local environment of the biological taxa.

Recommendations for the preservation of environmental samples in diatom metabarcoding studies

Implementation of DNA metabarcoding for diatoms for environmental monitoring is now moving from a research to an operational phase, requiring rigorous guidelines and standards. In particular, the first steps of the diatom metabarcoding process, which consist of sampling and storage, have been addressed in various ways in scientific and pilot studies and now need to be rationalised. The objective of this study was to compare three currently applied preservation protocols through different storage durations (ranging from one day to one year) for phytobenthos and phytoplankton samples intended for diatom DNA metabarcoding analysis. The experimental design used samples from four freshwater and two marine sites of diverse ecological characteristics. The impact of the sample preservation and storage duration was assessed through diatom metabarcoding endpoints: DNA quality and quantity, diversity and richness, diatom assemblage composition and ecological index values (for freshwater samples). The yield and quality of extracted DNA only decreased for freshwater phytobenthos samples preserved with ethanol. Diatom diversity was not affected and their taxonomic composition predominantly reflected the site origin. Only rare taxa (< 100 reads) differed among preservation methods and storage durations. For biomonitoring purposes, freshwater ecological index values were not affected by the preservation method and storage duration tested (including ethanol preservation), all treatments returning the same ecological status for a site. This study contributes to consolidating diatom metabarcoding. Thus, accompanied by operational standards, the method will be ready to be confidently deployed and prescribed in future regulatory monitoring.

Benthic invertebrates in Svalbard fjords-when metabarcoding does not outperform traditional biodiversity assessment

To protect and restore ecosystems and biodiversity is one of the 10 challenges identified by the United Nations's Decade of the Ocean Science. In this study we used eDNA from sediments collected in two fjords of the Svalbard archipelago and compared the taxonomic composition with traditional methods through metabarcoding, targeting mitochondrial CO1, to survey benthos. Clustering of 21.6 mill sequence reads with a d value of 13 in swarm, returned about 25 K OTU reads. An identification search with the BOLD database returned 12,000 taxonomy annotated sequences spanning a similarity range of 50% to 100%. Using an acceptance filter of minimum 90% similarity to the CO1 reference sequence, we found that 74% of the ca 100 taxon identified sequence reads were Polychaeta and 22% Nematoda. Relatively few other benthic invertebrate species were detected. Many of the identified sequence reads were extra-organismal DNA from terrestrial, planktonic, and photic zone sources. For the species rich Polychaeta, we found that, on average, only 20.6% of the species identified from morphology were also detected with DNA. This discrepancy was not due to missing reference sequences in the search database, because 90-100% (mean 96.7%) of the visually identified species at each station were represented with barcodes in Boldsystems. The volume of DNA samples is small compared with the volume searched in visual sorting, and the replicate DNA-samples in sum covered only about 2% of the surface area of a grab. This may considerably reduce the detection rate of species that are not uniformly distributed in the sediments. Along with PCR amplification bias and primer mismatch, this may be an important reason for the limited congruence of species identified with the two approaches. However, metabarcoding also identified 69 additional species that are usually overlooked in visual sample sorting, demonstrating how metabarcoding can complement traditional methodology by detecting additional, less conspicuous groups of organisms.

A manager's guide to using eDNA metabarcoding in marine ecosystems

Environmental DNA (eDNA) metabarcoding is a powerful tool that can enhance marine ecosystem/biodiversity monitoring programs. Here we outline five important steps managers and researchers should consider when developing eDNA monitoring program: (1) select genes and primers to target taxa; (2) assemble or develop comprehensive barcode reference databases; (3) apply rigorous site occupancy based decontamination pipelines; (4) conduct pilot studies to define spatial and temporal variance of eDNA; and (5) archive samples, extracts, and raw sequence data. We demonstrate the importance of each of these considerations using a case study of eDNA metabarcoding in the Ports of Los Angeles and Long Beach. eDNA metabarcoding approaches detected 94.1% (16/17) of species observed in paired trawl surveys while identifying an additional 55 native fishes, providing more comprehensive biodiversity inventories. Rigorous benchmarking of eDNA metabarcoding results improved ecological interpretation and confidence in species detections while providing archived genetic resources for future analyses. Well designed and validated eDNA metabarcoding approaches are ideally suited for biomonitoring applications that rely on the detection of species, including mapping invasive species fronts and endangered species habitats as well as tracking range shifts in response to climate change. Incorporating these considerations will enhance the utility and efficacy of eDNA metabarcoding for routine biomonitoring applications.

Discovering marine biodiversity in the 21st century

We review the current knowledge of the biodiversity of the ocean as well as the levels of decline and threat for species and habitats. The lack of understanding of the distribution of life in the ocean is identified as a significant barrier to restoring its biodiversity and health. We explore why the science of taxonomy has failed to deliver knowledge of what species are present in the ocean, how they are distributed and how they are responding to global and regional to local anthropogenic pressures. This failure prevents nations from meeting their international commitments to conserve marine biodiversity with the results that investment in taxonomy has declined in many countries. We explore a range of new technologies and approaches for discovery of marine species and their detection and monitoring. These include: imaging methods, molecular approaches, active and passive acoustics, the use of interconnected databases and citizen science. Whilst no one method is suitable for discovering or detecting all groups of organisms many are complementary and have been combined to give a more complete picture of biodiversity in marine ecosystems. We conclude that integrated approaches represent the best way forwards for accelerating species discovery, description and biodiversity assessment. Examples of integrated taxonomic approaches are identified from terrestrial ecosystems. Such integrated taxonomic approaches require the adoption of cybertaxonomy approaches and will be boosted by new autonomous sampling platforms and development of machine-speed exchange of digital information between databases.

Environmental DNA-based profiling of benthic bacterial and eukaryote communities along a crude oil spill gradient in a coral reef in the Persian Gulf

Coral reef ecosystems in the Persian Gulf are frequently exposed to crude oil spills. We investigated benthic bacterial and eukaryote community structures at such coral reef sites subjected to different degrees of polycyclic aromatic hydrocarbon (PAH) pollution using environmental DNA (eDNA) metabarcoding. Both bacterial and eukaryote communities responded with pronounced shifts to crude oil pollution and distinguished control sites, moderately and heavily impacted sites with significant confidentiality. The observed community patterns were predominantly driven by Alphaproteobacteria and metazoans. Among these, we identified individual genera that were previously linked to oil spill stress, but also taxa, for which a link to hydrocarbon still remains to be established. Considering the lack of an early-warning system for the environmental status of coral reef ecosystems exposed to frequent crude-oil spills, our results encourage further research towards the development of an eDNA-based biomonitoring tool that exploits benthic bacterial and eukaryote communities as bioindicators.

Environmental nucleic acids: A field-based comparison for monitoring freshwater habitats using eDNA and eRNA

Nucleic acids released by organisms and isolated from environmental substrates are increasingly being used for molecular biomonitoring. While environmental DNA (eDNA) has received much attention, the potential of environmental RNA as a biomonitoring tool remains under-explored. Several recent studies using paired DNA and RNA metabarcoding of bulk samples suggest that RNA might better reflect "metabolically active" parts of the community. However, such studies mainly capture organismal eDNA and eRNA. For larger eukaryotes, isolation of extra-organismal RNA will be important, but viability needs to be examined in a field-based setting. In this study we evaluate (a) whether extra-organismal eRNA release from macroeukaryotes can be detected given its supposedly rapid degradation, and (b) if the same field collection methods for eDNA can be applied to eRNA. We collected eDNA and eRNA from water in lakes where fish community composition is well documented, enabling a comparison between the two nucleic acids in two different seasons with monitoring using conventional methods. We found that eRNA is released from macroeukaryotes and can be filtered from water and metabarcoded in a similar manner as eDNA to reliably provide species composition information. eRNA had a small but significantly greater true positive rate than eDNA, indicating that it correctly detects more species known to exist in the lakes. Given relatively small differences between the two molecules in describing fish community composition, we conclude that if eRNA provides significant advantages in terms of lability, it is a strong candidate to add to the suite of molecular monitoring tools.

Total and active benthic foraminiferal community and their response to heavy metals revealed by high throughput DNA and RNA sequencing in the Zhejiang coastal waters, East China Sea

Benthic foraminifera, large protists abundant in marine environments, have been widely used as bioindicators of environmental conditions. In this study, high-throughput sequencing based on small subunit rDNA and rRNA amplifications was used to investigate total and active benthic foraminifera community composition and diversity from nineteen and twelve superficial marine sediment samples in the Zhejiang coastal waters, respectively. The results showed that the dominant taxa of total foraminifera changed from Buliminellidae (hyaline) to Saccamminidae (agglutinated) from north to south along the coastal waters of Zhejiang Province. According to our survey, heavy metal contamination was moderate in Zhejiang coastal waters, and the potential ecological risks posed by Cd and Hg were higher. The contamination level of heavy metals at Yueqing Bay was the highest, followed by those at Sanmen Bay and Hangzhou Bay. Cd, Cu and grain size may be key factors affecting the distribution and composition of active foraminiferal communities.