RNA-Seq or whole transcriptome shotgun sequencing has been adopted in the last years as a reference technique to determine the presence and the quantity of different species of RNA in determined biological samples, thanks to it allows the identification every single RNA species transcribed from a reference genome. Meta-profiling takes advantage of the public availability of an increasing set of RNA-Seq data produced by different laboratories to summarize the expression levels of the different RNA species of many samples according to their biological context, giving the opportunity to perform comparisons on the gene expression profiles of different tissues by integrating data derived from a high number of studies. By using Genevestigator™; a platform which integrates RNA-Seq data into meta-profiles, we have performed a comparison between the gene expression profiles of bone, cartilage, muscle tendon and skin by means of interrogating its database with different gene sets and families with relevance to the function of the tissues of the musculoskeletal system. The collagen gene family and genes coding for proteoglycans, matrix metalloproteinases and tissue inhibitors of metalloproteinases, mechanotransduction-related proteins and signalling pathways involved in tissue development and differentiation have been analysed. Hierarchical clustering for every gene set was performed for the understanding the differences and similarities between the different tissues included in the analyses. The results of this study will help to improve our understanding of the musculoskeletal system, and will help to identify new biomarkers and signalling pathways of specific relevance for the bone, cartilage, muscle and tendon.

Heterotopic ossification (HO) is defined as aberrant bone formation in extraskeletal locations. In this process, local stromal cells of mesenchymal origin abnormally differentiate, resulting in pathologic cartilage and bone matrix deposition. However, the specific cell type and mechanisms beyond this process are not well understood, in part due to the heterogeneity of progenitor cells involved. Here, a combination of single cell RNA sequencing (scRNA-Seq) and lineage tracing, defined the extent to which synovial / tendon sheath progenitor cells contribute to HO. For this purpose, a Tppp3 (tubulin polymerization-promoting protein family member 3) inducible reporter model was used, in combination with either Scx (Scleraxis) or Pdgfra (Platelet derived growth factor receptor alpha) reporter animals. Both arthroplasty-induced and tendon injury-mouse experimental HO models were utilized. ScRNA-Seq of tendon-induced traumatic HO suggested that Tppp3 is a progenitor cell marker for either osteochondral or tendon or cells. After HO induction, Tppp3 reporter+ cell population expanded in number and contributed to cartilage and bone formation in tendon and joint-associated HO. Using double reporter animals, we found that both Pdgfra+Tppp3+ and Pdgfra+Tppp3- progenitor cells produced HO-associated cartilage. Finally, the examination of human samples showed a significant population of TPPP3+ cells overlapping with osteogenic markers in areas of HO. Overall, these results provide novel observations that peritenon and synovial progenitor cells undergo abnormal osteochondral differentiation and contribute to heterotopic bone formation after trauma

Stem cell therapy is an effective means to address the repair of large segmental bone defects. However, the intense inflammatory response triggered by the implants severely impairs stem cell differentiation and tissue regeneration. High-dose transforming growth factor β1 (TGF-β1), the most locally expressed cytokine in implants, inhibits osteogenic differentiation of bone marrow mesenchymal stem cells (BMSCs) and promotes tissue fibrosis, severely compromising the efficacy of stem cell therapy. Small molecule inhibitors of TGF-β1 can be used to ameliorate the osteogenic disorders caused by high concentrations of TGF-β1, but systemic inhibition of TGF-β1 function will cause strong adverse effects. How to find safe and reliable molecular targets to antagonize TGF-β1 remains to be elucidated. Orphan nuclear receptor Nr4a1, an endogenous inhibitory molecule of TGF-β1, suppresses tissue fibrosis, but its role in BMSC osteogenesis is unclear. We found that TGF-β1 inhibited Nr4a1 expression through HDAC4. Overexpression of Nr4a1 in BMSCs reversed osteogenic differentiation inhibited by high levels of TGF- β1. Mechanistically, RNA sequencing showed that Nr4a1 activated the ECM-receptor interaction and Hippo signaling pathway, which in turn promoted BMSC osteogenesis. In bone defect repair and fracture healing models, transplantation of Nr4a1-overexpressing BMSCs into C57BL/6J mice or treatment with the Nr4a1 agonist Csn-B significantly ameliorated inflammation-induced bone regeneration disorders. In summary, our findings confirm the endogenous inhibitory effect of Nr4a1 on TGF- β1 and uncover the effectiveness of Nr4a1 agonists as a therapeutic tool to improve bone regeneration, which provides a new solution strategy for the treatment of clinical bone defects and inflammatory skeletal diseases

To address the current challenge of anterior cruciate ligament (ACL) reconstruction, this study is the first to fabricate a braided collagen rope (BCR) which mimics native hamstring for ACL reconstruction. The study aims to evaluate the biological and biomechanical properties of BCR both in vivo and vitro. Rabbit ACL reconstruction model using collagen rope and autograft (hamstring tendon) was conducted. The histological and biomechanical evaluations were conducted at 6-, 12-, 18, 26-week post-operation. In vitro study included cell morphology analysis, cell function evaluation and RNA sequencing of the tenocytes cultured on BCR. A cadaver study was also conducted to verify the feasibility of BCR for ACL reconstruction. BCR displays satisfactory mechanical strength similar to hamstring graft for ACL reconstruction in rabbit. Histological assessment showed BCR restore ACL morphology at 26 weeks similar to native ACL. The superior dynamic ligamentization in BCR over autograft group was evidenced by assessment of cell and collagen morphology and orientation. The in vitro study showed that the natural collagen fibres within BCR enables to signal the morphology adaptation and orientation of human tenocytes in bioreactor. BCR enables to enhance cell proliferation and tenogenic expression of tenocytes as compared to hydrolysed collagen. We performed an RNA-Sequencing (RNA-seq) experiment where RNA was extracted from tenocyte seeded with BCR. Analysis of enriched pathways of the up-regulated genes revealed that the most enriched pathways were the Hypoxia-inducible factor 1-alpha (HIF1A) regulated networks, implicating the possible mechanism BCR induced ACL regeneration. The subsequent cadaver study was conducted to proof the feasibility of BCR for ACL reconstruction. This study demonstrated the proof-of-concept of bio-textile braided collagen rope for ACL reconstruction, and the mechanism by which BCR induces natural collagen fibres that positively regulate morphology and function of tenocytes

Despite osteoarthritis (OA) representing a large burden for healthcare systems, there remains no effective intervention capable of regenerating the damaged cartilage in OA. Mesenchymal stromal cells (MSCs) are adult-derived, multipotent cells which are a candidate for musculoskeletal cell therapy. However, their precise mechanism of action remains poorly understood. The effects of an intra-articular injection of human bone-marrow derived MSCs into a knee osteochondral injury model were investigated in C57Bl/6 mice. The cell therapy was retrieved at different time points and single cell RNA sequencing was performed to elucidate the transcriptomic changes relevant to driving tissue repair. Mass cytometry was also used to study changes in the mouse immune cell populations during repair. Histological assessment reveals that MSC treatment is associated with improved tissue repair in C57Bl/6 mice. Single cell analysis of retrieved human MSCs showed spatial and temporal transcriptional heterogeneity between the repair tissue (in the epiphysis) and synovial tissue. A transcriptomic map has emerged of some of the distinct genes and pathways enriched in human MSCs isolated from different tissues following osteochondral injury. Several MSC subpopulations have been identified, including proliferative and reparative subpopulations at both 7 days and 28 days after injury. Supported by the mass cytometry results, the immunomodulatory role of MSCs was further emphasised, as MSC therapy was associated with the induction of increased numbers of regulatory T cells correlating with enhanced repair in the mouse knee. The transcriptomes of a retrieved MSC therapy were studied for the first time. An important barrier to the translation of MSC therapies is a lack of understanding of their heterogeneity, and the consequent lack of precision in its use. MSC subpopulations with different functional roles may be implicated in the different phases of tissue repair and this work offers further insights into repair process

During aging, tendons demonstrate substantial disruptions in homeostasis, leading to impairments in structure-function. Impaired tendon function contributes to substantial declines quality of life during aging. Aged tendons are more likely to undergo spontaneous rupture, and the healing response following injury is impaired in aged tendons. Thus, there is a need to develop strategies to maintain tendon homeostasis and healing capacity through the lifespan. Tendon cell density sharply declines by ∼12 months of age in mice, and this low cell density is retained in geriatric tendons. Our data suggests that this decline in cellularity initiates a degenerative cascade due to insufficient production of the extracellular matrix (ECM) components needed to maintain tendon homeostasis. Thus, preventing this decline in tendon cellularity has great potential for maintaining tendon health. Single cell RNA sequencing analysis identifies two changes in the aged tendon cell environment. First, aged tendons primarily lose tenocytes that are associated with ECM biosynthesis functions. Second, the tenocytes that remain in aged tendons have disruptions in proteostasis and an increased pro-inflammatory phenotype, with these changes collectively termed ‘programmatic skewing'. To determine which of these changes drives homeostatic disruption, we developed a model of tenocyte depletion in young animals. This model decreases tendon cellularity to that of an aged tendon, including decreased biosynthetic tenocyte function, while age-related programmatic skewing is absent. Loss of biosynthetic tenocyte function in young tendons was sufficient to induce homeostatic disruption comparable to natural aging, including deficits in ECM organization, composition, and material quality, suggesting loss biosynthetic tenocytes as an initiator of tendon degeneration. In contrast, our data suggest that programmatic skewing underpins impaired healing in aged tendons. Indeed, despite similar declines in the tenocyte environment, middle-aged and young-depleted tendons mount a physiological healing response characterized by robust ECM synthesis and remodeling, while aged tendons heal with insufficient ECM

Joint tissues release extracellular vesicles (EVs) that potentially sustain joint homeostasis and contribute to osteoarthritis (OA) pathogenesis. EVs are putative novel therapeutics for OA, and transport biologically active molecules (including small non-coding RNAs (SNCRNAs)) between cells. This study identified altering SNCRNA cargo in EVs in OA which may act as early diagnostic markers and treatment targets. OA was surgically induced in four skeletally mature Standardbred horses using an osteochondral fragment model in the left middle carpal joint. The right joint underwent sham surgery. Synovial fluid (SF) and plasma were obtained weekly throughout the 70-day study. EVs were isolated using size exclusion chromatography and characterised using nanoparticle tracking (Nanosight), and exosome fluorescence detection and tetraspanin phenotyping (Exoview). RNA was extracted from EVs derived from SF (sham and OA joints) and plasma collected at days 10, 35, 42, 49, 56, 63, and subjected to small RNA sequencing on a NovaSeq SP100 flow cell (Illumina). Nanosight-derived EV characteristics of size and concentration were not significantly different following disease induction. The diameter of the temporal population of plasma and SF-derived exosomes changed significantly for CD9 and CD81 following OA induction with significant temporal, and disease-related changes in CD63 and CD81 protein expressin in plasma and SF. In SF and plasma-derived EVs snoRNAs, snRNAs, tRNAs, lncRNA, y-RNA, piRNAs and scRNA were found. Following pairwise analysis of all-time points we identified 27 miRs DE in plasma and 45 DE miRs in SF. Seven were DE in plasma and SF; miR-451, miR-25, miR-215, miR-92a, miR-let-7c, miR-486-5p, miR-23a. In plasma and SF 35 and 21 snoRNAs were DE with four DE in plasma and SF; U3, snord15, snord46, snord58. This work has identified alterations to OA EV sncRNAs in plasma and SF providing a greater understanding of the role of EVs in early OA

We found that adipose stem cells are poorly differentiated into bone and that their ability to differentiate into bone varies from cell line to cell line. The osteogenic differentiation ability of the adipose stem cell lines was distinguished through Alzarin Red Staining, and the cell lines that performed well and those that did not were subjected to RNA-seq analysis. The selected gene GSTT1 (glutathione S-transferase theta-1) gene is a member of a protein superfamily that catalyzes the conjugation of reduced glutathione to a variety of hydrophilic and hydrophobic compounds. The purpose of this study is to treat avascular necrosis and bone defect by improving bone regeneration with adipose stem cells introduced with a new GSTT1 gene related to osteogenic differentiation of adipose stem cells. In addition, the GSTT1 gene has the potential as a genetic marker that can select a specific cell line in the development of an adipose stem cell bone regeneration drug. Total RNA was extracted from each sample using the TRIzol reagent. Its concentration and purity were determined based on A260 and A260/A280, respectively, using a spectrophotometer. RNA sequencing library of each sample was prepared using a TruSeq RNA Library Prep Kit. RNA-seq experiments were performed for hADSCs. Cells were transfected with either GSTT1 at 100 nM or siControl (scramble control) by electroporation using a 1050 pulse voltage for 30 ms with 2 pulses using a 10 μl pipette tip. The purpose of this study is to discover genetic markers that can promote osteogenic differentiation of adipose stem cells (hADSCs) through mRNA-seq gene analysis. The selected GSTT1 gene was found to be associated with the enhancement of osteogenic differentiation of adipose stem cells. siRNA against GSTT1 reduced osteogenic differentiation of hADSCs, whereas GSTT1 overexpression enhanced osteogenic differentiation of hADSCs under osteogenic conditions. In this study, GSTT1 transgenic adipose stem cells could be used in regenerative medicine to improve bone differentiation. In addition, the GSTT1 gene has important significance as a marker for selecting adipose stem cells with potential for bone differentiation in the development of a therapeutic agent for bone regeneration cells

Though dentin matrix protein 1 (Dmp1) is known to play critical role in mediating bone mineralization, it has also been validated to be expressed in brain and helps maintain blood brain barrier (BBB). Our study aims to clarify the expression pattern of Dmp1 in mouse brain and explore whether intercellular mitochondrial transfer occurs between Dmp1 positive astrocytes (DPAs) and endothelial cells, and thus acting as a mechanism in maintaining BBB during aging. Single cell RNA sequencing (scRNAseq) of 1 month, 6 month, and 20 month old mice brain (n=1, respectively) was employed to identify Dmp1 positive cell types. Dmp1. Cre. -mGmT and Dmp1. Cre. -COX8a fluorescent mice were generated to visualize DPAs and investigate their mitochondrial activities. A 3D noncontact coculture system and mitochondrial transplantation were applied to study the role of mitochondrial transfer between astrocytes and bEnd.3 endothelial cells. Dmp1. Cre. -Mfn2. f/f. mice were generated by depleting the ER-mitochondria tethering protein Mfn2 in DPAs. Dmp1 was mainly expressed in astrocytes at different ages. GO analysis revealed that cell projection and adhesion of DPAs were upregulated. Confocal imaging on Dmp1. Cre. -mGmT mice indicated that DPAs are a cluster of astrocytes that closely adhere to blood vessels (n=3). Bioinformatics analysis revealed that mitochondrial activity of DPAs were compromised during aging. Enriched scRNAseq of fluorescent cells from Dmp1. Cre. -COX8a mice (n=2) and immunofluorescent imaging (n=3) validated the acquisition of extrinsic mitochondria in endothelial cells. 3D coculture of astrocytes and bEnd.3 and direct mitochondrial transplantation revealed the rescue effect of mitochondrial transfer on damaged bEnd.3. BBB was impaired after depleting Mfn2 in DPAs, expressing a similar phenotype with aging brain. Astrocytes that express Dmp1 play a significant role in maintaining BBB via transferring mitochondria to vascular endothelial cells. Compromised mitochondrial transfer between DPAs and endothelial cells might be the potential mechanism of impaired BBB during aging

Abstract. Objectives. Mesenchymal stromal/stem cells (MSCs) are increasingly recognized as regulators of immune cells during disease or tissue repair. During these situations, the extracellular matrix (ECM) is very dynamic and therefore, our studies aim to understand how ECM influences the activity of MSCs. Methods. Human MSCs cultured on tissue culture plastic (TCP) and encapsulated within collagen type I, fibrin, or mixed Collagen-Fibrin were exposed to low dose TNFα and IFNɣ. Transcription profiles were examined using bulk RNA sequencing (RNAseq) after 24h of treatment. ELISA, Western blot, qPCR and immunofluorescence were employed to validate RNAseq results and to investigate the significance of transcriptional changes. Flow cytometry evaluated monocyte/macrophage phenotype. Results. Previously, we showed that human MSC expression of TNFAIP6 and CXCL10 in 3D environments is significantly upregulated in response to pro-inflammatory stimuli. Here, RNAseq revealed that there were 2,085 highly significant upregulated genes in 3D matrices compared to TCP. Notably, >90% of highly expressed genes (including FOSB, FOS and TNFAIP6) were shared in all hydrogels. Gene ontology confirmed the TNF signalling pathway among the most significantly represented. Protein-protein interaction predictions identified TNF-alpha/NF-kappa B and AP1 pathways as differentially influenced by the hydrogel environment. Using inhibitors to these pathways, NFkB, but not AP1, impacted on the upregulation of TNFAIP6 and CXCL10 in 3D culture. Conditioned media from these studies was added to cultures of human monocytes with distinct changes in the resulting macrophage phenotype. MSCs in a 3D environment promoted a greater acquisition of the M2 repair macrophage phenotype and impacted on the numbers of pro-inflammatory M1 macrophages. Conclusion. These data provide further evidence that the immunomodulatory action of human MSCs can be influenced by the surrounding structural environment. These observations have significance for understanding the events that following skeletal injury and the potential to be exploited in preconditioning MSCs for cell therapy