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.

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.

Structural hierarchy confers error tolerance in biological materials.

Structural hierarchy, in which materials possess distinct features on multiple length scales, is ubiquitous in nature. Diverse biological materials, such as bone, cellulose, and muscle, have as many as 10 hierarchical levels. Structural hierarchy confers many mechanical advantages, including improved toughness and economy of material. However, it also presents a problem: Each hierarchical level adds a new source of assembly errors and substantially increases the information required for proper assembly. This seems to conflict with the prevalence of naturally occurring hierarchical structures, suggesting that a common mechanical source of hierarchical robustness may exist. However, our ability to identify such a unifying phenomenon is limited by the lack of a general mechanical framework for structures exhibiting organization on disparate length scales. Here, we use simulations to substantiate a generalized model for the tensile stiffness of hierarchical filamentous networks with a nested, dilute triangular lattice structure. Following seminal work by Maxwell and others on criteria for stiff frames, we extend the concept of connectivity in network mechanics and find a similar dependence of material stiffness upon each hierarchical level. Using this model, we find that stiffness becomes less sensitive to errors in assembly with additional levels of hierarchy; although surprising, we show that this result is analytically predictable from first principles and thus potentially model independent. More broadly, this work helps account for the success of hierarchical, filamentous materials in biology and materials design and offers a heuristic for ensuring that desired material properties are achieved within the required tolerance.

One-pot system for synthesis, assembly, and display of functional single-span membrane proteins on oil-water interfaces.

Single-span membrane proteins (ssMPs) represent approximately one-half of all membrane proteins and play important roles in cellular communications. However, like all membrane proteins, ssMPs are prone to misfolding and aggregation because of the hydrophobicity of transmembrane helices, making them difficult to study using common aqueous solution-based approaches. Detergents and membrane mimetics can solubilize membrane proteins but do not always result in proper folding and functionality. Here, we use cell-free protein synthesis in the presence of oil drops to create a one-pot system for the synthesis, assembly, and display of functional ssMPs. Our studies suggest that oil drops prevent aggregation of some in vitro-synthesized ssMPs by allowing these ssMPs to localize on oil surfaces. We speculate that oil drops may provide a hydrophobic interior for cotranslational insertion of the transmembrane helices and a fluidic surface for proper assembly and display of the ectodomains. These functionalized oil drop surfaces could mimic cell surfaces and allow ssMPs to interact with cell surface receptors under an environment closest to cell-cell communication. Using this approach, we showed that apoptosis-inducing human transmembrane proteins, FasL and TRAIL, synthesized and displayed on oil drops induce apoptosis of cultured tumor cells. In addition, we take advantage of hydrophobic interactions of transmembrane helices to manipulate the assembly of ssMPs and create artificial clusters on oil drop surfaces. Thus, by coupling protein synthesis with self-assembly at the water-oil interface, we create a platform that can use recombinant ssMPs to communicate with cells.

Immotile Active Matter: Activity from Death and Reproduction.

Unlike equilibrium atomic solids, biofilms-soft solids composed of bacterial cells-do not experience significant thermal fluctuations at the constituent level. However, living cells stochastically reproduce and die, provoking a mechanical response. We investigate the mechanical consequences of cellular death and reproduction by measuring surface-height fluctuations of biofilms containing two mutually antagonistic strains of Vibrio cholerae that kill one another on contact via the type VI secretion system. While studies of active matter typically focus on activity via constituent mobility, here, activity is mediated by reproduction and death events in otherwise immobilized cells. Biofilm surface topography is measured in the nearly homeostatic limit via white light interferometry. Although biofilms are far from equilibrium systems, measured surface-height fluctuation spectra resemble the spectra of thermal permeable membranes but with an activity-mediated effective temperature, as predicted by Risler, Peilloux, and Prost [Phys. Rev. Lett. 115, 258104 (2015)PRLTAO0031-900710.1103/PhysRevLett.115.258104]. By comparing the activity of killer strains of V. cholerae with that of genetically modified strains that cannot kill each other and validating with individual-based simulations, we demonstrate that extracted effective temperatures increase with the amount of death and reproduction and that death and reproduction can fluidize biofilms. Together, these observations demonstrate the unique physical consequences of activity mediated by death and reproduction events.

Suppression of the coffee-ring effect by shape-dependent capillary interactions.

When a drop of liquid dries on a solid surface, its suspended particulate matter is deposited in ring-like fashion. This phenomenon, known as the coffee-ring effect, is familiar to anyone who has observed a drop of coffee dry. During the drying process, drop edges become pinned to the substrate, and capillary flow outward from the centre of the drop brings suspended particles to the edge as evaporation proceeds. After evaporation, suspended particles are left highly concentrated along the original drop edge. The coffee-ring effect is manifested in systems with diverse constituents, ranging from large colloids to nanoparticles and individual molecules. In fact--despite the many practical applications for uniform coatings in printing, biology and complex assembly-the ubiquitous nature of the effect has made it difficult to avoid. Here we show experimentally that the shape of the suspended particles is important and can be used to eliminate the coffee-ring effect: ellipsoidal particles are deposited uniformly during evaporation. The anisotropic shape of the particles significantly deforms interfaces, producing strong interparticle capillary interactions. Thus, after the ellipsoids are carried to the air-water interface by the same outward flow that causes the coffee-ring effect for spheres, strong long-ranged interparticle attractions between ellipsoids lead to the formation of loosely packed or arrested structures on the air-water interface. These structures prevent the suspended particles from reaching the drop edge and ensure uniform deposition. Interestingly, under appropriate conditions, suspensions of spheres mixed with a small number of ellipsoids also produce uniform deposition. Thus, particle shape provides a convenient parameter to control the deposition of particles, without modification of particle or solvent chemistry.

Measuring the Nonuniform Evaporation Dynamics of Sprayed Sessile Microdroplets with Quantitative Phase Imaging.

We demonstrate real-time quantitative phase imaging as a new optical approach for measuring the evaporation dynamics of sessile microdroplets. Quantitative phase images of various droplets were captured during evaporation. The images enabled us to generate time-resolved three-dimensional topographic profiles of droplet shape with nanometer accuracy and, without any assumptions about droplet geometry, to directly measure important physical parameters that characterize surface wetting processes. Specifically, the time-dependent variation of the droplet height, volume, contact radius, contact angle distribution along the droplet's perimeter, and mass flux density for two different surface preparations are reported. The studies clearly demonstrate three phases of evaporation reported previously: pinned, depinned, and drying modes; the studies also reveal instances of partial pinning. Finally, the apparatus is employed to investigate the cooperative evaporation of the sprayed droplets. We observe and explain the neighbor-induced reduction in evaporation rate, that is, as compared to predictions for isolated droplets. In the future, the new experimental methods should stimulate the exploration of colloidal particle dynamics on the gas-liquid-solid interface.

Physics in ordered and disordered colloidal matter composed of poly(N-isopropylacrylamide) microgel particles.

This review collects and describes experiments that employ colloidal suspensions to probe physics in ordered and disordered solids and related complex fluids. The unifying feature of this body of work is its clever usage of poly(N-isopropylacrylamide) (PNIPAM) microgel particles. These temperature-sensitive colloidal particles provide experimenters with a 'knob' for in situ control of particle size, particle interaction and particle packing fraction that, in turn, influence the structural and dynamical behavior of the complex fluids and solids. A brief summary of PNIPAM particle synthesis and properties is given, followed by a synopsis of current activity in the field. The latter discussion describes a variety of soft matter investigations including those that explore formation and melting of crystals and clusters, and those that probe structure, rearrangement and rheology of disordered (jammed/glassy) and partially ordered matter. The review, therefore, provides a snapshot of a broad range of physics phenomenology which benefits from the unique properties of responsive microgel particles.

Diffraction phase microscopy: monitoring nanoscale dynamics in materials science [invited].

Quantitative phase imaging (QPI) utilizes the fact that the phase of an imaging field is much more sensitive than its amplitude. As fields from the source interact with the specimen, local variations in the phase front are produced, which provide structural information about the sample and can be used to reconstruct its topography with nanometer accuracy. QPI techniques do not require staining or coating of the specimen and are therefore nondestructive. Diffraction phase microscopy (DPM) combines many of the best attributes of current QPI methods; its compact configuration uses a common-path off-axis geometry which realizes the benefits of both low noise and single-shot imaging. This unique collection of features enables the DPM system to monitor, at the nanoscale, a wide variety of phenomena in their natural environments. Over the past decade, QPI techniques have become ubiquitous in biological studies and a recent effort has been made to extend QPI to materials science applications. We briefly review several recent studies which include real-time monitoring of wet etching, photochemical etching, surface wetting and evaporation, dissolution of biodegradable electronic materials, and the expansion and deformation of thin-films. We also discuss recent advances in semiconductor wafer defect detection using QPI.

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.

Effects of particle shape on growth dynamics at edges of evaporating drops of colloidal suspensions.

We study the influence of particle shape on growth processes at the edges of evaporating drops. Aqueous suspensions of colloidal particles evaporate on glass slides, and convective flows during evaporation carry particles from drop center to drop edge, where they accumulate. The resulting particle deposits grow inhomogeneously from the edge in two dimensions, and the deposition front, or growth line, varies spatiotemporally. Measurements of the fluctuations of the deposition front during evaporation enable us to identify distinct growth processes that depend strongly on particle shape. Sphere deposition exhibits a classic Poisson-like growth process; deposition of slightly anisotropic particles, however, belongs to the Kardar-Parisi-Zhang universality class, and deposition of highly anisotropic ellipsoids appears to belong to a third universality class, characterized by Kardar-Parisi-Zhang fluctuations in the presence of quenched disorder.

Relationship between neighbor number and vibrational spectra in disordered colloidal clusters with attractive interactions.

We study connections between vibrational spectra and average nearest neighbor number in disordered clusters of colloidal particles with attractive interactions. Measurements of displacement covariances between particles in each cluster permit calculation of the stiffness matrix, which contains effective spring constants linking pairs of particles. From the cluster stiffness matrix, we derive vibrational properties of corresponding "shadow" glassy clusters, with the same geometric configuration and interactions as the "source" cluster but without damping. Here, we investigate the stiffness matrix to elucidate the origin of the correlations between the median frequency of cluster vibrational modes and average number of nearest neighbors in the cluster. We find that the mean confining stiffness of particles in a cluster, i.e., the ensemble-averaged sum of nearest neighbor spring constants, correlates strongly with average nearest neighbor number, and even more strongly with median frequency. Further, we find that the average oscillation frequency of an individual particle is set by the total stiffness of its nearest neighbor bonds; this average frequency increases as the square root of the nearest neighbor bond stiffness, in a manner similar to the simple harmonic oscillator.

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.

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.

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.

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.

Influence of particle shape on bending rigidity of colloidal monolayer membranes and particle deposition during droplet evaporation in confined geometries.

We investigate the influence of particle shape on the bending rigidity of colloidal monolayer membranes (CMMs) and on evaporative processes associated with these membranes. Aqueous suspensions of colloidal particles are confined between glass plates and allowed to evaporate. Confinement creates ribbonlike air-water interfaces and facilitates measurement and characterization of CMM geometry during drying. Interestingly, interfacial buckling events occur during evaporation. Extension of the description of buckled elastic membranes to our quasi-2D geometry enables the determination of the ratio of CMM bending rigidity to its Young's modulus. Bending rigidity increases with increasing particle anisotropy, and particle deposition during evaporation is strongly affected by membrane elastic properties. During drying, spheres are deposited heterogeneously, but ellipsoids are not. Apparently, increased bending rigidity reduces contact line bending and pinning and induces uniform deposition of ellipsoids. Surprisingly, suspensions of spheres doped with a small number of ellipsoids are also deposited uniformly.

Surfactant-induced Marangoni eddies alter the coffee-rings of evaporating colloidal drops.

The influence of the small ionic surfactant sodium dodecyl sulfate (SDS) on the evaporation of drying colloidal droplets is quantitatively investigated. The addition of SDS leads to a significantly more uniform deposition of colloidal particles after evaporation (i.e., the so-called "coffee-ring effect" is dramatically altered). We understand this phenomenon in the context of circulating radial Marangoni flows induced by the variation of SDS concentration along the air-water interface. Video microscopy permits the direct visualization of the colloidal particles involved in these flows, revealing a surprisingly stable "Marangoni eddy" that prevents particle deposition at the drop perimeter.

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.