Composition and spatial-temporal dynamics of phytoplankton community shaped by environmental selection and interactions in the Jiaozhou Bay

Biogeographic patterns of meio- and micro-eukaryotic communities in dam-induced river-reservoir systems

Although the Three Gorges Dam (TGD) is the world's largest hydroelectric dam, little is known about the spatial-temporal patterns and community assembly mechanisms of meio- and micro-eukaryotes and its two subtaxa (zooplankton and zoobenthos). This knowledge gap is particularly evident across various habitats and during different water-level periods, primarily arising from the annual regular dam regulation. To address this inquiry, we employed mitochondrial cytochrome c oxidase I (COI) gene-based environmental DNA (eDNA) metabarcoding technology to systematically analyze the biogeographic pattern of the three communities within the Three Gorges Reservoir (TGR). Our findings reveal distinct spatiotemporal characteristics and complementary patterns in the distribution of meio- and micro-eukaryotes. The three communities showed similar biogeographic patterns and assembly processes. Notably, the diversity of these three taxa gradually decreased along the river. Their communities were less shaped by stochastic processes, which gradually decreased along the longitudinal riverine-transition-lacustrine gradient. Hence, deterministic factors, such as seasonality, environmental, and spatial variables, along with species interactions, likely play a pivotal role in shaping these communities. Environmental factors primarily drive seasonal variations in these communities, while hydrological conditions, represented as spatial distance, predominantly influence spatial variations. These three communities followed the distance-decay pattern. In winter, compared to summer, both the decay and species interrelationships are more pronounced. Taken together, this study offers fresh insights into the composition and diversity patterns of meio- and micro-eukaryotes at the spatial-temporal level. It also uncovers the mechanisms behind community assembly in various environmental niches within the dam-induced river-reservoir systems. KEY POINTS: • Distribution and diversity of meio- and micro-eukaryotes exhibit distinct spatiotemporal patterns in the TGR. • Contribution of stochastic processes in community assembly gradually decreases along the river. • Deterministic factors and species interactions shape meio- and micro-eukaryotic community.

Intensified effect of nitrogen forms on dominant phytoplankton species succession by climate change

Nutrient proportion, light intensity, and temperature affect the succession of dominant phytoplankton species. Despite these insights, this transformation mechanism in highly turbid lakes remains a research gap, especially in response to climate change. To fill this gap, we investigated the mechanism by which multi-environmental factors influence the succession of dominant phytoplankton species in Lake Chagan. This investigation deployed the structural equation model (SEM) and the hydrodynamic-water quality-water ecology mechanism model. Results demonstrated that the dominant phytoplankton species in Lake Chagan transformed from diatom to cyanobacteria during 2012 and 2022. Notably, Microcystis was detected in 2022. SEM revealed the primary environment variables for this succession, including water temperature (Tw), nutrients (total nitrogen (TN), total phosphorus (TP), and ammonia nitrogen (NH4N)), and total suspended solids (TSS). Moreover, this event was not the consequence of zooplankton grazing. An integrated hydrodynamic-water quality-bloom mechanism model was built to explore the mechanism driving phytoplankton succession and its response to climate change. Nutrients determined the phytoplankton biomass and dominant species succession based on various proportions. High NH4N:NO3N ratios favored cyanobacteria and inhibited diatom under high TSS. Additionally, the biomass proportions of diatom (30.77 % vs. 22.28 %) and green (30.56 % vs. 23.30 %) decreased dramatically. In contrast, cyanobacteria abundance remarkably increased (35.78 % to 51.71 %) with the increasing NH4-N:NO3-N ratios. In addition, the proportion of non-nitrogen-fixing cyanobacteria was higher than that of the nitrogen-fixing cyanobacteria counterparts when TN:TP≥20 and NH4N:NO3N ≥ 10. Light-limitation phenotypes also experienced an increase with the rising NH4N:NO3N ratios. Notably, the cyanobacterial biomass reached 3-6 times that in the baseline scenario when the air temperature escalated by 3.0 °C until 2061 under the SSP585 scenario. We highlighted the effect of nitrogen forms on the succession of dominant phytoplankton species. Climate warming will increase nitrogen proportion, providing an insightful reference for controlling cyanobacterial blooms.

Composition and temporal dynamics of the phytoplankton community in Laizhou Bay revealed by microscopic observation and rbcL gene sequencing

Laizhou Bay, a major breeding ground for economic marine organisms in the northern waters of China, is facing rapid environmental degradation. In this study, field surveys in this area were conducted in the spring, summer, and autumn of 2020. Microscopic observation and RuBisCO large subunit (rbcL) gene analysis were employed to understand the community structure and temporal dynamics of phytoplankton. The phytoplankton community structures detected by the two methods showed significant differences. Microscopic observation revealed the dominance of dinoflagellates in spring that shifted to the dominance of diatoms in summer and autumn. However, rbcL gene sequencing consistently identified diatoms as dominant throughout all three seasons, with their relative abundance showing an increasing trend. Conversely, the relative abundance of the second- and third-most abundant taxa, namely, haptophytes and ochrophytes, decreased as the seasons transitioned. rbcL gene sequencing annotated more species than microscopy. It could detect haptophytes and cryptophytes, which were overlooked by microscopy. In addition, rbcL gene sequencing detected a remarkable amount of Thalassiosira profunda, which was previously unidentified in this sea area. However, it appeared to underestimate the contribution of dinoflagellates considerably, with most taxa being only identified through microscopic identification. The two methods jointly identified 28 harmful algal bloom taxa with similar detection quantities but substantial differences in species composition. Phytoplankton communities were influenced by temperature, salinity, and nutrients. The results of this work suggest that a combination of multiple techniques is necessary for a comprehensive understanding of phytoplankton.

Comparative assessment of the intragenomic variations of dinoflagellate Tripos species through single-cell sequencing

Tripos is a large dinoflagellate genus widely distributed in the world's oceans. Morphology-based species identification is inconclusive due to high morphological intraspecific variability. Metabarcoding analysis has been demonstrated to be effective for species identification and tracking their spatiotemporal dynamics. However, accumulating evidence suggests high levels of intragenomic variations (IGVs) are common in many algae, leading to concerns about overinterpretation of molecular diversity in metabarcoding studies. In this project, we evaluated and compared IGVs in Tripos species by conducting the first high-throughput sequencing (HTS) of 18S rDNA V4 of Tripos single cells. High numbers of haplotypes (19-172) were identified in each of the 30 Tripos cells. Each cell contained one dominant haplotype with high relative abundance and many haplotypes with lower abundances. Thus, the presence of multiple minor haplotypes substantially overestimate the molecular diversity identified in metabarcoding analysis, which encompass not only interspecific and intraspecific diversities, but high levels of IGVs.

Harmonized coexistence of intragenomic variations in diatom Skeletonema strains

Metabarcoding analysis has been demonstrated to be an effective technology for monitoring diversity and dynamics of phytoplankton including Skeletonema species. Although molecular diversity uncovered in metabarcoding projects has generally been interpreted as sum of interspecies diversity and intraspecies diversity, accumulating evidence suggests that it also harbors unprecedentedly high levels of intra-genomic variations (IGVs). As up to thousands of amplicon sequence variants (ASVs) identified in a typical metabarcoding project can be annotated to be Skeletonema species, we hypothesize that substantial portions of these ASVs are contributed by IGVs. Here, the nature of IGVs in Skeletonema species was quantitatively analyzed by carrying out single-strain metabarcoding analysis of 18S rDNA V4 in 49 strains belonging to seven Skeletonema species. Results showed that each Skeletonema strain harbored a high level of IGVs as expected. While many Skeletonema strains each contained one dominant ASV and a substantial number of ASVs displaying much lower relative abundance, other Skeletonema strains each contained multiple ASVs with comparable or nearly equally abundances. Thus the co-existence of multiple dominant ASVs in a single cell indicated a tug-of-war of these variants in evolution, which may eventually result in harmonized coexistence of multiple dominant ASVs. A total of nine dominant ASVs and 652 non-dominant ASVs were found in 49 strains of seven Skeletonema species, indicating rich interspecies and intraspecies variations, and complex evolution of IGVs in genus of Skeletonema. The results confirmed that the extensive degree of IGVs was the main contributor to the high molecular diversity revealed by metabarcoding analysis. This study highlights the importance of quantitative characterization of IGVs in Skeletonema species for accurate interpretation of species diversity in metabarcoding analysis.

Assessment of coastal eco-environmental sustainable development under multiple pressures: A case study of Jiaozhou Bay, China

The increasing demand for economic development in the coastal zone is exerting greater pressure on the eco-environment in this region. Therefore, it is imperative to conduct comprehensive research on the sustainable development of the eco-environment in coastal areas. Grounded in sustainable development theory and Integrated Coastal Zone Management (ICZM), this study assessed sustainable development of Jiaozhou Bay (JZB), China over the period from 2015 to 2019 under multiple pressures applying the Driver-Pressure-State-Impact-Response (DPSIR) framework. Subsequently, a SWOT (Strengths, Weaknesses, Opportunities and Threats) analysis was employed to propose specific countermeasures and recommendations. In terms of weight within the criterion layer, Response (R) and Pressure (P) were assigned higher weights. The comprehensive assessment value Z in the target layer indicated a slight decrease in 2016 followed by a gradual increase that eventually led to stabilization. The comprehensive assessment values of the indicators varied depending on different conditions, implying that higher values would contribute to sustainable development of the eco-environment. Based on a SWOT analysis, strategies for Qingdao were derived from a combination of research findings and future development policy. These strategies include preventing and controlling land and sea pollution, governing reclamation activities, optimizing industrial structure, and strengthening public opinion guidance. This study highlights the application of quantitative and qualitative analysis methods in ICZM, providing specific recommendations based on policy analysis for achieving sustainable development in JZB's eco-environment under multiple pressures.

Ecological impacts of the expansion of offshore wind farms on trophic level species of marine food chain

The global demand for renewable energy has resulted in a rapid expansion of offshore wind farms (OWFs) and increased attention to the ecological impacts of OWFs on the marine ecosystem. Previous reviews mainly focused on the OWFs' impacts on individual species like birds, bats, or mammals. This review collected numerous field-measured data and simulated results to summarize the ecological impacts on phytoplankton, zooplankton, zoobenthos, fishes, and mammals from each trophic level and also analyze their interactions in the marine food chain. Phytoplankton and zooplankton are positively or adversely affected by the 'wave effect', 'shading effect', oxygen depletion and predation pressure, leading to a ± 10% fluctuation of primary production. Although zoobenthos are threatened transiently by habitat destruction with a reduction of around 60% in biomass in the construction stage, their abundance exhibited an over 90% increase, dominated by sessile species, due to the 'reef effect' in the operation stage. Marine fishes and mammals are to endure the interferences of noise and electromagnetic, but they are also aggregated around OWFs by the 'reef effect' and 'reserve effect'. Furthermore, the complexity of marine ecosystem would increase with a promotion of the total system biomass by 40% through trophic cascade effects strengthen and resource partitioning alternation triggered by the proliferation of filter-feeders. The suitable site selection, long-term monitoring, and life-cycle-assessment of ecological impacts of OWFs that are lacking in current literature have been described in this review, as well as the carbon emission and deposition.

Chromosome-level genome assembly of marine diatom Skeletonema tropicum

Skeletonema tropicum is a marine diatom of the genus Skeletonema that also includes many well-known species including S. marinoi. S. tropicum is a high temperature preferring species thriving in tropical ocean regions or temperate ocean regions during summer-autumn. However, mechanisms of ecological adaptation of S. tropicum remain poorly understood due partially to the lack of a high-quality whole genome assembly. Here, we report the first high-quality chromosome-scale genome assembly for S. tropicum, using cutting-edge technologies including PacBio single molecular sequencing and high-throughput chromatin conformation capture. The assembled genome has a size of 78.78 Mb with a scaffold N50 of 3.17 Mb, anchored to 23 pseudo-chromosomes. In total, 20,613 protein-coding genes were predicted, of which 17,757 (86.14%) genes were functionally annotated. Collinearity analysis of the genomes of S. tropicum and S. marinoi revealed that these two genomes were highly homologous. This chromosome-level genome assembly of S. tropicum provides a valuable genomic platform for comparative analysis of mechanisms of ecological adaption.

Extensive intragenomic variations of the 18S rDNA V4 region in the toxigenic diatom species Pseudo-nitzschia multistriata revealed through high-throughput sequencing

Metabarcoding analysis is an effective technique for monitoring the domoic acid-producing Pseudo-nitzschia species in marine environments, uncovering high-levels of molecular diversity. However, such efforts may result in the overinterpretation of Pseudo-nitzschia species diversity, as molecular diversity not only encompasses interspecies and intraspecies diversities but also exhibits extensive intragenomic variations (IGVs). In this study, we analyzed the V4 region of the 18S rDNA of 30 strains of Pseudo-nitzschia multistriata collected from the coasts of China. The results showed that each P. multistriata strain harbored about a hundred of unique 18S rDNA V4 sequence varieties, of which each represented by a unique amplicon sequence variant (ASV). This study demonstrated the extensive degree of IGVs in P. multistriata strains, suggesting that IGVs may also present in other Pseudo-nitzschia species and other phytoplankton species. Understanding the scope and levels of IGVs is crucial for accurately interpreting the results of metabarcoding analysis.

A Pseudo-nitzschia metabarcoding approach with a calibrated ITS1 reference sequence database applied in the Taiwan Strait

Pseudo-nitzschia is a cosmopolitan phytoplankton genus of which some species can form blooms and produce the neurotoxin domoic acid (DA). Identification of Pseudo-nitzschia is generally based on field material or strains followed by morphological and/or molecular characterization. However, this process is time-consuming and laborious, and can not obtain a relatively complete and reliable profile of the Pseudo-nitzschia community, because species with low abundance in the field or potentially unavailable for culturing may easily be overlooked. In the present study, specific ITS primer sets were designed and evaluated using in silico matching. The primer set ITS-84F/456R involving the complete ITS1 region was found optimal. Based on matching with a Pseudo-nitzschia ITS1 reference sequence database carefully-calibrated in this study, a metabarcoding approach using annotated amplicon sequence variants (ASV) was applied in the Taiwan Strait of the East China Sea during two cruises in the spring and summer of 2019. In total, 48 Pseudo-nitzschia species/phylotypes including 36 known and 12 novel were uncovered, and verified by haplotype networks, ITS2 secondary structure comparisons and divergence analyses. Correlation analyses revealed that temperature was a key factor affecting the seasonal variation of the Pseudo-nitzschia community. This study provides an overview of the Pseudo-nitzschia community in the Taiwan Strait, with new insights into the diversity. The developed metabarcoding approach may be used elsewhere as a standard reference for accurate annotation of Pseudo-nitzschia.

Temperature and nutrients alter the relative importance of stochastic and deterministic processes in the coastal macroinvertebrates biodiversity assembly on long-time scales

Macroinvertebrates play a vital role in coastal ecosystems and are an important indicator of ecosystem quality. Both anthropogenic activity and environmental changes may lead to significant changes in the marine macroinvertebrate community. However, the assembly process of benthic biodiversity and its mechanism driven by environmental factors at large scales remains unclear. Here, using the benthic field survey data of 15 years at large spatial and temporal scales from the Yellow Sea Large Marine Ecosystem, we investigated the relative importance of environmental selection, dispersal processes, random-deterministic processes of macroinvertebrates community diversity assembly, and the responses of this relative importance driven by temperature and nutrients. Results showed that the macroinvertebrates community diversity is mainly affected by dispersal. Nitrogen and phosphorus are the most important negative factors among environmental variables, while geographical distance is the main limiting factor of β diversity. Within the range of 0.35-0.70 mg/L of nutrients, increasing nutrient concentration can significantly facilitate the contribution of the decay effect to β diversity. Within the temperature range studied (15.0-18.0°C), both warming and cooling can lead to a greater tendency for species diversity assembly processes to be dominated by deterministic processes. The analysis contributes to a better understanding of the assembly process of the diversity of coastal marine macroinvertebrates communities and how they adapt to global biogeochemical processes.

Remote sensing retrieval and driving analysis of phytoplankton density in the large storage freshwater lake: A study based on random forest and Landsat-8 OLI

Remote sensing monitoring of seasonal changes in phytoplankton density and analyses of the driving factors of phytoplankton densities are necessary for assessing the health of aquatic ecosystems, controlling lake eutrophication, and formulating ecological restoration policies. Building upon the satellite-ground synchronization experiment that involves the in situ aquatic ecological monitoring conducted in Nansi Lake, which is the largest storage lake situated along the eastern route of the South-to-North Water Diversion Project, we developed a phytoplankton density retrieval model utilizing the random forest (RF) method and Landsat-8 OLI data. On this basis, we mapped the seasonal fluctuations and spatial disparities in the phytoplankton densities from 2013 to 2023. Subsequently, we conducted a detailed analysis of the driving factors and considered both the natural and anthropogenic aspects. The results indicate that (1) the RF model, when utilizing three band combinations, yielded favorable results with R2, RMSE and MAE values of 0.67, 1.31 × 106 cells/L and 1.18 × 106 cells/L, respectively. (2) The phytoplankton densities exhibited both seasonal and spatial variations, with higher concentrations in summer and autumn than in spring and winter. Significantly, the northwestern region of Zhaoyang Lake and the southeastern region of Weishan Lake had substantially greater phytoplankton densities than did the other areas. Furthermore, overarching upward trends were observed from 2013 to 2023, reflecting an annual rate of increase of 3.32%. (3) An analysis of the causal factors indicated that temperatures and gross agricultural production levels are the primary drivers influencing the seasonal variations and distributions of phytoplankton densities. In the future, we will delve into the potential of deep learning and utilize various satellite sensors to explore the intricacies of phytoplankton monitoring, as well as the complex mechanisms that influence aquatic ecological health.

The high molecular diversity in Noctiluca scintillans is dominated by intra-genomic variations revealed by single cell high-throughput sequencing of 18S rDNA V4

The application of high-throughput sequencing (HTS) technologies has revolutionized research on phytoplankton biodiversity by generating an unprecedented amount of molecular data in marine ecosystem surveys. However, high-level of molecular diversity uncovered in HTS-based metabarcoding analyses may lead to overinterpretation of phytoplankton diversity due to excessive intra-genomic variations (IGVs). The aims in this study are to explore the nature of phytoplankton molecular diversity and to test the hypothesis. We carried out single-cell metabarcoding analysis of 18S rDNA V4 sequences obtained in single Noctiluca scintillans cells isolated from various sites in coastal waters of China. Results showed that each single N. scintillans cell harbored a high level of IGVs with about 100 amplicon sequence variants (ASVs). The large numbers of non-dominant ASVs identified in N. scintillans cells, which might correspond to the larger numbers of ASVs annotated as N. scintillans and showed similar temporal dynamics in metabarcoding analyses, could inflate the inter-species diversity or intra-species genetic diversity. In addition, there were large numbers of additional ASVs that were not annotated as N. scintillans. These non-N. scintillans ASVs might represent diverse preys for N. scintillans, consistent with previous reports that N. scintillans may act as chance predator of a broad-spectrum preys. This single-cell study has unambiguously demonstrated that the existence of high levels of IGVs in N. scintillans and most likely many other phytoplankton species, demonstrating that the majority of the molecular diversity revealed in metabarcoding analysis, which were generally interpreted as the sum of inter-species diversity and intra-species diversity, actually included high levels of IGVs and should be interpreted with caution.

Rich diversity and active spatial-temporal dynamics of Thalassiosira species revealed by time-series metabarcoding analysis

Thalassiosira is a species-rich genus in Bacillariophyta that not only contributes positively as primary producer, but also poses negative impacts on ecosystems by causing harmful algal blooms. Although taxonomical studies have identified a large number of Thalassiosira species, however, the composition of Thalassiosira species and their geographical distribution in marine ecosystems were not well understood due primarily to the lack of resolution of morphology-based approaches used previously in ecological expeditions. In this study, we systematically analyzed the composition and spatial-temporal dynamic distributions of Thalassiosira in the model marine ecosystem Jiaozhou Bay by applying metabarcoding analysis. Through analyzing samples collected monthly from 12 sampling sites, 14 Thalassiosira species were identified, including five species that were not previously reported in Jiaozhou Bay, demonstrating the resolution and effectiveness of metabarcoding analysis in ecological research. Many Thalassiosira species showed prominent temporal preferences in Jiaozhou Bay, with some displaying spring-winter preference represented by Thalassiosira tenera, while others displaying summer-autumn preference represented by Thalassiosira lundiana and Thalassiosira minuscula, indicating that the temperature is an important driving factor in the temporal dynamics. The application of metabarcoding analysis, equipped with appropriate molecular markers with high resolution and high specificity and databases of reference molecular marker sequences for potential all Thalassiosira species, will revolutionize ecological research of Thalassiosira species in Jiaozhou Bay and other marine ecosystems.

Stochasticity-driven weekly fluctuations distinguished the temporal pattern of particle-associated microorganisms from its free-living counterparts in temperate coastal seawater

Microbial community dynamics directly determine their ecosystem functioning. Despite the well-known annual recurrence pattern, little is known how different lifestyles affect the temporal variation and how community assembly mechanisms change over different temporal scales. Here, through a high-resolution observation of size fractionated samples over 60 consecutive weeks, we investigate the distinction in weekly distribution pattern and assembly mechanism between free-living (FL) and particle-associated (PA) communities in highly dynamic coastal environments. A clear pattern of annual recurrence was observed, which was more pronounced in FL compared to PA, resulting in higher temporal specificity in the former samples. Both the two size fractions displayed significant temporal distance-decay patterns, yet the PA community showed a higher magnitude of community variation between adjacent weeks, likely caused by sudden, drastic and long-lived blooms of heterotrophic bacteria. Generally, determinism (environmental selection) had a greater effect on the community assembly than stochasticity (random birth, death, and dispersal events), with significant contributions from temperature and inorganic nutrients. However, a clear shift in the temporal assembly pattern was observed, transitioning from a prevalence of stochastic processes driving short-term (within a month) fluctuations to a dominance of deterministic processes over longer time intervals. Between adjacent weeks, stochasticity was more important in the community assembly of PA than FL. This study revealed that stochastic processes can lead to rapid, dramatic and irregular PA community fluctuations, indicating weak resistance and resilience to disturbances, which considering the role of PA microbes in carbon processing would significantly affect the coastal carbon cycle. Our results provided a new insight into the microbial community assembly mechanisms in the temporal dimension.

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.

Environmental forcing of phytoplankton carbon-to-diversity ratio and carbon-to-chlorophyll ratio: A case study in Jiaozhou Bay, the Yellow Sea

The relationships between phytoplankton carbon (C) biomass and diversity (i.e., C-to-H' ratio) and chlorophyll a (i.e., C-to-Chl a ratio) are good indicators of marine ecosystem functioning and stability. Here we conducted four cruises spanning 2 years in Jiaozhou Bay to explore the dynamics of C-to-H' and C-to-Chl a ratios. The results showed that the phytoplankton C biomass and diversity were dominated by diatoms, followed by dinoflagellates. The average C-to-H' ratio ranged from 84.10 to 912.17, with high values occurring in the northern region of the bay. In contrast, the average C-to-Chl a ratio ranged between 15.55 and 89.47, and high values primarily appeared in the northern or northeastern part of the bay. In addition, the redundancy analysis showed that temperature and phosphate (DIP) were significantly correlated with both ratios in most cases, indicating that temperature and DIP may be key factors affecting the dynamics of C-to-H' and C-to-Chl a ratios.

Spatio-temporal dynamics of phytoplankton in a diversion reservoir and the major influencing factors: taxonomic versus functional groups classification

Identifying factors affecting phytoplankton dynamics is crucial to the management of aquatic ecosystems. A lot of scholars have conducted intensive studies on phytoplankton in lake or reservoirs, but not many studies have been conducted on diversion reservoirs. To explore the seasonal and spatial variation of phytoplankton communities and their relationship with environmental factors in the context of water diversion, a case study was carried out at XiKeng (XK) reservoir in South China. In this study, month-by-month water samples and phytoplankton were collected from this reservoir from December, 2021, to July, 2022. The results showed that the phytoplankton community was characterized by significant spatial and temporal variations. There were significant differences in phytoplankton abundance and structure in the reservoirs in terms of time. The abundance of phytoplankton cells and the proportion of Cyanobacteria in the reservoir showed a trend of increasing from autumn to spring and then decreasing from spring to summer, while the functional group evolved from S1 in autumn to SN in spring and summer. The abundance of phytoplankton was influenced by the dynamic water division and the characteristics of the reservoir itself, resulting in a spatial distribution characteristic of AIII > AII > AI. Water temperature (WT) and nutrients were the key factors driving the changes in phytoplankton abundance and community structure in the reservoir. These findings will deepen our understanding of the spatial and temporal dynamics of phytoplankton community structure in diversion reservoirs and provide a basis for freshwater water ecological management strategies.

Expansion of photoreception-related gene families may drive ecological adaptation of the dominant diatom species Skeletonema marinoi

Diatom species of the genus Skeletonema are dominant in global coastal waters with important roles in marine primary production and global biogeochemical cycling. Many Skeletonema species have been extensively studied also because they can cause harmful algae blooms (HABs) with negative impacts on marine ecosystems and aquaculture. In this study, the first chromosome-level assembly of the genome of Skeletonema marinoi was constructed. The genome size was 64.99 Mb with a contig N50 of 1.95 Mb. Up to 97.12 % of contigs were successfully anchored on 24 chromosomes. Analysis of the annotated genes revealed 28 large syntenic blocks with 2397 collinear gene pairs in the genome of S. marinoi, suggesting large-scale segmental duplication events in evolution. Substantial expansion of light-harvesting genes encoding fucoxanthin-chlorophyll a/c binding proteins, as well as expansion of photoreceptor gene families encoding aureochromes and cyptochromes (CRY) in S. marinoi were found, which may have shaped ecological adaptation of S. marinoi. In conclusion, the construction of the first high-quality Skeletonema genome assembly offers valuable clues on the ecological and evolutionary characteristics of this dominant coastal diatom species.

Harmful algal blooms in China: History, recent expansion, current status, and future prospects

The impacts of harmful algal blooms (HABs) on economies, public health, ecosystems, and aquaculture across the globe have all increased in recent decades, and this has been acutely the case in China. Here, we review the history of HABs and HABs research in China, as well as recent trends in HABs and future prospects of HAB science in China. The most updated analyses demonstrated that the number of HAB events, the number of HAB species, the aerial coverage of HABs, and the impacts of HABs in Chinese waters during the 21st century were all higher than that during the last two decades of the 20th century. The increase in the number of HABs in China has been significantly correlated with the increased discharge of ammonium and total phosphorus into coastal waters (p < 0.01 for both). Notable newly recognized events this century have included chronic HABs caused by Prorocentrum donghaiense and Karenia mikimotoi, a paralytic shellfish poisoning event caused by Gymnodinium catenatum that sickened 80 people, brown tides caused by Aureococcus anophagefferens, green tides caused by Ulva prolifera, golden tides caused by Sargassum horneri, and the disruption of a nuclear power plant caused by a bloom of Phaeocystis globosa. A series of key discoveries regarding HABs has been made this century including documentation of nearly all known HAB toxins in Chinese waters, discovery of novel cyst-formation and/or life stages of multiple HABs-causing species, identification of the chemical and physical oceanographic drivers of multiple HABs including those formed by P. donghaiense, K. mikimotoi, and U. prolifera, and the successful mitigation of HABs via the use of modified clay approaches. Future research prospects highlighted include the use of macroalgae as a means to prevent, mitigate, and control (PCM) HABs and the process by which multi-disciplinary studies involving molecular approaches (omics), remote in situ detection, artificial intelligence, and mega-data analyses might be used to develop refined and realistic HAB forecasting platforms. Collectively, this review demonstrates the significant evolution of HAB science since the 20th century in China and demonstrates that while HABs in China are complex and widespread, recent and on-going discoveries make the development of detailed understanding and effective measures to mitigate the negative effects of HABs a hopeful outcome in the coming years.

Phytoplankton community structure and environmental factors during the outbreak of Crown-of-Thorns Starfish in Xisha Islands, South China Sea

The "larval starvation hypothesis" proposed that the growing frequency of Crown-of-Thorns Starfish (CoTS) outbreaks could be attributed to increased availability of phytoplankton. However, comprehensive field investigation on the living environment of CoTS larvae and the availability of phytoplankton are still lacking. A cruise was conducted in June 2022 in Xisha Islands, South China Sea, to study the interaction between environmental conditions and phytoplankton communities during CoTS outbreak period. The average concentrations of dissolved inorganic phosphorus (0.05 ± 0.01 μmol L-1), dissolved inorganic nitrogen (0.66 ± 0.8 μmol L-1) and chlorophyll a (0.05 ± 0.05 μg L-1) suggested that phytoplankton may be limited for CoTS larvae in Xisha Islands. Microscopic observation and high-throughput sequencing were used to study the composition and structure of the phytoplankton communities. Bacillariophyta predominated in phytoplankton communities with the highest abundance and species richness. 29 dominant species, including 4 species with size-range preferred by CoTS larvae, were identified in Xisha Islands. The diversity index of all stations indicated a species-rich and structure-stable phytoplankton community in Xisha Islands during the period of CoTS outbreak, which may contribute to CoTS outbreak. These findings revealed the structure of phytoplankton community and environmental factors in the study area during CoTS outbreak, providing the groundwork for future research into the causes and processes of CoTS outbreak.

Phytoplankton Communities Miniaturization Driven by Extreme Weather in Subtropical Estuary under Climate Changes

Estuarine ecosystems are subject to extreme weather and climatic events, particularly global climate change. To characterize the effects of extreme weather, such as heavy precipitation and cold waves, on phytoplankton communities, four seasonal cruises were conducted in the subtropical Pearl River Estuary. Among the main phytoplankton taxa, small (pico- and nano-sized) phytoplankton accounted for approximately 50% and 70% in the upstream estuary. In spring, summer, and autumn, small-sized phytoplankton communities were dominated by Chlorophyta (20-45%), Ochrophyta (Bacillariophyceae; 25-40%), and Dinoflagellata (20-25%), associated with high phytoplankton network stability and river discharge. In winter, small phytoplankton communities were dominated by Cryptophyta (50%), Dinoflagellata (25%), and Ochrophyta (Bacillariophyceae; 10%), which were associated with low network stability and temperature. Low light and high nutrient conditions induced by large river discharge favored the growth of the smallest known brackish Chlorophyta, such as Ostreococcus tauri. Low temperatures provided a competitive advantage for the survival of the small-sized cryptophyte genus Teleaulax, especially in the upstream estuary during the 2020/2021 extreme cold wave period. Our findings highlight the impact of extreme weather induced by climatic events on the miniaturization of phytoplankton communities in subtropical estuaries by altering temperature, light availability, and nutrient dynamics.

Drivers of the spatiotemporal patterns of the mangrove crab metacommunity in a tropical bay

Revealing community patterns and driving forces is essential in community ecology and a prerequisite for effective management and conservation efforts. However, the mangrove ecosystem and its important fauna group such as the crabs, still lack multi-processes research under metacommunity framework, resulting in evidence and theorical application gaps. To fill these gaps, we selected China's most representative mangrove bay reserve in tropical zone as a stable experimental system and conducted a seasonal investigation (July 2020, October 2020, January 2021, and April 2021) of mangrove crabs. We performed a multi-approach analysis using both pattern-based and mechanistic method to distinguish the processes driving the mangrove crab metacommunity. Our results showed that the crab metacommunity exhibits a Clementsian pattern in the bay-wide mangrove ecosystem but is influenced by both local environmental heterogeneity and spatial processes, thus representing a combined paradigm of species sorting and mass effect. Moreover, the long-distance spatial constraints are more pronounced compared to the local environmental factors. This is reflected in the greater importance of the broad-scale Moran's Eigenvector Maps, the distance-decay pattern of similarity, and the difference in beta diversity dominated by the turnover component. This pattern changes throughout the year, mainly due to changes in dominant functional groups caused by the stress of changes in water salinity and temperature induced by air temperature and precipitation. This research provides multi-dimension research data and relevant analysis, offering clear evidence for understanding the patterns and related driving forces of crab metacommunity in tropical bay mangroves, and verifies the applicability of some general laws in the system. Future studies can address more diverse spatiotemporal scales, gaining a clearer understanding to serve the conservation of mangrove ecosystems and economically important fishery species.

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.

Diversity and Distribution of Harmful Algal Bloom Species from Seamount to Coastal Waters in the South China Sea

Mount Xianbei is one of the largest shallow seamounts located in the middle of the South China Sea (SCS), which might play a role in shaping the biodiversity of surrounding continental coastal waters, particularly the diversity of phytoplankton species causing frequent harmful algal blooms (HABs) in northern SCS. However, the diversity, composition, and distribution of phytoplankton species in the seamount regions of Xianbei remain largely unexplored. In this study, samples around and outside the seamount regions were collected during a late summer cruise of 2021 to test whether seamounts play a role in HAB species propagation. In total, we identified 19 HAB species across all samples using the ASV-based DNA metabarcoding approach, 6 of which had not been reported previously in the SCS, suggesting a diverse HAB species in the SCS. Specifically, 16 HAB species were found in the seamount region of Xianbei, and 5 of them were also found in the coastal waters, indicating a close connection between seamount and coastal waters. This study was the first attempt to explore HAB species' spatial diversity and vertical distribution in the seamount region of Xianbei at single-nucleotide resolution, which provides a novel explanation for the coastal HAB occurrence in the northern SCS. IMPORTANCE There are a number of seamounts under the water of the South China Sea (SCS). The seamounts might play a role in shaping the biodiversity of surrounding continental coastal waters. However, there is no direct evidence revealing the relationship of the biodiversity of phytoplankton between seamounts and coastal waters in the SCS, especially those species having the potential to form harmful algal blooms (HABs). Some HAB species might proliferate in certain geographic locations, while others may be broadly distributed across oceanic provinces. In this study, we provided a detailed analysis of phytoplankton composition and molecular detection of HAB species from seamount to coastal waters in the SCS, which suggested a strong interaction in the HAB species between the two areas. This finding provides new insights into the diversity and distribution of HABs in seamounts and their role in shaping the composition and the occurrence of HABs in coastal water.

Karenia selliformis bloom dynamics and growth rate estimation in the Sfax harbour (Tunisia), by using automated flow cytometry equipped with image in flow, during autumn 2019

A Karenia selliformis bloom event in the Gulf of Gabès (Mediterranean Sea), was monitored over 9 days at high frequency during fall 2019, by using an automated flow cytometer (Cytosense, Cytobuoy b.v.) with an image-in-flow attachment. The instrument recorded the shape of the optical signals that lead to the resolution of six cell groups of pico-, nano- and microphytoplankton, during the Harmful Algal Bloom (HAB). K. selliformis cell dimensions derived from the hourly records, enabled to estimate the daily division rate over the bloom period. Results revealed that K. selliformis was the only bloom-forming species and it reached its highest mean abundance the fourth day of the survey. A shift in the nutrient composition occurred with a potential P limitation during the bloom growth and N limitation during the bloom collapse. The co-inertia analysis revealed opposite patterns for K. selliformis and heterotrophic prokaryotes suggesting trophic interactions and possible mixotrophic behaviour of K. selliformis at the end of the bloom. K. selliformis exhibited low growth rates generally < 1 division day^-1, which could not explain the observed high abundance. The tide played a crucial role in the dynamics of K. selliformis at a semi-diurnal scale and at spring-neap tide scale and was probably enhancing K. selliformis accumulation.