The dynamic trophic architecture of open-ocean protist communities revealed through machine-guided metatranscriptomics.

Intricate networks of single-celled eukaryotes (protists) dominate carbon flow in the ocean. Their growth, demise, and interactions with other microorganisms drive the fluxes of biogeochemical elements through marine ecosystems. Mixotrophic protists are capable of both photosynthesis and ingestion of prey and are dominant components of open-ocean planktonic communities. Yet the role of mixotrophs in elemental cycling is obscured by their capacity to act as primary producers or heterotrophic consumers depending on factors that remain largely uncharacterized. Here, we develop and apply a machine learning model that predicts the in situ trophic mode of aquatic protists based on their patterns of gene expression. This approach leverages a public collection of protist transcriptomes as a training set to identify a subset of gene families whose transcriptional profiles predict trophic mode. We applied our model to nearly 100 metatranscriptomes obtained during two oceanographic cruises in the North Pacific and found community-level and population-specific evidence that abundant open-ocean mixotrophic populations shift their predominant mode of nutrient and carbon acquisition in response to natural gradients in nutrient supply and sea surface temperature. Metatranscriptomic data from ship-board incubation experiments revealed that abundant mixotrophic prymnesiophytes from the oligotrophic North Pacific subtropical gyre rapidly remodeled their transcriptome to enhance photosynthesis when supplied with limiting nutrients. Coupling this approach with experiments designed to reveal the mechanisms driving mixotroph physiology provides an avenue toward understanding the ecology of mixotrophy in the natural environment.

Diel transcriptional oscillations of light-sensitive regulatory elements in open-ocean eukaryotic plankton communities.

The 24-h cycle of light and darkness governs daily rhythms of complex behaviors across all domains of life. Intracellular photoreceptors sense specific wavelengths of light that can reset the internal circadian clock and/or elicit distinct phenotypic responses. In the surface ocean, microbial communities additionally modulate nonrhythmic changes in light quality and quantity as they are mixed to different depths. Here, we show that eukaryotic plankton in the North Pacific Subtropical Gyre transcribe genes encoding light-sensitive proteins that may serve as light-activated transcription factors, elicit light-driven electrical/chemical cascades, or initiate secondary messenger-signaling cascades. Overall, the protistan community relies on blue light-sensitive photoreceptors of the cryptochrome/photolyase family, and proteins containing the Light-Oxygen-Voltage (LOV) domain. The greatest diversification occurred within Haptophyta and photosynthetic stramenopiles where the LOV domain was combined with different DNA-binding domains and secondary signal-transduction motifs. Flagellated protists utilize green-light sensory rhodopsins and blue-light helmchromes, potentially underlying phototactic/photophobic and other behaviors toward specific wavelengths of light. Photoreceptors such as phytochromes appear to play minor roles in the North Pacific Subtropical Gyre. Transcript abundance of environmental light-sensitive protein-encoding genes that display diel patterns are found to primarily peak at dawn. The exceptions are the LOV-domain transcription factors with peaks in transcript abundances at different times and putative phototaxis photoreceptors transcribed throughout the day. Together, these data illustrate the diversity of light-sensitive proteins that may allow disparate groups of protists to respond to light and potentially synchronize patterns of growth, division, and mortality within the dynamic ocean environment.

Genome-wide diel growth state transitions in the diatom Thalassiosira pseudonana.

Marine diatoms are important primary producers that thrive in diverse and dynamic environments. They do so, in theory, by sensing changing conditions and adapting their physiology accordingly. Using the model species Thalassiosira pseudonana, we conducted a detailed physiological and transcriptomic survey to measure the recurrent transcriptional changes that characterize typical diatom growth in batch culture. Roughly 40% of the transcriptome varied significantly and recurrently, reflecting large, reproducible cell-state transitions between four principal states: (i) "dawn," following 12 h of darkness; (ii) "dusk," following 12 h of light; (iii) exponential growth and nutrient repletion; and (iv) stationary phase and nutrient depletion. Increases in expression of thousands of genes at the end of the reoccurring dark periods (dawn), including those involved in photosynthesis (e.g., ribulose-1,5-bisphosphate carboxylase oxygenase genes rbcS and rbcL), imply large-scale anticipatory circadian mechanisms at the level of gene regulation. Repeated shifts in the transcript levels of hundreds of genes encoding sensory, signaling, and regulatory functions accompanied the four cell-state transitions, providing a preliminary map of the highly coordinated gene regulatory program under varying conditions. Several putative light sensing and signaling proteins were associated with recurrent diel transitions, suggesting that these genes may be involved in light-sensitive and circadian regulation of cell state. These results begin to explain, in comprehensive detail, how the diatom gene regulatory program operates under varying environmental conditions. Detailed knowledge of this dynamic molecular process will be invaluable for new hypothesis generation and the interpretation of genetic, environmental, and metatranscriptomic data from field studies.

The Phaeodactylum genome reveals the evolutionary history of diatom genomes.

Diatoms are photosynthetic secondary endosymbionts found throughout marine and freshwater environments, and are believed to be responsible for around one-fifth of the primary productivity on Earth. The genome sequence of the marine centric diatom Thalassiosira pseudonana was recently reported, revealing a wealth of information about diatom biology. Here we report the complete genome sequence of the pennate diatom Phaeodactylum tricornutum and compare it with that of T. pseudonana to clarify evolutionary origins, functional significance and ubiquity of these features throughout diatoms. In spite of the fact that the pennate and centric lineages have only been diverging for 90 million years, their genome structures are dramatically different and a substantial fraction of genes ( approximately 40%) are not shared by these representatives of the two lineages. Analysis of molecular divergence compared with yeasts and metazoans reveals rapid rates of gene diversification in diatoms. Contributing factors include selective gene family expansions, differential losses and gains of genes and introns, and differential mobilization of transposable elements. Most significantly, we document the presence of hundreds of genes from bacteria. More than 300 of these gene transfers are found in both diatoms, attesting to their ancient origins, and many are likely to provide novel possibilities for metabolite management and for perception of environmental signals. These findings go a long way towards explaining the incredible diversity and success of the diatoms in contemporary oceans.

An atypical member of the light-harvesting complex stress-related protein family modulates diatom responses to light.

Diatoms are prominent phytoplanktonic organisms that contribute around 40% of carbon assimilation in the oceans. They grow and perform optimally in variable environments, being able to cope with unpredictable changes in the amount and quality of light. The molecular mechanisms regulating diatom light responses are, however, still obscure. Using knockdown Phaeodactylum tricornutum transgenic lines, we reveal the key function of a member of the light-harvesting complex stress-related (LHCSR) protein family, denoted LHCX1, in modulation of excess light energy dissipation. In contrast to green algae, this gene is already maximally expressed in nonstressful light conditions and encodes a protein required for efficient light responses and growth. LHCX1 also influences natural variability in photoresponse, as evidenced in ecotypes isolated from different latitudes that display different LHCX1 protein levels. We conclude, therefore, that this gene plays a pivotal role in managing light responses in diatoms.

Divergent functions of two clades of flavodoxin in diatoms mitigate oxidative stress and iron limitation.

Phytoplankton rely on diverse mechanisms to adapt to the decreased iron bioavailability and oxidative stress-inducing conditions of today's oxygenated oceans, including replacement of the iron-requiring ferredoxin electron shuttle protein with a less-efficient iron-free flavodoxin under iron-limiting conditions. Yet, diatoms transcribe flavodoxins in high-iron regions in contrast to other phytoplankton. Here, we show that the two clades of flavodoxins present within diatoms exhibit a functional divergence, with only clade II flavodoxins displaying the canonical role in acclimation to iron limitation. We created CRISPR/Cas9 knock-outs of the clade I flavodoxin from the model diatom Thalassiosira pseudonana and found that these cell lines are hypersensitive to oxidative stress, while maintaining a wild-type response to iron limitation. Within natural diatom communities, clade I flavodoxin transcript abundance is regulated over the diel cycle rather than in response to iron availability, whereas clade II transcript abundances increase either in iron-limiting regions or under artificially induced iron limitation. The observed functional specialization of two flavodoxin variants within diatoms reiterates two major stressors associated with contemporary oceans and illustrates diatom strategies to flourish in diverse aquatic ecosystems.

Flavobacterial exudates disrupt cell cycle progression and metabolism of the diatom Thalassiosira pseudonana.

Phytoplankton and bacteria form the base of marine ecosystems and their interactions drive global biogeochemical cycles. The effects of bacteria and bacteria-produced compounds on diatoms range from synergistic to pathogenic and can affect the physiology and transcriptional patterns of the interacting diatom. Here, we investigate physiological and transcriptional changes in the marine diatom Thalassiosira pseudonana induced by extracellular metabolites of a known antagonistic bacterium Croceibacter atlanticus. Mono-cultures of C. atlanticus released compounds that inhibited diatom cell division and elicited a distinctive morphology of enlarged cells with increased chloroplast content and enlarged nuclei, similar to what was previously observed when the diatom was co-cultured with live bacteria. The extracellular C. atlanticus metabolites induced transcriptional changes in diatom pathways that include recognition and signaling pathways, cell cycle regulation, carbohydrate and amino acid production, as well as cell wall stability. Phenotypic analysis showed a disruption in the diatom cell cycle progression and an increase in both intra- and extracellular carbohydrates in diatom cultures after bacterial exudate treatment. The transcriptional changes and corresponding phenotypes suggest that extracellular bacterial metabolites, produced independently of direct bacterial-diatom interaction, may modulate diatom metabolism in ways that support bacterial growth.

Complex marine microbial communities partition metabolism of scarce resources over the diel cycle.

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.

Diatom PtCPF1 is a new cryptochrome/photolyase family member with DNA repair and transcription regulation activity.

Members of the cryptochrome/photolyase family (CPF) are widely distributed throughout all kingdoms, and encode photosensitive proteins that typically show either photoreceptor or DNA repair activity. Animal and plant cryptochromes have lost DNA repair activity and now perform specialized photoperceptory functions, for example, plant cryptochromes regulate growth and circadian rhythms, whereas mammalian and insect cryptochromes act as transcriptional repressors that control the circadian clock. However, the functional differentiation between photolyases and cryptochromes is now being questioned. Here, we show that the PtCPF1 protein from the marine diatom Phaeodactylum tricornutum shows 6-4 photoproduct repair activity and can act as a transcriptional repressor of the circadian clock in a heterologous mammalian cell system. Conversely, it seems to have a wide role in blue-light-regulated gene expression in diatoms. The protein might therefore represent a missing link in the evolution of CPFs, and act as a novel ultraviolet/blue light sensor in marine environments.

Diel-Regulated Transcriptional Cascades of Microbial Eukaryotes in the North Pacific Subtropical Gyre.

Open-ocean surface waters host a diverse community of single-celled eukaryotic plankton (protists) consisting of phototrophs, heterotrophs, and mixotrophs. The productivity and biomass of these organisms oscillate over diel cycles, and yet the underlying transcriptional processes are known for few members of the community. Here, we examined a 4-day diel time series of transcriptional abundance profiles for the protist community (0.2-100 µm in cell size) in the North Pacific Subtropical Gyre near Station ALOHA. De novo assembly of poly-A+ selected metatranscriptomes yielded over 30 million contigs with taxonomic and functional annotations assigned to 54 and 25% of translated contigs, respectively. The completeness of the resulting environmental eukaryotic taxonomic bins was assessed, and 48 genera were further evaluated for diel patterns in transcript abundances. These environmental transcriptome bins maintained reproducible temporal partitioning of total gene family abundances, with haptophyte and ochrophyte genera generally showing the greatest diel partitioning of their transcriptomes. The haptophyte Phaeocystis demonstrated the highest proportion of transcript diel periodicity, while most other protists had intermediate levels of periodicity regardless of their trophic status. Dinoflagellates, except for the parasitoid genus Amoebophrya, exhibit the fewest diel oscillations of transcript abundances. Diel-regulated gene families were enriched in key metabolic pathways; photosynthesis, carbon fixation, and fatty acid biosynthesis gene families had peak times concentrated around dawn, while gene families involved in protein turnover (proteasome and protein processing) are most active during the high intensity daylight hours. TCA cycle, oxidative phosphorylation and fatty acid degradation predominantly peaked near dusk. We identified temporal pathway enrichments unique to certain taxa, including assimilatory sulfate reduction at dawn in dictyophytes and signaling pathways at early evening in haptophytes, pointing to possible taxon-specific channels of carbon and nutrients through the microbial community. These results illustrate the synchrony of transcriptional regulation to the diel cycle and how the protist community of the North Pacific Subtropical Gyre structures their transcriptomes to guide the daily flux of matter and energy through the gyre ecosystem.

Fe limitation decreases transcriptional regulation over the diel cycle in the model diatom Thalassiosira pseudonana.

Iron (Fe) is an important growth factor for diatoms and its availability is further restricted by changes in the carbonate chemistry of seawater. We investigated the physiological attributes and transcriptional profiles of the diatom Thalassiosira pseudonana grown on a day: night cycle under different CO2/pH and iron concentrations, that in combination generated available iron (Fe') concentrations of 1160, 233, 58 and 12 pM. We found the light-dark conditions to be the main driver of transcriptional patterns, followed by Fe' concentration and CO2 availability, respectively. At the highest Fe' (1160 pM), 55% of the transcribed genes were differentially expressed between day and night, whereas at the lowest Fe' (12 pM), only 28% of the transcribed genes displayed comparable patterns. While Fe limitation disrupts the diel expression patterns for genes in most central metabolism pathways, the diel expression of light- signaling molecules and glycolytic genes was relatively robust in response to reduced Fe'. Moreover, we identified a non-canonical splicing of transcripts encoding triose-phosphate isomerase, a key-enzyme of glycolysis, generating transcript isoforms that would encode proteins with and without an active site. Transcripts that encoded an active enzyme maintained a diel expression at low Fe', while transcripts that encoded the non-active enzyme lost the diel expression. This work illustrates the interplay between nutrient limitation and transcriptional regulation over the diel cycle. Considering that future ocean conditions will reduce the availability of Fe in many parts of the oceans, our work identifies some of the regulatory mechanisms that may shape future ecological communities.

Sulfonate-based networks between eukaryotic phytoplankton and heterotrophic bacteria in the surface ocean.

In the surface ocean, phytoplankton transform inorganic substrates into organic matter that fuels the activity of heterotrophic microorganisms, creating intricate metabolic networks that determine the extent of carbon recycling and storage in the ocean. Yet, the diversity of organic molecules and interacting organisms has hindered detection of specific relationships that mediate this large flux of energy and matter. Here, we show that a tightly coupled microbial network based on organic sulfur compounds (sulfonates) exists among key lineages of eukaryotic phytoplankton producers and heterotrophic bacterial consumers in the North Pacific Subtropical Gyre. We find that cultured eukaryotic phytoplankton taxa produce sulfonates, often at millimolar internal concentrations. These same phytoplankton-derived sulfonates support growth requirements of an open-ocean isolate of the SAR11 clade, the most abundant group of marine heterotrophic bacteria. Expression of putative sulfonate biosynthesis genes and sulfonate abundances in natural plankton communities over the diel cycle link sulfonate production to light availability. Contemporaneous expression of sulfonate catabolism genes in heterotrophic bacteria highlights active cycling of sulfonates in situ. Our study provides evidence that sulfonates serve as an ecologically important currency for nutrient and energy exchange between microbial autotrophs and heterotrophs, highlighting the importance of organic sulfur compounds in regulating ecosystem function.

Nutrient limitation is the main regulatory factor for carotenoid accumulation and for Psy and Pds steady state transcript levels in Dunaliella salina (Chlorophyta) exposed to high light and salt stress.

Dunaliella salina (Dunal) Teodoresco (1905) is a green unicellular alga able to withstand severe salt, light, and nutrient stress, adaptations necessary to grow in harsh environments such as salt ponds. In response to such growth conditions, this microalga accumulates high amounts of beta-carotene in its single chloroplast. In this study, we show that carotenoid accumulation is consistently inhibited in cells grown in nutrient-supplemented media and exposed either to high-light or medium-low-light conditions. Likewise, carotenogenesis in cells shifted to higher salinity (up to 27% NaCl) under medium-low-light conditions is inhibited by the presence of nutrients. The steady-state levels of transcripts encoding phytoene synthase and phytoene desaturase increased substantially in D. salina cells shifted to high light or high salt under nutrient-limiting conditions, whereas the presence of nutrients inhibited this response. The regulatory effect of nutrient availability on the accumulation of carotenoids and messenger RNA levels of the first two enzymes committed to carotenoid biosynthesis is discussed.

Molecular toolbox for studying diatom biology in Phaeodactylum tricornutum.

Research into diatom biology has now entered the post-genomics era, following the recent completion of the Thalassiosira pseudonana and Phaeodactylum tricornutum whole genome sequences and the establishment of Expressed Sequence Tag (EST) databases. The thorough exploitation of these resources will require the development of molecular tools to analyze and modulate the function of diatom genes in vivo. Towards this objective, we report here the identification of several reference genes that can be used as internal standards for gene expression studies by quantitative real-time PCR (qRT-PCR) in P. tricornutum cells grown over a diel cycle. In addition, we describe a series of diatom expression vectors based on Invitrogen Gateway technology for high-throughput protein tagging and overexpression studies in P. tricornutum. We demonstrate the utility of the diatom Destination vectors for determining the subcellular localization of a protein of interest and for immunodetection. The availability of these new resources significantly enriches the molecular toolbox for P. tricornutum and provides the diatom research community with well defined high-throughput methods for the analysis of diatom genes and proteins in vivo.

Evolutionary origins and functions of the carotenoid biosynthetic pathway in marine diatoms.

Carotenoids are produced by all photosynthetic organisms, where they play essential roles in light harvesting and photoprotection. The carotenoid biosynthetic pathway of diatoms is largely unstudied, but is of particular interest because these organisms have a very different evolutionary history with respect to the Plantae and are thought to be derived from an ancient secondary endosymbiosis between heterotrophic and autotrophic eukaryotes. Furthermore, diatoms have an additional xanthophyll-based cycle for dissipating excess light energy with respect to green algae and higher plants. To explore the origins and functions of the carotenoid pathway in diatoms we searched for genes encoding pathway components in the recently completed genome sequences of two marine diatoms. Consistent with the supplemental xanthophyll cycle in diatoms, we found more copies of the genes encoding violaxanthin de-epoxidase (VDE) and zeaxanthin epoxidase (ZEP) enzymes compared with other photosynthetic eukaryotes. However, the similarity of these enzymes with those of higher plants indicates that they had very probably diversified before the secondary endosymbiosis had occurred, implying that VDE and ZEP represent early eukaryotic innovations in the Plantae. Consequently, the diatom chromist lineage likely obtained all paralogues of ZEP and VDE genes during the process of secondary endosymbiosis by gene transfer from the nucleus of the algal endosymbiont to the host nucleus. Furthermore, the presence of a ZEP gene in Tetrahymena thermophila provides the first evidence for a secondary plastid gene encoded in a heterotrophic ciliate, providing support for the chromalveolate hypothesis. Protein domain structures and expression analyses in the pennate diatom Phaeodactylum tricornutum indicate diverse roles for the different ZEP and VDE isoforms and demonstrate that they are differentially regulated by light. These studies therefore reveal the ancient origins of several components of the carotenoid biosynthesis pathway in photosynthetic eukaryotes and provide information about how they have diversified and acquired new functions in the diatoms.

Isolation and characterization of a stress-inducible Dunaliella salina Lcy-beta gene encoding a functional lycopene beta-cyclase.

The halotolerant green alga Dunaliella salina accumulates large amounts of beta-carotene when exposed to various stress conditions. Although several studies concerning accumulation and biotechnological production of beta-carotene have been published, the molecular basis and regulation of the genes involved in carotenoid biosynthesis in D. salina are still poorly known. In this paper, we report the isolation and regulation of the lycopene beta-cyclase (Lcy-beta) gene by abiotic stress. The function of this gene was determined by heterologous genetic complementation in E. coli. Gene expression and physiological analyses revealed that D. salina Lcy-beta steady-state transcript and carotenoid levels were up-regulated in response to all stress conditions tested (salt, light and nutrient depletion). The results presented here suggest that nutrient availability is a key factor influencing carotenogenesis as well as carotenoid biosynthesis-related gene expression in D. salina.