Search

Time-series observations are critical to understand the structure, function, and dynamics of marine ecosystems. The Hawaii Ocean Time-series program has maintained near-monthly sampling at Station ALOHA (22°45′N, 158°00′W) in the oligotrophic North Pacific Subtropical Gyre (NPSG) since 1988 and has identified ecosystem variability over seasonal to interannual timescales. To further extend the temporal resolution of these near-monthly time-series observations, an extensive field campaign was conducted during July-September 2012 at Station ALOHA with near-daily sampling of upper water-column biogeochemistry, phytoplankton abundance, and activity. The resulting data set provided biogeochemical measurements at high temporal resolution and documents two important events at Station ALOHA: (1) a prolonged period of low productivity when net community production in the mixed layer shifted to a net heterotrophic state and (2) detection of a distinct sea-surface salinity minimum feature which was prominent in the upper water column (0-50 m) for a period of approximately 30 days. The shipboard observations during July-September 2012 were supplemented with in situ measurements provided by Seagliders, profiling floats, and remote satellite observations that together revealed the extent of the low productivity and the sea-surface salinity minimum feature in the NPSG.

Importance: Colonies of the cyanobacteria Trichodesmium act as a biological hotspot for the usage and recycling of key resources such as C, N, P, and Fe within an otherwise oligotrophic environment. While Trichodesmium colonies are known to interact and support a unique community of algae and particle-associated microbes, our understanding of the taxa that populate these colonies and the gene functions they encode is still limited. Characterizing the taxa and adaptive strategies that influence consortium physiology and its concomitant biogeochemistry is critical in a future ocean predicted to have increasingly resource-depleted regions., Competing Interests: The authors declare no conflict of interest.

Diatoms are important components of the marine food web and one of the most species-rich groups of phytoplankton. The diversity and composition of diatoms in eutrophic nearshore habitats have been well documented due to the outsized influence of diatoms on coastal ecosystem functioning. In contrast, patterns of both diatom diversity and community composition in offshore oligotrophic regions where diatom biomass is low have been poorly resolved. To compare the diatom diversity and community composition in oligotrophic and eutrophic waters, diatom communities were sampled along a 1,250 km transect from the oligotrophic Sargasso Sea to the coastal waters of the northeast US shelf. Diatom community composition was determined by amplifying and sequencing the 18S rDNA V4 region. Of the 301 amplicon sequence variants (ASVs) identified along the transect, the majority (70%) were sampled exclusively from oligotrophic waters of the Gulf Stream and Sargasso Sea and included the genera Bacteriastrum, Haslea, Hemiaulus, Pseudo-nitzschia, and Nitzschia. Diatom ASV richness did not vary along the transect, indicating that the oligotrophic Sargasso Sea and Gulf Stream are occupied by a diverse diatom community. Although ASV richness was similar between oligotrophic and coastal waters, diatom community composition in these regions differed significantly and was correlated with temperature and phosphate, two environmental variables known to influence diatom metabolism and geographic distribution. In sum, oligotrophic waters of the western North Atlantic harbor diverse diatom assemblages that are distinct from coastal regions, and these open ocean diatoms warrant additional study, as they may play critical roles in oligotrophic ecosystems., (© 2023 Phycological Society of America.)

One-quarter of photosynthesis-derived carbon on Earth rapidly cycles through a set of short-lived seawater metabolites that are generated from the activities of marine phytoplankton, bacteria, grazers and viruses. Here we discuss the sources of microbial metabolites in the surface ocean, their roles in ecology and biogeochemistry, and approaches that can be used to analyse them from chemistry, biology, modelling and data science. Although microbial-derived metabolites account for only a minor fraction of the total reservoir of marine dissolved organic carbon, their flux and fate underpins the central role of the ocean in sustaining life on Earth., (© 2022. Springer Nature Limited.)

The pelagophyte Aureococcus anophagefferens causes harmful brown tide blooms in marine embayments on three continents. Aureococcus anophagefferens was the first harmful algal bloom species to have its genome sequenced, an advance that evidenced genes important for adaptation to environmental conditions that prevail during brown tides. To expand the genomic tools available for this species, genomes for four strains were assembled, including three newly sequenced strains and one assembled from publicly available data. These genomes ranged from 57.11 to 73.62 Mb, encoding 13,191-17,404 potential proteins. All strains shared ~90% of their encoded proteins as determined by homology searches and shared most functional orthologs as determined by KEGG, although each strain also possessed coding sequences with unique functions. Like the original reference genome, the genomes assembled in this study possessed genes hypothesized to be important in bloom proliferation, including genes involved in organic compound metabolism and growth at low light. Cross-strain informatics and culture experiments suggest that the utilization of purines is a potentially important source of organic nitrogen for brown tides. Analyses of metatranscriptomes from a brown tide event demonstrated that use of a single genome yielded a lower read mapping percentage (~30% of library reads) as compared to a database generated from all available genomes (~43%), suggesting novel information about bloom ecology can be gained from expanding genomic space. This work demonstrates the continued need to sequence ecologically relevant algae to understand the genomic potential and their ecology in the environment., (© 2021 Phycological Society of America.)

Complex assemblages of microbes in the surface ocean are responsible for approximately half of global carbon fixation. The persistence of high taxonomic diversity despite competition for a small suite of relatively homogeneously distributed nutrients, that is, 'the paradox of the plankton', represents a long-standing challenge for ecological theory. Here we find evidence consistent with temporal niche partitioning of nitrogen assimilation processes over a diel cycle in the North Pacific Subtropical Gyre. We jointly analysed transcript abundances, lipids and metabolites and discovered that a small number of diel archetypes can explain pervasive periodic dynamics. Metabolic pathway analysis of identified diel signals revealed asynchronous timing in the transcription of nitrogen uptake and assimilation genes among different microbial groups-cyanobacteria, heterotrophic bacteria and eukaryotes. This temporal niche partitioning of nitrogen uptake emerged despite synchronous transcription of photosynthesis and central carbon metabolism genes and associated macromolecular abundances. Temporal niche partitioning may be a mechanism by which microorganisms in the open ocean mitigate competition for scarce resources, supporting community coexistence., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

The physical and biological dynamics that influence phytoplankton communities in the oligotrophic ocean are complex, changing across broad temporal and spatial scales. Eukaryotic phytoplankton (e.g., diatoms), despite their relatively low abundance in oligotrophic waters, are responsible for a large component of the organic matter flux to the ocean interior. Mesoscale eddies can impact both microbial community structure and function, enhancing primary production and carbon export, but the mechanisms that underpin these dynamics are still poorly understood. Here, mesoscale eddy influences on the taxonomic diversity and expressed functional profiles of surface communities of microeukaryotes and particle-associated heterotrophic bacteria from the North Pacific Subtropical Gyre were assessed over 2 years (spring 2016 and summer 2017). The taxonomic diversity of the microeukaryotes significantly differed by eddy polarity (cyclonic versus anticyclonic) and between sampling seasons/years and was significantly correlated with the taxonomic diversity of particle-associated heterotrophic bacteria. The expressed functional profile of these taxonomically distinct microeukaryotes varied consistently as a function of eddy polarity, with cyclones having a different expression pattern than anticyclones, and between sampling seasons/years. These data suggest that mesoscale forcing, and associated changes in biogeochemistry, could drive specific physiological responses in the resident microeukaryote community, independent of species composition., (© 2021 Society for Applied Microbiology and John Wiley & Sons Ltd.)

There is growing interest in the use of metatranscriptomics to study virus community dynamics. We used RNA samples collected from harmful brown tides caused by the eukaryotic alga Aureococcus anophagefferens within New York (United States) estuaries and in the process observed how preprocessing of libraries by either selection for polyadenylation or reduction in ribosomal RNA (rRNA) influenced virus community analyses. As expected, more reads mapped to the A. anophagefferens genome in polyadenylation-selected libraries compared to the rRNA-reduced libraries, with reads mapped in each sample correlating to one another regardless of preprocessing of libraries. Yet, this trend was not seen for reads mapping to the Aureococcus anophagefferens Virus (AaV), where significantly more reads (approximately two orders of magnitude) were mapped to the AaV genome in the rRNA-reduced libraries. In the rRNA-reduced libraries, there was a strong and significant correlation between reads mappings to AaV and A. anophagefferens . Overall, polyadenylation-selected libraries produced fewer viral contigs, fewer reads mapped to viral contigs, and different proportions across viral realms and families, compared to their rRNA-reduced pairs. This study provides evidence that libraries generated by rRNA reduction and not selected for polyadenylation are more appropriate for quantitative characterization of viral communities in aquatic ecosystems by metatranscriptomics., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Gann, Kang, Dyhrman, Gobler and Wilhelm.)

Sunlight is the most important environmental control on diel fluctuations in phytoplankton activity, and understanding diel microbial processes is essential to the study of oceanic biogeochemical cycles. Yet, little is known about the in situ temporal dynamics of phytoplankton metabolic activities and their coordination across different populations. We investigated diel orchestration of phytoplankton activity in photosynthesis, photoacclimation, and photoprotection by analyzing pigment and quinone distributions in combination with metatranscriptomes in surface waters of the North Pacific Subtropical Gyre (NPSG). We found diel cycles in pigment abundances resulting from the balance of their synthesis and consumption. These dynamics suggest that night represents a metabolic recovery phase, refilling cellular pigment stores, while photosystems are remodeled towards photoprotection during daytime. Transcript levels of genes involved in photosynthesis and pigment metabolism had synchronized diel expression patterns among all taxa, reflecting the driving force light imparts upon photosynthetic organisms in the ocean, while other environmental factors drive niche differentiation. For instance, observed decoupling of diel oscillations in transcripts and related pigments indicates that pigment abundances are modulated by environmental factors extending beyond gene expression/regulation reinforcing the need to combine metatranscriptomics with proteomics and metabolomics to fully understand the timing of these critical processes in situ.

The widespread coccolithophore Emiliania huxleyi is an abundant oceanic phytoplankton, impacting the global cycling of carbon through both photosynthesis and calcification. Here, we examined the transcriptional responses of populations of E. huxleyi in the North Pacific Subtropical Gyre to shifts in the nutrient environment. Using a metatranscriptomic approach, nutrient-amended microcosm studies were used to track the global metabolism of E. huxleyi. The addition of nitrate led to significant changes in transcript abundance for gene pathways involved in nitrogen and phosphorus metabolism, with a decrease in the abundance of genes involved in the acquisition of nitrogen (e.g. N-transporters) and an increase in the abundance of genes associated with phosphate acquisition (e.g. phosphatases). Simultaneously, after the addition of nitrate, genes associated with calcification and genes unique to the diploid life stages of E. huxleyi significantly increased. These results suggest that nitrogen is a major driver of the physiological ecology of E. huxleyi in this system and further suggest that the addition of nitrate drives shifts in the dominant life-stage of the population. Together, these results underscore the importance of phenotypic plasticity to the success of E. huxleyi, a characteristic that likely underpins its ability to thrive across a variety of marine environments., (© 2020 Society for Applied Microbiology and John Wiley & Sons Ltd.)

Emiliania huxleyi is a calcifying haptophyte, contributing to both the organic and inorganic marine carbon cycles. In marine ecosystems, light is a major driver of phytoplankton physiology and ultimately carbon flow through the ecosystem. Here, we analysed a Lagrangian time-series of metatranscriptomes collected in the North Pacific Subtropical Gyre (NPSG) to examine how in situ populations of E. huxleyi modulate gene expression over day-night transitions. Many E. huxleyi contigs had a diel expression pattern, with 61% of contigs clustering into modules with statistically significant diel periodicity. Contigs involved in processes that build up energy stores, like carbon fixation and lipid synthesis, peaked around dawn. In contrast, contigs involved in processes that released energy stores, like respiration and lipid degradation, peaked mid-day and towards dusk. These patterns suggest an orchestrated cycle of building, then consuming energy stores in E. huxleyi populations in the NPSG. Selected contigs related to the cell cycle also exhibited significant diel periodicity consistent with phased modulations of division observed in culture. Overall, these patterns of gene expression suggest a daily metabolic cascade that could contribute to both organic and inorganic carbon flow in this nutrient depleted ecosystem., (© 2019 Society for Applied Microbiology and John Wiley & Sons Ltd.)

Climate change is predicted to increase the severity and prevalence of harmful algal blooms (HABs). In the past twenty years, omics techniques such as genomics, transcriptomics, proteomics and metabolomics have transformed that data landscape of many fields including the study of HABs. Advances in technology have facilitated the creation of many publicly available omics datasets that are complementary and shed new light on the mechanisms of HAB formation and toxin production. Genomics have been used to reveal differences in toxicity and nutritional requirements, while transcriptomics and proteomics have been used to explore HAB species responses to environmental stressors, and metabolomics can reveal mechanisms of allelopathy and toxicity. In this review, we explore how omics data may be leveraged to improve predictions of how climate change will impact HAB dynamics. We also highlight important gaps in our knowledge of HAB prediction, which include swimming behaviors, microbial interactions and evolution that can be addressed by future studies with omics tools. Lastly, we discuss approaches to incorporate current omics datasets into predictive numerical models that may enhance HAB prediction in a changing world. With the ever-increasing omics databases, leveraging these data for understanding climate-driven HAB dynamics will be increasingly powerful., (Copyright © 2019 Elsevier B.V. All rights reserved.)

From June to August 2018, the eruption of Kīlauea volcano on the island of Hawai'i injected millions of cubic meters of molten lava into the nutrient-poor waters of the North Pacific Subtropical Gyre. The lava-impacted seawater was characterized by high concentrations of metals and nutrients that stimulated phytoplankton growth, resulting in an extensive plume of chlorophyll a that was detectable by satellite. Chemical and molecular evidence revealed that this biological response hinged on unexpectedly high concentrations of nitrate, despite the negligible quantities of nitrogen in basaltic lava. We hypothesize that the high nitrate was caused by buoyant plumes of nutrient-rich deep waters created by the substantial input of lava into the ocean. This large-scale ocean fertilization was therefore a unique perturbation event that revealed how marine ecosystems respond to exogenous inputs of nutrients., (Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

The N 2 fixing cyanobacterium Trichodesmium is a critically important organism in oligotrophic marine ecosystems, supplying "new" nitrogen (N) to the otherwise N-poor tropical and subtropical regions where it occurs. Low concentrations of phosphorus (P) in these regions can constrain Trichodesmium distribution and N 2 fixation rates. Physiological characterization of a single species in a mixed community can be challenging, and 'omic approaches are increasingly important tools for tracking nutritional physiology in a taxon-specific manner. As such, studies examining the dynamics of gene and protein markers of physiology (e.g., nutrient stress) are critical for the application and interpretation of such 'omic data in situ . Here we leveraged combined transcriptomics, proteomics, and enzyme activity assays to track the physiological response of Trichodesmium erythraeum IMS101 to P deficiency and subsequent P re-supply over 72 h of sampling. P deficiency resulted in differential gene expression, protein abundance, and enzyme activity that highlighted a synchronous shift in P physiology with increases in the transcripts and corresponding proteins for hydrolyzing organic phosphorus, taking up phosphate with higher affinity, and modulating intracellular P demand. After P deficiency was alleviated, gene expression of these biomarkers was reduced to replete levels within 4 h of P amendment. A number of these gene biomarkers were adjacent to putative pho boxes and their expression patterns were similar to a sphR response regulator. Protein products of the P deficiency biomarkers were slow to decline, with 84% of the original P deficient protein set still significantly differentially expressed after 72 h. Alkaline phosphatase activity tracked with proteins for this enzyme. With the rapid turnover time of transcripts, they appear to be good biomarkers of a P stress phenotype, whereas proteins, with a slower turnover time, may better reflect cellular activities. These results highlight the importance of validating and pairing transcriptome and proteome data that can be applied to physiological studies of key species in situ .

Harmful algal blooms (HABs) threaten ecosystems and human health worldwide. Controlling nitrogen inputs to coastal waters is a common HAB management strategy, as nutrient concentrations often suggest coastal blooms are nitrogen-limited. However, defining best nutrient management practices is a long-standing challenge: in part, because of difficulties in directly tracking the nutritional physiology of harmful species in mixed communities. Using metatranscriptome sequencing and incubation experiments, we addressed this challenge by assaying the in situ physiological ecology of the ecosystem destructive alga, Aureococcus anophagefferens . Here we show that gene markers of phosphorus deficiency were expressed in situ , and modulated by the enrichment of phosphorus, which was consistent with the observed growth rate responses. These data demonstrate the importance of phosphorus in controlling brown-tide dynamics, suggesting that phosphorus, in addition to nitrogen, should be evaluated in the management and mitigation of these blooms. Given that nutrient concentrations alone were suggestive of a nitrogen-limited ecosystem, this study underscores the value of directly assaying harmful algae in situ for the development of management strategies.

Heterosigma akashiwo is a raphidophyte known for forming ichthyotoxic blooms. In order to predict the potential impacts of rising CO 2 on H. akashiwo it is necessary to understand the factors influencing growth rates over a range of CO 2 concentrations. Here we examined the physiology and gene expression response of H. akashiwo to concentrations from 200 to 1000ppm CO 2 . Growth rate data were combined from this and previous studies and fit with a CO 2 limitation-inhibition model that revealed an apparent growth optimum around 600-800ppm CO 2 . Physiological changes included a significant increase in C:N ratio at ∼800ppm CO 2 and a significant decrease in hydrogen peroxide concentration at ∼1000ppm. Whole transcriptome sequencing of H. akashiwo revealed sharp distinctions in metabolic pathway gene expression between ∼600 and ∼800ppm CO 2 . Hierarchical clustering by co-expression identified groups of genes with significant correlations to CO 2 and growth rate. Genes with significant differential expression with CO 2 included carbon concentrating mechanism genes such as beta-carbonic anhydrases and a bicarbonate transporter, which may underpin shifts in physiology. Genes involved in cell motility were significantly changed by both elevated CO 2 and growth rate, suggesting that future ocean conditions could modify swimming behavior in this species., (Copyright © 2018 Elsevier GmbH. All rights reserved.)

In the surface ocean, light fuels photosynthetic carbon fixation of phytoplankton, playing a critical role in ecosystem processes including carbon export to the deep sea. In oligotrophic oceans, diatom-diazotroph associations (DDAs) play a keystone role in ecosystem function because diazotrophs can provide otherwise scarce biologically available nitrogen to the diatom host, fueling growth and subsequent carbon sequestration. Despite their importance, relatively little is known about the nature of these associations in situ. Here we used metatranscriptomic sequencing of surface samples from the North Pacific Subtropical Gyre (NPSG) to reconstruct patterns of gene expression for the diazotrophic symbiont Richelia and we examined how these patterns were integrated with those of the diatom host over day-night transitions. Richelia exhibited significant diel signals for genes related to photosynthesis, N 2 fixation, and resource acquisition, among other processes. N 2 fixation genes were significantly co-expressed with host nitrogen uptake and metabolism, as well as potential genes involved in carbon transport, which may underpin the exchange of nitrogen and carbon within this association. Patterns of expression suggested cell division was integrated between the host and symbiont across the diel cycle. Collectively these data suggest that symbiont-host physiological ecology is strongly interconnected in the NPSG.

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., (© 2018 Society for Applied Microbiology and John Wiley & Sons Ltd.)

The N 2 -fixing cyanobacterium Trichodesmium is intensely studied because of the control this organism exerts over the cycling of carbon and nitrogen in the low nutrient ocean gyres. Although iron (Fe) and phosphorus (P) bioavailability are thought to be major drivers of Trichodesmium distributions and activities, identifying resource controls on Trichodesmium is challenging, as Fe and P are often organically complexed and their bioavailability to a single species in a mixed community is difficult to constrain. Further, Fe and P geochemistries are linked through the activities of metalloenzymes, such as the alkaline phosphatases (APs) PhoX and PhoA, which are used by microbes to access dissolved organic P (DOP). Here we identified significant correlations between Trichodesmium-specific transcriptional patterns in the North Atlantic (NASG) and North Pacific Subtropical Gyres (NPSG) and patterns in Fe and P biogeochemistry, with the relative enrichment of Fe stress markers in the NPSG, and P stress markers in the NASG. We also observed the differential enrichment of Fe-requiring PhoX transcripts in the NASG and Fe-insensitive PhoA transcripts in the NPSG, suggesting that metalloenzyme switching may be used to mitigate Fe limitation of DOP metabolism in Trichodesmium. This trait may underpin Trichodesmium success across disparate ecosystems.

Trichodesmium is a widespread, N 2 fixing marine cyanobacterium that drives inputs of newly fixed nitrogen and carbon into the oligotrophic ecosystems where it occurs. Colonies of Trichodesmium ubiquitously occur with heterotrophic bacteria that make up a diverse microbiome, and interactions within this Trichodesmium holobiont could influence the fate of fixed carbon and nitrogen. Metatranscriptome sequencing was performed on Trichodesmium colonies collected during high-frequency Lagrangian sampling in the North Pacific Subtropical Gyre (NPSG) to identify possible interactions between the Trichodesmium host and microbiome over day-night cycles. Here we show significantly coordinated patterns of gene expression between host and microbiome, many of which had significant day-night periodicity. The functions of the co-expressed genes suggested a suite of interactions within the holobiont linked to key resources including nitrogen, carbon, and iron. Evidence of microbiome reliance on Trichodesmium-derived vitamin B12 was also detected in co-expression patterns, highlighting a dependency that could shape holobiont community structure. Collectively, these patterns of expression suggest that biotic interactions could influence colony cycling of resources like nitrogen and vitamin B12, and decouple activities, like N 2 fixation, from typical abiotic drivers of Trichodesmium physiological ecology.

Prochlorococcus is a globally important marine cyanobacterium that lacks the gene catalase and relies on 'helper' bacteria such as Alteromonas to remove reactive oxygen species. Increasing atmospheric CO 2 decreases the need for carbon concentrating mechanisms and photorespiration in phytoplankton, potentially altering their metabolism and microbial interactions even when carbon is not limiting growth. Here, Prochlorococcus (VOL4, MIT9312) was co-cultured with Alteromonas (strain EZ55) under ambient (400 p.p.m.) and elevated CO 2 (800 p.p.m.). Under elevated CO 2 , Prochlorococcus had a significantly longer lag phase and greater apparent die-offs after transfers suggesting an increase in oxidative stress. Whole-transcriptome analysis of Prochlorococcus revealed decreased expression of the carbon fixation operon, including carboxysome subunits, corresponding with significantly fewer carboxysome structures observed by electron microscopy. Prochlorococcus co-culture responsive gene 1 had significantly increased expression in elevated CO 2 , potentially indicating a shift in the microbial interaction. Transcriptome analysis of Alteromonas in co-culture with Prochlorococcus revealed decreased expression of the catalase gene, known to be critical in relieving oxidative stress in Prochlorococcus by removing hydrogen peroxide. The decrease in catalase gene expression was corroborated by a significant ~6-fold decrease in removal rates of hydrogen peroxide from co-cultures. These data suggest Prochlorococcus may be more vulnerable to oxidative stress under elevated CO 2 in part from a decrease in ecosystem services provided by heterotrophs like Alteromonas. This work highlights the importance of considering microbial interactions in the context of a changing ocean., (© 2018. International Society for Microbial Ecology.)

With rising atmospheric CO 2 , phytoplankton face shifts in ocean chemistry including increased dissolved CO 2 and acidification that will likely influence the relative competitive fitness of different phytoplankton taxa. Here we compared the physiological and gene expression responses of six species of phytoplankton including a diatom, a raphidophyte, two haptophytes, and two dinoflagellates to ambient (~400 ppm) and elevated (~800 ppm) CO 2 . Dinoflagellates had significantly slower growth rates and higher, yet variable, chlorophyll a per cell under elevated CO 2 . The other phytoplankton tended to have increased growth rates and/or decreased chlorophyll a per cell. Carbon and nitrogen partitioning of cells shifted under elevated CO 2 in some species, indicating potential changes in energy fluxes due to changes in carbon concentrating mechanisms (CCM) or photorespiration. Consistent with these phenotypic changes, gene set enrichment analyses revealed shifts in energy, carbon and nitrogen metabolic pathways, though with limited overlap between species in the genes and pathways involved. Similarly, gene expression responses across species revealed few conserved CO 2 -responsive genes within CCM and photorespiration categories, and a survey of available transcriptomes found high diversity in biophysical CCM and photorespiration expressed gene complements between and within the four phyla represented by these species. The few genes that displayed similar responses to CO 2 across phyla were from understudied gene families, making them targets for further research to uncover the mechanisms of phytoplankton acclimation to elevated CO 2 . These results underscore that eukaryotic phytoplankton have diverse gene complements and gene expression responses to CO 2 perturbations and highlight the value of cross-phyla comparisons for identifying gene families that respond to environmental change.

Trichodesmium is a genus of marine diazotrophic colonial cyanobacteria that exerts a profound influence on global biogeochemistry, by injecting 'new' nitrogen into the low nutrient systems where it occurs. Colonies of Trichodesmium ubiquitously contain a diverse assemblage of epibiotic microorganisms, constituting a microbiome on the Trichodesmium host. Metagenome sequences from Trichodesmium colonies were analyzed along a resource gradient in the western North Atlantic to examine microbiome community structure, functional diversity and metabolic contributions to the holobiont. Here we demonstrate the presence of a core Trichodesmium microbiome that is modulated to suit different ocean regions, and contributes over 10 times the metabolic potential of Trichodesmium to the holobiont. Given the ubiquitous nature of epibionts on colonies, the substantial functional diversity within the microbiome is likely an integral facet of Trichodesmium physiological ecology across the oligotrophic oceans where this biogeochemically significant diazotroph thrives.

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., (Copyright © 2017 Elsevier B.V. All rights reserved.)

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.

Establishing virus-host relationships has historically relied on culture-dependent approaches. Here we report on the use of marine metatranscriptomics to probe virus-host relationships. Statistical co-occurrence analyses of dsDNA, ssRNA and dsRNA viral markers of polyadenylation-selected RNA sequences from microbial communities dominated by Aureococcus anophagefferens (Quantuck Bay, NY), and diatoms (Narragansett Bay, RI) show active infections by diverse giant viruses (NCLDVs) associated with algal and nonalgal hosts. Ongoing infections of A. anophagefferens by a known Mimiviridae (AaV) occur during bloom peak and decline. Bloom decline is also accompanied by increased activity of viruses other than AaV, including (+) ssRNA viruses. In Narragansett Bay, increased temporal resolution reveals active NCLDVs with both 'boom-and-bust' and 'steady-state infection'-like ecologies that include known as well as novel virus-host interactions. Our approach offers a method for screening active viral infections and develops links between viruses and their potential hosts in situ. Our observations further demonstrate that previously unknown virus-host relationships in marine systems are abundant.

Protists, which are single-celled eukaryotes, critically influence the ecology and chemistry of marine ecosystems, but genome-based studies of these organisms have lagged behind those of other microorganisms. However, recent transcriptomic studies of cultured species, complemented by meta-omics analyses of natural communities, have increased the amount of genetic information available for poorly represented branches on the tree of eukaryotic life. This information is providing insights into the adaptations and interactions between protists and other microorganisms and macroorganisms, but many of the genes sequenced show no similarity to sequences currently available in public databases. A better understanding of these newly discovered genes will lead to a deeper appreciation of the functional diversity and metabolic processes in the ocean. In this Review, we summarize recent developments in our understanding of the ecology, physiology and evolution of protists, derived from transcriptomic studies of cultured strains and natural communities, and discuss how these novel large-scale genetic datasets will be used in the future.

Ciliates from the genus Mesodinium are globally distributed in marine and freshwater ecosystems and may possess either heterotrophic or mixotrophic nutritional modes. Members of the Mesodinium major/rubrum species complex photosynthesize by sequestering and maintaining organelles from cryptophyte prey, and under certain conditions form periodic or recurrent blooms (= red tides). Here, we present an analysis of the genetic diversity of Mesodinium and cryptophyte populations from 10 environmental samples (eight globally dispersed habitats including five Mesodinium blooms), using group-specific primers for Mesodinium partial 18S, ITS, and partial 28S rRNA genes as well as cryptophyte large subunit RuBisCO genes ( rbcL ). In addition, 22 new cryptophyte and four new M. rubrum cultures were used to extract DNA and sequence rbcL and 18S-ITS-28S genes, respectively, in order to provide a stronger phylogenetic context for our environmental sequences. Bloom samples were analyzed from coastal Brazil, Chile, two Northeastern locations in the United States, and the Pribilof Islands within the Bering Sea. Additionally, samples were also analyzed from the Baltic and Barents Seas and coastal California under non-bloom conditions. Most blooms were dominated by a single Mesodinium genotype, with coastal Brazil and Chile blooms composed of M. major and the Eastern USA blooms dominated by M. rubrum variant B. Sequences from all four blooms were dominated by Teleaulax amphioxeia -like cryptophytes. Non-bloom communities revealed more diverse assemblages of Mesodinium spp., including heterotrophic species and the mixotrophic Mesodinium chamaeleon . Similarly, cryptophyte diversity was also higher in non-bloom samples. Our results confirm that Mesodinium blooms may be caused by M. major , as well as multiple variants of M. rubrum , and further implicate T. amphioxeia as the key cryptophyte species linked to these phenomena in temperate and subtropical regions.

Populations of nitrogen-fixing cyanobacteria in the genus Trichodesmium are critical to ocean ecosystems, yet predicting patterns of Trichodesmium distribution and their role in ocean biogeochemistry is an ongoing challenge. This may, in part, be due to differences in the physiological ecology of Trichodesmium species, which are not typically considered independently in field studies. In this study, the abundance of the two dominant Trichodesmium clades (Clade I and Clade III) was investigated during a survey at Station ALOHA in the North Pacific Subtropical Gyre (NPSG) using a clade-specific qPCR approach. While Clade I dominated the Trichodesmium community, Clade III abundance was >50% in some NPSG samples, in contrast to the western North Atlantic where Clade III abundance was always <10%. Clade I populations were distributed down to depths >80 m, while Clade III populations were only observed in the mixed layer and found to be significantly correlated with depth and temperature. These data suggest active niche partitioning of Trichodesmium species from different clades, as has been observed in other cyanobacteria. Tracking the distribution and physiology of Trichodesmium spp. would contribute to better predictions of the physiological ecology of this biogeochemically important genus in the present and future ocean., (© 2016 Society for Applied Microbiology and John Wiley & Sons Ltd.)

Nitrogen-fixing cyanobacteria in the genus Trichodesmium play a critical role in the productivity of the tropical and subtropical oligotrophic oceans. The ecological success of these populations is likely associated with the diverse microbial interactions occurring within the Trichodesmium holobiont, especially between Trichodesmium and heterotrophic bacterial epibionts. Yet, the composition of the Trichodesmium holobiont and the processes governing microbial assemblage are not well documented. Here, we used high-resolution 16S rDNA amplicon sequencing to examine the diversity of Trichodesmium and associated epibionts across different ocean regions and colony morphologies (puffs and rafts). Trichodesmium Clade I (i.e., T. thiebautii-like) dominated the colonies in all ocean basins regardless of morphology, although the Trichodesmium community structure significantly varied between morphologies in some regions. On average, Alphaproteobacteria (i.e., Thalassobius), Gammaproteobacteria (i.e., Pseudoalteromonas), Sphingobacteria (i.e., Microscilla and Vibrio) and Flavobacteria dominated the epibiont communities, but community composition and structure significantly differed between regions. Epibionts from the two colony morphologies were taxonomically and functionally distinct in the North Atlantic and North Pacific. These findings suggest that the colony types might define two distinct niches and that epibiont assemblage might be driven in part by selective processes, where epibionts are selected according to their influence on colony metabolism., (© 2016 Society for Applied Microbiology and John Wiley & Sons Ltd.)

Dissolved organic matter (DOM) in the oceans is one of the largest pools of reduced carbon on Earth, comparable in size to the atmospheric CO2 reservoir. A vast number of compounds are present in DOM, and they play important roles in all major element cycles, contribute to the storage of atmospheric CO2 in the ocean, support marine ecosystems, and facilitate interactions between organisms. At the heart of the DOM cycle lie molecular-level relationships between the individual compounds in DOM and the members of the ocean microbiome that produce and consume them. In the past, these connections have eluded clear definition because of the sheer numerical complexity of both DOM molecules and microorganisms. Emerging tools in analytical chemistry, microbiology, and informatics are breaking down the barriers to a fuller appreciation of these connections. Here we highlight questions being addressed using recent methodological and technological developments in those fields and consider how these advances are transforming our understanding of some of the most important reactions of the marine carbon cycle.

A diverse microbial assemblage in the ocean is responsible for nearly half of global primary production. It has been hypothesized and experimentally demonstrated that nutrient loading can stimulate blooms of large eukaryotic phytoplankton in oligotrophic systems. Although central to balancing biogeochemical models, knowledge of the metabolic traits that govern the dynamics of these bloom-forming phytoplankton is limited. We used eukaryotic metatranscriptomic techniques to identify the metabolic basis of functional group-specific traits that may drive the shift between net heterotrophy and autotrophy in the oligotrophic ocean. Replicated blooms were simulated by deep seawater (DSW) addition to mimic nutrient loading in the North Pacific Subtropical Gyre, and the transcriptional responses of phytoplankton functional groups were assayed. Responses of the diatom, haptophyte, and dinoflagellate functional groups in simulated blooms were unique, with diatoms and haptophytes significantly (95% confidence) shifting their quantitative metabolic fingerprint from the in situ condition, whereas dinoflagellates showed little response. Significantly differentially abundant genes identified the importance of colimitation by nutrients, metals, and vitamins in eukaryotic phytoplankton metabolism and bloom formation in this system. The variable transcript allocation ratio, used to quantify transcript reallocation following DSW amendment, differed for diatoms and haptophytes, reflecting the long-standing paradigm of phytoplankton r- and K-type growth strategies. Although the underlying metabolic potential of the large eukaryotic phytoplankton was consistently present, the lack of a bloom during the study period suggests a crucial dependence on physical and biogeochemical forcing, which are susceptible to alteration with changing climate.

Trichodesmium abundance was elevated in certain cyclonic and anticyclonic eddiesEnhancement in cyclonic eddies could be driven by Ekman convergenceAnticyclonic eddies with elevated abundance were anomalously fresh.

Phosphorus in the +5 oxidation state (i.e., phosphate) is the most abundant form of phosphorus in the global ocean. An enigmatic pool of dissolved phosphonate molecules, with phosphorus in the +3 oxidation state, is also ubiquitous; however, cycling of phosphorus between oxidation states has remained poorly constrained. Using simple incubation and chromatography approaches, we measured the rate of the chemical reduction of phosphate to P(III) compounds in the western tropical North Atlantic Ocean. Colonial nitrogen-fixing cyanobacteria in surface waters played a critical role in phosphate reduction, but other classes of plankton, including potentially deep-water archaea, were also involved. These data are consistent with marine geochemical evidence and microbial genomic information, which together suggest the existence of a vast oceanic phosphorus redox cycle., (Copyright © 2015, American Association for the Advancement of Science.)

Diverse communities of marine phytoplankton carry out half of global primary production. The vast diversity of the phytoplankton has long perplexed ecologists because these organisms coexist in an isotropic environment while competing for the same basic resources (e.g., inorganic nutrients). Differential niche partitioning of resources is one hypothesis to explain this "paradox of the plankton," but it is difficult to quantify and track variation in phytoplankton metabolism in situ. Here, we use quantitative metatranscriptome analyses to examine pathways of nitrogen (N) and phosphorus (P) metabolism in diatoms that cooccur regularly in an estuary on the east coast of the United States (Narragansett Bay). Expression of known N and P metabolic pathways varied between diatoms, indicating apparent differences in resource utilization capacity that may prevent direct competition. Nutrient amendment incubations skewed N/P ratios, elucidating nutrient-responsive patterns of expression and facilitating a quantitative comparison between diatoms. The resource-responsive (RR) gene sets deviated in composition from the metabolic profile of the organism, being enriched in genes associated with N and P metabolism. Expression of the RR gene set varied over time and differed significantly between diatoms, resulting in opposite transcriptional responses to the same environment. Apparent differences in metabolic capacity and the expression of that capacity in the environment suggest that diatom-specific resource partitioning was occurring in Narragansett Bay. This high-resolution approach highlights the molecular underpinnings of diatom resource utilization and how cooccurring diatoms adjust their cellular physiology to partition their niche space.

Targeted gene expression using quantitative reverse transcription polymerase chain reaction (qRT-PCR) was employed to track patterns in the expression of genes indicative of nitrogen or phosphorus deficiency in the brown tide-forming alga Aureococcus anophagefferens. During culture experiments, a xanthine/uracil/vitamin C permease (XUV) was upregulated ∼20-fold under nitrogen-deficient conditions relative to a nitrogen-replete control and rapidly returned to nitrogen-replete levels after nitrogen-deficient cells were resupplied with nitrate or ammonium. It was not responsive to phosphorus deficiency. Expression of an inorganic phosphate transporter (PTA3) was enriched ∼10-fold under phosphorus-deficient conditions relative to a phosphorus-replete control, and this signal was rapidly lost upon phosphate resupply. PTA3 was not upregulated by nitrogen deficiency. Natural A. anophagefferens populations from a dense brown tide that occurred in Long Island, NY, in 2009 were assayed for XUV and PTA3 expression and compared with nutrient concentrations over the peak of a bloom. Patterns in XUV expression were consistent with nitrogen-replete growth, never reaching the values observed in N-deficient cultures. PTA3 expression was highest prior to peak bloom stages, reaching expression levels within the range of P-deficient cultures. These data highlight the value of molecular-level assessments of nutrient deficiency and suggest that phosphorus deficiency could play a role in the dynamics of destructive A. anophagefferens blooms., (© 2013 Society for Applied Microbiology and John Wiley & Sons Ltd.)

Transcriptome profiling was performed on the harmful algal bloom-forming pelagophyte Aureococcus anophagefferens strain CCMP 1850 to assess responses to common stressors for dense phytoplankton blooms: low inorganic nitrogen concentrations, low inorganic phosphorus concentrations, low light levels, and a replete control. The de novo assemblies of pooled reads from all treatments reconstructed ~54,000 transcripts using Trinity, and ~31,000 transcripts using ABySS. Comparison to the strain CCMP 1984 genome showed that the majority of the gene models were present in both de novo assemblies and that roughly 95% of contigs from both assemblies mapped to the genome, with Trinity capturing slightly more genome content. Sequence reads were mapped back to the de novo assemblies as well as the gene models and differential expression was analyzed using a Bayesian approach called Analysis of Sequence Counts (ASC). On average, 93% of significantly upregulated transcripts recovered by genome mapping were present in the significantly upregulated pool from both de novo assembly methods. Transcripts related to the transport and metabolism of nitrogen were upregulated in the low nitrogen treatment, transcripts encoding enzymes that hydrolyze organic phosphorus or relieve arsenic toxicity were upregulated in the low phosphorus treatment, and transcripts for enzymes that catabolize organic compounds, restructure lipid membranes, or are involved in sulfolipid biosynthesis were upregulated in the low light treatment. A comparison of this transcriptome to the nutrient regulated transcriptional response of CCMP 1984 identified conserved responses between these two strains. These analyses reveal the transcriptional underpinnings of physiological shifts that could contribute to the ecological success of this species in situ: organic matter processing, metal detoxification, lipid restructuring, and photosynthetic apparatus turnover.

Competing Interests: The authors have declared that no competing interests exist.

Phytoplankton alter their biochemical composition according to nutrient availability, such that their bulk elemental composition varies across oceanic provinces. However, the links between plankton biochemical composition and variation in biogeochemical cycling of nutrients remain largely unknown. In a survey of phytoplankton phosphorus stress in the western North Atlantic, we found that phytoplankton in the phosphorus-depleted subtropical Sargasso Sea were enriched in the biochemical polyphosphate (polyP) compared with nutrient-rich temperate waters, contradicting the canonical oceanographic view of polyP as a luxury phosphorus storage molecule. The enrichment in polyP coincided with enhanced alkaline phosphatase activity and substitution of sulfolipids for phospholipids, which are both indicators of phosphorus stress. Further, polyP appeared to be liberated preferentially over bulk phosphorus from sinking particles in the Sargasso Sea, thereby retaining phosphorus in shallow waters. Thus, polyP cycling may form a feedback loop that attenuates the export of phosphorus when it becomes scarce, contributes bioavailable P for primary production, and supports the export of carbon and nitrogen via sinking particles.

Coccolithophores have influenced the global climate for over 200 million years. These marine phytoplankton can account for 20 per cent of total carbon fixation in some systems. They form blooms that can occupy hundreds of thousands of square kilometres and are distinguished by their elegantly sculpted calcium carbonate exoskeletons (coccoliths), rendering them visible from space. Although coccolithophores export carbon in the form of organic matter and calcite to the sea floor, they also release CO2 in the calcification process. Hence, they have a complex influence on the carbon cycle, driving either CO2 production or uptake, sequestration and export to the deep ocean. Here we report the first haptophyte reference genome, from the coccolithophore Emiliania huxleyi strain CCMP1516, and sequences from 13 additional isolates. Our analyses reveal a pan genome (core genes plus genes distributed variably between strains) probably supported by an atypical complement of repetitive sequence in the genome. Comparisons across strains demonstrate that E. huxleyi, which has long been considered a single species, harbours extensive genome variability reflected in different metabolic repertoires. Genome variability within this species complex seems to underpin its capacity both to thrive in habitats ranging from the equator to the subarctic and to form large-scale episodic blooms under a wide variety of environmental conditions.

Genes that are constitutively expressed across multiple environmental stimuli are crucial to quantifying differentially expressed genes, particularly when employing quantitative reverse transcriptase polymerase chain reaction (RT-qPCR) assays. However, the identification of these potential reference genes in non-model organisms is challenging and is often guided by expression patterns in distantly related organisms. Here, transcriptome datasets from the diatom Thalassiosira pseudonana grown under replete, phosphorus-limited, iron-limited, and phosphorus and iron co-limited nutrient regimes were analyzed through literature-based searches for homologous reference genes, k-means clustering, and analysis of sequence counts (ASC) to identify putative reference genes. A total of 9759 genes were identified and screened for stable expression. Literature-based searches surveyed 18 generally accepted reference genes, revealing 101 homologs in T. pseudonana with variable expression and a wide range of mean tags per million. k-means analysis parsed the whole transcriptome into 15 clusters. The two most stable clusters contained 709 genes, but still had distinct patterns in expression. ASC analyses identified 179 genes that were stably expressed (posterior probability < 0.1 for 1.25 fold change). Genes known to have a stable expression pattern across the test treatments, like actin, were identified in this pool of 179 candidate genes. ASC can be employed on data without biological replicates and was more robust than the k-means approach in isolating genes with stable expression. The intersection of the genes identified through ASC with commonly used reference genes from the literature suggests that actin and ubiquitin ligase may be useful reference genes for T. pseudonana and potentially other diatoms. With the wealth of transcriptome sequence data becoming available, ASC can be easily applied to transcriptome datasets from other phytoplankton to identify reference genes.

Dissolved organic phosphorus (DOP) is increasingly recognized as an important phosphorus source to marine primary producers. Despite its importance, the production rate and fate of DOP is poorly understood. In this study, patterns of DOP production were evaluated by tracking the evolution of DOP during simulated phytoplankton blooms initiated with nutrient amended surface waters, relative to controls, from the Oregon (USA) coastal upwelling system. Nitrogen (N) and phosphorus (P) additions were used to decouple DOP production and hydrolysis by inducing or repressing, respectively, community alkaline phosphatase activity. In order to examine the progression of nutrient uptake and DOP production under upwelling versus relaxation conditions, two experiments were initiated with waters collected during upwelling events, and two with waters collected during relaxation events. Maximum [under (+P) conditions] and minimum [under (+N) conditions] DOP production rates were calculated and applied to in situ DOP levels to evaluate which end-member rate most closely approximates the in situ DOP production rate at the four study sites in this coastal system. Increases in DOP concentration occurred by day-5 in control treatments in all experiments. N treatments displayed increased chlorophyll a, increased alkaline phosphatase activity, and yielded lower net DOP production rates relative to controls, suggesting that DOP levels were depressed as a consequence of increased hydrolysis of bioavailable DOP substrates. Phosphorus additions resulted in a significant net production of DOP at all stations, but no increase in chlorophyll a relative to control treatments. The contrasting patterns in DOP production between treatments suggests that changes in the ambient dissolved inorganic nitrogen:dissolved inorganic phosphorus (DIN:DIP) ratio could exert profound control over DOP production rates in this system. Patterns of DOP production across the different experiments also suggest that bathymetry-driven differences in water residence times can influence DOP cycling. Taken together, these factors may impact the potential export of DOP to offshore ecosystems.

Colonies of the cyanobacterium Trichodesmium are abundant in the oligotrophic ocean, and through their ability to fix both CO(2) and N(2), have pivotal roles in the cycling of carbon and nitrogen in these highly nutrient-depleted environments. Trichodesmium colonies host complex consortia of epibiotic heterotrophic bacteria, and yet, the regulation of nutrient acquisition by these epibionts is poorly understood. We present evidence that epibiotic bacteria in Trichodesmium consortia use quorum sensing (QS) to regulate the activity of alkaline phosphatases (APases), enzymes used by epibionts in the acquisition of phosphate from dissolved-organic phosphorus molecules. A class of QS molecules, acylated homoserine lactones (AHLs), were produced by cultivated epibionts, and adding these AHLs to wild Trichodesmium colonies collected at sea led to a consistent doubling of APase activity. By contrast, amendments of (S)-4,5-dihydroxy-2,3-pentanedione (DPD)-the precursor to the autoinducer-2 (AI-2) family of universal interspecies signaling molecules-led to the attenuation of APase activity. In addition, colonies collected at sea were found by high performance liquid chromatography/mass spectrometry to contain both AHLs and AI-2. Both types of molecules turned over rapidly, an observation we ascribe to quorum quenching. Our results reveal a complex chemical interplay among epibionts using AHLs and AI-2 to control access to phosphate in dissolved-organic phosphorus.

Phosphorus (P) is a critical driver of phytoplankton growth and ecosystem function in the ocean. Diatoms are an abundant class of marine phytoplankton that are responsible for significant amounts of primary production. With the control they exert on the oceanic carbon cycle, there have been a number of studies focused on how diatoms respond to limiting macro and micronutrients such as iron and nitrogen. However, diatom physiological responses to P deficiency are poorly understood. Here, we couple deep sequencing of transcript tags and quantitative proteomics to analyze the diatom Thalassiosira pseudonana grown under P-replete and P-deficient conditions. A total of 318 transcripts were differentially regulated with a false discovery rate of <0.05, and a total of 136 proteins were differentially abundant (p<0.05). Significant changes in the abundance of transcripts and proteins were observed and coordinated for multiple biochemical pathways, including glycolysis and translation. Patterns in transcript and protein abundance were also linked to physiological changes in cellular P distributions, and enzyme activities. These data demonstrate that diatom P deficiency results in changes in cellular P allocation through polyphosphate production, increased P transport, a switch to utilization of dissolved organic P through increased production of metalloenzymes, and a remodeling of the cell surface through production of sulfolipids. Together, these findings reveal that T. pseudonana has evolved a sophisticated response to P deficiency involving multiple biochemical strategies that are likely critical to its ability to respond to variations in environmental P availability.

The toxic arsenate ion can behave as a phosphate analog, and this can result in arsenate toxicity especially in areas with elevated arsenate to phosphate ratios like the surface waters of the ocean gyres. In these systems, cellular arsenate resistance strategies would allow phytoplankton to ameliorate the effects of arsenate transport into the cell. Despite the potential coupling between arsenate and phosphate cycling in oligotrophic marine waters, relatively little is known about arsenate resistance in the nitrogen-fixing marine cyanobacteria that are key components of the microbial community in low nutrient systems. The unicellular diazotroph, Crocosphaera watsonii WH8501, was able to grow at reduced rates with arsenate additions up to 30 nM, and estimated arsenate to phosphate ratios of 6:1. The genome of strain WH8501 contains homologs for arsA, arsH, arsB, and arsC, allowing for the reduction of arsenate to arsenite and the pumping of arsenite out of the cell. The short-term addition of arsenate to the growth medium had no effect on nitrogen fixation. However, arsenate addition did result in the up-regulation of the arsB gene with increasing arsenate concentrations, indicating the induction of the arsenate detoxification response. The arsB gene was also up-regulated by phosphorus stress in concert with a gene encoding the high-affinity phosphate binding protein pstS. Both genes were down-regulated when phosphate was re-fed to phosphorus-stressed cells. A field survey of surface water from the low phosphate western North Atlantic detected expression of C. watsoniiarsB, suggestive of the potential importance of arsenate resistance strategies in this and perhaps other systems.

In low-phosphorus (P) marine systems, phytoplankton replace membrane phospholipids with non-phosphorus lipids, but it is not known how rapidly this substitution occurs. Here, when cells of the model diatom Thalassiosira pseudonana were transferred from P-replete medium to P-free medium, the phospholipid content of the cells rapidly declined within 48 h from 45±0.9 to 21±4.5% of the total membrane lipids; the difference was made up by non-phosphorus lipids. Conversely, when P-limited T. pseudonana were resupplied with P, cells reduced the percentage of their total membrane lipids contributed by a non-phosphorus lipid from 43±1.5 to 7.3±0.9% within 24 h, whereas the contribution by phospholipids rose from 2.2±0.1 to 44±3%. This dynamic phospholipid reservoir contained sufficient P to synthesize multiple haploid genomes, suggesting that phospholipid turnover could be an important P source for cells. Field observations of phytoplankton lipid content may thus reflect short-term changes in P supply and cellular physiology, rather than simply long-term adjustment to the environment.

Harmful algal blooms (HABs) cause significant economic and ecological damage worldwide. Despite considerable efforts, a comprehensive understanding of the factors that promote these blooms has been lacking, because the biochemical pathways that facilitate their dominance relative to other phytoplankton within specific environments have not been identified. Here, biogeochemical measurements showed that the harmful alga Aureococcus anophagefferens outcompeted co-occurring phytoplankton in estuaries with elevated levels of dissolved organic matter and turbidity and low levels of dissolved inorganic nitrogen. We subsequently sequenced the genome of A. anophagefferens and compared its gene complement with those of six competing phytoplankton species identified through metaproteomics. Using an ecogenomic approach, we specifically focused on gene sets that may facilitate dominance within the environmental conditions present during blooms. A. anophagefferens possesses a larger genome (56 Mbp) and has more genes involved in light harvesting, organic carbon and nitrogen use, and encoding selenium- and metal-requiring enzymes than competing phytoplankton. Genes for the synthesis of microbial deterrents likely permit the proliferation of this species, with reduced mortality losses during blooms. Collectively, these findings suggest that anthropogenic activities resulting in elevated levels of turbidity, organic matter, and metals have opened a niche within coastal ecosystems that ideally suits the unique genetic capacity of A. anophagefferens and thus, has facilitated the proliferation of this and potentially other HABs.