Co-occurrence patterns and assembly processes of microeukaryotic communities in an early-spring diatom bloom

Interaction patterns and assembly mechanisms of dinoflagellates and diatoms in a coastal bay suffering from long-term eutrophication

Dinoflagellates and diatoms are highly prevalent and ecologically important phytoplankton in coastal waters, greatly contributing to primary productivity in marine ecosystems. Although their composition and diversity have been extensively elucidated in the open ocean, their interaction patterns and community assembly in long-term eutrophic coastal waters remain poorly understood. This investigation aimed to elucidate the seasonal successional patterns of dinoflagellates and diatoms by 18S rRNA gene amplicon sequencing in a semi-enclosed bay. The results revealed that dinoflagellate and diatom communities have pronounced seasonal succession patterns, which are primarily associated with temperature. Furthermore, the most prevalent species throughout the year were Heterocapsa rotundata and Skeletonema costatum. Moreover, the assembly of dinoflagellate and diatom communities was mainly dominated by stochastic processes, with drift being the major factor. The co-occurrence of dinoflagellates and diatoms showed seasonal patterns, with the highest interactions observed in autumn. In addition, interactions of Syndiniales with dinoflagellates and diatoms highlighted the roles of parasites in eutrophic conditions. Flavobacteriaceae and Rhodobacteraceae are the bacterial taxa that most frequently interacted with dinoflagellates and diatoms, with interactions between dinoflagellates and bacteria being more complex than those between diatoms and bacteria. Overall, this study provides results that deepen our understanding of the phytoplankton dynamics in coastal eutrophic waters.IMPORTANCEDinoflagellates and diatoms are major phytoplankton groups in coastal waters. The composition and diversity of dinoflagellates and diatoms in the open ocean have been well documented; however, it remains uncertain to what extent their adaptation to long-term eutrophic conditions influences their response to environmental disturbances. Here, we investigated the interactions and assembly processes of dinoflagellates and diatoms in a eutrophic bay throughout the whole year. Our findings revealed that interactions between dinoflagellates and diatoms are primarily shaped by seasonal transitions, while prolonged eutrophic conditions tend to amplify stochastic processes in community assembly. These findings provide novel perspectives on the influence of long-term eutrophication on phytoplankton dynamics within eutrophic waters.

Time-series metagenomics reveals changing protistan ecology of a temperate dimictic lake

BACKGROUND

Protists, single-celled eukaryotic organisms, are critical to food web ecology, contributing to primary productivity and connecting small bacteria and archaea to higher trophic levels. Lake Mendota is a large, eutrophic natural lake that is a Long-Term Ecological Research site and among the world's best-studied freshwater systems. Metagenomic samples have been collected and shotgun sequenced from Lake Mendota for the last 20 years. Here, we analyze this comprehensive time series to infer changes to the structure and function of the protistan community and to hypothesize about their interactions with bacteria.

RESULTS

Based on small subunit rRNA genes extracted from the metagenomes and metagenome-assembled genomes of microeukaryotes, we identify shifts in the eukaryotic phytoplankton community over time, which we predict to be a consequence of reduced zooplankton grazing pressures after the invasion of a invasive predator (the spiny water flea) to the lake. The metagenomic data also reveal the presence of the spiny water flea and the zebra mussel, a second invasive species to Lake Mendota, prior to their visual identification during routine monitoring. Furthermore, we use species co-occurrence and co-abundance analysis to connect the protistan community with bacterial taxa. Correlation analysis suggests that protists and bacteria may interact or respond similarly to environmental conditions. Cryptophytes declined in the second decade of the timeseries, while many alveolate groups (e.g., ciliates and dinoflagellates) and diatoms increased in abundance, changes that have implications for food web efficiency in Lake Mendota.

CONCLUSIONS

We demonstrate that metagenomic sequence-based community analysis can complement existing efforts to monitor protists in Lake Mendota based on microscopy-based count surveys. We observed patterns of seasonal abundance in microeukaryotes in Lake Mendota that corroborated expectations from other systems, including high abundance of cryptophytes in winter and diatoms in fall and spring, but with much higher resolution than previous surveys. Our study identified long-term changes in the abundance of eukaryotic microbes and provided context for the known establishment of an invasive species that catalyzes a trophic cascade involving protists. Our findings are important for decoding potential long-term consequences of human interventions, including invasive species introduction. Video Abstract.

Diversity and distribution analysis of eukaryotic communities in the Xiangshan Bay, East China sea by metabarcoding approach

Eukaryotic communities play an important role in the coastal ecosystem of Xiangshan Bay, a narrow semi-closed bay famous for fisheries and marine farming. However, information on the diversity and composition of eukaryotic communities in Xiangshan Bay remains unclear. In this study, the metabarcoding approach was utilized to comprehensively investigate the eukaryotic plankton community structure and dominant taxa, particularly eukaryotic microalgae, in the Xiangshan Bay over a period of four months in 2018. The results showed that the three major phyla were Arthropoda, Chlorophyta, and Bacillariophyta. The richness indices revealed that species richness peaked in February and was at its lowest in May. Diversity indices showed that the samples collected in May had the lowest diversity. Centropages was detected in the samples of all months, however, its highest dominance was observed in the samples collected in February. In addition, compared to other months, a greater proportion of eukaryotic microalgae was witnessed in March. The three eukaryotic algae with highest abundances in March were Cyclotella, Prorocentrum, and Thalassiosira. Moreover, high diversity of pico-sized (0.2-2.0 μm) phytoplankton (which are often easily missed by microscopy) was discovered in this study by using metabarcoding approach. This study highlights the strength and significance of the metabarcoding approach to uncover a large number of eukaryotic species which remains undetectable during application of conventional approaches. The findings of this study reveals that the eukaryotic community structure varies noticeably in both time and space throughout sampling period, with temperature being the most important environmental factor influencing these changes. This study lays a solid foundation to understand eukaryotic plankton composition, temporal and spatial dynamics and the distribution mechanism of eukaryotic plankton community in Xiangshan Bay, providing theoretical reference for further studies related to marine ecology.

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.

Spatiotemporal assembly and functional composition of planktonic microeukaryotic communities along productivity gradients in a subtropical lake

Microeukaryotes play crucial roles in the microbial loop of freshwater ecosystems, functioning both as primary producers and bacterivorous consumers. However, understanding the assembly of microeukaryotic communities and their functional composition in freshwater lake ecosystems across diverse environmental gradients remains limited. Here, we utilized amplicon sequencing of 18S rRNA gene and multivariate statistical analyses to examine the spatiotemporal and biogeographical patterns of microeukaryotes in water columns (at depths of 0.5, 5, and 10 m) within a subtropical lake in eastern China, covering a 40 km distance during spring and autumn of 2022. Our results revealed that complex and diverse microeukaryotic communities were dominated by Chlorophyta (mainly Chlorophyceae), Fungi, Alveolata, Stramenopiles, and Cryptophyta lineages. Species richness was higher in autumn than in spring, forming significant hump-shaped relationships with chlorophyll a concentration (Chl-a, an indicator of phytoplankton biomass). Microeukaryotic communities exhibited significant seasonality and distance-decay patterns. By contrast, the effect of vertical depth was negligible. Stochastic processes mainly influenced the assembly of microeukaryotic communities, explaining 63, 67, and 55% of community variation for spring, autumn, and both seasons combined, respectively. Trait-based functional analysis revealed the prevalence of heterotrophic and phototrophic microeukaryotic plankton with a trade-off along N:P ratio, Chl-a, and dissolved oxygen (DO) gradients. Similarly, the mixotrophic proportions were significantly and positively correlated with Chl-a and DO concentrations. Overall, our findings may provide useful insights into the assembly patterns of microeukaryotes in lake ecosystem and how their functions respond to environmental changes.

Bacterial and Microeukaryotic Community Compositions and Their Assembly Processes in Lakes on the Eastern Qinghai-Tibet Plateau

Bacterial and microeukaryotic community compositions and their assembly processes have remained challenging and remained unclear in lake ecosystems on the Qinghai-Tibet Plateau (QTP). We revealed the diversity and community compositions, driving factors, ecological assembly processes, and co-occurrence networks of bacterial and microeukaryotic communities in water bodies of the eight lake ecosystems across the Eastern QTP. The results demonstrated that the predominant bacteria in most samples were Proteobacteria, with an average relative abundance of 41.78%, whereas the most abundant of microeukaryotes differed among the sample sites. The redundancy analysis revealed that latitude and pH were the most important driving factors in shaping the bacterial and microeukaryotic community compositions. Homogeneous selection (56.40%) was the dominant process in assembling the bacterial communities, whereas dispersal limitation (67.24%) was the major process in governing the microeukaryotic communities. Furthermore, dissolved organic carbon and salinity were the major factors mediating the balance of deterministic and stochastic assembly processes in the bacterial and microeukaryotic communities. Both the bacterial and microeukaryotic community co-occurrence networks exhibited topological features of modularity and non-random topological features. The results offer insights into the mechanisms underpinning bacterial and microeukaryotic diversities and communities in the lake ecosystems on the QTP.

Hydrological disturbances enhance stochastic assembly processes and decrease network stability of algae communities in a highland floodplain system

Floodplains are hotspots for biodiversity research and conservation worldwide. Hydrological disturbances can profoundly influence the ecological processes and functions of floodplain systems by altering key biological groups such as algae communities. However, the impacts of flood disturbance on the assembly processes and co-occurrence patterns of algae communities in floodplain ecosystems are still unclear. To ascertain the response patterns of algae communities to flood disturbance, we characterized planktonic and benthic algae communities in 144 water and sediment samples collected from the Tibetan floodplain during non-flood and flood periods based on 23S ribosomal RNA gene sequencing. Results showed that planktonic algae exhibited higher diversity and greater compositional variations compared with benthic communities after flood disturbance. Flooding promoted algae community homogenization at horizontal (rivers vs. oxbow lakes) and vertical levels (water vs. sediment). Stochastic processes governed the assembly of distinct algae communities, and their ecological impacts were enhanced in response to flooding. In the non-flood period, dispersal limitation (81.78 %) was the primary ecological process driving algae community assembly. In the flood period, the relative contribution of ecological drift (72.91 %) to algae community assembly markedly increased, with dispersal limitation (22.61 %) being less important. Flooding reduced the interactions among algae taxa, resulting in lower network complexity and stability. Compared with the planktonic algae subnetworks, the benthic subnetworks showed greater stability in the face of flooding. Findings of this study broaden our understanding of how algae communities respond to hydrological disturbances from an ecological perspective and could be useful for the management of highland floodplain ecosystems.

Temporal patterns of algae in different urban lakes and their correlations with environmental variables in Xi'an, China

Urban lakes were critical in aquatic ecology environments, but how environmental factors affected the distribution and change characteristics of algal communities in urban lakes of Xi'an city was not clearly. Here, we investigated the algal community structure of six urban lakes in Xi'an and evaluated the effects of water quality parameters on algae. The results indicated that the significant differences on physicochemical parameters existed in different urban lakes. The maximum concentration of total phosphorus in urban lakes was (0.18 ± 0.01) mg/L and there was a phenomenon of phosphorus limitation. In addition, 51 genera of algae were identified and Chlorella sp. was the dominant algal species, which was affiliated with Chlorophyta. Network analysis elucidated that each lake had a unique algal community network and the positive correlation was dominant in the interaction between algae species, illustrating that mature microbial communities existed or occupied similar niches. Redundancy analysis illustrated that environmental factors explained 47.35% variance of algal species-water quality correlation collectively, indicating that water quality conditions had a significant influence on the temporal variations of algae. Structural equation model further verified that algal community structure was directly or indirectly regulated by different water quality conditions. Our study shows that temporal patterns of algal communities can reveal the dynamics and interactions of different urban ecosystem types, providing a theoretical basis for assessing eutrophication levels and for water quality management.

Contrasting Mechanisms Determine the Microeukaryotic and Syndiniales Community Assembly in a Eutrophic bay

Syndiniales is a diverse parasitic group, increasingly gaining attention owing to its high taxonomic diversity in marine ecosystems and inhibitory effects on the dinoflagellate blooms. However, their seasonal dynamics, host interactions, and mechanisms of community assembly are largely unknown, particularly in eutrophic waters. Here, using 18S rRNA gene amplicon sequencing, we intended to elucidate the interactions between Syndiniales and microeukaryotes, as well as community assembly processes in a eutrophic bay. The results showed that Syndiniales group II was dominating throughout the year, with substantially higher abundance in the winter and spring, whereas Syndiniales group I was more abundant in the summer and autumn. Temperature and Dinoflagellata were the most important abiotic and biotic factors driving variations of the Syndiniales community, respectively. The assembly processes of microeukaryotes and Syndiniales were completely different, with the former being controlled by a balance between homogeneous selection and drift and the latter being solely governed by drift. Network analysis revealed that Syndiniales group II had the largest number of interactions with microeukaryotes, and they primarily associated with Dinoflagellata in the winter, while interactions with Chlorophyta and Bacillariophyta increased dramatically in summer and autumn. These findings provide significant insights in understanding the interactions and assembly processes of Syndiniales throughout the year, which is critical in revealing the roles of single-celled parasites in driving protist dynamics in eutrophic waters.

Seasonal assessment of the distribution, source apportionment, and risk of water-contaminated polycyclic aromatic hydrocarbons (PAHs)

The research aims to evaluate the seasonal differences in the distribution, source, and risks of water-contaminated PAHs. The PAHs were extracted by the liquid-liquid method and analyzed with GC-MS, and a total of eight PAHs were detected. There was a percentage increase in the average concentration of the PAHs from the wet to the dry season in the range of 20 (Anthracene)-350 (Pyrene)%. Total PAHs (∑PAHs) range from 0.31 to 1.23 mg/l in the wet period and from 0.42 to 1.96 mg/l in the dry period. The distribution of the average PAHs in mg/l showed that Fluoranthene ≤ Pyrene < Acenaphthene < Fluorene < Phenanthrene < Acenaphthylene < Anthracene < Naphthalene in wet period and while Fluoranthene < Acenaphthene < Pyrene < Fluorene < Phenanthrene < Acenaphthylene < Anthracene < Naphthalene in the dry period. The children were exposed to non-carcinogenic risk through non-dietary ingestion due to the accumulative effect (HI) of the PAHs in the dry period. Furthermore, the naphthalene was responsible for ecological and carcinogenic risk in the wet period, while the fluorene, phenanthrene, and anthracene were responsible for ecological and carcinogenic risk in the dry period. However, while adults and children are both susceptible to carcinogenic risk through the oral channel during the dry period, only children are susceptible to non-carcinogenic risk through this pathway. The multivariate statistical analysis revealed the influence of physicochemical parameters on the detected PAHs and also showed the PAHs' sources to be mainly combustion, pyrolysis, and vehicular emission.

Phytoplankton and Microzooplankton Community Structure and Assembly Mechanisms in Northwestern Pacific Ocean Estuaries with Environmental Heterogeneity and Geographic Segregation

Phytoplankton and microzooplankton are crucial players in marine ecosystems and first responders to environmental changes, but their community structures and how they are shaped by environmental conditions have rarely been studied simultaneously. In this study, we conducted an eDNA metabarcoding sequencing combined with multiple statistical methods to simultaneously analyze the phytoplankton and microzooplankton in Liaohe (LH) and Yalujiang (YLJ) estuaries. The major objective was to examine how plankton community structure and assembly mechanism may differ between two estuaries with similar latitudinal position and climate but geographical segregation and differential level of urbanization (more in LH). Clear differences in diversity and composition of phytoplankton and microzooplankton communities between LH and YLJ estuaries were observed. Richness of phytoplankton was significantly higher in LH than YLJ, while richness of microzooplankton was higher in YLJ. The magnitude of intrahabitat variations in phytoplankton communities was significantly stronger than that of microzooplankton. Some phytoplankton and microzooplankton taxa also showed interhabitat differences in their relative abundances. Phytoplankton showed a stronger geographic distance-decay of similarity than microzooplankton, while significant environmental distance-decay of similarity in microzooplankton was found in the less urbanized YLJ estuary. Community assembly of phytoplankton was, based on the neutral community models, driven primarily by stochastic processes, while deterministic processes contributed more for microzooplankton. Furthermore, we detected wider habitat niche breadths and stronger dispersal abilities in phytoplankton than in microzooplankton. These results suggest that passive dispersal shapes the phytoplankton community whereas environmental selection shapes the microzooplankton community. IMPORTANCE Understanding the underlying mechanisms shaping a metacommunity is useful to management for improving the ecosystem function. The research presented in the manuscript mainly tried to address the effects of habitat geography and environmental conditions on the phytoplankton and microzooplankton communities, and the underlying mechanisms of community assembly in temperate estuaries. In order to achieve this purpose, we developed a metabarcoding sequencing method based on 18S rRNA gene. The phytoplankton and microzooplankton communities from two estuaries with similar latitude and climatic conditions but obvious geographical segregation and significant environmental heterogeneity were investigated. The results of our study could lay a solid foundation for ascertaining phytoplankton and microzooplankton communities in estuaries with obvious environmental heterogeneity and geographic segregation and mechanisms underlying community assembly.

Novel insights in seasonal dynamics and co-existence patterns of phytoplankton and micro-eukaryotes in drinking water reservoir, Northwest China: DNA data and ecological model

Although associations between phytoplankton and micro-eukaryotes have been studied in aquatic ecosystems, there are still knowledge gaps in comprehending their dynamics and interactions in drinking water reservoirs. Here, the seasonal dynamics of phytoplankton and micro-eukaryotic diversities and their co-existence patterns were studied in a drinking water reservoir, Northwest China. The highest phytoplankton diversity was observed in summer, and Chlorella sp. that belongs to Chlorophyta was the most abundant genus. The highest eukaryotic diversity was also detected in summer, and Rimostrombidium sp. that belongs to Ciliophora was the most dominant genus. Mantel test showed that the phytoplankton diversity was significantly correlated with ammonia nitrogen (r = 0.561, p = 0.001) and dissolved organic carbon (r = 0.267, p = 0.017), while the eukaryotic diversity was significantly associated with ammonia nitrogen (r = 0.265, p = 0.034) and temperature (r = 0.208, p = 0.046). PLS-PM (Partial Least Squares Path Modeling) further revealed that nutrients (P < 0.01) significantly affected the phytoplankton diversity, while nutrients (P < 0.01) and temperature (P < 0.01) significantly influenced the eukaryotic diversity. Co-occurrence network displayed the primarily positive interactions (77.66% positive and 22.34% negative) between phytoplankton and micro-eukaryotes. These findings could deepen our understanding of interactions between phytoplankton and micro-eukaryotes and their driving factors under changing aquatic environments of drinking water reservoirs.

Change of the structure and assembly of bacterial and photosynthetic communities by the ecological engineering practices in Dianchi Lake

Cyanobacterial bloom challenges the aquatic ecosystem and ecological restoration is an effective approach for cyanobacterial bloom control, but the change of aquatic community after ecological restoration is still unclear. Dianchi Lake is an eutrophic lake with frequent cyanobacterial blooms in China, and recent ecological restoration projects in Caohai (north part) have a satisfactory performance. In this study, we collected 249 water samples at 23 sites from Dianchi Lake to explore the relationships between water physicochemical variables and aquatic microbial communities. Water physicochemical variables in Waihai (south part) intensively changed along time, whereas those in Caohai did not. Photoautotrophic communities were significantly divergent between Caohai and Waihai. Waihai had a lower diversity of photoautotrophic community, containing higher abundance of Cyanophyceae (89.9%) than Caohai (42.7%). Nutrient level and Cyanophyceae only exhibited strong correlations in Wahai (p < 0.05). Redundancy analysis and microbial ecological network suggested that microbial communities in Caohai had a higher stability. Deterministic process dominated the microbial assembly (50-80% for bacteria and >90% for photoautotrophs), and particularly in Caohai. Our results unraveled that the structure and assembly of bacterial and photoautotrophic communities significantly changed after ecological restoration, offering valuable suggestions that photosynthetic diversity should be focused for other ecological restoration projects.

Ecological Dynamics and Co-occurrences Among Prokaryotes and Microeukaryotes in a Diatom Bloom Process in Xiangshan Bay, China

Diatom blooms can significantly affect the succession of microbial communities, yet little is known about the assembly processes and interactions of microbial communities during autumn bloom events. In this study, we investigated the ecological effects of an autumn diatom bloom on prokaryotic communities (PCCs) and microeukaryotic communities (MECs), focusing on their assembly processes and interactions. The PCCs were largely dominated by Alphaproteobacteria, Gammaproteobacteria, Cyanobacteria, and Flavobacteria, while the MECs primarily included Diatomea, Dinoflagellata, and Chlorophyta. The succession of both PCCs and MECs was mainly driven by this diatom bloom and environmental factors, such as nitrate and silicate. Null modeling revealed that homogeneous selection had a more pronounced impact on the structure of PCCs compared with that of MECs. In particular, drift and dispersal limitation cannot be neglected in the assembly processes of MECs. Co-occurrence network analyses showed that Litorimicrobium, Cercozoa, Marine Group I (MGI), Cryptomonadales, Myrionecta, and Micromonas may affect the bloom process. In summary, these results elucidated the complex, robust interactions and obviously distinct assembly mechanisms of PCCs and MECs during a diatom bloom and extend our current comprehension of the ecological mechanisms and microbial interactions involved in an autumn diatom bloom process.

Water quality and periphyton functional response to input of dissolved manure-derived hydrochars (DHCs)

Dissolved organic matter (DOM) plays a critical role in the global carbon cycle and provides food and energy for aquatic organisms. Recently, hydrochar, as a solid carbonaceous substance derived from hydrothermal carbonization, has been increasingly used as a soil amendment. Upon entering the soil, dissolved components (DHCs) were released from hydrochar as exogenous DOM, finally entering the aquatic ecosystems by runoff, which participates in environmental geochemical processes. However, relevant reports revealing the response of the aquatic ecosystem to the input of DHCs remain insufficiently elucidated. For the first time, the fundamental features of DHCs and their influence on water quality and aquatic biological function were investigated in this study. DHCs at 260 °C (DHC260) had lower yields, a greater [C/N], worse biodegradability, and larger humic acid relative amounts than did DHCs at 180 °C (DHC180). The DHC structural alterations in periphyton-incubated aquatic ecosystems suggested that protein substances were more easily degraded or assimilated by periphyton, especially for DHC180, with rates of decrease of 34.5-63.5%. The increased chemical oxygen demand (COD) degradation in the DHC260 treatments was most likely due to humic acid substances with higher COD equivalents. Furthermore, DHC260 caused phosphorus to accumulate in periphyton, reducing aquatic phosphorus concentration. Notably, the abundances of Flavobacteria and Cyanobacteria associated with water blooms increased 12.7-25.5- and 1.3-8.3-fold, respectively; consequently, the promotional impact of DHCs on algal blooms should be considered. This result extends the nonnegligible role of DHCs in aquatic ecosystems and underlines the need to regulate the hydrochar application process.

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.

Salinity Is a Key Determinant for the Microeukaryotic Community in Lake Ecosystems of the Inner Mongolia Plateau, China

The arid and semiarid areas experienced remarkable lake shrinkage during recent decades due to intensive human activities and climate change, which would result in unprecedented changes of microeukaryotic communities. However, little is known about how climate change affects the structure and ecological mechanisms of microeukaryotic communities in this area. Here, we used an 18S rRNA gene-based high-throughput sequencing approach to explore the structure, interspecies interaction, and assembly processes of the microeukaryotic community in lake ecosystems of the Inner Mongolia Plateau. As a direct result of climate change, salinity has become the key determinant of the lacustrine microeukaryotic community in this region. The microeukaryotic community in this ecosystem can be divided into three groups: salt (Lake Daihai), brackish (Lake Dalinuoer) and freshwater lakes. Co-occurrence network analysis revealed that salinity shapes the interspecies interactions of the microeukaryotic community. This causes interspecies interactions to change from antagonistic to cooperative with an increase in salinity. Phylogenetic-based β-nearest taxon distance analyses revealed that stochastic processes mainly dominated the microeukaryotic community assembly in lake ecosystems of the Inner Mongolia Plateau, and salinity stress drove the assembly processes of the microeukaryotic community from stochastic to deterministic. Overall, these findings expand the current understanding of interspecies interactions and assembly processes of microeukaryotic communities during climate change in lake ecosystems of the Inner Mongolia Plateau.

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.

Changes in protist communities in drainages across the Pearl River Delta under anthropogenic influence

Drainages in the Pearl River Delta urban agglomeration (PRDUA) host vital aquatic ecosystems and face enormous pressures from human activities in one of the largest urban agglomerations in the world. Despite being crucial components of aquatic ecosystems, the interactions and assembly processes of the protistan community are rarely explored in areas with serious anthropogenic disturbance. To elucidate the mechanisms of these processes, we used environmental DNA sequencing of 18S rDNA to investigate the influence of environmental factors and species interactions on the protistan community and its assembly in drainages of the PRDUA during summer. The protistan community showed a high level of diversity and a marked spatial pattern in this region. Community assembly was driven primarily by stochastic processes based on the Sloan neutral community model, explaining 74.28%, 75.82%, 73.67%, 74.40% and 51.24% of community variations in the BJ (Beijiang), XJ (Xijiang), PRD (Pearl River Delta), PRE (Pearl River Estuary) areas and in total, respectively. Meanwhile, environmental variables including temperature, pH, dissolved oxygen, transparency, nutrients and land use were strongly correlated with the composition and assembly of the protistan community, explaining 40.40% of variation in the protistan community. Furthermore, the bacterial community was simultaneously analysed by the 16S rDNA sequencing. Co-occurrence network analysis revealed that species interactions within bacteria (81.41% positive) or protists (82.80% positive), and those between bacteria and protists (50% positive and 50% negative) impacted the protistan community assembly. In summary, stochastic processes dominated, whereas species interactions and environmental factors also played important roles in shaping the protistan communities in drainages across the PRDUA. This study provides insights into the ecological patterns, assembly processes and species interactions underlying protistan dynamics in urban aquatic ecosystems experiencing serious anthropogenic disturbance.

Community assembly mechanisms and co-occurrence patterns of nitrite-oxidizing bacteria communities in saline soils

Nitrite-oxidizing bacteria (NOB) catalyze the second step of nitrification by oxidizing nitrite to nitrate, which is a key process in the biogeochemical nitrogen cycling. However, little is known about the co-occurrence patterns and assembly processes of NOB communities in agricultural soils with different salinities. Here, we explored the effects of salinity on Nitrobacter and Nitrospira community using high-throughput sequencing and multivariate statistical analyses. Our results showed that high salinity significantly inhibited the nitrite oxidation rates and decreased the abundance of Nitrobacter and Nitrospira. Extreme salty conditions significantly altered the diversity and composition of Nitrospira community but had little effect on Nitrobacter community. Nitrobacter network in high salinity soils was more closely connected while the connectivity of Nitrospira network became weak. Nitrobacter and Nitrospira community exhibited distinct assembly processes at different salinity levels. Stochastic processes were dominant in the Nitrobacter community assembly in both low and high salinity soils. Interestingly, the assembly of Nitrospira community was governed by stochastic and deterministic processes in low and high salinity soils, respectively. To our knowledge, our study provides the first description of the co-occurrence patterns and assembly processes of NOB communities in agricultural soils with different salinities. These results can help us understand the NOB ecological roles and improve the nitrite oxidation activity in a high salinity environment.

Biogeographical patterns and mechanisms of microbial community assembly that underlie successional biocrusts across northern China

Biocrusts play critical eco-functions in many drylands, however it is challenging to explore their community assembly, particularly within patched successional types and across climate zones. Here, different successional biocrusts (alga, lichen, and moss-dominated biocrusts) were collected across the northern China, and assembly of biocrust microbial communities was investigated by high-throughput sequencing combined with measurements of soil properties and microclimate environments. Bacterial and eukaryotic communities showed that the maximum and minimum community variation occurred across longitude and latitude, respectively. In the regions where all three stages of biocrusts were involved, the highest community difference existed between successional stages, and decreased with distance. The community assembly was generally driven by dispersal limitation, although neutral processes have controlled the eukaryotic community assembly in hyperarid areas. Along the succession, bacterial community had no obvious patterns, but eukaryotic community showed increasing homogeneity, with increased species sorting and decreased dispersal limitation for community assembly. Compared to early successional biocrusts, there were higher microbial mutual exclusions and more complex networks at later stages, with distinct topological features. Correlation analysis further indicated that the balance between deterministic and stochastic processes might be mediated by aridity, salinity, and total phosphorus, although the mediations were opposite for bacteria and eukaryotes.

Community Compositions of Phytoplankton and Eukaryotes during the Mixing Periods of a Drinking Water Reservoir: Dynamics and Interactions

In deep drinking water reservoir ecosystems, the dynamics and interactions of community compositions of phytoplankton and eukaryotes during the mixing periods are still unclear. Here, morphological characteristics combined with high-throughput DNA sequencing (HTS) were used to investigate the variations of phytoplankton and the eukaryotic community in a large canyon-shaped, stratified reservoir located at the Heihe River in Shaanxi Province for three months. The results showed that Bacillariophyta and Chlorophyta were the dominant taxa of the phytoplankton community, accounting for more than 97% of total phytoplankton abundance, which mainly consisted of Melosira sp., Cyclotella sp., and Chlorella sp., respectively. Illumina Miseq sequencing suggested that the biodiversity of eukaryotes increased over time and thatspecies distribution was more even. Arthropoda (6.63% to 79.19%), Ochrophyta (5.60% to 35.16%), Ciliophora (1.81% to 10.93%) and Cryptomonadales (0.25% to 11.48%) were the keystone taxa in common, contributing over 50% of the total eukaryotic community. Cryptomycota as a unique fungus was observed to possess significant synchronization with algal density, reaching a maximum of 10.70% in December (when the algal density distinctly decreased) and suggesting that it might affect the growth of algae through parasitism. Co-occurrence network patterns revealed the complicated and diverse interactions between eukaryotes and phytoplankton, suggesting that eukaryotes respond to variations in dynamic structure of the phytoplankton community, although there might be antagonistic or mutualistic interactions between them. Redundancy analysis (RDA) results showed that environmental variables collectively explained a 96.7% variance of phytoplankton and 96.3% variance of eukaryotic microorganisms, indicating that the temporal variations of phytoplankton and eukaryotic microorganisms were significantly affected by environmental conditions. This study shows that potential interactions exist between phytoplankton and eukaryotic microorganism communities, andcould improve our understanding of the ecological roles of phytoplankton and eukaryotic microorganisms in changing aquatic ecosystems. However, long-term investigations are necessary in order to obtain comprehensive understandings of their complicated associations.