AP-1 family members act with Sox9 to promote chondrocyte hypertrophy

The role of NF-κB-SOX9 signalling pathway in osteoarthritis

The nuclear factor-κB (NF-κB) signalling pathway exists in a variety of cells and is involved in the gene regulation of various physiological and pathological processes such as inflammation, immunity, cell proliferation and apoptosis. It has been shown that this signaling pathway is also involved in numerous events associated with osteoarthritis, including chondrocyte catabolism, chondrocyte survival, and synovial inflammation. SRY-related high mobility group-box 9(SOX9) is the "master regulator" of chondrocytes and one of the key transcription factors that maintain chondrocyte phenotype and cartilage homeostasis. NF-κB can positively regulate the expression of SOX9 by directly binding to its promoter region, and play a role in the formation and development of chondrocytes. This article reviews the regulatory effect of the NF-κB-SOX9 signaling axis on osteoarthritis.

Transfer learning reveals sequence determinants of the quantitative response to transcription factor dosage

Deep learning approaches have made significant advances in predicting cell type-specific chromatin patterns from the identity and arrangement of transcription factor (TF) binding motifs. However, most models have been applied in unperturbed contexts, precluding a predictive understanding of how chromatin state responds to TF perturbation. Here, we used transfer learning to train and interpret deep learning models that use DNA sequence to predict, with accuracy approaching experimental reproducibility, how the concentration of two dosage-sensitive TFs (TWIST1, SOX9) affects regulatory element (RE) chromatin accessibility in facial progenitor cells. High-affinity motifs that allow for heterotypic TF co-binding and are concentrated at the center of REs buffer against quantitative changes in TF dosage and strongly predict unperturbed accessibility. In contrast, motifs with low-affinity or homotypic binding distributed throughout REs lead to sensitive responses with minimal contributions to unperturbed accessibility. Both buffering and sensitizing features show signatures of purifying selection. We validated these predictive sequence features using reporter assays and showed that a biophysical model of TF-nucleosome competition can explain the sensitizing effect of low-affinity motifs. Our approach of combining transfer learning and quantitative measurements of the chromatin response to TF dosage therefore represents a powerful method to reveal additional layers of the cis-regulatory code.

Runt-related Transcription Factors and Gene Regulatory Mechanisms in Skeletal Development and Diseases

PURPOSE OF REVIEW

Runt-related transcription factors (RUNX) play critical roles in skeletal development, metabolism, and diseases. In mammals, three RUNX members, namely RUNX1, RUNX2, and RUNX3, play distinct and redundant roles, although RUNX2 is a dominant factor in skeletal development and several skeletal diseases. This review is to provide an overview of the current understanding of RUNX-mediated transcriptional regulation in different skeletal cell types.

RECENT FINDINGS

Advances in chromatin immunoprecipitation and next-generation sequencing (ChIP-seq) have revealed genome-wide RUNX-mediated gene regulatory mechanisms, including their association with cis-regulatory elements and putative target genes. Further studies with genome-wide analysis and biochemical assays have shed light on RUNX-mediated pioneering action and involvements of RUNX2 in lipid-lipid phase separation. Emerging multi-layered mechanisms of RUNX-mediated gene regulations help us better understanding of skeletal development and diseases, which also provides clues to think how genome-wide studies can help develop therapeutic strategies for skeletal diseases.

Cartilage tissue from sites of weight bearing in patients with osteoarthritis exhibits a differential phenotype with distinct chondrocytes subests

OBJECTIVE

Osteoarthritis (OA) is a degenerative joint disease associated with excessive mechanical loading. The aim here was to elucidate whether different subpopulations of chondrocytes exhibit distinct phenotypes in response to variations in loading conditions. Furthermore, we seek to investigate the transcriptional switches and cell crosstalk among these chondrocytes subsets.

METHODS

Proteomic analysis was performed on cartilage tissues isolated from weight-bearing and non-weight-bearing regions. Additionally, single-cell RNA sequencing was employed to identify different subsets of chondrocytes. For disease-specific cells, in vitro differentiation induction was performed, and their presence was confirmed in human cartilage tissue sections using immunofluorescence. The molecular mechanisms underlying transcriptional changes in these cells were analysed through whole-transcriptome sequencing.

RESULTS

In the weight-bearing regions of OA cartilage tissue, a subpopulation of chondrocytes called OA hypertrophic chondrocytes (OAHCs) expressing the marker genes SLC39A14 and COL10A1 are present. These cells exhibit unique characteristics of active cellular interactions mediated by the TGFβ signalling pathway and express OA phenotypes, distinct from hypertrophic chondrocytes in healthy cartilage. OAHCs are mainly distributed in the superficial region of damaged cartilage in human OA tissue, and on TGFβ stimulation, exhibit activation of transcriptional expression of iron metabolism-related genes, along with enrichment of associated pathways.

CONCLUSION

This study identified and validated the existence of a subset of OAHCs in the weight-bearing area of OA cartilage tissue. Our findings provide a theoretical basis for targeting OAHCs to slow down the progression of OA and facilitate the repair of cartilage injuries.

Metabolic rewiring controlled by c-Fos governs cartilage integrity in osteoarthritis

OBJECTIVES

The activator protein-1 (AP-1) transcription factor component c-Fos regulates chondrocyte proliferation and differentiation, but its involvement in osteoarthritis (OA) has not been functionally assessed.

METHODS

c-Fos expression was evaluated by immunohistochemistry on articular cartilage sections from patients with OA and mice subjected to the destabilisation of the medial meniscus (DMM) model of OA. Cartilage-specific c-Fos knockout (c-FosΔCh) mice were generated by crossing c-fosfl/fl to Col2a1-CreERT mice. Articular cartilage was evaluated by histology, immunohistochemistry, RNA sequencing (RNA-seq), quantitative reverse transcription PCR (qRT-PCR) and in situ metabolic enzyme assays. The effect of dichloroacetic acid (DCA), an inhibitor of pyruvate dehydrogenase kinase (Pdk), was assessed in c-FosΔCh mice subjected to DMM.

RESULTS

FOS-positive chondrocytes were increased in human and murine OA cartilage during disease progression. Compared with c-FosWT mice, c-FosΔCh mice exhibited exacerbated DMM-induced cartilage destruction. Chondrocytes lacking c-Fos proliferate less, have shorter collagen fibres and reduced cartilage matrix. Comparative RNA-seq revealed a prominent anaerobic glycolysis gene expression signature. Consistently decreased pyruvate dehydrogenase (Pdh) and elevated lactate dehydrogenase (Ldh) enzymatic activities were measured in situ, which are likely due to higher expression of hypoxia-inducible factor-1α, Ldha, and Pdk1 in chondrocytes. In vivo treatment of c-FosΔCh mice with DCA restored Pdh/Ldh activity, chondrocyte proliferation, collagen biosynthesis and decreased cartilage damage after DMM, thereby reverting the deleterious effects of c-Fos inactivation.

CONCLUSIONS

c-Fos modulates cellular bioenergetics in chondrocytes by balancing pyruvate flux between anaerobic glycolysis and the tricarboxylic acid cycle in response to OA signals. We identify a novel metabolic adaptation of chondrocytes controlled by c-Fos-containing AP-1 dimers that could be therapeutically relevant.

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.

Hypomorphic and dominant-negative impact of truncated SOX9 dysregulates Hedgehog-Wnt signaling, causing campomelia

Haploinsufficiency for SOX9, the master chondrogenesis transcription factor, can underlie campomelic dysplasia (CD), an autosomal dominant skeletal malformation syndrome, because heterozygous Sox9 null mice recapitulate the bent limb (campomelia) and some other phenotypes associated with CD. However, in vitro cell assays suggest haploinsufficiency may not apply for certain mutations, notably those that truncate the protein, but in these cases in vivo evidence is lacking and underlying mechanisms are unknown. Here, using conditional mouse mutants, we compared the impact of a heterozygous Sox9 null mutation (Sox9+/-) with the Sox9+/Y440X CD mutation that truncates the C-terminal transactivation domain but spares the DNA-binding domain. While some Sox9+/Y440X mice survived, all Sox9+/- mice died perinatally. However, the skeletal defects were more severe and IHH signaling in developing limb cartilage was significantly enhanced in Sox9+/Y440X compared with Sox9+/-. Activating Sox9Y440X specifically in the chondrocyte-osteoblast lineage caused milder campomelia, and revealed cell- and noncell autonomous mechanisms acting on chondrocyte differentiation and osteogenesis in the perichondrium. Transcriptome analyses of developing Sox9+/Y440X limbs revealed dysregulated expression of genes for the extracellular matrix, as well as changes consistent with aberrant WNT and HH signaling. SOX9Y440X failed to interact with β-catenin and was unable to suppress transactivation of Ihh in cell-based assays. We propose enhanced HH signaling in the adjacent perichondrium induces asymmetrically localized excessive perichondrial osteogenesis resulting in campomelia. Our study implicates combined haploinsufficiency/hypomorphic and dominant-negative actions of SOX9Y440X, cell-autonomous and noncell autonomous mechanisms, and dysregulated WNT and HH signaling, as the cause of human campomelia.

SOX9 in organogenesis: shared and unique transcriptional functions

The transcription factor SOX9 is essential for the development of multiple organs including bone, testis, heart, lung, pancreas, intestine and nervous system. Mutations in the human SOX9 gene led to campomelic dysplasia, a haploinsufficiency disorder with several skeletal malformations frequently accompanied by 46, XY sex reversal. The mechanisms underlying the diverse SOX9 functions during organ development including its post-translational modifications, the availability of binding partners, and tissue-specific accessibility to target gene chromatin. Here we summarize the expression, activities, and downstream target genes of SOX9 in molecular genetic pathways essential for organ development, maintenance, and function. We also provide an insight into understanding the mechanisms that regulate the versatile roles of SOX9 in different organs.

Runx2 regulates chromatin accessibility to direct the osteoblast program at neonatal stages

The transcriptional regulator Runx2 (runt-related transcription factor 2) has essential but distinct roles in osteoblasts and chondrocytes in skeletal development. However, Runx2-mediated regulatory mechanisms underlying the distinctive programming of osteoblasts and chondrocytes are not well understood. Here, we perform an integrative analysis to investigate Runx2-DNA binding and chromatin accessibility ex vivo using neonatal osteoblasts and chondrocytes. We find that Runx2 engages with cell-type-distinct chromatin-accessible regions, potentially interacting with different combinations of transcriptional regulators, forming cell-type-specific hotspots, and potentiating chromatin accessibility. Genetic analysis and direct cellular reprogramming studies suggest that Runx2 is essential for establishment of chromatin accessibility in osteoblasts. Functional enhancer studies identify an Sp7 distal enhancer driven by Runx2-dependent binding and osteoblast-specific chromatin accessibility, contributing to normal osteoblast differentiation. Our findings provide a framework for understanding the regulatory landscape encompassing Runx2-mediated and cell-type-distinct enhancer networks that underlie the specification of osteoblasts.

Hojo et al. investigate the gene-regulatory landscape underlying specification of skeletal cell types in neonatal mice. Runx2, an osteoblast determinant, engages with cell-type-distinct chromatin-accessible regions and is essential for establishment of chromatin accessibility in osteoblasts. The study provides insights into enhancer networks in skeletal development.

Limb Mesoderm and Head Ectomesenchyme Both Express a Core Transcriptional Program During Chondrocyte Differentiation

To explain how cartilage appeared in different parts of the vertebrate body at discrete times during evolution, we hypothesize that different embryonic populations co-opted expression of a core gene regulatory network (GRN) driving chondrocyte differentiation. To test this hypothesis, laser-capture microdissection coupled with RNA-seq was used to reveal chondrocyte transcriptomes in the developing chick humerus and ceratobranchial, which are mesoderm- and neural crest-derived, respectively. During endochondral ossification, two general types of chondrocytes differentiate. Immature chondrocytes (IMM) represent the early stages of cartilage differentiation, while mature chondrocytes (MAT) undergo additional stages of differentiation, including hypertrophy and stimulating matrix mineralization and degradation. Venn diagram analyses generally revealed a high degree of conservation between chondrocyte transcriptomes of the limb and head, including SOX9, COL2A1, and ACAN expression. Typical maturation genes, such as COL10A1, IBSP, and SPP1, were upregulated in MAT compared to IMM in both limb and head chondrocytes. Gene co-expression network (GCN) analyses of limb and head chondrocyte transcriptomes estimated the core GRN governing cartilage differentiation. Two discrete portions of the GCN contained genes that were differentially expressed in limb or head chondrocytes, but these genes were enriched for biological processes related to limb/forelimb morphogenesis or neural crest-dependent processes, respectively, perhaps simply reflecting the embryonic origin of the cells. A core GRN driving cartilage differentiation in limb and head was revealed that included typical chondrocyte differentiation and maturation markers, as well as putative novel "chondrocyte" genes. Conservation of a core transcriptional program during chondrocyte differentiation in both the limb and head suggest that the same core GRN was co-opted when cartilage appeared in different regions of the skeleton during vertebrate evolution.

Identification of Transcription Factors Responsible for a Transforming Growth Factor-β-Driven Hypertrophy-like Phenotype in Human Osteoarthritic Chondrocytes

During osteoarthritis (OA), hypertrophy-like chondrocytes contribute to the disease process. TGF-β's signaling pathways can contribute to a hypertrophy(-like) phenotype in chondrocytes, especially at high doses of TGF-β. In this study, we examine which transcription factors (TFs) are activated and involved in TGF-β-dependent induction of a hypertrophy-like phenotype in human OA chondrocytes. We found that TGF-β, at levels found in synovial fluid in OA patients, induces hypertrophic differentiation, as characterized by increased expression of RUNX2, COL10A1, COL1A1, VEGFA and IHH. Using luciferase-based TF activity assays, we observed that the expression of these hypertrophy genes positively correlated to SMAD3:4, STAT3 and AP1 activity. Blocking these TFs using specific inhibitors for ALK-5-induced SMAD signaling (5 µM SB-505124), JAK-STAT signaling (1 µM Tofacitinib) and JNK signaling (10 µM SP-600125) led to the striking observation that only SB-505124 repressed the expression of hypertrophy factors in TGF-β-stimulated chondrocytes. Therefore, we conclude that ALK5 kinase activity is essential for TGF-β-induced expression of crucial hypertrophy factors in chondrocytes.

Cell-intrinsic Aryl Hydrocarbon Receptor signalling is required for the resolution of injury-induced colonic stem cells

The aryl hydrocarbon receptor (AHR) is an environmental sensor that integrates microbial and dietary cues to influence physiological processes within the intestinal microenvironment, protecting against colitis and colitis-associated colorectal cancer development. Rapid tissue regeneration upon injury is important for the reinstatement of barrier integrity and its dysregulation promotes malignant transformation. Here we show that AHR is important for the termination of the regenerative response and the reacquisition of mature epithelial cell identity post injury in vivo and in organoid cultures in vitro. Using an integrative multi-omics approach in colon organoids, we show that AHR is required for timely termination of the regenerative response through direct regulation of transcription factors involved in epithelial cell differentiation as well as restriction of chromatin accessibility to regeneration-associated Yap/Tead transcriptional targets. Safeguarding a regulated regenerative response places AHR at a pivotal position in the delicate balance between controlled regeneration and malignant transformation.

Rapid intestinal regeneration after injury is critical to maintain barrier integrity and homeostasis, but must be tightly controlled to prevent tumorigenesis. Here they show that the aryl hydrocarbon receptor is required to terminate the regenerative response after wound healing.

The Evolution of Biomineralization through the Co-Option of Organic Scaffold Forming Networks

Biomineralization is the process in which organisms use minerals to generate hard structures like teeth, skeletons and shells. Biomineralization is proposed to have evolved independently in different phyla through the co-option of pre-existing developmental programs. Comparing the gene regulatory networks (GRNs) that drive biomineralization in different species could illuminate the molecular evolution of biomineralization. Skeletogenesis in the sea urchin embryo was extensively studied and the underlying GRN shows high conservation within echinoderms, larval and adult skeletogenesis. The organic scaffold in which the calcite skeletal elements form in echinoderms is a tubular compartment generated by the syncytial skeletogenic cells. This is strictly different than the organic cartilaginous scaffold that vertebrates mineralize with hydroxyapatite to make their bones. Here I compare the GRNs that drive biomineralization and tubulogenesis in echinoderms and in vertebrates. The GRN that drives skeletogenesis in the sea urchin embryo shows little similarity to the GRN that drives bone formation and high resemblance to the GRN that drives vertebrates’ vascular tubulogenesis. On the other hand, vertebrates’ bone-GRNs show high similarity to the GRNs that operate in the cells that generate the cartilage-like tissues of basal chordate and invertebrates that do not produce mineralized tissue. These comparisons suggest that biomineralization in deuterostomes evolved through the phylum specific co-option of GRNs that control distinct organic scaffolds to mineralization.

Genome-scale actions of master regulators directing skeletal development

The mammalian skeleton develops through two distinct modes of ossification: intramembranous ossification and endochondral ossification. During the process of skeletal development, SRY-box containing gene 9 (Sox9), runt-related transcription factor 2 (Runx2), and Sp7 work as master transcription factors (TFs) or transcriptional regulators, underlying cell fate specification of the two distinct populations: bone-forming osteoblasts and cartilage-forming chondrocytes. In the past two decades, core transcriptional circuits underlying skeletal development have been identified mainly through mouse genetics and biochemical approaches. Recently emerging next-generation sequencer (NGS)-based studies have provided genome-scale views on the gene regulatory landscape programmed by the master TFs/transcriptional regulators. With particular focus on Sox9, Runx2, and Sp7, this review aims to discuss the gene regulatory landscape in skeletal development, which has been identified by genome-scale data, and provide future perspectives in this field.

Dmrt2 promotes transition of endochondral bone formation by linking Sox9 and Runx2

Endochondral bone formation is fundamental for skeletal development. During this process, chondrocytes undergo multiple steps of differentiation and coordinated transition from a proliferating to a hypertrophic stage, which is critical to advance skeletal development. Here, we identified the transcription factor Dmrt2 (double-sex and mab-3 related transcription factor 2) as a Sox9-inducible gene that promotes chondrocyte hypertrophy in pre-hypertrophic chondrocytes. Epigenetic analysis further demonstrated that Sox9 regulates Dmrt2 expression through an active enhancer located 18 kb upstream of the Dmrt2 gene and that this enhancer's chromatin status is progressively activated through chondrocyte differentiation. Dmrt2-knockout mice exhibited a dwarf phenotype with delayed initiation of chondrocyte hypertrophy. Dmrt2 augmented hypertrophic chondrocyte gene expression including Ihh through physical and functional interaction with Runx2. Furthermore, Dmrt2 deficiency reduced Runx2-dependent Ihh expression. Our findings suggest that Dmrt2 is critical for sequential chondrocyte differentiation during endochondral bone formation and coordinates the transcriptional network between Sox9 and Runx2.

SOX9 modulates cancer biomarker and cilia genes in pancreatic cancer

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive form of cancer with high mortality. The cellular origins of PDAC are largely unknown; however, ductal cells, especially centroacinar cells (CACs), have several characteristics in common with PDAC, such as expression of SOX9 and components of the Notch-signaling pathway. Mutations in KRAS and alterations to Notch signaling are common in PDAC, and both these pathways regulate the transcription factor SOX9. To identify genes regulated by SOX9, we performed siRNA knockdown of SOX9 followed by RNA-seq in PANC-1s, a human PDAC cell line. We report 93 differentially expressed (DE) genes, with convergence on alterations to Notch-signaling pathways and ciliogenesis. These results point to SOX9 and Notch activity being in a positive feedback loop and SOX9 regulating cilia production in PDAC. We additionally performed ChIP-seq in PANC-1s to identify direct targets of SOX9 binding and integrated these results with our DE gene list. Nine of the top 10 downregulated genes have evidence of direct SOX9 binding at their promoter regions. One of these targets was the cancer stem cell marker EpCAM. Using whole-mount in situ hybridization to detect epcam transcript in zebrafish larvae, we demonstrated that epcam is a CAC marker and that Sox9 regulation of epcam expression is conserved in zebrafish. Additionally, we generated an epcam null mutant and observed pronounced defects in ciliogenesis during development. Our results provide a link between SOX9, EpCAM and ciliary repression that can be exploited in improving our understanding of the cellular origins and mechanisms of PDAC.

Constitutive activation of mitogen-activated protein kinase kinase (MEK1) in ileal enterocytes leads to dysplasia and a predisposition to cancer

Activation of mitogen-activated protein kinases (MAPKs) is a key factor in the pathogenesis of cancer, although the specific role of mitogen-activated protein kinase kinase (MEK1) is not well understood. Villin promoter-driven Cre expression was used to excise a floxed stop cassette from a phosphomimetically constitutively activated MEK1 (caMEK1) expression construct in the intestine of C57BL/6 mice. Zygosity status of caMEK1 afforded assessment of the dose dependence of the effect. The expected mendelian distribution of genotypes and sex was observed in 443 progenies. Between 21 and 63 days of life, caMEK1 had no effect on body weight in male mice, but reduced body weight in female mice homozygous for caMEK1. At 10 wk of age, the ileum of caMEK1-expressing mice was characterized by the finding of dysplasia and profound changes in overall architecture. Paneth cells were nearly absent in caMEK1 homozygotes. Targeted proteomic profiling via reverse phase protein array analyses with confirmatory Western blotting revealed significant changes in protein and phosphoprotein expression, including upregulation of proteins downstream of MEK1, associated with enhanced markers of proliferation, diminished apoptosis, alterations in cell-fate determination, cell-cell interactions, and tight junctions. Long-term viability of caMEK1 homozygous mice was reduced with no survival beyond 1 yr. Invasive adenocarcinoma developed in three of ten older mice [15 wk (homozygous), 26 wk (homozygous), and 35 wk (heterozygous) of age]. Expression of caMEK1 in enterocytes leads to marked derangements in the intestinal epithelium, which is associated with a predisposition to the development of invasive cancer.NEW & NOTEWORTHY The ileum of mice with constitutive expression of activated MEK1 (via phosphomimetic changes) in enterocytes is markedly abnormal with architectural distortion and cytologic atypia, which evolves into an adenoma invasive carcinoma sequence. Phosphoproteomic analysis reveals upregulation of proteins downstream of MEK1, associated with enhanced markers of proliferation, diminished apoptosis, alterations in cell-fate determination, cell-cell interactions, and tight junctions. This novel model provides new insights into intestinal homeostasis and carcinogenesis.

Ultra-processed food targets bone quality via endochondral ossification

Ultra-processed foods have known negative implications for health; however, their effect on skeletal development has never been explored. Here, we show that young rats fed ultra-processed food rich in fat and sugar suffer from growth retardation due to lesions in their tibial growth plates. The bone mineral density decreases significantly, and the structural parameters of the bone deteriorate, presenting a sieve-like appearance in the cortices and poor trabecular parameters in long bones and vertebrae. This results in inferior mechanical performance of the entire bone with a high fracture risk. RNA sequence analysis of the growth plates demonstrated an imbalance in extracellular matrix formation and degradation and impairment of proliferation, differentiation and mineralization processes. Our findings highlight, for the first time, the severe impact of consuming ultra-processed foods on the growing skeleton. This pathology extends far beyond that explained by the known metabolic effects, highlighting bone as a new target for studies of modern diets.

Screening and verification of hub genes involved in osteoarthritis using bioinformatics

Osteoarthritis (OA) is one of the most common causes of disability and its development is associated with numerous factors. A major challenge in the treatment of OA is the lack of early diagnosis. In the present study, a bioinformatics method was employed to filter key genes that may be responsible for the pathogenesis of OA. From the Gene Expression Omnibus database, the datasets GSE55457, GSE12021 and GSE55325 were downloaded, which comprised 59 samples. Of these, 30 samples were from patients diagnosed with osteoarthritis and 29 were normal. Differentially expressed genes (DEGs) were obtained by downloading and analyzing the original data using bioinformatics. The Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathways were analyzed using the Database for Annotation, Visualization and Integrated Discovery online database. Protein-protein interaction network analysis was performed using the Search Tool for the Retrieval of Interacting Genes/proteins online database. BSCL2 lipid droplet biogenesis associated, seipin, FOS-like 2, activator protein-1 transcription factor subunit (FOSL2), cyclin-dependent kinase inhibitor 1A (CDKN1A) and kinectin 1 (KTN1) genes were identified as key genes by using Cytoscape software. Functional enrichment revealed that the DEGs were mainly accumulated in the ErbB, MAPK and PI3K-Akt pathways. Reverse transcription-quantitative PCR analysis confirmed a significant reduction in the expression levels of FOSL2, CDKN1A and KTN1 in OA samples. These genes have the potential to become novel diagnostic and therapeutic targets for OA.

ChIP-Seq Assays from Mammalian Cartilage and Chondrocytes

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is a powerful tool to identify binding profiles of transcriptional regulators and chromatin regulators as well as histone modification patterns in a genome-wide manner. ChIP-seq consists of five major steps: (1) preparation of cells and chromatin, (2) ChIP, (3) ChIP-seq library construction, (4) sequencing of ChIP DNA with a next-generation sequencer (NGS), and (5) computational analysis of sequence data. Recent ChIP-seq studies in skeletal tissues enable us to understand the modes of action of key skeletal regulators, functional interaction among the enhancers bound by the regulators, the complex nature of regulatory inputs, and thereby the gene regulatory landscape in skeletal development. Here we describe a ChIP-seq protocol that we have employed in our studies, with particular focus on chromatin preparation and subsequent ChIP in skeletal cells, including chondrocytes.

Transcription Factors in Cartilage Homeostasis and Osteoarthritis

Osteoarthritis (OA) is the most common degenerative joint disease, and it is characterized by articular cartilage loss. In part, OA is caused by aberrant anabolic and catabolic activities of the chondrocyte, the only cell type present in cartilage. These chondrocyte activities depend on the intra- and extracellular signals that the cell receives and integrates into gene expression. The key proteins for this integration are transcription factors. A large number of transcription factors exist, and a better understanding of the transcription factors activated by the various signaling pathways active during OA can help us to better understand the complex etiology of OA. In addition, establishing such a profile can help to stratify patients in different subtypes, which can be a very useful approach towards personalized therapy. In this review, we discuss crucial transcription factors for extracellular matrix metabolism, chondrocyte hypertrophy, chondrocyte senescence, and autophagy in chondrocytes. In addition, we discuss how insight into these factors can be used for treatment purposes.

Extracellular Matrix and Cellular Plasticity in Musculoskeletal Development

Cellular plasticity refers to the ability of cell fates to be reprogrammed given the proper signals, allowing for dedifferentiation or transdifferentiation into different cell fates. In vitro, this can be induced through direct activation of gene expression, however this process does not naturally occur in vivo. Instead, the microenvironment consisting of the extracellular matrix (ECM) and signaling factors, directs the signals presented to cells. Often the ECM is involved in regulating both biochemical and mechanical signals. In stem cell populations, this niche is necessary for maintenance and proper function of the stem cell pool. However, recent studies have demonstrated that differentiated or lineage restricted cells can exit their current state and transform into another state under different situations during development and regeneration. This may be achieved through (1) cells responding to a changing niche; (2) cells migrating and encountering a new niche; and (3) formation of a transitional niche followed by restoration of the homeostatic niche to sequentially guide cells along the regenerative process. This review focuses on examples in musculoskeletal biology, with the concept of ECM regulating cells and stem cells in development and regeneration, extending beyond the conventional concept of small population of progenitor cells, but under the right circumstances even “lineage-restricted” or differentiated cells can be reprogrammed to enter into a different fate.

Regulation and function of SOX9 during cartilage development and regeneration

Chondrogenesis is a highly coordinated event in embryo development, adult homeostasis, and repair of the vertebrate cartilage. Fate decisions and differentiation of chondrocytes accompany differential expression of genes critical for each step of chondrogenesis. SOX9 is a master transcription factor that participates in sequential events in chondrogenesis by regulating a series of downstream factors in a stage-specific manner. SOX9 either works alone or in combination with downstream SOX transcription factors, SOX5 and SOX6 as chondrogenic SOX Trio. SOX9 is reduced in the articular cartilage of patients with osteoarthritis while highly maintained during tumorigenesis of cartilage and bone. Gene therapy using viral and non-viral vectors accompanied by tissue engineering (scaffolds) is a promising tool to regenerate impaired cartilage. Delivery of SOX9 or chondrogenic SOX Trio into cells produces efficient therapeutic effects on chondrogenesis and this event is facilitated by scaffolds. Non-viral vector-guided delivery systems encapsulated or loaded in mechanically stable solid scaffolds are useful for the regeneration of articular cartilage. Here we review major milestones and most recent studies focusing on regulation and function of chondrogenic SOX Trio, during chondrogenesis and cartilage regeneration, and on the development of advanced technologies in gene delivery with tissue engineering to improve efficiency of cartilage repair process.

Histone ChIP-Seq identifies differential enhancer usage during chondrogenesis as critical for defining cell-type specificity

Epigenetic mechanisms are known to regulate gene expression during chondrogenesis. In this study, we have characterized the epigenome during the in vitro differentiation of human mesenchymal stem cells (hMSCs) into chondrocytes. Chromatin immunoprecipitation followed by next‐generation sequencing (ChIP‐seq) was used to assess a range of N‐terminal posttranscriptional modifications (marks) to histone H3 lysines (H3K4me3, H3K4me1, H3K27ac, H3K27me3, and H3K36me3) in both hMSCs and differentiated chondrocytes. Chromatin states were characterized using histone ChIP‐seq and cis‐regulatory elements were identified in chondrocytes. Chondrocyte enhancers were associated with chondrogenesis‐related gene ontology (GO) terms. In silico analysis and integration of DNA methylation data with chondrogenesis chromatin states revealed that enhancers marked by histone marks H3K4me1 and H3K27ac were de‐methylated during in vitro chondrogenesis. Similarity analysis between hMSC and chondrocyte chromatin states defined in this study with epigenomes of cell‐types defined by the Roadmap Epigenomics project revealed that enhancers are more distinct between cell‐types compared to other chromatin states. Motif analysis revealed that the transcription factor SOX9 is enriched in chondrocyte enhancers. Luciferase reporter assays confirmed that chondrocyte enhancers characterized in this study exhibited enhancer activity which may be modulated by DNA methylation and SOX9 overexpression. Altogether, these integrated data illustrate the cross‐talk between different epigenetic mechanisms during chondrocyte differentiation.

Sox100B Regulates Progenitor-Specific Gene Expression and Cell Differentiation in the Adult Drosophila Intestine

Robust production of terminally differentiated cells from self-renewing resident stem cells is essential to maintain proper tissue architecture and physiological functions, especially in high-turnover tissues. However, the transcriptional networks that precisely regulate cell transition and differentiation are poorly understood in most tissues. Here, we identified Sox100B, a Drosophila Sox E family transcription factor, as a critical regulator of adult intestinal stem cell differentiation. Sox100B is expressed in stem and progenitor cells and required for differentiation of enteroblast progenitors into absorptive enterocytes. Mechanistically, Sox100B regulates the expression of another critical stem cell differentiation factor, Sox21a. Supporting a direct control of Sox21a by Sox100B, we identified a Sox21a intronic enhancer that is active in all intestinal progenitors and directly regulated by Sox100B. Taken together, our results demonstrate that the activity and regulation of two Sox transcription factors are essential to coordinate stem cell differentiation and proliferation and maintain intestinal tissue homeostasis.

•

Sox100B is expressed in progenitor cells in the adult intestine

•

Sox100B is required for stem cell differentiation

•

Sox100B is required for Sox21a expression

•

Sox100B directly controls the activity of a Sox21a intronic enhancer

Transcriptional Regulation of Vih by Oct4 and Sox9 in Scylla paramamosain

Mud crab (Scylla paramamosain) is one of the most economically-important marine crabs in China. However, research on mechanisms of reproductive regulation is not sufficient. Vitellogenesis-inhibiting hormone (VIH) is a member of the crustacean hyperglycemia hormones (CHH) family, which plays an essential role in the regulation of gonadal development and maturation in crustaceans, and current studies on the regulation of Vih transcription in crabs are relatively rare. Our previous studies on the transcriptional regulation of mud crab Vih (SpVih) have proved that the binding site of Oct4/Sox9 transcription factor may be the key region for positively regulating the expression of SpVih. In this study, the electrophoretic mobility shift assay (EMSA) experiment confirmed that the nuclear protein extracted from the eyestalk could bind to the key region of SpVih promoter, and these specific bindings were dependent on the presence of Oct4/Sox9 binding sites. Two specific binding complex bands were detected in the supershift group of EMSA supershift experiments by Oct4 and Sox9 antibodies, further confirming the specific recognition of these two transcription factors on the key regulatory region of SpVih. In vitro, Oct4 and Sox9 gene overexpression vectors and SpVih core promoter fragment vector were constructed and co-transfected into HEK293T cells. As a result, SpVih activity increased with the concentration of transcription factors. In vivo, when Oct4 and Sox9 dsRNA were injected into the eyestalks of mud crab, respectively, the expression level of SpVih decreased significantly after interference with Oct4 or Sox9, and the expression level of SpVtg in the ovary and hepatopancreatic increased. Both in vitro and in vivo experiments showed that Oct4 and Sox9 had a positive regulatory effect on SpVih. The GST pull-down experiment was carried out by purified Oct4 and Sox9 proteins, and the results showed that there was an interaction between them. It was speculated that they regulated the expression of SpVih through the interaction.

Insights into Gene Regulatory Networks in Chondrocytes

Chondrogenesis is a key developmental process that molds the framework of our body and generates the skeletal tissues by coupling with osteogenesis. The developmental processes are well-coordinated by spatiotemporal gene expressions, which are hardwired with gene regulatory elements. Those elements exist as thousands of modules of DNA sequences on the genome. Transcription factors function as key regulatory proteins by binding to regulatory elements and recruiting cofactors. Over the past 30 years, extensive attempts have been made to identify gene regulatory mechanisms in chondrogenesis, mainly through biochemical approaches and genetics. More recently, newly developed next-generation sequencers (NGS) have identified thousands of gene regulatory elements on a genome scale, and provided novel insights into the multiple layers of gene regulatory mechanisms, including the modes of actions of transcription factors, post-translational histone modifications, chromatin accessibility, the concept of pioneer factors, and three-dimensional chromatin architecture. In this review, we summarize the studies that have improved our understanding of the gene regulatory mechanisms in chondrogenesis, from the historical studies to the more recent works using NGS. Finally, we consider the future perspectives, including efforts to improve our understanding of the gene regulatory landscape in chondrogenesis and potential applications to the treatment of chondrocyte-related diseases.

SOX9 in cartilage development and disease

SOX9 is a pivotal transcription factor in chondrocytes, a lineage essential in skeletogenesis. Its mandatory role in transactivating many cartilage-specific genes is well established, whereas its pioneer role in lineage specification, which along with transactivation defines master transcription factors, remains to be better defined. Abundant, but yet incomplete evidence exists that intricate molecular networks control SOX9 activity during the multi-step chondrogenesis pathway. They include a highly modular genetic regulation, post-transcriptional and post-translational modifications, and varying sets of functional partners. Fully uncovering SOX9 actions and regulation is fundamental to explain mechanisms underlying many diseases that directly or indirectly affect SOX9 activities and to design effective disease treatments. We here review current knowledge, highlight recent discoveries, and propose new research directions to answer remaining questions.

Roles and regulation of SOX transcription factors in skeletogenesis

SOX transcription factors participate in the specification, differentiation and activities of many cell types in development and beyond. The 20 mammalian family members are distributed into eight groups based on sequence identity, and while co-expressed same-group proteins often have redundant functions, different-group proteins typically have distinct functions. More than a handful of SOX proteins have pivotal roles in skeletogenesis. Heterozygous mutations in their genes cause human diseases, in which skeletal dysmorphism is a major feature, such as campomelic dysplasia (SOX9), or a minor feature, such as LAMSHF syndrome (SOX5) and Coffin-Siris-like syndromes (SOX4 and SOX11). Loss- and gain-of-function experiments in animal models have revealed that SOX4 and SOX11 (SOXC group) promote skeletal progenitor survival and control skeleton patterning and growth; SOX8 (SOXE group) delays the differentiation of osteoblast progenitors; SOX9 (SOXE group) is essential for chondrocyte fate maintenance and differentiation, and works in cooperation with SOX5 and SOX6 (SOXD group) and other types of transcription factors. These and other SOX proteins have also been proposed, mainly through in vitro experiments, to have key roles in other aspects of skeletogenesis, such as SOX2 in osteoblast stem cell self-renewal. We here review current knowledge of well-established and proposed skeletogenic roles of SOX proteins, their transcriptional and non-transcriptional actions, and their modes of regulation at the gene, RNA and protein levels. We also discuss gaps in knowledge and directions for future research to further decipher mechanisms underlying skeletogenesis in health and diseases and identify treatment options for skeletal malformation and degeneration diseases.

Redox control of chondrocyte differentiation and chondrogenesis

Chondrogenesis involves the recruitment and migration of mesenchymal cells, mesenchymal condensation, and chondrocyte differentiation and hypertrophy. Multiple factors precisely regulate chondrogenesis. Recent studies have demonstrated that the redox status of chondrocytes plays an essential role in the regulation of chondrocyte differentiation and chondrogenesis. Reactive oxygen species (ROS) and reactive nitrogen species (RNS) are important factors that change the intracellular redox status. Physiological levels of ROS/RNS act as intracellular signals in chondrocytes, and oxidative stress impairs the metabolism of chondrocytes. Under physiological conditions, the balance between ROS/RNS production and elimination ensures that redox-sensitive signalling proteins function correctly. The redox homeostasis of chondrocytes ensures that they respond appropriately to endogenous and exogenous stimuli. This review focuses on the redox regulation of key signalling pathways and transcription factors that control chondrogenesis and chondrocyte differentiation. Additionally, the mechanism by which ROS/RNS regulate signalling proteins and transcription factors in chondrocytes is also reviewed.

Stem and progenitor cells in skeletal development

Accumulating evidence supports the idea that stem and progenitor cells play important roles in skeletal development. Over the last decade, the definition of skeletal stem and progenitor cells has evolved from cells simply defined by their in vitro behaviors to cells fully defined by a combination of sophisticated approaches, including serial transplantation assays and in vivo lineage-tracing experiments. These approaches have led to better identification of the characteristics of skeletal stem cells residing in multiple sites, including the perichondrium of the fetal bone, the resting zone of the postnatal growth plate, the bone marrow space and the periosteum in adulthood. These diverse groups of skeletal stem cells appear to closely collaborate and achieve a number of important biological functions of bones, including not only bone development and growth, but also bone maintenance and repair. Although these are important findings, we are only beginning to understand the diversity and the nature of skeletal stem and progenitor cells, and how they actually behave in vivo.

Chromatin accessibility landscape of articular knee cartilage reveals aberrant enhancer regulation in osteoarthritis

Osteoarthritis (OA) is a common joint disorder with increasing impact in an aging society. While genetic and transcriptomic analyses have revealed some genes and non-coding loci associated to OA, the pathogenesis remains incompletely understood. Chromatin profiling, which provides insight into gene regulation, has not been reported in OA mainly due to technical difficulties. Here, we employed Assay for Transposase-Accessible Chromatin with high throughput sequencing (ATAC-seq) to map the accessible chromatin landscape in articular knee cartilage of OA patients. We identified 109,215 accessible chromatin regions for cartilages, of which 71% were annotated as enhancers. By overlaying them with genetic and DNA methylation data, we have determined potential OA-relevant enhancers and their putative target genes. Furthermore, through integration with RNA-seq data, we characterized genes that are altered both at epigenomic and transcriptomic levels in OA. These genes are enriched in pathways regulating ossification and mesenchymal stem cell (MSC) differentiation. Consistently, the differentially accessible regions in OA are enriched for MSC-specific enhancers and motifs of transcription factor families involved in osteoblast differentiation. In conclusion, we demonstrate how direct chromatin profiling of clinical tissues can provide comprehensive epigenetic information for a disease and suggest candidate genes and enhancers of translational potential.

TGF-β1-SOX9 axis-inducible COL10A1 promotes invasion and metastasis in gastric cancer via epithelial-to-mesenchymal transition

Molecular biomarkers that predict disease progression might promote drug development and therapeutic strategies in aggressive cancers, such as gastric cancer (GC). High-throughput mRNA sequencing (RNA-seq) revealed that collagen type X alpha 1 (COL10A1) is a disease progression-associated gene. Analysis of 103 GC patients showed that high COL10A1 mRNA expression was associated with GC metastasis and reduced survival. We analyzed the COL10A1 promoter using the UCSC genome website and JASPAR database, and we found potential SOX9 binding site. Here, we demonstrated that SOX9 and COL10A1 were both up-regulated in GC. We observed a positive correlation between the expression patterns of SOX9 and COL10A1 in GC cells and tissues. The results of electrophoretic mobility shift assay (EMSA), chromatin immunoprecipitation (ChIP) assay and promoter reporter indicated that SOX9 could directly bind to the COL10A1 gene promoter and activate its transcription. Biological function experiments showed that COL10A1 regulated the migration and invasion of GC cells. Knockdown COL10A1 inhibited lung and abdominal cavity metastasis in a nude mouse model. Moreover, transforming growth factor-β1 (TGF-β1) treatment up-regulated the phosphorylation of Smad2 and increased SOX9 and COL10A1 expression. COL10A1 was confirmed to be a potential inducer of epithelial-to-mesenchymal transition (EMT). SOX9 was essential for COL10A1-mediated EMT, and cell migration, invasion and metastasis. Co-expression of SOX9 and COL10A1 was associated with tumor progression and was strongly predictive of overall survival in GC patients. In summary, this study elucidated the mechanistic link between COL10A1 and the TGF-β1-SOX9 axis. These findings indicated that COL10A1 might play a crucial role in GC progression and serve as a potential biomarker and therapeutic target in GC patients.

Synergistic co-regulation and competition by a SOX9-GLI-FOXA phasic transcriptional network coordinate chondrocyte differentiation transitions

The growth plate mediates bone growth where SOX9 and GLI factors control chondrocyte proliferation, differentiation and entry into hypertrophy. FOXA factors regulate hypertrophic chondrocyte maturation. How these factors integrate into a Gene Regulatory Network (GRN) controlling these differentiation transitions is incompletely understood. We adopted a genome-wide whole tissue approach to establish a Growth Plate Differential Gene Expression Library (GP-DGEL) for fractionated proliferating, pre-hypertrophic, early and late hypertrophic chondrocytes, as an overarching resource for discovery of pathways and disease candidates. De novo motif discovery revealed the enrichment of SOX9 and GLI binding sites in the genes preferentially expressed in proliferating and prehypertrophic chondrocytes, suggesting the potential cooperation between SOX9 and GLI proteins. We integrated the analyses of the transcriptome, SOX9, GLI1 and GLI3 ChIP-seq datasets, with functional validation by transactivation assays and mouse mutants. We identified new SOX9 targets and showed SOX9-GLI directly and cooperatively regulate many genes such as Trps1, Sox9, Sox5, Sox6, Col2a1, Ptch1, Gli1 and Gli2. Further, FOXA2 competes with SOX9 for the transactivation of target genes. The data support a model of SOX9-GLI-FOXA phasic GRN in chondrocyte development. Together, SOX9-GLI auto-regulate and cooperate to activate and repress genes in proliferating chondrocytes. Upon hypertrophy, FOXA competes with SOX9, and control toward terminal differentiation passes to FOXA, RUNX, AP1 and MEF2 factors.

In the development of the mammalian growth plate, while several transcription factors are individually well known for their key roles in regulating phases of chondrocyte differentiation, there is little information on how they interact and cooperate with each other. We took an unbiased genome wide approach to identify the transcription factors and signaling pathways that play dominant roles in the chondrocyte differentiation cascade. We developed a searchable library of differentially expressed genes, GP-DGEL, which has fine spatial resolution and global transcriptomic coverage for discovery of processes, pathways and disease candidates. Our work identifies a novel regulatory mechanism that integrates the action of three transcription factors, SOX9, GLI and FOXA. SOX9-GLI auto-regulate and cooperate to activate and repress genes in proliferating chondrocytes. Upon entry into prehypertrophy, FOXA competes with SOX9, and control of hypertrophy passes to FOXA, RUNX, AP1 and MEF2 factors.

Transcriptional network systems in cartilage development and disease

Transcription factors play important roles in the regulation of cartilage development by controlling the expression of chondrogenic genes. Genetic studies have revealed that Sox9/Sox5/Sox6, Runx2/Runx3 and Osterix in particular are essential for the sequential steps of cartilage development. Importantly, these transcription factors form network systems that are also required for appropriate cartilage development. Molecular cloning approaches have largely contributed to the identification of several transcriptional partners for Sox9 and Runx2 during cartilage development. Although the importance of a negative-feedback loop between Indian hedgehog (Ihh) and parathyroid hormone-related protein (PTHrP) in chondrocyte hypertrophy has been well established, recent studies indicate that several transcription factors interact with the Ihh-PTHrP loop and demonstrated that Ihh has multiple functions in the regulation of cartilage development. The most common cartilage disorder, osteoarthritis, has been reported to result from the pathological action of several transcription factors, including Runx2, C/EBPβ and HIF-2α. On the other hand, NFAT family members appear to play roles in the protection of cartilage from osteoarthritis. It is also becoming important to understand the homeostasis and regulation of articular chondrocytes, because they have different cellular and molecular features from chondrocytes of the growth plate. This review summarizes the regulation and roles of transcriptional network systems in cartilage development and their pathological roles in osteoarthritis.

Transcriptional Network Controlling Endochondral Ossification

Endochondral ossification is the fundamental process of skeletal development in vertebrates. Chondrocytes undergo sequential steps of differentiation, including mesenchymal condensation, proliferation, hypertrophy, and mineralization. These steps, which are required for the morphological and functional changes in differentiating chondrocytes, are strictly regulated by a complex transcriptional network. Biochemical and mice genetic studies identified chondrogenic transcription factors critical for endochondral ossification. The transcription factor sex-determining region Y (SRY)-box 9 (Sox9) is essential for early chondrogenesis, and impaired Sox9 function causes severe chondrodysplasia in humans and mice. In addition, recent genome-wide chromatin immunoprecipitation-sequencing studies revealed the precise regulatory mechanism of Sox9 during early chondrogenesis. Runt-related transcription factor 2 promotes chondrocyte hypertrophy and terminal differentiation. Interestingly, endoplasmic reticulum (ER) stress-related transcription factors have recently emerged as novel regulators of chondrocyte differentiation. Here we review the transcriptional mechanisms that regulate endochondral ossification, with a focus on Sox9.

Identification of the gene-regulatory landscape in skeletal development and potential links to skeletal regeneration

A class of gene-regulatory elements called enhancers are the main mediators controlling quantitative, temporal and spatial gene expressions. In the course of evolution, the enhancer landscape of higher organisms such as mammals has become quite complex, exerting biological functions precisely and coordinately. In mammalian skeletal development, the master transcription factors Sox9, Runx2 and Sp7/Osterix function primarily through enhancers on the genome to achieve specification and differentiation of skeletal cells. Recently developed genome-scale analyses have shed light on multiple layers of gene regulations, uncovering not only the primary mode of actions of these transcription factors on skeletal enhancers, but also the relation of the epigenetic landscape to three-dimensional chromatin architecture. Here, we review findings on the emerging framework of gene-regulatory networks involved in skeletal development. We further discuss the power of genome-scale analyses to provide new insights into genetic diseases and regenerative medicine in skeletal tissues.

•

Skeletal development is coordinated by master transcription factors.

•

ChIP-seq analyses for the skeletal regulators identified their modes of actions.

•

Analyses of epigenetic landscape features distinct cell types in skeletal tissues.

•

Integrated analyses of the gene regulatory networks link to skeletal regeneration.

An Emerging Regulatory Landscape for Skeletal Development

Skeletal development creates the physical framework that shapes our body and its actions. In the past two decades, genetic studies have provided important insights into the molecular processes at play, including the roles of signaling pathways and transcriptional effectors that coordinate an orderly, progressive emergence and expansion of distinct cartilage and bone cell fates in an invariant temporal and spatial pattern for any given skeletal element within that specific vertebrate species. Genome-scale studies have provided additional layers of understanding, moving from individual genes to the gene regulatory landscape, integrating regulatory information through cis-regulatory modules into cell type-specific gene regulatory programs. This review discusses our current understanding of the transcriptional control of mammalian skeletal development, focusing on recent genome-scale studies.

Transcriptional control of chondrocyte specification and differentiation

A milestone in the evolutionary emergence of vertebrates was the invention of cartilage, a tissue that has key roles in modeling, protecting and complementing the bony skeleton. Cartilage is elaborated and maintained by chondrocytes. These cells derive from multipotent skeletal progenitors and they perform highly specialized functions as they proceed through sequential lineage commitment and differentiation steps. They form cartilage primordia, the primary skeleton of the embryo. They then transform these primordia either into cartilage growth plates, temporary drivers of skeletal elongation and endochondral ossification, or into permanent tissues, namely articular cartilage. Chondrocyte fate decisions and differentiated activities are controlled by numerous extrinsic and intrinsic cues, and they are implemented at the gene expression level by transcription factors. The latter are the focus of this review. Meritorious efforts from many research groups have led over the last two decades to the identification of dozens of key chondrogenic transcription factors. These regulators belong to all types of transcription factor families. Some have master roles at one or several differentiation steps. They include SOX9 and RUNX2/3. Others decisively assist or antagonize the activities of these masters. They include TWIST1, SOX5/6, and MEF2C/D. Many more have tissue-patterning roles and regulate cell survival, proliferation and the pace of cell differentiation. They include, but are not limited to, homeodomain-containing proteins and growth factor signaling mediators. We here review current knowledge of all these factors, one superclass, class, and family at a time. We then compile all knowledge into transcriptional networks. We also identify remaining gaps in knowledge and directions for future research to fill these gaps and thereby provide novel insights into cartilage disease mechanisms and treatment options.