PDF

Int J Med Sci 2019; 16(8):1057-1071. doi:10.7150/ijms.35611 This issue Cite

Research Paper

Identification of Novel Genes in Osteoarthritic Fibroblast-Like Synoviocytes Using Next-Generation Sequencing and Bioinformatics Approaches

Yi-Jen Chen^1,2, Wei-An Chang^1,3, Ling-Yu Wu¹, Ching-Fen Huang^1,2, Chia-Hsin Chen^2,4,5 , Po-Lin Kuo^1,6

1. Graduate Institute of Clinical Medicine, College of Medicine, Kaohsiung Medical University, Kaohsiung 807, Taiwan
2. Department of Physical Medicine and Rehabilitation, Kaohsiung Medical University Hospital, Kaohsiung 807, Taiwan
3. Division of Pulmonary and Critical Care Medicine, Kaohsiung Medical University Hospital, Kaohsiung 807, Taiwan
4. Department of Physical Medicine and Rehabilitation, School of Medicine, College of Medicine, Kaohsiung Medical University, Kaohsiung 807, Taiwan
5. Orthopaedic Research Center, Kaohsiung Medical University, Kaohsiung 807, Taiwan
6. Center for Cancer Research, Kaohsiung Medical University

✉ Corresponding author: Chia-Hsin Chen; chchenedu.tw; Tel.: +886-7-312-1101 ext. 5962 Po-Lin Kuo; kuopolinnet.tw; Tel.: +886-7-312-1101 ext. 2512-33

Abstract

Synovitis in osteoarthritis (OA) the consequence of low grade inflammatory process caused by cartilage breakdown products that stimulated the production of pro-inflammatory mediators by fibroblast-like synoviocytes (FLS). FLS participate in joint homeostasis and low grade inflammation in the joint microenvironment triggers FLS transformation. In the current study, we aimed to identify differentially expressed genes and potential miRNA regulations in human OA FLS through deep sequencing and bioinformatics approaches. The 245 differentially expressed genes in OA FLS were identified, and pathway analysis using various bioinformatics databases indicated their enrichment in functions related to altered extracellular matrix organization, cell adhesion and cellular movement. Moreover, among the 14 dysregulated genes with potential miRNA regulations identified, src kinase associated phosphoprotein 2 (SKAP2), adaptor related protein complex 1 sigma 2 subunit (AP1S2), PHD finger protein 21A (PHF21A), lipoma preferred partner (LPP), and transcription factor AP-2 alpha (TFAP2A) showed similar expression patterns in OA FLS and OA synovial tissue datasets in Gene Expression Omnibus database. Ingenuity Pathway Analysis identified the dysregulated LPP participated in cell migration and cell spreading of OA FLS, which was potentially regulated by miR-141-3p. The current findings suggested new perspectives into understanding the novel molecular signatures of FLS involved in the pathogenesis of OA, which may be potential therapeutic targets.

Keywords: osteoarthritis, synovitis, fibroblast-like synoviocytes, next-generation sequencing, messenger RNA, microRNA, bioinformatics

Introduction

Osteoarthritis (OA) is one of the common articular disorders that affect major weight bearing joints, causing joint pain and stiffness and lead to chronic disability [1]. The increasing prevalence of OA is likely due to increases in longevity and prevalence of obesity [2,3]. Clinically, the diagnosis of OA is mainly based on symptoms and radiographic findings, although discordance between pain and severity of radiographic joint pathology has been reported [4,5]. The major histopathological changes in OA joint are the cartilage destruction with hypertrophic differentiation of chondrocytes [6]. However, the contribution of low-grade inflammation and synovitis in OA progression has been appreciated, and OA is now considered a disease of the whole joint, not merely the cartilage [7,8].

The synovium forms the boundary between internal joint structure and adjacent soft tissues, and is essential for maintaining joint homeostasis. The major cellular components of this distinct tissue layer are the fibroblast-like and macrophage-like synoviocytes. The fibroblast-like synoviocytes (FLS) produce major constituents of synovial fluid that nourishes the chondrocytes through synovial vascular network, while the synovial macrophages help clearing debris from minor joint injuries [9]. Synovitis is a common feature of inflammatory arthritis, including rheumatoid arthritis (RA) and OA, and the degree of synovitis is associated with joint pain and structural progression [7,10]. The low-grade inflammatory OA joint microenvironment is caused by the cartilage breakdown products that provoke the release of proteolytic enzymes and increased production of pro-inflammatory mediators from FLS, followed by immune cell infiltration and vascular hyperplasia, leading to synovial inflammation [7,8]. This synovial change and overexpression of pro-inflammatory mediators can be observed in the early stage of OA, even before the presence of macroscopic cartilage degeneration [11].

OA being a disease of the whole joint involving the cartilage, synovium and subchondral bone, and synovitis associated with symptoms and progression of OA, the synovium may serve as a potential therapeutic target in the management of OA [8]. While several studies focusing on OA synovial fluid and FLS have proposed the role of microRNAs (miRNAs) in the pathogenesis of OA synovitis and disease progression [12-14], novel therapeutics targeting these small non-coding single-stranded RNAs through intra-articular injection may contribute to the maintenance of joint homeostasis, fine tuning downstream gene expressions related to inflammatory and catabolic pathways [14-16]. The transcriptome changes and novel molecular signatures between normal and arthritic pathologies can be efficiently identified using the high-throughput next-generation sequencing (NGS) technique [17,18], and the biological themes underlying the differentially expressed genomic profiling can be determined through the integrated analysis with bioinformatics approaches [19-21].

In the current study, the biological functions underlying the differentially expressed genes and potential miRNA regulations in OA FLS will be investigated using NGS and different bioinformatics databases, and validated in clinical OA synovium tissue data available in functional genomics data repository. We propose the findings will gain novel insights into understanding the role of FLS in the pathogenesis of OA and identify potential therapeutic targets in the management of OA.

Materials and Methods

Culturing Human Fibroblast-Like Synoviocytes (HFLS)

Human fibroblast-like synoviocytes isolated from adult normal (HFLS) and osteoarthritic synovial tissue (HFLS-OA) were obtained from Cell Applications, Inc. (San Diego, CA, USA). The isolated cells were cryopreserved at the first passage. The cryopreserved vials of HFLS and HFLS-OA were thawed and cultured in Synoviocyte Growth Medium (Cell Applications, Inc. San Diego, CA, USA) and incubated in a 37°C, 5% CO₂ humidified incubator until confluence. The cells were then harvested for total RNA extraction using Trizol Reagent (Invitrogen, Carlsbad, CA, USA). The quality of extracted RNAs were confirmed using ND-1000 spectrophotometer (Nanodrop Technology, Wilmington, DE, USA) for detection of OD₂₆₀/OD₂₈₀ absorbance ratio and Bioanalyzer 2100 (Agilent Technology, Santa Clara, CA, USA) for RNA integrity number (RIN) with RNA 6000 labchip kit (Agilent Technology, Santa Clara, CA, USA). The OD₂₆₀/OD₂₈₀ absorbance ratio was 1.95 for HFLS and 1.94 for HFLS-OA, while the RINs were 9.9 and 10 for HFLS and HFLS-OA, respectively, indicating good quality of the extracted RNA.

RNA Sequencing

The RNA and small RNA sequencing were carried out by Welgene Biotechnology Company (Welgene, Taipei, Taiwan). For RNA sequencing, all RNA samples were prepared according to the Illumina protocol. The Agilent's SureSelect Strand Specific RNA Library Preparation Kit was used for RNA library construction, followed by AMPure XP Beads size selection. The sequence was determined by sequencing-by-synthesis technology, with read length at 150 nucleotides pair-end. The sequence data was generated by Welgene's pipeline based on Illumina bcl2fastq v2.1.7. The raw reads were trimmed for qualified reads and remove lower quality bases using Trimmomatic (version 0.32), and the qualified reads were then aligned to reference human genome using HISAT2 alignment tool. The expression level of each aligned gene was normalized and expressed in fragments per kilobase of transcript per million mapped reads (FPKM). The differential expression between HFLS and HFLS-OA were analyzed based on Cuffdiff (Cufflinks version 2.2.1) with genome bias detection/correction and Welgene in-house programs. For small RNA sequencing, samples were prepared using Illumina sample preparation kit following the TruSeq Small RNA Sample Preparation Guide. The RNAs were reversed transcribed to cDNA, size-fractionated and purified to obtain bands with 18-40 nucleotides. The sequencing with read length at 75 nucleotides single-end was carried out on Illumina instrument and processed with Illumina software. The raw reads were trimmed for qualified reads and analyzed using miRDeep2 to clip 3' adaptor sequence before aligning to reference human genome from University of California, Santa Cruz (UCSC). The expression levels of known miRNAs were estimated using miRDeep2, normalized in reads per million (RPM). The selection criteria for differentially expressed mRNAs and miRNAs between HFLS and HFLS-OA were as following: fold change > 2.0, FPKM > 0.3 for mRNA and RPM > 1 for miRNA in at least one group.

Functional Enrichment Analysis Using Different Bioinformatics Tools

The gene lists of interest were uploaded into Database for Annotation, Visualization and Integrated Discovery (DAVID) bioinformatics resource [22] and Ingenuity Pathway Analysis (IPA) software (Ingenuity systems, Redwood City, CA, USA) [23] to perform integrated data mining and categorize large gene lists into different enriched biological functions and/or networks. The IPA software was also able to predict potential upstream regulators and downstream effectors of a given gene list. In the DAVID database, differentially expressed genes were uploaded for functional annotation analysis, setting the Expression Analysis Systematic Explorer (EASE) score at default cutoff value of 0.1, which represented the modified Fisher's exact p-value. In the IPA software, differentially expressed genes with fold changes between HFLS and HFLS-OA were uploaded for core analysis. The analytic results were obtained based on all direct and indirect relationships identified in all tissue types, and from either experimentally observed or moderate to highly predicted confidence.

Protein-Protein Interaction Network Analysis Using STRING Database

To identify the protein-protein interaction (PPI) network of differentially expressed genes, the STRING database (version 11.0) integrating functional interactions from known and predicted protein-protein association data was used [24]. For sub-network analysis, the Molecular Complex Detection (MCODE) plugin tool under Cytoscape software package was used to cluster the large PPI network into small networks [25].

MiRNA Target Prediction

For those identified differentially expressed miRNAs between HFLS and HFLS-OA, the putative targets were predicted using the miRmap database (miRmap version 1.0), an open-source software library that was developed using a comprehensive approach to predict the repression strength of a miRNA to specific genes [26]. Higher miRmap scores indicated higher repression strength. In the current study, 83 differentially expressed miRNAs were analyzed for their putative targets, and those putative targets with miRmap scores higher than 99.0 were selected. In addition, those potential miRNA-mRNA interactions of interest were further validated in other two miRNA prediction databases, including TargetScan [27] and miRDB [28].

Functional Genomics Data Repository -- Gene Expression Omnibus (GEO) Database

To assess the expression patterns of candidate genes of interest in clinical OA synovial tissue samples, we searched in the GEO database [29] for related high-throughput genomic datasets on synovial tissues from normal and OA patients. The genes of interest with their expression values could be obtained for further between-group comparison. In the current study, we assessed the expression patterns of candidate genes in five datasets of normal and OA synovial tissue samples (GSE55235, GSE55457, GSE82107, GSE1919 and GSE29746) and one dataset comparing non-inflammatory and inflammatory OA synovial tissues (GSE46750).

Statistical Analysis

The between-group difference of target gene expression values identified from selected GEO datasets were analyzed using non-parametric Mann-Whitney U test with SPSS Statistics software (version 19, IBM Corp., Armonk, NY, USA). A p-value < 0.05 was considered statistically significant.

Results

Identification of Differential Expression Profile between Normal and Osteoarthritic Human Fibroblast-Like Synoviocytes

The transcriptomic profile of adult HFLS and HFLS-OA cells were obtained from NGS results and the FPKM performance between two samples were displayed in density plot, as shown in Figure 1A. The differentially expressed genes between HFLS and HFLS-OA were screened for according to the following selection criteria: expression higher than 0.3 FPKM in either sample, at least two-fold change between HFLS and HFLS-OA, and significant differential expression with p-value < 0.05. The distribution of differential expression genes between HFLS and HFLS-OA were displayed in volcano plot (Figure 1B). The selection criteria yielded a total of 118 significantly up-regulated genes and 127 significantly down-regulated genes in HFLS-OA cells.

The Differentially Expressed Genes were Enriched in Functions Related to Extracellular Matrix Organization, Cell Adhesion and Cellular Movement

All 245 differentially expressed genes were uploaded into DAVID database for terms of biological process in Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway. In addition, these differentially expressed genes were also input into FunRich database for functional enrichment analysis. The functionally enriched biological processes, KEGG pathways and biological pathways with their p-values were shown in Figure 2. The top enriched functions were related to extracellular matrix (ECM) organization (p = 9.92x10^-6) and cellular movement such cell adhesion (p = 0.007) and epithelial-to-mesenchymal transition (p = 0.002). Five genes were also found to be associated with “response to mechanical stimulus” from the DAVID database (p = 0.007), including COL3A1, CHI3L1, POSTN, ASNS and CITED2, which were all down-regulated in HFLS-OA. Moreover, genes with differential expression values and fold-changes were also uploaded into IPA for core analysis. The results showed “cellular movement” was the top enriched molecular and cellular function, with 29 related molecules involved. Besides, the function annotation of “cell spreading” (p = 0.00238, z-score = -2.121) was predicted to have decreased activation, with the following molecules involved: CAP1, CDH11, LPP, MYH10, SERPINE1, SMAD4, SPARC, TGFBI.

Identification of Enriched Functions in Differentially Expressed Gene Clusters of HFLS

To identify gene clusters among the 245 differentially expressed genes in HFLS-OA and their associated biological functions, the list of differentially expressed genes were input into the STRING database to obtain a large PPI network. The sub-cluster analysis was performed under the Cytoscape software with plug-in tool MCODE. The sub-clusters of networks from MCODE were listed in Table 1, with cluster 1 and cluster 2 having higher scores. The two clusters of sub-networks were drawn in the Cytoscape software, as shown in Figure 3. To understand the biological functions of these two clusters of genes, the two clusters were separately input into DAVID database for functional annotation analysis. The top enriched biological functions in terms of biological process and KEGG pathway were listed in Table 2. Genes in cluster 1 were associated with RNA and protein processing, while genes in cluster 2 were associated with ECM organization and cell focal adhesion.

 Figure 1 

(A) The gene expression from next-generation sequencing in fragments per kilobase of transcript per million mapped reads (FPKM) performance of normal (HFLS) and osteoarthritic (HFLS-OA) human fibroblast-like synoviocytes were displayed in density plot. (B) The differential expression patterns between HFLS and HFLS-OA were plotted in volcano plot. The red dots represented up-regulated genes and the green dots represented down-regulated genes in HFLS-OA. Those genes with fold changes > 2.0 and p value < 0.05 were selected as significantly dysregulated genes.

 Figure 2 

The top enriched (A) biological processes in Gene Ontology terms and (B) Kyoto Encyclopedia of Genes and Genomes pathways in differentially expressed genes of HFLS-OA were identified from the Database for Annotation, Visualization and Integrated Discovery bioinformatics resource. The color scale indicated the corresponding p values and the x-axis indicated the gene counts of each biological function. (C) The enriched biological pathway in differentially expressed genes of HFLS-OA were identified from the FunRich database, where the percentage of genes and -log(p-value) of each biological pathway were indicated.

 Figure 3 

The potential interaction networks of (A) cluster 1 containing 13 molecules and (B) cluster 2 containing 12 molecules identified from Molecular Complex Detection (MCODE) were indicated. The two sub-networks were drawn from the Cytoscape software.

Identification of Differentially Expressed miRNAs and Potential miRNA-mRNA Interactions in HFLS-OA Cells

The differential miRNA expression profile between HFLS and HFLS-OA were simultaneously investigated with small RNA sequencing. The selection criteria for differentially expressed miRNAs in HFLS-OA were as following: normalized read counts > 1 RPM, at least 2.0-fold-change between HFLS and HFLS-OA. The result yielded 43 up-regulated and 40 down-regulated miRNAs in HFLS-OA. To obtain putative targets of dysregulated miRNAs, miRmap database, a miRNA target prediction database, was used, and those predicted targets with miRmap scores of at least 99.0 were selected. There were 956 putative targets of 43 up-regulated miRNAs and 1282 putative targets of 40 down-regulated miRNAs identified. These putative targets of up- and down-regulated miRNAs were matched to our differential expression mRNA profiles of 127 down- and 118 up-regulated genes in HFLS-OA. The heatmaps of differentially expressed miRNAs and mRNAs in z-score and the Venn diagram were shown in Figure 4. A total of 14 target genes with potential miRNA regulations were selected. The detailed gene names and their expression values in FPKM were listed in Table 3.

 Figure 4 

The differentially expressed miRNAs and mRNAs in HFLS and HFLS-OA displayed in heatmaps were indicated in left and right panels, respectively. Putative targets of dysregulated miRNAs were predicted from the miRmap database, selecting those with miRmap scores of ≥ 99.0 indicating high repression strength. The putative targets were matched to differentially expressed mRNAs in HFLS, and the Venn diagram was displayed in the middle panel. A total of 11 up-regulated genes and 3 down-regulated genes with potential miRNA regulations were identified.

Analysis of Target Genes with Potential miRNA-mRNA Interactions in Osteoarthritic Synovial Tissues and Identification of Potential Molecular Signatures in Osteoarthritic Synovium

To validate the expression patterns of these 14 target genes in clinical OA synovial tissues, we searched in the GEO database for OA synovial tissue datasets to further analysis of the expression patterns. Those datasets containing both normal and OA synovial tissue samples were selected for expression analysis. There were four OA synovial tissue datasets (GSE55235, GSE55457, GSE82107 and GSE1919) and one OA synovial fibroblast dataset (GSE29746) found in the database. In addition, one dataset comparing non-inflammatory and inflammatory OA synovial tissue expression profile (GSE46750) was also selected for analysis. The expression levels of the 14 target genes were analyzed in these 6 datasets to search for similar expression patterns found in our HFLS-OA data. The expression patterns of these target genes in the 6 datasets were summarized in Table 4. The more consistently dysregulated expression patterns in SKAP2, AP1S2, PHF21A and LPP were found in our HFLS-OA dataset and in at least two of the four OA synovial tissue datasets. Moreover, the down-regulated LPP was also observed in synovial fibroblast dataset. Additionally, the up-regulated TFAP2A was also found up-regulated in inflammatory OA synovial tissue samples. The expression patterns of the 14 target genes in one of the representative datasets (GSE55235) was shown in Figure 5.

 Figure 5 

The expression patterns of 14 target genes, including (A) 11 up-regulated and (B) 3 down-regulated genes were analyzed in one of the OA synovial tissue datasets extracted from GEO database (GSE55235). The significantly up-regulated expressions in SKAP2, KIF1B, TFAP2A, CACNA2D1 and AP1S2, and significantly down-regulated expressions in PHF21A and LPP in OA knee synovial tissues were in similar expression patterns to our HFLS-OA data. * indicated p < 0.05, ** indicated p < 0.01, *** indicated p < 0.001, and n.s. indicated no statistical significance between normal and OA groups.

Determination of Potential miRNA-mRNA Interactions in HFLS-OA

Since the expression patterns of SKAP2, AP1S2, PHF21A, LPP, and TFAP2A were more consistently observed in our HFLS-OA NGS dataset and OA synovial tissue datasets from GEO database, we further analyzed in the miRmap database for potential miRNA regulations of these candidate genes. Those potential miRNA-mRNA interactions with miRmap scores higher than 99.0 were selected and matched to our HFLS dataset of differentially expressed miRNAs. A total of 11 potential miRNA-mRNA interactions were identified. We then validated the putative 3'UTR binding sites and sequences of these potential miRNA-mRNA interactions in TargetScan and miRDB miRNA prediction databases. The results were listed in Table 5, and there were four potential miRNA regulations consistently validated in all three miRNA target prediction databases, including hsa-miR-450b-5p-SKAP2, hsa-miR-204-5p-AP1S2, hsa-miR-766-3p-PHF21A and hsa-miR-141-3p-LPP.

To understand the association of these miRNA targets among the two main clusters of differentially expressed genes in HFLS-OA, we uploaded these five target genes, SKAP2, AP1S2, PHF21A, LPP, and TFAP2A along with the two clusters of genes into the STRING database for interaction network analysis. The interaction network drawn from the STRING database was shown in Figure 6, where LPP was the only molecule having close association with ACTN1, one of the molecules in cluster 2.

LPP was Potentially Associated with Cell Migration and Cell Spreading

To understand the potential biological themes among the 14 target genes with potential miRNA regulations, these target genes were uploaded to the IPA software for functional enrichment analysis. The top networks from the IPA core analysis result indicated 8 of the 14 target genes were clustered in the top scored network related to diseases and functions of cellular development, cellular growth and proliferation, and neurological disease, in which LPP was one of the molecules in this top network. The detailed networks of these target genes were listed in Table 6. The interactions between the molecules in network 1 was shown in Figure 7, where tumor protein 53 (TP53) was the central hub of the network. Further overlay diseases and functions analysis indicated TP53, PAK3, MYC, LPP, and CYR61 (marked in purple frame in Figure 7) were associated with “migration of fibroblasts”. Along with previous finding from IPA that LPP was one of the molecules predicted to be involved in cell spreading, it is suggested that LPP potentially regulated by miR-141-3p was involved in functions of cell migration and cell spreading.

 Figure 6 

The five putative targets with potential miRNA regulations, including SKAP2, AP1S2, PHF21A, LPP, and TFAP2A, along with the two clusters of genes previously identified were input into STRING database for potential interaction network. Among the 5 putative targets, we found that LPP (indicated in black arrow) was the only molecule having direct interaction with ACTN1 in cluster 2, and indirect interaction with LMO7 in cluster 1.

 Figure 7 

The top scored network of the 14 target genes identified from the IPA software indicated 8 of the 14 genes were group in this network related to cellular development, cellular growth and proliferation, and neurological disease. The overlay diseases and functions analysis indicated TP53, PAK3, MYC, LPP and CYR61 (marked in purple frames) were associated with migration of fibroblasts. Molecules in red indicated up-regulated expression and molecules in green indicated down-regulated expression in HFLS-OA data. The color scales indicated the relative expression values of HFLS-OA to HFLS. The numbers indicated below each colored molecule indicated fold-changes and log₂(ratio) of HFLS-OA to HFLS expression.

 Figure 8 

Schematic summary of novel molecular signatures in OA fibroblast-like synoviocytes

Discussion

Through NGS and bioinformatics analysis, our current study identified the differentially expressed genes in HFLS-OA were enriched in functions related to ECM organization, cell adhesion and cellular movement. Moreover, the dysregulated LPP was suggested to participate in altered cellular movement of OA synoviocytes, potentially regulated by miR-141-3p, systematically validated through different miRNA prediction databases. The schematic summary of the novel molecular signatures identified in our current study is shown in Figure 8.

Synovitis is one of the macroscopic structural changes of OA, affecting joint integrity [30]. FLS transformation and synovitis have been reported to contribute to inflammatory arthritis [9,31]. In response to inflammation, FLS transformation occurs through activation of several signaling pathways, for instance, toll-like receptor signaling and epigenetic control [9]. Transformed FLS become aggressive in behavior and lead to persistent synovitis; they can proliferate rapidly and are highly migratory, and attach to the articular cartilage and produce matrix degrading enzymes to degrade and invade the cartilage [9]. The aggressive behavior of FLS in synovitis of RA and OA may differ in extend. In RA, diffuse synovial inflammation is observed, while synovitis in OA is usually in patchy distribution, confined to sites adjacent to cartilage damage [30]. The increased migratory and invasive ability of FLS have been more extensively recognized in RA than in OA [9,32], and several transcription factors are reported to be potentially responsible for the FLS transformation in inflammatory arthritic diseases [33,34]. In OA, several studies suggested that activation of receptors for advanced glycation end products and up-regulated expression of chemokine receptor 3 and cadherin-11 could increase the catabolic activity and migratory/invasive capacity of FLS [35-37]. In our HFLS-OA data, bioinformatics analysis suggested top enriched functions related to ECM organization and altered cellular movement, while the IPA results predicted the decreased activation in the cellular function of “cell spreading” (p = 0.00238, z-score = -2.121). The aggressiveness of transformed FLS may be affected by the severity of inflammatory arthritis [30,36]. The complexity of cell migration in synovial joint lies in the multifaceted biological factors dependent on cellular and ECM properties [38]. The location of synovial tissue in the affected OA joint may vary in histopathological changes [39,40], and the severity and degree of synovitis in OA may also influence the expression levels of cell adhesion and cellular movement related genes in FLS [41]. Of the 8 genes associated with function of cell spreading, including CDH11, LPP, SERPINE1, SMAD4, SPARC, TGFBI, CAP1, and MYH10, the dysregulated LPP is the candidate gene deducted from the differentially expressed mRNA results.

Through systematic bioinformatics analysis, we deducted the novel gene, LPP, 6.08-fold down-regulated in HFLS-OA potentially participated in altered cellular movement of OA synoviocytes. LPP (lipoma preferred partner) gene encodes LPP protein, a member of the zyxin family of LIM domain proteins that localizes at sites of cell adhesion and cell-cell contacts [42,43]. Besides, it also interacts with α-actinin (ACTN1) to participate in diverse cellular processes such as cell adhesion, spreading and migration [44]. LPP may partner with different molecules and reflect its various functions, therefore, LPP may exert as potential oncogene or oncosuppressor gene in different cancer cell lines [45,46]. The role of LPP in arthritic joint tissues has not been reported. However, ACTN1 was demonstrated to have increased expression in synovial tissues of RA compared to OA, and may participate in the signaling pathway initiated by TNF-α, thereby inducing cell proliferation and spreading of RA fibroblast-like synoviocytes to unaffected joints [47,48]. In our current study, OA FLS were obtained from commercialized human OA primary FLS. Whether the expression level of LPP can represent graded changes in different stages of OA and degree of synovitis and its interaction at sites of focal adhesion in OA FLS merits future investigation and validation in clinical synovial tissues of different OA stage subgroups.

The 2.37-fold up-regulated miR-141-3p in HLFS-OA was the miRNA potentially regulating LPP through systematic validation. The role of miR-141-3p in arthritis is still less investigated. One recent study by Park et al. reported the dysregulated miR-141-3p through ATP‐binding cassette transporter in human OA chondrocytes participated in altered lipid metabolism and suggested a novel insight into understanding the pathogenesis of OA [49]. In view of inflammatory conditions, miR-141 seemed to play dual roles in regulating inflammation-related signaling pathways. In diet-induced steatohepatitis mice, level of miR-141 in liver was significantly increased, and miR-141 knockout diminished hepatic inflammation and steatosis [50]. However, a recent report suggested that miR-141 decreased LPS-induced inflammation in human fibroblasts through p38 MAPK and NF-κB pathways [51]. As to the effect of miR-141-3p on cell function of fibroblasts, Feng et al. observed down-regulated miR-141-3p in keloid skin fibroblasts, and miR-141-3p exerted inhibitory effect on fibroblast proliferation and migration in keloids [52]. The results may be in line with our current findings of up-regulated miR-141-3p in HFLS-OA and functionally enriched biological functions of cellular movement in differentially expressed genes of HFLS-OA, particularly the predicted decreased activation in cell spreading from the IPA results. An interesting finding from literature review showed that miR-141 triggers senescence of human fibroblasts through the regulation of BMI1 expression [53]. In our HFLS-OA dataset, we also found the significantly down-regulated BMI1 in HFLS-OA, although it is not identified as one of the putative targets of miR-141-3p for its lower repression strength. Decreased autophagy in chondrocytes during aging leads to cellular senescence and ultimately OA progression [54]. Synoviocytes also serve to maintain joint homeostasis; however less is discussed about the senescence of synoviocytes and its link to OA progression. Further investigation on the role of miR-141-3p in OA FLS may have its clinical significance.

SKAP2 encodes src kinase associated phosphoprotein 2 and participates in different physiological processes, including integrin signaling, cell migration and cancer progression [55]. Alenghat et al. reported that SKAP2 is required for macrophage migration, chemotaxis and spreading, and may transmit signals required for cell migration [56], while others also claimed the importance of SKAP2 in integrin activation and neutrophil recruitment [57]. The evidence provides insight into further understanding of the essential role of SKAP2 in inflammatory disorders, including arthritis. Additionally, SKAP2 was detected in stromal cells infiltrating human gastric cancer tissues, and promoted tumor-associated macrophage infiltration and facilitated metastatic progression [58]. In contrary, Shimamura et al. demonstrated knockdown of SKAP2 in fibroblasts accelerated cell migration, and suggested the negative regulation of SKAP2 on cell migration of fibroblasts and tumor invasion of glioblastoma cells [59]. The evidence suggested that SKAP2 may exert distinct regulatory effect in different cell types. Our NGS results of HFLS showed the expression of SKAP2 in normal FLS was undetectable, while in OA FLS the expression of SKAP2 was 19.76 FPKM. The role of SKAP2 in OA FLS remains uncertain, but the expression pattern of SKAP2 in OA synovium may be associated with degree of inflammation in the OA joint microenvironment, and we proposed higher expression of SKAP2 may represent enhanced migratory effect of OA FLS that lead to synovial inflammation and predict cartilage breakdown in OA joint.

AP1S2 is a gene encoding sigma 2 subunit of the adaptor protein complex 1, and mediates the assembly of clathrin and the recognition of sorting signals of membrane proteins [60]. Studies related to AP1S2 have mainly focused on its mutation in neurodevelopmental disorders [61-63]. One recent report suggested miR-204, an effector of melanoma target therapy drug, exerts anti-migratory activity through targeting AP1S2 [64]. The potential miR-204-5p-AP1S2 interaction has been observed in our HFLS-OA data, with a 9.4-fold increase in AP1S2 expression and 3.22-fold decrease in miR-204-5p expression in HFLS-OA. To our understanding, there has been no related studies claiming the association of AP1S2 with inflammatory arthritis; however, expression of miR-204 were found decreased in T cells of RA patients [65], and decreased expression of miR-204 may enhance cartilage tissue destruction among OA patients through targeting IL-1β [66]. Whether miR-204-5p-AP1S2 interaction participates in OA synovitis or transformation of OA FLS merits future investigation.

PHF21A, also known as BHC80, is the molecule that recognizes unmethylated H3K4, and essential for neuronal and craniofacial development, and PHF21A loss of function may result in dysregulation of histone methylation and drive neurodevelopmental disorders [67,68]. In retinal endothelial cells, the expression of PHF21A was decreased under high glucose conditions, which was suggestive of the potential role in the development of diabetic retinopathy [69]. In mouse fibroblasts, PHF21A was one of the nuclear cofactors that modulate proinflammatory cytokine target genes, implicating its regulatory role in inflammatory status [70]. A 6.38-fold decrease in PHF21A in HFLS-OA was observed in our NGS result. PHF21A may be one of the potential inflammatory markers of OA, however, there is still lack of clinical or experimental evidence to link between PHF21A and OA synovitis in OA FLS.

Comparing the expression patterns of target genes with potential miRNA regulations identified from our HFLS-OA data in clinical samples of OA synovial tissues from GEO database, TFAP2A was the only gene significantly up-regulated in inflammatory OA synovium compared to that of non-inflammatory OA synovium. TFAP2A is a gene encoding activating enhancer-binding protein-2α (AP-2α) transcription factor, one of the members of AP-2 transcription factors known to play critical roles in cell cycle regulation and cell survival [71]. Previous studies suggested induction of transcription factor AP-2 expression during acute inflammation of rat primary afferent neurons [72] and in response to inflammatory cytokines in human keratinocytes [73], which are evident of possible AP-2 regulation through complex cytokine system and initiation of inflammatory process. Similar results were also observed in human lung fibroblasts and bronchial epithelial cells, in which induction of NF-κB and AP-2 led to IL-8 production and subsequent neutrophil infiltration and inflammation during bacterial infection [74]. Little is known about the role of AP-2 in arthritis. In OA chondrocytes, transfection of transcription factor AP-2ε led to increased expression of CXCL1, one of the chemokines inducing chondrocyte hypertrophy and apoptosis in OA [75]. Two recent studies also suggested TFAP2A is potentially one of the novel transcription regulators in psoriatic skin [76] and RA synoviocytes [77] through sequencing and bioinformatics approach. Taken together, we proposed the up-regulated TFAP2A in OA synovial tissue may represent higher degree of inflammation in OA, and TFAP2A may be a novel target for anti-inflammatory therapy in arthritis.

Conclusions

The current study identified differentially expressed genes in OA FLS were enriched in functions related to altered ECM, cell adhesion and cellular movement. In addition, miR-141-3p-LPP interaction potentially participated in cell migration and cell spreading of OA FLS, validated systematically in different miRNA target prediction databases. The current findings suggested novel insights into understanding the molecular signatures of FLS involved in the pathogenesis of OA, which may be potential therapeutic targets for the management of OA.

Acknowledgements

The authors thank the Center for Research Resources and Development of Kaohsiung Medical University.

Funding

This research was funded by Ministry of Science and Technology (MOST107-2320-B-037-011-MY3; MOST 107-2635-B-037-001; MOST 105-2314-B-037-012), Kaohsiung Medical University Hospital (KMUHS10701; KMUHS10712; KMUH107-7R82), and Kaohsiung Medical University (KMU-DK108003; KMU-TC108A04).

Competing Interests

The authors have declared that no competing interest exists.

References

1. Litwic A, Edwards MH, Dennison EM. et al. Epidemiology and burden of osteoarthritis. Br Med Bull. 2013;105:185-99

2. Wallace IJ, Worthington S, Felson DT. et al. Knee osteoarthritis has doubled in prevalence since the mid-20th century. Proc Natl Acad Sci. 2017;114:9332-6

3. Plotnikoff R, Karunamuni N, Lytvyak E. et al. Osteoarthritis prevalence and modifiable factors: a population study. BMC Public Health. 2015;15:1195

4. Martel-Pelletier J, Maheu E, Pelletier JP. et al. A new decision tree for diagnosis of osteoarthritis in primary care: international consensus of experts. Aging Clin Exp Res. 2019;31:19-30

5. Finan PH, Buenaver LF, Bounds SC. et al. Discordance between pain and radiographic severity in knee osteoarthritis: findings from quantitative sensory testing of central sensitization. Arthritis Rheum. 2013;65:363-72

6. Zhong L, Huang X, Karperien M. et al. The Regulatory Role of Signaling Crosstalk in Hypertrophy of MSCs and Human Articular Chondrocytes. Int J Mol Sci. 2015;16:19225-47

7. Mathiessen A, Conaghan PG. Synovitis in osteoarthritis: current understanding with therapeutic implications. Arthritis Res Ther. 2017;19:18

8. Attur M, Samuels J, Krasnokutsky S. et al. Targeting the synovial tissue for treating osteoarthritis (OA): where is the evidence? Best Pract Res Clin Rheumatol. 2010;24:71-9

9. Bhattaram P, Chandrasekharan U. The joint synovium: A critical determinant of articular cartilage fate in inflammatory joint diseases. Semin Cell Dev Biol. 2017;62:86-93

10. Silverstein AM, Stefani RM, Sobczak E. et al. Toward understanding the role of cartilage particulates in synovial inflammation. Osteoarthritis Cartilage. 2017;25:1353-61

11. Benito MJ, Veale DJ, FitzGerald O. et al. Synovial tissue inflammation in early and late osteoarthritis. Ann Rheum Dis. 2005;64:1263-7

12. Liu SC, Chuang SM, Hsu CJ. et al. CTGF increases vascular endothelial growth factor-dependent angiogenesis in human synovial fibroblasts by increasing miR-210 expression. Cell Death Dis. 2014;5:e1485

13. Kuo SJ, Yang WH, Liu SC. et al. Transforming growth factor beta1 enhances heme oxygenase 1 expression in human synovial fibroblasts by inhibiting microRNA 519b synthesis. PloS One. 2017;12:e0176052

14. Ko JY, Lee MS, Lian WS. et al. MicroRNA-29a Counteracts Synovitis in Knee Osteoarthritis Pathogenesis by Targeting VEGF. Sci Rep. 2017;7:3584

15. Wang JH, Shih KS, Wu YW. et al. Histone deacetylase inhibitors increase microRNA-146a expression and enhance negative regulation of interleukin-1beta signaling in osteoarthritis fibroblast-like synoviocytes. Osteoarthritis Cartilage. 2013;21:1987-96

16. Endisha H, Rockel J, Jurisica I. et al. The complex landscape of microRNAs in articular cartilage: biology, pathology, and therapeutic targets. JCI Insight. 2018:3

17. Tonge DP, Pearson MJ, Jones SW. The hallmarks of osteoarthritis and the potential to develop personalised disease-modifying pharmacological therapeutics. Osteoarthritis Cartilage. 2014;22:609-21

18. Bratus-Neuenschwander A, Castro-Giner F, Frank-Bertoncelj M. et al. Pain-Associated Transcriptome Changes in Synovium of Knee Osteoarthritis Patients. Genes. 2018:9

19. Gligorijevic V, Malod-Dognin N, Przulj N. Integrative methods for analyzing big data in precision medicine. Proteomics. 2016;16:741-58

20. Garcia-Campos MA, Espinal-Enriquez J, Hernandez-Lemus E. Pathway Analysis: State of the Art. Front Physiol. 2015;6:383

21. Kouskoumvekaki I, Shublaq N, Brunak S. Facilitating the use of large-scale biological data and tools in the era of translational bioinformatics. Brief Bioinform. 2014;15:942-52

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

23. Kramer A, Green J, Pollard J Jr. et al. Causal analysis approaches in Ingenuity Pathway Analysis. Bioinformatics. 2014;30:523-30

24. Szklarczyk D, Morris JH, Cook H. et al. The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 2017;45:D362-d8

25. Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 2003;4:2

26. Vejnar CE, Zdobnov EM. MiRmap: comprehensive prediction of microRNA target repression strength. Nucleic Acids Res. 2012;40:11673-83

27. Agarwal V, Bell GW, Nam JW. et al. Predicting effective microRNA target sites in mammalian mRNAs. Elife. 2015:4

28. Wong N, Wang X. miRDB: an online resource for microRNA target prediction and functional annotations. Nucleic Acids Res. 2015;43:D146-52

29. Clough E, Barrett T. The Gene Expression Omnibus Database. Methods Mol Biol. 2016;1418:93-110

30. Sellam J, Berenbaum F. The role of synovitis in pathophysiology and clinical symptoms of osteoarthritis. Nat Rev Rheumatol. 2010;6:625-35

31. Bottini N, Firestein GS. Duality of fibroblast-like synoviocytes in RA: passive responders and imprinted aggressors. Nat Rev Rheumatol. 2013;9:24-33

32. Lefevre S, Knedla A, Tennie C. et al. Synovial fibroblasts spread rheumatoid arthritis to unaffected joints. Nat Med. 2009;15:1414-20

33. Bhattaram P, Jones K. Regulation of fibroblast-like synoviocyte transformation by transcription factors in arthritic diseases. Biochem Pharmacol. 2019:10 1016/j.bcp.2019.03.018

34. Bhattaram P, Muschler G, Wixler V. et al. Inflammatory Cytokines Stabilize SOXC Transcription Factors to Mediate the Transformation of Fibroblast-Like Synoviocytes in Arthritic Disease. Arthritis Rheumatol. 2018;70:371-82

35. Steenvoorden MM, Huizinga TW, Verzijl N. et al. Activation of receptor for advanced glycation end products in osteoarthritis leads to increased stimulation of chondrocytes and synoviocytes. Arthritis Rheum. 2006;54:253-63

36. Ding X, Zhang Y, Huang Y. et al. Cadherin-11 involves in synovitis and increases the migratory and invasive capacity of fibroblast-like synoviocytes of osteoarthritis. Int Immunopharmacol. 2015;26:153-61

37. Chang X, Shen J, Yang H. et al. Upregulated expression of CCR3 in osteoarthritis and CCR3 mediated activation of fibroblast-like synoviocytes. Cytokine. 2016;77:211-9

38. Qu F, Guilak F, Mauck RL. Cell migration: implications for repair and regeneration in joint disease. Nat Rev Rheumatol. 2019;15:167-79

39. Shibakawa A, Aoki H, Masuko-Hongo K. et al. Presence of pannus-like tissue on osteoarthritic cartilage and its histological character. OsteoarthritisCartilage. 2003;11:133-40

40. Hutton CW, Hinton C, Dieppe PA. Intra-articular variation of synovial changes in knee arthritis: biopsy study comparing changes in patellofemoral synovium and the medial tibiofemoral synovium. Br J Rheumatol. 1987;26:5-8

41. Scanzello CR. Chemokines and inflammation in osteoarthritis: Insights from patients and animal models. J Orthop Res. 2017;35:735-9

42. Petit MM, Meulemans SM, Alen P. et al. The tumor suppressor Scrib interacts with the zyxin-related protein LPP, which shuttles between cell adhesion sites and the nucleus. BMC Cell Biol. 2005;6:1

43. Petit MM, Meulemans SM, Van de Ven WJ. The focal adhesion and nuclear targeting capacity of the LIM-containing lipoma-preferred partner (LPP) protein. J Biol Chem. 2003;278:2157-68

44. Li B, Zhuang L, Reinhard M. et al. The lipoma preferred partner LPP interacts with alpha-actinin. J Cell Sci. 2003;116:1359-66

45. Kuriyama S, Yoshida M, Yano S. et al. LPP inhibits collective cell migration during lung cancer dissemination. Oncogene. 2016;35:952-64

46. Colas E, Muinelo-Romay L, Alonso-Alconada L. et al. ETV5 cooperates with LPP as a sensor of extracellular signals and promotes EMT in endometrial carcinomas. Oncogene. 2012;31:4778-88

47. Huang RY, Huang QC, Burgering BM. Novel insight into the role of alpha-actinin-1 in rheumatoid arthritis. Discov Med. 2014;17:75-80

48. Qingchun H, Runyue H, LiGang J. et al. Comparison of the expression profile of apoptosis-associated genes in rheumatoid arthritis and osteoarthritis. Rheumatol Int. 2008;28:697-701

49. Park S, Oh J, Kim YI. et al. Suppression of ABCD2 dysregulates lipid metabolism via dysregulation of miR-141:ACSL4 in human osteoarthritis. Cell Biochem Funct. 2018;36:366-76

50. Tran M, Lee SM, Shin DJ. et al. Loss of miR-141/200c ameliorates hepatic steatosis and inflammation by reprogramming multiple signaling pathways in NASH. JCI Insight. 2017:2

51. Quan B, Zhang H, Xue R. miR-141 alleviates LPS-induced inflammation injury in WI-38 fibroblasts by up-regulation of NOX2. Life Sci. 2019;216:271-8

52. Feng J, Xue S, Pang Q. et al. miR-141-3p inhibits fibroblast proliferation and migration by targeting GAB1 in keloids. Biochem Biophys Res Commun. 2017;490:302-8

53. Dimri M, Carroll JD, Cho JH. et al. microRNA-141 regulates BMI1 expression and induces senescence in human diploid fibroblasts. Cell Cycle. 2013;12:3537-46

54. Vinatier C, Dominguez E, Guicheux J. et al. Role of the Inflammation-Autophagy-Senescence Integrative Network in Osteoarthritis. Front Physiol. 2018;9:706

55. Bureau JF, Cassonnet P, Grange L. et al. The SRC-family tyrosine kinase HCK shapes the landscape of SKAP2 interactome. Oncotarget. 2018;9:13102-15

56. Alenghat FJ, Baca QJ, Rubin NT. et al. Macrophages require Skap2 and Sirpalpha for integrin-stimulated cytoskeletal rearrangement. J Cell Sci. 2012;125:5535-45

57. Boras M, Volmering S, Bokemeyer A. et al. Skap2 is required for beta2 integrin-mediated neutrophil recruitment and functions. J Exp Med. 2017;214:851-74

58. Tanaka M, Shimamura S, Kuriyama S. et al. SKAP2 Promotes Podosome Formation to Facilitate Tumor-Associated Macrophage Infiltration and Metastatic Progression. Cancer Res. 2016;76:358-69

59. Shimamura S, Sasaki K, Tanaka M. The Src substrate SKAP2 regulates actin assembly by interacting with WAVE2 and cortactin proteins. J Biol Chem. 2013;288:1171-83

60. Kongsvik TL, Honing S, Bakke O. et al. Mechanism of interaction between leucine-based sorting signals from the invariant chain and clathrin-associated adaptor protein complexes AP1 and AP2. J Biol Chem. 2002;277:16484-8

61. Shaheen R, Sebai MA, Patel N. et al. The genetic landscape of familial congenital hydrocephalus. Ann Neurol. 2017;81:890-7

62. Cacciagli P, Desvignes JP, Girard N. et al. AP1S2 is mutated in X-linked Dandy-Walker malformation with intellectual disability, basal ganglia disease and seizures (Pettigrew syndrome). Eur J Hum Genet. 2014;22:363-8

63. Borck G, Molla-Herman A, Boddaert N. et al. Clinical, cellular, and neuropathological consequences of AP1S2 mutations: further delineation of a recognizable X-linked mental retardation syndrome. Hum Mutat. 2008;29:966-74

64. Vitiello M, Tuccoli A, D'Aurizio R. et al. Context-dependent miR-204 and miR-211 affect the biological properties of amelanotic and melanotic melanoma cells. Oncotarget. 2017;8:25395-417

65. Lai NS, Yu HC, Tung CH. et al. The role of aberrant expression of T cell miRNAs affected by TNF-alpha in the immunopathogenesis of rheumatoid arthritis. Arthritis Res Ther. 2017;19:261

66. Song X, Zhu M, Sun Y. et al. MiR-204 enhances the progression of osteoarthritis by suppressing the production of IL-1beta. Die Pharmazie. 2017;72:587-92

67. Garay PM, Wallner MA, Iwase S. Yin-yang actions of histone methylation regulatory complexes in the brain. Epigenomics. 2016;8:1689-708

68. Vallianatos CN, Iwase S. Disrupted intricacy of histone H3K4 methylation in neurodevelopmental disorders. Epigenomics. 2015;7:503-19

69. Ung C, Sanchez AV, Shen L. et al. Whole exome sequencing identification of novel candidate genes in patients with proliferative diabetic retinopathy. Vision Res. 2017;139:168-76

70. Meier-Soelch J, Jurida L, Weber A. et al. RNAi-Based Identification of Gene-Specific Nuclear Cofactor Networks Regulating Interleukin-1 Target Genes. Front Immunol. 2018;9:775

71. Eckert D, Buhl S, Weber S. et al. The AP-2 family of transcription factors. Genome Biol. 2005;6:246

72. Donaldson LF, McQueen DS, Seckl JR. Induction of transcription factor AP2 mRNA expression in rat primary afferent neurons during acute inflammation. Neurosci Lett. 1995;196:181-4

73. Oyama N, Iwatsuki K, Homma Y. et al. Induction of transcription factor AP-2 by inflammatory cytokines in human keratinocytes. J Invest Dermatol. 1999;113:600-6

74. Smith RS, Fedyk ER, Springer TA. et al. IL-8 production in human lung fibroblasts and epithelial cells activated by the Pseudomonas autoinducer N-3-oxododecanoyl homoserine lactone is transcriptionally regulated by NF-kappa B and activator protein-2. J Immunol. 2001;167:366-74

75. Wenke AK, Niebler S, Grassel S. et al. The transcription factor AP-2varepsilon regulates CXCL1 during cartilage development and in osteoarthritis. Osteoarthritis Cartilage. 2011;19:206-12

76. Zolotarenko A, Chekalin E, Mesentsev A. et al. Integrated computational approach to the analysis of RNA-seq data reveals new transcriptional regulators of psoriasis. Exp Mol Med. 2016;48:e268

77. Ham S, Bae JB, Lee S. et al. Epigenetic analysis in rheumatoid arthritis synoviocytes. Exp Mol Med. 2019;51:22

Author contact

Corresponding author: Chia-Hsin Chen; chchenedu.tw; Tel.: +886-7-312-1101 ext. 5962 Po-Lin Kuo; kuopolinnet.tw; Tel.: +886-7-312-1101 ext. 2512-33

Received 2019-4-9
Accepted 2019-7-5
Published 2019-7-21

Citation styles

©2024 Ivyspring International Publisher. Terms of use