Background: The nucleus pulposus is a constituent structure of the human intervertebral disc, and its degeneration can cause intervertebral disc degeneration (IDD). However, the cellular and molecular mechanisms involved remain elusive.Methods: Through bioinformatics analysis, the single-cell transcriptome sequencing expression profiles of human normal nucleus pulposus (NNP) cells and human degenerative nucleus pulposus (DNP) cells were compared to clarify the transcriptome differential expression profiles of human NNP and DNP. The single-cell sequencing results of the two samples were analyzed using bioinformatics methods to compare the differences in histiocytosis between human NNP and DNP, map the histiocytes of NNP and DNP, perform cell differentiation trajectories for the cell populations of interest and predict cell function, and explore their heterogeneity by pathway analysis and Gene Ontology analysis.Results: Nine cell types were identified, which were chondrocyte 1, chondrocyte 2, chondrocyte 3, chondrocyte 4, chondrocyte 5, endothelial, macrophage, neutrophil, and T cells. Analysis of the proportion of chondrocytes in different tissues revealed that chondrocyte 1 accounted for a higher proportion of NNP cells and highly expressed COL2A1 compared with DNP cells; chondrocyte 2, chondrocyte 3, chondrocyte 4, and chondrocyte 5 accounted for a higher proportion of DNP cells compared with NNP cells. Among them, chondrocyte 2 was an inhibitory calcified chondrocyte with high expression of MGP, chondrocytes 3 were fibrochondrocytes with high expression of COL1A1, chondrocytes 4 were chondrocytes that highly express pain inflammatory genes such as PTGES, and chondrocytes 5 were calcified chondrocytes with high expression of FN1 (chondrocytes 4 and chondrocytes 5 were found for the first time in a study of single-cell transcriptome sequencing of disc tissue). Cell trajectory analysis revealed that chondrocyte 1 was at the beginning of the trajectory and chondrocyte 3 was at the end of the trajectory, while chondrocyte 5 appeared first in the trajectory relative to chondrocyte 2 and chondrocyte 4.Conclusion: After functional identification of the specifically expressed genes in five chondrocytes, it was found that chondrocyte 1 was a chondrocyte with high expression of COL2A1, COL9A2, COL11A2, and CHRDL2 in a high proportion of NNP cells, and chondrocyte 3 was a fibrochondrocyte with high expression of COL1A1, COL6A3, COL1A2, COL3A1, AQP1, and COL15A1 in an increased proportion during nucleus pulposus cell degeneration. Through cell trajectory analysis, it was found that chondrocytes 5 specifically expressing FN1, SESN2, and GDF15 may be the key cells leading to degeneration of nucleus pulposus cells. Chondrocytes 2 expressing MGP, MT1G, and GPX3 may play a role in reversing calcification and degeneration, and chondrocytes 4 expressing PTGES, TREM1, and TIMP1 may play a role in disc degeneration pain and inflammation.

Background: Single-cell RNA sequencing (scRNA-seq) enables specific analysis of cell populations at single-cell resolution; however, there is still a lack of single-cell-level studies to characterize the dynamic and complex interactions between osteoporotic vertebral compression fractures (OVCFs) and Kümmell’s disease (KD) in the osteoimmune microenvironment. In this study, we used scRNA-seq analysis to investigate the osteoimmune microenvironment and cellular composition in OVCFs and KD.Methods: ScRNA-seq was used to perform analysis of fractured vertebral bone tissues from one OVCF and one KD patients, and a total of 8,741 single cells were captured for single-cell transcriptomic analysis. The cellularity of human vertebral bone tissue was further analyzed using uniform manifold approximation and projection. Pseudo-time analysis and gene enrichment analysis revealed the biological function of cell fate and its counterparts. CellphoneDB was used to identify the interactions between bone cells and immune cells in the osteoimmune microenvironment of human vertebral bone tissue and their potential functions.Results: A cellular profile of the osteoimmune microenvironment of human vertebral bone tissue was established, including mesenchymal stem cells (MSCs), pericytes, myofibroblasts, fibroblasts, chondrocytes, endothelial cells (ECs), granulocytes, monocytes, T cells, B cells, plasma cells, mast cells, and early erythrocytes. MSCs play an immunoregulatory function and mediate osteogenic differentiation and cell proliferation. The differentiation trajectory of osteoclasts in human vertebral bone tissue was also revealed. In addition, ECs actively participate in inflammatory infiltration and coupling with bone cells. T and B cells actively participate in regulating bone homeostasis. Finally, by identifying the interaction of ligand–receptor pairs, we found that immune cells and osteoclasts have bidirectional regulatory characteristics, have the effects of regulating bone resorption by osteoclasts and promoting bone formation, and are essential for bone homeostasis. It is also highlighted that CD8-TEM cells and osteoclasts might crosstalk via CD160–TNFRSF14 ligand–receptor interaction.Conclusion: Our analysis reveals a differential landscape of molecular pathways, population composition, and cell–cell interactions during OVCF development into KD. OVCFs exhibit a higher osteogenic differentiation capacity, owing to abundant immune cells. Conversely, KD results in greater bone resorption than bone formation due to depletion of MSCs and a relatively suppressed immune system, and this immune imbalance eventually leads to vertebral avascular necrosis. The site of action between immune cells and osteoclasts is expected to be a new therapeutic target, and these results may accelerate mechanistic and functional studies of osteoimmune cell types and specific gene action in vertebral avascular necrosis and pathological bone loss diseases, paving the way for drug discovery.

Background: Giant cell tumor of bone (GCTB) presents considerable complexity in tumor microenvironment (TME) because of its intricate intercellular heterogeneity and the presence of an immunosuppressive milieu. In order to understand the complex gene expression patterns and cell interactions in GCTB, we carried out a thorough investigation using single‐cell RNA sequencing (scRNA‐seq). Methods: We examined scRNA‐seq data from 7091 cells that were collected after surgical removal of GCTB. Following the initial quality control process, 10 separate groups of cells were distinguished, which consisted of dendritic cells, endothelial cells, macrophage cells, mast cells, monocyte cells, neutrophil cells, tumor cells, osteoclast cells, pericyte cells, and T cells. Additional analysis uncovered distinct categories within tumor‐associated macrophages (TAMs), CD8+ T cells, and CD4+ T cells. The differentiation mechanisms of TAMs, CD8+ T cells, and CD4+ T cells were explored using pseudo‐time trajectory analysis. The CellPhoneDB study revealed the interactions between various cell types within the TME of GCTB. Results: TAMs have been identified as the main infiltrating cells in GCTB. These TAMs exhibit several subtypes that are characterized by specific marker genes and functional states. The identification of several subgroups within CD8+ T cells that are involved in regulating immunological checkpoints underscores the difficulties encountered when attempting to employ immune checkpoint blockade therapy for GCTB. T cell exhaustion poses a major barrier to the efficacy of antitumor immune responses. Research suggests a strong correlation between TAMs and exhausted T cells (Texs) in the TME. The high number of regulatory T cells (Tregs) highlights the immunosuppressive nature of the immunological environment in GCTB. Significant interactions have been observed between TAMs and tumor cells, highlighting their crucial involvement in immune evasion strategies. Conclusion: This scRNA‐seq study provides a general overview of the different cellular compositions and immune interactions within GCTB. The identified subtypes and communication networks provide valuable information about the immunosuppressive environment of GCTB, laying the foundation for prospective therapeutic approaches targeting specific cell types or interactions.

Background: Osteoarthritis (OA), a refractory disease, is one of the leading contributors for disability worldwide. Since chondrocyte is the only resident cell in cartilage, this study aims to explore the roles of miR-129-3p and CPEB1 in chondrocyte apoptosis in knee joint fracture-induced OA. Methods: Cartilage was collected from 20 OA patients who underwent total knee replacement (OA group) and 20 patients with knee contusion (normal group). Then, miR-129-3p and CPEB1 levels in the cartilage were quantified by qRT-PCR. Primary rat chondrocytes in the knee were isolated and identified by toluidine blue staining and immunofluorescent staining of type II collagen. OA cellular models were induced by TNF-α treatment, in which miR-129-3p and CPEB1 expressions were assessed. Subsequently, cell viability, apoptosis, and the expression levels of apoptotic protein and caspase-3 were measured. Dual luciferase reporter assay identified the interaction between miR-129-3p and CPEB1. Results: Patients in the OA group had decreased miR-129-3p expression and increased CPEB1 expression than those in the normal group. TNF-α treatment successfully induced the OA cellular model. Downregulated miR-129-3p and upregulated CPEB1 expressions were found in OA-treated chondrocytes. miR-129-3p overexpression or CPEB1 knockdown improved chondrocyte viability and attenuated apoptosis, and vice versa. miR-129-3p negatively regulated CPEB1, thus ameliorating apoptosis and enhancing cell viability. Conclusion: miR-129-3p negatively targeted CPEB1 to facilitate chondrocyte viability and hamper apoptosis.

Worldwide, ovarian cancer (OC) is the seventh common cancer and the second most common cause of cancer death in women. Due to high rates of relapse, there is an urgent need for the identification of new targets for OC treatment. The far-upstream element binding protein 1 (FBP1) and enhancer of zeste homolog 2 (EZH2) are emerging proto-oncogenes that regulate cell proliferation and metastasis. In the present study, Oncomine data analysis demonstrated that FBP1 was closely associated with the development of OC, and The Cancer Genome Atlas (TCGA) data analysis indicated that there was a positive correlation between FBP1 and EZH2 in ovarian tissues. Moreover, we found that FBP1 knockdown suppressed tumor formation in nude mice and cisplatin resistance of OC cells, but the role of FBP1 in the cisplatin resistance of OC cells remained unclear. In addition, we verified physical binding between FBP1 and EZH2 in OC cells, and we demonstrated that FBP1 knockdown enhanced cisplatin cytotoxicity in OC cells and down-regulated EZH2 expression and trimethylation of H3K27. These results suggested that FBP1 increases cisplatin resistance of OC cells by up-regulating EZH2/H3K27me3. Thus, FBP1 is a prospective novel target for the development of OC treatment.

Corrigendum: Single-cell RNA sequencing reveals the difference in human normal and degenerative nucleus pulposus tissue profiles and cellular interactions