2019 Publications

  • Generalized size scaling of metabolic rates based on single-cell measurements with freshwater phytoplankton

    Zaoli S., Giometto A., Maranon E., Escrig S., Meibom A., Ahluwalia A., Stocker R., Maritan A. and Rinaldo A.

    , 2019, PNAS, 116, 17323-17329

    Kleiber’s law describes the scaling of metabolic rate with body size across several orders of magnitude in size and across taxa and is widely regarded as a fundamental law in biology. The physiological origins of Kleiber’s law are still debated and generalizations of the law accounting for deviations from the scaling behavior have been proposed. Most theoretical and experimental studies of Kleiber’s law, however, have focused on the relationship between the average body size of a species and its mean metabolic rate, neglecting intraspecific variation of these 2 traits. Here, we propose a theoretical characterization of such variation and report on proof-of-concept experiments with freshwater phytoplankton supporting such framework. We performed joint measurements at the single-cell level of cell volume and nitrogen/carbon uptake rates, as proxies of metabolic rates, of 3 phytoplankton species using nanoscale secondary ion mass spectrometry (NanoSIMS) and stable isotope labeling. Common scaling features of the distribution of nutrient uptake rates and cell volume are found to hold across 3 orders of magnitude in cell size. Once individual measurements of cell volume and nutrient uptake rate within a species are appropriately rescaled by a function of the average cell volume within each species, we find that intraspecific distributions of cell volume and metabolic rates collapse onto a universal curve. Based on the experimental results, this work provides the building blocks for a generalized form of Kleiber’s law incorporating intraspecific, correlated variations of nutrient-uptake rates and body sizes.

  • Systematic variation in marine dissolved organic matter stoichiometry and remineralization ratios as a function of lability

    Zakem, E. J. and N. M. Levine

    , 2019, Global Biogeochemical Cycles 33 (11): 1389-1407

    Remineralization of organic matter by heterotrophic organisms regulates the biological sequestration of carbon, thereby mediating atmospheric CO2. While surface nutrient supply impacts the elemental ratios of primary production, stoichiometric control by remineralization remains unclear. Here we develop a mechanistic description of remineralization and its stoichiometry in a marine microbial ecosystem model. The model simulates the observed elemental plasticity of phytoplankton and the relatively constant, lower C:N of heterotrophic biomass. In addition, the model captures the observed decreases in the C:N of more labile dissolved organic matter (DOM) and the C:N of its remineralization with depth, which are driven by a switch in the dominant source of DOM from phytoplankton to heterotrophic biomass. Only a model version with targeted remineralization of N‐rich components is able to simulate the observed profiles of preferential remineralization of N relative to C and the elevated C:N of bulk DOM. The model suggests that more labile substrates are associated with C‐limited heterotrophic growth and not with preferential remineralization, while less labile substrates are associated with growth limited by processing rates and with preferential remineralization. The resulting patterns of variable remineralization stoichiometry mediate the extent to which a proportional increase in carbon production resulting from changes in phytoplankton stoichiometry can increase the efficiency of the biological pump. Results emphasize the importance of understanding the physiology of both phytoplankton and heterotrophs for anticipating changes in biologically driven ocean carbon storage.

  • Not just going with the flow: the effects of fluid flow on bacteria and plankton

    Wheeler J.D., Secchi E., Rusconi R. and Stocker R.

    , 2019, Annual Review of Cell and Developmental Biology, 35, 213-237

    Microorganisms often live in habitats characterized by fluid flow, from lakes and oceans to soil and the human body. Bacteria and plankton experience a broad range of flows, from the chaotic motion characteristic of turbulence to smooth flows at boundaries and in confined environments. Flow creates forces and torques that affect the movement, behavior, and spatial distribution of microorganisms and shapes the chemical landscape on which they rely for nutrient acquisition and communication. Methodological advances and closer interactions between physicists and biologists have begun to reveal the importance of flow–microorganism interactions and the adaptations of microorganisms to flow. Here we review selected examples of such interactions from bacteria, phytoplankton, larvae, and zooplankton. We hope that this article will serve as a blueprint for a more in-depth consideration of the effects of flow in the biology of microorganisms and that this discussion will stimulate further multidisciplinary effort in understanding this important component of microorganism habitats.

  • Bacterial chemotaxis in a microfluidic T-maze reveals strong phenotypic heterogeneity in chemotactic sensitivity

    Salek M.M., Carrara F., Fernandez V.I., Guasto J.S. and Stocker R.

    , 2019, Nature Communications, 10, 1877

    Many microorganisms have evolved chemotactic strategies to exploit the microscale heterogeneity that frequently characterizes microbial habitats. Chemotaxis has been primarily studied as an average characteristic of a population, with little regard for variability among individuals. Here, we adopt a classic tool from animal ecology – the T-maze – and implement it at the microscale by using microfluidics to expose bacteria to a sequence of decisions, each consisting of migration up or down a chemical gradient. Single-cell observations of clonal Escherichia coli in the maze, coupled with a mathematical model, reveal that strong heterogeneity in the chemotactic sensitivity coefficient exists even within clonal populations of bacteria. A comparison of different potential sources of heterogeneity reveals that heterogeneity in the T-maze originates primarily from the chemotactic sensitivity coefficient, arising from a distribution of pathway gains. This heterogeneity may have a functional role, for example in the context of migratory bet-hedging strategies.

  • Bacterial maze runners reveal hidden diversity in chemotactic performance

    Salek M.M., Carrara F., Fernandez V.I. and Stocker R.

    , 2019, Microbial Cell, 6, 370-372

    Chemotaxis allows microorganisms to exploit gradients in chemical stimuli to find nutrient resources and hosts or escape noxious substances. Thus, the life of individual microbes in their natural environments is a continual sequence of decisions based on the perceived chemical gradients. However, it has remained unclear to what extent the chemotaxis properties vary among cells of one species, and whether there is a spectrum of different ‘decision makers’ within populations of bacteria. In our recent study (Salek, Carrara et al., Nature Communications 10 (1), 1877), we combine microfluidic experiments with mathematical modeling to demonstrate that even in clonal populations, bacteria are individuals with different abilities to climb chemical gradients.

  • The role of microbial motility and chemotaxis in symbiosis

    Raina J.B., Fernandez V.I., Lambert B., Stocker R. and Seymour J.R.

    , 2019, Nature Reviews Microbiology, 17, 284-294

    Many symbiotic relationships rely on the acquisition of microbial partners from the environment. However, the mechanisms by which microbial symbionts find and colonize their hosts are often unknown. We propose that the acquisition of environmental symbionts often necessitates active migration and colonization by the symbionts through motility and chemotaxis. The pivotal role of these behaviours in the onset and maintenance of symbiotic interactions is well established in a small number of model systems but remains largely overlooked for the many symbioses that involve the recruitment of microbial partners from the environment. In this Review, we highlight when, where and how chemotaxis and motility can enable symbiont recruitment and propose that these symbiont behaviours are important across a wide range of hosts and environments.

  • The molecular clock in the evolution of protein structures

    Pascual-García, A. and Arenas, M. and Bastolla, U.

    , 2019, Systematic Biology, 68(6),987-1002

    The molecular clock hypothesis, which states that substitutions accumulate in protein sequences at a constant rate, plays a fundamental role in molecular evolution but it is violated when selective or mutational processes vary with time. Such violations of the molecular clock have been widely investigated for protein sequences, but not yet for protein structures. Here, we introduce a novel statistical test (Significant Clock Violations) and perform a large scale assessment of the molecular clock in the evolution of both protein sequences and structures in three large superfamilies. After validating our method with computer simulations, we find that clock violations are generally consistent in sequence and structure evolution, but they tend to be larger and more significant in structure evolution. Moreover, changes of function assessed through Gene Ontology and InterPro terms are associated with large and significant clock violations in structure evolution. We found that almost one third of significant clock violations are significant in structure evolution but not in sequence evolution, highlighting the advantage to use structure information for assessing accelerated evolution and gathering hints of positive selection. Clock violations between closely related pairs are frequently significant in sequence evolution, consistent with the observed time dependence of the substitution rate attributed to segregation of neutral and slightly deleterious polymorphisms, but not in structure evolution, suggesting that these substitutions do not affect protein structure although they may affect stability. These results are consistent with the view that natural selection, both negative and positive, constrains more strongly protein structures than protein sequences. Our code for computing clock violations is freely available at https://github.com/ugobas/Molecular_clock.

  • Chain formation can enhance the vertical migration of phytoplankton through turbulence

    Lovecchio S., Climent E., Stocker R. and Durham W.M.

    , 2019, Science Advances, 5, eaaw7879.

    Many species of motile phytoplankton can actively form long multicellular chains by remaining attached to one another after cell division. While chains swim more rapidly than single cells of the same species, chain formation also markedly reduces phytoplankton’s ability to maintain their bearing. This suggests that turbulence, which acts to randomize swimming direction, could sharply attenuate a chain’s ability to migrate between well-lit surface waters during the day and deeper nutrient-rich waters at night. Here, we use numerical models to investigate how chain formation affects the migration of phytoplankton through a turbulent water column. Unexpectedly, we find that the elongated shape of chains helps them travel through weak to moderate turbulence much more effectively than single cells, and isolate the physical processes that confer chains this ability. Our findings provide a new mechanistic understanding of how turbulence can select for phytoplankton with elongated morphologies and may help explain why turbulence triggers chain formation.

  • A foraging mandala for aquatic microorganisms

    Fernandez V.I., Yawata Y. and Stocker R.

    , 2019, ISME Journal, 13, 563-575.

    Aquatic environments harbor a great diversity of microorganisms, which interact with the same patchy, particulate, or diffuse resources by means of a broad array of physiological and behavioral adaptations, resulting in substantially different life histories and ecological success. To date, efforts to uncover and understand this diversity have not been matched by equivalent efforts to identify unifying frameworks that can provide a degree of generality and thus serve as a stepping stone to scale up microscale dynamics to predict their ecosystem-level consequences. In particular, evaluating the ecological consequences of different resource landscapes and of different microbial adaptations has remained a major challenge in aquatic microbial ecology. Here, inspired by Ramon Margalef’s mandala for phytoplankton, we propose a foraging mandala for microorganisms in aquatic environments, which accounts for both the local environment and individual adaptations. This biophysical framework distills resource acquisition into two fundamental parameters: the search time for a new resource and the growth return obtained from encounter with a resource. We illustrate the foraging mandala by considering a broad range of microbial adaptations and environmental characteristics. The broad applicability of the foraging mandala suggests that it could be a useful framework to compare disparate microbial strategies in aquatic environments and to reduce the vast complexity of microbe-environment interactions into a minimal number of fundamental parameters.

  • Modular Assembly of Polysaccharide-Degrading Marine Microbial Communities

    Enke, T. N., M.S. Datta, J. Schwartzman, N. Cermak, D. Schmitz, J. Barrere, A.O. Pascual Garcia, O.X. Cordero.

    , 2019, Curr Biol 29.9 (2019): 1528-1535.

    Understanding the principles that govern the assembly of microbial communities across earth’s biomes is a major challenge in modern microbial ecology. This pursuit is complicated by the difficulties of mapping functional roles and interactions onto communities with immense taxonomic diversity and of identifying the scale at which microbes interact. To address this challenge, here, we focused on the bacterial communities that colonize and degrade particulate organic matter in the ocean. We show that the assembly of these communities can be simplified as a linear combination of functional modules. Using synthetic polysaccharide particles immersed in natural bacterioplankton assemblages, we showed that successional particle colonization dynamics are driven by the interaction of two types of modules: a first type made of narrowly specialized primary degraders, whose dynamics are controlled by particle polysaccharide composition, and a second type containing substrate-independent taxa whose dynamics are controlled by interspecific interactions—in particular, cross-feeding via organic acids, amino acids, and other metabolic byproducts. We show that, as a consequence of this trophic structure, communities can assemble modularly—i.e., by a simple sum of substrate-specific primary degrader modules, one for each complex polysaccharide in the particle, connected to a single broad-niche range consumer module. Consistent with this model, a linear combination of the communities on single-polysaccharide particles accurately predicts community composition on mixed-polysaccharide particles. Our results suggest that the assembly of heterotrophic communities that degrade complex organic materials follows simple design principles that could be exploited to engineer heterotrophic microbiomes.

  • Cooperation and spatial self-organization determine rate and efficiency of particulate organic matter degradation in marine bacteria

    Ebrahimi, A., Schwartzman, J. and Cordero, O. X.

    , 2019, Proceedings of the National Academy of Science, November 12, 2019 116 (46) 23309-23316

    The recycling of particulate organic matter (POM) by microbes is a key part of the global carbon cycle. This process is mediated by the extracellular hydrolysis of polysaccharides, which can trigger social behaviors in bacteria resulting from the production of public goods. Despite the potential importance of public good-mediated interactions, their relevance in the environment remains unclear. In this study, we developed a computational and experimental model system to address this challenge and studied how the POM depolymerization rate and its uptake efficiency (2 main ecosystem function parameters) depended on social interactions and spatial self-organization on particle surfaces. We found an emergent trade-off between rate and efficiency resulting from the competition between oligosaccharide diffusion and cellular uptake, with low rate and high efficiency being achieved through cell-to-cell cooperation between degraders. Bacteria cooperated by aggregating in cell clusters of ∼10 to 20 µm, in which cells were able to share public goods. This phenomenon, which was independent of any explicit group-level regulation, led to the emergence of critical cell concentrations below which degradation did not occur, despite all resources being available in excess. In contrast, when particles were labile and turnover rates were high, aggregation promoted competition and decreased the efficiency of carbon use. Our study shows how social interactions and cell aggregation determine the rate and efficiency of particulate carbon turnover in environmentally relevant scenarios.

  • Multicellular behaviour enables cooperation in microbial cell aggregates

    Ebrahimi, A., Schwartzman, J. and Cordero, O. X.

    , 2019, Phil. Trans. R. Soc.

    Oligosaccharides produced from the extracellular hydrolysis of biological materials can act as common goods that promote cooperative growth in microbial populations, whereby cell–cell aggregation increases both the per capita availability of resources and the per-cell growth rate. However, aggregation can also have detrimental consequences for growth, as gradients form within aggregates limiting the resource accessibility. We built a computational model, which predicts cooperation is restricted in dense cell aggregates larger than 10 µm because of the emergence of polymer and oligomer counter gradients. We compared these predictions to experiments performed with two well-studied alginate-degrading strains of Vibrio splendidus, which varied in their ability to secrete alginate lyase. We observed that both strains can form large aggregates (less than 50 µm), overcoming diffusion limitation by rearranging their internal structure. The stronger enzyme producer grew non-cooperatively and formed aggregates with internal channels that allowed exchange between the bulk environment and the aggregate, whereas the weak enzyme producer showed strongly cooperative growth and formed dense aggregates in which cells near the core mixed by active swimming. Our simulations suggest that the mixing and channelling reduce diffusion limitation and allow cells to uniformly grow in aggregates. Together, these data demonstrate that bacterial behaviour can help overcome competition imposed by resource gradients within cell aggregates.

  • NanoSIMS imaging reveals metabolic stratification within current-producing biofilms

    Chadwick, G.L., Otero, F.J., Gralnick, J.A., Bond, D.R. and Orphan, V.J.

    , 2019, Proceedings of the National Academy of Sciences, 116(41), pp.20716-20724.

    Metal-reducing bacteria direct electrons to their outer surfaces, where insoluble metal oxides or electrodes act as terminal electron acceptors, generating electrical current from anaerobic respiration. Geobacter sulfurreducens is a commonly enriched electricity-producing organism, forming thick conductive biofilms that magnify total activity by supporting respiration of cells not in direct contact with electrodes. Hypotheses explaining why these biofilms fail to produce higher current densities suggest inhibition by formation of pH, nutrient, or redox potential gradients; but these explanations are often contradictory, and a lack of direct measurements of cellular growth within biofilms prevents discrimination between these models. To address this fundamental question, we measured the anabolic activity of G. sulfurreducens biofilms using stable isotope probing coupled to nanoscale secondary ion mass spectrometry (nanoSIMS). Our results demonstrate that the most active cells are at the anode surface, and that this activity decreases with distance, reaching a minimum 10 µm from the electrode. Cells nearest the electrode continue to grow at their maximum rate in weeks-old biofilms 80-µm-thick, indicating nutrient or buffer diffusion into the biofilm is not rate-limiting. This pattern, where highest activity occurs at the electrode and declines with each cell layer, is present in thin biofilms (<5 µm) and fully grown biofilms (>20 µm), and at different anode redox potentials. These results suggest a growth penalty is associated with respiring insoluble electron acceptors at micron distances, which has important implications for improving microbial electrochemical devices as well as our understanding of syntrophic associations harnessing the phenomenon of microbial conductivity.

  • Bacteria push the limits of chemotactic precision to navigate dynamic chemical gradients

    Brumley D.R., Carrara F., Hein A.M., Yawata Y., Levin S.A. and Stocker R.

    , 2019, PNAS, 116, 10792-10797.

    Ephemeral aggregations of bacteria are ubiquitous in the environment, where they serve as hotbeds of metabolic activity, nutrient cycling, and horizontal gene transfer. In many cases, these regions of high bacterial concentration are thought to form when motile cells use chemotaxis to navigate to chemical hotspots. However, what governs the dynamics of bacterial aggregations is unclear. Here, we use an experimental platform to create realistic submillimeter-scale nutrient pulses with controlled nutrient concentrations. By combining experiments, mathematical theory, and agent-based simulations, we show that individual Vibrio ordalii bacteria begin chemotaxis toward hotspots of dissolved organic matter (DOM) when the magnitude of the chemical gradient rises sufficiently far above the sensory noise that is generated by stochastic encounters with chemoattractant molecules. Each DOM hotspot is surrounded by a dynamic ring of chemotaxing cells, which congregate in regions of high DOM concentration before dispersing as DOM diffuses and gradients become too noisy for cells to respond to. We demonstrate that V. ordaliioperates close to the theoretical limits on chemotactic precision. Numerical simulations of chemotactic bacteria, in which molecule counting noise is explicitly taken into account, point at a tradeoff between nutrient acquisition and the cost of chemotactic precision. More generally, our results illustrate how limits on sensory precision can be used to understand the location, spatial extent, and lifespan of bacterial behavioral responses in ecologically relevant environments.