Skip to main content

Long-term dynamic compression enhancement TGF-β3-induced chondrogenesis in bovine stem cells: a gene expression analysis

Abstract

Background

Bioengineering has demonstrated the potential of utilising mesenchymal stem cells (MSCs), growth factors, and mechanical stimuli to treat cartilage defects. However, the underlying genes and pathways are largely unclear. This is the first study on screening and identifying the hub genes involved in mechanically enhanced chondrogenesis and their potential molecular mechanisms.

Methods

The datasets were downloaded from the Gene Expression Omnibus (GEO) database and contain six transforming growth factor-beta-3 (TGF-β3) induced bovine bone marrow-derived MSCs specimens and six TGF-β3/dynamic-compression-induced specimens at day 42. Screening differentially expressed genes (DEGs) was performed and then analysed via bioinformatics methods. The Database for Annotation, Visualisation, and Integrated Discovery (DAVID) online analysis was utilised to obtain the Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway enrichment. The protein-protein interaction (PPI) network of the DEGs was constructed based on data from the STRING database and visualised through the Cytoscape software. The functional modules were extracted from the PPI network for further analysis.

Results

The top 10 hub genes ranked by their connection degrees were IL6, UBE2C, TOP2A, MCM4, PLK2, SMC2, BMP2, LMO7, TRIM36, and MAPK8. Multiple signalling pathways (including the PI3K-Akt signalling pathway, the toll-like receptor signalling pathway, the TNF signalling pathway, and the MAPK pathway) may impact the sensation, transduction, and reaction of external mechanical stimuli.

Conclusions

This study provides a theoretical finding showing that gene UBE2C, IL6, and MAPK8, and multiple signalling pathways may play pivotal roles in dynamic compression-enhanced chondrogenesis.

Background

Cartilage has little self-renewable ability due to its instinctive physiologies [1, 2], which include an avascular, aneural and non-lymphatic system [3], low cellularity in adult tissue, and a dense hydrated extracellular matrix, hampering resident chondrocytes or progenitor cells migration to the defect site to secrete a reparative matrix [2]. Mesenchymal stem cells (MSCs) are promising cell sources for osteochondral engineering. Numerous studies have demonstrated successful induction of chondrogenesis in various biomaterials. This strategy shows remarkable potential in repairing cartilage defects caused by osteoarthritis and athletic injuries [4, 5]. The most commonly used chondrogenic medium contains the TGF-β superfamily, which is a crucial mediator of MSCs chondrogenesis. Literature shows that TGF-β has proven a success in inducing chondrogenesis in vitro [6, 7]. However, TGF-β-induced chondrocytes alone were then witnessed a hypertrophic phenotype [8], which is not an ideal cell phenotype. Thus, inhibiting hypertrophy during the chondrogenic process in vitro and maintaining a stable cartilaginous phenotype need to be overcome.

Inspired by the physiology of native articular cartilage subjected to the dynamic joint environment, mainly under compression and shearing conditions [9], the significance of biomechanical stimuli has been well-established in the case of cartilage. Previous studies have shown that the ligand-integrin-cytoskeleton complex is the major mechanosensing component of the cell. The dynamic load and integrin activate the focal adhesion kinase (FAK) and mitogen-activated protein kinase (MAPK) pathways, increase the intracellular calcium, and induce further cell processes [10, 11]. Additionally, there are other pathways that do not rely on calcium. Dynamic compression is the most highly used physical condition to promote chondrogenesis [12]. Dynamic compression has been proven not only to enhance the efficiency of growth factors, but also play an important role in maintaining chondrocytes phenotypes and inhibiting hypertrophy. Despite increasing research on the impact of mechanical stimuli on chondrogenesis, there is no comprehensive understanding of underlying genes, while signal pathways remain elusive. Hence, in order to develop an optimal chondrogenic differentiation strategy, there is a pressing need to identify the key genes and signal pathways involved.

Microarray technology provides a powerful tool for exploring the gene regulation pattern and molecular mechanisms involved in mechanical-enhanced chondrogenesis. It enables to investigate thousands of gene expression patterns [13]. The microarray data can be uploaded and shared through open-source databases such as the Gene Expression Omnibus (GEO) database [14]. Huang et al. [15] provided the first study of how long-term (21 days) dynamic compression affected chondrogenesis. They briefly displayed a preliminary microarray screen for the genome expression profiles with chondrogenic induction and long-term dynamic compression. With limited data currently available on this topic, this study was conducted based on selected Huang’s data on gene expression patterns affected by dynamic compression after a sustained TGF-β3 chondrogenic induction of MSCs, and further analyses shed more in-depth understanding of the underlying mechanisms. Datasets were downloaded containing genes expression data between TGF-β3-induced and TGF-β3/dynamic-compression-induced chondrogenesis of bovine MSCs from the GEO. Differentially expressed genes (DEGs) were screened, and Gene Ontology (GO), Kyoto Encyclopaedia of Genes and Genomes (KEGG), and protein-protein interaction (PPI) network analyses were performed to explore the hub genes and key modules involved. To summarise, 236 DEGs and 10 hub genes were identified, which may be key candidates for responding to dynamic compression during chondrogenic differentiation of MSCs.

Results

Data pre-processing and identification of DEGs

Figure 1 displays the gene expression data of two groups containing 12 samples after normalisation. Medians show good alignment, indicating a high data quality after normalisation and suitability for the following analyses. The Volcano plot (Fig. 2) demonstrates the differential expression status of all detected genes highlighting DEGs beyond the set cut-off criterion. A total of 236 DEGs were obtained, of which 178 (75.42%) were up-regulated genes, and 58 (24.58%) were down-regulated genes in TGF-β3/dynamic-compression-induced MSCs compared to TGF-β3-induced MSCs. The cluster heatmap of DGEs is displayed in Fig. 3. The Euclidean distance was adopted to cluster the genes and produce the dendrograms. The red and green colours distinguish relatively higher or lower gene expression in each sample. Significant differences in DEGs expression patterns can be observed between these two groups (with and without dynamic compressive stimulation), indicating that the DEGs are reliable and eligible for the following analyses. The top 10 most significantly up-regulated and down-regulated genes are shown in Table 1.

Fig. 1
figure1

Box-plot of normalised data. Black lines in the boxes represent medians

Fig. 2
figure2

Volcano plot of all genes detected in the microarray. Each dot represents a gene. Dashed lines divide areas of down- and up-regulated genes. The X-axis is log2-base fold change, and Y-axis is −log10-base adjusted P-value

Fig. 3
figure3

Cluster heatmap demonstrates hierarchical clustering analysis results according to DEGs. Each row represents a DEG, and each column represents a sample. The colour displays the relative gene expression level. Green indicates lower values in gene expression, and red indicates higher values

Table 1 The top 10 most significantly up-regulated and down-regulated DEGs

GO and pathway enrichment analyses

GO enrichment and KEGG pathway enrichment analyses were performed to identify the biological function of DEGs. In GO terms, a negative regulation of angiogenesis, in utero embryonic development, and inflammatory responses provided the most significant enrichment in the biological process. The most significant enrichment in the cellular component was created through the cytoplasm, transcription elongation factor complex, and cortical actin cytoskeleton. Haemoglobin binding and ATP binding represented the most significant enrichment in the molecular function. A full list of enriched GO terms is shown in Table 2. In the KEGG pathway enrichment analysis, after screening and removing obviously irrelevant disease clusters, the PI3K-Akt signalling pathway, the toll-like receptor signalling pathway, and the TNF signalling pathway were remarkably enriched in dynamic compression-enhanced chondrogenesis (see Fig. 4 and Table 3).

Table 2 Significantly enriched GO terms of DEGs
Fig. 4
figure4

KEGG pathway enrichment analysis. The gradual colour stands for −log10-base adjusted P-value, red indicates a higher adjusted P-value, and green indicates a lower adjusted P-value. Dots size stands for gene count number. The X-axis represents the gene percentage ratio, and the Y-axis lays out pathway names

Table 3 Signalling pathway enrichment analysis of DEGs

PPI network construction

The PPI network of all DEGs (Fig. 5), constructed through the STRING database, includes 113 nodes and 185 edges. Among them, DEGs, IL6, UBE2C, TOP2A, MCM4, PLK2, SMC2, BMP2, LMO7, TRIM36, and MAPK8 were screened as hub genes, according to their connection degrees (Table 4). IL6 displayed the highest degree (= 14), followed by UBE2C (= 13). The deletion of IL6 and UBE2C remarkably loosens the structure of the PPI network as it reduces the interaction between proteins. Therefore, IL6 and UBE2C are the core nodes of PPI, suggesting that IL6 and UBE2C play an important role in responding to dynamic compression.

Fig. 5
figure5

PPI network of all DEGs. Red nodes with mesh patterns represent hub genes analysed by the cytoHubba. Node sizes reflect the connection degree, the higher degree is, the larger node size is

Table 4 The top 10 hub genes

Functional module analysis

The MCODE generated five sub-clusters, which reflect the high modularisation of a gene network. The top three amongst five modules contain nine of ten hub genes and are shown in Fig. 6. Module 1 consists of 14 nodes and 49 edges, and scores 7.54. Module 2 consists of 5 nodes and 10 edges, and scores 5.00. Module 3 consists of 4 nodes and 6 edges, and scores 4.00. As for annotation, this study focussed on Modules 1 and 3, which had the engagement of hub genes. Genes in Module 1 were mostly classified into GO terms of protein polyubiquitination, nuclear chromosome, and ATP binding, while genes in Module 3 were mainly classified into GO terms of defence responses to the virus, nucleus and cytokine activity (see Table 5). After screening and removing obvious irrelevant disease clusters, genes in Module 1 were mainly enriched through the ubiquitin-mediated proteolysis pathway, while the toll-like receptor signalling pathway, NOD-like receptor signalling pathway, cytosolic DNA-sensing pathway, and RIG-I-like receptor signalling pathway were identified for genes in Module 3 (see Table 6).

Fig. 6
figure6

The top three most significant modules. Red nodes with mesh patterns represent hub genes analysed by the cytoHubba. Node sizes reflect the connection degree. The higher connection degree is, the larger node size is

Table 5 Top 5 significantly enriched GO terms of module 1 and 3
Table 6 Signalling pathway enrichment analysis of module 1 and 3

Discussion

Chondrocytes respond to mechanical stimuli through regulating gene expression, proliferation, and metabolic functions. However, little is known about the key genes, signalling pathways, and proteins. Chondrocytes have been considered a post-mitotic tissue with nearly no cellular turnover. They are surrounded by an extracellular matrix comprised of glycosaminoglycan (GAG) and collagen and are subjected to daily dynamic compression. During the in vitro culture, growth factors such as bone morphogenetic protein (BMP) and the TGF-β superfamily are indispensable for the chondrogenic differentiation of MSCs [16]. However, compared to native cartilage, cartilage induced by TGF-β alone showed inferior mechanical properties [17]. Dynamic compression was proved to stabilise the chondrogenic phenotype by inhibiting hypertrophy in the presence of TGF-β3 [18]. To sum up, dynamic compression is essential for inducing non-hypertrophic chondrogenesis of MSCs.

Furthermore, in Huang’s [15] original study, the results revealed that the timing of applying dynamic compression was important. The loading initiated soon after MSCbeing encapsulated into agarose, led to reduced mechanical properties. In contrast, loading initiated after chondrogenic induction and ECM elaboration in the presence of TGF-β3, enhanced the mechanical properties of MSC-seeded constructs. This may be attributed to different mechanotransduction pathways between differentiated and undifferentiated MSCs. Following a shift from the 2% agarose to a denser, cartilage-like construct, the stresses induction was higher. The microarray analysis of the original study showed that several genes from the MMP/TIMP family were significantly modulated. However, the original microarray analysis merely took the fold change of genes into consideration when evaluated the gene importance. This may lead to an inadequate revelation of actual hub genes, as the fold change of genes is not always reliable and proportional to the actual influence on cells. Considering the availability of original data, and the fact that dynamic loading with TGF-β3 is the proven condition that promoted a stable chondrogenic phenotype, this study was built up on one of Huang’s series experiments for further bioinformatics analysis. It explored how compressive stimuli influence the gene expression after chondrogenic induction using TGF-β3, to shed important insight on the mechanism behind. Although the study was initially intended to collect a series of datasets at different time points, the uploaded datasets involving mechanical loading were only available at the time point of day 42. As consequence, a possible loss of some gene information at the initial time point might become inevitable, nevertheless, the long-term gene modulation data at the ending time point was indispensable for analysis. New understanding resulting from the data excavation may contribute towards developing a better strategy to enhance chondrogenic efficiency, quality, and stability.

The high-throughput microarray technology combined with bioinformatics analysis has been widely used in providing new insight into gene expression changes and molecular mechanisms. In the present study, the GEO database was utilised to obtain microarray raw data. A total of 236 DEGs were identified between TGF-β3-induced and TGF-β3/dynamic-compression-induced MSCs, including 178 up-regulated genes and 58 down-regulated genes. After that, the DEGs were analysed by GO functional enrichment analysis and classified into three groups, which were subsequently further clustered, based on functions and signalling pathways.

The results of GO functional enrichment analysis showed that the DEGs were mainly enriched in the GO terms of inflammatory response, in utero embryonic development and negative regulation of angiogenesis. This conforms to previous studies showing that the inflammatory response was involved in chondrogenic regulation. Inflammatory factors have been recognised as an important driving force leading to cartilage breakdown, and their down-regulation is vital for constructing initial collagen networks. A previous animal study revealed that the three-day cyclic compression of 0.5 MPa at 0.5 Hz on bovine chondrocytes counteracted the cartilage degradation induced by IL-1 [19]. Therefore, dynamic loading is not only a stimulator for chondrogenesis, but also an anti-inflammatory factor against pro-inflammatory cytokines. In this study, there were two other GO terms – GO:0001701 (in utero embryonic development) and GO:0016525 (negative regulation of angiogenesis) – that were significantly abnormal between the TGF-β3-induced and TGF-β3/dynamic-compression-induced MSCs. This demonstrates that dynamic compression may affect the anatomical structure development of chondrogenesis. During the early embryogenesis and cartilage maturation, various mechanical stimuli in the microenvironment promote chondrogenesis and limb formation and are responsible for adult chondrocyte phenotype maintenance [20]. Generally, biomechanics has been widely regarded as a promoter of angiogenesis and osteogenesis [21, 22]. On the other hand, cartilage is an avascular system [3], however, the understanding regarding how cartilage maintains avascularity under a mechanical load is limited in the literature, and the underlying biomechanics have not yet been fully established. This study suggests that appropriate mechanical stimuli are vital for inducing less angiogenesis.

Moreover, KEGG pathways enrichment analysis was performed. Because the KEGG database integrates data on genomes, chemical molecules and biochemical systems, including pathways, drug, disease, gene sequences, and genomes, some irrelevant disease clusters might be unexpectedly enriched. These disease-related clusters were screened and removed from the results and discussion. The KEGG pathway enrichment of DEGs and module analysis showed that the PI3K-Akt signalling pathway, toll-like receptor signalling pathway and TNF signalling pathway were highly enriched. Studies have demonstrated that the activation of the PI3K-Akt pathway promotes the terminal differentiation of chondrocytes and inhibits the hypertrophic differentiation of chondrocytes [23, 24]. The toll-like receptors mainly use MyD88-dependent signalling to activate NF-κB to transcript pro-inflammatory cytokines. Moreover, the activation of the toll-like receptor-2 induces the chondrogenic differentiation of MSCs [25, 26]. On the other hand, the mechanical load may promote chondrogenesis by inhibiting the TNF signalling pathway to reduce cartilage degradation. Further investigation is desired to support these findings. In brief, the findings of identified GO terms and the KEGG pathways may provide a theoretical basis on how dynamic compression regulates chondrogenesis.

The PPI network was constructed to predict the connections of proteins encoded by DEGs. The top 10 hub genes were screened according to connection degree as follows: IL6, UBE2C, TOP2A, MCM4, PLK2, SMC2, BMP2, LMO7, TRIM36, and MAPK8. Nine of them functioned in two of the top three most significant modules, suggesting that these genes play a more important role in chondrogenesis and are enhanced by dynamic compression. The Modules 1 and 3 were extracted from the PPI network. UBE2C, TOP2A, MCM4, PLK2, SMC2 LMO7, and TRIM36 were contained in Module 1, which were mainly enriched in GO terms related to the cellular metabolic process. These genes have closed relationships with the cell cycle and proliferation, and some of them were found overexpressed in various tumours. Moreover, UBE2C [27], TOP2A [28] and MCM4 [29] were identified as DEGs in OA. However, to the best of our knowledge, there is as yet no study on how these genes function in MSCs differential regulation were enhanced by mechanical load. This needs further investigation.

It was reported that the downregulation of PLK2 inhibited the degree of inflammation of knee joint synovial tissue and inhibited the cartilage collagen destruction in rats [30]. In recent years, studies have revealed that the SMC family might regulate bone development via mitogenic signals and the Wnt pathway, which is a central pathway in the bone and cartilage differentiation [31]. However, little is known on the specific function of SMC2 in response to mechanical stimuli, which requires further study. LMO7 and TRIM36 are both cell cycle-related genes. The overexpression of TRIM36 decelerates the cell cycle and attenuates cell growth [32], however, their functions in chondrogenesis have not been identified. The IL6 and MAPK8 showed vital roles in Module 3, which GO terms were mainly enriched in response to stimuli and the immune system. The pro-inflammatory cytokine IL6 constitutes an important factor involved in inflammation, immunoregulation, haematopoiesis and tumorigenesis. Its function in chondrogenesis remains controversial. Some studies reported that IL6 inhibited the chondrogenic differentiation [33, 34], while others demonstrated that activating the IL6/STAT3 signalling pathway promoted homeostasis maintenance and cartilage regeneration [35]. It is speculated that mechanical stimulus within the appropriate range of intensity, duration, and frequency may function as a potent anti-inflammatory signal and impose a positive influence on chondrogenesis, while overloading and unloading may lead to cartilage degradation. MAPK8 belongs to the c-Jun N-terminal kinase (JNK), a family which is one of the three main categories of MAPK families. JNK activation represents a protective response to external stimuli. Mechanical stress may activate the JNK pathway by phosphorylating ERK1/2, p38 MAPK, and SAPK/ERK kinase-1 (SEK1), resulting in chondrogenic differentiation [36] and apoptosis regulation [37]. Collectively, the comprehensive findings from this study show that UBE2C, IL6, and MAPK8 may play more important roles in dynamic compression enhanced chondrogenesis, unlike the original study which suggested the MMP/TIMP family might be the key genes (15).

Conclusions

This study analysed the gene expression profiles between TGF-β3-induced and TGF-β3/dynamic-compression-induced MSCs using a bioinformatics approach. 236 DEGs were found and annotated into GO terms and KEGG pathways, followed by constructing a PPI network and module mining. To our knowledge, this is the first time that genes, including UBE2C, IL6, and MAPK8, are identified to play a pivotal role in dynamic compression enhanced chondrogenesis via regulating proliferation, apoptosis and inflammatory response. Multiple signalling pathways, including the PI3K-Akt signalling pathway, toll-like receptor signalling pathway, TNF signalling pathway, and MAPK pathway, may be involved in sensation, transduction, and reaction of external mechanical stimuli. Although this is the first study giving a comprehensive genetic perspective on the interaction between mechanical stress and chondrogenesis, more experimental evidences are required to verify these findings. Further experimental studies are planned confirm these analysis results, which will be featured in the near future.

Methods

Microarray data information

The gene expression profiles of GSE18879 were downloaded from a public functional genomics data repository GEO database (https://www.ncbi.nlm.nih.gov/geo) [14] with the platform GPL2112 [Bovine] Affymetrix Bovine Genome Array (Affymetrix Inc., Santa Clara, CA, USA) [15]. This dataset includes negative control, TGF-β3-induced and TGF-β3/dynamic-compression-induced bovine bone marrow-derived MSCs specimens at three time points – day 3, 21 and 42 (repeated six times for each one). For specific groups, 10 ng/mL TGF-β3 was applied throughout 42 days, and the 10% strain dynamic compression at 1 Hz for 4 h daily began from day 21 onwards. Among them, the arrays of TGF-β3-induced and TGF-β3/dynamic-compression-induced specimens at day 42 were selected for analysis.

Data pre-processing

The CEL format files of raw data were converted into probe expression matrix, then underwent background adjustment, quantile normalisation, and ssummarisation using the Robust Multichip Average (RMA) in the RMAExpress software (version 1.2.0) [38]. Then, a log2 transformation was performed on the gene expression levels when the expression matrix was exported. After that, the probe serial numbers were transformed into official gene symbols.

Identification of DEGs

The up-regulated and down-regulated DEGs between TGF-β3-induced MSCs specimens and TGF-β3/dynamic-compression-induced MSCs specimens were identified through the Limma package on the NetworkAnalyst 3.0 web tool (https://www.networkanalyst.ca), which is a visual analytics platform for comprehensive gene expression profiling and meta-analysis [39]. Moreover, the p-value was corrected using the Benjamini-Hochberg test. Finally, the cut-off criterion of DEGs was set at the log2 fold change |log2FC| > 1.5 and adjusted as P < 0.05.

GO and pathway enrichment analyses

The Database for Annotation, Visualisation, and Integrated Discovery (DAVID, https://david.ncifcrf.gov) is an online functional enrichment analysis web tool that provides systematic annotation information for the biological function of large-scale gene list [40, 41]. In this study, GO enrichment and KEGG pathway enrichment analyses of DEGs were performed using DAVID with a cut-off criterion of gene count > 2 and P < 0.05. The GO analysis comprises of biological processes (BP), cellular components (CC), and molecular functions (MF). Irrelevant disease clusters in the KEGG pathway enrichment analysis were screened and removed before analysis and discussion.

PPI network construction

In order to understand the molecule mechanism and to study the interactions between dynamic compression and chondrogenesis, and between proteins encoded by DEGs and different proteins, the STRING (https://string-db.org) database [42] was utilised to recover the predicted associations between proteins encoded by DEGs and other proteins. The confidence score of > 0.4 was defined as significant. The results of the interaction data were then imported into the Cytoscape software (version 3.8.0) to visualise the PPI network. The degree distribution was established by counting the number of connections between different proteins in the network. The plug-in cytoHubba was utilised to screen the top 10 hub genes, ranked by degree.

Functional module analysis

Another built-in APP Molecular Complex Detection (MCODE) was utilised to detect the dense functionally connected sub-clusters within the large PPI network. The parameters of network scoring and cluster finding were set as follows: degree cutoff = 2, node score cutoff = 0.2, k-core = 2, and max depth = 100. The top three sub-clusters identified by modularity analysis were then selected for GO and pathway enrichment analysis via DAVID (gene count > 2 and P < 0.05). Similarly, irrelevant disease clusters in the KEGG pathway enrichment analysis were screened and removed before analysis and discussion.

Availability of data and materials

The datasets generated and/or analysed during the current study are available in the GEO repository, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE18879.

Abbreviations

MSCs:

Mesenchymal Stem Cells

GEO:

Gene Expression Omnibus

TGF-β3:

Transforming Growth Factor-Beta-3

DEGs:

Differentially Expressed Genes

DAVID:

Database for Annotation, Visualisation, and Integrated Discovery

GO:

Gene Ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

PPI:

Protein-Protein Interaction

FAK:

Focal Adhesion Kinase

MAPK:

Mitogen-Activated Protein Kinase

GEO:

Gene Expression Omnibus.

GAG:

Glycosaminoglycan.

BMP:

Bone Morphogenetic Protein.

JNK:

c-Jun N-terminal Kinase.

SEK1:

SAPK/ERK Kinase-1.

YAP:

Yes-Associated Protein.

TAZ:

Transcriptional Coactivator with PDZ-binding Domain.

PUU:

Poly (urea-urethane).

3D-TIPS:

3-Dimension Printing-guided Thermally Induced Phase Separation.

References

  1. 1.

    Vinatier C, Bouffi C, Merceron C, Gordeladze J, Brondello JM, Jorgensen C, Weiss P, Guicheux J, Noel D. Cartilage tissue engineering: towards a biomaterial-assisted mesenchymal stem cell therapy. Curr Stem Cell Res Ther. 2009;4(4):318–29. https://doi.org/10.2174/157488809789649205.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  2. 2.

    Roato I, Belisario DC, Compagno M, Lena A, Bistolfi A, Maccari L, Mussano F, Genova T, Godio L, Perale G, Formica M, Cambieri I, Castagnoli C, Robba T, Felli L, Ferracini R. Concentrated adipose tissue infusion for the treatment of knee osteoarthritis: clinical and histological observations. Int Orthop. 2019;43(1):15–23. https://doi.org/10.1007/s00264-018-4192-4.

    Article  PubMed  Google Scholar 

  3. 3.

    Ham O, Lee CY, Kim R, Lee J, Oh S, Lee MY, Kim J, Hwang KC, Maeng LS, Chang W. Therapeutic potential of differentiated Mesenchymal stem cells for treatment of osteoarthritis. Int J Mol Sci. 2015;16(7):14961–78. https://doi.org/10.3390/ijms160714961.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Zhang R, Ma J, Han J, Zhang W, Ma J. Mesenchymal stem cell related therapies for cartilage lesions and osteoarthritis. Am J Transl Res. 2019;11(10):6275–89.

    CAS  PubMed  PubMed Central  Google Scholar 

  5. 5.

    Zhou M, Lozano N, Wychowaniec JK, Hodgkinson T, Richardson SM, Kostarelos K, Hoyland JA. Graphene oxide: a growth factor delivery carrier to enhance chondrogenic differentiation of human mesenchymal stem cells in 3D hydrogels. Acta Biomater. 2019;96:271–80. https://doi.org/10.1016/j.actbio.2019.07.027.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Johnstone B, Hering TM, Caplan AI, Goldberg VM, Yoo JU. In vitro chondrogenesis of bone marrow-derived mesenchymal progenitor cells. Exp Cell Res. 1998;238(1):265–72. https://doi.org/10.1006/excr.1997.3858.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Tuli R, Tuli S, Nandi S, Huang X, Manner PA, Hozack WJ, Danielson KG, Hall DJ, Tuan RS. Transforming growth factor-beta-mediated chondrogenesis of human mesenchymal progenitor cells involves N-cadherin and mitogen-activated protein kinase and Wnt signaling cross-talk. J Biol Chem. 2003;278(42):41227–36. https://doi.org/10.1074/jbc.M305312200.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Pfeifer CG, Karl A, Kerschbaum M, Berner A, Lang S, Schupfner R, Koch M, Angele P, Nerlich M, Mueller MB. TGF-beta Signalling is suppressed under pro-hypertrophic conditions in MSC Chondrogenesis due to TGF-beta receptor Downregulation. Int J Stem Cells. 2019;12(1):139–50. https://doi.org/10.15283/ijsc18088.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Schatti O, Grad S, Goldhahn J, Salzmann G, Li Z, Alini M, Stoddart MJ. A combination of shear and dynamic compression leads to mechanically induced chondrogenesis of human mesenchymal stem cells. Eur Cell Mater. 2011;22:214–25. https://doi.org/10.22203/eCM.v022a17.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Sjaastad MD, Lewis RS, Nelson WJ. Mechanisms of integrin-mediated calcium signaling in MDCK cells: regulation of adhesion by IP3- and store-independent calcium influx. Mol Biol Cell. 1996;7(7):1025–41. https://doi.org/10.1091/mbc.7.7.1025.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Pommerenke H, Schmidt C, Durr F, Nebe B, Luthen F, Muller P, Rychly J. The mode of mechanical integrin stressing controls intracellular signaling in osteoblasts. J Bone Miner Res. 2002;17(4):603–11. https://doi.org/10.1359/jbmr.2002.17.4.603.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Grad S, Eglin D, Alini M, Stoddart MJ. Physical stimulation of chondrogenic cells in vitro: a review. Clin Orthop Relat Res. 2011;469(10):2764–72. https://doi.org/10.1007/s11999-011-1819-9.

    Article  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Lockhart DJ, Winzeler EA. Genomics, gene expression and DNA arrays. Nature. 2000;405(6788):827–36. https://doi.org/10.1038/35015701.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Clough E, Barrett T. The gene expression omnibus database. Methods Mol Biol. 2016;1418:93–110. https://doi.org/10.1007/978-1-4939-3578-9_5.

    Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Huang AH, Farrell MJ, Kim M, Mauck RL. Long-term dynamic loading improves the mechanical properties of chondrogenic mesenchymal stem cell-laden hydrogel. Eur Cell Mater. 2010;19:72–85. https://doi.org/10.22203/eCM.v019a08.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Wang W, Rigueur D, Lyons KM. TGFbeta signaling in cartilage development and maintenance. Birth Defects Res C Embryo Today. 2014;102(1):37–51. https://doi.org/10.1002/bdrc.21058.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Chiang H, Hsieh CH, Lin YH, Lin S, Tsai-Wu JJ, Jiang CC. Differences between chondrocytes and bone marrow-derived chondrogenic cells. Tissue Eng Part A. 2011;17(23–24):2919–29. https://doi.org/10.1089/ten.tea.2010.0732.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Aisenbrey EA, Bilousova G, Payne K, Bryant SJ. Dynamic mechanical loading and growth factors influence chondrogenesis of induced pluripotent mesenchymal progenitor cells in a cartilage-mimetic hydrogel. Biomater Sci. 2019;7(12):5388–403. https://doi.org/10.1039/C9BM01081E.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Torzilli PA, Bhargava M, Chen CT. Mechanical loading of articular cartilage reduces IL-1-induced enzyme expression. Cartilage. 2011;2(4):364–73. https://doi.org/10.1177/1947603511407484.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Responte DJ, Lee JK, Hu JC, Athanasiou KA. Biomechanics-driven chondrogenesis: from embryo to adult. FASEB J. 2012;26(9):3614–24. https://doi.org/10.1096/fj.12-207241.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Hankenson KD, Dishowitz M, Gray C, Schenker M. Angiogenesis in bone regeneration. Injury. 2011;42(6):556–61. https://doi.org/10.1016/j.injury.2011.03.035.

    Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Geris L, Vandamme K, Naert I, Vander Sloten J, Van Oosterwyck H, Duyck J. Mechanical loading affects angiogenesis and osteogenesis in an in vivo bone chamber: a modeling study. Tissue Eng Part A. 2010;16(11):3353–61. https://doi.org/10.1089/ten.tea.2010.0130.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Bai Y, Gong X, Dou C, Cao Z, Dong S. Redox control of chondrocyte differentiation and chondrogenesis. Free Radic Biol Med. 2019;132:83–9. https://doi.org/10.1016/j.freeradbiomed.2018.10.443.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Ulici V, Hoenselaar KD, Gillespie JR, Beier F. The PI3K pathway regulates endochondral bone growth through control of hypertrophic chondrocyte differentiation. BMC Dev Biol. 2008;8(1):40. https://doi.org/10.1186/1471-213X-8-40.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Najar M, Krayem M, Meuleman N, Bron D, Lagneaux L. Mesenchymal stromal cells and toll-like receptor priming: a critical review. Immune Netw. 2017;17(2):89–102. https://doi.org/10.4110/in.2017.17.2.89.

    Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Shirjang S, Mansoori B, Solali S, Hagh MF, Shamsasenjan K. Toll-like receptors as a key regulator of mesenchymal stem cell function: an up-to-date review. Cell Immunol. 2017;315:1–10. https://doi.org/10.1016/j.cellimm.2016.12.005.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Huan X, Jinhe Y, Rongzong Z. Identification of pivotal genes and pathways in osteoarthritic degenerative meniscal lesions via bioinformatics analysis of the GSE52042 dataset. Med Sci Monit. 2019;25:8891–904. https://doi.org/10.12659/MSM.920636.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Dunn SL, Soul J, Anand S, Schwartz JM, Boot-Handford RP, Hardingham TE. Gene expression changes in damaged osteoarthritic cartilage identify a signature of non-chondrogenic and mechanical responses. Osteoarthr Cartil. 2016;24(8):1431–40. https://doi.org/10.1016/j.joca.2016.03.007.

    CAS  Article  PubMed Central  Google Scholar 

  29. 29.

    Yang J, Wang N. Genome-wide expression and methylation profiles reveal candidate genes and biological processes underlying synovial inflammatory tissue of patients with osteoarthritis. Int J Rheum Dis. 2015;18(7):783–90. https://doi.org/10.1111/1756-185X.12643.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Liu W, Zha Z, Wang H. Upregulation of microRNA-27a inhibits synovial angiogenesis and chondrocyte apoptosis in knee osteoarthritis rats through the inhibition of PLK2. J Cell Physiol. 2019;234(12):22972–84. https://doi.org/10.1002/jcp.28858.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Cohen-Zinder M, Karasik D, Onn I. Structural maintenance of chromosome complexes and bone development: the beginning of a wonderful relationship? Bonekey Rep. 2013;2:388.

    Article  Google Scholar 

  32. 32.

    Miyajima N, Maruyama S, Nonomura K, Hatakeyama S. TRIM36 interacts with the kinetochore protein CENP-H and delays cell cycle progression. Biochem Biophys Res Commun. 2009;381(3):383–7. https://doi.org/10.1016/j.bbrc.2009.02.059.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Wei H, Shen G, Deng X, Lou D, Sun B, Wu H, Long L, Ding T, Zhao J. The role of IL-6 in bone marrow (BM)-derived mesenchymal stem cells (MSCs) proliferation and chondrogenesis. Cell Tissue Bank. 2013;14(4):699–706. https://doi.org/10.1007/s10561-012-9354-9.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Nakajima S, Naruto T, Miyamae T, Imagawa T, Mori M, Nishimaki S, Yokota S. Interleukin-6 inhibits early differentiation of ATDC5 chondrogenic progenitor cells. Cytokine. 2009;47(2):91–7. https://doi.org/10.1016/j.cyto.2009.05.002.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    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. https://doi.org/10.1002/art.39036.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Zhang Y, Pizzute T, Pei M. A review of crosstalk between MAPK and Wnt signals and its impact on cartilage regeneration. Cell Tissue Res. 2014;358(3):633–49. https://doi.org/10.1007/s00441-014-2010-x.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Davis RJ. Signal transduction by the JNK group of MAP kinases. Cell. 2000;103(2):239–52. https://doi.org/10.1016/S0092-8674(00)00116-1.

    CAS  Article  PubMed  Google Scholar 

  38. 38.

    Barash Y, Dehan E, Krupsky M, Franklin W, Geraci M, Friedman N, Kaminski N. Comparative analysis of algorithms for signal quantitation from oligonucleotide microarrays. Bioinformatics. 2004;20(6):839–46. https://doi.org/10.1093/bioinformatics/btg487.

    CAS  Article  PubMed  Google Scholar 

  39. 39.

    Zhou G, Soufan O, Ewald J, Hancock REW, Basu N, Xia J. NetworkAnalyst 3.0: a visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic Acids Res. 2019;47(W1):W234–41. https://doi.org/10.1093/nar/gkz240.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Huang da W, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res 2009, 37(1):1–13.

  41. 41.

    Huang da W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 2009, 4(1):44–57.

  42. 42.

    Jensen LJ, Kuhn M, Stark M, Chaffron S, Creevey C, Muller J, Doerks T, Julien P, Roth A, Simonovic M, Bork P, von Mering C. STRING 8--a global view on proteins and their functional interactions in 630 organisms. Nucleic Acids Res. 2009;37(Database issue):D412–6. https://doi.org/10.1093/nar/gkn760.

    CAS  Article  PubMed  Google Scholar 

Download references

Acknowledgments

The authors acknowledge financial support by the Engineering and Physical Science Research Council, the United Kingdom (EPSRC grant no. EP/L020904/1, EP/M026884/1 and EP/R02961X/1).

Funding

Financial support by the Engineering and Physical Science Research Council, the United Kingdom (EPSRC grant no. EP/L020904/1, EP/M026884/1 and EP/R02961X/1).

Author information

Affiliations

Authors

Contributions

JC had made substantial contributions to the data the acquisition, analysis, and was a major contributor in writing the manuscript. LC and WS contributed to the original conception and ensure correct data analysis. JH and WS had substantively revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Wenhui Song.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Appendix

Appendix

Table 7 Softwares and websites used in this paper

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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) 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

Chen, J., Chen, L., Hua, J. et al. Long-term dynamic compression enhancement TGF-β3-induced chondrogenesis in bovine stem cells: a gene expression analysis. BMC Genom Data 22, 13 (2021). https://doi.org/10.1186/s12863-021-00967-2

Download citation

Keywords

  • Bioinformatics
  • Chondrogenesis
  • Enrichment analysis
  • Mechanical stimulation
  • Mesenchymal stem cells