Research article | Open | Published:
miR-155 in the progression of lung fibrosis in systemic sclerosis
Arthritis Research & Therapyvolume 18, Article number: 155 (2016)
MicroRNA (miRNA) control key elements of mRNA stability and likely contribute to the dysregulated lung gene expression observed in systemic sclerosis associated interstitial lung disease (SSc-ILD). We analyzed the miRNA gene expression of tissue and cells from patients with SSc-ILD. A chronic lung fibrotic murine model was used.
RNA was isolated from lung tissue of 12 patients with SSc-ILD and 5 controls. High-resolution computed tomography (HRCT) was performed at baseline and 2–3 years after treatment. Lung fibroblasts and peripheral blood mononuclear cells (PBMC) were isolated from healthy controls and patients with SSc-ILD. miRNA and mRNA were analyzed by microarray, quantitative polymerase chain reaction, and/or Nanostring; pathway analysis was performed by DNA Intelligent Analysis (DIANA)-miRPath v2.0 software. Wild-type and miR-155 deficient (miR-155ko) mice were exposed to bleomycin.
Lung miRNA microarray data distinguished patients with SSc-ILD from healthy controls with 185 miRNA differentially expressed (q < 0.25). DIANA-miRPath revealed 57 Kyoto Encyclopedia of Genes and Genomes pathways related to the most dysregulated miRNA. miR-155 and miR-143 were strongly correlated with progression of the HRCT score. Lung fibroblasts only mildly expressed miR-155/miR-21 after several stimuli. miR-155 PBMC expression strongly correlated with lung function tests in SSc-ILD. miR-155ko mice developed milder lung fibrosis, survived longer, and weaker lung induction of several genes after bleomycin exposure compared to wild-type mice.
miRNA are dysregulated in the lungs and PBMC of patients with SSc-ILD. Based on mRNA-miRNA interaction analysis and pathway tools, miRNA may play a role in the progression of the disease. Our findings suggest that targeting miR-155 might provide a novel therapeutic strategy for SSc-ILD.
Systemic sclerosis (SSc) is a chronic autoimmune disease in which pulmonary involvement, particularly interstitial lung disease (SSc-ILD), is known to be the leading cause of mortality [1–3]. The prevalence of SSc-ILD is high and approximately 15–30 % of patients will progress to severe lung fibrosis. However, distinguishing those who will progress from those who will have slow or stable disease remains challenging [4, 5]. Using mRNA expression analysis of lung tissue in a prospective cohort of patients with SSc-ILD, we recently described immune activation pathways and upregulated transforming growth factor-β (TGFβ)-induced signature genes associated with progressive lung fibrosis .
MicroRNA (miRNA) are small RNA that are of vital importance in regulating gene expression within cells, having critical roles in innate and adaptive immunity in general . In addition, there is increasing evidence they play a role in fibrosis and are strongly associated with the pathogenesis of a wide range of human diseases [8, 9]. miRNA are remarkably stable and reliable for detection in the blood and urine and bring new hope to the development of clinical biomarkers. More importantly, disease-associated miRNA represent a new class of therapeutic targets.
Using high-throughput technology, we show that many miRNA are dysregulated in the lungs of patients with SSc-ILD. In addition, selected miRNA are associated with progressive lung fibrosis based on image lung score and lung function tests, with miR-155 as one of the best markers. miR-155-deficient mice survive longer and develop milder lung fibrosis. Therefore, miR-155 could be a potential novel therapeutic target for lung fibrosis in patients with SSc.
Between January 2002 and July 2004, 28 consecutive patients with SSc-ILD, identified on HRCT as having associated respiratory symptoms and/or reduced performance on a pulmonary function test (PFT) underwent an open lung biopsy and were followed for at least 3 years. Twelve subjects involved in the previous study  were selected (sixteen samples were not included due to lack of lung material for miRNA isolation). The patients included in the current study were classified as having diffuse SSc (dSSc) (n = 7) and limited SSc (lcSSc) (n = 5), according to diagnostic  and subtype criteria . Controls were five samples of normal histological healthy lung tissue obtained during lung cancer surgical resection. Lung explants from controls and patients with SSc-ILD were obtained from lung transplants at the University of Pittsburgh Medical Center (healthy controls, n = 6; SSc-ILD, n = 6) and used for fibroblast culture. Peripheral blood mononuclear cells (PBMC) were collected from healthy controls (n = 6) and patients with SSc-ILD (n = 13) from Boston University Medical Center. Informed consent was obtained from all patients and healthy subjects involved in the study. The Boston University Medical Center, the University of Sao Paulo, and the University of Pittsburgh Institutional Review Boards reviewed and approved the conduct of this study.
Two pathologists (VLC and ERP) classified the lung specimens as previously described .
High-resolution computed tomography
Two pulmonologists (CRC and RC) performed blinded evaluation of the images. As described previously  each lung was divided into three zones: upper (lung apex to aortic arch); middle (aortic arch to inferior pulmonary veins); and lower (inferior pulmonary veins to lung bases). For baseline and follow up HRCT, the extent of the pulmonary abnormality was scored as previously described . The sum of all scores on the lung zones is referred to as FibMax. A delta FibMax, which is the difference between the follow-up FibMax and the baseline FibMax, was used to evaluate the progression of the lung fibrosis.
Open lung biopsy was performed by formal thoracotomy, avoiding honeycombing areas .
Patients with evidence of non-specific interstitial pneumonia (NSIP) on lung biopsy were treated with monthly intravenous cyclophosphamide (0.5–1.0 g/m2) for at least one year (see previous report for details ). HRCT and PFT were performed before biopsy and repeated 2–3 years after the one-year treatment (Fig. 1a).
RNA and microRNA isolation and microarray analysis hybridization
Total RNA and miRNA from human lung biopsy tissue, lung fibroblasts, and PBMC were extracted using microRNeasy Mini Kits (Qiagen). RNA samples were stored at –80 °C.
MicroRNA and messenger RNA quantitative PCR (qPCR)
Complementary DNA (cDNA) was transcribed using the TaqMan MicroRNA Reverse Transcription Kit (Life Technologies) and Reverse Transcriptase (Gibco BRL) and primers using 200 ng of total RNA and according to the manufacturer’s protocol. Samples were diluted 1:15 and followed the TaqMan Universal PCR Master Mix (Life Technologies) protocol. Specific miRNA primers were used for miR-155, miR-21, miR-182, miR-15b, miR-193a, and RNU6B (Life Technologies). Specific mRNA primers were used for collagen type 3 alpha-1 (Col3a1), MS4A4A, periostin (POSTN), and 18S (Life Technologies). miRNA and mRNA qPCR was performed using the TaqMan assay and StepOne Plus Real Time PCR system (Applied Biosystems) in a 20-ul reaction for 40 cycles. Gene expression was analyzed using the difference in cycle threshold (ΔCt) method. The Ct values of the miRNAs were normalized to RNU6B and of the mRNA were normalized to 18S as an internal control using the equation:
ΔCt = Ct reference – Ct target and expressed as ΔCt.
Microarray data clustering
All procedures were performed at Boston University Microarray Resource Facility as described in the Affymetrix GeneChip 3′IVT Express user manual (Affymetrix, Santa Clara, CA, USA; www.affymetrix.com). For mRNA arrays biotin-labeled amplified total RNA was purified, fragmented, and hybridized to Affymetrix U133A 2.0 microarrays. The MAS5 algorithm with global scaling normalization was used to generate gene-level expression data. For miRNA arrays miRNA was purified, fragmented and hybridized to Affymetrix GenChip miRNA 3.0 microarrays. The signal of the samples was amplified and microarrays were immediately scanned using Affymetrix GeneArray Scanner 3000 7G Plus (Affymetrix, Santa Clara, CA, USA). The resulting CEL files were summarized using the Affymetrix Expression Console (current version 1.1). Clustering was performed using Cluster 3.0 software. Both mRNA and miRNA Affymetrix data [GEO:GSE81294] and SubSeries linked to GSE81294 [GEO:GSE81292; GEO:GSE81293] are available at the public repository Gene Expression Omnibus. Additional file 1: Table S1 and Additional file 2: Table S2 contain all the miRNA and mRNA, respectively, that were significantly different in patients with SSc-ILD and controls.
Lung tissue was minced and primary fibroblasts cultured using the outgrowth method as previously described . Lung fibroblasts were cultured in DMEM supplemented with 10 % fetal bovine serum and penicillin/streptomycin and utilized at passages 2–4. Fibroblasts (100 % confluent) were incubated in serum-free DMEM overnight prior to stimulation with TGFß (R&D System; 2.5 ng/ml), recombinant human IL-13 (R&D Systems, 20 ng/ml), or interferon-alpha (IFN) (R&D Systems; 500 U/ml) for 18 hours. Total RNA from fibroblasts was transferred in Qiazol buffer and purified using the miRNease mini kit protocol (Qiagen). RNA samples were stored at –80 °C.
Peripheral blood mononuclear cells
Blood was collected in cell preparation tubes (CPTs) designed for one-step cell separation (Becton Dickinson) from healthy controls and patients within 3 months of the date when the lung function tests were performed. The sample was immediately mixed and centrifuged at 1800 g at ambient temperature for 30 minutes. Total RNA from PBMC was extracted using the miRNease mini kit protocol (Qiagen). RNA samples were stored at –80 °C.
MicroRNA Nanostring technology
RNA (100 ng from each sample) was used for the miRNA analysis. The miRNA library contains 800 of the most relevant miRNA described in the literature. Normalization was performed using the average of expression of the top 20 most-expressed miRNA in all samples. The minimal threshold for the detection of miRNA was considered as 50 counts.
RNA Nanostring technology
RNA (100 ng of from each mouse) was used for RNA analysis. The lung gene panel was built based on a microarray analysis of the lungs of mice exposed to a bleomycin-pump compared to PBS in previous experiments (unpublished results). Normalization was performed using nine housekeeping genes that were not affected by the bleomycin-pump model.
DNA Intelligent Analysis (DIANA) DIANA-miRPath v2.0 utilizes miRNA targets based on DIANA-microT-CDS, which predicts miRNA-gene interaction including the binding region, position, and type, and/or miRNA-gene interaction experimentally validated from Tarbase v7.0 . Combined analysis results for miRNA and pathway-related information was obtained from miRBase 18  and the Kyoto Encyclopedia of Genes and Genomes (KEGG) v58 . The selected top 30 upregulated and downregulated miRNA were used (Additional file 1: Table S1).
Bleomycin lung fibrosis model
C57BL/6 J (B6) and B6.Cg-Mir155tm1.Rsky/J mice were purchased from The Jackson Laboratory (Bar Harbor, ME, USA). Age-matched and sex-matched mice were used in the experiments at 8 to 10 weeks. Mice were kept in pathogen-free conditions with food and water. The in vivo study was approved by the Animal Ethics Committee at the Boston University School of Medicine. Mice were chronically exposed to PBS or Bleomycin (90 U/Kg) by subcutaneous osmotic pump injections for 7 days. Mice were observed and sacrificed on day 28. The survival curve was calculated. Lungs were harvested for histology and gene expression. H&E and Masson’s Trichrome staining were performed. The Ashcroft lung score was measured blinded by one investigator (RBC).
Data analysis was performed using GraphPad Prism 6 software, version 6.0a (GraphPad Software, San Diego, CA, USA). Multiple group comparisons were analyzed by one-way analysis of variance (ANOVA) with Bonferroni adjustment for multiple comparisons. Two-group comparisons were analyzed by Student’s t test or the Mann-Whitney test. Correlation was tested and presented as Pearson’s correlation coefficient (r). The survival fraction was calculated by the Kaplan-Meier test with 95 % confidence intervals computed by the asymmetrical method using the log-rank test.
Baseline characteristics of patients with systemic sclerosis interstitial lung disease
A timeline of the study is summarized in Fig. 1a. We have previously described all the demographic and FibMax data of these patients with SSc-ILD before and after treatment . Note that despite treatment, most patients with SSc-ILD developed progressive disease based on the FibMax, with a mean baseline score of 6.88 ± 4.47 compared to 12.76 ± 6.59 after treatment; p = 0.006. The delta FibMax was used to gauge the progression of lung fibrosis as previously described . Healthy control lungs were obtained from cancer-free resected lungs.
Patients with systemic sclerosis interstitial lung disease have distinct microRNA gene expression patterns
We analyzed lung tissue samples from healthy controls (n = 5) and patients with SSc-ILD (n = 12) using miRNA microarrays. Both middle and lower lobe samples were analyzed in five of the patients with SSc-ILD, while only the lower lobes were analyzed in seven patients. One patient was excluded from the miRNA analysis on the basis of poor quality metrics. To identify genes with differential expression associated with SSc-ILD, Student’s two-sample t test was performed to compare SSc-ILD (expression averaged for each gene from both lung biopsy sites) and control samples (p < 0.05). A false discovery rate (FDR)-corrected p value (q < 0.25) was computed after removing miRNA present in ≤25 % of the samples, to eliminate miRNA with low overall expression. A heatmap of the most informative miRNA (n = 185 genes) was created containing the miRNA with the most significant differential expression between SSc-ILD and control samples. Most of the miRNA were downregulated in patients compared to the healthy controls (Fig. 1b). The complete list of significantly different miRNA can be found in Additional file 1: Table S1.
We validated several miRNA that were upregulated in patients with SSc-ILD compared to the controls, by qPCR. miR-182 was highly expressed in the lung tissue from patients with SSc-ILD (mean fold-change 17.5 ± 14.8) compared to controls (mean fold-change 1.1 ± 0.6; p = 0.001. Fig. 2a). miR-155 was upregulated in patients with SSc-ILD (mean fold-change 7.6 ± 9.2) compared to controls (mean fold-change 1.0 ± 0.5; p = 0.003. Fig. 2b), and miR-21 was upregulated in patients with SSc-ILD (mean fold-change 3.4 ± 2.6) compared to controls (mean fold-change 1.1 ± 0.7; p = 0.06), although not statistically significant (Fig. 2c).
Validation of relevant gene signatures in systemic sclerosis interstitial lung disease
In order to better understand the potential mRNA targeted by altered lung miRNA expression, we analyzed all samples for mRNA expression by microarray. We identified genes with differential expression associated with SSc-ILD, using Student’s two-sample t test (p < 0.05) corrected by FDR (q < 0.25). These microarray data confirmed the altered expression of hundreds of genes in SSc-ILD that we reported previously , such as chemokine C-C motif ligand 18 (CCL18), several collagens (types I and III), macrophage activation markers (for example, CD163 and MS4A4A), secreted phosphoprotein 1/osteopontin (SPP1), and others (Additional file 2: Table S2). Using qPCR, we validated several of the top upregulated genes in lung tissue from patients with SSc-ILD, such as COL3a1 (p = 0.004), MS4A4A (p = 0.02), and POSTN (p = 0.006) Fig. 2d-f.
Messenger RNA and microRNA: a complex network
Taking advantage of our simultaneous analysis of mRNA and miRNA in the same cohort of patients with SSc-ILD and healthy controls we further investigated interactions between mRNA and miRNA lung expression. Pearson’s correlation was tested between selected upregulated and downregulated miRNA (see Additional file 1: Table S1) and the most highly upregulated and downregulated mRNA in SSc-ILD (see Additional file 2: Table S2) by microarray analysis, based on average fold-change gene expression compared to controls. miR-155, miR-21, and miR-4459 were the miRNA that were most strongly correlated and their interactions are shown in Table 1.
RNA and microRNA: pathway analysis
We selected the top upregulated and downregulated miRNA in patients with SSc-ILD compared to the controls (see Additional file 1: Table S1) for pathway analysis. First, we applied the DIANA-Tarbase v7.0 database in order to identify miRNA target genes experimentally validated in the literature (see “Methods”). Only a few miRNA from the selected group had validated target genes: miR-155, miR-21, miR-125a, miR-31, miR-205, miR-20a, let-7f, miR-182, let-7e, miR-199b, miR-199a, miR-663, and miR-193a. miR-155 (n = 812) and miR-21 (n = 581) had the largest list of target genes already validated.
To gain further insight into the function of the top dysregulated miRNA and their targets, DIANA-microT-CDS miRPath was applied to identify the significant KEGG pathways (see “Methods”): 57 KEGG biological processes were significantly enriched (p < 0.05, FDR-corrected) in SSc-ILD (Additional file 3: Figure S1). The 10 most highly implicated processes identified using this approach, with their related miRNA, can be found in Table 2.
miR-155 lung expression correlates with progressive systemic sclerosis interstitial lung disease
There is an urgent need for a biomarker that predicts the progression of lung fibrosis in patients with SSc. miRNA are known to regulate a complex network; therefore, we tested the correlation between miRNA microarray lung gene expression and changes in the lung fibrosis score FibMax, the sum of all fibrotic scores in the six lung zones (see “Methods”). The delta FibMax of each patient, when available, was correlated with corresponding patient lung miRNA microarray fold-change gene expression (see Additional file 1: Table S1). Surprisingly, only four miRNA correlated positively with the delta FibMax (r > 0.4) with miR-155 having the strongest correlation (r = 0.65, p < 0.001, Fig. 3a), followed by miR-182 (r = 0.49, p = 0.06), miR-27a (r = 0.49, p = 0.06), and miR-21 (r = 0.47, p = 0.07). Five miRNA correlated negatively with the delta FibMax, with miR-143 having the strongest negative correlation (r = –0.64, p < 0.01), followed by miR-4270 (r = –0.59, p = 0.01), miR-4530 (r = –0.54, p = 0.03), miR-19b (r = –0.51, p = 0.05), and miR-4459 (r = –0.50, p = 0.05).
miR-155 expression by peripheral blood mononuclear cells is correlated with progressive systemic sclerosis interstitial lung disease
PBMC gene expression has been used as a surrogate for rarely available lung tissue to further investigate lung complications in SSc, such as SSc-ILD and SSc pulmonary arterial hypertension (SSc-PAH) . Therefore, we analyzed miRNA expression in SSc-ILD (n = 13) and control PBMC (n = 5) by Nanostring; demographic features are shown in Additional file 4: Table S3. Among the miRNA that were upregulated in the lungs of patients with SSc-ILD (Additional file 1: Table S1) let-7d (p = 0.01), miR-15b (p = 0.05), and miR-21 (p = 0.1) were the most expressed in PBMC from SSc-ILD compared to controls (Additional file 5: Table S4) based on fold-change expression.
More importantly, we observed that miR-155 expression on SSc-ILD PBMC was strongly negatively correlated with both %FVC (r = –0.60; p = 0.01) and %DLCO (r = –0.58; p = 0.02). let-7d was also strongly negatively correlated with %FVC (r = –0.65; p = 0.007) and with %DLCO (r = –0.56; p = 0.02). miR-21 PBMC expression did not correlate with lung function tests (%FVC r = 0.14; %DLCO r = 0.18; both p > 0.05).
Activated lung fibroblasts have minor influence on the microRNA signature in systemic sclerosis interstitial lung disease
In order to elucidate whether fibroblasts might be responsible for the miRNA expression in SSc-ILD we stimulated healthy lung fibroblasts (n = 6) and fibroblasts from the lungs of patients with SSc-ILD (n = 6) with cytokines known to be involved in the progression of lung disease: TGFβ, IL-13, and IFN . Fibroblasts were stimulated for 18 hours and miRNA expression was analyzed by Nanostring technology (Additional file 6: Table S5) and confirmed in a larger number of samples by qPCR (see “Methods”).
Expression of miR-21 was in general slightly increased after IFN and TGFβ stimulation of fibroblasts in SSc-ILD, although not statistically significant (Fig. 3d, p > 0.05). On the other hand, miR-155 expression was reduced overall after all three stimuli of fibroblasts in SSc-ILD, although not statistically significant (Fig. 3e, p > 0.05). There was mild expression of miR-193a and miR-15b in lung fibroblasts in both groups after all three stimuli (Additional file 7: Figure S2A and B, respectively).
The absence of miR-155 protects mice from bleomycin-induced lung fibrosis
We observed that miR-155 expression in the lungs and in the blood of patients with SSc-ILD are strongly correlated with progressive lung disease identified by image and by lung function tests, respectively. Therefore, we tested the bleomycin model of chronic lung fibrosis in mice deficient in miR-155 (miR-155 knockout (KO)) compared to wild-type mice (WT). miR-155 KO mice had a 100 % survival rate over the course of 28 days after starting bleomycin compared to a 45 % survival rate in WT mice; p < 0.001, log-rank test (Fig. 4a). WT (n = 7) and miR-155 KO (n = 6) mice exposed to PBS were used as controls. WT (n = 11) and miR-155 KO (n = 11) mice were exposed to bleomycin and three WT mice were found dead on days 8, 13, and 17, respectively. One WT mouse exposed to bleomycin was killed on day 20 due to poor conditions and was included in the lung Nanostring analysis. miR-155 KO mice (n = 11) developed milder fibrosis in the lungs as measured by the Ashcroft lung score (Fig. 4b, p < 0.001) compared to the WT (n = 7) exposed to bleomycin on day 28.
We tested the effect of miR-155 deletion on bleomycin-induced gene expression using a custom-designed Nanostring that we developed to detect proinflammatory and profibrotic genes upregulated in the lungs of bleomycin-treated lung tissue and lung tissue from patients with SSc-ILD. Overall, bleomycin-induced lung gene expression was lower in miR-155 KO mice (Additional file 8: Figure S3). Arginase-1 lung gene expression was strongly induced in the WT mice and correlated positively with the Ashcroft score (r = 0.70, p < 0.001). Arginase-1, a marker of M2 macrophages (Fig. 4c), and tissue inhibitor of metalloproteinase-1 (TIMP1) lung expression (Fig. 4d) were induced significantly less after bleomycin in miR155 KO mice compared to WT mice. Expression of other bleomycin-induced genes was blunted in miR-155 KO lungs compared to WT mice and clustered with Arginase-1 and TIMP1 (Fig. 4e). On the other hand, expression of CD68, a general marker of macrophages, was induced similarly in both groups exposed to bleomycin (p > 0.05).
Dysregulation of the immune system is critical for the development of skin and lung involvement in SSc. In this regard, we have recently shown that this inappropriate inflammatory response in the lungs is strongly correlated with progressive lung fibrosis , along with a profibrotic TGFβ signature. miRNA are a new class of regulators that tightly control the overall inflammatory response. We show here that miRNA are intensely dysregulated in the lungs and in PBMC in patients with SSc-ILD. Moreover, we found that dysregulated gene expression of only a few miRNA, including miR-155, are strongly associated with progressive lung disease and that the absence of miR-155 is protective for the bleomycin-lung fibrotic model.
miR-155 controls differentiation of CD4+ T cells into its several subtypes to regulate B-cell differentiation and antibody production, and is highly expressed in activated monocytes/macrophages . miR-155 is upregulated in several autoimmune disorders, such as rheumatoid arthritis , multiple sclerosis , systemic lupus erythematosus [20, 21], and ulcerative colitis . In our cohort, we showed that miR-155 expression correlates highly with profibrotic gene expression, such as SPP1 and POSTN. We also observed strong correlation between miR-155 PBMC expression in patients with SSc-ILD and lung function tests. On the other hand, lung fibroblasts seem not to be the most relevant cell type driving the strong expression of miR-155 in the lungs in patients with SSc-ILD.
We showed the functional relevance of miR-155 in lung fibrotic disease where miR-155 KO mice survived longer and developed substantially less aggressive lung fibrosis based on lung phenotype and on a comprehensive lung gene expression analysis. This protective effect might be underestimated as three mice in the WT bleomycin group did not survive until day 28 and neither the lung score nor the lung gene expression were considered in the analysis. In addition, the absence of miR-155 in the model blocked the alternative activation of lung macrophages, which correlates strongly with progressive lung fibrosis in patients with SSc-ILD  and with the Ashcroft lung fibrotic score in our murine model. Yan et al.  recently showed protection of bleomycin cutaneous fibrosis in miR-155-deficient mice and reduced skin thickness after topical antagomiR-155 treatment in mice primed with subcutaneous bleomycin . In the intratracheal bleomycin murine lung model miR-155 was upregulated and correlated with the degree of lung fibrosis . Altogether, our observations reinforce the hypothesis of a strong immune activation in SSc-ILD with miR-155 as a key regulator.
Our observations that miR-21 is upregulated in the lungs of patients with SSc-ILD suggest that modulating miR-21 might be a future therapeutic strategy for SSc-ILD. miR-21 has been extensively studied in fibrosis, with increased expression observed in the lungs of bleomycin-treated mice and in the lungs of patients with idiopathic pulmonary fibrosis (IPF) in whom it is mainly localized in myofibroblasts . More importantly, restoring miR-21 to normal levels inhibits bleomycin-induced pulmonary fibrosis . miR-21 is also upregulated in the skin and fibroblasts from patients with both diffuse and limited subtypes of SSc . miR-21 is induced by TGFβ in stellate hepatic cells and in a feed-forward loop it amplifies its signal, eventually promoting fibrosis . Lung fibroblasts might not be the main cell type responsible for the miRNA signature in SSc-ILD, as we observed only a slight increase in miR-21 expression after TGFβ stimulation. Although our data on lung tissues in SSc-ILD do not permit a direct assessment of intracellular pathways, we showed that several upregulated profibrotic genes, such as Col3a1 and POSTN, correlate positively with miR-21 lung expression. In addition to being well-established as a profibrotic-miR, miR-21 is also central to many inflammatory pathways, and most importantly, in controlling the Toll-like receptor signaling with strong connection with miR-155 .
We also showed for the first time that miR-182 is upregulated in the lungs of patients with SSc-ILD. miR-182 has been implicated in several types of cancer and it is known to regulate proliferation, invasion, and migration of cancer cells . Interestingly enough, some of the recently identified miR-182 targets are strong inhibitors of extracellular matrix degradation . miR-4484 and miR-4459 were significantly downregulated in our cohort of patients with SSc-ILD and neither of them have validated target genes based on Tarbase v7.0 .
Our study has limitations. Our cohort was small and our results should ideally be validated in a larger group. The suggested mRNA-miRNA interactions and pathways are mainly predicted, although our agnostic and unbiased analysis using algorithmic tools guided us in similar directions, which support the validity of our results.
Our unique approach of analyzing the relevance of miRNA in SSc-ILD using gene expression, pathway analysis, and functional studies, revealed a selected group of miRNA that might mediate the development and/or the progression of SSc-ILD. miR-155 and miR-21 were the most relevant miRNA having in common a close involvement in altered pathways observed in fibrotic diseases. mir-155 regulates the development of lung fibrosis in mice, thus, opening a new drug development opportunity.
%DLCO, percentage of diffusing capacity of the lungs for carbon monoxide; %FVC, percentage of forced vital capacity; ANOVA, analysis of variance; CCL18, chemokine C-C motif ligand 18; cDNA, complementary DNA; Col3a1, collagen type 3 alpha-1; Ct, cycle threshold; DIANA, DNA Intelligent Analysis; DMEM, Dulbecco’s modified Eagle’s medium; dSSc, diffuse systemic sclerosis; H&E, hematoxylin and eosin; FDR, false discovery rate; FiBMax, high-resolution computed tomography lung score; HRCT, high-resolution computed tomography; IFN, interferon alpha; IL-13, interleukin 13; KEGG, Kyoto Encyclopedia of Genes and Genomes; KO, knockout; lSSc, limited systemic sclerosis; mRNA, messenger RNA; miRNA, micro RNA; MS4A4A, membrane-spanning 4-domains, subfamily A, member 4A; NSIP, non-specific interstitial pneumonia; PBMC, peripheral blood mononuclear cells; PBS, phosphate-buffered saline; PFT, pulmonary function tests; POSTN, periostin; SPP1, osteopontin-1; SSc, systemic sclerosis; SSc-ILD, systemic sclerosis interstitial lung disease; TGF-beta, transforming growth factor beta; TIMP-1, tissue inhibitor of metalloproteinases type-1; WT, wild-type; IPF, idiopathic pulmonary fibrosis
Ioannidis JP, Vlachoyiannopoulos PG, Haidich AB, Medsger Jr TA, Lucas M, Michet CJ, et al. Mortality in systemic sclerosis: an international meta-analysis of individual patient data. Am J Med. 2005;118(1):2–10.
Tyndall AJ, Bannert B, Vonk M, Airo P, Cozzi F, Carreira PE, et al. Causes and risk factors for death in systemic sclerosis: a study from the EULAR Scleroderma Trials and Research (EUSTAR) database. Ann Rheum Dis. 2010;69(10):1809–15
Sampaio-Barros PD, Bortoluzzo AB, Marangoni RG, Rocha LF, Del Rio AP, Samara AM, et al. Survival, causes of death, and prognostic factors in systemic sclerosis: analysis of 947 Brazilian patients. J Rheumatol. 2012;39(10):1971–78.
Roth MD, Tseng CH, Clements PJ, Furst DE, Tashkin DP, Goldin JG, et al. Predicting treatment outcomes and responder subsets in scleroderma-related interstitial lung disease. Arthritis Rheum. 2011;63(9):2797–08.
Steen VD, Medsger TA. Changes in causes of death in systemic sclerosis, 1972-2002. Ann Rheum Dis. 2007;66(7):940–4.
Christmann RB, Sampaio-Barros P, Stifano G, Borges CL, de Carvalho CR, Kairalla R, et al. Association of Interferon- and transforming growth factor beta-regulated genes and macrophage activation with systemic sclerosis-related progressive lung fibrosis. Arthritis Rheumatol. 2014; 66(3):714–25.
Baumjohann D, Ansel KM. MicroRNA-mediated regulation of T helper cell differentiation and plasticity. Nat Rev Immunol. 2013;13(9):666–78.
Pandit KV, Milosevic J, Kaminski N. MicroRNAs in idiopathic pulmonary fibrosis. Transl Res. 2011;157(4):191–99.
Booton R, Lindsay MA. Emerging role of MicroRNAs and long noncoding RNAs in respiratory disease. Chest. 2014;146(1):193–04.
Preliminary criteria for the classification of systemic sclerosis (scleroderma). Subcommittee for scleroderma criteria of the American Rheumatism Association Diagnostic and Therapeutic Criteria Committee. Arthritis Rheum. 1980;23(5):581-90.
Furst DE, Clements PJ, Steen VD, Medsger Jr TA, Masi AT, D’Angelo WA, et al. The modified Rodnan skin score is an accurate reflection of skin biopsy thickness in systemic sclerosis. J Rheumatol. 1998;25(1):84–8.
Hsu E, Shi H, Jordan RM, Lyons-Weiler J, Pilewski JM, Feghali-Bostwick CA. Lung tissues in patients with systemic sclerosis have gene expression patterns unique to pulmonary fibrosis and pulmonary hypertension. Arthritis Rheum. 2011;63(3):783–94.
Vlachos IS, Paraskevopoulou MD, Karagkouni D, Georgakilas G, Vergoulis T, Kanellos I, et al. DIANA-TarBase v7.0: indexing more than half a million experimentally supported miRNA:mRNA interactions. Nucleic Acids Res. 2015;43(Database issue):D153–159.
Kozomara A, Griffiths-Jones S. miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res. 2011;39(Database issue):D152–7.
Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012;40(Database issue):D109–14.
Christmann RB, Hayes E, Pendergrass S, Padilla C, Farina G, Affandi AJ, et al. Interferon and alternative activation of monocyte/macrophages in systemic sclerosis-associated pulmonary arterial hypertension. Arthritis Rheum. 2011;63(6):1718–28.
Elton TS, Selemon H, Elton SM, Parinandi NL. Regulation of the MIR155 host gene in physiological and pathological processes. Gene. 2013;532(1):1–12.
Leah E. Rheumatoid arthritis: miR-155 mediates inflammation. Nat Rev Rheumatol. 2011;7(8):437.
Moore CS, Rao VT, Durafourt BA, Bedell BJ, Ludwin SK, Bar-Or A, et al. miR-155 as a multiple sclerosis-relevant regulator of myeloid cell polarization. Ann Neurol. 2013;74(5):709–20.
Wang G, Tam LS, Li EK, Kwan BC, Chow KM, Luk CC, et al. Serum and urinary cell-free MiR-146a and MiR-155 in patients with systemic lupus erythematosus. J Rheumatol. 2010;37(12):2516–22.
Zhou S, Wang Y, Meng Y, Xiao C, Liu Z, Brohawn P, et al. In Vivo Therapeutic Success of MicroRNA-155 (miR-155) Antagomir in a Mouse Model of Lupus Alveolar Hemorrhage. Arthritis Rheumatol. 2015. doi:10.1002/art.39485 [Epub ahead of print].
Takagi T, Naito Y, Mizushima K, Hirata I, Yagi N, Tomatsuri N, et al. Increased expression of microRNA in the inflamed colonic mucosa of patients with active ulcerative colitis. J Gastroenterol Hepatol. 2010;25 Suppl 1:S129–133.
Yan Q, Chen J, Li W, Bao C, Fu Q. Targeting miR-155 to treat experimental scleroderma. Sci Rep. 2016;6:20314.
Pottier N, Maurin T, Chevalier B, Puissegur MP, Lebrigand K, Robbe-Sermesant K, et al. Identification of keratinocyte growth factor as a target of microRNA-155 in lung fibroblasts: implication in epithelial-mesenchymal interactions. PLoS One. 2009;4(8), e6718.
Liu G, Friggeri A, Yang Y, Milosevic J, Ding Q, Thannickal VJ, et al. miR-21 mediates fibrogenic activation of pulmonary fibroblasts and lung fibrosis. J Exp Med. 2010;207(8):1589–97.
Zhu H, Li Y, Qu S, Luo H, Zhou Y, Wang Y, et al. MicroRNA expression abnormalities in limited cutaneous scleroderma and diffuse cutaneous scleroderma. J Clin Immunol. 2012;32(3):514–22.
Zhang Z, Zha Y, Hu W, Huang Z, Gao Z, Zang Y, et al. The autoregulatory feedback loop of microRNA-21/programmed cell death protein 4/activation protein-1 (MiR-21/PDCD4/AP-1) as a driving force for hepatic fibrosis development. J Biol Chem. 2013;288(52):37082–93.
Quinn SR, O’Neill LA. A trio of microRNAs that control Toll-like receptor signalling. Int Immunol. 2011;23(7):421–25.
Kouri FM, Hurley LA, Daniel WL, Day ES, Hua Y, Hao L, et al. miR-182 integrates apoptosis, growth, and differentiation programs in glioblastoma. Genes Dev. 2015;29(7):732–45.
Sachdeva M, Mito JK, Lee CL, Zhang M, Li Z, Dodd RD, et al. MicroRNA-182 drives metastasis of primary sarcomas by targeting multiple genes. J Clin Invest. 2014;124(10):4305–19.
We thank Yuriy Alekseyev, Ph.D. (Department of Pathology and Laboratory Medicine) and Adam Gower, Ph.D. (Clinical and Translational Science Institute) of Boston University Microarray Resource Core Facility for supervising all microarray-based experiments and help with data analysis. CTSI award U54-TR001012. This work was supported by 1P50AR060780-01 (RL) and 2R01AR051089-06A1 NIH/NIAMS grants (RL and RS), by the Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) 2013/14277-4 (VLC), by Universal National Council for Scientific and Technological Development (CNPq) 483005/2012-6 (VLC), and by Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), Foundation Ministry of Education of Brazil (TLM). This work was also supported in part by 5T32AI007309-25 NIAID training grant entitled “Research Training in Immunology” (AW), by K-08 Mentored Clinical Scientist Research Career Development Award, 2014 #1K08AR065507 (RBC), by Department of Medicine Career Investment Award, Boston University School of Medicine (RBC), by P30 AR061271 (CFB, SG, and RL), and by R01HL126990 (DJK).
All authors participated in the preparation of the manuscript in a significant way. RBC participated in the design and the coordination of the study, carried out material collection, executed the experiments, performed data analysis, and drafted the manuscript. AW and TLM participated in the design of the study and sample collection. JT and DA executed the lung fibroblast culture experiments. YM carried out the murine model experiments and sample collection. SG executed the lung staining. PSB, CLB, CRRC, and RAK participated in the design of the study and in the collection of the clinical data from baseline and follow up of the patients with SSc-ILD. VLC and EP participated in the design of the study and in the lung histologic diagnosis. JZ and GS carried out material collection and coordinated the study. CFB, DK, and MR provided the lung fibroblasts and participated in the design of the study. RS and RL participated in the design and coordination of the study. All authors read and approved the manuscript after revising it critically for relevant scientific content.
The authors declare that they have no competing interests.
The 185 differential miRNA expression (present in ≥25 % of samples, q < 0.25) in average fold-change observed in SSc-ILD compared to healthy controls. In sequence: transcript ID, fraction of samples. Top 30 upregulated and downregulated miRNA (present in >80 % of samples, p < 0.05) are highlighted in bold. (XLSX 55 kb)
List of mRNA gene expression in SSc-ILD compared to healthy controls. In sequence: gene symbol, gene name, and average fold-change expression in SSc-ILD compared to healthy controls. (XLSX 15 kb)
List of the 57 KEGG pathways related to the most expressed miRNA in the lungs of patients with SSc-ILD compared to controls. In sequence: pathway names, significance levels based on p value, number of targeted genes, number of miRNA targeting each pathway. (PDF 1456 kb)
Demographic data of patients with SSc-ILD and healthy controls for the miRNA analysis in PBMC. %FVC percentage forced vital capacity. %DLCO percentage diffusing capacity of the lung for carbon monoxide, LSSc limited systemic sclerosis, DSSc diffuse systemic sclerosis. (XLSX 35 kb)
miRNA expression in PBMC by Nanostring analysis. Data are expressed as counts ofgenes and fold-change of the average on PBMC in SSc-ILD compared to the average of the controls (column U). (XLSX 167 kb)
miRNA expression in lung fibroblasts by Nanostring analysis. Healthy control lung fibroblasts (HC) and SSc-ILD lung fibroblasts were stimulated with media, interleukin-13 (IL-13), interferon alpha (IFN), and TGF-beta for 18 hours. Data are expressed as counts of gene expression. (XLSX 56 kb)
miRNA expression of miR-193a (A) and miR-15b (B) on lung fibroblasts from cell lines of healthy controls (HC) and patients with SSc-ILD stimulated for 18 hours with media, IL-13, TGF-beta, and IFN-alpha. Data are expressed as fold-change compared to HC samples or media control. ANOVA, p > 0.05 for both miRNAs. (TIF 1112 kb)
Cluster of genes analyzed on lungs from wild-type and miR-155 KO mice exposed to PBS or bleomycin (see “Methods”). Two main clusters of genes upregulated (A) and downregulated (B) on lungs from wild-type mice exposed to bleomycin on day 28 of the model. Each miRNA was z-score-normalized across all samples and scaled to red and blue (≥2 or ≤ –2; respectively) and white indicating a z-score of zero. (TIF 3423 kb)