Region- and degeneration dependent stiffness distribution in intervertebral discs derived by shear wave elastography

Information on the local stiffness characteristics of the intervertebral disc (IVD) is crucial for the understanding of its structure-function properties in health and disease and may improve numerical modeling. Previous studies have attempted to map local tissue stiffness by sectioning the disc and performing mechanical testing on these discrete tissue units, which is technically challenging and may bias the results. Shear wave elastography (SWE) represents a nondestructive alternative that can provide spatially continuous elasticity estimates. We investigated the feasibility of SWE for human intervertebral disc elasticity mapping in a laboratory setting. To this end, global spinal segment mechanical behavior was determined in 6 loading directions and served as ground truth data for the validation of the approach. Subsequently, the cranial spinal vertebra was removed and shear wave elastographic scans of the IVD were acquired. SWE-measurements were reconstructed into three-dimensional elastographic maps, discretized into distinct IVD regions and correlated with global segment mechanical parameters. SWE-derived Young's modulus estimates were compared among different regions and as a function of their state of degeneration. We found annulus shear wave speed to be moderately correlated with segment mechanical behavior irrespective of the loading direction whereas shear wave speed in the nucleus pulposus showed a very weak association (mean (SD) absolute Pearson correlation coefficients: 0.51 (0.14) and 0.17 (0.12), respectively). Young's modulus mapping of the intervertebral disc revealed stiffness to be highest in the ventral annulus with a stiffness decrease both circumferentially towards the dorsal aspect as well as towards the center of the disc. SWE hence provides a valid alternative to disc sectioning and piecewise mechanical testing.