Interplay between activity, elasticity, and liquid transport in self-contractile biopolymer gels.

Active gels play an important role in biology and in inspiring biomimetic active materials, due to their ability to change shape, size, and create their own morphology. We study a particular class of active gels, generated by polymerizing actin in the presence of cross-linkers and clusters of myosin as molecular motors, which exhibit large contractions. The relevant mechanics for these highly swollen gels is the result of the interplay between activity and liquid flow: gel activity yields a structural reorganization of the gel network and produces a flow of liquid that eventually exits from the gel boundary. This dynamics inherits lengthscales that are typical of the liquid flow processes. The analyses we present provide insights into the contraction dynamics, and they focus on the effects of the geometry on both gel velocity and fluid flow.

Geometry Does Impact on the Plane Strain Directions of the Human Left Ventricle, Irrespective of Disease.

The directions of primary strain lines of local deformation in Epicardial and Endocardial layers have been the subject of debate in recent years. Different methods led to different conclusions and a complete assessment of strain direction patterns in large and variable (in terms of pathology) cohorts of healthy and diseased patients is still lacking. Here, we use local deformation tensors in order to evaluate the angle of strain lines with respect to the horizontal circumferential direction in both Epi- and Endo-layers. We evaluated this on a large group of 193 subjects including 82 healthy control and 111 patients belonging to a great variety of pathological conditions. We found that Epicardial strain lines obliquely directed while those of Endocardium are almost circumferential. This result occurs irrespective of pathological condition. We propose that the geometric vinculum characterizing Endocardium and Epicardium in terms of different lever arm length and orientation of muscular fibers during contraction inescapably requires Endocardial strain lines to be circumferentially oriented and this is corroborated by experimental results. Further investigations on transmural structure of myocytes could couple results presented here in order to furnish additional experimental explanations.

Local and global energies for shape analysis in medical imaging.

In a previous contribution, a new Riemannian shape space, named TPS space, was introduced to perform statistics on shape data. This space was endowed with a Riemannian metric and a flat connection, with torsion, compatible with the given metric. This connection allows the definition of a Parallel Transport of the deformation compatible with the three-fold decomposition in spherical, deviatoric, and non-affine components. Such a parallel transport also conserves the Γ-energy, strictly related to the total elastic strain energy stored by the body in the original deformation. A new approach is here presented in order to calculate the bending energy on the body alone (body bending energy) and to restrict it exclusively within physical boundaries of objects involved in the deformation analysis. The novelty of this new procedure resides in the fact that we propose a new metric to be preserved during the TPS direct transport. This allows transporting the shape change more coherently with the mechanical meaning of the deformation. The geometry of the TPS space is then further discussed in order to better represent the relationship between the Γ-energy, the strain energy, and the so-called bending energy densities.

Swelling and shrinking in prestressed polymer gels: an incremental stress-diffusion analysis.

Polymer gels are porous fluid-saturated materials which can swell or shrink triggered by various stimuli. The swelling/shrinking-induced deformation can generate large stresses which may lead to the failure of the material. In the present research, a nonlinear stress-diffusion model is employed to investigate the stress and the deformation state arising in hydrated constrained polymer gels when subject to a varying chemical potential. Two different constraint configurations are taken into account: (i) elastic constraint along the thickness direction and (ii) plane elastic constraint. The first step entirely defines a compressed/tensed configuration. From there, an incremental chemo-mechanical analysis is presented. The derived model extends the classical linear poroelastic theory with respect to a prestressed configuration. Finally, the comparison between the analytical results obtained by the proposed model and a particular problem already discussed in literature for a stress-free gel membrane (one-dimensional test case) will highlight the relevance of the derived model.

Non-invasive prediction of genotype positive-phenotype negative in hypertrophic cardiomyopathy by 3D modern shape analysis.

NEW FINDINGS: What is the central question of this study? Can impaired deformational indicators for genotype positive for hypertrophic cardiomyopathy in subjects that do not exhibit a left-ventricular wall hypertrophy condition (G+LVH-) be determined using non-invasive 3D echocardiography? What is the main finding and its importance? Using 3D-STE and modern shape analysis, peculiar deformational impairments can be detected in G+LVH- subjects that can be classified with good accuracy. Moreover, the patterns of impairment are located mainly on the apical region in agreement with other evidence coming from previous biomechanical investigations. ABSTRACT: We propose a non-invasive procedure for predicting genotype positive for hypertrophic cardiomyopathy (HCM) in subjects that do not exhibit a left-ventricular wall hypertrophy condition (G+LVH-); the procedure is based on the enhanced analysis of medical imaging from 3D speckle tracking echocardiography (3D-STE). 3D-STE, due to its low quality images, has not been used so far to detect effectively the G+LVH- condition. Here, we post-processed echocardiographic images exploiting the tools of modern shape analysis, and we studied the motion of the left ventricle (LV) during an entire cycle. We enrolled 82 controls, 21 HCM patients and 11 G+LVH- subjects. We followed two steps: (i) we selected the most impaired regions of the LV by analysing its strains; and (ii) we used shape analysis on these regions to classify the subjects. The G+LVH- subjects showed different trajectories and deformational attributes. We found high classification performance in terms of area under the receiver operating characteristic curve (∼90), sensitivity (∼78) and specificity (∼79). Our results showed that (i) G+LVH- subjects present important deformational impairments relative to healthy controls and (ii) modern shape analysis can efficiently predict genotype by means of a non-invasive and inexpensive technique such as 3D-STE.

Torque-induced reorientation in active fibre-reinforced materials.

We introduce a continuum model for a fibre reinforced material in which the reference orientation of the fibre may evolve with time, under the influence of external stimuli. The model is formulated in the framework of large strain hyperelasticity and the kinematics of the continuum is described by both a position vector and by a remodelling tensor which, in the present context, is an orthogonal tensor representing the fibre reorientation process. By imposing suitable thermodynamical restrictions on the constitutive equation, we obtain an evolution equation of the remodelling tensor governed by the Eshelby torque, whose stationary solutions are studied in absence of any external source terms. It is shown that the fibres reorient themselves in a configuration that minimises the elastic energy and get aligned along a direction that may or may not be of principal strain. The explicit analysis of the Hessian of the strain energy density allows us to discriminate among the stationary solutions, which ones are stable. Examples are given for passive reorientation processes driven by applied strains or external boundary tractions. Applications of the proposed theory to biological tissues, nematic or magneto-electro active elastomers are foreseen.

The decomposition of deformation: New metrics to enhance shape analysis in medical imaging.

In landmarks-based Shape Analysis size is measured, in most cases, with Centroid Size. Changes in shape are decomposed in affine and non affine components. Furthermore the non affine component can be in turn decomposed in a series of local deformations (partial warps). If the extent of deformation between two shapes is small, the difference between Centroid Size and m-Volume increment is barely appreciable. In medical imaging applied to soft tissues bodies can undergo very large deformations, involving large changes in size. The cardiac example, analyzed in the present paper, shows changes in m-Volume that can reach the 60%. We show here that standard Geometric Morphometrics tools (landmarks, Thin Plate Spline, and related decomposition of the deformation) can be generalized to better describe the very large deformations of biological tissues, without losing a synthetic description. In particular, the classical decomposition of the space tangent to the shape space in affine and non affine components is enriched to include also the change in size, in order to give a complete description of the tangent space to the size-and-shape space. The proposed generalization is formulated by means of a new Riemannian metric describing the change in size as change in m-Volume rather than change in Centroid Size. This leads to a redefinition of some aspects of the Kendall's size-and-shape space without losing Kendall's original formulation. This new formulation is discussed by means of simulated examples using 2D and 3D platonic shapes as well as a real example from clinical 3D echocardiographic data. We demonstrate that our decomposition based approaches discriminate very effectively healthy subjects from patients affected by Hypertrophic Cardiomyopathy.

Driving water cavitation in a hydrogel cavity.

We study the dynamics of the dehydration process of a hydrogel with a cavity filled with water. We identify two transient phases: the first one dominated by an inflatable-balloon deformation mode, and the second by a suction effect, determining highly not homogeneous deformation modes of the hydrogel walls. This last phase triggers negative pressures into the cavity up to the typical values of water cavitation. An analysis of the factors allowing cavitation pressure to form inside the cavity is proposed, to allow for precise tuning of the key geometrical and material parameters.

Homeostatic Left Heart integration and disintegration links atrio-ventricular covariation's dyshomeostasis in Hypertrophic Cardiomyopathy.

Left ventricle and left atrium are and have been practically always analyzed separately in common clinically and non-clinically oriented cardiovascular investigations. Both classic and speckle tracking echocardiographic data contributed to the knowledge about deformational impairments occurring in systo-diastolic differences. Recently new trajectory based approaches allowed a greater awareness about the entire left ventricle or left atrium revolution and on their deficiencies that take place in presence of hypertrophic cardiomyopathy. However, surprisingly, the concomitant function of the two left heart chambers has not been analyzed for their geometrical/mechanical relationship. For the first time we study here, by acquiring left ventricle and left atrial geometries on the same heartbeat, the trajectory attributes of the entire left heart treated as a whole shape and the shape covariation of its two subunits. We contrasted healthy subjects with patients affected by hypertrophic cardiomyopathy. We found impaired left heart trajectory mainly in terms of orientation and size. More importantly, we found profound differences in the direction of morphological covariation of left ventricle and left atrium. These findings open to new perspectives in pathophysiological evaluation of different diseases by allowing the appreciation of concomitant functioning of both left heart whole geometry and of its two chambers.

Left Atrial trajectory impairment in Hypertrophic Cardiomyopathy disclosed by Geometric Morphometrics and Parallel Transport.

The analysis of full Left Atrium (LA) deformation and whole LA deformational trajectory in time has been poorly investigated and, to the best of our knowledge, seldom discussed in patients with Hypertrophic Cardiomyopathy. Therefore, we considered 22 patients with Hypertrophic Cardiomyopathy (HCM) and 46 healthy subjects, investigated them by three-dimensional Speckle Tracking Echocardiography, and studied the derived landmark clouds via Geometric Morphometrics with Parallel Transport. Trajectory shape and trajectory size were different in Controls versus HCM and their classification powers had high AUC (Area Under the Receiving Operator Characteristic Curve) and accuracy. The two trajectories were much different at the transition between LA conduit and booster pump functions. Full shape and deformation analyses with trajectory analysis enabled a straightforward perception of pathophysiological consequences of HCM condition on LA functioning. It might be worthwhile to apply these techniques to look for novel pathophysiological approaches that may better define atrio-ventricular interaction.

Geometry and mechanics of thin growing bilayers.

We investigate how thin sheets of arbitrary shapes morph under the isotropic in-plane expansion of their top surface, which may represent several stimuli such as nonuniform heating, local swelling and differential growth. Inspired by geometry, an analytical model is presented that rationalizes how the shape of the disk influences morphing, from the initial spherical bending to the final isometric limit. We introduce a new measure of slenderness that describes a sheet in terms of both thickness and plate shape. We find that the mean curvature of the isometric state is three fourths the natural curvature, which we verify by numerics and experiments. We finally investigate the emergence of a preferred direction of bending in the isometric state, guided by numerical analyses. The scalability of our model suggests that it is suitable to describe the morphing of sheets spanning several orders of magnitude.

Structure of tracheae and the functional implications for collapse in the American cockroach.

The tracheal tubes of insects are complex and heterogeneous composites with a microstructural organization that affects their function as pumps, valves, or static conduits within the respiratory system. In this study, we examined the microstructure of the primary thoracic tracheae of the American cockroach (Periplaneta americana) using a combination of scanning electron microscopy and light microscopy. The organization of the taenidia, which represents the primary source of structural reinforcement of the tracheae, was analyzed. We found that the taenidia were more disorganized in the regions of highest curvature of the tracheal tube. We also used a simple finite element model to explore the effect of cross-sectional shape and distribution of taenidia on the collapsibility of the tracheae. The eccentricity of the tracheal cross-section had a stronger effect on the collapse properties than did the distribution of taenidia. The combination of the macro-scale geometry, meso-scale heterogeneity, and microscale organization likely enables rhythmic tracheal compression during respiration, ultimately driving oxygen-rich air to cells and tissues throughout the insect body. The material design principles of these natural composites could potentially aid in the development of new bio-inspired microfluidic systems based on the differential collapse of tracheae-like networks.

Morphing of geometric composites via residual swelling.

Understanding and controlling the shape of thin, soft objects has been the focus of significant research efforts among physicists, biologists, and engineers in the last decade. These studies aim to utilize advanced materials in novel, adaptive ways such as fabricating smart actuators or mimicking living tissues. Here, we present the controlled growth-like morphing of 2D sheets into 3D shapes by preparing geometric composite structures that deform by residual swelling. The morphing of these geometric composites is dictated by both swelling and geometry, with diffusion controlling the swelling-induced actuation, and geometric confinement dictating the structure's deformed shape. Building on a simple mechanical analog, we present an analytical model that quantitatively describes how the Gaussian and mean curvatures of a thin disk are affected by the interplay among geometry, mechanics, and swelling. This model is in excellent agreement with our experiments and numerics. We show that the dynamics of residual swelling is dictated by a competition between two characteristic diffusive length scales governed by geometry. Our results provide the first 2D analog of Timoshenko's classical formula for the thermal bending of bimetallic beams - our generalization explains how the Gaussian curvature of a 2D geometric composite is affected by geometry and elasticity. The understanding conferred by these results suggests that the controlled shaping of geometric composites may provide a simple complement to traditional manufacturing techniques.

A new 4D trajectory-based approach unveils abnormal LV revolution dynamics in hypertrophic cardiomyopathy.

The assessment of left ventricular shape changes during cardiac revolution may be a new step in clinical cardiology to ease early diagnosis and treatment. To quantify these changes, only point registration was adopted and neither Generalized Procrustes Analysis nor Principal Component Analysis were applied as we did previously to study a group of healthy subjects. Here, we extend to patients affected by hypertrophic cardiomyopathy the original approach and preliminarily include genotype positive/phenotype negative individuals to explore the potential that incumbent pathology might also be detected. Using 3D Speckle Tracking Echocardiography, we recorded left ventricular shape of 48 healthy subjects, 24 patients affected by hypertrophic cardiomyopathy and 3 genotype positive/phenotype negative individuals. We then applied Generalized Procrustes Analysis and Principal Component Analysis and inter-individual differences were cleaned by Parallel Transport performed on the tangent space, along the horizontal geodesic, between the per-subject consensuses and the grand mean. Endocardial and epicardial layers were evaluated separately, different from many ecocardiographic applications. Under a common Principal Component Analysis, we then evaluated left ventricle morphological changes (at both layers) explained by first Principal Component scores. Trajectories' shape and orientation were investigated and contrasted. Logistic regression and Receiver Operating Characteristic curves were used to compare these morphometric indicators with traditional 3D Speckle Tracking Echocardiography global parameters. Geometric morphometrics indicators performed better than 3D Speckle Tracking Echocardiography global parameters in recognizing pathology both in systole and diastole. Genotype positive/phenotype negative individuals clustered with patients affected by hypertrophic cardiomyopathy during diastole, suggesting that incumbent pathology may indeed be foreseen by these methods. Left ventricle deformation in patients affected by hypertrophic cardiomyopathy compared to healthy subjects may be assessed by modern shape analysis better than by traditional 3D Speckle Tracking Echocardiography global parameters. Hypertrophic cardiomyopathy pathophysiology was unveiled in a new manner whereby also diastolic phase abnormalities are evident which is more difficult to investigate by traditional ecocardiographic techniques.

Anisotropic swelling of thin gel sheets.

We describe the anisotropic swelling within the Flory-Rehner thermodynamic model through an extension of the elastic component of the free-energy, which takes into account the oriented hampering of the swelling-induced deformations due to the presence of stiffer fibers. We also characterize the homogeneous free-swelling solutions of the corresponding anisotropic stress-diffusion problem, and discuss an asymptotic approximation of the key equations, which allows us to explicitly derive the anisotropic solution of the problem. We propose a proof-of-concept of our model, realizing thin bilayered gel sheets with layers having different anisotropic structures. In particular, for seedpod-like sheets, we observe and quantitatively measure the helicoid versus ribbon transition determined by the aspect ratio of the composite sheet.

Buckling dynamics of a solvent-stimulated stretched elastomeric sheet.

When stretched uniaxially, a thin elastic sheet may exhibit buckling. The occurrence of buckling depends on the geometrical properties of the sheet and the magnitude of the applied strain. Here we show that an elastomeric sheet initially stable under uniaxial stretching can destabilize when exposed to a solvent that swells the elastomer. We demonstrate experimentally and computationally that the features of the buckling pattern depend on the magnitude of stretching, and this observation offers a new way for controlling the shape of a swollen homogeneous thin sheet.

4D-analysis of left ventricular heart cycle using procrustes motion analysis.

The aim of this study is to investigate human left ventricular heart morphological changes in time among 17 healthy subjects. Preliminarily, 2 patients with volumetric overload due to aortic insufficiency were added to our analyses. We propose a special strategy to compare the shape, orientation and size of cardiac cycle's morphological trajectories in time. We used 3D data obtained by Speckle Tracking Echocardiography in order to detect semi-automated and homologous landmarks clouds as proxies of left ventricular heart morphology. An extended Geometric Morphometrics toolkit in order to distinguish between intra- and inter-individual shape variations was used. Shape of trajectories with inter-individual variation were compared under the assumption that trajectories attributes, estimated at electrophysiologically homologous times are expressions of left ventricular heart function. We found that shape analysis as commonly applied in Geometric Morphometrics studies fails in identifying a proper morpho-space to compare the shape of morphological trajectories in time. To overcome this problem, we performed a special type of Riemannian Parallel Transport, called "linear shift". Whereas the two patients with aortic insufficiency were not differentiated in the static shape analysis from the healthy subjects, they set apart significantly in the analyses of motion trajectory's shape and orientation. We found that in healthy subjects, the variations due to inter-individual morphological differences were not related to shape and orientation of morphological trajectories. Principal Component Analysis showed that volumetric contraction, torsion and twist are differently distributed on different axes. Moreover, global shape change appeared to be more correlated with endocardial shape change than with the epicardial one. Finally, the total shape variation occurring among different subjects was significantly larger than that observable across properly defined morphological trajectories.