Single-cell genomics of multiple uncultured stramenopiles reveals underestimated functional diversity across oceans

How marine are Marine Stramenopiles (MAST)? A cross-system evaluation

Marine Stramenopiles (MAST) were first described two decades ago through ribosomal RNA gene (rRNA gene) sequences from marine surveys of microbial eukaryotes. MAST comprise several independent lineages at the base of the Stramenopiles. Despite their prevalence in the ocean, the majority of MAST diversity remains uncultured. Previous studies, mainly in marine environments, have explored MAST's cell morphology, distribution, trophic strategies, and genomics using culturing-independent methods. In comparison, less is known about their presence outside marine habitats. Here, we analyse the extensive EukBank dataset to assess the extent to which MAST can be considered marine protists. Additionally, by incorporating newly available rRNA gene sequences, we update Stramenopiles phylogeny, identifying three novel MAST lineages. Our results indicate that MAST are primarily marine with notable exceptions within MAST-2 and MAST-12, where certain subclades are prevalent in freshwater and soil habitats. In the marine water column, only a few MAST species, particularly within clades -1, -3, -4, and -7, dominate and exhibit clear latitudinal distribution patterns. Overall, the massive sequencing dataset analysed in our study confirms and partially expands the previously described diversity of MASTs groups and underscores the predominantly marine nature of most of these uncultured lineages.

Unveiling the microbial realm with VEBA 2.0: a modular bioinformatics suite for end-to-end genome-resolved prokaryotic, (micro)eukaryotic and viral multi-omics from either short- or long-read sequencing

The microbiome is a complex community of microorganisms, encompassing prokaryotic (bacterial and archaeal), eukaryotic, and viral entities. This microbial ensemble plays a pivotal role in influencing the health and productivity of diverse ecosystems while shaping the web of life. However, many software suites developed to study microbiomes analyze only the prokaryotic community and provide limited to no support for viruses and microeukaryotes. Previously, we introduced the Viral Eukaryotic Bacterial Archaeal (VEBA) open-source software suite to address this critical gap in microbiome research by extending genome-resolved analysis beyond prokaryotes to encompass the understudied realms of eukaryotes and viruses. Here we present VEBA 2.0 with key updates including a comprehensive clustered microeukaryotic protein database, rapid genome/protein-level clustering, bioprospecting, non-coding/organelle gene modeling, genome-resolved taxonomic/pathway profiling, long-read support, and containerization. We demonstrate VEBA's versatile application through the analysis of diverse case studies including marine water, Siberian permafrost, and white-tailed deer lung tissues with the latter showcasing how to identify integrated viruses. VEBA represents a crucial advancement in microbiome research, offering a powerful and accessible software suite that bridges the gap between genomics and biotechnological solutions.

Spatio-temporal changes of small protist and free-living bacterial communities in a temperate dimictic lake: insights from metabarcoding and machine learning

Microbial communities, which include prokaryotes and protists, play an important role in aquatic ecosystems and influence ecological processes. To understand these communities, metabarcoding provides a powerful tool to assess their taxonomic composition and track spatio-temporal dynamics in both marine and freshwater environments. While marine ecosystems have been extensively studied, there is a notable research gap in understanding eukaryotic microbial communities in temperate lakes. Our study addresses this gap by investigating the free-living bacteria and small protist communities in Lake Roś (Poland), a dimictic temperate lake. Metabarcoding analysis revealed that both the bacterial and protist communities exhibit distinct seasonal patterns that are not necessarily shaped by dominant taxa. Furthermore, machine learning and statistical methods identified crucial amplicon sequence variants (ASVs) specific to each season. In addition, we identified a distinct community in the anoxic hypolimnion. We have also shown that the key factors shaping the composition of analysed community are temperature, oxygen, and silicon concentration. Understanding these community structures and the underlying factors is important in the context of climate change potentially impacting mixing patterns and leading to prolonged stratification.

New insights into functional divergence and adaptive evolution of uncultured bacteria in anammox community by complete genome-centric analysis

Anaerobic ammonium-oxidation (anammox) bacteria play a crucial role in global nitrogen cycling and wastewater nitrogen removal, but they share symbiotic relationships with various other microorganisms. Functional divergence and adaptive evolution of uncultured bacteria in anammox community remain underexplored. Although shotgun metagenomics based on short reads has been widely used in anammox research, metagenome-assembled genomes (MAGs) are often discontinuous and highly contaminated, which limits in-depth analyses of anammox communities. Here, for the first time, we performed Pacific Biosciences high-fidelity (HiFi) long-read sequencing on the anammox granule sludge sample from a lab-scale bioreactor, and obtained 30 accurate and complete metagenome-assembled genomes (cMAGs). These cMAGs were obtained by selecting high-quality circular contigs from initial assemblies of long reads generated by HiFi sequencing, eliminating the need for Illumina short reads, binning, and reassembly. One new anammox species affiliated with Candidatus Jettenia and three species affiliated with novel families were found in this anammox community. cMAG-centric analysis revealed functional divergence in general and nitrogen metabolism among the anammox community members, and they might adopt a cross-feeding strategy in organic matter, cofactors, and vitamins. Furthermore, we identified 63 mobile genetic elements (MGEs) and 50 putative horizontal gene transfer (HGT) events within these cMAGs. The results suggest that HGT events and MGEs related to phage and integration or excision, particularly transposons containing tnpA in anammox bacteria, might play important roles in the adaptive evolution of this anammox community. The cMAGs generated in the present study could be used to establish of a comprehensive database for anammox bacteria and associated microorganisms. These findings highlight the advantages of HiFi sequencing for the studies of complex mixed cultures and advance the understanding of anammox communities.

Unveiling the Microbial Realm with VEBA 2.0: A modular bioinformatics suite for end-to-end genome-resolved prokaryotic, (micro)eukaryotic, and viral multi-omics from either short- or long-read sequencing

The microbiome is a complex community of microorganisms, encompassing prokaryotic (bacterial and archaeal), eukaryotic, and viral entities. This microbial ensemble plays a pivotal role in influencing the health and productivity of diverse ecosystems while shaping the web of life. However, many software suites developed to study microbiomes analyze only the prokaryotic community and provide limited to no support for viruses and microeukaryotes. Previously, we introduced the Viral Eukaryotic Bacterial Archaeal (VEBA) open-source software suite to address this critical gap in microbiome research by extending genome-resolved analysis beyond prokaryotes to encompass the understudied realms of eukaryotes and viruses. Here we present VEBA 2.0 with key updates including a comprehensive clustered microeukaryotic protein database, rapid genome/protein-level clustering, bioprospecting, non-coding/organelle gene modeling, genome-resolved taxonomic/pathway profiling, long-read support, and containerization. We demonstrate VEBA's versatile application through the analysis of diverse case studies including marine water, Siberian permafrost, and white-tailed deer lung tissues with the latter showcasing how to identify integrated viruses. VEBA represents a crucial advancement in microbiome research, offering a powerful and accessible platform that bridges the gap between genomics and biotechnological solutions.

Integrated overview of stramenopile ecology, taxonomy, and heterotrophic origin

Stramenopiles represent a significant proportion of aquatic and terrestrial biota. Most biologists can name a few, but these are limited to the phototrophic (e.g. diatoms and kelp) or parasitic species (e.g. oomycetes, Blastocystis), with free-living heterotrophs largely overlooked. Though our attention is slowly turning towards heterotrophs, we have only a limited understanding of their biology due to a lack of cultured models. Recent metagenomic and single-cell investigations have revealed the species richness and ecological importance of stramenopiles-especially heterotrophs. However, our lack of knowledge of the cell biology and behaviour of these organisms leads to our inability to match species to their particular ecological functions. Because photosynthetic stramenopiles are studied independently of their heterotrophic relatives, they are often treated separately in the literature. Here, we present stramenopiles as a unified group with shared synapomorphies and evolutionary history. We introduce the main lineages, describe their important biological and ecological traits, and provide a concise update on the origin of the ochrophyte plastid. We highlight the crucial role of heterotrophs and mixotrophs in our understanding of stramenopiles with the goal of inspiring future investigations in taxonomy and life history. To understand each of the many diversifications within stramenopiles-towards autotrophy, osmotrophy, or parasitism-we must understand the ancestral heterotrophic flagellate from which they each evolved. We hope the following will serve as a primer for new stramenopile researchers or as an integrative refresher to those already in the field.

MarFERReT, an open-source, version-controlled reference library of marine microbial eukaryote functional genes

Metatranscriptomics generates large volumes of sequence data about transcribed genes in natural environments. Taxonomic annotation of these datasets depends on availability of curated reference sequences. For marine microbial eukaryotes, current reference libraries are limited by gaps in sequenced organism diversity and barriers to updating libraries with new sequence data, resulting in taxonomic annotation of about half of eukaryotic environmental transcripts. Here, we introduce Marine Functional EukaRyotic Reference Taxa (MarFERReT), a marine microbial eukaryotic sequence library designed for use with taxonomic annotation of eukaryotic metatranscriptomes. We gathered 902 publicly accessible marine eukaryote genomes and transcriptomes and assessed their sequence quality and cross-contamination issues, selecting 800 validated entries for inclusion in MarFERReT. Version 1.1 of MarFERReT contains reference sequences from 800 marine eukaryotic genomes and transcriptomes, covering 453 species- and strain-level taxa, totaling nearly 28 million protein sequences with associated NCBI and PR2 Taxonomy identifiers and Pfam functional annotations. The MarFERReT project repository hosts containerized build scripts, documentation on installation and use case examples, and information on new versions of MarFERReT.

A k-mer-Based Approach for Phylogenetic Classification of Taxa in Environmental Genomic Data

In the age of genome sequencing, whole-genome data is readily and frequently generated, leading to a wealth of new information that can be used to advance various fields of research. New approaches, such as alignment-free phylogenetic methods that utilize k-mer-based distance scoring, are becoming increasingly popular given their ability to rapidly generate phylogenetic information from whole-genome data. However, these methods have not yet been tested using environmental data, which often tends to be highly fragmented and incomplete. Here, we compare the results of one alignment-free approach (which utilizes the D2 statistic) to traditional multi-gene maximum likelihood trees in 3 algal groups that have high-quality genome data available. In addition, we simulate lower-quality, fragmented genome data using these algae to test method robustness to genome quality and completeness. Finally, we apply the alignment-free approach to environmental metagenome assembled genome data of unclassified Saccharibacteria and Trebouxiophyte algae, and single-cell amplified data from uncultured marine stramenopiles to demonstrate its utility with real datasets. We find that in all instances, the alignment-free method produces phylogenies that are comparable, and often more informative, than those created using the traditional multi-gene approach. The k-mer-based method performs well even when there are significant missing data that include marker genes traditionally used for tree reconstruction. Our results demonstrate the value of alignment-free approaches for classifying novel, often cryptic or rare, species, that may not be culturable or are difficult to access using single-cell methods, but fill important gaps in the tree of life.

Holistic view of the seascape dynamics and environment impact on macro-scale genetic connectivity of marine plankton populations

BACKGROUND

Plankton seascape genomics studies have revealed different trends from large-scale weak differentiation to microscale structures. Previous studies have underlined the influence of the environment and seascape on species differentiation and adaptation. However, these studies have generally focused on a few single species, sparse molecular markers, or local scales. Here, we investigated the genomic differentiation of plankton at the macro-scale in a holistic approach using Tara Oceans metagenomic data together with a reference-free computational method.

RESULTS

We reconstructed the FST-based genomic differentiation of 113 marine planktonic taxa occurring in the North and South Atlantic Oceans, Southern Ocean, and Mediterranean Sea. These taxa belong to various taxonomic clades spanning Metazoa, Chromista, Chlorophyta, Bacteria, and viruses. Globally, population genetic connectivity was significantly higher within oceanic basins and lower in bacteria and unicellular eukaryotes than in zooplankton. Using mixed linear models, we tested six abiotic factors influencing connectivity, including Lagrangian travel time, as proxies of oceanic current effects. We found that oceanic currents were the main population genetic connectivity drivers, together with temperature and salinity. Finally, we classified the 113 taxa into parameter-driven groups and showed that plankton taxa belonging to the same taxonomic rank such as phylum, class or order presented genomic differentiation driven by different environmental factors.

CONCLUSION

Our results validate the isolation-by-current hypothesis for a non-negligible proportion of taxa and highlight the role of other physicochemical parameters in large-scale plankton genetic connectivity. The reference-free approach used in this study offers a new systematic framework to analyse the population genomics of non-model and undocumented marine organisms from a large-scale and holistic point of view.

Gene expression dynamics of natural assemblages of heterotrophic flagellates during bacterivory

BACKGROUND

Marine heterotrophic flagellates (HF) are dominant bacterivores in the ocean, where they represent the trophic link between bacteria and higher trophic levels and participate in the recycling of inorganic nutrients for regenerated primary production. Studying their activity and function in the ecosystem is challenging since most of the HFs in the ocean are still uncultured. In the present work, we investigated gene expression of natural HF communities during bacterivory in four unamended seawater incubations.

RESULTS

The most abundant species growing in our incubations belonged to the taxonomic groups MAST-4, MAST-7, Chrysophyceae, and Telonemia. Gene expression dynamics were similar between incubations and could be divided into three states based on microbial counts, each state displaying distinct expression patterns. The analysis of samples where HF growth was highest revealed some highly expressed genes that could be related to bacterivory. Using available genomic and transcriptomic references, we identified 25 species growing in our incubations and used those to compare the expression levels of these specific genes. Video Abstract CONCLUSIONS: Our results indicate that several peptidases, together with some glycoside hydrolases and glycosyltransferases, are more expressed in phagotrophic than in phototrophic species, and thus could be used to infer the process of bacterivory in natural assemblages.

Genomic differentiation of three pico-phytoplankton species in the Mediterranean Sea

For more than a decade, high-throughput sequencing has transformed the study of marine planktonic communities and has highlighted the extent of protist diversity in these ecosystems. Nevertheless, little is known relative to their genomic diversity at the species-scale as well as their major speciation mechanisms. An increasing number of data obtained from global scale sampling campaigns is becoming publicly available, and we postulate that metagenomic data could contribute to deciphering the processes shaping protist genomic differentiation in the marine realm. As a proof of concept, we developed a findable, accessible, interoperable and reusable (FAIR) pipeline and focused on the Mediterranean Sea to study three a priori abundant protist species: Bathycoccus prasinos, Pelagomonas calceolata and Phaeocystis cordata. We compared the genomic differentiation of each species in light of geographic, environmental and oceanographic distances. We highlighted that isolation-by-environment shapes the genomic differentiation of B. prasinos, whereas P. cordata is impacted by geographic distance (i.e. isolation-by-distance). At present time, the use of metagenomics to accurately estimate the genomic differentiation of protists remains challenging since coverages are lower compared to traditional population surveys. However, our approach sheds light on ecological and evolutionary processes occurring within natural marine populations and paves the way for future protist population metagenomic studies.

Discovering marine biodiversity in the 21st century

We review the current knowledge of the biodiversity of the ocean as well as the levels of decline and threat for species and habitats. The lack of understanding of the distribution of life in the ocean is identified as a significant barrier to restoring its biodiversity and health. We explore why the science of taxonomy has failed to deliver knowledge of what species are present in the ocean, how they are distributed and how they are responding to global and regional to local anthropogenic pressures. This failure prevents nations from meeting their international commitments to conserve marine biodiversity with the results that investment in taxonomy has declined in many countries. We explore a range of new technologies and approaches for discovery of marine species and their detection and monitoring. These include: imaging methods, molecular approaches, active and passive acoustics, the use of interconnected databases and citizen science. Whilst no one method is suitable for discovering or detecting all groups of organisms many are complementary and have been combined to give a more complete picture of biodiversity in marine ecosystems. We conclude that integrated approaches represent the best way forwards for accelerating species discovery, description and biodiversity assessment. Examples of integrated taxonomic approaches are identified from terrestrial ecosystems. Such integrated taxonomic approaches require the adoption of cybertaxonomy approaches and will be boosted by new autonomous sampling platforms and development of machine-speed exchange of digital information between databases.

Genome-centric metagenomics insights into functional divergence and horizontal gene transfer of denitrifying bacteria in anammox consortia

Denitrifying bacteria with high abundances in anammox communities play crucial roles in achieving stable anammox-based systems. Despite the relative constant composition of denitrifying bacteria, their functional diversity remains to be explored in anammox communities. Herein, a total of 77 high-quality metagenome-assembled genomes (MAGs) of denitrifying bacteria were recovered from the anammox community in a full-scale swine wastewater treatment plant. Among these microbes, a total of 26 MAGs were affiliated with the seven dominant denitrifying genera that have total abundances higher than 1%. A meta-analysis of these species suggested that external organics reduced the abundances of genus Ignavibacterium and species MAG.305 of UTPRO2 in anammox communities. Comparative genome analysis revealed functional divergence across different denitrifying bacteria, largely owing to their distinct capabilities for carbohydrate (including endogenous and exogenous) utilization and vitamin (e.g., pantothenate and thiamine) biosynthesis. Serval microbes in this system contained fewer genes encoding biotin, pantothenate and methionine biosynthesis compared with their related species from other habitats. In addition, the genes encoding energy production and conversion (73 genes) and inorganic ion transport (53 genes) putatively transferred from other species to denitrifying bacteria, while these denitrifying bacteria (especially genera UTPRO2 and SCN-69-89) likely donated the genes encoding nutrients (e.g., inorganic ion and amino acid) transporter (64 genes) for other members to utilize new metabolites. Collectively, these findings highlighted the functional divergence of these denitrifying bacteria and speculated that the genetic interactions within anammox communities through horizontal gene transfer may be one of the reasons for their functional divergence.

Genomic evidence for global ocean plankton biogeography shaped by large-scale current systems

Biogeographical studies have traditionally focused on readily visible organisms, but recent technological advances are enabling analyses of the large-scale distribution of microscopic organisms, whose biogeographical patterns have long been debated. Here we assessed the global structure of plankton geography and its relation to the biological, chemical, and physical context of the ocean (the 'seascape') by analyzing metagenomes of plankton communities sampled across oceans during the Tara Oceans expedition, in light of environmental data and ocean current transport. Using a consistent approach across organismal sizes that provides unprecedented resolution to measure changes in genomic composition between communities, we report a pan-ocean, size-dependent plankton biogeography overlying regional heterogeneity. We found robust evidence for a basin-scale impact of transport by ocean currents on plankton biogeography, and on a characteristic timescale of community dynamics going beyond simple seasonality or life history transitions of plankton.

Insights into the microbial community of treated wastewater, its year-round variability and impact on the receiver, using cultivation, microscopy and amplicon-based methods

Apart from chemical constituents, wastewater treatment plant (WWTP) effluents also release microorganisms that can be important to the receiving water bodies either from a sanitary point of view, or taking to the account the biogeochemical potential of the recipients. However, little is known about the treated wastewater microbial community, its composition, seasonal changes, functions and fate in the waters of the receiver. Thus, this study presents a synergistic approach coupling new and traditional methods: analytical chemistry, classical microbiology (cultivation- and microscopy-based methods), as well as Next Generation Sequencing and a quantitative real-time polymerase chain reaction (qPCR). The results show that in terms of bacterial community composition, treated wastewater differed from the environmental samples, irrespectively if they were related or unrelated to the WWTP effluent discharge. The canonical correspondence analysis (CCA) taking into account chemical parameters and taxonomical biodiversity indirectly confirmed the seasonal deterioration of the treated wastewater quality as a result of temperature-driven change of activated sludge community structure and biomass washout (observed also by DAPI staining). Despite seasonal fluctuations of total suspended solids and inter-related parameters (such as COD, BOD, TN, TP), the treated wastewater quality remained within current discharge limits. It was due to treatment processes intensively adjusted by WWTP operators, particularly those necessary to maintain an appropriate rate of autotrophic processes of nitrification and to support biological phosphorus removal. This can explain the observed microbiome composition similarity among WWTP effluents at high taxonomic levels. Obtained data also suggest that besides wastewater treatment efficiency, WWTP effluents are still sources of both human-related microorganisms as well as bacteria equipped in genes involved in N-cycling. Their potential of participation in nutrients cycling in the receivers is widely unknown and require critical attention and better understanding.

Notable predominant morphology of the smallest most abundant protozoa of the open ocean revealed by electron microscopy

In the microbe-driven ecosystems of the open ocean, the small heterotrophic flagellates (sHF) are the chief microbial predators and recyclers of essential nutrients to phototrophic microbes. Even with intensive molecular phylogenetic studies of the sHF, the origins of their feeding success remain obscure because of limited understanding of their morphological adaptations to feeding. Here, we examined the sHF morphologies in the largest, most oligotrophic South Pacific and Atlantic (sub)tropical gyres and adjacent mesotrophic waters. On four research cruises, the sHF cells were flow cytometrically sorted from bacterioplankton and phytoplankton for electron microscopy. The sorted sHF comprised chiefly heterokont (HK) biflagellates and unikont choanoflagellates numerically at around 10-to-1 ratio. Of the four differentiated morphological types of HK omnipresent in the open ocean, the short-tinsel heterokont (stHK), whose tinsel flagellum is too short to propagate a complete wave, is predominant and a likely candidate to be the most abundant predator on Earth. Modeling shows that the described stHK propulsion is effective in feeding on bacterioplankton cells at low concentrations; however, owing to general prey scarcity in the oligotrophic ocean, selective feeding is unsustainable and omnivory is equally obligatory for the seven examined sHF types irrespective of their mode of propulsion.

How Communities of Marine Stramenopiles Varied with Environmental and Biological Variables in the Subtropical Northwestern Pacific Ocean

MArine STramenopiles (MASTs) have been recognized as parts of heterotrophic protists and contribute substantially to protist abundances in the ocean. However, little is known about their spatiotemporal variations with respect to environmental and biological factors. The objectives of this study are to use canonical correspondence analysis to investigate how MASTs communities are shaped by environmental variables, and co-occurrence networks to examine their potential interactions with prokaryotic communities. Our dataset came from the southern East China Sea (sECS) in the subtropical northwestern Pacific, and involved 14 cruises along a coastal-oceanic transect, each of which sampled surface water from 4 to 7 stations. MASTs communities were revealed by metabarcoding of 18S rDNA V4 region. Most notably, MAST-9 had a high representation in warm waters in terms of read number and diversity. Subclades of MAST-9C and -9D showed slightly different niches, with MAST-9D dominating in more coastal waters where concentrations of nitrite and Synechococcus were higher. MAST-1C was a common component of colder water during spring. Overall, canonical correspondence analysis showed that MASTs communities were significantly influenced by temperature, nitrite and Synechococcus concentrations. The co-occurrence networks showed that certain other minor prokaryotic taxa can influence MAST communities. This study provides insight into how MASTs communities varied with environmental and biological variables.

Evolution of Phytoplankton as Estimated from Genetic Diversity

Phytoplankton are photosynthetic, single-celled organisms producing almost half of all oxygen on Earth and play a central role as prey for higher organisms, making them irreplaceable in the marine food web. As Global Change proceeds, imposing rapidly intensifying selection pressures, phytoplankton are forced to undergo evolution, local extinction, or redistribution, with potentially cascading effects throughout the marine ecosystem. Recent results from the field of population genetics display high levels of standing genetic diversity in natural phytoplankton populations, providing ample ‘evolutionary options’ and implying high adaptive potential to changing conditions. This potential for adaptive evolution is realized in several studies of experimental evolution, even though most of these studies investigate the evolution of only single strains. This, however, shows that phytoplankton not only evolve from standing genetic diversity, but also rely on de novo mutations. Recent global sampling campaigns show that the immense intraspecific diversity of phytoplankton in the marine ecosystem has been significantly underestimated, meaning we are only studying a minor portion of the relevant variability in the context of Global Change and evolution. An increased understanding of genomic diversity is primarily hampered by the low number of ecologically representative reference genomes of eukaryotic phytoplankton and the functional annotation of these. However, emerging technologies relying on metagenome and transcriptome data may offer a more realistic understanding of phytoplankton diversity.

Freshwater protists: unveiling the unexplored in a large floodplain system

Protists play a fundamental role in all ecosystems, but we are still far from estimating the total diversity of many lineages, in particular in highly diverse environments, such as freshwater. Here, we survey the protist diversity of the Paraná River using metabarcoding, and we applied an approach that includes sequence similarity and phylogeny to evaluate the degree of genetic novelty of the protists' communities against the sequences described in the reference database PR² . We observed that ~28% of the amplicon sequence variants were classified as novel according to their similarity with sequences from the reference database; most of them were related to heterotrophic groups traditionally overlooked in freshwater systems. This lack of knowledge extended to those groups within the green algae (Archaeplastida) that are well documented such as Mamiellophyceae, and also to the less studied Pedinophyceae, for which we found sequences representing novel deep-branching clusters. Among the groups with potential novel protists, Bicosoecida (Stramenopiles) were the best represented, followed by Codosiga (Opisthokonta), and the Perkinsea (Alveolata). This illustrates the lack of knowledge on freshwater planktonic protists and also the need for isolation and/or cultivation of new organisms to better understand their role in ecosystem functioning.

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.

Niche adaptation promoted the evolutionary diversification of tiny ocean predators

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

Genome sequencing and de novo assembly of the giant unicellular alga Acetabularia acetabulum using droplet MDA

The macroscopic single-celled green alga Acetabularia acetabulum has been a model system in cell biology for more than a century. However, no genomic information is available from this species. Since the alga has a long life cycle, is difficult to grow in dense cultures, and has an estimated diploid genome size of almost 2 Gb, obtaining sufficient genomic material for genome sequencing is challenging. Here, we have attempted to overcome these challenges by amplifying genomic DNA using multiple displacement amplification (MDA) combined with microfluidics technology to distribute the amplification reactions across thousands of microscopic droplets. By amplifying and sequencing DNA from five single cells we were able to recover an estimated ~ 7–11% of the total genome, providing the first draft of the A. acetabulum genome. We highlight challenges associated with genome recovery and assembly of MDA data due to biases arising during genome amplification, and hope that our study can serve as a reference for future attempts on sequencing the genome from non-model eukaryotes.

Comparative genomics reveals new functional insights in uncultured MAST species

Heterotrophic lineages of stramenopiles exhibit enormous diversity in morphology, lifestyle, and habitat. Among them, the marine stramenopiles (MASTs) represent numerous independent lineages that are only known from environmental sequences retrieved from marine samples. The core energy metabolism characterizing these unicellular eukaryotes is poorly understood. Here, we used single-cell genomics to retrieve, annotate, and compare the genomes of 15 MAST species, obtained by coassembling sequences from 140 individual cells sampled from the marine surface plankton. Functional annotations from their gene repertoires are compatible with all of them being phagocytotic. The unique presence of rhodopsin genes in MAST species, together with their widespread expression in oceanic waters, supports the idea that MASTs may be capable of using sunlight to thrive in the photic ocean. Additional subsets of genes used in phagocytosis, such as proton pumps for vacuole acidification and peptidases for prey digestion, did not reveal particular trends in MAST genomes as compared with nonphagocytotic stramenopiles, except a larger presence and diversity of V-PPase genes. Our analysis reflects the complexity of phagocytosis machinery in microbial eukaryotes, which contrasts with the well-defined set of genes for photosynthesis. These new genomic data provide the essential framework to study ecophysiology of uncultured species and to gain better understanding of the function of rhodopsins and related carotenoids in stramenopiles.

Current State and Future Directions of Genetics and Genomics of Endophytic Fungi for Bioprospecting Efforts

The bioprospecting of secondary metabolites from endophytic fungi received great attention in the 1990s and 2000s, when the controversy around taxol production from Taxus spp. endophytes was at its height. Since then, hundreds of reports have described the isolation and characterization of putative secondary metabolites from endophytic fungi. However, only very few studies also report the genetic basis for these phenotypic observations. With low sequencing cost and fast sample turnaround, genetics- and genomics-based approaches have risen to become comprehensive approaches to study natural products from a wide-range of organisms, especially to elucidate underlying biosynthetic pathways. However, in the field of fungal endophyte biology, elucidation of biosynthetic pathways is still a major challenge. As a relatively poorly investigated group of microorganisms, even in the light of recent efforts to sequence more fungal genomes, such as the 1000 Fungal Genomes Project at the Joint Genome Institute (JGI), the basis for bioprospecting of enzymes and pathways from endophytic fungi is still rather slim. In this review we want to discuss the current approaches and tools used to associate phenotype and genotype to elucidate biosynthetic pathways of secondary metabolites in endophytic fungi through the lens of bioprospecting. This review will point out the reported successes and shortcomings, and discuss future directions in sampling, and genetics and genomics of endophytic fungi. Identifying responsible biosynthetic genes for the numerous secondary metabolites isolated from endophytic fungi opens the opportunity to explore the genetic potential of producer strains to discover novel secondary metabolites and enhance secondary metabolite production by metabolic engineering resulting in novel and more affordable medicines and food additives.

CARD-FISH in the Sequencing Era: Opening a New Universe of Protistan Ecology

Phagotrophic protists are key players in aquatic food webs. Although sequencing-based studies have revealed their enormous diversity, ecological information on in situ abundance, feeding modes, grazing preferences, and growth rates of specific lineages can be reliably obtained only using microscopy-based molecular methods, such as Catalyzed Reporter Deposition-Fluorescence in situ Hybridization (CARD-FISH). CARD-FISH is commonly applied to study prokaryotes, but less so to microbial eukaryotes. Application of this technique revealed that Paraphysomonas or Spumella-like chrysophytes, considered to be among the most prominent members of protistan communities in pelagic environments, are omnipresent but actually less abundant than expected, in contrast to little known groups such as heterotrophic cryptophyte lineages (e.g., CRY1), cercozoans, katablepharids, or the MAST lineages. Combination of CARD-FISH with tracer techniques and application of double CARD-FISH allow visualization of food vacuole contents of specific flagellate groups, thus considerably challenging our current, simplistic view that they are predominantly bacterivores. Experimental manipulations with natural communities revealed that larger flagellates are actually omnivores ingesting both prokaryotes and other protists. These new findings justify our proposition of an updated model of microbial food webs in pelagic environments, reflecting more authentically the complex trophic interactions and specific roles of flagellated protists, with inclusion of at least two additional trophic levels in the nanoplankton size fraction. Moreover, we provide a detailed CARD-FISH protocol for protists, exemplified on mixo- and heterotrophic nanoplanktonic flagellates, together with tips on probe design, a troubleshooting guide addressing most frequent obstacles, and an exhaustive list of published probes targeting protists.

DMSP synthesis genes distinguish two types of DMSP producer phenotypes

Dimethylsulfoniopropionate (DMSP) is an important organic carbon and sulfur source in the surface ocean that fuels microbial activity and significantly impacts Earth's climate. After three decades of research, the cellular role(s) of DMSP and environmental drivers of production remain enigmatic. Recent work suggests that cellular DMSP concentrations, and changes in these concentrations in response to environmental stressors, define two major groups of DMSP producers: high DMSP producers that contain ≥ 50 mM intracellular DMSP and low DMSP producers that contain < 50 mM. Here we show that two recently described DMSP synthesis genes (DSYB and TpMT2) may differentiate these two DMSP phenotypes. A survey of prokaryotic and eukaryotic isolates found a significant correlation between the presence of DSYB and TpMT2 genes and previous measurements of high and low DMSP concentrations, respectively. Phylogenetic analysis demonstrated that DSYB and TpMT2 form two distinct clades. DSYB and TpMT2 were also found to be globally abundant in in situ surface communities, and their taxonomic annotations were similar to those observed for isolates. The strong correlation of the DSYB and TpMT2 synthesis genes with high and low producer phenotypes establishes a foundation for direct quantification of DMSP producers, enabling significantly improved predictions of DMSP in situ.

Benchmarking of single-virus genomics: a new tool for uncovering the virosphere

Metagenomics and single-cell genomics have enabled the discovery of relevant uncultured microbes. Recently, single-virus genomics (SVG), although still in an incipient stage, has opened new avenues in viral ecology by allowing the sequencing of one single virus at a time. The investigation of methodological alternatives and optimization of existing procedures for SVG is paramount to deliver high-quality genomic data. We report a sequencing dataset of viral single-amplified genomes (vSAGs) from cultured and uncultured viruses obtained by applying different conditions in each SVG step, from viral preservation and novel whole-genome amplification (WGA) to sequencing platforms and genome assembly. Sequencing data showed that cryopreservation and mild fixation were compatible with WGA, although fresh samples delivered better genome quality data. The novel TruPrime WGA, based on primase-polymerase features, and WGA-X employing a thermostable phi29 polymerase, were proven to be with sufficient sensitivity in SVG. The Oxford Nanopore (ON) sequencing platform did not provide a significant improvement of vSAG assembly compared to Illumina alone. Finally, the SPAdes assembler performed the best. Overall, our results represent a valuable genomic dataset that will help to standardized and advance new tools in viral ecology.

Protist taxonomic and functional diversity in soil, freshwater and marine ecosystems

Protists dominate eukaryotic diversity and play key functional roles in all ecosystems, particularly by catalyzing carbon and nutrient cycling. To date, however, a comparative analysis of their taxonomic and functional diversity that compares the major ecosystems on Earth (soil, freshwater and marine systems) is missing. Here, we present a comparison of protist diversity based on standardized high throughput 18S rRNA gene sequencing of soil, freshwater and marine environmental DNA. Soil and freshwater protist communities were more similar to each other than to marine protist communities, with virtually no overlap of Operational Taxonomic Units (OTUs) between terrestrial and marine habitats. Soil protists showed higher γ diversity than aquatic samples. Differences in taxonomic composition of the communities led to changes in a functional diversity among ecosystems, as expressed in relative abundance of consumers, phototrophs and parasites. Phototrophs (eukaryotic algae) dominated freshwater systems (49% of the sequences) and consumers soil and marine ecosystems (59% and 48%, respectively). The individual functional groups were composed of ecosystem- specific taxonomic groups. Parasites were equally common in all ecosystems, yet, terrestrial systems hosted more OTUs assigned to parasites of macro-organisms while aquatic systems contained mostly microbial parasitoids. Together, we show biogeographic patterns of protist diversity across major ecosystems on Earth, preparing the way for more focused studies that will help understanding the multiple roles of protists in the biosphere.

Divergent Microbiota Dynamics along the Coastal Marine Ecosystem of Puerto Rico

Understanding the different factors shaping the spatial and temporal distribution of marine microorganisms is fundamental in predicting their responses to future environmental disturbances. There has been, however, little effort to characterize the microbial diversity including the microbiome dynamics among regions in the Caribbean Sea. Toward this end, this study was designed to gain some critical insights into microbial diversity within the coastal marine ecosystem off the coast of Puerto Rico. Using Illumina MiSeq, the V4 region of the 16S rRNA gene was sequenced with the goal of characterizing the microbial diversity representative of different coastal sites around the island of Puerto Rico. This study provided valuable insights in terms of the local bacterial taxonomic abundance, α and β diversity, and the environmental factors shaping microbial community composition and structure. The most dominant phyla across all 11 sampling sites were the Proteobacteria, Bacteroidetes, and Planctomycetes, while the least dominant taxonomic groups were the NKB19, Tenericutes, OP3, Lentisphaerae, and SAR406. The geographical area (Caribbean and Atlantic seas) and salinity gradients were the main drivers shaping the marine microbial community around the island. Despite stable physical and chemical features of the different sites, a highly dynamic microbiome was observed. This highlights Caribbean waters as one of the richest marine sources for a microbial biodiversity hotspot. The data presented here provide a basis for further temporal evaluations aiming at deciphering microbial taxonomic diversity around the island, while determining how microbes adapt to changes in the climate.

Prey capture in protists utilizing microtubule filled processes and surface motility

Surface motility, which can be visualized by the movement of live prey organisms, polystyrene microspheres or other inert particles, has been shown to occur in a wide variety of microtubule-filled extensions of the protistan cell surface, although the associated functions remain enigmatic. This article integrates an extensive but poorly known body of literature showing that surface motility, associated with microtubule-filled cell extensions such as flagella, axopodia, actinopodia, reticulopodia, and haptonema, plays a crucial role in protistan prey capture. Surface motility has been most extensively studied in Chlamydomonas where it is responsible for flagella-dependent whole cell gliding motility. The force transduction machinery for gliding motility in Chlamydomonas is intraflagellar transport. Other than in Chlamydomonas, this field has not moved far beyond the descriptive to the mechanistic because of technical challenges associated with many of the protistan organisms that utilize surface motility for prey capture. The purpose of this article is to rekindle interest in the protistan systems that utilize surface motility for prey capture at a time when newly emerging molecular tools for working with protists are poised to reinvigorate a field that has been quiescent too long.

Single-virus genomics and beyond

Viruses are extremely diverse and modulate important biological and ecological processes globally. However, much of viral diversity remains uncultured and yet to be discovered. Several powerful culture-independent tools, in particular metagenomics, have substantially advanced virus discovery. Among those tools is single-virus genomics, which yields sequenced reference genomes from individual sorted virus particles without the need for cultivation. This new method complements virus culturing and metagenomic approaches and its advantages include targeted investigation of specific virus groups and investigation of genomic microdiversity within viral populations. In this Review, we provide a brief history of single-virus genomics, outline how this emergent method has facilitated advances in virus ecology and discuss its current limitations and future potential. Finally, we address how this method may synergistically intersect with other single-virus and single-cell approaches.

Tara Oceans: towards global ocean ecosystems biology

A planetary-scale understanding of the ocean ecosystem, particularly in light of climate change, is crucial. Here, we review the work of Tara Oceans, an international, multidisciplinary project to assess the complexity of ocean life across comprehensive taxonomic and spatial scales. Using a modified sailing boat, the team sampled plankton at 210 globally distributed sites at depths down to 1,000 m. We describe publicly available resources of molecular, morphological and environmental data, and discuss how an ecosystems biology approach has expanded our understanding of plankton diversity and ecology in the ocean as a planetary, interconnected ecosystem. These efforts illustrate how global-scale concepts and data can help to integrate biological complexity into models and serve as a baseline for assessing ecosystem changes and the future habitability of our planet in the Anthropocene epoch.

Reactivity, fate and functional roles of dissolved organic matter in anoxic inland waters

The transit of organic matter (OM) through the aquatic compartment of its global cycle has been intensively studied, traditionally with a focus on the processing and degradation of its dissolved fraction (dissolved organic matter, DOM). Because this is so intimately related to oxidation, the notion tenaciously persists that where oxygen is absent, DOM turnover is markedly slowed. In this Opinion Piece, we outline how diverse processes shape, transform and degrade DOM also in anoxic aquatic environments, and we focus here on inland waters as a particular case study. A suite of biogeochemical DOM functions that have received comparatively little attention may only be expressed in anoxic conditions and may result in enhanced biogeochemical roles of these deoxygenated habitats on a network scale.

Four high-quality draft genome assemblies of the marine heterotrophic nanoflagellate Cafeteria roenbergensis

The heterotrophic stramenopile Cafeteria roenbergensis is a globally distributed marine bacterivorous protist. This unicellular flagellate is host to the giant DNA virus CroV and the virophage mavirus. We sequenced the genomes of four cultured C. roenbergensis strains and generated 23.53 Gb of Illumina MiSeq data (99–282 × coverage per strain) and 5.09 Gb of PacBio RSII data (13–45 × coverage). Using the Canu assembler and customized curation procedures, we obtained high-quality draft genome assemblies with a total length of 34–36 Mbp per strain and contig N50 lengths of 148 kbp to 464 kbp. The C. roenbergensis genome has a GC content of ~70%, a repeat content of ~28%, and is predicted to contain approximately 7857–8483 protein-coding genes based on a combination of de novo, homology-based and transcriptome-supported annotation. These first high-quality genome assemblies of a bicosoecid fill an important gap in sequenced stramenopile representatives and enable a more detailed evolutionary analysis of heterotrophic protists.

Measurement(s)	DNA • mitochondrial_DNA • sequence_assembly • sequence feature annotation
Technology Type(s)	DNA sequencing • genome assembly • sequence annotation
Sample Characteristic - Organism	Cafeteria roenbergensis
Sample Characteristic - Environment	marine water body
Sample Characteristic - Location	Northwest Atlantic Ocean • Carribean Sea • Southeast Pacific Ocean • North East Pacific Ocean

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.11419098

Genome Resolved Biogeography of Mamiellales

Among marine phytoplankton, Mamiellales encompass several species from the genera Micromonas, Ostreococcus and Bathycoccus, which are important contributors to primary production. Previous studies based on single gene markers described their wide geographical distribution but led to discussion because of the uneven taxonomic resolution of the method. Here, we leverage genome sequences for six Mamiellales species, two from each genus Micromonas, Ostreococcus and Bathycoccus, to investigate their distribution across 133 stations sampled during the Tara Oceans expedition. Our study confirms the cosmopolitan distribution of Mamiellales and further suggests non-random distribution of species, with two triplets of co-occurring genomes associated with different temperatures: Ostreococcuslucimarinus, Bathycoccusprasinos and Micromonaspusilla were found in colder waters, whereas Ostreococcus spp. RCC809, Bathycoccus spp. TOSAG39-1 and Micromonascommoda were more abundant in warmer conditions. We also report the distribution of the two candidate mating-types of Ostreococcus for which the frequency of sexual reproduction was previously assumed to be very low. Indeed, both mating types were systematically detected together in agreement with either frequent sexual reproduction or the high prevalence of a diploid stage. Altogether, these analyses provide novel insights into Mamiellales' biogeography and raise novel testable hypotheses about their life cycle and ecology.

Phytoplankton in the Tara Ocean

Photosynthesis evolved in the ocean more than 2 billion years ago and is now performed by a wide range of evolutionarily distinct organisms, including both prokaryotes and eukaryotes. Our appreciation of their abundance, distributions, and contributions to primary production in the ocean has been increasing since they were first discovered in the seventeenth century and has now been enhanced by data emerging from the Tara Oceans project, which performed a comprehensive worldwide sampling of plankton in the upper layers of the ocean between 2009 and 2013. Largely using recent data from Tara Oceans, here we review the geographic distributions of phytoplankton in the global ocean and their diversity, abundance, and standing stock biomass. We also discuss how omics-based information can be incorporated into studies of photosynthesis in the ocean and show the likely importance of mixotrophs and photosymbionts.

Charting the Complexity of the Marine Microbiome through Single-Cell Genomics

Marine bacteria and archaea play key roles in global biogeochemistry. To improve our understanding of this complex microbiome, we employed single-cell genomics and a randomized, hypothesis-agnostic cell selection strategy to recover 12,715 partial genomes from the tropical and subtropical euphotic ocean. A substantial fraction of known prokaryoplankton coding potential was recovered from a single, 0.4 mL ocean sample, which indicates that genomic information disperses effectively across the globe. Yet, we found each genome to be unique, implying limited clonality within prokaryoplankton populations. Light harvesting and secondary metabolite biosynthetic pathways were numerous across lineages, highlighting the value of single-cell genomics to advance the identification of ecological roles and biotechnology potential of uncultured microbial groups. This genome collection enabled functional annotation and genus-level taxonomic assignments for >80% of individual metagenome reads from the tropical and subtropical surface ocean, thus offering a model to improve reference genome databases for complex microbiomes.

Unexpected mitochondrial genome diversity revealed by targeted single-cell genomics of heterotrophic flagellated protists

Most eukaryotic microbial diversity is uncultivated, under-studied and lacks nuclear genome data. Mitochondrial genome sampling is more comprehensive, but many phylogenetically important groups remain unsampled. Here, using a single-cell sorting approach combining tubulin-specific labelling with photopigment exclusion, we sorted flagellated heterotrophic unicellular eukaryotes from Pacific Ocean samples. We recovered 206 single amplified genomes, predominantly from underrepresented branches on the tree of life. Seventy single amplified genomes contained unique mitochondrial contigs, including 21 complete or near-complete mitochondrial genomes from formerly under-sampled phylogenetic branches, including telonemids, katablepharids, cercozoans and marine stramenopiles, effectively doubling the number of available samples of heterotrophic flagellate mitochondrial genomes. Collectively, these data identify a dynamic history of mitochondrial genome evolution including intron gain and loss, extensive patterns of genetic code variation and complex patterns of gene loss. Surprisingly, we found that stramenopile mitochondrial content is highly plastic, resembling patterns of variation previously observed only in plants.

Annual Protist Community Dynamics in a Freshwater Ecosystem Undergoing Contrasted Climatic Conditions: The Saint-Charles River (Canada)

Protists are key stone components of aquatic ecosystems, sustaining primary productivity and aquatic food webs. However, their diversity, ecology and structuring factors shaping their temporal distribution remain strongly misunderstood in freshwaters. Using high-throughput sequencing on water samples collected over 16 different months (including two summer and two winter periods), combined with geochemical measurements and climate monitoring, we comprehensively determined the pico- and nanoeukaryotic community composition and dynamics in a Canadian river undergoing prolonged ice-cover winters. Our analysis revealed a large protist diversity in this fluctuating ecosystem and clear seasonal patterns demonstrating a direct and/or indirect selective role of abiotic factors, such as water temperature or nitrogen concentrations, in structuring the eukaryotic microbial community. Nonetheless, our results also revealed that primary productivity, predatory as well as parasitism lifestyles, inferred from fine phylogenetic placements, remained potentially present over the annual cycle, despite the large seasonal fluctuations and the remodeling of the community composition under ice. In addition, potential interplays with the bacterial community composition were identified supporting a possible contribution of the bacterial community to the temporal dynamics of the protist community structure. Our results illustrate the complexity of the eukaryotic microbial community and provide a substantive and useful dataset to better understand the global freshwater ecosystem functioning.

The New Tree of Eukaryotes

For 15 years, the eukaryote Tree of Life (eToL) has been divided into five to eight major groupings, known as 'supergroups'. However, the tree has been profoundly rearranged during this time. The new eToL results from the widespread application of phylogenomics and numerous discoveries of major lineages of eukaryotes, mostly free-living heterotrophic protists. The evidence that supports the tree has transitioned from a synthesis of molecular phylogenetics and biological characters to purely molecular phylogenetics. Most current supergroups lack defining morphological or cell-biological characteristics, making the supergroup label even more arbitrary than before. Going forward, the combination of traditional culturing with maturing culture-free approaches and phylogenomics should accelerate the process of completing and resolving the eToL at its deepest levels.

A single-cell genome perspective on studying intracellular associations in unicellular eukaryotes

Single-cell genomics (SCG) methods provide a unique opportunity to analyse whole genome information at the resolution of an individual cell. While SCG has been extensively used to investigate bacterial and archaeal genomes, the technique has been rarely used to access the genetic makeup of uncultivated microbial eukaryotes. In this regard, the use of SCG can provide a wealth of information; not only do the methods allow exploration of the genome, they can also help elucidate the relationship between the cell and intracellular entities extant in nearly all eukaryotes. SCG enables the study of total eukaryotic cellular DNA, which in turn allows us to better understand the evolutionary history and diversity of life, and the physiological interactions that define complex organisms.

This article is part of a discussion meeting issue ‘Single cell ecology’.

Single cell ecology

Cells are the building blocks of life, from single-celled microbes through to multi-cellular organisms. To understand a multitude of biological processes we need to understand how cells behave, how they interact with each other and how they respond to their environment. The use of new methodologies is changing the way we study cells allowing us to study them on minute scales and in unprecedented detail. These same methods are allowing researchers to begin to sample the vast diversity of microbes that dominate natural environments. The aim of this special issue is to bring together research and perspectives on the application of new approaches to understand the biological properties of cells, including how they interact with other biological entities. This article is part of a discussion meeting issue 'Single cell ecology'.

Single cell ecogenomics reveals mating types of individual cells and ssDNA viral infections in the smallest photosynthetic eukaryotes

Planktonic photosynthetic organisms of the class Mamiellophyceae include the smallest eukaryotes (less than 2 µm), are globally distributed and form the basis of coastal marine ecosystems. Eight complete fully annotated 13-22 Mb genomes from three genera, Ostreococcus, Bathycoccus and Micromonas, are available from previously isolated clonal cultured strains and provide an ideal resource to explore the scope and challenges of analysing single cell amplified genomes (SAGs) isolated from a natural environment. We assembled data from 12 SAGs sampled during the Tara Oceans expedition to gain biological insights about their in situ ecology, which might be lost by isolation and strain culture. Although the assembled nuclear genomes were incomplete, they were large enough to infer the mating types of four Ostreococcus SAGs. The systematic occurrence of sequences from the mitochondria and chloroplast, representing less than 3% of the total cell's DNA, intimates that SAGs provide suitable substrates for detection of non-target sequences, such as those of virions. Analysis of the non-Mamiellophyceae assemblies, following filtering out cross-contaminations during the sequencing process, revealed two novel 1.6 and 1.8 kb circular DNA viruses, and the presence of specific Bacterial and Oomycete sequences suggests that these organisms might co-occur with the Mamiellales. This article is part of a discussion meeting issue 'Single cell ecology'.

Combining morphology, behaviour and genomics to understand the evolution and ecology of microbial eukaryotes

Microbial eukaryotes (protists) are structurally, developmentally and behaviourally more complex than their prokaryotic cousins. This complexity makes it more difficult to translate genomic and metagenomic data into accurate functional inferences about systems ranging all the way from molecular and cellular levels to global ecological networks. This problem can be traced back to the advent of the cytoskeleton and endomembrane systems at the origin of eukaryotes, which endowed them with a range of complex structures and behaviours that still largely dominate how they evolve and interact within microbial communities. But unlike the diverse metabolic properties that evolved within prokaryotes, the structural and behavioural characteristics that strongly define how protists function in the environment cannot readily be inferred from genomic data, since there is generally no simple correlation between a gene and a discrete activity or function. A deeper understanding of protists at both cellular and ecological levels, therefore, requires not only high-throughput genomics but also linking such data to direct observations of natural history and cell biology. This is challenging since these observations typically require cultivation, which is lacking for most protists. Potential remedies with current technology include developing a more phylogenetically diverse range of model systems to better represent the diversity, as well as combining high-throughput, single-cell genomics with microscopic documentation of the subject cells to link sequence with structure and behaviour. This article is part of a discussion meeting issue 'Single cell ecology'.

Phylogenetic Estimation of Community Composition and Novel Eukaryotic Lineages in Base Mine Lake: An Oil Sands Tailings Reclamation Site in Northern Alberta

Reclamation of anthropogenically impacted environments is a critical issue worldwide. In the oil sands extraction industry of Alberta, reclamation of mining-impacted areas, especially areas affected by tailings waste, is an important aspect of the mining life cycle. A reclamation technique currently under study is water-capping, where tailings are capped by water to create an end-pit lake (EPL). Base Mine Lake (BML) is the first full-scale end-pit lake in the Alberta oil sands region. In this study, we sequenced eukaryotic 18S rRNA genes recovered from 92 samples of Base Mine Lake water in a comprehensive sampling programme covering the ice-free period of 2015. The 565 operational taxonomic units (OTUs) generated revealed a dynamic and diverse community including abundant Microsporidia, Ciliata and Cercozoa, though 41% of OTUs were not classifiable below the phylum level by comparison to 18S rRNA databases. Phylogenetic analysis of five heterotrophic phyla (Cercozoa, Fungi, Ciliata, Amoebozoa and Excavata) revealed substantial novel diversity, with many clusters of OTUs that were more similar to each other than to any reference sequence. All of these groups are entirely or mostly heterotrophic, as a relatively small number of definitively photosynthetic clades were amplified from the BML samples.

Morphology, Ultrastructure, and Mitochondrial Genome of the Marine Non-Photosynthetic Bicosoecid Cafileria marina Gen. et sp. nov

In this paper, we describe a novel bacteriophagous biflagellate, Cafileria marina with two smooth flagellae, isolated from material collected from a rock surface in the Kvernesfjorden (Norway). This flagellate was characterized by scanning and transmission electron microscopy, fluorescence, and light microscopy. The sequence of the small subunit ribosomal RNA gene (18S) was used as a molecular marker for determining the phylogenetic position of this organism. Apart from the nuclear ribosomal gene, the whole mitochondrial genome was sequenced, assembled, and annotated. Morphological observations show that the newly described flagellate shares key ultrastructural characters with representatives of the family Bicosoecida (Heterokonta). Intriguingly, mitochondria of C. marina frequently associate with its nucleus through an electron-dense disc at the boundary of the two compartments. The function of this association remains unclear. Phylogenetic analyses corroborate the morphological data and place C. marina with other sequence data of representatives from the family Bicosoecida. We describe C. marina as a new species from a new genus in this family.

Single cell genomics yields a wide diversity of small planktonic protists across major ocean ecosystems

Marine planktonic protists are critical components of ocean ecosystems and are highly diverse. Molecular sequencing methods are being used to describe this diversity and reveal new associations and metabolisms that are important to how these ecosystems function. We describe here the use of the single cell genomics approach to sample and interrogate the diversity of the smaller (pico- and nano-sized) protists from a range of oceanic samples. We created over 900 single amplified genomes (SAGs) from 8 Tara Ocean samples across the Indian Ocean and the Mediterranean Sea. We show that flow cytometric sorting of single cells effectively distinguishes plastidic and aplastidic cell types that agree with our understanding of protist phylogeny. Yields of genomic DNA with PCR-identifiable 18S rRNA gene sequence from single cells was low (15% of aplastidic cell sorts, and 7% of plastidic sorts) and tests with alternate primers and comparisons to metabarcoding did not reveal phylogenetic bias in the major protist groups. There was little evidence of significant bias against or in favor of any phylogenetic group expected or known to be present. The four open ocean stations in the Indian Ocean had similar communities, despite ranging from 14°N to 20°S latitude, and they differed from the Mediterranean station. Single cell genomics of protists suggests that the taxonomic diversity of the dominant taxa found in only several hundreds of microliters of surface seawater is similar to that found in molecular surveys where liters of sample are filtered.