Estuarine gradients dictate spatiotemporal variations of microbiome networks in the Chesapeake Bay

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

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

Planktonic eukaryotes in the Chesapeake Bay: contrasting responses of abundant and rare taxa to estuarine gradients

Phytoplankton are important drivers of aquatic ecosystem function and environmental health. Their community compositions and distributions are directly impacted by environmental processes and human activities, including in the largest estuary in North America, the Chesapeake Bay. It is crucial to uncover how planktonic eukaryotes play fundamental roles as primary producers and trophic links and sustain estuarine ecosystems. In this study, we investigated the detailed community structure and spatiotemporal variations of planktonic eukaryotes in the Chesapeake Bay across space and time for three consecutive years. A clear seasonal and spatial shift of total, abundant, and rare planktonic eukaryotes was evident, and the pattern recurred interannually. Multiple harmful algal species have been identified in the Bay with varied distribution patterns, such as Karlodinium, Heterosigma akashiwo, Protoperidinium sp., etc. Compared to abundant taxa, rare subcommunities were more sensitive to environmental disturbance in terms of richness, diversity, and distribution. The combined effects of temporal variation (13.3%), nutrient availability (10.0%), and spatial gradients (8.8%) structured the distribution of eukaryotic microbial communities in the Bay. Similar spatiotemporal patterns between planktonic prokaryotes and eukaryotes suggest common mechanisms of adjustment, replacement, and species interaction for planktonic microbiomes under strong estuarine gradients. To our best knowledge, this work represents the first systematic study on planktonic eukaryotes in the Bay. A comprehensive view of the distribution of planktonic microbiomes and their interactions with environmental processes is critical in understanding the underlying microbial mechanisms involved in maintaining the stability, function, and environmental health of estuarine ecosystems.

IMPORTANCE

Deep sequencing analysis of planktonic eukaryotes in the Chesapeake Bay reveals high community diversity with many newly recognized phytoplankton taxa. The Chesapeake Bay planktonic eukaryotes show distinct seasonal and spatial variability, with recurring annual patterns of total, abundant, and rare groups. Rare taxa mainly contribute to eukaryotic diversity compared to abundant groups, and they are more sensitive to spatiotemporal variations and environmental filtering. Temporal variations, nutrient availability, and spatial gradients significantly affect the distribution of eukaryotic microbial communities. Similar spatiotemporal patterns in prokaryotes and eukaryotes suggest common mechanisms of adjustment, substitution, and species interactions in planktonic microbiomes under strong estuarine gradients. Interannually recurring patterns demonstrate that diverse eukaryotic taxa have well adapted to the estuarine environment with a long residence time. Further investigations of how human activities impact estuarine planktonic eukaryotes are critical in understanding their essential ecosystem roles and in maintaining environmental safety and public health.

Evaluating top-down, bottom-up, and environmental drivers of pelagic food web dynamics along an estuarine gradient

Identification of the key biotic and abiotic drivers within food webs is important for understanding species abundance changes in ecosystems, particularly across ecotones where there may be strong variation in interaction strengths. Using structural equation models (SEMs) and four decades of integrated data from the San Francisco Estuary, we investigated the relative effects of top-down, bottom-up, and environmental drivers on multiple trophic levels of the pelagic food web along an estuarine salinity gradient and at both annual and monthly temporal resolutions. We found that interactions varied across the estuarine gradient and that the detectability of different interactions depended on timescale. For example, for zooplankton and estuarine fishes, bottom-up effects appeared to be stronger in the freshwater upstream regions, while top-down effects were stronger in the brackish downstream regions. Some relationships (e.g., bottom-up effects of phytoplankton on zooplankton) were seen primarily at annual timescales, whereas others (e.g., temperature effects) were only observed at monthly timescales. We also found that the net effect of environmental drivers was similar to or greater than bottom-up and top-down effects for all food web components. These findings can help identify which trophic levels or environmental factors could be targeted by management actions to have the greatest impact on estuarine forage fishes and the spatial and temporal scale at which responses might be observed. More broadly, this study highlights how environmental gradients can structure community interactions and how long-term data sets can be leveraged to generate insights across multiple scales.

The Microbial Ecology of Estuarine Ecosystems

Human civilization relies on estuaries, and many estuarine ecosystem services are provided by microbial communities. These services include high rates of primary production that nourish harvests of commercially valuable species through fisheries and aquaculture, the transformation of terrestrial and anthropogenic materials to help ensure the water quality necessary to support recreation and tourism, and mutualisms that maintain blue carbon accumulation and storage. Research on the ecology that underlies microbial ecosystem services in estuaries has expanded greatly across a range of estuarine environments, including water, sediment, biofilms, biological reefs, and stands of seagrasses, marshes, and mangroves. Moreover, the application of new molecular tools has improved our understanding of the diversity and genomic functions of estuarine microbes. This review synthesizes recent research on microbial habitats in estuaries and the contributions of microbes to estuarine food webs, elemental cycling, and interactions with plants and animals, and highlights novel insights provided by recent advances in genomics.

Depth shapes microbiome assembly and network stability in the Mariana Trench

IMPORTANCE

Exploring microbial interactions and their stability/resilience from the surface to the hadal ocean is critical for further understanding of the microbiome structure and ecosystem function in the Mariana Trench. Vertical gradients did not destabilize microbial communities after long-term evolution and adaption. The uniform niche breadth, diversity, community complexity, and stability of microbiomes in both upper bathypelagic and hadal waters suggest the consistent roles of microbiomes in elemental cycling and adaptive strategies to overcome extreme environmental conditions. Compared with microeukaryotes, bacteria and archaea play a pivotal role in shaping the stability of the hadal microbiome. The consistent co-occurrence stability of microbiomes across vertical gradients was observed in the Mariana Trench. These results illuminate a key principle of microbiomes inhabiting the deepest trench: although distinct microbial communities occupy specific habitats, the interactions within microbial communities remain consistently stable from the upper bathypelagic to the hadal waters.

Environmental gradients shape microbiome assembly and stability in the East China sea

Microbiomes play a key role in marine ecosystem functioning and sustainability. Their organization and stability in coastal areas, particularly in anthropogenic-influenced regions, however, remains unclear compared with an understanding of how microbial community shifts respond to marine environmental gradients. Here, the assembly and community associations across vertical and horizontal gradients in the East China Sea are systematically researched. The seawater microbial communities possessed higher robustness and lower fragmentation and vulnerability compared to the sediment microbiomes. Spatial gradients act as a deterministic filtering factor for microbiome organization. Microbial communities had lower phylogenetic distance and higher niche breadth in the nearshore and offshore areas compared to intermediate areas. The phylogenetic distance of microbiomes decreased from the surface to the bottom but the niche breadth was enhanced in surface and bottom environments. Vertical gradients destabilized microbial associations, while the community diversity was enhanced. Multivariate regression tree analysis and canonical correspondence analysis indicated that depth, distance from shore, nutrient availability, temperature, salinity, and chlorophyll a, affected the distribution and co-occurrence of microbial groups. Our results highlight the crucial roles of environmental gradients in determining microbiome association and stability. These results improve our understanding of the survival strategies/adaptive mechanisms of microbial communities in response to environmental variation and provide new insights for protecting the ecosystems and maintaining the sustainability of ecological functions.

Deciphering characterization of seasonal variations in microbial communities of marine ranching: Diversity, co-occurrence network patterns, and assembly processes

Offshore coastal marine ranching ecosystems are one of the most productive ecosystems. The results showed that the composition and structure of the microbial communities varied considerably with the season. Co-occurrence network analysis demonstrated that the microbial network was more complex in summer and positively correlated links (cooperative or symbiotic) were dominated in autumn and winter. Null model indicated that the ecological processes of the bacterial communities were mainly governed by deterministic processes (mainly homogeneous selection) in summer. For microeukaryotic communities, assembly processes were more regulated by stochastic processes in all seasons. For rare taxa, assembly processes were regulated by stochastic processes and were not affected by seasonality. Changes in water temperature due to seasonal variations were the main, but not the only, environmental factor driving changes in microbial communities. This study will improve the understanding of offshore coastal ecosystems through the perspective of microbial ecology.

Distinct patterns of distribution, community assembly and cross-domain co-occurrence of planktonic archaea in four major estuaries of China

BACKGROUND

Archaea are key mediators of estuarine biogeochemical cycles, but comprehensive studies comparing archaeal communities among multiple estuaries with unified experimental protocols during the same sampling periods are scarce. Here, we investigated the distribution, community assembly, and cross-domain microbial co-occurrence of archaea in surface waters across four major estuaries (Yellow River, Yangtze River, Qiantang River, and Pearl River) of China cross climatic zones (~ 1,800 km) during the winter and summer cruises.

RESULTS

The relative abundance of archaea in the prokaryotic community and archaeal community composition varied with estuaries, seasons, and stations (reflecting local environmental changes such as salinity). Archaeal communities in four estuaries were overall predominated by ammonia-oxidizing archaea (AOA) (aka. Marine Group (MG) I; primarily Nitrosopumilus), while the genus Poseidonia of Poseidoniales (aka. MGII) was occasionally predominant in Pearl River estuary. The cross-estuary dispersal of archaea was largely limited and the assembly mechanism of archaea varied with estuaries in the winter cruise, while selection governed archaeal assembly in all estuaries in the summer cruise. Although the majority of archaea taxa in microbial networks were peripherals and/or connectors, extensive and distinct cross-domain associations of archaea with bacteria were found across the estuaries, with AOA as the most crucial archaeal group. Furthermore, the expanded associations of MGII taxa with heterotrophic bacteria were observed, speculatively indicating the endogenous demand for co-processing high amount and diversity of organic matters in the estuarine ecosystem highly impacted by terrestrial/anthropogenic input, which is worthy of further study.

CONCLUSIONS

Our results highlight the lack of common patterns in the dynamics of estuarine archaeal communities along the geographic gradient, expanding the understanding of roles of archaea in microbial networks of this highly dynamic ecosystem.

Response of planktonic microbial assemblages to disturbance in an urban sub-tropical estuary

Microbes are sensitive indicators of estuarine processes because they respond rapidly to dynamic disturbance events. As most of the world's population lives in urban areas and climate change-related disturbance events are becoming more frequent, estuaries bounded by cities are experiencing increasing stressors, at the same time that their ecosystem services are required more than ever. Here, using a multidisciplinary approach, we determined the response of planktonic microbial assemblages in response to seasonality and a rainfall disturbance in an urban estuary bounded by Australia's largest city, Sydney. We used molecular barcoding (16S, 18S V4 rRNA) and microscopy-based identification to compare microbial assemblages at locations with differing characteristics and urbanisation histories. Across 142 samples, we identified 8,496 unique free-living bacterial zOTUs, 8,175 unique particle associated bacterial zOTUs, and 1,920 unique microbial eukaryotic zOTUs. Using microscopy, we identified only the top <10% abundant, larger eukaryotic taxa (>10 µm), however quantification was possible. The site with the greater history of anthropogenic impact showed a more even community of associated bacteria and eukaryotes, and a significant increase in dissolved inorganic nitrogen following rainfall, when compared to the more buffered site. This coincided with a reduced proportional abundance of Actinomarina and Synechococcus spp., a change in SAR 11 clades, and an increase in the eukaryotic microbial groups Dinophyceae, Mediophyceae and Bathyoccocaceae, including a temporary dominance of the harmful algal bloom dinoflagellate Prorocentrum cordatum (syn. P. minimum). Finally, a validated hydrodynamic model of the estuary supported these results, showing that the more highly urbanised and upstream location consistently experienced a higher magnitude of salinity reduction in response to rainfall events during the study period. The best abiotic variables to explain community dissimilarities between locations were TDP, PN, modelled temperature and salinity (r = 0.73) for the free living bacteria, TP for the associated bacteria (r = 0.43), and modelled temperature (r = 0.28) for the microbial eukaryotic communities. Overall, these results show that a minor disturbance such as a brief rainfall event can significantly shift the microbial assemblage of an anthropogenically impacted area within an urban estuary to a greater degree than a seasonal change, but may result in a lesser response to the same disturbance at a buffered, more oceanic influenced location. Fine scale research into the factors driving the response of microbial communities in urban estuaries to climate related disturbances will be necessary to understand and implement changes to maintain future estuarine ecosystem services.

Distinct community assembly processes and habitat specialization driving the biogeographic patterns of abundant and rare bacterioplankton in a brackish coastal lagoon

The ecological diversity patterns and community assembly processes along spatio-temporal scales are least studied in the bacterioplankton sub-communities of brackish coastal lagoons. We examined the biogeographic patterns and relative influences of different assembly processes in structuring the abundant and rare bacterioplankton sub-communities of Chilika, the largest brackish water coastal lagoon of India. Rare taxa demonstrated significantly higher α- and β-diversity and biogeochemical functions than abundant taxa in the high-throughput 16S rRNA gene sequence dataset. The majority of the abundant taxa (91.4 %) were habitat generalists with a wider niche breadth (niche breadth index, B = 11.5), whereas most of the rare taxa (95.2 %) were habitat specialists with a narrow niche breadth (B = 8.9). Abundant taxa exhibited a stronger distance-decay relationship and higher spatial turnover rate than rare taxa. β-diversity partitioning revealed that the contribution of species turnover (72.2-97.8 %) was greater than nestedness (2.2-27.8 %) in causing the spatial variation in both abundant and rare taxa. Null model analyses revealed that the distribution of abundant taxa was mostly structured by stochastic processes (62.8 %), whereas deterministic processes (54.1 %) played a greater role in the rare taxa. However, the balance of these two processes varied across spatio-temporal scales in the lagoon. Salinity was the key deterministic factor controlling the variation of both abundant and rare taxa. Potential interaction networks showed a higher influence of negative interactions, indicating that species exclusion and top-down processes played a greater role in the community assembly. Notably, abundant taxa emerged as keystone taxa across spatio-temporal scales, suggesting their greater influences on other bacterial co-occurrences and network stability. Overall, this study provided detailed mechanistic insights into biogeographic patterns and underlying community assembly processes of the abundant and rare bacterioplankton over spatio-temporal scales in a brackish lagoon.

The underappreciated role of biotic factors in controlling the bloom ecology of potentially harmful microalgae in the Hudson-Raritan Bay

Despite widespread distribution of harmful algal blooms (HABs) and new and improved methods for detecting and quantifying them, no unifying ecological explanation has been found. Improved understanding depends upon local, ecological studies that include analysis of phytoplankton species data in relation to both abiotic and biotic parameters. Ecological network analysis was used to detect co-occurrence patterns among abiotic and biotic parameters in a long-term monitoring dataset (i.e., 2010-2021) from the eutrophic Hudson-Raritan Estuary (HRE) between the states of New York and New Jersey. The regular co-occurrence of potentially harmful bloom-forming species with companion species observed through microscopy was supported by the results of ecological network analysis, which showed that there were far more associations between HAB species and biotic parameters (∼95%) than abiotic parameters (∼5%). Temperature was the environmental variable that was most associated with HAB species throughout the estuary. The numerous network associations of HAB species with one another and with diatoms, dinoflagellates, and zooplankton highlight the complexity of planktonic food webs and interactions. Results also suggest that some taxa may play a central role in structuring the HRE plankton communities. These findings demonstrate that biotic associations play an underappreciated role in plankton structure and the value of examining the ecology of HAB species within the breadth of their biological communities. While network analysis does not fully explain and confirm complex associations among species, it does provide fresh insights and testable hypotheses to strengthen understanding and improve prediction.

Assembly and Network Stability of Planktonic Microorganisms under the Influence of Salinity Gradient: an Arctic Case Study from the Lena River Estuary to the Laptev Sea

The diversity and primary productivity in the Arctic ecosystem are rapidly changing due to global warming. Microorganisms play a vital role in biogeochemical cycling. However, the diversity of planktonic microorganism communities in the Laptev Sea, one of the most important marginal seas of the Western Arctic Ocean, have not been studied sufficiently in depth. The diversity and community structure of the planktonic microorganisms in the surface water were investigated at 20 stations on the Lena River flowing into the Laptev Sea. Multivariate statistical analyses demonstrated clear spatial patterns in the α diversity and community structure for microorganisms under different salinity levels. Co-occurrence networks of microbial communities revealed that spatial variation promoted differentiation of the characteristics and stability of microbial networks in the Laptev Sea. Contrary to expectations, abundant taxa were found to not have a large influence on the stability and resilience of microbial interactions in the region. On the contrary, less-abundant taxa were found to have far greater influence. The stability and resilience of the prokaryotic and microeukaryotic networks in the Lena River estuary and the continental shelf provided valuable insights into the impact of freshwater and land inflow disturbances on microbial assemblage. Overall, these results enhance our understanding of the composition of microbial communities and provide insights into how spatial changes of abundant versus rare species alter the nature and stability of microbial networks from the Lena River estuary to the Laptev Sea. In addition, this study explored microbial interactions and their ability to resist future disturbances. IMPORTANCE The regime of the Laptev Sea depends closely on the runoff of the Lena River. Microorganisms are essential components of aquatic food webs and play a significant role in polar ecosystems. In this study, we provided a basic microbial data set as well as new insights into the microbial networks from the Lena River estuary to the Laptev Sea, while exploring their potential to resist future disturbances. A comprehensive and systematic study of the community structure and function of the planktonic microorganisms in the Laptev Sea would greatly enhance our understanding of how polar microbial communities respond to the salinity gradient under climate warming.

The biogeography and co-occurrence network patterns of bacteria and microeukaryotes in the estuarine and coastal waters

Community and diversity shifts of bacteria and microeukaryotes with strong environmental and spatial variations have been unveiled in the Pearl River Estuary (PRE) and northern coastal part of South China Sea (SCS). However, it is not clear what the determining factors shape the microbial community and how the biotic interactions respond to the estuarine and oceanic environment. Here, we established the multiple regression models (MRM) and co-occurrence networks on microbial communities in PRE and SCS habitats. The results showed that there were significant differences of the abiotic factors affecting the bacterial and microeukaryotic communities between PRE and SCS habitats. Salinity explained the largest variations to the microbial community dissimilarities in PRE. Whereas spatial and environmental factors determined the microbial community dissimilarities in SCS. Positive relations between parasitic lineages (e.g. Perkinsea and Cercozoa) and algal taxa (Dinophyceae, Cryptophyta, Chlorophyta and Ochrophyta) dominated in the PRE network. While parasites Syndiniales positively correlated with other Syndiniales and protists in SCS. Strong positive associations among autotrophic and heterotrophic groups were revealed in both niches. Therefore, the biotic interactions are also important and may be responsible for the unexplained variations of the abiotic factors from MRM models. Microbial network in the PRE estuarine water had weakened resistance to environmental disturbances, while the SCS network had greater capacity to maintain network stability. This study shed light on the different mechanisms of abiotic and biotic factors in shaping the compositions of bacteria and microeukaryotes between PRE and SCS niches, and highlights the weakening effect of environmental disturbances on the microbial network stability.

The role of microbial communities on primary producers in aquatic ecosystems: Implications in turbidity stress resistance

Expanding the stress tolerance and adaptation potential of primary producers is of importance for the restoration and management of aquatic ecosystems. Microorganisms have been reported to mediate improved resistance toward different abiotic stresses of primary producers in terrestrial and marine ecosystems. However, it is not clear about the role of microbial communities in the turbidity resistance of primary producers, when aquatic ecosystems are under turbidity pressure. In this study, key microbes and the action path which enhance turbidity tolerance of primary producers were recognized by mesocosm and various multivariate statistical methods. Remarkable decrease of the biomass of primary producers was found with the increase of turbidity. Significant differences in microbial community under different turbidity pressure were recognized and key microbes which may expand the turbidity tolerance of primary producers were further identified. Rhodobacter and Rhodoferax were selected as key microbes by the investigation of keystone species in the microbial ecological network and significant discriminant taxa under different turbidity stress. The action path for microbial communities to help primary producers cope with turbidity pressure was found through structural equation model, and in which the increase of key microbes may expand the turbidity tolerance of primary producers through enhancing the microbial loop. The results may provide a new insight for aquatic ecosystems to resist turbidity stress, and provide a theoretical basis for the management and restoration of aquatic ecosystems.

TiO2 nanoparticles exert an adverse effect on aquatic microbial communities

Titanium dioxide nanoparticle (n-TiO₂) is a widely used nanomaterial, which is inevitably released as a residue into aquatic ecosystems during material production and usage. However, the effects of n-TiO₂ on aquatic microbial communities have not been completely elucidated. This study examined the toxic effects of n-TiO₂ on eukaryotic and prokaryotic microbial communities in freshwater environments. We determined that n-TiO₂ had a greater inhibitory effect on the growth of eukaryotic algae than cyanobacteria in monocultures. A similar phenomenon was observed in a microcosm experiment, revealing that n-TiO₂ slightly reduced the content of chlorophyll-a but evidently increased the phycocyanin content. Moreover, the alpha diversity of the eukaryotic community was not affected, whereas its beta diversity increased with exposure to n-TiO₂. Although n-TiO₂ altered the composition of freshwater microbial communities, it did not change the functions of the prokaryotic community, which might be attributed to the functional redundancy of microbiota. Co-occurrence network analysis indicated that n-TiO₂ destabilized the freshwater community, especially the eukaryotic community, and potentially disturbed the aquatic ecosystem. Our study revealed that the ecological risk of n-TiO₂ on aquatic microbial communities is complex; hence, rational utilization of n-TiO₂ should be emphasized.

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.

Contrasting Community Composition and Co-Occurrence Relationships of the Active Pico-Sized Haptophytes in the Surface and Subsurface Chlorophyll Maximum Layers of the Arctic Ocean in Summer

Haptophytes (Hacrobia: Haptophyta), which can perform phototrophic, phagotrophic, or mixotrophic nutritional modes, are critical for element cycling in a variety of aquatic ecosystems. However, their diversity, particularly in the changing Arctic Ocean (AO), remains largely unknown. In the present study, the biodiversity, community composition, and co-occurrence networks of pico-sized haptophytes in the surface water and subsurface chlorophyll maximum (SCM) layer of the AO were explored. Our results found higher alpha diversity estimates in the surface water compared with in the SCM based on high-throughput sequencing of haptophyte specific 18S rRNA. The community composition of the surface water was significantly different from that of the SCM, and water temperature was identified as the primary factor shaping the community compositions. Prymnesiales (mostly Chrysochromulina), uncultured Prymnesiophyceae, and Phaeocystis dominated the surface water communities, whereas Phaeocystis dominated the SCM communities, followed by Chrysochromulina, uncultured Prymnesiophyceae, and the remaining taxa. The communities of the surface water and SCM layer developed relatively independent modules in the metacommunity network. Nodes in the surface water were more closely connected to one another than those in the SCM. Network stability analysis revealed that surface water networks were more stable than SCM networks. These findings suggest that SCM communities are more susceptible to environmental fluctuations than those in surface water and that future global changes (e.g., global warming) may profoundly influence the development, persistence, and service of SCM in the AO.