Intra-genomic rRNA gene variability of Nassellaria and Spumellaria (Rhizaria, Radiolaria) assessed by Sanger, MinION and Illumina sequencing

Harmonized coexistence of intragenomic variations in diatom Skeletonema strains

Metabarcoding analysis has been demonstrated to be an effective technology for monitoring diversity and dynamics of phytoplankton including Skeletonema species. Although molecular diversity uncovered in metabarcoding projects has generally been interpreted as sum of interspecies diversity and intraspecies diversity, accumulating evidence suggests that it also harbors unprecedentedly high levels of intra-genomic variations (IGVs). As up to thousands of amplicon sequence variants (ASVs) identified in a typical metabarcoding project can be annotated to be Skeletonema species, we hypothesize that substantial portions of these ASVs are contributed by IGVs. Here, the nature of IGVs in Skeletonema species was quantitatively analyzed by carrying out single-strain metabarcoding analysis of 18S rDNA V4 in 49 strains belonging to seven Skeletonema species. Results showed that each Skeletonema strain harbored a high level of IGVs as expected. While many Skeletonema strains each contained one dominant ASV and a substantial number of ASVs displaying much lower relative abundance, other Skeletonema strains each contained multiple ASVs with comparable or nearly equally abundances. Thus the co-existence of multiple dominant ASVs in a single cell indicated a tug-of-war of these variants in evolution, which may eventually result in harmonized coexistence of multiple dominant ASVs. A total of nine dominant ASVs and 652 non-dominant ASVs were found in 49 strains of seven Skeletonema species, indicating rich interspecies and intraspecies variations, and complex evolution of IGVs in genus of Skeletonema. The results confirmed that the extensive degree of IGVs was the main contributor to the high molecular diversity revealed by metabarcoding analysis. This study highlights the importance of quantitative characterization of IGVs in Skeletonema species for accurate interpretation of species diversity in metabarcoding analysis.

Environmental selection and advective transport shape the distribution of two cyst-forming Acantharia clades in the Canadian Arctic

Anthropogenic induced climate perturbations are seen in changes in oceanic circulation patterns, and Arctic water masses defined by salinity are vulnerable to change. Biogeography of marine microbial eukaryotes is expected to be impacted by changes in local environmental conditions and advective processes, but tracking the extent of plankton distribution requires understanding routes for both active and passive tracers. To identify such tracers, we focused on samples collected in the western (Canada Basin) and eastern (Nares Strait); extremes of the Canadian High Arctic that are connected by an east flowing current north of Canada. Sequencing of the V4 region of 18S rRNA revealed that Acantharia, a taxonomically and functionally diverse group of large planktonic protists, were particularly common. Arctic acantharians in our study were dominated by two clades belonging to cyst-forming groups. The distribution of one clade suggested successful advective transport from the Pacific sourced water in the Beaufort Gyre to southern Nares Strait, with cells transported along the northern shelf of the Canadian Arctic. A second clade appeared to be a resident taxon of the Canada Basin whose distribution correlated to local environmental conditions, and detection in deeper samples would be consistent with swarmer formation enabling reestablishment the following year.

Identification of phytoplankton isolates from the eastern Canadian waters using long-read sequencing

As important primary producers and key components of biogeochemical processes, phytoplankton communities are highly diverse and continually adapt to changes in the environment, impacting the entire marine ecosystem. Consequently, it remains important to isolate, culture and properly describe new phytoplankton strains to provide relevant model organisms for laboratory research and accurate reference sequences for identification. Here, 73 phytoplankton strains from the eastern Canadian waters were isolated and genetically characterized using a long rRNA fragment (~4000 bp) covering the 18S, ITS and 28S rRNA regions generated with long-read sequencing technology. While most strains (66%) were accurately identified using the partial 18S rRNA gene (~1200 bp-V4-V9), this study broadens the current 28S rRNA reference database by incorporating 41 distinct strains including 30 unique species. In addition, previously unpublished long-read reference sequences were generated for a few genera, including Picochlorum and Droopiella. Species that were previously poorly described in the eastern Canadian waters were also genetically characterized, including a Chaetoceros similis strain. Overall, this work expands the total number of long-read references, an essential resource for accurate identification of phytoplankton and environmental sequences, thereby advancing the taxonomic resolution that could lead to a better understanding of environmental microbial diversity.

The global genetic diversity of planktonic foraminifera reveals the structure of cryptic speciation in plankton

The nature and extent of diversity in the plankton has fascinated scientists for over a century. Initially, the discovery of many new species in the remarkably uniform and unstructured pelagic environment appeared to challenge the concept of ecological niches. Later, it became obvious that only a fraction of plankton diversity had been formally described, because plankton assemblages are dominated by understudied eukaryotic lineages with small size that lack clearly distinguishable morphological features. The high diversity of the plankton has been confirmed by comprehensive metabarcoding surveys, but interpretation of the underlying molecular taxonomies is hindered by insufficient integration of genetic diversity with morphological taxonomy and ecological observations. Here we use planktonic foraminifera as a study model and reveal the full extent of their genetic diversity and investigate geographical and ecological patterns in their distribution. To this end, we assembled a global data set of ~7600 ribosomal DNA sequences obtained from morphologically characterised individual foraminifera, established a robust molecular taxonomic framework for the observed diversity, and used it to query a global metabarcoding data set covering ~1700 samples with ~2.48 billion reads. This allowed us to extract and assign 1 million reads, enabling characterisation of the structure of the genetic diversity of the group across ~1100 oceanic stations worldwide. Our sampling revealed the existence of, at most, 94 distinct molecular operational taxonomic units (MOTUs) at a level of divergence indicative of biological species. The genetic diversity only doubles the number of formally described species identified by morphological features. Furthermore, we observed that the allocation of genetic diversity to morphospecies is uneven. Only 16 morphospecies disguise evolutionarily significant genetic diversity, and the proportion of morphospecies that show genetic diversity increases poleward. Finally, we observe that MOTUs have a narrower geographic distribution than morphospecies and that in some cases the MOTUs belonging to the same morphospecies (cryptic species) have different environmental preferences. Overall, our analysis reveals that even in the light of global genetic sampling, planktonic foraminifera diversity is modest and finite. However, the extent and structure of the cryptic diversity reveals that genetic diversification is decoupled from morphological diversification, hinting at different mechanisms acting at different levels of divergence.

Intragenomic diversity of the V9 hypervariable domain in eukaryotes has little effect on metabarcoding

Metabarcoding revolutionized our understanding of diversity and ecology of microorganisms in different habitats. However, it is also associated with several inherent biases, one of which is associated with intragenomic diversity of a molecular barcode. Here, we compare intragenomic variability of the V9 region of the 18S rRNA gene in 19 eukaryotic phyla abundant in marine plankton. The level of intragenomic variability is comparable across all the phyla, and in most genomes and transcriptomes one V9 sequence and one OTU is predominant. However, most of the variability observed at the barcode level is probably caused by sequencing errors and is mitigated by using a denoising tool, DADA2. The SWARM algorithm commonly used in metabarcoding studies is not optimal for collapsing genuine and erroneous sequences into a single OTU, leading to an overestimation of diversity in metabarcoding data. For an unknown reason, SWARM inflates diversity of eupelagonemids more than that of other eukaryotes.

Macroevolutionary patterns in intragenomic rDNA variability among planktonic foraminifera

Ribosomal intragenomic variability in prokaryotes and eukaryotes is a genomic feature commonly studied for its inflationary impact on molecular diversity assessments. However, the evolutionary mechanisms and distribution of this phenomenon within a microbial group are rarely explored. Here, we investigate the intragenomic variability in 33 species of planktonic foraminifera, calcifying marine protists, by inspecting 2,403 partial SSU sequences obtained from single-cell clone libraries. Our analyses show that polymorphisms are common among planktonic foraminifera species, but the number of polymorphic sites significantly differs among clades. With our molecular simulations, we could assess that most of these mutations are located in paired regions that do not affect the secondary structure of the SSU fragment. Finally, by mapping the number of polymorphic sites on the phylogeny of the clades, we were able to discuss the evolution and potential sources of intragenomic variability in planktonic foraminifera, linking this trait to the distinctive nuclear and genomic dynamics of this microbial group.

Does filter pore size introduce bias in DNA sequence-based plankton community studies?

The cell size of microbial eukaryotic plankton normally ranges from 0.2 to 200 μm. During the past decade, high-throughput sequencing of DNA has been revolutionizing their study on an unprecedented scale. Nonetheless, it is currently unclear whether we can accurately, effectively, and quantitatively depict the microbial eukaryotic plankton community using size-fractionated filtration combined with environmental DNA (eDNA) molecular methods. Here we assessed the microbial eukaryotic plankton communities with two filtering strategies from two subtropical reservoirs, that is one-step filtration (0.2–200 μm) and size-fractionated filtration (0.2–3 and 3–200 μm). The difference of 18S rRNA gene copy abundance between the two filtering treatments was less than 50% of the 0.2–200 μm microbial eukaryotic community for 95% of the total samples. Although the microbial eukaryotic plankton communities within the 0.2–200 μm and the 0.2–3 and 3–200 μm size fractions had approximately identical 18S rRNA gene copies, there were significant differences in their community composition. Furthermore, our results demonstrate that the systemic bias introduced by size-fractionation filtration has more influence on unique OTUs than shared OTUs, and the significant differences in abundance between the two eukaryotic plankton communities largely occurred in low-abundance OTUs in specific seasons. This work provides new insights into the use of size-fractionation in molecular studies of microbial eukaryotes populating the plankton.