Hong Zheng

Domoic acid (DA) production by diatom genus Pseudo-nitzschia causes recurrent harmful algal blooms with major consequences for marine ecosystems, fisheries, and human health. Predicting DA risk remains challenging because toxin production varies across morphologically similar and cryptic lineages. Consequently, morphology or single marker genes often fail to resolve toxigenic potential. Here, we expand genomic resources for Pseudo-nitzschia by generating and compiling high-quality genomes spanning diverse species and biogeographic origins, establishing a genus-wide phylogenomic framework to investigate the distribution and evolution of the DA biosynthesis ( dab ) gene cluster. Among chemically characterized strains, the dab pathway was detected in species only with confirmed DA production. While the cluster exhibits a conserved core structure, several lineages display signatures of ongoing evolution, including intronization, locus rearrangement, and pseudogenization. One striking example of dab variation occurs in toxigenic P. australis strains, in which its biochemically validated isodomoic acid A synthase-encoding gene dabC is located outside the dab locus. Metatranscriptomic surveys reveal widespread environmental dab expression across size classes and species, implicating taxa not typically associated with environmental DA events, such as P . cf. obtusa and P . cf. plurisecta . To connect dab expression with chemical output, nutrient limitation experiments in the historically inconsistent DA producer P. pseudodelicatissima combined transcriptional profiling and nontargeted metabolomics demonstrating coordinated dab expression across distinct co-expression modules and DA accumulation in response to silicon limitation. Together, these results refine the phylogenetic scope of DA production and provide a genome-resolved basis for improved molecular monitoring of toxigenic diversity.

The McMurdo Dry Valleys contain a mosaic of specialized microbial habitats structured by strong physical and chemical gradients. Blood Falls, located at the terminus of the Taylor Glacier, is a red, iron-rich outflow of deep Antarctic subglacial brine with geochemical and isotopic evidence supporting a potential ancient marine origin. Previous studies suggest the brine formed when seawater inundated the Taylor Valley over warm climatic intervals before becoming isolated beneath the advancing glacier. Here, across 167 aquatic, sediment and aeolian samples from the McMurdo Dry Valleys and marine reference sites, we identified a distinct marine micro-eukaryotic assemblage restricted to red-hued ice, mud and sediment at the Taylor Glacier terminus. Marine indicator species were identified across several phylogenetic groups including diatoms, haptophytes, dinoflagellates and ciliates. Metatranscriptomic profiles revealed transcriptionally active phototrophs with enriched pathways for photosynthesis, osmotic stress responses and cellular repair. Haplotype networks showed lineage-specific divergence between Taylor Glacier terminus and McMurdo Sound diatoms, supporting geographic isolation. These findings indicate that the subglacial brine-fed system at the Taylor Glacier terminus retains marine-derived biological signatures long after physical separation from the ocean, linking contemporary Antarctic microbial assemblages to past climatic transitions.

Global overturning circulation partitions the deep ocean into regions, each with different physicochemical characteristics, but the extent to which these water masses represent distinct ecosystems remains unknown. In this work, we integrate extensive genomic information with hydrography and water mass age to delineate microbial taxonomic and functional boundaries across the South Pacific. Prokaryotic richness steeply increases with depth in the surface ocean, which forms a so-called phylocline, below which, richness is consistently high, dipping slightly in highly aged water. Reconstructed genomes self-organize into six spatially distinct taxonomic cohorts and 10 functionally distinct biomes that are primarily structured by wind-driven circulation at the surface and density-driven circulation at depth. Overall, water physicochemistry, modulated at depth by water age, drives microbial diversity patterns and functional potential in the pelagic ocean.

Global overturning circulation partitions the deep ocean into regions with unique physicochemical characteristics, but the extent to which these water masses represent distinct ecosystems remains unknown. Here, we integrate extensive genomic information with hydrography and water mass age to delineate microbial taxonomic and functional boundaries across the South Pacific. Prokaryotic richness steeply increases with depth in the surface ocean, forming a "phylocline", below which richness is consistently high, dipping slightly in highly aged water. Reconstructed genomes self-organize into six spatially-distinct taxonomic cohorts and ten functionally-distinct biomes that are primarily structured by wind-driven circulation at the surface and density-driven circulation at depth. Overall, water physicochemistry, modulated at depth by water age, drives microbial diversity and functional potential in the pelagic ocean.

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.

Ocean microbes are the foundation of marine food webs, regulating carbon cycling and ecosystem dynamics. How they proliferate, die, move, and interact is regulated by physical, chemical, and biological factors that are dynamic and challenging to quantify in the natural environment. A significant limitation in many marine field studies is the inability to continuously sample the ever-changing ocean environment over space and time. In this study, we integrated spatiotemporal and multi-omic sample collection in an intensive sampling effort of phytoplankton ecology in Monterey Bay, California during the spring of 2021. Sampling methods coupled: (1) manual shipboard CTD sampling, (2) autonomous sampling using a Long-Range Autonomous Underwater Vehicle (LRAUV) equipped with an Environmental Sampling Processor (ESP), and (3) high-resolution physical measurements by an autonomous vertical profiler (Wirewalker). Sampling occurred as upwelling waned alongside declining domoic acid (DA) and low abundances of toxigenic Pseudo-nitzschia . Conditions needed to spark a widespread and toxic Pseudo-nitzschia bloom were absent, yet low-level DA was driven by similar mechanisms to those causing elevated DA. Three DA biosynthetic intermediate molecules were reported in the environment for the first time. Both shipboard and ESP sampling approaches identified DA biosynthetic gene expression at frontal zones. DA and expression of dabA , the gene encoding the first committed step of DA biosynthesis, were higher in association with recently upwelled water that supplied nutrients for growth and DA biosynthesis. Detection of subtle variations in dab gene expression in response to environmental variation provide a window into the ecological dynamics underpinning major toxic events. Graphical

Genetic inactivation of diatom nitrate reductase abolishes nitrate assimilation but not nitrate uptake and results in transcriptional activation of nitrate stor