Metabolic response to a heterologous poly-3-hydroxybutyrate (PHB) pathway in Phaeodactylum tricornutum

The marine diatom Phaeodactylum tricornutum is an emerging host for metabolic engineering, but little is known about how introduced pathways are integrated into the existing metabolic framework of the host or influence transgene expression. In this study, we expressed the heterologous poly-3-hydroxybutyrate (PHB) pathway using episomal expression, which draws on the precursor acetyl coenzyme-A (AcCoA). By experimentally perturbing cultivation conditions, we gained insight into the regulation of the endogenous metabolism in transgenic lines under various environmental scenarios, as well as on alterations in AcCoA flux within the host cell. Biosynthesis of PHB led to distinct shifts in the metabolome of the host, and further analysis revealed a condition-dependent relationship between endogenous and transgenic metabolic pathways. Under N limitation, which induced a significant increase in neutral lipid content, both metabolic and transcriptomic data suggest that AcCoA was preferably shunted into the endogenous pathway for lipid biosynthesis over the transgenic PHB pathway. In contrast, supply of organic carbon in the form of glycerol supported both fatty acid and PHB biosynthesis, suggesting cross-talk between cytosolic and plastidial AcCoA precursors. This is the first study to investigate the transcriptomic and metabolomic response of diatom cell lines expressing a heterologous multi-gene pathway under different environmental conditions, providing useful insights for future engineering attempts for pathways based on the precursor AcCoA. KEY POINTS: • PHB expression had minimal effects on transcription of adjacent pathways. • N limitation favoured native lipid rather than transgenic PHB synthesis. • Glycerol addition allowed simultaneous lipid and PHB accumulation.

Unveiling the role of novel carbohydrate-binding modules in laminarin interaction of multimodular proteins from marine Bacteroidota during phytoplankton blooms

Laminarin, a β(1,3)-glucan, serves as a storage polysaccharide in marine microalgae such as diatoms. Its abundance, water solubility and simple structure make it an appealing substrate for marine bacteria. Consequently, many marine bacteria have evolved strategies to scavenge and decompose laminarin, employing carbohydrate-binding modules (CBMs) as crucial components. In this study, we characterized two previously unassigned domains as laminarin-binding CBMs in multimodular proteins from the marine bacterium Christiangramia forsetii KT0803T, thereby introducing the new laminarin-binding CBM families CBM102 and CBM103. We identified four CBM102s in a surface glycan-binding protein (SGBP) and a single CBM103 linked to a glycoside hydrolase module from family 16 (GH16_3). Our analysis revealed that both modular proteins have an elongated shape, with GH16_3 exhibiting greater flexibility than SGBP. This flexibility may aid in the recognition and/or degradation of laminarin, while the constraints in SGBP could facilitate the docking of laminarin onto the bacterial surface. Exploration of bacterial metagenome-assembled genomes (MAGs) from phytoplankton blooms in the North Sea showed that both laminarin-binding CBM families are widespread among marine Bacteroidota. The high protein abundance of CBM102- and CBM103-containing proteins during phytoplankton blooms further emphasizes their significance in marine laminarin utilization.

Metabolic engineering of omega-3 long chain polyunsaturated fatty acids in plants using different ∆6- and ∆5-desaturases co-expressed with LPCAT from the marine diatom Phaeodactylum tricornutum

Continuous research on obtaining an even more efficient production of very long-chain polyunsaturated fatty acids (VLC-PUFAs) in plants remains one of the main challenges of scientists working on plant lipids. Since crops are not able to produce these fatty acids due to the lack of necessary enzymes, genes encoding them must be introduced exogenously from native organisms producing VLC-PUFAs. In this study we reported, in tobacco leaves, the characterization of three distinct ∆6-desaturases from diatom Phaeodactylum tricornutum, fungi Rhizopus stolonifer and microalge Osterococcus tauri and two different ∆5-desaturases from P. tricornutum and single-celled saprotrophic eukaryotes Thraustochytrium sp. The in planta agroinfiltration of essential ∆6-desaturases, ∆6-elongases and ∆5-desaturases allowed for successful introduction of eicosapentaenoic acid (20:5∆5,8,11,14,17) biosynthesis pathway. However, despite the desired, targeted production of ω3-fatty acids we detected the presence of ω6-fatty acids, indicating and confirming previous results that all tested desaturases are not specifically restricted to neither ω3- nor ω6-pathway. Nevertheless, the additional co-expression of acyl-CoA:lysophosphatidylcholine acyltransferase (LPCAT) from Phaeodactylum tricornutum boosted the proportion of ω3-fatty acids in newly synthesized fatty acid pools. For the most promising genes combinations the EPA content reached at maximum 1.4% of total lipid content and 4.5% of all fatty acids accumulated in the TAG pool. Our results for the first time describe the role of LPCAT enzyme and its effectiveness in alleviating a bottleneck called 'substrate dichotomy' for improving the transgenic production of VLC-PUFAs in plants.

Biochemical characterization of acyl-CoA:diacylglycerol acyltransferase2 from the diatom Phaeodactylum tricornutum and its potential effect on LC-PUFAs biosynthesis in planta

BACKGROUND

Eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA), belonging to ω-3 long-chain polyunsaturated fatty acids (ω3-LC-PUFAs), are essential components of human diet. They are mainly supplemented by marine fish consumption, although their native producers are oleaginous microalgae. Currently, increasing demand for fish oils is insufficient to meet the entire global needs, which puts pressure on searching for the alternative solutions. One possibility may be metabolic engineering of plants with an introduced enzymatic pathway producing ω3-LC-PUFAs.

RESULT

In this study we focused on the acyl-CoA:diacylglycerol acyltransferase2b (PtDGAT2b) from the diatom Phaeodactylum tricornutum, an enzyme responsible for triacylglycerol (TAG) biosynthesis via acyl-CoA-dependent pathway. Gene encoding PtDGAT2b, incorporated into TAG-deficient yeast strain H1246, was used to confirm its activity and conduct biochemical characterization. PtDGAT2b exhibited a broad acyl-CoA preference with both di-16:0-DAG and di-18:1-DAG, whereas di-18:1-DAG was favored. The highest preference for acyl donors was observed for 16:1-, 10:0- and 12:0-CoA. PtDGAT2b also very efficiently utilized CoA-conjugated ω-3 LC-PUFAs (stearidonic acid, eicosatetraenoic acid and EPA). Additionally, verification of the potential role of PtDGAT2b in planta, through its transient expression in tobacco leaves, indicated increased TAG production with its relative amount increasing to 8%. Its co-expression with the gene combinations aimed at EPA biosynthesis led to, beside elevated TAG accumulation, efficient accumulation of EPA which constituted even 25.1% of synthesized non-native fatty acids (9.2% of all fatty acids in TAG pool).

CONCLUSIONS

This set of experiments provides a comprehensive biochemical characterization of DGAT enzyme from marine microalgae. Additionally, this study elucidates that PtDGAT2b can be used successfully in metabolic engineering of plants designed to obtain a boosted TAG level, enriched not only in ω-3 LC-PUFAs but also in medium-chain and ω-7 fatty acids.

Chromosome-level genome assembly of marine diatom Skeletonema tropicum

Skeletonema tropicum is a marine diatom of the genus Skeletonema that also includes many well-known species including S. marinoi. S. tropicum is a high temperature preferring species thriving in tropical ocean regions or temperate ocean regions during summer-autumn. However, mechanisms of ecological adaptation of S. tropicum remain poorly understood due partially to the lack of a high-quality whole genome assembly. Here, we report the first high-quality chromosome-scale genome assembly for S. tropicum, using cutting-edge technologies including PacBio single molecular sequencing and high-throughput chromatin conformation capture. The assembled genome has a size of 78.78 Mb with a scaffold N50 of 3.17 Mb, anchored to 23 pseudo-chromosomes. In total, 20,613 protein-coding genes were predicted, of which 17,757 (86.14%) genes were functionally annotated. Collinearity analysis of the genomes of S. tropicum and S. marinoi revealed that these two genomes were highly homologous. This chromosome-level genome assembly of S. tropicum provides a valuable genomic platform for comparative analysis of mechanisms of ecological adaption.

Zeaxanthin epoxidase 3 Knockout Mutants of the Model Diatom Phaeodactylum tricornutum Enable Commercial Production of the Bioactive Carotenoid Diatoxanthin

Carotenoids are pigments that have a range of functions in human health. The carotenoid diatoxanthin is suggested to have antioxidant, anti-inflammatory and chemo-preventive properties. Diatoxanthin is only produced by a few groups of microalgae, where it functions in photoprotection. Its large-scale production in microalgae is currently not feasible. In fact, rapid conversion into the inactive pigment diadinoxanthin is triggered when cells are removed from a high-intensity light source, which is the case during large-scale harvesting of microalgae biomass. Zeaxanthin epoxidase (ZEP) 2 and/or ZEP3 have been suggested to be responsible for the back-conversion of high-light accumulated diatoxanthin to diadinoxanthin in low-light in diatoms. Using CRISPR/Cas9 gene editing technology, we knocked out the ZEP2 and ZEP3 genes in the marine diatom Phaeodactylum tricornutum to investigate their role in the diadinoxanthin-diatoxanthin cycle and determine if one of the mutant strains could function as a diatoxanthin production line. Light-shift experiments proved that ZEP3 encodes the enzyme converting diatoxanthin to diadinoxanthin in low light. Loss of ZEP3 caused the high-light-accumulated diatoxanthin to be stable for several hours after the cultures had been returned to low light, suggesting that zep3 mutant strains could be suitable as commercial production lines of diatoxanthin.

Marine particle microbiomes during a spring diatom bloom contain active sulfate-reducing bacteria

Phytoplankton blooms fuel marine food webs with labile dissolved carbon and also lead to the formation of particulate organic matter composed of living and dead algal cells. These particles contribute to carbon sequestration and are sites of intense algal-bacterial interactions, providing diverse niches for microbes to thrive. We analyzed 16S and 18S ribosomal RNA gene amplicon sequences obtained from 51 time points and metaproteomes from 3 time points during a spring phytoplankton bloom in a shallow location (6-10 m depth) in the North Sea. Particulate fractions larger than 10 µm diameter were collected at near daily intervals between early March and late May in 2018. Network analysis identified two major modules representing bacteria co-occurring with diatoms and with dinoflagellates, respectively. The diatom network module included known sulfate-reducing Desulfobacterota as well as potentially sulfur-oxidizing Ectothiorhodospiraceae. Metaproteome analyses confirmed presence of key enzymes involved in dissimilatory sulfate reduction, a process known to occur in sinking particles at greater depths and in sediments. Our results indicate the presence of sufficiently anoxic niches in the particle fraction of an active phytoplankton bloom to sustain sulfate reduction, and an important role of benthic-pelagic coupling for microbiomes in shallow environments. Our findings may have implications for the understanding of algal-bacterial interactions and carbon export during blooms in shallow-water coastal areas.

Extensive intragenomic variations of the 18S rDNA V4 region in the toxigenic diatom species Pseudo-nitzschia multistriata revealed through high-throughput sequencing

Metabarcoding analysis is an effective technique for monitoring the domoic acid-producing Pseudo-nitzschia species in marine environments, uncovering high-levels of molecular diversity. However, such efforts may result in the overinterpretation of Pseudo-nitzschia species diversity, as molecular diversity not only encompasses interspecies and intraspecies diversities but also exhibits extensive intragenomic variations (IGVs). In this study, we analyzed the V4 region of the 18S rDNA of 30 strains of Pseudo-nitzschia multistriata collected from the coasts of China. The results showed that each P. multistriata strain harbored about a hundred of unique 18S rDNA V4 sequence varieties, of which each represented by a unique amplicon sequence variant (ASV). This study demonstrated the extensive degree of IGVs in P. multistriata strains, suggesting that IGVs may also present in other Pseudo-nitzschia species and other phytoplankton species. Understanding the scope and levels of IGVs is crucial for accurately interpreting the results of metabarcoding analysis.

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.

The chloroplast and mitochondrial genomes of two Gomphonema parvulum (Bacillariophyta) environmental isolates from South Carolina (United States) and Virginia (United States)

Gomphonema parvulum is a cosmopolitan freshwater diatom that is used as an indicator in water quality biomonitoring. In this study, we report the culturing of two geographically separated isolates from southeastern North America, their morphology, and the sequencing and assembly of their mitochondrial and chloroplast genomes. Morphologically, both strains fit G. parvulum sensu lato, but the frustules from a protected habitat in South Carolina were smaller than those cited in the historic data of this species from the same location as well as a second culture from Virginia. Phylogenetic analyses using the rbcL gene placed both within a clade with G. parvulum. Genetic markers, including full chloroplast and mitochondrial genomes and the nuclear small subunit rRNA gene region were assembled from each isolate. The organellar genomes of the two strains varied slightly in size due to small differences in intergenic regions with chloroplast genomes of 121,035 bp and 121,482 bp and mitochondrial genomes of 34,639 bp and 34,654 bp. The intraspecific pairwise identities of the chloroplast and mitochondrial genomes of these two isolates were 97.9% and 95.4%, respectively. Multigene phylogenetic analysis demonstrated a close relationship between G. parvulum, Gomphoneis minuta, and Didymosphenia geminata.

Succession of diversity, assembly mechanisms, and activities of the microeukaryotic community throughout Scrippsiella acuminata (Dinophyceae) bloom phases

Harmful algal bloom (HAB) is a rapidly expanding marine ecological hazard. Although numerous studies have been carried out about the ecological impact and the ecological mechanism of HAB outbreaks, few studies have comprehensively addressed the shifts of species composition, metabolic activity level, driving factors and community assembly mechanisms of microeukaryotic plankton in the course of the bloom event. To fill the gap of research, we conducted 18S ribosomal DNA and RNA sequencing during the initiation, development, sustenance and decline stages of a Scrippsiella acuminata (S. acuminata) bloom at the coastal sea of Fujian Province, China. We found that the bloom event caused a decrease in microeukaryotic plankton species diversity and increase in community homogeneity. Our results revealed that the RNA- and DNA-inferred communities were similar, but α-diversity was more dynamic in RNA- than in DNA-inferred communities. The main taxa with high projected metabolic activity (with RNA:DNA ratio as the proxy) during the bloom included dinoflagellates, Cercozoa, Chlorophyta, Protalveolata, and diatoms. The role of deterministic processes in microeukaryotic plankton community assembly increased during the bloom, but stochastic processes were always the dominant assembly mechanism throughout the bloom process. Our findings improve the understanding of temporal patterns, driving factors and assembly mechanisms underlying the microeukarytic plankton community in a dinoflagellate bloom.

Molecular physiology of Antarctic diatom natural assemblages and bloom event reveal insights into strategies contributing to their ecological success

The continental shelf of the Western Antarctic Peninsula (WAP) is a highly variable system characterized by strong cross-shelf gradients, rapid regional change, and large blooms of phytoplankton, notably diatoms. Rapid environmental changes coincide with shifts in plankton community composition and productivity, food web dynamics, and biogeochemistry. Despite the progress in identifying important environmental factors influencing plankton community composition in the WAP, the molecular basis for their survival in this oceanic region, as well as variations in species abundance, metabolism, and distribution, remains largely unresolved. Across a gradient of physicochemical parameters, we analyzed the metabolic profiles of phytoplankton as assessed through metatranscriptomic sequencing. Distinct phytoplankton communities and metabolisms closely mirrored the strong gradients in oceanographic parameters that existed from coastal to offshore regions. Diatoms were abundant in coastal, southern regions, where colder and fresher waters were conducive to a bloom of the centric diatom, Actinocyclus. Members of this genus invested heavily in growth and energy production; carbohydrate, amino acid, and nucleotide biosynthesis pathways; and coping with oxidative stress, resulting in uniquely expressed metabolic profiles compared to other diatoms. We observed strong molecular evidence for iron limitation in shelf and slope regions of the WAP, where diatoms in these regions employed iron-starvation induced proteins, a geranylgeranyl reductase, aquaporins, and urease, among other strategies, while limiting the use of iron-containing proteins. The metatranscriptomic survey performed here reveals functional differences in diatom communities and provides further insight into the environmental factors influencing the growth of diatoms and their predicted response to changes in ocean conditions.

IMPORTANCE

In the Southern Ocean, phytoplankton must cope with harsh environmental conditions such as low light and growth-limiting concentrations of the micronutrient iron. Using metratranscriptomics, we assessed the influence of oceanographic variables on the diversity of the phytoplankton community composition and on the metabolic strategies of diatoms along the Western Antarctic Peninsula, a region undergoing rapid climate change. We found that cross-shelf differences in oceanographic parameters such as temperature and variable nutrient concentrations account for most of the differences in phytoplankton community composition and metabolism. We opportunistically characterized the metabolic underpinnings of a large bloom of the centric diatom Actinocyclus in coastal waters of the WAP. Our results indicate that physicochemical differences from onshore to offshore are stronger than between southern and northern regions of the WAP; however, these trends could change in the future, resulting in poleward shifts in functional differences in diatom communities and phytoplankton blooms.

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.

Assessing the relevance of DNA metabarcoding compared to morphological identification for lake phytoplankton monitoring

Phytoplankton is a key biological group used to assess the ecological status of lakes. The classical monitoring approach relies on microscopic identification and counting of phytoplankton species, which is time-consuming and requires high taxonomic expertise. High-throughput sequencing, combined with metabarcoding, has recently demonstrated its potential as an alternative approach for plankton surveys. Several studies have confirmed the relevance of the diatom metabarcoding approach to calculate biotic indices based on species ecology. However, phytoplankton communities have not yet benefited from such validation. Here, by comparing the results obtained with the two methods (molecular and microscopic counting), we evaluated the relevance of metabarcoding approach for phytoplankton monitoring by considering different metrics: alpha diversity, taxonomic composition, community structure and a phytoplankton biotic index used to assess the trophic level of lakes. For this purpose, 55 samples were collected in four large alpine lakes (Aiguebelette, Annecy, Bourget, Geneva) during the year 2021. For each sample, a metabarcoding analysis based on two genetic markers (16S and 23S rRNA) was performed, in addition to the microscopic count. Regarding the trophic level of lakes, significant differences were found between index values obtained with the two approaches. The main hypothesis to explain these differences comes from the incompleteness, particularly at the species level, of the barcode reference library for the two genetic markers. It is therefore necessary to complete reference libraries for using such species-based biotic indices with metabarcoding data. Besides this, species richness and diversity were higher in the molecular inventories than in the microscopic ones. Moreover, despite differences in taxonomic composition of the floristic lists obtained by the two approaches, their community structures were similar. These results support the possibility of using metabarcoding for phytoplankton monitoring but in a different way. We suggest exploring alternative approaches to index development, such as a taxonomy-free approach.

Expanding the DNA editing toolbox: Novel lambda integrase variants targeting microalgal and human genome sequences

Recombinase enzymes are extremely efficient at integrating very large DNA fragments into target genomes. However, intrinsic sequence specificities curtail their use to DNA sequences with sufficient homology to endogenous target motifs. Extensive engineering is therefore required to broaden applicability and robustness. Here, we describe the directed evolution of novel lambda integrase variants capable of editing exogenous target sequences identified in the diatom Phaeodactylum tricornutum and the algae Nannochloropsis oceanica. These microorganisms hold great promise as conduits for green biomanufacturing and carbon sequestration. The evolved enzyme variants show >1000-fold switch in specificity towards the non-natural target sites when assayed in vitro. A single-copy target motif in the human genome with homology to the Nannochloropsis oceanica site can also be efficiently targeted using an engineered integrase, both in vitro and in human cells. The developed integrase variants represent useful additions to the DNA editing toolbox, with particular application for targeted genomic insertion of large DNA cargos.

Flexible B12 ecophysiology of Phaeocystis antarctica due to a fusion B12-independent methionine synthase with widespread homologues

Coastal Antarctic marine ecosystems are significant in carbon cycling because of their intense seasonal phytoplankton blooms. Southern Ocean algae are primarily limited by light and iron (Fe) and can be co-limited by cobalamin (vitamin B12). Micronutrient limitation controls productivity and shapes the composition of blooms which are typically dominated by either diatoms or the haptophyte Phaeocystis antarctica. However, the vitamin requirements and ecophysiology of the keystone species P. antarctica remain poorly characterized. Using cultures, physiological analysis, and comparative omics, we examined the response of P. antarctica to a matrix of Fe-B12 conditions. We show that P. antarctica is not auxotrophic for B12, as previously suggested, and identify mechanisms underlying its B12 response in cultures of predominantly solitary and colonial cells. A combination of proteomics and proteogenomics reveals a B12-independent methionine synthase fusion protein (MetE-fusion) that is expressed under vitamin limitation and interreplaced with the B12-dependent isoform under replete conditions. Database searches return homologues of the MetE-fusion protein in multiple Phaeocystis species and in a wide range of marine microbes, including other photosynthetic eukaryotes with polymorphic life cycles as well as bacterioplankton. Furthermore, we find MetE-fusion homologues expressed in metaproteomic and metatranscriptomic field samples in polar and more geographically widespread regions. As climate change impacts micronutrient availability in the coastal Southern Ocean, our finding that P. antarctica has a flexible B12 metabolism has implications for its relative fitness compared to B12-auxotrophic diatoms and for the detection of B12-stress in a more diverse set of marine microbes.

Mediator subunit MED8 interacts with heat shock transcription factor HSF3 to promote fucoxanthin synthesis in the diatom Phaeodactylum tricornutum

Fucoxanthin, a natural carotenoid that has substantial pharmaceutical value due to its anticancer, antioxidant, antiobesity, and antidiabetic properties, is biosynthesized from glyceraldehyde-3-phosphate (G3P) via a series of enzymatic reactions. However, our understanding of the transcriptional mechanisms involved in fucoxanthin biosynthesis remains limited. Using reverse genetics, the med8 mutant was identified based on its phenotype of reduced fucoxanthin content, and the biological functions of MED8 in fucoxanthin synthesis were characterized using approaches such as gene expression, protein subcellular localization, protein-protein interaction and chromatin immunoprecipitation assay. Gene-editing mutants of MED8 exhibited decreased fucoxanthin content as well as reduced expression levels of six key genes involved in fucoxanthin synthesis, namely DXS, PSY1, ZDS-like, CRTISO5, ZEP1, and ZEP3, when compared to the wild-type (WT) strain. Furthermore, we showed that MED8 interacts with HSF3, and genetic analysis revealed their shared involvement in the genetic pathway governing fucoxanthin synthesis. Additionally, HSF3 was required for MED8 association with the promoters of the six fucoxanthin synthesis genes. In conclusion, MED8 and HSF3 are involved in fucoxanthin synthesis by modulating the expression of the fucoxanthin synthesis genes. Our results increase the understanding of the molecular regulation mechanisms underlying fucoxanthin synthesis in the diatom P. tricornutum.

Coordinated phosphate uptake by extracellular alkaline phosphatase and solute carrier transporters in marine diatoms

Marine diatoms express genes encoding potential phosphate transporter and alkaline phosphatase (APase) under phosphate-limited (-P) condition. This indicates that diatoms use high-affinity phosphate uptake system with organic phosphate hydration. The function of molecules playing roles for Pi uptake was determined in this study. Pi uptake and APase activity of two marine diatoms, Phaeodactylum tricornutum and Thalassiosira pseudonana, were monitored during acclimation to -P condition. The transcript levels of Pi transporter were analyzed, and Pi transporters were localized with GFP tagging in diatom cells. KO mutants of plasma membrane solute carrier proteins (SLC34s) or APase were established, and their phenotype was evaluated. Some Na+ /Pi transporter candidates, SLC34s in P. tricornutum and T. pseudonana, increased transcript under -P condition. Whole-cell Pi transport was specifically stimulated by sodium ion but independent of potassium, lithium, or proton. Genome-editing KO of PtSLC34-5 and APase (Pt49678) in P. tricornutum was highly inhibitory for Pi uptake, and KO of TpSLC34-2 was also highly inhibitory for Pi uptake in T. pseudonana. SLC34s and APase were co-expressed under -P conditions in marine diatoms. SLC34s play a major role in the initial acclimation stage of diatom cells to -P condition and APase plays an increasing role in the prolonged Pi-starved condition.

Conserved cobalamin acquisition protein 1 is essential for vitamin B12 uptake in both Chlamydomonas and Phaeodactylum

Microalgae play an essential role in global net primary productivity and global biogeochemical cycling. Despite their phototrophic lifestyle, over half of algal species depend for growth on acquiring an external supply of the corrinoid vitamin B12 (cobalamin), a micronutrient produced only by a subset of prokaryotic organisms. Previous studies have identified protein components involved in vitamin B12 uptake in bacterial species and humans. However, little is known about its uptake in algae. Here, we demonstrate the essential role of a protein, cobalamin acquisition protein 1 (CBA1), in B12 uptake in Phaeodactylum tricornutum using CRISPR-Cas9 to generate targeted knockouts and in Chlamydomonas reinhardtii by insertional mutagenesis. In both cases, CBA1 knockout lines could not take up exogenous vitamin B12. Complementation of the C. reinhardtii mutants with the wild-type CBA1 gene restored B12 uptake, and regulation of CBA1 expression via a riboswitch element enabled control of the phenotype. When visualized by confocal microscopy, a YFP-fusion with C. reinhardtii CBA1 showed association with membranes. Bioinformatics analysis found that CBA1-like sequences are present in all major eukaryotic phyla. In algal taxa, the majority that encoded CBA1 also had genes for B12-dependent enzymes, suggesting CBA1 plays a conserved role. Our results thus provide insight into the molecular basis of algal B12 acquisition, a process that likely underpins many interactions in aquatic microbial communities.

Allele-dependent expression and functionality of lipid enzyme phospholipid:diacylglycerol acyltransferase affect diatom carbon storage and growth

In the acyl-CoA-independent pathway of triacylglycerol (TAG) synthesis unique to plants, fungi, and algae, TAG formation is catalyzed by the enzyme phospholipid:diacylglycerol acyltransferase (PDAT). The unique PDAT gene of the model diatom Phaeodactylum tricornutum strain CCMP2561 boasts 47 single nucleotide variants within protein coding regions of the alleles. To deepen our understanding of TAG synthesis, we observed the allele-specific expression of PDAT by the analysis of 87 published RNA-sequencing (RNA-seq) data and experimental validation. The transcription of one of the two PDAT alleles, Allele 2, could be specifically induced by decreasing nitrogen concentrations. Overexpression of Allele 2 in P. tricornutum substantially enhanced the accumulation of TAG by 44% to 74% under nutrient stress; however, overexpression of Allele 1 resulted in little increase of TAG accumulation. Interestingly, a more serious growth inhibition was observed in the PDAT Allele 1 overexpression strains compared with Allele 2 counterparts. Heterologous expression in yeast (Saccharomyces cerevisiae) showed that enzymes encoded by PDAT Allele 2 but not Allele 1 had TAG biosynthetic activity, and 7 N-terminal and 3 C-terminal amino acid variants between the 2 allele-encoded proteins substantially affected enzymatic activity. P. tricornutum PDAT, localized in the innermost chloroplast membrane, used monogalactosyldiacylglycerol and phosphatidylcholine as acyl donors as demonstrated by the increase of the 2 lipids in PDAT knockout lines, which indicated a common origin in evolution with green algal PDATs. Our study reveals unequal roles among allele-encoded PDATs in mediating carbon storage and growth in response to nitrogen stress and suggests an unsuspected strategy toward lipid and biomass improvement for biotechnological purposes.

Overexpression of PtVDL1 in Phaeodactylum tricornutum Increases Fucoxanthin Content under Red Light

Phaeodactylum tricornutum is a model diatom with significant biotechnological applications, including enhancing biomass, biofuel, and carotenoid production. Specifically, owing to the capacity of this organism to serve as a valuable source of essential raw materials for pharmaceuticals and nutraceuticals, ongoing research is actively focused on enhancing its productivity. One of the genes involved in various stages of fucoxanthin (Fx) biosynthesis, violaxanthin de-epoxidase like 1 (VDL1), has recently been identified. To validate the intracellular function of this gene and boost Fx production through overexpression, we established and examined three transgenic P. tricornutum lines characterized by elevated P. tricortunum VDL1 ( PtVDL1) expression and evaluate their cell growth and Fx productivity. These transgenic lines exhibited substantially increased PtVDL1 mRNA and protein levels compared to the wild type (WT). Notably, the enzyme substrate violaxanthin was entirely depleted and could not be detected in the transformants, whereas it remained at constant levels in the WT. Interestingly, under standard white light conditions, Fx productivity in the transformants remained unchanged; however, but after 48 h of exposure to red light, it increased by up to 15%. These results indicate that PtVDL1-overexpressing P. tricornutum has industrial potential, particularly for enhancing Fx production under red light conditions.

Declines in ice cover are accompanied by light limitation responses and community change in freshwater diatoms

The rediscovery of diatom blooms embedded within and beneath the Lake Erie ice cover (2007-2012) ignited interest in psychrophilic adaptations and winter limnology. Subsequent studies determined the vital role ice plays in winter diatom ecophysiology as diatoms partition to the underside of ice, thereby fixing their location within the photic zone. Yet, climate change has led to widespread ice decline across the Great Lakes, with Lake Erie presenting a nearly "ice-free" state in several recent winters. It has been hypothesized that the resultant turbid, isothermal water column induces light limitation amongst winter diatoms and thus serves as a competitive disadvantage. To investigate this hypothesis, we conducted a physiochemical and metatranscriptomic survey that spanned spatial, temporal, and climatic gradients of the winter Lake Erie water column (2019-2020). Our results suggest that ice-free conditions decreased planktonic diatom bloom magnitude and altered diatom community composition. Diatoms increased their expression of various photosynthetic genes and iron transporters, which suggests that the diatoms are attempting to increase their quantity of photosystems and light-harvesting components (a well-defined indicator of light limitation). We identified two gene families which serve to increase diatom fitness in the turbid ice-free water column: proton-pumping rhodopsins (a potential second means of light-driven energy acquisition) and fasciclins (a means to "raft" together to increase buoyancy and co-locate to the surface to optimize light acquisition). With large-scale climatic changes already underway, our observations provide insight into how diatoms respond to the dynamic ice conditions of today and shed light on how they will fare in a climatically altered tomorrow.

Plankton community changes during the last 124 000 years in the subarctic Bering Sea derived from sedimentary ancient DNA

Current global warming results in rising sea-water temperatures, and the loss of sea ice in Arctic and subarctic oceans impacts the community composition of primary producers with cascading effects on the food web and potentially on carbon export rates. This study analyzes metagenomic shotgun and diatom rbcL amplicon sequencing data from sedimentary ancient DNA of the subarctic western Bering Sea that records phyto- and zooplankton community changes over the last glacial-interglacial cycles, including the last interglacial period (Eemian). Our data show that interglacial and glacial plankton communities differ, with distinct Eemian and Holocene plankton communities. The generally warm Holocene period is dominated by picosized cyanobacteria and bacteria-feeding heterotrophic protists, while the Eemian period is dominated by eukaryotic picosized chlorophytes and Triparmaceae. By contrast, the glacial period is characterized by microsized phototrophic protists, including sea ice-associated diatoms in the family Bacillariaceae and co-occurring diatom-feeding crustaceous zooplankton. Our deep-time record of plankton community changes reveals a long-term decrease in phytoplankton cell size coeval with increasing temperatures, resembling community changes in the currently warming Bering Sea. The phytoplankton community in the warmer-than-present Eemian period is distinct from modern communities and limits the use of the Eemian as an analog for future climate scenarios. However, under enhanced future warming, the expected shift toward the dominance of small-sized phytoplankton and heterotrophic protists might result in an increased productivity, whereas the community's potential of carbon export will be decreased, thereby weakening the subarctic Bering Sea's function as an effective carbon sink.

Global Distribution and Diversity of Marine Parmales

Parmales (Bolidophyceae) is a minor eukaryotic phytoplankton group, sister to diatoms, which exists as two distinct forms of unicellular organisms: silicified cells and naked flagellates. Since their discovery, many field studies on Parmales have been performed; however, their global distribution has not yet been examined in detail. We herein compiled more than 3,000 marine DNA metabarcoding datasets targeting the V4 region of the 18S rRNA gene from the EukBank database. By linking this large dataset with the latest morphological and genetic information, we provide updated estimates on the diversity and distribution of Parmales in the global ocean at a fine taxonomic resolution. Parmalean amplicon sequence variants (ASVs) were detected in nearly 90% of the samples analyzed. However, the relative abundance of parmaleans in the eukaryotic community was less than 0.2% on average, and the estimated true richness of parmalean ASVs was approximately 316 ASVs, confirming their low abundance and diversity. A phylogenetic ana-lysis divided these algae into four clades, and three known morphotypes of silicified cells were classified into three different clades. The abundance of Parmales is generally high in the poles and decreases towards the tropics, and individual clades/subclades show further distinctions in their distribution. Collectively, the present results suggest clade/subclade-specific adaptation to different ecological niches.

Kaonashia insperata gen. et sp. nov., a eukaryotrophic flagellate, represents a novel major lineage of heterotrophic stramenopiles

Eukaryotrophic protists are ecologically significant and possess characteristics key to understanding the evolution of eukaryotes; however, they remain poorly studied, due partly to the complexities of maintaining predator-prey cultures. Kaonashia insperata, gen. nov., et sp. nov., is a free-swimming biflagellated eukaryotroph with a conspicuous ventral groove, a trait observed in distantly related lineages across eukaryote diversity. Di-eukaryotic (predator-prey) cultures of K. insperata with three marine algae (Isochrysis galbana, Guillardia theta, and Phaeodactylum tricornutum) were established by single-cell isolation. Growth trials showed that the studied K. insperata clone grew particularly well on G. theta, reaching a peak abundance of 1.0 × 105  ± 4.0 × 104 cells ml-1 . Small-subunit ribosomal DNA phylogenies infer that K. insperata is a stramenopile with moderate support; however, it does not fall within any well-defined phylogenetic group, including environmental sequence clades (e.g. MASTs), and its specific placement remains unresolved. Electron microscopy shows traits consistent with stramenopile affinity, including mastigonemes on the anterior flagellum and tubular mitochondrial cristae. Kaonashia insperata may represent a novel major lineage within stramenopiles, and be important for understanding the evolutionary history of the group. While heterotrophic stramenopile flagellates are considered to be predominantly bacterivorous, eukaryotrophy may be relatively widespread amongst this assemblage.

UV-C LED-induced cyclobutane pyrimidine dimer formation, lesion repair and mutagenesis in the biofilm-forming diatom, Navicula incerta

The use of ultraviolet-C (UV-C) irradiation in marine biofouling control is a relatively new and potentially disruptive technology. This study examined effects of UV-C exposure on the biofilm-forming diatom, Navicula incerta. UV-C-induced mutations were identified via Illumina HiSeq. A de novo genome was assembled from control sequences and reads from UV-C-exposed treatments were mapped to this genome, with a quantitative estimate of mutagenesis then derived from the frequency of single nucleotide polymorphisms. UV-C exposure increased cyclobutane pyrimidine dimer (CPD) abundance with a direct correlation between lesion formation and fluency. Cellular repair mechanisms gradually reduced CPDs over time, with the highest UV-C fluence treatments having the fastest repair rates. Mutation abundances were, however, negatively correlated with CPD abundance suggesting that UV-C exposure may influence lesion repair. The threshold fluence for CPD formation exceeding CPD repair was >1.27 J cm-2. Fluences >2.54 J cm-2 were predicted to inhibit repair mechanisms. While UV-C holds considerable promise for marine antifouling, diatoms are just one, albeit an important, component of marine biofouling communities. Determining fluence thresholds for other representative taxa, highlighting the most resistant, would allow UV-C treatments to be specifically tuned to target biofouling organisms, whilst limiting environmental effects and the power requirement.

Cytochrome c oxidase subunit I gene in Thalassiosira nordenskioeldii strains inhabiting in cold and warm sea waters

From the biota beneath the sea ice in Lake Saroma, which is adjacent to Sea of Okhotsk, a diatom culture of Saroma 16 was isolated. Strutted processes and a labiate process in Saroma 16 were characteristic of those in Thalassiosira nordenskioeldii. Similarity search analysis showed that the 826-bp rbcL-3P region sequence of this strain was 100% identical to multiple sequences registered as T. nordenskioeldii in a public database. The 4305-bp PCR-amplified mitochondrial cytochrome c oxidase subunit I (COI) gene (COI)-5P region of Saroma 16 included a 1060-bp open reading frame (ORF), which was interrupted by 934-bp and 2311-bp introns that included frame-shifted ORFs encoding reverse-transcriptase (RTase)-like proteins. Previous reports showed that a strain of the same species, CNS00052, originating from the East China Sea included no introns in the COI, whereas North Atlantic Ocean strains of the same species, such as CCMP992, CCMP993, and CCMP997, included a 2.3-kb intron in the same position as Saroma 16.

Targeted Gene Editing of Nuclear-Encoded Plastid Proteins in Phaeodactylum tricornutum via CRISPR/Cas9

Genome modifications in microalgae have emerged as a crucial and indispensable tool for research in fundamental and applied biology. In particular, CRISPR/Cas9 has gained significant recognition as a highly effective method for genome engineering in these photosynthetic organisms, enabling the targeted induction of mutations in specific regions of the genome. Here, we present a comprehensive protocol for generating knock-out mutants in the model diatom Phaeodactylum tricornutum using CRISPR/Cas9 by both biolistic transformation and bacterial conjugation. Our protocol outlines the step-by-step procedures and experimental conditions required to achieve successful genome editing, including the design and construction of guide RNAs, the delivery of CRISPR/Cas9 components into the algae cells, and the selection of the generated knockout mutants. Through the implementation of this protocol, researchers can harness the potential of CRISPR/Cas9 in P. tricornutum to advance the understanding of diatom biology and explore their potential applications in various fields.

Random mutagenesis of Phaeodactylum tricornutum using ultraviolet, chemical, and X-radiation demonstrates the need for temporal analysis of phenotype stability

We investigated two non-ionising mutagens in the form of ultraviolet radiation (UV) and ethyl methanosulfonate (EMS) and an ionising mutagen (X-ray) as methods to increase fucoxanthin content in the model diatom Phaeodactylum tricornutum. We implemented an ultra-high throughput method using fluorescence-activated cell sorting (FACS) and live culture spectral deconvolution for isolation and screening of potential pigment mutants, and assessed phenotype stability by measuring pigment content over 6 months using high-performance liquid chromatography (HPLC) to investigate the viability of long-term mutants. Both UV and EMS resulted in significantly higher fucoxanthin within the 6 month period after treatment, likely as a result of phenotype instability. A maximum fucoxanthin content of 135 ± 10% wild-type found in the EMS strain, a 35% increase. We found mutants generated using all methods underwent reversion to the wild-type phenotype within a 6 month time period. X-ray treatments produced a consistently unstable phenotype even at the maximum treatment of 1000 Grays, while a UV mutant and an EMS mutant reverted to wild-type after 4 months and 6 months, respectively, despite showing previously higher fucoxanthin than wild-type. This work provides new insights into key areas of microalgal biotechnology, by (i) demonstrating the use of an ionising mutagen (X-ray) on a biotechnologically relevant microalga, and by (ii) introducing temporal analysis of mutants which has substantial implications for strain creation and utility for industrial applications.

Evolutionary analysis and quality assessment of ζ-carbonic anhydrase sequences from environmental microbiome

Carbonic anhydrase (CA) is one of the most vital enzymes in living cells. This study has been performed due to the significance of this metalloenzyme for life and the novelty of some CA families like ζ-CA to evaluate evolutionary processes and quality check their sequences. In this study, bioinformatics methods revealed the presence of ζ-CA in some eukaryotic and prokaryotic microorganisms. Notably, it has not been previously reported in prokaryotes. The coexistence of β- and ζ-CAs in some microorganisms is also a novel finding as well. Also, our analysis identified several CA proteins with 6-14 amino acid intervals between histidine and cysteine in the second highly conserved motif, which can be classified as the novel ζ-CA subfamily members that emerged under the Zn deficiency of aquatic ecosystems and selection pressure in these environments. There is also a possibility that the achieved results are rooted in the contamination of samples from the environmental microbiome genome with genomes of diatom species and the occurrence of errors was observed in the DNA sequencing outcomes. Combining of all results from evolutionary analysis to quality control of ζ-CA DNA sequences is the incentive motivation to explore more the hidden aspects of ζ-CAs.

Oligotrophic waters of the Northwest Atlantic support taxonomically diverse diatom communities that are distinct from coastal waters

Diatoms are important components of the marine food web and one of the most species-rich groups of phytoplankton. The diversity and composition of diatoms in eutrophic nearshore habitats have been well documented due to the outsized influence of diatoms on coastal ecosystem functioning. In contrast, patterns of both diatom diversity and community composition in offshore oligotrophic regions where diatom biomass is low have been poorly resolved. To compare the diatom diversity and community composition in oligotrophic and eutrophic waters, diatom communities were sampled along a 1,250 km transect from the oligotrophic Sargasso Sea to the coastal waters of the northeast US shelf. Diatom community composition was determined by amplifying and sequencing the 18S rDNA V4 region. Of the 301 amplicon sequence variants (ASVs) identified along the transect, the majority (70%) were sampled exclusively from oligotrophic waters of the Gulf Stream and Sargasso Sea and included the genera Bacteriastrum, Haslea, Hemiaulus, Pseudo-nitzschia, and Nitzschia. Diatom ASV richness did not vary along the transect, indicating that the oligotrophic Sargasso Sea and Gulf Stream are occupied by a diverse diatom community. Although ASV richness was similar between oligotrophic and coastal waters, diatom community composition in these regions differed significantly and was correlated with temperature and phosphate, two environmental variables known to influence diatom metabolism and geographic distribution. In sum, oligotrophic waters of the western North Atlantic harbor diverse diatom assemblages that are distinct from coastal regions, and these open ocean diatoms warrant additional study, as they may play critical roles in oligotrophic ecosystems.

Phaeodactylum tricornutum: An established model species for diatom molecular research and an emerging chassis for algal synthetic biology

Diatoms are prominent and highly diverse microalgae in aquatic environments. Compared with other diatom species, Phaeodactylum tricornutum is an "atypical diatom" displaying three different morphotypes and lacking the usual silica shell. Despite being of limited ecological relevance, its ease of growth in the laboratory and well-known physiology, alongside the steady increase in genome-enabled information coupled with effective tools for manipulating gene expression, have meant it has gained increased recognition as a powerful experimental model for molecular research on diatoms. We here present a brief overview of how over the last 25 years P. tricornutum has contributed to the unveiling of fundamental aspects of diatom biology, while also emerging as a new tool for algal process engineering and synthetic biology.

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.

Resolving the dynamics of chrysolaminarin regulation in a marine diatom: A physiological and transcriptomic study

Diatom containing different active biological macromolecules are thought to be an excellent microbial cell factory. Phaeodactylum tricornutum, a model diatom, is a superb chassis organism accumulating chrysolaminarin with important bioactivities. However, the characteristic of chrysolaminarin accumulation and molecular mechanism of the fluctuated chrysolaminarin in diatom are still unknown. In this study, physiological data and transcriptomic analysis were carried out to clarify the mechanism involved in chrysolaminarin fluctuation. The results showed that chrysolaminarin content fluctuated, from 7.41 % dry weight (DW) to 40.01 % DW during one light/dark cycle, increase by day and decrease by night. The similar fluctuated characteristic was also observed in neutral lipid content. Genes related to the biosynthesis of chrysolaminarin and neutral lipid were up-regulated at the beginning of light-phase, explaining the accumulation of these biological macromolecules. Furthermore, genes involved in carbohydrate degradation, cell cycle, DNA replication and mitochondria-localized β-oxidation were up-regulated at the end of light phase and at the beginning of dark phase hinting an energy transition of carbohydrate to cell division during the dark period. Totally, our findings provide important information for the regulatory mechanism in the diurnal fluctuation of chrysolaminarin. It would also be of great help for the mass production of economical chrysolaminarin in marine diatom.

Genomic architecture constrains macromolecular allocation in dinoflagellates

Dinoflagellate genomes have a unique architecture that may constrain their physiological and biochemical responsiveness to environmental stressors. Here we quantified how nitrogen (N) starvation influenced macromolecular allocation and C:N:P of three photosynthetic marine dinoflagellates, representing different taxonomic classes and genome sizes. Dinoflagellates respond to nitrogen starvation by decreasing cellular nitrogen, protein and RNA content, but unlike many other eukaryotic phytoplankton examined RNA:protein is invariant. Additionally, 2 of the 3 species exhibit increases in cellular phosphorus and very little change in cellular carbon with N-starvation. As a consequence, N starvation induces moderate increases in C:N, but extreme decreases in N:P and C:P, relative to diatoms. Dinoflagellate DNA content relative to total C, N and P is much higher than similar sized diatoms, but similar to very small photosynthetic picoeukaryotes such as Ostreococcus. In aggregate these results indicate the accumulation of phosphate stores may be an important strategy employed by dinoflagellates to meet P requirements associated with the maintenance and replication of their large genomes.

Gene Duplication, Shifting Selection, and Dosage Balance of Silicon Transporter Proteins in Marine and Freshwater Diatoms

Numerous factors shape the evolution of protein-coding genes, including shifts in the strength or type of selection following gene duplications or changes in the environment. Diatoms and other silicifying organisms use a family of silicon transporters (SITs) to import dissolved silicon from the environment. Freshwaters contain higher silicon levels than oceans, and marine diatoms have more efficient uptake kinetics and less silicon in their cell walls, making them better competitors for a scarce resource. We compiled SITs from 37 diatom genomes to characterize shifts in selection following gene duplications and marine-freshwater transitions. A deep gene duplication, which coincided with a whole-genome duplication, gave rise to two gene lineages. One of them (SIT1-2) is present in multiple copies in most species and is known to actively import silicon. These SITs have evolved under strong purifying selection that was relaxed in freshwater taxa. Episodic diversifying selection was detected but not associated with gene duplications or habitat shifts. In contrast, genes in the second SIT lineage (SIT3) were present in just half the species, the result of multiple losses. Despite conservation of SIT3 in some lineages for the past 90-100 million years, repeated losses, relaxed selection, and low expression highlighted the dispensability of SIT3, consistent with a model of deterioration and eventual loss due to relaxed selection on SIT3 expression. The extensive but relatively balanced history of duplications and losses, together with paralog-specific expression patterns, suggest diatoms continuously balance gene dosage and expression dynamics to optimize silicon transport across major environmental gradients.

In situ diets of the bloom-forming dinoflagellate Noctiluca scintillans in Daya Bay

Red Noctiluca scintillans is a common heterotrophic dinoflagellate that forms blooms in temperate, subtropical, and tropical coastal ecosystems. The diet of this species plays an important role in its cell growth, development, and reproduction. Because limited gene diversity data are available regarding prey of this species, its diet in Daya Bay during a boreal winter bloom is reported using an integrated approach involving light microscopy, single cell isolation and plastid 16S rDNA cloning, and 18S rDNA V4 and V9 region amplification using isolated cells and environmental DNA as templates with high-throughput sequencing. While conventional light microscopy reveals the diet of this species to comprise Coscinodiscus sp. and Stephanopyxis turris (diatoms), copepod eggs, and detritus, plastid gene diversity identifies a diet comprising diatoms, cyanobacteria, and bacteria, and 18S rDNA high-throughput sequencing reveals a diet comprising 36 eukaryote families (primarily copepods, as well as diatoms, dinoflagellates, Ochrophyta, Haptophytes, Chordata, Cercozoans, Chlorophyta, Polychaeta, and ciliates). Dietary staples include copepods, diatoms, dinoflagellates, Ochrophyta, and Synechococcus. High copepod abundance in prey may reflect their relatively high abundance in environmental seawater. Thus, N. scintillans is generally omnivorous but prefers dominant phytoplankton taxa, including Rhizosoleniaceae, Leptocylindraceae, and Cymatosiraceae (diatoms), as well as Gonyaulacaceae (dinoflagellates). An integrated multi-disciplinary approach provides a more comprehensive picture of N. scintillans diet in Daya Bay, and an improved understanding of this species' ecological niche and trophic role in marine ecosystems.

Bioengineering of the Marine Diatom Phaeodactylum tricornutum with Cannabis Genes Enables the Production of the Cannabinoid Precursor, Olivetolic Acid

The increasing demand for novel natural compounds has prompted the exploration of innovative approaches in bioengineering. This study investigates the bioengineering potential of the marine diatom Phaeodactylum tricornutum through the introduction of cannabis genes, specifically, tetraketide synthase (TKS), and olivetolic acid cyclase (OAC), for the production of the cannabinoid precursor, olivetolic acid (OA). P. tricornutum is a promising biotechnological platform due to its fast growth rate, amenability to genetic manipulation, and ability to produce valuable compounds. Through genetic engineering techniques, we successfully integrated the cannabis genes TKS and OAC into the diatom. P. tricornutum transconjugants expressing these genes showed the production of the recombinant TKS and OAC enzymes, detected via Western blot analysis, and the production of cannabinoids precursor (OA) detected using the HPLC/UV spectrum when compared to the wild-type strain. Quantitative analysis revealed significant olivetolic acid accumulation (0.6-2.6 mg/L), demonstrating the successful integration and functionality of the heterologous genes. Furthermore, the introduction of TKS and OAC genes led to the synthesis of novel molecules, potentially expanding the repertoire of bioactive compounds accessible through diatom-based biotechnology. This study demonstrates the successful bioengineering of P. tricornutum with cannabis genes, enabling the production of OA as a precursor for cannabinoid production and the synthesis of novel molecules with potential pharmaceutical applications.

Validating a Promoter Library for Application in Plasmid-Based Diatom Genetic Engineering

While diatoms are promising synthetic biology platforms, there currently exists a limited number of validated genetic regulatory parts available for genetic engineering. The standard method for diatom transformation, nonspecific introduction of DNA into chromosomes via biolistic particle bombardment, is low throughput and suffers from clonal variability and epigenetic effects. Recent developments in diatom engineering have demonstrated that autonomously replicating episomal plasmids serve as stable expression platforms for diverse gene expression technologies. These plasmids are delivered via bacterial conjugation and, when combined with modular DNA assembly technologies, provide a flexibility and speed not possible with biolistic-mediated strain generation. In order to expand the current toolbox for plasmid-based engineering in the diatom Phaeodactylum tricornutum, a conjugation-based forward genetics screen for promoter discovery was developed, and application to a diatom genomic DNA library defined 252 P. tricornutum promoter elements. From this library, 40 promoter/terminator pairs were delivered via conjugation on episomal plasmids, characterized in vivo, and ranked across 4 orders of magnitude difference in reporter gene expression levels.

Extraction of Diatom DNA from Water Samples and Tissues

Diatom testing is an essential auxiliary means in forensic practice to determine whether the corpse drowned in water and to infer the drowning location. Diatom testing is also an important research content in the field of the environment and plankton. The diatom molecular biology testing technology, which focuses on diatom DNA as the primary research object, is a new method of diatom testing. Diatom DNA extraction is the basis of diatom molecular testing. At present, the kits commonly used for diatom DNA extraction are expensive, which increases the cost of carrying out related research. Our laboratory improved the general whole blood genomic DNA rapid extraction kit and obtained a satisfactory diatom DNA extraction effect, thus providing an alternative economical and affordable DNA extraction solution based on glass beads for related research. The diatom DNA extracted using this protocol could satisfy many downstream applications, such as PCR and sequencing.

Mapping of subcellular local pH in the marine diatom Phaeodactylum tricornutum

Diatoms are one of the most important phytoplankton on Earth. They comprise at least ten thousand species and contribute to up to 20% of the global primary production. Because of serial endosymbiotic events and horizontal gene transfers, diatoms have developed a "secondary plastid" bounded by four membranes containing a large phase-separated compartment, termed the pyrenoid. However, the physiological significance of this unique chloroplast morphology is poorly understood. Characterization of fundamental physiological parameters such as local pH in various subcellular compartments should facilitate a greater understanding of the physiological roles of the unique structure of the secondary plastid. A promising method to estimate local pH is the in situ expression of the pH-sensitive green fluorescent protein. Here, we first developed the molecular tool for the mapping of in situ local pH in the diatom Phaeodactylum tricornutum by heterologously expressing pHluorin2 in the cytosol, periplastidal compartment (PPC; the space in between two sets of outer and inner chloroplast envelopes), chloroplast stroma, and the pyrenoid matrix. Our data suggested that PPC and the pyrenoid matrix are more acidic than the adjacent areas, the cytosol and the chloroplast stroma. Finally, absolute pH values at each compartment were estimated from the ratiometric fluorescence of a recombinant pHluorin2 protein, giving pH values of approximately 7.9, 6.8, 8.0, and 7.5 respectively, for the cytosol, PPC, stroma, and pyrenoid of the P. tricornutum cells, indicating the occurrence of pH gradients and the associated electrochemical potentials at their boundary.

Seasonal dynamics of eukaryotic microbial diversity in hypersaline Tuz Lake characterized by 18S rDNA sequencing

Microbial diversity found in hypersaline ecosystems is structurally unique and essential in many microbiological and ecological processes. Tuz Lake, the second biggest lake in Türkiye, is a talassohaline (over 32% [w/v]) lake with near-neutral pH. The aim of study was to investigate the composition of the eukaryotic microbial community in Tuz Lake by 18S rDNA amplicon sequencing, as well as its relationship and change with environmental factors during 1-year period. Next-generation sequencing and bioinformatic analysis were applied to describe the eukaryotic microbial community in Tuz Lake. As a result of bioinformatics analysis, Archaeplastida (39%) and Stramenopiles, Alveolata, Rhizaria (SAR) (51%) were the most abundant taxa represented in the dataset. The Archaeplastida phylum showed a significant difference between winter and summer and higher abundance in summer in contrast to the SAR group, which represented higher abundance in winter. Genus level assessment showed that the most abundant genera were Navicula, Chlorophyta;unclassified_taxa, Dunaliella, Cladosporium, Paraphelidium, Scuticociliates;unclassified_taxa, and Chlamydomonadales;unclassified_taxa. Navicula abundance was significantly different and overwhelmingly dominant in winter. On the other hand, Cladosporium and Chlorophyta; unclassified_taxa represented a significant difference between seasons and high abundance in summer. Furthermore, Dunaliella populations were not detected in midsummer and early fall when the temperature increased and water volume in the lake decreased.

A taxonomy-free diatom eDNA-based technique for assessing lake trophic level using lake sediments

Anthropogenic eutrophication is one of the most pressing issues facing lakes globally. Our ability to manage lake eutrophication is hampered by the limited spatial and temporal extents of monitoring records, stemming from the time-consuming and expensive nature of physiochemical and biological monitoring. Diatom-based biomonitoring presents an alternative to traditional eutrophication monitoring, yet it is restricted by the high degree of taxonomic expertise required. Environmental DNA metabarcoding, while providing a promising substitute for diatom community enumeration, is plagued by inadequate taxonomic coverage of reference databases and methodological bias, limiting its use for biomonitoring. Here we show that taxonomy-free diatom-biomonitoring, in which environmental DNA metabarcoding data is utilised but not assigned to specific taxonomic classes, presents an accurate, fast, and relatively automated alternative to taxonomically assigned eutrophication biomonitoring. Our taxonomy-free index accounted for 85% of trophic level variability across 89 lakes and had the lowest average prediction error of the three approaches tested. By not relying on taxonomic identification or metabarcoding reference databases, taxonomy-free biomonitoring maintains diatom diversity that is lost in taxonomic assignment using molecular approaches. Furthermore, by utilising lake sediments, the approach outlined here presents a time-integrated estimation of lake trophic level and thus does not require time-consuming seasonal sampling. Taxonomy-free biomonitoring addresses the limitations of traditional physicochemical eutrophication monitoring and taxonomic biomonitoring alternatives and can be used to extend the spatial and temporal extents of eutrophication monitoring.

Expansion of photoreception-related gene families may drive ecological adaptation of the dominant diatom species Skeletonema marinoi

Diatom species of the genus Skeletonema are dominant in global coastal waters with important roles in marine primary production and global biogeochemical cycling. Many Skeletonema species have been extensively studied also because they can cause harmful algae blooms (HABs) with negative impacts on marine ecosystems and aquaculture. In this study, the first chromosome-level assembly of the genome of Skeletonema marinoi was constructed. The genome size was 64.99 Mb with a contig N50 of 1.95 Mb. Up to 97.12 % of contigs were successfully anchored on 24 chromosomes. Analysis of the annotated genes revealed 28 large syntenic blocks with 2397 collinear gene pairs in the genome of S. marinoi, suggesting large-scale segmental duplication events in evolution. Substantial expansion of light-harvesting genes encoding fucoxanthin-chlorophyll a/c binding proteins, as well as expansion of photoreceptor gene families encoding aureochromes and cyptochromes (CRY) in S. marinoi were found, which may have shaped ecological adaptation of S. marinoi. In conclusion, the construction of the first high-quality Skeletonema genome assembly offers valuable clues on the ecological and evolutionary characteristics of this dominant coastal diatom species.

Resolving Marine-Freshwater Transitions by Diatoms Through a Fog of Gene Tree Discordance

Despite the obstacles facing marine colonists, most lineages of aquatic organisms have colonized and diversified in freshwaters repeatedly. These transitions can trigger rapid morphological or physiological change and, on longer timescales, lead to increased rates of speciation and extinction. Diatoms are a lineage of ancestrally marine microalgae that have diversified throughout freshwater habitats worldwide. We generated a phylogenomic data set of genomes and transcriptomes for 59 diatom taxa to resolve freshwater transitions in one lineage, the Thalassiosirales. Although most parts of the species tree were consistently resolved with strong support, we had difficulties resolving a Paleocene radiation, which affected the placement of one freshwater lineage. This and other parts of the tree were characterized by high levels of gene tree discordance caused by incomplete lineage sorting and low phylogenetic signal. Despite differences in species trees inferred from concatenation versus summary methods and codons versus amino acids, traditional methods of ancestral state reconstruction supported six transitions into freshwaters, two of which led to subsequent species diversification. Evidence from gene trees, protein alignments, and diatom life history together suggest that habitat transitions were largely the product of homoplasy rather than hemiplasy, a condition where transitions occur on branches in gene trees not shared with the species tree. Nevertheless, we identified a set of putatively hemiplasious genes, many of which have been associated with shifts to low salinity, indicating that hemiplasy played a small but potentially important role in freshwater adaptation. Accounting for differences in evolutionary outcomes, in which some taxa became locked into freshwaters while others were able to return to the ocean or become salinity generalists, might help further distinguish different sources of adaptive mutation in freshwater diatoms.

Mapping harmful microalgal species by eDNA monitoring: A large-scale survey across the southwestern South China Sea

A large-scale sampling was undertaken during a research cruise across the South China Sea in August 2016, covering an area of about 100,000 km2 to investigate the molecular diversity and distributions of micro-eukaryotic protists, with a focus on the potentially harmful microalgal (HAB) species along the east coast of Peninsular Malaysia. Environmental DNAs from 30 stations were extracted and DNA metabarcoding targeting the V4 and V9 markers in the 18S rDNA was performed. Many protistan molecular units, including previously unreported HAB taxa, were discovered for the first time in the water. Our findings also revealed interesting spatial distribution patterns, with a marked signal of compositional turnover between latitudinal regimes of water masses, where dinophytes and diatom compositions were among the most strongly enhanced at the fronts, leading to distinct niches. Our results further confirmed the widespread distribution of HAB species, such as the toxigenic Alexandrium tamiyavaichii and Pseudo-nitzschia species, and the fish-killing Margalefidinium polykrikoides and Karlodinium veneficum. The molecular information obtained from this study provides an updated HAB species inventory and a toolset that could facilitate existing HAB monitoring schemes in the region to better inform management decisions.

Transcriptional responses of the marine diatom Chaetoceros tenuissimus to phosphate deficiency

The planktonic diatom Chaetoceros tenuissimus sometimes forms blooms in coastal surface waters where dissolved inorganic phosphorus (P) is typically deficient. To understand the molecular mechanisms for survival under P-deficient conditions, we compared whole transcripts and metabolites with P-sufficient conditions using stationary growth cells. Under P-deficient conditions, cell numbers and photosynthetic activities decreased as cells entered the stationary growth phase, with downregulation of transcripts related to the Calvin cycle and glycolysis/gluconeogenesis. Therefore, metabolites varied across nutritional conditions. Alkaline phosphatase, phosphodiesterase, phytase, phosphate transporter, and transcription factor genes were drastically upregulated under dissolved inorganic P deficiency. Genes related to phospholipid degradation and nonphospholipid synthesis were also upregulated. These results indicate that C. tenuissimus rearranges its membrane composition from phospholipids to nonphospholipids to conserve phosphate. To endure in P-deficient conditions, C. tenuissimus modifies its gene responses, suggesting a potential survival strategy in nature.

A chlorophyll c synthase widely co-opted by phytoplankton

Marine and terrestrial photosynthesis exhibit a schism in the accessory chlorophyll (Chl) that complements the function of Chl a: Chl b for green plants versus Chl c for most eukaryotic phytoplankton. The enzymes that mediate Chl c biosynthesis have long remained elusive. In this work, we identified the CHLC dioxygenase (Phatr3_J43737) from the marine diatom Phaeodactylum tricornutum as the Chl c synthase. The chlc mutants lacked Chl c, instead accumulating its precursors, and exhibited growth defects. In vitro, recombinant CHLC protein converted these precursors into Chl c, thereby confirming its identity. Phylogenetic evidence demonstrates conserved use of CHLC across phyla but also the existence of distinct Chl c synthases in different algal groups. Our study addresses a long-outstanding question with implications for both contemporary and ancient marine photosynthesis.

Genus-Wide Transcriptional Landscapes Reveal Correlated Gene Networks Underlying Microevolutionary Divergence in Diatoms

Marine microbes like diatoms make up the base of marine food webs and drive global nutrient cycles. Despite their key roles in ecology, biogeochemistry, and biotechnology, we have limited empirical data on how forces other than adaptation may drive diatom diversification, especially in the absence of environmental change. One key feature of diatom populations is frequent extreme reductions in population size, which can occur both in situ and ex situ as part of bloom-and-bust growth dynamics. This can drive divergence between closely related lineages, even in the absence of environmental differences. Here, we combine experimental evolution and transcriptome landscapes (t-scapes) to reveal repeated evolutionary divergence within several species of diatoms in a constant environment. We show that most of the transcriptional divergence can be captured on a reduced set of axes, and that repeatable evolution can occur along a single major axis of variation defined by core ortholog expression comprising common metabolic pathways. Previous work has associated specific transcriptional changes in gene networks with environmental factors. Here, we find that these same gene networks diverge in the absence of environmental change, suggesting these pathways may be central in generating phenotypic diversity as a result of both selective and random evolutionary forces. If this is the case, these genes and the functions they encode may represent universal axes of variation. Such axes that capture suites of interacting transcriptional changes during diversification improve our understanding of both global patterns in local adaptation and microdiversity, as well as evolutionary forces shaping algal cultivation.

Transcriptomics-guided identification of an algicidal protease of the marine bacterium Kordia algicida OT-1

In recent years, interest in algicidal bacteria has risen due to their ecological importance and their potential as biotic regulators of harmful algal blooms. Algicidal bacteria shape the plankton communities of the oceans by inhibiting or lysing microalgae and by consuming the released nutrients. Kordia algicida strain OT-1 is a model marine algicidal bacterium that was isolated from a bloom of the diatom Skeletonema costatum. Previous work has suggested that algicidal activity is mediated by secreted proteases. Here, we utilize a transcriptomics-guided approach to identify the serine protease gene KAOT1_RS09515, hereby named alpA1 as a key element in the algicidal activity of K. algicida. The protease AlpA1 was expressed and purified from a heterologous host and used in in vitro bioassays to validate its activity. We also show that K. algicida is the only algicidal species within a group of four members of the Kordia genus. The identification of this algicidal protease opens the possibility of real-time monitoring of the ecological impact of algicidal bacteria in natural phytoplankton blooms.