Genome Resolved Biogeography of Mamiellales

The tiny giant of the sea, Ostreococcus's unique adaptations

Ostreococcus spp. are unicellular organisms with one of the simplest cellular organizations. The sequencing of the genomes of different Ostreococcus species has reinforced this status since Ostreococcus tauri has one most compact nuclear genomes among eukaryotic organisms. Despite this, it has retained a number of genes, setting it apart from other organisms with similar small genomes. Ostreococcus spp. feature a substantial number of selenocysteine-containing proteins, which, due to their higher catalytic activity compared to their selenium-lacking counterparts, may require a reduced quantity of proteins. Notably, O. tauri encodes several ammonium transporter genes, that may provide it with a competitive edge for acquiring nitrogen (N). This characteristic makes it an intriguing model for studying the efficient use of N in eukaryotes. Under conditions of low N availability, O. tauri utilizes N from abundant proteins or amino acids, such as L-arginine, similar to higher plants. However, the presence of a nitric oxide synthase (L-arg substrate) sheds light on a new metabolic pathway for L-arg in algae. The metabolic adaptations of O. tauri to day and night cycles offer valuable insights into carbon and iron metabolic configuration. O. tauri has evolved novel strategies to optimize iron uptake, lacking the classic components of the iron absorption mechanism. Overall, the cellular and genetic characteristics of Ostreococcus contribute to its evolutionary success, making it an excellent model for studying the physiological and genetic aspects of how green algae have adapted to the marine environment. Furthermore, given its potential for lipid accumulation and its marine habitat, it may represent a promising avenue for third-generation biofuels.

Genomic characterisation and ecological distribution of Mantoniella tinhauana: a novel Mamiellophycean green alga from the Western Pacific

Mamiellophyceae are dominant marine algae in much of the ocean, the most prevalent genera belonging to the order Mamiellales: Micromonas, Ostreococcus and Bathycoccus, whose genetics and global distributions have been extensively studied. Conversely, the genus Mantoniella, despite its potential ecological importance, remains relatively under-characterised. In this study, we isolated and characterised a novel species of Mamiellophyceae, Mantoniella tinhauana, from subtropical coastal waters in the South China Sea. Morphologically, it resembles other Mantoniella species; however, a comparative analysis of the 18S and ITS2 marker genes revealed its genetic distinctiveness. Furthermore, we sequenced and assembled the first genome of Mantoniella tinhauana, uncovering significant differences from previously studied Mamiellophyceae species. Notably, the genome lacked any detectable outlier chromosomes and exhibited numerous unique orthogroups. We explored gene groups associated with meiosis, scale and flagella formation, shedding light on species divergence, yet further investigation is warranted. To elucidate the biogeography of Mantoniella tinhauana, we conducted a comprehensive analysis using global metagenomic read mapping to the newly sequenced genome. Our findings indicate this species exhibits a cosmopolitan distribution with a low-level prevalence worldwide. Understanding the intricate dynamics between Mamiellophyceae and the environment is crucial for comprehending their impact on the ocean ecosystem and accurately predicting their response to forthcoming environmental changes.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Comparative experimental evolution reveals species-specific idiosyncrasies in marine phytoplankton adaptation to warming

A number of experimental studies have demonstrated that phytoplankton can display rapid thermal adaptation in response to warmed environments. While these studies provide insight into the evolutionary responses of single species, they tend to employ different experimental techniques. Consequently, our ability to compare the potential for thermal adaptation across different, ecologically relevant, species remains limited. Here, we address this limitation by conducting simultaneous long-term warming experiments with the same experimental design on clonal isolates of three phylogenetically diverse species of marine phytoplankton; the cyanobacterium Synechococcus sp., the prasinophyte Ostreococcus tauri and the diatom Phaeodoactylum tricornutum. Over the same experimental time period, we observed differing levels of thermal adaptation in response to stressful supra-optimal temperatures. Synechococcus sp. displayed the greatest improvement in fitness (i.e., growth rate) and thermal tolerance (i.e., temperature limits of growth). Ostreococcus tauri was able to improve fitness and thermal tolerance, but to a lesser extent. Finally, Phaeodoactylum tricornutum showed no signs of adaptation. These findings could help us understand how the structure of phytoplankton communities may change in response to warming, and possible biogeochemical implications, as some species show relatively more rapid adaptive shifts in their thermal tolerance.

Structural Characterization of the Chlorophyllide a Oxygenase (CAO) Enzyme Through an In Silico Approach

Chlorophyllide a oxygenase (CAO) is responsible for converting chlorophyll a to chlorophyll b in a two-step oxygenation reaction. CAO belongs to the family of Rieske-mononuclear iron oxygenases. Although the structure and reaction mechanism of other Rieske monooxygenases have been described, a member of plant Rieske non-heme iron-dependent monooxygenase has not been structurally characterized. The enzymes in this family usually form a trimeric structure and electrons are transferred between the non-heme iron site and the Rieske center of the adjoining subunits. CAO is supposed to form a similar structural arrangement. However, in Mamiellales such as Micromonas and Ostreococcus, CAO is encoded by two genes where non-heme iron site and Rieske cluster localize on the distinct polypeptides. It is not clear if they can form a similar structural organization to achieve the enzymatic activity. In this study, the tertiary structures of CAO from the model plant Arabidopsis thaliana and the Prasinophyte Micromonas pusilla were predicted by deep learning-based methods, followed by energy minimization and subsequent stereochemical quality assessment of the predicted models. Furthermore, the chlorophyll a binding cavity and the interaction of ferredoxin, which is the electron donor, on the surface of Micromonas CAO were predicted. The electron transfer pathway was predicted in Micromonas CAO and the overall structure of the CAO active site was conserved even though it forms a heterodimeric complex. The structures presented in this study will serve as a basis for understanding the reaction mechanism and regulation of the plant monooxygenase family to which CAO belongs.

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

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

Whole-genome scanning reveals environmental selection mechanisms that shape diversity in populations of the epipelagic diatom Chaetoceros

Diatoms form a diverse and abundant group of photosynthetic protists that are essential players in marine ecosystems. However, the microevolutionary structure of their populations remains poorly understood, particularly in polar regions. Exploring how closely related diatoms adapt to different environments is essential given their short generation times, which may allow rapid adaptations, and their prevalence in marine regions dramatically impacted by climate change, such as the Arctic and Southern Oceans. Here, we address genetic diversity patterns in Chaetoceros, the most abundant diatom genus and one of the most diverse, using 11 metagenome-assembled genomes (MAGs) reconstructed from Tara Oceans metagenomes. Genome-resolved metagenomics on these MAGs confirmed a prevalent distribution of Chaetoceros in the Arctic Ocean with lower dispersal in the Pacific and Southern Oceans as well as in the Mediterranean Sea. Single-nucleotide variants identified within the different MAG populations allowed us to draw a landscape of Chaetoceros genetic diversity and revealed an elevated genetic structure in some Arctic Ocean populations. Gene flow patterns of closely related Chaetoceros populations seemed to correlate with distinct abiotic factors rather than with geographic distance. We found clear positive selection of genes involved in nutrient availability responses, in particular for iron (e.g., ISIP2a, flavodoxin), silicate, and phosphate (e.g., polyamine synthase), that were further supported by analysis of Chaetoceros transcriptomes. Altogether, these results highlight the importance of environmental selection in shaping diatom diversity patterns and provide new insights into their metapopulation genomics through the integration of metagenomic and environmental data.

Long-Term Stability of Bacterial Associations in a Microcosm of Ostreococcus tauri (Chlorophyta, Mamiellophyceae)

Phytoplankton-bacteria interactions rule over carbon fixation in the sunlit ocean, yet only a handful of phytoplanktonic-bacteria interactions have been experimentally characterized. In this study, we investigated the effect of three bacterial strains isolated from a long-term microcosm experiment with one Ostreococcus strain (Chlorophyta, Mamiellophyceae). We provided evidence that two Roseovarius strains (Alphaproteobacteria) had a beneficial effect on the long-term survival of the microalgae whereas one Winogradskyella strain (Flavobacteriia) led to the collapse of the microalga culture. Co-cultivation of the beneficial and the antagonistic strains also led to the loss of the microalga cells. Metagenomic analysis of the microcosm is consistent with vitamin B12 synthesis by the Roseovarius strains and unveiled two additional species affiliated to Balneola (Balneolia) and Muricauda (Flavobacteriia), which represent less than 4% of the reads, whereas Roseovarius and Winogradskyella recruit 57 and 39% of the reads, respectively. These results suggest that the low-frequency bacterial species may antagonize the algicidal effect of Winogradskyella in the microbiome of Ostreococcus tauri and thus stabilize the microalga persistence in the microcosm. Altogether, these results open novel perspectives into long-term stability of phytoplankton cultures.

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.

Evolutionary Genomics of Sex-Related Chromosomes at the Base of the Green Lineage

Although sex is now accepted as a ubiquitous and ancestral feature of eukaryotes, direct observation of sex is still lacking in most unicellular eukaryotic lineages. Evidence of sex is frequently indirect and inferred from the identification of genes involved in meiosis from whole genome data and/or the detection of recombination signatures from genetic diversity in natural populations. In haploid unicellular eukaryotes, sex-related chromosomes are named mating-type (MTs) chromosomes and generally carry large genomic regions where recombination is suppressed. These regions have been characterized in Fungi and Chlorophyta and determine gamete compatibility and fusion. Two candidate MT+ and MT− alleles, spanning 450–650 kb, have recently been described in Ostreococcus tauri, a marine phytoplanktonic alga from the Mamiellophyceae class, an early diverging branch in the green lineage. Here, we investigate the architecture and evolution of these candidate MT+ and MT− alleles. We analyzed the phylogenetic profile and GC content of MT gene families in eight different genomes whose divergence has been previously estimated at up to 640 Myr, and found evidence that the divergence of the two MT alleles predates speciation in the Ostreococcus genus. Phylogenetic profiles of MT trans-specific polymorphisms in gametologs disclosed candidate MTs in two additional species, and possibly a third. These Mamiellales MT candidates are likely to be the oldest mating-type loci described to date, which makes them fascinating models to investigate the evolutionary mechanisms of haploid sex determination in eukaryotes.

Tiara: deep learning-based classification system for eukaryotic sequences

MOTIVATION

With a large number of metagenomic datasets becoming available, eukaryotic metagenomics emerged as a new challenge. The proper classification of eukaryotic nuclear and organellar genomes is an essential step toward a better understanding of eukaryotic diversity.

RESULTS

We developed Tiara, a deep-learning-based approach for the identification of eukaryotic sequences in the metagenomic datasets. Its two-step classification process enables the classification of nuclear and organellar eukaryotic fractions and subsequently divides organellar sequences into plastidial and mitochondrial. Using the test dataset, we have shown that Tiara performed similarly to EukRep for prokaryotes classification and outperformed it for eukaryotes classification with lower calculation time. In the tests on the real data, Tiara performed better than EukRep in analyzing the small dataset representing eukaryotic cell microbiome and large dataset from the pelagic zone of oceans. Tiara is also the only available tool correctly classifying organellar sequences, which was confirmed by the recovery of nearly complete plastid and mitochondrial genomes from the test data and real metagenomic data.

AVAILABILITY AND IMPLEMENTATION

Tiara is implemented in python 3.8, available at https://github.com/ibe-uw/tiara and tested on Unix-based systems. It is released under an open-source MIT license and documentation is available at https://ibe-uw.github.io/tiara. Version 1.0.1 of Tiara has been used for all benchmarks.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Combining Nanopore and Illumina Sequencing Permits Detailed Analysis of Insertion Mutations and Structural Variations Produced by PEG-Mediated Transformation in Ostreococcus tauri

Ostreococcus tauri is a simple unicellular green alga representing an ecologically important group of phytoplankton in oceans worldwide. Modern molecular techniques must be developed in order to understand the mechanisms that permit adaptation of microalgae to their environment. We present for the first time in O. tauri a detailed characterization of individual genomic integration events of foreign DNA of plasmid origin after PEG-mediated transformation. Vector integration occurred randomly at a single locus in the genome and mainly as a single copy. Thus, we confirmed the utility of this technique for insertional mutagenesis. While the mechanism of double-stranded DNA repair in the O. tauri model remains to be elucidated, we clearly demonstrate by genome resequencing that the integration of the vector leads to frequent structural variations (deletions/insertions and duplications) and some chromosomal rearrangements in the genome at the insertion loci. Furthermore, we often observed variations in the vector sequence itself. From these observations, we speculate that a nonhomologous end-joining-like mechanism is employed during random insertion events, as described in plants and other freshwater algal models. PEG-mediated transformation is therefore a promising molecular biology tool, not only for functional genomic studies, but also for biotechnological research in this ecologically important marine alga.