Disentangling the mechanisms shaping the surface ocean microbiota

Response of particle-attached and free-living bacterial communities to Microcystis blooms

The massive proliferation of Microcystis threatens freshwater ecosystems and degrades water quality globally. Understanding the mechanisms that contribute to Microcystis growth is crucial for managing Microcystis blooms. The lifestyles of bacteria can be classified generally into two groups: particle-attached (PA; > 3 µm) and free-living (FL; 0.2-3.0 µm). However, little is known about the response of PA and FL bacteria to Microcystis blooms. Using 16S rRNA gene high-throughput sequencing, we investigated the stability, assembly process, and co-occurrence patterns of PA and FL bacterial communities during distinct bloom stages. PA bacteria were phylogenetically different from their FL counterparts. Microcystis blooms substantially influenced bacterial communities. The time decay relationship model revealed that Microcystis blooms might increase the stability of both PA and FL bacterial communities. A contrasting community assembly mechanism was observed between the PA and FL bacterial communities. Throughout Microcystis blooms, homogeneous selection was the major assembly process that impacted the PA bacterial community, whereas drift explained much of the turnover of the FL bacterial community. Both PA and FL bacterial communities could be separated into modules related to different phases of Microcystis blooms. Microcystis blooms altered the assembly process of PA and FL bacterial communities. PA bacterial community appeared to be more responsive to Microcystis blooms than FL bacteria. Decomposition of Microcystis blooms may enhance cooperation among bacteria. Our findings highlight the importance of studying bacterial lifestyles to understand their functions in regulating Microcystis blooms. KEY POINTS: • Microcystis blooms alter the assembly process of PA and FL bacterial communities • Microcystis blooms increase the stability of both PA and FL bacterial communities • PA bacteria seem to be more responsive to Microcystis blooms than FL bacteria.

A fungi hotspot deep in the ocean: explaining the presence of Gjaerumia minor in equatorial Pacific bathypelagic waters

A plant parasite associated with the white haze disease in apples, the Basidiomycota Gjaerumia minor, has been found in most samples of the global bathypelagic ocean. An analysis of environmental 18S rDNA sequences on 12 vertical profiles of the Malaspina 2010 expedition shows that the relative abundance of this cultured species increases with depth while its distribution is remarkably different between the deep waters of the Pacific and Atlantic oceans, being present in higher concentrations in the former. This is evident from sequence analysis and a microscopic survey with a species-specific newly designed TSA-FISH probe. Several hints point to the hypothesis that G. minor is transported to the deep ocean attached to particles, and the absence of G. minor in bathypelagic Atlantic waters could then be explained by the absence of this organism in surface waters of the equatorial Atlantic. The good correlation of G. minor biomass with Apparent Oxygen Utilization, recalcitrant carbon and free-living prokaryotic biomass in South Pacific waters, together with the identification of the observed cells as yeasts and not as resting spores (teliospores), point to the possibility that once arrived at deep layers this species keeps on growing and thriving.

Response of bacterial and micro-eukaryotic communities to spatio-temporal fluctuations of wastewater in full scale constructed wetlands

How microbial communities respond to wastewater fluctuations is poorly understood. Full-scale surface flow constructed wetlands (SFCWs) were constructed for investigating microbial communities. Results showed that influent wastewater changed sediment bacterial community composition seasonally, indicating that a single bacterial taxonomic group had low resistance (especially, Actinobacteriota and Gammaproteobacteria). However, copy numbers of 16S rRNA, ammonia oxidizing archaea, ammonia oxidizing bacteria, nirS and nirK in the first stage SFCWs were 2.49 × 1010, 3.48 × 109, 5.76 × 106, 8.77 × 108 and 9.06 × 108 g-1 dry sediment, respectively, which remained stable between seasons. Moreover, decreases in the nitrogen concentration in wastewater, changed microbial system state from heterotrophic to autotrophic. Micro-eukaryotic communities were more sensitive to wastewater fluctuations than bacterial communities. Overall, results revealed that microbial communities responded to spatio-temporal fluctuations in wastewater through state changes and species asynchrony. This highlighted complex processes of wastewater treatment by microbial components in SFCWs.

A generic approach to infer community-level fitness of microbial genes

The gene content in a metagenomic pool defines the function potential of a microbial community. Natural selection, operating on the level of genomes or genes, shapes the evolution of community functions by enriching some genes while depriving the others. Despite the importance of microbiomes in the environment and health, a general metric to evaluate the community-wide fitness of microbial genes remains lacking. In this work, we adapt the classic neutral model of species and use it to predict how the abundances of different genes will be shaped by selection, regardless of at which level the selection acts. We establish a simple metric that quantitatively infers the average survival capability of each gene in a microbiome. We then experimentally validate the predictions using synthetic communities of barcoded Escherichia coli strains undergoing neutral assembly and competition. We further show that this approach can be applied to publicly available metagenomic datasets to gain insights into the environment-function interplay of natural microbiomes.

Wallace's line structures seagrass microbiota and is a potential barrier to the dispersal of marine bacteria

BACKGROUND

The processes that shape microbial biogeography are not well understood, and concepts that apply to macroorganisms, like dispersal barriers, may not affect microorganisms in the same predictable ways. To better understand how known macro-scale biogeographic processes can be applied at micro-scales, we examined seagrass associated microbiota on either side of Wallace's line to determine the influence of this cryptic dispersal boundary on the community structure of microorganisms. Communities were examined from twelve locations throughout Indonesia on either side of this theoretical line.

RESULTS

We found significant differences in microbial community structure on either side of this boundary (R2 = 0.09; P = 0.001), and identified seven microbial genera as differentially abundant on either side of the line, six of these were more abundant in the West, with the other more strongly associated with the East. Genera found to be differentially abundant had significantly smaller minimum cell dimensions (GLM: t923 = 59.50, P < 0.001) than the overall community.

CONCLUSION

Despite the assumed excellent dispersal ability of microbes, we were able to detect significant differences in community structure on either side of this cryptic biogeographic boundary. Samples from the two closest islands on opposite sides of the line, Bali and Komodo, were more different from each other than either was to its most distant island on the same side. We suggest that limited dispersal across this barrier coupled with habitat differences are primarily responsible for the patterns observed. The cryptic processes that drive macroorganism community divergence across this region may also play a role in the bigeographic patterns of microbiota.

Decoding populations in the ocean microbiome

Understanding the characteristics and structure of populations is fundamental to comprehending ecosystem processes and evolutionary adaptations. While the study of animal and plant populations has spanned a few centuries, microbial populations have been under scientific scrutiny for a considerably shorter period. In the ocean, analyzing the genetic composition of microbial populations and their adaptations to multiple niches can yield important insights into ecosystem function and the microbiome's response to global change. However, microbial populations have remained elusive to the scientific community due to the challenges associated with isolating microorganisms in the laboratory. Today, advancements in large-scale metagenomics and metatranscriptomics facilitate the investigation of populations from many uncultured microbial species directly from their habitats. The knowledge acquired thus far reveals substantial genetic diversity among various microbial species, showcasing distinct patterns of population differentiation and adaptations, and highlighting the significant role of selection in structuring populations. In the coming years, population genomics is expected to significantly increase our understanding of the architecture and functioning of the ocean microbiome, providing insights into its vulnerability or resilience in the face of ongoing global change. Video Abstract.

Spatiotemporal assembly and functional composition of planktonic microeukaryotic communities along productivity gradients in a subtropical lake

Microeukaryotes play crucial roles in the microbial loop of freshwater ecosystems, functioning both as primary producers and bacterivorous consumers. However, understanding the assembly of microeukaryotic communities and their functional composition in freshwater lake ecosystems across diverse environmental gradients remains limited. Here, we utilized amplicon sequencing of 18S rRNA gene and multivariate statistical analyses to examine the spatiotemporal and biogeographical patterns of microeukaryotes in water columns (at depths of 0.5, 5, and 10 m) within a subtropical lake in eastern China, covering a 40 km distance during spring and autumn of 2022. Our results revealed that complex and diverse microeukaryotic communities were dominated by Chlorophyta (mainly Chlorophyceae), Fungi, Alveolata, Stramenopiles, and Cryptophyta lineages. Species richness was higher in autumn than in spring, forming significant hump-shaped relationships with chlorophyll a concentration (Chl-a, an indicator of phytoplankton biomass). Microeukaryotic communities exhibited significant seasonality and distance-decay patterns. By contrast, the effect of vertical depth was negligible. Stochastic processes mainly influenced the assembly of microeukaryotic communities, explaining 63, 67, and 55% of community variation for spring, autumn, and both seasons combined, respectively. Trait-based functional analysis revealed the prevalence of heterotrophic and phototrophic microeukaryotic plankton with a trade-off along N:P ratio, Chl-a, and dissolved oxygen (DO) gradients. Similarly, the mixotrophic proportions were significantly and positively correlated with Chl-a and DO concentrations. Overall, our findings may provide useful insights into the assembly patterns of microeukaryotes in lake ecosystem and how their functions respond to environmental changes.

Microbial community and network responses across strong environmental gradients: How do they compare with macroorganisms?

The way strong environmental gradients shape multispecific assemblages has allowed us to examine a suite of ecological and evolutionary hypotheses about structure, regulation and community responses to fluctuating environments. But whether the highly diverse co-occurring microorganisms are shaped in similar ways as macroscopic organisms across the same gradients has yet to be addressed in most ecosystems. Here, we characterize intertidal biofilm bacteria communities, comparing zonation at both the "species" and community levels, as well as network attributes, with co-occurring macroalgae and invertebrates in the same rocky shore system. The results revealed that the desiccation gradient has a more significant impact on smaller communities, while both desiccation and submersion gradients (surge) affect the larger, macroscopic communities. At the community level, we also confirmed the existence of distinct communities within each intertidal zone for microorganisms, similar to what has been previously described for macroorganisms. But our results indicated that dominant microbial organisms along the same environmental gradient exhibited less differentiation across tidal levels than their macroscopic counterparts. However, despite the substantial differences in richness, size and attributes of co-occurrence networks, both macro- and micro-communities respond to stress gradients, leading to the formation of similar zonation patterns in the intertidal rocky shore.

Subtropical coastal microbiome variations due to massive river runoff after a cyclonic event

BACKGROUND

Coastal ecosystem variability at tropical latitudes is dependent on climatic conditions. During the wet, rainy season, extreme climatic events such as cyclones, precipitation, and winds can be intense over a short period and may have a significant impact on the entire land‒sea continuum. This study focused on the effect of river runoff across the southwest coral lagoon ecosystem of Grand Terre Island of New Caledonia (South Pacific) after a cyclonic event, which is considered a pulse disturbance at our study site. The variability of coastal microbiomes, studied by the metabarcoding of V4 18S (protists) and V4-V5 16S (bacteria) rDNA genes, after the cyclone passage was associated with key environmental parameters describing the runoff impact (salinity, organic matter proxies, terrestrial rock origin metals) and compared to community structures observed during the dry season.

RESULTS

Microbiome biodiversity patterns of the dry season were destructured because of the runoff impact, and land-origin taxa were observed in the coastal areas. After the rainy event, different daily community dynamics were observed locally, with specific microbial taxa explaining these variabilities. Plume dispersal modeling revealed the extent of low salinity areas up to the coral reef area (16 km offshore), but a rapid (< 6 days) recovery to typical steady conditions of the lagoon's hydrology was observed. Conversely, during the same time, some biological components (microbial communities, Chl a) and biogeochemical components (particulate nickel, terrigenous organic matter) of the ecosystem did not recover to values observed during the dry season conditions.

CONCLUSION

The ecosystem resilience of subtropical ecosystems must be evaluated from a multidisciplinary, holistic perspective and over the long term. This allows evaluating the risk associated with a potential continued and long-term disequilibrium of the ecosystem, triggered by the change in the frequency and intensity of extreme climatic events in the era of planetary climatic changes.

Discrepant assembly processes of prokaryotic communities between the abyssal and hadal sediments in Yap Trench

Abyssal and hadal sediments represent two of the most type ecosystems on Earth and have the potential interactions with geochemistry. However, little is known about the prokaryotic community assembly and the response of prokaryotic communities to metal(loid)s in trench sediments due to the lack of adequate and appropriate samples. In this study, a systematic investigation combined the assembly mechanisms and co-occurrence patterns of prokaryotic communities between the hadal and abyssal sediments across the Yap Trench. The results revealed that the hadal prokaryotes had less species diversity, but more abundant function than the abyssal prokaryotes. The prokaryotic communities in the abyssal sediments had more core taxa than the hadal sediments. Twenty-one biomarkers mostly affiliated with Nitrosopumilaceae were detected using Random-Forests machine learning algorithm. Furthermore, stochasticity was dominant in the prokaryotic community assembly processes of the Yap Trench sediments. Meanwhile, homogeneous selection (32.6%-52.9%) belonging to deterministic processes governed the prokaryotic community assembly in hadal sediments with increasing of sediment depth. In addition to total nitrogen and total organic carbon, more metal(loid)s were significantly correlated with the prokaryotic community in the hadal sediments than that in the abyssal sediments. The hadal prokaryotic communities was most positively related to bismuth (r = 0.31, p < 0.01), followed by calcium, chromium, cerium, potassium, plumbum, scandium, titanium, and vanadium. Finally, co-occurrence networks revealed two potential dominant prokaryotic modules in Yap Trench sediments covaried across oceanographic zonation. By contrast, the hadal network had relatively more complexity, more bacterial taxa, and more associations among prokaryotic taxa, relative to the abyssal network. This study reveals potentially metal variables and community assembly mechanisms of the prokaryotic community in abyssal and hadal sediments and provides a better understanding on the prokaryotic diversity and ecology in trench sediment ecosystems.

Different Microeukaryotic Trophic Groups Show Different Latitudinal Spatial Scale Dependences in Assembly Processes across the Continental Shelves of China

The relative role of stochasticity versus determinism is critically dependent on the spatial scale over which communities are studied. However, only a few studies have attempted to reveal how spatial scales influence the balance of different assembly processes. In this study, we investigated the latitudinal spatial scale dependences in assembly processes of microeukaryotic communities in surface water and sediment along the continental shelves of China. It was hypothesized that different microeukaryotic trophic groups (i.e., autotroph, heterotroph, mixotroph, and parasite) showed different latitudinal scale dependences in their assembly processes. Our results disclosed that the relative importance of different assembly processes depended on a latitudinal space scale for planktonic microeukaryotes. In surface water, as latitudinal difference increased, the relative contributions of homogenous selection and homogenizing dispersal decreased for the entire community, while those of heterogeneous selection and drift increased. The planktonic autotrophic and heterotrophic groups shifted from stochasticity-dominated processes to heterogeneous selection as latitudinal differences surpassed thresholds of 8° and 16°, respectively. For mixotrophic and parasitic groups, however, the assembly processes were always dominated by drift across different spatial scales. The balance of different assembly processes for the autotrophic group was mainly driven by temperature, whereas that of the heterotrophic group was driven by salinity and geographical distance. In sediment, neither the entire microeukaryotic community nor the four trophic groups showed remarkable spatial scale dependences in assembly processes; they were always overwhelmingly dominated by the drift. This work provides a deeper understanding of the distribution mechanisms of microeukaryotes along the continental shelves of China from the perspective of trophic groups.

Disentangling microbial networks across pelagic zones in the tropical and subtropical global ocean

Microbial interactions are vital in maintaining ocean ecosystem function, yet their dynamic nature and complexity remain largely unexplored. Here, we use association networks to investigate possible ecological interactions in the marine microbiome among archaea, bacteria, and picoeukaryotes throughout different depths and geographical regions of the tropical and subtropical global ocean. Our findings reveal that potential microbial interactions change with depth and geographical scale, exhibiting highly heterogeneous distributions. A few potential interactions were global, meaning they occurred across regions at the same depth, while 11-36% were regional within specific depths. The bathypelagic zone had the lowest proportion of global associations, and regional associations increased with depth. Moreover, we observed that most surface water associations do not persist in deeper ocean layers despite microbial vertical dispersal. Our work contributes to a deeper understanding of the tropical and subtropical global ocean interactome, which is essential for addressing the challenges posed by global change.

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

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

Global Distribution and Diversity of Marine Parmales

Parmales (Bolidophyceae) is a minor eukaryotic phytoplankton group, sister to diatoms, which exists as two distinct forms of unicellular organisms: silicified cells and naked flagellates. Since their discovery, many field studies on Parmales have been performed; however, their global distribution has not yet been examined in detail. We herein compiled more than 3,000 marine DNA metabarcoding datasets targeting the V4 region of the 18S rRNA gene from the EukBank database. By linking this large dataset with the latest morphological and genetic information, we provide updated estimates on the diversity and distribution of Parmales in the global ocean at a fine taxonomic resolution. Parmalean amplicon sequence variants (ASVs) were detected in nearly 90% of the samples analyzed. However, the relative abundance of parmaleans in the eukaryotic community was less than 0.2% on average, and the estimated true richness of parmalean ASVs was approximately 316 ASVs, confirming their low abundance and diversity. A phylogenetic ana-lysis divided these algae into four clades, and three known morphotypes of silicified cells were classified into three different clades. The abundance of Parmales is generally high in the poles and decreases towards the tropics, and individual clades/subclades show further distinctions in their distribution. Collectively, the present results suggest clade/subclade-specific adaptation to different ecological niches.

Contrasting drivers of abundant phage and prokaryotic communities revealed in diverse coastal ecosystems

Phages (viruses of bacteria and archaea) are a ubiquitous top-down control on microbial communities by selectively infecting and killing cells. As obligate parasites, phages are inherently linked to processes that impact their hosts' distribution and physiology, but phages can also be impacted by external, environmental factors, such as UV radiation degrading their virions. To better understand these complex links of phages to their hosts and the environment, we leverage the unique ecological context of the Isthmus of Panama, which narrowly disconnects the productive Tropical Eastern Pacific (EP) and nutrient-poor Tropical Western Atlantic (WA) provinces. We could thus compare patterns of phage and prokaryotic communities at both global scales (between oceans) and local-scales (between habitats within an ocean). Although both phage and prokaryotic communities differed sharply between the oceans, phage community composition did not significantly differ between mangroves and reefs of the WA, while prokaryotic communities were distinct. These results suggest phages are more shaped by dispersal processes than local conditions regardless of spatial scale, while prokaryotes tend to be shaped by local conditions at smaller spatial scales. Collectively, we provide a framework for addressing the co-variability between phages and prokaryotes in marine systems and identifying factors that drive consistent versus disparate trends in community shifts, essential to informing models of biogeochemical cycles that include these interactions.

Spatial metatranscriptomics resolves host-bacteria-fungi interactomes

The interactions of microorganisms among themselves and with their multicellular host take place at the microscale, forming complex networks and spatial patterns. Existing technology does not allow the simultaneous investigation of spatial interactions between a host and the multitude of its colonizing microorganisms, which limits our understanding of host-microorganism interactions within a plant or animal tissue. Here we present spatial metatranscriptomics (SmT), a sequencing-based approach that leverages 16S/18S/ITS/poly-d(T) multimodal arrays for simultaneous host transcriptome- and microbiome-wide characterization of tissues at 55-µm resolution. We showcase SmT in outdoor-grown Arabidopsis thaliana leaves as a model system, and find tissue-scale bacterial and fungal hotspots. By network analysis, we study inter- and intrakingdom spatial interactions among microorganisms, as well as the host response to microbial hotspots. SmT provides an approach for answering fundamental questions on host-microbiome interplay.

Global biogeography of the smallest plankton across ocean depths

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

Molecular and morphological characterization of four new ancyromonad genera and proposal for an updated taxonomy of the Ancyromonadida

Ancyromonads are small biflagellated protists with a bean-shaped morphology. They are cosmopolitan in marine, freshwater, and soil environments, where they attach to surfaces while feeding on bacteria. These poorly known grazers stand out by their uncertain phylogenetic position in the tree of eukaryotes, forming a deep-branching "orphan" lineage that is considered key to a better understanding of the early evolution of eukaryotes. Despite their ecological and evolutionary interest, only limited knowledge exists about their true diversity. Here, we aimed to characterize ancyromonads better by integrating environmental surveys with behavioral observation and description of cell morphology, for which sample isolation and culturing are indispensable. We studied 18 ancyromonad strains, including 14 new isolates and seven new species. We described three new and genetically divergent genera: Caraotamonas, Nyramonas, and Olneymonas, together encompassing four species. The remaining three new species belong to the already-known genera Fabomonas and Ancyromonas. We also raised Striomonas, formerly a subgenus of Nutomonas, to full genus status, on morphological and phylogenetic grounds. We studied the morphology of diverse ancyromonads under light and electron microscopy and carried out molecular phylogenetic analyses, also including 18S rRNA gene sequences from several environmental surveys. Based on these analyses, we have updated the taxonomy of Ancyromonadida.

The Bay of Bengal exposes abundant photosynthetic picoplankton and newfound diversity along salinity-driven gradients

The Bay of Bengal (BoB) is a 2,600,000 km2 expanse in the Indian Ocean upon which many humans rely. However, the primary producers underpinning food chains here remain poorly characterized. We examined phytoplankton abundance and diversity along strong BoB latitudinal and vertical salinity gradients-which have low temperature variation (27-29°C) between the surface and subsurface chlorophyll maximum (SCM). In surface waters, Prochlorococcus averaged 11.7 ± 4.4 × 104  cells ml-1 , predominantly HLII, whereas LLII and 'rare' ecotypes, HLVI and LLVII, dominated in the SCM. Synechococcus averaged 8.4 ± 2.3 × 104  cells ml-1 in the surface, declined rapidly with depth, and population structure of dominant Clade II differed between surface and SCM; Clade X was notable at both depths. Across all sites, Ostreococcus Clade OII dominated SCM eukaryotes whereas communities differentiated strongly moving from Arabian Sea-influenced high salinity (southerly; prasinophytes) to freshwater-influenced low salinity (northerly; stramenopiles, specifically, diatoms, pelagophytes, and dictyochophytes, plus the prasinophyte Micromonas) surface waters. Eukaryotic phytoplankton peaked in the south (1.9 × 104  cells ml-1 , surface) where a novel Ostreococcus was revealed, named here Ostreococcus bengalensis. We expose dominance of a single picoeukaryote and hitherto 'rare' picocyanobacteria at depth in this complex ecosystem where studies suggest picoplankton are replacing larger phytoplankton due to climate change.

Soil Micro-eukaryotic Diversity Patterns Along Elevation Gradient Are Best Estimated by Increasing the Number of Elevation Steps Rather than Within Elevation Band Replication

The development of high-throughput sequencing (HTS) of environmental DNA (eDNA) has stimulated the study of soil microbial diversity patterns and drivers at all scales. However, given the heterogeneity of soils, a challenge is to define effective and efficient sampling protocols that allow sound comparison with other records, especially vegetation. In studies of elevational diversity pattern, a trade-off is choosing between replication within elevation bands vs. sampling more elevation bands. We addressed this question for soil protists along an elevation gradient on Mt. Asahi, Hokkaido, Japan. We compared two sampling approaches: (1) the replicate strategy (five replicates at six elevational bands, total = 30) and (2) the transect strategy (one sample in each of 16 different elevational bands). Despite a nearly twofold lower sampling effort, the transect strategy yielded congruent results compared to the replicate strategy for the estimation of elevational alpha diversity pattern: the regression coefficients between diversity indices and elevation did not differ between the two options. Furthermore, for a given total number of samples, gamma diversity estimated across the entire transect was higher when sampling more elevational bands as compared to replication from fewer elevational bands. Beta diversity (community composition turnover) was lower within a given elevational band than between adjacent bands and increased with elevation distance. In redundancy analyses, soil organic matter-related variable (the first principal component of soil organic matter, water content, total organic carbon, and nitrogen by whom were highly correlated) and elevation best explained elevational beta diversity pattern for both sampling approaches. Taken together, our results suggest that sampling a single plot per elevation band will be sufficient to obtain a good estimate of soil micro-eukaryotic diversity patterns along elevation gradients. This study demonstrated the effectiveness of the transect strategy in estimating diversity patterns along elevation gradients which is instructive for future environmental or even experimental studies. While not advocating for completely replacing replication-based sampling practices, it is important to note that both replicate and transect strategies have their merits and can be employed based on specific research goals and resource limitations.

Mechanisms shaping dissolved organic matter and microbial community in lake ecosystems

Lakes are active components of the global carbon cycle and host a range of processes that degrade and modify dissolved organic matter (DOM). Through the degradation of DOM molecules and the synthesis of new compounds, microbes in aquatic environments strongly and continuously influence chemodiversity, which can feedback to influence microbial diversity. Developing a better understanding of the biodiversity patterns that emerge along spatial and environmental gradients is one of the key objectives of community ecology. A changing climate may affect ecological feedback, including those that affect microbial communities. To maintain the function of a lake ecosystem and predict carbon cycling in the environment, it is increasingly important to understand the coupling between microbial and DOM diversity. To unravel the biotic and abiotic mechanisms that control the structure and patterns of DOM and microbial communities in lakes, we combined high-throughput sequencing and ultra-high resolution mass spectrometry together with a null modeling approach. The advantage of null models is their ability to evaluate the relative influences of stochastic and deterministic assembly processes in both DOM and microbial community assemblages. The present study includes spatiotemporal signatures of DOM and the microbial community in six temperate lakes contrasting continental and Mediterranean climates during the productive season. Different environmental conditions and nutrient sources characterized the studied lakes. Our results have shown high covariance between molecular-level DOM diversity and the diversity of individual microbial communities especially with diversity of microeukaryotes and free-living bacteria indicating their dynamic feedback. We found that the differences between lakes and climatic regions were mainly reflected in the diversity of DOM at the molecular formula-level and the microeukaryota community. Furthermore, using null models the DOM assembly was governed by deterministic variable selection operating consistently and strongly within and among lakes. In contrast, microbial community assembly processes were highly variable across lakes with different trophic status and climatic regions. Difference in the processes governing DOM and microbial composition does not indicate weak coupling between these components, rather it suggests that distinct factors may be influencing microbial communities and DOM assemblages separately. Further understanding of the DOM-microbe coupling (or lack thereof) is key to formulating predictive models of future lake ecology and function.

Ecological mechanisms and current systems shape the modular structure of the global oceans' prokaryotic seascape

Major biogeographic features of the microbial seascape in the oceans have been established and their underlying ecological mechanisms in the (sub)tropical oceans and the Pacific Ocean identified. However, we still lack a unifying understanding of how prokaryotic communities and biogeographic patterns are affected by large-scale current systems in distinct ocean basins and how they are globally shaped in line with ecological mechanisms. Here we show that prokaryotic communities in the epipelagic Pacific and Atlantic Ocean, in the southern Indian Ocean, and the Mediterranean Sea are composed of modules of co-occurring taxa with similar environmental preferences. The relative partitioning of these modules varies along latitudinal and longitudinal gradients and are related to different hydrographic and biotic conditions. Homogeneous selection and dispersal limitation were identified as the major ecological mechanisms shaping these communities and their free-living (FL) and particle-associated (PA) fractions. Large-scale current systems govern the dispersal of prokaryotic modules leading to the highest diversity near subtropical fronts.

A Metagenomic and Amplicon Sequencing Combined Approach Reveals the Best Primers to Study Marine Aerobic Anoxygenic Phototrophs

Studies based on protein-coding genes are essential to describe the diversity within bacterial functional groups. In the case of aerobic anoxygenic phototrophic (AAP) bacteria, the pufM gene has been established as the genetic marker for this particular functional group, although available primers are known to have amplification biases. We review here the existing primers for pufM gene amplification, design new ones, and evaluate their phylogenetic coverage. We then use samples from contrasting marine environments to evaluate their performance. By comparing the taxonomic composition of communities retrieved with metagenomics and with different amplicon approaches, we show that the commonly used PCR primers are biased towards the Gammaproteobacteria phylum and some Alphaproteobacteria clades. The metagenomic approach, as well as the use of other combinations of the existing and newly designed primers, show that these groups are in fact less abundant than previously observed, and that a great proportion of pufM sequences are affiliated to uncultured representatives, particularly in the open ocean. Altogether, the framework developed here becomes a better alternative for future studies based on the pufM gene and, additionally, serves as a reference for primer evaluation of other functional genes.

High Microeukaryotic Diversity in the Cold-Seep Sediment

Microeukaryotic diversity, community structure, and their regulating mechanisms remain largely unclear in chemosynthetic ecosystems. Here, using high-throughput sequencing data of 18S rRNA genes, we explored microeukaryotic communities from the Haima cold seep in the northern South China Sea. We compared three distinct habitats: active, less active, and non-seep regions, with vertical layers (0-25 cm) from sediment cores. The results showed that seep regions harbored more abundant and diverse parasitic microeukaryotes (e.g., Apicomplexa and Syndiniales) as indicator species, compared to nearby non-seep region. Microeukaryotic community heterogeneity was larger between habitats than within habitat, and greatly increased when considering molecular phylogeny, suggesting the local diversification in cold-seep sediments. Microeukaryotic α-diversity at cold seeps was positively increased by metazoan richness and dispersal rate of microeukaryotes, while its β-diversity was promoted by heterogeneous selection mainly from metazoan communities (as potential hosts). Their combined effects led to the significant higher γ-diversity (i.e., total diversity in a region) at cold seeps than non-seep regions, suggesting cold-seep sediment as a hotspot for microeukaryotic diversity. Our study highlights the importance of microeukaryotic parasitism in cold-seep sediment and has implications for the roles of cold seep in maintaining and promoting marine biodiversity.

Contrasting Mechanisms Determine the Microeukaryotic and Syndiniales Community Assembly in a Eutrophic bay

Syndiniales is a diverse parasitic group, increasingly gaining attention owing to its high taxonomic diversity in marine ecosystems and inhibitory effects on the dinoflagellate blooms. However, their seasonal dynamics, host interactions, and mechanisms of community assembly are largely unknown, particularly in eutrophic waters. Here, using 18S rRNA gene amplicon sequencing, we intended to elucidate the interactions between Syndiniales and microeukaryotes, as well as community assembly processes in a eutrophic bay. The results showed that Syndiniales group II was dominating throughout the year, with substantially higher abundance in the winter and spring, whereas Syndiniales group I was more abundant in the summer and autumn. Temperature and Dinoflagellata were the most important abiotic and biotic factors driving variations of the Syndiniales community, respectively. The assembly processes of microeukaryotes and Syndiniales were completely different, with the former being controlled by a balance between homogeneous selection and drift and the latter being solely governed by drift. Network analysis revealed that Syndiniales group II had the largest number of interactions with microeukaryotes, and they primarily associated with Dinoflagellata in the winter, while interactions with Chlorophyta and Bacillariophyta increased dramatically in summer and autumn. These findings provide significant insights in understanding the interactions and assembly processes of Syndiniales throughout the year, which is critical in revealing the roles of single-celled parasites in driving protist dynamics in eutrophic waters.

Development and validation of a real-time PCR assay protocol for the specific detection and quantification of pelagiphages in seawater samples

Earth is inhabited by numerous adaptations of cellular forms shaped by the persistent scrutiny of natural selection. Thus, as natural selection has fixed beneficial adaptations of functional traits, cellular life has conquered almost all environmental niches on our planet. However, cellular life succumbs in number and genetic diversity to viruses. Among all viruses, phages are highly prevalent in diverse environments, and due to their vast genetic diversity and abundance, their relevant role as significant players in several ecological processes is now fully recognized. Pelagiphages, bacteriophages infecting bacteria of the SAR11 clade, are the most abundant viruses in the oceans. However, the ecological contribution of pelagiphages on populations of Pelagibacterales remains largely underestimated. An essential aspect of estimating the impact of bacteriophages is their absolute and precise quantification, which provides relevant information about the host-virus interactions and the structure of viral assemblages. Consequently, due to its abundance and claimed influence in the biogeochemical cycling of elements, the accurate quantification of pelagiphages results in an essential task. This study describes the development and validation of a sensitive, specific, accurate and reproducible qPCR platform targeting pelagiphages. Moreover, this method allowed the detection and quantification of pelagiphages in the Gulf of California for the first time.

Predicting global distributions of eukaryotic plankton communities from satellite data

Satellite remote sensing is a powerful tool to monitor the global dynamics of marine plankton. Previous research has focused on developing models to predict the size or taxonomic groups of phytoplankton. Here, we present an approach to identify community types from a global plankton network that includes phytoplankton and heterotrophic protists and to predict their biogeography using global satellite observations. Six plankton community types were identified from a co-occurrence network inferred using a novel rDNA 18 S V4 planetary-scale eukaryotic metabarcoding dataset. Machine learning techniques were then applied to construct a model that predicted these community types from satellite data. The model showed an overall 67% accuracy in the prediction of the community types. The prediction using 17 satellite-derived parameters showed better performance than that using only temperature and/or the concentration of chlorophyll a. The constructed model predicted the global spatiotemporal distribution of community types over 19 years. The predicted distributions exhibited strong seasonal changes in community types in the subarctic-subtropical boundary regions, which were consistent with previous field observations. The model also identified the long-term trends in the distribution of community types, which suggested responses to ocean warming.

Spatiotemporal dynamics of high and low nucleic acid-content bacterial communities in Chinese coastal seawater: assembly process, co-occurrence relationship and the ecological functions

Studies of high nucleic acid-content (HNA) and low nucleic acid-content (LNA) bacterial communities are updating our view of their distributions and taxonomic composition. However, there are still large gaps in our knowledge of the composition, assembly processes, co-occurrence relationships and ecological functions of HNA and LNA bacterial communities. Here, using 16S rRNA gene amplicon sequencing, we investigated the spatiotemporal dynamics, assembly processes, co-occurrence relationships and ecological functions of HNA and LNA bacterial communities in the samples collected in summer and winter in Chinese coastal seas. The communities of HNA and LNA bacteria had clear spatiotemporal patterns and LNA bacteria was phylogenetically less diverse than HNA bacteria in both seasons. The distribution of HNA and LNA bacteria were significantly affected by the environmental factors and a significant seasonal-consistent distance-decay patterns were found in HNA and LNA bacteria. Furthermore, a quantitative assessment of ecological processes revealed that dispersal limitation, homogeneous selection exerted important roles in the community assembly of HNA and LNA bacteria. More importantly, we observed seasonality in the co-occurrence relationships: closer inter-taxa connections of HNA bacterial communities in winter than in summer and the opposite is true in the LNA bacterial communities. Some ecological functions, such as: phototrophy, photoautotrophy, oxygenic photoautotrophy, were different between HNA and LNA bacteria. These results provide a better understanding of spatiotemporal patterns, processes, and the ecological functions of HNA and LNA bacterial communities in Chinese coastal seawater.

Identification of the heterotrophic nanoflagellate Bilabrum latius in the southern Adriatic (Mediterranean Sea)

Heterotrophic flagellates (HF) represent an important protist group in marine ecosystem functioning. Characterised by high taxonomic diversity, identification and classification of HF is often difficult using classical methods of light microscopy (LM). Complementing LM with molecular methods, such as barcoding, enables reliable taxonomic identification of even small size nanoflagellates that share similar or unnoticeable morphological features. The order Bicosoecida is a group of heterotrophic nanoflagellates that are important part of protist plankton and benthic communities of the world oceans. Recently, on the basis of high-resolution light microscopy and barcoding, a new bicosoecid genus, Bilabrum, was described with B. latius sp. as a type species. Our study reports on identification of B. latius from co-culture with the diatom species Chaetoceros affinis isolated from fresh plankton samples collected in the southern Adriatic. This detection of the Adriatic B.latius represents first record of this species outside itś up to now known and described habitat (deep-sea sediment of the South - East Atlantic Ocean) and in diatom co-culture.

Assembly processes and functional diversity of marine protists and their rare biosphere

BACKGROUND

The mechanisms shaping the rare microbial biosphere and its role in ecosystems remain unclear. We developed an approach to study ecological patterns in the rare biosphere and use it on a vast collection of marine microbiomes, sampled in coastal ecosystems at a regional scale. We study the assembly processes, and the ecological strategies constituting the rare protistan biosphere. Using the phylogeny and morpho-trophic traits of these protists, we also explore their functional potential.

RESULTS

Taxonomic community composition remained stable along rank abundance curves. Conditionally rare taxa, driven by selection processes, and transiently rare taxa, with stochastic distributions, were evidenced along the rank abundance curves of all size-fractions. Specific taxa within the divisions Sagenista, Picozoa, Telonemia, and Choanoflagellida were rare across time and space. The distribution of traits along rank abundance curves outlined a high functional redundancy between rare and abundant protists. Nevertheless, trophic traits illustrated an interplay between the trophic groups of different size-fractions.

CONCLUSIONS

Our results suggest that rare and abundant protists are evolutionary closely related, most notably due to the high microdiversity found in the rare biosphere. We evidenced a succession of assembly processes and strategies of rarity along rank abundance curves that we hypothesize to be common to most microbiomes at the regional scale. Despite high functional redundancy in the rare protistan biosphere, permanently rare protists were evidenced, and they could play critical functions as bacterivores and decomposers from within the rare biosphere. Finally, changes in the composition of the rare protistan biosphere could be influenced by the trophic regime of aquatic ecosystems. Our work contributes to understanding the role of rare protists in microbiomes.

Biogeographical and Biodiversity Patterns of Marine Planktonic Bacteria Spanning from the South China Sea across the Gulf of Bengal to the Northern Arabian Sea

Understanding the biogeographical and biodiversity patterns of bacterial communities is essential in unraveling their responses to future environmental changes. However, the relationships between marine planktonic bacterial biodiversity and seawater chlorophyll a are largely understudied. Here, we used high-throughput sequencing to study the biodiversity patterns of marine planktonic bacteria across a broad chlorophyll a gradient spanning from the South China Sea across the Gulf of Bengal to the northern Arabian Sea. We found that the biogeographical patterns of marine planktonic bacteria complied with the scenario of homogeneous selection, with chlorophyll a concentration being the key environmental selecting variable of bacteria taxa. The relative abundance of Prochlorococcus, the SAR11 clade, the SAR116 clade, and the SAR86 clade significantly decreased in habitats with high chlorophyll a concentrations (>0.5 μg/L). Free-living bacteria (FLB) and particle-associated bacteria (PAB) displayed contrasting alpha diversity and chlorophyll a relationships with a positive linear correlation for FLB but a negative correlation for PAB. We further found that PAB had a narrower niche breadth of chlorophyll a than did FLB, with far fewer bacterial taxa being favored at higher chlorophyll a concentrations. Higher chlorophyll a concentrations were linked to the enhanced stochastic drift and reduced beta diversity of PAB but to the weakened homogeneous selection, enhanced dispersal limitation, and increased beta diversity of FLB. Taken together, our findings might broaden our knowledge about the biogeography of marine planktonic bacteria and advance the understanding of bacterial roles in predicting ecosystem functioning under future environmental changes that are derived from eutrophication. IMPORTANCE One of the long-standing interests of biogeography is to explore diversity patterns and uncover their underlying mechanisms. Despite intensive studies on the responses of eukaryotic communities to chlorophyll a concentrations, we know little about how changes in seawater chlorophyll a concentrations affect free-living bacteria (FLB) and particle-associated bacteria (PAB) diversity patterns in natural systems. Our biogeography study demonstrated that marine FLB and PAB displayed contrasting diversity and chlorophyll a relationships and exhibited completely different assembly mechanisms. Our findings broaden our knowledge about the biogeographical and biodiversity patterns of marine planktonic bacteria in nature systems and suggest that PAB and FLB should be considered independently in predicting marine ecosystem functioning under future frequent eutrophication.

Diversity and interaction of bacterial and microeukaryotic communities in sediments planted with different submerged macrophytes: Responses to decabromodiphenyl ether

Although various persistent organic pollutants (POPs) can affect microbial communities and functions in aquatic ecosystems, little is known about how bacteria and microeukaryotes respond to the POPs in sediments planted with different submerged macrophytes. Here, a 60-day microcosm experiment was carried out to investigate the changes in the diversity and interaction of bacterial and microeukaryotic communities in sediments collected from Taihu lake, either with decabromodiphenyl ether (BDE-209) own or combined with two common submerged macrophyte species (Vallisneria natans and Hydrilla verticillate). The results showed that BDE-209 significantly decreased the bacterial α-diversity but increased the microeukaryotic one. In sediments planted with submerged macrophytes, the negative effect of BDE-209 on bacterial diversity was weakened, and its positive effect on microeukaryotic one was strengthened. Co-occurrence network analysis revealed that the negative relationship was dominant in bacterial and microeukaryotic communities, while the cooperative relationship between microbial species was increased in planted sediments. Among nine keystone species, one belonging to bacterial family Thermoanaerobaculaceae was enriched by BDE-209, and others were inhibited. Notably, such inhibition was weakened, and the stimulation was enhanced in planted sediments. Together, these observations indicate that the responses of bacteria and microeukaryotes to BDE-209 are different, and their communities under BDE-209 contamination are more stable in sediments planted with submerged macrophytes. Moreover, the effects of plant species on the microbial responses to BDE-209 need to be explored by more specific field studies in the future.

A meta-analysis reveals that geographical factors drive the bacterial community variation in Chinese lakes

The biogeographical distribution of plants and animals has been extensively studied, however, the biogeographical patterns and the factors that influence bacterial communities in lakes over large scales are yet to be fully understood, even though they play critical roles in biogeochemical cycles. Here, bacterial community compositional data, geographic information, and environmental factors were integrated for 326 Chinese lakes based on previously published studies to determine the underlying factors that shape bacterial diversity among Chinese lakes. The composition of bacterial communities significantly varied among the three primary climatic regions of China (Northern China, NC; Southern China, SC; and the Tibetan Plateau, TIP), and across two different lake habitats (waters and sediments). Sediment bacterial communities exhibited significantly higher alpha-diversity and distance-decay relationships compared to water communities. The results indicate that the "scale-dependent patterns" of controlling factors, primarily influenced by geographical factors, become increasingly pronounced as the spatial scale increases. At a national scale, geographical factors exerted a dominant influence on both the water and sediment communities across all lakes, as geographical barriers restrict the dispersal of individuals. At smaller spatial scales, temperature-driven selection effects played a greater role in shaping water bacterial community variation in the NC, SC, and TIP, while geographical factors had a stronger association with sediment bacterial community variation in the lakes of the three regions. This synthesis offers novel insights into the ecological factors that determine the distribution of bacteria in Chinese lakes.

Disentangling temporal associations in marine microbial networks

BACKGROUND

Microbial interactions are fundamental for Earth's ecosystem functioning and biogeochemical cycling. Nevertheless, they are challenging to identify and remain barely known. Omics-based censuses are helpful in predicting microbial interactions through the statistical inference of single (static) association networks. Yet, microbial interactions are dynamic and we have limited knowledge of how they change over time. Here, we investigate the dynamics of microbial associations in a 10-year marine time series in the Mediterranean Sea using an approach inferring a time-resolved (temporal) network from a single static network.

RESULTS

A single static network including microbial eukaryotes and bacteria was built using metabarcoding data derived from 120 monthly samples. For the decade, we aimed to identify persistent, seasonal, and temporary microbial associations by determining a temporal network that captures the interactome of each individual sample. We found that the temporal network appears to follow an annual cycle, collapsing, and reassembling when transiting between colder and warmer waters. We observed higher association repeatability in colder than in warmer months. Only 16 associations could be validated using observations reported in literature, underlining our knowledge gap in marine microbial ecological interactions.

CONCLUSIONS

Our results indicate that marine microbial associations follow recurrent temporal dynamics in temperate zones, which need to be accounted for to better understand the functioning of the ocean microbiome. The constructed marine temporal network may serve as a resource for testing season-specific microbial interaction hypotheses. The applied approach can be transferred to microbiome studies in other ecosystems. Video Abstract.

Community assemblages and species coexistence of prokaryotes controlled by local environmental heterogeneity in a cold seep water column

The distribution and heterogeneity characteristics of microbial communities in cold seep water columns are significant factors governing the efficiency of methane filtering and carbon turnover. However, this process is poorly understood. The diversity of vertically stratified microbial communities and the factors controlling the community assemblage process in the water column above the Haima cold seep were investigated in this study. The prokaryotic community diversities varied distinctly with vertical changes in hydrochemistry. Cyanobacteria dominated the light-transmitting layers and Proteobacteria dominated the deeper layers. With respect to microbial community assemblages and co-occurrence networks, stochastic processes were particularly important in shaping prokaryotic communities. In the shallow (≥85 m) and mesopelagic water columns (600-800 m), microbial community characteristics were affected by deterministic processes, reduced network connectivity, and modularity. Microbial community diversities and assemblage processes along a vertical profile were influenced by the vertical variations in pH, temperature, DIC, and nutrients. Stochastic processes may have facilitated the formation of complex co-occurrence networks. Briefly, the distribution of local environmental heterogeneity along the vertical dimension could drive unique microbial community assemblage and species coexistence patterns. This study provides new perspectives on how microorganisms adapt to the environment and build communities, and how species coexist in shared habitats.

Distinct Ecological Processes Mediate Domain-Level Differentiation in Microbial Spatial Scaling

The spatial scaling of biodiversity, such as the taxa-area relationship (TAR) and distance-decay relationship (DDR), is a typical ecological pattern that is followed by both microbes and macrobes in natural ecosystems. Previous studies focusing on microbes mainly aimed to address whether and how different types of microbial taxa differ in spatial scaling patterns, leaving the underlying mechanisms largely untouched. In this study, the spatial scaling of different microbial domains and their associated ecological processes in an intertidal zone were comparatively investigated. The significant spatial scaling of biodiversity could be observed across all microbial domains, including archaea, bacteria, fungi, and protists. Among them, archaea and fungi were found with much stronger DDR slopes than those observed in bacteria and protists. For both TAR and DDR, rare subcommunities were mainly responsible for the observed spatial scaling patterns, except for the DDR of protists and bacteria. This was also evidenced by extending the TAR and DDR diversity metrics to Hill numbers. Further statistical analyses demonstrated that different microbial domains were influenced by different environmental factors and harbored distinct local community assembly processes. Of these, drift was mainly responsible for the compositional variations of bacteria and protists. Archaea were shaped by strong homogeneous selection, whereas fungi were more affected by dispersal limitation. Such differing ecological processes resulted in the domain-level differentiation of microbial spatial scaling. This study links ecological processes with microbial spatial scaling and provides novel mechanistic insights into the diversity patterns of microbes that belong to different trophic levels. IMPORTANCE As the most diverse and numerous life form on Earth, microorganisms play indispensable roles in natural ecological processes. Revealing their diversity patterns across space and through time is of essential importance to better understand the underlying ecological mechanisms controlling the distribution and assembly of microbial communities. However, the diversity patterns and their underlying ecological mechanisms for different microbial domains and/or trophic levels require further exploration. In this study, the spatial scaling of different microbial domains and their associated ecological processes in a mudflat intertidal zone were investigated. The results showed different spatial scaling patterns for different microbial domains. Different ecological processes underlie the domain-level differentiation of microbial spatial scaling. This study links ecological processes with microbial spatial scaling to provide novel mechanistic insights into the diversity patterns of microorganisms that belong to different trophic levels.

Beyond the limits of the unassigned protist microbiome: inferring large-scale spatio-temporal patterns of Syndiniales marine parasites

Marine protists are major components of the oceanic microbiome that remain largely unrepresented in culture collections and genomic reference databases. The exploration of this uncharted protist diversity in oceanic communities relies essentially on studying genetic markers from the environment as taxonomic barcodes. Here we report that across 6 large scale spatio-temporal planktonic surveys, half of the genetic barcodes remain taxonomically unassigned at the genus level, preventing a fine ecological understanding for numerous protist lineages. Among them, parasitic Syndiniales (Dinoflagellata) appear as the least described protist group. We have developed a computational workflow, integrating diverse 18S rDNA gene metabarcoding datasets, in order to infer large-scale ecological patterns at 100% similarity of the genetic marker, overcoming the limitation of taxonomic assignment. From a spatial perspective, we identified 2171 unassigned clusters, i.e., Syndiniales sequences with 100% similarity, exclusively shared between the Tropical/Subtropical Ocean and the Mediterranean Sea among all Syndiniales orders and 25 ubiquitous clusters shared within all the studied marine regions. From a temporal perspective, over 3 time-series, we highlighted 39 unassigned clusters that follow rhythmic patterns of recurrence and are the best indicators of parasite community's variation. These clusters withhold potential as ecosystem change indicators, mirroring their associated host community responses. Our results underline the importance of Syndiniales in structuring planktonic communities through space and time, raising questions regarding host-parasite association specificity and the trophic mode of persistent Syndiniales, while providing an innovative framework for prioritizing unassigned protist taxa for further description.

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

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

Environmental DNA sequencing reveals the regional difference in diversity and community assembly mechanisms of eukaryotic plankton in coastal waters

The diversity and community assembly mechanisms of eukaryotic plankton in coastal waters is so far not clear. In this study, we selected the coastal waters of Guangdong-Hong Kong-Macao Greater Bay Area, which is a highly developed region in China, as the research area. By use of high-throughput sequencing technologies, the diversity and community assembly mechanisms of eukaryotic marine plankton were studied in which a total of 7,295 OTUs were obtained, and 2,307 species were annotated by doing environmental DNA survey of 17 sites consist of surface and bottom layer. Ultimately, the analysis reveals that the species abundance of bottom layer is, by and large, higher than that in the surface layer. In the bottom, Arthropoda is the first largest group, accounting for more than 20% while Arthropoda and Bacillariophyta are dominant groups in surface waters accounting for more than 40%. It is significant of the variance in alpha-diversity between sampling sites, and the difference of alpha-diversity between bottom sites is greater than that of surface sites. The result suggests that the environmental factors that have significant influence on alpha-diversity are total alkalinity and offshore distance for surface sites, and water depth and turbidity for bottom sites. Likewise, the plankton communities obey the typical distance-decay pattern. Analysis about community assembly mechanisms reveals that, overall, dispersal limitation is the major pattern of community formation, which accounts for more than 83% of the community formation processes, suggesting that stochastic processes are the crucial assembly mechanism of the eukaryotic plankton community in the study area.

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

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

Microbial food webs share similar biogeographic patterns and driving mechanisms with depths in oligotrophic tropical western Pacific Ocean

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.

Geographic Scale Influences the Interactivities Between Determinism and Stochasticity in the Assembly of Sedimentary Microbial Communities on the South China Sea Shelf

Determinism and stochasticity in microbial community composition decisions have attracted wide attention. However, there is no consensus on their interrelationships and relative importance, and the mechanism controlling the interaction between the two ecological processes remains to be revealed. The interaction of the two ecological processes on the continental shelf of the South China Sea was studied by performing 16S rRNA gene amplicon sequencing on 90 sediments at multiple depths in five sites. Three nearshore sites have higher microbial diversity than those two close to the shelf margin. Different microbial composition was observed between sites and microbial composition of nearshore sites was positively correlated with total nitrogen, total sulfur, total organic carbon, and dissolved oxygen, while that of offshore was positively correlated with total carbon, salinity, and photosynthetically active radiation. The null model test showed that the community composition among layers of the same site and between nearby sites was mainly dominated by the homogeneous selection, while that between distant sites was mainly affected by dispersal limitation, which indicates that geographic scale influences the interactivities of determinism and stochasticity. Our research indicates that the balance of these two ecological processes along the geographic scale is mainly determined by the dispersal ability of microbes and environmental heterogeneity between areas. The study provides new insights into how deterministic and stochastic processes shape microbial community composition on the continental shelf.

Stochastic and Deterministic Processes Regulate Phytoplankton Assemblages in a Temperate Coastal Ecosystem

Assessing the relative contributions of the interacting deterministic and stochastic ecological processes for phytoplankton community assembly is crucial in understanding and predicting community organization and succession at different temporal and spatial scales. In this study, we hypothesized that deterministic and stochastic ecological processes regulating phytoplankton, present seasonal and repeating patterns. This hypothesis was explored during a 5-year survey (287 samples) conducted at a small spatial scale (~15km) in a temperate coastal ecosystem (eastern English Channel). Microscopy and flow cytometry quantified phytoplankton abundance and biomass, while metabarcoding data allowed an extended evaluation of diversity and the exploration of the ecological processes regulating phytoplankton using null model analysis. Alpha diversity of phytoplankton was governed by the effect of environmental conditions (environmental filtering). Temporal community turnover (beta diversity) evidenced a consistent interannual pattern that determined the phytoplankton seasonal structure. In winter and early spring (from January to March), determinism (homogeneous selection) was the major process in the phytoplankton community assembly. The overall mean in the year was 38%. Stochastic processes (ecological drift) prevailed during the rest of the year from April to December, where the overall mean for the year was 55%. The maximum values were recorded in late spring and summer, which often presented recurrent and transient monospecific phytoplankton peaks. Overall, the prevalence of stochastic processes rendered less predictable seasonal dynamics of phytoplankton communities to future environmental change. IMPORTANCE While ecological deterministic processes are conducive to modeling, stochastic ones are far less predictable. Understanding the overall assembly processes of phytoplankton is critical in tracking and predicting future changes. The novelty of this study was that it addressed a long-posed question, on a pluriannual scale. Was seasonal phytoplankton succession influenced by deterministic processes (e.g., abiotic environment) or by stochastic ones (e.g., dispersal, or ecological drift)? Our results provided strong support for a seasonal and repeating pattern with stochastic processes (drift) prevailing during most of the year and periods with monospecific phytoplankton peaks.

Selection, drift and community interactions shape microbial biogeographic patterns in the Pacific Ocean

Despite accumulating data on microbial biogeographic patterns in terrestrial and aquatic environments, we still lack a comprehensive understanding of how these patterns establish, in particular in ocean basins. Here we show the relative significance of the ecological mechanisms selection, dispersal and drift for shaping the composition of microbial communities in the Pacific Ocean over a transect of 12,400 km between subantarctic and subarctic regions. In the epipelagic, homogeneous selection contributes 50-60% and drift least to the three mechanism for the assembly of prokaryotic communities whereas in the upper mesopelagic, drift is relatively most important for the particle-associated subcommunities. Temperature is important for the relative significance of homogeneous selection and dispersal limitation for community assembly. The relative significance of both mechanisms was inverted with increasing temperature difference along the transect. For eukaryotes >8 µm, homogeneous selection is also the most important mechanisms at two epipelagic depths whereas at all other depths drift is predominant. As species interactions are essential for structuring microbial communities we further analyzed co-occurrence-based community metrics to assess biogeographic patterns over the transect. These interaction-adjusted indices explained much better variations in microbial community composition as a function of abiotic and biotic variables than compositional or phylogenetic distance measures like Bray-Curtis or UniFrac. Our analyses are important to better understand assembly processes of microbial communities in the upper layers of the largest ocean and how they adapt to effectively perform in global biogeochemical processes. Similar principles presumably act upon microbial community assembly in other ocean basins.

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.

Response of microbiota to exogenous inoculation improved the enzymatic activities of medium-temperature Daqu

To explore the potential mechanism of improving enzymatic activities in medium-temperature Daqu (MTD) by inoculation functional isolates, we inoculated a single strain of Bacillus licheniformis, and the microbiota composed of Bacillus velezensis and Bacillus subtilis in MTD to investigate the association between the response of the functional microbiota and the enzymatic activity. The results showed that the bacterial community of MTD might be more sensitive to bioturbation than the fungal community, and the indigenous microbiota responded to the single strain more than to the microbiota. Moreover, the differential microorganisms mainly included Lactobacillales, Bacillales, and Saccharomycetales between the conventional and fortified samples. Notably, the composition of functional microbiota related to liquefying activity (LA) and saccharifying activity (SA) were significantly different, changing from Lactobacillus and Rhizomucor to Bacillus, Weissella, and Hyphopichia. That might be closely related to the effect of the bioturbation on LA (31.33%) and SA (43.54%) associated microorganisms was more tellingly. Furthermore, the relative abundance changes of bioturbation-sensitive modules in the co-occurrence network might also lead to the difference in enzymatic activities. Therefore, the LA and SA of MTD were improved by bioturbation significantly. These results provide diverse insights into the exogenous functional isolates to regulate the MTD microbiota and improve enzymatic activities.

metaPR2 : A database of eukaryotic 18S rRNA metabarcodes with an emphasis on protists

In recent years, metabarcoding has become the method of choice for investigating the composition and assembly of microbial eukaryotic communities. The number of environmental data sets published has increased very rapidly. Although unprocessed sequence files are often publicly available, processed data, in particular clustered sequences, are rarely available in a usable format. Clustered sequences are reported as operational taxonomic units (OTUs) with different similarity levels or more recently as amplicon sequence variants (ASVs). This hampers comparative studies between different environments and data sets, for example examining the biogeographical patterns of specific groups/species, as well analysing the genetic microdiversity within these groups. Here, we present a newly-assembled database of processed 18S rRNA metabarcodes that are annotated with the PR² reference sequence database. This database, called metaPR² , contains 41 data sets corresponding to more than 4000 samples and 90,000 ASVs. The database, which is accessible through both a web-based interface (https://shiny.metapr2.org) and an R package, should prove very useful to all researchers working on protist diversity in a variety of systems.

Unravelling the interplay of ecological processes structuring the bacterial rare biosphere

Most ecological communities harbor many rare species (i.e., the rare biosphere), however, relatively little is known about how distinct ecological processes structure their existence. Here, we used spatiotemporal data on soil bacterial communities along a natural ecosystem gradient to model the relative influences of assembly processes structuring the rare and common biospheres. We found a greater influence of homogeneous selection (i.e., imposed by spatiotemporally constant variables) mediating the assembly of the rare biosphere, whereas the common biosphere was mostly governed by variable selection (i.e., imposed by spatial and/or temporal fluctuating variables). By partitioning the different types of rarity, we found homogeneous selection to explain the prevalence of permanently rare taxa, thus suggesting their persistence at low abundances to be restrained by physiological traits. Conversely, the dynamics of conditionally rare taxa were mostly structured by variable selection, which aligns with the ability of these taxa to switch between rarity and commonness as responses to environmental spatiotemporal variations. Taken together, our study contributes to the establishment of a link between conceptual and empirical developments in the ecology of the soil microbial rare biosphere. Besides, this study provides a framework to better understand, model, and predict the existence and dynamics of microbial rare biospheres across divergent systems and scales.

Ecological Dynamics and Co-occurrences Among Prokaryotes and Microeukaryotes in a Diatom Bloom Process in Xiangshan Bay, China

Diatom blooms can significantly affect the succession of microbial communities, yet little is known about the assembly processes and interactions of microbial communities during autumn bloom events. In this study, we investigated the ecological effects of an autumn diatom bloom on prokaryotic communities (PCCs) and microeukaryotic communities (MECs), focusing on their assembly processes and interactions. The PCCs were largely dominated by Alphaproteobacteria, Gammaproteobacteria, Cyanobacteria, and Flavobacteria, while the MECs primarily included Diatomea, Dinoflagellata, and Chlorophyta. The succession of both PCCs and MECs was mainly driven by this diatom bloom and environmental factors, such as nitrate and silicate. Null modeling revealed that homogeneous selection had a more pronounced impact on the structure of PCCs compared with that of MECs. In particular, drift and dispersal limitation cannot be neglected in the assembly processes of MECs. Co-occurrence network analyses showed that Litorimicrobium, Cercozoa, Marine Group I (MGI), Cryptomonadales, Myrionecta, and Micromonas may affect the bloom process. In summary, these results elucidated the complex, robust interactions and obviously distinct assembly mechanisms of PCCs and MECs during a diatom bloom and extend our current comprehension of the ecological mechanisms and microbial interactions involved in an autumn diatom bloom process.