Phylogenetic Responses of Marine Free-Living Bacterial Community to Phaeocystis globosa Bloom in Beibu Gulf, China

The influence of temperature and river runoff on phytoplankton community diversity in Beibu Gulf: insight from 18 S rDNA metabarcoding analysis

BACKGROUND

Sanniang Bay (SNB) and Dafeng River (DFR), located in the northern Beibu Gulf, is well-known as one of the eight habitats for humpback dolphins in China. This region is representative of typical estuarine and bay ecosystems and produce complex hydrodynamic seawater conditions. Moreover, anthropogenic pressure, such as eutrophication and large-scale infrastructure projects, have caused ongoing habitat deterioration and loss. It is urgent to know the phytoplankton community and their relationships with environmental factors in this region.

RESULTS

In this study, we assessed the diversity and assembly mechanisms of phytoplankton communities, as well as their relationship with the physicochemical characteristics of seawater in SNB and DFR region using 18 S rDNA metabarcoding analysis. The results showed that seasonal changes markedly impacted the alpha diversity of the phytoplankton community. From March to July, with the average temperature increasing from 25.2℃ to 28.1℃,the Shannon or Species Richness were negatively correlated with temperature. During hot season (in Sep, average temperature 32.1℃), phytoplankton diversity was negatively correlated with nutrients (NH4 +, NO3-, PO43-, TN). Additionally, during the rainy season, the Bray-Curtis similarity of the phytoplankton community was significantly lower than during the dry season. In March, the distance among the sampling sites was most strongly and positively correlated with the Bray-Curtis dissimilarity. Stochastic processes, specifically dispersal limitation and ecological drift, are the primary drivers of community assembly, while deterministic assembly processes (mainly heterogeneous selection) contribute a relatively minor portion (< 17%).

CONCLUSIONS

Rising temperature diminished the diversity of phytoplankton in SNB and DFR, and nutrient inputs and eutrophication in estuarine areas will aggravate the loss of phytoplankton diversity.

The micro-ecological feature of colonies is a potential strategy for Phaeocystis globosa bloom formation

Phaeocystis globosa is among the dominant microalgae associated with harmful algal blooms. P. globosa has a polymorphic life cycle and its ecological success has been attributed to algal colony formation, however, few studies have assessed differences in microbial communities and their functional profiles between intra- and extra-colonies during P. globosa blooms. To address this, environmental and metagenomics tools were used to conduct a time-series analysis of the bacterial composition and metabolic characteristics of intra- and extra-colonies during a natural P. globosa bloom. The results show that bacterial composition, biodiversity, and network interactions differed significantly between intra- and extra-colonies. Dominant extra-colonial bacteria were Bacteroidia and Saccharimonadis, while dominant intra-colonial bacteria included Alphaproteobacteria and Gammaproteobacteria. Despite the lower richness and diversity observed in the intra-colonial bacterial community, relative to extra-colonies, the complexity and interconnectedness of the intra-colonial networks were higher. Regarding bacterial function, more functional genes were enriched in substance metabolism (polysaccharides, iron element and dimethylsulfoniopropionate) and signal communication (quorum sensing, indoleacetic acid-IAA) pathways in intra- than in extra-colonies. Conceptual model construction showed that microbial cooperative synthesis of ammonium, vitamin B12, IAA, and siderophores were strongly related to the P. globosa bloom, particularly in the intra-colonial environment. Overall, our data highlight the differences in bacterial structure and functions within and outside the colony during P. globosa blooms. These findings represent fundamental information indicating that phenotypic heterogeneity is a selective strategy that improves microbial population competitiveness and environmental adaptation, benefiting P. globosa bloom formation and persistence.

Seasonal and anthropogenic influences on bacterioplankton communities: ecological impacts in the coastal waters of Qinhuangdao, Northern China

Marine bacterioplankton play a crucial role in the cycling of carbon, nitrogen, and phosphorus in coastal waters. And the impact of environmental factors on bacterial community structure and ecological functions is a dynamic ongoing process. To systematically assess the relationship between environmental changes and bacterioplankton communities, this study delved into the spatiotemporal distribution and predicted metabolic characteristics of bacterioplankton communities at two estuarine beaches in Northern China. Coastal water samples were collected regularly in spring, summer, and autumn, and were analyzed in combination with environmental parameters and bacterioplankton community. Results indicated significant seasonal variations in bacterioplankton communities as Bacteroidetes and Actinobacteria were enriched in spring, Cyanobacteria proliferated in summer. While Pseudomonadota and microorganisms associated with organic matter decomposition prevailed in autumn, closely linked to seasonal variation of temperature, light and nutrients such as nitrogen and phosphorus. Particularly in summer, increased tourism activities and riverine inputs significantly raised nutrient levels, promoting the proliferation of specific photosynthetic microorganisms, potentially linked to the occurrence of phytoplankton blooms. Spearman correlation analysis further revealed significant correlations between bacterioplankton communities and environmental factors such as salinity, chlorophyll a, and total dissolved phosphorus (TDP). Additionally, the metabolic features of the spring bacterioplankton community were primarily characterized by enhanced activities in the prokaryotic carbon fixation pathways, reflecting rapid adaptation to increased light and temperature, as well as significant contributions to primary productivity. In summer, the bacterial communities were involved in enhanced glycolysis and biosynthetic pathways, reflecting high energy metabolism and responses to increased light and biomass. In autumn, microorganisms adapted to the accelerated decomposition of organic matter and the seasonal changes in environmental conditions through enhanced amino acid metabolism and material cycling pathways. These findings demonstrate that seasonal changes and human activities significantly influence the structure and function of bacterioplankton communities by altering nutrient dynamics and physical environmental conditions. This study provides important scientific insights into the marine biological responses under global change.

Rare subcommunity maintains the stability of ecosystem multifunctionality by deterministic assembly processes in subtropical estuaries

Microorganisms, especially rare microbial species, are crucial in estuarine ecosystems for driving biogeochemical processes and preserving biodiversity. However, the understanding of the links between ecosystem multifunctionality (EMF) and the diversity of rare bacterial taxa in estuary ecosystems remains limited. Employing high-throughput sequencing and a variety of statistical methods, we assessed the diversities and assembly process of abundant and rare bacterioplankton and their contributions to EMF in a subtropical estuary. Taxonomic analysis revealed Proteobacteria as the predominant phylum among both abundant and rare bacterial taxa. Notably, rare taxa demonstrated significantly higher taxonomic diversity and a larger species pool than abundant taxa. Additionally, our findings highlighted that deterministic assembly processes predominantly shape microbial communities, with heterogeneous selection exerting a stronger influence on rare taxa. Further analysis reveals that rare bacterial beta-diversity significantly impacts to EMF, whereas alpha diversity did not. The partial least squares path modeling (PLS-PM) analysis demonstrated that the beta diversity of abundant and rare taxa, as the main biotic factor, directly affected EMF, while temperature and total organic carbon (TOC) were additional key factors to determine the relationship between beta diversity and EMF. These findings advance our understanding of the distribution features and ecological knowledge of the abundant and rare taxa in EMF in subtropical estuaries, and provide a reference for exploring the multifunctionality of different biospheres in aquatic environments.

Responses of microeukaryotic community structure to a Phaeocystis globosa bloom in a semi-enclosed subtropical bay

The occurrence of Phaeocystis globosa, a harmful algal bloom species in Chinese coastal waters, has significant impacts on marine organisms and poses a threat to the safety of coastal nuclear power plants. Although previous studies have established a close association between P. globosa blooms and the bacterial community, the relationship between the microeukaryotic community and P. globosa blooms remains poorly understood. In this study, the variations in the microeukaryotic community resulting from a P. globosa bloom were analyzed using 18S rRNA gene amplicon sequencing. The results indicated that the diversity of the microeukaryotic community during the bloom phase was significantly higher than that during the dissipation phase. The microeukaryotic community compositions varied significantly between the two phases of the P. globosa bloom. During the bloom phase, the dominant microeukaryotic was Viridiplantae, which was then replaced by Dinoflagellata during the dissipation phase. Co-occurrence network analysis showed that the relationship among the microeukaryotic community during the bloom phase was more complex than that during the dissipation phase, and the keystone taxa varied as the bloom progressed. Additionally, microeukaryotic community assembly was primarily driven by stochastic processes during the bloom phase based on the β-nearest taxon distance, whereas it was driven by both deterministic processes and stochastic processes during the dissipation phase. Overall, our findings provide novel insight into the mechanisms and interactions involved in microeukaryotic community dynamics in environments disturbed by P. globosa blooms.

Microbial community composition and metabolic potential during a succession of algal blooms from Skeletonema sp. to Phaeocystis sp

Elucidating the interactions between algal and microbial communities is essential for understanding the dynamic mechanisms regulating algal blooms in the marine environment. Shifts in bacterial communities when a single species dominates algal blooms have been extensively investigated. However, bacterioplankton community dynamics during bloom succession when one algal species shift to another is still poorly understood. In this study, we used metagenomic analysis to investigate the bacterial community composition and function during algal bloom succession from Skeletonema sp. to Phaeocystis sp. The results revealed that bacterial community structure and function shifted with bloom succession. The dominant group in the Skeletonema bloom was Alphaproteobacteria, while Bacteroidia and Gammaproteobacteria dominated the Phaeocystis bloom. The most noticeable feature during the successions was the change from Rhodobacteraceae to Flavobacteriaceae in the bacterial communities. The Shannon diversity indices were significantly higher in the transitional phase of the two blooms. Metabolic reconstruction of the metagenome-assembled genomes (MAGs) showed that dominant bacteria exhibited some environmental adaptability in both blooms, capable of metabolizing the main organic compounds, and possibly providing inorganic sulfur to the host algae. Moreover, we identified specific metabolic capabilities of cofactor biosynthesis (e.g., B vitamins) in MAGs in the two algal blooms. In the Skeletonema bloom, Rhodobacteraceae family members might participate in synthesizing vitamin B₁ and B₁₂ to the host, whereas in the Phaeocystis bloom, Flavobacteriaceae was the potential contributor for synthesizing vitamin B₇ to the host. In addition, signal communication (quorum sensing and indole-3-acetic acid molecules) might have also participated in the bacterial response to bloom succession. Bloom-associated microorganisms showed a noticeable response in composition and function to algal succession. The changes in bacterial community structure and function might be an internal driving factor for the bloom succession.

Bacteria Associated With Phaeocystis globosa and Their Influence on Colony Formation

Phaeocystis globosa (P. globosa) is one of the dominant algae during harmful algal blooms (HABs) in coastal regions of Southern China. P. globosa exhibits complex heteromorphic life cycles that could switch between solitary cells and colonies. The ecological success of P. globosa has been attributed to its colony formation, although underlying mechanisms remain unknown. Here, we investigated different bacterial communities associated with P. globosa colonies and their influence on colony formation of two P. globosa strains isolated from coastal waters of Guangxi (GX) and Shantou (ST). Eight operational taxonomic units (OTUs) were observed in ST co-cultures and were identified as biomarkers based on Linear discriminant analysis Effect Size (LEfSe) analysis, while seven biomarkers were identified in P. globosa GX co-cultures. Bacterial communities associated with the P. globosa GX were more diverse than those of the ST strain. The most dominant phylum in the two co-cultures was Proteobacteria, within which Marinobacter was the most abundant genus in both GX and ST co-cultures. Bacteroidota were only observed in the GX co-cultures and Planctomycetota were only observed in the ST co-cultures. Co-culture experiments revealed that P. globosa colony formation was not influenced by low and medium cell densities of Marinobacter sp. GS7, but was inhibited by high cell densities of Marinobacter sp. GS7. Overall, these results indicated that the associated bacteria are selected by different P. globosa strains, which may affect the colony formation and development of P. globosa.

Trophic states regulate assembly processes and network structures of small chromophytic phytoplankton communities in estuarine and coastal ecosystem

Small chromophytic phytoplankton (SCP) are anticipated to be more important for a significant proportion of primary production in estuarine-coastal ecosystems. However, responses of SCP community to coastal eutrophication are still unclear. In this study, we investigated diversity, co-occurrence and assembly features of SCP communities, as well as relationship with environmental factors in subtropical Beibu Gulf. The results exhibited that the alpha diversity and beta diversity of SCP communities were significantly different among eutrophic states. Co-occurrence network revealed a complex interaction that most amplicon sequence variants (ASVs) in modules of the network were specific to trophic states. Further, phylogenetic based β-nearest taxon distance analyses revealed that stochastic processes mainly provided 69.26% contribution to SCP community assembly, whereas deterministic processes dominated community assembly in heavy eutrophic state. Overall, our findings elucidate the mechanism of diversity and assembly in SCP community and promote the understanding of SCP ecology related to subtropical coastal eutrophication.

Co-occurrence of chromophytic phytoplankton and the Vibrio community during Phaeocystis globosa blooms in the Beibu Gulf

Accumulating research evidence has revealed that harmful algal blooms (HABs) can substantially affect the community structures of phytoplankton and heterotrophic bacteria in marine ecosystems. However, little is known about their species-specific interactions between phytoplankton and heterotrophic bacteria during the HABs period and about their interaction shifts in response to blooms. From this perspective, we investigated the co-occurrence of chromophytic phytoplankton and Vibrio during Phaeocystis globosa blooms in the Beibu Gulf. The results showed that Vibrio communities were distinct during the blooms, and P. globosa blooms resulted in a decline in phytoplankton alpha diversity, revealing that the blooms could affect their community compositions. The regression lines between the Shannon indices and Bray-Curtis distances of phytoplankton and Vibrio showed positive correlations with each other (p < 0.001), suggesting that they may have intrageneric symbiotic interactions overall. In addition, network analysis further demonstrated that relationships between phytoplankton and Vibrio were dominated by positive correlations, and more interaction modules were observed during the blooms, revealing that the blooms intensified synergistic association and mutual symbiotic interactions between them. Environmental factors (SiO₃^2-, NH₄⁺, NO₃^- and TN,) and P. globosa density more deeply affected network interactions between phytoplankton and Vibrio during the periods of P. globosa blooms than those before the blooms and after the blooms. This study provided new insight to elucidate community structure and interaction relationships between phytoplankton and Vibrio in response to P. globosa blooms and their ecological effects in marine ecosystems.

The Eco-Evo Mandala: Simplifying Bacterioplankton Complexity into Ecohealth Signatures

The microbiome emits informative signals of biological organization and environmental pressure that aid ecosystem monitoring and prediction. Are the many signals reducible to a habitat-specific portfolio that characterizes ecosystem health? Does an optimally structured microbiome imply a resilient microbiome? To answer these questions, we applied our novel Eco-Evo Mandala to bacterioplankton data from four habitats within the Great Barrier Reef, to explore how patterns in community structure, function and genetics signal habitat-specific organization and departures from theoretical optimality. The Mandala revealed communities departing from optimality in habitat-specific ways, mostly along structural and functional traits related to bacterioplankton abundance and interaction distributions (reflected by ϵ and λ as power law and exponential distribution parameters), which are not linearly associated with each other. River and reef communities were similar in their relatively low abundance and interaction disorganization (low ϵ and λ) due to their protective structured habitats. On the contrary, lagoon and estuarine inshore reefs appeared the most disorganized due to the ocean temperature and biogeochemical stress. Phylogenetic distances (D) were minimally informative in characterizing bacterioplankton organization. However, dominant populations, such as Proteobacteria, Bacteroidetes, and Cyanobacteria, were largely responsible for community patterns, being generalists with a large functional gene repertoire (high D) that increases resilience. The relative balance of these populations was found to be habitat-specific and likely related to systemic environmental stress. The position on the Mandala along the three fundamental traits, as well as fluctuations in this ecological state, conveys information about the microbiome's health (and likely ecosystem health considering bacteria-based multitrophic dependencies) as divergence from the expected relative optimality. The Eco-Evo Mandala emphasizes how habitat and the microbiome's interaction network topology are first- and second-order factors for ecosystem health evaluation over taxonomic species richness. Unhealthy microbiome communities and unbalanced microbes are identified not by macroecological indicators but by mapping their impact on the collective proportion and distribution of interactions, which regulates the microbiome's ecosystem function.

Spatiotemporal dynamics of marine microbial communities following a Phaeocystis bloom: biogeography and co-occurrence patterns

Marine microbes play important roles in the development of phytoplankton blooms. The diversity and composition of free living (FL) and particle attached (PA) microbial communities have been well studied, while little is known about their geographic and co-occurrence patterns, especially during the subsiding process of Phaeocystis globosa blooms. Herein, the beta-diversity of FL and PA microbial communities in both the surface and bottom layers of different habitats were comprehensively examined during succession of a P. globosa bloom event. The results showed that microbial communities from bloom and non-bloom sites exhibited distinct community compositions. Among the different sampling sites, the community similarities decreased with spatial distance, in which the FL communities' similarity in bottom waters was more influenced by spatial variation. The variation of microbial communities was mostly attributed to environmental selection, spatial distance, and the abundance of P. globosa successively. The co-occurrence networks of microbial communities in bloom and non-bloom waters differed in terms of structure and composition, and the bloom network had more links and closer relationships between genera than the non-bloom network. The correlation among genera and modules suggested that the bloom microbes were likely driven by high environmental selection and low competitive effect between each other.