Succession of Microbial Community in a Small Water Body within the Alluvial Aquifer of a Large River

Nitrogen is one of the essential elements limiting growth in aquatic environments. Being primarily of anthropogenic origin, it exerts negative impacts on freshwater ecosystems. The present study was carried out at the nitrate-vulnerable zone within the alluvial aquifer of the large lowland Drava River. The main aim was to investigate the ecosystem’s functionality by characterizing the bacterial and phytoplankton diversity of a small inactive gravel pit by using interdisciplinary approaches. The phytoplankton community was investigated via traditional microscopy analyses and environmental DNA (eDNA) metabarcoding, while the bacterial community was investigated by a molecular approach (eDNA). Variations in the algal and bacterial community structure indicated a strong correlation with nitrogen compounds. Summer samples were characterized by a high abundance of bloom-forming Cyanobacteria. Following the cyanobacterial breakdown in the colder winter period, Bacillariophyceae and Actinobacteriota became dominant groups. Changes in microbial composition indicated a strong correlation between N forms and algal and bacterial communities. According to the nitrogen dynamics in the alluvial aquifer, we emphasize the importance of small water bodies as potential buffer zones to anthropogenic nitrogen pressures and sentinels of the disturbances displayed as algal blooms within larger freshwater systems.


Introduction
Various aspects of nutrient dynamics in freshwater ecosystems are of paramount importance for understanding how the productivity of surface waters is controlled and provide the opportunity to analyse the current and future impacts of anthropogenic activities on freshwater ecosystems. In such environments, a large part of the primary production may depend on the recycling of nutrients such as nitrogen compounds [1]. Nitrogen is an essential element that often limits growth in aquatic ecosystems, and a key compound in many biochemical processes that are important for life, but can be harmful in high concentrations [2][3][4]. Nowadays, anthropogenic activities such as fertilizer synthesis and its widespread application on arable areas, as well as the burning of fossil fuels, significantly increase the N fluxes across different environmental compartments [3][4][5]. Its environmental 200 km 2 ( Figure 1). On the NW side, the alluvial aquifer is adjoined by the Varaždin Lake, an artificial reservoir of the hydroelectric power plant Varaždin, of which the Drava River watercourse constitutes a natural border in the NE direction. The Plitvica Stream flows at the S-SE edge of the study area, while in the middle there are several active and inactive gravel pits. All gravel pits represent exposed groundwater and vulnerable areas where the contamination of groundwater can occur faster from surface contaminants. In some cases, inactive (abandoned) gravel pits are used as waste (industrial or urban) disposal sites and are becoming a threat to groundwater quality [37][38][39].
(NW Croatia), belongs to the Black sea catchment and covers an area of approximately 200 km 2 ( Figure 1). On the NW side, the alluvial aquifer is adjoined by the Varaždin Lake, an artificial reservoir of the hydroelectric power plant Varaždin, of which the Drava River watercourse constitutes a natural border in the NE direction. The Plitvica Stream flows at the S-SE edge of the study area, while in the middle there are several active and inactive gravel pits. All gravel pits represent exposed groundwater and vulnerable areas where the contamination of groundwater can occur faster from surface contaminants. In some cases, inactive (abandoned) gravel pits are used as waste (industrial or urban) disposal sites and are becoming a threat to groundwater quality [37][38][39].
The study area is densely populated, with industrial and intensive agricultural production. The most common type of crops grown are corn, cabbage, potatoes, and vegetables. Extensive poultry farming is present, especially the fattening of chickens, quails, and pheasants, and the breeding of hens [40]. The Varaždin pumping site, one of two in the area, was shut down due to high nitrogen concentrations in groundwater caused by significant anthropogenic activities [41]. Nevertheless, demands for drinking and industrial water rise because of the growing production in the area. The gravel pit in the village of Šijanec was chosen because of its inactivity and accessibility. It is a small pit covering an area of approximately 12,000 m 2 ( Figure 1).  The study area is densely populated, with industrial and intensive agricultural production. The most common type of crops grown are corn, cabbage, potatoes, and vegetables. Extensive poultry farming is present, especially the fattening of chickens, quails, and pheasants, and the breeding of hens [40]. The Varaždin pumping site, one of two in the area, was shut down due to high nitrogen concentrations in groundwater caused by significant anthropogenic activities [41]. Nevertheless, demands for drinking and industrial water rise because of the growing production in the area. The gravel pit in the village of Šijanec was chosen because of its inactivity and accessibility. It is a small pit covering an area of approximately 12,000 m 2 ( Figure 1).

Sampling and Laboratory Analysis
Phytoplankton and bacterial samples were taken monthly from June 2017 until March 2018 on the deepest point of the gravel pit using the vertical Hydro-Bios water sampler (Hydro-Bios Apparatebau GmbH, Altenholz, Germany). Samples for chemical analysis were taken simultaneously with biological samples and transported in a portable freezer at 4 • C to the laboratory for further analysis. Before sampling, physico-chemical parameters of water, including electrical conductivity (EC), pH, temperature (T), and dissolved oxygen concentration (DO), were measured with a portable WTW Multi 3630 multimeter (Xylem Analytics, Weilheim, Germany). Water transparency (Z SD ) was estimated using a Secchi Disc. Total alkalinity was measured by titration with 1.6 N H 2 SO 4 using phenolphthalein and bromocresol green-methyl red as indicators. Dissolved cations (Na + , K + , Mg 2+ and Ca 2+ ) and anions (SO 4 2− , NO 3 − , Cl − ) were analyzed by ion chromatography using a Dionex ICS-6000 (Thermo Fisher Scientific, Waltham, MA, USA), while NH 4 + , NO 2 − , and PO 4 3− -P were measured spectrophotometrically. Dissolved inorganic and organic carbon (DIC, DOC) and total inorganic and organic carbon (TIC, TOC) were analysed using a HACH QbD1200 TOC analyser (Hach Company, Frederick, MD, USA). The analytical precision of the measurements of cations and anions, indicated by the ionic balance error (IBE), was computed on the basis of ions expressed in meqL −1 . The IBE value was observed to be within a limit of <±5% [42,43]. The PHREEQC geochemical code [44] was used to calculate the charge balance and pCO 2 pressure and to study the saturation state of the mineral phases. The samples for Chlorophyll a were filtered on Whatman GF/F glass filters (Whatman International Ltd, Kent, UK), extracted in 96% ethanol, and measured for Chl a using a UV-VIS spectrophotometer according to compliance monitoring standards [45].

Microbial Community Analysis
Phytoplankton samples were collected for both morphological and molecular analyses. Samples for morphological analysis were fixed with acid Lugol solution and stored in 250 mL volume glass bottles in the dark at 4 • C. The morphological analysis included a qualitative and quantitative community characterization according to the Utermöhl method [46], using a Zeiss Axiovert 200 inverted microscope (Carl Zeiss, Oberkochen, Germany).
Samples for molecular analysis were collected in sterile plastic bottles and preserved on ice during the transport to the laboratory. In the laboratory, they were immediately filtered on Nucleopore Track-Etch membrane filters (25 mm diameter, 0.4 µm pore size; Whatman International Ltd, Kent, UK) and stored at −20 • C until further processing. DNA was extracted from the frozen filters using DNeasy PowerWaterKit (Qiagen, Hilden, Germany). The manufacturer's instructions were followed with a slight modification in the final step, where 60 µL of sterile DNA-free PCR Grade water was added instead of Qiagen's C6 Solution. The quality of the extracted DNA was assessed with the Shimadzu BioSpec-nano spectrophotometer (Shimadzu Corporation, Kyoto, Japan).

PCR of the Phytoplankton (Eukaryotic) Community
The hypervariable V9-region of the SSU rRNA gene (ca. 150 base pairs) was amplified from environmental DNA using the universal eukaryotic primer pair [47,48]. The forward and reverse primers were 1391F (5 -GTAC ACACCGCCCGTC-3 ) and EukB (5 -TGATCCTTCTGCAGGTTCACCTAC-3 ), designed by Amaral-Zettler and colleagues [49]. Polymerase chain reactions (PCR) contained 1 U of Hot Start Taq DNA Polymerase (New England Biolabs, Ipswich, MA, USA), and for V9 amplification an initial activation step at 95 • C was employed for 5 min, followed by 30 three-step cycles consisting of 94 • C for 30 s, 57 • C for 45 s, and 72 • C for 1 min, followed by a final 2 min extension at 72 • C. PCR products were assessed by visualizing on a 1% agarose gel. Sequencing libraries were constructed using the NEB Next ® Ultra™ DNA Library Prep Kit for Illumina (New England Biolabs, Ipswich, MA, USA). Libraries were sequenced on an Illumina MiSeq platform, generating 250-bp paired-end reads (SeqIT GmbH & Co. KG, Kaiserslautern, Germany).

PCR for the Bacterial Community
The V3-4 region of bacterial rRNA genes was amplified using the forward primer 341 F 5 -CCTACGGGNGGCWGCAG-3 and reverse primer 805 R 5 -GACTACHVGGGTATCTAA TCC-3 [50]. PCR and sequencing were performed in the LGC Genomics GmbH laboratory (Berlin, Germany). The libraries were sequenced on an Illumina MiSeq platform, generating 300-bp paired-end reads.

Bioinformatic Analysis of the Phytoplankton (Eukaryotic) Community
Sequence reads were analysed using the programs implemented in the OBITools package, as described in De Barba et al. [51]. The quality of the reads was assessed using FastQC. Paired-end reads were aligned using Illumina paired-end, and alignments with quality scores >32 were kept. The aligned data set was demultiplexed using the ngsfilter command, which identified primers and tags and assigned the sequences to each sample. For dereplication, we used the obiuniq command for clustering together strictly identical sequences and keeping the information about their distribution among samples. Sequences shorter than 10 bp, or containing ambiguous nucleotides, or with occurrence lower or equal to 10 were excluded using the obigrep command. The obiclean command was then run to assign the status of "head", "internal", or "singleton" to each sequence within a PCR product. All sequences labeled "internal", corresponding to PCR errors, were discarded. Finally, the taxonomic assignment of OTUs was performed with the ecotag command, combining two reference databases, filtered according to target taxa from NCBI taxonomy and the EMBL database, after running the ecoPCR program [52,53]. Only sequences with a 98% match to the reference sequence were kept. Single-read OTUs were removed from the samples to avoid potential false positives. The final filtered file with taxonomically assigned OTUs of eukaryotic algae groups was used as a basis for all downstream analyses. The DNA sequencing reaction on two samples (September and October 2017) did not yield valid results. A list of commands with related parameters are presented in the Supplementary Materials (S1). Raw demultiplexed reads were deposited at the ENA's Sequence Read Archive and are publicly available under project number PRJEB40961.

Bioinformatic Analysis of the Bacterial Community
The quality of the reads was assessed using FastQC. Paired-end reads were qualitytrimmed using the bbduk function and merged using the bbmerge function of the BBMap package 34.48 (Lawrence Berkeley National Lab., Berkeley, CA, USA) [54]. Merged reads were quality-filtered using QIIME v1.8.0 [55]. Reads with exact barcodes and primers, unambiguous nucleotides, and a minimum length of 250 base pairs were retained. A Chimera check was done using UCHIME [56]. Non-chimeric reads were clustered with SWARM v3.0.0 [57] with default settings into Operational Taxonomic Units (OTUs). The bacterial reads were blasted against the SILVA database release 138 (Max Planck Institute for Marine Microbiology and Jacobs University, Bremen, Germany) using blastn (BLAST v2.9.0) [58]. Nontarget OTUs (chloroplasts, mitochondria), as well as singletons and doubletons, were excluded. The resulting OTUs were filtered by the quality of the blast results (≥98% identity). The DNA sequencing reaction on two samples (February and March 2018) did not yield valid results. Raw demultiplexed reads were deposited at the ENA's Sequence Read Archive and are publicly available under project number PRJEB40962.

Statistical Analysis
Statistical analyses were conducted in R v. 4.0.2 [59] using the program packages "vegan", "fossil", "factoextra", "devtools", and "ggbiplot", as well as "ggplot2" for all graphical representations. To access the comparability of morphological and molecular methodologies in the phytoplankton community, the taxa lists derived from both approaches were compared with regard to the presence or absence of taxa and community composition. The bacterial community composition was analysed by using the molecular approach. The results for downstream analysis were combined into a single data set for each approach and for each community. The molecular results were transformed into relative abundances to normalize the OTU database [60]. Biomass data obtained by microscopy for the phytoplankton community were transformed following the logarithmic scale [61].
The Shannon, Simpson, and richness indices were calculated for both approaches and both communities as measures of alpha diversity using program packages "Vegan v. 2.5.6" [62].
To test the statistical significance of the environmental parameters and which parameter was singled out depending on the month studied, a principal component analysis (PCA) was performed using the R package "vegan" [62].
Canonical correspondence analyses (CCA) were performed on phytoplankton morphological data and bacterial molecular data to estimate variance in environmental variables for both communities. ANOVA test was applied to test the statistical significance of all axes, and forward selections were used to evaluate the importance of each variable. The logarithm function was used to transform environmental parameters and both community datasets for statistical analysis.

Analysis of Environmental Parameters
The environmental variables of the investigated gravel pit are indicated in Table 1. The highest value of nitrates (NO 3 − ) concentration was measured in March (38.4 mg L −1 ) and the lowest in June (0.62 mg L −1 ), whereas the maximum concentrations of ammonium (NH 4 + ) and nitrites (NO 2 − ) were recorded in July (2.75 mg L −1 and 0.17 mg L −1 , respectively). The highest value of pH was detected in August (9.48), indicating that the water was alkaline. The highest temperature (T) value was recorded in June (24.2 • C) and the lowest in February (1.2 • C), respectively. The maximum value of dissolved oxygen was measured in August (17.1 mg L −1 ), and the minimum value was in July (7.1 mg L −1 ). For electrical conductivity (EC) and bicarbonates (HCO 3 − ), the maximum values were recorded in March (497 µS cm −1 and 249 mg L −1 , respectively) and the minimum ones in October (252 µS cm −1 and 107 mg L −1 , respectively). Silicon dioxide (SiO 2 ) concentrations were high during the warmer period, and the maximum was in September (26.8 mg L −1 ). During the colder period the concentrations were much lower, and the minimum was in December (7.2 mg L −1 ). In addition, a change was also observed in the concentrations of calcium, bicarbonates, and logpCO 2 pressure (Table 1), with the lowest concentrations detected in the summer period and the highest concentrations during winter.

PCA Analysis
The principal component analysis (PCA) performed for the 24 environmental variables explained 71.1% of the total variance on the first two PC axes. The overall strength of correlations between the samples and environmental parameters are summarized in Table 2. The most important variables for the PC1 axis were TIC, Ca 2+ , and EC (intra-set correlations: 0.271, 0.269, and 0.268, respectively). Regarding the PC2 axis, NH 4 + and DO (intra-set correlations: −0.355 and 0.407, respectively) were the variables that weighted most for the ordination. PCA arranged the samples into three groups (Figure 2): the first group consisted of samples from a warmer period of investigation (June, August, September, and October), the second group included a sample from July, while the third one comprised all samples from the colder period of investigation (November, December, January, February, and March).

Phytoplankton Succession
A total of 38 phytoplankton species were recorded by a morphological approach within the 10 samples collected during the investigated period. A total of 47,130 reads clustered into 88 OTUs were detected in the remaining eight samples. OTUs were taxonomically assigned into 30 eukaryotic algal taxa.
Phytoplankton abundance ranged between 1.9 × 10 6 cells L −1 in January 2018 to 2.8 × 10 8     A CCA was performed for the phytoplankton samples, and nine constrained environmental variables ( Figure 5) indicated eigenvalues for the first two axes of 0.4328 and 0.3210, respectively, explaining 38.5% of the total variance on the first two axes. A Pearson environment-species permutation for the two significant axes indicated a significant correlation between abiotic constrained values and phytoplankton functional variables. According to the ANOVA permutation test, the ordination of both axes for environmental variables was statistically significant (p < 0.05). Canonical coefficients and intra-set correlations on the phytoplankton samples showed that NO 3 − and EC were the most important variables for the ordination axis 1 (intra-set correlation coefficients: 0.6332 and 0.5467, respectively). Regarding axis 2, NO 2 − and NH 4 + (intra-set correlations −0.8155 and −0.6837, respectively) were the variables that weighted most for the ordination. At the positive end of both axes, phytoplankton samples were associated with EC, DO, HCO 3 − , and NO 3 − . At the negative end of both axes, phytoplankton samples were associated with pH, T, SiO 2 , NH 4 + , and NO 2 − . Considering the environmental pressure to phytoplankton, the CCA analysis showed the separation of samples into three groups ( Figure 5). The first group, comprised of summer and autumn samples (July to October 2017), correlated with high concentrations of NH 4 + (2.75 mg L −1 ), NO 2 − (0.173 mg L −1 ), and SiO 2 (26.8 mg L −1 ) and high values of pH (9.48) and T (24 • C). The most common species of the group were cyanobacteria M. aeruginosa and M. wesenbergii, green alga Scenedesmus quadricauda and cryptophyte Cryptomonas sp. According to the morphological approach, the samples collected in July were characterized by the highest Chl a concentration, the dominance of S. quadricauda, and the highest concentration of NH 4 + (2.75 mg L −1 ). A pronounced increase in the total phytoplankton biomass was recorded in August as a result of Microcystis spp. bloom. M. aeruginosa and Cryptomonas sp. were the descriptive species of the phytoplankton community in September and October, with the continuing decrease of the phytoplankton biomass. The molecular approach did not confirm the same composition, but rather detected the cyanobacterial predominance and a higher contribution of OTUs taxonomically assigned to the Cryptomonas genera. According to the morphological approach, the second group, consisting of an outlying sample from June 2017, was characterized by the predominance of filamentous cyanobacterium Limnotrix redekei (Goor) Meffert and centric diatom Aulacoseira muzzanensis (Meister) Krammer. The molecular approach confirmed the dominance of OTUs taxonomically assigned to the Aulacoseira genera as well, but also detected a high number of reads of OTUs taxonomically assigned to genus Parvodinium. The third group, composed of winter samples (November 2017 to March 2018), correlated with low T and higher concentrations of NO 3 − (38.4 mg L −1 ) and HCO 3 − (249 mg L −1 ) and with mostly constant values of DO (11.9 to 14.7 mg L −1 ) and EC (497 µS cm −1 ). As confirmed by both morphological and molecular approaches, the most dominant species/taxonomically assigned OTUs during the colder period were pennate diatoms (Fragilariaceae) of the genus Ulnaria. Dinoflagellate Peridinium sp. was recorded only by the morphological approach as the subdominant species in November and January. The higher contribution of the genus Dinobryon was confirmed by both analyses during February and March.
As confirmed by both the morphological and molecular approaches, the most dominant species during the colder period were pennate diatoms (Fragilariaceae) of the genus Ulnaria with dinoflagellate Peridinium sp. as the subdominant species.

Bacterial Community Composition
For the eight samples, a total of 582,900 reads clustered into 1743 OTUs were recorded. Sequence reads were taxonomically assigned into 52 bacterial phyla. The bacterial composition showed a succession of differences between months. The variations of alpha diversity in the rarefied richness, Shannon, and Simpson indices are shown in Figure 6. The highest bacterial richness was recorded in July 2017, and the lowest in October 2017. The Shannon and Simpson indices showed similar results, with the highest values in January 2018 and the lowest in August 2017. Comparing the composition of bacteria and eukaryotic algae inferred from the molecular approach, low values in richness and both indices during the summer period were noted.
The CCA analysis performed on the bacterial community and seven constrained environmental variables (Figure 7) indicated eigenvalues for the first two axes of 0.7131 and 0.6563, respectively, explaining 47.2% of the total variance on the first two axes. A Pearson environment-bacterial community permutation for the two significant axes indicated a significant correlation between abiotic constrained values and bacterial functional variables. According to the ANOVA permutation test, the ordination of both axes for environmental variables was statistically significant (p < 0.05). Canonical coefficients and intra-set correlations carried out on the bacterial community samples showed that pH and DO were the most important variables for the ordination axis 1 (intra-set correlation coefficients: 0.6140 and 0.3770, respectively). Regarding axis 2, NO 2 − and NH 4 + (intra-set correlations: 0.9229 and 0.8706, respectively) were the variables that weighted most for ordination. At the positive end of both axes, the bacterial community's samples were associated with T, pH, and NO 2 − . At the negative end of both axes, the bacterial community's samples were associated with EC and HCO 3 − . Based on the position of the samples related to the CCA1 axis, the bacterial community separated into two groups (Figure 7). The first group, comprised of summer samples (June to September 2017), correlated with a high concentration of NH 4 + and NO 2 − , and high values of pH and T. According to the percentage of classified OTUs, the most common phyla in the first group were Cyanobacteria, with the most dominant families being Pseudanabaenaceae (June 2017), Microcystaceae (July 2017), and Synechococcaceae and Microcystaceae (August 2017). This was also confirmed by a morphological approach, except for the picocyanobacterium of the genus Cyanobium, which is hard to detect under light microscopy. In June, Planctomycetota were codominant with Cyanobacteria. The sample from July was singled out as a result of the codominance of Cyanobacteria and Bacteroidota, which correlated with the concentration of NH 4 + . The second group, comprising samples from the colder months (October 2017 to January 2018), correlated with low T and higher concentrations of HCO 3 − , DO, and EC. Actinobacteriota was the dominant phylum in October 2017 and January 2018, whilst samples from November and December were dominated by Proteobacteria.

Discussion
Data on ecology and the importance of small water bodies in alluvial lowlands are still quite scarce, as those systems are not included in the national strategies for the protection of water resources. Within larger freshwater systems, small water bodies act as biochemical reactors because of their potential for supporting high metabolic rates that are often paired with naturally high concentrations of nutrients and trophic conditions [33,34,63]. Their role is modulating nutrient retention and recycling along the hydrological pathways [63]. Even though usually related to eutrophic or hypertrophic conditions, small standing water bodies collectively support a very diverse biodiversity, often with species-richer communities more adapted to eutrophication conditions and to a broad range of physical and chemical conditions than the communities in larger water bodies [38,64]. Eutrophication has been described as a major stressor for the freshwater biodiversity of both large and small water bodies [65,66]. Rosset et al. [32] suggested that the eutrophication management of lowland small water bodies should be regulated differently than for larger freshwater systems, with the conservation efforts focused on the protection of small water bodies representing a mosaic of trophic conditions (and acting as hosts of regional biodiversity).
Seasonal changes with complex dynamic phases govern the high rates of biodiversity in small lowland water bodies altered by high anthropogenic pressure and climate-related impacts, such as the Šijanec gravel pit [67,68]. Previous studies on the Drava River lowland did not consider the importance of small water bodies within the whole alluvial system [69,70]. Due to hydromorphological characteristics of the catchment, Šijanec receives a high nitrogen input from the Drava River aquifer via groundwater recharging [41]. The high concentration of nitrates in the Drava River groundwater system, with an average of 60.9 mg L −1 , was observed by Marković et al. [41]. Nitrogen compounds can easily percolate through the soil into the groundwater either from direct terrestrial runoff or with rainfall or irrigation water [71]. A higher nitrogen concentration was noted during the colder seasons as a consequence of the recorded decrease in the phytoplankton biomass and abundance, as well as the rise of the groundwater level due to an increase in precipitation [72,73], whereas nitrogen concentrations in the groundwater dropped during the warmer periods following the decrease of precipitation, as confirmed by the PCA analysis. Since nitrogen fertilizers are widely used in agriculture to increase crop production, the cropping practices and soil texture have been found to influence the extent of nitrate leaching [74].
The nutrient-based indication of eutrophic conditions in the gravel pit was further supported by high values of phytoplankton biomass and bacterial density throughout the investigated period, especially during the summer months, when the lowest NO 3 − concentration (June) and the highest concentrations of NH 4 + and NO 2 − (July) were detected. This occurrence of elevated nitrogen is presumably occurring as an effect of the cropping season, e.g., fertilizer application through intense irrigation [75]. This is particularly true for nitrates, which are normally assumed not to be absorbed by soil particles and are therefore easily leached, in which case the nitrate distribution should follow the wetting front [76,77]. As found by Paredes et al. [78], high intra-and interannual hydrological fluctuations influenced nitrate occurrence in freshwater streams and ponds, with the main source of nitrate linked to agricultural practices and the use of both organic and synthetic fertilizers. Since it can be rapidly oxygenized, the concentration of NO 2 − is usually deficient [79], as was exhibited in this study. Most of the nitrogen uptake in shallow eutrophic systems is the form of nitrate, which has a positive effect on the growth of phytoplankton biomass [80,81]. A sample from June had the lowest concentrations of NO 3 − , TN, TIC, and logpCO 2 , but the highest temperature and Si values. These conditions were characterized by the dominance of cyanobacteria Limnotrix redekei and centric diatom Aulacoseira muzzanensis. Limnotrix redekei, a species characteristic of eutrophic shallow water bodies used for recreation and fishing, often shows domination in spring and summer periods with co-occurrence of centric diatoms [82]. This species is known for its ability to adapt to low-light, cold conditions and capability to overwinter in considerable densities under the long-term ice-and snow-cover [83][84][85][86]. Aulacoseira muzzanensis, a species adapted to live in turbid and nutrient-rich waters [87], was described in a hyper-eutrophic lake (Lago di Muzzano) located in the Tessin region of Switzerland [88], which suffers from periodic Microcystis algal blooms [89][90][91][92]. During lower light conditions, both taxa can occur concomitantly with the Microcystis species in quantities capable of eliminating other phytoplankton taxa [93]. A sample from July in the grouped composed of summer and autumn samples was characterized by maximum concentrations of nitrites and ammonium. This suggested enhanced growth conditions for specific algal groups under a higher NH 4 + supply, which was also evident by a recorded high concentration of the chlorophyll a. Certain phytoplankton taxa, such as cyanobacteria and chlorophytes, prefer a high supply of energetically favorable NH 4 + [94], as was also noted in our investigation. These conditions can also inhibit NO 3 − uptake for other taxa, such as large diatoms [94][95][96]. The most abundant chlorophyte was Scenedesmus quadricauda, a small coenobium-forming and ammonium-tolerant species [97][98][99] with higher uptake abilities for ammonium under nitrogen limitation than species of the genus Microcystis and with competitive superiority in the large-pulse, low-frequency nutrient recharging [100], as was present in our study area. With regard to size, small algae can be more competitive than larger species, because they have high surface area to volume ratios, resulting in greater specific uptake and growth under low nutrient concentrations, when even slight NH 4 + additions can be enough to promote the growth of small algae [101]. Secondly, the most dominant species during the summer period was Microcystis spp., with the culmination of dominance (surface bloom) in August, which was connected with the maximum of phytoplankton biomass and number of cells per liter. Microcystis is a non-N 2 -fixing cyanobacteria, which dominates in highly eutrophicated, stratified ponds, rivers, or lakes receiving elevated N loadings, since its growth is dependent on nitrogen supplies [6]. A strong bloom of Microcystis spp. was associated with the warmer period, low nitrogen concentrations, and the subsequent water level drawdown [102]. Cyanobacteria have the ability to adapt to different environments by adjusting their light harvesting and carbon fixation mechanisms. Adversely, a high rate of photosynthesis induced by the Microcystis bloom can considerably reduce dissolved CO 2 concentrations and drive up the pH value [103]. Furthermore, Microcystis favors more alkaline conditions as a competitive advantage over other phytoplankton groups [104,105]. Along with the Microcystis bloom, the molecular approach detected the cyanobacterium of the genus Cyanobium, which was dominant in August. Cyanobium is a picocyanobacterial genus with a presumably significant role in the functioning of aquatic ecosystems, but rather hard to detect microscopically due to its size and taxonomical obscurity. Both Microcystis and Cyanobium genera are less demanding on nutrients and generally demonstrate summer peaks when the concentrations of nitrogen compounds are usually lower [106]. As detected by the molecular approach, the phylum Planctomycetota was subdominant in the bacterial community in June. Members of the phylum Planctomycetota have been found in a variety of environments, including freshwater [107], and are known for their role in the anaerobic oxidation of ammonium (anammox), as part of the biogeochemical nitrogen cycle [108][109][110][111][112]. In oxygen-limited systems, such as biofilm aggregates, the planctomycete anammox bacteria live closely associated with aerobic ammonium oxidizers. The aerobic ammonium-oxidizing bacteria consume oxygen at the outside of the biofilm, thus keeping the inside anoxic for the anammox bacteria. Together, they create conditions in which they can convert ammonium directly into dinitrogen gas. Anammox bacteria can contribute significantly to the loss of fixed nitrogen in both natural and anthropogenicinfluenced ecosystems [110,112,113]. The members of the phylum Bacteroidota recorded in July are typical for freshwater environments [114]. Their dominance correlated with OTUs taxonomically assigned in eukaryotic biflagellate cryptophytes from Cryptomonas genera, and with chlorophyte Scenedesmus quadricauda, as recorded by the morphological approach. The increase in their abundance correlated with higher algal concentrations, presumably due to their ability to establish a mutualistic relationship on the algal cell surface [115]. High densities of the Cryptomonas genus can occur following the period of nutrient enrichment [6]. During the colder period, the composition of the phytoplankton switched from cyanobacteria to a diatom-dominated community characterized by the genus Ulnaria, whose winter blooms require both biogenic silica for the formation of their outer cell wall structures (frustules) and lower basic pH conditions, as not to corrode them [116] and decrease the growth rate [117]. This resulted in a threefold drop in the SiO 2 concentration which was presumably consumed in the building of their frustules [118,119]. Diatoms are characterized as effective nitrate utilizers with high preferences to NO 3 − uptake [120]. Also subdominant in the colder period, together with diatoms, were dinoflagellate species, which were reported to have significant NO 3 − uptake rates [121]. A dense bloom of colonial chrysophyte Dinobryon spp. characterized the phytoplankton community in February and March. As adaptations to the lower temperature, the winter blooms of Dinobryon could indicate enhanced nutrient loading [6], but also the ability to obtain nutrients from bacteria by mixotrophy [122,123]. Members of Proteobacteria, as the most dominant bacterial phylum in the samples from November and December, are ubiquitous in freshwater environments, specifically in eutrophic conditions with high phytoplankton numbers [114,124]. Actinobacteriota were present throughout the investigated period, with increased abundance in October and January during the low abundance of cyanobacteria. This is plausibly correlated with the sensitivity of Actinobacteriota to conditions prevailing during the cyanobacterial blooms, such as high organic matter, inorganic nutrient availability, and high temperatures, under which the highly streamlined actinobacterial cells cannot compete [125]. Actinobacteriota can thrive under oligotrophic conditions due to their small size, high surface area-to-volume ratio, and enhanced capacity for efficient nutrient acquisition through high-affinity, broad-specific uptake systems [126]. Another important role of Actinobacteriota in freshwater habitats is connected to the proton-pumping rhodopsins (actinorhodopsins), thus revealing a photoheterotrophic lifestyle [127]. All these traits suggest that Actinobacteriota might serve as sentinels of impending ecological damage and have the potential to become standards of ecological freshwater quality [125]. They were all detected in the winter samples, except for Asulcocephalium miricentonis and Wislouchiella planctonica, which were detected in June. Besides type locality, no further ecological data were available on Asulcocephalium miricentonis, a species described in a temperate freshwater artificial pond in northeastern Japan [128]. Therefore, this record presents a contribution to the ecological conditions in which the species likely occurs. Wislouchiella planctonica is associated with man-made reservoirs and lentic freshwater habitats with eutrophic conditions [129,130], which is in line with our findings. Uroglenopsis americana was detected in February together with Dinobryon spp., due to its ability to compete for nutrients during the colder period, unlike algal species in eutrophic conditions [131]. The picoplanktonic species Meyerella planktonica is a major component of aquatic systems and a significant primary producer regularly occurring during winter [132,133]. Some of the species from the Chloromonas genera were found in snow samples [134], as was the case with the species Chloromonas subdivisa detected during snowy winter conditions in Šijanec. Actinochloris sphaerica is a cosmopolitan species mostly recorded in soil cultures and puddles [135]. Since very scarce information is available on all these species, the presented results also contribute to elucidating their ecological preferences.
Both approaches showed variations in diversity richness. Based on the morphological approach, the maximum value of alpha diversity was recorded in June, presumably due to the favorable conditions for algal growth. The minimum value of alpha diversity was detected in February, due to the lower number of algal species, whereas by using the molecular approach, higher values of alpha diversity were detected in the colder period, which can be associated with two possible causes: (1) when cell abundances of specific taxa in the water sample drop below a specific threshold, they can still be detectable with the molecular approach, but may not be found by microcopy; and (2) the resting stages of some algal species cannot be identified and assigned correctly by microscopy, but might be more easily recorded by the molecular approach [136]. Moreover, the lowest value of alpha diversity detected in August with the molecular approach could be correlated with the cyanobacterial bloom. Some cyanobacteria are known to produce toxins and cyclic peptides which can inhibit regulatory enzymes in eukaryotic cells, thus causing PCR inhibition [137,138]. Similar variations were detected in bacterial alpha diversity, with a maximum in July and low values throughout August, September, and October. This finding also corresponds to cyanobacterial bloom, which can inhibit the stabilization of microbial diversity [139]. Surprisingly, the morphological approach, as a basic descriptive method, succeeded in recognizing a higher microalgal diversity in Šijanec than the molecular approach, which is commonly considered a more powerful detector tool [140]. Events such as cyanobacterial blooms or discrepancies in the various DNA extraction methods can also have a discernible impact on the certainty of the community analysis via the molecular approach. Eukaryotic algae have a large range of cell wall structures, thus imposing challenges to the unbiased, uniform, and universal extraction of nucleic acids from such communities [138]. In spite of this, the molecular approach proved far more effective in discerning small-sized eukaryotic algae and cyanobacterial taxa, which are generally hard to detect with microscopy due to the scarcity of taxonomic knowledge and limitations of resolving power. Molecular methods showed that they can be successfully used to complement the morphological approach for assessing the microbial community's structure [138,141], especially in these kinds of eutrophic environments.

Conclusions
Despite the size, small water bodies like the Šijanec gravel pit play a key role as buffer zones within alluvial areas of large rivers. Due to the hydromorphological characteristics of the catchment, Šijanec receives a high nitrate input directly from the groundwater recharging. Nitrogen compounds likely control the phytoplankton biomass, thereby influencing the complete microbial community's structure. The integration of morphological and molecular approaches facilitates the comprehensive assessment of the microbial community's structure. Interdisciplinary approaches can be successfully used to elucidate the ecological preferences of microbial species and the prediction and prevention of algal blooms. The study emphasizes the importance of small water bodies in maintaining the state of water ecosystems and stresses the need for their enlistment in biomonitoring actions.