The distribution and diversity of eukaryotic phytoplankton in the Icelandic marine environment.

Phytoplankton play a crucial role in the marine food web and are sensitive indicators of environmental change. Iceland is at the center of a contrasting hydrography, with cold Arctic water coming in from the north and warmer Atlantic water from the south, making this geographical location very sensitive to climate change. We used DNA metabarcoding to determine the biogeography of phytoplankton in this area of accelerating change. Seawater samples were collected in spring (2012-2018), summer (2017) and winter (2018) together with corresponding physico-chemical metadata around Iceland. Amplicon sequencing of the V4 region of the 18S rRNA gene indicates that eukaryotic phytoplankton community composition is different between the northern and southern water masses, with some genera completely absent from Polar Water masses. Emiliania was more dominant in the Atlantic-influenced waters and in summer, and Phaeocystis was more dominant in the colder, northern waters and in winter. The Chlorophyta picophytoplankton genus, Micromonas, was similarly dominant to the dominant diatom genus, Chaetoceros. This study presents an extensive dataset which can be linked with other 18s rRNA datasets for further investigation into the diversity and biogeography of marine protists in the North Atlantic.

Global patterns and rates of habitat transitions across the eukaryotic tree of life.

The successful colonization of new habitats has played a fundamental role during the evolution of life. Salinity is one of the strongest barriers for organisms to cross, which has resulted in the evolution of distinct marine and non-marine (including both freshwater and soil) communities. Although microbes represent by far the vast majority of eukaryote diversity, the role of the salt barrier in shaping the diversity across the eukaryotic tree is poorly known. Traditional views suggest rare and ancient marine/non-marine transitions but this view is being challenged by the discovery of several recently transitioned lineages. Here, we investigate habitat evolution across the tree of eukaryotes using a unique set of taxon-rich phylogenies inferred from a combination of long-read and short-read environmental metabarcoding data spanning the ribosomal DNA operon. Our results show that, overall, marine and non-marine microbial communities are phylogenetically distinct but transitions have occurred in both directions in almost all major eukaryotic lineages, with hundreds of transition events detected. Some groups have experienced relatively high rates of transitions, most notably fungi for which crossing the salt barrier has probably been an important aspect of their successful diversification. At the deepest phylogenetic levels, ancestral habitat reconstruction analyses suggest that eukaryotes may have first evolved in non-marine habitats and that the two largest known eukaryotic assemblages (TSAR and Amorphea) arose in different habitats. Overall, our findings indicate that the salt barrier has played an important role during eukaryote evolution and provide a global perspective on habitat transitions in this domain of life.

metaPR2 : A database of eukaryotic 18S rRNA metabarcodes with an emphasis on protists.

In recent years, metabarcoding has become the method of choice for investigating the composition and assembly of microbial eukaryotic communities. The number of environmental data sets published has increased very rapidly. Although unprocessed sequence files are often publicly available, processed data, in particular clustered sequences, are rarely available in a usable format. Clustered sequences are reported as operational taxonomic units (OTUs) with different similarity levels or more recently as amplicon sequence variants (ASVs). This hampers comparative studies between different environments and data sets, for example examining the biogeographical patterns of specific groups/species, as well analysing the genetic microdiversity within these groups. Here, we present a newly-assembled database of processed 18S rRNA metabarcodes that are annotated with the PR2 reference sequence database. This database, called metaPR2 , contains 41 data sets corresponding to more than 4000 samples and 90,000 ASVs. The database, which is accessible through both a web-based interface (https://shiny.metapr2.org) and an R package, should prove very useful to all researchers working on protist diversity in a variety of systems.

pr2-primers: An 18S rRNA primer database for protists.

Metabarcoding of microbial eukaryotes (collectively known as protists) has developed tremendously in the last decade, almost solely relying on the 18S rRNA gene. As microbial eukaryotes are extremely diverse, many primers and primer pairs have been developed. To cover a relevant and representative fraction of the protist community in a given study system, an informed primer choice is necessary, as no primer pair can target all protists equally well. As such, a smart primer choice is very difficult even for experts and there are very few online resources available to list existing primers. We built a database listing 285 primers and 83 unique primer pairs that have been used for eukaryotic 18S rRNA gene metabarcoding. In silico performance of primer pairs was tested against two sequence databases: PR2 version 4.12.0 for eukaryotes and a subset of silva version 132 for bacteria and archaea. We developed an R-based web application enabling browsing of the database, visualization of the taxonomic distribution of the amplified sequences with the number of mismatches, and testing any user-defined primer or primer set (https://app.pr2-primers.org). Taxonomic specificity of primer pairs, amplicon size and location of mismatches can also be determined. We identified universal primer sets that matched the largest number of sequences and analysed the specificity of some primer sets designed to target certain groups. This tool enables guided primer choices that will help a wide range of researchers to include protists as part of their investigations.

Hemiselmis aquamarina sp. nov. (Cryptomonadales, Cryptophyceae), A Cryptophyte with A Novel Phycobiliprotein Type (Cr-PC 564).

Cryptophytes are a small group of photosynthetic biflagellate organisms distributed worldwide in fresh, brackish and marine waters. Although members of this class are easily distinguished from other groups, species identification is difficult and studies concerning their diversity are scarce. Two strains of an undescribed Hemiselmis species were isolated from the marine waters off Brazil and Japan. Analyses of morphology, phycobiliprotein spectral characterization, molecular phylogeny and ITS2 secondary structure comparisons were performed to assist the identification. The morphological features of Hemiselmis aquamarina sp. nov. matches that of other species from the same genus, but it has a new type of phycocyanin. Molecular phylogeny and ITS2 secondary structure support H. aquamarina as a distinct species. Furthermore, phylogenetic inferences indicate H. aquamarina as closely related to H. tepida, H. andersenii and H. rufescens. Currently, all Hemiselmis species have been described from the Northern Hemisphere and most from the subtropical region. H. aquamarina is the first species of this genus described from the South Atlantic.

Annual phytoplankton dynamics in coastal waters from Fildes Bay, Western Antarctic Peninsula.

Year-round reports of phytoplankton dynamics in the West Antarctic Peninsula are rare and mainly limited to microscopy and/or pigment-based studies. We analyzed the phytoplankton community from coastal waters of Fildes Bay in the West Antarctic Peninsula between January 2014 and 2015 using metabarcoding of the nuclear and plastidial 18/16S rRNA gene from both size-fractionated and flow cytometry sorted samples. Overall 14 classes of photosynthetic eukaryotes were present in our samples with the following dominating: Bacillariophyta (diatoms), Pelagophyceae and Dictyochophyceae for division Ochrophyta, Mamiellophyceae and Pyramimonadophyceae for division Chlorophyta, Haptophyta and Cryptophyta. Each metabarcoding approach yielded a different image of the phytoplankton community with for example Prymnesiophyceae more prevalent in plastidial metabarcodes and Mamiellophyceae in nuclear ones. Diatoms were dominant in the larger size fractions and during summer, while Prymnesiophyceae and Cryptophyceae were dominant in colder seasons. Pelagophyceae were particularly abundant towards the end of autumn (May). In addition of Micromonas polaris and Micromonas sp. clade B3, both previously reported in Arctic waters, we detected a new Micromonas 18S rRNA sequence signature, close to, but clearly distinct from M. polaris, which potentially represents a new clade specific of the Antarctic. These results highlight the need for complementary strategies as well as the importance of year-round monitoring for a comprehensive description of phytoplankton communities in Antarctic coastal waters.

No evidence of Phago-mixotropy in Micromonas polaris (Mamiellophyceae), the Dominant Picophytoplankton Species in the Arctic.

In the Arctic Ocean, the small green alga Micromonas polaris dominates picophytoplankton during the summer months but is also present in winter. It has been previously hypothesized to be phago-mixotrophic (capable of bacteria ingestion) based on laboratory and field experiments. Prey uptake was analyzed in several M. polaris strains isolated from different regions and depths of the Arctic Ocean and in Ochromonas triangulata, a known phago-mixotroph used as a control. Measuring ingestion of either fluorescent beads or fluorescently labeled bacteria by flow cytometry, we found no evidence of phago-mixotrophy in any M. polaris strain while O. triangulata was ingesting both beads and bacteria. In addition, in silico predictions revealed that members of the genus Micromonas lack a genetic signature of phagocytotic capacity.

Author Correction: Chloropicophyceae, a new class of picophytoplanktonic prasinophytes.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Mantoniella beaufortii and Mantoniella baffinensis sp. nov. (Mamiellales, Mamiellophyceae), two new green algal species from the high arctic1.

Members of the class Mamiellophyceae comprise species that can dominate picophytoplankton diversity in polar waters. Yet, polar species are often morphologically indistinguishable from temperate species, although clearly separated by molecular features. Here we examine four Mamiellophyceae strains from the Canadian Arctic. The 18S rRNA and Internal Transcribed Spacer 2 (ITS2) gene phylogeny place these strains within the family Mamiellaceae (Mamiellales, Mamiellophyceae) in two separate clades of the genus Mantoniella. ITS2 synapomorphies support their placement as two new species, Mantoniella beaufortii and Mantoniella baffinensis. Both species have round green cells with diameter between 3 and 5 µm, one long flagellum and a short flagellum (~1 µm) and are covered by spiderweb-like scales, making both species similar to other Mantoniella species. Morphologically, M. beaufortii and M. baffinensis are most similar to the cosmopolitan M. squamata with only minor differences in scale structure distinguishing them. Screening of global marine metabarcoding data sets indicates M. beaufortii has only been recorded in seawater and sea ice samples from the Arctic, while no environmental barcode matches M. baffinensis. Like other Mamiellophyceae genera that have distinct polar and temperate species, the polar distribution of these new species suggests they are cold or ice-adapted Mantoniella species.

Temporal and spatial dynamics of Bacteria, Archaea and protists in equatorial coastal waters.

Singapore, an equatorial island in South East Asia, is influenced by a bi-annual reversal of wind directions which defines two monsoon seasons. We characterized the dynamics of the microbial communities of Singapore coastal waters by collecting monthly samples between February 2017 and July 2018 at four sites located across two straits with different trophic status, and sequencing the V6-V8 region of the small sub-unit ribosomal RNA gene (rRNA gene) of Bacteria, Archaea, and Eukaryota. Johor Strait, which is subjected to wider environmental fluctuations from anthropogenic activities, presented a higher abundance of copiotrophic microbes, including Cellvibrionales and Rhodobacterales. The mesotrophic Singapore Strait, where the seasonal variability is caused by changes in the oceanographic conditions, harboured a higher proportion of typically marine microbe groups such as Synechococcales, Nitrosupumilales, SAR11, SAR86, Marine Group II Archaea and Radiolaria. In addition, we observed seasonal variability of the microbial communities in the Singapore Strait, which was possibly influenced by the alternating monsoon regime, while no seasonal pattern was detected in the Johor Strait.

Novel diversity within marine Mamiellophyceae (Chlorophyta) unveiled by metabarcoding.

Mamiellophyceae (unicellular green algae) are a key phytoplankton group in coastal waters. Although extensively studied over the last 20 years, the overall oceanic distribution of the major species/clades is still poorly known. To address this problem, we analyzed the 2014 Ocean Sampling Day (OSD) metabarcoding dataset providing sequences from the V4 hypervariable region of the 18S rRNA gene for 157 samples collected at 143 mostly coastal stations. Mamiellophyceae were found at nearly all OSD stations and represented 55% of the green microalgae (Chlorophyta) reads. We performed phylogenetic analyses of unique OSD metabarcodes (amplicon single variants, ASVs) and GenBank reference sequences from cultures and from the environment, focusing on the four most represented genera: Ostreococcus (45% of the Mamiellophyceae reads), Micromonas (34%), Bathycoccus (10%) and Mantoniella (8.7%). These analyses uncovered novel diversity within each genus except Bathycoccus. In Ostreococcus, a new clade (E) was the second most represented clade after Ostreococcus "lucimarinus". Micromonas could be separated into nine clades, exceeding the six species and candidate species already described. Finally, we found two new environmental clades within Mantoniella. Each Mamiellophyceae clade had a specific distribution in the OSD dataset suggesting that they are adapted to different ecological niches.

High contribution of Rhizaria (Radiolaria) to vertical export in the California Current Ecosystem revealed by DNA metabarcoding.

Passive sinking of particulate organic matter (POM) is the main mechanism through which the biological pump transports surface primary production to the ocean interior. However, the contribution and variability of different biological sources to vertical export is not fully understood. Here, we use DNA metabarcoding of the 18S rRNA gene and particle interceptor traps (PITs) to characterize the taxonomic composition of particles sinking out of the photic layer in the California Current Ecosystem (CCE), a productive system with high export potential. The PITs included formalin-fixed and 'live' traps to investigate eukaryotic communities involved in the export and remineralization of sinking particles. Sequences affiliated with Radiolaria dominated the eukaryotic assemblage in fixed traps (90%), with Dinophyta and Metazoa making minor contributions. The prominence of Radiolaria decreased drastically in live traps, possibly due to selective consumption by copepods, heterotrophic nanoflagellates, and phaeodarians that were heavily enriched in these traps. These patterns were consistent across the water masses surveyed extending from the coast to offshore, despite major differences in productivity and trophic structure of the epipelagic plankton community. Our findings identify Radiolaria as major actors in export fluxes in the CCE.

Rhythmicity of coastal marine picoeukaryotes, bacteria and archaea despite irregular environmental perturbations.

Seasonality in marine microorganisms has been classically observed in phytoplankton blooms, and more recently studied at the community level in prokaryotes, but rarely investigated at the scale of individual microbial taxa. Here we test if specific marine eukaryotic phytoplankton, bacterial and archaeal taxa display yearly rhythms at a coastal site impacted by irregular environmental perturbations. Our seven-year study in the Bay of Banyuls (North Western Mediterranean Sea) shows that despite some fluctuating environmental conditions, many microbial taxa displayed significant yearly rhythms. The robust rhythmicity was found in both autotrophs (picoeukaryotes and cyanobacteria) and heterotrophic prokaryotes. Sporadic meteorological events and irregular nutrient supplies did, however, trigger the appearance of less common non-rhythmic taxa. Among the environmental parameters that were measured, the main drivers of rhythmicity were temperature and day length. Seasonal autotrophs may thus be setting the pace for rhythmic heterotrophs. Similar environmental niches may be driving seasonality as well. The observed strong association between Micromonas and SAR11, which both need thiamine precursors for growth, could be a first indication that shared nutritional niches may explain some rhythmic patterns of co-occurrence.

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

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

A novel species of the marine cyanobacterium Acaryochloris with a unique pigment content and lifestyle.

All characterized members of the ubiquitous genus Acaryochloris share the unique property of containing large amounts of chlorophyll (Chl) d, a pigment exhibiting a red absorption maximum strongly shifted towards infrared compared to Chl a. Chl d is the major pigment in these organisms and is notably bound to antenna proteins structurally similar to those of Prochloron, Prochlorothrix and Prochlorococcus, the only three cyanobacteria known so far to contain mono- or divinyl-Chl a and b as major pigments and to lack phycobilisomes. Here, we describe RCC1774, a strain isolated from the foreshore near Roscoff (France). It is phylogenetically related to members of the Acaryochloris genus but completely lacks Chl d. Instead, it possesses monovinyl-Chl a and b at a b/a molar ratio of 0.16, similar to that in Prochloron and Prochlorothrix. It differs from the latter by the presence of phycocyanin and a vestigial allophycocyanin energetically coupled to photosystems. Genome sequencing confirmed the presence of phycobiliprotein and Chl b synthesis genes. Based on its phylogeny, ultrastructural characteristics and unique pigment suite, we describe RCC1774 as a novel species that we name Acaryochloris thomasi. Its very unusual pigment content compared to other Acaryochloris spp. is likely related to its specific lifestyle.

dinoref: A curated dinoflagellate (Dinophyceae) reference database for the 18S rRNA gene.

Dinoflagellates are a heterogeneous group of protists present in all aquatic ecosystems where they occupy various ecological niches. They play a major role as primary producers, but many species are mixotrophic or heterotrophic. Environmental metabarcoding based on high-throughput sequencing is increasingly applied to assess diversity and abundance of planktonic organisms, and reference databases are definitely needed to taxonomically assign the huge number of sequences. We provide an updated 18S rRNA reference database of dinoflagellates: dinoref. Sequences were downloaded from genbank and filtered based on stringent quality criteria. All sequences were taxonomically curated, classified taking into account classical morphotaxonomic studies and molecular phylogenies, and linked to a series of metadata. dinoref includes 1,671 sequences representing 149 genera and 422 species. The taxonomic assignation of 468 sequences was revised. The largest number of sequences belongs to Gonyaulacales and Suessiales that include toxic and symbiotic species. dinoref provides an opportunity to test the level of taxonomic resolution of different 18S barcode markers based on a large number of sequences and species. As an example, when only the V4 region is considered, 374 of the 422 species included in dinoref can still be unambiguously identified. Clustering the V4 sequences at 98% similarity, a threshold that is commonly applied in metabarcoding studies, resulted in a considerable underestimation of species diversity.

Small eukaryotic phytoplankton communities in tropical waters off Brazil are dominated by symbioses between Haptophyta and nitrogen-fixing cyanobacteria.

Symbioses between eukaryotic algae and nitrogen-fixing cyanobacteria have been recognized in recent years as a key source of new nitrogen in the oceans. We investigated the composition of the small photosynthetic eukaryote communities associated with nitrogen-fixing cyanobacteria in the Brazilian South Atlantic Bight using a combination of flow cytometry sorting and high throughput sequencing of two genes: the V4 region of 18S rRNA and nifH. Two distinct eukaryotic communities were often encountered, one dominated by the Mamiellophyceae Bathycoccus and Ostreococcus, and one dominated by a prymnesiophyte known to live in symbiosis with the UCYN-A1 nitrogen-fixing cyanobacterium. Among nifH sequences, those from UCYN-A1 were most abundant but three other UCYN-A clades (A2, A3, A4) were also found. Network analysis confirmed the relation between A1 and A2 clades and their hypothesized hosts and pointed out to the potential association between novel clade A4 with Braarudosphaera bigelowii, previously hypothesized to host A2.

Comparison of coastal phytoplankton composition estimated from the V4 and V9 regions of the 18S rRNA gene with a focus on photosynthetic groups and especially Chlorophyta.

We compared the composition of eukaryotic communities using two genetic markers (18S rRNA V4 and V9 regions) at 27 sites sampled during Ocean Sampling Day 2014, with a focus on photosynthetic groups and, more specifically green algae (Chlorophyta). Globally, the V4 and V9 regions of the 18S rRNA gene provided similar images of alpha diversity and ecological patterns. However, V9 provided 20% more OTUs built at 97% identity than V4. 34% of the genera were found with both markers and, of the remnant, 22% were found only with V4 and 44% only with V9. For photosynthetic groups, V4 and V9 performed equally well to describe global communities at different taxonomic levels from the division to the genus and provided similar Chlorophyta distribution patterns. However, at lower taxonomic level, the V9 dataset failed for example to describe the diversity of Dolichomastigales (Chlorophyta, Mamiellophyceae) emphasizing the lack of V9 sequences for this group and the importance of the reference database for metabarcode analysis. We conclude that in order to address questions regarding specific groups (e.g., a given genus), it is necessary to choose the marker based not only on the genetic divergence within this group but also on the existence of reference sequences in databases.

Chloropicophyceae, a new class of picophytoplanktonic prasinophytes.

Prasinophytes are a paraphyletic group of nine lineages of green microalgae that are currently classified either at the class or order level or as clades without formal taxonomic description. Prasinophyte clade VII comprises picoplanktonic algae that are important components of marine phytoplankton communities, particularly in moderately oligotrophic waters. Despite first being cultured in the 1960s, this clade has yet to be formally described. Previous phylogenetic analyses using the 18S rRNA gene divided prasinophyte clade VII into three lineages, termed A, B and C, the latter formed by a single species, Picocystis salinarum, that to date has only been found in saline lakes. Strains from lineages A and B cannot be distinguished by light microscopy and have very similar photosynthetic pigment profiles corresponding to the prasino-2A pigment group. We obtained phenotypic and genetic data on a large set of prasinophyte clade VII culture strains that allowed us to clarify the taxonomy of this important marine group. We describe two novel classes, the Picocystophyceae and the Chloropicophyceae, the latter containing two novel genera, Chloropicon and Chloroparvula, and eight new species of marine picoplanktonic green algae.

Revision of the Genus Micromonas Manton et Parke (Chlorophyta, Mamiellophyceae), of the Type Species M. pusilla (Butcher) Manton & Parke and of the Species M. commoda van Baren, Bachy and Worden and Description of Two New Species Based on the Genetic and Phenotypic Characterization of Cultured Isolates.

The green picoalgal genus Micromonas is broadly distributed in estuaries, coastal marine habitats and open oceans, from the equator to the poles. Phylogenetic, ecological and genomic analyses of culture strains and natural populations have suggested that this cosmopolitan genus is composed of several cryptic species corresponding to genetic lineages. We performed a detailed analysis of variations in morphology, pigment content, and sequences of the nuclear-encoded small-subunit rRNA gene and the second internal transcribed spacer (ITS2) from strains isolated worldwide. A new morphological feature of the genus, the presence of tip hairs at the extremity of the hair point, was discovered and subtle differences in hair point length were detected between clades. Clear non-homoplasious synapomorphies were identified in the small-subunit rRNA gene and ITS2 spacer sequences of five genetic lineages. These findings lead us to provide emended descriptions of the genus Micromonas, of the type species M. pusilla, and of the recently described species M. commoda, as well as to describe 2 new species, M. bravo and M. polaris. By clarifying the status of the genetic lineages identified within Micromonas, these formal descriptions will facilitate further interpretations of large-scale analyses investigating ecological trends in time and space for this widespread picoplankter.