Thermal Niche Differentiation in the Benthic Diatom Cylindrotheca closterium (Bacillariophyceae) Complex

Mixotrophic growth of a ubiquitous marine diatom

Diatoms are major players in the global carbon cycle, and their metabolism is affected by ocean conditions. Understanding the impact of changing inorganic nutrients in the oceans on diatoms is crucial, given the changes in global carbon dioxide levels. Here, we present a genome-scale metabolic model (iMK1961) for Cylindrotheca closterium, an in silico resource to understand uncharacterized metabolic functions in this ubiquitous diatom. iMK1961 represents the largest diatom metabolic model to date, comprising 1961 open reading frames and 6718 reactions. With iMK1961, we identified the metabolic response signature to cope with drastic changes in growth conditions. Comparing model predictions with Tara Oceans transcriptomics data unraveled C. closterium's metabolism in situ. Unexpectedly, the diatom only grows photoautotrophically in 21% of the sunlit ocean samples, while the majority of the samples indicate a mixotrophic (71%) or, in some cases, even a heterotrophic (8%) lifestyle in the light. Our findings highlight C. closterium's metabolic flexibility and its potential role in global carbon cycling.

Transcriptional chronology reveals conserved genes involved in pennate diatom sexual reproduction

Sexual reproduction is a major driver of adaptation and speciation in eukaryotes. In diatoms, siliceous microalgae with a unique cell size reduction-restitution life cycle and among the world's most prolific primary producers, sex also acts as the main mechanism for cell size restoration through the formation of an expanding auxospore. However, the molecular regulators of the different stages of sexual reproduction and size restoration are poorly explored. Here, we combined RNA sequencing with the assembly of a 55 Mbp reference genome for Cylindrotheca closterium to identify patterns of gene expression during different stages of sexual reproduction. These were compared with a corresponding transcriptomic time series of Seminavis robusta to assess the degree of expression conservation. Integrative orthology analysis revealed 138 one-to-one orthologues that are upregulated during sex in both species, among which 56 genes consistently upregulated during cell pairing and gametogenesis, and 11 genes induced when auxospores are present. Several early, sex-specific transcription factors and B-type cyclins were also upregulated during sex in other pennate and centric diatoms, pointing towards a conserved core regulatory machinery for meiosis and gametogenesis across diatoms. Furthermore, we find molecular evidence that the pheromone-induced cell cycle arrest is short-lived in benthic diatoms, which may be linked to their active mode of mate finding through gliding. Finally, we exploit the temporal resolution of our comparative analysis to report the first marker genes for auxospore identity called AAE1-3 ("Auxospore-Associated Expression"). Altogether, we introduce a multi-species model of the transcriptional dynamics during size restoration in diatoms and highlight conserved gene expression dynamics during different stages of sexual reproduction.

Comparison between p-distance and single-locus species delimitation models for delineating reproductively tested strains of pennate diatoms (Bacillariophyceae) using cox1, rbcL and ITS

Several automated molecular methods have emerged for distinguishing eukaryote species based on DNA sequence data. However, there are knowledge gaps around which of these single-locus methods is more accurate for the identification of microalgal species, such as the highly diverse and ecologically relevant diatoms. We applied genetic divergence, Automatic Barcode Gap Discovery for primary species delimitation (ABGD), Assemble Species by Automatic Partitioning (ASAP), Statistical Parsimony Network Analysis (SPNA), Generalized Mixed Yule Coalescent (GMYC) and Poisson Tree Processes (PTP) using partial cox1, rbcL, 5.8S + ITS2, ITS1 + 5.8S + ITS2 markers to delineate species and compare to published polyphasic identification data (morphological features, phylogeny and sexual reproductive isolation) to test the resolution of these methods. ASAP, ABGD, SPNA and PTP models resolved species of Eunotia, Seminavis, Nitzschia, Sellaphora and Pseudo-nitzschia corresponding to previous polyphasic identification, including reproductive isolation studies. In most cases, these models identified diatom species in similar ways, regardless of sequence fragment length. GMYC model presented smallest number of results that agreed with previous published identification. Following the recommendations for proper use of each model presented in the present study, these models can be useful tools to identify cryptic or closely related species of diatoms, even when the datasets have relatively few sequences.

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.

Investigation of the role of cell hydrophobicity and EPS production in the aggregation of the marine diatom Cylindrotheca closterium under hypo-saline conditions

Aggregation of diatoms is of global importance to understand settling of particulate organic carbon in aquatic systems. In this study, we investigate the aggregation of the marine diatom Cylindrotheca closterium during the exponential growth phase under hypo-saline conditions. The results of the flocculation/flotation experiments show that the aggregation of the diatom depends on the salinity. In favorable growth conditions for marine diatoms (salinity of 35), the highest aggregation is achieved. To explain these observations, we used a surface approach combining atomic force microscopy (AFM) and electrochemical methods to characterize both the cell surface properties and the structure of the extracellular polymeric substances (EPS) cell produce, and to quantify the amount of surface-active organic matter released. At a salinity of 35, the results showed that diatoms are soft, hydrophobic and release only small amounts of EPS organized into individual short fibrils. In contrast, diatoms adapt to a salinity of 5 by becoming much stiffer and more hydrophilic, producing larger amounts of EPS that structurally form an EPS network. Both adaptation responses of diatoms, the hydrophobic properties of diatoms and the release of EPS, appear to play an important role in diatom aggregation and explain the behavior observed at different salinities. This biophysical study provides important evidence allowing to get a deep insight into diatom interactions at the nanoscale, which may contribute to a better understanding of large-scale aggregation phenomena in aquatic systems.

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.

Selection constrains lottery assembly in the microbiomes of closely related diatom species

It is generally recognised that interactions between microalgae and bacteria play an important role in the functioning of marine ecosystems. In this context, increasing attention is paid to the processes that shape microalga-associated microbiomes. In recent years, conflicting evidence has been reported with respect to the relative importance of selective vs neutral processes in the assembly process. Whereas some studies report strong selection imposed by the host, others propose a more neutral, lottery-like assembly model according to which the chance of bacteria becoming part of the microbiome is proportional to their abundance in the environment and not driven by the selectional pressure created by the host. In the present study, we investigated to what degree selective vs neutral assembly processes constrain taxonomic, phylogenetic and functional variation within and between microbiomes associated with 69 isolates belonging to the Cylindrotheca closterium benthic marine diatom complex. The diatom cultures were initiated from non-axenic clonal isolates from different marine environments and geographic locations, and were then reared in a common garden (lab) environment. An important environmental imprint, likely due to in situ lottery dynamics, was apparent in the diatom microbiomes. However, microbiome assembly was also phylogenetically and functionally constrained through selective filtering related to the host microhabitat. Randomised microbiome assembly simulations revealed evidence for phylogenetic overdispersion in the observed microbiomes, reflecting an important role in the assembly process for competition between bacteria on the one hand and predominantly genetically driven differences between the hosts on the other hand. Our study thus shows that even between closely related diatom strains, host selection affects microbiome assembly, superimposing the predominantly stochastically driven recruitment process.

Pheromone Mediated Sexual Reproduction of Pennate Diatom Cylindrotheca closterium

Benthic diatoms dominate primary production in marine subtidal and intertidal environments. Their extraordinary species diversity and ecological success is thought to be linked with their predominantly heterothallic sexual reproduction. Little is known about pheromone involvement during mating of pennate diatoms. Here we describe pheromone guided mating in the coastal raphid diatom Cylindrotheca closterium. We show that the two mating types (mt⁺ and mt⁻) have distinct functions. Similar to other benthic diatoms, mt⁺ cells are searching for the mt⁻ cells to pair. To enhance mating efficiency mt⁻ exudes an attraction pheromone which we proved by establishing a novel capillary assay. Further, two more pheromones produced by mt⁻ promote the sexual events. One arrests the cell cycle progression of mt⁺ while the other induces gametogenesis of mt⁺. We suggest that C. closterium shares a functionally similar pheromone system with other pennate diatoms like Seminavis robusta and Pseudostaurosira trainorii which synchronize sexual events and mate attraction. Remarkably, we found no evidence of mt⁺ producing pheromones, which differentiates C. closterium from other pennates and suggests a less complex pheromone system in C. closterium.

Host specificity in diatom-bacteria interactions alleviates antagonistic effects

While different microalgae tend to be associated with different bacteria, it remains unclear whether such specific associations are beneficial for the microalgae. We assessed the impact of bacterial isolates, derived from various marine benthic diatoms, on the growth of several strains belonging to the Cylindrotheca closterium diatom species complex. We first tested the effect of 35 different bacterial isolates on the growth of a single C. closterium strain, and then evaluated the impact of 8 of these isolates on the growth of 6 C. closterium strains and 1 Cylindrotheca fusiformis strain. Surprisingly, most interactions were neutral to antagonistic. The interactions were highly specific, with diatom growth in the presence of specific bacteria differing between Cylindrotheca strains and species, and closely related bacteria eliciting contrasting diatom growth responses. These differences could be related to the origin of the bacterial isolates, as only isolates from foreign diatom hosts significantly reduced diatom growth, implying coadaptation between different Cylindrotheca strains and their associated bacteria. Interestingly, the antagonistic effect of a Marinobacter strain was alleviated by the presence of a microbial inoculum that was native to the diatom host, suggesting that coadapted bacteria might also benefit their host indirectly by preventing the establishment of harmful bacteria.