Single-cell RNA-sequencing atlas of bovine caudal intervertebral discs: Discovery of heterogeneous cell populations with distinct roles in homeostasis

Osteoporosis treatments for intervertebral disc degeneration and back pain: a perspective

Low back pain derived from intervertebral disc (IVD) degeneration is a debilitating spinal condition that, despite its prevalence, does not have any intermediary guidelines for pharmacological treatment between palliative care and invasive surgery. The development of treatments for the IVD is complicated by the variety of resident cell types needed to maintain the regionally distinct structural properties of the IVD that permit the safe, complex motions of the spine. Osteoporosis of the spine increases the risk of vertebral bone fracture that can increase the incidence of back pain. Fortunately, there are a variety of pharmacological treatments for osteoporosis that target osteoblasts, osteoclasts and/or osteocytes to build bone and prevent vertebral fracture. Of particular note, clinical and preclinical studies suggest that commonly prescribed osteoporosis drugs like bisphosphonates, intermittent parathyroid hormone, anti-sclerostin antibody, selective estrogen receptor modulators and anti-receptor activator of nuclear factor-kappa B ligand inhibitor denosumab may also relieve back pain. Here, we cite clinical and preclinical studies and include unpublished data to support the argument that a subset of these therapeutics for osteoporosis may alleviate low back pain by also targeting the IVD.

Characterization of the Nucleus Pulposus Progenitor Cells via Spatial Transcriptomics

Loss of refreshment in nucleus pulposus (NP) cellularity leads to intervertebral disc (IVD) degeneration. Nevertheless, the cellular sequence of NP cell differentiation remains unclear, although an increasing body of literature has identified markers of NP progenitor cells (NPPCs). Notably, due to their fragility, the physical enrichment of NP-derived cells has limited conventional transcriptomic approaches in multiple studies. To overcome this limitation, a spatially resolved transcriptional atlas of the mouse IVD is generated via the 10x Genomics Visium platform dividing NP spots into two clusters. Based on this, most reported NPPC-markers, including Cathepsin K (Ctsk), are rare and predominantly located within the NP-outer subset. Cell lineage tracing further evidence that a small number of Ctsk-expressing cells generate the entire adult NP tissue. In contrast, Tie2, which has long suggested labeling NPPCs, is actually neither expressed in NP subsets nor labels NPPCs and their descendants in mouse models; consistent with this, an in situ sequencing (ISS) analysis validated the absence of Tie2 in NP tissue. Similarly, no Tie2-cre-mediated labeling of NPPCs is observed in an IVD degenerative mouse model. Altogether, in this study, the first spatial transcriptomic map of the IVD is established, thereby providing a public resource for bone biology.

In vitro and in vivo models define a molecular signature reference for human embryonic notochordal cells

Understanding the emergence of human notochordal cells (NC) is essential for the development of regenerative approaches. We present a comprehensive investigation into the specification and generation of bona fide NC using a straightforward pluripotent stem cell (PSC)-based system benchmarked with human fetal notochord. By integrating in vitro and in vivo transcriptomic data at single-cell resolution, we establish an extended molecular signature and overcome the limitations associated with studying human notochordal lineage at early developmental stages. We show that TGF-β inhibition enhances the yield and homogeneity of notochordal lineage commitment in vitro. Furthermore, this study characterizes regulators of cell-fate decision and matrisome enriched in the notochordal niche. Importantly, we identify specific cell-surface markers opening avenues for differentiation refinement, NC purification, and functional studies. Altogether, this study provides a human notochord transcriptomic reference that will serve as a resource for notochord identification in human systems, diseased-tissues modeling, and facilitating future biomedical research.

Senescent response in inner annulus fibrosus cells in response to TNFα, H2O2, and TNFα-induced nucleus pulposus senescent secretome

Senescence, particularly in the nucleus pulposus (NP) cells, has been implicated in the pathogenesis of disc degeneration, however, the mechanism(s) of annulus fibrosus (AF) cell senescence is still not well understood. Both TNFα and H2O2, have been implicated as contributors to the senescence pathways, and their levels are increased in degenerated discs when compared to healthy discs. Thus, the objective of this study is to identify factor(s) that induces inner AF (iAF) cell senescence. Under TNFα exposure, at a concentration previously shown to induce senescence in NP cells, bovine iAF cells did not undergo senescence, indicated by their ability to continue to proliferate as demonstrated by Ki67 staining and growth curves and lack of expression of the senescent markers, p16 and p21. The lack of senescent response occurred even though iAF express higher levels of TNFR1 than NP cells. Interestingly, iAF cells showed no increase in intracellular ROS or secreted H2O2 in response to TNFα which contrasted to NP cells that did. Following TNFα treatment, only iAF cells had increased expression of the superoxide scavengers SOD1 and SOD2 whereas NP cells had increased NOX4 gene expression, an enzyme that can generate H2O2. Treating iAF cells with low dose H2O2 (50 μM) induced senescence, however unlike TNFα, H2O2 did not induce degenerative-like changes as there was no difference in COL2, ACAN, MMP13, or IL6 gene expression or number of COL2 and ACAN immunopositive cells compared to untreated controls. The latter result suggests that iAF cells may have distinct degenerative and senescent phenotypes. To evaluate paracrine signalling by senescent NP cells, iAF and TNFα-treated NP cells were co-cultured. In contact co-culture the NP cells induced iAF senescence. Thus, senescent NP cells may secrete soluble factors that induce degenerative and senescent changes within the iAF. This may contribute to a positive feedback loop of disc degeneration. It is possible these factors may include H2O2 and cytokines (such as TNFα). Further studies will investigate if human disc cells respond similarly.

Transcriptional profiling of human cartilage endplate cells identifies novel genes and cell clusters underlying degenerated and non-degenerated phenotypes

BACKGROUND

Low back pain is a leading cause of disability worldwide and is frequently attributed to intervertebral disc (IVD) degeneration. Though the contributions of the adjacent cartilage endplates (CEP) to IVD degeneration are well documented, the phenotype and functions of the resident CEP cells are critically understudied. To better characterize CEP cell phenotype and possible mechanisms of CEP degeneration, bulk and single-cell RNA sequencing of non-degenerated and degenerated CEP cells were performed.

METHODS

Human lumbar CEP cells from degenerated (Thompson grade ≥ 4) and non-degenerated (Thompson grade ≤ 2) discs were expanded for bulk (N=4 non-degenerated, N=4 degenerated) and single-cell (N=1 non-degenerated, N=1 degenerated) RNA sequencing. Genes identified from bulk RNA sequencing were categorized by function and their expression in non-degenerated and degenerated CEP cells were compared. A PubMed literature review was also performed to determine which genes were previously identified and studied in the CEP, IVD, and other cartilaginous tissues. For single-cell RNA sequencing, different cell clusters were resolved using unsupervised clustering and functional annotation. Differential gene expression analysis and Gene Ontology, respectively, were used to compare gene expression and functional enrichment between cell clusters, as well as between non-degenerated and degenerated CEP samples.

RESULTS

Bulk RNA sequencing revealed 38 genes were significantly upregulated and 15 genes were significantly downregulated in degenerated CEP cells relative to non-degenerated cells (|fold change| ≥ 1.5). Of these, only 2 genes were previously studied in CEP cells, and 31 were previously studied in the IVD and other cartilaginous tissues. Single-cell RNA sequencing revealed 11 unique cell clusters, including multiple chondrocyte and progenitor subpopulations with distinct gene expression and functional profiles. Analysis of genes in the bulk RNA sequencing dataset showed that progenitor cell clusters from both samples were enriched in "non-degenerated" genes but not "degenerated" genes. For both bulk- and single-cell analyses, gene expression and pathway enrichment analyses highlighted several pathways that may regulate CEP degeneration, including transcriptional regulation, translational regulation, intracellular transport, and mitochondrial dysfunction.

CONCLUSIONS

This thorough analysis using RNA sequencing methods highlighted numerous differences between non-degenerated and degenerated CEP cells, the phenotypic heterogeneity of CEP cells, and several pathways of interest that may be relevant in CEP degeneration.

Serglycin secreted by late-stage nucleus pulposus cells is a biomarker of intervertebral disc degeneration

Intervertebral disc degeneration is a natural process during aging and a leading cause of lower back pain. Here, we generate a comprehensive atlas of nucleus pulposus cells using single-cell RNA-seq analysis of human nucleus pulposus tissues (three males and four females, age 41.14 ± 18.01 years). We identify fibrotic late-stage nucleus pulposus cells characterized by upregulation of serglycin expression which facilitate the local inflammatory response by promoting the infiltration of inflammatory cytokines and macrophages. Finally, we discover that daphnetin, a potential serglycin ligand, substantially mitigates the local inflammatory response by downregulating serglycin expression in an in vivo mouse model, thus alleviating intervertebral disc degeneration. Taken together, we identify late-stage nucleus pulposus cells and confirm the potential mechanism by which serglycin regulates intervertebral disc degeneration. Our findings indicate that serglycin is a latent biomarker of intervertebral disc degeneration and may contribute to development of diagnostic and therapeutic strategies.

Intervertebral disc human nucleus pulposus cells associated with back pain trigger neurite outgrowth in vitro and pain behaviors in rats

Low back pain (LBP) is often associated with the degeneration of human intervertebral discs (IVDs). However, the pain-inducing mechanism in degenerating discs remains to be elucidated. Here, we identified a subtype of locally residing human nucleus pulposus cells (NPCs), generated by certain conditions in degenerating discs, that was associated with the onset of discogenic back pain. Single-cell transcriptomic analysis of human tissues showed a strong correlation between a specific cell subtype and the pain condition associated with the human degenerated disc, suggesting that they are pain-triggering. The application of IVD degeneration-associated exogenous stimuli to healthy NPCs in vitro recreated a pain-associated phenotype. These stimulated NPCs activated functional human iPSC-derived sensory neuron responses in an in vitro organ-chip model. Injection of stimulated NPCs into the healthy rat IVD induced local inflammatory responses and increased cold sensitivity and mechanical hypersensitivity. Our findings reveal a previously uncharacterized pain-inducing mechanism mediated by NPCs in degenerating IVDs. These findings could aid in the development of NPC-targeted therapeutic strategies for the clinically unmet need to attenuate discogenic LBP.

Comparative histopathological analysis of age-associated intervertebral disc degeneration in CD-1 and C57BL/6 mice: Anatomical and sex-based differences

Background

Intervertebral disc (IVD) degeneration is a major contributor to back pain and disability. The cause of IVD degeneration is multifactorial, with no disease-modifying treatments. Mouse models are commonly used to study IVD degeneration; however, the effects of anatomical location, strain, and sex on the progression of age-associated degeneration are poorly understood.

Methods

A longitudinal study was conducted to characterize age-, anatomical-, and sex-specific differences in IVD degeneration in two commonly used strains of mice, C57BL/6 and CD-1. Histopathological evaluation of the cervical, thoracic, lumbar, and caudal regions of mice at 6, 12, 20, and 24 months of age was conducted by two blinded observers at each IVD for the nucleus pulposus (NP), annulus fibrosus (AF), and the NP/AF boundary compartments, enabling analysis of scores by tissue compartment, summed scores for each IVD, or averaged scores for each anatomical region.

Results

C57BL/6 mice displayed mild IVD degeneration until 24 months of age; at this point, the lumbar spine demonstrated the most degeneration compared to other regions. Degeneration was detected earlier in the CD-1 mice (20 months of age) in both the thoracic and lumbar spine. In CD-1 mice, moderate to severe degeneration was noted in the cervical spine at all time points assessed. In both strains, age-associated IVD degeneration in the thoracic and lumbar spine was associated with increased histopathological scores in all IVD compartments. In both strains, minimal degeneration was detected in caudal IVDs out to 24 months of age. Both C57BL/6 and CD-1 mice displayed sex-specific differences in the presentation and progression of age-associated IVD degeneration.

Conclusions

These results showed that the progression and severity of age-associated degeneration in mouse models is associated with marked differences based on anatomical region, sex, and strain. This information provides a fundamental baseline characterization for users of mouse models to enable effective and appropriate experimental design, interpretation, and comparison between studies.

Spatially multicellular variability of intervertebral disc degeneration by comparative single-cell analysis

Previous studies have revealed cellular heterogeneity in intervertebral discs (IVDs). However, the cellular and molecular alteration patterns of cell populations during degenerative progression remain to be fully elucidated. To illustrate the cellular and molecular alteration of cell populations in intervertebral disc degeneration (IDD), we perform single cell RNA sequencing on cells from four anatomic sites of healthy and degenerative goat IVDs. EGLN3+ StressCs, TGFBR3+ HomCs and GPRC5A+ RegCs exhibit the characteristics associated with resistance to stress, maintaining homeostasis and repairing, respectively. The frequencies and signatures of these cell clusters fluctuate with IDD. Notably, the chondrogenic differentiation programme of PROCR+ progenitor cells is altered by IDD, while notochord cells turn to stemness exhaustion. In addition, we characterise CAV1+ endothelial cells that communicate with chondrocytes through multiple signalling pathways in degenerative IVDs. Our comprehensive analysis identifies the variability of key cell clusters and critical regulatory networks responding to IDD, which will facilitate in-depth investigation of therapeutic strategies for IDD.

Recommendations for intervertebral disc notochordal cell investigation: From isolation to characterization

Background

Lineage-tracing experiments have established that the central region of the mature intervertebral disc, the nucleus pulposus (NP), develops from the embryonic structure called "the notochord". However, changes in the cells derived from the notochord which form the NP (i.e., notochordal cells [NCs]), in terms of their phenotype and functional identity from early developmental stages to skeletal maturation are less understood. These key issues require further investigation to better comprehend the role of NCs in homeostasis and degeneration as well as their potential for regeneration. Progress in utilizing NCs is currently hampered due to poor consistency and lack of consensus methodology for in vitro NC extraction, manipulation, and characterization.

Methods

Here, an international group has come together to provide key recommendations and methodologies for NC isolation within key species, numeration, in vitro manipulation and culture, and characterization.

Results

Recommeded protocols are provided for isolation and culture of NCs. Experimental testing provided recommended methodology for numeration of NCs. The issues of cryopreservation are demonstrated, and a pannel of immunohistochemical markers are provided to inform NC characterization.

Conclusions

Together we hope this article provides a road map for in vitro studies of NCs to support advances in research into NC physiology and their potential in regenerative therapies.

Chemonucleolysis combined with dynamic loading for inducing degeneration in bovine caudal intervertebral discs

Chemonucleolysis has become an established method of producing whole organ culture models of intervertebral disc (IVD) degeneration. However, the field needs more side-by-side comparisons of the degenerative effects of the major enzymes used in chemonucleolysis towards gaining a greater understanding of how these organ culture models mimic the wide spectrum of characteristics observed in human degeneration. In the current work we induced chemonucleolysis in bovine coccygeal IVDs with 100 µL of papain (65 U/mL), chondroitinase ABC (chABC, 5 U/mL), or collagenase II (col'ase, 0.5 U/mL). Each enzyme was applied in a concentration projected to produce moderate levels of degeneration. After 7 days of culture with daily dynamic physiological loading (0.02-0.2 MPa, 0.2 Hz, 2 h), the cellular, biochemical and histological properties of the IVDs were evaluated in comparison to a PBS-injected control. Papain and collagenase, but not chABC, produced macroscopic voids in the tissues. Compared to day 0 intact IVDs, papain induced the greatest magnitude glycosaminoglycan (GAG) loss compared to chABC and col'ase. Papain also induced the greatest height loss (3%), compared to 0.7%, 1.2% and 0.4% for chABC, col'ase, and PBS, respectively. Cell viability in the region adjacent to papain and PBS-injection remained at nearly 100% over the 7-day culture period, whereas it was reduced to 60%-70% by chABC and col'ase. Generally, enzyme treatment tended to downregulate gene expression for major ECM markers, type I collagen (COL1), type II collagen (COL2), and aggrecan (ACAN) in the tissue adjacent to injection. However, chABC treatment induced an increase in COL2 gene expression, which was significant compared to the papain treated group. In general, papain and col'ase treatment tended to recapitulate aspects of advanced IVD degeneration, whereas chABC treatment captured aspects of early-stage degeneration. Chemonucleolysis of whole bovine IVDs is a useful tool providing researchers with a robust spectrum of degenerative changes and can be utilized for examination of therapeutic interventions.

Development and Degeneration of the Intervertebral Disc-Insights from Across Species

Back pain caused by intervertebral disc (IVD) degeneration has a major socio-economic impact in humans, yet historically has received minimal attention in species other than humans, mice and dogs. However, a general growing interest in this unique organ prompted the expansion of IVD research in rats, rabbits, cats, horses, monkeys, and cows, further illuminating the complex nature of the organ in both healthy and degenerative states. Application of recent biotechnological advancements, including single cell RNA sequencing and complex data analysis methods has begun to explain the shifting inflammatory signaling, variation in cellular subpopulations, differential gene expression, mechanical loading, and metabolic stresses which contribute to age and stress related degeneration of the IVD. This increase in IVD research across species introduces a need for chronicling IVD advancements and tissue biomarkers both within and between species. Here we provide a comprehensive review of recent single cell RNA sequencing data alongside existing case reports and histo/morphological data to highlight the cellular complexity and metabolic challenges of this unique organ that is of structural importance for all vertebrates.

Alterations of bovine nucleus pulposus cells with aging

Aging is one of the major etiological factors driving intervertebral disc (IVD) degeneration, the main cause of low back pain. The nucleus pulposus (NP) includes a heterogeneous cell population, which is still poorly characterized. Here, we aimed to uncover main alterations in NP cells with aging. For that, bovine coccygeal discs from young (12 months) and old (10-16 years old) animals were dissected and primary NP cells were isolated. Gene expression and proteomics of fresh NP cells were performed. NP cells were labelled with propidium iodide and analysed by flow cytometry for the expression of CD29, CD44, CD45, CD146, GD2, Tie2, CD34 and Stro-1. Morphological cell features were also dissected by imaging flow cytometry. Elder NP cells (up-regulated bIL-6 and bMMP1 gene expression) presented lower percentages of CD29+, CD44+, CD45+ and Tie2+ cells compared with young NP cells (upregulated bIL-8, bCOL2A1 and bACAN gene expression), while GD2, CD146, Stro-1 and CD34 expression were maintained with age. NP cellulome showed an upregulation of proteins related to endoplasmic reticulum (ER) and melanosome independently of age, whereas proteins upregulated in elder NP cells were also associated with glycosylation and disulfide bonds. Flow cytometry analysis of NP cells disclosed the existence of 4 subpopulations with distinct auto-fluorescence and size with different dynamics along aging. Regarding cell morphology, aging increases NP cell area, diameter and vesicles. These results contribute to a better understanding of NP cells aging and highlighting potential anti-aging targets that can help to mitigate age-related disc disease.

Loss of function mutation in Ank causes aberrant mineralization and acquisition of osteoblast-like-phenotype by the cells of the intervertebral disc

Pathological mineralization of intervertebral disc is debilitating and painful and linked to disc degeneration in a subset of human patients. An adenosine triphosphate efflux transporter, progressive ankylosis (ANK) is a regulator of extracellular inorganic pyrophosphate levels and plays an important role in tissue mineralization. However, the function of ANK in intervertebral disc has not been fully explored. Herein we analyzed the spinal phenotype of Ank mutant mice (ank/ank) with attenuated ANK function. Micro-computed tomography and histological analysis showed that loss of ANK function results in the aberrant annulus fibrosus mineralization and peripheral disc fusions with cranial to caudal progression in the spine. Vertebrae in ank mice exhibit elevated cortical bone mass and increased tissue non-specific alkaline phosphatase-positive endplate chondrocytes with decreased subchondral endplate porosity. The acellular dystrophic mineral inclusions in the annulus fibrosus were localized adjacent to apoptotic cells and cells that acquired osteoblast-like phenotype. Fourier transform infrared spectral imaging showed that the apatite mineral in the outer annulus fibrosus had similar chemical composition to that of vertebral bone. Transcriptomic analysis of annulus fibrosus and nucleus pulposus tissues showed changes in several biological themes with a prominent dysregulation of BMAL1/CLOCK circadian regulation. The present study provides new insights into the role of ANK in the disc tissue compartments and highlights the importance of local inorganic pyrophosphate metabolism in inhibiting the mineralization of this important connective tissue.

Intervertebral Disc Progenitors: Lessons Learned from Single-Cell RNA Sequencing and the Role in Intervertebral Disc Regeneration

The tremendous personal and economic burden worldwide caused by low back pain (LBP) has been surging in recent years. While intervertebral disc degeneration (IVDD) is the leading cause of LBP and vast efforts have been made to develop effective therapies, this problem is far from being resolved, as most treatments, such as painkillers and surgeries, mainly focus on relieving the symptoms rather than reversing the cause of IVDD. However, as stem/progenitor cells possess the potential to regenerate IVD, a deeper understanding of the early development and role of these cells could help to improve the effectiveness of stem/progenitor cell therapy in treating LBP. Single-cell RNA sequencing results provide fresh insights into the heterogeneity and development patterns of IVD progenitors; additionally, we compare mesenchymal stromal cells and IVD progenitors to provide a clearer view of the optimal cell source proposed for IVD regeneration.

The Breast Cancer Single-Cell Atlas: Defining cellular heterogeneity within model cell lines and primary tumors to inform disease subtype, stemness, and treatment options

PURPOSE

Breast Cancer (BC) is the most diagnosed cancer in women; however, through significant research, relative survival rates have significantly improved. Despite progress, there remains a gap in our understanding of BC subtypes and personalized treatments. This manuscript characterized cellular heterogeneity in BC cell lines through scRNAseq to resolve variability in subtyping, disease modeling potential, and therapeutic targeting predictions.

METHODS

We generated a Breast Cancer Single-Cell Cell Line Atlas (BSCLA) to help inform future BC research. We sequenced over 36,195 cells composed of 13 cell lines spanning the spectrum of clinical BC subtypes and leveraged publicly available data comprising 39,214 cells from 26 primary tumors.

RESULTS

Unsupervised clustering identified 49 subpopulations within the cell line dataset. We resolve ambiguity in subtype annotation comparing expression of Estrogen Receptor, Progesterone Receptor, and Human Epidermal Growth Factor Receptor 2 genes. Gene correlations with disease subtype highlighted S100A7 and MUCL1 overexpression in HER2 + cells as possible cell motility and localization drivers. We also present genes driving populational drifts to generate novel gene vectors characterizing each subpopulation. A global Cancer Stem Cell (CSC) scoring vector was used to identify stemness potential for subpopulations and model multi-potency. Finally, we overlay the BSCLA dataset with FDA-approved targets to identify to predict the efficacy of subpopulation-specific therapies.

CONCLUSION

The BSCLA defines the heterogeneity within BC cell lines, enhancing our overall understanding of BC cellular diversity to guide future BC research, including model cell line selection, unintended sample source effects, stemness factors between cell lines, and cell type-specific treatment response.

Comparison of transformations for single-cell RNA-seq data

The count table, a numeric matrix of genes × cells, is the basic input data structure in the analysis of single-cell RNA-sequencing data. A common preprocessing step is to adjust the counts for variable sampling efficiency and to transform them so that the variance is similar across the dynamic range. These steps are intended to make subsequent application of generic statistical methods more palatable. Here, we describe four transformation approaches based on the delta method, model residuals, inferred latent expression state and factor analysis. We compare their strengths and weaknesses and find that the latter three have appealing theoretical properties; however, in benchmarks using simulated and real-world data, it turns out that a rather simple approach, namely, the logarithm with a pseudo-count followed by principal-component analysis, performs as well or better than the more sophisticated alternatives. This result highlights limitations of current theoretical analysis as assessed by bottom-line performance benchmarks.

Single-cell RNA sequencing in orthopedic research

Although previous RNA sequencing methods have been widely used in orthopedic research and have provided ideas for therapeutic strategies, the specific mechanisms of some orthopedic disorders, including osteoarthritis, lumbar disc herniation, rheumatoid arthritis, fractures, tendon injuries, spinal cord injury, heterotopic ossification, and osteosarcoma, require further elucidation. The emergence of the single-cell RNA sequencing (scRNA-seq) technique has introduced a new era of research on these topics, as this method provides information regarding cellular heterogeneity, new cell subtypes, functions of novel subclusters, potential molecular mechanisms, cell-fate transitions, and cell‒cell interactions that are involved in the development of orthopedic diseases. Here, we summarize the cell subpopulations, genes, and underlying mechanisms involved in the development of orthopedic diseases identified by scRNA-seq, improving our understanding of the pathology of these diseases and providing new insights into therapeutic approaches.

Macrophages and Intervertebral Disc Degeneration

The intervertebral disc (IVD) aids in motion and acts to absorb energy transmitted to the spine. With little inherent regenerative capacity, degeneration of the intervertebral disc results in intervertebral disc disease, which contributes to low back pain and significant disability in many individuals. Increasing evidence suggests that IVD degeneration is a disease of the whole joint that is associated with significant inflammation. Moreover, studies show elevated macrophage accumulation within the IVD with increasing levels of disease severity; however, we still need to understand the roles, be they causative or consequential, of macrophages during the degenerative process. In this narrative review, we discuss hallmarks of IVD degeneration, showcase evidence of macrophage involvement during disc degeneration, and explore burgeoning research aimed at understanding the molecular pathways regulating macrophage functions during intervertebral disc degeneration.

Single-cell RNA sequencing reveals resident progenitor and vascularization-associated cell subpopulations in rat annulus fibrosus

Background

One of the main causes of low back pain is intervertebral disc degeneration (IDD). Annulus fibrosus (AF) is important for the integrity and functions of the intervertebral disc (IVD). However, the resident functional cell components such as progenitors and vascularization-associated cells in AF are yet to be fully identified.

Purpose

Identification of functional AF cell subpopulations including resident progenitors and vascularization-associated cells.

Methods

In this study, the single-cell RNA sequencing data of rat IVDs from a public database were analyzed using Seurat for cell clustering, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) for functional analysis, StemID for stem cell identification, Monocle and RNA velocity for pseudotime differentiation trajectory validation, single-cell regulatory network inference and clustering (SCENIC) for gene regulatory network (GRN) analysis, and CellChat for cell-cell interaction analysis. Immunostaining on normal and degenerated rat IVDs, as well as human AF, was used for validations.

Results

From the data analysis, seven AF cell clusters were identified, including two newly discovered functional clusters, the Grem1 ⁺ subpopulation and the Lum ⁺ ​subpopulation. The Grem1 ⁺ subpopulation had progenitor characteristics, while the Lum ⁺ ​subpopulation was associated with vascularization during IDD. The GRN analysis showed that Sox9 and Id1 were among the key regulators in the Grem1 ⁺ subpopulation, and Nr2f2 and Creb5 could be responsible for the vascularization function in the Lum ⁺ ​subpopulation. Cell-cell interaction analysis revealed highly regulated cellular communications between these cells, and multiple signaling networks including PDGF and MIF signaling pathways were involved in the interactions.

Conclusions

Our results revealed two new functional AF cell subpopulations, with stemness and vascularization induction potential, respectively.

The Translational potential of this article

These findings complement our knowledge about IVDs, especially the AF, and in return provide potential cell source and regulation targets for IDD treatment and tissue repair. The existence of the cell subpopulations was also validated in human AF, which strengthen the clinical relevance of the findings.

Single Cell RNA-Sequence Analyses Reveal Uniquely Expressed Genes and Heterogeneous Immune Cell Involvement in the Rat Model of Intervertebral Disc Degeneration

Intervertebral disc (IVD) degeneration is characterized by a loss of cellularity, and changes in cell-mediated activity that drives anatomic changes to IVD structure. In this study, we used single-cell RNA-sequencing analysis of degenerating tissues of the rat IVD following lumbar disc puncture. Two control, uninjured IVDs (L2-3, L3-4) and two degenerated, injured IVDs (L4-5, L5-6) from each animal were examined either at the two- or eight-week post-operative time points. The cells from these IVDs were extracted and transcriptionally profiled at the single-cell resolution. Unsupervised cluster analysis revealed the presence of four known cell types in both non-degenerative and degenerated IVDs based on previously established gene markers: IVD cells, endothelial cells, myeloid cells, and lymphoid cells. As a majority of cells were associated with the IVD cell cluster, sub-clustering was used to further identify the cell populations of the nucleus pulposus, inner and outer annulus fibrosus. The most notable difference between control and degenerated IVDs was the increase of myeloid and lymphoid cells in degenerated samples at two- and eight-weeks post-surgery. Differential gene expression analysis revealed multiple distinct cell types from the myeloid and lymphoid lineages, most notably macrophages and B lymphocytes, and demonstrated a high degree of immune specificity during degeneration. In addition to the heterogenous infiltrating immune cell populations in the degenerating IVD, the increased number of cells in the AF sub-cluster expressing Ngf and Ngfr, encoding for p75NTR, suggest that NGF signaling may be one of the key mediators of the IVD crosstalk between immune and neuronal cell populations. These findings provide the basis for future work to understand the involvement of select subsets of non-resident cells in IVD degeneration.

Genipin-crosslinked fibrin seeded with oxidized alginate microbeads as a novel composite biomaterial strategy for intervertebral disc cell therapy

Discectomy procedures alleviate disability caused by intervertebral disc (IVD) herniation, but do not repair herniation-induced annulus fibrosus (AF) defects. Cell therapy shows promise for IVD repair, yet cell delivery biomaterials capable of sealing AF defects and restoring biomechanical function have poor biological performance. To balance the biomechanical and biological demands of IVD cell delivery biomaterials, we engineered an injectable composite biomaterial using cell-laden, degradable oxidized alginate (OxAlg) microbeads (MBs) to deliver AF cells within high-modulus genipin-crosslinked fibrin (FibGen) hydrogels (FibGen + MB composites). Conceptually, the high-modulus FibGen would immediately stabilize injured IVDs, while OxAlg MBs would protect and release cells required for long-term healing. We first showed that AF cells microencapsulated in OxAlg MBs maintained high viability and, upon release, displayed phenotypic AF cell morphology and gene expression. Next, we created cell-laden FibGen + MB composites and demonstrated that OxAlg MBs functionalized with RGD peptides (MB-RGD) minimized AF cell apoptosis and retained phenotypic gene expression. Further, we showed that cell-laden FibGen + MB composites are biomechanically stable and promote extracellular matrix (ECM) synthesis in long-term in vitro culture. Lastly, we evaluated cell-laden FibGen + MB-RGD composites in a long-term bovine caudal IVD organ culture bioreactor and found that composites had low herniation risk, provided superior biomechanical and biological repair to discectomy controls, and retained anabolic cells within the IVD injury space. This novel injectable composite hydrogel strategy shows promise as an IVD cell delivery sealant with potentially broad applications for its capacity to balance biomechanical and biological performance.

Single-cell atlas unveils cellular heterogeneity and novel markers in human neonatal and adult intervertebral discs

The origin, composition, distribution, and function of cells in the human intervertebral disc (IVD) have not been fully understood. Here, cell atlases of both human neonatal and adult IVDs have been generated and further assessed by gene ontology pathway enrichment, pseudo-time trajectory, histology, and immunofluorescence. Comparison of cell atlases revealed the presence of two subpopulations of notochordal cells (NCs) and their associated markers in both the neonatal and adult IVDs. Developmental trajectories predicted 7 different cell states that describe the developmental process from neonatal to adult cells in IVD and analyzed the NC’s role in the IVD development. A high heterogeneity and gradual transition of annulus fibrosus cells (AFCs) in the neonatal IVD was detected and their potential relevance in IVD development assessed. Collectively, comparing single-cell atlases between neonatal and adult IVDs delineates the landscape of IVD cell biology and may help discover novel therapeutic targets for IVD degeneration.

•

Compared scRNA-seq between human neonatal and adult IVD

•

Identified two notochordal cell populations in adults and their novel markers

•

Notochordal cells preserved their identity and functions into adulthood

•

Unveiled heterogeneity of nucleus pulposus and annulus fibrosus cells in human IVD

IGFBP-6: At the Crossroads of Immunity, Tissue Repair and Fibrosis

Insulin-like growth factors binding protein-6 (IGFBP-6) is involved in a relevant number of cellular activities and represents an important factor in the immune response, particularly in human dendritic cells (DCs). Over the past several years, significant insights into the IGF-independent effects of IGFBP-6 were discovered, such as the induction of chemotaxis, capacity to increase oxidative burst and neutrophils degranulation, ability to induce metabolic changes in DCs, and, more recently, the regulation of the Sonic Hedgehog (SHH) signaling pathway during fibrosis. IGFBP-6 has been implicated in different human diseases, and it plays a rather controversial role in the biology of tumors. Notably, well established relationships between immunity, stroma activity, and fibrosis are prognostic and predictive of response to cancer immunotherapy. This review aims at describing the current understanding of mechanisms that link IGFBP-6 and fibrosis development and at highlighting the multiple roles of IGFBP-6 to provide an insight into evolutionarily conserved mechanisms that can be relevant for inflammation, tumor immunity, and immunological diseases.

Development, Pathogenesis, and Regeneration of the Intervertebral Disc: Current and Future Insights Spanning Traditional to Omics Methods

The intervertebral disc (IVD) is the fibrocartilaginous joint located between each vertebral body that confers flexibility and weight bearing capabilities to the spine. The IVD plays an important role in absorbing shock and stress applied to the spine, which helps to protect not only the vertebral bones, but also the brain and the rest of the central nervous system. Degeneration of the IVD is correlated with back pain, which can be debilitating and severely affects quality of life. Indeed, back pain results in substantial socioeconomic losses and healthcare costs globally each year, with about 85% of the world population experiencing back pain at some point in their lifetimes. Currently, therapeutic strategies for treating IVD degeneration are limited, and as such, there is great interest in advancing treatments for back pain. Ideally, treatments for back pain would restore native structure and thereby function to the degenerated IVD. However, the complex developmental origin and tissue composition of the IVD along with the avascular nature of the mature disc makes regeneration of the IVD a uniquely challenging task. Investigators across the field of IVD research have been working to elucidate the mechanisms behind the formation of this multifaceted structure, which may identify new therapeutic targets and inform development of novel regenerative strategies. This review summarizes current knowledge base on IVD development, degeneration, and regenerative strategies taken from traditional genetic approaches and omics studies and discusses the future landscape of investigations in IVD research and advancement of clinical therapies.

Decoding the annulus fibrosus cell atlas by scRNA-seq to develop an inducible composite hydrogel: A novel strategy for disc reconstruction

Low back pain is one of the most serious public health problems worldwide and the major clinical manifestation of intervertebral disc degeneration (IVDD). The key pathological change during IVDD is dysfunction of the annulus fibrosus (AF). However, due to the lack of an in-depth understanding of AF biology, the methods to reconstruct the AF are very limited. In this study, the mice AF cell atlas were decoded by single-cell RNA sequencing to provide a guide for AF reconstruction. The results first identify a new population of AF cells, fibrochondrocyte-like AF cells, which synthesize both collagen I and collagen II and are potential functional cells for AF reconstruction. According to the dual features of the AF extracellular matrix, a composite hydrogel based on the acylation of methacrylated silk fibroin with methacrylated hyaluronic acid was produced. To obtain the ability to stimulate differentiation, the composite hydrogels were combined with a fibrochondrocyte-inducing supplement. Finally, reconstruction of the AF defects, by the novel AF stem cell-loaded composite hydrogel, could be observed, its amount of chondroid matrices recovered to 31.7% of AF aera which is significantly higher than that in other control groups. In summary, this study decodes the AF cell atlas, based on which a novel strategy for AF reconstruction is proposed.

•

There are 10 populations of cells in the annulus fibrosus (AF), as decoded by single cell RNA sequencing.

•

Lineage tracing shows the route of migration and differentiation of annulus fibrosus-derived stem cells (AFSCs).

•

A new population of AF cells, fibrochondrocyte-like AF cells, was identified.

•

Both the fibrinoid and chondroid matrices of AF are reconstructed by the novel AFSCs-loaded composite hydrogel.