Molecular Forecasting of Domoic Acid during a Pervasive Toxic Diatom Bloom

Abstract 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.