Elevated pCO2 Impedes Succession of Phytoplankton Community From Diatoms to Dinoflagellates Along With Increased Abundance of Viruses and Bacteria

Eutrophic coastal regions are highly productive and greatly influenced by human activities. Primary production supporting the coastal ecosystems is supposed to be affected by progressive ocean acidification driven by increasing CO2 emissions. In order to investigate the effects of high pCO2 (HC) on eutrophic plankton community structure and ecological functions, we employed 9 mesocosms and carried out an experiment under ambient (∼410 ppmv) and future high (1000 ppmv) atmospheric pCO2 conditions, using in situ plankton community in Wuyuan Bay, East China Sea. Our results showed that HC along with natural seawater temperature rise significantly boosted biomass of diatoms with decreased abundance of dinoflagellates in the late stage of the experiment, demonstrating that HC repressed the succession from diatoms to dinoflagellates, a phenomenon observed during algal blooms in the East China Sea. HC did not significantly influence the primary production or biogenic silica contents of the phytoplankton assemblages. However, the HC treatments increased the abundance of viruses and heterotrophic bacteria, reflecting a refueling of nutrients for phytoplankton growth from virus-mediated cell lysis and bacterial degradation of organic matters. Conclusively, our results suggest that increasing CO2 concentrations can modulate plankton structure including the succession of phytoplankton community and the abundance of viruses and bacteria in eutrophic coastal waters, which may lead to altered biogeochemical cycles of carbon and nutrients.


INTRODUCTION
The ocean CO 2 sink has increased from 1.7 ± 0.4 Pg C yr −1 in 1980s to 2.5 ± 0.6 Pg C yr −1 in 2010s (Friedlingstein et al., 2020), leading to the pH drop of open ocean surface by 0.017-0.027 units per decade (IPCC, 2019), an environmental problem known as ocean acidification (OA). It is predicted that further reduction of around 0.3 pH units will occur by the end of this century under "business as usual" scenario (IPCC, 2019). Considering the huge changes, it is needed to investigate how marine life will respond directly or indirectly and the consequent changes in biogeochemical cycles. A number of studies have been conducted ranging from single organism to artificial community then to natural community levels, using laboratory, mesocosm and field investigation approaches (see the review by Gao et al., 2019 and literatures therein).
Field mesocosm experiments using natural communities are becoming effective approaches to explore the responses of marine planktonic ecosystems to climate change. Two mesocosm studies deployed in the Raunefjorden, western Norway, showed that Emiliania huxleyi, the most abundant coccolithophore species, reduced its growth rates under 710 µatm of pCO 2 during May-June 2001 (Engel et al., 2005), and even lost blooming capacity under 1000-3000 µatm conditions during May-June 2011 (Riebesell et al., 2016). In another mesocosm experiment carried out in the Raunefjorden during May 2006, the dominant group of prasinophyte community changed from Bathycoccus-like phylotypes to Micromonas-like phylotypes, as pCO 2 increased from current to predicted future levels (Meakin and Wyman, 2011). These results suggest that the OA effects on community structure can vary temporally in the same regions. In addition, highly variable effects of elevated pCO 2 on plankton also attribute to different regions. Under HC conditions, the composition of the community from the temperate Bering Sea shelf and offshore shifted away from diatoms toward nanophytoplankton (Hare et al., 2007). In the temperate Gullmar Fjord at the Swedish west coast, the total diatom biomass were unaffected by pCO 2 during bloom processes, but displayed a positive response to pCO 2 during post-bloom processes . Off the subtropical eastern coast of Gran Canaria, Spain, the total diatom biomass remained unaffected during oligotrophic conditions, but was significantly positively affected by high CO 2 after nutrient enrichment (Bach et al., 2019). Other plankton also showed varied responses to OA among different regions. The development of zooplankton community was unaffected significantly in the temperate Kongsfjorden Fjord of northwestern Norway (Aberle et al., 2013;Niehoff et al., 2013), but was prevented by increased pCO 2 off the subtropical eastern coast of Gran Canaria, Spain (Riebesell et al., 2018). In the mesocosm experiment conducted in the Kongsfjorden Fjord, bacterioplankton diversity differed significantly between high and low pCO 2 treatments (Zhang et al., 2013). However, pCO 2 did not affect bacterioplankton diversity, in our previous mesocosm experiment conducted in the subtropical eutrophic Wuyuan Bay (Lin et al., 2018). The shifts in plankton community structure will further influence the ecological functions, such as the production, consumption and export of organic matter, with different magnitudes and even directions (de Kluijver et al., 2013;Spilling et al., 2016;Taucher et al., 2020). By integrating and analyzing the results of mesocosm experiments, we noted that planktonic ecosystem level responses to elevated pCO 2 can vary spatiotemporally, which may depend on the initial community structure and other abiotic environmental factors such as nutrient levels.
To date, most mesocosm experiments have been conducted in oligotrophic or non-eutrophicated regions. In the mesocosm experiment employed in the Kongsfjorden during June-July 2010, increase in pCO 2 did not significantly influence chlorophyll a (Chl a) concentration under oligotrophic conditions, but increased Chl a concentration after nutrient fertilization . Increased pCO 2 similarly affected Chl a concentration after nutrient addition, in the mesocosm experiment employed in Raunefjord during May-June 2011 (Riebesell et al., 2016). In our previous mesocosm experiments which have been carried out in the highly eutrophic Wuyuan Bay, we found that HC treatment did not significantly affect Chl a concentration . The plankton community composition responses to increased pCO 2 also varied in different phases of the mesocosm experiment conducted in the Gullmar Fjord, resulting from plankton communities supported by regenerated nutrients were more sensitive to OA than those having access to higher availability of inorganic nutrients (Bach et al., 2017;Bach et al., 2019). These results suggested that planktonic ecosystem level responses to elevated pCO 2 can be influenced by nutrient levels.
Coastal waters are highly productive that support marine ecosystem and human activities. In these regions, the accumulative effects of natural and anthropogenic climate change are continuing and intensifying the course of eutrophication (Rabalais et al., 2009). With progressive OA, eutrophication is supposed to enhance the CO 2 -driven acidification in coastal waters, by an additional drop of 0.13 pH unit relative to open oceans at the end of this century (Cai et al., 2011). Such enhanced OA in pH fluctuating coastal waters can disrupt ecosystem service due to increased energetic demand for ecosystems to sustain. To better understand the consequences of OA in eutrophic coastal regions, further studies should be expanded to look into the ecological responses in different regions and at different seasons under fluctuating and changing environmental conditions. Here, we carried out a mesocosm experiment in the eutrophic coastal region Wuyuan Bay, using in situ community during April-May 2018. Apart from monitoring environmental and biogenic changes, we particularly focused on the phytoplankton community succession and the abundance of viruses and bacteria, and the results showed notable correlation between autotrophic assimilation and heterotrophic dissimilation.

Mesocosms Setup
The mesocosm experiment was conducted in the Facility for the Study of Ocean Acidification Impacts of Xiamen University, located in the subtropical coastal region Wuyuan Bay (24 • 31 48 N, 118 • 10 47 E), during 9th April-15th May 2018. Nine cylindrical mesocosm bags with 3 m in depth and 1.5 m in diameter, made of transparent thermoplastic polyurethane (TPU), were fixed within steel frames and sheltered by cone lids. The in situ seawater were filtered (pore size of 0.01 µm) by water purifier (MU801-4T, Midea, China) and pumped simultaneously into 9 mesocosm bags for 3000 L within 36 h. Bag 2, 4, 6, and 8 were set as control, aerated with ambient air pCO 2 (denoted as AC). To adjust seawater to projected 1000 µatm CO 2 in 2100, bag 1, 3, 5, 7, and 9 were added with approximately 11 L CO 2 saturation seawater (denoted as HC) . Subsequently, AC and HC mesocosm bags were aerated (5 L min −1 ) with ambient air of ∼410 ppmv CO 2 and premixed air-CO 2 of 1000 ppmv CO 2 , respectively, till the end of the experiment. After aerating to homogenize the carbonate system for about 2-3 h, 720 L of in situ seawater filtered by 180 µm mash was added simultaneously into 9 mesocosm bags, so that each bag was inoculated with 80 L in situ seawater containing a natural microbe community. Subsequently, the mesocosm samples were collected from 0.5 m depths of each bag at 10:00 a.m. every 1-3 days, for physical, chemical and biological analysis.

Physical and Chemical Parameters Measurements
The intensity of solar irradiance was measured with a realtime solar irradiance monitoring device (EKO, Japan). Salinity and temperature were measured using a salinometer and digital thermometer, respectively. The pH was determined instantly using a pH meter (Orion Star A211, Thermo Scientific), which was calibrated against NIST-traceable pH buffers. Samples for dissolved inorganic carbon (DIC) measurements were poisoned with final concentration of 1 saturated HgCl 2 solution and stored in 40 mL borosilicate glass vials until further analysis. The DIC was determined by acidifying 0.5 mL samples and then measuring CO 2 concentration using a DIC analyzer (AS-C3, Apollo, United States). The other parameters of carbonate system were calculated from DIC and pH NBS using CO2SYS program (Lewis and Wallace, 1998) based on known levels of nutrients, temperature, and salinity.
Samples for nutrient measurements were filtered through 0.45 mm cellulose acetate membrane by suction filter. The filtrate was divided into 2 of 125 mL high density polyethylene bottles: one was stored at −20 • C and measured concentrations of NO 3 − , NO 2 − , PO 4 3− , and SiO 3 2− using an auto-analyzer (AA3, Seal, Germany) at room temperature, and another added with 1 chloroform was stored at 4 • C and measured concentration of NH 4 + with indophenol blue spectrophotometry (Tri-223, Spectrum, China) at 25 • C.

Chlorophyll a Measurement
Water samples of 250-1000 mL were filtered onto GF/F membranes by suction filter with low vacuum pressure (<0.02 MPa), and soaked in pure methanol at 40 • C for 1 h. Then the extracts were centrifuged at 6000 × g for 10 min, and the absorption spectra of supernatants from 400 to 800 nm were measured using a UV-VIS Spectrophotometer (DU800, Beckman, United States). The Chl a concentration was calculated according to the equations from Ritchie (2006).

Viruses and Bacteria Abundance Measurement
Samples filtered by 20 µm mash were collected into 2 mL centrifuge tubes and fixed with glutaraldehyde to a final concentration of 0.5% at room temperature for 15 min, then frozen with liquid nitrogen and stored at −80 • C for further measurement. Viruses and bacteria were counted according to forward scatter, side scatter, and fluorescence intensity after samples were stained with SYBR Green I (Molecular Probe, United States), using the Epics Altra II flow cytometer (Beckman Coulter, United States) and Accuri C6 flow cytometer (Becton and Dickinson, United States), respectively (Chen et al., 2019). For the abundance measurement of both viruses and bacteria, fluorescent beads (Molecular Probes) with a diameter of 1 µm were added as an internal standard.

Phytoplankton Community Structure Determination
Water samples of 500-2000 mL were collected into polyethylene bottles and fixed by adding with 1% lugol's iodine. Then the samples were static placed, and concentrated to 25 mL using siphons within 3 days. The concentrated samples were observed using microscopy and plankton chamber to analyze phytoplankton abundance and diversity based on morphological characteristics.

Measurements of Primary Productivity and Respiratory Carbon Loss
Water samples of 120 mL were collected and transferred immediately into 6 glass scintillation vials (3 vials were used for 12 h inoculation and 3 vials were used for 24 h inoculation) before sunrise. The samples were inoculated with 100 µL of 5 µCi (0.185 MBq) NaH 14 CO 3 solution (ICN Radiochemicals, United States) for 12 h and 24 h with similar levels of sunlight and temperature as in the mesocosms. After 12 and 24 h incubation, the cells were filtered onto glass-fiber filters (25 mm, Whatman GF/F, United States), and the filters were stored at −20 • C. Later on, the filters were thawed and exposed to HCl fumes overnight and dried at 60 • C for 6 h to remove unincorporated labeled carbon, then 3 ml of scintillation cocktail (Hisafe 3, Perkin-Elmer) was added. The incorporated radioactivity was measured by liquid scintillation counting (LS 6500, Beckman Coulter, United States). The respiratory carbon loss during night was calculated as the difference between carbon fixation of 12 h (day-time primary productivity) and 24 h (daily net primary productivity).

Biogenic Silica Measurement
Water samples of 150-500 mL for biogenic silica measurement were collected onto 2 µm polycarbonate membranes, and rinsed with filtered seawater. Each membrane sample was put into 15 mL centrifuge tube, dried and stored in dryers until analysis in the laboratory. The biogenic silica contents were measured by molybdenum-blue-ness spectrophotometry (Brzezinski and Nelson, 1995). Briefly, the samples were placed in a boiling water bath for 40 min, after mixed with 4 mL of 0.2 mol L −1 NaOH. Then the samples were cooled on ice, and neutralized with 1 mL of 1 mol L −1 HCl. The supernatant, obtained by centrifuging samples at 6000 × g for 10 min, were diluted with MilliQ water, and then added with 2 mL of 0.006 mol L −1 ammonium molybdate and stood for 10 min. Finally, the samples were chromogenically reacted with reducing reagent for 3 h. The absorbance of samples was measured at 810 nm using a UV-VIS Spectrophotometer (DU800, Beckman, United States). For the calculation of biogenic silica contents, 0-30 µmol L −1 Na 2 SiF 6 were used as standards.

Statistical Analyses
The statistical analysis was performed using SPSS Statistics, R and Origin. Independent-samples t-test was conducted to check the significant effects of increased pCO 2 at the level of p < 0.1, because the variations were larger within ecological replicates than biological replicates. To further detect the pCO 2 effects on Chl a concentration, the abundance of viruses and bacteria, BSi concentration, primary productivity and night respiration, the generalized additive mixed-effect modeling (GAMM, R packages "mgcv") was applied to analyze the data with the following procedures (Bach et al., 2017). Three GAMM models (Model I: pCO 2 did not affect the dependent variable, Model II: pCO 2 influenced the constant offset in temporal trends, Model III: pCO 2 changed trend shape) were fitted to each independent variable as a function of time, in which pCO 2 and mesocosm number were set as explanatory and random variables, respectively. To ensure model assumptions were satisfied, we considered heteroscedasticity and temporal autocorrelation of residuals, as well as autoregressive structures. The model with the highest coefficient of correlation (R 2 ) was defined as the best that described the pCO 2 effects.

Abiotic Environment of Mesocosms
Throughout the mesocosm experiment, most of the days were sunny, with daily mean PAR (photosynthetically active radiation) ranged 690-1075 µmol photons m −2 s −1 (Supplementary Figure 1). The average salinity was 30.4 , and temperature increased gradually from 21.9 ± 0.1 to 26.1 ± 0.1 • C during the period of 9th April-15th May, 2018 ( Figure 1A).
The differences of pH NBS , pCO 2 and DIC concentration between AC and HC treatments were maintained, respectively, by 0.13-0.31, 230-600 µatm, and 58-130 µmol kg −1 throughout the experiment (p < 0.1, Figures 1B-D and Supplementary Table 1). The pH NBS of HC and AC treatments dropped, respectively, from the initial values of 7.79 ± 0.03 and 8.01 ± 0.01, to 7.70 ± 0.05 and 7.87 ± 0.03 on day 4. Correspondingly, the DIC concentration and pCO 2 increased by 26 µmol kg −1 and 206 µatm in HC treatment, and by 41 µmol kg −1 and 221 µatm in AC treatments, respectively. Then the pH NBS increased to peak values of 7.94 ± 0.13 in HC treatments on day 10 and 8.13 ± 0.03 in AC treatments on day 12, and the corresponding DIC concentration (pCO 2 ) decreased to 1920 ± 52 µmol kg −1 (735 ± 246 µatm) in HC treatments and to 1824 ± 8 µmol kg −1 (427 ± 30 µatm) in AC treatments. Subsequently, the pH NBS decreased by 0.21 units in HC treatments and by 0.18 units in AC treatments on day 28, and the corresponding DIC concentration (pCO 2 ) increased by 36 µmol kg −1 (501 µatm) in HC treatments and by 44 µmol kg −1 (263 µatm) in AC treatments. Afterward, the carbonate system of both HC and AC treatments sustained relatively stable levels till the end of the experiment.
The nutrient concentrations experienced drastic changes in the early stage of the mesocosm experiment, and maintained stable afterward (Figure 2). The initial concentrations of NO 3 − , NO 2 − , NH 4 + , PO 4 3− , and Si(OH) 4 in all mesocosm bags were 6.73-10.65, 0.49-0.79, 6.07-8.71, 0.18-0.47, and 5.42-9.50 µmol kg −1 , respectively. The concentrations of NH 4 + and PO 4 3− decreased rapidly to nearly 0 in both HC and AC treatments on day 2. However, the concentrations of NO 3 − , NO 2 − , and Si(OH) 4 increased at first and then decreased. The NO 3 − concentration increased to 20.62 ± 14.87 µmol kg −1 in HC treatments and to 12.48 ± 2.77 µmol kg −1 in AC treatments on day 2, and then decreased dramatically to 0.10-4.73 µmol kg −1 on day 4. The NO 2 − concentration increased up to day 2 and then maintained relatively stable till day 6, with the values of 1.99 ± 1.68 µmol kg −1 in HC treatments and 1.20 ± 0.38 µmol kg −1 in AC treatments. Subsequently, the NO 2 − concentration decreased to 0.39 ± 0.48 µmol kg −1 in HC and to 0.44 ± 0.61 µmol kg −1 in AC treatments on day 8. The Si(OH) 4 concentration increased to 10.60 ± 2.50 µmol kg −1 in HC treatments and to 12.68 ± 3.03 µmol kg −1 in AC treatments on day 2, then decreased progressively to 1.08 ± 1.00 µmol kg −1 in HC and to 0.64 ± 1.28 µmol kg −1 in AC treatments on day 10. After reaching the very low values, the concentrations of all nutrients were stable till the end of the mesocosm experiment. No significant differences in nutrient concentrations were observed between HC and AC treatments during the whole mesocosm experiment (p > 0.1), except that the average concentrations of both NO 3 − and NO 2 − were higher (p = 0. 289 and 0.095, respectively) in HC treatments than in AC treatments during the initial stage of the experiment.

Chlorophyll a Concentration
As shown in Figure 3, the temporal changes in autotrophic biomass were indicated by Chl a concentration, during the whole mesocosm experiment. The Chl a, with initial concentration of 0.17 µg L −1 , increased to 7.91 ± 4.02 µg L −1 on day 12 in AC treatments and 6.41 ± 2.80 µg L −1 on day 11 in HC treatments. Then it decreased to 0.91 ± 0.31 and to 1.03 ± 0.73 µg L −1 on day 22 in AC and HC treatments, respectively. Thereafter, Chl a concentration maintained stable till day 25 in AC treatments and to day 28 in HC treatments. Subsequently, it increased slightly to 1.34 ± 1.42 µg L −1 in AC treatments and to 1.00 ± 0.75 µg L −1 in HC treatments till the end of the experiment. Based on the ln scale of Chl a concentration, phytoplankton growth kinetics in both AC and HC treatments was classified into four phases: exponential phase (Phase I) from days 0 to 8, stationary phase (Phase II) from days 8 to 14, decline phase (Phase III) from days 12 to 25 (AC treatment) or day 28 (HC treatment), and exponential phase (Phase IV) of another growth cycle from day 25 (AC treatment) or day 28 (HC treatment) to day 31. The average value of Chl a concentration was slightly higher in HC treatments than in AC treatments during FIGURE 1 | Temporal variations of temperature (A), pH NBS (B), DIC concentration (C), and pCO 2 (D) in HC (1000 ppmv) and AC (∼410 ppmv) treatments. The pCO 2 was estimated from the measured pH and DIC concentration using CO2SYS program. Data are means ± SD of replicates for HC treatments and for AC treatments. days 8-16, although the GAMM results showed that HC did not influence the development of Chl a concentration ( Table 1).

Viruses and Bacteria Abundance
The abundance of viruses and bacteria underwent dynamics with different patterns, in both HC and AC treatments (Figure 4). The viruses abundance with initial values of 2.01-4.31 × 10 6 particles mL −1 increased, with the growth of autotrophic plankton, to 1.65 ± 0.20 × 10 8 particles mL −1 in HC treatments on day 10 and to 1.00 × 10 8 ± 0.34 × 10 8 particles mL −1 in AC treatments on day 11, respectively. After reaching the peak, the viruses abundance decreased to 3.22 ± 0.94 × 10 7 particles mL −1 in HC treatments and to 2.70 ± 0.56 × 10 7 particles mL −1 in HC treatments on day 16. Subsequently, the viruses abundance sustained relatively stable values till the end of the experiment in the AC treatments. In the HC treatments, however, it did not change significantly during days 16-22 but increased to 5.46 ± 2.75 × 10 7 particles mL −1 during days 25-31. The GAMM analyses revealed that elevated pCO 2 influenced the viruses trend shape (Table 1 and Supplementary Figure 2), with higher abundance in HC during days 6-14 and days 25-31, relative to AC.
The initial abundance of bacteria was 1.78 ± 0.25 × 10 6 cells mL −1 in HC treatments or 1.85 ± 0.96 × 10 6 cells mL −1 in AC treatments. In the early stage of the experiment, the bacteria abundance increased rapidly to 5.94 ± 0.94 × 10 6 cells mL −1 in HC treatments on day 2, which was higher than that in AC treatments of 4.44 ± 0.54 × 10 6 cells mL −1 (p = 0.059), and then the bacteria abundance decreased till day 4 under both HC and AC conditions. With the exponential growth of autotrophs (days 4-8), the bacteria abundance increased, however, it decreased in the stationary phase of phytoplankton (days 8-12). A constant offset in bacteria abundance existed between AC and HC, with higher bacteria abundance in HC throughout the experimental period, indicated by GAMM analyses results (Table 1 and Supplementary Figure 2).

Phytoplankton Community Composition and Succession
Using microscopy, a total of 34 genera were identified, in which 26 genera belonged to diatoms and 8 genera to dinoflagellates. From all of the mesocosm bags, the dominant genera were Chaetoceros, Thalassiosira, Rhizosolenis, Eucampia, and Prorocentrum, and the dominant species were Chaetoceros curvisetus, Rhizosolenis fragilissima, Eucampia zoodiacu, Cerataulina pelagica, and Guinardia delicatula.
Phytoplankton community structure underwent dynamic succession at phylum levels ( Figure 5). Diatoms accounted for a proportion of above 90% in both treatments during the period of days 0-16, subsequently decreased to 79.2 ± 13.7% in HC treatments on day 20, being higher than that of 52.2 ± 26.9% in AC treatments, although no significant differences were observed (p = 0.138). The relative abundance of diatoms appeared to decrease to 39.8 ± 28.8% in AC on day 28, but maintained stable by 81.8 ± 14.5% in HC treatments, resulting in more significant differences in relative abundance of diatoms between HC and AC treatments (p = 0.054). Correspondingly, dinoflagellates dominated by Prorocentrum minimum and Peridinium conicoides, maintained relatively stable abundance with variations less than 10% in HC treatments, but increased to 60.2 ± 28.8% in AC treatments.
The dynamic changes of community structure were more complicated at genus and species levels than at phylum level (Figure 6 and Supplementary Figure 3). The diatom Guinardia,  The pCO 2 effects were classified into three types (type A: pCO 2 did not affect the dependent variable, type B: pCO 2 influenced the constant offset in temporal trends, and type C: pCO 2 changed trend shape).
dominated by G. delicatula, was the most abundant with relative abundance of 53.4% in terms of total diatom + dinoflagellate species, with cell density of about 9.9 × 10 4 cells mL −1 , after in situ community was inoculated into mesocosm bags. However, it decreased rapidly to average relative abundance of 3.2% (density of 2.6 × 10 3 cells mL −1 ) and to 1.0% (density of 1.6 × 10 3 cells mL −1 ) on day 3 in AC and HC treatments, respectively, and no longer existed from day 6. The diatom Eucampia in both HC and AC treatments, mainly composed of E. zoodiacu, maintained density of approximately 6 × 10 4 -9 × 10 4 cells mL −1 during days 0-16, then dropped to 0.2 × 10 4 -3.7 × 10 4 cells mL −1 on day 20 and maintained stable afterward. The relative abundance of Eucampia increased from initial value of 29% to 61-76% and to 71-75% during days 3-6, and subsequently decreased to 7-16% and to 1-13% on day 10, in HC and AC treatments, respectively. The diatom Chaetoceros grew exponentially from 1.4 × 10 4 cells mL −1 , and reached the peak on day 10 with average concentration of 2.9 × 10 6 cells mL −1 in HC and 3.4 × 10 6 cells mL −1 in AC treatments, which accounted for 65.3 and 71.9% of diatom + dinoflagellate assemblages, respectively. Subsequently, the dominant genus in both HC and AC treatments became the diatom Thalassiosira and Rhizosolenis with the average relative abundance of 53 and 34% on day 16, respectively.  Then the abundance of Thalassiosira and Rhizosolenis decreased progressively. During the late stage of these experiment, HC stimulated the growth of the diatom Nitzschia on day 20, but AC stimulated the growth of the dinoflagellate Propocentrum on day 28. As a consequence, the diatom Nitzschia had significantly higher relative abundance during days 20-28 (p = 0.019), with the dinoflagellate Propocentrum accounting for lower proportion on day 28 (p = 0.078), in HC relative to AC treatments.

Primary Productivity and Night-Period Respiration
The temporal variations of day-time primary productivity, daily net primary productivity and night respiratory carbon loss normalized to water volume showed similar changing patterns with autotrophic biomass levels indicated by Chl a concentration (Figure 7). The day-time primary productivity and daily net primary productivity increased with the growth of phytoplankton, from initial values of 24-77 and 21-63 µg C L −1 , to the peak values of 573 and 392 µg C L −1 on average, respectively, as Chl a concentration reached the maximum. Subsequently, the day-time primary productivity and daily net primary productivity decreased progressively, with the decay of algal blooms, with lowered values of 6.5-204 and 6.3-111 µg C L −1 after day 22. Respiration carbon loss of phytoplankton during night period increased from 15 µg C L −1 on average during days 2-8 to 150 µg C L −1 on average on day 9. Subsequently, the night period respiration of phytoplankton maintained a rate of 140-180 µg C L −1 during the stationary phase of autotrophs, and decreased progressively with the decline of autotrophic blooms. Finally, the respiration stayed 7-9 µg C L −1 during days 22-28, and increased slightly at the end of the mesocosm experiment. The daily net primary productivity and day-time primary productivity normalized to Chl a increased immediately from the initial values of 17.3 ± 4.0 and 26.8 ± 9.6 µg C (µg Chl a) −1 , to 88.7 ± 23.9 and 100.1 ± 30.3 µg C (µg Chl a) −1 on day 6, and then sustained relatively stable levels with the average values of 48.2 and 67.9 µg C (µg Chl a) −1 . The night-period respiration normalized to Chl a, with lower values of 12.3 µg C (µg Chl a) −1 during days 0-9, increased drastically and maintained higher values of 12.3 µg C (µg Chl a) −1 . Neither primary productivity nor night respiration of phytoplankton was influenced significantly by HC (Table 1), normalized to water volume or to Chl a, at any time points examined during the experiment.

Biogenic Silica Contents
Throughout the mesocosm experiment, the differences in temporal development of biogenic silica (BSi) contents were not observed between HC and AC treatments ( Table 1 and Figure 8). With development of phytoplankton blooms, the BSi concentration in the water increased and reached peak values of 11.3 ± 0.5 µmol BSi L −1 on day 12. Correspondingly, the BSi contents normalized to Chl a increased from the initial values of approximately 0.05-2.1 ± 0.8 µmol BSi (µg Chl a) −1 on day 4, and then maintained relative stable till day 14. Afterward, the BSi concentration in the water decreased, while the BSi contents normalized to Chl a increased, with the decay of phytoplankton blooms till the end of the experiment.

DISCUSSION
The results from our mesocosm experiment during the period of late spring and early summer indicate that elevated CO 2 concentration to the level projected for the end of this century, along with natural rise of surface seawater temperature, altered phytoplankton community structure, with impeded succession from diatoms to dinoflagellates and enhanced abundance of freeliving viruses and heterotrophic bacteria. Such an ecological response may lead to significant changes in biogeochemical cycles of carbon and nutrients. The succession from diatoms to dinoflagellates occurred after diatom blooms declined (Figures 5, 6). Under the ambient CO 2 level, the dominant diatom Nitzschia closterium was replaced by FIGURE 7 | The changes of day-time primary productivity (A,D), night period respiration (B,E), and daily net primary productivity (C,F) normalized to water volume (left panels) and to Chl a (right panels) levels, in HC (1000 ppmv) and AC (∼410 ppmv) treatments. Data are means ± SD of replicates for HC treatments and for AC treatments.
the mixotrophic dinoflagellates P. minimum and P. conicoides. The phenomenon of succession from diatoms to dinoflagellates have been commonly observed in eutrophic coastal regions (Spilling et al., 2018), which appears to be due to different physiological strategies between the two groups in coping with environmental changes, including CO 2 and temperature changes. It is most likely that mixotrophic dinoflagellates take over the dominance by taking up organic matters after diatoms depleted the inorganic nutrients (Carlsson and Granéli, 1998). Alternatively, dinoflagellates may possess positive interaction with surrounding bacteria consortia to exchange and/or consume matters (Gong et al., 2017). In addition, the heterotrophic metabolism and flagellate motility may boost with increasing temperature. In the present work, water temperature increased by 4.2 • C from 21.9 ± 0.1 • C at the beginning to 26.1 ± 0.1 • C at the end of the experiment ( Figure 1A). Notably, the competitive capacity of dinoflagellates might become lower in the mesocosms relative to the bay, resulting from low abundance and diversity of dinoflagellates in the bay when the mesocosm experiment was initiated, and from suboptimal growth conditions for dinoflagellates in the mesocosms such as aeration disturbance (Lin et al., 2012;Wohlrab et al., 2020). We also noted that the composition of diatoms changed, with larger cells tending to account for higher proportion (Figure 6, Olenina, 2006;Bach et al., 2019). This observation seemed contradict to the strategy of diatoms in coping with nutrients limitation and increasing temperature, which usually results in smaller cells become the dominant (Finkel et al., 2010;Barton et al., 2013). However, larger diatoms have the advantages of storing nutrient for slower growth (Raven, 1987), and heterotrophic refueling of nutrients may also contribute the increased proportion of larger diatoms.
Under the elevated CO 2 concentration, however, the succession of phytoplankton from diatoms to dinoflagellates appeared to occur, but to much less extent, which was not enough for the latter to replace the former for the dominant group (Figures 5, 6). Higher relative abundance of diatoms under HC conditions was also observed in KOSMOS mesocosm experiment conducted in the Gulmar Fjord, Sweden (Bach et al., 2017), as well as at Shikine Island CO 2 seep, Japan Wada et al., 2021). These may attribute to competitive advantages of CCMs in diatoms, under HC conditions. Increased CO 2 availability associated with ocean acidification can fuel phytoplankton biomass when CO 2 becomes limited (Hansen et al., 2007;Berge et al., 2010;Raven et al., 2020). While increased CO 2 levels projected to the end of this century have been shown to down-regulate CCMs in diatoms (see the review by Gao and Campbell, 2014 and literatures therein) and in dinoflagellates (Rost et al., 2006), and the saved energy from operation of CCMs may be beneficial or harmful to phytoplankton growth under limiting and inhibiting light levels, respectively (Gao et al., 2012). Therefore, difference in CCMs efficiency between the two groups might be responsible for the suppressed succession from diatoms to dinoflagellates under the HC. The CO 2 affinity of ribulose 1,5-diphosphate carboxylase/oxidase (Rubisco) for carboxylation is higher in diatoms than dinoflagellates (Reinfelder, 2011). Dinoflagellates possess a high plasticity toward increased pCO 2 , while diatoms do not ( Van de Waal et al., 2019). Considering the higher CO 2 affinity and lack of CCMs plasticity, diatoms may gain a competitive advantage in DIC uptake under elevated pCO 2 conditions. In addition, OA effects on phytoplankton can be modified by other abiotic factors such as temperature and nutrients, depending on taxa and latitude or weather conditions, and therefore would influence plankton community structure (Sommer et al., 2015;Bach et al., 2016). Diatoms migrated poleward more slowly relative to dinoflagellates, implying that diatoms are more resilient to ocean warming (Chivers et al., 2017), and that diatoms would gain more advantages under OA conditions with increasing temperature. Algal responses to nutrient enrichment may differ at small and intermediate mesocosm volumes but are similar in larger volume mesocosms (>20 L), which are typically used in field (Spivak et al., 2011). Therefore, our experiment with a natural community in >3000 L mesocosm had been performed under different levels of CO 2 could provide basic understanding of the OA impacts on natural plankton communities in eutrophic regions. A previous study showed that plankton communities supported by regenerated nutrients may be more responsive to OA than those having access to higher availability of inorganic nutrients . Under HC conditions, the dominant diatoms Nitzschia spp. have a benefit in nutrient uptake (Qu et al., 2017), and less resources were necessary for diatom inorganic carbon acquisition (Bach et al., 2017). Therefore, elevated pCO 2 along with depletion of nutrients hindered the takeover by dinoflagellates, though a meta-analysis study showed that elevated CO 2 appears to favor dinoflagellates more than diatoms in terms of growth (Seifert et al., 2020). On the other hand, increased bacteria and viruses abundance under HC (Figure 4), along with boosted bacterial activity by high temperature (Spilling et al., 2018), could release more nutrients and CO 2 due to enhanced degradation of organic matters, though the inorganic macro-nutrients did not differ obviously (Figure 2).
Throughout processes of phytoplankton blooms, higher autotrophic biomass was observed in HC relative to AC, but the difference was not significant (Figure 3). The rapid increase in phytoplankton biomass led to the quick declines in nutrients from days 2 to 4. Phytoplankton biomass continued to increase even when the levels of NO 3 − and PO 4 3− were very low. This process may be supported by the degradation of organic matters driven by bacteria, as indicated by increased concentrations of NO x − and Si(OH) 4 , decreased NH 4 + concentration and increased bacteria abundance on day 2 (Figure 2). At the early stage, higher NO x − concentration and higher bacteria abundance in HC relative to AC treatments, indicated that OA stimulated more nutrient refueling by bacterial nitrification to support higher phytoplankton biomass. However, the abundances of viruses and bacteria were increased by HC during the growing and declining phase of phytoplankton, respectively, implying that the enhancement of cell lysis by viruses and the negative interaction of bacteria with diatoms (Vincent and Bowler, 2020) may offset the positive effects of OA on autotrophic biomass. In addition, phytoplankton biomass in HC reached the peak one day earlier than AC, which could be attributed to faster exhaustion of nutrients with higher biomass in HC.
While phytoplankton community structure and abundances of viruses and bacteria responded to the elevated pCO 2 , no significant changes were observed in primary productivity of phytoplankton assemblages over time, either normalized to water volume or to Chl a (Figure 7). High pCO 2 also did not significantly influence the production or degradation of BSi (Figure 8). Changes in the milieu pH could affect intracellular acid-base homeostasis, which in turn influences silicon metabolism, expansion of vesicle expansion and valve formation (Hervé et al., 2012). In this study, the initial pCO 2 of seawater in AC cultures was 610 µatm, which was higher than that equilibrated with ambient air pCO 2 of about 410 ppmv. This indicated that eutrophication contributed strongly to the CO 2 enrichment and acidification, with a pH drop of 0.15 units. In addition, phytoplankton in this region have adapted to fluctuating pH with 0.11-0.18 units of daily pH, due to high rates of day-time photosynthetic CO 2 removal and nighttime respiratory CO 2 release. Non-significant effects of HC on primary productivity and BSi contents may result from high resilience of phytoplankton physiological performance to OA and/or balanced effects of elevated pCO 2 and increased acidic stress (Gao et al., 2019). Although differences in primary productivity or BSi contents were not observed between HC and AC treatments, OA could influence biogeochemical cycles of carbon and nutrients because of varied plankton community structure. BSi contents normalized to Chl a maintained stable during the processes of bloom development, but increased with the decline of phytoplankton blooms, reflecting lower rates of BSi remineralizatin by bacteria relative to phytoplankton mortality rates.
In conclusion, our mesocosm experiment demonstrated a community level response to the elevated CO 2 in the eutrophic bay in the coastal waters of the East China Sea, showing repressed species succession of phytoplankton from diatoms to dinoflagellates and enhanced abundance of viruses and heterotrophic bacteria. This finding mirrors observed reduction in biodiversity near shallow CO 2 seeps (Cigliano et al., 2010;Hall-Spencer and Harvey, 2019;Harvey et al., 2021), since acidification would eventually result in "winners" and "losers" (Gattuso et al., 2015). Although we only found that bacteria increased their abundance, it is most likely that their community structure changed, since OA enhanced nitrification at the early stage. The altered bacterial communities by OA were also shown in the waters surrounding a CO 2 seep (Harvey et al., 2020). Increased viruses abundance suggested that OA may enhance viral activity, as shown previously that elevated CO 2 exacerbates viruses attack to its host, a red-tide microalga (Chen et al., 2015). Integratively, OA may alter phytoplankton succession and enhance viral lysis and bacterial degradation, consequently alter biogeochemical processes in eutrophic coastal ecosystems.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
KG, XL, and GG designed the experiment. TW, XY, DZ, JS, SC, HL, NB, RH, LQ, and XJ preserved normal operation of the mesocosm platform. JS, RH, and SC measured pH, salinity, temperature, primary productivity, and biogenic silica contents. XJ and LQ measured chlorophyll a concentration. ZW and YG provided phytoplankton community structure. YY, RZ, and NJ provided viruses and bacteria abundance. XS provided the concentration of nutrient and dissolved inorganic carbon. RH wrote the manuscript. All authors revised the manuscript.

FUNDING
This study was supported by National Natural Science Foundation of China (41720104005, 41721005, 41890803, and 41911530190) and National Key Research and Development Program of China (2016YFA0601400 and 2016YFA0601300).

ACKNOWLEDGMENTS
We sincerely thank Wenyan Zhao, Xianglan Zeng, and Liting Peng for their kind assistance in operation of mesocosm experiment.