A systematic dissection of human primary osteoblasts in vivo at single-cell resolution

Advancing skeletal health and disease research with single-cell RNA sequencing

Orthopedic conditions have emerged as global health concerns, impacting approximately 1.7 billion individuals worldwide. However, the limited understanding of the underlying pathological processes at the cellular and molecular level has hindered the development of comprehensive treatment options for these disorders. The advent of single-cell RNA sequencing (scRNA-seq) technology has revolutionized biomedical research by enabling detailed examination of cellular and molecular diversity. Nevertheless, investigating mechanisms at the single-cell level in highly mineralized skeletal tissue poses technical challenges. In this comprehensive review, we present a streamlined approach to obtaining high-quality single cells from skeletal tissue and provide an overview of existing scRNA-seq technologies employed in skeletal studies along with practical bioinformatic analysis pipelines. By utilizing these methodologies, crucial insights into the developmental dynamics, maintenance of homeostasis, and pathological processes involved in spine, joint, bone, muscle, and tendon disorders have been uncovered. Specifically focusing on the joint diseases of degenerative disc disease, osteoarthritis, and rheumatoid arthritis using scRNA-seq has provided novel insights and a more nuanced comprehension. These findings have paved the way for discovering novel therapeutic targets that offer potential benefits to patients suffering from diverse skeletal disorders.

Epigenetic regulations of cellular senescence in osteoporosis

Osteoporosis (OP) is a prevalent age-related disease that is characterized by a decrease in bone mineral density (BMD) and systemic bone microarchitectural disorders. With age, senescent cells accumulate and exhibit the senescence-associated secretory phenotype (SASP) in bone tissue, leading to the imbalance of bone homeostasis, osteopenia, changes in trabecular bone structure, and increased bone fragility. Cellular senescence in the bone microenvironment involves osteoblasts, osteoclasts, and bone marrow mesenchymal stem cells (BMSCs), whose effects on bone homeostasis are regulated by epigenetics. Therefore, the epigenetic regulatory mechanisms of cellular senescence have received considerable attention as potential targets for preventing and treating osteoporosis. In this paper, we systematically review the mechanisms of aging-associated epigenetic regulation in osteoporosis, emphasizing the impact of epigenetics on cellular senescence, and summarize three current methods of targeting cellular senescence, which is helpful better to understand the pathogenic mechanisms of cellular senescence in osteoporosis and provides strategies for the development of epigenetic drugs for the treatment of osteoporosis.

Characterization of adult human skeletal cells in different tissues reveals a CD90+CD34+ periosteal stem/progenitor population

The periosteum plays a crucial role in bone healing and is an important source of skeletal stem and progenitor cells. Recent studies in mice indicate that diverse populations of skeletal progenitors contribute to growth, homeostasis and healing. Information about the in vivo identity and diversity of skeletal stem and progenitor cells in different compartments of the adult human skeleton is limited. In this study, we compared non-hematopoietic populations in matched tissues from the femoral head and neck of 21 human participants using spectral flow cytometry of freshly isolated cells. High-dimensional clustering analysis indicated significant differences in marker distribution between periosteum, articular cartilage, endosteum and bone marrow populations, and identified populations that were highly enriched or unique to specific tissues. Periosteum-enriched markers included CD90 and CD34. Articular cartilage, which has very poor regenerative potential, showed enrichment of multiple markers, including the PDPN+CD73+CD164+CD146- population previously reported to represent human skeletal stem cells. We further characterized periosteal populations by combining CD90 with other strongly expressed markers. CD90+CD34+ cells sorted directly from periosteum showed significant colony-forming unit fibroblasts (CFU-F) enrichment, rapid expansion, and consistent multi-lineage differentiation of clonal populations in vitro. In situ, CD90+CD34+ cells include a perivascular population in the outer layer of the periosteum and non-perivascular cells closer to the bone surface. CD90+ cells are also highly enriched for CFU-F in bone marrow and endosteum, but not articular cartilage. In conclusion, our study indicates considerable diversity in the non-hematopoietic cell populations in different tissue compartments within the adult human skeleton, and suggests that periosteal progenitor cells reside within the CD90+CD34+ population.

Profiling joint tissues at single-cell resolution: advances and insights

Advances in the profiling of human joint tissues at single-cell resolution have provided unique insights into the organization and function of these tissues in health and disease. Data generated by various single-cell technologies, including single-cell RNA sequencing and cytometry by time-of-flight, have identified the distinct subpopulations that constitute these tissues. These timely studies have provided the building blocks for the construction of single-cell atlases of joint tissues including cartilage, bone and synovium, leading to the identification of developmental trajectories, deciphering of crosstalk between cells and discovery of rare populations such as stem and progenitor cells. In addition, these studies have revealed unique pathogenetic populations that are potential therapeutic targets. The use of these approaches in synovial tissues has helped to identify how distinct cell subpopulations can orchestrate disease initiation and progression and be responsible for distinct pathological outcomes. Additionally, repair of tissues such as cartilage and meniscus remains an unmet medical need, and single-cell methodologies can be invaluable in providing a blueprint for both effective tissue-engineering strategies and therapeutic interventions for chronic joint diseases such as osteoarthritis and rheumatoid arthritis.

Importance of Metalloproteinase Enzyme Group in Selected Skeletal System Diseases

Bone tissue is a dynamic structure that is involved in maintaining the homeostasis of the body due to its multidirectional functions, such as its protective, endocrine, or immunological role. Specialized cells and the extracellular matrix (ECM) are responsible for the remodeling of specific bone structures, which alters the biomechanical properties of the tissue. Imbalances in bone-forming elements lead to the formation and progression of bone diseases. The most important family of enzymes responsible for bone ECM remodeling are matrix metalloproteinases (MMPs)-enzymes physiologically present in the body's tissues and cells. The activity of MMPs is maintained in a state of balance; disruption of their activity is associated with the progression of many groups of diseases, including those of the skeletal system. This review summarizes the current understanding of the role of MMPs in bone physiology and the pathophysiology of bone tissue and describes their role in specific skeletal disorders. Additionally, this work collects data on the potential of MMPs as bio-markers for specific skeletal diseases.

Dissection of Cellular Communication between Human Primary Osteoblasts and Bone Marrow Mesenchymal Stem Cells in Osteoarthritis at Single-Cell Resolution

Background and Objectives

Osteoblasts are derived from bone marrow mesenchymal stem cells (BMMSCs) and play important role in bone remodeling. While our previous studies have investigated the cell subtypes and heterogeneity in osteoblasts and BMMSCs separately, cell-to-cell communications between osteoblasts and BMMSCs in vivo in humans have not been characterized. The aim of this study was to investigate the cellular communication between human primary osteoblasts and bone marrow mesenchymal stem cells.

Methods and Results

To investigate the cell-to-cell communications between osteoblasts and BMMSCs and identify new cell subtypes, we performed a systematic integration analysis with our single-cell RNA sequencing (scRNA-seq) transcriptomes data from BMMSCs and osteoblasts. We successfully identified a novel preosteoblasts subtype which highly expressed ATF3, CCL2, CXCL2 and IRF1. Biological functional annotations of the transcriptomes suggested that the novel preosteoblasts subtype may inhibit osteoblasts differentiation, maintain cells to a less differentiated status and recruit osteoclasts. Ligand-receptor interaction analysis showed strong interaction between mature osteoblasts and BMMSCs. Meanwhile, we found FZD1 was highly expressed in BMMSCs of osteogenic differentiation direction. WIF1 and SFRP4, which were highly expressed in mature osteoblasts were reported to inhibit osteogenic differentiation. We speculated that WIF1 and sFRP4 expressed in mature osteoblasts inhibited the binding of FZD1 to Wnt ligand in BMMSCs, thereby further inhibiting osteogenic differentiation of BMMSCs.

Conclusions

Our study provided a more systematic and comprehensive understanding of the heterogeneity of osteogenic cells. At the single cell level, this study provided insights into the cell-to-cell communications between BMMSCs and osteoblasts and mature osteoblasts may mediate negative feedback regulation of osteogenesis process.

Carboxypeptidase E is a prognostic biomarker co-expressed with osteoblastic genes in osteosarcoma

Osteosarcoma (OS) is a rare primary malignant bone tumor in adolescents and children with a poor prognosis. The identification of prognostic genes lags far behind advancements in treatment. In this study, we identified differential genes using mRNA microarray analysis of five paired OS tissues. Hub genes, gene set enrichment analysis, and pathway analysis were performed to gain insight into the pathway alterations of OS. Prognostic genes were screened using the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) dataset, then overlapped with the differential gene dataset. The carboxypeptidase E (CPE) gene, found to be an independent risk factor, was further validated using RT-PCR and Gene Expression Omnibus (GEO) datasets. Additionally, we explored the specific expression of CPE in OS tissues by reanalyzing single-cell genomics. Interestingly, CPE was found to be co-expressed with osteoblast lineage cell clusters that expressed RUNX2, SP7, SPP1, and IBSP marker genes in OS. These results suggest that CPE could serve as a prognostic factor in osteoblastic OS and should be further investigated as a potential therapeutic target.

Regulon active landscape reveals cell development and functional state changes of human primary osteoblasts in vivo

BACKGROUND

While transcription factor (TF) regulation is known to play an important role in osteoblast development, differentiation, and bone metabolism, the molecular features of TFs in human osteoblasts at the single-cell resolution level have not yet been characterized. Here, we identified modules (regulons) of co-regulated genes by applying single-cell regulatory network inference and clustering to the single-cell RNA sequencing profiles of human osteoblasts. We also performed cell-specific network (CSN) analysis, reconstructed regulon activity-based osteoblast development trajectories, and validated the functions of important regulons both in vivo and in vitro.

RESULTS

We identified four cell clusters: preosteoblast-S1, preosteoblast-S2, intermediate osteoblasts, and mature osteoblasts. CSN analysis results and regulon activity-based osteoblast development trajectories revealed cell development and functional state changes of osteoblasts. CREM and FOSL2 regulons were mainly active in preosteoblast-S1, FOXC2 regulons were mainly active in intermediate osteoblast, and RUNX2 and CREB3L1 regulons were most active in mature osteoblasts.

CONCLUSIONS

This is the first study to describe the unique features of human osteoblasts in vivo based on cellular regulon active landscapes. Functional state changes of CREM, FOSL2, FOXC2, RUNX2, and CREB3L1 regulons regarding immunity, cell proliferation, and differentiation identified the important cell stages or subtypes that may be predominantly affected by bone metabolism disorders. These findings may lead to a deeper understanding of the mechanisms underlying bone metabolism and associated diseases.

High Bone Mass Disorders: New Insights From Connecting the Clinic and the Bench

Monogenic high bone mass (HBM) disorders are characterized by an increased amount of bone in general, or at specific sites in the skeleton. Here, we describe 59 HBM disorders with 50 known disease-causing genes from the literature, and we provide an overview of the signaling pathways and mechanisms involved in the pathogenesis of these disorders. Based on this, we classify the known HBM genes into HBM (sub)groups according to uniform Gene Ontology (GO) terminology. This classification system may aid in hypothesis generation, for both wet lab experimental design and clinical genetic screening strategies. We discuss how functional genomics can shape discovery of novel HBM genes and/or mechanisms in the future, through implementation of omics assessments in existing and future model systems. Finally, we address strategies to improve gene identification in unsolved HBM cases and highlight the importance for cross-laboratory collaborations encompassing multidisciplinary efforts to transfer knowledge generated at the bench to the clinic. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Conditional Loss of Nmp4 in Mesenchymal Stem Progenitor Cells Enhances PTH-Induced Bone Formation

Activation of bone anabolic pathways is a fruitful approach for treating severe osteoporosis, yet FDA-approved osteoanabolics, eg, parathyroid hormone (PTH), have limited efficacy. Improving their potency is a promising strategy for maximizing bone anabolic output. Nmp4 (Nuclear Matrix Protein 4) global knockout mice exhibit enhanced PTH-induced increases in trabecular bone but display no overt baseline skeletal phenotype. Nmp4 is expressed in all tissues; therefore, to determine which cell type is responsible for driving the beneficial effects of Nmp4 inhibition, we conditionally removed this gene from cells at distinct stages of osteogenic differentiation. Nmp4-floxed (Nmp4fl/fl ) mice were crossed with mice bearing one of three Cre drivers including (i) Prx1Cre+  to remove Nmp4 from mesenchymal stem/progenitor cells (MSPCs) in long bones; (ii) BglapCre+  targeting mature osteoblasts, and (iii) Dmp1Cre+  to disable Nmp4 in osteocytes. Virgin female Cre+  and Cre- mice (10 weeks of age) were sorted into cohorts by weight and genotype. Mice were administered daily injections of either human PTH 1-34 at 30 μg/kg or vehicle for 4 weeks or 7 weeks. Skeletal response was assessed using dual-energy X-ray absorptiometry, micro-computed tomography, bone histomorphometry, and serum analysis for remodeling markers. Nmp4fl/fl ;Prx1Cre+  mice virtually phenocopied the global Nmp4-/- skeleton in the femur, ie, a mild baseline phenotype but significantly enhanced PTH-induced increase in femur trabecular bone volume/total volume (BV/TV) compared with their Nmp4fl/fl ;Prx1Cre- controls. This was not observed in the spine, where Prrx1 is not expressed. Heightened response to PTH was coincident with enhanced bone formation. Conditional loss of Nmp4 from the mature osteoblasts (Nmp4fl/fl ;BglapCre+ ) failed to increase BV/TV or enhance PTH response. However, conditional disabling of Nmp4 in osteocytes (Nmp4fl/fl ;Dmp1Cre+ ) increased BV/TV without boosting response to hormone under our experimental regimen. We conclude that Nmp4-/- Prx1-expressing MSPCs drive the improved response to PTH therapy and that this gene has stage-specific effects on osteoanabolism. © 2022 American Society for Bone and Mineral Research (ASBMR).

Single-cell RNA sequencing reveals in vivo osteoimmunology interactions between the immune and skeletal systems

BACKGROUND

While osteoimmunology interactions between the immune and skeletal systems are known to play an important role in osteoblast development, differentiation and bone metabolism related disease like osteoporosis, such interactions in either bone microenvironment or peripheral circulation in vivo at the single-cell resolution have not yet been characterized.

METHODS

We explored the osteoimmunology communications between immune cells and osteoblastic lineage cells (OBCs) by performing CellphoneDB and CellChat analyses with single-cell RNA sequencing (scRNA-seq) data from human femoral head. We also explored the osteoimmunology effects of immune cells in peripheral circulation on skeletal phenotypes. We used a scRNA-seq dataset of peripheral blood monocytes (PBMs) to perform deconvolution analysis. Then weighted gene co-expression network analysis (WGCNA) was used to identify monocyte subtype-specific subnetworks. We next used cell-specific network (CSN) and the least absolute shrinkage and selection operator (LASSO) to analyze the correlation of a gene subnetwork identified by WGCNA with bone mineral density (BMD).

RESULTS

We constructed immune cell and OBC communication networks and further identified L-R genes, such as JAG1 and NOTCH1/2, with ossification related functions. We also found a Mono4 related subnetwork that may relate to BMD variation in both older males and postmenopausal female subjects.

CONCLUSIONS

This is the first study to identify numerous ligand-receptor pairs that likely mediate signals between immune cells and osteoblastic lineage cells. This establishes a foundation to reveal advanced and in-depth osteoimmunology interactions to better understand the relationship between local bone microenvironment and immune cells in peripheral blood and the impact on bone phenotypes.

Imputation of Human Primary Osteoblast Single Cell RNA-Seq Data Identified Three Novel Osteoblastic Subtypes

BACKGROUND

Recently, single-cell RNA sequencing (scRNA-seq) technology was increasingly used to study transcriptomics at a single-cell resolution, scRNA-seq analysis was complicated by the "dropout", where the data only captures a small fraction of the transcriptome. This phenomenon can lead to the fact that the actual expressed transcript may not be detected. We previously performed osteoblast subtypes classification and dissection on freshly isolated human osteoblasts.

MATERIALS AND METHODS

Here, we used the scImpute method to impute the missing values of dropout genes from a scRNA-seq dataset generated on freshly isolated human osteoblasts.

RESULTS

Based on the imputed gene expression patterns, we discovered three new osteoblast subtypes. Specifically, these newfound osteoblast subtypes are osteoblast progenitors, and two undetermined osteoblasts. Osteoblast progenitors showed significantly high expression of proliferation related genes (FOS, JUN, JUNB and JUND). Analysis of each subtype showed that in addition to bone formation, these undetermined osteoblasts may involve osteoclast and adipocyte differentiation and have the potential function of regulate immune activation.

CONCLUSIONS

Our findings provided a new perspective for studying the osteoblast heterogeneity and potential biological functions of these freshly isolated human osteoblasts at the single-cell level, which provides further insight into osteoblasts subtypes under various (pathological) physiological conditions.

Multi-omics molecular biomarkers and database of osteoarthritis

Osteoarthritis (OA) is the most common form of arthritis in the adult population and is a leading cause of disability. OA-related genetic loci may play an important role in clinical diagnosis and disease progression. With the rapid development of diverse technologies and omics methods, many OA-related public data sets have been accumulated. Here, we retrieved a diverse set of omics experimental results from 159 publications, including genome-wide association study, differentially expressed genes and differential methylation regions, and 2405 classified OA-related gene markers. Meanwhile, based on recent single-cell RNA-seq data from different joints, 5459 cell-type gene markers of joints were collected. The information has been integrated into an online database named OAomics and molecular biomarkers (OAOB). The database (http://ibi.zju.edu.cn/oaobdb/) provides a web server for OA marker genes, omics features and so on. To our knowledge, this is the first database of molecular biomarkers for OA.

Single-cell transcriptional profiling reveals heterogeneity and developmental trajectories of Ewing sarcoma

PURPOSE

Ewing sarcoma (EwS) is an aggressive malignant neoplasm composed of small round cells. The heterogeneity and developmental trajectories of EwS are uncertain.

METHODS

Single-cell RNA sequencing was performed on 4 EwS tumor tissue samples, and 3 transcriptional atlases were generated. K-nearest neighbor algorithm was used to predict the origin of tumor cells at single-cell resolution. Monocle2 package was used to perform pseudotime trajectory analysis in tumor cells. Differentially expressed genes were compared against those in all other clusters via the FindMarkers function, and then they were subjected to GO analysis using clusterProfiler package.

RESULTS

Combined with the results of k-nearest neighbor algorithm and pseudotime trajectory analysis in tumor cells, we thought meningeal EwS originated from neural crest cells during epithelial to mesenchymal transition and simulated the process of neural crest cell lineage differentiation. But for perirenal EwS and spinal EwS, we hypothesized that after the neural crest cell lineage mutated into them, the tumor cells did not maintain the differentiation trajectory of neural crest cell lineage, and the development trajectory of tumor cells became chaotic. GO analysis results showed that interferon signaling pathway-related biological processes play an essential role in the tumorigenesis and tumor progression process of EwS, and among these biological processes genes, JAK1 gene up-regulated most significantly and highly expressed in all tumor cells. Ruxolitinib was used to explore the function of JAK1. Targeting JAK1 can promote apoptosis of EwS tumor cells, inhibit the migration and invasion of EwS tumor cells, and inhibit cell proliferation by inducing cell cycle S phase arrest.

CONCLUSION

EwS was derived from neural crest cell lineage with variable developmental timing of oncogenic conversion, and the JAK1 might be a candidate for therapeutic targets of EwS.

Insight into the potential pathogenesis of human osteoarthritis via single-cell RNA sequencing data on osteoblasts

Osteoarthritis (OA) is the most common degenerative joint disease caused by osteoblastic lineage cells. However, a comprehensive molecular program for osteoblasts in human OA remains underdeveloped. The single-cell gene expression of osteoblasts and microRNA array data were from human. After processing the single-cell RNA sequencing (scRNA-seq) data, it was subjected to principal component analysis (PCA) and T-Stochastic neighbor embedding analysis (TSNE). Differential expression analysis was aimed to find marker genes. Gene-ontology (GO) enrichment, Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis and Gene set enrichment analysis (GSEA) were applied to characterize the molecular function of osteoblasts with marker genes. Protein-protein interaction (PPI) networks and core module were established for marker genes by using the STRING database and Cytoscape software. All nodes in the core module were considered to be hub genes. Subsequently, we predicted the potential miRNA of hub genes through the miRWalk, miRDB and TargetScan database and experimentally verified the miRNA by GSE105027. Finally, miRNA-mRNA regulatory network was constructed using the Cytoscape software. We characterized the single-cell expression profiling of 4387 osteoblasts from normal and OA sample. The proportion of osteoblasts subpopulations changed dramatically in the OA, with 70.42% of the pre-osteoblasts. 117 marker genes were included and the results of GO analysis show that up-regulated marker genes enriched in collagen-containing extracellular matrix were highly expressed in the pre-osteoblasts cluster. Both KEGG and GSEA analyses results indicated that IL-17 and NOD-like receptor signaling pathways were enriched in down-regulated marker genes. We visualize the weight of marker genes and constructed the core module in PPI network. In potential mRNA-miRNA regulatory network, hsa-miR-449a and hsa-miR-218-5p may be involved in the development of OA. Our study found that alterations in osteoblasts state and cellular molecular function in the subchondral bone region may be involved in the pathogenesis of osteoarthritis.

A systematic dissection of human primary osteoblasts in vivo at single-cell resolution

Human osteoblasts are multifunctional bone cells, which play essential roles in bone formation, angiogenesis regulation, as well as maintenance of hematopoiesis. However, the categorization of primary osteoblast subtypes in vivo in humans has not yet been achieved. Here, we used single-cell RNA sequencing (scRNA-seq) to perform a systematic cellular taxonomy dissection of freshly isolated human osteoblasts from one 31-year-old male with osteoarthritis and osteopenia after hip replacement. Based on the gene expression patterns and cell lineage reconstruction, we identified three distinct cell clusters including preosteoblasts, mature osteoblasts, and an undetermined rare osteoblast subpopulation. This novel subtype was found to be the major source of the nuclear receptor subfamily 4 group A member 1 and 2 (NR4A1 and NR4A2) in primary osteoblasts, and the expression of NR4A1 was confirmed by immunofluorescence staining on mouse osteoblasts in vivo. Trajectory inference analysis suggested that the undetermined cluster, together with the preosteoblasts, are involved in the regulation of osteoblastogenesis and also give rise to mature osteoblasts. Investigation of the biological processes and signaling pathways enriched in each subpopulation revealed that in addition to bone formation, preosteoblasts and undetermined osteoblasts may also regulate both angiogenesis and hemopoiesis. Finally, we demonstrated that there are systematic differences between the transcriptional profiles of human and mouse osteoblasts, highlighting the necessity for studying bone physiological processes in humans rather than solely relying on mouse models. Our findings provide novel insights into the cellular heterogeneity and potential biological functions of human primary osteoblasts at the single-cell level.