- 1Key Laboratory of Marine Ecosystem Dynamics, Second Institute of Oceanography, Ministry of Natural Resources, Hangzhou, China
- 2Ningbo Institute of Oceanography, Ningbo, China
- 3Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Zhuhai, China
- 4College of Oceanography, Hohai University, Nanjing, China
- 5School of Oceanography, Shanghai Jiao Tong University, Shanghai, China
Microbial food web (MFW) dominates the energy flow in oligotrophic tropical open ocean pelagic ecosystems. Understanding biogeographic patterns and driving mechanisms of key components of the MFW is one of the central topics in current marine ecology. Investigations were conducted along an 1,100-km horizontal gradient and in the full-water column vertical gradient of the oligotrophic tropical western Pacific Ocean. High-throughput sequencing and association networking methods were used to analyze the community structure and interspecies interactions of MFW. The structure of MFW significantly differed with depths, but not across horizontal gradients. Bacteria and microeukaryotes were interconnected and had more predominantly positive and negative linkages in the aphotic layers. Key components of MFW exhibited similar biogeographic patterns and driving mechanisms. Geographic distance exerted minimal effects on the distribution patterns of the microbial food web, while environmental factors played more important roles, especially for temperature and inorganic nutrients. Stochastic processes were more important in the microbial food webs of the 5–200 m layer than the >500 m layer, and drift explained the majority of stochastic processes. Moreover, only a weak but not significant driving force for North Equatorial Current on the east–west connectivity of the microbial food web was found in the upper layers. This knowledge is a critical fundamental data for future planning of marine protected areas targeting the protection of tuna fishing in the western Pacific Ocean.
1. Introduction
The significance of microbial food web (MFW) in oligotrophic ocean was originally proposed by Pomeroy (1974). Azam et al. (1983) formalized the concept of the microbial loop; then, Sherr and Sherr (1988) broadened the definition of the microbial loop to include direct trophic links between phytoplankton and phagotrophic protozoa, in addition to the multi-step detritus-bacteria flagellate-ciliate pathway. Currently, MFW mainly consists of bacteria, archaea, Synechococcus, Prochlorococcus, picoeukaryotes, microzooplankton (heterotrophic and mixotrophic nanoflagellates and ciliates) and planktonic viruses, which are key groups that control nutrient cycling and biomass production in aquatic ecosystems, especially oligotrophic oceans (Seymour et al., 2017). In the MFW, different species are directly and indirectly linked through biogeochemical cycles and food web interactions. Heterotrophic bacteria and phytoplankton are interconnected by a two-way flux: phytoplankton excretion or cell lysis is a source of organic matter for heterotrophic bacteria and the mineralization of this organic matter by heterotrophic bacteria in turn provides nutrients to primary producers (Sarmento et al., 2010). Protists (e.g., nanoflagellates, ciliates) are grazers in the MFW, which are strictly heterotrophic or mixotrophic, and incorporate carbon by ingesting small primary producers and/or heterotrophic bacteria (Zubkov and Tarran, 2008). A correlation network analysis revealed that sequencing data can be used to hypothesize positive (mutualism and symbiosis) and negative (competition and parasitism) ecological interactions between various species (Faust and Raes, 2012). To date, only a few studies have investigated the MFW in oligotrophic oceans (Calbet and Landry, 1999).
The tropical western Pacific Ocean is one of the world’s largest warm oligotrophic regions, which has been largely unexplored with regard to the diversity and biogeography of planktonic microbiome, based on Tara Oceans data (Pesant et al., 2015). Only a few studies have used micro-optical based techniques, such as microscopy and flow cytometry, to study the community structure and biomass of phytoplankton (Dai et al., 2020; Wei et al., 2020), picoplankton, and virioplankton in the surface layer and vertical changes in the MFW around seamounts in the oligotrophic tropical western Pacific Ocean (Liang et al., 2017; Zhao et al., 2020; Wei et al., 2021). The development of high-throughput sequencing technology over the last decade has greatly improved our understanding of planktonic prokaryotes and eukaryotes in the oligotrophic tropical western Pacific Ocean (Wang et al., 2020; Wu et al., 2020; Zhao et al., 2020; Shao et al., 2022). Compared to optical based methods, it can be used to identify variations in less abundant species and picoplankton populations of the MFW. The tropical western Pacific Ocean is important tuna fishing ground (Worm et al., 2005), but primary production cannot be passed up to these large predators without the transfer of MFW (Calbet, 2008; Zeldis and Décima, 2020). Knowledge about the geographical patterns and driving mechanisms of MFW is therefore a critical fundamental data for future planning of marine protected areas targeting the protection of these large predators. In addition, we rely heavily on this knowledge to assess the impacts of climate change on pelagic ecosystems in open ocean (Pinsky et al., 2020).
One of central topics in current marine ecology is to elucidate the biogeographic patterns and driving mechanisms of marine microorganisms. Knowledge of the underlying mechanisms that shape microbial biogeographic patterns is necessary for understanding the community structure and to exercise best practices in the stewardship of marine resources (Zhou and Ning, 2017). However, most studies have focused on eutrophic and mesotrophic offshore and coastal environments (Djurhuus et al., 2020; Richter et al., 2020); only a few have investigated oligotrophic oceans. Previous studies demonstrated that environmental factors (temperature, light, and nutrients) and spatial factors significantly correlated with the distribution of bacteria and microeukaryotes in the oligotrophic tropical western Pacific Ocean (Wang et al., 2020; Wu et al., 2020; Zhao et al., 2020; Shao et al., 2022), indicating that environmental selection-related processes are more important than dispersal-related processes in shaping the distribution of bacteria and microeukaryotes. However, the proportion of bacterial and microeukaryotic community variations explained by both processes is small, which implies that the complex distribution patterns within communities and other drivers (currents, water mass, and biotic interactions) may significantly contribute to shaping the distribution of MFW in the oligotrophic tropical western Pacific Ocean. Null model analyses can be used to assess the importance of different drivers of bacterial and microeukaryotic community assembly in the oligotrophic tropical western Pacific Ocean (Kong et al., 2022; Xu et al., 2022), including dispersal limitation, homogenizing dispersal, heterogeneous selection, homogeneous selection, and drift (Stegen et al., 2013). The strong Pacific North Equatorial Current (NEC) is the main influencing current of tropical western Pacific Ocean, providing favorable conditions for east–west connectivity in the upper (0–200 m) western Pacific between 8 and 18°N (Qiu and Chen, 2012; Hsu et al., 2017). However, the importance of the NEC in shaping the MFW in the oligotrophic tropical western Pacific Ocean remains poorly understood.
In this study, the spatial variation and driving factors of bacterial and microeukaryotic communities in the oligotrophic tropical western Pacific Ocean were analyzed to enhance our understanding of the biogeographic patterns and assembly mechanisms of the MFW. In this extremely oligotrophic environment, we hypothesized that: (i) bacteria and microeukaryotes have different predominant linkages with depths; (ii) bacterial and microeukaryotic communities share similar biogeographic patterns and assembly mechanisms; and (iii) the NEC enhances the east–west connectivity of the MFW in the upper tropical western Pacific Ocean.
2. Materials and methods
2.1. Sampling site and collection
Sampling sites and sample collection were processed and analyzed following previously described methods (Shao et al., 2022), which were consistent between microeukaryotes and bacteria (Figure 1). Depths of eight stations ranged from 1,629 to 5,648 m, and we collected samples from six different depths (5 m, 75 m, deep chlorophyll maximum (DCM), 200 m, 500 m, and 3,000 m) at each station. The annual mean chlorophyll a concentration (mg/m3) was downloaded from the E.U. Copernicus Marine Service, which indicated that the sampling area was a typical oligotrophic sea. Eight sampling sites, except S1, were located in the influence area of the Pacific NEC (130–160°E, 8–18°N).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1098264/fmicb-14-1098264-HTML/image_m/fmicb-14-1098264-g001.jpg)
Figure 1. Location of eight water column sampling sites (black circles) and schematic illustration of the influence range of North Equatorial Current (NEC, yellow arrow) in the western tropical Pacific Ocean. Yellow boxes indicated two sections in the NEC influence area. The background was the annual mean chlorophyll a concentration (mg/m3).
2.2. DNA extraction, PCR amplification, and Illumina sequencing
DNA extraction was conducted using a DNeasy Power Soil Kit (QIAGEN, Valencia, CA, United States) according to the manufactures instructions. The concentration and purity of the total DNA were measured by a NanoDrop ND-100 spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, United States). Primers 547F (5’-CCAGCASCYGCGGTAATTCC-3′) and V4R (5’-ACTTTCGTTCTTGATYRA-3′) were used to amplify the V4 region of 18S rRNA (Stoeck et al., 2010). Adapter sequences were linked to the barcodes at the 5′ end of each primer. Microeukaryote PCR amplification, purification, and quantification were conducted following the same methods as for bacteria (Shao et al., 2022). High-throughput gene sequencing was performed on an Illumina Novaseq-PE250 platform (paired-end reads, 2 × 300 bp) by Personal Biotechnology, Co., Ltd., (Shanghai, China).
2.3. Sequence assembly, clustering, and annotation
The quality filter and assembly of raw microeukaryote sequence were the same as for bacteria, which followed previously described methods (Shao et al., 2022). High-quality representative sequences for amplicon sequence variants (ASVs) were assigned with 100% sequence identity using UCLUST (Edgar, 2010). Annotation was carried out using UCLUST and the National Center for Biotechnology Information (NCBI) gene database. The ASVs affiliated with chloroplasts, mitochondria, unclassified groups, and singletons were removed from the dataset. To correct possible errors induced by unequal sequencing efforts in all samples, the ASV table was randomly subsampled to ensure an equal number of sequences per sample (16,704) using MOTHUR v1.33.3 (Schloss et al., 2009). The α-diversity indices, including Pielou’s evenness, Faith’s phylogenetic diversity, and Shannon-Wiener and Simpson indices, were calculated using QIIME v2.0.0 (Caporaso et al., 2010). Good’s coverage was calculated using MOTHUR v1.33.3 (Schloss et al., 2009).
2.4. Statistical analyses
All statistical analyses were performed using R v4.1.1 software,1 unless noted otherwise. Prior to statistical analysis, the ASV table was Hellinger-transformed and the environmental variables were normalized using the “vegan” R package to improve normality and homoscedasticity (Oksanen et al., 2007).
The group compositions of the microeukaryotic communities with relative abundances >1% were visualized at the phylum and order levels. Bacterial taxonomic groups with relative abundances >1% were visualized at the order level. The α-diversity indices were visualized using the “ggplot2” R package (Gómez-Rubio, 2017). Wilcox tests were used to evaluate differences in the α-diversity indices between the six depths using the “ggpubr” R package (Kassambara, 2020). Venn diagrams were used to show the microeukaryotic relationships among the six depths using the “VennDiagram” R package (Chen and Boutros, 2011). The β-diversity was calculated using Bray-Curtis distances and visualized with a principal coordinated analysis (PCoA) ordination. Linear discriminant analysis effect size (LEfSe) was used to screen potential statistically significant microeukaryotic taxa with total abundances >10 among the six depths. The non-parametric factorial Kruskal-Wallis sum-rank test (alpha value = 0.05) was used to detect taxa with significantly different abundances between groups. A linear discriminant analysis (threshold score = 4.0) was used to estimate the effect sizes of each abundant taxa (Segata et al., 2011).
A network analysis was calculated and visualized using Cytoscape v3.9.1. First, the orders of bacteria and microeukaryotes with average relative abundance of less than 1% were removed, and a Spearman rank correlation matrix was created using the remaining data. Then, co-occurrence patterns were determined if the Spearman’s correlation coefficient was >0.4 (positive) or < −0.4 (negative) and p < 0.05 (Zhou et al., 2018). To test the effects of the NEC on the MFW, subsequent analyses were divided into 5–200 m and > 500 m layers of the bacterial and microeukaryotic communities. To study distance-decay relationships in the MFW, the relationship between the Bray-Curtis dissimilarities of bacterial and microeukaryotic communities and geographic distances was analyzed based on Spearman’s rank correlations. A permutational multivariate analysis of variation (PERMANOVA) was applied to partition the variations of all, 5–200 m, and > 500 m layers of bacterial and microeukaryotic community compositions among different depths, stations, and sections. Sections were classified by different latitudes in the NEC influence area (Figure 1; Supplementary Figures S2–S4 were one section, and Supplementary Figures S5–S8 were another section). A set of spatial factors was generated using the procedure coordinates of neighbor matrices (PCNMs) method based on the longitudes and latitudes of the sampling stations (Borcard and Legendre, 2002). Mantel tests were conducted to investigate the sources of bacterial and microeukaryotic variation while considering both spatial and environmental factors (Diniz-Filho et al., 2013). Redundancy analysis and Spearman’s rank correlation were used to explore the effects of environmental variables on the MFW (Oksanen et al., 2007). Then, we used mantel tests and variation partitioning analysis to quantify the relative effects of environmental and spatial factors in shaping the MFW following previously described methods (Shao et al., 2022).
To evaluate the relative effects of stochastic and deterministic processes, phylogenetic structures of the bacterial and microeukaryotic communities were used to obtain insights into the drivers of community assembly. We assessed the relative importance of heterogeneous selection, dispersal limitation, drift, homogenizing dispersal, and homogeneous selection using a three-step framework (Stegen et al., 2013). First, the intercommunity phylogenetic turnover was calculated using the mean nearest taxon distance (βMNTD). Then, the β-nearest taxon index (βNTI) was estimated, which represents the difference between the observed βMNTD and mean null expectation in units of standard deviation (999 randomizations). βNTI > +2 indicated heterogeneous selection, and βNTI < −2 indicated homogeneous selection. Finally, the Bray-Curtis-based Raup-Crick metric (RCbray) was calculated, which provided insight into the contribution of stochastic processes to community assembly, when −2 < βNTI < +2. RCbray < −0.95 indicated homogenizing dispersal, −0.95 < RCbray < +0.95 indicated drift, and RCbray > +0.95 indicated dispersal limitation. Randomization was performed 999 times. The averages (%) were subsequently calculated, as they indicated the relative importance of different ecological processes.
3. Results
3.1. Environmental factors of the water columns
Supplementary Figure S1 showed the highest temperatures were in surface water and decreased gradually to the lowest depth at 3,000 m. Salinity was highest in the DCM and lowest at 5 m. Dissolved oxygen increased from 5 m to 75 m, and then volatility decreased at 500 m and 3,000 m. The concentrations of chlorophyll an increased from 5 m to the DCM, then sharply decreased at 200 m and were reduced to zero in the deeper layers. The concentrations of nutrients (total inorganic nitrogen, dissolved inorganic phosphate and silicate) continuously increased from 5 m to 3,000 m depth. More details were described in Shao et al. (2022).
3.2. Community structure of MFW with depths
High-quality sequence numbers of each microeukaryotic sample ranged from 21,172 (S2-200) to 89,829 (S7-3000), with an average of 60,637. Good’s coverage ranged from 97.17% (S8-500) to 99.77% (S7-3000 and S8-3000), with an average of 98.98%. The 16,007 ASVs from 43 samples were obtained, and the ASV number in each sample ranged from 176 (S4-3000) to 1,396 (S8-500), with an average of 699 (Supplementary Table S1). Regarding the distribution patterns of all ASVs in different layers (Supplementary Figure S2), 63 ASVs (0.39%) occurred in all depths, while the total unique ASVs were 2,552, 2,045, 1,561, 2,569, 3,524, and 1,052 in the 5 m, 75 m, DCM, 200 m, 500 m, and 3,000 m layers, respectively.
The relative abundances of the microeukaryotic communities in different depths exhibited different patterns (Figure 2). Dinophyceae was the most abundant microeukaryotic phylum in the investigation area, accounting for more than half of the proportion of the community in almost all samples, except S7-3000. Dinophyceae mainly consisted of the orders Gymnodiniales and Syndiniales. Gymnodiniales contributed more to relative abundance in the 5–200 m layer, while Syndiniales showed the opposite trend. The other top microeukaryotic phyla included Ascomycota (3.99%), Ciliophora (3.91%), Chlorophyta (3.08%), Acantharea (2.80%), and Ichthyodinium (1.78%), which had relative abundances >1%. The other top microeukaryotic orders included Peridiniales (4.69%), Saccharomycetales (3.87%), and Dinophyceae sp. GD1590bp27 (2.29%), which had relative abundances >2%.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1098264/fmicb-14-1098264-HTML/image_m/fmicb-14-1098264-g002.jpg)
Figure 2. Relative abundance of microeukaryotic taxonomic groups at the phylum (top 10 phyla and others) and the order (relative abundance of > 1%) level across 43 samples.
The α-diversity indices of the microeukaryotes in the photic layers showed little fluctuation, which were higher in the 200-and 500-m layers than in the 300-m layer among the aphotic layers. Based on the Wilcox tests, the α-diversity indices of samples in the DCM significantly differed from that of samples in the other layers, while no significant differences were detected between the other layers (Supplementary Figure S3). As per the β-diversity analysis, the PCoA ordinations based on the ASV distributions across the 43 samples indicated that the microeukaryotic community shifted across the six depths (Supplementary Figure S4). LefSe analysis was used to identify differentially abundant taxa among the six depths categorized by microeukaryote successional dynamics (Supplementary Figure S5). The abundances of order Gonyaulacales, and genus Xiphacantha (including species Xiphacanthaalata), were higher in the 5 m layer. The classes, Polycystinea and Spirotrichea, significantly differed in the 200 m layer, while classes, Acantharea and Conidasida (including order Eucoccidiorida), had high abundances in the 500 m layer. The abundance of Oligohymenophorea in the 3,000 m layer was significantly higher than that in the other layers.
The taxonomic compositions of bacteria and their relative abundances are shown in Supplementary Figure S6. At the order level (relative abundances >2%), bacterial sequences were mainly assigned to Betaproteobacteriales (17.48%), Nitrosopumilales (12.81%), Clostridiales (10.73%), Synechococcales (9.30%), SAR11 (7.34%), Bacteroidales (4.11%), Lactobacillales (3.72%), Rhodospirillales (2.85%), Caulobacterales (2.44%), Microtrichales (2.35%), Flavobacteriales (2.31%), and SAR202 (2.18%). Synechococcales were mainly found in the photic zone. In contrast, Nitrosopumilales were mainly found in the aphotic zone. More details of bacterial community structure and α-diversity with depths were described in Shao et al. (2022), which found the 200 m layer had the highest Pieloou’s evenness, Faith’s phylogenetic diversity, Shannon–Wiener and Simpson indices, followed by the DCM layer and 500 m layer, and lowest in the 5-and 75-m layers.
3.3. Co-occurrence networks of bacteria and microeukaryotes in the MFW
The Spearman correlation analysis revealed strong correlations between 16 bacteria and 11 microeukaryote orders, except Sessilida (Figure 3). All microeukaryotic groups correlated with bacteria, except Vestibuliferida. Some bacteria specifically associated with other bacteria. For example, the bacteria SAR 324, SAR 406, SAR202, and Marine Group II highly positively correlated with each other (R > 0.87), as well as positively correlated with the microeukaryotes Syndiniales (R > 0.81) and Arthracanthida (R > 0.64), but negatively correlated with the bacteria Rhodobacterales (R < −0.75), Flavobacteriales (R < −0.74), and Synechococcales (R < −0.71) and the microeukaryotes Suessiales (R < −0.69), Chloropicales (R < −0.64), Thoracosphaerales (R < −0.63), and Dinophyceae sp. GD1590bp27 (R < −0.60). Moreover, Rhodobacterales, Flavobacteriales, and Synechococcales highly positively correlated with each other (R > 0.81) and positively correlated with the microeukaryotes Suessiales (R > 0.68), Thoracosphaerales (R > 0.67), and Gymnodiniales (R > 0.60).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1098264/fmicb-14-1098264-HTML/image_m/fmicb-14-1098264-g003.jpg)
Figure 3. Correlation relationship with the order (relative abundance of >1%) comprising bacteria and microeukaryotes across 43 samples. The network shows strong (R = 0.4) and significant (p < 0.05) Spearman rank correlation coefficients. Red circles represent microeukaryotic orders, blue circles represent bacterial orders. Numbers in the circles are the relative abundance of bacterial/microeukaryotic orders. Red lines indicate positive correlations, blue lines indicate negative correlations.
Significant correlations were detected between different microeukaryotic groups. For example, Syndiniales highly positively correlated with Arthracanthida (R = 0.81). Both microeukaryotes negatively correlated with Gymnodiniales, Suessiales, and Chloropicales (R < −0.64). Suessiales and Chloropicales positively correlated with Gymnodiniales and Thoracosphaerales (R > 0.64). Vestibuliferida only positively correlated with Saccharomycetales.
3.4. Distance-decay relationships in MFW
Spearman correlations comparing Bray-Cutis community dissimilarities with geographic distances between samples revealed no significant positive correlations for microeukaryotic or bacterial communities in 5–200 m and > 500 m layers (Figure 4A). The PERMANOVA tests showed that all, 5–200 m, and > 500 m bacteria and microeukaryotes significantly differed among the six depths. No significant differences were detected across the eight stations and two sections analyzed herein (Supplementary Table S2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1098264/fmicb-14-1098264-HTML/image_m/fmicb-14-1098264-g004.jpg)
Figure 4. (A) Spearman’s rank correlations between the Bray–Curtis dissimilarity of microeukaryotic (blue dots and line) and bacterial (red dots and line) communities and geographic distance in the 5–200 m and > 500 m layers, the n is the number of comparisons. (B) Relative importance of ecological processes in shaping the bacterial and microeukaryotic communities across six depths.
Null model analyses detected higher levels of stochastic processes (dispersal limitation, homogenizing dispersal, and drift) in the 5–200 m layer than in the >500 m layer (Figure 4B). Drift explained the majority of stochastic processes in bacterial (18.67% on average) and microeukaryotic (31.67% on average) community turnover, while homogenizing dispersal only affected (20.00%) bacterial community turnover in the 5 m layer. Among the deterministic processes, homogeneous selection explained 65.67% of bacterial community turnover on average, and it was the most important ecological process, followed by heterogeneous selection, which explained 12.33% of bacterial community turnover. Compared to bacterial community turnover, the relative importance of heterogeneous selection increased in microeukaryotic community turnover (40.83% on average), while homogeneous selection was lower (27.50% on average).
3.5. Driving mechanisms of MFW
Environmental and spatial factors associated with bacterial community were processed and analyzed following previously described methods (Shao et al., 2022). For the microeukaryotic community, the redundancy analysis showed that temperature, salinity, chlorophyll a concentration, and dissolved oxygen concentration were the controlling factors associated with population variation in the 5–200 m layer, while the strongest determinants of community structure in the >500 m layer were silicate and total inorganic nitrogen (TIN), including nitrite, nitrate, and ammonium (Supplementary Figure S7). The Mantel tests indicated that the bacterial and microeukaryotic community structures significantly correlated with almost all environmental factors, especially temperature, TIN, dissolved inorganic phosphate (DIP), and silicate, which exhibited almost no significant effects with five spatial factors (PCNMs, nos. 1–5), except for the 5–200 m bacterial and microeukaryotic communities (Supplementary Table S3).
The heatmap results showed that the bacteria Nitrosopumilales, Synechococcales, Flavobacteriales, SAR202, SAR324, Marine_Group_II, SAR406, and Rhodobacterales and the microeukaryotes Gymnodiniales, Syndiniales, Thoracosphacrales, Suessiales, Arthracanthida, and Chloropicales significantly correlated with temperature, dissolved oxygen, chlorophyll a concentration, DIP, TIN, and silicate (Figure 5). The bacteria Rhodospirillales and the microeukaryotes Gymnodiniales, Peridiniales, Saccharomycetales, Ichthyodinium, Arthracanthida, and Chloropicales significantly correlated with salinity, while the bacteria Betaproteobacteriales and the microeukaryote Saccharomycetales significantly correlated with longitude and latitude. Some microorganisms, like the bacteria Clostridiales, SAR11, Bacteroidales, Lactobacillales, Caulobacterales, and Microtrichales and the microeukaryotes Sessilida and Vestibuliferida, did not significantly correlate with the environmental factors examined in this study.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1098264/fmicb-14-1098264-HTML/image_m/fmicb-14-1098264-g005.jpg)
Figure 5. Heatmap depicting Spearman’s rank correlations between the relative abundances of main bacterial (yellow box) and microeukaryotic (green box) orders and environmental factors. Significant correlation coefficient values were marked in the figure. DIP, dissolved inorganic phosphate; TIN, total inorganic nitrogen.
The variation partitioning analysis showed that environmental effects had a relatively higher contribution on the bacterial and microeukaryotic communities than spatial factors (Figure 6). Comparatively, bacteria were more affected by environmental factors, while microeukaryotes were more affected by spatial factors. For the 5–200 m layer, spatial and environmental factors collectively explained 74.59% of bacterial community variation and 59.07% of microeukaryotic community variation, and remaining variations were unexplained. Spatial factors explained 13.56% of bacterial community variation and 17.83% of microeukaryotic community variation, while environmental factors explained 59.72% of bacterial community variation and 40.95% of microeukaryotic community variation. For the >500 m layer, environmental and spatial explanatory factors for the bacterial and microeukaryotic community variations both decreased; more variations remain unexplained.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1098264/fmicb-14-1098264-HTML/image_m/fmicb-14-1098264-g006.jpg)
Figure 6. Variation partitioning analysis of all, 5–200 m, and > 500 m bacterial and microeukaryotic community composition between environmental and spatial factors.
4. Discussion
4.1. Community diversity and distribution of MFW
Community diversity and distribution of bacteria had been discussed in Shao et al. (2022). In microeukaryotes, Dinophyceae were the most abundant group detected in this study (Figure 2), which was consistent with a previous study on the global ocean scale (de Vargas et al., 2015). Zhao et al. (2020) and Wu et al. (2020) also reported that dinoflagellate was the dominant group in the oligotrophic tropical western Pacific Ocean, while its operational taxonomic unit (OTU) proportion was much lower in aphotic communities than in photic communities. Dinoflagellates comprise both photosynthetic and heterotrophic flagellates, of which, 58% are heterotrophic species in the ocean (Gómez, 2012). Their adaptation to a wide range of environments is reflected by tremendous morphological and trophic diversity (Tailor, 1987). Of these, Gymnodiniales had a higher relative abundance in the 5–200 m layer, while a higher proportion of parasitic taxa, such as Syndiniales, was mainly detected in the aphotic layers (Figure 2). Gymnodiniales are mixotrophic non-thecate dinoflagellates that are generally surface-dwelling. Their appearance in deep waters is likely a result of the amplification of resting cysts, which is as an adaptive strategy for dinoflagellates to survive adverse conditions (Anderson, 1984). Syndiniales are ubiquitous marine parasite taxa that are widely distributed throughout the ocean (Anderson and Harvey, 2020). Our finding supports the hypothesis that microeukaryotic parasites may play key roles in the deep-sea food web (Zoccarato et al., 2016).
In this study, distinct vertical patterns of the microeukaryotic communities along the water column were observed in the oligotrophic tropical western Pacific Ocean (Supplementary Figure S4). Similar patterns were found for the bacterial communities (Shao et al., 2022). Both bacteria and microeukaryotes had their highest α-diversity indices detected in the 200-and 500-m layers (Supplementary Figure S3), wherein dramatic transitions occurred due to temperature, salinity, dissolved oxygen, chlorophyll a concentration, and major nutrients (Supplementary Figure S1). Different taxonomic groups exhibited various vertical patterns, which may be due to their biological characteristics and feeding habits. For instance, the diversity of Radiolarians, including Polycystinea and Acantharea, was higher in the aphotic layers than in the photic layers (Figure 2), and were biomarkers for the 200-and 500-m layers, respectively (Supplementary Figure S5). This may be partially attributed to their intracellular siliceous (Polycystinea) or celestite (Acantharea) skeletons consisting of strontium sulphate (SrSO4), which can thus benefit from transporting and preserving these organisms in deeper layers (Burki et al., 2010). A previous study found that the high diversity of Radiolarians in deep layers may include both local and resident populations, as well as taxa settling from the overlying water, which have been laterally advected from elsewhere (Xu et al., 2017). Ciliophora, mainly including the orders Sessilida and Vestibuliferida, appeared abundantly in the 500-and 3,000-m layers (Figure 2). The class Oligohymenophorea was the biomarker for the 3,000 m layer (Supplementary Figure S5). Sessilida are peritrichs of the class Oligohymenophorea, which are among the most speciose and commonly observed ciliates with >50 genera representing >1,000 species (Li et al., 2008). Vestibuliferida of the class Litostomatea have been reported to be found in the guts of several marine surgeonfish in the oligotrophic tropical western Pacific Ocean, as observed via microscopy (Grim et al., 2002). Some members of Oligohymenophorea and Litostomatea have been suggested to possess hydrogrnosome, belonging to obligatory anaerobes, and other members are facultative anaerobes that have been reportedly found in different suboxic and anoxic conditions (Lynn, 2008), which may explain their distributions in this study.
4.2. Interspecies interactions in the MFW with depths
Recent studies have indicated that interspecies interactions, such as mutualism, cross-feeding, competition, parasitism, predation, and allelopathy, in the MFW are responsible for community structure (Lima-Mendez et al., 2015). In our study, dinoflagellates positively correlated with specific bacterial groups and had the widest selection, where Vestibuliferida (Ciliophora) had almost no correlation with other bacteria or microeukaryotes (Figure 3). Most ciliophora are free-living, some are epi-and endo-commensals, parasites of other unicellular organisms and metazoans, and serve as hosts for bacterial epi-and endo-symbionts, fungi, algae, and other protists (Dziallas et al., 2012), like the ciliophoran Vestibuliferida and the fungus Saccharomycetales.
A previous study demonstrated that phytoplankton had a stronger relationship with bacteria than environmental factors (Pearman et al., 2016). Specific interactions between phytoplankton and bacteria have been documented, many of which are based on the exchange of energy sources and metabolites, including various types of chemical signaling (Ramanan et al., 2016). In photic layers, the bacteria Synechococcales, Flavobacteriales, and Rhodobacterale were highly positively correlated with the dinoflagellates Suessiales (R > 0.68), Thoracosphaerales (R > 0.67), and Gymnodiniales (R > 0.60; Figure 3). The bacterial communities that assemble around marine microphytoplankton are predictably dominated by Rhodobacterales, Flavobacteriales, and families within Gammaproteobacteria (Ferrer-González et al., 2021). A genome analysis of Flavobacteriales isolates identified several genes that encode CAZyme glycoside hydrolases (Fernández-Gomez et al., 2013). Taylor et al. (2014) proposed that Flavobacteriales were degrading and utilizing phytoplankton-derived transparent exopolymer particles. Rhodobacterales are primarily members of the Marine Roseobacter Clade (MRC), which is a metabolically diverse and opportunistic bacterioplankton group that is known to form close associations with phytoplankton blooms (Newton et al., 2010). Synechococcales (phylum Cyanobacteria) are ubiquitously distributed and are the most abundant photosynthetic organisms on Earth (Flombaum et al., 2013). Cyanobacterial symbionts (cyanobionts) have been found in numerous protist groups, including dinoflagellates, tintinnids, radiolarians, amoebae, diatoms, and haptophytes (Foster and Zehr, 2019). The relationships between cyanobionts and protistan hosts are particularly notable, as some nitrogen-fixing cyanobacteria (diazotrophs) play important roles in primary production, especially in nitrogen-limited oligotrophic oceans (Falkowski et al., 1998). However, the nature of cyanobiont-dinoflagellate host consortia remains poorly understood, our findings offer insights into these vacancies in the oligotrophic tropical western Pacific Ocean.
In the aphotic layers, the parasitic dinoflagellate Syndiniales highly positive correlated with the bacteria SAR324, SAR406, and SAR202, and Marine_Group_II (Figure 3; R > 0.81). Positive associations between parasites and heterotrophic bacteria may indicate that cell lysis from eukaryotic parasitism or predation leads to the release of dissolved organic matter, which in turn drives bacterial activity, similar to viral shunt (Jephcott et al., 2016). Members of SAR 324, SAR406, and SAR202 are typical free-living bacteria that serve extremely important roles in the mesopelagic and bathypelagic zones of all oceans (Salazar et al., 2016). Marine_Group_II Euryarchaeota are commonly connected with Syndiniales in oxygen-depleted water columns (Suter et al., 2021). However, each of these bacteria are described only through their ecological distributions or metagenome-assembled genomes and lack cultured representatives, making their functional characterization a challenge. Moreover, Arthracanthida (an order of Radiolarians) also highly positively correlated with the dinoflagellate Syndiniales (R = 0.81; Figure 3), which may result in the highly positive correlations with the bacteria SAR324, SAR406, SAR202, and Marine_Group_II (Figure 3; R > 0.64). Spero and Angel (1991) reported the mutualistic symbionts associated with the plankton rhizarian (Acantharia, Foraminifera, and Radiolaria) and the dinoflagellate genera Amphidinium, Aureodinium, Gloeodinium, Gymnodinium, Gyrodinium, Prorocentrum, Pyrocystis, and Scrippsiella. Our results expand this range of symbionts.
4.3. Factors that regulate the MFW
Generally, dispersal limitation and environmental heterogeneity are the main factors that shape the distributions of microorganisms, which generate a negative correlation between community similarity and geographic distance (Hanson et al., 2012). However, our study showed that geographic distance had little effect on the dispersal of bacterial and microeukaryotic communities over a spatial scale of 1,100 km in the oligotrophic tropical western Pacific Ocean (Figure 4A). The Tara Oceans dataset reported that community dissimilarity significantly positively correlated with geographic distance for all organismal size fractions over a transect of ~6,000 km (De Vargas et al., 2015). Several conceptual frameworks and modeling studies have demonstrated that the dispersal limitation of microbes may lead to a distance-decay relationship, while high dispersal ability weakens the relationship (Hanson et al., 2012). In our study region, Hsu et al. (2017) reported that the core of the NEC (defined as where the westward velocity is >15 cm s−1) occupied a region from 8 to 12.5°N and from 0 to 200 m. To quantitatively analyze the effects of the NEC on the dispersal activity of microorganisms in the six depths, we calculated the contribution of ecological processes to the assembly of MFW in the 5–200 and > 500 m layers. Higher levels of stochastic processes (dispersal limitation, homogenizing dispersal, and drift) were detected in the 5–200 m layer than in the >500 m layer (Figure 4B), which may weakly support our hypothesis that the NEC enhanced the dispersal ability and east–west connectivity of the MFW in the photic layers. Dispersal limitation had a weak effect on the MFW (Figure 4B), which was supported by the variation partitioning analysis, that is, spatial processes explained ~10% of the variance in all, 5–200 m, and > 500 m bacterial and microeukaryotic community composition (Figure 6). Drift explained the majority of stochastic processes (Figure 4B), which was the prevalent community-structuring mechanism in unicellular eukaryotes, and was lower important in structuring prokaryotic organisms (Logares et al., 2018). Some reports found that drift played a moderate role in the assembly of picoeukaryotes and prokaryotes over large spatial scales, such as those present in tropical and subtropical surface oceans (Logares et al., 2020). For deterministic processes, bacterial communities were mainly explained by homogeneous selection, while microeukaryotic communities were mainly explained by heterogeneous and homogeneous selection (Figure 4B). The shallow distance-decay relationships supported the dominance of homogeneous selection (Figure 4A), considering that homogeneous selection contributed to high community similarities (Dini-Andreote et al., 2015). Our results were in accordance with the findings of Logares et al. (2020) and Kong et al. (2022), who found that heterogeneous selection was more important in shaping microeukaryotic communities than bacterial communities, especially in the 3,000 m layer (Figure 4B). They speculated that the differences between the effects of selection may be explained by different modes of adaptation of bacteria and microeukaryotes. Moreover, the dispersal abilities of aquatic organisms and the distances they travel are poorly understood for entire biotas (Bohonak and Jenkins, 2003). In our study, only the bacteria Betaproteobacteriales and the microeukaryote Saccharomycetales significantly correlated with longitude and latitude (Figure 5), which indicated that these microorganisms are weak dispersing organisms (Heino et al., 2015). The fungi Saccharomycetales mainly appeared around stations S7 and S8 (Figure 2), which indicated the presence of organic matter (Korajkic et al., 2015).
In our study, environmental factors contributed more than spatial factors (Figure 6). Almost all environmental factors significantly correlated with the bacterial and microeukaryotic communities (Supplementary Table S3), which indicated that environmental selection played a more important role in shaping the MFW than dispersal limitation in the oligotrophic tropical western Pacific Ocean. Among these factors, temperature had the greatest effect on community, followed by inorganic nutrients (Supplementary Figure S7; Supplementary Table S3). Temperature is an important factor that alters community composition and diversity and is a stronger driver than other environmental factors in shaping microbial community composition, as per a global investigation (Sunagawa et al., 2015). Inorganic nutrients are essential for the growth and development of microorganisms, and different microorganisms adapt to their optimal growth concentrations (Follows and Dutkiewicz, 2011). Due to ocean warming, it is likely that the oligotrophic oceans will be more stratified (Sarmento et al., 2010). Vertical exchange of nutrient between nutrient-deficient upper waters and nutrient-rich deep waters will be reduced as a result of this enhanced stratification in the oligotrophic tropical western Pacific Ocean (Abirami et al., 2021). For different taxonomic groups (Figure 5), the bacteria Nitrosopumilales, SAR202, SAR324, Marine_Group_II, and SAR406 and the microeukaryotes Syndiniales and Arthracanthida negatively correlated with temperature, dissolved oxygen, and the chlorophyll a concentration and positively correlated with inorganic nutrient concentrations, while the bacteria Synechococcales, Flavobacteriales, and Rhodobacterales and the microeukaryotes Gymnodiniales, Thoracosphacrales, Suessiales, and Chloropicales exhibited opposing trends for environmental factors. These results were consistent with interspecies interactions in the MFW (Figure 3), which indicated that the adaptive capability of different taxonomic groups to environmental factors is an important factor that determines the community distribution.
Major bacterial and microeukaryotic groups were more or less regulated by environmental factors in our study (Figure 5), except for the bacteria Clostridiales, SAR11, Bacteroidales, Lactobacillales, Caulobacterales, and Microtrichales and the ciliophoran Sessilida and Vestibuliferida. Members of the SAR11 clade are the most abundant and ubiquitous bacterioplankton found in the oceans worldwide (Yooseph et al., 2010). The relative abundances of Bacteroidales, Clostridiales, and Lactobacillales were determined by the total ammonia and free ammonia concentrations (De Vrieze et al., 2015), which concentrations were not measured in our study. Ciliophora are a free-living group with cilium that are distributed in various habitats and are capable of adapting to different environments (Foissner et al., 2007). Broad niche widths and high ecophysiological tolerances of ciliates are reflected by a low limitation of dispersal (Grossmann et al., 2016). Overall, the responses of each group to environmental factors differed due to their specific adaptation strategies.
5. Conclusion
The interspecies interactions and biogeography of MFW have been focused on a long-standing debate. This research tested three hypotheses of MFW in the oligotrophic tropical western Pacific Ocean. First, we found predominant linkages of MFW were different with depths. Second, we revealed biogeographic patterns of microeukaryotic communities largely matched with bacterial communities. Third, we suggested local environmental factors were more important than spatial factors in shaping the biodiversity and biogeographic of MFW, while NEC had a weak effect. Such knowledge is critical fundamental data for forecasting MFW changes in the face of future global change and future planning of marine protected areas.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/, PRJNA791001 https://www.ncbi.nlm.nih.gov/, PRJNA879911.
Author contributions
QS and DS conceptualized the study and designed the experiment. QS performed the statistical analysis and wrote the manuscript with contributions from DS and CW. DS, CF, and YF participated in the research cruise and collected data. DS and CW approved the final version of the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China (42076122) and the China Ocean Mineral Resources Research and Development Association Program (DY-XZ-02; DY135-E2-3-04).
Acknowledgments
We thank Lei Wang (TIO, MNR) for sharing the environmental data with us and the members of research vessel Dayanghao for their great help during the investigation. We also thank the editor and two reviewers for their valuable comments on earlier versions. We thank LetPub (www.letpub.com) for its linguistic assistance during the preparation of this manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1098264/full#supplementary-material
Footnotes
References
Abirami, B., Radhakrishnan, M., Kumaran, S., and Wilson, A. (2021). Impacts of global warming on marine microbial communities. Sci. Total Environ. 791:147905. doi: 10.1016/j.scitotenv.2021.147905
Anderson, D. M. (1984). “The roles of dormant cysts in toxic dinoflagellate blooms and shellfish toxicity” in Seafood Toxins. American Chemical Society Symposium Series (Washington, United States: American Chemical Society)
Anderson, S. R., and Harvey, E. L. (2020). Temporal variability and ecological interactions of parasitic marine Syndiniales in coastal protist communities. Msphere 5, e00209–e00220. doi: 10.1128/mSphere.00209-20
Azam, F., Fenchel, T., Field, J. G., Gray, J. S., Meyer-Reil, L., and Thingstad, F. (1983). The ecological role of water-column microbes in the sea. Mar. Ecol. Prog. Ser. 10, 257–263. doi: 10.3354/meps010257
Bohonak, A. J., and Jenkins, D. G. (2003). Ecological and evolutionary significance of dispersal by freshwater invertebrates. Ecol. Lett. 6, 783–796. doi: 10.1046/j.1461-0248.2003.00486.x
Borcard, D., and Legendre, P. (2002). All-scale spatial analysis of ecological data by means of principal coordinates of neighbour matrices. Ecol. Model. 153, 51–68. doi: 10.1016/S0304-3800(01)00501-4
Burki, F., Kudryavtsev, A., Matz, M. V., Aglyamova, G. V., Bulman, S., Fiers, M., et al. (2010). Evolution of Rhizaria: new insights from phylogenomic analysis of uncultivated protists. BMC Evol. Biol. 10, 377–318. doi: 10.1186/1471-2148-10-377
Calbet, A. (2008). The trophic roles of microzooplankton in marine systems. ICES J. Mar. Sci. 65, 325–331. doi: 10.1093/icesjms/fsn013
Calbet, A., and Landry, M. R. (1999). Mesozooplankton influences on the microbial food web: direct and indirect trophic interactions in the oligotrophic open ocean. Limnol. Oceanogr. 44, 1370–1380. doi: 10.4319/lo.1999.44.6.1370
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Chen, H., and Boutros, P. C. (2011). VennDiagram: a package for the generation of highly-customizable Venn and Euler diagrams in R. BMC Bioinformatics 12, 1–7. doi: 10.1186/1471-2105-12-35
Dai, S., Zhao, Y. F., Li, X. G., Wang, Z. Y., Zhu, M. L., Liang, J. H., et al. (2020). The seamount effect on phytoplankton in the tropical western Pacific. Mar. Environ. Res. 162:105094. doi: 10.1016/j.marenvres.2020.105094
De Vargas, C., Audic, S., Henry, N., Decelle, J., Mahé, F., Logares, R., et al. (2015). Eukaryotic plankton diversity in the sunlit ocean. Science 348:1261605. doi: 10.1126/science.1261605
De Vrieze, J., Saunders, A. M., He, Y., Fang, J., Nielsen, P. H., Verstraete, W., et al. (2015). Ammonia and temperature determine potential clustering in the anaerobic digestion microbiome. Water Res. 75, 312–323. doi: 10.1016/j.watres.2015.02.025
Dini-Andreote, F., Stegen, J. C., Van Elsas, J. D., and Salles, J. F. (2015). Disentangling mechanisms that mediate the balance between stochastic and deterministic processes in microbial succession. Proc. Natl. Acad. Sci. 112, E1326–E1332. doi: 10.1073/pnas.1414261112
Diniz-Filho, J. A. F., Soares, T. N., Lima, J. S., Dobrovolski, R., Landeiro, V. L., Telles, M. P. D. C., et al. (2013). Mantel test in population genetics. Genet. Mol. Biol. 36, 475–485. doi: 10.1590/S1415-47572013000400002
Djurhuus, A., Closek, C. J., Kelly, R. P., Pitz, K. J., Michisaki, R. P., Starks, H. A., et al. (2020). Environmental DNA reveals seasonal shifts and potential interactions in a marine community. Nat. Commun. 11, 254–259. doi: 10.1038/s41467-019-14105-1
Dziallas, C., Allgaier, M., Monaghan, M. T., and Grossart, H. P. (2012). Act together—implications of symbioses in aquatic ciliates. Front. Microbiol. 3:288. doi: 10.3389/fmicb.2012.00288
Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461. doi: 10.1093/bioinformatics/btq461
Falkowski, P. G., Barber, R. T., and Smetacek, V. (1998). Biogeochemical controls and feedbacks on ocean primary production. Science 281, 200–206. doi: 10.1126/science.281.5374.200
Faust, K., and Raes, J. (2012). Microbial interactions: from networks to models. Nat. Rev. Microbiol. 10, 538–550. doi: 10.1038/nrmicro2832
Fernández-Gomez, B., Richter, M., Schüler, M., Pinhassi, J., Acinas, S. G., González, J. M., et al. (2013). Ecology of marine Bacteroidetes: a comparative genomics approach. ISME J. 7, 1026–1037. doi: 10.1038/ismej.2012.169
Ferrer-González, F. X., Widner, B., Holderman, N. R., Glushka, J., Edison, A. S., Kujawinski, E. B., et al. (2021). Resource partitioning of phytoplankton metabolites that support bacterial heterotrophy. ISME J. 15, 762–773. doi: 10.1038/s41396-020-00811-y
Flombaum, P., José, L. G., Gordillo, R. A., José, R., and Martiny, A. C. (2013). Present and future global distributions of the marine cyanobacteria Prochlorococcus and Synechococcus. Proc. Natl. Acad. Sci. 110, 9824–9829. doi: 10.1073/pnas.1307701110
Foissner, W., Chao, A., and Katz, L. A. (2007). “Diversity and geographic distribution of ciliates (Protista: Ciliophora)” in Protist Diversity and Geographical Distribution, ed. Foissner, W. (Berlin: Springer), 111–129.
Follows, M. J., and Dutkiewicz, S. (2011). Modeling diverse communities of marine microbes. Annu. Rev. Mar. Sci. 3, 427–451. doi: 10.1146/annurev-marine-120709-142848
Foster, R. A., and Zehr, J. P. (2019). Diversity, genomics, and distribution of phytoplankton-cyanobacterium single-cell symbiotic associations. Annu. Rev. Microbiol. 73, 435–456. doi: 10.1146/annurev-micro-090817-062650
Gómez, F. (2012). A quantitative review of the lifestyle, habitat and trophic diversity of dinoflagellates (Dinoflagellata, Alveolata). Syst. Biodivers. 10, 267–275. doi: 10.1080/14772000.2012.721021
Gómez-Rubio, V. (2017). Ggplot2-elegant graphics for data analysis. J. Stat. Softw. 77, 1–3. doi: 10.18637/jss.v077.b02
Grim, J. N., Clements, K. D., and Byfield, T. (2002). New species of Balantidium and Pamcichttdotherus (Ciliophora) inhabiting the intestines of four surgeonfish species from the Tuvalu Islands, Pacific Ocean. J. Eukaryot. Microbiol. 49, 146–153. doi: 10.1111/j.1550-7408.2002.tb00359.x
Grossmann, L., Jensen, M., Heider, D., Jost, S., Glücksman, E., Hartikainen, H., et al. (2016). Protistan community analysis: key findings of a large-scale molecular sampling. ISME J. 10, 2269–2279. doi: 10.1038/ismej.2016.10
Hanson, C. A., Fuhrman, J. A., Horner-Devine, M. C., and Martiny, J. B. H. (2012). Beyond biogeographic patterns: processes shaping the microbial landscape. Nat. Rev. Microbiol. 10, 497–506. doi: 10.1038/nrmicro2795
Heino, J., Melo, A. S., Siqueira, T., Soininen, J., Valanko, S., and Bini, L. M. (2015). Metacommunity organisation, spatial extent and dispersal in aquatic systems: patterns, processes and prospects. Freshw. Biol. 60, 845–869. doi: 10.1111/fwb.12533
Hsu, A. C., Xue, H. J., Chai, F., Xiu, P., and Han, Y. S. (2017). Variability of the Pacific north equatorial current and its implications on Japanese eel (Anguilla japonica) larval migration. Fish. Oceanogr. 26, 251–267. doi: 10.1111/fog.12189
Jephcott, T. G., Alves-De-Souza, C., Gleason, F. H., Van Ogtrop, F. F., Sime-Ngando, T., Karpov, S. A., et al. (2016). Ecological impacts of parasitic chytrids, syndiniales and perkinsids on populations of marine photosynthetic dinoflagellates. Fungal Ecol. 19, 47–58. doi: 10.1016/j.funeco.2015.03.007
Kassambara, A. (2020). Ggpubr: “Ggplot2” Based Publication Ready Plots. http://cran.rapporter.net/web/
Kong, J., Wang, L., Lin, C., Kuang, F. F., Zhou, X. W., Laws, E. A., et al. (2022). Contrasting community assembly mechanisms underlie similar biogeographic patterns of surface microbiota in the tropical North Pacific Ocean. Microbiol. Spectr. 10, e0079821–e0000721. doi: 10.1128/spectrum.00798-21
Korajkic, A., Parfrey, L. W., Mcminn, B. R., Baeza, Y. V., Vanteuren, W., Knight, R., et al. (2015). Changes in bacterial and eukaryotic communities during sewage decomposition in Mississippi river water. Water Res. 69, 30–39. doi: 10.1016/j.watres.2014.11.003
Li, L. F., Song, W. B., Warren, A., Shin, M. K., Chen, Z. G., Ji, D. D., et al. (2008). Reconsideration of the phylogenetic positions of five peritrich genera, vorticella, Pseudovorticella, Zoothamnopsis, Zoothamnium, and Epicarchesium (Ciliophora, Peritrichia, Sessilida), based on small subunit rRNA gene sequences. J. Eukaryot. Microbiol. 55, 448–456. doi: 10.1111/j.1550-7408.2008.00351.x
Liang, Y. T., Zhang, Y. Y., Zhang, Y., Luo, T. W., Rivkin, R. B., and Jiao, N. Z. (2017). Distributions and relationships of virio-and picoplankton in the epi-, meso-and bathypelagic zones of the Western Pacific Ocean. FEMS Microbiol. Ecol. 93:941323. doi: 10.1093/femsec/fiw238
Lima-Mendez, G., Faust, K., Henry, N., Decelle, J., Colin, S., Carcillo, F., et al. (2015). Determinants of community structure in the global plankton interactome. Science 348:1262073. doi: 10.1126/science.1262073
Logares, R., Deutschmann, I. M., Junger, P. C., Giner, C. R., Krabberød, A. K., Schmidt, T. S., et al. (2020). Disentangling the mechanisms shaping the surface ocean microbiota. Microbiome 8, 55–17. doi: 10.1186/s40168-020-00827-8
Logares, R., Tesson, S. V. M., Canbäck, B., Pontarp, M., Hedlund, K., and Rengefors, K. (2018). Contrasting prevalence of selection and drift in the community structuring of bacteria and microbial eukaryotes. Environ. Microbiol. 20, 2231–2240. doi: 10.1111/1462-2920.14265
Lynn, D. (2008). The Ciliated Protozoa: Characterization, Classification, and Guide to the Literature. Berlin: Springer Science & Business Media.
Newton, R. J., Griffin, L. E., Bowles, K. M., Meile, C., Gifford, S., Givens, C. E., et al. (2010). Genome characteristics of a generalist marine bacterial lineage. ISME J. 4, 784–798. doi: 10.1038/ismej.2009.150
Oksanen, J., Kindt, R., Legendre, P., O’hara, B., Stevens, M. H. H., Oksanen, M. J., et al. (2007). The vegan package. Community Ecol. Packag. 10:719. https://www.researchgate.net/
Pearman, J. K., Casas, L., Merle, T., Michell, C., and Irigoien, X. (2016). Bacterial and protist community changes during a phytoplankton bloom. Limnol. Oceanogr. 61, 198–213. doi: 10.1002/lno.10212
Pesant, S., Not, F., Picheral, M., Kandels-Lewis, S., Le Bescot, N., Gorsky, G., et al. (2015). Open science resources for the discovery and analysis of Tara oceans data. Sci. Data 2, 150023–150016. doi: 10.1038/sdata.2015.23
Pinsky, M. L., Selden, R. L., and Kitchel, Z. J. (2020). Climate-driven shifts in marine species ranges: scaling from organisms to communities. Annu. Rev. Mar. Sci. 12, 153–179. doi: 10.1146/annurev-marine-010419-010916
Pomeroy, L. R. (1974). The ocean's food web, a changing paradigm. Bioscience 24, 499–504. doi: 10.2307/1296885
Qiu, B., and Chen, S. M. (2012). Multidecadal Sea level and gyre circulation variability in the northwestern tropical Pacific Ocean. J. Phys. Oceanogr. 42, 193–206. doi: 10.1175/JPO-D-11-061.1
Ramanan, R., Kim, B. H., Cho, D. H., Oh, H. M., and Kim, H. S. (2016). Algae–bacteria interactions: evolution, ecology and emerging applications. Biotechnol. Adv. 34, 14–29. doi: 10.1016/j.biotechadv.2015.12.003
Richter, D. J., Watteaux, R., Vannier, T., Leconte, J., Frémont, P., Reygondeau, G., et al. (2020). Genomic evidence for global ocean plankton biogeography shaped by large-scale current systems. BioRxiv 11:e78129. doi: 10.7554/eLife.78129
Salazar, G., Cornejo-Castillo, F. M., Benítez-Barrios, V., Fraile-Nuez, E., Álvarez-Salgado, X. A., Duarte, C. M., et al. (2016). Global diversity and biogeography of deep-sea pelagic prokaryotes. ISME J. 10, 596–608. doi: 10.1038/ismej.2015.137
Sarmento, H., Montoya, J. M., Vázquez-Domínguez, E., Vaqué, D., and Gasol, J. M. (2010). Warming effects on marine microbial food web processes: how far can we go when it comes to predictions? Philos. Trans. R. Soc. B Biol. Sci. 365, 2137–2149. doi: 10.1098/rstb.2010.0045
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09
Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12, R60–R18. doi: 10.1186/gb-2011-12-6-r60
Seymour, J. R., Amin, S. A., Raina, J. B., and Stocker, R. (2017). Zooming in on the phycosphere: the ecological interface for phytoplankton–bacteria relationships. Nat. Microbiol. 2, 1–12. doi: 10.1038/nmicrobiol.2017.65
Shao, Q. W., Sun, D., Fang, C., Feng, Y. Z., and Wang, C. S. (2022). Biodiversity and biogeography of abundant and rare microbial assemblages in the Western subtropical Pacific Ocean. Front. Microbiol. 13:839562. doi: 10.3389/fmicb.2022.839562
Sherr, E., and Sherr, B. F. (1988). Role of microbes in pelagic food webs: a revised concept. Limnol. Oceanogr. 33, 1225–1227. doi: 10.4319/lo.1988.33.5.1225
Spero, H. J., and Angel, D. L. (1991). Planktonic sarcodines: microhabitat for oceanic dinoflagellates. J. Phycol. 27, 187–195. doi: 10.1111/j.0022-3646.1991.00187.x
Stegen, J. C., Lin, X. J., Fredrickson, J. K., Chen, X. Y., Kennedy, D. W., Murray, C. J., et al. (2013). Quantifying community assembly processes and identifying features that impose them. ISME J. 7, 2069–2079. doi: 10.1038/ismej.2013.93
Stoeck, T., Bass, D., Nebel, M., Christen, R., Jones, M. D. M., Breiner, H. W., et al. (2010). Multiple marker parallel tag environmental DNA sequencing reveals a highly complex eukaryotic community in marine anoxic water. Mol. Ecol. 19, 21–31. doi: 10.1111/j.1365-294X.2009.04480.x
Sunagawa, S., Coelho, L. P., Chaffron, S., Kultima, J. R., Labadie, K., Salazar, G., et al. (2015). Structure and function of the global ocean microbiome. Science 348:1261359. doi: 10.1126/science.1261359
Suter, E. A., Pachiadaki, M., Taylor, G. T., and Edgcomb, V. P. (2021). Eukaryotic parasites are integral to a productive microbial food web in oxygen-depleted waters. Front. Microbiol. 12:764605. doi: 10.3389/fmicb.2021.764605
Tailor, F. J. R. (1987). The Biology of Dinoflagellates, Botanical Monographs. 21, Malden, MA: Blackwell Science.
Taylor, J. D., Cottingham, S. D., Billinge, J., and Cunliffe, M. (2014). Seasonal microbial community dynamics correlate with phytoplankton-derived polysaccharides in surface coastal waters. ISME J. 8, 245–248. doi: 10.1038/ismej.2013.178
Wang, M. M., Ma, Y. Y., Feng, C. H., Cai, L., and Li, W. (2020). Diversity of pelagic and benthic bacterial assemblages in the western Pacific Ocean. Front. Microbiol. 11:1730. doi: 10.3389/fmicb.2020.01730
Wei, Y. Q., Chen, Z. G., Guo, C. C., Zhong, Q., Wu, C., and Sun, J. (2020). Physiological and ecological responses of photosynthetic processes to oceanic properties and phytoplankton communities in the oligotrophic western Pacific Ocean. Front. Microbiol. 11:1774. doi: 10.3389/fmicb.2020.01774
Wei, Y. Q., Sun, J., Chen, Z. G., Zhang, Z. Y., Zhang, G. C., and Liu, X. F. (2021). Significant contribution of picoplankton size fraction to biogenic silica standing stocks in the Western Pacific Ocean. Prog. Oceanogr. 192:102516. doi: 10.1016/j.pocean.2021.102516
Worm, B., Sandow, M., Oschlies, A., Lotze, H. K., and Myers, R. A. (2005). Global patterns of predator diversity in the open oceans. Science 309, 1365–1369. doi: 10.1126/science.1113399
Wu, P. F., Li, D. X., Kong, L. F., Li, Y. Y., Zhang, H., Xie, Z. X., et al. (2020). The diversity and biogeography of microeukaryotes in the euphotic zone of the northwestern Pacific Ocean. Sci. Total Environ. 698:134289. doi: 10.1016/j.scitotenv.2019.134289
Xu, Z. M., Cheung, S. Y., Endo, H., Xia, X. M., Wu, W. X., Chen, B. Z., et al. (2022). Disentangling the ecological processes shaping the latitudinal pattern of phytoplankton communities in the Pacific Ocean. Msystems 7, e0120321–e0101221. doi: 10.1128/msystems.01203-21
Xu, D. P., Li, R., Hu, C., Sun, P., Jiao, N. Z., and Warren, A. (2017). Microbial eukaryote diversity and activity in the water column of the South China Sea based on DNA and RNA high throughput sequencing. Front. Microbiol. 8:1121. doi: 10.3389/fmicb.2017.01121
Yooseph, S., Nealson, K. H., Rusch, D. B., Mccrow, J. P., Dupont, C. L., Kim, M., et al. (2010). Genomic and functional adaptation in surface ocean planktonic prokaryotes. Nature 468, 60–66. doi: 10.1038/nature09530
Zeldis, J. R., and Décima, M. (2020). Mesozooplankton connect the microbial food web to higher trophic levels and vertical export in the New Zealand subtropical convergence zone. Deep-Sea Res. I Oceanogr. Res. Pap. 155:103146. doi: 10.1016/j.dsr.2019.103146
Zhao, F., Filker, S., Xu, K. D., Huang, P. P., and Zheng, S. (2020). Microeukaryote communities exhibit phyla-specific distance-decay patterns and an intimate link between seawater and sediment habitats in the Western Pacific Ocean. Deep-Sea Res. I Oceanogr. Res. Pap. 160:103279. doi: 10.1016/j.dsr.2020.103279
Zhou, J. Z., and Ning, D. L. (2017). Stochastic community assembly: does it matter in microbial ecology? Microbiol. Mol. Biol. Rev. 81, e00002–e00017. doi: 10.1128/MMBR.00002-17
Zhou, J., Richlen, M. L., Sehein, T. R., Kulis, D. M., Anderson, D. M., and Cai, Z. (2018). Microbial community structure and associations during a marine dinoflagellate bloom. Front. Microbiol. 9:1201. doi: 10.3389/fmicb.2018.01201
Zoccarato, L., Pallavicini, A., Cerino, F., Umani, S. F., and Celussi, M. (2016). Water mass dynamics shape Ross Sea protist communities in mesopelagic and bathypelagic layers. Prog. Oceanogr. 149, 16–26. doi: 10.1016/j.pocean.2016.10.003
Keywords: tropical western Pacific Ocean, microbial food web, interspecies interaction, biogeographic pattern, driving mechanism, North Equatorial Current
Citation: Shao Q, Sun D, Fang C, Feng Y and Wang C (2023) Microbial food webs share similar biogeographic patterns and driving mechanisms with depths in oligotrophic tropical western Pacific Ocean. Front. Microbiol. 14:1098264. doi: 10.3389/fmicb.2023.1098264
Edited by:
Michaela M. Salcher, Biology Centre of the Czech Academy of Sciences, CzechiaReviewed by:
Yu Wang, Xiamen University, ChinaKasia Piwosz, National Marine Fisheries Research Institute, Poland
Copyright © 2023 Shao, Sun, Fang, Feng and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dong Sun, ✉ sund@sio.org.cn