REVIEW article

Front. Cell Dev. Biol., 29 June 2022

Sec. Molecular and Cellular Pathology

Volume 10 - 2022 | https://doi.org/10.3389/fcell.2022.924692

Immuno-Modulatory Effects of Intervertebral Disc Cells

  • 1. Tissue Engineering for Orthopaedics and Mechanobiology, Bone and Joint Program, Department for BioMedical Research (DBMR), Faculty of Medicine, University of Bern, Bern, Switzerland

  • 2. Department of Orthopaedic Surgery and Traumatology, Inselspital, Bern University Hospital, Medical Faculty, University of Bern, Bern, Switzerland

  • 3. BCN MedTech, Universitat Pompeu Fabra, Barcelona, Spain

  • 4. Biomolecular Sciences Research Centre, Sheffield Hallam University, Sheffield, United Kingdom

  • 5. ProtATonce Ltd., Athens, Greece

  • 6. School of Mechanical Engineering, National Technical University of Athens, Zografou, Greece

  • 7. GIGA In Silico Medicine, University of Liège, Liège, Belgium

  • 8. Skeletal Biology and Engineering Research Center, KU Leuven, Leuven, Belgium

  • 9. Twin Research and Genetic Epidemiology, St Thomas’ Hospital, King’s College London, London, United Kingdom

  • 10. Institut Hospital Del Mar D’Investigacions Mèdiques (IMIM), Barcelona, Spain

  • 11. Biomechanics Research Unit, KU Leuven, Leuven, Belgium

  • 12. Department of Biomedical Engineering, Rochester Institute of Technology, Rochester, NY, United States

  • 13. Spine Center, Schön Klinik München Harlaching Academic Teaching Hospital and Spine Research Institute of the Paracelsus Private Medical University Salzburg (Austria), Munich, Germany

Article metrics

View details

78

Citations

12,1k

Views

4,1k

Downloads

Abstract

Low back pain is a highly prevalent, chronic, and costly medical condition predominantly triggered by intervertebral disc degeneration (IDD). IDD is often caused by structural and biochemical changes in intervertebral discs (IVD) that prompt a pathologic shift from an anabolic to catabolic state, affecting extracellular matrix (ECM) production, enzyme generation, cytokine and chemokine production, neurotrophic and angiogenic factor production. The IVD is an immune-privileged organ. However, during degeneration immune cells and inflammatory factors can infiltrate through defects in the cartilage endplate and annulus fibrosus fissures, further accelerating the catabolic environment. Remarkably, though, catabolic ECM disruption also occurs in the absence of immune cell infiltration, largely due to native disc cell production of catabolic enzymes and cytokines. An unbalanced metabolism could be induced by many different factors, including a harsh microenvironment, biomechanical cues, genetics, and infection. The complex, multifactorial nature of IDD brings the challenge of identifying key factors which initiate the degenerative cascade, eventually leading to back pain. These factors are often investigated through methods including animal models, 3D cell culture, bioreactors, and computational models. However, the crosstalk between the IVD, immune system, and shifted metabolism is frequently misconstrued, often with the assumption that the presence of cytokines and chemokines is synonymous to inflammation or an immune response, which is not true for the intact disc. Therefore, this review will tackle immunomodulatory and IVD cell roles in IDD, clarifying the differences between cellular involvements and implications for therapeutic development and assessing models used to explore inflammatory or catabolic IVD environments.

Introduction

Epidemiology of Intervertebral Disc Degeneration

Low back pain (LBP) is the largest cause of morbidity worldwide, affecting approximately 80% of people from Western countries during their lifetime and resulting in 5 million disability-adjusted life-years in young adults (GBD 2017 Disease and Injury Incidence and Prevalence Collaborators, 2018). Lower intervertebral disc degeneration (IDD) is the cause of around half of all LBP cases in young adults; however not all cases of IDD result in LBP (Baumgartner et al., 2021c). Although IDD prevalence increases progressively with age, IDD is common in subjects younger than 30 years old, conveying those various other factors besides age, such as excessive or uneven mechanical load, obesity, genetics, nutrition, trauma, and gender are involved (Hoogendoorn et al., 2000; Paassilta et al., 2001; Pincus et al., 2002; Cheung et al., 2009; Samartzis et al., 2012; Teraguchi et al., 2014; Parenteau et al., 2021). For example, studies have shown that women experience greater pain and disability than men when they are treated for IDD (MacLean et al., 2020). Additionally, LBP prevalence in females after menopause further increases in comparison to men at comparable ages (Wáng et al., 2016). Further, it is unclear whether occupation-related heavy physical loading is an important risk factor for IDD, as studies have contradictory conclusions (Bongers et al., 1990; Videman and Battié, 1999). Some studies have found IDD is significantly more common in athletes compared to the general population (Swärd et al., 1991). However, various twin studies that have been conducted suggest that occupation or sport related risk factors have only a minor role in IDD, while genetic influences were found to play a greater role in predicting degeneration (Battié et al., 2004). On the other hand, obesity is associated with increased IDD severity and extent, likely due to altered biomechanical and/or biological processes such as those driven by adipokines (Samartzis et al., 2012; Li W. et al., 2022). Due to the complexity and multifactorial nature of IDD, the initiating and risk factors are poorly understood, which critically hampers proper LBP patient stratification and limits the development of personalized therapies.

The Structure of the Intervertebral Disc

The intervertebral disc (IVD) is the largest avascular organ in the human body with blood vessels only present in the outer annulus fibrosus (AF) and boney end plates, with all nutrient and waste exchange taking place via diffusion through the dense extracellular matrix of the disc (Urban 2002). Located between the vertebrae within the spine, the IVD consists of three highly hydrated, major tissues: 1) the nucleus pulposus (NP), 2) the annulus fibrosus (AF), and 3) the cartilage endplate (CEP). The central and proteoglycan-rich NP lies between the cranial and caudal CEPs and is laterally constrained by the peripheral and fiber-reinforced AF (Figure 1). This specialized composition and structure of the IVD ensures both trunk movements and resistance to high mechanical loads. The normal human IVD contains nucleus pulposus cells and annulus fibrosus cells within the NP and AF, respectively, with AF cells becoming more elongated and fibroblast-like towards the periphery. Cells occupy 1% volume of the disc, though are crucial in maintaining the balance between anabolic activity such as the production of proteoglycans and collagens type I and II, and the procatabolic effects of factors involved in ECM turnover, including metalloproteinases, prostaglandins, and nitric oxide (Kang et al., 1997). Furthermore, mechanical loads are thought to influence ECM homeostasis, where both excessive and insufficient loads lead to catabolism (Vergroesen et al., 2015). Due to the avascularity of the IVD, the environment is hypoxic, where the oxygen tension in the IVD is considered between 1 and 5% (Yao et al., 2016; Yao et al., 2017).

FIGURE 1

In comparison to NP and AF tissues, the CEP receives far less attention in the literature; however, it is a vital tissue when discussing LBP. Lakstins et al. (2021) demonstrated that imperfections and weakness in the CEP can be a better anticipator of pain than IVD degeneration because chemical changes to the CEP are directly related to intervertebral disc degeneration (IDD) (Yao et al., 2016). The CEP is rich in collagen type II (Yao et al., 2016) and performs both mechanical and chemical functions (Roberts et al., 1989; Lakstins et al., 2021). Mechanically, the CEP acts as a physical filter preventing macromolecules from escaping the disc through the subchondral bone and is considered important in controlling the hydration of the disc under mechanical loads (Roberts et al., 1989; Moore 2000; Ruiz Wills et al., 2018). Chemically, the CEP allows metabolites, small molecules and waste to travel between the IVD and neighboring blood vessels in the bony endplates (Roberts et al., 1989; Turgut et al., 2003; Yao et al., 2016; Ruiz Wills et al., 2018; Zuo et al., 2019; Sun et al., 2020; Lakstins et al., 2021).

The diffusivity of solutes through the CEP and towards the IVD depends greatly on their size and ionic charge. The healthy IVD is negatively charged due to the high proteoglycan concentration (Moore 2000; Pfannkuche et al., 2020). Therefore, only small, neutrally charged solutes such as glucose and oxygen, as well as cations such as sodium and calcium can penetrate the disc, but small anions such as sulphate and chloride ions can only cross through the CEP. In turn, large, neutrally charged solutes such as antibodies and enzymes usually cannot penetrate the healthy IVD (Moore 2000).

Intervertebral Disc Degeneration

Regarding disc morphology, as the IVD degenerates, it becomes more difficult to distinguish the boundaries between the AF and the NP. This loss of a distinct boundary worsens with age, as the nucleus loses its gel-like quality and becomes more fibrotic (Buckwalter 1995) which was seen as a common degenerative feature across all species (Dahia et al., 2021). Another significant biochemical change during disc degeneration is the loss of proteoglycans, which are necessary to provide the osmotic resistance for the IVD to withstand compressive loads and keep the disc hydrated (Knudson and Knudson, 2001). Such significant changes (loss of water content (Lyons et al., 1981) and disc height (Frobin et al., 2001)) in disc behavior strongly influence other spinal structures and may negatively impact their function and predispose them to injury.

During IDD, the CEP becomes thinner and fissured, with lower collagen and glycosaminoglycan (GAG) content (Turgut et al., 2003; Sun et al., 2020). This change in morphology affects the physiology and the performance of the CEP (Hamilton et al., 2006; Roberts et al., 2006) altering its permeability (Roberts et al., 2006; Yao et al., 2016). Furthermore, the CEP can lose its connection to the vasculature (Moon et al., 2013), which immunohistochemistry has shown leads to blood vessel and nerve fiber infiltration into the IVD through the CEP and subchondral bone and through fissures in the AF (Freemont et al., 1997; Roberts et al., 2006; Binch et al., 2015b; Yao et al., 2016; Lama et al., 2018; Sun et al., 2020). Moreover, the crosstalk between IVD and the bone marrow is facilitated due to the CEP damage (Dudli et al., 2016), causing possible adjacent “Modic discs”. Modic changes (MC) are defined as magnetic resonance imaging (MRI) signal alterations in the vertebral bone marrow close to a degenerated disc. There are several different types of MC, with MC1 fibrotic lesion having the highest association with pain, followed by MC2. MC3 are rare and often asymptomatic. MC1 and MC2 are commonly accompanied by persistent inflammatory stimulus. In addition, MC related pain could be related to the neovascularization and neurogenesis due to the increase in growth factor expression by blood vessels and disc cells and inflammatory cytokines (Rätsep et al., 2013; Sun et al., 2013; Dudli et al., 2018) which lead to increased expression of neurotrophic factors (Freemont et al., 2002; Purmessur et al., 2008; Binch et al., 2014) (Figure 1).

During disc degeneration, the balance between anabolism and catabolism is dysregulated, showing decreased synthesis of normal matrix, of collagen type II and aggrecan and increased presence of matrix degrading enzymes, reviewed previously by Baumgartner et al. (2021c). Moreover, several studies have reported decreased NP cell proliferation under catabolic cytokine stimulation (Wang et al., 2013; Li et al., 2019; Lin and Lin, 2020). Similarly, during the shift from anabolic to catabolic cell activity in the disc, the presence of these cytokines is also related to NP and AF cell apoptosis, (Hu et al., 2017; Yu et al., 2018; Zhang J. et al., 2019; Zhang S. et al., 2019).

These changes have been shown, at least in part, to be modulated by pro-catabolic cytokines in numerous studies, which are often referred to as inflammatory features, in the literature. However, since these factors are produced by native disc cells (NP, AF and CEP) in intact discs, this catabolic response can be easily misconstrued as an inflammatory response. Therefore, the aim of this review is to tackle immunomodulatory and IVD cell roles in IDD and clarifying the differences between cellular involvements. Furthermore, different in-silico, in-vivo and in-vitro models used to explore inflammatory or catabolic IVD environments will be discussed.

Cross-Talk Between the Immune System and IVD in IDD

As mentioned, the IVD is the largest avascular organ with blood vessels only present in the outer AF and bony end plates, so all metabolite exchange takes place via diffusion through the dense extracellular matrix (ECM) of the IVD. The dense ECM of the IVD inhibits blood vessel ingrowth both mechanically by having a high physical pressure, and chemically through high proteoglycan concentration (Johnson et al., 2002; Johnson et al., 2005), which combined with secretory inhibitors prevent nerve and blood vessel ingrowth in non-degenerate discs (Tolofari et al., 2010; Binch et al., 2015a). The AF and the CEP, along with the secretory inhibitors of angiogenesis, are defined as the blood-NP barrier (BNB), which strongly isolates the NP from the circulation and thus the host immune system (Sun et al., 2020).

Where both AF and CEP are intact, the IVD has been described as an immuno-privileged tissue (Sun et al., 2020) with a lack of immune cells (Figure 2). However, this is often confused as the native cells of the disc (i.e., the NP, AF and CEP cells) have been shown to take on roles and markers classically expressed by immune cells (Le Maitre et al., 2005; Jones et al., 2008; Phillips et al., 2013a; Risbud and Shapiro, 2013), and thus have been described by some authors as immune responses or inflammation. However, such activity of native IVD cells is not true inflammation. Therefore, distinguishing which cases of IDD involve an immune response is important as different clinical interventions and treatments would be required.

FIGURE 2

Native disc cells produce a plethora of cytokines and chemokines which are upregulated during disc degeneration and have been shown to drive many catabolic events in the IVD (Weiler et al., 2005; Le Maitre et al., 2005; Hoyland et al., 2008; Phillips et al., 2013a; Phillips et al., 2015). A shift to catabolism is at least in part driven by the increased production of cytokines in the disc by the native cells (in an intact disc) and a combination of inflammatory cells and native disc cells following CEP and AF rupture. Phillips et al. (2015) demonstrated that NP cells express a number of cytokine and chemokine receptors and are thus able to respond in a paracrine and autocrine manner (Figure 2). Caused by different, yet not fully understood mechanisms, disc cells upregulate the expression of inflammatory cytokines such as IL-1, TNFα, IL-6, IL-8 and IL-17 amongst others creating a cytokine rich catabolic environment. IL-1 has been shown to drive the catabolic events during disc degeneration (Le Maitre et al., 2005; Weiler et al., 2005; Hoyland et al., 2008; Phillips et al., 2013a; Phillips et al., 2015). Whilst other cytokines and chemokines (e.g., MCP-1, TNFα, IL-8) produced in the disc appear to have limited effects on the native disc cells due to lack of receptor expression in vivo (Le Maitre et al., 2007; Phillips et al., 2015), they undoubtedly diffuse to the surrounding tissues leading to increased inflammation in local tissues, and drive cellular infiltration following AF and CEP rupture and increased sensitisation of nerves (Ye et al., 2022). Such cytokines can stimulate specific intracellular signaling pathways that further enhance the degenerative process (Daniels et al., 2017; Suyama et al., 2018) and upregulate matrix-degrading enzymes known as matrix metalloproteinases (MMPs) and a disintegrant and metalloproteinase with thrombospondin motifs (ADAMTS), specifically MMP- 1, 2, 3, 9, 13 and ADAMTS-4, 5 (Baumgartner et al., 2021c). In later phases of IDD, these cytokines can upregulate neurotrophic and angiogenic factors, which could lead to further nerve and blood vessel ingrowth (Purmessur et al., 2008; Lee et al., 2011; Binch et al., 2014; Krock et al., 2014).

Remarkably, some of these cytokines, such as IL-1, have also been shown to be expressed in cells from non-degenerate discs and display roles in maintaining normal homeostasis (Le Maitre et al., 2005; Phillips et al., 2015). Indeed, if the IL-1 agonists are knocked out during development, IDD can be induced (Gorth et al., 2019). Thus, IL-1 plays a role as a normal regulatory mechanism during IVD homeostasis, which becomes imbalanced during IDD (Le Maitre et al., 2005) (Figure 2). Native NP cells have also been shown to take on other roles normally associated with immune cells. such as phagocytosis: Jones et al. (2008) observed the capacity of bovine NP cells to ingest latex beads at least as efficiently as phagocytic cells and ingested apoptotic cells. This capability could be of great physiological relevance to maintain a healthy disc, as it may prevent inflammation triggered by the release of toxic or immunogenic intracellular content by apoptotic cells (Fadok, 1999). Clearly, the suggestion from some reviews that cytokine production within the disc is solely from immune cells is inaccurate (Ye et al., 2022). However, when the AF or CEP becomes ruptured or fissures occur during injury or disc degeneration this provides a route for blood vessel ingrowth and migration of immune cells into the intervertebral disc. Within these “non-intact” IVDs, immune cells will migrate including T cells (CD4+, CD8+), B cells, macrophages, neutrophils and mast cells (Risbud and Shapiro, 2013) (Figure 2).

These immune cells then contribute to an inflamed environment in the disc, leading to further increases in cytokine and chemokine expression (Phillips et al., 2015). This leads to a viscous circle of inflammatory driven catabolism which acts synergistically with the native IVD cells to cause accelerated ECM breakdown (Figure 3) (Risbud and Shapiro, 2013). These cytokines and chemokines play a number of roles within this disc, including direct actions on NP, AF and CEP cells where their receptors are present (Le Maitre et al., 2005; Le Maitre et al., 2007; Phillips et al., 2015). They will also likely diffuse out of the IVD leading to increased cellular migration to the disc (Pattappa et al., 2014), or sensitization of local nerve roots (Leung and Cahill, 2010; Johnson et al., 2015).

FIGURE 3

The Shift From Anabolism to Catabolism

A Complex Interplay of Microenvironment, Biomechanics, Genetics and Epigenetics, Bacterial Infection and IVD Cells

As discussed during disc degeneration, there is a shift in metabolism from anabolism (matrix synthesis) to catabolism (matrix degradation), and this shift to catabolism is accompanied by increased production of neurotrophic and angiogenic factors which lead to nerve and blood vessel ingrowth leading to inflammation and increased pain sensation in the disc (Baumgartner et al., 2021c). There remains a poor understanding as to the initiating factors involved in this switch from anabolism to catabolism in disc degeneration, which is likely due to multifactorial processes including the disc microenvironment, biomechanics, genetics and epigenetics, and even bacterial infection of the disc and the gut microbiome (Figure 3) (Li W. et al., 2022).

Altered Disc Microenvironment

The IVD microenvironment is commonly described as harsh due to its limited nutrition (glucose and oxygen), low pH, and large changes in tissue osmolarity (Urban 2002). These factors not only impair the success of cell therapies (e.g., mesenchymal stromal cell injection) (Loibl et al., 2019; Williams et al., 2021) but can also negatively affect resident IVD cells and thus contribute to the catabolic-inflammatory shift observed during degeneration.

Limited nutrition and tissue acidity are a result of the avascular nature of the IVD (Hukins, 1988). Glucose and oxygen transport into the IVD, as well as the removal of cellular waste products such as lactic acid (which contributes to the drop in tissue pH), are hence dependent on diffusion via the capillaries in the endplate. Degeneration-associated calcification as well as a reduction in the density of the bone marrow contact channels in the endplates might further impair these transport mechanisms (Benneker et al., 2005; Chen et al., 2014). In the centre of the IVD, glucose levels, oxygen and pH can hence drop to 0.5 mM, 1% oxygen, and 6.5, respectively, although even lower values have also been reported (Urban et al., 2004). It is not clear whether glucose deprivation can activate the expression of pro-catabolic factors in NPC (Baumgartner et al., 2021b). However, multiple studies have demonstrated that glucose deprivation impairs NP cell proliferation and survival as well as proteoglycan synthesis/degradation and collagen synthesis (Bibby and Urban, 2004; Johnson et al., 2008; Chang et al., 2017; Saggese et al., 2020), and recent evidence highlights that non-coding RNAs (e.g., circ_0075062) may be involved in these processes (Chang et al., 2021). In contrast, oxygen deprivation alone was shown numerous times to have little effect on IVD cell proliferation or survival and mostly contributes to cell impairment prompted by other microenvironmental factors (Bibby and Urban, 2004; Johnson et al., 2008). Indeed, IVD cells are likely unaffected by hypoxia because of their robust and constitutive expression of hypoxia-inducible factor (HIF) 1, and more importantly, the inducible subunit HIF-1α (Merceron et al., 2014). Even more, hypoxia and HIF-1α were recently shown to attenuate the endoplasmic reticulum (ER) stress responses in NP cells (Novais et al., 2021). However, research related to rheumatoid arthritis demonstrated that the expression of Toll-like receptor 4 (TLR4) and TNF-α, but also of IL-10, are HIF-1-dependent processes in macrophages, indicating that IVD-infiltrating immune cells may be more affected by the hypoxic microenvironment than IVD cells themselves (Guo and Chen, 2020). Similar to low glucose concentrations, high lactic acid concentrations and the resulting drop in pH promotes cell death and a catabolic shift in mRNA expression (Horner and Urban, 2001; Bibby and Urban, 2004; Neidlinger-Wilke et al., 2012; Gilbert et al., 2016), likely via acid-sensing ion channels (Gilbert et al., 2016). Importantly, this response was more pronounced when cells were simultaneously exposed to low glucose levels (Bibby and Urban, 2004). Therefore, low glucose and high lactic acid levels, but not hypoxia, contribute to the catabolic shift observed during IVD degeneration. In general, there is a clear need to increase the number of experimental studies where different microenvironmental factors are combined, including both nutritional and pro-inflammatory cues (Baumgartner et al., 2021b).

Aside from low glucose, oxygen, and pH, the osmolality of the IVD is often considered a fourth harsh microenvironmental factor. This mostly refers to the relatively high osmolarity in the IVD, where 400 mOsm is considered iso-osmotic. However, it is important to note that IVD cells can be exposed to a wide range of tissue osmolarity, and these changes are more likely to affect IVD cell behaviour than the iso-osmotic condition. It can drop as low as 300 mOsm with the loss of proteoglycans during degeneration, and increase to approximately 500 mOsm during high mechanical loading (Sadowska et al., 2018). A reduction in tissue osmolarity leads to cell swelling (up to 20%) by the solubility-diffusion water transport across the cell membrane (Sadowska et al., 2018; Snuggs et al., 2021). Ample publications have shown that this hypoosmotic shift can activate and/or interplay with pro-inflammatory factors and catabolic responses and, hence, promote IVD inflammation and degeneration (Chen et al., 2002; Wuertz et al., 2007; Walter et al., 2016; Sadowska et al., 2020). Although the underlying mechanisms have not yet been identified, Transient Receptor Potential (TRP) channels and aquaporins may be involved (Sadowska et al., 2018; Sadowska et al., 2020; Snuggs et al., 2021). The hyper-osmotic shift in the IVD microenvironment leads to activation of the robustly expressed osmo-sensitive transcription factor TonEBP (tonicity-responsive enhancer binding protein) (Sadowska et al., 2018; Baumgartner et al., 2021c), which protects IVD cell viability under hyperosmotic stress (Tsai et al., 2006; Choi et al., 2018) and can also be regulated by cytokines (Johnson et al., 2017). These studies on IVD cells highlight that hyperosmolarity is likely not a main contributor to the catabolic shift in the IVD, whereas hypo-osmolarity seems to have detrimental effects on IVD cells. However, although no studies have specifically investigated the effect of osmolarity on IVD-infiltrating immune cells yet, research on other tissues (e.g., renal medulla, skin, lung epithelium) indicates that increased osmolarity activated macrophage inflammatory responses, which is at least partially TonEBP-dependent (with a threshold at approximately 360–380 mOsm) (Aramburu and López-Rodríguez, 2019). More research will hence be needed to better understand the role of the IVD microenvironment on macrophage polarization and the behaviour of other infiltrating immune cells.

Unbalanced Biomechanics and Mechanobiology in Catabolism

Biomechanics is another key contributor in the shift from anabolism to catabolism (Adams et al., 2000). The IVD experiences various forces throughout everyday life, which are necessary to maintain the health of the disc. For instance, the average intradiscal pressure in a healthy IVD ranges from 0.1 MPa (lying prone) to 0.5 MPa (standing flexed forward) (Wilke et al., 1999). However, damage occurs in the disc when it encounters abnormal or excessive forces, leading to catabolism, including increased cytokine production, and matrix degradation (Walter et al., 2011; Vergroesen et al., 2015; Fearing et al., 2018). This damage is believed to cause microinjuries within the disc, which gradually build up over time (Adams and Roughley, 2006; Baumgartner et al., 2021c), and is likely to contribute to the infiltration of immune cells because of the chemo-attraction effect of the pro-inflammatory cytokines released by the native IVD native cells. As the tissue degenerates, the size and composition of the IVD changes, leading to impaired response to any mechanical loading placed on the disc and causing further damage, possibly leading to disc herniation or endplate defects.

As a highly hydrated tissue, the NP provides protection to compressive forces imposed on the IVD (Adams and Roughley, 2006) while the more fibrous, surrounding AF confines the NP swelling pressure and helps the IVD to resist shear and tensile forces (Chu et al., 2018). When the NP loses hydration, the compressive load is transferred to the AF (Adams et al., 1996). Whereas the healthy AF, as a whole, is highly resistant mechanically, aberrant loading can further contribute to fissure formation where the tissue is already weakened by altered turnover, which allows for associated blood vessel growth and immune cell infiltration as discussed earlier (Lama et al., 2018).

At the cell level, specific biomechanical cues have been shown to impair IVD cell response. For example, shear stress has been found to lead to increased nitric oxide, causing downstream reduction in proteoglycan synthesis and increased apoptosis in IVD cells (Liu et al., 2001). Interactions with biochemical signalling was further demonstrated. For instance, it was found that AF and NP cells from a degenerated IVD respond differently to those from a healthy disc, suggesting that mechano-transduction pathways are altered through degeneration (Le Maitre et al., 2008; Le Maitre et al., 2009a; Chu et al., 2018) and can be modulated by cytokines such as IL-1 and IL-4 (Elfervig et al., 2001a; Gilbert et al., 2011). Additionally, a pro-inflammatory environment has been shown to change the mechanobiology of IVD cells. Treatment with inflammatory stimuli, specifically liposaccharide (LPS) or TNF-α, before osmotic loading was shown to increase hydraulic permeability and cell size, disrupt the F-actin cytoskeleton, and increase aquaporin-1, which is a main water channel in NP cells (Maidhof et al., 2014). Recently, Hernandez et al. demonstrated that inhibiting actomyosin contractility in NP cells caused a similar response as TNF-α induced inflammation, while increasing contractility protected the cells against TNF-α. Actomyosin contractility was also shown to regulate nuclear factor kappa-B (NF-κB) and downstream ECM degradation, conveying that mechano-transduction and inflammatory pathways are connected and the cross-talk could play an important role in IDD (Hernandez et al., 2020b). Thus, altered biomechanics can lead to mechanobiology alterations promoting matrix degradation and impacting the capacity of the disc to sense loads normally, leading to increased catabolism and IDD development.

Genetics and Epigenetics in IDD

Among the different causes for IDD, genetic susceptibility plays a crucial role. High heritability (over 70%) has been systematically reported for IDD (Battié and Videman, 2006; Kepler et al., 2013), as well as specific traits such as herniation (Sambrook et al., 1999) and endplate defects (Munir et al., 2018) and their progression (Williams et al., 2011b). Genetic burden, in such elevated polygenicity presented by IDD, is suggested to carry a larger effect than environmental factors, with the exception of body mass index (BMI). However, BMI itself has high heritability and polygenicity (Robinson et al., 2017), which is partially overlapped with IDD (Zhou et al., 2018).

Genetic associations for IDD have been mostly researched with candidate gene studies [see focalized reviews in Mayer et al. (2013), Feng et al. (2016), Kawaguchi (2018)]. The most representative functional group consists of genes of structural proteins and those regulating its turnover (Table 1).

TABLE 1

Structural/Regulatory componentFunctionGeneGene variantMolecular levelContribution to IDDReferernce
Collagen IXCartilage anabolic markerCOL9A2rs7533552-Associated with greater disc bulging (L1-L4)Näkki et al. (2011)
May contribute to reduced collagen crosslinkingMay contribute to disc instability and eventually prolapse in the elderlyWrocklage et al. (2000)
COL9A3Trp3 allele in IL1B 3954 C/T variantCOL9A3 gene on IDD might be modified by the IL 1β gene polymorphismInfluence MRI signal intensity in NP in the absence of the IL1β 3954 C/T alleleSolovieva et al. (2006)
Collagen XIAnabolic markerCOL11A1rs1676486Lower COL11A1 expressionHigh risk of herniationMio et al. (2007)
COL11A2rs2076311-Association with (i) disc signal intensity (ii) disc bulgingVideman et al. (2009)
Collagen IAF anabolic markerCOL1A1rs1800012-Not associated with IDD (taken as a single factor)Anjankar et al. (2015)
-Risk factor related to IDD in older peoplePluijm et al. (2004)
-Strong association with LDD in young maleTilkeridis et al. (2005)
AggrecanIVD anabolic markerACANACAN VNTR polymorphisms-Increased risk of LDD of shorter allelesKawaguchi et al. (1999), Gu et al. (2013)
-Aggrecan allele with 26 repeats is associated with dark NP MRI intensitySolovieva et al. (2007)
Cartilage Intermediate Layer ProteinCartilage-like catabolic markerCILPrs2073711TGF-β1 inhibition mediated induction of ECM proteins through direct interaction with TGF-β1Association between IDD and CILP rs2073711 variant in womenKelempisioti et al. (2011)
1184T/C-The CILP SNP 1184T/C is a risk factor for male collegiate athletesMin et al. (2010)
-Upregulation of CILP in intervertebral discs increased disc degeneration progressedSeki et al. (2005)
MetalloproteinaseCatabolic markerMMP3Combination of the T-C haplotype of IL 1α and the MMP3 minor 5A alleleIL1 promoted cartilage degradation through the induction of the matrix-degrading enzymes such as MMP1, MMP3, and MMP13Association between a combination of IL 1 and MMP3 gene variations and type II Modic changes among middle-aged Finnish menKarppinen et al. (2008)
promoter 5A/6AEnhanced the degeneration of IVD associated with environmental conditions resulting from the induction of a higher level of MMP3 expression in response to such conditionsaccelerate IVD degeneration in the elderlyTakahashi et al. (2001)
Intron 4 C/T-Associated with radiographic progression of LDDValdes et al. (2005)
MMP21306C/T-Correlation with more severe grades of disc degeneration and thus may be a genetic risk factor related to LDD susceptibility in the young adult populationDong et al. (2007)
MMP91562 C/T-Associated with a high risk of degenerative disc disease in the young adult population in North ChinaSun et al. (2009)
InterleukinCatabolic markerIL 1β3954 C/TCOL9A3 gene polymorphism on IDD might be modified by the IL-1β gene polymorphismAssociation between collagen gene polymorphisms and disc degeneration of the lumbar spine is modified or negatively confounded by the IL1β (C3954-T) polymorphism in middle-aged working menSolovieva et al. (2006)
IL 1αCombination of the T-C haplotype of IL1A and the MMP3 minor 5A alleleIL1 promotedcartilage degradation through the induction of matrix-degrading enzymes such as MMP1, MMP3, and MMP13Association between a combination of IL1 and MMP3 gene variations and type II Modic changes among middle-aged Finnish menKarppinen et al. (2008)
889C/T-IL1 gene cluster polymorphisms have an effect on the risk of disc degeneration, particularly TT genotype of the IL-1α gene promotes higher risk of disc bulgesSolovieva et al. (2004)
Significantly increased the transcriptional activity of the IL1A gene and IL-1β proteinIL 1α −889T represented a significant risk factor for the IDD-phenotypeVirtanen et al. (2007)
IL 6rs1800797, rs1800796 and rs1800795-IL 6 variants are associated with moderate IDD in a sample population of young adultsKelempisioti et al. (2011)
597G/A, 174G/C and 15T/A-association analysis provided support for a link between the IL 6 sequence variants and IDDNoponen-Hietala et al. (2005)
IL 18rs1420100-Association with severe degenerationOmair et al. (2013)
ThrombospondinECM regulationTHBS2rs9406328lower affinity for MMP binding and thus reduces MMP degradationRegulation of Intervertebral disc ECM metabolism by the THBS2-MMP system plays an essential role in the etiology and pathogenesis of lumbar disc herniationHirose et al. (2008)
A disintegrin and metalloproteinase with thrombospondin motifs (ADAMTS)Catabolic markerADAMTS5rs151058, rs229052, and rs162502Decreased binding affinity with LRP1 (protein that regulates its degradation by endocytosis)Genetic polymorphisms of ADAMTS 5 may be associated with susceptibility to LDDWu et al. (2014)
Growth differentiation factor 5Pro-chondrogenic factorsGDF5rs143383-5 population cohorts from Northern Europe indicate that a variant in the GDF5 gene is a risk factor for LDD in womenWilliams et al. (2011a)
tSKTKIAA121711 KIAA1217 variants in Exon 2, 3, 6, 7, 13, 14, 17 and 19-Strong causative candidates for the Vertebral Malformation phenotypesAl Dhaheri et al. (2020)
rs16924573-Association with lumbar disc herniationKarasugi et al. (2009)
-Association with lumbar disc herniationKelempisioti et al. (2011)
FAS receptor and ligandCell apoptosis factorsFAS and FASLrs2234767(FAS) and rs763110(FASL)-FAS and FASL may be associated with the presence and severity of LDDZhu et al. (2011)
Caspase-9CASP-91263A/G-Risk factors in the incidence of LBP in Chinese male soldiersMu et al. (2013)
rs1052576-Associated with lumbar disc herniation and disc degeneration in the Han population of northern ChinaSun et al. (2011)
Tumor necrosis factor related apoptosis-inducing ligandTRAIL1525 G/A and 1595 C/T-Associated with the susceptibility and severity of LDD in the Chinese Han populationDu et al. (2015)
Death receptor 4DR4rs4871857-Associated with the risk and severity of LDD in the Chinese Han populationTan et al. (2012)

Gene variants of structural/regulatory components of IVD associated with IDD by candidate gene approach.

Among structural components, collagen variants have been extensively assessed, and several collagen types and variants have been associated with IDD. Collagen type IX polymorphisms in alpha 2 and 3 (COL9A2, COL9A3) chains have been found to influence MRI signal intensity in NP (Wrocklage et al., 2000; Solovieva et al., 2006; Zhang et al., 2008; Näkki et al., 2011). Further, Solovieva et al. (2006) stated that the effect of the Trp3 allele in COL9A3 is dependent on an IL-1B polymorphism, reflecting the effect of immune-modulators/catabolic factors on IVD degeneration. Nevertheless, the pathophysiology of this interaction is not yet described. Polymorphisms of collagen type XI (COL11) have been associated with higher risk of herniation (Mio et al., 2007; Videman et al., 2009) and other degenerative traits (Noponen-Hietala et al., 2003; Solovieva et al., 2006; Videman et al., 2009; Kalb et al., 2012). Mio et al. (2007) stated that the variant rs1676486, which falls in cis elements region lowers COL11A1 expression due to decreased stability of its transcripts/mRNAs.

Similarly, a polymorphism found in an intronic region of the collagen 1 gene (COL1A1) that corresponds to a binding site of Specificity protein 1 (Sp1) has been shown to increase the risk of IDD, but the mechanism is not reported yet (Pluijm et al., 2004; Tilkeridis et al., 2005; Anjankar et al., 2015). However, it has been demonstrated that Sp1 downregulates pro-inflammatory cytokine-induced catabolic gene expression in disc cells (Xu et al., 2016). Additionally, Sp1-dependent mechanisms have been reported to modulate mechanically-induced apoptosis and autophagy in IDD (Li et al., 2020). Nonetheless, Sp1 also affects processes in other tissues including differentiation, angiogenesis and chromatin remodeling (Tan and Khachigian, 2009), but its potential effects on disc cells remains yet to be identified. Interestingly, Sp1 expression is inhibited by NF-κB (Tapias et al., 2008), which has been shown to be able to initiate a pro-inflammatory cascade in other tissues, often as a reaction to extracellular stimuli (Hunter and De Plaen, 2014). And although it is possible that this dual regulation confounds the potential catabolic effect of Sp1 deprivation with an inflammatory-like response, further studies are needed to distinguish the catabolic and inflammatory responses.

Another interesting gene, Aggrecan gene (ACAN), presents tandem repeat polymorphisms in the CS1 domain. Several studies have reported that lower repeat number can lead to lower chondroitin sulfate (CS), thus linking aggrecan with IVD degeneration (Kawaguchi et al., 1999; Solovieva et al., 2007; Kim et al., 2010; Gu et al., 2013). Low CS reduces the amount of water accumulated to withstand compression loadings, reducing disc mechanical properties. Further, it is possible that lower ACAN and CS reduce the IVD’s capability to recover after acute catabolic processes (Kim et al., 2010).

Another protein, cartilage intermediate layer protein (CILP), whose expression is restricted to different cartilage(-like) tissues including IVD, inhibits TGFβ, therefore preventing the ECM anabolism and cell proliferation promoted by TGFβ in IVD (Liu et al., 2021). A polymorphism in the interaction region between CILP and TGFβ has been shown to change their binding affinity, consequently, identifying it as a risk factor for IVD degeneration (Seki et al., 2005; Min et al., 2010; Kelempisioti et al., 2011). Additionally, CILP is able to inhibit Insulin-like growth factor-1 receptor (IGFR1), acting as an antagonist of Insulin-like growth factor 1 (IGF1), a factor that mediates chondrocyte anabolism and proliferation (Liu et al., 2021). Similarly, tandem repeat polymorphisms in Asporin gene inhibits TGFβ-induced anabolism with likely synergic effects with the CILP variant (Song et al., 2008; Min et al., 2010).

As stated before, pro-inflammatory cytokines are the key factors that start the catabolic shift through increase of matrix-degrading enzymes expression, with IL-1β being one of the most influential cytokines produced by the native IVD cells and immune cells following IVD rupture (Le Maitre et al., 2005; Millward-Sadler et al., 2009; Phillips et al., 2015). As mentioned earlier, combinations of specific IL-1 β and COL9A3 polymorphisms constitute a risk factor for IVD degeneration. Similarly, a combination of MMP3 and IL-1 β polymorphisms also presents a higher risk of IVD degeneration (Karppinen et al., 2008). Additionally, other MMP-3 (Takahashi et al., 2001; Valdes et al., 2005), MMP-2 (Dong et al., 2007) and MMP-9 (Sun et al., 2009) polymorphisms have shown greater risk. Likewise, different single nucleotide polymorphisms (SNPs) of IL-1 β and their combinations present increased risk of IDD, possibly due to overactivation under mechanical stress (Solovieva et al., 2004; Virtanen et al., 2007; Karppinen et al., 2008). Other interleukins such as IL-6 (Noponen-Hietala et al., 2005; Kelempisioti et al., 2011) and IL-18 (Omair et al., 2013) have also been reported to increase catabolic processes in the IVD, similar to IL-1 β.

Matrix-degrading enzymes can also be regulated through their degradation. A SNP in the Thrombospondin-2 gene (THBS2), which regulates degradation of MMPs through endocytosis, has shown lower affinity for MMP binding which in turn reduces MMP degradation and increases the risk of IDD (Hirose et al., 2008). In similar fashion, different variants of diverse ADAMTS family members, which degrade aggrecan, were identified to increase risk of degeneration (Rajasekaran et al., 2014; Wu et al., 2014; Liu et al., 2016) and are increased during IVD degeneration (Pockert et al., 2009). A SNP in ADAMTS identified by Wu et al. (2014) decreases binding affinity with the protein that regulates its degradation by endocytosis, LRP1, therefore increasing its catabolic activity.

A shift toward catabolic processes can also be provoked by perturbing cell differentiation and viability. This has identified different variants of growth factors to be associated with IDD, such as Growth Differentiation Factor 5 (GDF5) (Williams et al., 2011a), which is a key regulator of matrix synthesis in the disc (Le Maitre et al., 2009b); SKT gene (KIAA1217) (Karasugi et al., 2009; Kelempisioti et al., 2011; Al Dhaheri et al., 2020), vascular endothelial growth factor (VEGF), and endothelial nitric oxide synthase (eNOS) (Han et al., 2013). Variants in factors that modulate cell apoptosis, such as FAS receptor and its ligand (FASL) (Zhu et al., 2011), Caspase-9 (Sun et al., 2011; Mu et al., 2013), tumor necrosis factor related apoptosis-inducing ligand (TRAIL) (Du et al., 2015), and Death receptor 4 (DR4) (Tan et al., 2012) have been found to be associated with IDD. However, exact mechanisms of how such variants affect cell fate and IDD are still unclear and require further investigation.

In addition to the candidate gene approaches, where the set of genes tested is preselected, a few Genome-wide Association Studies (GWAS) have been performed. GWAS is an agnostic method that tests variants covering the “whole” genome (Duncan et al., 2019). A GWAS performed by FMK Williams et al. found a variant of Parkinson protein 2, E3 ubiquitin protein ligase (PARK2) and one of Proteasome 20S Subunit Beta 9 (PSMB9) to be associated with degenerative discs (Williams et al., 2013). Those genes encode for proteins that aim to tag and degrade unwanted proteins, providing another method to remove matrix-degrading enzymes, altering the metabolic balance. Another GWAS identified a variant of Carbohydrate Sulfo-Transferase 3 (CHST3), a catabolic enzyme that catalyses proteoglycan sulfation, as a susceptibility gene for IDD (Song et al., 2013). The authors suggest that this enzyme interacts with a micro-RNA (miRNA) that targets proteins with important regulatory functions in cell-mediated immune responses, but further analysis is needed to confirm such hypothesis.

Bacterial Infection and Disc Microbiome in IDD

Bacterial contamination has also been proposed as an important regulator of disc cell inflammation and catabolism, particularly in association to Modic changes (Gorth et al., 2015; Granville Smith et al., 2021). Despite the detection of various bacteria within isolated disc tissue, the presence of an IVD microbiome is still controversial as it has been traditionally considered as a sterile, immune privileged structure (Rajasekaran et al., 2020).

In Stirling et al. (2001), first reported the presence of anaerobic bacteria, particularly of Cutibacterium acnes (Gilchrist, 1900), previously known as Propionibacterium acnes, within the disc tissue of 43 out of 140 patients with sciatica. Recently, Granville Smith et al. (2021) performed a PRISMA systematic review identifying 36 articles from 34 research studies investigating bacteria in human IVDs. Bacteria were identified in 27 studies, whereas nine attributed bacterial presence to contamination. C. acnes, a Gram-positive anaerobic bacterium that is part of the natural skin microbiome, was the most abundant. Which is also associated with prosthetic joint infection (Jauregui et al., 2021) and was shown to be able to interact with bone cells (Aubin et al., 2017) and recently disc cells (Capoor et al., 2021). Coagulase-negative (CoNS) bacteria of the genus Staphylococcus Rosenbach 1884 were the second most abundant (Granville Smith et al., 2021). Inconsistencies between the identified bacteria and the prevalence of different bacteria can be partly administered to differences in tissue source (intact or herniated tissues), culture conditions (anaerobic vs. aerobic, culture time, culture media), differences in the methods used to detect bacteria, and differences in the administration of antibiotics. To date, there are few quantitative studies investigating bacterial infection to show whether bacteria are present in vivo or represent operative contamination.

Albert et al. (2013a) hypothesized that type 1 Modic changes in the adjacent vertebrae of herniated discs may be due to infection of the disc, highlighting the need to investigate bacteria presence in the disc. Treatment of chronic LBP and Modic changes with antibiotics has generated great controversy. It has been shown that in a certain subset of patients, antibiotic treatment was effective to reduce pain as well as disability (Albert et al., 2013a; Gilligan et al., 2021). However, this result has failed to be replicated in subsequent studies (Bråten et al., 2020).

Furthermore, the role of potential bacteria within the disc is unknown. However, previous studies have shown that LPS, a main component of Gram-negative bacteria, induces upregulation and production of various proinflammatory cytokines and matrix degrading enzymes in the NP over-activation of the NF-κB pathway (Li et al., 2016). However, most bacteria detected in the disc to date are Gram-positive and the potential influence of those bacteria on the disc remains poorly understood. Recently, Capoor et al. (2021) stimulated human IVD cells with C. acnes demonstrating an induction of catabolic cytokine expression by native NP cells, suggesting that at least in some individuals the increased catabolic cytokines during disc degeneration could be triggered by bacterial infection. Further work is required to understand whether bacteria are present within the disc and whether bacteria could act as a trigger to the catabolic stimuli seen during disc degeneration, and whether the gut microbiome could influence disc degeneration (Li W. et al., 2022).

Translating Knowledge of Initiating Factors of IDD to Clinics

Understanding the roles and interactions of each of these initiating factors is essential in order to diagnose IDD early and identify suitable treatments. Current treatments alleviate pain but do not regenerate the disc, therefore regenerative strategies are urgently needed in clinics (van Uden et al., 2017).

Despite extensive research, tissue engineering strategies have had limited success in translating from preclinical models to beneficial treatments in patients as they fail to address pain (Isa et al., 2022). Precision medicine appears promising to use multiomics profiling to elucidate the pathology of IDD in each patient and prescribe an individualized treatment plan and determine which therapies could be most effective (Isa et al., 2022). However, a more comprehensive understanding of the disc microenvironment, biomechanics, genetics and epigenetics, and even bacterial infection throughout the different stages of IDD and low back pain would be necessary to evaluate which therapies or combination of therapies could be effective in a patient. In particular, there is a lack of studies on cartilage endplate regeneration therapies, despite its importance in the nutrient supply of the IVD (van Uden et al., 2017). Overall, experimental and computational models of IDD, which will be explored in the next chapter, remain critical toward developing novel treatments and regenerative therapies for IDD.

Research Methods for Exploring the Inflammatory or Catabolic Environment of the IVD

Experimental Models in IVD Research

For many years different experimental models have been developed for IVD research to mimic IDD. Several approaches have been used to replicate the physiological state of the IVD as closely as possible, including 3D cell and organ culture models, bioreactors and animal studies (Figure 4B). These approaches are crucial in elucidating the causes and progression of IDD, as well as in developing and testing novel therapies. Nevertheless, the best strategy to investigate IDD remains unclear, and each culture system or animal model offers different advantages and disadvantages that should be considered when planning an experiment (Table 2).

FIGURE 4

TABLE 2

ModelCulture system or methodsRepresentative studiesScaleParameters that can be probedAdvantagesDisadvantagesContribution to IDD
Experimental3D cell cultureAlginate-based hydrogelsLe Maitre et al. (2005), Hernandez et al. (2020a), Naqvi and Buckley, (2014), Le Maitre et al. (2008), Zhang et al. (2014), Wang et al. (2014), Sun et al. (2015), Öztürk et al. (2016), Thorpe et al. (2018), Lazarus et al. (2021), Rosiak et al. (2021)CellDynamic loading and catabolism induction, Stimulation with cytokines and catabolic factorsInexpensive, non-toxic and excellent cell phenotype maintenance, human cells are availableNot recommended for AF cell cultureIVD phenotype maintenance and recovery from catabolism
Agarose carriersKelly et al. (2009), Smith et al. (2011)Dynamic loading and catabolism inductionlong-term 3D culture, human cells are availableCannot replicate macroscale forces the IVD experiencesCatabolism induction in IVD cells
Reinforced hydrogels (Silk)Frauchiger et al. (2017), Wöltje and Böbel, (2017)Porosity, coating and surface area of scaffoldBiodegradable and resorbable biomaterial with high cytocompatibilityThe outer layer can cause immune responseNovel therapeutic approach in IVD repair
Pellet culture systemsThorpe et al. (2018), Wangler et al. (2019), Hingert et al. (2019)Hydrostatic loading and nutrient perturbationSimple, inexpensive, human cells are availableChondrogenic phenotype induction and lack of ECMEffects of hydrostatic loading and nutrient deprivation in IDD
Ex vivoBioreactors for mechanical loadingGantenbein et al. (2006), Vergroesen et al., 2014; Costi et al., 2008, Walter et al. (2011), Paul et al. (2012), Chan et al. (2013), Paul et al. (2013), Salzer et al. (2021), Croft et al. (2021)Organ/TissueMechanical loading (compression, torsion, bending, flexion, extension, and asymmetric), environmental control (nutrition, pH, temperature, oxygen level), frequency and duration (static, dynamic, or diurnal)Mimics physiological conditions, possibility for automation, reproducible, in line with 3R principles (“Replacement, Reduction, and refinement")Expensive and difficult to build, culture time limited to ∼1 month, no connection to vasculature or immune system, most cannot test large sample sizes at once, human tissue is limitedDetermined physiological and catabolic ranges of mechanical loading regimens and test the mechanical properties of suitable biomaterials for IVD replacement
Perfusion bioreactors, microfluidics, “disc-on-a-chip"Elfervig et al. (2001b), Chou et al. (2016), de Bournonville et al. (2019), Dai et al. (2019), Hwang et al. (2020), Kim et al. (2021), Mainardi et al. (2021)Diffusion, shear stress, fluidic pattern, electrical impulses, environmental control (nutrition, pH, temperature, oxygen level)Possibility for automation, reproducible, extends culture time, in line with 3R principles (“Replacement, Reduction, and refinement")Cannot replicate macroscale forces the IVD experiences, difficult to design complex systemsIncreased cell viability in culture provided platform to investigate cellular response to shear stress and interactions with inflammatory and neurotrophic factors, and test treatments such as electrical stimulation
In vivoSpontaneous degenerationDaly et al. (2016)Whole body/Organ/tissue/Protein/GeneticDegeneration progression, therapies/treatmentsOccurs naturally (therefore more ethical), immune system and pain responseLong and unpredictable time course, inherent biological and biochemical differences to humans, expensive and complex, ethical considerationsChondrodystrophic dogs and sand rats have similar pathological changes to human in IDD and have useful in testing cellular therapies and other clinical treatments
Altered mechanical loadingLindblom (1957), Iatridis et al. (1999), Phillips et al. (2002), Ching et al. (2003)Magnitude, duration, and frequency of loadingImmune system and pain response, repeatable induction of IDD at specific time pointInherent biological and biochemical differences to humans, expensive and complex, ethical considerationsModels of bending, compression, and spinal fusion have shown that loading changes the mechanical properties of the IVD
Structural models (physical injury or chemical injection)Lü et al. (1997), Holm et al. (2004), Ulrich et al. (2007), Yoon et al. (2008), Dudli et al. (2011), Dudli et al. (2014), Alkhatib et al. (2014), Zhang et al. (2020), Zhou et al. (2021)Degeneration progression, proteoglycan degradation, therapies/treatmentsImmune system and pain response, repeatable induction of IDD at specific time point, useful in preclinical trialsInherent biological and biochemical differences to humans, expensive and complex, ethical considerations, fail to capture pathogenesis of human IDD, viability of native IVD cells preservedCritical in understanding IDD and developing/testing novel therapies for clinical application
Transgenic modelsMillecamps et al. (2011), Millecamps et al., 2012; , Li et al. (2009), Wang et al. (2012), phillips et al., 2013, Phillips et al., (2013a,b)Miyagi et al. (2014), Gorth et al. (2018)Genes and gene pathwaysImmune system and pain response, target specific pathways of interestInherent biological and biochemical differences to humansChanges in SPARC, Tg197, CCN2, IL-1rn, cAct, and SMAD3 genes have been identified to contribute to IDD
ComputationalFinite elementJones and Wilcox, (2008), Newell et al. (2017), Ghezelbash et al. (2020), Galbusera et al. (2011b), Malandrino et al. (2011), Volz et al. (2022)Organ/TissueLoading, environmental perturbations and catabolism inductionValuable prediction of altered mechanics and transport at the tissue level of the IVD, 3RsChallenging comprehensive validation and the cellular and sub-cellular level is not contemplated. Computationally intensive and time consumingPredictions about metabolic rates, oxygen and lactate transport, osmotic behaviour
Agent-based2D or 3DBaumgartner et al. (2021a), Baumgartner et al. (2021b)Tissue/CellularCell behavior and interactions, microenvironment, time, Cell types, environmental perturbationsDynamic, ability to model heterogenous populations, flexible, stochastic, reveals emergent phenomenaCan be computationally intensive, only as good as the rules inputtedVisual predictions of NP cells expressing TNF-α, IL-1β, or both TNF-α & IL-1β
NetworkKnowledge based or data drivenShannon et al. (2003), Hu et al. (2004), Batagelj and Andrej, (2004), Kashtan et al. (2004), Yu et al. (2004), Schreiber and Schwöbbermeyer (2005), Wernicke and Rasche (2006), Milenković et al. (2008), Baumgartner et al. (2021a), Baumgartner et al., 2021b, Xu et al. (2021)Protein/GeneticProtein-protein interaction (PPI) and transcriptomic/proteomic analysisIntuitive way to investigate, characterize, and understand interactions between biological components under micro-environmental stimuliDifficult to construct from available IVD data which has high sample variation, different stages of IVD and type of disc tissue, and the varying methods of analysisCapture interactions between transcriptome, proteins and their pathways in to understand the critical biochemical factors in IVD regulation. Reveal complex dynamics behind unbalanced metabolism
Genetic analysisCandidate gene studies, GWASWilliams et al. (2013), Song et al. (2013), Feng et al. (2016), Kawaguchi (2018), Mayer et al. (2013)GeneticGenes and gene pathwaysEffective in identifying genes implicated in IDDIncapable of explaining high heritability in complex diseases due to high polygenicity and unmet “common disease, common variants” hypothesis, and due to other heritable properties as epigeneticsIdentified large number of risk genomic loci involved in IDD (Table 1)
Machine learning/AI/Deep learningClassification of discsRim (2016), Niemeyer et al. (2021)Whole body/Organ/Tissue/Protein/GeneticImaging, clinical categories, compound structures, gene sequence, protein/RNA dataLink seemingly unrelated entities of complex/diverse biological dataNeed for algorithm creation and learning. Very subjective score systemDeep learning model for the classification of discs based on MRI with an average sensitivity of 90%
Simplifying or coupling complex modelsPfaff et al. (2020)Highly accurate surrogate models, significantly less computational resources and less time-consuming

Experimental and computational approaches applied for IDD research.

3D Cell Culture Systems

Over the last decades, three-dimensional (3D) cell culture models have been widely accepted due to the considerable improvements they possess in comparison to two-dimensional (2D) culture, including improved phenotypic retention to that seen in vivo, including: cell shape preservation; proliferation rates, and gene and protein phenotypic marker and matrix expression (Jensen and Teng, 2020). Conventional 2D monolayer culture systems lack the spatial architecture of the tissue, inducing a loss of cell phenotype and cell-ECM interactions. In contrast, 3D cultures environments promote extracellular matrix (ECM) deposition, a key factor for the maintenance of NP cell phenotype (Guerrero et al., 2020). Likewise, previous studies using notochordal cells have reported the negative effects of 2D culture (Rastogi et al., 2009) and the necessity of 3D culture system, preferably in hypoxia, and raised osmolality to maintain the phenotype (Gantenbein et al., 2014).

In the last decade many different biomaterials have been used in 3D cell culture of IVD cells, which will be listed according to the best outcome. Alginate-based hydrogels are commonly used because they are inexpensive, non-toxic and an easy 3D hydrogel model whilst maintaining excellent cell phenotype (Hernandez et al., 2020a). Notably, previous studies have reported IVD phenotype maintenance and recovery from catabolism after 3D alginate culture (Le Maitre et al., 2005; Naqvi and Buckley, 2014; Wang et al., 2014; Zhang et al., 2014; Sun et al., 2015). Additionally, alginate has been used with dynamic loading systems (Le Maitre et al., 2008) and for inducing catabolism (Le Maitre et al., 2005; Le Maitre et al., 2008). Moreover, modified alginates have shown novel properties and applications in biomedical research (Rosiak et al., 2021). For example, the sulfation of alginate hydrogel has been reported to preserve the phenotype of chondrocytes (Öztürk et al., 2016; Lazarus et al., 2021). Thus, alginate-based hydrogel systems are considered as promising biological constructs for NP cell culture (Thorpe et al., 2018). In contrast, the consistency of alginate-based hydrogels is not recommended for AF cell culture due to the lack of fibrotic structure. In addition, other materials such as agarose have been reported as long-term 3D culture models for inducing catabolism in IVD cells (Smith et al., 2011) and chondrocytes (Kelly et al., 2009). In contrast, other 3D hydrogel models such as fibrin-clots are not a realistic option for IVD 3D culture due to their lack of stiffness, despite the easy modification of this natural hydrogel. However, reinforced hydrogels with resorbable biomaterials, for instance silk, are a novel therapeutic approach in IVD repair due to high cytocompatibility (Frauchiger et al., 2017). Notably, the outer layer of silk filaments can cause immune response due to the presence of sericin (Wöltje and Böbel, 2017). 3D environment features can also be achieved without biomaterials through pellet culture. Notably, pellet culture systems have been utilized to investigate IVD degeneration, specifically examining the effects of hydrostatic loading and nutrient deprivation (Hingert et al., 2019; Wangler et al., 2019), however they induce a more chondrogenic phenotype rather than NP phenotype (Thorpe et al., 2018). Furthermore, the lack of ECM after the formation of the pellets could influence the response of hydrostatic loading (Zeiter et al., 2009), and they fail to mimic the cell-ECM connections and cell density seen in the IVD, therefore other materials would be better utilised.

In contrast, 3D constructs aim to offer a more physiological interaction between cells with ECM components without the presence of vasculature and innervation as well as the interaction with the immune system. Nevertheless, although biomaterial-based therapies have been developed in the last decades to prevent IDD, only a small number of bioengineered therapies are currently undergoing clinical trials (NCT02338271, NCT01290367, NCT01290367 and NCT02412735). Notably, pain relief does not correlate adequately with functional and structural IVD restoration. Overall, the main clinical challenge is to use clinical signs, patient pain, and disability history alongside advanced imaging techniques to design a sufficient biomaterial approach (Huang et al., 2018; Isa et al., 2022). However, ex vivo culture systems seem to be an appealing alternative resulting in a more representative model. Particularly, IVD explants gained attention due to a higher control of the degeneration state, sample geometry and loading (Salzer et al., 2021). Similarly, organ culture bioreactor models allow the presence of the native tissue microenvironment together with a loading system promoting a bridge between in-vitro and in-vivo models (Gantenbein et al., 2015). Such ex vivo culture systems are good model systems investigating intact IVDs. However, they do not enable connection with the vasculature and immune system, although co-culture systems could be developed to model these interactions.

Bioreactors and Microfluidic Devices

Bioreactors are widely accepted as pre-clinically relevant devices that simulate the microenvironment and offer a platform to evaluate the effects of limited nutrition and incorporate more complex parameters, such as mechanical loading and fluid flow, into in vitro and ex vivo experiments (Haglund et al., 2011; Illien-Jünger et al., 2014; Walter et al., 2014; Gantenbein et al., 2015; Gantenbein et al., 2019; Pfannkuche et al., 2020).

Various materials or organ culture are used for different bioreactor systems. Freshly isolated IVDs from bovine tails are often used in organ culture studies due to well established operating procedures, as well as their similarities to human discs (Chan and Gantenbein-Ritter, 2012; Saravi et al., 2021). Bioreactors are also used in dynamic 3D cell culture, however the material used must be able to withstand the imposed forces. Consequently, cells are often seeded into hydrogels, such as agarose or alginate, to offer more protection against mechanical loading (Fernando et al., 2011; Cambria et al., 2020).

Initial approaches were only capable of static loading; however, over time, bioreactors have gradually evolved to become more complex, incorporating diurnal loading (Gantenbein et al., 2006) and dynamic compression (Paul et al., 2012; Paul et al., 2013; Vergroesen et al., 2014). More recently, bioreactors have advanced past compression to include two- and six-dimensional degrees-of-freedom, allowing for the analysis of torsion, bending, flexion, and extension (Costi et al., 2008; Chan et al., 2013; Croft et al., 2021). Additionally, asymmetrical complex loading has been proposed as a model to study the effects of scoliosis on disc mechanobiology (Walter et al., 2011). These increasingly complex loading devices are crucial to better understanding the effects of mechanical loading on IDD and how aberrant mechanical loading contributes to the shift to catabolism. Further, these devices are highly clinically relevant as they can test the mechanical viability of novel regenerative therapies that aim to replace or regenerate the NP, AF, and/or CEP.

In addition to mechanical loading, bioreactors are useful in simulating fluid flow to the IVD. Perfusion bioreactors have been developed to allow for in vitro perfusion culture of scaffold-based tissue engineering constructs, offering the ability to monitor and control key parameters such as temperature, pH, and fluidic pattern (de Bournonville et al., 2019). Microfluidic, or “organ-on-a-chip”, platforms have also been explored to study IDD and have been reviewed recently by Mainardi et al. (2021). In 2019, one of the first microfluidic “disc-on-a-chip” devices was developed by Dai et al. (2019), permitting continuous media flow to mimic the disc microenvironment, and demonstrating higher cell viability than cells in static culture, allowing for the possibility of long-term organ culture to examine chronic disc degeneration. Microfluidic devices have also been used to investigate mechanical loading in AF cells through fluid-induced shear stress (Chou et al., 2016). Studies have found that AF cells had a greater response to shear stress when stimulated with IL-1β, suggesting that disc cells are more sensitive to shear during catabolic or inflammatory conditions, possibly affecting IDD development (Elfervig et al., 2001b). More recently, electrical stimulation was tested as treatment to modulate IL-1β-mediated catabolism in NP cells (Kim et al., 2021). Additionally, a microfluidic platform was used in a co-culture system of AF, NP, and endothelial cells to investigate IDD development from inflammatory and neurotrophic factors, which could be further developed to examine pain mechanisms in IDD (Hwang et al., 2020; Mainardi et al., 2021). The ability to evaluate pain in culture systems of IVD is currently lacking, which is a major issue in translating tissue engineering strategies successfully to clinics (Isa et al., 2022). Therefore, development of a microfluidic platform that could do this would be an immense step forward toward evaluating new therapies and treatments in vitro. Further, microfluidic devices have been proven valuable in testing drug delivery and improving screening strategies (Damiati et al., 2018). Although these systems have been around for less time and are therefore less validated than bioreactors that offer mechanical loading, perfusion and microfluidic systems offer a promising platform to probe inflammatory and catabolic parameters and test new treatments in vitro.

In Vivo Animal Models of Disc Degeneration

While in vitro and ex vivo models of the IVD are highly beneficial and provide insights on components of IDD, they do not offer the same level of complexity as in vivo studies, which may better examine the multifactorial nature of IDD and can include an immune system and pain response (Daly et al., 2016). Many animal models have been used to investigate the IVD, and the advantages and disadvantages have been reviewed extensively before (Lotz 2004; Alini et al., 2008; Showalter et al., 2012; Daly et al., 2016; Jin et al., 2018). However, no perfect model of disc degeneration currently exists, because there are many biological and biochemical differences between discs from animal species and those from humans (Oshima et al., 1993).

One major difference between human discs and animal discs is the presence of notochordal cells in the NP. In humans, notochordal cells are present at birth, but rapidly decrease and are gone by adulthood. In most other species, notochordal cells are present in adulthood. However, cows and sheep and some species of dog, classified as “chondrodystrophic”, lose their notochordal cells rapidly, similarly to humans. Although notochordal cells are not well understood, they are often considered as progenitor cells, and therefore their presence in degenerative animal models may lead to results that are minimally relevant to understanding human LDD (Alini et al., 2008). Other biochemical parameters should also be considered, such as water, GAG, and collagen content, as well as how these factors change with degeneration and age in animals versus humans, which has been reviewed previously (O'Connell et al., 2007; Beckstein et al., 2008; Miyazaki et al., 2009; Showalter et al., 2012). It should also be noted that rodents have distinctly different aggrecan proteins and do not express the same MMPs as humans, which is considered important in catabolism and tissue remodeling (Barry et al., 1994; Flannery et al., 1998). However, while rodent models may not be suitable for translational research and testing new therapies because of the major differences to humans, they offer a useful platform to elucidate the genetic basis of IDD and catabolic changes due to aging (Masuda and Lotz, 2010; Mainardi et al., 2021). Similarly, larger animal models are not a perfect match toward human IDD (Alini et al., 2008; Gullbrand et al., 2016). Nevertheless, they are furthermore suitable for initial tests of regenerative therapies as they offer valuable information on the changes in mechanical loading, whether an immune response is initiated, and possibly whether any pain is resolved.

Spontaneous degeneration occurs in mice, sand rats, chondrodystrophic dogs, and baboons; however, these models are unpredictable and often time-consuming (Daly et al., 2016). Therefore, there are various methods, categorized under mechanical or structural, that have been used to induce degeneration in animals.

In rat tails and rabbits, degeneration has been induced through altered mechanical loading, such as bending (Lindblom 1957), compression (Iatridis et al., 1999), or spinal fusion (Phillips et al., 2002; Oswald et al., 2021). In compression, the magnitude, duration, and frequency of loading cause significant changes in IVD mechanical properties, and static loading produces greater changes than cyclic loading (Ching et al., 2003).

Structural models involve a physical injury or chemical injection to the CEP, AF, or NP (Lotz 2004). Physical injuries are done using either a drill bit, scalpel, or needle. Annular injuries are commonly used and have been shown to cause decreased disc height, higher Pfirrman degeneration scores, and decreased NP GAG content (Yoon et al., 2008). Research has also shown that repetitive injury causes different inflammatory responses in the IVD. Ulrich et al. (2007) found that while a single stab injury in a rat led to localized, short-term pro-inflammatory response, while multiple stab injuries cause a prolonged upregulation of proinflammatory cytokines TNF-α, IL1-b, and IL-8 for up to 28 days after injury. CEP injuries have also been demonstrated to lead to disc degeneration similar to that of humans, characterized by decreased NP proteoglycan content and intradiscal pressure (Holm et al., 2004; Dudli et al., 2014; Zhou et al., 2021), as well as increased catabolic enzyme production and pro-inflammatory gene expression seen following CEP fracture (Dudli et al., 2011; Alkhatib et al., 2014). However injurious degeneration models fail to recapitulate the pathogenesis of human IDD and enable infiltration of inflammatory cells at a much earlier time frame than seen if at all in humans. Chemical injections with papain and chondroitinase ABC or papain are commonly used methods to induce degeneration through degrading proteoglycans in the disc (Daly et al., 2016). Although both cause catabolism, chymopapain was shown to cause greater destruction of the NP and AF proteoglycans, as well as greater spinal instability and disc space narrowing (Lü et al., 1997). However, chondroitinase ABC induced a similar catabolic shift to that seen in human IDD in the IVD of goats (Zhang et al., 2020).

In addition, groups have used transgenic animal models to represent IDD (Jin et al., 2018). The SPARC (secrete protein, acidic, rich in cysteine)-null transgenic mouse has been shown to develop behavioral signs consistent with chronic low back pain due to IDD, such as hypersensitivity to cold, axial discomfort, and motor impairment (Millecamps et al., 2011; Millecamps et al., 2012). Further, the SPARC-null mouse showed age-dependent increased innervation by sensory nerve fibers near the IVD (Miyagi et al., 2014). Gorth et al. (2018) also used Tg197 mice, a TNF-α transgenic line, to investigate the effects of systemic over-expression of TNF-α on IDD, finding that the experienced an increase in annulus tears and herniation with higher vascularity and immune cell infiltration. However, they found that intact IVDs remained healthy despite the elevated inflammation. Additionally, knockout technology has been used to create models of notochord-specific CCN2-null mice (Bedore et al., 2013), IL-1 receptor antagonist knockout mice (Phillips et al., 2013a), and β-catenin conditional activation (cAct) mice to examine the signaling pathway roles in disc degeneration (Wang et al., 2012).

Finally, emerging strategies such as tissue-engineered replacement discs have gained substantial attention in the IVD regeneration field. In terms of animal models, significant technical challenges must be addressed including cell source, construct size, culture strategies, and translational models (Gullbrand et al., 2018). Nevertheless, several studies in disc replacement, including prospective randomized comparative trials, have demonstrate advantages such as short-term superiority to spinal fusion (Hellum et al., 2012; Vital and Boissière, 2014) or at least non-inferiority to anterior spinal interbody fusion (Blumenthal et al., 2005; McAfee et al., 2005; Vital and Boissière, 2014).

Computational Modeling of the IVD

While experimental studies are valuable in determining cell sensitivity to biochemical and mechanical cues, they are not sufficient to capture the full complexity of cell response and interactions with the microenvironment, which is crucial for understanding the transition to catabolism and the initiation of an immune response. In addition, experiments are often expensive and time-consuming. In-silico models can use published literature and experimental data to predict multifactorial tissue and ECM regulation at multiple scales (Figure 4). Further, in silico models can offer the possibility of exploring patient-specific IDD and predicting the effects and risks of available therapies prior to being treated (Rijsbergen et al., 2018). Finite element models (FEM) are useful in determining the effect of mechanical loading at the tissue and organ level, while agent-based models (ABMs) are effective in predicting tissue and cellular level changes due to the microenvironment. At the subcellular level, network modelling provides further insight into the effects of cell signaling pathways and gene variants. Machine learning and deep learning are valuable tools to analyze and classify clinical images and predict the current and future status of a patient. Each of these in-silico tools offers a novel way to explore catabolism or inflammation in IDD, which will be further explained in the following paragraphs.

Finite Element Models in IDD

Finite element models (FEM) have been extensively used to represent the intervertebral disc and to simulate structural changes due to mechanical loading, providing a deeper understanding of each component’s role than what can be tested through experiments. The IVD is inhomogeneous, anisotropic, and porous, making it a highly complex structure (Newell et al., 2017). Material properties for each of IVD component, the NP, AF and CEP, are defined and validated through experimental measurements and clinical observations, however comprehensive validation of FE analysis is challenging due to the complex structure and interactions (Ghezelbash et al., 2020). The IVD components are generally modeled as a biphasic material, with an incompressible fluid phase and an elastic solid phase. Other reviews have already been written regarding FEM studies of the IVD (Jones and Wilcox, 2008; Newell et al., 2017; Ghezelbash et al., 2020), so here we will focus on how FEM has been applied to study degeneration and catabolism.

FEMs are a valuable tool to examine the vicious cycle of disc degeneration and aberrant loading. Once a disc is degenerative, the tissue biomechanics are altered, leading to a catabolic environment and causing further damage over time. However, degeneration is not uniform throughout subjects, and there are various ways previous groups have simulated degenerative discs. These include geometrical changes such as decreased height and reduced NP area (Ghezelbash et al., 2020), but also changes in the material properties, including reduced water content, calcified and thinner CEP, and a stiffer NP characterized by a decreased bulk modulus (Galbusera et al., 2011a). These studies predicted that a degenerative IVD experiences higher forces during axial rotation, as well as lower fluid flow and recovery of intradiscal pressure after loading. Models have also shown that as the NP loses fluid, it carries less load under compression, as well as with bending and shear (Ghezelbash et al., 2020). Investigation of the geometry of the IVD has concluded that simplified geometry is less stiff and does not capture the same strain distribution as FEMs based on more complex geometry obtained through segmentation of MRIs, conveying those accurate geometries are essential (Du et al., 2021).

However, there is a lack of studies that measure the effects of different patient-specific morphologies, either to observe the mechanical effects of deformation or their implications on nutrient transport. Recently, a coupled and patient-specific mechanoregulated model was developed to predict the effects of spinal fusion on disc degeneration and bone density, demonstrating how FEMs can be used by surgeons to provide insight into which patients could possibly benefit from spinal fusion treatments (Rijsbergen et al., 2018). Similarly, future models could aim to use available clinical data to help develop models that aid doctors in predicting which treatments and surgical interventions would have the best outcome.

Many FEMs of the IVD also simulate osmotic behavior, and it has been shown that a swelling model with strain-dependent osmotic pressure most accurately represents the IVD, and could be applied to investigate crack opening and fissure propagation (Galbusera et al., 2011b). A mechano-transport FEM of the IVD developed by Ruiz Wills et al. (2018) found that CEP permeability increases with aging and degeneration, and that CEP degeneration could be a cause of NP dehydration and play a key role in IDD. Other groups have simulated cell metabolism and nutrient levels in the IVD, predicting that higher cell metabolic rates lead to nutrient depletion and that application of mechanical loading led to decreased glucose levels throughout the IVD (Volz et al., 2022). Additionally, simulations of compression on oxygen and lactate transport within the IVD suggested that degenerative changes including disc height, fluid content, nucleus pressure, and cell density reductions significantly affected transport (Malandrino et al., 2011).

Overall, FEMs have proven valuable in predicting the effects of altered mechanics and transport due to degeneration at the tissue level of the IVD. However, FEMs fail to take into account what is happening at the cellular and sub-cellular level.

Agent-Based Models

Agent-based models (ABMs) are widely used across different spatial scales and research areas. Hence, agents might reflect human beings for socioeconomic studies (Alvarez-Galvez and Suarez-Lledo, 2019) or (sub-)cellular entities in cancer research (Metzcar et al., 2019). They are particularly useful for studying complex biological processes, such as inflammation and tissue degeneration, that are dynamic, spatially heterogeneous, and stochastic. ABMs can represent individual biological cells as computational agents and can simulate how collections of cells within a tissue will respond emergently to literature-derived rules. Previously, ABMs have been shown as valuable in simulating tissue degeneration and inflammation in musculoskeletal and cardiac tissue, spanning many cell types including immune cells, fibroblasts, and stem cells (Virgilio et al., 2018; Rikard et al., 2019). Thus, ABMs offer much potential to simulate cell dynamics in cartilage tissue (Pearce et al., 2020).

In IVD research, however, ABMs have been used only recently, though initial studies have demonstrated their value in predicting IVD cell responses in a pro-inflammatory environment. Baumgartner et al. (2021c) coupled an ABM with mathematical network models (see Section “Network Modeling”) to investigate the relative mRNA expression of proteins and proteases in NP cells of different pro-inflammatory cell states, i.e., immunopositive for TNF-α, IL-1β, both, or none (Baumgartner et al., 2021a; Baumgartner et al., 2021b). The ABM was used to visualize cell states within the NP through predicting how immunopositive cells could be arranged within a 3D environment. Thus, it was assumed that immunopositive cells were organized in clusters based on experimental data of autocrine and paracrine stimulation (Phillips et al., 2013b; Phillips et al., 2015), short half-lives of cytokines according to distantly related studies (Kudo et al., 1990; Oliver et al., 1993; Larson et al., 2006) and diffusion of pro-inflammatory cytokines. While validation is still limited due to the lack of experimental information on the arrangement of immunopositive cells, the ABM presented a novel projection of how those cells could be spatially distributed within the NP. In this regard, ABMs can be useful in identifying novel parameters and interactions implicated in IDD and therefore guiding future experiments. Future work could extend the ABM to simulate AF and CEP cells in addition to NP cells, or to simulate advanced stages of IDD by including cell migration and adding in immune cells to further investigate the intercellular interactions in the IVD during an immune response.

Network Modeling

Modeling biological networks provides a holistic and intuitive way to investigate, characterize, and understand the complex interactions between biological components. It is a static diagram represented by nodes (molecules) connected by lines (physical or functional interactions between nodes). The nodes are the stimuli or responses of the network, while the lines indicate either inhibition or activation between nodes directly or indirectly through other signaling pathways. The most frequently studied networks are protein-protein interaction (PPI) networks and the most commonly used software are: Cytoscape (Shannon et al., 2003), VisANT (Hu et al., 2004), TopNet (Yu et al., 2004), MAVisto (Schreiber and Schwöbbermeyer, 2005), FANMOD (Wernicke and Rasche, 2006), Pajek (Batagelj and Andrej, 2004), Mfinder (Kashtan et al., 2004), and GraphCrunch (Milenković et al., 2008).

Biological network modeling usually relies on “bottom-up” approaches, where intracellular interactions are simulated to estimate a final cellular response. Two methods can be used to build a network, either by gathering literature information regarding the pathways and the mechanisms that take part in the IVD degeneration (knowledge-based), or directly from experimental data (data-driven). In IVD degeneration, network modeling tries to capture the interactions between complex sets of proteins and their pathways and reveal the complex dynamics behind the imbalance between anabolic and catabolic processes. Identification of known NP cell high level cell regulatory factors is very important in order to integrate all the single stimuli into an IVD regulatory network model (RNM) for cell regulation, through which hypothesis and testing can be explored.

Regulatory network models can highlight the molecular signatures of the underlying pathological mechanisms that drive a condition. The reductionist view, one gene to one disease, is not applicable to a highly multifactorial condition such as IDD, thus cell signaling pathway analysis is of high importance in order to understand the system as a whole. A mechanistic understanding of the condition could pave the way for a mechanism-based biomarker selection for the effective and personalized treatment of IDD (Baumgartner et al., 2021b).

A common data-driven approach for modelling regulatory networks starts with the acquisition of -omics data from a web-based repository or by generating them. Differentially expressed genes or proteins between healthy and diseased samples are found and a functional enrichment analysis is performed in order to identify the most statistically significant pathways that are present in the condition. The final step includes experimental verification of the targets that were identified in-silico. However, constructing regulatory networks from IDD samples is a challenging task due to high sample variation, stage of IVD, type of disc tissue investigated and the chosen method of analysis which could be MassSpec, Microarrays or Next Generation RNA sequencing.

Xu et al. (2021) created a regulatory network behind IDD by combining transcriptomic and proteomic analysis. They hypothesized that post-transcriptional regulation could have an effect on protein content, thus, if a gene presents elevated mRNA and protein levels, it could be implicated in IDD. Their results identified six genes with these characteristics (CHI3L1, KRT19, COL6A2, DPT, TNFAIP6 and COL11A2), two of which were identified as important IDD markers in independent studies. Another group used transcriptomic data collected from lumbar-degenerated IVDs to build gene regulatory networks, finding differentially expressed genes in chemotactic signaling and matrix-degrading pathways that could later be used to help develop novel pharmacological approaches for IDD treatment (Zamanian et al., 2022). Li H. et al. (2022) constructed a protein-interaction network as well as a disease-gene interaction network that identified two potential therapeutic drugs, entrectinib and larotrectinib, demonstrating how emerging network models can be leveraged to identify novel treatments.

Recently, a top-down network modeling approach was presented to approximate cell responses of NP cells, where the cell is considered as a “black-box” (Baumgartner et al., 2021a; Baumgartner et al., 2021b). Approximations of cell responses were obtained by directly linking key relevant micro-environmental stimuli with cell responses of interest. Therefore, experimentally obtained data was systematically translated into parameters suitable for systems biology approaches. With this novel approach, interrelated results between NP cells of different pro-inflammatory states, i.e., immunopositive for TNF-α, IL-1β, or both; TNF-α and IL-1β, could be obtained for user-defined stimulus environments. This high-level network modeling methodology was embedded within an ABM (see Section “Agent-Based Models”) to visualize a proinflammatory environment and estimate the percentage of cells immunopositive for more than one proinflammatory cytokine, specifically TNF-α and IL-1β. Considering crucial nutritional and biochemical stimuli, in-silico results suggest that pro-inflammatory cytokines are important contributors in catabolic shifts in NP cell responses (Baumgartner et al., 2021b).

Top-down approaches appear promising to tackle highly complex multicellular multifactorial environments, as found in IVD tissues. Amongst others, focus might be set on the integration of more critical stimuli and cell responses in the network model.

Methods for Genetic Analysis in IDD

In the past decades, candidate genes and Genome-wide Association Studies (GWAS) have been implemented in the discovery of genetics underpinnings of complex disorders. The former strategy involves testing the association between a particular gene variant and a trait. Therefore, the selection of the studied gene is led by a priori knowledge of the biological pathways that are involved in the etiology of the disease. However, the high specificity of candidate gene approaches does not reflect the polygenicity in which multiple genomic loci are involved in the development of the disease (Tabor et al., 2002).

On the other hand, the aim of GWAS is to investigate relationships between genetic variants and traits spanning the whole genome in order to give an unbiased and comprehensive view on the allelic architecture underlying complex traits (Goddard et al., 2016). Despite GWAS having identified a large number of risk genomic loci and provided valuable outcomes on the agnostic genetic discovery for complex diseases such as IDD (Uffelmann et al., 2021), there are still gaps that have to be filled. GWAS are not capable of explaining the high heritability of complex diseases. Several reasons could contribute to this limitation, such as epigenetics or epistasis, which is the phenomenon for which the effect of a genetic variant is dependent on the presence of other variants, known as genetic background. Another hypothesis is the “common variants common disease”, which considers that the genetic contribution to a disease would come from an elevated number of SNPs, each one with a very small contribution difficult to identify (Manolio et al., 2009; Gibson 2010). Moreover, the interpretation of how a specific variant affects the downstream biological pathways is very challenging. First, due to linkage disequilibrium (LD), the phenomenon for which variants close to each other are inherited together, associated SNPs are often correlated with other neighboring variants. Thus, checking LD is essential to include all potential causal variants. Further, the majority of associations discovered through GWAS (90%) fall in non-coding regions, hindering the interpretation of how such SNPs affect the phenotype (Cano-Gamez and Trynka, 2020).

For all of these reasons, when interpreting GWAS results, one should consider several factors, including the number of different associations that exist at a given locus and their LD correlation. Then it is possible to pinpoint which could be the causal variant and establish the mechanistic effects on the downstream processes (Cannon and Mohlke, 2018). Many approaches have been carried out to tackle the challenge of understanding how many signals are present at a locus, such as applying a threshold of LD.

Another approach for gathering independent variants, as conceived by Yang et al. (2012), is to perform summary-level statistics conditional analysis. Here, the effect of a lead SNP is tested for association with all the other SNPs at a locus to determine the degree of association between them and detect the independent signals (Schizophrenia Working Group of the Psychiatric Genomics Consortium 2014).

After different signals at a locus have been defined, a subsequent analysis, referred to as “fine mapping”, is performed in order to identify the potential causal variant(s). As reviewed by Schaid et al. (2018), several approaches can be adopted to perform this “fine mapping” analysis: The heuristic LD approach considers all of the SNPs that are related to the main signal with a value of LD higher than a fixed threshold. Other methods that rely more on statistics consist of jointly analyzing the correlation amongst all of the SNPs at a given locus through regression methods. Due to high correlation between variants, penalized regression models have been shown to be the most effective strategy. Recently, Bayesian approaches have been implemented (Jiang et al., 2019) to assess the probability of SNP causality at risk loci with success. When a subset of potential causal SNPs is obtained, further analyses procure the annotation of the variant biological effect, often referred as “Variant to Function” (V2F) analysis (Sun et al., 2021). So far, to relate genes to non-coding SNPs, the closest gene is considered. However, recent works (Mountjoy et al., 2021; Forgetta et al., 2022) aimed to improve this method by integrating various data sources and statistical methods to reach a better interpretation of the effect of non-coding variants on complex traits such as IDD. For this, multiple tools such as Variant effect Predictor (VEP) (McLaren et al., 2016) and Functional Mapping and Annotation of Genome-Wide Association Studies (FUMA) (Watanabe et al., 2017) can be implemented to suggest the potential molecular alterations caused by a SNP.

Additionally, an increasing number of data repositories that have collected information regarding gene expression, regulatory elements, and epigenetics, such as the Encyclopedia of DNA Elements (ENCODE) (de Souza 2012), the roadmap Epigenomics Project (Romanoski et al., 2015), and Genotype-Tissue Expression (GTEx) (The GTEX consortium, 2020) can be queried in order to obtain meaningful insights on the functional elements that could be affected by a candidate variant. In this way, one could retrieve information about the effect of the variant on the expression of a given gene. A gene-variant pair in which the variant is correlated to the gene expression is called expression quantitative trait locus (eQTL). Nicolae et al. (2010) demonstrated that common SNPs associated with complex diseases are significantly more likely to be eQTL, in comparison to rare SNPs acting on the phenotype by altering the expression of the gene rather than modifying the gene itself.

For instance, a study on ischemic stroke by Amini et al. (2020) showed that in blood, the SNP rs78046578 was significantly correlated to the expression of C-X-C motif chemokine ligand 10 (CXCL10), a small protein whose increased levels in the blood serum have been correlated to patients with IDD and LBP as well as with Pfirmann grades (Vazirinejad et al., 2014; Yang et al., 2022).

Moreover, Raine et al. (2014) showed that in cartilage, the polymorphism rs8031440 is correlated with the expression of SMAD3. Specifically, carrying the G allele caused a reduction in the expression of the gene. SMAD3 is a key component of the TGF-β signaling pathway, which is highly involved in the anabolism of the extracellular matrix through enhancing the expression of type I collagen (Verrecchia and Mauviel, 2004). It was reported that SMAD3 knockout mice were smaller, had malformed and kyphotic spines, and had reduced levels of collagen and proteoglycans in the disc (Li et al., 2009). Although further evidence is needed, it appears that carrying a variant that leads to a dysregulation of SMAD3 expression could lead to altered development of the disc and the spine itself.

Machine Learning/AI/Deep Learning

Advances in computer science, computer programming, statistics, mathematics, and modelling allowed the creation of new algorithms that can “learn” and make predictions on new data. Machine learning (ML), AI, and deep learning are relatively new approaches used to tackle the complexity of a disease in order to identify biomarkers, or therapeutic interventions. The complexity and diversity of biological data (i.e., imaging, clinical categories, compound structures, gene sequence, protein/RNA data) is ideal for ML approaches to link seemingly unrelated entities.

A notable usage of deep learning in IDD is found in the classification of disc degeneration based on MRI images using the Pfirrman score. Although widely used, the Pfirrman score is very subjective and different observers sometimes classify the same image with differing scores (Rim 2016). Consistency in grading is essential for a clinician in order to have a clear idea of the patient’s condition, which has led to the development of deep learning models. Niemeyer et al. (2021) have successfully managed to develop a deep learning model for the classification of discs based on MRI data that has an average sensitivity of 90%.

Additionally, ML can be useful in studying disc degeneration through simplifying complex models to decrease computational demand, or through coupling models across several scales to offer a more holistic view of IDD. For example, while patient-specific FEMs, which were explained previously in this review, are useful in studying IVD biomechanics, they usually require complex procedures to set up and long computing times to obtain final simulation results. This therefore prevents prompt feedback to clinicians, resulting in studies with minimal sample sizes and severely hindering its suitability for time-sensitive clinical applications. As a response, neural networks have increasingly been employed in complex dynamical systems, resulting in highly accurate surrogate models that can be evaluated with significantly less computational resources and several orders of magnitude faster than conventional finite element solvers (Pfaff et al., 2020).

Conclusion

IDD is considered a complex multifactorial and pathological disease which alters biomechanical and biochemical aspects of the IVD, resulting in a shifted metabolism associated with increased cytokine and chemokine production by native disc cells driving catabolism, this together with abnormal mechanical loading can result in disc rupture. Following rupture, immune cells are able to invade the disc, allowing crosstalk between the IVD and the immune system. Therefore, a disrupted or herniated disc is needed in order to obtain an entry point for the immune system into the disc. Nevertheless, although catabolic and inflammatory features are different, since IVD cells (NP, AF and CEP) share classical immune cell’s roles and markers (Le Maitre et al., 2005; Jones et al., 2008; Risbud and Shapiro, 2013; Phillips et al., 2013) the resulting phenotype in each case remains controversial.

In terms of terminology, immunomodulatory or inflammatory terms should only be used when immune cells are present within the IVD. Consequently, the production of chemokines and cytokines without disc rupture are contributed by the native IVD cells (NP, AF and CEPs) and should be termed catabolic cytokines and chemokines. Hence, catabolic phenotypes which are related to production of inflammatory cytokines and chemokines are incorrectly commonly attributed to inflammatory processes and easily misconstrued in the literature. Further, both processes could also simultaneously appear during IDD, making the recognition of these phenotypes more challenging.

Regarding the shift from anabolism to catabolism, it is expected that many factors are at play, including changes in the microenvironment, biomechanics, genetics, and metabolism. Within the microenvironment, low glucose and high lactic acid levels contribute to a catabolic shift in IDD, while hypo-osmolarity can activate pro-inflammatory and catabolic factors. Altered biomechanics also contributes to this catabolic shift, and aberrant loading can lead to CEP fractures and AF fissures that allow for immune cell infiltration. IDD has shown high inheritability, and gene variants in genes of structural proteins and their turnover as well as cytokines such as IL-1β have been shown to provoke a catabolic shift. While bacterial presence in the IVD is still controversial, some studies indicate that in some cases of IDD, the increased catabolic cytokines could be due to bacterial infection and treatment with antibiotics could be effective to reduce pain although the results are varied. Overall, IDD is highly multifactorial, and each of these factors discussed play a role in the shift to catabolism within an intact disc and possible immune cell infiltration following AF or CEP rupture.

However, for the best clinical treatment, early diagnosis is crucial, which means that data analysis must be streamlined and the disc pathology must be classified correctly. Therefore, it is necessary to understand how different methodologies have been used to study the different features of IDD and assign the correct terminology. For that purpose, a wide range of in-silico, in-vivo and in-vitro models have been discussed in this review to select the best approaches for future IDD studies and thus provide the best clinical output. While 3D cell culture is effective in investigating individual parameters within IDD, bioreactors and microfluidics studies offer another level of complexity through the addition of mechanics and/or fluid flow. Further, animal models provide even more sophistication as they include interactions between the disc and other tissues, however there are still many biochemical and biological differences in comparison to IDD in humans. This is in part how in-silico studies can be useful, as they can predict changes in the IVD based on prior research, without harming humans. Additionally, computational modeling can offer insights into IDD that are difficult or expensive to obtain through experiments. FEMs are useful to determine the biomechanical effects on the disc, which are expected to simulate specific patient models and observe the effects of CEP shape and its implications on nutrient and water transport, as well as on the different NP morphologies. ABMs can offer visual and spatial predictions, and network models provide insight into complex interactions at the protein or genetic level. Moreover, methodologies using candidate genes and GWAS have identified influential gene variants in IDD. Machine learning can then be a useful tool to simplify these models and methodologies, or to streamline and reduce the bias in the classification of IDD.

Hence, this review summarizes the recent advances of cross-disciplinary approaches to identify the mechanisms of the shift of anabolism to catabolism in the progress of IDD and compared them with immunomodulatory features. It demonstrates our current knowledge of the interplay of the immune system, metabolism, genetics, epigenetics, physiology, and mechanics, as well as computational and experimental models used to investigate catabolism and inflammation in the IVD.

Statements

Author contributions

PB-L and KC conceptionalisation, major literature research and wrote the main text, AN, KC, PB-L, CM, RC and FG and contributed with writing and reviewing and visualisation, ST, EK, AA, LB and EM-M writing and reviewing, LA, LG, KW-K, CM, JN, LG and BG writing and reviewing conceptionalisation, and funding. All authors approved the final version of the manuscript.

Funding

Financial support was received from the Marie Skłodowska Curie International Training Network (ITN) “disc4all” (https://disc4all.upf.edu, accessed on 10 May 2022) grant agreement #955735 (https://cordis.europa.eu/project/id/955735, accessed on 10 May 2022).

Acknowledgments

We thank the Marie Skłodowska Curie International Training Network (ITN) “disc4all” for financial support. Furthermore, we thank the entire disc4all consortium for enabling this review.

Conflict of interest

EK and LK were employed by ProtATonce Ltd.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Glossary

  • ABM

    Agent-based model

  • ADAMTS

    a disintegrant and metalloproteinase with thrombospondin motifs

  • AF

    Annulus fibrosus

  • AI

    Artificial intelligence

  • BNB

    Blood NP barrier

  • BMI

    Body mass index

  • CEP

    Cartilage endplate

  • CLIP

    Cartilage Intermediate layer protein

  • CoNS

    Coagulase-negative Staphylococci

  • CS

    Chondroitin sulfate

  • CXCL10

    C-X-C motif chemokine ligand 10

  • DR

    Death receptor

  • ECM

    Extracellular Matrix

  • ENCODE

    Encyclopedia of DNA Elements

  • ENOS

    Endothelial nitric oxide synthase

  • eQTL

    expression quantitative trait locus

  • ER

    Endoplasmic Reticulum

  • FEM

    Finite element model

  • FUMA

    Functional Mapping and Annotation of Genome-Wide Association Studies

  • GAG

    Glycosaminoglycan

  • GCN

    Graph convolutional networks

  • GDF

    Growth differentiation factor

  • GTEx

    Genotype-Tissue Expression

  • GWAS

    Genome-wide Association Studies

  • HIF

    Hypoxia-inducible factor

  • IDD

    Intervertebral disc degeneration

  • IGF

    Insulin growth factor

  • IGFR

    Insulin growth factor receptor

  • IL

    Interleukin

  • IVD

    Intervertebral disc

  • LBP

    Low Back Pain

  • LD

    Linkage disequilibrium

  • LDD

    Lumbar Disc Degeneration

  • LPS

    Lipopolysaccharide

  • miRNA

    micro-RNA

  • MC

    Modic change

  • ML

    Machine learning

  • MMP

    matrix metalloproteinase

  • MRI

    Magnetic resonance imaging

  • NFκ-B

    nuclear factor kappa-B

  • NP

    Nucleus pulposus

  • SNP

    Single nucleotide polymorphism

  • SP1

    Specificity Protein 1 (syn. Sp1 transcription factor)

  • SPARC

    Secrete protein, acidic, rich in cysteine

  • THBS2

    Thrombospondin-2 gene

  • TonEBP

    Tonicity-responsive enhancer binding protein

  • TGF

    Tumor growth factor

  • TNF-α

    Tumor necrosis factor alpha

  • TRAIL

    TNF-α-related apoptosis-inducing ligand

  • TRP

    Transient receptor potential

  • V2F

    Variant to Function

  • VEP

    Variant effect Predictor

  • VEGF

    Vascular endothelial growth factor

Citation

Bermudez-Lekerika P, Crump KB, Tseranidou S, Nüesch A, Kanelis E, Alminnawi A, Baumgartner L, Muñoz-Moya E, Compte R, Gualdi F, Alexopoulos LG, Geris L, Wuertz-Kozak K, Le Maitre CL, Noailly J and Gantenbein B (2022) Immuno-Modulatory Effects of Intervertebral Disc Cells. Front. Cell Dev. Biol. 10:924692. doi: 10.3389/fcell.2022.924692

References

  • 1

    AdamsM. A.FreemanB. J. C.MorrisonH. P.NelsonI. W.DolanP. (2000). Mechanical Initiation of Intervertebral Disc Degeneration. Spine25, 16251636. 10.1097/00007632-200007010-00005

  • 2

    AdamsM. A.McMillanD. W.GreenT. P.DolanP. (1996). Sustained Loading Generates Stress Concentrations in Lumbar Intervertebral Discs. Spine21, 434438. 10.1097/00007632-199602150-00006

  • 3

    AdamsM. A.RoughleyP. J. (2006). What Is Intervertebral Disc Degeneration, and what Causes it?Spine31, 21512161. 10.1097/01.brs.0000231761.73859.2c

  • 4

    Al DhaheriN.WuN.ZhaoS.WuZ.BlankR. D.ZhangJ.et al (2020). KIAA1217 : A Novel Candidate Gene Associated with Isolated and Syndromic Vertebral Malformations. Am. J. Med. Genet.182, 16641672. 10.1002/ajmg.a.61607

  • 5

    AlbertH. B.LambertP.RollasonJ.SorensenJ. S.WorthingtonT.PedersenM. B.et al (2013a). Does Nuclear Tissue Infected with Bacteria Following Disc Herniations Lead to Modic Changes in the Adjacent Vertebrae?Eur. Spine J.22, 690696. 10.1007/s00586-013-2674-z

  • 6

    AlbertH. B.SorensenJ. S.ChristensenB. S.MannicheC. (2013b). Antibiotic Treatment in Patients with Chronic Low Back Pain and Vertebral Bone Edema (Modic Type 1 Changes): a Double-Blind Randomized Clinical Controlled Trial of Efficacy. Eur. Spine J.22, 697707. 10.1007/s00586-013-2675-y

  • 7

    AliniM.EisensteinS. M.ItoK.LittleC.KettlerA. A.MasudaK.et al (2008). Are Animal Models Useful for Studying Human Disc Disorders/degeneration?Eur. Spine J.17, 219. 10.1007/s00586-007-0414-y

  • 8

    AlkhatibB.RosenzweigD. H.RosenzweigD.KrockE.RoughleyP.BeckmanL.et al (2014). Acute Mechanical Injury of the Human Intervertebral Disc: Link to Degeneration and Pain. eCM28, 98111. ; discussion 110-1. 10.22203/ecm.v028a08

  • 9

    Alvarez-GalvezJ.Suarez-LledoV. (2019). Using Agent-Based Modeling to Understand the Emergence and Reproduction of Social Inequalities in Health. Proceedings44, 2. 10.3390/IECEHS-2-06372

  • 10

    AminiH.ShroffN.StamovaB.FerinoE.Carmona-MoraP.ZhanX.et al (2020). Genetic Variation Contributes to Gene Expression Response in Ischemic Stroke: an eQTL Study. Ann. Clin. Transl. Neurol.7, 16481660. 10.1002/acn3.51154

  • 11

    AnjankarS. D.PoornimaS.RajuS.JaleelM. A.BhiladvalaD.HasanQ. (2015). Degenerated Intervertebral Disc Prolapse and its Association of Collagen I Alpha 1 Spl Gene Polymorphism. Ijoo49, 589594. 10.4103/0019-5413.168765

  • 12

    AramburuJ.López-RodríguezC. (2019). Regulation of Inflammatory Functions of Macrophages and T Lymphocytes by NFAT5. Front. Immunol.10, 535. 10.3389/fimmu.2019.00535

  • 13

    AubinG. G.Baud’huinM.LavigneJ.-P.BrionR.GouinF.LepelletierD.et al (2017). Interaction of Cutibacterium ( Formerly Propionibacterium) Acnes with Bone Cells: a Step toward Understanding Bone and Joint Infection Development. Sci. Rep.7, 42918. 10.1038/srep42918

  • 14

    BarryF. P.NeameP. J.SasseJ.PearsonD. (1994). Length Variation in the Keratan Sulfate Domain of Mammalian Aggrecan. Matrix Biol.14, 323328. 10.1016/0945-053x(94)90198-8

  • 15

    BatageljV.MrvarA. (2004). “Pajek - Analysis and Visualization of Large Networks,” in Graph Drawing Software. Editor JüngerM. M. P. (Berlin, Heidelberg: Springer Berlin Heidelberg), 77103. 10.1007/978-3-642-18638-7_4

  • 16

    BattiéM. C.VidemanT. (2006). Lumbar Disc Degeneration: Epidemiology and Genetics. JBJS88, 39. 10.2106/JBJS.E.01313

  • 17

    BattiéM. C.VidemanT.ParentE. (2004). Lumbar Disc Degeneration. Spine29, 26792690. 10.1097/01.brs.0000146457.83240.eb

  • 18

    BaumgartnerL.ReaghJ. J.González BallesterM. A.NoaillyJ. (2021a). Simulating Intervertebral Disc Cell Behaviour within 3D Multifactorial Environments. Bioinformatics37, 12461253. 10.1093/bioinformatics/btaa939

  • 19

    BaumgartnerL.SadowskaA.TíoL.González BallesterM. A.Wuertz-KozakK.NoaillyJ. (2021b). Evidence-Based Network Modelling to Simulate Nucleus Pulposus Multicellular Activity in Different Nutritional and Pro-inflammatory Environments. Front. Bioeng. Biotechnol.9, 734258. 10.3389/fbioe.2021.734258

  • 20

    BaumgartnerL.Wuertz-KozakK.Le MaitreC. L.WignallF.RichardsonS. M.HoylandJ.et al (2021c). Multiscale Regulation of the Intervertebral Disc: Achievements in Experimental, In Silico, and Regenerative Research. Ijms22, 703. 10.3390/ijms22020703

  • 21

    BecksteinJ. C.SenS.SchaerT. P.VresilovicE. J.ElliottD. M. (2008). Comparison of Animal Discs Used in Disc Research to Human Lumbar Disc. Spine33, E166E173. 10.1097/BRS.0b013e318166e001

  • 22

    BedoreJ.ShaW.McCannM. R.LiuS.LeaskA.SéguinC. A. (2013). Loss of Notochord-Derived CCN2 Results in Impaired Intervertebral Disc Development and Premature Disc Degeneration. Arthritis & Rheumatism65, an. 10.1002/art.38075

  • 23

    BennekerL. M.HeiniP. F.AliniM.AndersonS. E.ItoK. (2005). 2004 Young Investigator Award Winner: Vertebral Endplate Marrow Contact Channel Occlusions and Intervertebral Disc Degeneration. Spine30, 167173. 10.1097/01.brs.0000150833.93248.09

  • 24

    BibbyS. R. S.UrbanJ. P. G. (2004). Effect of Nutrient Deprivation on the Viability of Intervertebral Disc Cells. Eur. Spine J.13, 695701. 10.1007/s00586-003-0616-x

  • 25

    BinchA. L. A.ColeA. A.BreakwellL. M.MichaelA. L. R.ChivertonN.CreemersL. B.et al (2015a). Class 3 Semaphorins Expression and Association with Innervation and Angiogenesis within the Degenerate Human Intervertebral Disc. Oncotarget6, 1833818354. 10.18632/oncotarget.4274

  • 26

    BinchA. L. A.ColeA. A.BreakwellL. M.MichaelA. L. R.ChivertonN.CreemersL. B.et al (2015b). Nerves Are More Abundant Than Blood Vessels in the Degenerate Human Intervertebral Disc. Arthritis Res. Ther.17, 370. 10.1186/s13075-015-0889-6

  • 27

    BlumenthalS.McAfeeP. C.GuyerR. D.HochschulerS. H.GeislerF. H.HoltR. T.et al (2005). A Prospective, Randomized, Multicenter Food and Drug Administration Investigational Device Exemptions Study of Lumbar Total Disc Replacement with the CHARITÉ Artificial Disc versus Lumbar Fusion. Spine30, 15651575. discussion E387-91. 10.1097/01.brs.0000170587.32676.0e

  • 28

    BinchA. L.ColeA. A.BreakwellL. M.MichaelA. L.ChivertonN.CrossA. K.et al (2014). Expression and Regulation of Neurotrophic and Angiogenic Factors during Human Intervertebral Disc Degeneration. Arthritis Res. Ther.16, 416. 10.1186/s13075-014-0416-1

  • 29

    BongersP. M.HulshofC. T. J.DlJKSTRAL.BoshuizenH. C.GroenhoutH. J. M.ValkenE. (1990). Back Pain and Exposure to Whole Body Vibration in Helicopter Pilots. Ergonomics33, 10071026. 10.1080/00140139008925309

  • 30

    BournonvilleS.LambrechtsT.VanhulstJ.LuytenF. P.PapantoniouI.GerisL. (2019). Towards Self-Regulated Bioprocessing: A Compact Benchtop Bioreactor System for Monitored and Controlled 3D Cell and Tissue Culture. Biotechnol. J.14, 1800545. 10.1002/biot.201800545

  • 31

    BråtenL. C. H.GrøvleL.EspelandA.PrippA. H.GrotleM.HelllumC.et al (2020). Clinical Effect Modifiers of Antibiotic Treatment in Patients with Chronic Low Back Pain and Modic Changes - Secondary Analyses of a Randomised, Placebo-Controlled Trial (The AIM Study). Bmc. Musculoskelet. Disord.21, 111. 10.1186/s12891-020-03422-y

  • 32

    BuckwalterJ. A. (1995). Aging and Degeneration of the Human Intervertebral Disc. Spine20, 13071314. 10.1097/00007632-199506000-00022

  • 33

    CambriaE.BrunnerS.HeusserS.FischP.HitzlW.FergusonS. J.et al (2020). Cell-Laden Agarose-Collagen Composite Hydrogels for Mechanotransduction Studies. Front. Bioeng. Biotechnol.8, 346. 10.3389/fbioe.2020.00346

  • 34

    CannonM. E.MohlkeK. L. (2018). Deciphering the Emerging Complexities of Molecular Mechanisms at GWAS Loci. Am. J. Hum. Genet.103, 637653. 10.1016/j.ajhg.2018.10.001

  • 35

    Cano-GamezE.TrynkaG. (2020). From GWAS to Function: Using Functional Genomics to Identify the Mechanisms Underlying Complex Diseases. Front. Genet.11, 424. 10.3389/fgene.2020.00424

  • 36

    CapoorM. N.KoniecznaA.McDowellA.RuzickaF.SmrckaM.JancalekR.et al (2021). Pro-Inflammatory and Neurotrophic Factor Responses of Cells Derived from Degenerative Human Intervertebral Discs to the Opportunistic Pathogen Cutibacterium Acnes. Ijms22, 2347. 10.3390/ijms22052347

  • 37

    ChanS. C. W.Gantenbein-RitterB. (2012). Preparation of Intact Bovine Tail Intervertebral Discs for Organ Culture. JoVE60, e3490. 10.3791/3490

  • 38

    ChanS. C. W.WalserJ.KäppeliP.ShamsollahiM. J.FergusonS. J.Gantenbein-RitterB. (2013). Region Specific Response of Intervertebral Disc Cells to Complex Dynamic Loading: An Organ Culture Study Using a Dynamic Torsion-Compression Bioreactor. PLoS. ONE.8, e72489. 10.1371/journal.pone.0072489

  • 39

    ChangH.CaiF.ZhangY.XueM.LiuL.YangA.et al (2017). Early-stage Autophagy Protects Nucleus Pulposus Cells from Glucose Deprivation-Induced Degeneration via the P -eIF2α/ATF4 Pathway. Biomed. Pharmacother.89, 529535. 10.1016/j.biopha.2017.02.074

  • 40

    ChangH.WangH.YangX.YouK.JiangM.CaiF.et al (2021). Comprehensive Profile Analysis of Differentially Expressed circRNAs in Glucose Deprivation-Induced Human Nucleus Pulposus Cell Degeneration. BioMed Res. Int.2021, 114. 10.1155/2021/4770792

  • 41

    ChenJ.-W.LiB.YangY.-H.JiangS.-D.JiangL.-S. (2014). Significance of Hypoxia in the Physiological Function of Intervertebral Disc Cells. Crit. Rev. Eukaryot. Gene. Expr.24, 193204. 10.1615/critreveukaryotgeneexpr.2014010485

  • 42

    ChenJ.BaerA. E.PaikP. Y.YanW.SettonL. A. (2002). Matrix Protein Gene Expression in Intervertebral Disc Cells Subjected to Altered Osmolarity. Biochem. Biophysical Res. Commun.293, 932938. 10.1016/S0006-291X(02)00314-5

  • 43

    CheungK. M. C.KarppinenJ.ChanD.HoD. W. H.SongY.-Q.ShamP.et al (2009). Prevalence and Pattern of Lumbar Magnetic Resonance Imaging Changes in a Population Study of One Thousand Forty-Three Individuals. Spine34, 934940. 10.1097/BRS.0b013e3181a01b3f

  • 44

    ChingC. T. S.ChowD. H. K.YaoF. Y. D.HolmesA. D. (2003). The Effect of Cyclic Compression on the Mechanical Properties of the Inter-vertebral Disc: an In Vivo Study in a Rat Tail Model. Clin. Biomech.18, 182189. 10.1016/s0268-0033(02)00188-2

  • 45

    ChoiH.ChaiyamongkolW.DoolittleA. C.JohnsonZ. I.GogateS. S.SchoepflinZ. R.et al (2018). COX-2 Expression Mediated by Calcium-TonEBP Signaling axis under Hyperosmotic Conditions Serves Osmoprotective Function in Nucleus Pulposus Cells. J. Biol. Chem.293, 89698981. 10.1074/jbc.RA117.001167

  • 46

    ChouP.-H.WangS.-T.YenM.-H.LiuC.-L.ChangM.-C.LeeO. K.-S. (2016). Fluid-induced, Shear Stress-Regulated Extracellular Matrix and Matrix Metalloproteinase Genes Expression on Human Annulus Fibrosus Cells. Stem. Cell. Res. Ther.7, 34. 10.1186/s13287-016-0292-5

  • 47

    ChuG.ShiC.WangH.ZhangW.YangH.LiB. (2018). Strategies for Annulus Fibrosus Regeneration: From Biological Therapies to Tissue Engineering. Front. Bioeng. Biotechnol.6, 90. 10.3389/fbioe.2018.00090

  • 48

    CostiJ. J.StokesI. A.Gardner-MorseM. G.IatridisJ. C. (2008). Frequency-Dependent Behavior of the Intervertebral Disc in Response to Each of Six Degree of Freedom Dynamic Loading. Spine33, 17311738. 10.1097/BRS.0b013e31817bb116

  • 49

    CroftA. S.RothY.OswaldK. A. C.ĆorlukaS.Bermudez-LekerikaP.GantenbeinB. (2021). In Situ Cell Signalling of the Hippo-YAP/TAZ Pathway in Reaction to Complex Dynamic Loading in an Intervertebral Disc Organ Culture. Ijms22, 13641. 10.3390/ijms222413641

  • 50

    DahiaC. L.EngilesJ. B.GullbrandS. E.IatridisJ. C.LaiA.Le MaitreC. L.et al (2021). A Perspective on the ORS Spine Section Initiative to Develop a Multi-Species JOR Spine Histopathology Series. JOR. Spine4, e1165. 10.1002/jsp2.1165

  • 51

    DaiJ.XingY.XiaoL.LiJ.CaoR.HeY.et al (2019). Microfluidic Disc-On-A-Chip Device for Mouse Intervertebral Disc-Pitching a Next-Generation Research Platform to Study Disc Degeneration. ACS Biomater. Sci. Eng.5, 20412051. 10.1021/acsbiomaterials.8b01522

  • 52

    DalyC.GhoshP.JenkinG.OehmeD.GoldschlagerT. (2016). A Review of Animal Models of Intervertebral Disc Degeneration: Pathophysiology, Regeneration, and Translation to the Clinic. BioMed Res. Int.2016, 114. 10.1155/2016/5952165

  • 53

    DamiatiS.KompellaU.DamiatiS.KodziusR. (2018). Microfluidic Devices for Drug Delivery Systems and Drug Screening. Genes9, 103. 10.3390/genes9020103

  • 54

    DanielsJ.BinchA. A. L.Le MaitreC. L. (2017). Inhibiting IL-1 Signaling Pathways to Inhibit Catabolic Processes in Disc Degeneration. J. Orthop. Res.35, 7485. 10.1002/jor.23363

  • 55

    de SouzaN. (2012). The ENCODE Project. Nat. Methods.9, 1046. 10.1038/nmeth.2238

  • 56

    DongD. M.YaoM.LiuB.SunC. Y.JiangY. Q.WangY. S. (2007). Association between the -1306C/T Polymorphism of Matrix Metalloproteinase-2 Gene and Lumbar Disc Disease in Chinese Young Adults. Eur. Spine J.16, 19581961. 10.1007/s00586-007-0454-3

  • 57

    DuH.BaiB.QiuY.YinS.BianW. (2015). Association between TRAIL Gene Polymorphisms and the Susceptibility and Severity of Lumbar Disc Degeneration. Int. J. Clin. Exp. Pathol.8, 74157420.

  • 58

    DuY.TavanaS.RahmanT.BaxanN.HansenU. N.NewellN. (2021). Sensitivity of Intervertebral Disc Finite Element Models to Internal Geometric and Non-geometric Parameters. Front. Bioeng. Biotechnol.9, 660013. 10.3389/fbioe.2021.660013

  • 59

    DudliS.FieldsA. J.SamartzisD.KarppinenJ.LotzJ. C. (2016). Pathobiology of Modic Changes. Eur. Spine J.25, 37233734. 10.1007/s00586-016-4459-7

  • 60

    DudliS.HaschtmannD.FergusonS. J. (2011). Fracture of the Vertebral Endplates, but Not Equienergetic Impact Load, Promotes Disc Degeneration In Vitro. J. Orthop. Res.30, 809816. 10.1002/jor.21573

  • 61

    DudliS.HaschtmannD.FergusonS. J. (2014). Persistent Degenerative Changes in the Intervertebral Disc after Burst Fracture in an In Vitro Model Mimicking Physiological Post-traumatic Conditions. Eur. Spine J.24, 19011908. 10.1007/s00586-014-3301-3

  • 62

    DudliS.LiebenbergE.MagnitskyS.LuB.LauricellaM.LotzJ. C. (2018). Modic Type 1 Change Is an Autoimmune Response that Requires a Proinflammatory Milieu provided by the 'Modic Disc'. Spine J.18, 831844. 10.1016/j.spinee.2017.12.004

  • 63

    DuncanL. E.OstacherM.BallonJ. (2019). How Genome-wide Association Studies (GWAS) Made Traditional Candidate Gene Studies Obsolete. Neuropsychopharmacol.44, 15181523. 10.1038/s41386-019-0389-5

  • 64

    ElfervigM. K.MinchewJ. T.FranckeE.TsuzakiM.BanesA. J. (2001a). IL-1? Sensitizes Intervertebral Disc Annulus Cells to Fluid-Induced Shear Stress. J. Cell. Biochem.82, 290298. 10.1002/jcb.1153

  • 65

    ElfervigM. K.MinchewJ. T.FranckeE.TsuzakiM.BanesA. J. (2001b). IL-1? Sensitizes Intervertebral Disc Annulus Cells to Fluid-Induced Shear Stress. J. Cell. Biochem.82, 290298. 10.1002/jcb.1153

  • 66

    FadokV. A. (1999). Clearance: The Last and Often Forgotten Stage of Apoptosis. J. Mamm. Gl. Biol. Neoplasia4, 203211. 10.1023/A:1011384009787

  • 67

    FearingB. V.HernandezP. A.SettonL. A.ChahineN. O. (2018). Mechanotransduction and Cell Biomechanics of the Intervertebral Disc. JOR. Spine1, e1026. 10.1002/jsp2.1026

  • 68

    FengY.EganB.WangJ. (2016). Genetic Factors in Intervertebral Disc Degeneration. Genes & Dis.3, 178185. 10.1016/j.gendis.2016.04.005

  • 69

    FernandoH. N.CzamanskiJ.YuanT.-Y.GuW.SalahadinA.HuangC.-Y. C. (2011). Mechanical Loading Affects the Energy Metabolism of Intervertebral Disc Cells. J. Orthop. Res.29, 16341641. 10.1002/jor.21430

  • 70

    FlanneryC. R.LittleC. B.CatersonB. (1998). Molecular Cloning and Sequence Analysis of the Aggrecan Interglobular Domain from Porcine, Equine, Bovine and Ovine Cartilage: Comparison of Proteinase-Susceptible Regions and Sites of Keratan Sulfate Substitution. Matrix Biol.16, 507511. 10.1016/s0945-053x(98)90021-x

  • 71

    ForgettaV.JiangL.VulpescuN. A.HoganM. S.ChenS.MorrisJ. A.et al (2022). An Effector Index to Predict Target Genes at GWAS Loci. Hum. Genet.[Epub ahead of print]. 10.1007/s00439-022-02434-z

  • 72

    FrauchigerD.TekariA.TekariA.WöltjeM.FortunatoG.BennekerL.et al (2017). A Review of the Application of Reinforced Hydrogels and Silk as Biomaterials for Intervertebral Disc Repair. eCM34, 271290. 10.22203/eCM.v034a17

  • 73

    FreemontA. J.WatkinsA.Le MaitreC.BairdP.JeziorskaM.KnightM. T. N.et al (2002). Nerve Growth Factor Expression and Innervation of the Painful Intervertebral Disc. J. Pathol.197, 286292. 10.1002/path.1108

  • 74

    FreemontA.PeacockT.GoupilleP.HoylandJ.O'BrienJ.JaysonM. (1997). Nerve Ingrowth into Diseased Intervertebral Disc in Chronic Back Pain. Lancet350, 178181. 10.1016/s0140-6736(97)02135-1

  • 75

    FrobinW.BrinckmannP.KramerM.HartwigE. (2001). Height of Lumbar Discs Measured from Radiographs Compared with Degeneration and Height Classified from MR Images. Eur. Radiol.11, 263269. 10.1007/s003300000556

  • 76

    GalbuseraF.SchmidtH.Neidlinger-WilkeC.WilkeH.-J. (2011a). The Effect of Degenerative Morphological Changes of the Intervertebral Disc on the Lumbar Spine Biomechanics: a Poroelastic Finite Element Investigation. Comput. Methods Biomechanics Biomed. Eng.14, 729739. 10.1080/10255842.2010.493522

  • 77

    GalbuseraF.SchmidtH.NoaillyJ.MalandrinoA.LacroixD.WilkeH.-J.et al (2011b). Comparison of Four Methods to Simulate Swelling in Poroelastic Finite Element Models of Intervertebral Discs. J. Mech. Behav. Biomed. Mater.4, 12341241. 10.1016/j.jmbbm.2011.04.008

  • 78

    GantenbeinB.FrauchigerD. A.MayR. D.BakirciE.UrsR.SibylleG. (2019). “Developing Bioreactors to Host Joint-Derived Tissues that Require Mechanical Stimulation,” in Encyclopedia of Tissue Engineering and Regenerative Medicine. Editors ReisR. L.GomesM. E. (NY: Elsevier), 261280. 10.1016/b978-0-12-801238-3.65611-8

  • 79

    GantenbeinB.CalandrielloE.Wuertz-KozakK.BennekerL. M.KeelM. J.ChanS. C. (2014). Activation of Intervertebral Disc Cells by Co-culture with Notochordal Cells, Conditioned Medium and Hypoxia. BMC Musculoskelet. Disord.15, 422. 10.1186/1471-2474-15-422

  • 80

    GantenbeinB.GrünhagenT.LeeC. R.van DonkelaarC. C.AliniM.ItoK. (2006). An In Vitro Organ Culturing System for Intervertebral Disc Explants with Vertebral Endplates. Spine31, 26652673. 10.1097/01.brs.0000244620.15386.df

  • 81

    GantenbeinB.Illien-JüngerS.ChanS.WalserJ.HaglundL.FergusonS.et al (2015). Organ Culture Bioreactors - Platforms to Study Human Intervertebral Disc Degeneration and Regenerative Therapy. Cscr10, 339352. 10.2174/1574888X10666150312102948

  • 82

    GBD 2017 Disease and Injury Incidence and Prevalence Collaborators (2018). Global, Regional, and National Incidence, Prevalence, and Years Lived with Disability for 354 Diseases and Injuries for 195 Countries and Territories, 1990-2017: A Systematic Analysis for the Global Burden of Disease Study 2017. Lancet392, 17891858. 10.1016/S0140-6736(18)32279-7

  • 83

    GhezelbashF.SchmidtH.Shirazi-AdlA.El-RichM. (2020). Internal Load-Sharing in the Human Passive Lumbar Spine: Review of In Vitro and Finite Element Model Studies. J. Biomechanics102, 109441. 10.1016/j.jbiomech.2019.109441

  • 84

    GibsonG. (2010). Hints of Hidden Heritability in GWAS. Nat. Genet.42, 558560. 10.1038/ng0710-558

  • 85

    GilbertH. T.HoylandJ. A.FreemontA. J.Millward-SadlerS. J. (2011). The Involvement of Interleukin-1 and Interleukin-4 in the Response of Human Annulus Fibrosus Cells to Cyclic Tensile Strain: an Altered Mechanotransduction Pathway with Degeneration. Arthritis Res. Ther.13, R8. 10.1186/ar3229

  • 86

    GilbertH. T. J.HodsonN.BairdP.RichardsonS. M.HoylandJ. A. (2016). Acidic pH Promotes Intervertebral Disc Degeneration: Acid-Sensing Ion Channel -3 as a Potential Therapeutic Target. Sci. Rep.6, 37360. 10.1038/srep37360

  • 87

    GilliganC. J.CohenS. P.FischettiV. A.HirschJ. A.CzaplewskiL. G. (2021). Chronic Low Back Pain, Bacterial Infection and Treatment with Antibiotics. Spine J.21, 903914. 10.1016/j.spinee.2021.02.013

  • 88

    GoddardM. E.KemperK. E.MacLeodI. M.ChamberlainA. J.HayesB. J. (2016). Genetics of Complex Traits: Prediction of Phenotype, Identification of Causal Polymorphisms and Genetic Architecture. Proc. R. Soc. B283, 20160569. 10.1098/rspb.2016.0569

  • 89

    GorthD. J.ShapiroI. M.RisbudM. V. (2015). Discovery of the Drivers of Inflammation Induced Chronic Low Back Pain: from Bacteria to Diabetes. Discov. Med.20, 177184.

  • 90

    GorthD. J.ShapiroI. M.RisbudM. V. (2019). A New Understanding of the Role of IL-1 in Age-Related Intervertebral Disc Degeneration in a Murine Model. J. Bone. Min. Res.34, 15311542. 10.1002/jbmr.3714

  • 91

    GorthD. J.ShapiroI. M.RisbudM. V. (2018). Transgenic Mice Overexpressing Human TNF-α Experience Early Onset Spontaneous Intervertebral Disc Herniation in the Absence of Overt Degeneration. Cell. death. Dis.10, 114. 10.1038/s41419-018-1246-x

  • 92

    Granville SmithI.DanckertN. P.FreidinM. B.WellsP.MarchesiJ. R.WilliamsF. M. K. (2021). Evidence for Infection in Intervertebral Disc Degeneration: a Systematic Review. Eur. Spine J.31, 414430. 10.1007/s00586-021-07062-1

  • 93

    GuJ.GuanF.GuanG.XuG.WangX.ZhaoW.et al (2013). Aggrecan Variable Number of Tandem Repeat Polymorphism and Lumbar Disc Degeneration. Spine38, E1600E1607. 10.1097/BRS.0000000000000012

  • 94

    GuerreroJ.HäckelS.CroftA. S.AlbersC. E.GantenbeinB. (2020). The Effects of 3D Culture on the Expansion and Maintenance of Nucleus Pulposus Progenitor Cell Multipotency. JOR Spine4, e1131. 10.1002/jsp2.1131

  • 95

    GullbrandS. E.MalhotraN. R.SchaerT. P.ZawackiZ.MartinJ. T.BendigoJ. R.et al (2017). A Large Animal Model that Recapitulates the Spectrum of Human Intervertebral Disc Degeneration. Osteoarthr. Cartil.25, 146156. 10.1016/j.joca.2016.08.006

  • 96

    GullbrandS. E.SmithL. J.SmithH. E.MauckR. L. (2018). Promise, Progress, and Problems in Whole Disc Tissue Engineering. JOR. Spine1, e1015. 10.1002/jsp2.1015

  • 97

    GuoX.ChenG. (2020). Hypoxia-Inducible Factor Is Critical for Pathogenesis and Regulation of Immune Cell Functions in Rheumatoid Arthritis. Front. Immunol.11, 1668. 10.3389/fimmu.2020.01668

  • 98

    HaglundL.MoirJ.BeckmanL.MulliganK. R.JimB.OuelletJ. A.et al (2011). Development of a Bioreactor for Axially Loaded Intervertebral Disc Organ Culture. Tissue Eng. Part C. Methods17, 10111019. 10.1089/ten.TEC.2011.0025

  • 99

    HamiltonD. J.SéguinC. A.WangJ.PilliarR. M.KandelR. A. (2006). BioEngineering of Skeletal Tissues TeamFormation of a Nucleus Pulposus-Cartilage Endplate Construct In Vitro. Biomaterials27, 397405. 10.1016/j.biomaterials.2005.07.007

  • 100

    HanI. B.RopperA. E.TengY. D.ShinD. A.JeonY. J.ParkH. M.et al (2013). Association between VEGF and eNOS Gene Polymorphisms and Lumbar Disc Degeneration in a Young Korean Population. Genet. Mol. Res.12, 22942305. 10.4238/2013.July.8.10

  • 101

    HellumC.JohnsenL. G.JohnsenBergL. G.GjertsenØ.BergL.NeckelmannG.et al (2012). Predictors of Outcome after Surgery with Disc Prosthesis and Rehabilitation in Patients with Chronic Low Back Pain and Degenerative Disc: 2-year Follow-Up. Eur. Spine J.21, 681690. 10.1007/s00586-011-2145-3

  • 102

    HernandezP. A.JacobsenT. D.BaratiZ.ChahineN. O. (2020a). Confocal Scanning of Intervertebral Disc Cells in 3D : Inside Alginate Beads and in Native Microenvironment. JOR Spine3, e1106. 10.1002/jsp2.1106

  • 103

    HernandezP. A.JacobsenT. D.ChahineN. O. (2020b). Actomyosin Contractility Confers Mechanoprotection against TNFα-Induced Disruption of the Intervertebral Disc. Sci. Adv.6, eaba2368. 10.1126/sciadv.aba2368

  • 104

    HingertD.NawilaijaroenP.AldridgeJ.BarantoA.BrisbyH. (2019). Investigation of the Effect of Secreted Factors from Mesenchymal Stem Cells on Disc Cells from Degenerated Discs. Cells. Tissues. Organs208, 7688. 10.1159/000506350

  • 105

    HiroseY.ChibaK.KarasugiT.NakajimaM.KawaguchiY.MikamiY.et al (2008). A Functional Polymorphism in THBS2 that Affects Alternative Splicing and MMP Binding Is Associated with Lumbar-Disc Herniation. Am. J. Hum. Genet.82, 11221129. 10.1016/j.ajhg.2008.03.013

  • 106

    HolmS.HolmA. K.EkströmL.KarladaniA.HanssonT. (2004). Experimental Disc Degeneration Due to Endplate Injury. J. Spinal Disord. Tech.17, 6471. 10.1097/00024720-200402000-00012

  • 107

    HoogendoornW. E.van PoppelM. N. M.BongersP. M.KoesB. W.BouterL. M. (2000). Systematic Review of Psychosocial Factors at Work and Private Life as Risk Factors for Back Pain. Spine25, 21142125. 10.1097/00007632-200008150-00017

  • 108

    HornerH. A.UrbanJ. P. G. (20012001). 2001 Volvo Award Winner in Basic Science Studies: Effect of Nutrient Supply on the Viability of Cells from the Nucleus Pulposus of the Intervertebral Disc. Spine26, 25432549. 10.1097/00007632-200112010-00006

  • 109

    HoylandJ. A.Le MaitreC.FreemontA. J. (2008). Investigation of the Role of IL-1 and TNF in Matrix Degradation in the Intervertebral Disc. Rheumatology47, 809814. 10.1093/rheumatology/ken056

  • 110

    HuJ.YanQ.ShiC.TianY.CaoP.YuanW. (2017). BMSC Paracrine Activity Attenuates Interleukin-1β-Induced Inflammation and Apoptosis in Rat AF Cells via Inhibiting Relative NF-Κb Signaling and the Mitochondrial Pathway. Am. J. Transl. Res.9, 7989.

  • 111

    HuZ.MellorJ.WuJ.DeLisiC. (2004). VisANT: an Online Visualization and Analysis Tool for Biological Interaction Data. Bmc. Bioinforma.5, 17. 10.1186/1471-2105-5-17

  • 112

    HuangY.-C.HuY.LiZ.LukK. D. K. (2018). Biomaterials for Intervertebral Disc Regeneration: Current Status and Looming Challenges. J. Tissue. Eng. Regen. Med.12, 21882202. 10.1002/term.2750

  • 113

    HunterC. J.De PlaenI. G. (2014). Inflammatory Signaling in NEC: Role of NF-Κb, Cytokines and Other Inflammatory Mediators. Pathophysiology21, 5565. 10.1016/j.pathophys.2013.11.010

  • 114

    HunterC. J. (1988). The Biology of the Intervertebral Disc. Boca Raton, Florida, US: CRC Press.

  • 115

    HwangM. H.SonH.-G.KimJ.ChoiH. (2020). In Vitro model of Distinct Catabolic and Inflammatory Response Patterns of Endothelial Cells to Intervertebral Disc Cell Degeneration. Sci. Rep.10, 20596. 10.1038/s41598-020-77785-6

  • 116

    IatridisJ. C.MenteP. L.StokesI. A. F.AronssonD. D.AliniM. (1999). Compression-induced Changes in Intervertebral Disc Properties in a Rat Tail Model. Spine24, 9961002. 10.1097/00007632-199905150-00013

  • 117

    Illien-JüngerS.WalterB. A.MayerJ. E.HechtA. C.IatridisJ. C. (2014). “Intervertebral Disc Culture Models and Their Applications to Study Pathogenesis and Repair,” in The Intervertebral Disc (Vienna: Springer), 353371. 10.1007/978-3-7091-1535-0_22

  • 118

    IsaI. L. M.MokhtarS. A.AbbahS. A.FauziM. B.DevittA.PanditA. (2022). Intervertebral Disc Degeneration: Biomaterials and Tissue Engineering Strategies toward Precision Medicine. Adv. Healthc. Mater., e2102530. 10.1002/adhm.202102530

  • 119

    JaureguiJ. J.TranA.KaveeshwarS.NadarajahV.ChaudhriM. W.HennR. F.et al (2021). Diagnosing a Periprosthetic Shoulder Infection: A Systematic Review. J. Orthop.26, 5866. 10.1016/j.jor.2021.07.012

  • 120

    JensenC.TengY. (2020). Is it Time to Start Transitioning from 2D to 3D Cell Culture?Front. Mol. Biosci.7, 33. 10.3389/fmolb.2020.00033

  • 121

    JiangJ.ColeJ. B.FreebernE.DaY.VanRadenP. M.MaL. (2019). Functional Annotation and Bayesian Fine-Mapping Reveals Candidate Genes for Important Agronomic Traits in Holstein Bulls. Commun. Biol.2, 112. 10.1038/s42003-019-0454-y

  • 122

    JinL.BalianG.LiX. J. (2018). Animal Models for Disc Degeneration-An Update. Histol. Histopathol.33, 543554. 10.14670/HH-11-910

  • 123

    JohnsonW. E. B.CatersonB.EisensteinS. M.HyndsD. L.SnowD. M.RobertsS. (2002). Human Intervertebral Disc Aggrecan Inhibits Nerve Growth In Vitro. Arthritis & Rheumatism46, 26582664. 10.1002/art.10585

  • 124

    JohnsonW. E. B.CatersonB.EisensteinS. M.RobertsS. (2005). Human Intervertebral Disc Aggrecan Inhibits Endothelial Cell Adhesion and Cell Migration In Vitro. Spine30, 11391147. 10.1097/01.brs.0000162624.95262.73

  • 125

    JohnsonW. E. B.StephanS.RobertsS. (2008). The Influence of Serum, Glucose and Oxygen on Intervertebral Disc Cell Growth In Vitro: Implications for Degenerative Disc Disease. Arthritis Res. Ther.10, R46. 10.1186/ar2405

  • 126

    JohnsonZ. I.DoolittleA. C.SnuggsJ. W.ShapiroI. M.Le MaitreC. L.RisbudM. V. (2017). TNF-α Promotes Nuclear Enrichment of the Transcription Factor TonEBP/NFAT5 to Selectively Control Inflammatory but Not Osmoregulatory Responses in Nucleus Pulposus Cells. J. Biol. Chem.292, 1756117575. 10.1074/jbc.M117.790378

  • 127

    JohnsonZ.SchoepflinZ. R.SchoepflinZ.ChoiH.ShapiroI.RisbudM. (2015). Disc in Flames: Roles of TNF-α and IL-1β in Intervertebral Disc Degeneration. eCM30, 104117. 10.22203/ecm.v030a08

  • 128

    JonesA. C.WilcoxR. K. (2008). Finite Element Analysis of the Spine: towards a Framework of Verification, Validation and Sensitivity Analysis. Med. Eng. Phys.30, 12871304. 10.1016/j.medengphy.2008.09.006

  • 129

    JonesP.GardnerL.MenageJ.WilliamsG. T.RobertsS. (2008). Intervertebral Disc Cells as Competent Phagocytes In Vitro: Implications for Cell Death in Disc Degeneration. Arthritis Res. Ther.10, R86. 10.1186/ar2466

  • 130

    KalbS.MartirosyanN. L.KalaniM. Y. S.BrocG. G.TheodoreN. (2012). Genetics of the Degenerated Intervertebral Disc. World Neurosurg.77, 491501. 10.1016/j.wneu.2011.07.014

  • 131

    KangJ. D.Stefanovic-RacicM.McIntyreL. A.GeorgescuH. I.EvansC. H. (1997). Toward a Biochemical Understanding of Human Intervertebral Disc Degeneration and Herniation. Spine22, 10651073. 10.1097/00007632-199705150-00003

  • 132

    KarasugiT.SembaK.HiroseY.KelempisiotiA.NakajimaM.MiyakeA.et al (2009). Association of the Tag SNPs in the HumanSKTGene (KIAA1217) with Lumbar Disc Herniation. J. Bone Mineral Res.24, 15371543. 10.1359/jbmr.090314

  • 133

    KarppinenJ.DaavittilaI.SolovievaS.KuismaM.TaimelaS.NatriA.et al (2008). Genetic Factors Are Associated with Modic Changes in Endplates of Lumbar Vertebral Bodies. Spine33, 12361241. 10.1097/BRS.0b013e318170fd0e

  • 134

    KashtanN.ItzkovitzS.MiloR.AlonU. (2004). Efficient Sampling Algorithm for Estimating Subgraph Concentrations and Detecting Network Motifs. Bioinformatics20, 17461758. 10.1093/bioinformatics/bth163

  • 135

    KawaguchiY. (2018). Genetic Background of Degenerative Disc Disease in the Lumbar Spine. Spine Surg. Relat. Res.2, 98112. 10.22603/ssrr.2017-0007

  • 136

    KawaguchiY.OsadaR.KanamoriM.IshiharaH.OhmoriK.MatsuiH.et al (1999). Association between an Aggrecan Gene Polymorphism and Lumbar Disc Degeneration. Spine24, 24562460. 10.1097/00007632-199912010-00006

  • 137

    KelempisiotiA.EskolaP. J.OkuloffA.KarjalainenU.TakataloJ.DaavittilaI.et al (2011). Genetic Susceptibility of Intervertebral Disc Degeneration Among Young Finnish Adults. Bmc. Med. Genet.12, 12153. 10.1186/1471-2350-12-153

  • 138

    KellyT.-A. N.NgK. W.AteshianG. A.HungC. T. (2009). Analysis of Radial Variations in Material Properties and Matrix Composition of Chondrocyte-Seeded Agarose Hydrogel Constructs. Osteoarthr. Cartil.17, 7382. 10.1016/j.joca.2008.05.019

  • 139

    KeplerC. K.PonnappanR. K.TannouryC. A.RisbudM. V.AndersonD. G. (2013). The Molecular Basis of Intervertebral Disc Degeneration. Spine J.13, 318330. 10.1016/j.spinee.2012.12.003

  • 140

    KimA.-G.KimT.-W.KwonW.-K.LeeK.-H.JeongS.HwangM.-H.et al (2021). Microfluidic Chip with Low Constant-Current Stimulation (LCCS) Platform: Human Nucleus Pulposus Degeneration In Vitro Model for Symptomatic Intervertebral Disc. Micromachines12, 1291. 10.3390/mi12111291

  • 141

    KimN. K.ShinD. A.HanI. B.YooE. H.KimS. H.ChungS. S. (2010). The Association of Aggrecan Gene Polymorphism with the Risk of Intervertebral Disc Degeneration. Acta Neurochir.153, 129133. 10.1007/s00701-010-0831-2

  • 142

    KnudsonC. B.KnudsonW. (2001). Cartilage Proteoglycans. Seminars Cell & Dev. Biol.12, 6978. 10.1006/scdb.2000.0243

  • 143

    KrockE.RosenzweigD. H.Chabot-DoréA. J.JarzemP.WeberM. H.OuelletJ. A.et al (2014). Painful, Degenerating Intervertebral Discs Up-Regulate Neurite Sprouting andCGRPthrough Nociceptive Factors. J. Cell. Mol. Med.18, 12131225. 10.1111/jcmm.12268

  • 144

    KudoS.MizunoK.HiraiY.ShimizuT. (1990). Clearance and Tissue Distribution of Recombinant Human Interleukin 1 Beta in Rats. Cancer. Res.50, 57515755.

  • 145

    LakstinsK.ArnoldL.GunschG.FlaniganD.KhanS.GaddeN.et al (2021). Characterization of the Human Intervertebral Disc Cartilage Endplate at the Molecular, Cell, and Tissue Levels. J. Orthop. Res.39, 18981907. 10.1002/jor.24854

  • 146

    LamaP.Le MaitreC. L.HardingI. J.DolanP.AdamsM. A. (2018). Nerves and Blood Vessels in Degenerated Intervertebral Discs Are Confined to Physically Disrupted Tissue. J. Anat.233, 8697. 10.1111/joa.12817

  • 147

    LarsonJ. W.LevicoffE. A.GilbertsonL. G.KangJ. D. (2006). Biologic Modification of Animal Models of Intervertebral Disc Degeneration. J. Bone. Jt. Surg. Am.88 (Suppl. 2), 8387. 10.2106/JBJS.F.00043

  • 148

    LazarusE.Bermudez-LekerikaP.FarchioneD.SchofieldT.HowardS.MambetkadyrovI.et al (2021). Sulfated Hydrogels in Intervertebral Disc and Cartilage Research. Cells10, 3568. 10.3390/cells10123568

  • 149

    Le MaitreC.FreemontA. J.HoylandJ. (2005). The Role of Interleukin-1 in the Pathogenesis of Human Intervertebral Disc Degeneration. Arthritis Res. Ther.7, R732R745. 10.1186/ar1732

  • 150

    Le MaitreC.HoylandJ.FreemontA. J. (2007). Catabolic Cytokine Expression in Degenerate and Herniated Human Intervertebral Discs: IL-1β and TNFα Expression Profile. Arthritis Res. Ther.9, R77. 10.1186/ar2275

  • 151

    Le MaitreC. L.FrainJ.FotheringhamA. P.FreemontA. J.HoylandJ. A. (2008). Human Cells Derived from Degenerate Intervertebral Discs Respond Differently to Those Derived from Non-degenerate Intervertebral Discs Following Application of Dynamic Hydrostatic Pressure. Biorheology45, 563575. 10.3233/BIR-2008-0498

  • 152

    Le MaitreC. L.FrainJ.Millward-SadlerJ.FotheringhamA. P.FreemontA. J.HoylandJ. A. (2009a). Altered Integrin Mechanotransduction in Human Nucleus Pulposus Cells Derived from Degenerated Discs. Arthritis. Rheum.60(2), 460469. 10.1002/art.24248

  • 153

    Le MaitreC. L.FreemontA. J.HoylandJ. A. (2009b). Expression of Cartilage-Derived Morphogenetic Protein in Human Intervertebral Discs and its Effect on Matrix Synthesis in Degenerate Human Nucleus Pulposus Cells. Arthritis Res. Ther.11(5), R137. 10.1186/ar2808

  • 154

    LeeJ. M.SongJ. Y.BaekM.JungH.-Y.KangH.HanI. B.et al (2011). Interleukin-1β Induces Angiogenesis and Innervation in Human Intervertebral Disc Degeneration. J. Orthop. Res.29, 265269. 10.1002/jor.21210

  • 155

    LeungL.CahillC. M. (2010). TNF-α and Neuropathic Pain - a Review. J. Neuroinflammation7, 2711. 10.1186/1742-2094-7-27

  • 156

    LiC.-G.LiangQ.-Q.ZhouQ.MengaE.CuiX.-J.ShuB.et al (2009). A Continuous Observation of the Degenerative Process in the Intervertebral Disc of Smad3 Gene Knock-Out Mice. Spine34, 13631369. 10.1097/BRS.0b013e3181a3c7c7

  • 157

    LiH.LiW.ZhangL.HeJ.TangL.LiZ.et al (2022). Comprehensive Network Analysis Identified SIRT7, NTRK2, and CHI3L1 as New Potential Markers for Intervertebral Disc Degeneration. J. Oncol.2022, 117. 10.1155/2022/4407541

  • 158

    LiW.LaiK.ChopraN.ZhengZ.DasA.DiwanA. D. (2022). Gut-disc axis: A Cause of Intervertebral Disc Degeneration and Low Back Pain?Eur. Spine J.31(4), 917992. 10.1007/s00586-022-07152-8

  • 159

    LiX.LinF.WuY.LiuN.WangJ.ChenR.et al (2019). Resveratrol Attenuates Inflammation Environment-Induced Nucleus Pulposus Cell Senescence In Vitro. Biosci. Rep.39, 30962260. 10.1042/BSR20190126

  • 160

    LiY.LiK.MaoL.HanX.ZhangK.ZhaoC.et al (2016). Cordycepin Inhibits LPS-Induced Inflammatory and Matrix Degradation in the Intervertebral Disc. PeerJ4, e1992. 10.7717/peerj.1992

  • 161

    LiZ.ShaoZ.ChenS.HuangD.PengY.ChenS.et al (2020). TIGAR Impedes Compression-Induced Intervertebral Disc Degeneration by Suppressing Nucleus Pulposus Cell Apoptosis and Autophagy. J. Cell. Physiology235, 17801794. 10.1002/jcp.29097

  • 162

    LinX.LinQ. (2020). MiRNA-495-3p Attenuates TNF-α Induced Apoptosis and Inflammation in Human Nucleus Pulposus Cells by Targeting IL5RA. Inflammation43, 17971805. 10.1007/s10753-020-01254-5

  • 163

    LindblomK. (1957). Intervertebral-disc Degeneration Considered as a Pressure Atrophy. J. Bone & Jt. Surg.39, 933945. 10.2106/00004623-195739040-00016

  • 164

    LiuG.-Z.IshiharaH.OsadaR.KimuraT.TsujiH. (2001). Nitric Oxide Mediates the Change of Proteoglycan Synthesis in the Human Lumbar Intervertebral Disc in Response to Hydrostatic Pressure. Spine26, 134141. 10.1097/00007632-200101150-00005

  • 165

    LiuL.HeJ.LiuC.YangM.FuJ.YiJ.et al (2021). Cartilage Intermediate Layer Protein Affects the Progression of Intervertebral Disc Degeneration by Regulating the Extracellular Microenvironment (Review). Int. J. Mol. Med.47, 475484. 10.3892/ijmm.2020.4832

  • 166

    LiuS.WuN.LiuJ.LiuH.SuX.LiuZ.et al (2016). Association betweenADAMTS-4gene Polymorphism and Lumbar Disc Degeneration in Chinese Han Population. J. Orthop. Res.34, 860864. 10.1002/jor.23081

  • 167

    LoiblM.Wuertz-KozakK.VadalaG.LangS.FairbankJ.UrbanJ. P. (2019). Controversies in Regenerative Medicine: Should Intervertebral Disc Degeneration Be Treated with Mesenchymal Stem Cells?JOR Spine2, e1043. 10.1002/jsp2.1043

  • 168

    LotzJ. C. (2004). Animal Models of Intervertebral Disc Degeneration. Spine29, 27422750. 10.1097/01.brs.0000146498.04628.f9

  • 169

    D.-S.ShonoY.OdaI.AbumiK.KanedaK. (1997). Effects of Chondroitinase ABC and Chymopapain on Spinal Motion Segment Biomechanics. Spine22, 18281834. ; discussion 1834-5. 10.1097/00007632-199708150-00006

  • 170

    LyonsG.EisensteinS. M.SweetM. B. E. (1981). Biochemical Changes in Intervertebral Disc Degeneration. Biochimica Biophysica Acta (BBA) - General Subj.673, 443453. 10.1016/0304-4165(81)90476-1

  • 171

    MacLeanM. A.TouchetteC. J.HanJ. H.ChristieS. D.PickettG. E. (2020). Gender Differences in the Surgical Management of Lumbar Degenerative Disease: a Scoping Review. Spine32, 799816. 10.3171/2019.11.SPINE19896

  • 172

    MaidhofR.JacobsenT.PapatheodorouA.ChahineN. O. (2014). Inflammation Induces Irreversible Biophysical Changes in Isolated Nucleus Pulposus Cells. PLoS ONE9, e99621. 10.1371/journal.pone.0099621

  • 173

    MainardiA.CambriaE.OcchettaP.MartinI.BarberoA.SchärenS.et al (2021). Intervertebral Disc-On-A-Chip as Advanced In Vitro Model for Mechanobiology Research and Drug Testing: A Review and Perspective. Front. Bioeng. Biotechnol.9, 826867. 10.3389/fbioe.2021.826867

  • 174

    MalandrinoA.NoaillyJ.LacroixD. (2011). The Effect of Sustained Compression on Oxygen Metabolic Transport in the Intervertebral Disc Decreases with Degenerative Changes. PLoS. Comput. Biol.7, e1002112. 10.1371/journal.pcbi.1002112

  • 175

    ManolioT. A.CollinsF. S.CoxN. J.GoldsteinD. B.HindorffL. A.HunterD. J.et al (2009). Finding the Missing Heritability of Complex Diseases. Nature461, 747753. 10.1038/nature08494

  • 176

    MasudaK.LotzJ. C. (2010). New Challenges for Intervertebral Disc Treatment Using Regenerative Medicine. Tissue Eng. Part B Rev.16, 147158. 10.1089/ten.TEB.2009.0451

  • 177

    MayerJ. E.IatridisJ. C.ChanD.QureshiS. A.GottesmanO.HechtA. C. (2013). Genetic Polymorphisms Associated with Intervertebral Disc Degeneration. Spine J.13, 299317. 10.1016/j.spinee.2013.01.041

  • 178

    McAfeeP. C.CunninghamB.HolsappleG.AdamsK.BlumenthalS.GuyerR. D.et al (2005). A Prospective, Randomized, Multicenter Food and Drug Administration Investigational Device Exemption Study of Lumbar Total Disc Replacement with the CHARITÉ Artificial Disc versus Lumbar Fusion. Spine30, 15761583. 10.1097/01.brs.0000170561.25636.1c

  • 179

    McLarenW.GilL.HuntS. E.RiatH. S.RitchieG. R. S.ThormannA.et al (2016). The Ensembl Variant Effect Predictor. Genome. Biol.17, 114. 10.1186/s13059-016-0974-4

  • 180

    MerceronC.MangiaviniL.RoblingA.WilsonT. L.GiacciaA. J.ShapiroI. M.et al (2014). Loss of HIF-1α in the Notochord Results in Cell Death and Complete Disappearance of the Nucleus Pulposus. PLoS. ONE9, e110768. 10.1371/journal.pone.0110768

  • 181

    MetzcarJ.WangY.HeilandR.MacklinP. (2019). A Review of Cell-Based Computational Modeling in Cancer Biology. JCO Clin. cancer Inf.3, 113. 10.1200/CCI.18.00069

  • 182

    MilenkovićT.LaiJ.PržuljN. (2008). GraphCrunch: a Tool for Large Network Analyses. Bmc. Bioinforma.9, 70. 10.1186/1471-2105-9-70

  • 183

    MillecampsI.TajerianM.SageE. H.StoneL. S. (2011). Behavioral Signs of Chronic Back Pain in the SPARC-Null Mouse. Spine36, 95102. 10.1097/BRS.0b013e3181cd9d75

  • 184

    MillecampsM.TajerianM.NasoL.SageH. E.StoneL. S. (2012). Lumbar Intervertebral Disc Degeneration Associated with Axial and Radiating Low Back Pain in Ageing SPARC-Null Mice. Pain153, 11671179. 10.1016/j.pain.2012.01.027

  • 185

    Millward-SadlerS. J.CostelloP. W.FreemontA. J.HoylandJ. A. (2009). Regulation of Catabolic Gene Expression in Normal and Degenerate Human Intervertebral Disc Cells: Implications for the Pathogenesis of Intervertebral Disc Degeneration. Arthritis Res. Ther.11, R65. 10.1186/ar2693

  • 186

    MinS.-K.NakazatoK.YamamotoY.GushikenK.FujimotoH.FujishiroH.et al (2010). Cartilage Intermediate Layer Protein Gene Is Associated with Lumbar Disc Degeneration in Male, but Not Female, Collegiate Athletes. Am. J. Sports. Med.38, 25522557. 10.1177/0363546510376714

  • 187

    MioF.ChibaK.HiroseY.KawaguchiY.MikamiY.OyaT.et al (2007). A Functional Polymorphism in COL11A1, Which Encodes the α1 Chain of Type XI Collagen, Is Associated with Susceptibility to Lumbar Disc Herniation. Am. J. Hum. Genet.81, 12711277. 10.1086/522377

  • 188

    MiyagiM.MillecampsM.DancoA. T.OhtoriS.TakahashiK.StoneL. S. (2014). ISSLS Prize Winner. Spine39, 13451354. 10.1097/BRS.0000000000000334

  • 189

    MiyazakiT.KobayashiS.TakenoK.MeirA.UrbanJ.BabaH. (2009). A Phenotypic Comparison of Proteoglycan Production of Intervertebral Disc Cells Isolated from Rats, Rabbits, and Bovine Tails; Which Animal Model Is Most Suitable to Study Tissue Engineering and Biological Repair of Human Disc Disorders?Tissue Eng. Part A15, 38353846. 10.1089/ten.tea.2009.0250

  • 190

    MoonS. M.YoderJ. H.WrightA. C.SmithL. J.VresilovicE. J.ElliottD. M. (2013). Evaluation of Intervertebral Disc Cartilaginous Endplate Structure Using Magnetic Resonance Imaging. Eur. Spine J.22, 18201828. 10.1007/s00586-013-2798-1

  • 191

    MooreR. J. (2000). The Vertebral End-Plate: what Do We Know?Eur. Spine J.9, 9296. 10.1007/s005860050217

  • 192

    MountjoyE.SchmidtE. M.CarmonaM.SchwartzentruberJ.PeatG.MirandaA.et al (2021). An Open Approach to Systematically Prioritize Causal Variants and Genes at All Published Human GWAS Trait-Associated Loci. Nat. Genet.53, 15271533. 10.1038/s41588-021-00945-5

  • 193

    MuJ.GeW.ZuoX.ChenY.HuangC. (2013). Analysis of Association between IL-1β, CASP-9, and GDF5 Variants and Low-Back Pain in Chinese Male Soldiers. Spi19, 243247. 10.3171/2013.4.SPINE12782

  • 194

    MunirS.FreidinM. B.RadeM.MäättäJ.LivshitsG.WilliamsF. M. K. (2018). Endplate Defect Is Heritable, Associated with Low Back Pain and Triggers Intervertebral Disc Degeneration. Spine43, 14961501. 10.1097/BRS.0000000000002721

  • 195

    NäkkiA.VidemanT.KujalaU. M.SuhonenM.MännikköM.PeltonenL.et al (2011). Candidate Gene Association Study of Magnetic Resonance Imaging-Based Hip Osteoarthritis (OA): Evidence for COL9A2 Gene as a Common Predisposing Factor for Hip OA and Lumbar Disc Degeneration. J. Rheumatol.38, 747752. 10.3899/jrheum.100080

  • 196

    NaqviS. M.BuckleyC. T. (2015). Differential Response of Encapsulated Nucleus Pulposus and Bone Marrow Stem Cells in Isolation and Coculture in Alginate and Chitosan Hydrogels. Tissue Eng. Part A21, 288299. 10.1089/ten.TEA.2013.0719

  • 197

    Neidlinger-WilkeC.MietschA.RinklerC.WilkeH.-J.IgnatiusA.UrbanJ. (2012). Interactions of Environmental Conditions and Mechanical Loads Have Influence on Matrix Turnover by Nucleus Pulposus Cells. J. Orthop. Res.30, 112121. 10.1002/jor.21481

  • 198

    NewellN.LittleJ.ChristouA.AdamsM.AdamC.MasourosS. (2017). Biomechanics of the Human Intervertebral Disc: A Review of Testing Techniques and Results. J. Mech. Behav. Biomed. Mater.69, 420434. 10.1016/j.jmbbm.2017.01.037

  • 199

    NicolaeD. L.GamazonE.ZhangW.DuanS.DolanM. E.CoxN. J. (2010). Trait-associated SNPs Are More Likely to Be eQTLs: Annotation to Enhance Discovery from GWAS. PLoS Genet.6, e1000888. 10.1371/journal.pgen.1000888

  • 200

    NiemeyerF.GalbuseraF.TaoY.KienleA.BeerM.WilkeH.-J. (2021). A Deep Learning Model for the Accurate and Reliable Classification of Disc Degeneration Based on MRI Data. Invest. Radiol.56, 7885. 10.1097/RLI.0000000000000709

  • 201

    Noponen-HietalaN.KyllönenE.MännikköM.IlkkoE.KarppinenJ.OttJ.et al (2003). Sequence Variations in the Collagen IX and XI Genes Are Associated with Degenerative Lumbar Spinal Stenosis. Ann. Rheumatic Dis.62, 12081214. 10.1136/ard.2003.008334

  • 202

    Noponen-HietalaN.VirtanenI.KarttunenR.SchwenkeS.JakkulaE.LiH.et al (2005). Genetic Variations in IL6 Associate with Intervertebral Disc Disease Characterized by Sciatica. Pain114, 186194. 10.1016/j.pain.2004.12.015

  • 203

    NovaisE. J.ChoiH.MadhuV.SuyamaK.AnjoS. I.ManadasB.et al (2021). Hypoxia and Hypoxia-Inducible Factor-1α Regulate Endoplasmic Reticulum Stress in Nucleus Pulposus Cells. Am. J. Pathology191, 487502. 10.1016/j.ajpath.2020.11.012

  • 204

    O′ConnellG. D.VresilovicE. J.ElliottD. M. (2007). Comparison of Animals Used in Disc Research to Human Lumbar Disc Geometry. Spine32, 328333. 10.1097/01.brs.0000253961.40910.c1

  • 205

    OliverJ. C.BlandL. A.OettingerC. W.ArduinoM. J.McAllisterS. K.AgueroS. M.et al (1993). Cytokine Kinetics in an In Vitro Whole Blood Model Following an Endotoxin Challenge. Lymphokine Cytokine Res.12, 115120.

  • 206

    OmairA.HoldenM.LieB. A.ReikerasO.BroxJ. I. (2013). Treatment Outcome of Chronic Low Back Pain and Radiographic Lumbar Disc Degeneration Are Associated with Inflammatory and Matrix Degrading Gene Variants: a Prospective Genetic Association Study. BMC Musculoskelet. Disord.14, 19. 10.1186/1471-2474-14-105

  • 207

    OshimaH.IshiharaH.UrbanJ. P. G.TsujiH. (1993). The Use of Coccygeal Discs to Study Intervertebral Disc Metabolism. J. Orthop. Res.11, 332338. 10.1002/jor.1100110304

  • 208

    OswaldK. A. C.BigdonS. F.CroftA. S.Bermudez-LekerikaP.BergadanoA.GantenbeinB.et al (2021). Establishment of a Novel Method for Spinal Discectomy Surgery in Elderly Rats in an In Vivo Spinal Fusion Model. MPs4, 79. 10.3390/mps4040079

  • 209

    ÖztürkE.ArlovAkselØ.AkselS.LiL.OrnitzD. M.Skjåk-BraekG.et al (2016). Sulfated Hydrogel Matrices Direct Mitogenicity and Maintenance of Chondrocyte Phenotype through Activation of FGF Signaling. Adv. Funct. Mat.26, 36493662. 10.1002/adfm.201600092

  • 210

    PaassiltaP.LohinivaJ.GöringH. H.PeräläM.RäinäS. S.KarppinenJ.et al (2001). Identification of a Novel Common Genetic Risk Factor for Lumbar Disk Disease. JAMA285, 18431849. 10.1001/jama.285.14.1843

  • 211

    ParenteauC. S.LauE. C.CampbellI. C.CourtneyA. (2021). Prevalence of Spine Degeneration Diagnosis by Type, Age, Gender, and Obesity Using Medicare Data. Sci. Rep.11, 5389. 10.1038/s41598-021-84724-6

  • 212

    PattappaG.PeroglioM.PeroglioM.SakaiD.MochidaJ.BennekerL.et al (2014). CCL5/RANTES Is a Key Chemoattractant Released by Degenerative Intervertebral Discs in Organ Culture. eCM27, 124136. ; discussion 136. 10.22203/ecm.v027a10

  • 213

    PaulC. P. L.SchoorlT.ZuiderbaanH. A.Zandieh DoulabiB.van der VeenA. J.van de VenP. M.et al (2013). Dynamic and Static Overloading Induce Early Degenerative Processes in Caprine Lumbar Intervertebral Discs. PLoS ONE8, e62411. 10.1371/journal.pone.0062411

  • 214

    PaulC. P. L.ZuiderbaanH. A.Zandieh DoulabiB.van der VeenA. J.van de VenP. M.SmitT. H.et al (2012). Simulated-physiological Loading Conditions Preserve Biological and Mechanical Properties of Caprine Lumbar Intervertebral Discs in Ex Vivo Culture. PLoS ONE7, e33147. 10.1371/journal.pone.0033147

  • 215

    PearceD.FischerS.HudaF.VahdatiA. (2020). Applications of Computer Modeling and Simulation in Cartilage Tissue Engineering. Tissue. Eng. Regen. Med.17, 113. 10.1007/s13770-019-00216-9

  • 216

    PfaffT.FortunatoM.Sanchez-GonzalezA.BattagliaP. W. (2020). Learning Mesh-Based Simulation with Graph Networks. 10.48550/ARXIV.2010.03409

  • 217

    PfannkucheJ.-J.GuoW.CuiS.MaJ.LangG.PeroglioM.et al (2020). Intervertebral Disc Organ Culture for the Investigation of Disc Pathology and Regeneration - Benefits, Limitations, and Future Directions of Bioreactors. Connect. Tissue Res.61, 304321. 10.1080/03008207.2019.1665652

  • 218

    PhillipsF. M.ReubenJ.WetzelF. T. (2002). Intervertebral Disc Degeneration Adjacent to a Lumbar Fusion. An Experimental Rabbit Model. J. Bone Jt. Surg. Br. volume84, 289294. 10.1302/0301-620x.84b2.11937

  • 219

    PhillipsK. L. E.ChivertonN.MichaelA. L.ColeA. A.BreakwellL. M.HaddockG.et al (2013a). The Cytokine and Chemokine Expression Profile of Nucleus Pulposus Cells: Implications for Degeneration and Regeneration of the Intervertebral Disc. Arthritis Res. Ther.15, R213. 10.1186/ar4408

  • 220

    PhillipsK. L. E.CullenK.ChivertonN.MichaelA. L. R.ColeA. A.BreakwellL. M.et al (2015). Potential Roles of Cytokines and Chemokines in Human Intervertebral Disc Degeneration: Interleukin-1 Is a Master Regulator of Catabolic Processes. Osteoarthr. Cartil.23, 11651177. 10.1016/j.joca.2015.02.017

  • 221

    PhillipsK. L. E.Jordan-MahyN.NicklinM. J. H.Le MaitreC. L. (2013b). Interleukin-1 Receptor Antagonist Deficient Mice Provide Insights into Pathogenesis of Human Intervertebral Disc Degeneration. Ann. Rheum. Dis.72(11), 18601867. 10.1136/annrheumdis-2012-202266

  • 222

    PincusT.BurtonA. K.VogelS.FieldA. P. (2002). A Systematic Review of Psychological Factors as Predictors of Chronicity/disability in Prospective Cohorts of Low Back Pain. Spine27, E109E120. 10.1097/00007632-200203010-00017

  • 223

    PluijmS. M. F.HwV. E.BravenboerN.UitterlindenA. G.SmitJ. H.PolsH. A.et al (2004). Collagen Type I 1 Sp1 Polymorphism, Osteoporosis, and Intervertebral Disc Degeneration in Older Men and Women. Ann. Rheumatic Dis.63, 7177. 10.1136/ard.2002.002287

  • 224

    PockertA. J.RichardsonS. M.Le MaitreC. L.LyonM.DeakinJ. A.ButtleD. J.et al (2009). Modified Expression of the ADAMTS Enzymes and Tissue Inhibitor of Metalloproteinases 3 during Human Intervertebral Disc Degeneration. Arthritis Rheum.60, 482491. 10.1002/art.24291

  • 225

    PurmessurD.FreemontA. J.HoylandJ. A. (2008). Expression and Regulation of Neurotrophins in the Nondegenerate and Degenerate Human Intervertebral Disc. Arthritis Res. Ther.10, R99. 10.1186/ar2487

  • 226

    RaineE. V. A.ReynardL. N.van de LaarI. M. B. H.Bertoli-AvellaA. M.LoughlinJ. (2014). Identification and Analysis of a SMAD3 Cis-Acting eQTL Operating in Primary Osteoarthritis and in the Aneurysms and Osteoarthritis Syndrome. Osteoarthr. Cartil.22, 698705. 10.1016/j.joca.2014.02.931

  • 227

    RajasekaranS.KannaR. M.SenthilN.RaveendranM.RanjaniV.CheungK. M. C.et al (2014). Genetic Susceptibility of Lumbar Degenerative Disc Disease in Young Indian Adults. Eur. Spine J.24, 19691975. 10.1007/s00586-014-3687-y

  • 228

    RajasekaranS.SoundararajanD. C. R.TangavelC.MuthurajanR.Sri Vijay AnandK. S.MatchadoM. S.et al (2020). Human Intervertebral Discs Harbour a Unique Microbiome and Dysbiosis Determines Health and Disease. Eur. Spine J.29, 16211640. 10.1007/s00586-020-06446-z

  • 229

    RastogiA.ThakoreP.LeungA.BenavidesM.MachadoM.MorschauserM. A.et al (2009). Environmental Regulation of Notochordal Gene Expression in Nucleus Pulposus Cells. J. Cell. Physiol.220, 698705. 10.1002/jcp.21816

  • 230

    RätsepT.MinajevaA.AsserT. (2013). Relationship between Neovascularization and Degenerative Changes in Herniated Lumbar Intervertebral Discs. Eur. Spine J.22, 24742480. 10.1007/s00586-013-2842-1

  • 231

    RijsbergenM. V.van RietbergenB.BarthelemyV.EltesP.LazáryÁ.LacroixD.et al (2018). Comparison of Patient-specific Computational Models vs. Clinical Follow-Up, for Adjacent Segment Disc Degeneration and Bone Remodelling after Spinal Fusion. PLoS ONE13, e0200899. 10.1371/journal.pone.0200899

  • 232

    RikardS. M.AtheyT. L.NelsonA. R.ChristiansenS. L. M.LeeJ.-J.HolmesJ. W.et al (2019). Multiscale Coupling of an Agent-Based Model of Tissue Fibrosis and a Logic-Based Model of Intracellular Signaling. Front. Physiol.10, 1481. 10.3389/fphys.2019.01481

  • 233

    RimD. C. (2016). Quantitative Pfirrmann Disc Degeneration Grading System to Overcome the Limitation of Pfirrmann Disc Degeneration Grade. Korean. J. Spine13, 18. 10.14245/kjs.2016.13.1.1

  • 234

    RisbudM. V.ShapiroI. M. (2013). Role of Cytokines in Intervertebral Disc Degeneration: Pain and Disc Content. Nat. Rev. Rheumatol.10, 4456. 10.1038/nrrheum.2013.160

  • 235

    RobertsS.EvansH.TrivediJ.MenageJ. (2006). Histology and Pathology of the Human Intervertebral Disc. J. Bone. Jt. Surg. Am.88 (Suppl. 2), 1014. 10.2106/JBJS.F.00019

  • 236

    RobertsS.MenageJ.UrbanJ. P. G. (1989). Biochemical and Structural Properties of the Cartilage End-Plate and its Relation to the Intervertebral Disc. Spine14, 166174. 10.1097/00007632-198902000-00005

  • 237

    RobinsonM. R.EnglishG.EnglishG.MoserG.Lloyd-JonesL. R.TriplettM. A.et al (2017). Genotype-covariate Interaction Effects and the Heritability of Adult Body Mass Index. Nat. Genet.49, 11741181. 10.1038/ng.3912

  • 238

    RomanoskiC. E.GlassC. K.StunnenbergH. G.WilsonL.AlmouzniG. (2015). Roadmap for Regulation. Nature518, 314316. 10.1038/518314a

  • 239

    RosiakP.LatanskaI.PaulP.SujkaW.KolesinskaB. (2021). Modification of Alginates to Modulate Their Physic-Chemical Properties and Obtain Biomaterials with Different Functional Properties. Molecules26, 7264. 10.3390/molecules26237264

  • 240

    Ruiz WillsC.FoataB.González BallesterM. Á.KarppinenJ.NoaillyJ. (2018). Theoretical Explorations Generate New Hypotheses about the Role of the Cartilage Endplate in Early Intervertebral Disk Degeneration. Front. Physiol.9, 1210. 10.3389/fphys.2018.01210

  • 241

    SadowskaA.AltinayB.HitzlW.FergusonS. J.Wuertz-KozakK. (2020). Hypo-Osmotic Loading Induces Expression of IL-6 in Nucleus Pulposus Cells of the Intervertebral Disc Independent of TRPV4 and TRPM7. Front. Pharmacol.11, 952. 10.3389/fphar.2020.00952

  • 242

    SadowskaA.KamedaT.KamedaT.KrupkovaO.Wuertz-KozakK. (2018). Osmosensing, Osmosignalling and Inflammation: How Intervertebral Disc Cells Respond to Altered Osmolarity. eCM36, 231250. 10.22203/eCM.v036a17

  • 243

    SaggeseT.ThambyahA.WadeK.McGlashanS. R. (2020). Differential Response of Bovine Mature Nucleus Pulposus and Notochordal Cells to Hydrostatic Pressure and Glucose Restriction. Cartilage11, 221233. 10.1177/1947603518775795

  • 244

    SalzerE.MouserV. H. M.TryfonidouM. A.ItoK. (2021). A Bovine Nucleus Pulposus Explant Culture Model. J. Orthop. Res[Epub ahead of print]. 10.1002/jor.25226

  • 245

    SamartzisD.KarppinenJ.ChanD.LukK. D. K.CheungK. M. C. (2012). The Association of Lumbar Intervertebral Disc Degeneration on Magnetic Resonance Imaging with Body Mass Index in Overweight and Obese Adults: a Population-Based Study. Arthritis & Rheumatism64, 14881496. 10.1002/art.33462

  • 246

    SambrookP. N.MacGregorA. J.SpectorT. D. (1999). Genetic Influences on Cervical and Lumbar Disc Degeneration: A Magnetic Resonance Imaging Study in Twins. Arthritis & Rheumatism42, 366372. 10.1002/1529-0131(199902)42:2<366::aid-anr20>3.0.co;2-6

  • 247

    SaraviB.LangG.GradS.AliniM.RichardsR. G.SchmalH.et al (2021). A Proinflammatory, Degenerative Organ Culture Model to Simulate Early-Stage Intervertebral Disc Disease. JoVE168. 10.3791/62100

  • 248

    SchaidD. J.ChenW.LarsonN. B. (2018). From Genome-wide Associations to Candidate Causal Variants by Statistical Fine-Mapping. Nat. Rev. Genet.19, 491504. 10.1038/s41576-018-0016-z

  • 249

    Schizophrenia Working Group of the Psychiatric Genomics Consortium (2014). Biological Insights from 108 Schizophrenia-Associated Genetic Loci. Nature511, 421427. 10.1038/nature13595

  • 250

    SchreiberF.SchwöbbermeyerH. (2005). MAVisto: a Tool for the Exploration of Network Motifs. Bioinformatics21, 35723574. 10.1093/bioinformatics/bti556

  • 251

    SekiS.KawaguchiY.ChibaK.MikamiY.KizawaH.OyaT.et al (2005). A Functional SNP in CILP, Encoding Cartilage Intermediate Layer Protein, Is Associated with Susceptibility to Lumbar Disc Disease. Nat. Genet.37, 607612. 10.1038/ng1557

  • 252

    ShannonP.MarkielA.OzierO.BaligaN. S.WangJ. T.RamageD.et al (2003). Cytoscape: a Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res.13, 24982504. 10.1101/gr.1239303

  • 253

    ShowalterB. L.BecksteinJ. C.MartinJ. T.BeattieE. E.OríasA. A. E.SchaerT. P.et al (2012). Comparison of Animal Discs Used in Disc Research to Human Lumbar Disc. Spine37, E900E907. 10.1097/BRS.0b013e31824d911c

  • 254

    SmithL.ChiaroJ. A.ChiaroJ.NerurkarN.CortesD.HoravaS.et al (2011). Nucleus Pulposus Cells Synthesize a Functional Extracellular Matrix and Respond to Inflammatory Cytokine Challenge Following Long-Term Agarose Culture. eCM22, 291301. 10.22203/ecm.v022a22

  • 255

    SnuggsJ.BunningR. A.BunningR.Le MaitreC. (2021). Osmotic Adaptation of Nucleus Pulposus Cells: the Role of Aquaporin 1, Aquaporin 4 and Transient Receptor Potential Vanilloid 4. eCM41, 121141. 10.22203/eCM.v041a09

  • 256

    SolovievaS.KouhiaS.Leino-ArjasP.Ala-KokkoL.LuomaK.RaininkoR.et al (2004). Interleukin 1 Polymorphisms and Intervertebral Disc Degeneration. Epidemiology15, 626633. 10.1097/01.ede.0000135179.04563.35

  • 257

    SolovievaS.LohinivaJ.Leino-ArjasP.RaininkoR.LuomaK.Ala-KokkoL.et al (2006). Intervertebral Disc Degeneration in Relation to the COL9A3 and the IL-1ß Gene Polymorphisms. Eur. Spine J.15, 613619. 10.1007/s00586-005-0988-1

  • 258

    SolovievaS.NoponenN.MännikköM.Leino-ArjasP.LuomaK.RaininkoR.et al (2007). Association between the Aggrecan Gene Variable Number of Tandem Repeats Polymorphism and Intervertebral Disc Degeneration. Spine32, 17001705. 10.1097/BRS.0b013e3180b9ed51

  • 259

    SongY.-Q.CheungK. M. C.HoD. W. H.PoonS. C. S.ChibaK.KawaguchiY.et al (2008). Association of the Asporin D14 Allele with Lumbar-Disc Degeneration in Asians. Am. J. Hum. Genet.82, 744747. 10.1016/j.ajhg.2007.12.017

  • 260

    SongY.-Q.KarasugiT.CheungK. M. C.ChibaK.HoD. W. H.MiyakeA.et al (2013). Lumbar Disc Degeneration Is Linked to a Carbohydrate Sulfotransferase 3 Variant. J. Clin. Invest.123, 49094917. 10.1172/JCI69277

  • 261

    StirlingA.WorthingtonT.RafiqM.LambertP. A.ElliottT. S. (2001). Association between Sciatica and Propionibacterium Acnes. Lancet357, 20242025. 10.1016/S0140-6736(00)05109-6

  • 262

    SunQ.CrowleyC. A.HuangL.WenJ.ChenJ.BaoE. L.et al (2022). From GWAS Variant to Function: A Study of ∼148,000 Variants for Blood Cell Traits. Hum. Genet. Genomics Adv.3, 100063. 10.1016/j.xhgg.2021.100063

  • 263

    SunY.LvM.ZhouL.TamV.LvF.ChanD.et al (2015). Enrichment of Committed Human Nucleus Pulposus Cells Expressing Chondroitin Sulfate Proteoglycans under Alginate Encapsulation. Osteoarthr. Cartil.23, 11941203. 10.1016/j.joca.2015.02.166

  • 264

    SunZ.ZhangM.ZhaoX. H.LiuZ. H.GaoY.SamartzisD.et al (2013). Immune Cascades in Human Intervertebral Disc: the Pros and Cons. Int. J. Clin. Exp. Pathol.6, 10091014.

  • 265

    SunZ.-M.LingM.HuoY.ChangY.LiY.QinH.et al (2011). Caspase 9 Gene Polymorphism and Susceptibility to Lumbar Disc Disease in the Han Population in Northern China. Connect. Tissue Res.52, 198202. 10.3109/03008207.2010.510914

  • 266

    SunZ.-m.MiaoL.ZhangY.-g.MingL. (2009). Association between the -1562 C/T Polymorphism of Matrix Metalloproteinase-9 Gene and Lumbar Disc Disease in the Young Adult Population in North China. Connect. Tissue Res.50, 181185. 10.1080/03008200802585630

  • 267

    SunZ.LiuB.LuoZ.-J. (2020). The Immune Privilege of the Intervertebral Disc: Implications for Intervertebral Disc Degeneration Treatment. Int. J. Med. Sci.17, 685692. 10.7150/ijms.42238

  • 268

    SuyamaK.SakaiD.HirayamaN.NakamuraY.MatsushitaE.TerayamaH.et al (2018). Effects of interleukin-17A in Nucleus Pulposus Cells and its Small-Molecule Inhibitors for Intervertebral Disc Disease. J. Cell. Mol. Med.22, 55395551. 10.1111/jcmm.13828

  • 269

    SwärdL.HellströmM.JacobssonB. O.NymanR.PetersonL. (1991). Disc Degeneration and Associated Abnormalities of the Spine in Elite Gymnasts. Spine16, 437443. 10.1097/00007632-199104000-00009

  • 270

    TaborH. K.RischN. J.MyersR. M. (2002). Candidate-gene Approaches for Studying Complex Genetic Traits: Practical Considerations. Nat. Rev. Genet.3, 391397. 10.1038/nrg796

  • 271

    TakahashiM.HaroH.WakabayashiY.Kawa-uchiT.KomoriH.ShinomiyaK. (2001). The Association of Degeneration of the Intervertebral Disc with 5a/6a Polymorphism in the Promoter of the Human Matrix Metalloproteinase-3 Gene. J. Bone Jt. Surg. Br. volume83, 491495. 10.1302/0301-620x.83b4.11617

  • 272

    TanH.ZhaoJ.JiangJ.RenY. (2012). Association of the Polymorphism of DR4 with the Risk and Severity of Lumbar Disc Degeneration in the Chinese Han Population. Scand. J. Clin. Laboratory Investigation72, 576579. 10.3109/00365513.2012.713176

  • 273

    TanN. Y.KhachigianL. M. (2009). Sp1 Phosphorylation and its Regulation of Gene Transcription. Mol. Cell. Biol.29, 24832488. 10.1128/MCB.01828-08

  • 274

    TapiasA.CiudadC. J.RoninsonI. B.NoéV. (2008). Regulation of Sp1 by Cell Cycle Related Proteins. Cell. Cycle.7, 28562867. 10.4161/cc.7.18.6671

  • 275

    TeraguchiM.YoshimuraN.HashizumeH.MurakiS.YamadaH.MinamideA.et al (2014). Prevalence and Distribution of Intervertebral Disc Degeneration over the Entire Spine in a Population-Based Cohort: the Wakayama Spine Study. Osteoarthr. Cartil.22, 104110. 10.1016/j.joca.2013.10.019

  • 276

    The GTEX consortium (2020). The GTEx Consortium Atlas of Genetic Regulatory Effects across Human Tissues. Science369, 13181330. 10.1126/science.aaz1776

  • 277

    ThorpeA. A.BachF. C.TryfonidouM. A.Le MaitreC. L.MwaleF.DiwanA. D.et al (2018). Leaping the Hurdles in Developing Regenerative Treatments for the Intervertebral Disc from Preclinical to Clinical. JOR. Spine1, e1027. 10.1002/jsp2.1027

  • 278

    TilkeridisC.BeiT.GarantziotisS.StratakisC. A. (2005). Association of a COL1A1 Polymorphism with Lumbar Disc Disease in Young Military Recruits. J. Med. Genet.42, e44. 10.1136/jmg.2005.033225

  • 279

    TolofariS. K.RichardsonS. M.FreemontA. J.HoylandJ. A. (2010). Expression of Semaphorin 3A and its Receptors in the Human Intervertebral Disc: Potential Role in Regulating Neural Ingrowth in the Degenerate Intervertebral Disc. Arthritis Res. Ther.12, R1R11. 10.1186/ar2898

  • 280

    TsaiT.-T.DanielsonK. G.GuttapalliA.OguzE.AlbertT. J.ShapiroI. M.et al (2006). TonEBP/OREBP Is a Regulator of Nucleus Pulposus Cell Function and Survival in the Intervertebral Disc. J. Biol. Chem.281, 2541625424. 10.1074/jbc.M601969200

  • 281

    TurgutM.UsluS.UysalA. l.YurtsevenM. E.UstünH. (2003). Changes in Vascularity of Cartilage Endplate of Degenerated Intervertebral Discs in Response to Melatonin Administration in Rats. Neurosurg. Rev.26, 133138. 10.1007/s10143-003-0259-8

  • 282

    UffelmannE.HuangQ. Q.MunungN. S.de VriesJ.OkadaY.MartinA. R.et al (2021). Genome-wide Association Studies. Nat. Rev. Methods Prim.1, 121. 10.1038/s43586-021-00056-9

  • 283

    UlrichJ. A.LiebenbergE. C.ThuillierD. U.LotzJ. C. (2007). ISSLS Prize Winner: Repeated Disc Injury Causes Persistent Inflammation. Spine32, 28122819. 10.1097/BRS.0b013e31815b9850

  • 284

    UrbanJ. P. G.SmithS.FairbankJ. C. T. (2004). Nutrition of the Intervertebral Disc. Spine29, 27002709. 10.1097/01.brs.0000146499.97948.52

  • 285

    UrbanJ. P. G. (2002). The Role of the Physicochemical Environment in Determining Disc Cell Behaviour. Biochem. Soc. Trans.30, 858863. 10.1042/bst0300858

  • 286

    ValdesA. M.HassettG.HartD. J.SpectorT. D. (2005). Radiographic Progression of Lumbar Spine Disc Degeneration Is Influenced by Variation at Inflammatory Genes. Spine30, 24452451. 10.1097/01.brs.0000184369.79744.a5

  • 287

    van UdenS.Silva-CorreiaJ.OliveiraJ. M.ReisR. L. (2017). Current Strategies for Treatment of Intervertebral Disc Degeneration: Substitution and Regeneration Possibilities. Biomater. Res.21, 22. 10.1186/s40824-017-0106-6

  • 288

    VazirinejadR.AhmadiZ.Kazemi ArababadiM.HassanshahiG.KennedyD. (2014). The Biological Functions, Structure and Sources of CXCL10 and its Outstanding Part in the Pathophysiology of Multiple Sclerosis. Neuroimmunomodulation21, 322330. 10.1159/000357780

  • 289

    VergroesenP.-P. A.KingmaI.EmanuelK. S.HoogendoornR. J. W.WeltingT. J.van RoyenB. J.et al (2015). Mechanics and Biology in Intervertebral Disc Degeneration: a Vicious Circle. Osteoarthr. Cartil.23, 10571070. 10.1016/j.joca.2015.03.028

  • 290

    VergroesenP.-P. A.van der VeenA. J.van RoyenB. J.KingmaI.SmitT. H. (2014). Intradiscal Pressure Depends on Recent Loading and Correlates with Disc Height and Compressive Stiffness. Eur. Spine J.23, 23592368. 10.1007/s00586-014-3450-4

  • 291

    VerrecchiaF.MauvielA. (2004). TGF-β and TNF-α: Antagonistic Cytokines Controlling Type I Collagen Gene Expression. Cell. Signal.16, 873880. 10.1016/j.cellsig.2004.02.007

  • 292

    VidemanT.BattiéM. C. (1999). Spine Update. Spine24, 11641168. 10.1097/00007632-199906010-00020

  • 293

    VidemanT.SaarelaJ.KaprioJ.NäkkiA.LevälahtiE.GillK.et al (2009). Associations of 25 Structural, Degradative, and Inflammatory Candidate Genes with Lumbar Disc Desiccation, Bulging, and Height Narrowing. Arthritis Rheum.60, 470481. 10.1002/art.24268

  • 294

    VirgilioK. M.MartinK. S.PeirceS. M.BlemkerS. S. (2018). Agent-based Model Illustrates the Role of the Microenvironment in Regeneration in Healthy and Mdx Skeletal Muscle. J. Appl. Physiology125, 14241439. 10.1152/japplphysiol.00379.2018

  • 295

    VirtanenI. M.KarppinenJ.TaimelaS.OttJ.BarralS.KaikkonenK.et al (2007). Occupational and Genetic Risk Factors Associated with Intervertebral Disc Disease. Spine32, 11291134. 10.1097/01.brs.0000261473.03274.5c

  • 296

    VitalJ.-M.BoissièreL. (2014). Total Disc Replacement. Orthop. Traumatology Surg. Res.100, S1S14. 10.1016/j.otsr.2013.06.018

  • 297

    VolzM.ElmasryS.JacksonA. R.TravascioF. (2022). Computational Modeling Intervertebral Disc Pathophysiology: A Review. Front. Physiol.12, 750668. 10.3389/fphys.2021.750668

  • 298

    WalterB. A.Illien-JüngerS.NasserP. R.HechtA. C.IatridisJ. C. (2014). Development and Validation of a Bioreactor System for Dynamic Loading and Mechanical Characterization of Whole Human Intervertebral Discs in Organ Culture. J. Biomechanics47, 20952101. 10.1016/j.jbiomech.2014.03.015

  • 299

    WalterB. A.KoreckiC. L.PurmessurD.RoughleyP. J.MichalekA. J.IatridisJ. C. (2011). Complex Loading Affects Intervertebral Disc Mechanics and Biology. Osteoarthr. Cartil.19, 10111018. 10.1016/j.joca.2011.04.005

  • 300

    WalterB.PurmessurD.PurmessurD.MoonA.OcchiogrossoJ.LaudierD.et al (2016). Reduced Tissue Osmolarity Increases TRPV4 Expression and Pro-inflammatory Cytokines in Intervertebral Disc Cells. eCM32, 123136. 10.22203/ecm.v032a08

  • 301

    WangH.ZhouY.HuangB.LiuL.-T.LiuM.-H.WangJ.et al (2014). Utilization of Stem Cells in Alginate for Nucleus Pulposus Tissue Engineering. Tissue Eng. Part A20, 908920. 10.1089/ten.TEA.2012.0703

  • 302

    WangM.TangD.ShuB.WangB.JinH.HaoS.et al (2012). Conditional Activation of β-catenin Signaling in Mice Leads to Severe Defects in Intervertebral Disc Tissue. Arthritis & Rheumatism64, 26112623. 10.1002/art.34469

  • 303

    WangS.-L.YuY.-L.TangC.-L.LvF.-Z. (2013). Effects of TGF-Β1 and IL-1β on Expression of ADAMTS Enzymes and TIMP-3 in Human Intervertebral Disc Degeneration. Exp. Ther. Med.6, 15221526. 10.3892/etm.2013.1348

  • 304

    WángY. X. J.WángJ.-Q.KáplárZ. (2016). Increased Low Back Pain Prevalence in Females Than in Males after Menopause Age: Evidences Based on Synthetic Literature Review. Quant. Imaging Med. Surg.6, 199206. 10.21037/qims.2016.04.06

  • 305

    WanglerS.MenzelU.LiZ.MaJ.HoppeS.BennekerL. M.et al (2019). CD146/MCAM Distinguishes Stem Cell Subpopulations with Distinct Migration and Regenerative Potential in Degenerative Intervertebral Discs. Osteoarthr. Cartil.27, 10941105. 10.1016/j.joca.2019.04.002

  • 306

    WatanabeK.TaskesenE.van BochovenA.PosthumaD. (2017). Functional Mapping and Annotation of Genetic Associations with FUMA. Nat. Commun.8, 111. 10.1038/s41467-017-01261-5

  • 307

    WeilerC.NerlichA. G.BachmeierB. E.BoosN. (2005). Expression and Distribution of Tumor Necrosis Factor Alpha in Human Lumbar Intervertebral Discs: A Study in Surgical Specimen and Autopsy Controls. Spine30, 4453. 10.1097/01.brs.0000149186.63457.20

  • 308

    WernickeS.RascheF. (2006). FANMOD: a Tool for Fast Network Motif Detection. Bioinformatics22, 11521153. 10.1093/bioinformatics/btl038

  • 309

    WilkeH. J.NeefP.CaimiM.HooglandT.ClaesL. E. (1999). New In Vivo Measurements of Pressures in the Intervertebral Disc in Daily Life. Spine24, 755762. 10.1097/00007632-199904150-00005

  • 310

    WilliamsF. M. K.BansalA. T.van MeursJ. B.BellJ. T.MeulenbeltI.SuriP.et al (2013). Novel Genetic Variants Associated with Lumbar Disc Degeneration in Northern Europeans: a Meta-Analysis of 4600 Subjects. Ann. Rheum. Dis.72, 11411148. 10.1136/annrheumdis-2012-201551

  • 311

    WilliamsF. M. K.PophamM.HartD. J.de SchepperE.Bierma-ZeinstraS.HofmanA.et al (2011a). GDF5 Single-Nucleotide Polymorphism Rs143383 Is Associated with Lumbar Disc Degeneration in Northern European Women. Arthritis & Rheumatism63, 708712. 10.1002/art.30169

  • 312

    WilliamsF. M. K.PophamM.SambrookP. N.JonesA. F.SpectorT. D.MacGregorA. J. (2011b). Progression of Lumbar Disc Degeneration over a Decade: a Heritability Study. Ann. Rheumatic Dis.70, 12031207. 10.1136/ard.2010.146001

  • 313

    WilliamsR. J.TryfonidouM. A.SnuggsJ. W.Le MaitreC. L. (2021). Cell Sources Proposed for Nucleus Pulposus Regeneration. JOR. Spine4, e1175. 10.1002/jsp2.1175

  • 314

    WöltjeM.BöbelM. (2017). “Natural Biodegradable Medical Polymers,” in Science and Principles of Biodegradable and Bioresorbable Medical Polymers. Editor ZhangX. (Amsterdam, The Netherlands: Woodhead Publishing), 351376. 10.1016/B978-0-08-100372-5.00012-X

  • 315

    WrocklageC.WassmannH.PaulusW. (2000). COL9A2 Allelotypes in Intervertebral Disc Disease. Biochem. Biophysical Res. Commun.279, 398400. 10.1006/bbrc.2000.3967

  • 316

    WuN.ChenJ.LiuH.ZhaoL.LiuS.LiuJ.et al (2014). The Involvement of ADAMTS-5 Genetic Polymorphisms in Predisposition and Diffusion Tensor Imaging Alterations of Lumbar Disc Degeneration. J. Orthop. Res.32, 686694. 10.1002/jor.22582

  • 317

    WuertzK.UrbanJ. P. G.KlasenJ.IgnatiusA.WilkeH.-J.ClaesL.et al (2007). Influence of Extracellular Osmolarity and Mechanical Stimulation on Gene Expression of Intervertebral Disc Cells. J. Orthop. Res.25, 15131522. 10.1002/jor.20436

  • 318

    XuC.LuoS.WeiL.WuH.GuW.ZhouW.et al (2021). Integrated Transcriptome and Proteome Analyses Identify Novel Regulatory Network of Nucleus Pulposus Cells in Intervertebral Disc Degeneration. Bmc. Med. Genomics.14, 112. 10.1186/s12920-021-00889-z

  • 319

    XuK.WangX.ZhangQ.LiangA.ZhuH.HuangD.et al (2016). Sp1 Downregulates Proinflammatory Cytokine-Induced Catabolic Gene Expression in Nucleus Pulposus Cells. Mol. Med. Rep.14, 39613968. 10.3892/mmr.2016.5730

  • 320

    YangJ.-E.ZhaoK.-H.QuY.ZouY.-C. (2022). Increased Serum CXCL10 Levels Are Associated with Clinical Severity and Radiographic Progression in Patients with Lumbar Disc Degeneration. Clin. Chim. Acta525, 1522. 10.1016/j.cca.2021.12.006

  • 321

    YangJ.FerreiraT.FerreiraT.MorrisA. P.MedlandS. E.MaddenP. A. F.et al (2012). Conditional and Joint Multiple-SNP Analysis of GWAS Summary Statistics Identifies Additional Variants Influencing Complex Traits. Nat. Genet.44, 369375. 10.1038/ng.2213

  • 322

    YaoY.DengQ.SongW.ZhangH.LiY.YangY.et al (2016). MIF Plays a Key Role in Regulating Tissue-specific Chondro-Osteogenic Differentiation Fate of Human Cartilage Endplate Stem Cells under Hypoxia. Stem. Cell. Rep.7, 249262. 10.1016/j.stemcr.2016.07.003

  • 323

    YaoY.SongW.DengQ.ZhangH.WangJ.LiuH.et al (2017). General Regulatory Effects of Hypoxia on Human Cartilage Endplate-Derived Stem Cells: A Genome-wide Analysis of Differential Gene Expression and Alternative Splicing Events. Mol. Med. Rep.16, 30013009. 10.3892/mmr.2017.6907

  • 324

    YeF.LyuF. J.WangH.ZhengZ. (2022). The Involvement of Immune System in Intervertebral Disc Herniation and Degeneration. JOR. Spine5, e1196. 10.1002/jsp2.1196

  • 325

    YoonS.-H.MiyazakiM.HongS.-W.TowB.MorishitaY.HuM.et al (2008). A Porcine Model of Intervertebral Disc Degeneration Induced by Annular Injury Characterized with Magnetic Resonance Imaging and Histopathological Findings. Spi8, 450457. 10.3171/SPI/2008/8/5/450

  • 326

    YuH.ZhuX.GreenbaumD.KarroJ.GersteinM. (2004). TopNet: a Tool for Comparing Biological Sub-networks, Correlating Protein Properties with Topological Statistics. Nucleic Acids Res.32, 328337. 10.1093/nar/gkh164

  • 327

    YuW.FuJ.LiuY.WuY.JiangD. (2018). Osteogenic Protein-1 Inhibits Nucleus Pulposus Cell Apoptosis through Regulating the NF-Κb/ROS Pathway in an Inflammation Environment. Biosci. Rep.38, BSR20181530. 10.1042/BSR20181530

  • 328

    ZamanianC.BhandarkarA. R.MonieD. D.MoinuddinF.VileR. G.Quiñones-HinojosaA.et al (2022). Systems Neuroimmunology: a Review of Multiomics Methodologies to Characterize Neuroimmunological Interactions in Spinal and Cranial Diseases. Neurosurg. Focus.52, E9. 10.3171/2021.11.FOCUS21571

  • 329

    ZeiterS.der WerfM. V.ItoK. (2009). The Fate of Bovine Bone Marrow Stromal Cells in Hydrogels: a Comparison to Nucleus Pulposus Cells and Articular Chondrocytes. J. Tissue. Eng. Regen. Med.3, 310320. 10.1002/term.167

  • 330

    ZhangC.GullbrandS. E.SchaerT. P.LauY. K.JiangZ.DodgeG. R.et al (2020). Inflammatory Cytokine and Catabolic Enzyme Expression in a Goat Model of Intervertebral Disc Degeneration. J. Orthop. Res.38, 25212531. 10.1002/jor.24639

  • 331

    ZhangJ.WangX.LiuH.LiZ.ChenF.WangH.et al (2019). TNF-α Enhances Apoptosis by Promoting Chop Expression in Nucleus Pulposus Cells: Role of the MAPK and NF-Κb Pathways. J. Orthop. Res.37, 697705. 10.1002/jor.24204

  • 332

    ZhangS.DingL.LiX.FanH. (2019). Identification of Biomarkers Associated with the Recurrence of Osteosarcoma Using ceRNA Regulatory Network Analysis. Int. J. Mol. Med.43, 17231733. 10.3892/ijmm.2019.4108

  • 333

    ZhangY. g.SunZ.LiuJ.GuoX. (2008). Advances in Susceptibility Genetics of Intervertebral Degenerative Disc Disease. Int. J. Biol. Sci.4, 283290. 10.7150/ijbs.4.283

  • 334

    ZhangZ.LiF.TianH.GuanK.ZhaoG.ShanJ.et al (2014). Differentiation of Adipose-Derived Stem Cells toward Nucleus Pulposus-like Cells Induced by Hypoxia and a Three-Dimensional Chitosan-Alginate Gel Scaffold In Vitro. Chin. Med. J. Engl.127, 314321. 10.3760/cma.j.issn.0366-6999.20130584

  • 335

    ZhouX.CheungC.-L.KarasugiT.KarppinenJ.SamartzisD.HsuY.-H.et al (2018). Trans-Ethnic Polygenic Analysis Supports Genetic Overlaps of Lumbar Disc Degeneration with Height, Body Mass Index, and Bone Mineral Density. Front. Genet.9, 267. 10.3389/fgene.2018.00267

  • 336

    ZhouZ.CuiS.DuJ.RichardsR. G.AliniM.GradS.et al (2021). One Strike Loading Organ Culture Model to Investigate the Post-traumatic Disc Degenerative Condition. J. Orthop. Transl.26, 141150. 10.1016/j.jot.2020.08.003

  • 337

    ZhuG.-b.JiangX.-r.XiaC.-l.SunY.-j.ZengQ.-s.WuX.-m.et al (2011). Association of FAS and FAS Ligand Polymorphisms with the Susceptibility and Severity of Lumbar Disc Degeneration in Chinese Han Population. Biomarkers16, 485490. 10.3109/1354750X.2011.598563

  • 338

    ZuoR.WangY.LiJ.WuJ.WangW.LiB.et al (2019). Rapamycin Induced Autophagy Inhibits Inflammation-Mediated Endplate Degeneration by Enhancing Nrf2/Keap1 Signaling of Cartilage Endplate Stem Cells. Stem Cells37, 828840. 10.1002/stem.2999

Summary

Keywords

intervertebral disc degeneration, low back pain, inflammation, catabolism, immune-privileged microenvironment, GWAS, artificial intelligence–AI, agent-based model (ABM)

Citation

Bermudez-Lekerika P, Crump KB, Tseranidou S, Nüesch A, Kanelis E, Alminnawi A, Baumgartner L, Muñoz-Moya E, Compte R, Gualdi F, Alexopoulos LG, Geris L, Wuertz-Kozak K, Le Maitre CL, Noailly J and Gantenbein B (2022) Immuno-Modulatory Effects of Intervertebral Disc Cells. Front. Cell Dev. Biol. 10:924692. doi: 10.3389/fcell.2022.924692

Received

20 April 2022

Accepted

20 May 2022

Published

29 June 2022

Volume

10 - 2022

Edited by

Arjun Singh, Memorial Sloan Kettering Cancer Center, United States

Reviewed by

Sumedha Bagga, Boston University, United States

Ipsita Mohanty, Children’s Hospital of Philadelphia, United States

Updates

Copyright

*Correspondence: Benjamin Gantenbein,

†ORCID: Paola Bermudez-Lekerika, orcid.org/0000-0002-6858-8213; Katherine Crump, orcid.org/0000-0001-5328-667X; Sofia Tseranidou, orcid.org/0000-0003-1459-5650; Andrea Nüesch, orcid.org/0000-0002-0962-8705; Exarchos Kanelis, orcid.org/0000-0002-2059-1480; Ahmad Alminnawi, orcid.org/0000-0002-9998-3275; Laura Baumgartner, orcid.org/0000-0002-9801-8318; Estefano Muñoz-Moya, orcid.org/0000-0001-5222-4071; Roger Compte, orcid.org/0000-0002-7034-1092; Francesco Gualdi, orcid.org/0000-0003-0449-9884; Leonidas G. Alexopoulos, orcid.org/0000-0003-0425-166X; Liesbet Geris, orcid.org/0000-0002-8180-1445; Karin Wuertz-Kozak, orcid.org/0000-0002-3281-4629; Christine Le Maitre, orcid.org/0000-0003-4489-7107; Jérôme Noailly, orcid.org/0000-0003-3446-7621; Benjamin Gantenbein, orcid.org/0000-0002-9005-0655

‡These authors have contributed equally to this work and share first authorship

This article was submitted to Molecular and Cellular Pathology, a section of the journal Frontiers in Cell and Developmental Biology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics