Water masses and their associated temperature and cross-domain biotic factors co-shape upwelling microbial communities

An unignorable human health risk posed by antibiotic resistome and microbiome in urban rivers: Insights from Beijing, China

Urban rivers are the main water bodies humans frequently come into contact with, so the risks posed are closely monitored. Antibiotic resistance genes (ARGs) residues in reclaimed water pose serious risks to human health. There are urgent needs to improve the understanding of distribution of and risks posed by ARGs in urban rivers. In this study, shotgun metagenomic approach was used to characterize ARGs, mobile genetic elements (MGEs), and virulence factors (VFs) in water and sediment from Xinfeng River in Beijing and to identify microbes, potential antibiotic resistant bacteria, and human pathogens (HPs). MGE, microbial community, VF, and ARG co-occurrences were used to assess the environmental risks posed by ARGs. The results indicated that quinolone was the most abundant ARG type and that tufA and fusA were the two dominant ARG subtypes. Wetland effluent increased ARG abundance in the river, and the effect was detected even 50 m downstream. ARG abundances and distribution in the river had difference in different seasons. The dominant bacteria in the river were Proteobacteria, Bacteroidetes, and Actinobacteria, and 59 HPs were detected. In total, 69 MGEs and 19 VFs were found. Co-occurrence networks indicated that potential antibiotic resistant bacteria, MGEs, VFs, and ARGs in the river significantly correlated, indicating the potential risks posed by ARGs. The results improve our understanding of ARG distribution and environmental risks in urban river water. More attention should be paid to controlling environmental risks posed by ARGs in urban river and reclaimed water.

Disentangling the role of biotic and abiotic factors in shaping microbial biogeographic patterns in the large spatial scale ocean

Across large oceanic spatial scales, the mechanisms governing microbial community composition remain poorly understood, particularly regarding the influence of biotic factors. In this study, samples were collected over a 3772.55-km transect in the Western Pacific surface waters, with 16S and 18S ribosomal RNA gene sequencing employed to elucidate the assembly processes and drivers of microbial communities. Both eukaryotic (EM) and prokaryotic microbial (PM) communities exhibited a significant distance-decay relationship throughout the survey, although their assembly mechanisms differed. Biotic factors played a pivotal role in shaping the biogeographical patterns of both EM and PM communities, indicating that cross-domain microbial interactions significantly affect each other's distribution. Furthermore, the stability of the co-occurrence network was more strongly influenced by cross-domain microbial relationships than by keystone taxa. These findings advance our understanding of the mechanisms and biogeographical patterns sustaining microbial diversity across large oceanic spatial scales.

Microbial responses on upwelling signature across reversal wind pattern in tropical coastal environment off Mumbai, India

Upwelling promotes marine productivity through water column mixing. The process disturbs the ecosystem, causing oxygen depletion and thermal variability. This study analyses effect of upwelling processes on microbial signature in coastal waters off Mumbai. The coastal environment with seasonal reversal winds was analysed using data during ten cruises. Coastal metocean processes are examined using water quality parameters and the Ekman approximation with wind stress. This analysis explains oxygen depletion and coastal upwelling, influenced by seasonal reversal wind pattern. The study connects hypoxia in the coastal water column to wind-induced upwelling. Concurrently, microbial structure is assessed through metrics such as Total Viable Count, Total Bacterial Count, Sulfate Reducing Bacteria (SRB), and denitrifiers. Notably, high levels of SRB are observed during hypoxia associated with coastal upwelling. This study investigates microbial level with combined result of physical processes and water quality parameters.

Flooding promotes the coalescence of microbial community in estuarine habitats

Microbial community coalescence describes the mixing of microbial communities and their integration with the surrounding environment, which is common in natural ecosystems and has potential impacts on ecological processes. However, few studies have focused on microbial community coalescence between different habitats in estuarine regions. In this study, we comprehensively investigated the environmental characteristics and bacterial community changes of different habitats (water body (Water), subtidal sediments (SS) and intertidal salt marsh sediments (SM)) in Luanhe estuary during flood and normal flow periods. The results showed that flood event significantly reduced the salinity of the estuarine habitats, changed the nutrient structure and intensified the eutrophication of estuarine water. By calculating the proportion of overlapping groups and applying the 'FEAST' algorithm, we revealed that flood event facilitated the migration of bacterial communities along alternative pathways across habitats, markedly enhanced the cross-habitat mobility of bacterial communities, which underscores the pivotal role of flood event in driving bacterial community coalescence. Flood-induced community coalescence not only increased the α-diversity of bacterial communities within habitats, but also increased the proportion of overlapped species between habitats, ultimately leading to homogenization between habitats. Canonical correlation analysis combined co-occurrence network analysis revealed that flood event attenuated the role of environmental filtration in microbial assembly, while increased the impact of dispersal processes and intensified interspecific competition among microorganisms, led to the change of keystone species and reduced the complexity and stability of bacterial communities. In conclusion, this study demonstrates the complex effects of flood events on estuarine microbial communities from the perspective of multi-habitat interactions in the estuary, and emphasizes the key role of river hydrodynamic conditions in facilitating the coalescence of estuarine microbial communities. We look forward to further attention and research on estuarine microbial coalescence, which will provide new insights into assessing the stability and resilience of estuarine ecosystems under flood challenges and the sustainable management of estuarine wetlands.

Microbial distribution in Mudbank regions off Alappuzha, South-West coast of India

The coastal waters of Kerala, in the South Eastern Arabian Sea (SEAS), are unique during the Southwest monsoon season due to the concurrent occurrence of two physical processes, the upwelling and Mudbanks. However, little is known about the viral ecology and activity in a system where upwelling and Mudbanks coexist, though it is generally recognized that microbial assemblages play a vital role in the food web dynamics of marine systems, particularly in upwelling. Water samples were taken from monsoon and pre-monsoon period from three locations, (M1, M2, and M3) off Alappuzha, on the southwest coast of India to examine the viral activity and distribution. The dissolved oxygen levels showed the incursion of hypoxic waters in all the stations during the peak upwelling period. Upwelling signals were prominent in all the stations, but Mudbank and upwelling co-occurred at M2 alone during monsoon. The abundance of viruses ranged from 0.86 to 15.68 × 106 Viral-like Particles (VLPs mL-1) and prokaryotic abundance ranged from 2.73-16.26 × 105 cells mL-1. Viral and prokaryotic abundance was significantly higher in the monsoon compared to pre and late-monsoon. Based on Transmission electron microscopy (TEM) results, the non-tailed viruses constituted the major (43%) proportion of the total viruses in this study region. However, the viral production rates and viral-mediated bacterial mortality were high in the pre-monsoon compared to the monsoon and late-monsoon periods. There was no obvious effect of Mudbanks on viral dynamics and the observed variations in virological and hydrological features were governed mainly by coastal upwelling.

Biogeography and dynamics of prokaryotic and microeukaryotic community assembly across 2600 km in the coastal and shelf ecosystems of the China Seas

Marine prokaryotes and microeukaryotes are essential components of microbial food webs, and drive the biogeochemical cycling. However, the underlying ecological mechanisms driving prokaryotic and microeukaryotic community assembly in large-scale coastal ecosystems remain unclear. In this study, we studied biogeographic patterns of prokaryotic and microeukaryotic communities in the coastal and shelf ecosystem of the China Seas. Results showed that prokaryotic richness was the highest in the Yangtze River Plume, whereas microeukaryotic richness decreased from south to north. Prokaryotic-microeukaryotic co-occurrence networks display greater complexity in the Yangtze River Plume compared to other regions, potentially indicating higher environmental heterogeneity. Furthermore, the cross-domain networks revealed that prokaryotes were more interconnected with each other than with microeukaryotes or between microeukaryotes, and all hub nodes were bacterial taxa, suggesting that prokaryotes may be more important for sustaining the stability and multifunctionality of coastal ecosystem than microeukaryotes. Variation Partitioning Analysis revealed that approximately equal proportions of environmental, biotic and spatial factors contribute to variations in microbial community composition. Temperature was the primary environmental driver of both prokaryotic and microeukaryotic communities across the China Seas. Additionally, stochastic processes (dispersal limitation) and deterministic processes (homogeneous selection) were two major ecological factors in shaping microeukaryotic and prokaryotic assemblages, respectively, suggesting their different environmental plasticity and evolutionary mechanisms. Overall, these results demonstrate both prokaryotic and microeukaryotic communities displayed a latitude-driven distribution pattern and different assembly mechanisms, improving our understanding of microbial biogeography patterns under global change and anthropogenic activity driven habitat diversification in the coastal and shelf ecosystem.

Identification of keystone taxa in rhizosphere microbial communities using different methods and their effects on compounds of the host Cinnamomum migao

Exploring keystone taxa affecting microbial community stability and host function is crucial for understanding ecosystem functions. However, identifying keystone taxa from humongous microbial communities remains challenging. We collected 344 rhizosphere and bulk soil samples from the endangered plant C. migao for 2 years consecutively. Used high-throughput sequencing 16S rDNA and ITS to obtain the composition of bacterial and fungal communities. We explored keystone taxa and the applicability and limitations of five methods (SPEC-OCCU, Zi-Pi, Subnetwork, Betweenness, and Module), as well as the impact of microbial community domain, time series, and rhizosphere boundary on the identification of keystone taxa in the communities. Our results showed that the five methods, identified abundant keystone taxa in rhizosphere and bulk soil microbial communities. However, the keystone taxa shared by the rhizosphere and bulk soil microbial communities over time decreased rapidly decrease in the five methods. Among five methods on the identification of keystone taxa in the rhizosphere community, Module identified 113 taxa, SPEC-OCCU identified 17 taxa, Betweenness identified 3 taxa, Subnetwork identified 3 taxa, and Zi-Pi identified 4 taxa. The keystone taxa are mainly conditionally rare taxa, and their ecological functions include chemoheterotrophy, aerobic chemoheterotrophy, nitrate reduction, and anaerobic photoautotrophy. The results of the random forest model and structural equation model predict that keystone taxa Mortierella and Ellin6513 may have an effects on the accumulation of 1, 4, 7, - Cycloundecatriene, 1, 5, 9, 9-tetramethyl-, Z, Z, Z-, beta-copaene, bicyclogermacrene, 1,8-Cineole in C. migao fruits, but their effects still need further evidence. Our study evidence an unstable microbial community in the bulk soil, and the definition of microbial boundary and ecologically functional affected the identification of keystone taxa in the community. Subnetwork and Module are more in line with the definition of keystone taxa in microbial ecosystems in terms of maintaining community stability and hosting function.

Biogeographical and biodiversity patterns of planktonic microeukaryotes along the tropical western to eastern Pacific Ocean transect revealed by metabarcoding

Microeukaryotic plankton (0.2-200 µm), which are morphologically and genetically highly diverse, play a crucial role in ocean productivity and carbon consumption. The Pacific Ocean (PO), one of the world's largest oligotrophic regions, remains largely unexplored in terms of the biogeography and biodiversity of microeukaryotes based on large-scale sampling. We investigated the horizontal distribution of microeukaryotes along a 16,000 km transect from the west to the east of the PO. The alpha diversity indices showed a distinct decreasing trend from west to east, which was highly correlated with water temperature. The microeukaryotic community, which was clustered into the western, central, and eastern PO groups, displayed a significant distance-decay relationship. Syndiniales, a lineage of parasitic dinoflagellates, was ubiquitously distributed along the transect and dominated the community in terms of both sequence and zero-radius operational taxonomic unit (ZOTU) proportions. The prevailing dominance of Syndiniales-affiliated ZOTUs and their close associations with dinoflagellates, diatoms, and radiolarians, as revealed by SparCC correlation analysis, suggested that parasitism may be an important trophic strategy in the surface waters of the PO. Geographical distance and temperature were the most important environmental factors that significantly correlated with community structure. Overall, our study sheds more light on the distribution pattern of both alpha and beta diversities of microeukaryotic communities and highlighted the importance of parasitisms by Syndiniales across the tropical PO.IMPORTANCEUnderstanding the biogeographical and biodiversity patterns of microeukaryotic communities is essential to comprehending their roles in biogeochemical cycling. In this study, planktonic microeukaryotes were collected along a west-to-east Pacific Ocean transect (ca. 16,000 km). Our study revealed that the alpha diversity indices were highly correlated with water temperature, and the microeukaryotic communities displayed a distinct geographical distance-driven pattern. The predominance of the parasitic dinoflagellate lineage Syndiniales and their close relationship with other microeukaryotic groups suggest that parasitism may be a crucial survival strategy for microeukaryotes in the surface waters of the Pacific Ocean. Our findings expand our understanding of the biodiversity and biogeographical pattern of microeukaryotes and highlight the significance of parasitic Syndiniales in the surface ocean.

Climate and local environment co-mediate the taxonomic and functional diversity of bacteria and archaea in the Qinghai-Tibet Plateau rivers

Understanding the environmental response patterns of riverine microbiota is essential for predicting the potential impact of future environmental change on river ecosystems. Vulnerable plateau ecosystems are particularly sensitive to climate and local environmental changes, however, the environmental response patterns of the taxonomic and functional diversity of riverine microbiota remain unclear. Here, we conducted a systematic investigation of the taxonomic and functional diversity of bacteria and archaea from riparian soils, sediments, and water across the elevation of 1800- 4800 m in the Qinghai-Tibet Plateau rivers. We found that within the elevation range of 1800 to 3800 m, riparian soils and sediments exhibited similarities and stabilities in microbial taxonomic and functional diversity, and water microbiomes were more sensitive with great fluctuations in microbial diversity. Beyond the elevation of 3800 m, microbial diversity declined across all riverine matrixes. Local environmental conditions can influence the sensitivity of microbiomes to climate change. The combination of critical climate and local environmental factors, including total nitrogen, total organic carbon, as well as climate variables associated with temperature and precipitation, provided better explanations for microbial diversity than single-factor analyses. Under the extremely adverse scenario of high greenhouse gas emission concentrations (SSP585), we anticipate that by the end of this century, the bacterial, archaeal, and microbial functional diversity across the river network of the Yangtze and Yellow source basin would potentially change by -16.9- 5.2 %, -16.1- 5.7 %, and -9.3- 6.4 %, respectively. Overall, climate and local environments jointly shaped the microbial diversity in plateau river ecosystems, and water microbiomes would provide early signs of environmental changes. Our study provides effective theoretical foundations for the conservation of river biodiversity and functional stability under environmental changes.

Microbial interactions strengthen deterministic processes during community assembly in a subtropical estuary

Systematic studies on the assembly process and driving mechanisms of microbial communities in estuaries with diverse seasonal and spatial scales are still limited. In this study, high-throughput sequencing, and microbial network analysis were combined to decipher the impact of environmental changes and biological interactions on the maintenance of microbial diversity patterns in the Jiulong River Estuary (JRE). The results showed that overall, stochastic processes dominated the bacterioplankton community assembly in the estuary, accounting for 49.66-74.78 % of the total. Additionally, bacterioplankton community diversity varied significantly across seasons and subzones. Specifically, the concentration of soluble reactive phosphorus (SRP) in the estuary steadily reduced from winter to summer, and the corresponding bacterioplankton community interactions gradually shifted from the weakest interaction in winter to the strongest in summer. The deterministic processes contributed more than half (50.34 %) to microbial assembly in the summer, but only 25.22 % in winter. Deterministic processes prevailed in the seaward with low SRP concentrations and strong bacterioplankton community interactions, while stochastic processes contributed 70.14 % to the assembly of microbial communities riverward. Biotic and abiotic factors, such as nutrients and microbial interactions, jointly drove the seasonal and spatial patterns of bacterioplankton community assembly, but overall, nutrients played a dominant role. Nevertheless, the contributions of nutrients and microbial interactions were equivalent in spatial assembly processes, albeit nutrients were the primary seasonal driver of the bacterioplankton community assembly process. This study emphasizes the significance of microbial interactions in the bacterioplankton community assemblage. These findings provide new and comprehensive insights into the microbial communities' organization in estuaries.

Selective cell lysis pressure on rare and abundant prokaryotic taxa across a shelf-to-slope continuum in the Northern South China Sea

IMPORTANCE

Virus-induced host lysis contributes up to 40% of total prokaryotic mortality and plays crucial roles in shaping microbial composition and diversity in the ocean. Nonetheless, what taxon-specific cell lysis is caused by viruses remains to be studied. The present study, therefore, examined the taxon-specific cell lysis and estimated its contribution to the variations in the rare and abundant microbial taxa. The results demonstrate that taxon-specific mortality differed in surface and bottom of the coastal environment. In addition, active rare taxa are more susceptible to heightened lytic pressure and suggested the importance of viral lysis in regulating the microbial community composition. These results improve our understanding of bottom-up (abiotic environmental variables) and top-down (viral lysis) controls contributing to microbial community assembly in the ocean.

Global biogeography of the smallest plankton across ocean depths

Tiny ocean plankton (picoplankton) are fundamental for the functioning of the biosphere, but the ecological mechanisms shaping their biogeography were partially understood. Comprehending whether these microorganisms are structured by niche versus neutral processes is relevant in the context of global change. We investigate the ecological processes (selection, dispersal, and drift) structuring global-ocean picoplanktonic communities inhabiting the epipelagic (0 to 200 meters), mesopelagic (200 to 1000 meters), and bathypelagic (1000 to 4000 meters) zones. We found that selection decreased, while dispersal limitation increased with depth, possibly due to differences in habitat heterogeneity and dispersal barriers such as water masses and bottom topography. Picoplankton β-diversity positively correlated with environmental heterogeneity and water mass variability, but this relationship tended to be weaker for eukaryotes than for prokaryotes. Community patterns were more pronounced in the Mediterranean Sea, probably because of its cross-basin environmental heterogeneity and deep-water isolation. We conclude that different combinations of ecological mechanisms shape the biogeography of the ocean microbiome across depths.

Community structures of mangrove endophytic and rhizosphere bacteria in Zhangjiangkou National Mangrove Nature Reserve

Bacterial communities play an important role in mangrove ecosystems. In order to gain information on the bacterial communities in mangrove species and rhizospheres grown in Zhangjiangkou National Mangrove Nature Reserve, this study collected root, branch, and leaf samples from five mangrove species as well as rhizosphere and non-rhizosphere samples and analyzed the community structure of endophytic bacteria and bacteria in rhizosphere and non-rhizosphere using Illumina high-throughput sequencing technique. Bacteria in 52 phyla, 64 classes, 152 orders, 295 families, and 794 genera were identified, which mainly belonged to Proteobacteria, Cyanobacteria, Actinobacteria, Firmicutes, Bacteroidetes, Fusobacteria, and Nitrospirota. At each taxonomic level, the community structure of the rhizosphere bacteria varied slightly with mangrove species, but endophytic bacteria differed greatly with plant species. The diversity indices of endophytic bacteria in branch and leaf samples of Acanthus ilicifolius were significantly lower, and endophytic bacteria in the plant tissues had higher abundance in the replication/repair and translation Clusters of Orthologous Genes functional categories but lower abundance in the carbohydrate metabolism category. This study helps to understand the community structure and diversity characteristics of endophytic and rhizosphere bacteria in different mangrove plants. Provide a theoretical basis for in-depth research on the functions of mangrove ecosystems.

Dissemination of antibiotic resistance genes from the Pearl River Estuary to adjacent coastal areas

The spread of antibiotic resistance genes (ARGs) is a growing concern over the world's various environments. Coastal environments may receive pollutants from land runoffs via estuaries. However, the impact of ARG contamination from estuarine regions to coastal areas is rarely reported. This study used high-throughput quantitative PCR to examine the diversity and abundance of ARGs in Pearl River Estuary (PRE) and adjacent coastal areas. We found that the distribution of ARGs in seawater exhibited the distance-decay phenomenon from the estuary to coastal areas, while the sediment samples did not exhibit an obvious distribution pattern. The estuarine water was found to be the hotspot of ARGs, with 74 ARG species detected and absolute abundance being 5.93 × 10⁵ copies per mL, on average, while less species and lower abundance of ARGs were detected in coastal waters. Ordination analysis showed that estuarine ARG communities were significantly different from coastal ARG communities for water samples. SourceTracker analysis revealed that ARGs from the estuarine environment contributed only a minor fraction of ARG contamination to downstream coastal areas (1.5%-7.4% for water samples, and 0.7-1.8% for sediment samples), indicating the strong dilution effect of seawater. Mantel tests, redundancy analysis and random forest model analysis identified salinity, nutrients, microbial community structure and mobile genetic elements (MGEs) as important factors influencing ARG distribution. Partial least squares-path model revealed that, among all environmental factors, MGEs directly affected the distribution of ARGs, while other factors indirectly contributed by affecting the MGEs assemblage. Our study provides insight into the dissemination of ARGs from the PRE to adjacent coastal areas.

From the Sunlit to the Aphotic Zone: Assembly Mechanisms and Co-Occurrence Patterns of Protistan-Bacterial Microbiotas in the Western Pacific Ocean

We know little about the assembly processes and association patterns of microbial communities below the photic zone. In marine pelagic systems, there are insufficient observational data regarding why and how the microbial assemblies and associations vary from photic to aphotic zones. In this study, we investigated size-fractionated oceanic microbiotas, specifically free-living (FL; 0.22 to 3 μm) and particle-associated (PA; >3 μm) bacteria and protists (0.22 to 200 μm) collected from the surface to 2,000 m in the western Pacific Ocean, to see how assembly mechanisms and association patterns changed from photic to aphotic zones. Taxonomic analysis revealed a distinct community composition between photic and aphotic zones that was largely driven by biotic associations rather than abiotic factors. Aphotic community co-occurrence was less widespread and robust than its photic counterparts, and biotic associations were crucial in microbial co-occurrence, having a higher influence on photic than aphotic co-occurrences. The decrease in biotic associations and the increase in dispersal limitation from the photic to the aphotic zone affect the deterministic-stochastic balance, leading to a more stochastic-process-driven community assembly for all three microbial groups in the aphotic zone. Our findings significantly contribute to our understanding of how and why microbial assembly and co-occurrence vary from photic to aphotic zones, offering insight into the dynamics of the protistan-bacterial microbiota in the western Pacific's photic and aphotic zones. IMPORTANCE We know little about the assembly processes and association patterns of microbial communities below the photic zone in marine pelagic systems. We discovered that community assembly processes differed between photic and aphotic zones, with all three microbial groups studied (protists and FL and PA bacteria) being more influenced by stochastic processes than in the photic zone. The decrease in organismic associations and the increase in dispersal limitation from the photic to the aphotic zone both have an impact on the deterministic-stochastic balance, resulting in a more stochastic process-driven community assembly for all three microbial groups in the aphotic zone. Our findings significantly contribute to the understanding of how and why microbial assembly and co-occurrence change between photic and aphotic zones, offering insight into the dynamics of the protist-bacteria microbiota in the western Pacific oceans.

Freshwater trophic status mediates microbial community assembly and interdomain network complexity

Freshwater microorganisms and their interactions are important drivers of nutrient cycling that are in turn affected by nutrient status, causing shifts in microbial community diversity, composition, and interactions. However, the impact of water trophic status on bacterial-archaeal interdomain interactions remains poorly understood. This study focused on the impact of trophic status, as characterized by trophic state index (TSI), on the interdomain interactions of freshwater microbial communities from 45 ponds in Hangzhou. Our results showed that the mesotrophic wetland bordering on lightly eutrophic (Hemu: TSI of 49; lightly eutrophic is defined as 50 ≤ TSI <60) harbored a much more complex bacterial-archaeal interdomain network, which showed significantly (P < 0.05) higher connectivity than the wetlands with lower (TSI of 38) or higher (TSI of 57) trophic levels. Notably, light eutrophication strengthened the network modules' negative associations with organic carbon through some network hubs, which could trigger carbon loss in wetlands. We also detected a non-linear response of interdomain network complexity to the increasing of nutrients with a turning point of approximately TSI 50. Quantitative estimates of community assembly processes and structural equation modelling analysis indicated that chlorophyll-a, total nitrogen, and total phosphorus could regulate interdomain network complexity (50% of the variation explanation rate) by driving microbial community assembly. This study demonstrates that microbial interdomain network complexity could be used as a bioindicator for ecological changes, which would helpful for improving ecological assessment of the freshwater eutrophication.

Different Responses of Bacteria and Microeukaryote to Assembly Processes and Co-occurrence Pattern in the Coastal Upwelling

Upwelling may generate unique hydrological and environmental heterogeneity, leading to enhanced diffusion to reshape microbial communities. However, it remains largely unknown how different microbial taxa respond to highly complex and dynamic upwelling systems. In the present study, geographic patterns and co-occurrence network of different microbial communities in response to upwelling were examined. Our results showed that coastal upwelling shaped prokaryotic and eukaryotic microbial community and decreased their diversity. In addition, bacteria and microeukaryote had similar biogeographical patterns with distinct assembly mechanisms. The impact of stochastic processes on bacteria was significantly stronger compared with microeukaryote in upwelling. Lower network complexity but more frequent interaction was found in upwelling microbial co-occurrence. However, the upwelling environment increased the robustness and modularity of bacterial network, while eukaryotic network was just the opposite. Co-occurrence networks of bacteria and microeukaryote showed significant distance-decay patterns, while the bacterial network had a stronger spatial variation. Temperature and salinity were the strongest environmental factors affecting microbial coexistence, whereas the topological characteristics of bacterial and eukaryotic networks had different responses to the upwelling environment. These findings expanded our understanding of biogeographic patterns of microbial community and ecological network and the underlying mechanisms of different microbial taxa in upwelling.