Modeling endosymbioses: Insights and hypotheses from theoretical approaches.

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.

Metabolic compatibility and the rarity of prokaryote endosymbioses.

The evolution of the mitochondria was a significant event that gave rise to the eukaryotic lineage and most large complex life. Central to the origins of the mitochondria was an endosymbiosis between prokaryotes. Yet, despite the potential benefits that can stem from a prokaryotic endosymbiosis, their modern occurrence is exceptionally rare. While many factors may contribute to their rarity, we lack methods for estimating the extent to which they constrain the appearance of a prokaryotic endosymbiosis. Here, we address this knowledge gap by examining the role of metabolic compatibility between a prokaryotic host and endosymbiont. We use genome-scale metabolic flux models from three different collections (AGORA, KBase, and CarveMe) to assess the viability, fitness, and evolvability of potential prokaryotic endosymbioses. We find that while more than half of host-endosymbiont pairings are metabolically viable, the resulting endosymbioses have reduced growth rates compared to their ancestral metabolisms and are unlikely to gain mutations to overcome these fitness differences. In spite of these challenges, we do find that they may be more robust in the face of environmental perturbations at least in comparison with the ancestral host metabolism lineages. Our results provide a critical set of null models and expectations for understanding the forces that shape the structure of prokaryotic life.

Minor variations in multicellular life cycles have major effects on adaptation.

Multicellularity has evolved several independent times over the past hundreds of millions of years and given rise to a wide diversity of complex life. Recent studies have found that large differences in the fundamental structure of early multicellular life cycles can affect fitness and influence multicellular adaptation. Yet, there is an underlying assumption that at some scale or categorization multicellular life cycles are similar in terms of their adaptive potential. Here, we consider this possibility by exploring adaptation in a class of simple multicellular life cycles of filamentous organisms that only differ in one respect, how many daughter filaments are produced. We use mathematical models and evolutionary simulations to show that despite the similarities, qualitatively different mutations fix. In particular, we find that mutations with a tradeoff between cell growth and group survival, i.e. "selfish" or "altruistic" traits, spread differently. Specifically, altruistic mutations more readily spread in life cycles that produce few daughters while in life cycles producing many daughters either type of mutation can spread depending on the environment. Our results show that subtle changes in multicellular life cycles can fundamentally alter adaptation.

Overexpression of the plastidial pseudo-protease AtFtsHi3 enhances drought tolerance while sustaining plant growth.

With climate change, droughts are expected to be more frequent and severe, severely impacting plant biomass and quality. Here, we show that overexpressing the Arabidopsis gene AtFtsHi3 (FtsHi3OE) enhances drought-tolerant phenotypes without compromising plant growth. AtFtsHi3 encodes a chloroplast envelope pseudo-protease; knock-down mutants (ftshi3-1) are found to be drought tolerant but exhibit stunted growth. Altered AtFtsHi3 expression therefore leads to drought tolerance, while only diminished expression of this gene leads to growth retardation. To understand the underlying mechanisms of the enhanced drought tolerance, we compared the proteomes of ftshi3-1 and pFtsHi3-FtsHi3OE (pFtsHi3-OE) to wild-type plants under well-watered and drought conditions. Drought-related processes like osmotic stress, water transport, and abscisic acid response were enriched in pFtsHi3-OE and ftshi3-1 mutants following their enhanced drought response compared to wild-type. The knock-down mutant ftshi3-1 showed an increased abundance of HSP90, HSP93, and TIC110 proteins, hinting at a potential downstream role of AtFtsHi3 in chloroplast pre-protein import. Mathematical modeling was performed to understand how variation in the transcript abundance of AtFtsHi3 can, on the one hand, lead to drought tolerance in both overexpression and knock-down lines, yet, on the other hand, affect plant growth so differently. The results led us to hypothesize that AtFtsHi3 may form complexes with at least two other protease subunits, either as homo- or heteromeric structures. Enriched amounts of AtFtsH7/9, AtFtsH11, AtFtsH12, and AtFtsHi4 in ftshi3-1 suggest a possible compensation mechanism for these proteases in the hexamer.

Syntrophy emerges spontaneously in complex metabolic systems.

Syntrophy allows a microbial community as a whole to survive in an environment, even though individual microbes cannot. The metabolic interdependence typical of syntrophy is thought to arise from the accumulation of degenerative mutations during the sustained co-evolution of initially self-sufficient organisms. An alternative and underexplored possibility is that syntrophy can emerge spontaneously in communities of organisms that did not co-evolve. Here, we study this de novo origin of syntrophy using experimentally validated computational techniques to predict an organism's viability from its metabolic reactions. We show that pairs of metabolisms that are randomly sampled from a large space of possible metabolism and viable on specific primary carbon sources often become viable on new carbon sources by exchanging metabolites. The same biochemical reactions that are required for viability on primary carbon sources also confer viability on novel carbon sources. Our observations highlight a new and important avenue for the emergence of metabolic adaptations and novel ecological interactions.

Vitamin D receptor regulates autophagy in the normal mammary gland and in luminal breast cancer cells.

Women in North America have a one in eight lifetime risk of developing breast cancer (BC), and a significant proportion of these individuals will develop recurrent BC and will eventually succumb to the disease. Metastatic, therapy-resistant BC cells are refractory to cell death induced by multiple stresses. Here, we document that the vitamin D receptor (VDR) acts as a master transcriptional regulator of autophagy. Activation of the VDR by vitamin D induces autophagy and an autophagic transcriptional signature in BC cells that correlates with increased survival in patients; strikingly, this signature is present in the normal mammary gland and is progressively lost in patients with metastatic BC. A number of epidemiological studies have shown that sufficient vitamin D serum levels might be protective against BC. We observed that dietary vitamin D supplementation in mice increases basal levels of autophagy in the normal mammary gland, highlighting the potential of vitamin D as a cancer-preventive agent. These findings point to a role of vitamin D and the VDR in modulating autophagy and cell death in both the normal mammary gland and BC cells.

Shortsighted Evolution Constrains the Efficacy of Long-Term Bet Hedging.

To survive unpredictable environmental change, many organisms adopt bet-hedging strategies that are initially costly but provide a long-term fitness benefit. The temporal extent of these deferred fitness benefits determines whether bet-hedging organisms can survive long enough to realize them. In this article, we examine a model of microbial bet hedging in which there are two paths to extinction: unpredictable environmental change and demographic stochasticity. In temporally correlated environments, these drivers of extinction select for different switching strategies. Rapid phenotype switching ensures survival in the face of unpredictable environmental change, while slower-switching organisms become extinct. However, when both switching strategies are present in the same population, then demographic stochasticity-enforced by a limited population size-leads to extinction of the faster-switching organism. As a result, we find a novel form of evolutionary suicide whereby selection in a fluctuating environment can favor bet-hedging strategies that ultimately increase the risk of extinction. Population structures with multiple subpopulations and dispersal can reduce the risk of extinction from unpredictable environmental change and shift the balance so as to facilitate the evolution of slower-switching organisms.

Bistability in a metabolic network underpins the de novo evolution of colony switching in Pseudomonas fluorescens.

Phenotype switching is commonly observed in nature. This prevalence has allowed the elucidation of a number of underlying molecular mechanisms. However, little is known about how phenotypic switches arise and function in their early evolutionary stages. The first opportunity to provide empirical insight was delivered by an experiment in which populations of the bacterium Pseudomonas fluorescens SBW25 evolved, de novo, the ability to switch between two colony phenotypes. Here we unravel the molecular mechanism behind colony switching, revealing how a single nucleotide change in a gene enmeshed in central metabolism (carB) generates such a striking phenotype. We show that colony switching is underpinned by ON/OFF expression of capsules consisting of a colanic acid-like polymer. We use molecular genetics, biochemical analyses, and experimental evolution to establish that capsule switching results from perturbation of the pyrimidine biosynthetic pathway. Of central importance is a bifurcation point at which uracil triphosphate is partitioned towards either nucleotide metabolism or polymer production. This bifurcation marks a cell-fate decision point whereby cells with relatively high pyrimidine levels favour nucleotide metabolism (capsule OFF), while cells with lower pyrimidine levels divert resources towards polymer biosynthesis (capsule ON). This decision point is present and functional in the wild-type strain. Finally, we present a simple mathematical model demonstrating that the molecular components of the decision point are capable of producing switching. Despite its simple mutational cause, the connection between genotype and phenotype is complex and multidimensional, offering a rare glimpse of how noise in regulatory networks can provide opportunity for evolution.

Geometry shapes evolution of early multicellularity.

Organisms have increased in complexity through a series of major evolutionary transitions, in which formerly autonomous entities become parts of a novel higher-level entity. One intriguing feature of the higher-level entity after some major transitions is a division of reproductive labor among its lower-level units in which reproduction is the sole responsibility of a subset of units. Although it can have clear benefits once established, it is unknown how such reproductive division of labor originates. We consider a recent evolution experiment on the yeast Saccharomyces cerevisiae as a unique platform to address the issue of reproductive differentiation during an evolutionary transition in individuality. In the experiment, independent yeast lineages evolved a multicellular "snowflake-like" cluster formed in response to gravity selection. Shortly after the evolution of clusters, the yeast evolved higher rates of cell death. While cell death enables clusters to split apart and form new groups, it also reduces their performance in the face of gravity selection. To understand the selective value of increased cell death, we create a mathematical model of the cellular arrangement within snowflake yeast clusters. The model reveals that the mechanism of cell death and the geometry of the snowflake interact in complex, evolutionarily important ways. We find that the organization of snowflake yeast imposes powerful limitations on the available space for new cell growth. By dying more frequently, cells in clusters avoid encountering space limitations, and, paradoxically, reach higher numbers. In addition, selection for particular group sizes can explain the increased rate of apoptosis both in terms of total cell number and total numbers of collectives. Thus, by considering the geometry of a primitive multicellular organism we can gain insight into the initial emergence of reproductive division of labor during an evolutionary transition in individuality.

The contribution of private and public information in foraging by Australasian gannets.

Predators that forage on foods with temporally and spatially patchy distributions may rely on private or public sources of information to enhance their chances of foraging success. Using GPS tracking, field observations, and videography, we examined potential sites and mechanisms of information acquisition in departures for foraging trips by colonially breeding Australasian gannets (Morus serrator). Analyses of the bill-fencing ceremony between mated pairs of breeding gannets did not detect correlations between parameters of this reciprocal behavior and foraging trips, as would have been predicted if gannets used this behavior as a source of private information. Instead, 60% of the departing birds flew directly to join water rafts of other conspecific en route to the feeding grounds. The departure of solitary birds from the water rafts was synchronized (within 60 s) with the arrival of incoming foragers and also among departing birds. Furthermore, solitary departing birds from the rafts left in the same directional quadrant (90º slices) as the prior arriving (67%) and also prior departing forager (79%). When associated plunge dives of conspecific were visible from the colony, providing a public source of information, gannets more often departed from the water rafts in groups. Our study thus provides evidence for the use of water rafts, but not the nest site, as locations of information transfer, and also confirms the use of local enhancement as a strategy for foraging flights by Australasian gannets.

Emergence and maintenance of stable coexistence during a long-term multicellular evolution experiment.

The evolution of multicellular life spurred evolutionary radiations, fundamentally changing many of Earth's ecosystems. Yet little is known about how early steps in the evolution of multicellularity affect eco-evolutionary dynamics. Through long-term experimental evolution, we observed niche partitioning and the adaptive divergence of two specialized lineages from a single multicellular ancestor. Over 715 daily transfers, snowflake yeast were subjected to selection for rapid growth, followed by selection favouring larger group size. Small and large cluster-forming lineages evolved from a monomorphic ancestor, coexisting for over ~4,300 generations, specializing on divergent aspects of a trade-off between growth rate and survival. Through modelling and experimentation, we demonstrate that coexistence is maintained by a trade-off between organismal size and competitiveness for dissolved oxygen. Taken together, this work shows how the evolution of a new level of biological individuality can rapidly drive adaptive diversification and the expansion of a nascent multicellular niche, one of the most historically impactful emergent properties of this evolutionary transition.

A conceptual framework for the evolutionary origins of multicellularity.

The evolution of multicellular organisms from unicellular counterparts involved a transition in Darwinian individuality from single cells to groups. A particular challenge is to understand the nature of the earliest groups, the causes of their evolution, and the opportunities for emergence of Darwinian properties. Here we outline a conceptual framework based on a logical set of possible pathways for evolution of the simplest self-replicating groups. Central to these pathways is the recognition of a finite number of routes by which genetic information can be transmitted between individual cells and groups. We describe the form and organization of each primordial group state and consider factors affecting persistence and evolution of the nascent multicellular forms. Implications arising from our conceptual framework become apparent when attempting to partition fitness effects at individual and group levels. These are discussed with reference to the evolutionary emergence of individuality and its manifestation in extant multicellular life-including those of marginal Darwinian status.

Evolutionary consequences of nascent multicellular life cycles.

A key step in the evolutionary transition to multicellularity is the origin of multicellular groups as biological individuals capable of adaptation. Comparative work, supported by theory, suggests clonal development should facilitate this transition, although this hypothesis has never been tested in a single model system. We evolved 20 replicate populations of otherwise isogenic clonally reproducing 'snowflake' yeast (Δace2/∆ace2) and aggregative 'floc' yeast (GAL1p::FLO1 /GAL1p::FLO1) with daily selection for rapid growth in liquid media, which favors faster cell division, followed by selection for rapid sedimentation, which favors larger multicellular groups. While both genotypes adapted to this regime, growing faster and having higher survival during the group-selection phase, there was a stark difference in evolutionary dynamics. Aggregative floc yeast obtained nearly all their increased fitness from faster growth, not improved group survival; indicating that selection acted primarily at the level of cells. In contrast, clonal snowflake yeast mainly benefited from higher group-dependent fitness, indicating a shift in the level of Darwinian individuality from cells to groups. Through genome sequencing and mathematical modeling, we show that the genetic bottlenecks in a clonal life cycle also drive much higher rates of genetic drift-a result with complex implications for this evolutionary transition. Our results highlight the central role that early multicellular life cycles play in the process of multicellular adaptation.

Emergence and maintenance of stable coexistence during a long-term multicellular evolution experiment.

The evolution of multicellular life spurred evolutionary radiations, fundamentally changing many of Earth’s ecosystems. Yet little is known about how early steps in the evolution of multicellularity transform eco-evolutionary dynamics, e.g., via niche expansion processes that may facilitate coexistence. Using long-term experimental evolution in the snowflake yeast model system, we show that the evolution of multicellularity drove niche partitioning and the adaptive divergence of two distinct, specialized lineages from a single multicellular ancestor. Over 715 daily transfers, snowflake yeast were subject to selection for rapid growth in rich media, followed by selection favoring larger group size. Both small and large cluster-forming lineages evolved from a monomorphic ancestor, coexisting for over ~4,300 generations. These small and large sized snowflake yeast lineages specialized on divergent aspects of a trade-off between growth rate and survival, mirroring predictions from ecological theory. Through modeling and experimentation, we demonstrate that coexistence is maintained by a trade-off between organismal size and competitiveness for dissolved oxygen. Taken together, this work shows how the evolution of a new level of biological individuality can rapidly drive adaptive diversification and the expansion of a nascent multicellular niche, one of the most historically-impactful emergent properties of this evolutionary transition.

Exclusion rules, bottlenecks and the evolution of stochastic phenotype switching.

Stochastic phenotype switching--often considered a bet hedging or risk-reducing strategy--can enhance the probability of survival in fluctuating environments. A recent experiment provided direct evidence for an adaptive origin by showing the de novo evolution of switching in bacterial populations propagated under a selective regime that captured essential features of the host immune response. The regime involved strong frequency-dependent selection realized via dual imposition of an exclusion rule and population bottleneck. Applied at the point of transfer between environments, the phenotype common in the current environment was assigned a fitness of zero and was thus excluded from participating in the next round (the exclusion rule). In addition, also at the point of transfer, and so as to found the next bout of selection, a single phenotypically distinct type was selected at random from among the survivors (the bottleneck). Motivated by this experiment, we develop a mathematical model to explore the broader significance of key features of the selective regime. Through a combination of analytical and numerical results, we show that exclusion rules and population bottlenecks act in tandem as potent selective agents for stochastic phenotype switching, such that even when initially rare, and when switching engenders a cost in Malthusian fitness, organisms with the capacity to switch can invade non-switching populations and replace non-switching types. Simulations demonstrate the robustness of our findings to alterations in switching rate, fidelity of exclusion, bottleneck size, duration of environmental state and growth rate. We also demonstrate the relevance of our model to a range of biological scenarios such as bacterial persistence and the evolution of sex.

The evolutionary emergence of stochastic phenotype switching in bacteria.

Stochastic phenotype switching - or bet hedging - is a pervasive feature of living systems and common in bacteria that experience fluctuating (unpredictable) environmental conditions. Under such conditions, the capacity to generate variable offspring spreads the risk of being maladapted in the present environment, against offspring likely to have some chance of survival in the future. While a rich subject for theoretical studies, little is known about the selective causes responsible for the evolutionary emergence of stochastic phenotype switching. Here we review recent work - both theoretical and experimental - that sheds light on ecological factors that favour switching types over non-switching types. Of particular relevance is an experiment that provided evidence for an adaptive origin of stochastic phenotype switching by subjecting bacterial populations to a selective regime that mimicked essential features of the host immune response. Central to the emergence of switching types was frequent imposition of 'exclusion rules' and 'population bottlenecks' - two complementary faces of frequency dependent selection. While features of the immune response, exclusion rules and bottlenecks are likely to operate in many natural environments. Together these factors define a set of selective conditions relevant to the evolution of stochastic switching, including antigenic variation and bacterial persistence.

The Consequences of Budding versus Binary Fission on Adaptation and Aging in Primitive Multicellularity.

Early multicellular organisms must gain adaptations to outcompete their unicellular ancestors, as well as other multicellular lineages. The tempo and mode of multicellular adaptation is influenced by many factors including the traits of individual cells. We consider how a fundamental aspect of cells, whether they reproduce via binary fission or budding, can affect the rate of adaptation in primitive multicellularity. We use mathematical models to study the spread of beneficial, growth rate mutations in unicellular populations and populations of multicellular filaments reproducing via binary fission or budding. Comparing populations once they reach carrying capacity, we find that the spread of mutations in multicellular budding populations is qualitatively distinct from the other populations and in general slower. Since budding and binary fission distribute age-accumulated damage differently, we consider the effects of cellular senescence. When growth rate decreases with cell age, we find that beneficial mutations can spread significantly faster in a multicellular budding population than its corresponding unicellular population or a population reproducing via binary fission. Our results demonstrate that basic aspects of the cell cycle can give rise to different rates of adaptation in multicellular organisms.

Oxygen suppression of macroscopic multicellularity.

Atmospheric oxygen is thought to have played a vital role in the evolution of large, complex multicellular organisms. Challenging the prevailing theory, we show that the transition from an anaerobic to an aerobic world can strongly suppress the evolution of macroscopic multicellularity. Here we select for increased size in multicellular 'snowflake' yeast across a range of metabolically-available O2 levels. While yeast under anaerobic and high-O2 conditions evolved to be considerably larger, intermediate O2 constrained the evolution of large size. Through sequencing and synthetic strain construction, we confirm that this is due to O2-mediated divergent selection acting on organism size. We show via mathematical modeling that our results stem from nearly universal evolutionary and biophysical trade-offs, and thus should apply broadly. These results highlight the fact that oxygen is a double-edged sword: while it provides significant metabolic advantages, selection for efficient use of this resource may paradoxically suppress the evolution of macroscopic multicellular organisms.

The Grayness of the Origin of Life.

In the search for life beyond Earth, distinguishing the living from the non-living is paramount. However, this distinction is often elusive, as the origin of life is likely a stepwise evolutionary process, not a singular event. Regardless of the favored origin of life model, an inherent "grayness" blurs the theorized threshold defining life. Here, we explore the ambiguities between the biotic and the abiotic at the origin of life. The role of grayness extends into later transitions as well. By recognizing the limitations posed by grayness, life detection researchers will be better able to develop methods sensitive to prebiotic chemical systems and life with alternative biochemistries.

Ecological Advantages and Evolutionary Limitations of Aggregative Multicellular Development.

All multicellular organisms develop through one of two basic routes: they either aggregate from free-living cells, creating potentially chimeric multicellular collectives, or they develop clonally via mother-daughter cellular adhesion. Although evolutionary theory makes clear predictions about trade-offs between these developmental modes, these have never been experimentally tested in otherwise genetically identical organisms. We engineered unicellular baker's yeast (Saccharomyces cerevisiae) to develop either clonally ("snowflake"; Δace2) or aggregatively ("floc"; GAL1p::FLO1) and examined their fitness in a fluctuating environment characterized by periods of growth and selection for rapid sedimentation. When cultured independently, aggregation was far superior to clonal development, providing a 35% advantage during growth and a 2.5-fold advantage during settling selection. Yet when competed directly, clonally developing snowflake yeast rapidly displaced aggregative floc. This was due to unexpected social exploitation: snowflake yeast, which do not produce adhesive FLO1, nonetheless become incorporated into flocs at a higher frequency than floc cells themselves. Populations of chimeric clusters settle much faster than floc alone, providing snowflake yeast with a fitness advantage during competition. Mathematical modeling suggests that such developmental cheating may be difficult to circumvent; hypothetical "choosy floc" that avoid exploitation by maintaining clonality pay an ecological cost when rare, often leading to their extinction. Our results highlight the conflict at the heart of aggregative development: non-specific cellular binding provides a strong ecological advantage-the ability to quickly form groups-but this very feature leads to its exploitation.