Bacterioplankton Community as a Biological Element for Reservoirs Water Quality Assessment

: The Water Framework Directive (WFD) is used to evaluate the water quality of aquatic ecosystems. Phytoplankton is the only biological element considered in the reservoirs water quality assessment. In this study, we aimed to assess the use of the bacterioplankton community as an indicator of water quality, using a culture-independent assay (denaturing gradient gel electrophoresis, DGGE). Four Portuguese reservoirs (Miranda, Pocinho, Aguieira and Alqueva) were analysed in four periods (autumn 2018, spring and autumn 2019, and spring 2020). Bacterial total abundance had similar values for Miranda, Pocinho and Aguieira, and generally lower values for Alqueva. Diversity and richness values did not show a clear trend. Negative correlations were observed between some nutrients and the bacterial community. Overall, members of Actinobacteria , Bacteroidetes , Cyanobacteria , Nitrospirae and Proteobacteria were identiﬁed in all sampling sites. In Alqueva, no spatial, temporal or water body quality relationships with bacterial community were observed, which may be due to its higher size, low water velocity rate and higher residence times. However, in Miranda, Pocinho and Aguieira, a strong spatial and temporal bacterial community dynamic was observed. Furthermore, the presence of some species (e.g., Acinetobacter sp.) may reﬂect the poor water quality that was not detected by the WFD approach.


Introduction
Superficial freshwaters, especially reservoirs, are subjected to eutrophication, mainly due to the excessive charge of nutrients, in particular, in the area upstream of the reservoir [1,2]. The occurrence of this phenomenon leads to excessive growth of primary producers, which can disrupt the dynamics of the biological equilibrium in the aquatic ecosystem [3].
The Water Framework Directive 2000/60/EC (WFD), created in 2000 and applied by the European Union, aims to assess the ecological quality of water bodies [4] which, in the case of heavily modified water bodies (reservoirs), is denominated as "ecological potential". This classification allows determining which parameters are in good conditions (similar to reference values) and which parameters are changed and therefore responsible for water degradation. As far as reservoirs are concerned, this directive is based on general physical and chemical parameters, on a specific list of specific pollutants and priority substances and on hydromorphological parameters. However, as biological elements only, the phytoplankton community is considered for this water typology.
Lentic freshwater environments have been heavily studied to understand their dynamics, focusing namely on physical and chemical characteristics. Relative to the community composition, phytoplankton and ichthyofauna are the most studied [5]. Of the freshwater ecosystems, reservoirs are the less studied and the most similar to natural lakes, with small differences, such as shorter retention times and water level fluctuations, which makes the Four Portuguese reservoirs, Miranda (M), Pocinho (P), Aguieira (Ag) and Alqueva (Al), were selected to conduct this study (Figure 1).
Miranda and Pocinho reservoirs are located in the north of Portugal, in the international and national sections of Douro River, respectively. Miranda reservoir has been in operation since 1960 to produce energy, and it is the smallest reservoir of this study, with an area of 1 km 2 [17]. Pocinho reservoir started operating in 1982, is used for navigation as well as hydroelectric exploration and occupies an area of 8 km 2 [18]. The selection of these reservoirs is related to the fact that Miranda and Pocinho are two ancient reservoirs, with small dimensions and located in a main course river in the northern region. Aguieira reservoir is located in the centre of Portugal, in the Mondego River, and started operating in 1981 for irrigation and hydroelectric exploration [19]. This reservoir occupies an area of 20 km 2 and was included in the WFD inter-calibration exercise [20]. Aguieira was selected because it is one of the reference reservoirs for the WFD calibration and is located in the centre of Portugal on a main course river. Alqueva reservoir is situated in the Alentejo, in the Guadiana River, being the largest artificial water reservoir in Western Europe, with an area of 250 km 2 [21]. Since 2004, it has been operating with the objective of irrigation, water supply, navigation and electricity production [22,23]. The Alqueva was selected for being a recent reservoir, having high dimensions and being located in the south of Portugal on one of the main course rivers.
Aguieira was selected because it is one of the reference reservoirs for the WFD cal and is located in the centre of Portugal on a main course river. Alqueva rese situated in the Alentejo, in the Guadiana River, being the largest artificial water r in Western Europe, with an area of 250 km 2 [21]. Since 2004, it has been operating objective of irrigation, water supply, navigation and electricity production [22, Alqueva was selected for being a recent reservoir, having high dimensions an located in the south of Portugal on one of the main course rivers. Due to the small area of Miranda and Pocinho reservoirs and the fact that th reservoir has the same characteristics and pressures, only one sampling site in these reservoirs was selected, (M and P respectively; Figure 1). In the Aguieira re four sampling sites (Ag1 to Ag4; Figure 1), and in Alqueva reservoir, five sampl (Al1 to Al5; Figure 1) were selected. Each reservoir was sampled in 4 periods, 2018 (Aut18), spring 2019 (Spr19), autumn 2019 (Aut19) and spring 2020 (Spr20) in each sampling site, several abiotic parameters were measured sub-superficial m depth): pH, conductivity (Cond, μS/cm), temperature (Temp, °C) and dissolved (O2, mg/L and %), with a multiparameter probe (Multi 3630 IDS SET F). To exam bacterial community, a water sample of 500 mL was collected in each site. Addi for the quantification of other chemical and biological parameters, 5 L of water were collected in plastic bottles. All samples were transported to the laboratory and under dark conditions for further analyses.

Laboratory Procedures
The concentration of the following compounds was determined using chromatographic methods. Nitrate concentration (NO3 − , mg/L) was quantified b chromatography of ions [24]. Turbidity (Turb, m −1 ) was determined spectrophotometer method using the protocol of Brower et al. [25]. Total suspende (TSS, mg/L) were performed according to APHA [26]. Total nitrogen concentratio mg/L) was determined using the Kjeldahl nitrogen method [27]. Total pho concentration (Ptotal, mg/L) was determined by the application of inductively plasma mass spectrometry [28]. Ammonium concentration (NH4 + , mg/L) was det  Due to the small area of Miranda and Pocinho reservoirs and the fact that the entire reservoir has the same characteristics and pressures, only one sampling site in each of these reservoirs was selected, (M and P respectively; Figure 1). In the Aguieira reservoir, four sampling sites (Ag1 to Ag4; Figure 1), and in Alqueva reservoir, five sampling sites (Al1 to Al5; Figure 1) were selected. Each reservoir was sampled in 4 periods, autumn 2018 (Aut18), spring 2019 (Spr19), autumn 2019 (Aut19) and spring 2020 (Spr20). In situ in each sampling site, several abiotic parameters were measured sub-superficially (<0.50 m depth): pH, conductivity (Cond, µS/cm), temperature (Temp, • C) and dissolved oxygen (O 2 , mg/L and %), with a multiparameter probe (Multi 3630 IDS SET F). To examine the bacterial community, a water sample of 500 mL was collected in each site. Additionally, for the quantification of other chemical and biological parameters, 5 L of water sample were collected in plastic bottles. All samples were transported to the laboratory at 4 • C and under dark conditions for further analyses.

Laboratory Procedures
The concentration of the following compounds was determined using specific chromatographic methods. Nitrate concentration (NO 3 − , mg/L) was quantified by liquid chromatography of ions [24]. Turbidity (Turb, m −1 ) was determined in a spectrophotometer method using the protocol of Brower et al. [25]. Total suspended solids (TSS, mg/L) were performed according to APHA [26]. Total nitrogen concentration (N total , mg/L) was determined using the Kjeldahl nitrogen method [27]. Total phosphorus concentration (P total , mg/L) was determined by the application of inductively coupled plasma mass spectrometry [28]. Ammonium concentration (NH 4 + , mg/L) was determined by spectrophotometric and turbidimetric analyses [29]. The determination of the content in the chlorophyll a (Chl a, µg/L) was followed according to the protocol of the Agência Portuguesa do Ambiente [30] and Instituto da Água [31], and three replicates were performed for each site.

DNA Extraction and PCR Amplification
In the laboratory, the water samples collected in a Schott glass flask was fully filtered with a 0.22 µm polycarbonate sterile filter. Half of the filter was used for the extraction of the metagenomic microbial DNA. DNA extraction was performed using the commercial kit DNeasy ® PowerSoil ® Kit (QIAGEN Group, Hilden, Germany) according to the manufacturer's instructions. The extracted DNA was used for PCR amplification of the 16S rRNA gene, using the universal primers GC-358F (5 -CCT ACG GGA GGC AGC AG-3 ; with a GC-clamp at the 5 end) and 907r (5 -CCG TCA ATT CMT TTG AGT TT-3 ) in 50 µL of PCR mixture (25 µL of NZYTaq 2× Green Master Mix; 0.5 µM of each primer and 10 µL of DNA as template). The PCR reaction was performed in a MyCycler™ Thermo Cycler (Bio-Rad, Hercules, CA, USA), with an initial denaturation step for 5 min at 94 • C; followed by 10 cycles of 1 min at 94 • C, 1 min at decreasing temperature (each cycle starting at 65 • C and ending at 55 • C) and 3 min at 72 • C; subsequent 20 cycles of 1 min at 94 • C, 1 min at 55 • C and 3 min at 72 • C; and a final extension of 10 min at 72 • C.

DGGE
The products of PCR amplification were loaded on an 8% polyacrylamide DGGE gel, using a denaturing gradient of 30-70% 7 M urea/40% formamide. DGGE electrophoresis was performed in a DCode™ universal mutation detection system (Bio-Rad Laboratories Lda, Amadora, Portugal) at 60 • C for 16 h in 1× TAE buffer with a constant voltage of 65 volts. The gels were stained for 1 h in 1× TAE buffer with SYBR ® Gold Nucleic Acid Gel Stain (Thermo Fisher Scientific, Waltham, MA, USA) and subsequently visualized in a Gel Doc EZ system (Bio-Rad Laboratories Lda, Amadora, Portugal) with Image Lab Software v4.0.1 (Bio-Rad, Hercules, CA, USA) using UV light. DGGE gels were scanned and analysed with the QuantityOne software v4.6.9 (Bio-Rad, Hercules, CA, USA) for further statistical analysis.

DNA Sequencing of Excised DGGE Bands
Thirty bands of different sizes and positions were excised from the DGGE gel according to the following criteria: (i) common to all or almost all sites; (ii) unique to a specific site; (iii) common to a specific reservoir or reservoirs; and (iv) seasonal presence. Each band was eluted in 20 mL ultrapure water prior to DNA re-amplification using the same oligonucleotide primers as above, but without the GC clamp. Amplicons were purified using the Illustra GFX PCR DNA and Gel Band Purification Kit (GE Healthcare, Chicago, IL, USA) and directly sequenced by GATC Biotech, to identify microbiota OTUs (Operational Taxonomic Units). The sequences were assembled and corrected using GENEIOUS PRIME (version 2020.0.5). Alignment of all consensus sequences was performed using MEGA 7 (Molecular Evolutionary Genetics Analysis) software, which permits to infer over time the molecular evolution between genes, genomes and species [31].

Statistical Analysis
Bacterioplankton DGGE banding patterns were transformed into presence/absence matrices and band intensities were measured using QuantityOne 1-D Analysis Software v4.6.9 (Bio-Rad Laboratories Lda., Amadora, Portugal). Relative similarities between samples and replicates were calculated using the Primer v7.0.11 (PRIMER-E, Ivybridge, UK). Besides the total number of OTUs (operational taxonomic units), the Shannon-Wiener index was used to estimate microbial diversity, while the Margalef index was used to estimate the microbial richness of each sample. To test the existence of relations between the environmental water parameters and communities variations, a Pearson correlation coefficient was calculated. The level of significance considered was α = 0.05.
For each hydrographic basin, clustering patterns were achieved by the construction of complete linkage dendrograms, based on the Bray-Curtis coefficient, using Primer software v7.0.11. The construction of the phylogenetic tree from the DGGE bands DNA sequences was performed using the calculation method maximum likelihood-ML in MEGA 7, applying general time reversible model and gamma distributed with invariant sites (G + I). The sequences were compared to the GenBank database using the Basic Local Alignment Search Tool (BLAST). Table 1 presents the results of physical, chemical and biological (Chl a as a proxy of phytoplankton biomass) parameters measured for each site over the sampling period. According to the physical and chemical parameters used in the WFD, for heavily modified and artificial water bodies, only the pH, O 2 , NO 3 − and P total have environmental quality standard (EQS) values established for a good ecological potential (EP) ( Table 1). In general, the pH presented higher values in the spring samples for all reservoirs. Similarly, the O 2 showed also higher values in the spring samples in all reservoirs. According to the WFD, the O 2 concentration was below the minimum EQS limit in M, Ag1 and Ag2 in Aut19. The values of NO 3 − did not show any pattern in Miranda, Pocinho and Alqueva reservoirs, and in Aguieira, the highest values were generally recorded in autumn samplings. In all reservoirs, the NO 3 − concentrations were always below the maximum established value. P total presented the highest concentrations, namely in the upstream sites of Aguieira and Alqueva reservoirs (Ag3, Ag4 and Al5). High concentrations of P total were also observed in Miranda and Pocinho in Aut18 and Spr20. This parameter was found several times above the EQS, namely in Miranda and Pocinho in Aut18 and Spr20; in Aguieira during all sampling periods at site Ag3; and in Alqueva at Al4_Aut18 and in site Al5 in Aut18, Spr19 and Spr20.

Environmental Variables and Nutrient Concentrations
Regarding the remaining physical and chemical parameters that have no reference value (Table 1), the Cond presented the highest values in Alqueva reservoir (>500 µS/cm), followed by Miranda and Pocinho reservoirs (>250 µS/cm) and Aguieira reservoir presented the lowest values observed (almost always <100 µS/cm). Turb and TSS appear to be related, as sites with higher turbidity also have higher TSS content (e.g., Ag3_Aut18 and Al5 in all seasons). The water temperature fluctuated considerably throughout the sampling periods (9 to 33 • C), revealing seasonality in all reservoirs. N total and NH 4 + showed several values below the detection level during the sampling periods; however, when they were quantifiable, they presented higher concentrations in the most upstream sites of the reservoirs (Ag3, Ag4 and Al5).
Concerning the analysed biological parameter, Chl a, was mostly above the permitted limit (Table 1). For Miranda, Pocinho and Alqueva reservoirs, values were found above the limit for good EP (EQS ≤ 9.50 µg/L) in all spring samples of Miranda and Pocinho and in Al5 in all sampling periods. Regarding the Aguieira reservoir, all sites in spring samples, Ag2 and Ag3 of Aut18 and in Ag3 of Aut19 presented values above the limit (EQS ≤ 7.90 µg/L). According to the previously analysed parameters, the water quality of reservoirs tends to show low water quality, namely in Miranda and Pocinho, almost all sites of Aguieira and in the Al5 site of Alqueva. Table 1. Results of physical, chemical and biological parameters: pH, conductivity (Cond), temperature (Temp), dissolved oxygen (O 2 ), turbidity (Turb), total suspended solids (TSS), nitrates (NO 3 − ), total Kjeldahl nitrogen (N total ), total phosphorus (P total ), ammonium (NH 4 + ) and chlorophyll a (Chl a). Sampling sites: Miranda-M, Pocinho-P, Aguieira-Ag1 to Ag4 and Alqueva-Al1 to Al5; sampling periods: autumn of 2018 (Aut18), spring of 2019 (Spr19), autumn of 2019 (Aut19) and spring of 2020 (Spr20). The bold values represent the values outside the environmental quality standards. * stands for physical and chemical EQS of north type. ** stands for physical and chemical EQS of south type; (1) [30].

Relations between Bacterial Communities and Environmental Parameters
Correlations between environmental parameters and the bacterioplankton communities are presented in Table 3. Miranda and Pocinho reservoirs did not show correlations between the analysed data. Regarding Aguieira reservoir, the bacterial community showed    (Table 1). Within the autumn samples, it was also possible to verify the separation of samples by reservoir (M/P); however, in the spring samples, this separation did not occur.  Figure 3 represents (a) the Bray-Curtis dendrogram and (b) the profile of the DGGE bands for each sampling site of Aguieira reservoir. The analysis of bacterioplankton for this reservoir revealed different groups separated by the season (Aut/Spr) with the exception of Ag3_Aut19, which was separated from all groups. In addition, it was possible to observe the greater similarity between sites that are geographically close to the dam (Ag1 and Ag2) as well as for sites that are geographically furthest from the dam (Ag3 and  of Ag3_Aut19, which was separated from all groups. In addition, it was possible to observe the greater similarity between sites that are geographically close to the dam (Ag1 and Ag2) as well as for sites that are geographically furthest from the dam (Ag3 and Ag4).  Figure 3 represents (a) the Bray-Curtis dendrogram and (b) the profile of the DGGE bands for each sampling site of Aguieira reservoir. The analysis of bacterioplankton for this reservoir revealed different groups separated by the season (Aut/Spr) with the exception of Ag3_Aut19, which was separated from all groups. In addition, it was possible to observe the greater similarity between sites that are geographically close to the dam (Ag1 and Ag2) as well as for sites that are geographically furthest from the dam (Ag3 and Ag4).

Microbial Community Phylogeny
Of the 30 DGGE excised bands, 24 bands were phylogenetically identified based on the analysis of the 16S rRNA gene. These belonged to five different bacterial phyla: Actinobacteria, Bacteroidetes, Cyanobacteria, Nitrospirae and Proteobacteria (Table S1 and Figure  5). Although the PCR was performed with bacterial specific primers, eukaryotic DNA was also detected, which belong to the Cyclopidae family (Macrocyclops albidus) ( Table S1). The 24 different excised bands identified seventeen different strains. In four of the bands, it

Microbial Community Phylogeny
Of the 30 DGGE excised bands, 24 bands were phylogenetically identified based on the analysis of the 16S rRNA gene. These belonged to five different bacterial phyla: Acti-nobacteria, Bacteroidetes, Cyanobacteria, Nitrospirae and Proteobacteria (Table S1 and Figure 5). Although the PCR was performed with bacterial specific primers, eukaryotic DNA was also detected, which belong to the Cyclopidae family (Macrocyclops albidus) ( Table S1). The 24 different excised bands identified seventeen different strains. In four of the bands, it was not possible to make an identification beyond the phylum.

Microbial Community Phylogeny
Of the 30 DGGE excised bands, 24 bands were phylogenetically identified based on the analysis of the 16S rRNA gene. These belonged to five different bacterial phyla: Actinobacteria, Bacteroidetes, Cyanobacteria, Nitrospirae and Proteobacteria (Table S1 and Figure  5). Although the PCR was performed with bacterial specific primers, eukaryotic DNA was also detected, which belong to the Cyclopidae family (Macrocyclops albidus) ( Table S1). The 24 different excised bands identified seventeen different strains. In four of the bands, it was not possible to make an identification beyond the phylum.  From Actinobacteria phylum, Actinobacteria bacterium IMCC19121 was identified in Miranda reservoir only in Spr19, in Pocinho reservoir in both autumns and in all sites of Aguieira only in Aut18. Actinobacterium SCGC AAA043-I11 was identified in all autumn samples from Pocinho and Aguieira reservoirs and in almost all sites of autumn samples from Alqueva reservoir. Candidatus Nanopelagicus limnes MMS-21-122 was found in all the reservoirs and in all seasons. Candidatus Planktophila vernalis MMS-IIA-15 was identified in Miranda and Pocinho except in the Spr19 and in the Aguieira in both spring samplings. The Alphaproteobacteria, Candidatus Fonsibacter ubiquis LSUCC0530 was identified in all reservoirs and seasons except in Aguieira in Spr19. The Bacteroidetes, Bacteroidetes bacterium SCGC AAA023-M10 and Bacteroidetes bacterium zo35 were only retrieved in Alqueva reservoir, except in Spr19 for the first species and in both autumns for the second species. The Pedobacter sp. strain AR-3-17 was found in Miranda and Pocinho reservoir, only in the Aut18. In the Betaproteobacteria, Curvibacter sp. W2.09-301r was only identified in some sites of autumn samples (Miranda and Pocinho of Aut18 as well as in Ag1_Aut18 and Ag3_Aut19). Polynucleobacter sp. strain LSUCC0497 was identified in Ag1 and Ag2 of Aut19 and Spr20 and in Pocinho of Aut18 and Spr19. Vogesella perlucida strain WS10 and Vogesella sp. strain UTW7-1 was only observed in Al1_Spr19. The Cyanobacteria, Calothrix elsteri strain CCALA 953 was observed in almost all sites of Alqueva reservoir. The other cyanobacteria identified in this study, Microcystis aeruginosa FD4, was only identified in Ag3_Aut18. The Gammaprotobacteria, Acinetobacter sp. St01 was only observed in spring samples from Pocinho and Aguieira. Arenimonas maotaiensis strain YT8 was only identified in Al5 of autumn samples and in Al4_Spr19. The Nitropirae, Nitrospira japonica strain NJ11 was only found in Ag1 and Ag2 of Aut18.

Discussion
The water quality of the reservoirs observed in the present study are mainly due to the high concentrations of P total and Chl a. Over the last years, all the studied reservoirs were characterized as eutrophic, with high P total [17,[32][33][34] and Chl a concentrations [33,35]. Recently, Pinto et al. [36] also described high Chl a concentration in Aguieira reservoir. In the present study, the authors observed a bloom of Cyanobacteria in Ag3 in the autumn of 2018. The study area presents constant exposure to fertilizers from agricultural runoff and from urban and livestock effluent discharges, promoting the accumulation of nutrients such as P [21,37]. Regarding O 2 levels for all reservoirs, Bordalo et al. [17] and SNIRH: Sistema Nacional de Informação de Recursos Hídricos [38] reported occasionally low O 2 levels, with 2 mg/L being the lowest concentration observed in the Miranda reservoir and 3 mg/L in Pocinho, Aguieira and Alqueva reservoirs, values in accordance with the results obtained in the present study.
The bacterioplankton showed sensitivity to the characteristics and changes that occurred in the studied reservoirs. Luo et al. [39] described that the bacterial community of reservoirs has a strong spatial and temporal dynamics but referred that this community is scarcely studied in these ecosystems. In addition, Tamaki et al. [40], in a study of bacterial diversity in freshwater sediment of a shallow eutrophic lake, also referred to the existence of spatial and temporal dynamics and a close relation of the bacterioplankton to the water quality and possible presence of pathogens species. In our study, some of this evidence could also be observed, such as bacterial temporal dynamics in Miranda, Pocinho and Aguieira, bacterial spatial dynamics in Aguieira, and the presence of Acinetobacter sp. St01, a species related to membrane bioreactors of wastewater treatment.
With the exception of Nitrospirae, the bacterial phyla identified in the studied reservoirs correspond to the most common phyla recovered from lake communities: Actinobacteria, Proteobacteria, Bacteroidetes and Cyanobacteria. Accordingly, the taxonomic taxa found have been previously described in freshwater environments [41]. The most abundant taxa were Proteobacteria (7 distinct taxa) followed by Actinobacteria (four distinct taxa), three different Bacteroidetes, two Cyanobacteria and one Nitrospirae.
The phylum Actinobacteria is known to constitute 10% to 60% of the bacterial community in freshwater habitats [42][43][44][45], with acI clade being the most abundant microbes in freshwater systems [41,46]. These bacteria are ubiquitous in these environments and are able to live in different ranges of pH, oxygen and salinity [47]. In accordance with the literature, we retrieved genetic material of acI clade (Candidatus Nanopelagicus limnes isolate MMS-21-122, Actinobacteria bacterium IMCC19121, Candidatus Planktophila vernalis isolate MMS-IIA-15) from all the reservoirs and seasons, reinforcing the expressive presence of this clade in different ecosystems. acI clade has been described as to suffer small changes with seasons, being the most constant freshwater bacteria [41]. This ubiquity can be explained by alternative energy metabolisms through rhodopsins, small cell size, cell wall composition and/or UV stress resistance. Genetic material of actinobacterial strains was found in all reservoirs. However, Actinobacteria bacterium IMCC19121 and Actinobacterium SCGC AAA043-I11 were observed only in autumn samplings, with the exception of M_Spr20. Previous studies already described that pH and carbon sources could explain fluctuations in Actinobacteria freshwater populations [41]. As physical and chemical parameters of the studied reservoirs suffer slight changes through the seasons and locations, Actinobacteria communities are probably not affected by those.
Regarding classes within the Proteobacteria observed in these reservoirs, Betaproteobacteria prevailed, followed by Alphaproteobacteria and Gamaproteobacteria, as already described in other studies [44,48,49]. Betaproteobacteria is the most studied bacterial class in freshwater lakes, and it is thought to be selected by depth, pH, carbon availability and seasons [41]. In this study, it was not possible to establish any correlation between physical and chemical factors and the identified betaproteobacteria. Contrary to Actinobacteria or Proteobacteria, Bacteroidetes were present only in some samples (e.g., M_Aut18, Ag1 in almost all samples and Al1_Spr19), but no seasonal or biographic pattern could be established.
The Alphaproteobacteria Candidatus Fonsibacter ubiquis strain LSUCC0530 was present in all reservoirs in almost all samples. This bacterium belongs to the LD12 subclade within the order Pelagibacterales (SAR11), which has evolved as a unique freshwater lineage and is the most dominant freshwater bacterioplankton [50].
The bacterioplankton community of Miranda and Pocinho reservoirs showed temporal and spatial dynamics, but no significant correlation with the physical and chemical parameters analysed could be established. Furthermore, in the spring samples, where no spatial pattern existed, the sites grouped according to water ecological quality (Miranda and Pocinho in Spr20 as moderate and in Spr19 as good quality), which indicates the possibility of characterizing these waterbodies. In spring samples, we observed the presence of Acinetobacter sp. St01, a species related to membrane bioreactors of wastewater treatment [51], which may be the result of non-treated discharges by wastewater treatment plants (WWTP), already described for Miranda reservoir and confirmed by the large amounts of faecal coliforms (10,000 CFU/100 mL) and total coliforms (30,000 CFU/100 mL) verified in 2014 [38]. Furthermore, according to the SNIRH: Sistema Nacional de Informação de Recursos Hídricos [38], the amount of Escherichia coli showed high values such as 500 MPN/100 mL in 2019 and 1000 MPN/100 mL in 2020, possibly indicating that this water body is polluted, with low quality. This water degradation can be problematic since this water body is also used for human consumption [52]. In the case of the Pocinho reservoir, according to the SNIRH: Sistema Nacional de Informação de Recursos Hídricos monitoring stations in this area [38], the quantification of faecal and total coliforms was always found in low quantities (<400 CFU/100 mL and <1000 CFU/100 mL, respectively, in 2014) in the last 20 years; thus apparently, the presence of Acinetobacter sp. St01 is not due to the same disturbances as observed in Miranda. However, the existence of pollution sources from wastewater is recognized [53]. Furthermore, Acinetobacter sp., is known for its function as degrader of polycyclic aromatic hydrocarbons (PAHs) [54], and following the WFD guidelines regarding these type of compounds, the PAHs analysed were not detected in the present work [18,36]. Curvibacter sp. W2.09-301r was observed in Miranda and Pocinho reservoirs in Aut18 and has the function of trichloroethylene degradation [55], a compound not included in the WFD analysis. However, it is impossible to chemically analyse all compounds that may be present in the surface waters. Conversely, Huang et al. [54] observed that this species, when in the presence of petroleum-hydrocarbon, tends to disappear or not be detected in DGGE. This decrease of its abundance in the spring samples and autumn of 2019 may reflect the existence of disturbances, which may indicate stress conditions; yet, these disturbances have not been mentioned for both reservoirs. Additionally, both reservoirs are part of several cruise routes, where in each one, there is a boat dock [56].
Regarding the Aguieira reservoir, the temporal dynamic may be due to the positive correlation observed with Temp, since warmer periods are associated with greater abundance, diversity and evenness of bacterial species [57]. Gurung et al. [57] in Lake Biwa, in Japan, attributed the positive correlation with temperature to the fact that high nutrients concentrations are not a limiting factor. This observation is consistent with our results since Aguieira reservoir presented high nutrients concentration. Conversely, researchers suggest that an increase in bacterial abundance is associated with an increase in the trophic status of the waters, namely with Chl a [58]. However, in the present work, no correlation was observed with this element, although it tends to present higher values in spring, when temperatures also tend to be higher. Similar to Miranda and Pocinho, Acinetobacter sp. St01 was also observed at all spring samples. According to the SNIRH: Sistema Nacional de Informação de Recursos Hídricos [38], there is only information on faecal and total coliforms for Ag1 and Ag4. In Ag1 the values observed are always low, whereas in Ag4, these values reach 400 CFU/100 mL for faecal coliforms and 600 CFU/100 mL for total coliforms in 2013, reflecting the low quality of this water body and the possible contamination by wastewater. According to Agência Portuguesa do Ambiente [59], this reservoir presents significant pressures, namely agricultural, livestock, urban and industrial, that can justify the presence of Curvibacter sp. W2.09-301r. Nitrospira japonica is not usually described as one of the major phyla in freshwater ecosystems; however, it is known to play an important role in nitrification and can colonize a wide range of habitats and to sustain shifts in environmental conditions such as changing in oxygen concentrations [60]. However, we do not observe any relation with NO 3 − (Table 3). Nevertheless, the Nitrospira japonica strain NJ11 [61] was first described in activated sludges of a wastewater treatment plant, and this reservoir presents urban pressure which can justify the presence of this organism. Cyanobacteria is also one of the most described bacterial phyla in freshwater ecosystems. These bacteria are well known for their ability to perform oxygenic photosynthesis, relying only on Chl a [62]. We were able to identify Microcystis aeruginosa, but only in Ag3_Aut18, although this organism has already been reported for the entire reservoir [34]. Ag3 is one of the sampling sites with a higher value of Chl a (Table 1), indicating a possible Cyanobacteria bloom. Pinto et al. [36] had already reported, for this site, the existence of a bloom of cyanobacteria, namely of Microcystis sp. The presence of this organism and the existence of blooms reflect the poor quality of the water body, namely the enrichment of nutrients [63]. Furthermore, this microorganism is known for the production of toxic toxins that negatively affect the entire dynamics of the ecosystem [64]. Therefore, the fact that this reservoir is used to capture water for consumption and for recreational activities [59] is worrying as there are species that are dangerous to humans health.
Gannon et al. [65] observed that the high-size, low-rate of velocity as well as the high residence times of a body water promote bacterial die-off as well as bacterial sedimentation. These characteristics can be framed in the Alqueva reservoir, where we can observe the same results regarding the bacterioplankton community. He et al. [66] reported that high levels of conductivity promote high salt concentrations, which can inhibit bacterial growth or damage microorganisms, which is also consistent with our results. The bacterial community in Alqueva tends to be different from the other reservoirs. Bacteroidetes bacterium SCGC AAA023-M10, Bacteroidetes bacterium zo35, Vogesella perlucida strain WS10, Vogesella sp. strain UTW7-1 and Arenimonas maotaiensis strain YT8 were observed only and occasionally in this reservoir, and little is known about their functions in the ecosystem.

Conclusions
In conclusion, from a general point of view, the bacterioplankton community analysed in the present study showed local sensitivity. The presence or absence as well as the existence of seasonality of some organisms may indicate different degrees of water quality as well as punctual and local disturbances. However, this study represents only an initial and general approach to the biological parameter analysed. Using the DGGE in the approach to assess reservoir water bodies allows us to make an initial screening of the reservoir status and find sites with potential disturbance. Furthermore, it is a standardized methodology, of quick and easy performance, and is economical, which would allow the execution of this approach in a larger number of water samples and their comparison. An important aspect missing from this study would be information about the total structure of the bacterial community in each site. Future work with different approaches (e.g., number of sites, number of samples, different laboratorial methodologies) as well as highthroughput analyses are important to determine how this microbial community is present and active in these ecosystems. Furthermore, this work showed us the existence of many organisms that are not yet properly characterized, which makes it difficult to correlate with the reservoirs environmental characteristics and water quality.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/w13202836/s1, Table S1: Presence or absence for all phylogenetically identified species for all sampling sites and periods. Funding: This research was funded by National Funds (through the FCT-Foundation for Science and Technology) and by the European Regional Development Fund (through COMPETE2020 and PT2020) through the research project ReDEFine (POCI-01-0145-FEDER-029368) and the strategic program UIDB/04423/2020 and UIDP/04423/2020. Sara Antunes was hired through the Regulamento do Emprego Científico e Tecnológico-RJEC from the Portuguese Foundation for Science and Technology (FCT) program (CEECIND/01756/2017).