Plastid proteome prediction for diatoms and other algae with secondary plastids of the red lineage

New plastids, old proteins: repeated endosymbiotic acquisitions in kareniacean dinoflagellates

Dinoflagellates are a diverse group of ecologically significant micro-eukaryotes that can serve as a model system for plastid symbiogenesis due to their susceptibility to plastid loss and replacement via serial endosymbiosis. Kareniaceae harbor fucoxanthin-pigmented plastids instead of the ancestral peridinin-pigmented ones and support them with a diverse range of nucleus-encoded plastid-targeted proteins originating from the haptophyte endosymbiont, dinoflagellate host, and/or lateral gene transfers (LGT). Here, we present predicted plastid proteomes from seven distantly related kareniaceans in three genera (Karenia, Karlodinium, and Takayama) and analyze their evolutionary patterns using automated tree building and sorting. We project a relatively limited ( ~ 10%) haptophyte signal pointing towards a shared origin in the family Chrysochromulinaceae. Our data establish significant variations in the functional distributions of these signals, emphasizing the importance of micro-evolutionary processes in shaping the chimeric proteomes. Analysis of plastid genome sequences recontextualizes these results by a striking finding the extant kareniacean plastids are in fact not all of the same origin, as two of the studied species (Karlodinium armiger, Takayama helix) possess plastids from different haptophyte orders than the rest.

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.

Proteomes reveal the lipid metabolic network in the complex plastid of Phaeodactylum tricornutum

Phaeodactylum tricornutum plastid is surrounded by four membranes, and its protein composition and function remain mysterious. In this study, the P. tricornutum plastid-enriched fraction was obtained and 2850 proteins were identified, including 92 plastid-encoded proteins, through label-free quantitative proteomic technology. Among them, 839 nuclear-encoded proteins were further determined to be plastidial proteins based on the BLAST alignments within Plant Proteome DataBase and subcellular localization prediction, in spite of the strong contamination by mitochondria-encoded proteins and putative plasma membrane proteins. According to our proteomic data, we reconstructed the metabolic pathways and highlighted the hybrid nature of this diatom plastid. Triacylglycerol (TAG) hydrolysis and glycolysis, as well as photosynthesis, glycan metabolism, and tocopherol and triterpene biosynthesis, occur in the plastid. In addition, the synthesis of long-chain acyl-CoAs, elongation, and desaturation of fatty acids (FAs), and synthesis of lipids including TAG are confined in the four-layered-membrane plastid based on the proteomic and GFP-fusion localization data. The whole process of generation of docosahexaenoic acid (22:6) from palmitic acid (16:0), via elongation and desaturation of FAs, occurs in the chloroplast endoplasmic reticulum membrane, the outermost membrane of the plastid. Desaturation that generates 16:4 from 16:0 occurs in the plastid stroma and outer envelope membrane. Quantitative analysis of glycerolipids between whole cells and isolated plastids shows similar composition, and the FA profile of TAG was not different. This study shows that the diatom plastid combines functions usually separated in photosynthetic eukaryotes, and differs from green alga and plant chloroplasts by undertaking the whole process of lipid biosynthesis.

Complex Endosymbioses I: From Primary to Complex Plastids, Serial Endosymbiotic Events

A considerable part of the diversity of eukaryotic phototrophs consists of algae with plastids that evolved from endosymbioses between two eukaryotes. These complex plastids are characterized by a high number of envelope membranes (more than two) and some of them contain a residual nucleus of the endosymbiotic alga called a nucleomorph. Complex plastid-bearing algae are thus chimeric cell assemblies, eukaryotic symbionts living in a eukaryotic host. In contrast, the primary plastids of the Archaeplastida (plants, green algae, red algae, and glaucophytes) possibly evolved from a single endosymbiosis with a cyanobacterium and are surrounded by two membranes. Complex plastids have been acquired several times by unrelated groups of eukaryotic heterotrophic hosts, suggesting that complex plastids are somewhat easier to obtain than primary plastids. Evidence suggests that complex plastids arose twice independently in the green lineage (euglenophytes and chlorarachniophytes) through secondary endosymbiosis, and four times in the red lineage, first through secondary endosymbiosis in cryptophytes, then by higher-order events in stramenopiles, alveolates, and haptophytes. Engulfment of primary and complex plastid-containing algae by eukaryotic hosts (secondary, tertiary, and higher-order endosymbioses) is also responsible for numerous plastid replacements in dinoflagellates. Plastid endosymbiosis is accompanied by massive gene transfer from the endosymbiont to the host nucleus and cell adaptation of both endosymbiotic partners, which is related to the trophic switch to phototrophy and loss of autonomy of the endosymbiont. Such a process is essential for the metabolic integration and division control of the endosymbiont in the host. Although photosynthesis is the main advantage of acquiring plastids, loss of photosynthesis often occurs in algae with complex plastids. This chapter summarizes the essential knowledge of the acquisition, evolution, and function of complex plastids.

Circadian rhythms and circadian clock gene homologs of complex alga Chromera velia

Most organisms on Earth are affected by periodic changes in their environment. The circadian clock is an endogenous device that synchronizes behavior, physiology, or biochemical processes to an approximately 24-hour cycle, allowing organisms to anticipate the periodic changes of day and night. Although circadian clocks are widespread in organisms, the actual molecular components differ remarkably among the clocks of plants, animals, fungi, and prokaryotes. Chromera velia is the closest known photosynthetic relative of apicomplexan parasites. Formation of its motile stage, zoospores, has been described as associated with the light part of the day. We examined the effects on the periodic release of the zoospores under different light conditions and investigated the influence of the spectral composition on zoosporogenesis. We performed a genomic search for homologs of known circadian clock genes. Our results demonstrate the presence of an almost 24-hour free-running cycle of zoosporogenesis. We also identified the blue light spectra as the essential compound for zoosporogenesis. Further, we developed a new and effective method for zoospore separation from the culture and estimated the average motility speed and lifespan of the C. velia zoospores. Our genomic search identified six cryptochrome-like genes, two genes possibly related to Arabidopsis thaliana CCA/LHY, whereas no homolog of an animal, cyanobacterial, or fungal circadian clock gene was found. Our results suggest that C. velia has a functional circadian clock, probably based mainly on a yet undefined mechanism.

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.

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.

Subcellular proteomics for determining iron-limited remodeling of plastids in the model diatom Thalassiosira pseudonana (Bacillariophyta)

Diatoms are important primary producers in the world's oceans, yet their growth is constrained in large regions by low bioavailable iron (Fe). Low-Fe stress-induced limitation of primary production is due to requirements for Fe in components of essential metabolic pathways including photosynthesis and other chloroplast plastid functions. Studies have shown that under low-Fe stress, diatoms alter plastid-specific processes, including components of electron transport. These physiological changes suggest changes of protein content and in protein abundances within the diatom plastid. While in silico predictions provide putative information on plastid-localized proteins, knowledge of diatom plastid proteins remains limited in comparison to well-studied model photosynthetic organisms. To address this, we employed shotgun proteomics to investigate the proteome of subcellular plastid-enriched fractions from Thalassiosira pseudonana to gain a better understanding of how the plastid proteome is remodeled in response to Fe limitation. Using mass spectrometry-based peptide identification and quantification, we analyzed T. pseudonana grown under Fe-replete and -limiting conditions. Through these analyses, we inferred the relative quantities of each protein, revealing that Fe limitation regulates major metabolic pathways in the plastid, including the Calvin cycle. Additionally, we observed changes in the expression of light-harvesting proteins. In silico localization predictions of proteins identified in this plastid-enriched proteome allowed for an in-depth comparison of theoretical versus observed plastid-localization, providing evidence for the potential of additional protein import pathways into the diatom plastid.

Genome analysis of Parmales, the sister group of diatoms, reveals the evolutionary specialization of diatoms from phago-mixotrophs to photoautotrophs

The order Parmales (class Bolidophyceae) is a minor group of pico-sized eukaryotic marine phytoplankton that contains species with cells surrounded by silica plates. Previous studies revealed that Parmales is a member of ochrophytes and sister to diatoms (phylum Bacillariophyta), the most successful phytoplankton group in the modern ocean. Therefore, parmalean genomes can serve as a reference to elucidate both the evolutionary events that differentiated these two lineages and the genomic basis for the ecological success of diatoms vs. the more cryptic lifestyle of parmaleans. Here, we compare the genomes of eight parmaleans and five diatoms to explore their physiological and evolutionary differences. Parmaleans are predicted to be phago-mixotrophs. By contrast, diatoms have lost genes related to phagocytosis, indicating the ecological specialization from phago-mixotrophy to photoautotrophy in their early evolution. Furthermore, diatoms show significant enrichment in gene sets involved in nutrient uptake and metabolism, including iron and silica, in comparison with parmaleans. Overall, our results suggest a strong evolutionary link between the loss of phago-mixotrophy and specialization to a silicified photoautotrophic life stage early in diatom evolution after diverging from the Parmales lineage.

Characterization of oil body-associated proteins obtained from oil bodies with different sizes in oleaginous diatom Fistulifera solaris

Oil body-associated proteins from the oleaginous diatom Fistulifera solaris were identified by proteomic analysis of oil bodies of various sizes (small, middle, and large) by time-dependent culturing upon nutrient-starvation at 36, 96 and 168 h. This diatom strain has the capability to accumulate neutral lipids and triacylglycerol. Liquid chromatography-tandem mass spectrometry analysis revealed 662 proteins in all oil body sizes. Among these, 132 proteins were predicted to be localized to the endoplasmic reticulum. Seventeen proteins that exhibited a positive correlation with gene expression and the oil body size were selected as novel candidates for oil body-associated proteins. Among the 17 protein candidates, two proteins encoded by fso:g8246 and fso:g10200 were confirmed to be localized on the surface of the oil body and endoplasmic reticulum. A protein encoded by fso:g2514, which is involved in sterol biosynthesis, was also identified. This protein was likely to localize to mitochondria; however, inhibitor assays suggested that it might play a role in lipid degradation. Our work provides new insights into the proteomics of microalgae and provides a valuable strategy for boosting lipid productivity in microalgae.

Bioinformatics-Based Screening Approach for the Identification and Characterization of Lipolytic Enzymes from the Marine Diatom Phaeodactylum tricornutum

Oleaginous diatoms accumulate lipids of biotechnological interest when exposed to nutrient stress conditions such as nitrogen starvation. While accumulation mechanisms are well-known and have been engineered to improve lipid production, degradation mechanisms remain poorly investigated in diatoms. Identifying lipid-degrading enzymes is the initial step to understanding the catabolic processes. In this study, an in silico screening of the genome of Phaeodactylum tricornutum led to the identification of 57 putative triacylglycerol lipases (EC 3.1.1.3) grouped in 4 families. Further analysis revealed the presence of conserved domains and catalytic residues of lipases. Physico-chemical characteristics and subcellular localization predictions highlighted that a majority of these putative proteins are hydrophilic and cytosolic, suggesting they could be recruited to lipid droplets directly from the cytosol. Among the 57 identified putative proteins, three lipases were identified as possibly involved in lipophagy due to a potential vacuolar localization. The expression of the mRNA corresponding to the 57 proteins was then searched in 3 transcriptomic datasets obtained under nitrogen starvation. Nine genes were highly regulated and were considered as encoding enzymes with a probable important function in lipid catabolism. A tertiary structure prediction of these nine candidates yielded eight functional 3D models. Among those, two downregulated enzymes, Phatr3_J54974 and Phatr3_EG00720, were highlighted as good targets for future functional genomics and purification studies to investigate their role in lipid degradation.

In vivo localization of iron starvation induced proteins under variable iron supplementation regimes in Phaeodactylum tricornutum

The model pennate diatom Phaeodactylum tricornutum is able to assimilate a range of iron sources. It therefore provides a platform to study different mechanisms of iron processing concomitantly in the same cell. In this study, we follow the localization of three iron starvation induced proteins (ISIPs) in vivo, driven by their native promoters and tagged by fluorophores in an engineered line of P. tricornutum. We find that the localization patterns of ISIPs are dynamic and variable depending on the overall iron status of the cell and the source of iron it is exposed to. Notwithstanding, a shared destination of the three ISIPs both under ferric iron and siderophore-bound iron supplementation is a globular compartment in the vicinity of the chloroplast. In a proteomic analysis, we identify that the cell engages endocytosis machinery involved in the vesicular trafficking as a response to siderophore molecules, even when these are not bound to iron. Our results suggest that there may be a direct vesicle traffic connection between the diatom cell membrane and the periplastidial compartment (PPC) that co-opts clathrin-mediated endocytosis and the "cytoplasm to vacuole" (Cvt) pathway, for proteins involved in iron assimilation. Proteomics data are available via ProteomeXchange with identifier PXD021172.

HIGHLIGHT

The marine diatom P. tricornutum engages a vesicular network to traffic siderophores and phytotransferrin from the cell membrane directly to a putative iron processing site in the vicinity of the chloroplast.

Expression of nuclear-encoded, haptophyte-derived ftsH genes support extremely rapid PSII repair and high-light photoacclimation in Karenia brevis (Dinophyceae)

Karenia brevis, a neurotoxic dinoflagellate that produces brevetoxins, is endemic to the Gulf of Mexico and can grow at high irradiances typical of surface waters found there. To build upon a growing number of studies addressing high-light tolerance in K. brevis, specific photobiology and molecular mechanisms underlying this capacity were evaluated in culture. Since photosystem II (PSII) repair cycle activity can be crucial to high light tolerance in plants and algae, the present study assessed this capacity in K. brevis and characterized the ftsH-like genes which are fundamental to this process. Compared with cultures grown in low-light, cultures grown in high-light showed a 65-fold increase in PSII photoinactivation, a ∼50-fold increase in PSII repair, enhanced nonphotochemical quenching (NPQ), and depressed F_v/F_m. Repair rates were among the fastest reported in phytoplankton. Publicly available K. brevis transcriptomes (MMETSP) were queried for ftsH-like sequences and refined with additional sequencing from two K. brevis strains. The genes were phylogenetically related to haptophyte orthologs, implicating acquisition during tertiary endosymbiosis. RT-qPCR of three of the four ftsH-like homologs revealed that poly-A tails predominated in all homologs, and that the most highly expressed homolog had a 5' splice leader and amino-acid motifs characteristic of chloroplast targeting, indicating nuclear encoding for this plastid-targeted gene. High-light cultures showed a ∼1.5-fold upregulation in mRNA expression of the thylakoid-associated genes. Overall, in conjunction with NPQ mechanisms, rapid PSII repair mediated by a haptophyte-derived ftsH prevents chronic photoinhibition in K. brevis. Our findings continue to build the case that high-light photobiology-supported by the acquisition and maintenance of tertiary endosymbiotic genes-is critical to the success of K. brevis in the Gulf of Mexico.

A metabolic, phylogenomic and environmental atlas of diatom plastid transporters from the model species Phaeodactylum

Diatoms are an important group of algae, contributing nearly 40% of total marine photosynthetic activity. However, the specific molecular agents and transporters underpinning the metabolic efficiency of the diatom plastid remain to be revealed. We performed in silico analyses of 70 predicted plastid transporters identified by genome-wide searches of Phaeodactylum tricornutum. We considered similarity with Arabidopsis thaliana plastid transporters, transcriptional co-regulation with genes encoding core plastid metabolic pathways and with genes encoded in the mitochondrial genomes, inferred evolutionary histories using single-gene phylogeny, and environmental expression trends using Tara Oceans meta-transcriptomics and meta-genomes data. Our data reveal diatoms conserve some of the ion, nucleotide and sugar plastid transporters associated with plants, such as non-specific triose phosphate transporters implicated in the transport of phosphorylated sugars, NTP/NDP and cation exchange transporters. However, our data also highlight the presence of diatom-specific transporter functions, such as carbon and amino acid transporters implicated in intricate plastid-mitochondria crosstalk events. These confirm previous observations that substrate non-specific triose phosphate transporters (TPT) may exist as principal transporters of phosphorylated sugars into and out of the diatom plastid, alongside suggesting probable agents of NTP exchange. Carbon and amino acid transport may be related to intricate metabolic plastid-mitochondria crosstalk. We additionally provide evidence from environmental meta-transcriptomic/meta- genomic data that plastid transporters may underpin diatom sensitivity to ocean warming, and identify a diatom plastid transporter (J43171) whose expression may be positively correlated with temperature.

Trypsin is a coordinate regulator of N and P nutrients in marine phytoplankton

Trypsin is best known as a digestive enzyme in animals, but remains unexplored in phytoplankton, the major primary producers in the ocean. Here we report the prevalence of trypsin genes in global ocean phytoplankton and significant influences of environmental nitrogen (N) and phosphorus (P) on their expression. Using CRISPR/Cas9 mediated-knockout and overexpression analyses, we further reveal that a trypsin in Phaeodactylum tricornutum (PtTryp2) functions to repress N acquisition, but its expression decreases under N-deficiency to promote N acquisition. On the contrary, PtTryp2 promotes phosphate uptake per se, and its expression increases under P-deficiency to further reinforce P acquisition. Furthermore, PtTryp2 knockout led to amplitude magnification of the nitrate and phosphate uptake ‘seesaw’, whereas PtTryp2 overexpression dampened it, linking PtTryp2 to stabilizing N:P stoichiometry. Our data demonstrate that PtTryp2 is a coordinate regulator of N:P stoichiometric homeostasis. The study opens a window for deciphering how phytoplankton adapt to nutrient-variable marine environments.

Using CRISPR-Cas9 mediated-knockout and overexpression analyses, this study shows that a trypsin in the diatom Phaeodactylum tricornutum promotes phosphorus uptake and inhibits nitrogen uptake but its expression is downregulated under nitrogen stress and upregulated under phosphorus stress. Together, the findings suggest this trypsin is a coordinate regulator of nutrient homeostasis.

Strain-specific transcriptional responses overshadow salinity effects in a marine diatom sampled along the Baltic Sea salinity cline

The salinity gradient separating marine and freshwater environments represents a major ecological divide for microbiota, yet the mechanisms by which marine microbes have adapted to and ultimately diversified in freshwater environments are poorly understood. Here, we take advantage of a natural evolutionary experiment: the colonization of the brackish Baltic Sea by the ancestrally marine diatom Skeletonema marinoi. To understand how diatoms respond to low salinity, we characterized transcriptomic responses of acclimated S. marinoi grown in a common garden. Our experiment included eight strains from source populations spanning the Baltic Sea salinity cline. Gene expression analysis revealed that low salinities induced changes in the cellular metabolism of S. marinoi, including upregulation of photosynthesis and storage compound biosynthesis, increased nutrient demand, and a complex response to oxidative stress. However, the strain effect overshadowed the salinity effect, as strains differed significantly in their response, both regarding the strength and the strategy (direction of gene expression) of their response. The high degree of intraspecific variation in gene expression observed here highlights an important but often overlooked source of biological variation associated with how diatoms respond to environmental change.

Molecular Insights into Lipoxygenases in Diatoms Based on Structure Prediction: a Pioneering Study on Lipoxygenases Found in Pseudo-nitzschia arenysensis and Fragilariopsis cylindrus

Diatoms produce a variety of oxylipins which are oxygenated polyunsaturated fatty acids and are involved in chemical defense and intercellular communication, among other roles. Although the chemistry of diatom oxylipins has long been studied, the enzymes involved in their production, in particular lipoxygenase (LOX), which catalyzes the initial reaction of the synthesis, have not been discovered in diatom genomes. Recently, diatom LOXs were found in two species, Pseudo-nitzschia arenysensis (PaLOX) and Fragilariopsis cylindrus (FcLOX); however, the enzymology of these LOXs is largely unknown. In this review article, we discuss the potential functions of the diatom LOXs based on previously reported structures of LOXs derived from various organisms other than diatoms. Since the structures of PaLOX and FcLOX have not yet been solved, we discussed their functions, such as regio- and stereospecificities, on the basis of their structures predicted using a computational tool based on deep learning technology. Both diatom LOXs were predicted to conserve common core domains with relatively wide substrate-binding pockets. The stereo-determinant residues in PaLOX and FcLOX suggest S specificity. We assume that the highly conserved common core domain can be a clue to reveal unidentified lox genes from the accumulated diatom genome information with the aid of high-throughput structure prediction tools and structure-based alignment tools in the near future.

Rhodopsin-bestrophin fusion proteins from unicellular algae form gigantic pentameric ion channels

Many organisms sense light using rhodopsins, photoreceptive proteins containing a retinal chromophore. Here we report the discovery, structure and biophysical characterization of bestrhodopsins, a microbial rhodopsin subfamily from marine unicellular algae, in which one rhodopsin domain of eight transmembrane helices or, more often, two such domains in tandem, are C-terminally fused to a bestrophin channel. Cryo-EM analysis of a rhodopsin-rhodopsin-bestrophin fusion revealed that it forms a pentameric megacomplex (~700 kDa) with five rhodopsin pseudodimers surrounding the channel in the center. Bestrhodopsins are metastable and undergo photoconversion between red- and green-absorbing or green- and UVA-absorbing forms in the different variants. The retinal chromophore, in a unique binding pocket, photoisomerizes from all-trans to 11-cis form. Heterologously expressed bestrhodopsin behaves as a light-modulated anion channel.

Shedding light on silica biomineralization by comparative analysis of the silica-associated proteomes from three diatom species

Morphogenesis of the intricate patterns of diatom silica cell walls is a protein-guided process, yet to date only very few such silica biomineralization proteins have been identified. Therefore, it is currently unknown whether all diatoms share conserved proteins of a basal silica forming machinery, and whether unique proteins are responsible for the morphogenesis of species-specific silica patterns. To answer these questions, we extracted proteins from the silica of three diatom species (Thalassiosira pseudonana, Thalassiosira oceanica, and Cyclotella cryptica) by complete demineralization of the cell walls. Liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS) analysis of the extracts identified 92 proteins that we name 'soluble silicome proteins' (SSPs). Surprisingly, no SSPs are common to all three species, and most SSPs showed very low similarity to one another in sequence alignments. In-depth bioinformatics analyses revealed that SSPs could be grouped into distinct classes based on short unconventional sequence motifs whose functions are yet unknown. The results from the in vivo localization of selected SSPs indicates that proteins, which lack sequence homology but share unconventional sequence motifs may exert similar functions in the morphogenesis of the diatom silica cell wall.

Genome evolution of a nonparasitic secondary heterotroph, the diatom Nitzschia putrida

Secondary loss of photosynthesis is observed across almost all plastid-bearing branches of the eukaryotic tree of life. However, genome-based insights into the transition from a phototroph into a secondary heterotroph have so far only been revealed for parasitic species. Free-living organisms can yield unique insights into the evolutionary consequence of the loss of photosynthesis, as the parasitic lifestyle requires specific adaptations to host environments. Here, we report on the diploid genome of the free-living diatom Nitzschia putrida (35 Mbp), a nonphotosynthetic osmotroph whose photosynthetic relatives contribute ca. 40% of net oceanic primary production. Comparative analyses with photosynthetic diatoms and heterotrophic algae with parasitic lifestyle revealed that a combination of gene loss, the accumulation of genes involved in organic carbon degradation, a unique secretome, and the rapid divergence of conserved gene families involved in cell wall and extracellular metabolism appear to have facilitated the lifestyle of a free-living secondary heterotroph.

Impaired photoprotection in Phaeodactylum tricornutum KEA3 mutants reveals the proton regulatory circuit of diatoms light acclimation

Diatoms are successful phytoplankton clades able to acclimate to changing environmental conditions, including e.g. variable light intensity. Diatoms are outstanding at dissipating light energy exceeding the maximum photosynthetic electron transfer (PET) capacity via the nonphotochemical quenching (NPQ) process. While the molecular effectors of NPQ as well as the involvement of the proton motive force (PMF) in its regulation are known, the regulators of the PET/PMF relationship remain unidentified in diatoms. We generated mutants of the H⁺ /K⁺ antiporter KEA3 in the model diatom Phaeodactylum tricornutum. Loss of KEA3 activity affects the PET/PMF coupling and NPQ responses at the onset of illumination, during transients and in steady-state conditions. Thus, this antiporter is a main regulator of the PET/PMF coupling. Consistent with this conclusion, a parsimonious model including only two free components, KEA3 and the diadinoxanthin de-epoxidase, describes most of the feedback loops between PET and NPQ. This simple regulatory system allows for efficient responses to fast (minutes) or slow (e.g. diel) changes in light environment, thanks to the presence of a regulatory calcium ion (Ca²⁺ )-binding domain in KEA3 modulating its activity. This circuit is likely tuned by the NPQ-effector proteins, LHCXs, providing diatoms with the required flexibility to thrive in different ocean provinces.

An enigmatic stramenopile sheds light on early evolution in Ochrophyta plastid organellogenesis

Ochrophyta is an algal group belonging to the Stramenopiles and comprises diverse lineages of algae which contribute significantly to the oceanic ecosystems as primary producers. However, early evolution of the plastid organelle in Ochrophyta is not fully understood. In this study, we provide a well-supported tree of the Stramenopiles inferred by the large-scale phylogenomic analysis that unveils the eukaryvorous (nonphotosynthetic) protist Actinophrys sol (Actinophryidae) is closely related to Ochrophyta. We used genomic and transcriptomic data generated from A. sol to detect molecular traits of its plastid and we found no evidence of plastid genome and plastid-mediated biosynthesis, consistent with previous ultrastructural studies that did not identify any plastids in Actinophryidae. Moreover, our phylogenetic analyses of particular biosynthetic pathways provide no evidence of a current and past plastid in A. sol. However, we found more than a dozen organellar aminoacyl-tRNA synthases (aaRSs) that are of algal origin. Close relationships between aaRS from A. sol and their ochrophyte homologs document gene transfer of algal genes that happened before the divergence of Actinophryidae and Ochrophyta lineages. We further showed experimentally that organellar aaRSs of A. sol are targeted exclusively to mitochondria, although organellar aaRSs in Ochrophyta are dually targeted to mitochondria and plastids. Together, our findings suggested that the last common ancestor of Actinophryidae and Ochrophyta had not yet completed the establishment of host–plastid partnership as seen in the current Ochrophyta species, but acquired at least certain nuclear-encoded genes for the plastid functions.

Silencing of a Pseudo-nitzschia arenysensis lipoxygenase transcript leads to reduced oxylipin production and impaired growth

Because of their importance as chemical mediators, the presence of a rich and varied family of lipoxygenase (LOX) products, collectively named oxylipins, has been investigated thoroughly in diatoms, and the involvement of these products in important processes such as bloom regulation has been postulated. Nevertheless, little information is available on the enzymes and pathways operating in these protists. Exploiting transcriptome data, we identified and characterized a LOX gene, PaLOX, in Pseudo-nitzschia arenysensis, a marine diatom known to produce different species of oxylipins by stereo- and regio-selective oxidation of eicosapentaenoic acid (EPA) at C12 and C15. PaLOX RNA interference correlated with a decrease of the lipid-peroxidizing activity and oxylipin synthesis, as well as with a reduction of growth of P. arenysensis. In addition, sequence analysis and structure models of the C-terminal part of the predicted protein closely fitted with the data for established LOXs from other organisms. The presence in the genome of a single LOX gene, whose downregulation impairs both 12- and 15-oxylipins synthesis, together with the in silico 3D protein modelling suggest that PaLOX encodes for a 12/15S-LOX with a dual specificity, and provides additional support to the correlation between cell growth and oxylipin biosynthesis in diatoms.

Proteomic analysis of metabolic pathways supports chloroplast-mitochondria cross-talk in a Cu-limited diatom

Diatoms are one of the most successful phytoplankton groups in our oceans, being responsible for over 20% of the Earth's photosynthetic productivity. Their chimeric genomes have genes derived from red algae, green algae, bacteria, and heterotrophs, resulting in multiple isoenzymes targeted to different cellular compartments with the potential for differential regulation under nutrient limitation. The resulting interactions between metabolic pathways are not yet fully understood. We previously showed how acclimation to Cu limitation enhanced susceptibility to overreduction of the photosynthetic electron transport chain and its reorganization to favor photoprotection over light harvesting in the oceanic diatom Thalassiosira oceanica (Hippmann et al., 2017, 10.1371/journal.pone.0181753). In order to gain a better understanding of the overall metabolic changes that help alleviate the stress of Cu limitation, we have further analyzed the comprehensive proteomic datasets generated in that study to identify differentially expressed proteins involved in carbon, nitrogen, and oxidative stress-related metabolic pathways. Metabolic pathway analysis showed integrated responses to Cu limitation. The upregulation of ferredoxin (Fdx) was correlated with upregulation of plastidial Fdx-dependent isoenzymes involved in nitrogen assimilation as well as enzymes involved in glutathione synthesis, thus suggesting an integration of nitrogen uptake and metabolism with photosynthesis and oxidative stress resistance. The differential expression of glycolytic isoenzymes located in the chloroplast and mitochondria may enable them to channel both excess electrons and/or ATP between these compartments. An additional support for chloroplast-mitochondrial cross-talk is the increased expression of chloroplast and mitochondrial proteins involved in the proposed malate shunt under Cu limitation.

Genome-Scale Metabolic Modeling Enables In-Depth Understanding of Big Data

Genome-scale metabolic models (GEMs) enable the mathematical simulation of the metabolism of archaea, bacteria, and eukaryotic organisms. GEMs quantitatively define a relationship between genotype and phenotype by contextualizing different types of Big Data (e.g., genomics, metabolomics, and transcriptomics). In this review, we analyze the available Big Data useful for metabolic modeling and compile the available GEM reconstruction tools that integrate Big Data. We also discuss recent applications in industry and research that include predicting phenotypes, elucidating metabolic pathways, producing industry-relevant chemicals, identifying drug targets, and generating knowledge to better understand host-associated diseases. In addition to the up-to-date review of GEMs currently available, we assessed a plethora of tools for developing new GEMs that include macromolecular expression and dynamic resolution. Finally, we provide a perspective in emerging areas, such as annotation, data managing, and machine learning, in which GEMs will play a key role in the further utilization of Big Data.

Identification of sequence motifs in Lhcx proteins that confer qE-based photoprotection in the diatom Phaeodactylum tricornutum

Photosynthetic organisms in nature often experience light fluctuations. While low light conditions limit the energy uptake by algae, light absorption exceeding the maximal rate of photosynthesis may go along with enhanced formation of potentially toxic reactive oxygen species. To preempt high light-induced photodamage, photosynthetic organisms evolved numerous photoprotective mechanisms. Among these, energy-dependent fluorescence quenching (qE) provides a rapid mechanism to dissipate thermally the excessively absorbed energy. Diatoms thrive in all aquatic environments and thus belong to the most important primary producers on earth. qE in diatoms is provided by a concerted action of Lhcx proteins and the xanthophyll cycle pigment diatoxanthin. While the exact Lhcx activation mechanism of diatom qE is unknown, two lumen-exposed acidic amino acids within Lhcx proteins were proposed to function as regulatory switches upon light-induced lumenal acidification. By introducing a modified Lhcx1 lacking these amino acids into a Phaeodactylum tricornutum Lhcx1-null qE knockout line, we demonstrate that qE is unaffected by these two amino acids. Based on sequence comparisons with Lhcx4, being incapable of providing qE, we perform domain swap experiments of Lhcx4 with Lhcx1 and identify two peptide motifs involved in conferring qE. Within one of these motifs, we identify a tryptophan residue with a major influence on qE establishment. This tryptophan residue is located in close proximity to the diadinoxanthin/diatoxanthin-binding site based on the recently revealed diatom Lhc crystal structure. Our findings provide a structural explanation for the intimate link of Lhcx and diatoxanthin in providing qE in diatoms.

The genome of a nonphotosynthetic diatom provides insights into the metabolic shift to heterotrophy and constraints on the loss of photosynthesis

Although most of the tens of thousands of diatom species are photoautotrophs, a small number of heterotrophic species no longer photosynthesize. We sequenced the genome of a nonphotosynthetic diatom, Nitzschia Nitz4, to determine how carbon metabolism was altered in the wake of this trophic shift. Nitzschia Nitz4 has retained its plastid and plastid genome, but changes associated with the transition to heterotrophy were cellular-wide and included losses of photosynthesis-related genes from the nuclear and plastid genomes, elimination of isoprenoid biosynthesis in the plastid, and remodeling of mitochondrial glycolysis to maximize adenosine triphosphte (ATP) yield. The genome contains a β-ketoadipate pathway that may allow Nitzschia Nitz4 to metabolize lignin-derived compounds. Diatom plastids lack an oxidative pentose phosphate pathway (oPPP), leaving photosynthesis as the primary source of NADPH to support essential biosynthetic pathways in the plastid and, by extension, limiting available sources of NADPH in nonphotosynthetic plastids. The genome revealed similarities between nonphotosynthetic diatoms and apicomplexan parasites for provisioning NADPH in their plastids and highlighted the ancestral absence of a plastid oPPP as a potentially important constraint on loss of photosynthesis, a hypothesis supported by the higher frequency of transitions to parasitism or heterotrophy in lineages that have a plastid oPPP.

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

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

Fatty acid photodecarboxylase is an ancient photoenzyme that forms hydrocarbons in the thylakoids of algae

Fatty acid photodecarboxylase (FAP) is one of the few enzymes that require light for their catalytic cycle (photoenzymes). FAP was first identified in the microalga Chlorella variabilis NC64A, and belongs to an algae-specific subgroup of the glucose-methanol-choline oxidoreductase family. While the FAP from C. variabilis and its Chlamydomonas reinhardtii homolog CrFAP have demonstrated in vitro activities, their activities and physiological functions have not been studied in vivo. Furthermore, the conservation of FAP activity beyond green microalgae remains hypothetical. Here, using a C. reinhardtii FAP knockout line (fap), we showed that CrFAP is responsible for the formation of 7-heptadecene, the only hydrocarbon of this alga. We further showed that CrFAP was predominantly membrane-associated and that >90% of 7-heptadecene was recovered in the thylakoid fraction. In the fap mutant, photosynthetic activity was not affected under standard growth conditions, but was reduced after cold acclimation when light intensity varied. A phylogenetic analysis that included sequences from Tara Ocean identified almost 200 putative FAPs and indicated that FAP was acquired early after primary endosymbiosis. Within Bikonta, FAP was retained in secondary photosynthetic endosymbiosis lineages but absent from those that lost the plastid. Characterization of recombinant FAPs from various algal genera (Nannochloropsis, Ectocarpus, Galdieria, Chondrus) provided experimental evidence that FAP photochemical activity was present in red and brown algae, and was not limited to unicellular species. These results thus indicate that FAP was conserved during the evolution of most algal lineages where photosynthesis was retained, and suggest that its function is linked to photosynthetic membranes.

Diurnal transcript profiling of the diatom Seminavis robusta reveals adaptations to a benthic lifestyle

Coastal regions contribute an estimated 20% of annual gross primary production in the oceans, despite occupying only 0.03% of their surface area. Diatoms frequently dominate coastal sediments, where they experience large variations in light regime resulting from the interplay of diurnal and tidal cycles. Here, we report on an extensive diurnal transcript profiling experiment of the motile benthic diatom Seminavis robusta. Nearly 90% (23 328) of expressed protein-coding genes and 66.9% (1124) of expressed long intergenic non-coding RNAs showed significant expression oscillations and are predominantly phasing at night with a periodicity of 24 h. Phylostratigraphic analysis found that rhythmic genes are enriched in highly conserved genes, while diatom-specific genes are predominantly associated with midnight expression. Integration of genetic and physiological cell cycle markers with silica depletion data revealed potential new silica cell wall-associated gene families specific to diatoms. Additionally, we observed 1752 genes with a remarkable semidiurnal (12-h) periodicity, while the expansion of putative circadian transcription factors may reflect adaptations to cope with highly unpredictable external conditions. Taken together, our results provide new insights into the adaptations of diatoms to the benthic environment and serve as a valuable resource for the study of diurnal regulation in photosynthetic eukaryotes.

Using Diatom and Apicomplexan Models to Study the Heme Pathway of Chromera velia

Heme biosynthesis is essential for almost all living organisms. Despite its conserved function, the pathway’s enzymes can be located in a remarkable diversity of cellular compartments in different organisms. This location does not always reflect their evolutionary origins, as might be expected from the history of their acquisition through endosymbiosis. Instead, the final subcellular localization of the enzyme reflects multiple factors, including evolutionary origin, demand for the product, availability of the substrate, and mechanism of pathway regulation. The biosynthesis of heme in the apicomonad Chromera velia follows a chimeric pathway combining heme elements from the ancient algal symbiont and the host. Computational analyses using different algorithms predict complex targeting patterns, placing enzymes in the mitochondrion, plastid, endoplasmic reticulum, or the cytoplasm. We employed heterologous reporter gene expression in the apicomplexan parasite Toxoplasma gondii and the diatom Phaeodactylum tricornutum to experimentally test these predictions. 5-aminolevulinate synthase was located in the mitochondria in both transfection systems. In T. gondii, the two 5-aminolevulinate dehydratases were located in the cytosol, uroporphyrinogen synthase in the mitochondrion, and the two ferrochelatases in the plastid. In P. tricornutum, all remaining enzymes, from ALA-dehydratase to ferrochelatase, were placed either in the endoplasmic reticulum or in the periplastidial space.

Phylogenomic fingerprinting of tempo and functions of horizontal gene transfer within ochrophytes

Horizontal gene transfer (HGT) is an important source of novelty in eukaryotic genomes. This is particularly true for the ochrophytes, a diverse and important group of algae. Previous studies have shown that ochrophytes possess a mosaic of genes derived from bacteria and eukaryotic algae, acquired through chloroplast endosymbiosis and from HGTs, although understanding of the time points and mechanisms underpinning these transfers has been restricted by the depth of taxonomic sampling possible. We harness an expanded set of ochrophyte sequence libraries, alongside automated and manual phylogenetic annotation, in silico modeling, and experimental techniques, to assess the frequency and functions of HGT across this lineage. Through manual annotation of thousands of single-gene trees, we identify continuous bacterial HGT as the predominant source of recently arrived genes in the model diatom Phaeodactylum tricornutum Using a large-scale automated dataset, a multigene ochrophyte reference tree, and mathematical reconciliation of gene trees, we note a probable elevation of bacterial HGTs at foundational points in diatom evolution, following their divergence from other ochrophytes. Finally, we demonstrate that throughout ochrophyte evolutionary history, bacterial HGTs have been enriched in genes encoding secreted proteins. Our study provides insights into the sources and frequency of HGTs, and functional contributions that HGT has made to algal evolution.

Plastidial acyl carrier protein Δ9-desaturase modulates eicosapentaenoic acid biosynthesis and triacylglycerol accumulation in Phaeodactylum tricornutum

The unicellular marine diatom Phaeodactylum tricornutum accumulates up to 35% eicosapentaenoic acid (EPA, 20:5n3) and has been used as a model organism to study long chain polyunsaturated fatty acids (LC-PUFA) biosynthesis due to an excellent annotated genome sequence and established transformation system. In P. tricornutum, the majority of EPA accumulates in polar lipids, particularly in galactolipids such as mono- and di-galactosyldiacylglycerol. LC-PUFA biosynthesis is considered to start from oleic acid (18:1n9). EPA can be synthesized via a series of desaturation and elongation steps occurring at the endoplasmic reticulum and newly synthesized EPA is then imported into the plastids for incorporation into galactolipids via an unknown route. The basis for the flux of EPA is fundamental to understanding LC-PUFA biosynthesis in diatoms. We used P. tricornutum to study acyl modifying activities, upstream of 18:1n9, on subsequent LC-PUFA biosynthesis. We identified the gene coding for the plastidial acyl carrier protein Δ9-desaturase, a key enzyme in fatty acid modification and analyzed the impact of overexpression and knock out of this gene on glycerolipid metabolism. This revealed a previously unknown role of this soluble desaturase in EPA synthesis and production of triacylglycerol. This study provides further insight into the distinctive nature of lipid metabolism in the marine diatom P. tricornutum and suggests additional approaches for tailoring oil composition in microalgae.

Transcriptomes of an oceanic diatom reveal the initial and final stages of acclimation to copper deficiency

Copper (Cu) concentration is greatly reduced in the open sea so that phytoplankton must adjust their uptake systems and acclimate to sustain growth. Acclimation to low Cu involves changes to the photosynthetic apparatus and specific biochemical reactions that use Cu, but little is known how Cu affects cellular metabolic networks. Here we report results of whole transcriptome analysis of a plastocyanin-containing diatom, Thalassiosira oceanica 1005, during its initial stages of acclimation and after long-term adaptation in Cu-deficient seawater. Gene expression profiles, used to identify Cu-regulated metabolic pathways, show downregulation of anabolic and energy-yielding reactions in Cu-limited cells. These include the light reactions of photosynthesis, carbon fixation, nitrogen assimilation and glycolysis. Reduction of these pathways is consistent with reduced growth requirements for C and N caused by slower rates of photosynthetic electron transport. Upregulation of oxidative stress defence systems persists in adapted cells, suggesting cellular damage by increased reactive oxygen species (ROS) occurs even after acclimation. Copper deficiency also alters fatty acid metabolism, possibly in response to an increase in lipid peroxidation and membrane damage driven by ROS. During the initial stages of Cu-limitation the majority of differentially regulated genes are associated with photosynthetic metabolism, highlighting the chloroplast as the primary target of low Cu availability. The results provide insights into the mechanisms of acclimation and adaptation of T. oceanica to Cu deficiency.

ppGpp influences protein protection, growth and photosynthesis in Phaeodactylum tricornutum

Chloroplasts retain elements of a bacterial stress response pathway that is mediated by the signalling nucleotides guanosine penta- and tetraphosphate ((p)ppGpp). In the model flowering plant Arabidopsis, ppGpp acts as a potent regulator of plastid gene expression and influences photosynthesis, plant growth and development. However, little is known about ppGpp metabolism or its evolution in other photosynthetic eukaryotes. Here, we studied the function of ppGpp in the diatom Phaeodactylum tricornutum using transgenic lines containing an inducible system for ppGpp accumulation. We used these lines to investigate the effects of ppGpp on growth, photosynthesis, lipid metabolism and protein expression. We demonstrate that ppGpp accumulation reduces photosynthetic capacity and promotes a quiescent-like state with reduced proliferation and ageing. Strikingly, using nontargeted proteomics, we discovered that ppGpp accumulation also leads to the coordinated upregulation of a protein protection response in multiple cellular compartments. Our findings highlight the importance of ppGpp as a fundamental regulator of chloroplast function across different domains of life, and lead to new questions about the molecular mechanisms and roles of (p)ppGpp signalling in photosynthetic eukaryotes.

Comprehensive and Functional Analysis of Horizontal Gene Transfer Events in Diatoms

Diatoms are a diverse group of mainly photosynthetic algae, responsible for 20% of worldwide oxygen production, which can rapidly respond to favorable conditions and often outcompete other phytoplankton. We investigated the contribution of horizontal gene transfer (HGT) to its ecological success. A large-scale phylogeny-based prokaryotic HGT detection procedure across nine sequenced diatoms showed that 3-5% of their proteome has a horizontal origin and a large influx occurred at the ancestor of diatoms. More than 90% of HGT genes are expressed, and species-specific HGT genes in Phaeodactylum tricornutum undergo strong purifying selection. Genes derived from HGT are implicated in several processes including environmental sensing and expand the metabolic toolbox. Cobalamin (vitamin B12) is an essential cofactor for roughly half of the diatoms and is only produced by bacteria. Five consecutive genes involved in the final synthesis of the cobalamin biosynthetic pathway, which could function as scavenging and repair genes, were detected as HGT. The full suite of these genes was detected in the cold-adapted diatom Fragilariopsis cylindrus. This might give diatoms originating from the Southern Ocean, a region typically depleted in cobalamin, a competitive advantage. Overall, we show that HGT is a prevalent mechanism that is actively used in diatoms to expand its adaptive capabilities.

Integrative analysis of chloroplast DNA methylation in a marine alga-Saccharina japonica

KEY MESSAGE

We applied an integrative approach using multiple methods to verify cytosine methylation in the chloroplast DNA of the multicellular brown alga Saccharina japonica. Cytosine DNA methylation is a heritable process which plays important roles in regulating development throughout the life cycle of an organism. Although methylation of nuclear DNA has been studied extensively, little is known about the state and role of DNA methylation in chloroplast genomes, especially in marine algae. Here, we have applied an integrated approach encompassing whole-genome bisulfite sequencing, methylated DNA immunoprecipitation, gene co-expression networks and photophysiological analyses to provide evidence for the role of chloroplast DNA methylation in a marine alga, the multicellular brown alga Saccharina japonica. Although the overall methylation level was relatively low in the chloroplast genome of S. japonica, gametophytes exhibited higher methylation levels than sporophytes. Gene-specific bisulfite-cloning sequencing provided additional evidence for the methylation of key photosynthetic genes. Many of them were highly expressed in sporophytes whereas genes involved in transcription, translation and biosynthesis were strongly expressed in gametophytes. Nucleus-encoded photosynthesis genes were co-expressed with their chloroplast-encoded counterparts potentially contributing to the higher photosynthetic performance in sporophytes compared to gametophytes where these co-expression networks were less pronounced. A nucleus-encoded DNA methyltransferase of the DNMT2 family is assumed to be responsible for the methylation of the chloroplast genome because it is predicted to possess a plastid transit peptide.

Genome-scale metabolic model of the diatom Thalassiosira pseudonana highlights the importance of nitrogen and sulfur metabolism in redox balance

Diatoms are unicellular photosynthetic algae known to secrete organic matter that fuels secondary production in the ocean, though our knowledge of how their physiology impacts the composition of dissolved organic matter remains limited. Like all photosynthetic organisms, their use of light for energy and reducing power creates the challenge of avoiding cellular damage. To better understand the interplay between redox balance and organic matter secretion, we reconstructed a genome-scale metabolic model of Thalassiosira pseudonana strain CCMP 1335, a model for diatom molecular biology and physiology, with a 60-year history of studies. The model simulates the metabolic activities of 1,432 genes via a network of 2,792 metabolites produced through 6,079 reactions distributed across six subcellular compartments. Growth was simulated under different steady-state light conditions (5–200 μmol photons m^-2 s^-1) and in a batch culture progressing from exponential growth to nitrate-limitation and nitrogen-starvation. We used the model to examine the dissipation of reductants generated through light-dependent processes and found that when available, nitrate assimilation is an important means of dissipating reductants in the plastid; under nitrate-limiting conditions, sulfate assimilation plays a similar role. The use of either nitrate or sulfate uptake to balance redox reactions leads to the secretion of distinct organic nitrogen and sulfur compounds. Such compounds can be accessed by bacteria in the surface ocean. The model of the diatom Thalassiosira pseudonana provides a mechanistic explanation for the production of ecologically and climatologically relevant compounds that may serve as the basis for intricate, cross-kingdom microbial networks. Diatom metabolism has an important influence on global biogeochemistry; metabolic models of marine microorganisms link genes to ecosystems and may be key to integrating molecular data with models of ocean biogeochemistry.

Iron metabolism strategies in diatoms

Diatoms are one of the most successful group of photosynthetic eukaryotes in the contemporary ocean. They are ubiquitously distributed and are the most abundant primary producers in polar waters. Equally remarkable is their ability to tolerate iron deprivation and respond to periodic iron fertilization. Despite their relatively large cell sizes, diatoms tolerate iron limitation and frequently dominate iron-stimulated phytoplankton blooms, both natural and artificial. Here, we review the main iron use strategies of diatoms, including their ability to assimilate and store a range of iron sources, and the adaptations of their photosynthetic machinery and architecture to iron deprivation. Our synthesis relies on published literature and is complemented by a search of 82 diatom transcriptomes, including information collected from seven representatives of the most abundant diatom genera in the world's oceans.

Proximity proteomics in a marine diatom reveals a putative cell surface-to-chloroplast iron trafficking pathway

Iron is a biochemically critical metal cofactor in enzymes involved in photosynthesis, cellular respiration, nitrate assimilation, nitrogen fixation, and reactive oxygen species defense. Marine microeukaryotes have evolved a phytotransferrin-based iron uptake system to cope with iron scarcity, a major factor limiting primary productivity in the global ocean. Diatom phytotransferrin is endocytosed; however, proteins downstream of this environmentally ubiquitous iron receptor are unknown. We applied engineered ascorbate peroxidase APEX2-based subcellular proteomics to catalog proximal proteins of phytotransferrin in the model marine diatom Phaeodactylum tricornutum. Proteins encoded by poorly characterized iron-sensitive genes were identified including three that are expressed from a chromosomal gene cluster. Two of them showed unambiguous colocalization with phytotransferrin adjacent to the chloroplast. Further phylogenetic, domain, and biochemical analyses suggest their involvement in intracellular iron processing. Proximity proteomics holds enormous potential to glean new insights into iron acquisition pathways and beyond in these evolutionarily, ecologically, and biotechnologically important microalgae.

Genome-wide identification of chitinase genes in Thalassiosira pseudonana and analysis of their expression under abiotic stresses

BACKGROUND

The nitrogen-containing polysaccharide chitin is the second most abundant biopolymer on earth and is found in the cell walls of diatoms, where it serves as a scaffold for biosilica deposition. Diatom chitin is an important source of carbon and nitrogen in the marine environment, but surprisingly little is known about basic chitinase metabolism in diatoms.

RESULTS

Here, we identify and fully characterize 24 chitinase genes from the model centric diatom Thalassiosira pseudonana. We demonstrate that their expression is broadly upregulated under abiotic stresses, despite the fact that chitinase activity itself remains unchanged, and we discuss several explanations for this result. We also examine the potential transcriptional complexity of the intron-rich T. pseudonana chitinase genes and provide evidence for two separate tandem duplication events during their evolution.

CONCLUSIONS

Given the many applications of chitin and chitin derivatives in suture production, wound healing, drug delivery, and other processes, new insight into diatom chitin metabolism has both theoretical and practical value.

PhycoCosm, a comparative algal genomics resource

Algae are a diverse, polyphyletic group of photosynthetic eukaryotes spanning nearly all eukaryotic lineages of life and collectively responsible for ∼50% of photosynthesis on Earth. Sequenced algal genomes, critical to understanding their complex biology, are growing in number and require efficient tools for analysis. PhycoCosm (https://phycocosm.jgi.doe.gov) is an algal multi-omics portal, developed by the US Department of Energy Joint Genome Institute to support analysis and distribution of algal genome sequences and other ‘omics’ data. PhycoCosm provides integration of genome sequence and annotation for >100 algal genomes with available multi-omics data and interactive web-based tools to enable algal research in bioenergy and the environment, encouraging community engagement and data exchange, and fostering new sequencing projects that will further these research goals.

Characterization of a Novel Acyl-ACP Δ9 Desaturase Gene Responsible for Palmitoleic Acid Accumulation in a Diatom Phaeodactylum tricornutum

Palmitoleic acid (16:1Δ⁹) possesses a double bond at the seventh carbon atom from methyl end of the acyl chain and belongs to unusual ω-7 monounsaturated fatty acids with broad applications in food, pharmaceuticals, cosmetics, biofuel, and other industries. This high-value fatty acid accumulates up to >40% of total lipid in the marine diatom Phaeodactylum tricornutum. The present study was conducted to determine the key gene responsible for 16:1Δ⁹ biosynthesis in this unicellular alga. A new full-length cDNA and genomic DNA encoding acyl-ACP Δ⁹ desaturase (PtAAD) were isolated from P. tricornutum cells. Expression levels of PtAAD gene under normal and stress culture conditions were both positively correlated with 16:1Δ⁹ accumulation, implying its potential role for fatty acid determination. Functional complementation assay of a yeast mutant strain BY4839 evidenced that PtAAD could restore the synthesis of unsaturated fatty acid, especially generating high levels of 16:1Δ⁹. Further transient expression of PtAAD gene in Nicotiana benthamiana leaves was accompanied by the accumulation of 16:1Δ⁹, which was absent from control groups. Three-dimensional structure modeling studies showed that functional domain of PtAAD contained three variant amino acids (F160, A223, and L156), which may narrow the space shape of substrate-binding cavity to ensure the entry of 16:0-ACP. Consistent with this prediction, the mutated version of PtAAD gene (F160L, A223T, and L156M) in N. benthamiana systems failed to accumulate 16:1Δ⁹, but increased levels of 18:1Δ⁹. Taken together, PtAAD exhibits a strong enzymatic activity and substrate preference for 16:0-ACP, acting as the key player for high biosynthesis and accumulation of 16:1Δ⁹ in this alga. These findings provide new insights for better understanding the palmitoleic acid and oil biosynthetic mechanism in P. tricornutum, indicating that PtAAD gene may have practical applications for enriching palmitoleic acid and oil yield in other commercial oleaginous algae and crops.

Transcriptional Response of Osmolyte Synthetic Pathways and Membrane Transporters in a Euryhaline Diatom During Long-term Acclimation to a Salinity Gradient

How diatoms respond to fluctuations in osmotic pressure is important from both ecological and applied perspectives. It is well known that osmotic stress affects photosynthesis and can result in the accumulation of compounds desirable in pharmaceutical and alternative fuel industries. Gene expression responses to osmotic stress have been studied in short-term trials, but it is unclear whether the same mechanisms are recruited during long-term acclimation. We used RNA-seq to study the genome-wide transcription patterns in the euryhaline diatom, Cyclotella cryptica, following long-term acclimation to salinity that spanned the natural range of fresh to oceanic water. Long-term acclimated C. cryptica exhibited induced synthesis or repressed degradation of the osmolytes glycine betaine, taurine and dimethylsulfoniopropionate (DMSP). Although changes in proline concentration is one of the main responses in short-term osmotic stress, we did not detect a transcriptional change in proline biosynthetic pathways in our long-term experiment. Expression of membrane transporters showed a general tendency for increased import of potassium and export of sodium, consistent with the electrochemical gradients and dependence on co-transported molecules. Our results show substantial between-genotype differences in growth and gene expression reaction norms and suggest that the regulation of proline synthesis important in short-term osmotic stress might not be maintained in long-term acclimation. Further examination using time-course gene expression experiments, metabolomics and genetic validation of gene functions would reinforce patterns inferred from RNA-seq data.

The genome of Ectocarpus subulatus - A highly stress-tolerant brown alga

Brown algae are multicellular photosynthetic stramenopiles that colonize marine rocky shores worldwide. Ectocarpus sp. Ec32 has been established as a genomic model for brown algae. Here we present the genome and metabolic network of the closely related species, Ectocarpus subulatus Kützing, which is characterized by high abiotic stress tolerance. Since their separation, both strains show new traces of viral sequences and the activity of large retrotransposons, which may also be related to the expansion of a family of chlorophyll-binding proteins. Further features suspected to contribute to stress tolerance include an expanded family of heat shock proteins, the reduction of genes involved in the production of halogenated defence compounds, and the presence of fewer cell wall polysaccharide-modifying enzymes. Overall, E. subulatus has mainly lost members of gene families down-regulated in low salinities, and conserved those that were up-regulated in the same condition. However, 96% of genes that differed between the two examined Ectocarpus species, as well as all genes under positive selection, were found to encode proteins of unknown function. This underlines the uniqueness of brown algal stress tolerance mechanisms as well as the significance of establishing E. subulatus as a comparative model for future functional studies.

Fatty Acid Biosynthesis in Chromerids

Fatty acids are essential components of biological membranes, important for the maintenance of cellular structures, especially in organisms with complex life cycles like protozoan parasites. Apicomplexans are obligate parasites responsible for various deadly diseases of humans and livestock. We analyzed the fatty acids produced by the closest phototrophic relatives of parasitic apicomplexans, the chromerids Chromera velia and Vitrella brassicaformis, and investigated the genes coding for enzymes involved in fatty acids biosynthesis in chromerids, in comparison to their parasitic relatives. Based on evidence from genomic and metabolomic data, we propose a model of fatty acid synthesis in chromerids: the plastid-localized FAS-II pathway is responsible for the de novo synthesis of fatty acids reaching the maximum length of 18 carbon units. Short saturated fatty acids (C14:0-C18:0) originate from the plastid are then elongated and desaturated in the cytosol and the endoplasmic reticulum. We identified giant FAS I-like multi-modular enzymes in both chromerids, which seem to be involved in polyketide synthesis and fatty acid elongation. This full-scale description of the biosynthesis of fatty acids and their derivatives provides important insights into the reductive evolutionary transition of a phototropic algal ancestor to obligate parasites.

Alanine to serine substitutions drive thermal adaptation in a psychrophilic diatom cytochrome c6

In this study, we investigate the thermodynamic mechanisms by which electron transfer proteins adapt to environmental temperature by directly comparing the redox properties and folding stability of a psychrophilic cytochrome c and a mesophilic homolog. Our model system consists of two cytochrome c₆ proteins from diatoms: one adapted specifically to polar environments, the other adapted generally to surface ocean environments. Direct electrochemistry shows that the midpoint potential for the mesophilic homolog is slightly higher at all temperatures measured. Cytochrome c₆ from the psychrophilic diatom unfolds with a melting temperature 10.4 °C lower than the homologous mesophilic cytochrome c₆. Changes in free energy upon unfolding are identical, within error, for the psychrophilic and mesophilic protein; however, the chemical unfolding transition of the psychrophilic cytochrome c₆ is more cooperative than for the mesophilic cytochrome c₆. Substituting alanine residues found in the mesophile with serine found in corresponding positions of the psychrophile demonstrates that burial of the polar serine both decreases the thermal stability and decreases the midpoint potential. The mutagenesis data, combined with differences in the m-value of chemical denaturation, suggest that differences in solvent accessibility of the hydrophobic core underlie the adaptation of cytochrome c₆ to differing environmental temperature.

Extrachromosomal Genetic Engineering of the Marine Diatom Phaeodactylum tricornutum Enables the Heterologous Production of Monoterpenoids

Geraniol is a commercially relevant plant-derived monoterpenoid that is a main component of rose essential oil and used as insect repellent. Geraniol is also a key intermediate compound in the biosynthesis of the monoterpenoid indole alkaloids (MIAs), a group of over 2000 compounds that include high-value pharmaceuticals. As plants naturally produce extremely small amounts of these molecules and their chemical synthesis is complex, industrially sourcing these compounds is costly and inefficient. Hence, microbial hosts suitable to produce MIA precursors through synthetic biology and metabolic engineering are currently being sought. Here, we evaluated the suitability of a eukaryotic microalga, the marine diatom Phaeodactylum tricornutum, for the heterologous production of monoterpenoids. Profiling of endogenous metabolism revealed that P. tricornutum, unlike other microbes employed for industrial production of terpenoids, accumulates free pools of the precursor geranyl diphosphate. To evaluate the potential for larger synthetic biology applications, we engineered P. tricornutum through extrachromosomal, episome-based expression, for the heterologous biosynthesis of the MIA intermediate geraniol. By profiling the production of geraniol resulting from various genetic and cultivation arrangements, P. tricornutum reached the maximum geraniol titer of 0.309 mg/L in phototrophic conditions. This work provides (i) a detailed analysis of P. tricornutum endogenous terpenoid metabolism, (ii) a successful demonstration of extrachromosomal expression for metabolic pathway engineering with potential gene-stacking applications, and (iii) a convincing proof-of-concept of the suitability of P. tricornutum as a novel production platform for heterologous monoterpenoids, with potential for complex pathway engineering aimed at the heterologous production of MIAs.