Assessment of species diversity and distribution of an ancient diatom lineage using a DNA metabarcoding approach

A comparison of two gene regions for assessing community composition of eukaryotic marine microalgae from coastal ecosystems

Two gene regions commonly used to characterise the diversity of eukaryotic communities using metabarcoding are the 18S ribosomal DNA V4 and V9 gene regions. We assessed the effectiveness of these two regions for characterising diverisity of coastal eukaryotic microalgae communities (EMCs) from tropical and temperate sites. We binned amplicon sequence variants (ASVs) into the high level taxonomic groups: dinoflagellates, pennate diatoms, radial centric diatoms, polar centric diatoms, chlorophytes, haptophytes and 'other microalgae'. When V4 and V9 generated ASV abundances were compared, the V9 region generated a higher number of raw reads, captured more diversity from all high level taxonomic groups and was more closely aligned with the community composition determined using light microscopy. The V4 region did resolve more ASVs to a deeper taxonomic resolution within the dinoflagellates, but did not effectively resolve other major taxonomic divisions. When characterising these communities via metabarcoding, the use of multiple gene regions is recommended, but the V9 gene region can be used in isolation to provide high-level community biodiversity to reflect relative abundances within groups. This approach reduces the cost of sequencing multiple gene regions whilst still providing important baseline ecosystem function information.

Metabarcoding reveals high genetic diversity of harmful algae in the coastal waters of Texas, Gulf of Mexico

Environmental-DNA (eDNA) for metabarcoding is a rapid and effective means to investigate microplankton community composition and species diversity. The objective of this study was to examine the genetic diversity of the phytoplankton community in the Gulf of Mexico, with particular emphasis on harmful algal bloom species. Samples were collected at stations along the coast of Texas in September-October 2017 that were inundated by low salinity waters in the aftermath of Hurricane Harvey. Metabarcodes were generated from the eDNA targeting both the V4 and V8-V9 regions of the 18S rDNA gene. Evaluation of the metabarcodes revealed an unexpectedly high number of harmful algal species during this short period, including five that had not been documented in this region previously. A total of 36 harmful algal species could be differentiated based on V4 and V8-V9 metabarcode markers. Using a phylogenetic approach, the taxonomic resolution of each marker differed and not all species could be differentiated using solely one marker. The V4 region resolved species within some genera (e.g., Heterocapsa), while the V8-V9 marker was necessary to resolve species within other genera (e.g., Chattonella). In other cases, species differentiation within a genus required a combination of both markers (e.g., Prorocentrum, Karenia), or another marker will be needed to resolve all species (e.g., Alexandrium, Dinophysis). We conclude that no single marker can delineate all species, so it is recommended HAB monitoring programs use more than one marker. Overall, the observed diversity of HAB species along the Texas coast using metabarcoding exceeded reports from other parts of the world.

Resolving cryptic species complexes in marine protists: phylogenetic haplotype networks meet global DNA metabarcoding datasets

Marine protists have traditionally been assumed to be lowly diverse and cosmopolitan. Yet, several recent studies have shown that many protist species actually consist of cryptic complexes of species whose members are often restricted to particular biogeographic regions. Nonetheless, detection of cryptic species is usually hampered by sampling coverage and application of methods (e.g. phylogenetic trees) that are not well suited to identify relatively recent divergence and ongoing gene flow. In this paper, we show how these issues can be overcome by inferring phylogenetic haplotype networks from global metabarcoding datasets. We use the Chaetoceros curvisetus (Bacillariophyta) species complex as study case. Using two complementary metabarcoding datasets (Ocean Sampling Day and Tara Oceans), we equally resolve the cryptic complex in terms of number of inferred species. We detect new hypothetical species in both datasets. Gene flow between most of species is absent, but no barcoding gap exists. Some species have restricted distribution patterns whereas others are widely distributed. Closely related taxa occupy contrasting biogeographic regions, suggesting that geographic and ecological differentiation drive speciation. In conclusion, we show the potential of the analysis of metabarcoding data with evolutionary approaches for systematic and phylogeographic studies of marine protists.

Biases in bulk: DNA metabarcoding of marine communities and the methodology involved

With the growing anthropogenic pressure on marine ecosystems, the need for efficient monitoring of biodiversity grows stronger. DNA metabarcoding of bulk samples is increasingly being implemented in ecosystem assessments and is more cost-efficient and less time-consuming than monitoring based on morphology. However, before raw sequences are obtained from bulk samples, a profound number of methodological choices must be made. Here, we critically review the recent methods used for metabarcoding of marine bulk samples (including benthic, plankton and diet samples) and indicate how potential biases can be introduced throughout sampling, preprocessing, DNA extraction, marker and primer selection, PCR amplification and sequencing. From a total of 64 studies evaluated, our recommendations for best practices include to (a) consider DESS as a fixative instead of ethanol, (b) use the DNeasy PowerSoil kit for any samples containing traces of sediment, (c) not limit the marker selection to COI only, but preferably include multiple markers for higher taxonomic resolution, (d) avoid touchdown PCR profiles, (e) use a fixed annealing temperature for each primer pair when comparing across studies or institutes, (f) use a minimum of three PCR replicates, and (g) include both negative and positive controls. Although the implementation of DNA metabarcoding still faces several technical complexities, we foresee wide-ranging advances in the near future, including improved bioinformatics for taxonomic assignment, sequencing of longer fragments and the use of whole-genome information. Despite the bulk of biases involved in metabarcoding of bulk samples, if appropriate controls are included along the data generation process, it is clear that DNA metabarcoding provides a valuable tool in ecosystem assessments.

Genome-enabled phylogenetic and functional reconstruction of an araphid pennate diatom Plagiostriata sp. CCMP470, previously assigned as a radial centric diatom, and its bacterial commensal

Diatoms are an ecologically fundamental and highly diverse group of algae, dominating marine primary production in both open-water and coastal communities. The diatoms include both centric species, which may have radial or polar symmetry, and the pennates, which include raphid and araphid species and arose within the centric lineage. Here, we use combined microscopic and molecular information to reclassify a diatom strain CCMP470, previously annotated as a radial centric species related to Leptocylindrus danicus, as an araphid pennate species in the staurosiroid lineage, within the genus Plagiostriata. CCMP470 shares key ultrastructural features with Plagiostriata taxa, such as the presence of a sternum with parallel striae, and the presence of a highly reduced labiate process on its valve; and this evolutionary position is robustly supported by multigene phylogenetic analysis. We additionally present a draft genome of CCMP470, which is the first genome available for a staurosiroid lineage. 270 Pfams (19%) found in the CCMP470 genome are not known in other diatom genomes, which otherwise does not hold big novelties compared to genomes of non-staurosiroid diatoms. Notably, our DNA library contains the genome of a bacterium within the Rhodobacterales, an alpha-proteobacterial lineage known frequently to associate with algae. We demonstrate the presence of commensal alpha-proteobacterial sequences in other published algal genome and transcriptome datasets, which may indicate widespread and persistent co-occurrence.

Lake-depth related pattern of genetic and morphological diatom diversity in boreal Lake Bolshoe Toko, Eastern Siberia

Large, old and heterogenous lake systems are valuable sources of biodiversity. The analysis of current spatial variability within such lakes increases our understanding of the origin and establishment of biodiversity. The environmental sensitivity and the high taxonomic richness of diatoms make them ideal organisms to investigate intra-lake variability. We investigated modern intra-lake diatom diversity in the large and old sub-arctic Lake Bolshoe Toko in Siberia. Our study uses diatom-specific metabarcoding, applying a short rbcL marker combined with next-generation sequencing and morphological identification to analyse the diatom diversity in modern sediment samples of 17 intra-lake sites. We analysed abundance-based compositional taxonomic diversity and generic phylogenetic diversity to investigate the relationship of diatom diversity changes with water depth. The two approaches show differences in taxonomic identification and alpha diversity, revealing a generally higher diversity with the genetic approach. With respect to beta diversity and ordination analyses, both approaches result in similar patterns. Water depth or related lake environmental conditions are significant factors influencing intra-lake diatom patterns, showing many significant negative correlations between alpha and beta diversity and water depth. Further, one near-shore and two lagoon lake sites characterized by low (0-10m) and medium (10-30m) water depth are unusual with unique taxonomic compositions. At deeper (>30m) water sites we identified strongest phylogenetic clustering in Aulacoseira, but generally much less in Staurosira, which supports that water depth is a strong environmental filter on the Aulacoseira communities. Our study demonstrates the utility of combining analyses of genetic and morphological as well as phylogenetic diversity to decipher compositional and generic phylogenetic patterns, which are relevant in understanding intra-lake heterogeneity as a source of biodiversity in the sub-arctic glacial Lake Bolshoe Toko.

Intraspecific Diversity in the Cold Stress Response of Transposable Elements in the Diatom Leptocylindrus aporus

Transposable elements (TEs), activated as a response to unfavorable conditions, have been proposed to contribute to the generation of genetic and phenotypic diversity in diatoms. Here we explore the transcriptome of three warm water strains of the diatom Leptocylindrus aporus, and the possible involvement of TEs in their response to changing temperature conditions. At low temperature (13 °C) several stress response proteins were overexpressed, confirming low temperature to be unfavorable for L. aporus, while TE-related transcripts of the LTR retrotransposon superfamily were the most enriched transcripts. Their expression levels, as well as most of the stress-related proteins, were found to vary significantly among strains, and even within the same strains analysed at different times. The lack of overexpression after many months of culturing suggests a possible role of physiological plasticity in response to growth under controlled laboratory conditions. While further investigation on the possible central role of TEs in the diatom stress response is warranted, the strain-specific responses and possible role of in-culture evolution draw attention to the interplay between the high intraspecific variability and the physiological plasticity of diatoms, which can both contribute to the adaptation of a species to a wide range of conditions in the marine environment.

Time series and beyond: multifaceted plankton research at a marine Mediterranean LTER site

Plankton are a pivotal component of the diversity and functioning of coastal marine ecosystems. A long time-series of observations is the best tool to trace their patterns and variability over multiple scales, ultimately providing a sound foundation for assessing, modelling and predicting the effects of anthropogenic and natural environmental changes on pelagic communities. At the same time, a long time-series constitutes a formidable asset for different kinds of research on specific questions that emerge from the observations, whereby the results of these complementary studies provide precious interpretative tools that augment the informative value of the data collected. In this paper, we review more than 140 studies that have been developed around a Mediterranean plankton time series gathered in the Gulf of Naples at the station LTER-MC since 1984. These studies have addressed different topics concerning marine plankton, which have included: i) seasonal patterns and trends; ii) taxonomic diversity, with a focus on key or harmful algal species and the discovery of many new taxa; iii) molecular diversity of selected species, groups of species or the whole planktonic community; iv) life cycles of several phyto- and zooplankton species; and v) interactions among species through trophic relationships, parasites and viruses. Overall, the products of this research demonstrate the great value of time series besides the record of fluctuations and trends, and highlight their primary role in the development of the scientific knowledge of plankton much beyond the local scale.

Diatom diversity through HTS-metabarcoding in coastal European seas

Diatoms constitute a diverse lineage of unicellular organisms abundant and ecologically important in aquatic ecosystems. Compared to other protists, their biology and taxonomy are well-studied, offering the opportunity to combine traditional approaches and new technologies. We examined a dataset of diatom 18S rRNA- and rDNA- (V4 region) reads from different plankton size-fractions and sediments from six European coastal marine sites, with the aim of identifying peculiarities and commonalities with respect to the whole protistan community. Almost all metabarcodes (99.6%) were assigned to known genera (121) and species (236), the most abundant of which were those already known from classic studies and coincided with those seen in light microscopy. rDNA and rRNA showed comparable patterns for the dominant taxa, but rRNA revealed a much higher diversity particularly in the sediment communities. Peculiar to diatoms is a tight bentho-pelagic coupling, with many benthic or planktonic species colonizing both water column and sediments and the dominance of planktonic species in both habitats. Overall metabarcoding results reflected the marked specificity of diatoms compared to other protistan groups in terms of morphological and ecological characteristics, at the same time confirming their great potential in the description of protist communities.

Tackling critical parameters in metazoan meta-barcoding experiments: a preliminary study based on coxI DNA barcode

Nowadays DNA meta-barcoding is a powerful instrument capable of quickly discovering the biodiversity of an environmental sample by integrating the DNA barcoding approach with High Throughput Sequencing technologies. It mainly consists of the parallel reading of informative genomic fragment/s able to discriminate living entities. Although this approach has been widely studied, it still needs optimization in some necessary steps requested in its advanced accomplishment. A fundamental element concerns the standardization of bioinformatic analyses pipelines. The aim of the present study was to underline a number of critical parameters of laboratory material preparation and taxonomic assignment pipelines in DNA meta-barcoding experiments using the cytochrome oxidase subunit-I (coxI) barcode region, known as a suitable molecular marker for animal species identification. We compared nine taxonomic assignment pipelines, including a custom in-house method, based on Hidden Markov Models. Moreover, we evaluated the potential influence of universal primers amplification bias in qPCR, as well as the correlation between GC content with taxonomic assignment results. The pipelines were tested on a community of known terrestrial invertebrates collected by pitfall traps from a chestnut forest in Italy. Although the present analysis was not exhaustive and needs additional investigation, our results suggest some potential improvements in laboratory material preparation and the introduction of additional parameters in taxonomic assignment pipelines. These include the correct setup of OTU clustering threshold, the calibration of GC content affecting sequencing quality and taxonomic classification, as well as the evaluation of PCR primers amplification bias on the final biodiversity pattern. Thus, careful attention and further validation/optimization of the above-mentioned variables would be required in a DNA meta-barcoding experimental routine.

Diatom flagellar genes and their expression during sexual reproduction in Leptocylindrus danicus

Background

Flagella have been lost in the vegetative phase of the diatom life cycle, but they are still present in male gametes of centric species, thereby representing a hallmark of sexual reproduction. This process, besides maintaining and creating new genetic diversity, in diatoms is also fundamental to restore the maximum cell size following its reduction during vegetative division. Nevertheless, sexual reproduction has been demonstrated in a limited number of diatom species, while our understanding of its different phases and of their genetic control is scarce.

Results

In the transcriptome of Leptocylindrus danicus, a centric diatom widespread in the world’s seas, we identified 22 transcripts related to the flagella development and confirmed synchronous overexpression of 6 flagellum-related genes during the male gamete formation process. These transcripts were mostly absent in the closely related species L. aporus, which does not have sexual reproduction. Among the 22 transcripts, L. danicus showed proteins that belong to the Intra Flagellar Transport (IFT) subcomplex B as well as IFT-A proteins, the latter previously thought to be absent in diatoms. The presence of flagellum-related proteins was also traced in the transcriptomes of several other centric species. Finally, phylogenetic reconstruction of the IFT172 and IFT88 proteins showed that their sequences are conserved across protist species and have evolved similarly to other phylogenetic marker genes.

Conclusion

Our analysis describes for the first time the diatom flagellar gene set, which appears to be more complete and functional than previously reported based on the genome sequence of the model centric diatom, Thalassiosira pseudonana. This first recognition of the whole set of diatom flagellar genes and of their activation pattern paves the way to a wider recognition of the relevance of sexual reproduction in individual species and in the natural environment.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4210-8) contains supplementary material, which is available to authorized users.

Eukaryotic diversity in late Pleistocene marine sediments around a shallow methane hydrate deposit in the Japan Sea

Marine sediments contain eukaryotic DNA deposited from overlying water columns. However, a large proportion of deposited eukaryotic DNA is aerobically biodegraded in shallow marine sediments. Cold seep sediments are often anaerobic near the sediment-water interface, so eukaryotic DNA in such sediments is expected to be preserved. We investigated deeply buried marine sediments in the Japan Sea, where a methane hydrate deposit is associated with cold seeps. Quantitative PCR analysis revealed the reproducible recovery of eukaryotic DNA in marine sediments at depths up to 31.0 m in the vicinity of the methane hydrate deposit. In contrast, the reproducible recovery of eukaryotic DNA was limited to a shallow depth (8.3 m) in marine sediments not adjacent to the methane hydrate deposit in the same area. Pyrosequencing of an 18S rRNA gene variable region generated 1,276-3,307 reads per sample, which was sufficient to cover the biodiversity based on rarefaction curves. Phylogenetic analysis revealed that most of the eukaryotic DNA originated from radiolarian genera of the class Chaunacanthida, which have SrSO4 skeletons, the sea grass genus Zostera, and the seaweed genus Sargassum. Eukaryotic DNA originating from other planktonic fauna and land plants was also detected. Diatom sequences closely related to Thalassiosira spp., indicative of cold climates, were obtained from sediments deposited during the last glacial period (MIS-2). Plant sequences of the genera Alnus, Micromonas, and Ulmus were found in sediments deposited during the warm interstadial period (MIS-3). These results suggest the long-term persistence of eukaryotic DNA from terrestrial and aquatic sources in marine sediments associated with cold seeps, and that the genetic information from eukaryotic DNA from deeply buried marine sediments associated with cold seeps can be used to reconstruct environments and ecosystems from the past.

A high resolution melting method for the molecular identification of the potentially toxic diatom Pseudo-nitzschia spp. in the Mediterranean Sea

The aim of this study was to develop and validate a high resolution melting (HRM) method for the rapid, accurate identification of the various harmful diatom Pseudo-nitzschia species in marine environments. Pseudo-nitzschia has a worldwide distribution and some species are toxic, producing the potent domoic acid toxin, which poses a threat to both human and animal health. Hence, it is important to identify toxic Pseudo-nitzschia species. A pair of primers targeting the LSU rDNA of the genus Pseudo-nitzschia was designed for the development of the assay and its specificity was validated using 22 control DNAs of the P. calliantha, P. delicatissima/P. arenysensis complex and P. pungens. The post-PCR HRM assay was applied to numerous unidentified Pseudo-nitzschia strains isolated from the northwestern Adriatic Sea (Mediterranean Sea), and it was able to detect and discriminate three distinct Pseudo-nitzschia taxa from unidentified samples. Moreover, the species-specific identification of Pseudo-nitzschia isolates by the HRM assay was consistent with phylogenetic analyses. The HRM assay was specific, robust and rapid when applied to high numbers of cultured samples in order to taxonomically identify Pseudo-nitzschia isolates recovered from environmental samples.

Strengths and Biases of High-Throughput Sequencing Data in the Characterization of Freshwater Ciliate Microbiomes

Molecular surveys of eukaryotic microbial communities employing high-throughput sequencing (HTS) techniques are rapidly supplanting traditional morphological approaches due to their larger data output and reduced bench work time. Here, we directly compare morphological and Illumina data obtained from the same samples, in an effort to characterize ciliate faunas from sediments in freshwater environments. We show how in silico processing affects the final outcome of our HTS analysis, providing evidence that quality filtering protocols strongly impact the number of predicted taxa, but not downstream conclusions such as biogeography patterns. We determine the abundance distribution of ciliates, showing that a small fraction of abundant taxa dominates read counts. At the same time, we advance reasons to believe that biases affecting HTS abundances may be significant enough to blur part of the underlying biological picture. We confirmed that the HTS approach detects many more taxa than morphological inspections, and highlight how the difference varies among taxonomic groups. Finally, we hypothesize that the two datasets actually correspond to different conceptions of "diversity," and consequently that neither is entirely superior to the other when investigating environmental protists.

Benthic protists: the under-charted majority

Marine protist diversity inventories have largely focused on planktonic environments, while benthic protists have received relatively little attention. We therefore hypothesize that current diversity surveys have only skimmed the surface of protist diversity in marine sediments, which may harbor greater diversity than planktonic environments. We tested this by analyzing sequences of the hypervariable V4 18S rRNA from benthic and planktonic protist communities sampled in European coastal regions. Despite a similar number of OTUs in both realms, richness estimations indicated that we recovered at least 70% of the diversity in planktonic protist communities, but only 33% in benthic communities. There was also little overlap of OTUs between planktonic and benthic communities, as well as between separate benthic communities. We argue that these patterns reflect the heterogeneity and diversity of benthic habitats. A comparison of all OTUs against the Protist Ribosomal Reference database showed that a higher proportion of benthic than planktonic protist diversity is missing from public databases; similar results were obtained by comparing all OTUs against environmental references from NCBI's Short Read Archive. We suggest that the benthic realm may therefore be the world's largest reservoir of marine protist diversity, with most taxa at present undescribed.

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.

The diatom molecular toolkit to handle nitrogen uptake

Nutrient concentrations in the oceans display significant temporal and spatial variability, which strongly affects growth, distribution and survival of phytoplankton. Nitrogen (N) in particular is often considered a limiting resource for prominent marine microalgae, such as diatoms. Diatoms possess a suite of N-related transporters and enzymes and utilize a variety of inorganic (e.g., nitrate, NO3(-); ammonium, NH4(+)) and organic (e.g., urea; amino acids) N sources for growth. However, the molecular mechanisms allowing diatoms to cope efficiently with N oscillations by controlling uptake capacities and signaling pathways involved in the perception of external and internal clues remain largely unknown. Data reported in the literature suggest that the regulation and the characteristic of the genes, and their products, involved in N metabolism are often diatom-specific, which correlates with the peculiar physiology of these organisms for what N utilization concerns. Our study reveals that diatoms host a larger suite of N transporters than one would expected for a unicellular organism, which may warrant flexible responses to variable conditions, possibly also correlated to the phases of life cycle of the cells. All this makes N transporters a crucial key to reveal the balance between proximate and ultimate factors in diatom life.