A single cell transcriptional atlas of early synovial joint development

Synergistic effects of biological stimuli and flexion induce microcavities promote hypertrophy and inhibit chondrogenesis during in vitro culture of human mesenchymal stem cell aggregates

Interzone/cavitation are key steps in early stage joint formation that have not been successfully developed in vitro. Further, current models of endochondral ossification, an important step in early bone formation, lack key morphology morphological structures such as microcavities found during development in vivo. This is possibly due to the lack of appropriate strategies for incorporating chemical and mechanical stimuli that are thought to be involved in joint development. We designed a bioreactor system and investigated the synergic effect of chemical stimuli (chondrogenesis-inducing [CIM] and hypertrophy-inducing medium [HIM]) and mechanical stimuli (flexion) on the growth of human mesenchymal stem cells (hMSCs) based linear aggregates under different conditions over 4 weeks of perfusion culture. Computational studies were used to evaluate tissue stress qualitatively. After harvesting, both Safranin-O and hematoxylin & eosin (H&E) staining histology demonstrated microcavity structures and void structures in the region of higher stresses for tissue aggregates cultured only in HIM under flexion. In comparison to either HIM treatment or flexion only, increased glycosaminoglycan (GAG) content in the extracellular matrix (ECM) at this region indicates the morphological change resembles the early stage of joint cavitation; while decreased type II collagen (Col II), and increased type X collagen (Col X) and vascular endothelial growth factor (VEGF) with a clear boundary in the staining section indicates it resembles the early stage of ossification. Further, cell alignment analysis indicated that cells were mostly oriented toward the direction of flexion in high-stress region only in HIM under flexion, resembling cell morphology in both joint cavitation and hypertrophic cartilage in growth plate. Collectively, our results suggest that flexion and HIM inhibit chondrogenesis and promote hypertrophy and development of microcavities that resemble the early stage of joint cavitation and endochondral ossification. We believe the tissue model described in this work can be used to develop in vitro models of joint tissue for applications such as pathophysiology and drug discovery.

Time- and cell-specific activation of BMP signaling restrains chondrocyte hypertrophy

Stem cell therapies for degenerative cartilage disease are limited by an incomplete understanding of hyaline cartilage formation and maintenance. Human bone marrow stromal cells/skeletal stem cells (hBMSCs/SSCs) produce stable hyaline cartilage when attached to hyaluronic acid-coated fibrin microbeads (HyA-FMBs), yet the mechanism remains unclear. In vitro, hBMSC/SSC/HyA-FMB organoids exhibited reduced BMP signaling early in chondrogenic differentiation, followed by restoration of BMP signaling in chondrogenic IGFBP5 + /MGP + cells. Subsequently, human-induced pluripotent stem cell (hiPSC)-derived sclerotome cells were established (BMP inhibition) and then treated with transforming growth factor β (TGF-β) -/+ BMP2 and growth differentiation factor 5 (GDF5) (BMP signaling activation). TGF-β alone elicited a weak chondrogenic response, but TGF-β/BMP2/GDF5 led to delamination of SOX9 + aggregates (chondrospheroids) with high expression of COL2A1, ACAN, and PRG4 and minimal expression of COL10A1 and ALP in vitro. While transplanted hBMSCs/SSCs/HyA-FMBs did not heal articular cartilage defects in immunocompromised rodents, chondrospheroid-derived cells/HyA-FMBs formed non-hypertrophic cartilage that persisted until at least 5 months in vivo.

Cell cycle expression heterogeneity predicts degree of differentiation

Methods that predict fate potential or degree of differentiation from transcriptomic data have identified rare progenitor populations and uncovered developmental regulatory mechanisms. However, some state-of-the-art methods are too computationally burdensome for emerging large-scale data and all methods make inaccurate predictions in certain biological systems. We developed a method in R (stemFinder) that predicts single cell differentiation time based on heterogeneity in cell cycle gene expression. Our method is computationally tractable and is as good as or superior to competitors. As part of our benchmarking, we implemented four different performance metrics to assist potential users in selecting the tool that is most apt for their application. Finally, we explore the relationship between differentiation time and cell fate potential by analyzing a lineage tracing dataset with clonally labelled hematopoietic cells, revealing that metrics of differentiation time are correlated with the number of downstream lineages.

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.

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.

Sox, Fox, and Lmx1b binding sites differentially regulate a Gdf5-Associated regulatory region during elbow development

Introduction: The articulating ends of limb bones have precise morphology and asymmetry that ensures proper joint function. Growth differentiation factor 5 (Gdf5) is a secreted morphogen involved in cartilage and bone development that contributes to the architecture of developing joints. Dysregulation of Gdf5 results in joint dysmorphogenesis often leading to progressive joint degeneration or osteoarthritis (OA). The transcription factors and cis-regulatory modules (CRMs) that regulate Gdf5 expression are not well characterized. We previously identified a Gdf5-associated regulatory region (GARR) that contains predicted binding sites for Lmx1b, Osr2, Fox, and the Sox transcription factors. These transcription factors are recognized factors involved in joint morphogenesis and skeletal development. Methods: We used in situ hybridization to Gdf5, Col2A1, and the transcription factors of interest in developing chicken limbs to determine potential overlap in expression. We further analyzed scRNA-seq data derived from limbs and knees in published mouse and chicken datasets, identifying cells with coexpression of Gdf5 and the transcription factors of interest. We also performed site-directed mutatgenesis of the predicted transcription factor binding sites in a GARR-reporter construct and determined any change in activity using targeted regional electroporation (TREP) in micromass and embryonic chicken wing bioassays. Results: Gdf5 expression overlapped the expression of these transcription factors during joint development both by in situ hybridization (ISH) and scRNA-seq analyses. Within the GARR CRM, mutation of two binding sites common to Fox and Sox transcripstion factors reduced enhancer activity to background levels in micromass cultures and in ovo embryonic chicken wing bioassays, whereas mutation of two Sox-only binding sites caused a significant increase in activity. These results indicate that the Fox/Sox binding sites are required for activity, while the Sox-only sites are involved in repression of activity. Mutation of Lmx1b binding sites in GARR caused an overall reduction in enhancer activity in vitro and a dorsal reduction in ovo. Despite a recognized role for Osr2 in joint development, disruption of the predicted Osr2 site did not alter GARR activity. Conclusion: Taken together, our data indicates that GARR integrates positive, repressive, and asymmetrical inputs to fine-tune the expression of Gdf5 during elbow joint development.

Single-cell RNA sequencing in orthopedic research

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

Application of Single-Cell and Spatial Omics in Musculoskeletal Disorder Research

Musculoskeletal disorders, including fractures, scoliosis, heterotopic ossification, osteoporosis, osteoarthritis, disc degeneration, and muscular injury, etc., can occur at any stage of human life. Understanding the occurrence and development mechanism of musculoskeletal disorders, as well as the changes in tissues and cells during therapy, might help us find targeted treatment methods. Single-cell techniques provide excellent tools for studying alterations at the cellular level of disorders. However, the application of these techniques in research on musculoskeletal disorders is still limited. This review summarizes the current single-cell and spatial omics used in musculoskeletal disorders. Cell isolation, experimental methods, and feasible experimental designs for single-cell studies of musculoskeletal system diseases have been reviewed based on tissue characteristics. Then, the paper summarizes the latest findings of single-cell studies in musculoskeletal disorders from three aspects: bone and ossification, joint, and muscle and tendon disorders. Recent discoveries about the cell populations involved in these diseases are highlighted. Furthermore, the therapeutic responses of musculoskeletal disorders, especially single-cell changes after the treatments of implants, stem cell therapies, and drugs are described. Finally, the application potential and future development directions of single-cell and spatial omics in research on musculoskeletal diseases are discussed.

Deciphering postnatal limb development at single-cell resolution

The early postnatal limb developmental progression bridges embryonic and mature stages and mirrors the pathological remodeling of articular cartilage. However, compared with multitudinous research on embryonic limb development, the early postnatal stage seems relatively unnoticed. Here, a systematic work to portray the postnatal limb developmental landscape was carried out by characterization of 19,952 single cells from murine hindlimbs at 4 postnatal stages using single-cell RNA sequencing technique. By delineation of cell heterogeneity, the candidate progenitor sub-clusters marked by Cd34 and Ly6e were discovered in articular cartilage and enthesis, and three cellular developmental branches marked by Col10a1, Spp1, and Tnni2 were reflected in growth plate. The representative transcriptomes and developmental patterns were intensively explored, and the key regulation mechanisms as well as evolvement in osteoarthritis were discussed. Above all, these results expand horizons of postnatal limb developmental biology and reach the interconnections between limb development, remodeling, and regeneration.

Sirt3 Promotes Chondrogenesis, Chondrocyte Mitochondrial Respiration and the Development of High-Fat Diet-Induced Osteoarthritis in Mice

Understanding how obesity-induced metabolic stress contributes to synovial joint tissue damage is difficult because of the complex role of metabolism in joint development, maintenance, and repair. Chondrocyte mitochondrial dysfunction is implicated in osteoarthritis (OA) pathology, which motivated us to study the mitochondrial deacetylase enzyme sirtuin 3 (Sirt3). We hypothesized that combining high-fat-diet (HFD)-induced obesity and cartilage Sirt3 loss at a young age would impair chondrocyte mitochondrial function, leading to cellular stress and accelerated OA. Instead, we unexpectedly found that depleting cartilage Sirt3 at 5 weeks of age using Sirt3-flox and Acan-Cre^ERT2 mice protected against the development of cartilage degeneration and synovial hyperplasia following 20 weeks of HFD. This protection was associated with increased cartilage glycolysis proteins and reduced mitochondrial fatty acid metabolism proteins. Seahorse-based assays supported a mitochondrial-to-glycolytic shift in chondrocyte metabolism with Sirt3 deletion. Additional studies with primary murine juvenile chondrocytes under hypoxic and inflammatory conditions showed an increased expression of hypoxia-inducible factor (HIF-1) target genes with Sirt3 deletion. However, Sirt3 deletion impaired chondrogenesis using a murine bone marrow stem/stromal cell pellet model, suggesting a context-dependent role of Sirt3 in cartilage homeostasis. Overall, our data indicate that Sirt3 coordinates HFD-induced changes in mature chondrocyte metabolism that promote OA. © 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).

Identification of candidate enhancers controlling the transcriptome during the formation of interphalangeal joints

The formation of the synovial joint begins with the visible emergence of a stripe of densely packed mesenchymal cells located between distal ends of the developing skeletal anlagen called the interzone. Recently the transcriptome of the early synovial joint was reported. Knowledge about enhancers would complement these data and lead to a better understanding of the control of gene transcription at the onset of joint development. Using ChIP-sequencing we have mapped the H3-signatures H3K27ac and H3K4me1 to locate regulatory elements specific for the interzone and adjacent phalange, respectively. This one-stage atlas of candidate enhancers (CEs) was used to map the association between these respective joint tissue specific CEs and biological processes. Subsequently, integrative analysis of transcriptomic data and CEs identified new putative regulatory elements of genes expressed in interzone (e.g., GDF5, BMP2 and DACT2) and phalange (e.g., MATN1, HAPLN1 and SNAI1). We also linked such CEs to genes known as crucial in synovial joint hypermobility and osteoarthritis, as well as phalange malformations. These analyses show that the CE atlas can serve as resource for identifying, and as starting point for experimentally validating, putative disease-causing genomic regulatory regions in patients with synovial joint dysfunctions and/or phalange disorders, and enhancer-controlled synovial joint and phalange formation.

Single-cell multimodal analysis identifies common regulatory programs in synovial fibroblasts of rheumatoid arthritis patients and modeled TNF-driven arthritis

BACKGROUND

Synovial fibroblasts (SFs) are specialized cells of the synovium that provide nutrients and lubricants for the proper function of diarthrodial joints. Recent evidence appreciates the contribution of SF heterogeneity in arthritic pathologies. However, the normal SF profiles and the molecular networks that govern the transition from homeostatic to arthritic SF heterogeneity remain poorly defined.

METHODS

We applied a combined analysis of single-cell (sc) transcriptomes and epigenomes (scRNA-seq and scATAC-seq) to SFs derived from naïve and hTNFtg mice (mice that overexpress human TNF, a murine model for rheumatoid arthritis), by employing the Seurat and ArchR packages. To identify the cellular differentiation lineages, we conducted velocity and trajectory analysis by combining state-of-the-art algorithms including scVelo, Slingshot, and PAGA. We integrated the transcriptomic and epigenomic data to infer gene regulatory networks using ArchR and custom-implemented algorithms. We performed a canonical correlation analysis-based integration of murine data with publicly available datasets from SFs of rheumatoid arthritis patients and sought to identify conserved gene regulatory networks by utilizing the SCENIC algorithm in the human arthritic scRNA-seq atlas.

RESULTS

By comparing SFs from healthy and hTNFtg mice, we revealed seven homeostatic and two disease-specific subsets of SFs. In healthy synovium, SFs function towards chondro- and osteogenesis, tissue repair, and immune surveillance. The development of arthritis leads to shrinkage of homeostatic SFs and favors the emergence of SF profiles marked by Dkk3 and Lrrc15 expression, functioning towards enhanced inflammatory responses and matrix catabolic processes. Lineage inference analysis indicated that specific Thy1+ SFs at the root of trajectories lead to the intermediate Thy1+/Dkk3+/Lrrc15+ SF states and culminate in a destructive and inflammatory Thy1- SF identity. We further uncovered epigenetically primed gene programs driving the expansion of these arthritic SFs, regulated by NFkB and new candidates, such as Runx1. Cross-species analysis of human/mouse arthritic SF data determined conserved regulatory and transcriptional networks.

CONCLUSIONS

We revealed a dynamic SF landscape from health to arthritis providing a functional genomic blueprint to understand the joint pathophysiology and highlight the fibroblast-oriented therapeutic targets for combating chronic inflammatory and destructive arthritic disease.

Regeneration of Jaw Joint Cartilage in Adult Zebrafish

The poor intrinsic repair capacity of mammalian joint cartilage likely contributes to the high incidence of arthritis worldwide. Adult zebrafish can regenerate many structures that show limited or no healing capacity in mammals, including the jawbone. To test whether zebrafish can also regenerate damaged joints, we developed a surgical injury model in which the zebrafish jaw joint is destabilized via transection of the major jaw joint ligament, the interopercular-mandibular (IOM). Unilateral transection of the IOM ligament in 1-year-old fish resulted in an initial reduction of jaw joint cartilage by 14 days, with full regeneration of joint cartilage by 28 days. Joint cartilage regeneration involves the re-entry of articular chondrocytes into the cell cycle and the upregulated expression of sox10, a marker of developing chondrocytes in the embryo that becomes restricted to a subset of joint chondrocytes in adults. Genetic ablation of these sox10-expressing chondrocytes shows that they are essential for joint cartilage regeneration. To uncover the potential source of new chondrocytes during joint regeneration, we performed single-cell RNA sequencing of the uninjured adult jaw joint and identified multiple skeletal, connective tissue, and fibroblast subtypes. In particular, we uncovered a joint-specific periosteal population expressing coch and grem1a, with the jaw joint chondrocytes marked by grem1a expression during regeneration. Our findings demonstrate the capacity of zebrafish to regenerate adult joint cartilage and identify candidate cell types that can be tested for their roles in regenerative response.

Dcx expression defines a subpopulation of Gdf5 + cells with chondrogenic potentials in E12.5 mouse embryonic limbs

Growth differentiation factor 5 (Gdf5) and doublecortin (Dcx) genes are both expressed in joint interzone cells during synovial joint development. In this study, we re-analyzed the single cell RNA-sequencing data (Gene Expression Omnibus GSE151985) generated from Gdf5 ⁺ cells of mouse knee joints at embryonic stages of E12.5, E13.5, E14.5, and E15.5, with a new focus on Dcx. We found that Dcx expression was enriched in clusters of Gdf5 ⁺ cells, with high expression levels of pro-chondrogenic genes including sex determining region Y-box transcription factor 5 (Sox5), Sox6, Sox9, Gdf5, versican, matrilin 4, collagen type II α 1 chain (Col2a1), Col9a1, Col9a2, and Col9a3 at E12.5. Dcx ⁺ and Dcx ^- cells had differential gene expression profiles. The up-regulated genes in Dcx ⁺ vs. Dcx ^- cells at E12.5 and E13.5 were enriched in chondrocyte differentiation and cartilage development, whereas those genes up-regulated at E14.5 and E15.5 were enriched in RNA splicing, protein stability, cell proliferation, and cell growth. Gene expression profiles in Dcx ⁺ cells showed rapid daily changes from E12.5 to E15.5, with limited number of genes shared across the time period. Expression of Gdf5, Sox5, Sox6, melanoma inhibitory activity, noggin, odd-skipped related transcription factor 2, matrilin 4, and versican was positively correlated with Dcx expression. Our results demonstrate that Dcx expression defines a subpopulation of Gdf5 ⁺ cells with chondrogenic potentials in E12.5 mouse embryonic limbs.

Single-Cell RNA-Seq Reveals Transcriptomic Heterogeneity and Post-Traumatic Osteoarthritis-Associated Early Molecular Changes in Mouse Articular Chondrocytes

Articular cartilage is a connective tissue lining the surfaces of synovial joints. When the cartilage severely wears down, it leads to osteoarthritis (OA), a debilitating disease that affects millions of people globally. The articular cartilage is composed of a dense extracellular matrix (ECM) with a sparse distribution of chondrocytes with varying morphology and potentially different functions. Elucidating the molecular and functional profiles of various chondrocyte subtypes and understanding the interplay between these chondrocyte subtypes and other cell types in the joint will greatly expand our understanding of joint biology and OA pathology. Although recent advances in high-throughput OMICS technologies have enabled molecular-level characterization of tissues and organs at an unprecedented resolution, thorough molecular profiling of articular chondrocytes has not yet been undertaken, which may be in part due to the technical difficulties in isolating chondrocytes from dense cartilage ECM. In this study, we profiled articular cartilage from healthy and injured mouse knee joints at a single-cell resolution and identified nine chondrocyte subtypes with distinct molecular profiles and injury-induced early molecular changes in these chondrocytes. We also compared mouse chondrocyte subpopulations to human chondrocytes and evaluated the extent of molecular similarity between mice and humans. This work expands our view of chondrocyte heterogeneity and rapid molecular changes in chondrocyte populations in response to joint trauma and highlights potential mechanisms that trigger cartilage degeneration.

Turning Nature's own processes into design strategies for living bone implant biomanufacturing: a decade of Developmental Engineering

A decade after the term developmental engineering (DE) was coined to indicate the use of developmental processes as blueprints for the design and development of engineered living implants, a myriad of proof-of-concept studies demonstrate the potential of this approach in small animal models. This review provides an overview of DE work, focusing on applications in bone regeneration. Enabling technologies allow to quantify the distance between in vitro processes and their developmental counterpart, as well as to design strategies to reduce that distance. By embedding Nature's robust mechanisms of action in engineered constructs, predictive large animal data and subsequent positive clinical outcomes can be gradually achieved. To this end, the development of next generation biofabrication technologies should provide the necessary scale and precision for robust living bone implant biomanufacturing.