Induced pluripotent stem cells (iPSCs) are potential cell sources for regenerative medi-cine. The iPSCs exhibit a preference for lineage differentiation to the donor cell type indicating the existence of memory of origin. Although the intrinsic effect of the donor cell type on differentiation of iPSCs is well recognized, whether disease-specific factors of donor cells influence the differentiation capacity of iPSC remains unknown. Using viral based reprogramming, we demonstrated the generation of iPSCs from chondrocytes isolated from healthy (AC-iPSCs) and osteoarthritis cartilage (OA-iPSCs). These reprogrammed cells acquired markers of pluripotency and differentiated into uncommitted mesenchymal-like progenitors. Interestingly, AC-iPSCs exhibited enhanced chondrogenic potential as compared OA-iPSCs and showed increased expression of chondrogenic genes. Pan-transcriptome analysis showed that chondrocytes derived from AC-iPSCs were enriched in molecular pathways related to energy metabolism and epigenetic regulation, together with distinct expression signature that distin-guishes them from OA-iPSCs. Our molecular tracing data demonstrated that dysregulation of epigenetic and metabolic factors seen in OA chondrocytes relative to healthy chondrocytes persisted following iPSC reprogramming and differentiation toward mesenchymal progenitors. Our results suggest that the epigenetic and metabolic memory of disease may predispose OA-iPSCs for their reduced chondrogenic differentiation and thus regulation at epigenetic and metabolic level may be an effective strategy for controlling the chondrogenic potential of iPSCs.

Intervertebral disc degeneration is the main contributor to low back pain, now the leading cause of disability worldwide. Gene transfer, either in a therapeutic attempt or in basic research to understand the mechanisms of disc degeneration, is a fascinating and promising tool to manipulate the complex physiology of the disc. Viral vectors based on the adeno-associated virus (AAV) have emerged as powerful transgene delivery vehicles yet a systematic investigation into their respective tropism, transduction efficiency, and relative toxicity have not yet been performed in the disc in vivo. Herein, we used in vivo bioluminescence imaging to systematically compare multiple AAV serotypes, injection volumes, titers, promoters, and luciferase reporters to determine which result in high transduction efficiency of murine nucleus pulposus (NP) cells in vivo. We find that AAV6 using a CAG promoter to drive transgene expression, delivered into the NP of murine caudal discs at a titer of 1011 GC/mL, provides excellent transduction efficiency/kinetics and low toxicity in vivo. We also show, for the first time, that the transduction of NP cells can be significantly boosted in vivo by the use of small cell permeabilization peptides. Finally, to our knowledge, we are the first to demonstrate the use of optical tissue clearing and three-dimensional lightsheet microscopy in the disc, which was used to visualize fine details of tissue and cell architecture in whole intact discs following AAV6 delivery. Taken together, these data will contribute to the success of using AAV-mediated gene delivery for basic and translational studies of the IVD.

SOX4 belongs to the group C of the SOX transcription factor family. It is a critical mediator of tumor necrosis factor alpha (TNF)-induced transformation of fibroblast-like s-ynoviocytes (FLS) in arthritis. In this study we investigated the genome wide association between the DNA binding and transcriptional activities of SOX4 and the NF-kappaB signaling transcription factor RELA/p65 downstream of TNF signaling. We used ChIP-seq assays in mouse FLS to compare the global DNA binding profiles of SOX4 and RELA. RNA-seq of TNF-induced wildtype and SoxC-knockout FLS was used to identify the SOX4-dependent and independent aspects of the TNF-regulated transcriptome. We found that SOX4 and RELA physically interact with each other on the chromatin. Interestingly, ChIP-seq assays revealed that 70.4% of SOX4 peak summits were within 50bp of the RELA peak summits suggesting that both proteins bind in close-proximity on regulatory sequences, enabling them to co-operatively regulate gene expression. By integrating the ChIP-seq results with RNA-seq from SoxC-knockout FLS we identified a set of TNF-responsive genes that are targets of the RELA-SOX4 transcriptional complex. These TNF-responsive and RELA-SOX4-depenedent genes included inflammation mediators, histone remodeling enzymes and components of the AP-1 signaling pathway. We also identified an autoregulatory mode of SoxC gene expression that involves a TNF-mediated switch from RELA binding to SOX4 binding in the 3’ UTR of Sox4 and Sox11 genes. In conclusion, our results show that SOX4 and RELA together orchestrate a multimodal regulation of gene expression downstream of TNF signaling. Their interdependent activities play a pivotal role in the transformation of FLS in arthritis and in the inflammatory pathology of diverse tissues where RELA and SOX4 are co-expressed.

Intervertebral disc degeneration (IDD) is a public health dilemma as it is associated with low back and neck pain, a frequent reason for patients to visit the physician. During IDD, nucleus pulposus (NP), the central compartment of intervertebral disc (IVD) undergo degeneration. Stem cells have been adopted as a promising biological source to regenerate the IVD and restore its function. Here, we describe a simple, two-step differentiation strategy using a cocktail of four factors (LDN, AGN, FGF, and CHIR) for efficient derivation of notochordal cells from human embryonic stem cells (hESCs). We employed a CRISPR/Cas9 based genome-editing approach to knock-in the mCherry reporter vector upstream of the 3′ untranslated region of the Noto gene in H9-hESCs and monitored notochordal cell differentiation. Our data show that treatment of H9-hESCs with the above-mentioned four factors for 6 days successfully resulted in notochordal cells. These cells were characterized by morphology, immunostaining, and gene and protein expression analyses for established notochordal cell markers including FoxA2, SHH, and Brachyury. Additionally, pan-genomic high-throughput single cell RNA-sequencing revealed an efficient and robust notochordal differentiation. We further identified a key regulatory network consisting of eight candidate genes encoding transcription factors including PAX6, GDF3, FOXD3, TDGF1, and SOX5, which are considered as potential drivers of notochordal differentiation. This is the first single cell transcriptomic analysis of notochordal cells derived from hESCs. The ability to efficiently obtain notochordal cells from pluripotent stem cells provides an additional tool to develop new cell-based therapies for the treatment of IDD.

Craniofacial bone loss is a complex clinical problem with limited regenerative solutions. Currently, BMP2 is used as a bone-regenerative therapy in adults, but in pediatric cases of bone loss, it is not FDA-approved due to concerns of life-threatening inflammation and cancer. Development of a bone-regenerative therapy for children will transform our ability to reduce the morbidity associated with current autologous bone grafting techniques. We discovered that JAGGED1 (JAG1) induces cranial neural crest (CNC) cell osteoblast commitment during craniofacial intramembranous ossification, suggesting that exogenous JAG1 delivery is a potential craniofacial bone-regenerative approach. In this study, we found that JAG1 delivery using synthetic hydrogels containing O9–1 cells, a CNC cell line, into critical-sized calvarial defects in C57BL/6 mice provided robust bone-regeneration. Since JAG1 signals through canonical (Hes1/Hey1) and non-canonical (JAK2) NOTCH pathways in CNC cells, we used RNAseq to analyze transcriptional pathways activated in CNC cells treated with JAG1 ± DAPT, a NOTCH-canonical pathway inhibitor. JAG1 upregulated expression of multiple NOTCH canonical pathway genes (Hes1), which were downregulated in the presence of DAPT. JAG1 also induced bone chemokines (Cxcl1), regulators of cytoskeletal organization and cell migration (Rhou), signaling targets (STAT5), promoters of early osteoblast cell proliferation (Prl2c2, Smurf1 and Esrra), and, inhibitors of osteoclasts (Id1). In the presence of DAPT, expression levels of Hes1 and Cxcl1 were decreased, whereas, Prl2c2, Smurf1, Esrra, Rhou and Id1 remain elevated, suggesting that JAG1 induces osteoblast proliferation through these non-canonical genes. Pathway analysis of JAG1 + DAPT-treated CNC cells revealed significant upregulation of multiple non-canonical pathways, including the cell cycle, tubulin pathway, regulators of Runx2 initiation and phosphorylation of STAT5 pathway. In total, our data show that JAG1 upregulates multiple pathways involved in osteogenesis, independent of the NOTCH canonical pathway. Moreover, our findings suggest that JAG1 delivery using a synthetic hydrogel, is a bone-regenerative approach with powerful translational potential.

During craniofacial development, cranial neural crest (CNC) cells migrate into the developing face and form bone through intramembranous ossification. Loss of JAGGED1 (JAG1) signaling in the CNC cells is associated with maxillary hypoplasia or maxillary bone deficiency (MBD) in mice and recapitulates the MBD seen in humans with Alagille syndrome. JAGGED1, a membrane-bound NOTCH ligand, is required for normal craniofacial development, and Jagged1 mutations in humans are known to cause Alagille Syndrome, which is associated with cardiac, biliary, and bone phenotypes and these children experience increased bony fractures. Previously, we demonstrated deficient maxillary osteogenesis in Wnt1-cre;Jagged1f/f (Jag1CKO) mice by conditional deletion of Jagged1 in maxillary CNC cells. In this study, we investigated the JAG1 signaling pathways in a CNC cell line. Treatment with JAG1 induced osteoblast differentiation and maturation markers, Runx2 and Ocn, respectively, Alkaline Phosphatase (ALP) production, as well as classic NOTCH1 targets, Hes1 and Hey1. While JAG1-induced Hes1 and Hey1 expression levels were predictably decreased after DAPT (NOTCH inhibitor) treatment, JAG1-induced Runx2 and Ocn levels were surprisingly constant in the presence of DAPT, indicating that JAG1 effects in the CNC cells are independent of the canonical NOTCH pathway. JAG1 treatment of CNC cells increased Janus Kinase 2 (JAK2) phosphorylation, which was refractory to DAPT treatment, highlighting the importance of the non-canonical NOTCH pathway during CNC cells osteoblast commitment. Pharmacologic inhibition of JAK2 phosphorylation, with and without DAPT treatment, upon JAG1 induction reduced ALP production and, Runx2 and Ocn gene expression. Collectively, these data suggest that JAK2 is an essential component downstream of a non-canonical JAG1-NOTCH1 pathway through which JAG1 stimulates expression of osteoblast-specific gene targets in CNC cells that contribute to osteoblast differentiation and bone mineralization.

Craniosynostosis (CS), the premature fusion of one or more cranial sutures, is a relatively common congenital anomaly, occurring in 3–5 per 10,000 live births. Nonsyndromic CS (NCS) accounts for up to 80% of all CS cases, yet the genetic factors contributing to the disorder remain largely unknown. The RUNX2 gene, encoding a transcription factor critical for bone and skull development, is a well known CS candidate gene, as copy number variations of this gene locus have been found in patients with syndromic craniosynostosis. In the present study, we aimed to characterize RUNX2 to better understand its role in the genetic etiology and in the molecular mechanisms underlying midline suture ossification in NCS. We report four nonsynonymous variants, one intronic variant and one 18 bp in-frame deletion in RUNX2 not found in our study control population. Significant difference in allele frequency (AF) for the deletion variant RUNX2 p.Ala84-Ala89del (ClinVar 257,095; dbSNP rs11498192) was observed in our sagittal NCS cohort when compared to the general population (P = 1.28 × 10−6), suggesting a possible role in the etiology of NCS. Dual-luciferase assays showed that three of four tested RUNX2 variants conferred a gain-of-function effect on RUNX2, further suggesting their putative pathogenicity in the tested NCS cases. Downregulation of RUNX2 expression was observed in prematurely ossified midline sutures. Metopic sites showed significant downregulation of promoter 1-specific isoforms compared to sagittal sites. Suture-derived mesenchymal stromal cells showed an increased expression of RUNX2 over matched unfused suture derived cells. This demonstrates that RUNX2, and particularly the distal promoter 1-isoform group, are overexpressed in the osteogenic precursors within the pathological suture sites.

Fibroblast-like synoviocytes (FLS) play a critical role in the pathogenesis of rheumatoid arthritis (RA). Chronic inflammation induces transcriptomic and epigenetic modifications that imparts a persistent catabolic phenotype to the FLS, despite their dissociation from the inflammatory environment. We analyzed high throughput gene expression and chromatin accessibility data from human and mouse FLS from our and other studies available on public repositories, with the goal of identifying the persistently reprogrammed signaling pathways driven by chronic inflammation. We found that the gene expression changes induced by short-term tumor necrosis factor-alpha (TNF) treatment were largely sustained in the FLS exposed to chronic inflammation. These changes that included both activation and repression of gene expression, were accompanied by the remodeling of chromatin accessibility. The sustained activated genes (SAGs) included established pro-inflammatory signaling components known to act at multiple levels of NF-kappaB, STAT and AP-1 signaling cascades. Interestingly, the sustained repressed genes (SRGs) included critical mediators and targets of the BMP signaling pathway. We thus identified sustained repression of BMP signaling as a unique constituent of the long-term inflammatory memory induced by chronic inflammation. We postulate that simultaneous targeting of these activated and repressed signaling pathways may be necessary to combat RA persistence.

Intervertebral disc (IVD) disease (IDD) is a complex, multifactorial disease. While various aspects of IDD progression have been reported, the underlying molecular pathways and transcriptional networks that govern the maintenance of healthy nucleus pulposus (NP) and annulus fibrosus (AF) have not been fully elucidated. We defined the transcriptome map of healthy human IVD by performing single-cell RNA-sequencing (scRNA-seq) in primary AF and NP cells isolated from non-degenerated lumbar disc. Our systematic and comprehensive analyses revealed distinct genetic architecture of human NP and AF compartments and identified 2,196 differentially expressed genes. Gene enrichment analysis showed that SFRP1, BIRC5, CYTL1, ESM1 and CCNB2 genes were highly expressed in the AF cells; whereas, COL2A1, DSC3, COL9A3, COL11A1, and ANGPTL7 were mostly expressed in the NP cells. Further, functional annotation clustering analysis revealed the enrichment of receptor signaling pathways genes in AF cells, while NP cells showed high expression of genes related to the protein synthesis machinery. Subsequent interaction network analysis revealed a structured network of extracellular matrix genes in NP compartments. Our regulatory network analysis identified FOXM1 and KDM4E as signature transcription factor of AF and NP respectively, which might be involved in the regulation of core genes of AF and NP transcriptome.

The bone healing environment in the posterolateral spine following arthrodesis surgery is one of the most challenging in all of orthopedics and our understanding of the molecular signaling pathways mediating osteogenesis during spinal fusion is limited. In this study, the spatial and temporal expression pattern of Wnt signaling factors and inhibitors during spinal fusion was assessed for the first time. Bilateral posterolateral spine arthrodesis with autologous iliac crest bone graft was performed on 21 New Zealand White rabbits. At 1-, 2-, 3-, 4-, and 6-weeks, the expression of sclerostin and a variety of canonical and noncanonical Wnts signaling factors was measured by qRT-PCR from tissue separately collected from the transverse processes, the Outer and Inner Zones of the fusion mass, and the adjancent paraspinal muscle. Immunohistochemistry for sclerostin protein was also performed. Sclerostin and many Wnt factors, especially Wnt3a and Wnt5a, were found to have distinct spatial and temporal expression patterns. For example, harvesting ICBG caused a significant increase in sclerostin expression. Furthermore, the paraspinal muscle immediately adjacent to the transplanted ICBG also had significant increases in sclerostin expression at 3 weeks, suggesting new potential mechanisms for pseudarthroses following spinal arthrodesis. The presented work is the first description of the spatial and temporal expression of sclerostin and Wnt signaling factors in the developing spine fusion, filling an important knowledge gap in the basic biology of spinal fusion and potentially aiding in the development of novel biologics to increase spinal fusion rates.