Microbial Eukaryotes Associated With Sediments in Deep-Sea Cold Seeps

Microbial eukaryotes are key components of the marine food web, but their distribution in deep-sea chemosynthetic ecosystems has not been well studied. Here, high-throughput sequencing of the 18S rRNA gene and network analysis were applied to investigate the diversity, distribution and potential relationships between microbial eukaryotes in samples collected from two cold seeps and one trough in the northern South China Sea. SAR (i.e., Stramenopiles, Alveolata, and Rhizaria) was the predominant group in all the samples, and it was highly affiliated to genotypes with potential symbiotic and parasitic strategies identified from other deep-sea extreme environments (e.g., oxygen deficient zones, bathypelagic waters, and hydrothermal vents). Our findings indicated that specialized lineages of deep-sea microbial eukaryotes exist in chemosynthetic cold seeps, where microbial eukaryotes affiliated with parasitic/symbiotic taxa were prevalent in the community. The biogeographic pattern of the total community was best represented by the intermediate operational taxonomic unit (OTU) category, whose relative abundance ranged 0.01–1% within a sample, and the communities of the two cold seeps were distinct from the trough, which suggests that geographical proximity has no critical impact on the distribution of deep-sea microbial eukaryotes. Overall, this study has laid the foundations for future investigations regarding the ecological function and in situ trophic relationships of microbial eukaryotes in deep-sea ecosystems.

Genomic diversity and biosynthetic capabilities of sponge-associated chlamydiae

Sponge microbiomes contribute to host health, nutrition, and defense through the production of secondary metabolites. Chlamydiae, a phylum of obligate intracellular bacteria ranging from animal pathogens to endosymbionts of microbial eukaryotes, are frequently found associated with sponges. However, sponge-associated chlamydial diversity has not yet been investigated at the genomic level and host-interactions remain thus far unexplored. Here, we sequenced the microbiomes of three sponge species and found high, though variable, Chlamydiae relative abundances of up to 21.2% of bacterial diversity. Using genome-resolved metagenomics 18 high-quality sponge-associated chlamydial genomes were reconstructed, covering four chlamydial families. Among these, Sorochlamydiaceae shares a common ancestor with Chlamydiaceae animal pathogens, suggesting long-term co-evolution with animals. Sponge-associated chlamydiae genomes mostly resembled environmental chlamydial endosymbionts, but not pathogens, and encoded genes for degrading diverse compounds associated with sponges, such as taurine. Unexpectedly, we identified widespread genetic potential for secondary metabolite biosynthesis across Chlamydiae, which may represent an explored reservoir of novel natural products. This finding suggests that chlamydiae may partake in defensive symbioses and that secondary metabolites play a wider role in mediating intracellular interactions. Furthermore, sponge-associated chlamydiae relatives were found in other marine invertebrates, pointing towards wider impacts of this phylum on marine ecosystems.

The Diversity Patterns of Rare to Abundant Microbial Eukaryotes Across a Broad Range of Salinities in a Solar Saltern

Solar salterns are excellent artificial systems for examining species diversity and succession along salinity gradients. Here, the eukaryotic community in surface water of a Korean solar saltern (30 to 380 practical salinity units) was investigated from April 2019 to October 2020 using Illumina sequencing targeting the V4 and V9 regions of 18S rDNA. A total of 926 operational taxonomic units (OTUs) and 1,999 OTUs were obtained with the V4 and V9 regions, respectively. Notably, most of the OTUs were microbial eukaryotes, and the high-abundance groups (> 5% relative abundance (RA), Alveolata, Stramenopila, Archaeplastida, and Opisthokonta) usually accounted for > 90% of the total cumulative read counts and > 80% of all OTUs. Moreover, the high-abundance Alveolata (larger forms) and Stramenopila (smaller forms) groups displayed a significant inverse relationship, probably due to predator–prey interactions. Most of the low-abundance (0.1–5% RA) and rare (< 0.1% RA) groups remained small portion during the field surveys. Taxonomic novelty (at < 90% sequence identity) was high in the Amoebozoa, Cryptista, Haptista, Rhizaria, and Stramenopila groups (69.8% of all novel OTUs), suggesting the presence of a large number of hidden species in hypersaline environments. Remarkably, the high-abundance groups had little overlap with the other groups, implying the weakness of rare-to-prevalent community dynamics. The low-abundance Discoba group alone temporarily became the high-abundance group, suggesting that it is an opportunistic group. Overall, the composition and diversity of the eukaryotic community in hypersaline environments may be persistently stabilized, despite diverse disturbance events.

Linking extreme seasonality and gene expression in arctic marine protists

At high latitudes, strong seasonal differences in light availability affect marine organisms and restrict the timing of ecosystem processes. Marine protists are key players in Arctic aquatic ecosystems, yet little is known about their ecological roles over yearly cycles. This is especially true for the dark polar night period, which up until recently was assumed to be devoid of biological activity. A 12 million transcripts catalogue was built from 0.45-10 μm protist assemblages sampled over 13 months in a time series station in an arctic fjord in Svalbard. Community gene expression was correlated with seasonality, with light as the main driving factor. Transcript diversity and evenness were higher during polar night compared to polar day. Light-dependent functions had higher relative expression during polar day, except phototransduction. 64% of the most expressed genes could not be functionally annotated, yet up to 78% were identified in arctic samples from Tara Oceans, suggesting that arctic marine assemblages are distinct from those from other oceans. Our study increases understanding of the links between extreme seasonality and biological processes in pico- and nanoplanktonic protists. Our results set the ground for future monitoring studies investigating the seasonal impact of climate change on the communities of microbial eukaryotes in the High Arctic.

Seasonal Variability of Photosynthetic Microbial Eukaryotes (<3 µm) in the Kara Sea Revealed by 18S rDNA Metabarcoding of Sediment Trap Fluxes

This survey is the first to explore the seasonal cycle of microbial eukaryote diversity (<3 µm) using the NGS method and a 10-month sediment trap (2018–2019). The long-term trap was deployed from September to June in the northwestern part of the Kara Sea. A water sample collected before the sediment trap was deployed and also analyzed. The taxonomic composition of microbial eukaryotes in the water sample significantly differed from sediment trap samples, characterized by a high abundance of Ciliophora reads and low abundance of Fungi while trap samples contained an order of magnitude less Ciliophora sequences and high contribution of Fungi. Photosynthetic eukaryotes (PEs) accounting for about 34% of total protists reads were assigned to five major divisions: Chlorophyta, Cryptophyta, Dinoflagellata, Haptophyta, and Ochrophyta. The domination of phototrophic algae was revealed in late autumn. Mamiellophyceae and Trebouxiophyceae were the predominant PEs in mostly all of the studied seasons. Micromonas polaris was constantly present throughout the September–June period in the PE community. The obtained results determine the seasonal dynamics of picoplankton in order to improve our understanding of their role in polar ecosystems.

Short- and long-read metabarcoding of the eukaryotic rRNA operon: evaluation of primers and comparison to shotgun metagenomics sequencing

High-throughput sequencing for analysis of environmental microbial diversity has evolved vastly over the last decade. Currently the go-to method for microbial eukaryotes is short-read metabarcoding of variable regions of the 18S rRNA gene with <500 bp amplicons. However, there is a growing interest in long-read sequencing of amplicons covering the rRNA operon for improving taxonomic resolution. For both methods, the choice of primers is crucial. It determines if community members are covered, if they can be identified at a satisfactory taxonomic level, and if the obtained community profile is representative. Here, we designed new primers targeting 18S and 28S rRNA based on 177,934 and 21,072 database sequences, respectively. The primers were evaluated in silico along with published primers on reference sequence databases and marine metagenomics datasets. We further evaluated a subset of the primers for short- and long-read sequencing on environmental samples in vitro and compared the obtained community profile with primer-unbiased metagenomic sequencing. Of the short-read pairs, a new V6-V8 pair and the V4_Balzano pair used with a simplified PCR protocol provided good results in silico and in vitro. Fewer differences were observed between the long-read primer pairs. The long-read amplicons and ITS1 alone provided higher taxonomic resolution than V4. Together, our results represent a reference and guide for selection of robust primers for research on and environmental monitoring of microbial eukaryotes.

Estimating the maximal growth rates of eukaryotic microbes from cultures and metagenomes via codon usage patterns

Microbial eukaryotes are ubiquitous in the environment and play important roles in key ecosystem processes, including accounting for a significant portion of global primary production. Yet, our tools for assessing the functional capabilities of eukaryotic microbes in the environment are quite limited because many microbes have yet to be grown in culture. Maximum growth rate is a fundamental parameter of microbial lifestyle that reveals important information about an organism's functional role in a community. We developed and validated a genomic estimator of maximum growth rate for eukaryotic microbes, enabling the assessment of growth potential for both cultivated and yet-to-be-cultivated organisms. We produced a database of over 700 growth predictions from genomes, transcriptomes, and metagenome-assembled genomes, and found that closely related and/or functionally similar organisms tended to have similar maximal growth rates. By comparing the maximal growth rates of existing culture collections with environmentally-derived genomes we found that, unlike for prokaryotes, culture collections of microbial eukaryotes are only minimally biased in terms of growth potential. We then extended our tool to make community-wide estimates of growth potential from over 500 marine metagenomes, mapping growth potential across the global oceans. We found that prokaryotic and eukaryotic communities have highly correlated growth potentials near the ocean surface, but that this relationship disappears deeper in the water column. This suggests that fast growing eukaryotes and prokaryotes thrive under similar conditions at the ocean surface, but that there is a decoupling of these communities as resources become scarce deeper in the water column.

Protist.Guru: A Comparative Transcriptomics Database for the Protist Kingdom

Summary: During the last few decades, the study of microbial ecology has been enabled by molecular and genomic data. DNA sequencing has revealed the surprising extent of microbial diversity and how microbial processes run global ecosystems. However, significant gaps in our understanding of the microbial world remain, and one example is that microbial eukaryotes, or protists, are still neglected. To address this gap, we used gene expression data from 15 distinct protist species to create protist.guru: an online database equipped with tools for identifying functional co-expression networks, gene families, and enriched gene clusters. Here, we show how our database can be used to reveal genes involved in essential pathways, such as the synthesis of secondary carotenoids in Haematococcus lacustris. We expect protist.guru to serve as a valuable resource for protistologists, as well as a catalyst for discoveries and new insights into the biological processes of microbial eukaryotes. Availability: The database and co-expression networks are freely available from http://protist.guru/. The expression matrices and sample annotations are found in the supplementary data.