Comparative Transcriptomics to Identify RNA Writers and Erasers in Microalgae

Epitranscriptomics is considered as a new regulatory step in eukaryotes for developmental processes and stress responses. The aim of this study was, for the first time, to identify RNA methyltransferase (writers) and demethylase (erasers) in four investigated species, i.e., the dinoflagellates Alexandrium tamutum and Amphidinium carterae, the diatom Cylindrotheca closterium, and the green alga Tetraselmis suecica. As query sequences for the enzymatic classes of interest, we selected those ones that were previously detected in marine plants, evaluating their expression upon nutrient starvation stress exposure. The hypothesis was that upon stress exposure, the activation/deactivation of specific writers and erasers may occur. In microalgae, we found almost all plant writers and erasers (ALKBH9B, ALKBH10B, MTB, and FIP37), except for three writers (MTA, VIRILIZER, and HAKAI). A sequence similarity search by scanning the corresponding genomes confirmed their presence. Thus, we concluded that the three writer sequences were lacking from the studied transcriptomes probably because they were not expressed in those experimental conditions, rather than a real lack of these genes from their genomes. This study showed that some of them were expressed only in specific culturing conditions. We also investigated their expression in other culturing conditions (i.e., nitrogen depletion, phosphate depletion, and Zinc addition at two different concentrations) in A. carterae, giving new insights into their possible roles in regulating gene expression upon stress.

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.

First regional reference database of northern Adriatic diatom transcriptomes

Marine microbial communities form the basis for the functioning of marine ecosystems and the conservation of biodiversity. With the application of metagenomics and metatranscriptomics in marine environmental studies, significant progress has been made in analysing the functioning of microbial communities as a whole. These molecular techniques are highly dependent on reliable, well-characterised, comprehensive and taxonomically diverse sequenced reference transcriptomes of microbial organisms. Here we present a set of 12 individual transcriptome assemblies derived from 6 representative diatom species from the northern Adriatic Sea grown under 2 environmentally relevant growth conditions (phosphate replete vs. phosphate deprived). After filtering the reads and assembly, an average number of 64,932 transcripts per assembly was obtained, of which an average of 8856 were assigned to functionally known proteins. Of all assigned transcripts, an average of 6483 proteins were taxonomically assigned to diatoms (Bacillariophyta). On average, a higher number of assigned proteins was detected in the transcriptome assemblies of diatoms grown under replete media condition. On average, 50% of the mapped proteins were shared between the two growth conditions. All recorded proteins in the dataset were classified into 24 COG categories, with approximately 25% belonging to the unknown function and the remaining 75% belonging to all other categories. The resulting diatom reference database for the northern Adriatic, focussing on the response to nutrient limitation as characteristic for the region and predicted for the future world oceans, provides a valuable resource for analysing environmental metatranscriptome and metagenome data. Each northern Adriatic transcriptome can also be used by itself as a reference database for the (meta)transcriptomes and gene expression studies of the associated species that will be generated in the future.

Effects of turbulence on diatoms of the genus Pseudo-nitzschia spp. and associated bacteria

Turbulence is one of the least investigated environmental factors impacting the ecophysiology of phytoplankton, both at the community and individual species level. Here, we investigated, for the first time, the effect of a turbulence gradient (Reynolds number, from Reλ = 0 to Reλ = 360) on two species of the marine diatom Pseudo-nitzschia and their associated bacterial communities under laboratory conditions. Cell abundance, domoic acid (DA) production, chain formation, and Chl a content of P. fraudulenta and P. multiseries were higher for intermediate turbulence (Reλ = 160 or 240). DA was detectable only in P. multiseries samples. These observations were supported by transcriptomic analyses results, which suggested the turbulence related induction of the expression of the DA production locus, with a linkage to an increased photosynthetic activity of the total metatranscriptome. This study also highlighted a higher richness of the bacterial community associated with the nontoxic strain of P. fraudulenta in comparison to the toxic strain of P. multiseries. Bacillus was an important genus in P. multiseries cultures (relative abundance 15.5%) and its highest abundances coincided with the highest DA levels. However, associated bacterial communities of both Pseudo-nitzschia species did not show clear patterns relative to turbulence intensity.

Enhancing Acetate Utilization in Phaeodactylum tricornutum through the Introduction of Acetate Transport Protein

The diatom Phaeodactylum tricornutum, known for its high triacylglycerol (TAG) content and significant levels of n-3 long chain polyunsaturated fatty acids (LC-PUFAs), such as eicosapentaenoic acid (EPA), has a limited ability to utilize exogenous organic matter. This study investigates the enhancement of acetate utilization in P. tricornutum by introducing an exogenous acetate transport protein. The acetate transporter gene ADY2 from Saccharomyces cerevisiae endowed the organism with the capability to assimilate acetate and accelerating its growth. The transformants exhibited superior growth rates at an optimal NaAc concentration of 0.01 M, with a 1.7- to 2.0-fold increase compared to the wild-type. The analysis of pigments and photosynthetic activities demonstrated a decline in photosynthetic efficiency and maximum electron transport rate. This decline is speculated to result from the over-reduction of the electron transport components between photosystems due to acetate utilization. Furthermore, the study assessed the impact of acetate on the crude lipid content and fatty acid composition, revealing an increase in the crude lipid content and alterations in fatty acid profiles, particularly an increase in C16:1n-7 at the expense of EPA and a decrease in the unsaturation index. The findings provide insights into guiding the biomass and biologically active products production of P. tricornutum through metabolic engineering.

Responses to phytoplankton community succession and expression of key functional genes in plateau lakes under 17β-estradiol interference

Steroid estrogens (SEs) have garnered global attention because of their potential hazards to human health and aquatic organisms at low concentrations (ng/L). The ecosystems of plateau freshwater lakes are fragile, the water lag time is long, and pollutants easily accumulate, making them more vulnerable to the impact of SEs. However, the knowledge of the impact of SEs on the growth and decomposition of phytoplankton communities in plateau lakes and the eutrophication process is limited. This study investigated the effects and mechanisms of SEs exposure on dominant algal communities and the expression of typical algal functional genes in Erhai Lake using indoor simulations and molecular biological methods. The results showed that phytoplankton were sensitive to 17β-estradiol (E2β) pollution, with a concentration of 50, and 100 ng/L E2β exposure promoting the growth of cyanophyta and chlorophyta in the short term; this poses an ecological risk of inducing algal blooms. E2β of 1000 ng/L exposure led to cross-effects of estrogenic effects and toxicity, with most phytoplankton being inhibited. However, small filamentous cyanobacteria and diatoms exhibited greater tolerance; Melosira sp. even exhibited "low inhibition, high promotion" behavior. Exposure to E2β reduced the Shannon-Wiener diversity index (H'), Pielou index (J), and the number of dominant algal species (S) in phytoplankton communities, leading to instability in community succession. E2β of 50 ng/L enhanced the expression levels of relevant functional genes, such as ftsH, psaB, atpB, and prx, related to Microcystis aeruginosa. E2β of 50 ng/L and 5 mg/L can promote the transcription of Microcystis toxins (MC) related genes (mcyA), leading to more MC production by algal cells.

A rapid aureochrome opto-switch enables diatom acclimation to dynamic light

Diatoms often outnumber other eukaryotic algae in the oceans, especially in coastal environments characterized by frequent fluctuations in light intensity. The identities and operational mechanisms of regulatory factors governing diatom acclimation to high light stress remain largely elusive. Here, we identified the AUREO1c protein from the coastal diatom Phaeodactylum tricornutum as a crucial regulator of non-photochemical quenching (NPQ), a photoprotective mechanism that dissipates excess energy as heat. AUREO1c detects light stress using a light-oxygen-voltage (LOV) domain and directly activates the expression of target genes, including LI818 genes that encode NPQ effector proteins, via its bZIP DNA-binding domain. In comparison to a kinase-mediated pathway reported in the freshwater green alga Chlamydomonas reinhardtii, the AUREO1c pathway exhibits a faster response and enables accumulation of LI818 transcript and protein levels to comparable degrees between continuous high-light and fluctuating-light treatments. We propose that the AUREO1c-LI818 pathway contributes to the resilience of diatoms under dynamic light conditions.

The characteristics of PtHSP40 gene family in Phaeodactylum tricornutum and its response to environmental stresses

Diatom has evolved response mechanisms to cope with multiple environmental stresses. Heat shock protein 40 (HSP40) plays a key role in these response mechanisms. HSP40 gene family in higher plants has been well-studied. However, the HSP40 gene family has not been systematically investigated in marine diatom. In this study, the bioinformatic characteristics, phylogenetic relationship, conserved motifs, gene structure, chromosome distribution and the transcriptional response of PtHSP40 to different environmental stresses were analyzed in the diatom Phaeodactylum tricornutum, and quantitative real-time PCR was conducted. Totally, 55 putative PtHSP40 genes are distributed to 21 chromosomes. All PtHSP40 proteins can be divided into four groups based on their evolutionary relationship, and 54 of them contain a conserved HPD (histidine-proline-aspartic acid tripeptide) motif. Additionally, six, eleven, ten and four PtHSP40 genes were significantly upregulated under the treatments of nitrogen starvation, phosphorus deprivation, 2,2',4,4'-tetrabrominated biphenyl ether (BDE-47) and ocean acidification, respectively. More interestingly, the expression level of 9 PtHSP40 genes was obviously upregulated in response to nickel stress, suggesting the sensitive to metal stress. The different expression models of PtHSP40 genes to environmental stresses imply the specificity of PtHSP40 proteins under different stresses. This study provides a systematic understanding of the PtHSP40 gene family in P. tricornutum and a comprehensive cognition in its functions and response mechanisms to environmental stresses.

The role of cis-zeatin in enhancing high-temperature resistance and fucoxanthin biosynthesis in Phaeodactylum tricornutum

Phaeodactylum tricornutum a prominent source of industrial fucoxanthin production, faces challenges in its application due to its tolerance to high-temperature environments. This study investigates the physiological responses of P. tricornutum to high-temperature stress and its impact on fucoxanthin content, with a specific focus on the role of cis-zeatin. The results reveal that high-temperature stress inhibits P. tricornutum's growth and photosynthetic activity, leading to a decrease in fucoxanthin content. Transcriptome analysis shows that high temperature suppresses the expression of genes related to photosynthesis (e.g., psbO, psbQ, and OEC) and fucoxanthin biosynthesis (e.g., PYS, PDS1, and PSD2), underscoring the negative effects of high temperature on P. tricornutum. Interestingly, genes associated with cis-zeatin biosynthesis and cytokinesis signaling pathways exhibited increased expression under high-temperature conditions, indicating a potential role of cis-zeatin signaling in response to elevated temperatures. Content measurements confirm that high temperature enhances cis-zeatin content. Furthermore, the exogenous addition of cytokinesis mimetics or inhibitors significantly affected P. tricornutum's high-temperature resistance. Overexpression of the cis-zeatin biosynthetic enzyme gene tRNA DMATase enhanced P. tricornutum's resistance to high-temperature stress, while genetic knockout of tRNA DMATase reduced its resistance to high temperatures. Therefore, this research not only uncovers a novel mechanism for high-temperature resistance in P. tricornutum but also offers a possible alga species that can withstand high temperatures for the industrial production of fucoxanthin, offering valuable insights for practical utilization.IMPORTANCEThis study delves into Phaeodactylum tricornutum's response to high-temperature stress, specifically focusing on cis-zeatin. We uncover inhibited growth, reduced fucoxanthin, and significant cis-zeatin-related gene expression under high temperatures, highlighting potential signaling mechanisms. Crucially, genetic engineering and exogenous addition experiments confirm that the change in cis-zeatin levels could influence P. tricornutum's resistance to high-temperature stress. This breakthrough deepens our understanding of microalgae adaptation to high temperatures and offers an innovative angle for industrial fucoxanthin production. This research is a pivotal step toward developing heat-resistant microalgae for industrial use.

The draft genome of Nitzschia closterium f. minutissima and transcriptome analysis reveals novel insights into diatom biosilicification

BACKGROUND

Nitzschia closterium f. minutissima is a commonly available diatom that plays important roles in marine aquaculture. It was originally classified as Nitzschia (Bacillariaceae, Bacillariophyta) but is currently regarded as a heterotypic synonym of Phaeodactylum tricornutum. The aim of this study was to obtain the draft genome of the marine microalga N. closterium f. minutissima to understand its phylogenetic placement and evolutionary specialization. Given that the ornate hierarchical silicified cell walls (frustules) of diatoms have immense applications in nanotechnology for biomedical fields, biosensors and optoelectric devices, transcriptomic data were generated by using reference genome-based read mapping to identify significantly differentially expressed genes and elucidate the molecular processes involved in diatom biosilicification.

RESULTS

In this study, we generated 13.81 Gb of pass reads from the PromethION sequencer. The draft genome of N. closterium f. minutissima has a total length of 29.28 Mb, and contains 28 contigs with an N50 value of 1.23 Mb. The GC content was 48.55%, and approximately 18.36% of the genome assembly contained repeat sequences. Gene annotation revealed 9,132 protein-coding genes. The results of comparative genomic analysis showed that N. closterium f. minutissima was clustered as a sister lineage of Phaeodactylum tricornutum and the divergence time between them was estimated to be approximately 17.2 million years ago (Mya). CAFF analysis demonstrated that 220 gene families that significantly changed were unique to N. closterium f. minutissima and that 154 were specific to P. tricornutum, moreover, only 26 gene families overlapped between these two species. A total of 818 DEGs in response to silicon were identified in N. closterium f. minutissima through RNA sequencing, these genes are involved in various molecular processes such as transcription regulator activity. Several genes encoding proteins, including silicon transporters, heat shock factors, methyltransferases, ankyrin repeat domains, cGMP-mediated signaling pathways-related proteins, cytoskeleton-associated proteins, polyamines, glycoproteins and saturated fatty acids may contribute to the formation of frustules in N. closterium f. minutissima.

CONCLUSIONS

Here, we described a draft genome of N. closterium f. minutissima and compared it with those of eight other diatoms, which provided new insight into its evolutionary features. Transcriptome analysis to identify DEGs in response to silicon will help to elucidate the underlying molecular mechanism of diatom biosilicification in N. closterium f. minutissima.

HSP70A promotes the photosynthetic activity of marine diatom Phaeodactylum tricornutum under high temperature

With global climate change, the high-temperature environment has severely impacted the community structure and phenotype of marine diatoms. Phaeodactylum tricornutum, a model species of marine diatom, is sensitive to high temperature, which grow slowly under high temperature. However, the regulatory mechanism of P. tricornutum in response to high-temperature is still unclear. In this study, we found that the expression level of the HSP70A in the wild type (WT) increased 28 times when exposed to high temperature (26°C) for 1 h, indicating that HSP70A plays a role in high temperature in P. tricornutum. Furthermore, overexpression and interference of HSP70A have great impact on the exponential growth phase of P. tricornutum under 26°C. Moreover, the results of Co-immunoprecipitation (Co-IP) suggested that HSP70A potentially involved in the correct folding of the photosynthetic system-related proteins (D1/D2), preventing aggregation. The photosynthetic activity results demonstrated that overexpression of HSP70A improves non-photochemical quenching (NPQ) activity under high-temperature stress. These results reveal that HSP70A regulates the photosynthetic activity of P. tricornutum under high temperatures. This study not only helps us to understand the photosynthetic activity of marine diatoms to high temperature but also provides a molecular mechanism for HSP70A in P. tricornutum under high-temperature stress.

Low-CO2-inducible bestrophins outside the pyrenoid sustain high photosynthetic efficacy in diatoms

Anion transporters sustain a variety of physiological states in cells. Bestrophins (BSTs) belong to a Cl- and/or HCO3- transporter family conserved in bacteria, animals, algae, and plants. Recently, putative BSTs were found in the green alga Chlamydomonas reinhardtii, where they are upregulated under low CO2 (LC) conditions and play an essential role in the CO2-concentrating mechanism (CCM). The putative BST orthologs are also conserved in diatoms, secondary endosymbiotic algae harboring red-type plastids, but their physiological functions are unknown. Here, we characterized the subcellular localization and expression profile of BSTs in the marine diatoms Phaeodactylum tricornutum (PtBST1 to 4) and Thalassiosira pseudonana (TpBST1 and 2). PtBST1, PtBST2, and PtBST4 were localized at the stroma thylakoid membrane outside of the pyrenoid, and PtBST3 was localized in the pyrenoid. Contrarily, TpBST1 and TpBST2 were both localized in the pyrenoid. These BST proteins accumulated in cells grown in LC but not in 1% CO2 (high CO2 [HC]). To assess the physiological functions, we generated knockout mutants for the PtBST1 gene by genome editing. The lack of PtBST1 decreased photosynthetic affinity for dissolved inorganic carbon to the level comparable with the HC-grown wild type. Furthermore, non-photochemical quenching in LC-grown cells was 1.5 to 2.0 times higher in the mutants than in the wild type. These data suggest that HCO3- transport at the stroma thylakoid membranes by PtBST1 is a critical part of the CO2-evolving machinery of the pyrenoid in the fully induced CCM and that PtBST1 may modulate photoprotection under CO2-limited environments in P. tricornutum.

De novo transcriptome assembly of a lipoxygenase knock-down strain in the diatom Pseudo-nitzschia arenysensis

Diatoms are microalgae that live in marine and freshwater environments and are responsible for about 20% of the world's carbon fixation. Population dynamics of these cells is finely regulated by intricate signal transduction systems, in which oxylipins are thought to play a relevant role. These are oxygenated fatty acids whose biosynthesis is initiated by a lipoxygenase enzyme (LOX) and are widely distributed in all phyla, including diatoms. Here, we present a de novo transcriptome obtained from the RNA-seq performed in the diatom species Pseudo-nitzschia arenysensis, using both a wild-type and a LOX-silenced strain, which will represent a reliable reference for comparative analyses within the Pseudo-nitzschia genus and at a broader taxonomic scale. Moreover, the RNA-seq data can be interrogated to go deeper into the oxylipins metabolic pathways.

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.