Whole-genome duplication in the Multicellularity Long Term Evolution Experiment.

Whole-genome duplication (WGD) is widespread across eukaryotes and can promote adaptive evolution1-4. However, given the instability of newly-formed polyploid genomes5-7, understanding how WGDs arise in a population, persist, and underpin adaptations remains a challenge. Using our ongoing Multicellularity Long Term Evolution Experiment (MuLTEE)8, we show that diploid snowflake yeast (Saccharomyces cerevisiae) under selection for larger multicellular size rapidly undergo spontaneous WGD. From its origin within the first 50 days of the experiment, tetraploids persist for the next 950 days (nearly 5,000 generations, the current leading edge of our experiment) in ten replicate populations, despite being genomically unstable. Using synthetic reconstruction, biophysical modeling, and counter-selection experiments, we found that tetraploidy evolved because it confers immediate fitness benefits in this environment, by producing larger, longer cells that yield larger clusters. The same selective benefit also maintained tetraploidy over long evolutionary timescales, inhibiting the reversion to diploidy that is typically seen in laboratory evolution experiments. Once established, tetraploidy facilitated novel genetic routes for adaptation, playing a key role in the evolution of macroscopic multicellular size via the origin of evolutionarily conserved aneuploidy. These results provide unique empirical insights into the evolutionary dynamics and impacts of WGD, showing how it can initially arise due to its immediate adaptive benefits, be maintained by selection, and fuel long-term innovations by creating additional dimensions of heritable genetic variation.

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.

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.

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.

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.

Geometry, packing, and evolutionary paths to increased multicellular size.

The evolutionary transition to multicellularity transformed life on earth, heralding the evolution of large, complex organisms. Recent experiments demonstrated that laboratory-evolved multicellular "snowflake yeast" readily overcome the physical barriers that limit cluster size by modifying cellular geometry [Jacobeen et al., Nat. Phys. 14, 286 (2018)10.1038/s41567-017-0002-y]. However, it is unclear why this route to large size is observed, rather than an evolved increase in intercellular bond strength. Here, we use a geometric model of the snowflake yeast growth form to examine the geometric efficiency of increasing size by modifying geometry and bond strength. We find that changing geometry is a far more efficient route to large size than evolving increased intercellular adhesion. In fact, increasing cellular aspect ratio is on average ∼13 times more effective than increasing bond strength at increasing the number of cells in a cluster. Modifying other geometric parameters, such as the geometric arrangement of mother and daughter cells, also had larger effects on cluster size than increasing bond strength. Simulations reveal that as cells reproduce, internal stress in the cluster increases rapidly; thus, increasing bond strength provides diminishing returns in cluster size. Conversely, as cells become more elongated, cellular packing density within the cluster decreases, which substantially decreases the rate of internal stress accumulation. This suggests that geometrically imposed physical constraints may have been a key early selective force guiding the emergence of multicellular complexity.

Zwitterionic, cationic, and anionic fluorinated chemicals in aqueous film forming foam formulations and groundwater from U.S. military bases by nonaqueous large-volume injection HPLC-MS/MS.

A new analytical method was developed to quantify 26 newly-identified and 21 legacy (e.g. perfluoroalkyl carboxylates, perfluoroalkyl sulfonates, and fluorotelomer sulfonates) per and polyfluorinated alkyl substances (PFAS) in groundwater and aqueous film forming foam (AFFF) formulations. Prior to analysis, AFFF formulations were diluted into methanol and PFAS in groundwater were micro liquid-liquid extracted. Methanolic dilutions of AFFF formulations and groundwater extracts were analyzed by large-volume injection (900 µL) high-performance liquid chromatography tandem mass spectrometry. Orthogonal chromatography was performed using cation exchange (silica) and anion exchange (propylamine) guard columns connected in series to a reverse-phase (C18) analytical column. Method detection limits for PFAS in groundwater ranged from 0.71 ng/L to 67 ng/L, and whole-method accuracy ranged from 96% to 106% for analytes for which matched authentic analytical standards were available. For analytes without authentic analytical standards, whole-method accuracy ranged from 78 % to 144 %, and whole-method precision was less than 15 % relative standard deviation for all analytes. A demonstration of the method on groundwater samples from five military bases revealed eight of the 26 newly-identified PFAS present at concentrations up to 6900 ng/L. The newly-identified PFAS represent a minor fraction of the fluorinated chemicals in groundwater relative to legacy PFAS. The profiles of PFAS in groundwater differ from those found in fluorotelomer- and electrofluorination-based AFFF formulations, which potentially indicates environmental transformation of PFAS.

Multiple-site concerted proton-electron transfer reactions of hydrogen-bonded phenols are nonadiabatic and well described by semiclassical Marcus theory.

Photo-oxidations of hydrogen-bonded phenols using excited-state polyarenes are described to derive fundamental understanding of multiple-site concerted proton-electron transfer reactions (MS-CPET). Experiments have examined phenol bases having -CPh(2)NH(2), -Py, and -CH(2)Py groups ortho to the phenol hydroxyl group and tert-butyl groups in the 4,6-positions for stability (HOAr-NH(2), HOAr-Py, and HOAr-CH(2)Py, respectively; Py = pyridyl; Ph = phenyl). The photo-oxidations proceed by intramolecular proton transfer from the phenol to the pendent base concerted with electron transfer to the excited polyarene. For comparison, 2,4,6-(t)Bu(3)C(6)H(2)OH, a phenol without a pendent base and tert-butyl groups in the 2,4,6-positions, has also been examined. Many of these bimolecular reactions are fast, with rate constants near the diffusion limit. Combining the photochemical k(CPET) values with those from prior thermal stopped-flow kinetic studies gives data sets for the oxidations of HOAr-NH(2) and HOAr-CH(2)Py that span over 10(7) in k(CPET) and nearly 0.9 eV in driving force (ΔG(o)'). Plots of log(k(CPET)) vs ΔG(o)', including both excited-state anthracenes and ground state aminium radical cations, define a single Marcus parabola in each case. These two data sets are thus well described by semiclassical Marcus theory, providing a strong validation of the use of this theory for MS-CPET. The parabolas give λ(CPET) ≅ 1.15-1.2 eV and H(ab) ≅ 20-30 cm(-1). These experiments represent the most direct measurements of H(ab) for MS-CPET reactions to date. Although rate constants are available only up to the diffusion limit, the parabolas clearly peak well below the adiabatic limit of ca. 6 × 10(12) s(-1). Thus, this is a very clear demonstration that the reactions are nonadiabatic. The nonadiabatic character slows the reactions by a factor of ~45. Results for the oxidation of HOAr-Py, in which the phenol and base are conjugated, and for oxidation of 2,4,6-(t)Bu(3)C(6)H(2)OH, which lacks a base, show that both have substantially lower λ and larger pre-exponential terms. The implications of these results for MS-CPET reactions are discussed.