Survey of the green picoalga Bathycoccus genomes in the global ocean

An INDEL genomic approach to explore population diversity of phytoplankton

BACKGROUND

Although metabarcoding and metagenomic approaches have generated large datasets on worldwide phytoplankton species diversity, the intraspecific genetic diversity underlying the genetic adaptation of marine phytoplankton to specific environmental niches remains largely unexplored. This is mainly due to the lack of biological resources and tools for monitoring the dynamics of this diversity in space and time.

RESULTS

To gain insight into population diversity, a novel method based on INDEL markers was developed on Bathycoccus prasinos (Mamiellophyceae), an abundant and cosmopolitan species with strong seasonal patterns. Long read sequencing was first used to characterize structural variants among the genomes of six B. prasinos strains sampled from geographically distinct regions in the world ocean. Markers derived from identified insertions/deletions were validated by PCR then used to genotype 55 B. prasinos strains isolated during the winter bloom 2018-2019 in the bay of Banyuls-sur-Mer (Mediterranean Sea, France). This led to their classification into eight multi-loci genotypes and the sequencing of strains representative of local diversity, further improving the available genetic diversity of B. prasinos. Finally, selected markers were directly tracked on environmental DNA sampled during 3 successive blooms from 2018 to 2021, showcasing a fast and cost-effective approach to follow local population dynamics.

CONCLUSIONS

This method, which involves (i) pre-identifying the genetic diversity of B. prasinos in environmental samples by PCR, (ii) isolating cells from selected environmental samples and (iii) identifying genotypes representative of B. prasinos diversity for sequencing, can be used to comprehensively describe the diversity and population dynamics not only in B. prasinos but also potentially in other generalist phytoplankton species.

Diverse patterns of correspondence between protist metabarcodes and protist metagenome-assembled genomes

Two common approaches to study the composition of environmental protist communities are metabarcoding and metagenomics. Raw metabarcoding data are usually processed into Operational Taxonomic Units (OTUs) or amplicon sequence variants (ASVs) through clustering or denoising approaches, respectively. Analogous approaches are used to assemble metagenomic reads into metagenome-assembled genomes (MAGs). Understanding the correspondence between the data produced by these two approaches can help to integrate information between the datasets and to explain how metabarcoding OTUs and MAGs are related with the underlying biological entities they are hypothesised to represent. MAGs do not contain the commonly used barcoding loci, therefore sequence homology approaches cannot be used to match OTUs and MAGs. We made an attempt to match V9 metabarcoding OTUs from the 18S rRNA gene (V9 OTUs) and MAGs from the Tara Oceans expedition based on the correspondence of their relative abundances across the same set of samples. We evaluated several metrics for detecting correspondence between features in these two datasets and developed controls to filter artefacts of data structure and processing. After selecting the best-performing metrics, ranking the V9 OTU/MAG matches by their proportionality/correlation coefficients and applying a set of selection criteria, we identified candidate matches between V9 OTUs and MAGs. In some cases, V9 OTUs and MAGs could be matched with a one-to-one correspondence, implying that they likely represent the same underlying biological entity. More generally, matches we observed could be classified into 4 scenarios: one V9 OTU matches many MAGs; many V9 OTUs match many MAGs; many V9 OTUs match one MAG; one V9 OTU matches one MAG. Notably, we found some instances in which different OTU-MAG matches from the same taxonomic group were not classified in the same scenario, with all four scenarios possible even within the same taxonomic group, illustrating that factors beyond taxonomic lineage influence the relationship between OTUs and MAGs. Overall, each scenario produces a different interpretation of V9 OTUs, MAGs and how they compare in terms of the genomic and ecological diversity they represent.

Decoding populations in the ocean microbiome

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

Hidden genomic diversity drives niche partitioning in a cosmopolitan eukaryotic picophytoplankton

Marine eukaryotic phytoplankton are fundamental to the marine food web, yet the lack of reference genomes or just a single genome representing a taxon has led to an underestimation of their taxonomic, adaptive, and functional diversity. Here, we integrated strain isolation with metagenomic binning to recover genomes from the cosmopolitan picophytoplankton genus Bathycoccus, traditionally considered monospecific. Our recovery and analysis of 37 Bathycoccus genomes delineated their global genomic diversity and established four evolutionary clades (BI, BII, BIII, BIV). Our metagenomic abundance survey revealed well-differentiated ecological niches and distinct biogeographic distributions for each clade, predominantly shaped by temperature, salinity, and nutrient availability. Comparative genomics analyses further revealed clade-specific genomic traits that underpin niche adaptation and contribute to the global prevalence of Bathycoccus. Our findings underscore temperature as a major driver of genome diversification in this genus, with clade divergences coinciding with major paleoclimatic events that influenced their contemporary thermal niches. Moreover, the unique enrichment of C2H2 zinc finger and ankyrin repeat gene families in polar-adapted clades suggests previously unrecognized cold-adaptation mechanisms in marine eukaryotic phytoplankton. Our study offers a comprehensive genomic landscape of this crucial eukaryotic picophytoplankton, providing insights into their microdiversity and adaptive evolution in response to changing environments.

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.

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.

Functional repertoire convergence of distantly related eukaryotic plankton lineages abundant in the sunlit ocean

Marine planktonic eukaryotes play critical roles in global biogeochemical cycles and climate. However, their poor representation in culture collections limits our understanding of the evolutionary history and genomic underpinnings of planktonic ecosystems. Here, we used 280 billion Tara Oceans metagenomic reads from polar, temperate, and tropical sunlit oceans to reconstruct and manually curate more than 700 abundant and widespread eukaryotic environmental genomes ranging from 10 Mbp to 1.3 Gbp. This genomic resource covers a wide range of poorly characterized eukaryotic lineages that complement long-standing contributions from culture collections while better representing plankton in the upper layer of the oceans. We performed the first, to our knowledge, comprehensive genome-wide functional classification of abundant unicellular eukaryotic plankton, revealing four major groups connecting distantly related lineages. Neither trophic modes of plankton nor its vertical evolutionary history could completely explain the functional repertoire convergence of major eukaryotic lineages that coexisted within oceanic currents for millions of years.

Can Aggregate-Associated Organisms Influence the Fouling in a SWRO Desalination Plant?

This pilot study investigates the formation of aggregates within a desalination plant, before and after pre-treatment, as well as their potential impact on fouling. The objective is to provide an understanding of the biofouling potential of the feed water within a seawater reverse osmosis (SWRO) desalination plant, due to the limited removal of fouling precursors. The 16S and 18S rRNA was extracted from the water samples, and the aggregates and sequenced. Pre-treatment systems, within the plant remove < 5 µm precursors and organisms; however, smaller size particles progress through the plant, allowing for the formation of aggregates. These become hot spots for microbes, due to their nutrient gradients, facilitating the formation of niche environments, supporting the proliferation of those organisms. Aggregate-associated organisms are consistent with those identified on fouled SWRO membranes. This study examines, for the first time, the factors supporting the formation of aggregates within a desalination system, as well as their microbial communities and biofouling potential.

Diversity and Evolution of Mamiellophyceae: Early-Diverging Phytoplanktonic Green Algae Containing Many Cosmopolitan Species

The genomic revolution has bridged a gap in our knowledge about the diversity, biology and evolution of unicellular photosynthetic eukaryotes, which bear very few discriminating morphological features among species from the same genus. The high-quality genome resources available in the class Mamiellophyceae (Chlorophyta) have been paramount to estimate species diversity and screen available metagenomic data to assess the biogeography and ecological niches of different species on a global scale. Here we review the current knowledge about the diversity, ecology and evolution of the Mamiellophyceae and the large double-stranded DNA prasinoviruses infecting them, brought by the combination of genomic and metagenomic analyses, including 26 metabarcoding environmental studies, as well as the pan-oceanic GOS and the Tara Oceans expeditions.

Size-Fractionated Microbiome Structure in Subarctic Rivers and a Coastal Plume Across DOC and Salinity Gradients

Little is known about the microbial diversity of rivers that flow across the changing subarctic landscape. Using amplicon sequencing (rRNA and rRNA genes) combined with HPLC pigment analysis and physicochemical measurements, we investigated the diversity of two size fractions of planktonic Bacteria, Archaea and microbial eukaryotes along environmental gradients in the Great Whale River (GWR), Canada. This large subarctic river drains an extensive watershed that includes areas of thawing permafrost, and discharges into southeastern Hudson Bay as an extensive plume that gradually mixes with the coastal marine waters. The microbial communities differed by size-fraction (separated with a 3-μm filter), and clustered into three distinct environmental groups: (1) the GWR sites throughout a 150-km sampling transect; (2) the GWR plume in Hudson Bay; and (3) small rivers that flow through degraded permafrost landscapes. There was a downstream increase in taxonomic richness along the GWR, suggesting that sub-catchment inputs influence microbial community structure in the absence of sharp environmental gradients. Microbial community structure shifted across the salinity gradient within the plume, with changes in taxonomic composition and diversity. Rivers flowing through degraded permafrost had distinct physicochemical and microbiome characteristics, with allochthonous dissolved organic carbon explaining part of the variation in community structure. Finally, our analyses of the core microbiome indicated that while a substantial part of all communities consisted of generalists, most taxa had a more limited environmental range and may therefore be sensitive to ongoing change.

Viruses infecting a warm water picoeukaryote shed light on spatial co-occurrence dynamics of marine viruses and their hosts

The marine picoeukaryote Bathycoccus prasinos has been considered a cosmopolitan alga, although recent studies indicate two ecotypes exist, Clade BI (B. prasinos) and Clade BII. Viruses that infect Bathycoccus Clade BI are known (BpVs), but not that infect BII. We isolated three dsDNA prasinoviruses from the Sargasso Sea against Clade BII isolate RCC716. The BII-Vs do not infect BI, and two (BII-V2 and BII-V3) have larger genomes (~210 kb) than BI-Viruses and BII-V1. BII-Vs share ~90% of their proteins, and between 65% to 83% of their proteins with sequenced BpVs. Phylogenomic reconstructions and PolB analyses establish close-relatedness of BII-V2 and BII-V3, yet BII-V2 has 10-fold higher infectivity and induces greater mortality on host isolate RCC716. BII-V1 is more distant, has a shorter latent period, and infects both available BII isolates, RCC716 and RCC715, while BII-V2 and BII-V3 do not exhibit productive infection of the latter in our experiments. Global metagenome analyses show Clade BI and BII algal relative abundances correlate positively with their respective viruses. The distributions delineate BI/BpVs as occupying lower temperature mesotrophic and coastal systems, whereas BII/BII-Vs occupy warmer temperature, higher salinity ecosystems. Accordingly, with molecular diagnostic support, we name Clade BII Bathycoccus calidus sp. nov. and propose that molecular diversity within this new species likely connects to the differentiated host-virus dynamics observed in our time course experiments. Overall, the tightly linked biogeography of Bathycoccus host and virus clades observed herein supports species-level host specificity, with strain-level variations in infection parameters.

The pico-sized Mamiellophyceae and a novel Bathycoccus clade from the summer plankton of Russian Arctic Seas and adjacent waters

Global climate changes and anthropogenic activity greatly impact Arctic marine biodiversity including phytoplankton which contribute greatly to atmospheric oxygen production. Thus the study of microalgae has rising topicality. Class Mamiellophyceae is an important component of phototrophic picoplankton. To gain more knowledge about Mamiellophyceae distribution and diversity special studies were performed in such remote areas as the Russian Arctic seas. A metabarcoding of pico-sized Mamiellophyceae was undertaken by high-throughput sequencing of the 18S rRNA gene sequence V4 region from samples collected in July-September 2017 in the Barents, Kara and Laptev seas, and in the adjacent waters of the Norwegian Sea. Our study is the first to show that Mamiellophyceae among the summer picoplankton of Russian Arctic seas are diverse and represented by 16 algae species/phylotypes. We discovered a new candidate species of Bathycoccus assigned to a new Bathycoccus clade A-uncultured Bathycoccus Kara 2017. It was found that several Micromonas species can co-exist, with Micromonas polaris dominating north of 72°N. The presence of Ostreococcus tauri, Ostreococcus lucimarinus and Ostreococcus mediterraneus at high latitudes beyond 65°N was documented for the first time, similar to findings for some other taxa. Our results will be important for obtaining a global view of Mamiellophyceae community dynamics.

No evidence of Phago-mixotropy in Micromonas polaris (Mamiellophyceae), the Dominant Picophytoplankton Species in the Arctic

In the Arctic Ocean, the small green alga Micromonas polaris dominates picophytoplankton during the summer months but is also present in winter. It has been previously hypothesized to be phago-mixotrophic (capable of bacteria ingestion) based on laboratory and field experiments. Prey uptake was analyzed in several M. polaris strains isolated from different regions and depths of the Arctic Ocean and in Ochromonas triangulata, a known phago-mixotroph used as a control. Measuring ingestion of either fluorescent beads or fluorescently labeled bacteria by flow cytometry, we found no evidence of phago-mixotrophy in any M. polaris strain while O. triangulata was ingesting both beads and bacteria. In addition, in silico predictions revealed that members of the genus Micromonas lack a genetic signature of phagocytotic capacity.

Annual phytoplankton dynamics in coastal waters from Fildes Bay, Western Antarctic Peninsula

Year-round reports of phytoplankton dynamics in the West Antarctic Peninsula are rare and mainly limited to microscopy and/or pigment-based studies. We analyzed the phytoplankton community from coastal waters of Fildes Bay in the West Antarctic Peninsula between January 2014 and 2015 using metabarcoding of the nuclear and plastidial 18/16S rRNA gene from both size-fractionated and flow cytometry sorted samples. Overall 14 classes of photosynthetic eukaryotes were present in our samples with the following dominating: Bacillariophyta (diatoms), Pelagophyceae and Dictyochophyceae for division Ochrophyta, Mamiellophyceae and Pyramimonadophyceae for division Chlorophyta, Haptophyta and Cryptophyta. Each metabarcoding approach yielded a different image of the phytoplankton community with for example Prymnesiophyceae more prevalent in plastidial metabarcodes and Mamiellophyceae in nuclear ones. Diatoms were dominant in the larger size fractions and during summer, while Prymnesiophyceae and Cryptophyceae were dominant in colder seasons. Pelagophyceae were particularly abundant towards the end of autumn (May). In addition of Micromonas polaris and Micromonas sp. clade B3, both previously reported in Arctic waters, we detected a new Micromonas 18S rRNA sequence signature, close to, but clearly distinct from M. polaris, which potentially represents a new clade specific of the Antarctic. These results highlight the need for complementary strategies as well as the importance of year-round monitoring for a comprehensive description of phytoplankton communities in Antarctic coastal waters.

Estimating the quality of eukaryotic genomes recovered from metagenomic analysis with EukCC

Microbial eukaryotes constitute a significant fraction of biodiversity and have recently gained more attention, but the recovery of high-quality metagenomic assembled eukaryotic genomes is limited by the current availability of tools. To help address this, we have developed EukCC, a tool for estimating the quality of eukaryotic genomes based on the automated dynamic selection of single copy marker gene sets. We demonstrate that our method outperforms current genome quality estimators, particularly for estimating contamination, and have applied EukCC to datasets derived from two different environments to enable the identification of novel eukaryote genomes, including one from the human skin.

Diversity and evolution of bacterial bioluminescence genes in the global ocean

Although bioluminescent bacteria are the most abundant and widely distributed of all light-emitting organisms, the biological role and evolutionary history of bacterial luminescence are still shrouded in mystery. Bioluminescence has so far been observed in the genomes of three families of Gammaproteobacteria in the form of canonical lux operons that adopt the CDAB(F)E(G) gene order. LuxA and luxB encode the two subunits of bacterial luciferase responsible for light-emission. Our deep exploration of public marine environmental databases considerably expands this view by providing a catalog of new lux homolog sequences, including 401 previously unknown luciferase-related genes. It also reveals a broader diversity of the lux operon organization, which we observed in previously undescribed configurations such as CEDA, CAED and AxxCE. This expanded operon diversity provides clues for deciphering lux operon evolution and propagation within the bacterial domain. Leveraging quantitative tracking of marine bacterial genes afforded by planetary scale metagenomic sampling, our study also reveals that the novel lux genes and operons described herein are more abundant in the global ocean than the canonical CDAB(F)E(G) operon.

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.

Transcriptome reconstruction and functional analysis of eukaryotic marine plankton communities via high-throughput metagenomics and metatranscriptomics

Large-scale metagenomic and metatranscriptomic data analyses are often restricted by their gene-centric approach, limiting the ability to understand organismal and community biology. De novo assembly of large and mosaic eukaryotic genomes from complex meta-omics data remains a challenging task, especially in comparison with more straightforward bacterial and archaeal systems. Here, we use a transcriptome reconstruction method based on clustering co-abundant genes across a series of metagenomic samples. We investigated the co-abundance patterns of ∼37 million eukaryotic unigenes across 365 metagenomic samples collected during the Tara Oceans expeditions to assess the diversity and functional profiles of marine plankton. We identified ∼12,000 co-abundant gene groups (CAGs), encompassing ∼7 million unigenes, including 924 metagenomics-based transcriptomes (MGTs, CAGs larger than 500 unigenes). We demonstrated the biological validity of the MGT collection by comparing individual MGTs with available references. We identified several key eukaryotic organisms involved in dimethylsulfoniopropionate (DMSP) biosynthesis and catabolism in different oceanic provinces, thus demonstrating the potential of the MGT collection to provide functional insights on eukaryotic plankton. We established the ability of the MGT approach to capture interspecies associations through the analysis of a nitrogen-fixing haptophyte-cyanobacterial symbiotic association. This MGT collection provides a valuable resource for analyses of eukaryotic plankton in the open ocean by giving access to the genomic content and functional potential of many ecologically relevant eukaryotic species.

Photosynthetic Picoeukaryotes Diversity in the Underlying Ice Waters of the White Sea, Russia

The White Sea is a unique basin combining features of temperate and arctic seas. The current state of its biocenoses can serve as a reference point in assessing the expected desalination of the ocean as a result of climate change. A metagenomic study of under-ice ice photosynthetic picoeukaryotes (PPEs) was undertaken by Illumina high-throughput sequencing of the 18S rDNA V4 region from probes collected in March 2013 and 2014. The PPE biomass in samples was 0.03–0.17 µg C·L−1 and their abundance varied from 10 cells·mL−1 to 140 cells·mL−1. There were representatives of 16 algae genera from seven classes and three supergroups, but Chlorophyta, especially Mamiellophyceae, dominated. The most represented genera were Micromonas and Mantoniella. For the first time, the predominance of Mantoniella (in four samples) and Bolidophyceae (in one sample) was observed in under-ice water. It can be assumed that a change in environmental conditions will lead to a considerable change in the structure of arctic PPE communities.

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.

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'.

A streamlined and predominantly diploid genome in the tiny marine green alga Chloropicon primus

Tiny marine green algae issued from two deep branches of the Chlorophyta, the Mamiellophyceae and Chloropicophyceae, dominate different regions of the oceans and play key roles in planktonic communities. Considering that the Mamiellophyceae is the sole lineage of prasinophyte algae that has been intensively investigated, the extent to which these two algal groups differ in their metabolic capacities and cellular processes is currently unknown. To address this gap of knowledge, we investigate here the nuclear genome sequence of a member of the Chloropicophyceae, Chloropicon primus. Among the main biological insights that emerge from this 17.4 Mb genome, we find an unexpected diploid structure for most chromosomes and a propionate detoxification pathway in green algae. Our results support the notion that separate events of genome minimization, which entailed differential losses of genes/pathways, have occurred in the Chloropicophyceae and Mamiellophyceae, suggesting different strategies of adaptation to oceanic environments.

Simplified Transformation of Ostreococcus tauri Using Polyethylene Glycol

Ostreococcustauri is an easily cultured representative of unicellular algae (class Mamiellophyceae) that abound in oceans worldwide. Eight complete 13-22 Mb genomes of phylogenetically divergent species within this class are available, and their DNA sequences are nearly always present in metagenomic data produced from marine samples. Here we describe a simplified and robust transformation protocol for the smallest of these algae (O. tauri). Polyethylene glycol (PEG) treatment was much more efficient than the previously described electroporation protocol. Short (2 min or less) incubation times in PEG gave >104 transformants per microgram DNA. The time of cell recovery after transformation could be reduced to a few hours, permitting the experiment to be done in a day rather than overnight as used in previous protocols. DNA was randomly inserted in the O. tauri genome. In our hands PEG was 20-40-fold more efficient than electroporation for the transformation of O. tauri, and this improvement will facilitate mutagenesis of all of the dispensable genes present in the tiny O. tauri genome.

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.

Novel diversity within marine Mamiellophyceae (Chlorophyta) unveiled by metabarcoding

Mamiellophyceae (unicellular green algae) are a key phytoplankton group in coastal waters. Although extensively studied over the last 20 years, the overall oceanic distribution of the major species/clades is still poorly known. To address this problem, we analyzed the 2014 Ocean Sampling Day (OSD) metabarcoding dataset providing sequences from the V4 hypervariable region of the 18S rRNA gene for 157 samples collected at 143 mostly coastal stations. Mamiellophyceae were found at nearly all OSD stations and represented 55% of the green microalgae (Chlorophyta) reads. We performed phylogenetic analyses of unique OSD metabarcodes (amplicon single variants, ASVs) and GenBank reference sequences from cultures and from the environment, focusing on the four most represented genera: Ostreococcus (45% of the Mamiellophyceae reads), Micromonas (34%), Bathycoccus (10%) and Mantoniella (8.7%). These analyses uncovered novel diversity within each genus except Bathycoccus. In Ostreococcus, a new clade (E) was the second most represented clade after Ostreococcus "lucimarinus". Micromonas could be separated into nine clades, exceeding the six species and candidate species already described. Finally, we found two new environmental clades within Mantoniella. Each Mamiellophyceae clade had a specific distribution in the OSD dataset suggesting that they are adapted to different ecological niches.

Live cell analysis at sea reveals divergent thermal performance between photosynthetic ocean microbial eukaryote populations

Experimentation at sea provides insight into which traits of ocean microbes are linked to performance in situ. Here we show distinct patterns in thermal tolerance of microbial phototrophs from adjacent water masses sampled in the south-west Pacific Ocean, determined using a fluorescent marker for reactive oxygen species (ROS). ROS content of pico-eukaryotes was assessed after 1, 5 and 25 h of incubation along a temperature gradient (15.6-32.1 °C). Pico-eukaryotes from the East Australian Current (EAC) had relatively constant ROS and showed greatest mortality after 25 h at 7 °C below ambient, whereas those from the Tasman Sea had elevated ROS in both warm and cool temperature extremes and greatest mortality at temperatures 6-10 °C above ambient, interpreted as the outcome of thermal stress. Tracking of water masses within an oceanographic circulation model showed populations had distinct thermal histories, with EAC pico-eukaryotes experiencing higher average temperatures for at least 1 week prior to sampling. While acclimatization and community assembly could both influence biological responses, this study clearly demonstrates that phenotypic divergence occurs along planktonic drift trajectories.

Rapidity of Genomic Adaptations to Prasinovirus Infection in a Marine Microalga

Prasinoviruses are large dsDNA viruses commonly found in aquatic systems worldwide, where they can infect and lyse unicellular prasinophyte algae such as Ostreococcus. Host susceptibility is virus strain-specific, but resistance of susceptible Ostreococcus tauri strains to a virulent virus arises frequently. In clonal resistant lines that re-grow, viruses are usually present for many generations, and genes clustered on chromosome 19 show physical rearrangements and differential expression. Here, we investigated changes occurring during the first two weeks after inoculation of the prasinovirus OtV5. By serial dilutions of cultures at the time of inoculation, we estimated the frequency of resistant cells arising in virus-challenged O. tauri cultures to be 10^-3⁻10^-4 of the inoculated population. Re-growing resistant cells were detectable by flow cytometry 3 days post-inoculation (dpi), visible re-greening of cultures occurred by 6 dpi, and karyotypic changes were visually detectable at 8 dpi. Resistant cell lines showed a modified spectrum of host-virus specificities and much lower levels of OtV5 adsorption.

Small eukaryotic phytoplankton communities in tropical waters off Brazil are dominated by symbioses between Haptophyta and nitrogen-fixing cyanobacteria

Symbioses between eukaryotic algae and nitrogen-fixing cyanobacteria have been recognized in recent years as a key source of new nitrogen in the oceans. We investigated the composition of the small photosynthetic eukaryote communities associated with nitrogen-fixing cyanobacteria in the Brazilian South Atlantic Bight using a combination of flow cytometry sorting and high throughput sequencing of two genes: the V4 region of 18S rRNA and nifH. Two distinct eukaryotic communities were often encountered, one dominated by the Mamiellophyceae Bathycoccus and Ostreococcus, and one dominated by a prymnesiophyte known to live in symbiosis with the UCYN-A1 nitrogen-fixing cyanobacterium. Among nifH sequences, those from UCYN-A1 were most abundant but three other UCYN-A clades (A2, A3, A4) were also found. Network analysis confirmed the relation between A1 and A2 clades and their hypothesized hosts and pointed out to the potential association between novel clade A4 with Braarudosphaera bigelowii, previously hypothesized to host A2.

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

Single-celled eukaryotes (protists) are critical players in global biogeochemical cycling of nutrients and energy in the oceans. While their roles as primary producers and grazers are well appreciated, other aspects of their life histories remain obscure due to challenges in culturing and sequencing their natural diversity. Here, we exploit single-cell genomics and metagenomics data from the circumglobal Tara Oceans expedition to analyze the genome content and apparent oceanic distribution of seven prevalent lineages of uncultured heterotrophic stramenopiles. Based on the available data, each sequenced genome or genotype appears to have a specific oceanic distribution, principally correlated with water temperature and depth. The genome content provides hypotheses for specialization in terms of cell motility, food spectra, and trophic stages, including the potential impact on their lifestyles of horizontal gene transfer from prokaryotes. Our results support the idea that prominent heterotrophic marine protists perform diverse functions in ocean ecology.

The biology of many marine protists, such as stramenopiles, remains obscure. Here, the authors exploit single-cell genomics and metagenomics to analyze the genome content and apparent oceanic distribution of seven prevalent lineages of uncultured heterotrophic stramenopiles.

Tools for Genomic and Transcriptomic Analysis of Microbes at Single-Cell Level

Microbiologists traditionally study population rather than individual cells, as it is generally assumed that the status of individual cells will be similar to that observed in the population. However, the recent studies have shown that the individual behavior of each single cell could be quite different from that of the whole population, suggesting the importance of extending traditional microbiology studies to single-cell level. With recent technological advances, such as flow cytometry, next-generation sequencing (NGS), and microspectroscopy, single-cell microbiology has greatly enhanced the understanding of individuality and heterogeneity of microbes in many biological systems. Notably, the application of multiple ‘omics’ in single-cell analysis has shed light on how individual cells perceive, respond, and adapt to the environment, how heterogeneity arises under external stress and finally determines the fate of the whole population, and how microbes survive under natural conditions. As single-cell analysis involves no axenic cultivation of target microorganism, it has also been demonstrated as a valuable tool for dissecting the microbial ‘dark matter.’ In this review, current state-of-the-art tools and methods for genomic and transcriptomic analysis of microbes at single-cell level were critically summarized, including single-cell isolation methods and experimental strategies of single-cell analysis with NGS. In addition, perspectives on the future trends of technology development in the field of single-cell analysis was also presented.

Revision of the Genus Micromonas Manton et Parke (Chlorophyta, Mamiellophyceae), of the Type Species M. pusilla (Butcher) Manton & Parke and of the Species M. commoda van Baren, Bachy and Worden and Description of Two New Species Based on the Genetic and Phenotypic Characterization of Cultured Isolates

The green picoalgal genus Micromonas is broadly distributed in estuaries, coastal marine habitats and open oceans, from the equator to the poles. Phylogenetic, ecological and genomic analyses of culture strains and natural populations have suggested that this cosmopolitan genus is composed of several cryptic species corresponding to genetic lineages. We performed a detailed analysis of variations in morphology, pigment content, and sequences of the nuclear-encoded small-subunit rRNA gene and the second internal transcribed spacer (ITS2) from strains isolated worldwide. A new morphological feature of the genus, the presence of tip hairs at the extremity of the hair point, was discovered and subtle differences in hair point length were detected between clades. Clear non-homoplasious synapomorphies were identified in the small-subunit rRNA gene and ITS2 spacer sequences of five genetic lineages. These findings lead us to provide emended descriptions of the genus Micromonas, of the type species M. pusilla, and of the recently described species M. commoda, as well as to describe 2 new species, M. bravo and M. polaris. By clarifying the status of the genetic lineages identified within Micromonas, these formal descriptions will facilitate further interpretations of large-scale analyses investigating ecological trends in time and space for this widespread picoplankter.

Quantitative biogeography of picoprasinophytes establishes ecotype distributions and significant contributions to marine phytoplankton

Bathycoccus and Ostreococcus are broadly distributed marine picoprasinophyte algae. We enumerated small phytoplankton using flow cytometry and qPCR assays for phylogenetically distinct Bathycoccus clades BI and BII and Ostreococcus clades OI and OII. Among 259 photic-zone samples from transects and time-series, Ostreococcus maxima occurred in the North Pacific coastal upwelling for OI (36 713 ± 1485 copies ml^-1 ) and the Kuroshio Front for OII (50 189 ± 561 copies ml^-1 ) and the two overlapped only in frontal regions. The Bathycoccus overlapped more often with maxima along Line-P for BI (10 667 ± 1299 copies ml^-1 ) and the tropical Atlantic for BII (4125 ± 339 copies ml^-1 ). Only BII and OII were detected at warm oligotrophic sites, accounting for 34 ± 13% of 1589 ± 448 eukaryotic phytoplankton cells ml^-1 (annual average) at Station ALOHA's deep chlorophyll maximum. Significant distributional and molecular differences lead us to propose that Bathycoccus clade BII represents a separate species which tolerates higher temperature oceanic conditions than Bathycoccus prasinos (BI). Morphological differences were not evident, but quick-freeze deep-etch electron microscopy provided insight into Bathycoccus scale formation. Our results highlight the importance of quantitative seasonal abundance data for inferring ecological distributions and demonstrate significant, differential picoprasinophyte contributions in mesotrophic and open-ocean waters.

Spontaneous Mutation Rate in the Smallest Photosynthetic Eukaryotes

Mutation is the ultimate source of genetic variation, and knowledge of mutation rates is fundamental for our understanding of all evolutionary processes. High throughput sequencing of mutation accumulation lines has provided genome wide spontaneous mutation rates in a dozen model species, but estimates from nonmodel organisms from much of the diversity of life are very limited. Here, we report mutation rates in four haploid marine bacterial-sized photosynthetic eukaryotic algae; Bathycoccus prasinos, Ostreococcus tauri, Ostreococcus mediterraneus, and Micromonas pusilla. The spontaneous mutation rate between species varies from μ = 4.4 × 10-10 to 9.8 × 10-10 mutations per nucleotide per generation. Within genomes, there is a two-fold increase of the mutation rate in intergenic regions, consistent with an optimization of mismatch and transcription-coupled DNA repair in coding sequences. Additionally, we show that deviation from the equilibrium GC content increases the mutation rate by ∼2% to ∼12% because of a GC bias in coding sequences. More generally, the difference between the observed and equilibrium GC content of genomes explains some of the inter-specific variation in mutation rates.