Niche adaptation promoted the evolutionary diversification of tiny ocean predators.

Unicellular eukaryotic predators play a crucial role in the functioning of the ocean ecosystem by recycling nutrients and energy that are channeled to upper trophic levels. Traditionally, these evolutionarily diverse organisms have been combined into a single functional group (heterotrophic flagellates), overlooking their organismal differences. Here, we investigated four evolutionarily related species belonging to one cosmopolitan group of uncultured marine picoeukaryotic predators: marine stramenopiles (MAST)-4 (species A, B, C, and E). Co-occurrence and distribution analyses in the global surface ocean indicated contrasting patterns in MAST-4A and C, suggesting adaptation to different temperatures. We then investigated whether these spatial distribution patterns were mirrored by MAST-4 genomic content using single-cell genomics. Analyses of 69 single cells recovered 66 to 83% of the MAST-4A/B/C/E genomes, which displayed substantial interspecies divergence. MAST-4 genomes were similar in terms of broad gene functional categories, but they differed in enzymes of ecological relevance, such as glycoside hydrolases (GHs), which are part of the food degradation machinery in MAST-4. Interestingly, MAST-4 species featuring a similar GH composition (A and C) coexcluded each other in the surface global ocean, while species with a different set of GHs (B and C) appeared to be able to coexist, suggesting further niche diversification associated with prey digestion. We propose that differential niche adaptation to temperature and prey type has promoted adaptive evolutionary diversification in MAST-4. We show that minute ocean predators from the same phylogenetic group may have different biogeography and genomic content, which needs to be accounted for to better comprehend marine food webs.

Functional redundancy of seasonal vitamin B12 biosynthesis pathways in coastal marine microbial communities.

Vitamin B12 (cobalamin) is a major cofactor required by most marine microbes, but only produced by a few prokaryotes in the ocean, which is globally B12 -depleted. Despite the ecological importance of B12 , the seasonality of B12 metabolisms and the organisms involved in its synthesis in the ocean remain poorly known. Here we use metagenomics to assess the monthly dynamics of B12 -related pathways and the functional diversity of associated microbial communities in the coastal NW Mediterranean Sea over 7 years. We show that genes related to potential B12 metabolisms were characterized by an annual succession of different organisms carrying distinct production pathways. During the most productive winter months, archaea (Nitrosopumilus and Nitrosopelagicus) were the main contributors to B12 synthesis potential through the anaerobic pathway (cbi genes). In turn, Alphaproteobacteria (HIMB11, UBA8309, Puniceispirillum) contributed to B12 synthesis potential in spring and summer through the aerobic pathway (cob genes). Cyanobacteria could produce pseudo-cobalamin from spring to autumn. Finally, we show that during years with environmental perturbations, the organisms usually carrying B12 synthesis genes were replaced by others having the same gene, thus maintaining the potential for B12 production. Such ecological insurance could contribute to the long-term functional resilience of marine microbial communities exposed to contrasting inter-annual environmental conditions.

Seasonality of biogeochemically relevant microbial genes in a coastal ocean microbiome.

Microbes drive the biogeochemical cycles of marine ecosystems through their vast metabolic diversity. While we have a fairly good understanding of the spatial distribution of these metabolic processes in various ecosystems, less is known about their seasonal dynamics. We investigated the annual patterns of 21 biogeochemical relevant functions in an oligotrophic coastal ocean site by analysing the presence of key genes, analysing high-rank gene taxonomy and the dynamics of nucleotide variants. Most genes presented seasonality: photoheterotrophic processes were enriched during spring, phosphorous-related genes were dominant during summer, coinciding with potential phosphate limitation, and assimilatory nitrate reductases appeared mostly during summer and autumn, correlating negatively with nitrate availability. Additionally, we identified the main taxa driving each function at each season and described the role of underrecognized taxa such as Litoricolaceae in carbon fixation (rbcL), urea degradation (ureC), and CO oxidation (coxL). Finally, the seasonality of single variants of some families presented a decoupling between the taxonomic abundance patterns and the functional gene patterns, implying functional specialization of the different genera. Our study unveils the seasonality of key biogeochemical functions and the main taxonomic groups that harbour these relevant functions in a coastal ocean ecosystem.

Role of Bacterial Community Composition as a Driver of the Small-Sized Phytoplankton Community Structure in a Productive Coastal System.

We present here the first detailed description of the seasonal patterns in bacterial community composition (BCC) in shelf waters off the Ría de Vigo (Spain), based on monthly samplings during 2 years. Moreover, we studied the relationship between bacterial and small-sized eukaryotic community composition to identify potential biotic interactions among components of these two communities. Bacterial operational taxonomic unit (OTU) richness and diversity systematically peaked in autumn-winter, likely related to low resource availability during this period. BCC showed seasonal and vertical patterns, with Rhodobacteraceae and Flavobacteriaceae families dominating in surface waters, and SAR11 clade dominating at the base of the photic zone (30 m depth). BCC variability was significantly explained by environmental variables (e.g., temperature of water, solar radiation, or dissolved organic matter). Interestingly, a strong and significant correlation was found between BCC and small-sized eukaryotic community composition (ECC), which suggests that biotic interactions may play a major role as structuring factors of the microbial plankton in this productive area. In addition, co-occurrence network analyses revealed strong and significant, mostly positive, associations between bacteria and small-sized phytoplankton. Positive associations likely result from mutualistic relationships (e.g., between Dinophyceae and Rhodobacteraceae), while some negative correlations suggest antagonistic interactions (e.g., between Pseudo-nitzchia sp. and SAR11). These results support the key role of biotic interactions as structuring factors of the small-sized eukaryotic community, mostly driven by positive associations between small-sized phytoplankton and bacteria.

18S rDNA gene metabarcoding of microeukaryotes and epi-endophytes in the holobiome of seven species of large brown algae.

Brown algae (Phaeophyceae) are habitat-forming species in coastal ecosystems and include kelp forests and seaweed beds that support a wide diversity of marine life. Host-associated microbial communities are an integral part of phaeophyte biology, and whereas the bacterial microbial partners have received considerable attention, the microbial eukaryotes associated with brown algae have hardly been studied. Here, we used broadly targeted "pan-eukaryotic" primers (metabarcoding) to investigate brown algal-associated eukaryotes (the eukaryome). Using this approach, we aimed to investigate the eukaryome of seven large brown algae that are important and common species in coastal ecosystems. We also aimed to assess whether these macroalgae harbor novel eukaryotic diversity and to ascribe putative functional roles to the host-associated eukaryome based on taxonomic affiliation and phylogenetic placement. We detected a significant diversity of microeukaryotic and algal lineages associated with the brown algal species investigated. The operational taxonomic units (OTUs) were taxonomically assigned to 10 of the eukaryotic major supergroups, including taxonomic groups known to be associated with seaweeds as epibionts, endobionts, parasites, and commensals. Additionally, we revealed previously unrecorded sequence types, including novel phaeophyte OTUs, particularly in the Fucus spp. samples, that may represent fucoid genomic variants, sequencing artifacts, or undescribed epi-/endophytes. Our results provide baseline data and technical insights that will be useful for more comprehensive seaweed eukaryome studies investigating the evidently lineage-rich and functionally diverse symbionts of brown algae.

Diversity of eukaryotic and prokaryotic microbiota revealed by metabarcoding in Neotropical floodplain lakes.

The diversity of eukaryotic and prokaryotic communities has been assessed by morphological and genetic approaches, which are used to characterize the microbiota in different environments. Here, planktonic prokaryotic and eukaryotic communities of the Araguaia River, located in the Central region of Brazil, were analyzed based on metabarcoding analysis of rRNA genes to evaluate the diversity of these groups in tropical floodplain lakes. Also, we tested their spatial concordance throughout the Araguaia river. Water samples were collected from 8 floodplain lakes in Araguaia River. The 16S and 18S rRNA genes were amplified and sequenced using Illumina MiSeq. For eukaryotes, 34,242 merged reads were obtained and 225 distinct OTUs were delineated, of which 106 OTUs were taxonomically classified. For prokaryotes, 26,426 sequences were obtained and 351 OTUs were detected. Of them, 231 were classified in at least one taxonomic category. The most representative eukaryotes belonged to Ciliophora, Chlorophyta and Charophyta. The prokaryotic phylum with the most OTUs classified were Proteobacteria, Actinobacteria and Bacteroidetes. The lakes did not show spatial concordance when comparing the similarity between their microbiota. The knowledge of freshwater biodiversity using DNA sequencing for important rivers, such as Araguaia River, can improve microbiota inventories of tropical biodiversity hotspots.

Shared and contrasting associations in the dynamic nano- and picoplankton communities of two close but contrasting sites from the Bay of Biscay.

Pico- and nanoplankton are key players in the marine ecosystems due to their implication in the biogeochemical cycles, nutrient recycling and the pelagic food webs. However, the specific dynamics and niches of most bacterial, archaeal and eukaryotic plankton remain unknown, as well as the interactions between them. Better characterization of these is critical for understanding and predicting ecosystem functioning under anthropogenic pressures. We used environmental DNA metabarcoding across a 6-year time series to explore the structure and seasonality of pico- and nanoplankton communities in two sites of the Bay of Biscay, one coastal and one offshore, and construct association networks to reveal potential keystone and connector taxa. Temporal trends in alpha diversity were similar between the two sites, and concurrent communities more similar than within the same site at different times. However, we found differences between the network topologies of the two sites, with both shared and site-specific keystones and connectors. For example, Micromonas, with lower abundance in the offshore site is a keystone here, indicating a stronger effect of associations such as resource competition. This study provides an example of how time series and association network analysis can reveal how similar communities may function differently despite being geographically close.

Microbial population genomes from the Amazon River reveal possible modulation of the organic matter degradation process in tropical freshwaters.

Rivers connect the carbon cycle in land with that in aquatic ecosystems by transporting and transforming terrestrial organic matter (TeOM). The Amazon River receives huge loads of TeOM from the surrounding rainforest, promoting a substantial microbial heterotrophic activity and consequently, CO2 outgassing. In the Amazon River, microbes degrade up to 55% of the lignin present in the TeOM. Yet, the main microbial genomes involved in TeOM degradation were unknown. Here, we characterize 51 population genomes (PGs) representing some of the most abundant microbes in the Amazon River deriving from 106 metagenomes. The 51 reconstructed PGs are among the most abundant microbes in the Amazon River, and 53% of them are not able to degrade TeOM. Among the PGs capable of degrading TeOM, 20% were exclusively cellulolytic, while the others could also oxidize lignin. The transport and consumption of lignin oxidation byproducts seemed to be decoupled from the oxidation process, being apparently performed by different groups of microorganisms. By connecting the genomic features of abundant microbes in the Amazon River with the degradation machinery of TeOM, we suggest that a complex microbial consortium could explain the quick turnover of TeOM previously observed in this ecosystem.

Sinking particles promote vertical connectivity in the ocean microbiome.

The sinking of organic particles formed in the photic layer is a main vector of carbon export into the deep ocean. Although sinking particles are heavily colonized by microbes, so far it has not been explored whether this process plays a role in transferring prokaryotic diversity from surface to deep oceanic layers. Using Illumina sequencing of the 16S rRNA gene, we explore here the vertical connectivity of the ocean microbiome by characterizing marine prokaryotic communities associated with five different size fractions and examining their compositional variability from surface down to 4,000 m across eight stations sampled in the Atlantic, Pacific, and Indian Oceans during the Malaspina 2010 Expedition. Our results show that the most abundant prokaryotes in the deep ocean are also present in surface waters. This vertical community connectivity seems to occur predominantly through the largest particles because communities in the largest size fractions showed the highest taxonomic similarity throughout the water column, whereas free-living communities were more isolated vertically. Our results further suggest that particle colonization processes occurring in surface waters determine to some extent the composition and biogeography of bathypelagic communities. Overall, we postulate that sinking particles function as vectors that inoculate viable particle-attached surface microbes into the deep-sea realm, determining to a considerable extent the structure, functioning, and biogeography of deep ocean communities.

Perspectives from Ten Years of Protist Studies by High-Throughput Metabarcoding.

During the last decade, high-throughput metabarcoding became routine for analyzing protistan diversity and distributions in nature. Amid a multitude of exciting findings, scientists have also identified and addressed technical and biological limitations, although problems still exist for inference of meaningful taxonomic and ecological knowledge based on short DNA sequences. Given the extensive use of this approach, it is critical to settle our understanding on its strengths and weaknesses and to synthesize up-to-date methodological and conceptual trends. This article summarizes key scientific and technical findings, and identifies current and future directions in protist research that uses metabarcoding.

Impact of grazing, resource availability and light on prokaryotic growth and diversity in the oligotrophic surface global ocean.

The impact of grazing, resource competition and light on prokaryotic growth and taxonomic composition in subtropical and tropical surface waters were studied through 10 microcosm experiments conducted between 30°N and 30°S in the Atlantic, Pacific and Indian oceans. Under natural sunlight conditions, significant changes in taxonomic composition were only observed after the reduction of grazing by sample filtration in combination with a decrease in resource competition by sample dilution. Sunlight exposure significantly reduced prokaryote growth (11 ± 6%) and community richness (14 ± 4%) compared to continuous darkness but did not significantly change community composition. The largest growth inhibition after sunlight exposure occurred at locations showing deep mixed layers. The reduction of grazing had an expected and significant positive effect on growth, but caused a significant decrease in community richness (16 ± 6%), suggesting that the coexistence of many different OTUs is partly promoted by the presence of predators. Dilution of the grazer-free prokaryotic community significantly enhanced growth at the level of community, but consistently and sharply reduced the abundance of Prochlorococcus and SAR11 populations. The decline of these oligotrophic bacterial taxa following an increase in resource availability is consistent with their high specialization for exploiting the limited resources available in the oligotrophic warm ocean.

Primer Design for an Accurate View of Picocyanobacterial Community Structure by Using High-Throughput Sequencing.

High-throughput sequencing (HTS) of the 16S rRNA gene has been used successfully to describe the structure and dynamics of microbial communities. Picocyanobacteria are important members of bacterioplankton communities, and, so far, they have predominantly been targeted using universal bacterial primers, providing a limited resolution of the picocyanobacterial community structure and dynamics. To increase such resolution, the study of a particular target group is best approached with the use of specific primers. Here, we aimed to design and evaluate specific primers for aquatic picocyanobacterial genera to be used with high-throughput sequencing. Since the various regions of the 16S rRNA gene have different degrees of conservation in different bacterial groups, we therefore first determined which hypervariable region of the 16S rRNA gene provides the highest taxonomic and phylogenetic resolution for the genera Synechococcus, Prochlorococcus, and Cyanobium An in silico analysis showed that the V5, V6, and V7 hypervariable regions appear to be the most informative for this group. We then designed primers flanking these hypervariable regions and tested them in natural marine and freshwater communities. We successfully detected that most (97%) of the obtained reads could be assigned to picocyanobacterial genera. We defined operational taxonomic units as exact sequence variants (zero-radius operational taxonomic units [zOTUs]), which allowed us to detect higher genetic diversity and infer ecologically relevant information about picocyanobacterial community composition and dynamics in different aquatic systems. Our results open the door to future studies investigating picocyanobacterial diversity in aquatic systems.IMPORTANCE The molecular diversity of the aquatic picocyanobacterial community cannot be accurately described using only the available universal 16S rRNA gene primers that target the whole bacterial and archaeal community. We show that the hypervariable regions V5, V6, and V7 of the 16S rRNA gene are better suited to study the diversity, community structure, and dynamics of picocyanobacterial communities at a fine scale using Illumina MiSeq sequencing. Due to its variability, it allows reconstructing phylogenies featuring topologies comparable to those generated when using the complete 16S rRNA gene sequence. Further, we successfully designed a new set of primers flanking the V5 to V7 region whose specificity for picocyanobacterial genera was tested in silico and validated in several freshwater and marine aquatic communities. This work represents a step forward for understanding the diversity and ecology of aquatic picocyanobacteria and sets the path for future studies on picocyanobacterial diversity.

Quantifying long-term recurrence in planktonic microbial eukaryotes.

How much temporal recurrence is present in microbial assemblages is still an unanswered ecological question. Even though marked seasonal changes have been reported for whole microbial communities, less is known on the dynamics and seasonality of individual taxa. Here, we aim at understanding microbial recurrence at three different levels: community, taxonomic group and operational taxonomic units (OTUs). For that, we focused on a model microbial eukaryotic community populating a long-term marine microbial observatory using 18S rRNA gene data from two organismal size fractions: the picoplankton (0.2-3 µm) and the nanoplankton (3-20 µm). We have developed an index to quantify recurrence in particular taxa. We found that community structure oscillated systematically between two main configurations corresponding to winter and summer over the 10 years studied. A few taxonomic groups such as Mamiellophyceae or MALV-III presented clear recurrence (i.e., seasonality), whereas 13%-19% of the OTUs in both size fractions, accounting for ~40% of the relative abundance, featured recurrent dynamics. Altogether, our work links long-term whole community dynamics with that of individual OTUs and taxonomic groups, indicating that recurrent and non-recurrent changes characterize the dynamics of microbial assemblages.

Assessing the viral content of uncultured picoeukaryotes in the global-ocean by single cell genomics.

Viruses are the most abundant biological entities on Earth and have fundamental ecological roles in controlling microbial communities. Yet, although their diversity is being increasingly explored, little is known about the extent of viral interactions with their protist hosts as most studies are limited to a few cultivated species. Here, we exploit the potential of single-cell genomics to unveil viral associations in 65 individual cells of 11 essentially uncultured stramenopiles lineages sampled during the Tara Oceans expedition. We identified viral signals in 57% of the cells, covering nearly every lineage and with narrow host specificity signal. Only seven out of the 64 detected viruses displayed homologies to known viral sequences. A search for our viral sequences in global ocean metagenomes showed that they were preferentially found at the DCM and within the 0.2-3 µm size fraction. Some of the viral signals were widely distributed, while others geographically constrained. Among the viral signals we detected an endogenous mavirus virophage potentially integrated within the nuclear genome of two distant uncultured stramenopiles. Virophages have been previously reported as a cell's defence mechanism against other viruses, and may therefore play an important ecological role in regulating protist populations. Our results point to single-cell genomics as a powerful tool to investigate viral associations in uncultured protists, suggesting a wide distribution of these relationships, and providing new insights into the global viral diversity.

Higher contribution of globally rare bacterial taxa reflects environmental transitions across the surface ocean.

Microbial taxa range from being ubiquitous and abundant across space to extremely rare and endemic, depending on their ecophysiology and on different processes acting locally or regionally. However, little is known about how cosmopolitan or rare taxa combine to constitute communities and whether environmental variations promote changes in their relative abundances. Here we identified the Spatial Abundance Distribution (SpAD) of individual prokaryotic taxa (16S rDNA-defined Operational Taxonomic Units, OTUs) across 108 globally-distributed surface ocean stations. We grouped taxa based on their SpAD shape ("normal-like"- abundant and ubiquitous; "logistic"- globally rare, present in few sites; and "bimodal"- abundant only in certain oceanic regions), and investigated how the abundance of these three categories relates to environmental gradients. Most surface assemblages were numerically dominated by a few cosmopolitan "normal-like" OTUs, yet there was a gradual shift towards assemblages dominated by "logistic" taxa in specific areas with productivity and temperature differing the most from the average conditions in the sampled stations. When we performed the SpAD categorization including additional habitats (deeper layers and particles of varying sizes), the SpAD of many OTUs changed towards fewer "normal-like" shapes, and OTUs categorized as globally rare in the surface ocean became abundant. This suggests that understanding the mechanisms behind microbial rarity and dominance requires expanding the context of study beyond local communities and single habitats. We show that marine bacterial communities comprise taxa displaying a continuum of SpADs, and that variations in their abundances can be linked to habitat transitions or barriers that delimit the distribution of community members.

Protist Diversity and Seasonal Dynamics in Skagerrak Plankton Communities as Revealed by Metabarcoding and Microscopy.

Protist community composition and seasonal dynamics are of major importance for the production of higher trophic levels, such as zooplankton and fish. Our aim was to reveal how the protist community in the Skagerrak changes through the seasons by combining high-throughput sequencing and microscopy of plankton collected monthly over two years. The V4 region of the 18S rRNA gene was amplified by eukaryote universal primers from the total RNA/cDNA. We found a strong seasonal variation in protist composition and proportional abundances, and a difference between two depths within the euphotic zone. Highest protist richness was found in late summer-early autumn, and lowest in winter. Temperature was the abiotic factor explaining most of the variation in diversity. Dinoflagellates was the most abundant and diverse group followed by ciliates and diatoms. We found about 70 new taxa recorded for the first time in the Skagerrak. The seasonal pattern in relative read abundance of major phytoplankton groups was well in accordance with microscopical biovolumes. This is the first metabarcoding study of the protist plankton community of all taxonomic groups and through seasons in the Skagerrak, which may serve as a baseline for future surveys to reveal effects of climate and environmental changes.

Contrasting prevalence of selection and drift in the community structuring of bacteria and microbial eukaryotes.

Whether or not communities of microbial eukaryotes are structured in the same way as bacteria is a general and poorly explored question in ecology. Here, we investigated this question in a set of planktonic lake microbiotas in Eastern Antarctica that represent a natural community ecology experiment. Most of the analysed lakes emerged from the sea during the last 6000 years, giving rise to waterbodies that originally contained marine microbiotas and that subsequently evolved into habitats ranging from freshwater to hypersaline. We show that habitat diversification has promoted selection driven by the salinity gradient in bacterial communities (explaining ∼ 72% of taxa turnover), while microeukaryotic counterparts were predominantly structured by ecological drift (∼72% of the turnover). Nevertheless, we also detected a number of microeukaryotes with specific responses to salinity, indicating that albeit minor, selection has had a role in the structuring of specific members of their communities. In sum, we conclude that microeukaryotes and bacteria inhabiting the same communities can be structured predominantly by different processes. This should be considered in future studies aiming to understand the mechanisms that shape microbial assemblages.

Seasonal succession of small planktonic eukaryotes inhabiting surface waters of a coastal upwelling system.

Small eukaryotes (0.2-20 µm cell-size) represent a significant fraction of the microbial plankton community in shelf waters of NW-Spain. The community composition of small eukaryotes living at the surface and at the base of the photic zone was analysed by means of 18S rDNA high-throughput sequencing on a circa-monthly basis over a 23 months period. Ostreococcus was the most abundant taxon in surface waters, showing marked peaks in read abundance in spring and late summer, while Syndiniales dominated at the base of the photic zone. A well-defined seasonal pattern of community composition, linked to the succession of the dominant taxa, was found in surface waters. Seasonality was less apparent at the base of the euphotic zone. Temporal changes in abiotic factors significantly correlated with changes in community composition in surface (r = 0.71) and at the base of the photic zone (r = 0.38). Changes in community composition significantly correlated with changes in community function-related variables (including biomass, primary production and respiration) only in surface water (r = 0.36). Co-occurrence network analyses revealed 45 significant interspecies associations among the 50 most abundant taxa with highly connected OTUs belonging to cryptophyceans. The network topology, with small-world characteristics, suggests a stabilizing role of biotic interactions to environmental disturbance.

Delineating closely related dinoflagellate lineages using phylotranscriptomics.

Recently radiated dinoflagellates Apocalathium aciculiferum (collected in Lake Erken, Sweden), Apocalathium malmogiense (Baltic Sea) and Apocalathium aff. malmogiense (Highway Lake, Antarctica) represent a lineage with an unresolved phylogeny. We determined their phylogenetic relationships using phylotranscriptomics based on 792 amino acid sequences. Our results showed that A. aciculiferum diverged from the other two closely related lineages, consistent with their different morphologies in cell size, relative cell length and presence of spines. We hypothesized that A. aff. malmogiense and A. malmogiense, which inhabit different hemispheres, are evolutionarily more closely related because they diverged from a marine common ancestor, adapting to a wide salinity range, while A. aciculiferum colonized a freshwater habitat, by acquiring adaptations to this environment, in particular, salinity intolerance. We show that phylotranscriptomics can resolve the phylogeny of recently diverged protists. This has broad relevance, given that many phytoplankton species are morphologically very similar, and single genes sometimes lack the information to determine species' relationships.

Abundant and rare picoeukaryotic sub-communities present contrasting patterns in the epipelagic waters of marginal seas in the northwestern Pacific Ocean.

In this work, they compared patterns of abundant and rare picoeukaryotic sub-communities in the epipelagic waters (surface and 40-75 m depth subsurface layers) of the East and South China Seas across seasons via 454 pyrosequencing of the V4 region of 18S rDNA. They also examined the relative effects of environmental filtering, dispersal limitations and seasonality on community assembly. Their results indicated that (i) in the surface layer, abundant taxa are primarily influenced by dispersal limitations and rare taxa are primarily influenced by environmental filtering, whereas (ii) in the subsurface layer, both abundant and rare sub-communities are only weakly influenced by environmental filtering but are strongly influenced by dispersal limitations. Moreover, (iii) abundant taxa exhibit stronger temporal variability than rare taxa. They also found that abundant and rare sub-communities display similar spatial richness patterns that are negatively correlated with latitude and chlorophyll a and positively correlated with temperature. In summary, environmental filtering and dispersal limitations have different effects on abundant and rare picoeukaryotic sub-communities in different layers. Thus, depth appears as an essential variable that governs the structuring patterns of picoeukaryotic communities in the oceans and should be thoroughly considered to develop a more comprehensive understanding of oceanic microbial assemblages.