Endosymbiotic relationships are pervasive across diverse taxa of life, offering key avenues for eco-evolutionary dynamics. Although a variety of experimental and empirical frameworks have shed light on critical aspects of endosymbiosis, theoretical frameworks (mathematical models) are especially well-suited for certain tasks. Mathematical models can integrate multiple factors to determine the net outcome of endosymbiotic relationships, identify broad patterns that connect endosymbioses with other systems, simplify biological complexity, generate hypotheses for underlying mechanisms, evaluate different hypotheses, identify constraints that limit certain biological interactions, and open new lines of inquiry. This Essay highlights the utility of mathematical models in endosymbiosis research, particularly in generating relevant hypotheses. Despite their limitations, mathematical models can be used to address known unknowns and discover unknown unknowns.

Multicellular life forms have evolved many times on our planet, suggesting that this is a common evolutionary innovation. Multiple advantages have been proposed for the emergence of multicellularity (MC). In this paper, we address the problem of how the first precondition for MC, namely 'stay together', might have occurred under spatially limited resources exploited by a population of unicellular agents. Using a minimal model of evolved cell-cell adhesion among growing and dividing cells that exploit a localized resource with a given size, we show that a transition occurs at a critical resource size separating a phase of evolved multicellular aggregates from a phase where unicellularity (UC) is favoured. The two phases are separated by an intermediate domain where both UC and MC can be selected by evolution. This model provides a minimal approach to the early stages that were required to transition from individuality to cohesive groups of cells associated with a physical cooperative effect: when resources are present only in a localized portion of the habitat, MC is a desirable property as it helps cells to keep close to the available local nutrients.

Changing environments pose a challenge to living organisms. Cells need to gather and process incoming information, adapting to changes in predictable ways. This requires in particular the presence of memory, which allows different internal states to be stored. Biological memory can be stored by switches that retain information on past and present events. Synthetic biologists have implemented a number of memory devices for biological applications, mostly in single cells. It has been shown that the use of multicellular consortia provides interesting advantages to implement biological circuits. Here we show how to build a synthetic biological memory switch using an eukaryotic consortium. We engineered yeast cells that can communicate and retain memory of changes in the extracellular environment. These cells were able to produce and secrete a pheromone and sense a different pheromone following NOT logic. When the two strains were cocultured, they behaved as a double-negative-feedback motif with memory. In addition, we showed that memory can be effectively changed by the use of external inputs. Further optimization of these modules and addition of other cells could lead to new multicellular circuits that exhibit memory over a broad range of biological inputs.

The rise of multicellularity in the early evolution of life represents a major challenge for evolutionary biology. Guidance for finding answers has emerged from disparate fields, from phylogenetics to modelling and synthetic biology, but little is known about the potential origins of multicellular aggregates before genetic programmes took full control of developmental processes. Such aggregates should involve spatial organization of differentiated cells and the modification of flows and concentrations of metabolites within well-defined boundaries. Here, we show that, in an environment where limited nutrients and toxic metabolites are introduced, a population of cells capable of stochastic differentiation and differential adhesion can develop into multicellular aggregates with conflict mediation mechanisms and a complex internal structure. The morphospace of possible patterns is shown to be very rich, including proto-organisms that display a high degree of organizational complexity, far beyond simple heterogeneous populations of cells. Our findings reveal that there is a potentially enormous richness of organismal complexity between simple mixed cooperators and embodied living organisms.

Spatial self-organization emerges in distributed systems exhibiting local interactions when nonlinearities and the appropriate propagation of signals are at work. These kinds of phenomena can be modeled with different frameworks, typically cellular automata or reaction-diffusion systems. A different class of dynamical processes involves the correlated movement of agents over space, which can be mediated through chemotactic movement or minimization of cell-cell interaction energy. A classic example of the latter is given by the formation of spatially segregated assemblies when cells display differential adhesion. Here, we consider a new class of dynamical models, involving cell adhesion among two stochastically exchangeable cell states as a minimal model capable of exhibiting well-defined, ordered spatial patterns. Our results suggest that a whole space of pattern-forming rules is hosted by the combination of physical differential adhesion and the value of probabilities modulating cell phenotypic switching, showing that Turing-like patterns can be obtained without resorting to reaction-diffusion processes. If the model is expanded allowing cells to proliferate and die in an environment where diffusible nutrient and toxic waste are at play, different phases are observed, characterized by regularly spaced patterns. The analysis of the parameter space reveals that certain phases reach higher population levels than other modes of organization. A detailed exploration of the mean-field theory is also presented. Finally, we let populations of cells with different adhesion matrices compete for reproduction, showing that, in our model, structural organization can improve the fitness of a given cell population. The implications of these results for ecological and evolutionary models of pattern formation and the emergence of multicellularity are outlined.

Cells are complex machines capable of processing information by means of an entangled network of molecular interactions. A crucial component of these decision-making systems is the presence of memory and this is also a specially relevant target of engineered synthetic systems. A classic example of memory devices is a 1-bit memory element known as the flip-flop. Such system can be in principle designed using a single-cell implementation, but a direct mapping between standard circuit design and a living circuit can be cumbersome. Here we present a novel computational implementation of a 1-bit memory device using a reliable multicellular design able to behave as a set-reset flip-flop that could be implemented in yeast cells. The dynamics of the proposed synthetic circuit is investigated with a mathematical model using biologically-meaningful parameters. The circuit is shown to behave as a flip-flop in a wide range of parameter values. The repression strength for the NOT logics is shown to be crucial to obtain a good flip-flop signal. Our model also shows that the circuit can be externally tuned to achieve different memory states and dynamics, such as persistent and transient memory. We have characterized the parameter domains for robust memory storage and retrieval as well as the corresponding time response dynamics.