Pattern selection in growing tubular tissues

Hox gene activity directs physical forces to differentially shape chick small and large intestinal epithelia

Hox transcription factors play crucial roles in organizing developmental patterning across metazoa, but how these factors trigger regional morphogenesis has largely remained a mystery. In the developing gut, Hox genes help demarcate identities of intestinal subregions early in embryogenesis, which ultimately leads to their specialization in both form and function. Although the midgut forms villi, the hindgut develops sulci that resolve into heterogeneous outgrowths. Combining mechanical measurements of the embryonic chick intestine and mathematical modeling, we demonstrate that the posterior Hox gene HOXD13 regulates biophysical phenomena that shape the hindgut lumen. We further show that HOXD13 acts through the transforming growth factor β (TGF-β) pathway to thicken, stiffen, and promote isotropic growth of the subepithelial mesenchyme-together, these features lead to hindgut-specific surface buckling. TGF-β, in turn, promotes collagen deposition to affect mesenchymal geometry and growth. We thus identify a cascade of events downstream of positional identity that direct posterior intestinal morphogenesis.

Comparing finite viscoelastic constitutive relations and variational principles in modeling gastrointestinal soft tissue deformation

The mechanical attributes of soft tissues within the gastrointestinal (GI) tract are crucial for the effective operation of the GI system, and alterations in these properties may play a role in motility-related disorders. Various constitutive modeling approaches have been suggested to comprehend the response of soft tissues to diverse loading conditions. Among these, hyperelastic constitutive models based on finite elasticity have gained popularity. However, these models fall short in capturing rate- and time-dependent tissue properties. In contrast, finite viscoelastic models offer a solution to overcome these limitations. Nevertheless, the development of a suitable finite viscoelastic model, coupled with a variational formulation for efficient finite element (FE) implementation, remains an ongoing challenge. This study aims to address this gap by developing diverse finite viscoelastic constitutive relations and applying them to characterize soft tissue. Furthermore, the research explores the creation of compressible, nearly incompressible, and incompressible versions of viscoelastic constitutive relations, along with their variational formulation, to facilitate efficient FE implementation. The proposed model demonstrates remarkable accuracy in replicating experimental results, achieving an R2 value exceeding 0.99.

Elastic Instability behind Brittle Fracture

We argue that nucleation of brittle cracks in initially flawless soft elastic solids is preceded by a nonlinear elastic instability, which cannot be captured without accounting for geometrically precise description of finite elastic deformation. As a prototypical problem we consider a homogeneous elastic body subjected to tension and assume that it is weakened by the presence of a free surface which then serves as a location of cracks nucleation. We show that in this maximally simplified setting, brittle fracture emerges from a symmetry breaking elastic instability activated by softening and involving large elastic rotations. The implied bifurcation of the homogeneous elastic equilibrium is highly unconventional for nonlinear elasticity as it exhibits strong sensitivity to geometry, reminiscent of the transition to turbulence in fluids. We trace the postbifurcational development of this instability beyond the limits of applicability of scale-free continuum elasticity and use a phase-field approach to capture the scale dependent subcontinuum strain localization, signaling the formation of actual cracks.

Curvature-Regulated Multiphase Patterns in Tori

Biological functions in living systems are closely related to their geometries and morphologies. Toroidal structures, which widely exist in nature, present interesting features containing positive, zero, and negative Gaussian curvatures within one system. Such varying curvatures would significantly affect the growing or dehydrating morphogenesis, as observed in various intricate patterns in abundant biological structures. To understand the underlying morphoelastic mechanism and to determine the crucial factors that govern the patterning in toroidal structures, we develop a core-shell model and derive a scaling law to characterize growth- or dehydration-induced instability patterns. We find that the eventual patterns are mainly determined by two dimensionless parameters that are composed of stiffness and curvature of the system. Moreover, we construct a phase diagram showing the multiphase wrinkling pattern selection in various toroidal structures in terms of these two parameters, which is confirmed by our experimental observations. Physical insights into the multiphase transitions and existence of bistable modes are further provided by identifying hysteresis loops and the Maxwell equal-energy conditions. The universal law for morphology selection on core shell structures with varying curvatures can fundamentally explain and precisely predict wrinkling patterns of diverse toroidal structures, which may also provide a platform to design morphology-related functional surfaces.

Chiral topographic instability in shrinking spheres

Many biological structures exhibit intriguing morphological patterns adapted to environmental cues, which contribute to their important biological functions and also inspire material designs. Here, we report a chiral wrinkling topography in shrinking core-shell spheres, as observed in excessively dehydrated passion fruit and experimentally demonstrated in silicon core-shells under air extraction. Upon shrinkage deformation, the surface initially buckles into a buckyball pattern (periodic hexagons and pentagons) and then transforms into a chiral mode. The neighbouring chiral cellular patterns can further interact with each other, resulting in secondary symmetry breaking and the formation of two types of topological network. We develop a core-shell model and derive a universal scaling law to understand the underlying morphoelastic mechanism and to effectively describe and predict such chiral symmetry breaking far beyond the critical instability threshold. Moreover, we show experimentally that the chiral characteristic adapted to local perturbation can be harnessed to effectively and stably grasp small-sized objects of various shapes and made of different stiff and soft materials. Our results not only reveal chiral instability topographies, providing fundamental insights into the surface morphogenesis of the deformed core-shell spheres that are ubiquitous in the real world, but also demonstrate potential applications of adaptive grasping based on delicate chiral localization.

Pattern selection of directionally oriented chitosan tubes

The growth of viscoelastic curved materials, inspired by biological systems, may give rise to various complex structures. One of the simplest ways to control the pattern formation is to vary the orientation of the reaction vessel while keeping all other experimental conditions constant. Here, we report the self-organization of soft chitosan tubes by injecting acidic chitosan sol into a pool of sodium hydroxide solution, where the adhesive force between the gel and container keeps the tubules on the bottom of the reactor. The horizontal growth of the tubular structure undergoes spontaneous symmetry breaking, where instabilities develop on the surface of the chitosan tubules. Transformation of folds into wrinkles and finally to a smooth tube takes place by varying the orientation of the container. In addition to characterizing the evolving structures, we have also shown that the linear growth rate of the tube scales with the tilt angle of the container from the horizontal.

Mechanical regulation of tissues that reproduces wrinkle patterns of gastrointestinal tracts

Gastrointestinal tracts exhibit a number of surface morphologies including zigzags, labyrinths, protrusions, and invaginations which are associated with digestive functions and are suggested to be formed by mechanical mechanisms. In this study, we investigate loading conditions and mechanical properties of tissues that reproduce different wrinkle patterning of gastrointestinal tracts on cell culture platforms. Numerical simulations of wrinkling dynamics are performed for a layered model consisting of an anisotropic epithelial layer resting on a bimodular soft substrate, which in turn adheres to a rigid foundation. Motivated by the patterning of intestinal villi of chicks and mice, we examine two-step compression, where the epithelial layer is subjected to uniaxial compression followed by biaxial compression, and one-step compression, where the epithelial layer is compressed in biaxial directions. Under different mechanical conditions of tissues, a wide variety of surface patterns are displayed that reproduce luminal patterns of digestive tracts. These results suggest possible conditions for mechanical regulation of tissues to duplicate gastrointestinal surface patterns in vitro and provide insight into mechanistic understandings of biological tissues.

Influence of the growth gradient on surface wrinkling and pattern transition in growing tubular tissues

Growth-induced pattern formations in curved film-substrate structures have attracted extensive attention recently. In most existing literature, the growth tensor is assumed to be homogeneous or piecewise homogeneous. In this paper, we aim at clarifying the influence of a growth gradient on pattern formation and pattern evolution in bilayered tubular tissues under plane-strain deformation. In the framework of finite elasticity, a bifurcation condition is derived for a general material model and a generic growth function. Then we suppose that both layers are composed of neo-Hookean materials. In particular, the growth function is assumed to decay linearly either from the inner surface or from the outer surface. It is found that a gradient in the growth has a weak effect on the critical state, compared with the homogeneous growth type where both layers share the same growth factor. Furthermore, a finite-element model is built to validate the theoretical model and to investigate the post-buckling behaviours. It is found that the associated pattern transition is not controlled by the growth gradient but by the ratio of the shear modulus between two layers. Different morphologies can occur when the modulus ratio is varied. The current analysis could provide useful insight into the influence of a growth gradient on surface instabilities and suggests that a homogeneous growth field may provide a good approximation on interpreting complicated morphological formations in multiple systems.

Active elasticity drives the formation of periodic beading in damaged axons

In several pathological conditions, such as coronavirus infections, multiple sclerosis, Alzheimer's and Parkinson's diseases, the physiological shape of axons is altered and a periodic sequence of bulges appears. Experimental evidences suggest that such morphological changes are caused by the disruption of the microtubules composing the cytoskeleton of the axon. In this paper, we develop a mathematical model of damaged axons based on the theory of continuum mechanics and nonlinear elasticity. The axon is described as a cylinder composed of an inner passive part, called axoplasm, and an outer active cortex, composed mainly of F-actin and able to contract thanks to myosin-II motors. Through a linear stability analysis we show that, as the shear modulus of the axoplasm diminishes due to the disruption of the cytoskeleton, the active contraction of the cortex makes the cylindrical configuration unstable to axisymmetric perturbations, leading to a beading pattern. Finally, the nonlinear evolution of the bifurcated branches is investigated through finite element simulations.

Gas-Flow-Assisted Wrinkle-Free Transfer of a Centimeter-Scale Ultrathin Alumina Membrane onto Arbitrary Substrates

The transfer of an ultrathin membrane onto arbitrary substrates is important in different practical fields. Conventional wet-transfer methods inevitably induce wrinkle defects as a result of the large contact angle of the trapped droplet between the membrane and the substrate. Here, we demonstrate a gas flow-assisted method (GFAM) to transfer centimeter (cm)-scale ultrathin membranes onto arbitrary substrates (including a curved substrate) without wrinkles. GFAM makes use of contact angle hysteresis to bulge the trapped droplet between the substrate and the ultrathin membrane and simultaneously stretch the ultrathin membrane during rapid dewetting driven by gas flow. Moreover, GFAM can be easily fulfilled by using compressed air for seconds. Compared with conventional hydrophilic treatments or organic liquid wetting, this method has no durability concern and does not change the surface nature of substrates. Taking a widely used ultrathin anodic aluminum oxide (AAO) membrane as an example, we successfully demonstrate the application of a large-area wrinkle-free ultrathin AAO membrane to defect-free ordered nanostructure array fabrication and investigate the micro-scale details of macro-scale wrinkles generated by the conventional ways. In addition, its corresponding superiority over the defective counterpart is further studied in optical sensing. This method is highly valuable for promoting the simplicity of large-area ultrathin membrane transfer in practice.

Morphologies of compressed active epithelial monolayers

Using a three-dimensional active vertex model, we numerically study the shapes of strained unsupported epithelial monolayers subject to active junctional noise due to stochastic binding and unbinding of myosin. We find that while uniaxial, biaxial, and isotropic in-plane compressive strains do lead to the formation of longitudinal, herringbone pattern, and labyrinthine folds, respectively, the villus morphology characteristic of, e.g., the small intestine appears only if junctional tension fluctuations are strong enough to fluidize the tissue. Moreover, the fluidized epithelium features villi even in the absence of compressive strain provided that the apico-basal differential surface tension is large enough. We analyze several details of the different epithelial forms including the role of strain rate and the modulation of tissue thickness across folds. Our results show that even unsupported, non-patterned epithelia can form nontrivial morphologies.

Wrinkling instabilities for biologically relevant fiber-reinforced composite materials with a case study of Neo-Hookean/Ogden-Gasser-Holzapfel bilayer

Wrinkling is a ubiquitous surface phenomenon in many biological tissues and is believed to play an important role in arterial health. As arteries are highly nonlinear, anisotropic, multilayered composite systems, it is necessary to investigate wrinkling incorporating these material characteristics. Several studies have examined surface wrinkling mechanisms with nonlinear isotropic material relationships. Nevertheless, wrinkling associated with anisotropic constitutive models such as Ogden-Gasser-Holzapfel (OGH), which is suitable for soft biological tissues, and in particular arteries, still requires investigation. Here, the effects of OGH parameters such as fibers' orientation, stiffness, and dispersion on the onset of wrinkling, wrinkle wavelength and amplitude are elucidated through analysis of a bilayer system composed of a thin, stiff neo-Hookean membrane and a soft OGH substrate subjected to compression. Critical contractile strain at which wrinkles occur is predicted using both finite element analysis and analytical linear perturbation approach. Results suggest that besides stiffness mismatch, anisotropic features associated with fiber stiffness and distribution might be used in natural layered systems to adjust wrinkling and subsequent folding behaviors. Further analysis of a bilayer system with fibers in the (x-y) plane subjected to compression in the x direction shows a complex dependence of wrinkling strain and wavelength on fiber angle, stiffness, and dispersion. This behavior is captured by an approximation utilizing the linearized anisotropic properties derived from OGH model. Such understanding of wrinkling in this artery wall-like system will help identify the role of wrinkling mechanisms in biological artery in addition to the design of its synthetic counterparts.

Nonuniform growth and surface friction determine bacterial biofilm morphology on soft substrates

During development, organisms acquire three-dimensional (3D) shapes with important physiological consequences. While basic mechanisms underlying morphogenesis are known in eukaryotes, it is often difficult to manipulate them in vivo. To circumvent this issue, here we present a study of developing Vibrio cholerae biofilms grown on agar substrates in which the spatiotemporal morphological patterns were altered by varying the agar concentration. Expanding biofilms are initially flat but later undergo a mechanical instability and become wrinkled. To gain mechanistic insights into this dynamic pattern-formation process, we developed a model that considers diffusion of nutrients and their uptake by bacteria, bacterial growth/biofilm matrix production, mechanical deformation of both the biofilm and the substrate, and the friction between them. Our model shows quantitative agreement with experimental measurements of biofilm expansion dynamics, and it accurately predicts two distinct spatiotemporal patterns observed in the experiments-the wrinkles initially appear either in the peripheral region and propagate inward (soft substrate/low friction) or in the central region and propagate outward (stiff substrate/high friction). Our results, which establish that nonuniform growth and friction are fundamental determinants of stress anisotropy and hence biofilm morphology, are broadly applicable to bacterial biofilms with similar morphologies and also provide insight into how other bacterial biofilms form distinct wrinkle patterns. We discuss the implications of forming undulated biofilm morphologies, which may enhance the availability of nutrients and signaling molecules and serve as a "bet hedging" strategy.

Mechanics of human brain organoids

Organoids are prototypes of human organs derived from cultured human stem cells. They provide a reliable and accurate experimental model to study the physical mechanisms underlying the early developmental stages of human organs and, in particular, the early morphogenesis of the cortex. Here we propose a mathematical model to elucidate the role played by two mechanisms which have been experimentally proven to be crucial in shaping human brain organoids: the contraction of the inner core of the organoid and the microstructural remodeling of its outer cortex. Our results show that both mechanisms are crucial for the final shape of the organoid and that perturbing those mechanisms can lead to pathological morphologies which are reminiscent of those associated with lissencephaly (smooth brain).

Prescribing patterns in growing tubular soft matter by initial residual stress

Initial residual stress is omnipresent in biological tissues and soft matter, and can affect growth-induced pattern selection significantly. Here we demonstrate this effect experimentally by letting soft tubes grow in the presence or absence of initial residual stress and by observing different growth pattern evolutions. These experiments motivate us to model the mechanisms at play when a growing bilayer tubular organ spontaneously displays buckling patterns on its inner surface. We demonstrate that not only differential growth, geometry and elasticity, but also initial residual stress distribution, exert a notable influence on these pattern phenomena. Prescribing an initial residual stress distribution offers an alternative or a more effective way to implement pattern selection for growable bio-tissues or soft matter. The results also show promise for the design of 4D bio-mimic printing protocols or for controlling hydrogel actuators.

On Rate Boundary Conditions for Soft Tissue Bifurcation Analysis

Mechanical instability of soft tissues can either risk their normal function or alternatively trigger patterning mechanisms during growth and morphogenesis processes. Unlike standard stability analysis of linear elastic bodies, for soft tissues undergoing large deformations it is imperative to account for the nonlinearities induced by the coupling between load and surface changes at onset of instability. The related issue of boundary conditions, in context of soft tissues, has hardly been addressed in the literature, with most of available research employing dead-load conditions. This paper is concerned with the influence of imposed homogeneous rate (incremental) surface data on critical loads and associated modes in soft tissues, within the context of linear bifurcation analysis. Material behavior is modeled by compressible isotropic hyperelastic strain energy functions (SEFs), with experimentally validated material parameters for the Fung-Demiray SEF, over a range of constitutive response (including brain and liver tissues). For simplicity, we examine benchmark problems of basic spherical patterns: full sphere, spherical cavity, and thick spherical shell. Limiting the analysis to primary hydrostatic states we arrive at universal closed-form solutions, thus providing insight on the role of imposed boundary data. Influence of selected rate boundary conditions (RBCs) like dead-load and fluid-pressure (FP), coupled with constitutive parameters, on the existence and levels of bifurcation loads is compared and discussed. It is argued that the selection of the appropriate type of homogeneous RBC can have a critical effect on the level of bifurcation loads and even exclude the emergence of bifurcation instabilities.

Growth and remodelling of living tissues: perspectives, challenges and opportunities

One of the most remarkable differences between classical engineering materials and living matter is the ability of the latter to grow and remodel in response to diverse stimuli. The mechanical behaviour of living matter is governed not only by an elastic or viscoelastic response to loading on short time scales up to several minutes, but also by often crucial growth and remodelling responses on time scales from hours to months. Phenomena of growth and remodelling play important roles, for example during morphogenesis in early life as well as in homeostasis and pathogenesis in adult tissues, which often adapt to changes in their chemo-mechanical environment as a result of ageing, diseases, injury or surgical intervention. Mechano-regulated growth and remodelling are observed in various soft tissues, ranging from tendons and arteries to the eye and brain, but also in bone, lower organisms and plants. Understanding and predicting growth and remodelling of living systems is one of the most important challenges in biomechanics and mechanobiology. This article reviews the current state of growth and remodelling as it applies primarily to soft tissues, and provides a perspective on critical challenges and future directions.

Non-uniform curvature and anisotropic deformation control wrinkling patterns on tori

We investigate wrinkling patterns in a tri-layer torus consisting of an expanding thin outer layer, an intermediate soft layer and an inner core with a tunable shear modulus, inspired by pattern formation in developmental biology, such as follicle pattern formation during the development of chicken embryos. We show from large-scale finite element simulations that hexagonal wrinkling patterns form for stiff cores whereas stripe wrinkling patterns develop for soft cores. Hexagons and stripes co-exist to form hybrid patterns for cores with intermediate stiffness. The governing mechanism for the pattern transition is that the stiffness of the inner core controls the degree to which the major radius of the torus expands - this has a greater effect on deformation in the long direction as compared to the short direction of the torus. This anisotropic deformation alters stress states in the outer layer which change from biaxial (preferred hexagons) to uniaxial (preferred stripes) compression as the core stiffness is reduced. As the outer layer continues to expand, stripe and hexagon patterns will evolve into zigzags and segmented labyrinths, respectively. Stripe wrinkles are observed to initiate at the inner surface of the torus while hexagon wrinkles start from the outer surface as a result of curvature-dependent stresses in the torus. We further discuss the effects of elasticities and geometries of the torus on the wrinkling patterns.

Influence of Initial Residual Stress on Growth and Pattern Creation for a Layered Aorta

Residual stress is ubiquitous and indispensable in most biological and artificial materials, where it sustains and optimizes many biological and functional mechanisms. The theory of volume growth, starting from a stress-free initial state, is widely used to explain the creation and evolution of growth-induced residual stress and the resulting changes in shape, and to model how growing bio-tissues such as arteries and solid tumors develop a strategy of pattern creation according to geometrical and material parameters. This modelling provides promising avenues for designing and directing some appropriate morphology of a given tissue or organ and achieve some targeted biomedical function. In this paper, we rely on a modified, augmented theory to reveal how we can obtain growth-induced residual stress and pattern evolution of a layered artery by starting from an existing, non-zero initial residual stress state. We use experimentally determined residual stress distributions of aged bi-layered human aortas and quantify their influence by a magnitude factor. Our results show that initial residual stress has a more significant impact on residual stress accumulation and the subsequent evolution of patterns than geometry and material parameters. Additionally, we provide an essential explanation for growth-induced patterns driven by differential growth coupled to an initial residual stress. Finally, we show that initial residual stress is a readily available way to control growth-induced pattern creation for tissues and thus may provide a promising inspiration for biomedical engineering.

On the existence of elastic minimizers for initially stressed materials

A soft solid is said to be initially stressed if it is subjected to a state of internal stress in its unloaded reference configuration. In physical terms, its stored elastic energy may not vanish in the absence of an elastic deformation, being also dependent on the spatial distribution of the underlying material inhomogeneities. Developing a sound mathematical framework to model initially stressed solids in nonlinear elasticity is key for many applications in engineering and biology. This work investigates the links between the existence of elastic minimizers and the constitutive restrictions for initially stressed materials subjected to finite deformations. In particular, we consider a subclass of constitutive responses in which the strain energy density is taken as a scalar-valued function of both the deformation gradient and the initial stress tensor. The main advantage of this approach is that the initial stress tensor belongs to the group of divergence-free symmetric tensors satisfying the boundary conditions in any given reference configuration. However, it is still unclear which physical restrictions must be imposed for the well-posedness of this elastic problem. Assuming that the constitutive response depends on the choice of the reference configuration only through the initial stress tensor, under given conditions we prove the local existence of a relaxed state given by an implicit tensor function of the initial stress distribution. This tensor function is generally not unique, and can be transformed according to the symmetry group of the material at fixed initial stresses. These results allow one to extend Ball's existence theorem of elastic minimizers for the proposed constitutive choice of initially stressed materials. This article is part of the theme issue 'Rivlin's legacy in continuum mechanics and applied mathematics'.

Chick midgut morphogenesis

The gastrointestinal tract is an essential system of organs required for nutrient absorption. As a simple tube early in development, the primitive gut is patterned along its anterior-posterior axis into discrete compartments with unique morphologies relevant to their functions in the digestive process. These morphologies are acquired gradually through development as the gut is patterned by tissue interactions, both molecular and mechanical in nature, involving all three germ layers. With a focus on midgut morphogenesis, we review work in the chick embryo demonstrating how these molecular signals and mechanical forces sculpt the developing gut tube into its mature form. In particular, we highlight two mechanisms by which the midgut increases its absorptive surface area: looping and villification. Additionally, we review the differentiation and patterning of the intestinal mesoderm into the layers of smooth muscle that mechanically drive peristalsis and the villification process itself. Where relevant, we discuss the mechanisms of chick midgut morphogenesis in the context of experimental data from other model systems.

Understanding the mechanical link between oriented cell division and cerebellar morphogenesis

The cerebellum is a tightly folded structure located at the back of the head. Unlike the folds of the cerebrum, the folds of the cerebellum are aligned such that the external surface appears to be covered in parallel grooves. Experiments have shown that anchoring center initiation drives cerebellar foliation. However, the mechanism guiding the location of these anchoring centers, and subsequently cerebellar morphology, remains poorly understood. In particular, there is no definitive mechanistic explanation for the preferential emergence of parallel folds instead of an irregular folding pattern like in the cerebral cortex. Here we use mechanical modeling on the cellular and tissue scales to show that the oriented granule cell division observed in the experimental setting leads to the characteristic parallel folding pattern of the cerebellum. Specifically, we propose an agent-based model of cell clones, a strategy for propagating information from our in silico cell clones to the tissue scale, and an analytical solution backed by numerical results to understand how differential growth between the cerebellar layers drives geometric instability in three dimensional space on the tissue scale. This proposed mechanical model provides further insight into the process of anchoring center initiation and establishes a framework for future multiscale mechanical analysis of developing organs.

Photocontrollable Wrinkle Morphology Evolution on Azo-Based Multilayers for Hierarchical Surface Micropatterns Fabrication

Inspired by nature, comprehensive understanding and ingenious utilization of the self-organized wrinkling behaviors of the sandwiched multilayer bonded on substrates are important for engineering and/or functional laminated devices design. Herein, we report a facile and effective strategy to regulate the wrinkles morphology evolution and the resultant hierarchical surface micropatterns on azobenzene-based laminated multilayers by visible-light irradiation. Revealed by systematic experiments, the photocontrolled dynamic wrinkle evolutions are triggered by the reversible photoisomerization of azobenzene in the top azopolymer film and are strongly dependent on the intermediate photoinert layers (e.g., polystyrene and oxygen plasma-induced SiO _x layer) with the wrinkle-reinforcing effect or the stress relaxation acceleration effect. Interestingly, large-area well-defined hierarchical surface wrinkle patterns could be fabricated on the multilayers upon selective exposure. In the unexposed region, the wrinkles evolved into highly oriented patterns, whereas in the exposed region, they were fully erased or evolved into smaller-wavelength wrinkles. This study not only sheds light on the morphological evolution of the wrinkling laminated composites in engineering and nature but also paves a new avenue to conveniently and controllably realize the hierarchical stimulus-responsive surface patterns.

Pattern Transitions in a Soft Cylindrical Shell

Instability patterns of rolling up a sleeve appear more intricate than the ones of walking over a rug on floor, both characterized as systems of uniaxially compressed soft film on stiff substrate. This can be explained by curvature effects. To investigate pattern transitions on a curved surface, we study a soft shell sliding on a rigid cylinder by experiments, computations and theoretical analyses. We reveal a novel postbuckling phenomenon involving multiple successive bifurcations: smooth-wrinkle-ridge-sagging transitions. The shell initially buckles into periodic axisymmetric wrinkles at the threshold and then a wrinkle-to-ridge transition occurs upon further axial compression. When the load increases to the third bifurcation, the amplitude of the ridge reaches its limit and the symmetry is broken with the ridge sagging into a recumbent fold. It is identified that hysteresis loops and the Maxwell equal-energy conditions are associated with the coexistence of wrinkle-ridge or ridge-sagging patterns. Such a bifurcation scenario is inherently general and independent of material constitutive models.

Mechanical morphogenesis and the development of neocortical organisation

The development of complex neocortical organisations is thought to result from the interaction of genetic and activity-dependent processes. We propose that a third type of process - mechanical morphogenesis - may also play an important role. We review theoretical and experimental results in physics showing how even homogeneous growth can produce a variety of forms, in particular neocortical folding. The mechanical instabilities that produce these forms induce heterogeneous patterns of stress at the scale of the organ. We review the evidence showing how these stresses can influence cell proliferation, migration and apoptosis, cell differentiation and shape, migration and axonal guidance, and could thus be able to influence regional neocortical identity and connectivity.

Wrinkling patterns in soft shells

Curvature plays an important role in the morphological evolution of soft shells under stretch. Here, through a combination of experiment, theory and simulation, we investigate the behavior of a hemispherical soft shell subject to an increasing outward point force at its pole. In contrast to an inward point force inducing a polygonal pattern of buckling in the shell, we observe a four-stage morphological transition and symmetry breaking under an increasing outward point force. The shell undergoes axisymmetric deformation around its pole and then buckles into a non-axisymmetric shape with a number of shallow wrinkles emanating from the pole, followed by the emergence of crater-like deep crumples and ultimately a transformation into a wrinkled pseudocone. Our theoretical analysis and numerical simulations yield the critical conditions for the morphological transitions at each stage of deformation and reveal the underlying interplays between elastic bending and stretching energies and the curvature of the shell.

Matched asymptotic solution for crease nucleation in soft solids

A soft solid subjected to a large compression develops sharp self-contacting folds at its free surface, known as creases. Creasing is physically different from structural elastic instabilities, like buckling or wrinkling. Indeed, it is a fully nonlinear material instability, similar to a phase-transformation. This work provides theoretical insights of the physics behind crease nucleation. Creasing is proved to occur after a global bifurcation allowing the co-existence of an outer deformation and an inner solution with localised self-contact at the free surface. The most fundamental result here is the analytic prediction of the nucleation threshold, in excellent agreement with experiments and numerical simulations. A matched asymptotic solution is given within the intermediate region between the two co-existing states. The self-contact acts like the point-wise disturbance in the Oseen's correction for the Stokes flow past a circle. Analytic expressions of the matching solution and its range of validity are also derived.

Quantitative predictions of diverse wrinkling patterns in film/substrate systems

A basic characteristic of stiff film/soft substrate systems is their ability to experience large deformation under compressive stresses, which inevitably leads to formation of patterns on the surface. Such pattern formation is the result of loss of stability and symmetry breaking. Knowledge on how such instabilities arise and evolve is essential to describe, understand, predict, and ultimately to design complex functional materials and structures, for example the fabrication of stretchable electronic devices and micro/nano-scale surface patterning control. In this paper, quantitative predictions of various instability pattern formations and evolutions, which involve highly nonlinear deformation and multiple bifurcations, will be presented based on advanced mechanical models and methods, from planar to curved geometry. The results can provide further insight into fundamental understanding in a whole view of a variety of surface patterning morphology and imply a potential way to facilitate the design of functional materials and structures by quantitatively harnessing surface instabilities.

Rayleigh-Taylor instability in soft elastic layers

This work investigates the morphological stability of a soft body composed of two heavy elastic layers attached to a rigid surface and subjected only to the bulk gravity force. Using theoretical and computational tools, we characterize the selection of different patterns as well as their nonlinear evolution, unveiling the interplay between elastic and geometric effects for their formation. Unlike similar gravity-induced shape transitions in fluids, such as the Rayleigh-Taylor instability, we prove that the nonlinear elastic effects saturate the dynamic instability of the bifurcated solutions, displaying a rich morphological diagram where both digitations and stable wrinkling can emerge. The results of this work provide important guidelines for the design of novel soft systems with tunable shapes, with several applications in engineering sciences.This article is part of the themed issue 'Patterning through instabilities in complex media: theory and applications.'

Wrinkling instabilities in soft bilayered systems

Wrinkling phenomena control the surface morphology of many technical and biological systems. While primary wrinkling has been extensively studied, experimentally, analytically and computationally, higher-order instabilities remain insufficiently understood, especially in systems with stiffness contrasts well below 100. Here, we use the model system of an elastomeric bilayer to experimentally characterize primary and secondary wrinkling at moderate stiffness contrasts. We systematically vary the film thickness and substrate prestretch to explore which parameters modulate the emergence of secondary instabilities, including period-doubling, period-tripling and wrinkle-to-fold transitions. Our experiments suggest that period-doubling is the favourable secondary instability mode and that period-tripling can emerge under disturbed boundary conditions. High substrate prestretch can suppress period-doubling and primary wrinkles immediately transform into folds. We combine analytical models with computational simulations to predict the onset of primary wrinkling, the post-buckling behaviour, secondary bifurcations and the wrinkle-to-fold transition. Understanding the mechanisms of pattern selection and identifying the critical control parameters of wrinkling will allow us to fabricate smart surfaces with tunable properties and to control undesired surface patterns like in the asthmatic airway.This article is part of the themed issue 'Patterning through instabilities in complex media: theory and applications.'

Mechanics Reveals the Biological Trigger in Wrinkly Fingers

Fingertips wrinkle due to long exposure to water. The biological reason for this morphological change is unclear and still not fully understood. There are two main hypotheses for the underlying mechanism of fingertip wrinkling: the 'shrink' model (in which the wrinkling is driven by the contraction of the lower layers of skin, associated with the shrinking of the underlying vasculature), and the 'swell' model (in which the wrinkling is driven by the swelling of the upper layers of the skin, associated with osmosis). In reality, contraction of the lower layers of the skin and swelling of the upper layers will happen simultaneously. However, the relative importance of these two mechanisms to drive fingertip wrinkling also remains unclear. Simulating the swelling in the upper layers of skin alone, which is associated with neurological disorders, we found that wrinkles appeared above an increase of volume of [Formula: see text] Therefore, the upper layers can not exceed this swelling level in order to not contradict in vivo observations in patients with such neurological disorders. Simulating the contraction of the lower layers of the skin alone, we found that the volume have to decrease a [Formula: see text] to observe wrinkles. Furthermore, we found that the combined effect of both mechanisms leads to pronounced wrinkles even at low levels of swelling and contraction when individually they do not. This latter results indicates that the collaborative effect of both hypothesis are needed to induce wrinkles in the fingertips. Our results demonstrate how models from continuum mechanics can be successfully applied to testing hypotheses for the mechanisms that underly fingertip wrinkling, and how these effects can be quantified.

Inhomogeneous distribution of Alzheimer pathology along the isocortical relief. Are cortical convolutions an Achilles heel of evolution?

Alzheimer's disease (AD) is neuropathologically characterized by neuritic plaques and neurofibrillary tangles. Progression of both plaques and tangles throughout the brain follows a hierarchical distribution which is defined by intrinsic cytoarchitectonic features and extrinsic connectivity patterns. What has less well been studied is how cortical convolutions influence the distribution of AD pathology. Here, the distribution of both plaques and tangles within subsulcal gyral components (fundi) to components forming their top regions at the subarachnoidal brain surface (crowns) by stereological methods in seven different cortical areas was systematically compared. Further, principle differences in cytoarchitectonic organization of cortical crowns and fundi that might provide the background for regionally selective vulnerability were attempted to identify. It was shown that both plaques and tangles were more prominent in sulcal fundi than gyri crowns. The differential distribution of pathology along convolutions corresponds to subgyral differences in the vascular network, GFAP-positive astrocytes and intracortical and subcortical connectivity. While the precise mechanisms accounting for these differences remain open, the presence of systematic inhomogeneities in the distribution of AD pathology along cortical convolutions indicates that the phylogenetic shaping of the cortex is associated with features that render the human brain vulnerable to AD pathology.

Tri-layer wrinkling as a mechanism for anchoring center initiation in the developing cerebellum

During cerebellar development, anchoring centers form at the base of each fissure and remain fixed in place while the rest of the cerebellum grows outward. Cerebellar foliation has been extensively studied; yet, the mechanisms that control anchoring center initiation and position remain insufficiently understood. Here we show that a tri-layer model can predict surface wrinkling as a potential mechanism to explain anchoring center initiation and position. Motivated by the cerebellar microstructure, we model the developing cerebellum as a tri-layer system with an external molecular layer and an internal granular layer of similar stiffness and a significantly softer intermediate Purkinje cell layer. Including a weak intermediate layer proves key to predicting surface morphogenesis, even at low stiffness contrasts between the top and bottom layers. The proposed tri-layer model provides insight into the hierarchical formation of anchoring centers and establishes an essential missing link between gene expression and evolution of shape.

Elastosis during airway wall remodeling explains multiple co-existing instability patterns

Living structures can undergo morphological changes in response to growth and alterations in microstructural properties in response to remodeling. From a biological perspective, airway wall inflammation and airway elastosis are classical hallmarks of growth and remodeling during chronic lung disease. From a mechanical point of view, growth and remodeling trigger mechanical instabilities that result in inward folding and airway obstruction. While previous analytical and computational studies have focused on identifying the critical parameters at the onset of folding, few have considered the post-buckling behavior. All prior studies assume constant microstructural properties during the folding process; yet, clinical studies now reveal progressive airway elastosis, the degeneration of elastic fibers associated with a gradual stiffening of the inner layer. Here, we explore the influence of temporally evolving material properties on the post-bifurcation behavior of the airway wall. We show that a growing and stiffening inner layer triggers an additional subsequent bifurcation after the first instability occurs. Evolving material stiffnesses provoke failure modes with multiple co-existing wavelengths, associated with the superposition of larger folds evolving on top of the initial smaller folds. This phenomenon is exclusive to material stiffening and conceptually different from the phenomenon of period doubling observed in constant-stiffness growth. Our study suggests that the clinically observed multiple wavelengths in diseased airways are a result of gradual airway wall stiffening. While our evolving material properties are inspired by the clinical phenomenon of airway elastosis, the underlying concept is broadly applicable to other types of remodeling including aneurysm formation or brain folding.

Quantitative Morphology of Epithelial Folds

The shape of spatially modulated epithelial morphologies such as villi and crypts is usually associated with the epithelium-stroma area mismatch leading to buckling. We propose an alternative mechanical model based on intraepithelial stresses generated by differential tensions of apical, lateral, and basal sides of cells as well as on the elasticity of the basement membrane. We use it to theoretically study longitudinal folds in simple epithelia and we identify four types of corrugated morphologies: compact, invaginated, evaginated, and wavy. The obtained tissue contours and thickness profiles are compared to epithelial folds observed in invertebrates and vertebrates, and for most samples, the agreement is within the estimated experimental error. Our model establishes the groove-crest modulation of tissue thickness as a morphometric parameter that can, together with the curvature profile, be used to estimate the relative differential apicobasal tension in the epithelium.

Modulation of apical constriction by Wnt signaling is required for lung epithelial shape transition

In lung development the apically constricted columnar epithelium forms numerous buds during the pseudoglandular stage and subsequently changes the shape into flat or cuboidal pneumocytes that compose the air sacs during the canalicular and saccular (canalicular-saccular) stages, yet the impact of cell shapes on tissue morphogenesis remains unclear. The expression of Wnt components were decreased in the canalicular-saccular stages, and genetically constitutive activation of Wnt signaling impaired air sac formation by inducing apical constriction in the epithelium as seen in the pseudoglandular stage. Organ culture models also demonstrated that Wnt signaling induces apical constriction through the apical actomyosin cytoskeletal organization. Mathematical modeling revealed that apical constriction induces bud formation and loss of apical constriction is required for the formation of an air sac-like structure. MAP/Microtubule affinity-regulating kinase (MARK1) was identified as a downstream molecule of Wnt signaling and required for the apical cytoskeletal organization and bud formation. These results suggest that Wnt signaling is required for bud formation by inducing apical constriction during the pseudoglandular stage, while loss of Wnt signaling is for air sac formation in the canalicular-saccular stages.

Systems biology and mechanics of growth

In contrast to inert systems, living biological systems have the advantage to adapt to their environment through growth and evolution. This transfiguration is evident during embryonic development, when the predisposed need to grow allows form to follow function. Alterations in the equilibrium state of biological systems breed disease and mutation in response to environmental triggers. The need to characterize the growth of biological systems to better understand these phenomena has motivated the continuum theory of growth and stimulated the development of computational tools in systems biology. Biological growth in development and disease is increasingly studied using the framework of morphoelasticity. Here, we demonstrate the potential for morphoelastic simulations through examples of volume, area, and length growth, inspired by tumor expansion, chronic bronchitis, brain development, intestine formation, plant shape, and myopia. We review the systems biology of living systems in light of biochemical and optical stimuli and classify different types of growth to facilitate the design of growth models for various biological systems within this generic framework. Exploring the systems biology of growth introduces a new venue to control and manipulate embryonic development, disease progression, and clinical intervention.

The effect of mechanical-driven volumetric change on instability patterns of bilayered soft solids

If a soft solid is compressible, its volume changes with imposed loading. The extent of the volume change depends on its Poisson's ratio. Here, we study the effect of the mechanical-driven volumetric change on buckling and post-buckling behaviors of a hard thin film perfectly bound on a compliant substrate through the theoretical analysis and finite element method. Poisson's ratio of the substrate has been chosen to be in the range of -1 to 0.5, allowing its volume change during deformation. We find that Poisson's ratio cannot only shift the critical strain for the onset of buckling, but also affect the buckling modes. When Poisson's ratio of the substrate is close to -1, the surface instabilities of the thin film can be suppressed and delayed to large deformation. The present study demonstrates a new way to control surface instabilities of a bilayered system by changing Poisson's ratio of the material.

Hierarchical wrinkling in a confined permeable biogel

Confined thin surfaces may wrinkle as a result of the growth of excess material. Elasticity or gravity usually sets the wavelength. We explore new selection mechanisms based on hydrodynamics. First, inspired by yoghurt-making processes, we use caseins (a family of milk proteins) as pH-responsive building blocks and the acidulent glucono-δ-lactone to design a porous biogel film immersed in a confined buoyancy-matched viscous medium. Under specific boundary conditions yet without any external stimulus, the biogel film spontaneously wrinkles in cascade. Second, using a combination of titration, rheology, light microscopy, and confocal microscopy, we demonstrate that, during continuous acidification, the gel first shrinks and then swells, inducing wrinkling. Third, taking into account both Darcy flow through the gel and Poiseuille flow in the surrounding solvent, we develop a model that correctly predicts the wrinkling wavelength. Our results should be universal for acid-induced protein gels because they are based on pH-induced charge stabilization/destabilization and therefore could set a benchmark to gain fundamental insights into wrinkled biological tissues, to texture food, or to design surfaces for optical purposes.

Physical biology of human brain development

Neurodevelopment is a complex, dynamic process that involves a precisely orchestrated sequence of genetic, environmental, biochemical, and physical events. Developmental biology and genetics have shaped our understanding of the molecular and cellular mechanisms during neurodevelopment. Recent studies suggest that physical forces play a central role in translating these cellular mechanisms into the complex surface morphology of the human brain. However, the precise impact of neuronal differentiation, migration, and connection on the physical forces during cortical folding remains unknown. Here we review the cellular mechanisms of neurodevelopment with a view toward surface morphogenesis, pattern selection, and evolution of shape. We revisit cortical folding as the instability problem of constrained differential growth in a multi-layered system. To identify the contributing factors of differential growth, we map out the timeline of neurodevelopment in humans and highlight the cellular events associated with extreme radial and tangential expansion. We demonstrate how computational modeling of differential growth can bridge the scales-from phenomena on the cellular level toward form and function on the organ level-to make quantitative, personalized predictions. Physics-based models can quantify cortical stresses, identify critical folding conditions, rationalize pattern selection, and predict gyral wavelengths and gyrification indices. We illustrate that physical forces can explain cortical malformations as emergent properties of developmental disorders. Combining biology and physics holds promise to advance our understanding of human brain development and enable early diagnostics of cortical malformations with the ultimate goal to improve treatment of neurodevelopmental disorders including epilepsy, autism spectrum disorders, and schizophrenia.

Patient-Specific Airway Wall Remodeling in Chronic Lung Disease

Chronic lung disease affects more than a quarter of the adult population; yet, the mechanics of the airways are poorly understood. The pathophysiology of chronic lung disease is commonly characterized by mucosal growth and smooth muscle contraction of the airways, which initiate an inward folding of the mucosal layer and progressive airflow obstruction. Since the degree of obstruction is closely correlated with the number of folds, mucosal folding has been extensively studied in idealized circular cross sections. However, airflow obstruction has never been studied in real airway geometries; the behavior of imperfect, non-cylindrical, continuously branching airways remains unknown. Here we model the effects of chronic lung disease using the nonlinear field theories of mechanics supplemented by the theory of finite growth. We perform finite element analysis of patient-specific Y-branch segments created from magnetic resonance images. We demonstrate that the mucosal folding pattern is insensitive to the specific airway geometry, but that it critically depends on the mucosal and submucosal stiffness, thickness, and loading mechanism. Our results suggests that patient-specific airway models with inherent geometric imperfections are more sensitive to obstruction than idealized circular models. Our models help to explain the pathophysiology of airway obstruction in chronic lung disease and hold promise to improve the diagnostics and treatment of asthma, bronchitis, chronic obstructive pulmonary disease, and respiratory failure.

Period-doubling and period-tripling in growing bilayered systems

Growing layers on elastic substrates are capable of creating a wide variety of surface morphologies. Moderate growth generates a regular pattern of sinusoidal wrinkles with a homogeneous energy distribution. While the critical conditions for periodic wrinkling have been extensively studied, the rich pattern formation beyond this first instability point remains poorly understood. Here we show that upon continuing growth, the energy progressively localizes and new complex morphologies emerge. Previous studies have often overlooked these secondary bifurcations; they have focused on large stiffness ratios between layer and substrate, where primary instabilities occur early, long before secondary instabilities emerge. We demonstrate that secondary bifurcations are particularly critical in the low stiffness ratio regime, where the critical conditions for primary and secondary instabilities move closer together. Amongst all possible secondary bifurcations, the mode of period-doubling plays a central role - it is energetically favorable over all other modes. Yet, we can numerically suppress period-doubling, by choosing boundary conditions, which favor alternative higher order modes. Our results suggest that in the low stiffness regime, pattern formation is highly sensitive to small imperfections: surface morphologies emerge rapidly, change spontaneously, and quickly become immensely complex. This is a common paradigm in developmental biology. Our results have significantly applications in the morphogenesis of living systems where growth is progressive and stiffness ratios are low.

Morphological patterns of a growing biological tube in a confined environment with contacting boundary

Morphological transition from creasing to detachment of a growing biological tube in a confined environment with the contacting boundary.