Bioassessment of a Drinking Water Reservoir Using Plankton: High Throughput Sequencing vs. Traditional Morphological Method

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.

Dynamics of Bacterioplankton Communities during Wet and Dry Seasons in the Danjiangkou Reservoir in Hubei, China

Water quality is directly linked to drinking water safety for millions of people receiving the water. The Danjiangkou Reservoir is the main water source for the Middle Route of the South-to-North Water Diversion Project (MR-SNWDP), located in the vicinity of Henan and Hubei provinces in China. Aquatic microorganisms are key indicators of biologically assessing and monitoring the water quality of the reservoir as they are sensitive to environmental and water quality changes. This study aimed to investigate the spatiotemporal variations in bacterioplankton communities during wet (April) and dry (October) seasons at eight monitoring points in Hanku reservoir and five monitoring points in Danku reservoir. Each time point had three replicates, labeled as wet season Hanku (WH), wet season Danku (WD), dry season Hanku (DH), and dry season Danku (DD) of Danjiangkou Reservoir in 2021. High-throughput sequencing (Illumina PE250) of the 16S rRNA gene was performed, and alpha (ACE and Shannon) and beta (PCoA and NDMS) diversity indices were analyzed. The results showed that the dry season (DH and DD) had more diverse bacterioplankton communities compared to the wet season (WH and WD). Proteobacteria, Actinobacteria, and Firmicutes were the most abundant phyla, and Acinetobacter, Exiguobacterium, and Planomicrobium were abundant in the wet season, while polynucleobacter was abundant in the dry season. The functional prediction of metabolic pathways revealed six major functions including carbohydrate metabolism, membrane transport, amino acid metabolism, signal transduction, and energy metabolism. Redundancy analysis showed that environmental parameters greatly affected bacterioplankton diversity during the dry season compared to the wet season. The findings suggest that seasonality has a significant impact on bacterioplankton communities, and the dry season has more diverse communities influenced by environmental parameters. Further, the relatively high abundance of certain bacteria such as Acinetobacter deteriorated the water quality during the wet season compared to the dry season. Our findings have significant implications for water resource management in China, and other countries facing similar challenges. However, further investigations are required to elucidate the role of environmental parameters in influencing bacterioplankton diversity in order to devise potential strategies for improving water quality management in the reservoir.

Dissolved Oxygen and Water Temperature Drive Vertical Spatiotemporal Variation of Phytoplankton Community: Evidence from the Largest Diversion Water Source Area

In order to study the vertical distribution characteristics of phytoplankton in the Danjiangkou Reservoir, the water source of the Middle Route Project of the South-North Water Diversion, seven sampling sites were set up in the Reservoir for quarterly sampling from 2017 to 2019, and water environment surveys were conducted simultaneously. The results showed that 157 species (including varieties) were identified, belonging to 9 phyla and 88 genera. In terms of species richness, Chlorophyta had the largest number of species, accounting for 39.49% of the total species. The Bacillariophyta and Cyanobacteria accounted for 28.03% and 13.38% of the total species, respectively. From the whole Danjiangkou Reservoir, the total phytoplankton abundance varied from 0.09 × 10² to 20.01 × 10⁶ cells/L. In the vertical distribution, phytoplankton were mainly observed in the surface-thermospheric layer (I-II layer) and the bottom layer, while the Shannon-Wiener index showed a trend of gradually decreasing from the I-V layer. The Surfer model analysis showed that there was no significant stratification in the Q site's water temperature (WT) and dissolved oxygen (DO) levels in the water diversion area during the dynamic water diversion process. A canonical correspondence analysis (CCA) showed that DO, WT, pH, electrical conductivity (Cond), chemical oxygen demand (COD_Mn), total phosphorus (TP), ammonia nitrogen (NH₄⁺-N), and total nitrogen (TN) had significant effects on the vertical distribution of phytoplankton (p < 0.05). A partial Mantel analysis showed that the vertical distribution of the phytoplankton community was related to WT, and the phytoplankton community structure at the other sites, except for Heijizui (H) and Langhekou (L), was affected by DO. This study has positive significance for exploring the vertical distribution characteristics of a phytoplankton community in a deep-water dynamic water diversion reservoir.

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.

Comparative analysis of freshwater phytoplankton communities in two lakes of Burabay National Park using morphological and molecular approaches

We analyzed phytoplankton assemblages' variations in oligo-mesotrophic Shchuchie and Burabay lakes using traditional morphological and next-generation sequencing (NGS) approaches. The total phytoplankton biodiversity and abundance estimated by both microscopy and NGS were significantly higher in Lake Burabay than in Lake Shchuchie. NGS of 16S and 18S rRNA amplicons adequately identify phytoplankton taxa only on the genera level, while species composition obtained by microscopic examination was significantly larger. The limitations of NGS analysis could be related to insufficient coverage of freshwater lakes phytoplankton by existing databases, short algal sequences available from current instrumentation, and high homology of chloroplast genes in eukaryotic cells. However, utilization of NGS, together with microscopy allowed us to perform a complete taxonomic characterization of phytoplankton lake communities including picocyanobacteria, often overlooked by traditional microscopy. We demonstrate the high potential of an integrated morphological and molecular approach in understanding the processes of organization in aquatic ecosystem assemblages.

Changes in the Community Structure of Under-Ice and Open-Water Microbiomes in Urban Lakes Exposed to Road Salts

Salinization of freshwater is increasingly observed in regions where chloride de-icing salts are applied to the roads in winter, but little is known about the effects on microbial communities. In this study, we analyzed the planktonic microbiomes of four lakes that differed in degree of urbanization, eutrophication and salinization, from an oligotrophic reference lake with no surrounding roads, to a eutrophic, salinized lake receiving runoff from a highway. We tested the hypothesis that an influence of road salts would be superimposed on the effects of season and trophic status. We evaluated the microbial community structure by 16S rRNA sequencing for Bacteria, and by four methods for eukaryotes: 16S rRNA chloroplast analysis, 18S rRNA sequencing, photosynthetic pigment analysis and microscopy. Consistent with our hypothesis, chloride and total nitrogen concentrations were among the most important statistical factors explaining the differences in taxonomic composition. These factors were positively correlated with the abundance of cryptophytes, haptophytes, and cyanobacteria. Ice-cover was also a major structuring factor, with clear differences between the winter communities and those of the open-water period. Nitrifying and methane oxidizing bacteria were more abundant in winter, suggesting the importance of anaerobic sediment processes and release of reduced compounds into the ice-covered water columns. The four methods for eukaryotic analysis provided complementary information. The 18S rRNA observations were strongly influenced by the presence of ribosome-rich ciliates, but revealed a much higher degree of taxonomic richness and greater separation of lakes, seasonal changes and potential salinity effects than the other methods.

Algal and cyanobacterial diversity in saline rivers of the Elton Lake Basin (Russia) studied via light microscopy and next-generation sequencing

Naturally saline rivers are known in various regions of the world. Saline rivers with a salinity gradient from the source to the mouth are particularly interesting, because the range of salinity is the structure-forming factor of the hydrobiont assemblage. Such rivers are represented by saline rivers of the Elton Lake Basin in Volgograd region of Russia (the Bolshaya Samoroda River and the Malaya Samoroda River). Herein, we analyzed taxonomic structure and species diversity of microalgae and Cyanobacteria of the saline rivers flowing into the Elton Lake by light microscopy and next-generation sequencing. The differences and possible causes of inconsistencies in the results obtained by these methods are discussed. In total, 91 taxa of microorganisms were identified by integrated approach in the assemblages of microalgae and Cyanobacteria in the middle course of the Bolshaya Samoroda River, and 60 taxa – in the river mouth. The species diversity of those assemblages in the hypersaline Malaya Samoroda River was lower: 27 taxa from the middle course and 23 taxa from the mouth. Next-generation sequencing allowed us to refine and expand the list of microalgae taxa in the studied saline rivers due to detection of species which were hard to identify, low-abundance taxa, as well as extremely small-cell forms. Some discrepancies between the data obtained by light microscopy and next-generation sequencing indicate the advantage of simultaneous use of both methods for study of the algae communities. Such a comprehensive approach provides the most accurate and correct list of taxa added with the morphological descriptions and 18S rRNA and 16S rRNA partial sequences. Generally, 18 taxa have been recorded for the first time in the Bolshaya Samoroda River, belonging to the phyla Chlorophyta (Borodinellopsis sp., Chlorochytrium lemnae Cohn, Caespitella sp., Halochlorococcum sp., Tetraselmis cordiformis (H. J. Carter) F. Stein), Ochrophyta (Pseudocharaciopsis ovalis (Chodat) D. J. Hibberd, Characiopsis sp., Poterioochromonas stipitata Scherffel, Chrysolepidomonas sp.), Euglenozoa (Euglena bucharica I. Kisselev, Lepocinclis tripteris (Dujardin) B. Marin & Melkonian, Phacus orbicularis K. Hübner, P. parvulus G. A. Klebs), Cryptophyta (Hemiselmis cryptochromatica C. E. Lane & J. M. Archibald, Rhodomonas sp., Hanusia phi J. A. Deane), Haptophyta (Pavlova sp.), Cyanobacteria (Johanseninema constrictum (Szafer) Hasler, Dvorák & Poulícková). Seven taxa have been detected for the first time in the algal and cyanobacterial assemblages of the Malaya Samoroda River from the phyla Chlorophyta (Tetraselmis cordiformis, T. arnoldii (Proschkina-Lavrenko) R. E. Norris, Hori & Chihara, T. tetrathele (West) Butcher, Pyrobotrys elongatus Korshikov), Cryptophyta (Hanusia phi), and Cyanobacteria (Synechococcus elongatus (Nägeli) Nägeli, Oscillatoria simplicissima Gomont).

Structural Characteristics and Driving Factors of the Planktonic Eukaryotic Community in the Danjiangkou Reservoir, China

Planktonic eukaryotes are widespread in aquatic ecosystems, and the study of their community composition and driving factors is of great significance to protecting and maintaining the balance of these ecosystems. This study evaluates five typical ecological sites in the Danjiangkou Reservoir—the water source for the project. This was done to comprehensively understand the composition of Danjiangkou Reservoir planktonic eukaryotes, and ensure the ecological balance of the water source for the South-to-North Water Diversion Project. The diversity of the planktonic eukaryotes in surface water and the factors driving changes in their abundance are analyzed with an 18S ribosomal DNA sequencing approach. Monitoring shows that the Danjiangkou Reservoir has good water quality. The Danjiangkou Reservoir planktonic eukaryote community is mainly composed of 11 phyla, of which Cryptomonadales is dominant, accounting for an average percentage of 65.19% of the community (47.2–84.90%). LEFSe analysis shows significant differences among samples in the abundances of 13 phyla, 20 classes, 23 orders, 26 families, and 27 genera, and there are also significant differences in the diversity of planktonic eukaryotes at different temporal and spatial scales. Redundancy analysis (RDA) show that water temperature, DO, SD, TN, and Chla are significant factors that affect the composition of the planktonic eukaryote community. Spearman rank correlation analysis combined with taxonomic difference analysis shows that Kathablepharidae and Choanoflagellida are not sensitive to environmental or physicochemical factors and that the interannual variations in their abundance are not significant. Network analysis shows that Protalveolata, Basidiomycota, P1-31, Bicosoecida, and Ochrophyta represent important nodes in the single-factor network, while Chytridiomycota, P1-31, Cryptomycota, Ochrophyta, Ichthyosporea, Bicosoecida, Protalveolata, and physicochemical factors (ORP, TN, WT, DO, SD, NH3-N, and NO3-N) represent important nodes in the two-factor network.

Bacterial and Eukaryotic Small-Subunit Amplicon Data Do Not Provide a Quantitative Picture of Microbial Communities, but They Are Reliable in the Context of Ecological Interpretations

High-throughput sequencing (HTS) of amplified fragments of rRNA genes provides unprecedented insight into the diversity of prokaryotic and eukaryotic microorganisms. Unfortunately, HTS data are prone to quantitative biases, which may lead to an erroneous picture of microbial community composition and thwart efforts to advance its understanding. These concerns motivated us to investigate how accurately HTS data characterize the variability of microbial communities, the relative abundances of specific phylotypes, and their relationships with environmental factors in comparison to an established microscopy-based method. We compared results obtained by HTS and catalyzed reporter deposition-fluorescence in situ hybridization (CARD-FISH) from three independent aquatic time series for both prokaryotic and eukaryotic microorganisms (almost 900 data points, the largest obtained with both methods so far). HTS and CARD-FISH data disagree with regard to relative abundances of bacterial and eukaryotic phylotypes but identify similar environmental drivers shaping bacterial and eukaryotic communities.

High-throughput sequencing (HTS) of gene amplicons is a preferred method of assessing microbial community composition, because it rapidly provides information from a large number of samples at high taxonomic resolution and low costs. However, mock community studies show that HTS data poorly reflect the actual relative abundances of individual phylotypes, casting doubt on the reliability of subsequent statistical analysis and data interpretation. We investigated how accurately HTS data reflect the variability of bacterial and eukaryotic community composition and their relationship with environmental factors in natural samples. For this, we compared results of HTS from three independent aquatic time series (n = 883) with those from an established, quantitative microscopic method (catalyzed reporter deposition-fluorescence in situ hybridization [CARD-FISH]). Relative abundances obtained by CARD-FISH and HTS disagreed for most bacterial and eukaryotic phylotypes. Nevertheless, the two methods identified the same environmental drivers to shape bacterial and eukaryotic communities. Our results show that amplicon data do provide reliable information for their ecological interpretations. Yet, when studying specific phylogenetic groups, it is advisable to combine HTS with quantification using microscopy and/or the addition of internal standards.

IMPORTANCE High-throughput sequencing (HTS) of amplified fragments of rRNA genes provides unprecedented insight into the diversity of prokaryotic and eukaryotic microorganisms. Unfortunately, HTS data are prone to quantitative biases, which may lead to an erroneous picture of microbial community composition and thwart efforts to advance its understanding. These concerns motivated us to investigate how accurately HTS data characterize the variability of microbial communities, the relative abundances of specific phylotypes, and their relationships with environmental factors in comparison to an established microscopy-based method. We compared results obtained by HTS and catalyzed reporter deposition-fluorescence in situ hybridization (CARD-FISH) from three independent aquatic time series for both prokaryotic and eukaryotic microorganisms (almost 900 data points, the largest obtained with both methods so far). HTS and CARD-FISH data disagree with regard to relative abundances of bacterial and eukaryotic phylotypes but identify similar environmental drivers shaping bacterial and eukaryotic communities.

Rhizosphere Bacterial Community Structure and Predicted Functional Analysis in the Water-Level Fluctuation Zone of the Danjiangkou Reservoir in China During the Dry Period

The water-level fluctuation zone (WLFZ) is a transitional zone between terrestrial and aquatic ecosystems. Plant communities that are constructed artificially in the WLFZ can absorb and retain nutrients such as nitrogen (N) and phosphorus (P). However, the microbial community composition and function associated with this process have not been elucidated. In this study, four artificially constructed plant communities, including those of herbs (Cynodon dactylon and Chrysopogon zizanioides), trees (Metasequoia glyptostroboides), and shrubs (Salix matsudana) from the newly formed WLFZ of the Danjiangkou Reservoir were evaluated. The bacterial community compositions were analyzed by 16S rRNA gene sequencing using a MiSeq platform, and the functions of these communities were assessed via Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) analysis. The results showed that the bacterial communities primarily comprised 362 genera from 24 phyla, such as Proteobacteria, Acidobacteria, Actinobacteria, and Gemmatimonadetes, showing the richness of the community composition. Planting altered the bacterial community composition, with varying effects observed among the different plant types. The bacterial community functional analysis revealed that these bacteria were primarily associated with six biological metabolic pathway categories (e.g., metabolism, genetic information processing, and environmental information processing) with 34 subfunctions, showing the richness of community functions. The planting of M. glyptostroboides, S. matsudana, and C. dactylon improved the metabolic capabilities of bacterial communities. N- and P-cycling gene analysis showed that planting altered the N- and P-cycling metabolic capacities of soil bacteria. The overall N- and P-metabolic capacity was highly similar between C. dactylon and C. zizanioides samples and between S. matsudana and M. glyptostroboides samples. The results of this study provide a preliminary analysis of soil bacterial community structure and function in the WLFZ of the Danjiangkou Reservoir and provides a reference for vegetation construction in this zone.

Illumina MiSeq sequencing and network analysis the distribution and co-occurrence of bacterioplankton in Danjiangkou Reservoir, China

Network analysis has contributed to studies of the interactions of microorganisms and the identification of key populations. However, such analysis has rarely been conducted in the study of reservoir bacterioplankton communities. This study investigated the bacterioplankton community composition in the surface water of the Danjiangkou Reservoir using the Illumina MiSeq sequencing platform. We observed that the bacterioplankton community primarily consisted of 27 phyla and 336 genera, including Actinobacteria, Proteobacteria, and Bacteroidetes, demonstrating the richness of the community composition. Redundancy analysis of the bacterioplankton communities and environmental variables showed that the total nitrogen (TN), pH, chemical oxygen demand (COD), and permanganate index (COD_Mn) were important factors affecting the bacterioplankton distribution. Network analysis was performed using the relative abundances of bacterioplankton based on the phylogenetic molecular ecological network (pMEN) method. The connectivity of node i within modules (Zi), the connectivity of node i among modules (Pi), and the number of key bacteria were high at the Taizishan and Heijizui sites, which were associated with higher TN contents than at the other sites. Among the physicochemical properties of water, TN, ammonia nitrogen (NH₄-N), pH, COD, and dissolved oxygen (DO) might have great influences on the functional units of the bacterial communities in bacterioplankton molecular networks. This study improves the understanding of the structure and function of bacterioplankton communities in the Danjiangkou Reservoir.

Planktonic fungal community structures and their relationship to water quality in the Danjiangkou Reservoir, China

Planktonic fungi are important components of aquatic ecosystems, and analyses of their community composition and function have far-reaching significance for the ecological management and maintenance of reservoir environments. However, few studies have investigated the composition, distribution, and function of planktonic fungi in reservoir ecosystems and their relationship with water quality. Here, the composition of the planktonic fungal community in the surface water layer of the Danjiangkou Reservoir is investigated using Illumina MiSeq sequencing. According to the results, the reservoir community is primarily composed of 7 phyla, including Ascomycota, Rozellomycota, Basidiomycota, Chytridiomycota, and Zygomycota, comprising 294 genera, demonstrating the rich diversity of this community. Redundancy analysis (RDA) of the planktonic fungal community and environmental factors showed dissolved oxygen (DO), chemical oxygen demand (COD), total nitrogen (TN), chlorophyll a (Chl a), and permanganate (CODMn) to be important factors influencing the distribution of planktonic fungi. Spearman correlation analysis of the planktonic fungal community composition and diversity indices with physical and chemical water quality parameters showed that the impacts of TN, COD and DO were the most significant. The results of this study on the planktonic fungal community in the Danjiangkou Reservoir area using high-throughput sequencing revealed that the community is sensitive to water quality parameters. This result provides a reference for studying the composition and distribution of the planktonic fungal community in Danjiangkou Reservoir and its role in the biogeochemical cycle.