Next Article in Journal
Child Nutrition Trends Over the Past Two Decades and Challenges for Achieving Nutrition SDGs and National Targets in China
Next Article in Special Issue
Mechanism of Membrane Fouling Control by HMBR: Effect of Microbial Community on EPS
Previous Article in Journal
HGBEnviroScreen: Enabling Community Action through Data Integration in the Houston–Galveston–Brazoria Region
Previous Article in Special Issue
Affect Path to Flood Protective Coping Behaviors Using SEM Based on a Survey in Shenzhen, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Community Compositions of Phytoplankton and Eukaryotes during the Mixing Periods of a Drinking Water Reservoir: Dynamics and Interactions

1
Shaanxi Key Laboratory of Environmental Engineering, Key Laboratory of Northwest Water Resource, Environment and Ecology, MOE, Xi’an University of Architecture and Technology, Xi’an 710055, China
2
Guangdong Key Laboratory of Integrated Agro-Environmental Pollution Control and Management, Guangdong Institute of Eco-Environmental Science and Technology, Guangzhou 510650, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2020, 17(4), 1128; https://doi.org/10.3390/ijerph17041128
Submission received: 3 December 2019 / Revised: 5 February 2020 / Accepted: 7 February 2020 / Published: 11 February 2020
(This article belongs to the Special Issue Water Resources and Hydro-Ecology)

Abstract

:
In deep drinking water reservoir ecosystems, the dynamics and interactions of community compositions of phytoplankton and eukaryotes during the mixing periods are still unclear. Here, morphological characteristics combined with high-throughput DNA sequencing (HTS) were used to investigate the variations of phytoplankton and the eukaryotic community in a large canyon-shaped, stratified reservoir located at the Heihe River in Shaanxi Province for three months. The results showed that Bacillariophyta and Chlorophyta were the dominant taxa of the phytoplankton community, accounting for more than 97% of total phytoplankton abundance, which mainly consisted of Melosira sp., Cyclotella sp., and Chlorella sp., respectively. Illumina Miseq sequencing suggested that the biodiversity of eukaryotes increased over time and that species distribution was more even. Arthropoda (6.63% to 79.19%), Ochrophyta (5.60% to 35.16%), Ciliophora (1.81% to 10.93%) and Cryptomonadales (0.25% to 11.48%) were the keystone taxa in common, contributing over 50% of the total eukaryotic community. Cryptomycota as a unique fungus was observed to possess significant synchronization with algal density, reaching a maximum of 10.70% in December (when the algal density distinctly decreased) and suggesting that it might affect the growth of algae through parasitism. Co-occurrence network patterns revealed the complicated and diverse interactions between eukaryotes and phytoplankton, suggesting that eukaryotes respond to variations in dynamic structure of the phytoplankton community, although there might be antagonistic or mutualistic interactions between them. Redundancy analysis (RDA) results showed that environmental variables collectively explained a 96.7% variance of phytoplankton and 96.3% variance of eukaryotic microorganisms, indicating that the temporal variations of phytoplankton and eukaryotic microorganisms were significantly affected by environmental conditions. This study shows that potential interactions exist between phytoplankton and eukaryotic microorganism communities, andcould improve our understanding of the ecological roles of phytoplankton and eukaryotic microorganisms in changing aquatic ecosystems. However, long-term investigations are necessary in order to obtain comprehensive understandings of their complicated associations.

Graphical Abstract

1. Introduction

Reservoirs have many important ecosystem functions such as being a major source of drinking water for local inhabitants and providing agriculture irrigation, flood control, and power generation [1,2]. Reservoir-ecosystem state shifts could threaten biodiversity and cause substantial losses of sustainable ecosystem goods and essential ecosystem services for human beings [3,4]. In recent years, eutrophications, and especially the formation of algal blooms, have become increasingly severe in reservoirs, and have caused periodic-rising-water quality deterioration [5,6,7]. Therefore, adequate management is essential and of utmost importance in reservoirs. A reservoir is an ideal ecological-ecosystem model to explore the physical, chemical, and biological characteristics of aquatic environments [3]. Vertical profiles of water quality are influenced by physical thermal stratification [2,3], and the drivers (e.g., temperature, light, nutrients) of algal blooms in oligotrophic drinking water reservoir shave been investigated in the past few decades. However, fewer studies have focused on the dynamics of phytoplankton and eukaryotes during mixing periods of drinking water reservoirs.
Algal blooms are widespread, and are an ongoing global problem that can affect aquatic organisms and alter water quality in aquatic ecosystems [8,9]. The dynamics of the phytoplankton population are traditionally considered to be largely influenced by abiotic variables, including nutrients concentrations, water temperature, adequate light intensity, pH values, and CO2 [10,11]. The biological interactions of regulating algal blooms are only just beginning to be identified. With the continuously deepening studies of freshwater-system ecological evaluation and the outbreak of mechanism of algal blooms, it is increasingly recognized that the formation of algal blooms is closely related to aquatic microbes such as bacterial and fungal communities [1,10]. These microorganisms, as a major contributor to primary productivity, play a critical role in aquatic food webs, nutrient cycling, global biogeochemical cycles, and oxygen release [12,13]. However, our knowledge about biotic factors is far from holistic, and lacks an understanding of the correlation between phytoplankton and microbial eukaryotic community composition, especially instratified drinking water reservoirs [14,15].
Phytoplankton and water eukaryotes occupy an important position in the structure and function of freshwater ecosystems, which are essential components of aquatic food webs [16]. Phytoplankton supplies oxygen and nutrients for grazers, whereas microbes decompose organic matter, recycle nutrients, and maintain the balance of ecosystems [17,18,19]. In recent years, the mutual interaction between microbes and algal bloom species has been recognized as a vital factor regulating the dynamics of microbes and algal populations [1,8,9,19,20,21,22,23,24,25,26,27,28,29]. Many interaction processes, including mutualism, algicidal effects, allelopathy, and parasitism between microorganisms and algae, can affect the formation of microbial communities [25]. Over the past few years, extensive studies have investigated bacterial [1,8,9,12,13,16,18,19,20,21,22,23,25,26,27,28,30], fungal [1,11,24,29], and zooplankton [31] community structures associated with phytoplankton dynamics using different molecular approaches, including denaturing gradient gel electrophoresis, terminal restriction fragment length polymorphism, and the Illumina Miseq DNA sequencing technique in various aquatic ecosystems such as lakes [8,9,13,16,18,19,20,31], rivers [12,28], oceans [23,24,26,30], and reservoirs [1,27]. They found that phytoplanktonis strongly linked with symbiotic microbes, which area biological disturbance to microbial communities that alter their surrounding environmental characteristics. Several recent studies have demonstrated a strong effect of phytoplankton blooms on eukaryotic plankton, which can remarkably alter the composition of eukaryotic microorganism communities [14,15,25]. Compared with bacterial, the community compositions and dynamics of eukaryotic plankton have relatively seldom been investigated in drinking water reservoirs. As primary producers and consumers, eukaryotic plankton plays an essential role in microbial food webs, which are crucial for the recycling of nutrients [14]. Therefore, their response to environmental fluctuations associated with algae may directly affect the function of aquatic ecosystems [15,32].
The ecological-network approach, a powerful bioinformatics-analysis technology, has been widely used to explore the coexistence and interaction of microbial communities [33,34] harbored in soils [35,36], sludge [37], compost [38], urban lakes [13,21,39] and marine environments [40]. In recent years, this approach has been developed and widely applied to explore the complicated mutual correlations among microbial communities [13,15,16,25], providing a unique and deep perspective on microbial interactions and ecological rules. The network contains a set of nodes or vertices that are connected in pairs by lines or edges, and it has the advantage of being able to simplify molecular ecological networks [13]. Recently, Zhang et al. [37] used Illumina Miseq DNA sequencing to reveal nirS-type denitrifying bacterial communities from different wastewater treatment plants, and to determine the coexistence and interaction of dominant and rare populations using network. Zhao et al. [13] investigated and compared the environmental niches of bacterioplankton in six urban lakes using network analysis, which demonstrated that co-occurrence networks varied in different seasons and that the ecological roles of cyanobacteria differed from season to season. They believed that strong co-occurrence patterns and appropriate living conditions could be critical reasons for the dominance of cyanobacteria in summer. However, from an ecological viewpoint, the dynamics and interactions of the community compositions of phytoplankton and eukaryotes during the mixing periods of drinking water reservoirs are still not comprehensively understood.
To this end, we used a high-throughput sequencing technique (Illumina Miseq DNA sequence) to investigate the composition of the water-eukaryote community in Jinpen drinking water reservoir. A traditional microscopic method was employed to analyze the abundance and composition of the phytoplankton community. Moreover, correlation-based network analysis was utilized to investigate the specific associations between these two communities. The main objectives of our present study were as follows: (1) explore the temporal variation of the phytoplankton community from October to December 2018; (2) evaluate the eukaryotic plankton community composition and structure through a high-throughput DNA sequencing technique; and (3) assess the coexistence and interaction of phytoplankton and eukaryotic microbe compositions, as well as their relationships with water-quality parameters. This study will enhance our understanding of the potential role of microbial diversity as an important driver for algal blooms in freshwater ecosystems.

2. Experiments

2.1. Study Sites and Sampling

Field sampling was conducted at Jinpen reservoir (34°42′–34°13′ N; 107°43′–108°24′ E) [2], located in the city of Xi’an, Shaanxi province, Northwest China. It is a canyon-shaped reservoir with strong thermal stratification (during summer and autumn), and serves as an important drinking water source with a daily water-supply capacity of 8.0 × 105 m3, which accounts for 76% of the total water supply for about 8.58 million inhabitants of Xi’an [2,7,41]. The total storage capacity of Jinpen reservoir is nearly 2.0 × 108 m3 [2,7,41]. It has a water surface area of approximately 4.55 km2 and a catchment area of 1418 km2, with a mean depth of 40 m and a maximum depth of 95 m [2,7]. Annual average rainfall, evaporation and runoff volume of Jinpen reservoir are 898 mm, 948.5 mm, and 545.22 million m3, respectively [2].
During the mixing periods (Figure 1), water samples were collected on a monthly basis between October and December in 2018 from different depths (0, 2, 5, 10 m) of three stations (coded as sites A, B, and C). Sites A (34°02′46″ N; 108°12′24″ E) and B (34°02′24″ N; 108°11′48″ E) are located in the main reservoir area; site A is near the intake tower for the drinking water plants, while site C (34°00′51″ N; 108°10′39″ E) is located in the upstream river area. Thus, the three sampling sites could represent the majority of cases in the entire reservoir according to the results of our research group survey [2,7]. The specific information on water level, average daily inflow and outflow, rainfall, sampling times, and depth in the Jinpen drinking water reservoir from October to December 2018 are shown in Figure 1.
All samples were collected in triplicate by a sterilized water-sampling device [2,41], stored in sterile polyethylene containers (2.0 L), kept at 4 °C, and immediately transferred to the laboratory within 24 h for further processing. Each of the samples was 2.0 L and used for chemical parameter, algal biomass, and water microbial community analyses, respectively. For microbial community analysis, 500 mL water from each sample was filtered through a 0.22 μm pore-size polycarbonate filter (Millipore, USA) and filters containing microbes were frozen at −20 °C immediately until deoxyribonucleic acid (DNA) extraction.

2.2. Measurement of Water-Quality Parameters

The physical characteristics of the samples, namely water temperature (WT), dissolved-oxygen (DO) concentration, pH, water turbidity, electrical conductivity (Cond.), and chlorophylla (Chl-a) were simultaneously recorded in situ by a Multi-Parameter-Water Quality Sonde (Hydro-lab DS5, HACH Company, Loveland, CO, USA). The chemical analyses of the water samples, total nitrogen (TN), ammonium nitrogen (NH4+-N), nitrate nitrogen (NO3-N), and total phosphorus (TP) were measured in accordance with standard methods [42,43,44]. Briefly, concentrations of NH4+-N and NO3-N were determined following Nessler’s reagent spectrophotometry and ultraviolet spectrophotometric methods, respectively. TN and TP concentrations were quantified according to the standard UV spectrophotometric method after digestion (121 °C, 30 min; DR5000, HACH, USA) [7]. The concentrations of dissolved organic carbon (DOC) were measured with a TOC analyzer (TOC-L CPN, Japan). Concentrations of iron (Fe) and manganese (Mn) were examined using flame atomic absorption spectrometry (FAAS, AA 6800, Shimadzu, Japan) [1,2,42]. All the physicochemical analyses were performed within 48 h, and all parameters of this study were the averages of duplicate sample measurements.

2.3. Phytoplankton Identification and Counting

For phytoplankton analysis, 500 mL of the water sample was collected, preserved with 1% acidic Lugol’s solution, and settled for 48 h. After sedimentation for 48 h, samples were concentrated to 10 mL in order to determine algal-cell (AC) concentrations. Then, 100 μL of the concentrated sample was enumerated using a counting chamber under 400× microscope magnification after gently shaking, and recorded in units of ×104 cells per liter. Each sample was generally counted twice (if the results of the two counts differed from the average by more than 15%, a third enumeration was required), then averaged to obtain a mean value for the sample [43,44]. Phytoplankton identification was performed according to Hu et al. [43], supplemented with other references [5,16,18], and identified to the lowest possible taxonomic level. All phytoplankton were enumerated and identified using a Nikon microscope (50i, Nikon, Japan).

2.4. Water Microbial DNA Extraction

For each site, total DNA from water microbes on the filters was extracted and purified using a Microbial DNA Kit (Omega Bio-Tek, CA, USA) according to the manufacturer’s protocol. The concentration and quality of genomic DNA were determined through spectrophotometric analysis using a Nano-Drop ND-2000 spectrophotometer (Thermo Scientific, Waltham, MA, USA) [45,46]. Then, the extracted DNA samples were stored at −20 °C prior to further analysis [19].

2.5. Illumina MiSeq Sequencing and Sequence Analysis

The Illumina Miseq sequencing technique was conducted to assess the water eukaryotic microbe communities. The 18S rRNA V9 region, a hypervariable area of the eukaryotic 18S rRNA gene, was amplified for analysis of the eukaryotic microorganism community using eukaryote-specific primers 1380F (5′-CCCTGCCHTTTGTACACAC-3′) and 1510R (5′-CCTTCYGCAGGTTCACCTAC-3′) [15]. All plankton DNA samples were PCR-amplified in triplicate with a PCR thermal cycler (ABI GeneAmp 9700, CA, USA). The PCR reaction mixture (30 μL) consisted of 15 μL of Phusion High-Fidelity PCR Master Mix (New England Biolabs, Beverly, MA, USA), 0.2 μM each of the forward and reverse primers, and about 10 ng of template DNA. The PCR reaction consisted of an initial denaturation for 1 min at 98 °C, followed by 30 cycles of 10 s denaturation at 98 °C, annealing at 50 °C for 30 s, elongation at 72 °C for 30 s, and a final 5 min of extension at 72 °C. The triplicate PCR products of each sample were pooled in equal quantity for gel purification, and paired-end sequencing (2 × 300 bp) was conducted on the Illumina HiSeq platform (Illumina Inc., San Diego, CA, USA) [47] at Shanghai Biozeron Bio-Technology Co., Ltd. (Shanghai, China).

2.6. Nucleotide-Sequence Accession Numbers

All Illumina MiSeq DNA sequence data in this study were deposited in the public NCBI Sequence Read Archive (SRA) database (http://www.ncbi.nlm.nih.gov/) under Bio-Project number: PRJNA531813 and accession numberSRP192061.

2.7. Statistical Analysis

Two-way ANOVA was employed with SPSS software (version 17.0, SPSS Inc., Chicago, IL, USA) using a Tukey’s HSD posthoc test (p-value < 0.05) to assess statistical significance among water-quality parameters, algal-cell concentration, and algal biomass under different spatiotemporal conditions.
To better understand the interactions between phytoplankton and eukaryotes, we utilized network analysis methods to analyze those two communities. In a network model, the nodes represent tax awhile the connecting edges represent correlation over time. The correlation network indicates the interactions between the abundant genera with significant (p < 0.05) and significantly strong (|R| > 0.8) Spearman correlations. The size of each node indicates the number of correlations. The thickness of each edge indicates the strength of the Spearman’s correlation coefficient. A network clustering diagram was generated by Cytoscape Software (version 2.8.0, http://cytoscape.org). The raw Illumina MiSeq DNA sequences were processed and checked using Quantitative Insights into Microbial Ecology (QIIME2, version 1.9.1, https://qiime2.org/). Water eukaryotic microorganism community sequence reads were filtered by the QIIME pipeline, and trimmed when they were shorter than 50 bp [15]. After assembling and quality filtering, sequences with ≥97% similarity were assigned to the same operational taxonomic units (OTUs) that were used to compute the Chao richness estimator (Chao1), Shannon diversity (H’), and Simpson diversity (D) index using the MOTHUR package (version 1.22.2, http://www.mothur.org) [1,15]. The most representative sequences for each OTU were selected, and we then used the Ribosomal Database Project (RDP) classifier (RDP Release 11.5, https://rdp.cme.msu.edu/) to annotate taxonomic information for each representative sequence [1]. Heat-map profiles were created by R Software (version 3.2.3, R core team, Vienna, Austria) [19]. To identify of the drivers of dominance in phytoplankton and eukaryote communities, redundancy analysis (RDA) was implemented on the basis of CANOCO software (version 4.5, Wageningen, The Netherlands) to explore the relationships that were consistent between eukaryotic-plankton communities and environmental variables. Prior to analysis, all water quality indicators were log(x+1)-transformed except for pH to meet normality and homoscedasticity using a Monte Carlo permutation test (p-value < 0.05).

3. Results and Discussion

3.1. Variation of Water Quality Parameters

The trends of physicochemical variables in Jinpen reservoir from October to December 2018 exhibited distinct types of temporal variation, as shown in Table 1. During the studied periods, the highest temperature was observed in October, with a maximum of 17.11 ± 0.10 °C. DO concentration declined from October to November, after which it significantly increased in December (p < 0.001), while pH was observed to have the same variation over the studied periods (p < 0.01). Generally, DOC concentration (p < 0.001) and conductivity (p < 0.001) gradually increased and peaked in November (at a depth of 5 m) and December (at a depth of 0 m) alone, with a maximum of 3.30 ± 0.14 mg/L and 157.33 ± 4.16 S/cm, respectively. The highest concentrations of Fe and Mn were both observed in December, with values of 0.038 ± 0.00 and 0.024 ± 0.00 mg/L (Table 1). A previous study conducted by Landa et al. [48] indicated that a higher water iron concentration could contribute to the proliferation of algal blooms in oceans. However, recently, Zhang et al. [1] have found that the concentration of Fe was at its lowest before the algal blooms, and that it increased markedly during the decline periods (p < 0.001). Compared with DOC and conductivity, concentrations of TN, NH4+-N and TP exhibited inverse trends, with the highest concentrations (as well as the highest concentrations of algal density) mainly observed in October. Recently, Su et al. [49] investigated the effect of nutrient limitation on phytoplankton growth, and found that the combination of N and P exhibited a more marked and positive influence on phytoplankton proliferation than N or P addition did independently in Spring Lake, USA. NO3-N concentrations were in the range of 0.87 ± 0.00–0.94 ± 0.06 mg/L, and were not significantly different throughout the whole periods. Furthermore, trends of physicochemical factors were mainly consistent across all four depths, with the exception of turbidity. It was revealed that the spatial variation of all water quality parameters from October to December in Jinpen reservoir was not significant, resulting from Jinpen reservoir being in a mixed state during that period, with a similar vertical stratification of its physical and chemical factors [50]. Recently, Zhou et al. [51] conducted field research onthe impact of a water-lifting aeration system on water quality and microbial community structure in a drinking water reservoir, and found that reservoir stratification was successfully eliminated due to the artificial mixing, which strongly affected environmental variables and obviously diminished differences of water quality between the bottom and surface. Our results indicate that water quality parameters varied a lot with the changing environment, which also significantly altered the community structure of phytoplankton and eukaryotes in Jinpen reservoir. However, changing environmental conditions are in complete extrinsic factors for predicting the dynamics of phytoplankton and eukaryote communities.

3.2. Shifts in Algal Cell Concentration and Chlorophyll a

Variations of algal-cell concentrations and chlorophyll a from October to December in 2018 are shown in Figure 2. The peaks of chlorophyll a (3.49 μg/L) and algal-cell concentrations (303.5 × 104 cells/L) both appeared in October over the sampling period (p < 0.05). After that, cell concentration first decreased from 303.5 × 104 to 212.5 × 104 cells/L (p < 0.05) and finally exhibited relatively stable trends with concentrations ranging from 200.3 × 104 to 212.5× 104 cells/L (Figure 2A), which were generally consistent with changes in chlorophyll a (Figure 2B). Cell concentrations and chlorophyll a also exhibited significant differences (p < 0.05) from the depth of 0 to 10 m in October, which gradually declined (Figure 2). In contrast, both metrics exhibited less vertical variation in November and December than in October, and tended to be consistent. Results of this research indicate that phytoplankton growth in winter was not as active as in other seasonson the basis of an observed low algal biomass. Additionally, algal-cell concentrations and chlorophyll a were consistent with change of depth from November to December. Some explanations are as follows. On the one hand, water temperature is the most important environmental factor affecting growth, development, community structure, quantity change, and the horizontal and vertical distribution of microbes in aquatic ecosystems [11]. On the other hand, the reservoir was at the natural mixed state during the studied period, which has an extreme influence on the entire functioning of the aquatic ecosystem [50,52]. Thus, phytoplankton in reservoirs are affected when the temperature of winter is the lowest of the whole year. Bullerjahn et al. [53] investigated the dynamics and functional roles of algal and bacterial communities during winter in two seasonally ice-covered central European Great Lakes, observing that the algal biomasses of the two lakes during winter were lower than its means during a long-term summer. Diatoms, green algae, cryptophytes and chrysoflagellates were the dominant groups of algal communities, accounting for the majority of algal biomass.

3.3. Phytoplankton Dynamics and Community Structure

In the present study, a total of 31 phytoplankton taxa were detected on the basis of morphological identification by Nikon microscopy (Figure S1). They were derived from Bacillariophyta, Chlorophyta, Cryptophyta, Cyanophyta, Dinophyta, Euglenophyta and Chrysophyta. The phytoplankton community was mainly dominated by species belonging to Bacillariophyta and Chlorophyta, accounting for more than 97% of the total phytoplankton abundance based on algal-cell concentration throughout the studied periods (Figure 3). With time, Chlorophyta decreased in relative abundance from 62.50% to 39.28% of the total phytoplankton community, whereas Bacillariophyta followed the completely inverse pattern, with a relative abundance increasing from 39.86% to 62.80%. The reservoirs showed an apparent switch from a Chlorophyta-dominated state to a Bacillariophyta taxa-dominated state. In winter, picoeukaryotes, such as common members of diatoms, green algae, and small flagellates, usually dominate the phytoplankton community [53]. Diatoms are common species in all kinds of fresh-water bodies, and widely exist in rivers, reservoirs, and lakes [54]. Previous studies have found that diatoms are very sensitive to changes in water temperature, nutrients, flow rate, predation by zooplankton, and other environmental fluctuations, which have led to corresponding variations in diatom species and densities [54]. Moreover, with advantages in turbidity and mixing environments, they are the main phytoplankton phylum found in aquatic ecosystems, and are often used as indicators for assessing ecological states such as water-quality degradation and hydrological and climate changes [12,54]. Thus, phytoplankton communities and densities were undergoing a directional change responding to environmental fluctuations during the mixing period in Jinpen reservoir.
With respect to vertical variation within the phytoplankton community composition, Bacillariophyta and Chlorophyta were most dominant at depths of 5, 10, and 0, 2 m. Due to the lack of a self-regulating mechanism, Chlorophyta prefer living at the optimal light intensity in the reservoir, while diatoms are mainly distributed in the lower layer as a result of disliking light [55]. Hence, the vertical distribution of phytoplankton is closely correlated with their living habits. The Bacillariophyta and Chlorophyta communities were primarily composed of Melosira sp., Cyclotella sp., and Chlorella sp., respectively, and all other algal divisions accounted for less than 5% during the study period (Figure 4). However, in Lough Neagh, Bacillariophyta were mainly composed by Fragilaria construens, Amphora ovalis, and others during the winter periods, which is inconsistent with our results due to the differences of geographical distribution and water quality [56]. Seasonal variation in the phytoplankton community composition is significant. Recently, Escalas et al. [55] investigated the impact of dominant phytoplankton across four summer campaigns from 50 contrasted environmental waterbodies in the Île-de-France region, and found that Chlorophyta and Cyanobacteria were widespread in summer communities, representing most of the cases of dominance with a relative abundance of 35.5–40.6% and 30.3–36.5%, respectively. However, by investigating the effects of environmental changes on dominant phytoplankton composition, diversity, and community stability in the Římov reservoir (Czech Republic) over 32 years, Znachor et al. [57] found that the phytoplankton of the reservoir underwent a substantial compositional shift towards diatoms, and that the most significant predictors of phytoplankton-composition structure driving the overall phytoplankton assemblage were TN, inflow rate, and surface level.
A co-occurrence network that was developed and has been widely applied to reveal associations among different microbial functional groups or populations can visually show the interaction between microbial communities [13,25]. As shown in Table 2 and Figure 4, the resulting networks consisted of 17 nodes linked by 15 edges in October, 18 nodes linked by 15 edges in November, and 16 nodes linked by 14 edges in December. Overall, all nodes from different sampling times were clearly parsed into six different modules, respectively (Figure 4). Among these modules, modules I and II were the largest co-occurrence modules in the whole networks, and each module revealed different composing characteristics in a different co-occurrence network. For each co-occurrence network, most modules were relatively independent of the others, which were observed to have fewer relationships with other modules. Conversely, nodes were linked more frequently between module I and II in December. The majority of phytoplankton genera (nodes) in the same modules exhibited close positive correlations, while only a few showed negative. In total, 35 had positive correlations and 9 had negative among these significant correlations. Interspecific interactions exist widely in phytoplankton communities. Among these correlations, positive correlations in particular may indicate that microbes thrive and perform similar ecological functions under the same environment [58]. The obvious differences of the three different networks indicated that correlations between phytoplankton could adjust with time. The main reasons for the seasonal variation of the networks are that the phytoplankton community composition varies over time, and that the dynamics of the environmental parameters can affect the relationships between algal species. For examples, species that share similar ecological niches probably exhibit competitive or mutually exclusive roles in resource shortages, but may reveal positive interactions in resource-rich environments [13]. Some studies proposed that the community structure of phytoplankton is potentially regulated conjointly by both extrinsic and intrinsic drivers deriving from the surrounding environment and interspecific interactions between two species. Recently, Yang et al. [59] explored the mechanisms of phytoplankton community dynamics and demonstrated that interspecific interactions were more important than extrinsic factors in shaping phytoplankton community structure, with a maximum value of 48% in random-forest-based models. Furthermore, they found the key stone species of Oscillatoria granulate were positively correlated with Euglena pisciformis, Trachelomonas and Phacus, but negatively correlated with Chlorella.

3.4. Microbial Eukaryotic Community Composition

Rarefaction curves of OTU number sat a 97% similarity box plot revealed that many reads were sampled (Figure 5). On the basis of the Illumina Miseq high-throughput sequencing, a total of 503,848 high-quality sequences were yielded from 12 samples, with an average length of 134 bp (triplicates) after quality trimming (Table 3 and Figure 5). Miseq sequencing obtained a total of 3610 operational taxonomic units (OTUs) with 97% similarity. The coverage of eukaryotic microbial communities from 12 samples reached 99%, indicating that the sequencing depth of eukaryotes had reached a high level with high data reliability. There were 203,565,159,021,141, and 262 read numbers from October to December, and the average read numbers for each month were 50,891,39,755, and 35,316, respectively. The highest OTUs were observed at a depth of 5 m in November, and the lowest were at a depth of 0 m in December (415 and 132, respectively). Similarly, on the basis of the Illumina Miseq sequencing technique, Zhang et al. [41] investigated the vertical distribution of bacterial communities during thermal stratification in a drinking water reservoir and observed that the highest and lowest OTUs were found in 5 and 0.5 m, respectively.
To explore the dynamic of water eukaryotic community structures, we utilized species-richness estimates (Chao 1) and diversity indices (Shannon, Simpson). Diversity estimators Chao 1, Shannon, and Simpson clearly varied between months and depths (Table 3). Overall, Chao 1 indices were lowest in October and highest from November to December, ranging from 143 to 476. Moreover, the Chao 1 richness index was higher than the OTU value, indicating that more unidentified sequences are still to be studied in the reservoir [50]. Variations of Shannon diversity (H’) were analogous to the Chao 1 index. As shown in Table 3, the highest Shannon diversity indices (H’) and the lowest Simpson diversity indexes (D) were both observed in November at a depth of 2 m, with values of 4.3 and 0.026, respectively, indicating the highest diversity and most even taxon distribution of eukaryotic communities during this period. Spatial variations of eukaryotic community diversity generally also exhibited surface water that was lower than the other layers. Previously, Zhang et al. [41] reported that the water microbial community Chao 1 index ranged from 407 to 506, changing with depth. Additionally, Shannon diversity (H’) was lower in the surface layer and higher in the bottom. By investigating the impacts of thermal stratification on the bacterial community with the methods of BIOLOG and 454 pyrosequencing, Yang et al. [50] also found that bacterial community diversity in the hypolimnion was higher than that in the epilimnion, indicating that stratification significantly affected microbial communities.
The eukaryotic community at the phylum level distinctly varied during the studied periods in Jinpen reservoir. Figure 6 shows an overview of the eukaryotic community composition. The 3610 operational taxonomic units (OTUs) with 97% similarity were mainly classified into 20 phyla: Arthropoda, Ochrophyta, Ciliophora, Cryptomonadales, Nematoda, Cercozoa, Phragmoplastophyta, Cryptomycota, Vertebrata, Kathablepharidae, Chytridiomycota, Choanoflagellida, Euglenozoa, Dinoflagellata, Ascomycota, Basidiomycota, Blastocladiomycota, Protalveolata, Unclassified and Others (Figure 6). Different from the variation of algal community composition and abundance, eukaryote biodiversity increased with time, and species distribution was more even (Figure 6). From October to December, Arthropoda (6.63% to 79.19%), Ochrophyta (5.60% to 35.16%), Ciliophora (1.81% to 10.93%), and Cryptomonadales (0.25% to 11.48%) were the most dominant eukaryotes in common, which contributed over 50% of the total eukaryotic community (except for the point of 0 m in November). From those taxa, arthropod abundance declined significantly over time, especially at the point of 0 to 5 m in November, with an average decrease of 55.55% (Figure 6), while algal density was also markedly reduced during this period. Arthropoda, as a kind of zooplankton, grow by feeding on algae, organic particles, and bacteria, and some even have algae as their sole source of nutrition. Thus, to a certain extent, arthropod growth might be affected by the decrease of their food sources [31]. Recently, research observed that a relative abundance of arthropods declined over time in a subtropical reservoir, and that they were one of the major contributors to eukaryotic plankton community rearrangements [15]. In addition, Nematoda peaked in November (at a depth of 0 m), with a maximal relative abundance of 60.87% (Figure 6), becoming the most abundant and unique species during this period, which was rare in other periods. At the same time, the algal-cell concentration sharply decreased, with a maximal reduction of 91 × 104 cells/L. Special variation of Nematoda in November might be associated with the sharp decrease of algal cell concentration at a depth of 0 m, which meant that a marked decrease of algae might boost the active growth of Nematoda. Under the same condition, differences between Arthropod and Nematoda dynamics possibly resulted from two aspects: on the one hand, there are different resource-utilization characteristics and abilities of these two eukaryotes to adapt to particular ecosystems; on the other hand, there are interrelationships between different eukaryotes in terms of predation, competition, and parasitism [11]. Cercozoa dominated from November to December, with a maximal relative abundance of 12.50% in November (at a depth of 2 m). Similarly, as the common dominant taxon, Ciliophora were also observed with a maximal relative abundance (10.93%) in November (at a depth of 2 m). Ciliophora and Cercozoa are widely distributed and ubiquitous in lakes [11,14,32,60,61] and reservoirs [15,62], and they are important consumers of algae and bacteria, which occupy a pivotal position in aquatic food webs by transferring nutrients to higher trophic levels [61]. Nematoda, Ciliophora, and Cercozoa are also affiliated with zooplankton, inducing a decrease of algal density by predation and altering the structure of phytoplankton communities by selecting specific algae [31]. Several studies have reported that variations in zooplankton species distribution and abundance may alter ecosystems via cascading effects, and found that climate variability can be indicated by zooplankton [55,56,63]. Freshwater fungi are ubiquitous and diverse, occupying an important position as saprotrophs or parasites in the function of freshwater ecosystems [24]. In this study, five fungus phyla were observed in the eukaryotic community: Cryptomycota, Chytridiomycota, Ascomycota, Basidiomycota, and Blastocladiomycota. Compared with other fungi, the abundance of Cryptomycota significantly increased in November and December, reaching a maximum of 10.70% (at a depth of 5 m) in November, whereas algal density decreased during this period. Cryptomycota, a newly discovered type of fungus, is widely found in water environments and usually parasitizes algae or other fungi, absorbing nutrients through phagocytosis and thereby directly affecting the growth of hosts [64,65]. Hence, results in our study partly confirmed that there might be a saprophytic interaction between Cryptomycota and algal decomposition. Recently, Cryptomycota were also identified by a study conducted by Zhang et al. [66] exploring the relationships between axenic Microcystis aeruginosa and aquatic microbes by using a laboratory co-culture system. They found that abundance of Cryptomycota and other eukaryotes decreased after co-culturing.
To further investigate the eukaryotic community, a heat-map fingerprint with 79 representative eukaryotic species was generated at the genus level, as shown in Figure 7. Consistent with the variation of phylum level, the eukaryotic community composition in Jinpen reservoir exhibited distinct patterns at the genus level, which became more diverse and different over time. Aulacoseira and Thalassiosira were the common genera in the entire sampling time, and contributed to a higher relative abundance for the eukaryotic community with a sum of relative abundance ranging from 17.93% to 49.56%. In October, Tetraselmis and Cryptomonas were the most abundant genera in addition to the common genera. From these two taxa, Tetraselmis reached its maximal relative abundance of 25.82% at a depth of 0 m, while Cryptomonas was absolutely dominant at depths from 2 to 10 m (16.24% to 20.29%). Strombidium was another abundant genus in October appearing at a depth of 0 m witha relative abundance of 8.28%, which was rare atother sampling times. In November, Sargassum contributed to the highest relative proportions in the community at a depth of 0 m, with a relative abundance of 18.86%, and was closely followed by Entodinium (14.20%). Bodo, Rimostrombidium, and Halteria dominated at depths of 2, 5, and 10 m and accounted for 8.92%, 8.55%, and 7.57%, respectively (with the exception of Aulacoseira and Thalassiosira). Consistent with previous observations, Alveolata comprising the specific genera Halteriagrandinella and Rimostrombidiumlacustris was found as the second most important taxon, dominating in the day-4 clone library during a Mesocosm experiment [14]. In December, Teleaulax and Dinobryon were other abundant genera after Aulacoseira and Thalassiosira, with an average relative abundance of 11.33% and 11.17%, respectively. Teleaulax belongs to Cryptophyta, which was also detected by Mikhailov et al. [60] during spring algal blooms in Lake Baikal. It contributed the highest proportions of 15% for total eukaryotic reads, and was positively correlated with Cryptococcus. Mikhailov et al. also observed Halteria ranging from 0.1% to 23%, with a higher betweenness centrality of 0.105 in the co-occurrence network of microbial eukaryotes.

3.5. Relationships between Phytoplankton and Eukaryotic Communities Explored by the Network Approach

Co-occurrence networks provide synthetic and clear representations of the microbial community for a better understanding of its structure and function [25,55]. Here, we assessed the co-occurrence patterns of eukaryotic organisms at the genus level from October to December to explore interactions between different eukaryotes. Co-occurring eukaryotic genus networks showed that community composition significantly variedin the three months. As shown in Table 4 and Figure 8, Figure 9 and Figure 10, the resulting networks consisted of 67 nodes linked by 147 edges in October, 72 nodes linked by 144 edges in November, and 65 nodes linked by 204 edges in December. The observed link numbers increased, reaching their maximum in December, whereas node numbers were similar in the three months. Network density was highest in December, while the October and December networks exhibited lower and similar density. The indices of average clustering coefficient and average path length represent the size of the networks [55], which decreased and increased from October to December, respectively, and reached their minimum and maximum in December (Table 4).
In October, all nodes were clearly parsed into eight different modules (Figure 8). Among these modules, modules I to VI, as the largest co-occurrence modules, dominated in the whole network, containing 12, 11, 11, and 9 nodes (genera), respectively, witheach module exhibiting different composing characteristics. The connected nodes with a larger size (e.g., Pleiochaeta) had higher significant negative or positive associations with other taxa (genera). In module I, Hyphochytrium was observed to have significant negative correlations with Poteriospumella, Phanerochaete, and Korotnevella, while significant positive correlations existed among Pleiochaeta, Poteriospumella, Kazachstania, Phanerochaete, Korotnevella, and Aulacoseira. Modules VI and VII and others were independent, and observed to have fewer relationships with other modules. Eukaryotes from these three modules revealed a significant positive correlation with one another. However, eukaryotes from module IV mostly showed significant a negative correlation with one another, and only three significant positive correlations were observed among Naganishia, Obertrumia, and Paramicrosporidium.
In November, all nodes were clearly parsed into seven different modules, of which modules I–III dominated in the network, accounting for 36.11%, 20.83%, and 13.89% of the whole network, respectively (Figure 9). In module I, notably, Penicillium was observed to have more significant correlations with other eukaryotic genera, with the highest degree (number of connections) of 19. Among these correlations, 18 were negative and only one was positive (with Bromeliophrya), indicating that Penicillium negatively affected the growth of the majority of other genera, or that there might be an antagonistic and competitive [36], or saprophytic [14], association between Penicillium and these genera. Correlations from modules II and V–VII were all positive, and positive correlations from modules V–VII were especially significant, indicating that the eukaryotic genera from these modules might arise from similar niches or mutualism [60,61,66]. Additionally, in module II, Phanerochaete was observed to have the highest positive correlations with other eukaryotic genera (adegree of 20), followed closely by Pleiochaeta with a degree of 13. Naganishia, Dinobryon, Nuclearia, Ascochyta, and Euglena from module III were significantly and positively correlated with each other, which is consistent with the significant association among Grammatophora, Strombidium, and Heterocapsa from module IV.
In December, all nodes were clearly parsed into four different modules, which was the lowest during the studied period (Figure 10). Modules I–III comprised 22, 21, and 18 nodes (genera), respectively, when the module contained only four nodes. In module I, significant negative correlations with Aureobasidium and Paramicrosporidium were observed for Paraphysomonas, while significant positive correlations were observed between Aureobasidium and Paramicrosporidium. Likewise, in module II, Monodus was significantly and positively linked with Tintinnidium and Stauridium, whereas a distinct negative association existed between Ascochyta and Gymnophrys. In module III, Malassezia had the highest degree (18), but was negatively correlated with numerous eukaryotic genera. In contrast, Cryothecomonas was observed to have the highest positive correlations with other eukaryotic genera (a degree of 15). It was followed closely by Chrysoxys with a degree of 12. Module VI was relatively independent while other modules were linked more frequently, and genera from module VI all exhibited significant positive correlations with one another.
To summarize, the characteristics of co-occurring eukaryotic genera networks further showed that community composition strongly varied over time. Eukaryotes exhibited significant and predominantly positive or negative correlations between each other, which were specific to time. Positive correlations might arise from similar niches or mutualism while negative correlations might reflect a non-overlapping niche or antagonism [61]. Analysis of macro-ecological networks shows that co-occurrences may be a better indicator of niche preferences than biotic interactions [11]. Recently, Xue et al. [15] employed co-occurrence network analysis to investigate the relationships of abundant and rare planktonic eukaryotes in a subtropical reservoir, and found that dominant taxa were mainly derived from rare species that may play vital roles in regulating the ecological function in aquatic systems as keystone species. The authors also observed that rare and non-rare taxa exhibited synergistic effects that might contribute to maintaining the stability and resilience of eukaryotic communities and ecological function. Mikhailov et al. [60] found that species from the microbial metacommunity occupied positive correlations and were highly interconnected with one another within a specific domain. Positive associations between bacteria and eukaryotes may be used as evidence of mutualistic interactions, whereas negative associations may indicate competition.

3.6. Effect of the Environmental Factors on Communities

To better understand the effects of water quality on phytoplankton and eukaryote community diversity, multivariable statistics from redundancy analysis (RDA) were employed to investigate the relationships between water-quality parameters and compositions of algal and eukaryotic communities on the basis of microscopy analysis and Illunina MiSeq sequencing. RDA results revealed that distinct phytoplankton and eukaryote communities were detected during the studied periods in the Jinpen drinking water reservoir. As shown in Figure 11, environmental parameters explained 96.6% variance of phytoplankton communities and 96.3% variance of eukaryote communities, respectively. Phytoplankton community structure was significantly influenced by water temperature, NH4+-N, TN, TP, DOC, conductivity, Mn, Fe, DO, and pH. The most important determinants of eukaryotic community structure were TN, NH4+-N, water temperature, DOC, Mn, conductivity, and pH among all the physic-chemical parameters. Moreover, TN, NH4+-N, and water temperature exhibited a reverse correlation relative to DOC, Mn, conductivity, and pH. This suggests that water quality significantly affects the variability of phytoplankton and eukaryote community compositions, particularly TN, NH4+-N, water temperature, DOC, Mn, conductivity, and pH, which were common factors in controlling population variation. However, DO, TP, and Fe were the only influential factors on the structure of the phytoplankton communities. A similar study conducted by Chen et al. [14] indicated that the microbial eukaryotic community composition was mostly correlated with changes of pH, DO, and DOC concentrations, which is generally consistent with our results.
By investigating the variations of phytoplankton community structure over two years (a wet year and a dry year), Qiu et al. [45] observed that nitrogen fluctuations caused by rainfall were more significant than other environmental parameters in the Zhoucun reservoir, which is one of the crucial factors responsible for the dynamics of phytoplankton growth. In addition, they found that water temperature, thermal stratification, and light were considerable driving factors for shaping the phytoplankton community. Recently, Yang et al. [4] found that combined disturbance events such as water-level fluctuations, human activities, climate change, and rainfall exert a complex and diverse impact on phytoplankton communities across space and time compared to single disturbances. They observed that phytoplankton community dynamics in two subtropical reservoirs over the past six years reflected distinct and reversible patterns (from a cyanobacteria-dominated state to non-cyanobacterial-taxa-dominated state) under multiple environmental-disturbance conditions. From those disturbances, water temperature and rainfall strongly affected cyanobacterium biomass, while nutrients and water level, as the main driving factors, had great impacts onthephytoplankton community structure. Moreover, the impact of the disturbances varied in the taxon and organizational levels. However, Bullerjahn et al. [53] recently observed that winter phytoplankton communities in ice fractions were photo-synthetically active. Additionally, the similarities of ice-derived microbial communities in these two totally different lakes suggest that the conditions for ice formation chose common community compositions. Capo et al. [32] demonstrated that modules comprising Dinophyceae and unclassified Alveolata were significantly related to air temperature in two lakes that were exposed to varying degrees of phosphorus enrichment. Compared with our study, these inconsistent results indicate that microbial community compositions in an aquatic system are affected by a variety of environmental factors and are altered by changes in local conditions, leading to totally different community structure variation trajectories.
In conclusion, our results demonstrated that phytoplankton and eukaryotes in the water column were both affected by changing environmental conditions, and that they responded accordingly. At the same time, such extrinsic factors deriving from the surrounding environment are incomplete for explaining the dynamics of phytoplankton and eukaryote communities. Most importantly, intrinsic feedback mechanisms and interspecific interactions play critical roles in a community’s response to variations in external factors [10]. Our results revealed that relationships between phytoplankton and eukaryotes were complicated and important. The mutual interactions existed widely and exerted a strong impact on shaping phytoplankton and eukaryote assemblage. For example, Ciliophora and Cercozoa, which are important consumers of algae and bacteria, were dominant in the eukaryotic community. Phytoplankton blooms can boost the growth rate of microbes [15]. Thus, the increase of these taxa most likely resulted from the increase of their prey or food. Results in our study also exhibited proof that the growth of phytoplankton is closely related to fungus. For example, Cryptomycota had a saprophytic interaction with phytoplankton. Furthermore, Penicillium exhibited higher connections with other eukaryotic genera that were negatively related to others such as Cryptomonas. These negative correlations indicated that there might be antagonistic or saprophytic relationships between Penicillium and these genera. Conversely, our co-occurrence network also showed some negative associations among eukaryotes. Phanerochaete, Pleiochaeta, Cryothecomonas, and Chrysoxys positively correlated with other genera. Our results provide a deeper insight into the dynamics of phytoplankton and eukaryote communities in reservoirs from a new perspective. Given our observations, we recommend that future studies focus on understanding the correlations of phytoplankton and eukaryote communities through long-term and extensive investigations in order to fully understand their intrinsic mechanisms and relative importance in affecting microbial community dynamics.

4. Conclusions

In summary, results of this study proved that variations of phytoplankton and water eukaryotes were closely associated with each other during winter in the Jinpen drinking water reservoir. Significant spatial-temporal changes were revealed in the composition of the eukaryotic and phytoplankton communities. The co-occurrence of phytoplankton indicated that the community structure varied remarkably over time. Moreover, Bacillariophyta and Chlorophyta were the most abundant taxa, with a total relative abundance of more than 97% throughout the studied periods, which were primarily composed of Melosira sp., Cyclotella sp., and Chlorella sp., respectively. Illumina Miseq sequencing suggested that eukaryote biodiversity gradually increased and species distribution became more even over time. Arthropoda (6.63% to 79.19%), Ochrophyta (5.60% to 35.16%), Ciliophora (1.81% to 10.93%), and Cryptomonadales (0.25% to 11.48%) were the most dominant eukaryotes in common, contributing over 50% of the total eukaryotic community. The abundance of arthropod was significantly reduced at the point of 0 to 5 m from October to November, with an average decrease of 55.55%, while algal density significantly declined at the same time, indicating a negative association between arthropod and phytoplankton. Nematoda and Cercozoa peaked in November with maximal relative abundances of 60.87% and 12.50%, respectively. Cryptomycota, a newly discovered type of fungus, was special among the most dominant eukaryotes; when the number of algae decreased, the abundance of Cryptomycota increased, suggesting that there might be a saprophytic interaction between Cryptomycota and phytoplankton decomposition. The co-occurrence network analysis showed that there existed significant correlations between eukaryotic and phytoplankton community composition, suggesting that eukaryotes respond to variations in the dynamic structure of phytoplankton community. For instance, Penicillium was observed in the majority of significant negative correlations with other eukaryotic genera, which indicates that Penicillium negatively affects the growth of other genera and that there might be an antagonistic and competitive (or saprophytic) association between them. In addition, RDA results demonstrated that total nitrogen, ammonia nitrogen, water temperature, DOC, Mn, conductivity, and pH contributed the most to variations of phytoplankton and eukaryotic microorganism communities. Furthermore, the DO, TP, and Fe factors can also significantly alter the structures of phytoplankton communities. These conclusions confirm that intricate and vital mutual interactions exist between phytoplankton and eukaryotes in freshwater ecosystems. On this basis, more studies are needed to further investigate the ecological functions of water microbes in drinking water reservoirs.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1660-4601/17/4/1128/s1. Figure S1: Images of phytoplankton under microscopy.

Author Contributions

S.C. and T.H. conceived the initial idea for the present study; M.Y. determined all the environmental parameters, analyzed the data and wrote the first draft of the manuscript; N.L. wrote the manuscript with critical revision among all co-authors; B.L. constructed networks. K.L., R.Z. and Y.M. partly give suggestions on the early versions of the manuscript. All authors have read and approved the final manuscript.

Funding

This research was jointly supported by the National Natural Science Foundation of China (Nos. 51978561; 51979217), the International Science and TechnologyCooperation Program in Shaanxi Province (No. 2018KW-011), the Natural Science Basic Research Plan in Shaanxi Province of China (No. 2018JM5125), and the Shaanxi Provincial Key Research and Development Project (No. 2018ZDXM-SF-020). The authors would like to thank the editor and anonymous reviewers for their constructive comments that improved the earlier version of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, H.H.; Jia, J.Y.; Chen, S.N.; Huang, T.L.; Wang, Y.; Zhao, Z.F.; Feng, J.; Hao, H.H.; Li, S.L.; Ma, X. Dynamics of bacterial and fungal communities during the outbreak and decline of an algal bloom in a drinking water reservoir. Int. J. Environ. Res. Public Health 2018, 15, 361. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Huang, T.L.; Li, X.; Rijnaarts, H.; Grotenhuis, T.; Ma, W.X.; Sun, X.; Xu, J. Effects of storm runoff on the thermal regime and water quality of a deep, stratified reservoir in a temperate monsoon zone, in Northwest China. Sci. Total Environ. 2014, 485–486, 820–827. [Google Scholar] [CrossRef] [PubMed]
  3. Zhang, H.; Zhou, X.; Wang, L.B.; Wang, W.D.; Xu, J.L. Concentrations and potential health risks of strontium in drinking water from Xi’an, Northwest China. Ecotoxicol. Environ. Saf. 2018, 164, 181–188. [Google Scholar] [CrossRef] [PubMed]
  4. Yang, J.R.; Lv, H.; Isabwe, A.; Liu, L.M.; Yu, X.Q.; Chen, H.H.; Yang, J. Disturbance-induced phytoplankton regime shifts and recovery of cyanobacteria dominance in two subtropical reservoirs. Water Res. 2017, 120, 52–63. [Google Scholar] [CrossRef] [PubMed]
  5. Xiao, L.J.; Wang, T.; Hu, R.; Han, B.P.; Wang, S.; Qian, X.; Padisák, J. Succession of phytoplankton functional groups regulated by monsoonal hydrology in a large canyon-shaped reservoir. Water Res. 2011, 45, 5099–5109. [Google Scholar] [CrossRef]
  6. Wang, S.; Qian, X.; Han, B.P.; Luo, L.C.; Hamilton, D.P. Effects of local climate and hydrological conditions on the thermal regime of a reservoir at Tropic of Cancer, in southern China. Water Res. 2012, 46, 2591–2604. [Google Scholar] [CrossRef]
  7. Ma, W.X.; Huang, T.L.; Li, X.; Zhang, H.H.; Ju, T. Impact of short-term climate variation and hydrology change on thermal structure and water quality of a canyon-shaped, stratified reservoir. Environ. Sci. Pollut. Res. 2015, 22, 18372–18380. [Google Scholar] [CrossRef]
  8. Wilhelm, S.W.; Farnsley, S.E.; LeCleir, G.R.; Layton, A.C.; Satchwell, M.F.; DeBruyn, J.M.; Boyer, G.L.; Zhu, G.; Paerl, H.W. The relationships between nutrients, cyanobacterial toxins and the microbial community in Taihu (Lake Tai), China. Harmful Algae 2011, 10, 207–215. [Google Scholar] [CrossRef]
  9. Li, H.; Xing, P.; Chen, M.; Bian, Y.; Wu, Q.L. Short-term bacterial community composition dynamics in response to accumulation and breakdown of Microcystis blooms. Water Res. 2011, 45, 1702–1710. [Google Scholar] [CrossRef]
  10. Van Wichelen, J.; Vanormelingen, P.; Codd, G.A.; Vyverman, W. The common bloom-forming cyanobacterium Microcystis is prone to a wide array of microbial antagonists. Harmful Algae 2016, 55, 97–111. [Google Scholar] [CrossRef]
  11. Shi, X.L.; Li, S.N.; Zhang, M.; Liu, C.Q.; Wu, Q.L. Temperature mainly determines the temporal succession of the photosynthetic picoeukaryote community in Lake Chaohu, a highly eutrophic shallow lake. Sci. Total Environ. 2020, 702, 134803. [Google Scholar] [CrossRef]
  12. Isabwe, A.; Yang, J.R.; Wang, Y.; Liu, L.; Chen, H.; Yang, J. Community assembly processes underlying phytoplankton and bacterioplankton across a hydrologic change in a human-impacted river. Sci. Total Environ. 2018, 630, 658–667. [Google Scholar] [CrossRef] [PubMed]
  13. Zhao, D.Y.; Shen, F.; Zeng, J.; Huang, R.; Yu, Z.; Wu, Q.L. Network analysis reveals seasonal variation of co-occurrence correlations between Cyanobacteria and other bacterioplankton. Sci. Total Environ. 2016, 573, 817–825. [Google Scholar] [CrossRef]
  14. Chen, M.J.; Chen, F.Z.; Xing, P.; Li, H.B.; Wu, Q.L. Microbial eukaryotic community in response to Microcystis spp. bloom, as assessed by an enclosure experiment in Lake Taihu, China. FEMS Microbiol. Ecol. 2010, 74, 19–31. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Xue, Y.Y.; Chen, H.H.; Yang, J.R.; Liu, M.; Huang, B.Q.; Yang, J. Distinct patterns and processes of abundant and rare eukaryotic plankton communities following a reservoir cyanobacterial bloom. ISME J. 2018, 12, 2263–2277. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Su, X.M.; Steinman, A.D.; Xue, Q.J.; Zhao, Y.Y.; Tang, X.M.; Xie, L.Q. Temporal patterns of phyto- and bacterioplankton and their relationships with environmental factors in Lake Taihu, China. Chemosphere 2017, 184, 299–308. [Google Scholar] [CrossRef] [PubMed]
  17. Naeem, S.; Hahn, D.R.; Schuurman, G. Producer-decomposer co-dependency influences biodiversity effects. Nature 2000, 403, 762–764. [Google Scholar] [CrossRef]
  18. Niu, Y.; Shen, H.; Chen, J.; Xie, P.; Yang, X.; Tao, M.; Ma, Z.; Qi, M. Phytoplankton community succession shaping bacterioplankton community composition in Lake Taihu, China. Water Res. 2011, 45, 4169–4182. [Google Scholar] [CrossRef]
  19. Su, X.M.; Steinman, A.D.; Tang, X.M.; Xue, Q.J.; Zhao, Y.Y.; Xie, L.Q. Response of bacterial communities to cyanobacterial harmful algal blooms in Lake Taihu, China. Harmful Algae 2017, 68, 168–177. [Google Scholar] [CrossRef]
  20. Berry, M.A.; Davis, T.W.; Cory, R.M.; Duhaime, M.B.; Johengen, T.H.; Kling, G.W.; Marino, J.A.; Den Uyl, P.A.; Gossiaux, D.; Dick, G.J.; et al. Cyanobacterial harmful algal blooms are a biological disturbance to Western Lake Erie bacterial communities. Environ. Microbiol. 2017, 19, 1149–1162. [Google Scholar] [CrossRef]
  21. Chen, S.N.; He, H.Y.; Zong, R.R.; Liu, K.W.; Yan, M.M.; Xu, L.; Miao, Y.T.; Huang, X. Geographical patterns of algal communities associated with different urban lakes in China. Int. J. Environ. Res. Public Health 2020, 17, 1009. [Google Scholar] [CrossRef] [Green Version]
  22. Scherer, P.I.; Millard, A.D.; Miller, A.; Schoen, R.; Raeder, U.; Geist, J.; Zwirglmaier, K. Temporal dynamics of the microbial community composition with a focus on toxic cyanobacteria and toxin presence during harmful algal blooms in two South German lakes. Front. Microbiol. 2017, 8, 2387. [Google Scholar] [CrossRef] [PubMed]
  23. Liu, M.; Xiao, T.; Sun, J.; Wei, H.; Wu, Y.; Zhao, Y.; Zhang, W.C. Bacterial community structures associated with a natural spring phytoplankton bloom in the Yellow Sea, China. Deep Sea Res. Part II Top. Stud. Oceanogr. 2013, 97, 85–92. [Google Scholar] [CrossRef]
  24. Sun, J.Y.; Song, Y.; Ma, Z.P.; Zhang, H.J.; Yang, Z.D.; Cai, Z.H.; Zhou, J. Fungal community dynamics during a marine dinoflagellate (Noctiluca scintillans) bloom. Mar. Environ. Res. 2017, 131, 183–194. [Google Scholar] [CrossRef]
  25. Tan, S.; Zhou, J.; Zhu, X.; Yu, S.; Zhan, W.; Wang, B.; Cai, Z. An association network analysis among microeukaryotes and bacterioplankton reveals algal bloom dynamics. J. Phycol. 2015, 51, 120–132. [Google Scholar] [CrossRef]
  26. Teeling, H.; Fuchs, B.M.; Becher, D.; Klockow, C.; Gardebrecht, A.; Bennke, C.M.; Kassabgy, M.; Huang, S.; Mann, A.J.; Waldmann, J.; et al. Substrate-controlled succession of marine bacterioplankton populations induced by a phytoplankton bloom. Science 2012, 336, 608–611. [Google Scholar] [CrossRef]
  27. Guo, Y.; Liu, M.; Liu, L.; Liu, X.; Chen, H.; Yang, J. The antibiotic resistome of free-living and particle-attached bacteria under a reservoir cyanobacterial bloom. Environ. Int. 2018, 117, 107–115. [Google Scholar] [CrossRef] [PubMed]
  28. Xiao, X.; Pei, M.; Liu, X.; Zhao, Y.; Liang, Y. Planktonic algal bloom significantly alters sediment bacterial community structure. J. Soils Sediments 2017, 17, 2547–2556. [Google Scholar] [CrossRef]
  29. Jia, Y.; Han, G.; Wang, C.; Guo, P.; Jiang, W.; Li, X.; Tian, X. The efficacy and mechanisms of fungal suppression of freshwater harmful algal bloom species. J. Hazard. Mater. 2010, 183, 176–181. [Google Scholar] [CrossRef]
  30. Rooney-Varga, J.N.; Giewat, M.W.; Savin, M.C.; Sood, S.; LeGresley, M.; Martin, J.L. Links between Phytoplankton and Bacterial Community Dynamics in a Coastal Marine Environment. Microb. Ecol. 2005, 49, 163–175. [Google Scholar] [CrossRef]
  31. Li, C.C.; Feng, W.Y.; Chen, H.Y.; Li, X.F.; Song, F.H.; Guo, W.J.; Giesy, J.P.; Sun, F.H. Temporal variation in zooplankton and phytoplankton community species composition and the affecting factors in Lake Taihu-a large freshwater lake in China. Environ. Pollut. 2019, 245, 1050–1057. [Google Scholar] [CrossRef]
  32. Capo, E.; Debroas, D.; Arnaud, F.; Perga, M.E.; Chardon, C.; Domaizon, I. Tracking a century of changes in microbial eukaryotic diversity in lakes driven by nutrient enrichment and climate warming. Environ. Microbiol. 2017, 19, 2873–2892. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Faust, K.; Raes, J. Microbial interactions: From networks to models. Nat. Rev. Microbiol. 2012, 10, 538–550. [Google Scholar] [CrossRef] [PubMed]
  34. Larsen, P.E.; Field, D.; Gilbert, J.A. Predicting bacterial community assemblages using an artificial neural network approach. Nat. Methods 2012, 9, 621–625. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Sun, W.; Xiao, E.; Häggblom, M.; Krumins, V.; Dong, Y.; Sun, X.; Li, F.; Wang, Q.; Li, B.; Yan, B. Bacterial survival strategies in an alkaline tailing site and the physiological mechanisms of dominant phylotypes as revealed by metagenomic analyses. Environ. Sci. Technol. 2018, 52, 13370–13380. [Google Scholar] [CrossRef]
  36. Sun, W.M.; Xiao, E.Z.; Krumins, V.; Häggblom, M.M.; Dong, Y.R.; Pu, Z.L.; Li, B.Q.; Wang, Q.; Xiao, T.F.; Li, F.B. Rhizosphere microbial response to multiple metal(loid)s in different contaminated arable soils indicates crop-specific metal-microbe interactions. Appl. Environ. Microbiol. 2018. [Google Scholar] [CrossRef] [Green Version]
  37. Zhang, H.H.; Feng, J.; Chen, S.N.; Zhao, Z.F.; Li, B.Q.; Wang, Y.; Jia, J.Y.; Li, S.L.; Wang, Y.; Yan, M.M.; et al. Geographical patterns of nirS gene abundance and nirS-Type denitrifying bacterial community associated with activated sludge from different wastewater treatment plants. Microb. Ecol. 2019, 77, 304–316. [Google Scholar] [CrossRef] [PubMed]
  38. Yin, Y.; Gu, J.; Wang, X.; Zhang, Y.; Zheng, W.; Chen, R.; Wang, X. Effects of rhamnolipid and Tween-80 on cellulase activities and metabolic functions of the bacterial community during chicken manure composting. Bioresour. Technol. 2019, 288, 121507. [Google Scholar] [CrossRef]
  39. Zhang, H.H.; Zhao, Z.F.; Chen, S.N.; Wang, Y.; Feng, J.; Jia, J.Y.; Kang, P.L.; Li, S.L. Geographical patterns of denitrifying bacterial communities associated with different urban lakes. RSC Adv. 2018, 8, 17079–17090. [Google Scholar] [CrossRef] [Green Version]
  40. Wang, Y.; Zhang, R.; Zheng, Q.; Deng, Y.; Van Nostrand, J.D.; Zhou, J.; Jiao, N. Bacterioplankton community resilience to ocean acidification: Evidence from microbial network analysis. ICES J. Mar. Sci. 2015, 73, 865–875. [Google Scholar] [CrossRef]
  41. Zhang, H.H.; Chen, S.N.; Huang, T.L.; Ma, W.X.; Xu, J.L.; Sun, X. Vertical distribution of bacterial community diversity and water quality during the Reservoir thermal stratification. Int. J. Environ. Res. Public Health 2015, 12, 6933–6945. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Greenberg, A.E.; Clesceri, L.S.; Eaton, A.D. Standard Methods for the Examination of Water and Wastewater, 18th ed.; American Public Health Association: Washington, DC, USA, 1992. [Google Scholar]
  43. Hu, H.J.; Wei, Y.X. The Freshwater Algae of China; Science and Technology Press: Beijing, China, 2006. [Google Scholar]
  44. Ministry of Environmental Protection of China (MEPC). Monitoring and Analysis Methods of Water and Wastewater, 4th ed.; China Environmental Science Press: Beijing, China, 2002.
  45. Qiu, X.P.; Huang, T.L.; Zeng, M.Z. Differences in phytoplankton dynamics and community structure between a wet year and dry year in the Zhoucun Reservoir. J. Freshw. Ecol. 2016, 31, 377–391. [Google Scholar] [CrossRef] [Green Version]
  46. Parulekar, N.N.; Kolekar, P.; Jenkins, A.; Kleiven, S.; Utkilen, H.; Johansen, A.; Sawant, S.; Kulkarni-Kale, U.; Kale, M.; Sæbø, M. Characterization of bacterial community associated with phytoplankton bloom in a eutrophic lake in South Norway using 16S rRNA gene amplicon sequence analysis. PLoS ONE 2017, 12, e0173408. [Google Scholar] [CrossRef] [PubMed]
  47. Liu, L.M.; Liu, M.; Wilkinson, D.M.; Chen, H.H.; Yu, X.Q.; Yang, J. DNA metabarcoding reveals that 200-μm-size-fractionated filtering is unable to discriminate between planktonic microbial and large eukaryotes. Mol. Ecol. Resour. 2017, 17, 991–1002. [Google Scholar] [CrossRef]
  48. Landa, M.; Blain, S.; Christaki, U.; Monchy, S.; Obernosterer, I. Shifts in bacterial community composition associated with increased carbon cycling in a mosaic of phytoplankton blooms. ISME J. 2015, 10, 39–50. [Google Scholar] [CrossRef]
  49. Su, X.M.; Steinman, A.D.; Oudsema, M.; Hassett, M.; Xie, L.Q. The influence of nutrients limitation on phytoplankton growth and microcystins production in Spring Lake, USA. Chemosphere 2019, 234, 34–42. [Google Scholar] [CrossRef]
  50. Yang, X.; Huang, T.L.; Zhang, H.H. Effects of seasonal thermal stratification on the functional diversity and composition of the microbial community in a drinking water reservoir. Water 2015, 7, 5525–5546. [Google Scholar] [CrossRef]
  51. Zhou, Z.Z.; Huang, T.L.; Gong, W.J.; Li, Y.; Liu, Y.; Zhou, S.L. Field research on mixing aeration in a drinking water reservoir: Performance and microbial community structure. Int. J. Environ. Res. Public Health 2019, 16, 4221. [Google Scholar] [CrossRef] [Green Version]
  52. Yu, Z.; Yang, J.; Amalfitano, S.; Yu, X.Q.; Liu, L.M. Effects of water stratification and mixing on microbial community structure in a subtropical deep reservoir. Sci. Rep. 2014, 4, 5821. [Google Scholar] [CrossRef] [Green Version]
  53. Bullerjahn, G.S.; McKay, R.M.L.; Bernát, G.; Prášil, O.; Vörös, L.; Pálffy, K.; Tugyi, N.; Somogyi, B. Community dynamics and function of algae and bacteria during winter in central European great lakes. J. Great Lakes Res. 2019. [Google Scholar] [CrossRef]
  54. Lei, Y.; Du, X.; Wang, Y.; Chen, Q.; Tang, H.; Jiang, S. Diatom succession dynamics controlled by multiple forces in a subtropical reservoir in southern China. Quat. Int. 2018, 493, 227–244. [Google Scholar] [CrossRef]
  55. Escalas, A.; Catherine, A.; Maloufi, S.; Cellamare, M.; Hamlaoui, S.; Yéprémian, C.; Louvard, C.; Troussellier, M.; Bernard, C. Drivers and ecological consequences of dominance in periurban phytoplankton communities using networks approaches. Water Res. 2019, 163, 114893. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Dokulil, M.T.; Herzig, A.; Somogyi, B.; Vörös, L.; Donabaum, K.; May, L.; Nõges, T. Winter conditions in six European shallow lakes: A comparative synopsis. Est. J. Ecol. 2014, 63, 111–129. [Google Scholar] [CrossRef] [Green Version]
  57. Znachor, P.; Nedoma, J.; Hejzlar, J.; Seďa, J.; Komárková, J.; Kolář, V.; Mrkvička, T.; Boukal, D.S. Changing environmental conditions underpin long-term patterns of phytoplankton in a freshwater reservoir. Sci. Total Environ. 2019, 710, 135626. [Google Scholar] [CrossRef]
  58. Lima-Mendez, G.; Faust, K.; Henry, N.; Decelle, J.; Colin, S.; Carcillo, F.; Chaffron, S.; Ignacio-Espinosa, J.C.; Roux, S.; Vincent, F.; et al. Determinants of community structure in the global plankton interactome. Science 2015, 348, 1262073. [Google Scholar] [CrossRef] [Green Version]
  59. Yang, W.; Zheng, Z.M.; Zheng, C.; Lu, K.H.; Ding, D.W.; Zhu, J.Y. Temporal variations in a phytoplankton community in a subtropical reservoir: An interplay of extrinsic and intrinsic community effects. Sci. Total Environ. 2018, 612, 720–727. [Google Scholar] [CrossRef]
  60. Mikhailov, I.S.; Zakharova, Y.R.; Bukin, Y.S.; Galachyants, Y.P.; Petrova, D.P.; Sakirko, M.V.; Likhoshway, Y.V. Co-occurrence networks among bacteria and microbial eukaryotes of lake Baikal during a spring phytoplankton bloom. Microb. Ecol. 2019, 77, 96–109. [Google Scholar] [CrossRef]
  61. Hou, F.R.; Zhang, H.J.; Xie, W.J.; Zhou, X.Y.; Zhu, X.Y.; Zhang, D.M. Co-occurrence patterns and assembly processes of microeukaryotic communities in an early-spring diatom bloom. Sci. Total Environ. 2019, 134624. [Google Scholar] [CrossRef]
  62. Liu, L.M.; Chen, H.H.; Liu, M.; Yang, J.R.; Xiao, P.; Wilkinson, D.M.; Yang, J. Response of the eukaryotic plankton community to the cyanobacterial biomass cycle over 6 years in two subtropical reservoirs. ISME J. 2019. [Google Scholar] [CrossRef] [Green Version]
  63. Kwasniewski, S.; Gluchowska, M.; Walkusz, W.; Karnovsky, N.J.; Jakubas, D.; Wojczulanis-Jakubas, K.; Harding, A.M.A.; Goszczko, I.; Cisek, M.; Beszczynska-Möller, A.; et al. Interannual changes in zooplankton on the West Spitsbergen Shelf in relation to hydrography and their consequences for the diet of planktivorous seabirds. ICES J. Mar. Sci. 2012, 69, 890–901. [Google Scholar] [CrossRef] [Green Version]
  64. James, T.Y.; Pelin, A.; Bonen, L.; Ahrendt, S.; Sain, D.; Corradi, N.; Stajich, J.E. Shared signatures of parasitism and phylogenomics unite Cryptomycota and Microsporidia. Curr. Biol. Cb 2013, 23, 1548–1553. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Jones, M.D.; Forn, I.; Gadelha, C.; Egan, M.J.; Bass, D.; Massana, R.; Richards, T.A. Discovery of novel intermediate forms redefines the fungal tree of life. Nature 2011, 474, 200–203. [Google Scholar] [CrossRef] [PubMed]
  66. Zhang, M.; Lu, T.; Paerl, H.W.; Chen, Y.L.; Zhang, Z.Y.; Zhou, Z.G.; Qian, H.F. Feedback regulation between aquatic microorganisms and the bloom-forming cyanobacterium Microcystis aeruginosa. Appl. Environ. Microbiol. 2019. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Water level, average daily inflow and outflow, rainfall, sampling times and depth in Jinpen reservoir from October to December 2018.
Figure 1. Water level, average daily inflow and outflow, rainfall, sampling times and depth in Jinpen reservoir from October to December 2018.
Ijerph 17 01128 g001
Figure 2. Algal cell (A) and chlorophyll a concentrations (B) in Jinpen reservoir from October to December, 2018, the data for each month are the means of all stations and the error bars represent the standard deviations (n = 3). * p < 0.05, ** p < 0.01.
Figure 2. Algal cell (A) and chlorophyll a concentrations (B) in Jinpen reservoir from October to December, 2018, the data for each month are the means of all stations and the error bars represent the standard deviations (n = 3). * p < 0.05, ** p < 0.01.
Ijerph 17 01128 g002
Figure 3. Phytoplankton community’s temporal-spatial variation at a phylum level from October to December in 2018.
Figure 3. Phytoplankton community’s temporal-spatial variation at a phylum level from October to December in 2018.
Ijerph 17 01128 g003
Figure 4. Network of co-occurring phytoplankton communities at the genus level based on a Spearman’s correlation significant analysis (p < 0.01) from October to December in 2018. (AC) represents the October, November, and December sampling times, respectively. The nodes were colored according to their modularity class. Each node represents one individual algal species. The larger the node is, the more connections it has to other nodes. The green lines signify positive correlations, while the red lines signify negative correlations.
Figure 4. Network of co-occurring phytoplankton communities at the genus level based on a Spearman’s correlation significant analysis (p < 0.01) from October to December in 2018. (AC) represents the October, November, and December sampling times, respectively. The nodes were colored according to their modularity class. Each node represents one individual algal species. The larger the node is, the more connections it has to other nodes. The green lines signify positive correlations, while the red lines signify negative correlations.
Ijerph 17 01128 g004aIjerph 17 01128 g004b
Figure 5. Rarefaction curves of the operational taxonomic unit (OTU) numbers at 97% similarity for temporal-spatial distribution of the eukaryotic community in Jinpen reservoir from October to December 2018.
Figure 5. Rarefaction curves of the operational taxonomic unit (OTU) numbers at 97% similarity for temporal-spatial distribution of the eukaryotic community in Jinpen reservoir from October to December 2018.
Ijerph 17 01128 g005
Figure 6. Eukaryotic microorganism community compositions at phylum levels (top 18) in samples from October to December in 2018.
Figure 6. Eukaryotic microorganism community compositions at phylum levels (top 18) in samples from October to December in 2018.
Ijerph 17 01128 g006
Figure 7. A color-scale heat map showing 79 predominant 18S rRNA gene-based eukaryotic sequences classified at the genus level. Red colors indicate higher abundance; blue colors indicate lower abundance.
Figure 7. A color-scale heat map showing 79 predominant 18S rRNA gene-based eukaryotic sequences classified at the genus level. Red colors indicate higher abundance; blue colors indicate lower abundance.
Ijerph 17 01128 g007
Figure 8. Coexistence and interaction of eukaryote community compositions at the genus level based on the 18S rRNA sequence and associated with a water sample in October 2018. The nodes were colored according to modularity class. Each node represents one individual eukaryotic species. The larger the node is, the more connections it has to other nodes. The green lines signify positive correlations, while the red lines signify negative correlations.
Figure 8. Coexistence and interaction of eukaryote community compositions at the genus level based on the 18S rRNA sequence and associated with a water sample in October 2018. The nodes were colored according to modularity class. Each node represents one individual eukaryotic species. The larger the node is, the more connections it has to other nodes. The green lines signify positive correlations, while the red lines signify negative correlations.
Ijerph 17 01128 g008
Figure 9. Coexistence and interaction of eukaryote community compositions at the genus level based on the 18S rRNA sequence and associated with a water sample in November 2018. The nodes were colored according to modularity class. Each node represents one individual eukaryotic species. The larger the node is, the more connections it has to other nodes. The green lines signify positive correlations, while the red lines signify negative correlations.
Figure 9. Coexistence and interaction of eukaryote community compositions at the genus level based on the 18S rRNA sequence and associated with a water sample in November 2018. The nodes were colored according to modularity class. Each node represents one individual eukaryotic species. The larger the node is, the more connections it has to other nodes. The green lines signify positive correlations, while the red lines signify negative correlations.
Ijerph 17 01128 g009
Figure 10. Coexistence and interaction of eukaryote community compositions at the genus level based on the 18S sRNA sequence and associated with water sample in December 2018. The nodes were colored according to modularity class. Each node represents one individual eukaryotic species. The larger the node is, the more connections it has to other nodes. The green lines signify positive correlations, while the red lines signify negative correlations.
Figure 10. Coexistence and interaction of eukaryote community compositions at the genus level based on the 18S sRNA sequence and associated with water sample in December 2018. The nodes were colored according to modularity class. Each node represents one individual eukaryotic species. The larger the node is, the more connections it has to other nodes. The green lines signify positive correlations, while the red lines signify negative correlations.
Ijerph 17 01128 g010
Figure 11. Redundancy analysis (RDA) of algal (A) and eukaryotic (B) communities in Jinpen reservoir from October to December 2018. Stars represent sampling points and 1,2, 3, 4, 5, 6, 7, 8, 9, 10, 11, and 12 represent Oct.—0 m, Oct.—2 m, Oct.—5 m, Oct.—10 m, Nov.—0 m, Nov.—2 m, Nov.—5 m, Nov.—10 m, Dec.—0 m, Dec.—2 m, Dec.—5 m, and Dec.—10 m, respectively. For the algal community, RDA1 explains 86.9%, and RDA2 explains 9.7% of the total variance. For the eukaryotic community, RDA1 explains 66.8%, and RDA 2 explains 29.5% of the total variance. The primary factors for the variables of the water quality data are represented by arrows. DO = dissolved oxygen; TN = total nitrogen; TP = total phosphorus; DOC = dissolved organic carbon.
Figure 11. Redundancy analysis (RDA) of algal (A) and eukaryotic (B) communities in Jinpen reservoir from October to December 2018. Stars represent sampling points and 1,2, 3, 4, 5, 6, 7, 8, 9, 10, 11, and 12 represent Oct.—0 m, Oct.—2 m, Oct.—5 m, Oct.—10 m, Nov.—0 m, Nov.—2 m, Nov.—5 m, Nov.—10 m, Dec.—0 m, Dec.—2 m, Dec.—5 m, and Dec.—10 m, respectively. For the algal community, RDA1 explains 86.9%, and RDA2 explains 9.7% of the total variance. For the eukaryotic community, RDA1 explains 66.8%, and RDA 2 explains 29.5% of the total variance. The primary factors for the variables of the water quality data are represented by arrows. DO = dissolved oxygen; TN = total nitrogen; TP = total phosphorus; DOC = dissolved organic carbon.
Ijerph 17 01128 g011
Table 1. Water quality parameters of Jinpen reservoir measured from October to December in 2018.
Table 1. Water quality parameters of Jinpen reservoir measured from October to December in 2018.
MonthDepth (m)T (°C)DO (mg/L)Turbidity (NTU)Conductivity (S/cm)pHDOC (mg/L)TN (mg/L)NO3N (mg/L)NH4+-N (mg/L)TP (mg/L)Fe (mg/L)Mn (mg/L)
October017.11 ± 0.108.57 ± 0.100.23 ± 0.06142.67 ± 2.527.69 ± 0.022.09 ± 0.091.33 ± 0.040.90 ± 0.020.09 ± 0.030.02 ± 0.000.013 ± 0.000.005 ± 0.00
216.93 ± 0.218.49 ± 0.110.47 ± 0.07142.67 ± 2.527.67 ± 0.031.92 ± 0.071.27 ± 0.010.88 ± 0.010.09 ± 0.020.02 ± 0.000.014 ± 0.000.003 ± 0.00
516.78 ± 0.248.42 ± 0.180.70 ± 0.00142.33 ± 2.087.65 ± 0.042.20 ± 0.301.36 ± 0.030.88 ± 0.010.10 ± 0.020.03 ± 0.000.013 ± 0.000.004 ± 0.00
1016.78 ± 0.268.42 ± 0.211.03 ± 0.00142.33 ± 2.087.63 ± 0.041.98 ± 1.331.39 ± 0.090.89 ± 0.000.10 ± 0.020.02 ± 0.000.012 ± 0.000.004 ± 0.00
November014.44 ± 0.078.06 ± 0.180.62 ± 0.00146.33 ± 0.587.55 ± 0.052.99 ± 0.201.18 ± 0.090.88 ± 0.020.06 ± 0.010.03 ± 0.000.008 ± 0.000.008 ± 0.00
214.45 ± 0.068.05 ± 0.220.80 ± 0.06146.00 ± 1.007.61 ± 0.042.69 ± 0.501.13 ± 0.030.87 ± 0.010.06 ± 0.010.03 ± 0.000.008 ± 0.000.008 ± 0.00
514.45 ± 0.068.03 ± 0.251.13 ± 0.07146.33 ± 0.587.60 ± 0.033.30 ± 0.141.18 ± 0.010.87 ± 0.000.07 ± 0.010.03 ± 0.000.009 ± 0.000.010 ± 0.00
1014.46 ± 0.078.01 ± 0.221.57 ± 0.14146.00 ± 1.007.57 ± 0.023.28 ± 0.111.20 ± 0.060.90 ± 0.030.06 ± 0.010.03 ± 0.000.013 ± 0.000.010 ± 0.00
December011.07 ± 0.208.99 ± 0.610.17 ± 0.12157.33 ± 4.167.72 ± 0.053.06 ± 0.391.21 ± 0.100.94 ± 0.060.06 ± 0.010.01 ± 0.000.038 ± 0.000.016 ± 0.00
211.02 ± 0.169.09 ± 0.530.53 ± 0.07155.33 ± 1.157.73 ± 0.083.01 ± 0.301.14 ± 0.030.86 ± 0.010.06 ± 0.010.01 ± 0.000.036 ± 0.000.022 ± 0.00
511.01 ± 0.169.30 ± 0.731.00 ± 0.07155.33 ± 1.157.69 ± 0.072.96 ± 0.211.17 ± 0.030.86 ± 0.000.06 ± 0.010.02 ± 0.000.034 ± 0.000.023 ± 0.00
1011.00 ± 0.169.03 ± 0.941.57 ± 0.06156.33 ± 2.087.70 ± 0.052.71 ± 0.411.14 ± 0.000.86 ± 0.010.06 ± 0.010.02 ± 0.000.032 ± 0.000.024 ± 0.00
Two-way ANOVAMonth******NS************************
DepthNSNS*NSNSNSNSNSNSNSNSNS
Month × DepthNS***NS**NSNSNS*NSNSNSNS
T = water temperature; DO = dissolved oxygen; TN = total nitrogen; TP = total phosphorus. Data showed are means ± standard deviations (n = 3). * p < 0.05; ** p < 0.01; *** p < 0.001 represent statistical significance using two-way ANOVA. NS means no statistical significance.
Table 2. Three individual network properties of phytoplankton communities at the genus level.
Table 2. Three individual network properties of phytoplankton communities at the genus level.
ParametersOctoberNovemberDecember
Average degree0.8240.7780.875
Average weighted degree0.8050.7580.843
Network diameter323
Graph density0.0510.0460.058
Modularity0.730.7960.597
Connected components665
Average clustering coefficient0.1180.1670.062
Average path length1.4351.2221.5
Nodes171816
Edges151514
Table 3. Eukaryotic community diversity and richness estimators of Jinpen reservoir from October to December in 2018.
Table 3. Eukaryotic community diversity and richness estimators of Jinpen reservoir from October to December in 2018.
MonthWater Depth (m)Reads NumberOTUs0.97 Level
Chao1Shannon Diversity (H’)Simpson Diversity (D)Coverage (%)
October048,084365425 (398, 497)2.92 (2.9, 2.94)0.215 (0.211, 0.219)99
248,874218228 (221, 254)2.38 (2.36, 2.4)0.311 (0.306, 0.315)99
558,482256267 (259, 291)3.26 (3.24, 3.28)0.165 (0.163, 0.168)99
1048,125222236 (225, 278)3.15 (3.13, 3.17)0.181 (0.178, 0.185)99
November036,313132143 (135, 176)2.12 (2.1, 2.15)0.381 (0.375, 0.387)99
244,037304308 (305, 322)4.3 (4.29, 4.31)0.026 (0.025, 0.026)99
531,597415476 (450, 523)4.18 (4.16, 4.2)0.034 (0.033, 0.035)99
1047,074378406 (392, 435)3.37 (3.35, 3.39)0.151 (0.147, 0.154)99
December030,176169184 (174, 218)3.8 (3.79, 3.82)0.052 (0.050, 0.053)99
233,106393444 (421, 487)4.17 (4.16, 4.19)0.033 (0.032, 0.034)99
534,724365429 (401, 478)3.58 (3.56, 3.6)0.098 (0.096, 0.100)99
1043,256393443 (421, 484)3.98 (3.97, 4)0.047 (0.046, 0.048)99
Table 4. Network parameters of eukaryote community compositions based on the 18S DNA sequence at the genus level (R = 0.8).
Table 4. Network parameters of eukaryote community compositions based on the 18S DNA sequence at the genus level (R = 0.8).
ParametersOctoberNovemberDecember
Average degree2.19423.138
Average weighted degree2.1151.9282.99
Network diameter14612
Graph density0.0330.0280.049
Modularity0.7430.7890.584
Connected components762
Average clustering coefficient0.2240.2840.11
Average path length3.7772.2284.881
Nodes677265
Edges147144204

Share and Cite

MDPI and ACS Style

Yan, M.; Chen, S.; Huang, T.; Li, B.; Li, N.; Liu, K.; Zong, R.; Miao, Y.; Huang, X. Community Compositions of Phytoplankton and Eukaryotes during the Mixing Periods of a Drinking Water Reservoir: Dynamics and Interactions. Int. J. Environ. Res. Public Health 2020, 17, 1128. https://doi.org/10.3390/ijerph17041128

AMA Style

Yan M, Chen S, Huang T, Li B, Li N, Liu K, Zong R, Miao Y, Huang X. Community Compositions of Phytoplankton and Eukaryotes during the Mixing Periods of a Drinking Water Reservoir: Dynamics and Interactions. International Journal of Environmental Research and Public Health. 2020; 17(4):1128. https://doi.org/10.3390/ijerph17041128

Chicago/Turabian Style

Yan, Miaomiao, Shengnan Chen, Tinglin Huang, Baoqin Li, Nan Li, Kaiwen Liu, Rongrong Zong, Yutian Miao, and Xin Huang. 2020. "Community Compositions of Phytoplankton and Eukaryotes during the Mixing Periods of a Drinking Water Reservoir: Dynamics and Interactions" International Journal of Environmental Research and Public Health 17, no. 4: 1128. https://doi.org/10.3390/ijerph17041128

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop