Diatom DNA metabarcoding for ecological assessment: Comparison among bioinformatics pipelines used in six European countries reveals the need for standardization

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.

A pile of pipelines: An overview of the bioinformatics software for metabarcoding data analyses

Environmental DNA (eDNA) metabarcoding has gained growing attention as a strategy for monitoring biodiversity in ecology. However, taxa identifications produced through metabarcoding require sophisticated processing of high-throughput sequencing data from taxonomically informative DNA barcodes. Various sets of universal and taxon-specific primers have been developed, extending the usability of metabarcoding across archaea, bacteria and eukaryotes. Accordingly, a multitude of metabarcoding data analysis tools and pipelines have also been developed. Often, several developed workflows are designed to process the same amplicon sequencing data, making it somewhat puzzling to choose one among the plethora of existing pipelines. However, each pipeline has its own specific philosophy, strengths and limitations, which should be considered depending on the aims of any specific study, as well as the bioinformatics expertise of the user. In this review, we outline the input data requirements, supported operating systems and particular attributes of thirty-two amplicon processing pipelines with the goal of helping users to select a pipeline for their metabarcoding projects.

A comparison of two gene regions for assessing community composition of eukaryotic marine microalgae from coastal ecosystems

Two gene regions commonly used to characterise the diversity of eukaryotic communities using metabarcoding are the 18S ribosomal DNA V4 and V9 gene regions. We assessed the effectiveness of these two regions for characterising diverisity of coastal eukaryotic microalgae communities (EMCs) from tropical and temperate sites. We binned amplicon sequence variants (ASVs) into the high level taxonomic groups: dinoflagellates, pennate diatoms, radial centric diatoms, polar centric diatoms, chlorophytes, haptophytes and 'other microalgae'. When V4 and V9 generated ASV abundances were compared, the V9 region generated a higher number of raw reads, captured more diversity from all high level taxonomic groups and was more closely aligned with the community composition determined using light microscopy. The V4 region did resolve more ASVs to a deeper taxonomic resolution within the dinoflagellates, but did not effectively resolve other major taxonomic divisions. When characterising these communities via metabarcoding, the use of multiple gene regions is recommended, but the V9 gene region can be used in isolation to provide high-level community biodiversity to reflect relative abundances within groups. This approach reduces the cost of sequencing multiple gene regions whilst still providing important baseline ecosystem function information.

Maximising environmental pressure-response relationship signals from diatom-based metabarcoding in rivers

DNA metabarcoding has been performed on a large number of river phytobenthos samples collected from the UK, using rbcL primers optimised for diatoms. Within this dataset the composition of non-diatom sequence reads was studied and the effect of including these in models for evaluating the nutrient gradient was assessed. Whilst many non-diatom taxonomic groups were detected, few contained the full diversity expected in riverine environments. This may be due to the performance of the current primers in characterising the wider phytobenthic community and influenced by the sampling method employed, as both were developed specifically for diatoms. Nevertheless, the study identified considerable diversity in some groups, e.g. Eustigmatophyceae and a wider distribution than previously thought for freshwater Phaeophyceae. These results offer a strong case for the benefits of metabarcoding for expanding knowledge of aquatic biodiversity in the UK and elsewhere. Many of the ASVs associated with non-diatoms showed significant pressure responses; however, models that included non-diatoms had similar predictive strength to those based on diatoms alone. Whilst limitations of the primers for assessing non-diatoms may play a role in explaining these results, the diatoms provide a strong signal along the nutrient gradient and other algae, therefore, add little unique information. We recommend that future developments should use ASVs to calculate metrics, with links to reference databases made as a final step to generate lists of taxa to support interpretation. Any further exploration of the potential of non-diatoms would benefit from access to a well-curated reference database, similar to diat.barcode. Such a database does not yet exist, and we caution against the indiscriminate use of NCBI GenBank as a taxonomic resource as many rbcL sequences deposited have not been curated.

A taxonomy-free diatom eDNA-based technique for assessing lake trophic level using lake sediments

Anthropogenic eutrophication is one of the most pressing issues facing lakes globally. Our ability to manage lake eutrophication is hampered by the limited spatial and temporal extents of monitoring records, stemming from the time-consuming and expensive nature of physiochemical and biological monitoring. Diatom-based biomonitoring presents an alternative to traditional eutrophication monitoring, yet it is restricted by the high degree of taxonomic expertise required. Environmental DNA metabarcoding, while providing a promising substitute for diatom community enumeration, is plagued by inadequate taxonomic coverage of reference databases and methodological bias, limiting its use for biomonitoring. Here we show that taxonomy-free diatom-biomonitoring, in which environmental DNA metabarcoding data is utilised but not assigned to specific taxonomic classes, presents an accurate, fast, and relatively automated alternative to taxonomically assigned eutrophication biomonitoring. Our taxonomy-free index accounted for 85% of trophic level variability across 89 lakes and had the lowest average prediction error of the three approaches tested. By not relying on taxonomic identification or metabarcoding reference databases, taxonomy-free biomonitoring maintains diatom diversity that is lost in taxonomic assignment using molecular approaches. Furthermore, by utilising lake sediments, the approach outlined here presents a time-integrated estimation of lake trophic level and thus does not require time-consuming seasonal sampling. Taxonomy-free biomonitoring addresses the limitations of traditional physicochemical eutrophication monitoring and taxonomic biomonitoring alternatives and can be used to extend the spatial and temporal extents of eutrophication monitoring.

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.

Diatom DNA metabarcoding to assess the effect of natural radioactivity in mineral springs on ASV of benthic diatom communities

Little is still known about the low dose effects of radiation on the microbial communities in the environment. Mineral springs are ecosystems than can be affected by natural radioactivity. These extreme environments are, therefore, observatories for studying the influence of chronic radioactivity on the natural biota. In these ecosystems we find diatoms, unicellular microalgae, playing an essential role in the food chain. The present study aimed to investigate, using DNA metabarcoding, the effect of natural radioactivity in two environmental compartments (i.e. spring sediments and water) on the genetic richness, diversity and structure of diatom communities in 16 mineral springs in the Massif Central, France. Diatom biofilms were collected during October 2019, and a 312 bp region of the chloroplast gene rbcL (coding for the Ribulose Bisphosphate Carboxylase) used as a barcode for taxonomic assignation. A total of 565 amplicon sequence variants (ASV) were found. The dominant ASV were associated with Navicula sanctamargaritae, Gedaniella sp., Planothidium frequentissimum, Navicula veneta, Diploneis vacillans, Amphora copulata, Pinnularia brebissonii, Halamphora coffeaeformis, Gomphonema saprophilum, and Nitzschia vitrea, but some of the ASVs could not be assigned at the species level. Pearson correlation failed to show a correlation between ASV' richness and radioactivity parameters. Non-parametric MANOVA analysis based on ASVs occurrence or abundances revealed that geographical location was the main factor influencing ASVs distribution. Interestingly, ²³⁸U was the second factor that explained diatom ASV structure. Among the ASVs in the mineral springs monitored, ASV associated with one of the genetic variants of Planothidium frequentissimum was well represented in the springs and with higher levels of ²³⁸U, suggesting its high tolerance to this particular radionuclide. This diatom species may therefore represent a bio-indicator of high natural levels of uranium.

Comparative analysis of fish environmental DNA reveals higher sensitivity achieved through targeted sequence-based metabarcoding

Environmental DNA (eDNA)-based methods of species detection are enabling various applications in ecology and conservation including large-scale biomonitoring efforts. qPCR is widely used as the standard approach for species-specific detection, often targeting a fish species of interest from aquatic eDNA. However, DNA metabarcoding has the potential to displace qPCR in certain eDNA applications. In this study, we compare the sensitivity of the latest Illumina NovaSeq 6000 NGS platform to qPCR TaqMan assays by measuring limits of detection and by analysing eDNA from water samples collected from Churchill River and Lake Melville, NL, Canada. Species-specific, targeted next generation sequencing (NGS) assays had significantly higher sensitivity than qPCR, with limits of detection 14- to 29-fold lower. For example, when analysing eDNA, qPCR detected Gadus ogac (Greenland cod) in 21% of samples, but targeted NGS detected this species in 29% of samples. General NGS assays were as sensitive as qPCR, while simultaneously detecting 15 fish species from eDNA samples. With over 34,000 fish species on the planet, parallel and sensitive methods such as NGS will be required to support effective biomonitoring at both regional and global scales.

First insights into the gut microbiomes and the diet of the Littorina snail ecotypes, a recently emerged marine evolutionary model

Microbes can play a prominent role in the evolution of their hosts, facilitating adaptation to various environments and promoting ecological divergence. The Wave and Crab ecotypes of the intertidal snail Littorina saxatilis is an evolutionary model of rapid and repeated adaptation to environmental gradients. While patterns of genomic divergence of the Littorina ecotypes along the shore gradients have been extensively studied, their microbiomes have been so far overlooked. The aim of the present study is to start filling this gap by comparing gut microbiome composition of the Wave and Crab ecotypes using metabarcoding approach. Since Littorina snails are micro-grazers feeding on the intertidal biofilm, we also compare biofilm composition (i.e. typical snail diet) in the crab and wave habitats. In the results, we found that bacterial and eukaryotic biofilm composition varies between the typical habitats of the ecotypes. Further, the snail gut bacteriome was different from outer environments, being dominated by Gammaproteobacteria, Fusobacteria, Bacteroidia and Alphaproteobacteria. There were clear differences in the gut bacterial communities between the Crab and the Wave ecotypes as well as between the Wave ecotype snails from the low and high shores. These differences were both observed in the abundances and in the presence of different bacteria, as well as at different taxonomic level, from bacterial OTU's to families. Altogether, our first insights show that Littorina snails and their associated bacteria are a promising marine system to study co-evolution of the microbes and their hosts, which can help us to predict the future for wild species in the face of rapidly changing marine environments.

Environmental DNA study on aquatic ecosystem monitoring and management: Recent advances and prospects

Environmental DNA (eDNA) is organismal DNA that can be detected in the environment and is derived from cellular material of organisms shed into aquatic or terrestrial environments. It can be sampled and monitored using molecular methods, which is important for the early detection of invasive and native species as well as the discovery of rare and cryptic species. While few reviews have summarized the latest findings on eDNA for most aquatic animal categories in the aquatic ecosystem, especially for aquatic eDNA processing and application. In the present review, we first performed a bibliometric network analysis of eDNA studies on aquatic animals. Subsequently, we summarized the abiotic and biotic factors affecting aquatic eDNA occurrence. We also systematically discussed the relevant experiments and analyses of aquatic eDNA from various aquatic organisms, including fish, molluscans, crustaceans, amphibians, and reptiles. Subsequently, we discussed the major achievements of eDNA application in studies on the aquatic ecosystem and environment. The application of eDNA will provide an entirely new paradigm for biodiversity conservation, environment monitoring, and aquatic species management at a global scale.

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.

Evaluation of two short overlapping rbcL markers for diatom metabarcoding of environmental samples: Effects on biomonitoring assessment and species resolution

Two short diatom rbcL barcodes, 331 bp and 263 bp in length, have frequently been used in diatom metabarcoding studies. They overlap in a common 263-bp region but differ in the presence or absence of a 68-bp tail at the 5' end. Though the effectiveness of both has been demonstrated in separate biomonitoring and diversity studies, the impact of the 68-bp non-shared region has not been evaluated. Here we compare the two barcodes in terms of the values of a biotic index (IPS) and the ecological status classes derived from their application to an extensive metabarcoding dataset from United Kingdom rivers; this comprised 1703 samples and was produced using the 331-bp primers. In addition, we assess the effectiveness of each barcode for discrimination of genetic variants around and below the species level. The strong correlation found in IPS values between barcodes (Pearson's R = 0.98) indicates that the choice of the barcode does not have major implications for current WFD ecological assessments, although a very few sites (55: 3.23% of those analysed) were downgraded from an acceptable WFD class ("Good") to an unacceptable one ("Moderate"). Analyses of the taxonomic resolution of the two barcodes indicate that for many ASVs, the use of either marker - 263-bp and 331-bp - gives unambiguous assignations at species level though with differences in bootstrap confidence values. Such differences are caused by the stochasticity involved in the naïve Bayesian classifier used and by the fact that genetic distance, regarding closely related species, is increased when using the 331-bp barcode. However, in three cases, species differentiation fails with the shorter marker, leading to underestimates of species diversity. Finally, two ASVs from Nitzschia species evidenced that the use of the shorter marker can sometimes lead to false positives when the extent and nature of infraspecific variation are poorly known.

Phytoplankton community and HAB species in the South China Sea detected by morphological and metabarcoding approaches

The southern Chinese coast is one of the most developed regions in China and is an area where harmful algal blooms (HABs) have occurred frequently. In this study, differences in the phytoplankton community between microscopic observations and 18S rDNA metabarcoding were compared in 89 surface water samples collected from the southern Chinese coast and the western South China Sea (SCS). This is the first report investigating the phytoplankton community and HAB species using a combination of morphological and metabarcoding approaches in this sea area. There were substantial differences in phytoplankton community structure detected by the two methods. Microscopic observation revealed diatom predominance in the phytoplankton community, while metabarcoding indicated dinoflagellate dominance. The phytoplankton community structure obtained by microscopic observation better reflects the real situation in the water column. Metabarcoding annotated more species than morphospecies observed by microscopy. Haptophyta and Cryptophyta were the specific phyla detected in metabarcoding but were missed in microscopy due to their small size. Conversely, some taxa were found in microscopic analysis alone, such as species in Dinophysis, Prorocentrum, and Scrippsiella, suggesting some biases during metabarcoding and gaps in sequence databases. Metabarcoding is superior for detecting morphologically cryptic, small-sized and HAB taxa, such as unarmored dinoflagellates, nanosized hatophytes and chlorophytes, as well as multiple species in Alexandrium, Pseudonitzschia, and Chaetoceros in our study. A total of 62 HAB taxa were identified in this study, including blooming and potentially toxic species. Diatom abundances generally decreased southward, while those of dinoflagellates and haptophytes showed the opposite trend. Chlorophytes were mainly distributed in coastal waters, especially in the Pearl River Estuary. Phytoplankton community structures were shaped by nutrients and salinity, and phosphorus was the most limiting factor for phytoplankton growth. The phytoplankton community in the western SCS showed unique characteristics away from those in the coastal sea areas. The results suggest that the combination of morphological and metabarcoding approaches comprehensively reveals the phytoplankton community structure and diversity of HAB species.

Diatom eDNA metabarcoding and morphological methods for bioassessment of karstic river

Karst ecosystems play a unique role as exceptional natural habitats in sustaining biodiversity. This study focuses on diatoms, a diverse group of microeukaryotes in the periphytic community of a karstic river. In a multi-microhabitat study along the Krka River (Croatia), our goal was to obtain a detailed overview of diatom diversity and community structure using morphological and molecular approaches, and to assess the applicability of eDNA metabarcoding as a reliable tool for biomonitoring assessment. The results revealed a relatively low agreement in the diatom community composition between the two approaches, but also provided complementary information, with no differences in beta diversity detected between microhabitats. The SIMPER analysis underlined the importance of the molecular approach in identifying diatom community composition, due to errors in distinguishing between deposited diatom cells that occurred in the morphological analysis. In contrast, the morphological approach indicated a clear diatom community separation along the river with a strong location effect. Despite certain differences, both approaches provided a feasible assessment of the ecological status according to the relationship to environmental pressures, classifying the Krka River as High (morphological approach) or Good (molecular approach) throughout the most of its course. Moreover, diatom diversity based on both approaches provides a reliable dataset applicable in routine monitoring assessment and offers a deeper understanding of the presented ecological status. The incompleteness of a reference database presents one major drawback of the molecular approach, which needs further updating in order to improve routine diatom metabarcoding.

Comparing microscopy and DNA metabarcoding techniques for identifying cyanobacteria assemblages across hundreds of lakes

Accurately identifying the species present in an ecosystem is vital to lake managers and successful bioassessment programs. This is particularly important when monitoring cyanobacteria, as numerous taxa produce toxins and can have major negative impacts on aquatic ecosystems. Increasingly, DNA-based techniques such as metabarcoding are being used for measuring aquatic biodiversity, as they could accelerate processing time, decrease costs and reduce some of the biases associated with traditional light microscopy. Despite the continuing use of traditional microscopy and the growing use of DNA metabarcoding to identify cyanobacteria assemblages, methodological comparisons between the two approaches have rarely been reported from a wide suite of lake types. Here, we compare planktonic cyanobacteria assemblages generated by inverted light microscopy and DNA metabarcoding from a 379-lake dataset spanning a longitudinal and trophic gradient. We found moderate levels of congruence between methods at the broadest taxonomic levels (i.e., Order, RV=0.40, p < 0.0001). This comparison revealed distinct cyanobacteria communities from lakes of different trophic states, with Microcystis, Aphanizomenon and Dolichospermum dominating with both methods in eutrophic and hypereutrophic sites. This finding supports the use of either method when monitoring eutrophication in lake surface waters. The biggest difference between the two methods was the detection of picocyanobacteria, which are typically underestimated by light microscopy. This reveals that the communities generated by each method currently are complementary as opposed to identical and promotes a combined-method strategy when monitoring a range of trophic systems. For example, microscopy can provide measures of cyanobacteria biomass, which are critical data in managing lakes. Going forward, we believe that molecular genetic methods will be increasingly adopted as reference databases are routinely updated with more representative sequences and will improve as cyanobacteria taxonomy is resolved with the increase in available genetic information.

Meta-analysis shows both congruence and complementarity of DNA and eDNA metabarcoding to traditional methods for biological community assessment

DNA metabarcoding is increasingly used for the assessment of aquatic communities, and numerous studies have investigated the consistency of this technique with traditional morpho‐taxonomic approaches. These individual studies have used DNA metabarcoding to assess diversity and community structure of aquatic organisms both in marine and freshwater systems globally over the last decade. However, a systematic analysis of the comparability and effectiveness of DNA‐based community assessment across all of these studies has hitherto been lacking. Here, we performed the first meta‐analysis of available studies comparing traditional methods and DNA metabarcoding to measure and assess biological diversity of key aquatic groups, including plankton, microphytobentos, macroinvertebrates, and fish. Across 215 data sets, we found that DNA metabarcoding provides richness estimates that are globally consistent to those obtained using traditional methods, both at local and regional scale. DNA metabarcoding also generates species inventories that are highly congruent with traditional methods for fish. Contrastingly, species inventories of plankton, microphytobenthos and macroinvertebrates obtained by DNA metabarcoding showed pronounced differences to traditional methods, missing some taxa but at the same time detecting otherwise overseen diversity. The method is generally sufficiently advanced to study the composition of fish communities and replace more invasive traditional methods. For smaller organisms, like macroinvertebrates, plankton and microphytobenthos, DNA metabarcoding may continue to give complementary rather than identical estimates compared to traditional approaches. Systematic and comparable data collection will increase the understanding of different aspects of this complementarity, and increase the effectiveness of the method and adequate interpretation of the results.

New and poorly known "araphid" diatom species (Bacillariophyta) from regions near Lake Titicaca, South America and a discussion on the continued use of morphological characters in "araphid" diatom taxonomy

Based on two Andean Altiplano samples and on light and scanning electron microscopy analyses, we present six new species of "araphid" diatoms in the genus Pseudostaurosira, P.aedes sp. nov., P.frankenae sp. nov., P.heteropolaris sp. nov., P.oblonga sp. nov., P.occulta sp. nov., and P.pulchra sp. nov. Additional data are provided for four other known taxa, Nanofrustulumcataractarum, N.rarissimum, P.sajamaensis and P.vulpina, the latter species corresponding to a stat. nov. based on a variety of P.laucensis. Each taxon is described morphologically and compared with closely related published taxa, using characters such as axial area, virgae, vimines, areolar shape, volae, internal striae depositions, spines, flaps and apical pore fields, which are not usually used for species distinction within the genus. It is our intention that the detailed morphological descriptions of each taxon and the elaborate comparative tables we provide serve as a basis for correction of neo and paleo-databases for the Altiplano to produce a better account of autecological data and ecological change in the region. Some arguments for our continued use of a morphologically based approach are given in the context of rapid environmental degradation in the Andes and the difficulties in applying molecular approaches in countries such as Bolivia.

DNA metabarcoding reveals differences in distribution patterns and ecological preferences among genetic variants within some key freshwater diatom species

Our study evaluates differences in the distribution and ecology of genetic variants within several ecologically important diatom species that are also key for Water Framework Directive monitoring of European rivers: Fistulifera saprophila (FSAP), Achnanthidium minutissimum (ADMI), Nitzschia inconspicua (NINC) and Nitzschia soratensis (NSTS). We used DADA2 to infer amplicon sequence variants (ASVs) of a short rbcL barcode in 531 environmental samples from biomonitoring campaigns in Catalonia and France. ASVs within each species showed different distribution patterns. Threshold Indicator Taxa ANalysis revealed three ecological groupings of ASVs in both ADMI and FSAP. Two of these in each species were separated by opposite responses to calcium and conductivity. Boosted regression trees additionally showed that both variables greatly influenced the occurrence of these groupings. A third grouping in FSAP was characterized by a negative response to total organic carbon and hence was better represented in waters with higher ecological status than the other FSAP ASVs, contrasting with what is generally assumed for the species. In the two Nitzschia species, our analyses confirmed earlier studies: NINC preferred higher levels of calcium and conductivity. Our findings suggest that the broad ecological tolerance of some diatom species results from overlapping preferences among genetic variants, which individually show much more restricted preferences and distributions. This work shows the importance of studying the ecological preferences of genetic variants within species complexes, now possible with DNA metabarcoding. The results will help reveal and understand biogeographical distributions and facilitate the development of more accurate biological indexes for biomonitoring programmes.

Assessing changes in stream macroinvertebrate communities across ecological gradients using morphological versus DNA metabarcoding approaches

Freshwater macroinvertebrates provide valuable indicators for biomonitoring ecosystem change in relation to natural and anthropogenic drivers. DNA metabarcoding is an efficient approach for estimating such indicators, but its results may differ from morphotaxonomic approaches traditionally used in biomonitoring. Here we test the hypothesis that despite differences in the number and identity of taxa recorded, both approaches may retrieve comparable patterns of community change, and detect similar ecological gradients influencing such changes. We compared results obtained with morphological identification at family level of macroinvertebrates collected at 80 streams under a Water Framework Directive biomonitoring program in Portugal, with results obtained with metabarcoding from the ethanol preserving the bulk samples, using either single (COI-M19BR2, 16S-Inse01, 18S-Euka02) or multiple markers. Metabarcoding recorded less families and different communities compared to morphotaxonomy, but community sensitivities to disturbance estimated with the IASPT index were more similar across approaches. Spatial variation in local community metrics and the factors influencing such variation were significantly correlated between morphotaxonomy and metabarcoding. After reducing random noise in the dissimilarity matrices, the spatial variation in community composition was also significantly correlated across methods. A dominant gradient of community change was consistently retrieved, and all methods identified a largely similar set of anthropogenic stressors strongly influencing such gradient. Overall, results confirm our initial hypothesis, suggesting that morphotaxonomy and metabarcoding can estimate consistent spatial patterns of community variation and their main drivers. These results are encouraging for macroinvertebrate biomonitoring using metabarcoding approaches, suggesting that they can be intercalibrated with morphotaxonomic approaches to recover equivalent spatial and temporal gradients of ecological change.

Paleo-diatom composition from Santa Barbara Basin deep-sea sediments: a comparison of 18S-V9 and diat-rbcL metabarcoding vs shotgun metagenomics

Sedimentary ancient DNA (sedaDNA) analyses are increasingly used to reconstruct marine ecosystems. The majority of marine sedaDNA studies use a metabarcoding approach (extraction and analysis of specific DNA fragments of a defined length), targeting short taxonomic marker genes. Promising examples are 18S-V9 rRNA (~121-130 base pairs, bp) and diat-rbcL (76 bp), targeting eukaryotes and diatoms, respectively. However, it remains unknown how 18S-V9 and diat-rbcL derived compositional profiles compare to metagenomic shotgun data, the preferred method for ancient DNA analyses as amplification biases are minimised. We extracted DNA from five Santa Barbara Basin sediment samples (up to ~11 000 years old) and applied both a metabarcoding (18S-V9 rRNA, diat-rbcL) and a metagenomic shotgun approach to (i) compare eukaryote, especially diatom, composition, and (ii) assess sequence length and database related biases. Eukaryote composition differed considerably between shotgun and metabarcoding data, which was related to differences in read lengths (~112 and ~161 bp, respectively), and overamplification of short reads in metabarcoding data. Diatom composition was influenced by reference bias that was exacerbated in metabarcoding data and characterised by increased representation of Chaetoceros, Thalassiosira and Pseudo-nitzschia. Our results are relevant to sedaDNA studies aiming to accurately characterise paleo-ecosystems from either metabarcoding or metagenomic data.

DNA metabarcoding and morphological methods show complementary patterns in the metacommunity organization of lentic epiphytic diatoms

Diatoms are important organisms in freshwater ecosystems due to their position as primary producers and therefore, analyzing their assemblages provides relevant information on ecosystem functioning. Diatoms have historically been identified based on morphological traits, which is time-consuming and requires well-trained specialists. Nevertheless, DNA barcoding offers an alternative approach to overcome some limitations of the morphological method. Here, we assess if both approaches are comparable methods to study patterns and mechanisms (including environmental filtering and dispersal limitation) of epiphytic diatom metacommunities using a comprehensive dataset from 22 Mediterranean ponds at different taxonomic resolutions. We used a fragment of rbcL barcode gene combined with High-Throughput Sequencing to infer diatom community composition. The overall degree of correspondence between both approaches was assessed by Procrustean rotation analysis and Procrustean randomization tests, whereas the role of local environmental variables and geographical distances was studied using a comprehensive combination of BIOENV, Mantel tests and distance-based redundancy analysis. Our results showed a relatively poor correspondence in the compositional variation of diatom metacommunity between both approaches. We speculate that the incompleteness of the reference database and the bioinformatics processing are the biases most likely affecting the molecular approach, whereas the limited counting effort and the presence of cryptic species are presumably the major biases related with the morphological method. On the other hand, variation in diatom community composition detected with both approaches was strongly related to the environmental template, which may be related with the narrow community-environment relationships in diatoms. Nevertheless, we found no significant relationship between compositional variation and geographical distances. Overall, our work shows the complementary nature of both approaches and highlights the importance of DNA metabarcoding to address empirical research questions of community ecology in freshwaters, especially once the reference databases include most genotypes of occurring taxa and bioinformatics biases are overcome.

How can integrated morphotaxonomy- and metabarcoding-based diatom assemblage analyses best contribute to the ecological assessment of streams?

Environmental conditions, such as nutrient concentrations, salinity, elevation etc., shape diatom assemblages of periphytic biofilms. These assemblages respond rapidly to environmental changes, a fact which makes diatoms valuable bioindicators. Hence, freshwater biomonitoring programmes currently use diatom indices (e.g. EU Water Framework Directive - WFD). To date, microscopy-based assessments require high taxonomic expertise for diatom identification at the species level. High-throughput technologies now provide cost-effective identification approaches that are promising, complementary or alternative tools for bioassessment. The suitability of the metabarcoding method is evaluated for the first time in the Cyprus streams WFD monitoring network, an eastern Mediterranean country with many endemic species and results are compared to the results acquired from the morphotaxonomic analysis. Morphotaxonomic identification was conducted microscopically, using the most updated taxonomic concepts, literature and online resources. At the same time, DNA metabarcoding involved the use of the rbcL 312 bp barcode, high-throughput sequencing and bioinformatic analysis. The ecological status was calculated using the IPS Index. Results show a positive correlation between morpho-taxonomic and molecular IPS scores. Discrepancies between the two methodologies are related to the limitations of both techniques. This study confirmed that Fistulifera saprophila can have a crucial role in key differences observed, as it negatively influences IPS scores and microscopy methods frequently overlook it. Importantly, gaps in the DNA barcoding reference databases lead to a positive overestimation in IPS scores. Overall, we conclude that DNA metabarcoding offsets the morphotaxonomic methodology for the ecological quality assessment of freshwaters.

Towards reproducible metabarcoding data: Lessons from an international cross-laboratory experiment

Advances in high-throughput sequencing (HTS) are revolutionizing monitoring in marine environments by enabling rapid, accurate and holistic detection of species within complex biological samples. Research institutions worldwide increasingly employ HTS methods for biodiversity assessments. However, variance in laboratory procedures, analytical workflows and bioinformatic pipelines impede the transferability and comparability of results across research groups. An international experiment was conducted to assess the consistency of metabarcoding results derived from identical samples and primer sets using varying laboratory procedures. Homogenized biofouling samples collected from four coastal locations (Australia, Canada, New Zealand and the USA) were distributed to 12 independent laboratories. Participants were asked to follow one of two HTS library preparation workflows. While DNA extraction, primers and bioinformatic analyses were purposefully standardized to allow comparison, many other technical variables were allowed to vary among laboratories (amplification protocols, type of instrument used, etc.). Despite substantial variation observed in raw results, the primary signal in the data was consistent, with the samples grouping strongly by geographical origin for all data sets. Simple post hoc data clean-up by removing low-quality samples gave the best improvement in sample classification for nuclear 18S rRNA gene data, with an overall 92.81% correct group attribution. For mitochondrial COI gene data, the best classification result (95.58%) was achieved after correction for contamination errors. The identified critical methodological factors that introduced the greatest variability (preservation buffer, sample defrosting, template concentration, DNA polymerase, PCR enhancer) should be of great assistance in standardizing future biodiversity studies using metabarcoding.

Potential for cross-contamination of diatom DNA samples when using toothbrushes

The use of toothbrushes and similar devices for sampling diatoms from hard surfaces is a well-established approach. Toothbrushes are routinely cleaned and reused when sampling for analysis by light microscopy. This paper looks at the scale of contamination encountered when this technique is used to sample diatoms for metabarcoding analyses, as well as at the scale of contamination to be expected if stream, rather than distilled water, is used to wash diatoms from stones. Although some contamination attributable to toothbrushes was detected, read numbers were low and had no effect on index calculation or ecological status estimates. However, if the primary focus of a study is to thoroughly document diversity in a sample, then even this small level of contamination may be unacceptable and more stringent measures may be required.

Epiphytic Diatom-Based Biomonitoring in Mediterranean Ponds: Traditional Microscopy versus Metabarcoding Approaches

Benthic diatoms have traditionally been used as bioindicators of aquatic ecosystems. Because diatom-based monitoring of water quality is required by European legislation, molecular-based methods had emerged as useful alternatives to classical methods based on morphological identification using light microscopy. The aim of this study was to test the reliability of DNA metabarcoding combined with High-Throughput Sequencing (HTS) techniques in the bioassessment of the trophic status of 22 Mediterranean shallow ponds in NW Spain. For each pond, the Trophic Diatom Index (TDI) was calculated from inventories obtained by identification using light microscopy (LM) followed by high-throughput sequencing (HTS) at the molecular level. Ponds were subsequently classified into five water quality classes. The results showed a good correspondence between both methods, especially after applying a correction factor that depended on the biovolume of the cells. This correspondence led to the assignment to the same quality class in 59% of the ponds. The determination and quantification of valves or DNA sequences was one of the main pitfalls, which mainly included those related to the variability in the relative abundances of some species. Accordingly, ponds with similar relative abundances for the dominant species were assigned to the same quality class. Moreover, other difficulties leading the discrepancies were the misidentification of some species due to the presence of semi-cryptic taxa, the incompleteness of the reference database and the bioinformatic protocol. Thus, the validation of DNA-based methods for the identification of freshwater diatoms represents an important goal, as an alternative to using traditional methods in Mediterranean shallow ponds.

Resources and Practices to Improve Diatom Data Quality

Environmental programs in the United States face technical challenges that inhibit the ability to use diatoms in water quality monitoring and assessment projects. Specifically, inconsistent taxonomy can obscure diatom responses to environmental variables. Problems are the result of (1) limited access to a common set of taxonomic references, especially those that are geographically relevant, (2) inefficient enumeration protocols, (3) lack of complete and transparent documentation of taxa, and (4) limited opportunities for continued education, training, and knowledge sharing. However, robust resources and practices are available to improve diatom data quality and interpretation. Several resources improve diatom data quality, including a publicly accessible taxonomic reference (diatoms.org) and recommended practices. These practices include adoption of the voucher floras, random sample assignment, replicate microscope slides, and improved quality control. Finally, the Society for Freshwater Science Diatom Taxonomic Certification Committee is developing educational materials and certification exams to support practitioner training and to increase the diatom research knowledge base. The resources and practices in this article are broadly applicable to improving basic and applied research on diatoms worldwide.

Expanding ecological assessment by integrating microorganisms into routine freshwater biomonitoring

Bioindication has become an indispensable part of water quality monitoring in most countries of the world, with the presence and abundance of bioindicator taxa, mostly multicellular eukaryotes, used for biotic indices. In contrast, microbes (bacteria, archaea and protists) are seldom used as bioindicators in routine assessments, although they have been recognized for their importance in environmental processes. Recently, the use of molecular methods has revealed unexpected diversity within known functional groups and novel metabolic pathways that are particularly important in energy and nutrient cycling. In various habitats, microbial communities respond to eutrophication, metals, and natural or anthropogenic organic pollutants through changes in diversity and function. In this review, we evaluated the common trends in these changes, documenting that they have value as bioindicators and can be used not only for monitoring but also for improving our understanding of the major processes in lotic and lentic environments. Current knowledge provides a solid foundation for exploiting microbial taxa, community structures and diversity, as well as functional genes, in novel monitoring programs. These microbial community measures can also be combined into biotic indices, improving the resolution of individual bioindicators. Here, we assess particular molecular approaches complemented by advanced bioinformatic analysis, as these are the most promising with respect to detailed bioindication value. We conclude that microbial community dynamics are a missing link important for our understanding of rapid changes in the structure and function of aquatic ecosystems, and should be addressed in the future environmental monitoring of freshwater ecosystems.

Lake Sedimentary DNA Research on Past Terrestrial and Aquatic Biodiversity: Overview and Recommendations

The use of lake sedimentary DNA to track the long-term changes in both terrestrial and aquatic biota is a rapidly advancing field in paleoecological research. Although largely applied nowadays, knowledge gaps remain in this field and there is therefore still research to be conducted to ensure the reliability of the sedimentary DNA signal. Building on the most recent literature and seven original case studies, we synthesize the state-of-the-art analytical procedures for effective sampling, extraction, amplification, quantification and/or generation of DNA inventories from sedimentary ancient DNA (sedaDNA) via high-throughput sequencing technologies. We provide recommendations based on current knowledge and best practises.

Host-parasitoid associations in marine planktonic time series: Can metabarcoding help reveal them?

In this study, we created a dataset of a continuous three-year 18S metabarcoding survey to identify eukaryotic parasitoids, and potential connections to hosts at the Long-Term Ecological Research station Helgoland Roads. The importance of parasites and parasitoids for food web dynamics has previously been recognized mostly in terrestrial and freshwater systems, while marine planktonic parasitoids have been understudied in comparison to those. Therefore, the occurrence and role of parasites and parasitoids remains mostly unconsidered in the marine environment. We observed high abundances and diversity of parasitoid operational taxonomic units in our dataset all year round. While some parasitoid groups were present throughout the year and merely fluctuated in abundances, we also detected a succession of parasitoid groups with peaks of individual species only during certain seasons. Using co-occurrence and patterns of seasonal occurrence, we were able to identify known host-parasitoid dynamics, however identification of new potential host-parasitoid interactions was not possible due to their high dynamics and variability in the dataset.