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.

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

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.

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

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.

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.

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.

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.

Shotgun mass spectrometry was used to detect proteins in the harmful alga, Aureococcus anophagefferens, and monitor their relative abundance across nutrient replete (control), phosphate-deficient (-P) and -P refed with phosphate (P-refed) conditions. Spectral counting techniques identified differentially abundant proteins and demonstrated that under phosphate deficiency, A. anophagefferens increases proteins involved in both inorganic and organic phosphorus (P) scavenging, including a phosphate transporter, 5'-nucleotidase, and alkaline phosphatase. Additionally, an increase in abundance of a sulfolipid biosynthesis protein was detected in -P and P-refed conditions. Analysis of the polar membrane lipids showed that cellular concentrations of the sulfolipid sulphoquinovosyldiacylglycerol (SQDG) were nearly two-fold greater in the -P condition versus the control condition, while cellular phospholipids were approximately 8-fold less. Transcript and protein abundances were more tightly coupled for gene products involved in P metabolism compared to those involved in a range of other metabolic functions. Comparison of protein abundances between the -P and P-refed conditions identified differences in the timing of protein degradation and turnover. This suggests that culture studies examining nutrient starvation responses will be valuable in interpreting protein abundance patterns for cellular nutritional status and history in metaproteomic datasets.

Background: Recent technological advancements have made high throughput sequencing an increasingly popular approach for transcriptome analysis. Advantages of sequencing-based transcriptional profiling over microarrays have been reported, including lower technical variability. However, advances in technology do not remove biological variation between replicates and this variation is often neglected in many analyses., Results: We propose an empirical Bayes method, titled Analysis of Sequence Counts (ASC), to detect differential expression based on sequencing technology. ASC borrows information across sequences to establish prior distribution of sample variation, so that biological variation can be accounted for even when replicates are not available. Compared to current approaches that simply tests for equality of proportions in two samples, ASC is less biased towards highly expressed sequences and can identify more genes with a greater log fold change at lower overall abundance., Conclusions: ASC unifies the biological and statistical significance of differential expression by estimating the posterior mean of log fold change and estimating false discovery rates based on the posterior mean. The implementation in R is available at http://www.stat.brown.edu/Zwu/research.aspx.

Alexandrium catenella is widespread in western North America and produces a suite of potent neurotoxins that cause paralytic shellfish poisoning (PSP) in humans and have deleterious impacts on public health and economic resources. There are seasonal PSP-related closures of recreational and commercial shellfisheries in the Puget Sound, but the factors that influence cell distribution, abundance, and relationship to paralytic shellfish toxins (PSTs) in this system are poorly described. Here, a quantitative PCR assay was used to detect A. catenella cells in parallel with state shellfish toxicity testing during the 2006 bloom season at 41 sites from April through October. Over 500,000 A. catenella cells liter(-1) were detected at several stations, with two main pulses of cells driving cell distribution, one in June and the other in August. PSTs over the closure limit of 80 mug of PST 100 per g of shellfish tissue were detected at 26 of the 41 sites. Comparison of cell numbers and PST data shows that shellfish toxicity is preceded by an increase in A. catenella cells in 71% of cases. However, cells were also observed in the absence of PSTs in shellfish, highlighting the complex relationship between A. catenella and the resulting shellfish toxicity. These data provide important information on the dynamics of A. catenella cells in the Puget Sound and are a first step toward assessing the utility of plankton monitoring to augment shellfish toxicity testing in this system.

Trichodesmium are responsible for a large fraction of open ocean nitrogen fixation, and are often found in complex consortia of other microorganisms, including viruses, prokaryotes, microbial eukaryotes and metazoa. We applied a community gene expression (metatranscriptomic) approach to study the patterns of microbial gene utilization within colonies of Trichodesmium collected during a bloom in the Southwest Pacific Ocean in April 2007. The survey generated 5711-day and 5385-night putative mRNA reads. The majority of mRNAs were from the co-occurring microorganisms and not Trichodesmium, including other cyanobacteria, heterotrophic bacteria, eukaryotes and phage. Most transcripts did not share homology with proteins from cultivated microorganisms, but were similar to shotgun sequences and unannotated proteins from open ocean metagenomic surveys. Trichodesmium transcripts were mostly expressed photosynthesis, N(2) fixation and S-metabolism genes, whereas those in the co-occurring microorganisms were mostly involved in genetic information storage and processing. Detection of Trichodesmium genes involved in P uptake and As detoxification suggest that local enrichment of N through N(2) fixation may lead to a P-stress response. Although containing similar dominant transcripts to open ocean metatranscriptomes, the overall pattern of gene expression in Trichodesmium colonies was distinct from free-living pelagic assemblages. The identifiable genes expressed by Trichodesmium and closely associated microorganisms reflect the constraints of life in well-lit and nutrient-poor waters, with biosynthetic investment in nutrient acquisition and cell maintenance, which is in contrast to gene transcription by soil and coastal seawater microbial assemblages. The results provide insight into aggregate microbial communities in contrast to planktonic free-living assemblages that are the focus of other studies.

Phosphorus is an obligate requirement for the growth of all organisms; major biochemical reservoirs of phosphorus in marine plankton include nucleic acids and phospholipids. However, eukaryotic phytoplankton and cyanobacteria (that is, 'phytoplankton' collectively) have the ability to decrease their cellular phosphorus content when phosphorus in their environment is scarce. The biochemical mechanisms that allow phytoplankton to limit their phosphorus demand and still maintain growth are largely unknown. Here we show that phytoplankton, in regions of oligotrophic ocean where phosphate is scarce, reduce their cellular phosphorus requirements by substituting non-phosphorus membrane lipids for phospholipids. In the Sargasso Sea, where phosphate concentrations were less than 10 nmol l-1, we found that only 1.3 +/- 0.6% of phosphate uptake was used for phospholipid synthesis; in contrast, in the South Pacific subtropical gyre, where phosphate was greater than 100 nmol l-1, plankton used 17 6% (ref. 6). Examination of the planktonic membrane lipids at these two locations showed that classes of sulphur- and nitrogen-containing membrane lipids, which are devoid of phosphorus, were more abundant in the Sargasso Sea than in the South Pacific. Furthermore, these non-phosphorus, 'substitute lipids' were dominant in phosphorus-limited cultures of all of the phytoplankton species we examined. In contrast, the marine heterotrophic bacteria we examined contained no substitute lipids and only phospholipids. Thus heterotrophic bacteria, which compete with phytoplankton for nutrients in oligotrophic regions like the Sargasso Sea, appear to have a biochemical phosphorus requirement that phytoplankton avoid by using substitute lipids. Our results suggest that phospholipid substitutions are fundamental biochemical mechanisms that allow phytoplankton to maintain growth in the face of phosphorus limitation.

It is widely appreciated that eukaryotic marine phytoplankton can hydrolyze a variety of compounds within the dissolved organic matter (DOM) pool in marine environments. Herein, cultures and field populations of marine phytoplankton were assayed for beta-N-acetylglucosaminidase activity, a terminal enzyme of chitin degradation. A traditional bulk assay, which can assess hydrolytic rate, but is not cell-specific, was complemented with a cell-specific assay that images the activity associated with single cells using an enzyme labeled fluorescence (ELF) substrate. beta-N-acetylglucosaminidase activity was widespread across various taxa of marine phytoplankton, and activity was observed both under controlled culture conditions and in field populations. The number of cells with enzyme activity varied with the nutritional physiology of the test species in three of the 17 cultures tested. In these three cases the number of cells with activity in the low nutrient medium was higher than in nutrient replete medium. Taken together, these data suggest that a broad group of marine phytoplankton may be a relevant part of chitin-like DOM degradation and should be incorporated into conceptual models of chitin cycling in marine systems.

Through the fixation of atmospheric nitrogen and photosynthesis, marine diazotrophs play a critical role in the global cycling of nitrogen and carbon. Crocosphaera watsonii is a recently described unicellular diazotroph that may significantly contribute to marine nitrogen fixation in tropical environments. One of the many factors that can constrain the growth and nitrogen fixation rates of marine diazotrophs is phosphorus bioavailability. Using genomic and physiological approaches, we examined phosphorus scavenging mechanisms in strains of C. watsonii from both the Atlantic and the Pacific. Observations from the C. watsonii WH8501 genome suggest that this organism has the capacity for high-affinity phosphate transport (e.g., homologs of pstSCAB) in low-phosphate, oligotrophic systems. The pstS gene (high-affinity phosphate binding) is present in strains isolated from both the Atlantic and the Pacific, and its expression was regulated by the exogenous phosphate supply in strain WH8501. Genomic observation also indicated a broad capacity for phosphomonoester hydrolysis (e.g., a putative alkaline phosphatase). In contrast, no clear homologs of genes for phosphonate transport and hydrolysis could be identified. Consistent with these genomic observations, C. watsonii WH8501 is able to grow on phosphomonoesters as a sole source of added phosphorus but not on the phosphonates tested to date. Taken together these data suggest that C. watsonii has a robust capacity for scavenging phosphorus in oligotrophic systems, although this capacity differs from that of other marine cyanobacterial genera, such as Synechococcus, Prochlorococcus, and Trichodesmium.

The factors that control the growth and nitrogen fixation rates of marine diazotrophs such as Trichodesmium have been intensively studied because of the role that these processes have in the global cycling of carbon and nitrogen, and in the sequestration of carbon to the deep sea. Because the phosphate concentrations of many ocean gyres are low, the bioavailability of the larger, chemically heterogeneous pool of dissolved organic phosphorus could markedly influence Trichodesmium physiology. Here we describe the induction, by phosphorus stress, of genes from the Trichodesmium erythraeum IMS101 genome that are predicted to encode proteins associated with the high-affinity transport and hydrolysis of phosphonate compounds by a carbon-phosphorus lyase pathway. We show the importance of these genes through expression analyses with T. erythraeum from the Sargasso Sea. Phosphonates are known to be present in oligotrophic marine systems, but have not previously been considered to be bioavailable to marine diazotrophs. The apparent absence of genes encoding a carbon-phosphorus lyase pathway in the other marine cyanobacterial genomes suggests that, relative to other phytoplankton, Trichodesmium is uniquely adapted for scavenging phosphorus from organic sources. This adaptation may help to explain the prevalence of Trichodesmium in low phosphate, oligotrophic systems.

The abundant and widespread coccolithophore Emiliania huxleyi plays an important role in mediating CO2 exchange between the ocean and the atmosphere through its impact on marine photosynthesis and calcification. Here, we use long serial analysis of gene expression (SAGE) to identify E. huxleyi genes responsive to nitrogen (N) or phosphorus (P) starvation. Long SAGE is an elegant approach for examining quantitative and comprehensive gene expression patterns without a priori knowledge of gene sequences via the detection of 21-bp nucleotide sequence tags. E. huxleyi appears to have a robust transcriptional-level response to macronutrient deficiency, with 42 tags uniquely present or up-regulated twofold or greater in the N-starved library and 128 tags uniquely present or up-regulated twofold or greater in the P-starved library. The expression patterns of several tags were validated with reverse transcriptase PCR. Roughly 48% of these differentially expressed tags could be mapped to publicly available genomic or expressed sequence tag (EST) sequence data. For example, in the P-starved library a number of the tags mapped to genes with a role in P scavenging, including a putative phosphate-repressible permease and a putative polyphosphate synthetase. In short, the long SAGE analyses have (i) identified many new differentially regulated gene sequences, (ii) assigned regulation data to EST sequences with no database homology and unknown function, and (iii) highlighted previously uncharacterized aspects of E. huxleyi N and P physiology. To this end, our long SAGE libraries provide a new public resource for gene discovery and transcriptional analysis in this biogeochemically important marine organism.

Alkaline phosphatase activity is a common marker of phosphate stress in many phytoplankton, but it has been difficult to attribute alkaline phosphatase activity to specific organisms or groups of phytoplankton in the field with traditional biochemical procedures. A new alkaline phosphatase substrate, ELF-97 (enzyme-labeled fluorescence), shows promise in this regard. When a phosphate group is cleaved from the ELF-97 reagent, the remaining molecule precipitates near the site of enzyme activity, thus fluorescently tagging cells with alkaline phosphatase activity. We characterized ELF-97 labeling in axenic cultures of a common dinoflagellate, Prorocentrum minimum, in order to understand ELF-97 labeling dynamics when phosphate nutrition varies. Enzyme activity, as detected by ELF-97 labeling, appears to be induced in late-log- or early-stationary-phase cultures if cells are grown in low-phosphate media and is lost when phosphate-stressed cells are refed with phosphate. ELF-97 appears to label an inducible intracellular alkaline phosphatase in P. minimum based on confocal microscopy studies. This may limit the use of this reagent to organisms that lack high levels of constitutive intracellular phosphatases. After laboratory cultures were characterized, ELF-97 was used to assay field populations of P. minimum in Narragansett Bay during two 1-week periods, and 12 to 100% of the P. minimum cells were labeled. The level of cell labeling was reduced by 3 days of incubation with added inorganic phosphate. Our results indicate that ELF-97 is an excellent new tool for monitoring phytoplankton phosphate stress in the environment when the data are supported by appropriate laboratory studies.