Franck Lejzerowicz

Ultra-deep sequencing of foraminiferal microbarcodes unveils hidden richness of early monothalamous lineages in deep-sea sediments

Deep-sea floors represent one of the largest and most complex ecosystems on Earth but remain essentially unexplored. The vastness and remoteness of this ecosystem make deep-sea sampling difficult, hampering traditional taxonomic observations and diversity assessment. This problem is particularly true in the case of the deep-sea meiofauna, which largely comprises small-sized, fragile, and difficult-to-identify metazoans and protists. Here, we introduce an ultra-deep sequencing-based metagenetic approach to examine the richness of benthic foraminifera, a principal component of deep-sea meiofauna. We used Illumina sequencing technology to assess foraminiferal richness in 31 unsieved deep-sea sediment samples from five distinct oceanic regions. We sequenced an extremely short fragment (36 bases) of the small subunit ribosomal DNA hypervariable region 37f, which has been shown to accurately distinguish foraminiferal species. In total, we obtained 495,978 unique sequences that were grouped into 1,643 operational taxonomic units, of which about half (841) could be reliably assigned to foraminifera. The vast majority of the operational taxonomic units (nearly 90%) were either assigned to early (ancient) lineages of soft-walled, single-chambered (monothalamous) foraminifera or remained undetermined and yet possibly belong to unknown early lineages. Contrasting with the classical view of multichambered taxa dominating foraminiferal assemblages, our work reflects an unexpected diversity of monothalamous lineages that are as yet unknown using conventional micropaleontological observations. Although we can only speculate about their morphology, the immense richness of deep-sea phylotypes revealed by this study suggests that ultra-deep sequencing can improve understanding of deep-sea benthic diversity considered until now as unknowable based on a traditional taxonomic approach.

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.

Environmental monitoring through protist next-generation sequencing metabarcoding: assessing the impact of fish farming on benthic foraminifera communities

The measurement of species diversity represents a powerful tool for assessing the impacts of human activities on marine ecosystems. Traditionally, the impact of fish farming on the coastal environment is evaluated by monitoring the dynamics of macrobenthic infaunal populations. However, taxonomic sorting and morphology‐based identification of the macrobenthos demand highly trained specialists and are extremely time‐consuming and costly, making it unsuitable for large‐scale biomonitoring efforts involving numerous samples. Here, we propose to alleviate this laborious task by developing protist metabarcoding tools based on next‐generation sequencing (NGS) of environmental DNA and RNA extracted from sediment samples. In this study, we analysed the response of benthic foraminiferal communities to the variation of environmental gradients associated with salmon farms in Scotland. We investigated the foraminiferal diversity based on ribosomal minibarcode sequences generated by the Illumina NGS technology. We compared the molecular data with morphospecies counts and with environmental gradients, including distance to cages and redox used as a proxy for sediment oxygenation. Our study revealed high variations between foraminiferal communities collected in the vicinity of fish farms and at distant locations. We found evidence for species richness decrease in impacted sites, especially visible in the RNA data. We also detected some candidate bioindicator foraminiferal species. Based on this proof‐of‐concept study, we conclude that NGS metabarcoding using foraminifera and other protists has potential to become a new tool for surveying the impact of aquaculture and other industrial activities in the marine environment.

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.

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.

Protist metabarcoding and environmental biomonitoring: Time for change

High-throughput amplicon sequencing of environmental DNA and/or RNA proved to be a powerful tool to describe protist diversity. This new approach called also the metabarcoding has totally transformed our view of protist diversity, revealing a large number of novel lineages and expanding the range of protist phylogenetic diversity at almost every taxonomic level. However, until now the objectives of the vast majority of metabarcoding studies were purely academic. Practical applications of protist metabarcoding are surprisingly scarce, despite the fact that several groups of protists are commonly used as bioindicators of environmental impacts in freshwater or marine ecosystems. Here, we are reviewing studies that examine the ecological applications of metabarcoding for two groups of well-known protist bioindicators: diatoms and foraminifera. The results of these studies show that despite some biological and technical biases, molecular data quite faithfully reflect the morphology-based biotic indices and provide a similar assessment of ecosystem status. In view of these results, protist metabarcoding appears as a rapid and accurate tool for the evaluation of the quality of aquatic ecosystems. Hence, we plead for integration of protist metabarcoding in future biomonitoring projects as a complement of traditional methods and a source of new biosensors for environmental impact assessment.

Molecular evidence for widespread occurrence of Foraminifera in soils.

Environmental SSU rDNA-based surveys are contributing to the dramatic revision of eukaryotic high-level diversity and phylogeny as the number of sequence data increases. This ongoing revolution gives the opportunity to test for the presence of some eukaryotic taxa in environments where they have not been found using classical microscopic observations. Here, we test whether the foraminifera, a group of single-celled eukaryotes, considered generally as typical for the marine ecosystems are present in soil. We performed foraminiferal-specific nested PCR on 20 soil DNA samples collected in contrasted environments. Unexpectedly, we found that the majority of the samples contain foraminiferal SSU rDNA sequences. In total, we obtained 49 sequences from 17 localities. Phylogenetic analysis clusters them in four groups branching among the radiation of early foraminiferal lineages. Three of these groups also include sequences originated from previous freshwater surveys, suggesting that there were up to four independent colonization events of terrestrial and/or freshwater ecosystems by ancestral foraminifera. As shown by our data, foraminifera are a widespread and diverse component of soil microbial communities. Yet, identification of terrestrial foraminiferal species and understanding of their ecological role represent an exciting challenge for future research.

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.

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.

Predicting the Ecological Quality Status of Marine Environments from eDNA Metabarcoding Data Using Supervised Machine Learning

Monitoring biodiversity is essential to assess the impacts of increasing anthropogenic activities in marine environments. Traditionally, marine biomonitoring involves the sorting and morphological identification of benthic macro-invertebrates, which is time-consuming and taxonomic-expertise demanding. High-throughput amplicon sequencing of environmental DNA (eDNA metabarcoding) represents a promising alternative for benthic monitoring. However, an important fraction of eDNA sequences remains unassigned or belong to taxa of unknown ecology, which prevent their use for assessing the ecological quality status. Here, we show that supervised machine learning (SML) can be used to build robust predictive models for benthic monitoring, regardless of the taxonomic assignment of eDNA sequences. We tested three SML approaches to assess the environmental impact of marine aquaculture using benthic foraminifera eDNA, a group of unicellular eukaryotes known to be good bioindicators, as features to infer macro-invertebrates based biotic indices. We found similar ecological status as obtained from macro-invertebrates inventories. We argue that SML approaches could overcome and even bypass the cost and time-demanding morpho-taxonomic approaches in future biomonitoring.