Diat.barcode, an open-access curated barcode library for diatoms

Biogeography and genetic diversity of freshwater diatoms: The potential of large combined rbcL metabarcoding datasets

Aiming to gain a general picture of rbcL diversity within freshwater diatom species, this study assembles and analyzes multiple metabarcoding datasets spanning various geographical regions. From these datasets, we inferred >10,000 amplicon sequence variants (ASVs) of 263-bp length. More than half of the 1000 most abundant ASVs were recorded in both Eurasia and N America and there was only limited evidence for continent-specific lineages. The geographical range was extended for some species, illustrating the potential of metabarcoding datasets for such checks. For detailed analysis of intraspecific diversity, 73 freshwater species were selected, corresponding to 360 ASVs assigned phylogenetically. We found notable variation, some species being represented by only one or a few ASVs, while others were represented by a higher number. Furthermore, within species, ASVs exhibited different dominance and distribution patterns, in some cases with a head-tail pattern, in others a more equal spread of abundance or unresolved reticulate relationships. Except for Ulnaria ulna, no geographical structure among species' ASVs was detectable in haplotype networks using the 263-bp rbcL marker. Observed heterogeneity within species was categorized by computing several metrics of genetic variation and classified into three groups, reflecting optimal sampling strategies based on the patterns of intraspecific variation in the 73 target species There was a significant relationship between intraspecific diversity and the traditional separation between 'centric' and 'pennate' diatoms, with centric species exhibiting significantly fewer variants than pennates, possibly because of different plastid inheritance patterns.

Enhancing microbial predator-prey detection with network and trait-based analyses

BACKGROUND

Network analyses are often applied to microbial communities using sequencing survey datasets. However, associations in such networks do not necessarily indicate actual biotic interactions, and even if they do, the nature of the interactions commonly remains unclear. While network analyses are valuable for generating hypotheses, the inferred hypotheses are rarely experimentally confirmed.

RESULTS

We employed cross-kingdom network analyses, applied trait-based functions to the microorganisms, and subsequently experimentally investigated the found putative predator-prey interactions to evaluate whether, and to what extent, correlations indicate actual predator-prey relationships. For this, we investigated algae and their protistan predators in biocrusts of three distinct polar regions, i.e., Svalbard, the Antarctic Peninsula, and Continental Antarctica. Network analyses using FlashWeave indicated that 89, 138, and 51 correlations occurred between predatory protists and algae, respectively. However, trait assignment revealed that only 4.7-9.3% of said correlations link predators to actually suitable prey. We further confirmed these results with HMSC modeling, which resulted in similar numbers of 7.5% and 4.8% linking predators to suitable prey for full co-occurrence and abundance models, respectively. The combination of network analyses and trait assignment increased confidence in the prediction of predator-prey interactions, as we show that 82% of all experimentally investigated correlations could be verified. Furthermore, we found that more vicious predators, i.e., predators with the highest growth rate in co-culture with their prey, exhibit higher stress and betweenness centrality - giving rise to the future possibility of determining important predators from their network statistics.

CONCLUSIONS

Our results support the idea of using network analyses for inferring predator-prey interactions, but at the same time call for cautionary consideration of the results, by combining them with trait-based approaches to increase confidence in the prediction of biological interactions. Video Abstract.

Sexual reproduction during diatom bloom

Phytoplankton supports food webs in all aquatic ecosystems. Ecological studies highlighted the links between environmental variables and species successions in situ. However, the role of life cycle characteristics on phytoplankton community dynamics remains poorly characterized. In diatoms, sexual reproduction creates new genetic combinations and prevents excessive cell size miniaturization. It has been extensively studied in vitro but seldom in the natural environment. Here, analyzing metatranscriptomic data in the light of the expression patterns previously characterized in vitro, we identified a synchronized and transient sexual reproduction event during a bloom of the toxic diatom species Pseudo-nitzschia australis. Despite the complexity of environmental conditions encountered in situ, sexual reproduction appeared to be the strongest differential gene expression signal that occurred during the bloom. The potential link between environmental conditions and the initiation of sexual reproduction remain to be determined, but sexual reproduction probably had a major impact on the bloom dynamic.

Could landscape ecology principles apply at the microscale? A metabarcoding approach on Trichoptera larvae-associated microbial diversity

Landscape heterogeneity is known as a major factor of community structure and composition. Whether this effect of the landscape extends at different scales and particularly at the relevant scale for microorganisms remained to be determined. We used the cases produced by aquatic larvae of Trichoptera, which assemble organic or mineral particles, as naturally replicated experimental systems representing structured substrates to determine the effect of landscape structuration on microbial communities. A metabarcoding approach was used to characterise fungal, bacterial and diatom communities on cases produced by six Trichoptera species and related unstructured organic and mineral substrates. The structuration of the particles constituting the cases was also determined as a measure of microscale landscape. Structured substrates harboured communities of diatoms, fungi and bacteria that differed from those found on unstructured substrates. Microbial communities also differed between organic and mineral substrates. We found a higher microbial diversity on structured substrates than on unstructured substrates. The heterogeneity of the microscale landscape also affected bacterial and fungal communities within cases. These results highlight the importance of microscale landscape structuration for microbial diversity and demonstrate that approaches of landscape ecology could be downscaled to the microscale.

The importance of integrating phycological research, teaching, outreach, and engagement in a changing world

The ecological, evolutionary, economic, and cultural importance of algae necessitates a continued integration of phycological research, education, outreach, and engagement. Here, we comment on several topics discussed during a networking workshop-Algae and the Environment-that brought together phycological researchers from a variety of institutions and career stages. We share some of our perspectives on the state of phycology by examining gaps in teaching and research. We identify action areas where we urge the phycological community to prepare itself to embrace the rapidly changing world. We emphasize the need for more trained taxonomists as well as integration with molecular techniques, which may be expensive and complicated but are important. An essential benefit of these integrative studies is the creation of high-quality algal reference barcoding libraries augmented with morphological, physiological, and ecological data that are important for studies of systematics and crucial for the accuracy of the metabarcoding bioassessment. We highlight different teaching approaches for engaging undergraduate students in algal studies and the importance of algal field courses, forays, and professional phycological societies in supporting the algal training of students, professionals, and citizen scientists.

BFR2: a curated ribosomal reference dataset for benthic foraminifera

Benthic foraminifera are one of the major groups of marine protists that also occur in freshwater and terrestrial habitats. They are widely used to monitor current and past environmental conditions. Over the last three decades, thousands of DNA sequences have been obtained from benthic foraminiferal isolates. The results of this long-term effort are compiled here in the form of the first curated benthic foraminiferal ribosomal reference dataset (BFR2). The present dataset contains over 5000 sequences of a fragment of the 18S rDNA gene, which is recognized as the DNA barcode of foraminifera. The sequences represent 279 species and 204 genera belonging to 91 families. Thirteen percent of these sequences have not been assigned to any morphologically described group and may represent species new to science. Furthermore, forty-five percent of the sequences have not been previously published. The BFR2 dataset aims to collect all DNA barcodes of benthic foraminifera and to provide a much-needed reference dataset for the rapidly developing field of molecular foraminiferal studies.

The Diatom Diversity and Ecological Status of a Tufa-Depositing River through eDNA Metabarcoding vs. a Morphological Approach-A Case Study of the Una River (Bosnia and Herzegovina)

Tufa deposits in karst rivers are unique habitats created by mutual interactions between specific environmental and biotope features and inhabited by diatoms as a highly abundant and diverse algal group. This pilot study aimed to investigate the diversity of diatom communities on tufa depositing habitats and assess the Una River's ecological status using a comparative molecular and morphological approach for diatom identification. The 312 base pairs of the rbcL gene were barcoded and analyzed using MiSeq reads and amplicon sequence variants (ASVs) obtained by the DADA2 pipeline. The reference database Diat.barcode v7 was used for taxonomic assignment. The morphological identification of the diatoms was carried out in parallel. In total, the combined dataset revealed 46 taxa identified at genus rank, 125 on the subgenus, and 145 on combined taxonomy rank. The metabarcoding approach mostly leads to a lower number of identified taxa at species rank (58 in molecular vs. 119 in optical inventory), resulting in higher values of beta diversity and heterogeneity in diatom assemblages in samples obtained by morphological approach. Despite the high percentage of taxonomically not assigned diatom ASVs to the species rank, high Shannon diversity index values and a similar number of taxa per locations compared to the morphological approach were obtained. Taxa Achnanthidium minutissimum (Kützing) Czarnecki, Achnanthidium pyrenaicum (Hustedt) H.Kobayasi, Amphora pediculus (Kützing) Grunow, Diatoma vulgaris Bory, Navicula cryptotenella Lange-Bertalot, and Navicula tripunctata (O.F.Müller) Bory were identified at all locations in both inventories. Although limited consistency in the diatom abundances between the two inventory datasets was found, a similar grouping of samples was observed connected to the river's longitudinal gradient. The data obtained using molecular approach in most sites indicated a mostly lower ecological status (good or moderate) compared to the data obtained from the morphological approach (high, good, and moderate). The potential of environmental DNA (eDNA) diatom metabarcoding for water monitoring and diversity studies is undeniable, but to fully realize the benefits of these methods in the future, it is essential to standardize protocols and expand the reference database for species found in specific habitats, such as tufa deposits.

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.

First regional reference database of northern Adriatic diatom transcriptomes

Marine microbial communities form the basis for the functioning of marine ecosystems and the conservation of biodiversity. With the application of metagenomics and metatranscriptomics in marine environmental studies, significant progress has been made in analysing the functioning of microbial communities as a whole. These molecular techniques are highly dependent on reliable, well-characterised, comprehensive and taxonomically diverse sequenced reference transcriptomes of microbial organisms. Here we present a set of 12 individual transcriptome assemblies derived from 6 representative diatom species from the northern Adriatic Sea grown under 2 environmentally relevant growth conditions (phosphate replete vs. phosphate deprived). After filtering the reads and assembly, an average number of 64,932 transcripts per assembly was obtained, of which an average of 8856 were assigned to functionally known proteins. Of all assigned transcripts, an average of 6483 proteins were taxonomically assigned to diatoms (Bacillariophyta). On average, a higher number of assigned proteins was detected in the transcriptome assemblies of diatoms grown under replete media condition. On average, 50% of the mapped proteins were shared between the two growth conditions. All recorded proteins in the dataset were classified into 24 COG categories, with approximately 25% belonging to the unknown function and the remaining 75% belonging to all other categories. The resulting diatom reference database for the northern Adriatic, focussing on the response to nutrient limitation as characteristic for the region and predicted for the future world oceans, provides a valuable resource for analysing environmental metatranscriptome and metagenome data. Each northern Adriatic transcriptome can also be used by itself as a reference database for the (meta)transcriptomes and gene expression studies of the associated species that will be generated in the future.

Distinct multitrophic biodiversity composition and community organization in a freshwater lake and a hypersaline lake on the Tibetan Plateau

Alpine lakes play pivotal roles in plateau hydrological processes but are highly sensitive to climate change, yet we lack comprehensive knowledge of their multitrophic biodiversity patterns. Here, we compared the biodiversity characteristics of diverse taxonomic groups across water depths and in surface sediments from a freshwater lake and a hypersaline lake on the northwestern Tibetan Plateau. Using multi-marker environmental DNA metabarcoding, we detected 134 cyanobacteria, 443 diatom, 1,519 invertebrate, and 28 vertebrate taxa. Each group had a substantially different community composition in the two lakes, and differences were also found between water and sediments within each lake. Cooccurrence network analysis revealed higher network complexity, lower modularity, and fewer negative cohesions in the hypersaline lake, suggesting that high salinity may destabilize ecological networks. Our results provide the first holistic view of Tibetan lake biodiversity under contrasting salinity levels and reveal structural differences in the ecological networks that may impact ecosystem resilience.

Successive accumulation of biotic assemblages at a fine spatial scale along glacier-fed waters

Glacier-fed waters create strong environmental filtering for biota, whereby different organisms may assume distinct distribution patterns. By using environmental DNA-based metabarcoding, we investigated the multi-group biodiversity distribution patterns of the Parlung No. 4 Glacier, on the Tibetan Plateau. Altogether, 642 taxa were identified from the meltwater stream and the downstream Ranwu Lake, including 125 cyanobacteria, 316 diatom, 183 invertebrate, and 18 vertebrate taxa. As the distance increased from the glacier terminus, community complexity increased via sequential occurrences of cyanobacteria, diatoms, invertebrates, and vertebrates, as well as increasing taxa numbers. The stream and lake showed different community compositions and distinct taxa. Furthermore, the correlations with environmental factors and community assembly mechanisms showed group- and habitat-specific patterns. Our results reveal the rapid spatial succession and increasing community complexity along glacial flowpaths and highlight the varying adaptivity of different organisms, while also providing insight into the ecosystem responses to global change.

The chloroplast and mitochondrial genomes of two Gomphonema parvulum (Bacillariophyta) environmental isolates from South Carolina (United States) and Virginia (United States)

Gomphonema parvulum is a cosmopolitan freshwater diatom that is used as an indicator in water quality biomonitoring. In this study, we report the culturing of two geographically separated isolates from southeastern North America, their morphology, and the sequencing and assembly of their mitochondrial and chloroplast genomes. Morphologically, both strains fit G. parvulum sensu lato, but the frustules from a protected habitat in South Carolina were smaller than those cited in the historic data of this species from the same location as well as a second culture from Virginia. Phylogenetic analyses using the rbcL gene placed both within a clade with G. parvulum. Genetic markers, including full chloroplast and mitochondrial genomes and the nuclear small subunit rRNA gene region were assembled from each isolate. The organellar genomes of the two strains varied slightly in size due to small differences in intergenic regions with chloroplast genomes of 121,035 bp and 121,482 bp and mitochondrial genomes of 34,639 bp and 34,654 bp. The intraspecific pairwise identities of the chloroplast and mitochondrial genomes of these two isolates were 97.9% and 95.4%, respectively. Multigene phylogenetic analysis demonstrated a close relationship between G. parvulum, Gomphoneis minuta, and Didymosphenia geminata.

Assessing the relevance of DNA metabarcoding compared to morphological identification for lake phytoplankton monitoring

Phytoplankton is a key biological group used to assess the ecological status of lakes. The classical monitoring approach relies on microscopic identification and counting of phytoplankton species, which is time-consuming and requires high taxonomic expertise. High-throughput sequencing, combined with metabarcoding, has recently demonstrated its potential as an alternative approach for plankton surveys. Several studies have confirmed the relevance of the diatom metabarcoding approach to calculate biotic indices based on species ecology. However, phytoplankton communities have not yet benefited from such validation. Here, by comparing the results obtained with the two methods (molecular and microscopic counting), we evaluated the relevance of metabarcoding approach for phytoplankton monitoring by considering different metrics: alpha diversity, taxonomic composition, community structure and a phytoplankton biotic index used to assess the trophic level of lakes. For this purpose, 55 samples were collected in four large alpine lakes (Aiguebelette, Annecy, Bourget, Geneva) during the year 2021. For each sample, a metabarcoding analysis based on two genetic markers (16S and 23S rRNA) was performed, in addition to the microscopic count. Regarding the trophic level of lakes, significant differences were found between index values obtained with the two approaches. The main hypothesis to explain these differences comes from the incompleteness, particularly at the species level, of the barcode reference library for the two genetic markers. It is therefore necessary to complete reference libraries for using such species-based biotic indices with metabarcoding data. Besides this, species richness and diversity were higher in the molecular inventories than in the microscopic ones. Moreover, despite differences in taxonomic composition of the floristic lists obtained by the two approaches, their community structures were similar. These results support the possibility of using metabarcoding for phytoplankton monitoring but in a different way. We suggest exploring alternative approaches to index development, such as a taxonomy-free approach.

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.

Metabarcoding reveals unique microbial mat communities and evidence of biogeographic influence in low-oxygen, high-sulfur sinkholes and springs

High-sulfur, low-oxygen environments formed by underwater sinkholes and springs create unique habitats populated by microbial mat communities. To explore the diversity and biogeography of these mats, samples were collected from three sites in Alpena, Michigan, one site in Monroe, Michigan, and one site in Palm Coast, Florida. Our study investigated previously undescribed eukaryotic diversity in these habitats and further explored their bacterial communities. Mat samples and water parameters were collected from sulfur spring sites during the spring, summer, and fall of 2022. Cyanobacteria and diatoms were cultured from mat subsamples to create a culture-based DNA reference library. Remaining mat samples were used for metabarcoding of the 16S and rbcL regions to explore bacterial and diatom diversity, respectively. Analyses of water chemistry, alpha diversity, and beta diversity articulated a range of high-sulfur, low-oxygen habitats, each with distinct microbial communities. Conductivity, pH, dissolved oxygen, temperature, sulfate, and chloride had significant influences on community composition but did not describe the differences between communities well. Chloride concentration had the strongest correlation with microbial community structure. Mantel tests revealed that biogeography contributed to differences between communities as well. Our results provide novel information on microbial mat composition and present evidence that both local conditions and biogeography influence these unique communities.

Alien Algae Species Invasions in Humic Rivers within Weakly Human Impact Basin

Increasing anthropogenic influence and climate change are leading to significant transformations in living conditions for indigenous representatives of aquatic communities. The problem of alien species invasions is actively discussed in the example of large rivers and water reservoirs, but medium and small rivers with weak anthropogenic influence have been insufficiently studied in this aspect. With the help of analysis of literary data and the results of our own long-term observations of phytoplankton using morphological and molecular genetic methods in some left-bank Volga tributaries, we identified six invasive species of different taxonomic groups of algae, with a predominance of diatoms. The relevance of using both traditional and modern approaches to identifying invasive algae species is revealed. Such invasive species as Thalassiosira incerta, T. faurii, Skeletonema subsalsum, Unruhdinium kevei, and Gonyostomum semen were part of planktonic communities; the benthic species Plagiotropis lepidoptera var. proboscidea sometimes reached a significant level of development (up 6 to 44% from total biomass) in plankton. It was demonstrated that some algae species have firmly taken the position of dominants and subdominants in planktonic algae communities. The expansion of alien representatives of microphytobenthos was noted in the Volga River basin for the first time. For Gonyostomum semen, its European origin was revealed, for plankton and benthic diatom-Ponto-Caspian. Our study showed that the processes of invasion and subsequent development of alien species take place in habitats with weak anthropogenic influence, which is likely determined by the hydrological, hydrochemical, and climatic changes in river basins and the high adaptive capabilities of invasive species.

100 years of anthropogenic impact causes changes in freshwater functional biodiversity

Despite efforts from scientists and regulators, biodiversity is declining at an alarming rate. Unless we find transformative solutions to preserve biodiversity, future generations may not be able to enjoy nature's services. We have developed a conceptual framework that establishes the links between biodiversity dynamics and abiotic change through time and space using artificial intelligence. Here, we apply this framework to a freshwater ecosystem with a known history of human impact and study 100 years of community-level biodiversity, climate change and chemical pollution trends. We apply explainable network models with multimodal learning to community-level functional biodiversity measured with multilocus metabarcoding, to establish correlations with biocides and climate change records. We observed that the freshwater community assemblage and functionality changed over time without returning to its original state, even if the lake partially recovered in recent times. Insecticides and fungicides, combined with extreme temperature events and precipitation, explained up to 90% of the functional biodiversity changes. The community-level biodiversity approach used here reliably explained freshwater ecosystem shifts. These shifts were not observed when using traditional quality indices (e.g. Trophic Diatom Index). Our study advocates the use of high-throughput systemic approaches on long-term trends over species-focused ecological surveys to identify the environmental factors that cause loss of biodiversity and disrupt ecosystem functions.

eDNA-based diversity and multitrophic network reveal the effects of land use and pollutants on the subtropical Dongjiang River systems

Biodiversity and its constituted multitrophic network in rivers are accelerating change under human land use and pollutants. However, due to the lack of complete datasets across taxa limited by traditional morphological biomonitoring, the change patterns of biodiversity and multitrophic networks are still unclear. Here, we used the eDNA approach to capture multitrophic communities (including fish, aquatic insects, protozoa, diatom and bacteria) in the Dongjiang River, a typical subtropical river in southeast China, and analyzed the changing patterns of biodiversity and multitrophic networks in relation to land use and water pollution. First, our data showed that the eDNA approach provided a snapshot of the multitrophic communities in the Dongjiang River, and the monitored 5833 OTUs were annotated to 55 phyla, 144 classes, 329 orders, 521 families, 945 genera and 406 species. Second, the multitrophic diversity index had similar patterns on the longitudinal scale of rivers, with significant decreases from the upstream to the downstream, while individual taxonomic groups exhibited variable spatial patterns. While there were similar spatial patterns between network metrics and diversity index, the former had stronger relationships with the spatial distance. Third, the multitrophic diversity and networks were significantly negatively correlated with land use and water pollution (e.g., CODMn), and network structures often had stronger and non-linear responses. Overall, this study highlights that eDNA biomonitoring of multitrophic communities and networks can provide deeper insights into ecosystem changes and help develop more targeted management strategies.

A metabarcode based (species) inventory of the northern Adriatic phytoplankton

Background

The northern Adriatic is characterised as the coldest and most productive marine area of the Mediterranean, which is due to high nutrient levels introduced by river discharges, the largest of which is the Italian Po River (at the same time also the largest freshwater input into the Mediterranean). The northern Adriatic is a very shallow marine ecosystem with ocean current patterns that result in long retention times of plankton in the area. The northern Adriatic phytoplankton biodiversity and abundance are well-studied, through many scientific and long-term monitoring reports. These datasets were based on phytoplankton morphological traits traditionally obtained with light microscopy. The most recent comprehensive eastern Adriatic phytoplankton checklist was published more than 20 years ago and is still valuable today. Since phytoplankton taxonomy and systematics are constantly being reviewed (partly also due to new molecular methods of species identification that complement classical methodologies), checklists need to be updated and complemented. Today, metabarcoding of molecular markers gains more and more importance in biodiversity research and monitoring. Here, we report the use of high throughput sequencing methods to re-examine taxonomic richness and provide updated knowledge of phytoplankton diversity in the eastern northern Adriatic to complement the standardised light microscopy method.

New information

This study aimed to report an up-to-date list of the phytoplankton taxonomic richness and phylogenetic relationships in the eastern northern Adriatic, based on sequence variability of barcoding genes resolved with advanced molecular tools, namely metabarcoding. Here, metabarcoding is used to complement standardised light microscopy to advance conventional monitoring and research of phytoplankton communities for the purpose of assessing biodiversity and the status of the marine environments. Monthly two-year net sampling targeted six phytoplankton groups including Bacillariophyceae (diatoms) and Chrysophyceae (golden algae) belonging to Ochrophyta, Dinophyceae (dinoflagellates), Cryptophyceae (cryptophytes), Haptophyta (mostly coccolithophorids) and Chlorophyta with Prasinophyceae (prasinophytes) and Chlorophyceae (protist green algae). Generated sequence data were taxonomically assigned and redistributed in two kingdoms, five classes, 32 orders, 49 families and 67 genera. The most diverse group were dinoflagellates, comprising of 34 found genera (48.3%), following by diatoms with 23 (35.4%) and coccolithophorids with three genera (4.0%). In terms of genetic diversity, results were a bit different: a great majority of sequences with one nucleotide tolerance (ASVs, Amplicon sequence variants) assigned to species or genus level were dinoflagellates (83.8%), 13.7% diatoms and 1.6% Chlorophyta, respectively. Although many taxa have not been detected that have been considered as common in this area, metabarcoding revealed five diatoms and 20 dinoflagellate genera that were not reported in previous checklists, along with a few species from other targeted groups that have been reported previously. We here describe the first comprehensive 18S metabarcode inventory for the northern Adriatic Sea.

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 endemism and taxonomic turnover: Assessment in high-altitude alpine lakes covering a large geographical range

Diatoms are widely used as ecological indicators and show various degrees of endemism. Many studies that support the idea of endemic species integrate several climate zones, a variety of ecosystem types, and often focus on a global scale. Here, we investigated whether endemism could be detected when considering a homogeneous type of ecosystem in a single climate zone. We sampled stone biofilms at 40-50 cm depth in high-altitude lakes in the Alpine climate zone. A total of 149 samples were obtained from the French and Georgian mountains, two areas separated by ∼3000 km. Using Amplicon Sequence Variants derived from DNA metabarcoding, we assessed taxonomic turnover and Zeta-diversity (a measure of endemism). We ran haplotype networks and phylogenetic tests to measure geographical signal in the phylogenies of dominant taxa. The French and Georgian communities shared 51 % of species. Species that were not shared across both regions were mostly rare, and often not characteristic of lakes but of neighboring habitats instead. In contrast, at the sub-species level, 87 % of the genotypes showed restricted distributions. Whereas endemism was the rule at sub-species level, most species were shared across both French and Georgian lakes, suggesting that geographic barriers strongly limited dispersal at the sub-species level but not species level. Dominant species hosted higher levels of sub-specific diversity than rare species. In contrast to global-scale studies, we did not find any significant geographical structuring in the phylogeny of the investigated species. This could indicate ongoing dispersal at a frequency fast enough to prevent allopatric divergence, yet slow enough to prevent sharing most haplotypes between France and Georgia. These results have implications for biomonitoring: depending on the taxonomic level chosen, robust generic tools (species level) or tools dedicated to a region able to discriminate fine pressures differences (sub-species level) may be developed.

The Roles of Microbes in Stream Restorations

The goods and services provided by riverine systems are critical to humanity, and our reliance increases with our growing population and demands. As our activities expand, these systems continue to degrade throughout the world even as we try to restore them, and many efforts have not met expectations. One way to increase restoration effectiveness could be to explicitly design restorations to promote microbial communities, which are responsible for much of the organic matter breakdown, nutrient removal or transformation, pollutant removal, and biomass production in river ecosystems. In this paper, we discuss several design concepts that purposefully create conditions for these various microbial goods and services, and allow microbes to act as ecological restoration engineers. Focusing on microbial diversity and function could improve restoration effectiveness and overall ecosystem resilience to the stressors that caused the need for the restoration. Advances in next-generation sequencing now allow the use of microbial 'omics techniques (e.g., metagenomics, metatranscriptomics) to assess stream ecological conditions in similar fashion to fish and benthic macroinvertebrates. Using representative microbial communities from stream sediments, biofilms, and the water column may greatly advance assessment capabilities. Microbes can assess restorations and ecosystem function where animals may not currently be present, and thus may serve as diagnostics for the suitability of animal reintroductions. Emerging applications such as ecological metatranscriptomics may further advance our understanding of the roles of specific restoration designs towards ecological services as well as assess restoration effectiveness.

50%, not great, not terrible: Pan-European gap-analysis shows the real status of the DNA barcode reference libraries in two aquatic invertebrate groups and points the way ahead

The essential key to routine molecular species identification (DNA barcoding/metabarcoding) is the existence of an error-free DNA barcode reference library providing full coverage of all species. Published studies generally state the need to produce more barcodes, and control their quality, but unfortunately, the number of barcoded species is still low. However, to initiate real progress, we need to know where the gaps lie, how big they are and why they persist. Our aims were to draw and understand the current state of knowledge regarding species diversity, distribution, and barcode coverage, and offer solutions for improvement. In this study, we used two groups of aquatic insects, beetles and true bugs. We have compiled and critically evaluated an essentially complete and up-to-date European list, containing 1527 species. The list served as a basis for the barcode gap analyses in the Barcode-of-Life-Data-System (BOLD) conducted in three subsequent years (2020-2022). The overall barcode coverage of the pan-European fauna was around 50 % in both groups. The lowest coverage was in the Mediterranean, the Balkans and South-eastern Europe. The coverage in each country depended significantly on the local diversity, the number of rare, endemic species and the similarity of its fauna to that of the most active barcoding European countries. Gap analyses showed a very small increase in species coverage (<1 % in European aquatic beetles) despite an ~25 % increase in the number of barcodes. Hence, it is clear that future barcoding campaigns must prioritise quality over quantity. To visibly improve reference libraries, we need to increase the involvement of taxonomic experts and focus on targeted studies and underexplored but biodiversity-rich areas.

Multi-group biodiversity distributions and drivers of metacommunity organization along a glacial-fluvial-limnic pathway on the Tibetan plateau

Extensive global glacial retreats are threatening cryosphere ecosystem functioning and the associated biota in glacier-fed water systems. Understanding multi-group biodiversity distributions and compositional variation across diverse but hydrologically linked habitats under varying glacial influences will help explain the mechanisms underlying glacial community organization and ecosystem processes. However, such data are generally lacking due to the difficulty of obtaining biodiversity information across wide taxonomic ranges. Here, we used a multi-marker environmental DNA metabarcoding approach to simultaneously investigate the spatial patterns of community compositions and assembly mechanisms of four taxonomic groups (cyanobacteria, diatoms, invertebrates, and vertebrates) along the flowpaths of a tributary of Lake Nam Co on the Tibetan Plateau-from its glacier headwaters, through its downstream river and wetlands, to its estuary. We detected 869 operational taxonomic units: 119 cyanobacterial, 395 diatom, 269 invertebrate, and 86 vertebrate. Taxonomic richnesses consistently increased from upstream to downstream, and although all groups showed community similarity distance decay patterns, the trend for vertebrates was the weakest. Cyanobacteria, diatom, and invertebrate community compositions were significantly correlated with several environmental factors, while the vertebrate community was only correlated with waterway width. Variation partitioning analysis indicated that varying extents of environmental conditions and spatial factors affected community organizations for different groups. Furthermore, stochastic processes contributed prominently to the microorganisms' community assembly (Sloan's neutral model R² = 0.77 for cyanobacteria and 0.73 for diatoms) but were less important for macroorganisms (R² = 0.21 for invertebrates and 0.15 for vertebrates). That trend was further substantiated by modified stochasticity ratio analyses. This study provides the first holistic picture of the diverse biotic communities residing in a series of hydrologically connected glacier-influenced habitats. Our results both uncovered the distinct mechanisms that underlie the metacommunity organizations of different glacial organisms and helped comprehensively predict the ecological impacts of the world's melting glaciers.

Delimitation of Some Taxa of Ulnaria and Fragilaria (Bacillariophyceae) Based on Genetic, Morphological Data and Mating Compatibility

Fragilaria and Ulnaria are two closely related diatom genera for which the delimitation and circumscription of several species is unclear. We studied strains isolated from Lake Baikal and compared them with the species from freshwater reservoirs in Europe and Asia using phylogenetic and species delimitation methods, microscopy and interclonal crossing experiments. The results of the phylogenetic analyses of the fragments of rbcL and 18S rRNA genes revealed that baikalian F. radians clade was independent from the representatives of the genus from other localities. Among Ulnaria we found the following 18S rRNA phylogenetic tree groups at species level: U. acus, U. ulna and U. danica. Genetic distance between genera varied between 3.9–10.2% substitutions in rbcL gene and 3.2–11.5% in 18S rRNA. The boundary between intraspecies and interspecies polymorphism for studied species of Ulnaria and Fragilaria in these marker genes was around 0.8% substitutions. Morphometric characters of individual strains showed their variability and division into F. radians, U. acus and U. ulna together with U. danica. Strains of U. acus and U. danica from different localities of Europe and Asia were sexually compatible inside the species. Sexual reproduction has never been observed in monoclonal cultures, either between this species or with strains of the Fragilaria.

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 (&lt; 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.

More than just hitchhikers: a survey of bacterial communities associated with diatoms originating from sea turtles

Diatoms and bacteria are known for being the first colonizers of submerged surfaces including the skin of marine reptiles. Sea turtle carapace and skin harbor diverse prokaryotic and eukaryotic microbes, including several epizoic diatoms. However, the importance of diatom-bacteria associations is hardly investigated in biofilms associated with animal hosts. This study provides an inventory of diatoms, bacteria and diatom-associated bacteria originating from loggerhead sea turtles using both metabarcoding and culturing approaches. Amplicon sequencing of the carapace and skin samples chloroplast gene rbcL and 16S rRNA gene detected, in total, 634 diatom amplicon sequence variants (ASVs) and 3661 bacterial ASVs, indicating high diversity. Cultures of putative epizoic and non-epizoic diatoms contained 458 bacterial ASVs and their bacterial assemblages reflected those of their host. Diatom strains allowed for enrichment and isolation of bacterial families rarely observed on turtles, such as Marinobacteraceae, Alteromonadaceae and Alcanivoracaceae. When accounting for phylogenetic relationships between bacterial ASVs, we observed that related diatom genera might retain similar microbial taxa in culture, regardless of the turtle's skin or carapace source. These data provide deeper insights into the sea turtle-associated microbial communities, and reveal the potential of epizoic biofilms as a source of novel microbes and possibly important diatom-bacteria associations.

Metabarcoding as an effective complement of microscopic studies in revealing the composition of the diatom community – a case study of an oxbow lake of Tisza River (Hungary) with the description of a new Mayamaea species

Diatoms are valuable bioindicators and their traditional classification and identification are mainly based on the morphological characteristics of their frustules. However, in recent years, DNA-based methods have been proposed and are rapidly growing in the scientific literature as a complementary tool to assess the ecological status of freshwaters. Diatom-based ecological status assessment uses indices calculated from sensitivity and tolerance values as well as relative abundance of species. Correct assessment requires an accurate identification of species. In the present study, diatom assemblages of an oxbow lake were investigated using light and scanning electron microscopy as well as metabarcoding using rbcL marker, and the identification results were compared, intending to match barcode sequences of species that are currently missing in the diatom reference database. The investigated oxbow is an important wetland for bird conservation, although it is impacted by land use. Taxon lists based on morphology and metabarcoding considerably differed when bioinformatics analysis involved DADA2 pipeline with Diat.barcode database. Previously unknown sequence variants of four pennate species were found with additional BLAST search. Using phylogeny and p-distance calculations sequences could be matched to three small-celled naviculoid species that were found under a microscope. One of them was found to be a new species of the genus Mayamaea and was described as a new species, Mayamaea ectorii. Additionally, spatial distribution maps for several small-celled naviculoid species are provided for the Hungarian territory.

MetaWorks: A flexible, scalable bioinformatic pipeline for high-throughput multi-marker biodiversity assessments

Multi-marker metabarcoding is increasingly being used to generate biodiversity information across different domains of life from microbes to fungi to animals such as for molecular ecology and biomonitoring applications in different sectors from academic research to regulatory agencies and industry. Current popular bioinformatic pipelines support microbial and fungal marker analysis, while ad hoc methods are often used to process animal metabarcode markers from the same study. MetaWorks provides a harmonized processing environment, pipeline, and taxonomic assignment approach for demultiplexed Illumina reads for all biota using a wide range of metabarcoding markers such as 16S, ITS, and COI. A Conda environment is provided to quickly gather most of the programs and dependencies for the pipeline. Several workflows are provided such as: taxonomically assigning exact sequence variants, provides an option to generate operational taxonomic units, and facilitates single-read processing. Pipelines are automated using Snakemake to minimize user intervention and facilitate scalability. All pipelines use the RDP classifier to provide taxonomic assignments with confidence measures. We extend the functionality of the RDP classifier for taxonomically assigning 16S (bacteria), ITS (fungi), and 28S (fungi), to also support COI (eukaryotes), rbcL (eukaryotes, land plants, diatoms), 12S (fish, vertebrates), 18S (eukaryotes, diatoms) and ITS (fungi, plants). MetaWorks properly handles ITS by trimming flanking conserved rRNA gene regions as well as protein coding genes by providing two options for removing obvious pseudogenes. MetaWorks can be downloaded from https://github.com/terrimporter/MetaWorks and quickstart instructions, pipeline details, and a tutorial for new users can be found at https://terrimporter.github.io/MetaWorksSite.

Molecular and Morphological Investigations of Two Giant Diatom Cymbella Species from the Transbaikal Area (Russia, Siberia) with Comments on Their Distributions

For the first time, a giant diatom species of the genus Cymbella from Lake Baikal was studied using molecular methods. Molecular and morphological investigations allowed to us to described one new species, Cymbella baicalaspera Glushchenko, Kulikovskiy and Kociolek sp. nov. This species is both morphologically similar and phylogenetically close to a second giant Cymbella species that we investigated here, identified by us as Cymbella himalaspera Jüttner and Van de Vijver in Jüttner et al. 2010. This species was first described from Nepal on the basis of a morphological investigation. Small morphological differences exist between the type population and specimens from Lake Baikal, but otherwise the two are identical. These very interesting results show that some Baikalian diatoms can be distributed more widely and are not only endemic to this ancient lake. Similarity between Cymbella baicalaspera sp. nov. and Cymbella himalaspera on the basis of both morphological features and their close phylogenetic relationships suggested by molecular data indicate they are sister species and an example of sympatric speciation. These results also suggest an early development of a species flock. This species group warrants additional research in terms of. their diversification and biogeography.

Diatoms as indicators in running waters: trends of studies on biological assessment and monitoring

Despite the biodiversity and ecosystem services provided by lotic ecosystems, they are strongly affected by anthropogenic activities. Therefore, biological monitoring and assessment strategies are crucial in helping maintain these ecosystems and developing mitigation policies. We provide a global overview of the use of benthic diatoms as bioindicators in lotic environments, by analyzing 764 articles published in the past 20 years. We analyzed the influence of substrate type on samplings, which species have been highlighted as indicators and for which type of impacts, which anthropogenic impacts have been most commonly evaluated, and which metrics have been commonly used in studies using diatoms to assess and monitor the quality of lotic environments. We found that the most studied anthropogenic impact is artificial eutrophication and that some species, especially Nitzschia palea, have been thoroughly mentioned as indicators of this impact. Indicator species related to other types of impact are less common, demonstrating the need for studies on this issue. Moreover, we verified that traditional taxonomic metrics, such as diversity and diatom indices, have been widely used. Some alternative metrics have been used recently, such as those based on teratological valves, lipid bodies, valve size, and DNA metabarcoding. The number of biomonitoring and assessment studies based on diatoms has increased considerably in the past 20 years. Nonetheless, the demand for natural resources and consequently the degradation of lotic ecosystems have accelerated significantly. Thus, the development of low-cost and time-efficient biological assessment and monitoring strategies is essential for evaluating the health of lotic environments.

DNA sequence and taxonomic gap analyses to quantify the coverage of aquatic cyanobacteria and eukaryotic microalgae in reference databases: Results of a survey in the Alpine region

The taxonomic identification of organisms based on the amplification of specific genetic markers (metabarcoding) implicitly requires adequate discriminatory information and taxonomic coverage of environmental DNA sequences in taxonomic databases. These requirements were quantitatively examined by comparing the determination of cyanobacteria and microalgae obtained by metabarcoding and light microscopy. We used planktic and biofilm samples collected in 37 lakes and 22 rivers across the Alpine region. We focused on two of the most used and best represented genetic markers in the reference databases, namely the 16S rRNA and 18S rRNA genes. A sequence gap analysis using blastn showed that, in the identity range of 99-100%, approximately 30% (plankton) and 60% (biofilm) of the sequences did not find any close counterpart in the reference databases (NCBI GenBank). Similarly, a taxonomic gap analysis showed that approximately 50% of the cyanobacterial and eukaryotic microalgal species identified by light microscopy were not represented in the reference databases. In both cases, the magnitude of the gaps differed between the major taxonomic groups. Even considering the species determined under the microscope and represented in the reference databases, 22% and 26% were still not included in the results obtained by the blastn at percentage levels of identity ≥95% and ≥97%, respectively. The main causes were the absence of matching sequences due to amplification and/or sequencing failure and potential misidentification in the microscopy step. Our results quantitatively demonstrated that in metabarcoding the main obstacles in the classification of 16S rRNA and 18S rRNA sequences and interpretation of high-throughput sequencing biomonitoring data were due to the existence of important gaps in the taxonomic completeness of the reference databases and the short length of reads. The study focused on the Alpine region, but the extent of the gaps could be much greater in other less investigated geographic areas.

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.

Two new species of Navicula (Bacillariophyta) from Southeast Asia

We present light and scanning electron microscopical observations on two new species of Navicula Bory sensu stricto from Southeast Asia. Naviculawinoniformis Chudaev, Glushchenko, Kulikovskiy & Kociolek, sp. nov. differs from taxa with similar outline and size by the combination of simple drop-like external raphe endings deflected to the primary valve side, presence of well-developed external longitudinal grooves and relatively low lineolae density. Naviculasparsilineolata Chudaev, Glushchenko, Kulikovskiy & Kociolek, sp. nov. can be discriminated from the taxa of similar valve shape and size by consistently lower lineolae density, and from the majority of them also by the structure of external proximal raphe endings with small projections in proximal parts and larger triangular insertions in distal parts. Some remarks on Navicula species diversity and its distribution in the Southeast Asia are given.

Assessment of marine benthic diatom communities: insights from a combined morphological-metabarcoding approach in Mediterranean shallow coastal waters

We investigated the advantages and disadvantages of light microscope (LM)-based identifications and DNA metabarcoding, based on a 312-bp rbcL marker, for examining benthic diatom communities from Mediterranean shallow coastal environments. For this, we used biofilm samples collected from different substrata in the Ebro delta bays. We show that 1) Ebro delta bays harbour high-diversity diatom communities [LM identified 249 taxa] and 2) DNA metabarcoding effectively reflects this diversity at genus- but not species level, because of the incompleteness of the DNA reference library. Nevertheless, DNA metabarcoding offers new opportunities for detecting small, delicate and rare diatom species missed by LM and diatoms that lack silica frustules. The primers used, though designed for diatoms, successfully amplified rarely reported members of other stramenopile groups. Combining LM and DNA approaches offers stronger support for ecological studies of benthic microalgal communities in shallow coastal environments than using either approach on its own.

Phytool, a ShinyApp to homogenise taxonomy of freshwater microalgae from DNA barcodes and microscopic observations

Methods for biomonitoring of freshwater phytoplankton are evolving rapidly with eDNA-based methods, offering great complementarity with microscopy. Metabarcoding approaches have been more commonly used over the last years, with a continuous increase in the amount of data generated. Depending on the researchers and the way they assigned barcodes to species (bioinformatic pipelines and molecular reference databases), the taxonomic assignment obtained for HTS DNA reads might vary. This is also true for traditional taxonomic studies by microscopy with regular adjustments of the classification and taxonomy.

For those reasons (leading to non-homogeneous taxonomies), gap-analyses and comparisons between studies become even more challenging and the curation processes to find potential consensus names are time-consuming. Here, we present a web-based application (Phytool), developed with ShinyApp (Rstudio), that aims to make the harmonisation of taxonomy easier and in a more efficient way, using a complete and up-to-date taxonomy reference database for freshwater microalgae. Phytool allows users to homogenise and update freshwater phytoplankton taxonomical names from sequence files and data tables directly uploaded in the application. It also gathers barcodes from curated references in a user-friendly way in which it is possible to search for specific organisms. All the data provided are downloadable with the possibility to apply filters in order to select only the required taxa and fields (e.g. specific taxonomic ranks). The main goal is to make accessible to a broad range of users the connection between microscopy and molecular biology and taxonomy through different ready-to-use functions. This study estimates that only 25% of species of freshwater phytoplankton in Phytobs are associated with a barcode. We plead for an increased effort to enrich reference databases by coupling taxonomy and molecular methods. Phytool should make this crucial work more efficient.

The application is available at https://caninuzzo.shinyapps.io/phytool_v1/

Fish eDNA metabarcoding from aquatic biofilm samples: Methodological aspects

Fish eDNA metabarcoding is usually performed from filtered water samples. The volume of filtered water depends on the study scope and can rapidly become time consuming according to the number of samples that have to be processed. To avoid time allocated to filtration, passive DNA samplers have been used to recover fish eDNA from marine environments faster. In freshwater ecosystems, aquatic biofilms were used to catch eDNA from macroinvertebrates. Here, we test the capacity of aquatic biofilms to entrap fish eDNA in a large lake and, therefore, the possibility to perform fish eDNA metabarcoding from this matrix compared to the traditional fish eDNA approach from filtered water samples. Methodological aspects of the use of aquatic biofilms for fish eDNA metabarcoding (e.g. PCR replicates, biological replicates, bioinformatics pipeline, reference database and taxonomic assignment) were validated against a mock community. When using biofilms from habitats sheltered from wind and waves, biofilm and water approach provided similar inventories. Richness and diversity were comparable between both approaches. Approaches differed only for rare taxa. Our results illustrate the capacity of aquatic biofilms to act as passive eDNA samplers of fish eDNA and, therefore, the possibility to use biofilms to monitor fish communities efficiently from biofilms. Furthermore, our results open up avenues of research to study a diversity of biological groups (among which bioindicators as diatoms, macroinvertebrates and fish) from eDNA isolated from a single environmental matrix reducing sampling efforts, analysis time and costs.

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.

Herbicide-Induced Shifts in the Periphyton Community Composition Indirectly Affect Feeding Activity and Physiology of the Gastropod Grazer Physella acuta

Herbicides are well known for unintended effects on freshwater periphyton communities. Large knowledge gaps, however, exist regarding indirect herbicide impacts on primary consumers through changes in the quality of periphyton as a food source (i.e., diet-related effects). To address this gap, the grazer Physella acuta (Gastropoda) was fed for 21 days with periphyton that grew for 15 days in the presence or absence of the herbicide diuron (8 μg/L) to quantify changes in the feeding rate, growth rate, and energy storage (neutral lipid fatty acids; NLFAs) of P. acuta. Periphyton biomass, cell viability, community structure, and FAs served as proxies for food quality that support a mechanistic interpretation of the grazers' responses. Diuron changed the algae periphyton community and fatty acid profiles, indicating alterations in the food quality, which could explain differences in the snails' feeding rate compared to the control. While the snails' growth rate was, despite an effect size of 55%, not statistically significantly changed, NLFA profiles of P. acuta were altered. These results indicate that herbicides can change the food quality of periphyton by shifts in the algae composition, which may affect the physiology of grazers.

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.

Dumpster diving for diatom plastid 16S rRNA genes

High throughput sequencing is improving the efficiency of monitoring diatoms, which inhabit and support aquatic ecosystems across the globe. In this study, we explored the potential of a standard V4 515F-806RB primer pair in recovering diatom plastid 16S rRNA sequences. We used PhytoREF to classify the 16S reads from our freshwater biofilm field sampling from three stream segments across two streams in south-eastern Australia and retrieved diatom community data from other, publicly deposited, Australian 16S amplicon datasets. When these diatom operational taxonomic units (OTUs) were traced using the default RDPII and NCBI databases, 68% were characterized as uncultured cyanobacteria. We analysed the 16S rRNA sequences from 72 stream biofilm samples, separated the chloroplast OTUs, and classified them using the PhytoREF database. After filtering the reads attributed to Bacillariophyta (relative abundance >1%), 71 diatom OTUs comprising more than 90% of the diatom reads in each stream biofilm sample were identified. Beta-diversity analyses demonstrated significantly different diatom assemblages and discrimination among river segments. To further test the approach, the diatom OTUs from our biofilm sampling were used as reference sequences to identify diatom reads from other Australian 16S rRNA datasets in the NCBI-SRA database. Across the three selected public datasets, 67 of our 71 diatom OTUs were detected in other Australian ecosystems. Our results show that diatom plastid 16S rRNA genes are readily amplified with existing 515F-806RB primer sets. Therefore, the volume of existing 16S rRNA amplicon datasets initially generated for microbial community profiling can also be used to detect, characterize, and map diatom distribution to inform phylogeny and ecological health assessments, and can be extended into a range of ecological and industrial applications. To our knowledge, this study represents the first attempt to classify freshwater samples using this approach and the first application of PhytoREF in Australia.

Cyprus' diatom diversity and the association of environmental and anthropogenic influences for ecological assessment of rivers using DNA metabarcoding

Human activities are the leading cause of environmental impairments. Appropriate biomonitoring of ecosystems is needed to assess these activities effectively. In freshwater ecosystems, periphytic and epilithic biofilms have diatom assemblages. These assemblages respond rapidly to environmental changes, making diatoms valuable bioindicators. For this reason, freshwater biomonitoring programs are currently using diatoms (e.g., Water Framework Directive). In the past ten years, DNA metabarcoding coupled with next-generation sequencing and bioinformatics represents a complementary approach for diatom biomonitoring. In this study, this approach is used for the first time in Cyprus by considering the association of environmental and anthropogenic pressures to diatom assemblages. Statistical analysis was then applied to identify the environmental (i.e., river types, geo-morphological) and anthropogenic (i.e., physicochemical, human land-use pressures) variables' role in the observed diatom diversity. Results indicate differences in diatom assemblages between intermittent and perennial rivers. Achnanthidium minutissimum was more abundant in intermittent rivers; whereas Amphora pediculus and Planothidium caputium in perennial ones. Additionally, we could demonstrate the correlation between nutrients (e.g., nitrogen, phosphorus), stations' local characteristics (e.g., elevation), and land use activities on the observed differences in diatom diversity. Finally, we conclude that multi-stressors and anthropogenic pressures together as multiple stressors have a significant statistical relationship to the observed diatom diversity and play a pivotal role in determining Cyprus' rivers' ecological status.

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.

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.

Metadata standards and practical guidelines for specimen and DNA curation when building barcode reference libraries for aquatic life

DNA barcoding and metabarcoding is increasingly used to effectively and precisely assess and monitor biodiversity in aquatic ecosystems. As these methods rely on data availability and quality of barcode reference libraries, it is important to develop and follow best practices to ensure optimal quality and traceability of the metadata associated with the reference barcodes used for identification. Sufficient metadata, as well as vouchers, corresponding to each reference barcode must be available to ensure reliable barcode library curation and, thereby, provide trustworthy baselines for downstream molecular species identification. This document (1) specifies the data and metadata required to ensure the relevance, the accessibility and traceability of DNA barcodes and (2) specifies the recommendations for DNA harvesting and for the storage of both voucher specimens/samples and barcode data.

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

Ecological assessment of lakes and rivers using benthic diatom assemblages currently requires considerable taxonomic expertise to identify species using light microscopy. This traditional approach is also time-consuming. Diatom metabarcoding is a promising alternative and there is increasing interest in using this approach for routine assessment. However, until now, analysis protocols for diatom metabarcoding have been developed and optimised by research groups working in isolation. The diversity of existing bioinformatics methods highlights the need for an assessment of the performance and comparability of results of different methods. The aim of this study was to test the correspondence of outputs from six bioinformatics pipelines currently in use for diatom metabarcoding in different European countries. Raw sequence data from 29 biofilm samples were treated by each of the bioinformatics pipelines, five of them using the same curated reference database. The outputs of the pipelines were compared in terms of sequence unit assemblages, taxonomic assignment, biotic index score and ecological assessment outcomes. The three last components were also compared to outputs from traditional light microscopy, which is currently accepted for ecological assessment of phytobenthos, as required by the Water Framework Directive. We also tested the performance of the pipelines on the two DNA markers (rbcL and 18S-V4) that are currently used by the working groups participating in this study. The sequence unit assemblages produced by different pipelines showed significant differences in terms of assigned and unassigned read numbers and sequence unit numbers. When comparing the taxonomic assignments at genus and species level, correspondence of the taxonomic assemblages between pipelines was weak. Most discrepancies were linked to differential detection or quantification of taxa, despite the use of the same reference database. Subsequent calculation of biotic index scores also showed significant differences between approaches, which were reflected in the final ecological assessment. Use of the rbcL marker always resulted in better correlation among molecular datasets and also in results closer to these generated using traditional microscopy. This study shows that decisions made in pipeline design have implications for the dataset's structure and the taxonomic assemblage, which in turn may affect biotic index calculation and ecological assessment. There is a need to define best-practice bioinformatics parameters in order to ensure the best representation of diatom assemblages. Only the use of similar parameters will ensure the compatibility of data from different working groups. The future of diatom metabarcoding for ecological assessment may also lie in the development of new metrics using, for example, presence/absence instead of relative abundance data.