Changes in gene expression of Prymnesium parvum induced by nitrogen and phosphorus limitation

Long-read genome sequencing provides novel insights into the harmful algal bloom species Prymnesium parvum

Prymnesium parvum is a toxin-producing haptophyte that causes harmful algal blooms worldwide, which are often associated with massive fish-kills and subsequent economic losses. In here, we present nuclear and plastid genome assemblies using PacBio HiFi long reads and DNBseq short reads for the two P. parvum strains UTEX 2797 and CCMP 3037, representing producers of type A prymnesins. Our results show that the P. parvum strains have a moderate haptophyte genome size of 97.56 and 107.32 Mb. The genome assemblies present one of highest contiguous assembled contig sequences to date consisting of 463 and 362 contigs with a contig N50 of 596.99 kb and 968.39 kb for strain UTEX 2797 and CCMP 3037, respectively. The assembled contigs of UTEX 2797 and CCMP 3037 were anchored to 34 scaffolds, with a scaffold N50 of 5.35 Mb and 3.61 Mb, respectively, accounting for 93.2 % and 97.9 % of the total length. Each plastid genome comprises a circular contig. A total of 20,578 and 19,426 protein-coding genes were annotated for UTEX 2797 and CCMP 3037. The expanded gene family analysis showed that starch and sucrose metabolism, sulfur metabolism, energy metabolism and ABC transporters are involved in the evolution of P. parvum. Polyketide synthase (PKS) genes responsible for the production of secondary metabolites such as prymnesins displayed different expression patterns under nutrient limitation. Overlap with repeats and horizontal gene transfer may be two contributing factors to the high number of PKS genes found in this species. The two high quality P. parvum genomes will serve as valuable resources for ecological, genetic, and toxicological studies of haptophytes that can be used to monitor and potentially manage harmful blooms of ichthyotoxic P. parvum in the future.

Phosphate Limitation Responses in Marine Green Algae Are Linked to Reprogramming of the tRNA Epitranscriptome and Codon Usage Bias

Marine algae are central to global carbon fixation, and their productivity is dictated largely by resource availability. Reduced nutrient availability is predicted for vast oceanic regions as an outcome of climate change; however, there is much to learn regarding response mechanisms of the tiny picoplankton that thrive in these environments, especially eukaryotic phytoplankton. Here, we investigate responses of the picoeukaryote Micromonas commoda, a green alga found throughout subtropical and tropical oceans. Under shifting phosphate availability scenarios, transcriptomic analyses revealed altered expression of transfer RNA modification enzymes and biased codon usage of transcripts more abundant during phosphate-limiting versus phosphate-replete conditions, consistent with the role of transfer RNA modifications in regulating codon recognition. To associate the observed shift in the expression of the transfer RNA modification enzyme complement with the transfer RNAs encoded by M. commoda, we also determined the transfer RNA repertoire of this alga revealing potential targets of the modification enzymes. Codon usage bias was particularly pronounced in transcripts encoding proteins with direct roles in managing phosphate limitation and photosystem-associated proteins that have ill-characterized putative functions in "light stress." The observed codon usage bias corresponds to a proposed stress response mechanism in which the interplay between stress-induced changes in transfer RNA modifications and skewed codon usage in certain essential response genes drives preferential translation of the encoded proteins. Collectively, we expose a potential underlying mechanism for achieving growth under enhanced nutrient limitation that extends beyond the catalog of up- or downregulated protein-encoding genes to the cell biological controls that underpin acclimation to changing environmental conditions.

Microcystis spp. and phosphorus in aquatic environments: A comprehensive review on their physiological and ecological interactions

Cyanobacterial blooms caused by eutrophication have become a major environmental problem in aquatic ecosystems worldwide over the last few decades. Phosphorus is a limiting nutrient that affects the growth of cyanobacteria and plays a role in dynamic changes in algal density and the formation of cyanobacterial blooms. Therefore, identifying the association between phosphorus sources and Microcystis, which is the most representative and harmful cyanobacteria, is essential for building an understanding of the ecological risks of cyanobacterial blooms. However, systematic reviews summarizing the relationships between Microcystis and phosphorus in aquatic environments are rare. Thus, this study provides a comprehensive overview of the physiological and ecological interactions between phosphorus sources and Microcystis in aquatic environments from the following perspectives: (i) the effects of phosphorus source and concentration on Microcystis growth, (ii) the impacts of phosphorus on the environmental behaviors of Microcystis, (iii) mechanisms of phosphorus-related metabolism in Microcystis, and (iv) role of Microcystis in the distribution of phosphorus sources within aquatic environments. In addition, relevant unsolved issues and essential future investigations (e.g., secondary ecological risks) have been highlighted and discussed. This review provides deeper insights into the relationship between phosphorus sources and Microcystis and can serve as a reference for the evaluation, monitoring, and effective control of cyanobacterial blooms.

Insights into Toxic Prymnesium parvum Blooms as a Cause of the Ecological Disaster on the Odra River

In 2022, Poland and Germany experienced a prolonged and extensive mass fish kill in the Odra River. During the period from the end of July to the beginning of September 2022, a high level of incidental disease and mortality was observed in various fish species (dozens of different species were found dead). The fish mortality affected five Polish provinces (Silesia, Opole, Lower Silesia, Lubuskie, and Western Pomerania) and involved reservoir systems covering most of the river (the Odra River is 854 km long, of which 742 km are in Poland). Fatal cases were investigated using toxicological, anatomopathological, and histopathological tests. Water samples were collected to determine nutrient status in the water column, phytoplankton biomass, and community composition. High nutrient concentrations indicated high phytoplankton productivity, with favorable conditions for golden algal blooms. The harmful toxins (prymnesins secreted by Prymnesium parvum habitats) had not been found in Poland before, but it was only a matter of time, especially in the Odra River, whose waters are permanently saline and still used for navigation. The observed fish mortality resulted in a 50% decrease in the fish population in the river and affected mainly cold-blooded species. Histopathological examinations of fish showed acute damage to the most perfused organs (gills, spleen, kidneys). The disruption to hematopoietic processes and damage to the gills were due to the action of hemolytic toxins (prymnesins). An evaluation of the collected hydrological, meteorological, biological, and physico-chemical data on the observed spatio-temporal course of the catastrophe, as well as the detection of three compounds from the group of B-type prymnesins in the analyzed material (the presence of prymnesins was confirmed using an analysis of the fragmentation spectrum and the accurate tandem mass spectrometry (MS/MS) measurement, in combination with high-performance liquid chromatography-tandem mass spectrometry (LC-MS/MS), allowed the formulation and subsequent testing of the hypothesis for a direct link between the observed fish mortality and the presence of prymnesins in the Odra River. This article systematizes what is known about the causes of the fish kill in the Odra River in 2022, based on official government reports (one Polish and one German) and the EU technical report by the Joint Research Centre. A review and critical analysis of government findings (Polish and German) on this disaster were conducted in the context of what is known to date about similar cases of mass fish kills.

The Effect of pH and Salinity on the Toxicity and Growth of the Golden Alga, Prymnesium parvum

Bioassays using cultures of the toxic haptophyte Prymnesium parvum and the ciliate Cyclidium sp. as prey were conducted to test the effect of pH (range = 6.5 - 8.5), salinity (range = 1.50 - 7.50‰), and a combination of pH and salinity on the toxicity of P. parvum. pH had a significant effect on P. parvum toxicity. Toxicity was rapidly (within 24 hr) induced by increasing pH of the medium, or reduced by lowering pH. Conversely, lowering salinity reduced toxicity, albeit less effectively compared to pH, and P. parvum cells remained toxic at the lowest values tested (1.50‰ at pH 7.5). An additional effect between pH and salinity was also observed: low salinity combined with low pH led to not only decreased toxicity, but also resulted in lower P. parvum growth rates. Such effects of pH and salinity on P. parvum growth and toxicity provide insight into the environmental factors supporting community dominance and toxic blooms of the alga.

Mixotrophy in the bloom-forming genus Phaeocystis and other haptophytes

Phaeocystis is a globally widespread marine phytoplankton genus, best known for its colony-forming species that can form large blooms and odorous foam during bloom decline. In the North Sea, Phaeocystis globosa typically becomes abundant towards the end of the spring bloom, when nutrients are depleted and the share of mixotrophic protists increases. Although mixotrophy is widespread across the eukaryotic tree of life and is also found amongst haptophytes, a mixotrophic nutrition has not yet been demonstrated in Phaeocystis. Here, we sampled two consecutive Phaeocystis globosa spring blooms in the coastal North Sea. In both years, bacterial cells were observed inside 0.6 - 2% of P. globosa cells using double CARD-FISH hybridizations in combination with laser scanning confocal microscopy. Incubation experiments manipulating light and nutrient availability showed a trend towards higher occurrence of intracellular bacteria under P-deplete conditions. Based on counts of bacteria inside P. globosa cells in combination with theoretical values of prey digestion times, maximum ingestion rates of up to 0.08 bacteria cell^-1 h^-1 were estimated. In addition, a gene-based predictive model was applied to the transcriptome assemblies of seven Phaeocystis strains and 24 other haptophytes to assess their trophic mode. This model predicted a phago-mixotrophic feeding strategy in several (but not all) strains of P. globosa, P. antarctica and other haptophytes that were previously assumed to be autotrophic. The observation of bacterial cells inside P. globosa and the gene-based model predictions strongly suggest that the phago-mixotrophic feeding strategy is widespread among members of the Phaeocystis genus and other haptophytes, and might contribute to their remarkable success to form nuisance blooms under nutrient-limiting conditions.

Transcriptomic response of the harmful algae Heterosigma akashiwo to polyphosphate utilization and phosphate stress

Phosphorus (P) is one of the major macronutrients necessary for phytoplankton growth. In some parts of the ocean, however, P is frequently scarce, hence, there is limited phytoplankton growth. To cope with P deficiency, phytoplankton evolved a variety of strategies, including, utilization of different P sources. Polyphosphate (polyP) is ubiquitously present and serves an essential function in aquatic environments, but it is unclear if and how this polymer is utilized by phytoplankton. Here, we examined the physiological and molecular responses of the widely present harmful algal bloom (HAB) species, Heterosigma akashiwo in polyP utilization, and in coping with P-deficiency. Our results revealed that two forms of inorganic polyP, namely, sodium tripolyphosphate and sodium hexametaphosphate, support H. akashiwo growth as efficiently as orthophosphate. However, few genes involved in polyP utilization have been identified. Under P-deficient conditions, genes associated with P transport, dissolved organic P utilization, sulfolipid synthesis, and energy production, were markedly elevated. In summary, our results indicate that polyP is bioavailable to H. akashiwo, and this HAB species have evolved a comprehensive strategy to cope with P deficiency.

Transcriptomic responses to phosphorus in an invasive cyanobacterium, Raphidiopsis raciborskii: Implications for nutrient management

Phosphorus (P) is a vital macronutrient associated with the growth and proliferation of Raphidiopsis raciborskii, an invasive and notorious bloom-forming cyanobacterium. However, the molecular mechanisms involved in P acclimation remain largely unexplored for Raphidiopsis raciborskii. Here, transcriptome sequencing of Raphidiopsis raciborskii was conducted to reveal multifaceted mechanisms involved in mimicking dipotassium phosphate (DIP), β-glycerol phosphate (Gly), 2-aminoethylphosphonic acid (AEP), and P-free conditions (NP). Chlorophyll a fluorescence parameters showed significant differences in the NP and AEP groups compared with the DIP and Gly-groups. Expression levels of genes related to phosphate transportation and uptake, organic P utilization, nitrogen metabolism, urea cycling, carbon fixation, amino acid metabolism, environmental information, the ATP-synthesis process in glycolysis, the tricarboxylic acid (TCA) cycle, and the pentose phosphate pathway were remarkably upregulated, while those related to photosynthesis, phycobiliproteins, respiration, oxidative phosphorylation, sulfur metabolism, and genetic information were markedly downregulated in the NP group relative to the DIP group. However, the expression of genes involved in organic P utilization, the urea cycle, and genetic information in the Gly-group, and carbon-phosphorus lyase, genetic information and environmental information in the AEP group were significantly increased compared to the DIP group. Together, these results indicate that Raphidiopsis raciborskii exhibits the evolution of coordination of multiple metabolic pathways and certain key genes to adapt to ambient P changes, which implies that if P is reduced to control Raphidiopsis raciborskii bloom, there is a risk that external nutrients (such as nitrogen, amino acids, and urea) will stimulate the growth or metabolism of Raphidiopsis.

No evidence of Phago-mixotropy in Micromonas polaris (Mamiellophyceae), the Dominant Picophytoplankton Species in the Arctic

In the Arctic Ocean, the small green alga Micromonas polaris dominates picophytoplankton during the summer months but is also present in winter. It has been previously hypothesized to be phago-mixotrophic (capable of bacteria ingestion) based on laboratory and field experiments. Prey uptake was analyzed in several M. polaris strains isolated from different regions and depths of the Arctic Ocean and in Ochromonas triangulata, a known phago-mixotroph used as a control. Measuring ingestion of either fluorescent beads or fluorescently labeled bacteria by flow cytometry, we found no evidence of phago-mixotrophy in any M. polaris strain while O. triangulata was ingesting both beads and bacteria. In addition, in silico predictions revealed that members of the genus Micromonas lack a genetic signature of phagocytotic capacity.

In-Situ Metatranscriptomic Analyses Reveal the Metabolic Flexibility of the Thermophilic Anoxygenic Photosynthetic Bacterium Chloroflexus aggregans in a Hot Spring Cyanobacteria-Dominated Microbial Mat

Chloroflexus aggregans is a metabolically versatile, thermophilic, anoxygenic phototrophic member of the phylum Chloroflexota (formerly Chloroflexi), which can grow photoheterotrophically, photoautotrophically, chemoheterotrophically, and chemoautotrophically. In hot spring-associated microbial mats, C. aggregans co-exists with oxygenic cyanobacteria under dynamic micro-environmental conditions. To elucidate the predominant growth modes of C. aggregans, relative transcription levels of energy metabolism- and CO₂ fixation-related genes were studied in Nakabusa Hot Springs microbial mats over a diel cycle and correlated with microscale in situ measurements of O₂ and light. Metatranscriptomic analyses indicated two periods with different modes of energy metabolism of C. aggregans: (1) phototrophy around midday and (2) chemotrophy in the early morning hours. During midday, C. aggregans mainly employed photoheterotrophy when the microbial mats were hyperoxic (400–800 µmol L⁻¹ O₂). In the early morning hours, relative transcription peaks of genes encoding uptake hydrogenase, key enzymes for carbon fixation, respiratory complexes as well as enzymes for TCA cycle and acetate uptake suggest an aerobic chemomixotrophic lifestyle. This is the first in situ study of the versatile energy metabolism of C. aggregans based on gene transcription patterns. The results provide novel insights into the metabolic flexibility of these filamentous anoxygenic phototrophs that thrive under dynamic environmental conditions.

Transcriptome responses of the dinoflagellate Karenia mikimotoi driven by nitrogen deficiency

The availability of ambient N nutrient is often correlated with the occurrences of harmful algal bloom formed by certain dinoflagellates, making it important to understand how these species might be responding to such conditions. Here, transcriptome sequencing of Karenia mikimotoi was conducted to understand the underlying molecular mechanisms by which this dinoflagellate copes with nitrogen (N) deficiency. Transcriptomic analysis revealed 8802 unigenes (3.56%) that were differentially expressed with ≥ 2-fold change. Under N-depleted conditions, genes involved in glycolysis, fatty acid metabolism, and the tricarboxylic acid (TCA) cycle as well as lipid accumulation were significantly upregulated. The elevated expression of enzymes used in protein degradation and turnover suggests possible metabolic reconfiguration towards accelerated N recycling. Moreover, a significant increase in urea transporter was observed, indicating increased assimilation of organic nitrogen resources as an alternative in N-depleted cultures of K. mikimotoi. The down-regulated glutamate synthase genes were also identified under N deficiency, suggesting suppression of primary amino acid synthesis to save N resource. Taken together, results of this study show enhanced multiple N resource acquisition and reuse of multiple N resources constitute a comprehensive strategy to cope with N deficiency in a dinoflagellate.

Integrated omics unveil the secondary metabolic landscape of a basal dinoflagellate

BACKGROUND

Some dinoflagellates cause harmful algal blooms, releasing toxic secondary metabolites, to the detriment of marine ecosystems and human health. Our understanding of dinoflagellate toxin biosynthesis has been hampered by their unusually large genomes. To overcome this challenge, for the first time, we sequenced the genome, microRNAs, and mRNA isoforms of a basal dinoflagellate, Amphidinium gibbosum, and employed an integrated omics approach to understand its secondary metabolite biosynthesis.

RESULTS

We assembled the ~ 6.4-Gb A. gibbosum genome, and by probing decoded dinoflagellate genomes and transcriptomes, we identified the non-ribosomal peptide synthetase adenylation domain as essential for generation of specialized metabolites. Upon starving the cells of phosphate and nitrogen, we observed pronounced shifts in metabolite biosynthesis, suggestive of post-transcriptional regulation by microRNAs. Using Iso-Seq and RNA-seq data, we found that alternative splicing and polycistronic expression generate different transcripts for secondary metabolism.

CONCLUSIONS

Our genomic findings suggest intricate integration of various metabolic enzymes that function iteratively to synthesize metabolites, providing mechanistic insights into how dinoflagellates synthesize secondary metabolites, depending upon nutrient availability. This study provides insights into toxin production associated with dinoflagellate blooms. The genome of this basal dinoflagellate provides important clues about dinoflagellate evolution and overcomes the large genome size, which has been a challenge previously.

Utilization of various forms of nitrogen and expression regulation of transporters in the harmful alga Heterosigma akashiwo (Raphidophyceae)

Nitrogen (N) is an essential nutrient for phytoplankton growth. There is ample evidence that N enrichment promotes harmful algae blooms (HABs) but molecular mechanisms regulating N-nutrient uptake and metabolism are not so clear, especially for the raphidophyte Heterosigma akashiwo, which forms ichthyotoxic HABs in many coastal waters. In this study, the utilization of three different chemical forms of N (nitrate, ammonium, and urea) by H. akashiwo CCMA 369 was investigated in batch culture conditions. Results showed that H. akashiwo grew well on all three N compounds, and the highest cell yield occurred in the NH₄⁺ culture group. Reverse transcription quantitative PCR analysis revealed that the expression of high-affinity NO₃^- transporter (NIT), NH₄⁺ transporter (AMT) and high-affinity urea active transporter (DUR3), were significantly up-regulated under N-limitation compared to the N-replete control. The mRNA levels of AMT and DUR3 also displayed a clear diel rhythm, with high levels at midnight. In addition, NH₄⁺ addition (5 μM) did not depress the transcript abundance of any of the three N transporters. Compared with the co-occurring immobile diatom Skeletonema costatum, the high expression of AMT in dark period in H. akashiwo is consistent with its diel vertical migration behavior, which may promote N-nutrient acquisition from deeper layers and give advantages for H. akashiwo to form blooms.

Transcriptional and physiological responses to inorganic nutrition in a tropical Pacific strain of Alexandrium minutum: Implications for nutrient uptakes and assimilation

The marine dinoflagellate Alexandrium minutum is known to produce saxitoxins that cause paralytic shellfish poisoning in human worldwide through consumption of the contaminated shellfish mollusks. Despite numerous studies on the growth physiology and saxitoxin production of this species, the knowledge on the molecular basis of nutrient uptakes in relation to toxin production in this species is limited. In this study, relative expressions of the high-affinity transporter genes of nitrate, ammonium, and phosphate (AmNrt2, AmAmt1 and AmPiPT1) and the assimilation genes, nitrate reductase (AmNas), glutamine synthase (AmGSIII) and carbamoyl phosphate synthase (AmCPSII) from A. minutum were studied in batch clonal culture condition with two nitrogen sources (nitrate: NO₃^- or ammonium: NH₄⁺) under different N:P ratios (high-P: N:P of 14 and 16, and low-P: N:P of 155). The expression of AmAmt1 was suppressed in excess NH₄⁺-grown condition but was not observed in AmNrt2 and AmNas. Expressions of AmAmt1, AmNrt2, AmNas, AmGSIII, AmCPSII, and AmPiPT1 were high in P-deficient condition, showing that A. minutum is likely to take up nutrients for growth under P-stress condition. Conversely, relative expression of AmCPSII was incongruent with cell growth, but was well correlated with toxin quota, suggesting that the gene might involve in arginine metabolism and related toxin production pathway. The expression of AmGSIII is found coincided with higher toxin production and is believed to involve in mechanism to detoxify the cells from excess ammonium stress. The gene regulation observed in this study has provided better insights into the ecophysiology of A. minutum in relation to its adaptive strategies in unfavorable environments.

Temporal phosphate gradients reveal diverse acclimation responses in phytoplankton phosphate uptake

Phytoplankton face environmental nutrient variations that occur in the dynamic upper layers of the ocean. Phytoplankton cells are able to rapidly acclimate to nutrient fluctuations by adjusting their nutrient-uptake system and metabolism. Disentangling these acclimation responses is a critical step in bridging the gap between phytoplankton cellular physiology and community ecology. Here, we analyzed the dynamics of phosphate (P) uptake acclimation responses along different P temporal gradients by using batch cultures of the diatom Phaeodactylum tricornutum. We employed a multidisciplinary approach that combined nutrient-uptake bioassays, transcriptomic analysis, and mathematical models. Our results indicated that cells increase their maximum nutrient-uptake rate (V_max) both in response to P pulses and strong phosphorus limitation. The upregulation of three genes coding for different P transporters in cells experiencing low intracellular phosphorus levels supported some of the observed V_max variations. In addition, our mathematical model reproduced the empirical V_max patterns by including two types of P transporters upregulated at medium-high environmental and low intracellular phosphorus levels, respectively. Our results highlight the existence of a sequence of acclimation stages along the phosphate continuum that can be understood as a succession of acclimation responses. We provide a novel conceptual framework that can contribute to integrating and understanding the dynamics and wide diversity of acclimation responses developed by phytoplankton.

Effects of phosphorus sources on volatile organic compound emissions from Microcystis flos-aquae and their toxic effects on Chlamydomonas reinhardtii

There is diverse phosphorus (P) in eutrophicated waters, but it is considered as a crucial nutrient for cyanobacteria growth due to its easy precipitation as insoluble salts. To uncover the effects of complex P nutrients on the emission of volatile organic compounds (VOCs) from cyanobacteria and their toxic effects on other algae, the VOCs from Microcystis flos-aquae supplied with different types and amount of P nutrients were analyzed, and the effects of VOCs and their two main compounds on Chlamydomonas reinhardtii growth were investigated. When M. flos-aquae cells were supplied with K₂HPO₄, sodium pyrophosphate and sodium hexametaphosphate as the sole P source, 27, 23 and 29 compounds were found, respectively, including furans, sulfocompounds, terpenoids, benzenes, aldehydes, hydrocarbons and esters. With K₂HPO₄ as the sole P source, the VOC emission increased with reducing P amount, and the maximum emission was found under Non-P condition. In the treatments of M. flos-aquae VOCs under Non-P condition and two main terpenoids (eucalyptol and limonene) in the VOCs, remarkable decreases were found in C. reinhardtii cell growth, photosynthetic pigment content and photosynthetic abilities. Therefore, we deduce that multiple P nutrients in eutrophicated waters induce different VOC emissions from cyanobacteria, and P amount reduction caused by natural precipitation and algal massive growth results in more VOC emissions. These VOCs play toxic roles in cyanobacteria becoming dominant species, and eucalyptol and limonene are two toxic agents.

Shifting metabolic priorities among key protistan taxa within and below the euphotic zone

A metatranscriptome study targeting the protistan community was conducted off the coast of Southern California, at the San Pedro Ocean Time-series station at the surface, 150 m (oxycline), and 890 m to link putative metabolic patterns to distinct protistan lineages. Comparison of relative transcript abundances revealed depth-related shifts in the nutritional modes of key taxonomic groups. Eukaryotic gene expression in the sunlit surface environment was dominated by phototrophs, such as diatoms and chlorophytes, and high abundances of transcripts associated with synthesis pathways (e.g., photosynthesis, carbon fixation, fatty acid synthesis). Sub-euphotic depths (150 and 890 m) exhibited strong contributions from dinoflagellates and ciliates, and were characterized by transcripts relating to digestion or intracellular nutrient recycling (e.g., breakdown of fatty acids and V-type ATPases). These transcriptional patterns underlie the distinct nutritional modes of ecologically important protistan lineages that drive marine food webs, and provide a framework to investigate trophic dynamics across diverse protistan communities.

Eukaryotic phytoplankton community spatiotemporal dynamics as identified through gene expression within a eutrophic estuary

Over the span of a year, we investigated the interactions between biotic and abiotic factors within the eutrophic Neuse River Estuary (NRE). Through metatranscriptomic sequencing in combination with water quality measurements, we show that there are different metabolic strategies deployed along the NRE. In the upper estuary, taxonomically resolved phytoplankton groups express more transcripts of genes for synthesis of cellular components and carbon metabolism whereas in the lower estuary, transcripts allocated to nutrient metabolism and transport were more highly expressed. Metabolisms for polysaccharide synthesis and transportation were elevated in the lower estuary and could be reflective of unbalanced growth and/or interactions with their surrounding microbial consortia. Our results indicate phytoplankton have high metabolic activity, suggestive of increased growth rates in the upper estuary and display patterns reflective of nutrient limitation in the lower estuary. Among all the environmental parameters varying along the NRE, nitrogen availability is found to be the main driving factor for the observed spatial divergence.

Utilization of urea and expression profiles of related genes in the dinoflagellate Prorocentrum donghaiense

Urea has been shown to contribute more than half of total nitrogen (N) required by phytoplankton in some estuaries and coastal waters and to provide a substantial portion of the N demand for many harmful algal blooms (HABs) of dinoflagellates. In this study, we investigated the physiological and transcriptional responses in Prorocentrum donghaiense to changes in nitrate and urea availability. We found that this species could efficiently utilize urea as sole N source and achieve comparable growth rate and photosynthesis capability as it did under nitrate. These physiological parameters were markedly lower in cultures grown under nitrate- or urea-limited conditions. P. donghaiense N content was similarly low under nitrate- or urea-limited culture condition, but was markedly higher under urea-replete condition than under nitrate-replete condition. Carbon (C) content was consistently elevated under N-limited condition. Consequently, the C:N ratio was as high as 21:1 under nitrate- or urea-limitation, but 7:1 under urea-replete condition and 9:1 to 10:1 under nitrate-replete condition. Using quantitative reverse transcription PCR, we investigated the expression pattern for four genes involved in N transport and assimilation. The results indicated that genes encoding nitrate transport, urea hydrolysis, and nickel transporter gene were sensitive to changes in general N nutrient availability whereas the urea transporter gene responded much more strongly to changes in urea concentration. Taken together, our study shows the high bioavailability of urea, its impact on C:N stoichiometry, and the sensitivity of urea transporter gene expression to urea availability.

Transcriptional response of the harmful raphidophyte Heterosigma akashiwo to nitrate and phosphate stress

The marine eukaryotic alga Heterosigma akashiwo (Raphidophyceae) is known for forming ichthyotoxic harmful algal blooms (HABs). In the past 50 years, H. akashiwo blooms have increased, occurring globally in highly eutrophic coastal and estuarine systems. These systems often incur dramatic physicochemical changes, including macronutrient (nitrogen and phosphorus) enrichment and depletion, on short timescales. Here, H. akashiwo cultures grown under nutrient replete, low N and low P growth conditions were examined for changes in biochemical and physiological characteristics in concert with transcriptome sequencing to provide a mechanistic perspective on the metabolic processes involved in responding to N and P stress. There was a marked difference in the overall transcriptional pattern between low N and low P transcriptomes. Both nutrient stresses led to significant changes in the abundance of thousands of contigs related to a wide diversity of metabolic pathways, with limited overlap between the transcriptomic responses to low N and low P. Enriched contigs under low N included many related to nitrogen metabolism, acquisition, and transport. In addition, metabolic modules like photosynthesis and carbohydrate metabolism changed significantly under low N, coincident with treatment-specific changes in photosynthetic efficiency and particulate carbohydrate content. P-specific contigs responsible for P transport and organic P use were more enriched in the low P treatment than in the replete control and low N treatment. These results provide new insight into the genetic mechanisms that distinguish how this HAB species responds to these two common nutrient stresses, and the results can inform future field studies, linking transcriptional patterns to the physiological ecology of H. akashiwo in situ.

Conserved Transcriptional Responses to Nutrient Stress in Bloom-Forming Algae

The concentration and composition of bioavailable nitrogen (N) and phosphorus (P) in the upper ocean shape eukaryotic phytoplankton communities and influence their physiological responses. Phytoplankton are known to exhibit similar physiological responses to limiting N and P conditions such as decreased growth rates, chlorosis, and increased assimilation of N and P. Are these responses similar at the molecular level across multiple species? To interrogate this question, five species from biogeochemically important, bloom-forming taxa (Bacillariophyta, Dinophyta, and Haptophyta) were grown under similar low N, low P, and replete nutrient conditions to identify transcriptional patterns and associated changes in biochemical pools related to N and P stress. Metabolic profiles, revealed through the transcriptomes of these taxa, clustered together based on species rather than nutrient stressor, suggesting that the global metabolic response to nutrient stresses was largely, but not exclusively, species-specific. Nutrient stress led to few transcriptional changes in the two dinoflagellates, consistent with other research. An orthologous group analysis examined functionally conserved (i.e., similarly changed) responses to nutrient stress and therefore focused on the diatom and haptophytes. Most conserved ortholog changes were specific to a single nutrient treatment, but a small number of orthologs were similarly changed under both N and P stress in 2 or more species. Many of these orthologs were related to photosynthesis and may represent generalized stress responses. A greater number of orthologs were conserved across more than one species under low P compared to low N. Screening the conserved orthologs for functions related to N and P metabolism revealed increased relative abundance of orthologs for nitrate, nitrite, ammonium, and amino acid transporters under N stress, and increased relative abundance of orthologs related to acquisition of inorganic and organic P substrates under P stress. Although the global transcriptional responses were dominated by species-specific changes, the analysis of conserved responses revealed functional similarities in resource acquisition pathways among different phytoplankton taxa. This overlap in nutrient stress responses observed among species may be useful for tracking the physiological ecology of phytoplankton field populations.

Gene expression characterizes different nutritional strategies among three mixotrophic protists

Mixotrophic protists, i.e. protists that can carry out both phototrophy and heterotrophy, are a group of organisms with a wide range of nutritional strategies. The ecological and biogeochemical importance of these species has recently been recognized. In this study, we investigated and compared the gene expression of three mixotrophic protists, Prymnesium parvum, Dinobyron sp. and Ochromonas sp. under light and dark conditions in the presence of prey using RNA-Seq. Gene expression of the obligately phototrophic P. parvum and Dinobryon sp. changed significantly between light and dark treatments, while that of primarily heterotrophic Ochromonas sp. was largely unchanged. Gene expression of P. parvum and Dinobryon sp. shared many similarities, especially in the expression patterns of genes related to reproduction. However, key genes involved in central carbon metabolism and phagotrophy had different expression patterns between these two species, suggesting differences in prey consumption and heterotrophic nutrition in the dark. Transcriptomic data also offered clues to other physiological traits of these organisms such as preference of nitrogen sources and photo-oxidative stress. These results provide potential target genes for further exploration of the mechanisms of mixotrophic physiology and demonstrate the potential usefulness of molecular approaches in characterizing the nutritional modes of mixotrophic protists.