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

Nitrogen source and availability regulate plastic population dynamics in the marine diatom Thalassiosira weissflogii

Variation in nitrogen (N) availability significantly influences population dynamics and the productivity of marine phytoplankton. As N availability in the ocean is conditioned by the N source, it is important to understand the capacity of phytoplankton organisms to adjust their physiology and dynamics under different N conditions. We investigated the growth dynamics of Thalassiosira weissflogii, a coastal diatom, in response to different N sources (Nitrate, NO3-; Ammonium, NH4+; urea, CH4N2O) and availabilities (45 and 5 μM). Our findings demonstrate that T. weissflogii can display plastic adjustments in population dynamics to different N sources. These responses evidenced a greater preference for NH4+ and urea than NO3-, particularly under high N availability. The relative growth rate (μ) is higher (1.18 ± 0.01) under NH4+-high treatment compared to NO3--high (1.01 ± 0.01). The carrying capacity (K) varied only among concentrations, indicating equal N utilization efficiency for biomass production. No effects of N source were detected under the low concentration, suggesting that the preference for NH₄⁺ and urea was diminished by limited nitrogen supply due to potential interactions. These results provide valuable insights into the physiological flexibility of T. weissflogii to varying N conditions, shedding light on the ecological success and resilience of this species in highly variable coastal environments.

Nitrogen uptake rates and phytoplankton composition across contrasting North Atlantic Ocean coastal regimes north and south of Cape Hatteras

Understanding nitrogen (N) uptake rates respect to nutrient availability and the biogeography of phytoplankton communities is crucial for untangling the complexities of marine ecosystems and the physical, biological, and chemical forces shaping them. In the summer of 2016, we conducted measurements of bulk microbial uptake rates for six 15N-labeled substrates: nitrate, nitrite, ammonium, urea, cyanate, and dissolve free amino acids across distinct marine provinces, including the continental shelf of the Mid-and South Atlantic Bights (MAB and SAB), the Slope Sea, and the Gulf Stream, marking the first instance of simultaneously measuring six different N uptake rates in this dynamic region. Total measured N uptake rates were lowest in the Gulf Stream followed by the SAB. Notably, the MAB exhibited significantly higher N uptake rates compared to the SAB, likely due to the excess levels of pre-existing phosphorus present in the MAB. Together, urea and nitrate uptake contributed approximately 50% of the total N uptake across the study region. Although cyanate uptake rates were consistently low, they accounted for up to 11% of the total measured N uptake at some Gulf Stream stations. Phytoplankton groups were identified based on specific pigment markers, revealing a dominance of diatoms in the shelf community, while Synechococcus, Prochlorococcus, and pico-eukaryotes dominated in oligotrophic Gulf Stream waters. The reported uptake rates in this study were mostly in agreement with previous studies conducted in coastal waters of the North Atlantic Ocean. This study suggests there are distinct regional patterns of N uptake in this physically dynamic region, correlating with nutrient availability and phytoplankton community composition. These findings contribute valuable insights into the intricate interplay of biological and chemical factors shaping N dynamics in disparate marine ecosystems.

Phytoplankton thermal trait parameterization alters community structure and biogeochemical processes in a modeled ocean

Phytoplankton exhibit diverse physiological responses to temperature which influence their fitness in the environment and consequently alter their community structure. Here, we explored the sensitivity of phytoplankton community structure to thermal response parameterization in a modelled marine phytoplankton community. Using published empirical data, we evaluated the maximum thermal growth rates (μmax ) and temperature coefficients (Q10 ; the rate at which growth scales with temperature) of six key Phytoplankton Functional Types (PFTs): coccolithophores, cyanobacteria, diatoms, diazotrophs, dinoflagellates, and green algae. Following three well-documented methods, PFTs were either assumed to have (1) the same μmax and the same Q10 (as in to Eppley, 1972), (2) a unique μmax but the same Q10 (similar to Kremer et al., 2017), or (3) a unique μmax and a unique Q10 (following Anderson et al., 2021). These trait values were then implemented within the Massachusetts Institute of Technology biogeochemistry and ecosystem model (called Darwin) for each PFT under a control and climate change scenario. Our results suggest that applying a μmax and Q10 universally across PFTs (as in Eppley, 1972) leads to unrealistic phytoplankton communities, which lack diatoms globally. Additionally, we find that accounting for differences in the Q10 between PFTs can significantly impact each PFT's competitive ability, especially at high latitudes, leading to altered modeled phytoplankton community structures in our control and climate change simulations. This then impacts estimates of biogeochemical processes, with, for example, estimates of export production varying by ~10% in the Southern Ocean depending on the parameterization. Our results indicate that the diversity of thermal response traits in phytoplankton not only shape community composition in the historical and future, warmer ocean, but that these traits have significant feedbacks on global biogeochemical cycles.

Eukaryotic genomes from a global metagenomic data set illuminate trophic modes and biogeography of ocean plankton

IMPORTANCE

Single-celled eukaryotes play ecologically significant roles in the marine environment, yet fundamental questions about their biodiversity, ecological function, and interactions remain. Environmental sequencing enables researchers to document naturally occurring protistan communities, without culturing bias, yet metagenomic and metatranscriptomic sequencing approaches cannot separate individual species from communities. To more completely capture the genomic content of mixed protistan populations, we can create bins of sequences that represent the same organism (metagenome-assembled genomes [MAGs]). We developed the EukHeist pipeline, which automates the binning of population-level eukaryotic and prokaryotic genomes from metagenomic reads. We show exciting insight into what protistan communities are present and their trophic roles in the ocean. Scalable computational tools, like EukHeist, may accelerate the identification of meaningful genetic signatures from large data sets and complement researchers' efforts to leverage MAG databases for addressing ecological questions, resolving evolutionary relationships, and discovering potentially novel biodiversity.

Linking extreme seasonality and gene expression in Arctic marine protists

At high latitudes, strong seasonal differences in light availability affect marine organisms and regulate the timing of ecosystem processes. Marine protists are key players in Arctic aquatic ecosystems, yet little is known about their ecological roles over yearly cycles. This is especially true for the dark polar night period, which up until recently was assumed to be devoid of biological activity. A 12 million transcripts catalogue was built from 0.45 to 10 μm protist assemblages sampled over 13 months in a time series station in an Arctic fjord in Svalbard. Community gene expression was correlated with seasonality, with light as the main driving factor. Transcript diversity and evenness were higher during polar night compared to polar day. Light-dependent functions had higher relative expression during polar day, except phototransduction. 64% of the most expressed genes could not be functionally annotated, yet up to 78% were identified in Arctic samples from Tara Oceans, suggesting that Arctic marine assemblages are distinct from those from other oceans. Our study increases understanding of the links between extreme seasonality and biological processes in pico- and nanoplanktonic protists. Our results set the ground for future monitoring studies investigating the seasonal impact of climate change on the communities of microbial eukaryotes in the High Arctic.

High Growth Rate of Diatoms Explained by Reduced Carbon Requirement and Low Energy Cost of Silica Deposition

The rapid growth of diatoms makes them one of the most pervasive and productive types of plankton in the world's ocean, but the physiological basis for their high growth rates remains poorly understood. Here, we evaluate the factors that elevate diatom growth rates, relative to other plankton, using a steady-state metabolic flux model that computes the photosynthetic C source from intracellular light attenuation and the carbon cost of growth from empirical cell C quotas, across a wide range of cell sizes. For both diatoms and other phytoplankton, growth rates decline with increased cell volume, consistent with observations, because the C cost of division increases with size faster than photosynthesis. However, the model predicts overall higher growth rates for diatoms due to reduced C requirements and the low energetic cost of Si deposition. The C savings from the silica frustule are supported by metatranscriptomic data from Tara Oceans, which show that the abundance of transcripts for cytoskeleton components in diatoms is lower than in other phytoplankton. Our results highlight the importance of understanding the origins of phylogenetic differences in cellular C quotas, and suggest that the evolution of silica frustules may play a critical role in the global dominance of marine diatoms. IMPORTANCE This study addresses a longstanding issue regarding diatoms, namely, their fast growth. Diatoms, which broadly are phytoplankton with silica frustules, are the world's most productive microorganisms and dominate in polar and upwelling regions. Their dominance is largely supported by their high growth rate, but the physiological reasoning behind that characteristic has been obscure. In this study, we combine a quantitative model and metatranscriptomic approaches and show that diatoms' low carbon requirements and low energy costs for silica frustule production are the key factors supporting their fast growth. Our study suggests that the effective use of energy-efficient silica as a cellular structure, instead of carbon, enables diatoms to be the most productive organisms in the global ocean.

Variability in phytoplankton shape dominance in marine pelagic systems: prevalence of different adaptive strategies in the eastern Arabian Sea during the winter monsoon

Phytoplankton communities from pelagic systems were assessed to explore the potential of using commonly used traits (such as cell geometry and taxa) as ecological function indicators from the data generated during the winter monsoon in the eastern Arabian Sea (AS). Altogether, data from two oceanic, i.e., convective mixing influenced non-oligotrophic northeastern-AS (NEAS-O) and Rossby wave-influenced oligotrophic southeastern-AS (SEAS-O) and one coastal (NEAS-C) cruises were utilized to decipher the ecological inferences. Overall phytoplankton shapes showed a high level of redundancy by selecting only a few dominant shapes (5 of 22 shapes), though taxonomic diversity was rich (164 species). The taxonomic and morphological approach adopted revealed high species and shape diversity in NEAS-O than in high-abundance NEAS-C and low-abundance SEAS-O. Also, the shape diversity and dominant shapes (cylinder, elliptic-prism, and prism-on-parallelogram) remained the same in oceans than NEAS-C where combined (cylinder + 2 half-sphere) and simple (elliptic-prism) shapes dominated. Additionally, the Rossby-wave front and its reminiscence in SEAS-O and sea-surface-temperature fronts in NEAS-C favored simple and combine shaped phytoplankton, respectively. The morphological properties assessment revealed that the dominant shapes adapted the strategy to conserve the optimal surface-to-volume ratio (S:V) irrespective of changes in greatest-axial-linear-dimension (GALD) in NEAS-O and SEAS-O but not in NEAS-C. However, the dominant shapes in the NEAS-O and SEAS-O opted for high S:V with low GALD and low S:V with high GALD, respectively, while high S:V with no relation with GALD in NEAS-C suggests the prevalence of different adaptive strategies to cope with the respective hydrographic conditions, particularly nutrient availability.

Using metatranscriptomics to better understand the role of microbial nitrogen cycling in coastal sediment benthic flux denitrification efficiency

Spatial and temporal variability in benthic flux denitrification efficiency occurs across Port Phillip Bay, Australia. Here, we assess the capacity for untargeted metatranscriptomics to resolve spatiotemporal differences in the microbial contribution to benthic nitrogen cycling. The most abundant sediment transcripts assembled were associated with the archaeal nitrifier Nitrosopumilus. In sediments close to external inputs of organic nitrogen, the dominant transcripts were associated with Nitrosopumilus nitric oxide nitrite reduction (nirK). The environmental conditions close to organic nitrogen inputs that select for increased transcription in Nitrosopumilus (amoCAB, nirK, nirS, nmo, hcp) additionally selected for increased transcription of bacterial nitrite reduction (nxrB) and transcripts associated with anammox (hzo) but not denitrification (bacterial nirS/nirk). In sediments that are more isolated from external inputs of organic nitrogen dominant transcripts were associated with nitrous oxide reduction (nosZ) and changes in nosZ transcript abundance were uncoupled from transcriptional profiles associated with archaeal nitrification. Coordinated transcription of coupled community-level nitrification-denitrification was not well supported by metatranscriptomics. In comparison, the abundance of archaeal nirK transcripts were site- and season-specific. This study indicates that the transcription of archaeal nirK in response to changing environmental conditions may be an important and overlooked feature of coastal sediment nitrogen cycling.

To assemble or not to assemble: metagenomic profiling of microbially mediated biogeochemical pathways in complex communities

High-throughput profiling of microbial functional traits involved in various biogeochemical cycling pathways using shotgun metagenomic sequencing has been routinely applied in microbial ecology and environmental science. Multiple bioinformatics data processing approaches are available, including assembly-based (single-sample assembly and multi-sample assembly) and read-based (merged reads and raw data). However, it remains not clear how these different approaches may differ in data analyses and affect result interpretation. In this study, using two typical shotgun metagenome datasets recovered from geographically distant coastal sediments, the performance of different data processing approaches was comparatively investigated from both technical and biological/ecological perspectives. Microbially mediated biogeochemical cycling pathways, including nitrogen cycling, sulfur cycling and B12 biosynthesis, were analyzed. As a result, multi-sample assembly provided the most amount of usable information for targeted functional traits, at a high cost of computational resources and running time. Single-sample assembly and read-based analysis were comparable in obtaining usable information, but the former was much more time- and resource-consuming. Critically, different approaches introduced much stronger variations in microbial profiles than biological differences. However, community-level differences between the two sampling sites could be consistently observed despite the approaches being used. In choosing an appropriate approach, researchers shall balance the trade-offs between multiple factors, including the scientific question, the amount of usable information, computational resources and time cost. This study is expected to provide valuable technical insights and guidelines for the various approaches used for metagenomic data analysis.

Haptophyte communities along the Kuroshio current reveal their geographical sources and ecological traits

Haptophytes are one of the most ecologically successful phytoplankton groups in the modern ocean and tend to maintain balanced and stable communities across various environments. However, little is known about the mechanisms that enable community stability and ecological success. To reveal the community characteristics and interactions among haptophytes, we conducted comprehensive observations from the upstream to downstream regions of the Kuroshio Current. Haptophyte abundance and taxonomy were assessed using quantitative polymerase chain reaction and metabarcoding of 18S rRNA sequences, respectively. The haptophyte community structure changed abruptly at sites on the shelf-slope of the East China Sea, indicating the strong influence of shelf waters with high phytoplankton biomass on downstream communities. Correlation network analysis combined with the phylogeny suggested that haptophytes can coexist with their close relatives, possibly owing to their nutritional flexibility, thereby escaping from resource competition. Consistently, some noncalcifying haptophyte genera with high mixotrophic capacities such as Chrysochromulina constituted a major component of the co-occurrence network, whereas coccolithophores such as Emiliania/Gephyrocapsa were rarely observed. Our study findings suggest that noncalcifying haptophytes play crucial roles in community diversity and stability, and in sustaining the food web structure in the Kuroshio ecosystems.

Strong population genomic structure of the toxic dinoflagellate Alexandrium minutum inferred from meta-transcriptome samples

Despite theoretical expectations, marine microeukaryote population are often highly structured and the mechanisms behind such patterns remain to be elucidated. These organisms display huge census population sizes, yet genotyping usually requires clonal strains originating from single cells, hindering proper population sampling. Estimating allelic frequency directly from population wide samples, without any isolation step, offers an interesting alternative. Here, we validate the use of meta-transcriptome environmental samples to determine the population genetic structure of the dinoflagellate Alexandrium minutum. Strain and meta-transcriptome based results both indicated a strong genetic structure for A. minutum in Western Europe, to the level expected between cryptic species. The presence of numerous private alleles, and even fixed polymorphism, would indicate ancient divergence and absence of gene flow between populations. Single nucleotide polymorphisms (SNPs) displaying strong allele frequency differences were distributed throughout the genome, which might indicate pervasive selection from standing genetic variation (soft selective sweeps). However, a few genomic regions displayed extremely low diversity that could result from the fixation of adaptive de novo mutations (hard selective sweeps) within the populations.

Daily dynamics of contrasting spring algal blooms in Santa Monica Bay (central Southern California Bight)

Protistan algae (phytoplankton) dominate coastal upwelling ecosystems where they form massive blooms that support the world's most important fisheries and constitute an important sink for atmospheric CO₂ . Bloom initiation is well understood, but the biotic and abiotic forces that shape short-term dynamics in community composition are still poorly characterized. Here, high-frequency (daily) changes in relative abundance dynamics of the metabolically active protistan community were followed via expressed 18S V4 rRNA genes (RNA) throughout two algal blooms during the spring of 2018 and 2019 in Santa Monica Bay (central Southern California Bight). A diatom bloom formed after wind-driven, nutrient upwelling events in both years, but different taxa dominated each year. Whereas diatoms bloomed following elevated nutrients and declined after depletion each year, a massive dinoflagellate bloom manifested under relatively low inorganic nitrogen conditions following diatom bloom senescence in 2019 but not 2018. Network analysis revealed associations between diatoms and cercozoan putative parasitic taxa and syndinean parasites during 2019 that may have influenced the demise of the diatoms, and the transition to a dinoflagellate-dominated bloom.

Phytoplankton Responses to Bacterially Regenerated Iron in a Southern Ocean Eddy

In the Subantarctic sector of the Southern Ocean, vertical entrainment of iron (Fe) triggers the seasonal productivity cycle but diminishing physical supply during the spring to summer transition forces microbial assemblages to rapidly acclimate. Here, we tested how phytoplankton and bacteria within an isolated eddy respond to different dissolved Fe (DFe)/ligand inputs. We used three treatments: one that mimicked the entrainment of new DFe (Fe-NEW), another in which DFe was supplied from bacterial regeneration of particles (Fe-REG), and a control with no addition of DFe (Fe-NO). After 6 days, 3.5 (Fe-NO, Fe-NEW) to 5-fold (Fe-REG) increases in Chlorophyll a were observed. These responses of the phytoplankton community were best explained by the differences between the treatments in the amount of DFe recycled during the incubation (Fe-REG, 15% recycled c.f. 40% Fe-NEW, 60% Fe-NO). This additional recycling was more likely mediated by bacteria. By day 6, bacterial production was comparable between Fe-NO and Fe-NEW but was approximately two-fold higher in Fe-REG. A preferential response of phytoplankton (haptophyte-dominated) relative to high nucleic acid (HNA) bacteria was also found in the Fe-REG treatment while the relative proportion of diatoms increased faster in the Fe-NEW and Fe-NO treatments. Comparisons between light and dark incubations further confirmed the competition between picophytoplankton and HNA for DFe. Overall, our results demonstrate great versatility by microorganisms to use different Fe sources that results in highly efficient Fe recycling within surface waters. This study also encourages future research to further investigate the interactions between functional groups of microbes (e.g. HNA and cyanobacteria) to better constraint modeling in Fe and carbon biogeochemical cycles.

Reduced forms of nitrogen control the spatial distribution of phytoplankton communities: The functional winner, dinoflagellates in an anthropogenically polluted estuary

The effects of reduced forms of nitrogen (NH₄⁺ and dissolved organic nitrogen (DON)) on the spatial distribution of diatoms and dinoflagellates in an estuarine-coastal water continuum were investigated from 2015 to 2019. The proportion of non-DIN in total nitrogen was utilized as an indicator of DON along with direct measurements of DON. While NO₃^- originated from Seomjin River, the abundant NH₄⁺ and DON occurred from Gwangyang Bay through Namhae. Diatoms were mostly confined to the upper estuarine system and dinoflagellates dominated in the regions with high levels of NH₄⁺ and DON. Generalized additive models also presented the different responses of diatoms and dinoflagellates to increases in NH₄⁺ and DON. Thus, our results highlight that diatoms dominate in NO₃^--replete water with full access to the source and dinoflagellates take over the ecologically open niche in an anthropogenically polluted estuary with full access to reduced forms of nitrogen.

Proteomic traits vary across taxa in a coastal Antarctic phytoplankton bloom

Production and use of proteins is under strong selection in microbes, but it is unclear how proteome-level traits relate to ecological strategies. We identified and quantified proteomic traits of eukaryotic microbes and bacteria through an Antarctic phytoplankton bloom using in situ metaproteomics. Different taxa, rather than different environmental conditions, formed distinct clusters based on their ribosomal and photosynthetic proteomic proportions, and we propose that these characteristics relate to ecological differences. We defined and used a proteomic proxy for regulatory cost, which showed that SAR11 had the lowest regulatory cost of any taxa we observed at our summertime Southern Ocean study site. Haptophytes had lower regulatory cost than diatoms, which may underpin haptophyte-to-diatom bloom progression in the Ross Sea. We were able to make these proteomic trait inferences by assessing various sources of bias in metaproteomics, providing practical recommendations for researchers in the field. We have quantified several proteomic traits (ribosomal and photosynthetic proteomic proportions, regulatory cost) in eukaryotic and bacterial taxa, which can then be incorporated into trait-based models of microbial communities that reflect resource allocation strategies.

Contrasting Community Composition and Co-Occurrence Relationships of the Active Pico-Sized Haptophytes in the Surface and Subsurface Chlorophyll Maximum Layers of the Arctic Ocean in Summer

Haptophytes (Hacrobia: Haptophyta), which can perform phototrophic, phagotrophic, or mixotrophic nutritional modes, are critical for element cycling in a variety of aquatic ecosystems. However, their diversity, particularly in the changing Arctic Ocean (AO), remains largely unknown. In the present study, the biodiversity, community composition, and co-occurrence networks of pico-sized haptophytes in the surface water and subsurface chlorophyll maximum (SCM) layer of the AO were explored. Our results found higher alpha diversity estimates in the surface water compared with in the SCM based on high-throughput sequencing of haptophyte specific 18S rRNA. The community composition of the surface water was significantly different from that of the SCM, and water temperature was identified as the primary factor shaping the community compositions. Prymnesiales (mostly Chrysochromulina), uncultured Prymnesiophyceae, and Phaeocystis dominated the surface water communities, whereas Phaeocystis dominated the SCM communities, followed by Chrysochromulina, uncultured Prymnesiophyceae, and the remaining taxa. The communities of the surface water and SCM layer developed relatively independent modules in the metacommunity network. Nodes in the surface water were more closely connected to one another than those in the SCM. Network stability analysis revealed that surface water networks were more stable than SCM networks. These findings suggest that SCM communities are more susceptible to environmental fluctuations than those in surface water and that future global changes (e.g., global warming) may profoundly influence the development, persistence, and service of SCM in the AO.

Disentangling the Ecological Processes Shaping the Latitudinal Pattern of Phytoplankton Communities in the Pacific Ocean

Phytoplankton diversity and community compositions vary across spaces and are fundamentally affected by several deterministic (e.g., environmental selection) and stochastic (e.g., ecological drift) processes. How this suite of different processes regulates the biogeography of phytoplankton remains to be comprehensively explored. Using high-throughput sequencing data and null model analysis, we revealed the ecological processes shaping the latitudinal community structure of three major phytoplankton groups (i.e., diatoms, Synechococcus, and haptophytes) across the Pacific Ocean (70°N, 170°W to 35°S, 170°W). At the basin scale, heterogeneous selection (selection under heterogeneous environmental conditions) dominated the assembly processes of all phytoplankton groups; however, its relative importance varied greatly at the climatic zonal scale, explaining the distinct latitudinal α- and β-diversity among phytoplankton groups. Assembly processes in Synechococcus and haptophyte communities were mainly controlled by physical and nutrient factors, respectively. High temperature drove Synechococcus communities to be more deterministic with higher diversity, while haptophyte communities were less environmentally selected at low latitudes due to their wide niche breadth and mixotrophic lifestyle. Diatom communities were overwhelmingly dominated by the selection process but with low correlation of measured environmental factors to their community compositions. This could be attributed to the high growth rate of diatoms, as indicated by their lower site occupation frequency than predicted in the neutral community model. Our study showed that heterogeneous selection is the main force that shaped the biogeography of three key phytoplankton groups in the Pacific Ocean, with a latitudinal variation of relative importance due to the distinct traits among phytoplankton.

IMPORTANCE Phytoplankton are diverse and abundant as primary producers in the ocean, with diversity and community compositions varying spatially. How fundamental processes (e.g., selection, dispersal, and drift) regulate their global biogeography remains to be comprehensively explored. In this study, we disentangled the ecological processes of three key phytoplankton groups (i.e., diatoms, Synechococcus, and haptophytes) along the same latitudinal gradients in the Pacific Ocean. Heterogeneous selection, by promoting species richness and reducing similarity between communities, was the dominant process shaping the communities of each phytoplankton group at the basin scale. However, its relative importance varied greatly among different phytoplankton groups in different climate zones, explaining the uneven latitudinal α- and β-diversity. We also highlight the importance of identifying key factors mediating the relative importance of assembly processes in phytoplankton communities, which will enhance our understanding of their biogeography in the ocean and future patterns under climate changes.

Marine phytoplankton functional types exhibit diverse responses to thermal change

Marine phytoplankton generate half of global primary production, making them essential to ecosystem functioning and biogeochemical cycling. Though phytoplankton are phylogenetically diverse, studies rarely designate unique thermal traits to different taxa, resulting in coarse representations of phytoplankton thermal responses. Here we assessed phytoplankton functional responses to temperature using empirically derived thermal growth rates from four principal contributors to marine productivity: diatoms, dinoflagellates, cyanobacteria, and coccolithophores. Using modeled sea surface temperatures for 1950-1970 and 2080-2100, we explored potential alterations to each group's growth rates and geographical distribution under a future climate change scenario. Contrary to the commonly applied Eppley formulation, our data suggest phytoplankton functional types may be characterized by different temperature coefficients (Q₁₀), growth maxima thermal dependencies, and thermal ranges which would drive dissimilar responses to each degree of temperature change. These differences, when applied in response to global simulations of future temperature, result in taxon-specific projections of growth and geographic distribution, with low-latitude coccolithophores facing considerable decreases and cyanobacteria substantial increases in growth rates. These results suggest that the singular effect of changing temperature may alter phytoplankton global community structure, owing to the significant variability in thermal response between phytoplankton functional types.

Elucidating gene expression adaptation of phylogenetically divergent coral holobionts under heat stress

As coral reefs struggle to survive under climate change, it is crucial to know whether they have the capacity to withstand changing conditions, particularly increasing seawater temperatures. Thermal tolerance requires the integrative response of the different components of the coral holobiont (coral host, algal photosymbiont, and associated microbiome). Here, using a controlled thermal stress experiment across three divergent Caribbean coral species, we attempt to dissect holobiont member metatranscriptome responses from coral taxa with different sensitivities to heat stress and use phylogenetic ANOVA to study the evolution of gene expression adaptation. We show that coral response to heat stress is a complex trait derived from multiple interactions among holobiont members. We identify host and photosymbiont genes that exhibit lineage-specific expression level adaptation and uncover potential roles for bacterial associates in supplementing the metabolic needs of the coral-photosymbiont duo during heat stress. Our results stress the importance of integrative and comparative approaches across a wide range of species to better understand coral survival under the predicted rise in sea surface temperatures.

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.

Microbial community transcriptional patterns vary in response to mesoscale forcing in the North Pacific Subtropical Gyre

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

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

Most organisms coordinate key biological events to coincide with the day/night cycle. These diel oscillations are entrained through the activity of light-sensitive photoreceptors that allow organisms to respond rapidly to changes in light exposure. In the ocean, the plankton community must additionally contend with dramatic changes in the quantity and quality of light over depth. Here, we show that the predominantly blue-light field in the open-ocean environment may have driven expansion of blue light-sensitive regulatory elements in open-ocean eukaryotic plankton derived from secondary and tertiary endosymbiosis. The diel transcription of genes encoding light-sensitive elements indicate that photosynthetic and heterotrophic marine protists respond to and anticipate fluctuating light conditions in the dynamic marine environment.

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.

New paradigm in diatom omics and genetic manipulation

Diatoms are one of the most heterogeneous eukaryotic plankton known for regulating earth's biogeochemical cycles and maintaining the marine ecosystems ever since the late Eocene epoch. The advent of multidisciplinary omics approach has both epitomized and revolutionized the nature of their chimeric genetic toolkit, ecophysiology, and metabolic adaptability as well as their interaction with other communities. In addition, advanced functional annotation of transcriptomic and proteomic data using cutting edge bioinformatics tools together with high-resolution genome-scale mathematical modeling has effectively proven as the catapult in solving genetic bottlenecks in microbial as well as diatom exploration. In this review, a corroborative summation of the robust work done in manipulating, engineering, and sequencing of the diatom genomes besides underpinning the holistic application of omics in transcription and translation has been discussed in order to shrewd their multifarious novel potential in the field of biotechnology and provide an insight into their dynamic evolutionary relevance.

Bistability in oxidative stress response determines the migration behavior of phytoplankton in turbulence

Turbulence has long been known to drive phytoplankton fitness and species succession: motile species dominate in calmer environments and non-motile species in turbulent conditions. Yet a mechanistic understanding of the effect of turbulence on phytoplankton migratory behavior and physiology is lacking. By combining a method to generate turbulent cues, quantification of stress accumulation and physiology, and a mathematical model of stress dynamics, we show that motile phytoplankton use their mechanical stability to sense the intensity of turbulent cues and integrate these cues in time via stress signaling to trigger switches in migratory behavior. The stress-mediated warning strategy we discovered provides a paradigm for how phytoplankton cope with turbulence, thereby potentially governing which species will be successful in a changing ocean.

Turbulence is an important determinant of phytoplankton physiology, often leading to cell stress and damage. Turbulence affects phytoplankton migration both by transporting cells and by triggering switches in migratory behavior, whereby vertically migrating cells can actively invert their direction of migration upon exposure to turbulent cues. However, a mechanistic link between single-cell physiology and vertical migration of phytoplankton in turbulence is currently missing. Here, by combining physiological and behavioral experiments with a mathematical model of stress accumulation and dissipation, we show that the mechanism responsible for the switch in the direction of migration in the marine raphidophyte Heterosigma akashiwo is the integration of reactive oxygen species (ROS) signaling generated by turbulent cues. Within timescales as short as tens of seconds, the emergent downward-migrating subpopulation exhibited a twofold increase in ROS, an indicator of stress, 15% lower photosynthetic efficiency, and 35% lower growth rate over multiple generations compared to the upward-migrating subpopulation. The origin of the behavioral split as a result of a bistable oxidative stress response is corroborated by the observation that exposure of cells to exogenous stressors (H₂O₂, UV-A radiation, or high irradiance), in lieu of turbulence, caused comparable ROS accumulation and an equivalent split into the two subpopulations. By providing a mechanistic link between the single-cell mechanics of swimming and physiology on the one side and the emergent population-scale migratory response and impact on fitness on the other, the ROS-mediated early warning response we discovered contributes to our understanding of phytoplankton community composition in future ocean conditions.

Combined pigment and metatranscriptomic analysis reveals highly synchronized diel patterns of phenotypic light response across domains in the open oligotrophic ocean

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.

Linking regional shifts in microbial genome adaptation with surface ocean biogeochemistry

Linking 'omics measurements with biogeochemical cycles is a widespread challenge in microbial community ecology. Here, we propose applying genomic adaptation as 'biosensors' for microbial investments to overcome nutrient stress. We then integrate this genomic information with a trait-based model to predict regional shifts in the elemental composition of marine plankton communities. We evaluated this approach using metagenomic and particulate organic matter samples from the Atlantic, Indian and Pacific Oceans. We find that our genome-based trait model significantly improves our prediction of particulate C : P (carbon : phosphorus) across ocean regions. Furthermore, we detect previously unrecognized ocean areas of iron, nitrogen and phosphorus stress. In many ecosystems, it can be very challenging to quantify microbial stress. Thus, a carefully calibrated genomic approach could become a widespread tool for understanding microbial responses to environmental changes and the biogeochemical outcomes. This article is part of the theme issue 'Conceptual challenges in microbial community ecology'.

Functional Genomics Differentiate Inherent and Environmentally Influenced Traits in Dinoflagellate and Diatom Communities

Dinoflagellates and diatoms are among the most prominent microeukaryotic plankton groups, and they have evolved different functional traits reflecting their roles within ecosystems. However, links between their metabolic processes and functional traits within different environmental contexts warrant further study. The functional biodiversity of dinoflagellates and diatoms was accessed with metatranscriptomics using Pfam protein domains as proxies for functional processes. Despite the overall geographic similarity of functional responses, abiotic (i.e., temperature and salinity; ~800 Pfam domains) and biotic (i.e., taxonomic group; ~1500 Pfam domains) factors influencing particular functional responses were identified. Salinity and temperature were identified as the main drivers of community composition. Higher temperatures were associated with an increase of Pfam domains involved in energy metabolism and a decrease of processes associated with translation and the sulfur cycle. Salinity changes were correlated with the biosynthesis of secondary metabolites (e.g., terpenoids and polyketides) and signal transduction processes, indicating an overall strong effect on the biota. The abundance of dinoflagellates was positively correlated with nitrogen metabolism, vesicular transport and signal transduction, highlighting their link to biotic interactions (more so than diatoms) and suggesting the central role of species interactions in the evolution of dinoflagellates. Diatoms were associated with metabolites (e.g., isoprenoids and carotenoids), as well as lysine degradation, which highlights their ecological role as important primary producers and indicates the physiological importance of these metabolic pathways for diatoms in their natural environment. These approaches and gathered information will support ecological questions concerning the marine ecosystem state and metabolic interactions in the marine environment.

Species specific gene expression dynamics during harmful algal blooms

Harmful algal blooms are caused by specific members of microbial communities. Understanding the dynamics of these events requires comparing the strategies developed by the problematic species to cope with environmental fluctuations to the ones developed by the other members of the community. During three consecutive years, the meta-transcriptome of micro-eukaryote communities was sequenced during blooms of the toxic dinoflagellate Alexandrium minutum. The dataset was analyzed to investigate species specific gene expression dynamics. Major shifts in gene expression were explained by the succession of different species within the community. Although expression patterns were strongly correlated with fluctuation of the abiotic environment, and more specifically with nutrient concentration, transcripts specifically involved in nutrient uptake and metabolism did not display extensive changes in gene expression. Compared to the other members of the community, A. minutum displayed a very specific expression pattern, with lower expression of photosynthesis transcripts and central metabolism genes (TCA cycle, glucose metabolism, glycolysis…) and contrasting expression pattern of ion transporters across environmental conditions. These results suggest the importance of mixotrophy, cell motility and cell-to-cell interactions during A. minutum blooms.

Diatoms Are Selective Segregators in Global Ocean Planktonic Communities

Diatoms are key phytoplankton in the modern ocean that are involved in numerous biotic interactions, ranging from symbiosis to predation and viral infection, which have considerable effects on global biogeochemical cycles. However, despite recent large-scale studies of plankton, we are still lacking a comprehensive picture of the diversity of diatom biotic interactions in the marine microbial community. Through the ecological interpretation of both inferred microbial association networks and available knowledge on diatom interactions compiled in an open-access database, we propose an ecosystems approach for exploring diatom interactions in the ocean.

Diatoms are a major component of phytoplankton, believed to be responsible for around 20% of the annual primary production on Earth. As abundant and ubiquitous organisms, they are known to establish biotic interactions with many other members of plankton. Through analyses of cooccurrence networks derived from the Tara Oceans expedition that take into account both biotic and abiotic factors in shaping the spatial distributions of species, we show that only 13% of diatom pairwise associations are driven by environmental conditions; the vast majority are independent of abiotic factors. In contrast to most other plankton groups, on a global scale, diatoms display a much higher proportion of negative correlations with other organisms, particularly toward potential predators and parasites, suggesting that their biogeography is constrained by top-down pressure. Genus-level analyses indicate that abundant diatoms are not necessarily the most connected and that species-specific abundance distribution patterns lead to negative associations with other organisms. In order to move forward in the biological interpretation of cooccurrence networks, an open-access extensive literature survey of diatom biotic interactions was compiled, of which 18.5% were recovered in the computed network. This result reveals the extent of what likely remains to be discovered in the field of planktonic biotic interactions, even for one of the best-known organismal groups.

IMPORTANCE Diatoms are key phytoplankton in the modern ocean that are involved in numerous biotic interactions, ranging from symbiosis to predation and viral infection, which have considerable effects on global biogeochemical cycles. However, despite recent large-scale studies of plankton, we are still lacking a comprehensive picture of the diversity of diatom biotic interactions in the marine microbial community. Through the ecological interpretation of both inferred microbial association networks and available knowledge on diatom interactions compiled in an open-access database, we propose an ecosystems approach for exploring diatom interactions in the ocean.

Comparative metatranscriptomic profiling and microRNA sequencing to reveal active metabolic pathways associated with a dinoflagellate bloom

Harmful algal blooms (HABs) have increased as a result of global climate and environmental changes, exerting increasing impacts on the aquatic ecosystem, coastal economy, and human health. Despite great research efforts, our understanding on the drivers of HABs is still limited in part because HAB species' physiology is difficult to probe in situ. Here, we used molecular ecological analyses to characterize a dinoflagellate bloom at Xiamen Harbor, China. Prorocentrum donghaiense was identified as the culprit, which nutrient bioassays showed were not nutrient-limited. Metatranscriptome profiling revealed that P. donghaiense highly expressed genes related to N- and P-nutrient uptake, phagotrophy, energy metabolism (photosynthesis, oxidative phophorylation, and rhodopsin) and carbohydrate metabolism (glycolysis/gluconeogenesis, TCA cycle and pentose phosphate) during the bloom. Many genes in P. donghaiense were up-regulated at night, including phagotrophy and environmental communication genes, and showed active expression in mitosis. Eight microbial defense genes were up-regulated in the bloom compared with previously analyzed laboratory cultures. Furthermore, 76 P. donghaiense microRNA were identified from the bloom, and their target genes exhibited marked differences in amino acid metabolism between the bloom and cultures and the potential of up-regulated antibiotic and cell communication capabilities. These findings, consistent with and complementary to recent reports, reveal major metabolic processes in P. donghaiense potentially important for bloom formation and provide a gene repertoire for developing bloom markers in future research.

Progress and promise of omics for predicting the impacts of climate change on harmful algal blooms

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

Into the bloom: Molecular response of pelagic tunicates to fluctuating food availability

The planktonic tunicates appendicularians and thaliaceans are highly efficient filter feeders on a wide range of prey size including bacteria and have shorter generation times than any other marine grazers. These traits allow some tunicate species to reach high population densities and ensure their success in a favourable environment. However, there are still few studies focusing on which genes and gene pathways are associated with responses of pelagic tunicates to environmental variability. Herein, we present the effect of food availability increase on tunicate community and gene expression at the Marquesas Islands (South-East Pacific Ocean). By using data from the Tara Oceans expedition, we show that changes in phytoplankton density and composition trigger the success of a dominant larvacean species (an undescribed appendicularian). Transcriptional signature to the autotroph bloom suggests key functions in specific physiological processes, i.e., energy metabolism, muscle contraction, membrane trafficking, and proteostasis. The relative abundance of reverse transcription-related Pfams was lower at bloom conditions, suggesting a link with adaptive genetic diversity in tunicates in natural ecosystems. Downstream of the bloom, pelagic tunicates were outcompeted by copepods. Our work represents the first metaomics study of the biological effects of phytoplankton bloom on a key zooplankton taxon.

Using single-cell transcriptomics to understand functional states and interactions in microbial eukaryotes

Understanding the diversity and evolution of eukaryotic microorganisms remains one of the major challenges of modern biology. In recent years, we have advanced in the discovery and phylogenetic placement of new eukaryotic species and lineages, which in turn completely transformed our view on the eukaryotic tree of life. But we remain ignorant of the life cycles, physiology and cellular states of most of these microbial eukaryotes, as well as of their interactions with other organisms. Here, we discuss how high-throughput genome-wide gene expression analysis of eukaryotic single cells can shed light on protist biology. First, we review different single-cell transcriptomics methodologies with particular focus on microbial eukaryote applications. Then, we discuss single-cell gene expression analysis of protists in culture and what can be learnt from these approaches. Finally, we envision the application of single-cell transcriptomics to protist communities to interrogate not only community components, but also the gene expression signatures of distinct cellular and physiological states, as well as the transcriptional dynamics of interspecific interactions. Overall, we argue that single-cell transcriptomics can significantly contribute to our understanding of the biology of microbial eukaryotes. This article is part of a discussion meeting issue 'Single cell ecology'.

Functional Differences in the Blooming Phytoplankton Heterosigma akashiwo and Prorocentrum donghaiense Revealed by Comparative Metaproteomics

Phytoplankton blooms are natural phenomena in the ocean, which are the results of rapid cell growth of some phytoplankton species in a unique environment. However, little is known about the molecular events occurring during the bloom. Here, we compared metaproteomes of two phytoplankton Heterosigma akashiwo and Prorocentrum donghaiense in the coastal East China Sea. H. akashiwo and P. donghaiense accounted for 7.82% and 4.74% of the phytoplankton community protein abundances in the nonbloom sample, whereas they contributed to 60.13% and 78.09%, respectively, in their individual blooming samples. Compared with P. donghaiense, H. akashiwo possessed a significantly higher abundance of light-harvesting complex proteins, carbonic anhydrasem and RuBisCO. The blooming H. akashiwo cells expressed more proteins related to external nutrient acquisition, such as bicarbonate transporter SLC4, ammonium transporter, nitrite transporter, and alkaline phosphatase, while the blooming P. donghaiense cells highly expressed proteins related to extra- and intracellular organic nutrient utilization, such as amino acid transporter, 5'-nucleotidase, acid phosphatase, and tripeptidyl-peptidase. The strong capabilities of light harvesting, as well as acquisition and assimilation of inorganic carbon, nitrogen, and phosphorus, facilitated the formation of the H. akashiwo bloom under the high turbidity and inorganic nutrient-sufficient condition, whereas the competitive advantages in organic nutrient acquisition and reallocation guaranteed the occurrence of the P. donghaiense bloom under the inorganic nutrient-insufficient condition. This study highlights the power of metaproteomics for revealing the underlying molecular behaviors of different coexisting phytoplankton species and advances our knowledge on the formation of phytoplankton blooms.IMPORTANCE A deep understanding of the mechanisms driving bloom formation is a prerequisite for effective bloom management. Metaproteomics was applied in this study to reveal the adaptive and responsive strategies of two coexisting phytoplankton species, H. akashiwo and P. donghaiense, during their bloom periods. Metabolic features and niche divergence in light harvesting, as well as carbon, nitrogen, and phosphorus acquisition and assimilation likely promoted the bloom occurrence under different environments. The molecular behaviors of coexisting bloom-causing species will give clues for bloom monitoring and management in the oceans.

Transcriptome sequencing of a toxic dinoflagellate, Karenia mikimotoi subjected to stress from solar ultraviolet radiation

Solar ultraviolet radiation (UVR) is a stress factor in aquatic environments and may act directly or indirectly on orgnisms in the upper layers of the water column. However, UVR effects are usually species-specific and difficult to extrapolate. Here we use the HAB-forming, toxic dinoflagellate Karenia mikimotoi (which was found to be relatively resistant in previous studies) to investigate its transcriptional responses to a one-week UVR exposure. For this, batch cultures of K. mikimotoi were grown with and without UVR, and their transcriptomes (generated via RNAseq technology) were compared. RNA-seq generated 45.31 million reads, which were further assembled to 202600 unigenes (>300bp). Among these, ca. 61% were annotated with NCBI, NR, GO, KOG, PFAM, Swiss-Prot, and KEGG database. Transcriptomic analysis revealed 722 differentially expressed unigenes (DEGs, defined as being within a |log2 fold change| ≥ 2 and padj < 0.05) responding to solar UVR, which were only 0.36% of all unigenes. 716 unigenes were down-regulated, and only 6 unigenes were up-regulated in the UVR compared to non-UVR treatment. KEGG pathway further analysis revealed DEGs were involved in the different pathway; genes involved in the ribosome, endocytosis and steroid biosynthesis pathways were highly down-regulated, but this was not the case for those involved in the energy metabolisms (including photosynthesis, oxidative phosphorylation) which may contribute to the sustainable growth observed in UVR treatment. The up-regulated expression of both zinc-finger proteins (ZFPs) and ribosomal protein L11 (RPL11) may be one of the acclimated mechanisms against UVR. In addition, this work identified down-regulated genes involved in fatty acid degradation and the hydrophobic branched chain amino acids (e.g., Valine, leucine, and isoleucine), which act as structural components of cell membranes modulating lipid homeostasis or turnover. In conclusion, the present study suggests that the toxic dinoflagellate K. mikimotoi has limited transcriptomic regulation but confirms that it appears as a tolerant species in response to solar UVR. These findings expand current knowledge of gene expression in HAB-forming species in response to natural environment factors such as solar radiation.

Strategies among phytoplankton in response to alleviation of nutrient stress in a subtropical gyre

Despite generally low primary productivity and diatom abundances in oligotrophic subtropical gyres, the North Atlantic Subtropical Gyre (NASG) exhibits significant diatom-driven carbon export on an annual basis. Subsurface pulses of nutrients likely fuel brief episodes of diatom growth, but the exact mechanisms utilized by diatoms in response to these nutrient injections remain understudied within near-natural settings. Here we simulated delivery of subsurface nutrients and compare the response among eukaryotic phytoplankton using a combination of physiological techniques and metatranscriptomics. We show that eukaryotic phytoplankton groups exhibit differing levels of transcriptional responsiveness and expression of orthologous genes in response to release from nutrient limitation. In particular, strategies for use of newly delivered nutrients are distinct among phytoplankton groups. Diatoms channel new nitrate to growth-related strategies while physiological measurements and gene expression patterns of other groups suggest alternative strategies. The gene expression patterns displayed here provide insights into the cellular mechanisms that underlie diatom subsistence during chronic nitrogen-depleted conditions and growth upon nutrient delivery that can enhance carbon export from the surface ocean.

A kleptoplastidic dinoflagellate and the tipping point between transient and fully integrated plastid endosymbiosis

Plastid endosymbiosis has been a major force in the evolution of eukaryotic cellular complexity, but how endosymbionts are integrated is still poorly understood at a mechanistic level. Dinoflagellates, an ecologically important protist lineage, represent a unique model to study this process because dinoflagellate plastids have repeatedly been reduced, lost, and replaced by new plastids, leading to a spectrum of ages and integration levels. Here we describe deep-transcriptomic analyses of the Antarctic Ross Sea dinoflagellate (RSD), which harbors long-term but temporary kleptoplasts stolen from haptophyte prey, and is closely related to dinoflagellates with fully integrated plastids derived from different haptophytes. In some members of this lineage, called the Kareniaceae, their tertiary haptophyte plastids have crossed a tipping point to stable integration, but RSD has not, and may therefore reveal the order of events leading up to endosymbiotic integration. We show that RSD has retained its ancestral secondary plastid and has partitioned functions between this plastid and the kleptoplast. It has also obtained genes for kleptoplast-targeted proteins via horizontal gene transfer (HGT) that are not derived from the kleptoplast lineage. Importantly, many of these HGTs are also found in the related species with fully integrated plastids, which provides direct evidence that genetic integration preceded organelle fixation. Finally, we find that expression of kleptoplast-targeted genes is unaffected by environmental parameters, unlike prey-encoded homologs, suggesting that kleptoplast-targeted HGTs have adapted to posttranscriptional regulation mechanisms of the host.

Diel transcriptional response of a California Current plankton microbiome to light, low iron, and enduring viral infection

Phytoplankton and associated microbial communities provide organic carbon to oceanic food webs and drive ecosystem dynamics. However, capturing those dynamics is challenging. Here, an in situ, semi-Lagrangian, robotic sampler profiled pelagic microbes at 4 h intervals over ~2.6 days in North Pacific high-nutrient, low-chlorophyll waters. We report on the community structure and transcriptional dynamics of microbes in an operationally large size class (>5 μm) predominantly populated by dinoflagellates, ciliates, haptophytes, pelagophytes, diatoms, cyanobacteria (chiefly Synechococcus), prasinophytes (chiefly Ostreococcus), fungi, archaea, and proteobacteria. Apart from fungi and archaea, all groups exhibited 24-h periodicity in some transcripts, but larger portions of the transcriptome oscillated in phototrophs. Periodic photosynthesis-related transcripts exhibited a temporal cascade across the morning hours, conserved across diverse phototrophic lineages. Pronounced silica:nitrate drawdown, a high flavodoxin to ferredoxin transcript ratio, and elevated expression of other Fe-stress markers indicated Fe-limitation. Fe-stress markers peaked during a photoperiodically adaptive time window that could modulate phytoplankton response to seasonal Fe-limitation. Remarkably, we observed viruses that infect the majority of abundant taxa, often with total transcriptional activity synchronized with putative hosts. Taken together, these data reveal a microbial plankton community that is shaped by recycled production and tightly controlled by Fe-limitation and viral activity.

Metatranscriptomic Signatures Associated With Phytoplankton Regime Shift From Diatom Dominance to a Dinoflagellate Bloom

Diatoms and dinoflagellates dominate coastal marine phytoplankton communities as major players of marine biogeochemical cycles and their seasonal succession often leads to harmful algal blooms (HABs). What regulates their respective dominances and the development of the HABs remains elusive. Here we conducted time-sequential metatranscriptomic profiling on a natural assemblage that evolved from diatom dominance to a dinoflagellate bloom to interrogate the underlying major metabolic and ecological drivers. Data reveals similarity between diatoms and dinoflagellates in exhibiting high capacities of energy production, nutrient acquisition, and stress protection in their respective dominance stages. The diatom-to-dinoflagellate succession coincided with an increase in turbidity and sharp declines in silicate and phosphate availability, concomitant with the transcriptomic shift from expression of silicate uptake and urea utilization genes in diatoms to that of genes for light harvesting, diversified phosphorus acquisition and autophagy-based internal nutrient recycling in dinoflagellates. Furthermore, the diatom-dominant community featured strong potential to carbohydrate metabolism and a strikingly high expression of trypsin potentially promoting frustule building. In contrast, the dinoflagellate bloom featured elevated expression of xanthorhodopsin, and antimicrobial defensin genes, indicating potential importance of energy harnessing and microbial defense in bloom development. This study sheds light on mechanisms potentially governing diatom- and dinoflagellate-dominance and regulating bloom development in the natural environment and raises new questions to be addressed in future studies.

Stormwater influences phytoplankton assemblages within the diverse, but impacted Sydney Harbour estuary

Sydney Harbour is subjected to persistent stress associated with anthropogenic activity and global climate change, but is particularly subjected to pulse stress events associated with stormwater input during episodic periods of high rainfall. Photosynthetic microbes underpin metazoan diversity within estuarine systems and are therefore important bioindicators of ecosystem health; yet how stormwater input affects their occurrence and distribution in Sydney Harbour remains poorly understood. We utilised molecular tools (16S/18S rRNA and petB genes) to examine how the phytoplankton community structure (both prokaryotes and eukaryotes) within Sydney Harbour varies between high and low rainfall periods. The relative proportion of phytoplankton sequences was more abundant during the high rainfall period, comprising mainly of diatoms, an important functional group supporting increased productivity within estuarine systems, together with cyanobacteria. Increased spatial variability in the phytoplankton community composition was observed, potentially driven by the steepened physico-chemical gradients associated with stormwater inflow. Conversely, during a low rainfall period, the proportion of planktonic photosynthetic microbes was significantly lower and the persistent phytoplankton were predominantly represented by chlorophyte and dinoflagellate sequences, with lower overall diversity. Differences in phytoplankton composition between the high and low rainfall periods were correlated with temperature, salinity, total nitrogen and silicate. These results suggest that increased frequency of high-rainfall events may change the composition, productivity and health of the estuary. Our study begins to populate the knowledge gap in the phytoplankton community structure and substantial changes associated with transient environmental perturbations, an essential step towards unravelling the dynamics of primary production in a highly urbanised estuarine ecosystem in response to climate change and other anthropogenic stressors.

Daily changes in phytoplankton lipidomes reveal mechanisms of energy storage in the open ocean

Sunlight is the dominant control on phytoplankton biosynthetic activity, and darkness deprives them of their primary external energy source. Changes in the biochemical composition of phytoplankton communities over diel light cycles and attendant consequences for carbon and energy flux in environments remain poorly elucidated. Here we use lipidomic data from the North Pacific subtropical gyre to show that biosynthesis of energy-rich triacylglycerols (TAGs) by eukaryotic nanophytoplankton during the day and their subsequent consumption at night drives a large and previously uncharacterized daily carbon cycle. Diel oscillations in TAG concentration comprise 23 ± 11% of primary production by eukaryotic nanophytoplankton representing a global flux of about 2.4 Pg C yr⁻¹. Metatranscriptomic analyses of genes required for TAG biosynthesis indicate that haptophytes and dinoflagellates are active members in TAG production. Estimates suggest that these organisms could contain as much as 40% more calories at sunset than at sunrise due to TAG production.

Day-night cycles in the biochemical composition of phytoplankton remain poorly understood. Here, Becker et al. use lipidomic and transcriptomic data from the North Pacific subtropical gyre to describe a daily cycle of production and consumption of energy-rich lipids by eukaryotic phytoplankton.

Effects of nutrient enrichment on surface microbial community gene expression in the oligotrophic North Pacific Subtropical Gyre

Marine microbial communities are critical for biogeochemical cycles and the productivity of ocean ecosystems. Primary productivity in the surface ocean is constrained by nutrients which are supplied, in part, by mixing with deeper water. Little is known about the time scales, frequency, or impact of mixing on microbial communities. We combined in situ sampling using the Environmental Sample Processor and a small-scale mixing experiment with lower euphotic zone water to determine how individual populations respond to mixing. Transcriptional responses were measured using the MicroTOOLs (Microbiological Targets for Ocean Observing Laboratories) microarray, which targets all three domains of life and viruses. The experiment showed that mixing substantially affects photosynthetic taxa as expected, but surprisingly also showed that populations respond differently to unfiltered deep water which contains particles (organisms and detritus) compared to filtered deep water that only contains nutrients and viruses, pointing to the impact of biological interactions associated with these events. Comparison between experimental and in situ population transcription patterns indicated that manipulated populations can serve as analogs for natural populations, and that natural populations may be frequently or continuously responding to nutrients from deeper waters. Finally, this study also shows that the microarray approach, which is complementary to metatranscriptomic sequencing, is useful for determining the physiological status of in situ microbial communities.