Metabolic trade-offs constrain the cell size ratio in a nitrogen-fixing symbiosis.

Biological dinitrogen (N2) fixation is a key metabolic process exclusively performed by prokaryotes, some of which are symbiotic with eukaryotes. Species of the marine haptophyte algae Braarudosphaera bigelowii harbor the N2-fixing endosymbiotic cyanobacteria UCYN-A, which might be evolving organelle-like characteristics. We found that the size ratio between UCYN-A and their hosts is strikingly conserved across sublineages/species, which is consistent with the size relationships of organelles in this symbiosis and other species. Metabolic modeling showed that this size relationship maximizes the coordinated growth rate based on trade-offs between resource acquisition and exchange. Our findings show that the size relationships of N2-fixing endosymbionts and organelles in unicellular eukaryotes are constrained by predictable metabolic underpinnings and that UCYN-A is, in many regards, functioning like a hypothetical N2-fixing organelle (or nitroplast).

Flexible B12 ecophysiology of Phaeocystis antarctica due to a fusion B12-independent methionine synthase with widespread homologues.

Coastal Antarctic marine ecosystems are significant in carbon cycling because of their intense seasonal phytoplankton blooms. Southern Ocean algae are primarily limited by light and iron (Fe) and can be co-limited by cobalamin (vitamin B12). Micronutrient limitation controls productivity and shapes the composition of blooms which are typically dominated by either diatoms or the haptophyte Phaeocystis antarctica. However, the vitamin requirements and ecophysiology of the keystone species P. antarctica remain poorly characterized. Using cultures, physiological analysis, and comparative omics, we examined the response of P. antarctica to a matrix of Fe-B12 conditions. We show that P. antarctica is not auxotrophic for B12, as previously suggested, and identify mechanisms underlying its B12 response in cultures of predominantly solitary and colonial cells. A combination of proteomics and proteogenomics reveals a B12-independent methionine synthase fusion protein (MetE-fusion) that is expressed under vitamin limitation and interreplaced with the B12-dependent isoform under replete conditions. Database searches return homologues of the MetE-fusion protein in multiple Phaeocystis species and in a wide range of marine microbes, including other photosynthetic eukaryotes with polymorphic life cycles as well as bacterioplankton. Furthermore, we find MetE-fusion homologues expressed in metaproteomic and metatranscriptomic field samples in polar and more geographically widespread regions. As climate change impacts micronutrient availability in the coastal Southern Ocean, our finding that P. antarctica has a flexible B12 metabolism has implications for its relative fitness compared to B12-auxotrophic diatoms and for the detection of B12-stress in a more diverse set of marine microbes.

Trophic interactions with heterotrophic bacteria limit the range of Prochlorococcus.

Prochlorococcus is both the smallest and numerically most abundant photosynthesizing organism on the planet. While thriving in the warm oligotrophic gyres, Prochlorococcus concentrations drop rapidly in higher-latitude regions. Transect data from the North Pacific show the collapse occurring at a wide range of temperatures and latitudes (temperature is often hypothesized to cause this shift), suggesting an ecological mechanism may be at play. An often used size-based theory of phytoplankton community structure that has been incorporated into computational models correctly predicts the dominance of Prochlorococcus in the gyres, and the relative dominance of larger cells at high latitudes. However, both theory and computational models fail to explain the poleward collapse. When heterotrophic bacteria and predators that prey nonspecifically on both Prochlorococcus and bacteria are included in the theoretical framework, the collapse of Prochlorococcus occurs with increasing nutrient supplies. The poleward collapse of Prochlorococcus populations then naturally emerges when this mechanism of "shared predation" is implemented in a complex global ecosystem model. Additionally, the theory correctly predicts trends in both the abundance and mean size of the heterotrophic bacteria. These results suggest that ecological controls need to be considered to understand the biogeography of Prochlorococcus and predict its changes under future ocean conditions. Indirect interactions within a microbial network can be essential in setting community structure.

A nutrient relay sustains subtropical ocean productivity.

The expansive gyres of the subtropical ocean account for a significant fraction of global organic carbon export from the upper ocean. In the gyre interior, vertical mixing and the heaving of nutrient-rich waters into the euphotic layer sustain local productivity, in turn depleting the layers below. However, the nutrient pathways by which these subeuphotic layers are themselves replenished remain unclear. Using a global, eddy-permitting simulation of ocean physics and biogeochemistry, we quantify nutrient resupply mechanisms along and across density surfaces, including the contribution of eddy-scale motions that are challenging to observe. We find that mesoscale eddies (10 to 100 km) flux nutrients from the shallow flanks of the gyre into the recirculating interior, through time-varying motions along density surfaces. The subeuphotic layers are ultimately replenished in approximately equal contributions by this mesoscale eddy transport and the remineralization of sinking particles. The mesoscale eddy resupply is most important in the lower thermocline for the whole subtropical region but is dominant at all depths within the gyre interior. Subtropical gyre productivity may therefore be sustained by a nutrient relay, where the lateral transport resupplies nutrients to the thermocline and allows vertical exchanges to maintain surface biological production and carbon export.

A model of time-dependent macromolecular and elemental composition of phytoplankton.

Phytoplankton Chl:C:N:P ratios are important from both an ecological and a biogeochemical perspective. We show that these elemental ratios can be represented by a phytoplankton physiological model of low complexity that includes major cellular macromolecular pools. In particular, our model resolves time-dependent intracellular pools of chlorophyll, proteins, nucleic acids, carbohydrates/lipids, and N and P storage. Batch culture data for two diatom and two prasinophyte species are used to constrain parameters that represent specific allocation traits and strategies. A key novelty is the simultaneous estimation of physiological parameters for two phytoplankton groups of such different sizes. The number of free parameters is reduced by assuming (i) allometric scaling for maximum uptake rates, (ii) shared half-saturation constants for synthesis of functional macromolecules, (iii) shared exudation rates of functional macromolecules across the species. The rationale behind this assumption is that across the different species, the same or similar processes, enzymes, and metabolites play a role in key physiological processes. For the turnover numbers of macromolecular synthesis and storage exudation rates, differences between diatoms and prasinophytes need to be taken into account to obtain a good fit. Our model fits suggest that the parameters related to storage dynamics dominate the differences in the C:N:P ratios between the different phytoplankton groups. Since descriptions of storage dynamics are still incomplete and imprecise, predictions of C:N:P ratios by phytoplankton models likely have a large uncertainty.

Reverse engineering environmental metatranscriptomes clarifies best practices for eukaryotic assembly.

BACKGROUND: Diverse communities of microbial eukaryotes in the global ocean provide a variety of essential ecosystem services, from primary production and carbon flow through trophic transfer to cooperation via symbioses. Increasingly, these communities are being understood through the lens of omics tools, which enable high-throughput processing of diverse communities. Metatranscriptomics offers an understanding of near real-time gene expression in microbial eukaryotic communities, providing a window into community metabolic activity. RESULTS: Here we present a workflow for eukaryotic metatranscriptome assembly, and validate the ability of the pipeline to recapitulate real and manufactured eukaryotic community-level expression data. We also include an open-source tool for simulating environmental metatranscriptomes for testing and validation purposes. We reanalyze previously published metatranscriptomic datasets using our metatranscriptome analysis approach. CONCLUSION: We determined that a multi-assembler approach improves eukaryotic metatranscriptome assembly based on recapitulated taxonomic and functional annotations from an in-silico mock community. The systematic validation of metatranscriptome assembly and annotation methods provided here is a necessary step to assess the fidelity of our community composition measurements and functional content assignments from eukaryotic metatranscriptomes.

Microbial feedbacks optimize ocean iron availability.

Iron is the limiting factor for biological production over a large fraction of the surface ocean because free iron is rapidly scavenged or precipitated under aerobic conditions. Standing stocks of dissolved iron are maintained by association with organic molecules (ligands) produced by biological processes. We hypothesize a positive feedback between iron cycling, microbial activity, and ligand abundance: External iron input fuels microbial production, creating organic ligands that support more iron in seawater, leading to further macronutrient consumption until other microbial requirements such as macronutrients or light become limiting, and additional iron no longer increases productivity. This feedback emerges in numerical simulations of the coupled marine cycles of macronutrients and iron that resolve the dynamic microbial production and loss of iron-chelating ligands. The model solutions resemble modern nutrient distributions only over a finite range of prescribed ligand source/sink ratios where the model ocean is driven to global-scale colimitation by micronutrients and macronutrients and global production is maximized. We hypothesize that a global-scale selection for microbial ligand cycling may have occurred to maintain "just enough" iron in the ocean.

Anthropogenic Asian aerosols provide Fe to the North Pacific Ocean.

Fossil-fuel emissions may impact phytoplankton primary productivity and carbon cycling by supplying bioavailable Fe to remote areas of the ocean via atmospheric aerosols. However, this pathway has not been confirmed by field observations of anthropogenic Fe in seawater. Here we present high-resolution trace-metal concentrations across the North Pacific Ocean (158°W from 25°to 42°N). A dissolved Fe maximum was observed around 35°N, coincident with high dissolved Pb and Pb isotope ratios matching Asian industrial sources and confirming recent aerosol deposition. Iron-stable isotopes reveal in situ evidence of anthropogenic Fe in seawater, with low δ56Fe (-0.23‰ > δ56Fe > -0.65‰) observed in the region that is most influenced by aerosol deposition. An isotope mass balance suggests that anthropogenic Fe contributes 21-59% of dissolved Fe measured between 35° and 40°N. Thus, anthropogenic aerosol Fe is likely to be an important Fe source to the North Pacific Ocean.

A steady-state model of microbial acclimation to substrate limitation.

Microbes acclimate to changes in substrate availability by altering the number of transporters on the cell surface, however there is some disagreement on just how. We revisit the physics of substrate uptake and consider the steady-state scenario whereby cells have acclimated to maximize fitness. Flux balance analysis of a stoichiometric model of Escherichia coli was used in conjunction with quantitative proteomics data and molecular modeling of membrane transporters to reconcile these opposing views. An emergent feature of the proposed model is a critical substrate concentration S*, which delineates two rate limits. At concentrations above S*, transporter abundance can be regulated to maintain uptake rates as demanded by maximal growth rates, whereas below S*, uptake rates are strictly diffusion limited. In certain scenarios, the proposed model can take on a qualitatively different shape from the familiar hyperbolic kinetics curves, instead resembling the long-forgotten Blackman kinetics.

Moving ecological and biogeochemical transitions across the North Pacific.

In the North Pacific Ocean, nutrient rich surface waters flow south from the subpolar gyre through a transitional region and into the subtropics. Along the way, nutrients are used, recycled, and exported, leading to lower biomass and a commensurate change in ecosystem structure moving southward. We focus on the region between the two gyres (the Transition Zone) using a coupled biophysical ocean model, remote sensing, floats, and cruise data to explore the nature of the physical, biogeochemical, and ecological fields in this region. Nonlinear interactions between biological processes and the meridional gradient in nutrient supply lead to sharp shifts across this zone. These transitions between a southern region with more uniform biological and biogeochemical properties and steep meridional gradients to the north are diagnosed from extrema in the first derivative of the properties with latitude. Some transitions like that for chlorophyll a (the transition zone chlorophyll front [TZCF]) experience large seasonal excursions while the location of the transitions in other properties moves very little. The seasonal shifts are not caused by changes in the horizontal flow field, but rather by the interaction of seasonal, depth related, forcing with the mean latitudinal gradients. Focusing on the TZCF as a case study, we express its phase velocity in terms of vertical nutrient flux and internal ecosystem processes, demonstrating their nearly equal influence on its motion. This framework of propagating biogeochemical transitions can be systematically expanded to better understand the processes that structure ecosystems and biogeochemistry in the North Pacific and beyond.

Generalized Stoichiometry and Biogeochemistry for Astrobiological Applications.

A central need in the field of astrobiology is generalized perspectives on life that make it possible to differentiate abiotic and biotic chemical systems McKay (2008). A key component of many past and future astrobiological measurements is the elemental ratio of various samples. Classic work on Earth's oceans has shown that life displays a striking regularity in the ratio of elements as originally characterized by Redfield (Redfield 1958; Geider and La Roche 2002; Eighty years of Redfield 2014). The body of work since the original observations has connected this ratio with basic ecological dynamics and cell physiology, while also documenting the range of elemental ratios found in a variety of environments. Several key questions remain in considering how to best apply this knowledge to astrobiological contexts: How can the observed variation of the elemental ratios be more formally systematized using basic biological physiology and ecological or environmental dynamics? How can these elemental ratios be generalized beyond the life that we have observed on our own planet? Here, we expand recently developed generalized physiological models (Kempes et al. 2012, 2016, 2017, 2019) to create a simple framework for predicting the variation of elemental ratios found in various environments. We then discuss further generalizing the physiology for astrobiological applications. Much of our theoretical treatment is designed for in situ measurements applicable to future planetary missions. We imagine scenarios where three measurements can be made-particle/cell sizes, particle/cell stoichiometry, and fluid or environmental stoichiometry-and develop our theory in connection with these often deployed measurements.

Metabolic evolution and the self-organization of ecosystems.

Metabolism mediates the flow of matter and energy through the biosphere. We examined how metabolic evolution shapes ecosystems by reconstructing it in the globally abundant oceanic phytoplankter Prochlorococcus To understand what drove observed evolutionary patterns, we interpreted them in the context of its population dynamics, growth rate, and light adaptation, and the size and macromolecular and elemental composition of cells. This multilevel view suggests that, over the course of evolution, there was a steady increase in Prochlorococcus' metabolic rate and excretion of organic carbon. We derived a mathematical framework that suggests these adaptations lower the minimal subsistence nutrient concentration of cells, which results in a drawdown of nutrients in oceanic surface waters. This, in turn, increases total ecosystem biomass and promotes the coevolution of all cells in the ecosystem. Additional reconstructions suggest that Prochlorococcus and the dominant cooccurring heterotrophic bacterium SAR11 form a coevolved mutualism that maximizes their collective metabolic rate by recycling organic carbon through complementary excretion and uptake pathways. Moreover, the metabolic codependencies of Prochlorococcus and SAR11 are highly similar to those of chloroplasts and mitochondria within plant cells. These observations lead us to propose a general theory relating metabolic evolution to the self-amplification and self-organization of the biosphere. We discuss the implications of this framework for the evolution of Earth's biogeochemical cycles and the rise of atmospheric oxygen.

An empirical model of carbon flow through marine viruses and microzooplankton grazers.

Viruses and microzooplankton grazers represent major sources of mortality for marine phytoplankton and bacteria, redirecting the flow of organic material throughout the world's oceans. Here, we investigate the use of nonlinear population models of interactions between phytoplankton, viruses and grazers as a means to quantitatively constrain the flow of carbon through marine microbial ecosystems. We augment population models with a synthesis of laboratory-based estimates of prey, predator and viral life history traits that constrain transfer efficiencies. We then apply the model framework to estimate loss rates in the California Current Ecosystem (CCE). With our empirically parameterized model, we estimate that, of the total losses mediated by viruses and microzooplankton grazing at the focal CCE site, 22 ± 3%, 46 ± 27%, 3 ± 2% and 29 ± 20% were directed to grazers, sloppy feeding (as well as excretion and respiration), viruses and viral lysate respectively. We identify opportunities to leverage ecosystem models and conventional mortality assays to further constrain the quantitative rates of critical ecosystem processes.

Biochemical diversity of glycosphingolipid biosynthesis as a driver of Coccolithovirus competitive ecology.

Coccolithoviruses (EhVs) are large, double-stranded DNA-containing viruses that infect the single-celled, marine coccolithophore Emiliania huxleyi. Given the cosmopolitan nature and global importance of E. huxleyi as a bloom-forming, calcifying, photoautotroph, E. huxleyi-EhV interactions play a key role in oceanic carbon biogeochemistry. Virally-encoded glycosphingolipids (vGSLs) are virulence factors that are produced by the activity of virus-encoded serine palmitoyltransferase (SPT). Here, we characterize the dynamics, diversity and catalytic production of vGSLs in an array of EhV strains in relation to their SPT sequence composition and explore the hypothesis that they are a determinant of infectivity and host demise. vGSL production and diversity was positively correlated with increased virulence, virus replication rate and lytic infection dynamics in laboratory experiments, but they do not explain the success of less-virulent EhVs in natural EhV communities. The majority of EhV-derived SPT amplicon sequences associated with infected cells in the North Atlantic derived from slower infecting, less virulent EhVs. Our lab-, field- and mathematical model-based data and simulations support ecological scenarios whereby slow-infecting, less-virulent EhVs successfully compete in North Atlantic populations of E. huxleyi, through either the preferential removal of fast-infecting, virulent EhVs during active infection or by having access to a broader host range.

Light regulation of coccolithophore host-virus interactions.

Viruses that infect photoautotrophs have a fundamental relationship with light, given the need for host resources. We investigated the role of light on Coccolithovirus (EhV) infection of the globally distributed coccolithophore, Emiliania huxleyi. Light was required for EhV adsorption, and viral production was highest when host cultures were maintained in continuous light or at irradiance levels of 150-300 µmol m-2  s-1 . During the early stages of infection, photosynthetic electron transport remained high, while RuBisCO expression decreased concomitant with an induction of the pentose phosphate pathway, the primary source of de novo nucleotides. A mathematical model developed and fitted to the laboratory data supported the hypothesis that EhV replication was controlled by a trade-off between host nucleotide recycling and de novo synthesis, and that photoperiod and photon flux could toggle this switch. Laboratory results supported field observations that light was the most robust driver of EhV replication within E. huxleyi populations collected across a 2000 nautical mile transect in the North Atlantic. Collectively, these findings demonstrate that light can drive host-virus interactions through a mechanistic interplay between host metabolic processes, which serve to structure infection and phytoplankton mortality in the upper ocean.

Marine mixotrophy increases trophic transfer efficiency, mean organism size, and vertical carbon flux.

Mixotrophic plankton, which combine the uptake of inorganic resources and the ingestion of living prey, are ubiquitous in marine ecosystems, but their integrated biogeochemical impacts remain unclear. We address this issue by removing the strict distinction between phytoplankton and zooplankton from a global model of the marine plankton food web. This simplification allows the emergence of a realistic trophic network with increased fidelity to empirical estimates of plankton community structure and elemental stoichiometry, relative to a system in which autotrophy and heterotrophy are mutually exclusive. Mixotrophy enhances the transfer of biomass to larger sizes classes further up the food chain, leading to an approximately threefold increase in global mean organism size and an ∼35% increase in sinking carbon flux.

Morphological optimization for access to dual oxidants in biofilms.

A major theme driving research in biology is the relationship between form and function. In particular, a longstanding goal has been to understand how the evolution of multicellularity conferred fitness advantages. Here we show that biofilms of the bacterium Pseudomonas aeruginosa produce structures that maximize cellular reproduction. Specifically, we develop a mathematical model of resource availability and metabolic response within colony features. This analysis accurately predicts the measured distribution of two types of electron acceptors: oxygen, which is available from the atmosphere, and phenazines, redox-active antibiotics produced by the bacterium. Using this model, we demonstrate that the geometry of colony structures is optimal with respect to growth efficiency. Because our model is based on resource dynamics, we also can anticipate shifts in feature geometry based on changes to the availability of electron acceptors, including variations in the external availability of oxygen and genetic manipulation that renders the cells incapable of phenazine production.

A bottom-up perspective on ecosystem change in Mesozoic oceans.

Mesozoic and Early Cenozoic marine animals across multiple phyla record secular trends in morphology, environmental distribution, and inferred behaviour that are parsimoniously explained in terms of increased selection pressure from durophagous predators. Another systemic change in Mesozoic marine ecosystems, less widely appreciated than the first, may help to explain the observed animal record. Fossils, biomarker molecules, and molecular clocks indicate a major shift in phytoplankton composition, as mixotrophic dinoflagellates, coccolithophorids and, later, diatoms radiated across shelves. Models originally developed to probe the ecology and biogeography of modern phytoplankton enable us to evaluate the ecosystem consequences of these phytoplankton radiations. In particular, our models suggest that the radiation of mixotrophic dinoflagellates and the subsequent diversification of marine diatoms would have accelerated the transfer of primary production upward into larger size classes and higher trophic levels. Thus, phytoplankton evolution provides a mechanism capable of facilitating the observed evolutionary shift in Mesozoic marine animals.

Growth, metabolic partitioning, and the size of microorganisms.

Population growth rate is a fundamental ecological and evolutionary characteristic of living organisms, but individuals must balance the metabolism devoted to biosynthesis and reproduction against the maintenance of existing structure and other functionality. Here we present a mathematical model that relates metabolic partitioning to the form of growth. The model captures the observed growth trajectory of single cells and individuals for a variety of species and taxa spanning prokaryotes, unicellular eukaryotes, and small multicellular eukaryotes. Our analysis suggests that the per-unit costs of biosynthesis and maintenance are conserved across prokaryotes and eukaryotes. However, the relative metabolic expenditure on growth and maintenance of whole organisms clearly differentiates taxa: prokaryotes spend an increasing fraction of their entire metabolism on growth with increasing cell size, whereas eukaryotes devote a diminishing fraction. These differences allow us to predict the minimum and maximum size for each taxonomic group, anticipating observed evolutionary life-history transitions. The framework provides energetic insights into taxonomic tradeoffs related to growth and metabolism and constrains traits that are important for size-structured modeling of microbial communities and their ecological and biogeochemical effects.