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

Molecular Forecasting of Domoic Acid during a Pervasive Toxic Diatom Bloom

In 2015, the largest recorded harmful algal bloom (HAB) occurred in the Northeast Pacific, causing nearly 100 million dollars in damages to fisheries and killing many protected marine mammals. Dominated by the toxic diatom Pseudo-nitzschia australis , this bloom produced high levels of the neurotoxin domoic acid (DA). Through molecular and transcriptional characterization of 52 near-weekly phytoplankton net-tow samples collected at a bloom hotspot in Monterey Bay, California, we identified active transcription of known DA biosynthesis ( dab ) genes from the three identified toxigenic species, including P. australis as the primary origin of toxicity. Elevated expression of silicon transporters ( sit1 ) during the bloom supports the previously hypothesized role of dissolved silica (Si) exhaustion in contributing to bloom physiology and toxicity. We find that co-expression of the dabA and sit1 genes serves as a robust predictor of DA one week in advance, potentially enabling the forecasting of DA-producing HABs. We additionally present evidence that low levels of iron could have co-limited the diatom population along with low Si. Iron limitation represents a previously unrecognized driver of both toxin production and ecological success of the low iron adapted Pseudo-nitzschia genus during the 2015 bloom, and increasing pervasiveness of iron limitation may fuel the escalating magnitude and frequency of toxic Pseudo-nitzschia blooms globally. Our results advance understanding of bloom physiology underlying toxin production, bloom prediction, and the impact of global change on toxic blooms.

Significance

Pseudo-nitzschia diatoms form oceanic harmful algal blooms that threaten human health through production of the neurotoxin domoic acid (DA). DA biosynthetic gene expression is hypothesized to control DA production in the environment, yet what regulates expression of these genes is yet to be discovered. In this study, we uncovered expression of DA biosynthesis genes by multiple toxigenic Pseudo-nitzschia species during an economically impactful bloom along the North American West Coast, and identified genes that predict DA in advance of its production. We discovered that iron and silica co-limitation restrained the bloom and likely promoted toxin production. This work suggests that increasing iron limitation due to global change may play a previously unrecognized role in driving bloom frequency and toxicity.

Mechanistic constraints on the trade-off between photosynthesis and respiration in response to warming

Phytoplankton are responsible for half of all oxygen production and drive the ocean carbon cycle. Metabolic theory predicts that increasing global temperatures will cause phytoplankton to become more heterotrophic and smaller. Here, we uncover the metabolic trade-offs between cellular space, energy, and stress management driving phytoplankton thermal acclimation and how these might be overcome through evolutionary adaptation. We show that the observed relationships between traits such as chlorophyll, lipid content, C:N, and size can be predicted on the basis of the metabolic demands of the cell, the thermal dependency of transporters, and changes in membrane lipids. We suggest that many of the observed relationships are not fixed physiological constraints but rather can be altered through adaptation. For example, the evolution of lipid metabolism can favor larger cells with higher lipid content to mitigate oxidative stress. These results have implications for rates of carbon sequestration and export in a warmer ocean.

The cellular response to ocean warming in Emiliania huxleyi

Marine phytoplankton contribute substantially to the global flux of carbon from the atmosphere to the deep ocean. Sea surface temperatures will inevitably increase in line with global climate change, altering the performance of marine phytoplankton. Differing sensitivities of photosynthesis and respiration to temperature, will likely shift the strength of the future oceanic carbon sink. To further clarify the molecular mechanisms driving these alterations in phytoplankton function, shotgun proteomic analysis was carried out on the globally-occurring coccolithophore Emiliania huxleyi exposed to moderate- (23°C) and elevated- (28°C) warming. Compared to the control (17°C), growth of E. huxleyi increased under elevated temperatures, with higher rates recorded under moderate- relative to elevated- warming. Proteomic analysis revealed a significant modification of the E. huxleyi cellular proteome as temperatures increased: at lower temperature, ribosomal proteins and photosynthetic machinery appeared abundant, as rates of protein translation and photosynthetic performance are restricted by low temperatures. As temperatures increased, evidence of heat stress was observed in the photosystem, characterized by a relative down-regulation of the Photosystem II oxygen evolving complex and ATP synthase. Acclimation to elevated warming (28°C) revealed a substantial alteration to carbon metabolism. Here, E. huxleyi made use of the glyoxylate cycle and succinate metabolism to optimize carbon use, maintain growth and maximize ATP production in heat-damaged mitochondria, enabling cultures to maintain growth at levels significantly higher than those recorded in the control (17°C). Based on the metabolic changes observed, we can predict that warming may benefit photosynthetic carbon fixation by E. huxleyi in the sub-optimal to optimal thermal range. Past the thermal optima, increasing rates of respiration and costs of repair will likely constrain growth, causing a possible decline in the contribution of this species to the oceanic carbon sink depending on the evolvability of these temperature thresholds.

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.