Ripe for reassessment: A synthesis of available molecular data for the speciose diatom family Bacillariaceae

Maximising environmental pressure-response relationship signals from diatom-based metabarcoding in rivers

DNA metabarcoding has been performed on a large number of river phytobenthos samples collected from the UK, using rbcL primers optimised for diatoms. Within this dataset the composition of non-diatom sequence reads was studied and the effect of including these in models for evaluating the nutrient gradient was assessed. Whilst many non-diatom taxonomic groups were detected, few contained the full diversity expected in riverine environments. This may be due to the performance of the current primers in characterising the wider phytobenthic community and influenced by the sampling method employed, as both were developed specifically for diatoms. Nevertheless, the study identified considerable diversity in some groups, e.g. Eustigmatophyceae and a wider distribution than previously thought for freshwater Phaeophyceae. These results offer a strong case for the benefits of metabarcoding for expanding knowledge of aquatic biodiversity in the UK and elsewhere. Many of the ASVs associated with non-diatoms showed significant pressure responses; however, models that included non-diatoms had similar predictive strength to those based on diatoms alone. Whilst limitations of the primers for assessing non-diatoms may play a role in explaining these results, the diatoms provide a strong signal along the nutrient gradient and other algae, therefore, add little unique information. We recommend that future developments should use ASVs to calculate metrics, with links to reference databases made as a final step to generate lists of taxa to support interpretation. Any further exploration of the potential of non-diatoms would benefit from access to a well-curated reference database, similar to diat.barcode. Such a database does not yet exist, and we caution against the indiscriminate use of NCBI GenBank as a taxonomic resource as many rbcL sequences deposited have not been curated.

Nuclear Transformation of the Marine Pennate Diatom Nitzschia sp. Strain NIES-4635 by Multi-Pulse Electroporation

Nitzschia is one of the largest genera of diatoms found in a range of aquatic environments, from freshwater to seawater. This genus contains evolutionarily and ecologically unique species, such as those that have lost photosynthetic capacity or those that live symbiotically in dinoflagellates. Several Nitzschia species have been used as indicators of water pollution. Recently, Nitzschia species have attracted considerable attention in the field of biotechnology. In this study, a transformation method for the marine pennate diatom Nitzschia sp. strain NIES-4635, isolated from the coastal Seto Inland Sea, was established. Plasmids containing the promoter/terminator of the fucoxanthin chlorophyll a/c binding protein gene (fcp, or Lhcf) derived from Nitzschia palea were constructed and introduced into cells by multi-pulse electroporation, resulting in 500 μg/mL nourseothricin-resistant transformants with transformation frequencies of up to 365 colonies per 108 cells. In addition, when transformation was performed using a new plasmid containing a promoter derived from a diatom-infecting virus upstream of the green fluorescent protein gene (gfp), 44% of the nourseothricin-resistant clones exhibited GFP fluorescence. The integration of the genes introduced into the genomes of the transformants was confirmed by Southern blotting. The Nitzschia transformation method established in this study will enable the transformation this species, thus allowing the functional analysis of genes from the genus Nitzschia, which are important species for environmental and biotechnological development.

Diversity and toxicity of the planktonic diatom genus Pseudo-nitzschia from coastal and offshore waters of the Southeast Pacific, including Pseudo-nitzschia dampieri sp. nov

To expand knowledge of Pseudo-nitzschia species in the Southeast Pacific, we isolated specimens from coastal waters of central Chile (36°S-30°S), the Gulf of Corcovado, and the oceanic Robinson Crusoe Island (700 km offshore) and grew them into monoclonal strains. A total of 123 Pseudo-nitzschia strains were identified to 11 species based on sequencing of the ITS region of the nuclear rDNA and on ultrastructural and morphometric analyses of the frustule in selected representatives of each clade: P. australis, P. bucculenta, P. cf. chiniana, P. cf. decipiens, P. fraudulenta, P. hasleana, P. multistriata, P. plurisecta, P. cf. sabit, the new species P. dampieri sp. nov., and one undescribed species. Partial 18S and 28S rDNA sequences, including the hypervariable V4 and D1-D3 regions used for barcoding, were gathered from representative strains of each species to facilitate future metabarcoding studies. Results showed different levels of genetic, and at times ultrastructural, diversity among the above-mentioned entities, suggesting morphological variants (P. bucculenta), rapidly radiating complexes with ill-defined species boundaries (P. cf. decipiens and P. cf. sabit), and the presence of new species (P. dampieri sp. nov., Pseudo-nitzschia sp. 1, and probably P. cf. chiniana). Domoic acid (DA) was detected in 18 out of 82 strains tested, including those of P. australis, P. plurisecta, and P. multistriata. Toxicity varied among species mostly corresponding to expectations from previous reports, with the prominent exception of P. fraudulenta; DA was not detected in any of its 10 strains tested. In conclusion, a high diversity of Pseudo-nitzschia exists in Chilean waters, particularly offshore.

Prey preference in a kleptoplastic dinoflagellate is linked to photosynthetic performance

Dinoflagellates of the family Kryptoperidiniaceae, known as "dinotoms", possess diatom-derived endosymbionts and contain individuals at three successive evolutionary stages: a transiently maintained kleptoplastic stage; a stage containing multiple permanently maintained diatom endosymbionts; and a further permanent stage containing a single diatom endosymbiont. Kleptoplastic dinotoms were discovered only recently, in Durinskia capensis; until now it has not been investigated kleptoplastic behavior and the metabolic and genetic integration of host and prey. Here, we show D. capensis is able to use various diatom species as kleptoplastids and exhibits different photosynthetic capacities depending on the diatom species. This is in contrast with the prey diatoms in their free-living stage, as there are no differences in their photosynthetic capacities. Complete photosynthesis including both the light reactions and the Calvin cycle remain active only when D. capensis feeds on its habitual associate, the "essential" diatom Nitzschia captiva. The organelles of another edible diatom, N. inconspicua, are preserved intact after ingestion by D. capensis and expresses the psbC gene of the photosynthetic light reaction, while RuBisCO gene expression is lost. Our results indicate that edible but non-essential, "supplemental" diatoms are used by D. capensis for producing ATP and NADPH, but not for carbon fixation. D. capensis has established a species-specifically designed metabolic system allowing carbon fixation to be performed only by its essential diatoms. The ability of D. capensis to ingest supplemental diatoms as kleptoplastids may be a flexible ecological strategy, to use these diatoms as "emergency supplies" while no essential diatoms are available.

Luticola tenera sp. nov. (Diadesmidaceae, Naviculales)-A New Diatom from the Soil of the State Nature Reserve "Bastak" (Jewish Autonomous Region, Russia)

Diatoms inhabit different aquatic and non-aquatic environments, including soils. The naviculoid genus Luticola is widespread in various habitats and accounts for 264 species that are only based on morphological and morphometric characteristics. These parameters can greatly vary during the life cycle, making the species very similar to each other and complicating their unambiguous identification. During a study on soil algal diversity in the Russian Far East (Jewish Autonomous Region), we isolated a strain of naviculoid diatom and examined it using an integrative approach (phylogenetic, morphological, ultrastructural data, and life cycle). Phylogenetic analyses, based on chloroplast rbcL gene data, showed affinity of the new strain with the genus Luticola. Our alga shares morphological features typical of the genus members but differs from them by having valves with a larger width and hook-shaped external proximal raphe ends deflected to the side opposite the stigma. It was revealed that the strain reproduces via two types of sexual reproduction-isogamy and cis-anisogamy. Based on these phenotypic traits, we described the new isolate as Luticola tenera sp. nov.

Craspedostaurosnazmii sp. nov., a new diatom species (Bacillariophyta) from the Turkish Coast of the Black Sea

Craspedostauros E.J. Cox is a diatom genus comprising 17 taxa reported from various regions of the world. While many species of Craspedostauros are epibiontic, the taxa have variable ecological preferences. In this study we formally describe Craspedostaurosnazmiisp. nov., an epilithic species discovered along the Turkish Black Sea Coast, based on light and scanning electron microscopy. Craspedostaurosnazmiisp. nov. is characterized by valves that are lanceolate to narrowly lanceolate, slightly constricted near the apices with uniseriate, parallel throughout the whole valve, transapical striae and and the presence of an apical silica flap. The areolae are distributed over the valve face and the mantle. The differences and similarities between C.nazmiisp. nov. and established species of Craspedostauros are discussed. Based on shape and morphometrics, the most similar species is Craspedostauroscapensis, but it is easily distinguished from C.nazmiisp. nov. by its lack of an apical silica flap.

Evolutionary dynamics of plastomes in coscinodiscophycean diatoms revealed by comparative genomics

To understand the evolution of coscinodiscophycean diatoms, plastome sequences of six coscinodiscophycean diatom species were constructed and analyzed in this study, doubling the number of constructed plastome sequences in Coscinodiscophyceae (radial centrics). The platome sizes varied substantially in Coscinodiscophyceae, ranging from 119.1 kb of Actinocyclus subtilis to 135.8 kb of Stephanopyxis turris. Plastomes in Paraliales and Stephanopyxales tended to be larger than those in Rhizosoleniales and Coscinodiacales, which were due to the expansion of the inverted repeats (IRs) and to the marked increase of the large single copy (LSC). Phylogenomic analysis indicated that Paralia and Stephanopyxis clustered tightly to form the Paraliales-Stephanopyxales complex, which was sister to the Rhizosoleniales-Coscinodiscales complex. The divergence time between Paraliales and Stephanopyxales was estimated at 85 MYA in the middle Upper Cretaceous, indicating that Paraliales and Stephanopyxales appeared later than Coscinodiacales and Rhizosoleniales according to their phylogenetic relationships. Frequent losses of housekeeping protein-coding genes (PCGs) were observed in these coscinodiscophycean plastomes, indicating that diatom plastomes showed an ongoing reduction in gene content during evolution. Two acpP genes (acpP1 and acpP2) detected in diatom plastomes were found to be originated from an early gene duplication event occurred in the common progenitor after diatom emergence, rather than multiple independent gene duplications occurring in different lineages of diatoms. The IRs in Stephanopyxis turris and Rhizosolenia fallax-imbricata exhibited a similar trend of large expansion to the small single copy (SSC) and slightly small contraction from the LSC, which eventually led to the conspicuous increase in IR size. Gene order was highly conserved in Coscinodiacales, while multiple rearrangements were observed in Rhizosoleniales and between Paraliales and Stephanopyxales. Our results greatly expanded the phylogenetic breadth in Coscinodiscophyceae and gained novel insights into the evolution of plastomes in diatoms.

Spatial fragmentation in the distribution of diatom endosymbionts from the taxonomically clarified dinophyte Kryptoperidinium triquetrum (= Kryptoperidinium foliaceum, Peridiniales)

Among the photosynthetically active dinophytes, the Kryptoperidiniaceae are unique in having a diatom as endosymbiont instead of the widely present peridinin chloroplast. Phylogenetically, it is unresolved at present how the endosymbionts are inherited, and the taxonomic identities of two iconic dinophyte names, Kryptoperidinium foliaceum and Kryptoperidinium triquetrum, are also unclear. Multiple strains were newly established from the type locality in the German Baltic Sea off Wismar and inspected using microscopy as well as molecular sequence diagnostics of both host and endosymbiont. All strains were bi-nucleate, shared the same plate formula (i.e., po, X, 4', 2a, 7'', 5c, 7s, 5''', 2'''') and exhibited a narrow and characteristically L-shaped precingular plate 7''. Within the molecular phylogeny of Bacillariaceae, endosymbionts were scattered over the tree in a highly polyphyletic pattern, even if they were gained from different strains of a single species, namely K. triquetrum. Notably, endosymbionts from the Baltic Sea show molecular sequences distinct from the Atlantic and the Mediterranean Sea, which is the first report of such a spatial fragmentation in a planktonic species of dinophytes. The two names K. foliaceum and K. triquetrum are taxonomically clarified by epitypification, with K. triquetrum having priority over its synonym K. foliaceum. Our study underlines the need of stable taxonomy for central questions in evolutionary biology.

Biology, Genetic Diversity, and Ecology of Nitzschia acidoclinata Lange-Bertalot (Bacillariophyta)

The diatom Nitzschia acidoclinata is a widespread eurybiontic alga. There is little information on its life cycle properties and cardinal points. To fill this gap, we analyzed six N. acidoclinata clones from a range of habitats in Asiatic Russia regarding their genetic diversity, morphology, morphometry, geography, and ecology. A comparison of 15 N. acidoclinata rbcL sequences sampled across its relatively wide distribution area and contrasting habitats revealed no distinct genotypes in the species. We demonstrated that the valve morphology, their length, and the sexual activity of the investigated clones varied depending on the phase of their life cycle. In this species, abrupt size reduction was observed. It was revealed that N. acidoclinata reproduced by pedogamy, and its auxosporulation was season-dependent and observed in spring and autumn only. The mating activity in our clones was detected only when the cell size was reduced to 9–22 µm in length. The available data on sexual reproduction in the genus Nitzschia suggest that neither clades nor subclades comprise pedogamous or anisogamous taxa at the same time. However, isogamy could occur in the same clade with either pedogamy or anisogamy. These data provide a fundamental basis for the development of N. acidoclinata mass cultivation and long-term maintenance in culture technologies.

Phylotranscriptomics reveals the reticulate evolutionary history of a widespread diatom species complex

In contrast to surveys based on a few genes that often provide limited taxonomic resolution, transcriptomes provide a wealth of genomic loci that can resolve relationships among taxonomically challenging lineages. Diatoms are a diverse group of aquatic microalgae that includes important bioindicator species and many such lineages. One example is Nitzschia palea, a widespread species complex with several morphologically defined taxonomic varieties, some of which are critical pollution indicators. Morphological differences among the varieties are subtle and phylogenetic studies based on a few genes fail to resolve their evolutionary relationships. We conducted morphometric and transcriptome analyses of 10 Nitzschia palea strains to resolve the relationships among strains and taxonomic varieties. Nitzschia palea was resolved into three clades, one of which corresponds to a group of strains with narrow linear-lanceolate valves. The other morphological group recovered in the shape outline analysis was not monophyletic and consisted of two clades. Gene-tree concordance analyses and phylogenetic network estimations revealed patterns of incomplete lineage sorting and gene flow between intraspecific lineages. We detected reticulated evolutionary patterns among lineages with different morphologies, resulting in a putative recent hybrid. Our study shows that phylogenomic analyses of unlinked nuclear loci, complemented with morphometrics, can resolve complex evolutionary histories of recently diverged species complexes.

Comparative Analysis of Bacillariophyceae Chloroplast Genomes Uncovers Extensive Genome Rearrangements Associated with Speciation

The Bacillariophyceae is a species-rich, ecologically significant class of Bacillariophyta. Despite their critical importance in marine ecosystems as primary producers and in the development of harmful algal blooms (HABs), taxonomic research on Bacillariophyceae species has been hindered because of their limited morphological features, plasticity of morphologies, and the low resolution of common molecular markers. Hence molecular markers with improved resolution are urgently needed. Organelle genomes, which can be constructed efficiently with the recent development of high throughput DNA sequencing technologies and the advancement of bioinformatics tools, have been proposed as super barcodes for their higher resolution for distinguishing different species and intra-species genomic variations. In this study, we tested the value of full-length chloroplast genomes (cpDNAs) as super barcodes for distinguishing diatom species, by constructing cpDNAs of 11 strains of the class Bacillariophyceae, including Nitzschia ovalis, Nitzschia traheaformis, Cylindrotheca spp., Psammodictyon constrictum, Bacillaria paxillifer, two strains of Haslea tsukamotoi, Haslea avium, Navicula arenaria, and Pleurosigma sp. Comparative analysis of cpDNAs revealed that cpDNAs were not only adequate for resolving different species, but also for enabling recognition of high levels of genome rearrangements between cpDNAs of different species, especially for species of the genera Nitzschia, Cylindrotheca, Navicula and Haslea. Additionally, comparative analysis suggested that the positioning of species in the genus Haslea should be transferred to the genus Navicula. Chloroplast genome-based evolutionary analysis suggested that the Bacillariophyceae species first appeared during the Cretaceous period and the diversity of species rose after the mass extinction about 65 Mya. This study highlighted the value of cpDNAs in research on the biodiversity and evolution of Bacillariophyceae species, and, with the construction of more cpDNAs representing additional genera, deeper insight into the biodiversity and evolutionary relationships of Bacillariophyceae species will be gained.

Diatom Biodiversity and Speciation Revealed by Comparative Analysis of Mitochondrial Genomes

Diatoms (Bacillariophyta) constitute one of the most diverse and ecologically significant groups of phytoplankton, comprising 100,000-200,000 species in three classes Bacillariophyceae, Mediophyceae, and Coscinodiscophyceae. However, due to the limited resolution of common molecular markers including 18S rDNA, 28S rDNA, ITS, rbcL, and cox1, diatom biodiversity has not been adequately ascertained. Organelle genomes including mitochondrial genomes (mtDNAs) have been proposed to be "super barcodes" for distinguishing diatom species because of their rich genomic content, and the rapid progress of DNA sequencing technologies that has made it possible to construct mtDNAs with increasing throughout and decreasing cost. Here, we constructed complete mtDNAs of 15 diatom species including five Coscinodiscophyceae species (Guinardia delicatula, Guinardia striata, Stephanopyxis turris, Paralia sulcata, and Actinocyclus sp.), four Mediophyceae species (Hemiaulus sinensis, Odontella aurita var. minima, Lithodesmioides sp., and Helicotheca tamesis), and six Bacillariophyceae species (Nitzschia ovalis, Nitzschia sp., Nitzschia traheaformis, Cylindrotheca closterium, Haslea tsukamotoi, and Pleurosigma sp.) to test the practicality of using mtDNAs as super barcodes. We found that mtDNAs have much higher resolution compared to common molecular markers as expected. Comparative analysis of mtDNAs also suggested that mtDNAs are valuable in evolutionary studies by revealing extensive genome rearrangement events with gene duplications, gene losses, and gains and losses of introns. Synteny analyses of mtDNAs uncovered high conservation among species within an order, but extensive rearrangements including translocations and/or inversions between species of different orders within a class. Duplication of cox1 was discovered for the first time in diatoms in Nitzschia traheaformis and Haslea tsukamotoi. Molecular dating analysis revealed that the three diatom classes split 100 Mya and many diatom species appeared since 50 Mya. In conclusion, more diatom mtDNAs representing different orders will play great dividends to explore biodiversity and speciation of diatoms in different ecological regions.

Domoic acid biosynthesis in the red alga Chondria armata suggests a complex evolutionary history for toxin production

Originally isolated from the red alga Chondria armata, domoic acid (DA) is best known as a potent marine neurotoxin produced by oceanic harmful algal blooms of planktonic diatoms. Sequencing efforts to date of kainoid-producing red algae have focused exclusively on a closely related molecule, kainic acid, leaving a gap in the understanding of DA biosynthesis in red algae and its evolutionary linkage to diatoms. Here, we present the phylogenetic and biochemical investigation of DA biosynthesis in C. armata. This work demonstrates the high synteny of DA biosynthetic genes between relatively distant taxonomic groups of algae and suggests a complex evolutionary history for DA biosynthesis involving gene transfer and neofunctionalization.

Domoic acid (DA), the causative agent of amnesic shellfish poisoning, is produced by select organisms within two distantly related algal clades: planktonic diatoms and red macroalgae. The biosynthetic pathway to isodomoic acid A was recently solved in the harmful algal bloom–forming diatom Pseudonitzschia multiseries, establishing the genetic basis for the global production of this potent neurotoxin. Herein, we sequenced the 507-Mb genome of Chondria armata, the red macroalgal seaweed from which DA was first isolated in the 1950s, identifying several copies of the red algal DA (rad) biosynthetic gene cluster. The rad genes are organized similarly to the diatom DA biosynthesis cluster in terms of gene synteny, including a cytochrome P450 (CYP450) enzyme critical to DA production that is notably absent in red algae that produce the simpler kainoid neurochemical, kainic acid. The biochemical characterization of the N-prenyltransferase (RadA) and kainoid synthase (RadC) enzymes support a slightly altered DA biosynthetic model in C. armata via the congener isodomoic acid B, with RadC behaving more like the homologous diatom enzyme despite higher amino acid similarity to red algal kainic acid synthesis enzymes. A phylogenetic analysis of the rad genes suggests unique origins for the red macroalgal and diatom genes in their respective hosts, with native eukaryotic CYP450 neofunctionalization combining with the horizontal gene transfer of N-prenyltransferases and kainoid synthases to establish DA production within the algal lineages.

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.

Nitzschia anatoliensis sp. nov., a cryptic diatom species from the highly alkaline Van Lake (Turkey)

In this article we describe Nitzschia anatoliensis Górecka, Gastineau & Solak sp. nov., an example of a diatom species inhabiting extreme habitats. The new species has been isolated and successfully grown from the highly alkaline Van Lake in East Turkey. The description is based on morphology (light and scanning electron microscopy), the sequencing of its organellar genomes and several molecular phylogenies. This species could easily be overlooked because of its extreme similarity to Nitzschia aurariae but molecular phylogenies indicate that they are only distantly related. Furthermore, molecular data suggest that N. anatoliensis may occur in several alkaline lakes of Asia Minor and Siberia, but was previously misidentified as Nitzschia communis. It also revealed the very close genetic proximity between N. anatoliensis and the endosymbiont of the dinotom Kryptoperidinium foliaceum, providing additional clues on what might have been the original species of diatoms to enter symbiosis.

Diploid genomic architecture of Nitzschia inconspicua, an elite biomass production diatom

A near-complete diploid nuclear genome and accompanying circular mitochondrial and chloroplast genomes have been assembled from the elite commercial diatom species Nitzschia inconspicua. The 50 Mbp haploid size of the nuclear genome is nearly double that of model diatom Phaeodactylum tricornutum, but 30% smaller than closer relative Fragilariopsis cylindrus. Diploid assembly, which was facilitated by low levels of allelic heterozygosity (2.7%), included 14 candidate chromosome pairs composed of long, syntenic contigs, covering 93% of the total assembly. Telomeric ends were capped with an unusual 12-mer, G-rich, degenerate repeat sequence. Predicted proteins were highly enriched in strain-specific marker domains associated with cell-surface adhesion, biofilm formation, and raphe system gliding motility. Expanded species-specific families of carbonic anhydrases suggest potential enhancement of carbon concentration efficiency, and duplicated glycolysis and fatty acid synthesis pathways across cytosolic and organellar compartments may enhance peak metabolic output, contributing to competitive success over other organisms in mixed cultures. The N. inconspicua genome delivers a robust new reference for future functional and transcriptomic studies to illuminate the physiology of benthic pennate diatoms and harness their unique adaptations to support commercial algae biomass and bioproduct production.