Search

Microbial communities are ubiquitous and play crucial roles in many natural processes. Despite their importance for the environment, industry and human health, there are still many aspects of microbial community dynamics that we do not understand quantitatively. Recent experiments have shown that the structure and composition of microbial communities are intertwined with the metabolism of the species that inhabit them, suggesting that properties at the intracellular level such as the allocation of cellular proteomic resources must be taken into account when describing microbial communities with a population dynamics approach. In this work, we reconsider one of the theoretical frameworks most commonly used to model population dynamics in competitive ecosystems, MacArthur’s consumer-resource model, in light of experimental evidence showing how proteome allocation affects microbial growth. This new framework allows us to describe community dynamics at an intermediate level of complexity between classical consumer-resource models and biochemical models of microbial metabolism, accounting for temporally-varying proteome allocation subject to constraints on growth and protein synthesis in the presence of multiple resources, while preserving analytical insight into the dynamics of the system. We first show with a simple experiment that proteome allocation needs to be accounted for to properly understand the dynamics of even the simplest microbial community, i.e. two bacterial strains competing for one common resource. Then, we study our consumer-proteome-resource model analytically and numerically to determine the conditions that allow multiple species to coexist in systems with arbitrary numbers of species and resources.

Antagonistic interactions are widespread in the microbial world and affect microbial evolutionary dynamics. Natural microbial communities often display spatial structure, which affects biological interactions, but much of what we know about microbial antagonism comes from laboratory studies of well-mixed communities. To overcome this limitation, we manipulated two killer strains of the budding yeast Saccharomyces cerevisiae, expressing different toxins, to independently control the rate at which they released their toxins. We developed mathematical models that predict the experimental dynamics of competition between toxin-producing strains in both well-mixed and spatially structured populations. In both situations, we experimentally verified theory’s prediction that a stronger antagonist can invade a weaker one only if the initial invading population exceeds a critical frequency or size. Finally, we found that toxin-resistant cells and weaker killers arose in spatially structured competitions between toxin-producing strains, suggesting that adaptive evolution can affect the outcome of microbial antagonism in spatial settings.

Microbial communities are ubiquitous and play crucial roles in many natural processes. Despite their importance for the environment, industry and human health, there are still many aspects of microbial community dynamics that we do not understand quantitatively. Recent experiments have shown that the metabolism of species in a community is intertwined with its composition, suggesting that properties at the intracellular level such as the allocation of cellular proteomic resources must be taken into account when describing microbial communities with a population dynamics approach. In this work we reconsider one of the theoretical frameworks most commonly used to model population dynamics in competitive ecosystems, MacArthur’s consumer-resource model, in light of experimental evidence showing how pro-teome allocation affects microbial growth. This new framework allows us to describe community dynamics at an intermediate level of complexity between classical consumer-resource models and biochemical models of microbial metabolism, accounting for temporally-varying proteome allocation subject to constraints on growth and protein synthesis in the presence of multiple resources, while preserving analytical insight into the dynamics of the system. We first show experimentally that proteome allocation needs to be accounted for to properly understand the dynamics of even the simplest microbial community, i.e. two bacterial strains competing for one common resource. We then study the model analytically and numerically to determine the conditions that allow multiple species to coexist in systems with arbitrary numbers of species and resources.

Significance Microbes often live in dense populations such as colonies and biofilms. We show that the success and extinction of yeast strains within a growing colony are determined by a combination of their relative fitness and the forces exerted by proliferating cells on their neighbors. These physical interactions prolong the survival of less-fit strains at the growing frontier of the colony and slow down the colony’s takeover by fitter strains. This reduction in the power of natural selection favors the maintenance of genetic diversity in environments in which the strains’ relative growth rates vary with time. Growth-induced physical interactions may thus favor the maintenance of drug-resistant cells, which are typically less fit than nonresistant cells, within dense microbial populations., Microbial populations often assemble in dense populations in which proliferating individuals exert mechanical forces on the nearby cells. Here, we use yeast strains whose doubling times depend differently on temperature to show that physical interactions among cells affect the competition between different genotypes in growing yeast colonies. Our experiments demonstrate that these physical interactions have two related effects: they cause the prolonged survival of slower-growing strains at the actively-growing frontier of the colony and cause faster-growing strains to increase their frequency more slowly than expected in the absence of physical interactions. These effects also promote the survival of slower-growing strains and the maintenance of genetic diversity in colonies grown in time-varying environments. A continuum model inspired by overdamped hydrodynamics reproduces the experiments and predicts that the strength of natural selection depends on the width of the actively growing layer at the colony frontier. We verify these predictions experimentally. The reduced power of natural selection observed here may favor the maintenance of drug-resistant cells in microbial populations and could explain the apparent neutrality of interclone competition within tumors.

Phototaxis, the process through which motile organisms direct their swimming toward or away from light, is implicated in key ecological phenomena (including algal blooms and diel vertical migration) that shape the distribution, diversity, and productivity of phytoplankton and thus energy transfer to higher trophic levels in aquatic ecosystems. Phototaxis also finds important applications in biofuel reactors and microbiopropellers and is argued to serve as a benchmark for the study of biological invasions in heterogeneous environments owing to the ease of generating stochastic light fields. Despite its ecological and technological relevance, an experimentally tested, general theoretical model of phototaxis seems unavailable to date. Here, we present accurate measurements of the behavior of the alga Euglena gracilis when exposed to controlled light fields. Analysis of E. gracilis' phototactic accumulation dynamics over a broad range of light intensities proves that the classic Keller-Segel mathematical framework for taxis provides an accurate description of both positive and negative phototaxis only when phototactic sensitivity is modeled by a generalized "receptor law," a specific nonlinear response function to light intensity that drives algae toward beneficial light conditions and away from harmful ones. The proposed phototactic model captures the temporal dynamics of both cells' accumulation toward light sources and their dispersion upon light cessation. The model could thus be of use in integrating models of vertical phytoplankton migrations in marine and freshwater ecosystems, and in the design of bioreactors.

Classical models of biological invasions assess species spread in homogeneous landscapes by assuming constant growth rates and random local movement. Mounting evidence suggests, however, that demographic stochasticity, environmental heterogeneity and non-random movement of individuals affect considerably the spread dynamics. Here, we show that the dynamics of biological invasions are controlled by the spatial heterogeneity of the resource distribution. We show theoretically that increasing the landscape resource autocorrelation length causes a reduction in the average speed of species spread. Demographic stochasticity plays a key role in the slowdown, which is streghtened when individuals can actively move towards resources. The reduction in the front propagation speed is verified in laboratory microcosm experiments with the flagellated protistEuglena gracilisby comparing spread in habitats characterized by different resource heterogeneity. Our theoretical and experimental findings highlight the need to account for the intrinsic stochasticity of population dynamics to describe spread in spatially extended landscapes, which are inevitably characterized by heterogeneous spatial distributions of resources controlling vital rates. Our work identifies the resource autocorrelation length as a key modulator and a simple measure of landscape susceptibility to biological invasions, with implications for predicting the characters of biological invasions within naturally heterogeneous environmental corridors.

1. Natural communities commonly contain many different species and functional groups, and multiple types of species interactions act simultaneously, such as competition, predation, commensalism or mutualism. However, experimental and theoretical investigations have generally been limited by focusing on one type of interaction at a time or by a lack of a common methodological and conceptual approach to measure species interactions. 2. We compared four methods to measure and express species interactions. These approaches are, with increasing degree of model complexity, an extinction-based model, a relative yield model and two generalized Lotka-Volterra (LV) models. All four approaches have been individually applied in different fields of community ecology, but rarely integrated. We provide an overview of the definitions, assumptions and data needed for the specific methods and apply them to empirical data by experimentally deriving the interaction matrices among 11 protist and rotifer species, belonging to three functional groups. Furthermore, we compare their advantages and limitations to predict multispecies community dynamics and ecosystem functioning. 3. The relative yield method is, in terms of final biomass production, the best method in predicting the 11-species community dynamics from the pairwise competition experiments. The LV model, which is considering equilibrium among the species, suffers from experimental constraints given the strict equilibrium assumption, and this may be rarely satisfied in ecological communities. 4. We show how simulations of a LV stochastic community model, derived from an empirical interaction matrix, can be used to predict multispecies community dynamics across multiple functional groups. 5.Our work unites available tools to measure species interactions under one framework. This improves our ability to make management-oriented predictions of species coexistence/extinction and to compare ecosystem processes across study systems.

The size of an organism matters for its metabolic, growth, mortality, and other vital rates. Scale-free community size spectra (i.e., size distributions regardless of species) are routinely observed in natural ecosystems and are the product of intra- and interspecies regulation of the relative abundance of organisms of different sizes. Intra- and interspecies distributions of body sizes are thus major determinants of ecosystems’ structure and function. We show experimentally that single-species mass distributions of unicellular eukaryotes covering different phyla exhibit both characteristic sizes and universal features over more than four orders of magnitude in mass. Remarkably, we find that the mean size of a species is sufficient to characterize its size distribution fully and that the latter has a universal form across all species. We show that an analytical physiological model accounts for the observed universality, which can be synthesized in a log-normal form for the intraspecies size distributions. We also propose how ecological and physiological processes should interact to produce scale-invariant community size spectra and discuss the implications of our results on allometric scaling laws involving body mass.

Unveiling the mechanisms that promote coexistence in biological communities is a fundamental problem in ecology. Stable coexistence of many species is commonly observed in natural communities. Most of these natural communities, however, are composed of species from multiple trophic and functional groups, while theory and experiments on coexistence have been focusing on functionally similar species. Here, we investigated how functional diversity affects the stability of species coexistence and productivity in multispecies communities by characterizing experimentally all pairwise species interactions in a pool of 11 species of eukaryotes (10 protists and one rotifer) belonging to three different functional groups. Species within the same functional group showed stronger competitive interactions compared to among-functional group interactions. This often led to competitive exclusion between species that had higher functional relatedness, but only at low levels of species richness. Communities with higher functional diversity resulted in increased species coexistence and community biomass production. Our experimental findings and the results of a stochastic model tailored to the experimental interaction matrix suggest the emergence of strong stabilizing forces when species from different functional groups interact in a homogeneous environment. By combining theoretical analysis with experiments we could also disentangle the relationship between species richness and functional diversity, showing that functional diversity per se is a crucial driver of productivity and stability in multispecies community.