The Combination of RNA and Protein Profiling Reveals the Response to Nitrogen Depletion in Thalassiosira pseudonana

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.

Heterotrophic bacteria trigger transcriptome remodelling in the photosynthetic picoeukaryote Micromonas commoda

Marine biogeochemical cycles are built on interactions between surface ocean microbes, particularly those connecting phytoplankton primary producers to heterotrophic bacteria. Details of these associations are not well understood, especially in the case of direct influences of bacteria on phytoplankton physiology. Here we catalogue how the presence of three marine bacteria (Ruegeria pomeroyi DSS-3, Stenotrophomonas sp. SKA14 and Polaribacter dokdonensis MED152) individually and uniquely impact gene expression of the picoeukaryotic alga Micromonas commoda RCC 299. We find a dramatic transcriptomic remodelling by M. commoda after 8 h in co-culture, followed by an increase in cell numbers by 56 h compared with the axenic cultures. Some aspects of the algal transcriptomic response are conserved across all three bacterial co-cultures, including an unexpected reduction in relative expression of photosynthesis and carbon fixation pathways. Expression differences restricted to a single bacterium are also observed, with the Flavobacteriia P. dokdonensis uniquely eliciting changes in relative expression of algal genes involved in biotin biosynthesis and the acquisition and assimilation of nitrogen. This study reveals that M. commoda has rapid and extensive responses to heterotrophic bacteria in ways that are generalizable, as well as in a taxon specific manner, with implications for the diversity of phytoplankton-bacteria interactions ongoing in the surface ocean.

De novo transcriptome assembly and molecular response mechanism analysis of a diatom Cyclotella meneghiniana Kützing exposed to cadmium

Cadmium is a persistent heavy metal commonly found in aquatic ecosystems and has a strong toxic effect on organisms. The sensitivity of phytoplankton to environmental changes and its role as an indicator of aquatic ecosystem health have been well-established. However, the mechanisms by which phytoplankton respond to cadmium remain incompletely understood. In this study, we chose the typical planktonic diatom Cyclotella meneghiniana Kützing, by integrating physiological-biochemical data and transcriptome analysis, to reveal the molecular mechanisms of C. meneghiniana responing to cadmium. Under cadmium stress, the cell density and chlorophyll-a content of C. meneghiniana significantly decreased, while MDA content and SOD activity gradually increased. At 72 h of cadmium stress, we found that at this time point, cell abundance and physiological variation were very significant, therefore we selected 72 h for subsequent analysis. To better understand the cadmium stress response mechanisms of C. meneghiniana, a de novo transcriptome method was used to analyse C. meneghiniana under cadmium stress for 72 h, and 1704 (M vs. CK) and 4788 (H vs. CK) differentially expressed genes were found. Our results showed that the changes in gene expression were closely correlated to the physiological-biochemical changes. Although cadmium stress could promote the nitrogen metabolism pathway, ROS scavenging system, and photosynthesis. While, C. meneghiniana under medium and high concentrations of cadmium can also limit various intracellular metabolic pathways, such as the MAPK pathway and phosphatidylinositol metabolic pathway, and the degree of inhibition increases with the increase of stress concentration. In present study, the complete molecular mechanism of the planktonic diatom response to cadmium has been established, which provided important information for further studies on heavy metal pollutants and the multiple functional genes responsible for cadmium sensitivity and tolerance in planktonic diatoms.

Thalassiosira pseudonana (Cyclotella nana) (Hustedt) Hasle et Heimdal (Bacillariophyceae): A genetically tractable model organism for studying diatom biology, including biological silica formation

In 2004, Thalassiosira pseudonana was the first eukaryotic marine alga to have its genome sequenced. Since then, this species has quickly emerged as a valuable model species for investigating the molecular underpinnings of essentially all aspects of diatom life, particularly bio-morphogenesis of the cell wall. An important prerequisite for the model status of T. pseudonana is the ongoing development of increasingly precise tools to study the function of gene networks and their encoded proteins in vivo. Here, we briefly review the current toolbox for genetic manipulation, highlight specific examples of its application in studying diatom metabolism, and provide a peek into the role of diatoms in the emerging field of silica biotechnology.

Genome Sequences of Two Strains of Prototheca wickerhamii Provide Insight Into the Protothecosis Evolution

The Prototheca alga is the only chlorophyte known to be involved in a series of clinically relevant opportunistic infections in humans and animals, namely, protothecosis. Most pathogenic cases in humans are caused by Prototheca wickerhamii. In order to investigate the evolution of Prototheca and the genetic basis for its pathogenicity, the genomes of two P. wickerhamii strains S1 and S931 were sequenced using Nanopore long-read and Illumina short-read technologies. The mitochondrial, plastid, and nuclear genomes were assembled and annotated including a transcriptomic data set. The assembled nuclear genome size was 17.57 Mb with 19 contigs and 17.45 Mb with 26 contigs for strains S1 and S931, respectively. The number of predicted protein-coding genes was approximately 5,700, and more than 96% of the genes could be annotated with a gene function. A total of 2,798 gene families were shared between the five currently available Prototheca genomes. According to the phylogenetic analysis, the genus of Prototheca was classified in the same clade with A. protothecoides and diverged from Chlorella ~500 million years ago (Mya). A total of 134 expanded genes were enriched in several pathways, mostly in metabolic pathways, followed by biosynthesis of secondary metabolites and RNA transport. Comparative analysis demonstrated more than 96% consistency between the two herein sequenced strains. At present, due to the lack of sufficient understanding of the Prototheca biology and pathogenicity, the diagnosis rate of protothecosis is much lower than the actual infection rate. This study provides an in-depth insight into the genome sequences of two strains of P. wickerhamii isolated from the clinic to contribute to the basic understanding of this alga and explore future prevention and treatment strategies.

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.

Toxic effects and mechanisms of PFOA and its substitute GenX on the photosynthesis of Chlorella pyrenoidosa

Perfluorooctanoic acid (PFOA) and its substitute GenX are toxic chemicals that are widespread in the aquatic environment. However, there is little information about their toxicity mechanisms to aquatic organisms. In this study, Chlorella pyrenoidosa (C. pyrenoidosa) was treated with two concentrations (100 ng L^-1 and 100 μg L^-1) of PFOA or GenX for 12 days. The results showed that these two concentrations of PFOA and GenX began to inhibit the growth of algae after 6 days of treatment, and the Chlorophyll content and photosynthetic activity of C. pyrenoidosa were also negatively affected by these two chemicals. The transcriptomic results indicated that most of the genes related to the photosynthetic metabolism of C. pyrenoidosa were down-regulated (in 100 ng L^-1 treatment groups) on the 12th day. Besides, GenX and PFOA showed similar effects on algae photosynthesis including physical damage and metabolic disorders. According to this study, GenX might not be an ideal substitute for PFOA, and more attention should be paid on the management of emerging perfluoroalkyl substances.

Quantitative Proteomic Profiling of Marine Diatom Skeletonema dohrnii in Response to Temperature and Silicate Induced Environmental Stress

Global warming is expected to reduce the nutrient concentration in the upper ocean and affect the physiology of marine diatoms, but the underlying molecular mechanisms controlling these physiological changes are currently unknown. To understand these mechanisms, here we investigated iTRAQ based proteomic profiling of diatom Skeletonema dohrnii in a multifactorial experimental with a combining change of temperature and silicate concentrations. In total, 3369 differently abundant proteins were detected in four different environmental conditions, and the function of all proteins was identified using Gene Ontology and KEGG pathway analysis. For discriminating the proteome variation among samples, multivariate statistical analysis (PCA, PLS-DA) was performed by comparing the protein ratio differences. Further, performing pathway analysis on diatom proteomes, we here demonstrated downregulation of photosynthesis, carbon metabolism, and ribosome biogenesis in the cellular process that leads to decrease the oxidoreductase activity and affects the cell cycle of the diatom. Using PLS-DA VIP score plot analysis, we identified 15 protein biomarkers for discriminating studied samples. Of these, five proteins or gene (rbcL, PRK, atpB, DNA-binding, and signal transduction) identified as key biomarkers, induced by temperature and silicate stress in diatom metabolism. Our results show that proteomic finger-printing of S. dohrnii with different environmental conditions adds biological information that strengthens marine phytoplankton proteome analysis.

The molecular response mechanisms of a diatom Thalassiosira pseudonana to the toxicity of BDE-47 based on whole transcriptome analysis

Polybrominated diphenyl ethers (PBDEs) are ubiquitously distributed persistent organic pollutants (POPs) in marine environments. Phytoplankton are the entrance of PBDEs entering to biotic environments from abiotic environments, while the responding mechanisms of phytoplankton to PBDEs have not been full established. Therefore, we chose the model diatom Thalassiosira pseudonana in this study, by integrating whole transcriptome analysis with physiological-biochemical data, to reveal the molecular responding mechanisms of T. pseudonana to the toxicity of BDE-47. Our results indicated the changes of genes expressions correlated to the physiological-biochemical changes, and there were multiple molecular mechanisms of T. pseudonana responding to the toxicity of BDE-47: Gene expressions evidence explained the suppression of light reaction and proved the occurrence of cellular oxidative stress; In the meanwhile, up-regulations of genes in pathways involving carbon metabolisms happened, including the Calvin cycle, glycolysis, TCA cycle, fatty acid synthesis, and triacylglycerol synthesis; Lastly, DNA damage was found and three outcome including DNA repair, cell cycle arrest and programmed cell death (PCD) happened, which could finally inhibit the cell division and population growth of T. pseudonana. This study presented the most complete molecular responding mechanisms of phytoplankton cells to PBDEs, and provided valuable information of various PBDEs-sensitive genes with multiple functions for further research involving organic pollutants and phytoplankton.

Bioengineering of Microalgae: Recent Advances, Perspectives, and Regulatory Challenges for Industrial Application

Microalgae, due to their complex metabolic capacity, are being continuously explored for nutraceuticals, pharmaceuticals, and other industrially important bioactives. However, suboptimal yield and productivity of the bioactive of interest in local and robust wild-type strains are of perennial concerns for their industrial applications. To overcome such limitations, strain improvement through genetic engineering could play a decisive role. Though the advanced tools for genetic engineering have emerged at a greater pace, they still remain underused for microalgae as compared to other microorganisms. Pertaining to this, we reviewed the progress made so far in the development of molecular tools and techniques, and their deployment for microalgae strain improvement through genetic engineering. The recent availability of genome sequences and other omics datasets form diverse microalgae species have remarkable potential to guide strategic momentum in microalgae strain improvement program. This review focuses on the recent and significant improvements in the omics resources, mutant libraries, and high throughput screening methodologies helpful to augment research in the model and non-model microalgae. Authors have also summarized the case studies on genetically engineered microalgae and highlight the opportunities and challenges that are emerging from the current progress in the application of genome-editing to facilitate microalgal strain improvement. Toward the end, the regulatory and biosafety issues in the use of genetically engineered microalgae in commercial applications are described.

Comparative Proteomic Analysis Reveals New Insights Into the Common and Specific Metabolic Regulation of the Diatom Skeletonema dohrnii to the Silicate and Temperature Availability

Silicate (Si) and temperature are essential drivers for diatom growth and development in the ocean. Response of diatoms to these particular stress has been investigated; however, their common and specific responses to regulate intracellular development and growth are not known. Here, we investigated the combination of physiological characteristics and comparative proteomics of the diatom Skeletonema dohrnii grown in silicate- and temperature-limited conditions. Results show that cell carbon and lipid quotas were higher at lower-temperature cells, whereas cellular phosphate was higher in cells grown with lower Si. In silicate-limited cells, nitrate transporters were downregulated and resulted in lower nitrate assimilation, whereas the phosphate transporters and its assimilation were reduced in lower-temperature conditions. In photosynthesis, lower silicate caused impact in the linear electron flow and NADPH production, whereas cycling electron transport and ATP production were affected by the lower temperature. Concerning cell cycle, imbalances in the translation process were observed in lower-silicate cells, whereas impact in the transcription mechanism was observed in lower-temperature cells. However, proteins associated with carbon fixation and photorespiration were downregulated in both stress conditions, while the carbohydrate and lipid synthesis proteins were upregulated. Our results showed new insights into the common and specific responses on the proteome and physiology of S. dohrnii to silicate and temperature limitation, providing particular nutrient (Si)- and temperature-dependent mechanisms in diatoms.

Quantitative Proteomic Analysis Reveals Novel Insights into Intracellular Silicate Stress-Responsive Mechanisms in the Diatom Skeletonema dohrnii

Diatoms are a successful group of marine phytoplankton that often thrives under adverse environmental stress conditions. Members of the Skeletonema genus are ecologically important which may subsist during silicate stress and form a dense bloom following higher silicate concentration. However, our understanding of diatoms' underlying molecular mechanism involved in these intracellular silicate stress-responses are limited. Here an iTRAQ-based proteomic method was coupled with multiple physiological techniques to explore distinct cellular responses associated with oxidative stress in the diatom Skeletonema dohrnii to the silicate limitation. In total, 1768 proteins were detected; 594 proteins were identified as differentially expressed (greater than a two-fold change; p < 0.05). In Si-limited cells, downregulated proteins were mainly related to photosynthesis metabolism, light-harvesting complex, and oxidative phosphorylation, corresponding to inducing oxidative stress, and ROS accumulation. None of these responses were identified in Si-limited cells; in comparing with other literature, Si-stress cells showed that ATP-limited diatoms are unable to rely on photosynthesis, which will break down and reshuffle carbon metabolism to compensate for photosynthetic carbon fixation losses. Our findings have a good correlation with earlier reports and provides a new molecular level insight into the systematic intracellular responses employed by diatoms in response to silicate stress in the marine environment.

Physiological and transcriptomic analyses to determine the responses to phosphorus utilization in Nostoc sp

Phosphorus (P) is an important factor driving algal growth in aquatic ecosystems. In the present study, the growth, P uptake and utilization, photosynthesis, and transcriptome profile of Nostoc sp. were measured when Nostoc sp. cultured in media containing β-glycerol phosphate (β-gly, containing COP bonds), 2-aminoethylphosphonic acid (2-amin, containing CP bonds), or orthophosphate (K₂HPO₄), and in P-free (NP) medium. The results revealed that NP treatment adversely affected the growth and photosynthesis of Nostoc sp. and significantly down-regulated the expression of genes related to nutrient transport and material metabolism. Furthermore, 2-amin treatment reduced the growth of Nostoc sp. but did not significantly reduce photosynthesis, and the treatments of NP and 2-amin up-regulated the expressions of genes related antioxidation and stress. Additionally, there were no obvious differences in growth, photosynthesis, and phosphorus utilization between the β-gly and K₂HPO₄ treatments. These results suggested that Nostoc had a flexible ability to utilize P, which might play an important role in its widespread distribution in the environment.

Molecular adaptations to phosphorus deprivation and comparison with nitrogen deprivation responses in the diatom Phaeodactylum tricornutum

Phosphorus, an essential element for all living organisms, is a limiting nutrient in many regions of the ocean due to its fast recycling. Changes in phosphate (Pi) availability in aquatic systems affect diatom growth and productivity. We investigated the early adaptive mechanisms in the marine diatom Phaeodactylum tricornutum to P deprivation using a combination of transcriptomics, metabolomics, physiological and biochemical experiments. Our analysis revealed strong induction of gene expression for proteins involved in phosphate acquisition and scavenging, and down-regulation of processes such as photosynthesis, nitrogen assimilation and nucleic acid and ribosome biosynthesis. P deprivation resulted in alterations of carbon allocation through the induction of the pentose phosphate pathway and cytosolic gluconeogenesis, along with repression of the Calvin cycle. Reorganization of cellular lipids was indicated by coordinated induced expression of phospholipases, sulfolipid biosynthesis enzymes and a putative betaine lipid biosynthesis enzyme. A comparative analysis of nitrogen- and phosphorus-deprived P. tricornutum revealed both common and distinct regulation patterns in response to phosphate and nitrate stress. Regulation of central carbon metabolism and amino acid metabolism was similar, whereas unique responses were found in nitrogen assimilation and phosphorus scavenging in nitrogen-deprived and phosphorus-deprived cells, respectively.