Spatial and temporal dynamics of a freshwater eukaryotic plankton community revealed via 18S rRNA gene metabarcoding

Metabarcoding reveals potentially mixotrophic flagellates and picophytoplankton as key groups of phytoplankton in the Elbe estuary

In estuaries, phytoplankton are faced with strong environmental forcing (e.g. high turbidity, salinity gradients). Taxa that appear under such conditions may play a critical role in maintaining food webs and biological carbon pumping, but knowledge about estuarine biota remains limited. This is also the case in the Elbe estuary where the lower 70 km of the water body are largely unexplored. In the present study, we investigated the phytoplankton composition in the Elbe estuary via metabarcoding. Our aim was to identify key taxa in the unmonitored reaches of this ecosystem and compare our results from the monitored area with available microscopy data. Phytoplankton communities followed distinct seasonal and spatial patterns. Community composition was similar across methods. Contributions of key classes and genera were correlated to each other (p < 0.05) when obtained from reads and biovolume (R2 = 0.59 and 0.33, respectively). Centric diatoms (e.g. Stephanodiscus) were the dominant group - comprising on average 55 % of the reads and 66-69 % of the biovolume. However, results from metabarcoding imply that microscopy underestimates the prevalence of picophytoplankton and flagellates with a potential for mixotrophy (e.g. cryptophytes). This might be due to their small size and sensitivity to fixation agents. We argue that mixotrophic flagellates are ecologically relevant in the mid to lower estuary, where, e.g., high turbidity render living conditions rather unfavorable, and skills such as phagotrophy provide fundamental advantages. Nevertheless, further findings - e.g. important taxa missing from the metabarcoding dataset - emphasize potential limitations of this method and quantitative biases can result from varying numbers of gene copies in different taxa. Further research should address these methodological issues but also shed light on the causal relationship of taxa with the environmental conditions, also with respect to active mixotrophic behavior.

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.

Diet Diversity of the Fluviatile Masu Salmon, Oncorhynchus masou (Brevoort 1856) Revealed via Gastrointestinal Environmental DNA Metabarcoding and Morphological Identification of Contents

Masu salmon, Oncorhynchus masou (Brevoort 1856), a commercially important fish species endemic to the North Pacific Ocean, attained national second-level protected animal status in China in 2021. Despite this recognition, knowledge about the trophic ecology of this fish remains limited. This study investigated the diet diversity of fluviatile Masu salmon in the Mijiang River, China, utilizing the gastrointestinal tract environmental DNA (GITeDNA) metabarcoding and morphological identification. The results revealed a diverse prey composition, ranging from terrestrial and aquatic invertebrates to small fishes. The fluviatile Masu salmon in general consumed noteworthily more aquatic prey than terrestrial prey. There were much more prey taxa and a higher diet diversity detected by GITeDNA metabarcoding than by morphological identification. GITeDNA metabarcoding showed that larger and older Masu salmon consumed significantly more terrestrial insects than aquatic prey species did, with 7366 verses 5012 sequences in the group of ≥20 cm, 9098 verses 4743 sequences in the group of ≥100 g and 11,540 verses 729 sequences in the group of age 3+. GITeDNA metabarcoding also showed size- and age-related diet diversity, indicating that the dietary niche breadth and trophic diversity of larger and older Masu salmon increased with food resources expanding to more terrestrial prey. Terrestrial invertebrates of riparian habitats play a vital role in the diet of fluviatile Masu salmon, especially larger individuals, highlighting their importance in connecting aquatic and terrestrial food webs. Conservation plans should prioritize the protection and restoration of riparian habitats. This study advocates the combined use of GITeDNA metabarcoding and morphological observation for a comprehensive understanding of fish diet diversity.

Efficient modulation of cellular phosphorus components in response to phosphorus deficiency in the dinoflagellate Karenia mikimotoi

IMPORTANCE

Dinoflagellates are the most common phytoplankton group and account for more than 75% of harmful algal blooms in coastal waters. In recent decades, dinoflagellates seem to prevail in phosphate-depleted waters. However, the underlying acclimation mechanisms and competitive strategies of dinoflagellates in response to phosphorus deficiency are poorly understood, especially in terms of intracellular phosphorus modulation and recycling. Here, we focused on the response of intracellular phosphorus metabolism to phosphorus deficiency in the model dinoflagellate Karenia mikimotoi. Our work reveals the strong capability of K. mikimotoi to efficiently regulate intracellular phosphorus resources, particularly through membrane phospholipid remodeling and miRNA regulation of energy metabolism. Our research improved the understanding of intracellular phosphorus metabolism in marine phytoplankton and underscored the advantageous strategies of dinoflagellates in the efficient modulation of internal phosphorus resources to maintain active physiological activity and growth under unsuitable phosphorus conditions, which help them outcompete other species in coastal phosphate-depleted environments.

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

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

Does filter pore size introduce bias in DNA sequence-based plankton community studies?

The cell size of microbial eukaryotic plankton normally ranges from 0.2 to 200 μm. During the past decade, high-throughput sequencing of DNA has been revolutionizing their study on an unprecedented scale. Nonetheless, it is currently unclear whether we can accurately, effectively, and quantitatively depict the microbial eukaryotic plankton community using size-fractionated filtration combined with environmental DNA (eDNA) molecular methods. Here we assessed the microbial eukaryotic plankton communities with two filtering strategies from two subtropical reservoirs, that is one-step filtration (0.2–200 μm) and size-fractionated filtration (0.2–3 and 3–200 μm). The difference of 18S rRNA gene copy abundance between the two filtering treatments was less than 50% of the 0.2–200 μm microbial eukaryotic community for 95% of the total samples. Although the microbial eukaryotic plankton communities within the 0.2–200 μm and the 0.2–3 and 3–200 μm size fractions had approximately identical 18S rRNA gene copies, there were significant differences in their community composition. Furthermore, our results demonstrate that the systemic bias introduced by size-fractionation filtration has more influence on unique OTUs than shared OTUs, and the significant differences in abundance between the two eukaryotic plankton communities largely occurred in low-abundance OTUs in specific seasons. This work provides new insights into the use of size-fractionation in molecular studies of microbial eukaryotes populating the plankton.

Monitoring harmful microalgal species and their appearance in Tokyo Bay, Japan, using metabarcoding

During the recent decade, high-throughput sequencing (HTS) techniques, in particular, DNA metabarcoding, have facilitated increased detection of biodiversity, including harmful algal bloom (HAB) species. In this study, the presence of HAB species and their appearance patterns were investigated by employing molecular and light microscopy-based monitoring in Tokyo Bay, Japan. The potential co-appearance patterns between the HAB species, as well as with other eukaryotes and prokaryotes were investigated using correlation and association rule-based time-series analysis. In total, 40 unique HAB species were detected, including 12 toxin-producing HAB species previously not reported from the area. More than half of the HAB species were present throughout the sampling season (summer to autumn) and no structuring or succession patterns associated with the environmental conditions could be detected. Statistically significant (p &lt; 0.05, rS ranging from −0.88 to 0.90) associations were found amongst the HAB species and other eukaryotic and prokaryotic species, including genera containing growth-limiting bacteria. However, significant correlations between species differed amongst the years, indicating that variability in environmental conditions between the years may have a stronger influence on the microalgal community structure and interspecies interactions than the variability during the sampling season. The association rule-based time-series analysis allowed the detection of a previously reported negative relationship between Synechococcus sp. and Skeletonema sp. in nature. Overall, the results support the applicability of metabarcoding and HTS-based microalgae monitoring, as it facilitates more precise species identification compared to light microscopy, as well as provides input for investigating potential interactions amongst different species/groups through simultaneous detection of multiple species/genera.

Airborne Prokaryotic, Fungal and Eukaryotic Communities of an Urban Environment in the UK

Bioaerosols often contain human pathogens and allergens affecting public health. However, relatively little attention has been given to bioaerosols compared with non-biological aerosols. In this study, we aimed to identify bioaerosol compositions in Manchester, UK by applying high throughput sequencing methods and to find potential sources. Samples were collected at Manchester Air Quality Super Site at the Firs Environmental Research Station in November 2019 and in February 2020. Total DNA has been extracted and sequenced targeting the 16S rRNA gene of prokaryotes, ITS region of fungal DNA and 18S rRNA gene of eukaryotes. We found marine environment-associated bacteria and archaea were relatively more abundant in the February 2020 samples compared with the November 2019 samples, consistent with the North West marine origin based on wind back-trajectory analysis. In contrast, an OTU belonging to Methylobacterium, which includes many species resistant to heavy metals, was relatively more abundant in November 2019 when there were higher metal concentrations. Fungal taxa that fruit all year were relatively more abundant in the February 2020 samples while autumn fruiting species generally had higher relative abundance in the November 2019 samples. There were higher relative abundances of land plants and algae in the February 2020 samples based on 18S rRNA gene sequencing. One of the OTUs belonging to the coniferous yew genus Taxus was more abundant in the February 2020 samples agreeing with the usual pollen season of yews in the UK which is from mid-January until late April. The result from this study suggests a potential application of bioaerosol profiling for tracing the source of atmospheric particles.

Multi-marker DNA metabarcoding detects suites of environmental gradients from an urban harbour

There is increasing need for biodiversity monitoring, especially in places where potential anthropogenic disturbance may significantly impact ecosystem health. We employed a combination of traditional morphological and bulk macroinvertebrate metabarcoding analyses to benthic samples collected from Toronto Harbour (Ontario, Canada) to compare taxonomic and functional diversity of macroinvertebrates and their responses to environmental gradients. At the species rank, sites assessed using COI metabarcoding showed more variation than sites assessed using morphological methods. Depending on the assessment method, we detected gradients in magnesium (morphological taxa), ammonia (morphological taxa, COI sequence variants), pH (18S sequence variants) as well as gradients in contaminants such as metals (COI & 18S sequence variants) and organochlorines (COI sequence variants). Observed responses to contaminants such as aromatic hydrocarbons and metals align with known patchy distributions in harbour sediments. We determined that the morphological approach may limit the detection of macroinvertebrate responses to lake environmental conditions due to the effort needed to obtain fine level taxonomic assignments necessary to investigate responses. DNA metabarcoding, however, need not be limited to macroinvertebrates, can be automated, and taxonomic assignments are associated with a certain level of accuracy from sequence variants to named taxonomic groups. The capacity to detect change using a scalable approach such as metabarcoding is critical for addressing challenges associated with biodiversity monitoring and ecological investigations.

High-Throughput Sequencing of Diatom Community, Its Spatial and Temporal Variation and Interrelationships with Physicochemical Factors in Danjiangkou Reservoir, China

Diatoms constitute an important part of the phytoplankton community in lakes and reservoirs and play a significant role in regulating ecological balance. Danjiangkou Reservoir is the water source area of the middle route of China’s South-to-North Water Diversion project. In order to explore the spatial and temporal distribution and know the governing factors of the diatom community, 18srRNA sequencing was carried out from seven sampling sites of the reservoir. At the same time, the concentration of nutrients present in the collected sample water was also determined. The results showed that a total of 51 genera and 96 species were thriving the community of diatoms in Danjiangkou Reservoir. Discostella was dominant in summer and autumn, accounting for 98.84% and 62.71% of the diatom abundance, respectively. Aulacoseira was dominant in spring and winter, accounting for 60.62% and 60.90%, respectively. Discostella and Aulacoseira showed significant differences in seasonal variation (p < 0.05). The colinear network of diatoms changed significantly with the seasons, mainly consisting of Aulacoseira, Discostella, and Stephanodiscus. RDA redundancy analysis showed that water temperature (WT), total nitrogen (TN), NH4+-N, pH, and electrical conductivity (Cond) were the main environmental factors driving the changes in diatom community structure.

Determination of the direct and indirect effects of bend on the urban river ecological heterogeneity

The ecological heterogeneity created by river bends benefits the diversity of microorganisms, which is vital for the pollutant degradation and overall river health. However, quantitative tools capable of determining the interactions among different trophic levels and species are lacking, and research regarding ecological heterogeneity has been limited to a few species. By integrating the multi-species-based index of biotic integrity (Mt-IBI) and the structure equation model (SEM), an interactions-based prediction modeling framework was established. Based on DNA metabarcoding, a multi-species (i.e., bacteria, protozoans, and metazoans) based index of biotic integrity including 309 candidate metrics was developed. After a three-step screening process, eight core metrics were obtained to assess the ecological heterogeneity, quantitatively. The Mt-IBI value, which ranged from 2.08 to 7.17, was calculated as the sum of each single core metric value. The Mt-IBI revealed that the ecological heterogeneity of concave banks was higher than other sites. According to the result of the SEM, D₉₀ was the controlling factor (r = -0.779) of the ecological heterogeneity under the influence of the river bends. The bend-induced redistribution of sediment particle further influenced the concentrations of carbon, nitrogen, and sulphur. The nitrogen group (r = 0.668) also played an essential role in determining the ecological heterogeneity, follow by carbon group (r = 0.455). Furthermore, the alteration of niches would make a difference on the ecological heterogeneity. This multi-species interactions-based prediction modeling framework proposed a novel method to quantify ecological heterogeneity and provided insight into the enhancement of ecological heterogeneity in river bends.

The Effects of Water Level Fluctuation on Zooplankton Communities in Shahu Lake Based on DNA Metabarcoding and Morphological Methods

Background: The water level of Poyang Lake (China) fluctuates seasonally. Shahu Lake, a smaller body of water connected to Poyang Lake during the wet season, is separated in the dry season. Due to a special fishing method termed ‘lake enclosed in autumn’, the water level is lowered and reaches its lowest point in January, which is <0.5 m deep in the middle of the lake. Our research investigated the effect of water level changes on the zooplankton community composition in Shahu Lake. Methods: We used both DNA metabarcoding method (MBC) (18S rRNA gene V4 region) and morphological method (MOI) to track the zooplankton community structure over four seasons in Shahu Lake (China). Results: Totals of 90 and 98 species of zooplankton were detected by MOI and MBC, respectively, with rotifers being the main zooplankton component. The α-diversity index of both methods increased from spring to summer and decreased from summer to autumn, reaching the lowest value in winter. NMDS and a cluster analysis showed that all zooplankton communities detected by MOI and MBC were significantly separated by season. The zooplankton community in winter was separated from that of the other three seasons, but the summer and autumn communities were more similar. Conclusions: Changes in the water level had significant effects on the zooplankton community composition. We found that MBC was more able to detect the differences in the zooplankton composition than MOI. MBC also had more advantages in copepod recognition. In our study, 37 species of copepods were detected by MBC, but only 11 species were detected by MOI. We concluded that MBC should be used to research the seasonal variations of zooplankton.

Seasonal Diversity of Microeukaryotes in the Han River, Korea Through 18S rRNA Gene Metabarcoding

Freshwater ecosystems contain a large diversity of microeukaryotes that play important roles in maintaining their structure. Microeukaryote communities vary in composition and abundance on the basis of temporal and environmental variables and may serve as useful bioindicators of environmental changes. In the present study, 18S rRNA metabarcoding was employed to investigate the seasonal diversity of microeukaryote communities during four seasons in the Han River, Korea. In total, 882 unique operational taxonomic units (OTUs) were detected, including various diatoms, metazoans (e.g., arthropods and rotifers), chlorophytes, and fungi. Although alpha diversity revealed insignificant differences based on seasons, beta diversity exhibited a prominent variation in the community composition as per seasons. The analysis revealed that the diversity of microeukaryotes was primarily driven by seasonal changes in the prevailing conditions of environmental water temperature and dissolved oxygen. Moreover, potential indicator OTUs belonging to diatoms and chlorophytes were associated with seasonal and environmental factors. This analysis was a preliminary study that established a continuous monitoring system using metabarcoding. This approach could be an effective tool to manage the Han River along with other freshwater ecosystems in Korea.

Plankton biodiversity and species co-occurrence based on environmental DNA – a multiple marker study

Metabarcoding in combination with high-throughput sequencing (HTS) allows simultaneous detection of multiple taxa by targeting single or several taxonomically informative gene regions from environmental DNA samples. In this study, a multiple-marker HTS approach was applied to investigate the plankton diversity and seasonal succession in the Baltic Sea from winter to autumn. Four different markers targeting the 16S, 18S, and 28S ribosomal RNA genes were employed, including a marker for more efficient dinoflagellate detection. Typical seasonal changes were observed in phyto- and bacterioplankton communities. In phytoplankton, the appearance patterns of selected common, dominant, or harmful species followed the patterns also confirmed based on 20 years of phytoplankton monitoring data. In the case of zooplankton, both macro- and microzooplankton species were detected. However, no seasonal patterns were detected in their appearance. In total, 15 and 2 new zoo- and phytoplankton species were detected from the Baltic Sea. HTS approach was especially useful for detecting microzooplankton species as well as for investigating the co-occurrence and potential interactions of different taxa. The results of this study further exemplify the efficiency of metabarcoding for biodiversity monitoring and the advantage of employing multiple markers through the detection of species not identifiable based on a single marker survey and/or by traditional morphology-based methods.

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.

Ecosystems monitoring powered by environmental genomics: A review of current strategies with an implementation roadmap

A decade after environmental scientists integrated high-throughput sequencing technologies in their toolbox, the genomics-based monitoring of anthropogenic impacts on the biodiversity and functioning of ecosystems is yet to be implemented by regulatory frameworks. Despite the broadly acknowledged potential of environmental genomics to this end, technical limitations and conceptual issues still stand in the way of its broad application by end-users. In addition, the multiplicity of potential implementation strategies may contribute to a perception that the routine application of this methodology is premature or "in development", hence restraining regulators from binding these tools into legal frameworks. Here, we review recent implementations of environmental genomics-based methods, applied to the biomonitoring of ecosystems. By taking a general overview, without narrowing our perspective to particular habitats or groups of organisms, this paper aims to compare, review and discuss the strengths and limitations of four general implementation strategies of environmental genomics for monitoring: (a) Taxonomy-based analyses focused on identification of known bioindicators or described taxa; (b) De novo bioindicator analyses; (c) Structural community metrics including inferred ecological networks; and (d) Functional community metrics (metagenomics or metatranscriptomics). We emphasise the utility of the three latter strategies to integrate meiofauna and microorganisms that are not traditionally utilised in biomonitoring because of difficult taxonomic identification. Finally, we propose a roadmap for the implementation of environmental genomics into routine monitoring programmes that leverage recent analytical advancements, while pointing out current limitations and future research needs.

Bacterial diversity and predicted enzymatic function in a multipurpose surface water system - from wastewater effluent discharges to drinking water production

BACKGROUND

Rivers and lakes are used for multiple purposes such as for drinking water (DW) production, recreation, and as recipients of wastewater from various sources. The deterioration of surface water quality with wastewater is well-known, but less is known about the bacterial community dynamics in the affected surface waters. Understanding the bacterial community characteristics -from the source of contamination, through the watershed to the DW production process-may help safeguard human health and the environment.

RESULTS

The spatial and seasonal dynamics of bacterial communities, their predicted functions, and potential health-related bacterial (PHRB) reads within the Kokemäenjoki River watershed in southwest Finland were analyzed with the 16S rRNA-gene amplicon sequencing method. Water samples were collected from various sampling points of the watershed, from its major pollution sources (sewage influent and effluent, industrial effluent, mine runoff) and different stages of the DW treatment process (pre-treatment, groundwater observation well, DW production well) by using the river water as raw water with an artificial groundwater recharge (AGR). The beta-diversity analysis revealed that bacterial communities were highly varied among sample groups (R = 0.92, p <  0.001, ANOSIM). The species richness and evenness indices were highest in surface water (Chao1; 920 ± 10) among sample groups and gradually decreased during the DW treatment process (DW production well; Chao1: 320 ± 20). Although the phylum Proteobacteria was omnipresent, its relative abundance was higher in sewage and industrial effluents (66-80%) than in surface water (55%). Phyla Firmicutes and Fusobacteria were only detected in sewage samples. Actinobacteria was more abundant in the surface water (≥13%) than in other groups (≤3%). Acidobacteria was more abundant in the DW treatment process (≥13%) than in others (≤2%). In total, the share of PHRB reads was higher in sewage and surface water than in the DW treatment samples. The seasonal effect in bacterial communities was observed only on surface water samples, with the lowest diversity during summer.

CONCLUSIONS

The low bacterial diversity and absence of PHRB read in the DW samples indicate AGR can produce biologically stable and microbiologically safe drinking water. Furthermore, the significantly different bacterial communities at the pollution sources compared to surface water and DW samples highlight the importance of effective wastewater treatment for protecting the environment and human health.

Evaluation of V2 18S rDNA barcode marker and assessment of sample collection and DNA extraction methods for metabarcoding of autotrophic euglenids

Even though the interest in metabarcoding in environmental studies is growing, euglenids are still underrepresented in both sea and freshwater bodies researches. The reason for this situation could be the unsuitability of universal eukaryotic DNA barcodes and primers as well as the lack of a verified protocol, suitable to assess euglenid diversity. In this study, using specific primers for the V2 hypervariable region of 18S rDNA for metabarcoding resulted in obtaining a high fraction (85%) of euglenid reads and species‐level identification of almost 90% of them. Fifty species were detected by the metabarcoding method, including almost all species observed using a light microscope. We investigated three biomass harvesting methods (filtering, centrifugation and scraping the side of a collection vessel) and determined that centrifugation and filtration outperformed scrapes, but the choice between them is not crucial for the reliability of the analysis. In addition, eight DNA extraction methods were evaluated. We compared five commercially available DNA isolation kits, two CTAB‐based protocols and a chelating resin. For this purpose, the efficiency of extraction, quality of obtained DNA, preparation time and generated costs were taken into consideration. After examination of the aforementioned criteria, we chose the GeneMATRIX Soil DNA Purification Kit as the most suitable for DNA isolation.

Highlighting the promise of qPCR-based environmental monitoring: response of the ribosomal RNA:DNA ratio of calanoid copepods to toxic cyanobacteria

Calanoid copepods are integral to aquatic food webs and may drive the bioaccumulation of toxins and heavy metals, spread of infectious diseases, and occurrence of toxic cyanobacterial harmful algal blooms (HABs) in freshwater aquatic systems. However, interrelationships between copepod and cyanobacterial population dynamics and ecophysiology remain unclear. Insights into these relationships are important to aquatic resource management, as they may help guide mitigation efforts. We developed a calanoid copepod qPCR assay to investigate how copepod abundance and physiological status relate to the abundance of cyanobacteria and the concentration of total microcystin in a HAB-prone freshwater multi-use eutrophic lake. Through in silico and in vitro validation of primers and analyses of time series, we demonstrate that our assay can be used as a reliable tool for environmental monitoring. Importantly, copepod RNA:DNA ratios on and shortly after the day when microcystin concentration was at its highest within the lake were not significantly lower (or higher) than before or after this period, suggesting that copepods may have been tolerant of microcystin levels observed and capable of perpetuating bloom events by consuming competitors of toxic cyanobacteria.

DNA metabarcoding of zooplankton communities: species diversity and seasonal variation revealed by 18S rRNA and COI

Background

Zooplankton is an important component of aquatic organisms and has important biological and economical significance in freshwater ecosystems. However, traditional methods that rely on morphology to classify zooplankton require expert taxonomic skills. Moreover, traditional classification methods are time-consuming and labor-intensive, which is not practical for the design of conservation measures and ecological management tools based on zooplankton diversity assessment.

Methods

We used DNA metabarcoding technology with two different markers: the nuclear small subunit ribosomal RNA (18S rRNA) and mitochondrial cytochrome c oxidase (COI), to analyze 72 zooplankton samples collected in 4 seasons and 9 locations from the Sanmenxia Reservoir. We investigated seasonal changes in the zooplankton community and their relationship with water environmental factors.

Results

A total of 190 species of zooplankton were found, belonging to 12 phyla, 24 classes, 61 orders, 111 families, and 174 genera. Protozoa, especially ciliates, were the most diverse taxa. Richness and relative abundance of zooplankton showed significant seasonal changes. Both alpha and beta diversity showed seasonal trends: the diversity in summer and autumn was higher than that in winter and spring. The zooplankton diversity was most similar in winter and spring. By correlating metabarcoding data and water environmental factors, we proved that water temperature, chemical oxygen demand, total nitrogen and ammoniacal nitrogen were the main environmental factors driving the seasonal changes in zooplankton in the Sanmenxia Reservoir. Water temperature, followed by total nitrogen, were the most influential factors. This study highlights the advantages and some limitations of zooplankton molecular biodiversity assessment using two molecular markers.

Development of an absolute quantification method for ribosomal RNA gene copy numbers per eukaryotic single cell by digital PCR

Recent increase of Harmful Algal Blooms (HAB) causes world-wide ecological, economical, and health issues, and more attention is paid to frequent coastal monitoring for the early detection of HAB species to prevent or reduce such impacts. Use of molecular tools in addition to traditional microscopy-based observation has become one of the promising methodologies for coastal monitoring. However, as ribosomal RNA (rRNA) genes are commonly targeted in molecular studies, variability in the rRNA gene copy number within and between species must be considered to provide quantitative information in quantitative PCR (qPCR), digital PCR (dPCR), and metabarcoding analyses. Currently, this information is only available for a limited number of species. The present study utilized a dPCR technology to quantify copy numbers of rRNA genes per single cell in 16 phytoplankton species, the majority of which are toxin-producers, using a newly developed universal primer set accompanied by a labeled probe with a fluorophore and a double-quencher. In silico PCR using the newly developed primers allowed the detection of taxa from 8 supergroups, demonstrating universality and broad coverage of the primer set. Chelex buffer was found to be suitable for DNA extraction to obtain DNA fragments with suitable size to avoid underestimation of the copy numbers. The study successfully demonstrated the first comparison of absolute quantification of 18S rRNA copy numbers per cell from 16 phytoplankton species by the dPCR technology.

A strategy for power generation from bilgewater using a photosynthetic microalgal fuel cell (MAFC)

Microbial fuel cells (MFCs) have recently been applied to generate electricity from oily wastewater. Although MFCs that utilize microalgae to provide a self-supporting oxygen (O₂) supply at the cathode have been well discussed, those with microalgae at the anode as an active biomass for treating wastewater and producing electrons are still poorly studied and understood. Here, we demonstrated a bilgewater treatment using single- and double-chamber microalgal fuel cells (SMAFC and DMAFC) capable of generating energy with a novel microalgal strain (Chlorella sorokiniana) that was initially isolated from oily wastewater. Compared to previous MFC studies using green algae, relatively high voltage output (151.3-160.1 mV, 71.3-83.4 mV m^-2 of power density) was observed in the SMAFC under O₂ controlled systems (i.e., acetate addition or light/dark cycle). It was assumed that, under the O₂ depletion, alternative electron acceptors such as bicarbonate may be utilized for power generation. A DMAFC showed better power density (up to 23.9%) compared to the SMAFC due to the separated cathode chamber which fully utilizes O₂ as an electron acceptor. Both SMAFC and DMAFC removed 67.2-77.4% of soluble chemical oxygen demands (SCOD) from the synthetic bilgewater. This study demonstrates that the application of algae-based MFCs is a feasible strategy to treat oil-in-water emulsion while generating electricity.

Aquatic food webs in deep temperate lakes: Key species establish through their autecological versatility

Microbial planktonic communities are the basis of food webs in aquatic ecosystems since they contribute substantially to primary production and nutrient recycling. Network analyses of DNA metabarcoding data sets emerged as a powerful tool to untangle the complex ecological relationships among the key players in food webs. In this study, we evaluated co-occurrence networks constructed from time-series metabarcoding data sets (12 months, biweekly sampling) of protistan plankton communities in surface layers (epilimnion) and bottom waters (hypolimnion) of two temperate deep lakes, Lake Mondsee (Austria) and Lake Zurich (Switzerland). Lake Zurich plankton communities were less tightly connected, more fragmented and had a higher susceptibility to a species extinction scenario compared to Lake Mondsee communities. We interpret these results as a lower robustness of Lake Zurich protistan plankton to environmental stressors, especially stressors resulting from climate change. In all networks, the phylum Ciliophora contributed the highest number of nodes, among them several in key positions of the networks. Associations in ciliate-specific subnetworks resembled autecological species-specific traits that indicate adaptions to specific environmental conditions. We demonstrate the strength of co-occurrence network analyses to deepen our understanding of plankton community dynamics in lakes and indicate biotic relationships, which resulted in new hypotheses that may guide future research in climate-stressed ecosystems.

Protocols for Monitoring Harmful Algal Blooms for Sustainable Aquaculture and Coastal Fisheries in Chile

Harmful algae blooms (HABs) cause acute effects on marine ecosystems due to their production of endogenous toxins or their enormous biomass, leading to significant impacts on local economies and public health. Although HAB monitoring has been intensively performed at spatiotemporal scales in coastal areas of the world over the last decades, procedures have not yet been standardized. HAB monitoring procedures are complicated and consist of many methodologies, including physical, chemical, and biological water sample measurements. Each monitoring program currently uses different combinations of methodologies depending on site specific purposes, and many prior programs refer to the procedures in quotations. HAB monitoring programs in Chile have adopted the traditional microscopic and toxin analyses but not molecular biology and bacterial assemblage approaches. Here we select and optimize the HAB monitoring methodologies suitable for Chilean geography, emphasizing on metabarcoding analyses accompanied by the classical tools with considerations including cost, materials and instrument availability, and easiness and efficiency of performance. We present results from a pilot study using the standardized stepwise protocols, demonstrating feasibility and plausibility for sampling and analysis for the HAB monitoring. Such specific instructions in the standardized protocol are critical obtaining quality data under various research environments involving multiple stations, different analysts, various time-points, and long HAB monitoring duration.

The case for research integration, from genomics to remote sensing, to understand biodiversity change and functional dynamics in the world's lakes

Freshwater ecosystems are heavily impacted by multiple stressors, and a freshwater biodiversity crisis is underway. This realization has prompted calls to integrate global freshwater ecosystem data, including traditional taxonomic and newer types of data (e.g., eDNA, remote sensing), to more comprehensively assess change among systems, regions, and organism groups. We argue that data integration should be done, not only with the important purpose of filling gaps in spatial, temporal, and organismal representation, but also with a more ambitious goal: to study fundamental cross-scale biological phenomena. Such knowledge is critical for discerning and projecting ecosystem functional dynamics, a realm of study where generalizations may be more tractable than those relying on taxonomic specificity. Integration could take us beyond cataloging biodiversity losses, and toward predicting ecosystem change more broadly. Fundamental biology questions should be central to integrative, interdisciplinary research on causal ecological mechanisms, combining traditional measures and more novel methods at the leading edge of the biological sciences. We propose a conceptual framework supporting this vision, identifying key questions and uncertainties associated with realizing this research potential. Our framework includes five interdisciplinary "complementarities." First, research approaches may provide comparative complementarity when they offer separate realizations of the same focal phenomenon. Second, for translational complementarity, data from one research approach is used to translate that from another, facilitating new inferences. Thirdly, causal complementarity arises when combining approaches allows us to "fill in" cause-effect relationships. Fourth, contextual complementarity is realized when together research methodologies establish the wider ecological and spatiotemporal context within which focal biological responses occur. Finally, integration may allow us to cross inferential scales through scaling complementarity. Explicitly identifying the modes and purposes of integrating research approaches, and reaching across disciplines to establish appropriate collaboration will allow researchers to address major biological questions that are more than the sum of the parts.

A Review and Perspective of eDNA Application to Eutrophication and HAB Control in Freshwater and Marine Ecosystems

Changing ecological communities in response to anthropogenic activities and climate change has become a worldwide problem. The eutrophication of waterbodies in freshwater and seawater caused by the effects of human activities and nutrient inputs could result in harmful algae blooms (HABs), decreases water quality, reductions in biodiversity and threats to human health. Rapid and accurate monitoring and assessment of aquatic ecosystems are imperative. Environmental DNA (eDNA) analysis using high-throughput sequencing has been demonstrated to be an effective and sensitive assay for detecting and monitoring single or multiple species in different samples. In this study, we review the potential applications of eDNA approaches in controlling and mitigating eutrophication and HABs in freshwater and marine ecosystems. We use recent studies to highlight how eDNA methods have been shown to be a useful tool for providing comprehensive data in studies of eutrophic freshwater and marine environments. We also provide perspectives on using eDNA techniques to reveal molecular mechanisms in biological processes and mitigate eutrophication and HABs in aquatic ecosystems. Finally, we discuss the feasible applications of eDNA for monitoring biodiversity, surveying species communities and providing instructions for the conservation and management of the environment by integration with traditional methods and other advanced techniques.

Evaluating putative ecological drivers of microcystin spatiotemporal dynamics using metabarcoding and environmental data

Microcystin is a cyanobacterial hepatotoxin of global concern. Understanding the environmental factors that cause high concentrations of microcystin is crucial to the development of lake management strategies that minimize harmful exposures. While the literature is replete with studies linking cyanobacterial production of microcystin to changes in various nutrients, abiotic stressors, grazers, and competitors, no single biotic or abiotic factor has been shown to be reliably predictive of microcystin concentrations in complex ecosystems. We performed random forest regression analyses with 16S and 18S rRNA gene sequencing data and environmental data to determine which putative ecological drivers best explained spatiotemporal variation in total microcystin and several individual congeners in a eutrophic freshwater reservoir. Model performance was best for predicting concentrations of the congener MC-LR, with ca. 88% of spatiotemporal variance explained. Most of the variance was associated with changes in the relative abundance of the cyanobacterial genus Microcystis. Follow-up RF regression analyses revealed that factors that were the most important in predicting MC-LR were also the most important in predicting Microcystis population dynamics. We discuss how these results relate to prevailing ecological hypotheses regarding the function of microcystin.

Seasonality of Planktonic Freshwater Ciliates: Are Analyses Based on V9 Regions of the 18S rRNA Gene Correlated With Morphospecies Counts?

Ciliates represent central nodes in freshwater planktonic food webs, and many species show pronounced seasonality, with short-lived maxima of a few dominant taxa while many being rare or ephemeral. These observations are primarily based on morphospecies counting methods, which, however, have limitations concerning the amount and volume of samples that can be processed. For high sampling frequencies at large scales, high throughput sequencing (HTS) of freshwater ciliates seems to be a promising tool. However, several studies reported large discrepancy between species abundance determinations by molecular compared to morphological means. Therefore, we compared ciliate DNA metabarcodes (V9 regions of the 18S rRNA gene) with morphospecies counts for a 3-year study (Lake Zurich, Switzerland; biweekly sampling, n = 74). In addition, we isolated, cultivated and sequenced the 18S rRNA gene of twelve selected ciliate species that served as seeds for HTS analyses. This workflow allowed for a detailed comparison of V9 data with microscopic analyses by quantitative protargol staining (QPS). The dynamics of V9 read abundances over the seasonal cycle corroborated well with morphospecies population patterns. Annual successions of rare and ephemeral species were more adequately characterized by V9 reads than by QPS. However, numbers of species specific sequence reads only partly reflected rank orders seen by counts. In contrast, biomass-based assemblage compositions showed higher similarity to V9 read numbers, probably indicating a relation between cell sizes and numbers / sizes of macronuclei (or 18S rRNA operons). Full-length 18S rRNA sequences of ciliates assigned to certain morphospecies are urgently needed for barcoding approaches as planktonic taxa are still poorly represented in public databases and the interpretation of HTS data depends on profound reference sequences. Through linking operational taxonomic units (OTUs) with known morphospecies, we can use the deep knowledge about the autecology of these species.

Diverse microbial communities hosted by the model carnivorous pitcher plant Sarracenia purpurea: analysis of both bacterial and eukaryotic composition across distinct host plant populations

BACKGROUND

The pitcher plant Sarracenia purpurea supplements nutrient acquisition through carnivory, capturing insect prey which are digested by a food web community of eukaryotes and bacteria. While the food web invertebrates are well studied, and some recent studies have characterized bacteria, detailed genetic analysis of eukaryotic diversity is lacking. This study aimed to compare eukaryotic and bacterial composition and diversity of pitcher communities within and between populations of host plants in nearby but distinct wetland habitats, and to characterize microbial functions across populations and in comparison with another freshwater community.

METHODS

Pitcher fluid was sampled from the two wetlands, Cedarburg and Sapa Bogs, community DNA was extracted, and 16S and 18S rRNA amplicons were sequenced and data processed for community-level comparisons.

RESULTS AND CONCLUSIONS

Bacterial diversity in the small pitcher volume rivaled that of larger aquatic communities. Between pitcher plant populations, several bacterial families (Kiloniellaceae, Acetobacteraceae, Xanthobacteraceae, Sanguibacteraceae, Oligoflexaceae, Nitrosomonadaceae, Chromatiaceae, Saprospiraceae) were significantly higher in one population. However, although predicted pitcher bacterial functions were distinct from other freshwater communities, especially for some amino acid metabolism, functions were similar across all the pitchers in the two populations. This suggests some functional redundancy among bacterial taxa, and that functions converge to achieve similar food web processes. The sequencing identified a previously under-appreciated high diversity of ciliates, Acari mites, fungi and flagellates in pitcher communities; the most abundant sequences from eukaryotic taxa were Oligohymenophorea ciliates, millipedes and Ichthyosporea flagellates. Two thirds of taxa were identified as food web inhabitants and less than one third as prey organisms. Although eukaryotic composition was not significantly different between populations, there were different species of core taxonomic groups present in different pitchers-these differences may be driven by wetland habitats providing different populations to colonize new pitchers. Eukaryotic composition was more variable than bacterial composition, and there was a poor relationship between bacterial and eukaryotic composition within individual pitchers, suggesting that colonization by eukaryotes may be more stochastic than for bacteria, and bacterial recruitment to pitchers may involve factors other than prey capture and colonization by eukaryotic food web inhabitants.

High-throughput environmental DNA analysis informs a biological assessment of an urban stream

There is growing interest in the use of DNA barcoding and metabarcoding approaches to aid biological assessments and monitoring of waterbodies. While biodiversity measured by morphology and by DNA often has been found correlated, few studies have compared DNA data to established measures of impairment such as multimetric pollution tolerance indices used by many bioassessment programs. We incorporated environmental DNA (eDNA) metabarcoding of seston into a rigorous watershed-scale biological assessment of an urban stream to examine the extent to which eDNA richness and diversity patterns were correlated with multimetric indices and ecological impairment status designations. We also evaluated different filtering approaches and taxonomic classifications to identify best practices for environmental assessments. Seston eDNA revealed a wide diversity of eukaryotic taxa but was dominated by diatoms (36%). Differentiation among sites in alpha and beta diversity was greater when operational taxonomic units (OTUs) were classified taxonomically, but coarse resolution taxonomy (kingdom) was more informative than finer resolution taxonomy (family, genus). Correlations of DNA richness and diversity with multimetric indices for fish and macroinvertebrates were generally weak, possibly because Metazoa were not highly represented in our DNA dataset. Nonetheless, sites could be differentiated based on ecological impairment status, with more impaired sites having lower eDNA diversity as measured by the Shannon index, but higher taxonomic richness. Significant environmental drivers of community structure, as inferred from constrained ordination analyses, differed among kingdoms within the eDNA dataset, as well as from fish and macrobenthos, suggesting that eDNA provides novel environmental information. These results suggest that even a simple seston eDNA filtering protocol can provide biodiversity information of value to stream bioassessment programs. The approach bears further investigation as a potentially useful rapid assessment protocol to supplement more intensive field sampling efforts.