New findings confirm the viscoelastic behaviour of the inter-lamellar matrix of the disc annulus fibrosus in radial and circumferential directions of loading

[Biomechanics of the intervertebral disc : Consequences of degenerative changes]

The intervertebral disc represents the flexible connection between two adjacent vertebral bodies. Intervertebral discs, therefore, give the spine its enormous range of motion. At the same time, intervertebral discs distribute the load evenly over the bony vertebral bodies to ensure load transfer from the upper body to the pelvis, provide sufficient stability, and absorb shocks during everyday movements as well as under extreme loads. The two tissue regions of the intervertebral disc, the central fluid-rich core (nucleus pulposus) and the lamellae of the outer fibrous ring (annulus fibrosus) experience different stresses under load and movement. This article summarizes current knowledge on the biomechanical properties of the intervertebral disc and explains how these are altered by degeneration and which surgical treatment options for degenerated and herniated discs are advisable from a biomechanical perspective.

Stress relaxation behavior of the transition zone in the intervertebral disc

The stress relaxation of the TZ region, located at the interface of the Annulus Fibrosus (AF) and Nucleus Pulposus (NP) of the disc, and how its stress is relaxed compared to the adjacent regions is unknown. The current study aimed to identify the TZ stress relaxation properties under different strain magnitudes (0.2, 0.4, and 0.6 mm/mm) and compared the TZ stress relaxation characteristics to the NP and inner AF (IAF) regions at a specific strain magnitude (0.6 mm/mm). The results of the current study revealed that the TZ region exhibited different stress relaxation properties under various strain magnitudes with significantly higher initial (p < 0.008) and reduced stresses (marginally; p = 0.06) at higher strains. Our experimental stress relaxation data revealed a significantly higher equilibrium stress for the IAF compared to the TZ and NP regions (p < 0.001) but not between the TZ and NP regions (p = 0.7). We found that NP radial stress relaxed significantly faster (p < 0.04) than the TZ and NP. Additionally, the current study proposed a simple mathematical model and identified that, consistent with experimental data, the overall effect of region on both the level of decayed stress and the rate at which stress is relaxed was significant (p < 0.006). The current study found a similar stress relaxation characteristic between the NP and TZ regions, while IAF exhibited different stress relaxation properties. It is possible that this mismatch in stress relaxation acts as a shape transformation mechanism triggered by viscoelastic behavior. STATEMENT OF SIGNIFICANCE: Our understanding of the biomechanical properties of the transition zone (TZ) in the IVD, a region at the interface of the Nucleus Pulposus (NP) and Annulus Fibrosus (AF), is sparse. Unfortunately, there are no current studies that investigate the TZ stress relaxation properties and how stress is relaxed in the TZ compared to the adjacent regions. For the first time, the current study characterized the stress relaxation properties of the TZ and described how the TZ stress is relaxed compared to its adjacent regions.

Mechanical consequences to the annulus fibrosus following rapid internal pressurization and endplate fracture under restrained-expansion conditions

Intervertebral disc herniation is not a common injury in the adolescent population, but the correlation between trauma and herniation warrants concern. Previous research demonstrated the capacity for rapid internal pressurization to reduce the mechanical integrity of the intervertebral disc's annulus fibrosus, even in the absence of fracture. The purpose of this study was to modify previous internal pressurization procedures towards a more transferable injury model, then investigate the capacity for these procedures to damage the mechanical integrity of the annulus fibrosus. Porcine cervical motion segments with intact facet joints were confined between a vice and force plate under 300 N of static compression, then a single, manual, rapid internal pressurization was delivered. Posterolateral annulus samples were extracted and situated in a 180° peel test configuration, exposing the interlamellar matrix of samples to separations of 0.5 mm/s, until complete separation of the sample occurred. Multilayer tensile testing was performed on superficial and mid-span samples of annulus by applying uniaxial tension of 1 %/s to 50 % strain. Compared to unpressurized controls, rapid pressurization causing fracture resulted in reduced lamellar adhesion and increased toe-region stress and strain properties in the annulus. Morphological assessment reported similar fracture patterns between endplate fractures achieved in the present experiment and endplate fractures documented in human patients. Mechanical plus morphological results suggest that rapid internal pressurization resulting in endplate fracture may represent a potent mechanism for subsequent damage to the intervertebral disc.

Comprehensive modeling of annulus fibrosus: From biphasic refined characterization to damage accumulation under viscous loading

The annulus fibrosus (AF), a permeable, hydrated, and fiber-reinforced soft tissue, exhibits complex responses influenced by fluid pressure, osmotic pressure, and structural mechanics. Existing models struggle to comprehensively represent these intricate interactions and the heterogeneous solid responses within the AF. Additionally, the mechanisms driving differential damage accumulation between non-degenerative and degenerative intervertebral discs remain poorly understood. In this study, we introduce a biphasic-swelling damage model for the AF. We conceptually develop and rigorously validate this model through tissue-level tests employing various loading modes, consistently aligning model predictions with experimental data. Leveraging parametric geometric algorithms and custom Python scripts, we construct models simulating both non-degenerative and degenerative discs. Following calibration, we subject these models to viscous loading protocols. Our findings reveal the posterior AF's susceptibility to damage, contingent upon loading rate and water content. We elucidate the underlying mechanisms by examining the temporal evolution of fluid pressure, osmotic pressure, and the regionally dependent fiber network. This research presents a highly accurate model of the AF, providing valuable insights into disc damage. Future research endeavors should expand this model to incorporate ionic transport and diffusion, enabling a more profound exploration of intervertebral disc mechanobiology. This comprehensive model contributes to a better understanding of AF behavior and may inform therapeutic strategies for disc-related pathologies. STATEMENT OF SIGNIFICANCE: This research presents a comprehensive model of the annulus fibrosus (AF), a crucial component of the intervertebral disc that provides structural support and resists deformation. The study introduces a biphasic-swelling damage model for the AF and validates it through tissue-level tests. The model accounts for fluid pressure, osmotic pressure, and matrix mechanics, providing a more accurate representation of the AF's behavior. The study also investigates the differential damage accumulation between non-degenerative and degenerative discs, shedding light on the mechanisms driving disc degeneration. The findings have significant implications for medical treatments and interventions, as they highlight the posterior AF's susceptibility to damage. This research is of great interest to readers interested in biomechanics, tissue engineering, and medical treatments for disc degeneration.

Biomechanical response of lumbar intervertebral disc in daily sitting postures: a poroelastic finite element analysis

This study aims to establish and validate a poroelastic L4-L5 finite element model to evaluate the effect of different sitting postures and their durations on the mechanical responses of the disc. During the sustained loading conditions, the height loss, fluid loss and von-Mises stress gradually increased, but the intradiscal pressure decreased. The varying rates of aforementioned parameters were more significant at the initial loading stage and less so at the end. The predicted values in the flexed sitting posture were significantly greater than other postures. The extended sitting posture caused an obvious von-Mises stress concentration in the posterior region of the inter-lamellar matrix. From the biomechanical perspective, prolonged sitting may pose a high risk of lumbar disc degeneration, and therefore adjusting the posture properly in the early stage of sitting time may be useful to mitigate that. Additionally, upright sitting is a safer posture, while flexed sitting posture is more harmful.

On the identification of the ultra-structural organization of elastic fibers and their effects on the integrity of annulus fibrosus

Due to the complicated structure of the elastic fiber network in annulus fibrosus, existing in-silico studies either simplified or just overlooked its distribution pattern. Nonetheless, experimental and simulation results have proven that elastic fibers are of great importance to maintaining the structural integrity of annulus fibrosus and therefore to ensuring the load-bearing ability of intervertebral discs. Such needs call for a fine model. This work aims at developing a biphasic annulus fibrosus model by incorporating the accurate distribution pattern of collagen and elastic fibers. Both the structural parameters and intrinsic mechanical parameters were successfully identified using single lamella and inter-lamella microscopy anatomy and micromechanical testing data. The proposed model was then used to implement finite element simulations on various anterior and posterolateral multi-lamellae annulus fibrosus specimens. In general, simulation results agree well with available experimental and simulation data. On this basis, the effects of elastic fibers on the integrity of annulus fibrosus were further investigated. It was found that elastic fibers significantly influence the free swelling, radial stretching and circumferential shear performances of annulus fibrosus. Nonetheless, no significant effects were found for the circumferential stretching capability. The proposed biphasic model considers for the first time the distribution characteristics of elastic fibers at two scales, including both the principal orientations of all fiber families and the detailed distribution pattern within each family. Better understandings on the functions of collagen and elastic fibers can therefore be realized. To further enhance its prediction capability, the current model can be extended in the future by taking the fiber-matrix interaction as well as progressive damages into consideration.

Emerging tissue engineering strategies for annulus fibrosus therapy

Low back pain is a major public health concern experienced by 80% of the world's population during their lifetime, which is closely associated with intervertebral disc (IVD) herniation. IVD herniation manifests as the nucleus pulposus (NP) protruding beyond the boundaries of the intervertebral disc due to disruption of the annulus fibrosus (AF). With a deepening understanding of the importance of the AF structure in the pathogenesis of intervertebral disc degeneration, numerous advanced therapeutic strategies for AF based on tissue engineering, cellular regeneration, and gene therapy have emerged. However, there is still no consensus concerning the optimal approach for AF regeneration. In this review, we summarized strategies in the field of AF repair and highlighted ideal cell types and pro-differentiation targeting approaches for AF repair, and discussed the prospects and difficulties of implant systems combining cells and biomaterials to guide future research directions. STATEMENT OF SIGNIFICANCE: Low back pain is a major public health concern experienced by 80% of the world's population during their lifetime, which is closely associated with intervertebral disc (IVD) herniation. However, there is still no consensus concerning the optimal approach for annulus fibrosus (AF) regeneration. In this review, we summarized strategies in the field of AF repair and highlighted ideal cell types and pro-differentiation targeting approaches for AF repair, and discussed the prospects and difficulties of implant systems combining cells and biomaterials to guide future research directions.

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.

A Microstructure-Based Mechanistic Approach to Detect Degeneration Effects on Potential Damage Zones and Morphology of Young and Old Human Intervertebral Discs

There is an increasing demand to develop predictive medicine through the creation of predictive models and digital twins of the different body organs. To obtain accurate predictions, real local microstructure, morphology changes and their accompanying physiological degenerative effects must be taken into account. In this article, we present a numerical model to estimate the long-term aging effect on the human intervertebral disc response by means of a microstructure-based mechanistic approach. It allows to monitor in-silico the variations in disc geometry and local mechanical fields induced by age-dependent long-term microstructure changes. Both lamellar and interlamellar zones of the disc annulus fibrosus are constitutively represented by considering the main underlying microstructure features in terms of proteoglycans network viscoelasticity, collagen network elasticity (along with content and orientation) and chemical-induced fluid transfer. With age, a noticeable increase in shear strain is especially observed in the posterior and lateral posterior regions of the annulus which is in correlation with the high vulnerability of elderly people to back problems and posterior disc hernia. Important insights about the relation between age-dependent microstructure features, disc mechanics and disc damage are revealed using the present approach. These numerical observations are hardly obtainable using current experimental technologies which makes our numerical tool useful for patient-specific long-term predictions.

Elastic Fibers in the Intervertebral Disc: From Form to Function and toward Regeneration

Despite extensive efforts over the past 40 years, there is still a significant gap in knowledge of the characteristics of elastic fibers in the intervertebral disc (IVD). More studies are required to clarify the potential contribution of elastic fibers to the IVD (healthy and diseased) function and recommend critical areas for future investigations. On the other hand, current IVD in-vitro models are not true reflections of the complex biological IVD tissue and the role of elastic fibers has often been ignored in developing relevant tissue-engineered scaffolds and realistic computational models. This has affected the progress of IVD studies (tissue engineering solutions, biomechanics, fundamental biology) and translation into clinical practice. Motivated by the current gap, the current review paper presents a comprehensive study (from the early 1980s to 2022) that explores the current understanding of structural (multi-scale hierarchy), biological (development and aging, elastin content, and cell-fiber interaction), and biomechanical properties of the IVD elastic fibers, and provides new insights into future investigations in this domain.

Delamination of the Annulus Fibrosus of the Intervertebral Disc: Using a Bovine Tail Model to Examine Effect of Separation Rate

The intervertebral disc (IVD) is a complex structure, and recent evidence suggests that separations or delamination between layers of the annulus may contribute to degeneration development, a common cause of low back pain The purpose of the present experiment was to quantify the mechanical response of the layer-adjoining interlamellar matrix at different rates of separation. Understanding the rate-dependency of the interlamellar matrix, or the adhesion between adjacent layers of the disc, is important as the spine experiences various loading velocities during activities of daily living. Twelve discs were dissected from four bovine tails (three extracts per tail). Two multi-layered annulus samples were collected from each IVD (total = 24, mean bond width = 3.82 ± 0.96 mm) and randomly assigned to a 180° peel test at one of three delamination rates; 0.05 mm/s, 0.5 mm/s, or 5 mm/s. Annulus extracts were found to have similar maximal adhesion strengths (p = 0.39) and stiffness (p = 0.97) across all rate conditions. However, a significant difference in lamellar adhesion strength variability was observed between the 5 mm/s condition (0.96 N/mm ± 0.31) when compared to the 0.5 mm/s (0.50 N/mm ± 0.19) and 0.05 mm/s (0.37 N/mm ± 0.13) conditions (p < 0.05). Increased variability may be indicative of non-uniform strength due to inconsistent adhesion throughout the interlamellar matrix, which is exacerbated by increased rates of loading. The observed non-uniform strength could possibly lead to a scenario more favourable to the development of microtrauma, and eventual delamination.

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.

Modeling multiaxial damage regional variation in human annulus fibrosus

In the present article, a fully three-dimensional human annulus fibrosus model is developed by considering the regional variation of the complex structural organization of collagen network at different scales to predict the regional anisotropic multiaxial damage of the intervertebral disc. The model parameters are identified using experimental data considering as elementary structural unit, the single annulus lamellae stretched till failure along the micro-sized collagen fibers. The multi-layered lamellar/inter-lamellar annulus model is constructed by considering the effective interactions between adjacent layers and the chemical-induced volumetric strain. The regional dependent model predictions are analyzed under various loading modes and compared to experimental data when available. The stretching along the circumferential and radial directions till failure serves to check the predictive capacities of the annulus model. Model results under simple shear, biaxial stretching and plane-strain compression are further presented and discussed. Finally, a full disc model is constructed using the regional annulus model and simulations are presented to assess the most likely failed areas under disc axial compression. STATEMENT OF SIGNIFICANCE: The damage in annulus soft tissues is a complex multiscale phenomenon due to a complex structural arrangement of collagen network at different scales of hierarchical organization. A fully three-dimensional constitutive representation that considers the regional variation of the structural complexity to estimate annulus multiaxial mechanics till failure has not yet been developed. Here, a model is developed to predict deformation-induced damage and failure of annulus under multiaxial loading histories considering as time-dependent physical process both chemical-induced volumetric effects and damage accumulation. After model identification using single lamellae extracted from different disc regions, the model predictability is verified for various multiaxial elementary loading modes representative of the spine movement. The heterogeneous mechanics of a full human disc model is finally presented.

Angle-ply scaffold supports annulus fibrosus matrix expression and remodeling by mesenchymal stromal and annulus fibrosus cells

The angle-ply multilaminate structure of the annulus fibrosus is not reestablished following discectomy which leads to reherniation of the intervertebral disc (IVD). Biomimetic scaffolds developed to repair these defects should be evaluated for their ability to support tissue regeneration by endogenous and exogenous cells. Herein a collagen-based, angle-ply multilaminate patch designed to repair the outer annulus fibrosus was assessed for its ability to support mesenchymal stromal and annulus fibrosus cell viability, elongation, alignment, extracellular matrix gene expression, and scaffold remodeling. Results demonstrated that the cells remained viable, elongated, and aligned along the collagen fiber preferred direction of the scaffold, upregulated genes associated with annulus fibrosus matrix and produced collagen on the scaffold yielding biaxial mechanical properties that resembled native annulus fibrosus tissue. In conclusion, these scaffolds have demonstrated their potential to promote a living repair of defects in the annulus fibrosus and thus may be used to prevent recurrent IVD herniations.

A microstructure-based model for a full lamellar-interlamellar displacement and shear strain mapping inside human intervertebral disc core

The determinant role of the annulus fibrosus interlamellar zones in the intervertebral disc transversal and volumetric responses and hence on their corresponding three-dimensional conducts have been only revealed and appreciated recently. Their consideration in disc modeling strategies has been proven to be essential for the reproduction of correct local strain and displacement fields inside the disc especially in the unconstrained directions of the disc. In addition, these zones are known to be the starting areas of annulus fibrosus circumferential tears and disc delamination failure mode, which is often judged as one of the most dangerous disc failure modes that could evolve with time leading to disc hernia. For this latter reason, the main goal of the current contribution is to incorporate physically for the first time, the interlamellar zones, at the scale of a complete human lumbar intervertebral disc, in order to allow a correct local vision and replication of the different lamellar-interlamellar interactions and an identification of the interlamellar critical zones. By means of a fully tridimensional chemo-viscoelastic constitutive model, which we implemented into a finite element code, the physical, mechanical and chemical contribution of the interlamellar zones is added to the disc. The chemical-induced volumetric response is accounted by the model for both the interlamellar zones and the lamellae using experimentally-based fluid kinetics. Computational simulations are performed and critically discussed upon different simple and complex physiological movements. The disc core and the interlamellar zones are numerically accessed, allowing the observation of the displacement and shear strain fields that are compared to direct MRI experiments from the literature. Important conclusions about the correct lamellar-interlamellar-nucleus interactions are provided thanks to the developed model. The critical interlamellar spots with the highest delamination potentials are defined, analyzed and related to the local kinetics and microstructure.

The effect of compressive loading rate on annulus fibrosus strength following endplate fracture

Intervertebral disc degeneration poses a considerable healthcare challenge, although the process is not well understood. Endplate fracture marks severe biomechanical compromise in a segment and may be correlated with degeneration of the disc. The purpose of this experiment was to investigate the relationship between endplate fracture velocity and damage to the annulus fibrosus. Following overnight-thawing, 27 frozen porcine cervical spines were dissected into motion segments (vertebra-disc-vertebra) and compressed until fracture at one of three loading rates (fast=15 mm/s, medium=1.5 mm/s, and slow=0.15 mm/s), or remained unfractured (control). Two annular samples were extracted and mechanically tested from each segment: 1) Bilayer samples underwent uniaxial tension to a stretch-ratio of 1.5; 2) Multilayer samples were delaminated with a 180° peel test configuration. All three rates of compression resulted in specimen fracture observed in the endplate and/or vertebra with varying degree of severity. Significant differences were detected in compressive strength and stiffness of motion segments when loaded at different rates of compression; interestingly these differences were not observed in the mechanical properties of the annulus fibrosus suggesting that at slow rates of loading, fracture of the endplate precedes destruction of the annulus fibrosus. In corroboration of these findings, gross and histological analysis reported no signs of annular disruption, strengthening assertions that annular damage did not occur.

Modeling of human intervertebral disc annulus fibrosus with complex multi-fiber networks

Collagen fibers within the annulus fibrosus (AF) lamellae are unidirectionally aligned with alternating orientations between adjacent layers. AF constitutive models often combine two adjacent lamellae into a single equivalent layer containing two fiber networks with a crisscross pattern. Additionally, AF models overlook the inter-lamellar matrix (ILM) as well as elastic fiber networks in between lamellae. We developed a nonhomogenous micromechanical model as well as two coarser homogenous hyperelastic and microplane models of the human AF, and compared their performances against measurements (tissue level uniaxial and biaxial tests as well as whole disc experiments) and seven published hyperelastic models. The micromechanical model had a realistic non-homogenous distribution of collagen fiber networks within each lamella and elastic fiber network in the ILM. For small matrix linear moduli (<0.2 MPa), the ILM showed substantial anisotropy (>10%) due to the elastic fiber network. However, at moduli >0.2 MPa, the effects of the elastic fiber network on differences in stress-strain responses at different directions disappeared (<10%). Variations in sample geometry and boundary conditions (due to uncertainty) markedly affected stress-strain responses of the tissue in uniaxial and biaxial tests (up to 16 times). In tissue level tests, therefore, simulations should represent testing conditions (e.g., boundary conditions, specimen geometry, preloads) as closely as possible. Stress/strain fields estimated from the single equivalent layer approach (conventional method) yielded different results from those predicted by the anatomically more accurate apparoach (i.e., layerwise). In addition, in a disc under a compressive force (symmetric loading), asymmetric stress-strain distributions were computed when using a layerwise simulation. Although all developed and selected published AF models predicted gross compression-displacement responses of the whole disc within the range of measured data, some showed excessively stiff or compliant responses under tissue-level uniaxial/biaxial tests. This study emphasizes, when constructing and validating constitutive models of AF, the importance of the proper simulation of individual lamellae as distinct layers, and testing parameters (sample geometric dimensions/loading/boundary conditions).

Interlamellar matrix governs human annulus fibrosus multiaxial behavior

Establishing accurate structure–property relationships for intervertebral disc annulus fibrosus tissue is a fundamental task for a reliable computer simulation of the human spine but needs excessive theoretical-numerical-experimental works. The difficulty emanates from multiaxiality and anisotropy of the tissue response along with regional dependency of a complex hierarchic structure interacting with the surrounding environment. We present a new and simple hybrid microstructure-based experimental/modeling strategy allowing adaptation of animal disc model to human one. The trans-species strategy requires solely the basic knowledge of the uniaxial circumferential response of two different animal disc regions to predict the multiaxial response of any human disc region. This work demonstrates for the first time the determining role of the interlamellar matrix connecting the fibers-reinforced lamellae in the disc multiaxial response. Our approach shows encouraging multiaxial predictive capabilities making it a promising tool for human spine long-term prediction.

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).

Elastic fibers: The missing key to improve engineering concepts for reconstruction of the Nucleus Pulposus in the intervertebral disc

The increasing prevalence of low back pain has imposed a heavy economic burden on global healthcare systems. Intense research activities have been performed for the regeneration of the Nucleus Pulposus (NP) of the IVD; however, tissue-engineered scaffolds have failed to capture the multi-scale structural hierarchy of the native tissue. The current study revealed for the first time, that elastic fibers form a network across the NP consisting of straight and thick parallel fibers that were interconnected by wavy fine fibers and strands. Both straight fibers and twisted strands were regularly merged or branched to form a fine elastic network across the NP. As a key structural feature, ultrathin (53 ± 7 nm), thin (215 ± 20 nm), and thick (890 ± 12 nm) elastic fibers were observed in the NP. While our quantitative analysis for measurement of the thickness of elastic fibers revealed no significant differences (p < 0.633), the preferential orientation of fibers was found to be significantly different (p < 0.001) across the NP. The distribution of orientation for the elastic fibers in the NP represented one major organized angle of orientation except for the central NP. We found that the distribution of elastic fibers in the central NP was different from those located in the peripheral regions representing two symmetrically organized major peaks (±45⁰). No significant differences in the maximum fiber count at the major angles of orientation (±45⁰) were observed for both peripheral (p = 0.427) and central NP (p = 0.788). Based on these new findings a structural model for the elastic fibers in the NP was proposed. The geometrical presentation, along with the distribution of elastic fibers orientation, resulting from the present study identifies the ultrastructural organization of elastic fibers in the NP important towards understanding their mechanical role which is still under investigation. Given the results of this new geometrical analysis, more-accurate multiscale finite element models can now be developed, which will provide new insights into the mechanobiology of the IVD. In addition, the results of this study can potentially be used for the fabrication of bio-inspired tissue-engineered scaffolds and IVD models to truly capture the multi-scale structural hierarchy of IVDs. STATEMENT OF SIGNIFICANCE: Visualization of elastic fibers in the nucleus of the intervertebral disk under high magnification was not reported before. The present research utilized extracellular matrix partial digestion to address significant gaps in understanding of nucleus microstructure that can potentially be used for the fabrication of bio-inspired tissue-engineered scaffolds and disk models to truly capture the multi-scale structural hierarchy of discs.

Advanced Strategies for the Regeneration of Lumbar Disc Annulus Fibrosus

Damage to the annulus fibrosus (AF), the outer region of the intervertebral disc (IVD), results in an undesirable condition that may accelerate IVD degeneration causing low back pain. Despite intense research interest, attempts to regenerate the IVD have failed so far and no effective strategy has translated into a successful clinical outcome. Of particular significance, the failure of strategies to repair the AF has been a major drawback in the regeneration of IVD and nucleus replacement. It is unlikely to secure regenerative mediators (cells, genes, and biomolecules) and artificial nucleus materials after injection with an unsealed AF, as IVD is exposed to significant load and large deformation during daily activities. The AF defects strongly change the mechanical properties of the IVD and activate catabolic routes that are responsible for accelerating IVD degeneration. Therefore, there is a strong need to develop effective therapeutic strategies to prevent or reconstruct AF damage to support operational IVD regenerative strategies and nucleus replacement. By the way of this review, repair and regenerative strategies for AF reconstruction, their current status, challenges ahead, and future outlooks were discussed.

Osmo-inelastic response of the intervertebral disc annulus fibrosus tissue

The aim of this article is to provide some insights on the osmo-inelastic response under stretching of annulus fibrosus of the intervertebral disc. Circumferentially oriented specimens of square cross section, extracted from different regions of bovine cervical discs (ventral-lateral and dorsal-lateral), are tested under different strain-rates and saline concentrations within normal range of strains. An accurate optical strain measuring technique, based upon digital image correlation, is used in order to determine the full-field displacements in the lamellae and fibers planes of the layered soft tissue. Annulus stress-stretch relationships are measured along with full-field transversal strains in the two planes. The mechanical response is found hysteretic, rate-dependent and osmolarity-dependent with a Poisson's ratio higher than 0.5 in the fibers plane and negative (auxeticity) in the lamellae plane. While the stiffness presents a regional-dependency due to variations in collagen fibers content/orientation, the strain-rate sensitivity of the response is found independent on the region. A significant osmotic effect is found on both the auxetic response in the lamellae plane and the stiffness rate-sensitivity. These local experimental observations will result in more accurate chemo-mechanical modeling of the disc annulus and a clearer multi-scale understanding of the disc intervertebral function.

How pre-strain affects the chemo-torsional response of the intervertebral disc

BACKGROUND

The role of the axial pre-strain on the torsional response of the intervertebral disc remains largely undefined. Moreover, the chemo-mechanical interactions in disc tissues are still unclear and corresponding data are rare in the literature. The paper deals with an in-vitro study of the pre-strain effect on the chemical sensitivity of the disc torsional response.

METHODS

Fifteen non-frozen 'motion segments' (two vertebrae and the intervening soft tissues) were extracted from the cervical spines of mature sheep. The motion segments were loaded in torsion at various saline concentrations and axial pre-strain levels in order to modulate the intradiscal pressure. After preconditioning with successive low-strain compressions at a magnitude of 0.1 mm (10 cycles at 0.05 mm/s), the motion segment was subjected to a cyclic torsion until a twisting level of 2 deg. at 0.05 deg./s while a constant axial pre-strain (in compression or in tension) is maintained, the saline concentration of the surrounding fluid bath being changed from hypo-osmotic condition to hyper-osmotic condition.

FINDINGS

Analysis of variance shows that the saline concentration influences the torsional response only when the motion segments are pre-compressed (p < .001) with significant differences between hypo-osmotic condition and hyper-osmotic condition.

INTERPRETATION

The combination of a compressive pre-strain with twisting amplifies the nucleus hydrostatic pressure on the annulus and the annulus collagen fibers tensions. The proteoglycans density increases with the compressive pre-strain and leads to higher chemical imbalances, which would explain the increase in chemical sensitivity of the disc torsional response.

Interlamellar-induced time-dependent response of intervertebral disc annulus: A microstructure-based chemo-viscoelastic model

The annulus fibrosus of the intervertebral disc exhibits an unusual transversal behavior for which a constitutive representation that considers as well regional effect, chemical sensitivity and time-dependency has not yet been developed, and it is hence the aim of the present contribution. A physically-based model is proposed by introducing a free energy function that takes into account the actual disc annulus structure in relation with the surrounding biochemical environment. The response is assumed to be dominated by the viscoelastic contribution of the extracellular matrix, the elastic contribution of the oriented collagen fibers and the osmo-induced volumetric contribution of the internal fluid content variation. The regional dependence of the disc annulus response due to variation in fibers content/orientation allows a micromechanical treatment of the soft tissue. A finite element model of the annulus specimen is designed while taking into consideration the 'interlamellar' ground substance zone between lamellae of the layered soft tissue. The kinetics is designed using full-field strain measurements performed on specimens extracted from two disc annulus regions and tested under different osmotic conditions. The time-dependency of the tissue response is reported on stress-free volumetric changes, on hysteretic stress and transversal strains during quasi-static stretching at different strain-rates and on their temporal changes during an interrupted stretching. Considering the effective contributions of the internal fluid transfer and the extracellular matrix viscosity, the microstructure-based chemo-mechanical model is found able to successfully reproduce the significant features of the macro-response and the unusual transversal behavior including the strong regional dependency from inner to outer parts of the disc: Poisson's ratio lesser than 0 (auxetic) in lamellae plane, higher than 0.5 in fibers plane, and their temporal changes towards usual values (between 0 and 0.5) at chemo-mechanical equilibrium. The underlying time-dependent mechanisms occurring in the tissue are analyzed via the local numerical fields and important insights about the effective role of the interlamellar zone are revealed for the different disc localizations. STATEMENT OF SIGNIFICANCE: The structural complexity of the annulus fibrosus has only been appreciated through recent experimental contributions and a constitutive representation that considers as well regional effect, chemical sensitivity and time-dependency of the unusual transversal behavior has not yet been developed. Here, a microstructure-based chemo-viscoelastic model is developed to highlight the interlamellar-induced time-dependent response by means of a two-scale strategy. The model provides important insights about the origin of the time-dependent phenomena in disc annulus along with regional dependency, essential for understanding disc functionality.

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.

Crystallinity dependency of the time-dependent mechanical response of polyethylene: application in total disc replacement

Degeneration of the intervertebral disc (IVD) is a leading source of chronic low back pain or neck pain, and represents the main cause of long-term disability worldwide. In the aim to relieve pain, total disc replacement (TDR) is a valuable surgical treatment option, but the expected benefit strongly depends on the prosthesis itself. The present contribution is focused on the synthetic mimic of the native IVD in the aim to optimally restore its functional anatomy and biomechanics, and especially its time-dependency. Semi-crystalline polyethylene (PE) materials covering a wide spectrum of the crystallinity are used to propose new designs of TDR. The influence of the crystallinity on various features of the time-dependent mechanical response of the PE materials is reported over a large strain range by means of dynamic mechanical thermo-analysis and video-controlled tensile mechanical tests. The connection of the stiffness and the yield strength with the microstructure is reported in the aim to propose a model predicting the crystallinity dependency of the response variation with the frequency. New designs of TDR are proposed and implemented into an accurate computational model of a cervical spine segment in order to simulate the biomechanical response under physiological conditions. Predicted in-silico motions are found in excellent agreement with experimental data extracted from published in-vitro studies under compression and different neck movements, namely, rotation, flexion/extension and lateral bending. The simulation results are also criticized by analyzing the local stresses and the predicted biomechanical responses provided by the different prosthetic solutions in terms of time-dependency manifested by the hysteretic behavior under a cyclic movement and the frequency effect.

Osmo-inelastic response of the intervertebral disc

The intervertebral disc exhibits a complex inelastic response characterized by relaxation, hysteresis during cyclic loading and rate dependency. All these inelastic phenomena depend on osmotic interactions between disc tissues and their surrounding chemical environment. Coupling between osmotic and inelastic effects is not fully understood, so this article aimed to study the influence of chemical conditions on the inelastic behaviour of the intervertebral disc in response to different modes of loading. A total of 18 non-frozen ‘motion segments’ (two vertebrae and the intervening soft tissues) were dissected from the cervical spines of mature sheep. The motion segments were loaded in tension, compression and torsion at various loading rates and saline concentrations. Analysis of variance showed that saline concentration significantly influenced inelastic effects in tension and especially in compression (p < 0.05), but not in torsion. Opposite effects were seen in tension and compression. An interpretation of the underlying osmo-inelastic mechanisms is proposed in which two sources of inelastic effects are identified, that is, extracellular matrix rearrangements and fluid exchange created by osmosis.

New insights into the viscoelastic and failure mechanical properties of the elastic fiber network of the inter-lamellar matrix in the annulus fibrosus of the disc

The mechanical role of elastic fibers in the inter-lamellar matrix (ILM) is unknown; however, it has been suggested that they play a role in providing structural integrity to the annulus fibrosus (AF). Therefore, the aim of this study was to measure the viscoelastic and failure properties of the elastic fiber network in the ILM of ovine discs under both tension and shear directions of loading. Utilizing a technique, isolated elastic fibers within the ILM from ovine discs were stretched to 40% of their initial length at three strain rates of 0.1% s^-1 (slow), 1% s^-1 (medium) and 10% s^-1 (fast), followed by a ramp test to failure at 10% s^-1. A significant strain-rate dependent response was found, particularly at the fastest rate for phase angle and normalized stiffness (p < 0.001). The elastic fibers in the ILM demonstrated a significantly higher capability for energy absorption at slow compared to medium and fast strain rates (p < 0.001). These finding suggests that the elastic fiber network of the ILM exhibits nonlinear elastic behavior. When tested to failure, a significantly higher normalized failure force was found in tension compared to shear loading (p = 0.011), which is consistent with the orthotropic structure of elastic fibers in the ILM. The results of this study confirmed the mechanical contribution of the elastic fiber network to the ILM and the structural integrity of the AF. This research serves as a foundation for future studies to investigate the relationship between degeneration and ILM mechanical properties.

STATEMENT OF SIGNIFICANCE

The mechanical role of elastic fibres in the inter-lamellar matrix (ILM) of the disc is unknown. The viscoelastic and failure properties of the elastic fibre network in the ILM in both tension and shear directions of loading was measured for the first time. We found a strain-rate dependent response for the elastic fibres in the ILM. The elastic fibres in the ILM demonstrated a significantly higher capability for energy absorption at slow compared to medium and fast strain rates. When tested to failure, a significantly higher normalized failure force was found in tension compared to shear loading, which is consistent with the orthotropic structure of elastic fibres in the ILM.