Yunlan Yang
Fei Zhang
Xiaowei Chen
Huifang Li
Nianzhi Jiao
Rui Zhang- 1College of the Environment and Ecology, Xiamen University, Xiamen, China
- 2State Key Laboratory of Marine Environmental Science, Fujian Key Laboratory of Marine Carbon Sequestration, Xiamen University, Xiamen, China
- 3Laboratory of Marine Biology and Ecology, Ministry of Natural Resources, Third Institute of Oceanography, Xiamen, China
- 4College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
Ocean acidification, as one of the major consequences of global climate change, markedly affects multiple ecosystem functions in disparate marine environments from coastal habitats to the deep ocean. Evaluation of the responses of marine microbial community to the increasing partial pressure of CO2 (pCO2) is crucial to explore the microbe-driven biogeochemical processes in the future ocean. In this study, a microcosm incubation of eutrophic coastal water from Xiamen Bay under elevated pCO2 (about 1,000 μatm) and control (ambient air, about 380–410 μatm) conditions was conducted to investigate the effect of ocean acidification on the natural bacterioplankton community. During the 5-day incubation period, the chlorophyll a concentration and bacterioplankton abundance were not significantly affected by increased pCO2. Hierarchical clustering and non-metric multidimensional scaling analysis based on Bray-Curtis similarity among the bacterioplankton community derived from the 16S rRNA genes revealed an inconspicuous impact of elevated pCO2 on the bacterial community. During the incubation period, Proteobacteria, Bacteroidetes, Actinobacteria, Cyanobacteria, and Epsilonbacteraeota were predominant in all microcosms. Despite the distinct temporal variation in the composition of the bacterioplankton community during the experimental period, statistical analyses showed that no significant difference was found on bacterioplankton taxa between elevated pCO2 and control, indicating that the bacterioplankton at the population-level were also insensitive to elevated pCO2. Our results therefore suggest that the bacterioplankton communities in the fluctuating and eutrophic coastal ecosystems appear to be adaptable to the short-term elevated pCO2.
Introduction
Human activities have triggered substantial changes in global climate systems with a preternatural rate over the past two centuries, leading to massive CO2 absorption by the world’s oceans and a reduction in the pH of seawater which is known as ocean acidification (Caldeira and Wickett, 2003; Orr et al., 2005). The partial pressure of atmospheric CO2 (pCO2) has increased by nearly 40% from the preindustrial period to the present day (about 400 μatm) and is predicted to reach approximately 1,000 μatm by the end of this century (Caldeira and Wickett, 2003; IPCC, 2013). Ocean uptake of CO2 changes the equilibrium of the carbonate system, and the continued release of anthropogenic CO2 may lead to another 0.3–0.4units decline in seawater pH globally by 2100 (Orr et al., 2005; IPCC, 2013). As the most complicated and productive ecosystems, coastal oceans generally consist of diversiform but tightly connected aquatic environments, such as rivers, estuaries, tidal wetlands, and sea margins, all of which are strongly influenced by climatic and anthropogenic factors (Cai et al., 2011; Melzner et al., 2012; Wallace et al., 2014). Recent syntheses of the air-sea CO2 fluxes in coastal waters suggest that CO2 uptake in coastal ecosystems has reached to 0.22–0.45 Pg C yr.−1, which is expected to affect the global carbon flux (Cai et al., 2006; Bauer et al., 2013; Regnier et al., 2013). In addition to a large CO2 sink for the atmosphere, ocean acidification in coastal habitats was also reported to be amplified by eutrophication and hypoxia, and these regions might be grimmer by future climate change than previously thought (Cai et al., 2011; Melzner et al., 2012; Wallace et al., 2014). Therefore, the subsequent effects of ocean acidification on coastal life have become one of the most important issues.
Microorganisms exist everywhere and are abundant in density and genetic diversity (Azam and Worden, 2004; Azam and Malfatti, 2007). It has been estimated that microorganisms account for more than two-thirds of marine biomass, despite their tiny size (Bar-On and Milo, 2019). Furthermore, they are key components of the marine food web and play crucial roles in marine ecosystem function and carbon cycling (Azam, 1998; Jiao et al., 2010). The responses of microorganisms to climate change (e.g., ocean acidification) will be pivotal for the marine food web and the biogeochemical cycle in the future ocean. The responses of cyanobacteria, as important primary producers in the ocean, to elevated pCO2 have been investigated in terms of their growth rates, photosynthesis, carbon concentration mechanisms, and cellular affinities for inorganic carbon (Fu et al., 2007; Flombaum et al., 2013; Dutkiewicz et al., 2015; Hutchins and Fu, 2017; Ma and Wang, 2021). With a higher abundance than cyanobacteria by 1–2 orders of magnitude, heterotrophic bacteria play a major role in recycling dissolved organic carbon and nutrients through the microbial loop; however, far less studies were conducted on relationships between heterotrophic bacteria and ocean acidification than phytoplankton (Azam et al., 1983; Azam, 1998; Robinson and Ramaiah, 2011). Previous studies have shown that responses of bacterioplankton communities to elevated pCO2 are diverse and complex and even conflicting. Bacterioplankton communities were found to be susceptible to changes in pCO2 or pH in some studies that were reflected in abundance, diversity, or composition (Krause et al., 2012; Zhang et al., 2013; Xia et al., 2019; Aguayo et al., 2020; Crummett, 2020), whereas other studies reported negligible effects (Roy et al., 2013; Oliver et al., 2014; Wang et al., 2015). Take coastal ecosystems for example, microbial community from coastal ecosystem with naturally low pH (average=7.8) was still sensitive to acidification (Crummett, 2020). However, recent evidence demonstrated that bacterioplankton community from variable coastal ecosystem was relatively stable under elevated pCO2 condition (Allen et al., 2020). The relationships between bacterioplankton communities and elevated pCO2 were more likely to be associated with environmental characteristics, such as nutrients and temperature (Bergen et al., 2016; Sala and Aparicio, 2016; Allen et al., 2020; Wang et al., 2021). Nevertheless, most studies on this topic have been performed in mesotrophic high-latitude regions, especially the Arctic Ocean, while less investigations were conducted in eutrophic low- and middle-latitude coastal ecosystems.
As a subtropical coastal region, Xiamen Bay is characterized by the input of nutrient-rich freshwater from Jiulong River, intrusion of saltwater from the South China Sea, and inflow of artificial wastewater, as well as being affected by intensive human activities (Chen et al., 2013; Cai et al., 2016). The Xiamen coastal ecosystem has showed nutrient-enhanced eutrophication since mid-1990s because of its hydrographical setting (Chen et al., 2013; Cai et al., 2016; Fu et al., 2016). Analyses of the long-term variations in the concentrations of nutrients in Xiamen coastal seawater suggest that dissolved inorganic nitrogen and phosphate have increased by several fold over recent decades (Cai et al., 2016). A time series sampling investigation revealed that the concentration of dissolved organic matter has even in Xiamen coastal ecosystem (Chen et al., 2021b). And estuarine input contributed a mass of dissolved organic matter, including humic-like fluorescent dissolved organic matter, S-containing, and N-containing organic molecules, to the coastal ecosystem (Chen et al., 2021b). In addition to eutrophication, interactions between anthropogenic CO2 emissions and local drivers in coastal ecosystems (such as eutrophication, hypoxia, and biological activities) result in complex regulation of the pH and carbon cycle in coastal waters (Cai et al., 2011; Melzner et al., 2012; Wallace et al., 2014). Analyses of the annual variations in water quality indices in Xiamen Bay from 1986 to 2007 revealed that the mean pH values ranged from 7.86 to 8.21 and continuously decreased (Cai et al., 2016). To investigate the effects of elevated pCO2 on bacterioplankton community at such a changing environment, a 5-day microcosm incubation experiment with eutrophic coastal water from Xiamen Bay under elevated pCO2 and ambient conditions was conducted. During the incubation, the response of bacterioplankton community structure to elevated CO2 and the population size of bacterioplankton were continuously detected.
Materials and Methods
Location and Experimental Setup
The microcosm experiment was carried out in March 2013 using surface seawater collected from a site east of Xiamen, China (24°29′47″N, 118°14′12″E; Figure 1A). Environmental parameters, such as pH, temperature, salinity, and dissolved oxygen, were measured in situ using a YSI Professional Plus multiparameter meter (YSI Incorporated, Yellow Springs, OH, United States). For nutrients analysis, water samples were filtered through a 0.45μm cellulose acetate filter and stored at −20°C. The nitrite (NO2−), nitrate (NO3−), phosphate (PO43−), and silicate (SiO32−) concentrations were analyzed according to colorimetric method with a Technicon AA3 Auto-Analyzer (Bran+Lube, GmbH, Norderstedt, Germany; Dai et al., 2008; Han et al., 2012). To determine microbial abundance, samples (1.98ml for each one) were fixed with a final concentration of 0.5% glutaraldehyde (Sangon Biotech, Shanghai, China) for 15–30min at room temperature (about 25°C), flash-frozen in liquid nitrogen, and then kept at −80°C until analysis (Marie et al., 1997).
Figure 1. Location of sampling station (A) and setup of microcosm experiment (B). The mesocosms were treated with two levels of pCO2, HC (about 1,000 μatm) and control (about 380–410 μatm), and incubated under dark (0 μE m−2 s−1) and light (64.4 μE m−2 s−1) cycles at approximately in situ temperature (about 17.5°C). Map is colored according to the pH values of seawater around Xiamen Island. Map was constructed using Ocean Data View software (version 5.3.0).
The microcosm experimental design is shown in Figure 1B. Briefly, a total of 80L of surface seawater were collected and pre-filtered through a 20μm mesh to remove the large-size fraction. The pre-filtered seawater was then distributed into 20L polycarbonate bottles (Nalgene, United States), with approximately 60–70% of light transmittance, that served as the experimental microcosms (Ha et al., 2013). To simulate atmospheric CO2 concentrations currently and by the end of this century, two levels of pCO2 (ambient air was set as control, about 380–410 μatm; elevated pCO2 was set as HC treatment, about 1,000 μatm) were obtained by adjusting ambient air with CO2 using an enrichment device (Wuhan Ruihua Instrument and Equipment, Wuhan, China). For equilibration of the carbonate system, bubbling was continued throughout the experiment. The pH was monitored using a pH detector (Thermo Scientific, Waltham, MA, United States) during the whole incubation period (Supplementary Figure 1). The microcosms were incubated under a 12h light (64.4 μE m−2 s−1) and 12h dark (0 μE m−2 s−1) cycle at approximately in situ temperature (about 17.5°C), using a light incubator equipped with fluorescent lamps of 400–700nm wavelength (PGX-450B-HM, Ningbo Saifu, China). Two replications were established for the HC treatment and the control. During 5-days of incubation, subsamples were collected from all microcosms every other day to determine the chlorophyll a (Chl a) concentration and the microbial abundance and community structure.
Measurements of Chlorophyll a Concentration and Microbial Abundance
For Chl a determination, about 100ml water samples were filtered onto GF/F filters (25mm, Whatman, Sigma-Aldrich, MO, United States) and extracted overnight at 4°C in absolute methanol. Then, each sample was centrifuged at 5,000× g for 10min to remove particulates, and the absorbance of the supernatant was determined using a UV-VIS spectrophotometer (Beckman Coulter, Brea, CA, United States). The Chl a content was calculated according to the formulae reported by Porra (2002).
Total bacterioplankton (including autotrophic and heterotrophic bacteria), virioplankton, and autotrophic picoplankton were counted with a flow cytometer (Epics Altra II, Beckman Coulter, United States; Marie et al., 1997; Jiao et al., 2002; Brussaard, 2004; Brussaard et al., 2010). After thawing, bacterioplankton samples were diluted in Tris-EDTA buffer (pH=8, Sigma-Aldrich) and then stained with SYBR Green I (10,000× final concentration in DMSO, Molecular Probes, Invitrogen, Carlsbad, CA, United States) for 15min in the dark at room temperature (about 25°C). Samples for virus counting were dyed with SYBR Green I for 10min at 80°C in the dark before analyzed. Abundance of autotrophic picoplankton (picoeukaryote, Synechococcus, and Prochlorococcus) could be directly detected by flow cytometry without dyeing. As an internal standard, 10μl of 1μm diameter fluorescent microspheres (Molecular Probes Inc.) was added to all samples before flow cytometry analysis. The data were obtained and analyzed with EXPOTM32 MultiCOMP software (Beckman Coulter) and FCM Express software (De Novo Software, version 3). Green fluorescence and side scatter were recorded and used as discriminators of bacterioplankton and virioplankton, while autotrophic picoplankton was identified in the plots of side scatter vs. red fluorescence and orange fluorescence vs. red fluorescence.
DNA Extraction, PCR, and Sequencing
Microbial cells were collected from a 1 L sample of each microcosm for bacterial community structure analysis using membrane filtration (0.22μm-pore-size Isopore membrane, Millipore, Billerica, MA, United States). To avoid nucleic acid degradation, samples were flash-frozen in liquid nitrogen and stored at −80°C until DNA extraction. Parallel samples were mixed, and total genomic DNA was extracted using a bacterial DNA extraction kit (Tiangen DP302, Beijing, China) following the manufacturer’s instructions. Amplification, library construction, and sequencing of the extracted DNA were conducted by the Shanghai Personal Biotechnology Co., Ltd. (Shanghai, China). The V4 region of the bacterial 16S rRNA gene was amplified with the primers 520F (5'– GCACCTAAYTGGGYDTAAAGNG–3') and 802R (5'–TACNVGGGTATCTAATCC–3'; Lu et al., 2018; Zhu et al., 2018). The DNA libraries were constructed using a TruSeq Nano DNA LT Library Prep Kit (Illumina, San Diego, CA, United States) following the preparation guide, and the amplicons were sequenced on the Illumina Miseq PE300 Platform.
Sequence Assignment and Data Analysis
All sequence analyses were performed using fast length adjustment of SHort reads version 1.2.7 (FLASH) and quantitative insights into microbial ecology version 1.8.0 (QIIME; Caporaso et al., 2010; Magoc and Salzberg, 2011). A total of 291,315 raw sequences were obtained, and high-quality sequences (quality score≥Q20 and without poly-N strings) longer than 150bp were reserved. The unique operational taxonomic units (OTUs) were clustered at 97% sequence similarity. The OTUs with only one sequence (singleton) were eliminated, and the same number of sequences from each sample were subsampled for further analysis. Classification was carried out using the SILVA database (release 132) with 80% cutoff (Quast et al., 2013). Rarefaction curves based on the identified OTUs were estimated by PAST (version 3.18; Supplementary Figure 2). The Good’s coverage, richness (Chao1), and diversity (Simpson’s and Shannon’s) indexes were calculated using QIIME. Clustering and non-metric multidimensional scaling (NMDS) analysis based on Bray-Curtis similarity of OTUs relative abundance were carried out using PRIMER 6 (version 6.1.16). To determine significant differences in bacterioplankton communities between HC and control and among different sampling days, the analysis of similarity (ANOSIM) was tested using PRIMER 6.
Statistical Analyses
The significance of differences in Chl a content, total bacterioplankton abundance, viral abundance, and Good’s coverage, Chao1, Shannon’s, and Simpson’s indexes between HC and the control were assessed by the paired-samples t-test using PASW statistical software (version 18.0.0). To determine the significance of differences in bacterioplankton communities between HC and the control and among different sampling times, statistical analysis of similarities was conducted using PRIMER 6. Significant differences in bacterioplankton community composition between HC and control at the phylum, class, order, family, and genus levels were detected using the statistical analysis of metagenomic profiles (STAMP) software package (Parks et al., 2014). To explore the relationship between bacterioplankton communities and environmental variables (pH, Chl a, bacterial abundance, and viral abundance), distance-based multivariate regression analysis (DistLM) was carried out using forward selection in Primer 6 with the PERMANOVA+ add on package (version 1.0.6).
Results
Environmental Conditions
The abiotic and biotic characteristics of in situ seawater are shown in Supplementary Table 1. The initial pH value of seawater was 7.55, while the mean pH value of global seawater is estimated to be about 8.1. The temperature of the seawater collection was 17.5°C, and this temperature was maintained during the incubation period. The sampling station was located at a subtropical coastal ecosystem that was affected by saline water from the South China Sea. Therefore, the salinity at this location (31.55) was higher than that at other areas around Xiamen (about 23.86-30.00, unpublished data), and this result was consistent with prior studies (Liu et al., 2017; Wang et al., 2019). The collected seawater had an initial dissolved oxygen concentration of 8.68mgl−1. Similar to the results of other studies, the seawater was eutrophic and the nutrient concentrations of NO2−, NO3−, PO43−, and SiO32− were 7.11, 51.66, 1.38, and 34.40μmoll−1, respectively (Chen et al., 2019; Wang et al., 2019). The abundance of in situ bacterioplankton (1.43±0.07×106 cells ml−1) at the sampling station was consistent with published values for Xiamen Bay (Wang et al., 2019). The abundance of Synechococcus and picoeukaryotes was 2.40±0.15×103 cells ml−1 and 1.09±0.00×104 cells ml−1, respectively. In addition, no Prochlorococcus was detected at the sampling site.
Dynamics of Bacterioplankton Abundance
During the 5-day incubation period, phytoplankton abundance was assessed by the Chl a concentration, which increased from 1.39±0.00μgl−1 to 26.70±3.29μgl−1 in HC and from 0.93±0.00μgl−1 to 31.10±4.27μgl−1 in the control in the first 3-days and declined afterward (Figure 2A), and no significant difference in Chl a between HC and the control was found (t-test, p>0.05). The abundance of total bacterioplankton increased on day 1, dropped to the lowest values of 1.75±0.11×106 cells ml−1 in HC and 1.72±0.22×106 cells ml−1 in the control on day 3, and then slightly increased toward the end of the incubation period (Figure 2B). In addition, bacterial abundance was insensitive to elevated pCO2, although it was slightly higher in HC than in the control on day 1 of incubation. The bacterioplankton grew preferentially on the first day but not keep increasing concomitant with phytoplankton in the following days, and nutrients competition may be one of the reasons for the different growth patterns of bacterioplankton and phytoplankton (Thingstad et al., 2008; Zhang et al., 2017; Huang et al., 2018). During the incubation period, viral abundance in both HC and the control remained nearly constant with a range of 1.16–1.91×107 particles ml−1, and no significant difference between HC and the control was observed (t-test, p>0.05).
Figure 2. Changes in Chl a concentration (A) and total bacterioplankton abundance (B) during incubation under elevated pCO2 (HC, about 1,000 μatm) or ambient pCO2 (Control, about 380–410 μatm) conditions in microcosm experiment. Error bars indicated SD of replicate samples.
Bacterioplankton Community Composition
After screening and quality control, a total of 230,562 high-quality sequences were obtained from all samples. Subsamples with 8,832 sequences from each sample were clustered into 403 to 686 OTUs at the 97% sequence similarity level (Supplementary Table 2). The richness of the bacterioplankton community was found to be invariable in the first 3-days but decreased to 71.33% in HC and 56.53% in the control on day 5 of the incubation period. Although the bacterioplankton community richness did change during the incubation period, it did not significantly differ between HC and the control (t-test, p>0.05; Figure 3). Likewise, bacterioplankton community diversity, as indicated by Shannon’s and Simpson’s indexes, decreased slightly on the last day and was not affected by elevated pCO2 (t-test, p>0.05; Figure 3). The results of the clustering and NMDS analyses showed that the bacterioplankton community composition changed during the incubation and was significantly different on day 5 of the incubation period than on other days (Figure 4). Further analyses revealed that the temporal shift in bacterioplankton community composition was significant (ANOSIM global test: global R=1, p=0.01). However, there was no evidence that elevated pCO2 affect the bacterioplankton community composition (ANOSIM global test: global R=−0.222, p>0.05).
Figure 3. Bacterioplankton community richness (A) and diversity (B,C) in elevated pCO2 (HC, about 1,000 μatm) and ambient pCO2 (Control, about 380–410 μatm) conditions. Changes in Chao1 index, Simpson’s index, and Shannon’s index over time are shown on the left. Boxplots on the right show responses of Chao1 index, Simpson’s index, and Shannon’s index to elevated pCO2 throughout the whole experiment.
Figure 4. Clustering (A) and non-metric multidimensional scaling (NMDS; B) analyses showing relationships between the bacterioplankton communities in elevated pCO2 (HC, about 1,000 μatm, solid circles) and ambient pCO2 (Control, about 380–410 μatm, hollow circles) conditions during microcosm experiment. Similarities were based on Bray-Curtis similarity. Samples were clustered at a similarity level of 70% in NMDS analysis.
Overall, bacterioplankton communities were allocated to 26, 44, 125, 198, and 372 groups at the phylum, class, order, family, and genus levels, respectively (Supplementary Table 3). Bacterioplankton community compositions at the phylum, class, order, family, and genus levels (with relative abundance higher than 1%) are shown in Figure 5. Proteobacteria, which mainly consisted of Alphaproteobacteria and Gammaproteobacteria, were predominant in both HC and the control throughout the incubation period, accounting for about 57.17% of the bacterioplankton community (Figure 5). Bacteroidetes, Actinobacteria, Cyanobacteria, and Epsilonbacteraeota were also abundant in all microcosms and were not affected by elevated pCO2 (Figures 5, 6). The relative abundance of Bacteroidetes increased by about three-fold during the incubation period, while the proportion of Actinobacteria decreased to 16.67% (Supplementary Figure 3). As expected, Cyanobacteria reached to the highest abundance after 3-day incubation and were insensitive to elevated pCO2, consistent with the Chl a concentration (Supplementary Figure 3; Figure 6). Likewise, bacterial groups at finer levels varied over time but did not differ significantly between HC and the control (STAMP analysis, p>0.05; Figures 5, 6). These results suggested that dominant groups constantly changed over time resulting in a prominent temporal shift in the bacterioplankton community. However, contrary to our expectation, the relative abundance of these taxa did not differ significantly between HC and the control (t-test, p>0.05), indicating that the bacterioplankton community from the Xiamen coastal ecosystem remained stable under elevated pCO2 conditions upon a short-term incubation.
Figure 5. Composition of bacterioplankton community at the phylum (A), class (B), order (C), family (D), and genus (E) levels in elevated pCO2 (HC, about 1,000 μatm) and ambient pCO2 (Control, about 380–410 μatm) conditions during incubation. Charts show relative abundance of taxa present at >1%.
Figure 6. STAMP analysis of relative abundance of bacterioplankton communities at the phylum (A), class (B), order (C), family (D), and genus (E) levels between elevated pCO2 (HC, about 1,000 μatm) and ambient pCO2 (Control) conditions.
To identify the potential drivers of changes in the bacterioplankton community, DistLM-forward analysis was carried out to explore the relationships between environmental variables and the bacterioplankton community. However, the bacterioplankton community was not significantly correlated with any of the variables of the microcosms in this study, indicating that these variables could not explain the changes in the bacterioplankton community during the incubation period (Table 1).
Table 1. Relationship between bacterioplankton communities and environmental variables (pH, Chl a, bacterial abundance and viral abundance) in microcosms as determined by distance-based multivariate regression analysis with forward selection (DistLM-forward).
Discussion
Our results showed that both the population size and community structure of bacterioplankton were not significantly affected by elevated pCO2 during 5-days of incubation, indicating that the bacterioplankton community in the coastal Xiamen Bay ecosystem was adaptable to the short-term elevated pCO2. Similar results were also observed in earlier investigations (Grossart et al., 2006; Newbold et al., 2012; Oliver et al., 2014). For instance, elevated pCO2 did not significantly affect the total bacterial cell count distributions in a marine picoplankton community under phytoplankton pre-bloom and post-bloom conditions (Newbold et al., 2012). Bacterioplankton communities were found to be highly resistant to short-term catastrophic pCO2 perturbation in a mesocosm experiment, and no significant differences in community abundance, structure, or composition were observed (Oliver et al., 2014).
Recently, a hypothesis was re-proposed that environmental stability might influence the sensitivity of the bacterioplankton community to climate change (Liu et al., 2010; Joint et al., 2011; Wang et al., 2021). In other words, many bacterioplankton communities are already adapted to changing environments due to long-term exposure to variable environmental conditions and subsequent influence. The evidence to date suggests that coastal communities might be more resistant or flexible under changing environmental conditions than communities in more stable environments such as open ocean gyres (Wang et al., 2021). In addition, many studies have detected minor effects of elevated pCO2 on bacterial abundance, while only few studies showed statistical significant responses which were mainly conducted in oligotrophic oceans, implying that nutrients might be an important influence on the relationship between bacterial population size and elevated pCO2 (Grossart et al., 2006; Allgaier et al., 2008; Newbold et al., 2012; Maas et al., 2013; Bergen et al., 2016; Lin et al., 2018; James et al., 2019; Xia et al., 2019; Allen et al., 2020; Crummett, 2020; Hu et al., 2021). Similarly, significant impacts of elevated pCO2 on bacterioplankton community diversity and composition also have been reported to be associated with nutrient regimes (Roy et al., 2013; Sala and Aparicio, 2016; Allen et al., 2020). By way of illustration, bacterioplankton community composition changed consistently in response to elevated pCO2 at the ultra-oligotrophic center of the South Pacific gyre, while no significant pCO2 treatment effect was found at the mesotrophic fringe of the South Pacific gyre (Allen et al., 2020). In general, therefore, it seems that bacterioplankton communities from coastal ecosystems were more stable in response to elevated pCO2 and these ecosystems were always characterized as rapidly changing and eutrophic.
The Xiamen coastal ecosystem is a typical subtropical coastal ecosystem that is subjected to complex geographical and environmental influences. One of the consequences caused by anthropogenic activities and hydrological factors was the drastic pH fluctuation in Xiamen coastal ecosystem. The highest pH value of Xiamen coastal water has exceeded 8.5, and the lowest was under 7.5 for the last few decades (Cai et al., 2016). In addition, the probability of acid rain in Xiamen reached up to 68.8%, and this might also contribute to the low pH of the seawater.1 Therefore, the bacterioplankton community in the Xiamen coastal region has already experienced the average surface ocean pH predicted to occur at the end of the century or even lower. Previous studies showed that the bacterioplankton community composition in an elevated pCO2 mesocosm would be more conserved through time and resistant to CO2 perturbation (Oliver et al., 2014). In addition, the variable pH in coastal ecosystems should also consider the effects of biological activity although the changes might be contrary to expectation in the future ocean. Biological driven diel fluctuations in pH could reach to 0.3–0.5 pH units in coastal ecosystems and even exceed 0.5 pH units during phytoplankton blooms or red tides (Joint et al., 2011; Hendriks et al., 2015). According to the records of red tide outbreak in the Xiamen coastal ecosystem, a total of 53 red tide events were recorded from 1986 to 2017 (Chen et al., 2021a). All these evidence indicated that the Xiamen coastal seawater has undergo drastic pH fluctuations on daily, seasonal, and even inter-annual scales. The coastal bacterioplankton communities in this area have probably adapted to these changes through processes including physiological acclimation and evolution (Evans and Hofmann, 2012). Moreover, as a result of the complexity and fluctuation in coastal habitats, bacterioplankton communities in coastal ecosystems were suggested to be highly variable reflecting the heterogeneity. And this heterogeneity might play a role in community stability (Zinger et al., 2011; Shade et al., 2012).
Prior studies have also noted the importance of trophic states in the response of the bacterioplankton community to ocean acidification. Bacterioplankton communities were suggested to be more resistant to ocean acidification in nutrient-rich waters (Roy et al., 2013; Sala and Aparicio, 2016; Allen et al., 2020). There was an experimental demonstration of the trophic effect in response of bacterioplankton to elevated pCO2 revealed that more pronounced pH homeostasis genes were aroused to cope with pH stress in oligotrophic marine environments compared with high-nutrient conditions (Bunse et al., 2016). Additionally, the high expression levels of pH homeostasis genes in some bacterial groups are at the expense of growth, and this can ultimately affect the composition and diversity of the bacterioplankton community (Bunse et al., 2016; Allen et al., 2020). However, the energy cost of pH homeostasis expression was not necessary for bacterial cells in eutrophic oceans. The evidence thus far supports the idea that physiological acclimation of the bacterioplankton community to elevated pCO2 is highly possible in eutrophic and highly changeable primitive environments.
Our results differ from those of a previous study on phytoplankton from the Xiamen nearshore, which showed that CO2 enrichment enhanced the relative abundance of Flavobacteria during the early stage of a phytoplankton bloom (Lin et al., 2018). Notably, we used in situ bacterioplankton communities in our study, while the previous study introduced an artificial phytoplankton community into mesocosms system and conducted the incubation for a longer period (Lin et al., 2018). Therefore, one explanation for the differences in results might be that our 5-day incubation was too short to detect the long-term responses of the bacterioplankton community to seawater acidification, since Flavobacteria only showed increased relative abundance in the HC treatments at day 10 in the study of Lin et al. (2018). Another possible explanation is that artificial phytoplankton inoculated in the mesocosms influenced the competitive ability of Flavobacteria group at high pCO2 level. In addition, the possible interference of bacterial community (including Flavobacteria) of the inoculated phytoplankton cultures could not be ruled out although the authors thought natural bacterioplankton was the determiner of responses to different CO2 concentrations.
Conclusion
Bacterioplankton communities response to elevated pCO2 in coastal regions are supposed to be foresight and important. Our results suggest the bacterioplankton community in the coastal region of Xiamen appears to be adaptable to the short-term elevated pCO2 on account of the eutrophic and changeable habitat. To understand the ecological processes and mechanism underlined these phenomena, better experimental setup (e.g., >3 replicates), more comprehensive analysis of relevant environmental parameters (such as dynamics of organic and inorganic nutrients) and including other ecological components (e.g., heterotrophic nanoflagellates and phytoplankton) are required. In addition, given the influences of long-term environmental exposures on microbial phenotypic plasticity, acclimation, and evolutionary adaptation, this study cannot rule out the long-term effects of ocean acidification on coastal bacterioplankton communities. Thus, further experimentation at multiple temporal scales are needed to address issues related to acclimation and adaptation. Considering the diversification of coastal marine ecosystems caused by specific hydrogeological conditions and anthropogenic activities, predicting how coastal bacterioplankton communities will respond to elevated pCO2 requires more investigations in more coastal ecosystems. Overall, the findings of this study contribute to our knowledge of bacterioplankton community responses to ocean acidification in coastal area and highlight the need for further research toward to understanding the long-term effects of ocean acidification on dynamic coastal ecosystems.
Data Availability Statement
The data presented in the study are deposited in the national center for biotechnology information (NCBI) sequence read archive (SRA) repository, accession numbers SRR14766467–SRR14766473 (BioProject accession number PRJNA736025; BioSample accession numbers SAMN19606215–SAMN19606221). The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/sra/PRJNA736025, PRJNA736025.
Author Contributions
RZ and NJ supervised the project and revised manuscript. FZ and YY performed the experiments. YY, FZ, XC, HL, and RZ analyzed data and wrote the manuscript. All authors interpreted the data and gave comments on the manuscript.
Funding
This study was supported by the National Key Research and Development Program of China (2020YFA0608300, 2021YFE0193000), National Natural Science Foundation (41861144018), China Postdoctoral Science Foundation (2019 M662237), the Senior User Project of RV KEXUE (KEXUE2020G10) from Center for Ocean Mega-Science, Chinese Academy of Sciences, and China Scholarship Council.
Conflict of Interest
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.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors thank Wenfang Lin, Jia Sun, Jianning Wang and other colleagues for assistance with in situ sampling and parameters determination. The authors also thank Hong Chen and Yu Wang for their suggestions on data analysis.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.730377/full#supplementary-material.
Footnotes
References
Aguayo, P., Campos, V. L., Henriquez, C., Olivares, F., De Ia Iglesia, R., Ulloa, O., et al. (2020). The influence of pCO2-driven ocean acidification on open ocean bacterial communities during a short-term microcosm experiment in the eastern tropical South Pacific (ETSP) off northern Chile. Microorganisms. 8:1924. doi: 10.3390/microorganisms8121924
Allen, R., Hoffmann, L. J., Law, C. S., and Summerfield, T. C. (2020). Subtle bacterioplankton community responses to elevated CO2 and warming in the oligotrophic South Pacific gyre. Environ. Microbiol. Rep. 12, 377–386. doi: 10.1111/1758-2229.12844
Allgaier, M., Riebesell, U., Vogt, M., Thyrhaug, R., and Grossart, H.-P. (2008). Coupling of heterotrophic bacteria to phytoplankton bloom development at different pCO2 levels: a mesocosm study. Biogeosciences 5, 1007–1022. doi: 10.5194/bg-5-1007-2008
Azam, F. (1998). Microbial control of oceanic carbon flux the plot thickens. Science 280, 694–696. doi: 10.1126/science.280.5364.694
Azam, F., Fenchel, T., Field, J. G., Gray, J. S., Meyer-Reil, L. A., and Thingstad, F. (1983). The ecological role of water-column microbes in the sea. Mar. Ecol. Prog. Ser. 10, 257–263. doi: 10.3354/meps010257
Azam, F., and Malfatti, F. (2007). Microbial structuring of marine ecosystems. Nat. Rev. Microbiol. 5, 782–791. doi: 10.1038/nrmicro1747
Azam, F., and Worden, A. Z. (2004). Microbes, molecules, and marine ecosystems. Science 303, 1622–1624. doi: 10.1126/science.1093892
Bar-On, Y. M., and Milo, R. (2019). The biomass composition of the oceans: A blueprint of our blue planet. Cell 179, 1451–1454. doi: 10.1016/j.cell.2019.11.018
Bauer, J. E., Cai, W. J., Raymond, P. A., Bianchi, T. S., Hopkinson, C. S., and Regnier, P. A. (2013). The changing carbon cycle of the coastal ocean. Nature 504, 61–70. doi: 10.1038/nature12857
Bergen, B., Endres, S., Engel, A., Zark, M., Dittmar, T., Sommer, U., et al. (2016). Acidification and warming affect prominent bacteria in two seasonal phytoplankton bloom mesocosms. Environ. Microbiol. 18, 4579–4595. doi: 10.1111/1462-2920.13549
Brussaard, C. P. D. (2004). Optimization of procedures for counting viruses by flow cytometry. Appl. Environ. Microbiol. 70, 1506–1513. doi: 10.1128/aem.70.3.1506-1513.2004
Brussaard, C. P. D., Payet, J. P., Winter, C., and Weinbauer, M. G. (2010). Quantification of aquatic viruses by flow cytometry. Manual Aquantic Viral Ecology. 11, 102–109. doi: 10.4319/mave.2010.978-0-9845591-0-7.102
Bunse, C., Lundin, D., Karlsson, C. M. G., Akram, N., Vila-Costa, M., Palovaara, J., et al. (2016). Response of marine bacterioplankton pH homeostasis gene expression to elevated CO2. Nat. Clim. Chang. 6, 483–487. doi: 10.1038/nclimate2914
Cai, M., Liu, Y., Chen, K., Huang, D., and Yang, S. (2016). Quantitative analysis of anthropogenic influences on coastal water – A new perspective. Ecol. Indic. 67, 673–683. doi: 10.1016/j.ecolind.2016.03.037
Cai, W. J., Dai, M., and Wang, Y. (2006). Air-sea exchange of carbon dioxide in ocean margins: A province-based synthesis. Geophys. Res. Lett. 33, 1–4. doi: 10.1029/2006GL026219
Cai, W. J., Lohrenz, S. E., Chou, W. C., Zhai, W., Hollibaugh, J. T., Wang, Y., et al. (2011). Acidification of subsurface coastal waters enhanced by eutrophication. Nat. Geosci. 4, 766–770. doi: 10.1038/ngeo1297
Caldeira, K., and Wickett, M. E. (2003). Anthropogenic carbon and ocean pH. Nature 425:365. doi: 10.1038/425365a
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Chen, B., Wang, K., Dong, X., and Lin, H. (2021a). Long-term changes in red tide outbreaks in Xiamen Bay in China from 1986 to 2017. Estuar. Coast. Shelf Sci. 249:107095. doi: 10.1016/j.ecss.2020.107095
Chen, N., Peng, B., Hong, H., Turyaheebwa, N., Cui, S., and Mo, X. (2013). Nutrient enrichment and N:P ratio decline in a coastal bay-river system in Southeast China: The need for a dual nutrient (N and P) management strategy. Ocean Coast. Manage. 81, 7–13. doi: 10.1016/j.ocecoaman.2012.07.013
Chen, Q., Chen, F., Gonsior, M., Li, Y., Wang, Y., He, C., et al. (2021b). Correspondence between DOM molecules and microbial community in a subtropical coastal estuary on a spatiotemporal scale. Environ. Int. 154:106558. doi: 10.1016/j.envint.2021.106558
Chen, X., Ma, R., Yang, Y., Jiao, N., and Zhang, R. (2019). Viral regulation on bacterial community impacted by lysis-lysogeny switch: A microcosm experiment in eutrophic coastal waters. Front. Microbiol. 10:1763. doi: 10.3389/fmicb.2019.01763
Crummett, L. T. (2020). Acidification decreases microbial community diversity in the Salish Sea, a region with naturally high pCO2. PLoS One 15:e0241183. doi: 10.1371/journal.pone.0241183
Dai, M., Wang, L., Guo, X., Zhai, W., Li, Q., He, B., et al. (2008). Nitrification and inorganic nitrogen distribution in a large perturbed river/estuarine system: the Pearl River estuary, China. Biogeosciences 5, 1227–1244. doi: 10.5194/bg-5-1227-2008
Dutkiewicz, S., Morris, J. J., Follows, M. J., Scott, J., Levitan, O., Dyhrman, S. T., et al. (2015). Impact of ocean acidification on the structure of future phytoplankton communities. Nat. Clim. Chang. 5, 1002–1006. doi: 10.1038/nclimate2722
Evans, T. G., and Hofmann, G. E. (2012). Defining the limits of physiological plasticity: how gene expression can assess and predict the consequences of ocean change. Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci. 367, 1733–1745. doi: 10.1098/rstb.2012.0019
Flombaum, P., Gallegos, J. L., Gordillo, R. A., Rincón, J., Zabala, L. L., Jiao, N., et al. (2013). Present and future global distributions of the marine cyanobacteria Prochlorococcus and Synechococcus. Proc. Natl. Acad. Sci. U. S. A. 110, 9824–9829. doi: 10.1073/pnas.1307701110
Fu, F. X., Warner, M. E., Zhang, Y., Feng, Y., and Hutchins, D. A. (2007). Effects of increased temperature and CO2 on photosynthesis, growth, and elemental ratios in marine Synechococcus and Prochlorococcus (cyanobacteria). J. Phycol. 43, 485–496. doi: 10.1111/j.1529-8817.2007.00355.x
Fu, T., Chen, B., Ji, W., Chen, H., Chen, W., Dong, X., et al. (2016). Size structure of phytoplankton community and its response to environmental factors in Xiamen Bay China. Environ. Earth Sci. 75:734. doi: 10.1007/s12665-016-5552-2
Grossart, H.-P., Allgaier, M., Passow, U., and Riebesell, U. (2006). Testing the effect of CO2 concentration on the dynamics of marine heterotrophic. Limnol. Oceanogr. 51, 1–11. doi: 10.4319/lo.2006.51.1.0001
Ha, S.-Y., Joo, H.-M., Kang, S.-H., Ahn, I.-Y., and Shin, K.-H. (2013). Effect of ultraviolet irradiation on the production and composition of fatty acids in plankton in a sub-Antarctic environment. J. Oceanogr. 70, 1–10. doi: 10.1007/s10872-013-0207-3
Han, A., Dai, M., Kao, S.-J., Gan, J., Li, Q., Wang, L., et al. (2012). Nutrient dynamics and biological consumption in a large continental shelf system under the influence of both a river plume and coastal upwelling. Limnol. Oceanogr. 57, 486–502. doi: 10.4319/lo.2012.57.2.048E
Hendriks, I. E., Duarte, C. M., Olsen, Y. S., Steckbauer, A., Ramajo, L., Moore, T. S., et al. (2015). Biological mechanisms supporting adaptation to ocean acidification in coastal ecosystems. Estuar. Coast. Shelf Sci. 152, A1–A8. doi: 10.1016/j.ecss.2014.07.019
Hu, C., Li, X., He, M., Jiang, P., Long, A., and Xu, J. (2021). Effect of ocean acidification on bacterial metabolic activity and community composition in oligotrophic oceans, inferred from short-term bioassays. Front. Microbiol. 12:583982. doi: 10.3389/fmicb.2021.583982
Huang, Y., Liu, X., Laws, E. A., Chen, B., Li, Y., Xie, Y., et al. (2018). Effects of increasing atmospheric CO2 on the marine phytoplankton and bacterial metabolism during a bloom: A coastal mesocosm study. Sci. Total Environ. 633, 618–629. doi: 10.1016/j.scitotenv.2018.03.222
Hutchins, D. A., and Fu, F. X. (2017). Microorganisms and ocean global change. Nat. Microbiol. 2:17058. doi: 10.1038/nmicrobiol.2017.58
James, A. K., Kelly, L. W., Nelson, C. E., Wilbanks, E. G., and Carlson, C. A. (2019). Elevated pCO2 alters marine heterotrophic bacterial community composition and metabolic potential in response to a pulse of phytoplankton organic matter. Environ. Microbiol. 21, 541–556. doi: 10.1111/1462-2920.14484
Jiao, N., Herndl, G. J., Hansell, D. A., Benner, R., Kattner, G., Wilhelm, S. W., et al. (2010). Microbial production of recalcitrant dissolved organic matter: long-term carbon storage in the global ocean. Nat. Rev. Microbiol. 8, 593–599. doi: 10.1038/nrmicro2386
Jiao, N., Yang, Y., Koshikawa, H., and Watanabe, M. (2002). Influence of hydrographic conditions on picoplankton distribution in the East China Sea. Aquat. Microb. Ecol. 30, 37–48. doi: 10.3354/ame030037
Joint, I., Doney, S. C., and Karl, D. M. (2011). Will ocean acidification affect marine microbes? ISME J. 5, 1–7. doi: 10.1038/ismej.2010.79
Krause, E., Wichels, A., Gimenez, L., Lunau, M., Schilhabel, M. B., and Gerdts, G. (2012). Small changes in pH have direct effects on marine bacterial community composition: a microcosm approach. PLoS One 7:e47035. doi: 10.1371/journal.pone.0047035
Lin, X., Huang, R., Li, Y., Li, F., Wu, Y., Hutchins, D. A., et al. (2018). Interactive network configuration maintains bacterioplankton community structure under elevated CO2 in a eutrophic coastal mesocosm experiment. Biogeosciences 15, 551–565. doi: 10.5194/bg-15-551-2018
Liu, J., Weinbauer, M. G., Maier, C., Dai, M., and Gattuso, J. P. (2010). Effect of ocean acidification on microbial diversity and on microbe-driven biogeochemistry and ecosystem functioning. Aquat. Microb. Ecol. 61, 291–305. doi: 10.3354/ame01446
Liu, L., Cai, L., and Zhang, R. (2017). Co-existence of freshwater and marine T4-like myoviruses in a typical subtropical estuary. FEMS Microbiol. Ecol. 93:fix119. doi: 10.1093/femsec/fix119
Lu, Y., Feng, L., Yang, G., Yang, Q., Zhang, X., and Mu, J. (2018). Intensification and microbial pathways of simultaneous nitrification-denitrification in a sequencing batch biofilm reactor for seawater-based saline wastewater treatment. J. Chem. Technol. Biotechnol. 93, 2766–2773. doi: 10.1002/jctb.5652
Ma, J., and Wang, P. (2021). Effects of rising atmospheric CO2 levels on physiological response of cyanobacteria and cyanobacterial bloom development: A review. Sci. Total Environ. 754:141889. doi: 10.1016/j.scitotenv.2020.141889
Maas, E. W., Law, C. S., Hall, J. A., Pickmere, S., Currie, K. I., Chang, F. H., et al. (2013). Effect of ocean acidification on bacterial abundance, activity and diversity in the Ross Sea. Antarctica. Aquat. Microb. Ecol. 70, 1–15. doi: 10.3354/ame01633
Magoc, T., and Salzberg, S. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507
Marie, D., Partensky, F., Jacquet, S., and Vaulot, D. (1997). Enumeration and cell cycle analysis of natural populations of marine picoplankton by flow cytometry using the nucleic acid stain SYBR green I. Appl. Environ. Microbiol. 63, 186–193. doi: 10.1128/aem.63.1.186-193.1997
Melzner, F., Thomsen, J., Koeve, W., Oschlies, A., Gutowska, M. A., Bange, H. W., et al. (2012). Future Ocean acidification will be amplified by hypoxia in coastal habitats. Mar. Biol. 160, 1875–1888. doi: 10.1007/s00227-012-1954-1
Newbold, L. K., Oliver, A. E., Booth, T., Tiwari, B., DeSantis, T., Maguire, M., et al. (2012). The response of marine picoplankton to ocean acidification. Environ. Microbiol. 14, 2293–2307. doi: 10.1111/j.1462-2920.2012.02762.x
Oliver, A. E., Newbold, L. K., Whiteley, A. S., and van der Gast, C. J. (2014). Marine bacterial communities are resistant to elevated carbon dioxide levels. Environ. Microbiol. Rep. 6, 574–582. doi: 10.1111/1758-2229.12159
Orr, J. C., Fabry, V. J., Aumont, O., Bopp, L., Doney, S. C., Feely, R. A., et al. (2005). Anthropogenic Ocean acidification over the twenty-first century and its impact on calcifying organisms. Nature 437, 681–686. doi: 10.1038/nature04095
Parks, D. H., Tyson, G. W., Hugenholtz, P., and Beiko, R. G. (2014). STAMP: statistical analysis of taxonomic and functional profiles. Bioinformatics 30, 3123–3124. doi: 10.1093/bioinformatics/btu494
Porra, R. J. (2002). The chequered history of the development and use of simultaneous equations for the accurate determination of chlorophylls a and b. Photosynthesis Res. 73, 149–156. doi: 10.1023/A:1020470224740
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219
Regnier, P., Friedlingstein, P., Ciais, P., Mackenzie, F. T., Gruber, N., Janssens, I. A., et al. (2013). Anthropogenic perturbation of the carbon fluxes from land to ocean. Nat. Geosci. 6, 597–607. doi: 10.1038/ngeo1830
Robinson, C., and Ramaiah, N. (2011). “Microbial heterotrophic metabolic rates constrain the microbial carbon pump,” in Microbial Carbon Pump in the Ocean. eds. N. Jiao, F. Azam, and S. Sanders (Washington: The American Association for the Advancement of Science), 52–53.
Roy, A. S., Gibbons, S. M., Schunck, H., Owens, S., Caporaso, J. G., Sperling, M., et al. (2013). Ocean acidification shows negligible impacts on high-latitude bacterial community structure in coastal pelagic mesocosms. Biogeosciences 10, 555–566. doi: 10.5194/bg-10-555-2013
Sala, M. M., and Aparicio, F. L. (2016). Contrasting effects of ocean acidification on the microbial food web under different trophic conditions. ICES J. Mar. Sci. 73, 670–679. doi: 10.1093/icesjms/fsv130
Shade, A., Peter, H., Allison, S. D., Baho, D. L., Berga, M., Burgmann, H., et al. (2012). Fundamentals of microbial community resistance and resilience. Front. Microbiol. 3:417. doi: 10.3389/fmicb.2012.00417
Thingstad, T. F., Bellerby, R. G., Bratbak, G., Borsheim, K. Y., Egge, J. K., Heldal, M., et al. (2008). Counterintuitive carbon-to-nutrient coupling in an Arctic pelagic ecosystem. Nature 455, 387–390. doi: 10.1038/nature07235
Wallace, R. B., Baumann, H., Grear, J. S., Aller, R. C., and Gobler, C. J. (2014). Coastal Ocean acidification: The other eutrophication problem. Estuar. Coast. Shelf Sci. 148, 1–13. doi: 10.1016/j.ecss.2014.05.027
Wang, Y., Liu, Y., Wang, J., Luo, T., Zhang, R., Sun, J., et al. (2019). Seasonal dynamics of bacterial communities in the surface seawater around subtropical Xiamen Island, China, as determined by 16S rRNA gene profiling. Mar. Pollut. Bull. 142, 135–144. doi: 10.1016/j.marpolbul.2019.03.035
Wang, Y., Zhang, R., Zheng, Q., Deng, Y., Nostrand, J. D. V., Zhou, J., et al. (2015). Bacterioplankton community resilience to ocean acidification: evidence from microbial network analysis. ICES J. Mar. Sci. 73, 865–875. doi: 10.1093/icesjms/fsv187
Wang, Z., Tsementzi, D., Williams, T. C., Juarez, D. L., Blinebry, S. K., Garcia, N. S., et al. (2021). Environmental stability impacts the differential sensitivity of marine microbiomes to increases in temperature and acidity. ISME J. 15, 19–28. doi: 10.1038/s41396-020-00748-2
Xia, X., Wang, Y., Yang, Y., Luo, T., Van Nostrand, J. D., Zhou, J., et al. (2019). Ocean acidification regulates the activity, community structure and functional potential of heterotrophic bacterioplankton in an oligotrophic gyre. J. Geopgys. Res. 124, 1001–1017. doi: 10.1029/2018jg004707
Zhang, R., Jiang, T., Tian, Y., Xie, S., Zhou, L., Li, Q., et al. (2017). Volcanic ash stimulates growth of marine autotrophic and heterotrophic microorganisms. Geology 45, 679–682. doi: 10.1130/G38833.1
Zhang, R., Xia, X., Lau, S. C. K., Motegi, C., Weinbauer, M. G., and Jiao, N. (2013). Response of bacterioplankton community structure to an artificial gradient of CO2 in the Arctic Ocean. Biogeosciences 10, 3679–3689. doi: 10.5194/bg-10-3679-2013
Zhu, W., Wang, C., Sun, F., Zhao, L., Dou, W., Mao, Z., et al. (2018). Overall bacterial community composition and abundance of nitrifiers and denitrifiers in a typical macrotidal estuary. Mar. Pollut. Bull. 126, 540–548. doi: 10.1016/j.marpolbul.2017.09.062
Keywords: elevated pCO2, bacterioplankton community, abundance, community composition, eutrophic coastal ecosystem
Citation: Yang Y, Zhang F, Chen X, Li H, Jiao N and Zhang R (2021) Insignificant Response of Bacterioplankton Community to Elevated pCO2 During a Short-Term Microcosm Experiment in a Subtropical Eutrophic Coastal Ecosystem. Front. Microbiol. 12:730377. doi: 10.3389/fmicb.2021.730377
Edited by:
Feixue Fu, University of Southern California, United StatesReviewed by:
Juntian Xu, Jiangsu Ocean University, ChinaYantao Liang, Ocean University of China, China
Copyright © 2021 Yang, Zhang, Chen, Li, Jiao and Zhang. 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(s) 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: Nianzhi Jiao, amlhb0B4bXUuZWR1LmNu; Rui Zhang, cnVpemhhbmdAeG11LmVkdS5jbg==
†These authors have contributed equally to this work and share first authorship