Micro-ribonucleic acid profiles from microarray in ankylosing spondylitis.
Micro-ribonucleic acids (miRNAs) are noncoding RNA oligonucleotides that have been highly conserved during evolution and have recently emerged as potent regulators of gene expression of posttranscriptional regulators. (5) miRNAs have been identified and appear to be critical for fine-tuning many biological processes, in addition to offering the prospect of multiple targets. (6,7) Emerging data suggests that single miRNA species can profoundly alter the phenotype and outcome of immune responses. (8-10) Although dysregulated expression profiles of miRNAs have been identified in the blood of patients with immune disease, (11) only limited data is available regarding miRNA with regard to AS. (12-14)
Therefore, in this study, we aimed to detect candidate miRNAs from microarray within peripheral blood mononuclear cells (PBMCs) and synovial fluid mononuclear cells (SFMCs) in patients with AS.
PATIENTS AND METHODS
Samples from three AS patients (3 males, mean age 37.3 [+ or -] 2.5 years; range 35 to 40 years) with peripheral arthritis and three healthy controls (HCs) (3 males, mean age 39.0 [+ or -] 2.6 years; range 37 to 42 years) were obtained between December 2011 and January 2012 from Chonnam National University Hospital for miRNA microarray. To validate miRNA expression, samples from 25 HCs, and paired peripheral blood and synovial fluid from 21 patients with AS were obtained. The study was approved by the Institutional Review Board of Chonnam National University Hospital, and written informed consent was obtained from all participants.
Peripheral venous blood and synovial fluid samples were collected in heparin-containing tubes. PBMCs and SFMCs were isolated by density-gradient centrifugation using Ficoll-Paque Plus solution (Amersham Biosciences, Uppsala, Sweden). Freshly isolated PBMCs were suspended in a complete medium, consisting of RPMI 1640, 2 mM L-glutamine, 100 units/mL of penicillin, and 100 pg/mL of streptomycin, supplemented with 10% FBS (Gibco BRL, Grand Island, NY, USA).
Total RNA was isolated using the miRNeasy kit (Qiagen Inc., Valencia, CA, USA). The presence of miRNA in the total RNA was assayed using the Agilent 2100 Bioanalyzer Small RNA assay (Agilent Technologies, CA, USA). The microarray experiment proceeded only when the quality of total RNAs were considered as "passed" and their integrity was good by total RNA quality control using Agilent Bioanalyzer 2100. 100 ng of each RNA sample were hybridized to Agilent Human microRNA Microarray v2.0 (G4470B, Agilent Technologies). MicroRNA labeling, hybridization, and washing were carried out following Agilent's instructions. Images of hybridized microarrays were acquired with a microarray scanner (Agilent G2565BA), and features were extracted using the agilent feature extraction image analysis tool version A.18.104.22.168. To understand the impact of the storage condition on the miRNA expression profiles, hierarchical clustering was performed. Microarray data were quantile-normalized and analyzed using GeneSpring software (Agilent).
The miScript Reverse Transcription Kit (Qiagen GmbH, Hilden, Germany) was used for complementary deoxyribonucleic acid (cDNA) preparation. MiScript primer assays (Qiagen GmbH, Hilden, Germany) were used for semiquantitative determination of the expression of human miRNAs. The expression of U6B snRNA was used as an endogenous control. A threshold cycle (CT) was observed in the exponential phase of amplification, and quantification of relative expression levels were performed using standard curves for target genes and the endogenous control. Geometric means were used to calculate the [DELTA][DELTA]CT (delta-delta CT) values and were expressed as 2-[DELTA][DELTA]CT. The value of each control sample was set at one and used to calculate the fold-change of the target genes.
The Kruskal-Wallis test was used for the comparison of gene expression among groups. Dunn's test was used post-hoc to compare all pairs of groups. A p-value <0.05 was considered statistically significant. All statistical analyses were performed using PASW version 17.0 software (SPSS Inc., Chicago, IL, USA).
A total of 887 miRNAs were screened by microarray among groups (Figure 1). After normalization of the raw data, we noted that the expression of five miRNAs, including hsa-miR-136, hsa-miR-29c*, hsa-miR-376a, hsamiR-376c, hsa-miR-377, were significantly lower (fold change [less than or equal to]0.5 and p[less than or equal to]0.05) and only hsa-miR424-5p was significantly higher in AS PBMCs (fold change [greater than or equal to]2 and p<0.05, Table 1). In AS SFMCs, it was identified that expressions of 16 miRNAs were significantly down regulated: hsa-miR-125b, hsa-miR-126, hsa-miR-126*, hsa-miR-136, hsa-miR-144, hsa-miR-151-3p, hsa-miR-199a-5p, hsa-miR-29c*, hsa-miR-30b, hsa-miR-326, hsa-miR-376a, hsa-miR-376c, hsa-miR-377, hsa-miR-451, hsa-miR-486-5p, hsa-miR-766 (fold change [less than or equal to]0.5 and p[less than or equal to]0.05) whereas only hsa-miR-424-5p was significantly upregulated (fold change [greater than or equal to]2 and p[less than or equal to]0.05, Table 1). All above-mentioned miRNAs were re-evaluated for further validation. Finally, a significantly elevated hsa-miR-424-5p in SFMCs from AS patients (8.95 [+ or -] 0.25) compared with HCs (1.31 [+ or -] 0.08, p<0.01) (Figure 2a) was notified. As shown in Figure 2b, hsa-miR-377 was significantly decreased in SFMCs from AS patients (0.01 [+ or -] 0.00) compared with HCs (0.97 [+ or -] 0.04, p<0.01). There were no significant differences of gene expression in other miRNAs among groups (data not shown). Based on target prediction programs (miRTarBase, TargetScan, miRanda, and MiRBase) and published papers, we found highly evidenced candidate target genes (hypoxia-inducible factor 1A and Cullin 2 genes for hypoxia and angiogenesis; cyclin-dependent kinase 6 (CDK6), Cyclin D1 (CCND1), Cyclin E1 (CCNE1), and CCND3 genes for cell cycle) of hsa-miR-424-5p. Meanwhile, protein phosphatase magnesium-dependent 1A (PPM1A) and p21-activated kinase 1 (PAK1) were predicted as candidate genes of hsa-miR-377 (Table 2).
[FIGURE 1 OMITTED]
In the present study, miRNAs in the PBMCs and SFMCs of AS patients were comprehensively isolated and analyzed by using miRNA microarray analysis. Through this analysis, we discovered significantly high expression of miR-424-5p and lower expression of miR-377 in SFMCs from patients with AS, compared to HCs.
Based on target prediction programs and published papers, miR-424-5p was highly linked to pathways of hypoxia and the cell cycle. Inflammatory arthritis is a chronic progressive autoimmune disease characterized by the proliferation of synovial tissue. The role of hypoxia in the pathogenesis of inflammatory arthritis has been demonstrated in several in vitro and in vivo studies. (15-17) Hypoxia is a key driving force for angiogenesis and is recognized as an important event in inflammatory arthritis. Several studies have demonstrated hypoxia-inducible factor-1[alpha] expression in synovial tissue. It has been revealed in synovial cell cultures, that hypoxia induces key angiogenic growth factors, chemokines, and matrix metalloproteinases. (18,19) Various stresses, including hypoxia, regulate the expression and function of miRNA. (20) For example, a subset of miRNAs induced by hypoxia contributes to the regulation of the broad spectrum of genes regulated by hypoxia. It was reported that miR-424-5p is a critical mediator of oxygen-dependent changes in endothelial cells and is physiologically upregulated in tissues undergoing vascular remodelin g and angiogenesis. (21) Cyclin D1, the product of CCND1 gene, is an important regulator of G1 to S phase progression in many different cell types. Cyclin D1, together with its binding partners CDK4 and CDK6, forms active complexes that promote cell cycle progression by phosphorylating and inactivating the retinoblastoma protein. (22) According to our data, miR-424-5p might inhibit cell proliferation, arrest cell cycle progression, and induce cell apoptosis in AS.
[FIGURE 2 OMITTED]
Recently, it was published that PPM1A was highly expressed in AS synovial tissue. Moreover, PPM1A detecting antibodies were correlated positively with Bath Ankylosing Spondylitis Disease Activity Index score. (23) Given the result from our data, decreased expression of miR-377 in SFMCs might contribute to increase levels of PPM1A in AS synovium. The expression of PAK1 was increased in rheumatoid arthritis fibroblast-like synoviocytes fibroblast-like synoviocytes. PAK1 plays an important role in regulating the migration, invasion and production and activity of matrix metalloproteinases in fibroblast-like synoviocytes, which is mediated by the c-Jun N-terminal kinase pathway. (24) This might suggest that hsa-miR-377 has a role on PAK1 expression in arthritis of AS, although there is no report of PAK1 in AS.
In conclusion, miR-424-5p was increased and miR-377 was decreased in SFMCs from patients with AS. These two miRs might have functional roles in patients with arthritis via different pathways. However, further research is necessary to ascertain whether miR-424-5p and miR-377 can truly have functional roles in AS patients with peripheral arthritis. Due to the small sample size, further validation of other candidate miRNAs from this microarray assay should be carried out. Nevertheless, data from this microarray may provide important information for research on miRNA in AS.
Declaration of conflicting interests
The authors declared no conflicts of interest with respect to the authorship and/or publication of this article.
This study was supported by grants from the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science, and Technology (grant no. 2011-0008867), by the Chonnam National University Hospital Research Institute of Clinical Medicine (grant no. CRI 11074-21).
Received: August 12, 2015 Accepted: October 27, 2015 Published online: December 12, 2015
Correspondence: Tae-Jong Kim, MD. Department of Rheumatology, Research Institute of Medical Sciences, Chonnam National University Medical School and Hospital, 501-757 Gwangju, South Korea. Tel: 07468079257 e-mail: email@example.com
(1.) Khan MA, van der Linden SM. A wider spectrum of spondyloarthropathies. Semin Arthritis Rheum 1990;20:107-13.
(2.) Kim TJ, Lee S, Joo KB, Park DJ, Park YW, Lee SS, et al. The presence of peripheral arthritis delays spinal radiographic progression in ankylosing spondylitis: Observation Study of the Korean Spondyloarthropathy Registry. Rheumatology (Oxford) 2014;53:1404-8.
(3.) Kim TJ, Kim TH. Clinical spectrum of ankylosing spondylitis in Korea. Joint Bone Spine 2010;77:235-40.
(4.) Heuft-Dorenbosch L, van Tubergen A, Spoorenberg A, Landewe R, Dougados M, Mielants H, et al. The influence of peripheral arthritis on disease activity in ankylosing spondylitis patients as measured with the Bath Ankylosing Spondylitis Disease Activity Index. Arthritis Rheum 2004;51:154-9.
(5.) Bartel DR MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 2004;116:281-97.
(6.) O'Connell RM, Rao DS, Chaudhuri AA, Baltimore D. Physiological and pathological roles for microRNAs in the immune system. Nat Rev Immunol 2010;10:111-22.
(7.) Pillai RS, Bhattacharyya SN, Filipowicz W. Repression of protein synthesis by miRNAs: how many mechanisms? Trends Cell Biol 2007;17:118-26.
(8.) Jin HM, Kim TJ, Choi JH, Kim MJ, Cho YN, Nam KI, et al. MicroRNA-155 as a proinflammatory regulator via SHIP-1 down-regulation in acute gouty arthritis. Arthritis Res Ther 2014;16:88.
(9.) Rodriguez A, Vigorito E, Clare S, Warren MV, Couttet P, Soond DR, et al. Rodriguez A, Vigorito E, Clare S, Warren MV, Couttet P, Soond DR, et al. Science 2007;316:608-11.
(10.) Xiao C, Srinivasan L, Calado DP, Patterson HC, Zhang B, Wang J, et al. Lymphoproliferative disease and autoimmunity in mice with increased miR17-92 expression in lymphocytes. Nat Immunol 2008;9:405-14.
(11.) Murata K, Yoshitomi H, Tanida S, Ishikawa M, Nishitani K, Ito H, et al. Plasma and synovial fluid microRNAs as potential biomarkers of rheumatoid arthritis and osteoarthritis. Arthritis Res Ther 2010;12:86.
(12.) Huang J, Song G, Yin Z, Luo X, Ye Z. Elevated miR-29a expression is not correlated with disease activity index in PBMCs of patients with ankylosing spondylitis. Mod Rheumatol 2014;24:331-4.
(13.) Lai NS, Yu HC, Chen HC, Yu CL, Huang HB, Lu MC. Aberrant expression of microRNAs in T cells from patients with ankylosing spondylitis contributes to the immunopathogenesis. Clin Exp Immunol 2013;173:47-57.
(14.) Huang CH, Wei JC, Chang WC, Chiou SY, Chou CH, Lin YJ, et al. Higher expression of whole blood microRNA-21 in patients with ankylosing spondylitis associated with programmed cell death 4 mRNA expression and collagen cross-linked C-telopeptide concentration. J Rheumatol 2014;41:1104-11.
(15.) Konisti S, Kiriakidis S, Paleolog EM. Hypoxia--a key regulator of angiogenesis and inflammation in rheumatoid arthritis. Nat Rev Rheumatol 2012;8:153-62.
(16.) Del Rey MJ, Izquierdo E, Usategui A, Gonzalo E, Blanco FJ, Acquadro F, et al. The transcriptional response of normal and rheumatoid arthritis synovial fibroblasts to hypoxia. Arthritis Rheum 2010;62:3584-94.
(17.) Brouwer E, Gouw AS, Posthumus MD, van Leeuwen MA, Boerboom AL, Bijzet J, et al. Hypoxia inducible factor-1-alpha (HIF-1alpha) is related to both angiogenesis and inflammation in rheumatoid arthritis. Clin Exp Rheumatol 2009;27:945-51.
(18.) Akhavani MA, Madden L, Buysschaert I, Sivakumar B, Kang N, Paleolog EM. Hypoxia upregulates angiogenesis and synovial cell migration in rheumatoid arthritis. Arthritis Res Ther 2009;11:64.
(19.) Muz B, Khan MN, Kiriakidis S, Paleolog EM. Hypoxia. The role of hypoxia and HIF-dependent signalling events in rheumatoid arthritis. Arthritis Res Ther 2009;11:201.
(20.) Kulshreshtha R, Ferracin M, Wojcik SE, Garzon R, Alder H, Agosto-Perez FJ, et al. A microRNA signature of hypoxia. Mol Cell Biol 2007;27:1859-67.
(21.) Ghosh G, Subramanian IV, Adhikari N, Zhang X, Joshi HP, Basi D, et al. Hypoxia-induced microRNA-424 expression in human endothelial cells regulates HIF-a isoforms and promotes angiogenesis. J Clin Invest 2010;120:4141-54.
(22.) Alao JP. The regulation of cyclin D1 degradation: roles in cancer development and the potential for therapeutic invention. Mol Cancer 2007;6:24.
(23.) Kim YG, Sohn DH, Zhao X, Sokolove J, Lindstrom TM, Yoo B, et al. Role of protein phosphatase magnesium-dependent 1A and antiprotein phosphatase magnesium-dependent 1A autoantibodies in ankylosing spondylitis. Arthritis Rheumatol 2014;66:2793-803.
(24.) Fu D, Yang Y, Xiao Y, Lin H, Ye Y, Zhan Z, et al. Role of p21-activated kinase 1 in regulating the migration and invasion of fibroblast-like synoviocytes from rheumatoid arthritis patients. Rheumatology (Oxford) 2012;51:1170-80.
Hye-Mi JIN,  Young-Nan CHO,  Seung-Jung KEE,  Shin-Seok LEE,  Yong-Wook PARK,  Tae-Jong KIM 
 Department of Rheumatology, Research Institute of Medical Sciences, Chonnam National University Medical School and Hospital, Gwangju, South Korea
 Department of Laboratory Medicine, Research Institute of Medical Sciences, Chonnam National University Medical School and Hospital, Gwangju, South Korea
Table 1. Comparison of micro-ribonucleic acid level of peripheral blood mononuclear cells and synovial fluid mononuclear cells in patients with ankylosing spondylitis PBMC (Fold change) Name of gene Downregulated Upregulated hsa-miR-125b 0.523497 hsa-miR-126 0.989866 hsa-miR-126 * 0.849571 hsa-miR-136 0.099969 hsa-miR-144 1.833485 hsa-miR-151-3p 1.311719 hsa-miR-199a-5p 0.994883 hsa-miR-29c * 0.219023 hsa-miR-30b 0.783106 hsa-miR-326 0.543503 hsa-miR-376a 0.184359 hsa-miR-376c 0.250179 hsa-miR-377 0.219750 hsa-miR-424-5p 2.158072 hsa-miR-451 1.921773 hsa-miR-486-5p 1.091827 hsa-miR-766 0.698829 SFMC (Fold change) Name of gene Downregulated Upregulated hsa-miR-125b 0.01 hsa-miR-126 0.01 hsa-miR-126 * 0.01 hsa-miR-136 0.01 hsa-miR-144 0.01 hsa-miR-151-3p 0.01 hsa-miR-199a-5p 0.01 hsa-miR-29c * 0.01 hsa-miR-30b 0.445761 hsa-miR-326 0.01 hsa-miR-376a 0.01 hsa-miR-376c 0.01 hsa-miR-377 0.01 hsa-miR-424-5p 3.522999 hsa-miR-451 0.051164 hsa-miR-486-5p 0.01 hsa-miR-766 0.01 PBMC: Peripheral blood mononuclear cell; SFMC: Synovial fluid mononuclear cell; Fold [greater than or equal to] 2 or [less than or equal to] 0.5 indicates the difference is significant; hsa-miR-126 * : It is the original gene name of mature sequence hsa-mir-126-5p; hsa-miR-29c*: It is the original gene name of mature sequence hsa-mir-29c-5p. Table 2. Candidate target genes of hsa-miR-424-5p and hsa-miR-377 Gene name Main pathway Target Accession gene ID hsa-miR-424-5p Hypoxia and angiogenesis HIF1A MIRT005926 CUL2 MIRT005927 Cell cycle CDK6 MIRT000938 CCND1 MIRT000941 CCNE1 MIRT000936 CCND3 MIRT000937 hsa-miR-377 MAPK pathway PPM1A MIRT000990 PI3K-akt pathway PAK1 MIRT000991 HIF1A: Hypoxia-inducible factor 1A; CUL2: Cullin 2; CDK: Cyclin-dependent kinase; CCND1: Cyclin D1; CCNE1: Cyclin E1; CCND3: Cyclin D3; MAPK: Mitogen-activated protein kinase; PPM1A: Protein phosphatase magnesium-dependent 1A; PI3K: Phosphatidyl inositol 3'-kinase; PAK1: p21-activated kinase 1.
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||ORIGINAL ARTICLE|
|Author:||Jin, Hye-Mi; Cho, Young-Nan; Kee, Seung-Jung; Lee, Shin-Seok; Park, Yong-Wook; Kim, Tae-Jong|
|Publication:||Turkish Journal of Rheumatology|
|Date:||Jun 1, 2016|
|Previous Article:||Glutathione-S-transferase variants are not associated with increased carotid intima-media thickness in Turkish Familial Mediterranean fever patients.|
|Next Article:||Are we able to suppress disease activity adequately in patients with established rheumatoid arthritis? An observational and cross-sectional study.|