Iame A. Guedes1
Caio T. C. C. Rachid2
Luciana M. Rangel3
Lúcia H. S. Silva3
Paulo M. Bisch1
Sandra M. F. O. Azevedo1
Ana B. F. Pacheco1*- 1Instituto de Biofísica Carlos Chagas Filho, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil
- 2Instituto de Microbiologia Professor Paulo de Góes, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil
- 3Departamento de Botânica, Museu Nacional, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil
Cyanobacteria tend to become the dominant phytoplankton component in eutrophic freshwater environments during warmer seasons. However, general observations of cyanobacterial adaptive advantages in these circumstances are insufficient to explain the prevalence of one species over another in a bloom period, which may be related to particular strategies and interactions with other components of the plankton community. In this study, we present an integrative view of a mixed cyanobacterial bloom occurring during a warm, rainy period in a tropical hydropower reservoir. We used high-throughput sequencing to follow temporal shifts in the dominance of cyanobacterial genera and shifts in the associated heterotrophic bacteria community. The bloom occurred during late spring-summer and included two distinct periods. The first period corresponded to Microcystis aeruginosa complex (MAC) dominance with a contribution from Dolichospermum circinale; this pattern coincided with high water retention time and low transparency. The second period corresponded to Cylindrospermopsis raciborskii and Synechococcus spp. dominance, and the reservoir presented lower water retention time and higher water transparency. The major bacterial phyla were primarily Cyanobacteria and Proteobacteria, followed by Actinobacteria, Bacteroidetes, Verrucomicrobia, and Planctomycetes. Temporal shifts in the dominance of cyanobacterial genera were not only associated with physical features of the water but also with shifts in the associated heterotrophic bacteria. The MAC bloom was associated with a high abundance of Bacteroidetes, particularly Cytophagales. In the second bloom period, Planctomycetes increased in relative abundance, five Planctomycetes OTUs were positively correlated with Synechococcus or C. raciborskii OTUs. Our results suggest specific interactions of the main cyanobacterial genera with certain groups of the heterotrophic bacterial community. Thus, considering biotic interactions may lead to a better understanding of the shifts in cyanobacterial dominance.
Introduction
Cyanobacterial blooms occur in freshwater environments around the world, mainly as a result of eutrophication (Rigosi et al., 2014). These events cause deleterious environmental and socioeconomic effects, impacting the ecosystem due to disruption of the food webs and thus potentially decreasing biodiversity. Cyanobacterial blooms also impair the use of water by human populations (Paerl and Paul, 2012). During the last few decades, an expansion of cyanobacterial blooms has been recorded, and global warming is expected to further exacerbate this situation (Paerl, 2008; Visser et al., 2016).
Potentially harmful cyanobacteria tend to become the dominant phytoplankton component in eutrophic freshwater environments, especially during warmer seasons (Paerl and Huisman, 2008). Indeed, the main reported abiotic drivers of cyanobacterial blooms are increased nutrient loading (nitrogen and phosphorus) and rising temperatures (Kosten and Huszar, 2012; Lürling et al., 2017). Some bloom-forming cyanobacteria have efficient nutrient uptake and storage abilities and can use both inorganic and organic N and P pools (O'Neil et al., 2012; Paerl and Paul, 2012; Harke and Gobler, 2013). Generally, cyanobacteria will grow faster at higher temperatures compared to other phytoplankton groups, although growth rates can be differentially affected depending on the species considered (Paerl, 2008; Paerl and Paul, 2012; Lürling et al., 2013). Rising temperatures also intensify vertical stratification, and in this situation some bloom-forming cyanobacteria can be benefited due to their buoyancy ability (Kruk et al., 2010). Cells can accumulate at the surface and shade underlying layers, thus outcompeting other phytoplankton groups through light limitation (Reynolds, 2006; Harke et al., 2016). Additionally, buoyant cyanobacteria can access nutrients from deeper waters when epilimnion concentrations are diminished (Graham et al., 2016). These general observations, however, are insufficient to explain the prevalence of one species over another in a bloom period, which may be related to particular adaptive strategies and interactions with other components of the plankton community (Woodhouse et al., 2016; Wood et al., 2017).
In contrast to abiotic factors, the role of biotic interactions on cyanobacterial bloom dynamics has been less explored. Diverse biotic factors such as grazing, predation, parasitism and mutualism influence cyanobacterial biomass through interactions with other plankton components such as protozoans, zooplankton, bacteria, and viruses (Paerl and Otten, 2013; Gerphagnon et al., 2015; Ger et al., 2016; Steffen et al., 2017). Interactions between cyanobacteria and heterotrophic bacteria can be positive due to the exchange of nutrients and oxygen, which can benefit both microorganisms (Bagatini et al., 2014; Gerphagnon et al., 2015). Interactions may also be negative, as for cyanolytic bacteria (Van Wichelen et al., 2016; Osman et al., 2017). During cyanobacterial blooms, bacteria can be found directly attached to cyanobacterial cells or adjacent to them, occupying the cyanobacterial phycosphere (Louati et al., 2015). In many cases, these associations were reported under laboratory conditions, but several recent studies in natural environments have provided evidence of the close relation between cyanobacterial biomass and associated heterotrophic bacteria (Eiler and Bertilsson, 2004; Wu et al., 2007; Kormas et al., 2010; Cheng et al., 2011; Dziallas and Grossart, 2011; Li et al., 2011, 2015; Wilhelm et al., 2011; Steffen et al., 2012; Cai et al., 2013; Parveen et al., 2013; Louati et al., 2015; Woodhouse et al., 2016; Parulekar et al., 2017; Salmaso et al., 2017). Moreover, specific associations between some cyanobacterial genera and heterotrophic bacteria have recently been reported (Bagatini et al., 2014; Louati et al., 2015), pointing to a possible connection of those bacteria attached to cyanobacteria and their participation in cyanobacterial bloom dynamics. Some authors extend this view and suggest that microbial communities of distinct taxonomic composition can play similar functional roles in bloom events (Steffen et al., 2012). Clearly, the complexity of these microbial interactions is still little explored, and their impact in the ecophysiology of cyanobacteria is underestimated.
High-throughput sequencing has revolutionized microbial ecology in recent years, particularly through the characterization of metagenomes. By following temporal variation in the composition of a community, it is possible to reveal correlations and to infer interactions among taxa (Stubbendieck et al., 2016). This application can be explored to better understand cyanobacterial species dynamics during a bloom, and the accompanying bacterial community can be seen as an integral and essential part of these events.
In this study, we present an integrative view of a mixed cyanobacterial bloom occurring during a warm, rainy period in a tropical hydropower reservoir; we used high-throughput sequencing to follow temporal shifts in the dominance of cyanobacterial genera and the associated heterotrophic bacterial community. We also looked for correlations between temporal shifts in the bacterial community and abiotic factors.
Materials and Methods
Sampling
Water samples were collected in the Funil Reservoir (22°30′S, 44°45′W), Rio de Janeiro, Brazil (Supplementary Figure 1), from October 2013 to March 2014. This is a eutrophic reservoir with an area of 40 km2, a volume of 8.9 × 106 m3, maximum and medium depth of 70 and 22 m, respectively, and average residence time of 41.5 days (Soares et al., 2009; Rangel et al., 2012). The region in which the reservoir is located typically presents warm, rainy conditions during summer and cold, dry conditions during winter. The sampling period encompassed the end of spring and the summer.
Samples were collected from two locations: one in the central part of the reservoir (point 1) and the other near the dam (point 2). From one to three samples per month were obtained from the integrated euphotic zone (determined as 2.7 times the Secchi disk depth; Cole, 1994) for a total of 22 samples. Water temperature and pH were measured using a Yellow Spring multiparametric probe (model 600 QS), and the water transparency was determined using a Secchi disk. Water temperature and pH values were measured at 0.5-m intervals to the end of the euphotic zone, and average values are presented. Volumes of 0.3–1 L of water (depending on the phytoplankton density) from the integrated euphotic zone were filtered (Whatman GF/F, 0.7 μm) to collect cells and then stored at −20°C for DNA extraction.
Nutrient Analysis
Aliquots of water for total and dissolved nutrient analysis were collected in polypropylene tubes and stored at −20°C. For dissolved nutrient analysis, aliquots were filtered (Whatman GF/F, 0.7 μm) before storage. The soluble reactive phosphorus (SRP), ammonium, nitrate, nitrite, and total phosphorus (TP) were measured using flow injection analysis (FIAlab 2500) according to the manufacturer's instructions (FIALab Instruments Inc., Seattle, Washington).
Quantitative Analysis of Phytoplankton
Aliquots of the integrated water samples were stored in amber glass vials with Lugol's solution. Phytoplankton abundance was determined by the Utermöhl method (Utermöhl, 1958) using an inverted optical microscope (Olympus BX-51). The biovolume (mm3 L−1) was estimated by multiplying the density of each species by the average volume of its cells (Hillebrand et al., 1999).
DNA Extraction and 16SrDNA Amplification
DNA was extracted from cells collected in filters (Whatman GF/F, 0.7 μm) using the Power Soil DNA Isolation Kit (Mo Bio) according to the manufacturer's instructions. DNA samples were quantified using a fluorimeter (Qubit, Thermo Fisher Scientific). Amplification of the v3-v4 region of 16S rDNA genes was performed with the primers S-D-Bact-0341-b-S-17F (5′-CCTACGGGNGGCWGCAG-3′) and S-D-Bact-0785-a-A-21R (5′-GACTACHVGGGTATCTAATCC-3′) (Klindworth et al., 2013) containing the appropriate adaptors for sequencing in the Illumina platform. Amplifications were performed in a 25 μl reaction mixture containing 12.5 μL of HiFi HotStart ReadyMix (KAPA Biosystems), 0.2 μM of each primer and 12.5 μg of DNA. The PCR program included an initial denaturation at 95°C for 3 min followed by 25 cycles of amplification (95°C for 30 s, 55°C for 30 s, and 72°C for 30 s) and a final step of 72°C for 5 min. Products were purified using magnetic beads (Agencourt AMPure XP, Beckman Coulter) and subjected to a second PCR to incorporate dual indices, as described in the 16S Metagenomic Sequencing Library Preparation Protocol for the Illumina MiSeq System. The size and quality of DNA in the final libraries were verified on a Bioanalyzer 2100 (Agilent) using a Bioanalyzer DNA 1000 chip (Agilent Technologies). After quantitative PCR with a KAPA Library Quantification Kit for Illumina (KAPA Biosystems), samples were normalized and pooled for sequencing.
DNA Sequencing and Data Analysis
Sequencing was performed in a MiSeq platform (Illumina) using the MiSeq Reagent Kit v3 (2 × 300 base pairs) according to the manufacturer's instructions. Files were recovered (.fastq), and paired-end reads were joined using mothur v.1.35.1 (Schloss et al., 2009). Sequences are available for download via the NCBI short read archive under BioProject PRJNA406945. The following criteria were used to eliminate low-quality reads: average quality (window size = 50) <30, length 460 base pairs, presence of ambiguous characters (“N”), homopolymer <8. The remaining reads were aligned using the SILVA database, trimmed and filtered. Then, sequences were preclustered with diff = 4. Chimeras were detected using UCHIME (Edgar et al., 2011) and excluded. Taxonomic classification was carried out using the RDP database (Release 11) with a confidence threshold of 80%. Sequences not assigned as Bacteria or classified as Chloroplast or Mitochondria were discharged. Singletons and doubletons were removed, and the number of sequences in the 22 samples was normalized to the same number of sequences. Sequences were then used as input to generate a distance matrix and clustered into operational taxonomic units (OTUs) at the sequence similarity cutoff of 97%. Species richness and Shannon diversity index were calculated in mothur. Taxonomic assignment of OTUs was performed using Greengenes (version 13_5).
The OTU relative abundances in the samples were used to generate an ordination plot by nMDS (non-metric multidimensional scaling) based on Bray-Curtis similarity coefficients. The limnological parameters and cyanobacterial cell counts with significant differences between the periods (t-test, p < 0.001) were plotted together with the nMDS. Statistical significance of sample grouping was tested with PERMANOVA. These analyses were performed on package Past3 (Hammer et al., 2001). To identify the major OTU contributors to grouping differentiation (periods), we used a similarity percentage analysis (SIMPER) (Clarke, 1993). Spearman correlation was used to test the degree of association among limnological variables with cyanobacterial cell counts and cyanobacterial OTUs and to test the association between specific cyanobacterial OTUs and heterotrophic bacterial OTUs (considering only those that contributed at least 0.2% to the total of sequences). We considered relevant correlations those with p < 0.001 and r > 0.6. The visualization of these correlations was made with the Cytoscape package version 3.4.0 (available at: www.cytoscape.org) using the plugin CoNet. The r-value (>0.6) was selected to support the generation of an edge-weighted spring-embedded network (Assenov et al., 2008). The network included the classified OTUs which contributed with at least 0.2% of the total.
Results
Limnological Characterization
During the study period, Funil Reservoir was characterized by turbid, dynamic, and slightly eutrophic conditions (Table 1). Water transparency varied between 0.25 and 2 m, with lower values between October and early January. Water retention time ranged between 23 and 38 days, with lower values beginning in mid-January, coinciding with the increase in precipitation. The pH was generally alkaline, ranging between 7.1 and 10.9, with highest values on January and February. Nitrogen concentrations were always above the potentially limiting concentration for phytoplankton growth (>100 μg L−1), while SRP concentrations were below the detection limit (<2 μg L−1) in almost all samples, suggesting strong phosphorus limitation during most of the studied period (Table 1).
Table 1. Limnological variables and abiotic parameters associated with the collected samples.
Cyanobacterial Dynamics
Microscopic analysis revealed that cyanobacteria were the dominant (96.2–99.8%) component of phytoplankton all over the sampling period; 22 cyanobacterial species were identified (Supplementary Table 1). The most abundant species were Microcystis spp. (referred herein as Microcystis aeruginosa complex, MAC), Dolichospermum circinale, Cylindrospermopsis raciborskii, Pseudanabaena mucicola, Synechococcus nidulans, and Synechocystis aquatilis (Figure 1). An MAC bloom was apparent from October to mid-January with the co-occurrence of D. circinale and P. mucicola. Beginning in mid-January, C. raciborskii dominated the phytoplankton community together with D. circinale, S. nidulans, MAC and other cyanobacteria. Although the cyanobacterial biomass differed between the two sampling points during the period, the relative contribution of the cyanobacterial species was similar in the two locations (Figure 1).
Figure 1. Variation of cyanobacteria density from October 2013 to March 2014 in two locations of Funil reservoir (A) point 1, central part, (B) point 2, near the dam.
The cyanobacterial community composition assessed by 16S rDNA sequencing revealed a somewhat different pattern for the mixed bloom. The main cyanobacteria genera identified were Microcystis, Dolichospermum, Cylindrospermopsis, Pseudanabaena, and Synechococcus (Figure 2). Microcystis dominated from October to end of January, as also observed with microscopy, but the dominance of C. raciborskii was not as evident. Instead, a Synechococcus predominance followed the decay of the Microcystis bloom. This pattern was evident in the two sampling points despite some differences in specific dates.
Figure 2. Variation of cyanobacterial community accessed by 16S rDNA sequencing from October 2013 to March 2014 in two sampling points of Funil reservoir (A) point 1, central part, (B) point 2, near the dam.
Microbial Community Structure
The metagenomic sequencing of the bacterioplankton community resulted in a total of 2,547,075 sequences (101,883 per sample) that clustered into 6,452 OTUs (3% dissimilarity). Rarefaction curves showed an excellent coverage (Supplementary Figure 2) for all samples. The bacterial richness varied from 613 to 2,038 OTUs, and Shannon indices varied from 3.2 to 5.0 (Table 1).
Analysis of the whole microbial community showed that Cyanobacteria and Proteobacteria were the dominant phyla followed by Actinobacteria, Bacteroidetes, Verrucomicrobia, and Planctomycetes (Figure 3). The relationship between all samples was assessed using an NMDS based on the bacterial community distribution at OTUs level (3% dissimilarity). The ordination showed clustering of the samples according to the two periods. The difference between these periods was supported by PERMANOVA analysis (F = 3.19, p = 0.0001). Period 1 corresponded to samples from October to mid-January (Microcystis dominance), and period 2 from mid-January to March (Synechococcus dominance; Figure 4). For spatial variability, there was no significant difference between the sampling points (F = 1.13, P = 0.29). The distinction between these two periods was also supported by some limnological parameters and cyanobacterial microscopic counts. Water retention time and transparency showed significant differences between periods one and two [t(20) = 2.8, p < 0.01 and t(20) = −4.9, p < 0.01, respectively]. Microcystis spp. and C. raciborskii cell counts were significantly different between periods [t(16) = 3.2, p < 0.01 and t(16) = −3.29, p < 0.01, respectively; Figure 4].
Figure 3. Relative abundance of OTUs classified as Order and Phyla across the sampling period. The area of the bubbles represents the relative abundance of OTUs (average values of the two sampling points). The color of the bubbles indicates the Phylum to which the OTUs were assigned.
Figure 4. Non-metric multidimensional scaling ordination based on Bray-Curtis similarity of data from OTU abundance in the samples from the two sampling points. Squares correspond to samples from October 2013 to mid January 2014, and triangles correspond to samples from the end of January 2014 to March. Vectors are environmental variables and cyanobacterial microscopic counts that were significantly different between the two defined periods (p < 0.01). RT (retention time) WT (water transparency).
Regarding the main contributors to the distinction between the two periods, SIMPER analysis revealed that OTUs assigned as Synechococcus and Microcystis accounted for 12.75 and 10.13% of the variability, respectively. Considering other OTUs with over 1% relative abundance, 11 orders also contributed to the variability of the two periods, although their combined impact was <2% (Figure 5).
Figure 5. Average relative contribution of OTUs in the two defined periods (Period 1 from October to mid January and Period 2 from mid January to March). The selected OTUs contributed to at least 2% for the differentiation between the periods (SIMPER analysis).
Associations of Biotic and Abiotic Factors during the Bloom
The data pointed to a clear distinction of a Microcystis dominated bloom followed by a Synechococcus dominated period. The 16S rDNA sequencing and microscopy data revealed the correlations among the cyanobacterial taxa. Microcystis was positively correlated with Dolichospermum (both for microscopy counts and 16SrDNA data). C. raciborskii microscopy counts were positively correlated with Synechococcus (16S rDNA), while Pseudanabaena, which was a contributor to the transition between periods 1 and 2, was negatively correlated with Microcystis and positively with C. raciborskii (Table 2).
Table 2. Values of Spearman correlations (r) among cyanobacterial taxa considering both 16S rDNA sequencing (OTU) and microscopy data and limnological parameters (*p < 0.05, **p < 0.001).
Considering abiotic factors, Microcystis microscopic counts (MO) were positively correlated with water retention time and were negatively correlated with water transparency; C. raciborskii (MO) was negatively correlated with water transparency; Dolichospermum (MO) was negatively correlated with DIN and water transparency; Synechococcus (OTU) was positively correlated with temperature and retention time and was negatively correlated with water transparency (Table 2).
Several significant correlations were demonstrated between cyanobacterial taxa and heterotrophic bacteria taxa (Figure 6, Table 3). Regarding the dominant cyanobacteria, a few correlations, mostly negative, were observed between Microcystis and other bacterial OTUs, particularly with OTUs from Proteobacteria and Chloroflexi. Cylindrospermopsis, the least abundant among the known main Cyanobacteria, showed an opposite pattern with most of the correlations being positive, mainly with OTUs from Proteobacteria and with Planctomycetes. Synechococcus showed a higher number of interactions compared to Microcystis, and Cylindrospermopsis and those were more balanced between positive and negative correlations. Synechococcus showed numerous correlations with Proteobacteria, particularly with taxa from the family Comamonadaceae. The Dolichospermum OTU showed only one positive correlation with an OTU from Comamonadaceae. Microcystis and Cylindrospermopsis did not correlate with any other Cyanobacteria, while Synechococcus and Pseudanabaena showed positive correlation with each other.
Figure 6. Edge-weighted Spring-embedded network with significant correlations (p < 0.001 and r > 0.6) between cyanobacterial and heterotrophic bacterial OTUs. Node size is proportional to the OTU abundance and the colors indicate the Phylum to which the OTUs were assigned. Line colors are indicative of the Spearman correlation coefficient (green = positive and red = negative). The network included the OTUs which contributed with at least 0.2% of the total. Centroid sequences of the dominant cyanobacterial OTUs are listed on Supplementary File 1.
Table 3. Significant (p < 0.001) Spearman correlations (r > 0.6) among the main cyanobacterial OTUs and heterotrophic bacteria.
Discussion
In this study, we investigated the bacterial community coupled with a mixed cyanobacterial bloom (initially dominated by Microcystis and Dolichospermum, followed by C. raciborskii and Synechococcus) occurring in a tropical reservoir during late spring-summer, combining microscopy with metagenomics. Temporal shifts in the dominance of bloom-forming cyanobacterial genera were associated with physical features of the water and with shifts in the associated heterotrophic bacteria. Changes in abiotic factors during this period probably influenced both the cyanobacterial and the non-cyanobacterial microbial community. To gain insights about possible biotic associations during a cyanobacterial bloom in this tropical system we have explored specific interactions of the main cyanobacterial genera with some components of the heterotrophic bacterial community.
Cyanobacterial Dynamics
For the last several decades the Funil Reservoir has experienced cyanobacterial blooms, mainly during the summer, dominated by Microcystis, Dolichospermum, and Cylindrospermopsis (Soares et al., 2013; Guedes et al., 2014; Rangel et al., 2016), which are also the major bloom-forming genera globally (O'Neil et al., 2012; Antunes et al., 2015; Harke et al., 2016). The succession of a bloom dominated by Microcystis and Dolichospermum by the dominance of C. raciborskii during the summer has been reported for several years in this system, and it has been more tightly associated with temperature and physical variables than with nutrient availability (Soares et al., 2012). Thus, the stratification of the water column and the reservoir residence time are important abiotic factors associated with the bloom dynamics (Soares et al., 2009; Rangel et al., 2016). In the present study, the Microcystis bloom was associated with high retention time and low water transparency. Similarly, previous studies in this reservoir have associated the dominance of Microcystis with a more stable and prolonged period of thermal stratification (Soares et al., 2009, 2012; Rangel et al., 2016). M. aeruginosa is considered a species adapted to high light intensities (Robarts and Zohary, 1987); it can thus form blooms on the surface. As a consequence, Microcystis can outcompete other species by reducing available light for nonbuoyant phytoplankton competitors (Harke et al., 2016). Indeed, in the initial period of the bloom, Microcystis accounted for 89% of the total cyanobacterial density. In addition, this species is sensitive to mixing and turbulence (Reynolds, 2006), which can explain the decline of its abundance in the middle of January (Period 1), when precipitation increased and retention time decreased.
After the decay of the Microcystis bloom, C. raciborskii and Synechococcus were the dominant cyanobacteria. Opposite to the former species, the latter two species were positively correlated with lower retention time and higher water transparency. Previous investigations in this reservoir have associated C. raciborskii blooms with periods of water column mixing (Soares et al., 2009, 2012). This species is considered to tolerate both stratified and mixed conditions and can dominate in mixed periods (Berger et al., 2006; Bonilla et al., 2012; Soares et al., 2013). This is likely due to the ability of cells to photoadapt to dark and fluctuating light conditions (O'Brien et al., 2009).
The period of C. raciborskii dominance was coupled with the high relative abundance of Synechococcus (shown by 16S rDNA sequences). This has never been reported for this reservoir, since previous studies characterized the phytoplankton community by microscopy only. Picoplanktonic cells (size 0.2–2 μm) seem to be generally subquantified by the Utermöhl method, while other techniques such as fluorescence microcopy and flow cell cytometry are more efficient for their quantification (Callieri, 2007). The high relative abundance of Synechococcus in conditions of higher water transparency and low retention time is in accordance with the high light requirement reported for this genus (Reynolds et al., 2002), which together with its high growth rate can explain its rise just after the senescence of the Microcystis bloom.
Microbial Community Dynamics
In the present study, we found a clear distinction between the heterotrophic bacteria associated with the two periods of the cyanobacterial bloom. The distinction of these periods was based on limnological parameters (water retention time and transparency) and cyanobacterial community composition (Microcystis spp. and C. raciborskii/Synechococcus dominance). While the influence of the changes in abiotic environmental parameters on the composition of the bacterioplankton community cannot be disregarded, we discuss the overall community composition shifts emphasizing the links between the dominant cyanobacterial genera and the associated bacterioplankton in this reservoir.
Changes in the composition of the bacterioplankton community can be linked to phytoplankton blooms (Xing et al., 2007; Li et al., 2015; Woodhouse et al., 2016). Recent studies noted specific associations between heterotrophic bacteria taxa and cyanobacteria, which has been attributed to different dissolved organic matter qualities and quantities produced by different cyanobacteria (Bagatini et al., 2014; Louati et al., 2015). In our case, we did not observe a correlation between cyanobacterial abundance and diversity of the microbial community.
In this study, the major observed bacteria phyla were Cyanobacteria and Proteobacteria followed by Actinobacteria, Bacteroidetes, Verrucomicrobia, and Planctomycetes (Figure 3). Actinobacteria is one of the most frequent phyla in freshwater environments and is reported as highly abundant in different lakes from oligotrophic to eutrophic systems (Newton et al., 2011). This phylum has not been found in physical association with Cyanobacteria but can be a major part of the bacterioplankton community during phytoplankton blooms (Kormas et al., 2010; Newton et al., 2011; Steffen et al., 2017). On the other hand, Bacteroidetes, Verrucomicrobia, and Planctomycetes frequently reach high relative abundance during cyanobacterial blooms (Li et al., 2015; Woodhouse et al., 2016; Parulekar et al., 2017; Steffen et al., 2017). In this study, Verrucomicrobia ranged from 1 to 18% of relative abundance. Bacteria from this phylum are described as able to degrade algal polysaccharides and organic matter, and so can be favored by high cyanobacterial abundances (Bagatini et al., 2014; Woodhouse et al., 2016; Parulekar et al., 2017). The relative abundance of Bacteroidetes also increased during the bloom, and other studies have shown that some taxa in this phylum, such as Sphingobacteria and Flavobacteria, are often found in high abundance after phytoplankton bloom decay either adjacent or attached to phytoplankton (Li et al., 2011; Bagatini et al., 2014). In the present study, the most abundant Bacteroidetes taxa were Cytophagales and Saprospirales, and higher abundances were associated with the Microcystis bloom. Planctomycetes can be abundant in nutrient-enriched waters (Woebken et al., 2007), which frequently have a high abundance of phytoplankton. Other studies also reported higher densities of Planctomycetes after diatom (Morris et al., 2006) or cyanobacterial blooms (Eiler and Bertilsson, 2004), suggesting a possible association of this phylum with phytoplankton. For example, it has already been observed that Planctomycetes prefer to remain attached rather than be free-living (Allgaier and Grossart, 2006). The functional roles of these associations between cyanobacteria and heterotrophic bacterial taxa are still unclear but likely will be revealed from studies exploring the microbial community functional response during a bloom and the links to environmental conditions (Steffen et al., 2012, 2017).
Previous studies that evaluated an entire year in this reservoir defined the cyanobacteria bloom period as October–March (Soares et al., 2009; Guedes et al., 2014; Rangel et al., 2016). Thus, in the period evaluated here, cyanobacterial cell density was always high, so we could not distinguish between situations with and without a bloom. The temporal analysis of the variation in the composition of the bacterial community was related to the abundance of the main cyanobacterial genera that dominated the phytoplankton community. In this late spring-summer bloom, two periods could be distinguished; the first period was characterized by higher retention time, low transparency and dominance of Microcystis. This period was also characterized by the high abundance of Bacteroidetes OTUs, particularly Cytophagales. Some bacteria of the genus Cytophaga are described as predatory agents in freshwater and marine systems (Daft and Stewart, 1971; Imai et al., 1993; Rashidan and Bird, 2001; Kirchman, 2002). In a previous study, Daft and Stewart (1971) first described a Cytophaga strain whose dynamics closely followed Microcystis dynamics; the species was able to terminate a bloom. Rashidan and Bird (2001) found a close relationship between the abundance of Anabaena sp. and a lytic strain of Cytophaga; they suggested that one reason for the decay of the Anabaena bloom was the lysis induced by this Cytophaga strain. Moreover, in a broad study that followed the microbial community of Lake Champlain by for 8 years using 16S rDNA sequencing, Tromas et al. (2017) defined Cytophaga as a bloom biomarker. In the present study, Cytophagales OTUs co-occurred with Microcystis in the first period of the bloom. Thus, we can speculate that the lytic capacity of these bacteria affects other cyanobacteria but not Microcystis, that the lytic capacity differentially affects various strains of Microcystis, or that such predatory activity is present but is not sufficient to decrease this population of cyanobacteria. Altogether, these observations indicate that Cytophaga can be an important biotic factor contributing to the prevalence of the Microcystis bloom in the Funil Reservoir.
Among the numerous negative correlations observed for Microcystis were those with the genera Planctomyces, Limnohabitans, and Hyphomicrobium. These genera have already been associated with Microcystis by others, but establishing positive correlations. Cai et al. (2013) observed an increased abundance of Planctomycetes associated with Microcystis colonies and pointed to a possible role of Planctomycetes in the degradation of sulfated polysaccharides produced by cyanobacteria. Comparing the particle-associated and free-living bacteria community from a lake during a M. aeruginosa bloom, Yang et al. reported an increase in the relative abundance of Limnohabitans and Limnobacter, genera known for utilizing algal-derived dissolved organic matter (Yang et al., 2017). Finally, Hyphomicrobium has also been identified in bacterial communities associated to M. aeruginosa, involved in cyanopeptide degradation (Briand et al., 2016; Tsao et al., 2017). In the present study, rare positive correlations were found between Microcystis and heterotrophic bacteria, among them with the order Sphingomonadales which includes strains from Sphingopyxis and Sphingomonas that can degrade microcystin (Kormas and Lymperopoulou, 2013).
The second period of the summer bloom was characterized by the dominance of Synechococcus (16S rDNA) and C. raciborskii (MO), which coincided with lower retention time and less turbidity. In this period, Planctomycetes increased in relative abundance, and five Planctomycetes OTUs were positively correlated with Synechococcus or C. raciborskii. Specifically, we observed a strong correlation between Planctomyces and Synechococcus. Ruber et al. (2017) also observed high Planctomycetes relative abundance (up to 25%) in a lake dominated by Synechococcus. Planctomyces was also identified as the most abundant genus within the phylum Planctomycetes in a bacterioplankton community of a deep alpine lake in which Synechococcus accounted for circa 30% of cyanobacteria (Salmaso et al., 2017). We also identified positive correlations between Synechococcus and other genera such as Prosthecobacter, Luteolibacter, Mycobacterium, Pseudanabaena, and Fluviicola. The co-occurrence of Synechococcus with all these taxa has been reported in natural communities (Bertos-Fortis et al., 2016; Salmaso et al., 2017) and members of Flavobacteriales have been identified as dominant associated heterotrophic bacteria in cultures of Synechococcus isolates, (Zheng et al., 2017).
In contrast to what has been reported for Microcystis and Synechococcus, little is known about bacterial groups directly associated with C. raciborskii. Here, we observed positive correlations between C. raciborskii and the genera Gemmatimonas and Roseococcus. Gemmatimonas related OTUs have been identified in bacterial communities associated with C. raciborskii cultures (Bagatini et al., 2014).
It is important to highlight that any disagreement observed in the correlation dynamics from the present study to others may represent real differences in biological dynamics, but also can be derived from the use of different methodologies. The present study has the advantage of presenting a very comprehensive survey of the microbial community in the reservoir, but as many studies using the next generation sequencing, it lacks reliable taxonomical resolution at the species level, which limits the description of ecologically relevant correlations with reliability.
While other studies have analyzed the bacterial community associated with blooms of cyanobacteria by considering factors associated with the rise and decay of blooms, some investigations, including the present study, focused on variations that occur during a bloom period; this allowed us to detail not only the influence of cyanobacteria biomass but also the possible favoring of different bacterial communities by different cyanobacterial species (or vice versa). This approach is a first step to understanding the mechanisms associated with the shift of the dominance of cyanobacteria genera, which until now have been based almost solely on interactions that involve cyanobacteria alone such as direct competition or allelopathy. Future studies under laboratory conditions can simulate the interactions among the cyanobacterial and bacterial genera described here, evaluating possible synergistic or antagonist relationships, with the potential to develop biocontrol tools for cyanobacterial blooms.
Conclusion
Our study has provided insights on the bacterial communities associated with bloom-forming cyanobacteria in a tropical system. Temporal shifts in the dominance of cyanobacteria genera were associated not only with physical features of the water (retention time and transparency) but also with shifts in the associated heterotrophic bacteria. Our results suggest specific interactions of the main cyanobacterial genera with certain groups of the heterotrophic bacterial community, but further studies exploring the microbial community functional aspects and environmental conditions are needed to better understand the ecophysiological role of these associations.
Author Contributions
IG: Conceived the study design, conducted the sampling, performed sequencing and data analysis and wrote the paper; CR: Conducted the bioinformatic analysis; LR: Conducted sampling, phytoplankton quantification, and analyzed data; LS: Analyzed data and edited the manuscript; PB: Provided counseling in sequencing and edited the manuscript; SA: Conceived the study design, analyzed data and edited the manuscript; AP: Conceived the study design, performed sequencing and data analysis, and wrote the paper. All authors approved the final submitted manuscript.
Funding
This work was supported by CNPq (446216/2015-1) and FAPERJ (E-26/210.146/2014). This study was conducted under the auspices of the CAPES (Brazil) NUFFIC (The Netherlands) project 045 2012. IG was supported by a fellowship from the National Council for Scientific and Technological Development (141668/2013-0).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We are thankful to Marcelo M. Manzi for the nutrient analysis. We also thank Marcelo M. Manzi, Vera Huszar and the Center for Environmental Studies and Sustainable Development of Universidade Estadual do Rio de Janeiro in Ilha Grande for providing an intense and pleasant writing week for preparing this manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.00424/full#supplementary-material
Supplementary Figure 1. Map of the Funil Reservoir showing geographical location and sampling stations (Modified from Soares et al., 2009).
Supplementary Figure 2. Rarefaction curves of 16SrDNA sequences from all samples.
Supplementary Table 1. Cyanobacterial species identified by microscopic analysis.
Supplementary File 1. Dominant cyanobacterial OTUs as represented in Figure 6. OTU identification, number of retrieved sequences for each OTU, genus classification and representative sequence centroid.
References
Allgaier, M., and Grossart, H. P. (2006). Seasonal dynamics and phylogenetic diversity of free-living and particle-associated bacterial communities in four lakes in northeastern Germany. Aquat. Microb. Ecol. 45, 115–128. doi: 10.3354/ame045115
Antunes, J. T., Leão, P. N., and Vasconcelos, V. M. (2015). Cylindrospermopsis raciborskii: review of the distribution, phylogeography, and ecophysiology of a global invasive species. Front. Microbiol. 6:473. doi: 10.3389/fmicb.2015.00473
Assenov, Y., Ramírez, F., Schelhorn, S. E., Lengauer, T., and Albrecht, M. (2008). Computing topological parameters of biological networks. Bioinformatics 24, 282–284. doi: 10.1093/bioinformatics/btm554
Bagatini, I. L., Eiler, A., Bertilsson, S., Klaveness, D., Tessarolli, L. P., and Vieira, A. A. (2014). Host-specificity and dynamics in bacterial communities associated with bloom-forming freshwater phytoplankton. PLoS ONE 9:e85950. doi: 10.1371/journal.pone.0085950
Berger, C., Ba, N., Gugger, M., Bouvy, M., Rusconi, F., Couté, A., et al. (2006). Seasonal dynamics and toxicity of Cylindrospermopsis raciborskii in Lake Guiers (Senegal, West Africa). FEMS Microbiol. Ecol. 57, 355–66. doi: 10.1111/j.1574-6941.2006.00141.x
Bertos-Fortis, M., Farnelid, H. M., Lindh, M. V., Casini, M., Andersson, A., Pinhassi, J., et al. (2016). Unscrambling cyanobacteria community dynamics related to environmental factors. Front. Microbiol. 7:625. doi: 10.3389/fmicb.2016.00625
Bonilla, S., Aubriot, L., Soares, M. C., González-Piana, M., Fabre, A., Huszar, V. L., et al. (2012). What drives the distribution of the bloom-forming cyanobacteria Planktothrix agardhii and Cylindrospermopsis raciborskii? FEMS Microbiol. Ecol. 79, 594–607. doi: 10.1111/j.1574-6941.2011.01242.x
Briand, E., Humbert, J. F., Tambosco, K., Bormans, M., and Gerwick, W. H. (2016). Role of bacteria in the production and degradation of Microcystis cyanopeptides. Microbiologyopen 5, 469–478. doi: 10.1002/mbo3.343
Cai, H., Yan, Z., and Wang, A. (2013). Analysis of the attached microbial community on mucilaginous cyanobacterial aggregates in the eutrophic lake taihu reveals the importance of Planctomycetes 66, 73–83. doi: 10.1007/s00248-013-0224-1
Callieri, C. (2007). Picophytoplankton in freshwater ecosystems : the importance of small-sized phototrophs. Freshw. Rev. 1, 1–28. doi: 10.1608/FRJ-1.1.1
Cheng, C., Zaichao, Z., Aizhong, D., Jiayan, W., Jingfa, X., and Yujiao, S. (2011). Bar-Coded pyrosequencing reveals the bacterial community during microcystis water bloom in guanting reservoir, Beijing. Proc. Eng. 18, 341–346. doi: 10.1016/j.proeng.2011.11.054
Clarke, K. (1993). Non-parametric multivariate analyses of changes in community structure. Aust. J. Ecol. 18, 117–143. doi: 10.1111/j.1442-9993.1993.tb00438.x
Daft, M., and Stewart, W. D. (1971). Bacterial pathogens of freshwater blue-green algae. New Phytol. 7, 819–829. doi: 10.1111/j.1469-8137.1971.tb02582.x
Dziallas, C., and Grossart, H. (2011). Temperature and biotic factors influence bacterial communities associated with the cyanobacterium Microcystis sp. Environ. Microbiol. 13, 1632–1641. doi: 10.1111/j.1462-2920.2011.02479.x
Edgar, R. C., Haas, B., Clemente, J., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381
Eiler, A., and Bertilsson, S. (2004). Composition of freshwater bacterial communities associated with cyanobacterial blooms in four Swedish lakes. Environ. Microbiol. 6, 1228–43. doi: 10.1111/j.1462-2920.2004.00657.x
Ger, K. A., Urrutia-Cordero, P., Frost, P. C., Hansson, L. A., Sarnelle, O., Wilson, A. E., et al. (2016). The interaction between cyanobacteria and zooplankton in a more eutrophic world. Harmful Algae 54, 128–144. doi: 10.1016/j.hal.2015.12.005
Gerphagnon, M., Macarthur, D. J., Latour, D., Gachon, C. M., Van Ogtrop, F., Gleason, F. H., et al. (2015). Microbial players involved in the decline of filamentous and colonial cyanobacterial blooms with a focus on fungal parasitism. Environ. Microbiol. 17, 2573–2587. doi: 10.1111/1462-2920.12860
Graham, L. E., Graham, J. M., Wilcox, L. W., and Cook, M. E. (2016). Algae, 3rd Edn. LJLM Press. Available online at: http://www.ljlmpress.com/algae.html
Guedes, I. A., da Costa Leite, D. M., Manhães, L. A., Bisch, P. M., Azevedo, S. M. F. O., Pacheco, A. B. F., et al. (2014). Fluctuations in microcystin concentrations, potentially toxic Microcystis and genotype diversity in a cyanobacterial community from a tropical reservoir. Harmful Algae 39, 303–309. doi: 10.1016/j.hal.2014.09.001
Hammer, Ø., Harper, D. A. T., and Ryan, P. D. (2001). PAST: Paleontological statistics software package for education and data analysis. Palaeontol. Electron. 4:9. Available online at: http://palaeo-electronica.org/2001_1/past/issue1_01.htm
Harke, M. J., and Gobler, C. J. (2013). Global transcriptional responses of the toxic cyanobacterium, Microcystis aeruginosa, to nitrogen stress, phosphorus stress, and growth on organic matter. PLoS ONE 8:e69834. doi: 10.1371/journal.pone.0069834
Harke, M. J., Steffen, M. M., Gobler, C. J., Otten, T. G., Wilhelm, S. W., Wood, S. A., et al. (2016). A review of the global ecology, genomics, and biogeography of the toxic cyanobacterium, Microcystis spp. Harmful Algae 54, 4–20. doi: 10.1016/j.hal.2015.12.007
Hillebrand, H., Dürselen, C. D., Kirschtel, D., Pollingher, U., Zohary, T., Durselen, C. D., et al. (1999). Biovolume calculation for pelagic and benthic microalgae. J. Phycol. 35, 403–424. doi: 10.1046/j.1529-8817.1999.3520403.x
Imai, I., Ishida, Y., and Hata, Y. (1993). Killing of marine phytoplankton by a gliding bacterium Cytophaga sp., isolated from the coastal sea of Japan. Mar. Biol. Int. J. Life Ocean Coast. Waters 116, 527–532.
Kirchman, D. L. (2002). The ecology of Cytophaga-Flavobacteria in aquatic environments. FEMS Microbiol. Ecol. 39, 91–100. doi: 10.1016/S0168-6496(01)00206-9
Klindworth, A., Pruesse, E., Schweer, T., Peplies, J., Quast, C., Horn, M., et al. (2013). Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 41, 1–11. doi: 10.1093/nar/gks808
Kormas, K. A., Vardaka, E., Moustaka-Gouni, M., Kontoyanni, V., Petridou, E., Gkelis, S., et al. (2010) Molecular detection of potentially toxic cyanobacteria their associated bacteria in lake water column sediment. World J. Microbiol. Biotechnol. 26, 1473–1482 doi: 10.1007/s11274-010-0322-x
Kormas, K., and Lymperopoulou, D. S. (2013) Cyanobacterial toxin degrading bacteria: who are they? Biomed. Res. Int. 2013:463894. doi: 10.1155/2013/463894
Kosten, S., Huszar, V. L. M., Bécares, E., Costa, L. S., van Donk, E., Hansson, L.-A., et al. (2012). Warmer climates boost cyanobacterial dominance in shallow lakes. Glob. Change Biol. 18, 118–126. doi: 10.1111/j.1365-2486.2011.02488.x
Kruk, C., Huszar, V. L. M., Peeters, E. T. H. M., Bonilla, S., Costa, L., Lürling, M., et al. (2010). A morphological classification capturing functional variation in phytoplankton. Freshw. Biol. 55, 614–627. doi: 10.1111/j.1365-2427.2009.02298.x
Li, J., Zhang, J., Liu, L., Fan, Y., Li, L., Yang, Y., et al. (2015). Annual periodicity in planktonic bacterial and archaeal community composition of eutrophic Lake Taihu. Sci. Rep. 5:15488. doi: 10.1038/srep15488
Li, N., Zhang, L., Li, F., Wang, Y., Zhu, Y., Kang, H., et al. (2011). Metagenome of microorganisms associated with the toxic Cyanobacteria Microcystis aeruginosa analyzed using the 454 sequencing platform. Chinese J. Oceanol. Limnol. 29, 505–513. doi: 10.1007/s00343-011-0056-0
Louati, I., Pascault, N., Debroas, D., and Bernard, C. (2015). Structural diversity of bacterial communities associated with bloom-forming freshwater cyanobacteria differs according to the cyanobacterial Genus. PLoS ONE 10:e0140614. doi: 10.1371/journal.pone.0140614
Lürling, M., Eshetu, F., Faassen, E. J., Kosten, S., and Huszar, V. L. M. (2013). Comparison of cyanobacterial and green algal growth rates at different temperatures. Freshw. Biol. 58, 552–559. doi: 10.1111/j.1365-2427.2012.02866.x
Lürling, M., van Oosterhout, F., and Faassen, E. (2017). Eutrophication and warming boost cyanobacterial biomass and microcystins. Toxins 9:64. doi: 10.3390/toxins9020064
Morris, R. M., Longnecker, K., and Giovannoni, S. J. (2006). Pirellula and OM43 are among the dominant lineages identified in an Oregon coast diatom bloom. Environ. Microbiol. 8, 1361–1370. doi: 10.1111/j.1462-2920.2006.01029.x
Newton, R. J., Jones, S. E., Eiler, A., McMahon, K. D., and Bertilsson, S. (2011). A guide to the natural history of freshwater lake bacteria. Microbiol. Mol. Biol. Rev. 75, 14–49. doi: 10.1128/MMBR.00028-10
O'Brien, K., Burford, M., and Brookes, J. D. (2009). Effects of light history on primary productivity in a Cylindrospermopsis raciborskii dominated reservoir. Freshw. Biol. 54, 272–282. doi: 10.1111/j.1365-2427.2008.02106.x
O'Neil, J. M., Davis, T. W., Burford, M., and Gobler, C. J. (2012). The rise of harmful cyanobacteria blooms: the potential roles of eutrophication and climate change. Harmful Algae 14, 313–334. doi: 10.1016/j.hal.2011.10.027
Osman, O. A., Beier, S., Grabherr, M., and Bertilsson, S. (2017). Interactions of freshwater cyanobacteria with bacterial antagonists. Appl. Environ. Microbiol. 83, 1–18. doi: 10.1128/AEM.02634-16
Paerl, H. (2008). Nutrient and other environmental controls of harmful cyanobacterial blooms along the freshwater-marine continuum. Adv. Exp. Med. Biol. 619, 217–237. doi: 10.1007/978-0-387-75865-7_10
Paerl, H. W., and Huisman, J. (2008). Blooms like it hot. Science 320, 57–58. doi: 10.1126/science.1155398
Paerl, H. W., and Otten, T. G. (2013). Environmental science. blooms bite the hand that feeds them. Science 342, 433–434. doi: 10.1126/science.1245276
Paerl, H. W., and Paul, V. J. (2012). Climate change: links to global expansion of harmful cyanobacteria. Water Res. 46, 1349–1363. doi: 10.1016/j.watres.2011.08.002
Parulekar, N. N., Kolekar, P., Jenkins, A., Kleiven, S., Utkilen, H., Johansen, A., et al. (2017). Characterization of bacterial community associated with phytoplankton bloom in a eutrophic lake in South Norway using 16S rRNA gene amplicon sequence analysis. PLoS ONE 12:e0173408. doi: 10.1371/journal.pone.0173408
Parveen, B., Ravet, V., Djediat, C., Mary, I., Quiblier, C., Debroas, D., et al. (2013). Bacterial communities associated with Microcystis colonies differ from free-living communities living in the same ecosystem. Environ. Microbiol. Rep. 5, 716–724. doi: 10.1111/1758-2229.12071
Rangel, L. M., Carolina, M., Soares, S., Paiva, R., Helena, L., and Silva, S. (2016). Morphology-based functional groups as effective indicators of phytoplankton dynamics in a tropical cyanobacteria-dominated transitional river – reservoir system. Ecol. Indic. 64, 217–227. doi: 10.1016/j.ecolind.2015.12.041
Rangel, L. M., Silva, L. H. S., Rosa, P., Roland, F., and Huszar, V. L. M. (2012). Phytoplankton biomass is mainly controlled by hydrology and phosphorus concentrations in tropical hydroelectric reservoirs. Hydrobiologia 693, 13–28. doi: 10.1007/s10750-012-1083-3
Rashidan, K. K., and Bird, D. F. (2001). Role of predatory bacteria in the termination of a cyanobacterial bloom. Microb. Ecol. 41, 97–105. doi: 10.1007/s002480000074
Reynolds, C. S. (2006). The Ecology of Phytoplankton - Ecology, Biodiversity and Conservation. Cambridge: Cambridge University Press.
Reynolds, C. S., Huszar, V., Kruk, C., Naselli-Flores, L., and Melo, S. (2002). Towards a functional classification of the freshwater phytoplankton. J. Plankton Res. 24, 417–428. doi: 10.1093/plankt/24.5.417
Rigosi, A., Carey, C. C., Ibelings, B. W., and Brookes, J. D. (2014). The interaction between climate warming and eutrophication to promote cyanobacteria is dependent on trophic state and varies among taxa. Limnol. Oceanogr. 59, 99–114. doi: 10.4319/lo.2014.59.1.0099
Robarts, R. D., and Zohary, T. (1987). Temperature effects on photosynthetic capacity, respiration, and growth rates of bloom-forming cyanobacteria. New Zeal. J. Mar. Fresh. 21, 391–399. doi: 10.1080/00288330.1987.9516235
Ruber, J., Bauer, F. R., Millard, A. D., Raeder, U., Geist, J., and Zwirglmaier, K. (2017). Synechococcus diversity along a trophic gradient in the Osterseen Lake District, Bavaria. Microbiology 162, 2053–2063. doi: 10.1099/mic.0.000389
Salmaso, N., Albanese, D., Capelli, C., Boscaini, A., Pindo, M., and Donati, C. (2017). Diversity and cyclical seasonal transitions in the bacterial community in a large and deep perialpine lake. Microb. Ecol. 1–19. doi: 10.1007/s00248-017-1120-x
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open source, platform independent, community supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09
Soares, M., de A Rocha, M., Marinho, M., Azevedo, S., Branco, C., and Huszar, V. (2009). Changes in species composition during annual cyanobacterial dominance in a tropical reservoir: physical factors, nutrients and grazing effects. Aquat. Microb. Ecol. 57, 137–149. doi: 10.3354/ame01336
Soares, M. C. S., Huszar, V. L. M., Miranda, M. N., Mello, M. M., Roland, F., and Lürling, M. (2013). Cyanobacterial dominance in Brazil: distribution and environmental preferences. Hydrobiologia 717, 1–12. doi: 10.1007/s10750-013-1562-1
Soares, M. C. S., Marinho, M. M., Azevedo, S. M. O. F., Branco, C. W. C., and Huszar, V. L. M. (2012). Eutrophication and retention time affecting spatial heterogeneity in a tropical reservoir. Limnologica 42, 197–203. doi: 10.1016/j.limno.2011.11.002
Steffen, M. M., Davis, T. W., McKay, R. M. L., Bullerjahn, G. S., Krausfeldt, L. E., Stough, J. M. A., et al. (2017). Ecophysiological examination of the lake erie microcystis bloom in 2014: linkages between biology and the water supply shutdown of Toledo, OH. Environ. Sci. Technol. 51, 6745–6755. doi: 10.1021/acs.est.7b00856
Steffen, M. M., Li, Z., Effler, T. C., Hauser, L. J., Boyer, G. L., and Wilhelm, S. W. (2012). Comparative metagenomics of toxic freshwater cyanobacteria bloom communities on two continents. PLoS ONE 7:e44002. doi: 10.1371/journal.pone.0044002
Stubbendieck, R. M., Vargas-Bautista, C., and Straight, P. D. (2016). Bacterial communities: interactions to scale. Front. Microbiol. 7:1234. doi: 10.3389/fmicb.2016.01234
Tromas, N., Fortin, N., Bedrani, L., Terrat, Y., Cardoso, P., Bird, D., et al. (2017). Characterizing and predicting cyanobacterial blooms in an 8-year amplicon sequencing time-course. ISME J. 11, 1746–1763. doi: 10.1038/ismej.2017.58
Tsao, S., Wei, D. J., Chang, Y. T., and Lee, J. F. (2017). Aerobic biodegradation of microcystin-LR by an indigenous bacterial mixed culture isolated in Taiwan. Int. Biodeterior. Biodegr. 124, 101–108. doi: 10.1016/j.ibiod.2017.04.011
Utermöhl, H. (1958). Zur Vervollkommnung der quantitativen Phytoplankton-Methodik: Mit 1 Tabelle und 15 abbildungen im Text und auf 1 Tafel. Int. Ver. Theor. Angew. Limnol. 9, 1–38. doi: 10.1080/05384680.1958.11904091
Van Wichelen, J., Vanormelingen, P., Codd, G. A., and Vyverman, W. (2016). The common bloom-forming cyanobacterium Microcystis is prone to a wide array of microbial antagonists. Harmful Algae 55, 97–111. doi: 10.1016/j.hal.2016.02.009
Visser, P. M., Verspagen, J. M. H., Sandrini, G., Stal, L. J., Matthijs, H. C. P., Davis, T. W., et al. (2016). How rising CO2 and global warming may stimulate harmful cyanobacterial blooms. Harmful Algae 54, 145–159. doi: 10.1016/j.hal.2015.12.006
Wilhelm, S. W., Farnsley, S. E., LeCleir, G. R., Layton, A. C., Satchwell, M. F., DeBruyn, J. M., et al. (2011). The relationships between nutrients, cyanobacterial toxins and the microbial community in Taihu (Lake Tai), China. Harmful Algae 10, 207–215. doi: 10.1016/j.hal.2010.10.001
Woebken, D., Teeling, H., Wecker, P., Dumitriu, A., Kostadinov, I., DeLong, E. F., et al. (2007). Fosmids of novel marine Planctomycetes from the Namibian and Oregon coast upwelling systems and their cross-comparison with planctomycete genomes. ISME J. 1, 419–435. doi: 10.1038/ismej.2007.63
Wood, S. A., Borges, H., Puddick, J., Biessy, L., Atalah, J., Hawes, I., et al. (2017). Contrasting cyanobacterial communities and microcystin concentrations in summers with extreme weather events: insights into potential effects of climate change. Hydrobiologia 785, 71–89. doi: 10.1007/s10750-016-2904-6
Woodhouse, J. N., Kinsela, A. S., Collins, R. N., Bowling, L. C., Honeyman, G. L., Holliday, J. K., et al. (2016). Microbial communities reflect temporal changes in cyanobacterial composition in a shallow ephemeral freshwater lake. ISME J. 10, 1337–1351. doi: 10.1038/ismej.2015.218
Wu, X., Xi, W., Ye, W., and Yang, H. (2007). Bacterial community composition of a shallow hypertrophic freshwater lake in China, revealed by16S rRNA gene sequences. FEMS Microbiol. Ecol. 61, 85–96. doi: 10.1111/j.1574-6941.2007.00326.x
Xing, P., Kong, F., Cao, H., Zhang, M., and Tan, X. (2007). Variations of bacterioplankton community composition during Microcystis spp. blooms in a shallow eutrophic lake. J. Freshw. Ecol. 22, 61–67. doi: 10.1080/02705060.2007.9664146
Yang, C., Wang, Q., Simon, P. N., Liu, J., Liu, L., Dai, X., et al. (2017). Distinct network interactions in particle-associated and free-living bacterial communities during a Microcystis aeruginosa bloom in a Plateau lake. Front. Microbiol. 8:1202. doi: 10.3389/fmicb.2017.01202
Keywords: 16S rDNA, Microcystis, Synechococcus, Cylindrospermopsis, cyanobacterial bloom, microbial community, Illumina
Citation: Guedes IA, Rachid CTCC, Rangel LM, Silva LHS, Bisch PM, Azevedo SMFO and Pacheco ABF (2018) Close Link Between Harmful Cyanobacterial Dominance and Associated Bacterioplankton in a Tropical Eutrophic Reservoir. Front. Microbiol. 9:424. doi: 10.3389/fmicb.2018.00424
Received: 19 September 2017; Accepted: 22 February 2018;
Published: 12 March 2018.
Edited by:
Petra M. Visser, University of Amsterdam, NetherlandsReviewed by:
Konstantinos Ar. Kormas, University of Thessaly, GreeceJanelle Renee Thompson, Massachusetts Institute of Technology, United States
Copyright © 2018 Guedes, Rachid, Rangel, Silva, Bisch, Azevedo and Pacheco. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ana B. F. Pacheco, biafp@biof.ufrj.br; biafpac@gmail.com