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

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. ABSTRACT 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.