Skip to main content

Identifying the key genes regulating mesenchymal stem cells chondrogenic differentiation: an in vitro study



Mesenchymal stem cells (MSCs) possess the potential to differentiate into chondrocytes, which makes them an ideal source for healing cartilage defects. Here, we seek to identify the essential genes participating in MSCs chondrogenesis.


Human MSCs were induced for chondrogenesis for 7, 14, and 21 days using a high-density micromass culture system, and RNA was extracted for RNA-seq.


A total of 6247 differentially expressed genes (DEGs) were identified on day 7, and 85 DEGs were identified on day 14. However, no significant DEGs was identified on day 21. The top 30 DEGs at day 7, including COL9A3, COL10A1, and CILP2, are closely related to extracellular matrix organization. While the top 30 DEGs at day 14 revealed that inflammation-related genes were enriched, including CXCL8, TLR2, and CCL20. We also conducted protein–protein interaction (PPI) networks analysis using the search tool for the retrieval of interacting genes (STRING) database and identified key hub genes, including CXCL8, TLR2, CCL20, and MMP3. The transcriptional factors were also analyzed, identifying the top 5 TFs: LEF1, FOXO1, RORA, BHLHE41, and SOX5. We demonstrated one particular TF, RORA, in promoting early MSCs chondrogenesis.


Taken together, our results suggested that these DEGs may have a complex effect on MSCs chondrogenesis both synergistically and solitarily.

Peer Review reports


Cartilage tissues are responsible for attenuating force and providing lubrication for joint movement. Due to aging, sports injuries, or repeated use, cartilage defect is a major challenge in orthopedic surgery. Articular cartilage has limited regenerative ability, and cartilage injury can progress to osteoarthritis if poorly treated [1]. Stem cells combined with scaffolds have been used clinically as cartilage defect therapy [2]. The crucial process for these cell-based scaffolds to heal the cartilage defect requires cartilage differentiation (or chondrogenesis), but the regulatory factors involved are still not fully understood.

During cartilage development, the mesenchymal stem cells are recruited and undergo mesenchymal condensation and chondrogenic differentiation [3, 4]. Chondrogenesis is a tightly regulated process in which mesenchymal progenitors differentiate into mature chondrocytes and secrete extracellular matrix (ECM) components. The chondrocytes in the articular cartilage maintain the resting state and constitute hyaline cartilage. In pathologic conditions such as osteoarthritis's early stage, the chondrocytes undergo hypertrophic differentiation [5]. During chondrocyte hypertrophy, the ECM secreted by chondrocytes was changed from type II to type X collagen, leading to osteoarthritis progression [6]. However, the gene expression profile during the transition from resting to hypertrophic chondrocyte remained unclear.

Adult mesenchymal stem cells (MSCs) are pluripotent stem cells and can differentiate into multiple cell types, including osteocytes, chondrocytes, tenocytes, and adipocytes [6, 7]. Besides, MSCs exhibit other exciting properties, such as anti-inflammation, anti-senescence, and anti-allergic [8, 9]. Because of these properties and their simplicity of obtaining, MSCs injection has been used as a treatment for multiple clinical conditions [10]. Previous studies have shown that multiple genes and pathways are crucial to the chondrogenesis of MSCs, such as SOX9 and BMP signaling pathways [11]. However, despite these well-established genes and pathways that regulate chondrogenesis, emerging evidence has shown that chondrogenesis is regulated by a complex network of other participants, such as non-coding RNA and growth factors [12]. Therefore, further studies are still needed to analyze the dynamic change during MSC chondrogenesis.

To investigate the vital genes in regulating MSCs chondrogenesis, it is important to characterize the gene expression pattern during chondrogenesis. In this study, MSCs chondrogenesis was induced by a high-density micromass culture system. We then performed RNA-seq and subsequent analysis at different time points during MSCs chondrogenic differentiation. We hope this study may enhance our understanding of the gene profile during chondrogenesis. More importantly, identifying the key transcriptional factors that initiate and maintain MSCs chondrogenesis in vitro may enhance our comprehension on the mechanism of induced chondrogenesis in vivo.


MSC isolation and culture

The Ethical Committee of Sun Yat-sen University approved this study. Every participant enrolled in the study signed a formal consent form for collecting MSCs. The MSCs from three male participants with an average age of 37.6 ± 4.2 years were used for chondrogenic differentiation. Briefly, the bone marrow aspirate was collected from the posterior iliac crest. Bone marrow aspirate was collected using a 14-gauge needle connected to a 20-mL syringe. Approximately 2 mL of bone marrow was collected from each puncture in the posterior iliac crest. The bone marrow aspirate was immediately transferred to a vacutainer tube containing the anticoagulant. MSCs were separated from bone marrow aspirates as described previously [13, 14]. Briefly, the bone marrow aspirates were diluted with sterile phosphate buffer saline (PBS) and then fractionated with density gradient by 500 g for 30 min. The mononuclear cells layer was carefully collected and resuspended in low-glucose Dulbecco’s modified Eagle medium (DMEM; Gibco, Waltham, MA, USA). The MSCs were incubated at 37 ℃ incubators at 5% CO2. The medium was replaced every 3 days and was digested with 0.25% trypsin (Gibco, Waltham, MA, USA) for passage.

Chondrogenic differentiation

As previously described, the chondrogenic differentiation of MSCs was achieved using a high-density micromass culture system [14, 15]. Only MSCs from passages 3–5 were subjected to chondrogenic differentiation. The procedure of chondrogenic differentiation was described as follows: MSCs were separated with 0.25% trypsin and washed with PBS. MSCs were suspended at a cell number of 2×107 Cells/ml in Mesenchymal Stem Cell Chondrogenic Differentiation Medium (MCDM, ScienCell) supplemented with 10ug/ml TGF-β3 (PEPROTECH, USA, Cat No: 100-36E). Each MSCs suspension droplet (12.5 μl) was carefully placed on a 24-well plate. The droplets were then placed in the 37 ℃ incubators for 2 h to evaporate liquid and form aggregate. Then 500 μl chondrogenic medium was added to each well. The medium was replaced every 3 days, and the pellets were harvested after 7, 14, and 21 days of induction. The CHO group was chondrocytes isolated from an age-paired patient who underwent total knee arthroplasty (TKA) surgery. The relatively intact part of the knee cartilage was dissected, the chondrocytes were digested, and total RNA was extracted.

RNA extraction and RNA-seq

Total RNA was extracted from the cartilage pellets with TRIzol reagent (Invitrogen). A total amount of 1 µg RNA per sample was used as input material for the RNA sample preparations. Cartilage pellets were sent to Novogene (Beijing, China) for RNA-seq analysis. Briefly, RNA integrity number (RIN) was evaluated by Agilent 2100 Bioanalyzer. The sequencing libraries were generated and sequenced by an Illumina HiSeq X Ten sequencer (Illumina).

Data processing and the differentially expressed genes (DEGs) analysis

We calculated the RPKM (reads per thousand base exons per million mapped reads) value to evaluate gene expression. The R software package FlashClust was used to cluster the genes expressed in different samples. The DESeq package was used for differential gene expression analysis at different time-point. The false discovery rate (FDR) was controlled by adjusting the p-value using the Benjamin-Hochberg algorithm. The significant DEGs were identified as log2FoldChange > 1 or < -1, P < 0.05. The R package ggPlot2 and pheatmap were used to visualize the DEGs.

Protein–protein interaction (PPI) network analysis

The DEGs in the form of gene entrez was uploaded to the STRING database for PPI network analysis ( Briefly, the selection criteria were as follows: organism: Homo sapiens; network: full network; interaction score: high confidence. Each hub gene is clustered by kmeans clustering. The networks were further analyzed by Cytoscape version 3.8.2 ( to identify hub genes and their relationships. The inserted app cytoHubba and MCODE was used to assist the analysis.

Statistical analysis

The numeric data were presented as the mean ± SD of three replicates. SPSS 23.0 was used to perform the statistical analysis. GraphPad Prism 8.0 software was used to visualize the results. Differences between the two grouped data were analyzed with the student's t-test.


Identification of DEGs during MSCs chondrogenesis

The schematic illustration of this experiment was shown in Fig. 1A. We performed immunohistochemistry(IHC) staining to validate the expression of the chondrogenic marker COL2A1. The successful induction of MSCs chondrogenesis was confirmed (Fig. 1B). The samples were then collected and subjected to RNA sequencing, described in the method section. Pearson’s correlation between different groups was presented as R2 value. The in-group difference was more negligible than the inter-group (Fig. 1C). As shown in Fig. 1D, there was a significant difference in the gene expression profile between different time points.

Fig. 1
figure 1

Study design and heatmap showing all DEGs. A The expreimental procedure of this study. B The expression of chondrogenic marker COL2A1 in chondrocyte pellet after induced for 7, 14, and 21 days were assessed by IHC. C The correlation coefficient between different biological replicates and indicated time points. MSCs = undiferentiated MSCs, C7 = induced chondrogenesis for 7 days, C14 = induced chondrogenesis for 14 days, C21 = induced chondrogenesis for 21 days, CHO = mature chondrocytes. D Heatmap showing all DEGs. The expression values are |log2 fold changes|≥ 1 and P ≤ 0.05

The pellet harvested on day 7 was identified 6247 DEGs compared to undifferentiated MSCs. Among the 6247 DEGs, 3333 genes were up-regulated, and 2914 genes were down-regulated (Fig. 2A). Compared to the pellet harvested on day 7, a total of 85 DEGs was identified in the pellet harvested on day 14. Among the 85 DEGs, 24 were up-regulated, and 61 were down-regulated (Fig. 2B). The top 30 DEGs at 7 days after induced chondrogenesis was shown in Table 1, and the top 30 DEGs at 14 days compared to 7 days were shown in Table 2. However, there were no statistically significant DEGs between day 21 and day 14 pellet, suggesting the induced chondrogenesis was halted (Fig. 2C & D).

Fig. 2
figure 2

Volcano plot showing DEGs. A DEGs between MSC and C7 group. B DEGs between C7 and C14 group. C DEGs between C14 and C21 group. D DEGs between C21 and CHO group

Table 1 Top 30 markedly differentially expressed genes between MSC and 7 days after induced chondrogenesis
Table 2 Top 30 markedly differentially expressed genes between 7 and 14 days after induced chondrogenesis

Functional GO terms and pathway enrichment analysis

The GO functional annotation analysis of the 6247 DEGs at day 7 revealed that extracellular matrix organization (P = 9.13E-11), extracellular structure organization (P = 1.22E-10), cell-substrate adhesion (P = 2.65E-09), cell junction organization (P = 5.54E-09), and ossification (P = 3.89E-08) were the top 5 altered biological processes (Fig. 3A). The KEGG functional analysis of these DEGs at day 7 revealed that focal adhesion (P = 9.88E-09), arrhythmogenic right ventricular cardiomyopathy (P = 7.96E-06), regulation of actin cytoskeleton (P = 9.07E-06), hepatocellular carcinoma (P = 3.87E-05) and pathways in cancer (P = 4.74E-05) were the top 5 altered pathways (Fig. 3B).

Fig. 3
figure 3

GO and KEGG analysis of DEGs in day 7 and day 14. A & B The GO and KEGG analysis of the 6247 DEGs between MSC and C7 group. C & D The GO and KEGG analysis of the 85 DEGs between C7 and C14 group

The GO functional annotation analysis of 85 DEGs on day 14 was also analyzed. The top 5 terms of biological processes included neutrophil chemotaxis (P = 8.11E-08), neutrophil migration (P = 1.55E-07), granulocyte chemotaxis (P = 4.12E-07), granulocyte migration (P = 8.41E-07), and response to bacterium (P = 1.47E-06) (Fig. 3C). These biological processes are related to common genes, including CXCL8, CCL20, TREM1, RAC2, SAA1, RIPOR2, and ITGB2. The KEGG functional analysis of these DEGs at day 14 identified two pathways of significance: rheumatoid arthritis (P = 1.38E-04) and legionellosis (P = 4.35E-04) (Fig. 3D).

PPI network analysis

In order to identify the key DEGs during the whole process of chondrogenesis, the common DEGs in the 7 days chondrogenesis group and 14 days chondrogenesis group were identified and selected. A total of 65 DEGs were identified and subjected to STRING for PPI analysis. The results of the PPI network are shown (Fig. 4A). To eliminate unrelated genes and to identify hub node genes in the process, the results were further analyzed with Cytoscape and cytoHubba. The top 10 hub nodes were identified, including CXCL8, TLR2, CCL20, MMP3, BCL2A1, SAA1, NT5E, RAC2, IL24, and MEFV (Fig. 4B).

Fig. 4
figure 4

PPI network of co-DEGs between day 7 and day 14. A PPI network of 65 co-DEGs in C7 and C14 group. B The top 10 hub nodes of the PPI network as analyzed by cytoHubba

Identification of key transcription factors during chondrogenesis

Transcriptional factors (TFs) are crucial in initiating the chondrogenic process. Therefore, we screened the top 5 markedly differentially expressed TFs on day 7 (Table 3). Top 5 differentially expressed TFs includes LEF1 (P = 3.94E-48), FOXO1 (P = 6.59E-26), RORA (P = 4.46E-24), BHLHE41 (P = 6.68E-24), SOX5 (P = 3.45E-20). The TF differentially expressed on day 14 was also analyzed. FOXQ1 (P = 1.23E-07) and TBX3 (P = 1.77E-05) were identified as markedly differentially expressed TFs on day 14 (Table 4).

Table 3 Top 5 markedly differentially expressed transcriptional factors 7 days after induced chondrogenesis
Table 4 Markedly differentially expressed transcriptional factors 14 days after induced chondrogenesis

Validation on transcriptional factors involved in chondrogenesis

During MSCs chondrogenesis, the initial stage (before day 7) was regarded as the most critical stage in determining MSCs' fate [16]. Therefore, we focused on the TFs that are differentially expressed on day 7. The role of LEF1, FOXO1, and SOX5 in MSCs chondrogenesis has been identified [17,18,19]. Because our previous studies show that melatonin promotes chondrogenic differentiation [15, 20], we are particularly interested in RORA, the nuclear receptor for melatonin. Treatment throughout the induction process with an antagonist for RORA SR3335 (1 μM) attenuated the effect of induced chondrogenesis. While continued treatment with 1 μM SR1078, an agonist for RORA, promoted chondrogenesis (Fig. 5A & B). The expression of chondrogenesis marker genes COL2A1 and ACAN were upregulated upon SR1078 treatment after 7 and 14 days of induction (Fig. 5C). When treated with SR3335 for 7 days, the expression level of COL2A1 was suppressed, but no significant difference between the vehicle and SR3335 in 14 days (Fig. 5D).

Fig. 5
figure 5

RORA promotes MSCs chondrogenesis. A Cartilage pellets collected 7 and 14 days after inducing MSCs chondrogenesis. Indicated treatment was added to the culture media. The cartilage pellets were observed with stereoscopic microscope. B The diameter (mm) of cartilage pellet after 7 and 14 days of induced chondrogenesis, 6 pellets were used for measurement in each group. C RT-qPCR analysis was conducted to assess the mRNA expression of COL2A1 and ACAN in chondrocyte pellet with the addition of indicated treatments for 7 days. D RT-qPCR analysis was conducted to assess the mRNA expression of COL2A1 and ACAN in chondrocyte pellet with the addition of indicated treatments for 14 days.* p ≤ 0.05 and ** p ≤ 0.01 compared to the vehicle treated group


In this study, we identified the transcriptional profile during chondrogenesis of human MSCs. Recent studies have demonstrated that MSCs are promising therapies for treating musculoskeletal disorders [21, 22]. Inducing MSCs chondrogenesis is an especially promising strategy to treat cartilage defects because of its low immunogenicity and capability to repair cartilage defects in situ. However, the mechanism underlying different stages during chondrogenesis remained relatively unknown. Therefore, identifying key genes and pathways may provide new insight into optimizing conditions of chondrogenesis both in vitro and in vivo.

Cartilage regeneration is a promising disease-modifying OA therapy [23]. The cartilage regeneration process involves stem cell differentiation and ECM synthesis. However, the therapeutic application of induced chondrogenesis is limited, partially due to the declined differentiation capacity of MSCs. Osteoarthritis and cartilage defects create an inflammatory microenvironment, which makes MSCs aberrantly differentiate to the fibroblast phenotype [24]. Therefore, understanding the key regulatory factors can contribute to optimizing the chondrogenic method and microenvironment. For example, activating the pro-inflammatory IL-6/STAT-3 pathway promotes chondrogenesis, revealing that inflammation may play a dual role in regulating chondrogenesis [25]. However, overproduction and accumulation of the pro-inflammatory factor IL-1β hamper the chondrogenic differentiation of MSCs [14].In this study, we found that neutrophil-related cytokines are significantly upregulated during chondrogenesis. Similar to our findings, another study found that immune response-related DEGs, including CXCL8 and CCL20, are upregulated during chondrogenesis [26].

The PPI analysis suggested some important hub genes in the network, as shown in Fig. 3B, including CXCL8, TLR2, CCL20, MMP3, BCL2A1, SAA1, RAC2, and IL24. These genes are closely related to regulating inflammation and response to inflammatory cytokines. Interestingly, top DEGs and hub genes reveal some common genes, including CXCL8, CCL20, SAA1, and RAC2. CXCL8, a member of the CXC chemokine family, encodes the interleukin-8 (IL-8) protein. The role of IL-8 in chondrogenesis remained controversial. IL-8 stimulates chondrocyte growth and differentiation and may contribute to cartilage defect repair [27, 28]. In contrast, other studies suggested that IL-8 may promote cartilage destruction during the pathogenesis of osteoarthritis [29, 30]. Therefore, the role of CXCL8 in chondrogenesis remains unclear and requires further exploration. As for TLR2 (encoding Toll-like receptor 2), a study by Lee et al. demonstrated that TLR2 promotes chondrogenic differentiation and consequent calcification in vascular smooth muscle cells [31]. Moreover, studies found that TLR2 is expressed in MSCs, and treatment with TLR2 agonists promotes chondrogenesis [32]. Interestingly, stimulation with TLR2 agonists led to increased production of IL-8 in MSCs in a dose-dependent manner [33]. Therefore, TLR2 and CXCL8 may represent a new pathway in regulating chondrogenesis and are worth further exploration.

We also identified several critical TFs that participated in chondrogenesis, including LEF1, FOXO1, RORA, BHLHE41, and SOX5. LEF1 forms a transcriptional complex with β-catenin and acts as an important downstream effector of the Wnt/β-catenin signaling pathway. A study showed that LEF1 cooperates with WNT5A and WNT11, directs MSCs chondrogenesis, and prevents chondrocytes from hypertrophic differentiation [17]. Besides, FOXO1 knockdown attenuated chondrogenesis by regulating SOX9 expression [18]. The role of RORA in MSCs chondrogenesis has not been investigated. We found that RORA pharmacological inhibition attenuated MSCs chondrogenesis while activating RORA promoted chondrogenesis. Previous studies suggested that RORA plays a contributory role in the pathogenesis of osteoarthritis by regulating cholesterol metabolism and IL6/STAT3 pathway [34, 35]. RORA also regulates chondrocytes hypertrophy in tibia organ cultures isolated from E15.5 mice [36]. These results suggested that different from LEF1, RORA may promote chondrogenesis and chondrocyte hypertrophy simutaneously.

In conclusion, this study provided insight into the gene expression profile during MSCs chondrogenesis. The hub genes and key transcriptional factors were identified, which may assist the research for cartilage regeneration and the development of osteoarthritis in future studies.

Availability of data and materials

The datasets generated and/or analysed during the current study are available in the GEO repository, accession number GSE210984 (



Mesenchymal stem cells


Transcription factor


Fold change


Differentially expressed genes


Kyoto Encyclopedia of Genes and Genomes


Gene Ontology


Protein–protein interaction


RNA integrity number


  1. Correa D, Lietman SA. Articular cartilage repair: Current needs, methods and research directions. SEMIN CELL DEV BIOL. 2017;62:67–77.

    Article  PubMed  Google Scholar 

  2. Makris EA, Gomoll AH, Malizos KN, Hu JC, Athanasiou KA. Repair and tissue engineering techniques for articular cartilage. NAT REV RHEUMATOL. 2015;11(1):21–34.

    Article  PubMed  CAS  Google Scholar 

  3. Saito T, Fukai A, Mabuchi A, Ikeda T, Yano F, Ohba S, Nishida N, Akune T, Yoshimura N, Nakagawa T, et al. Transcriptional regulation of endochondral ossification by HIF-2alpha during skeletal growth and osteoarthritis development. NAT MED. 2010;16(6):678–86.

    Article  PubMed  CAS  Google Scholar 

  4. Allas L, Boumediene K, Bauge C. Epigenetic dynamic during endochondral ossification and articular cartilage development. Bone. 2019;120:523–32.

    Article  PubMed  CAS  Google Scholar 

  5. Rim YA, Nam Y, Ju JH. The Role of Chondrocyte Hypertrophy and Senescence in Osteoarthritis Initiation and Progression. Int J Mol Sci. 2020;21(7):2358.

    Article  PubMed Central  CAS  Google Scholar 

  6. Lian C, Wang X, Qiu X, Wu Z, Gao B, Liu L, Liang G, Zhou H, Yang X, Peng Y, et al. Collagen type II suppresses articular chondrocyte hypertrophy and osteoarthritis progression by promoting integrin beta1-SMAD1 interaction. BONE RES. 2019;7:8.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Song J, Li J, Yang F, Ning G, Zhen L, Wu L, Zheng Y, Zhang Q, Lin D, Xie C, et al. Nicotinamide mononucleotide promotes osteogenesis and reduces adipogenesis by regulating mesenchymal stromal cells via the SIRT1 pathway in aged bone marrow. CELL DEATH DIS. 2019;10(5):336.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Hyun SY, Lee JH, Kang KJ, Jang YJ. Effect of FGF-2, TGF-beta-1, and BMPs on Teno/Ligamentogenesis and Osteo/Cementogenesis of Human Periodontal Ligament Stem Cells. MOL CELLS. 2017;40(8):550–7.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  9. Park IS, Kim JH, Bae JS, Kim DK, Mo JH. The Supernatant of Tonsil-Derived Mesenchymal Stem Cell Has Antiallergic Effects in Allergic Rhinitis Mouse Model. Mediators Inflamm. 2020;2020:6982438.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Zhao Y, Xie L. An Update on Mesenchymal Stem Cell-Centered Therapies in Temporomandibular Joint Osteoarthritis. Stem Cells Int. 2021;2021:6619527.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Murphy MP, Koepke LS, Lopez MT, Tong X, Ambrosi TH, Gulati GS, Marecic O, Wang Y, Ransom RC, Hoover MY, et al. Articular cartilage regeneration by activated skeletal stem cells. Nat Med. 2020;26(10):1583–92.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Lin H, Sohn J, Shen H, Langhans MT, Tuan RS. Bone marrow mesenchymal stem cells: Aging and tissue engineering applications to enhance bone healing. Biomaterials. 2019;203:96–110.

    Article  PubMed  CAS  Google Scholar 

  13. Mushahary D, Spittler A, Kasper C, Weber V, Charwat V. Isolation, cultivation, and characterization of human mesenchymal stem cells. Cytometry A. 2018;93(1):19–31.

    Article  PubMed  CAS  Google Scholar 

  14. Gao B, Gao W, Wu Z, Zhou T, Qiu X, Wang X, Lian C, Peng Y, Liang A, Qiu J, et al. Melatonin rescued interleukin 1beta-impaired chondrogenesis of human mesenchymal stem cells. Stem Cell Res Ther. 2018;9(1):162.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Wu Z, Qiu X, Gao B, Lian C, Peng Y, Liang A, Xu C, Gao W, Zhang L, Su P, et al. Melatonin-mediated miR-526b-3p and miR-590-5p upregulation promotes chondrogenic differentiation of human mesenchymal stem cells. J Pineal Res. 2018;65(1): e12483.

    Article  PubMed  Google Scholar 

  16. Almalki SG, Agrawal DK. Key transcription factors in the differentiation of mesenchymal stem cells. Differentiation. 2016;92(1–2):41–51.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  17. Diederichs S, Tonnier V, Marz M, Dreher SI, Geisbusch A, Richter W. Regulation of WNT5A and WNT11 during MSC in vitro chondrogenesis: WNT inhibition lowers BMP and hedgehog activity, and reduces hypertrophy. Cell Mol Life Sci. 2019;76(19):3875–89.

    Article  PubMed  CAS  Google Scholar 

  18. Kurakazu I, Akasaki Y, Hayashida M, Tsushima H, Goto N, Sueishi T, Toya M, Kuwahara M, Okazaki K, Duffy T, et al. FOXO1 transcription factor regulates chondrogenic differentiation through transforming growth factor beta1 signaling. J Biol Chem. 2019;294(46):17555–69.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Park JS, Yi SW, Kim HJ, Kim SM, Kim JH, Park KH. Construction of PLGA Nanoparticles Coated with Polycistronic SOX5, SOX6, and SOX9 Genes for Chondrogenesis of Human Mesenchymal Stem Cells. ACS Appl Mater Interfaces. 2017;9(2):1361–72.

    Article  PubMed  CAS  Google Scholar 

  20. Gao W, Lin M, Liang A, Zhang L, Chen C, Liang G, Xu C, Peng Y, Chen C, Huang D, et al. Melatonin enhances chondrogenic differentiation of human mesenchymal stem cells. J PINEAL RES. 2014;56(1):62–70.

    Article  PubMed  CAS  Google Scholar 

  21. Liu F, Xu H, Huang H. A novel kartogenin-platelet-rich plasma gel enhances chondrogenesis of bone marrow mesenchymal stem cells in vitro and promotes wounded meniscus healing in vivo. Stem Cell Res Ther. 2019;10(1):201.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Richardson SM, Kalamegam G, Pushparaj PN, Matta C, Memic A, Khademhosseini A, Mobasheri R, Poletti FL, Hoyland JA, Mobasheri A. Mesenchymal stem cells in regenerative medicine: Focus on articular cartilage and intervertebral disc regeneration. Methods. 2016;99:69–80.

    Article  PubMed  CAS  Google Scholar 

  23. Cho Y, Jeong S, Kim H, Kang D, Lee J, Kang SB, Kim JH. Disease-modifying therapeutic strategies in osteoarthritis: current status and future directions. Exp Mol Med. 2021;53(11):1689–96.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Li M, Yin H, Yan Z, Li H, Wu J, Wang Y, Wei F, Tian G, Ning C, Li H, et al. The immune microenvironment in cartilage injury and repair. Acta Biomater. 2022;140:23–42.

    Article  PubMed  CAS  Google Scholar 

  25. Kondo M, Yamaoka K, Sakata K, Sonomoto K, Lin L, Nakano K, Tanaka Y. Contribution of the Interleukin-6/STAT-3 Signaling Pathway to Chondrogenic Differentiation of Human Mesenchymal Stem Cells. Arthritis Rheumatol. 2015;67(5):1250–60.

    Article  PubMed  CAS  Google Scholar 

  26. Zhou J, Li C, Yu A, Jie S, Du X, Liu T, Wang W, Luo Y. Bioinformatics analysis of differentially expressed genes involved in human developmental chondrogenesis. Medicine (Baltimore). 2019;98(27):e16240.

    Article  Google Scholar 

  27. Brittberg M, Sjogren-Jansson E, Thornemo M, Faber B, Tarkowski A, Peterson L, Lindahl A. Clonal growth of human articular cartilage and the functional role of the periosteum in chondrogenesis. Osteoarthritis Cartilage. 2005;13(2):146–53.

    Article  PubMed  CAS  Google Scholar 

  28. Yoon DS, Lee KM, Kim SH, Kim SH, Jung Y, Kim SH, Park KH, Choi Y, Ryu HA, Choi WJ, et al. Synergistic Action of IL-8 and Bone Marrow Concentrate on Cartilage Regeneration Through Upregulation of Chondrogenic Transcription Factors. Tissue Eng Part A. 2016;22(3–4):363–74.

    Article  PubMed  CAS  Google Scholar 

  29. Nees TA, Rosshirt N, Zhang JA, Reiner T, Sorbi R, Tripel E, Walker T, Schiltenwolf M, Hagmann S, Moradi B. Synovial Cytokines Significantly Correlate with Osteoarthritis-Related Knee Pain and Disability: Inflammatory Mediators of Potential Clinical Relevance. J Clin Med. 2019;8(9):1343.

    Article  PubMed Central  CAS  Google Scholar 

  30. van Geffen EW, van Caam AP, van Beuningen HM, Vitters EL, Schreurs W, van de Loo FA, van Lent PL, Koenders MI, Blaney DE, van der Kraan PM. IL37 dampens the IL1beta-induced catabolic status of human OA chondrocytes. Rheumatology (Oxford). 2017;56(3):351–61.

    Google Scholar 

  31. Lee GL, Yeh CC, Wu JY, Lin HC, Wang YF, Kuo YY, Hsieh YT, Hsu YJ, Kuo CC. TLR2 Promotes Vascular Smooth Muscle Cell Chondrogenic Differentiation and Consequent Calcification via the Concerted Actions of Osteoprotegerin Suppression and IL-6-Mediated RANKL Induction. Arterioscler Thromb Vasc Biol. 2019;39(3):432–45.

    Article  PubMed  CAS  Google Scholar 

  32. Pevsner-Fischer M, Morad V, Cohen-Sfady M, Rousso-Noori L, Zanin-Zhorov A, Cohen S, Cohen IR, Zipori D. Toll-like receptors and their ligands control mesenchymal stem cell functions. Blood. 2007;109(4):1422–32.

    Article  PubMed  CAS  Google Scholar 

  33. Kim HS, Shin TH, Yang SR, Seo MS, Kim DJ, Kang SK, Park JH, Kang KS. Implication of NOD1 and NOD2 for the differentiation of multipotent mesenchymal stem cells derived from human umbilical cord blood. PLoS ONE. 2010;5(10):e15369.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Choi WS, Lee G, Song WH, Koh JT, Yang J, Kwak JS, Kim HE, Kim SK, Son YO, Nam H, et al. The CH25H-CYP7B1-RORalpha axis of cholesterol metabolism regulates osteoarthritis. Nature. 2019;566(7743):254–8.

    Article  PubMed  CAS  Google Scholar 

  35. Liang T, Chen T, Qiu J, Gao W, Qiu X, Zhu Y, Wang X, Chen Y, Zhou H, Deng Z, et al. Inhibition of nuclear receptor RORalpha attenuates cartilage damage in osteoarthritis by modulating IL-6/STAT3 pathway. CELL DEATH DIS. 2021;12(10):886.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. Woods A, James CG, Wang G, Dupuis H, Beier F. Control of chondrocyte gene expression by actin dynamics: a novel role of cholesterol/Ror-alpha signalling in endochondral bone growth. J Cell Mol Med. 2009;13(9B):3497–516.

    Article  PubMed  Google Scholar 

Download references


We thank all clinical staff of Department of Orthopedics of Sun Yat-sen Memorial Hospital for assisting this research.


This research was supported by grants from the Guangdong Natural Science Foundation (No. 2019A1515110122 & No. 2021A1515012619), National Natural Science Foundation of China (No. 82072473 & No. 82202653).

Author information

Authors and Affiliations



Tongzhou Liang, Pengfei Li performed induced chondrogenesis induction. Anjing Liang, Yuanxin Zhu, Xianjian Qiu performed RNA extraction and sequencing. Jincheng Qiu, Yan Peng and Dongsheng Huang analyzed the sequencing data. Tongzhou Liang, Dongsheng Huang, Wenjie Gao and Bo Gao revised and finalized the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Dongsheng Huang, Wenjie Gao or Bo Gao.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Ethical Committee of Sun Yat-sen University, and written informed consent was obtained from all subjects included in the study. All methods were carried out in accordance with relevant guidelines and regulations.

Consent for publication

Not applicable.

Competing interests

The authors declare no potential conflicts of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Liang, T., Li, P., Liang, A. et al. Identifying the key genes regulating mesenchymal stem cells chondrogenic differentiation: an in vitro study. BMC Musculoskelet Disord 23, 985 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Mesenchymal Stem Cells
  • Chondrogenesis
  • Transcriptome Analysis
  • Transcriptional Factors