Phytoplankton diversity explained by connectivity across a mesoscale frontal system in the open ocean

Phytoplankton community composition is important in establishing ecosystem structure and function. Intuitively, we recognize that water movements must be important for modifying spatial gradients and plankton diversity. However, identifying boundaries and exchange between habitats in the open ocean is not straightforward. Here, we use the abundance of nine phytoplankton species closely sampled in a mesoscale frontal system in the northeastern North Sea as a proxy for community composition and explore the relationship between phytoplankton biogeography and transport patterns. Subsurface community distributions could be related to modeled patterns in water movement. A methodology for analyzing pelagic diversity that includes a representation of plankton community composition and an Eulerian connectivity tracer was developed, and the relative importance of connectivity and geographical distance for phytoplankton species composition analyzed. The connectivity tracer identifies timescales and dispersal barriers in the open ocean. Connectivity was found to be superior in explaining pelagic plankton diversity and found to be a prerequisite for understanding the pelagic phytoplankton composition. This approach is a valuable tool for establishing the link between ocean transports, ecosystem structure and biodiversity and for informing the placement of marine protected areas.

Trophic strategies explain the ocean niches of small eukaryotic phytoplankton

A large fraction of marine primary production is performed by diverse small protists, and many of these phytoplankton are phagotrophic mixotrophs that vary widely in their capacity to consume bacterial prey. Prior analyses suggest that mixotrophic protists as a group vary in importance across ocean environments, but the mechanisms leading to broad functional diversity among mixotrophs, and the biogeochemical consequences of this, are less clear. Here we use isolates from seven major taxa to demonstrate a tradeoff between phototrophic performance (growth in the absence of prey) and phagotrophic performance (clearance rate when consuming Prochlorococcus). We then show that trophic strategy along the autotrophy-mixotrophy spectrum correlates strongly with global niche differences, across depths and across gradients of stratification and chlorophyll a. A model of competition shows that community shifts can be explained by greater fitness of faster-grazing mixotrophs when nutrients are scarce and light is plentiful. Our results illustrate how basic physiological constraints and principles of resource competition can organize complexity in the surface ocean ecosystem.

Genomic evidence for global ocean plankton biogeography shaped by large-scale current systems

Biogeographical studies have traditionally focused on readily visible organisms, but recent technological advances are enabling analyses of the large-scale distribution of microscopic organisms, whose biogeographical patterns have long been debated. Here we assessed the global structure of plankton geography and its relation to the biological, chemical, and physical context of the ocean (the ‘seascape’) by analyzing metagenomes of plankton communities sampled across oceans during the Tara Oceans expedition, in light of environmental data and ocean current transport. Using a consistent approach across organismal sizes that provides unprecedented resolution to measure changes in genomic composition between communities, we report a pan-ocean, size-dependent plankton biogeography overlying regional heterogeneity. We found robust evidence for a basin-scale impact of transport by ocean currents on plankton biogeography, and on a characteristic timescale of community dynamics going beyond simple seasonality or life history transitions of plankton.

Oceans are brimming with life invisible to our eyes, a myriad of species of bacteria, viruses and other microscopic organisms essential for the health of the planet. These ‘marine plankton’ are unable to swim against currents and should therefore be constantly on the move, yet previous studies have suggested that distinct species of plankton may in fact inhabit different oceanic regions. However, proving this theory has been challenging; collecting plankton is logistically difficult, and it is often impossible to distinguish between species simply by examining them under a microscope. However, within the last decade, a research schooner called Tara has travelled the globe to gather thousands of plankton samples. At the same time, advances in genomics have made it possible to identify species based only on fragments of their DNA sequence.

To understand the hidden geography of plankton communities in Earth’s oceans, Richter et al. pored over DNA from the Tara Oceans expedition. This revealed that, despite being unable to resist the flow of water, various planktonic species which live close to the surface manage to occupy distinct, stable provinces shaped by currents. Different sizes of plankton are distributed in different sized provinces, with the smallest organisms tending to inhabit the smallest areas. Comparing DNA similarities and speeds of currents at the ocean surface revealed how these might stretch and mix plankton communities.

Plankton play a critical role in the health of the ocean and the chemical cycles of planet Earth. These results could allow deeper investigation by marine modellers, ecologists, and evolutionary biologists. Meanwhile, work is already underway to investigate how climate change might impact this hidden geography.

Global drivers of eukaryotic plankton biogeography in the sunlit ocean

[Figure: see text].

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.

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.

Stable aerobic and anaerobic coexistence in anoxic marine zones

Mechanistic description of the transition from aerobic to anaerobic metabolism is necessary for diagnostic and predictive modeling of fixed nitrogen loss in anoxic marine zones (AMZs). In a metabolic model where diverse oxygen- and nitrogen-cycling microbial metabolisms are described by underlying redox chemical reactions, we predict a transition from strictly aerobic to predominantly anaerobic regimes as the outcome of ecological interactions along an oxygen gradient, obviating the need for prescribed critical oxygen concentrations. Competing aerobic and anaerobic metabolisms can coexist in anoxic conditions whether these metabolisms represent obligate or facultative populations. In the coexistence regime, relative rates of aerobic and anaerobic activity are determined by the ratio of oxygen to electron donor supply. The model simulates key characteristics of AMZs, such as the accumulation of nitrite and the sustainability of anammox at higher oxygen concentrations than denitrification, and articulates how microbial biomass concentrations relate to associated water column transformation rates as a function of redox stoichiometry and energetics. Incorporating the metabolic model into an idealized two-dimensional ocean circulation results in a simulated AMZ, in which a secondary chlorophyll maximum emerges from oxygen-limited grazing, and where vertical mixing and dispersal in the oxycline also contribute to metabolic co-occurrence. The modeling approach is mechanistic yet computationally economical and suitable for global change applications.

The role of submesoscale currents in structuring marine ecosystems

From microbes to large predators, there is increasing evidence that marine life is shaped by short-lived submesoscales currents that are difficult to observe, model, and explain theoretically. Whether and how these intense three-dimensional currents structure the productivity and diversity of marine ecosystems is a subject of active debate. Our synthesis of observations and models suggests that the shallow penetration of submesoscale vertical currents might limit their impact on productivity, though ecological interactions at the submesoscale may be important in structuring oceanic biodiversity.

Short-lived three-dimensional submesoscale currents, responsible for swirling ocean color chlorophyll filaments, have long been thought to affect productivity. Current research suggests they may not be effective in enhancing phytoplankton growth, but may have important contributions to biodiversity.

Ecological control of nitrite in the upper ocean

Microorganisms oxidize organic nitrogen to nitrate in a series of steps. Nitrite, an intermediate product, accumulates at the base of the sunlit layer in the subtropical ocean, forming a primary nitrite maximum, but can accumulate throughout the sunlit layer at higher latitudes. We model nitrifying chemoautotrophs in a marine ecosystem and demonstrate that microbial community interactions can explain the nitrite distributions. Our theoretical framework proposes that nitrite can accumulate to a higher concentration than ammonium because of differences in underlying redox chemistry and cell size between ammonia- and nitrite-oxidizing chemoautotrophs. Using ocean circulation models, we demonstrate that nitrifying microorganisms are excluded in the sunlit layer when phytoplankton are nitrogen-limited, but thrive at depth when phytoplankton become light-limited, resulting in nitrite accumulation there. However, nitrifying microorganisms may coexist in the sunlit layer when phytoplankton are iron- or light-limited (often in higher latitudes). These results improve understanding of the controls on nitrification, and provide a framework for representing chemoautotrophs and their biogeochemical effects in ocean models.

Nitrite tends to peak at the base of the sunlit zone in the ocean, but the ecological drivers of the local and global distributions of nitrite are not known. Here, Zakem et al. use a marine ecosystem model to show how the interactions of nitrifying microbes mediate nitrite accumulation.

Large-scale ocean connectivity and planktonic body size

Global patterns of planktonic diversity are mainly determined by the dispersal of propagules with ocean currents. However, the role that abundance and body size play in determining spatial patterns of diversity remains unclear. Here we analyse spatial community structure - β-diversity - for several planktonic and nektonic organisms from prokaryotes to small mesopelagic fishes collected during the Malaspina 2010 Expedition. β-diversity was compared to surface ocean transit times derived from a global circulation model, revealing a significant negative relationship that is stronger than environmental differences. Estimated dispersal scales for different groups show a negative correlation with body size, where less abundant large-bodied communities have significantly shorter dispersal scales and larger species spatial turnover rates than more abundant small-bodied plankton. Our results confirm that the dispersal scale of planktonic and micro-nektonic organisms is determined by local abundance, which scales with body size, ultimately setting global spatial patterns of diversity.

Body size is hypothesised to be a major determinant of β-diversity in passively-dispersing marine organisms. Here, Villarino et al. show that plankton body size determines rates of dispersal along marine currents, with shorter dispersal and higher species spatial turnover in larger organisms.

Dispersal similarly shapes both population genetics and community patterns in the marine realm

Dispersal plays a key role to connect populations and, if limited, is one of the main processes to maintain and generate regional biodiversity. According to neutral theories of molecular evolution and biodiversity, dispersal limitation of propagules and population stochasticity are integral to shaping both genetic and community structure. We conducted a parallel analysis of biological connectivity at genetic and community levels in marine groups with different dispersal traits. We compiled large data sets of population genetic structure (98 benthic macroinvertebrate and 35 planktonic species) and biogeographic data (2193 benthic macroinvertebrate and 734 planktonic species). We estimated dispersal distances from population genetic data (i.e., FST vs. geographic distance) and from β-diversity at the community level. Dispersal distances ranked the biological groups in the same order at both genetic and community levels, as predicted by organism dispersal ability and seascape connectivity: macrozoobenthic species without dispersing larvae, followed by macrozoobenthic species with dispersing larvae and plankton (phyto- and zooplankton). This ranking order is associated with constraints to the movement of macrozoobenthos within the seabed compared with the pelagic habitat. We showed that dispersal limitation similarly determines the connectivity degree of communities and populations, supporting the predictions of neutral theories in marine biodiversity patterns.

Estimating planktonic diversity through spatial dominance patterns in a model ocean

In the open ocean, the observation and quantification of biodiversity patterns is challenging. Marine ecosystems are indeed largely composed by microbial planktonic communities whose niches are affected by highly dynamical physico-chemical conditions, and whose observation requires advanced methods for morphological and molecular classification. Optical remote sensing offers an appealing complement to these in-situ techniques. Global-scale coverage at high spatiotemporal resolution is however achieved at the cost of restrained information on the local assemblage. Here, we use a coupled physical and ecological model ocean simulation to explore one possible metrics for comparing measures performed on such different scales. We show that a large part of the local diversity of the virtual plankton ecosystem - corresponding to what accessible by genomic methods - can be inferred from crude, but spatially extended, information - as conveyed by remote sensing. Shannon diversity of the local community is indeed highly correlated to a 'seascape' index, which quantifies the surrounding spatial heterogeneity of the most abundant functional group. The error implied in drastically reducing the resolution of the plankton community is shown to be smaller in frontal regions as well as in regions of intermediate turbulent energy. On the spatial scale of hundreds of kms, patterns of virtual plankton diversity are thus largely sustained by mixing communities that occupy adjacent niches. We provide a proof of principle that in the open ocean information on spatial variability of communities can compensate for limited local knowledge, suggesting the possibility of integrating in-situ and satellite observations to monitor biodiversity distribution at the global scale.

Drift in ocean currents impacts intergenerational microbial exposure to temperature

Microbes are the foundation of marine ecosystems [Falkowski PG, Fenchel T, Delong EF (2008) Science 320(5879):1034-1039]. Until now, the analytical framework for understanding the implications of ocean warming on microbes has not considered thermal exposure during transport in dynamic seascapes, implying that our current view of change for these critical organisms may be inaccurate. Here we show that upper-ocean microbes experience along-trajectory temperature variability up to 10 °C greater than seasonal fluctuations estimated in a static frame, and that this variability depends strongly on location. These findings demonstrate that drift in ocean currents can increase the thermal exposure of microbes and suggests that microbial populations with broad thermal tolerance will survive transport to distant regions of the ocean and invade new habitats. Our findings also suggest that advection has the capacity to influence microbial community assemblies, such that regions with strong currents and large thermal fluctuations select for communities with greatest plasticity and evolvability, and communities with narrow thermal performance are found where ocean currents are weak or along-trajectory temperature variation is low. Given that fluctuating environments select for individual plasticity in microbial lineages, and that physiological plasticity of ancestors can predict the magnitude of evolutionary responses of subsequent generations to environmental change [Schaum CE, Collins S (2014) Proc Biol Soc 281(1793):20141486], our findings suggest that microbial populations in the sub-Antarctic (∼40°S), North Pacific, and North Atlantic will have the most capacity to adapt to contemporary ocean warming.

Anthropogenic climate change drives shift and shuffle in North Atlantic phytoplankton communities

Anthropogenic climate change has shifted the biogeography and phenology of many terrestrial and marine species. Marine phytoplankton communities appear sensitive to climate change, yet understanding of how individual species may respond to anthropogenic climate change remains limited. Here, using historical environmental and phytoplankton observations, we characterize the realized ecological niches for 87 North Atlantic diatom and dinoflagellate taxa and project changes in species biogeography between mean historical (1951-2000) and future (2051-2100) ocean conditions. We find that the central positions of the core range of 74% of taxa shift poleward at a median rate of 12.9 km per decade (km⋅dec(-1)), and 90% of taxa shift eastward at a median rate of 42.7 km⋅dec(-1) The poleward shift is faster than previously reported for marine taxa, and the predominance of longitudinal shifts is driven by dynamic changes in multiple environmental drivers, rather than a strictly poleward, temperature-driven redistribution of ocean habitats. A century of climate change significantly shuffles community composition by a basin-wide median value of 16%, compared with seasonal variations of 46%. The North Atlantic phytoplankton community appears poised for marked shift and shuffle, which may have broad effects on food webs and biogeochemical cycles.

The Impact of Submesoscale Physics on Primary Productivity of Plankton

Life in the ocean relies on the photosynthetic production of phytoplankton, which is influenced by the availability of light and nutrients that are modulated by a host of physical processes. Submesoscale processes are particularly relevant to phytoplankton productivity because the timescales on which they act are similar to those of phytoplankton growth. Their dynamics are associated with strong vorticity and strain rates that occur on lateral scales of 0.1-10 km. They can support vertical velocities as large as 100 m d(-1) and play a crucial role in transporting nutrients into the sunlit ocean for phytoplankton production. In regimes with deep surface mixed layers, submesoscale instabilities can cause stratification within days, thereby increasing light exposure for phytoplankton trapped close to the surface. These instabilities help to create and maintain localized environments that favor the growth of phytoplankton, contribute to productivity, and cause enormous heterogeneity in the abundance of phytoplankton, which has implications for interactions within the ecosystem.

The dynamical landscape of marine phytoplankton diversity

Observations suggest that the landscape of marine phytoplankton assemblage might be strongly heterogeneous at the dynamical mesoscale and submesoscale (10-100 km, days to months), with potential consequences in terms of global diversity and carbon export. But these variations are not well documented as synoptic taxonomic data are difficult to acquire. Here, we examine how phytoplankton assemblage and diversity vary between mesoscale eddies and submesoscale fronts. We use a multi-phytoplankton numerical model embedded in a mesoscale flow representative of the North Atlantic. Our model results suggest that the mesoscale flow dynamically distorts the niches predefined by environmental contrasts at the basin scale and that the phytoplankton diversity landscape varies over temporal and spatial scales that are one order of magnitude smaller than those of the basin-scale environmental conditions. We find that any assemblage and any level of diversity can occur in eddies and fronts. However, on a statistical level, the results suggest a tendency for larger diversity and more fast-growing types at fronts, where nutrient supplies are larger and where populations of adjacent water masses are constantly brought into contact; and lower diversity in the core of eddies, where water masses are kept isolated long enough to enable competitive exclusion.

Mechanisms of Physical-Biological-Biogeochemical Interaction at the Oceanic Mesoscale

Mesoscale phenomena are ubiquitous and highly energetic features of ocean circulation. Their influence on biological and biogeochemical processes varies widely, stemming not only from advective transport but also from the generation of variations in the environment that affect biological and chemical rates. The ephemeral nature of mesoscale features in the ocean makes it difficult to elucidate the attendant mechanisms of physical-biological-biogeochemical interaction, necessitating the use of multidisciplinary approaches involving in situ observations, remote sensing, and modeling. All three aspects are woven through this review in an attempt to synthesize current understanding of the topic, with particular emphasis on novel developments in recent years.

How turbulence regulates biodiversity in systems with cyclic competition

Cyclic, nonhierarchical interactions among biological species represent a general mechanism by which ecosystems are able to maintain high levels of biodiversity. However, species coexistence is often possible only in spatially extended systems with a limited range of dispersal, whereas in well-mixed environments models for cyclic competition often lead to a loss of biodiversity. Here we consider the dispersal of biological species in a fluid environment, where mixing is achieved by a combination of advection and diffusion. In particular, we perform a detailed numerical analysis of a model composed of turbulent advection, diffusive transport, and cyclic interactions among biological species in two spatial dimensions and discuss the circumstances under which biodiversity is maintained when external environmental conditions, such as resource supply, are uniform in space. Cyclic interactions are represented by a model with three competitors, resembling the children's game of rock-paper-scissors, whereas the flow field is obtained from a direct numerical simulation of two-dimensional turbulence with hyperviscosity. It is shown that the space-averaged dynamics undergoes bifurcations as the relative strengths of advection and diffusion compared to biological interactions are varied.