Creb5 establishes the competence for Prg4 expression in articular cartilage

Cell and transcriptomic diversity of infrapatellar fat pad during knee osteoarthritis

OBJECTIVES

In this study, we employ a multiomic approach to identify major cell types and subsets, and their transcriptomic profiles within the infrapatellar fat pad (IFP), and to determine differences in the IFP based on knee osteoarthritis (KOA), sex and obesity status.

METHODS

Single-nucleus RNA sequencing of 82 924 nuclei from 21 IFPs (n=6 healthy control and n=15 KOA donors), spatial transcriptomics and bioinformatic analyses were used to identify contributions of the IFP to KOA. We mapped cell subclusters from other white adipose tissues using publicly available literature. The diversity of fibroblasts within the IFP was investigated by bioinformatic analyses, comparing by KOA, sex and obesity status. Metabolomics was used to further explore differences in fibroblasts by obesity status.

RESULTS

We identified multiple subclusters of fibroblasts, macrophages, adipocytes and endothelial cells with unique transcriptomic profiles. Using spatial transcriptomics, we resolved distributions of cell types and their transcriptomic profiles and computationally identified putative cell-cell communication networks. Furthermore, we identified transcriptomic differences in fibroblasts from KOA versus healthy control donor IFPs, female versus male KOA-IFPs and obese versus normal body mass index (BMI) KOA-IFPs. Finally, using metabolomics, we defined differences in metabolite levels in supernatants of naïve, profibrotic stimuli-treated and proinflammatory stimuli-treated fibroblasts from obese compared to normal BMI KOA-IFPs.

CONCLUSIONS

Overall, by employing a multiomic approach, this study provides the first comprehensive map of the cellular and transcriptomic diversity of human IFP and identifies IFP fibroblasts as key cells contributing to transcriptomic and metabolic differences related to KOA disease, sex or obesity.

Functional genomics of human skeletal development and the patterning of height heritability

Underlying variation in height are regulatory changes to chondrocytes, cartilage cells comprising long-bone growth plates. Currently, we lack knowledge on epigenetic regulation and gene expression of chondrocytes sampled across the human skeleton, and therefore we cannot understand basic regulatory mechanisms controlling height biology. We first rectify this issue by generating extensive epigenetic and transcriptomic maps from chondrocytes sampled from different growth plates across developing human skeletons, discovering novel regulatory networks shaping human bone/joint development. Next, using these maps in tandem with height genome-wide association study (GWAS) signals, we disentangle the regulatory impacts that skeletal element-specific versus global-acting variants have on skeletal growth, revealing the prime importance of regulatory pleiotropy in controlling height variation. Finally, as height is highly heritable, and thus often the test case for complex-trait genetics, we leverage these datasets within a testable omnigenic model framework to discover novel chondrocyte developmental modules and peripheral-acting factors shaping height biology and skeletal growth.

Serinc5 Regulates Sequential Chondrocyte Differentiation by Inhibiting Sox9 Function in Pre-Hypertrophic Chondrocytes

The growth plate is the primary site of longitudinal bone growth with chondrocytes playing a pivotal role in endochondral bone development. Chondrocytes undergo a series of differentiation steps, resulting in the formation of a unique hierarchical columnar structure comprising round, proliferating, pre-hypertrophic, and hypertrophic chondrocytes. Pre-hypertrophic chondrocytes, which exist in the transitional stage between proliferating and hypertrophic stages, are a critical cell population in the growth plate. However, the molecular basis of pre-hypertrophic chondrocytes remains largely undefined. Here, we employed scRNA-seq analysis on fluorescently labeled growth plate chondrocytes for their molecular characterization. Serine incorporator 5 (Serinc5) was identified as a marker gene for pre-hypertrophic chondrocytes. Histological analysis revealed that Serinc5 is specifically expressed in pre-hypertrophic chondrocytes, overlapping with Indian hedgehog (Ihh). Serinc5 represses cell proliferation and Col2a1 and Acan expression by inhibiting the transcriptional activity of Sox9 in primary chondrocytes. Chromatin profiling using ChIP-seq and ATAC-seq revealed an active enhancer of Serinc5 located in intron 1, with its chromatin status progressively activated during chondrocyte differentiation. Collectively, our findings suggest that Serinc5 regulates sequential chondrocyte differentiation from proliferation to hypertrophy by inhibiting Sox9 function in pre-hypertrophic chondrocytes, providing novel insights into the mechanisms underlying chondrocyte differentiation in growth plates.

A conserved transcription factor regulatory program promotes tendon fate

Tendons, which transmit force from muscles to bones, are highly prone to injury. Understanding the mechanisms driving tendon fate would impact efforts to improve tendon healing, yet this knowledge is limited. To find direct regulators of tendon progenitor emergence, we performed a zebrafish high-throughput chemical screen. We established forskolin as a tenogenic inducer across vertebrates, functioning through Creb1a, which is required and sufficient for tendon fate. Putative enhancers containing cyclic AMP (cAMP) response elements (CREs) in humans, mice, and fish drove specific expression in zebrafish cranial and fin tendons. Analysis of these genomic regions identified motifs for early B cell factor (Ebf/EBF) transcription factors. Mutation of CRE or Ebf/EBF motifs significantly disrupted enhancer activity and specificity in tendons. Zebrafish ebf1a/ebf3a mutants displayed defects in tendon formation. Notably, Creb1a/CREB1 and Ebf1a/Ebf3a/EBF1 overexpression facilitated tenogenic induction in zebrafish and human pluripotent stem cells. Together, our work identifies the functional conservation of two transcription factors in promoting tendon fate.

Single cell analysis of Idh mutant growth plates identifies cell populations responsible for longitudinal bone growth and enchondroma formation

Enchondromas are a common tumor in bone that can occur as multiple lesions in enchondromatosis, which is associated with deformity of the affected bone. These lesions harbor somatic mutations in IDH and driving expression of a mutant Idh1 in Col2 expressing cells in mice causes an enchondromatosis phenotype. Here we compared growth plates from E18.5 mice expressing a mutant Idh1 with control littermates using single cell RNA sequencing. Data from Col2 expressing cells were analysed using UMAP and RNA pseudo-time analyses. A unique cluster of cells was identified in the mutant growth plates that expressed genes known to be upregulated in enchondromas. There was also a cluster of cells that was underrepresented in the mutant growth plates that expressed genes known to be important in longitudinal bone growth. Immunofluorescence showed that the genes from the unique cluster identified in the mutant growth plates were expressed in multiple growth plate anatomic zones, and pseudo-time analysis also suggested these cells could arise from multiple growth plate chondrocyte subpopulations. This data supports the notion that a subpopulation of chondrocytes become enchondromas at the expense of contributing to longitudinal growth.

Decoding cellular plasticity and niche regulation of limbal stem cells during corneal wound healing

BACKGROUND

Dysfunction or deficiency of corneal epithelium results in vision impairment or blindness in severe cases. The rapid and effective regeneration of corneal epithelial cells relies on the limbal stem cells (LSCs). However, the molecular and functional responses of LSCs and their niche cells to injury remain elusive.

METHODS

Single-cell RNA sequencing was performed on corneal tissues from normal mice and corneal epithelium defect models. Bioinformatics analysis was performed to confirm the distinct characteristics and cell fates of LSCs. Knockdown of Creb5 and OSM treatment experiment were performed to determine their roles of in corneal epithelial wound healing.

RESULTS

Our data defined the molecular signatures of LSCs and reconstructed the pseudotime trajectory of corneal epithelial cells. Gene network analyses characterized transcriptional landmarks that potentially regulate LSC dynamics, and identified a transcription factor Creb5, that was expressed in LSCs and significantly upregulated after injury. Loss-of-function experiments revealed that silencing Creb5 delayed the corneal epithelial healing and LSC mobilization. Through cell-cell communication analysis, we identified 609 candidate regeneration-associated ligand-receptor interaction pairs between LSCs and distinct niche cells, and discovered a unique subset of Arg1+ macrophages infiltrated after injury, which were present as the source of Oncostatin M (OSM), an IL-6 family cytokine, that were demonstrated to effectively accelerate the corneal epithelial wound healing.

CONCLUSIONS

This research provides a valuable single-cell resource and reference for the discovery of mechanisms and potential clinical interventions aimed at ocular surface reconstruction.

Leveraging single cell multiomic analyses to identify factors that drive human chondrocyte cell fate

Cartilage plays a crucial role in skeletal development and function, and abnormal development contributes to genetic and age-related skeletal disease. To better understand how human cartilage develops in vivo , we jointly profiled the transcriptome and open chromatin regions in individual nuclei recovered from distal femurs at 2 fetal timepoints. We used these multiomic data to identify transcription factors expressed in distinct chondrocyte subtypes, link accessible regulatory elements with gene expression, and predict transcription factor-based regulatory networks that are important for growth plate or epiphyseal chondrocyte differentiation. We developed a human pluripotent stem cell platform for interrogating the function of predicted transcription factors during chondrocyte differentiation and used it to test NFATC2 . We expect new regulatory networks we uncovered using multiomic data to be important for promoting cartilage health and treating disease, and our platform to be a useful tool for studying cartilage development in vitro .

Statement of Significance

The identity and integrity of the articular cartilage lining our joints are crucial to pain-free activities of daily living. Here we identified a gene regulatory landscape of human chondrogenesis at single cell resolution, which is expected to open new avenues of research aimed at mitigating cartilage diseases that affect hundreds of millions of individuals world-wide.

Single-cell transcriptomic analyses of mouse idh1 mutant growth plate chondrocytes reveal distinct cell populations responsible for longitudinal growth and enchondroma formation

Enchondromas are a common tumor in bone that can occur as multiple lesions in enchondromatosis, which is associated with deformity of the effected bone. These lesions harbor mutations in IDH and driving expression of a mutant Idh1 in Col2 expressing cells in mice causes an enchondromatosis phenotype. In this study we compared growth plates from E18.5 mice expressing a mutant Idh1 with control littermates using single cell RNA sequencing. Data from Col2 expressing cells were analyzed using UMAP and RNA pseudo-time analyses. A unique cluster of cells was identified in the mutant growth plates that expressed genes known to be upregulated in enchondromas. There was also a cluster of cells that was underrepresented in the mutant growth plates that expressed genes known to be important in longitudinal bone growth. Immunofluorescence showed that the genes from the unique cluster identified in the mutant growth plates were expressed in multiple growth plate anatomic zones, and pseudo-time analysis also suggested these cells could arise from multiple growth plate chondrocyte subpopulations. This data identifies subpopulations of cells in control and mutant growth plates, and supports the notion that a mutant Idh1 alters the subpopulations of growth plate chondrocytes, resulting a subpopulation of cells that become enchondromas at the expense of other populations that contribute to longitudinal growth.

Prg4-Expressing Chondroprogenitor Cells in the Superficial Zone of Articular Cartilage

Joint-resident chondrogenic precursor cells have become a significant therapeutic option due to the lack of regenerative capacity in articular cartilage. Progenitor cells are located in the superficial zone of the articular cartilage, producing lubricin/Prg4 to decrease friction of cartilage surfaces during joint movement. Prg4-positive progenitors are crucial in maintaining the joint's structure and functionality. The disappearance of progenitor cells leads to changes in articular hyaline cartilage over time, subchondral bone abnormalities, and the formation of ectopic ossification. Genetic labeling cell technology has been the main tool used to characterize Prg4-expressing progenitor cells of articular cartilage in vivo through drug injection at different time points. This technology allows for the determination of the origin of progenitor cells and the tracking of their progeny during joint development and cartilage damage. We endeavored to highlight the currently known information about the Prg4-producing cell population in the joint to underline the significance of the role of these cells in the development of articular cartilage and its homeostasis. This review focuses on superficial progenitors in the joint, how they contribute to postnatal articular cartilage formation, their capacity for regeneration, and the consequences of Prg4 deficiency in these cells. We have accumulated information about the Prg4+ cell population of articular cartilage obtained through various elegantly designed experiments using transgenic technologies to identify potential opportunities for further research.

Widespread variation in molecular interactions and regulatory properties among transcription factor isoforms

Most human Transcription factors (TFs) genes encode multiple protein isoforms differing in DNA binding domains, effector domains, or other protein regions. The global extent to which this results in functional differences between isoforms remains unknown. Here, we systematically compared 693 isoforms of 246 TF genes, assessing DNA binding, protein binding, transcriptional activation, subcellular localization, and condensate formation. Relative to reference isoforms, two-thirds of alternative TF isoforms exhibit differences in one or more molecular activities, which often could not be predicted from sequence. We observed two primary categories of alternative TF isoforms: "rewirers" and "negative regulators", both of which were associated with differentiation and cancer. Our results support a model wherein the relative expression levels of, and interactions involving, TF isoforms add an understudied layer of complexity to gene regulatory networks, demonstrating the importance of isoform-aware characterization of TF functions and providing a rich resource for further studies.

Synovium and infrapatellar fat pad share common mesenchymal progenitors and undergo coordinated changes in osteoarthritis

Osteoarthritis (OA) affects multiple tissues in the knee joint, including the synovium and intra-articular adipose tissue (IAAT) that are attached to each other. However, whether these two tissues share the same progenitor cells and hence function as a single unit in joint homeostasis and diseases is largely unknown. Single-cell transcriptomic profiling of synovium and infrapatellar fat pad (IFP), the largest IAAT, from control and OA mice revealed five mesenchymal clusters and predicted mesenchymal progenitor cells (MPCs) as the common progenitors for other cells: synovial lining fibroblasts (SLFs), myofibroblasts (MFs), and preadipocytes 1 and 2. Histologic examination of joints in reporter mice having Dpp4-CreER and Prg4-CreER that label MPCs and SLFs, respectively, demonstrated that Dpp4+ MPCs reside in the synovial sublining layer and give rise to Prg4+ SLFs and Perilipin+ adipocytes during growth and OA progression. After OA injury, both MPCs and SLFs gave rise to MFs, which remained in the thickened synovium at later stages of OA. In culture, Dpp4+ MPCs possessed mesenchymal progenitor properties, such as proliferation and multilineage differentiation. In contrast, Prg4+ SLFs did not contribute to adipocytes in IFP and Prg4+ cells barely grew in vitro. Taken together, we demonstrate that the synovium and joint fat pad are one integrated functional tissue sharing common mesenchymal progenitors and undergoing coordinated changes during OA progression.

The Actin Cytoskeleton as a Regulator of Proteoglycan 4

OBJECTIVE

The superficial zone (SZ) of articular cartilage is responsible for distributing shear forces for optimal cartilage loading and contributes to joint lubrication through the production of proteoglycan 4 (PRG4). PRG4 plays a critical role in joint homeostasis and is chondroprotective. Normal PRG4 production is affected by inflammation and irregular mechanical loading in post-traumatic osteoarthritis (PTOA). THe SZ chondrocyte (SZC) phenotype, including PRG4 expression, is regulated by the actin cytoskeleton in vitro. There remains a limited understanding of the regulation of PRG4 by the actin cytoskeleton in native articular chondrocytes. The filamentous (F)-actin cytoskeleton is a potential node in crosstalk between mechanical stimulation and cytokine activation and the regulation of PRG4 in SZCs, therefore developing insights in the regulation of PRG4 by actin may identify molecular targets for novel PTOA therapies.

MATERIALS AND METHODS

A comprehensive literature search on PRG4 and the regulation of the SZC phenotype by actin organization was performed.

RESULTS

PRG4 is strongly regulated by the actin cytoskeleton in isolated SZCs in vitro. Biochemical and mechanical stimuli have been characterized to regulate PRG4 and may converge upon actin cytoskeleton signaling.

CONCLUSION

Actin-based regulation of PRG4 in native SZCs is not fully understood and requires further elucidation. Understanding the regulation of PRG4 by actin in SZCs requires an in vivo context to further potential of leveraging actin arrangement to arthritic therapeutics.

Altered Sox9 and FGF signaling gene expression in Aga2 OI mice negatively affects linear growth

Osteogenesis imperfecta (OI), or brittle bone disease, is a disorder characterized by bone fragility and increased fracture incidence. All forms of OI also feature short stature, implying an effect on endochondral ossification. Using the Aga2+/- mouse, which has a mutation in type I collagen, we show an affected growth plate primarily due to a shortened proliferative zone. We used single-cell RNA-Seq analysis of tibial and femoral growth plate tissues to understand transcriptional consequences on growth plate cell types. We show that perichondrial cells, which express abundant type I procollagen, and growth plate chondrocytes, which were found to express low amounts of type I procollagen, had ER stress and dysregulation of the same unfolded protein response pathway as previously demonstrated in osteoblasts. Aga2+/- proliferating chondrocytes showed increased FGF and MAPK signaling, findings consistent with accelerated differentiation. There was also increased Sox9 expression throughout the growth plate, which is expected to accelerate early chondrocyte differentiation but reduce late hypertrophic differentiation. These data reveal that mutant type I collagen expression in OI has an impact on the cartilage growth plate. These effects on endochondral ossification indicate that OI is a biologically complex phenotype going beyond its known impacts on bone to negatively affect linear growth.

Single-cell, whole-embryo phenotyping of mammalian developmental disorders

Mouse models are a critical tool for studying human diseases, particularly developmental disorders1. However, conventional approaches for phenotyping may fail to detect subtle defects throughout the developing mouse2. Here we set out to establish single-cell RNA sequencing of the whole embryo as a scalable platform for the systematic phenotyping of mouse genetic models. We applied combinatorial indexing-based single-cell RNA sequencing3 to profile 101 embryos of 22 mutant and 4 wild-type genotypes at embryonic day 13.5, altogether profiling more than 1.6 million nuclei. The 22 mutants represent a range of anticipated phenotypic severities, from established multisystem disorders to deletions of individual regulatory regions4,5. We developed and applied several analytical frameworks for detecting differences in composition and/or gene expression across 52 cell types or trajectories. Some mutants exhibit changes in dozens of trajectories whereas others exhibit changes in only a few cell types. We also identify differences between widely used wild-type strains, compare phenotyping of gain- versus loss-of-function mutants and characterize deletions of topological associating domain boundaries. Notably, some changes are shared among mutants, suggesting that developmental pleiotropy might be 'decomposable' through further scaling of this approach. Overall, our findings show how single-cell profiling of whole embryos can enable the systematic molecular and cellular phenotypic characterization of mouse mutants with unprecedented breadth and resolution.

CaSR modulates proliferation of the superficial zone cells in temporomandibular joint cartilage via the PTHrP nuclear localization sequence

The superficial zone cells in mandibular condylar cartilage are proliferative. The present purpose was to delineate the relation of calcium-sensing receptor (CaSR) and parathyroid hormone-related peptide nuclear localization sequence (PTHrP87-139 ), and their role in the proliferation behaviors of the superficial zone cells. A gain- and loss-of-function strategy were used in an in vitro fluid flow shear stress (FFSS) model and an in vivo bilateral elevation bite model which showed mandibular condylar cartilage thickening. CaSR and PTHrP87-139 were modulated through treating the isolated superficial zone cells with activator/SiRNA and via deleting CaSR or parathyroid hormone-related peptide (PTHrP) gene in mice with the promoter gene of proteoglycan 4 (Prg4-CreERT2 ) in the tamoxifen-inducible pattern with or without additional injection of Cinacalcet, the CaSR agonist, or PTHrP87-139 peptide. FFSS stimulated CaSR and PTHrP expression, and accelerated proliferation of the Prg4-expressing superficial zone cells, in which process CaSR acted as an up-streamer of PTHrP. Proteoglycan 4 specific knockout of CaSR or PTHrP reduced the cartilage thickness, suppressed the proliferation and early differentiation of the superficial zone cells, and inhibited cartilage thickening and matrix production promoted by bilateral elevation bite. Injections of CaSR agonist Cinacalcet could not improve the phenotype caused by PTHrP mutation. Injections of PTHrP87-139 peptide rescued the cartilage from knockout of CaSR gene. CaSR modulates proliferation of the superficial zone cells in mandibular condylar cartilage through activation of PTHrP nuclear localization sequence. Our data support the therapeutic target of CaSR in promoting PTHrP production in superficial zone cartilage.

Lineage-specific differences and regulatory networks governing human chondrocyte development

To address large gaps in our understanding of the molecular regulation of articular and growth plate cartilage development in humans, we used our directed differentiation approach to generate these distinct cartilage tissues from human embryonic stem cells. The resulting transcriptomic profiles of hESC-derived articular and growth plate chondrocytes were similar to fetal epiphyseal and growth plate chondrocytes, with respect to genes both known and previously unknown to cartilage biology. With the goal to characterize the regulatory landscapes accompanying these respective transcriptomes, we mapped chromatin accessibility in hESC-derived chondrocyte lineages, and mouse embryonic chondrocytes, using ATAC-sequencing. Integration of the expression dataset with the differentially accessible genomic regions revealed lineage-specific gene regulatory networks. We validated functional interactions of two transcription factors (TFs) (RUNX2 in growth plate chondrocytes and RELA in articular chondrocytes) with their predicted genomic targets. The maps we provide thus represent a framework for probing regulatory interactions governing chondrocyte differentiation. This work constitutes a substantial step towards comprehensive and comparative molecular characterizations of distinct chondrogenic lineages and sheds new light on human cartilage development and biology.

Digits in a dish: An in vitro system to assess the molecular genetics of hand/foot development at single-cell resolution

In vitro models allow for the study of developmental processes outside of the embryo. To gain access to the cells mediating digit and joint development, we identified a unique property of undifferentiated mesenchyme isolated from the distal early autopod to autonomously re-assemble forming multiple autopod structures including: digits, interdigital tissues, joints, muscles and tendons. Single-cell transcriptomic analysis of these developing structures revealed distinct cell clusters that express canonical markers of distal limb development including: Col2a1, Col10a1, and Sp7 (phalanx formation), Thbs2 and Col1a1 (perichondrium), Gdf5, Wnt5a, and Jun (joint interzone), Aldh1a2 and Msx1 (interdigital tissues), Myod1 (muscle progenitors), Prg4 (articular perichondrium/articular cartilage), and Scx and Tnmd (tenocytes/tendons). Analysis of the gene expression patterns for these signature genes indicates that developmental timing and tissue-specific localization were also recapitulated in a manner similar to the initiation and maturation of the developing murine autopod. Finally, the in vitro digit system also recapitulates congenital malformations associated with genetic mutations as in vitro cultures of Hoxa13 mutant mesenchyme produced defects present in Hoxa13 mutant autopods including digit fusions, reduced phalangeal segment numbers, and poor mesenchymal condensation. These findings demonstrate the robustness of the in vitro digit system to recapitulate digit and joint development. As an in vitro model of murine digit and joint development, this innovative system will provide access to the developing limb tissues facilitating studies to discern how digit and articular joint formation is initiated and how undifferentiated mesenchyme is patterned to establish individual digit morphologies. The in vitro digit system also provides a platform to rapidly evaluate treatments aimed at stimulating the repair or regeneration of mammalian digits impacted by congenital malformation, injury, or disease.

Taxonomy of fibroblasts and progenitors in the synovial joint at single-cell resolution

OBJECTIVES

Fibroblasts in synovium include fibroblast-like synoviocytes (FLS) in the lining and Thy1+ connective-tissue fibroblasts in the sublining. We aimed to investigate their developmental origin and relationship with adult progenitors.

METHODS

To discriminate between Gdf5-lineage cells deriving from the embryonic joint interzone and other Pdgfrα-expressing fibroblasts and progenitors, adult Gdf5-Cre;Tom;Pdgfrα-H2BGFP mice were used and cartilage injury was induced to activate progenitors. Cells were isolated from knees, fibroblasts and progenitors were sorted by fluorescence-activated cell-sorting based on developmental origin, and analysed by single-cell RNA-sequencing. Flow cytometry and immunohistochemistry were used for validation. Clonal-lineage mapping was performed using Gdf5-Cre;Confetti mice.

RESULTS

In steady state, Thy1+ sublining fibroblasts were of mixed ontogeny. In contrast, Thy1-Prg4+ lining fibroblasts predominantly derived from the embryonic joint interzone and included Prg4-expressing progenitors distinct from molecularly defined FLS. Clonal-lineage tracing revealed compartmentalisation of Gdf5-lineage fibroblasts between lining and sublining. Following injury, lining hyperplasia resulted from proliferation and differentiation of Prg4-expressing progenitors, with additional recruitment of non-Gdf5-lineage cells, into FLS. Consistent with this, a second population of proliferating cells, enriched near blood vessels in the sublining, supplied activated multipotent cells predicted to give rise to Thy1+ fibroblasts, and to feed into the FLS differentiation trajectory. Transcriptional programmes regulating fibroblast differentiation trajectories were uncovered, identifying Sox5 and Foxo1 as key FLS transcription factors in mice and humans.

CONCLUSIONS

Our findings blueprint a cell atlas of mouse synovial fibroblasts and progenitors in healthy and injured knees, and provide novel insights into the cellular and molecular principles governing the organisation and maintenance of adult synovial joints.

Engraftment of allogeneic iPS cell-derived cartilage organoid in a primate model of articular cartilage defect

Induced pluripotent stem cells (iPSCs) are a promising resource for allogeneic cartilage transplantation to treat articular cartilage defects that do not heal spontaneously and often progress to debilitating conditions, such as osteoarthritis. However, to the best of our knowledge, allogeneic cartilage transplantation into primate models has never been assessed. Here, we show that allogeneic iPSC-derived cartilage organoids survive and integrate as well as are remodeled as articular cartilage in a primate model of chondral defects in the knee joints. Histological analysis revealed that allogeneic iPSC-derived cartilage organoids in chondral defects elicited no immune reaction and directly contributed to tissue repair for at least four months. iPSC-derived cartilage organoids integrated with the host native articular cartilage and prevented degeneration of the surrounding cartilage. Single-cell RNA-sequence analysis indicated that iPSC-derived cartilage organoids differentiated after transplantation, acquiring expression of PRG4 crucial for joint lubrication. Pathway analysis suggested the involvement of SIK3 inactivation. Our study outcomes suggest that allogeneic transplantation of iPSC-derived cartilage organoids may be clinically applicable for the treatment of patients with chondral defects of the articular cartilage; however further assessment of functional recovery long term after load bearing injuries is required.

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.

TGF-β signaling and Creb5 cooperatively regulate Fgf18 to control pharyngeal muscle development

The communication between myogenic cells and their surrounding connective tissues is indispensable for muscle morphogenesis. During late embryonic development in mice, myogenic progenitors migrate to discrete sites to form individual muscles. The detailed mechanism of this process remains unclear. Using mouse levator veli palatini (LVP) development as a model, we systematically investigated how a distinct connective tissue subpopulation, perimysial fibroblasts, communicates with myogenic cells to regulate mouse pharyngeal myogenesis. Using single-cell RNAseq data analysis, we identified that TGF-β signaling is a key regulator for the perimysial fibroblasts. Loss of TGF-β signaling in the neural crest-derived palatal mesenchyme leads to defects in perimysial fibroblasts and muscle malformation in the soft palate in Osr2Cre;Tgfbr1fl/fl mice. In particular, Creb5, a transcription factor expressed in the perimysial fibroblasts, cooperates with TGF-β signaling to activate expression of Fgf18. Moreover, Fgf18 supports pharyngeal muscle development in vivo and exogenous Fgf18 can partially rescue myogenic cell numbers in Osr2Cre;Tgfbr1fl/fl samples, illustrating that TGF-β-regulated Fgf18 signaling is required for LVP development. Collectively, our findings reveal the mechanism by which TGF-β signaling achieves its functional specificity in defining the perimysial-to-myogenic signals for pharyngeal myogenesis.

Creb5 coordinates synovial joint formation with the genesis of articular cartilage

While prior work has established that articular cartilage arises from Prg4-expressing perichondrial cells, it is not clear how this process is specifically restricted to the perichondrium of synovial joints. We document that the transcription factor Creb5 is necessary to initiate the expression of signaling molecules that both direct the formation of synovial joints and guide perichondrial tissue to form articular cartilage instead of bone. Creb5 promotes the generation of articular chondrocytes from perichondrial precursors in part by inducing expression of signaling molecules that block a Wnt5a autoregulatory loop in the perichondrium. Postnatal deletion of Creb5 in the articular cartilage leads to loss of both flat superficial zone articular chondrocytes coupled with a loss of both Prg4 and Wif1 expression in the articular cartilage; and a non-cell autonomous up-regulation of Ctgf. Our findings indicate that Creb5 promotes joint formation and the subsequent development of articular chondrocytes by driving the expression of signaling molecules that both specify the joint interzone and simultaneously inhibit a Wnt5a positive-feedback loop in the perichondrium.

RSPO2 defines a distinct undifferentiated progenitor in the tendon/ligament and suppresses ectopic ossification

Ectopic endochondral ossification in the tendon/ligament is caused by repetitive mechanical overload or inflammation. Tendon stem/progenitor cells (TSPCs) contribute to tissue repair, and some express lubricin [proteoglycan 4 (PRG4)]. However, the mechanisms of ectopic ossification and association of TSPCs are not yet known. Here, we investigated the characteristics of Prg4-positive (+) cells and identified that R-spondin 2 (RSPO2), a WNT activator, is specifically expressed in a distinct Prg4+ TSPC cluster. The Rspo2+ cluster was characterized as mostly undifferentiated, and RSPO2 overexpression suppressed ectopic ossification in a mouse Achilles tendon puncture model via chondrogenic differentiation suppression. RSPO2 expression levels in patients with ossification of the posterior longitudinal ligament were lower than those in spondylosis patients, and RSPO2 protein suppressed chondrogenic differentiation of human ligament cells. RSPO2 was induced by inflammatory stimulation and mechanical loading via nuclear factor κB. Rspo2+ cells may contribute to tendon/ligament homeostasis under pathogenic conditions.

Roles of Cartilage-Resident Stem/Progenitor Cells in Cartilage Physiology, Development, Repair and Osteoarthritis

Osteoarthritis (OA) is a degenerative disease that causes irreversible destruction of articular cartilage for which there is no effective treatment at present. Although articular cartilage lacks intrinsic reparative capacity, numerous studies have confirmed the existence of cartilage-resident stem/progenitor cells (CSPCs) in the superficial zone (SFZ) of articular cartilage. CSPCs are characterized by the expression of mesenchymal stromal cell (MSC)-related surface markers, multilineage differentiation ability, colony formation ability, and migration ability in response to injury. In contrast to MSCs and chondrocytes, CSPCs exhibit extensive proliferative and chondrogenic potential with no signs of hypertrophic differentiation, highlighting them as suitable cell sources for cartilage repair. In this review, we focus on the organizational distribution, markers, cytological features and roles of CSPCs in cartilage development, homeostasis and repair, and the application potential of CSPCs in cartilage repair and OA therapies.

Krüppel-like factor-4 and Krüppel-like factor-2 are important regulators of joint tissue cells and protect against tissue destruction and inflammation in osteoarthritis

OBJECTIVES

Analysing expression patterns of Krüppel-like factor (KLF) transcription factors in normal and osteoarthritis (OA) human cartilage, and determining functions and mechanisms of KLF4 and KLF2 in joint homoeostasis and OA pathogenesis.

METHODS

Experimental approaches included human joint tissues cells, transgenic mice and mouse OA model with viral KLF4 gene delivery to demonstrate therapeutic benefit in structure and pain improvement. Mechanistic studies applied global gene expression analysis and chromatin immunoprecipitation sequencing (ChIP-seq).

RESULTS

Several KLF genes were significantly decreased in OA cartilage. Among them, KLF4 and KLF2 were strong inducers of cartilage collagen genes and Proteoglycan-4. Cartilage-specific deletion of Klf2 in mature mice aggravated severity of experimental OA. Transduction of human chondrocytes with Adenovirus (Ad) expressing KLF4 or KLF2 enhanced expression of major cartilage extracellular matrix (ECM) genes and SRY-box transcription factor-9, and suppressed mediators of inflammation and ECM-degrading enzymes. Ad-KLF4 and Ad-KLF2 enhanced similar protective functions in meniscus cells and synoviocytes, and promoted chondrocytic differentiation of human mesenchymal stem cells. Viral KLF4 delivery into mouse knees reduced severity of OA-associated changes in cartilage, meniscus and synovium, and improved pain behaviours. ChIP-seq analysis suggested that KLF4 directly bound cartilage signature genes. Ras-related protein-1 signalling was the most enriched pathway in KLF4-transduced cells, and its signalling axis was involved in upregulating cartilage ECM genes by KLF4 and KLF2.

CONCLUSIONS

KLF4 and KLF2 may be central transcription factors that increase protective and regenerative functions in joint tissue cells, suggesting that KLF gene transfer or molecules upregulating KLFs are therapeutic candidates for OA.

A novel CREB5/TOP1MT axis confers cisplatin resistance through inhibiting mitochondrial apoptosis in head and neck squamous cell carcinoma

BACKGROUND

Cisplatin resistance is one of the main causes of treatment failure and death in head and neck squamous cell carcinoma (HNSCC). A more comprehensive understanding of the cisplatin resistance mechanism and the development of effective treatment strategies are urgent.

METHODS

RNA sequencing, RT-PCR, and immunoblotting were used to identify differentially expressed genes associated with cisplatin resistance. Gain- and loss-of-function experiments were performed to detect the effect of CREB5 on cisplatin resistance and mitochondrial apoptosis in HNSCC. Chromatin immunoprecipitation (ChIP) assay, dual-luciferase reporter assay, and immunoblotting experiments were performed to explore the underlying mechanisms of CREB5.

RESULTS

CREB5 was significantly upregulated in cisplatin-resistant HNSCC (CR-HNSCC) patients, which was correlated with poor prognosis. CREB5 overexpression strikingly facilitated the cisplatin resistance of HNSCC cells in vitro and in vivo, while CREB5 knockdown enhanced cisplatin sensitivity in CR-HNSCC cells. Interestingly, the activation of AKT signaling induced by cisplatin promoted nucleus translocation of CREB5 in CR-HNSCC cells. Furthermore, CREB5 transcriptionally activated TOP1MT expression depending on the canonical motif. Moreover, CREB5 silencing could trigger mitochondrial apoptosis and overcome cisplatin resistance in CR-HNSCC cells, which could be reversed by TOP1MT overexpression. Additionally, double-targeting of CREB5 and TOP1MT could combat cisplatin resistance of HNSCC in vivo.

CONCLUSIONS

Our findings reveal a novel CREB5/TOP1MT axis conferring cisplatin resistance in HNSCC, which provides a new basis to develop effective strategies for overcoming cisplatin resistance.

Proteoglycan 4 (PRG4) treatment enhances wound closure and tissue regeneration

The wound healing response is one of most primitive and conserved physiological responses in the animal kingdom, as restoring tissue integrity/homeostasis can be the difference between life and death. Wound healing in mammals is mediated by immune cells and inflammatory signaling molecules that regulate tissue resident cells, including local progenitor cells, to mediate closure of the wound through formation of a scar. Proteoglycan 4 (PRG4), a protein found throughout the animal kingdom from fish to elephants, is best known as a glycoprotein that reduces friction between articulating surfaces (e.g. cartilage). Previously, PRG4 was also shown to regulate the inflammatory and fibrotic response. Based on this, we asked whether PRG4 plays a role in the wound healing response. Using an ear wound model, topical application of exogenous recombinant human (rh)PRG4 hastened wound closure and enhanced tissue regeneration. Our results also suggest that rhPRG4 may impact the fibrotic response, angiogenesis/blood flow to the injury site, macrophage inflammatory dynamics, recruitment of immune and increased proliferation of adult mesenchymal progenitor cells (MPCs) and promoting chondrogenic differentiation of MPCs to form the auricular cartilage scaffold of the injured ear. These results suggest that PRG4 has the potential to suppress scar formation while enhancing connective tissue regeneration post-injury by modulating aspects of each wound healing stage (blood clotting, inflammation, tissue generation and tissue remodeling). Therefore, we propose that rhPRG4 may represent a potential therapy to mitigate scar and improve wound healing.

Primary Cilia Direct Murine Articular Cartilage Tidemark Patterning Through Hedgehog Signaling and Ambulatory Load

Articular cartilage (AC) is essential for body movement but is highly susceptible to degenerative diseases and has poor self-repair capacity. To improve current subpar regenerative treatments, developmental mechanisms of AC should be clarified and, specifically, how its postnatal multizone organization is acquired. Primary cilia are cell surface organelles crucial for mammalian tissue morphogenesis. Although their importance for chondrocyte function is appreciated, their specific roles in postnatal AC morphogenesis remain unclear. To explore these mechanisms, we used a murine conditional loss-of-function approach (Ift88-flox) targeting joint-lineage progenitors (Gdf5Cre) and monitored postnatal knee AC development. Joint formation and growth up to juvenile stages were largely unaffected. However, mature AC (aged 2 months) exhibited disorganized extracellular matrix, decreased aggrecan and collagen II due to reduced gene expression (not increased catabolism), and marked reduction of AC modulus by 30%-50%. In addition, and unexpectedly, we discovered that tidemark patterning was severely disrupted, as was hedgehog signaling, and exhibited specificity based on regional load-bearing functions of AC. Interestingly, Prg4 expression was markedly increased in highly loaded sites in mutants. Together, our data provide evidence that primary cilia orchestrate postnatal AC morphogenesis including tidemark topography, zonal matrix composition, and ambulation load responses. © 2022 American Society for Bone and Mineral Research (ASBMR).

The superficial zone of articular cartilage

The superficial zone of articular cartilage contributes to smooth joint motion through the production of proteoglycan 4 (PRG4), also known as lubricin. Recent studies indicate novel effects of PRG4 as a signaling molecule, other than a simple extracellular matrix protein. Additionally, the accumulating evidence displays that various molecules and signaling pathways are involved in regulating the superficial zone and PRG4 expression. In addition, Prg4-expressing cells include a progenitor population of articular chondrocytes. Several non-clinical and clinical studies have shown that PRG4 and related molecules are promising candidates for disease-modifying drugs for treating osteoarthritis. Since PRG4 is also expressed in the synovium, tendons, and ligaments, further studies of PRG4-related pathways and PRG4-positive cells may elucidate the mechanisms underlying joint homeostasis.

Regulatory Mechanisms of Prg4 and Gdf5 Expression in Articular Cartilage and Functions in Osteoarthritis

Owing to the rapid aging of society, the numbers of patients with joint disease continue to increase. Accordingly, a large number of patients require appropriate treatment for osteoarthritis (OA), the most frequent bone and joint disease. Thought to be caused by the degeneration and destruction of articular cartilage following persistent and excessive mechanical stimulation of the joints, OA can significantly impair patient quality of life with symptoms such as knee pain, lower limb muscle weakness, or difficulty walking. Because articular cartilage has a low self-repair ability and an extremely low proliferative capacity, healing of damaged articular cartilage has not been achieved to date. The current pharmaceutical treatment of OA is limited to the slight alleviation of symptoms (e.g., local injection of hyaluronic acid or non-steroidal anti-inflammatory drugs); hence, the development of effective drugs and regenerative therapies for OA is highly desirable. This review article summarizes findings indicating that proteoglycan 4 (Prg4)/lubricin, which is specifically expressed in the superficial zone of articular cartilage and synovium, functions in a protective manner against OA, and covers the transcriptional regulation of Prg4 in articular chondrocytes. We also focused on growth differentiation factor 5 (Gdf5), which is specifically expressed on the surface layer of articular cartilage, particularly in the developmental stage, describing its regulatory mechanisms and functions in joint formation and OA pathogenesis. Because several genetic studies in humans and mice indicate the involvement of these genes in the maintenance of articular cartilage homeostasis and the presentation of OA, molecular targeting of Prg4 and Gdf5 is expected to provide new insights into the aetiology, pathogenesis, and potential treatment of OA.

RGD v2.0: a major update of the ruminant functional and evolutionary genomics database

Ruminant Genome Database (RGD; http://animal.nwsuaf.edu.cn/RGD) provides visualization and analysis tools for ruminant comparative genomics and functional annotations. As more high-quality ruminant genome assemblies have become available, we have redesigned the user interface, integrated and expanded multi-omics data, and developed novel features to improve the database. The new version, RGD v2.0, houses 78 ruminant genomes; 110-species synteny alignments for major livestock (including cattle, sheep, goat) and wild ungulates; 21 012 orthologous gene clusters with Gene Ontology and pathway annotation; ∼8 600 000 conserved elements; and ∼1 000 000 cis-regulatory elements by utilizing 1053 epigenomic data sets. The transcriptome data in RGD v2.0 has nearly doubled, currently with 1936 RNA-seq data sets, and 155 174 phenotypic data sets have been newly added. New and updated features include: (i) The UCSC Genome Browser, BLAT, BLAST and Table Browser tools were updated for six available ruminant livestock species. (ii) The LiftOver tool was newly introduced into our browser to allow coordinate conversion between different ruminant assemblies. And (iii) tissue specificity index, tau, was calculated to facilitate batch screening of specifically expressed genes. The enhanced genome annotations and improved functionality in RGD v2.0 will be useful for study of genome evolution, environmental adaption, livestock breeding and biomedicine.