Coastal Antarctic marine ecosystems are significant in carbon cycling because of their intense seasonal phytoplankton blooms. Southern Ocean algae are primarily limited by light and iron (Fe) and can be co-limited by cobalamin (vitamin B12). Micronutrient limitation controls productivity and shapes the composition of blooms which are typically dominated by either diatoms or the haptophyte Phaeocystis antarctica. However, the vitamin requirements and ecophysiology of the keystone species P. antarctica remain poorly characterized. Using cultures, physiological analysis, and comparative omics, we examined the response of P. antarctica to a matrix of Fe-B12 conditions. We show that P. antarctica is not auxotrophic for B12, as previously suggested, and identify mechanisms underlying its B12 response in cultures of predominantly solitary and colonial cells. A combination of proteomics and proteogenomics reveals a B12-independent methionine synthase fusion protein (MetE-fusion) that is expressed under vitamin limitation and interreplaced with the B12-dependent isoform under replete conditions. Database searches return homologues of the MetE-fusion protein in multiple Phaeocystis species and in a wide range of marine microbes, including other photosynthetic eukaryotes with polymorphic life cycles as well as bacterioplankton. Furthermore, we find MetE-fusion homologues expressed in metaproteomic and metatranscriptomic field samples in polar and more geographically widespread regions. As climate change impacts micronutrient availability in the coastal Southern Ocean, our finding that P. antarctica has a flexible B12 metabolism has implications for its relative fitness compared to B12-auxotrophic diatoms and for the detection of B12-stress in a more diverse set of marine microbes.

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. closteriums 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. closteriums metabolic flexibility and its potential role in global carbon cycling.

Diatoms outcompete other phytoplankton for nitrate, yet little is known about the mechanisms underpinning this ability. Genomes and genome-enabled studies have shown that diatoms possess unique features of nitrogen metabolism however, the implications for nutrient utilization and growth are poorly understood. Using a combination of transcriptomics, proteomics, metabolomics, fluxomics, and flux balance analysis to examine short-term shifts in nitrogen utilization in the model pennate diatom in Phaeodactylum tricornutum, we obtained a systems-level understanding of assimilation and intracellular distribution of nitrogen. Chloroplasts and mitochondria are energetically integrated at the critical intersection of carbon and nitrogen metabolism in diatoms. Pathways involved in this integration are organelle-localized GS-GOGAT cycles, aspartate and alanine systems for amino moiety exchange, and a split-organelle arginine biosynthesis pathway that clarifies the role of the diatom urea cycle. This unique configuration allows diatoms to efficiently adjust to changing nitrogen status, conferring an ecological advantage over other phytoplankton taxa.

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 coexpression 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 colimited the diatom population along with low Si. Iron limitation represents an overlooked 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.

Microbial communities in the world ocean are affected strongly by oceanic circulation, creating characteristic marine biomes. The high connectivity of most of the ocean makes it difficult to disentangle selective retention of colonizing genotypes (with traits suited to biome specific conditions) from evolutionary selection, which would act on founder genotypes over time. The Arctic Ocean is exceptional with limited exchange with other oceans and ice covered since the last ice age. To test whether Arctic microalgal lineages evolved apart from algae in the global ocean, we sequenced four lineages of microalgae isolated from Arctic waters and sea ice. Here we show convergent evolution and highlight geographically limited HGT as an ecological adaptive force in the form of PFAM complements and horizontal acquisition of key adaptive genes. Notably, ice-binding proteins were acquired and horizontally transferred among Arctic strains. A comparison with Tara Oceans metagenomes and metatranscriptomes confirmed mostly Arctic distributions of these IBPs. The phylogeny of Arctic-specific genes indicated that these events were independent of bacterial-sourced HGTs in Antarctic Southern Ocean microalgae.