Development of a molecular-based index for assessing iron status in bloom-forming pennate diatoms

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.

Whole-genome scanning reveals environmental selection mechanisms that shape diversity in populations of the epipelagic diatom Chaetoceros

Diatoms form a diverse and abundant group of photosynthetic protists that are essential players in marine ecosystems. However, the microevolutionary structure of their populations remains poorly understood, particularly in polar regions. Exploring how closely related diatoms adapt to different environments is essential given their short generation times, which may allow rapid adaptations, and their prevalence in marine regions dramatically impacted by climate change, such as the Arctic and Southern Oceans. Here, we address genetic diversity patterns in Chaetoceros, the most abundant diatom genus and one of the most diverse, using 11 metagenome-assembled genomes (MAGs) reconstructed from Tara Oceans metagenomes. Genome-resolved metagenomics on these MAGs confirmed a prevalent distribution of Chaetoceros in the Arctic Ocean with lower dispersal in the Pacific and Southern Oceans as well as in the Mediterranean Sea. Single-nucleotide variants identified within the different MAG populations allowed us to draw a landscape of Chaetoceros genetic diversity and revealed an elevated genetic structure in some Arctic Ocean populations. Gene flow patterns of closely related Chaetoceros populations seemed to correlate with distinct abiotic factors rather than with geographic distance. We found clear positive selection of genes involved in nutrient availability responses, in particular for iron (e.g., ISIP2a, flavodoxin), silicate, and phosphate (e.g., polyamine synthase), that were further supported by analysis of Chaetoceros transcriptomes. Altogether, these results highlight the importance of environmental selection in shaping diatom diversity patterns and provide new insights into their metapopulation genomics through the integration of metagenomic and environmental data.

Nitrogen and Iron Availability Drive Metabolic Remodeling and Natural Selection of Diverse Phytoplankton during Experimental Upwelling

Nearly half of carbon fixation and primary production originates from marine phytoplankton, and much of it occurs in episodic blooms in upwelling regimes. Here, we simulated blooms limited by nitrogen and iron by incubating Monterey Bay surface waters with subnutricline waters and inorganic nutrients and measured the whole-community transcriptomic response during mid- and late-bloom conditions. Cell counts revealed that centric and pennate diatoms (largely Pseudo-nitzschia and Chaetoceros spp.) were the major blooming taxa, but dinoflagellates, prasinophytes, and prymnesiophytes also increased. Viral mRNA significantly increased in late bloom and likely played a role in the bloom's demise. We observed conserved shifts in the genetic similarity of phytoplankton populations to cultivated strains, indicating adaptive population-level changes in community composition. Additionally, the density of single nucleotide variants (SNVs) declined in late-bloom samples for most taxa, indicating a loss of intraspecific diversity as a result of competition and a selective sweep of adaptive alleles. We noted differences between mid- and late-bloom metabolism and differential regulation of light-harvesting complexes (LHCs) under nutrient stress. While most LHCs are diminished under nutrient stress, we showed that diverse taxa upregulated specialized, energy-dissipating LHCs in low iron. We also suggest the relative expression of NRT2 compared to the expression of GSII as a marker of cellular nitrogen status and the relative expression of iron starvation-induced protein genes (ISIP1, ISIP2, and ISIP3) compared to the expression of the thiamine biosynthesis gene (thiC) as a marker of iron status in natural diatom communities. IMPORTANCE Iron and nitrogen are the nutrients that most commonly limit phytoplankton growth in the world's oceans. The utilization of these resources by phytoplankton sets the biomass available to marine systems and is of particular interest in high-nutrient, low-chlorophyll (HNLC) coastal fisheries. Previous research has described the biogeography of phytoplankton in HNLC regions and the transcriptional responses of representative taxa to nutrient limitation. However, the differential transcriptional responses of whole phytoplankton communities to iron and nitrogen limitation has not been previously described, nor has the selective pressure that these competitive bloom environments exert on major players. In addition to describing changes in the physiology of diverse phytoplankton, we suggest practical indicators of cellular nitrogen and iron status for future monitoring.

Diminished carbon and nitrate assimilation drive changes in diatom elemental stoichiometry independent of silicification in an iron-limited assemblage

In the California Current Ecosystem, upwelled water low in dissolved iron (Fe) can limit phytoplankton growth, altering the elemental stoichiometry of the particulate matter and dissolved macronutrients. Iron-limited diatoms can increase biogenic silica (bSi) content >2-fold relative to that of particulate organic carbon (C) and nitrogen (N), which has implications for carbon export efficiency given the ballasted nature of the silica-based diatom cell wall. Understanding the molecular and physiological drivers of this altered cellular stoichiometry would foster a predictive understanding of how low Fe affects diatom carbon export. In an artificial upwelling experiment, water from 96 m depth was incubated shipboard and left untreated or amended with dissolved Fe or the Fe-binding siderophore desferrioxamine-B (+DFB) to induce Fe-limitation. After 120 h, diatoms dominated the communities in all treatments and displayed hallmark signatures of Fe-limitation in the +DFB treatment, including elevated particulate Si:C and Si:N ratios. Single-cell, taxon-resolved measurements revealed no increase in bSi content during Fe-limitation despite higher transcript abundance of silicon transporters and silicanin-1. Based on these findings we posit that the observed increase in bSi relative to C and N was primarily due to reductions in C fixation and N assimilation, driven by lower transcript expression of key Fe-dependent genes.

In contrast to diatoms, cryptophytes are susceptible to iron limitation, but not to ocean acidification

Previous field studies in the Southern Ocean (SO) indicated an increased occurrence and dominance of cryptophytes over diatoms due to climate change. To gain a better mechanistic understanding of how the two ecologically important SO phytoplankton groups cope with ocean acidification (OA) and iron (Fe) availability, we chose two common representatives of Antarctic waters, the cryptophyte Geminigera cryophila and the diatom Pseudo-nitzschia subcurvata. Both species were grown at 2°C under different pCO₂ (400 vs. 900 μatm) and Fe (0.6 vs. 1.2 nM) conditions. For P. subcurvata, an additional high pCO₂ level was applied (1400 μatm). At ambient pCO₂ under low Fe supply, growth of G. cryophila almost stopped while it remained unaffected in P. subcurvata. Under high Fe conditions, OA was not beneficial for P. subcurvata, but stimulated growth and carbon production of G. cryophila. Under low Fe supply, P. subcurvata coped much better with OA than the cryptophyte, but invested more energy into photoacclimation. Our study reveals that Fe limitation was detrimental for the growth of G. cryophila and suppressed the positive OA effect. The diatom was efficient in coping with low Fe, but was stressed by OA while both factors together strongly impacted its growth. The distinct physiological response of both species to OA and Fe limitation explains their occurrence in the field. Based on our results, Fe availability is an important modulator of OA effects on SO phytoplankton, with different implications on the occurrence of cryptophytes and diatoms in the future.

Dinoflagellates alter their carbon and nutrient metabolic strategies across environmental gradients in the central Pacific Ocean

Marine microeukaryotes play a fundamental role in biogeochemical cycling through the transfer of energy to higher trophic levels and vertical carbon transport. Despite their global importance, microeukaryote physiology, nutrient metabolism and contributions to carbon cycling across offshore ecosystems are poorly characterized. Here, we observed the prevalence of dinoflagellates along a 4,600-km meridional transect extending across the central Pacific Ocean, where oligotrophic gyres meet equatorial upwelling waters rich in macronutrients yet low in dissolved iron. A combined multi-omics and geochemical analysis provided a window into dinoflagellate metabolism across the transect, indicating a continuous taxonomic dinoflagellate community that shifted its functional transcriptome and proteome as it extended from the euphotic to the mesopelagic zone. In euphotic waters, multi-omics data suggested that a combination of trophic modes were utilized, while mesopelagic metabolism was marked by cytoskeletal investments and nutrient recycling. Rearrangement in nutrient metabolism was evident in response to variable nitrogen and iron regimes across the gradient, with no associated change in community assemblage. Total dinoflagellate proteins scaled with particulate carbon export, with both elevated in equatorial waters, suggesting a link between dinoflagellate abundance and total carbon flux. Dinoflagellates employ numerous metabolic strategies that enable broad occupation of central Pacific ecosystems and play a dual role in carbon transformation through both photosynthetic fixation in the euphotic zone and remineralization in the mesopelagic zone.

Genomic mosaicism underlies the adaptation of marine Synechococcus ecotypes to distinct oceanic iron niches

Phytoplankton are limited by iron (Fe) in ~40% of the world's oceans including high-nutrient low-chlorophyll (HNLC) regions. While low-Fe adaptation has been well-studied in large eukaryotic diatoms, less is known for small, prokaryotic marine picocyanobacteria. This study reveals key physiological and genomic differences underlying Fe adaptation in marine picocyanobacteria. HNLC ecotype CRD1 strains have greater physiological tolerance to low Fe congruent with their expanded repertoire of Fe transporter, storage and regulatory genes compared to other ecotypes. From metagenomic analysis, genes encoding ferritin, flavodoxin, Fe transporters and siderophore uptake genes were more abundant in low-Fe waters, mirroring paradigms of low-Fe adaptation in diatoms. Distinct Fe-related gene repertories of HNLC ecotypes CRD1 and CRD2 also highlight how coexisting ecotypes have evolved independent approaches to life in low-Fe habitats. Synechococcus and Prochlorococcus HNLC ecotypes likewise exhibit independent, genome-wide reductions of predicted Fe-requiring genes. HNLC ecotype CRD1 interestingly was most similar to coastal ecotype I in Fe physiology and Fe-related gene content, suggesting populations from these different biomes experience similar Fe-selective conditions. This work supports an improved perspective that phytoplankton are shaped by more nuanced Fe niches in the oceans than previously implied from mostly binary comparisons of low- versus high-Fe habitats and populations.

Different iron storage strategies among bloom-forming diatoms

Diatoms are prominent eukaryotic phytoplankton despite being limited by the micronutrient iron in vast expanses of the ocean. As iron inputs are often sporadic, diatoms have evolved mechanisms such as the ability to store iron that enable them to bloom when iron is resupplied and then persist when low iron levels are reinstated. Two iron storage mechanisms have been previously described: the protein ferritin and vacuolar storage. To investigate the ecological role of these mechanisms among diatoms, iron addition and removal incubations were conducted using natural phytoplankton communities from varying iron environments. We show that among the predominant diatoms, Pseudo-nitzschia were favored by iron removal and displayed unique ferritin expression consistent with a long-term storage function. Meanwhile, Chaetoceros and Thalassiosira gene expression aligned with vacuolar storage mechanisms. Pseudo-nitzschia also showed exceptionally high iron storage under steady-state high and low iron conditions, as well as following iron resupply to iron-limited cells. We propose that bloom-forming diatoms use different iron storage mechanisms and that ferritin utilization may provide an advantage in areas of prolonged iron limitation with pulsed iron inputs. As iron distributions and availability change, this speculated ferritin-linked advantage may result in shifts in diatom community composition that can alter marine ecosystems and biogeochemical cycles.

Pervasive iron limitation at subsurface chlorophyll maxima of the California Current

The vertical distribution of phytoplankton cells and chlorophyll concentrations throughout the sunlit water column is rarely uniform. In many ocean regions, chlorophyll concentrations peak in distinct and persistent layers deep below the surface called subsurface chlorophyll maximum layers (SCMLs). SCML formation is hypothesized to reflect the consequences of phytoplankton light/macronutrient colimitation, behavior, and/or photoacclimation. We discovered unexpectedly persistent and widespread phytoplankton iron limitation and iron/light colimitation in SCMLs of the California Current and at the edge of the North Pacific Subtropical Gyre using shipboard incubations, metatranscriptomics, and biogeochemical proxies. These results suggest that interactions and feedbacks between iron and light availability play an important and previously unrecognized role in controlling the productivity and biogeochemical dynamics of SCMLs.

Subsurface chlorophyll maximum layers (SCMLs) are nearly ubiquitous in stratified water columns and exist at horizontal scales ranging from the submesoscale to the extent of oligotrophic gyres. These layers of heightened chlorophyll and/or phytoplankton concentrations are generally thought to be a consequence of a balance between light energy from above and a limiting nutrient flux from below, typically nitrate (NO₃). Here we present multiple lines of evidence demonstrating that iron (Fe) limits or with light colimits phytoplankton communities in SCMLs along a primary productivity gradient from coastal to oligotrophic offshore waters in the southern California Current ecosystem. SCML phytoplankton responded markedly to added Fe or Fe/light in experimental incubations and transcripts of diatom and picoeukaryote Fe stress genes were strikingly abundant in SCML metatranscriptomes. Using a biogeochemical proxy with data from a 40-y time series, we find that diatoms growing in California Current SCMLs are persistently Fe deficient during the spring and summer growing season. We also find that the spatial extent of Fe deficiency within California Current SCMLs has significantly increased over the last 25 y in line with a regional climate index. Finally, we show that diatom Fe deficiency may be common in the subsurface of major upwelling zones worldwide. Our results have important implications for our understanding of the biogeochemical consequences of marine SCML formation and maintenance.

Pseudo-nitzschia, Nitzschia, and domoic acid: New research since 2011

Some diatoms of the genera Pseudo-nitzschia and Nitzschia produce the neurotoxin domoic acid (DA), a compound that caused amnesic shellfish poisoning (ASP) in humans just over 30 years ago (December 1987) in eastern Canada. This review covers new information since two previous reviews in 2012. Nitzschia bizertensis was subsequently discovered to be toxigenic in Tunisian waters. The known distribution of N. navis-varingica has expanded from Vietnam to Malaysia, Indonesia, the Philippines and Australia. Furthermore, 15 new species (and one new variety) of Pseudo-nitzschia have been discovered, bringing the total to 52. Seven new species were found to produce DA, bringing the total of toxigenic species to 26. We list all Pseudo-nitzschia species, their ability to produce DA, and show their global distribution. A consequence of the extended distribution and increased number of toxigenic species worldwide is that DA is now found more pervasively in the food web, contaminating new marine organisms (especially marine mammals), affecting their physiology and disrupting ecosystems. Recent findings highlight how zooplankton grazers can induce DA production in Pseudo-nitzschia and how bacteria interact with Pseudo-nitzschia. Since 2012, new discoveries have been reported on physiological controls of Pseudo-nitzschia growth and DA production, its sexual reproduction, and infection by an oomycete parasitoid. Many advances are the result of applying molecular approaches to discovering new species, and to understanding the population genetic structure of Pseudo-nitzschia and mechanisms used to cope with iron limitation. The availability of genomes from three Pseudo-nitzschia species, coupled with a comparative transcriptomic approach, has allowed advances in our understanding of the sexual reproduction of Pseudo-nitzschia, its signaling pathways, its interactions with bacteria, and genes involved in iron and vitamin B₁₂ and B₇ metabolism. Although there have been no new confirmed cases of ASP since 1987 because of monitoring efforts, new blooms have occurred. A massive toxic Pseudo-nitzschia bloom affected the entire west coast of North America during 2015-2016, and was linked to a 'warm blob' of ocean water. Other smaller toxic blooms occurred in the Gulf of Mexico and east coast of North America. Knowledge gaps remain, including how and why DA and its isomers are produced, the world distribution of potentially toxigenic Nitzschia species, the prevalence of DA isomers, and molecular markers to discriminate between toxigenic and non-toxigenic species and to discover sexually reproducing populations in the field.

Endocytosis-mediated siderophore uptake as a strategy for Fe acquisition in diatoms

Phytoplankton growth is limited in vast oceanic regions by the low bioavailability of iron. Iron fertilization often results in diatom blooms, yet the physiological underpinnings for how diatoms survive in chronically iron-limited waters and outcompete other phytoplankton when iron becomes available are unresolved. We show that some diatoms can use siderophore-bound iron, and exhibit a species-specific recognition for siderophore types. In Phaeodactylum tricornutum, hydroxamate siderophores are taken up without previous reduction by a high-affinity mechanism that involves binding to the cell surface followed by endocytosis-mediated uptake and delivery to the chloroplast. The affinity recorded is the highest ever described for an iron transport system in any eukaryotic cell. Collectively, our observations suggest that there are likely a variety of iron uptake mechanisms in diatoms besides the well-established reductive mechanism. We show that iron starvation-induced protein 1 (ISIP1) plays an important role in the uptake of siderophores, and through bioinformatics analyses we deduce that this protein is largely diatom-specific. We quantify expression of ISIP1 in the global ocean by querying the Tara Oceans atlas of eukaryotic genes and show a link between the abundance and distribution of diatom-associated ISIP1 with ocean provinces defined by chronic iron starvation.