Computational and Mathematical Methods in Medicine
Volume 2017 (2017), Article ID 9426280, 9 pages
https://doi.org/10.1155/2017/9426280

Research Article

Construction of an miRNA-Regulated Pathway Network Reveals Candidate Biomarkers for Postmenopausal Osteoporosis

Min Shao

Department of Orthopedics, Third Affiliated Hospital of Guangzhou University of Chinese Medicine, Guangzhou, Guangdong Province 510240, China

Correspondence should be addressed to Min Shao; shaomin6980@163.com

Received 14 June 2017; Revised 15 August 2017; Accepted 21 August 2017; Published 11 October 2017

Academic Editor: Nadia A. Chuzhanova

Copyright © 2017 Min Shao. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We aimed to identify risk pathways for postmenopausal osteoporosis (PMOP) via establishing an microRNAs- (miRNA-) regulated pathway network (MRPN). Firstly, we identified differential pathways through calculating gene- and pathway-level statistics based on the accumulated normal samples using the individual pathway aberrance score (iPAS). Significant pathways based on differentially expressed genes (DEGs) using DAVID were extracted, followed by identifying the common pathways between iPAS and DAVID methods. Next, miRNAs prediction was implemented via calculating TargetScore values with precomputed input (log fold change (FC), TargetScan context score (TSCS), and probabilities of conserved targeting (PCT)). An MRPN construction was constructed using the common genes in the common pathways and the predicted miRNAs. Using false discovery rate (FDR) < 0.05, 279 differential pathways were identified. Using the criteria of FDR < 0.05 and , 39 DEGs were retrieved, and these DEGs were enriched in 64 significant pathways identified by DAVID. Overall, 27 pathways were the common ones between two methods. Importantly, MAPK signaling pathway and PI3K-Akt signaling pathway were the first and second significantly enriched ones, respectively. These 27 common pathways separated PMOP from controls with the accuracy of 0.912. MAPK signaling pathway and PI3K/Akt signaling pathway might play crucial roles in PMOP.

1. Introduction

Estrogen plays a crucial role in keeping healthy bones. Estrogen deficiency has significantly negative effect on bone cell functions [1]. Postmenopausal osteoporosis (PMOP) is believed to result directly from the decreased endogenous estrogen in menopausal women [2, 3]. PMOP is characterized by the severe loss of bone mass in the vertebrae and long bones. Globally, thousands of women over 50 years of age are influenced by PMOP which causes an increase of economic and societal burden recently. Therefore, studying PMOP-related molecular mechanisms is urgently needed.

In recent years, previous studies have demonstrated that there are powerful genetic effects on the etiology of osteoporosis (OP) in postmenopausal women [4–6]. For example, polymorphisms in a set of genes have been found to be related to PMOP women in China, including ESR2 [7], HLA-B [8], and OPG [9]. Nevertheless, the etiology of PMOP still remains limited. MicroRNAs (miRNAs), as a kind of small noncoding RNAs of 22–25 nucleotides, have been implicated to exert indispensable regulatory functions in gene expression in animal kingdoms [10]. miRNAs have also been reported to play important roles in the regulation of biological processes, such as cell apoptosis, proliferation, and differentiation [11, 12]. Growing evidence has demonstrated that miRNAs have an important role in regulating bone mass. In recent years, the role of miRNAs in the progress of OP has gained interest. Numerous studies have demonstrated that the dysregulation of miRNAs is associated with OP. miRNA-2861 has been found to affect osteoblast differentiation, contributing to OP in the femur via its effects on osteoblasts [13]. Moreover, miRNA-34c has been implicated to enhance the osteoclastogenesis through regulating the genes involved in the Notch signaling pathway [14]. Moreover, Chen et al. [15] have demonstrated that miRNA-125b is significantly overexpressed in PMOP, and miRNA-125b may be a potentially noninvasive biosignature for PMOP. Similarly, Xiao et al. [16] have suggested the role of miRNA-129-5p in osteoblast differentiation and bone homeostasis. These findings suggest that miRNAs have potential effects on the pathways in the progression of OP and can thereby affect the treatment of OP. However, most of the studies have focused on only one or a few miRNAs/pathways using limited samples. Of note, little is known about specific pathways regulated by miRNAs during the development of PMOP.

Hence, we aimed to detect the risk pathways regulated by PMOP-related miRNAs through establishing an miRNA-regulated pathway network (MRPN). Then, we analyzed the topological characteristics of the MRPN, and we also determined how the miRNAs regulated PMOP risk pathways. The results of our study may provide a novel viewpoint with respect to the mechanism and treatment of PMOP.

2. Methods

The current analysis was comprised of the following steps-data selection (microarray profile, pathway data): pathway identification through calculating gene-level statistics and pathway-level statistics based on the accumulated normal samples using individual pathway aberrance score (iPAS), finding differentially expressed genes (DEGs) and significant pathways based on DEGs using DAVID software, identification of the common pathways of these two methods, miRNAs prediction based on the common genes via calculating TargetScores, MRPN construction, and the topological properties for MRPN. The detailed information of each step was described in Figure 1.

Figure 1: The flowchart showing the workflow step-by-step.

2.1. Data Availability

To reveal the molecular mechanisms of PMOP, microarray analyses of monocytes were performed using Affymetrix HG-133A arrays in 80 females, including 40 pre- and 40 postmenopausal subjects. These data were deposited in the E-GEOD-56815 of the ArrayExpress database [14]. Probe IDs owning concentrated expression level were transformed into human gene symbols. Duplicated genes in matrix were eliminated. Overall, we obtained 8450 genes.

2.2. Collection of Pathway Data

All biological pathways of humans were retrieved from REACTOME database (http://www.reactome.org/) [17]. Then, several pathways were filtered out when the gene set size of these pathways was larger than 100. In addition, the pathways with the intersection of 0 between genes enriched in the pathways and 8450 genes identified in the present analysis were also discarded. Ultimately, 1078 pathways were reserved for further analysis.

2.3. Gene-Level Statistics

Herein, accumulated premenopausal samples were used as a reference and the expression level of genes was computed by comparing one postmenopause sample with many accumulated premenopause samples according to robust multichip average [18]. In detail, we firstly normalized the genes in all the premenopause subjects for the reference one by one to the reference, and then we calculated the mean value as well as standard deviation (SD) of the expression level. For a single postmenopause sample, quantile normalization was conducted after integrating the individual postmenopause sample with all reference samples. Assuming the genes having several probes, gene expression value was computed via averaging probe expression value. Gene-level statistics of each gene in an individual postmenopause subject was standardized according to the mean value and SD of the reference genes.

2.4. Pathway-Level Statistics

We used the average method to calculate the pathway activity. Specifically, the gene-level statistics of all genes in each pathway were measured and summarized, following by transforming the mean value of gene-level into the pathway statistics of this pathway. The pathway statistics was calculated using the following formula:

In this formula, iPAS is the expression status of a pathway, represents the gene-level statistics of th gene, and denotes the count of genes enriched in this given pathway.

2.5. Identification of Differential Pathways

With the goal of analyzing the changes of pathways in PMOP and control samples, Gitools [19] was utilized to construct the cluster heatmap of pathways. -test was used to measure the pathway statistics of each pathway in all samples of the two groups. All the premenopause samples for the reference were compared one by one with reference to output the null distribution of pathway statistics. On the basis of the comparison between an individual postmenopause sample and the null distribution, we obtained the raw value. Subsequently, the values were corrected for multiple hypothesis testing based on false discovery rate (FDR) control using Benjamini and Hochberg (BH). Pathways with a FDR value of less than 0.05 were considered significantly differential.

2.6. Identification of DEGs

In our study, LIMMA package [20] of R language was used to calculate the expression level of each gene between two groups to further identify DEGs. Next, all expression scores were converted into fold-changes (FC) with base (herein, this means ). The for each gene was determined as log(postmenopausal) − log(premenopausal), and the distribution of the value of each gene expression was investigated. was on behalf of the differential expression degree. Then, we used multiple test to adjust the original values by means of Benjamini et al. [21] method using FDR. DEGs were finally extracted when the threshold was set at FDR < 0.05 and .

2.7. Pathway Analyses for DEGs

To further reveal the biological roles of DEGs, DAVID [22] was used to perform the pathway analysis using the Expression Analysis Systematic Explorer (EASE) test [23]. The EASE score was utilized to extract the significant categories. In our work, the pathways having < 0.05 were regarded to be significant pathways.

2.8. Comparison of iPAS and DAVID

To evaluate whether the iPAS method was feasible, the method proposed in the current study was compared with traditional DAVID software. The common pathways between the two methods were identified. In an attempt to evaluate the classification capacity of the common pathways, the hierarchical clustering analysis was implemented for the common pathways using Gitools [19].

2.9. TargetScore Calculation and Prediction of Potential miRNAs

In detail, TargetScan context score (TSCS) is a sequence-based score for the individual target site, which was computed using TargetScan [24]. Probabilities of conserved targeting (PCT) are the probability of conserved targeting for single target site [25]. TargetScore is a flexible package which takes values and sequence scores as inputs to calculate the probability of each gene being the target of a given miRNA. The function “TargetScore” offers a convenient way to acquire TargetScore with precomputed input (, TSCS, and PCT), and the range of TargetScore is from 0 to 1 [26]. The greater the TargetScore is, the higher the accuracy in detecting known targets is. Therefore, relying on the , TSCS, and PCT values, we counted the TargetScore for the genes enriched in the significant pathways to further identify the potential miRNAs by means of Variational Bayesian-Gaussian Mixture Model (VB-GMM) [27]. The TargetScore distribution for validated and nonvalidated targets of all miRNA-mRNA interactions each having at least 1 validated targets was evaluated. In the current study, the predefined was determined as the threshold when little overlap between the two distributions when TargetScore = κ existed. Furthermore, the miRNA-regulated pathways were achieved using their target gene sets. Next, we constructed an MRPN consisting of miRNAs, genes, and their common pathways. Cytoscape 2.8.3 was used to visualize the MRPN, and the topological properties of the MRPN were analyzed based on the Network Analysis plugin [28].

3. Results

3.1. Identification of Differential Pathways

After the gene-level statistics of all genes was transformed into the pathway-level statistics value of each pathway, -test was used for measuring the pathway statistics of each pathway in control and PMOP samples. Moreover, FDR approach was applied to correct significance levels ( values) for multiple hypothesis testing. Using the criteria of FDR < 0.05, a total of 279 pathways were extracted between the two groups. Cluster analysis was conducted to explore the changes of these 279 pathways. The heatmap of the 279 pathways was shown in Figure 2.

Figure 2: Cluster analysis of using average approach as the iPAS on the dataset of postmenopausal osteoporosis (PMOP) through Gitools method. Pathways ( = 279) and samples are clustered based on iPAS. The color scale represents the relative levels of pathway aberrance. Horizontal axis is samples; vertical coordinate is differential pathways. iPAS means individual pathway aberrance score.

3.2. Identifying DEGs

Before pathway analysis, DEGs between the two groups were firstly identified. Using the criteria of FDR < 0.05 and , a total of 39 DEGs were retrieved. Specific information about DEGs was exhibited in Table 1. The most 5 significant DEGs were PTGES2, HAB1, SCT, ASNSD1, and TUT1.

Table 1: List of differentially expressed genes (DEGs).

Genes
Log (fold change)
False discovery rate (FDR)

PTGES2

HAB1

SCT

ASNSD1

TUT1

SEPT5-GP1BB

NXT2

ACIN1

RPRD1A

GMFB

GNL1

MORN1

SH3GLB2

NEK4

NUPR1

AES

DEF6

LOC100507630

ABO

KCNG2

ITGB7

TTC33

KCNF1

F7

KIAA1551

DSPP

RNASE4

MAP1S

DYNLT3

YIPF5

ITM2C

HSPB1

ZFYVE16

UAP1L1

CBR4

TSC22D1

HAMP

ZFP36L2

3.3. Comparison of iPAS and DAVID

To evaluate whether the iPAS method was feasible, the method proposed in the current study was compared with traditional DAVID software. It was observed that a total of 27 pathways were the common ones obtained from the iPAS approach and DAVID software. These common pathways were shown in Table 2. More importantly, hsa0410 (MAPK signaling pathway) was the first significantly enriched pathway. Of note, the pathway quantity selected from iPAS was increased, relative to the DAVID method (a total of 279 differential pathways were obtained by means of iPAS method, while only 64 significant pathways were extracted based on the traditional DAVID). Thus, our study demonstrated that the iPAS method showed more efficiency to identify significant pathways, compared with DAVID.

Table 2: List of the common pathways between the two methods.

Pathways
DEGs

MAPK signaling pathway [PATH: hsa04010]
HSPB1

PI3K-Akt signaling pathway [PATH: hsa04151]
ITGB7

VEGF signaling pathway [PATH: hsa04370]
HSPB1

TGF-beta signaling pathway [PATH: hsa04350]
ZFYVE16

Transcriptional misregulation in cancers [PATH: hsa05202]
NUPR1, ITGB7, CEBPE

mRNA surveillance pathway [PATH: hsa03015]
NXT2, ACIN1

RNA transport [PATH: hsa03013]
NXT2, ACIN1

Glycosphingolipid biosynthesis-lacto and neolacto series [PATH: hsa00601]
ABO

Endocytosis [PATH: hsa04144]
SH3GLB2, ZFYVE16

Amino sugar and nucleotide sugar metabolism [PATH: hsa00520]
UAP1L1

Intestinal immune network for IgA production [PATH: hsa04672]
ITGB7

Arachidonic acid metabolism [PATH: hsa00590]
PTGES2

Complement and coagulation cascades [PATH: hsa04610]
F7

Arrhythmogenic right ventricular cardiomyopathy (ARVC) [PATH: hsa05412]
ITGB7

TGF-beta signaling pathway [PATH: hsa04350]
ZFYVE16

Ribosome biogenesis in eukaryotes [PATH: hsa03008]
NXT2

Hypertrophic cardiomyopathy (HCM) [PATH: hsa05410]
ITGB7

ECM-receptor interaction [PATH: hsa04512]
ITGB7

Dilated cardiomyopathy (DCM) [PATH: hsa05414]
ITGB7

Pancreatic secretion [PATH: hsa04972]
SCT

Amoebiasis [PATH: hsa05146]
HSPB1

Spliceosome [PATH: hsa03040]
ACIN1

Cell adhesion molecules (CAMs) [PATH: hsa04514]
ITGB7

Influenza A [PATH: hsa05164]
NXT2

Epstein-Barr virus infection [PATH: hsa05169]
HSPB1

Focal adhesion [PATH: hsa04510]
ITGB7

Regulation of actin cytoskeleton [PATH: hsa04810]
ITGB7

Neuroactive ligand-receptor interaction [PATH: hsa04080]
SCT

The clustering analysis was carried out to evaluate the classification capacity for microarray data samples using the common 27 pathways detected by iPAS and DAVID method (Figure 3). The classification efficiency for microarray data samples was computed based on accuracy. These pathways separated PMOP from controls with the accuracy of 0.912. In light of this, we further infer that these pathways could be used to be marker pathways for PMOP with a high accuracy.

Figure 3: The cluster heatmap of the common pathways based on Gitools method. The color scale stood for the iPAS level; horizontal axis denoted samples; vertical coordinate represented pathways.

3.4. Construction of MRPN for the Common Pathways Based on the Common Genes of miRNAs and Pathways

TargetScores were firstly calculated based on the combination of values, and TSCS as well as PCT scores. Then, the distribution of TargetScores was obtained. According to the distribution of TargetScores, little overlap between the two distributions at TargetScore = 0.4 was found. Hence, TargetScore > 0.4 was selected as the threshold to filter out less confidence miRNAs. Using the TargetScore value > 0.4, a total of 39 genes, 50 miRNAs, and 145 miRNA-gene interactions were identified, and specific relationship was displayed in Figure 4.

Figure 4: Relationships between the common genes in the common pathways and miRNAs.

We then constructed an MRPN for PMOP, which was shown in Figure 5. In this network, the most significant common pathway (MAPK signaling pathway) was regulated by has-miRNA-9. Furthermore, has-miRNA-34b-5p and has-miRNA-530f indirectly regulated 10 significantly differential pathways; among these 10 pathways, hsa-04151 (PI3K-Akt signaling pathway) was the second most significant common pathway.

Figure 5: MiRNA-regulated pathway network (MRPN) in PMOP. Network organization of pathway, gene, and miRNA associations. Red nodes represented miRNAs, and blue squares as well as blue nodes indicated genes, pathways, respectively. Yellow ones were the hub nodes.

3.5. Topological Features of the MRPN

We determined the topological characteristics of the network, including the degree distribution. The degree distribution followed a power law distribution using all of the nodes in the MRPN. Based on the degree distribution, we selected 5 nodes as hub nodes, including NXT2, has-miR-1285, GMFB, has-miR-530f, and RPRD1A.

4. Discussion

In our study, a novel method of extracting pathway-based biomarkers in PMOP was developed, by investigating the interactions among genes, miRNAs, and pathways associated with pathogenesis through constructing an MRPN. The results of our work identified 279 differential pathways between PMOP and normal samples, including MAPK signaling pathway and PI3K/Akt signaling pathway.

MAPKs play crucial roles in cellular response to growth factors or cytokines [29]. The musculoskeletal system is frequently observed in several of the Ras/MAPK pathway syndromes, suggesting that activation of the Ras/MAPK pathway impacts cells regulating bone development and homeostasis [30]. OP is the systemic bone remodeling disease caused by imbalance between bone formation and bone resorption. Significantly, MAPK has been implicated to trigger the induction of c-Fos [31], and c-Fos is a key transcription factor of osteoclastogenesis [32]. Moreover, osteoclastogenesis exerts important functions in the development of OP [33]. MAPK signaling pathway has been demonstrated to be differentially activated in MSCs derived from osteoporotic postmenopausal women [34, 35]. Therefore, we infer that MAPK is important for bone formation and may represent a target to prevent bone loss associated with PMOP. Of note, this pathway was regulated by has-miRNA-9. Has-miRNA-9 was a proinflammatory miRNA [36]. A former study has demonstrated that has-miRNA-9 plays important roles in suppressing tumor growth via mediating the downstream p38 MAPK pathway [37]. Accordingly, has-miRNA-9 might play important roles in the PMOP progression, partially through regulating the MAPK signaling pathway.

Previous studies have suggested that PI3K-Akt signaling pathway is implicated in the regulation of cell growth, proliferation, and survival [38]. Osteoblastic bone formation is stimulated through activating the PI3K-Akt signaling pathway [39]. Bone remodeling is an essential physiological process regulating bone mass and strength; an imbalance in bone resorption and formation results in bone loss during aging and in OP. Moreover, PI3K/Akt pathway is very important in osteoclastogenesis and bone mass in mice, which may offer a potential target in OP [40]. Remarkably, osteoclastogenesis exerts important functions in OP progression [33]. Has-miRNA-530f as a hub node in the MDPN regulated this pathway of PI3K/Akt signaling pathway in our study. According to the literatures, no studies have reported the association of has-miRNA-530f and PMOP. However, based on our results, we speculate that PMOP might be related to has-miRNA-530f-regulated PI3K/Akt signaling pathway.

5. Conclusion

In conclusion, we comprehensively collected the PMOP risk pathways and the miRNAs of PMOP to characterize the molecular mechanism and treatment strategy for PMOP. We then constructed an MRPN to better reveal the relationship among genes, pathways, and miRNAs and to develop promising drug candidates for PMOP. Based on the results, differential pathways, including MAPK signaling pathway and PI3K/Akt signaling pathway, were successfully identified. Moreover, these pathways-related miRNAs may be involved in the pathogenic process of PMOP. The findings of our study may be applied clinically for the diagnosis and treatment of PMOP patients.

Conflicts of Interest

The author declares no conflicts of interest.

Acknowledgments

This study was supported by a grant from the National Natural Science Foundation of China (no. 81373654).

References

S. Khosla, M. J. Oursler, and D. G. Monroe, “Estrogen and the skeleton,” Trends in Endocrinology and Metabolism, vol. 23, no. 11, pp. 576–581, 2012.

R. Marcus, “Post-menopausal osteoporosis,” Best Practice and Research: Clinical Obstetrics and Gynaecology, vol. 16, no. 3, pp. 309–327, 2002.

M. J. Seibel, C. R. Dunstan, H. Zhou, C. M. Allan, and D. J. Handelsman, “Sex steroids, not FSH, influence bone mass,” Cell, vol. 127, no. 6, pp. 1080-1081, 2006.

K. Michaëlsson, H. Melhus, H. Ferm, A. Ahlbom, and N. L. Pedersen, “Genetic liability to fractures in the elderly,” Archives of Internal Medicine, vol. 165, no. 16, pp. 1825–1830, 2005.

Y. Guo, S.-S. Dong, X.-F. Chen et al., “Integrating Epigenomic Elements and GWASs Identifies BDNF Gene Affecting Bone Mineral Density and Osteoporotic Fracture Risk,” Scientific Reports, vol. 6, Article ID 30558, 2016.

C. Wang, Z. Zhang, H. Zhang et al., “Susceptibility genes for osteoporotic fracture in postmenopausal Chinese women,” Journal of Bone and Mineral Research, vol. 27, pp. 2582–2591, 2012.

G. Li, Z. Yao, H. Yang, J. Luo, and L. Han, “Association of CA Repeat Polymorphism,” Journal of Genetics Genomics, vol. 34, pp. 868–876, 2007.

S. M. Li, H. J. Yang, D. X. Zhou, and M. Y. Liu, “Associations between polymorphisms of HLA-B gene and postmenopausal osteoporosis in Chinese Han population,” International Journal of Immunogenetics, vol. 41, pp. 324–329.

J. F. Song, Z. Z. Jing, W. Hu, and Y. X. Su, “Association between single nucleotide polymorphisms of the osteoprotegerin gene and postmenopausal osteoporosis in Chinese women,” Genetics and Molecular Research, vol. 12, no. 3, pp. 3279–3285, 2013.

L. He and G. J. Hannon, “MicroRNAs: small RNAs with a big role in gene regulation,” Nature Reviews Genetics, vol. 5, no. 7, pp. 522–531, 2004.

G. Stefani and F. J. Slack, “Small non-coding RNAs in animal development,” Nature Reviews Molecular Cell Biology, vol. 9, no. 3, pp. 219–230, 2008.

W. Filipowicz, S. N. Bhattacharyya, and N. Sonenberg, “Mechanisms of post-transcriptional regulation by microRNAs: are the answers in sight?” Nature Reviews Genetics, vol. 9, no. 2, pp. 102–114, 2008.

H. Li, H. Xie, W. Liu et al., “A novel microRNA targeting HDAC5 regulates osteoblast differentiation in mice and contributes to primary osteoporosis in humans,” The Journal of Clinical Investigation, vol. 119, no. 12, pp. 3666–3677, 2009.

Y. Bae, T. Yang, H.-C. Zeng et al., “miRNA-34c regulates Notch signaling during bone development,” Human Molecular Genetics, vol. 21, no. 13, Article ID dds129, pp. 2991–3000, 2012.

H. Chen, H. Jiang, D. Can, H. Xu, K. Zhang, and S. Guo, “Evaluation of microRNA 125b as a potential biomarker for postmenopausal osteoporosis,” Tropical Journal of Pharmaceutical Research, vol. 16, no. 3, pp. 641–647, 2017.

W. Z. Xiao, X. C. Gu, B. Hu, X. W. Liu, Y. Zi, and M. Li, “Role of microRNA-129-5p in osteoblast differentiation from bone marrow mesenchymal stem cells,” Cell Mol Biol (Noisy-le-grand), vol. 62, pp. 95–99, 2016.

D. Croft, G. OKelly, G. Wu et al., “Reactome: a database of reactions, pathways and biological processes,” Nucleic acids research, 2010.

R. A. Irizarry, B. Hobbs, F. Collin et al., “Exploration, normalization, and summaries of high density oligonucleotide array probe level data,” Biostatistics, vol. 4, no. 2, pp. 249–264, 2003.

P. L. Christian and L. B. Nuria, “Gitools: Analysis and Visualisation of Genomic Data Using Interactive Heat-Maps,” Plos One, vol. 6, Article ID e19541, 2011.

G. K. Smyth, “Limma: linear models for microarray data,” in Bioinformatics and computational biology solutions using R and Bioconductor, pp. 397–420, Springer, 2005.

Y. Benjamini, D. Drai, G. Elmer, N. Kafkafi, and I. Golani, “Controlling the false discovery rate in behavior genetics research,” Behavioural Brain Research, vol. 125, no. 1-2, pp. 279–284, 2001.

D. W. Huang, B. T. Sherman, and R. A. Lempicki, “Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources,” Nature Protocols, vol. 4, no. 1, pp. 44–57, 2009.

G. Ford, Z. Xu, A. Gates, J. Jiang, and B. D. Ford, “Expression Analysis Systematic Explorer (EASE) analysis reveals differential gene expression in permanent and transient focal stroke rat models,” Brain Research, vol. 1071, no. 1, pp. 226–236, 2006.

D. M. Garcia, D. Baek, C. Shin, G. W. Bell, A. Grimson, and D. P. Bartel, “Weak seed-pairing stability and high target-site abundance decrease the proficiency of lsy-6 and other microRNAs,” Nature Structural and Molecular Biology, vol. 18, no. 10, pp. 1139–1146, 2011.

R. C. Friedman, K. K. Farh, C. B. Burge, and D. P. Bartel, “Most mammalian mRNAs are conserved targets of microRNAs,” Genome Research, vol. 19, no. 1, pp. 92–105, 2009.

S. D. Wilson and P. A. Keddy, “Species competitive ability and position along a natural stress/disturbance gradient.,” Ecology, vol. 67, no. 5, pp. 1236–1242, 1986.

A. A. Khan, D. Betel, M. L. Miller, C. Sander, C. S. Leslie, and D. S. Marks, “Transfection of small RNAs globally perturbs gene regulation by endogenous microRNAs,” Nature Biotechnology, vol. 27, pp. 549–555, 2009.

M. E. Smoot, K. Ono, J. Ruscheinski, P. L. Wang, and T. Ideker, “Cytoscape 2.8: new features for data integration and network visualization,” Bioinformatics, vol. 27, no. 3, pp. 431-432, 2011.

Z. Wei and H. T. Liu, “MAPK signal pathways in the regulation of cell proliferation in mammalian cells,” Cell Research, vol. 12, no. 1, pp. 9–18, 2002.

D. Stevenson, E. Schwarz, J. Carey et al., “Bone resorption in syndromes of the Ras/MAPK pathway,” Clinical Genetics, vol. 80, no. 6, pp. 566–573, 2011.

J. Shan, W. Donelan, J. N. Hayner, F. Zhang, E. E. Dudenhausen, and M. S. Kilberg, “MAPK signaling triggers transcriptional induction of cFOS during amino acid limitation of HepG2 cells,” Biochimica et Biophysica Acta, vol. 1853, no. 3, pp. 539–548, 2015.

Y. Omata, T. Yasui, J. Hirose et al., “Genomewide comprehensive analysis reveals critical cooperation between smad and c-Fos in RANKL-induced osteoclastogenesis,” Journal of Bone and Mineral Research, vol. 30, no. 5, pp. 869–877, 2015.

T.-H. Kuo, T.-H. Lin, R.-S. Yang, S.-C. Kuo, W.-M. Fu, and H.-Y. Hung, “Novel Pyrazole Derivatives Effectively Inhibit Osteoclastogenesis, a Potential Target for Treating Osteoporosis,” Journal of Medicinal Chemistry, vol. 58, no. 12, pp. 4954–4963, 2015.

J. P. Rodríguez, S. Ríos, M. Fernández, and J. F. Santibañez, “Differential activation of ERK1,2 MAP kinase signaling pathway in mesenchymal stem cell from control and osteoporotic postmenopausal women,” Journal of Cellular Biochemistry, vol. 92, no. 4, pp. 745–754, 2004.

G. Xia, Z. Yu, Y. Zhao et al., “Sialoglycoproteins isolated from the eggs of Carassius auratus prevents osteoporosis by suppressing the activation of osteoclastogenesis related NF-κB and MAPK pathways,” Journal of Functional Foods, vol. 17, pp. 491–503, 2015.

M. F. Tolba, S. S. Azab, A. E. Khalifa, S. Z. Abdel-Rahman, and A. B. Abdel-Naim, “Caffeic acid phenethyl ester, a promising component of propolis with a plethora of biological activities: a review on its anti-inflammatory, neuroprotective, hepatoprotective, and cardioprotective effects,” IUBMB Life, vol. 65, no. 8, pp. 699–709, 2013.

J. Lu, H. Luo, X. Liu et al., “miR-9 targets CXCR4 and functions as a potential tumor suppressor in nasopharyngeal carcinoma,” Carcinogenesis, vol. 35, no. 3, Article ID bgt354, pp. 554–563, 2014.

M. A. Krasilnikov, “Phosphatidylinositol-3 kinase dependent pathways: the role in control of cell growth, survival, and malignant transformation,” Biochemistry, vol. 65, no. 1, pp. 59–67, 2000.

Y. Zhang, X. Zeng, L. Zhang, and X. Zheng, “Stimulatory effect of puerarin on bone formation through activation of PI3K/Akt pathway in rat calvaria osteoblasts,” Planta Medica, vol. 73, no. 4, pp. 341–347, 2007.

H. Kang, W. Chang, M. Hurley, A. Vignery, and D. Wu, “Important roles of PI3Kγ in osteoclastogenesis and bone homeostasis,” Proceedings of the National Academy of Sciences of the United States of America, vol. 107, no. 29, pp. 12901–12906, 2010.