In vivo intervertebral disc deformation: intratissue strain patterns within adjacent discs during flexion-extension

Intervertebral Disc Elastography to Relate Shear Modulus and Relaxometry in Compression and Bending

Intervertebral disc degeneration is the most recognized cause of low back pain, characterized by the decline of tissue structure and mechanics. Image-based mechanical parameters (e.g., strain, stiffness) may provide an ideal assessment of disc function that is lost with degeneration but unfortunately remains underdeveloped. Moreover, it is unknown whether strain or stiffness of the disc may be predicted by MRI relaxometry (e.g. T1 or T2), an increasingly accepted quantitative measure of disc structure. In this study, we quantified T1 and T2 relaxation times and in-plane strains using displacement-encoded MRI within the disc under physiological levels of compression and bending. We then estimated shear modulus in orthogonal image planes and compared these values to relaxation times and strains within regions of the disc. Intratissue strain depended on the loading mode, and shear modulus in the nucleus pulposus was typically an order of magnitude lower than the annulus fibrosis, except in bending, where the apparent stiffness depended on the loading. Relative shear moduli estimated from strain data derived under compression generally did not correspond with those from bending experiments, with no correlations in the sagittal plane and only 4 of 15 regions correlated in the coronal plane, suggesting that future inverse models should incorporate multiple loading conditions. Strain imaging and strain-based estimation of material properties may serve as imaging biomarkers to distinguish healthy and diseased discs. Additionally, image-based elastography and relaxometry may be viewed as complementary measures of disc structure and function to assess degeneration in longitudinal studies.

Multi-frame biomechanical and relaxometry analysis during in vivo loading of the human knee by spiral dualMRI and compressed sensing

PURPOSE

Knee cartilage experiences repetitive loading during physical activities, which is altered during the pathogenesis of diseases like osteoarthritis. Analyzing the biomechanics during motion provides a clear understanding of the dynamics of cartilage deformation and may establish essential imaging biomarkers of early-stage disease. However, in vivo biomechanical analysis of cartilage during rapid motion is not well established.

METHODS

We used spiral displacement encoding with stimulated echoes (DENSE) MRI on in vivo human tibiofemoral cartilage during cyclic varus loading (0.5 Hz) and used compressed sensing on the k-space data. The applied compressive load was set for each participant at 0.5 times body weight on the medial condyle. Relaxometry methods were measured on the cartilage before (T1ρ , T2 ) and after (T1ρ ) varus load.

RESULTS

Displacement and strain maps showed a gradual shift of displacement and strain in time. Compressive strain was observed in the medial condyle cartilage and shear strain was roughly half of the compressive strain. Male participants had more displacement in the loading direction compared to females, and T1ρ values did not change after cyclic varus load. Compressed sensing reduced the scanning time up to 25% to 40% when comparing the displacement maps and substantially lowered the noise levels.

CONCLUSION

These results demonstrated the ease of which spiral DENSE MRI could be applied to clinical studies because of the shortened imaging time, while quantifying realistic cartilage deformations that occur through daily activities and that could serve as biomarkers of early osteoarthritis.

Physiological and degenerative loading of bovine intervertebral disc in a bioreactor: A finite element study of complex motions

Intervertebral disc (IVD) degeneration and regenerative therapies are commonly studied in organ-culture experiments with uniaxial compressive loading. Recently, in our laboratory, we established a bioreactor system capable of applying loads in six degrees-of-freedom (DOF) to bovine IVDs, which replicates more closely the complex multi-axial loading of the IVD in vivo. However, the magnitudes of loading that are physiological (able to maintain cell viability) or mechanically degenerative are unknown for load cases combining several DOFs. This study aimed to establish physiological and degenerative levels of maximum principal strains and stresses in the bovine IVD tissue and to investigate how they are achieved under complex load cases related to common daily activities. The physiological and degenerative levels of maximum principal strains and stresses were determined via finite element (FE) analysis of bovine IVD subjected to experimentally established physiological and degenerative compressive loading protocols. Then, complex load cases, such as a combination of compression + flexion + torsion, were applied on the FE-model with increasing magnitudes of loading to discover when physiological and degenerative tissue strains and stresses were reached. When applying 0.1 MPa of compression and ±2-3° of flexion and ±1-2° of torsion the investigated mechanical parameters remained at physiological levels, but with ±6-8° of flexion in combination with ±2-4° of torsion, the stresses in the outer annulus fibrosus (OAF) exceeded degenerative levels. In the case of compression + flexion + torsion, the mechanical degeneration likely initiates at the OAF when loading magnitudes are high enough. The physiological and degenerative magnitudes can be used as guidelines for bioreactor experiments with bovine IVDs.

The Digital Twin: A Potential Solution for the Personalized Diagnosis and Treatment of Musculoskeletal System Diseases

Due to the high prevalence and rates of disability associated with musculoskeletal system diseases, more thorough research into diagnosis, pathogenesis, and treatments is required. One of the key contributors to the emergence of diseases of the musculoskeletal system is thought to be changes in the biomechanics of the human musculoskeletal system. However, there are some defects concerning personal analysis or dynamic responses in current biomechanical research methodologies. Digital twin (DT) was initially an engineering concept that reflected the mirror image of a physical entity. With the application of medical image analysis and artificial intelligence (AI), it entered our lives and showed its potential to be further applied in the medical field. Consequently, we believe that DT can take a step towards personalized healthcare by guiding the design of industrial personalized healthcare systems. In this perspective article, we discuss the limitations of traditional biomechanical methods and the initial exploration of DT in musculoskeletal system diseases. We provide a new opinion that DT could be an effective solution for musculoskeletal system diseases in the future, which will help us analyze the real-time biomechanical properties of the musculoskeletal system and achieve personalized medicine.

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.

High frame rate deformation analysis of knee cartilage by spiral dualMRI and relaxation mapping

PURPOSE

Daily activities including walking impose high-frequency cyclic forces on cartilage and repetitive compressive deformation. Analyzing cartilage deformation during walking would provide spatial maps of displacement and strain and enable viscoelastic characterization, which may serve as imaging biomarkers for early cartilage degeneration when the damage is still reversible. However, the time-dependent biomechanics of cartilage is not well described, and how defects in the joint impact the viscoelastic response is unclear.

METHODS

We used spiral acquisition with displacement-encoding MRI to quantify displacement and strain maps at a high frame rate (25 frames/s) in tibiofemoral joints. We also employed relaxometry methods (T1 , T1ρ , T2 , T2 *) on the cartilage.

RESULTS

Normal and shear strains were concentrated on the bovine tibiofemoral contact area during loading, and the defected joint exhibited larger compressive strains. We also determined a positive correlation between the change of T1ρ in cartilage after cyclic loading and increased compressive strain on the defected joint. Viscoelastic behavior was quantified by the time-dependent displacement, where the damaged joint showed increased creep behavior compared to the intact joint. This technique was also successfully demonstrated on an in vivo human knee showing the gradual change of displacement during varus load.

CONCLUSION

Our results indicate that spiral scanning with displacement encoding can quantitatively differentiate the damaged from intact joint using the strain and creep response. The viscoelastic response identified with this methodology could serve as biomarkers to detect defects in joints in vivo and facilitate the early diagnosis of joint diseases such as osteoarthritis.

In vivo intervertebral disc mechanical deformation following a treadmill walking "stress test" is inversely related to T1rho relaxation time

OBJECTIVE

To assess the in vivo relationship between the mechanical response of intervertebral discs (IVDs) to dynamic activity and IVD biochemical composition assessed via T1rho relaxation imaging.

DESIGN

Eighteen asymptomatic participants with no history of low back pain (LBP), injury, or surgery underwent magnetic resonance (MR) imaging of their lumbar spine prior to and immediately following a treadmill walking "stress test." Anatomic (SPACE, FLASH) MR images were obtained pre- and post-exercise and utilized to measure IVD mechanical deformation. Quantitative (T1rho) imaging was performed pre-exercise to reflect IVD composition. Pre-exercise anatomic images were also utilized to assess IVD degenerative status based on the modified Pfirrmann scale. To quantify mechanical response, 3D surface models of the L1-L2-L5-S1 IVDs were created from manual segmentations of pre- and post-exercise anatomic images and utilized to assess changes in IVD height. IVD strain (%) was defined as change in IVD height normalized to pre-activity height. Linear mixed models were used to assess the relationships between IVD mechanical deformation (strain), composition (T1rho relaxation time), and degenerative status (Pfirrmann grade).

RESULTS

Increased compressive IVD strain was associated with lower T1rho relaxation times in the nucleus pulposus (NP) of the disc (βT1rho=5.07,CI:[1.52,7.77],Rmarg2=0.52,p=0.005). Thus, an inverse relationship between IVD strain and NP T1rho relaxation time was observed.

CONCLUSION

The in vivo mechanical response of the IVD to the "stress test" was sensitive to differences in NP composition. The results of this study suggest that quantification of in vivo IVD mechanical function and composition may provide insight into IVD health.

Ogden material calibration via magnetic resonance cartography, parameter sensitivity and variational system identification

Contemporary material characterization techniques that leverage deformation fields and the weak form of the equilibrium equations face challenges in the numerical solution procedure of the inverse characterization problem. As material models and descriptions differ, so too must the approaches for identifying parameters and their corresponding mechanisms. The widely used Ogden material model can be comprised of a chosen number of terms of the same mathematical form, which presents challenges of parsimonious representation, interpretability and stability. Robust techniques for system identification of any material model are important to assess and improve experimental design, in addition to their centrality to forward computations. Using fully three-dimensional displacement fields acquired in silicone elastomers with our recently developed magnetic resonance cartography (MR-u) technique on the order of greater than [Formula: see text], we leverage partial differential equation-constrained optimization as the basis of variational system identification of our material parameters. We incorporate the statistical F-test to maintain parsimony of representation. Using a new, local deformation decomposition locally into mixtures of biaxial and uniaxial tensile states, we evaluate experiments based on an analytical sensitivity metric and discuss the implications for experimental design. This article is part of the theme issue 'The Ogden model of rubber mechanics: Fifty years of impact on nonlinear elasticity'.

A BIOMECHANICAL STUDY OF THE EFFECTS OF FLEXION ANGLE ON THE INDUCTION MECHANISM OF CERVICAL SPONDYLOSIS

Lesions in facet joints such as bone hyperplasia and degenerative changes in the intervertebral discs, can compress nerve roots and the spinal cord, leading to cervical spondylosis (CS). Lesions in these parts of the spine are commonly related to abnormal loads caused by bad posture of the cervical spine. This study aimed to understand the potential mechanical effects of load amplitude on cervical spine motion to provide a theoretical basis for the biomechanical causes of CS, and to provide a reference for preventing of the condition. In this study, a finite element model of the normal human cervical spine (C1-C7) was established and validated using an infrared motion capture system to analyze the effects of flexion angle on the stresses experienced by intervertebral discs, the anterior edge of the vertebral body, the pedicle, uncinate and facet joints. Our analysis indicated that the intervertebral disc load increased by at least 70% during the 20∘ to 45∘ flexion of the neck with 121% load increase in the vertebrae. In the intervertebral discs, the stress was largest at C4-C5, and the stress was moderate at C5-C6. These results are consistent with clinical CS prone site research. According to Wolff’s law, when bones are placed under large stresses, hyperplasia can result to allow adaptation to large loads. Increased cervical spine flexion angles caused the proliferation of bone in the above-mentioned parts of the spine and can accelerate accelerating the appearance of CS.

Tissue-Engineering of Human Intervertebral Disc: A Concise Review

Intervertebral disc (IVD) represents a structure of crucial structural and functional importance for human spine. Pathology of IVD institutes a frequently encountered condition in current clinical practice. Degenerative Disc Disease (DDD), the principal clinical representative of IVD pathology, constitutes an increasingly diagnosed spinal disorder associated with substantial morbidity and mortality in recent years. Despite the considerable incidence and socioeconomic burden of DDD, existing treatment modalities including conservative and surgical methods have been demonstrated to provide a limited therapeutic effect, being not capable of interrupting or reversing natural progress of underlying disease. These limitations underline the requirement for development of novel, innovative and more effective therapeutic strategies for DDD management. Within this literature framework, compromised IVD replacement with a viable IVD construct manufactured with Tissue-Engineering (TE) methods has been recommended as a promising therapeutic strategy for DDD. Existing preliminary preclinical data demonstrate that proper combination of cells from various sources, different scaffold materials and appropriate signaling molecules renders manufacturing of whole-IVD tissue-engineered constructs a technically feasible process. Aim of this narrative review is to critically summarize current published evidence regarding particular aspects of IVD-TE, primarily emphasizing in providing researchers in this field with practicable knowledge in order to enhance clinical translatability of their research and informing clinical practitioners about the features and capabilities of innovative TE science in the field of IVD-TE.