Bioinformatic Analysis of MicroRNAs and Signaling Networks that Orchestrates the hMSC Osteogenic Differentiation

doi:10.21203/rs.3.rs-957426/v1

Background

The drive to understand how microRNAs regulate the hMSCs fate and their osteogenic differentiation has inspired stem cells biology and regenerative biology field development. Numerous experimental data shows crucial involvement of miRs in skeletal development in embryos, osteogenic differentiation, and maturation. However, molecular mechanisms of miRs action, in other words their target signaling pathways and transcriptional factors that specific drives osteogenic differentiation is far from being understood.

Results

With meta-analysis we identified miRs significantly involved in hMSCs osteogenic differentiation. Statistical analysis revealed significant trend of upregulation of let-7a, mir-21, mir-26a, mir-29b, mir-101, mir-143, mir-218 during hMSCs differentiation into osteoblast. And the opposite trend was shown for mir-17, mir-31, mir-138 and mir-222: their content was significantly lower during osteogenic differentiation. Using bioinformatics approaches we identified predictable genes-target for each miRs, analyzed signaling networks and biological process enriched by these genes.

Conclusions

Our date shows that miRNAs specifically involved in hMSC transition into osteogenic differentiation via microenvironment formation (i.e. let-7a, mir-17, mir-21, mir-29b and mir-101), TGF-β/BMP–SMAD dependent pathway (i.e. let-7a, mir-17, mir-21, mir-26a mir-101) and MAPK signaling pathway (i.e. let-7a, mir-21, mir-26a, mir-29b, mir-143 and mir17). Yap-dependent expression of osteogenic TF modulated by let-7a, mir-31mir-101, mir-138 and mir-222. We predicted that mir-17, mir-26a, mir-29b, mir-101, mir-138 and mir-222 are specifically involved in canonical Wnt signaling dependent osteogenesis as well as in osteoblast maturation together with let-7a, mir-29b and mir-218 which modulate AMPK signaling. Additionally, identified mir-101 is likely involved into osteoblast homeostasis via Hedgehog signaling.

Our data expands knowledge in the field of hMSCs fate and osteogenesis orchestration by miRs, points to pro-osteogenic and anti-osteogenic miRs and their potential molecular pathways.

Molecular Biology

Molecular Genetics

Bioinformatics

osteogenesis

miR

microRNA

differentiation

MSC

osteogenic lineage

cells signaling

Mesenchymal stem cells (MSCs) are characterized by high plasticity and their ability to self-renewal and differentiation into various types of cells. Both, bone mesenchymal stromal cells (BMSCs), as well as adipose-derived mesenchymal stem cells (ASCs), can become osteoblasts, chondrocytes, and adipocytes. However, this ability is regulated by many microenvironment factors, including growth factors and hormones, extracellular matrix [1, 2] and microRNAs [3-5].

MicroRNAs (miRs/miRNAs) are important posttranscriptional regulators of gene expression. Through the RNA interference pathway, miRs regulate negatively or positively the endogenous gene expression [6,7]. It is interesting that each miRNA has hundreds of targets, and vice versa a single mRNA may be regulated by diverse miRNAs. During last decades, the crucial function of some miRNAs in osteogenesis and bone formation was shown [8, 38]. Moreover, deregulation of these miRNAs has been associated with skeletal and bone malformations.

The osteogenic differentiation is a highly organized and complex process, which is orchestrated by several transcription factors including RUNX2, OSTERIX, BMP, and others [9-11]. The expression of these transcription factors is negatively or positively regulated by canonical Wnt, JAK/STAT, MAPK and BMP signaling pathways [10,12-15]. However, miRs are able to regulate mentioned signaling pathways activity or osteogenic transcriptional factors expression. Indeed, miR-125 targets VDR, ErbB2 and Osterix in mice, miR-133 recognize RANX2 and this way inhibit differentiation of osteoprogenitors [16-18]. Diminishing of osteogenic differentiation of ASCs by miR-140-5p was shown in rats, moreover authors demonstrate that miR-140-5p effects osteogenesis via its direct targets: TLR4 and BMP2 [19]. On the contrary, promotion of osteogenic differentiation was shown for miR-27a and miR-128, where both miRs activate Wnt signaling pathway through DKK2 targeting in rats BMSCs [20, 21]. Regulation of humans BMSCs differentiation into osteoblast via targeting of canonical Wnt signaling was demonstrated for miR-381 and miR-889 [22, 23] Moreover, authors found higher level of miR-889 in osteoporosis patients [23]. Thus, our understanding of miRs regulation of osteogenesis continues to expand, however role of miRs in hMSCs fate is far from understood. Thus, all published works are focusing on different miRs function and targets analysis in humans, rats, or mice mesenchymal stem cells. Considering the complex nature of miRs regulation (i.e., many targets for each known miR), identification of the molecular mechanisms of miRs in osteogenic signaling pathways and transcriptional factors regulation can be beneficial for the future clinical applications of MSCs.

In this study, we used publicly available experimental data. First, we conducted a meta-analysis and identified a list of miRNAs that significantly change during hMSCs differentiation into osteoblast. Using bioinformatic approaches, we predict the signaling pathways and genes targeted by identified miRNAs during osteogenic differentiation of human mesenchymal stem cells.

First, we have analyzed all available experimental data as described in the Materials and Methods section (Fig. 1). Basing on results of this meta-analysis, we identified a line of microRNAs that are highly likely to be up- or downregulated during osteogenic differentiation of hMSCs (Fig. 2). Further statistical analysis revealed significant trend to upregulation of let-7a, mir-21, mir-26a, mir-29b, mir-101, mir-143, mir-218 during hMSCs differentiation into osteoblast. At the same time, the opposite trend was shown for mir-17, mir-31, mir-138 and mir-222: their content in such conditions was significantly lower.

Then, using three miR target gene prediction databases (miRanda, miRDB, and Targetscan) we verified the gene-targets for each identified miRs. As a result, we identified 523 target genes for let-7a; 151 for mir-21; 510 for mir-26a; 557 for mir-29b; 464 for mir-101; 195 for mir-143; 502 for mir-218; 797 for mir-17; 169 for mir-31; 250 for mir-138 and 234 for mir-222, and presented them at Venn diagrams (Fig. 3a, 4a).

The results of Gene Ontology (GO) enrichment analysis of identified miRs target genes revealed the overrepresentation of genes involved into regulation of transcription from RNA polymerase II promoter, transcription DNA templated biological processes (Fig. 3b, 4b). However, as listed in Fig. 3b and 4b identified miRs were enriched at the same time in more than one biological processes, so, the overlapping has taken place.

Mentioned above stem cell-specific biological processes - stem cell differentiation, stem cell division, positive regulation of mesenchymal cell proliferation, somatic stem cell population maintenance and mesodermal cell differentiation, - were enriched with gene targets of let-7a, mir-17, mir-21, mir-29b, mir-101, mir-222 and mir-218, whereas in case of mir-21 the results of GO data analysis prove the negative effect on stem cell differentiation processes. Specific biological process for cells division, proliferation and apoptosis are greatly enriched by gene targets of mir-17, mir-26a, mir-29b, mir-101, mir-138 and mir-222. However, analysis of experimental data proves positive regulation of cell proliferation by mir-222, whilst there was no unambiguous effect shown for mir-17. Then, data analysis confirms the hypothesis of apoptosis affection related to high mir-26a expression, but mir-138 and mir-222 are both seem to regulate the apoptotic process negatively.

Further analysis reveals specific biological processes that correspond to bone development and morphogenesis, namely osteoblast differentiation, bone morphogenesis, bone development, endochondral ossification, and bone mineralization to be enriched by mir-218, mir-143, mir-101, mir-21 and mir-17 target genes. The obtained results shown mir-101 and mir-218 to affect almost all mentioned process of bone development and morphogenesis. Morphogenesis process that is extremely close to bone development also is enriched by target genes of analyzed miRs. GO enrichment analysis points that chondrocyte development and differentiation is regulated by mir-26a and mir-21, where latter has a pro-chondrogenic effect. Cartilage development process is enriched by mir-17 and mir-138 gene targets. Cartilage condensation and cartilage morphogenesis are driven by mir-29b and mir-21 targets genes, respectively. Hence, it is not surprising that the majority of these miRs involves in embryonic skeletal system development and limb formation (i.e., mir-29b, mir-101, mir-218 and mir-17). Additionaly, let-7a and mir-222 take part in embryonic skeletal system development and embryonic forelimb morphogenesis regulation.

It is worth noting that such biological process as extracellular matrix organization, collagen biosynthesis and fibril organization as well as collagen catabolic process is orchestrated by identified by us miRs that are expressed on higher level during hMSCs osteogenic differentiation (i.e., let-7a, mir-29b). The same set of miRs together with mir-101 was found to be associated with platelet-derived growth factor (PDGF) receptor signaling pathway regulation. Gathering of following Gene Ontology: beta-catenin binding (i.e. mir-101, mir-29, mir-26a, mir-222 and mir-17); beta-catenin-TCF complex assembly (mir-101); Wnt signaling pathway including calcium modulating pathway (mir-222) and canonical Wnt regulation (i.e. mir-26a, mir-17 and mir-138); SMAD binding (i.e. let-7a, mir-26a, mir-101 and mir-17); pathway-restricted SMAD protein phosphorylation (i.e. let-7a); SMAD protein signal transduction (mir-21 and mir-26a); BMP signaling pathway (mir-26a, mir-21 and mir-17) additionally supported the crucial function of analyzed miRs in hMSCs osteogenic differentiation. Moreover, we found another important signaling pathways enriched by analyzed miRs targets what points their involvement in wider number of biological process including hMSCs self-renewal, differentiation, division, chondro- and osteogenesis and bone tissue repair. Indeed, we found that specific biological processes such as wound healing is regulated by mir-26a, mir-17 and mir-138. Additionaly, some of miRs level are linked with regulation of macrophage activation and macrophage differentiation; cytokine and interleukin biosynthetic; inflammatory response regulation (Fig. 3b, 4b).

KEGG pathway analysis revealed 17 potential pathways significantly enriched by analyzed miRs target genes: TGF-beta signaling pathway, Wnt signaling pathway, mTOR signaling pathway, Jak-STAT signaling pathway, Signaling pathways regulating pluripotency of stem cells, Rap1 signaling pathway, FoxO signaling pathway, Insulin resistance and signaling, p53 signaling pathway, Cell cycle, Ras signaling pathway, Hippo signaling pathway, ErbB signaling pathway and MAPK signaling pathway (Fig. 3с, 4с). It is interesting that, cGMP-PKG signaling pathway, Calcium signaling pathway, Chemokine signaling pathway, Hedgehog signaling pathway, Platelet activation, ECM-receptor interaction pathway, VEGF signaling pathway, T cell receptor signaling pathway and B cell receptor signaling are greatly enriched by up-regulated miRs targets (Fig. 3с). Further analysis show autophagy, endocytosis, estrogen signaling pathway and osteoclast differentiation pathways were found to be targeted by drop expressed miRs during hMSC differentiation (Fig. 4с).

Differentiation of osteogenic lineage cells is highly organized biological process which underlies in skeletal elements developing, in embryo, bone tissue remodeling and healing in adult. Many signaling networks and transcription factors as well as co-factors gathered MCSs differentiation into mature osteocyte [26]. Numerous research also pointed out that mircoRNAs has crucial function in osteogenesis and bone tissue remodeling targeting the key signaling pathways and transcriptional factors (Tabl. 1). However, different paper show that different miR or set of miRs are driving the hMSCs osteogenic differentiation. On one hand, a different biological model has been used to explain this variety of results, on the other – nature of genes expression regulation by miRs. Indeed, bioinformatic analysis revealed that a single miRNA can regulate up to several hundred target genes and vice versa, the same 3′-UTR of a gene can be targeted by many different miRs [30, 31]. In our opinion, it is depending on miRs or its targets tissue specific expression, hence on cells type.

Table 1

MicroRNAs regulates osteogenic differentiation via transcription factors and signaling networks.

miRNAs	Targets	Expression level	Reference
Pro-osteogenic differentiation miRNAs
miR-27a	PPARγ, GREM1	Up	Rev [38]
miR-21	P-Akt, HIF-1α	Up	Rev [38]
miR-217	DKK1	Up	Rev [38]
miR-26a	Runx2, OC, GSK3β	Up	Rev [38]
miR-148a	IGF1	Down	Rev [38]
miR-200b	Cx43, VEGF-A	Down	Rev [38]
miR-335-5p	DKK1	Up	Rev [38]
miR-92a	Smad6 3′-UTR	Up	Rev [38]
miR-9	Runx2, ERK	Up	Rev [38]
miR-199b-5p	GSK3β	Up	Rev [38]
miR-26a/b miR-29b	CDK6, HDAC4	Up	[39]
miR-1249-5p	PDX1/PI3K/Akt	Up	[40]
Anti-osteogenic differentiation miRNAs
miR-124	Sp7	Up	Rev [38]
miR-23a	CXCL12, LRP5, Runx2 3′UTR	Up	Rev [38]
miR-23a-5p	MAPK13	Up	Rev [38]
miR-23b	Runx2 3′UTR	Up	Rev [38]
miR-214	BMP2 3′UTR, FGF, p-JNK, p-p38	Up	Rev [38]
miR-214-5p	COL4A1	Up	Rev [38]
miR-144-3p	Smad4	Up	Rev [38]
miR-383	Satb2	Up	Rev [38]
miR-96	Osterix	Up	Rev [38]
miR-183-5p	Hmox1	Up	Rev [38]
miR-217	Runx 2 3′UTR	Up	Rev [38]
miR-133a-3p	MEG3	Up	Rev [38]
miR-10	Runx 2 3′UTR	Up	Rev [38]
miR-503-5p	Runx2, ALP	Up	Rev [38]
miR-205	SATB2, Runx2	Up	Rev [38]
miR-141, miR-200a	SVCT2 3′UTR	Up	Rev [38]
miR-204	Runx2 3′UTR	Up	Rev [38]
miR-138	ALP, RUNX2	Up	Rev [38]
miR-338-3p	Runx2, Fgfr2	Up	Rev [38]
miR-31	Satb2	Up	Rev [38]
miR-138	FAK, ERK1/2, RUNX2	Up	Rev [38]
miR-27a	Sp7	Up	Rev [38]
miR-381	WNT5A, FZD3	Up	[22]
miR-889	WNT7A	Up	[23]
miR204	BMP2	Up	[41]

Here we provide meta-analysis of experimental data and identify the line of miRs which are significantly involved in hMCSs osteogenic differentiation. The results of our analysis demonstrate complex involvement of different miRs in osteogenic hMSCs differentiation and show that each step of osteoblasts differentiation is orchestrated by a few miRs via similar molecular mechanisms (Fig. 5). Vice versa, the same miRs regulate different signaling pathways, regulating osteogenesis.

With KEGG pathway and GO analysis we found that let-7a, mir-17, mir-21 and mir-101 in hMSC are likely to be involved in transforming growth factor β (TGF-β) signaling where mir-17 is predicted to inhibit the TGF-β cytoplasmic mediator activity. At the same time, mir-29b, let-7a and mir-101 were predicted to regulate the platelet-derived growth factor (PDGF). Aadditionally, mir-101 was predicted to target the fibroblasts growth factor (FGF). Mir-29b was predicted to activate the vascular endothelial growth factor (VEGF) and maintains the cellular response to TGF-β and FGF. Together up-regulated PDGF, TGF-β, FGF and VEGF create appropriate microenvironment for hMSCs to enter in osteogenic lineage differentiation, namely promote collagen synthesis and fibril organization. GO analysis additionally shows important regulatory role of let-7a and mir-29b in collagen synthesis and fibril organization. Having summarized our data, we suppose that balance between high level of let-7a, mir-21, mir-29b, mir-101 and low level of mir-17 is necessary for microenvironment formation and hMSCs turn on into osteogenic lineage (Fig. 5).

Widely known TGF-β/BMP activates SMAD-dependent and non-SMAD-dependent signaling which results in main osteogenic transcriptional factor Runx2 expression [12]. TGF-β–SMAD signaling raises MSCs proliferation, chemotaxis, and early differentiation of osteoprogenitors, where BMP–SMAD signaling regulate almost every stage of osteoblast differentiation and maturation [12]. During osteoprogenitors development BMP–SMAD also activates the Runx2 and promotes hMSCs differentiation. We suggest that let-7a is mostly involves in TGF-β modulation, where mir-26a leads to all steps of osteoblast development and maturation via BMP–SMAD targeting. Additionaly, TGF-β/BMP–SMAD dependent pathway is supposed to be regulated by lowering in mir-17 which targets SMAD7, SMAD6 and SMAD1. But, surprisingly, overexpressed miRs mir-21 and mir-101 also targets SMAD proteins (SMAD7 and SMAD6) inhibiting SMAD dependent transcriptions.

However, non-SMAD-dependent pathway (JNK, ERK, p38 MAPK, Pi3K-Akt signaling) also promote osteogenesis and main transcriptional factors (Runx2, Osx, Dlx5 and other) expression [12-15, 26]. Mentioned signaling can activate the Runx2 expression in hMSCs. Moreover, all these signaling together with mentioned above BMP pathway regulates the Osterix expression Runx2 independently. With DAVID Bioinformatics Resources we found that non-SMAD-dependent pathway activity is orchestrated by interplaying between different up- and downregulated miRs. As we mentioned above let-7a targets TGF-β–SMAD signaling and is crucial for entering of hMSCs into osteogenesis, however this microRNA regulates the TGF-β–non-SMAD-dependent signaling including MAPK, PI3K-Akt, cGMP-PKG, Jak-STAT, Hippo and AMPK positively affected their activity. Our data showed that let-7a is versatile regulator of osteogenesis and high level of its expression maintain osteoprogenitors differentiation as well as osteoblast maturation and metabolism (via AMPK pathway).

Our data show that not only let-7a, but MAPK, ErbB/MAPK and Pi3K-Akt signaling are predictably regulated by mir-29b, mir-101, mir-143 and mir-218 that are also expressed at higher level during osteogenic differentiation. Additionaly MAPK and Pi3K-Akt signaling supposed to be regulated by let-7a and mir-26a via cGMP-PKG signaling pathway and mir-21, mir-101, mir-29b and mir-218 through Ras signaling pathway. However, GO analysis revealed that mir-21 can affect negatively the ERK1/1 and Ras signaling pathway via inhibiting of MAPK activity. As shown, Ras/MAPK signaling branch activity is associated with osteogenesis inhibiting [32], so mir-21 may inactivate Ras/MAPK signaling promoting in such way osteogenic differentiation. We found that mir-101 and mir-218 also are possible inactivators of MAPK-cascade and Jak-STAT signaling. Thus, our research indicates that high content of let-7a, mir-21, mir-26a, mir-29b, mir-143 can be specifically involved in non-SMAD dependent signaling and promote osteogenic TF expression (Fig. 5). On the other side non-SMAD dependent signaling modulated by miRs shows significantly lower expression in osteogenesis. Namely, mir-17, mir-222, mir-138 and mir-31 are supposed to be involved in non-SMAD-dependent pathway (i.e., MAPK, Pi3K-Akt, and canonical Wnt signaling) regulation. However, GO analysis indicates that mir-17 positively regulates Wnt and p38-MAPK signaling cascades where inhibits TGF-β signaling. Hence, we suppose that mir-17 level diminishing is necessary for TGF-β/non-SMAD dependent signaling activity. Our data show that mir-31 also positively regulates MAPK-cascade, however low level of its content points that this microRNA does not require for MAPK-dependent expression of osteogenic TF.

GO analysis reviled that mir-222 and mir-138 can modulate canonical Wnt via noncanonical calcium and Hippo signaling, respectively. Moreover mir-138, based on KEGG pathway analysis, is supposed to negatively regulate canonical Wnt-dependent transcription through Notch signaling pathway. Hence, the lower level of mir-222 and mir-138 in hMSCs maintains the canonical Wnt signaling up-regulation, where the latter one regulates other osteogenic TF expression (Runx2, Osx, Dlx5 and ATF4), strongly stimulates osteoblastogenesis and inhibits adipogenic TF (C/ERB and PPARγ) [10, 13]. Moreover, canonical Wnt involves in osteoblast maturation and metabolism [33]. GO analysis revealed that canonical Wnt activity can be positively regulated by mir-26a, mir-101 and mir-29b which expresses on higher level during osteogenesis, where mir-101 and mir-29b regulate β-catenin binding and β-catenin-TCF complex assembly. Summarizing these data we came to the conclusion that canonical Wnt signaling activity in osteogenesis depends on low content of mir-222, mir-17 and mir-138 and high expression of mir-26a, mir-101 and mir-29b.

Other crucial for osteogenesis signaling is Hippo/Yap driven by a few miRs. Notably, Yap is a co-regulator for main osteogenic TF and act via both SMAD-dependent and non-SMAD-dependent molecular mechanisms, so it is up-regulation is crucial for osteogenic lineage commitment as well as osteoblasts maturation [34,35]. GO and KEGG pathway analysis revealed that Yap dependent transcription can be modulated by up-regulated let-7a, mir-101 and downregulated mir-31, mir-138 and mir-222. Our finding points that mir-138 and mir-222 can regulate Yap/TCF/Lef dependent transcription. Mir-31 and mir-222 are supposed to inhibit Yap/Taz and prevent its transcription activity, where mir-17 modulate TGF-β–SMAD/Yap, BMP–SMAD/Yap and Yap/TCF/Lef dependent signaling. Summarizing we assume that lower level of mir-31, miR-138, mir-222 are important for Yap dependent transcription of main osteogenic TF in hMSC during all steps of osteogenesis (Fig. 5).

Osteoblast differentiation relates to metabolic changes, hence AMPK signaling pathway is one of crucial factors which controls different aspects of metabolism in cells including osteoprogenitors and osteoblast. Moreover, AMPK signaling is required for hMSC to enter the adipogenic and osteogenic lineages through Runx2 and PPARγ and to be able to maintain osteogenesis and osteoblast homeostasis [36]. With bioinformatic analysis we found that AMPK signaling during osteogenesis can be positively modulated by let-7a, mir-29b and mir-218 which are significantly up-regulated (Fig. 5).

Other important signaling involved in osteoblast homeostasis as well as differentiation is Hedgehog signaling pathway, which play dual role in bone homeostasis and required in both bone resorption and formation process. Additionally, Hedgehog signaling induce osteoblast differentiation in BM-MSC through Runx2 and Osx expression [37]. We found that during osteogenesis mir-101 can promote Hedgehog signaling activity.

Osteogenic differentiation, osteoblast maturation and homeostasis are coordinated by several signaling pathways which has many have many crosstalks and common targets genes. In this study, we used bioinformatic approaches to identify specific association between cells signaling and microRNAs up- and downregulated during osteogenic lineage differentiation. We suppose mentioned below miRNA to be specifically involved in hMSC transition into osteogenic differentiation via microenvironment formation (i.e. let-7a, mir-17, mir-21, mir-29b and mir-101) and TGF-β/BMP–SMAD dependent pathway (i.e. let-7a, mir-17, mir-21, mir-26a mir-101) as well as MAPK signaling pathway (i.e. let-7a, mir-21, mir-26a, mir-29b, mir-143 and mir17), while Yap-dependent expression of osteogenic TF is modulated by let-7a, mir-31mir-101, mir-138 and mir-222. We identify mir-17, mir-26a, mir-29b, mir-101, mir-138 and mir-222 as those that can be specifically involved in canonical Wnt signaling dependent osteogenesis as well as in osteoblast maturation and metabolism maintainance together with let-7a, mir-29b and mir-218 which modulate AMPK signaling. Additionally, identified mir-101 is likely involved in osteoblast homeostasis via Hedgehog signaling.

Our data expands knowledge in the field of hMSCs fate and osteogenesis orchestration by miRs, demostrates pro-osteogenic and anti-osteogenic miRs and their potential molecular pathways.

Search strategy and selection criteria

We carried out a comprehensive literature search in the PubMed, Medline and Google Scholar databases using the following keywords: miR, microRNA, osteogenesis, osteogenic differentiation, MSC. The articles were not limited by the publication date, but their language was restricted to English. A total of 632 papers were identified with the initial search. The inclusion criterion for article selection was the presence of original research data in the field of miRNA expression in osteogenesis or osteogenic differentiation. The exclusion criteria were: 1) the paper did not contain original research data; 2) the paper did not contain numeric data or any formalized miRNA expression data; 3) data provided in the paper could not be compared with other research data.

As a result, 87 papers were included in this meta-analysis. (Fig. 1). The detailed information of relevant citations is listed in Supplement 1.

Data extraction and study assessment

Two researchers independently extracted data and reviewed the essence of the articles to determine whether they met the criteria for inclusion. All disagreements were discussed.

Gathered data were analyzed by the following criteria: 1) miRNA expression has been measured; 2) process has been studied. In total 2480 different miRNA have been analyzed in included papers, but only miRNA mentioned by 3 different authors have been studied.

Target gene analysis

To identify potential gene targets for each miRNA we used specialized databases such as Targetscan (http://www.targetscan.org/vert_72/), miRDB (http://mirdb.org/), and miRanda (http://www.microrna.org/microrna/home.do). For common gene targets visualization, the InteractiVenn (http://www.interactivenn.net/) instrument was employed [24]. Then, for identified gene analysis, DAVID Bioinformatics Resources 6.8 (https://david.ncifcrf.gov/) was used, namely we performed Gene Ontology (GO) function verification and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis.

Statistical analysis

The Aggregate Data approach was used. First, we calculated the Pearson correlation coefficient using the standard procedures. Then we performed one-way ANOVA for the study of changes in levels of individual miRNAs during processes of osteogenesis, osteogenic differentiation, and in osteopathological conditions. For statistical analysis of the received data, the simple regression method was used.

We did not employ data weighting methodology, as data used for the purposes of this study did not include any clinical trial data and all experimental groups in the studies were rather small.

The effect size was calculated by Hedges and Olkin.

Both types of publications – those that demonstrated a significant changes of miRNA expression and/or experimental modulation of processes have been studied, and those that did not - were included into analysis; the fail-safe number was calculated by Rosenthal method.

All calculations were performed with Origin 8.1.

Author Contributions:

LLM, and OOP planned work, discussed the results, and wrote the manuscript. DSA, LLM extracted and analyzed data. LLM performed the statistical analysis. Bioinformatic analysis was done by OOP. All authors participated in the design of the study, interpretation, and analysis of the data, reviewed of the manuscript.

Funding:

This work was funded by the Ukrainian State Budget Program
"Support for the Development of Priority Areas of Scientific Research" (Code: 6541230).

Availability of data and materials

The datasets used and analyzed during the current study are available from the corresponding author on request.

Conflicts of Interest:

The author declares no conflict of interest.

Zhang Y, Gordon A, Qian W, Chen W. Engineering Nanoscale Stem Cell Niche: Direct Stem Cell Behavior at Cell-Matrix Interface. Advanced Healthcare Materials. 2015; 4(13): 1900–1914. doi:10.1002/adhm.201500351
Discher D, Mooney D, Zandstra P. Growth factors, matrices, and forces combine and control stem cells Science. 2009; 324(5935): 1673–1677. doi:10.1126/science.1171643.
Landrier J-F, Derghal A, Mounien L. MicroRNAs in Obesity and Related Metabolic Disorders. Cells. 2019; 8(8):859
Sera S, zur Nieden S. MicroRNA Regulation of Skeletal Development. Curr Osteoporos Rep. 2017: 15(4): 353–366. doi:10.1007/s11914-017-0379-7.
Yeh YT, Wei J, Thorossian S, Nguyen K, Hoffman C, del Álamo CJ, et al. MiR-145 mediates cell morphology-regulated mesenchymal stem cell differentiation to smooth muscle cells. Biomaterials. 2019; 204:59-69. https://doi.org/10.1016/j.biomaterials.2019.03.003
Gebauer F, Hentze MW. Molecular mechanisms of translational control. Nat Rev Mol Cell Biol. 2004; 5(10):827–835
Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004; 116(2):281–297
Jing D, Hao J, Shen Y, et al. The role of microRNAs in bone remodeling. Int J Oral Sci. 2015; 7: 131–143. https://doi.org/10.1038/ijos.2015.22
Bruderer M, Richards RG, Alini M, Stoddart MJ. Role and regulation of RUNX2 in osteogenesis. European Cells and Materials. 2014; 28: 269-286. DOI: 10.22203/eCM.v028a19
Felber K, Elks PM, Lecca M, Roehl HH. Expression of osterix Is Regulated by FGF and Wnt/β-Catenin Signalling during Osteoblast Differentiation. PLoS ONE. 2015; 10(12): e0144982. https://doi.org/10.1371/journal.pone.0144982
Zhang. Transcriptional regulation of bone formation by the osteoblast-specific transcription factor Osx. Journal of Orthopaedic Surgery and Research. 2010; 5:37. doi: 10.1186/1749-799X-5-37
Wu M, Chen G, Li YP. TGF-β and BMP signaling in osteoblast, skeletal development, and bone formation, homeostasis and disease. Bone Res. 2016; 4:16009. https://doi.org/10.1038/boneres.2016.9
Houschyar KS, Tapking C, Borrelli MR, Popp D, Duscher D, Maan ZN, et al. Wnt Pathway in Bone Repair and Regeneration – What Do We Know So Far. Frontiers in Cell and Developmental Biology. 2019; 6:170. DOI=10.3389/fcell.2018.00170
Rodríguez-Carballo E, Gámez B, Ventura F. p38 MAPK Signaling in Osteoblast Differentiation. Frontiers in Cell and Developmental Biology. 2016; 4:40. DOI=10.3389/fcell.2016.00040
Sanpaolo ER, Rotondo C, Cici D, et al. JAK/STAT pathway and molecular mechanism in bone remodeling. Mol Biol Rep. 2020;47:9087–9096. https://doi.org/10.1007/s11033-020-05910-9
Mohri T, Nakajima M, Takagi S, Komagata S, Yokoi T. MicroRNA regulates human vitamin D receptor. Int J Cancer. 2009; 125 (6): 1328-1333. 10.1002/ijc.24459.
Mizuno Y, Yagi K, Tokuzawa Y, Kanesaki-Yatsuka Y, Suda T, Katagiri T, et al. miR-125b inhibits osteoblastic differentiation by down-regulation of cell proliferation. Biochem Biophys Res Commun. 2008; 368 (2): 267-272. 10.1016/j.bbrc.2008.01.073.
Peng H, Lu S-L, Bai Y, Fang X, Huang H, Zhuang H-Q. MiR-133a inhibits fracture healing via targeting RUNX2/BMP2. Eur Rev Med Pharmacol Sci. 2018; 22(9):2519-2526.
DOI: 10.26355/eurrev_201805_14914
Guo P-Y, Wu L-F, Xiao Z-Y, Huang T-L, Li X. Knockdown of MiR-140-5 promotes osteogenesis of adipose-derived mesenchymal stem cells by targeting TLR4 and BMP2 and promoting fracture healing in the atrophic nonunion rat model. Eur Rev Med Pharmacol Sci. 2019;23(5):2112-2124. DOI: 10.26355/eurrev_201903_17255
Wu X, Gu Q, Chen X, Mi W, Wu T, Huang H. MiR-27a targets DKK2 and SFRP1 to promote reosseointegration in the regenerative treatment of peri-implantitis. J. Bone Miner Res.2018;34: 123–157.
Wang C, Qiao X, Zhang Z, Li C. MiR-128 promotes osteogenic differentiation of bone marrow mesenchymal stem cells in rat by targeting DKK2. Biosci Rep. 2020;40(2):BSR20182121. doi: 10.1042/BSR20182121.
Long H, Zhu Y, Lin Z, et al. miR-381 modulates human bone mesenchymal stromal cells (BMSCs) osteogenesis via suppressing Wnt signaling pathway during atrophic nonunion development. Cell Death Dis. 2019; 10:470. https://doi.org/10.1038/s41419-019-1693-z
Xu1 G, Ding Z, Shi H-F. The mechanism of miR-889 regulates osteogenesis in human bone marrow mesenchymal stem cells. Journal of Orthopaedic Surgery and Research. 2019; 14:366 https://doi.org/10.1186/s13018-019-1399-z
Heberle H, Meirelles G V, da Silva F R, Telles GP, Minghim R. InteractiVenn: a web-based tool for the analysis of sets through Venn diagrams. BMC Bioinformatics. 2015; 16:169. DOI: 10.1186/s12859-015-0611-3
Jeggari A, Alekseenko Z, Dias J, Ericson J, Alexeyenko A. EviNet: a web platform for network enrichment analysis with flexible definition of gene sets. Nucleic Acids Res. 2018;46(1):W163-W170. doi: 10.1002/1878-0261.12350.
Huang W, Yang S, Shao J, Li Y-P. Signaling and transcriptional regulation in osteoblast commitment and differentiation. Front Biosci. 2007;12: 3068–3092.
Peng S, Gao D, Gao C, Wei P, Niu M, Shuai C. MicroRNAs regulate signaling pathways in osteogenic differentiation of mesenchymal stem cells (Review). Mol Med Rep. 2016;14: 623-629.
Fröhlich LF. Micrornas at the Interface between Osteogenesis and Angiogenesis as Targets for Bone Regeneration. Cells. 2019; 8(2): doi:10.3390/cells8020121
van Wijnen A J, van de Peppel J, van Leeuwen J P, Lian J B, Stein G S, Westendorf J J, et al. MicroRNA functions in osteogenesis and dysfunctions in osteoporosis. Current osteoporosis reports. 2013; 11(2):72–82. https://doi.org/10.1007/s11914-013-0143-6
Grimson A, Farh K K, Johnston W K, Garrett-Engele P, Lim L P, Bartel D P. MicroRNA targeting specificity in mammals: Determinants beyond seed pairing. Mol. Cell. 2007;27:91–105. doi: 10.1016/j.molcel.2007.06.017.
Betel D, Wilson M, Gabow A, Marks D S, Sander C. The microRNA.org resource: Targets and expression. Nucleic Acids Res. 2008;36:D149–D153. doi: 10.1093/nar/gkm995.
Schindeler A, Little D G. Ras‐MAPK Signaling in Osteogenic Differentiation: Friend or Foe? JBMR Volume. 2006; 21(9):1331-1338.
Karner C M, Long F. Wnt signaling and cellular metabolism in osteoblasts. Cellular and molecular life sciences: CMLS. 2017;74(9):1649–1657. https://doi.org/10.1007/s00018-016-2425-5
Pan JX, Xiong L, Zhao K, et al. YAP promotes osteogenesis and suppresses adipogenic differentiation by regulating β-catenin signaling. Bone Res. 2018;6:18.
Brandaõ A S, Bensimon-Brito A, Lourenço R, Borbinha1 J, Soares A R, Mateus R, Jacinto A. Yap induces osteoblast differentiation by modulating Bmp signalling during zebrafish caudal fin regeneration. Journal of Cell Science. 2019; 132: jcs231993. doi:10.1242/jcs.231993
Jeyabalan J, Shah M, Viollet B, Chenu C. AMP-activated protein kinase pathway and bone metabolism. Journal of Endocrinology. 2012;212(3):277-290.
Lv W T, Du D H, Gao R J, Yu C W, Jia Y, Jia Z F, Wang C J. Regulation of Hedgehog signaling Offers A Novel Perspective for Bone Homeostasis Disorder Treatment. International journal of molecular sciences. 2019;20(16):3981. https://doi.org/10.3390/ijms20163981
Wang J, Liu S, Li J, et al. Roles for miRNAs in osteogenic differentiation of bone marrow mesenchymal stem cells. Stem Cell Res Ther. 2019; 10: 197. https://doi.org/10.1186/s13287-019-1309-7
Trompeter et al.: MicroRNAs miR-26a, miR-26b, and miR-29b accelerate osteogenic differentiation of unrestricted somatic stem cells from human cord blood. BMC Genomics. 2013; 14:111.
Yang XM, Song YQ, Li L, et al. miR-1249-5p regulates the osteogenic differentiation of ADSCs by targeting PDX1. J Orthop Surg Res. 2012;16:10. https://doi.org/10.1186/s13018-020-02147-x
Jiang X, Zhang Z, Peng T, Wang G, Xu Q, Li G. miR204 inhibits the osteogenic differentiation of mesenchymal stem cells by targeting bone morphogenetic protein 2. Mol Med Rep. 2020; 21: 43-50.

No competing interests reported.

supplementary.xls

Bioinformatic Analysis of MicroRNAs and Signaling Networks that Orchestrates the hMSC Osteogenic Differentiation

Status:

Version 1

Abstract

Figures

Introduction

Results

Discussion

Materials And Methods