Novel Lineage-Tracing System to Identify Site-Specific Ectopic Bone Precursor Cells

Vegfc-expressing cells form heterotopic bone after musculoskeletal injury

Heterotopic ossification (HO) is a challenging condition that occurs after musculoskeletal injury and is characterized by the formation of bone in non-skeletal tissues. While the effect of HO on blood vessels is well established, little is known about its impact on lymphatic vessels. Here, we use a mouse model of traumatic HO to investigate the relationship between HO and lymphatic vessels. We show that injury triggers lymphangiogenesis at the injury site, which is associated with elevated vascular endothelial growth factor C (VEGF-C) levels. Through single-cell transcriptomic analyses, we identify mesenchymal progenitor cells and tenocytes as sources of Vegfc. We demonstrate by lineage tracing that Vegfc-expressing cells undergo osteochondral differentiation and contribute to the formation of HO. Last, we show that Vegfc haploinsufficiency results in a nearly 50% reduction in lymphangiogenesis and HO formation. These findings shed light on the complex mechanisms underlying HO formation and its impact on lymphatic vessels.

Cutaneous burn injury represents a major risk factor for the development of traumatic ectopic bone formation following blast-related extremity injury

Blast-related traumatic heterotopic ossification (tHO) impacts clinical outcomes in combat-injured patients, leading to delayed wound healing, inflammatory complications, and reduced quality of life. Blast injured patients often have significant burns. This study investigated whether a partial thickness thermal burn injury exacerbates blast-related tHO in a clinically relevant polytrauma animal model. Adult male Sprague Dawley rats were subjected to an established model involving a whole-body blast overpressure exposure (BOP), complex extremity trauma followed by hind limb amputation (CET) followed by the addition of a 10 % total body surface area (TBSA) second degree thermal burn (BU). Micro-CT scans on post-operative day 56 showed a significant increase in HO volume in the CET + BU as compared to the CET alone injury group (p < .0001; 22.83 ± 3.41 mm3 vs 4.84 ± 5.77 mm3). Additionally, CET + BU concomitant with BOP significantly increased HO (p < .0001; 34.95 ± 7.71 mm3) as compared to CET + BU alone, confirming BOP has a further synergistic effect. No HO was detectable in rats in the absence of CET. Serum analysis revealed similar significant elevated (p < .0001) levels of pro-inflammatory markers (Cxcl1 and Il6) at 6 h post-injury (hpi) in the CET + BU and BOP + CET + BU injury groups as compared to naïve baseline values. Real-time qPCR demonstrated similar levels of chondrogenic and osteogenic gene expression in muscle tissue at the site of injury at 168 hpi in both the CET + BU and BOP+CET + BU injury groups. These results support the hypothesis that a 10 % TBSA thermal burn markedly enhances tHO following acute musculoskeletal extremity injury in the presence and absence of blast overpressure. Furthermore, the influence of BOP on tHO cannot be accounted for either in regards to systemic inflammation induced from remote injury or inflammatory-osteo-chondrogenic expression changes local to the musculoskeletal trauma, suggesting that another mechanism beyond BOP and BU synergistic effects are at play. Therefore, these findings warrant future investigations to explore other mechanisms by which blast and burn influence tHO, and testing prophylactic measures to mitigate the local and systemic inflammatory effects of these injuries on development of HO.

The HIF-1α/PLOD2 axis integrates extracellular matrix organization and cell metabolism leading to aberrant musculoskeletal repair

While hypoxic signaling has been shown to play a role in many cellular processes, its role in metabolism-linked extracellular matrix (ECM) organization and downstream processes of cell fate after musculoskeletal injury remains to be determined. Heterotopic ossification (HO) is a debilitating condition where abnormal bone formation occurs within extra-skeletal tissues. Hypoxia and hypoxia-inducible factor 1α (HIF-1α) activation have been shown to promote HO. However, the underlying molecular mechanisms by which the HIF-1α pathway in mesenchymal progenitor cells (MPCs) contributes to pathologic bone formation remain to be elucidated. Here, we used a proven mouse injury-induced HO model to investigate the role of HIF-1α on aberrant cell fate. Using single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics analyses of the HO site, we found that collagen ECM organization is the most highly up-regulated biological process in MPCs. Zeugopod mesenchymal cell-specific deletion of Hif1α (Hoxa11-CreERT2; Hif1afl/fl) significantly mitigated HO in vivo. ScRNA-seq analysis of these Hoxa11-CreERT2; Hif1afl/fl mice identified the PLOD2/LOX pathway for collagen cross-linking as downstream of the HIF-1α regulation of HO. Importantly, our scRNA-seq data and mechanistic studies further uncovered that glucose metabolism in MPCs is most highly impacted by HIF-1α deletion. From a translational aspect, a pan-LOX inhibitor significantly decreased HO. A newly screened compound revealed that the inhibition of PLOD2 activity in MPCs significantly decreased osteogenic differentiation and glycolytic metabolism. This suggests that the HIF-1α/PLOD2/LOX axis linked to metabolism regulates HO-forming MPC fate. These results suggest that the HIF-1α/PLOD2/LOX pathway represents a promising strategy to mitigate HO formation.

Neutrophil and NETosis Modulation in Traumatic Heterotopic Ossification

OBJECTIVE

To characterize the role of neutrophil extracellular traps (NETs) in heterotopic ossification (HO) formation and progression and to use mechanical and pharmacological methods to decrease NETosis and mitigate HO formation.

BACKGROUND

Traumatic HO is the aberrant osteochondral differentiation of mesenchymal progenitor cells after traumatic injury, burns, or surgery. While the innate immune response has been shown to be necessary for HO formation, the specific immune cell phenotype and function remain unknown. Neutrophils, one of the earliest immune cells to respond after HO-inducing injuries, can extrude DNA, forming highly inflammatory NETs. We hypothesized that neutrophils and NETs would be diagnostic biomarkers and therapeutic targets for the detection and mitigation of HO.

METHODS

C57BL6J mice underwent burn/tenotomy (a well-established mouse model of HO) or a non-HO-forming sham injury. These mice were either (1) ambulated ad libitum, (2) ambulated ad libitum with daily intraperitoneal hydroxychloroquine, ODN-2088 (both known to affect NETosis pathways), or control injections, or (3) had the injured hind limb immobilized. Single-cell analysis was performed to analyze neutrophils, NETosis, and downstream signaling after the HO-forming injury. Immunofluorescence microscopy was used to visualize NETosis at the HO site and neutrophils were identified using flow cytometry. Serum and cell lysates from HO sites were analyzed using enzyme-linked immunosorbent assay for myeloperoxidase-DNA and ELA2-DNA complexes to identify NETosis. Micro-computerized tomography was performed on all groups to analyze the HO volume.

RESULTS

Molecular and transcriptional analyses revealed the presence of NETs within the HO injury site, which peaked in the early phases after injury. These NETs were highly restricted to the HO site, with gene signatures derived from both in vitro NET induction and clinical neutrophil characterizations showing a high degree of NET "priming" at the site of injury, but not in neutrophils in the blood or bone marrow. Cell-cell communication analyses revealed that this localized NET formation coincided with high levels of toll-like receptor signaling specific to neutrophils at the injury site. Reducing the overall neutrophil abundance within the injury site, either pharmacologically through treatment with hydroxychloroquine, the toll-like receptor 9 inhibitor OPN-2088, or mechanical treatment with limb offloading, results in the mitigation of HO formation.

CONCLUSIONS

These data provide a further understanding of the ability of neutrophils to form NETs at the injury site, clarify the role of neutrophils in HO, and identify potential diagnostic and therapeutic targets for HO mitigation.

Itaconate-producing neutrophils regulate local and systemic inflammation following trauma

Modulation of the immune response to initiate and halt the inflammatory process occurs both at the site of injury as well as systemically. Due to the evolving role of cellular metabolism in regulating cell fate and function, tendon injuries that undergo normal and aberrant repair were evaluated by metabolic profiling to determine its impact on healing outcomes. Metabolomics revealed an increasing abundance of the immunomodulatory metabolite itaconate within the injury site. Subsequent single-cell RNA-Seq and molecular and metabolomic validation identified a highly mature neutrophil subtype, not macrophages, as the primary producers of itaconate following trauma. These mature itaconate-producing neutrophils were highly inflammatory, producing cytokines that promote local injury fibrosis before cycling back to the bone marrow. In the bone marrow, itaconate was shown to alter hematopoiesis, skewing progenitor cells down myeloid lineages, thereby regulating systemic inflammation. Therapeutically, exogenous itaconate was found to reduce injury-site inflammation, promoting tenogenic differentiation and impairing aberrant vascularization with disease-ameliorating effects. These results present an intriguing role for cycling neutrophils as a sensor of inflammation induced by injury - potentially regulating immune cell production in the bone marrow through delivery of endogenously produced itaconate - and demonstrate a therapeutic potential for exogenous itaconate following tendon injury.

Genetic models for lineage tracing in musculoskeletal development, injury, and healing

Musculoskeletal development and later post-natal homeostasis are highly dynamic processes, marked by rapid structural and functional changes across very short periods of time. Adult anatomy and physiology are derived from pre-existing cellular and biochemical states. Consequently, these early developmental states guide and predict the future of the system as a whole. Tools have been developed to mark, trace, and follow specific cells and their progeny either from one developmental state to the next or between circumstances of health and disease. There are now many such technologies alongside a library of molecular markers which may be utilized in conjunction to allow for precise development of unique cell 'lineages'. In this review, we first describe the development of the musculoskeletal system beginning as an embryonic germ layer and at each of the key developmental stages that follow. We then discuss these structures in the context of adult tissues during homeostasis, injury, and repair. Special focus is given in each of these sections to the key genes involved which may serve as markers of lineage or later in post-natal tissues. We then finish with a technical assessment of lineage tracing and the techniques and technologies currently used to mark cells, tissues, and structures within the musculoskeletal system.

Tppp3+ synovial/tendon sheath progenitor cells contribute to heterotopic bone after trauma

Heterotopic ossification (HO) is a pathological process resulting in aberrant bone formation and often involves synovial lined tissues. During this process, mesenchymal progenitor cells undergo endochondral ossification. Nonetheless, the specific cell phenotypes and mechanisms driving this process are not well understood, in part due to the high degree of heterogeneity of the progenitor cells involved. Here, using a combination of lineage tracing and single-cell RNA sequencing (scRNA-seq), we investigated the extent to which synovial/tendon sheath progenitor cells contribute to heterotopic bone formation. For this purpose, Tppp3 (tubulin polymerization-promoting protein family member 3)-inducible reporter mice were used in combination with either Scx (Scleraxis) or Pdgfra (platelet derived growth factor receptor alpha) reporter mice. Both tendon injury- and arthroplasty-induced mouse experimental HO models were utilized. ScRNA-seq of tendon-associated traumatic HO suggested that Tppp3 is an early progenitor cell marker for either tendon or osteochondral cells. Upon HO induction, Tppp3 reporter⁺ cells expanded in number and partially contributed to cartilage and bone formation in either tendon- or joint-associated HO. In double reporter animals, both Pdgfra⁺Tppp3⁺ and Pdgfra⁺Tppp3^- progenitor cells gave rise to HO-associated cartilage. Finally, analysis of human samples showed a substantial population of TPPP3-expressing cells overlapping with osteogenic markers in areas of heterotopic bone. Overall, these data demonstrate that synovial/tendon sheath progenitor cells undergo aberrant osteochondral differentiation and contribute to HO after trauma.

Discoidin domain receptor 2 regulates aberrant mesenchymal lineage cell fate and matrix organization

Extracellular matrix (ECM) interactions regulate both the cell transcriptome and proteome, thereby determining cell fate. Traumatic heterotopic ossification (HO) is a disorder characterized by aberrant mesenchymal lineage (MLin) cell differentiation, forming bone within soft tissues of the musculoskeletal system following traumatic injury. Recent work has shown that HO is influenced by ECM-MLin cell receptor signaling, but how ECM binding affects cellular outcomes remains unclear. Using time course transcriptomic and proteomic analyses, we identified discoidin domain receptor 2 (DDR2), a cell surface receptor for fibrillar collagen, as a key MLin cell regulator in HO formation. Inhibition of DDR2 signaling, through either constitutive or conditional Ddr2 deletion or pharmaceutical inhibition, reduced HO formation in mice. Mechanistically, DDR2 perturbation alters focal adhesion orientation and subsequent matrix organization, modulating Focal Adhesion Kinase (FAK) and Yes1 Associated Transcriptional Regulator and WW Domain Containing Transcription Regulator 1 (YAP/TAZ)-mediated MLin cell signaling. Hence, ECM-DDR2 interactions are critical in driving HO and could serve as a previously unknown therapeutic target for treating this disease process.

Single-cell RNA sequencing unravels heterogeneity of skeletal progenitors and cell–cell interactions underlying the bone repair process

Introduction

Activation of skeletal progenitors upon tissue injury and the subsequent cell fate specification are tightly coordinated in the bone repair process. Although known osteoimmunological signaling networks play important roles in the microenvironment of the bone defect sites, the molecular mechanism underlying the bone repair process has not been fully understood.

Methods

To better understand the behavior of the skeletal progenitors and the heterogeneity of the cells during bone repair at the microenvironmental level, we performed a combinatorial analysis consisting of lineage tracing for skeletal progenitors using the Sox9-CreERT2;R26R^tdTomato mouse line followed by single-cell RNA sequencing (scRNA-seq) analysis using a mouse model of calvarial bone repair. To identify a therapeutic target for bone regeneration, further computational analysis was performed focusing on the identification of the cell–cell interactions, followed by pharmacological assessments with a critical-size calvarial bone defect mouse model.

Results

Lineage tracing analysis showed that skeletal progenitors marked by Sox9 were activated upon bone injury and contributed to bone repair by differentiating into osteoblasts. The scRNA-seq analysis characterized heterogeneous cell populations at the bone defect sites; the computational analysis predicted a bifurcated lineage from skeletal progenitors toward osteogenic and adipogenic lineages. Chemokine C–C motif ligand 9 (Ccl9) was identified as a signaling molecule that regulates bone regeneration in the mouse model, possibly through the regulation of adipogenic differentiation at the bone defect site.

Conclusion

Multipotential skeletal progenitors and the direction of the cell differentiation were characterized at single cell resolution in a mouse bone repair model. The Ccl9 signaling pathway may be a key factor directing osteogenesis from the progenitors in the model and may be a therapeutic target for bone regeneration.

•

Sox9-positive skeletal progenitors contributed to the calvaria bone repair process.

•

scRNA-seq analysis revealed a heterogeneous cell population at bone defect sites.

•

Skeletal progenitors had a bifurcated lineages of osteogenesis and adipogenesis.

•

Ccl9 was identified as an important signaling molecule regulating bone regeneration.

Single-cell RNA sequencing unravels heterogeneity of skeletal progenitors and cell-cell interactions underlying the bone repair process

Introduction

Activation of skeletal progenitors upon tissue injury and the subsequent cell fate specification are tightly coordinated in the bone repair process. Although known osteoimmunological signaling networks play important roles in the microenvironment of the bone defect sites, the molecular mechanism underlying the bone repair process has not been fully understood.

Methods

To better understand the behavior of the skeletal progenitors and the heterogeneity of the cells during bone repair at the microenvironmental level, we performed a combinatorial analysis consisting of lineage tracing for skeletal progenitors using the Sox9-CreERT2;R26R ^tdTomato mouse line followed by single-cell RNA sequencing (scRNA-seq) analysis using a mouse model of calvarial bone repair. To identify a therapeutic target for bone regeneration, further computational analysis was performed focusing on the identification of the cell-cell interactions, followed by pharmacological assessments with a critical-size calvarial bone defect mouse model.

Results

Lineage tracing analysis showed that skeletal progenitors marked by Sox9 were activated upon bone injury and contributed to bone repair by differentiating into osteoblasts. The scRNA-seq analysis characterized heterogeneous cell populations at the bone defect sites; the computational analysis predicted a bifurcated lineage from skeletal progenitors toward osteogenic and adipogenic lineages. Chemokine C-C motif ligand 9 (Ccl9) was identified as a signaling molecule that regulates bone regeneration in the mouse model, possibly through the regulation of adipogenic differentiation at the bone defect site.

Conclusion

Multipotential skeletal progenitors and the direction of the cell differentiation were characterized at single cell resolution in a mouse bone repair model. The Ccl9 signaling pathway may be a key factor directing osteogenesis from the progenitors in the model and may be a therapeutic target for bone regeneration.

Single-cell RNA sequencing unravels heterogeneity of skeletal progenitors and cell-cell interactions underlying the bone repair process

Introduction

Activation of skeletal progenitors upon tissue injury and the subsequent cell fate specification are tightly coordinated in the bone repair process. Although known osteoimmunological signaling networks play important roles in the microenvironment of the bone defect sites, the molecular mechanism underlying the bone repair process has not been fully understood.

Methods

To better understand the behavior of the skeletal progenitors and the heterogeneity of the cells during bone repair at the microenvironmental level, we performed a combinatorial analysis consisting of lineage tracing for skeletal progenitors using the Sox9-CreERT2;R26R ^tdTomato mouse line followed by single-cell RNA sequencing (scRNA-seq) analysis using a mouse model of calvarial bone repair. To identify a therapeutic target for bone regeneration, further computational analysis was performed focusing on the identification of the cell-cell interactions, followed by pharmacological assessments with a critical-size calvarial bone defect mouse model.

Results

Lineage tracing analysis showed that skeletal progenitors marked by Sox9 were activated upon bone injury and contributed to bone repair by differentiating into osteoblasts. The scRNA-seq analysis characterized heterogeneous cell populations at the bone defect sites; the computational analysis predicted a bifurcated lineage from skeletal progenitors toward osteogenic and adipogenic lineages. Chemokine C-C motif ligand 9 (Ccl9) was identified as a signaling molecule that regulates bone regeneration in the mouse model, possibly through the regulation of adipogenic differentiation at the bone defect site.

Conclusion

Multipotential skeletal progenitors and the direction of the cell differentiation were characterized at single cell resolution in a mouse bone repair model. The Ccl9 signaling pathway may be a key factor directing osteogenesis from the progenitors in the model and may be a therapeutic target for bone regeneration.

Single-cell RNA sequencing unravels heterogeneity of skeletal progenitors and cell-cell interactions underlying the bone repair process

Introduction

Activation of skeletal progenitors upon tissue injury and the subsequent cell fate specification are tightly coordinated in the bone repair process. Although known osteoimmunological signaling networks play important roles in the microenvironment of the bone defect sites, the molecular mechanism underlying the bone repair process has not been fully understood.

Methods

To better understand the behavior of the skeletal progenitors and the heterogeneity of the cells during bone repair at the microenvironmental level, we performed a combinatorial analysis consisting of lineage tracing for skeletal progenitors using the Sox9-CreERT2;R26R ^tdTomato mouse line followed by single-cell RNA sequencing (scRNA-seq) analysis using a mouse model of calvarial bone repair. To identify a therapeutic target for bone regeneration, further computational analysis was performed focusing on the identification of the cell-cell interactions, followed by pharmacological assessments with a critical-size calvarial bone defect mouse model.

Results

Lineage tracing analysis showed that skeletal progenitors marked by Sox9 were activated upon bone injury and contributed to bone repair by differentiating into osteoblasts. The scRNA-seq analysis characterized heterogeneous cell populations at the bone defect sites; the computational analysis predicted a bifurcated lineage from skeletal progenitors toward osteogenic and adipogenic lineages. Chemokine C-C motif ligand 9 (Ccl9) was identified as a signaling molecule that regulates bone regeneration in the mouse model, possibly through the regulation of adipogenic differentiation at the bone defect site.

Conclusion

Multipotential skeletal progenitors and the direction of the cell differentiation were characterized at single cell resolution in a mouse bone repair model. The Ccl9 signaling pathway may be a key factor directing osteogenesis from the progenitors in the model and may be a therapeutic target for bone regeneration.

Single-cell RNA sequencing unravels heterogeneity of skeletal progenitors and cell-cell interactions underlying the bone repair process

Introduction

Activation of skeletal progenitors upon tissue injury and the subsequent cell fate specification are tightly coordinated in the bone repair process. Although known osteoimmunological signaling networks play important roles in the microenvironment of the bone defect sites, the molecular mechanism underlying the bone repair process has not been fully understood.

Methods

To better understand the behavior of the skeletal progenitors and the heterogeneity of the cells during bone repair at the microenvironmental level, we performed a combinatorial analysis consisting of lineage tracing for skeletal progenitors using the Sox9-CreERT2;R26R ^tdTomato mouse line followed by single-cell RNA sequencing (scRNA-seq) analysis using a mouse model of calvarial bone repair. To identify a therapeutic target for bone regeneration, further computational analysis was performed focusing on the identification of the cell-cell interactions, followed by pharmacological assessments with a critical-size calvarial bone defect mouse model.

Results

Lineage tracing analysis showed that skeletal progenitors marked by Sox9 were activated upon bone injury and contributed to bone repair by differentiating into osteoblasts. The scRNA-seq analysis characterized heterogeneous cell populations at the bone defect sites; the computational analysis predicted a bifurcated lineage from skeletal progenitors toward osteogenic and adipogenic lineages. Chemokine C-C motif ligand 9 (Ccl9) was identified as a signaling molecule that regulates bone regeneration in the mouse model, possibly through the regulation of adipogenic differentiation at the bone defect site.

Conclusion

Multipotential skeletal progenitors and the direction of the cell differentiation were characterized at single cell resolution in a mouse bone repair model. The Ccl9 signaling pathway may be a key factor directing osteogenesis from the progenitors in the model and may be a therapeutic target for bone regeneration.

Single-cell RNA sequencing unravels heterogeneity of skeletal progenitors and cell-cell interactions underlying the bone repair process

Introduction

Activation of skeletal progenitors upon tissue injury and the subsequent cell fate specification are tightly coordinated in the bone repair process. Although known osteoimmunological signaling networks play important roles in the microenvironment of the bone defect sites, the molecular mechanism underlying the bone repair process has not been fully understood.

Methods

To better understand the behavior of the skeletal progenitors and the heterogeneity of the cells during bone repair at the microenvironmental level, we performed a combinatorial analysis consisting of lineage tracing for skeletal progenitors using the Sox9-CreERT2;R26R ^tdTomato mouse line followed by single-cell RNA sequencing (scRNA-seq) analysis using a mouse model of calvarial bone repair. To identify a therapeutic target for bone regeneration, further computational analysis was performed focusing on the identification of the cell-cell interactions, followed by pharmacological assessments with a critical-size calvarial bone defect mouse model.

Results

Lineage tracing analysis showed that skeletal progenitors marked by Sox9 were activated upon bone injury and contributed to bone repair by differentiating into osteoblasts. The scRNA-seq analysis characterized heterogeneous cell populations at the bone defect sites; the computational analysis predicted a bifurcated lineage from skeletal progenitors toward osteogenic and adipogenic lineages. Chemokine C-C motif ligand 9 (Ccl9) was identified as a signaling molecule that regulates bone regeneration in the mouse model, possibly through the regulation of adipogenic differentiation at the bone defect site.

Conclusion

Multipotential skeletal progenitors and the direction of the cell differentiation were characterized at single cell resolution in a mouse bone repair model. The Ccl9 signaling pathway may be a key factor directing osteogenesis from the progenitors in the model and may be a therapeutic target for bone regeneration.

CTR9 drives osteochondral lineage differentiation of human mesenchymal stem cells via epigenetic regulation of BMP-2 signaling

Cell fate determination of human mesenchymal stem/stromal cells (hMSCs) is precisely regulated by lineage-specific transcription factors and epigenetic enzymes. We found that CTR9, a key scaffold subunit of polymerase-associated factor complex (PAFc), selectively regulates hMSC differentiation to osteoblasts and chondrocytes, but not to adipocytes. An in vivo ectopic osteogenesis assay confirmed the essentiality of CTR9 in hMSC-derived bone formation. CTR9 counteracts the activity of Enhancer Of Zeste 2 (EZH2), the epigenetic enzyme that deposits H3K27me3, in hMSCs. Accordingly, CTR9 knockdown (KD) hMSCs gain H3K27me3 mark, and the osteogenic differentiation defects of CTR9 KD hMSCs can be partially rescued by treatment with EZH2 inhibitors. Transcriptome analyses identified bone morphology protein-2 (BMP-2) as a downstream effector of CTR9. BMP-2 secretion, membrane anchorage, and the BMP-SMAD pathway were impaired in CTR9 KD MSCs, and the effects were rescued by BMP-2 supplementation. This study uncovers an epigenetic mechanism engaging the CTR9-H3K27me3-BMP-2 axis to regulate the osteochondral lineage differentiation of hMSCs.

Functional Heterogeneity of Bone Marrow Mesenchymal Stem Cell Subpopulations in Physiology and Pathology

Bone marrow mesenchymal stem cells (BMSCs) are multi-potent cell populations and are capable of maintaining bone and body homeostasis. The stemness and potential therapeutic effect of BMSCs have been explored extensively in recent years. However, diverse cell surface antigens and complex gene expression of BMSCs have indicated that BMSCs represent heterogeneous populations, and the natural characteristics of BMSCs make it difficult to identify the specific subpopulations in pathological processes which are often obscured by bulk analysis of the total BMSCs. Meanwhile, the therapeutic effect of total BMSCs is often less effective partly due to their heterogeneity. Therefore, understanding the functional heterogeneity of the BMSC subpopulations under different physiological and pathological conditions could have major ramifications for global health. Here, we summarize the recent progress of functional heterogeneity of BMSC subpopulations in physiology and pathology. Targeting tissue-resident single BMSC subpopulation offers a potentially innovative therapeutic strategy and improves BMSC effectiveness in clinical application.

Macrophage-mediated PDGF Activation Correlates with Regenerative Outcomes Following Musculoskeletal Trauma

OBJECTIVE

Our objective was to identify macrophage subpopulations and gene signatures associated with regenerative or fibrotic healing across different musculoskeletal injury types.

BACKGROUND

Subpopulations of macrophages are hypothesized to fine tune the immune response after damage, promoting either normal regenerative, or aberrant fibrotic healing.

METHODS

Mouse single-cell RNA sequencing data before and after injury were assembled from models of musculoskeletal injury, including regenerative and fibrotic mouse volumetric muscle loss (VML), regenerative digit tip amputation (DTA), and fibrotic heterotopic ossification (HO). R packages Harmony , MacSpectrum and Seurat were used for data integration, analysis and visualizations.

RESULTS

There was substantial overlap between macrophages from the regenerative VML (2 mm injury) and regenerative bone (DTA) models, as well as a separate overlap between the fibrotic VML (3 mm injury) and fibrotic bone (HO) models. We identified 2 fibrotic-like (FL 1 and FL 2) along with 3 regenerative-like (RL 1, RL 2, and RL 3) subpopulations of macrophages, each of which was transcriptionally distinct. We found that regenerative and fibrotic conditions had similar compositions of pro-inflammatory and anti-inflammatory macrophages, suggesting that macrophage polarization state did not correlate with healing outcomes. Receptor/ligand analysis of macrophage-to-mesenchymal progenitor cell (MPC) crosstalk showed enhanced transforming growth factor beta (TGF-β) in fibrotic conditions and enhanced platelet derived growth factor (PDGF) signaling in regenerative conditions.

CONCLUSION

Characterization of macrophage subtypes could be used to predict fibrotic responses following injury and provide a therapeutic target to tune the healing microenvironment towards more regenerative conditions.

Single-cell mapping of regenerative and fibrotic healing responses after musculoskeletal injury

After injury, a cascade of events repairs the damaged tissue, including expansion and differentiation of the progenitor pool and redeposition of matrix. To guide future wound regeneration strategies, we compared single-cell sequencing of regenerative (third phalangeal element [P3]) and fibrotic (second phalangeal element [P2]) digit tip amputation (DTA) models as well as traumatic heterotopic ossification (HO; aberrant). Analyses point to a common initial response to injury, including expansion of progenitors, redeposition of matrix, and activation of transforming growth factor β (TGF-β) and WNT pathways. Surprisingly, fibrotic P2 DTA showed greater transcriptional similarity to HO than to regenerative P3 DTA, suggesting that gene expression more strongly correlates with healing outcome than with injury type or cell origin. Differential analysis and immunostaining revealed altered activation of inflammatory pathways, such as the complement pathway, in the progenitor cells. These data suggests that common pathways are activated in response to damage but are fine tuned within each injury. Modulating these pathways may shift the balance toward regenerative outcomes.

•

Regenerative and fibrotic injuries share common early response mechanisms

•

Transcriptomes correlate with healing outcome more than injury type or cell source

•

Matrix composition after injury-induced tissue repair is highly injury type dependent

•

Inflammatory cascades are activated in immune and mesenchymal cells

Single-Cell Analyses of Heterotopic Ossification: Characteristics of Injury-Related Senescent Fibroblasts

Introduction

Injury-related cellular senescence may be involved in heterotopic ossification, and no research has been performed about this before.

Methods and Results

The study utilized integrated single-cell RNA-sequencing (scRNA-seq) data from heterotopic ossification samples. The number of senescent cells increased from day 3 and reached the highest level at day 21. However, the expression level of Cyclin Dependent Kinase Inhibitor 2A (Cdkn2a) has no such tendency as the change of cell amount, indicating that the expression level of Cdkn2a may be different in different types of senescent cells or the same time of senescent cell at different time points. The expression level of SASPs (senescence associated secret phenotypes) was also different in different types of senescent cells or at different time points. The GO (gene ontology) analysis revealed that the senescent cells were significantly correlated with the ossification processes, like ECM organization, cell adhesion, ossification, cartilage development, etc. Trajectory analysis showed that injury-related senescent fibroblasts (day 7 and 21) and age-related senescent fibroblasts (day 0 and 42) were in different branches. GO analysis demonstrated that injury-related senescent fibroblasts were mainly related to ossification and ECM remodeling. The KEGG (Kyoto Encyclopedia of Genes and Genomes) results revealed that the ossification was significantly corrected with protein processing in PI3K-Akt signaling, MAPK signaling, focal adhesion, etc.

Conclusion

Consequently, we demonstrated that, unlike age-related senescence, the injury-related senescence demonstrated significantly different SASP phenotypes. The injury-related senescence of fibroblasts is associated with heterotopic ossification formation and may act through PI3K/Akt-induced SASPs.

A new mouse model of post-traumatic joint injury allows to identify the contribution of Gli1+ mesenchymal progenitors in arthrofibrosis and acquired heterotopic endochondral ossification

Complex injury and open reconstructive surgeries of the knee often lead to joint dysfunction that may alter the normal biomechanics of the joint. Two major complications that often arise are excessive deposition of fibrotic tissue and acquired heterotopic endochondral ossification. Knee arthrofibrosis is a fibrotic joint disorder where aberrant buildup of scar tissue and adhesions develop around the joint. Heterotopic ossification is ectopic bone formation around the periarticular tissues. Even though arthrofibrosis and heterotopic ossification pose an immense clinical problem, limited studies focus on their cellular and molecular mechanisms. Effective cell-targeted therapeutics are needed, but the cellular origin of both knee disorders remains elusive. Moreover, all the current animal models of knee arthrofibrosis and stiffness are developed in rats and rabbits, limiting genetic experiments that would allow us to explore the contribution of specific cellular targets to these knee pathologies. Here, we present a novel mouse model where surgically induced injury and hyperextension of the knee lead to excessive deposition of disorganized collagen in the meniscus, synovium, and joint capsule in addition to formation of extra-skeletal bone in muscle and soft tissues within the joint capsule. As a functional outcome, arthrofibrosis and acquired heterotopic endochondral ossification coupled with a significant increase in total joint stiffness were observed. By employing this injury model and genetic lineage tracing, we also demonstrate that Gli1+ mesenchymal progenitors proliferate after joint injury and contribute to the pool of fibrotic cells in the synovium and ectopic osteoblasts within the joint capsule. These findings demonstrate that Gli1+ cells are a major cellular contributor to knee arthrofibrosis and acquired heterotopic ossification that manifest after knee injury. Our data demonstrate that genetic manipulation of Gli1+ cells in mice may offer a platform for identification of novel therapeutic targets to prevent knee joint dysfunction after chronic injury.

Contemporary perspectives on heterotopic ossification

Heterotopic ossification (HO) is the formation of ectopic bone that is primarily genetically driven (fibrodysplasia ossificans progressiva [FOP]) or acquired in the setting of trauma (tHO). HO has undergone intense investigation, especially over the last 50 years, as awareness has increased around improving clinical technologies and incidence, such as with ongoing wartime conflicts. Current treatments for tHO and FOP remain prophylactic and include NSAIDs and glucocorticoids, respectively, whereas other proposed therapeutic modalities exhibit prohibitive risk profiles. Contemporary studies have elucidated mechanisms behind tHO and FOP and have described new distinct niches independent of inflammation that regulate ectopic bone formation. These investigations have propagated a paradigm shift in the approach to treatment and management of a historically difficult surgical problem, with ongoing clinical trials and promising new targets.

Pathophysiology and Emerging Molecular Therapeutic Targets in Heterotopic Ossification

The term heterotopic ossification (HO) describes bone formation in tissues where bone is normally not present. Musculoskeletal trauma induces signalling events that in turn trigger cells, probably of mesenchymal origin, to differentiate into bone. The aetiology of HO includes extremely rare but severe, generalised and fatal monogenic forms of the disease; and as a common complex disorder in response to musculoskeletal, neurological or burn trauma. The resulting bone forms through a combination of endochondral and intramembranous ossification, depending on the aetiology, initiating stimulus and affected tissue. Given the heterogeneity of the disease, many cell types and biological pathways have been studied in efforts to find effective therapeutic strategies for the disorder. Cells of mesenchymal, haematopoietic and neuroectodermal lineages have all been implicated in the pathogenesis of HO, and the emerging dominant signalling pathways are thought to occur through the bone morphogenetic proteins (BMP), mammalian target of rapamycin (mTOR), and retinoic acid receptor pathways. Increased understanding of these disease mechanisms has resulted in the emergence of several novel investigational therapeutic avenues, including palovarotene and other retinoic acid receptor agonists and activin A inhibitors that target both canonical and non-canonical signalling downstream of the BMP type 1 receptor. In this article we aim to illustrate the key cellular and molecular mechanisms involved in the pathogenesis of HO and outline recent advances in emerging molecular therapies to treat and prevent HO that have had early success in the monogenic disease and are currently being explored in the common complex forms of HO.

Single-cell RNA sequencing unravels heterogeneity of skeletal progenitors and cell-cell interactions underlying the bone repair process

Introduction

Activation of skeletal progenitors upon tissue injury and the subsequent cell fate specification are tightly coordinated in the bone repair process. Although known osteoimmunological signaling networks play important roles in the microenvironment of the bone defect sites, the molecular mechanism underlying the bone repair process has not been fully understood.

Methods

To better understand the behavior of the skeletal progenitors and the heterogeneity of the cells during bone repair at the microenvironmental level, we performed a combinatorial analysis consisting of lineage tracing for skeletal progenitors using the Sox9-CreERT2;R26R ^tdTomato mouse line followed by single-cell RNA sequencing (scRNA-seq) analysis using a mouse model of calvarial bone repair. To identify a therapeutic target for bone regeneration, further computational analysis was performed focusing on the identification of the cell-cell interactions, followed by pharmacological assessments with a critical-size calvarial bone defect mouse model.

Results

Lineage tracing analysis showed that skeletal progenitors marked by Sox9 were activated upon bone injury and contributed to bone repair by differentiating into osteoblasts. The scRNA-seq analysis characterized heterogeneous cell populations at the bone defect sites; the computational analysis predicted a bifurcated lineage from skeletal progenitors toward osteogenic and adipogenic lineages. Chemokine C-C motif ligand 9 (Ccl9) was identified as a signaling molecule that regulates bone regeneration in the mouse model, possibly through the regulation of adipogenic differentiation at the bone defect site.

Conclusion

Multipotential skeletal progenitors and the direction of the cell differentiation were characterized at single cell resolution in a mouse bone repair model. The Ccl9 signaling pathway may be a key factor directing osteogenesis from the progenitors in the model and may be a therapeutic target for bone regeneration.

Cell autonomous TGFβ signaling is essential for stem/progenitor cell recruitment into degenerative tendons

Understanding cell recruitment in damaged tendons is critical for improvements in regenerative therapy. We recently reported that targeted disruption of transforming growth factor beta (TGFβ) type II receptor in the tendon cell lineage (Tgfbr2^ScxCre) resulted in resident tenocyte dedifferentiation and tendon deterioration in postnatal stages. Here we extend the analysis and identify direct recruitment of stem/progenitor cells into the degenerative mutant tendons. Cre-mediated lineage tracing indicates that these cells are not derived from tendon-ensheathing tissues or from a Scleraxis-expressing lineage, and they turned on tendon markers only upon entering the mutant tendons. Through immunohistochemistry and inducible gene deletion, we further find that the recruited cells originated from a Sox9-expressing lineage and their recruitment was dependent on cell autonomous TGFβ signaling. The cells identified in this study thus differ from previous reports of cell recruitment into injured tendons and suggest a critical role for TGFβ signaling in cell recruitment, providing insights that may support improvements in tendon repair.

•

Targeted deletion of TGFβ signaling led to degenerative changes in mouse tendons

•

Stem/progenitor cells were recruited into the degenerative mutant tendons

•

The recruited cells are different from the ones so far reported in tendon injury

•

Recruitment was dependent on cell autonomous TGFβ signaling in the recruited cells

High Frequency Spectral Ultrasound Imaging Detects Early Heterotopic Ossification in Rodents

Heterotopic ossification (HO) is a devastating condition in which ectopic bone forms inappropriately in soft tissues following traumatic injuries and orthopedic surgeries as a result of aberrant mesenchymal progenitor cell (MPC) differentiation. HO leads to chronic pain, decreased range of motion, and an overall decrease in quality of life. While several treatments have shown promise in animal models, all must be given during early stages of formation. Methods for early determination of whether and where endochondral ossification/soft tissue mineralization (HO anlagen) develop are lacking. At-risk patients are not identified sufficiently early in the process of MPC differentiation and soft tissue endochondral ossification for potential treatments to be effective. Hence, a critical need exists to develop technologies capable of detecting HO anlagen soon after trauma, when treatments are most effective. In this study, we investigate high frequency spectral ultrasound imaging (SUSI) as a noninvasive strategy to identify HO anlagen at early time points after injury. We show that by determining quantitative parameters based on tissue organization and structure, SUSI identifies HO anlagen as early as 1-week postinjury in a mouse model of burn/tenotomy and 3 days postinjury in a rat model of blast/amputation. We analyze single cell RNA sequencing profiles of the MPCs responsible for HO formation and show that the early tissue changes detected by SUSI match chondrogenic and osteogenic gene expression in this population. SUSI identifies sites of soft tissue endochondral ossification at early stages of HO formation so that effective intervention can be targeted when and where it is needed following trauma-induced injury. Furthermore, we characterize the chondrogenic to osteogenic transition that occurs in the MPCs during HO formation and correlate gene expression to SUSI detection of the HO anlagen.