Valentin Vasselon

R-Syst::diatom: an open-access and curated barcode database for diatoms and freshwater monitoring

Diatoms are micro-algal indicators of freshwater pollution. Current standardized methodologies are based on microscopic determinations, which is time consuming and prone to identification uncertainties. The use of DNA-barcoding has been proposed as a way to avoid these flaws. Combining barcoding with next-generation sequencing enables collection of a large quantity of barcodes from natural samples. These barcodes are identified as certain diatom taxa by comparing the sequences to a reference barcoding library using algorithms. Proof of concept was recently demonstrated for synthetic and natural communities and underlined the importance of the quality of this reference library. We present an open-access and curated reference barcoding database for diatoms, called R-Syst::diatom, developed in the framework of R-Syst, the network of systematic supported by INRA (French National Institute for Agricultural Research), see http://www.rsyst.inra.fr/en. R-Syst::diatom links DNA-barcodes to their taxonomical identifications, and is dedicated to identify barcodes from natural samples. The data come from two sources, a culture collection of freshwater algae maintained in INRA in which new strains are regularly deposited and barcoded and from the NCBI (National Center for Biotechnology Information) nucleotide database. Two kinds of barcodes were chosen to support the database: 18S (18S ribosomal RNA) and rbcL (Ribulose-1,5-bisphosphate carboxylase/oxygenase), because of their efficiency. Data are curated using innovative (Declic) and classical bioinformatic tools (Blast, classical phylogenies) and up-to-date taxonomy (Catalogues and peer reviewed papers). Every 6 months R-Syst::diatom is updated. The database is available through the R-Syst microalgae website (http://www.rsyst.inra.fr/) and a platform dedicated to next-generation sequencing data analysis, virtual_BiodiversityL@b (https://galaxy-pgtp.pierroton.inra.fr/). We present here the content of the library regarding the number of barcodes and diatom taxa. In addition to these information, morphological features (e.g. biovolumes, chloroplasts…), life-forms (mobility, colony-type) or ecological features (taxa preferenda to pollution) are indicated in R-Syst::diatom. Database URL: http://www.rsyst.inra.fr/

Application of high-throughput sequencing (HTS) metabarcoding to diatom biomonitoring: do DNA extraction methods matter?

Current freshwater biomonitoring with diatoms is based on microscopic examination of the morphology of their silica skeleton. This standardized approach is time consuming and requires a high degree of taxonomic expertise. Metabarcoding combined with high-throughput sequencing (HTS) has great potential for next-generation biomonitoring applications but requires standardization. Molecular inventories are strongly influenced by the DNA extraction method used, but the effect of extraction protocols has not been tested to enable selection of the best DNA extraction method for HTS metabarcoding. We used 5 DNA extraction methods combining various types of cell lysis and DNA purification to extract DNA from 8 pure diatom cultures and 8 samples from streams and lakes with differing water quality. We compared the methods based on: 1) quality and purity of the extracted DNA, 2) community inventories obtained from HTS targeting the ribulose-1, 5-bisphosphate carboxylase (rbcL) barcode, and 3) similarity between molecular and microscopy-based inventories of community composition and the Specific Pollution-sensitivity Index [SPI]. A method based on GenElute™-LPA had higher extraction efficiency than the 4 commercial kits but had the highest polymerase chain reaction inhibition level. All 5 methods were efficient for HTS, and method did not affect operational taxonomic unit richness. We observed variations in the relative abundance of some taxa within Nitzschia, Amphora, Encyonema, Gomphonema, and Navicula between 2 of the 5 methods, but method did not affect global diatom community composition or SPI values. SPI values calculated from microscopy-based inventories and molecular inventories based on all 5 extraction methods were strongly correlated. For convenience purposes (high DNA quantity and low cost), we encourage standardization of HTS diatom biomonitoring based on the SA-Gen method.

Avoiding quantification bias in metabarcoding: Application of a cell biovolume correction factor in diatom molecular biomonitoring

Handling Editor: Andrew Mahon Abstract 1. In recent years, remarkable progress has been made in developing environmental DNA metabarcoding. However, its ability to quantify species relative abundance remains uncertain, limiting its application for biomonitoring. In diatoms, although the rbcL gene appears to be a suitable barcode for diatoms, providing relevant qualitative data to describe taxonomic composition, improvement of species quantification is still required. 2. Here, we hypothesized that rbcL copy number is correlated with diatom cell biovolume (as previously described for the 18S gene) and that a correction factor (CF) based on cell biovolume should be applied to improve taxa quantification. We carried out a laboratory experiment using pure cultures of eight diatom species with contrasted cell biovolumes in order to (1) verify the relationship between rbcL copy numbers (estimated by qPCR) and diatom cell biovolumes and (2) define a potential CF. In order to evaluate CF efficiency, five mock communities were created by mixing different amounts of DNA from the eight species, and were sequenced using HTS and targeting the same rbcL barcode. 3. As expected, the correction of DNA reads proportions by the CF improved the congruence between morphological and molecular inventories. Final validation of the CF was obtained on environmental samples (metabarcoding data from 80 benthic biofilms) for which the application of CF allowed differences between molecular and morphological water quality indices to be reduced by 47%. 4. Overall, our results highlight the usefulness of applying a CF factor, which is effective in reducing over-estimation of high biovolume species, correcting quantitative biases in diatom metabarcoding studies and improving final water quality assessment.

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.

Do we similarly assess diversity with microscopy and high-throughput sequencing? Case of microalgae in lakes

Diatoms are a species-diverse phylum of microalgae often presenting high biomass in aquatic habitats. This makes them excellent ecological indicators in rivers and lakes. They are routinely used to assess ecological quality of rivers and lakes using microscopy, which is time consuming. An alternative is to determine species in samples based on short DNA barcodes and high-throughput sequencing (HTS). Former studies showed that community structure and water quality assessments based on diatoms deliver similar results with both methods. But, none evaluated if diversities were assessed in the same way despite the importance of this ecological metric. Based on littoral benthic samplings carried out in 56 pristine alpine lakes, we compared different diversity indices measured with microscopy and metabarcoding. Each lake was sampled in three different places of its littoral. We showed that α (diversity measured in a single sampling site of a given lake) and ϒ (total diversity in a lake where three independent samples were considered) diversities obtained with HTS were higher than those obtained with microscopy. This may be explained by the capacity of HTS to detect morphologically cryptic species and to better detect rare taxa. On the other hand, β diversity obtained with HTS was smaller, which may be explained by the capacity of HTS to detect very rare species and free-floating extracellular DNA. Nevertheless, diversity indices obtained with both methodologies were well correlated each other. This study validates the possibility to assess diatom diversity with HTS in a comparable way to the classical microscopic analysis.

DNA metabarcoding and microscopic analyses of sea turtles biofilms: Complementary to understand turtle behavior

Sea turtles are distributed in tropical and subtropical seas worldwide. They play several ecological roles and are considered important indicators of the health of marine ecosystems. Studying epibiotic diatoms living on turtle shells suggestively has great potential in the study of turtle behavior because diatoms are always there. However, diatom identification at the species level is time consuming, requires well-trained specialists, and there is a high probability of finding new taxa growing on turtle shells, which makes identification tricky. An alternative approach based on DNA barcoding and high throughput sequencing (HTS), metabarcoding, has been developed in recent years to identify species at the community level by using a DNA reference library. The suitabilities of morphological and molecular approaches were compared. Diatom assemblages were sampled from seven juvenile green turtles (Chelonia mydas) from Mayotte Island, France. The structures of the epibiotic diatom assemblages differed between both approaches. This resulted in different clustering of the turtles based on their diatom communities. Metabarcoding allowed better discrimination between turtles based on their epibiotic diatom assemblages and put into evidence the presence of a cryptic diatom diversity. Microscopy, for its part, provided more ecological information of sea turtles based on historical bibliographical data and the abundances of ecological guilds of the diatom species present in the samples. This study shows the complementary nature of these two methods for studying turtle behavior.

Boosting DNA metabarcoding for biomonitoring with phylogenetic estimation of operational taxonomic units’ ecological profiles

DNA metabarcoding has been introduced as a revolutionary way to identify organisms and monitor ecosystems. However, the potential of this approach for biomonitoring remains partially unfulfilled because a significant part of the sampled DNA cannot be affiliated to species due to incomplete reference libraries. Thus, biotic indices, which are based on the estimated abundances of species in a community and their ecological profiles, can be inaccurate. We propose to compute biotic indices using phylogenetic imputation of operational taxonomic units (OTUs’) ecological profiles (OTU‐PITI approach). First, OTUs sequences are inserted within a reference phylogeny. Second, OTUs’ ecological profiles are estimated on the basis of their phylogenetic relationships with reference species whose ecology is known. Based on these ecological profiles, biotic indices can be computed using all available OTUs. Using freshwater diatoms as a case study, we show that short DNA barcodes can be placed accurately within a phylogeny and their ecological preferences estimated with a satisfactory level of precision. In the light of these results, we tested the approach with a data set of 139 environmental samples of benthic river diatoms for which the same biotic index (specific sensitivity index) was calculated using (a) traditional microscopy, (b) OTUs with taxonomic assignment approach, (c) OTUs with phylogenetic estimation of ecological profiles (OTU‐PITI) and (d) OTU with taxonomic assignment completed by the phylogenetic approach (OTU‐PITI) for unclassified OTUs. Using traditional microscopy as a reference, we found that the combination of the OTUs’ taxonomic assignment completed by the phylogenetic method performed satisfactorily and substantially better than the other methods tested.