PhytoREF: a reference database of the plastidial 16S rRNA gene of photosynthetic eukaryotes with curated taxonomy

Photosynthetic eukaryotes have a critical role as the main producers in most ecosystems of the biosphere. The ongoing environmental metabarcoding revolution opens the perspective for holistic ecosystems biological studies of these organisms, in particular the unicellular microalgae that often lack distinctive morphological characters and have complex life cycles. To interpret environmental sequences, metabarcoding necessarily relies on taxonomically curated databases containing reference sequences of the targeted gene (or barcode) from identified organisms. To date, no such reference framework exists for photosynthetic eukaryotes. In this study, we built the PhytoREF database that contains 6490 plastidial 16S rDNA reference sequences that originate from a large diversity of eukaryotes representing all known major photosynthetic lineages. We compiled 3333 amplicon sequences available from public databases and 879 sequences extracted from plastidial genomes, and generated 411 novel sequences from cultured marine microalgal strains belonging to different eukaryotic lineages. A total of 1867 environmental Sanger 16S rDNA sequences were also included in the database. Stringent quality filtering and a phylogeny-based taxonomic classification were applied for each 16S rDNA sequence. The database mainly focuses on marine microalgae, but sequences from land plants (representing half of the PhytoREF sequences) and freshwater taxa were also included to broaden the applicability of PhytoREF to different aquatic and terrestrial habitats. PhytoREF, accessible via a web interface (http://phytoref.fr), is a new resource in molecular ecology to foster the discovery, assessment and monitoring of the diversity of photosynthetic eukaryotes using high-throughput sequencing.

PhytoREF: a reference database of the plastidial 16S rRNA gene of photosynthetic eukaryotes with curated taxonomy

Photosynthetic eukaryotes have a critical role as the main producers in most ecosystems of the biosphere. The ongoing environmental metabarcoding revolution opens the perspective for holistic ecosystems biological studies of these organisms, in particular the unicellular microalgae that often lack distinctive morphological characters and have complex life cycles. To interpret environmental sequences, metabarcoding necessarily relies on taxonomically curated databases containing reference sequences of the targeted gene (or barcode) from identified organisms. To date, no such reference framework exists for photosynthetic eukaryotes. In this study, we built the PhytoREF database that contains 6490 plastidial 16S rDNA reference sequences that originate from a large diversity of eukaryotes representing all known major photosynthetic lineages. We compiled 3333 amplicon sequences available from public databases and 879 sequences extracted from plastidial genomes, and generated 411 novel sequences from cultured marine microalgal strains belonging to different eukaryotic lineages. A total of 1867 environmental Sanger 16S rDNA sequences were also included in the database. Stringent quality filtering and a phylogeny-based taxonomic classification were applied for each 16S rDNA sequence. The database mainly focuses on marine microalgae, but sequences from land plants (representing half of the PhytoREF sequences) and freshwater taxa were also included to broaden the applicability of PhytoREF to different aquatic and terrestrial habitats. PhytoREF, accessible via a web interface (http://phytoref.fr), is a new resource in molecular ecology to foster the discovery, assessment and monitoring of the diversity of photosynthetic eukaryotes using high-throughput sequencing.

Recent Reticulate Evolution in the Ecologically Dominant Lineage of Coccolithophores

The coccolithophore family Noëlaerhabdaceae contains a number of taxa that are very abundant in modern oceans, including the cosmopolitan bloom-forming Emiliania huxleyi. Introgressive hybridization has been suggested to account for incongruences between nuclear, mitochondrial and plastidial phylogenies of morphospecies within this lineage, but the number of species cultured to date remains rather limited. Here, we present the characterization of 5 new Noëlaerhabdaceae culture strains isolated from samples collected in the south-east Pacific Ocean. These were analyzed morphologically using scanning electron microscopy and phylogenetically by sequencing 5 marker genes (nuclear 18S and 28S rDNA, plastidial tufA, and mitochondrial cox1 and cox3 genes). Morphologically, one of these strains corresponded to Gephyrocapsa ericsonii and the four others to Reticulofenestra parvula. Ribosomal gene sequences were near identical between these new strains, but divergent from G. oceanica, G. muellerae, and E. huxleyi. In contrast to the clear distinction in ribosomal phylogenies, sequences from other genomic compartments clustered with those of E. huxleyi strains with which they share an ecological range (i.e., warm temperate to tropical waters). These data provide strong support for the hypothesis of past (and potentially ongoing) introgressive hybridization within this ecologically important lineage and for the transfer of R. parvula to Gephyrocapsa. These results have important implications for understanding the role of hybridization in speciation in vast ocean meta-populations of phytoplankton.

Recent Reticulate Evolution in the Ecologically Dominant Lineage of Coccolithophores

The coccolithophore family Noëlaerhabdaceae contains a number of taxa that are very abundant in modern oceans, including the cosmopolitan bloom-forming Emiliania huxleyi. Introgressive hybridization has been suggested to account for incongruences between nuclear, mitochondrial and plastidial phylogenies of morphospecies within this lineage, but the number of species cultured to date remains rather limited. Here, we present the characterization of 5 new Noëlaerhabdaceae culture strains isolated from samples collected in the south-east Pacific Ocean. These were analyzed morphologically using scanning electron microscopy and phylogenetically by sequencing 5 marker genes (nuclear 18S and 28S rDNA, plastidial tufA, and mitochondrial cox1 and cox3 genes). Morphologically, one of these strains corresponded to Gephyrocapsa ericsonii and the four others to Reticulofenestra parvula. Ribosomal gene sequences were near identical between these new strains, but divergent from G. oceanica, G. muellerae, and E. huxleyi. In contrast to the clear distinction in ribosomal phylogenies, sequences from other genomic compartments clustered with those of E. huxleyi strains with which they share an ecological range (i.e., warm temperate to tropical waters). These data provide strong support for the hypothesis of past (and potentially ongoing) introgressive hybridization within this ecologically important lineage and for the transfer of R. parvula to Gephyrocapsa. These results have important implications for understanding the role of hybridization in speciation in vast ocean meta-populations of phytoplankton.

Morphological and Phylogenetic Characterization of New Gephyrocapsa Isolates Suggests Introgressive Hybridization in the Emiliania/Gephyrocapsa Complex (Haptophyta)

The coccolithophore genus Gephyrocapsa contains a cosmopolitan assemblage of pelagic species, including the bloom-forming Gephyrocapsa oceanica, and is closely related to the emblematic coccolithophore Emiliania huxleyi within the Noëlaerhabdaceae. These two species have been extensively studied and are well represented in culture collections, whereas cultures of other species of this family are lacking. We report on three new strains of Gephyrocapsa isolated into culture from samples from the Chilean coastal upwelling zone using a novel flow cytometric single-cell sorting technique. The strains were characterized by morphological analysis using scanning electron microscopy and phylogenetic analysis of 6 genes (nuclear 18S and 28S rDNA, plastidial 16S and tufA, and mitochondrial cox1 and cox3 genes). Morphometric features of the coccoliths indicate that these isolates are distinct from G. oceanica and best correspond to G. muellerae. Surprisingly, both plastidial and mitochondrial gene phylogenies placed these strains within the E. huxleyi clade and well separated from G. oceanica isolates, making Emiliania appear polyphyletic. The only nuclear sequence difference, 1 bp in the 28S rDNA region, also grouped E. huxleyi with the new Gephyrocapsa isolates and apart from G. oceanica. Specifically, the G. muellerae morphotype strains clustered with the mitochondrial β clade of E. huxleyi, which, like G. muellerae, has been associated with cold (temperate and sub-polar) waters. Among putative evolutionary scenarios that could explain these results we discuss the possibility that E. huxleyi is not a valid taxonomic unit, or, alternatively the possibility of past hybridization and introgression between each E. huxleyi clade and older Gephyrocapsa clades. In either case, the results support the transfer of Emiliania to Gephyrocapsa. These results have important implications for relating morphological species concepts to ecological and evolutionary units of diversity.

Morphological and Phylogenetic Characterization of New Gephyrocapsa Isolates Suggests Introgressive Hybridization in the Emiliania/Gephyrocapsa Complex (Haptophyta)

The coccolithophore genus Gephyrocapsa contains a cosmopolitan assemblage of pelagic species, including the bloom-forming Gephyrocapsa oceanica, and is closely related to the emblematic coccolithophore Emiliania huxleyi within the Noëlaerhabdaceae. These two species have been extensively studied and are well represented in culture collections, whereas cultures of other species of this family are lacking. We report on three new strains of Gephyrocapsa isolated into culture from samples from the Chilean coastal upwelling zone using a novel flow cytometric single-cell sorting technique. The strains were characterized by morphological analysis using scanning electron microscopy and phylogenetic analysis of 6 genes (nuclear 18S and 28S rDNA, plastidial 16S and tufA, and mitochondrial cox1 and cox3 genes). Morphometric features of the coccoliths indicate that these isolates are distinct from G. oceanica and best correspond to G. muellerae. Surprisingly, both plastidial and mitochondrial gene phylogenies placed these strains within the E. huxleyi clade and well separated from G. oceanica isolates, making Emiliania appear polyphyletic. The only nuclear sequence difference, 1 bp in the 28S rDNA region, also grouped E. huxleyi with the new Gephyrocapsa isolates and apart from G. oceanica. Specifically, the G. muellerae morphotype strains clustered with the mitochondrial β clade of E. huxleyi, which, like G. muellerae, has been associated with cold (temperate and sub-polar) waters. Among putative evolutionary scenarios that could explain these results we discuss the possibility that E. huxleyi is not a valid taxonomic unit, or, alternatively the possibility of past hybridization and introgression between each E. huxleyi clade and older Gephyrocapsa clades. In either case, the results support the transfer of Emiliania to Gephyrocapsa. These results have important implications for relating morphological species concepts to ecological and evolutionary units of diversity.