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

Recurring seasonality exposes dominant species and niche partitioning strategies of open ocean picoeukaryotic algae

Ocean spring phytoplankton blooms are dynamic periods important to global primary production. We document vertical patterns of a diverse suite of eukaryotic algae, the prasinophytes, in the North Atlantic Subtropical Gyre with monthly sampling over four years at the Bermuda Atlantic Time-series Study site. Water column structure was used to delineate seasonal stability periods more ecologically relevant than seasons defined by calendar dates. During winter mixing, tiny prasinophytes dominated by Class II comprise 46  ±  24% of eukaryotic algal (plastid-derived) 16S rRNA V1-V2 amplicons, specifically Ostreococcus Clade OII, Micromonas commoda, and Bathycoccus calidus. In contrast, Class VII are rare and Classes I and VI peak during warm stratified periods when surface eukaryotic phytoplankton abundances are low. Seasonality underpins a reservoir of genetic diversity from multiple prasinophyte classes during warm periods that harbor ephemeral taxa. Persistent Class II sub-species dominating the winter/spring bloom period retreat to the deep chlorophyll maximum in summer, poised to seed the mixed layer upon winter convection, exposing a mechanism for initiating high abundances at bloom onset. Comparisons to tropical oceans reveal broad distributions of the dominant sub-species herein. This unparalleled window into temporal and spatial niche partitioning of picoeukaryotic primary producers demonstrates how key prasinophytes prevail in warm oceans.

Gradients of bacteria in the oceanic water column reveal finely-resolved vertical distributions

Bacterial communities directly influence ecological processes in the ocean, and depth has a major influence due to the changeover in primary energy sources between the sunlit photic zone and dark ocean. Here, we examine the abundance and diversity of bacteria in Monterey Bay depth profiles collected from the surface to just above the sediments (e.g., 2000 m). Bacterial abundance in these Pacific Ocean samples decreased by >1 order of magnitude, from 1.22 ±0.69 ×106 cells ml-1 in the variable photic zone to 1.44 ± 0.25 ×105 and 6.71 ± 1.23 ×104 cells ml-1 in the mesopelagic and bathypelagic, respectively. V1-V2 16S rRNA gene profiling showed diversity increased sharply between the photic and mesopelagic zones. Weighted Gene Correlation Network Analysis clustered co-occurring bacterial amplicon sequence variants (ASVs) into seven subnetwork modules, of which five strongly correlated with depth-related factors. Within surface-associated modules there was a clear distinction between a 'copiotrophic' module, correlating with chlorophyll and dominated by e.g., Flavobacteriales and Rhodobacteraceae, and an 'oligotrophic' module dominated by diverse Oceanospirillales (such as uncultured JL-ETNP-Y6, SAR86) and Pelagibacterales. Phylogenetic reconstructions of Pelagibacterales and SAR324 using full-length 16S rRNA gene data revealed several additional subclades, expanding known microdiversity within these abundant lineages, including new Pelagibacterales subclades Ia.B, Id, and IIc, which comprised 4-10% of amplicons depending on the subclade and depth zone. SAR324 and Oceanospirillales dominated in the mesopelagic, with SAR324 clade II exhibiting its highest relative abundances (17±4%) in the lower mesopelagic (300-750 m). The two newly-identified SAR324 clades showed highest relative abundances in the photic zone (clade III), while clade IV was extremely low in relative abundance, but present across dark ocean depths. Hierarchical clustering placed microbial communities from 900 m samples with those from the bathypelagic, where Marinimicrobia was distinctively relatively abundant. The patterns resolved herein, through high resolution and statistical replication, establish baselines for marine bacterial abundance and taxonomic distributions across the Monterey Bay water column, against which future change can be assessed.

Genetic and physiological responses to light quality in a deep ocean ecotype of Ostreococcus, an ecologically important photosynthetic picoeukaryote

Phytoplankton are exposed to dramatic variations in light quality when cells are carried by upwelling or downwelling currents or encounter sediment. We investigated the potential impact of light quality changes in Ostreococcus, a key marine photosynthetic picoeukaryote, by analysing changes in its transcriptome, pigment content, and photophysiology after acclimation to monochromatic red, green, or blue light. The clade B species RCC809, isolated from the deep euphotic zone of the tropical Atlantic Ocean, responded to blue light by accelerating cell division at the expense of storage reserves and by increasing the relative level of blue-light-absorbing pigments. It responded to red and green light by increasing its potential for photoprotection. In contrast, the clade A species OTTH0595, which originated from a shallow water environment, showed no difference in photosynthetic properties and minor differences in carotenoid contents between light qualities. This was associated with the loss of candidate light-quality responsive promoter motifs identified in RCC809 genes. These results demonstrate that light quality can have a major influence on the physiology of eukaryotic phytoplankton and suggest that different light quality environments can drive selection for diverse patterns of responsiveness and environmental niche partitioning.

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

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

The Land-Sea Connection: Insights Into the Plant Lineage from a Green Algal Perspective

The colonization of land by plants generated opportunities for the rise of new heterotrophic life forms, including humankind. A unique event underpinned this massive change to earth ecosystems-the advent of eukaryotic green algae. Today, an abundant marine green algal group, the prasinophytes, alongside prasinodermophytes and nonmarine chlorophyte algae, is facilitating insights into plant developments. Genome-level data allow identification of conserved proteins and protein families with extensive modifications, losses, or gains and expansion patterns that connect to niche specialization and diversification. Here, we contextualize attributes according to Viridiplantae evolutionary relationships, starting with orthologous protein families, and then focusing on key elements with marked differentiation, resulting in patchy distributions across green algae and plants. We place attention on peptidoglycan biosynthesis, important for plastid division and walls; phytochrome photosensors that are master regulators in plants; and carbohydrate-active enzymes, essential to all manner of carbohydratebiotransformations. Together with advances in algal model systems, these areas are ripe for discovering molecular roles and innovations within and across plant and algal lineages.

Phytoplankton Surveys in the Arctic Fram Strait Demonstrate the Tiny Eukaryotic Alga Micromonas and Other Picoprasinophytes Contribute to Deep Sea Export

Critical questions exist regarding the abundance and, especially, the export of picophytoplankton (≤2 µm diameter) in the Arctic. These organisms can dominate chlorophyll concentrations in Arctic regions, which are subject to rapid change. The picoeukaryotic prasinophyte Micromonas grows in polar environments and appears to constitute a large, but variable, proportion of the phytoplankton in these waters. Here, we analyze 81 samples from the upper 100 m of the water column from the Fram Strait collected over multiple years (2009−2015). We also analyze sediment trap samples to examine picophytoplankton contributions to export, using both 18S rRNA gene qPCR and V1-V2 16S rRNA Illumina amplicon sequencing to assess the Micromonas abundance within the broader diversity of photosynthetic eukaryotes based on the phylogenetic placement of plastid-derived 16S amplicons. The material sequenced from the sediment traps in July and September 2010 showed that 11.2 ± 12.4% of plastid-derived amplicons are from picoplanktonic prasinophyte algae and other green lineage (Viridiplantae) members. In the traps, Micromonas dominated (83.6 ± 21.3%) in terms of the overall relative abundance of Viridiplantae amplicons, specifically the species Micromonas polaris. Temporal variations in Micromonas abundances quantified by qPCR were also observed, with higher abundances in the late-July traps and deeper traps. In the photic zone samples, four prasinophyte classes were detected in the amplicon data, with Micromonas again being the dominant prasinophyte, based on the relative abundance (89.4 ± 8.0%), but with two species (M. polaris and M. commoda-like) present. The quantitative PCR assessments showed that the photic zone samples with higher Micromonas abundances (>1000 gene copies per mL) had significantly lower standing stocks of phosphate and nitrate, and a shallower average depth (20 m) than those with fewer Micromonas. This study shows that despite their size, prasinophyte picophytoplankton are exported to the deep sea, and that Micromonas is particularly important within this size fraction in Arctic marine ecosystems.

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.

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

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

Small pigmented eukaryote assemblages of the western tropical North Atlantic around the Amazon River plume during spring discharge

Small pigmented eukaryotes (⩽ 5 µm) are an important, but overlooked component of global marine phytoplankton. The Amazon River plume delivers nutrients into the oligotrophic western tropical North Atlantic, shades the deeper waters, and drives the structure of microphytoplankton (> 20 µm) communities. For small pigmented eukaryotes, however, diversity and distribution in the region remain unknown, despite their significant contribution to open ocean primary production and other biogeochemical processes. To investigate how habitats created by the Amazon river plume shape small pigmented eukaryote communities, we used high-throughput sequencing of the 18S ribosomal RNA genes from up to five distinct small pigmented eukaryote cell populations, identified and sorted by flow cytometry. Small pigmented eukaryotes dominated small phytoplankton biomass across all habitat types, but the population abundances varied among stations resulting in a random distribution. Small pigmented eukaryote communities were consistently dominated by Chloropicophyceae (0.8-2 µm) and Bacillariophyceae (0.8-3.5 µm), accompanied by MOCH-5 at the surface or by Dinophyceae at the chlorophyll maximum. Taxonomic composition only displayed differences in the old plume core and at one of the plume margin stations. Such results reflect the dynamic interactions of the plume and offshore oceanic waters and suggest that the resident small pigmented eukaryote diversity was not strongly affected by habitat types at this time of the year.

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.

Temperature Acclimation of the Picoalga Ostreococcus tauri Triggers Early Fatty-Acid Variations and Involves a Plastidial ω3-Desaturase

Alteration of fatty-acid unsaturation is a universal response to temperature changes. Marine microalgae display the largest diversity of polyunsaturated fatty-acid (PUFA) whose content notably varies according to temperature. The physiological relevance and the molecular mechanisms underlying these changes are however, still poorly understood. The ancestral green picoalga Ostreococcus tauri displays original lipidic features that combines PUFAs from two distinctive microalgal lineages (Chlorophyceae, Chromista kingdom). In this study, optimized conditions were implemented to unveil early fatty-acid and desaturase transcriptional variations upon chilling and warming. We further functionally characterized the O. tauri ω3-desaturase which is closely related to ω3-desaturases from Chromista species. Our results show that the overall omega-3 to omega-6 ratio is swiftly and reversibly regulated by temperature variations. The proportion of the peculiar 18:5 fatty-acid and temperature are highly and inversely correlated pinpointing the importance of 18:5 temperature-dependent variations across kingdoms. Chilling rapidly and sustainably up-regulated most desaturase genes. Desaturases involved in the regulation of the C18-PUFA pool as well as the Δ5-desaturase appear to be major transcriptional targets. The only ω3-desaturase candidate, related to ω3-desaturases from Chromista species, is localized at chloroplasts in Nicotiana benthamiana and efficiently performs ω3-desaturation of C18-PUFAs in Synechocystis sp. PCC6803. Overexpression in the native host further unveils a broad impact on plastidial and non-plastidial glycerolipids illustrated by the alteration of omega-3/omega-6 ratio in C16-PUFA and VLC-PUFA pools. Global glycerolipid features of the overexpressor recall those of chilling acclimated cells.

Seasonal and Geographical Transitions in Eukaryotic Phytoplankton Community Structure in the Atlantic and Pacific Oceans

Much is known about how broad eukaryotic phytoplankton groups vary according to nutrient availability in marine ecosystems. However, genus- and species-level dynamics are generally unknown, although important given that adaptation and acclimation processes differentiate at these levels. We examined phytoplankton communities across seasonal cycles in the North Atlantic (BATS) and under different trophic conditions in the eastern North Pacific (ENP), using phylogenetic classification of plastid-encoded 16S rRNA amplicon sequence variants (ASVs) and other methodologies, including flow cytometric cell sorting. Prasinophytes dominated eukaryotic phytoplankton amplicons during the nutrient-rich deep-mixing winter period at BATS. During stratification (‘summer’) uncultured dictyochophytes formed ∼35 ± 10% of all surface plastid amplicons and dominated those from stramenopile algae, whereas diatoms showed only minor, ephemeral contributions over the entire year. Uncultured dictyochophytes also comprised a major fraction of plastid amplicons in the oligotrophic ENP. Phylogenetic reconstructions of near-full length 16S rRNA sequences established 11 uncultured Dictyochophyte Environmental Clades (DEC). DEC-I and DEC-VI dominated surface dictyochophytes under stratification at BATS and in the ENP, and DEC-IV was also important in the latter. Additionally, although less common at BATS, Florenciella-related clades (FC) were prominent at depth in the ENP. In both ecosystems, pelagophytes contributed notably at depth, with PEC-VIII (Pelagophyte Environmental Clade) and (cultured) Pelagomonas calceolata being most important. Q-PCR confirmed the near absence of P. calceolata at the surface of the same oligotrophic sites where it reached ∼1,500 18S rRNA gene copies ml^–1 at the DCM. To further characterize phytoplankton present in our samples, we performed staining and at-sea single-cell sorting experiments. Sequencing results from these indicated several uncultured dictyochophyte clades are comprised of predatory mixotrophs. From an evolutionary perspective, these cells showed both conserved and unique features in the chloroplast genome. In ENP metatranscriptomes we observed high expression of multiple chloroplast genes as well as expression of a selfish element (group II intron) in the psaA gene. Comparative analyses across the Pacific and Atlantic sites support the conclusion that predatory dictyochophytes thrive under low nutrient conditions. The observations that several uncultured dictyochophyte lineages are seemingly capable of photosynthesis and predation, raises questions about potential shifts in phytoplankton trophic roles associated with seasonality and long-term ocean change.

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.

A distinct lineage of giant viruses brings a rhodopsin photosystem to unicellular marine predators

Giant viruses are remarkable for their large genomes, often rivaling those of small bacteria, and for having genes thought exclusive to cellular life. Most isolated to date infect nonmarine protists, leaving their strategies and prevalence in marine environments largely unknown. Using eukaryotic single-cell metagenomics in the Pacific, we discovered a Mimiviridae lineage of giant viruses, which infects choanoflagellates, widespread protistan predators related to metazoans. The ChoanoVirus genomes are the largest yet from pelagic ecosystems, with 442 of 862 predicted proteins lacking known homologs. They are enriched in enzymes for modifying organic compounds, including degradation of chitin, an abundant polysaccharide in oceans, and they encode 3 divergent type-1 rhodopsins (VirR) with distinct evolutionary histories from those that capture sunlight in cellular organisms. One (VirR_DTS) is similar to the only other putative rhodopsin from a virus (PgV) with a known host (a marine alga). Unlike the algal virus, ChoanoViruses encode the entire pigment biosynthesis pathway and cleavage enzyme for producing the required chromophore, retinal. We demonstrate that the rhodopsin shared by ChoanoViruses and PgV binds retinal and pumps protons. Moreover, our 1.65-Å resolved VirR_DTS crystal structure and mutational analyses exposed differences from previously characterized type-1 rhodopsins, all of which come from cellular organisms. Multiple VirR types are present in metagenomes from across surface oceans, where they are correlated with and nearly as abundant as a canonical marker gene from Mimiviridae Our findings indicate that light-dependent energy transfer systems are likely common components of giant viruses of photosynthetic and phagotrophic unicellular marine eukaryotes.

Diel transcriptional response of a California Current plankton microbiome to light, low iron, and enduring viral infection

Phytoplankton and associated microbial communities provide organic carbon to oceanic food webs and drive ecosystem dynamics. However, capturing those dynamics is challenging. Here, an in situ, semi-Lagrangian, robotic sampler profiled pelagic microbes at 4 h intervals over ~2.6 days in North Pacific high-nutrient, low-chlorophyll waters. We report on the community structure and transcriptional dynamics of microbes in an operationally large size class (>5 μm) predominantly populated by dinoflagellates, ciliates, haptophytes, pelagophytes, diatoms, cyanobacteria (chiefly Synechococcus), prasinophytes (chiefly Ostreococcus), fungi, archaea, and proteobacteria. Apart from fungi and archaea, all groups exhibited 24-h periodicity in some transcripts, but larger portions of the transcriptome oscillated in phototrophs. Periodic photosynthesis-related transcripts exhibited a temporal cascade across the morning hours, conserved across diverse phototrophic lineages. Pronounced silica:nitrate drawdown, a high flavodoxin to ferredoxin transcript ratio, and elevated expression of other Fe-stress markers indicated Fe-limitation. Fe-stress markers peaked during a photoperiodically adaptive time window that could modulate phytoplankton response to seasonal Fe-limitation. Remarkably, we observed viruses that infect the majority of abundant taxa, often with total transcriptional activity synchronized with putative hosts. Taken together, these data reveal a microbial plankton community that is shaped by recycled production and tightly controlled by Fe-limitation and viral activity.

Closely related viruses of the marine picoeukaryotic alga Ostreococcus lucimarinus exhibit different ecological strategies

In marine ecosystems, viruses are major disrupters of the direct flow of carbon and nutrients to higher trophic levels. Although the genetic diversity of several eukaryotic phytoplankton virus groups has been characterized, their infection dynamics are less understood, such that the physiological and ecological implications of their diversity remain unclear. We compared genomes and infection phenotypes of the two most closely related cultured phycodnaviruses infecting the widespread picoprasinophyte Ostreococcus lucimarinus under standard- (1.3 divisions per day) and limited-light (0.41 divisions per day) nutrient replete conditions. OlV7 infection caused early arrest of the host cell cycle, coinciding with a significantly higher proportion of infected cells than OlV1-amended treatments, regardless of host growth rate. OlV7 treatments showed a near-50-fold increase of progeny virions at the higher host growth rate, contrasting with OlV1's 16-fold increase. However, production of OlV7 virions was more sensitive than OlV1 production to reduced host growth rate, suggesting fitness trade-offs between infection efficiency and resilience to host physiology. Moreover, although organic matter released from OlV1- and OlV7-infected hosts had broadly similar chemical composition, some distinct molecular signatures were observed. Collectively, these results suggest that current views on viral relatedness through marker and core gene analyses underplay operational divergence and consequences for host ecology.

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.

Green microalgae in marine coastal waters: The Ocean Sampling Day (OSD) dataset

The ecology and distribution of green phytoplankton (Chlorophyta) in the ocean is poorly known because most studies have focused on groups with large cell size such as diatoms or dinoflagellates that are easily recognized by traditional techniques such as microscopy. The Ocean Sampling Day (OSD) project sampled surface waters quasi-simultaneously at 141 marine locations, mostly in coastal waters. The analysis of the 18S V4 region OSD metabarcoding dataset reveals that Chlorophyta are ubiquitous and can be locally dominant in coastal waters. Chlorophyta represented 29% of the global photosynthetic reads (Dinoflagellates excluded) and their contribution was especially high at oligotrophic stations (up to 94%) and along the European Atlantic coast. Mamiellophyceae dominated most coastal stations. At some coastal stations, they were replaced by Chlorodendrophyceae, Ulvophyceae, Trebouxiophyceae or Chlorophyceae as the dominating group, while oligotrophic stations were dominated either by Chloropicophyceae or the uncultured prasinophytes clade IX. Several Chlorophyta classes showed preferences in terms of nitrate concentration, distance to the coast, temperature and salinity. For example, Chlorophyceae preferred cold and low salinity coastal waters, and prasinophytes clade IX warm, high salinity, oligotrophic oceanic waters.

Specialized proteomic responses and an ancient photoprotection mechanism sustain marine green algal growth during phosphate limitation

Marine algae perform approximately half of global carbon fixation, but their growth is often limited by the availability of phosphate or other nutrients^1,2. As oceans warm, the area of phosphate-limited surface waters is predicted to increase, resulting in ocean desertification^3,4. Understanding the responses of key eukaryotic phytoplankton to nutrient limitation is therefore critical^5,6. We used advanced photo-bioreactors to investigate how the widespread marine green alga Micromonas commoda grows under transitions from replete nutrients to chronic phosphate limitation and subsequent relief, analysing photosystem changes and broad cellular responses using proteomics, transcriptomics and biophysical measurements. We find that physiological and protein expression responses previously attributed to stress are critical to supporting stable exponential growth when phosphate is limiting. Unexpectedly, the abundance of most proteins involved in light harvesting does not change, but an ancient light-harvesting-related protein, LHCSR, is induced and dissipates damaging excess absorbed light as heat throughout phosphate limitation. Concurrently, a suite of uncharacterized proteins with narrow phylogenetic distributions increase multifold. Notably, of the proteins that exhibit significant changes, 70% are not differentially expressed at the mRNA transcript level, highlighting the importance of post-transcriptional processes in microbial eukaryotes. Nevertheless, transcript-protein pairs with concordant changes were identified that will enable more robust interpretation of eukaryotic phytoplankton responses in the field from metatranscriptomic studies. Our results show that P-limited Micromonas responds quickly to a fresh pulse of phosphate by rapidly increasing replication, and that the protein network associated with this ability is composed of both conserved and phylogenetically recent proteome systems that promote dynamic phosphate homeostasis. That an ancient mechanism for mitigating light stress is central to sustaining growth during extended phosphate limitation highlights the possibility of interactive effects arising from combined stressors under ocean change, which could reduce the efficacy of algal strategies for optimizing marine photosynthesis.

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.