Metatranscriptomes reveal functional variation in diatom communities from the Antarctic Peninsula

Molecular physiology of Antarctic diatom natural assemblages and bloom event reveal insights into strategies contributing to their ecological success

The continental shelf of the Western Antarctic Peninsula (WAP) is a highly variable system characterized by strong cross-shelf gradients, rapid regional change, and large blooms of phytoplankton, notably diatoms. Rapid environmental changes coincide with shifts in plankton community composition and productivity, food web dynamics, and biogeochemistry. Despite the progress in identifying important environmental factors influencing plankton community composition in the WAP, the molecular basis for their survival in this oceanic region, as well as variations in species abundance, metabolism, and distribution, remains largely unresolved. Across a gradient of physicochemical parameters, we analyzed the metabolic profiles of phytoplankton as assessed through metatranscriptomic sequencing. Distinct phytoplankton communities and metabolisms closely mirrored the strong gradients in oceanographic parameters that existed from coastal to offshore regions. Diatoms were abundant in coastal, southern regions, where colder and fresher waters were conducive to a bloom of the centric diatom, Actinocyclus. Members of this genus invested heavily in growth and energy production; carbohydrate, amino acid, and nucleotide biosynthesis pathways; and coping with oxidative stress, resulting in uniquely expressed metabolic profiles compared to other diatoms. We observed strong molecular evidence for iron limitation in shelf and slope regions of the WAP, where diatoms in these regions employed iron-starvation induced proteins, a geranylgeranyl reductase, aquaporins, and urease, among other strategies, while limiting the use of iron-containing proteins. The metatranscriptomic survey performed here reveals functional differences in diatom communities and provides further insight into the environmental factors influencing the growth of diatoms and their predicted response to changes in ocean conditions.

IMPORTANCE

In the Southern Ocean, phytoplankton must cope with harsh environmental conditions such as low light and growth-limiting concentrations of the micronutrient iron. Using metratranscriptomics, we assessed the influence of oceanographic variables on the diversity of the phytoplankton community composition and on the metabolic strategies of diatoms along the Western Antarctic Peninsula, a region undergoing rapid climate change. We found that cross-shelf differences in oceanographic parameters such as temperature and variable nutrient concentrations account for most of the differences in phytoplankton community composition and metabolism. We opportunistically characterized the metabolic underpinnings of a large bloom of the centric diatom Actinocyclus in coastal waters of the WAP. Our results indicate that physicochemical differences from onshore to offshore are stronger than between southern and northern regions of the WAP; however, these trends could change in the future, resulting in poleward shifts in functional differences in diatom communities and phytoplankton blooms.

Active metabolism and biomass dynamics of biocrusts are shaped by variation in their successional state and seasonal energy sources

Seasonal growth and changes in biomass within communities are the core of ecosystem dynamics. Biocrusts play a prominent role as pioneers in dryland soils. However, the seasonal dynamics of biocrusts remain poorly resolved. In this study, we collected biocrusts across a successional gradient (cyanobacteria, cyanolichen, chlorolichen, and moss-dominated) from southeastern Tengger Desert (China) during the summer and autumn seasons, and explored seasonal changes in metabolism and biomass using multi-omics approaches. We found that Cyanobacteria and Ascomycota were the dominant active taxa and both exhibited higher abundances in autumn. We also found that the dominant primary producers in biocrusts strongly affected community-wide characteristics of metabolism. Along with seasonal differences in light energy utilization, utilization of inorganic energy sources exhibited higher expression in the summer while for organic sources, in the autumn. We found that overall metabolism was significantly regulated by the ratio of intracellular to extracellular polymer degradation, and affected by NO₃^-, PO₄^3- and EC (in the summer)/NO₂^- (in the autumn). In summary, biocrust growth varied with seasonal variation in light energy utilization and complementary chemical energy sources, with the most suitable season varying with biocrust successional type.

Diatoms and Their Microbiomes in Complex and Changing Polar Oceans

Diatoms, a key group of polar marine microbes, support highly productive ocean ecosystems. Like all life on earth, diatoms do not live in isolation, and they are therefore under constant biotic and abiotic pressures which directly influence their evolution through natural selection. Despite their importance in polar ecosystems, polar diatoms are understudied compared to temperate species. The observed rapid change in the polar climate, especially warming, has created increased research interest to discover the underlying causes and potential consequences on single species to entire ecosystems. Next-Generation Sequencing (NGS) technologies have greatly expanded our knowledge by revealing the molecular underpinnings of physiological adaptations to polar environmental conditions. Their genomes, transcriptomes, and proteomes together with the first eukaryotic meta-omics data of surface ocean polar microbiomes reflect the environmental pressures through adaptive responses such as the expansion of protein families over time as a consequence of selection. Polar regions and their microbiomes are inherently connected to climate cycles and their feedback loops. An integrated understanding built on "omics" resources centered around diatoms as key primary producers will enable us to reveal unifying concepts of microbial co-evolution and adaptation in polar oceans. This knowledge, which aims to relate past environmental changes to specific adaptations, will be required to improve climate prediction models for polar ecosystems because it provides a unifying framework of how interacting and co-evolving biological communities might respond to future environmental change.

TRAPID 2.0: a web application for taxonomic and functional analysis of de novo transcriptomes

Advances in high-throughput sequencing have resulted in a massive increase of RNA-Seq transcriptome data. However, the promise of rapid gene expression profiling in a specific tissue, condition, unicellular organism or microbial community comes with new computational challenges. Owing to the limited availability of well-resolved reference genomes, de novo assembled (meta)transcriptomes have emerged as popular tools for investigating the gene repertoire of previously uncharacterized organisms. Yet, despite their potential, these datasets often contain fragmented or contaminant sequences, and their analysis remains difficult. To alleviate some of these challenges, we developed TRAPID 2.0, a web application for the fast and efficient processing of assembled transcriptome data. The initial processing phase performs a global characterization of the input data, providing each transcript with several layers of annotation, comprising structural, functional, and taxonomic information. The exploratory phase enables downstream analyses from the web application. Available analyses include the assessment of gene space completeness, the functional analysis and comparison of transcript subsets, and the study of transcripts in an evolutionary context. A comparison with similar tools highlights TRAPID’s unique features. Finally, analyses performed within TRAPID 2.0 are complemented by interactive data visualizations, facilitating the extraction of new biological insights, as demonstrated with diatom community metatranscriptomes.

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.

Microbial activity during a coastal phytoplankton bloom on the Western Antarctic Peninsula in late summer

Phytoplankton biomass during the austral summer is influenced by freezing and melting cycles as well as oceanographic processes that enable nutrient redistribution in the West Antarctic Peninsula (WAP). Microbial functional capabilities, metagenomic and metatranscriptomic activities as well as inorganic 13C- and 15N-assimilation rates were studied in the surface waters of Chile Bay during two contrasting summer periods in 2014. Concentrations of Chlorophyll a (Chla) varied from 0.3 mg m-3 in February to a maximum of 2.5 mg m-3 in March, together with a decrease in nutrients; however, nutrients were never depleted. The microbial community composition remained similar throughout both sampling periods; however, microbial abundance and activity changed with Chla levels. An increased biomass of Bacillariophyta, Haptophyceae and Cryptophyceae was observed along with night-grazing activity of Dinophyceae and ciliates (Alveolates). During high Chla conditions, HCO3- uptake rates during daytime incubations increased 5-fold (>2516 nmol C L-1 d-1), and increased photosynthetic transcript numbers that were mainly associated with cryptophytes; meanwhile night time NO3- (>706 nmol N L-1 d-1) and NH4+ (41.7 nmol N L-1 d-1) uptake rates were 2- and 3-fold higher, respectively, due to activity from Alpha-/Gammaproteobacteria and Bacteroidetes (Flavobacteriia). Due to a projected acceleration in climate change in the WAP, this information is valuable for predicting the composition and functional changes in Antarctic microbial communities.

Multi-metal tolerance of von Willebrand factor type D domain isolated from metal contaminated site by metatranscriptomics approach

Environmental pollution through heavy metals is an upcoming universal problem that relentlessly endangers human health, biodiversity and ecosystems. Hence remediating these heavy metal pollutants from the environment by engineering soil microbiome through metatranscriptomics is befitting reply. In the present investigation, we have constructed size fractionated cDNA libraries from eukaryotic mRNA of cadmium (Cd) contaminated soil and screened for Cd tolerant genes by yeast complementation system by using Cd sensitive ycf1^Δ mutant. We are reporting one of the transformants PLCe10 (from library C, 1-4 kb) with potential tolerance towards Cd toxicity (40 μM-80 μM). Sequence analysis of PLCe10 transcript showed homology to von Willebrand factor type D domain (VWD) of vitellogenin-6 of Ascaris suum encoding 338 amino acids peptide. qPCR analysis revealed that PLCe10 induced in presence of Cd (32 fold) and also accumulated maximum amount of Cd at 60 μM Cd. This cDNA was further tested for its tolerance against other heavy metals like copper (Cu), zinc (Zn) and cobalt (Co). Heterologous complementation assays of cDNA PLCe10 showed a range of tolerance to Cu (150 μM-500 μM), Zn (10 mM-12 mM) and Co (2-4 mM). Results of the present study suggest that cDNA PLCe10 is one of the functional eukaryotic heavy metal tolerant genes present among the soil microbial community and could be exploited to rehabilitate metal contaminated sites.

Shifting metabolic priorities among key protistan taxa within and below the euphotic zone

A metatranscriptome study targeting the protistan community was conducted off the coast of Southern California, at the San Pedro Ocean Time-series station at the surface, 150 m (oxycline), and 890 m to link putative metabolic patterns to distinct protistan lineages. Comparison of relative transcript abundances revealed depth-related shifts in the nutritional modes of key taxonomic groups. Eukaryotic gene expression in the sunlit surface environment was dominated by phototrophs, such as diatoms and chlorophytes, and high abundances of transcripts associated with synthesis pathways (e.g., photosynthesis, carbon fixation, fatty acid synthesis). Sub-euphotic depths (150 and 890 m) exhibited strong contributions from dinoflagellates and ciliates, and were characterized by transcripts relating to digestion or intracellular nutrient recycling (e.g., breakdown of fatty acids and V-type ATPases). These transcriptional patterns underlie the distinct nutritional modes of ecologically important protistan lineages that drive marine food webs, and provide a framework to investigate trophic dynamics across diverse protistan communities.

Ion and metabolite transport in the chloroplast of algae: lessons from land plants

Chloroplasts are endosymbiotic organelles and play crucial roles in energy supply and metabolism of eukaryotic photosynthetic organisms (algae and land plants). They harbor channels and transporters in the envelope and thylakoid membranes, mediating the exchange of ions and metabolites with the cytosol and the chloroplast stroma and between the different chloroplast subcompartments. In secondarily evolved algae, three or four envelope membranes surround the chloroplast, making more complex the exchange of ions and metabolites. Despite the importance of transport proteins for the optimal functioning of the chloroplast in algae, and that many land plant homologues have been predicted, experimental evidence and molecular characterization are missing in most cases. Here, we provide an overview of the current knowledge about ion and metabolite transport in the chloroplast from algae. The main aspects reviewed are localization and activity of the transport proteins from algae and/or of homologues from other organisms including land plants. Most chloroplast transporters were identified in the green alga Chlamydomonas reinhardtii, reside in the envelope and participate in carbon acquisition and metabolism. Only a few identified algal transporters are located in the thylakoid membrane and play role in ion transport. The presence of genes for putative transporters in green algae, red algae, diatoms, glaucophytes and cryptophytes is discussed, and roles in the chloroplast are suggested. A deep knowledge in this field is required because algae represent a potential source of biomass and valuable metabolites for industry, medicine and agriculture.

Comparative Metatranscriptomics of Wheat Rhizosphere Microbiomes in Disease Suppressive and Non-suppressive Soils for Rhizoctonia solani AG8

The soilborne fungus Rhizoctonia solani anastomosis group (AG) 8 is a major pathogen of grain crops resulting in substantial production losses. In the absence of resistant cultivars of wheat or barley, a sustainable and enduring method for disease control may lie in the enhancement of biological disease suppression. Evidence of effective biological control of R. solani AG8 through disease suppression has been well documented at our study site in Avon, South Australia. A comparative metatranscriptomic approach was applied to assess the taxonomic and functional characteristics of the rhizosphere microbiome of wheat plants grown in adjacent fields which are suppressive and non-suppressive to the plant pathogen R. solani AG8. Analysis of 12 rhizosphere metatranscriptomes (six per field) was undertaken using two bioinformatic approaches involving unassembled and assembled reads. Differential expression analysis showed the dominant taxa in the rhizosphere based on mRNA annotation were Arthrobacter spp. and Pseudomonas spp. for non-suppressive samples and Stenotrophomonas spp. and Buttiauxella spp. for the suppressive samples. The assembled metatranscriptome analysis identified more differentially expressed genes than the unassembled analysis in the comparison of suppressive and non-suppressive samples. Suppressive samples showed greater expression of a polyketide cyclase, a terpenoid biosynthesis backbone gene (dxs) and many cold shock proteins (csp). Non-suppressive samples were characterised by greater expression of antibiotic genes such as non-heme chloroperoxidase (cpo) which is involved in pyrrolnitrin synthesis, and phenazine biosynthesis family protein F (phzF) and its transcriptional activator protein (phzR). A large number of genes involved in detoxifying reactive oxygen species (ROS) and superoxide radicals (sod, cat, ahp, bcp, gpx1, trx) were also expressed in the non-suppressive rhizosphere samples most likely in response to the infection of wheat roots by R. solani AG8. Together these results provide new insight into microbial gene expression in the rhizosphere of wheat in soils suppressive and non-suppressive to R. solani AG8. The approach taken and the genes involved in these functions provide direction for future studies to determine more precisely the molecular interplay of plant-microbe-pathogen interactions with the ultimate goal of the development of management options that promote beneficial rhizosphere microflora to reduce R. solani AG8 infection of crops.

Proton-pumping rhodopsins are abundantly expressed by microbial eukaryotes in a high-Arctic fjord

Proton-pumping rhodopsins provide an alternative pathway to photosynthesis by which solar energy can enter the marine food web. Rhodopsin genes are widely found in marine bacteria, also in the Arctic, and were recently reported from several eukaryotic lineages. So far, little is known about rhodopsin expression in Arctic eukaryotes. In this study, we used metatranscriptomics and 18S rDNA tag sequencing to examine the mid-summer function and composition of marine protists (size 0.45-10 µm) in the high-Arctic Billefjorden (Spitsbergen), especially focussing on the expression of microbial proton-pumping rhodopsins. Rhodopsin transcripts were highly abundant, at a level similar to that of genes involved in photosynthesis. Phylogenetic analyses placed the environmental rhodopsins within disparate eukaryotic lineages, including dinoflagellates, stramenopiles, haptophytes and cryptophytes. Sequence comparison indicated the presence of several functional types, including xanthorhodopsins and a eukaryotic clade of proteorhodopsin. Transcripts belonging to the proteorhodopsin clade were also abundant in published metatranscriptomes from other oceanic regions, suggesting a global distribution. The diversity and abundance of rhodopsins show that these light-driven proton pumps play an important role in Arctic microbial eukaryotes. Understanding this role is imperative to predicting the future of the Arctic marine ecosystem faced by a changing light climate due to diminishing sea-ice.

Evaluating trace element bioavailability and potential transfer into marine food chains using immobilised diatom model species Phaeodactylum tricornutum, on King George Island, Antarctica

In order to evaluate trace element bioavailability and potential transfer into marine food chains in human impacted areas of the Fildes Peninsula (King George Island, South Shetland Islands Archipelago), element levels (Cr, Ni, Cu, Zn, Cd, and Pb) were determined in water, sediments, phytoplankton, and in diatom Phaeodactylum tricornutum Bohlin (Bacillariophyceae) cells immobilised in alginate and exposed to water and sediments, from the Bellingshausen Dome (reference site) and Ardley Cove (human impacted area), during January 2014. High element concentrations in exposed P. tricornutum indicated element mobilisation from sediments into the water. Levels in exposed cells reflected the sediment element content pattern, comparable to those found in phytoplankton, supporting phytoplankton as an important path of trace element entry into marine food chains. This study clearly shows immobilised P. tricornutum as good proxy of phytoplankton concerning element accumulation efficiency, and an effective tool to monitor trace element contamination in polar coastal ecosystems.

High particulate iron(II) content in glacially sourced dusts enhances productivity of a model diatom

Little is known about the bioavailability of iron (Fe) in natural dusts and the impact of dust mineralogy on Fe utilization by photosynthetic organisms. Variation in the supply of bioavailable Fe to the ocean has the potential to influence the global carbon cycle by modulating primary production in the Southern Ocean. Much of the dust deposited across the Southern Ocean is sourced from South America, particularly Patagonia, where the waxing and waning of past and present glaciers generate fresh glaciogenic material that contrasts with aged and chemically weathered nonglaciogenic sediments. We show that these two potential sources of modern-day dust are mineralogically distinct, where glaciogenic dust sources contain mostly Fe(II)-rich primary silicate minerals, and nearby nonglaciogenic dust sources contain mostly Fe(III)-rich oxyhydroxide and Fe(III) silicate weathering products. In laboratory culture experiments, Phaeodactylum tricornutum, a well-studied coastal model diatom, grows more rapidly, and with higher photosynthetic efficiency, with input of glaciogenic particulates compared to that of nonglaciogenic particulates due to these differences in Fe mineralogy. Monod nutrient accessibility models fit to our data suggest that particulate Fe(II) content, rather than abiotic solubility, controls the Fe bioavailability in our Fe fertilization experiments. Thus, it is possible for this diatom to access particulate Fe in dusts by another mechanism besides uptake of unchelated Fe (Fe') dissolved from particles into the bulk solution. If this capability is widespread in the Southern Ocean, then dusts deposited to the Southern Ocean in cold glacial periods are likely more bioavailable than those deposited in warm interglacial periods.

Discovering Protein-Coding Genes from the Environment: Time for the Eukaryotes?

Eukaryotic microorganisms from diverse environments encompass a large number of taxa, many of them still unknown to science. One strategy to mine these organisms for genes of biotechnological relevance is to use a pool of eukaryotic mRNA directly extracted from environmental samples. Recent reports demonstrate that the resulting metatranscriptomic cDNA libraries can be screened by expression in yeast for a wide range of genes and functions from many of the different eukaryotic taxa. In combination with novel emerging high-throughput technologies, we anticipate that this approach should contribute to exploring the functional diversity of the eukaryotic microbiota.

Genome complexity of harmful microalgae

During the past decade, next generation sequencing (NGS) technologies have provided new insights into the diversity, dynamics, and metabolic pathways of natural microbial communities. But, these new techniques face challenges related to the genome size and level of genome complexity of the species under investigation. Moreover, the coverage depth and the short-read length achieved by NGS based approaches also represent a major challenge for assembly. These factors could limit the use of these high-throughput sequencing methods for species lacking a reference genome and characterized by a high level of complexity. In the present work, the evolutionary history, mainly consisting of gene transfer events from bacteria and unicellular eukaryotes to microalgae, including harmful species, is discussed and reviewed as it relates to NGS application in microbial communities, with a particular focus on harmful algal bloom species and dinoflagellates. In the context of genetic population studies, genotyping-by-sequencing (GBS), an NGS based approach, could be used for the discovery and analysis of single nucleotide polymorphisms (SNPs). The NGS technologies are still relatively new and require further improvement. Specifically, there is a need to develop and standardize tools and approaches to handle large data sets, which have to be used for the majority of HAB species characterized by evolutionary highly dynamic genomes.

High Prevalence of Gammaproteobacteria in the Sediments of Admiralty Bay and North Bransfield Basin, Northwestern Antarctic Peninsula

Microorganisms dominate most Antarctic marine ecosystems, in terms of biomass and taxonomic diversity, and play crucial role in ecosystem functioning due to their high metabolic plasticity. Admiralty Bay is the largest bay on King George Island (South Shetland Islands, Antarctic Peninsula) and a combination of hydro-oceanographic characteristics (bathymetry, sea ice and glacier melting, seasonal entrance of water masses, turbidity, vertical fluxes) create conditions favoring organic carbon deposition on the seafloor and microbial activities. We sampled surface sediments from 15 sites across Admiralty Bay (100–502 m total depth) and the adjacent North Bransfield Basin (693–1147 m), and used the amplicon 454-sequencing of 16S rRNA gene tags to compare the bacterial composition, diversity, and microbial community structure across environmental parameters (sediment grain size, pigments and organic nutrients) between the two areas. Marine sediments had a high abundance of heterotrophic Gammaproteobacteria (92.4% and 83.8% inside and outside the bay, respectively), followed by Alphaproteobacteria (2.5 and 5.5%), Firmicutes (1.5 and 1.6%), Bacteroidetes (1.1 and 1.7%), Deltaproteobacteria (0.8 and 2.5%) and Actinobacteria (0.7 and 1.3%). Differences in alpha-diversity and bacterial community structure were found between the two areas, reflecting the physical and chemical differences in the sediments, and the organic matter input.

Characterization of CRISPR RNA transcription by exploiting stranded metatranscriptomic data

CRISPR-Cas systems are bacterial adaptive immune systems, each typically composed of a locus of cas genes and a CRISPR array of spacers flanked by repeats. Processed transcripts of CRISPR arrays (crRNAs) play important roles in the interference process mediated by these systems, guiding targeted immunity. Here we developed computational approaches that allow us to characterize the expression of many CRISPRs in their natural environments, using community RNA-seq (metatranscriptomic) data. By exploiting public human gut metatranscriptomic data sets, we studied the expression of 56 repeat-sequence types of CRISPRs, revealing that most CRISPRs are transcribed in one direction (producing crRNAs). In rarer cases, including a type II system associated with Bacteroides fragilis, CRISPRs are transcribed in both directions. Type III CRISPR-Cas systems were found in the microbiomes, but metatranscriptomic reads were barely found for their CRISPRs. We observed individual-level variation of the crRNA transcription, and an even greater transcription of a CRISPR from the antisense strand than the crRNA strand in one sample. The orientations of CRISPR expression implicated by metatranscriptomic data are largely in agreement with prior predictions for CRISPRs, with exceptions. Our study shows the promise of exploiting community RNA-seq data for investigating the transcription of CRISPR-Cas systems.

Multisignal control of expression of the LHCX protein family in the marine diatom Phaeodactylum tricornutum

Diatoms are phytoplanktonic organisms that grow successfully in the ocean where light conditions are highly variable. Studies of the molecular mechanisms of light acclimation in the marine diatom Phaeodactylum tricornutum show that carotenoid de-epoxidation enzymes and LHCX1, a member of the light-harvesting protein family, both contribute to dissipate excess light energy through non-photochemical quenching (NPQ). In this study, we investigate the role of the other members of the LHCX family in diatom stress responses. Our analysis of available genomic data shows that the presence of multiple LHCX genes is a conserved feature of diatom species living in different ecological niches. Moreover, an analysis of the levels of four P. tricornutum LHCX transcripts in relation to protein expression and photosynthetic activity indicates that LHCXs are differentially regulated under different light intensities and nutrient starvation, mostly modulating NPQ capacity. We conclude that multiple abiotic stress signals converge to regulate the LHCX content of cells, providing a way to fine-tune light harvesting and photoprotection. Moreover, our data indicate that the expansion of the LHCX gene family reflects functional diversification of its members which could benefit cells responding to highly variable ocean environments.

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

A diverse microbial assemblage in the ocean is responsible for nearly half of global primary production. It has been hypothesized and experimentally demonstrated that nutrient loading can stimulate blooms of large eukaryotic phytoplankton in oligotrophic systems. Although central to balancing biogeochemical models, knowledge of the metabolic traits that govern the dynamics of these bloom-forming phytoplankton is limited. We used eukaryotic metatranscriptomic techniques to identify the metabolic basis of functional group-specific traits that may drive the shift between net heterotrophy and autotrophy in the oligotrophic ocean. Replicated blooms were simulated by deep seawater (DSW) addition to mimic nutrient loading in the North Pacific Subtropical Gyre, and the transcriptional responses of phytoplankton functional groups were assayed. Responses of the diatom, haptophyte, and dinoflagellate functional groups in simulated blooms were unique, with diatoms and haptophytes significantly (95% confidence) shifting their quantitative metabolic fingerprint from the in situ condition, whereas dinoflagellates showed little response. Significantly differentially abundant genes identified the importance of colimitation by nutrients, metals, and vitamins in eukaryotic phytoplankton metabolism and bloom formation in this system. The variable transcript allocation ratio, used to quantify transcript reallocation following DSW amendment, differed for diatoms and haptophytes, reflecting the long-standing paradigm of phytoplankton r- and K-type growth strategies. Although the underlying metabolic potential of the large eukaryotic phytoplankton was consistently present, the lack of a bloom during the study period suggests a crucial dependence on physical and biogeochemical forcing, which are susceptible to alteration with changing climate.