Prediction of MicroRNA and Gene Target in Synovium-Associated Pain of Knee Osteoarthritis Based on Canonical Correlation Analysis

Wang, Haiming; Hu, Yue; Xie, Yujie; Wang, Li; Wang, Jianxiong; Lei, Lei; Huang, Maomao; Zhang, Chi

doi:https://doi.org/10.1155/2019/4506876

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2019 | Article ID 4506876 | https://doi.org/10.1155/2019/4506876

Prediction of MicroRNA and Gene Target in Synovium-Associated Pain of Knee Osteoarthritis Based on Canonical Correlation Analysis

Haiming Wang,¹Yue Hu,^2,3Yujie Xie,^2,3Li Wang,^2,3Jianxiong Wang,^2,3Lei Lei,^2,3Maomao Huang,^2,3and Chi Zhang^2,3

Academic Editor: Przemko Tylzanowski

Received14 May 2019

Revised30 Jul 2019

Accepted26 Aug 2019

Published13 Oct 2019

Abstract

Inflammation plays a central role in knee osteoarthritis (OA) pathogenesis (C. R. Scanzello, 2017). The synovial membrane inflammation is associated with disease progression and represents a primary source of agony in knee OA (L. A. Stoppiello et al., 2014). Many inflammatory mediators may have biomarker utility. To identify synovium related to knee OA pain biomarkers, we used canonical correlation analysis to analyze the miRNA-mRNA dual expression profiling data and extracted the miRNAs and mRNAs. After identifying miRNAs and mRNAs, we built an interaction network by integrating miRWalk2.0. Then, we extended the network by increasing miRNA-mRNA pairs and identified five miRNAs and four genes (TGFBR2, DST, TBXAS1, and FHLI) through the Spearman rank correlation test. For miRNAs involved in the network, we further performed the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses, whereafter only those mRNAs overlapped with the Online Mendelian Inheritance in Man (OMIM) genetic database were analyzed. Receiver operating characteristic (ROC) curve and support vector machine (SVM) classification were taken into the analysis. The results demonstrated that all the recognized miRNAs and their gene targets in the network might be potential biomarkers for synovial-associated pain in knee OA. This study predicts the underlying risk biomarkers of synovium pain in knee OA.

1. Introduction

Pain is the main presenting symptom of knee osteoarthritis (OA) patients [3], but the mechanism is rarely known. Synovial inflammation is a major cause of pain in patients with knee OA [4]. Increasing evidence suggested that synovial inflammation plays a significant role in knee OA [5]. Synovium tissue in knee OA is regarded as a treatment target in symptomatic knee OA [6]. Synovial inflammation is promoted by some important mediators [7], and a study has identified some candidate genes associated with synovium pain in knee OA patients [8]. It is crucial to explore the pathogenesis of knee OA from theses biomarker.

Recent studies have established that microRNAs (miRNAs) are mediators of knee OA, which were responsible for the development and progression of OA. For example, circulating miRNA (23a-3p, 24-3p, 27a-3p, 27b-3p, 29c-3p, 34a-5p, and 186-5p) in synovial fluid (SF) has been identified as a potential predictor for knee OA [9]. Bijlsma and Knahr found that an intra-articular injection of miRNA-140 could alleviate knee OA progression by regulating extracellular matrix in rats [10]. Moreover, interacting levels of miRNA targeted by multiple miRNAs have demonstrated the altered expression of genes in knee OA [11]. Notably, targeting miRNA-mRNA interactions, such as miR-140 with ADAMTS5, MMP13, and IGFBP5 have a great influence on OA pathophysiology [12]. However, limited studies have reported the combinatorial analysis of miRNA-mRNA interaction based on miRNA-mRNA expression profiles in knee OA patients related to synovium pain. This study is dedicated to finding the potential risk biomarkers and predictors of synovium OA pain.

2. Materials and Methods

2.1. Data Source

In the current study, we used miRNA-mRNA dual expression profiling data (GSE99662) to implement our analysis. miRNAs were profiled in knee synovium tissues from knee OA patients with high knee pain (n = 5) and low knee pain (n = 5) by the next generation of sequencing from the study of Sprott H (Table 1). A number of five female normal knee synovium tissues provided miRNAs from the results of Lotz and coworkers [13]. The processing of raw miRNA and mRNA sequencing data was modified according to our in-house analysis pipeline. Significant differential expression mRNAs and miRNAs between patients and normal tissues were distinguished by the DESeq R package [14]. More rigorous criteria, as the cutoff, was used to filter differentially expressed miRNAs and mRNAs. As a result, 2485 differentially expressed mRNAs and 36 differentially expressed miRNAs were identified for further analysis.

2.2. Canonical Correlation Analysis (CCA)

CCA was performed to identify and measure the correlation between two sets of variables to find the association matrix between variables. Notably, if there were too many variables in the analysis, it would be inappropriate to use this analysis method. Only differential expression microRNA and mRNA were selected for target gene prediction, which was achieved by mirWalk2 [15]. In this analysis, miRNAs and mRNAs were taken as the first and second set of variables, respectively. Then, we performed RCC to search for potential associations between miRNAs and mRNAs. The leave-one-out criterion was applied before performing RCC on two sets of variables. R software package (http://www.r-Project.org) was used for CCA.

2.3. Network Construction

For the miRNAs and mRNAs, extracted from the CCA, we extended the network of miRNA and mRNA based on retaining the miRNA-targeting relationship obtained from the three miRNA-target sources. According to the Spearman correlation coefficient, we chose only the miRNA-mRNA interaction with value <0.05 to be significant. In this analysis, there were two reasons why multiple tests were not performed. First, it enabled us to retain a high degree of confidence in the interaction network building. Second, it avoided the loss of information. The constructed network predicted miRNA-targeting relationships that might be related to synovium-associated pain with KOA.

2.4. GO and KEGG Function Enrichment Analysis

For the five miRNAs in the constructed network, we performed Gene Ontology (GO) function and KEGG enrichment analysis on its target genes. ClusterProfiler software R package (http://www.r-project.org) was used for the analysis. Multiple test corrections were performed based on the Benjamini–Hochberg method [16]. A GO and KEGG term with was considered to be significant.

2.5. Verification of Potential Genes Related to Synovium Pain in Knee OA in the Network

2.5.1. Classification Performance Analysis of Potential Synovial-Associated Pain in Knee OA Related Genes Identified

All the ten patients were female, aged 62 to 78 years. OA pain severity status of patients was assessed a day before surgery by a clinical investigator using a 10-point VAS. VAS score of 0–3 indicates low pain and 7–10 high pain. Five patients with high pain, five patients with low pain, and no healthy donors were enrolled. The Local Ethical Committee (KEK-ZH-Nr. 2013-0210) approved the study. Normal human knee cartilage tissues were procured by tissue banks (approved by Scripps Institutional Review Board) from five women. All patients have written informed consent. All experiments were performed in accordance with relevant guidelines and regulations. Standardized surgical excision (excision of the capsule with synovia in the recesses suprapatellar) was performed. The biopsies (two per patient with a size of approximately 2 cm²) were fresh frozen (dry ice). RNA was extracted randomly from one biopsy using TRIzol® Reagent (Thermo Fisher Scientific, Waltham, MA, USA). Then, RNA was analyzed by Qubit® 1.0 Fluorometer (Thermo Fisher Scientific) and Agilent 2200 Tape Station (Agilent Technologies, Santa Clara, CA, USA). At last, RNA with RIN (RNA Integrity Number) values 8.3–9.1 was used for sequencing.

To find out whether the genes in the constructed network is related to synovitis-caused knee pain, we manually chose genes overlapping with online Mendelian Inheritance in Man (OMIM) genetic database [17]. The area under the ROC of overlapping genes was calculated to detect target genes. The expression of these overlapping genes was used to distinguish low pain and high pain by the support vector machine (SVM) method [18]. Besides, cross validation could limit problems such as overfitting to obtain more accurate model prediction performance estimates. Because there were not enough data to divide into “training set” and “test set”, it may affect important modeling or testing capabilities. To solve this problem, we carried out cross validation five times in the SVM program. So we divided all the samples into five groups. For each analysis, one group served as test data and the rest as train data.

2.5.2. Cluster Analysis of miRNA and mRNA Expression Profiles to Validate Possible Genes Related to Synovium Pain in Knee OA

To verify the potential genes involved in the OA-related network, we combined miRNA with the mRNA expression profiles of 10 knee OA patients and used the similar network fusion (SNF) method in the cluster analysis [19]. SNF gathered the similarities between different data types of knee OA patients through a cross-network diffusion process to build the similarity matrix of the fusion patient. Low-pain and high-pain patients with knee OA were divided into clusters, which intended to compare with normal female synovium tissue. The most significant differentially expressed genes between the clusters and the normal samples were extracted. Figure 1 illustrates our workflow flow chart for this analysis.

3. Results

3.1. CCA Results

A total of 6 differentially expressed miRNAs and 425 differentially expressed mRNAs were analyzed by CCA. Because the number of genes was larger than the number of samples, we took the CCA to identify miRNA-mRNA associations from OA synovium pain. Optimal values of two regularization parameters were 0.75 and 0.25 based on the leave-one-out criterion. We obtained four matrixes through regularized canonical correlation: the correlation with miRNA-miRNA canonical variables; the correlation with mRNA-miRNA canonical variables; the correlation with miRNA-mRNA canonical variables; and the correlation with mRNA-mRNA canonical variables. For these four matrixes, we selected the correlation coefficients of the original variable and the first canonical variable because the first canonical variable had the highest canonical correlation coefficient of 0.954 (Appendix 1). We identified 5 miRNAs (miR-133a-3p, miR-145-5p, miR-335-5p, miR-224-5p, and miR-215-5p) in Table 1 and 284 mRNAs, 251 genes, and 260 miRNA-mRNA target pairs with correlation coefficients >0.3 (Appendixes 2 and 3). These miRNAs and mRNAs were further used for network construction.

3.2. Network Construction

The Spearman rank correlation test was used to analyze 13 miRNAs and 26 mRNAs, which were identified by CCA (Appendix 4). Based on the criterion, , we found five miRNAs and 129 genes in the network (Figure 2). From the constructed network, we found that some target genes of miRNAs were associated with synovium OA pain. For example, leptin receptor (LEPR) gene has been found in the network and associated with susceptibility to knee OA [20]. Interesting, our network found that miR-335-5p could regulate LEPR. Transforming growth factor-b type II receptor (TGFBR2) has found to influence cartilage metablism [21], which can be regulated by miR-133a-3p, miR-145-5p, and miR-224-5p.

3.3. Gene Ontology (GO) Function Enrichment Analysis

GO enrichment analysis was performed for each miRNA target gene in the network. Figure 3 shows GO and KEGG function enrichment analysis of target mRNA of miR-133a-3p. Cooler colors indicate more significant GO terms (the values displayed in the graph were not adjusted). As shown in Figure 3, the target genes of miR-133a-3p were enriched on vasculogenesis by GO analysis, MAPK signaling pathway cytokine, or receptor interaction by KEGG analysis.

3.4. Classification Performance Analysis of Genes with Potential Synovium-Related OA Pain

We manually chose genes overlapping with synovial-associated pain in KOA-related genes in online database OMIM. Four important genes (TGFBR2, DST, TBXAS1, FHLI) have been found to be associated with synovium-related OA pain and the area under the curve (AUC) more than 0.8 (Figure 4(a)). The expression of these four genes is considered as the prediction. The classification accuracy rate of TGFBR2, DST, TBXAS1, and FHLI, respectively, was more significant than 0.8. TGFBR2 showed higher classification accuracy rate than other three genes based on SVM classification plot (Figure 4(b)). DST and FHL1 showed excellent classification performance in low pain and were distinguished well from high-pain patients (Figures 4(c) and 4(d)).

(a)

(b)

(c)

(d)

Figure 4

Classification performance analysis of four genes overlapped with synovial pain knee OA genes identified in the OMIM database. (a) ROC analysis of the four genes: TGFBR2, DST, TBXAS1, and FHL1. AUC was above 0.8. Classification accuracy was above 80%. (b) Classification accuracy of the four genes based on the SVM method. (c, d) The SVM classification diagram of TGFBR2, DST, TBXAS1, and FHL1. OMIM: Online Mendelian Inheritance in Man; ROC: receiver operating characteristic curve; AUC: area under the Curve; SVM: support vector machine.

3.5. Cluster Analysis of miRNA and mRNA Expression Profiles to Validate Possible Genes Involved Synovial-Related Knee OA Pain

SNF results suggested that low-pain and high-pain patients with knee OA be divided into three clusters and then compared with five normal female synovium tissues individually. A heat map in Figure 5(a) and a silhouette plot map in Figure 5(b) described the three clusters. Cluster 1 included four patients (P22L, P18H, P17H, and P21L). Cluster 2 included three patients (P14L, P15L, and P12H), and cluster 3 included three patients (P19H, P16H, and P13L). In this analysis, the silhouette score of cluster 2 was the highest (S1 = 0.31), which suggested a greater similarity among samples of the same cluster. We extracted 400 significant differentially expressed genes from the 3 clusters and the normal samples for each of the three clusters. Surprisingly, DST and TBXAS1 were the common genes that existed in the compartment results of all the individual clusters vs. normal. Hence, these hub genes in the network could be considered as potential synovium-related OA pain genes.

(a)

(b)

4. Discussion

This study predicted potential biomarkers of synovial pain in knee OA. At present, it has proved that age, sex, and body weight are risk factors for knee OA [22]. By performing integration data analysis, we can find the potential biomarkers of synovium-associated knee OA pain. Previous researches have shown that the metastasis of cancer was related to the abnormal regulation of miRNA-mRNA [23]. There are few studies to prove the relationship between miRNA and their targets in synovium-associated knee OA pain. Meanwhile, few studies have reported a correlation between miRNA-mRNA maladjustment and synovium-associated pain in knee OA. The reason might be that there is poor data availability on pre-double expression of miRNA-mRNA.

In this study, we performed CCA to find potential biomarkers of synovium-associated knee OA pain. miRNA-mRNA coexpression data were analyzed to ascertain the correlation between miRNA-mRNA and synovial pain in knee OA. CCA applies a nonlinear function to the original variables and extracts related components from two sets of variables [24]. CCA does not require a dataset to have the same dimension and can be used for more than two datasets at the same time. On the other hand, CCA does not assume the directionality of the relationships between datasets. That is contrary to specifying the regression method for both independent and dependent datasets. Besides, CCA not only quantifies the similarity between datasets, but also describes relationships between datasets in an interpretable way [25]. Meanwhile, GO enrichment analysis was performed on target genes in this research. GO is commonly used in the biomedical research community for analysis-histology and related data. It describes functions from three aspects: molecular functions, cellular components, and biological processes [26].

Our research found that TGFBR2, DST, TBXAS1, and FHL1 genes were closely linked to susceptibility to synovial pain in patients with knee OA. TGFBR2 gene can promote the bone formation of osteoblasts and affect the maintenance of bone mass homeostasis. Knockout of TGFBR2 will reduce the severity of OA [27]. Dystonin is a cytoskeletal crosslinker protein that binds to actin and tubulin. This protein dysfunction can lead to sensory neuropathy, known as the human muscle-atrophic sensory autonomic neuropathy [28]. According to our research, this is the first time that dystonin and synovium pain in KOA have correlated. The relationship between the two needs further study. Then, TbXAS1 is a member of cytochrome P450 superfamily, which can directly or indirectly provide angiogenic factors to platelets to promote angiogenesis. TBXAS1 can catalyze the conversion of prostaglandin H2 to thromboxane A2 and then induce platelet aggregation and smooth muscle contraction [29]. FHL1 gene is highly expressed in skeletal muscle and participates in the growth and differentiation of muscle cells. It may regulate cell functions mediated by integrins, such as cytoskeletal rearrangement [30]. We took the expression of these four genes as prediction variables and used a support vector machine to classify the samples. The classification accuracy of TGFBR2, DST, TBXAS1, and FHL1 was 100%, 92%, 80%, and 92%, respectively (Figure 4(a)). Among them, TGFBR2 showed a classification accuracy of 100%. Therefore, we used the SVM classification map to determine the classification performance of the three genes (TGFBR2, DST, and FHL1), which showed a higher classification accuracy than TBXAS1.

In this construction network, we observed five newly identified miRNAs has-miR-133a-3p, has-miR-145-5p, has-miR-335-5p, has-miR-224-5p, and has-miR-215-5p, which can regulate synovial pain in knee OA. Has-miR-133a-3p was involved in osteoclast formation, differentiation, apoptosis, and resorption, showing great promise as a potential therapeutic target for biomarkers and osteoporosis [31], and was found to regulate TGFBR2 and FHL1 expression in this study. Surprisingly, has-miR-145-5p was also detected to regulate TGFBR2 and FHL1, which has been discovered to be associated with changes in cartilage matrix production and degradation levels in patients with OA [32]. Furthermore, has-miR-335-5p can regulate DST and FHL1. Previous studies have shown that has-miR-335-5p expression is elevated in OA patients. MiR-335-5p can regulate the differentiation process of hMSCs chondrocytes induced by transforming growth factor β3. SoMiR-335-5p has a potential role in regulating genes to maintain more progenitor phenotypes [33]. Furthermore, has-miR-224-5p plays a significant role in the pathogenesis of femoral head necrosis and may be a potential biomarker for the prognosis of it [34]. Interestingly, we found that this microRNA also regulated the expression of DST, TBXAS1, and FHL1, suggesting that it may be related to synovial pain in KOA.

Meanwhile, there has been no evidence of a link between has-miR-215-5p and KOA. In this research, we studied that has-miR-215-5p regulated the expression of FHL1. However, these miRNA target relationships obtained require further validation in molecular biology experiments.

SNF shows its ability to integrate data when using different types of data measurements at the same time [35]. Fortunately, we used the SNF analysis to find that DST and TBXAS1 genes exist in all 3 clusters in our study. These results suggested that the two genes may play a key role in the progress of synovium pain in knee OA.

There are three limitations in our study. First, we have selected only verified targeted genes by experiment in the database. Some predicted miRNA target regulation had not been incorporated into the analysis. Second, some exciting findings were based on a small sample. We plan to increase the sample in the next study. Third, synovial fluid exosomal miRNA content has been found altered in patients with OA, and these changes were gender-specific [36]. The samples’ data we used were of only female, so the research result cannot be translated into a general OA population, especially male. The research data should be interpreted carefully.

These predictive analyses revealed that some miRNAs associated with knee OA were involved in synovium pain-related genes, especially DST and TBXAS1. Therefore, our consequences can be used as a reference for future explorers on the pathogenesis and genetic causes of synovial pain in knee OA. In future studies, it is necessary to explore the potential mechanism of controlling complex molecular regulation by integrating multiple data types. These studies will be helpful for further experimental verification and discovery of knee OA biomarkers.

Data Availability

Most of the data used to support the findings of this study are included within the supplementary information files. Other accurate data could available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

HMW and CZ conceptualized and designed the article. LW, YH, MMH, HMW, and YJX analyzed and interpreted the data. JXW, YH, and CZ drafted the article. HMW and JXW were responsible for critical revision of the article. HMW and CZ finally approved the article.

Acknowledgments

We are grateful to Chengdu Basebiotech Co. Ltd. for providing bioinformatic analysis assistance. This study was supported by the Natural Science Foundation of China (Grant no. 81802246) and Research Fund for the Doctoral Program of Affiliated Hospital of Southwest Medical University (Project no. 2016-QB-4).

Supplementary Materials

Appendix 1: four matrixes. We selected the correlation coefficients of the original variable and the first canonical variable because the first canonical variable had the highest canonical correlation coefficient of 0.954. Appendixes 2 and 3: the miRNAs and mRNAs selected in this study, which will be further used in network construction. Appendix 4: 13 miRNAs and 26 mRNAs identified by CCA, which were analyzed by the Spearman rank correlation test. (Supplementary Materials)

References

C. R. Scanzello, “Role of low-grade inflammation in osteoarthritis,” Current Opinion in Rheumatology, vol. 29, no. 1, pp. 79–85, 2017.
View at: Publisher Site | Google Scholar
L. A. Stoppiello, P. I. Mapp, D. Wilson, R. Hill, B. E. Scammell, and D. A. Walsh, “Structural associations of symptomatic knee osteoarthritis,” Arthritis & Rheumatology, vol. 66, no. 11, pp. 3018–3027, 2014.
View at: Publisher Site | Google Scholar
M. W. Lark, E. K. Bayne, J. Flanagan et al., “Aggrecan degradation in human cartilage. Evidence for both matrix metalloproteinase and aggrecanase activity in normal, osteoarthritic, and rheumatoid joints,” Journal of Clinical Investigation, vol. 100, no. 1, pp. 93–106, 1997.
View at: Publisher Site | Google Scholar
R. C. Billinghurst, L. Dahlberg, M. Ionescu et al., “Enhanced cleavage of type II collagen by collagenases in osteoarthritic articular cartilage,” Journal of Clinical Investigation, vol. 99, no. 7, pp. 1534–1545, 1997.
View at: Publisher Site | Google Scholar
M. R. Hussein, N. A. Fathi, A. M. E. El-Din et al., “Alterations of the CD4⁺, CD8⁺ T cell subsets, interleukins-1β, IL-10, IL-17, tumor necrosis factor-α and soluble intercellular adhesion molecule-1 in rheumatoid arthritis and osteoarthritis: preliminary observations,” Pathology & Oncology Research, vol. 14, no. 3, pp. 321–328, 2008.
View at: Publisher Site | Google Scholar
A. Fragoulis, J. Laufs, S. Müller et al., “Sulforaphane has opposing effects on TNF-alpha stimulated and unstimulated synoviocytes,” Arthritis Research & Therapy, vol. 14, no. 5, p. R220, 2012.
View at: Publisher Site | Google Scholar
R. A. Black, C. T. Rauch, C. J. Kozlosky et al., “A metalloproteinase disintegrin that releases tumour-necrosis factor-α from cells,” Nature, vol. 385, no. 6618, pp. 729–733, 1997.
View at: Publisher Site | Google Scholar
U. Pankiewicz, M. Sujka, M. Włodarczyk-Stasiak, A. Mazurek, and J. Jamroz, “Effect of pulse electric fields (PEF) on accumulation of magnesium and zinc ions in Saccharomyces cerevisiae cells,” Food Chemistry, vol. 157, no. 1, pp. 125–131, 2014.
View at: Publisher Site | Google Scholar
F. Linda, K. R. B. Hagen, J. W. J. Bijlsma et al., “EULAR recommendations for the non-pharmacological core management of hip and knee osteoarthritis,” Annals of the Rheumatic Diseases, vol. 72, no. 7, pp. 1125–1135, 2013.
View at: Publisher Site | Google Scholar
J. W. J. Bijlsma and K. Knahr, “Strategies for the prevention and management of osteoarthritis of the hip and knee,” Best Practice & Research Clinical Rheumatology, vol. 21, no. 1, pp. 59–76, 2007.
View at: Publisher Site | Google Scholar
M. Cross, E. Smith, D. Hoy et al., “The global burden of hip and knee osteoarthritis: estimates from the global burden of disease 2010 study,” Annals of the Rheumatic Diseases, vol. 73, no. 7, pp. 1323–1330, 2014.
View at: Publisher Site | Google Scholar
T. D. Spector and A. J. Macgregor, “Risk factors for osteoarthritis: genetics,” Osteoarthritis & Cartilage, vol. 12, pp. 39–44, 2004.
View at: Publisher Site | Google Scholar
K. M. Fisch, R. Gamini, O. Alvarez-Garcia et al., “Identification of transcription factors responsible for dysregulated networks in human osteoarthritis cartilage by global gene expression analysis,” Osteoarthritis and Cartilage, vol. 26, no. 11, pp. 1531–1538, 2018.
View at: Publisher Site | Google Scholar
Y.-D. Zhang, S.-C. Zhao, Z.-S. Zhu et al., “Cx43- and smad-mediated TGF-β/BMP signaling pathway promotes cartilage differentiation of bone marrow mesenchymal stem cells and inhibits osteoblast differentiation,” Cellular Physiology and Biochemistry, vol. 42, no. 4, pp. 1277–1293, 2017.
View at: Publisher Site | Google Scholar
T. Kawasaki, Y. Niki, T. Miyamoto et al., “The effect of timing in the administration of hepatocyte growth factor to modulate BMP-2-induced osteoblast differentiation,” Biomaterials, vol. 31, no. 6, pp. 1191–1198, 2010.
View at: Publisher Site | Google Scholar
V. B. Kraus, S. Feng, S. C. Wang et al., “Trabecular morphometry by fractal signal analysis is a novel marker of osteoarthritis progression,” Arthritis Rheum, vol. 60, no. 12, pp. 3711–3722, 2009.
View at: Publisher Site | Google Scholar
A. Hamosh, A. F. Scott, J. Amberger, D. Valle, and V. A. Mckusick, “Online Mendelian Inheritance in Man (OMIM),” Human Mutation, vol. 15, no. 1, pp. 57–61, 2015.
View at: Publisher Site | Google Scholar
L. Du, C. Zhang, Q. Liu, X. Zhang, and B. Yue, “Krait: an ultrafast tool for genome-wide survey of microsatellites and primer design,” Bioinformatics, vol. 34, no. 4, pp. 681–683, 2017.
View at: Publisher Site | Google Scholar
B. Wang, A. M. Mezlini, F. Demir et al., “Similarity network fusion for aggregating data types on a genomic scale,” Nature Methods, vol. 11, no. 3, pp. 333–337, 2014.
View at: Publisher Site | Google Scholar
J. M. Wilkins, J. Loughlin, and S. J. B. Snelling, “Osteoarthritis genetics: current status and future prospects,” Future Rheumatology, vol. 2, no. 6, pp. 607–620, 2007.
View at: Publisher Site | Google Scholar
S. Ikegawa, “New gene associations in osteoarthritis: what do they provide, and where are we going?” Current Opinion in Rheumatology, vol. 19, no. 5, pp. 429–434, 2007.
View at: Publisher Site | Google Scholar
S. Lee and S. J. Kim, “Prevalence of knee osteoarthritis, risk factors, and quality of life: the fifth Korean national health and nutrition examination survey,” International Journal of Rheumatic Diseases, vol. 20, no. 7, pp. 809–817, 2017.
View at: Publisher Site | Google Scholar
A. Jacobsen, J. Silber, G. Harinath, J. T. Huse, N. Schultz, and C. Sander, “Analysis of microRNA-target interactions across diverse cancer types,” Nature Structural & Molecular Biology, vol. 20, no. 11, pp. 1325–1332, 2013.
View at: Publisher Site | Google Scholar
W. W. Hsieh, “Nonlinear canonical correlation analysis by neural networks,” Neural Networks, vol. 13, no. 10, pp. 1095–1105, 2000.
View at: Publisher Site | Google Scholar
N. Y. Bilenko and J. L. Gallant, “Pyrcca: regularized kernel canonical correlation analysis in Python and its applications to neuroimaging,” Frontiers in Neuroinformatics, vol. 10, p. 49, 2016.
View at: Publisher Site | Google Scholar
T. Yang, J. Zhang, Y. Cao et al., “Wnt5a/Ror2 mediates temporomandibular joint subchondral bone remodeling,” Journal of Dental Research, vol. 94, no. 6, p. 803, 2015.
View at: Publisher Site | Google Scholar
G. Zhen, C. Wen, X. Jia et al., “Inhibition of TGF-beta signaling in mesenchymal stem cells of subchondral bone attenuates osteoarthritis,” Nature Medicine, vol. 19, no. 6, pp. 704–712, 2013.
View at: Publisher Site | Google Scholar
S. F. Kornfeld, A. Lynch-Godrei, S. R. Bonin, S. Gibeault, Y. De Repentigny, and R. Kothary, “Cytoskeletal linker protein dystonin is not critical to terminal oligodendrocyte differentiation or CNS myelination,” PLoS One, vol. 11, no. 2, Article ID e0149201, 2016.
View at: Publisher Site | Google Scholar
S. H. Oh, Y. H. Kim, S. M. Park et al., “Association analysis of thromboxane a synthase 1 gene polymorphisms with aspirin intolerance in asthmatic patients,” Pharmacogenomics, vol. 12, no. 3, pp. 351–363, 2011.
View at: Publisher Site | Google Scholar
H. Joos, W. Albrecht, S. Laufer, H. Reichel, and R. E. Brenner, “IL-1beta regulates FHL2 and other cytoskeleton-related genes in human chondrocytes,” Molecular Medicine, vol. 14, no. 3-4, pp. 150–159, 2008.
View at: Publisher Site | Google Scholar
X. Ji, X. Chen, and X. Yu, “MicroRNAs in osteoclastogenesis and function: potential therapeutic targets for osteoporosis,” International Journal of Molecular Sciences, vol. 17, no. 3, p. 349, 2016.
View at: Publisher Site | Google Scholar
L. A. Vonk, A. H. Kragten, W. J. Dhert, D. B. Saris, and L. B. Creemers, “Overexpression of hsa-miR-148a promotes cartilage production and inhibits cartilage degradation by osteoarthritic chondrocytes,” Osteoarthritis Cartilage, vol. 22, no. 1, pp. 145–153, 2014.
View at: Publisher Site | Google Scholar
M. Kopanska, D. Szala, J. Czech et al., “MiRNA expression in the cartilage of patients with osteoarthritis,” Journal of Orthopaedic Surgery and Research, vol. 12, no. 1, p. 51, 2017.
View at: Publisher Site | Google Scholar
J. Moreno-Rubio, G. Herrero-Beaumont, L. Tardio, M. A. Alvarez-Soria, and R. Largo, “Nonsteroidal antiinflammatory drugs and prostaglandin E₂ modulate the synthesis of osteoprotegerin and RANKL in the cartilage of patients with severe knee osteoarthritis,” Arthritis Rheum, vol. 62, no. 2, pp. 478–488, 2010.
View at: Publisher Site | Google Scholar
A. N. Zizzo, L. Erdman, B. M. Feldman, and A. Goldenberg, “Similarity network fusion: a novel application to making clinical diagnoses,” Rheumatic Disease Clinics of North America, vol. 44, no. 2, pp. 285–293, 2018.
View at: Publisher Site | Google Scholar
R. Kolhe, M. Hunter, S. Liu et al., “Gender-specific differential expression of exosomal miRNA in synovial fluid of patients with osteoarthritis,” Scientific Reports, vol. 7, no. 1, p. 2029, 2017.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2019 Haiming Wang et al. 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.

PDF Download Citation

Download other formats

Order printed copies