Diversity and Spatial Distribution of Chromophytic Phytoplankton in the Bay of Bengal Revealed by RuBisCO Genes (rbcL)

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.

Environmental DNA metabarcoding reveals diverse phytoplankton assemblages and potentially harmful algal distribution along the urban coasts of Türkiye

Marine phytoplankton are widely used to monitor the state of the water column due to their rapid changes in response to environmental conditions. In this study, we aimed to investigate the coastal phytoplankton assemblages, including bloom-forming species using high-throughput sequencing of 18S rRNA genes targeting the V4 region and their relationship with environmental variables along the Istanbul coasts of the Sea of Marmara. A total of 118 genera belonging to six phyla were detected. Among them, Dinoflagellata (36) and Bacillariophyta (26) were represented with the highest number of genera. According to the relative abundance of DNA reads, the most abundant taxa were Dinoflagellata_phylum (18.1%), Emiliania (8.4%), Biecheleria (8.4), and Noctiluca (8.1%). The ANOSIM test showed that there was a significant temporal difference in the assemblages, while the driving environmental factors were pH, water temperature, and salinity. According to the TRIX index, the trophic state of the coasts was highly mesotrophic and eutrophic. In addition, 45 bloom-forming and HAB taxa were detected and two species of Noctiluca and Emiliania, which frequently cause blooms in the area, were recorded in high abundance. Our results provide insight into the phytoplankton assemblages along the urbanized coastlines by analysing the V4 region of 18S rRNA. This data can support future studies that use both traditional methods and metabarcoding, employing various primers and targeting different genes and regions.

The Bay of Bengal exposes abundant photosynthetic picoplankton and newfound diversity along salinity-driven gradients

The Bay of Bengal (BoB) is a 2,600,000 km2 expanse in the Indian Ocean upon which many humans rely. However, the primary producers underpinning food chains here remain poorly characterized. We examined phytoplankton abundance and diversity along strong BoB latitudinal and vertical salinity gradients-which have low temperature variation (27-29°C) between the surface and subsurface chlorophyll maximum (SCM). In surface waters, Prochlorococcus averaged 11.7 ± 4.4 × 104  cells ml-1 , predominantly HLII, whereas LLII and 'rare' ecotypes, HLVI and LLVII, dominated in the SCM. Synechococcus averaged 8.4 ± 2.3 × 104  cells ml-1 in the surface, declined rapidly with depth, and population structure of dominant Clade II differed between surface and SCM; Clade X was notable at both depths. Across all sites, Ostreococcus Clade OII dominated SCM eukaryotes whereas communities differentiated strongly moving from Arabian Sea-influenced high salinity (southerly; prasinophytes) to freshwater-influenced low salinity (northerly; stramenopiles, specifically, diatoms, pelagophytes, and dictyochophytes, plus the prasinophyte Micromonas) surface waters. Eukaryotic phytoplankton peaked in the south (1.9 × 104  cells ml-1 , surface) where a novel Ostreococcus was revealed, named here Ostreococcus bengalensis. We expose dominance of a single picoeukaryote and hitherto 'rare' picocyanobacteria at depth in this complex ecosystem where studies suggest picoplankton are replacing larger phytoplankton due to climate change.

Spatial and temporal changes in the assembly mechanism and co-occurrence network of the chromophytic phytoplankton communities in coastal ecosystems under anthropogenic influences

As a typical semiclosed coastal sea area in China, the ecological environment of Bohai Bay has been significantly disturbed by human activities. As primary producers, the chromophytic phytoplankton are the basis of ecosystems, especially in coastal ecosystems, and changes in the chromophytic phytoplankton community can affect the stability of the entire ecosystem. In this study, we investigated the effects of the human activity-induced spatial and temporal environmental heterogeneity on the community composition, diversity, assembly mechanisms, and co-occurrence networks of chromophytic phytoplankton in Bohai Bay during the wet season and the dry season. The results showed that in both seasons, there was obvious environmental heterogeneity between the nearshore area and the offshore area, and the nearshore areas were more affected by human disturbance. Although higher diversity was supported by the abundance of nutrients in nearshore areas, co-occurrence network analysis revealed that the chromophytic phytoplankton were less closely connected to each other in nearshore areas than in offshore areas due to chemical oxygen demand (COD), eutrophication index (EI), and dissolved inorganic nitrogen (DIN). The nearshore network was less stable than the offshore co-occurrence network in both seasons, which was related to the concentration of dissolved oxygen and COD. Both stochastic and deterministic processes dominated the assembly of the chromophytic phytoplankton communities, with different importance rankings of stochastic and deterministic processes in the nearshore and offshore areas. Drift dominated the assembly of the communities in nearshore areas, while variable selection dominated the assembly of the communities in offshore areas. DIN, EI, and COD, rather than geographic distance, were the main environmental factors affecting the phylogenetic turnover of the chromophytic phytoplankton. Our study showed that environmental heterogeneity caused by human disturbance had a greater impact on the chromophytic phytoplankton communities in Bohai Bay than natural factors such as temperature and salinity.

Does environmental DNA reflect the actual phytoplankton diversity in the aquatic environment? Case study of marine mucilage in the Sea of Marmara

The present study was designed to assess the effectiveness of the eDNA metabarcoding approach to determine the phytoplankton composition in the marine environment with a special focus on mucilage episodes in the Sea of Marmara. For this purpose, the samples were collected from 5 different sites located in the Sea of Marmara and the northern Aegean Sea during the mucilage episode in June 2021. The phytoplankton diversity was analyzed morphologically and by 18S rRNA gene amplicon sequencing, and the dataset of both methods was compared, accordingly. The results showed significant differences between methods in terms of composition and the abundance of the phytoplankton groups. While Miozoa was the most abundant group by metabarcoding, light microscopy (LM) indicated a dominance of Bacillariophyta. Katablepharidophyta was found at lower abundances by the metabarcoding (representing < 1% of the community); the members of this phylum were not observed by a microscope. At the lower taxonomic levels, Chaetoceros was the only genus detected in all samples by both methods. Additionally, while mucilage-forming Gonyaulax fragilis, Cylindrotheca closterium, and Thalassiosira rotula were detected to species-level by LM, metabarcoding was able to determine these organisms at the genus level. On the other hand, the genus Arcocellulus was found in all metabarcoding datasets and not detected by microscopy. The results indicated that metabarcoding can detect a greater number of genera and reveal taxa that were overlooked by light microscopy but to develop a complete picture of phytoplankton diversity in the sample, microscopical observations still are in need.

Deciphering the diversity and distribution of chromophytic phytoplankton in the Bohai Sea and the Yellow Sea via RuBisCO genes (rbcL)

Present study investigated composition and distribution of chromophytic phytoplankton in the Bohai Sea (BS) and the Yellow Sea (YS) by using rbcL genes. Bacillariophyceae, Haptophyceae and Pelagophyceae were the most abundant phytoplankton groups. Distinct phytoplankton communities were observed in the BS and the YS: offshore stations were dominated by bloom forming genera Thalassiosira and Skeletonema, while brown tide-forming species including Chrysochromulina spp. and Aureococcus anophagefferens were commonly found in the nearshore areas. Redundancy analysis showed that phosphate, temperature and silicic acid play key roles in structuring chromophytic phytoplankton, such as phytoplankton at nearshore stations were affected by nutrient runoff from adjacent rivers (Yellow River). Anthropogenic activities in the Bohai Sea and seasonal circulation of ocean currents may also contribute to shaping chromophytic phytoplankton communities. This study provides data support and foundational observations of chromophytic phytoplankton in the BS and the YS, and their responses to environmental gradients and human activities.

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.

Meta-analysis shows both congruence and complementarity of DNA and eDNA metabarcoding to traditional methods for biological community assessment

DNA metabarcoding is increasingly used for the assessment of aquatic communities, and numerous studies have investigated the consistency of this technique with traditional morpho‐taxonomic approaches. These individual studies have used DNA metabarcoding to assess diversity and community structure of aquatic organisms both in marine and freshwater systems globally over the last decade. However, a systematic analysis of the comparability and effectiveness of DNA‐based community assessment across all of these studies has hitherto been lacking. Here, we performed the first meta‐analysis of available studies comparing traditional methods and DNA metabarcoding to measure and assess biological diversity of key aquatic groups, including plankton, microphytobentos, macroinvertebrates, and fish. Across 215 data sets, we found that DNA metabarcoding provides richness estimates that are globally consistent to those obtained using traditional methods, both at local and regional scale. DNA metabarcoding also generates species inventories that are highly congruent with traditional methods for fish. Contrastingly, species inventories of plankton, microphytobenthos and macroinvertebrates obtained by DNA metabarcoding showed pronounced differences to traditional methods, missing some taxa but at the same time detecting otherwise overseen diversity. The method is generally sufficiently advanced to study the composition of fish communities and replace more invasive traditional methods. For smaller organisms, like macroinvertebrates, plankton and microphytobenthos, DNA metabarcoding may continue to give complementary rather than identical estimates compared to traditional approaches. Systematic and comparable data collection will increase the understanding of different aspects of this complementarity, and increase the effectiveness of the method and adequate interpretation of the results.

Dynamics of bacterial communities during a seasonal hypoxia at the Bohai Sea: Coupling and response between abundant and rare populations

Marine bacterial community plays a vital role in the formation of the hypoxia zone in coastal oceans. Yet, their dynamics in the seasonal hypoxia zone of the Bohai Sea (BHS) are barely studied. Here, the 16S rRNA gene-based high-throughput sequencing was used to explore the dynamics of their diversity, structure, and function as well as driving factors during the gradual deoxygenation process in the BHS. Our results evinced that the bacterial community was dominated by Proteobacteria, followed by Bacteroidetes, Firmicutes, Actinobacteria, and Cyanobacteria, etc. The abundant subcommunity dominated in the number of sequences (49%) while the rare subcommunity dominated in the number of species (99.61%). Although abundant subcommunity accounted for most sequences, rare subcommunity possessed higher diversity, richness and their population dramatically changed (higher turnover) during the hypoxia transition. Further, co-occurrence network analysis proved the vital role of rare subcommunity in the process of community assembly. Additionally, beta diversity partition revealed that both subcommunities possessed a higher turnover component than nestedness and/or richness component, implying species replacement could explain a considerable percentage of community variation. This variation might be governed by both environmental selection and stochastic processes, and further, it influenced the nitrogen cycle (PICRUSt-based prediction) of the hypoxia zone. Overall, this study provides insight into the spatial-temporal heterogeneity of bacterial and their vital role in biogeochemical cycles in the hypoxia zone of the BHS. These findings will extend our horizons about the stabilization mechanism, feedback regulation, and interactive model inside the bacterial community under oxygen-depleted ecosystems.

Paleo-diatom composition from Santa Barbara Basin deep-sea sediments: a comparison of 18S-V9 and diat-rbcL metabarcoding vs shotgun metagenomics

Sedimentary ancient DNA (sedaDNA) analyses are increasingly used to reconstruct marine ecosystems. The majority of marine sedaDNA studies use a metabarcoding approach (extraction and analysis of specific DNA fragments of a defined length), targeting short taxonomic marker genes. Promising examples are 18S-V9 rRNA (~121-130 base pairs, bp) and diat-rbcL (76 bp), targeting eukaryotes and diatoms, respectively. However, it remains unknown how 18S-V9 and diat-rbcL derived compositional profiles compare to metagenomic shotgun data, the preferred method for ancient DNA analyses as amplification biases are minimised. We extracted DNA from five Santa Barbara Basin sediment samples (up to ~11 000 years old) and applied both a metabarcoding (18S-V9 rRNA, diat-rbcL) and a metagenomic shotgun approach to (i) compare eukaryote, especially diatom, composition, and (ii) assess sequence length and database related biases. Eukaryote composition differed considerably between shotgun and metabarcoding data, which was related to differences in read lengths (~112 and ~161 bp, respectively), and overamplification of short reads in metabarcoding data. Diatom composition was influenced by reference bias that was exacerbated in metabarcoding data and characterised by increased representation of Chaetoceros, Thalassiosira and Pseudo-nitzschia. Our results are relevant to sedaDNA studies aiming to accurately characterise paleo-ecosystems from either metabarcoding or metagenomic data.

Finding needles in a haystack-Extensive diversity in the eustigmatophyceae revealed by community metabarcode analysis targeting the rbcL gene using lineage-directed primers

Sequences from the Stramenopile class Eustigmatophyceae are rarely reported in metabarcoding studies, and when they have been reported, there are very few haplotypes. We hypothesized that the paucity of eustigmatophyte species detected in these studies may be a result of the metabarcoding techniques used, which have primarily employed universal ribosomal RNA gene regions. In this study, we examined environmental DNA samples from 22 sites in southwestern Virginia, some of which had previously been studied using ribosomal RNA analysis. We used metabarcoding techniques targeting the plastid rbcL gene with new primers designed to produce a 370 bp amplicon from all lineages of the Eustigmatophyceae in a reference collection. The amplicons were then analyzed with DADA2 to produce amplicon sequence variants (ASVs). Our results revealed 184 rbcL haplotypes that can be tentatively assigned to the Eustigmatophyceae from these sites, representing much higher diversity than has been detected by ribosomal DNA-based studies. The techniques employed can be used for future studies of population structure, ecology, distribution, and diversity of this class. With these techniques, it should be possible to make realistic estimates of the species-level diversity of the Eustigmatophyceae on local, regional, and perhaps global scales.

Assembly Processes and Co-occurrence Patterns of Abundant and Rare Bacterial Community in the Eastern Indian Ocean

Microbial communities are composed of many rare species and a few abundant species. Considering the disproportionate importance of rare species for ecosystem functioning, it is important to understand the mechanisms structuring the rare and abundant components of a diverse community in response to environmental changes. Here, we used a 16S ribosomal RNA gene sequencing approach to investigate the bacterial community diversity in the Eastern Indian Ocean (EIO) during the monsoon and intermonsoon. We employed a phylogenetic null model and network analysis to evaluate the assembly processes and co-occurrence pattern of the microbial community. We found that higher bacterial diversity was detected in the intermonsoon with high temperature and low Chlorophyll a concentrations and N/P ratios. The balance between ecological deterministic processes and stochastic processes varied with seasons in the EIO. Meanwhile, conditionally rare taxa (CRT) were more likely modulated by variable selection processes than always rare taxa (ART) and abundant taxa (AT) (CRT > ART > AT). By linking assembly process and species co-occurrence, we demonstrated that the microbial co-occurrence associations tended to be higher when deterministic processes (mainly variable selection) were weaker. This negative trend was observed in rare species rather than abundant species. The linkage could enhance our understanding of the underlying mechanisms underpinning the generation and maintenance of microbial community diversity.

Carbonate chemistry variability in the southern Yellow Sea and East China Sea during spring of 2017 and summer of 2018

Marginal seas are highly productive and disproportionately large contributors to global air-sea CO₂ fluxes. Due to complex physical and biogeochemical conditions, the southern Yellow-East China Sea is an ideal site for studying carbonate chemistry variability. The carbonate system was investigated in the area in spring of 2017 and summer of 2018. Dissolved inorganic carbon (DIC) and total alkalinity (TA) concentrations were higher in the SYS than the ECS due to material from carbonate weathering and erosion carried by the Yellow River. High pH and low DIC and TA were observed in the Zhe-Min Coastal Current in spring due to high primary productivity caused by Changjiang River input and the Taiwan Warm Current. Temperature and biological activity were the primary drivers controlling the partial pressure of CO₂ (pCO₂) in the SYS, pCO₂ was controlled by primary productivity related to nutrients carried by the Changjiang River and physical mixing in the Changjiang River plume and inner/middle shelves of the ECS, whereas temperature was the dominant factor determining pCO₂ distributions in the ECS outer shelf waters influenced by the Kuroshio Current. Overall, the entire study area shifted from a CO₂ sink (-4.18 ± 5.60 mmol m^-2 d^-1) to a weak source (1.02 ± 4.87 mmol m^-2 d^-1) from spring to summer. Specifically, the SYS and ECS offshore waters changed from CO₂ sinks in spring to sources in summer, while the Changjiang River plume was always a CO₂ sink.

pr2-primers: An 18S rRNA primer database for protists

Metabarcoding of microbial eukaryotes (collectively known as protists) has developed tremendously in the last decade, almost solely relying on the 18S rRNA gene. As microbial eukaryotes are extremely diverse, many primers and primer pairs have been developed. To cover a relevant and representative fraction of the protist community in a given study system, an informed primer choice is necessary, as no primer pair can target all protists equally well. As such, a smart primer choice is very difficult even for experts and there are very few online resources available to list existing primers. We built a database listing 285 primers and 83 unique primer pairs that have been used for eukaryotic 18S rRNA gene metabarcoding. In silico performance of primer pairs was tested against two sequence databases: PR² version 4.12.0 for eukaryotes and a subset of silva version 132 for bacteria and archaea. We developed an R-based web application enabling browsing of the database, visualization of the taxonomic distribution of the amplified sequences with the number of mismatches, and testing any user-defined primer or primer set (https://app.pr2-primers.org). Taxonomic specificity of primer pairs, amplicon size and location of mismatches can also be determined. We identified universal primer sets that matched the largest number of sequences and analysed the specificity of some primer sets designed to target certain groups. This tool enables guided primer choices that will help a wide range of researchers to include protists as part of their investigations.

The Composition and Primary Metabolic Potential of Microbial Communities Inhabiting the Surface Water in the Equatorial Eastern Indian Ocean

Simple Summary

Marine microbes are regarded as the most diverse organisms in the biosphere and drive biogeochemical cycles through their metabolism. It is essential to understand the structure and metabolic function of microbial communities. The Indian Ocean is the third largest ocean in the world, and it possesses unique hydrographical properties. So far, assessments of microbial diversity and metabolism need to be improved in the Indian Ocean. Therefore, we carried out a series of investigations in the equatorial eastern Indian Ocean in order to clarify the local microbial communities and detect the genetic potential for microbial functions. The obtained results suggested Cyanobacteria was the dominant microbial group, and predicted the Calvin cycle and the assimilatory nitrate and nitrite reduction played important role in the pathway of carbon fixation and nitrogen metabolism respectively. This study provides insights into microbial community structures as well as the metabolic potential that may be active in the local environment, and lays the groundwork for understanding the roles of microbes in energy and resource cycling in this habitat.

Abstract

Currently, there is scant information about the biodiversity and functional diversity of microbes in the eastern Indian Ocean (EIO). Here, we used a combination of high-throughput sequencing of 16S rRNA genes and a metagenomic approach to investigate the microbial population structure and its metabolic function in the equatorial EIO. Our results show that Cyanobacterial Prochlorococcus made up the majority of the population. Interestingly, there were fewer contributions from clades SAR11 (Alphaproteobacteria) and SAR86 (Gammaproteobacteria) to microbial communities than contributions from Prochlorococcus. Based on functional gene analysis, functional genes rbcL, narB, and nasA were relatively abundant among the relevant genes. The abundance of Prochlorococcus implies its typically ecological adaptation in the local ecosystem. The microbial metabolic potential shows that in addition to the main carbon fixation pathway Calvin cycle, the rTCA cycle and the 3-HP/4-HB cycle have potential alternative carbon fixation contributions to local ecosystems. For the nitrogen cycle, the assimilatory nitrate and nitrite reduction pathway is potentially the crucial form of nitrogen utilization; unexpectedly, nitrogen fixation activity was relatively weak. This study extends our knowledge of the roles of microbes in energy and resource cycling in the EIO and provides a foundation for revealing profound biogeochemical processes driven by the microbial community in the ocean.

Distribution of Chromophytic Phytoplankton in the Eddy-Induced Upwelling Region of the West Pacific Ocean Revealed Using rbcL Genes

Marine chromophytic phytoplankton are a diverse group of algae and contribute significantly to the total oceanic primary production. However, the spatial distribution of chromophytic phytoplankton is understudied in the West Pacific Ocean (WPO). In this study, we have investigated the community structure and spatial distribution of chromophytic phytoplankton using RuBisCO genes (Form ID rbcL). Our results showed that Haptophyceae, Pelagophyceae, Cyanophyceae, Xanthophyceae, and Bacillariophyceae were the dominant groups. Further, chromophytic phytoplankton can be distinguished between upwelling and non-upwelling zones of the WPO. Surface and 75 m depths of a non-upwelling area were dominated by Prochlorococcus strains, whereas chromophytic phytoplankton were homogenously distributed at the surface layer in the upwelling zone. Meanwhile, Pelagomonas-like sequences were dominant at DCM (75 m) and 150 m depths of the upwelling zone. Non-metric multidimensional scaling (NMDS) analysis did not differentiate between chromophytic phytoplankton in the upwelling and non-upwelling areas, however, it showed clear trends of them at different depths. Further, redundancy analysis (RDA) showed the influence of physicochemical parameters on the distribution of chromophytic phytoplankton. Along with phosphate (p < 0.01), temperature and other dissolved nutrients were important in driving community structure. The upwelling zone was impacted by a decrease in temperature, salinity, and re-supplement of nutrients, where Pelagomonas-like sequences outnumbered other chromophytic groups presented.

Influence of the cold bottom water on taxonomic and functional composition and complexity of microbial communities in the southern Yellow Sea during the summer

The formation and presence of the cold bottom water (Yellow Sea Cold Water Mass, YSCWM) is a striking hydrological phenomenon in the southern Yellow Sea during the summer and has important effects on the marine ecosystem. To better understand its influence on microbial community structure and function, we compared the bacterial, archaeal and microeukaryotic communities in the cold water mass area (CWMA) and the southern area (SA) during the summer using amplicon and metagenomic sequencings. The habitat environment in the deep waters of the CWMA was characterized by higher salinity/DO/PO₄-P, greater depth/distance to the coast, and lower levels of temperature/chlorophyll a/DIN/SiO₃-Si/N:P ratio compared to that of the SA. Pure depth or distance to the coast explained a small portion of the microbial community variance, while environment explained a significant fraction of the variance when partialling the effects of depth and distance to the coast. Oligotrophic taxa (e.g. SAR11 clade Ia, Nitrosopumilus, Chloropicophyceae) dominated the deep water communities in the CWMA, while the common coastal taxa (e.g. Roseobacter strain HIMB11, Bacillariophyta, Noctilucophyceae) were more dominant in the deep waters of the SA, suggesting the great impact of the oligotrophic condition in the YSCWM on microbial communities. The microbial co-occurrence networks in the CWMA were less complex but contained a higher proportion of mutual exclusion relationship among prokaryotes; the prokaryotic α-diversity in the CWMA was significantly lower than in the SA while the microeukaryotic α-diversity was significantly higher in the CWMA, implying that prokaryotes and microeukaryotes respond to the cold water mass differently and the competition among prokaryotes was intensified under the impact of the YSCWM. Genes that relate to replication and repair accounted for a significantly lower proportion in the CWMA, which was likely an adaptation to the low carbon environment.