Intron Invasions Trace Algal Speciation and Reveal Nearly Identical Arctic and Antarctic Micromonas Populations

Direct repeats found in the vicinity of intron splice sites

Four main classes of introns (group I, group II, spliceosomal, and archaeal) have been reported for all major types of RNA from nuclei and organelles of a wide range of taxa. When and how introns inserted within the genic regions of genomes, however, is often unclear. Introns were examined from Archaea, Bacteria, and Eukarya. Up to 80 bp surrounding each of the 5' and 3' intron/exon borders were compared to search for direct repeats (DRs). For each of the 213 introns examined, DNA sequence analysis revealed DRs at or near the intron/exon borders, ranging from 4 to 30 bp in length, with a mean of 11.4 bp. More than 80% of the repeats were within 10 bp of the intron/exon borders. The numbers of DRs 6-30 bp in length were greater than expected by chance. When a DNA segment moves into a new genomic location, the insertion involves a double-strand DNA break that must be repaired to maintain genome stability and often results in a pair of DRs that now flank the insert. This insertion process applies to both mobile genetic elements (MGEs), such as transposons, and to introns as reported here. The DNA break at the insertion site may be caused by transposon-like events or recombination. Thus, introns and transposons appear to be members of a group of parasitic MGEs that secondarily may benefit their host cell and have expanded greatly in eukaryotes from their prokaryotic ancestors.

Intron lariat spliceosomes convert lariats to true circles: implications for intron transposition

Rare, full-length circular intron RNAs distinct from lariats have been reported in several species, but their biogenesis is not understood. We envisioned and tested a hypothesis for their formation using Saccharomyces cerevisiae, documenting full-length and novel processed circular RNAs from multiple introns. Evidence implicates a previously undescribed catalytic activity of the intron lariat spliceosome (ILS) in which the 3'-OH of the lariat tail (with optional trimming and adenylation by the nuclear 3' processing machinery) attacks the branch, joining the intron 3' end to the 5' splice site in a 3'-5' linked circle. Human U2 and U12 spliceosomes produce analogous full-length and processed circles. Postsplicing catalytic activity of the spliceosome may promote intron transposition during eukaryotic genome evolution.

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.

Intron-lariat spliceosomes convert lariats to true circles: implications for intron transposition

Rare, full length circular intron RNAs distinct from lariats have been reported in several species, but their biogenesis is not understood. We envision and test a hypothesis for their formation using Saccharomyces cerevisiae, documenting full length and novel processed circular RNAs from multiple introns. Evidence implicates a previously undescribed catalytic activity of the intron-lariat spliceosome (ILS) in which the 3'-OH of the lariat tail (with optional trimming and adenylation by the nuclear 3' processing machinery) attacks the branch, joining the intron 3' end to the 5' splice site in a 3'-5' linked circle. Human U2 and U12 spliceosomes produce analogous full length and processed circles. Post-splicing catalytic activity of the spliceosome may promote intron transposition during eukaryotic genome evolution.

Transposable elements drive intron gain in diverse eukaryotes

There is massive variation in intron numbers across eukaryotic genomes, yet the major drivers of intron content during evolution remain elusive. Rapid intron loss and gain in some lineages contrast with long-term evolutionary stasis in others. Episodic intron gain could be explained by recently discovered specialized transposons called Introners, but so far Introners are only known from a handful of species. Here, we performed a systematic search across 3,325 eukaryotic genomes and identified 27,563 Introner-derived introns in 175 genomes (5.2%). Species with Introners span remarkable phylogenetic diversity, from animals to basal protists, representing lineages whose last common ancestor dates to over 1.7 billion years ago. Aquatic organisms were 6.5 times more likely to contain Introners than terrestrial organisms. Introners exhibit mechanistic diversity but most are consistent with DNA transposition, indicating that Introners have evolved convergently hundreds of times from nonautonomous transposable elements. Transposable elements and aquatic taxa are associated with high rates of horizontal gene transfer, suggesting that this combination of factors may explain the punctuated and biased diversity of species containing Introners. More generally, our data suggest that Introners may explain the episodic nature of intron gain across the eukaryotic tree of life. These results illuminate the major source of ongoing intron creation in eukaryotic genomes.

Intron losses and gains in the nematodes

BACKGROUND

The evolution of spliceosomal introns has been widely studied among various eukaryotic groups. Researchers nearly reached the consensuses on the pattern and the mechanisms of intron losses and gains across eukaryotes. However, according to previous studies that analyzed a few genes or genomes, Nematoda seems to be an eccentric group.

RESULTS

Taking advantage of the recent accumulation of sequenced genomes, we extensively analyzed the intron losses and gains using 104 nematode genomes across all the five Clades of the phylum. Nematodes have a wide range of intron density, from less than one to more than nine per kbp coding sequence. The rates of intron losses and gains exhibit significant heterogeneity both across different nematode lineages and across different evolutionary stages of the same lineage. The frequency of intron losses far exceeds that of intron gains. Five pieces of evidence supporting the model of cDNA-mediated intron loss have been observed in ten Caenorhabditis species, the dominance of the precise intron losses, frequent loss of adjacent introns, high-level expression of the intron-lost genes, preferential losses of short introns, and the preferential losses of introns close to 3'-ends of genes. Like studies in most eukaryotic groups, we cannot find the source sequences for the limited number of intron gains detected in the Caenorhabditis genomes.

CONCLUSIONS

These results indicate that nematodes are a typical eukaryotic group rather than an outlier in intron evolution.

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.

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.

Spatiotemporal Variations in Antarctic Protistan Communities Highlight Phytoplankton Diversity and Seasonal Dominance by a Novel Cryptophyte Lineage

The Andvord fjord in the West Antarctic Peninsula (WAP) is known for its productivity and abundant megafauna. Nevertheless, seasonal patterns of the molecular diversity and abundance of protistan community members underpinning WAP productivity remain poorly resolved. We performed spring and fall expeditions pursuing protistan diversity, abundance of photosynthetic taxa, and the connection to changing conditions. 18S rRNA amplicon sequence variant (ASV) profiles revealed diverse predatory protists spanning multiple eukaryotic supergroups, alongside enigmatic heterotrophs like the Picozoa. Among photosynthetic protists, cryptophyte contributions were notable. Analysis of plastid-derived 16S rRNA ASVs supported 18S ASV results, including a dichotomy between cryptophytes and diatom contributions previously reported in other Antarctic regions. We demonstrate that stramenopile and cryptophyte community structures have distinct attributes. Photosynthetic stramenopiles exhibit high diversity, with the polar diatom Fragilariopsis cylindrus, unidentified Chaetoceros species, and others being prominent. Conversely, ASV analyses followed by environmental full-length rRNA gene sequencing, electron microscopy, and flow cytometry revealed that a novel alga dominates the cryptophytes. Phylogenetic analyses established that TPG clade VII, as named here, is evolutionarily distinct from cultivated cryptophyte lineages. Additionally, cryptophyte cell abundance correlated with increased water temperature. Analyses of global data sets showed that clade VII dominates cryptophyte ASVs at Southern Ocean sites and appears to be endemic, whereas in the Arctic and elsewhere, Teleaulax amphioxeia and Plagioselmis prolonga dominate, although both were undetected in Antarctic waters. Collectively, our studies provide baseline data against which future change can be assessed, identify different diversification patterns between stramenopiles and cryptophytes, and highlight an evolutionarily distinct cryptophyte clade that thrives under conditions enhanced by warming.

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

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

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.

Internally Symmetrical Stwintrons and Related Canonical Introns in Hypoxylaceae Species

Spliceosomal introns are pervasive in eukaryotes. Intron gains and losses have occurred throughout evolution, but the origin of new introns is unclear. Stwintrons are complex intervening sequences where one of the sequence elements (5′-donor, lariat branch point element or 3′-acceptor) necessary for excision of a U2 intron (external intron) is itself interrupted by a second (internal) U2 intron. In Hypoxylaceae, a family of endophytic fungi, we uncovered scores of donor-disrupted stwintrons with striking sequence similarity among themselves and also with canonical introns. Intron–exon structure comparisons suggest that these stwintrons have proliferated within diverging taxa but also give rise to proliferating canonical introns in some genomes. The proliferated (stw)introns have integrated seamlessly at novel gene positions. The recently proliferated (stw)introns appear to originate from a conserved ancestral stwintron characterised by terminal inverted repeats (45–55 nucleotides), a highly symmetrical structure that may allow the formation of a double-stranded intron RNA molecule. No short tandem duplications flank the putatively inserted intervening sequences, which excludes a DNA transposition-based mechanism of proliferation. It is tempting to suggest that this highly symmetrical structure may have a role in intron proliferation by (an)other mechanism(s).

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.

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.

Controlled sampling of ribosomally active protistan diversity in sediment-surface layers identifies putative players in the marine carbon sink

Marine sediments are one of the largest carbon reservoir on Earth, yet the microbial communities, especially the eukaryotes, that drive these ecosystems are poorly characterised. Here, we report implementation of a sampling system that enables injection of reagents into sediments at depth, allowing for preservation of RNA in situ. Using the RNA templates recovered, we investigate the 'ribosomally active' eukaryotic diversity present in sediments close to the water/sediment interface. We demonstrate that in situ preservation leads to recovery of a significantly altered community profile. Using SSU rRNA amplicon sequencing, we investigated the community structure in these environments, demonstrating a wide diversity and high relative abundance of stramenopiles and alveolates, specifically: Bacillariophyta (diatoms), labyrinthulomycetes and ciliates. The identification of abundant diatom rRNA molecules is consistent with microscopy-based studies, but demonstrates that these algae can also be exported to the sediment as active cells as opposed to dead forms. We also observe many groups that include, or branch close to, osmotrophic-saprotrophic protists (e.g. labyrinthulomycetes and Pseudofungi), microbes likely to be important for detrital decomposition. The sequence data also included a diversity of abundant amplicon-types that branch close to the Fonticula slime moulds. Taken together, our data identifies additional roles for eukaryotic microbes in the marine carbon cycle; where putative osmotrophic-saprotrophic protists represent a significant active microbial-constituent of the upper sediment layer.

Rapidly evolving protointrons in Saccharomyces genomes revealed by a hungry spliceosome

Introns are a prevalent feature of eukaryotic genomes, yet their origins and contributions to genome function and evolution remain mysterious. In budding yeast, repression of the highly transcribed intron-containing ribosomal protein genes (RPGs) globally increases splicing of non-RPG transcripts through reduced competition for the spliceosome. We show that under these “hungry spliceosome” conditions, splicing occurs at more than 150 previously unannotated locations we call protointrons that do not overlap known introns. Protointrons use a less constrained set of splice sites and branchpoints than standard introns, including in one case AT-AC in place of GT-AG. Protointrons are not conserved in all closely related species, suggesting that most are not under positive selection and are fated to disappear. Some are found in non-coding RNAs (e. g. CUTs and SUTs), where they may contribute to the creation of new genes. Others are found across boundaries between noncoding and coding sequences, or within coding sequences, where they offer pathways to the creation of new protein variants, or new regulatory controls for existing genes. We define protointrons as (1) nonconserved intron-like sequences that are (2) infrequently spliced, and importantly (3) are not currently understood to contribute to gene expression or regulation in the way that standard introns function. A very few protointrons in S. cerevisiae challenge this classification by their increased splicing frequency and potential function, consistent with the proposed evolutionary process of “intronization”, whereby new standard introns are created. This snapshot of intron evolution highlights the important role of the spliceosome in the expansion of transcribed genomic sequence space, providing a pathway for the rare events that may lead to the birth of new eukaryotic genes and the refinement of existing gene function.

The protein coding information in eukaryotic genes is broken by intervening sequences called introns that are removed from RNA during transcription by a large protein-RNA complex called the spliceosome. Where introns come from and how the spliceosome contributes to genome evolution are open questions. In this study, we find more than 150 new places in the yeast genome that are recognized by the spliceosome and spliced out as introns. Since they appear to have arisen very recently in evolution by sequence drift and do not appear to contribute to gene expression or its regulation, we call these protointrons. Protointrons are found in both protein-coding and non-coding RNAs and are not efficiently removed by the splicing machinery. Although most protointrons are not conserved and will likely disappear as evolution proceeds, a few are spliced more efficiently, and are located where they might begin to play functional roles in gene expression, as predicted by the proposed process of intronization. The challenge now is to understand how spontaneously appearing splicing events like protointrons might contribute to the creation of new genes, new genetic controls, and new protein isoforms as genomes evolve.

Rapidly evolving protointrons in Saccharomyces genomes revealed by a hungry spliceosome

Introns are a prevalent feature of eukaryotic genomes, yet their origins and contributions to genome function and evolution remain mysterious. In budding yeast, repression of the highly transcribed intron-containing ribosomal protein genes (RPGs) globally increases splicing of non-RPG transcripts through reduced competition for the spliceosome. We show that under these "hungry spliceosome" conditions, splicing occurs at more than 150 previously unannotated locations we call protointrons that do not overlap known introns. Protointrons use a less constrained set of splice sites and branchpoints than standard introns, including in one case AT-AC in place of GT-AG. Protointrons are not conserved in all closely related species, suggesting that most are not under positive selection and are fated to disappear. Some are found in non-coding RNAs (e. g. CUTs and SUTs), where they may contribute to the creation of new genes. Others are found across boundaries between noncoding and coding sequences, or within coding sequences, where they offer pathways to the creation of new protein variants, or new regulatory controls for existing genes. We define protointrons as (1) nonconserved intron-like sequences that are (2) infrequently spliced, and importantly (3) are not currently understood to contribute to gene expression or regulation in the way that standard introns function. A very few protointrons in S. cerevisiae challenge this classification by their increased splicing frequency and potential function, consistent with the proposed evolutionary process of "intronization", whereby new standard introns are created. This snapshot of intron evolution highlights the important role of the spliceosome in the expansion of transcribed genomic sequence space, providing a pathway for the rare events that may lead to the birth of new eukaryotic genes and the refinement of existing gene function.

The evolution of the plastid phosphate translocator family

The plastid phosphate translocators evolved in algae but diversified into several groups, which adopted different physiological functions by extensive gene duplications and losses in Streptophyta. The plastid phosphate translocators (pPT) are a family of transporters involved in the exchange of metabolites and inorganic phosphate between stroma and cytosol. Based on their substrate specificities, they were divided into four subfamilies named TPT, PPT, GPT and XPT. To analyse the occurrence of these transporters in different algae and land plant species, we identified 652 pPT genes in 101 sequenced genomes for phylogenetic analysis. The first three subfamilies are found in all species and evolved before the split of red and green algae while the XPTs were derived from the duplication of a GPT gene at the base of Streptophyta. The analysis of the intron-exon structures of the pPTs corroborated these findings. While the number and positions of introns are conserved within each subfamily, they differ between the subfamilies suggesting an insertion of the introns shortly after the three subfamilies evolved. During angiosperm evolution, the subfamilies further split into different groups (TPT1-2, PPT1-3, GPT1-6). Angiosperm species differ significantly in the total number of pPTs, with many species having only a few, while several plants, especially crops, have a higher number, pointing to the importance of these transporters for improved source-sink strength and yield. The differences in the number of pPTs can be explained by several small-scale gene duplications and losses in plant families or single species, but also by whole genome duplications, for example, in grasses. This work could be the basis for a comprehensive analysis of the molecular and physiological functions of this important family of transporters.

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.

Genome streamlining via complete loss of introns has occurred multiple times in lichenized fungal mitochondria

Reductions in genome size and complexity are a hallmark of obligate symbioses. The mitochondrial genome displays clear examples of these reductions, with the ancestral alpha-proteobacterial genome size and gene number having been reduced by orders of magnitude in most descendent modern mitochondrial genomes. Here, we examine patterns of mitochondrial evolution specifically looking at intron size, number, and position across 58 species from 21 genera of lichenized Ascomycete fungi, representing a broad range of fungal diversity and niches. Our results show that the cox1gene always contained the highest number of introns out of all the mitochondrial protein-coding genes, that high intron sequence similarity (>90%) can be maintained between different genera, and that lichens have undergone at least two instances of complete, genome-wide intron loss consistent with evidence for genome streamlining via loss of parasitic, noncoding DNA, in Phlyctis boliviensisand Graphis lineola. Notably, however, lichenized fungi have not only undergone intron loss but in some instances have expanded considerably in size due to intron proliferation (e.g., Alectoria fallacina and Parmotrema neotropicum), even between closely related sister species (e.g., Cladonia). These results shed light on the highly dynamic mitochondrial evolution that is occurring in lichens and suggest that these obligate symbiotic organisms are in some cases undergoing recent, broad-scale genome streamlining via loss of protein-coding genes as well as noncoding, parasitic DNA elements.

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.

Transcriptional responses of the marine green alga Micromonas pusilla and an infecting prasinovirus under different phosphate conditions

Prasinophytes are widespread marine algae for which responses to nutrient limitation and viral infection are not well understood. We studied the picoprasinophyte, Micromonas pusilla, grown under phosphate-replete (0.65 ± 0.07 d^-1 ) and 10-fold lower (low)-phosphate (0.11 ± 0.04 d^-1 ) conditions, and infected by the phycodnavirus MpV-SP1. Expression of 17% of Micromonas genes in uninfected cells differed by >1.5-fold (q < 0.01) between nutrient conditions, with genes for P-metabolism and the uniquely-enriched Sel1-like repeat (SLR) family having higher relative transcript abundances, while phospholipid-synthesis genes were lower in low-P than P-replete. Approximately 70% (P-replete) and 30% (low-P) of cells were lysed 24 h post-infection, and expression of ≤5.8% of host genes changed relative to uninfected treatments. Host genes for CAZymes and glycolysis were activated by infection, supporting importance in viral production, which was significantly lower in slower growing (low-P) hosts. All MpV-SP1 genes were expressed, and our analyses suggest responses to differing host-phosphate backgrounds involve few viral genes, while the temporal program of infection involves many more, and is largely independent of host-phosphate background. Our study (i) identifies genes previously unassociated with nutrient acclimation or viral infection, (ii) provides insights into the temporal program of prasinovirus gene expression by hosts and (iii) establishes cell biological aspects of an ecologically important host-prasinovirus system that differ from other marine algal-virus systems.

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.

Origin and Spread of Spliceosomal Introns: Insights from the Fungal Clade Zymoseptoria

Spliceosomal introns are a key feature of eukaryote genome architecture and have been proposed to originate from selfish group II introns from an endosymbiotic bacterium, that is, the ancestor of mitochondria. However, the mechanisms underlying the wide spread of spliceosomal introns across eukaryotic genomes have been obscure. In this study, we characterize the dynamic evolution of spliceosomal introns in the fungal genus Zymoseptoria at different evolutionary scales, that is, within a genome, among conspecific strains within species, and between different species. Within the genome, spliceosomal introns can proliferate in unrelated genes and intergenic regions. Among conspecific strains, spliceosomal introns undergo rapid turnover (gains and losses) and frequent sequence exchange between geographically distinct strains. Furthermore, spliceosomal introns could undergo introgression between distinct species, which can further promote intron invasion and proliferation. The dynamic invasion and proliferation processes of spliceosomal introns resemble the life cycles of mobile selfish (group I/II) introns, and these intron movements, at least in part, account for the dramatic processes of intron gain and intron loss during eukaryotic evolution.

Photosynthetic Picoeukaryotes in the Land-Fast Ice of the White Sea, Russia

The White Sea is a unique marine environment combining features of temperate and Arctic seas. The composition and abundance of photosynthetic picoeukaryotes (PPEs) were investigated in the land-fast ice of the White Sea, Russia, in March 2013 and 2014. High-throughput tag sequencing (Illumina MiSeq system) of the V4 region of the 18S rRNA gene was used to reveal the diversity of PPE ice community. The integrated PPE abundance varied from 11 × 10⁶ cells/m² to 364 × 10⁶ cells/m²; the integrated biomass ranged from 0.02 to 0.26 mg С/m². The composition of sea-ice PPEs was represented by 16 algae genera belonging to eight classes and three super-groups. Chlorophyta, especially Mamiellophyceae, dominated among ice PPEs. The detailed analysis revealed the latent diversity of Micromonas and Mantоniella. Micromonas clade E2 revealed in the subarctic White Sea ice indicates that the area of distribution of this species is wider than previously thought. We suppose there exists a new Micromonas clade F. Micromonas clade C and Minutocellulus polymorphus were first discovered in the ice and extend the modern concept of sympagic communities' diversity generally and highlights the importance of further targeting subarctic sea ice for microbial study.

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.

Seasonal patterns in Arctic prasinophytes and inferred ecology of Bathycoccus unveiled in an Arctic winter metagenome

Prasinophytes occur in all oceans but rarely dominate phytoplankton populations. In contrast, a single ecotype of the prasinophyte Micromonas is frequently the most abundant photosynthetic taxon reported in the Arctic from summer through autumn. However, seasonal dynamics of prasinophytes outside of this period are little known. To address this, we analyzed high-throughput V4 18S rRNA amplicon data collected from November to July in the Amundsen Gulf Region, Beaufort Sea, Arctic. Surprisingly during polar sunset in November and December, we found a high proportion of reads from both DNA and RNA belonging to another prasinophyte, Bathycoccus. We then analyzed a metagenome from a December sample and the resulting Bathycoccus metagenome assembled genome (MAG) covered ~90% of the Bathycoccus Ban7 reference genome. In contrast, only ~20% of a reference Micromonas genome was found in the metagenome. Our phylogenetic analysis of marker genes placed the Arctic Bathycoccus in the B1 coastal clade. In addition, substitution rates of 129 coding DNA sequences were ~1.6% divergent between the Arctic MAG and coastal Chilean upwelling MAGs and 17.3% between it and a South East Atlantic open ocean MAG in the B2 Clade. The metagenomic analysis also revealed a winter viral community highly skewed toward viruses targeting Micromonas, with a much lower diversity of viruses targeting Bathycoccus. Overall a combination of Micromonas being relatively less able to maintain activity under dark winter conditions and viral suppression of Micromonas may have contributed to the success of Bathycoccus in the Amundsen Gulf during winter.

Responses of the picoprasinophyte Micromonas commoda to light and ultraviolet stress

Micromonas is a unicellular marine green alga that thrives from tropical to polar ecosystems. We investigated the growth and cellular characteristics of acclimated mid-exponential phase Micromonas commoda RCC299 over multiple light levels and over the diel cycle (14:10 hour light:dark). We also exposed the light:dark acclimated M. commoda to experimental shifts from moderate to high light (HL), and to HL plus ultraviolet radiation (HL+UV), 4.5 hours into the light period. Cellular responses of this prasinophyte were quantified by flow cytometry and changes in gene expression by qPCR and RNA-seq. While proxies for chlorophyll a content and cell size exhibited similar diel variations in HL and controls, with progressive increases during day and decreases at night, both parameters sharply decreased after the HL+UV shift. Two distinct transcriptional responses were observed among chloroplast genes in the light shift experiments: i) expression of transcription and translation-related genes decreased over the time course, and this transition occurred earlier in treatments than controls; ii) expression of several photosystem I and II genes increased in HL relative to controls, as did the growth rate within the same diel period. However, expression of these genes decreased in HL+UV, likely as a photoprotective mechanism. RNA-seq also revealed two genes in the chloroplast genome, ycf2-like and ycf1-like, that had not previously been reported. The latter encodes the second largest chloroplast protein in Micromonas and has weak homology to plant Ycf1, an essential component of the plant protein translocon. Analysis of several nuclear genes showed that the expression of LHCSR2, which is involved in non-photochemical quenching, and five light-harvesting-like genes, increased 30 to >50-fold in HL+UV, but was largely unchanged in HL and controls. Under HL alone, a gene encoding a novel nitrite reductase fusion protein (NIRFU) increased, possibly reflecting enhanced N-assimilation under the 625 μmol photons m^-2 s^-1 supplied in the HL treatment. NIRFU’s domain structure suggests it may have more efficient electron transfer than plant NIR proteins. Our analyses indicate that Micromonas can readily respond to abrupt environmental changes, such that strong photoinhibition was provoked by combined exposure to HL and UV, but a ca. 6-fold increase in light was stimulatory.

Survey of the green picoalga Bathycoccus genomes in the global ocean

Bathycoccus is a cosmopolitan green micro-alga belonging to the Mamiellophyceae, a class of picophytoplankton that contains important contributors to oceanic primary production. A single species of Bathycoccus has been described while the existence of two ecotypes has been proposed based on metagenomic data. A genome is available for one strain corresponding to the described phenotype. We report a second genome assembly obtained by a single cell genomics approach corresponding to the second ecotype. The two Bathycoccus genomes are divergent enough to be unambiguously distinguishable in whole DNA metagenomic data although they possess identical sequence of the 18S rRNA gene including in the V9 region. Analysis of 122 global ocean whole DNA metagenome samples from the Tara-Oceans expedition reveals that populations of Bathycoccus that were previously identified by 18S rRNA V9 metabarcodes are only composed of these two genomes. Bathycoccus is relatively abundant and widely distributed in nutrient rich waters. The two genomes rarely co-occur and occupy distinct oceanic niches in particular with respect to depth. Metatranscriptomic data provide evidence for gain or loss of highly expressed genes in some samples, suggesting that the gene repertoire is modulated by environmental conditions.

Mechanism for DNA transposons to generate introns on genomic scales

The discovery of introns four decades ago was one of the most unexpected findings in molecular biology. Introns are sequences interrupting genes that must be removed as part of messenger RNA production. Genome sequencing projects have shown that most eukaryotic genes contain at least one intron, and frequently many. Comparison of these genomes reveals a history of long evolutionary periods during which few introns were gained, punctuated by episodes of rapid, extensive gain. However, although several detailed mechanisms for such episodic intron generation have been proposed, none has been empirically supported on a genomic scale. Here we show how short, non-autonomous DNA transposons independently generated hundreds to thousands of introns in the prasinophyte Micromonas pusilla and the pelagophyte Aureococcus anophagefferens. Each transposon carries one splice site. The other splice site is co-opted from the gene sequence that is duplicated upon transposon insertion, allowing perfect splicing out of the RNA. The distributions of sequences that can be co-opted are biased with respect to codons, and phasing of transposon-generated introns is similarly biased. These transposons insert between pre-existing nucleosomes, so that multiple nearby insertions generate nucleosome-sized intervening segments. Thus, transposon insertion and sequence co-option may explain the intron phase biases and prevalence of nucleosome-sized exons observed in eukaryotes. Overall, the two independent examples of proliferating elements illustrate a general DNA transposon mechanism that can plausibly account for episodes of rapid, extensive intron gain during eukaryotic evolution.

Intron gain by tandem genomic duplication: a novel case in a potato gene encoding RNA-dependent RNA polymerase

The origin and subsequent accumulation of spliceosomal introns are prominent events in the evolution of eukaryotic gene structure. However, the mechanisms underlying intron gain remain unclear because there are few proven cases of recently gained introns. In an RNA-dependent RNA polymerase (RdRp) gene, we found that a tandem duplication occurred after the divergence of potato and its wild relatives among other Solanum plants. The duplicated sequence crosses the intron-exon boundary of the first intron and the second exon. A new intron was detected at this duplicated region, and it includes a small previously exonic segment of the upstream copy of the duplicated sequence and the intronic segment of the downstream copy of the duplicated sequence. The donor site of this new intron was directly obtained from the small previously exonic segment. Most of the splicing signals were inherited directly from the parental intron/exon structure, including a putative branch site, the polypyrimidine tract, the 3' splicing site, two putative exonic splicing enhancers, and the GC contents differed between the intron and exon. In the widely cited model of intron gain by tandem genomic duplication, the duplication of an AGGT-containing exonic segment provides the GT and AG splicing sites for the new intron. Our results illustrate that the tandem duplication model of intron gain should be diverse in terms of obtaining the proper splicing signals.

Phylogenetic diversity and biogeography of the Mamiellophyceae lineage of eukaryotic phytoplankton across the oceans

High-throughput diversity amplicon sequencing of marine microbial samples has revealed that members of the Mamiellophyceae lineage are successful phytoplankton in many oceanic habitats. Indeed, these eukaryotic green algae can dominate the picoplanktonic biomass, however, given the broad expanses of the oceans, their geographical distributions and the phylogenetic diversity of some groups remain poorly characterized. As these algae play a foundational role in marine food webs, it is crucial to assess their global distribution in order to better predict potential changes in abundance and community structure. To this end, we analyzed the V9-18S small subunit rDNA sequences deposited from the Tara Oceans expedition to evaluate the diversity and biogeography of these phytoplankton. Our results show that the phylogenetic composition of Mamiellophyceae communities is in part determined by geographical provenance, and do not appear to be influenced - in the samples recovered - by water depth, at least at the resolution possible with the V9-18S. Phylogenetic classification of Mamiellophyceae sequences revealed that the Dolichomastigales order encompasses more sequence diversity than other orders in this lineage. These results indicate that a large fraction of the Mamiellophyceae diversity has been hitherto overlooked, likely because of a combination of size fraction, sequencing and geographical limitations.

Spliceosomal intronogenesis

The presence of intervening sequences, termed introns, is a defining characteristic of eukaryotic nuclear genomes. Once transcribed into pre-mRNA, these introns must be removed within the spliceosome before export of the processed mRNA to the cytoplasm, where it is translated into protein. Although intron loss has been demonstrated experimentally, several mysteries remain regarding the origin and propagation of introns. Indeed, documented evidence of gain of an intron has only been suggested by phylogenetic analyses. We report the use of a strategy that detects selected intron gain and loss events. We have experimentally verified, to our knowledge, the first demonstrations of intron transposition in any organism. From our screen, we detected two separate intron gain events characterized by the perfect transposition of a reporter intron into the yeast genes RPL8B and ADH2, respectively. We show that the newly acquired introns are able to be removed from their respective pre-mRNAs by the spliceosome. Additionally, the novel allele, RPL8Bint, is functional when overexpressed within the genome in a strain lacking the Rpl8 paralogue RPL8A, demonstrating that the gene targeted for intronogenesis is functional.

Evidence-based green algal genomics reveals marine diversity and ancestral characteristics of land plants

Background

Prasinophytes are widespread marine green algae that are related to plants. Cellular abundance of the prasinophyte Micromonas has reportedly increased in the Arctic due to climate-induced changes. Thus, studies of these unicellular eukaryotes are important for marine ecology and for understanding Viridiplantae evolution and diversification.

Results

We generated evidence-based Micromonas gene models using proteomics and RNA-Seq to improve prasinophyte genomic resources. First, sequences of four chromosomes in the 22 Mb Micromonas pusilla (CCMP1545) genome were finished. Comparison with the finished 21 Mb genome of Micromonas commoda (RCC299; named herein) shows they share ≤8,141 of ~10,000 protein-encoding genes, depending on the analysis method. Unlike RCC299 and other sequenced eukaryotes, CCMP1545 has two abundant repetitive intron types and a high percent (26 %) GC splice donors. Micromonas has more genus-specific protein families (19 %) than other genome sequenced prasinophytes (11 %). Comparative analyses using predicted proteomes from other prasinophytes reveal proteins likely related to scale formation and ancestral photosynthesis. Our studies also indicate that peptidoglycan (PG) biosynthesis enzymes have been lost in multiple independent events in select prasinophytes and plants. However, CCMP1545, polar Micromonas CCMP2099 and prasinophytes from other classes retain the entire PG pathway, like moss and glaucophyte algae. Surprisingly, multiple vascular plants also have the PG pathway, except the Penicillin-Binding Protein, and share a unique bi-domain protein potentially associated with the pathway. Alongside Micromonas experiments using antibiotics that halt bacterial PG biosynthesis, the findings highlight unrecognized phylogenetic complexity in PG-pathway retention and implicate a role in chloroplast structure or division in several extant Viridiplantae lineages.

Conclusions

Extensive differences in gene loss and architecture between related prasinophytes underscore their divergence. PG biosynthesis genes from the cyanobacterial endosymbiont that became the plastid, have been selectively retained in multiple plants and algae, implying a biological function. Our studies provide robust genomic resources for emerging model algae, advancing knowledge of marine phytoplankton and plant evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2585-6) contains supplementary material, which is available to authorized users.

Abundance and Biogeography of Picoprasinophyte Ecotypes and Other Phytoplankton in the Eastern North Pacific Ocean

Eukaryotic algae within the picoplankton size class (≤2 μm in diameter) are important marine primary producers, but their spatial and ecological distributions are not well characterized. Here, we studied three picoeukaryotic prasinophyte genera and their cyanobacterial counterparts, Prochlorococcus and Synechococcus, during two cruises along a North Pacific transect characterized by different ecological regimes. Picoeukaryotes and Synechococcus reached maximum abundances of 1.44 × 10(5) and 3.37 × 10(5) cells · ml(-1), respectively, in mesotrophic waters, while Prochlorococcus reached 1.95 × 10(5) cells · ml(-1) in the oligotrophic ocean. Of the picoeukaryotes, Bathycoccus was present at all stations in both cruises, reaching 21,368 ± 327 18S rRNA gene copies · ml(-1). Micromonas and Ostreococcus clade OI were detected only in mesotrophic and coastal waters and Ostreococcus clade OII only in the oligotrophic ocean. To resolve proposed Bathycoccus ecotypes, we established genetic distances for 1,104 marker genes using targeted metagenomes and the Bathycoccus prasinos genome. The analysis was anchored in comparative genome analysis of three Ostreococcus species for which physiological and environmental data are available to facilitate data interpretation. We established that two Bathycoccus ecotypes exist, named here BI (represented by coastal isolate Bathycoccus prasinos) and BII. These share 82% ± 6% nucleotide identity across homologs, while the Ostreococcus spp. share 75% ± 8%. We developed and applied an analysis of ecomarkers to metatranscriptomes sequenced here and published -omics data from the same region. The results indicated that the Bathycoccus ecotypes cooccur more often than Ostreococcus clades OI and OII do. Exploratory analyses of relative transcript abundances suggest that Bathycoccus NRT2.1 and AMT2.2 are high-affinity NO3 (-) and low-affinity NH4 (+) transporters, respectively, with close homologs in multiple picoprasinophytes. Additionally, in the open ocean, where dissolved iron concentrations were low (0.08 nM), there appeared to be a shift to the use of nickel superoxide dismutases (SODs) from Mn/Fe/Cu SODs closer inshore. Our study documents the distribution of picophytoplankton along a North Pacific ecological gradient and offers new concepts and techniques for investigating their biogeography.