Effect of Small-Scale Turbulence on the Physiology and Morphology of Two Bloom-Forming Cyanobacteria.

The main goal of the present work is to test the hypothesis that small-scale turbulence affected physiological activities and the morphology of cyanobacteria in high turbulence environments. Using quantified turbulence in a stirring device, we conducted one set of experiments on cultures of two strains of cyanobacteria with different phenotypes; i.e., unicellular Microcystis flos-aquae and colonial Anabaena flos-aquae. The effect of small-scale turbulence examined varied from 0 to 8.01×10-2 m2s-3, covering the range of turbulence intensities experienced by cyanobacteria in the field. The results of photosynthesis activity and the cellular chlorophyll a in both strains did not change significantly among the turbulence levels, indicating that the potential indirect effects of a light regime under the gradient of turbulent mixing could be ignored. However, the experiments demonstrated that small-scale turbulence significantly modulated algal nutrient uptake and growth in comparison to the stagnant control. Cellular N and C of the two stains showed approximately the same responses, resulting in a similar pattern of C/N ratios. Moreover, the change in the phosphate uptake rate was similar to that of growth in two strains, which implied that growth characteristic responses to turbulence may be dependent on the P strategy, which was correlated with accumulation of polyphosphate. Additionally, our results also showed the filament length of A. flos-aquae decreased in response to high turbulence, which could favor enhancement of the nutrient uptake. These findings suggested that both M. flos-aquae and A. flos-aquae adjust their growth rates in response to turbulence levels in the ways of asynchronous cellular stoichiometry of C, N, and P, especially the phosphorus strategy, to improve the nutrient application efficiency. The fact that adaptation strategies of cyanobacteria diversely to turbulence depending on their physiological conditions presents a good example to understand the direct cause-effect relationship between hydrodynamic forces and algae.

Introduction

Accelerated eutrophication of freshpapan class="Chemical">n class="Chemical">water bodies oftenpan> leads to the excessive proliferation of phytoplanktonpan>, of which cyanobacteria tenpan>d to predominate. Rapid growth of cyanobacteria oftenpan> leads to blooms that can have serious adverse effects onpan> aquatic communpan>ities in lakes throughout the world [1]. Cyanobacteria bloom occurrenpan>ce is affected by a complex set of physical, chemical, biological, hydrological, and meteorological conpan>ditionpan>s, so it is difficult to isolate specific causative enpan>vironpan>menpan>tal factors [2]. Previous investigationpan>s suggested that interactionpan>s amonpan>g nutrienpan>t inputs, turbulenpan>t mixing, unpan>derwater light availability, and grazing pressure can strongly influence the components of phytoplankton species in the water column [3, 4].

Turbulence mixipan class="Chemical">ng has been widely accepted as a source of external energy to the planktonic environment, such as regulating light accessibility in the water column [5], bringing eutrophic water into contact with algal cells [6, 7], and affecting encounter possibilities between phytoplankton and predators [8]. In particular, small-scale turbulence could play an important role in determining which phytoplanktonic taxon dominates over other competitors in the ever-changing environment [9]. Phytoplankton in the water column are constantly exposed to small-scale turbulence, which may have a differential effect on the physiology of algae [10]. The positive effect of small-scale turbulence on phytoplankton has been reported to occur through changing diffusive sublayers and regulating nutrient fluxes of cells [11-13], while the negative impacts include mechanical damage, behavioral alteration and physiological impairment [10, 14–16]. However, to date, there is only a small body of literature on the direct effects of small-scale turbulence on cyanobacteria in the freshwater rivers and lakes.

Cyanobacteria have copan class="Chemical">nstantly been subjected to a complex field environment imposed by these small-scale turbulent conditions. Survival of different cyanobacteria species in the mixed water column depenpan>ds onpan> the combined effects of varying enpan>vironpan>menpan>tal factors [17]. Nutrienpan>ts and light are crucial, but the turbulenpan>ce and its inpan>tensity also affect growth strategies of cyanobacteria species in their competitionpan> for these basic commodities [18, 19]. Because of a variety of adaptive strategies, such as buoyancy [20], N2-fixationpan> [21], light tolerance [22], and efficienpan>t utilizationpan> of phosphorus and inorganic carbon [23, 24], cyanobacteria can dominate over other competitors throughout the life cycle. Under certain environmental stressors, cyanobacteria adapt habitat energy and nutrients through the adjustment of their metabolism while keeping cellular carbon, nitrogen and phosphorus in an acceptable range [25, 26]. Nevertheless, this plasticity of adaptation regarding nutrients has not been well reported when cyanobacteria species encounter turbulence mixing.

Many studies have reported that heterogepan class="Chemical">neity in response to turbulence is often observed among phytoplankton taxa [27], and certain species may exhibit specific adaptive strategies to turbulence. For instance, genus-specific shear effects have been illustrated for cyanobacteria such that the non-n class="Disease">heterocystous filamentous Oscillatoria appears to be more shear tolerant than the heterocystous Anabaenpan>a [28, 29]. Furthermore, the seasonpan>al field studies of mixed phytoplanktonpan> communpan>ities and enpan>closure experimenpan>ts also have shownpan> that shear has either increased or decreased/arrested cyanobacterial abunpan>dance or activity [28-30]. Some researchers suggested these differenpan>ces in senpan>sitivity induced by turbulenpan>ce depenpan>d onpan> cell size. It has beenpan> reported that the increase in the nutrienpan>t diffusionpan> rate into cells due to small-scale turbulenpan>ce varies from negligible to signpan>ificant as cell size increases [13].

Hard evidence of a direct cause—effect relatiopan class="Chemical">nships between turbulence and algae is not easy to find [8, 31]. Inpan> the presenpan>t study, we inpan>vestigated the effect of small-scale turbulence onpan> the physiology and morphology of two bloom-forming cyanobacteria with differenpan>t phenpan>otypes, i.e., unpan>icellular Microcystis flos-aquae and colonpan>ial Anabaenpan>a flos-aquae. To our knowledge, few experimenpan>ts have beenpan> performed employing quantified turbulenpan>ce levels to evaluate the effect of small-scale turbulenpan>ce onpan> the physiology of cyanobacteria. Thus we developed a model to quantify turbulenpan>ce exposure in the experimenpan>t and systematically investigate the change in photosynthesis and total cellular C, N and P of the two stains, and filamenpan>t lenpan>gth of A. flos-aquae under different turbulent dissipation rates. The data from the study were analyzed to illustrate the mechanisms of cyanobacteria with different morphology adaptation strategies under small-scale turbulence mixing.

Materials and Methods

Strains and culture conditions

Unicellular papan class="Chemical">n class="Species">M. flos-aquae FACHB1028 anpan>d colonial A. flos-aquae FACHB245 were obtained from the Freshwater Algae Culture Collection of the Institute of Hydrobiology (FACHB-Collection, Wuhan, China). Routine cultures of the two strains were maintained in BG11 medium with a phosphate concentration of 2 mg/L [32] under a cool-white fluorescent light intensity of 25 μmol m-2 s-1 with a 12:12 light:dark cycle at a temperature of 25±1°C, and they were manually shaken 3 times every day during incubation.

In the experimepan class="Chemical">nts, cells at their exponential growth phase were served with P free BG11 medium for a period of 10 days before the growing experimenpan>t. After that, they were inoculated into 1 L glass beakers with 2 mg/L npan> class="Chemical">phosphate concentration BG11 medium and exposed to different agitation intensities. A 1-L beaker was used as the reactor in the experiment. A self-designed impeller was located at the center of the beaker. Geometric dimensions are shown in Fig 1. The liquid height was 113 mm. Continuous stirring was performed at the start of the growing experiment at rotation speeds of 100, 200, 300, 400 and 500 rpm. The samples cultured in completely stagnant conditions were used as controls. The test was run in triplicate for two weeks under the conditions described above, and they were sampled every 2 days for the determination of varying parameters, as described below.

Fig 1

Geometric dimensions of the reactor (Unit: mm).

Modeling and calculation of turbulent dissipation rates

A two-phase flow Computational Fluid Dypan class="Chemical">namics (CFD) model was used to simulate hydrodynamics in the reactor under different rotation speed [33]. The whole reactor was considered as the computational domain and was discretized using unstructured grids, where finer grids were performed in impeller and shaft regions (maximum size = 1 mm). About 1,000,000 total computational grids (maximum size = 2 mm) were created using the tetra mesh optionpan> of ANSYS ICEM CFD (ANSYS Inpan>c.) in order to get the grid indepenpan>denpan>t solutionpan> for the flow. The simulationpan> of the process was conpan>sidered as a steady-state, and numerically solved by the commercial software package CFX 12.0 (ANSYS Inpan>c.). The conpan>vergenpan>ce criteria used in all the simulationpan>s was 1×10−4, which was a factor by which the initial mass flow residual reduced as the simulationpan> progressed. The simulationpan>s were carried out on the 6 nodes, 96 processor SUGON cluster. A commercial 2D particle image velocimetry (PIV) system (ILA, German), including a laser (Leamtech, Nd:YAG, 200Mj, 15Hz), a PIV camera (PCO2000, 2048×2048 pixels), a synchronizer and Vidpiv 4.7 software, was employed to validate the model. The setup and data processing were the same as described earlier [34]. Detail information of its validation results were reported in S1 File. The integrated results of turbulence mixing rate and Kolmogorov microscale in different levels of rotation speed were shown in Table 1.

Table 1

Turbulent dissipation rate and Kolmogorov microscale in different levels of rotation speed estimated by the CFD hydrodynamic model.

Turbulence
Rotation speed (rpm)	100	200	300	400	500
Turbulent dissipation rate (m2s-3)	1.51×10−3	6.63×10−3	2.26×10−2	5.06×10−2	8.01×10−2
Kolmogorov microscale (mm)	0.186	0.126	0.092	0.074	0.065

Determination of growth rates and chlorophyll a (Chla) of cyanobacteria

Two strains were batch cultured upan class="Chemical">nder different agitation intenpan>sities (0 rpm to 500 rpm) for two weeks anpan>d were sampled every 2 days for cell counting unpan>der a microscope (Olympus CX 43, Japan) using the appropriate magnpan>ificationpan> and a hemocytometer. Prior to cell counpan>ting, A. flos-aquae colonies were treated with ultrasonication for 30 s or 60 s at 20 kHz and 100 W (JY92-2D, SCIENTZ, China) to obtain single cells. The growth rate was measured during the exponential growth phase according to the exponential growth equation: μ = (lnBt2-lnBt1)/Δt [35], where μ is the specific growth rate, Bt2 and Bt1 represent the cell concentration at the end (t2) and beginning (t1) of the log phase, respectively, and Δt = t2-t1. For Chla measurements, the samples were centrifuged at 6,000 rpm for 15 min, and cell pellets were extracted with 80% acetone/20%water.

Determination of photosynthetic activity of cyanobacteria

In vivo papan class="Chemical">n class="Chemical">chlorophyll fluorescenpan>ce was measured with a phytoplankton analyzer (PHYTO-PAM, Walz GmbH, Germany). The samples were dark-adapted for at least 15 min before measuring the fluorescence parameters (photosystem II activity, PSII). The electron transport rate of PSII (ETR) was calculated according to [36]. Relative ETR = (((F- F) /F)×0.84×0.5×PAR (m-2s-1)). The maximum effective quantum yield of PSII was calculated as F/F = (F-F)/F [37]. F is the difference between F and F, which are the maximum and minimum fluorescence values of the dark adapted stage of PSII.

Analysis of the carbon, nitrogen and phosphorus contents of cyanobacteria

To measure the carbonpan>, nitrogen and phosphorus contents of the two strains, culture suspensions were taken from each sample, filtered through GF/F filters (Whatman, UK) and dried at 65°C for 24 h. Each filter was split for analysis of C, N, and total cellular phosphorus (TCP). C and N concentrations were measured using a Vario El cube elemental analyzer (Elementar, Germany). The filtrates were used for measuring soluble reactive phosphorus (SRP). The filtered extracts were diluted with distilled water, and TCP was determined by digestion [38, 39] followed by colorimetric analysis [40]. The polyphosphate (polyP) was extracted using the hot water method according to [41].

Uptake rate of phosphate

pan class="Chemical">Phosphate uptake rates were calculated by the chanpan>ge in phosphate concentration in the culture media during 72 h batch cultures at the same culture conditions as mentioned above. Phosphate levels were determined using ascorbic acid-molybdenum blue spectrophotometric methods [42]. The uptake rate was calculated according to [43]: phosphate uptake rate = (W-W) · V · N· t, where W is the initial concentration of phosphate, W is the concentration of the phosphate after t hours, N is the initial level of cell density and V is the volume of culture medium.

Measurement of the filament size of Anabaena

Anabaepan class="Chemical">na filaments were examined using an Olympus BX 53 light microscope with a digital camera (Olympus DP 71, Japan). The sizes of the filaments were measured in microphotographs using image analysis Olympus DP-Soft. One hundred filaments were chosen at random from each sample to capture the images, and then the surface area (S), width (d), and length (l) of individual filaments were obtained using Olympus DP-Soft, which was calibrated separately for each magnification.

Statistical analysis

The data in this study are presepan class="Chemical">nted as the mean ± standard deviation (SD). Significant differences between controls and treated samples were determined with an ANOVA. All statistical analyses were carried out with Origin 8.0 (OriginLab, USA). Differences were considered to be significant at P < 0.05.

Results

Effect of turbulence on the growth rate of cyanobacteria

Specific growth rates of unicellular papan class="Chemical">n class="Species">M. flos-aquae anpan>d colonial A. flos-aquae at different turbulent dissipation rates are shown in Table 2. The growth curve was plotted to roughly ascertain the period of the exponential phase of each strain, and thus the growth rates of individual strains were calculated in the exponential phase accordingly. For M. flos-aquae, the growth rate increased with the elevation of the turbulent dissipation rate. At a turbulent dissipation rate of 1.51×10−3 m2s-3, the growth rate of M. flos-aquae was approximately 3.4 times than that in the stagnation condition. When the turbulent dissipation rate was above 5.06×10−2 m2s-3, the growth rate of M. flos-aquae remained relatively unchanged (ANOVA, P > 0.05). On the other hand, the growth rate of A. flos-aquae increased dramatically (ANOVA, P < 0.05) at lower turbulent dissipation rates and reached the max at 2.26×10−2 m2s-3, while the growth rate decreased with increases in the turbulence intensity at higher levels. Overall, the growth rates of A. flos-aquae at all levels of turbulence intensity were lower than that of M. flos-aquae.

Table 2

The specific growth rates of M. flos-aquae and A. flos-aquae under different levels of turbulence.

Turbulent dissipation rate (m2s-3)
Strain	0	1.51×10−3	6.63×10−3	2.26×10−2	5.06×10−2	8.01×10−2
M. flos-aquae	0.039±0.003	0.133±0.002	0.144±0.001	0.223±0.005	0.261±0.005	0.286±0.004
A. flos-aquae	0.009±0.001	0.014±0.002	0.027±0.001	0.043±0.002	0.040±0.003	0.029±0.003

The data are the mean ± SD (n = 3)

The data are the mean ± SD (pan class="Chemical">n = 3)

Photosynthetic parameters and cellular Chla concentrations of cyanobacteria

After incubatiopan class="Chemical">n with varied levels of turbulent dissipation rates, the maximum electron transport rate (ETRmax) and the maximum quantum yield (F/F) of M. flos-aquae and A. flos-aquae showed considerable differences (Fig 2a and 2b). The parameters for the photosynthetic responses to changes in turbulence intensity of M. flos-aquae were significantly higher than in the control (non-turbulent) treatment (ANOVA, P < 0.05). Among all the turbulent mixing treatments, markedly low values of F/F of M. flos-aquae were found in cells grown at a turbulence intensity of 1.51×10−3 m2s-3. The values of ETRmax and F/F were not found to be significantly different among levels of turbulence intensities from 2.26×10−2 to 8.01×10−2 m2s-3 (ANOVA, P > 0.05). Coincident with the elevation in the ETRmax of M. flos-aquae, the changes in the ETRmax of A. flos-aquae showed similar trends to that of M. flos-aquae (ANOVA, P < 0.05). However, it appears that there was no significant difference in the F/F of A. flos-aquae between all the turbulent treatments and the controls (ANOVA, P > 0.05). The cellular Chla showed a similar pattern to the F/F ratios, although the cellular Chla was relatively higher at higher turbulent dissipation rates, probably due to the reduction in cell density at the highest turbulent dissipation rates (Fig 2c).

Fig 2

Changes in the maximum electron transfer rate (ETRmax) (a), F/F (b) and cellular chlorophyll a (Chla) (c) of M. flos-aquae and A. flos-aquae at different turbulent dissipation rate levels in the logarithmic growth phase. The data are the mean ± SD (n = 3).

Phosphate uptake of cyanobacteria

The time course of SRpan class="Chemical">P inpan> the culture medium and phosphate uptake rate of M. flos-aquae and A. flos-aquae cultured at different turbulent dissipation rate levels are shown in Fig 3, respectively. It was observed that the SRP concentrations of the two strains decreased significantly with culturing time at all turbulence intensity levels. The residual P concentrations in the M. flos-aquae cultures were affected by the turbulent dissipation rate. The concentration of SRP in the cultures treated with higher turbulent dissipation rates (8.01×10−2 m2s-3) were significantly lower than that in the turbulent dissipation rates of 1.51×10−3 and 6.63×10−3 m2s-3 at the end of the test (ANOVA, P < 0.05). However, the SRP concentration in the A. flos-aquae culture with a turbulence intensity of 2.26×10−2 m2s-3 was the greatest at the end of the test, whereas no significant differences were observed among the other treatment cultures (ANOVA, P > 0.05). The phosphate uptake rate was analyzed as a function of the external phosphate concentration. It showed similar trends to that of SRP (Fig 4). Under the turbulent dissipation rate of 2.26×10−2 m2s-3, the uptake rates of phosphate in the two strains increased with the elevation of turbulence intensity. When the turbulent dissipation rate was above 2.26×10−2 m2s-3, the phosphate uptake rate of M. flos-aquae remained relatively unchanged, while the phosphate uptake rate of A. flos-aquae decreased dramatically (ANOVA, P < 0.05).

Fig 3

The time course of soluble reactive phosphorus (SRP) in the culture medium of M. flos-aquae (a) and A. flos-aquae(b) at different turbulent dissipation rate levels.

The data are the mean ± SD (n = 3).

Fig 4

Uptake rate for phosphate and polyphosphate content of M. flos-aquae (a) and A. flos-aquae (b) cultured at different turbulent dissipation rate levels in the logarithmic growth phase.

The data are the mean ± SD (n = 3).

The data are the mean ± SD (pan class="Chemical">n = 3).

The data are the mean ± SD (pan class="Chemical">n = 3).

Effect of turbulence on the elemental composition of cyanobacteria

To further explore the cellular stoichiometry changes, cellular C, pan class="Chemical">N and P in response to different turbulent dissipation rates were analyzed (Table 3). There were significant differences in cellular C, N, and P among different levels of turbulent mixing (ANOVA, P < 0.05). For M. flos-aquae, cellular C was the highest unpan>der the turbulenpan>t dissipationpan> rate of 5.06×10−2 m2s-3, with an increase of 11% compared with the conpan>trol treatmenpan>t. However, the cellular N conpan>cenpan>trationpan> was signpan>ificantly decreased with increasing turbulenpan>ce intenpan>sity, except unpan>der the turbulenpan>t dissipationpan> rate of 1.51×10−3 m2s-3. Furthermore, as dissipationpan> rates increased, there was also anpan> increase in cellular P, and the maximum value was observed unpan>der the turbulenpan>t dissipationpan> rate of 2.26×10−2 m2s-3, which was 2.0-fold compared with that at the stagnpan>ationpan> conpan>ditionpan>. Additionpan>ally, the changes in cellular PolyP showed similar trends to that of cellular P (Fig 4). In A. flos-aquae, the highest concentrations of cellular C, N and P were found at the turbulent dissipation of 2.26×10−2 m2s-3. The cellular PolyP kept increasing along with increases in turbulence intensity at the rate of 0–2.26×10−2 m2s-3. Then, a decline of cellular PolyP was observed under the dissipation rate of 5.06×10−2 m2s-3 and above (ANOVA, P < 0.05).

Table 3

Cellular stoichiometry of M. flos-aquae and A. flos-aquae under different levels of turbulence in the logarithmic growth phase.

Turbulent dissipation rate (m2s-3)	M. flos-aquae					A. flos-aquae
	pg/cell			mole ratio		pg/cell			mole ratio
	C	N	P	C/N	N/P	C	N	P	C/N	N/P
0	6.14	1.58	0.29	4.53	12.06	6.13	1.38	0.25	5.18	12.22
1.51×10−3	6.70	1.89	0.31	4.13	13.50	7.27	1.48	0.19	5.73	17.24
6.63×10−3	6.73	1.38	0.26	5.68	11.75	6.38	1.31	0.28	5.68	10.35
2.26×10−2	6.36	1.49	0.59	4.97	5.59	8.17	1.55	0.33	6.14	10.40
5.06×10−2	6.89	1.16	0.43	6.92	5.97	6.81	1.38	0.21	5.75	14.55
8.01×10−2	5.52	0.91	0.33	7.07	6.11	7.22	1.49	0.2	5.65	16.49

Changes ipan class="Chemical">n molar ratios of C/N and N/P of M. flos-aquae and A. flos-aquae are also shown in Table 3. The results from our experiment showed plasticity of cellular ratios in M. flos-aquae and A. flos-aquae in response to small-scale turbulence. There was an increase in the C/N ratio and a decrease in the N/P ratio in M. flos-aquae at the dissipation rate of 6.63×10−3 m2s-3 and above, compared with the controls. Variations in the C/N ratio in A. flos-aquae were similar to that in M. flos-aquae. The maximum C/N ratio was recorded at the turbulent dissipation rate of 2.26×10−2 m2/s3. However, the N/P ratio in A. flos-aquae decreased only at the turbulent dissipation rates of 6.63×10−3 and 2.26×10−2 m2s-3, compared with the control treatment, which was due to the increase in cellular P.

The filament length of A. flos-aquae

The filament lepan class="Chemical">ngth of A. flos-aquae varied with differenpan>t turbulenpan>t dissipation rates (Fig 5a). It was observed that the geometric mean of the filamenpan>t lenpan>gths of A. flos-aquae decreased with the elevation of turbulence intensity. Compared with the stagnation condition, the filament length of A. flos-aquae decreased by 69% and 47% at the turbulent dissipation rates of 2.26×10−2 and 5.06×10−2 m2s-3, respectively. The frequency distribution of each filament length range under the different turbulence intensities is shown in Fig 5b. The result also demonstrated that the frequency of 0–50 μm lengths was increased with increasing dissipation rates and reached the maximum at the 2.26×10−2 m2s-3, while the large filament (>200 μm) was the lowest. The frequencies of 100–200 μm lengths decreased at the turbulent dissipation rates of 1.51×10−3 and 6.63×10−3 m2s-3, while for the high turbulence intensities of 8.01×10−2 m2s-3, the frequencies of large filaments (>200 μm) also decreased significantly.

Fig 5

Box and whisker plots (a) and frequency distribution (b) of A. flos-aquae filament length as a response to various turbulence intensities at the end of the test.

a, filament lengths are presented in box plots showing the 5, 10, 25, 50, 75, 90 and 95 percentiles; whiskers represent minimum and maximum values, and dots inside the boxes are the mean values.

Discussion

Turbulence is ubiquitous ipan class="Chemical">n aquatic systems and thus can potentially influence a wide scope of phytoplankton and eco-physiological processes. It is still often referred to as one of the unsolved problems in field studies due to the interactions with other variables, such as light penetration and supply of nutrients, and their consequences for organism physiology and behaviors [44]. Some of the current knowledge has been derived from laboratories with devices to generate controlled turbulence; however, these should correctly evaluate the response of phytoplankton to a certain level of turbulence in natural conditions [45]. A gradient of increasing levels of dissipation rates was used in our study to investigate the shape and character of potential relationships between small-scale turbulence, nutrient conditions and cyanobacteria. Turbulent dissipation in the device varied from 10−4 to 10−2 m2s-3. On the Kolmogorov microscale, it was in the higher range of turbulences in natural aquatic ecosystems. However, it covered the range of turbulence intensities experienced by cyanobacteria in the field [2].

In most studies described ipan class="Chemical">n the literature to date, the effect of turbulence on phytoplankton was determined by changes in growth and the viability of the organism [46, 47]. In this study, responses of Microcystis and Anabaena to turbulence intensities were quite different as inferred by growth. The turbulent regimes of 1.51×10−3 to 2.26×10−2 m2s-3 favored the growth of the two strains. However, in regard to the higher dissipation rate of 8.01×10−2 m2s-3, turbulenpan>ce had a negative impact onpan> the growth rate of A. flos-aquae but did not affect the growth of M. flos-aquae, suggesting that unicellular M. flos-aquae was better adapted to growing under high turbulent dissipation rates.

The photosynthetic system is a good ipan class="Chemical">ndicator system because it is very sensitive to any changes in living conditions [48]. F/F and ETRmax are conpan>sidered as the senpan>sitive photosynthetic parameters to enpan>vironpan>menpan>tal stress [49]. At the used range of turbulenpan>t mixing, the chanpan>ges in F/F and ETRmax in the exponential phase support the hypothesis that the photosynthetic activity of A. flos-aquae significantly increases when compared with stagnant conditions but might not change in contrast to the unimodal response of a growth rate to a gradient of turbulent mixing. Similar results were also discussed by [50], indicating that at saturating irradiances, photosynthesis was less sensitive than growth to changes in turbulence mixing. Additionally, the Chla content in both M. flos-aquae and A. flos-aquae increased with turbulence even under a higher level, indicating that the photosynthetic apparatus of algae showed little or no disruption in response to turbulence. Therefore, the potential indirect effects of a light regime under the gradient of turbulent mixing, e.g., the shading effect due to an increase in the population in the reactors, could be ignored.

A great majority of studies reported that turbulence may have bepan class="Chemical">neficial effects on phytoplankton due to an increase in the diffusive transport of nutrients toward the cell surface and, hence, an increase the nutrient uptake rate [14, 51]. Under turbulent mixing in our study, cellular N and C of the two stains showed approximately the same responses, resulting in a similar pattern of C/N ratios in response to turbulent mixing. Unlike cellular C and N, P uptake and assimilation in cells did not show a unimodal response. Turbulence is reported to affect the algal growth by changing their nutrient acquisition traits [52]. The fact that the change in phosphate uptake rates was similar to that of growth in both of the strains implied that differenpan>t growth characteristic responpan>ses to turbulenpan>ce may be depenpan>denpan>t onpan> the P strategy. This finding was conpan>sistenpan>t with previous research [47] in that the phosphate uptake rate of algae was approximately similar to the growth rate.

Formation apan class="Chemical">nd storage of PolyP is commonpan>ly regarded as a nutritionpan>al strategy of micronpan> class="Species">algae when encountering environmental stress [53]. According to Fig 4, it was shown that luxury uptake of phosphorus by cyanobacteria was mostly applied to accumulate PolyP, implying that accumulation of PolyP at high dissipation rates was the strategy for cyanobacteria to maintain phosphorus concentrations in unfavorable conditions. However, compared with M. flos-aquae, the concentration of PolyP in A. flos-aquae at dissipation rates of 2.26×10−2 to 5.06×10−2 m2s-3 were much higher than those with a dissipation rate of 1.51×10−3 m2s-3, indicating that the reallocation of cellular P in A. flos-aquae could be a response to the changed turbulent conditions, and the filament integrity in A. flos-aquae was more sensitive to turbulence than in unicellular M. flos-aquae. Furthermore, our results also showed that the filament length of A. flos-aquae decreased in response to turbulence. It seemed that turbulence negatively affects cyanobacterial morphology prior to physiological activity. We stress that this may be related to the nutrient acquisition strategies, which was depended on algal size. Theoretical predictions regarding the thickness of external diffusion boundary layer is equal to the equivalent spherical diameters of the organism [54]. Thus, larger organism increased the restriction on nutrient uptake from low external concentrations, probably because the greater thickness of the diffusion boundary layers around itself. As small filaments have higher rates of nutrient uptake and lower metabolic requirement, which allow them to grow at much lower resource concentrations than larger ones [55], we inferred that with the exhaustion of SRP in the culture, the decreased filaments of A. flos-aquae at higher turbulent dissipation rates could favor the resource gain. Our findings support what had been suggested by [47], such that, in the lotic diatoms where the cells were smaller, there was higher P affinity. Overall, our study also highlights the need for further investigation of the relationship between the boundary layer around filaments and turbulence mixing.

A gradient of increasing levels of dissipation rates is useful when attempting to investigate the shape and character of potential relationships between small-scale turbulence, nutrient conditions and planktonic organisms [56]. Small-scale shearing affects physiological activities and the morphology of phytoplankton in high turbulence environments. The variety of physiological and morphological properties of phytoplankton may be largely considered as adaptations to different scales and characteristics of turbulent motions in the aquatic environment [57]. In concert with earlier investigations, the results from the two types of cyanobacteria in our study illustrated that both M. flos-aquae and A. flos-aquae adjust their growth rates in response to turbulence levels in the ways of asynchronous cellular stoichiometry of C, N, and P, especially the phosphorus strategy, to improve the nutrient application efficiency. Additionally, the colonial A. flos-aquae could change the effective filament length to enhance the nutrient uptake under high turbulent conditions. Again, these evidences supported the hypothesis that hydrodynamic regimes may impact eco-physiology of cyanobacteria to form certain adaptive strategies, which should be carefully considered in elucidating mechanisms of cyanobacteria blooms. Clearly, future studies are warranted to assess the effects of turbulence on nutrient flux of cells including different morphological cyanobacteria species.

The validation results of CFD modeling and PIV measurement of flow pattern.

(DOCX)

Click here for additional data file.