A single-cell atlas of conventional central chondrosarcoma reveals the role of endoplasmic reticulum stress in malignant transformation.

The transformation of benign lesions to malignant tumours is a crucial aspect of understanding chondrosarcomas, which are malignant cartilage tumours that could develop from benign chondroid lesions. However, the process of malignant transformation for chondroid lesions remains poorly understood, and no reliable markers are available to aid clinical decision-making. To address this issue, we conducted a study analysing 11 primary cartilage tumours and controls using single-cell RNA sequencing. By creating a single-cell atlas, we were able to identify the role of endoplasmic reticulum (ER) stress in the malignant transformation of conventional central chondrosarcomas (CCCS). Our research revealed that lower levels of ER stress promote chondrosarcoma growth in a patient-derived xenograft mouse model, while intensive ER stress reduces primary chondrosarcoma cell viability. Furthermore, we discovered that the NF-κB pathway alleviates ER stress-induced apoptosis during chondrosarcoma progression. Our single-cell signatures and large public data support the use of key ER stress regulators, such as DNA Damage Inducible Transcript 3 (DDIT3; also known as CHOP), as malignant markers for overall patient survival. Ultimately, our study highlights the significant role that ER stress plays in the malignant transformation of cartilaginous tumours and provides a valuable resource for future diagnostic markers and therapeutic strategies.

MMP14 cleaves PTH1R in the chondrocyte-derived osteoblast lineage, curbing signaling intensity for proper bone anabolism.

Bone homeostasis is regulated by hormones such as parathyroid hormone (PTH). While PTH can stimulate osteo-progenitor expansion and bone synthesis, how the PTH-signaling intensity in progenitors is controlled is unclear. Endochondral bone osteoblasts arise from perichondrium-derived osteoprogenitors and hypertrophic chondrocytes (HC). We found, via single-cell transcriptomics, that HC-descendent cells activate membrane-type 1 metalloproteinase 14 (MMP14) and the PTH pathway as they transition to osteoblasts in neonatal and adult mice. Unlike Mmp14 global knockouts, postnatal day 10 (p10) HC lineage-specific Mmp14 null mutants (Mmp14ΔHC) produce more bone. Mechanistically, MMP14 cleaves the extracellular domain of PTH1R, dampening PTH signaling, and consistent with the implied regulatory role, in Mmp14ΔHC mutants, PTH signaling is enhanced. We found that HC-derived osteoblasts contribute ~50% of osteogenesis promoted by treatment with PTH 1-34, and this response was amplified in Mmp14ΔHC. MMP14 control of PTH signaling likely applies also to both HC- and non-HC-derived osteoblasts because their transcriptomes are highly similar. Our study identifies a novel paradigm of MMP14 activity-mediated modulation of PTH signaling in the osteoblast lineage, contributing new insights into bone metabolism with therapeutic significance for bone-wasting diseases.

Hedgehog signaling orchestrates cartilage-to-bone transition independently of Smoothened.

Although recent lineage studies strongly support a chondrocyte-to-osteoblast differentiation continuum, the biological significance and molecular basis remain undetermined. In silico analysis at a single-cell level indicates a transient shutdown of Hedgehog-related transcriptome during simulated cartilage-to-bone transition. Prompted by this, we genetically induce gain- and loss-of function to probe the role of Hedgehog signaling in cartilage-to-bone transition. Ablating Smo in hypertrophic chondrocytes (HCs) does not result in any phenotypic outcome, whereas deleting Ptch1 in HCs leads to disrupted formation of primary spongiosa and actively proliferating HCs-derived osteogenic cells that contribute to bony bulges seen in adult mutant mice. In HCs-derived osteoblasts, constitutive activation of Hedgehog signaling blocks their further differentiation to osteocytes. Moreover, ablation of both Smo and Ptch1 in HCs reverses neither persistent Hedgehog signaling nor bone overgrowths. These results establish a functional contribution of extended chondrocyte lineage to bone homeostasis and diseases, governed by an unanticipated mode of regulation for Hedgehog signaling independently of Smo.

Hypertrophic chondrocytes serve as a reservoir for marrow-associated skeletal stem and progenitor cells, osteoblasts, and adipocytes during skeletal development.

Hypertrophic chondrocytes give rise to osteoblasts during skeletal development; however, the process by which these non-mitotic cells make this transition is not well understood. Prior studies have also suggested that skeletal stem and progenitor cells (SSPCs) localize to the surrounding periosteum and serve as a major source of marrow-associated SSPCs, osteoblasts, osteocytes, and adipocytes during skeletal development. To further understand the cell transition process by which hypertrophic chondrocytes contribute to osteoblasts or other marrow associated cells, we utilized inducible and constitutive hypertrophic chondrocyte lineage tracing and reporter mouse models (Col10a1CreERT2; Rosa26fs-tdTomato and Col10a1Cre; Rosa26fs-tdTomato) in combination with a PDGFRaH2B-GFP transgenic line, single-cell RNA-sequencing, bulk RNA-sequencing, immunofluorescence staining, and cell transplantation assays. Our data demonstrate that hypertrophic chondrocytes undergo a process of dedifferentiation to generate marrow-associated SSPCs that serve as a primary source of osteoblasts during skeletal development. These hypertrophic chondrocyte-derived SSPCs commit to a CXCL12-abundant reticular (CAR) cell phenotype during skeletal development and demonstrate unique abilities to recruit vasculature and promote bone marrow establishment, while also contributing to the adipogenic lineage.

DIPPER, a spatiotemporal proteomics atlas of human intervertebral discs for exploring ageing and degeneration dynamics.

The spatiotemporal proteome of the intervertebral disc (IVD) underpins its integrity and function. We present DIPPER, a deep and comprehensive IVD proteomic resource comprising 94 genome-wide profiles from 17 individuals. To begin with, protein modules defining key directional trends spanning the lateral and anteroposterior axes were derived from high-resolution spatial proteomes of intact young cadaveric lumbar IVDs. They revealed novel region-specific profiles of regulatory activities and displayed potential paths of deconstruction in the level- and location-matched aged cadaveric discs. Machine learning methods predicted a 'hydration matrisome' that connects extracellular matrix with MRI intensity. Importantly, the static proteome used as point-references can be integrated with dynamic proteome (SILAC/degradome) and transcriptome data from multiple clinical samples, enhancing robustness and clinical relevance. The data, findings, and methodology, available on a web interface (http://www.sbms.hku.hk/dclab/DIPPER/), will be valuable references in the field of IVD biology and proteomic analytics.

The backbone of vertebrate animals consists of a series of bones called vertebrae that are joined together by disc-like structures that allow the back to move and distribute forces to protect it during daily activities. It is common for these intervertebral discs to degenerate with age, resulting in back pain and severely reducing quality of life. The mechanical features of intervertebral discs are the result of their proteins. These include extracellular matrix proteins, which form the external scaffolding that binds cells together in a tissue, and signaling proteins, which allow cells to communicate. However, how the levels of different proteins in each region of the disc vary with time has not been fully examined. To establish how protein composition changes with age, Tam, Chen et al. quantified the protein levels and gene activity (which leads to protein production) of intervertebral discs from young and old deceased individuals. They found that the position of different mixtures of proteins in the intervertebral disc changes with age, and that young people have high levels of extracellular matrix proteins and signaling proteins. Levels of these proteins decreased as people got older, as did the amount of proteins produced. To determine which region of the intervertebral disc different proteins were in, Tam, Chen et al. also performed magnetic resonance imaging (MRI) of the samples to correlate image intensity (which represents water content) with the corresponding protein signature. The data obtained provides a high-quality map of how the location of different proteins changes with age, and is available online under the name DIPPER. This database is an informative resource for research into skeletal biology, and it will likely advance the understanding of intervertebral disc degeneration in humans and animals, potentially leading to the development of new treatment strategies for this condition.

Lgr5 and Col22a1 Mark Progenitor Cells in the Lineage toward Juvenile Articular Chondrocytes.

The synovial joint forms from a pool of progenitor cells in the future region of the joint, the interzone. Expression of Gdf5 and Wnt9a has been used to mark the earliest cellular processes in the formation of the interzone and the progenitor cells. However, lineage specification and progression toward the different tissues of the joint are not well understood. Here, by lineage-tracing studies we identify a population of Lgr5+ interzone cells that contribute to the formation of cruciate ligaments, synovial membrane, and articular chondrocytes of the joint. This finding is supported by single-cell transcriptome analyses. We show that Col22a1, a marker of early articular chondrocytes, is co-expressed with Lgr5+ cells prior to cavitation as an important lineage marker specifying the progression toward articular chondrocytes. Lgr5+ cells contribute to the repair of a joint defect with the re-establishment of a Col22a1-expressing superficial layer.

Identification of Genes Associated With Hirschsprung Disease, Based on Whole-Genome Sequence Analysis, and Potential Effects on Enteric Nervous System Development.

BACKGROUND & AIMS: Hirschsprung disease, or congenital aganglionosis, is believed to be oligogenic-that is, caused by multiple genetic factors. We performed whole-genome sequence analyses of patients with Hirschsprung disease to identify genetic factors that contribute to disease development and analyzed the functional effects of these variants. METHODS: We performed whole-genome sequence analyses of 443 patients with short-segment disease, recruited from hospitals in China and Vietnam, and 493 ethnically matched individuals without Hirschsprung disease (controls). We performed genome-wide association analyses and gene-based rare-variant burden tests to identify rare and common disease-associated variants and study their interactions. We obtained induced pluripotent stem cell (iPSC) lines from 4 patients with Hirschsprung disease and 2 control individuals, and we used these to generate enteric neural crest cells for transcriptomic analyses. We assessed the neuronal lineage differentiation capability of iPSC-derived enteric neural crest cells using an in vitro differentiation assay. RESULTS: We identified 4 susceptibility loci, including 1 in the phospholipase D1 gene (PLD1) (P = 7.4 × 10-7). The patients had a significant excess of rare protein-altering variants in genes previously associated with Hirschsprung disease and in the ß-secretase 2 gene (BACE2) (P = 2.9 × 10-6). The epistatic effects of common and rare variants across these loci provided a sensitized background that increased risk for the disease. In studies of the iPSCs, we observed common and distinct pathways associated with variants in RET that affect risk. In functional assays, we found variants in BACE2 to protect enteric neurons from apoptosis. We propose that alterations in BACE1 signaling via amyloid ß precursor protein and BACE2 contribute to pathogenesis of Hirschsprung disease. CONCLUSIONS: In whole-genome sequence analyses of patients with Hirschsprung disease, we identified rare and common variants associated with disease risk. Using iPSC cells, we discovered some functional effects of these variants.

Trans-Ethnic Polygenic Analysis Supports Genetic Overlaps of Lumbar Disc Degeneration With Height, Body Mass Index, and Bone Mineral Density.

Lumbar disc degeneration (LDD) is age-related break-down in the fibrocartilaginous joints between lumbar vertebrae. It is a major cause of low back pain and is conventionally assessed by magnetic resonance imaging (MRI). Like most other complex traits, LDD is likely polygenic and influenced by both genetic and environmental factors. However, genome-wide association studies (GWASs) of LDD have uncovered few susceptibility loci due to the limited sample size. Previous epidemiology studies of LDD also reported multiple heritable risk factors, including height, body mass index (BMI), bone mineral density (BMD), lipid levels, etc. Genetics can help elucidate causality between traits and suggest loci with pleiotropic effects. One such approach is polygenic score (PGS) which summarizes the effect of multiple variants by the summation of alleles weighted by estimated effects from GWAS. To investigate genetic overlaps of LDD and related heritable risk factors, we calculated the PGS of height, BMI, BMD and lipid levels in a Chinese population-based cohort with spine MRI examination and a Japanese case-control cohort of lumbar disc herniation (LDH) requiring surgery. Because most large-scale GWASs were done in European populations, PGS of corresponding traits were created using weights from European GWASs. We calibrated their prediction performance in independent Chinese samples, then tested associations with MRI-derived LDD scores and LDH affection status. The PGS of height, BMI, BMD and lipid levels were strongly associated with respective phenotypes in Chinese, but phenotype variances explained were lower than in Europeans which would reduce the power to detect genetic overlaps. Despite of this, the PGS of BMI and lumbar spine BMD were significantly associated with LDD scores; and the PGS of height was associated with the increased the liability of LDH. Furthermore, linkage disequilibrium score regression suggested that, osteoarthritis, another degenerative disorder that shares common features with LDD, also showed genetic correlations with height, BMI and BMD. The findings suggest a common key contribution of biomechanical stress to the pathogenesis of LDD and will direct the future search for pleiotropic genes.

Inhibiting the integrated stress response pathway prevents aberrant chondrocyte differentiation thereby alleviating chondrodysplasia.

The integrated stress response (ISR) is activated by diverse forms of cellular stress, including endoplasmic reticulum (ER) stress, and is associated with diseases. However, the molecular mechanism(s) whereby the ISR impacts on differentiation is incompletely understood. Here, we exploited a mouse model of Metaphyseal Chondrodysplasia type Schmid (MCDS) to provide insight into the impact of the ISR on cell fate. We show the protein kinase RNA-like ER kinase (PERK) pathway that mediates preferential synthesis of ATF4 and CHOP, dominates in causing dysplasia by reverting chondrocyte differentiation via ATF4-directed transactivation of Sox9. Chondrocyte survival is enabled, cell autonomously, by CHOP and dual CHOP-ATF4 transactivation of Fgf21. Treatment of mutant mice with a chemical inhibitor of PERK signaling prevents the differentiation defects and ameliorates chondrodysplasia. By preventing aberrant differentiation, titrated inhibition of the ISR emerges as a rationale therapeutic strategy for stress-induced skeletal disorders.

Activating the unfolded protein response in osteocytes causes hyperostosis consistent with craniodiaphyseal dysplasia.

Bone remodeling is a balanced process between bone synthesis and degradation, maintaining homeostasis and a constant bone mass in adult life. Imbalance will lead to conditions such as osteoporosis or hyperostosis. Osteoblasts build bone, becoming embedded in bone matrix as mature osteocytes. Osteocytes have a role in sensing and translating mechanical loads into biochemical signals, regulating the differentiation and activity of osteoblasts residing at the bone surface through the secretion of Sclerostin (SOST), an inhibitor of WNT signaling. Excessive mechanical load can lead to activation of cellular stress responses altering cell behavior and differentiation. The unfolded protein response (UPR) is a shared pathway utilized by cells to cope with stress stimuli. We showed that in a transgenic mouse model, activation of the UPR in early differentiating osteocytes delays maturation, maintaining active bone synthesis. In addition, expression of SOST is delayed or suppressed; resulting in active WNT signaling and enhanced periosteal bone formation, and the combined outcome is generalized hyperostosis. A clear relationship between the activation of the unfolded protein response was established and the onset of hyperostosis that can be suppressed with a chemical chaperone, sodium 4-phenobutyrate (4-PBA). As the phenotype is highly consistent with craniodiaphyseal dysplasia (CDD; OMIM 122860), we propose activation of the UPR could be part of the disease mechanism for CDD patients as these patients are heterozygous for SOST mutations that impair protein folding and secretion. Thus, therapeutic agents ameliorating protein folding or the UPR can be considered as a potential therapeutic treatment.