Sec. Molecular and Cellular Pathology
Immuno-Modulatory Effects of Intervertebral Disc Cells
- 1Tissue Engineering for Orthopaedics and Mechanobiology, Bone and Joint Program, Department for BioMedical Research (DBMR), Faculty of Medicine, University of Bern, Bern, Switzerland
- 2Department of Orthopaedic Surgery and Traumatology, Inselspital, Bern University Hospital, Medical Faculty, University of Bern, Bern, Switzerland
- 3BCN MedTech, Universitat Pompeu Fabra, Barcelona, Spain
- 4Biomolecular Sciences Research Centre, Sheffield Hallam University, Sheffield, United Kingdom
- 5ProtATonce Ltd., Athens, Greece
- 6School of Mechanical Engineering, National Technical University of Athens, Zografou, Greece
- 7GIGA In Silico Medicine, University of Liège, Liège, Belgium
- 8Skeletal Biology and Engineering Research Center, KU Leuven, Leuven, Belgium
- 9Twin Research and Genetic Epidemiology, St Thomas’ Hospital, King’s College London, London, United Kingdom
- 10Institut Hospital Del Mar D’Investigacions Mèdiques (IMIM), Barcelona, Spain
- 11Biomechanics Research Unit, KU Leuven, Leuven, Belgium
- 12Department of Biomedical Engineering, Rochester Institute of Technology, Rochester, NY, United States
- 13Spine Center, Schön Klinik München Harlaching Academic Teaching Hospital and Spine Research Institute of the Paracelsus Private Medical University Salzburg (Austria), Munich, Germany
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.
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. Comparison of a healthy and a degenerated IVD disc (focused on ECM components). In the intact IVD, the NP matrix mostly contains proteoglycans (PG) and non-oriented collagen type II fibers. The proteoglycans contain negatively charged sulfated groups leading to an intradiscal osmotic pressure crucial for the basal hydration of the NP and the biomechanical function of the IVD. Within the degenerated disc, the total content of PG decreases. Small non-aggregating PGs are present. This drop-in PG content negatively affects the swelling capacity of the disc. Additionally, during disc degeneration, the production of catabolic cytokines, matrix-degrading enzymes, and neurotrophic as well as angiogenic factors occur due to cellular changes. This can lead to blood and nerve vessel ingrowth in the AF. The AF is composed of highly oriented concentric lamella of type I collagen whereas the cell density is higher in intact than in degenerated discs.
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. Comparison of a healthy, degenerated, and herniated IVD discs (focused on cellular involvement). (A) Intact IVD: Native disc cells produce a plethora of cytokines and chemokines expressing the corresponding receptors and maintaining homeostasis in a para and autocrine manner. The CEP is intact with blood vessels. The NP has a high number of proteoglycans. AF cells aligned. (B) Degenerated IVD: Shift to catabolic environment. Cytokines are expressed by disc cells themselves. The CEP has a higher amount of blood vessels than in the intact IVD. Proteoglycan number decreases in the NP. In the AF, there is a loss of alignment and support for AF cells. (C) Herniated IVD with crack in CEP: As soon as the AF or CEP is ruptured during injury or disc degeneration, a route for migration of immune cells into the IVD is provided. Immune cells, including T cells, B cells, macrophages, neutrophils and mast cells, contribute to an inflamed environment within the disc, further increasing the cytokine and chemokine expression and leading to a viscous circle of inflammatory driven catabolism. A crack in CEP allows blood vessels to grow into the AF and NP. The AF herniates/bulges, which is where blood vessel in-growth primarily occurs.
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. Schematic diagram of the different factors contributing to the metabolic shift from anabolism to catabolism in IDD, including genetics and epigenetics, biomechanics, microenvironment, presence of bacteria and other factors. All these contributors can promote a downstream biochemical effects (matrix breakdown and neurotrophins production) leading to structural and biomechanical alterations, nerve ingrowth and blood vessel formation. Thus, involvement of immune system could be achieved by chemotaxis losing the immuno-privileged state of the IVD.
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. 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. Flow chart of the intersection of experiments and computational modeling. (A) First, a literature review is necessary to determine the current state of the research. Then, the researcher can either perform additional experiments to fill gaps of knowledge in the literature, or use published data to create an in-silico model. (B) There are many options for experimental model design, including use of imaging modalities to view the state of the IVD, in-vivo animal studies which better examine the complexity of IDD, bioreactors and microfluidic devices that allow the investigation of mechanical loading or fluid flow in the IVD, and in-vitro/ex-vivo culture of the whole IVD or IVD cells from human or animal tissue. (C) In-silico models or methodologies can use published literature or additional experiments to provide deeper investigations into complex tissue (FEM), cell (ABM), protein (network modeling), and genetic responses (network modeling, GWAS) as well as explore interactions at multiple scales which would be difficult and expensive to do experimentally. These models can help identify novel parameters and interactions that should be validated or explored further through experiments.
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 (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.
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).
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.
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.
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).
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.
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.
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
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.
Adams, M. A., Freeman, B. J. C., Morrison, H. P., Nelson, I. W., and Dolan, P. (2000). Mechanical Initiation of Intervertebral Disc Degeneration. Spine 25, 1625–1636. doi:10.1097/00007632-200007010-00005
Adams, M. A., McMillan, D. W., Green, T. P., and Dolan, P. (1996). Sustained Loading Generates Stress Concentrations in Lumbar Intervertebral Discs. Spine 21, 434–438. doi:10.1097/00007632-199602150-00006
Al Dhaheri, N., Wu, N., Zhao, S., Wu, Z., Blank, R. D., Zhang, J., et al. (2020). KIAA1217 : A Novel Candidate Gene Associated with Isolated and Syndromic Vertebral Malformations. Am. J. Med. Genet. 182, 1664–1672. doi:10.1002/ajmg.a.61607
Albert, H. B., Lambert, P., Rollason, J., Sorensen, J. S., Worthington, T., Pedersen, M. 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, 690–696. doi:10.1007/s00586-013-2674-z
Albert, H. B., Sorensen, J. S., Christensen, B. S., and Manniche, C. (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, 697–707. doi:10.1007/s00586-013-2675-y
Alini, M., Eisenstein, S. M., Ito, K., Little, C., Kettler, A. A., Masuda, K., et al. (2008). Are Animal Models Useful for Studying Human Disc Disorders/degeneration? Eur. Spine J. 17, 2–19. doi:10.1007/s00586-007-0414-y
Alkhatib, B., Rosenzweig, D. H., Rosenzweig, D., Krock, E., Roughley, P., Beckman, L., et al. (2014). Acute Mechanical Injury of the Human Intervertebral Disc: Link to Degeneration and Pain. eCM 28, 98–111. ; discussion 110-1. doi:10.22203/ecm.v028a08
Alvarez-Galvez, J., and Suarez-Lledo, V. (2019). Using Agent-Based Modeling to Understand the Emergence and Reproduction of Social Inequalities in Health. Proceedings 44, 2. doi:10.3390/IECEHS-2-06372
Amini, H., Shroff, N., Stamova, B., Ferino, E., Carmona-Mora, P., Zhan, X., et al. (2020). Genetic Variation Contributes to Gene Expression Response in Ischemic Stroke: an eQTL Study. Ann. Clin. Transl. Neurol. 7, 1648–1660. doi:10.1002/acn3.51154
Anjankar, S. D., Poornima, S., Raju, S., Jaleel, M. A., Bhiladvala, D., and Hasan, Q. (2015). Degenerated Intervertebral Disc Prolapse and its Association of Collagen I Alpha 1 Spl Gene Polymorphism. Ijoo 49, 589–594. doi:10.4103/0019-5413.168765
Aubin, G. G., Baud’huin, M., Lavigne, J.-P., Brion, R., Gouin, F., Lepelletier, D., 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. doi:10.1038/srep42918
Batagelj, V., and Mrvar, A. (2004). “Pajek - Analysis and Visualization of Large Networks,” in Graph Drawing Software. Editor M. M. P. Jünger (Berlin, Heidelberg: Springer Berlin Heidelberg), 77–103. doi:10.1007/978-3-642-18638-7_4
Baumgartner, L., Reagh, J. J., González Ballester, M. A., and Noailly, J. (2021a). Simulating Intervertebral Disc Cell Behaviour within 3D Multifactorial Environments. Bioinformatics 37, 1246–1253. doi:10.1093/bioinformatics/btaa939
Baumgartner, L., Sadowska, A., Tío, L., González Ballester, M. A., Wuertz-Kozak, K., and Noailly, J. (2021b). Evidence-Based Network Modelling to Simulate Nucleus Pulposus Multicellular Activity in Different Nutritional and Pro-inflammatory Environments. Front. Bioeng. Biotechnol. 9, 734258. doi:10.3389/fbioe.2021.734258
Baumgartner, L., Wuertz-Kozak, K., Le Maitre, C. L., Wignall, F., Richardson, S. M., Hoyland, J., et al. (2021c). Multiscale Regulation of the Intervertebral Disc: Achievements in Experimental, In Silico, and Regenerative Research. Ijms 22, 703. doi:10.3390/ijms22020703
Beckstein, J. C., Sen, S., Schaer, T. P., Vresilovic, E. J., and Elliott, D. M. (2008). Comparison of Animal Discs Used in Disc Research to Human Lumbar Disc. Spine 33, E166–E173. doi:10.1097/BRS.0b013e318166e001
Bedore, J., Sha, W., McCann, M. R., Liu, S., Leask, A., and Séguin, C. A. (2013). Loss of Notochord-Derived CCN2 Results in Impaired Intervertebral Disc Development and Premature Disc Degeneration. Arthritis & Rheumatism 65, a–n. doi:10.1002/art.38075
Benneker, L. M., Heini, P. F., Alini, M., Anderson, S. E., and Ito, K. (2005). 2004 Young Investigator Award Winner: Vertebral Endplate Marrow Contact Channel Occlusions and Intervertebral Disc Degeneration. Spine 30, 167–173. doi:10.1097/01.brs.0000150833.93248.09
Binch, A. L. A., Cole, A. A., Breakwell, L. M., Michael, A. L. R., Chiverton, N., Creemers, L. B., et al. (2015a). Class 3 Semaphorins Expression and Association with Innervation and Angiogenesis within the Degenerate Human Intervertebral Disc. Oncotarget 6, 18338–18354. doi:10.18632/oncotarget.4274
Binch, A. L. A., Cole, A. A., Breakwell, L. M., Michael, A. L. R., Chiverton, N., Creemers, L. B., et al. (2015b). Nerves Are More Abundant Than Blood Vessels in the Degenerate Human Intervertebral Disc. Arthritis Res. Ther. 17, 370. doi:10.1186/s13075-015-0889-6
Blumenthal, S., McAfee, P. C., Guyer, R. D., Hochschuler, S. H., Geisler, F. H., Holt, R. 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. Spine 30, 1565–1575. discussion E387-91. doi:10.1097/01.brs.0000170587.32676.0e
Binch, A. L., Cole, A. A., Breakwell, L. M., Michael, A. L., Chiverton, N., Cross, A. K., et al. (2014). Expression and Regulation of Neurotrophic and Angiogenic Factors during Human Intervertebral Disc Degeneration. Arthritis Res. Ther. 16, 416. doi:10.1186/s13075-014-0416-1
Bongers, P. M., Hulshof, C. T. J., DlJKSTRA, L., Boshuizen, H. C., Groenhout, H. J. M., and Valken, E. (1990). Back Pain and Exposure to Whole Body Vibration in Helicopter Pilots. Ergonomics 33, 1007–1026. doi:10.1080/00140139008925309
Bournonville, S., Lambrechts, T., Vanhulst, J., Luyten, F. P., Papantoniou, I., and Geris, L. (2019). Towards Self-Regulated Bioprocessing: A Compact Benchtop Bioreactor System for Monitored and Controlled 3D Cell and Tissue Culture. Biotechnol. J. 14, 1800545. doi:10.1002/biot.201800545
Bråten, L. C. H., Grøvle, L., Espeland, A., Pripp, A. H., Grotle, M., Helllum, C., 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, 1–11. doi:10.1186/s12891-020-03422-y
Cambria, E., Brunner, S., Heusser, S., Fisch, P., Hitzl, W., Ferguson, S. J., et al. (2020). Cell-Laden Agarose-Collagen Composite Hydrogels for Mechanotransduction Studies. Front. Bioeng. Biotechnol. 8, 346. doi:10.3389/fbioe.2020.00346
Capoor, M. N., Konieczna, A., McDowell, A., Ruzicka, F., Smrcka, M., Jancalek, R., et al. (2021). Pro-Inflammatory and Neurotrophic Factor Responses of Cells Derived from Degenerative Human Intervertebral Discs to the Opportunistic Pathogen Cutibacterium Acnes. Ijms 22, 2347. doi:10.3390/ijms22052347
Chan, S. C. W., Walser, J., Käppeli, P., Shamsollahi, M. J., Ferguson, S. J., and Gantenbein-Ritter, B. (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. doi:10.1371/journal.pone.0072489
Chang, H., Cai, F., Zhang, Y., Xue, M., Liu, L., Yang, A., et al. (2017). Early-stage Autophagy Protects Nucleus Pulposus Cells from Glucose Deprivation-Induced Degeneration via the P -eIF2α/ATF4 Pathway. Biomed. Pharmacother. 89, 529–535. doi:10.1016/j.biopha.2017.02.074
Chang, H., Wang, H., Yang, X., You, K., Jiang, M., Cai, F., et al. (2021). Comprehensive Profile Analysis of Differentially Expressed circRNAs in Glucose Deprivation-Induced Human Nucleus Pulposus Cell Degeneration. BioMed Res. Int. 2021, 1–14. doi:10.1155/2021/4770792
Chen, J.-W., Li, B., Yang, Y.-H., Jiang, S.-D., and Jiang, L.-S. (2014). Significance of Hypoxia in the Physiological Function of Intervertebral Disc Cells. Crit. Rev. Eukaryot. Gene. Expr. 24, 193–204. doi:10.1615/critreveukaryotgeneexpr.2014010485
Chen, J., Baer, A. E., Paik, P. Y., Yan, W., and Setton, L. A. (2002). Matrix Protein Gene Expression in Intervertebral Disc Cells Subjected to Altered Osmolarity. Biochem. Biophysical Res. Commun. 293, 932–938. doi:10.1016/S0006-291X(02)00314-5
Cheung, K. M. C., Karppinen, J., Chan, D., Ho, D. W. H., Song, Y.-Q., Sham, P., et al. (2009). Prevalence and Pattern of Lumbar Magnetic Resonance Imaging Changes in a Population Study of One Thousand Forty-Three Individuals. Spine 34, 934–940. doi:10.1097/BRS.0b013e3181a01b3f
Ching, C. T. S., Chow, D. H. K., Yao, F. Y. D., and Holmes, A. 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, 182–189. doi:10.1016/s0268-0033(02)00188-2
Choi, H., Chaiyamongkol, W., Doolittle, A. C., Johnson, Z. I., Gogate, S. S., Schoepflin, Z. 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, 8969–8981. doi:10.1074/jbc.RA117.001167
Chou, P.-H., Wang, S.-T., Yen, M.-H., Liu, C.-L., Chang, M.-C., and Lee, O. 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. doi:10.1186/s13287-016-0292-5
Chu, G., Shi, C., Wang, H., Zhang, W., Yang, H., and Li, B. (2018). Strategies for Annulus Fibrosus Regeneration: From Biological Therapies to Tissue Engineering. Front. Bioeng. Biotechnol. 6, 90. doi:10.3389/fbioe.2018.00090
Costi, J. J., Stokes, I. A., Gardner-Morse, M. G., and Iatridis, J. C. (2008). Frequency-Dependent Behavior of the Intervertebral Disc in Response to Each of Six Degree of Freedom Dynamic Loading. Spine 33, 1731–1738. doi:10.1097/BRS.0b013e31817bb116
Croft, A. S., Roth, Y., Oswald, K. A. C., Ćorluka, S., Bermudez-Lekerika, P., and Gantenbein, B. (2021). In Situ Cell Signalling of the Hippo-YAP/TAZ Pathway in Reaction to Complex Dynamic Loading in an Intervertebral Disc Organ Culture. Ijms 22, 13641. doi:10.3390/ijms222413641
Dahia, C. L., Engiles, J. B., Gullbrand, S. E., Iatridis, J. C., Lai, A., Le Maitre, C. L., et al. (2021). A Perspective on the ORS Spine Section Initiative to Develop a Multi-Species JOR Spine Histopathology Series. JOR. Spine 4, e1165. doi:10.1002/jsp2.1165
Dai, J., Xing, Y., Xiao, L., Li, J., Cao, R., He, Y., 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, 2041–2051. doi:10.1021/acsbiomaterials.8b01522
Daly, C., Ghosh, P., Jenkin, G., Oehme, D., and Goldschlager, T. (2016). A Review of Animal Models of Intervertebral Disc Degeneration: Pathophysiology, Regeneration, and Translation to the Clinic. BioMed Res. Int. 2016, 1–14. doi:10.1155/2016/5952165
Dong, D. M., Yao, M., Liu, B., Sun, C. Y., Jiang, Y. Q., and Wang, Y. 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, 1958–1961. doi:10.1007/s00586-007-0454-3
Du, H., Bai, B., Qiu, Y., Yin, S., and Bian, W. (2015). Association between TRAIL Gene Polymorphisms and the Susceptibility and Severity of Lumbar Disc Degeneration. Int. J. Clin. Exp. Pathol. 8, 7415–7420.
Du, Y., Tavana, S., Rahman, T., Baxan, N., Hansen, U. N., and Newell, N. (2021). Sensitivity of Intervertebral Disc Finite Element Models to Internal Geometric and Non-geometric Parameters. Front. Bioeng. Biotechnol. 9, 660013. doi:10.3389/fbioe.2021.660013
Dudli, S., Haschtmann, D., and Ferguson, S. J. (2011). Fracture of the Vertebral Endplates, but Not Equienergetic Impact Load, Promotes Disc Degeneration In Vitro. J. Orthop. Res. 30, 809–816. doi:10.1002/jor.21573
Dudli, S., Haschtmann, D., and Ferguson, S. 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, 1901–1908. doi:10.1007/s00586-014-3301-3
Dudli, S., Liebenberg, E., Magnitsky, S., Lu, B., Lauricella, M., and Lotz, J. C. (2018). Modic Type 1 Change Is an Autoimmune Response that Requires a Proinflammatory Milieu provided by the 'Modic Disc'. Spine J. 18, 831–844. doi:10.1016/j.spinee.2017.12.004
Duncan, L. E., Ostacher, M., and Ballon, J. (2019). How Genome-wide Association Studies (GWAS) Made Traditional Candidate Gene Studies Obsolete. Neuropsychopharmacol. 44, 1518–1523. doi:10.1038/s41386-019-0389-5
Elfervig, M. K., Minchew, J. T., Francke, E., Tsuzaki, M., and Banes, A. J. (2001a). IL-1? Sensitizes Intervertebral Disc Annulus Cells to Fluid-Induced Shear Stress. J. Cell. Biochem. 82, 290–298. doi:10.1002/jcb.1153
Elfervig, M. K., Minchew, J. T., Francke, E., Tsuzaki, M., and Banes, A. J. (2001b). IL-1? Sensitizes Intervertebral Disc Annulus Cells to Fluid-Induced Shear Stress. J. Cell. Biochem. 82, 290–298. doi:10.1002/jcb.1153
Fernando, H. N., Czamanski, J., Yuan, T.-Y., Gu, W., Salahadin, A., and Huang, C.-Y. C. (2011). Mechanical Loading Affects the Energy Metabolism of Intervertebral Disc Cells. J. Orthop. Res. 29, 1634–1641. doi:10.1002/jor.21430
Flannery, C. R., Little, C. B., and Caterson, B. (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, 507–511. doi:10.1016/s0945-053x(98)90021-x
Forgetta, V., Jiang, L., Vulpescu, N. A., Hogan, M. S., Chen, S., Morris, J. A., et al. (2022). An Effector Index to Predict Target Genes at GWAS Loci. Hum. Genet. [Epub ahead of print]. doi:10.1007/s00439-022-02434-z
Frauchiger, D., Tekari, A., Tekari, A., Wöltje, M., Fortunato, G., Benneker, L., et al. (2017). A Review of the Application of Reinforced Hydrogels and Silk as Biomaterials for Intervertebral Disc Repair. eCM 34, 271–290. doi:10.22203/eCM.v034a17
Freemont, A. J., Watkins, A., Le Maitre, C., Baird, P., Jeziorska, M., Knight, M. T. N., et al. (2002). Nerve Growth Factor Expression and Innervation of the Painful Intervertebral Disc. J. Pathol. 197, 286–292. doi:10.1002/path.1108
Freemont, A., Peacock, T., Goupille, P., Hoyland, J., O'Brien, J., and Jayson, M. (1997). Nerve Ingrowth into Diseased Intervertebral Disc in Chronic Back Pain. Lancet 350, 178–181. doi:10.1016/s0140-6736(97)02135-1
Frobin, W., Brinckmann, P., Kramer, M., and Hartwig, E. (2001). Height of Lumbar Discs Measured from Radiographs Compared with Degeneration and Height Classified from MR Images. Eur. Radiol. 11, 263–269. doi:10.1007/s003300000556
Galbusera, F., Schmidt, H., Neidlinger-Wilke, C., and Wilke, H.-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, 729–739. doi:10.1080/10255842.2010.493522
Galbusera, F., Schmidt, H., Noailly, J., Malandrino, A., Lacroix, D., Wilke, H.-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, 1234–1241. doi:10.1016/j.jmbbm.2011.04.008
Gantenbein, B., Frauchiger, D. A., May, R. D., Bakirci, E., Urs, R., and Sibylle, G. (2019). “Developing Bioreactors to Host Joint-Derived Tissues that Require Mechanical Stimulation,” in Encyclopedia of Tissue Engineering and Regenerative Medicine. Editors R. L. Reis, and M. E. Gomes (NY: Elsevier), 261–280. doi:10.1016/b978-0-12-801238-3.65611-8
Gantenbein, B., Calandriello, E., Wuertz-Kozak, K., Benneker, L. M., Keel, M. J., and Chan, S. C. (2014). Activation of Intervertebral Disc Cells by Co-culture with Notochordal Cells, Conditioned Medium and Hypoxia. BMC Musculoskelet. Disord. 15, 422. doi:10.1186/1471-2474-15-422
Gantenbein, B., Grünhagen, T., Lee, C. R., van Donkelaar, C. C., Alini, M., and Ito, K. (2006). An In Vitro Organ Culturing System for Intervertebral Disc Explants with Vertebral Endplates. Spine 31, 2665–2673. doi:10.1097/01.brs.0000244620.15386.df
Gantenbein, B., Illien-Jünger, S., Chan, S., Walser, J., Haglund, L., Ferguson, S., et al. (2015). Organ Culture Bioreactors - Platforms to Study Human Intervertebral Disc Degeneration and Regenerative Therapy. Cscr 10, 339–352. doi:10.2174/1574888X10666150312102948
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. Lancet 392, 1789–1858. doi:10.1016/S0140-6736(18)32279-7
Ghezelbash, F., Schmidt, H., Shirazi-Adl, A., and El-Rich, M. (2020). Internal Load-Sharing in the Human Passive Lumbar Spine: Review of In Vitro and Finite Element Model Studies. J. Biomechanics 102, 109441. doi:10.1016/j.jbiomech.2019.109441
Gilbert, H. T., Hoyland, J. A., Freemont, A. J., and Millward-Sadler, S. 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. doi:10.1186/ar3229
Gilbert, H. T. J., Hodson, N., Baird, P., Richardson, S. M., and Hoyland, J. A. (2016). Acidic pH Promotes Intervertebral Disc Degeneration: Acid-Sensing Ion Channel -3 as a Potential Therapeutic Target. Sci. Rep. 6, 37360. doi:10.1038/srep37360
Gilligan, C. J., Cohen, S. P., Fischetti, V. A., Hirsch, J. A., and Czaplewski, L. G. (2021). Chronic Low Back Pain, Bacterial Infection and Treatment with Antibiotics. Spine J. 21, 903–914. doi:10.1016/j.spinee.2021.02.013
Goddard, M. E., Kemper, K. E., MacLeod, I. M., Chamberlain, A. J., and Hayes, B. J. (2016). Genetics of Complex Traits: Prediction of Phenotype, Identification of Causal Polymorphisms and Genetic Architecture. Proc. R. Soc. B 283, 20160569. doi:10.1098/rspb.2016.0569
Gorth, D. J., Shapiro, I. M., and Risbud, M. 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, 1531–1542. doi:10.1002/jbmr.3714
Gorth, D. J., Shapiro, I. M., and Risbud, M. V. (2018). Transgenic Mice Overexpressing Human TNF-α Experience Early Onset Spontaneous Intervertebral Disc Herniation in the Absence of Overt Degeneration. Cell. death. Dis. 10, 1–14. doi:10.1038/s41419-018-1246-x
Granville Smith, I., Danckert, N. P., Freidin, M. B., Wells, P., Marchesi, J. R., and Williams, F. M. K. (2021). Evidence for Infection in Intervertebral Disc Degeneration: a Systematic Review. Eur. Spine J. 31, 414–430. doi:10.1007/s00586-021-07062-1
Gu, J., Guan, F., Guan, G., Xu, G., Wang, X., Zhao, W., et al. (2013). Aggrecan Variable Number of Tandem Repeat Polymorphism and Lumbar Disc Degeneration. Spine 38, E1600–E1607. doi:10.1097/BRS.0000000000000012
Guerrero, J., Häckel, S., Croft, A. S., Albers, C. E., and Gantenbein, B. (2020). The Effects of 3D Culture on the Expansion and Maintenance of Nucleus Pulposus Progenitor Cell Multipotency. JOR Spine 4, e1131. doi:10.1002/jsp2.1131
Gullbrand, S. E., Malhotra, N. R., Schaer, T. P., Zawacki, Z., Martin, J. T., Bendigo, J. R., et al. (2017). A Large Animal Model that Recapitulates the Spectrum of Human Intervertebral Disc Degeneration. Osteoarthr. Cartil. 25, 146–156. doi:10.1016/j.joca.2016.08.006
Guo, X., and Chen, G. (2020). Hypoxia-Inducible Factor Is Critical for Pathogenesis and Regulation of Immune Cell Functions in Rheumatoid Arthritis. Front. Immunol. 11, 1668. doi:10.3389/fimmu.2020.01668
Haglund, L., Moir, J., Beckman, L., Mulligan, K. R., Jim, B., Ouellet, J. A., et al. (2011). Development of a Bioreactor for Axially Loaded Intervertebral Disc Organ Culture. Tissue Eng. Part C. Methods 17, 1011–1019. doi:10.1089/ten.TEC.2011.0025
Hamilton, D. J., Séguin, C. A., Wang, J., Pilliar, R. M., and Kandel, R. A. (2006). BioEngineering of Skeletal Tissues TeamFormation of a Nucleus Pulposus-Cartilage Endplate Construct In Vitro. Biomaterials 27, 397–405. doi:10.1016/j.biomaterials.2005.07.007
Han, I. B., Ropper, A. E., Teng, Y. D., Shin, D. A., Jeon, Y. J., Park, H. M., et al. (2013). Association between VEGF and eNOS Gene Polymorphisms and Lumbar Disc Degeneration in a Young Korean Population. Genet. Mol. Res. 12, 2294–2305. doi:10.4238/2013.July.8.10
Hellum, C., Johnsen, L. G., JohnsenBerg, L. G., Gjertsen, Ø., Berg, L., Neckelmann, G., 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, 681–690. doi:10.1007/s00586-011-2145-3
Hernandez, P. A., Jacobsen, T. D., Barati, Z., and Chahine, N. O. (2020a). Confocal Scanning of Intervertebral Disc Cells in 3D : Inside Alginate Beads and in Native Microenvironment. JOR Spine 3, e1106. doi:10.1002/jsp2.1106
Hernandez, P. A., Jacobsen, T. D., and Chahine, N. O. (2020b). Actomyosin Contractility Confers Mechanoprotection against TNFα-Induced Disruption of the Intervertebral Disc. Sci. Adv. 6, eaba2368. doi:10.1126/sciadv.aba2368
Hingert, D., Nawilaijaroen, P., Aldridge, J., Baranto, A., and Brisby, H. (2019). Investigation of the Effect of Secreted Factors from Mesenchymal Stem Cells on Disc Cells from Degenerated Discs. Cells. Tissues. Organs 208, 76–88. doi:10.1159/000506350
Hirose, Y., Chiba, K., Karasugi, T., Nakajima, M., Kawaguchi, Y., Mikami, Y., 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, 1122–1129. doi:10.1016/j.ajhg.2008.03.013
Holm, S., Holm, A. K., Ekström, L., Karladani, A., and Hansson, T. (2004). Experimental Disc Degeneration Due to Endplate Injury. J. Spinal Disord. Tech. 17, 64–71. doi:10.1097/00024720-200402000-00012
Hoogendoorn, W. E., van Poppel, M. N. M., Bongers, P. M., Koes, B. W., and Bouter, L. M. (2000). Systematic Review of Psychosocial Factors at Work and Private Life as Risk Factors for Back Pain. Spine 25, 2114–2125. doi:10.1097/00007632-200008150-00017
Horner, H. A., and Urban, J. 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. Spine 26, 2543–2549. doi:10.1097/00007632-200112010-00006
Hoyland, J. A., Le Maitre, C., and Freemont, A. J. (2008). Investigation of the Role of IL-1 and TNF in Matrix Degradation in the Intervertebral Disc. Rheumatology 47, 809–814. doi:10.1093/rheumatology/ken056
Hu, J., Yan, Q., Shi, C., Tian, Y., Cao, P., and Yuan, W. (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, 79–89.
Huang, Y.-C., Hu, Y., Li, Z., and Luk, K. D. K. (2018). Biomaterials for Intervertebral Disc Regeneration: Current Status and Looming Challenges. J. Tissue. Eng. Regen. Med. 12, 2188–2202. doi:10.1002/term.2750
Hwang, M. H., Son, H.-G., Kim, J., and Choi, H. (2020). In Vitro model of Distinct Catabolic and Inflammatory Response Patterns of Endothelial Cells to Intervertebral Disc Cell Degeneration. Sci. Rep. 10, 20596. doi:10.1038/s41598-020-77785-6
Iatridis, J. C., Mente, P. L., Stokes, I. A. F., Aronsson, D. D., and Alini, M. (1999). Compression-induced Changes in Intervertebral Disc Properties in a Rat Tail Model. Spine 24, 996–1002. doi:10.1097/00007632-199905150-00013
Illien-Jünger, S., Walter, B. A., Mayer, J. E., Hecht, A. C., and Iatridis, J. C. (2014). “Intervertebral Disc Culture Models and Their Applications to Study Pathogenesis and Repair,” in The Intervertebral Disc (Vienna: Springer), 353–371. doi:10.1007/978-3-7091-1535-0_22
Isa, I. L. M., Mokhtar, S. A., Abbah, S. A., Fauzi, M. B., Devitt, A., and Pandit, A. (2022). Intervertebral Disc Degeneration: Biomaterials and Tissue Engineering Strategies toward Precision Medicine. Adv. Healthc. Mater., e2102530. doi:10.1002/adhm.202102530
Jauregui, J. J., Tran, A., Kaveeshwar, S., Nadarajah, V., Chaudhri, M. W., Henn, R. F., et al. (2021). Diagnosing a Periprosthetic Shoulder Infection: A Systematic Review. J. Orthop. 26, 58–66. doi:10.1016/j.jor.2021.07.012
Jiang, J., Cole, J. B., Freebern, E., Da, Y., VanRaden, P. M., and Ma, L. (2019). Functional Annotation and Bayesian Fine-Mapping Reveals Candidate Genes for Important Agronomic Traits in Holstein Bulls. Commun. Biol. 2, 1–12. doi:10.1038/s42003-019-0454-y
Johnson, W. E. B., Caterson, B., Eisenstein, S. M., Hynds, D. L., Snow, D. M., and Roberts, S. (2002). Human Intervertebral Disc Aggrecan Inhibits Nerve Growth In Vitro. Arthritis & Rheumatism 46, 2658–2664. doi:10.1002/art.10585
Johnson, W. E. B., Caterson, B., Eisenstein, S. M., and Roberts, S. (2005). Human Intervertebral Disc Aggrecan Inhibits Endothelial Cell Adhesion and Cell Migration In Vitro. Spine 30, 1139–1147. doi:10.1097/01.brs.0000162624.95262.73
Johnson, W. E. B., Stephan, S., and Roberts, S. (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. doi:10.1186/ar2405
Johnson, Z. I., Doolittle, A. C., Snuggs, J. W., Shapiro, I. M., Le Maitre, C. L., and Risbud, M. 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, 17561–17575. doi:10.1074/jbc.M117.790378
Johnson, Z., Schoepflin, Z. R., Schoepflin, Z., Choi, H., Shapiro, I., and Risbud, M. (2015). Disc in Flames: Roles of TNF-α and IL-1β in Intervertebral Disc Degeneration. eCM 30, 104–117. doi:10.22203/ecm.v030a08
Jones, A. C., and Wilcox, R. K. (2008). Finite Element Analysis of the Spine: towards a Framework of Verification, Validation and Sensitivity Analysis. Med. Eng. Phys. 30, 1287–1304. doi:10.1016/j.medengphy.2008.09.006
Jones, P., Gardner, L., Menage, J., Williams, G. T., and Roberts, S. (2008). Intervertebral Disc Cells as Competent Phagocytes In Vitro: Implications for Cell Death in Disc Degeneration. Arthritis Res. Ther. 10, R86. doi:10.1186/ar2466
Kang, J. D., Stefanovic-Racic, M., McIntyre, L. A., Georgescu, H. I., and Evans, C. H. (1997). Toward a Biochemical Understanding of Human Intervertebral Disc Degeneration and Herniation. Spine 22, 1065–1073. doi:10.1097/00007632-199705150-00003
Karasugi, T., Semba, K., Hirose, Y., Kelempisioti, A., Nakajima, M., Miyake, A., et al. (2009). Association of the Tag SNPs in the HumanSKTGene (KIAA1217) with Lumbar Disc Herniation. J. Bone Mineral Res. 24, 1537–1543. doi:10.1359/jbmr.090314
Karppinen, J., Daavittila, I., Solovieva, S., Kuisma, M., Taimela, S., Natri, A., et al. (2008). Genetic Factors Are Associated with Modic Changes in Endplates of Lumbar Vertebral Bodies. Spine 33, 1236–1241. doi:10.1097/BRS.0b013e318170fd0e
Kashtan, N., Itzkovitz, S., Milo, R., and Alon, U. (2004). Efficient Sampling Algorithm for Estimating Subgraph Concentrations and Detecting Network Motifs. Bioinformatics 20, 1746–1758. doi:10.1093/bioinformatics/bth163
Kawaguchi, Y., Osada, R., Kanamori, M., Ishihara, H., Ohmori, K., Matsui, H., et al. (1999). Association between an Aggrecan Gene Polymorphism and Lumbar Disc Degeneration. Spine 24, 2456–2460. doi:10.1097/00007632-199912010-00006
Kelempisioti, A., Eskola, P. J., Okuloff, A., Karjalainen, U., Takatalo, J., Daavittila, I., et al. (2011). Genetic Susceptibility of Intervertebral Disc Degeneration Among Young Finnish Adults. Bmc. Med. Genet. 12, 12–153. doi:10.1186/1471-2350-12-153
Kelly, T.-A. N., Ng, K. W., Ateshian, G. A., and Hung, C. T. (2009). Analysis of Radial Variations in Material Properties and Matrix Composition of Chondrocyte-Seeded Agarose Hydrogel Constructs. Osteoarthr. Cartil. 17, 73–82. doi:10.1016/j.joca.2008.05.019
Kepler, C. K., Ponnappan, R. K., Tannoury, C. A., Risbud, M. V., and Anderson, D. G. (2013). The Molecular Basis of Intervertebral Disc Degeneration. Spine J. 13, 318–330. doi:10.1016/j.spinee.2012.12.003
Kim, A.-G., Kim, T.-W., Kwon, W.-K., Lee, K.-H., Jeong, S., Hwang, M.-H., et al. (2021). Microfluidic Chip with Low Constant-Current Stimulation (LCCS) Platform: Human Nucleus Pulposus Degeneration In Vitro Model for Symptomatic Intervertebral Disc. Micromachines 12, 1291. doi:10.3390/mi12111291
Kim, N. K., Shin, D. A., Han, I. B., Yoo, E. H., Kim, S. H., and Chung, S. S. (2010). The Association of Aggrecan Gene Polymorphism with the Risk of Intervertebral Disc Degeneration. Acta Neurochir. 153, 129–133. doi:10.1007/s00701-010-0831-2
Krock, E., Rosenzweig, D. H., Chabot-Doré, A. J., Jarzem, P., Weber, M. H., Ouellet, J. A., et al. (2014). Painful, Degenerating Intervertebral Discs Up-Regulate Neurite Sprouting andCGRPthrough Nociceptive Factors. J. Cell. Mol. Med. 18, 1213–1225. doi:10.1111/jcmm.12268
Lakstins, K., Arnold, L., Gunsch, G., Flanigan, D., Khan, S., Gadde, N., et al. (2021). Characterization of the Human Intervertebral Disc Cartilage Endplate at the Molecular, Cell, and Tissue Levels. J. Orthop. Res. 39, 1898–1907. doi:10.1002/jor.24854
Lama, P., Le Maitre, C. L., Harding, I. J., Dolan, P., and Adams, M. A. (2018). Nerves and Blood Vessels in Degenerated Intervertebral Discs Are Confined to Physically Disrupted Tissue. J. Anat. 233, 86–97. doi:10.1111/joa.12817
Larson, J. W., Levicoff, E. A., Gilbertson, L. G., and Kang, J. D. (2006). Biologic Modification of Animal Models of Intervertebral Disc Degeneration. J. Bone. Jt. Surg. Am. 88 (Suppl. 2), 83–87. doi:10.2106/JBJS.F.00043
Lazarus, E., Bermudez-Lekerika, P., Farchione, D., Schofield, T., Howard, S., Mambetkadyrov, I., et al. (2021). Sulfated Hydrogels in Intervertebral Disc and Cartilage Research. Cells 10, 3568. doi:10.3390/cells10123568
Le Maitre, C., Hoyland, J., and Freemont, A. J. (2007). Catabolic Cytokine Expression in Degenerate and Herniated Human Intervertebral Discs: IL-1β and TNFα Expression Profile. Arthritis Res. Ther. 9, R77. doi:10.1186/ar2275
Le Maitre, C. L., Frain, J., Fotheringham, A. P., Freemont, A. J., and Hoyland, J. 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. Biorheology 45, 563–575. doi:10.3233/BIR-2008-0498
Le Maitre, C. L., Frain, J., Millward-Sadler, J., Fotheringham, A. P., Freemont, A. J., and Hoyland, J. A. (2009a). Altered Integrin Mechanotransduction in Human Nucleus Pulposus Cells Derived from Degenerated Discs. Arthritis. Rheum. 60(2), 460–469. doi:10.1002/art.24248
Le Maitre, C. L., Freemont, A. J., and Hoyland, J. 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. doi:10.1186/ar2808
Lee, J. M., Song, J. Y., Baek, M., Jung, H.-Y., Kang, H., Han, I. B., et al. (2011). Interleukin-1β Induces Angiogenesis and Innervation in Human Intervertebral Disc Degeneration. J. Orthop. Res. 29, 265–269. doi:10.1002/jor.21210
Li, C.-G., Liang, Q.-Q., Zhou, Q., Menga, E., Cui, X.-J., Shu, B., et al. (2009). A Continuous Observation of the Degenerative Process in the Intervertebral Disc of Smad3 Gene Knock-Out Mice. Spine 34, 1363–1369. doi:10.1097/BRS.0b013e3181a3c7c7
Li, H., Li, W., Zhang, L., He, J., Tang, L., Li, Z., et al. (2022). Comprehensive Network Analysis Identified SIRT7, NTRK2, and CHI3L1 as New Potential Markers for Intervertebral Disc Degeneration. J. Oncol. 2022, 1–17. doi:10.1155/2022/4407541
Li, W., Lai, K., Chopra, N., Zheng, Z., Das, A., and Diwan, A. D. (2022). Gut-disc axis: A Cause of Intervertebral Disc Degeneration and Low Back Pain? Eur. Spine J. 31(4), 917–992. doi:10.1007/s00586-022-07152-8
Li, X., Lin, F., Wu, Y., Liu, N., Wang, J., Chen, R., et al. (2019). Resveratrol Attenuates Inflammation Environment-Induced Nucleus Pulposus Cell Senescence In Vitro. Biosci. Rep. 39, 30962260. doi:10.1042/BSR20190126
Li, Y., Li, K., Mao, L., Han, X., Zhang, K., Zhao, C., et al. (2016). Cordycepin Inhibits LPS-Induced Inflammatory and Matrix Degradation in the Intervertebral Disc. PeerJ 4, e1992. doi:10.7717/peerj.1992
Li, Z., Shao, Z., Chen, S., Huang, D., Peng, Y., Chen, S., et al. (2020). TIGAR Impedes Compression-Induced Intervertebral Disc Degeneration by Suppressing Nucleus Pulposus Cell Apoptosis and Autophagy. J. Cell. Physiology 235, 1780–1794. doi:10.1002/jcp.29097
Lin, X., and Lin, Q. (2020). MiRNA-495-3p Attenuates TNF-α Induced Apoptosis and Inflammation in Human Nucleus Pulposus Cells by Targeting IL5RA. Inflammation 43, 1797–1805. doi:10.1007/s10753-020-01254-5
Liu, G.-Z., Ishihara, H., Osada, R., Kimura, T., and Tsuji, H. (2001). Nitric Oxide Mediates the Change of Proteoglycan Synthesis in the Human Lumbar Intervertebral Disc in Response to Hydrostatic Pressure. Spine 26, 134–141. doi:10.1097/00007632-200101150-00005
Liu, L., He, J., Liu, C., Yang, M., Fu, J., Yi, J., 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, 475–484. doi:10.3892/ijmm.2020.4832
Liu, S., Wu, N., Liu, J., Liu, H., Su, X., Liu, Z., et al. (2016). Association betweenADAMTS-4gene Polymorphism and Lumbar Disc Degeneration in Chinese Han Population. J. Orthop. Res. 34, 860–864. doi:10.1002/jor.23081
Loibl, M., Wuertz-Kozak, K., Vadala, G., Lang, S., Fairbank, J., and Urban, J. P. (2019). Controversies in Regenerative Medicine: Should Intervertebral Disc Degeneration Be Treated with Mesenchymal Stem Cells? JOR Spine 2, e1043. doi:10.1002/jsp2.1043
Lü, D.-S., Shono, Y., Oda, I., Abumi, K., and Kaneda, K. (1997). Effects of Chondroitinase ABC and Chymopapain on Spinal Motion Segment Biomechanics. Spine 22, 1828–1834. ; discussion 1834-5. doi:10.1097/00007632-199708150-00006
Lyons, G., Eisenstein, S. M., and Sweet, M. B. E. (1981). Biochemical Changes in Intervertebral Disc Degeneration. Biochimica Biophysica Acta (BBA) - General Subj. 673, 443–453. doi:10.1016/0304-4165(81)90476-1
MacLean, M. A., Touchette, C. J., Han, J. H., Christie, S. D., and Pickett, G. E. (2020). Gender Differences in the Surgical Management of Lumbar Degenerative Disease: a Scoping Review. Spine 32, 799–816. doi:10.3171/2019.11.SPINE19896
Maidhof, R., Jacobsen, T., Papatheodorou, A., and Chahine, N. O. (2014). Inflammation Induces Irreversible Biophysical Changes in Isolated Nucleus Pulposus Cells. PLoS ONE 9, e99621. doi:10.1371/journal.pone.0099621
Mainardi, A., Cambria, E., Occhetta, P., Martin, I., Barbero, A., Schären, S., 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. doi:10.3389/fbioe.2021.826867
Malandrino, A., Noailly, J., and Lacroix, D. (2011). The Effect of Sustained Compression on Oxygen Metabolic Transport in the Intervertebral Disc Decreases with Degenerative Changes. PLoS. Comput. Biol. 7, e1002112. doi:10.1371/journal.pcbi.1002112
Manolio, T. A., Collins, F. S., Cox, N. J., Goldstein, D. B., Hindorff, L. A., Hunter, D. J., et al. (2009). Finding the Missing Heritability of Complex Diseases. Nature 461, 747–753. doi:10.1038/nature08494
Mayer, J. E., Iatridis, J. C., Chan, D., Qureshi, S. A., Gottesman, O., and Hecht, A. C. (2013). Genetic Polymorphisms Associated with Intervertebral Disc Degeneration. Spine J. 13, 299–317. doi:10.1016/j.spinee.2013.01.041
McAfee, P. C., Cunningham, B., Holsapple, G., Adams, K., Blumenthal, S., Guyer, R. 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. Spine 30, 1576–1583. doi:10.1097/01.brs.0000170561.25636.1c
Merceron, C., Mangiavini, L., Robling, A., Wilson, T. L., Giaccia, A. J., Shapiro, I. M., et al. (2014). Loss of HIF-1α in the Notochord Results in Cell Death and Complete Disappearance of the Nucleus Pulposus. PLoS. ONE 9, e110768. doi:10.1371/journal.pone.0110768
Millecamps, M., Tajerian, M., Naso, L., Sage, H. E., and Stone, L. S. (2012). Lumbar Intervertebral Disc Degeneration Associated with Axial and Radiating Low Back Pain in Ageing SPARC-Null Mice. Pain 153, 1167–1179. doi:10.1016/j.pain.2012.01.027
Millward-Sadler, S. J., Costello, P. W., Freemont, A. J., and Hoyland, J. 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. doi:10.1186/ar2693
Min, S.-K., Nakazato, K., Yamamoto, Y., Gushiken, K., Fujimoto, H., Fujishiro, H., 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, 2552–2557. doi:10.1177/0363546510376714
Mio, F., Chiba, K., Hirose, Y., Kawaguchi, Y., Mikami, Y., Oya, T., 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, 1271–1277. doi:10.1086/522377
Miyazaki, T., Kobayashi, S., Takeno, K., Meir, A., Urban, J., and Baba, H. (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 A 15, 3835–3846. doi:10.1089/ten.tea.2009.0250
Moon, S. M., Yoder, J. H., Wright, A. C., Smith, L. J., Vresilovic, E. J., and Elliott, D. M. (2013). Evaluation of Intervertebral Disc Cartilaginous Endplate Structure Using Magnetic Resonance Imaging. Eur. Spine J. 22, 1820–1828. doi:10.1007/s00586-013-2798-1
Mountjoy, E., Schmidt, E. M., Carmona, M., Schwartzentruber, J., Peat, G., Miranda, A., et al. (2021). An Open Approach to Systematically Prioritize Causal Variants and Genes at All Published Human GWAS Trait-Associated Loci. Nat. Genet. 53, 1527–1533. doi:10.1038/s41588-021-00945-5
Mu, J., Ge, W., Zuo, X., Chen, Y., and Huang, C. (2013). Analysis of Association between IL-1β, CASP-9, and GDF5 Variants and Low-Back Pain in Chinese Male Soldiers. Spi 19, 243–247. doi:10.3171/2013.4.SPINE12782
Munir, S., Freidin, M. B., Rade, M., Määttä, J., Livshits, G., and Williams, F. M. K. (2018). Endplate Defect Is Heritable, Associated with Low Back Pain and Triggers Intervertebral Disc Degeneration. Spine 43, 1496–1501. doi:10.1097/BRS.0000000000002721
Näkki, A., Videman, T., Kujala, U. M., Suhonen, M., Männikkö, M., Peltonen, L., 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, 747–752. doi:10.3899/jrheum.100080
Naqvi, S. M., and Buckley, C. 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 A 21, 288–299. doi:10.1089/ten.TEA.2013.0719
Neidlinger-Wilke, C., Mietsch, A., Rinkler, C., Wilke, H.-J., Ignatius, A., and Urban, J. (2012). Interactions of Environmental Conditions and Mechanical Loads Have Influence on Matrix Turnover by Nucleus Pulposus Cells. J. Orthop. Res. 30, 112–121. doi:10.1002/jor.21481
Newell, N., Little, J., Christou, A., Adams, M., Adam, C., and Masouros, S. (2017). Biomechanics of the Human Intervertebral Disc: A Review of Testing Techniques and Results. J. Mech. Behav. Biomed. Mater. 69, 420–434. doi:10.1016/j.jmbbm.2017.01.037
Nicolae, D. L., Gamazon, E., Zhang, W., Duan, S., Dolan, M. E., and Cox, N. J. (2010). Trait-associated SNPs Are More Likely to Be eQTLs: Annotation to Enhance Discovery from GWAS. PLoS Genet. 6, e1000888. doi:10.1371/journal.pgen.1000888
Niemeyer, F., Galbusera, F., Tao, Y., Kienle, A., Beer, M., and Wilke, H.-J. (2021). A Deep Learning Model for the Accurate and Reliable Classification of Disc Degeneration Based on MRI Data. Invest. Radiol. 56, 78–85. doi:10.1097/RLI.0000000000000709
Noponen-Hietala, N., Kyllönen, E., Männikkö, M., Ilkko, E., Karppinen, J., Ott, J., et al. (2003). Sequence Variations in the Collagen IX and XI Genes Are Associated with Degenerative Lumbar Spinal Stenosis. Ann. Rheumatic Dis. 62, 1208–1214. doi:10.1136/ard.2003.008334
Noponen-Hietala, N., Virtanen, I., Karttunen, R., Schwenke, S., Jakkula, E., Li, H., et al. (2005). Genetic Variations in IL6 Associate with Intervertebral Disc Disease Characterized by Sciatica. Pain 114, 186–194. doi:10.1016/j.pain.2004.12.015
Novais, E. J., Choi, H., Madhu, V., Suyama, K., Anjo, S. I., Manadas, B., et al. (2021). Hypoxia and Hypoxia-Inducible Factor-1α Regulate Endoplasmic Reticulum Stress in Nucleus Pulposus Cells. Am. J. Pathology 191, 487–502. doi:10.1016/j.ajpath.2020.11.012
Oliver, J. C., Bland, L. A., Oettinger, C. W., Arduino, M. J., McAllister, S. K., Aguero, S. M., et al. (1993). Cytokine Kinetics in an In Vitro Whole Blood Model Following an Endotoxin Challenge. Lymphokine Cytokine Res. 12, 115–120.
Omair, A., Holden, M., Lie, B. A., Reikeras, O., and Brox, J. 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, 1–9. doi:10.1186/1471-2474-14-105
Oswald, K. A. C., Bigdon, S. F., Croft, A. S., Bermudez-Lekerika, P., Bergadano, A., Gantenbein, B., et al. (2021). Establishment of a Novel Method for Spinal Discectomy Surgery in Elderly Rats in an In Vivo Spinal Fusion Model. MPs 4, 79. doi:10.3390/mps4040079
Öztürk, E., ArlovAksel, Ø., Aksel, S., Li, L., Ornitz, D. M., Skjåk-Braek, G., et al. (2016). Sulfated Hydrogel Matrices Direct Mitogenicity and Maintenance of Chondrocyte Phenotype through Activation of FGF Signaling. Adv. Funct. Mat. 26, 3649–3662. doi:10.1002/adfm.201600092
Paassilta, P., Lohiniva, J., Göring, H. H., Perälä, M., Räinä, S. S., Karppinen, J., et al. (2001). Identification of a Novel Common Genetic Risk Factor for Lumbar Disk Disease. JAMA 285, 1843–1849. doi:10.1001/jama.285.14.1843
Parenteau, C. S., Lau, E. C., Campbell, I. C., and Courtney, A. (2021). Prevalence of Spine Degeneration Diagnosis by Type, Age, Gender, and Obesity Using Medicare Data. Sci. Rep. 11, 5389. doi:10.1038/s41598-021-84724-6
Pattappa, G., Peroglio, M., Peroglio, M., Sakai, D., Mochida, J., Benneker, L., et al. (2014). CCL5/RANTES Is a Key Chemoattractant Released by Degenerative Intervertebral Discs in Organ Culture. eCM 27, 124–136. ; discussion 136. doi:10.22203/ecm.v027a10
Paul, C. P. L., Schoorl, T., Zuiderbaan, H. A., Zandieh Doulabi, B., van der Veen, A. J., van de Ven, P. M., et al. (2013). Dynamic and Static Overloading Induce Early Degenerative Processes in Caprine Lumbar Intervertebral Discs. PLoS ONE 8, e62411. doi:10.1371/journal.pone.0062411
Paul, C. P. L., Zuiderbaan, H. A., Zandieh Doulabi, B., van der Veen, A. J., van de Ven, P. M., Smit, T. H., et al. (2012). Simulated-physiological Loading Conditions Preserve Biological and Mechanical Properties of Caprine Lumbar Intervertebral Discs in Ex Vivo Culture. PLoS ONE 7, e33147. doi:10.1371/journal.pone.0033147
Pearce, D., Fischer, S., Huda, F., and Vahdati, A. (2020). Applications of Computer Modeling and Simulation in Cartilage Tissue Engineering. Tissue. Eng. Regen. Med. 17, 1–13. doi:10.1007/s13770-019-00216-9
Pfannkuche, J.-J., Guo, W., Cui, S., Ma, J., Lang, G., Peroglio, M., 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, 304–321. doi:10.1080/03008207.2019.1665652
Phillips, F. M., Reuben, J., and Wetzel, F. T. (2002). Intervertebral Disc Degeneration Adjacent to a Lumbar Fusion. An Experimental Rabbit Model. J. Bone Jt. Surg. Br. volume 84, 289–294. doi:10.1302/0301-620x.84b2.11937
Phillips, K. L. E., Chiverton, N., Michael, A. L., Cole, A. A., Breakwell, L. M., Haddock, G., 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. doi:10.1186/ar4408
Phillips, K. L. E., Cullen, K., Chiverton, N., Michael, A. L. R., Cole, A. A., Breakwell, L. 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, 1165–1177. doi:10.1016/j.joca.2015.02.017
Phillips, K. L. E., Jordan-Mahy, N., Nicklin, M. J. H., and Le Maitre, C. L. (2013b). Interleukin-1 Receptor Antagonist Deficient Mice Provide Insights into Pathogenesis of Human Intervertebral Disc Degeneration. Ann. Rheum. Dis. 72(11), 1860–1867. doi:10.1136/annrheumdis-2012-202266
Pincus, T., Burton, A. K., Vogel, S., and Field, A. P. (2002). A Systematic Review of Psychological Factors as Predictors of Chronicity/disability in Prospective Cohorts of Low Back Pain. Spine 27, E109–E120. doi:10.1097/00007632-200203010-00017
Pluijm, S. M. F., Hw, V. E., Bravenboer, N., Uitterlinden, A. G., Smit, J. H., Pols, H. A., et al. (2004). Collagen Type I 1 Sp1 Polymorphism, Osteoporosis, and Intervertebral Disc Degeneration in Older Men and Women. Ann. Rheumatic Dis. 63, 71–77. doi:10.1136/ard.2002.002287
Pockert, A. J., Richardson, S. M., Le Maitre, C. L., Lyon, M., Deakin, J. A., Buttle, D. J., et al. (2009). Modified Expression of the ADAMTS Enzymes and Tissue Inhibitor of Metalloproteinases 3 during Human Intervertebral Disc Degeneration. Arthritis Rheum. 60, 482–491. doi:10.1002/art.24291
Purmessur, D., Freemont, A. J., and Hoyland, J. A. (2008). Expression and Regulation of Neurotrophins in the Nondegenerate and Degenerate Human Intervertebral Disc. Arthritis Res. Ther. 10, R99. doi:10.1186/ar2487
Raine, E. V. A., Reynard, L. N., van de Laar, I. M. B. H., Bertoli-Avella, A. M., and Loughlin, J. (2014). Identification and Analysis of a SMAD3 Cis-Acting eQTL Operating in Primary Osteoarthritis and in the Aneurysms and Osteoarthritis Syndrome. Osteoarthr. Cartil. 22, 698–705. doi:10.1016/j.joca.2014.02.931
Rajasekaran, S., Kanna, R. M., Senthil, N., Raveendran, M., Ranjani, V., Cheung, K. M. C., et al. (2014). Genetic Susceptibility of Lumbar Degenerative Disc Disease in Young Indian Adults. Eur. Spine J. 24, 1969–1975. doi:10.1007/s00586-014-3687-y
Rajasekaran, S., Soundararajan, D. C. R., Tangavel, C., Muthurajan, R., Sri Vijay Anand, K. S., Matchado, M. S., et al. (2020). Human Intervertebral Discs Harbour a Unique Microbiome and Dysbiosis Determines Health and Disease. Eur. Spine J. 29, 1621–1640. doi:10.1007/s00586-020-06446-z
Rastogi, A., Thakore, P., Leung, A., Benavides, M., Machado, M., Morschauser, M. A., et al. (2009). Environmental Regulation of Notochordal Gene Expression in Nucleus Pulposus Cells. J. Cell. Physiol. 220, 698–705. doi:10.1002/jcp.21816
Rätsep, T., Minajeva, A., and Asser, T. (2013). Relationship between Neovascularization and Degenerative Changes in Herniated Lumbar Intervertebral Discs. Eur. Spine J. 22, 2474–2480. doi:10.1007/s00586-013-2842-1
Rijsbergen, M. V., van Rietbergen, B., Barthelemy, V., Eltes, P., Lazáry, Á., Lacroix, D., 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 ONE 13, e0200899. doi:10.1371/journal.pone.0200899
Rikard, S. M., Athey, T. L., Nelson, A. R., Christiansen, S. L. M., Lee, J.-J., Holmes, J. 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. doi:10.3389/fphys.2019.01481
Roberts, S., Menage, J., and Urban, J. P. G. (1989). Biochemical and Structural Properties of the Cartilage End-Plate and its Relation to the Intervertebral Disc. Spine 14, 166–174. doi:10.1097/00007632-198902000-00005
Robinson, M. R., English, G., English, G., Moser, G., Lloyd-Jones, L. R., Triplett, M. A., et al. (2017). Genotype-covariate Interaction Effects and the Heritability of Adult Body Mass Index. Nat. Genet. 49, 1174–1181. doi:10.1038/ng.3912
Rosiak, P., Latanska, I., Paul, P., Sujka, W., and Kolesinska, B. (2021). Modification of Alginates to Modulate Their Physic-Chemical Properties and Obtain Biomaterials with Different Functional Properties. Molecules 26, 7264. doi:10.3390/molecules26237264
Ruiz Wills, C., Foata, B., González Ballester, M. Á., Karppinen, J., and Noailly, J. (2018). Theoretical Explorations Generate New Hypotheses about the Role of the Cartilage Endplate in Early Intervertebral Disk Degeneration. Front. Physiol. 9, 1210. doi:10.3389/fphys.2018.01210
Sadowska, A., Altinay, B., Hitzl, W., Ferguson, S. J., and Wuertz-Kozak, K. (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. doi:10.3389/fphar.2020.00952
Sadowska, A., Kameda, T., Kameda, T., Krupkova, O., and Wuertz-Kozak, K. (2018). Osmosensing, Osmosignalling and Inflammation: How Intervertebral Disc Cells Respond to Altered Osmolarity. eCM 36, 231–250. doi:10.22203/eCM.v036a17
Saggese, T., Thambyah, A., Wade, K., and McGlashan, S. R. (2020). Differential Response of Bovine Mature Nucleus Pulposus and Notochordal Cells to Hydrostatic Pressure and Glucose Restriction. Cartilage 11, 221–233. doi:10.1177/1947603518775795
Samartzis, D., Karppinen, J., Chan, D., Luk, K. D. K., and Cheung, K. 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 & Rheumatism 64, 1488–1496. doi:10.1002/art.33462
Sambrook, P. N., MacGregor, A. J., and Spector, T. D. (1999). Genetic Influences on Cervical and Lumbar Disc Degeneration: A Magnetic Resonance Imaging Study in Twins. Arthritis & Rheumatism 42, 366–372. doi:10.1002/1529-0131(199902)42:2<366::aid-anr20>3.0.co;2-6
Saravi, B., Lang, G., Grad, S., Alini, M., Richards, R. G., Schmal, H., et al. (2021). A Proinflammatory, Degenerative Organ Culture Model to Simulate Early-Stage Intervertebral Disc Disease. JoVE 168. doi:10.3791/62100
Seki, S., Kawaguchi, Y., Chiba, K., Mikami, Y., Kizawa, H., Oya, T., et al. (2005). A Functional SNP in CILP, Encoding Cartilage Intermediate Layer Protein, Is Associated with Susceptibility to Lumbar Disc Disease. Nat. Genet. 37, 607–612. doi:10.1038/ng1557
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 13, 2498–2504. doi:10.1101/gr.1239303
Showalter, B. L., Beckstein, J. C., Martin, J. T., Beattie, E. E., Orías, A. A. E., Schaer, T. P., et al. (2012). Comparison of Animal Discs Used in Disc Research to Human Lumbar Disc. Spine 37, E900–E907. doi:10.1097/BRS.0b013e31824d911c
Smith, L., Chiaro, J. A., Chiaro, J., Nerurkar, N., Cortes, D., Horava, S., et al. (2011). Nucleus Pulposus Cells Synthesize a Functional Extracellular Matrix and Respond to Inflammatory Cytokine Challenge Following Long-Term Agarose Culture. eCM 22, 291–301. doi:10.22203/ecm.v022a22
Snuggs, J., Bunning, R. A., Bunning, R., and Le Maitre, C. (2021). Osmotic Adaptation of Nucleus Pulposus Cells: the Role of Aquaporin 1, Aquaporin 4 and Transient Receptor Potential Vanilloid 4. eCM 41, 121–141. doi:10.22203/eCM.v041a09
Solovieva, S., Kouhia, S., Leino-Arjas, P., Ala-Kokko, L., Luoma, K., Raininko, R., et al. (2004). Interleukin 1 Polymorphisms and Intervertebral Disc Degeneration. Epidemiology 15, 626–633. doi:10.1097/01.ede.0000135179.04563.35
Solovieva, S., Lohiniva, J., Leino-Arjas, P., Raininko, R., Luoma, K., Ala-Kokko, L., et al. (2006). Intervertebral Disc Degeneration in Relation to the COL9A3 and the IL-1ß Gene Polymorphisms. Eur. Spine J. 15, 613–619. doi:10.1007/s00586-005-0988-1
Solovieva, S., Noponen, N., Männikkö, M., Leino-Arjas, P., Luoma, K., Raininko, R., et al. (2007). Association between the Aggrecan Gene Variable Number of Tandem Repeats Polymorphism and Intervertebral Disc Degeneration. Spine 32, 1700–1705. doi:10.1097/BRS.0b013e3180b9ed51
Song, Y.-Q., Cheung, K. M. C., Ho, D. W. H., Poon, S. C. S., Chiba, K., Kawaguchi, Y., et al. (2008). Association of the Asporin D14 Allele with Lumbar-Disc Degeneration in Asians. Am. J. Hum. Genet. 82, 744–747. doi:10.1016/j.ajhg.2007.12.017
Song, Y.-Q., Karasugi, T., Cheung, K. M. C., Chiba, K., Ho, D. W. H., Miyake, A., et al. (2013). Lumbar Disc Degeneration Is Linked to a Carbohydrate Sulfotransferase 3 Variant. J. Clin. Invest. 123, 4909–4917. doi:10.1172/JCI69277
Sun, Q., Crowley, C. A., Huang, L., Wen, J., Chen, J., Bao, E. 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. doi:10.1016/j.xhgg.2021.100063
Sun, Y., Lv, M., Zhou, L., Tam, V., Lv, F., Chan, D., et al. (2015). Enrichment of Committed Human Nucleus Pulposus Cells Expressing Chondroitin Sulfate Proteoglycans under Alginate Encapsulation. Osteoarthr. Cartil. 23, 1194–1203. doi:10.1016/j.joca.2015.02.166
Sun, Z.-M., Ling, M., Huo, Y., Chang, Y., Li, Y., Qin, H., et al. (2011). Caspase 9 Gene Polymorphism and Susceptibility to Lumbar Disc Disease in the Han Population in Northern China. Connect. Tissue Res. 52, 198–202. doi:10.3109/03008207.2010.510914
Sun, Z.-m., Miao, L., Zhang, Y.-g., and Ming, L. (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, 181–185. doi:10.1080/03008200802585630
Sun, Z., Liu, B., and Luo, Z.-J. (2020). The Immune Privilege of the Intervertebral Disc: Implications for Intervertebral Disc Degeneration Treatment. Int. J. Med. Sci. 17, 685–692. doi:10.7150/ijms.42238
Suyama, K., Sakai, D., Hirayama, N., Nakamura, Y., Matsushita, E., Terayama, H., 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, 5539–5551. doi:10.1111/jcmm.13828
Swärd, L., Hellström, M., Jacobsson, B. O., Nyman, R., and Peterson, L. (1991). Disc Degeneration and Associated Abnormalities of the Spine in Elite Gymnasts. Spine 16, 437–443. doi:10.1097/00007632-199104000-00009
Takahashi, M., Haro, H., Wakabayashi, Y., Kawa-uchi, T., Komori, H., and Shinomiya, K. (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. volume 83, 491–495. doi:10.1302/0301-620x.83b4.11617
Tan, H., Zhao, J., Jiang, J., and Ren, Y. (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 Investigation 72, 576–579. doi:10.3109/00365513.2012.713176
Teraguchi, M., Yoshimura, N., Hashizume, H., Muraki, S., Yamada, H., Minamide, A., 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, 104–110. doi:10.1016/j.joca.2013.10.019
Thorpe, A. A., Bach, F. C., Tryfonidou, M. A., Le Maitre, C. L., Mwale, F., Diwan, A. D., et al. (2018). Leaping the Hurdles in Developing Regenerative Treatments for the Intervertebral Disc from Preclinical to Clinical. JOR. Spine 1, e1027. doi:10.1002/jsp2.1027
Tilkeridis, C., Bei, T., Garantziotis, S., and Stratakis, C. A. (2005). Association of a COL1A1 Polymorphism with Lumbar Disc Disease in Young Military Recruits. J. Med. Genet. 42, e44. doi:10.1136/jmg.2005.033225
Tolofari, S. K., Richardson, S. M., Freemont, A. J., and Hoyland, J. 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, R1–R11. doi:10.1186/ar2898
Tsai, T.-T., Danielson, K. G., Guttapalli, A., Oguz, E., Albert, T. J., Shapiro, I. M., et al. (2006). TonEBP/OREBP Is a Regulator of Nucleus Pulposus Cell Function and Survival in the Intervertebral Disc. J. Biol. Chem. 281, 25416–25424. doi:10.1074/jbc.M601969200
Turgut, M., Uslu, S., Uysal, A. l., Yurtseven, M. E., and Ustün, H. (2003). Changes in Vascularity of Cartilage Endplate of Degenerated Intervertebral Discs in Response to Melatonin Administration in Rats. Neurosurg. Rev. 26, 133–138. doi:10.1007/s10143-003-0259-8
Ulrich, J. A., Liebenberg, E. C., Thuillier, D. U., and Lotz, J. C. (2007). ISSLS Prize Winner: Repeated Disc Injury Causes Persistent Inflammation. Spine 32, 2812–2819. doi:10.1097/BRS.0b013e31815b9850
Valdes, A. M., Hassett, G., Hart, D. J., and Spector, T. D. (2005). Radiographic Progression of Lumbar Spine Disc Degeneration Is Influenced by Variation at Inflammatory Genes. Spine 30, 2445–2451. doi:10.1097/01.brs.0000184369.79744.a5
van Uden, S., Silva-Correia, J., Oliveira, J. M., and Reis, R. L. (2017). Current Strategies for Treatment of Intervertebral Disc Degeneration: Substitution and Regeneration Possibilities. Biomater. Res. 21, 22. doi:10.1186/s40824-017-0106-6
Vazirinejad, R., Ahmadi, Z., Kazemi Arababadi, M., Hassanshahi, G., and Kennedy, D. (2014). The Biological Functions, Structure and Sources of CXCL10 and its Outstanding Part in the Pathophysiology of Multiple Sclerosis. Neuroimmunomodulation 21, 322–330. doi:10.1159/000357780
Vergroesen, P.-P. A., Kingma, I., Emanuel, K. S., Hoogendoorn, R. J. W., Welting, T. J., van Royen, B. J., et al. (2015). Mechanics and Biology in Intervertebral Disc Degeneration: a Vicious Circle. Osteoarthr. Cartil. 23, 1057–1070. doi:10.1016/j.joca.2015.03.028
Vergroesen, P.-P. A., van der Veen, A. J., van Royen, B. J., Kingma, I., and Smit, T. H. (2014). Intradiscal Pressure Depends on Recent Loading and Correlates with Disc Height and Compressive Stiffness. Eur. Spine J. 23, 2359–2368. doi:10.1007/s00586-014-3450-4
Videman, T., Saarela, J., Kaprio, J., Näkki, A., Levälahti, E., Gill, K., et al. (2009). Associations of 25 Structural, Degradative, and Inflammatory Candidate Genes with Lumbar Disc Desiccation, Bulging, and Height Narrowing. Arthritis Rheum. 60, 470–481. doi:10.1002/art.24268
Virgilio, K. M., Martin, K. S., Peirce, S. M., and Blemker, S. S. (2018). Agent-based Model Illustrates the Role of the Microenvironment in Regeneration in Healthy and Mdx Skeletal Muscle. J. Appl. Physiology 125, 1424–1439. doi:10.1152/japplphysiol.00379.2018
Virtanen, I. M., Karppinen, J., Taimela, S., Ott, J., Barral, S., Kaikkonen, K., et al. (2007). Occupational and Genetic Risk Factors Associated with Intervertebral Disc Disease. Spine 32, 1129–1134. doi:10.1097/01.brs.0000261473.03274.5c
Walter, B. A., Illien-Jünger, S., Nasser, P. R., Hecht, A. C., and Iatridis, J. C. (2014). Development and Validation of a Bioreactor System for Dynamic Loading and Mechanical Characterization of Whole Human Intervertebral Discs in Organ Culture. J. Biomechanics 47, 2095–2101. doi:10.1016/j.jbiomech.2014.03.015
Walter, B. A., Korecki, C. L., Purmessur, D., Roughley, P. J., Michalek, A. J., and Iatridis, J. C. (2011). Complex Loading Affects Intervertebral Disc Mechanics and Biology. Osteoarthr. Cartil. 19, 1011–1018. doi:10.1016/j.joca.2011.04.005
Walter, B., Purmessur, D., Purmessur, D., Moon, A., Occhiogrosso, J., Laudier, D., et al. (2016). Reduced Tissue Osmolarity Increases TRPV4 Expression and Pro-inflammatory Cytokines in Intervertebral Disc Cells. eCM 32, 123–136. doi:10.22203/ecm.v032a08
Wang, H., Zhou, Y., Huang, B., Liu, L.-T., Liu, M.-H., Wang, J., et al. (2014). Utilization of Stem Cells in Alginate for Nucleus Pulposus Tissue Engineering. Tissue Eng. Part A 20, 908–920. doi:10.1089/ten.TEA.2012.0703
Wang, M., Tang, D., Shu, B., Wang, B., Jin, H., Hao, S., et al. (2012). Conditional Activation of β-catenin Signaling in Mice Leads to Severe Defects in Intervertebral Disc Tissue. Arthritis & Rheumatism 64, 2611–2623. doi:10.1002/art.34469
Wang, S.-L., Yu, Y.-L., Tang, C.-L., and Lv, F.-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, 1522–1526. doi:10.3892/etm.2013.1348
Wáng, Y. X. J., Wáng, J.-Q., and Káplár, Z. (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, 199–206. doi:10.21037/qims.2016.04.06
Wangler, S., Menzel, U., Li, Z., Ma, J., Hoppe, S., Benneker, L. M., et al. (2019). CD146/MCAM Distinguishes Stem Cell Subpopulations with Distinct Migration and Regenerative Potential in Degenerative Intervertebral Discs. Osteoarthr. Cartil. 27, 1094–1105. doi:10.1016/j.joca.2019.04.002
Weiler, C., Nerlich, A. G., Bachmeier, B. E., and Boos, N. (2005). Expression and Distribution of Tumor Necrosis Factor Alpha in Human Lumbar Intervertebral Discs: A Study in Surgical Specimen and Autopsy Controls. Spine 30, 44–53. doi:10.1097/01.brs.0000149186.63457.20
Wilke, H. J., Neef, P., Caimi, M., Hoogland, T., and Claes, L. E. (1999). New In Vivo Measurements of Pressures in the Intervertebral Disc in Daily Life. Spine 24, 755–762. doi:10.1097/00007632-199904150-00005
Williams, F. M. K., Bansal, A. T., van Meurs, J. B., Bell, J. T., Meulenbelt, I., Suri, P., et al. (2013). Novel Genetic Variants Associated with Lumbar Disc Degeneration in Northern Europeans: a Meta-Analysis of 4600 Subjects. Ann. Rheum. Dis. 72, 1141–1148. doi:10.1136/annrheumdis-2012-201551
Williams, F. M. K., Popham, M., Hart, D. J., de Schepper, E., Bierma-Zeinstra, S., Hofman, A., et al. (2011a). GDF5 Single-Nucleotide Polymorphism Rs143383 Is Associated with Lumbar Disc Degeneration in Northern European Women. Arthritis & Rheumatism 63, 708–712. doi:10.1002/art.30169
Williams, F. M. K., Popham, M., Sambrook, P. N., Jones, A. F., Spector, T. D., and MacGregor, A. J. (2011b). Progression of Lumbar Disc Degeneration over a Decade: a Heritability Study. Ann. Rheumatic Dis. 70, 1203–1207. doi:10.1136/ard.2010.146001
Wöltje, M., and Böbel, M. (2017). “Natural Biodegradable Medical Polymers,” in Science and Principles of Biodegradable and Bioresorbable Medical Polymers. Editor X. Zhang (Amsterdam, The Netherlands: Woodhead Publishing), 351–376. doi:10.1016/B978-0-08-100372-5.00012-X
Wu, N., Chen, J., Liu, H., Zhao, L., Liu, S., Liu, J., 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, 686–694. doi:10.1002/jor.22582
Wuertz, K., Urban, J. P. G., Klasen, J., Ignatius, A., Wilke, H.-J., Claes, L., et al. (2007). Influence of Extracellular Osmolarity and Mechanical Stimulation on Gene Expression of Intervertebral Disc Cells. J. Orthop. Res. 25, 1513–1522. doi:10.1002/jor.20436
Xu, C., Luo, S., Wei, L., Wu, H., Gu, W., Zhou, W., et al. (2021). Integrated Transcriptome and Proteome Analyses Identify Novel Regulatory Network of Nucleus Pulposus Cells in Intervertebral Disc Degeneration. Bmc. Med. Genomics. 14, 1–12. doi:10.1186/s12920-021-00889-z
Xu, K., Wang, X., Zhang, Q., Liang, A., Zhu, H., Huang, D., et al. (2016). Sp1 Downregulates Proinflammatory Cytokine-Induced Catabolic Gene Expression in Nucleus Pulposus Cells. Mol. Med. Rep. 14, 3961–3968. doi:10.3892/mmr.2016.5730
Yang, J.-E., Zhao, K.-H., Qu, Y., and Zou, Y.-C. (2022). Increased Serum CXCL10 Levels Are Associated with Clinical Severity and Radiographic Progression in Patients with Lumbar Disc Degeneration. Clin. Chim. Acta 525, 15–22. doi:10.1016/j.cca.2021.12.006
Yang, J., Ferreira, T., Ferreira, T., Morris, A. P., Medland, S. E., Madden, P. A. F., et al. (2012). Conditional and Joint Multiple-SNP Analysis of GWAS Summary Statistics Identifies Additional Variants Influencing Complex Traits. Nat. Genet. 44, 369–375. doi:10.1038/ng.2213
Yao, Y., Deng, Q., Song, W., Zhang, H., Li, Y., Yang, Y., 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, 249–262. doi:10.1016/j.stemcr.2016.07.003
Yao, Y., Song, W., Deng, Q., Zhang, H., Wang, J., Liu, H., 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, 3001–3009. doi:10.3892/mmr.2017.6907
Yoon, S.-H., Miyazaki, M., Hong, S.-W., Tow, B., Morishita, Y., Hu, M., et al. (2008). A Porcine Model of Intervertebral Disc Degeneration Induced by Annular Injury Characterized with Magnetic Resonance Imaging and Histopathological Findings. Spi 8, 450–457. doi:10.3171/SPI/2008/8/5/450
Yu, H., Zhu, X., Greenbaum, D., Karro, J., and Gerstein, M. (2004). TopNet: a Tool for Comparing Biological Sub-networks, Correlating Protein Properties with Topological Statistics. Nucleic Acids Res. 32, 328–337. doi:10.1093/nar/gkh164
Yu, W., Fu, J., Liu, Y., Wu, Y., and Jiang, D. (2018). Osteogenic Protein-1 Inhibits Nucleus Pulposus Cell Apoptosis through Regulating the NF-Κb/ROS Pathway in an Inflammation Environment. Biosci. Rep. 38, BSR20181530. doi:10.1042/BSR20181530
Zamanian, C., Bhandarkar, A. R., Monie, D. D., Moinuddin, F., Vile, R. G., Quiñones-Hinojosa, A., et al. (2022). Systems Neuroimmunology: a Review of Multiomics Methodologies to Characterize Neuroimmunological Interactions in Spinal and Cranial Diseases. Neurosurg. Focus. 52, E9. doi:10.3171/2021.11.FOCUS21571
Zeiter, S., der Werf, M. V., and Ito, K. (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, 310–320. doi:10.1002/term.167
Zhang, C., Gullbrand, S. E., Schaer, T. P., Lau, Y. K., Jiang, Z., Dodge, G. R., et al. (2020). Inflammatory Cytokine and Catabolic Enzyme Expression in a Goat Model of Intervertebral Disc Degeneration. J. Orthop. Res. 38, 2521–2531. doi:10.1002/jor.24639
Zhang, J., Wang, X., Liu, H., Li, Z., Chen, F., Wang, H., 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, 697–705. doi:10.1002/jor.24204
Zhang, S., Ding, L., Li, X., and Fan, H. (2019). Identification of Biomarkers Associated with the Recurrence of Osteosarcoma Using ceRNA Regulatory Network Analysis. Int. J. Mol. Med. 43, 1723–1733. doi:10.3892/ijmm.2019.4108
Zhang, Z., Li, F., Tian, H., Guan, K., Zhao, G., Shan, J., 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, 314–321. doi:10.3760/cma.j.issn.0366-6999.20130584
Zhou, X., Cheung, C.-L., Karasugi, T., Karppinen, J., Samartzis, D., Hsu, Y.-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. doi:10.3389/fgene.2018.00267
Zhou, Z., Cui, S., Du, J., Richards, R. G., Alini, M., Grad, S., et al. (2021). One Strike Loading Organ Culture Model to Investigate the Post-traumatic Disc Degenerative Condition. J. Orthop. Transl. 26, 141–150. doi:10.1016/j.jot.2020.08.003
Zhu, G.-b., Jiang, X.-r., Xia, C.-l., Sun, Y.-j., Zeng, Q.-s., Wu, X.-m., et al. (2011). Association of FAS and FAS Ligand Polymorphisms with the Susceptibility and Severity of Lumbar Disc Degeneration in Chinese Han Population. Biomarkers 16, 485–490. doi:10.3109/1354750X.2011.598563
Zuo, R., Wang, Y., Li, J., Wu, J., Wang, W., Li, B., et al. (2019). Rapamycin Induced Autophagy Inhibits Inflammation-Mediated Endplate Degeneration by Enhancing Nrf2/Keap1 Signaling of Cartilage Endplate Stem Cells. Stem Cells 37, 828–840. doi:10.1002/stem.2999
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
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
IVD Intervertebral disc
LBP Low Back Pain
LD Linkage disequilibrium
LDD Lumbar Disc Degeneration
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
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.
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
Copyright © 2022 Bermudez-Lekerika, Crump, Tseranidou, Nüesch, Kanelis, Alminnawi, Baumgartner, Muñoz-Moya, Compte, Gualdi, Alexopoulos, Geris, Wuertz-Kozak, Le Maitre, Noailly and Gantenbein. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Benjamin Gantenbein, Benjamin.Gantenbein@dbmr.unibe.ch
†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