- Research article
- Open Access
Circulating plasma microRNA profiling in patients with polymyositis/dermatomyositis before and after treatment: miRNA may be associated with polymyositis/dermatomyositis
Inflammation and Regenerationvolume 38, Article number: 1 (2018)
MicroRNAs (miRNAs) are involved in the regulation of key biological processes and have been implicated in various diseases, including autoimmune disorders. The pathogenesis of polymyositis (PM) and dermatomyositis (DM) is considered to be mediated by autoimmune reactions. To determine miRNA role in the development and progression of PM and DM, we performed plasma miRNA profiling in PM/DM patients before and after treatment.
Total RNA was isolated from plasma of 10 patients before and after treatment with prednisolone, or, in case of prednisolone resistance or complications, with the combination of calcineurin inhibitors (cyclosporine or tacrolims) and/or pulse intravenous cyclophosphamide. The expression of miRNAs was determined using miRNA microarray and validated by qRT-PCR.
More differentially expressed miRNAs were found in plasma of DM patients compared to PM patients before and after treatment, and their profiles were different. Among the differentially expressed plasma miRNA identified by microarray, the levels of hsa-miR-4442 were confirmed by qRT-PCR to be significantly decreased by treatment. In addition, plasma hsa-miR-4442 content in active PM/DM significantly exceeded that in other active autoimmune diseases such as rheumatoid arthritis and systemic lupus erythematosus, as well as in healthy individuals. The level of plasma hsa-miR-4442 was positively correlated with Skeletal Disease Activity in MITAX (Myositis Intention to Treat Activity Index).
This is the first report describing plasma miRNA expression profiles in PM/DM patients. The present data suggest that plasma levels of miRNAs may be associated with polymyositis/dermatomyositis and hsa-miR-4442 could be used as a biomarker for PM/DM diagnosis and/or disease activity.
Polymyositis (PM) is an idiopathic inflammatory disease characterized by the degeneration of the muscles, while dermatomyositis (DM) is polymyositis accompanied by skin inflammation manifested by heliotrope eyelids and Gottron’s signs which are red papules erupting on the finger joints. The pathogenesis of PM/DM is considered to be mediated by autoimmune response triggered by various causes, including cancer, and the disease can overlap with other autoimmune disorders such as systemic lupus erythematosus (SLE), rheumatoid arthritis (RA), and systemic sclerosis. However, the exact cause of PM/DM is still unclear.
Micro (mi)RNAs are small noncoding RNAs about 22 nucleotides in length that usually function as negative regulators of target mRNA translation by binding to mRNA 3′-untranslated region (UTR) . miRNAs are involved in the regulation of key biological processes, including immune response and cell differentiation, proliferation, and apoptosis. The pathogenic role of miRNAs has been intensely studied in malignant diseases as well as in autoimmune disorders [2, 3]. miRNAs have been found in plasma and other body fluids of humans, indicating their high stability in the extracellular environment and suggesting that they may control cell-to-cell communication in health and disease. Highly stable extracellular miRNAs circulating in blood can be delivered by various cargo molecules such as high-density lipoproteins, nucleophosmin, and Argonaute 2 protein, to recipient cells where they regulate the expression of key proteins involved in disease pathogenesis [4,5,6,7]. Therefore, extracellular miRNAs present in blood circulation have attracted attention as disease biomarkers . Furthermore, miRNA profiling of body fluids, including plasma and serum, has revealed significant differences in the spectrum and concentration of miRNAs between blood and tissues , suggesting that blood miRNA content may reflect specific disease-associated processes.
There are some reports about the changes of miRNA level in skin, muscle, or serum of PM/DM patients , but miRNA concentration in plasma has not been investigated. Based on recent reports, we also expected that more abundant microRNAs will be collected from exosomes and lipoproteins contained in plasma . In addition, our previous high throughput studies of patients with autoimmune conditions before and after treatment revealed that serum/plasma contained biomarkers of disease pathogenesis and progression [12,13,14], suggesting that plasma may also be a source of PM/DM markers, including miRNA.
To test this hypothesis, here, we examined the changes in plasma miRNA profiles of patients with PM/DM before and after treatment.
Patients and samples
Plasma samples were obtained from four patients with PM and four patients with DM before and after treatment (active and inactive phase, respectively) to use in microarray. In order to validate the expression by qRT-PCR, we added samples from one PM patient and one DM patient. Clinicopathological characteristics of the patients enrolled in this study are shown in Table 1. All the patients fulfilled the criteria of Bohan and Peter, which include symmetrical muscle weakness, myositis, increase in serum skeletal muscle enzymes, characteristic electromyogram, and typical rash [15, 16]. The patients were analyzed for malignant tumors by computerized tomography and gastrointestinal endoscopy, and they did not have malignant tumors. All the patients received therapy with prednisolone were treated by the combination of calcineurin inhibitors (cyclosporine or tacrolims) and/or intravenous cyclophosphamide pulse therapy. Post treatment samples were obtained from patients who improved muscle symptoms. We investigated the Myositis Intention to Treat Activity Index (MITAX), CK, ALD, LDH, ESR as disease activity evaluation. MITAX is modified from the BILAG approach to assess disease activity in lupus . MITAX assesses specific manifestations in seven organ systems (constitutional, cutaneous, skeletal, gastrointestinal, pulmonary, cardiac, and muscle systems). For MITAX assessment, three to nine criteria related to symptoms, physical findings, or laboratory abnormalities in each of the seven organ systems are assessed and converted into a score that range from zero to nine. The total disease activity score is calculated by summing the score of each organ system. Healthy individuals (seven women and three men, age 24–48 years old), RA patients (10 women, 51–82 years old), and SLE patients (11 women and 1 man, 17–54 years old) served as controls. RA patients were treated with conventional synthetic disease-modifying anti-rheumatic drugs (csDMARDs) or biologics in addition to csDMARDs, and the mean of DAS-ESR was 5.26. They were scheduled to introduce or to switch the biologics because they were ineffective. SLE patients were the onset or the recurrence, and the mean of SLEDAI scores was 17.4.
miRNA array analysis
Total RNA was extracted from 300 μl of plasma obtained from four PM and four DM patients using the 3D-Gene™ RNA extraction kit (Toray Industries Inc., Tokyo, Japan) and miRNeasy Mini Kit (Qiagen, Venlo, The Netherlands). The extracted RNA was labeled using the 3D-Gene miRNA labeling kit (Toray Industries Inc.), and the labeled targets were then hybridized to a 3D-Gene Human miRNA 4-plex chip (V21_V1.0.0, Toray Industries Inc.). Hybridization images were scanned in a GenePix4400A device (Molecular Devices, Sunnyvale, CA, USA), and miRNA expression was assessed based on signal intensity calculated as the median of foreground signals minus the mean of negative control signals + 2SD. As there was no definite internal control for plasma miRNA, the obtained data were normalized to the median intensity levels of miR-762, miR-3665, miR-3960, and miR-4516 using a per-chip 95th percentile method.
Quantitative reverse transcription polymerase chain reaction
Total RNA was extracted from 200 μl of individual plasma samples using ISOGEN II (Nippon Gene Co., Ltd., Toyama, Japan) and the miRNeasy Mini Kit (Qiagen). For normalization of sample-to-sample variation in the RNA isolation step, 25 fmol of Caenorhabditis elegans cel-miR-39-3p (mirVana miRNA mimic; Thermo Fisher Scientific, Waltham, MA, USA) was added to each denatured sample mixed with lysis buffer. Reverse transcription was carried out using the TaqMan MicroRNA Reverse Transcription Kit (Thermo Fisher Scientific) and specific primers (hsa-miR-4442, Assay ID: 463327_mat; hsa-miR-3187-3p, Assay ID: 245619_mat; cel-miR-39-3p, Assay ID: 000200; Thermo Fisher Scientific), and qRT-PCR was performed using THUNDERBIRD Probe qPCR Mix (TOYOBO Co., Ltd., Osaka, Japan) in an ABI Prism 7900HT system (Thermo Fisher Scientific) according to the manufacturer’s instructions. Cycle threshold (Ct) values were calculated using the SDS 1.4 software (Thermo Fisher Scientific), and miRNA expression was normalized to cel-miR-39-3p using the ∆∆Ct method.
Microarray data were analyzed using Microsoft Excel (Microsoft, Redmond, WA, USA). Welch’s t test was applied to compare the numbers of miRNAs, and paired t test was used to compare samples taken from patients before and after treatment. Statistical analysis of miRNA expression was performed with R (version 3.2.3; available from https://www.r-project.org/). The Wilcoxon signed-rank test was used to identify differentially expressed miRNAs in paired samples before and after treatment, and Steel’s test was used to compare the expression of miRNA in PM/DM patients with that in RA and SLE patients or healthy individuals. The Spearman rank correlation was used to analyze the correlation of the level of miRNA and the disease activities. P values of less than 0.05 were considered statistically significant.
miRNA levels in plasma of PM/DM patients before and after treatment
Analysis of circulating extracellular miRNA in patients with myositis revealed differences in miRNA plasma levels between DM and PM patients and showed that treatment significantly affected miRNA expression in plasma (Fig. 1). The numbers of differentially expressed miRNAs with at least twofold change caused by treatment indicated that there were significantly more miRNAs whose plasma levels were affected by treatment in DM compared to PM (P = 0.0093) (Table 2). The average numbers of miRNAs differentially expressed in plasma of DM and PM patients before and after treatment are presented as the Venn diagram (Fig. 2), which showed the difference between DM and PM regarding the treatment effect. There were 11 upregulated and 4 downregulated miRNAs common for DM and PM (Table 3); among them, the expression of hsa-miR-3187-3p was significantly upregulated and that of hsa-miR-4442 was significantly downregulated (P = 0.023 and P = 0.016, respectively).
Validation of hsa-miR-3187-3p and hsa-miR-4442 expression in plasma
To validate the differential expression of hsa-miR-3187-3p and hsa-miR-4442, we performed qRT-PCR which confirmed that the levels of plasma hsa-miR-4442 were markedly decreased by treatment (P = 0.047). However, there was no significant difference in the expression of hsa-miR-3187-3p (P = 0.28) (Fig. 3). These results indicate that hsa-miR-4442 may be associated with PM/DM.
hsa-miR-4442 expression in plasma from patients with active PM/DM, RA, and SLE, and healthy individuals
To verify the relationship between plasma hsa-miR-4442 levels and PM/DM, we compared the expression of this miRNA in PM/DM patients with that in patients with other autoimmune diseases (RA and SLE) and healthy individuals. Plasma expression of hsa-miR-4442 was significantly higher in active PM/DM compared with active RA (P = 0.00045), active SLE (P = 0.0011), or healthy control (P = 0.028) (Fig. 4a), further supporting the notion that circulating extracellular hsa-miR-4442 may be associated with PM/DM. We additionally analyzed the differential expression of hsa-miR-4442 in active/inactive SLE and RA patients using qRT-PCR (Fig. 4b, c). There was no statistically significant difference in either patient.
Correlation of hsa-miR-4442 level in plasma and disease activities
We examined the correlation between the plasma hsa-miR-4442 level and disease activities. The level of miRNA was positively correlated with Skeletal Disease Activity in MITAX (Table 4). There was a weak positive correlation with ESR. There was no correlation with MITAX total score, CK, ALD, and LDH.
Serum or plasma miRNAs have been reported to be associated with RA , SLE , and other autoimmune diseases . The current study was the first to compare plasma miRNA expression in PM/DM patients before and after treatment using miRNA microarray. We found a difference in plasma miRNA profiles between PM and DM patients and identified several miRNAs whose plasma levels were affected by treatment (Figs. 1 and 2). A recent study revealed differential expression of serum miRNAs between PM and DM, showing significant association of three miRNAs (hsa-miR-3676, hsa-miR-3907, and hsa-miR-877*) with disease activity . However, we did not observe changes in these miRNAs in our study, suggesting that miRNA spectra in serum and plasma of PM/DM patients may be different.
There were some reports about the miRNA level in skin and muscle in inflammatory myopathies, and we compare our plasma data with these reports. Inoue et al. showed that skin tissue miRNA array analysis demonstrated that the hsa-miR-223 level was markedly decreased in Gottron’s papules of DM and CADM . The transfection of a specific inhibitor of hsa-miR-223 in keratinocytes led to upregulation of the PKCε protein and resulted in abnormally increased cell proliferation. The serum hsa-miR-223 concentration was decreased in PM/DM patients, particularly in CADM (clinically amyopathic dermatomyositis) patients, compared with healthy controls. The expression of plasma hsa-miR-223 in our study seemed to be increased after treatment. However, the increase was not significant. This result could support that low level of hsa-miR-223 was associated with skin lesion and increased in parallel with treatment. hsa-miR-1, hsa-miR-133a, hsa-miR-133b, and hsa-miR-206 are critical regulators of myoblast-to-myocyte differentiation through regulation of multiple genes . Georgantas et al. observed decreased expression of hsa-miR-1, hsa-miR-133a, and hsa-miR-133b in PM/DM and inclusion body myositis (IBM), as well as decreased expression of miR-206 in DM . TNFα was significantly inversely correlated with decreased myogenic miRNA expression in the inflammatory myopathy subtypes. These miRNAs could not be detected in our miRNA array analysis. These miRNAs in muscle tissue may not be released into the blood.
Recently, the classification of PM/DM was being changed by the concept of viewing inflammatory myopathies as a spectrum, with muscle and skin involvement occurring to varying degrees . ADM (amyopathic dermatomyositis) is on the pure skin portion of the spectrum, while HDM (hypomyopathic dermatomyositis) is slightly closer to the muscle side. CDM (classic dermatomyositis) is in the middle of the spectrum and PM have purely muscle disease with a different pathologic process from ADM, HDM and CDM. Furthermore, the subsets of autoantibody status (such as anti-synthetase antibodies, Mi-2, SRP, and others) were increasingly being recognized as an improved way to phenotype patients with inflammatory myopathies . It is possible that the difference of miRNA profiles between PM and DM is involved in the appearance of these spectrum, phenotypes, and autoantibody status. In this study, since the number of available samples was not sufficient, only miRNAs commonly fluctuated in PM/DM were analyzed. We thought that these miRNAs might be involved in common pathology to PM/DM.
Our results indicate that the expression of plasma hsa-miR-4442 was significantly decreased after treatment (Fig. 3) and that it was significantly higher in PM/DM compared to RA, SLE, or healthy controls (Fig. 4a). The change of hsa-miR-4442 could not be denied the possibility of merely the influence of the drug; however, the fluctuations of SLE and RA seemed to be not as much as PM/DM (Fig. 4b, c). Therefore, we believe that the change may not be caused just by the drug. There has been no report on the function of hsa-miR-4442 or its association with a disease in general. Furthermore, we confirmed a positive correlation between plasma has-miR-4442 level and Skeletal Disease Activity in MITAX and a weak positive correlation with ESR (Table 4). However, we could not find any correlation between hsa-miR-4442 expression level and serum autoantibody profiles (Additional file 1: Table S1). Thus, this is the first study to reveal possible correlation of the plasma hsa-miR-4442 level with a pathological condition (myositis), suggesting that it may be a candidate biomarker for PM/DM diagnosis or disease activity.
However, this study had limitations. First, the sample size was small, which might have caused a type 2 error in statistical analysis. Therefore, the role of hsa-miR-4442 in PM/DM suggested here should be confirmed in studies using larger patient cohorts. Also, we did not investigate the cellular source of plasma hsa-miR-4442 or its functional role in the pathogenesis of PM/DM, as well as the mechanism underlying the decrease of plasma hsa-miR-4442 content after treatment. These important aspects should be investigated to clarify the significance of hsa-miR-4442 in the development and progression of autoimmune myositis.
In summary, we found differences in plasma miRNA profiles between PM and DM and between active and inactive PM/DM. Among the differentially expressed plasma miRNAs, has-miR-4442 levels were confirmed to be significantly lower after treatment, while being significantly higher in active PM/DM compared with RA, SLE, or healthy individuals. This is the first report about plasma miRNA profiling in PM/DM patients, laying a foundation for further studies to reveal the role of miRNAs in PM/DM.
Clinically amyopathic dermatomyositis
Myositis Intention to Treat Activity Index
Polymyositis and dermatomyositis
Quantitative reverse transcription polymerase chain reaction
Systemic lupus erythematosus
Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–97.
Jakymiw A, Ikeda K, Fritzler MJ, Reeves WH, Satoh M, Chan EK. Autoimmune targeting of key components of RNA interference. Arthritis Res Ther. 2006;8:R87.
Long H, Yin H, Wang L, Gershwin ME, Lu Q. The critical role of epigenetics in systemic lupus erythematosus and autoimmunity. J Autoimmun. 2016;74:118–38.
Vickers KC, Palmisano BT, Shoucri BM, Shamburek RD, Remaley AT. MicroRNAs are transported in plasma and delivered to recipient cells by high-density lipoproteins. Nat Cell Biol. 2011;13:423–33.
Wang K, Zhang S, Weber J, Baxter D, Galas DJ. Export of microRNAs and microRNA-protective protein by mammalian cells. Nucleic Acids Res. 2010;38:7248–59.
Arroyo JD, Chevillet JR, Kroh EM, Ruf IK, Pritchard CC, Gibson DF, Mitchell PS, Bennett CF, Pogosova-Agadjanyan EL, Stirewalt DL, et al. Argonaute2 complexes carry a population of circulating microRNAs independent of vesicles in human plasma. Proc Natl Acad Sci U S A. 2011;108:5003–8.
Ikeda K, Satoh M, Pauley KM, Fritzler MJ, Reeves WH, Chan EK. Detection of the argonaute protein Ago2 and microRNAs in the RNA induced silencing complex (RISC) using a monoclonal antibody. J Immunol Methods. 2006;317:38–44.
Pritchard CC, Kroh E, Wood B, Arroyo JD, Dougherty KJ, Miyaji MM, Tait JF, Tewari M. Blood cell origin of circulating microRNAs: a cautionary note for cancer biomarker studies. Cancer Prev Res (Phila). 2012;5:492–7.
He Y, Lin J, Kong D, Huang M, Xu C, Kim TK, Etheridge A, Luo Y, Ding Y, Wang K. Current state of circulating MicroRNAs as cancer biomarkers. Clin Chem. 2015;61:1138–55.
Parkes JE, Day PJ, Chinoy H, Lamb JA. The role of microRNAs in the idiopathic inflammatory myopathies. Curr Opin Rheumatol. 2015;27:608–15.
Huang X, Yuan T, Tschannen M, Sun Z, Jacob H, Du M, Liang M, Dittmar RL, Liu Y, Liang M, et al. Characterization of human plasma-derived exosomal RNAs by deep sequencing. BMC Genomics. 2013;14:319.
Sekigawa I, Yanagida M, Iwabuchi K, Kaneda K, Kaneko H, Takasaki Y, Jung G, Sone S, Tanaka Y, Ogawa H, Takamori K. Protein biomarker analysis by mass spectrometry in patients with rheumatoid arthritis receiving anti-tumor necrosis factor-alpha antibody therapy. Clin Exp Rheumatol. 2008;26:261–7.
Nozawa K, Fujishiro M, Kawasaki M, Yamaguchi A, Ikeda K, Morimoto S, Iwabuchi K, Yanagida M, Ichinose S, Morioka M, et al. Inhibition of connective tissue growth factor ameliorates disease in a murine model of rheumatoid arthritis. Arthritis Rheum. 2013;65:1477–86.
Kawasaki M, Fujishiro M, Yamaguchi A, Nozawa K, Kaneko H, Takasaki Y, Takamori K, Ogawa H, Sekigawa I. Possible role of the JAK/STAT pathways in the regulation of T cell-interferon related genes in systemic lupus erythematosus. Lupus. 2011;20:1231–9.
Bohan A, Peter JB. Polymyositis and dermatomyositis (first of two parts). N Engl J Med 1975;292:344-347.
Bohan A, Peter JB. Polymyositis and dermatomyositis (second of two parts). N Engl J Med 1975;292:403-407.
Hay EM, Bacon PA, Gordon C, Isenberg DA, Maddison P, Snaith ML, Symmons DP, Viner N, Zoma A. The BILAG index: a reliable and valid instrument for measuring clinical disease activity in systemic lupus erythematosus. Q J Med. 1993;86:447–58.
Castro-Villegas C, Perez-Sanchez C, Escudero A, Filipescu I, Verdu M, Ruiz-Limon P, Aguirre MA, Jimenez-Gomez Y, Font P, Rodriguez-Ariza A, et al. Circulating miRNAs as potential biomarkers of therapy effectiveness in rheumatoid arthritis patients treated with anti-TNFalpha. Arthritis Res Ther. 2015;17:49.
Carlsen AL, Schetter AJ, Nielsen CT, Lood C, Knudsen S, Voss A, Harris CC, Hellmark T, Segelmark M, Jacobsen S, et al. Circulating microRNA expression profiles associated with systemic lupus erythematosus. Arthritis Rheum. 2013;65:1324–34.
Steen SO, Iversen LV, Carlsen AL, Burton M, Nielsen CT, Jacobsen S, Heegaard NH. The circulating cell-free microRNA profile in systemic sclerosis is distinct from both healthy controls and systemic lupus erythematosus. J Rheumatol. 2015;42:214–21.
Misunova M, Salinas-Riester G, Luthin S, Pommerenke C, Husakova M, Zavada J, Klein M, Plestilova L, Svitalkova T, Cepek P, et al. Microarray analysis of circulating micro RNAs in the serum of patients with polymyositis and dermatomyositis reveals a distinct disease expression profile and is associated with disease activity. Clin Exp Rheumatol. 2016;34:17–24.
Inoue K, Jinnin M, Yamane K, Makino T, Kajihara I, Makino K, Honda N, Nakayama W, Fukushima S, Ihn H. Down-regulation of miR-223 contributes to the formation of Gottron’s papules in dermatomyositis via the induction of PKCvarepsilon. Eur J Dermatol. 2013;23:160–7.
van Rooij E, Liu N, Olson EN. MicroRNAs flex their muscles. Trends Genet. 2008;24:159–66.
Georgantas RW, Streicher K, Greenberg SA, Greenlees LM, Zhu W, Brohawn PZ, Higgs BW, Czapiga M, Morehouse CA, Amato A, et al. Inhibition of myogenic microRNAs 1, 133, and 206 by inflammatory cytokines links inflammation and muscle degeneration in adult inflammatory myopathies. Arthritis Rheumatol. 2014;66:1022–33.
Sontheimer RD. Skin manifestations of systemic autoimmune connective tissue disease: diagnostics and therapeutics. Best Pract Res Clin Rheumatol. 2004;18:429–62.
Ghirardello A, Bassi N, Palma L, Borella E, Domeneghetti M, Punzi L, Doria A. Autoantibodies in polymyositis and dermatomyositis. Curr Rheumatol Rep. 2013;15:335.
This study was supported in part by a Grant-in-Aide (S1311011) from MEXT (Ministry of Education, Culture, Sports, Science and Technology)-supported Program for the Strategic Research Foundation at Private University, and by JSPS (Japan Society for the Promotion of Science) KAKENHI Grant Number JP15K08107 (to KI) and JP26860265 (to KH).
Availability of data and materials
All microRNA array data are available at the National Centre for Biotechnology Information Gene Expression Omnibus site (www.ncbi.nlm.nih.gov/geo; accession no. GSE94535).
Ethics approval and consent to participate
Ethical approval for this study was granted by the institutional review board of the Juntendo University Urayasu Hospital. Written informed consent was obtained from all the patients according to the Declaration of Helsinki.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Serum autoantibody profiles of PM/DM patients in this study. (DOCX 51 kb)