# Genomic adaptation of the picoeukaryote *Pelagomonas calceolata* to iron-poor oceans revealed by a chromosome-scale genome sequence

Nina Guérin, Marta Ciccarella, Elisa Flamant, Paul Frémont, Sophie Mangenot, Benjamin Istace, Benjamin Noel, Caroline Belser, Laurie Bertrand, Karine Labadie, et al.

► To cite this version:

Nina Guérin, Marta Ciccarella, Elisa Flamant, Paul Frémont, Sophie Mangenot, et al.. Genomic adaptation of the picoeukaryote *Pelagomonas calceolata* to iron-poor oceans revealed by a chromosome-scale genome sequence. Communications Biology, 2022, 5 (1), pp.983. 10.1038/s42003-022-03939-z.
cea-03791710

HAL Id: cea-03791710

https://cea.hal.science/cea-03791710

Submitted on 29 Sep 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Genomic adaptation of the picoeukaryote *Pelagomonas calceolata* to iron-poor oceans revealed by a chromosome-scale genome sequence

Nina Guérin<fnref n="1" />,², Marta Ciccarella ID<fnref n="1" />, Elisa Flamant<fnref n="1" />,², Paul Frémont ID<fnref n="1" />,², Sophie Mangenot<fnref n="1" />,², Benjamin Istace<fnref n="1" />, Benjamin Noel ID<fnref n="1" />, Caroline Belser ID<fnref n="1" />, Laurie Bertrand<fnref n="1" />,², Karine Labadie ID²,³, Corinne Cruaud²,³, Sarah Romac ID⁴, Charles Bachy ID⁴,⁵, Martin Gachenot⁵, Eric Pelletier<fnref n="1" />,², Adriana Alberti<fnref n="1" />,²,⁶, Olivier Jaillon ID<fnref n="1" />,², Patrick Wincker ID<fnref n="1" />,², Jean-Marc Aury ID<fnref n="1" /> & Quentin Carradec ID<fnref n="1" />,²✉

The smallest phytoplankton species are key actors in oceans biogeochemical cycling and their abundance and distribution are affected with global environmental changes. Among them, algae of the Pelagophyceae class encompass coastal species causative of harmful algal blooms while others are cosmopolitan and abundant. The lack of genomic reference in this lineage is a main limitation to study its ecological importance. Here, we analysed *Pelagomonas calceolata* relative abundance, ecological niche and potential for the adaptation in all oceans using a complete chromosome-scale assembled genome sequence. Our results show that *P. calceolata* is one of the most abundant eukaryotic species in the oceans with a relative abundance favoured by high temperature, low-light and iron-poor conditions. Climate change projections based on its relative abundance suggest an extension of the *P. calceolata* habitat toward the poles at the end of this century. Finally, we observed a specific gene repertoire and expression level variations potentially explaining its ecological success in low-iron and low-nitrate environments. Collectively, these findings reveal the ecological importance of *P. calceolata* and lay the foundation for a global scale analysis of the adaptation and acclimation strategies of this small phytoplankton in a changing environment.

<footnote n="1">¹Génomique Métabolique, Genoscope, Institut François Jacob, CEA, CNRS, Univ Evry, Université Paris-Saclay, 91057 Evry, France. ²Research Federation for the Study of Global Ocean Systems Ecology and Evolution, R2022/Tara Oceans GO-SEE, 3 rue Michel-Ange, 75016 Paris, France. ³Genoscope, Institut François Jacob, CEA, Université Paris-Saclay, 2 Rue Gaston Crémieux, 91057 Evry, France. ⁴Sorbonne Université, CNRS, Station Biologique de Roscoff, AD2M, UMR7144, Place Georges Teissier, 29680 Roscoff, France. ⁵Sorbonne Université, CNRS, FR2424, Station Biologique de Roscoff, 29680 Roscoff, France. ⁶Université Paris-Saclay, CEA, CNRS, Institute for Integrative Biology of the Cell (I2BC), 91198 Gif-sur-Yvette, France. ✉email: qcarrade@genoscope.cns.fr</footnote>
Marine phytoplankton accounts for more than 45% of photosynthetic primary production on Earth and play an essential role in supplying organic matter to marine food webs¹. They are key global actors in CO₂ uptake and provide gaseous oxygen to the atmosphere. A global decline of phytoplankton biomass has been reported over the past century (1% of chlorophyll-a concentration per year) leading to a decrease of net primary production in many oceanic regions². This decline is probably a consequence of global ocean warming which drives water column stratification, reducing the nutrient supply to surface waters. Temperature-driven reductions in phytoplankton productivity in tropical and temperate regions are likely to have cascading effects on higher trophic levels and ecosystem functioning³.

Photosynthetic picoeukaryotes (PPEs), defined by a cell diameter <3 µm, belong to different phyla, including Chlorophyta, Cryptophyta, Haptophyta, and Stramenopiles⁴. Present in all oceans, PPEs are the dominant primary producers in warm and oligotrophic regions⁵. Ocean warming and expansion of oligotrophic regions in the next decades may extend the ecological niche of PPEs, and a global shift from large photosynthetic organisms toward smaller primary producers is expected³,⁶. For example, sea ice melting in the Canadian Arctic Basin has been associated with an increase in the abundance of PPEs such as *Micromonas* at the expense of larger algae⁷. In the laboratory, this alga has the capacity to change its optimum temperature for growth in only a few hundred generations, suggesting that it will be less affected by global warming than many larger organisms⁸. In addition, the larger cell surface-volume ratio of PPEs compared to larger phytoplankton cells is advantageous for resource acquisition and growth in nutrient-limited environments⁹,¹⁰.

Iron is one key compound required for the activity of the respiratory chain, photosynthesis and nitrogen fixation¹⁰. Because bioavailable iron is extremely low in more than one-third of the surface ocean, small phytoplankton has developed several strategies to optimize iron uptake and reduce iron needs¹¹. In diatoms, reductive and non-reductive iron uptake mechanisms involve many proteins, including phytotransferrins, transmembrane ferric reductases, iron permeates, and siderophore-binding proteins¹². The iron needs can be modulated by the variation of gene expression levels between iron-required proteins and their iron-free equivalent. These protein switches include electron transfer (flavodoxin/ferredoxin), gluconeogenesis (fructose-bisphosphate aldolase type I or type II) and superoxide dismutases (Mn/Fe-SOD, Cu/Zn-SOD or Ni-SOD)¹³⁻¹⁵.

PPE growth is also limited by nitrogen (N) availability in large portions of the global ocean¹⁶. Ammonium (NH₄⁺), nitrate (NO₃⁻) and nitrite (NO₂⁻) are the primary source of inorganic N for PPEs, however, several studies have shown that dissolved organic N, like urea, can be metabolized in N-limited environments¹⁷. For example, several membrane-localized urea transporters in the diatom *Phaeodactylum tricornutum* are maximally expressed in nitrogen-limited conditions¹⁸ and the harmful algal blooms of the pelagophyceae *Aureococcus anophagefferens* may be fueled by urea¹⁹.

Despite their large taxonomic distribution, most molecular studies on the ecological role of PPEs and their adaptation to the environment are restricted to a few species. PPEs are suspected of possessing highly developed acclimation/adaptation capacities, but the underlying molecular mechanisms remain poorly characterized due to the lack of reference genomic data.

Among PPEs, *Pelagomonas calceolata* was the first described member of the Pelagophyceae class²⁰. It has since been identified in many oceanic regions using its 18 S rRNA sequence and chloroplast genome²¹⁻²³. Several studies have demonstrated the capacity of *P. calceolata* to adapt to different environmental conditions. In the laboratory, *P. calceolata* has been shown to exhibit a high degree of

acclimation to light fluctuations with rapid activation of the photo-protective xanthophyll cycle and non-photochemical quenching²⁴. In the Marquesas archipelago, *P. calceolata* is one of the most responsive species to iron fertilization with upregulation of genes involved in photosynthesis, amino acid synthesis and nitrogen assimilation¹³. A global-scale analysis of pelagophyte genes revealed that they are adapted to low-iron conditions¹⁴. In the subtropical Pacific, *P. calceolata* expresses stress genes in surface samples and genes involved in nitrogen assimilation are overexpressed in the deep-chlorophyll maximum²⁵. A laboratory study suggests that *P. calceolata* also has the ability to increase the transcription levels of organic-nitrogenous compound cleavage enzymes (cathepsin, urease, arginase) under low nitrogen concentration²⁶. Thus, gene expression appears to be controlled according to the nitrogen source and quantity. Taken together, this apparent adaptive plasticity may explain the presence of *P. calceolata* in many different oceanic environments, however, an exhaustive analysis of the genetic capacity of this species and the in situ characterization of its ecological niche is lacking.

Here we sequenced, assembled and annotated the *Pelagomonas calceolata* genome, with a combination of long- and short reads. We examined its genomic structure and gene content relative to other unicellular phytoplankton. We used this genome to detect *P. calceolata* in environmental datasets of *Tara* expeditions across all oceans to characterize its ecological niche and to identify the environmental conditions controlling its relative abundance. Finally, environmental expression levels of genes involved in nitrogen compounds and iron uptake and metabolism were studied.

## Results

**Chromosome-scale assembly and annotation of the *P. calceolata* genome.** To measure the abundance of *P. calceolata* in the oceans and study its genetic capacity to grow in different environmental conditions, we sequenced and assembled the genome of *P. calceolata* RCC100 using long reads of Oxford Nanopore Technologies (ONT) and Illumina short reads. Using the k-mer distribution of short reads, the genome was estimated to be homozygous with a size of 31 Mb (Supplementary Fig. S1a). The ONT long reads were assembled with Flye into six nuclear contigs for a total of 32.4 Mb, 1 plastid circular contig (90 Kb) and 1 mitochondrial circular contig (39 Kb) (Fig. 1, Supplementary Fig. S1b, c, and Supplementary Data 1). Two large and highly similar duplicated regions (>99% of identity) were detected at the extremity of contig 1 and 5 (393 Kb) and at the extremities of contig 3 and 6 (192 Kb; Supplementary Note 1 and Supplementary Fig. S1d). (TTAGGG)n telomeric repeats were detected at both ends of contigs 2, 3, 4, and 6, indicating that these four contigs represent complete chromosomes (Supplementary Fig. S1d). For contig 1 and contig 5, telomeric sequences were identified at only one extremity, the other extremity ending in the duplicated region. We used the Hi-C long-range technology to validate the assembly of *P. calceolata* genome sequence. The interaction map revealed a high number of contacts within contigs and very few across contigs (Supplementary Fig. S2 and Supplementary Note 2). No chimeras or fragmentations were detected. This result confirms that the six contigs correspond to six chromosomes of *P. calceolata*.

A total of 16,667 genes were predicted on the *P. calceolata* genome (see “Methods”), which is a high number for a PPE (Table 1 and Supplementary Data 2). There was an average of 0.45 intron per gene, and the distribution of their lengths reveals a peak at around 210 bp, which is the characteristic length of Introner Elements described in *A. anophagefferens*²⁷ (Supplementary Fig. S3, Supplementary Data 3, and Supplementary Note 3). In all, 9812 *P. calceolata* predicted proteins (58%) are homologous with at least
<figure><img src="image_3.png" /><figcaption>**Fig. 1** *Pelagomonas calceolata* nuclear genome. Representation of the 6 nuclear contigs of *P. calceolata*. The top layer indicates the number of genes per 100 Kb (black bars), the middle layer represents the GC content in percentage over a window of 200 Kb and the bottom layer is the position of DNA repeats of more than 500 bases repeated at least five times over the entire genome. Red, orange, and yellow bars indicate three different repeats in low-GC regions present in at least three different contigs. Dashed red and blue rectangles are duplicated chromosomal regions.</figcaption></figure>

<table><caption>**Table 1** Genome characteristics of several unicellular photosynthetic eukaryotes.</caption><thead><tr><th>Phylum/class</th><th>Species</th><th>Genome size (Mb)</th><th>Number of chromosomes</th><th>Predicted genes</th><th>GC%</th><th>Cell size</th><th>References</th></tr></thead><tbody><tr><td>Pelagophyceae</td><td>Pelagomonas calceolata</td><td>32.4</td><td>6</td><td>16,667</td><td>63.6</td><td>2 µm</td><td>This study</td></tr><tr><td>Pelagophyceae</td><td>Aureococcus<br>anophagefferens</td><td>56.0</td><td>Unknown</td><td>11,520</td><td>67.4</td><td>2 µm</td><td>89</td></tr><tr><td>Eustigmatophyceae</td><td>Nannochloropsis oceanica</td><td>29.3</td><td>32</td><td>7730</td><td>54.0</td><td>3 µm</td><td>90</td></tr><tr><td>Diatom</td><td>Phaeoactylum tricornutum</td><td>27.0</td><td>33</td><td>10,402</td><td>48.8</td><td>11 µm</td><td>91</td></tr><tr><td>Diatom</td><td>Thalassiosira pseudonana</td><td>34.5</td><td>24</td><td>11,776</td><td>46.9</td><td>5 µm</td><td>92</td></tr><tr><td>Chlorophyta</td><td>Micromonas pusilla</td><td>21.9</td><td>17</td><td>10,575</td><td>65</td><td>≤2 µm</td><td>93</td></tr><tr><td>Chlorophyta</td><td>Ostreococcus lucimarinus</td><td>13.2</td><td>21</td><td>7651</td><td>60</td><td>1.3 µm</td><td>94</td></tr><tr><td>Chlorophyta</td><td>Bathyococcus prasinos</td><td>15.1</td><td>19</td><td>7847</td><td>48</td><td>1-2 µm</td><td>95</td></tr><tr><td>Haptophyta</td><td>Emiliania huxleyi</td><td>141.7</td><td>Unknown</td><td>38,549</td><td>64.5</td><td>4-5 µm</td><td>96</td></tr></tbody></table>

one gene in a stramenopile genome, including 2631 (16%) only shared with the pelagophyte *A. anophagefferens* (Supplementary Fig. S4). A conserved functional domain (Pfam, KO or InterProScan) was found in 11,240 proteins (67%). Even if gene completeness estimations are imprecise for species distant from model organisms, we obtained 94.0% of completeness with BUSCO<sup>28</sup> (88% single-copy and 6% duplicated genes), showing that our genome is near complete (Supplementary Data 4).

**GC content, centromeres, and meiosis in *P. calceolata*.** A remarkable feature in the *P. calceolata* genome is the distribution of GC content along *P. calceolata* chromosomes (Fig. 1). While the average GC content of the nuclear genome is 63%, one large region

in each contig (259 Kb on average) is 52% GC. These unique large troughs in GC content in each chromosome suggest that these regions encompass centromeres. The Hi-C result confirms centromere positions with the presence of contacts between low-GC regions across chromosomes, suggesting the physical proximity of these regions in the nucleus (Supplementary Fig. S2). Interestingly, we did not observe an accumulation of repeated elements or transposons in these low-GC regions and only a slight decrease of gene density. Genomic specificities and gene content of low-GC regions are detailed in Supplementary Data 5, S6 and Supplementary Note 4. Low-GC patterns in centromeres could be explained by the inhibition of recombination<sup>29,30</sup>, suggesting that *P. calceolata* is capable of meiosis and recombination. Among 23 meiosis-specific
genes characterized in other species³¹⁻³³, 18 homologs are present in the *P. calceolata* genome (Supplementary Data 7). These genes include the double-strand DNA break (DSB) initiator SPO11; RAD50, RAD52, and MRE11 to bind DSBs; HOP2, MND1, DMC1, and RAD51 to ensure pairing and invade the homologous strand; MSH2, MSH3 PMS1, and MSH6 genes involved in the synthesis-dependent strand annealing pathway and MUS81 necessary for non-interfering (class II) crossing over. The five missing genes are absent in many organisms capable of meiosis, suggesting that they are not required to perform genetic recombination (Supplementary Note 5). Taken together, the genomic structure and the genetic content of *P. calceolata* strongly suggests that this species performs meiosis.

**Relative abundance of *P. calceolata* across oceanic basins.** To estimate the relative abundance of *P. calceolata* across all oceans, we first used the abundance of the V9 region of the 18 S rRNA sequenced from all samples of the *Tara* Oceans expedition³⁴. The most abundant *P. calceolata* Operational Taxonomic Unit (OTU) in the 0.8–5 µm size fraction is on average 0.80% for the 104 surface samples and 1.23% for the 61 deep-chlorophyll maximum (DCM) samples (Supplementary Data 8). According to this abundance estimation method, *P. calceolata* is the third most abundant eukaryote OTU of the 0.8–5 µm size fraction after two Dinophyceae OTUs affiliated to *Ankistrodinium* and an unknown Gymnodiniaceae.

To estimate *P. calceolata* abundance independently from PCRs and 18 S rRNA copy number bias, we used the mapping of metagenomic reads on the *P. calceolata* genome. For the 0.8–5 µm size fraction, the percentage of sequenced reads aligned on the genome is 1.39% (n = 93, sd = 1.5) in surface samples and 2.67% (n = 55, sd = 1.6) in DCM samples. In the 0.8–2000 µm size fraction, *P. calceolata* represents 1.01% (n = 80, sd = 1.2) of all reads in surface samples and 1.56% (n = 39, sd = 1.3) in DCM samples. A maximal relative abundance of 6.7% in the 0.8–5 µm size fraction was observed in the North Indian Ocean (station TARA_38) at the DCM (Fig. 2a). In the Indian Ocean, Red Sea and Mediterranean Sea, *P. calceolata* is significantly more abundant in the DCM than at the surface (Fig. 2b). In cold waters (below 10 °C), *P. calceolata* is not detected above our threshold of 25% of horizontal genomic coverage. Important variations between and within each oceanic basin are observed, suggesting that many biotic or abiotic factors influence *P. calceolata* abundance.

Finally, we compared the two methods of abundance estimations (Supplementary Fig. S5). The metagenomic-based relative abundance is strongly correlated to the metabarcoding-based relative abundance (Pearson correlations of 0.91 and 0.70 for the 0.8–2000 and 0.8–5 µm size fractions, respectively). However, the metabarcoding-based abundance is on average 2.3 lower in the 0.8–5 µm size fraction and 3.1 lower in the 0.8–2000 size fraction compared to the metagenomic-based abundance (Supplementary Fig. S5).

**High relative abundance of *P. calceolata* in temperate, low-light, and iron-poor regions.** In order to identify factors controlling *P. calceolata* abundance in the oceans, we used physical-chemical parameters available for each oceanic station (see “Methods”). Principal component analysis revealed a positive relation between metagenomic-based *P. calceolata* abundance, the temperature, and the coast distance and a negative relation with iron concentration (Fig. 3a, b). This result was consistent over the 2 size fractions containing *P. calceolata* cells (0.8–5 µm and 0.8–2000 µm size fractions; Supplementary Fig. S6). Despite the numerous factors potentially influencing *P. calceolata* abundance, we observed a weak but significant Pearson's positive correlation with the

temperature, a negative correlation with Photosynthetically Active Radiation (PAR, mean of 30 days) and a negative correlation with iron concentrations (Table 2). In the 0.8–5 µm size fraction, the relative abundance of *P. calceolata* is higher in low-iron conditions (<0.2 nmol/l, 54 samples) with an average 2.3% of metagenomic reads than in high-iron environments (>0.2 nmol/l, 88 samples) with an average 1.7% of metagenomic reads (Wilcoxon test, *P* value = 0.02). In the 0.8–2000 µm size fraction, we observe the same tendency with a relative abundance of 1.9% of metagenomic reads on average in low-iron waters (49 samples) and a lower relative abundance of 0.78% of metagenomic reads on average in high-iron environments (59 samples) (Wilcoxon test, *P* value = 9.6e⁻⁷). In addition, *P. calceolata* relative abundance is weakly correlated with the 9'butanoyloxyfucoxanthin concentration, a signature pigment for pelagophytes (Pearson 0.22, *P* value = 0.02 and Pearson 0.41, *P* value = 4.82e⁻⁰⁵ in the 0.8–5 µm and 0.8–2000 µm size fraction, respectively). We used a general additive model to estimate the contribution of temperature, PAR and iron concentration to *P. calceolata* relative abundance (Table 2). The three factors explain 32.3% of the variations of *P. calceolata* abundance in the 0.8–5 µm size fraction and 56.8% in the 0.8–2000 µm size fraction.

Finally, we projected the relative abundance of *P. calceolata* at the end of the century following Frémont et al. methodology³⁵. We modeled the ecological niche of *P. calceolata* using the World Ocean Atlas (WOA18) datasets at the time and location of sampling or using the projected climatology in 2099 using the RPC8.5 scenario (see “Methods”). We used four machine-learning techniques: Generalized Additive Models (GAM), Neural Networks (nn), Random Forest (rf) and Gradient Boosted Trees (bt) and evaluated their performances with two parameters. The Pearson correlation coefficient indicates the correlation between the model and in situ measurements of *P. calceolata* abundance. The four machine-learning tools have similar performances based on Pearson's correlations (nn = 0.676; gam = 0.621; bt = 0.683; rf = 0.694). The second parameter is the root mean square error (rmse) and reflects the magnitude of the errors in the models (the number of standard deviations from the mean). Using this metric the GAM approach is less good (rmse = 1.04) than the three other tools (nn = 0.964; bt = 0.952; rf = 0.941). These results indicate that we have enough in situ data to capture the global trends on the relative abundance of *P. calceolata* but these models could be imprecise on the amplitude of abundance variations. In addition, the predictions in the tropical waters are uncertain because this environment in 2099 is out of the range of the training dataset. Because the performances of the four models are similar, we combined them to obtain the most accurate projection (Fig. 3c and Supplementary Fig. S7). Despite these limitations, we projected an increase of up to 1.12% of *P. calceolata* relative abundance from latitude 40° to latitude 50° in the North and South hemispheres and a decrease in temperate and tropical waters (−0.8% maximum).

**Genes related to iron uptake, storage and usage in *P. calceolata*.** Iron is a critical metal for all photosynthetic organisms, required for the photosynthesis, the nitrogen cycle and the protection against reactive oxygen species. Since *P. calceolata* seems to thrive in iron-poor waters, we identified *P. calceolata* genes coding for iron uptake and storage, then compared their expression levels in low (<0.2 nmol/l) versus high (>0.2 nmol/l) iron conditions using *Tara* Oceans metatranscriptomes. *P. calceolata* has five genes encoding the phytotransferrin *ISIP2A* involved in Fe³⁺ uptake via endosomal vesicles and 2 putative iron-storage protein *ISIP3*. In iron-poor environments, three *ISIP2A* and one *ISIP3* are over-expressed (Fig. 4a and Supplementary Data 9). This result indicates
<figure><img src="image_4.png" /><figcaption>**Fig. 2** Relative abundance and distribution of *P. calceolata* in the oceans. **a** World map of the relative abundance of *P. calceolata* metagenomic reads. The color code indicates the percentage of sequenced reads aligned on the genome. The DCM samples of size fractions 0.8–5 µm (circles) or 0.8–2000 µm (triangles) are shown. *P. calceolata* is considered to be absent when the horizontal coverage is below 25% of the genome (gray dots). **b** Boxplot of the relative abundance in each oceanic region in surface and DCM samples. Red stars indicate a significant difference between SUR and DCM samples (Wilcoxon test, *P* value <0.01).</figcaption></figure>

that similarly to diatoms, ISIP are upregulated following iron
starvation in *P. calceolata* potentially improving cell growth in
low-iron environments. Three genes encode the iron transporter
ferroportin but are not differentially expressed according to
the environment (Supplementary Fig. S8a). These proteins are
iron exporters in multicellular organisms but their function in
microalgae remains to be studied³⁶. Finally, we identified eight
Zinc/iron permeases potentially involved in iron uptake from the

environment in the *P. calceolata* genome. Among them, two are
overexpressed in high-iron and one in low-iron environments.
Interestingly, we note the absence of the iron permease *FTR1*,
the iron-storage ferritin and the starvation induced protein *ISIP1*
 invольved in endocytosis of siderophores in diatoms). In compar-
ison to *P. calceolata*, the coastal Pelagophyceae *A. anophagefferens*
do not have *ISIP3* gene and a lower number of Zinc-Iron
permeases.
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /><figcaption>**Fig. 3** Ecological niche of *P. calceolata*. **a** Principal component analysis of the metagenomic-based relative abundance of *P. calceolata* in the 0.8–5 µm size fraction. Percentages of variance explained by each axis are indicated on axis titles. Top panel: each dot represents a sample with a size proportional to the relative abundance of *P. calceolata* and the colors indicate the oceanic basins. Bottom panel: nine environmental parameters are represented as vectors alongside the relative abundance of *P. calceolata* (blue vector). **b** Bubble plot of the relative abundance of *P. calceolata* for the 0.8–5 µm size fraction according to the nine environmental parameters. **c** Delta of the modeled relative abundance of *P. calceolata* between 2010 and 2099. Green areas correspond to a decrease while purple areas correspond to an increase of *P. calceolata* relative abundance. Small stars indicate locations where at least one of the predictor drivers is out of range compared to the training dataset values.</figcaption></figure>
<table><caption>Table 2 Environmental parameters explaining *P. calceolata* relative abundance for the 0.8-5 µm (a) and the 0.8-2000 µm (b) size fractions.</caption><thead><tr><th rowspan="2">(a)<br>0.8-5 µm</th><th colspan="3">GAM model</th><th colspan="2">GAM verification</th><th colspan="2">Pearson correlations</th></tr><tr><th>edf</th><th>F value</th><th>P value</th><th>k-index</th><th>k P value</th><th>r</th><th>P value</th></tr></thead><tbody><tr><td>s(Temperature)</td><td>1</td><td>22.16</td><td>6.84e-6</td><td>1.11</td><td>0.87</td><td>0.23</td><td>0.001</td></tr><tr><td>s(Iron concentration)</td><td>1.257</td><td>13.12</td><td>2.26e-4</td><td>0.93</td><td>0.12</td><td>-0.25</td><td>0.001</td></tr><tr><td>s(PAR 30 days)</td><td>1.859</td><td>15.94</td><td>4.56e-7</td><td>1.01</td><td>0.44</td><td>-0.32</td><td>0.001</td></tr><tr><td>Adjusted R<sup>2</sup></td><td colspan="3">0.3</td><td colspan="3"></td></tr><tr><td>Deviance explained</td><td colspan="3">32.30%</td><td colspan="3"></td></tr></tbody></table><table><thead><tr><th rowspan="2">(b)<br>0.8-2000 µm</th><th colspan="3">GAM model</th><th colspan="2">GAM verification</th><th colspan="2">Pearson correlations</th></tr><tr><th>edf</th><th>F value</th><th>P value</th><th>k-index</th><th>k P value</th><th>r</th><th>P value</th></tr></thead><tbody><tr><td>s(Temperature)</td><td>3.628</td><td>9.442</td><td>1.71e-6</td><td>0.96</td><td>0.31</td><td>0.57</td><td>0.0001</td></tr><tr><td>s(Iron concentration)</td><td>1.54</td><td>1.225</td><td>2.42e-1</td><td>0.91</td><td>0.14</td><td>-0.47</td><td>0.0001</td></tr><tr><td>s(PAR 30 days)</td><td>2.454</td><td>6.962</td><td>0.00027</td><td>1.14</td><td>0.93</td><td>-0.051</td><td>0.6</td></tr><tr><td>Adjusted R<sup>2</sup></td><td colspan="3">0.53</td><td colspan="3"></td></tr><tr><td>Deviance explained</td><td colspan="3">56.80%</td><td colspan="3"></td></tr></tbody></table>

Several important ferrous proteins can be substituted by non-ferrous equivalents in iron-poor environments¹⁴. In the *P. calceolata* genome we identified 11 flavodoxin genes involved in electron transfer during photosynthesis, potentially replacing ferredoxins (17 genes; Supplementary Data 9). This number of genes is important compared to other algae, including *A. anophagefferens* (5 flavodoxins and 9 ferredoxins). Expression levels of these genes across the oceans revealed overexpression of flavodoxin genes in low-iron environments, replaced by ferredoxins in high-iron conditions (Fig. 4b). The Fructose-Bisphosphate Aldolase (FBA), necessary for gluconeogenesis and the Calvin cycle, is encoded by six genes in *P. calceolata*. Two genes are dependent on a divalent cation (FBA type II), and the four others are Zinc/Iron-independent (FBA type I). FBA type I is overexpressed in high-iron conditions and quasi-absent in low-iron conditions (Fig. 4c). Finally, all types of Superoxide dismutases (SOD) are found in the *P. calceolata* genome. Non-ferrous SODs (Cu/Zn and Ni) encoded by three genes and Mn/ Fe-SOD encoded by two genes are not differentially expressed according to iron concentrations (Supplementary Fig. S8b). The gene content and the transcriptomic flexibility suggest important capacities for *P. calceolata* growth in low-iron environments.

**Genes involved in nitrogen uptake, storage, and recycling.** Because *P. calceolata* could be an important player in the nitrogen (N) cycle in oceanic ecosystems²⁵, we explored its gene content and analyzed the expression levels of genes involved in nitrogen metabolism in the environment using *Tara* Oceans metatranscriptomes (Fig. 5 and Supplementary Data 9). The uptake of nitrogen-containing inorganic compounds is supported by 11 genes in *P. calceolata*. Three genes encode nitrate/nitrite transporters, three genes encode formate/nitrite transporters and five genes encode ammonium transporters. One formate/nitrite and two nitrate/nitrite transporter genes are significantly overexpressed in low-nitrate conditions (Supplementary Fig. S9). One ammonium transporter is overexpressed in low-ammonium environments (Supplementary Fig. S9). We identified one nitrate reductase and two nitrite reductases in the *P. calceolata* genome (Supplementary Data 9). The expression level is significantly higher in high-nitrate environments for the nitrate reductase and one nitrite reductase (Supplementary Fig. S10). But the second nitrite reductase is surprisingly less expressed in high-nitrate conditions. The number of enzymes incorporating ammonium into organic compounds (GS/GOGAT pathway) is higher in *P. calceolata* than in other species: five glutamine

synthetase (GS) and four glutamate synthase (GOGAT) genes are present in the *P. calceolata* genome. Two GS genes are more expressed in high-nitrate samples and 1 GOGAT gene is more expressed in low-nitrate samples (Supplementary Fig. S10).

We identified three genes carrying the nitrate and nitrite sensing (NIT) domain (IPR013587) in the *P. calceolata* genome. Using NCBI non-redundant proteins and marine genomic databases (see “Methods”), 60 homologous proteins of the NIT-sensing domains of *P. calceolata* were identified. These homologs are restricted to the Pelagophyceae class (16 transcriptomes), the Dictyochophyceae class (6 transcriptomes) and one putative cryptophyte transcriptome. The phylogenetic tree of this protein family shows three subfamilies diverging before the Dictyochophyceae/Pelagophyceae separation (Supplementary Fig. S11). One *P. calceolata* protein carries a NIT-sensing domain surrounded by two transmembrane domains suggesting a capacity for external nitrate/nitrite sensing while the two other NIT genes carry a protein-kinase domain (IPR000719) suggesting phosphorylation-based signal transduction dependent on intracellular nitrate or nitrite concentration (Fig. 5b). To investigate this possibility, we studied the expression levels of NIT-sensing genes. The potential external NIT-sensing gene is indeed significantly overexpressed in low-nitrate environments. In contrast, only one of the intracellular NIT-sensing genes is upregulated in nitrate-rich environments (Fig. 5c). This result suggests an important role of the NIT-sensing genes in the acclimation to environmental nitrate concentrations.

Finally, we identified genes involved in nitrogen recycling from organic compounds which are important in several species in case of inorganic nitrogen deprivation. One arginase and one cyanate lyase were detected in the *P. calceolata* genome but no gene encoding formamidase. In addition, the number of gene copies for enzymes involved in the urea cycle (carbamoyl-phosphate synthetase, ornithine carbamoyltransferase, argininosuccinate synthase and argininosuccinate lyase) is equal or slightly lower than in other algae (Fig. 5a and Supplementary Data 9). Among these genes, only the cyanate lyase is overexpressed in low-nitrate conditions suggesting cyanate is an important alternative source of organic nitrogen for *P. calceolata* (Supplementary Fig. S10).

Discussion

The essential roles of phytoplankton in oceanic ecosystems have been illustrated many times, however, numerous lineages are still poorly explored and model organisms are restricted to a few taxa (mainly diatoms, prasinophytes, and haptophytes) limiting the global understanding of phytoplankton activity. The *P. calceolata*
<figure><img src="image_9.png" /><figcaption>**Fig. 4** Expression of iron-related genes in *P. calceolata*. **a** Relative gene expression levels normalized in transcript per million (TPM) of five phytotransferrins (ISIP2A) and two putative iron storage (ISIP3) in low-iron (<0.2 nM) and high-iron (>0.2 nM) oceanic stations. *P* values of Wilcoxon statistical tests between low- and high-iron conditions are indicated for each gene. Significant *P* values (<0.01) are in bold. **b** Relative expression levels (TPM) of genes coding for ferredoxin (orange) and its non-ferrous equivalent flavodoxin (purple) in each *Tara Oceans* sample. Samples are sorted from low-iron (left) to high-iron (right) conditions. Iron concentrations are indicated in nM on the colored horizontal bar. **c** Same representation for genes coding for fructose-bisphosphate aldolase II (orange) and its non-ferrous equivalent fructose-bisphosphate aldolase I (purple).</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Fig. 5** Nitrogen sensing and metabolism in *P. calceolata*. **a** Schematic representation of N transport and assimilation in *P. calceolata* based on the gene content. The color code indicates if the number of gene copies for a specific function is overrepresented (green), equally represented (blue), underrepresented (orange) or absent (red) in *P. calceolata* genome compared to the mean of eight pico-nano photosynthetic eukaryotes. Gene copy number for each function is indicated in Supplementary Data 9. **b** Domain organization of NIT-sensing proteins in *P. calceolata*. Orange boxes are NIT-sensing domains (IPR13587), blue boxes are serine-threonine/tyrosine-kinase domain (IPR20635), and yellow rectangles are transmembrane domains. **c** Relative expression levels (TPM) of three NIT-sensing genes in low-nitrate (<2 µM) and high-nitrate (>2 µM) environments. P values of Mann-Whitney-Wilcoxon tests between low- and high-nitrate samples are indicated for each gene. Significant P values (<0.01) are in bold.</figcaption></figure>

genome assembled and annotated in this study reveals a pre-
viously underestimated high abundance of *P. calceolata* in the
oceans and brings new insights into specific genomic features of
this algae class related to its adaptation to specific environments.

We have shown in this study that *P. calceolata* is cosmopolitan
in oceanic samples above 10 °C with a relative abundance

generally >1% of all sequenced reads. In contrast to the coastal
Pelagophyceae *A. anophagefferens* that can present high peaks of
abundance³⁷, no *P. calceolata* blooms were reported, but *P. cal-
ceolata* is well-adapted to an extensive range of environmental
conditions as suggested by previous studies²¹,²³. Although the
abundance of an organism calculated from metabarcoding or
metagenomic data provides only an indirect and relative quantification of organism abundances, both methods suggest that *P. calceolata* is one of the most abundant pico-nano eukaryote in offshore data. The high relative abundance of *P. calceolata* measured with a metabarcoding approach has recently been confirmed with a qPCR method (average of 5 882 ± 2 855 rRNA gene copies mL⁻¹ on the surface of the eastern North Pacific)²¹. In addition, we have shown that the metabarcoding approach probably underestimates the relative abundance of *P. calceolata* compared to the metagenomic analysis owing to the low copy number of rRNAs in this organism. However, we cannot exclude that the large genome size of *P. calceolata* compared to bacterial genomes present in the 0.8–5 µm size fraction overestimates its relative abundance in metagenomic datasets. Further studies may combine microscopic and flow sorting approaches with genomic data to assess the number of cells and the biomass of this organism in the oceans. Our model analysis has revealed a probable increase of *P. calceolata* relative abundance at the end of the century in high latitudes where the seawater temperature is currently too low for this species. This result is coherent with previous studies suggesting a global increase of phytoplankton in subpolar regions³⁸,³⁹.

Iron is essential for growth, photosynthesis, primary production, nitrogen fixation and reduction for PPEs⁴⁰. Our results show that *P. calceolata* thrives in iron-poor waters and thus occupies a large ecological niche for a PPE. Two main strategies exist against iron deprivation: the optimization of iron uptake and the modulation of iron needs. Genes coding for iron chelators and ferritin are absent from the *P. calceolata* genome, and genes coding for passive iron transporters are under- or equally represented compared to other PPEs. In contrast, phytotransferrins (ISIP2A), putative iron-storage proteins (ISIP3) and ferroportins are overrepresented in the *P. calceolata* genome. The expression levels of ISIP genes are anti-correlated with iron concentration in *P. calceolata* showing an acclimation to low-iron conditions. Because phytotransferrins are dependent on carbonate ions, ocean acidification may reduce iron uptake efficiency in many species including *P. calceolata*⁴¹,⁴². Compared to *P. calceolata*, the low abundance of *A. anophagefferens* in open oceans where iron is limited is consistent with the absence of genes involved in iron uptake and storage, including iron permeases and ISIP3 genes. The presence of three ferroportin genes in *P. calceolata* is interesting since these transmembrane iron export proteins play a significant role in iron homeostasis in multicellular organisms³⁶. Ferroportin function in microalgae is unknown but could act to export iron from endosomes to the cytoplasm⁴³. In the green alga *Chlamydomonas reinhardtii*, a ferroportin gene is overexpressed under low Fe conditions⁴⁴. Even though we did not find significant changes in the expression of the 3 ferroportin genes of *P. calceolata* according to iron concentration, the function of these genes could be investigated to understand their role in variable iron concentrations.

Modulation of iron needs seems to be a major acclimation strategy for *P. calceolata* in low-iron environments. All known molecular switches between ferrous and non-ferrous proteins are present in the *P. calceolata* genome and the transcriptomic regulation of these genes suggests a central role of these proteins for its growth in low-iron environments.

Expressing more than 90% of all nitrate transporter transcripts, pelagophytes may dominate nitrate uptake and assimilation in the North Pacific Ocean²⁵. Indeed, *P. calceolata* contains a large collection of genes for inorganic nitrogen transporters. The main difference with the coastal *A. anophagefferens* is the reduced number of ammonium transporters in *P. calceolata*.

Organic nitrogen compounds could also be a major nitrogen source of for *P. calceolata*. We have shown that the cyanate lyase

and urease genes are expressed in many environments but only the cyanate lyase is overexpressed in low-nitrate conditions. These two genes, largely present among phytoplankton lineages, could be significant components of acclimation to low-nitrate environments⁴⁵. In addition, *A. anophagefferens* grow faster on the organic nitrogen substrates (urea and glutamic acid) than on nitrate or ammonium, suggesting that the dominance of Pelagophyceae in low-nitrate environments could be due to an optimized usage of these organic molecules⁴⁶.

One remarkable feature of the *P. calceolata* genome is the presence of three genes carrying NIT domains (PF08376). This NIT domain was first described in bacterial nitrite and nitrate sensor proteins⁴⁷. This sensor is an alpha-helical protein playing a signal transduction role in regulating gene expression, cell motility and enzyme activity in *Klebsiella oxytoca*⁴⁸. In pico-nano algae, NIT-sensing domains can be associated with a serine-threonine/tyrosine-kinase domain, suggesting signal transduction according to intracellular nitrate/nitrite concentration, or surrounded by two transmembrane domains suggesting extracellular sensing. Even though NIT-sensing domains can be found across various phyla, homologs of *P. calceolata* NIT proteins are restricted to pelagophytes and dictyochophytes. The NIT genes in *P. calceolata* are highly expressed in subtropical Pacific N-depleted waters, suggesting that these proteins have a role in transcription regulation according to nitrate availability²⁵. Our results suggest that NIT-sensing proteins respond differently to environmental nitrate depletion. We can hypothesize that the NIT-sensing protein overexpressed in nitrate-rich environments plays a role in the intracellular regulation of stored nitrogen, activating pathways when nitrate or nitrite stocks are sufficient. In contrast, putative extracellular NIT-sensing could be an environmental nitrate or nitrite sensor activated to regulate the expression of genes involved in the acclimation to low-nitrate conditions.

In summary, due to its widespread distribution and its high abundance in the open oceans, *Pelagomonas calceolata* can serve as an ecologically-relevant model to study marine photosynthetic protists. We used the chromosome-scale genome sequence, mostly telomere-to-telomere, generated in this study to estimate its abundance in environmental datasets. We have shown that the *P. calceolata* genome has specific genomic features potentially explaining its ecological success in open oceans. Compared to the coastal Pelagophyceae *A. anophagefferens*, the *P. calceolata* genome contains specific genes involved in the acclimation to low-iron conditions. The large repertoire of genes involved in nutrient acquisition from the environment is coherent with its widespread pattern of relative abundance distribution across different environments. The ecological niche of *P. calceolata* suggests that this alga will benefit from the global climate change with the extension of oligotrophic regions and global ocean warming. Future studies could use the *P. calceolata* genome to explore adaptation and acclimation processes controlling the distribution and abundance of this alga.

## Methods

**Pelagomonas culture.** *Pelagomonas calceolata* RCC100 culture was grown in 12:12-h light:dark photoperiod in K medium with natural seawater base at 20 °C. At the Roscoff Culture Collection, cells were kept at a light intensity of ~80 µmol photon m⁻² s⁻¹ and the volume of culture was ramped up to 1 litre in mid-exponential growth phase before harvesting. RCC100 culture was not axenic and grown in the presence of undefined bacterial microbiota.

**DNA extraction, library preparation, and sequencing.** We pelleted cells from 500 ml of culture by two successive centrifugations at 10,000 × g for 15 min at 4 °C. Genomic DNA was extracted using the NucleoSpin Plant II Mini kit according to the manufacturer's instructions (Macherey-Nagel, Hoerdt, France) with the following exception for the lysis step: 400 µL of lysis buffer PL1 and 25 µL of proteinase K 25 mg/mL were added to strain pellets, and lysates were incubated at 55 °C for 1 h at 900 rpm. DNA quantity and integrity were respectively evaluated
on a Qubit 2.0 spectrofluorometer (Invitrogen, Carlsbad, CA, USA) and a Nanodrop1000 spectrophotometer (Thermo Fisher Scientific, MA, USA).

For Illumina sequencing, DNA (1.5 µg) was sonicated using a Covaris E220 sonicator (Covaris, Woburn, MA, USA). Fragments were end-repaired, 3'-adenylated and Illumina adapters (Bioo Scientific, Austin, TX, USA) added using the Kapa Hyper Prep Kit (KapaBiosystems, Wilmington, MA, USA). Ligation products were purified with AMPure XP beads (Beckmann Coulter Genomics, Danvers, MA, USA). The library was then quantified by qPCR using the KAPA Library Quantification Kit for Illumina Libraries (KapaBiosystems), and the library profile was assessed using a High Sensitivity DNA kit on an Agilent Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). The library was sequenced on an Illumina NovaSeq instrument (Illumina, San Diego, CA, USA) using 150 base-length read chemistry in paired-end mode.

For ONT sequencing, the library was prepared using the 1D Native barcoding genomic DNA (with EXP-NBD104 and SQK-LSK109). Genomic DNA fragments (1 µg) were repaired and 3'-adenylated with the NEBNext FFPE DNA Repair Mix, and the NEBNext® Ultra™ II End Repair/dA-Tailing Module (New England Biolabs, Ipswich, MA, USA). Adapters with barcodes provided by ONT were then ligated using the NEB Blunt/T'A Ligase Master Mix (NEB). After purification with AMPure XP beads (Beckmann Coulter, Brea, CA, USA), the sequencing adapters (ONT) were added using the NEBNext Quick T4 DNA ligase (NEB). The library was purified with AMPure XP beads (Beckmann Coulter), then mixed with the Sequencing Buffer (ONT) and the Loading Band (ONT), and loaded on a MinION R9.4.1 flow cell. Reads were basecalled using Guppy 3.1.5.

For the Hi-C sequencing, the library was prepared using the Dovetail Omni-C kit (Dovetail Genomics, Scotts Valley, CA, USA). A *P. calceolata* RCC100 culture (60 mL corresponding approximately to $6 \times 10^7$ cells) was first centrifuged at 5000 × g for 10 min. The pellet was processed as mammalian cells, following the Mammalian Cell Protocol for Sample Preparation (version 1.4) without using DSG cross-linking reagent. Briefly, the chromatin was fixed with formaldehyde, randomly digested with DNase I and then extracted. Chromatin ends were repaired and ligated to a biotinylated bridge adapter, followed by proximity ligation of adapter-containing ends. After proximity ligation, crosslinks were reversed and DNA was purified. Purified DNA was treated to remove biotin that was not internal to ligated fragments, and a sequencing library was generated using NEBNext Ultra enzymes and Illumina-compatible adapters. Biotin-containing fragments were isolated using streptavidin beads before PCR enrichment of the library. The Dovetail Hi-C library quality was checked as described above and sequenced on an Illumina MiSeq instrument (Illumina, San Diego, CA, USA) in paired mode (2\*150 bp), producing 2,832,092 reads. The Hi-C raw reads were aligned against the assembly (-s none option) using Juicer (Juicer version 1.5.6 - Juicer Tools Version 1.9.9)$^{49}$. The contact map representation was generated with R version 4.1.1 using the merged nodups file.

**RNA extraction, library preparation, and sequencing.** When the cell concentration reached 10 million cells/mL in the mid-exponential growth phase, 160 mL of culture were collected by three successive filtrations on 1.2-µm poly-carbonate filters of 47 mm to avoid prokaryotic contamination. To preserve cell and RNA integrity, we kept filtration time and pressure below 10 min and 20 mmHg, respectively. Then filters were stored in 15-mL Falcon tubes with 3 mL of Trizol (Invitrogen, Carlsbad, CA, USA), mixed and flash-frozen in liquid nitrogen for further processing. RNA was extracted by incubation at 65 °C for 15 min, followed by chloroform extraction. The aqueous phase was purified using a Purelink RNA Isolation kit (Ambion Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instructions. DNA contamination was removed by digestion using the TURBO DNA-free™ Kit (Ambion Invitrogen) according to the manufacturer's DNase treatment protocol. After two rounds of 30-min incubation at 37 °C, the efficiency of DNase treatment was assessed by PCR. Quantity and quality of extracted RNA were analyzed with RNA-specific fluorimetric quantitation on a Qubit 2.0 Fluorometer using Qubit RNA HS Assay (Invitrogen). The qualities of total RNA were checked by capillary electrophoresis on an Agilent Bioanalyzer, using the RNA 6000 Nano LabChip kit (Agilent Technologies, Santa Clara, CA).

RNA-Seq library preparation was carried out from 1 µg total RNA using the TruSeq Stranded mRNA kit (Illumina, San Diego, CA, USA), allowing mRNA strand orientation. Briefly, poly(A) + RNAs were selected with oligo(dT) beads, chemically fragmented and converted into single-stranded cDNA using random hexamer priming. After second strand synthesis, double-stranded cDNA was 3'-adenylated and ligated to Illumina adapters. Ligation products were PCR-amplified following the manufacturer's recommendations. Finally, the ready-to-sequence Illumina library was quantified by qPCR using the KAPA Library Quantification Kit for Illumina libraries (KapaBiosystems, Wilmington, MA, USA), and evaluated with an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). The library was sequenced using 101 bp paired-end read chemistry on a HiSeq2000 Illumina sequencer. Low-quality nucleotides (Q < 20) from both ends of the reads were discarded. Illumina sequencing adapters and primer sequences were removed and reads shorter than 30 nucleotides after trimming were discarded. These trimming and cleaning steps were achieved using in-house-designed software based on the FastX package (https://www.genoscope.cns.fr/externe/fastxtend/). The last step identifies and discards read pairs mapped to the phage phiX genome, using the SOAP aligner$^{50}$ and the Enterobacteria phage PhiX174 reference sequence

(GenBank: NC_001422.1). This processing, described in ref. 51, resulted in high-quality data. Moreover, ribosomal RNA-like reads were excluded using SortMeRNA$^{52}$ 2.1 and SILVA databases.

**Long-read-based genome assembly.** Raw nanopore reads were used for genome assembly. The taxonomic assignation was performed using Centrifuge$^{53}$ version 1.0.3 to detect potential contamination. Genome size and heterozygosity rate were estimated using Genomescope$^{54}$ and Illumina short reads. For the genome assembly, we generated three sets of ONT reads: all the reads, 30× genome coverage with the longest reads and 30× genome coverage of the highest quality reads estimated by the Filtrong tool (https://github.com/rrwick/Filtrong). We used Filtrong filtering with default parameters using *Pelagomonas* Illumina short reads as a reference (ONT reads covered by Illumina reads have higher scores). We then applied four different assemblers with default settings, Smartdenovo, Redbean, Flye and Ra on these three sets of reads (Supplementary Data 1)$^{55–58}$. After the assembly phase, we selected the best assembly (Flye with all reads) based on the cumulative size and fragmentation. Indeed, the Wtdbg2 and Smartdenovo assembler generated fragmented assemblies with lower N90. Raven and Flye were very close, but only the Flye assembly with all ONT reads contained both the mitochondrial and chloroplastic circular contigs. To display connections that are not present in the contigs file, we used Bandage tool$^{59}$. The selected Flye assembly was polished three times using *Racon*$^{60}$ with ONT reads, and two times with HapoG$^{61}$ and Illumina reads. Gene completeness of the assembly was estimated using the single-copy orthologous gene analysis from BUSCO v5 with the stramenopile dataset version 10 containing 100 genes$^{28}$.

**Repeat masking and GC analyses.** Repetitive regions on the genome were masked using Tandem Repeat Finder tool$^{62}$. Dust tool to detect low-complexity regions$^{63}$ and RepeatMasker$^{64}$ to identify interspersed repeats based on homology search within the Stramenopile clade and other low-complexity sequences. The positions of detected repeats were merged and hard-masked on the genome, amounting to 8% of its length. Ab initio identification of repeat family sequences was performed using RepeatScout$^{65}$. The algorithm first calculates the frequency of all k-mers in the genome, then removes low-complexity regions and tandem repeats. In >80% of the cases, repeat families identified using ab initio approaches do not overlap with repetitive regions identified by homology search. GC content along the genome was calculated with Bedtools nuc version 2.29.2$^{66}$ and the coverage over a non-overlapping window of 2 Kb with Mosdepth version 0.2.8$^{67}$.

**Transcriptome assembly.** RNA sequencing reads from *P. calceolata* RCC100 were assembled using Velvet 1.2.07 and Oases 0.2.08 with a k-mer size of 63 bp$^{68,69}$. Reads were mapped back to the contigs with BWA-mem$^{70}$ and only consistent paired-end reads were kept. Uncovered regions were detected and used to identify chimeric contigs. In addition, open reading frames (ORF) and domains were searched using respectively TransDecoder (http://transdecoder.sourceforge.net) and CDDsearch$^{71}$. Contig extremities without predicted ORFs or functional domains were removed. Lastly, we used the read strand information to orient RNA contigs. We completed the RNA contigs dataset with the two transcriptome assemblies of the RCC100 strain of *P. calceolata* from the Marine Eukaryotes Transcriptomes database (METdb) (http://metdb.sb-roscoff.fr/metdb/) $^{72}$.

**Gene prediction.** Nuclear gene prediction was performed using 23,696 Pelagomonadales proteins (mainly *A. anophagefferens*) downloaded from the NCBI website. Proteins were aligned on the genome in a two-step strategy. First, BLAT$^{73}$ (version 36 with default parameters) was used to rapidly localize corresponding putative regions of these proteins on the genome. The best match and the matches with a score greater than or equal to 90% of the best match score were retained. Then, the regions with BLAT alignments were masked and we aligned the same set of proteins using BLAST$^{74}$, which can identify more divergent matches. Second, alignments were refined using Genewise$^{75}$ (version 2.2.0 default parameters, except the -splice model option to detect non-canonical splicing sites), which is more accurate for detecting intron/exon boundaries. Alignments were kept if more than 50% of the length of the protein is aligned on the genome. Additionally, the transcriptome assemblies of *P. calceolata* RCC969, RCC2362, RCC706 and RCC981 included in the METdb were translated into proteins and aligned to the genome using BLAT, a BLAT score > 50 % filter, and alignments refined with Genewise as previously described.

We selected alignments from the newly generated transcriptome assembly and the two assemblies available in METdb belonging to the *P. calceolata* RCC100 strain to build a training set for the AUGUSTUS ab initio gene predictor$^{76}$. Only gene models with complete coding DNA sequences were retained for training and 1000 genes were set aside for testing AUGUSTUS accuracy. Initial training produced exon and intron parameters for *P. calceolata* species. Parameters were optimized using successive steps of training and testing. We calculated gene prediction accuracy by running AUGUSTUS on the test set. At the exon level, AUGUSTUS performed well in terms of sensitivity (0.619) and specificity (0.669). We thus run AUGUSTUS on the masked genome based on trained parameters.

The ab initio prediction and all the transcriptomic and protein alignments were combined using Gmove, an easy-to-use predictor with no need for a pre-calibration
step<sup>77</sup>. Briefly, putative exons and introns, extracted from predictions and alignments, were used to build a graph, where nodes and edges represent exons and introns respectively. From this graph, Gmove extracts all paths and searches open reading frames (ORFs) consistent with the protein evidence. We trimmed untranslated transcribed regions that overlapped the coding part of a neighbor gene and renamed the genes following the standard nomenclature. Mono-exonic genes models encoding proteins of less than 200 amino acids without significant protein match (1006 genes) were excluded. Chloroplast and mitochondrial genes (contig 7 and 8) were predicted using previously published annotations for *P.calceolata*<sup>23,78</sup>. Following this pipeline, we predicted 16,667 genes with 0.45 intron per gene on average.

**Functional annotation.** Predicted gene models of *P. calceolata* nuclear genome (contig 1 to contig 6) were annotated for protein function using InterProScan v5.41-78.0<sup>79</sup>. A protein alignment against the NR database (01-12-2021 version) was performed with diamond v0.9.24<sup>80</sup>. The best protein match with a functional annotation and an *e*-value < 10<sup>-5</sup> was retained. KEGG Orthologues (KO) were identified with the HMM search tool KofamKoala v1.3.0 and KO annotations with an *e*-value < 10<sup>-5</sup> and a score above the HMM threshold were retained<sup>81</sup>. Finally, Gene Ontology (GO) terms and Enzyme commission (EC) numbers were recovered from the Interproscan and KO analysis respectively. Previously published chloroplast and mitochondrial gene names and functions were reported on the corresponding genes<sup>23,78</sup>. All functional gene annotations of *P. calceolata* are available in Supplementary Data 2.

In order to compare the functional annotation of *P. calceolata* with other small free-living photosynthetic eukaryotes, we applied the same analysis to the predicted proteins available for the following species: *Aureococcus anophagefferens*, *Thalassiosira pseudonana*, *Phaeodactylum tricornutum*, *Nannochloropsis oceanica*, *Bathycoccus prasinos*, *Micromonas pusilla*, *Ostreococcus lucimarinus*, and *Emiliania huxleyi* (references are indicated in Table 1).

We defined a list of 23 meiosis-specific genes using three previous studies<sup>31–33</sup>. KO annotations and Interproscan domains were used to recover these genes in the *P. calceolata* genome. Transcriptomic reads of *P. calceolata* were mapped back to the predicted genes with BWA-MEM v 2.2.1. Reads aligned over more than 80% of their length were retained to estimate the relative expression levels of meiosis genes.

Genes containing the NIT-sensing domain were identified based on Interproscan annotations. Eukaryotic homologs of the three *P. calceolata* NIT-sensing genes were retrieved with a BLASTP (*e*-value < 10<sup>-5</sup>, coverage > 100 aa) against 27.7 million proteins from NR, the METDb<sup>72</sup> transcriptome database, eukaryotic algal proteomes from the JGI database, *Tara Oceans* single-cell amplified genomes and metagenome assembled genomes (SMAGs)<sup>82</sup>. The 60 retrieved proteins and the 3 *P. calceolata* NIT-domain-containing proteins were then aligned with Mafft v7.0 (https://mafft.cbcrc.jp/alignment/server/) and a Maximum Likelihood phylogenetic tree (Jones–Taylor–Thornton substitution model and 100 bootstraps) was made with MEGAX software. Transmembrane regions in NIT-sensing domain-containing proteins were identified with TMHMM v 2.055.

**Estimation of *P. calceolata* relative abundance in environmental metagenomic reads.** We used metagenomic datasets of *Tara Oceans* and *Tara Polar Circle* expeditions to detect and estimate the relative abundance of the *P. calceolata* in the oceans. Datasets from water samples collected on the photic zone: surface (SUR) and deep-chlorophyll maximum (DCM) were analyzed. Size-fractionated water samples containing pico-nano algae (organisms < 5 µm in cell diameter) were selected: 0.2–3 µm (100 samples), 0.8–2000 µm (119 samples) and 0.8–5 µm (148 samples)<sup>51</sup>. Metagenomic reads were aligned on the *P. calceolata* genome assembled in this study with BWA-mem version 0.7.15 with default parameters<sup>83</sup>. Alignments with at least 90% of identity over 80% of read length were retained for further analysis. In the case of several possible best matches, a random one was picked. In order to remove putative PCR duplicates, multiple read pairs aligned at the same position on the *P. calceolata* genome were removed with samtools rmmpid version 1.10.2<sup>70</sup>. For the metagenomic abundance, we divided the total number of reads aligned on the *P. calceolata* genome assembled in this study by the total number of sequenced reads for each sample. For the metabarcoding abundance, we used the 18SV9 rRNA OTU table published in 2019 and available here https://zenodo.org/record/3768510#.YEX2S9zjJaQ<sup>34</sup>. Bacterial and archaea OTUs were removed for this analysis.

**P. calceolata** relative abundance models. Two *P. calceolata* relative abundance models were performed based on the in situ environmental conditions measured at the time of sampling or using the WOA18 datasets. The environmental parameters measured during the expedition are available in the Pangaea database (https://www.pangaea.de/) and are described in ref. 84. Iron concentrations are annual means derived from PISCES2 model<sup>85</sup> and described in ref. 13. Ammonium concentrations at the date and location of sampling are derived from the MITgcm Darwin model and available in the Pangaea database<sup>86</sup>. Environmental parameters for each sample are available in Supplementary Data 10. We consider oceanic samples as “low-iron” if they contain less than 0.2 nM of iron, “low-nitrate” if they

contain less than 2 µM of nitrate and “low-ammonium” if they contain less than 25 nM of ammonium. These thresholds were defined with the distribution of nutrient concentrations in the dataset and previous studies<sup>13,14</sup>. Pearson's correlations between the relative abundance of *P. calceolata* and all environmental parameters were calculated with the *cor* function the R package FactoMineR version 2.4 and the GGally package version 2.1.0. Principal component analysis (PCA) was performed with 9 parameters presenting significant Pearson's correlations. We used a Generalized Additive Model (GAM) for its ability to fit non-linear and non-monotonic functions and for its low sensitivity to extreme values to model the relative abundance of *P. calceolata* as a function of iron concentration, temperature and PAR light<sup>87</sup>. This function is implemented in the mgcv R package version 1.8–33.

The mean of several climatologies of the Earth System Models under the RCP8.5 scenario (ESM8.5) was used to define the environmental conditions at the end of the century following the method of ref. 35. *P. calceolata* relative abundance models based on the WOA18 at locations, depths, and months of *Tara Oceans* samples or ESM8.5 were obtained using four machine-learning techniques as described in ref. 35, with some differences. The four tools were trained in regression mode, we used the neuralnet R package<sup>88</sup> with the following parameters: hidden = 1, decay = 0.1, mxit = 1e6 and size = 5, 6 splines were selected for the GAM. We performed 100-fold random cross-validation for each model and evaluated their performance using the Pearson correlation coefficient and rmse. We used the ensemble model approaches<sup>35</sup> for final global-scale models of the relative abundance of *P. calceolata* (i.e., the mean projections of the validated machine-learning techniques). Figures 2a, 3c, and Supplementary Fig. S7 were generated using the world map of R package “maps”, “mapdata” and “ggmap”. World maps data are imported from the public domain Natural Earth project.

**Environmental metatranscriptomic reads mapping and filtering.** Metatranscriptomic reads from *Tara Oceans* datasets were aligned on predicted coding sequences of the *P. calceolata* genome with BWA-mem 2.2.1 using default parameters. We selected reads aligned with more than 95% of identity over 80% of the read length. We kept nuclear genes covered by at least ten reads in a minimum of ten samples and removed those detected in more than 90% of samples which probably aggregate metagenomic reads from other organisms (cross-mapping). To obtain robust relative expression values, we removed samples where less than 25% of *P. calceolata* genes were detected. We finally got the expression of 15,679 genes of *P. calceolata* across 167 samples. In all gene expression analyses, we normalized the gene expression levels in transcripts per kilobase million (TPM). Gene expression levels of *P. calceolata* in all *Tara Oceans* samples are available here: https://doi.org/10.5281/zenodo.6983365.

**Statistical analysis and reproducibility.** Sample sizes are indicated in the methods, in the figure captions or in the main text. Nonparametric Wilcoxon signed-rank tests were applied with the two-sided alternative hypothesis and not paired. Two sides Fisher statistical tests were applied for gene enrichments analysis. All statistical tests in the manuscript were generated with R version 4.0.3 and *P* values < 0.01 are considered significant.

**Reporting summary.** Further information on research design is available in the Nature Research Reporting Summary linked to this article.

## Data availability

*Pelagomonas calceolata* genomic and transcriptomic reads, genome assembly, and gene prediction are available at the ENA (EMBL-EBI) website under the accession number PRJEB47931. *P. calceolata* transcriptomes are available under the accession number PRJEB34158, runs ERR3497221 and ERR3497222. *Tara Oceans* and *Tara Polar Circle* metagenomic sequences are archived at the ENA under the following accession numbers: PRJEB9740, PRJEB9691, PRJEB4352, and PRJEB1787. All other data are available from the corresponding author upon reasonable request.

Received: 13 December 2021; Accepted: 2 September 2022;
Published online: 16 September 2022

## References

1. Field, C. B., Behrenfeld, M. J., Randerson, J. T. & Falkowski, P. Primary production of the biosphere: integrating terrestrial and oceanic components. *Science* **281**, 237–240 (1998).

2. Boyce, D. G., Lewis, M. R. & Worm, B. Global phytoplankton decline over the past century. *Nature* **466**, 591–596 (2010).

3. Henson, S. A., Cael, B. B., Allen, S. R. & Dutkiewicz, S. Future phytoplankton diversity in a changing climate. *Nat. Commun.* **12**, 5372 (2021).
4. Vaulot, D., Eikrem, W., Viprey, M. & Moreau, H. The diversity of small eukaryotic phytoplankton (≤3 µm) in marine ecosystems. *FEMS Microbiol. Rev.* **32**, 795–820 (2008).

5. Agawin, N. S. R., Duarte, C. M. & Agustí, S. Nutrient and temperature control of the contribution of picoplankton to phytoplankton biomass and production. *Limnol. Oceanogr.* **45**, 591–600 (2000).

6. Morán, X. A. G., López-Urrutia, Á., Calvo-Díaz, A. & Li, W. K. W. Increasing importance of small phytoplankton in a warmer ocean. *Glob. Change Biol.* **16**, 1137–1144 (2010).

7. Li, W. K. W., McLaughlin, F. A., Lovejoy, C. & Carmack, E. C. Smallest algae thrive as the arctic ocean freshens. *Science* **326** https://doi.org/10.1126/science.1179798 (2009).

8. Benner, I., Irwin, A. J. & Finkel, Z. V. Capacity of the common Arctic picoeukaryote Micromonas to adapt to a warming ocean. *Limnol. Oceanography Lett.* **5**, 221–227 (2020).

9. Sunda, W. G. & Huntsman, S. A. Iron uptake and growth limitation in oceanic and coastal phytoplankton. *Mar. Chem.* **50**, 189–206 (1995).

10. Raven, J. A. The twelfth Tansley Lecture. Small is beautiful: the picophytoplankton. *Funct. Ecol.* **12**, 503–513 (1998).

11. Morel, F. M. M. & Price, N. M. The biogeochemical cycles of trace metals in the oceans. *Science* **300**, 944–947 (2003).

12. Gao, X., Bowler, C. & Kazamia, E. Iron metabolism strategies in diatoms. *J. Exp. Bot.* **72**, 2165–2180 (2021).

13. Caputi, L. et al. Community-level responses to iron availability in open ocean plankton ecosystems. *Glob. Biogeochemical Cycles* **33**, 391–419 (2019).

14. Carradec, Q. et al. A global ocean atlas of eukaryotic genes. *Nat. Commun.* **9**, 373 (2018).

15. Morrissey, J. et al. A novel protein, ubiquitous in marine phytoplankton, concentrates iron at the cell surface and facilitates uptake. *Curr. Biol.* **25**, 364–371 (2015).

16. Moore, C. M. et al. Processes and patterns of oceanic nutrient limitation. *Nat. Geosci.* **6**, 701–710 (2013).

17. Kumar, A. & Bera, S. Revisiting nitrogen utilization in algae: a review on the process of regulation and assimilation. *Bioresour. Technol. Rep.* **12**, 100584 (2020).

18. Smith, S. R. et al. Evolution and regulation of nitrogen flux through compartmentalized metabolic networks in a marine diatom. *Nat. Commun.* **10**, 4552 (2019).

19. Berg, G. M., Glibert, P. M., Lomas, M. W. & Burford, M. A. Organic nitrogen uptake and growth by the chrysophyte *Aureococcus anophagefferens* during a brown tide event. *Mar. Biol.* **129**, 377–387 (1997).

20. Andersen, R. A., Saunders, G. W., Paskind, M. P. & Sexton, J. P. Ultrastructure and 18s rRNA gene sequence for *Pelagomonas calceolata* gen. et sp. nov. and the description of a new algal class, the pelagophyceae classis nov. *J. Phycol.* **29**, 701–715 (1993).

21. Choi, C. J. et al. Seasonal and geographical transitions in eukaryotic phytoplankton community structure in the Atlantic and Pacific Oceans. *Front. Microbiol.* **11**, 542372 (2020).

22. Duerschlag, J. et al. Niche partitioning by photosynthetic plankton as a driver of CO2-fixation across the oligotrophic South Pacific Subtropical Ocean. *ISME J* **1**–12 https://doi.org/10.1038/s41396-021-01072-z (2021).

23. Worden, A. Z. et al. Global distribution of a wild alga revealed by targeted metagenomics. *Curr. Biol.* **22**, R675–R677 (2012).

24. Dimier, C. é, Brunet, C., Geider, R. & Raven, J. Growth and photoregulation dynamics of the picoeukaryote *Pelagomonas calceolata* in fluctuating light. *Limnol. Oceanogr.* **54**, 823–836 (2009).

25. Dupont, C. L. et al. Genomes and gene expression across light and productivity gradients in eastern subtropical Pacific microbial communities. *ISME J.* **9**, 1076–1092 (2015).

26. Kang, Y. et al. Transcriptomic responses of four pelagophytes to nutrient (N, P) and light stress. *Front. Mar. Sci.* **8**, 636699 (2021).

27. Huff, J. T., Zilberman, D. & Roy, S. W. Mechanism for DNA transposons to generate introns on genomic scales. *Nature* **538**, 533–536 (2016).

28. Waterhouse, R. M. et al. BUSCO applications from quality assessments to gene prediction and phylogenomics. *Mol. Biol. Evol.* **35**, 543–548 (2018).

29. Nambiar, M. & Smith, G. R. Repression of harmful meiotic recombination in centromeric regions. *Semin Cell Dev. Biol.* **54**, 188–197 (2016).

30. Pessia, E. et al. Evidence for widespread GC-biased gene conversion in eukaryotes. *Genome Biol. Evol.* **4**, 675–682 (2012).

31. Chi, J., Mahé, F., Loidl, J., Logsdon, J. & Dunthorn, M. Meiosis gene inventory of four ciliates reveals the prevalence of a synaptonemal complex-independent crossover pathway. *Mol. Biol. Evol.* **31**, 660–672 (2014).

32. Ramesh, M. A., Malik, S.-B. & Logsdon, J. M. A phylogenomic inventory of meiotic genes; evidence for sex in Giardia and an early eukaryotic origin of meiosis. *Curr. Biol.* **15**, 185–191 (2005).

33. Schurko, A. M. & Logsdon, J. M. Using a meiosis detection toolkit to investigate ancient asexual 'scandals' and the evolution of sex. *Bioessays* **30**, 579–589 (2008).

34. Ibarbalz, F. M. et al. Global trends in marine plankton diversity across kingdoms of life. *Cell* **179**, 1084–1097.e21 (2019).

35. Frémont, P. et al. Restructuring of plankton genomic biogeography in the surface ocean under climate change. *Nat. Clim. Chang.* **12**, 393–401 (2022).

36. Ward, D. M. & Kaplan, J. Ferroportin-mediated iron transport: expression and regulation. *Biochim. Biophys. Acta* **1823**, 1426–1433 (2012).

37. Gobler, C. J., Lonsdale, D. J. & Boyer, G. L. A review of the causes, effects, and potential management of harmful brown tide blooms caused by *Aureococcus anophagefferens* (Hargraves et sieburth). *Estuaries* **28**, 726–749 (2005).

38. Agustí, S., Lubián, L. M., Moreno-Ostos, E., Estrada, M. & Duarte, C. M. Projected changes in photosynthetic picoplankton in a warmer subtropical ocean. *Front. Mar. Sci.* **5**, 506 (2019).

39. Anderson, S. I., Barton, A. D., Clayton, S., Duttiewicz, S. & Rynearson, T. A. Marine phytoplankton functional types exhibit diverse responses to thermal change. *Nat. Commun.* **12**, 6413 (2021).

40. Martin, J. H. et al. Testing the iron hypothesis in ecosystems of the equatorial Pacific Ocean. *Nature* **371**, 123–129 (1994).

41. Shi, D., Xu, Y., Hopkinson, B. M. & Morel, F. M. M. Effect of ocean acidification on iron availability to marine phytoplankton. *Science* **327**, 676–679 (2010).

42. McQuaid, J. B. et al. Carbonate-sensitive phytotransferrin controls high-affinity iron uptake in diatoms. *Nature* **555**, 534–537 (2018).

43. Turnšek, J. et al. Proximity proteomics in a marine diatom reveals a putative cell surface-to-chloroplast iron trafficking pathway. *eLife* **10**, e52770 (2021).

44. Urzica, E. I. et al. Systems and trans-system level analysis identifies conserved iron deficiency responses in the plant lineage[W][OA]. *Plant Cell* **24**, 3921–3948 (2012).

45. Mao, X. et al. Diversity, prevalence, and expression of cyanase genes (cynS) in planktonic marine microorganisms. *ISME J.* **16**, 602–605 (2022).

46. Ou, L., Cai, Y., Jin, W., Wang, Z. & Lu, S. Understanding the nitrogen uptake and assimilation of the Chinese strain of *Aureococcus anophagefferens* (Pelagophyceae). *Algal Res.* **34**, 182–190 (2018).

47. Shu, C. J., Ulrich, L. E. & Zhulin, I. B. The NIT domain: a predicted nitrate-responsive module in bacterial sensory receptors. *Trends Biochem Sci.* **28**, 121–124 (2003).

48. Wu, S. Q., Chai, W., Lin, J. T. & Stewart, V. General nitrogen regulation of nitrate assimilation regulatory gene nasR expression in Klebsiella oxytoca M5al. *J. Bacteriol.* **181**, 7274–7289 (1994).

49. Durand, N. C. et al. Juicer provides a one-click system for analyzing loop-resolution Hi-C experiments. *Cell Syst.* **3**, 95–98 (2016).

50. Li, R., Li, Y., Kristiansen, K. & Wang, J. SOAP: short oligonucleotide alignment program. *Bioinformatics* **24**, 713–714 (2008).

51. Alberti, A. et al. Viral to metazoa marine plankton nucleotide sequences from the Tara Oceans expedition. *Sci. Data* **4**, 170093 (2017).

52. Kopylova, E., Noé, L. & Touzet, H. SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data. *Bioinformatics* **28**, 3211–3217 (2012).

53. Kim, D., Song, L., Breitwieser, F. P. & Salzberg, S. L. Centrifuge: rapid and sensitive classification of metagenomic sequences. *Genome Res.* https://doi.org/10.1101/gr.210641.116 (2016).

54. Vulture, G. W. et al. GenomeScope: fast reference-free genome profiling from short reads. *Bioinformatics* **33**, 2202–2204 (2017).

55. Vaser, R. & Šikić, M. Yet another de novo genome assembler. *BioRxiv*. https://doi.org/10.1101/656306 (2019).

56. Liu, H. et al. SMARTdenovo: a de novo assembler using long noisy reads. *Gigabyte* **2021**, 1–9 (2021).

57. Kolmogorov, M., Yuan, J., Lin, Y. & Pevzner, P. A. Assembly of long, error-prone reads using repeat graphs. *Nat. Biotechnol.* **37**, 540–546 (2019).

58. Ruan, J. & Li, H. Fast and accurate long-read assembly with wtdbg2. *Nat. Methods* **17**, 155–158 (2020).

59. Wick, R. R., Schultz, M. B., Zobel, J. & Holt, K. E. Bandage: interactive visualization of de novo genome assemblies. *Bioinformatics* **31**, 3350–3352 (2015).

60. Vaser, R., Sović, I., Nagarajan, N. & Šikić, M. Fast and accurate de novo genome assembly from long uncorrected reads. *Genome Res* **27**, 737–746 (2017).

61. Aury, J.-M. & Istace, B. Hapo-G, haplotype-aware polishing of genome assemblies with accurate reads. *NAR Genomics Bioinform.* **3**, lqab034 (2021).

62. Benson, G. Tandem repeats finder: a program to analyze DNA sequences. *Nucleic Acids Res.* **27**, 573–580 (1999).

63. Morgulis, A., Gertz, E. M., Schäffer, A. A. & Agarwala, R. A fast and symmetric DUST implementation to mask low-complexity DNA sequences. *J. Comput Biol.* **13**, 1028–1040 (2006).

64. Smit, A. F. A., Hubley, R. & Green, P. RepeatMasker. http://repeatmasker.org/ (2013).

65. Price, A. L., Jones, N. C. & Pevzner, P. A. De novo identification of repeat families in large genomes. *Bioinformatics* **21**, i351–i358 (2005).
66. Quinlan, A. R. & Hall, I. M. BEDTools: a flexible suite of utilities for comparing genomic features. *Bioinformatics* **26**, 841–842 (2010).

67. Pedersen, B. S. & Quinlan, A. R. Mosdepth: quick coverage calculation for genomes and exomes. *Bioinformatics* **34**, 867–868 (2018).

68. Schulz, M. H., Zerbino, D. R., Vingron, M. & Birney, E. Oases: robust de novo RNA-seq assembly across the dynamic range of expression levels. *Bioinformatics* **28**, 1086–1092 (2012).

69. Zerbino, D. R. & Birney, E. Velvet: algorithms for de novo short read assembly using de Brujin graphs. *Genome Res.* **18**, 821–829 (2008).

70. Li, H. et al. The sequence Alignment/Map format and SAMtools. *Bioinformatics* **25**, 2078–2079 (2009).

71. Marchler-Bauer, A. et al. CDD: NCBI's conserved domain database. *Nucleic Acids Res.* **43**, D222–D226 (2015).

72. Niang, G. et al. METdb: A genomic reference database for marine species. *F1000Research*, https://doi.org/10.7490/f1000research.1118000.1 (2020).

73. Kent, W. J. BLAT—the BLAST-like alignment tool. *Genome Res.* **12**, 656–664 (2002).

74. Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic local alignment search tool. *J. Mol. Biol.* **215**, 403–410 (1990).

75. Birney, E., Clamp, M. & Durbin, R. GeneWise and genomewise. *Genome Res.* **14**, 988–995 (2004).

76. Stanke, M. et al. AUGUSTUS: ab initio prediction of alternative transcripts. *Nucleic Acids Res.* **34**, W435–W439 (2006).

77. Dubarry, M. et al. Gmove a tool for eukaryotic gene predictions using various evidences. *F1000Research*, https://doi.org/10.7490/f1000research.1111735.1 (2016).

78. Sibbald, S. J., Lawton, M. & Archibald, J. M. Mitochondrial genome evolution in pelagophyte algae. *Genome Biol. Evol.* **13**, evab018 (2021).

79. Quevillon, E. et al. InterProScan: protein domains identifier. *Nucleic Acids Res.* **33**, W116–W120 (2005).

80. Buchfink, B., Reuter, K. & Drost, H.-G. Sensitive protein alignments at tree-of-life scale using DIAMOND. *Nat. Methods* **18**, 366–368 (2021).

81. Aramaki, T. et al. KofamKOALA: KEGG Ortholog assignment based on profile HMM and adaptive score threshold. *Bioinformatics* **36**, 2251–2252 (2020).

82. Delmont, T. O. et al. Functional repertoire convergence of distantly related eukaryotic plankton lineages abundant in the sunlit ocean. *Cell Genomics* **2**, 100123 (2022).

83. Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. *Bioinformatics* **25**, 1754–1760 (2009).

84. Pesant, S. et al. Open science resources for the discovery and analysis of Tara Oceans data. *Sci. Data* **2**, 150023 (2015).

85. Aumont, O., Ethé, C., Tagliabue, A., Bopp, L. & Gehlen, M. PICES-v2: an ocean biogeochemical model for carbon and ecosystem studies. *Geoscientific Model Dev.* **8**, 2465–2513 (2015).

86. Clayton, S. et al. Biogeochemical versus ecological consequences of modeled ocean physics. *Biogeosciences* **14**, 2877–2889 (2017).

87. Ravindra, K., Rattan, P., Mor, S. & Aggarwal, A. N. Generalized additive models: building evidence of air pollution, climate change and human health. *Environ. Int.* **132**, 104987 (2019).

88. Günther, F. & Fritsch, S. neuralnet: training of neural networks. *R. J.* **2**, 30–38 (2010).

89. Gobler, C. J. et al. Niche of harmful alga *Aureococcus anophagefferens* revealed through ecogenomics. *Proc. Natl Acad. Sci. USA* **108**, 4352–4357 (2011).

90. Guo, L. et al. Genome assembly of *Nannochloropsis oceanica* provides evidence of host nucleus overthrow by the symbiont nucleus during speciation. *Commun. Biol.* **2**, 1–12 (2019).

91. Bowler, C. et al. The Phaeodactylum genome reveals the evolutionary history of diatom genomes. *Nature* **456**, 239–244 (2008).

92. Armbrust, E. V. et al. The genome of the diatom thalassiosira pseudonana: ecology, evolution, and metabolism. *Science* **306**, 79–86 (2004).

93. Worden, A. Z. et al. Green evolution and dynamic adaptations revealed by genomes of the marine picoeukaryotes micromonas. *Science* **324**, 268–272 (2009).

94. Palenik, B. et al. The tiny eukaryote Ostreococcus provides genomic insights into the paradox of plankton speciation. *PNAS* **104**, 7705–7710 (2007).

95. Moreau, H. et al. Gene functionalities and genome structure in Bathycoccus prasinos reflect cellular specializations at the base of the green lineage. *Genome Biol.* **13**, R74 (2012).

96. Read, B. A. et al. Pan genome of the phytoplankton *Emiliania underpins* its global distribution. *Nature* **499**, 209–213 (2013).

## Acknowledgements

We thank the commitment of the following people who made this work possible: the Genoscope/CEA, the CNRS (in particular the Federation de Recherche R2022/Tara Oceans GO-SEE), Marie-José Garet-Delmas from the Roscoff Culture Collection for growing the RCC100 strain, Claude Scarpelli for support in high-performance computing and Linda Sperling for language editing. Computations were performed using the cobalt HPC machine. We acknowledge the financial support of FRANCE GENOMIQUE (ANR-10-INBS-09-08) and Oceanomics (ANR-11-BTBR-0008). We also thank the Tara Expedition Foundation and their partners for the organization of marine scientific expeditions (http://oceans.taraexpeditions.org). This article is contribution number 139 of Tara Oceans.

## Author contributions

S.R., C.B., and M.G. performed *P. calceolata* cultures and DNA/RNA extractions. A.A., E.P., and C.C. coordinated DNA/RNA sequencing. B.I, B.N., and J.M.A. did the assembly and annotation of the genome. M.C., S.M., Q.C., and J.M.A. carried out genomic analysis. N.G., E.F., and Q.C. analyzed environmental datasets. P.F. and O.J. worked on *P. calceolata* models. Q.C. and N.G. wrote the paper with a strong support of J.M.A. and P.W. C.B, L.B., and K.L. did the Hi-C library preparation, sequencing and analysis, respectively. All authors contributed to the manuscript preparation and approved the final version of the paper.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s42003-022-03939-z.

**Correspondence** and requests for materials should be addressed to Quentin Carradec.

**Peer review information** *Communications Biology* thanks Senjie Lin, Ana Cabello and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Primary Handling Editors: Linn Hoffmann and Luke R. Grinham. Peer reviewer reports are available.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022