3D Hydrogel Encapsulation Regulates Nephrogenesis in Kidney Organoids

Stem cell-derived kidney organoids contain nephron segments that recapitulate morphological and functional aspects of the human kidney. However, directed differentiation protocols for kidney organoids are largely conducted using biochemical signals to control differentiation. Here, the hypothesis that mechanical signals regulate nephrogenesis is investigated in 3D culture by encapsulating kidney organoids within viscoelastic alginate hydrogels with varying rates of stress relaxation. Tubular nephron segments are significantly more convoluted in kidney organoids differentiated in encapsulating hydrogels when compared with those in suspension culture. Hydrogel viscoelasticity regulates the spatial distribution of nephron segments within the differentiating kidney organoids. Consistent with these observations, a particle-based computational model predicts that the extent of deformation of the hydrogel-organoid interface regulates the morphology of nephron segments. Elevated extracellular calcium levels in the culture medium, which can be impacted by the hydrogels, decrease the glomerulus-to-tubule ratio of nephron segments. These findings reveal that hydrogel encapsulation regulates nephron patterning and morphology and suggest that the mechanical microenvironment is an important design variable for kidney regenerative medicine.

A mechanochemical model recapitulates distinct vertebrate gastrulation modes

During vertebrate gastrulation, an embryo transforms from a layer of epithelial cells into a multilayered gastrula. This process requires the coordinated movements of hundreds to tens of thousands of cells, depending on the organism. In the chick embryo, patterns of actomyosin cables spanning several cells drive coordinated tissue flows. Here, we derive a minimal theoretical framework that couples actomyosin activity to global tissue flows. Our model predicts the onset and development of gastrulation flows in normal and experimentally perturbed chick embryos, mimicking different gastrulation modes as an active stress instability. Varying initial conditions and a parameter associated with active cell ingression, our model recapitulates distinct vertebrate gastrulation morphologies, consistent with recently published experiments in the chick embryo. Altogether, our results show how changes in the patterning of critical cell behaviors associated with different force-generating mechanisms contribute to distinct vertebrate gastrulation modes via a self-organizing mechanochemical process.

Learning to write with the fluid rope trick

Direct ink writing, a versatile method of 3D and 4D printing, requires the precise placement of a nozzle just above the print surface to prevent fluid instabilities that cause deviations from the prescribed print path. But what if one could harness the instability associated with the spontaneously folding or coiling of a thin stream of viscous fluid, i.e. use the "fluid rope trick" to write specified patterns on a substrate? Here we use Deep Reinforcement Learning to derive control strategies for the motion of the extruding nozzle and thus the fluid patterns that are deposited on the surface. The method proceeds by having a learner (nozzle) repeatedly interact with the environment (a viscous filament simulator), and improves its strategy using the results of this experience. We demonstrate the outcome of the learned control instructions using experiments to manipulate a falling viscous jet and create cursive writing patterns and Pollockian paintings on substrates.

Contractility-Induced Phase Separation in Active Solids

Experiments over many decades are suggestive that the combination of cellular contractility and active phase separation in cell-matrix composites can enable spatiotemporal patterning in multicellular tissues across scales. To characterize these phenomena, we provide a general theory that incorporates active cellular contractility into the classical Cahn-Hilliard-Larché model for phase separation in passive viscoelastic solids. Within this framework, we show how a homogeneous cell-matrix mixture can be destabilized by activity via either a pitchfork or Hopf bifurcation, resulting in stable phase separation and/or traveling waves. Numerical simulations of the full equations allow us to track the evolution of the resulting self-organized patterns in periodic and mechanically constrained domains, and in different geometries. Altogether, our study underscores the importance of integrating both cellular activity and mechanical phase separation in understanding patterning in soft, active biosolids in both in vivo and in vitro settings.

Beak morphometry and morphogenesis across avian radiations

Adaptive avian radiations associated with the diversification of bird beaks into a multitude of forms enabling different functions are exemplified by Darwin's finches and Hawaiian honeycreepers. To elucidate the nature of these radiations, we quantified beak shape and skull shape using a variety of geometric measures that allowed us to collapse the variability of beak shape into a minimal set of geometric parameters. Furthermore, we find that just two measures of beak shape-the ratio of the width to length and the normalized sharpening rate (increase in the transverse beak curvature near the tip relative to that at the base of the beak)-are strongly correlated with diet. Finally, by considering how transverse sections to the beak centreline evolve with distance from the tip, we show that a simple geometry-driven growth law termed 'modified mean curvature flow' captures the beak shapes of Darwin's finches and Hawaiian honeycreepers. A surprising consequence of the simple growth law is that beak shapes that are not allowed based on the developmental programme of the beak are also not observed in nature, suggesting a link between evolutionary morphology and development in terms of growth-driven developmental constraints.

Evaporation-Driven Cellular Patterns in Confined Hyperelastic Hydrogels

When a hyperelastic hydrogel confined between two parallel glass plates begins to dry from a lateral boundary, the volume lost by evaporation is accommodated by an inward displacement of the air-hydrogel interface that induces an elastic deformation of the hydrogel. Once a critical front displacement is reached, we observe intermittent fracture events initiated by a geometric instability resulting in localized bursts at the interface. These bursts relax the stresses and irreversibly form air cavities that lead to cellular networks. We show that the spatial extent of the strain field prior to a burst, influenced by the air-hydrogel interfacial tension and the confinement of the gel, determines the characteristic size of the cavities.

Statics and diffusive dynamics of surfaces driven by p-atic topological defects

Inspired by epithelial morphogenesis, we consider a minimal model for the shaping of a surface driven by p-atic topological defects. We show that a positive (negative) defect can dynamically generate a (hyperbolic) cone whose shape evolves diffusively, and predict that a defect of charge + 1/p leads to a final semi-cone angle β which satisfies the inequality . By exploiting the fact that for axisymmetric surfaces, the extrinsic geometry is tightly coupled to the intrinsic geometry, we further show that the resulting stationary shape of a membrane with negligible bending modulus and embedded polar order is a deformed lemon with two defects at antipodal points. Finally, we close by pointing out that our results may be relevant beyond epithelial morphogenesis in such contexts as shape transitions in macroscopic closed spheroidal surfaces such as pollen grains.

How to Grow a Flat Leaf

Growing a flat lamina such as a leaf is almost impossible without some feedback to stabilize long wavelength modes that are easy to trigger since they are energetically cheap. Here we combine the physics of thin elastic plates with feedback control theory to explore how a leaf can remain flat while growing. We investigate both in-plane (metric) and out-of-plane (curvature) growth variation and account for both local and nonlocal feedback laws. We show that a linearized feedback theory that accounts for both spatially nonlocal and temporally delayed effects suffices to suppress long wavelength fluctuations effectively and explains recently observed statistical features of growth in tobacco leaves. Our work provides a framework for understanding the regulation of the shape of leaves and other leaflike laminar objects.

Ethics Inside the Black Box: Integrating Science and Technology Studies into Engineering and Public Policy Curricula

There is growing need for hybrid curricula that integrate constructivist methods from Science and Technology Studies (STS) into both engineering and policy courses at the undergraduate and graduate levels. However, institutional and disciplinary barriers have made implementing such curricula difficult at many institutions. While several programs have recently been launched that mix technical training with consideration of "societal" or "ethical issues," these programs often lack a constructivist element, leaving newly-minted practitioners entering practical fields ill-equipped to unpack the politics of knowledge and technology or engage with skeptical publics. This paper presents a novel format for designing interdisciplinary coursework that combines conceptual content from STS with training in engineering and policy. Courses following this format would ideally be team taught by instructors with advanced training in diverse fields, and hence co-learning between instructors and disciplines is a key element of the format. Several instruments for facilitating both student and instructor collaborative learning are introduced. The format is also designed for versatility: in addition to being adaptable to both technical and policy training environments, topics are modularized around a conceptual core so that issues ranging from biotech to nuclear security can be incorporated to fit programmatic needs and resources.

An additive framework for kirigami design

We present an additive approach for the inverse design of kirigami-based mechanical metamaterials by focusing on the empty (negative) spaces instead of the solid tiles. By considering each negative space as a four-bar linkage, we identify a simple recursive relationship between adjacent linkages, yielding an efficient method for creating kirigami patterns. This allows us to solve the kirigami design problem using elementary linear algebra, with compatibility, reconfigurability and rigid-deployability encoded into an iterative procedure involving simple matrix multiplications. The resulting linear design strategy circumvents the solution of a non-convex global optimization problem and allows us to control the degrees of freedom in the deployment angle field, linkage offsets and boundary conditions. We demonstrate this by creating a large variety of rigid-deployable, compact, reconfigurable kirigami patterns. We then realize our kirigami designs physically using two simple but effective fabrication strategies with very different materials. Altogether, our additive approaches present routes for efficient mechanical metamaterial design and fabrication based on ori/kirigami art forms.

Direct force measurement and loading on developing tissues in intact avian embryos

Developmental morphogenesis is driven by tissue stresses acting on tissue rheology. Direct measurements of forces in small tissues (100 µm-1 mm) in situ, such as in early embryos, require high spatial precision and minimal invasiveness. Here, we introduce a control-based approach, tissue force microscopy (TiFM), that integrates a mechanical cantilever probe and live imaging with closed-loop feedback control of mechanical loading in early chicken embryos. By testing previously qualitatively characterized force-producing tissues in the elongating body axis, we show that TiFM quantitatively captures stress dynamics with high sensitivity. TiFM also provides the means to apply stable, minimally invasive and physiologically relevant loads to drive tissue deformation and to follow the resulting morphogenetic progression associated with large-scale cell movements. Together, TiFM allows us to control tissue force measurement and manipulation in small developing embryos, and promises to contribute to the quantitative understanding of complex multi-tissue mechanics during development.

Explosive rigidity percolation in kirigami

Controlling the connectivity and rigidity of kirigami, i.e. the process of cutting paper to deploy it into an articulated system, is critical in the manifestations of kirigami in art, science and technology, as it provides the resulting metamaterial with a range of mechanical and geometric properties. Here, we combine deterministic and stochastic approaches for the control of rigidity in kirigami using the power of k choices, an approach borrowed from the statistical mechanics of explosive percolation transitions. We show that several methods for rigidifying a kirigami system by incrementally changing either the connectivity or the rigidity of individual components allow us to control the nature of the explosive transition by a choice of selection rules. Our results suggest simple lessons for the design of mechanical metamaterials.

Combined measures of mimetic fidelity explain imperfect mimicry in a brood parasite-host system

The persistence of imperfect mimicry in nature presents a challenge to mimicry theory. Some hypotheses for the existence of imperfect mimicry make differing predictions depending on how mimetic fidelity is measured. Here, we measure mimetic fidelity in a brood parasite-host system using both trait-based and response-based measures of mimetic fidelity. Cuckoo finches Anomalospiza imberbis lay imperfectly mimetic eggs that lack the fine scribbling characteristic of eggs of the tawny-flanked prinia Prinia subflava, a common host species. A trait-based discriminant analysis based on Minkowski functionals-that use geometric and topological morphometric methods related to egg pattern shape and coverage-reflects this consistent difference between host and parasite eggs. These methods could be applied to quantify other phenotypes including stripes and waved patterns. Furthermore, by painting scribbles onto cuckoo finch eggs and testing their rate of rejection compared to control eggs (i.e. a response-based approach to quantify mimetic fidelity), we show that prinias do not discriminate between eggs based on the absence of scribbles. Overall, our results support relaxed selection on cuckoo finches to mimic scribbles, since prinias do not respond differently to eggs with and without scribbles, despite the existence of this consistent trait difference.

Matrix viscoelasticity controls spatiotemporal tissue organization

Biomolecular and physical cues of the extracellular matrix environment regulate collective cell dynamics and tissue patterning. Nonetheless, how the viscoelastic properties of the matrix regulate collective cell spatial and temporal organization is not fully understood. Here we show that the passive viscoelastic properties of the matrix encapsulating a spheroidal tissue of breast epithelial cells guide tissue proliferation in space and in time. Matrix viscoelasticity prompts symmetry breaking of the spheroid, leading to the formation of invading finger-like protrusions, YAP nuclear translocation and epithelial-to-mesenchymal transition both in vitro and in vivo in a Arp2/3-complex-dependent manner. Computational modelling of these observations allows us to establish a phase diagram relating morphological stability with matrix viscoelasticity, tissue viscosity, cell motility and cell division rate, which is experimentally validated by biochemical assays and in vitro experiments with an intestinal organoid. Altogether, this work highlights the role of stress relaxation mechanisms in tissue growth dynamics, a fundamental process in morphogenesis and oncogenesis.

Rheomergy : collective behaviour mediated by active flow-based recruitment

The physics of signal propagation in a collection of organisms that communicate with each other both enables and limits how active excitations at the individual level reach, recruit and lead to collective patterning. Inspired by the spatio-temporal patterns in a planar swarm of bees that use pheromones and fanning flows to recruit additional bees, we develop a theoretical framework for patterning via active flow-based recruitment. Our model generalizes the well-known Patlak–Keller–Segel model of diffusion dominated aggregation and leads to an enhanced phase space of patterns spanned by two dimensionless parameters that measure the scaled stimulus/activity and the scaled chemotactic response. Together these determine the efficacy of signal communication via fluid flow (which we dub rheomergy ) that leads to a variety of migration and aggregation patterns, consistent with observations.

An allometric prior enhances acoustic niche partitioning signal

The acoustic niche hypothesis suggests that vocal signals of sympatric animal species are structured so as to minimize acoustic interference and facilitate communication. Accordingly, each species attempts to establish its own acoustic bandwidth so that intra-species signals are not masked. Detecting a non-random partitioning of the frequency spectrum among sympatric species could constitute evidence for the existence of acoustic avoidance behaviour. However, results from previous studies have been mixed or inconclusive, possibly as a consequence of overlooking the importance of physiological and ecological constraints. Here we introduce an improved test that incorporates prior information on body mass to account for the allometric correlation between mass (size) and vocalization frequency. By correcting for the bias induced by this correlation, the new test uncovers evidence of acoustic niche partitioning as a function of frequency in several tropical bird communities that would not be detected under a more standard test. Separately, we introduce a spatial version of the acoustic partitioning test which, in theory, could prove effective when data are collected from multiple sites located in close spatial proximity.

A continuous statistical-geometric framework for normative and impaired gaits

A quantitative analysis of human gait patterns in space-time provides an opportunity to observe variability within and across individuals of varying motor capabilities. Impaired gait significantly affects independence and quality of life, and thus a large part of clinical research is dedicated to improving gait through rehabilitative therapies. Evaluation of these paradigms relies on understanding the characteristic differences in the kinematics and underlying biomechanics of impaired and unimpaired locomotion, which has motivated quantitative measurement and analysis of the gait cycle. Previous analysis has largely been limited to a statistical comparison of manually selected pointwise metrics identified through expert knowledge. Here, we use a recent statistical-geometric framework, elastic functional data analysis (FDA), to decompose kinematic data into continuous 'amplitude' (spatial) and 'phase' (temporal) components, which can then be integrated with established dimensionality reduction techniques. We demonstrate the utility of elastic FDA through two unsupervised applications to post-stroke gait datasets. First, we distinguish between unimpaired, paretic and non-paretic gait presentations. Then, we use FDA to reveal robust, interpretable groups of differential response to exosuit assistance. The proposed methods aim to benefit clinical practice for post-stroke gait rehabilitation, and more broadly, to automate the quantitative analysis of motion.

Active Nematic Defects and Epithelial Morphogenesis

Inspired by recent experiments that highlight the role of nematic defects in the morphogenesis of epithelial tissues, we develop a minimal framework to study the dynamics of an active curved surface driven by its nematic texture. Allowing the surface to evolve via relaxational dynamics leads to a theory linking nematic defect dynamics, cellular division rates, and Gaussian curvature. Regions of large positive (negative) curvature and positive (negative) growth are colocalized with the presence of positive (negative) defects. In an ex-vivo setting of cultured murine neural progenitor cells, we show that our framework is consistent with the observed cell accumulation at positive defects and depletion at negative defects. In an in-vivo setting, we show that the defect configuration consisting of a bound +1 defect state, which is stabilized by activity, surrounded by two -1/2 defects can create a stationary ring configuration of tentacles, consistent with observations of a basal marine invertebrate Hydra.

Modular representation and control of floppy networks

Geometric graph models of systems as diverse as proteins, DNA assemblies, architected materials and robot swarms are useful abstract representations of these objects that also unify ways to study their properties and control them in space and time. While much work has been done in the context of characterizing the behaviour of these networks close to critical points associated with bond and rigidity percolation, isostaticity, etc., much less is known about floppy, underconstrained networks that are far more common in nature and technology. Here, we combine geometric rigidity and algebraic sparsity to provide a framework for identifying the zero energy floppy modes via a representation that illuminates the underlying hierarchy and modularity of the network and thence the control of its nestedness and locality. Our framework allows us to demonstrate a range of applications of this approach that include robotic reaching tasks with motion primitives, and predicting the linear and nonlinear response of elastic networks based solely on infinitesimal rigidity and sparsity, which we test using physical experiments. Our approach is thus likely to be of use broadly in dissecting the geometrical properties of floppy networks using algebraic sparsity to optimize their function and performance.

Optimal policies for mitigating pandemic costs: a minimal model

There have been a number of pharmaceutical and non-pharmaceutical interventions associated with COVID-19 over the past two years. Of the various non-pharmaceutical interventions that were proposed and implemented to control the spread of the COVID-19 pandemic partial and complete lockdowns were used repeatedly in an attempt to minimize the costs associated with mortality, economic losses and social factors, while being subject to constraints such as finite hospital capacity. Here, we use a minimal model to understand the costs and benefits of these strategies that mitigate pandemic costs subject to constraints, we adopt the language of optimal control theory. This allows us to determine top-down policies for the nature and dynamics of social contact rates given an age-structured model for the dynamics of the disease. Depending on the relative weights allocated to mortality and socioeconomic losses, we see that the optimal strategies range from long-term social-distancing only for the most vulnerable, to partial lockdown to ensure not over-running hospitals, to alternating-shifts with significant reduction in mortality and/or socioeconomic losses. Crucially, commonly used strategies that involve long periods of broad lockdown are almost never optimal, as they are highly unstable to reopening {and entail high socioeconomic costs}. Using parameter estimates from data available for Germany and the USA early in the pandemic, we quantify these policies and use sensitivity analysis in the relevant model parameters and initial conditions to determine the range of robustness of our policies. Finally we also discuss how bottom-up behavioral changes affect the dynamics of the pandemic and show they can work in tandem with top-down control policies to mitigate pandemic costs even more effectively.

Safeguarding human–wildlife cooperation

Human–wildlife cooperation occurs when humans and free‐living wild animals actively coordinate their behavior to achieve a mutually beneficial outcome. These interactions provide important benefits to both the human and wildlife communities involved, have wider impacts on the local ecosystem, and represent a unique intersection of human and animal cultures. The remaining active forms are human–honeyguide and human–dolphin cooperation, but these are at risk of joining several inactive forms (including human–wolf and human–orca cooperation). Human–wildlife cooperation faces a unique set of conservation challenges, as it requires multiple components—a motivated human and wildlife partner, a suitable environment, and compatible interspecies knowledge—which face threats from ecological and cultural changes. To safeguard human–wildlife cooperation, we recommend: (i) establishing ethically sound conservation strategies together with the participating human communities; (ii) conserving opportunities for human and wildlife participation; (iii) protecting suitable environments; (iv) facilitating cultural transmission of traditional knowledge; (v) accessibly archiving Indigenous and scientific knowledge; and (vi) conducting long‐term empirical studies to better understand these interactions and identify threats. Tailored safeguarding plans are therefore necessary to protect these diverse and irreplaceable interactions. Broadly, our review highlights that efforts to conserve biological and cultural diversity should carefully consider interactions between human and animal cultures.

Please see AfricanHoneyguides.com/abstract‐translations for Kiswahili and Portuguese translations of the abstract.

Prestrain-induced contraction in one-dimensional random elastic chains

Prestrained elastic networks arise in a number of biological and technological systems ranging from the cytoskeleton of cells to tensegrity structures. Motivated by this observation, we here consider a minimal model in one dimension to set the stage for understanding the response of such networks as a function of the prestrain. To this end we consider a chain [one-dimensional (1D) network] of elastic springs upon which a random, zero mean, finite variance prestrain is imposed. Numerical simulations and analytical predictions quantify the magnitude of the contraction as a function of the variance of the prestrain, and show that the chain always shrinks. To test these predictions, we vary the topology of the chain, consider more complex connectivity and show that our results are relatively robust to these changes.

Combing a double helix

Combing hair involves brushing away the topological tangles in a collective curl, defined as a bundle of interacting elastic filaments. Using a combination of experiment and computation, we study this problem that naturally links topology, geometry and mechanics. Observations show that the dominant interactions in hair are those of a two-body nature, corresponding to a braided homochiral double helix. This minimal model allows us to study the detangling of an elastic double helix driven by a single stiff tine that moves along it and leaves two untangled filaments in its wake. Our results quantify how the mechanics of detangling correlates with the dynamics of a topological quantity, the link density, that propagates ahead of the tine and flows out the free end as a link current. This in turn provides a measure of the maximum characteristic length of a single combing stroke in the many-body problem on a head of hair, producing an optimal combing strategy that balances trade-offs between comfort, efficiency and speed of combing in hair curls of varying geometrical and topological complexity.

Rectified random cell motility as a mechanism for embryo elongation

The body of vertebrate embryos forms by posterior elongation from a terminal growth zone called the tail bud. The tail bud is a source of highly motile cells that eventually constitute the presomitic mesoderm (PSM), a tissue that plays an important role in elongation movements. PSM cells establish an anterior-posterior cell motility gradient that parallels a gradient associated with the degradation of a specific cellular signal (FGF) known to be implicated in cell motility. Here, we combine the electroporation of fluorescent reporters in the PSM with time-lapse imaging in the chicken embryo to quantify cell diffusive movements along the motility gradient. We show that a simple microscopic model for random cell motility induced by FGF activity along with geometric confinement leads to rectified tissue elongation consistent with our observations. A continuum analog of the microscopic model leads to a macroscopic mechano-chemical model for tissue extension that couples FGF activity-induced cell motility and tissue rheology, and is consistent with the experimentally observed speed and extent of elongation. Together, our experimental observations and theoretical models explain how the continuous addition of cells at the tail bud combined with lateral confinement can be converted into oriented movement and drive body elongation.

A buckling-sheet ring oscillator for electronics-free, multimodal locomotion

Locomotion of soft robots typically relies on control of multiple inflatable actuators by electronic computers and hard valves. Soft pneumatic oscillators can reduce the demand on controllers by generating complex movements required for locomotion from a single, constant input pressure, but either have been constrained to low rates of flow of air or have required complex fabrication processes. Here, we describe a pneumatic oscillator fabricated from flexible, but inextensible, sheets that provides high rates of airflow for practical locomotion by combining three instabilities: out-of-plane buckling of the sheets, kinking of tubing attached to the sheets, and a system-level instability resulting from connection of an odd number of pneumatic inverters made from these sheets in a loop. This device, which we call a "buckling-sheet ring oscillator" (BRO), directly generates movement from its own interaction with its surroundings and consists only of readily available materials assembled in a simple process-specifically, stacking acetate sheets, nylon film, and double-sided tape, and attaching an elastomeric tube. A device incorporating a BRO is capable of both translational and rotational motion over varied terrain (even without a tether) and can climb upward against gravity and downward against the buoyant force encountered under water.

Euclid's Random Walk: Developmental Changes in the Use of Simulation for Geometric Reasoning

Euclidean geometry has formed the foundation of architecture, science, and technology for millennia, yet the development of human's intuitive reasoning about Euclidean geometry is not well understood. The present study explores the cognitive processes and representations that support the development of humans' intuitive reasoning about Euclidean geometry. One-hundred-twenty-five 7- to 12-year-old children and 30 adults completed a localization task in which they visually extrapolated missing parts of fragmented planar triangles and a reasoning task in which they answered verbal questions about the general properties of planar triangles. While basic Euclidean principles guided even young children's visual extrapolations, only older children and adults reasoned about triangles in ways that were consistent with Euclidean geometry. Moreover, a relation beteen visual extrapolation and reasoning appeared only in older children and adults. Reasoning consistent with Euclidean geometry may thus emerge when children abandon incorrect, axiomatic-based reasoning strategies and come to reason using mental simulations of visual extrapolations.

Geometrical dynamics of edge-driven accretive surface growth

Accretion of mineralized thin wall-like structures via localized growth along their edges is observed in physical and biological systems ranging from molluscan and brachiopod shells to carbonate–silica composite precipitates. To understand the shape of these mineralized structures, we develop a mathematical framework that treats the thin-walled shells as a smooth surface left in the wake of the growth front that can be described as an evolving space curve. Our theory then takes an explicit geometric form for the prescription of the velocity of the growth front curve, along with compatibility relations and a closure equation related to the nature of surface curling. Solutions of these equations capture a range of geometric precipitate patterns seen in abiotic and biotic forms across scales. In addition to providing a framework for the growth and form of these thin-walled morphologies, our theory suggests a new class of dynamical systems involving moving space curves that are compatible with non-Euclidean embeddings of surfaces.

Dynamics of cooperative excavation in ant and robot collectives

The solution of complex problems by the collective action of simple agents in both biologically evolved and synthetically engineered systems involves cooperative action. Understanding the resulting emergent solutions requires integrating across the organismal behavior of many individuals. Here, we investigate an ecologically relevant collective task in black carpenter ants Camponotus pennsylvanicus: excavation of a soft, erodible confining corral. These ants show a transition from individual exploratory excavation at random locations to spatially localized collective exploitative excavation and escape from the corral. Agent-based simulations and a minimal continuum theory that coarse-grains over individual actions and considers their integrated influence on the environment leads to the emergence of an effective phase space of behaviors, characterized in terms of excavation strength and cooperation intensity. To test the theory over the range of both observed and predicted behaviors, we use custom-built robots (RAnts) that respond to stimuli to characterize the phase space of emergence (and failure) of cooperative excavation. Tuning the amount of cooperation between RAnts, allows us to vary the efficiency of excavation and synthetically generate the entire range of macroscopic phases predicted by our theory. Overall, our approach shows how the cooperative completion of tasks can arise from simple rules that involve the interaction of agents with a dynamically changing environment that serves as both an enabler and a modulator of behavior.

Thermoregulatory morphodynamics of honeybee swarm clusters

During reproductive swarming, honeybees clusters of more than 10,000 individuals that hang from structures in the environment (e.g., tree branches) are exposed to diurnal variations in ambient temperature for up to a week. Swarm clusters collectively modulate their morphology in response to these variations (i.e., expanding/contracting in response to heating/cooling) to maintain their internal temperature within a tolerable range and to avoid exhausting their honey stores prematurely. To understand the spatiotemporal aspects of thermoregulatory morphing, we measured the change in size, shape and internal temperature profiles of swarm clusters in response to dynamic temperature ramp perturbations. We see that swarm clusters show a two-fold variation in their volume/density when heated from 15°C to 30°C. However, they do not reach an equilibrium size or shape when held at 30°C for 5 hours, long after the core temperature of the cluster has stabilized. Furthermore, the changes in cluster shape and size are hysteretic, contracting in response to cooling faster than expanding in response to heating. Although the base contact diameter of the cluster increased continuously when the swarm is heated, the change in length of the swarm (base totip) over time is non-monotonic. Consequently, the aspect ratio of the swarm fluctuated continuously even when held at a constant temperature. Taken together, our results quantify the hysteretic and anisotropic morphological responses of swarm clusters to ambient temperature variations while suggesting that both mechanical constraints and heat transfer govern their thermoregulatory morphodynamics.

Instabilities and patterns in a submerged jelling jet

When a thin stream of aqueous sodium alginate is extruded into a reacting calcium chloride bath, it polymerizes into a soft elastic tube that spontaneously forms helical coils due to the ambient fluid drag. We quantify the onset of this drag-induced instability and its nonlinear evolution using experiments, and explain the results using a combination of scaling, theory and simulations. By co-extruding a second (internal) liquid within the aqueous sodium alginate jet and varying the diameter of the jet and the rates of the co-extrusion of the two liquids, we show that we can tune the local composition of the composite filament and the nature of the ensuing instabilities to create soft filaments of variable relative buoyancy, shape and mechanical properties. Altogether, by harnessing the fundamental varicose (jetting) and sinuous (buckling) instabilities associated with the extrusion of a submerged jelling filament, we show that it is possible to print complex three-dimensional filamentous structures in an ambient fluid.

Feedback control of protein aggregation

The self-assembly of peptides and proteins into amyloid fibrils plays a causative role in a wide range of increasingly common and currently incurable diseases. The molecular mechanisms underlying this process have recently been discovered, prompting the development of drugs that inhibit specific reaction steps as possible treatments for some of these disorders. A crucial part of treatment design is to determine how much drug to give and when to give it, informed by its efficacy and intrinsic toxicity. Since amyloid formation does not proceed at the same pace in different individuals, it is also important that treatment design is informed by local measurements of the extent of protein aggregation. Here, we use stochastic optimal control theory to determine treatment regimens for inhibitory drugs targeting several key reaction steps in protein aggregation, explicitly taking into account variability in the reaction kinetics. We demonstrate how these regimens may be updated "on the fly" as new measurements of the protein aggregate concentration become available, in principle, enabling treatments to be tailored to the individual. We find that treatment timing, duration, and drug dosage all depend strongly on the particular reaction step being targeted. Moreover, for some kinds of inhibitory drugs, the optimal regimen exhibits high sensitivity to stochastic fluctuations. Feedback controls tailored to the individual may therefore substantially increase the effectiveness of future treatments.

Wallpaper group kirigami

Kirigami, the art of paper cutting, has become a paradigm for mechanical metamaterials in recent years. The basic building blocks of any kirigami structures are repetitive deployable patterns that derive inspiration from geometric art forms and simple planar tilings. Here, we complement these approaches by directly linking kirigami patterns to the symmetry associated with the set of 17 repeating patterns that fully characterize the space of periodic tilings of the plane. We start by showing how to construct deployable kirigami patterns using any of the wallpaper groups, and then design symmetry-preserving cut patterns to achieve arbitrary size changes via deployment. We further prove that different symmetry changes can be achieved by controlling the shape and connectivity of the tiles and connect these results to the underlying kirigami-based lattice structures. All together, our work provides a systematic approach for creating a broad range of kirigami-based deployable structures with any prescribed size and symmetry properties.

Genetic architecture of floral traits in bee- and hummingbird-pollinated sister species of Aquilegia (columbine)

Interactions with animal pollinators have helped shape the stunning diversity of flower morphologies across the angiosperms. A common evolutionary consequence of these interactions is that some flowers have converged on suites of traits, or pollination syndromes, that attract and reward specific pollinator groups. Determining the genetic basis of these floral pollination syndromes can help us understand the processes that contributed to the diversification of the angiosperms. Here, we characterize the genetic architecture of a bee-to-hummingbird pollination shift in Aquilegia (columbine) using QTL mapping of 17 floral traits encompassing color, nectar composition, and organ morphology. In this system, we find that the genetic architectures underlying differences in floral color are quite complex, and we identify several likely candidate genes involved in anthocyanin and carotenoid floral pigmentation. Most morphological and nectar traits also have complex genetic underpinnings; however, one of the key floral morphological phenotypes, nectar spur curvature, is shaped by a single locus of large effect.

Puckering and wrinkling in a growing composite ring

Pattern formation driven by differential strain in constrained elastic systems is a common motif in many technological and biological systems. Here we introduce a biologically motivated case of elastic patterning that allows us to explore the conditions for the existence of local puckering and global wrinkling patterns: a soft growing composite ring adhered elastically to a constraining rigid ring. We explore how differential growth of the soft ring and the elastic resistance to shear and stretching deformations induced by soft adherence lead to a range of phenomena that include uniform aperture-like modes, localized puckers that are Nambu–Goldstone-like modes and global wrinkles in the system. Our analysis combines computer simulations of a discrete rod model with a nonlinear stability analysis of the differential equations in the continuum limit. We provide phase diagrams and scaling relations that reveal the nature and extent of the deformation patterns. Overall, our study reveals how geometry and mechanics conspire to yield a rich phenomenology that could serve as a guide to the design of programmable localized elastic deformations while being relevant for the mechanical basis of biological morphogenesis.

Static adhesion hysteresis in elastic structures

Adhesive interactions between elastic structures such as graphene sheets, carbon nanotubes, and microtubules have been shown to exhibit hysteresis due to irrecoverable energy loss associated with bond breakage, even in static (rate-independent) experiments. To understand this phenomenon, we start with a minimal theory for the peeling of a thin sheet from a substrate, coupling the local event of bond breaking to the nonlocal elastic relaxation of the sheet and show that this can drive static adhesion hysteresis over a bonding/debonding cycle. Using this model we quantify hysteresis in terms of the adhesion and elasticity parameters of the system. This allows us to derive a scaling relation that preserves hysteresis at different levels of granularity while resolving a seeming paradox of lattice trapping in the continuum limit of a discrete fracture process. Finally, to verify our theory, we use new experiments to demonstrate and measure adhesion hysteresis in bundled microtubules.

Models of benthic bipedalism

Walking is a common bipedal and quadrupedal gait and is often associated with terrestrial and aquatic organisms. Inspired by recent evidence of the neural underpinnings of primitive aquatic walking in the little skate Leucoraja erinacea, we introduce a theoretical model of aquatic walking that reveals robust and efficient gaits with modest requirements for body morphology and control. The model predicts undulatory behaviour of the system body with a regular foot placement pattern, which is also observed in the animal, and additionally predicts the existence of gait bistability between two states, one with a large energetic cost for locomotion and another associated with almost no energetic cost. We show that these can be discovered using a simple reinforcement learning scheme. To test these theoretical frameworks, we built a bipedal robot and show that its behaviours are similar to those of our minimal model: its gait is also periodic and exhibits bistability, with a low efficiency mode separated from a high efficiency mode by a 'jump' transition. Overall, our study highlights the physical constraints on the evolution of walking and provides a guide for the design of efficient biomimetic robots.

Mechanical basis for fibrillar bundle morphology

Understanding the morphology of self-assembled fibrillar bundles and aggregates is relevant to a range of problems in molecular biology, supramolecular chemistry and materials science. Here, we propose a coarse-grained approach that averages over specific molecular details and yields an effective mechanical theory for the spatial complexity of self-assembling fibrillar structures that arises due to the competing effects of (the bending and twisting) elasticity of individual filaments and the adhesive interactions between them. We show that our theoretical framework accounting for this allows us to capture a number of diverse fibril morphologies observed in natural and synthetic systems, ranging from Filopodia to multi-walled carbon nanotubes, and leads to a phase diagram of possible fibril shapes. We also show how the extreme sensitivity of these morphologies can lead to spatially chaotic structures. Together, these results suggest a common mechanical basis for mesoscale fibril morphology as a function of the nanoscale mechanical properties of its filamentous constituents.

Flow-Driven Branching in a Frangible Porous Medium

Channel formation and branching is widely seen in physical systems where movement of fluid through a porous structure causes the spatiotemporal evolution of the medium. We provide a simple theoretical framework that embodies this feedback mechanism in a multiphase model for flow through a frangible porous medium with a dynamic permeability. Numerical simulations of the model show the emergence of branched networks whose topology is determined by the geometry of external flow forcing. This allows us to delineate the conditions under which splitting and/or coalescing branched network formation is favored, with potential implications for both understanding and controlling branching in soft frangible media.

Coordinated crawling via reinforcement learning

Rectilinear crawling locomotion is a primitive and common mode of locomotion in slender soft-bodied animals. It requires coordinated contractions that propagate along a body that interacts frictionally with its environment. We propose a simple approach to understand how this coordination arises in a neuromechanical model of a segmented, soft-bodied crawler via an iterative process that might have both biological antecedents and technological relevance. Using a simple reinforcement learning algorithm, we show that an initial all-to-all neural coupling converges to a simple nearest-neighbour neural wiring that allows the crawler to move forward using a localized wave of contraction that is qualitatively similar to what is observed in Drosophila melanogaster larvae and used in many biomimetic solutions. The resulting solution is a function of how we weight gait regularization in the reward, with a trade-off between speed and robustness to proprioceptive noise. Overall, our results, which embed the brain-body-environment triad in a learning scheme, have relevance for soft robotics while shedding light on the evolution and development of locomotion.

Elastohydrodynamic Scaling Law for Heart Rates

Animal hearts are soft shells that actively pump blood to oxygenate tissues. Here, we propose an allometric scaling law for the heart rate based on the idea of elastohydrodynamic resonance of a fluid-loaded soft active elastic shell that buckles and contracts axially when twisted periodically. We show that this picture is consistent with numerical simulations of soft cylindrical shells that twist-buckle while pumping a viscous fluid, yielding optimum ejection fractions of 35%-40% when driven resonantly. Our scaling law is consistent with experimental measurements of heart rates over 2 orders of magnitude, and provides a mechanistic basis for how metabolism scales with organism size. In addition to providing a physical rationale for the heart rate and metabolism of an organism, our results suggest a simple design principle for soft fluidic pumps.

Suspension Jams in a Leaky Microfluidic Channel

Inspired by the jamming in leaky systems that arises in many physiological and industrial settings, we study the propagation of clogs in a leaky microfluidic channel. By driving a colloidal suspension through such a channel with a fluid-permeable wall adjoining a gutter, we follow the formation and propagation of jams and show that they move at a steady speed, in contrast with jams in channels that have impermeable walls. Furthermore, by varying the ratio of the resistance from the leaky wall and that of the gutter, we show that it is possible to control the shape of the propagating jam, which is typically wedge shaped. We complement our experiments with numerical simulations, where we implement an Euler-Lagrangian framework for the simultaneous evolution of both immersed colloidal particles and the carrier fluid. Finally, we show that the particle ordering in the clog can be tuned by adjusting the geometry of the leaky wall. Altogether, the leaky channel serves both as a filter and a shunt with the potential for a range of uses.

Self-Excited Motions of Volatile Drops on Swellable Sheets

When a volatile droplet is deposited on a floating swellable sheet, it becomes asymmetric, lobed and mobile. We describe and quantify this phenomena that involves nonequilibrium swelling, evaporation and motion, working together to realize a self-excitable spatially extended oscillator. Solvent penetration causes the film to swell locally and eventually buckle, changing its shape and the drop responds by moving. Simultaneously, solvent evaporation from the swollen film causes it to regain its shape once the droplet has moved away. The process repeats and leads to complex pulsatile spinning and/or sliding movements. We use a one-dimensional experiment to highlight the slow swelling of and evaporation from the film and the fast motion of the drop, a characteristic of excitable systems. Finally, we provide a phase diagram for droplet excitability as a function of drop size and film thickness and scaling laws for the motion of the droplet.

Genetic and Mechanical Regulation of Intestinal Smooth Muscle Development

The gastrointestinal tract is enveloped by concentric and orthogonally aligned layers of smooth muscle; however, an understanding of the mechanisms by which these muscles become patterned and aligned in the embryo has been lacking. We find that Hedgehog acts through Bmp to delineate the position of the circumferentially oriented inner muscle layer, whereas localized Bmp inhibition is critical for allowing formation of the later-forming, longitudinally oriented outer layer. Because the layers form at different developmental stages, the muscle cells are exposed to unique mechanical stimuli that direct their alignments. Differential growth within the early gut tube generates residual strains that orient the first layer circumferentially, and when formed, the spontaneous contractions of this layer align the second layer longitudinally. Our data link morphogen-based patterning to mechanically controlled smooth muscle cell alignment and provide a mechanistic context for potentially understanding smooth muscle organization in a wide variety of tubular organs.

Early warning signals in motion inference

The ability to infer intention lies at the basis of many social interactions played out via motor actions. We consider a simple paradigm of this ability in humans using data from experiments simulating an antagonistic game between an Attacker and a Blocker. Evidence shows early inference of an Attacker move by as much as 100ms but the nature of the informational cues signaling the impending move remains unknown. We show that the transition to action has the hallmark of a critical transition that is accompanied by early warning signals. These early warning signals occur as much as 130 ms before motion ensues-showing a sharp rise in motion autocorrelation at lag-1 and a sharp rise in the autocorrelation decay time. The early warning signals further correlate strongly with Blocker response times. We analyze the variance of the motion near the point of transition and find that it diverges in a manner consistent with the dynamics of a fold-transition. To test if humans can recognize and act upon these early warning signals, we simulate the dynamics of fold-transition events and ask people to recognize the onset of directional motion: participants react faster to fold-transition dynamics than to its uncorrelated counterpart. Together, our findings suggest that people can recognize the intent and onset of motion by inferring its early warning signals.

Spatio-temporal integration in plant tropisms

Tropisms, growth-driven responses to environmental stimuli, cause plant organs to respond in space and time and reorient themselves. Classical experiments from nearly a century ago reveal that plant shoots respond to the integrated history of light and gravity stimuli rather than just responding instantaneously. We introduce a temporally non-local response function for the dynamics of shoot growth formulated as an integro-differential equation whose solution allows us to qualitatively reproduce experimental observations associated with intermittent and unsteady stimuli. Furthermore, an analytic solution for the case of a pulse stimulus expresses the response function as a function of experimentally tractable variables, which we calculate for the case of the phototropic response of Arabidopsis hypocotyls. All together, our model enables us to predict tropic responses to time-varying stimuli, manifested in temporal integration phenomena, and sets the stage for the incorporation of additional effects such as multiple stimuli, gravitational sagging, etc.

Mechanics and kinetics of dynamic instability

During dynamic instability, self-assembling microtubules (MTs) stochastically alternate between phases of growth and shrinkage. This process is driven by the presence of two distinct states of MT subunits, GTP- and GDP-bound tubulin dimers, that have different structural properties. Here, we use a combination of analysis and computer simulations to study the mechanical and kinetic regulation of dynamic instability in three-dimensional (3D) self-assembling MTs. Our model quantifies how the 3D structure and kinetics of the distinct states of tubulin dimers determine the mechanical stability of MTs. We further show that dynamic instability is influenced by the presence of quenched disorder in the state of the tubulin subunit as reflected in the fraction of non-hydrolysed tubulin. Our results connect the 3D geometry, kinetics and statistical mechanics of these tubular assemblies within a single framework, and may be applicable to other self-assembled systems where these same processes are at play.

Rotation of a submerged finite cylinder moving down a soft incline

A submerged finite cylinder moving under its own weight along a soft incline lifts off and slides at a steady velocity while also spinning. Here, we experimentally quantify the steady spinning of the cylinder and show theoretically that it is due to a combination of an elastohydrodynamic torque generated by flow in the variable gap, and the viscous friction on the edges of the finite-length cylinder. The relative influence of the latter depends on the aspect ratio of the cylinder, the angle of the incline, and the deformability of the substrate, which we express in terms of a single scaled compliance parameter. By independently varying these quantities, we show that our experimental results are consistent with a transition from an edge-effect dominated regime for short cylinders to a gap-dominated elastohydrodynamic regime when the cylinder is very long.

Deterministic and stochastic control of kirigami topology

Kirigami, the creative art of paper cutting, is a promising paradigm for mechanical metamaterials. However, to make kirigami-inspired structures a reality requires controlling the topology of kirigami to achieve connectivity and rigidity. We address this question by deriving the maximum number of cuts (minimum number of links) that still allow us to preserve global rigidity and connectivity of the kirigami. A deterministic hierarchical construction method yields an efficient topological way to control both the number of connected pieces and the total degrees of freedom. A statistical approach to the control of rigidity and connectivity in kirigami with random cuts complements the deterministic pathway, and shows that both the number of connected pieces and the degrees of freedom show percolation transitions as a function of the density of cuts (links). Together, this provides a general framework for the control of rigidity and connectivity in planar kirigami.