Marked changes in diversity and relative activity of picoeukaryotes with depth in the world ocean

Integrating depth-dependent protist dynamics and microbial interactions in spring succession of a freshwater reservoir

BACKGROUND

Protists are essential contributors to eukaryotic diversity and exert profound influence on carbon fluxes and energy transfer in freshwaters. Despite their significance, there is a notable gap in research on protistan dynamics, particularly in the deeper strata of temperate lakes. This study aimed to address this gap by integrating protists into the well-described spring dynamics of Římov reservoir, Czech Republic. Over a 2-month period covering transition from mixing to established stratification, we collected water samples from three reservoir depths (0.5, 10 and 30 m) with a frequency of up to three times per week. Microbial eukaryotic and prokaryotic communities were analysed using SSU rRNA gene amplicon sequencing and dominant protistan groups were enumerated by Catalysed Reporter Deposition-Fluorescence in situ Hybridization (CARD-FISH). Additionally, we collected samples for water chemistry, phyto- and zooplankton composition analyses.

RESULTS

Following the rapid changes in environmental and biotic parameters during spring, protistan and bacterial communities displayed swift transitions from a homogeneous community to distinct strata-specific communities. A prevalence of auto- and mixotrophic protists dominated by cryptophytes was associated with spring algal bloom-specialized bacteria in the epilimnion. In contrast, the meta- and hypolimnion showcased a development of a protist community dominated by putative parasitic Perkinsozoa, detritus or particle-associated ciliates, cercozoans, telonemids and excavate protists (Kinetoplastida), co-occurring with bacteria associated with lake snow.

CONCLUSIONS

Our high-resolution sampling matching the typical doubling time of microbes along with the combined microscopic and molecular approach and inclusion of all main components of the microbial food web allowed us to unveil depth-specific populations' successions and interactions in a deep lentic ecosystem.

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.

Planktonic eukaryotes in the Chesapeake Bay: contrasting responses of abundant and rare taxa to estuarine gradients

Phytoplankton are important drivers of aquatic ecosystem function and environmental health. Their community compositions and distributions are directly impacted by environmental processes and human activities, including in the largest estuary in North America, the Chesapeake Bay. It is crucial to uncover how planktonic eukaryotes play fundamental roles as primary producers and trophic links and sustain estuarine ecosystems. In this study, we investigated the detailed community structure and spatiotemporal variations of planktonic eukaryotes in the Chesapeake Bay across space and time for three consecutive years. A clear seasonal and spatial shift of total, abundant, and rare planktonic eukaryotes was evident, and the pattern recurred interannually. Multiple harmful algal species have been identified in the Bay with varied distribution patterns, such as Karlodinium, Heterosigma akashiwo, Protoperidinium sp., etc. Compared to abundant taxa, rare subcommunities were more sensitive to environmental disturbance in terms of richness, diversity, and distribution. The combined effects of temporal variation (13.3%), nutrient availability (10.0%), and spatial gradients (8.8%) structured the distribution of eukaryotic microbial communities in the Bay. Similar spatiotemporal patterns between planktonic prokaryotes and eukaryotes suggest common mechanisms of adjustment, replacement, and species interaction for planktonic microbiomes under strong estuarine gradients. To our best knowledge, this work represents the first systematic study on planktonic eukaryotes in the Bay. A comprehensive view of the distribution of planktonic microbiomes and their interactions with environmental processes is critical in understanding the underlying microbial mechanisms involved in maintaining the stability, function, and environmental health of estuarine ecosystems.

IMPORTANCE

Deep sequencing analysis of planktonic eukaryotes in the Chesapeake Bay reveals high community diversity with many newly recognized phytoplankton taxa. The Chesapeake Bay planktonic eukaryotes show distinct seasonal and spatial variability, with recurring annual patterns of total, abundant, and rare groups. Rare taxa mainly contribute to eukaryotic diversity compared to abundant groups, and they are more sensitive to spatiotemporal variations and environmental filtering. Temporal variations, nutrient availability, and spatial gradients significantly affect the distribution of eukaryotic microbial communities. Similar spatiotemporal patterns in prokaryotes and eukaryotes suggest common mechanisms of adjustment, substitution, and species interactions in planktonic microbiomes under strong estuarine gradients. Interannually recurring patterns demonstrate that diverse eukaryotic taxa have well adapted to the estuarine environment with a long residence time. Further investigations of how human activities impact estuarine planktonic eukaryotes are critical in understanding their essential ecosystem roles and in maintaining environmental safety and public health.

Insight into planktonic protistan and fungal communities across the nutrient-depleted environment of the South Pacific Subtropical Gyre

Ocean microorganisms constitute ~70% of the marine biomass, contribute to ~50% of the Earth's primary production, and play a vital role in global biogeochemical cycles. The marine heterotrophic and mixotrophic protistan and fungal communities have often been overlooked mainly due to limitations in morphological species identification. Despite the accumulation of studies on biogeographic patterns observed in microbial communities, our understanding of the abundance and distribution patterns within the microbial community of the largest subtropical gyre, the South Pacific Gyre (SPG), remains incomplete. Here, we investigated the diversity and vertical composition of protistan and fungal communities in the water column of the ultra-oligotrophic SPG. Our results showed apparent differences in protistan community diversity in the photic and aphotic regions. The entire protistan community diversity was significantly affected by temperature, salinity, oxygen, and nutrient concentrations, while the parasitic community diversity was also affected by chlorophyll a concentration. The parasitic protists were assigned to the class Syndiniales accounting for over 98% of the total parasitic protists, exhibiting higher relative sequence abundance along the water depth and displaying consistent patterns among different sampling stations. In contrast to the protistan community, the fungal community along the SPG primarily clustered based on the sampling station and pelagic zones. In particular, our study reveals a significant presence of parasitic protists and functionally diverse fungi in SPG and their potential impact on carbon cycling in the gyre.IMPORTANCEOur findings carry important implications for understanding the distribution patterns of the previously unrecognized occurrence of parasitic protists and functionally diverse fungi in the nutrient-limited South Pacific Gyre. In particular, our study reveals a significant presence of parasitic Syndiniales, predominantly abundant in the upper 300 m of the aphotic zone in the gyre, and a distinct presence of fungal communities in the aphotic zone at the central part of the gyre. These findings strongly suggest that these communities play a substantial role in yet insufficiently described microbial food web. Moreover, our research enhances our understanding of their contribution to the dynamics of the food webs in oligotrophic gyres and is valuable for projecting the ecological consequences of future climate warming.

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.

Protistan community composition and metabolism in the North Pacific Subtropical Gyre: Influences of mesoscale eddies and depth

Marine protists and their metabolic activities are intricately tied to the cycling of nutrients and the flow of energy through microbial food webs. Physiochemical changes in the environment, such as those that result from mesoscale eddies, may impact protistan communities, but the effects that such changes have on protists are poorly known. A metatranscriptomic study was conducted to investigate how eddies affected protists at adjacent cyclonic and anticyclonic eddy sites in the oligotrophic ocean at four depths from 25 to 250 m. Eddy polarity impacted protists at all depths sampled, although the effects of eddy polarity were secondary to the impact of depth across the depth range. Eddy-induced vertical shifts in the water column yielded differences in the cyclonic and anticyclonic eddy protistan communities, and these differences were the most pronounced at and just below the deep chlorophyll maximum. An analysis of transcripts associated with protistan nutritional physiology at 150 m revealed that cyclonic eddies may support a more heterotrophic community, while anticyclonic eddies promote a more phototrophic community. The results of this study indicate that eddies alter the metabolism of protists particularly in the lower euphotic zone and may therefore impact carbon export from the euphotic zone.

The active free-living bathypelagic microbiome is largely dominated by rare surface taxa

A persistent microbial seed bank is postulated to sustain the marine biosphere, and recent findings show that prokaryotic taxa present in the ocean's surface dominate prokaryotic communities throughout the water column. Yet, environmental conditions exert a tight control on the activity of prokaryotes, and drastic changes in these conditions are known to occur from the surface to deep waters. The simultaneous characterization of the total (DNA) and active (i.e. with potential for protein synthesis, RNA) free-living communities in 13 stations distributed across the tropical and subtropical global ocean allowed us to assess their change in structure and diversity along the water column. We observed that active communities were surprisingly more similar along the vertical gradient than total communities. Looking at the vertical connectivity of the active vs. the total communities, we found that taxa detected in the surface sometimes accounted for more than 75% of the active microbiome of bathypelagic waters (50% on average). These active taxa were generally rare in the surface, representing a small fraction of all the surface taxa. Our findings show that the drastic vertical change in environmental conditions leads to the inactivation and disappearance of a large proportion of surface taxa, but some surface-rare taxa remain active (or with potential for protein synthesis) and dominate the bathypelagic active microbiome.

Increased contribution of parasites in microbial eukaryotic communities of different Aegean Sea coastal systems

Background-Aim

Protistan communities have a major contribution to biochemical processes and food webs in coastal ecosystems. However, related studies are scarce and usually limited in specific groups and/or sites. The present study examined the spatial structure of the entire protistan community in seven different gulfs and three different depths in a regional Mediterranean Sea, aiming to define taxa that are important for differences detected in the marine microbial network across the different gulfs studied as well as their trophic interactions.

Methods

Protistan community structure analysis was based on the diversity of the V2-V3 hypervariable region of the 18S rRNA gene. Operational taxonomic units (OTUs) were identified using a 97% sequence identity threshold and were characterized based on their taxonomy, trophic role, abundance and niche specialization level. The differentially abundant, between gulfs, OTUs were considered for all depths and interactions amongst them were calculated, with statistic and network analysis.

Results

It was shown that Dinophyceae, Bacillariophyta and Syndiniales were the most abundant groups, prevalent in all sites and depths. Gulfs separation was more striking at surface corroborating with changes in environmental factors, while it was less pronounced in higher depths. The study of differentially abundant, between gulfs, OTUs revealed that the strongest biotic interactions in all depths occurred between parasite species (mainly Syndiniales) and other trophic groups. Most of these species were generalists but not abundant highlighting the importance of rare species in protistan community assemblage.

Conclusion

Overall this study revealed the emergence of parasites as important contributors in protistan network regulation regardless of depth.

Globally-distributed microbial eukaryotes exhibit endemism at deep-sea hydrothermal vents

Single-celled microbial eukaryotes inhabit deep-sea hydrothermal vent environments and play critical ecological roles in the vent-associated microbial food web. 18S rRNA amplicon sequencing of diffuse venting fluids from four geographically- and geochemically-distinct hydrothermal vent fields was applied to investigate community diversity patterns among protistan assemblages. The four vent fields include Axial Seamount at the Juan de Fuca Ridge, Sea Cliff and Apollo at the Gorda Ridge, all in the NE Pacific Ocean, and Piccard and Von Damm at the Mid-Cayman Rise in the Caribbean Sea. We describe species diversity patterns with respect to hydrothermal vent field and sample type, identify putative vent endemic microbial eukaryotes, and test how vent fluid geochemistry may influence microbial community diversity. At a semi-global scale, microbial eukaryotic communities at deep-sea vents were composed of similar proportions of dinoflagellates, ciliates, Rhizaria, and stramenopiles. Individual vent fields supported distinct and highly diverse assemblages of protists that included potentially endemic or novel vent-associated strains. These findings represent a census of deep-sea hydrothermal vent protistan communities. Protistan diversity, which is shaped by the hydrothermal vent environment at a local scale, ultimately influences the vent-associated microbial food web and the broader deep-sea carbon cycle.

Global perspective of environmental distribution and diversity of Perkinsea (Alveolata) explored by a meta-analysis of eDNA surveys

Perkinsea constitutes a lineage within the Alveolata eukaryotic superphylum, mainly composed of parasitic organisms. Some described species represent significant ecological and economic threats due to their invasive ability and pathogenicity, which can lead to mortality events. However, the genetic diversity of these described species is just the tip of the iceberg. Environmental surveys targeting this lineage are still scarce and mainly limited to the Northern Hemisphere. Here, we aim to conduct an in depth exploration of the Perkinsea group, uncovering the diversity across a variety of environments, including those beyond freshwater and marine ecosystems. We seek to identify and describe putative novel organisms based on their genetic signatures. In this study, we conducted an extensive analysis of a metabarcoding dataset, focusing on the V4 region of the 18S rRNA gene (the EukBank dataset), to investigate the diversity, distribution and environmental preferences of the Perkinsea. Our results reveal a remarkable diversity within the Perkinsea, with 1568 Amplicon Sequence Variants (ASVs) identified across thousands of environmental samples. Surprisingly, we showed a substantial diversity of Perkinsea within soil samples (269 ASVs), challenging the previous assumption that this group is confined to marine and freshwater environments. In addition, we revealed that a notable proportion of Perkinsea ASVs (428 ASVs) could correspond to putative new organisms, encompassing the well-established taxonomic group Perkinsidae. Finally, our study shed light on previously unveiled taxonomic groups, including the Xcellidae, and revealed their environmental distribution. These findings demonstrate that Perkinsea exhibits far greater diversity than previously detected and surprisingly extends beyond marine and freshwater environments. The meta-analysis conducted in this study has unveiled the existence of previously unknown clusters within the Perkinsea lineage, solely identified based on their genetic signatures. Considering the ecological and economic importance of described Perkinsea species, these results suggest that Perkinsea may play a significant, yet previously unrecognized, role across a wide range of environments, spanning from soil environments to the abyssal zone of the open ocean with important implications for ecosystem functioning.

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.

Interpreting phylogenetic placements for taxonomic assignment of environmental DNA

Taxonomic assignment of operational taxonomic units (OTUs) is an important bioinformatics step in analyzing environmental sequencing data. Pairwise alignment and phylogenetic-placement methods represent two alternative approaches to taxonomic assignments, but their results can differ. Here we used available colpodean ciliate OTUs from forest soils to compare the taxonomic assignments of VSEARCH (which performs pairwise alignments) and EPA-ng (which performs phylogenetic placements). We showed that when there are differences in taxonomic assignments between pairwise alignments and phylogenetic placements at the subtaxon level, there is a low pairwise similarity of the OTUs to the reference database. We then showcase how the output of EPA-ng can be further evaluated using GAPPA to assess the taxonomic assignments when there exist multiple equally likely placements of an OTU, by taking into account the sum over the likelihood weights of the OTU placements within a subtaxon, and the branch distances between equally likely placement locations. We also inferred the evolutionary and ecological characteristics of the colpodean OTUs using their placements within subtaxa. This study demonstrates how to fully analyze the output of EPA-ng, by using GAPPA in conjunction with knowledge of the taxonomic diversity of the clade of interest.

Ocean-wide comparisons of mesopelagic planktonic community structures

For decades, marine plankton have been investigated for their capacity to modulate biogeochemical cycles and provide fishery resources. Between the sunlit (epipelagic) layer and the deep dark waters, lies a vast and heterogeneous part of the ocean: the mesopelagic zone. How plankton composition is shaped by environment has been well-explored in the epipelagic but much less in the mesopelagic ocean. Here, we conducted comparative analyses of trans-kingdom community assemblages thriving in the mesopelagic oxygen minimum zone (OMZ), mesopelagic oxic, and their epipelagic counterparts. We identified nine distinct types of intermediate water masses that correlate with variation in mesopelagic community composition. Furthermore, oxygen, NO3- and particle flux together appeared as the main drivers governing these communities. Novel taxonomic signatures emerged from OMZ while a global co-occurrence network analysis showed that about 70% of the abundance of mesopelagic plankton groups is organized into three community modules. One module gathers prokaryotes, pico-eukaryotes and Nucleo-Cytoplasmic Large DNA Viruses (NCLDV) from oxic regions, and the two other modules are enriched in OMZ prokaryotes and OMZ pico-eukaryotes, respectively. We hypothesize that OMZ conditions led to a diversification of ecological niches, and thus communities, due to selective pressure from limited resources. Our study further clarifies the interplay between environmental factors in the mesopelagic oxic and OMZ, and the compositional features of communities.

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.

Mechanistic insights into dissolved organic matter-driven protistan and bacterial community dynamics influenced by vegetation restoration

Vegetation restoration projects can not only improve water quality by absorbing and transferring pollutants and nutrients from non-vegetation sources, but also protect biodiversity by providing habitat for biological growth. However, the mechanism of the protistan and bacterial assembly processes in the vegetation restoration project were rarely explored. To address this, based on 18 S rRNA and 16 S rRNA high-throughput sequencing, we investigated the mechanism of protistan and bacterial community assembly processes, environmental conditions, and microbial interactions in the rivers with (out) vegetation restoration. The results indicated that the deterministic process dominated the protistan and bacterial community assembly (94.29% and 92.38%), influenced by biotic and abiotic factors. For biotic factors, microbial network connectivity was higher in the vegetation zone (average degree = 20.34) than in the bare zone (average degree = 11.00). For abiotic factors, the concentration of dissolved organic carbon ([DOC]) was the most important environmental factor affecting the microbial community composition. [DOC] was lower significantly in vegetation zone (18.65 ± 6.34 mg/L) than in the bare zone (28.22 ± 4.82 mg/L). In overlying water, vegetation restoration upregulated the protein-like fluorescence components (C1 and C2) by 1.26 and 1.01-folds and downregulated the terrestrial humic-like fluorescence components (C3 and C4) by 0.54 and 0.55-folds, respectively. The different DOM components guided bacteria and protists to select different interactive relationships. The protein-like DOM components led to bacterial competition, whereas the humus-like DOM components resulted in protistan competition. Finally, the structural equation model was established to explain that DOM components can affect protistan and bacterial diversity by providing substrates, facilitating microbial interactions, and promoting nutrient input. In general, our study provides insights into the responses of vegetation restored ecosystems to the dynamics and interactives in the anthropogenically influenced river and evaluates the ecological restoration performance of vegetation restoration from a molecular biology perspective.

RNA outperforms DNA-based metabarcoding in assessing the diversity and response of microeukaryotes to environmental variables in the Arctic Ocean

The Arctic Ocean (AO) has a harsh environment characterized by low temperatures, extensive ice coverage, and periodic freezing and melting of sea ice, which has provided diverse habitats for microorganisms. Prior studies primarily focused on microeukaryote communities in the upper water or sea ice based on environmental DNA, leaving the composition of active microeukaryotes in the diverse AO environments largely unknown. This study provided a vertical assessment of microeukaryote communities in the AO from snow and ice to sea water at a depth of 1670 m using high-throughput sequencing of co-extracted DNA and RNA. RNA extracts depicted microeukaryote community structure and intergroup correlations more accurately and responded more sensitively to environmental conditions than those derived from DNA. Using RNA:DNA ratios as a proxy for relative activity of major taxonomic groups, the metabolic activities of major microeukaryote groups were determined along depth. Analysis of co-occurrence networks showed that parasitism between Syndiniales and dinoflagellates/ciliates in the deep ocean may be significant. This study increased our knowledge of the diversity of active microeukaryote communities and highlighted the importance of using RNA-based sequencing over DNA-based sequencing to examine the relationship between microeukaryote assemblages and the responses of microeukaryotes to environmental variables in the AO.

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.

Water masses shape pico-nano eukaryotic communities of the Weddell Sea

Polar oceans belong to the most productive and rapidly changing environments, yet our understanding of this fragile ecosystem remains limited. Here we present an analysis of a unique set of DNA metabarcoding samples from the western Weddell Sea sampled throughout the whole water column and across five water masses with different characteristics and different origin. We focus on factors affecting the distribution of planktonic pico-nano eukaryotes and observe an ecological succession of eukaryotic communities as the water masses move away from the surface and as oxygen becomes depleted with time. At the beginning of this succession, in the photic zone, algae, bacteriovores, and predators of small eukaryotes dominate the community, while another community develops as the water sinks deeper, mostly composed of parasitoids (syndinians), mesoplankton predators (radiolarians), and diplonemids. The strongly correlated distribution of syndinians and diplonemids along the depth and oxygen gradients suggests their close ecological link and moves us closer to understanding the biological role of the latter group in the ocean ecosystem.

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.

Global patterns and rates of habitat transitions across the eukaryotic tree of life

The successful colonization of new habitats has played a fundamental role during the evolution of life. Salinity is one of the strongest barriers for organisms to cross, which has resulted in the evolution of distinct marine and non-marine (including both freshwater and soil) communities. Although microbes represent by far the vast majority of eukaryote diversity, the role of the salt barrier in shaping the diversity across the eukaryotic tree is poorly known. Traditional views suggest rare and ancient marine/non-marine transitions but this view is being challenged by the discovery of several recently transitioned lineages. Here, we investigate habitat evolution across the tree of eukaryotes using a unique set of taxon-rich phylogenies inferred from a combination of long-read and short-read environmental metabarcoding data spanning the ribosomal DNA operon. Our results show that, overall, marine and non-marine microbial communities are phylogenetically distinct but transitions have occurred in both directions in almost all major eukaryotic lineages, with hundreds of transition events detected. Some groups have experienced relatively high rates of transitions, most notably fungi for which crossing the salt barrier has probably been an important aspect of their successful diversification. At the deepest phylogenetic levels, ancestral habitat reconstruction analyses suggest that eukaryotes may have first evolved in non-marine habitats and that the two largest known eukaryotic assemblages (TSAR and Amorphea) arose in different habitats. Overall, our findings indicate that the salt barrier has played an important role during eukaryote evolution and provide a global perspective on habitat transitions in this domain of life.

Spatial dynamics of active microeukaryotes along a latitudinal gradient: Diversity, assembly process, and co-occurrence relationships

Recent global warming is profoundly and increasingly influencing the Arctic ecosystem. Understanding how microeukaryote communities respond to changes in the Arctic Ocean is crucial for understanding their roles in the biogeochemical cycles of nutrients and elements. Between July 22 and August 19, 2016, during cruise ARA07, seawater samples were collected along a latitudinal transect extending from the East Sea of Korea to the central Arctic Ocean. Environmental RNA was extracted and the V4 hypervariable regions of the reverse transcribed SSU rRNA were amplified. The sequences generated by high throughput sequencing were clustered into zero-radius OTUs (ZOTUs), and the taxonomic identities of each ZOTU were assigned using SINTAX against the PR2 database. Thus, the diversity, community composition, and co-occurrence networks of size fractionated microeukaryotes were revealed. The present study found: 1) the alpha diversity of pico- and nano-sized microeukaryotes showed a latitudinal diversity gradient; 2) three distinct communities were identified, i.e., the Leg-A, Leg-B surface, and Leg-B subsurface chlorophyll a maximum (SCM) groups; 3) distinct network structure and composition were found in the three groups; and 4) water temperature was identified as the primary factor driving both the alpha and beta diversities of microeukaryotes. This study conducted a comprehensive and systematic survey of active microeukaryotes along a latitudinal gradient, elucidated the diversity, community composition, co-occurrence relationships, and community assembly processes among major microeukaryote assemblages, and will help shed more light on our understanding of the responses of microeukaryote communities to the changing Arctic Ocean.

Protist Diversity and Metabolic Strategy in Freshwater Lakes Are Shaped by Trophic State and Watershed Land Use on a Continental Scale

Protists play key roles in aquatic food webs as primary producers, predators, nutrient recyclers, and symbionts. However, a comprehensive view of protist diversity in freshwaters has been challenged by the immense environmental heterogeneity among lakes worldwide. We assessed protist diversity in the surface waters of 366 freshwater lakes across a north temperate to subarctic range covering nearly 8.4 million km² of Canada. Sampled lakes represented broad gradients in size, trophic state, and watershed land use. Hypereutrophic lakes contained the least diverse and most distinct protist communities relative to nutrient-poor lakes. Greater taxonomic variation among eutrophic lakes was mainly a product of heterotroph and mixotroph diversity, whereas phototroph assemblages were more similar under high-nutrient conditions. Overall, local physicochemical factors, particularly ion and nutrient concentrations, elicited the strongest responses in community structure, far outweighing the effects of geographic gradients. Despite their contrasting distribution patterns, obligate phototroph and heterotroph turnover was predicted by an overlapping set of environmental factors, while the metabolic plasticity of mixotrophs may have made them less predictable. Notably, protist diversity was associated with variation in watershed soil pH and agricultural crop coverage, pointing to human impact on the land-water interface that has not been previously identified in studies on smaller scales. Our study exposes the importance of both within-lake and external watershed characteristics in explaining protist diversity and biogeography, critical information for further developing an understanding of how freshwater lakes and their watersheds are impacted by anthropogenic stressors.

IMPORTANCE Freshwater lakes are experiencing rapid changes under accelerated anthropogenic stress and a warming climate. Microorganisms underpin aquatic food webs, yet little is known about how freshwater microbial communities are responding to human impact. Here, we assessed the diversity of protists and their myriad ecological roles in lakes varying in size across watersheds experiencing a range of land use pressures by leveraging data from a continental-scale survey of Canadian lakes. We found evidence of human impact on protist assemblages through an association with lake trophic state and extending to agricultural activity and soil characteristics in the surrounding watershed. Furthermore, trophic state appeared to explain the distributions of phototrophic and heterotrophic protists in contrasting ways. Our findings highlight the vulnerability of lake ecosystems to increased land use and the importance of assessing terrestrial interfaces to elucidate freshwater ecosystem dynamics.

Exploring the Taxonomic and Functional Diversity of Marine Benthic Micro-Eukaryotes Along the Red Sea Coast of Jeddah City

Backgrounds

Diverse marine habitats along Jeddah's Red Sea coast support rich biodiversity. Few studies have been done on its diverse communities, especially its microbial counterparts. Metagenomic analysis of marine benthic micro-eukaryotic communities was performed for the first time on the Red Sea coast of Jeddah. This research looks into their community structure and metabolic potential.

Methods

Next-generation sequencing was used to examine the micro-eukaryotic communities of seven sedimentary soil samples from four Jeddah coast locations. After isolating DNA from seven benthic sedimentary soil samples, the 18S rDNA V4 regions were amplified and sequenced on the Illumina MiSeq. It was also verified using an Agilent Technologies 2100 Bioanalyzer with a DNA 1000 chip (Agilent Technologies, Fisher Scientific). A standard curve of fluorescence readings generated by qPCR quantification using the Illumina library was achieved using the GS FLX library. Metagenomic data analysis was used to evaluate the microbial communities' biochemical and enzymatic allocations in studied samples.

Results

Blast analysis showed that the top ten phyla were Annelida, Eukaryota, Diatomea, Porifera, Phragmoplastophyta, Arthropoda, Dinoflagellata, Xenacoelomorpha Nematoda, and uncultured. Annelida was also found in the highest percentage (93%), in the sample M followed by Porifera (64%), the most abundant in the control sample then Eukaryotes (61%), Phragmatoplastophyta (55%), Arthropoda, and Diatomea (the least common) (32%). community diversity analysis: using Shannon and inverse Simpson indices showed sediment composition to be effective. Also, PICRUST2 indicated that the most abundant pathways were pyruvate fermentation to isobutanol, pyrimidine deoxyribonucleotide phosphorylation, adenosine ribonucleotide de novo biosynthesis, guanosine ribonucleotide de novo biosynthesis, NAD salvage pathway I, the super pathway of glyoxylate bypass and aerobic respiration I (cytochrome c).

Conclusion

Results showed that high throughput metagenomics could reveal species diversity and estimate gene profiles. Environmental factors appear to be more important than geographic variation in determining the structure of these microbial communities. This study provides the first report of marine benthic micro-eukaryotic communities found on the Red Sea coast of Jeddah and will serve as a good platform for future research.

Diversity of Microbial Eukaryotes Along the West Antarctic Peninsula in Austral Spring

During a cruise from October to November 2019, along the West Antarctic Peninsula, between 64.32 and 68.37°S, we assessed the diversity and composition of the active microbial eukaryotic community within three size fractions: micro- (> 20 μm), nano- (20-5 μm), and pico-size fractions (5-0.2 μm). The communities and the environmental parameters displayed latitudinal gradients, and we observed a strong similarity in the microbial eukaryotic communities as well as the environmental parameters between the sub-surface and the deep chlorophyll maximum (DCM) depths. Chlorophyll concentrations were low, and the mixed layer was shallow for most of the 17 stations sampled. The richness of the microplankton was higher in Marguerite Bay (our southernmost stations), compared to more northern stations, while the diversity for the nano- and pico-plankton was relatively stable across latitude. The microplankton communities were dominated by autotrophs, mostly diatoms, while mixotrophs (phototrophs-consuming bacteria and kleptoplastidic ciliates, mostly alveolates, and cryptophytes) were the most abundant and active members of the nano- and picoplankton communities. While phototrophy was the dominant trophic mode, heterotrophy (mixotrophy, phagotrophy, and parasitism) tended to increase southward. The samples from Marguerite Bay showed a distinct community with a high diversity of nanoplankton predators, including spirotrich ciliates, and dinoflagellates, while cryptophytes were observed elsewhere. Some lineages were significantly related-either positively or negatively-to ice coverage (e.g., positive for Pelagophyceae, negative for Spirotrichea) and temperature (e.g., positive for Cryptophyceae, negative for Spirotrichea). This suggests that climate changes will have a strong impact on the microbial eukaryotic community.

Contrasting diversity patterns of prokaryotes and protists over time and depth at the San-Pedro Ocean Time series

Community dynamics are central in microbial ecology, yet we lack studies comparing diversity patterns among marine protists and prokaryotes over depth and multiple years. Here, we characterized microbes at the San-Pedro Ocean Time series (2005-2018), using SSU rRNA gene sequencing from two size fractions (0.2-1 and 1-80 μm), with a universal primer set that amplifies from both prokaryotes and eukaryotes, allowing direct comparisons of diversity patterns in a single set of analyses. The 16S + 18S rRNA gene composition in the small size fraction was mostly prokaryotic (>92%) as expected, but the large size fraction unexpectedly contained 46-93% prokaryotic 16S rRNA genes. Prokaryotes and protists showed opposite vertical diversity patterns; prokaryotic diversity peaked at mid-depth, protistan diversity at the surface. Temporal beta-diversity patterns indicated prokaryote communities were much more stable than protists. Although the prokaryotic communities changed monthly, the average community stayed remarkably steady over 14 years, showing high resilience. Additionally, particle-associated prokaryotes were more diverse than smaller free-living ones, especially at deeper depths, contributed unexpectedly by abundant and diverse SAR11 clade II. Eukaryotic diversity was strongly correlated with the diversity of particle-associated prokaryotes but not free-living ones, reflecting that physical associations result in the strongest interactions, including symbioses, parasitism, and decomposer relationships.

Eukaryotic Parasites Are Integral to a Productive Microbial Food Web in Oxygen-Depleted Waters

Oxygen-depleted water columns (ODWCs) host a diverse community of eukaryotic protists that change dramatically in composition over the oxic-anoxic gradient. In the permanently anoxic Cariaco Basin, peaks in eukaryotic diversity occurred in layers where dark microbial activity (chemoautotrophy and heterotrophy) were highest, suggesting a link between prokaryotic activity and trophic associations with protists. Using 18S rRNA gene sequencing, parasites and especially the obligate parasitic clade, Syndiniales, appear to be particularly abundant, suggesting parasitism is an important, but overlooked interaction in ODWC food webs. Syndiniales were also associated with certain prokaryotic groups that are often found in ODWCs, including Marinimicrobia and Marine Group II archaea, evocative of feedbacks between parasitic infection events, release of organic matter, and prokaryotic assimilative activity. In a network analysis that included all three domains of life, bacterial and archaeal taxa were putative bottleneck and hub species, while a large proportion of edges were connected to eukaryotic nodes. Inclusion of parasites resulted in a more complex network with longer path lengths between members. Together, these results suggest that protists, and especially protistan parasites, play an important role in maintaining microbial food web complexity, particularly in ODWCs, where protist diversity and microbial productivity are high, but energy resources are limited relative to euphotic waters.

How habitat heterogeneity shapes bacterial and protistan communities in temperate coastal areas near estuaries

In this study, we investigated microbial communities (bacteria and protist) in two coastal areas near the estuaries of the Liaohe (LH) River and Yalujiang (YLJ) River in the Northwestern Pacific Ocean. Due to the existence of Liaodong Peninsula and different levels of urbanization, geographical segregation and significant environmental heterogeneity were observed between these two areas. There were significantly different regional species pools and biogeographic patterns for both bacterial and protistan communities between LH and YLJ coastal areas. Species turnover was the main mechanism driving β-diversity patterns of both bacterial and protistan communities in each area. In addition, the contributed ratio of nestedness to the β-diversity patterns was significantly higher for protists compared to bacteria. Variation in regional species pools was found to be the dominant driver of differences of bacterial and protistan communities between the LH and YLJ coastal areas. For a single-studied area, local community assembly mechanisms, including heterogeneous selection and dispersal limitation, were found to shape the bacterial and protistan communities through calculation of the β-deviation index. Among them, the relative importance of heterogeneous selection and dispersal limitation on the community assembly varied according to microorganism type and habitat.

Comparison of Deep-Sea Picoeukaryotic Composition Estimated from the V4 and V9 Regions of 18S rRNA Gene with a Focus on the Hadal Zone of the Mariana Trench

Diversity of microbial eukaryotes is estimated largely based on sequencing analysis of the hypervariable regions of 18S rRNA genes. But the use of different regions of 18S rRNA genes as molecular markers may generate bias in diversity estimation. Here, we compared the differences between the two most widely used markers, V4 and V9 regions of the 18S rRNA gene, in describing the diversity of epipelagic, bathypelagic, and hadal picoeukaryotes in the Challenger Deep of the Mariana Trench, which is a unique and little explored environment. Generally, the V9 region identified more OTUs in deeper waters than V4, while the V4 region provided greater Shannon diversity than V9. In the epipelagic zone, where Alveolata was the dominant group, picoeukaryotic community compositions identified by V4 and V9 markers are similar at different taxonomic levels. However, in the deep waters, the results of the two datasets show clear differences. These differences were mainly contributed by Retaria, Fungi, and Bicosoecida. The primer targeting the V9 region has an advantage in amplifying Bicosoecids in the bathypelagic and hadal zone of the Mariana Trench, and its high abundance in V9 dataset pointed out the possibility of Bicosoecids as a dominant group in this environment. Chrysophyceae, Fungi, MALV-I, and Retaria were identified as the dominant picoeukaryotes in the bathypelagic and hadal zone and potentially play important roles in deep-sea microbial food webs and biogeochemical cycling by their phagotrophic, saprotrophic, and parasitic life styles. Overall, the use of different markers of 18S rRNA gene allows a better assessment and understanding of the picoeukaryotic diversity in deep-sea environments.

The polar night shift: seasonal dynamics and drivers of Arctic Ocean microbiomes revealed by autonomous sampling

The Arctic Ocean features extreme seasonal differences in daylight, temperature, ice cover, and mixed layer depth. However, the diversity and ecology of microbes across these contrasting environmental conditions remain enigmatic. Here, using autonomous samplers and sensors deployed at two mooring sites, we portray an annual cycle of microbial diversity, nutrient concentrations and physical oceanography in the major hydrographic regimes of the Fram Strait. The ice-free West Spitsbergen Current displayed a marked separation into a productive summer (dominated by diatoms and carbohydrate-degrading bacteria) and regenerative winter state (dominated by heterotrophic Syndiniales, radiolarians, chemoautotrophic bacteria, and archaea). The autumn post-bloom with maximal nutrient depletion featured Coscinodiscophyceae, Rhodobacteraceae (e.g. Amylibacter) and the SAR116 clade. Winter replenishment of nitrate, silicate and phosphate, linked to vertical mixing and a unique microbiome that included Magnetospiraceae and Dadabacteriales, fueled the following phytoplankton bloom. The spring-summer succession of Phaeocystis, Grammonema and Thalassiosira coincided with ephemeral peaks of Aurantivirga, Formosa, Polaribacter and NS lineages, indicating metabolic relationships. In the East Greenland Current, deeper sampling depth, ice cover and polar water masses concurred with weaker seasonality and a stronger heterotrophic signature. The ice-related winter microbiome comprised Bacillaria, Naviculales, Polarella, Chrysophyceae and Flavobacterium ASVs. Low ice cover and advection of Atlantic Water coincided with diminished abundances of chemoautotrophic bacteria while others such as Phaeocystis increased, suggesting that Atlantification alters microbiome structure and eventually the biological carbon pump. These insights promote the understanding of microbial seasonality and polar night ecology in the Arctic Ocean, a region severely affected by climate change.

Diversity and ecology of protists revealed by metabarcoding

Protists are the dominant eukaryotes in the biosphere where they play key functional roles. While protists have been studied for over a century, it is the high-throughput sequencing of molecular markers from environmental samples - the approach of metabarcoding - that has revealed just how diverse, and abundant, these small organisms are. Metabarcoding is now routine to survey environmental diversity, so data have rapidly accumulated from a multitude of environments and at different sampling scales. This mass of data has provided unprecedented opportunities to study the taxonomic and functional diversity of protists, and how this diversity is organised in space and time. Here, we use metabarcoding as a common thread to discuss the state of knowledge in protist diversity research, from technical considerations of the approach to important insights gained on diversity patterns and the processes that might have structured this diversity. In addition to these insights, we conclude that metabarcoding is on the verge of an exciting added dimension thanks to the maturation of high-throughput long-read sequencing, so that a robust eco-evolutionary framework of protist diversity is within reach.

Single-Cell Transcriptomics of Abedinium Reveals a New Early-Branching Dinoflagellate Lineage

Dinoflagellates possess many cellular characteristics with unresolved evolutionary histories. These include nuclei with greatly expanded genomes and chromatin packaged using histone-like proteins and dinoflagellate-viral nucleoproteins instead of histones, highly reduced mitochondrial genomes with extensive RNA editing, a mix of photosynthetic and cryptic secondary plastids, and tertiary plastids. Resolving the evolutionary origin of these traits requires understanding their ancestral states and early intermediates. Several early-branching dinoflagellate lineages are good candidates for such reconstruction, however these cells tend to be delicate and environmentally sparse, complicating such analyses. Here, we employ transcriptome sequencing from manually isolated and microscopically documented cells to resolve the placement of two cells of one such genus, Abedinium, collected by remotely operated vehicle in deep waters off the coast of Monterey Bay, CA. One cell corresponds to the only described species, Abedinium dasypus, whereas the second cell is distinct and formally described as Abedinium folium, sp. nov. Abedinium has classically been assigned to the early-branching dinoflagellate subgroup Noctilucales, which is weakly supported by phylogenetic analyses of small subunit ribosomal RNA, the single characterized gene from any member of the order. However, an analysis based on 221 proteins from the transcriptome places Abedinium as a distinct lineage, separate from and basal to Noctilucales and the rest of the core dinoflagellates. The transcriptome also contains evidence of a cryptic plastid functioning in the biosynthesis of isoprenoids, iron-sulfur clusters, and heme, a mitochondrial genome with all three expected protein-coding genes (cob, cox1, and cox3), and the presence of some but not all dinoflagellate-specific chromatin packaging proteins.

Protistan grazing impacts microbial communities and carbon cycling at deep-sea hydrothermal vents

Microbial eukaryotes (or protists) in marine ecosystems are a link between primary producers and all higher trophic levels, and the rate at which heterotrophic protistan grazers consume microbial prey is a key mechanism for carbon transport and recycling in microbial food webs. At deep-sea hydrothermal vents, chemosynthetic bacteria and archaea form the base of a food web that functions in the absence of sunlight, but the role of protistan grazers in these highly productive ecosystems is largely unexplored. Here, we pair grazing experiments with a molecular survey to quantify protistan grazing and to characterize the composition of vent-associated protists in low-temperature diffuse venting fluids from Gorda Ridge in the northeast Pacific Ocean. Results reveal protists exert higher predation pressure at vents compared to the surrounding deep seawater environment and may account for consuming 28 to 62% of the daily stock of prokaryotic biomass within discharging hydrothermal vent fluids. The vent-associated protistan community was more species rich relative to the background deep sea, and patterns in the distribution and co-occurrence of vent microbes provide additional insights into potential predator-prey interactions. Ciliates, followed by dinoflagellates, Syndiniales, rhizaria, and stramenopiles, dominated the vent protistan community and included bacterivorous species, species known to host symbionts, and parasites. Our findings provide an estimate of protistan grazing pressure within hydrothermal vent food webs, highlighting the important role that diverse protistan communities play in deep-sea carbon cycling.

Changes in protist communities in drainages across the Pearl River Delta under anthropogenic influence

Drainages in the Pearl River Delta urban agglomeration (PRDUA) host vital aquatic ecosystems and face enormous pressures from human activities in one of the largest urban agglomerations in the world. Despite being crucial components of aquatic ecosystems, the interactions and assembly processes of the protistan community are rarely explored in areas with serious anthropogenic disturbance. To elucidate the mechanisms of these processes, we used environmental DNA sequencing of 18S rDNA to investigate the influence of environmental factors and species interactions on the protistan community and its assembly in drainages of the PRDUA during summer. The protistan community showed a high level of diversity and a marked spatial pattern in this region. Community assembly was driven primarily by stochastic processes based on the Sloan neutral community model, explaining 74.28%, 75.82%, 73.67%, 74.40% and 51.24% of community variations in the BJ (Beijiang), XJ (Xijiang), PRD (Pearl River Delta), PRE (Pearl River Estuary) areas and in total, respectively. Meanwhile, environmental variables including temperature, pH, dissolved oxygen, transparency, nutrients and land use were strongly correlated with the composition and assembly of the protistan community, explaining 40.40% of variation in the protistan community. Furthermore, the bacterial community was simultaneously analysed by the 16S rDNA sequencing. Co-occurrence network analysis revealed that species interactions within bacteria (81.41% positive) or protists (82.80% positive), and those between bacteria and protists (50% positive and 50% negative) impacted the protistan community assembly. In summary, stochastic processes dominated, whereas species interactions and environmental factors also played important roles in shaping the protistan communities in drainages across the PRDUA. This study provides insights into the ecological patterns, assembly processes and species interactions underlying protistan dynamics in urban aquatic ecosystems experiencing serious anthropogenic disturbance.

Swarm v3: towards tera-scale amplicon clustering

MOTIVATION

Previously we presented swarm, an open-source amplicon clustering programme that produces fine-scale molecular operational taxonomic units (OTUs) that are free of arbitrary global clustering thresholds. Here, we present swarm v3 to address issues of contemporary datasets that are growing towards tera-byte sizes.

RESULTS

When compared with previous swarm versions, swarm v3 has modernized C++ source code, reduced memory footprint by up to 50%, optimized CPU-usage and multithreading (more than 7 times faster with default parameters), and it has been extensively tested for its robustness and logic.

AVAILABILITY AND IMPLEMENTATION

Source code and binaries are available at https://github.com/torognes/swarm.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A Comparison of DNA Metabarcoding and Microscopy Methodologies for the Study of Aquatic Microbial Eukaryotes

The procedures and methodologies employed to study microbial eukaryotic plankton have been thoroughly discussed. Two main schools exist—one insisting on classic microscopy methodologies and the other supporting modern high-throughput sequencing (DNA metabarcoding). However, few studies have attempted to combine both these approaches; most studies implement one method while ignoring the other. This work aims to contribute to this discussion and examine the advantages and disadvantages of each methodology by comparing marine plankton community results from microscopy and DNA metabarcoding. The results obtained by the two methodologies do not vary significantly for Bacillariophyta, although they do for Dinoflagellata and Ciliophora. The lower the taxonomic level, the higher the inconsistency between the two methodologies for all the studied groups. Considering the different characteristics of microscopy-based identification and DNA metabarcoding, this work underlines that each method should be chosen depending on the aims of the study. DNA metabarcoding provides a better estimate of the taxonomic richness of an ecosystem while microscopy provides more accurate quantitative results regarding abundance and biomass. In any case, the combined use of the two methods, if properly standardized, can provide much more reliable and accurate results for the study of marine microbial eukaryotes.

Optimal proteome allocation and the temperature dependence of microbial growth laws

Although the effect of temperature on microbial growth has been widely studied, the role of proteome allocation in bringing about temperature-induced changes remains elusive. To tackle this problem, we propose a coarse-grained model of microbial growth, including the processes of temperature-sensitive protein unfolding and chaperone-assisted (re)folding. We determine the proteome sector allocation that maximizes balanced growth rate as a function of nutrient limitation and temperature. Calibrated with quantitative proteomic data for Escherichia coli, the model allows us to clarify general principles of temperature-dependent proteome allocation and formulate generalized growth laws. The same activation energy for metabolic enzymes and ribosomes leads to an Arrhenius increase in growth rate at constant proteome composition over a large range of temperatures, whereas at extreme temperatures resources are diverted away from growth to chaperone-mediated stress responses. Our approach points at risks and possible remedies for the use of ribosome content to characterize complex ecosystems with temperature variation.

Bathymetric gradient shapes the community composition rather than the species richness of deep-sea benthic ciliates

The bathymetric gradient is one of the most important factors that regulate the distribution of life. However, community variations of benthic ciliates along bathymetric gradients in the deep sea remain rather unexplored. In this study, we hypothesize that in the deep sea, the bathymetric gradient shapes the benthic ciliate community composition rather than the species richness. Here, we evaluated the distribution patterns and drivers of benthic ciliate communities of an abyssal plain, a seamount, and a trench with water depths ranging from 800 m down to 6600 m by high throughput eDNA sequencing and statistical analyses. We observed no significant correlation between ciliate operated taxonomic unit (OTU) richness and water depth. A meta-analysis, which combined our previously published data from the neritic habitats, supports the notion that water depth exceeding 800 m has little effect on the richness of benthic ciliate species. In contrast, the composition of deep-sea ciliate communities was significantly distinct in different habitats along the bathymetric gradients. A SourceTracker analysis revealed extremely low connectivity among ciliate communities along the bathymetric gradients. More than 95% of the community dissimilarity in the deep-sea floor was attributed to species replacement, which might be caused by environmental sorting or historical constraints. Furthermore, the observed community variations could be ascribed more to water depth than to geographic distance. The findings imply that the strong force of environmental sorting along the bathymetric gradients and the low connectivity among the ciliate communities might lead to an isolated evolution. This could shape the community composition rather than the species richness, which is mainly determined by the limited nutrient availability and the extreme environmental conditions in the deep sea.

Shift and Metabolic Potentials of Microbial Eukaryotic Communities Across the Full Depths of the Mariana Trench

Microbial eukaryotes are widespread and play important roles in marine ecosystems. However, their ecological characteristics in the deep sea (>1,000 m), especially hadal trenches, were largely unknown. Here, we investigated the diversity and metabolic potentials of microbial eukaryotes along the whole water column of the Mariana Trench by metagenomics. Our results showed clear depth-related distribution of microbial eukaryotic community and associated metabolic potentials. Surface seawater was dominated by phototrophic/mixotrophic groups (e.g., Dinoflagellata) and genes involved in biosynthesis (photosynthesis and fatty acid biosynthesis), while deep (bathypelagic and/or hadal) seawaters were enriched with heterotrophic groups (e.g., Bicoecea) and genes related to digestion (lysosomal enzymes and V-type ATPase) and carbohydrate metabolism. Co-occurrence analysis revealed high intra-domain connectivity, indicating that microbial eukaryotic composition was more influenced by microbial eukaryotes themselves than bacteria. Increased abundance of genes associated with unsaturated fatty acid biosynthesis likely plays a role in resisting high hydrostatic pressure. Top1 and hupB genes, responsible for the formation and stabilization of DNA structure, were unique and abundant in the hadal zone and thus may be helpful to stabilize DNA structure in the deep sea. Overall, our results provide insights into the distribution and potential adaptability of microbial eukaryotes in the hadal zone.

Colpodean ciliate phylogeny and reference alignments for phylogenetic placements

The Colpodea form a major clade of ciliates that are often found in environmental DNA sequencing studies. They are united by similar somatic ciliature, but differentiated by complex oral structures. Although there are four well supported colpodean subclades, there is disagreement in molecular phylogenetic inferences about their branching order. Using available nuclear SSU-rRNA sequences, we evaluated if the bursariomorphids or the platyophryids are sister to the remaining colpodeans. We inferred the "platyophryids-early" topologies using different alignment and masking methods, but constrained analyses could not reject the "bursariomorphids-early" topology. Both bursariomorphids and platyophryids clades have a similar number of nucleotide positions shared with the outgroup, and both are interconnected with the outgroup in phylogenetic networks. Based on these discordant results, it is hard to determine which clade branched off first, although the "platyophryids-early topology" is also supported by mitochondrial SSU-rRNA data. We also offer different reference alignments that can be used to phylogenetically place short- and long-read data from environmental DNA sequencing studies, and we propose some tentative evolutionary and ecological interpretations of those placements.

The Windblown: Possible Explanations for Dinophyte DNA in Forest Soils

Dinophytes are widely distributed in marine- and fresh-waters, but have yet to be conclusively documented in terrestrial environments. Here, we evaluated the presence of these protists from an environmental DNA metabarcoding dataset of Neotropical rainforest soils. Using a phylogenetic placement approach with a reference alignment and tree, we showed that the numerous sequencing reads that were phylogenetically placed as dinophytes did not correlate with taxonomic assignment, environmental preference, nutritional mode, or dormancy. All the dinophytes in the soils are rather windblown dispersal units of aquatic species and are not biologically active residents of terrestrial environments.

Influence of Regional Oceanography and Hydrothermal Activity on Protist Diversity and Community Structure in the Okinawa Trough

Microbial eukaryotes (protists) contribute substantially to ecological functioning in marine ecosystems, but the relative importance of factors shaping protist diversity, such as environmental selection and dispersal, remains difficult to parse. Water masses of a back-arc basin with hydrothermal activity provide a unique opportunity for studying the effects of dispersal and environmental selection on protist communities. In this study, we used metabarcoding to characterize protist communities in the Okinawa Trough, a back-arc spreading basin containing at least twenty-five active hydrothermal vent fields. Water was sampled from four depths at fourteen stations spanning the length of the Okinawa Trough, including three sites influenced by nearby hydrothermal vent sites. While significant differences in community structure reflecting water depth were present, protist communities were mostly homogeneous horizontally. Protist communities in the bottom waters affected by hydrothermal activity were significantly different from communities in other bottom waters, suggesting that environmental factors can be especially important in shaping community composition under specific conditions. Amplicon sequence variants that were enriched in hydrothermally influenced bottom waters largely derived from cosmopolitan protists that were present, but rare, in other near-bottom samples, thus highlighting the importance of the rare biosphere.

Gene expression during bacterivorous growth of a widespread marine heterotrophic flagellate

Phagocytosis is a fundamental process in marine ecosystems by which prey organisms are consumed and their biomass incorporated in food webs or remineralized. However, studies searching for the genes underlying this key ecological process in free-living phagocytizing protists are still scarce, in part due to the lack of appropriate ecological models. Our reanalysis of recent molecular datasets revealed that the cultured heterotrophic flagellate Cafeteria burkhardae is widespread in the global oceans, which prompted us to design a transcriptomics study with this species, grown with the cultured flavobacterium Dokdonia sp. We compared the gene expression between exponential and stationary phases, which were complemented with three starvation by dilution phases that appeared as intermediate states. We found distinct expression profiles in each condition and identified 2056 differentially expressed genes between exponential and stationary samples. Upregulated genes at the exponential phase were related to DNA duplication, transcription and translational machinery, protein remodeling, respiration and phagocytosis, whereas upregulated genes in the stationary phase were involved in signal transduction, cell adhesion, and lipid metabolism. We identified a few highly expressed phagocytosis genes, like peptidases and proton pumps, which could be used to target this ecologically relevant process in marine ecosystems.

Ciliate diversity and distribution across horizontal and vertical scales in the open ocean

Ciliates are globally distributed eukaryotic organisms inhabiting virtually all environments on Earth. Although ciliates range from 10 µm to a few millimetres in cell size, they are repeatedly reported in the pico-sized fraction (<2-3 µm) of molecular surveys. Here, we used existing data sets (BioMarKs and Tara Oceans) with different size fractions to demonstrate that the ciliate pico-sized signal, probably derived from cell breakage during filtration, is informative and reliable to study marine ciliate biodiversity and biogeography. We then used sequences from the pico-eukaryotic fraction of two circumnavigation expeditions, Malaspina-2010 and Tara Oceans, to give insights into the taxonomic composition and horizontal and vertical distribution of ciliates in the global ocean. The results suggested a high homogeneity of ciliate communities along the ocean surface from temperate to tropical waters, with ciliate assemblages dominated by a few abundant and widely distributed taxa. Very few taxa were found in a single oceanic region, therefore suggesting a high level of ciliate cosmopolitanism in the global ocean. In vertical profiles, ciliates were detected up to 4,000 m depth, and a clear vertical community structuring was observed. Our results provided evidence supporting ciliates as deeply integrated organisms in the deep-sea trophic web, where they may play a relevant role as symbionts of metazoans and grazers of prokaryotes and small eukaryotes in the water column and in aggregates.

Does Intraspecific Variation in rDNA Copy Number Affect Analysis of Microbial Communities?

Amplicon sequencing of partial regions of the ribosomal RNA loci (rDNA) is widely used to profile microbial communities. However, the rDNA is dynamic and can exhibit substantial interspecific and intraspecific variation in copy number in prokaryotes and, especially, in microbial eukaryotes. As change in rDNA copy number is a common response to environmental change, rDNA copy number is not necessarily a property of a species. Variation in rDNA copy number, especially the capacity for large intraspecific changes driven by external cues, complicates analyses of rDNA amplicon sequence data. We highlight the need to (i) interpret amplicon sequence data in light of possible interspecific and intraspecific variation, and (ii) examine the potential plasticity in rDNA copy number as an important ecological factor to better understand how microbial communities are structured in heterogeneous environments.

Protist diversity and function in the dark ocean - Challenging the paradigms of deep-sea ecology with special emphasis on foraminiferans and naked protists

The dark ocean and the underlying deep seafloor together represent the largest environment on this planet, comprising about 80% of the oceanic volume and covering more than two-thirds of the Earth's surface, as well as hosting a major part of the total biosphere. Emerging evidence suggests that these vast pelagic and benthic habitats play a major role in ocean biogeochemistry and represent an "untapped reservoir" of high genetic and metabolic microbial diversity. Due to its huge volume, the water column of the dark ocean is the largest reservoir of organic carbon in the biosphere and likely plays a major role in the global carbon budget. The dark ocean and the seafloor beneath it are also home to a largely enigmatic food web comprising little-known and sometimes spectacular organisms, mainly prokaryotes and protists. This review considers the globally important role of pelagic and benthic protists across all protistan size classes in the deep-sea realm, with a focus on their taxonomy, diversity, and physiological properties, including their role in deep microbial food webs. We argue that, given the important contribution that protists must make to deep-sea biodiversity and ecosystem processes, they should not be overlooked in biological studies of the deep ocean.

Insights into spatial effects of ceria nanoparticles on oxygen mass transfer in wastewater biofilms: Interfacial microstructure, in-situ microbial activity and metabolism regulation mechanism

Growing international exploitation of ceria nanoparticles (CeO₂ NPs) for varied applications has increased their release into wastewater treatment plants. Mass transfer of oxygen (MTO) in wastewater biofilm is of considerable importance to influence the activity and purification ability of biofilm. Herein, we investigated the spatial distribution of oxygen in gas-liquid-biofilm phases, the microstructure of interfaces and the in-situ microbial activity to reveal the impacts of CeO₂ NPs on MTO in wastewater biofilm and the related mechanisms. After exposure to 1 and 10 mg/L CeO₂ NPs, the oxygen transfer coefficient (K_La) from gas to wastewater increased by 28.1% and 75.3% with a reduction of thickness in gas-liquid boundary layer, indicating the enhanced MTO in gas-liquid interface. In contrast, the MTO in liquid-biofilm interface was negatively affected and the thickness of liquid-biofilms boundary layer significantly increased, which was mainly attributed to the smoother surface and the decreased surface area difference of biofilm. Within biofilm, the microbial activity was inhibited by 10 mg/L CeO₂ NPs, whereas the production of extracellular polymeric substance (EPS) was significantly improved, leading to a decline of 35.0% in the internal effective diffusivity (D_B) and a 300-μm reduction of oxygen penetration depth. Moreover, the relative activities of key enzymes involved in glycometabolism indicated the transition of Embden-Meyerhof pathway to pentose phosphate pathway, which probably contributed to the enhanced EPS production and consequently increased mass transfer resistance in liquid-biofilm interface and inner biofilm. These results could potentially expand the knowledge on mass transfer of nutrients or pollutants in wastewater biofilm in response to NPs exposure.

High throughput sequencing of 18S rRNA and its gene to characterize a Prorocentrum shikokuense (Dinophyceae) bloom

Comparing phytoplankton non-bloom and bloom communities using rRNA and its coding gene can help understand the shift of dominant species and its driving processes (e.g., intrinsic growth or grazing). Here we conducted high-throughput sequencing of 18S rRNA and its coding gene for studying non-bloom and bloom plankton communities in East China Sea. The non-bloom community was dominated by diatoms whereas during the bloom it was dominated by the dinoflagellate Prorocentrum shikokuense (formerly P. donghaiense). P. shikokuense rRNA abundance and rRNA:rRNA gene ratio both increased markedly in the bloom, indicating that the bloom arose from active growth. In contrast, some non-bloom species showed low DNA abundances during the bloom albeit high rRNA:rRNA gene ratios, suggesting that cell loss processes such as grazing might have prevented these species from blooming or that these species might be at an early stage of bloom development. Furthermore, Pearson's correlation analysis showed that dinoflagellate abundance was positively correlated with temperature and negatively related to dissolved inorganic phosphate (DIP) concentrations, suggesting warm and DIP-poor environment as a niche space for P. shikokuense. Our results demonstrate the usefulness of combined analysis of rRNA and its gene in characterizing phytoplankton bloom development to shed light on the complex phytoplankton dynamics and regulating mechanisms in the course of bloom development.

A metagenomic assessment of microbial eukaryotic diversity in the global ocean

Surveying microbial diversity and function is accomplished by combining complementary molecular tools. Among them, metagenomics is a PCR free approach that contains all genetic information from microbial assemblages and is today performed at a relatively large scale and reasonable cost, mostly based on very short reads. Here, we investigated the potential of metagenomics to provide taxonomic reports of marine microbial eukaryotes. We prepared a curated database with reference sequences of the V4 region of 18S rDNA clustered at 97% similarity and used this database to extract and classify metagenomic reads. More than half of them were unambiguously affiliated to a unique reference whilst the rest could be assigned to a given taxonomic group. The overall diversity reported by metagenomics was similar to that obtained by amplicon sequencing of the V4 and V9 regions of the 18S rRNA gene, although either one or both of these amplicon surveys performed poorly for groups like Excavata, Amoebozoa, Fungi and Haptophyta. We then studied the diversity of picoeukaryotes and nanoeukaryotes using 91 metagenomes from surface down to bathypelagic layers in different oceans, unveiling a clear taxonomic separation between size fractions and depth layers. Finally, we retrieved long rDNA sequences from assembled metagenomes that improved phylogenetic reconstructions of particular groups. Overall, this study shows metagenomics as an excellent resource for taxonomic exploration of marine microbial eukaryotes.