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

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.


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.

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.

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%.

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.

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. NH 4 + concentrations varied from 3.25-28.70 µmol/L, while NO 2 − and NO 3 − concentrations were relatively low with values of 0.04-0.65 µmol/L and 0.59-9.57 µmol/L, respectively. PO 4 3− 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), NO 2 − concentration (p < 0.01) and NO 3 − concentration (p < 0.01) varied significantly between the outbreak and decline phases.

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).
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 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

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)

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 NO 2 − 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  (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 Flavobacte-riales (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). 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 PO4 3− 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).

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 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.

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 (NH 4 + , NO 2 − , NO 3 − and PO 4 3− ) 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.

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 cooccurrence 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.

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.