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

Genetic Markers for Metabarcoding of Freshwater Microalgae: Review

The metabarcoding methods for studying the diversity of freshwater microalgae and routine biomonitoring are actively used in modern research. A lot of experience has been accumulated already, and many methodological questions have been solved (such as the influence of the methods and time of sample conservation, DNA extraction and bioinformatical processing). The reproducibility of the method has been tested and confirmed. However, one of the main problems-choosing a genetic marker for the study-still lacks a clear answer. We analyzed 70 publications and found out that studies on eukaryotic freshwater microalgae use 12 markers (different nuclear regions 18S and ITS and plastids rbcL, 23S and 16S). Each marker has its peculiarities; they amplify differently and have various levels of efficiency (variability) in different groups of algae. The V4 and V9 18S and rbcL regions are used most often. We concentrated especially on the studies that compare the results of using different markers and microscopy. We summarize the data on the primers for each region and on how the choice of a marker affects the taxonomic composition of a community.

Environmental filtering and mass effect are two important processes driving lake benthic diatoms: Results of a DNA metabarcoding study in a large lake

Environmental filtering is often found to dominate assembly rules in diatoms. These microalgae are diverse, especially at subspecies level, and tend to exhibit a niche phylogenetic conservatism. Therefore, other rules, such as competition or mass effects, should be detectable when environmental gradients and dispersal barriers are limited. We used metabarcoding to analyse benthic littoral diatom communities in 153 sites in a large lake (Geneva) exhibiting weak geographical barriers and weak environmental gradients outside river estuaries. We assessed assembly rules using variance partitioning, phylogenetic and source tracking analyses. No phylogenetic over-dispersion of communities, indicative of exclusive competition, was detected. Instead, we found these communities to be phylogenetically over-clustered, indicating environmental filtering, which was even stronger near river estuaries where environmental gradients are stronger. Finally, using a Bayesian method (SourceTracker), we found that rivers flowing into the lake bring communities that settle, especially in sites close to estuaries. Rivers with the highest discharges are primarily responsible for immigration, explaining 27% of lake composition. Therefore, despite favourable conditions to observe other rules, our results support that diatom communities are prominently assembled by environmental filtering and immigration processes, in particular from rivers. However, this does not exclude that other assembly rules may be at play at a finer spatial, temporal and/or phylogenetic scale.

More than just hitchhikers: a survey of bacterial communities associated with diatoms originating from sea turtles

Diatoms and bacteria are known for being the first colonizers of submerged surfaces including the skin of marine reptiles. Sea turtle carapace and skin harbor diverse prokaryotic and eukaryotic microbes, including several epizoic diatoms. However, the importance of diatom-bacteria associations is hardly investigated in biofilms associated with animal hosts. This study provides an inventory of diatoms, bacteria and diatom-associated bacteria originating from loggerhead sea turtles using both metabarcoding and culturing approaches. Amplicon sequencing of the carapace and skin samples chloroplast gene rbcL and 16S rRNA gene detected, in total, 634 diatom amplicon sequence variants (ASVs) and 3661 bacterial ASVs, indicating high diversity. Cultures of putative epizoic and non-epizoic diatoms contained 458 bacterial ASVs and their bacterial assemblages reflected those of their host. Diatom strains allowed for enrichment and isolation of bacterial families rarely observed on turtles, such as Marinobacteraceae, Alteromonadaceae and Alcanivoracaceae. When accounting for phylogenetic relationships between bacterial ASVs, we observed that related diatom genera might retain similar microbial taxa in culture, regardless of the turtle's skin or carapace source. These data provide deeper insights into the sea turtle-associated microbial communities, and reveal the potential of epizoic biofilms as a source of novel microbes and possibly important diatom-bacteria associations.

Red, Gold and Green: Microbial Contribution of Rhodophyta and Other Algae to Green Turtle (Chelonia mydas) Gut Microbiome

The fitness of the endangered green sea turtle (Chelonia mydas) may be strongly affected by its gut microbiome, as microbes play important roles in host nutrition and health. This study aimed at establishing environmental microbial baselines that can be used to assess turtle health under altered future conditions. We characterized the microbiome associated with the gastrointestinal tract of green turtles from Guinea Bissau in different life stages and associated with their food items, using 16S rRNA metabarcoding. We found that the most abundant (% relative abundance) bacterial phyla across the gastrointestinal sections were Proteobacteria (68.1 ± 13.9% “amplicon sequence variants”, ASVs), Bacteroidetes (15.1 ± 10.1%) and Firmicutes (14.7 ± 21.7%). Additionally, we found the presence of two red algae bacterial indicator ASVs (the Alphaproteobacteria Brucella pinnipedialis with 75 ± 0% and a Gammaproteobacteria identified as methanotrophic endosymbiont of Bathymodiolus, with <1%) in cloacal compartments, along with six bacterial ASVs shared only between cloacal and local environmental red algae samples. We corroborate previous results demonstrating that green turtles fed on red algae (but, to a lower extent, also seagrass and brown algae), thus, acquiring microbial components that potentially aid them digest these food items. This study is a foundation for better understanding the microbial composition of sea turtle digestive tracts.

Cultivating epizoic diatoms provides insights into the evolution and ecology of both epibionts and hosts

Our understanding of the importance of microbiomes on large aquatic animals—such as whales, sea turtles and manatees—has advanced considerably in recent years. The latest observations indicate that epibiotic diatom communities constitute diverse, polyphyletic, and compositionally stable assemblages that include both putatively obligate epizoic and generalist species. Here, we outline a successful approach to culture putatively obligate epizoic diatoms without their hosts. That some taxa can be cultured independently from their epizoic habitat raises several questions about the nature of the interaction between these animals and their epibionts. This insight allows us to propose further applications and research avenues in this growing area of study. Analyzing the DNA sequences of these cultured strains, we found that several unique diatom taxa have evolved independently to occupy epibiotic habitats. We created a library of reference sequence data for use in metabarcoding surveys of sea turtle and manatee microbiomes that will further facilitate the use of environmental DNA for studying host specificity in epizoic diatoms and the utility of diatoms as indicators of host ecology and health. We encourage the interdisciplinary community working with marine megafauna to consider including diatom sampling and diatom analysis into their routine practices.

Surface microbiota of Mediterranean loggerhead sea turtles unraveled by 16S and 18S amplicon sequencing

The loggerhead sea turtle is considered a keystone species with a major ecological role in Mediterranean marine environment. As is the case with other wild reptiles, their outer microbiome is rarely studied. Although there are several studies on sea turtle’s macro-epibionts and endo-microbiota, there has been little research on epibiotic microbiota associated with turtle skin and carapace. Therefore we aimed to provide the identification of combined epibiotic eukaryotic, bacterial and archaeal microbiota on Mediterranean loggerhead sea turtles. In this study, we sampled skins and carapaces of 26 loggerheads from the Mediterranean Sea during 2018 and 2019. To investigate the overall microbial diversity and composition, amplicon sequencing of 16S and 18S rRNA genes was performed. We found that the Mediterranean loggerhead sea turtle epibiotic microbiota is a reservoir of a vast variety of microbial species. Microbial communities mostly varied by different locations and seas, while within bacterial communities’ significant difference was observed between sampled body sites (carapace vs. skin). In terms of relative abundance, Proteobacteria and Bacteroidota were the most represented phyla within prokaryotes, while Alveolata and Stramenopiles thrived among eukaryotes. This study, besides providing a first survey of microbial eukaryotes on loggerheads via metabarcoding, identifies fine differences within both bacterial and eukaryotic microbial communities that seem to reflect the host anatomy and habitat. Multi-domain epi-microbiome surveys provide additional layers of information that are complementary with previous morphological studies and enable better understanding of the biology and ecology of these vulnerable marine reptiles.

Novel Approach to Freshwater Diatom Profiling and Identification Using Raman Spectroscopy and Chemometric Analysis

(1) An approach with great potential for fast and cost-effective profiling and identification of diatoms in lake ecosystems is presented herein. This approach takes advantage of Raman spectroscopy. (2) The study was based on the analysis of 790 Raman spectra from 29 species, belonging to 15 genera, 12 families, 9 orders and 4 subclasses, which were analysed using chemometric methods. The Raman data were first analysed by a partial least squares regression discriminant analysis (PLS-DA) to characterise the diatom species. Furthermore, a method was developed to streamline the integrated interpretation of PLS-DA when a high number of significant components is extracted. Subsequently, an artificial neural network (ANN) was used for taxa identification from Raman data. (3) The PLS interpretation produced a Raman profile for each species reflecting its biochemical composition. The ANN models were useful to identify various taxa with high accuracy. (4) Compared to studies in the literature, involving huge datasets one to four orders of magnitude larger than ours, high sensitivity was found for the identification of Achnanthidium exiguum (67%), Fragilaria pararumpens (67%), Amphora pediculus (71%), Achnanthidium minutissimum (80%) and Melosira varians (82%).

Meta-analysis shows both congruence and complementarity of DNA and eDNA metabarcoding to traditional methods for biological community assessment

DNA metabarcoding is increasingly used for the assessment of aquatic communities, and numerous studies have investigated the consistency of this technique with traditional morpho‐taxonomic approaches. These individual studies have used DNA metabarcoding to assess diversity and community structure of aquatic organisms both in marine and freshwater systems globally over the last decade. However, a systematic analysis of the comparability and effectiveness of DNA‐based community assessment across all of these studies has hitherto been lacking. Here, we performed the first meta‐analysis of available studies comparing traditional methods and DNA metabarcoding to measure and assess biological diversity of key aquatic groups, including plankton, microphytobentos, macroinvertebrates, and fish. Across 215 data sets, we found that DNA metabarcoding provides richness estimates that are globally consistent to those obtained using traditional methods, both at local and regional scale. DNA metabarcoding also generates species inventories that are highly congruent with traditional methods for fish. Contrastingly, species inventories of plankton, microphytobenthos and macroinvertebrates obtained by DNA metabarcoding showed pronounced differences to traditional methods, missing some taxa but at the same time detecting otherwise overseen diversity. The method is generally sufficiently advanced to study the composition of fish communities and replace more invasive traditional methods. For smaller organisms, like macroinvertebrates, plankton and microphytobenthos, DNA metabarcoding may continue to give complementary rather than identical estimates compared to traditional approaches. Systematic and comparable data collection will increase the understanding of different aspects of this complementarity, and increase the effectiveness of the method and adequate interpretation of the results.

Assessment of marine benthic diatom communities: insights from a combined morphological-metabarcoding approach in Mediterranean shallow coastal waters

We investigated the advantages and disadvantages of light microscope (LM)-based identifications and DNA metabarcoding, based on a 312-bp rbcL marker, for examining benthic diatom communities from Mediterranean shallow coastal environments. For this, we used biofilm samples collected from different substrata in the Ebro delta bays. We show that 1) Ebro delta bays harbour high-diversity diatom communities [LM identified 249 taxa] and 2) DNA metabarcoding effectively reflects this diversity at genus- but not species level, because of the incompleteness of the DNA reference library. Nevertheless, DNA metabarcoding offers new opportunities for detecting small, delicate and rare diatom species missed by LM and diatoms that lack silica frustules. The primers used, though designed for diatoms, successfully amplified rarely reported members of other stramenopile groups. Combining LM and DNA approaches offers stronger support for ecological studies of benthic microalgal communities in shallow coastal environments than using either approach on its own.

DNA metabarcoding reveals differences in distribution patterns and ecological preferences among genetic variants within some key freshwater diatom species

Our study evaluates differences in the distribution and ecology of genetic variants within several ecologically important diatom species that are also key for Water Framework Directive monitoring of European rivers: Fistulifera saprophila (FSAP), Achnanthidium minutissimum (ADMI), Nitzschia inconspicua (NINC) and Nitzschia soratensis (NSTS). We used DADA2 to infer amplicon sequence variants (ASVs) of a short rbcL barcode in 531 environmental samples from biomonitoring campaigns in Catalonia and France. ASVs within each species showed different distribution patterns. Threshold Indicator Taxa ANalysis revealed three ecological groupings of ASVs in both ADMI and FSAP. Two of these in each species were separated by opposite responses to calcium and conductivity. Boosted regression trees additionally showed that both variables greatly influenced the occurrence of these groupings. A third grouping in FSAP was characterized by a negative response to total organic carbon and hence was better represented in waters with higher ecological status than the other FSAP ASVs, contrasting with what is generally assumed for the species. In the two Nitzschia species, our analyses confirmed earlier studies: NINC preferred higher levels of calcium and conductivity. Our findings suggest that the broad ecological tolerance of some diatom species results from overlapping preferences among genetic variants, which individually show much more restricted preferences and distributions. This work shows the importance of studying the ecological preferences of genetic variants within species complexes, now possible with DNA metabarcoding. The results will help reveal and understand biogeographical distributions and facilitate the development of more accurate biological indexes for biomonitoring programmes.

Epibiont Assemblages on Nesting Hawksbill Turtles Show Site-Specificity in the Persian Gulf

Sea turtle epibionts can provide insights into the hosts' habitat use. However, at present, there is a lack of information on sea turtle epibiont communities in many locations worldwide. Here, we describe the epibiont communities of 46 hawksbill turtles (Eretmochelys imbricata) in the Persian Gulf. Specifically, we sampled 28 turtles from the Dayyer-Nakhiloo National Park (DNNP) in the northern Gulf and 18 turtles from Shibderaz beach in the Strait of Hormuz. A total of 54 macro, meio, and micro-epibiont taxa were identified, including 46 taxa from Shibderaz and 29 taxa from DNNP. The barnacles Chelonibia testudinaria and Platylepas hexastylos, as well as harpacticoid copepods and Rotaliid foraminifers, had the highest frequency of occurrence found on almost all turtle individuals. Harpacticoids were the most abundant epizoic taxa (19.55 ± 3.9 ind. per 9 cm2) followed by forams (Quinqueloculina spp.: 6.25 ± 1.5 ind. per 9 cm2 and Rotaliids: 6.02 ± 1.3 ind. per 9 cm2). Our results showed significant differences between the study sites in the composition of micro and macro-epibiont communities found on hawksbill turtles. We speculate that the differences in epibiont communities were largely influenced by local environmental conditions.

DNA metabarcoding and morphological methods show complementary patterns in the metacommunity organization of lentic epiphytic diatoms

Diatoms are important organisms in freshwater ecosystems due to their position as primary producers and therefore, analyzing their assemblages provides relevant information on ecosystem functioning. Diatoms have historically been identified based on morphological traits, which is time-consuming and requires well-trained specialists. Nevertheless, DNA barcoding offers an alternative approach to overcome some limitations of the morphological method. Here, we assess if both approaches are comparable methods to study patterns and mechanisms (including environmental filtering and dispersal limitation) of epiphytic diatom metacommunities using a comprehensive dataset from 22 Mediterranean ponds at different taxonomic resolutions. We used a fragment of rbcL barcode gene combined with High-Throughput Sequencing to infer diatom community composition. The overall degree of correspondence between both approaches was assessed by Procrustean rotation analysis and Procrustean randomization tests, whereas the role of local environmental variables and geographical distances was studied using a comprehensive combination of BIOENV, Mantel tests and distance-based redundancy analysis. Our results showed a relatively poor correspondence in the compositional variation of diatom metacommunity between both approaches. We speculate that the incompleteness of the reference database and the bioinformatics processing are the biases most likely affecting the molecular approach, whereas the limited counting effort and the presence of cryptic species are presumably the major biases related with the morphological method. On the other hand, variation in diatom community composition detected with both approaches was strongly related to the environmental template, which may be related with the narrow community-environment relationships in diatoms. Nevertheless, we found no significant relationship between compositional variation and geographical distances. Overall, our work shows the complementary nature of both approaches and highlights the importance of DNA metabarcoding to address empirical research questions of community ecology in freshwaters, especially once the reference databases include most genotypes of occurring taxa and bioinformatics biases are overcome.

Dumpster diving for diatom plastid 16S rRNA genes

High throughput sequencing is improving the efficiency of monitoring diatoms, which inhabit and support aquatic ecosystems across the globe. In this study, we explored the potential of a standard V4 515F-806RB primer pair in recovering diatom plastid 16S rRNA sequences. We used PhytoREF to classify the 16S reads from our freshwater biofilm field sampling from three stream segments across two streams in south-eastern Australia and retrieved diatom community data from other, publicly deposited, Australian 16S amplicon datasets. When these diatom operational taxonomic units (OTUs) were traced using the default RDPII and NCBI databases, 68% were characterized as uncultured cyanobacteria. We analysed the 16S rRNA sequences from 72 stream biofilm samples, separated the chloroplast OTUs, and classified them using the PhytoREF database. After filtering the reads attributed to Bacillariophyta (relative abundance >1%), 71 diatom OTUs comprising more than 90% of the diatom reads in each stream biofilm sample were identified. Beta-diversity analyses demonstrated significantly different diatom assemblages and discrimination among river segments. To further test the approach, the diatom OTUs from our biofilm sampling were used as reference sequences to identify diatom reads from other Australian 16S rRNA datasets in the NCBI-SRA database. Across the three selected public datasets, 67 of our 71 diatom OTUs were detected in other Australian ecosystems. Our results show that diatom plastid 16S rRNA genes are readily amplified with existing 515F-806RB primer sets. Therefore, the volume of existing 16S rRNA amplicon datasets initially generated for microbial community profiling can also be used to detect, characterize, and map diatom distribution to inform phylogeny and ecological health assessments, and can be extended into a range of ecological and industrial applications. To our knowledge, this study represents the first attempt to classify freshwater samples using this approach and the first application of PhytoREF in Australia.

Algal and cyanobacterial diversity in saline rivers of the Elton Lake Basin (Russia) studied via light microscopy and next-generation sequencing

Naturally saline rivers are known in various regions of the world. Saline rivers with a salinity gradient from the source to the mouth are particularly interesting, because the range of salinity is the structure-forming factor of the hydrobiont assemblage. Such rivers are represented by saline rivers of the Elton Lake Basin in Volgograd region of Russia (the Bolshaya Samoroda River and the Malaya Samoroda River). Herein, we analyzed taxonomic structure and species diversity of microalgae and Cyanobacteria of the saline rivers flowing into the Elton Lake by light microscopy and next-generation sequencing. The differences and possible causes of inconsistencies in the results obtained by these methods are discussed. In total, 91 taxa of microorganisms were identified by integrated approach in the assemblages of microalgae and Cyanobacteria in the middle course of the Bolshaya Samoroda River, and 60 taxa – in the river mouth. The species diversity of those assemblages in the hypersaline Malaya Samoroda River was lower: 27 taxa from the middle course and 23 taxa from the mouth. Next-generation sequencing allowed us to refine and expand the list of microalgae taxa in the studied saline rivers due to detection of species which were hard to identify, low-abundance taxa, as well as extremely small-cell forms. Some discrepancies between the data obtained by light microscopy and next-generation sequencing indicate the advantage of simultaneous use of both methods for study of the algae communities. Such a comprehensive approach provides the most accurate and correct list of taxa added with the morphological descriptions and 18S rRNA and 16S rRNA partial sequences. Generally, 18 taxa have been recorded for the first time in the Bolshaya Samoroda River, belonging to the phyla Chlorophyta (Borodinellopsis sp., Chlorochytrium lemnae Cohn, Caespitella sp., Halochlorococcum sp., Tetraselmis cordiformis (H. J. Carter) F. Stein), Ochrophyta (Pseudocharaciopsis ovalis (Chodat) D. J. Hibberd, Characiopsis sp., Poterioochromonas stipitata Scherffel, Chrysolepidomonas sp.), Euglenozoa (Euglena bucharica I. Kisselev, Lepocinclis tripteris (Dujardin) B. Marin & Melkonian, Phacus orbicularis K. Hübner, P. parvulus G. A. Klebs), Cryptophyta (Hemiselmis cryptochromatica C. E. Lane & J. M. Archibald, Rhodomonas sp., Hanusia phi J. A. Deane), Haptophyta (Pavlova sp.), Cyanobacteria (Johanseninema constrictum (Szafer) Hasler, Dvorák & Poulícková). Seven taxa have been detected for the first time in the algal and cyanobacterial assemblages of the Malaya Samoroda River from the phyla Chlorophyta (Tetraselmis cordiformis, T. arnoldii (Proschkina-Lavrenko) R. E. Norris, Hori & Chihara, T. tetrathele (West) Butcher, Pyrobotrys elongatus Korshikov), Cryptophyta (Hanusia phi), and Cyanobacteria (Synechococcus elongatus (Nägeli) Nägeli, Oscillatoria simplicissima Gomont).

Benthic Diatoms in River Biomonitoring—Present and Future Perspectives within the Water Framework Directive

The European Water Framework Directive 2000/60/EC (WFD) has been implemented over the past 20 years, using physicochemical, biological and hydromorphological elements to assess the ecological status of surface waters. Benthic diatoms (i.e., phytobenthos) are one of the most common biological quality elements (BQEs) used in surface water monitoring and are particularly successful in detecting eutrophication, organic pollution and acidification. Herein, we reviewed their implementation in river biomonitoring for the purposes of the WFD, highlighting their advantages and disadvantages over other BQEs, and we discuss recent advances that could be applied in future biomonitoring. Until now, phytobenthos have been intercalibrated by the vast majority (26 out of 28) of EU Member States (MS) in 54% of the total water bodies assessed and was the most commonly used BQE after benthic invertebrates (85% of water bodies), followed by fish (53%), macrophytes (27%) and phytoplankton (4%). To meet the WFD demands, numerous taxonomy-based quality indices have been developed among MS, presenting, however, uncertainties possibly related to species biogeography. Recent development of different types of quality indices (trait-based, DNA sequencing and predictive modeling) could provide more accurate results in biomonitoring, but should be validated and intercalibrated among MS before their wide application in water quality assessments.

Consistent patterns in 16S and 18S microbial diversity from the shells of the common and widespread red-eared slider turtle (Trachemys scripta)

Microbial communities associated with freshwater aquatic habitats and resident species are both critical to and indicative of ecosystem status and organismal health. External surfaces of turtle shells readily accumulate microbial growth and could carry representation of habitat-wide microbial diversity, since they are in regular contact with multiple elements of freshwater environments. Yet, microbial diversity residing on freshwater turtle shells is poorly understood. We applied 16S and 18S metabarcoding to characterize microbiota associated with external shell surfaces of 20 red-eared slider (Trachemys scripta) turtles collected from varied habitats in central and western Oklahoma, and ranging to southeast Iowa. Shell-associated microbial communities were highly diverse, with samples dominated by Bacteroidia and alpha-/gamma-proteobacteria, and ciliophoran alveolates. Alpha diversity was lower on turtle shells compared to shallow-water-associated environmental samples, likely resulting from basking-drying behavior and seasonal scute shedding, while alpha diversity was higher on carapace than plastron surfaces. Beta diversity of turtle shells was similarly differentiated from environmental samples, although sampling site was consistently a significant factor. Deinococcus-Thermus bacteria and ciliophoran alveolates were recovered with significantly higher abundance on turtle shells versus environmental samples, while bacterial taxa known to include human-pathogenic species were variably more abundant between shell and environmental samples. Microbial communities from a single, shared-site collection of the ecologically similar river cooter (P. concinna) largely overlapped with those of T. scripta. These data add to a foundation for further characterization of turtle shell microbial communities across species and habitats, with implications for freshwater habitat assessment, microbial ecology and wildlife conservation efforts.

Geographical variation in the diatom communities associated with loggerhead sea turtles (Caretta caretta)

Epizoic diatoms form an important part of micro-epibiota of marine vertebrates such as whales and sea turtles. The present study explores and compares the diversity and biogeography of diatom communities growing on the skin and shell of loggerhead sea turtles (Caretta caretta) from four different localities: Adriatic Sea (Croatia), Ionian Sea (Greece), South Africa and Florida Bay (USA) using both light and scanning electron microscopy. We observed almost 400 diatom taxa belonging to more than 100 genera. Diatom communities from Greece and Croatia showed the highest similarity and were statistically different from those recorded from South Africa and Florida. Part of this variation could be attributed to differences in sampling techniques; however, we believe that geography had an important role. In general, contrary to several previous observations from sea turtles, the presumably exclusively epizoic diatoms contributed less than common benthic taxa to the total diatom flora, which might have been related to the loggerhead feeding behavior. Moreover, skin samples differed from carapace samples in having a distinct diatom composition with a higher proportion of the putative true epizoonts. Our results indicate that epizoic diatom communities differ according to loggerhead geographical location and substrate (skin vs. carapace). The relative abundances of common benthic diatoms and putative exclusive epizoic taxa may inform about sea turtle habitat use or behavior though detailed comparisons among different host species have yet to be performed.

Evaluation and sensitivity analysis of diatom DNA metabarcoding for WFD bioassessment of Mediterranean rivers

Our study of 164 diatom samples from Catalonia (NE Spain) is the first to evaluate the applicability of DNA metabarcoding, based on high throughput sequencing (HTS) using a 312-bp rbcL marker, for biomonitoring Mediterranean rivers. For this, we compared the values of a biotic index (IPS) and the ecological status classes derived from them, between light microscope-based (LM) and HTS methods. Very good correspondence between methods gives encouraging results concerning the applicability of DNA metabarcoding for Catalan rivers for the EU Water Framework Directive (WFD). However, in 10 sites, the ecological status class was downgraded from "Good"/"High" obtained by LM to "Moderate"/"Poor"/"Bad" by HTS; these "critical" sites are especially important, because the WFD requires remedial action by water managers for any river with Moderate or lower status. We investigated the contribution of each species to the IPS using a "leave-one-out" sensitivity analysis, paying special attention to critical sites. Discrepancies in IPS between LM and HTS were mainly due to the misidentification and overlooking in LM of a few species, which were better recovered by HTS. This bias was particularly important in the case of Fistulifera saprophila, whose clear underrepresentation in LM was important for explaining 8 out of the 10 critical sites and probably reflected destruction of weakly-silicified frustules during sample preparation. Differences between species in the rbcL copy number per cell affected the relative abundance obtained by HTS for Achnanthidium minutissimum, Nitzschia inconspicua and Ulnaria ulna, which were also identified by the sensitivity analysis as important for the WFD. Only minor IPS discrepancies were attributed to the incompleteness of the reference library, as most of the abundant and influential species (to the IPS) were well represented there. Finally, we propose that leave-one-out analysis is a good method for identifying priority species for isolation and barcoding.

Biodiversity of carapace epibiont diatoms in loggerhead sea turtles (Caretta caretta Linnaeus 1758) in the Aegean Sea Turkish coast

Background

The Aegean Sea coast of Turkey hosts one of the most important nesting grounds for loggerhead sea turtles (Caretta caretta) in the Mediterranean Sea. Previous studies have revealed that the sea turtle carapace provides favourable conditions for various epibiontic organisms. Epibionts occurring on the carapace have been examined from different locations in the oceans.

Methods

This is the first time such a high number (39) of samples collected from nesting turtles during such a long time period (extending from 2011 to 2018) has been used for the study of the diatom component of the microbiome on the turtle carapaces. A total of 33 samples were investigated in terms of light microscopy (LM) and scanning electron microscopy (SEM). Six unprocessed biofilm fragments were subject to SEM observations.

Results

A total of 457 epizoic diatom taxa belonging to 86 genera were identified. Epizoic forms, e.g., Achnanthes spp., Chelonicola spp. or Tripterion spp. (also identified by SEM observations of the undisturbed pieces of the microbiome) dominated in terms of relative abundance, but the highest numbers of taxa were ubiquitously represented by Navicula (79), Nitzschia (45), Amphora (40), Cocconeis (32), Diploneis (25) and Mastogloia (23). Navicula perminuta and Delphineis australis were the most frequent taxa, present in 65% of the samples, both with an average relative abundance of 10%. The results of our study revealed that diatoms are an essential component of the loggerhead sea turtles’ microbiome, in terms of high biodiversity and abundance. Although strict epibionts provide a signature of the turtle microbiome, the carapace as a solid substrate attracts numerous benthic diatom species which are considered opportunistic forms and can be found in the surrounding benthic habitats of the vast ocean littoral space.

Diat.barcode, an open-access curated barcode library for diatoms

Diatoms (Bacillariophyta) are ubiquitous microalgae which produce a siliceous exoskeleton and which make a major contribution to the productivity of oceans and freshwaters. They display a huge diversity, which makes them excellent ecological indicators of aquatic ecosystems. Usually, diatoms are identified using characteristics of their exoskeleton morphology. DNA-barcoding is an alternative to this and the use of High-Throughput-Sequencing enables the rapid analysis of many environmental samples at a lower cost than analyses under microscope. However, to identify environmental sequences correctly, an expertly curated reference library is needed. Several curated libraries for protists exists; none, however are dedicated to diatoms. Diat.barcode is an open-access library dedicated to diatoms which has been maintained since 2012. Data come from two sources (1) the NCBI nucleotide database and (2) unpublished sequencing data of culture collections. Since 2017, several experts have collaborated to curate this library for rbcL, a chloroplast marker suitable for species-level identification of diatoms. For the latest version of the database (version 7), 605 of the 3482 taxonomical names originally assigned by the authors of the rbcL sequences were modified after curation. The database is accessible at https://www6.inra.fr/carrtel-collection_eng/Barcoding-database.

DNA barcode reference libraries for the monitoring of aquatic biota in Europe: Gap-analysis and recommendations for future work

Effective identification of species using short DNA fragments (DNA barcoding and DNA metabarcoding) requires reliable sequence reference libraries of known taxa. Both taxonomically comprehensive coverage and content quality are important for sufficient accuracy. For aquatic ecosystems in Europe, reliable barcode reference libraries are particularly important if molecular identification tools are to be implemented in biomonitoring and reports in the context of the EU Water Framework Directive (WFD) and the Marine Strategy Framework Directive (MSFD). We analysed gaps in the two most important reference databases, Barcode of Life Data Systems (BOLD) and NCBI GenBank, with a focus on the taxa most frequently used in WFD and MSFD. Our analyses show that coverage varies strongly among taxonomic groups, and among geographic regions. In general, groups that were actively targeted in barcode projects (e.g. fish, true bugs, caddisflies and vascular plants) are well represented in the barcode libraries, while others have fewer records (e.g. marine molluscs, ascidians, and freshwater diatoms). We also found that species monitored in several countries often are represented by barcodes in reference libraries, while species monitored in a single country frequently lack sequence records. A large proportion of species (up to 50%) in several taxonomic groups are only represented by private data in BOLD. Our results have implications for the future strategy to fill existing gaps in barcode libraries, especially if DNA metabarcoding is to be used in the monitoring of European aquatic biota under the WFD and MSFD. For example, missing species relevant to monitoring in multiple countries should be prioritized for future collaborative programs. We also discuss why a strategy for quality control and quality assurance of barcode reference libraries is needed and recommend future steps to ensure full utilisation of metabarcoding in aquatic biomonitoring.

The impact of OTU sequence similarity threshold on diatom-based bioassessment: A case study of the rivers of Mayotte (France, Indian Ocean)

Extensive studies on the taxonomic resolution required for bioassessment purposes have determined that resolution above species level (genus, family) is sufficient for their use as indicators of relevant environmental pressures. The high-throughput sequencing (HTS) and meta-barcoding methods now used for bioassessment traditionally employ an arbitrary sequence similarity threshold (SST) around 95% or 97% to cluster sequences into operational taxonomic units, which is considered descriptive of species-level resolution. In this study, we analyzed the effect of the SST on the resulting diatom-based ecological quality index, which is based on OTU abundance distribution along a defined environmental gradient, ideally avoiding taxonomic assignments that could result in high rates of unclassified OTUs and biased final values. A total of 90 biofilm samples were collected in 2014 and 2015 from 51 stream sites on Mayotte Island in parallel with measures of relevant physical and chemical parameters. HTS sequencing was performed on the biofilms using the rbcL region as the genetic marker and diatom-specific primers. Hierarchical clustering was used to group sequences into OTUs using 20 experimental SST levels (80%-99%). An OTU-based quality index (IdxOTU) was developed based on a weighted average equation using the abundance profiles of the OTUs. The developed IdxOTU revealed significant correlations between the IdxOTU values and the reference pressure gradient, which reached maximal performance using an SST of 90% (well above species level delimitation). We observed an interesting and important trade-off with the power to discriminate between sampling sites and index stability that will greatly inform future applications of the index. Taken together, the results from this study detail a thoroughly optimized and validated approach to generating robust, reproducible, and complete indexes that will greatly facilitate effective and efficient environmental monitoring.