Single-cell RNA sequencing of human femoral head in vivo

In Vitro Cell Surface Marker Expression on Mesenchymal Stem Cell Cultures does not Reflect Their Ex Vivo Phenotype

Cell surface marker expression is one of the criteria for defining human mesenchymal stem or stromal cells (MSC) in vitro. However, it is unclear if expression of markers including CD73 and CD90 reflects the in vivo origin of cultured cells. We evaluated expression of 15 putative MSC markers in primary cultured cells from periosteum and cartilage to determine whether expression of these markers reflects either the differentiation state of cultured cells or the self-renewal of in vivo populations. Cultured cells had universal and consistent expression of various putative stem cell markers including > 95% expression CD73, CD90 and PDPN in both periosteal and cartilage cultures. Altering the culture surface with extracellular matrix coatings had minimal effect on cell surface marker expression. Osteogenic differentiation led to loss of CD106 and CD146 expression, however CD73 and CD90 were retained in > 90% of cells. We sorted freshly isolated periosteal populations capable of CFU-F formation on the basis of CD90 expression in combination with CD34, CD73 and CD26. All primary cultures universally expressed CD73 and CD90 and lacked CD34, irrespective of the expression of these markers ex vivo indicating phenotypic convergence in vitro. We conclude that markers including CD73 and CD90 are acquired in vitro in most 'mesenchymal' cells capable of expansion. Overall, we demonstrate that in vitro expression of many cell surface markers in plastic-adherent cultures is unrelated to their expression prior to culture.

Carboxypeptidase M modulates BMSCs osteogenesis-adipogenesis via the MAPK/ERK pathway: An integrated single-cell and bulk transcriptomic study

The pathogenesis of osteoporosis (OP) is closely associated with the disrupted balance between osteogenesis and adipogenesis in bone marrow-derived mesenchymal stem cells (BMSCs). We analyzed published single-cell RNA sequencing (scRNA-seq) data to dissect the transcriptomic profiles of bone marrow-derived cells in OP, reviewing 56 377 cells across eight scRNA-seq datasets from femoral heads (osteoporosis or osteopenia n = 5, osteoarthritis n = 3). Seventeen genes, including carboxypeptidase M (CPM), were identified as key osteogenesis-adipogenesis regulators through comprehensive gene set enrichment, differential expression, regulon activity, and pseudotime analyses. In vitro, CPM knockdown reduced osteogenesis and promoted adipogenesis in BMSCs, while adenovirus-mediated CPM overexpression had the reverse effects. In vivo, intraosseous injection of CPM-overexpressing BMSCs mitigated bone loss in ovariectomized mice. Integrated scRNA-seq and bulk RNA sequencing analyses provided insight into the MAPK/ERK pathway's role in the CPM-mediated regulation of BMSC osteogenesis and adipogenesis; specifically, CPM overexpression enhanced MAPK/ERK signaling and osteogenesis. In contrast, the ERK1/2 inhibitor binimetinib negated the effects of CPM overexpression. Overall, our findings identify CPM as a pivotal regulator of BMSC differentiation, which provides new clues for the mechanistic study of OP.

Long noncoding RNA Malat1 protects against osteoporosis and bone metastasis

MALAT1, one of the few highly conserved nuclear long noncoding RNAs (lncRNAs), is abundantly expressed in normal tissues. Previously, targeted inactivation and genetic rescue experiments identified MALAT1 as a suppressor of breast cancer lung metastasis. On the other hand, Malat1-knockout mice are viable and develop normally. On a quest to discover the fundamental roles of MALAT1 in physiological and pathological processes, we find that this lncRNA is downregulated during osteoclastogenesis in humans and mice. Remarkably, Malat1 deficiency in mice promotes osteoporosis and bone metastasis of melanoma and mammary tumor cells, which can be rescued by genetic add-back of Malat1. Mechanistically, Malat1 binds to Tead3 protein, a macrophage-osteoclast-specific Tead family member, blocking Tead3 from binding and activating Nfatc1, a master regulator of osteoclastogenesis, which results in the inhibition of Nfatc1-mediated gene transcription and osteoclast differentiation. Notably, single-cell transcriptome analysis of clinical bone samples reveals that reduced MALAT1 expression in pre-osteoclasts and osteoclasts is associated with osteoporosis and metastatic bone lesions. Altogether, these findings identify Malat1 as a lncRNA that protects against osteoporosis and bone metastasis.

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.

A roadmap for delivering a human musculoskeletal cell atlas

Advances in single-cell technologies have transformed the ability to identify the individual cell types present within tissues and organs. The musculoskeletal bionetwork, part of the wider Human Cell Atlas project, aims to create a detailed map of the healthy musculoskeletal system at a single-cell resolution throughout tissue development and across the human lifespan, with complementary generation of data from diseased tissues. Given the prevalence of musculoskeletal disorders, this detailed reference dataset will be critical to understanding normal musculoskeletal function in growth, homeostasis and ageing. The endeavour will also help to identify the cellular basis for disease and lay the foundations for novel therapeutic approaches to treating diseases of the joints, soft tissues and bone. Here, we present a Roadmap delineating the critical steps required to construct the first draft of a human musculoskeletal cell atlas. We describe the key challenges involved in mapping the extracellular matrix-rich, but cell-poor, tissues of the musculoskeletal system, outline early milestones that have been achieved and describe the vision and directions for a comprehensive musculoskeletal cell atlas. By embracing cutting-edge technologies, integrating diverse datasets and fostering international collaborations, this endeavour has the potential to drive transformative changes in musculoskeletal medicine.

Major vault protein (MVP) suppresses aging- and estrogen deficiency-related bone loss through Fas-mediated apoptosis in osteoclasts

Osteoclasts (OCs), derived from monocyte/macrophage lineage, are key orchestrators in bone remodeling. Targeting osteoclast apoptosis is a promising approach to cut down excessive osteoclast numbers, and thus slow down the rate of bone mass loss that inevitably occurs during aging. However, the therapeutic target of apoptosis in osteoclasts has not been fully studied. Our previous work generated Mvpf/fLyz2-Cre mice, conditionally depleting major vault protein (MVP) in monocyte lineage, and identified MVP as a bone protector for its negative role in osteoclastogenesis in vivo and in vitro. Here, we observed a notable decline of MVP in osteoclasts with aging in mice, encouraging us to further investigate the regulatory role of osteoclast MVP. Then, Mvpf/fLyz2-Cre mice were exploited in two osteoporosis contexts, aging and abrupt loss of estrogen, and we revealed that conditional knockout of MVP inhibited osteoclast apoptosis in vivo and in vitro. Moreover, we reported the interaction between MVP and death receptor Fas, and MVP-Fas signaling cascade was identified to positively regulate the apoptosis of osteoclasts, thus preventing osteoporosis. Collectively, our comprehensive discovery of MVP's regulatory role in osteoclasts provides new insight into osteoclast biology and therapeutic targets for osteoporosis.

Inflammation produced by senescent osteocytes mediates age-related bone loss

Purpose

The molecular mechanisms of age-related bone loss are unclear and without valid drugs yet. The aims of this study were to explore the molecular changes that occur in bone tissue during age-related bone loss, to further clarify the changes in function, and to predict potential therapeutic drugs.

Methods

We collected bone tissues from children, middle-aged individuals, and elderly people for protein sequencing and compared the three groups of proteins pairwise, and the differentially expressed proteins (DEPs) in each group were analyzed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). K-means cluster analysis was then used to screen out proteins that continuously increased/decreased with age. Canonical signaling pathways that were activated or inhibited in bone tissue along with increasing age were identified by Ingenuity Pathway Analysis (IPA). Prediction of potential drugs was performed using the Connectivity Map (CMap). Finally, DEPs from sequencing were verified by Western blot, and the drug treatment effect was verified by quantitative real-time PCR.

Results

The GO and KEGG analyses show that the DEPs were associated with inflammation and bone formation with aging, and the IPA analysis shows that pathways such as IL-8 signaling and acute-phase response signaling were activated, while glycolysis I and EIF2 signaling were inhibited. A total of nine potential drugs were predicted, with rapamycin ranking the highest. In cellular experiments, rapamycin reduced the senescence phenotype produced by the H₂O₂-stimulated osteocyte-like cell MLO-Y4.

Conclusion

With age, inflammatory pathways are activated in bone tissue, and signals that promote bone formation are inhibited. This study contributes to the understanding of the molecular changes that occur in bone tissue during age-related bone loss and provides evidence that rapamycin is a drug of potential clinical value for this disease. The therapeutic effects of the drug are to be further studied in animals.

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.

The role of immune cells in modulating chronic inflammation and osteonecrosis

Osteonecrosis occurs when, under continuous stimulation by adverse factors such as glucocorticoids or alcohol, the death of local bone and marrow cells leads to abnormal osteoimmune function. This creates a chronic inflammatory microenvironment, which interferes with bone regeneration and repair. In a variety of bone tissue diseases, innate immune cells and adaptive immune cells interact with bone cells, and their effects on bone metabolic homeostasis have attracted more and more attention, thus developing into a new discipline - osteoimmunology. Immune cells are the most important regulator of inflammation, and osteoimmune disorder may be an important cause of osteonecrosis. Elucidating the chronic inflammatory microenvironment regulated by abnormal osteoimmune may help develop potential treatments for osteonecrosis. This review summarizes the inflammatory regulation of bone immunity in osteonecrosis, explains the pathophysiological mechanism of osteonecrosis from the perspective of osteoimmunology, and provides new ideas for the treatment of osteonecrosis.

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.

Single-cell Transcriptomic Analysis Reveals the Cellular Heterogeneity of Mesenchymal Stem Cells

Ex vivo-expanded mesenchymal stem cells (MSCs) have been demonstrated to be a heterogeneous mixture of cells exhibiting varying proliferative, multipotential, and immunomodulatory capacities. However, the exact characteristics of MSCs remain largely unknown. By single-cell RNA sequencing of 61,296 MSCs derived from bone marrow and Wharton’s jelly, we revealed five distinct subpopulations. The developmental trajectory of these five MSC subpopulations was mapped, revealing a differentiation path from stem-like active proliferative cells (APCs) to multipotent progenitor cells, followed by branching into two paths: 1) unipotent preadipocytes or 2) bipotent prechondro-osteoblasts that were subsequently differentiated into unipotent prechondrocytes. The stem-like APCs, expressing the perivascular mesodermal progenitor markers CSPG4/MCAM/NES, uniquely exhibited strong proliferation and stemness signatures. Remarkably, the prechondrocyte subpopulation specifically expressed immunomodulatory genes and was able to suppress activated CD3⁺ T cell proliferation in vitro, supporting the role of this population in immunoregulation. In summary, our analysis mapped the heterogeneous subpopulations of MSCs and identified two subpopulations with potential functions in self-renewal and immunoregulation. Our findings advance the definition of MSCs by identifying the specific functions of their heterogeneous cellular composition, allowing for more specific and effective MSC application through the purification of their functional subpopulations.