Harriet Alexander

The Transcriptome and Proteome of the Diatom Thalassiosira pseudonana Reveal a Diverse Phosphorus Stress Response

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.

Metatranscriptome analyses indicate resource partitioning between diatoms in the field

Significance Nutrient availability plays a central role in driving the activities and large-scale distributions of phytoplankton, yet there are still fundamental gaps in understanding how phytoplankton metabolize nutrients, like nitrogen (N) and phosphorus (P), and how this metabolic potential is modulated in field populations. Here, we show that cooccurring diatoms in a dynamic coastal marine system have apparent differences in their metabolic capacity to use N and P. Further, bioinformatic approaches enabled the identification and species-specific comparison of resource-responsive (RR) genes. Variation of these RR gene sets highlights the disparate transcriptional responses these species have to the same environment, which likely reflects the role resource partitioning has in facilitating the vast diversity of the phytoplankton. 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.

Probing the evolution, ecology and physiology of marine protists using transcriptomics

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.

Functional group-specific traits drive phytoplankton dynamics in the oligotrophic ocean

Significance Blooms of phytoplankton can shift the ecosystem state of low-nutrient ocean regions between net heterotrophic and autotrophic conditions, increasing carbon sequestration and driving carbon export to the deep sea. Little is known about the traits that govern the timing and magnitude of these bloom events. We used metatranscriptomics to assess phytoplankton functional group-specific metabolic shifts experimentally during simulated blooms in the North Pacific Subtropical Gyre. The results indicate blooms form when phytoplankton are released from limitation by resources (nutrients, vitamins, and trace metals) and that the mechanistic basis for the success of one functional group over another may be driven by how efficiently the transcriptome is modulated following a nutrient pulse. 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.

Identifying reference genes with stable expression from high throughput sequence data.

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.

Virus-host relationships of marine single-celled eukaryotes resolved from metatranscriptomics

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.

Conserved Transcriptional Responses to Nutrient Stress in Bloom-Forming Algae

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.

Variable depth distribution of Trichodesmium clades in the North Pacific Ocean

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.

Transcriptional patterns identify resource controls on the diazotroph Trichodesmium in the Atlantic and Pacific oceans

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

Transcriptional response of the harmful raphidophyte Heterosigma akashiwo to nitrate and phosphate stress

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.