A full-length 18S ribosomal DNA metabarcoding approach for determining protist community diversity using Nanopore sequencing.

Protist diversity studies are frequently conducted using DNA metabarcoding methods. Currently, most studies have utilized short read sequences to assess protist diversity. One limitation of using short read sequences is the low resolution of the markers. For better taxonomic resolution longer sequences of the 18S rDNA are required because the full-length has both conserved and hypervariable regions. In this study, a new primer pair combination was used to amplify the full-length 18S rDNA and its efficacy was validated with a test community and then validated with field samples. Full-length sequences obtained with the Nanopore MinION for protist diversity from field samples were compared with Illumina MiSeq V4 and V8-V9 short reads. Sequences generated from the high-throughput sequencers are Amplicon Sequence Variants (ASVs). Metabarcoding results show high congruency among the long reads and short reads in taxonomic annotation at the major taxonomic group level; however, not all taxa could be successfully detected from sequences. Based on the criteria of ≥95% similarity and ≥1000 bp query length, 298 genera were identified by all markers in the field samples, 250 (84%) were detected by 18S, while only 226 (76%) by V4 and 213 (71%) by V8-V9. Of the total 85 dinoflagellate genera observed, 19 genera were not defined by 18S dinoflagellate ASVs compared to only three among the total 52 diatom genera. The discrepancy in this resolution is due to the lack of taxonomically available 18S reference sequences in particular for dinoflagellates. Overall, this preliminary investigation demonstrates that application of the full-length 18S rDNA approach can be successful in field studies.

De novo transcriptome assembly and gene annotation for the toxic dinoflagellate Dinophysis.

Species within the dinoflagellate genus Dinophysis can produce okadiac acid and dinophysistoxins leading to diarrhetic shellfish poisoning. Since the first report of D. ovum from the Gulf of Mexico in 2008, reports of other Dinophysis species across US have increased. Members of the D. cf. acuminata complex (D. acuminata, D. acuta, D. ovum, D. sacculus) are difficult to differentiate due to their morphological similarities. Dinophysis feeds on and steals the chloroplasts from the ciliate, Mesodinium rubrum, which in turn has fed on and captured the chloroplasts of its prey, the cryptophyte Teleaulax amphioxeia. The objective of this study was to generate de novo transcriptomes for new isolates of these mixotrophic organisms. The transcriptomes obtained will serve as a reference for future experiments to assess the effect of different abiotic and biotic conditions and will also provide a useful resource for screening potential marker genes to differentiate among the closely related species within the D. cf. acuminata-complex. The complete comprehensive detailed workflow and links to obtain the transcriptome data are provided.

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.

Application of a convolutional neural network to improve automated early warning of harmful algal blooms.

Continuous monitoring and early warning together represent an important mitigation strategy for harmful algal blooms (HAB). The coast of Texas experiences periodic blooms of three HAB dinoflagellates: Karenia brevis, Dinophysis ovum, and Prorocentrum texanum. A plankton image data set acquired by an Imaging FlowCytobot over a decade of operation was used to train and evaluate two new automated image classifiers. A 112 class, random forest classifier (RF_112) and a 112 class, convolutional neural network classifier (CNN_112) were developed and compared with an existing, 54 class, random forest classifier (RF_54) already in use as an early warning notification system. Both 112 class classifiers exhibited improved performance over the RF_54 classifier when tested on three different HAB species with the CNN_112 classifier producing fewer false positives and false negatives in most of the cases tested. For K. brevis and P. texanum, the current threshold of 2 cells.mL-1 was identified as the best threshold to minimize the number of false positives and false negatives. For D. ovum, a threshold of 1 cell.mL-1 was found to produce the best results with regard to the number of false positives/negatives. A lower threshold will result in earlier notification of an increase in cell concentration and will provide state health managers with increased lead time to prepare for an impending HAB.

Species detection and delineation in the marine planktonic diatoms Chaetoceros and Bacteriastrum through metabarcoding: making biological sense of haplotype diversity.

High-throughput sequencing (HTS) metabarcoding is commonly applied to assess phytoplankton diversity. Usually, haplotypes are grouped into operational taxonomic units (OTUs) through clustering, whereby the resulting number of OTUs depends on chosen similarity thresholds. We applied, instead, a phylogenetic approach to infer taxa among 18S rDNA V4-metabarcode haplotypes gathered from 48 time-series samples using the marine planktonic diatoms Chaetoceros and Bacteriastrum as test case. The 73 recovered taxa comprised both solitary haplotypes and polytomies, the latter composed each of a highly abundant, dominant haplotype and one to several minor, peripheral haplotypes. The solitary and dominant haplotypes usually matched reference sequences, enabling species assignation of taxa. We hypothesise that the super-abundance of reads in dominant haplotypes results from the homogenization effect of concerted evolution. Reads of populous peripheral haplotypes and dominant haplotypes show comparable distribution patterns over the sample dates, suggesting that they are part of the same population. Many taxa revealed marked seasonality, with closely related ones generally showing distinct periodicity, whereas others occur year-round. Phylogenies inferred from metabarcode haplotypes enable delineation of biologically meaningful taxa, whereas OTUs resulting from clustering algorithms often deviate markedly from such taxa.

Global distribution and diversity of Chaetoceros (Bacillariophyta, Mediophyceae): integration of classical and novel strategies.

Information on taxa distribution is a prerequisite for many research fields, and biological records are a major source of data contributing to biogeographic studies. The Global Biodiversity Information Facility (GBIF) and the Ocean Biogeographic Information System (OBIS) are important infrastructures facilitating free and open access to classical biological data from several sources in both temporal and spatial scales. Over the last ten years, high throughput sequencing (HTS) metabarcoding data have become available, which constitute a great source of detailed occurrence data. Among the global sampling projects that have contributed to such data are Tara Oceans and the Ocean Sampling Day (OSD). Integration of classical and metabarcoding data may aid a more comprehensive assessment of the geographic range of species, especially of microscopic ones such as protists. Rare, small and cryptic species are often ignored in surveys or mis-assigned with the classical approaches. Here we show how integration of data from various sources can contribute to insight in the biogeography and diversity at the genus- and species-level using Chaetoceros as study system, one of the most diverse and abundant genera among marine planktonic diatoms. Chaetoceros records were extracted from GBIF and OBIS and literature data were collected by means of a Google Scholar search. Chaetoceros references barcodes where mapped against the metabarcode datasets of Tara Oceans (210 sites) and OSD (144 sites). We compared the resolution of different data sources in determining the global distribution of the genus and provided examples, at the species level, of detection of cryptic species, endemism and cosmopolitan or restricted distributions. Our results highlighted at genus level a comparable picture from the different sources but a more complete assessment when data were integrated. Both the importance of the integration but also the challenges related to it were illustrated. Chaetoceros data collected in this study are organised and available in the form of tables and maps, providing a powerful tool and a baseline for further research in e.g., ecology, conservation and evolutionary biology.

Annotated 18S and 28S rDNA reference sequences of taxa in the planktonic diatom family Chaetocerotaceae.

The species-rich diatom family Chaetocerotaceae is common in the coastal marine phytoplankton worldwide where it is responsible for a substantial part of the primary production. Despite its relevance for the global cycling of carbon and silica, many species are still described only morphologically, and numerous specimens do not fit any described taxa. Nowadays, studies to assess plankton biodiversity deploy high throughput sequencing metabarcoding of the 18S rDNA V4 region, but to translate the gathered metabarcodes into biologically meaningful taxa, there is a need for reference barcodes. However, 18S reference barcodes for this important family are still relatively scarce. We provide 18S rDNA and partial 28S rDNA reference sequences of 443 morphologically characterized chaetocerotacean strains. We gathered 164 of the 216 18S sequences and 244 of the 413 28S sequences of strains from the Gulf of Naples, Atlantic France, and Chile. Inferred phylogenies showed 84 terminal taxa in seven principal clades. Two of these clades included terminal taxa whose rDNA sequences contained spliceosomal and Group IC1 introns. Regarding the commonly used metabarcode markers in planktonic diversity studies, all terminal taxa can be discriminated with the 18S V4 hypervariable region; its primers fit their targets in all but two species, and the V4-tree topology is similar to that of the 18S. Hence V4-metabarcodes of unknown Chaetocerotaceae are assignable to the family. Regarding the V9 hypervariable region, most terminal taxa can be discriminated, but several contain introns in their primer targets. Moreover, poor phylogenetic resolution of the V9 region affects placement of metabarcodes of putative but unknown chaetocerotacean taxa, and hence, uncertainty in taxonomic assignment, even of higher taxa.

Two new species in the Chaetoceros socialis complex (Bacillariophyta): C. sporotruncatus and C. dichatoensis, and characterization of its relatives, C. radicans and C. cinctus.

The diatom genus Chaetoceros is one of the most abundant and diverse phytoplankton in marine and brackish waters worldwide. Within this genus, Chaetoceros socialis has been cited as one of the most common species. However, recent studies from different geographic areas have shown the presence of pseudo-cryptic diversity within the C. socialis complex. Members of this complex are characterized by curved chains (primary colonies) aggregating into globular clusters, where one of the four setae of each cell curves toward the center of the cluster and the other three orient outwards. New light and electron microscopy observations as well as molecular data on marine planktonic diatoms from the coastal waters off Chile revealed the presence of two new species, Chaetoceros sporotruncatus sp. nov. and C. dichatoensis. sp. nov. belonging to the C. socialis complex. The two new species are similar to other members of the complex (i.e., C. socialis and C. gelidus) in the primary and secondary structure of the colony, the orientation pattern of the setae, and the valve ultrastructure. The only morphological characters that can be used to differentiate the species of this complex are aspects related to resting spore morphology. The two newly described species are closely related to each other and form a sister clade to C. gelidus in molecular phylogenies. We also provide a phylogenetic status along with the morphological characterization of C. radicans and C. cintus, which are genetically related to the C. socialis complex.

Diversity in the Globally Distributed Diatom Genus Chaetoceros (Bacillariophyceae): Three New Species from Warm-Temperate Waters.

Chaetoceros is one of the most species rich, widespread and abundant diatom genera in marine and brackish habitats worldwide. It therefore forms an excellent model for in-depth biodiversity studies, assessing morphological and genetic differentiation among groups of strains. The global Chaetoceros lorenzianus complex presently comprises three species known to science. However, our recent studies have shown that the group includes several previously unknown species. In this article, 50 strains, mainly from high latitudes and from warm-temperate waters, were examined morphologically and genetically and the results compared with those of field studies from elsewhere. The strains clustered into five groups, two of which are formed by C. decipiens Cleve and C. mitra (Bailey) Cleve, respectively. Their species descriptions are emended based on samples collected close to the type localities. The three other groups are formed by new species, C. elegans sp. nov., C. laevisporus sp. nov. and C. mannaii sp. nov. Characters used to distinguish each species are: orientation of setae, shape and size of the apertures, shape, size and density of the poroids on the setae and, at least in some species, characters of the resting spores. Our aim is to cover the global species diversity in this complex, as correct species delineation is the basis for exploring biodiversity, distribution of organisms, interactions in the food web and effects of environmental changes.

Daily variations in pathogenic bacterial populations in a monsoon influenced tropical environment.

Changing climatic conditions have influenced the monsoon pattern in recent years. Variations in bacterial population in one such tropical environment were observed everyday over two years and point out intra and inter annual changes driven by the intensity of rainfall. Vibrio spp. were abundant during the monsoon and so were faecal coliforms. Vibrio alginolyticus were negatively influenced by nitrate, whereas, silicate and rainfall positively influenced Vibrio parahaemolyticus numbers. It is also known that pathogenic bacteria are associated with the plankton. Changes in the abundance of plankton, which are governed mainly by environmental changes, could be responsible for variation in pathogenic bacterial abundance during monsoon, other than the land runoff due to precipitation and influx of fresh water.