Wide Occurrence of SSU rDNA Intragenomic Polymorphism in Foraminifera and its Implications for Molecular Species Identification

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.

Quantitative assessment of reef foraminifera community from metabarcoding data

Describing living community compositions is essential to monitor ecosystems in a rapidly changing world, but it is challenging to produce fast and accurate depiction of ecosystems due to methodological limitations. Morphological methods provide absolute abundances with limited throughput, whereas metabarcoding provides relative abundances of genes that may not correctly represent living communities from environmental DNA assessed with morphological methods. However, it has the potential to deliver fast descriptions of living communities provided that it is interpreted with validated species-specific calibrations and reference databases. Here, we developed a quantitative approach to retrieve from metabarcoding data the assemblages of living large benthic foraminifera (LBF), photosymbiotic calcifying protists, from Indonesian coral reefs that are under increasing anthropogenic pressure. To depict the diversity, we calculated taxon-specific correction factors to reduce biological biases by comparing surface area, biovolume and calcite volume, and the number of mitochondrial gene copies in seven common LBF species. To validate the approach, we compared calibrated datasets of morphological communities from mock samples with bulk reef sediment; both sample types were metabarcoded. The calibration of the data significantly improved the estimations of genus relative abundance, with a difference of ±5% on average, allowing for comparison of past morphological datasets with future molecular ones. Our results also highlight the application of our quantitative approach to support reef monitoring operations by capturing fine-scale processes, such as seasonal and pollution-driven dynamics, that require high-throughput sampling treatment.

Comparative assessment of the intragenomic variations of dinoflagellate Tripos species through single-cell sequencing

Tripos is a large dinoflagellate genus widely distributed in the world's oceans. Morphology-based species identification is inconclusive due to high morphological intraspecific variability. Metabarcoding analysis has been demonstrated to be effective for species identification and tracking their spatiotemporal dynamics. However, accumulating evidence suggests high levels of intragenomic variations (IGVs) are common in many algae, leading to concerns about overinterpretation of molecular diversity in metabarcoding studies. In this project, we evaluated and compared IGVs in Tripos species by conducting the first high-throughput sequencing (HTS) of 18S rDNA V4 of Tripos single cells. High numbers of haplotypes (19-172) were identified in each of the 30 Tripos cells. Each cell contained one dominant haplotype with high relative abundance and many haplotypes with lower abundances. Thus, the presence of multiple minor haplotypes substantially overestimate the molecular diversity identified in metabarcoding analysis, which encompass not only interspecific and intraspecific diversities, but high levels of IGVs.

Assigning the unassigned: A signature-based classification of rDNA metabarcodes reveals new deep-sea diversity

Environmental DNA metabarcoding reveals a vast genetic diversity of marine eukaryotes. Yet, most of the metabarcoding data remain unassigned due to the paucity of reference databases. This is particularly true for the deep-sea meiofauna and eukaryotic microbiota, whose hidden diversity is largely unexplored. Here, we tackle this issue by using unique DNA signatures to classify unknown metabarcodes assigned to deep-sea foraminifera. We analyzed metabarcoding data obtained from 311 deep-sea sediment samples collected in the Clarion-Clipperton Fracture Zone, an area of potential polymetallic nodule exploitation in the Eastern Pacific Ocean. Using the signatures designed in the 37F hypervariable region of the 18S rRNA gene, we were able to classify 802 unassigned metabarcodes into 61 novel lineages, which have been placed in 27 phylogenetic clades. The comparison of new lineages with other foraminiferal datasets shows that most novel lineages are widely distributed in the deep sea. Five lineages are also present in the shallow-water datasets; however, phylogenetic analysis of these lineages separates deep-sea and shallow-water metabarcodes except in one case. While the signature-based classification does not solve the problem of gaps in reference databases, this taxonomy-free approach provides insight into the distribution and ecology of deep-sea species represented by unassigned metabarcodes, which could be useful in future applications of metabarcoding for environmental monitoring.

Unlocking foraminiferal genetic diversity on estuarine mudflats with eDNA metabarcoding

Environmental biomonitoring is a prerequisite for efficient evaluation and remediation of ecosystem degradation due to anthropogenic pressure or climate change. Estuaries are key habitats subject to multiple anthropogenic and natural stressors. Due to these multiple stressors, the detection of anthropogenic pressure is challenging. The fact that abundant natural stressors often lead to negative quality assessments has been coined the "estuarine quality paradox". To solve this issue, the application of molecular approaches with successful bioindicators like foraminifera is promising. However, sampling protocols, molecular procedures and data analyses need to be validated before such tools can be routinely applied. We conducted an environmental DNA survey of estuarine mudflats along the French Atlantic coast, using a metabarcoding approach targeting foraminifera. Our results demonstrate that estuarine environments have only a few active OTUs dominating the community composition and a large stock of dormant or propagule stages. This last genetic diversity components constitute an important reservoir, with different species which can potentially develop in response to the temporal variability of the multiple stressors. In fact, different OTUs were dominant in the studied estuaries. Our statistical model shows that the physical and chemical characteristics of the sediment and the climatic conditions explain only 43 % of the community composition variance. This suggests that other, less easily quantifiable factors, such as the history and use of the estuaries or the ecological drift could play an important role as well. Environmental DNA biomonitoring opens new perspectives to better characterize the genetic diversity in estuaries.

Invasion success of a Lessepsian symbiont-bearing foraminifera linked to high dispersal ability, preadaptation and suppression of sexual reproduction

Among the most successful Lessepsian invaders is the symbiont-bearing benthic foraminifera Amphistegina lobifera. In its newly conquered habitat, this prolific calcifier and ecosystem engineer is exposed to environmental conditions that exceed the range of its native habitat. To disentangle which processes facilitated the invasion success of A. lobifera into the Mediterranean Sea we analyzed a ~ 1400 bp sequence fragment covering the SSU and ITS gene markers to compare the populations from its native regions and along the invasion gradient. The genetic variability was studied at four levels: intra-genomic, population, regional and geographical. We observed that the invasion is not associated with genetic differentiation, but the invasive populations show a distinct suppression of intra-genomic variability among the multiple copies of the rRNA gene. A reduced genetic diversity compared to the Indopacific is observed already in the Red Sea populations and their high dispersal potential into the Mediterranean appears consistent with a bridgehead effect resulting from the postglacial expansion from the Indian Ocean into the Red Sea. We conclude that the genetic structure of the invasive populations reflects two processes: high dispersal ability of the Red Sea source population pre-adapted to Mediterranean conditions and a likely suppression of sexual reproduction in the invader. This discovery provides a new perspective on the cost of invasion in marine protists: The success of the invasive A. lobifera in the Mediterranean Sea comes at the cost of abandonment of sexual reproduction.

Pilot Study on the Prevalence of Entamoeba gingivalis in Austria-Detection of a New Genetic Variant

Entamoeba gingivalis is a parasitic protist that resides in the oral cavity. Although E. gingivalis has been frequently detected in individuals with periodontitis, its precise role in this context remains to be established, since E. gingivalis is also regularly found in healthy individuals. Sequence data on E. gingivalis are still scarce, with only a limited number of sequences available in public databases. In this study, a diagnostic PCR protocol was established in order to obtain a first impression on the prevalence of E. gingivalis in Austria and enable a differentiation of isolates by targeting the variable internal transcribed spacer regions. In total, 59 voluntary participants were screened for E. gingivalis and almost 50% of the participants were positive, with a significantly higher prevalence of participants with self-reported gingivitis. Moreover, in addition to the established subtypes ST1 and ST2, a potentially new subtype was found, designated ST3. 18S DNA sequencing and phylogenetic analyses clearly supported a separate position of ST3. Interestingly, subtype-specific PCRs revealed that, in contrast to ST2, ST3 only occurred in association with ST1. ST2 and ST1/ST3 were more often associated with gingivitis; however, more data will be necessary to corroborate this observation.

Single-Cell Genomics Reveals the Divergent Mitochondrial Genomes of Retaria (Foraminifera and Radiolaria)

Mitochondria originated from an ancient bacterial endosymbiont that underwent reductive evolution by gene loss and endosymbiont gene transfer to the nuclear genome. The diversity of mitochondrial genomes published to date has revealed that gene loss and transfer processes are ongoing in many lineages. Most well-studied eukaryotic lineages are represented in mitochondrial genome databases, except for the superphylum Retaria-the lineage comprising Foraminifera and Radiolaria. Using single-cell approaches, we determined two complete mitochondrial genomes of Foraminifera and two nearly complete mitochondrial genomes of radiolarians. We report the complete coding content of an additional 14 foram species. We show that foraminiferan and radiolarian mitochondrial genomes contain a nearly fully overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. In contrast to animals and fungi, many protists encode a diverse set of proteins on their mitochondrial genomes, including several ribosomal genes; however, some aerobic eukaryotic lineages (euglenids, myzozoans, and chlamydomonas-like algae) have reduced mitochondrial gene content and lack all ribosomal genes. Similar to these reduced outliers, we show that retarian mitochondrial genomes lack ribosomal protein and tRNA genes, contain truncated and divergent small and large rRNA genes, and contain only 14 or 15 protein-coding genes, including nad1, -3, -4, -4L, -5, and -7, cob, cox1, -2, and -3, and atp1, -6, and -9, with forams and radiolarians additionally carrying nad2 and nad6, respectively. In radiolarian mitogenomes, a noncanonical genetic code was identified in which all three stop codons encode amino acids. Collectively, these results add to our understanding of mitochondrial genome evolution and fill in one of the last major gaps in mitochondrial sequence databases. IMPORTANCE We present the reduced mitochondrial genomes of Retaria, the rhizarian lineage comprising the phyla Foraminifera and Radiolaria. By applying single-cell genomic approaches, we found that foraminiferan and radiolarian mitochondrial genomes contain an overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. An alternative genetic code was identified in radiolarian mitogenomes in which all three stop codons encode amino acids. Collectively, these results shed light on the divergent nature of the mitochondrial genomes from an ecologically important group, warranting further questions into the biological underpinnings of gene content variability and genetic code variation between mitochondrial genomes.

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.

Mitochondrial cytochrome c oxidase subunit I (COI) metabarcoding of Foraminifera communities using taxon-specific primers

Foraminifera are a species-rich phylum of rhizarian protists that are highly abundant in most marine environments. Molecular methods such as metabarcoding have revealed a high, yet undescribed diversity of Foraminifera. However, so far only one molecular marker, the 18S ribosomal RNA, was available for metabarcoding studies on Foraminifera. Primers that allow amplification of foraminiferal mitochondrial cytochrome oxidase I (COI) and identification of Foraminifera species were recently published. Here we test the performance of these primers for the amplification of whole foraminiferal communities, and compare their performance to that of the highly degenerate LerayXT primers, which amplify the same COI region in a wide range of eukaryotes. We applied metabarcoding to 48 samples taken along three transects spanning a North Sea beach in the Netherlands from dunes to the low tide level, and analysed both sediment samples and meiofauna samples, which contained taxa between 42 µm and 1 mm in body size obtained by decantation from sand samples. We used single-cell metabarcoding (Girard et al., 2022) to generate a COI reference library containing 32 species of Foraminifera, and used this to taxonomically annotate our community metabarcoding data. Our analyses show that the highly degenerate LerayXT primers do not amplify Foraminifera, while the Foraminifera primers are highly Foraminifera- specific, with about 90% of reads assigned to Foraminifera and amplifying taxa from all major groups, i.e., monothalamids, Globothalamea, and Tubothalamea. We identified 176 Foraminifera ASVs and found a change in Foraminifera community composition along the beach transects from high tide to low tide level, and a dominance of single-chambered monothalamid Foraminifera. Our results highlight that COI metabarcoding can be a powerful tool for assessing Foraminiferal communities.

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

Ribosomal RNA (rRNA) genes are known to be valuable markers for the barcoding of eukaryotic life and its phylogenetic classification at various taxonomic levels. The large-scale exploration of environmental microbial diversity through metabarcoding approaches have been focused mainly on the V4 and V9 regions of the 18S rRNA gene. The accurate interpretation of such environmental surveys is hampered by technical (e.g., PCR and sequencing errors) and biological biases (e.g., intra-genomic variability). Here we explored the intra-genomic diversity of Nassellaria and Spumellaria specimens (Radiolaria) by comparing Sanger sequencing with Illumina and Oxford Nanopore Technologies (MinION). Our analysis determined that intra-genomic variability of Nassellaria and Spumellaria is generally low, yet some Spumellaria specimens showed two different copies of the V4 with <97% similarity. From the different sequencing methods, Illumina showed the highest number of contaminations (i.e., environmental DNA, cross-contamination, tag-jumping), revealed by its high sequencing depth; and MinION showed the highest sequencing rate error (~14%). Yet the long reads produced by MinION (~2900 bp) allowed accurate phylogenetic reconstruction studies. These results highlight the requirement for a careful interpretation of Illumina based metabarcoding studies, in particular regarding low abundant amplicons, and open future perspectives towards full-length rDNA environmental metabarcoding surveys.

Coupling between Ribotypic and Phenotypic Traits of Protists across Life Cycle Stages and Temperatures

Relationships between ribotypic and phenotypic traits of protists across life cycle stages remain largely unknown. Herein, we used single cells of two soil and two marine ciliate species to examine phenotypic and ribotypic traits and their relationships across lag, log, plateau, cystic stages and temperatures. We found that Colpoda inflata and Colpoda steinii demonstrated allometric relationships between 18S ribosomal DNA (rDNA) copy number per cell (CNPC), cell volume (CV), and macronuclear volume across all life cycle stages. Integrating previously reported data of Euplotes vannus and Strombidium sulcatum indicated taxon-dependent rDNA CNPC-CV functions. Ciliate and prokaryote data analysis revealed that the rRNA CNPC followed a unified power-law function only if the rRNA-deficient resting cysts were not considered. Hence, a theoretical framework was proposed to estimate the relative quantity of resting cysts in the protistan populations with total cellular rDNA and rRNA copy numbers. Using rDNA CNPC was a better predictor of growth rate at a given temperature than rRNA CNPC and CV, suggesting replication of redundant rDNA operons as a key factor that slows cell division. Single-cell high-throughput sequencing and analysis after correcting sequencing errors revealed multiple rDNA and rRNA variants per cell. Both encystment and temperature affected the number of rDNA and rRNA variants in several cases. The divergence of rDNA and rRNA sequence in a single cell ranged from 1% to 10% depending on species. These findings have important implications for inferring cell-based biological traits (e.g., species richness, abundance and biomass, activity, and community structure) of protists using molecular approaches. IMPORTANCE Based on phenotypic traits, traditional surveys usually characterize organismal richness, abundance, biomass, and growth potential to describe diversity, organization, and function of protistan populations and communities. The rRNA gene (rDNA) and its transcripts have been widely used as molecular markers in ecological studies of protists. Nevertheless, the manner in which these molecules relate to cellular (organismal) and physiological traits remains poorly understood, which could lead to misinterpretations of protistan diversity and ecology. The current research highlights the dynamic nature of cellular rDNA and rRNA contents, which tightly couple with multiple phenotypic traits in ciliated protists. We demonstrate that quantity of resting cysts and maximum growth rate of a population can be theoretically estimated using ribotypic trait-based models. The intraindividual sequence polymorphisms of rDNA and rRNA can be influenced by encystment and temperature, which should be considered when interpreting species-level diversity and community structure of microbial eukaryotes.

First report of mitochondrial COI in foraminifera and implications for DNA barcoding

Foraminifera are a species-rich phylum of rhizarian protists that are highly abundant in many marine environments and play a major role in global carbon cycling. Species recognition in Foraminifera is mainly based on morphological characters and nuclear 18S ribosomal RNA barcoding. The 18S rRNA contains variable sequence regions that allow for the identification of most foraminiferal species. Still, some species show limited variability, while others contain high levels of intragenomic polymorphisms, thereby complicating species identification. The use of additional, easily obtainable molecular markers other than 18S rRNA will enable more detailed investigation of evolutionary history, population genetics and speciation in Foraminifera. Here we present the first mitochondrial cytochrome c oxidase subunit 1 (COI) gene sequences ("barcodes") of Foraminifera. We applied shotgun sequencing to single foraminiferal specimens, assembled COI, and developed primers that allow amplification of COI in a wide range of foraminiferal species. We obtained COI sequences of 49 specimens from 17 species from the orders Rotaliida and Miliolida. Phylogenetic analysis showed that the COI tree is largely congruent with previously published 18S rRNA phylogenies. Furthermore, species delimitation with ASAP and ABGD algorithms showed that foraminiferal species can be identified based on COI barcodes.

Didymium pseudonivicola: A new myxomycete from the austral Andes emerges from broad-scale morphological and molecular analyses of D. nivicola collections

A new nivicolous myxomycete is described as a result of a comprehensive study of Didymium nivicola collections from the entire range of its occurrence. Statistical analysis of 12 morphological characters, phylogenetic analyses of nuc 18S rDNA and elongation factor 1-alpha gene (EF1A), and a delimitation method (automatic barcode gap diversity) have been applied to corroborate the identity of the new species. A preliminary morphological analysis of D. nivicola revealed high variability of South American populations where four types of spore ornamentation were noted. However, results of molecular study and statistical analysis of morphological characters did not support recognition of these four forms but the distinction of two morphotypes. Consequently, two species have been recognized: D. nivicola and the newly proposed D. pseudonivicola. The new species can be distinguished from D. nivicola by distinctly larger and mostly plasmodiocarpic sporophores, which are scattered to gregarious, paler spores, and by the paler, more delicate and more elastic capillitium. Spore ornamentation of D. pseudonivicola is uniform and can be described as distinctly spiny (pilate under scanning electron microscope [SEM]), whereas those of D. nivicola is more variable, where spines (pilae under SEM) are delicate, distinct, or conspicuous. Additionally, whereas D. nivicola is a species distributed worldwide, D. pseudonivicola occurs only in the austral Andes of Argentina and Chile.

Empirical evidence for concerted evolution in the 18S rDNA region of the planktonic diatom genus Chaetoceros

Concerted evolution is a process of homogenisation of repetitive sequences within a genome through unequal crossing over and gene conversion. This homogenisation is never fully achieved because mutations always create new variants. Classically, concerted evolution has been detected as "noise" in electropherograms and these variants have been characterised through cloning and sequencing of subsamples of amplified products. However, this approach limits the number of detectable variants and provides no information about the abundance of each variant. In this study, we investigated concerted evolution by using environmental time-series metabarcoding data, single strain high-throughput sequencing (HTS) and a collection of Sanger reference barcode sequences. We used six species of the marine planktonic diatom genus Chaetoceros as study system. Abundance plots obtained from environmental metabarcoding and single strain HTS showed the presence of a haplotype far more abundant than all the others (the "dominant" haplotype) and identical to the reference sequences of that species obtained with Sanger sequencing. This distribution fitted best with Zipf's law among the rank abundance/ dominance models tested. Furthermore, in each strain 99% of reads showed a similarity of 99% with the dominant haplotype, confirming the efficiency of the homogenisation mechanism of concerted evolution. We also demonstrated that minor haplotypes found in the environmental samples are not only technical artefacts, but mostly intragenomic variation generated by incomplete homogenisation. Finally, we showed that concerted evolution can be visualised inferring phylogenetic networks from environmental data. In conclusion, our study provides an important contribution to the understanding of concerted evolution and to the interpretation of DNA barcoding and metabarcoding data based on multigene family markers.

Integrating morphology and metagenomics to understand taxonomic variability of Amphisorus (Foraminifera, Miliolida) from Western Australia and Indonesia

Foraminifera are a group of mostly marine protists with high taxonomic diversity. Species identification is often complex, as both morphological and molecular approaches can be challenging due to a lack of unique characters and reference sequences. An integrative approach combining state of the art morphological and molecular tools is therefore promising. In this study, we analysed large benthic Foraminifera of the genus Amphisorus from Western Australia and Indonesia. Based on previous findings on high morphological variability observed in the Soritidae and the discontinuous distribution of Amphisorus along the coast of western Australia, we expected to find multiple morphologically and genetically unique Amphisorus types. In order to gain detailed insights into the diversity of Amphisorus, we applied micro CT scanning and shotgun metagenomic sequencing. We identified four distinct morphotypes of Amphisorus, two each in Australia and Indonesia, and showed that each morphotype is a distinct genotype. Furthermore, metagenomics revealed the presence of three dinoflagellate symbiont clades. The most common symbiont was Fugacium Fr5, and we could show that its genotypes were mostly specific to Amphisorus morphotypes. Finally, we assembled the microbial taxa associated with the two Western Australian morphotypes, and analysed their microbial community composition. Even though each Amphisorus morphotype harboured distinct bacterial communities, sampling location had a stronger influence on bacterial community composition, and we infer that the prokaryotic community is primarily shaped by the microhabitat rather than host identity. The integrated approach combining analyses of host morphology and genetics, dinoflagellate symbionts, and associated microbes leads to the conclusion that we identified distinct, yet undescribed taxa of Amphisorus. We argue that the combination of morphological and molecular methods provides unprecedented insights into the diversity of foraminifera, which paves the way for a deeper understanding of their biodiversity, and facilitates future taxonomic and ecological work.

Anthropogenic impact on the historical phytoplankton community of Lake Constance reconstructed by multimarker analysis of sediment-core environmental DNA

During the 20th century, many lakes in the Northern Hemisphere were affected by increasing human population and urbanization along their shorelines and catchment, resulting in aquatic eutrophication. Ecosystem monitoring commenced only after the changes became apparent, precluding any examination of timing and dynamics of initial community change in the past and comparison of pre- and postimpact communities. Peri-Alpine Lake Constance (Germany) underwent a mid-century period of eutrophication followed by re-oligotrophication since the 1980s and is now experiencing warm temperatures. We extended the period for which monitoring data of indicator organisms exist by analysing historical environmental DNA (eDNA) from a sediment core dating back some 110 years. Using three metabarcoding markers-for microbial eukaryotes, diatoms and cyanobacteria-we revealed two major breakpoints of community change, in the 1930s and the mid-1990s. In our core, the latest response was exhibited by diatoms, which are classically used as palaeo-bioindicators for the trophic state of lakes. Following re-oligotrophication, overall diversity values reverted to similar ones of the early 20th century, but multivariate analysis indicated that the present community is substantially dissimilar. Community changes of all three groups were strongly correlated to phosphorus concentration changes, whereas significant relationships to temperature were only observed when we did not account for temporal autocorrelation. Our results indicate that each microbial group analysed exhibited a unique response, highlighting the particular strength of multimarker analysis of eDNA, which is not limited to organisms with visible remains and can therefore discover yet unknown responses and abiotic-biotic relationships.

Diversity hotspot and unique community structure of foraminifera in the world's deepest marine blue hole - Sansha Yongle Blue Hole

Marine blue holes are precious geological heritages with high scientific research values. Their physical and chemical characteristics are unique because of the steep-walled structure and isolated water column which create isolated ecosystems in geographically restricted areas. The Sansha Yongle Blue Hole (SYBH) is the world's deepest marine blue hole. Here, we generated the first DNA metabarcoding dataset from SYBH sediment focusing on foraminifera, a group of protists that have colonized various marine environments. We collected sediment samples from SYBH along a depth gradient to characterize the foraminiferal diversity and compared them with the foraminiferal diversity of the costal Jiaozhou Bay (JZB) and the abyssal Northwest Pacific Ocean (NWP). We amplified the SSU rDNA of foraminifera and sequenced them with high-throughput sequencing. The results showed that the foraminiferal assemblages in SYBH were vertically structured in response to the abiotic gradients and diversity was higher than in JZB and NWP. This study illustrates the capacity of foraminifera to colonize hostile environments and shows that blue holes are natural laboratories to explore physiological innovation associated with anoxia.

Species detection and delineation in the marine planktonic diatoms Chaetoceros and Bacteriastrum through metabarcoding: making biological sense of haplotype diversity

High-throughput sequencing (HTS) metabarcoding is commonly applied to assess phytoplankton diversity. Usually, haplotypes are grouped into operational taxonomic units (OTUs) through clustering, whereby the resulting number of OTUs depends on chosen similarity thresholds. We applied, instead, a phylogenetic approach to infer taxa among 18S rDNA V4-metabarcode haplotypes gathered from 48 time-series samples using the marine planktonic diatoms Chaetoceros and Bacteriastrum as test case. The 73 recovered taxa comprised both solitary haplotypes and polytomies, the latter composed each of a highly abundant, dominant haplotype and one to several minor, peripheral haplotypes. The solitary and dominant haplotypes usually matched reference sequences, enabling species assignation of taxa. We hypothesise that the super-abundance of reads in dominant haplotypes results from the homogenization effect of concerted evolution. Reads of populous peripheral haplotypes and dominant haplotypes show comparable distribution patterns over the sample dates, suggesting that they are part of the same population. Many taxa revealed marked seasonality, with closely related ones generally showing distinct periodicity, whereas others occur year-round. Phylogenies inferred from metabarcode haplotypes enable delineation of biologically meaningful taxa, whereas OTUs resulting from clustering algorithms often deviate markedly from such taxa.

The global-scale distributions of soil protists and their contributions to belowground systems

Protists are ubiquitous in soil, where they are key contributors to nutrient cycling and energy transfer. However, protists have received far less attention than other components of the soil microbiome. We used amplicon sequencing of soils from 180 locations across six continents to investigate the ecological preferences of protists and their functional contributions to belowground systems. We complemented these analyses with shotgun metagenomic sequencing of 46 soils to validate the identities of the more abundant protist lineages. We found that most soils are dominated by consumers, although parasites and phototrophs are particularly abundant in tropical and arid ecosystems, respectively. The best predictors of protist composition (primarily annual precipitation) are fundamentally distinct from those shaping bacterial and archaeal communities (namely, soil pH). Some protists and bacteria co-occur globally, highlighting the potential importance of these largely undescribed belowground interactions. Together, this study allowed us to identify the most abundant and ubiquitous protists living in soil, with our work providing a cross-ecosystem perspective on the factors structuring soil protist communities and their likely contributions to soil functioning.

Evidence of Intra-individual SSU Polymorphisms in Dark-spored Myxomycetes (Amoebozoa)

The nuclear small subunit rRNA gene (SSU or 18S) is a marker frequently used in phylogenetic and barcoding studies in Amoebozoa, including Myxomycetes. Despite its common usage and the confirmed existence of divergent copies of ribosomal genes in other protists, the potential presence of intra-individual SSU variability in Myxomycetes has never been studied before. Here we investigated the pattern of nucleotide polymorphism in the 5' end fragment of SSU by cloning and sequencing a total of 238 variants from eight specimens, each representing a species of the dark-spored orders Stemonitidales and Physarales. After excluding singletons, a relatively low SSU intra-individual variability was found but our data indicate that this might be a widely distributed phenomenon in Myxomycetes as all samples analyzed possessed various ribotypes. To determine if the occurrence of multiple SSU variants within a single specimen has a negative effect on the circumscription of species boundaries, we conducted phylogenetic analyses that revealed that clone variation may be detrimental for inferring phylogenetic relationships among some of the specimens analyzed. Despite that intra-individual variability should be assessed in additional taxa, our results indicate that special care should be taken for species identification when working with closely related species.

The future of biotic indices in the ecogenomic era: Integrating (e)DNA metabarcoding in biological assessment of aquatic ecosystems

The bioassessment of aquatic ecosystems is currently based on various biotic indices that use the occurrence and/or abundance of selected taxonomic groups to define ecological status. These conventional indices have some limitations, often related to difficulties in morphological identification of bioindicator taxa. Recent development of DNA barcoding and metabarcoding could potentially alleviate some of these limitations, by using DNA sequences instead of morphology to identify organisms and to characterize a given ecosystem. In this paper, we review the structure of conventional biotic indices, and we present the results of pilot metabarcoding studies using environmental DNA to infer biotic indices. We discuss the main advantages and pitfalls of metabarcoding approaches to assess parameters such as richness, abundance, taxonomic composition and species ecological values, to be used for calculation of biotic indices. We present some future developments to fully exploit the potential of metabarcoding data and improve the accuracy and precision of their analysis. We also propose some recommendations for the future integration of DNA metabarcoding to routine biomonitoring programs.

Current view on phylogeny within the genus Flabellula Schaeffer, 1926 (Amoebozoa: Leptomyxida)

The molecular phylogeny of Flabellula Schaeffer, 1926 has been updated by analysing 18S rRNA and actin gene sequences of 19 new strains collected and characterised by the authors over the past ten years. The genus Flabellula Schaeffer, 1926 (Amoebozoa: Leptomyxida) is a taxon in which species delineation based on morphological data by themselves is insufficient or even misleading. The description of two novel species, F. schaefferi n. sp. and F. sawyeri n. sp., is justified by the congruence of morphological data with 18S rRNA and actin gene sequence phylogenies, in-silico secondary structure prediction of the V2 region in the 18S rRNA, and by recognition of species-specific sequential motifs within this region.

Surface ocean metabarcoding confirms limited diversity in planktonic foraminifera but reveals unknown hyper-abundant lineages

Since the advent of DNA metabarcoding surveys, the planktonic realm is considered a treasure trove of diversity, inhabited by a small number of abundant taxa, and a hugely diverse and taxonomically uncharacterized consortium of rare species. Here we assess if the apparent underestimation of plankton diversity applies universally. We target planktonic foraminifera, a group of protists whose known morphological diversity is limited, taxonomically resolved and linked to ribosomal DNA barcodes. We generated a pyrosequencing dataset of ~100,000 partial 18S rRNA foraminiferal sequences from 32 size fractioned photic-zone plankton samples collected at 8 stations in the Indian and Atlantic Oceans during the Tara Oceans expedition (2009–2012). We identified 69 genetic types belonging to 41 morphotaxa in our metabarcoding dataset. The diversity saturated at local and regional scale as well as in the three size fractions and the two depths sampled indicating that the diversity of foraminifera is modest and finite. The large majority of the newly discovered lineages occur in the small size fraction, neglected by classical taxonomy. These unknown lineages dominate the bulk [>0.8 µm] size fraction, implying that a considerable part of the planktonic foraminifera community biomass has its origin in unknown lineages.

Assessing SSU rDNA Barcodes in Foraminifera: A Case Study using Bolivina quadrata

The Small Subunit Ribosomal RNA gene (SSU rDNA) is a widely used tool to reconstruct phylogenetic relationships among foraminiferal species. Recently, the highly variable regions of this gene have been proposed as DNA barcodes to identify foraminiferal species. However, the resolution of these barcodes has not been well established, yet. In this study, we evaluate four SSU rDNA hypervariable regions (37/f, 41/f, 43/e, and 45/e) as DNA barcodes to distinguish among species of the genus Bolivina, with particular emphasis on Bolivina quadrata for which ten new sequences (KY468817-KY468826) were obtained during this study. Our analyses show that a single SSU rDNA hypervariable sequence is insufficient to resolve all Bolivina species and that some regions (37/f and 41/f) are more useful than others (43/e and 45/e) to distinguish among closely related species. In addition, polymorphism analyses reveal a high degree of variability. In the context of barcoding studies, these results emphasize the need to assess the range of intraspecific variability of DNA barcodes prior to their application to identify foraminiferal species in environmental samples; our results also highlight the possibility that a longer SSU rDNA region might be required to distinguish among species belonging to the same taxonomic group (i.e. genus).

Ciliate diversity and distribution patterns in the sediments of a seamount and adjacent abyssal plains in the tropical Western Pacific Ocean

Background

Benthic ciliates and the environmental factors shaping their distribution are far from being completely understood. Likewise, deep-sea systems are amongst the least understood ecosystems on Earth. In this study, using high-throughput DNA sequencing, we investigated the diversity and community composition of benthic ciliates in different sediment layers of a seamount and an adjacent abyssal plain in the tropical Western Pacific Ocean with water depths ranging between 813 m and 4566 m. Statistical analyses were used to assess shifts in ciliate communities across vertical sediment gradients and water depth.

Results

Nine out of 12 ciliate classes were detected in the different sediment samples, with Litostomatea accounting for the most diverse group, followed by Plagiopylea and Oligohymenophorea. The novelty of ciliate genetic diversity was extremely high, with a mean similarity of 93.25% to previously described sequences. On a sediment depth gradient, ciliate community structure was more similar within the upper sediment layers (0-1 and 9-10 cm) compared to the lower sediment layers (19-20 and 29-30 cm) at each site. Some unknown ciliate taxa which were absent from the surface sediments were found in deeper sediments layers. On a water depth gradient, the proportion of unique OTUs was between 42.2% and 54.3%, and that of OTUs shared by all sites around 14%. However, alpha diversity of the different ciliate communities was relatively stable in the surface layers along the water depth gradient, and about 78% of the ciliate OTUs retrieved from the surface layer of the shallowest site were shared with the surface layers of sites deeper than 3800 m. Correlation analyses did not reveal any significant effects of measured environmental factors on ciliate community composition and structure.

Conclusions

We revealed an obvious variation in ciliate community along a sediment depth gradient in the seamount and the adjacent abyssal plain and showed that water depth is a less important factor shaping ciliate distribution in deep-sea sediments unlike observed for benthic ciliates in shallow seafloors. Additionally, an extremely high genetic novelty of ciliate diversity was found in these habitats, which points to a hot spot for the discovery of new ciliate species.

Electronic supplementary material

The online version of this article (10.1186/s12866-017-1103-6) contains supplementary material, which is available to authorized users.

DNA barcoding in autotrophic euglenids: evaluation of COI and 18s rDNA

Autotrophic euglenids (Euglenophyceae) are a common and abundant group of microbial eukaryotes in freshwater habitats. They have a limited number of features, which can be observed using light microscopy, thus species identification is often problematic. Establishing a barcode for this group is therefore an important step toward the molecular identification of autotrophic euglenids. Based on the literature, we selected verified species and used a plethora of available methods to validate two molecular markers: COI and 18S rDNA (the whole sequence and three fragments separately) as potential DNA barcodes. Analyses of the COI gene were performed based on the data set of 43 sequences (42 obtained in this study) representing 24 species and the COI gene was discarded as a DNA barcode mainly due to a lack of universal primer sites. For 18S rDNA analyses we used a data set containing 263 sequences belonging to 86 taxonomically verified species. We demonstrated that the whole 18S rDNA is too long to be a useful marker, but from the three shorter analyzed variable regions we recommend variable regions V2V3 and V4 of 18S rDNA as autotrophic euglenid barcodes due to their high efficiency (above 95% and 90%, respectively).

Neovahlkampfia nana n. sp. Reinforcing an Underrepresented Subclade of Tetramitia, Heterolobosea

The study provides robust genetic evidence that a newly isolated naked ameba with morphological and ultrastructural features indicative of Heterolobosea is a new species. Neovahlkampfia nana n. sp. associates with the yet underrepresented subclade of Tetramitia I. Considerable differences found in 18S rRNA gene sequences of individual molecular clones derived from DNA of five clonal cultures, using a low fidelity DNA polymerase, raised the issue of intragenomic sequence variation, a phenomenon that has not been previously studied in Heterolobosea. However, as proved using a higher fidelity DNA polymerase, the sequence variability observed was introduced by PCR mediated by the low fidelity polymerase and fixed by molecular cloning. This points to the potentially dubious validity of some current nominal species of Heterolobosea that differ from one another in just one or two base positions.

Contrasting evolutionary patterns of 28S and ITS rRNA genes reveal high intragenomic variation in Cephalenchus (Nematoda): Implications for species delimitation

Concerted evolution is often assumed to be the evolutionary force driving multi-family genes, including those from ribosomal DNA (rDNA) repeat, to complete homogenization within a species, although cases of non-concerted evolution have been also documented. In this study, sequence variation of 28S and ITS ribosomal RNA (rRNA) genes in the genus Cephalenchus is assessed at three different levels, intragenomic, intraspecific, and interspecific. The findings suggest that not all Cephalenchus species undergo concerted evolution. High levels of intraspecific polymorphism, mostly due to intragenomic variation, are found in Cephalenchus sp1 (BRA-01). Secondary structure analyses of both rRNA genes and across different species show a similar substitution pattern, including mostly compensatory (CBC) and semi-compensatory (SBC) base changes, thus suggesting the functionality of these rRNA copies despite the variation found in some species. This view is also supported by low sequence variation in the 5.8S gene in relation to the flanking ITS-1 and ITS-2 as well as by the existence of conserved motifs in the former gene. It is suggested that potential cross-fertilization in some Cephalenchus species, based on inspection of female reproductive system, might contribute to both intragenomic and intraspecific polymorphism of their rRNA genes. These results reinforce the potential implications of intragenomic and intraspecific genetic diversity on species delimitation, especially in biodiversity studies based solely on metagenetic approaches. Knowledge of sequence variation will be crucial for accurate species diversity estimation using molecular methods.

Accurate assessment of the impact of salmon farming on benthic sediment enrichment using foraminiferal metabarcoding

Assessing the environmental impact of salmon farms on benthic systems is traditionally undertaken using biotic indices derived from microscopic analyses of macrobenthic infaunal (MI) communities. In this study, we tested the applicability of using foraminiferal-specific high-throughput sequencing (HTS) metabarcoding for monitoring these habitats. Sediment samples and physico-chemical data were collected along an enrichment gradient radiating out from three Chinook salmon (Oncorhynchus tshawytscha) farms in New Zealand. HTS of environmental DNA and RNA (eDNA/eRNA) resulted in 1,875,300 sequences that clustered into 349 Operational Taxonomic Units. Strong correlations were observed among various biotic indices calculated from MI data and normalized fourth-root transformed HTS data. Correlations were stronger using eRNA compared to eDNA data. Quantile regression spline analyses identified 12 key foraminiferal taxa that have potential to be used as bioindicator species. This study demonstrates the huge potential for using this method for biomonitoring of fish-farming and other marine industrial activities.

Next-generation environmental diversity surveys of foraminifera: preparing the future

Foraminifera are commonly defined as marine testate protists, and their diversity is mainly assessed on the basis of the morphology of their agglutinated or mineralized tests. Diversity surveys based on environmental DNA (eDNA) have dramatically changed this view by revealing an unexpected diversity of naked and organic-walled lineages as well as detecting foraminiferal lineages in soil and freshwater environments. Moreover, single-cell analyses have allowed discrimination among genetically distinctive types within almost every described morphospecies. In view of these studies, the foraminiferal diversity appeared to be largely underestimated, but its accurate estimation was impeded by the low speed and coverage of a cloning-based eDNA approach. With the advent of high-throughput sequencing (HTS) technologies, these limitations disappeared in favor of exhaustive descriptions of foraminiferal diversity in numerous samples. Yet, the biases and errors identified in early HTS studies raised some questions about the accuracy of HTS data and their biological interpretation. Among the most controversial issues affecting the reliability of HTS diversity estimates are (1) the impact of technical and biological biases, (2) the sensitivity and specificity of taxonomic sequence assignment, (3) the ability to distinguish rare species, and (4) the quantitative interpretation of HTS data. Here, we document the lessons learned from previous HTS surveys and present the current advances and applications focusing on foraminiferal eDNA. We discuss the problems associated with HTS approaches and predict the future trends and avenues that hold promises for surveying foraminiferal diversity accurately and efficiently.