Patterns of diversity in marine phytoplankton

Environmental controls on the seasonal variations of diatoms and dinoflagellates in the Qingdao coastal region, the Yellow Sea

Diatoms and dinoflagellates are two typical functional groups of phytoplankton assemblages, which play a crucial role in the structure and functioning of most marine ecosystems. To date, a novel challenge in ecology and biogeochemistry is to address the influences of environmental changes associated with climate change and human activities on the dynamics of diatoms and dinoflagellates. However, the knowledge of the key environmental factors controlling the diatom-dinoflagellate dynamics remains to be improved, particularly in the coastal ecosystems. Therefore, we conducted four cruises along the Qingdao coastline in spring, summer, autumn, and winter 2022 to explore how diatoms and dinoflagellates varied in response to regional environmental changes. The results showed that the phytoplankton communities were dominated by diatoms and dinoflagellates in terms of abundance and species diversity throughout the year in the study region. Yet, there were significant seasonal variability of diatoms and dinoflagellates across the four seasons. For example, diatom species was the most diverse during autumn, and the higher average abundance was observed in the fall and winter. In contrast, the average abundance of dinoflagellates was maximum during the summer and minimum in the autumn season. Moreover, the abundance and species ratios of diatoms/dinoflagellates (dia/dino) also showed significant seasonal variations in the region. The dia/dino abundance ratio was lowest in summer, while the dia/dino species ratio showed an increasing trend from spring to fall and a slight descending trend during winter. Based on the redundancy analysis, we revealed that diatoms and dinoflagellates responded differently to various environmental variables in different seasons, of which temperature and nutrients (especially dissolved inorganic nitrogen, DIN) had highly significant correlations with both the dia/dino abundance and species ratios. Thus, we suggested that temperature and DIN were the key factors controlling the seasonal dynamics of diatoms and dinoflagellates in the Qingdao coastal area.

Neutral Theory and Plankton Biodiversity

The biodiversity of the plankton has been interpreted largely through the monocle of competition. The spatial distancing of phytoplankton in nature is so large that cell boundary layers rarely overlap, undermining opportunities for resource-based competitive exclusion. Neutral theory accounts for biodiversity patterns based purely on random birth, death, immigration, and speciation events and has commonly served as a null hypothesis in terrestrial ecology but has received comparatively little attention in aquatic ecology. This review summarizes basic elements of neutral theory and explores its stand-alone utility for understanding phytoplankton diversity. A theoretical framework is described entailing a very nonneutral trophic exclusion principle melded with the concept of ecologically defined neutral niches. This perspective permits all phytoplankton size classes to coexist at any limiting resource level, predicts greater diversity than anticipated from readily identifiable environmental niches but less diversity than expected from pure neutral theory, and functions effectively in populations of distantly spaced individuals.

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.

Trichocoleus desertorum isolated from Negev desert petroglyphs: Characterization, adaptation and bioerosion potential

The Negev petroglyphs are considered valuable cultural heritage sites, but unfortunately, they are exposed to deterioration processes caused by anthropogenic and natural forces. Despite the many studies that have already pointed to the role of cyanobacteria in biogenic rock weathering, the knowledge involved in the process is still lacking. In this study, a cyanobacterial strain was isolated from the Negev Desert petroglyphs aiming to reveal its involvement in geochemical cycles and in the weathering processes of the rock substrate. The strain was characterized using morphological, molecular, and microscopic studies. The morphological research revealed a green-bluish, bundle-forming filamentous strain characterized by trichomes embedded in a common sheath. A combination of Nanopore and Illumina sequencing technologies facilitated the assembly of a near-complete genome containing 5,458,034 base pairs. A total of 5027 coding sequences were revealed by implementing PROKKA software. Annotation of five replicas of the 16S ribosomal RNA genes revealed that the Negev cyanobacteria isolate is closely (99.73 %) related to Trichocoleus desertorum LSB90_MW403957 isolated from the Sahara Desert, Algeria. The local strain was thus named Trichocoleus desertorum NBK24 CP116619. Several gene sequences that code for possible environmental adaptation mechanisms were found. Our study also identified genes for membrane transporters involved in the exchange of chemical elements, suggesting a possible interaction with rock minerals. Microscopic observations of T. desertorum NBK24 CP116619 infected onto calcareous stone slabs under laboratory conditions demonstrated the effect of the isolated cyanobacteria on stone surface degradation. In conclusion, the findings of this study further our understanding of terrestrial cyanobacterial genomes and functions and highlight the role of T. desertorum NBK24 CP116619 in stone weathering processes. This information may contribute to the creation of efficient restoration strategies for stone monuments affected by cyanobacteria.

Environmental forcing of phytoplankton carbon-to-diversity ratio and carbon-to-chlorophyll ratio: A case study in Jiaozhou Bay, the Yellow Sea

The relationships between phytoplankton carbon (C) biomass and diversity (i.e., C-to-H' ratio) and chlorophyll a (i.e., C-to-Chl a ratio) are good indicators of marine ecosystem functioning and stability. Here we conducted four cruises spanning 2 years in Jiaozhou Bay to explore the dynamics of C-to-H' and C-to-Chl a ratios. The results showed that the phytoplankton C biomass and diversity were dominated by diatoms, followed by dinoflagellates. The average C-to-H' ratio ranged from 84.10 to 912.17, with high values occurring in the northern region of the bay. In contrast, the average C-to-Chl a ratio ranged between 15.55 and 89.47, and high values primarily appeared in the northern or northeastern part of the bay. In addition, the redundancy analysis showed that temperature and phosphate (DIP) were significantly correlated with both ratios in most cases, indicating that temperature and DIP may be key factors affecting the dynamics of C-to-H' and C-to-Chl a ratios.

A Gulf Stream frontal eddy harbors a distinct microbiome compared to adjacent waters

Mesoscale oceanographic features, including eddies, have the potential to alter productivity and other biogeochemical rates in the ocean. Here, we examine the microbiome of a cyclonic, Gulf Stream frontal eddy, with a distinct origin and environmental parameters compared to surrounding waters, in order to better understand the processes dominating microbial community assembly in the dynamic coastal ocean. Our microbiome-based approach identified the eddy as distinct from the surround Gulf Stream waters. The eddy-associated microbial community occupied a larger area than identified by temperature and salinity alone, increasing the predicted extent of eddy-associated biogeochemical processes. While the eddy formed on the continental shelf, after two weeks both environmental parameters and microbiome composition of the eddy were most similar to the Gulf Stream, suggesting the effect of environmental filtering on community assembly or physical mixing with adjacent Gulf Stream waters. In spite of the potential for eddy-driven upwelling to introduce nutrients and stimulate primary production, eddy surface waters exhibit lower chlorophyll a along with a distinct and less even microbial community, compared to the Gulf Stream. At the population level, the eddy microbiome exhibited differences among the cyanobacteria (e.g. lower Trichodesmium and higher Prochlorococcus) and in the heterotrophic alpha Proteobacteria (e.g. lower relative abundances of specific SAR11 phylotypes) versus the Gulf Stream. However, better delineation of the relative roles of processes driving eddy community assembly will likely require following the eddy and surrounding waters since inception. Additionally, sampling throughout the water column could better clarify the contribution of these mesoscale features to primary production and carbon export in the oceans.

Periodic temporal environmental variations induce coexistence in resource competition models

Natural ecosystems, in particular on the microbial scale, are inhabited by a large number of species. The population size of each species is affected by interactions of individuals with each other and by spatial and temporal changes in environmental conditions, such as resource abundance. Here, we use a generic population dynamics model to study how, and under what conditions, a periodic temporal environmental variation can alter an ecosystem's composition and biodiversity. We demonstrate that using timescale separation allows one to qualitatively predict the long-term population dynamics of interacting species in varying environments. We show that the notion of Tilman's R* rule, a well-known principle that applies for constant environments, can be extended to periodically varying environments if the timescale of environmental changes (e.g., seasonal variations) is much faster than the timescale of population growth (doubling time in bacteria). When these timescales are similar, our analysis shows that a varying environment deters the system from reaching a steady state, and stable coexistence between multiple species becomes possible. Our results posit that biodiversity can in part be attributed to natural environmental variations.

Fungal Diversity and Community Composition across Ecosystems

Fungi have shaped the biosphere since the development of life on Earth. Despite fungi being present in all environments, most of the available fungal research has focused on soils. As a result, the role and composition of fungal communities in aquatic (marine and freshwater) environments remain largely unexplored. The use of different primers to characterise fungal communities has additionally complicated intercomparisons among studies. Consequently, we lack a basic global assessment of fungal diversity across major ecosystems. Here, we took advantage of a recently published 18S rRNA dataset comprising samples from major ecosystems (terrestrial, freshwater, and marine) to attempt a global assessment of fungal diversity and community composition. We found the highest fungal diversities for terrestrial > freshwater > marine environments, and pronounced gradients of fungal diversity along temperature, salinity, and latitude in all ecosystems. We also identified the most abundant taxa in each of these ecosystems, mostly dominated by Ascomycota and Basidiomycota, except in freshwater rivers where Chytridiomycota dominated. Collectively, our analysis provides a global analysis of fungal diversity across all major environmental ecosystems, highlighting the most distinct order and ASVs (amplicon sequencing variants) by ecosystem, and thus filling a critical gap in the study of the Earth's mycobiome.

Modeling present and future distribution of plankton populations in a coastal upwelling zone: the copepod Calanus chilensis as a study case

Predicting species distribution in the ocean has become a crucial task to assess marine ecosystem responses to ongoing climate change. In the Humboldt Current System (HCS), the endemic copepod Calanus chilensis is one of the key species bioindicator of productivity and water masses. Here we modeled the geographic distribution of Calanus chilensis for two bathymetric ranges, 0-200 and 200-400 m. For the 0-200 m layer, we used the Bayesian Additive Regression Trees (BART) method, whereas, for the 200-400 m layer, we used the Ensembles of Small Models (ESMs) method and then projected the models into two future scenarios to assess changes in geographic distribution patterns. The models were evaluated using the multi-metric approach. We identified that chlorophyll-a (0.34), Mixed Layer Depth (0.302) and salinity (0.36) explained the distribution of C. chilensis. The geographic prediction of the BART model revealed a continuous distribution from Ecuador to the southernmost area of South America for the 0-200 m depth range, whereas the ESM model indicated a discontinuous distribution with greater suitability for the coast of Chile for the 200-400 m depth range. A reduction of the distribution range of C. chilensis is projected in the future. Our study suggests that the distribution of C. chilensis is conditioned by productivity and mesoscale processes, with both processes closely related to upwelling intensity. These models serve as a tool for proposing indicators of changes in the ocean. We further propose that the species C. chilensis is a high productivity and low salinity indicator at the HCS. We recommend further examining multiple spatial and temporal scales for stronger inference.

Phytoplankton Community Response to Environmental Factors along a Salinity Gradient in a Seagoing River, Tianjin, China

A river-estuary ecosystem usually features a distinct salinity gradient and a complex water environment, so it is enormously valuable to study the response mechanism of living organisms to multiple abiotic factors under salinity stress. Phytoplankton, as an important part of aquatic microorganisms, has always been of concern for its crucial place in the aquatic ecosystem. In this study, phytoplankton data and 18 abiotic factors collected from 15 stations in Duliujian River, a seagoing river, were investigated in different seasons. The results showed that the river studied was of a Cyanophyta-dominant type. Salinity (SAL) was the key control factor for phytoplankton species richness, while water temperature (WT) was critical not only for species richness, but also community diversity, and the abundance and biomass of dominant species. Apart from WT, the abundance and biomass of dominant species were also driven by total nitrogen (TN), nitrate (NO₃^-), pH, and water transparency (SD). Moreover, total dissolved phosphorus (TDP), pH, and chemical oxygen demand (COD) were crucial for community diversity and evenness. The bloom of dominant species positively associated with TDP led to lower diversity and evenness in autumn. In addition, when available nitrogen was limited, Pseudoanabaena sp. could obtain a competitive advantage through the N₂ fixation function. Increased available nitrogen concentration could favor the abundance of Chlorella vulgaris to resist the negative effect of WT. The results show that Oscillatoria limosa could serve as an indicator of organic contamination, and nutrient-concentration control must be effective to inhibit Microcystis bloom. This could help managers to formulate conservation measures.

Fine-scale diversity of microbial communities due to satellite niches in boom and bust environments

Recent observations have revealed that closely related strains of the same microbial species can stably coexist in natural and laboratory settings subject to boom and bust dynamics and serial dilutions, respectively. However, the possible mechanisms enabling the coexistence of only a handful of strains, but not more, have thus far remained unknown. Here, using a consumer-resource model of microbial ecosystems, we propose that by differentiating along Monod parameters characterizing microbial growth rates in high and low nutrient conditions, strains can coexist in patterns similar to those observed. In our model, boom and bust environments create satellite niches due to resource concentrations varying in time. These satellite niches can be occupied by closely related strains, thereby enabling their coexistence. We demonstrate that this result is valid even in complex environments consisting of multiple resources and species. In these complex communities, each species partitions resources differently and creates separate sets of satellite niches for their own strains. While there is no theoretical limit to the number of coexisting strains, in our simulations, we always find between 1 and 3 strains coexisting, consistent with known experiments and observations.

Whole-genome scanning reveals environmental selection mechanisms that shape diversity in populations of the epipelagic diatom Chaetoceros

Diatoms form a diverse and abundant group of photosynthetic protists that are essential players in marine ecosystems. However, the microevolutionary structure of their populations remains poorly understood, particularly in polar regions. Exploring how closely related diatoms adapt to different environments is essential given their short generation times, which may allow rapid adaptations, and their prevalence in marine regions dramatically impacted by climate change, such as the Arctic and Southern Oceans. Here, we address genetic diversity patterns in Chaetoceros, the most abundant diatom genus and one of the most diverse, using 11 metagenome-assembled genomes (MAGs) reconstructed from Tara Oceans metagenomes. Genome-resolved metagenomics on these MAGs confirmed a prevalent distribution of Chaetoceros in the Arctic Ocean with lower dispersal in the Pacific and Southern Oceans as well as in the Mediterranean Sea. Single-nucleotide variants identified within the different MAG populations allowed us to draw a landscape of Chaetoceros genetic diversity and revealed an elevated genetic structure in some Arctic Ocean populations. Gene flow patterns of closely related Chaetoceros populations seemed to correlate with distinct abiotic factors rather than with geographic distance. We found clear positive selection of genes involved in nutrient availability responses, in particular for iron (e.g., ISIP2a, flavodoxin), silicate, and phosphate (e.g., polyamine synthase), that were further supported by analysis of Chaetoceros transcriptomes. Altogether, these results highlight the importance of environmental selection in shaping diatom diversity patterns and provide new insights into their metapopulation genomics through the integration of metagenomic and environmental data.

Differentiated Evolutionary Strategies of Genetic Diversification in Atlantic and Pacific Thaumarchaeal Populations

Some marine microbes are seemingly “ubiquitous,” thriving across a wide range of environmental conditions. While the increased depth in metagenomic sequencing has led to a growing body of research on within-population heterogeneity in environmental microbial populations, there have been fewer systematic comparisons and characterizations of population-level genetic diversity over broader expanses of time and space. Here, we investigated the factors that govern the diversification of ubiquitous microbial taxa found within and between ocean basins. Specifically, we use mapped metagenomic paired reads to examine the genetic diversity of ammonia-oxidizing archaeal (“Candidatus Nitrosopelagicus brevis”) populations in the Pacific (Hawaii Ocean Time-series [HOT]) and Atlantic (Bermuda Atlantic Time Series [BATS]) Oceans sampled over 2 years. We observed higher nucleotide diversity in “Ca. N. brevis” at HOT, driven by a higher rate of homologous recombination. In contrast, “Ca. N. brevis” at BATS featured a more open pangenome with a larger set of genes that were specific to BATS, suggesting a history of dynamic gene gain and loss events. Furthermore, we identified highly differentiated genes that were regulatory in function, some of which exhibited evidence of recent selective sweeps. These findings indicate that different modes of genetic diversification likely incur specific adaptive advantages depending on the selective pressures that they are under. Within-population diversity generated by the environment-specific strategies of genetic diversification is likely key to the ecological success of “Ca. N. brevis.”

IMPORTANCE Ammonia-oxidizing archaea (AOA) are one of the most abundant chemolithoautotrophic microbes in the marine water column and are major contributors to global carbon and nitrogen cycling. Despite their ecological importance and geographical pervasiveness, there have been limited systematic comparisons and characterizations of their population-level genetic diversity over time and space. Here, we use metagenomic time series from two ocean observatories to address the fundamental questions of how abiotic and biotic factors shape the population-level genetic diversity and how natural microbial populations adapt across diverse habitats. We show that the marine AOA “Candidatus Nitrosopelagicus brevis” in different ocean basins exhibits distinct modes of genetic diversification in response to their selective regimes shaped by nutrient availability and patterns of environmental fluctuations. Our findings specific to “Ca. N. brevis” have broader implications, particularly in understanding the population-level responses to the changing climate and predicting its impact on biogeochemical cycles.

Influence of nutrient supply on plankton microbiome biodiversity and distribution in a coastal upwelling region

The ecological and oceanographic processes that drive the response of pelagic ocean microbiomes to environmental changes remain poorly understood, particularly in coastal upwelling ecosystems. Here we show that seasonal and interannual variability in coastal upwelling predicts pelagic ocean microbiome diversity and community structure in the Southern California Current region. Ribosomal RNA gene sequencing, targeting prokaryotic and eukaryotic microbes, from samples collected seasonally during 2014-2020 indicate that nitracline depth is the most robust predictor of spatial microbial community structure and biodiversity in this region. Striking ecological changes occurred due to the transition from a warm anomaly during 2014-2016, characterized by intense stratification, to cooler conditions in 2017-2018, representative of more typical upwelling conditions, with photosynthetic eukaryotes, especially diatoms, changing most strongly. The regional slope of nitracline depth exerts strong control on the relative proportion of highly diverse offshore communities and low biodiversity, but highly productive nearshore communities.

Coastal upwelling sustains some of the most productive ocean regions. Here, the authors find that spatial patterns and temporal changes in nutrient supply explain marine microbial community structure and diversity in the Southern California Current region.

Highly-resolved interannual phytoplankton community dynamics of the coastal Northwest Atlantic

Microbial observatories can track phytoplankton at frequencies that resolve monthly, seasonal, and multiyear trends in environmental change from short-lived events. Using 4-years of weekly flow cytometry along with chloroplast and cyanobacterial 16S rRNA gene sequence data from a time-series station in the coastal Northwest Atlantic (Bedford Basin, Nova Scotia, Canada), we analyzed temporal observations for globally-relevant genera (e.g., Bolidomonas, Teleaulax, Minidiscus, Chaetoceros, Synechococcus, and Phaeocystis) in an oceanic region that has been recognized as a likely hotspot for phytoplankton diversity. Contemporaneous Scotian Shelf data also collected during our study established that the major phytoplankton within the Bedford Basin were important in the Scotian Shelf during spring and fall, therefore pointing to their broader significance within the coastal Northwest Atlantic (NWA). Temporal trends revealed a subset of indicator taxa along with their DNA signatures (e.g., Eutreptiella and Synechococcus), whose distribution patterns make them essential for timely detection of environmentally-driven shifts in the NWA. High-resolution sampling was key to identifying important community shifts towards smaller phytoplankton under anomalous environmental conditions, while further providing a detailed molecular view of community compositions underpinning general phytoplankton succession within the coastal NWA. Our study demonstrates the importance of accessible coastal time-series sites where high-frequency DNA sampling allows for the detection of shifting baselines in phytoplankton communities.

Functional Traits Resolve Mechanisms Governing the Assembly and Distribution of Nitrogen-Cycling Microbial Communities in the Global Ocean

Microorganisms drive much of the marine nitrogen (N) cycle, which jointly controls the primary production in the global ocean. However, our understanding of the microbial communities driving the global ocean N cycle remains fragmented. Focusing on “who is doing what, where, and how?”, this study draws a clear picture describing the global biogeography of marine N-cycling microbial communities by utilizing the Tara Oceans shotgun metagenomes. The marine N-cycling communities are highly variable taxonomically but relatively even at the functional trait level, showing clear functional redundancy properties. The functional traits and taxonomic groups are shaped by the same set of geo-environmental factors, among which, depth is the major factor impacting marine N-cycling communities, differentiating mesopelagic from epipelagic communities. Latitudinal diversity gradients and distance-decay relationships are observed for taxonomic groups, but rarely or weakly for functional traits. The composition of functional traits is strongly deterministic as revealed by null model analysis, while a higher degree of stochasticity is observed for taxonomic composition. Integrating multiple lines of evidence, in addition to drawing a biogeographic picture of marine N-cycling communities, this study also demonstrated an essential microbial ecological theory—determinism governs the assembly of microbial communities performing essential biogeochemical processes; the environment selects functional traits rather than taxonomic groups; functional redundancy underlies stochastic taxonomic community assembly.

Biodiversity exploration in autumn using environmental DNA in the South China sea

The South China Sea (SCS) is an important part of the Indo-Pacific convergence zone, with high biodiversity and abundant marine resources. Traditional methods are primarily used to monitor biodiversity. However, a few studies have used environmental DNA (eDNA) metabarcoding to research the assemblage structure of the SCS. This study used eDNA metabarcoding to survey the SCS assemblage and its relationship with environmental factors over a month-long time-series (August 30th to September 30th, 2020) of seawater samples from the central part of the SCS (9°-20°86' N, 113°-118°47' E). 32 stations were divided into six groups (A, B, C, D, E, F) according to longitude. We collected water samples, extracted eDNA, and amplified 18S rRNA gene V4 region (18S V4), 18S rRNA gene V9 region (18S V9), and 12S rRNA gene (12S). Krona diagrams were used to show species composition. We identified 192 phytoplankton, 104 invertebrate, and 61 fish species from 18S V4, 18S V9, and 12S, respectively. Generally, the three assemblage structures exhibited an increase in species diversity with increasing longitude. Group E had the highest fish diversity. Groups F and C had the highest phytoplankton and invertebrate diversity, respectively. Canonical correspondence analysis showed that four factors (chlorophyll a, depth, salinity, and temperature) were correlated with assemblage structure. Chlorophyll a was the main environmental factor that affected fish, phytoplankton, and invertebrate assemblage structures; salinity was strongly correlated with fish and invertebrate assemblage structures; temperature was a key factor that impacted fish and invertebrate assemblage structures; and depth was strongly correlated with invertebrate assemblage structure. Our results revealed that eDNA metabarcoding is a powerful tool for improving detection rate and using multiple markers is an effective approach for monitoring biodiversity. This study provided information that can be used to enhance biodiversity protection efforts in the SCS.

Prochlorococcus Exudate Stimulates Heterotrophic Bacterial Competition with Rival Phytoplankton for Available Nitrogen

The marine cyanobacterium Prochlorococcus numerically dominates the phytoplankton community of the nutrient-limited open ocean, establishing itself as the most abundant photosynthetic organism on Earth. This ecological success has been attributed to lower cell quotas for limiting nutrients, superior resource acquisition, and other advantages associated with cell size reduction and genome streamlining. In this study, we tested the prediction that Prochlorococcus outcompetes its rivals for scarce nutrients and that this advantage leads to its numerical success in nutrient-limited waters. Strains of Prochlorococcus and its sister genus Synechococcus grew well in both mono- and cocultures when nutrients were replete. However, in nitrogen-limited medium, Prochlorococcus outgrew Synechococcus but only when heterotrophic bacteria were also present. In the nitrogen-limited medium, the heterotroph Alteromonas macleodii outcompeted Synechococcus for nitrogen but only if stimulated by the exudate released by Prochlorococcus or if a proxy organic carbon source was provided. Genetic analysis of Alteromonas suggested that it outcompetes Synechococcus for nitrate and/or nitrite, during which cocultured Prochlorococcus grows on ammonia or other available nitrogen species. We propose that Prochlorococcus can stimulate antagonism between heterotrophic bacteria and potential phytoplankton competitors through a metabolic cross-feeding interaction, and this stimulation could contribute to the numerical success of Prochlorococcus in nutrient-limited regions of the ocean. IMPORTANCE In nutrient-poor habitats, competition for limited resources is thought to select for organisms with an enhanced ability to scavenge nutrients and utilize them efficiently. Such adaptations characterize the cyanobacterium Prochlorococcus, the most abundant photosynthetic organism in the nutrient-limited open ocean. In this study, the competitive superiority of Prochlorococcus over a rival cyanobacterium, Synechococcus, was captured in laboratory culture. Critically, this outcome was achieved only when key aspects of the open ocean were simulated: a limited supply of nitrogen and the presence of heterotrophic bacteria. The results indicate that Prochlorococcus promotes its numerical dominance over Synechococcus by energizing the heterotroph's ability to outcompete Synechococcus for available nitrogen. This study demonstrates how interactions between trophic groups can influence interactions within trophic groups and how these interactions likely contribute to the success of the most abundant photosynthetic microorganism.

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.

Higher-order effects, continuous species interactions, and trait evolution shape microbial spatial dynamics

Persistently diverse microbial communities are one of biology’s great puzzles. Using a modeling framework that accommodates high mutation rates and a continuum of species traits, we studied microbial communities in which antagonistic interactions occur via the production of, inhibition of, and vulnerability to toxins (e.g., antibiotics). Mutation size and mobility enhanced microbial diversity and temporal persistence to extraordinarily high levels. These findings—including the discovery that the duration of the transient phase in community assembly provides a guide to equilibrial diversity—highlight the potentially critical role that antagonistic interactions play in promoting the diversity of bacterial systems. Such interactions, together with resource-driven interactions and spatial structure, may drive the enigmatic levels of biodiversity seen in microbial systems.

The assembly and maintenance of microbial diversity in natural communities, despite the abundance of toxin-based antagonistic interactions, presents major challenges for biological understanding. A common framework for investigating such antagonistic interactions involves cyclic dominance games with pairwise interactions. The incorporation of higher-order interactions in such models permits increased levels of microbial diversity, especially in communities in which antibiotic-producing, sensitive, and resistant strains coexist. However, most such models involve a small number of discrete species, assume a notion of pure cyclic dominance, and focus on low mutation rate regimes, none of which well represent the highly interlinked, quickly evolving, and continuous nature of microbial phenotypic space. Here, we present an alternative vision of spatial dynamics for microbial communities based on antagonistic interactions—one in which a large number of species interact in continuous phenotypic space, are capable of rapid mutation, and engage in both direct and higher-order interactions mediated by production of and resistance to antibiotics. Focusing on toxin production, vulnerability, and inhibition among species, we observe highly divergent patterns of diversity and spatial community dynamics. We find that species interaction constraints (rather than mobility) best predict spatiotemporal disturbance regimes, whereas community formation time, mobility, and mutation size best explain patterns of diversity. We also report an intriguing relationship among community formation time, spatial disturbance regimes, and diversity dynamics. This relationship, which suggests that both higher-order interactions and rapid evolution are critical for the origin and maintenance of microbial diversity, has broad-ranging links to the maintenance of diversity in other systems.

Bottom-Heavy Trophic Pyramids Impair Methylmercury Biomagnification in the Marine Plankton Ecosystems

Methylmercury (CH₃Hg⁺, MMHg) in the phytoplankton and zooplankton, which form the bottom of marine food webs, is a good predictor of MMHg in top predators, including humans. Therefore, evaluating the potential exposure of MMHg to higher trophic levels (TLs) requires a better understanding of relationships between MMHg biomagnification and plankton dynamics. In this study, a coupled ecological/physical model with 366 plankton types of different sizes, biogeochemical functions, and temperature tolerance is used to simulate the relationships between MMHg biomagnification and the ecosystem structure. The study shows that the MMHg biomagnification becomes more significant with increasing TLs. Trophic magnification factors (TMFs) in the lowest two TLs show the opposite spatial pattern to TMFs in higher TLs. The low TMFs are usually associated with a short food-chain length. The less bottom-heavy trophic pyramids in the oligotrophic oceans enhance the MMHg trophic transfer. The global average TMF is increased from 2.3 to 2.8 in the warmer future with a medium climate sensitivity of 2.5 °C. Our study suggests that if there are no mitigation measures for Hg emission, MMHg in the high-trophic-level plankton is increased more dramatically in the warming future, indicating greater MMHg exposure for top predators such as humans.

Abrupt shifts in 21st-century plankton communities

Marine microbial communities sustain ocean food webs and mediate global elemental cycles. These communities will change with climate; these changes can be gradual or foreseeable but likely have much more substantial consequences when sudden and unpredictable. In a complex virtual marine microbial ecosystem, we find that climate change–driven shifts over the 21st century are often abrupt, large in amplitude and extent, and unpredictable using standard early warning signals. Phytoplankton with unique resource needs, especially fast-growing species such as diatoms, are more prone to abrupt shifts. Abrupt shifts in biomass, productivity, and community structure are concentrated in Atlantic and Pacific subtropics. Abrupt changes in environmental variables such as temperature and nutrients rarely precede these ecosystem shifts, indicating that rapid community restructuring can occur in response to gradual environmental changes, particularly in nutrient supply rate ratios.

Can Top-Down Controls Expand the Ecological Niche of Marine N2 Fixers?

The ability of marine diazotrophs to fix dinitrogen gas (N₂) is one of the most influential yet enigmatic processes in the ocean. With their activity diazotrophs support biological production by fixing about 100–200 Tg N/year and turning otherwise unavailable dinitrogen into bioavailable nitrogen (N), an essential limiting nutrient. Despite their important role, the factors that control the distribution of diazotrophs and their ability to fix N₂ are not fully elucidated. We discuss insights that can be gained from the emerging picture of a wide geographical distribution of marine diazotrophs and provide a critical assessment of environmental (bottom-up) versus trophic (top-down) controls. We expand a simplified theoretical framework to understand how top-down control affects competition for resources that determine ecological niches. Selective mortality, mediated by grazing or viral-lysis, on non-fixing phytoplankton is identified as a critical process that can broaden the ability of diazotrophs to compete for resources in top-down controlled systems and explain an expanded ecological niche for diazotrophs. Our simplified analysis predicts a larger importance of top-down control on competition patterns as resource levels increase. As grazing controls the faster growing phytoplankton, coexistence of the slower growing diazotrophs can be established. However, these predictions require corroboration by experimental and field data, together with the identification of specific traits of organisms and associated trade-offs related to selective top-down control. Elucidation of these factors could greatly improve our predictive capability for patterns and rates of marine N₂ fixation. The susceptibility of this key biogeochemical process to future changes may not only be determined by changes in environmental conditions but also via changes in the ecological interactions.

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.

Truncated bimodal latitudinal diversity gradient in early Paleozoic phytoplankton

The latitudinal diversity gradient (LDG)-the decline in species richness from the equator to the poles-is classically considered as the most pervasive macroecological pattern on Earth, but the timing of its establishment, its ubiquity in the geological past, and explanatory mechanisms remain uncertain. By combining empirical and modeling approaches, we show that the first representatives of marine phytoplankton exhibited an LDG from the beginning of the Cambrian, when most major phyla appeared. However, this LDG showed a single peak of diversity centered on the Southern Hemisphere, in contrast to the equatorial peak classically observed for most modern taxa. We find that this LDG most likely corresponds to a truncated bimodal gradient, which probably results from an uneven sediment preservation, smaller sampling effort, and/or lower initial diversity in the Northern Hemisphere. Variation of the documented LDG through time resulted primarily from fluctuations in annual sea-surface temperature and long-term climate changes.

Depth-Dependent Variables Shape Community Structure and Functionality in the Prince Edward Islands

Physicochemical variables limit and control the distribution of microbial communities in all environments. In the oceans, this may significantly influence functional processes such the consumption of dissolved organic material and nutrient sequestration. Yet, the relative contributions of physical factors, such as water mass variability and depth, on functional processes are underexplored. We assessed microbial community structure and functionality in the Prince Edward Islands (PEIs) using 16S rRNA gene amplicon analysis and extracellular enzymatic activity assays, respectively. We found that depth and nutrients substantially drive the structural patterns of bacteria and archaea in this region. Shifts from epipelagic to bathypelagic zones were linked to decreases in the activities of several extracellular enzymes. These extracellular enzymatic activities were positively correlated with several phyla including several Alphaproteobacteria (including members of the SAR 11 clade and order Rhodospirillales) and Cyanobacteria. We show that depth-dependent variables may be essential drivers of community structure and functionality in the PEIs.

Phytoplankton taxonomic and functional diversity patterns across a coastal tidal front

Oceanic physics at fine scale; e.g. eddies, fronts, filaments; are notoriously difficult to sample. However, an increasing number of theoretical approaches hypothesize that these processes affect phytoplankton diversity which have cascading effects on regional ecosystems. In 2015, we targeted the Iroise Sea (France) and evidenced the setting up of the Ushant tidal front from the beginning of spring to late summer. Seawater samples were taken during three sampling cruises and DNA-barcoding allowed us to investigate patterns of eukaryotic phytoplankton diversity across this front. First focusing on patterns of taxonomic richness, we evidenced that the front harbored a hotspot of eukaryotic phytoplankton diversity sustained throughout summer. We then detail the ecological processes leading to the formation of this hotspot by studying shifts in community composition across the Iroise Sea. Physical mixing mingled the communities surrounding the front, allowing the formation of a local ecotone, but it was cycles of disturbances and nutrient inputs over the front that allowed a decrease in competitive exclusion, which maintained a higher diversity of rare phytoplankton taxa. These processes did not select a specific ecological strategy as inferred by a trait approach coupled to our taxonomic approach. Instead the front favored higher richness within widespread strategies, resulting in functional redundancy. We detail how fine-scale ocean physics affect phytoplankton diversity and suppose that this interplay is a major control on regional ecosystems.

Microbial tropicalization driven by a strengthening western ocean boundary current

Western boundary currents (WBCs) redistribute heat and oligotrophic seawater from the tropics to temperate latitudes, with several displaying substantial climate change-driven intensification over the last century. Strengthening WBCs have been implicated in the poleward range expansion of marine macroflora and fauna, however, the impacts on the structure and function of temperate microbial communities are largely unknown. Here we show that the major subtropical WBC of the South Pacific Ocean, the East Australian Current (EAC), transports microbial assemblages that maintain tropical and oligotrophic (k-strategist) signatures, to seasonally displace more copiotrophic (r-strategist) temperate microbial populations within temperate latitudes of the Tasman Sea. We identified specific characteristics of EAC microbial assemblages compared with non-EAC assemblages, including strain transitions within the SAR11 clade, enrichment of Prochlorococcus, predicted smaller genome sizes and shifts in the importance of several functional genes, including those associated with cyanobacterial photosynthesis, secondary metabolism and fatty acid and lipid transport. At a temperate time-series site in the Tasman Sea, we observed significant reductions in standing stocks of total carbon and chlorophyll a, and a shift towards smaller phytoplankton and carnivorous copepods, associated with the seasonal impact of the EAC microbial assemblage. In light of the substantial shifts in microbial assemblage structure and function associated with the EAC, we conclude that climate-driven expansions of WBCs will expand the range of tropical oligotrophic microbes, and potentially profoundly impact the trophic status of temperate waters.

Exploring the role of spatial and stoichiometric heterogeneity in the top-down control in eutrophic planktonic ecosystems

Understanding the impact of eutrophication on the dynamics of aquatic food webs, remains a long-term challenge in ecology. Mathematical models generally predict the destabilisation of such webs, under increasing eutrophication levels, with large oscillations of species densities that can result in their extinction. This is at odds with a number of ecological observations that show stable dynamics even for high nutrient loads. The apparent discrepancy between theory and observations is known as the Rosenzweig's 'paradox of enrichment' and various solutions to the problem have been proposed over the years. In this study, we explore the stabilisation of dynamics of a tri-trophic plankton model in a eutrophic environment which occurs as a result of interplay of space heterogeneity, ecological stoichiometry, and food taxis of predators. We build a variety of models of increasing complexity, to explore various scenarios of phytoplankton growth, zooplankton food-dependent vertical movement, and different stoichiometric limitations of zooplankton. We show that the synergy among the vertical gradient in phytoplankton growth, phytoplankton structuring in terms of their stoichiometric ratio, and food-dependent vertical movement of zooplankton, would result in a postponing of destabilisation of eutrophic systems as compared to a well-mixed system. Our approach reveals a high complexity of the bifurcation structure of the system when key model parameters, such as the degree of eutrophication and light shading, are varied. We find coexistence of limit cycles and stable equilibria and that the possibility of multiple attractors in the system can result in hysteresis phenomena when the nutrient load is manipulated. These results are relevant and should be taken into account in lake restoration programs.

Global Trends in Marine Plankton Diversity across Kingdoms of Life

The ocean is home to myriad small planktonic organisms that underpin the functioning of marine ecosystems. However, their spatial patterns of diversity and the underlying drivers remain poorly known, precluding projections of their responses to global changes. Here we investigate the latitudinal gradients and global predictors of plankton diversity across archaea, bacteria, eukaryotes, and major virus clades using both molecular and imaging data from Tara Oceans. We show a decline of diversity for most planktonic groups toward the poles, mainly driven by decreasing ocean temperatures. Projections into the future suggest that severe warming of the surface ocean by the end of the 21^st century could lead to tropicalization of the diversity of most planktonic groups in temperate and polar regions. These changes may have multiple consequences for marine ecosystem functioning and services and are expected to be particularly significant in key areas for carbon sequestration, fisheries, and marine conservation.

Video Abstract

•

Most epipelagic planktonic groups exhibit a poleward decline of diversity

•

No latitudinal diversity gradient was observed below the photic zone

•

Temperature emerges as the best predictor of epipelagic plankton diversity

•

Global warming may increase plankton diversity, particularly at high latitudes

Biomagnification of Methylmercury in a Marine Plankton Ecosystem

Methylmercury is greatly bioconcentrated and biomagnified in marine plankton ecosystems, and these communities form the basis of marine food webs. Therefore, the evaluation of the potential exposure of methylmercury to higher trophic levels, including humans, requires a better understanding of its distribution in the ocean and the factors that control its biomagnification. In this study, a coupled physical/ecological model is used to simulate the trophic transfer of monomethylmercury (MMHg) in a marine plankton ecosystem. The model includes phytoplankton, a microbial community, herbivorous zooplankton (HZ), and carnivorous zooplankton (CZ). The model captures both shorter food chains in oligotrophic regions, with small HZ feeding on small phytoplankton, and longer chains in higher nutrient conditions, with larger HZ feeding on larger phytoplankton and larger CZ feeding on larger HZ. In the model, trophic dilution occurs in the food webs that involve small zooplankton, as the grazing fluxes of small zooplankton are insufficient to accumulate more MMHg in themselves than in their prey. The model suggests that biomagnification is more prominent in large zooplankton and that the microbial community plays an important role in the trophic transfer of MMHg. Sensitivity analyses show that with increasing body size, the sensitivity of the trophic magnification ratio to grazing, mortality rates, and food assimilation efficiency (AE_C) increases, while the sensitivity to excretion rates decreases. More predation or a longer zooplankton lifespan may lead to more prominent biomagnification, especially for large species. Because lower AE_C results in more predation, modeled ratios of MMHg concentrations between large plankton are doubled or even tripled when the AE_C decreases from 50% to 10%. This suggests that the biomagnification of large zooplankton is particularly sensitive to food assimilation efficiency.

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.

Small phytoplankton dominate western North Atlantic biomass

The North Atlantic phytoplankton spring bloom is the pinnacle in an annual cycle that is driven by physical, chemical, and biological seasonality. Despite its important contributions to the global carbon cycle, transitions in plankton community composition between the winter and spring have been scarcely examined in the North Atlantic. Phytoplankton composition in early winter was compared with latitudinal transects that captured the subsequent spring bloom climax. Amplicon sequence variants (ASVs), imaging flow cytometry, and flow-cytometry provided a synoptic view of phytoplankton diversity. Phytoplankton communities were not uniform across the sites studied, but rather mapped with apparent fidelity onto subpolar- and subtropical-influenced water masses of the North Atlantic. At most stations, cells < 20-µm diameter were the main contributors to phytoplankton biomass. Winter phytoplankton communities were dominated by cyanobacteria and pico-phytoeukaryotes. These transitioned to more diverse and dynamic spring communities in which pico- and nano-phytoeukaryotes, including many prasinophyte algae, dominated. Diatoms, which are often assumed to be the dominant phytoplankton in blooms, were contributors but not the major component of biomass. We show that diverse, small phytoplankton taxa are unexpectedly common in the western North Atlantic and that regional influences play a large role in modulating community transitions during the seasonal progression of blooms.

Microbial abundance, composition, and function in nectar are shaped by flower visitor identity

Microbial dispersal is essential for establishment in new habitats, but the role of vector identity is poorly understood in community assembly and function. Here, we compared microbial assembly and function in floral nectar visited by legitimate pollinators (hummingbirds) and nectar robbers (carpenter bees). We assessed effects of visitation on the abundance and composition of culturable bacteria and fungi and their taxonomy and function using shotgun metagenomics and nectar chemistry. We also compared metagenome-assembled genomes (MAGs) of Acinetobacter, a common and highly abundant nectar bacterium, among visitor treatments. Visitation increased microbial abundance, but robbing resulted in 10× higher microbial abundance than pollination. Microbial communities differed among visitor treatments: robbed flowers were characterized by predominant nectar specialists within Acetobacteraceae and Metschnikowiaceae, with a concurrent loss of rare taxa, and these resulting communities harbored genes relating to osmotic stress, saccharide metabolism and specialized transporters. Gene differences were mirrored in function: robbed nectar contained a higher percentage of monosaccharides. Draft genomes of Acinetobacter revealed distinct amino acid and saccharide utilization pathways in strains isolated from robbed versus pollinated flowers. Our results suggest an unrecognized cost of nectar robbing for pollination and distinct effects of visitor type on interactions between plants and pollinators. Overall, these results suggest vector identity is an underappreciated factor structuring microbial community assembly and function.

Combining Molecular Observations and Microbial Ecosystem Modeling: A Practical Guide

Advances in technologies for molecular observation are leading to novel types of data, including gene, transcript, protein, and metabolite levels, which are fundamentally different from the types traditionally compared with microbial ecosystem models, such as biomass (e.g., chlorophyll a) and nutrient concentrations. A grand challenge is to use these data to improve predictive models and use models to explain observed patterns. This article presents a framework that aligns observations and models along the dimension of abstraction or biological organization-from raw sequences to ecosystem patterns for observations, and from sequence simulators to ecological theory for models. It then reviews 16 studies that compared model results with molecular observations. Molecular data can and are being combined with microbial ecosystem models, but to keep up with and take advantage of the full scope of observations, models need to become more mechanistically detailed and complex, which is a technical and cultural challenge for the ecological modeling community.

Stable coexistence of equivalent nutrient competitors through niche differentiation in the light spectrum

Niche-based theories and the neutral theory of biodiversity differ in their predictions of how the species composition of natural communities will respond to changes in nutrient availability. This is an issue of major environmental relevance, as many ecosystems have experienced changes in nitrogen (N) and phosphorus (P) due to anthropogenic manipulation of nutrient loading. To understand how changes in N and P limitation may impact community structure, we conducted laboratory competition experiments using a multispecies phytoplankton community sampled from the North Sea. Results showed that picocyanobacteria (Cyanobium sp.) won the competition under N limitation, while picocyanobacteria and nonmotile nanophytoplankton (Nannochloropsis sp.) coexisted at equal abundances under P limitation. Additional experiments using isolated monocultures confirmed that Cyanobium sp. depleted N to lower levels than Nannochloropsis sp., but that both species had nearly identical P requirements, suggesting a potential for neutral coexistence under P-limited conditions. Pairwise competition experiments with the two isolates seemed to support the consistency of these results, but P limitation resulted in stable species coexistence irrespective of the initial conditions rather than the random drift of species abundances predicted by neutral theory. Comparison of the light absorption spectra indicates that coexistence of the two species was stabilized through differential use of the underwater light spectrum. Our results provide an interesting experimental example of modern coexistence theory, where species were equal competitors in one niche dimension but their competitive traits differed in other niche dimensions, thus enabling stable species coexistence on a single limiting nutrient through niche differentiation in the light spectrum.

Responses of marine phytoplankton communities to environmental changes: New insights from a niche classification scheme

Predicting changes of phytoplankton communities in response to global warming is one of the challenges of ecological forecasting. One of the constraints is the paucity of general principles applicable to community ecology. Based on a synecological analysis of a decadal-scale database, we created a niche habitat classification scheme relating nine phytoplankton groups to fifteen statistically refined realized niches comprised of three niche dimensions: temperature, irradiance, and nitrate concentrations. The niche scheme assigned the nine phytoplankton groups to three types of niches: a cold type, a warm type, and a type associated with high irradiance and high nitrate concentrations. The fact that phytoplankton groups in cold niches were governed by irradiance and those in warm niches by nitrate is consistent with general ecological theories, but the fact that diatoms were the only dominant group in high-irradiance, high-nitrate niches challenges the idea based on autecological studies that diatoms are generally better adapted to low-irradiance, high-nutrient conditions. When combined with an irradiance model, the niche scheme revealed that photoinhibition of Prochlorococcus, which is predicted from autecological studies, is a function of temperature. We used the niche scheme to predict the responses of phytoplankton communities to environmental changes due to seawater warming and eutrophication. The results of the study suggest that a synecological analysis of large databases from field studies facilitates identification of general principles of community ecology that can be used to forecast responses of biological communities to environmental changes.

Considering the Role of Adaptive Evolution in Models of the Ocean and Climate System

Numerical models have been highly successful in simulating global carbon and nutrient cycles in today's ocean, together with observed spatial and temporal patterns of chlorophyll and plankton biomass at the surface. With this success has come some confidence in projecting the century-scale response to continuing anthropogenic warming. There is also increasing interest in using such models to understand the role of plankton ecosystems in past oceans. However, today's marine environment is the product of billions of years of continual evolution-a process that continues today. In this paper, we address the questions of whether an assumption of species invariance is sufficient, and if not, under what circumstances current model projections might break down. To do this, we first identify the key timescales and questions asked of models. We then review how current marine ecosystem models work and what alternative approaches are available to account for evolution. We argue that for timescales of climate change overlapping with evolutionary timescales, accounting for evolution may to lead to very different projected outcomes regarding the timescales of ecosystem response and associated global biogeochemical cycling. This is particularly the case for past extinction events but may also be true in the future, depending on the eventual degree of anthropogenic disruption. The discipline of building new numerical models that incorporate evolution is also hugely beneficial in itself, as it forces us to question what we know about adaptive evolution, irrespective of its quantitative role in any specific event or environmental changes.

Current State of and Future Opportunities for Prediction in Microbiome Research: Report from the Mid-Atlantic Microbiome Meet-up in Baltimore on 9 January 2019

Accurate predictions across multiple fields of microbiome research have far-reaching benefits to society, but there are few widely accepted quantitative tools to make accurate predictions about microbial communities and their functions. More discussion is needed about the current state of microbiome analysis and the tools required to overcome the hurdles preventing development and implementation of predictive analyses. We summarize the ideas generated by participants of the Mid-Atlantic Microbiome Meet-up in January 2019. While it was clear from the presentations that most fields have advanced beyond simple associative and descriptive analyses, most fields lack essential elements needed for the development and application of accurate microbiome predictions. Participants stressed the need for standardization, reproducibility, and accessibility of quantitative tools as key to advancing predictions in microbiome analysis. We highlight hurdles that participants identified and propose directions for future efforts that will advance the use of prediction in microbiome research.

Diversity and biogeography of SAR11 bacteria from the Arctic Ocean

The Arctic Ocean is relatively isolated from other oceans and consists of strongly stratified water masses with distinct histories, nutrient, temperature, and salinity characteristics, therefore providing an optimal environment to investigate local adaptation. The globally distributed SAR11 bacterial group consists of multiple ecotypes that are associated with particular marine environments, yet relatively little is known about Arctic SAR11 diversity. Here, we examined SAR11 diversity using ITS analysis and metagenome-assembled genomes (MAGs). Arctic SAR11 assemblages were comprised of the S1a, S1b, S2, and S3 clades, and structured by water mass and depth. The fresher surface layer was dominated by an ecotype (S3-derived P3.2) previously associated with Arctic and brackish water. In contrast, deeper waters of Pacific origin were dominated by the P2.3 ecotype of the S2 clade, within which we identified a novel subdivision (P2.3s1) that was rare outside the Arctic Ocean. Arctic S2-derived SAR11 MAGs were restricted to high latitudes and included MAGs related to the recently defined S2b subclade, a finding consistent with bi-polar ecotypes and Arctic endemism. These results place the stratified Arctic Ocean into the SAR11 global biogeography and have identified SAR11 lineages for future investigation of adaptive evolution in the Arctic Ocean.

Scientists' warning to humanity: microorganisms and climate change

In the Anthropocene, in which we now live, climate change is impacting most life on Earth. Microorganisms support the existence of all higher trophic life forms. To understand how humans and other life forms on Earth (including those we are yet to discover) can withstand anthropogenic climate change, it is vital to incorporate knowledge of the microbial 'unseen majority'. We must learn not just how microorganisms affect climate change (including production and consumption of greenhouse gases) but also how they will be affected by climate change and other human activities. This Consensus Statement documents the central role and global importance of microorganisms in climate change biology. It also puts humanity on notice that the impact of climate change will depend heavily on responses of microorganisms, which are essential for achieving an environmentally sustainable future.

Biogeography and Diversity of Freshwater Bacteria on a River Catchment Scale

To illustrate how freshwater bacterial community changes with geographic gradient, we investigated the spatial changes of bacterial abundance and community structures from over 200 samples on a catchment scale in the Songhua River using heterotrophic plate counts, flow cytometry, denaturing gradient gel electrophoresis, and pyrosequencing analysis. The results showed that the mainstream had higher cultivable bacteria and total bacterial concentration than tributaries in the Songhua River catchment. Response model analysis demonstrated that the bacterial community exhibits a biogeographical signature even in an interconnected river network system, and the total bacterial concentration and biodiversity were significantly correlated to latitude (p < 0.001) and longitude (p < 0.001). Multivariate redundancy analysis indicated that temperature was the most important factor driving bacterial community structure in the Songhua River, which accounts for 35.30% variance of communities, then dissolved oxygen (17.60%), latitude (17.60%), longitude (11.80%), and pH (5.88%). High-throughput pyrosequencing revealed that at the phylum level, Proteobacteria was numerically dominant (89.6%) in river catchment, followed by Bacteroidetes (8.1%) and Cyanobacteria (1.2%). The overall results revealed that the bacterial community was driven by geographical distance regardless of the continuum of the river on a catchment scale.

Global pattern of phytoplankton diversity driven by temperature and environmental variability

Despite their importance to ocean productivity, global patterns of marine phytoplankton diversity remain poorly characterized. Although temperature is considered a key driver of general marine biodiversity, its specific role in phytoplankton diversity has remained unclear. We determined monthly phytoplankton species richness by using niche modeling and >540,000 global phytoplankton observations to predict biogeographic patterns of 536 phytoplankton species. Consistent with metabolic theory, phytoplankton richness in the tropics is about three times that in higher latitudes, with temperature being the most important driver. However, below 19°C, richness is lower than expected, with ~8°- 14°C waters (~35° to 60° latitude) showing the greatest divergence from theoretical predictions. Regions of reduced richness are characterized by maximal species turnover and environmental variability, suggesting that the latter reduces species richness directly, or through enhancing competitive exclusion. The nonmonotonic relationship between phytoplankton richness and temperature suggests unanticipated complexity in responses of marine biodiversity to ocean warming.

Coupling between taxonomic and functional diversity in protistan coastal communities

The study of protistan functional diversity is crucial to understand the dynamics of oceanic ecological processes. We combined the metabarcoding data of various coastal ecosystems and a newly developed trait-based approach to study the link between taxonomic and functional diversity across marine protistan communities of different size-classes. Environmental DNA was extracted and the V4 18S rDNA genomic region was amplified and sequenced. In parallel, we tried to annotate the operational taxonomic units (OTUs) from our metabarcoding dataset to 30 biological traits using published and accessible information on protists. We then developed a method to study trait correlations across protists (i.e. trade-offs) in order to build the best functional groups. Based on the annotated OTUs and our functional groups, we demonstrated that the functional diversity of marine protist communities varied in parallel with their taxonomic diversity. The coupling between functional and taxonomic diversity was conserved across different protist size classes. However, the smallest size-fraction was characterized by wider taxonomic and functional groups diversity, corroborating the idea that nanoplankton and picoplankton are part of a more stable ecological background on which larger protists and metazoans might develop.