Automatic Segmentation of Mechanically Inhomogeneous Tissues Based on Deformation Gradient Jump

Evaluation of an Inverse Method for Quantifying Spatially Variable Mechanics

Soft biological tissues often function as highly deformable membranes in vivo and exhibit impressive mechanical behavior effectively characterized by planar biaxial testing. The Generalized Anisotropic Inverse Mechanics (GAIM) method links full-field deformations and boundary forces from mechanical testing to quantify material properties of soft, anisotropic, heterogeneous tissues. In this study, we introduced an orthotropic constraint to GAIM to improve the quality and physical significance of its mechanical characterizations. We evaluated the updated GAIM method using simulated and experimental biaxial testing datasets obtained from soft tissue analogs (PDMS and TissueMend) with well-defined mechanical properties. GAIM produced stiffnesses (first Kelvin moduli, K1) that agreed well with previously published Young's moduli of PDMS samples. It also matched the stiffness moduli determined via uniaxial testing for TissueMend, a collagen-rich patch intended for tendon repair. We then conducted the first biaxial testing of TissueMend and confirmed that the sample was mechanically anisotropic via a relative anisotropy metric produced by GAIM. Next, we demonstrated the benefits of full-field laser micrometry in distinguishing between spatial variations in thickness and stiffness. Finally, we conducted an analysis to verify that results were independent of partitioning scheme. The success of the newly implemented constraints on GAIM suggests notable potential for applying this tool to soft tissues, particularly following the onset of pathologies that induce mechanical and structural heterogeneities.

Segmenting mechanically heterogeneous domains via unsupervised learning

From biological organs to soft robotics, highly deformable materials are essential components of natural and engineered systems. These highly deformable materials can have heterogeneous material properties, and can experience heterogeneous deformations with or without underlying material heterogeneity. Many recent works have established that computational modeling approaches are well suited for understanding and predicting the consequences of material heterogeneity and for interpreting observed heterogeneous strain fields. In particular, there has been significant work toward developing inverse analysis approaches that can convert observed kinematic quantities (e.g., displacement, strain) to material properties and mechanical state. Despite the success of these approaches, they are not necessarily generalizable and often rely on tight control and knowledge of boundary conditions. Here, we will build on the recent advances (and ubiquity) of machine learning approaches to explore alternative approaches to detect patterns in heterogeneous material properties and mechanical behavior. Specifically, we will explore unsupervised learning approaches to clustering and ensemble clustering to identify heterogeneous regions. Overall, we find that these approaches are effective, yet limited in their abilities. Through this initial exploration (where all data and code are published alongside this manuscript), we set the stage for future studies that more specifically adapt these methods to mechanical data.

The Lumbar Facet Capsular Ligament Becomes More Anisotropic and the Fibers Become Stiffer With Intervertebral Disc and Facet Joint Degeneration

Degeneration of the lumbar spine, and especially how that degeneration may lead to pain, remains poorly understood. In particular, the mechanics of the facet capsular ligament may contribute to low back pain, but the mechanical changes that occur in this ligament with spinal degeneration are unknown. Additionally, the highly nonlinear, heterogeneous, and anisotropic nature of the facet capsular ligament makes understanding mechanical changes more difficult. Clinically, magnetic resonance imaging (MRI)-based signs of degeneration in the facet joint and the intervertebral disc (IVD) correlate. Therefore, this study examined how the nonlinear, heterogeneous mechanics of the facet capsular ligament change with degeneration of the lumbar spine as characterized using MRI. Cadaveric human spines were imaged via MRI, and the L2-L5 facet joints and IVDs were scored using the Fujiwara and Pfirrmann grading systems. Then, the facet capsular ligament was isolated and biaxially loaded. The nonlinear mechanical properties of the ligament were obtained using a nonlinear generalized anisotropic inverse mechanics analysis (nGAIM). Then a Holzapfel-Gasser-Ogden (HGO) model was fit to the stress-strain data obtained from nGAIM. The facet capsular ligament is stiffer and more anisotropic at larger Pfirrmann grades and higher Fujiwara scores than at lower grades and scores. Analysis of ligament heterogeneity showed all tissues are highly heterogeneous, but no distinct spatial patterns of heterogeneity were found. These results show that degeneration of the lumbar spine including the facet capsular ligament appears to be occurring as a whole joint phenomenon and advance our understanding of lumbar spine degeneration.

Characterizing Tissue Remodeling and Mechanical Heterogeneity in Cerebral Aneurysms

Accurately assessing the complex tissue mechanics of cerebral aneurysms (CAs) is critical for elucidating how CAs grow and whether that growth will lead to rupture. The factors that have been implicated in CA progression - blood flow dynamics, immune infiltration, and extracellular matrix remodeling - all occur heterogeneously throughout the CA. Thus, it stands to reason that the mechanical properties of CAs are also spatially heterogeneous. Here, we present a new method for characterizing the mechanical heterogeneity of human CAs using generalized anisotropic inverse mechanics, which uses biaxial stretching experiments and inverse analyses to determine the local Kelvin moduli and principal alignments within the tissue. Using this approach, we find that there is significant mechanical heterogeneity within a single acquired human CA. These results were confirmed using second harmonic generation imaging of the CA's fiber architecture and a correlation was observed. This approach provides a single-step method for determining the complex heterogeneous mechanics of CAs, which has important implications for future identification of metrics that can improve accuracy in prediction risk of rupture.

A nonlinear anisotropic inverse method for computational dissection of inhomogeneous planar tissues

Quantification of the mechanical behavior of soft tissues is challenging due to their anisotropic, heterogeneous, and nonlinear nature. We present a method for the 'computational dissection' of a tissue, by which we mean the use of computational tools both to identify and to analyze regions within a tissue sample that have different mechanical properties. The approach employs an inverse technique applied to a series of planar biaxial experimental protocols. The aggregated data from multiple protocols provide the basis for (1) segmentation of the tissue into regions of similar properties, (2) linear analysis for the small-strain behavior, assuming uniform, linear, anisotropic behavior within each region, (3) subsequent nonlinear analysis following each individual experimental protocol path and using local linear properties, and (4) construction of a strain energy data set W(E) at every point in the material by integrating the differential stress-strain functions along each strain path. The approach has been applied to simulated data and captures not only the general nonlinear behavior but also the regional differences introduced into the simulated tissue sample.

Marker-Free Tracking of Facet Capsule Motion Using Polarization-Sensitive Optical Coherence Tomography

We proposed and tested a method by which surface strains of biological tissues can be captured without the use of fiducial markers by instead, utilizing the inherent structure of the tissue. We used polarization-sensitive optical coherence tomography (PS OCT) to obtain volumetric data through the thickness and across a partial surface of the lumbar facet capsular ligament during three cases of static bending. Reflectivity and phase retardance were calculated from two polarization channels, and a power spectrum analysis was performed on each a-line to extract the dominant banding frequency (a measure of degree of fiber alignment) through the maximum value of the power spectrum (maximum power). Maximum powers of all a-lines for each case were used to create 2D visualizations, which were subsequently tracked via digital image correlation. In-plane strains were calculated from measured 2D deformations and converted to 3D surface strains by including out-of-plane motion obtained from the PS OCT image. In-plane strains correlated with 3D strains (R(2) ≥ 0.95). Using PS OCT for marker-free motion tracking of biological tissues is a promising new technique because it relies on the structural characteristics of the tissue to monitor displacement instead of external fiducial markers.