Next Article in Journal
Squalamine and Its Aminosterol Derivatives: Overview of Biological Effects and Mechanisms of Action of Compounds with Multiple Therapeutic Applications
Previous Article in Journal
Landscape Composition and Soil Physical–Chemical Properties Drive the Assemblages of Bacteria and Fungi in Conventional Vegetable Fields
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Microorganisms
Volume 10
Issue 6
10.3390/microorganisms10061204
microorganisms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Variations in Marine Bacterial and Archaeal Communities during an Ulva prolifera Green Tide in Coastal Qingdao Areas

by
Guihua Zhao
1,
Hui He
1,*,
Hualong Wang
1,
Yantao Liang
1,
Cui Guo
1,
Hongbing Shao
1,
Yong Jiang
1 and
Min Wang
1,2,3,*
1
College of Marine Life Sciences, Institute of Evolution and Marine Biodiversity, Frontiers Science Center for Deep Ocean Multispheres and Earth System, Ocean University of China, Qingdao 266003, China
2
The Affiliated Hospital of Qingdao University, Qingdao 266000, China
3
OUC-UMT Joint Academic Centre for Marine Studies, Qingdao 266003, China
*
Authors to whom correspondence should be addressed.
Microorganisms 2022, 10(6), 1204; https://doi.org/10.3390/microorganisms10061204
Submission received: 13 May 2022 / Revised: 9 June 2022 / Accepted: 9 June 2022 / Published: 13 June 2022
(This article belongs to the Section Environmental Microbiology)
Browse Figures
Versions Notes

Abstract

:
Green tides caused by Ulva prolifera occur annually in the Yellow Sea, potentially influencing the marine microorganisms. Here, we focused on the variations in marine bacterial and archaeal communities during an U. prolifera green tide in coastal Qingdao areas with Illumina high-throughput sequencing analysis. Our results revealed that the diversity and structure of bacterial and archaeal communities, as well as the organization and structure of microbial co-occurrence networks, varied during the green tide. The decline phase may be favorable to the bacterial and archaeal diversity and richness. The bacterial community, as well as the archaeal community, showed clear variations between the outbreak and decline phases. A simpler and less connected microbial co-occurrence network was observed during the outbreak phase compared with the decline phase. Flavobacteriales and Rhodobacterales separately dominated the bacterial community during the outbreak and decline phase, and Marine Group II (MGII) dominated the archaeal community during the green tide. Combined with microbial co-occurrence network analysis, Flavobacteriales, Rhodobacterales and MGII may be important organisms during the green tide. Temperature, chlorophyll a content and salinity may have an important impact on the variations in bacterial and archaeal communities during the green tide.

1. Introduction

Green tide is a kind of ecological disaster caused by macroalgae mainly including genera Ulva, Chaetomorpha and Cladophora. Since 2007, green tide primarily driven by Ulva prolifera has occurred annually in the Yellow Sea, which is regarded as the largest green tide in the world [1,2]. U. prolifera generally floats on the surface of seawater, posing a serious threat to coastal tourism, fishery and environment [3,4]. The outbreak of U. prolifera green tide affects not only the nutrient contents of seawater but also the interactions among marine organisms, and further influences the marine biogeochemical cycles [5,6,7,8]. During the decomposing process of U. prolifera, it can consume dissolved oxygen resulting in water hypoxia or even anoxia, and also release hydrogen sulfide into seawater, which is hazardous to marine organisms [9,10]. In addition, decomposing U. prolifera can release nutrients into seawater, causing a brief eutrophication and possibly leading to secondary disasters such as jellyfish blooms [6,11].
Due to extremely high diversity and various ecological functions, microorganisms that are involved in biogeochemical cycles, primary production and microbial loop are the keystone of marine ecosystems [12,13,14]. Microorganisms, on the other hand, are regarded as sensitive indicators in characterizing the environmental changes [15]. Algal blooms have crucial effects on the microbial community that may differ depending on causative species, bloom phases and environmental factors [16,17]. For example, Actinobacteria and Betaproteobacteria dominated the bacterial community during a cyanobacteria bloom, while Flavobacteria and Alphaproteobacteria were the major taxa during a diatom bloom [16,18]. Zhou et al. [19] reported that Gammaproteobacteria and Bacteroidetes dominated the bacterial community during the initial phase of a dinoflagellate bloom, while Alphaproteobacteria, Cyanobacteria and Actinobacteria were the major bacterial taxa during the onset and termination phases. Variations in the archaeal community also occurred during a dinoflagellate bloom, with methanogens and ammonia-oxidizing archaea being the main taxa during the early and late phases, respectively [19].
It is reported that U. prolifera may settle in seawater southeast of the Shandong Peninsula, particularly along the coastal Qingdao areas [20]. Variations in bacterial community during the U. prolifera green tide have been reported successively [8,21,22,23,24,25]. For example, Zhang et al. [21] found that the bacterial community varied during an U. prolifera green tide, and the abundance of nitrogen-fixing bacteria dramatically increased during the outbreak phase. U. prolifera green tide in coastal Qingdao areas led to a reduction in bacterial diversity and abundance, and some functional groups, such as sulfate-reducing bacteria and Cytophaga-Flavobacter-Bacteroides, were also affected by the U. prolifera green tide [8]. In recent years, microbial interactions have been considered to be vital drivers in determining the algal blooms [26]. Co-occurrence network analysis has been used in a variety of environments including seawater, sediments and soils [27,28,29], and is an effective tool for the reflection of interactions between microorganisms, which is critical to the assembly and stability of microbial community [30,31,32]. By means of co-occurrence network analysis, Zhou et al. [33] found that time-shifted interdependencies were prevalent in the bacterioplankton community, and Rickettsiaceae, Saprospriaceae, Winogradskyella, Flavobacteriaceae, Chitinophagaceae and NS 11–12 were the core taxa during a dinoflagellate bloom. However, to our knowledge, microbial interactions during an U. prolifera green tide derived from the co-occurrence network analysis have not been reported yet. In this study, variations in marine bacterial and archaeal communities during an U. prolifera green tide in coastal Qingdao areas were explored through the Illumina high-throughput sequencing method based on the 16S rRNA gene. Our primary goals were: (i) to explore the variations in bacterial and archaeal communities during an U. prolifera green tide, (ii) to identify the shift of microbial interactions using co-occurrence network analysis during an U. prolifera green tide, and (iii) to discuss the potential vital environmental factors regulating the variations in bacterial and archaeal communities during an U. prolifera green tide.

2. Materials and Methods

2.1. Sampling

In this study, we chose three stations (XG, ZQ and MD) in coastal Qingdao areas to study the variations in marine bacterial and archaeal communities during an U. prolifera green tide (Figure 1). Among three sampling stations, XG (36.07° E, 120.30° N) was less likely impacted by the green tide, whereas ZQ (36.06° E, 120.31° N) and MD (36.05° E, 120.42° N) were greatly influenced by the green tide. Twenty liters of surface seawater at each station was collected from 13 June to 5 September 2019 (Table 1). After collection, the seawater was filtered through an 800-mesh plankton net to remove large zooplankton, then filtered through a 0.22 μm pore-size polycarbonate membrane. The obtained filters were immediately frozen in liquid nitrogen and stored at −80 °C until further molecular analysis.
Temperature, salinity, dissolved oxygen content (DO) and pH were recorded in situ using a YSI proplus multi-parameter water quality analyzer (YSI, Yellow Springs, OH, USA). Chlorophyll a (chla) content was determined using the spectrophotometric method after extraction using acetone. To determine the nutrient concentrations, the seawater was filtered through a 0.45 μm pore-size GF/F membrane. Then, concentrations of nitrate (NO3−), nitrite (NO2−), ammonium (NH4+) and phosphate (PO43−) were measured with a QuAAtro nutrient auto analyzer (Seal Analytical Ltd., King’s Lynn, UK).

2.2. DNA Extraction and Illumina High-Throughput Sequencing

Genomic DNA was extracted with FastDNA Spin kit for Soil (MP Biomedicals, Solon, OH, USA). The obtained DNA was checked through 1% agarose gel electrophoresis analysis, and quantified with NanoDrop 2000 UV-Vis spectrophotometer (Thermo Scientific, Waltham, MA, USA). The hypervariable regions (V3–V4) of bacterial and archaeal 16S rRNA genes were amplified with the primer pairs 338F (5′-ACTC CTA CGG GAG GCA GCA G-3′) and 806R (5′-GGA CTA CHV GGG TWT CTA AT-3′), 524F10extF (5′-TGY CAG CCG CCG CGG TAA-3′) and Arch958RmodR (5′-YCC GGC GTT GAV TCC AAT T-3′), respectively [34,35]. The amplification procedures were as follows: initial denaturation at 95 °C for 3 min, followed by 27 cycles of denaturation at 95 °C for 30 s, annealing at 55 °C for 30 s and extension at 72 °C for 45 s, and a final extension at 72 °C for 10 min. After purification and quantification, the amplicons were pooled in equimolar, then sequenced on Miseq PE300 platform (Illumina, San Diego, CA, USA) at Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China). The raw reads were submitted to National Centre for Biotechnology Information (NCBI) Sequence Read Archive (SRA) database under accession numbers PRJNA732997, PRJNA739445, PRJNA741737 and PRJNA742371.

2.3. Sequence Processing

Through Fast Length Adjustment of Short reads (FLASH, version 1.2.11, Magoč and Salzberg, Baltimore, MD, USA) [36], raw reads which meet all the following criteria: (i) overlap sequences longer than 10 bp, and (ii) mismatch ratio of overlap greater than 0.2, were merged. The merged sequences were quality-filtered using Quantitative Insights into Microbial Ecology (QIIME, version 1.9.1, Caporaso et al., Boulder, CO, USA) [37]. Briefly, sequences were discarded if they met any of the following criteria: (i) an average quality score below 20, (ii) shorter than 50 bp, and (iii) containing any ambiguous base. Then, operational taxonomic units (OTUs) were clustered with a 97% sequence similarity cutoff using UPARSE (version 7.1, Edgar, CA, USA) [38]. Chimeric sequences were identified and removed at the same time. The predominant sequence in each OTU was chosen as the representative sequence, and annotated with the Silva database (Silva SSU128) using a confidence threshold of 70%.

2.4. Statistical Analysis

To illustrate the alpha-diversity of bacterial and archaeal communities, Chao1 index, Shannon index and Good’s coverage were calculated using QIIME. Principal component analysis (PCA) at the OTU level was applied to evaluate the similarities or differences in community structure between the outbreak and decline phases. The top 15 taxa at the family level were compared between the outbreak and decline phases through Wilcoxon rank-sum test. Redundancy analysis (RDA) was employed to explain how environmental factors influenced the bacterial and archaeal communities during the green tide. Correlations of Chao1 index, Shannon index, and relative abundance of dominant taxa with chla content were estimated with SPSS (version 26, IBM Corporation, New York, NY, USA). One-way analysis of variance (ANOVA) was carried out to determine whether there were statistically significant differences among stations or phases.
Co-occurrence network analysis was conducted at the out level during the green tide, the outbreak phase and the decline phase, respectively. OTUs occurring in less than 50% of samples were filtered from total sequences. The network was performed with the package “fdrtool” and “igraph” in R statistical software (version 4.0.4, R Core Team, Vienna, Austria). Only the OTUs with statistically significant values (p < 0.01 and Q value < 0.05) and the Spearman’s coefficient > |0.6| were put into the further analysis. The co-occurrence network was visualized in Gephi (version 0.9.2, Bastian et al., Paris, France) [39]. The topology of network, including degree, modularity, average path length and average clustering coefficient, were also calculated with Gephi.

3. Results

3.1. Variations in Environmental Factors during the Green Tide

In this study, a total of 29 surface seawater samples were collected during the green tide. Variations in environmental factors in coastal Qingdao areas during the green tide were shown in Table 2. Over the sampling period, temperature, salinity, pH and DO ranged from 18.60–26.50 °C, 29.28–30.39, 7.58–8.37 and 3.92–12.86 mg/L, respectively. Chla contents ranged from 0.65–4.86 μg/L. NH4+ concentrations varied from 3.25–28.70 μmol/L, while NO2− and NO3− concentrations were relatively low with values of 0.04–0.65 μmol/L and 0.59–9.57 μmol/L, respectively. PO43− concentrations varied from 0.10–4.67 μmol/L.
Based on the chla content, all samples were divided into two phases, that is, the outbreak phase (13 June to 22 June 2019) and the decline phase (6 August to 5 September 2019). Temperature (p < 0.01), chla content (p < 0.01), NO2− concentration (p < 0.01) and NO3− concentration (p < 0.01) varied significantly between the outbreak and decline phases.

3.2. Variations in Bacterial and Archaeal Richness and Diversity during the Green Tide

In total, 1,590,865 and 1,270,196 high-quality bacterial and archaeal 16S rRNA gene sequences were generated in this study, with an average length of 415 bp and 428 bp, respectively. All high-quality bacterial and archaeal sequences separately yielded 1088 bacterial OTUs and 163 archaeal OTUs. Good’s coverage was greater than 99.51% in all cases, demonstrating that the sequences we obtained sufficiently covered the majority of bacterial and archaeal taxa in this study.
For the bacterial community, Chao1 index and Shannon index were in the range of 539.73–949.26 and 2.31–4.87, respectively (Table 3). Compared with the bacterial community, lower archaeal Chao1 index and lower archaeal Shannon index were observed, ranging from 56.00–123.00 and 1.61–2.95, respectively. The average value of bacterial Chao1 index and the average value of bacterial Shannon index were greater during the decline phase, and the average value of archaeal Chao 1 index and the average value of archaeal Shannon index showed similar trends. Thus, it is speculated that the decline phase is beneficial to the richness and diversity of bacterial and archaeal communities. In addition, we also observed that the diversity of bacterial community differed significantly among stations (p < 0.05), whereas the diversity of archaeal community varied markedly between the outbreak and decline phases (p < 0.05).

3.3. Variations in Bacterial and Archaeal Community Structure during the Green Tide

Thirty bacterial phyla and five archaeal phyla were observed during the green tide. Here, the bacterial and archaeal community composition was analyzed at the order level. Rhodobacterales (with a relative abundance ranging from 8.07–77.47%) and Flavobacteriales (4.93–32.76%) dominated the bacterial community across all samples (Figure 2a). The relative abundance of Flavobacteriales was significantly greater during the outbreak phase compared with the decline phase (p < 0.05). Rhodobacterales exhibited greater abundance during the decline phase, but not to an obvious extent (p > 0.05). Cellvibrionales (0.55–9.96%), Synechococcales (0.08–48.33%), SAR86 clade (0.01–13.11%), Oceanospirillales (0.39–13.72%), Alteromonadales (0.27–20.81%), Pseudomonadales (0.01–20.04%), Bacteroidales (0.04–29.39%) and Bacillales (0.00–12.84%) also occupied greater proportions in the bacterial community during the green tide. Cellvibrionales (p > 0.05), Synechococcales (p > 0.05), Alteromonadales (p < 0.05), Pseudomonadales (p > 0.05) and Bacillales (p > 0.05) were more abundant during the outbreak phase; in contrast, SAR86 clade (p > 0.05), Oceanospirillales (p > 0.05) and Bacteroidales (p > 0.05) exhibited greater abundance during the decline phase.
The predominant archaeal orders comprised of Marine Group II (MGII, 2.95–98.90%), and to a lesser extent, Nitrosopumilales (0.95–74.36%), unclassified Bathyarchaeia (0.03–6.31%), Nitrososphaerales (0.00–13.57%), Methanosarcinales (0.00–6.28%) and unclassified Thermoplasmata (0.00–5.04%, Figure 2b). It is noteworthy that MGII dominated the archaeal community during the green tide (except for MD0823 and MD0829), and differed markedly between the outbreak and decline phases (p < 0.01). Other dominant archaeal taxa, including Nitrosopumilales, Nitrososphaerales, Methanosarcinales and unclassified Thermoplasmata, presented greater values during the decline phase, but this was non-significant (p > 0.05).
Clear variations in the bacterial community structure, as well as the archaeal community structure, were observed between the outbreak and decline phases (Figure 3). Based on the analysis of similarities (ANOSIM), both the bacterial and archaeal communities revealed significant differences between the outbreak and decline phases (p < 0.01). To identify the differences in community compositions between the outbreak and decline phases, the top 15 taxa at the family level were analyzed using the Wilcoxon rank-sum test. The variations in the bacterial community during the green tide were caused by increases in Rhodobacteraceae, unclassified SAR86 clade, Cryomorphaceae, Ilumatobacteraceae, Actinomarinaceae, SAR116 clade, AEGEAN-169 Marine Group, Sphingomonadaceae and Nitrincolaceae, and decreases in Flavobacteriaceae, Cyanobiaceae, Halieaceae, Alteromonadaceae, Moraxellaceae and Cycloclasticaceae (Figure 4a). Decreases in unclassified MGII, unclassified Bathyarchaeia, Halomicrobiaceae, Methanosaetaceae, Methanocorpusculaceae, Haloferacaceae, unclassified Lokiarchaeia and Halococcaceae, and increases in Nitrosopumilaceae, Nitrososphaeraceae, Methanosarcinaceae, unclassified Thermoplasmata, Methanobacteriaceae, unclassified Woesearchaeia and Methanomicrobiaceae led to the variations in the archaeal community during the green tide (Figure 4b).

3.4. Correlation Analysis of Environmental Factors with Bacterial and Archaeal Communities

RDA was employed to explore the correlation of environmental factors with bacterial and archaeal communities. Temperature (p < 0.01, F = 7.29, 999 Monte Carlo permutations) and NO2− concentration (p < 0.05, F = 2.08, 999 Monte Carlo permutations) contributed most to the variations in bacterial community during the green tide (Figure 5a). The archaeal community during the green tide was strongly associated with temperature (p < 0.01, F = 8.77, 999 Monte Carlo permutations), pH (p < 0.01, F = 6.30, 999 Monte Carlo permutations), NO2− concentration (p < 0.01, F = 4.04, 999 Monte Carlo permutations) and salinity (p < 0.05, F = 2.33, 999 Monte Carlo permutations) (Figure 5d). Meanwhile, it was observed that the bacterial community during the outbreak phase was significantly related to salinity (p < 0.01, F = 3.44, 999 Monte Carlo permutations) and PO43− concentration (p < 0.01, F = 2.31, 999 Monte Carlo permutations) (Figure 5b), while temperature (p < 0.01, F = 16.86, 999 Monte Carlo permutations) exhibited marked correlations with archaeal community during the outbreak phase (Figure 5e). It is noteworthy that during the decline phase, only the archaeal community was found to be significantly impacted by pH (p < 0.01, F = 5.74, 999 Monte Carlo permutations) (Figure 5f).
As previously mentioned, chla content is a vital indicator of U. prolifera green tide. In order to illustrate how chla content affected the bacterial and archaeal communities, relationships of Chao1 index, Shannon index, relative abundance of dominant taxa with chla content were explored through Spearman correlation analysis (Figure 6, Table 4). It is demonstrated that chla content (r = 0.638, p < 0.01) plays a crucial role in the relative abundance of Flavobacteriales during the green tide. During the outbreak phase, the Shannon index of bacterial community (r = 0.857, p < 0.01) and the relative abundance of Flavobacteriales (r = 0.762, p < 0.05) was markedly associated with chla content. Nevertheless, chla content had no significant effect on the Chao1 index, Shannon index and relative abundance of MGII of archaeal community during any phases (p > 0.05).

3.5. Microbial Co-Occurrence Network Analysis during the Green Tide

Co-occurrence network was constructed to explore the microbial interactions during the green tide based on the Spearman correlation. Simultaneously, in order to elucidate the distinctions of microbial interactions between phases, we also constructed two separated co-occurrence networks during the outbreak and decline phases. Most of the correlations in the network of green tide were positive, suggesting that mutualism was the main type of interaction among microorganisms during the green tide (Figure 7a). The modularity index of network was 0.464, implying a modular structure during the green tide (Figure 7d). The nodes with relatively more abundance in each module were specific to a particular phase. For example, the majority of OTUs in module 1 and 3 showed greater relative abundance during the outbreak phase, whereas most of OTUs in module 4 had greater relative abundance during the decline phase.
The nodes and edges in the network of the outbreak phase (681 and 3497) were less than those in the network of the decline phase (809 and 6825), suggesting a more complex network during the decline phase (Figure 7b,c). The modularity index of network was 0.577 and 0.585 during the outbreak and the decline phase, respectively, suggesting a more modular structure during the decline phase (Figure 7e,f). In this study, the module with proportion to the whole network greater than 5% was recognized as a major module. During the outbreak phase, the network was attributed to six major modules, named module 1, 2, 3, 4, 5 and 6 accounting for 22.47%, 15.57%, 12.63%, 11.16%, 7.64% and 7.05% of the network, respectively. The network of the decline phase was comprised of five major modules, module 1, 2, 3, 4 and 5, which separately accounted for 31.64%, 21.38%, 12.98%, 11.00% and 10.75% of the network during the decline phase.
In this study, key species should meet the following criteria: (i) have a high degree ranking in the top 20%, and (ii) have a low betweenness centrality ranking in the bottom 20% among the above high degree nodes [40]. During the outbreak and the decline phase, 27 and 32 key species were identified, respectively, which differed between phases. For instance, most of the key species during the outbreak phase belonged to order Flavobacteriales within class Bacteroidia, whereas order Rhodobacterales within class Alphaproteobacteria and order MGII within class Thermoplasmata were observed as key species during the decline phase. Hence, Flavobacteriales may play a major part in maintaining the stability of microbial community during the outbreak phase, while Rhodobacterales and MGII may be crucial for the persistence of the microbial community during the decline phase.

4. Discussion

In this study, we explored the variations in marine bacterial and archaeal communities during an U. prolifera green tide in coastal Qingdao areas through Illumina high-throughput sequencing analysis. It was concluded that the diversity and structure of bacterial and archaeal communities, as well as the organization and structure of microbial co-occurrence networks, have changed throughout the green tide.
Notable differences on the richness and diversity of both bacterial and archaeal communities were observed between the outbreak and decline phases. In this study, the bacterial and archaeal communities presented greater richness and diversity during the decline phase (Table 3), implying that the decline phase may be beneficial to the bacterial and archaeal richness and diversity. Greater richness and diversity during the decline phase might mean a more complex microbial community at this phase, which was concordant with our results obtained from the microbial co-occurrence network analysis. During the decline phase, decomposing U. prolifera released large quantities of organic matter into surrounding seawater. Meanwhile, inorganic nutrient (NH4+, NO2−, NO3− and PO43−) concentrations were also observed to be greater during the decline phase compared with the outbreak phase (Table 2). Large amounts of organic matter and inorganic nutrient concentrations during the decline phase might explain the reason for the greater richness and diversity of bacterial and archaeal communities during this time.
Clear variation in bacterial community structure was observed during the green tide. Overall, Flavobacteriales and Rhodobacterales dominated the bacterial community during the outbreak and decline phase, respectively (Figure 2a), indicating that they play important roles during the green tide. Due to the rapid response to the enrichment of organic matter, Flavobacteriales was regarded as an essential component of the bacterial community during blooms [18,19,41]. Similar to our results, the observations of field studies and incubation experiments in laboratory all exhibited the dominant position of Flavobacteriales in the bacterial community during the outbreak phase of green tide [22,24,25]. Flavobacteirales were also observed to be more abundant at the late period of the outbreak phase, which may be related to their algicidal activity [42]. Rhodobacterales are related to the degradation of organic matter [22], and play a prominent part in various ecological processes, particularly sulfur oxidation and dimethylsulfoniopropionate demethylation [43]. In addition, Rhodobacterales have the capability to reduce the nutrient concentrations and hydrogen sulfide contents [43,44,45], which may in turn improve the water quality during the green tide.
It is worth noting that Alteromonadaceae exhibited a greater relative abundance during the outbreak phase (Figure 2a), which was in agreement with a previous report on the bacterial community during the bloom caused by Prorocentrum donghaiense [46]. Recent studies demonstrated that a certain genus within family Alteromonadaceae, such as Alteromonas, plays a major part in total bacterial production during blooms, owing to its rapid response to the increasing of organic matter [47]. As one of the first reported algicidal bacteria, Alteromonas can secrete algicidal substances to kill certain algae or algal dissolved cells [48,49]. Hence, Alteromonadaceae, together with Flavobacteriales, might secrete algicidal substances to kill U. prolifera and thus dissolve cells directly or indirectly, which may speed up the coming of decline phase [50].
During the U. prolifera green tide, the archaeal community was less complex compared with the bacterial community. MGII was the predominant archaeal taxon during the green tide (Figure 2b), which is consistent with the results of archaeal community composition during a marine dinoflagellate bloom [19]. MGII is the most common archaeal taxon across surface seawater, and it plays a crucial role in carbon and nitrogen cycles, as well as in the attachment and utilization of particulate organic matter in marine ecosystems [51,52,53,54]. Furthermore, we also observed a remarkable increase in the relative abundance of Nitrosopumilales during the decline phase (Figure 2b). Nitrosopumilales, formerly known as Marine Group I (MGI), is an important participant in ammonia oxidation and carbon fixation, and also one of the major contributors to marine primary productivity [55,56,57,58]. The lower abundance of Nitrosopumilus during the outbreak phase possibly implied a competition for nitrogen nutrients between rapidly growing U. prolifera and archaea (ammonia oxidation process) during this phase. Ammonia oxidation, coupled with denitrification, can remove nitrogen pollution [59], thus greater relative abundance of Nitrosopumilales during the decline phase may prevent the accumulation of toxic ammonium, reduce the eutrophic level and maintain the nitrogen balance in our studied areas.
As an important indicator of phytoplankton biomass and eutrophication level, chla content was a vital factor affecting the bacterial community during the green tide, mainly through its positive impact on the relative abundance of Flavobacterales (r = 0.638, p < 0.01) during the green tide, and Shannon index (r = 0.857, p < 0.01) and relative abundance of Flavobacterales (r = 0.762, p < 0.05) during the outbreak phase (Figure 6, Table 4). Nevertheless, chla content did not display any remarkable influences on the archaeal community during any phases, indicating that the bacterial community, rather than the archaeal community, was more sensitive to the outbreak of U. prolifera. Temperature may be the principal shaping factor of bacterial and archaeal communities during the green tide, which was shown previously [60,61,62]. For example, Bergen and his colleagues identified that temperature was the major force in structuring the bacterial community during both a diatom bloom in autumn and a mixed phytoplankton bloom in summer [62]; Lucas et al. [63] pointed out that temperature was the main driving factor for the short-term successions of bacterioplankton community during phytoplankton blooms in the German Bight. In this study, salinity also showed the potential to impact the variations in microbial community during the green tide, mainly affecting the archaea. Salinity was one of the most significant environmental factors impacting the distribution of MGII archaea [54]. The importance of salinity for the archaeal community may be due to the dominant position of MGII in this study.
In addition, complex microbial interactions elucidated by the co-occurrence network also play important roles in regulating the variations in bacterial and archaeal communities during blooms [26], while few studies have paid attention to the microbial interactions during the U. prolifera green tide so far. Here, two independent co-occurrence networks during the outbreak and decline phase were established to uncover the differences in microbial interactions between phases, and we found that the organization and structure of microbial co-occurrence networks varied during the green tide. Positive correlations manifest mutualistic interactions, while negative correlations may suggest competition or predation among microorganisms [27,64]. Complex positive or negative correlations among microorganisms were displayed in this study, and mutualistic interactions dominated among the microbial interactions during both the outbreak and decline phases (Figure 7b,c). The proportion of positive correlation during the decline phase (87.44%) was greater than that during the outbreak phase (75.86%), suggesting a greater synergic relationship among microorganisms, which may be in favor of the ecological restoration during the decline phase. The nodes and edges in the network of the outbreak phase were less than those in the network of decline phase, which is in line with the above-mentioned lower microbial richness and diversity during the outbreak phase, indicating that the outbreak of U. prolifera may reduce the interactions among microorganisms, as well as the complexity of network structure. In addition, compared with the network of the decline phase, the network of the outbreak phase exhibited lower average degree (4.454 versus 10.523), lower graph density (0.015 versus 0.021), lower clustering coefficient (0.386 versus 0.399) and longer path length (5.082 versus 3.938), indicating that there was a less connected microbial co-occurrence network during the outbreak phase, which may arise from a stronger competition of limited nutrients during this phase [65].
Key species play an exceptionally vital role in maintaining the stability of microbial community [66]. In this study, key species during the outbreak phase were affiliated to Flavobacteriales, while Rhodobacterales and MGII were recognized as the main players during the decline phase, suggesting that Flavobacteriales, Rhodobacterales and MGII may display vital roles during the green tide, and the lack of these organisms may cause disintegration of the structure of microbial co-occurrence networks. As a common bacterial group, Flavobacteriales generally degrade organic matter to obtain energy for its rapid growth and also own algicidal activity, which may speed up the arrival of the decline phase [42,50,67]. Rhodobacterales and MGII are known to participate in diverse biogeochemical cycles, particularly of sulfur, carbon and nitrogen [43,52]. Variations in key species might mean that the microbial potential ecological functions may change during the green tide, and further study is warranted.

5. Conclusions

In this study, variations in marine bacterial and archaeal communities during an U. prolifera green tide were examined in coastal Qingdao areas. The results revealed that the diversity and structure of bacterial and archaeal communities, as well as the organization and structure of microbial co-occurrence networks, varied during the green tide. The decline phase benefits the bacterial and archaeal richness and diversity. The bacterial community, as well as the archaeal community, exhibited clear variations between the outbreak and decline phases. Moreover, distinct differences in the microbial interactions were also observed between phases, that is, a simpler and less connected microbial co-occurrence network existed during the outbreak phase compared with the decline phase. Flavobacteriales, Rhodobacterales and MGII may be pivotal organisms playing irreplaceable roles during the green tide. In addition, temperature, chlorophyll a content and salinity may have an important influence on the variations in bacterial and archaeal community during the green tide.

Author Contributions

Conceptualization, H.H., H.W., Y.L., C.G., H.S., Y.J. and M.W.; methodology, G.Z., H.H., Y.L. and C.G.; formal analysis, G.Z.; investigation, G.Z.; data curation, G.Z.; writing—original draft preparation, G.Z.; writing—review and editing, H.H. and H.W.; supervision, H.H., H.W., Y.L., C.G., H.S., Y.J. and M.W.; funding acquisition, H.H. and M.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (41806131, 41976117, 42120104006, 42176111, 41676178) and Fundamental Research Funds for the Central Universities (202172002).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw reads were submitted to National Centre for Biotechnology Information (NCBI) Sequence Read Archive (SRA) database under accession numbers PRJNA732997, PRJNA739445, PRJNA741737 and PRJNA742371.

Acknowledgments

We are grateful to Ding Zhang, Jianhua Sun and Guiliang Zheng for their assistance in the measurement of environmental factors. We also appreciate the computing resources provided by IEMB-1, a high-performance computation cluster operated by the Institute of Evolution and Marine Biodiversity, Ocean University of China.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhao, J.; Jiang, P.; Liu, Z.; Wei, W.; Lin, H.; Li, F.; Wang, J.; Qin, S. The Yellow Sea green tides were dominated by one species, Ulva (Enteromorpha) prolifera, from 2007 to 2011. Chin. Sci. Bull. 2013, 58, 2298–2302. [Google Scholar] [CrossRef] [Green Version]
  2. Cao, Y.; Wu, Y.; Fang, Z.; Cui, X.; Liang, J.; Song, X. Spatiotemporal patterns and morphological characteristics of Ulva prolifera distribution in the Yellow Sea, China in 2016–2018. Remote Sens. 2019, 11, 445. [Google Scholar] [CrossRef] [Green Version]
  3. Zhang, Y.; He, P.; Li, H.; Li, G.; Liu, J.; Jiao, F.; Zhang, J.; Huo, Y.; Shi, X.; Su, R.; et al. Ulva prolifera green-tide outbreaks and their environmental impact in the Yellow Sea, China. Natl. Sci. Rev. 2019, 6, 825–838. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Liu, X.; Wang, Z. A review of the green tides in the Yellow Sea, China. Mar. Environ. Res. 2016, 119, 189–196. [Google Scholar] [CrossRef]
  5. Song, X.; Shi, Y.; Liu, A.; Xing, H.; Jiang, H.; Wang, W.; Zhang, L. The impact of green tide on the phytoplankton community in Yellow Sea. Appl. Mech. Mater. 2013, 260–261, 1130–1137. [Google Scholar] [CrossRef]
  6. Zhang, T.; Wang, X. Release and microbial degradation of dissolved organic matter (DOM) from the macroalgae Ulva prolifera. Mar. Pollut. Bull. 2017, 125, 192–198. [Google Scholar] [CrossRef]
  7. Wang, W.; Zhang, G.; Sun, X.; Zhang, F.; Zhang, X. Temporal variability in zooplankton community in the western Yellow Sea and its possible links to green tides. PeerJ 2019, 7, e6641. [Google Scholar] [CrossRef]
  8. Qu, T.; Zhao, X.; Hao, Y.; Zhong, Y.; Guan, C.; Hou, C.; Tang, X.; Wang, Y. Ecological effects of Ulva prolifera green tide on bacterial community structure in Qingdao offshore environment. Chemosphere 2020, 244, 125477. [Google Scholar] [CrossRef]
  9. O’Boyle, S.; McDermott, G.; Silke, J.; Cusack, C. Potential impact of an exceptional bloom of Karenia mikimotoi on dissolved oxygen levels in waters off western Ireland. Harmful Algae 2016, 53, 77–85. [Google Scholar] [CrossRef]
  10. Achá, D.; Guédron, S.; Amouroux, D.; Point, D.; Lazzaro, X.; Fernandez, P.E.; Sarret, G. Algal bloom exacerbates hydrogen sulfide and methylmercury contamination in the emblematic high-altitude Lake Titicaca. Geosciences 2018, 8, 438. [Google Scholar] [CrossRef] [Green Version]
  11. Wang, C.; Yu, R.C.; Zhou, M.J. Effects of the decomposing green macroalgae Ulva (Enteromorpha) prolifera on the growth of four red-tide species. Harmful Algae 2012, 16, 12–19. [Google Scholar] [CrossRef]
  12. Ducklow, H.W.; Purdie, D.A.; Williams, P.J.L.; Davies, J.M. Bacterioplankton: A sink for carbon in a coastal marine plankton community. Science 1986, 232, 865–867. [Google Scholar] [CrossRef] [PubMed]
  13. Hoppe, H.G.; Gocke, K.; Koppe, R.; Begler, C. Bacterial growth and primary production along a north-south transect of the Atlantic Ocean. Nature 2002, 416, 168–171. [Google Scholar] [CrossRef] [PubMed]
  14. Seymour, J.R.; Stocker, M.R. Resource patch formation and exploitation throughout the marine microbial food web. Am. Nat. 2009, 173, E15–E29. [Google Scholar] [CrossRef] [Green Version]
  15. Paerl, H.W.; Dyble, J.; Moisander, P.H.; Noble, R.T.; Piehler, M.F.; Pinckney, J.L.; Steppe, T.F.; Twomey, L.; Valdes, L.M. Microbial indicators of aquatic ecosystem change: Current applications to eutrophication studies. FEMS Microbiol. Ecol. 2003, 46, 233–246. [Google Scholar] [CrossRef]
  16. Su, X.; Steinman, A.D.; Tang, X.; Xue, Q.; Zhao, Y.; Xie, L. Response of bacterial communities to cyanobacterial harmful algal blooms in Lake Taihu, China. Harmful Algae 2017, 68, 168–177. [Google Scholar] [CrossRef]
  17. Zhang, H.; Hou, F.; Xie, W.; Wang, K.; Zhou, X.; Zhang, D.; Zhu, X. Interaction and assembly processes of abundant and rare microbial communities during a diatom bloom process. Mar. Microb. Ecol. 2020, 22, 1707–1719. [Google Scholar] [CrossRef]
  18. Klindworth, A.; Mann, A.J.; Huang, S.; Wichels, A.; Quast, C.; Waldmann, J.; Teeling, H.; Glöckner, F.O. Diversity and activity of marine bacterioplankton during a diatom bloom in the North Sea assessed by total RNA and pyrotag sequencing. Mar. Genom. 2014, 18, 185–192. [Google Scholar] [CrossRef]
  19. Zhou, J.; Richlen, M.L.; Sehein, T.R.; Kulis, D.M.; Anderson, D.M.; Cai, Z. Microbial community structure and associations during a marine dinoflagellate bloom. Front. Microbiol. 2018, 9, 1201. [Google Scholar] [CrossRef]
  20. Geng, H.; Yu, R.; Zhang, Q.; Yan, T.; Kong, F.; Zhou, M. Tracing the settlement region of massive floating green algae in the Yellow Sea. J. Oceanol. Limnol. 2019, 37, 1555–1565. [Google Scholar] [CrossRef]
  21. Zhang, X.; Song, Y.; Liu, D.; Keesing, J.K.; Gong, J. Macroalgal blooms favor heterotrophic diazotrophic bacteria in nitrogen-rich and phosphorus-limited coastal surface waters in the Yellow Sea. Estuar. Coast. Shelf Sci. 2015, 163, 75–81. [Google Scholar] [CrossRef]
  22. Lin, G.; Sun, F.; Wang, C.; Zhang, L.; Zhang, X. Assessment of the effect of Enteromorpha prolifera on bacterial community structures in aquaculture environment. PLoS ONE 2017, 12, e0179792. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Chen, J.; Li, H.; Zhang, Z.; He, C.; Shi, Q.; Jiao, N.; Zhang, Y. DOC dynamics and bacterial community succession during long-term degradation of Ulva prolifera and their implications for the legacy effect of green tides on refractory DOC pool in seawater. Water Res. 2020, 185, 116268. [Google Scholar] [CrossRef] [PubMed]
  24. Sun, F.; Wang, C.; Chen, H.; Zheng, Z. Metagenomic analysis of the effect of Enteromorpha prolifera on microbial community and function in aquaculture environment. Curr. Microbiol. 2020, 77, 816–825. [Google Scholar] [CrossRef]
  25. Liang, J.; Liu, J.; Zhan, Y.; Zhou, S.; Xue, C.X.; Sun, C.; Lin, Y.; Luo, C.; Wang, X.; Zhang, X.H. Succession of marine bacteria in response to Ulva prolifera-derived dissolved organic matter. Environ. Int. 2021, 155, 106687. [Google Scholar] [CrossRef]
  26. Needham, D.M.; Fuhrman, J.A. Pronounced daily succession of phytoplankton, archaea and bacteria following a spring bloom. Nat. Microbiol. 2016, 1, 16005. [Google Scholar] [CrossRef]
  27. Lupatini, M.; Suleiman, A.K.A.; Jacques, R.J.S.; Antoniolli, Z.I.; de Siqueira Ferreira, A.; Kuramae, E.E.; Roesch, L.F.W. Network topology reveals high connectance levels and few key microbial genera within soils. Front. Environ. Sci. 2014, 2, 10. [Google Scholar] [CrossRef] [Green Version]
  28. Fuhrman, J.A.; Cram, J.A.; Needham, D.M. Marine microbial community dynamics and their ecological interpretation. Nat. Rev. Microbiol. 2015, 13, 133–146. [Google Scholar] [CrossRef]
  29. Dai, T.; Zhang, Y.; Ning, D.; Su, Z.; Tang, Y.; Huang, B.; Mu, Q.; Wen, D. Dynamics of sediment microbial functional capacity and community interaction networks in an urbanized coastal estuary. Front. Microbiol. 2018, 9, 2731. [Google Scholar] [CrossRef]
  30. Chaffron, S.; Rehrauer, H.; Pernthaler, J.; von Mering, C. A global network of coexisting microbes from environmental and whole-genome sequence data. Genome Res. 2010, 20, 947–959. [Google Scholar] [CrossRef] [Green Version]
  31. Konopka, A.; Lindemann, S.; Fredrickson, J. Dynamics in microbial communities: Unraveling mechanisms to identify principles. ISME J. 2015, 9, 1488–1495. [Google Scholar] [CrossRef] [PubMed]
  32. Hou, F.; Zhang, H.; Xie, W.; Zhou, X.; Zhu, X.; Zhang, D. Co-occurrence patterns and assembly processes of microeukaryotic communities in an early-spring diatom bloom. Sci. Total Environ. 2020, 711, 134624. [Google Scholar] [CrossRef] [PubMed]
  33. Zhou, J.; Lao, Y.M.; Song, J.T.; Jin, H.; Zhu, J.M.; Cai, Z.H. Temporal heterogeneity of microbial communities and metabolic activities during a natural algal bloom. Water Res. 2020, 183, 116020. [Google Scholar] [CrossRef] [PubMed]
  34. Liu, C.; Li, H.; Zhang, Y.; Si, D.; Chen, Q. Evolution of microbial community along with increasing solid concentration during high-solids anaerobic digestion of sewage sludge. Bioresour. Technol. 2016, 216, 87–94. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Xu, N.; Tan, G.; Wang, H.; Gai, X. Effect of biochar additions to soil on nitrogen leaching, microbial biomass and bacterial community structure. Eur. J. Soil Biol. 2016, 74, 1–8. [Google Scholar] [CrossRef]
  36. Magoč, T.; Salzberg, S.L. FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics 2011, 27, 2957–2963. [Google Scholar] [CrossRef] [PubMed]
  37. Caporaso, J.G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F.D.; Costello, E.K.; Fierer, N.; Peña, A.G.; Goodrich, J.K.; Gordon, J.I. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 2010, 7, 335–336. [Google Scholar] [CrossRef] [Green Version]
  38. Edgar, R.C. UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 2013, 10, 996–998. [Google Scholar] [CrossRef]
  39. Bastian, M.; Heymann, S.; Jacomy, M. Gephi: An open source software for exploring and manipulating networks. In Proceedings of the International AAAI Conference on Weblogs and Social Media, San Jose, CA, USA, 17–20 May 2009. [Google Scholar]
  40. Zhang, Y.H.; Yao, P.; Sun, C.; Li, S.Z.; Shi, X.C.; Zhang, X.H.; Liu, J.W. Vertical diversity and association pattern of total, abundant and rare microbial communities in deep-sea sediments. Mol. Ecol. 2021, 30, 2800–2816. [Google Scholar] [CrossRef]
  41. Simon, H.; Lipsewers, Y.A.; Giebel, H.A.; Wiltshire, K.H.; Simon, M. Temperature effects on aggregation during a spring diatom bloom. Limnol. Oceanogr. 2014, 59, 2089–2100. [Google Scholar] [CrossRef]
  42. Skerratt, J.H.; Bowman, J.P.; Hallegraeff, G.; James, S.; Nichols, P.D. Algicidal bacteria associated with blooms of a toxic dinoflagellate in a temperate Australian estuary. Mar. Ecol. Prog. Ser. 2002, 244, 1–15. [Google Scholar] [CrossRef] [Green Version]
  43. Brinkhoff, T.; Giebel, H.A.; Simon, M. Diversity, ecology, and genomics of the Roseobacter clade: A short overview. Arch. Microbiol. 2008, 189, 531–539. [Google Scholar] [CrossRef] [PubMed]
  44. Lenk, S.; Moraru, C.; Hahnke, S.; Arnds, J.; Richter, M.; Kube, M.; Reinhardt, R.; Brinkhoff, T.; Harder, J.; Amann, R.; et al. Roseobacter clade bacteria are abundant in coastal sediments and encode a novel combination of sulfur oxidation genes. ISME J. 2012, 6, 2178–2187. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Voget, S.; Wemheuer, B.; Brinkhoff, T.; Vollmers, J.; Dietrich, S.; Giebel, H.A.; Beardsley, C.; Sardemann, C.; Bakenhus, I.; Billerbeck, S.; et al. Adaptation of an abundant Roseobacter RCA organism to pelagic systems revealed by genomic and transcriptomic analyses. ISME J. 2015, 9, 371–384. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Li, D.X.; Zhang, H.; Chen, X.H.; Xie, Z.X.; Zhang, Y.; Zhang, S.F.; Lin, L.; Chen, F.; Wang, D.Z. Metaproteomics reveals major microbial players and their metabolic activities during the blooming period of a marine dinoflagellate Prorocentrum donghaiense. Environ. Microbiol. 2018, 20, 632–644. [Google Scholar] [CrossRef]
  47. Tada, Y.; Taniguchi, A.; Nagao, I.; Miki, T.; Uematsu, M.; Tsuda, A.; Hamasaki, K. Differing growth responses of major phylogenetic groups of marine bacteria to natural phytoplankton blooms in the western north Pacific Ocean. Appl. Environ. Microbiol. 2011, 77, 4055–4065. [Google Scholar] [CrossRef] [Green Version]
  48. Imai, I.; Ishida, Y.; Sakaguchi, K.; Hata, Y. Algicidal marine bacteria isolated from northern Hiroshima Bay, Japan. Fish. Sci. 1995, 61, 628–636. [Google Scholar] [CrossRef] [Green Version]
  49. Umetsu, S.; Kanda, M.; Imai, I.; Sakai, R.; Fujita, M.J. Questiomycins, algicidal compounds produced by the marine bacterium Alteromonas sp. D and their production cue. Molecules 2019, 24, 4522. [Google Scholar] [CrossRef] [Green Version]
  50. Mayali, X.; Azamn, F. Algicidal bacteria in the sea and their impact on algal blooms. J. Eukaryot. Microbiol. 2004, 51, 139–144. [Google Scholar] [CrossRef]
  51. Orsi, W.D.; Smith, J.M.; Wilcox, H.M.; Swalwell, J.E.; Carini, P.; Worden, A.Z.; Santoro, A.E. Ecophysiology of uncultivated marine euryarchaea is linked to particulate organic matter. ISME J. 2015, 9, 1747–1763. [Google Scholar] [CrossRef] [Green Version]
  52. Zhang, C.L.; Xie, W.; Martin-Cuadrado, A.B.; Rodriguez-Valera, F. Marine Group II archaea, potentially important players in the global ocean carbon cycle. Front. Microbiol. 2015, 6, 1108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Liu, H.; Zhang, C.L.; Yang, C.; Chen, S.; Cao, Z.; Zhang, Z.; Tian, J. Marine Group II dominates planktonic archaea in water column of the northeastern South China Sea. Front. Microbiol. 2017, 8, 1098. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Xie, W.; Luo, H.; Murugapiran, S.K.; Dodsworth, J.A.; Chen, S.; Sun, Y.; Hedlund, B.P.; Wang, P.; Fang, H.; Deng, M.; et al. Localized high abundance of Marine Group II archaea in the subtropical Pearl River Estuary: Implications for their niche adaptation. Environ. Microbiol. 2018, 20, 734–754. [Google Scholar] [CrossRef]
  55. Christman, H.; Campbell, E.; Meeks, J. Global transcription profiles of the nitrogen stress response resulting in heterocyst or hormogonium development in Nostoc punctiforme. J. Bacteriol. 2011, 193, 6874–6886. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Morono, Y.; Terada, N.; Nishizawa, M.; Ito, M.; Hillion, F.; Takahata, N.; Sano, Y.; Inagaki, F. Carbon and nitrogen assimilation in deep subseafloor microbial cells. Proc. Natl. Acad. Sci. USA 2011, 108, 18295–18300. [Google Scholar] [CrossRef] [Green Version]
  57. Stahl, D.A.; de la Torre, J.R. Physiology and diversity of ammonia-oxidizing archaea. Annu. Rev. Microbiol. 2012, 66, 83–101. [Google Scholar] [CrossRef]
  58. Offre, P.; Spang, A.; Schleper, C. Archaea in biogeochemical cycles. Annu. Rev. Microbiol. 2013, 67, 437–457. [Google Scholar] [CrossRef] [Green Version]
  59. Bothe, H.; Jost, G.; Schloter, M.; Ward, B.B.; Witzel, K.P. Molecular analysis of ammonia oxidation and denitrification in natural environments. FEMS Microbiol. Rev. 2000, 24, 673–690. [Google Scholar] [CrossRef]
  60. Niu, Y.; Shen, H.; Chen, J.; Xie, P.; Yang, X.; Tao, M.; Ma, Z.; Qi, M. Phytoplankton community succession shaping bacterioplankton community composition in Lake Taihu, China. Water Res. 2011, 45, 4169–4182. [Google Scholar] [CrossRef]
  61. Xing, P.; Li, H.; Liu, Q.; Zhang, J. Composition of the archaeal community involved in methane production during the decomposition of Microcystis blooms in the laboratory. Can. J. Microbiol. 2012, 58, 1153–1158. [Google Scholar] [CrossRef]
  62. Bergen, B.; Endres, S.; Engel, A.; Zark, M.; Dittmar, T.; Sommer, U.; Jürgens, K. Acidification and warming affect prominent bacteria in two seasonal phytoplankton bloom mesocosms. Environ. Microbiol. 2016, 18, 4579–4595. [Google Scholar] [CrossRef] [PubMed]
  63. Lucas, J.; Wichels, A.; Teeling, H.; Chafee, M.; Scharfe, M.; Gerdts, G. Annual dynamics of North Sea bacterioplankton: Seasonal variability superimposes short-term variation. FEMS Microbiol. Ecol. 2015, 91, fiv099. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Steele, J.A.; Countway, P.D.; Xia, L.; Vigil, P.D.; Beman, J.M.; Kim, D.Y.; Chow, C.T.; Sachdeva, R.; Jones, A.C.; Schwalbach, M.S.; et al. Marine bacterial, archaeal and protistan association networks reveal ecological linkages. ISME J. 2011, 5, 1414–1425. [Google Scholar] [CrossRef] [PubMed]
  65. Wang, H.; Chen, F.; Zhang, C.; Wang, M.; Kan, J. Estuarine gradients dictate spatiotemporal variations of microbiome networks in the Chesapeake Bay. Environ. Microbiome 2021, 16, 22. [Google Scholar] [CrossRef]
  66. Herren, C.M.; McMahon, K.D. Keystone taxa predict compositional change in microbial communities. Environ. Microbiol. 2018, 20, 2207–2217. [Google Scholar] [CrossRef] [Green Version]
  67. Mayali, X.; Franks, P.J.S.; Azam, F. Cultivation and ecosystem role of a marine Roseobacter clade-affiliated cluster bacterium. Appl. Environ. Microbiol. 2008, 74, 2595–2603. [Google Scholar] [CrossRef] [Green Version]
Microorganisms 10 01204 g001 550
Figure 1. Locations of sampling stations in coastal Qingdao areas during the green tide.
Figure 1. Locations of sampling stations in coastal Qingdao areas during the green tide.
Microorganisms 10 01204 g001
Microorganisms 10 01204 g002 550
Figure 2. Relative abundance of dominant bacterial (a) and archaeal (b) taxa at the order level during the green tide.
Figure 2. Relative abundance of dominant bacterial (a) and archaeal (b) taxa at the order level during the green tide.
Microorganisms 10 01204 g002
Microorganisms 10 01204 g003 550
Figure 3. Principal component analysis (PCA) of bacterial (a) and archaeal (b) communities at the OTU level during the green tide.
Figure 3. Principal component analysis (PCA) of bacterial (a) and archaeal (b) communities at the OTU level during the green tide.
Microorganisms 10 01204 g003
Microorganisms 10 01204 g004 550
Figure 4. Wilcoxon rank-sum test of the top 15 bacterial (a) and archaeal (b) taxa at the family level during the green tide. * p < 0.05, ** p < 0.01.
Figure 4. Wilcoxon rank-sum test of the top 15 bacterial (a) and archaeal (b) taxa at the family level during the green tide. * p < 0.05, ** p < 0.01.
Microorganisms 10 01204 g004
Microorganisms 10 01204 g005 550
Figure 5. RDA analysis of environmental factors with bacterial community during the green tide (a), the outbreak (b) and decline (c) phase, and with archaeal community during the green tide (d), the outbreak (e) and decline (f) phase.
Figure 5. RDA analysis of environmental factors with bacterial community during the green tide (a), the outbreak (b) and decline (c) phase, and with archaeal community during the green tide (d), the outbreak (e) and decline (f) phase.
Microorganisms 10 01204 g005
Microorganisms 10 01204 g006 550
Figure 6. Correlations of Chao1 index, Shannon index, relative abundance of dominant taxa with chla content during the green tide (a), the outbreak (b) and decline (c) phase.
Figure 6. Correlations of Chao1 index, Shannon index, relative abundance of dominant taxa with chla content during the green tide (a), the outbreak (b) and decline (c) phase.
Microorganisms 10 01204 g006
Microorganisms 10 01204 g007 550
Figure 7. Co-occurrence network colored by taxonomy during the green tide (a), the outbreak (b) and decline phase (c), and colored by modularity during the green tide (d), the outbreak (e) and decline phase (f).
Figure 7. Co-occurrence network colored by taxonomy during the green tide (a), the outbreak (b) and decline phase (c), and colored by modularity during the green tide (d), the outbreak (e) and decline phase (f).
Microorganisms 10 01204 g007
Table
Table 1. Sampling information.
Table 1. Sampling information.
Sampling DateSample Name
XGZQMD
13 June 2019XG0613ZQ0613MD0613
27 June 2019XG0627ZQ0627MD0627
11 July 2019XG0711ZQ0711MD0711
18 July 2019XG0718ZQ0718MD0718
22 July 2019XG0722ZQ0722MD0722
8 August 2019XG0806ZQ0806MD0806
16 August 2019XG0816ZQ0816/
23 August 2019XG0823ZQ0823MD0823
29 August 2019XG0829ZQ0829MD0829
5 September 2019XG0905ZQ0905MD0905
Table
Table 2. Environmental factors in coastal Qingdao areas during the green tide.
Table 2. Environmental factors in coastal Qingdao areas during the green tide.
Temperature (°C)SalinitypHDO (mg/L)chla (μg/L)NH4+ (μmol/L)NO2− (μmol/L)NO3− (μmol/L)PO43− (μmol/L)
XG061318.7330.158.057.611.674.270.131.532.88
XG062720.9029.968.057.561.0724.350.141.674.31
XG071122.8029.837.987.311.185.370.081.150.30
XG071823.6029.788.006.74/23.470.172.300.26
XG072224.7729.497.997.342.9613.370.335.020.44
XG080625.8329.957.587.620.7512.660.262.900.61
XG081626.1029.287.757.480.8318.420.659.570.30
XG082326.1029.668.0210.711.5115.320.273.700.21
XG082926.0729.798.376.47/8.010.233.080.10
XG090525.3729.618.099.420.6528.700.526.843.92
ZQ061320.1030.088.119.392.9228.660.131.760.24
ZQ062722.0030.078.057.243.9118.780.111.370.89
ZQ071123.8030.018.045.593.375.690.070.73/
ZQ071823.7029.868.015.08/21.130.111.38/
ZQ072225.9029.868.047.944.865.710.081.06/
ZQ080626.0030.397.827.680.7811.170.202.092.68
ZQ081626.2030.008.029.101.5110.140.213.980.44
ZQ082326.1030.048.1110.201.1320.740.405.03/
ZQ082926.5030.178.2512.17/16.010.161.870.27
ZQ090525.2029.888.1012.860.9321.070.273.561.83
MD061318.6030.228.118.74/3.250.040.590.44
MD062720.8030.048.139.07/7.690.081.070.27
MD071122.9029.988.138.37/18.400.050.730.25
MD071823.0029.878.057.84/4.910.070.970.35
MD072224.3029.778.058.19/23.020.111.930.18
MD080625.8029.888.018.05/21.400.202.660.32
MD082326.2029.907.993.92/19.950.080.990.27
MD082926.4029.918.2111.05/22.230.314.760.22
MD090525.5029.748.088.35/23.030.385.454.67
/, Missing data.
Table
Table 3. Total OTUs, Chao1 index, Shannon index and Good’s coverage of bacterial and archaeal communities during the green tide.
Table 3. Total OTUs, Chao1 index, Shannon index and Good’s coverage of bacterial and archaeal communities during the green tide.
Bacterial CommunityArchaeal Community
OTUsChao1ShannonCoverageOTUsChao1ShannonCoverage
XG0613643750.364.440.99655765.271.640.9995
XG0627588705.503.800.99635360.501.650.9996
XG0711659824.283.430.99588283.002.470.9999
XG0718768871.214.610.9963101123.002.430.9996
XG0722735851.324.550.996289102.001.850.9995
XG0806671825.014.280.99565772.001.850.9994
XG0816750856.654.620.99615985.251.760.9994
XG0823756854.704.310.9963106115.232.550.9994
XG0829787890.154.870.996489114.002.450.9991
XG0905674765.884.690.99675362.002.330.9996
ZQ0613642738.043.840.66976465.201.920.9999
ZQ0627719784.324.630.99737379.071.610.9995
ZQ0711625748.044.040.99595969.912.110.9994
ZQ0718813932.724.760.99628392.072.130.9994
ZQ0722811903.244.790.996392117.501.960.9993
ZQ0806803935.374.680.995990106.872.190.9991
ZQ0816712789.114.550.9969102122.002.110.9994
ZQ0823650797.663.990.99566266.672.490.9997
ZQ0829837949.264.850.9962112120.572.600.9994
ZQ0905762915.254.410.99587693.102.390.9993
MD0613632753.153.890.99676767.751.950.9999
MD0627687781.614.450.99694956.001.610.9997
MD0711646769.103.340.99607075.602.420.9997
MD0718806916.574.520.996383113.002.180.9991
MD0722750830.534.230.99658284.502.370.9998
MD0806792945.384.370.99548698.672.190.9993
MD0823395539.733.430.99685764.502.950.9998
MD0829438631.802.310.99607580.602.360.9997
MD0905617812.633.830.99518387.582.390.9996
Table
Table 4. Spearman correlation analysis of Chao1 index, Shannon index, relative abundance of dominant taxa with chla content during the green tide, the outbreak and decline phase.
Table 4. Spearman correlation analysis of Chao1 index, Shannon index, relative abundance of dominant taxa with chla content during the green tide, the outbreak and decline phase.
PhasesBacterial CommunityArchaeal Community
Chao1ShannonRhodobacteralesFlavobacterialesChao1ShannonMGII
chla contentgreen tide−0.2120.0180.3260.638 **0.209−0.250−0.165
outbreak phase0.5480.857 **−0.0480.762 *0.6190.000−0.476
decline phase0.048−0.4760.6900.3330.6670.4760.143
* p < 0.05, ** p < 0.01.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhao, G.; He, H.; Wang, H.; Liang, Y.; Guo, C.; Shao, H.; Jiang, Y.; Wang, M. Variations in Marine Bacterial and Archaeal Communities during an Ulva prolifera Green Tide in Coastal Qingdao Areas. Microorganisms 2022, 10, 1204. https://doi.org/10.3390/microorganisms10061204

AMA Style

Zhao G, He H, Wang H, Liang Y, Guo C, Shao H, Jiang Y, Wang M. Variations in Marine Bacterial and Archaeal Communities during an Ulva prolifera Green Tide in Coastal Qingdao Areas. Microorganisms. 2022; 10(6):1204. https://doi.org/10.3390/microorganisms10061204

Chicago/Turabian Style

Zhao, Guihua, Hui He, Hualong Wang, Yantao Liang, Cui Guo, Hongbing Shao, Yong Jiang, and Min Wang. 2022. "Variations in Marine Bacterial and Archaeal Communities during an Ulva prolifera Green Tide in Coastal Qingdao Areas" Microorganisms 10, no. 6: 1204. https://doi.org/10.3390/microorganisms10061204

APA Style

Zhao, G., He, H., Wang, H., Liang, Y., Guo, C., Shao, H., Jiang, Y., & Wang, M. (2022). Variations in Marine Bacterial and Archaeal Communities during an Ulva prolifera Green Tide in Coastal Qingdao Areas. Microorganisms, 10(6), 1204. https://doi.org/10.3390/microorganisms10061204

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop