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.

Spatial constraints and stochastic seeding subvert microbial arms race.

Surface attached communities of microbes grow in a wide variety of environments. Often, the size of these microbial community is constrained by their physical surroundings. However, little is known about how size constraints of a colony impact the outcome of microbial competitions. Here, we use individual-based models to simulate contact killing between two bacterial strains with different killing rates in a wide range of community sizes. We found that community size has a substantial impact on outcomes; in fact, in some competitions the identity of the most fit strain differs in large and small environments. Specifically, when at a numerical disadvantage, the strain with the slow killing rate is more successful in smaller environments than in large environments. The improved performance in small spaces comes from finite size effects; stochastic fluctuations in the initial relative abundance of each strain in small environments lead to dramatically different outcomes. However, when the slow killing strain has a numerical advantage, it performs better in large spaces than in small spaces, where stochastic fluctuations now aid the fast killing strain in small communities. Finally, we experimentally validate these results by confining contact killing strains of Vibrio cholerae in transmission electron microscopy grids. The outcomes of these experiments are consistent with our simulations. When rare, the slow killing strain does better in small environments; when common, the slow killing strain does better in large environments. Together, this work demonstrates that finite size effects can substantially modify antagonistic competitions, suggesting that colony size may, at least in part, subvert the microbial arms race.

Examining the role of oxygen-binding proteins on the early evolution of multicellularity.

Oxygen availability is a key factor in the evolution of multicellularity, as larger and more sophisticated organisms often require mechanisms allowing efficient oxygen delivery to their tissues. One such mechanism is the presence of oxygen-binding proteins, such as globins and hemerythrins, which arose in the ancestor of bilaterian animals. Despite their importance, the precise mechanisms by which oxygen-binding proteins influenced the early stages of multicellular evolution under varying environmental oxygen levels are not yet clear. We addressed this knowledge gap by heterologously expressing the oxygen binding proteins myoglobin and myohemerythrin in snowflake yeast, a model system of simple, undifferentiated multicellularity. These proteins increased the depth and rate of oxygen diffusion, increasing the fitness of snowflake yeast growing aerobically. Experiments show that, paradoxically, oxygen-binding proteins confer a greater fitness benefit for larger organisms under high, not low, O2 conditions. We show via biophysical modeling that this is because facilitated diffusion is more efficient when oxygen is abundant, transporting a greater quantity of O2 which can be used for metabolism. By alleviating anatomical diffusion limitations to oxygen consumption, the evolution of O2-binding proteins in the oxygen-rich Neoproterozoic may have been a key breakthrough enabling the evolution of increasingly large, complex multicellular metazoan lineages.

The biophysical basis of bacterial colony growth.

Bacteria often attach to surfaces and grow densely-packed communities called biofilms. As biofilms grow, they expand across the surface, increasing their surface area and access to nutrients. Thus, the overall growth rate of a biofilm is directly dependent on its "range expansion" rate. One factor that limits the range expansion rate is vertical growth; at the biofilm edge there is a direct trade-off between horizontal and vertical growth-the more a biofilm grows up, the less it can grow out. Thus, the balance of horizontal and vertical growth impacts the range expansion rate and, crucially, the overall biofilm growth rate. However, the biophysical connection between horizontal and vertical growth remains poorly understood, due in large part to difficulty in resolving biofilm shape with sufficient spatial and temporal resolution from small length scales to macroscopic sizes. Here, we experimentally show that the horizontal expansion rate of bacterial colonies is controlled by the contact angle at the biofilm edge. Using white light interferometry, we measure the three-dimensional surface morphology of growing colonies, and find that small colonies are surprisingly well-described as spherical caps. At later times, nutrient diffusion and uptake prevent the tall colony center from growing exponentially. However, the colony edge always has a region short enough to grow exponentially; the size and shape of this region, characterized by its contact angle, along with cellular doubling time, determines the range expansion rate. We found that the geometry of the exponentially growing biofilm edge is well-described as a spherical-cap-napkin-ring, i.e., a spherical cap with a cylindrical hole in its center (where the biofilm is too tall to grow exponentially). We derive an exact expression for the spherical-cap-napkin-ring-based range expansion rate; further, to first order, the expansion rate only depends on the colony contact angle, the thickness of the exponentially growing region, and the cellular doubling time. We experimentally validate both of these expressions. In line with our theoretical predictions, we find that biofilms with long cellular doubling times and small contact angles do in fact grow faster than biofilms with short cellular doubling times and large contact angles. Accordingly, sensitivity analysis shows that biofilm growth rates are more sensitive to their contact angles than to their cellular growth rates. Thus, to understand the fitness of a growing biofilm, one must account for its shape, not just its cellular doubling time.

Spontaneous Emergence of Multicellular Heritability.

The major transitions in evolution include events and processes that result in the emergence of new levels of biological individuality. For collectives to undergo Darwinian evolution, their traits must be heritable, but the emergence of higher-level heritability is poorly understood and has long been considered a stumbling block for nascent evolutionary transitions. Using analytical models, synthetic biology, and biologically-informed simulations, we explored the emergence of trait heritability during the evolution of multicellularity. Prior work on the evolution of multicellularity has asserted that substantial collective-level trait heritability either emerges only late in the transition or requires some evolutionary change subsequent to the formation of clonal multicellular groups. In a prior analytical model, we showed that collective-level heritability not only exists but is usually more heritable than the underlying cell-level trait upon which it is based, as soon as multicellular groups form. Here, we show that key assumptions and predictions of that model are borne out in a real engineered biological system, with important implications for the emergence of collective-level heritability.

De novo evolution of macroscopic multicellularity.

While early multicellular lineages necessarily started out as relatively simple groups of cells, little is known about how they became Darwinian entities capable of sustained multicellular evolution1-3. Here we investigate this with a multicellularity long-term evolution experiment, selecting for larger group size in the snowflake yeast (Saccharomyces cerevisiae) model system. Given the historical importance of oxygen limitation4, our ongoing experiment consists of three metabolic treatments5-anaerobic, obligately aerobic and mixotrophic yeast. After 600 rounds of selection, snowflake yeast in the anaerobic treatment group evolved to be macroscopic, becoming around 2 × 104 times larger (approximately mm scale) and about 104-fold more biophysically tough, while retaining a clonal multicellular life cycle. This occurred through biophysical adaptation-evolution of increasingly elongate cells that initially reduced the strain of cellular packing and then facilitated branch entanglements that enabled groups of cells to stay together even after many cellular bonds fracture. By contrast, snowflake yeast competing for low oxygen5 remained microscopic, evolving to be only around sixfold larger, underscoring the critical role of oxygen levels in the evolution of multicellular size. Together, this research provides unique insights into an ongoing evolutionary transition in individuality, showing how simple groups of cells overcome fundamental biophysical limitations through gradual, yet sustained, multicellular evolution.

Vertical growth dynamics of biofilms.

During the biofilm life cycle, bacteria attach to a surface and then reproduce, forming crowded, growing communities. Many theoretical models of biofilm growth dynamics have been proposed; however, difficulties in accurately measuring biofilm height across relevant time and length scales have prevented testing these models, or their biophysical underpinnings, empirically. Using white light interferometry, we measure the heights of microbial colonies with nanometer precision from inoculation to their final equilibrium height, producing a detailed empirical characterization of vertical growth dynamics. We propose a heuristic model for vertical growth dynamics based on basic biophysical processes inside a biofilm: diffusion and consumption of nutrients and growth and decay of the colony. This model captures the vertical growth dynamics from short to long time scales (10 min to 14 d) of diverse microorganisms, including bacteria and fungi.

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.

Varied solutions to multicellularity: The biophysical and evolutionary consequences of diverse intercellular bonds.

The diversity of multicellular organisms is, in large part, due to the fact that multicellularity has independently evolved many times. Nonetheless, multicellular organisms all share a universal biophysical trait: cells are attached to each other. All mechanisms of cellular attachment belong to one of two broad classes; intercellular bonds are either reformable or they are not. Both classes of multicellular assembly are common in nature, having independently evolved dozens of times. In this review, we detail these varied mechanisms as they exist in multicellular organisms. We also discuss the evolutionary implications of different intercellular attachment mechanisms on nascent multicellular organisms. The type of intercellular bond present during early steps in the transition to multicellularity constrains future evolutionary and biophysical dynamics for the lineage, affecting the origin of multicellular life cycles, cell-cell communication, cellular differentiation, and multicellular morphogenesis. The types of intercellular bonds used by multicellular organisms may thus result in some of the most impactful historical constraints on the evolution of multicellularity.

Evolution of a cis-Acting SNP That Controls Type VI Secretion in Vibrio cholerae.

Mutations in regulatory mechanisms that control gene expression contribute to phenotypic diversity and thus facilitate the adaptation of microbes and other organisms to new niches. Comparative genomics can be used to infer rewiring of regulatory architecture based on large effect mutations like loss or acquisition of transcription factors but may be insufficient to identify small changes in noncoding, intergenic DNA sequence of regulatory elements that drive phenotypic divergence. In human-derived Vibrio cholerae, the response to distinct chemical cues triggers production of multiple transcription factors that can regulate the type VI secretion system (T6), a broadly distributed weapon for interbacterial competition. However, to date, the signaling network remains poorly understood because no regulatory element has been identified for the major T6 locus. Here we identify a conserved cis-acting single nucleotide polymorphism (SNP) controlling T6 transcription and activity. Sequence alignment of the T6 regulatory region from diverse V. cholerae strains revealed conservation of the SNP that we rewired to interconvert V. cholerae T6 activity between chitin-inducible and constitutive states. This study supports a model of pathogen evolution through a noncoding cis-regulatory mutation and preexisting, active transcription factors that confers a different fitness advantage to tightly regulated strains inside a human host and unfettered strains adapted to environmental niches. IMPORTANCE Organisms sense external cues with regulatory circuits that trigger the production of transcription factors, which bind specific DNA sequences at promoters ("cis" regulatory elements) to activate target genes. Mutations of transcription factors or their regulatory elements create phenotypic diversity, allowing exploitation of new niches. Waterborne pathogen Vibrio cholerae encodes the type VI secretion system "nanoweapon" to kill competitor cells when activated. Despite identification of several transcription factors, no regulatory element has been identified in the promoter of the major type VI locus, to date. Combining phenotypic, genetic, and genomic analysis of diverse V. cholerae strains, we discovered a single nucleotide polymorphism in the type VI promoter that switches its killing activity between a constitutive state beneficial outside hosts and an inducible state for constraint in a host. Our results support a role for noncoding DNA in adaptation of this pathogen.

Cellular organization in lab-evolved and extant multicellular species obeys a maximum entropy law.

The prevalence of multicellular organisms is due in part to their ability to form complex structures. How cells pack in these structures is a fundamental biophysical issue, underlying their functional properties. However, much remains unknown about how cell packing geometries arise, and how they are affected by random noise during growth - especially absent developmental programs. Here, we quantify the statistics of cellular neighborhoods of two different multicellular eukaryotes: lab-evolved 'snowflake' yeast and the green alga Volvox carteri. We find that despite large differences in cellular organization, the free space associated with individual cells in both organisms closely fits a modified gamma distribution, consistent with maximum entropy predictions originally developed for granular materials. This 'entropic' cellular packing ensures a degree of predictability despite noise, facilitating parent-offspring fidelity even in the absence of developmental regulation. Together with simulations of diverse growth morphologies, these results suggest that gamma-distributed cell neighborhood sizes are a general feature of multicellularity, arising from conserved statistics of cellular packing.

Biomechanics of pollen pellet removal by the honey bee.

Honey bees (Apis mellifera) carry pollen back to their hive by mixing it with nectar and forming it into a pellet. The pellet must be firmly attached to their legs during flight, but also easily removable when deposited in the hive. How does the honey bee achieve these contrary aims? In this experimental study, we film honey bees removing pollen pellets and find they peel them off at speeds 2-10 times slower than their typical grooming speeds. Using a self-built pollen scraper, we find that slow removal speeds reduce the force and work required to remove the pellet under shear stress. Creep tests on individual pollen pellets revealed that pollen pellets are viscoelastic materials characterized by a Maxwell model with long relaxation times. The relaxation time enables the pellet to remain a solid during both transport and removal. We hope that this work inspires further research into viscoelastic materials in nature.

A New Contact Killing Toxin Permeabilizes Cells and Belongs to a Broadly Distributed Protein Family.

Vibrio cholerae is an aquatic Gram-negative bacterium that causes severe diarrheal cholera disease when ingested by humans. To eliminate competitor cells in both the external environment and inside hosts, V. cholerae uses the type VI secretion system (T6SS). The T6SS is a macromolecular contact-dependent weapon employed by many Gram-negative bacteria to deliver cytotoxic proteins into adjacent cells. In addition to canonical T6SS gene clusters encoded by all sequenced V. cholerae isolates, strain BGT49 encodes another locus, which we named auxiliary (Aux) cluster 4. The Aux 4 cluster is located on a mobile genetic element and can be used by killer cells to eliminate both V. cholerae and Escherichia coli cells in a T6SS-dependent manner. A putative toxin encoded in the cluster, which we name TpeV (type VI permeabilizing effector Vibrio), shares no homology to known proteins and does not contain motifs or domains indicative of function. Ectopic expression of TpeV in the periplasm of E. coli permeabilizes cells and disrupts the membrane potential. Using confocal microscopy, we confirm that susceptible target cells become permeabilized when competed with killer cells harboring the Aux 4 cluster. We also determine that tpiV, the gene located immediately downstream of tpeV, encodes an immunity protein that neutralizes the toxicity of TpeV. Finally, we show that TpeV homologs are broadly distributed across important human, animal, and plant pathogens and are localized in proximity to other T6SS genes. Our results suggest that TpeV is a toxin that belongs to a large family of T6SS proteins. IMPORTANCE Bacteria live in polymicrobial communities where competition for resources and space is essential for survival. Proteobacteria use the T6SS to eliminate neighboring cells and cause disease. However, the mechanisms by which many T6SS toxins kill or inhibit susceptible target cells are poorly understood. The sequence of the TpeV toxin that we describe here is unlike any previously described protein. We demonstrate that it has antimicrobial activity by permeabilizing cells, eliminating membrane potentials, and causing severe cytotoxicity. TpeV homologs are found near known T6SS genes in human, animal, and plant bacterial pathogens, indicating that the toxin is a representative member of a broadly distributed protein family. We propose that TpeV-like toxins contribute to the fitness of many bacteria. Finally, since antibiotic resistance is a critical global health threat, the discovery of new antimicrobial mechanisms could lead to the development of new treatments against resistant strains.

Glucose confers protection to Escherichia coli against contact killing by Vibrio cholerae.

Evolutionary arms races are broadly prevalent among organisms including bacteria, which have evolved defensive strategies against various attackers. A common microbial aggression mechanism is the type VI secretion system (T6SS), a contact-dependent bacterial weapon used to deliver toxic effector proteins into adjacent target cells. Sibling cells constitutively express immunity proteins that neutralize effectors. However, less is known about factors that protect non-sibling bacteria from T6SS attacks independently of cognate immunity proteins. In this study, we observe that human Escherichia coli commensal strains sensitive to T6SS attacks from Vibrio cholerae are protected when co-cultured with glucose. We confirm that glucose does not impair V. cholerae T6SS activity. Instead, we find that cells lacking the cAMP receptor protein (CRP), which regulates expression of hundreds of genes in response to glucose, survive significantly better against V. cholerae T6SS attacks even in the absence of glucose. Finally, we show that the glucose-mediated T6SS protection varies with different targets and killers. Our findings highlight the first example of an extracellular small molecule modulating a genetically controlled response for protection against T6SS attacks. This discovery may have major implications for microbial interactions during pathogen-host colonization and survival of bacteria in environmental communities.

Accumulation of dead cells from contact killing facilitates coexistence in bacterial biofilms.

Bacterial communities are governed by a wide variety of social interactions, some of which are antagonistic with potential significance for bacterial warfare. Several antagonistic mechanisms, such as killing via the type VI secretion system (T6SS), require killer cells to directly contact target cells. The T6SS is hypothesized to be a highly potent weapon, capable of facilitating the invasion and defence of bacterial populations. However, we find that the efficacy of contact killing is severely limited by the material consequences of cell death. Through experiments with Vibrio cholerae strains that kill via the T6SS, we show that dead cell debris quickly accumulates at the interface that forms between competing strains, preventing physical contact and thus preventing killing. While previous experiments have shown that T6SS killing can reduce a population of target cells by as much as 106-fold, we find that, as a result of the formation of dead cell debris barriers, the impact of contact killing depends sensitively on the initial concentration of killer cells. Killer cells are incapable of invading or eliminating competitors on a community level. Instead, bacterial warfare itself can facilitate coexistence between nominally antagonistic strains. While a variety of defensive strategies against microbial warfare exist, the material consequences of cell death provide target cells with their first line of defence.

Topological constraints in early multicellularity favor reproductive division of labor.

Reproductive division of labor (e.g. germ-soma specialization) is a hallmark of the evolution of multicellularity, signifying the emergence of a new type of individual and facilitating the evolution of increased organismal complexity. A large body of work from evolutionary biology, economics, and ecology has shown that specialization is beneficial when further division of labor produces an accelerating increase in absolute productivity (i.e. productivity is a convex function of specialization). Here we show that reproductive specialization is qualitatively different from classical models of resource sharing, and can evolve even when the benefits of specialization are saturating (i.e. productivity is a concave function of specialization). Through analytical theory and evolutionary individual-based simulations, we demonstrate that reproductive specialization is strongly favored in sparse networks of cellular interactions that reflect the morphology of early, simple multicellular organisms, highlighting the importance of restricted social interactions in the evolution of reproductive specialization.

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.

Cyberphysical risks of hacked internet-connected vehicles.

The integration of automotive technology with internet connectivity promises to both dramatically improve transportation while simultaneously introducing the potential for new unknown risks. Internet-connected vehicles are like digital data because they can be targeted for malicious hacking. Unlike digital data, however, internet-connected vehicles are cyberphysical systems that physically interact with each other and their environment. As such, the extension of cybersecurity concerns into the cyberphysical domain introduces new possibilities for self-organized phenomena in traffic flow. Here we study a scenario envisioned by cybersecurity experts leading to a large number of internet-connected vehicles being suddenly and simultaneously disabled. We investigate posthack traffic using agent-based simulations and discover the critical relevance of percolation for probabilistically predicting the outcomes on a multilane road in the immediate aftermath of a vehicle-targeted cyberattack. We develop an analytic percolation-based model to rapidly assess road conditions given the density of disabled vehicles and apply it to study the street network of Manhattan (New York City, New York, USA) revealing the city's vulnerability to this particular cyberphysical attack. While a comprehensive investigation of city-scale traffic around hacked vehicles is an extremely complicated problem, we find that the statistical physics of percolation can provide an estimate of the number of vehicles that critically disrupts citywide traffic flow. Our upper-bound estimate represents a quantification of citywide traffic disruptions when multiple vehicles are hacked.

Black Soldier Fly Larvae Rearrange under Compression.

Thousands of black soldier larvae hatch simultaneously from eggs laid within rotting vegetation or animal carcasses. Over the next few weeks, they grow while compressed by both their surroundings and each other. When compressed, these larvae rearrange to reduce the forces upon them. How quickly can larvae rearrange, and what final state do they choose? In this experimental study, we use a universal testing machine to conduct creep tests on larvae, squeezing them to set volume fractions and measuring the time course of their reaction force. Live larvae come to equilibrium at a rate 10 times faster than dead larvae, indicating that their small movements can rearrange them faster than just settling. The relaxation of dead larvae is well described by stretched exponentials, which also characterize hierarchical self-avoiding materials such as polymers or balls of crumpled aluminum foil. The equilibrium pressures of live larvae are comparable to those of dead larvae, suggesting that such pressures are dictated by the physics of their bodies rather than their behavior. Live larvae perform fluctuations to actively maintain this equilibrium pressure. This ability to survive large pressures might have applications in the larvae-rearing industry, where both live and dead larvae are packed in containers for shipping.

Drivers of Spatial Structure in Social Microbial Communities.

Microbes are social organisms, interacting primarily through secreted biomolecules. Many traits have evolved based solely on their effects upon other community members, yet even individually beneficial traits often create social side effects that are mediated by spatial population structure. Predicting the evolution of many microbial traits thus requires a comprehensive understanding of their social consequences. In this review, we examine the critical role of population spatial structure in microbial social evolution. We briefly review key mechanisms structuring microbial communities, focusing primarily on the universal roles of cellular death and reproduction. Finally, we explain how spatial assortment can be efficiently calculated in two-dimensional, surface-attached populations.