Empirical bayes analysis of sequencing-based transcriptional profiling without replicates.

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.

Growth on ATP Elicits a P-Stress Response in the Picoeukaryote Micromonas pusilla.

The surface waters of oligotrophic oceans have chronically low phosphate (Pi) concentrations, which renders dissolved organic phosphorus (DOP) an important nutrient source. In the subtropical North Atlantic, cyanobacteria are often numerically dominant, but picoeukaryotes can dominate autotrophic biomass and productivity making them important contributors to the ocean carbon cycle. Despite their importance, little is known regarding the metabolic response of picoeukaryotes to changes in phosphorus (P) source and availability. To understand the molecular mechanisms that regulate P utilization in oligotrophic environments, we evaluated transcriptomes of the picoeukaryote Micromonas pusilla grown under Pi-replete and -deficient conditions, with an additional investigation of growth on DOP in replete conditions. Genes that function in sulfolipid substitution and Pi uptake increased in expression with Pi-deficiency, suggesting cells were reallocating cellular P and increasing P acquisition capabilities. Pi-deficient M. pusilla cells also increased alkaline phosphatase activity and reduced their cellular P content. Cells grown with DOP were able to maintain relatively high growth rates, however the transcriptomic response was more similar to the Pi-deficient response than that seen in cells grown under Pi-replete conditions. The results demonstrate that not all P sources are the same for growth; while M. pusilla, a model picoeukaryote, may grow well on DOP, the metabolic demand is greater than growth on Pi. These findings provide insight into the cellular strategies which may be used to support growth in a stratified future ocean predicted to favor picoeukaryotes.

Genetic indicators of iron limitation in wild populations of Thalassiosira oceanica from the northeast Pacific Ocean.

Assessing the iron (Fe) nutritional status of natural diatom populations has proven challenging as physiological and molecular responses can differ in diatoms of the same genus. We evaluated expression of genes encoding flavodoxin (FLDA1) and an Fe-starvation induced protein (ISIP3) as indicators of Fe limitation in the marine diatom Thalassiosira oceanica. The specificity of the response to Fe limitation was tested in cultures grown under Fe- and macronutrient-deficient conditions, as well as throughout the diurnal light cycle. Both genes showed a robust and specific response to Fe limitation in laboratory cultures and were detected in small volume samples collected from the northeast Pacific, demonstrating the sensitivity of this method. Overall, FLDA1 and ISIP3 expression was inversely related to Fe concentrations and offered insight into the Fe nutritional health of T. oceanica in the field. As T. oceanica is a species tolerant to low Fe, indications of Fe limitation in T. oceanica populations may serve as a proxy for severe Fe stress in the overall diatom community. At two shallow coastal locations, FLD1A and ISIP3 expression revealed Fe stress in areas where dissolved Fe concentrations were high, demonstrating that this approach may be powerful for identifying regions where Fe supply may not be biologically available.

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.

Characterization of putative iron responsive genes as species-specific indicators of iron stress in thalassiosiroid diatoms.

Iron (Fe) availability restricts diatom growth and primary production in large areas of the oceans. It is a challenge to assess the bulk Fe nutritional health of natural diatom populations, since species can differ in their physiological and molecular responses to Fe limitation. We assayed expression of selected genes in diatoms from the Thalassiosira genus to assess their potential utility as species-specific molecular markers to indicate Fe status in natural diatom assemblages. In this study, we compared the expression of the photosynthetic genes encoding ferredoxin (a Fe-requiring protein) and flavodoxin (a Fe-free protein) in culture experiments with Fe replete and Fe stressed Thalassiosira pseudonana (CCMP 1335) isolated from coastal waters and Thalassiosira weissflogii (CCMP 1010) isolated from the open ocean. In T. pseudonana, expression of flavodoxin and ferredoxin genes were not sensitive to Fe status but were found to display diel periodicities. In T. weissflogii, expression of flavodoxin was highly responsive to iron levels and was only detectable when cultures were Fe limited. Flavodoxin genes have been duplicated in most diatoms with available genome data and we show that T. pseudonana has lost its copy related to the Fe-responsive copy in T. weissflogii. We also examined the expression of genes for a putative high affinity, copper (Cu)-dependent Fe uptake system in T. pseudonana. Our results indicate that genes encoding putative Cu transporters, a multi-Cu oxidase, and a Fe reductase are not linked to Fe status. The expression of a second putative Fe reductase increased in Fe limited cultures, but this gene was also highly expressed in Fe replete cultures, indicating it may not be a useful marker in the field. Our findings highlight that Fe metabolism may differ among diatoms even within a genus and show a need to validate responses in different species as part of the development pipeline for genetic markers of Fe status in field populations.