DNA metabarcoding effectively quantifies diatom responses to nutrients in streams

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.

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.

The impact of flooding from the Minjiang River on the succession of harmful algal blooms (HABs) caused by diatoms in China's offshore waters

This study examines diatom assemblages in the Matsu Archipelago, an area influenced by Minjiang River runoff. It focuses on harmful algal blooms (HABs) that occurred between August 2021 and July 2022. Utilizing 18S rRNA metabarcoding and microscopic analysis, we observed a significant diatom bloom during early summer runoff, peaking at 5 × 105 cells L-1. The research reveals dynamic community changes during the runoff season, with dominant genera including Pseudo-nitzschia, Chaetoceros, and Skeletonema. Skeletonema cell density correlated with NO3 levels, Chaetoceros had a slight PO4 affinity, and Pseudo-nitzschia showed a negative correlation with Skeletonema. Pseudo-nitzschia, which prefers high light and pH conditions, had notably high concentrations in the flood season and in the autumn. In both, it was dominated by potential toxin-producing species - P. multistriata and P. pungens during the flooding, and P. cuspidate in the autumn. These findings highlight the intricate relationship between diatom dynamics and environmental factors, providing essential insights for managing HABs, especially Pseudo-nitzschia species, amidst environmental changes.

Applying water quality standards to pollution from diffuse sources

Water quality standards are instrumental in evaluating the status of water bodies, and in providing protective and restorative endpoints. To date, much of the infrastructure used to implement water quality standards has been directed towards remediating and managing pollution from point source discharges. However, pollution from diffuse sources is the leading cause of water quality impairment, especially by nutrients. Although the effects of nutrient enrichment on streams is well studied, and ecological thresholds identified, those thresholds have not been widely adopted as standards primarily because they are not attainable by point sources. Clearly, a framework for adopting and applying standards to manage pollution from diffuse sources needs to be decoupled from those intended for point sources. This paper argues for a relatively unstructured distributional approach to predict how ecological responses might shift in response to management of diffuse sources. The approach calls for first developing a deterministic model of stressor and response variables, followed by a reformulation as a Bayesian model. In the case here, a structural equation model was developed that linked nutrient enrichment, habitat quality, and chloride and manganese concentrations to an index of macroinvertebrate quality. Results from the Bayesian representation suggest that in landscapes where the drainage network has been highly modified for agricultural production, reduction in total phosphorus alone is expected to have a modest (but non-trivial) effect on macroinvertebrate condition, shifting the distribution of scores up by 1 point. The addition of habitat restoration is likely to shift the distribution upwards by 4 points, an effect size observed in Ohio, USA from other large-scale restoration efforts.

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.

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.

Understanding the effects of phosphorus on diatom richness in rivers and streams using taxon-environment relationships

Changes in phosphorus concentrations affect periphytic diatom composition in streams, yet we rarely observe strong relationships between diatom richness and phosphorus. In contrast, changes in conductivity are strongly associated with differences in both diatom composition and richness. We hypothesised that we could better understand the mechanisms that control the phosphorus-richness relationship by examining relationships between phosphorus and the occurrence of individual diatom taxa, comparing these with relationships between conductivity and taxon occurrence, and documenting how niche breadths of taxa affect richness patterns. We estimated relationships between phosphorus and taxon occurrence using DNA metabarcoding data of diatoms collected from 1,811 sites distributed across the conterminous U.S.A. and contrasted patterns in these relationships with those between conductivity and taxon occurrence. The distribution of taxon optima for phosphorus was bimodal, with most optima located at either the maximum or minimum observed phosphorus concentration. The distribution of taxon optima for conductivity was unimodal. Niche breadths of taxa for phosphorus and for conductivity both generally increased with optimum values. The distribution of conductivity optima gave rise to a prominent hump-shaped relationship between richness and conductivity. The relationship between richness and phosphorus was also slightly hump-shaped, but this relationship would not be expected from the bimodal distribution of optima. Instead, we determined that broad niche breadths caused the hump-shaped relationship between richness and phosphorus. Our results highlight the nuanced effects that increased P loadings exert on diatom assemblages in rivers and streams and identify reasons that weak relationships between taxon richness and increased phosphorus have been observed. These findings allow us to better describe how excess phosphorus and subsets of taxa and their niche breadths contribute to patterns of taxa richness in diatom assemblages, and to improve the tools used to manage phosphorus pollution.

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.

Developing Indicators of Nutrient Pollution in Streams Using 16S rRNA Gene Metabarcoding of Periphyton-Associated Bacteria

Indicators based on nutrient-biota relationships in streams can inform water quality restoration and protection programs. Bacterial assemblages could be particularly useful indicators of nutrient effects because they are species-rich, important contributors to ecosystem processes in streams, and responsive to rapidly changing conditions. Here, we sampled 25 streams weekly (12-14 times each) and used 16S rRNA gene metabarcoding of periphyton-associated bacteria to quantify the effects of total phosphorus (TP) and total nitrogen (TN). Threshold indicator taxa analysis identified assemblage-level changes and amplicon sequence variants (ASVs) that increased or decreased with increasing TP and TN concentrations (i.e., low P, high P, low N, and high N ASVs). Boosted regression trees confirmed that relative abundances of gene sequence reads for these four indicator groups were associated with nutrient concentrations. Gradient forest analysis complemented these results by using multiple predictors and random forest models for each ASV to identify portions of TP and TN gradients at which the greatest changes in assemblage structure occurred. Synthesized statistical results showed bacterial assemblage structure began changing at 24 μg TP/L with the greatest changes occurring from 110 to 195 μg/L. Changes in the bacterial assemblages associated with TN gradually occurred from 275 to 855 μg/L. Taxonomic and phylogenetic analyses showed that low nutrient ASVs were commonly Firmicutes, Verrucomicrobiota, Flavobacteriales, and Caulobacterales, Pseudomonadales, and Rhodobacterales of Proteobacteria, whereas other groups, such as Chitinophagales of Bacteroidota, and Burkholderiales, Rhizobiales, Sphingomonadales, and Steroidobacterales of Proteobacteria comprised the high nutrient ASVs. Overall, the responses of bacterial ASV indicators in this study highlight the utility of metabarcoding periphyton-associated bacteria for quantifying biotic responses to nutrient inputs in streams.

Characterizing temporal variability in streams supports nutrient indicator development using diatom and bacterial DNA metabarcoding

Interest in developing periphytic diatom and bacterial indicators of nutrient effects continues to grow in support of the assessment and management of stream ecosystems and their watersheds. However, temporal variability could confound relationships between indicators and nutrients, subsequently affecting assessment outcomes. To document how temporal variability affects measures of diatom and bacterial assemblages obtained from DNA metabarcoding, we conducted weekly periphyton and nutrient sampling from July to October 2016 in 25 streams in a 1293 km2 mixed land use watershed. Measures of both diatom and bacterial assemblages were strongly associated with the percent agriculture in upstream watersheds and total phosphorus (TP) and total nitrogen (TN) concentrations. Temporal variability in TP and TN concentrations increased with greater amounts of agriculture in watersheds, but overall diatom and bacterial assemblage variability within sites-measured as mean distance among samples to corresponding site centroids in ordination space-remained consistent. This consistency was due in part to offsets between decreasing variability in relative abundances of taxa typical of low nutrient conditions and increasing variability in those typical of high nutrient conditions as mean concentrations of TP and TN increased within sites. Weekly low and high nutrient diatom and bacterial metrics were more strongly correlated with site mean nutrient concentrations over the sampling period than with same day measurements and more strongly correlated with TP than with TN. Correlations with TP concentrations were consistently strong throughout the study except briefly following two major precipitation events. Following these events, biotic relationships with TP reestablished within one to three weeks. Collectively, these results can strengthen interpretations of survey results and inform monitoring strategies and decision making. These findings have direct applications for improving the use of diatoms and bacteria, and the use of DNA metabarcoding, in monitoring programs and stream site assessments.

Quantifying spatial and temporal relationships between diatoms and nutrients in streams strengthens evidence of nutrient effects from monitoring data

Observational data are frequently used to better understand the effects of changes in P and N on stream biota, but nutrient gradients in streams are usually associated with gradients in other environmental factors, a phenomenon that complicates efforts to accurately estimate the effects of nutrients. Here, we propose a new approach for analyzing observational data in which we compare the effects of changes in nutrient concentrations in time within individual sites and in space among many sites. Covarying relationships between other, potentially confounding environmental factors and nutrient concentrations are unlikely to be the same in both time and space, and, therefore, estimated effects of nutrients that are similar in time and space are more likely to be accurate. We applied this approach to diatom rbcL metabarcoding data collected from streams in the East Fork of the Little Miami River watershed, Ohio, USA. Changes in diatom assemblage composition were consistently associated with changes in the concentration of total reactive P in both time and space. In contrast, despite being associated with spatial differences in ammonia and urea concentrations, diatom assemblage composition was not associated with temporal changes in these nitrogen species. We suggest that the results of this analysis provide evidence of a causal effect of increased P on diatom assemblage composition. We further analyzed the effects of temporal variability in measurements of total reactive P and found that averaging periods greater than ~1 wk prior to sampling best represented the effects of P on the diatom assemblage. Comparisons of biological responses in space and time can sharpen insights beyond those that are based on analyses conducted on only 1 of the 2 dimensions.

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.

Finding needles in a haystack-Extensive diversity in the eustigmatophyceae revealed by community metabarcode analysis targeting the rbcL gene using lineage-directed primers

Sequences from the Stramenopile class Eustigmatophyceae are rarely reported in metabarcoding studies, and when they have been reported, there are very few haplotypes. We hypothesized that the paucity of eustigmatophyte species detected in these studies may be a result of the metabarcoding techniques used, which have primarily employed universal ribosomal RNA gene regions. In this study, we examined environmental DNA samples from 22 sites in southwestern Virginia, some of which had previously been studied using ribosomal RNA analysis. We used metabarcoding techniques targeting the plastid rbcL gene with new primers designed to produce a 370 bp amplicon from all lineages of the Eustigmatophyceae in a reference collection. The amplicons were then analyzed with DADA2 to produce amplicon sequence variants (ASVs). Our results revealed 184 rbcL haplotypes that can be tentatively assigned to the Eustigmatophyceae from these sites, representing much higher diversity than has been detected by ribosomal DNA-based studies. The techniques employed can be used for future studies of population structure, ecology, distribution, and diversity of this class. With these techniques, it should be possible to make realistic estimates of the species-level diversity of the Eustigmatophyceae on local, regional, and perhaps global scales.

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.