Quantifying deformations and strains in human intervertebral discs using Digital Volume Correlation combined with MRI (DVC-MRI)

Full-Field Strain Measurements of the Muscle-Tendon Junction Using X-ray Computed Tomography and Digital Volume Correlation

We report, for the first time, the full-field 3D strain distribution of the muscle-tendon junction (MTJ). Understanding the strain distribution at the junction is crucial for the treatment of injuries and to predict tear formation at this location. Three-dimensional full-field strain distribution of mouse MTJ was measured using X-ray computer tomography (XCT) combined with digital volume correlation (DVC) with the aim of understanding the mechanical behavior of the junction under tensile loading. The interface between the Achilles tendon and the gastrocnemius muscle was harvested from adult mice and stained using 1% phosphotungstic acid in 70% ethanol. In situ XCT combined with DVC was used to image and compute strain distribution at the MTJ under a tensile load (2.4 N). High strain measuring 120,000 µε, 160,000 µε, and 120,000 µε for the first principal stain (εp1), shear strain (γ), and von Mises strain (εVM), respectively, was measured at the MTJ and these values reduced into the body of the muscle or into the tendon. Strain is concentrated at the MTJ, which is at risk of being damaged in activities associated with excessive physical activity.

The influence of geometry on intervertebral disc stiffness

Geometry plays an important role in intervertebral disc (IVD) mechanics. Previous computational studies have found a link between IVD geometry and stiffness. However, few experimental studies have investigated this link, possibly due to difficulties in non-destructively quantifying internal geometric features. Recent advances in ultra-high resolution MRI provides the opportunity to visualise IVD features in unprecedented detail. This study aimed to quantify 3D human IVD geometries using 9.4 T MRIs and to investigate correlations between geometric variations and IVD stiffness. Thirty human lumbar motion segments (fourteen non-degenerate and sixteen degenerate) were scanned using a 9.4 T MRI and geometric parameters were measured. A 1kN compressive load was applied to each motion segment and stiffness was calculated. Degeneration caused a reduction (p < 0.05) in IVD height, a decreased nucleus-annulus area ratio, and a 1.6 ± 3.0 mm inward collapse of the inner annulus. The IVD height, anteroposterior (AP) width, lateral width, cross-sectional area, nucleus-annulus boundary curvature, and nucleus-annulus area ratio had a significant (p < 0.05) influence on IVD stiffness. Linear relationships (p < 0.05, r > 0.47) were observed between these geometric features and IVD compressive stiffness and a multivariate regression model was generated to enable stiffness to be predicted from features observable on clinical imaging (stiffness, N/mm = 6062 - (61.2 × AP width, mm) - (169.2 × IVD height, mm)). This study advances our understanding of disc structure-function relationships and how these change with degeneration, which can be used to both generate and validate more realistic computational models.

Quantifying internal intervertebral disc strains to assess nucleus replacement device designs: a digital volume correlation and ultra-high-resolution MRI study

Introduction: Nucleus replacement has been proposed as a treatment to restore biomechanics and relieve pain in degenerate intervertebral discs (IVDs). Multiple nucleus replacement devices (NRDs) have been developed, however, none are currently used routinely in clinic. A better understanding of the interactions between NRDs and surrounding tissues may provide insight into the causes of implant failure and provide target properties for future NRD designs. The aim of this study was to non-invasively quantify 3D strains within the IVD through three stages of nucleus replacement surgery: intact, post-nuclectomy, and post-treatment. Methods: Digital volume correlation (DVC) combined with 9.4T MRI was used to measure strains in seven human cadaveric specimens (42 ± 18 years) when axially compressed to 1 kN. Nucleus material was removed from each specimen creating a cavity that was filled with a hydrogel-based NRD. Results: Nucleus removal led to loss of disc height (12.6 ± 4.4%, p = 0.004) which was restored post-treatment (within 5.3 ± 3.1% of the intact state, p > 0.05). Nuclectomy led to increased circumferential strains in the lateral annulus region compared to the intact state (-4.0 ± 3.4% vs. 1.7 ± 6.0%, p = 0.013), and increased maximum shear strains in the posterior annulus region (14.6 ± 1.7% vs. 19.4 ± 2.6%, p = 0.021). In both cases, the NRD was able to restore these strain values to their intact levels (p ≥ 0.192). Discussion: The ability of the NRD to restore IVD biomechanics and some strain types to intact state levels supports nucleus replacement surgery as a viable treatment option. The DVC-MRI method used in the present study could serve as a useful tool to assess future NRD designs to help improve performance in future clinical trials.

A novel tool to quantify in vivo lumbar spine kinematics and 3D intervertebral disc strains using clinical MRI

Medical imaging modalities that calculate tissue morphology alone cannot provide direct information regarding the mechanical behaviour of load-bearing musculoskeletal organs. Accurate in vivo measurement of spine kinematics and intervertebral disc (IVD) strains can provide important information regarding the mechanical behaviour of the spine, help to investigate the effects of injuries on the mechanics of the spine, and assess the effectiveness of treatments. Additionally, strains can serve as a functional biomechanical marker for detecting normal and pathologic tissues. We hypothesised that combining digital volume correlation (DVC) with 3T clinical MRI can provide direct information regarding the mechanics of the spine. Here, we have developed a novel non-invasive tool for in vivo displacement and strain measurement within the human lumbar spine and we used this tool to calculate lumbar kinematics and IVD strains in six healthy subjects during lumbar extension. The proposed tool enabled spine kinematics and IVD strains to be measured with errors that did not exceed 0.17 mm and 0.5%, respectively. The findings of the kinematics study identified that during extension the lumbar spine of healthy subjects experiences total 3D translations ranging from 1 mm to 4.5 mm for different vertebral levels. The findings of strain analysis identified that the average of the maximum tensile, compressive, and shear strains for different lumbar levels during extension ranged from 3.5% to 7.2%. This tool can provide base-line data that can be used to describe the mechanical environment of healthy lumbar spine, which can help clinicians manage preventative treatments, define patient-specific treatments, and to monitor the effectiveness of surgical and non-surgical interventions.

Image quality and scan time optimisation for in situ phase contrast x-ray tomography of the intervertebral disc

In-line phase contrast synchrotron tomography combined with in situ mechanical loading enables the characterisation of soft tissue micromechanics via digital volume correlation (DVC) within whole organs. Optimising scan time is important for reducing radiation dose from multiple scans and to limit sample movement during acquisition. Also, although contrasted edges provided by in-line phase contrast tomography of soft tissues are useful for DVC, the effect of phase contrast imaging on its accuracy has yet to be investigated. Due to limited time at synchrotron facilities, scan parameters are often decided during imaging and their effect on DVC accuracy is not fully understood. Here, we used previously published data of intervertebral disc phase contrast tomography to evaluate the influence of i) fibrous image texture, ii) number of projections, iii) tomographic reconstruction method, and iv) phase contrast propagation distance on DVC results. A greater understanding of how image texture influences optimal DVC tracking was obtained by visualising objective function mapping, enabling tracking inaccuracies to be identified. When reducing the number of projections, DVC was minimally affected by image high frequency noise but with a compromise in accuracy. Iterative reconstruction methods improved image signal-to-noise and consequently significantly lowered DVC displacement uncertainty. Propagation distance was shown to affect DVC accuracy. Consistent DVC results were achieved within a propagation distance range which provided contrast to the smallest scale features, where; too short a distance provided insufficient features to track, whereas too long led to edge effect inconsistencies, particularly at greater deformations. Although limited to a single sample type and image setup, this study provides general guidelines for future investigations when optimising image quality and scan times for in situ phase contrast x-ray tomography of fibrous connective tissues.

A practical guide for in situ mechanical testing of musculoskeletal tissues using synchrotron tomography

Musculoskeletal tissues are complex hierarchical materials where mechanical response is linked to structural and material properties at different dimensional levels. Therefore, high-resolution three-dimensional tomography is very useful for assessing tissue properties at different scales. In particular, Synchrotron Radiation micro-Computed Tomography (SR-microCT) has been used in several applications to analyze the structure of bone and biomaterials. In the past decade the development of digital volume correlation (DVC) algorithms applied to SR-microCT images and its combination with in situ mechanical testing (four-dimensional imaging) have allowed researchers to visualise, for the first time, the deformation of musculoskeletal tissues and their interaction with biomaterials under different loading scenarios. However, there are several experimental challenges that make these measurements difficult and at high risk of failure. Challenges relate to sample preparation, imaging parameters, loading setup, accumulated tissue damage for multiple tomographic acquisitions, reconstruction methods and data processing. Considering that access to SR-microCT facilities is usually associated with bidding processes and long waiting times, the failure of these experiments could notably slow down the advancement of this research area and reduce its impact. Many of the experimental failures can be avoided with increased experience in performing the tests and better guidelines for preparation and execution of these complex experiments; publication of negative results could help interested researchers to avoid recurring mistakes. Therefore, the goal of this article is to highlight the potential and pitfalls in the design and execution of in situ SR-microCT experiments, involving multiple scans, of musculoskeletal tissues for the assessment of their structural and/or mechanical properties. The advice and guidelines that follow should improve the success rate of this type of experiment, allowing the community to reach higher impact more efficiently.

Digital volume correlation for the characterization of musculoskeletal tissues: Current challenges and future developments

Biological tissues are complex hierarchical materials, difficult to characterise due to the challenges associated to the separation of scale and heterogeneity of the mechanical properties at different dimensional levels. The Digital Volume Correlation approach is the only image-based experimental approach that can accurately measure internal strain field within biological tissues under complex loading scenarios. In this minireview examples of DVC applications to study the deformation of musculoskeletal tissues at different dimensional scales are reported, highlighting the potential and challenges of this relatively new technique. The manuscript aims at reporting the wide breath of DVC applications in the past 2 decades and discuss future perspective for this unique technique, including fast analysis, applications on soft tissues, high precision approaches, and clinical applications.

High resolution three-dimensional strain measurements in human articular cartilage

An unresolved challenge in osteoarthritis research is characterising the localised intra-tissue mechanical response of articular cartilage. The aim of this study was to explore whether laboratory micro-computed tomography (micro-CT) and digital volume correlation (DVC) permit non-destructive quantification of three-dimensional (3D) strain fields in human articular cartilage. Human articular cartilage specimens were harvested from the knee, mounted into a loading device and imaged in the unloaded and loaded states using a micro-CT scanner. Strain was measured throughout the cartilage volume using the micro-CT image data and DVC analysis. The volumetric DVC-measured strain was within 5% of the known applied strain. Variation in strain distribution between the superficial, middle and deep zones was observed, consistent with the different architecture of the material in these locations. These results indicate DVC method may be suitable for calculating strain in human articular cartilage.

Sensitivity of Intervertebral Disc Finite Element Models to Internal Geometric and Non-geometric Parameters

Finite element models are useful for investigating internal intervertebral disc (IVD) behaviours without using disruptive experimental techniques. Simplified geometries are commonly used to reduce computational time or because internal geometries cannot be acquired from CT scans. This study aimed to (1) investigate the effect of altered geometries both at endplates and the nucleus-anulus boundary on model response, and (2) to investigate model sensitivity to material and geometric inputs, and different modelling approaches (graduated or consistent fibre bundle angles and glued or cohesive inter-lamellar contact). Six models were developed from 9.4 T MRIs of bovine IVDs. Models had two variations of endplate geometry (a simple curved profile from the centre of the disc to the periphery, and precise geometry segmented from MRIs), and three variations of NP-AF boundary (linear, curved, and segmented). Models were subjected to axial compressive loading (to 0.86 mm at a strain rate of 0.1/s) and the effect on stiffness and strain distributions, and the sensitivity to modelling approaches was investigated. The model with the most complex geometry (segmented endplates, curved NP-AF boundary) was 3.1 times stiffer than the model with the simplest geometry (curved endplates, linear NP-AF boundary), although this difference may be exaggerated since segmenting the endplates in the complex geometry models resulted in a shorter average disc height. Peak strains were close to the endplates at locations of high curvature in the segmented endplate models which were not captured in the curved endplate models. Differences were also seen in sensitivity to material properties, graduated fibre angles, cohesive rather than glued inter-lamellar contact, and NP:AF ratios. These results show that FE modellers must take care to ensure geometries are realistic so that load is distributed and passes through IVDs accurately.

Experimental and Computational Comparison of Intervertebral Disc Bulge for Specimen-Specific Model Evaluation Based on Imaging

Finite element modelling of the spinal unit is a promising preclinical tool to assess the biomechanical outcome of emerging interventions. Currently, most models are calibrated and validated against range of motion and rarely directly against soft-tissue deformation. The aim of this contribution was to develop an in vitro methodology to measure disc bulge and assess the ability of different specimen-specific modelling approaches to predict disc bulge. Bovine bone-disc-bone sections (N = 6) were prepared with 40 glass markers on the intervertebral disc surface. These were initially magnetic resonance (MR)-imaged and then sequentially imaged using peripheral-qCT under axial compression of 1 mm increments. Specimen-specific finite-element models were developed from the CT data, using three different methods to represent the nucleus pulposus geometry with and without complementary use of the MR images. Both calibrated specimen-specific and averaged compressive material properties for the disc tissues were investigated. A successful methodology was developed to quantify the disc bulge in vitro, enabling observation of surface displacement on qCT. From the finite element model results, no clear advantage was found in using geometrical information from the MR images in terms of the models' ability to predict stiffness or disc bulge for bovine intervertebral disc.

Toward Patient Specific Models of Pediatric IVDs: A Parametric Study of IVD Mechanical Properties

Patient specific finite element (FE) modeling of the pediatric spine is an important challenge which offers to revolutionize the treatment of pediatric spinal pathologies, for example adolescent idiopathic scoliosis (AIS). In particular, modeling of the intervertebral disc (IVD) is a unique challenge due to its structural and mechanical complexity. This is compounded by limited ability to non-invasively interrogate key mechanical parameters of a patient's IVD. In this work, we seek to better understand the link between mechanical properties and mechanical behavior of patient specific FE models of the pediatric lumbar spine. A parametric study of IVD parameter was conducted, coupled with insights from current knowledge of the pediatric IVD. In particular, the combined effects of parameters was investigated. Recommendations are made toward areas of importance in patient specific FE modeling of the pediatric IVD. In particular, collagen fiber bundles of the IVD are found to dominate IVD mechanical behavior and are thus recommended as an area of primary focus for patient specific FE models. In addition, areas requiring further experimental research are identified. This work provides a valuable building block toward the development of patient specific models of the pediatric spine.

The Effect of Degeneration on Internal Strains and the Mechanism of Failure in Human Intervertebral Discs Analyzed Using Digital Volume Correlation (DVC) and Ultra-High Field MRI

The intervertebral disc (IVD) plays a main role in absorbing and transmitting loads within the spinal column. Degeneration alters the structural integrity of the IVDs and causes pain, especially in the lumbar region. The objective of this study was to investigate non-invasively the effect of degeneration on human 3D lumbar IVD strains (n = 8) and the mechanism of spinal failure (n = 10) under pure axial compression using digital volume correlation (DVC) and 9.4 Tesla magnetic resonance imaging (MRI). Degenerate IVDs had higher (p < 0.05) axial strains (58% higher), maximum 3D compressive strains (43% higher), and maximum 3D shear strains (41% higher), in comparison to the non-degenerate IVDs, particularly in the lateral and posterior annulus. In both degenerate and non-degenerate IVDs, peak tensile and shear strains were observed close to the endplates. Inward bulging of the inner annulus was observed in all degenerate IVDs causing an increase in the AF compressive, tensile, and shear strains at the site of inward bulge, which may predispose it to circumferential tears (delamination). The endplate is the spine's “weak link” in pure axial compression, and the mechanism of human vertebral fracture is associated with disc degeneration. In non-degenerate IVDs the locations of failure were close to the endplate centroid, whereas in degenerate IVDs they were in peripheral regions. These findings advance the state of knowledge on mechanical changes during degeneration of the IVD, which help reduce the risk of injury, optimize treatments, and improve spinal implant designs. Additionally, these new data can be used to validate computational models.

Correlating Local Volumetric Tissue Strains with Global Lung Mechanics Measurements

The mechanics of breathing is a fascinating and vital process. The lung has complexities and subtle heterogeneities in structure across length scales that influence mechanics and function. This study establishes an experimental pipeline for capturing alveolar deformations during a respiratory cycle using synchrotron radiation micro-computed tomography (SR-micro-CT). Rodent lungs were mechanically ventilated and imaged at various time points during the respiratory cycle. Pressure-Volume (P-V) characteristics were recorded to capture any changes in overall lung mechanical behaviour during the experiment. A sequence of tomograms was collected from the lungs within the intact thoracic cavity. Digital volume correlation (DVC) was used to compute the three-dimensional strain field at the alveolar level from the time sequence of reconstructed tomograms. Regional differences in ventilation were highlighted during the respiratory cycle, relating the local strains within the lung tissue to the global ventilation measurements. Strains locally reached approximately 150% compared to the averaged regional deformations of approximately 80-100%. Redistribution of air within the lungs was observed during cycling. Regions which were relatively poorly ventilated (low deformations compared to its neighbouring region) were deforming more uniformly at later stages of the experiment (consistent with its neighbouring region). Such heterogenous phenomena are common in everyday breathing. In pathological lungs, some of these non-uniformities in deformation behaviour can become exaggerated, leading to poor function or further damage. The technique presented can help characterize the multiscale biomechanical nature of a given pathology to improve patient management strategies, considering both the local and global lung mechanics.

In Vivo Deformation and Strain Measurements in Human Bone Using Digital Volume Correlation (DVC) and 3T Clinical MRI

Strains within bone play an important role in the remodelling process and the mechanisms of fracture. The ability to assess these strains in vivo can provide clinically relevant information regarding bone health, injury risk, and can also be used to optimise treatments. In vivo bone strains have been investigated using multiple experimental techniques, but none have quantified 3D strains using non-invasive techniques. Digital volume correlation based on clinical MRI (DVC-MRI) is a non-invasive technique that has the potential to achieve this. However, before it can be implemented, uncertainties associated with the measurements must be quantified. Here, DVC-MRI was evaluated to assess its potential to measure in vivo strains in the talus. A zero-strain test (two repeated unloaded scans) was conducted using three MRI sequences, and three DVC approaches to quantify errors and to establish optimal settings. With optimal settings, strains could be measured with a precision of 200 με and accuracy of 480 με for a spatial resolution of 7.5 mm, and a precision of 133 με and accuracy of 251 με for a spatial resolution of 10 mm. These results demonstrate that this technique has the potential to measure relevant levels of in vivo bone strain and to be used for a range of clinical applications.

Developing an Analogue Residual Limb for Comparative DVC Analysis of Transtibial Prosthetic Socket Designs

Personalised prosthetic sockets are fabricated by expert clinicians in a skill- and experience-based process, with research providing tools to support evidence-based practice. We propose that digital volume correlation (DVC) may offer a deeper understanding of load transfer from prosthetic sockets into the residual limb, and tissue injury risk. This study’s aim was to develop a transtibial amputated limb analogue for volumetric strain estimation using DVC, evaluating its ability to distinguish between socket designs. A soft tissue analogue material was developed, comprising silicone elastomer and sand particles as fiducial markers for image correlation. The material was cast to form an analogue residual limb informed by an MRI scan of a person with transtibial amputation, for whom two polymer check sockets were produced by an expert prosthetist. The model was micro-CT scanned according to (i) an unloaded noise study protocol and (ii) a case study comparison between the two socket designs, loaded to represent two-legged stance. The scans were reconstructed to give 108 µm voxels. The DVC noise study indicated a 64 vx subvolume and 50% overlap, giving better than 0.32% strain sensitivity, and ~3.5 mm spatial resolution of strain. Strain fields induced by the loaded sockets indicated tensile, compressive and shear strain magnitudes in the order of 10%, with a high signal:noise ratio enabling distinction between the two socket designs. DVC may not be applicable for socket design in the clinical setting, but does offer critical 3D strain information from which existing in vitro and in silico tools can be compared and validated to support the design and manufacture of prosthetic sockets, and enhance the biomechanical understanding of the load transfer between the limb and the prosthesis.

Quantifying 3D Strain in Scaffold Implants for Regenerative Medicine

Regenerative medicine solutions require thoughtful design to elicit the intended biological response. This includes the biomechanical stimulus to generate an appropriate strain in the scaffold and surrounding tissue to drive cell lineage to the desired tissue. To provide appropriate strain on a local level, new generations of scaffolds often involve anisotropic spatially graded mechanical properties that cannot be characterised with traditional materials testing equipment. Volumetric examination is possible with three-dimensional (3D) imaging, in situ loading and digital volume correlation (DVC). Micro-CT and DVC were utilised in this study on two sizes of 3D-printed inorganic/organic hybrid scaffolds (n = 2 and n = 4) with a repeating homogenous structure intended for cartilage regeneration. Deformation was observed with a spatial resolution of under 200 µm whilst maintaining displacement random errors of 0.97 µm, strain systematic errors of 0.17% and strain random errors of 0.031%. Digital image correlation (DIC) provided an analysis of the external surfaces whilst DVC enabled localised strain concentrations to be examined throughout the full 3D volume. Strain values derived using DVC correlated well against manually calculated ground-truth measurements (R2 = 0.98, n = 8). The technique ensures the full 3D micro-mechanical environment experienced by cells is intimately considered, enabling future studies to further examine scaffold designs for regenerative medicine.