High-throughput sequencing and morphology perform equally well for benthic monitoring of marine ecosystems

Environmental diversity surveys are crucial for the bioassessment of anthropogenic impacts on marine ecosystems. Traditional benthic monitoring relying on morphotaxonomic inventories of macrofaunal communities is expensive, time-consuming and expertise-demanding. High-throughput sequencing of environmental DNA barcodes (metabarcoding) offers an alternative to describe biological communities. However, whether the metabarcoding approach meets the quality standards of benthic monitoring remains to be tested. Here, we compared morphological and eDNA/RNA-based inventories of metazoans from samples collected at 10 stations around a fish farm in Scotland, including near-cage and distant zones. For each of 5 replicate samples per station, we sequenced the V4 region of the 18S rRNA gene using the Illumina technology. After filtering, we obtained 841,766 metazoan sequences clustered in 163 Operational Taxonomic Units (OTUs). We assigned the OTUs by combining local BLAST searches with phylogenetic analyses. We calculated two commonly used indices: the Infaunal Trophic Index and the AZTI Marine Biotic Index. We found that the molecular data faithfully reflect the morphology-based indices and provides an equivalent assessment of the impact associated with fish farms activities. We advocate that future benthic monitoring should integrate metabarcoding as a rapid and accurate tool for the evaluation of the quality of marine benthic ecosystems.

Marine protist diversity in European coastal waters and sediments as revealed by high-throughput sequencing

Although protists are critical components of marine ecosystems, they are still poorly characterized. Here we analysed the taxonomic diversity of planktonic and benthic protist communities collected in six distant European coastal sites. Environmental deoxyribonucleic acid (DNA) and ribonucleic acid (RNA) from three size fractions (pico-, nano- and micro/mesoplankton), as well as from dissolved DNA and surface sediments were used as templates for tag pyrosequencing of the V4 region of the 18S ribosomal DNA. Beta-diversity analyses split the protist community structure into three main clusters: picoplankton-nanoplankton-dissolved DNA, micro/mesoplankton and sediments. Within each cluster, protist communities from the same site and time clustered together, while communities from the same site but different seasons were unrelated. Both DNA and RNA-based surveys provided similar relative abundances for most class-level taxonomic groups. Yet, particular groups were overrepresented in one of the two templates, such as marine alveolates (MALV)-I and MALV-II that were much more abundant in DNA surveys. Overall, the groups displaying the highest relative contribution were Dinophyceae, Diatomea, Ciliophora and Acantharia. Also, well represented were Mamiellophyceae, Cryptomonadales, marine alveolates and marine stramenopiles in the picoplankton, and Monadofilosa and basal Fungi in sediments. Our extensive and systematic sequencing of geographically separated sites provides the most comprehensive molecular description of coastal marine protist diversity to date.

Environmental Monitoring: Inferring the Diatom Index from Next-Generation Sequencing Data

Diatoms are widely used as bioindicators for the assessment of water quality in rivers and streams. Classically, the diatom biotic indices are based on the relative abundance of morphologically identified species weighted by their autoecological value. Obtaining such indices is time-consuming, costly, and requires excellent taxonomic expertise, which is not always available. Here we tested the possibility to overcome these limitations using a next-generation sequencing (NGS) approach to identify and quantify diatoms found in environmental DNA and RNA samples. We analyzed 27 river sites in the Geneva area (Switzerland), in order to compare the values of the Swiss Diatom Index (DI-CH) computed either by microscopic quantification of diatom species or directly from NGS data. Despite gaps in the reference database and variations in relative abundance of analyzed species, the diatom index shows a significant correlation between morphological and molecular data indicating similar biological quality status for the majority of sites. This proof-of-concept study demonstrates the potential of the NGS approach for identification and quantification of diatoms in environmental samples, opening new avenues toward the routine application of genetic tools for bioassessment and biomonitoring of aquatic ecosystems.

Molecular barcoding of aquatic oligochaetes: implications for biomonitoring

Aquatic oligochaetes are well recognized bioindicators of quality of sediments and water in watercourses and lakes. However, the difficult taxonomic determination based on morphological features compromises their more common use in eco-diagnostic analyses. To overcome this limitation, we investigated molecular barcodes as identification tool for broad range of taxa of aquatic oligochaetes. We report 185 COI and 52 ITS2 rDNA sequences for specimens collected in Switzerland and belonging to the families Naididae, Lumbriculidae, Enchytraeidae and Lumbricidae. Phylogenetic analyses allowed distinguishing 41 lineages separated by more than 10 % divergence in COI sequences. The lineage distinction was confirmed by Automatic Barcode Gap Discovery (ABGD) method and by ITS2 data. Our results showed that morphological identification underestimates the oligochaete diversity. Only 26 of the lineages could be assigned to morphospecies, of which seven were sequenced for the first time. Several cryptic species were detected within common morphospecies. Many juvenile specimens that could not be assigned morphologically have found their home after genetic analysis. Our study showed that COI barcodes performed very well as species identifiers in aquatic oligochaetes. Their easy amplification and good taxonomic resolution might help promoting aquatic oligochaetes as bioindicators for next generation environmental DNA biomonitoring of aquatic ecosystems.

Transcriptome analyses to investigate symbiotic relationships between marine protists

Rhizaria are an important component of oceanic plankton communities worldwide. A number of species harbor eukaryotic microalgal symbionts, which are horizontally acquired in the environment at each generation. Although these photosymbioses are determinant for Rhizaria ability to thrive in oceanic ecosystems, the mechanisms for symbiotic interactions are unclear. Using high-throughput sequencing technology (i.e., 454), we generated large Expressed Sequence Tag (EST) datasets from four uncultured Rhizaria, an acantharian (Amphilonche elongata), two polycystines (Collozoum sp. and Spongosphaera streptacantha), and one phaeodarian (Aulacantha scolymantha). We assessed the main genetic features of the host/symbionts consortium (i.e., the holobiont) transcriptomes and found rRNA sequences affiliated to a wide range of bacteria and protists in all samples, suggesting that diverse microbial communities are associated with the holobionts. A particular focus was then carried out to search for genes potentially involved in symbiotic processes such as the presence of c-type lectins-coding genes, which are proteins that play a role in cell recognition among eukaryotes. Unigenes coding putative c-type lectin domains (CTLD) were found in the species bearing photosynthetic symbionts (A. elongata, Collozoum sp., and S. streptacantha) but not in the non-symbiotic one (A. scolymantha). More particularly, phylogenetic analyses group CTLDs from A. elongata and Collozoum sp. on a distinct branch from S. streptacantha CTLDs, which contained carbohydrate-binding motifs typically observed in other marine photosymbiosis. Our data suggest that similarly to other well-known marine photosymbiosis involving metazoans, the interactions of glycans with c-type lectins is likely involved in modulation of the host/symbiont specific recognition in Radiolaria.

Accurate multiplexing and filtering for high-throughput amplicon-sequencing

Tagging amplicons with tag sequences appended to PCR primers allow the multiplexing of numerous samples for high-throughput sequencing (HTS). This approach is routinely used in HTS-based diversity analyses, especially in microbial ecology and biomedical diagnostics. However, amplicon library preparation is subject to pervasive sample sequence cross-contaminations as a result of tag switching events referred to as mistagging. Here, we sequenced seven amplicon libraries prepared using various multiplexing designs in order to measure the magnitude of this phenomenon and its impact on diversity analyses. Up to 28.2% of the unique sequences correspond to undetectable (critical) mistags in single- or saturated double-tagging libraries. We show the advantage of multiplexing samples following Latin Square Designs in order to optimize the detection of mistags and maximize the information on their distribution across samples. We use this information in designs incorporating PCR replicates to filter the critical mistags and to recover the exact composition of mock community samples. Being parameter-free and data-driven, our approach can provide more accurate and reproducible HTS data sets, improving the reliability of their interpretations.

Molecular Phylogeny and Morphology of Leannia veloxifera n. gen. et sp. Unveils a New Lineage of Monothalamous Foraminifera

Monothalamous (single-chambered) foraminifera have long been considered as the "poor cousins" of multichambered species, which calcareous and agglutinated tests dominate in the fossil record. This view is currently changing with environmental DNA surveys showing that the monothalamids may be as diverse as hard-shelled foraminifera. Yet, the majority of numerous molecular lineages revealed by eDNA studies remain anonymous. Here, we describe a new monothalamous species and genus isolated from the sample of sea grass collected in Gulf of Eilat (Red Sea). This new species, named Leannia veloxifera, is characterized by a tiny ovoid theca (about 50-100 μm) composed of thin organic wall, with two opposite apertures. The examined individuals are multinucleated and show very active reticulopodial movement. Phylogenetic analyses of SSU rDNA, actin, and beta-tubulin (ß-tubulin) show that the species represents a novel lineage branching separately from other monothalamous foraminifera. Interestingly, the SSU rDNA sequence of the new species is very similar to an environmental foraminiferal sequence from Bahamas, suggesting that the novel lineage may represent a group of shallow-water tropical allogromiids, poorly studied until now.

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.

Dependence of interfacial film organization on lipid molecular structure

Combination of surface analytical techniques was employed to investigate the interfacial behavior of the two designed lipids-N-stearoylglycine (1) and its bulky neutral headgroup-containing derivative N-stearoylvaline ethyl ester (2)-at the air-solution interface and as transferred layers on different substrates. Formation of monolayers at the air-water interface was monitored on pure water and on aqueous solutions of different pH. Crystallization effects were visualized at pure water by recording the hystereses in the Langmuir-Blodgett (LB) isotherms and by transferring the layers onto mica, gold (111), and ITO (indium-tin oxide on glass) electrodes. Subphase pH affects the morphology and patch formation in monolayers of 1, as evidenced by BAM measurements. At pH 8.2, formation of well-ordered crystallites is observed, which upon compression elongate according to predominantly 1-D growth mechanism to form a dense layer of crystallites. This effect is not observed in monolayers of 2, whose headgroup is not protonated. The orientation of layers of 1 transferred to the solid supports is also pH dependent, and their stability can be related to formation of a hydrogen-bonded networks. AFM images of 1 exhibited platelets of multilayer phase. The IR spectra of the ITO substrates covered by 1 indicated formation of hydrogen bonds between the amide groups. The nature of the adsorption layer and its organization as a function of potential were studied in-depth by EC STM using Au(111) as the substrate. A model showing the arrangement of hydrogen bonds between adsorbed molecules is presented and related to the observed organization of the layer.

Ancient DNA sheds new light on the Svalbard foraminiferal fossil record of the last millennium

Recent palaeogenetic studies have demonstrated the occurrence of preserved ancient DNA (aDNA) in various types of fossilised material. Environmental aDNA sequences assigned to modern species have been recovered from marine sediments dating to the Pleistocene. However, the match between the aDNA and the fossil record still needs to be evaluated for the environmental DNA approaches to be fully exploited. Here, we focus on foraminifera in sediments up to one thousand years old retrieved from the Hornsund fjord (Svalbard). We compared the diversity of foraminiferal microfossil assemblages with the diversity of aDNA sequenced from subsurface sediment samples using both cloning and high-throughput sequencing (HTS). Our study shows that 57% of the species archived in the fossil record were also detected in the aDNA data. However, the relative abundance of aDNA sequence reads and fossil specimens differed considerably. We also found a limited match between the stratigraphic occurrence of some fossil species and their aDNA sequences, especially in the case of rare taxa. The aDNA data comprised a high proportion of non-fossilised monothalamous species, which are known to dominate in modern foraminiferal communities of the Svalbard region. Our results confirm the relevance of HTS for studying past micro-eukaryotic diversity and provide insight into its ability to reflect fossil assemblages. Palaeogenetic studies including aDNA analyses of non-fossilised groups expand the range of palaeoceanographical proxies and therefore may increase the accuracy of palaeoenvironmental reconstructions.

The Marine Microbial Eukaryote Transcriptome Sequencing Project (MMETSP): illuminating the functional diversity of eukaryotic life in the oceans through transcriptome sequencing

Current sampling of genomic sequence data from eukaryotes is relatively poor, biased, and inadequate to address important questions about their biology, evolution, and ecology; this Community Page describes a resource of 700 transcriptomes from marine microbial eukaryotes to help understand their role in the world's oceans.

Patterns of rare and abundant marine microbial eukaryotes

BACKGROUND

Biological communities are normally composed of a few abundant and many rare species. This pattern is particularly prominent in microbial communities, in which most constituent taxa are usually extremely rare. Although abundant and rare subcommunities may present intrinsic characteristics that could be crucial for understanding community dynamics and ecosystem functioning, microbiologists normally do not differentiate between them. Here, we investigate abundant and rare subcommunities of marine microbial eukaryotes, a crucial group of organisms that remains among the least-explored biodiversity components of the biosphere. We surveyed surface waters of six separate coastal locations in Europe, independently considering the picoplankton, nanoplankton, and microplankton/mesoplankton organismal size fractions.

RESULTS

Deep Illumina sequencing of the 18S rRNA indicated that the abundant regional community was mostly structured by organismal size fraction, whereas the rare regional community was mainly structured by geographic origin. However, some abundant and rare taxa presented similar biogeography, pointing to spatiotemporal structure in the rare microeukaryote biosphere. Abundant and rare subcommunities presented regular proportions across samples, indicating similar species-abundance distributions despite taxonomic compositional variation. Several taxa were abundant in one location and rare in other locations, suggesting large oscillations in abundance. The substantial amount of metabolically active lineages found in the rare biosphere suggests that this subcommunity constitutes a diversity reservoir that can respond rapidly to environmental change.

CONCLUSIONS

We propose that marine planktonic microeukaryote assemblages incorporate dynamic and metabolically active abundant and rare subcommunities, with contrasting structuring patterns but fairly regular proportions, across space and time.

Molecular evidence for Lessepsian invasion of soritids (larger symbiont bearing benthic foraminifera)

The Mediterranean Sea is considered as one of the hotspots of marine bioinvasions, largely due to the influx of tropical species migrating through the Suez Canal, so-called Lessepsian migrants. Several cases of Lessepsian migration have been documented recently, however, little is known about the ecological characteristics of the migrating species and their aptitude to colonize the new areas. This study focused on Red Sea soritids, larger symbiont-bearing benthic foraminifera (LBF) that are indicative of tropical and subtropical environments and were recently found in the Israeli coast of the Eastern Mediterranean. We combined molecular phylogenetic analyses of soritids and their algal symbionts as well as network analysis of Sorites orbiculus Forskål to compare populations from the Gulf of Elat (northern Red Sea) and from a known hotspot in Shikmona (northern Israel) that consists of a single population of S. orbiculus. Our phylogenetic analyses show that all specimens found in Shikmona are genetically identical to a population of S. orbiculus living on a similar shallow water pebbles habitat in the Gulf of Elat. Our analyses also show that the symbionts found in Shikmona and Elat soritids belong to the Symbiodinium clade F5, which is common in the Red Sea and also present in the Indian Ocean and Caribbean Sea. Our study therefore provides the first genetic and ecological evidences that indicate that modern population of soritids found on the Mediterranean coast of Israel is probably Lessepsian, and is less likely the descendant of a native ancient Mediterranean species.

Molecular evidence for β-tubulin neofunctionalization in Retaria (Foraminifera and radiolarians)

Foraminifera and radiolarians are closely related amoeboid protists (i.e., retarians) often characterized by their shells and pseudopodia. Previous studies hypothesized that the unusual "Type 2" β-tubulin (β2) is critically involved in forming helical filaments (HFs), a unique microtubule (MT) assembly/disassembly intermediate found in foraminiferan reticulopodia. Such noncanonical β-tubulin sequences have also been found in two radiolarian species and appear to be closely related to the foraminiferan β2. In this study, we report 119 new β-tubulin transcript sequences from six foraminiferans, four radiolarians, and a related non-retarian species. We found that foraminiferan and radiolarian β2-tubulins share some of the unusual substitutions in the structurally essential and usually conserved domains. In the β-tubulin phylogeny, retarian β2-tubulin forms a monophyletic clade, well separated from the canonical β-tubulin (β1) ubiquitous in eukaryotes. Furthermore, we found that foraminiferan and radiolarian β2-tubulin lineages were under positive selection, and used homology models for foraminiferan α- and β-tubulin hexamers to understand the structural effect of the positively selected substitutions. We suggest that the positively selected substitutions play key roles in the transition of MT to HF by altering the lateral and longitudinal interactions between α- and β-tubulin heterodimers. Our results indicate that the unusual β2-tubulin is a molecular synapomorphy of retarians, and the β-tubulin gene duplication occurred before the divergence of Foraminifera and radiolarians. The duplicates have likely been subjected to neofunctionalization responsible for the unique MT to HF assembly/disassembly dynamics, and/or other unknown physiological processes in retarian protists.

Ancient DNA complements microfossil record in deep-sea subsurface sediments

Deep-sea subsurface sediments are the most important archives of marine biodiversity. Until now, these archives were studied mainly using the microfossil record, disregarding large amounts of DNA accumulated on the deep-sea floor. Accessing ancient DNA (aDNA) molecules preserved down-core would offer unique insights into the history of marine biodiversity, including both fossilized and non-fossilized taxa. Here, we recover aDNA of eukaryotic origin across four cores collected at abyssal depths in the South Atlantic, in up to 32.5 thousand-year-old sediment layers. Our study focuses on Foraminifera and Radiolaria, two major groups of marine microfossils also comprising diverse non-fossilized taxa. We describe their assemblages in down-core sediment layers applying both micropalaeontological and environmental DNA sequencing approaches. Short fragments of the foraminiferal and radiolarian small subunit rRNA gene recovered from sedimentary DNA extracts provide evidence that eukaryotic aDNA is preserved in deep-sea sediments encompassing the last glacial maximum. Most aDNA were assigned to non-fossilized taxa that also dominate in molecular studies of modern environments. Our study reveals the potential of aDNA to better document the evolution of past marine ecosystems and opens new horizons for the development of deep-sea palaeogenomics.

The revised classification of eukaryotes

This revision of the classification of eukaryotes, which updates that of Adl et al. [J. Eukaryot. Microbiol. 52 (2005) 399], retains an emphasis on the protists and incorporates changes since 2005 that have resolved nodes and branches in phylogenetic trees. Whereas the previous revision was successful in re-introducing name stability to the classification, this revision provides a classification for lineages that were then still unresolved. The supergroups have withstood phylogenetic hypothesis testing with some modifications, but despite some progress, problematic nodes at the base of the eukaryotic tree still remain to be statistically resolved. Looking forward, subsequent transformations to our understanding of the diversity of life will be from the discovery of novel lineages in previously under-sampled areas and from environmental genomic information.