Surface and interfacial creases in a bilayer tubular soft tissue

Mechanobiological tissue instability induced by stress-modulated growth

The growth of biological tissues, which is regulated by a variety of factors, can induce stresses that may, in turn, destabilize the tissues into diverse patterns. In most previous studies, however, tissue growth was usually assumed as a prescribed parameter independent of stresses, limiting our understanding of the mechanobiological morphogenesis of real tissues. In this paper, we propose a theoretical model to investigate the mechanobiological response of soft tissues undergoing stress-modulated growth. Linear stability analysis is first performed to elucidate the surface instability mechanism induced by stress-modulated volumetric growth. We further conduct finite element simulations to validate the theoretical prediction and, particularly, to capture the post-buckling pattern evolution. Our results show that the non-uniform stresses, which evolve with the tissue growth and morphogenesis, exert mechanical feedback on the growth itself, producing up-down asymmetric surface morphologies as observed in, for example, the gyrification of human brains and brain organoids. It is also revealed that large residual stresses are unnecessary to cause mechanobiological instability and subsequent asymmetric patterning, which has long been believed to be driven by sufficiently high stresses. The present work could help us to understand the morphological changes of biological tissues under physiological and pathological conditions.

Ballooning, bulging, and necking: An exact solution for longitudinal phase separation in elastic systems near a critical point

Prominent examples of longitudinal phase separation in elastic systems include elastic necking, the propagation of a bulge in a cylindrical party balloon, and the beading of a gel fiber subject to surface tension. Here we demonstrate that if the parameters of such a system are tuned near a critical point (where the difference between the two phases vanishes), then the behavior of all systems is given by the minimization of a simple and universal elastic energy familiar from Ginzburg-Landau theory in an external field. We minimize this energy analytically, which yields not only the well known interfacial tanh solution, but also the complete set of stable and unstable solutions in both finite and infinite length systems, unveiling the elastic system's full shape evolution and hysteresis. Correspondingly, we also find analytic results for the the delay of onset, changes in criticality, and ultimate suppression of instability with diminishing system length, demonstrating that our simple near-critical theory captures much of the complexity and choreography of far-from-critical systems. Finally, we find critical points for the three prominent examples of phase separation given above, and demonstrate how each system then follows the universal set of solutions.

Mechanics of pulmonary airways: Linking structure to function through constitutive modeling, biochemistry, and histology

Breathing involves fluid-solid interactions in the lung; however, the lack of experimental data inhibits combining the mechanics of air flow to airway deformation, challenging the understanding of how biomaterial constituents contribute to tissue response. As such, lung mechanics research is increasingly focused on exploring the relationship between structure and function. To address these needs, we characterize mechanical properties of porcine airways using uniaxial tensile experiments, accounting for bronchial orientation- and location- dependency. Structurally-reinforced constitutive models are developed to incorporate the role of collagen and elastin fibers embedded within the extrafibrillar matrix. The strain-energy function combines a matrix description (evaluating six models: compressible NeoHookean, unconstrained Ogden, uncoupled Mooney-Rivlin, incompressible Ogden, incompressible Demiray and incompressible NeoHookean), superimposed with non-linear fibers (evaluating two models: exponential and polynomial). The best constitutive formulation representative of all bronchial regions is determined based on curve-fit results to experimental data, accounting for uniqueness and sensitivity. Glycosaminoglycan and collagen composition, alongside tissue architecture, indicate fiber form to be primarily responsible for observed airway anisotropy and heterogeneous mechanical behavior. To the authors' best knowledge, this study is the first to formulate a structurally-motivated constitutive model, augmented with biochemical analysis and microstructural observations, to investigate the mechanical function of proximal and distal bronchi. Our systematic pulmonary tissue characterization provides a necessary foundation for understanding pulmonary mechanics; furthermore, these results enable clinical translation through simulations of airway obstruction in disease, fluid-structure interaction insights during breathing, and potentially, predictive capabilities for medical interventions. STATEMENT OF SIGNIFICANCE: The advancement of pulmonary research relies on investigating the biomechanical response of the bronchial tree. Experiments demonstrating the non-linear, heterogeneous, and anisotropic material behavior of porcine airways are used to develop a structural constitutive model representative of proximal and distal bronchial behavior. Calibrated material parameters exhibit regional variation in biomaterial properties, initially hypothesized to originate from tissue constituents. Further exploration through biochemical and histological analysis indicates mechanical function is primarily governed by microstructural form. The results of this study can be directly used in finite element and fluid-structure interaction models to enable physiologically relevant and more accurate computational simulations aimed to help diagnose and monitor pulmonary disease.

Pattern selection when a layer buckles on a soft substrate

If a neo-Hookean elastic layer adhered to a neo-Hookean substrate grows equibiaxially, it will buckle into a topographic pattern. Here, we combine higher order perturbation theory and finite element numerics to predict the pattern formed just beyond the buckling threshold. More precisely, we construct a series of solutions corresponding to hexagonal, square and stripe patterns, and expand the elastic energy for each pattern as a Landau-like energy series in the topography amplitude. We see that, for square and stripe patterns, the elastic energy is invariant under topography inversion, making the instabilities supercritical. However, since patterns of hexagonal dents are physically different to patterns of hexagonal bumps, the hexagonal energy lacks this invariance. This lack introduces a cubic term which causes hexagonal patterns to be formed subcritically and are hence energetically favoured. Our analytic calculation of the cubic term allows us to determine that dents are favoured in incompressible systems, but bumps are favored in sufficiently compressible systems. Finally, we consider a stiff layer sandwiched between an identical substrate and superstrate pair. This system has topography inversion symmetry, so hexagons form supercritically, and square patterns are favoured. We use finite element calculations to verify our theoretical predictions for each pattern, and confirm which pattern is selected. Previous work has used a simplified elastic model (a plate & a linear elastic substrate) that possesses invariance under topography inversion, and hence incorrectly predicted square patterns. Our work demonstrates that large strain geometry is sufficient to break this symmetry and explain the hexagonal dent patterns observed in buckling experiments.

Denser Growing Fiber Connections Induce 3-hinge Gyral Folding

Recent studies have shown that quantitative description of gyral shape patterns offers a novel window to examine the relationship between brain structure and function. Along this research line, this paper examines a unique and interesting type of cortical gyral region where 3 different gyral crests meet, termed 3-hinge gyral region. We extracted 3-hinge gyral regions in macaque/chimpanzee/human brains, quantified and compared the relevant DTI-derived fiber densities in 3-hinge and 2-hinge gyral regions. Our observations consistently showed that DTI-derived fiber densities in 3-hinge regions are much higher than those in 2-hinge regions. Therefore, we hypothesize that besides the cortical expansion, denser fiber connections can induce the formation of 3-hinge gyri. To examine the biomechanical basis of this hypothesis, we constructed a series of 3-dimensional finite element soft tissue models based on continuum growth theory to investigate fundamental biomechanical mechanisms of consistent 3-hinge gyri formation. Our computational simulation results consistently showed that during gyrification gyral regions with higher concentrations of growing axonal fibers tend to form 3-hinge gyri. Our integrative approach combining neuroimaging data analysis and computational modeling appears effective in probing a plausible theory of 3-hinge gyri formation and providing new insights into structural and functional cortical architectures and their relationship.

Peristaltic Elastic Instability in an Inflated Cylindrical Channel

A long cylindrical cavity through a soft solid forms a soft microfluidic channel, or models a vascular capillary. We observe experimentally that, when such a channel bears a pressurized fluid, it first dilates homogeneously, but then becomes unstable to a peristaltic elastic instability. We combine theory and numerics to fully characterize the instability in a channel with initial radius a through an incompressible bulk neo-Hookean solid with shear modulus μ. We show instability occurs supercritically with wavelength 12.278…a when the cavity pressure exceeds 2.052…μ. In finite solids, the wavelength for peristalsis lengthens, with peristalsis ultimately being replaced by a long-wavelength bulging instability in thin-walled cylinders. Peristalsis persists in Gent strain-stiffening materials, provided the material can sustain extension by more than a factor of 6. Although naively a pressure driven failure mode of soft channels, the instability also offers a route to fabricate periodically undulating channels, producing, e.g., waveguides with photonic or phononic stop bands.

Mechanical properties of the airway tree: heterogeneous and anisotropic pseudoelastic and viscoelastic tissue responses

Airway obstruction and pulmonary mechanics remain understudied despite lung disease being the third cause of death in the United States. Lack of relevant data has led computational pulmonary models to infer mechanical properties from available material data for the trachea. Additionally, the time-dependent, viscoelastic behaviors of airways have been largely overlooked, despite their potential physiological relevance and utility as metrics of tissue remodeling and disease progression. Here, we address the clear need for airway-specific material characterization to inform biophysical studies of the bronchial tree. Specimens from three airway levels (trachea, large bronchi, and small bronchi) and two orientations (axial and circumferential) were prepared from five fresh pig lungs. Uniaxial tensile tests revealed substantial heterogeneity and anisotropy. Overall, the linear pseudoelastic modulus was significantly higher axially than circumferentially (30.5 ± 3.1 vs. 8.4 ± 1.1 kPa) and significantly higher among circumferential samples for small bronchi than for the trachea and large bronchi (12.5 ± 1.9 vs. 6.0 ± 0.6 and 6.6 ± 0.9 kPa). Circumferential samples exhibited greater percent stress relaxation over 300 s than their axial counterparts (38.0 ± 1.4 vs. 23.1 ± 1.5%). Axial and circumferential trachea samples displayed greater percent stress relaxation (26.4 ± 1.6 and 42.5 ± 1.7%) than corresponding large and small bronchi. This ex vivo pseudoelastic and viscoelastic characterization reveals novel anisotropic and heterogeneous behaviors and equips us to construct airway-specific constitutive relations. Our results establish necessary fundamentals for airway mechanics, laying the groundwork for future studies to extend to clinical questions surrounding lung injury, and further directly enables computational tools for lung disease obstruction predictions. NEW & NOTEWORTHY Understanding the mechanics of the lung is necessary for investigating disease progression. Trachea mechanics comprises the vast majority of ex vivo airway tissue characterization despite distal airways being the site of disease manifestation and occlusion. Furthermore, viscoelastic studies are scarce, whereas time-dependent behaviors could be potential physiological metrics of tissue remodeling. In this study, the critical need for airway-specific material properties is addressed, reporting bronchial tree anisotropic and heterogeneous material properties.

Radial Structure Scaffolds Convolution Patterns of Developing Cerebral Cortex

Commonly-preserved radial convolution is a prominent characteristic of the mammalian cerebral cortex. Endeavors from multiple disciplines have been devoted for decades to explore the causes for this enigmatic structure. However, the underlying mechanisms that lead to consistent cortical convolution patterns still remain poorly understood. In this work, inspired by prior studies, we propose and evaluate a plausible theory that radial convolution during the early development of the brain is sculptured by radial structures consisting of radial glial cells (RGCs) and maturing axons. Specifically, the regionally heterogeneous development and distribution of RGCs controlled by Trnp1 regulate the convex and concave convolution patterns (gyri and sulci) in the radial direction, while the interplay of RGCs' effects on convolution and axons regulates the convex (gyral) convolution patterns. This theory is assessed by observations and measurements in literature from multiple disciplines such as neurobiology, genetics, biomechanics, etc., at multiple scales to date. Particularly, this theory is further validated by multimodal imaging data analysis and computational simulations in this study. We offer a versatile and descriptive study model that can provide reasonable explanations of observations, experiments, and simulations of the characteristic mammalian cortical folding.

Effect of Cyclic Loading on Surface Instability of Silicone Rubber under Compression

This work combines experiments and finite element simulations to study the effect of pre-imposed cyclic loading on surface instability of silicon rubber under compression. We first fabricate cuboid blocks of silicon rubber and pinch them cyclicly a few times. Then, an in-house apparatus is set to apply uniaxial compression on the silicon rubber under exact plane strain conditions. Surprisingly, we find multiple creases on the surface of silicone rubber, significantly different from what have been observed on the samples without the cyclic pinching. To reveal the underlying physics for these experimentally observed multiple creases, we perform detailed nanoindentation experiments to measure the material properties at different locations of the silicon rubber. The modulus is found to be nonuniform and varies along the thickness direction after the cyclic pinching. According to these experimental results, three-layer and multilayer finite element models are built with different materials properties informed by experiments. The three-layer finite element model can excellently explain the nucleation and pattern of multiple surface creases on the surface of compressed silicone rubber, in good agreement with experiments. Counterintuitively, the multilayer model with gradient modulus cannot be used to explain the multiple creases observed in our experiments. According to these simulations, the experimentally observed multiple creases should be attributed to a thin and stiff layer formed on the surface of silicon rubber after the pre-imposed cyclic loading.

Mechanism of Consistent Gyrus Formation: an Experimental and Computational Study

As a significant type of cerebral cortical convolution pattern, the gyrus is widely preserved across species. Although many hypotheses have been proposed to study the underlying mechanisms of gyrus formation, it is currently still far from clear which factors contribute to the regulation of consistent gyrus formation. In this paper, we employ a joint analysis scheme of experimental data and computational modeling to investigate the fundamental mechanism of gyrus formation. Experimental data on mature human brains and fetal brains show that thicker cortices are consistently found in gyral regions and gyral cortices have higher growth rates. We hypothesize that gyral convolution patterns might stem from heterogeneous regional growth in the cortex. Our computational simulations show that gyral convex patterns may occur in locations where the cortical plate grows faster than the cortex of the brain. Global differential growth can only produce a random gyrification pattern, but it cannot guarantee gyrus formation at certain locations. Based on extensive computational modeling and simulations, it is suggested that a special area in the cerebral cortex with a relatively faster growth speed could consistently engender gyri.