Three-dimensional microstructural reconstruction of the ovine intervertebral disc using ultrahigh field MRI

Unveiling the mechanical role of radial fibers in meniscal tissue: Toward structural biomimetics

The meniscus tissue is crucial for knee joint biomechanics and is frequently susceptible to injuries resulting in early-onset osteoarthritis. Consequently, the need for meniscal substitutes spurs ongoing development. The meniscus is a composite tissue reinforced with circumferential and radial collagenous fibers; the mechanical role of the latter has yet to be fully unveiled. Here, we investigated the role of radial fibers using a synergistic methodology combining meniscal tissue structure imaging, a computational knee joint model, and the fabrication of simple biomimetic composite laminates. These laminates mimic the basic structural units of the meniscus, utilizing longitudinal and transverse fibers equivalent to the circumferential and radial fibers in meniscal tissue. In the computational model, the absence of radial fibers resulted in stress concentration within the meniscus matrix and up to 800 % greater area at the same stress level. Furthermore, the contact pressure on the tibial cartilage increased drastically, affecting up to 322 % larger areas. Conversely, in models with radial fibers, we observed up to 25 % lower peak contact pressures and width changes of less than 0.1 %. Correspondingly, biomimetic composite laminates containing transverse fibers exhibited minor transverse deformations and smaller Poisson's ratios. They demonstrated structural shielding ability, maintaining their mechanical performance with the reduced amount of fibers in the loading direction, similar to the ability of the torn meniscus to carry and transfer loads to some extent. These results indicate that radial fibers are essential to distribute contact pressure and tensile stresses and prevent excessive deformations, suggesting the importance of incorporating them in novel designs of meniscal substitutes. STATEMENT OF SIGNIFICANCE: The organization of the collagen fibers in the meniscus tissue is crucial to its biomechanical function. Radially oriented fibers are an important structural element of the meniscus and greatly affect its mechanical behavior. However, despite their importance to the meniscus mechanical function, radially oriented fibers receive minor attention in meniscal substitute designs. Here, we used a synergistic methodology that combines imaging of the meniscal tissue structure, a structural computational model of the knee joint, and the fabrication of simplistic biomimetic composite laminates that mimic the basic structural units of the meniscus. Our findings highlight the importance of the radially oriented fibers, their mechanical role in the meniscus tissue, and their importance as a crucial element in engineering novel meniscal substitutes.

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.

Structure-function characterization of the transition zone in the intervertebral disc

Understanding the structure-function relationship in the intervertebral disk (IVD) is crucial for the development of novel tissue engineering strategies to regenerate IVD and the establishment of accurate computational models for low back pain research. A large number of studies have improved our knowledge of the mechanical and structural properties of the nucleus pulposus (NP) and annulus fibrosus (AF), two of the main regions in the IVD. However, few studies have focused on the AF-NP interface (transition zone; TZ). Therefore, the current study aims to, for the first time, characterize the cyclic and failure mechanical properties of the TZ region under physiological loading (1, 3, and 5%s^-1 strain rates) and investigate the structural integration mechanisms between the NP, TZ, and AF regions. The results of the current study reveal significant effects of region (NP, TZ, and AF) and strain rates (1, 3, and 5%s^-1) on stiffness (p < 0.001). In addition, energy absorption is significantly higher for the AF compared to the TZ and NP (p <0.001) as well as between the TZ and NP (p <0.001). The current research finds adaptation, direct penetration, and entanglement between TZ and AF fibers as three common mechanisms for structural integration between the TZ and AF regions. STATEMENT OF SIGNIFICANCE: Despite a large number of studies that have mechanically, structurally, and biologically characterized the nucleus pulposus (NP) and annulus fibrosus (AF) regions, few studies have focused on the NP-AF interface region (known as Transition Zone; TZ) in the IVD; hence, our understanding of the TZ structure-function relationship is still incomplete. Of particular importance, the cyclic mechanical properties of the TZ, compared to the adjacent regions (NP and AF), are yet to be explored and the precise nature of the structural integration between the NP and AF via the TZ region is not yet known. The current study explores both the mechanical and structural properties of the TZ region to ultimately identify the mechanism of integration between the NP and AF.

An update on animal models of intervertebral disc degeneration and low back pain: Exploring the potential of artificial intelligence to improve research analysis and development of prospective therapeutics

Animal models have been invaluable in the identification of molecular events occurring in and contributing to intervertebral disc (IVD) degeneration and important therapeutic targets have been identified. Some outstanding animal models (murine, ovine, chondrodystrophoid canine) have been identified with their own strengths and weaknesses. The llama/alpaca, horse and kangaroo have emerged as new large species for IVD studies, and only time will tell if they will surpass the utility of existing models. The complexity of IVD degeneration poses difficulties in the selection of the most appropriate molecular target of many potential candidates, to focus on in the formulation of strategies to effect disc repair and regeneration. It may well be that many therapeutic objectives should be targeted simultaneously to effect a favorable outcome in human IVD degeneration. Use of animal models in isolation will not allow resolution of this complex issue and a paradigm shift and adoption of new methodologies is required to provide the next step forward in the determination of an effective repairative strategy for the IVD. AI has improved the accuracy and assessment of spinal imaging supporting clinical diagnostics and research efforts to better understand IVD degeneration and its treatment. Implementation of AI in the evaluation of histology data has improved the usefulness of a popular murine IVD model and could also be used in an ovine histopathological grading scheme that has been used to quantify degenerative IVD changes and stem cell mediated regeneration. These models are also attractive candidates for the evaluation of novel anti-oxidant compounds that counter inflammatory conditions in degenerate IVDs and promote IVD regeneration. Some of these compounds also have pain-relieving properties. AI has facilitated development of facial recognition pain assessment in animal IVD models offering the possibility of correlating the potential pain alleviating properties of some of these compounds with IVD regeneration.

Detailed mechanical characterization of the transition zone: New insight into the integration between the annulus and nucleus of the intervertebral disc

The Nucleus Pulposus (NP) and Annulus Fibrous (AF) are two primary regions of the intervertebral disc (IVD). The interface between the AF and NP, where the gradual transition in structure and type of fibers are observed, is known as the Transition Zone (TZ). Recent structural studies have shown that the TZ contains organized fibers that appear to connect the NP to the AF. However, the mechanical characteristics of the TZ are yet to be explored. The current study aimed to investigate the mechanical properties of the TZ at the anterolateral (AL) and posterolateral (PL) regions in both radial and circumferential directions of loading using ovine IVDs (N = 28). Young's and toe moduli, maximum stress, failure strain, strain at maximum stress, and toughness were calculated mechanical parameters. The findings from this study revealed that the mechanical properties of the TZ, including young's modulus (p = 0.001), failure strain (p < 0.001), strain at maximum stress (p = 0.002), toughness (p = 0.027), and toe modulus (p = 0.005), were significantly lower for the PL compared to the AL region. Maximum stress was not significantly different between the PL and AL regions (p = 0.164). We found that maximum stress (p = 0.002), failure strain (p < 0.001), and toughness (p = 0.001) were significantly different in different loading directions. No significant differences for modulus (young's; p = 0.169 and toe; p = 0.352) and strain at maximum stress (p = 0.727) were found between the radial and circumferential loading directions. STATEMENT OF SIGNIFICANCE: To date there has not been a study that has investigated the mechanical characterization of the annulus (AF)-nucleus (NP) interface (transition zone; TZ) in the intervertebral disc (IVD), nor is it known whether the posterolateral (PL) and anterolateral (AL) regions of the TZ exhibit different mechanical properties. Accordingly, the TZ mechanical properties have been rarely used in the development of computational IVD models and relevant tissue-engineered scaffolds. The current research reported the mechanical properties of the TZ region and revealed that its mechanical properties were significantly lower for the PL compared to the AL region. These new findings enhance our knowledge about the nature of AF-NP integration and may help to develop more realistic tissue-engineered or computational IVD models.

Regional variations in discrete collagen fibre mechanics within intact intervertebral disc resolved using synchrotron computed tomography and digital volume correlation

Many soft tissues, such as the intervertebral disc (IVD), have a hierarchical fibrous composite structure which suffers from regional damage. We hypothesise that these tissue regions have distinct, inherent fibre structure and structural response upon loading. Here we used synchrotron computed tomography (sCT) to resolve collagen fibre bundles (∼5μm width) in 3D throughout an intact native rat lumbar IVD under increasing compressive load. Using intact samples meant that tissue boundaries (such as endplate-disc or nucleus-annulus) and residual strain were preserved; this is vital for characterising both the inherent structure and structural changes upon loading in tissue regions functioning in a near-native environment. Nano-scale displacement measurements along >10,000 individual fibres were tracked, and fibre orientation, curvature and strain changes were compared between the posterior-lateral region and the anterior region. These methods can be widely applied to other soft tissues, to identify fibre structures which cause tissue regions to be more susceptible to injury and degeneration. Our results demonstrate for the first time that highly-localised changes in fibre orientation, curvature and strain indicate differences in regional strain transfer and mechanical function (e.g. tissue compliance). This included decreased fibre reorientation at higher loads, specific tissue morphology which reduced capacity for flexibility and high strain at the disc-endplate boundary.

Statement of significance

The analyses presented here are applicable to many collagenous soft tissues which suffer from regional damage. We aimed to investigate regional intervertebral disc (IVD) structural and functional differences by characterising collagen fibre architecture and linking specific fibre- and tissue-level deformation behaviours. Synchrotron CT provided the first demonstration of tracking discrete fibres in 3D within an intact IVD. Detailed analysis of regions was performed using over 200k points, spaced every 8 μm along 10k individual fibres. Such comprehensive structural characterisation is significant in informing future computational models. Morphological indicators of tissue compliance (change in fibre curvature and orientation) and fibre strain measurements revealed localised and regional differences in tissue behaviour.

Radial trend in murine annulus fibrosus fiber orientation is best explained by vertebral growth

PURPOSE

The intervertebral disc (IVD) annulus fibrosus (AF) is composed of concentric lamellae with alternating right- and left-handed helically oriented collagen fiber bundles. This arrangement results in anisotropic material properties, which depend on local fiber orientations. Prior measurements of fiber inclination angles in human lumbar and bovine caudal IVDs found a significantly higher inclination angle in the inner AF than outer, though it is currently unknown if this pattern is conserved in smaller mammalian species. Additionally, the physical mechanism behind this pattern remains un-determined.

METHODS

In this study, AF fiber angles were measured histologically in murine caudal IVDs and compared to previously published values from bovine caudal IVDs. Fiber angles were also predicted using three theoretical models, including two based on adaptation to internal swelling pressure and one based on vertebral body growth.

RESULTS

Fiber angle was found to significantly decrease from 49.5 ± 3.8° in the inner AF to 34.5 ± 6.6° in the outer AF. While steeper than in bovine discs at all locations, the trend with radial position was comparable between species. This trend was best fit by growth-based model and opposite of that predicted by the pressure vessel models.

CONCLUSION

Trends in AF fiber orientation are conserved between mammalian species. Modeling results suggest that the AF tissue microstructure is more likely to be driven by adjacent vertebral body growth than adapted for optimal mechanical performance.

The ultrastructural organization of elastic fibers at the interface of the nucleus and annulus of the intervertebral disk

There has been no study to describe the ultrastructural organization of elastic fibers at the interface of the nucleus pulposus and annulus fibrosus of the intervertebral disk (IVD), a region called the transition zone (TZ). A previously developed digestion technique was optimized to eliminate cells and non-elastin ECM components except for the elastic fibers from the anterolateral (AL) and posterolateral (PL) regions of the TZ in ovine IVDs. Not previously reported, the current study identified a complex elastic fiber network across the TZ for both AL and PL regions. In the AL region, this network consisted of major thick elastic fibers (≈ 1 µm) that were interconnected with delicate (< 200 nm) elastic fibers. While the same ultrastructural organization was observed in the PL region, interestingly the size of the elastic fibers was smaller (< 100 nm) compared to those that were located in the AL region. Quantitative analysis of the elastic fibers revealed significant differences in the size (p < 0.001) and the orientation of elastic fibers (p = 0.001) between the AL and PL regions, with a higher orientation and larger size of elastic fibers observed in the AL region. The gradual elimination of cells and non-elastin extracellular matrix components identified that elastic fibers in the TZ region in combination with the extracellular matrix created a honeycomb structure that was more compact at the AF interface compared to that located close to the NP. Three different symmetrically organized angles of rotation (0⁰ and ±90⁰) were detected for the honeycomb structure at both interfaces, and the structure was significantly orientated at the TZ-AF compared to the TZ-NP interface (p = 0.003).

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

Physical disruptions to intervertebral discs (IVDs) can cause mechanical changes that lead to degeneration and to low back pain which affects 75% of us in our lifetimes. Quantifying the effects of these changes on internal IVD strains may lead to better preventative strategies and treatments. Digital Volume Correlation (DVC) is a non-invasive technique that divides volumetric images into subsets, and measures strains by tracking the internal patterns within them under load. Applying DVC to MRIs may allow non-invasive strain measurements. However, DVC-MRI for strain measurements in IVDs has not been used previously. The purpose of this study was to quantify the strain and deformation errors associated with DVC-MRI for measurements in human IVDs. Eight human lumbar IVDs were MRI scanned (9.4 T) for a 'zero-strain study' (multiple unloaded scans to quantify noise within the system), and a loaded study (2 mm axial compression). Three DVC methodologies: Fast-Fourier transform (FFT), direct correlation (DC), and a combination of both FFT and DC approaches were compared with subset sizes ranging from 8 to 88 voxels to establish the optimal DVC methodology and settings which were then used in the loaded study. FFT + DC was the optimal method and a subset size of 56 voxels (2520 µm) was found to be a good compromise between errors and spatial resolution. Displacement and strain errors did not exceed 28 µm and 3000 microstrain, respectively. These findings demonstrate that DVC-MRI can quantify internal strains within IVDs non-invasively and accurately. The method has unique potential for assessing IVD strains within patients.

Synchrotron tomography of intervertebral disc deformation quantified by digital volume correlation reveals microstructural influence on strain patterns

The intervertebral disc (IVD) has a complex and multiscale extracellular matrix structure which provides unique mechanical properties to withstand physiological loading. Low back pain has been linked to degeneration of the disc but reparative treatments are not currently available. Characterising the disc's 3D microstructure and its response in a physiologically relevant loading environment is required to improve understanding of degeneration and to develop new reparative treatments. In this study, techniques for imaging the native IVD, measuring internal deformation and mapping volumetric strain were applied to an in situ compressed ex vivo rat lumbar spine segment. Synchrotron X-ray micro-tomography (synchrotron CT) was used to resolve IVD structures at microscale resolution. These image data enabled 3D quantification of collagen bundle orientation and measurement of local displacement in the annulus fibrosus between sequential scans using digital volume correlation (DVC). The volumetric strain mapped from synchrotron CT provided a detailed insight into the micromechanics of native IVD tissue. The DVC findings showed that there was no slipping at lamella boundaries, and local strain patterns were of a similar distribution to the previously reported elastic network with some heterogeneous areas and maximum strain direction aligned with bundle orientation, suggesting bundle stretching and sliding. This method has the potential to bridge the gap between measures of macro-mechanical properties and the local 3D micro-mechanical environment experienced by cells. This is the first evaluation of strain at the micro scale level in the intact IVD and provides a quantitative framework for future IVD degeneration mechanics studies and testing of tissue engineered IVD replacements. STATEMENT OF SIGNIFICANCE: Synchrotron in-line phase contrast X-ray tomography provided the first visualisation of native intact intervertebral disc microstructural deformation in 3D. For two annulus fibrosus volumes of interest, collagen bundle orientation was quantified and local displacement mapped as strain. Direct evidence of microstructural influence on strain patterns could be seen such as no slipping at lamellae boundaries and maximum strain direction aligned with collagen bundle orientation. Although disc elastic structures were not directly observed, the strain patterns had a similar distribution to the previously reported elastic network. This study presents technical advances and is a basis for future X-ray microscopy, structural quantification and digital volume correlation strain analysis of soft tissue.

Microstructural characterization of annulus fibrosus by ultrasonography: a feasibility study with an in vivo and in vitro approach

The main function of the intervertebral disc is biomechanical function, since it must resist repetitive high loadings, while giving the spine its flexibility and protecting the spinal cord from over-straining. It partially owes its mechanical characteristics to the lamellar architecture of its outer layer, the annulus fibrosus. Today, no non-invasive means exist to characterize annulus lamellar structure in vivo. The aim of this work was to test the feasibility of imaging annulus fibrosus microstructure in vivo with ultrasonography. Twenty-nine healthy adolescents were included. Ultrasonographies of L3-L4 disc were acquired with a frontal approach. Annulus fibrosus was segmented in the images to measure the thickness of the lamellae. To validate lamellar appearance in ultrasonographies, multimodality images of two cow tail discs were compared: ultrasonography, magnetic resonance and optical microscopy. In vivo average lamellar thickness was 229.7 ± 91.5 μm, and it correlated with patient body mass index and age. Lamellar appearance in the three imaging modalities in vitro was consistent. Lamellar measurement uncertainty was 7%, with good agreement between two operators. Feasibility of ultrasonography for the analysis of lumbar annulus fibrosus structure was confirmed. Further work should aim at validating measurement reliability, and to assess the relevance of the method to characterize annulus alterations, for instance in disc degeneration or scoliosis.

The Mechanical Role of the Radial Fiber Network Within the Annulus Fibrosus of the Lumbar Intervertebral Disc: A Finite Elements Study

The annulus fibrosus (AF) of the intervertebral disc (IVD) consists of a set of concentric layers composed of a primary circumferential collagen fibers arranged in an alternating oblique orientation. Moreover, there exists an additional secondary set of radial translamellar collagen fibers which connects the concentric layers, creating an interconnected fiber network. The aim of this study was to investigate the mechanical role of the radial fiber network. Toward that goal, a three-dimensional (3D) finite element model of the L3–L4 spinal segment was generated and calibrated to axial compression and pure moment loading. The AF model explicitly recognizes the two heterogeneous networks of fibers. The presence of radial fibers demonstrated a pronounced effect on the local disc responses under lateral bending, flexion, and extension modes. In these modes, the radial fibers were in a tensile state in the disc region that subjected to compression. In addition, the circumferential fibers, on the opposite side of the IVD, were also under tension. The local stress in the matrix was decreased in up to 9% in the radial fibers presence. This implies an active fiber network acting collectively to reduce the stresses and strains in the AF lamellae. Moreover, a reduction of 26.6% in the matrix sideways expansion was seen in the presence of the radial fibers near the neutral bending axis of the disc. The proposed biomechanical model provided a new insight into the mechanical role of the radial collagen fibers in the AF structure. This model can assist in the design of future IVD substitutes.