The biological basis of degenerative disc disease: proteomic and biomechanical analysis of the canine intervertebral disc

IntroductionIn the present study, we sought to quantify and contrast the secretome and biomechanical properties of the non-chondrodystrophic (NCD) and chondrodystrophic (CD) canine intervertebral disc (IVD) nucleus pulposus (NP).MethodsWe used iTRAQ proteomic methods to quantify the secretome of both CD and NCD NP. Differential levels of proteins detected were further verified using immunohistochemistry, Western blotting, and proteoglycan extraction in order to evaluate the integrity of the small leucine-rich proteoglycans (SLRPs) decorin and biglycan. Additionally, we used robotic biomechanical testing to evaluate the biomechanical properties of spinal motion segments from both CD and NCD canines.ResultsWe detected differential levels of decorin, biglycan, and fibronectin, as well as of other important extracellular matrix (ECM)-related proteins, such as fibromodulin and HAPLN1 in the IVD NP obtained from CD canines compared with NCD canines. The core proteins of the vital SLRPs decorin and biglycan were fragmented in CD NP but were intact in the NP of the NCD animals. CD and NCD vertebral motion segments demonstrated significant differences, with the CD segments having less stiffness and a more varied range of motion.ConclusionsThe CD NP recapitulates key elements of human degenerative disc disease. Our data suggest that at least some of the compromised biomechanical properties of the degenerative disc arise from fibrocartilaginous metaplasia of the NP secondary to fragmentation of SLRP core proteins and associated degenerative changes affecting the ECM. This study demonstrates that the degenerative changes that naturally occur within the CD NP make this animal a valuable animal model with which to study IVD degeneration and potential biological therapeutics.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-015-0733-z) contains supplementary material, which is available to authorized users.     

Accelerated cellular senescence in degenerate intervertebral discs: a possible role in the pathogenesis of intervertebral disc degeneration

Current evidence implicates intervertebral disc degeneration as a major cause of low back pain, although its pathogenesis is poorly understood. Numerous characteristic features of disc degeneration mimic those seen during ageing but appear to occur at an accelerated rate. We hypothesised that this is due to accelerated cellular senescence, which causes fundamental changes in the ability of disc cells to maintain the intervertebral disc (IVD) matrix, thus leading to IVD degeneration. Cells isolated from non-degenerate and degenerate human tissue were assessed for mean telomere length, senescence-associated β-galactosidase (SA-β-gal), and replicative potential. Expression of P16 ^INK4A (increased in cellular senescence) was also investigated in IVD tissue by means of immunohistochemistry. RNA from tissue and cultured cells was used for real-time polymerase chain reaction analysis for matrix metalloproteinase-13, ADAMTS 5 (a disintegrin and metalloprotease with thrombospondin motifs 5), and P16 ^INK4A . Mean telomere length decreased with age in cells from non-degenerate tissue and also decreased with progressive stages of degeneration. In non-degenerate discs, there was an age-related increase in cellular expression of P16 ^INK4A . Cells from degenerate discs (even from young patients) exhibited increased expression of P16 ^INK4A , increased SA-β-gal staining, and a decrease in replicative potential. Importantly, there was a positive correlation between P16 ^INK4A and matrix-degrading enzyme gene expression. Our findings indicate that disc cell senescence occurs in vivo and is accelerated in IVD degeneration. Furthermore, the senescent phenotype is associated with increased catabolism, implicating cellular senescence in the pathogenesis of IVD degeneration.     

Morphometric analysis of lumbar vertebrae in extinct Malagasy strepsirrhines

Previous research on subfossil lemurs has revealed much about the positional behavior of these extinct strepsirrhines, but a thorough quantitative analysis of their vertebral form and function has not been performed. In this study, 156 lumbar vertebrae of Pachylemur, Archaeolemur, Megaladapis, Mesopropithecus, Babakotia, and Palaeopropithecus (11 species in all) were compared to those of 26 species of extant strepsirrhines and haplorhines. Lumbar shape was compared among species, using a principal components analysis (PCA) in conjunction with selected vertebral indices. The first principal component revealed strong separation between Palaeopropithecus at one extreme, and Archaeolemur/Pachylemur at the other, with Babakotia, Mesopropithecus, and Megaladapis in an intermediate position. Palaeopropithecus has markedly shorter spinous processes and wider laminae than do the other subfossil taxa, consistent with sloth-like, inverted suspensory postures. The moderately reduced lumbar spinous processes of Babakotia, Mesopropithecus, and Megaladapis are convergent with those of lorisids and Pongo, reflecting antipronogrady, but a less specialized adaptation than that of Palaeopropithecus. Archaeolemur and Pachylemur share relatively elongated spinous processes, in conjunction with other features (e.g., transverse process orientation and relatively short vertebral bodies) indicative of pronograde, quadrupedal locomotion characterized by reduced agility. All subfossil taxa exhibit adaptations emphasizing lumbar spinal stability (e.g., relatively short vertebral bodies, and transverse processes that are not oriented ventrally); we believe this probably reflects convergent mechanical demands connected to large body size, irrespective of specific locomotor mode. Reconstructions of positional behavior in subfossil lemurs based on lumbar vertebrae are largely consistent with those based on other aspects of the postcrania.     

Accumulation of platinum in the intervertebral discs and vertebrae of ovarian tumor-bearing patients treated with cisplatin

Platinum was determined by the inductively coupled plasma mass spectrometry (ICP-MS) in the intervertebral discs and vertebrae of ovarian tumor bearing patients treated withcis-diamminedichloro-platinum (II) (cisplatin). Platinum was 0.05 ng/mL at the absolute detection limit, and platinum was undetectable in the intervertebral discs and vertebrae of human specimens without cisplatin treatments. On the other hand, platinum was detected in the intervertebral discs and vertebrae of patients administered cisplatin, and platinum concentration was at levels of 1.06–10.31 μg/g dry tissue in the intervertebral discs and 0.60–1.28 μg/g dry tissue in the vertebrae, respectively. The platinum level of intervertebral discs was 4.3-fold higher than that of the vertebrae. Thus, platinum accumulates greatly in the intervertebral discs and somewhat in the vertebrae after administering cisplatin to patients for therapy.     

Biomechanical comparison between fusion of two vertebrae and implantation of an artificial intervertebral disc

Surgical treatments for lower back pain can be distributed into two main groups: fusion (arthrodesis) and disc replacement (arthroplasty). The objective of this study was to compare, under severe loading conditions, the biomechanics of the lumbar spine treated either by fusion or total disc replacement (TDR). A three-dimensional model of a two-level ligamentous lumbar segment was created and simulated through static analyses with the finite-element method (FEM) software ABAQUS. The model was validated by comparing mobility, pressure on the facets, force in the ligaments, maximum stresses, disc bulge, and endplate deflection with measured data given in the literature. The FEM analysis predicted that the mobility of the model after arthrodesis on the upper level was reduced in all rotational degrees of freedom by an average of approximately 44%, relative to healthy normal discs. Conversely, the mobility of the model after TDR on the upper level was increased in all rotational degrees of freedom by an average of approximately 52%. The level implanted with the artificial disc showed excessive ligament tensions (greater than 500 N), high facet pressures (greater than 3 MPa), and a high risk of instability. The mobility and the stresses in the level adjacent to the arthroplasty were also increased. In conclusion, the model for an implanted movable artificial disc illustrated complications common to spinal arthroplasty and showed greater risk of instability and further degeneration than predicted for the fused model. This modeling technique provides an accurate means for assessing potential biomechanical risks and can be used to improve the design of future artificial intervertebral discs.     

Thermal analysis of the human intervertebral disc

Intervertebral disc (IVD) degeneration is one of the most common musculuskeletal disorders affecting western society. Degeneration alters the morphology and the mechanical properties of the discs. According to previous reports, DSC proved to be a suitable method for the demonstration of thermal consequences of local as well as global conformational changes in the structure of the human intervertebral discs. In the present study, a wide spectrum of degenerated IVD was examined by DSC. The results suggest that definitive differences exist between the stages of disc degeneration in calorimetric measures.     

Mechanical profiling of intervertebral discs

Despite recent advances in imaging diagnostic technology and additional treatment options our ability to prevent or inhibit discogenic back pain has not drastically improved. The challenge of linking early degenerative patterns to dysfunction and pain remains. Using a novel material testing device designated the tissue diagnostic instrument (TDI) we measured the local stiffness and strain energy absorption in the radial direction of 13 intact intervertebral discs; effectively generating a mechanical profile of each disc. Prior to measuring mechanical properties, an MR image was taken of each spine segment and the discs were radiologically scored according to the Pfirrmann scale. After testing, a sagittal portion of each L1–L2 disc was excised from each of four spines for histology. No significant correlations were found between Pfirrmann grade and mechanical data. However, polarized light microscopy images of disc sections indicated correlations between local tissue modulus measured with the TDI and the clarity and density of lamellar striations.     

Loads distributed in vivo among vertebrae,muscles, spinal ligaments,and intervertebral discs in a passively flexed lumbar spine

Biomechanics and Modeling in Mechanobiology - The load distribution among lumbar spinal structures—still an unanswered question—has been in the focus of this hybrid experimental and...     

A hyperelastic and almost incompressible material model as an approach to intervertebral disc analysis

The literature contains many references to intervertebral disc analysis using experimental techniques which represent a basis for analytical and numerical approaches. The scatter of experimental results may be one reason for difficulties in interpretation. Apart from analytical approaches, which often seen inadequate to deal with the peculiar complexities of the problem, numerical techniques are reliable and can lead to significant results. In this work, a formulation based on the finite element method is described, adopting a nonlinear model with hyperelastic material configuration. Particular attention is paid to modelling the material constituting the nucleus, involving incompressibility characteristics and avoiding simulation techniques. This approach allows the mechanical behaviour of the real configuration of the disc to be investigated and also provides a reliable analysis of disc degeneration phenomena. The theoretical and operational aspects of the formulation are reported. The results obtained are compared with responses from various numerical and experimental data.     

Gene expression profile analysis of human intervertebral disc degeneration

In this study, we used microarray analysis to investigate the biogenesis and progression of intervertebral disc degeneration. The gene expression profiles of 37 disc tissue samples obtained from patients with herniated discs and degenerative disc disease collected by the National Cancer Institute Cooperative Tissue Network were analyzed. Differentially expressed genes between more and less degenerated discs were identified by significant analysis of microarray. A total of 555 genes were significantly overexpressed in more degenerated discs with a false discovery rate of < 3%. Functional annotation showed that these genes were significantly associated with membrane-bound vesicles, calcium ion binding and extracellular matrix. Protein-protein interaction analysis showed that these genes, including previously reported genes such as fibronectin, COL2A1 and β-catenin, may play key roles in disc degeneration. Unsupervised clustering indicated that the widely used morphology-based Thompson grading system was only marginally associated with the molecular classification of intervertebral disc degeneration. These findings indicate that detailed, systematic gene analysis may be a useful way of studying the biology of intervertebral disc degeneration.     

Degeneration of the intervertebral disc

The intervertebral disc is a cartilaginous structure that resembles articular cartilage in its biochemistry, but morphologically it is clearly different. It shows degenerative and ageing changes earlier than does any other connective tissue in the body. It is believed to be important clinically because there is an association of disc degeneration with back pain. Current treatments are predominantly conservative or, less commonly, surgical; in many cases there is no clear diagnosis and therapy is considered inadequate. New developments, such as genetic and biological approaches, may allow better diagnosis and treatments in the future.     

Effect of intervertebral disc degeneration on disc cell viability: a numerical investigation

Degeneration of the intervertebral disc may be initiated and supported by impairment of the nutrition processes of the disc cells. The effects of degenerative changes on cell nutrition are, however, only partially understood. In this work, a finite volume model was used to investigate the effect of endplate calcification, water loss, reduction of disc height and cyclic mechanical loading on the sustainability of the disc cell population. Oxygen, lactate and glucose diffusion, production and consumption were modelled with non-linear coupled partial differential equations. Oxygen and glucose consumption and lactate production were expressed as a function of local oxygen concentration, pH and cell density. The cell viability criteria were based on local glucose concentration and pH. Considering a disc with normal water content, cell death was initiated in the centre of the nucleus for oxygen, glucose, and lactate diffusivities in the cartilaginous endplate below 20% of the physiological values. The initial cell population could not be sustained even in the non-calcified endplates when a reduction of diffusion inside the disc due to water loss was modelled. Alterations in the disc shape such as height loss, which shortens the transport route between the nutrient sources and the cells, and cyclic mechanical loads, could enhance cell nutrition processes.     

Diffusion of antibiotics in intervertebral disc

Delivering charged antibiotics to the intervertebral disc is challenging because of the avascular, negatively charged extracellular matrix (ECM) of the tissue. The purpose of this study was to measure the apparent diffusion coefficient of two clinically relevant, charged antibiotics, vancomycin (positively charged) and oxacillin (negatively charged) in IVD. A one-dimensional steady state diffusion experiment was employed to measure the apparent diffusion coefficient of the two antibiotics in bovine coccygeal annulus fibrosus (AF) tissue. The averaged apparent diffusion coefficient for vancomycin under 20% compressive strain was 7.94 ± 2.00 × 10⁻¹² m²/s (n = 10), while that of oxacillin was 2.26 ± 0.68 × 10⁻¹⁰ m²/s (n = 10). A student’s t-test showed that the diffusivity of vancomycin was significantly lower than that of oxacillin. This finding may be attributed to two factors: solute size and possible binding effects. Vancomycin is approximately 3 times larger in molecular weight than oxacillin, meaning that steric hindrance likely plays a role in the slower transport. Reversible binding between positive vancomycin and the negative ECM could also slow down the rate of diffusion. Therefore, more investigation is necessary to determine the specific relationship between net charge on antibiotic and diffusion coefficients in IVD. This study provides essential quantitative information regarding the transport rates of antibiotics in the IVD, which is critical in using computational modeling to design effective strategies to treat disc infection.