Dinoflagellates alter their carbon and nutrient metabolic strategies across environmental gradients in the central Pacific Ocean

Seasonal metabolic dynamics of microeukaryotic plankton: a year-long metatranscriptomic study in a temperate sea

Seasonal fluctuations profoundly affect marine microeukaryotic plankton composition and metabolism, but accurately tracking these changes has been a long-standing challenge. In this study, we present a year-long metatranscriptomic data set from the Southern Bight of the North Sea, shedding light on the seasonal dynamics in temperate plankton ecosystems. We observe distinct shifts in active plankton species and their metabolic processes in response to seasonal changes. We characterized the metabolic signatures of different seasonal phases in detail, thereby revealing the metabolic versatility of dinoflagellates, the heterotrophic dietary strategy of Phaeocystis during its late-stage blooms, and stark variations in summer and fall diatom abundance and metabolic activity across nearby sampling stations. Our data illuminate the varied contributions of microeukaryotic taxa to biomass production and nutrient cycling at different times of the year and allow delineation of their ecological niches.

IMPORTANCE

Ecosystem composition and metabolic functions of temperate marine microeukaryote plankton are strongly influenced by seasonal dynamics. Although monitoring of species composition of microeukaryotes has expanded recently, few methods also contain seasonally resolved information on ecosystem functioning. We generated a year-long spatially resolved metatranscriptomic data set to assess seasonal dynamics of microeukaryote species and their associated metabolic functions in the Southern Bight of the North Sea. Our study underscores the potential of metatranscriptomics as a powerful tool for advancing our understanding of marine ecosystem functionality and resilience in response to environmental changes, emphasizing its potential in continuous marine ecosystem monitoring to enhance our ecological understanding of the ocean's eukaryotic microbiome.

Global biogeography of microbes driving ocean ecological status under climate change

Microbial communities play a crucial role in ocean ecology and global biogeochemical processes. However, understanding the intricate interactions among diversity, taxonomical composition, functional traits, and how these factors respond to climate change remains a significant challenge. Here, we propose seven distinct ecological statuses by systematically considering the diversity, structure, and biogeochemical potential of the ocean microbiome to delineate their biogeography. Anthropogenic climate change is expected to alter the ecological status of the surface ocean by influencing environmental conditions, particularly nutrient and oxygen contents. Our predictive model, which utilizes machine learning, indicates that the ecological status of approximately 32.44% of the surface ocean may undergo changes from the present to the end of this century, assuming no policy interventions. These changes mainly include poleward shifts in the main taxa, increases in photosynthetic carbon fixation and decreases in nutrient metabolism. However, this proportion can decrease significantly with effective control of greenhouse gas emissions. Our study underscores the urgent necessity for implementing policies to mitigate climate change, particularly from an ecological perspective.

New plastids, old proteins: repeated endosymbiotic acquisitions in kareniacean dinoflagellates

Dinoflagellates are a diverse group of ecologically significant micro-eukaryotes that can serve as a model system for plastid symbiogenesis due to their susceptibility to plastid loss and replacement via serial endosymbiosis. Kareniaceae harbor fucoxanthin-pigmented plastids instead of the ancestral peridinin-pigmented ones and support them with a diverse range of nucleus-encoded plastid-targeted proteins originating from the haptophyte endosymbiont, dinoflagellate host, and/or lateral gene transfers (LGT). Here, we present predicted plastid proteomes from seven distantly related kareniaceans in three genera (Karenia, Karlodinium, and Takayama) and analyze their evolutionary patterns using automated tree building and sorting. We project a relatively limited ( ~ 10%) haptophyte signal pointing towards a shared origin in the family Chrysochromulinaceae. Our data establish significant variations in the functional distributions of these signals, emphasizing the importance of micro-evolutionary processes in shaping the chimeric proteomes. Analysis of plastid genome sequences recontextualizes these results by a striking finding the extant kareniacean plastids are in fact not all of the same origin, as two of the studied species (Karlodinium armiger, Takayama helix) possess plastids from different haptophyte orders than the rest.

eDNA metabarcoding reveals biodiversity and depth stratification patterns of dinoflagellate assemblages within the epipelagic zone of the western Coral Sea

BACKGROUND

Dinoflagellates play critical roles in the functioning of marine ecosystems but also may pose a hazard to human and ecosystem health by causing harmful algal blooms (HABs). The Coral Sea is a biodiversity hotspot, but its dinoflagellate assemblages in pelagic waters have not been studied by modern sequencing methods. We used metabarcoding of the 18 S rRNA V4 amplicon to assess the diversity and structure of dinoflagellate assemblages throughout the water column to a depth of 150 m at three stations in the Western Coral Sea. Additionally, at one station we compared metabarcoding with morphological methods to optimise identification and detection of dinoflagellates.

RESULTS

Stratification of dinoflagellate assemblages was evident in depth-specific relative abundances of taxonomic groups; the greatest difference was between the 5-30 m assemblages and the 130-150 m assemblages. The relative abundance of Dinophyceae (photosynthetic and heterotrophic) decreased with increasing depth, whereas that of Syndiniales (parasitic) increased with increasing depth. The composition of major taxonomic groups was similar among stations. Taxonomic richness and diversity of amplicon sequence variants (ASVs) were similar among depths and stations; however, the abundance of dominant taxa was highest within 0-30 m, and the abundance of rare taxa was highest within 130-150 m, indicating adaptations to specific depth strata. The number of unclassified ASVs at the family and species levels was very high, particularly for Syndinian representatives.

CONCLUSIONS

Dinoflagellate assemblages in open water of the Coral Sea are highly diverse and taxonomically stratified by depth; patterns of relative abundance along the depth gradient reflect environmental factors and ecological processes. Metabarcoding detects more species richness than does traditional microscopical methods of sample analysis, yet the methods are complementary, with morphological analysis revealing additional richness. The large number of unclassified dinoflagellate-ASVs indicates a need for improved taxonomic reference databases and suggests presence of dinoflagellate-crypto and-morphospecies.

Insight into planktonic protistan and fungal communities across the nutrient-depleted environment of the South Pacific Subtropical Gyre

Ocean microorganisms constitute ~70% of the marine biomass, contribute to ~50% of the Earth's primary production, and play a vital role in global biogeochemical cycles. The marine heterotrophic and mixotrophic protistan and fungal communities have often been overlooked mainly due to limitations in morphological species identification. Despite the accumulation of studies on biogeographic patterns observed in microbial communities, our understanding of the abundance and distribution patterns within the microbial community of the largest subtropical gyre, the South Pacific Gyre (SPG), remains incomplete. Here, we investigated the diversity and vertical composition of protistan and fungal communities in the water column of the ultra-oligotrophic SPG. Our results showed apparent differences in protistan community diversity in the photic and aphotic regions. The entire protistan community diversity was significantly affected by temperature, salinity, oxygen, and nutrient concentrations, while the parasitic community diversity was also affected by chlorophyll a concentration. The parasitic protists were assigned to the class Syndiniales accounting for over 98% of the total parasitic protists, exhibiting higher relative sequence abundance along the water depth and displaying consistent patterns among different sampling stations. In contrast to the protistan community, the fungal community along the SPG primarily clustered based on the sampling station and pelagic zones. In particular, our study reveals a significant presence of parasitic protists and functionally diverse fungi in SPG and their potential impact on carbon cycling in the gyre.IMPORTANCEOur findings carry important implications for understanding the distribution patterns of the previously unrecognized occurrence of parasitic protists and functionally diverse fungi in the nutrient-limited South Pacific Gyre. In particular, our study reveals a significant presence of parasitic Syndiniales, predominantly abundant in the upper 300 m of the aphotic zone in the gyre, and a distinct presence of fungal communities in the aphotic zone at the central part of the gyre. These findings strongly suggest that these communities play a substantial role in yet insufficiently described microbial food web. Moreover, our research enhances our understanding of their contribution to the dynamics of the food webs in oligotrophic gyres and is valuable for projecting the ecological consequences of future climate warming.

In situ community transcriptomics illuminates CO2-fixation potentials and supporting roles of phagotrophy and proton pump in plankton in a subtropical marginal sea

Lineage-wise physiological activities of plankton communities in the ocean are important but challenging to characterize. Here, we conducted whole-assemblage metatranscriptomic profiling at continental shelf and slope sites in the South China Sea to investigate carbon fixation potential in different lineages. RuBisCO expression, the proxy of Calvin carbon fixation (CCF) potential, was mainly contributed by Bacillariophyta, Chlorophyta, Cyanobacteria, and Haptophyta, which was differentially affected by environmental factors among lineages. CCF potential exhibited positive or negative correlations with phagotrophy gene expression, suggesting phagotrophy possibly enhances or complements CCF. Our data also reveal significant non-Calvin carbon fixation (NCF) potential, as indicated by the active expression of genes in all five currently recognized NCF pathways, mainly contributed by Flavobacteriales, Alteromonadales, and Oceanospirillales. Furthermore, in Flavobacteriales, Alteromonadales, Pelagibacterales, and Rhodobacterales, NCF potential was positively correlated with proton-pump rhodopsin (PPR) expression, suggesting that NCF might be energetically supported by PPR. The novel insights into the lineage-differential potential of carbon fixation, widespread mixotrophy, and PPR as an energy source for NCF lay a methodological and informational foundation for further research to understand carbon fixation and the trophic landscape in the ocean.IMPORTANCEMarine plankton plays an important role in global carbon cycling and climate regulation. Phytoplankton and cyanobacteria fix CO2 to produce organic compounds using solar energy and mainly by the Calvin cycle, whereas autotrophic bacteria and archaea may fix CO2 by non-Calvin cycle carbon fixation pathways. How active individual lineages are in carbon fixation and mixotrophy, and what energy source bacteria may employ in non-Calvin carbon fixation, in a natural plankton assemblage are poorly understood and underexplored. Using metatranscriptomics, we studied carbon fixation in marine plankton with lineage resolution in tropical marginal shelf and slope areas. Based on the sequencing results, we characterized the carbon fixation potential of different lineages and assessed Calvin- and non-Calvin- carbon fixation activities and energy sources. Data revealed a high number of unigenes (4.4 million), lineage-dependent differential potentials of Calvin carbon fixation and responses to environmental conditions, major contributors of non-Calvin carbon fixation, and their potential energy source.

Protistan community composition and metabolism in the North Pacific Subtropical Gyre: Influences of mesoscale eddies and depth

Marine protists and their metabolic activities are intricately tied to the cycling of nutrients and the flow of energy through microbial food webs. Physiochemical changes in the environment, such as those that result from mesoscale eddies, may impact protistan communities, but the effects that such changes have on protists are poorly known. A metatranscriptomic study was conducted to investigate how eddies affected protists at adjacent cyclonic and anticyclonic eddy sites in the oligotrophic ocean at four depths from 25 to 250 m. Eddy polarity impacted protists at all depths sampled, although the effects of eddy polarity were secondary to the impact of depth across the depth range. Eddy-induced vertical shifts in the water column yielded differences in the cyclonic and anticyclonic eddy protistan communities, and these differences were the most pronounced at and just below the deep chlorophyll maximum. An analysis of transcripts associated with protistan nutritional physiology at 150 m revealed that cyclonic eddies may support a more heterotrophic community, while anticyclonic eddies promote a more phototrophic community. The results of this study indicate that eddies alter the metabolism of protists particularly in the lower euphotic zone and may therefore impact carbon export from the euphotic zone.

Alexandrium spp.: From Toxicity to Potential Biotechnological Benefits

Many dinoflagellates of the genus Alexandrium are well known for being responsible for harmful algal blooms (HABs), producing potent toxins that cause damages to other marine organisms, aquaculture, fishery, tourism, as well as induce human intoxications and even death after consumption of contaminated shellfish or fish. In this review, we summarize potential bioprospecting associated to the genus Alexandrium, including which Alexandrium spp. produce metabolites with anticancer, antimicrobial, antiviral, as well as anti-Alzheimer applications. When available, we report their mechanisms of action and targets. We also discuss recent progress on the identification of secondary metabolites with biological properties favorable to human health and aquaculture. Altogether, this information highlights the importance of studying which culturing conditions induce the activation of enzymatic pathways responsible for the synthesis of bioactive metabolites. It also suggests considering and comparing clones collected in different locations for toxin monitoring and marine bioprospecting. This review can be of interest not only for the scientific community, but also for the entire population and industries.

Persistent equatorial Pacific iron limitation under ENSO forcing

Projected responses of ocean net primary productivity to climate change are highly uncertain1. Models suggest that the climate sensitivity of phytoplankton nutrient limitation in the low-latitude Pacific Ocean plays a crucial role1-3, but this is poorly constrained by observations4. Here we show that changes in physical forcing drove coherent fluctuations in the strength of equatorial Pacific iron limitation through multiple El Niño/Southern Oscillation (ENSO) cycles, but that this was overestimated twofold by a state-of-the-art climate model. Our assessment was enabled by first using a combination of field nutrient-addition experiments, proteomics and above-water hyperspectral radiometry to show that phytoplankton physiological responses to iron limitation led to approximately threefold changes in chlorophyll-normalized phytoplankton fluorescence. We then exploited the >18-year satellite fluorescence record to quantify climate-induced nutrient limitation variability. Such synoptic constraints provide a powerful approach for benchmarking the realism of model projections of net primary productivity to climate changes.

Global analysis of ocean phytoplankton nutrient limitation reveals high prevalence of co-limitation

Nutrient availability limits phytoplankton growth throughout much of the global ocean. Here we synthesize available experimental data to identify three dominant nutrient limitation regimes: nitrogen is limiting in the stratified subtropical gyres and in the summertime Arctic Ocean, iron is most commonly limiting in upwelling regions, and both nutrients are frequently co-limiting in regions in between the nitrogen and iron limited systems. Manganese can be co-limiting with iron in parts of the Southern Ocean, whilst phosphate and cobalt can be co-/serially limiting in some settings. Overall, an analysis of experimental responses showed that phytoplankton net growth can be significantly enhanced through increasing the number of different nutrients supplied, regardless of latitude, temperature, or trophic status, implying surface seawaters are often approaching nutrient co-limitation. Assessments of nutrient deficiency based on seawater nutrient concentrations and nutrient stress diagnosed via molecular biomarkers showed good agreement with experimentally-assessed nutrient limitation, validating conceptual and theoretical links between nutrient stoichiometry and microbial ecophysiology.

Phylogeny and biogeography of the algal DMS-releasing enzyme in the global ocean

Phytoplankton produce the volatile dimethyl sulfide (DMS), an important infochemical mediating microbial interactions, which is also emitted to the atmosphere and affecting the global climate. Albeit the enzymatic source for DMS in eukaryotes was elucidated, namely a DMSP lyase (DL) called Alma1, we still lack basic knowledge regarding its taxonomic distribution. We defined unique sequence motifs which enable the identification of DL homologs (DLHs) in model systems and environmental populations. We used these motifs to predict DLHs in diverse algae by analyzing hundreds of genomic and transcriptomic sequences from model systems under stress conditions and from environmental samples. Our findings show that the DL enzyme is more taxonomically widespread than previously thought, as it is encoded by known algal taxa as haptophytes and dinoflagellates, but also by chlorophytes, pelagophytes and diatoms, which were conventionally considered to lack the DL enzyme. By exploring the Tara Oceans database, we showed that DLHs are widespread across the oceans and are predominantly expressed by dinoflagellates. Certain dinoflagellate DLHs were differentially expressed between the euphotic and mesopelagic zones, suggesting a functional specialization and an involvement in the metabolic plasticity of mixotrophic dinoflagellates. In specific regions as the Southern Ocean, DLH expression by haptophytes and diatoms was correlated with environmental drivers such as nutrient availability. The expanded repertoire of putative DL enzymes from diverse microbial origins and geographic niches suggests new potential players in the marine sulfur cycle and provides a foundation to study the cellular function of the DL enzyme in marine microbes.

Mixoplankton and mixotrophy: future research priorities

Phago-mixotrophy, the combination of photoautotrophy and phagotrophy in mixoplankton, organisms that can combine both trophic strategies, have gained increasing attention over the past decade. It is now recognized that a substantial number of protistan plankton species engage in phago-mixotrophy to obtain nutrients for growth and reproduction under a range of environmental conditions. Unfortunately, our current understanding of mixoplankton in aquatic systems significantly lags behind our understanding of zooplankton and phytoplankton, limiting our ability to fully comprehend the role of mixoplankton (and phago-mixotrophy) in the plankton food web and biogeochemical cycling. Here, we put forward five research directions that we believe will lead to major advancement in the field: (i) evolution: understanding mixotrophy in the context of the evolutionary transition from phagotrophy to photoautotrophy; (ii) traits and trade-offs: identifying the key traits and trade-offs constraining mixotrophic metabolisms; (iii) biogeography: large-scale patterns of mixoplankton distribution; (iv) biogeochemistry and trophic transfer: understanding mixoplankton as conduits of nutrients and energy; and (v) in situ methods: improving the identification of in situ mixoplankton and their phago-mixotrophic activity.

Gene expression dynamics of natural assemblages of heterotrophic flagellates during bacterivory

BACKGROUND

Marine heterotrophic flagellates (HF) are dominant bacterivores in the ocean, where they represent the trophic link between bacteria and higher trophic levels and participate in the recycling of inorganic nutrients for regenerated primary production. Studying their activity and function in the ecosystem is challenging since most of the HFs in the ocean are still uncultured. In the present work, we investigated gene expression of natural HF communities during bacterivory in four unamended seawater incubations.

RESULTS

The most abundant species growing in our incubations belonged to the taxonomic groups MAST-4, MAST-7, Chrysophyceae, and Telonemia. Gene expression dynamics were similar between incubations and could be divided into three states based on microbial counts, each state displaying distinct expression patterns. The analysis of samples where HF growth was highest revealed some highly expressed genes that could be related to bacterivory. Using available genomic and transcriptomic references, we identified 25 species growing in our incubations and used those to compare the expression levels of these specific genes. Video Abstract CONCLUSIONS: Our results indicate that several peptidases, together with some glycoside hydrolases and glycosyltransferases, are more expressed in phagotrophic than in phototrophic species, and thus could be used to infer the process of bacterivory in natural assemblages.

Evaluation of Protein Extraction Methods for Metaproteomic Analyses of Root-Associated Microbes

Metaproteomics is a powerful tool for the characterization of metabolism, physiology, and functional interactions in microbial communities, including plant-associated microbiota. However, the metaproteomic methods that have been used to study plant-associated microbiota are very laborious and require large amounts of plant tissue, hindering wider application of these methods. We optimized and evaluated different protein extraction methods for metaproteomics of plant-associated microbiota in two different plant species (Arabidopsis and maize). Our main goal was to identify a method that would work with low amounts of input material (40 to 70 mg) and that would maximize the number of identified microbial proteins. We tested eight protocols, each comprising a different combination of physical lysis method, extraction buffer, and cell-enrichment method on roots from plants grown with synthetic microbial communities. We assessed the performance of the extraction protocols by liquid chromatography-tandem mass spectrometry-based metaproteomics and found that the optimal extraction method differed between the two species. For Arabidopsis roots, protein extraction by beating whole roots with small beads provided the greatest number of identified microbial proteins and improved the identification of proteins from gram-positive bacteria. For maize, vortexing root pieces in the presence of large glass beads yielded the greatest number of microbial proteins identified. Based on these data, we recommend the use of these two methods for metaproteomics with Arabidopsis and maize. Furthermore, detailed descriptions of the eight tested protocols will enable future optimization of protein extraction for metaproteomics in other dicot and monocot plants. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Temperature variability interacts with mean temperature to influence the predictability of microbial phenotypes

Despite their relatively high thermal optima (T_opt ), tropical taxa may be particularly vulnerable to a rising baseline and increased temperature variation because they live in relatively stable temperatures closer to their T_opt . We examined how microbial eukaryotes with differing thermal histories responded to temperature fluctuations of different amplitudes (0 control, ±2, ±4°C) around mean temperatures below or above their T_opt . Cosmopolitan dinoflagellates were selected based on their distinct thermal traits and included two species of the same genus (tropical and temperate Coolia spp.), and two strains of the same species maintained at different temperatures for >500 generations (tropical Amphidinium massartii control temperature and high temperature, CT and HT, respectively). There was a universal decline in population growth rate under temperature fluctuations, but strains with narrower thermal niche breadth (temperate Coolia and HT) showed ~10% greater reduction in growth. At suboptimal mean temperatures, cells in the cool phase of the fluctuation stopped dividing, fixed less carbon (C) and had enlarged cell volumes that scaled positively with elemental C, N, and P and C:Chlorophyll-a. However, at a supra-optimal mean temperature, fixed C was directed away from cell division and novel trait combinations developed, leading to greater phenotypic diversity. At the molecular level, heat-shock proteins, and chaperones, in addition to transcripts involving genome rearrangements, were upregulated in CT and HT during the warm phase of the supra-optimal fluctuation (30 ± 4°C), a stress response indicating protection. In contrast, the tropical Coolia species upregulated major energy pathways in the warm phase of its supra-optimal fluctuation (25 ± 4°C), indicating a broadscale shift in metabolism. Our results demonstrate divergent effects between taxa and that temporal variability in environmental conditions interacts with changes in the thermal mean to mediate microbial responses to global change, with implications for biogeochemical cycling.

The microbiome of the dinoflagellate Prorocentrum cordatum in laboratory culture and its changes at higher temperatures

In the ocean, phytoplankton are dependent on communities of bacteria living in the phycosphere, a hot spot of metabolic and genetic exchange. Many types of interactions between phytoplankton and phycosphere bacteria have been shown, but it is unclear if the microbial communities associated with microalgae strains in culture collections are beneficial or harmful to the host strain. Here, we studied the microbial communities associated with four strains of the dinoflagellate Prorocentrum cordatum that had been isolated from distant geographical locations and maintained in culture collection for hundreds of generations. Community composition was determined by 16S rRNA gene amplicon sequencing. The dinoflagellate host strain was the strongest parameter separating communities, while growth phase, lifestyle (particle-attached versus free-living) and temperature had only a modulating effect. Although the strains had been isolated from distant locations in the Atlantic and Pacific Ocean, 14 ASVs were shared among all strains, the most abundant ones being Gilvibacter, Marivita, uncultivated Rhodobacteraceae, Marinobacter, Hyphomonadaceae, Cupriavidus, Variovorax, and Paucibacter. Adaptation to higher temperatures resulted in specific changes in each phycosphere microbiome, including increased abundance of rare community members. We then compared the growth of the four xenic cultures to that of the axenic P. cordatum CCMP1329. At 20°C, growth of the xenic cultures was similar or slower than that of CCMP1329. At 26°C, all four xenic cultures experienced a death phase, while the axenic culture stably remained in the stationary phase. At 30°C, only two of the xenic cultures were able to grow. A shift of dinoflagellate metabolism from autotrophy to mixotrophy and competition between dinoflagellate and bacteria for limiting nutrients, including essential vitamins, may contribute to these differences in growth patterns.

Microbial functional diversity across biogeochemical provinces in the central Pacific Ocean

Enzymes catalyze key reactions within Earth's life-sustaining biogeochemical cycles. Here, we use metaproteomics to examine the enzymatic capabilities of the microbial community (0.2 to 3 µm) along a 5,000-km-long, 1-km-deep transect in the central Pacific Ocean. Eighty-five percent of total protein abundance was of bacterial origin, with Archaea contributing 1.6%. Over 2,000 functional KEGG Ontology (KO) groups were identified, yet only 25 KO groups contributed over half of the protein abundance, simultaneously indicating abundant key functions and a long tail of diverse functions. Vertical attenuation of individual proteins displayed stratification of nutrient transport, carbon utilization, and environmental stress. The microbial community also varied along horizontal scales, shaped by environmental features specific to the oligotrophic North Pacific Subtropical Gyre, the oxygen-depleted Eastern Tropical North Pacific, and nutrient-rich equatorial upwelling. Some of the most abundant proteins were associated with nitrification and C1 metabolisms, with observed interactions between these pathways. The oxidoreductases nitrite oxidoreductase (NxrAB), nitrite reductase (NirK), ammonia monooxygenase (AmoABC), manganese oxidase (MnxG), formate dehydrogenase (FdoGH and FDH), and carbon monoxide dehydrogenase (CoxLM) displayed distributions indicative of biogeochemical status such as oxidative or nutritional stress, with the potential to be more sensitive than chemical sensors. Enzymes that mediate transformations of atmospheric gases like CO, CO₂, NO, methanethiol, and methylamines were most abundant in the upwelling region. We identified hot spots of biochemical transformation in the central Pacific Ocean, highlighted previously understudied metabolic pathways in the environment, and provided rich empirical data for biogeochemical models critical for forecasting ecosystem response to climate change.

Rapid evolution allows coexistence of highly divergent lineages within the same niche

Marine microbial communities are extremely complex and diverse. The number of locally coexisting species often vastly exceeds the number of identifiable niches, and taxonomic composition often appears decoupled from local environmental conditions. This is contrary to the view that environmental conditions should select for a few locally well-adapted species. Here we use an individual-based eco-evolutionary model to show that virtually unlimited taxonomic diversity can be supported in highly evolving assemblages, even in the absence of niche separation. With a steady stream of heritable changes to phenotype, competitive exclusion may be weakened, allowing sustained coexistence of nearly neutral phenotypes with highly divergent lineages. This behaviour is robust even to abrupt environmental perturbations that might be expected to cause strong selection pressure and an associated loss of diversity. We, therefore, suggest that rapid evolution and individual-level variability are key drivers of species coexistence and maintenance of microbial biodiversity.

A robust approach to estimate relative phytoplankton cell abundances from metagenomes

Phytoplankton account for >45% of global primary production, and have an enormous impact on aquatic food webs and on the entire Earth System. Their members are found among prokaryotes (cyanobacteria) and multiple eukaryotic lineages containing chloroplasts. Genetic surveys of phytoplankton communities generally consist of PCR amplification of bacterial (16S), nuclear (18S) and/or chloroplastic (16S) rRNA marker genes from DNA extracted from environmental samples. However, our appreciation of phytoplankton abundance or biomass is limited by PCR-amplification biases, rRNA gene copy number variations across taxa, and the fact that rRNA genes do not provide insights into metabolic traits such as photosynthesis. Here, we targeted the photosynthetic gene psbO from metagenomes to circumvent these limitations: the method is PCR-free, and the gene is universally and exclusively present in photosynthetic prokaryotes and eukaryotes, mainly in one copy per genome. We applied and validated this new strategy with the size-fractionated marine samples collected by Tara Oceans, and showed improved correlations with flow cytometry and microscopy than when based on rRNA genes. Furthermore, we revealed unexpected features of the ecology of these ecosystems, such as the high abundance of picocyanobacterial aggregates and symbionts in the ocean, and the decrease in relative abundance of phototrophs towards the larger size classes of marine dinoflagellates. To facilitate the incorporation of psbO in molecular-based surveys, we compiled a curated database of >18,000 unique sequences. Overall, psbO appears to be a promising new gene marker for molecular-based evaluations of entire phytoplankton communities.

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.

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.

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.

Online Nanoflow Two-Dimension Comprehensive Active Modulation Reversed Phase-Reversed Phase Liquid Chromatography High-Resolution Mass Spectrometry for Metaproteomics of Environmental and Microbiome Samples

Metaproteomics is a powerful analytical approach that can assess the functional capabilities deployed by microbial communities in both environmental and biomedical microbiome settings. Yet, the mass spectra resulting from these mixed biological communities are challenging to obtain due to the high number of low intensity peak features. The use of multiple dimensions of chromatographic separation prior to mass spectrometry analyses has been applied to proteomics previously but can require increased sampling handling and instrument time. Here, we demonstrate an automated online comprehensive active modulation two-dimensional liquid chromatography method for metaproteome sample analysis. A high pH PLRP-S column was used in the first dimension followed by low pH separation in the second dimension using dual modulating C18 traps and a C18 column. This method increased the number of unique peptides found in ocean metaproteome samples by more than 50% when compared to a one-dimension separation while using the same amount of sample and instrument time.