Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Altered miRNA expression in pulmonary sarcoidosis

BMC Medical Genetics volume 17, Article number: 2 (2016) Cite this article

Abstract

Background

miRNAs control important cellular functions including angiogenesis/angiostasis or fibrosis and reveal altered expression during pathological processes in the lung.

Methods

The aim of the study was to investigate the expression of selected miRNAs (miR-let7f, miR-15b, miR-16, miR-20a, miR-27b, miR-128a, miR-130a, miR-192 miR-221, miR-222) in patients with pulmonary sarcoidosis (n = 94) and controls (n = 50). The expression was assessed by q-PCR in BALF cells and peripheral blood lymphocytes (PB lymphocytes). For statistical analysis, the Kruskal–Wallis test, Mann–Whitney U- test, Neuman–Keuls’ multiple comparison test, and Spearman’s rank correlation were used.

Results

In BALF cells, significantly higher expression of miR-192 and miR-221 and lower expression of miR-15b were found in patients than controls. MiR-27b, miR-192 and miR-221 expression was significantly higher in patients without parenchymal involvement (stages I) than those at stages II-IV. Patients with acute disease demonstrated significantly higher miR-27b, miR-192 and miR-221 expression than those with insidious onset. For PB lymphocytes, patients demonstrated significantly greater miR-15b, miR-27b, miR-192, miR-221 and miR-222 expression, but lower miR-let7f and miR-130a expression, than controls. Stage I patients demonstrated significantly higher miR-16 and miR-15b expression than those in stages II-IV, and patients with the acute form demonstrated higher miR-130a and miR-15b expression. In BALF cells, miR-16 and miR-20a expression was significantly higher in patients with lung volume restriction, and miR-let7f was higher in the PB lymphocytes in patients with obturation. Several correlations were observed between the pattern of miRNA expression, lung function parameters and selected laboratory markers.

Conclusion

The obtained results suggest that the studied miRNAs play a role in the pathogenesis of sarcoidosis, and that some of them might have negative prognostic value.

Peer Review reports

Background

Sarcoidosis is a chronic inflammatory granulomatous disorder of unknown etiology, characterized in the majority of patients by intrathoracic lymph node enlargement and parenchymal lung disease, with the possible involvement of many extrapulmonary organs.

Among many risk factors, genetic susceptibility seems to be crucial for its development in subjects exposed to unknown infectious or non-infectious environmental factors [1, 2]. This genetic predisposition to sarcoidosis is known to generally concern immunoregulatory genes [3], genes regulating various aspects of angiogenesis (e.g., cell proliferation, apoptosis, migration, morphogenesis) and angiostasis [4, 5] and genes involved in the pathogenesis of lung fibrosis [1, 6–8].

Recently, microRNAs (miRNAs) have attracted attention as being involved in many tissue and cell functions [5, 9]. Mature miRNAs act as small negative regulators of gene expression via reduction of protein translation or degradation of mRNA [10]; they are also involved in feedback interactions as crosstalk between many target genes.

Changes in miRNA expression levels may contribute to various lung pathologies and diseases, including cystic fibrosis, chronic obstructive pulmonary disease, asthma, idiopathic pulmonary fibrosis and lung cancer, by targeting many important signaling pathways, such as the TGF-β/SMAD signaling pathway, or by influencing HIF-1A related genes [9, 11–14].

The number of studies focused on the crosstalk between miRNAs and the genes involved in the pathogenesis of sarcoidosis remains limited [15–17]. Therefore, the present study evaluates the expression pattern of several miRNAs which target genes involved in proangiogenic or angiostatic functions and ECM/EMT remodeling. These are recognized as important processes in the pathogenesis of sarcoidosis, as they are involved in granulomatous formation, fibrosis, and the inhibition of the migration and proliferation of endothelial and epithelial cells [2, 4, 5, 18–21]. The study assesses the expression of miR-let7f, miR-15b, miR-16, miR-20a, miR-27b, miR-128a, miR-130a, miR-192, miR-221, and miR-222: miRNAs which target the VEGF, TFG-B, SMADs, and HIF-1A genes, as well as others under the control of hypoxia. The analyses were performed in BALF cells and peripheral blood (PB) lymphocytes of sarcoidosis patients, using the qPCR method.

Methods

All patients signed an individual consent form. The study was approved by the Medical University of Lodz Ethics Committee (RNN/141/10/KE).

Study group

A total of 94 patients with pulmonary sarcoidosis were recruited into the study. The patients were admitted to the Department of Pneumology and Allergy, Norbert Barlicki memorial University Hospital No. 1 in Lodz, during the years 2010-2014. The diagnosis was made based on current standards [1, 22]. For each patient, a consistent clinical and radiological picture of sarcoidosis was confirmed, with the presence of non-caseating granuloma indicated in tissue biopsy. The diagnosis was documented by EBUS-TBNA, bronchial mucosal biopsy, transbronchial peripheral lung biopsy, mediastinoscopy, or extrathoracic biopsy (skin, peripheral lymph nodes). No biopsy was obligatory only in patients with a typical clinico-radiological picture (bilateral hilar lymph nodes enlargement) and typical BAL results (increased percentage of lymphocytes with CD+4/CD+8 > 3.5). Based on their chest X-ray results, patients were classified into the following radiological subgroups: stage I (hilar lymph node enlargement without signs of parenchymal involvement), stage II (signs of parenchymal involvement in addition to hilar lymph node enlargement), stage III (parenchymal involvement without visible hilar lymph node enlargement) and IV (signs of irreversible extensive lung fibrosis). An independent comparison was performed between patients with acute disease onset (Löfgren syndrome with arthritis, erythema nodosum, elevated body temperature – with at least two symptoms present) and patients with insidious disease onset. The clinical and biological characteristics of the study group are presented Table 1.

Table 1 Clinical and biological characteristics of the study patients
Full size table

The control group consisted of healthy non-smoking persons, referred for bronchoscopy due to chronic cough or the presence of undefined changes on a chest X-ray. If the radiological signs were defined as clinically insignificant changes or artifacts, the patients underwent a thorough examination and were finally diagnosed either with idiopathic cough, or as healthy subjects. A group of 50 subjects was included for gene expression analysis in BALF cells while 20 subjects were included for gene expression analysis in PB lymphocytes.

Bronchoscopy and bronchoalveolar lavage fluid (BALF) collection

Bronchoscopy was performed with a flexible bronchoscope (Pentax, Tokyo, Japan) according to Polish Respiratory Society guidelines [23]. Patients optionally received midanium and atropine before the examination, and 2 % lidocaine was used as a topical anesthetic. BAL fluid (BALF) was collected from medial lobe by instillation and subsequent withdrawal of 4 x 50 ml of 0.9 % NaCl. The fluid recovery was 52.1 ± 1.2 %. The crude BALF was filtered through a gauze to clear the thick mucus and other contaminants, centrifuged, and the pellet was suspended in a phosphate buffer. The total number of non-epithelial cells (total cell count – TCC) was presented as n x 106. Cytospin slides were prepared and stained by May-Grünwald-Giemsa stain. The number of macrophages, lymphocytes, neutrophils, and eosinophils was calculated under a light microscope and presented as % of TCC. After the calculations, the fluid was centrifuged (10 min at 1200 rpm) and the BALF supernatant was suspended in around 350 μl RNAlater RNA Stabilization Reagent (QIAGEN, Hilden, Germany), and frozen (-80 °C) until further RNA isolation procedures could be performed. The percentages of the numbers of macrophages, lymphocytes, neutrophilis and eosinophils obtained from the BALF are shown in Fig. 1.

Fig. 1
figure 1

Graphs presenting the percentage of lymphocytes, neutrophils and eosinophils obtained from BALF

Full size image

Lung function

Spirometry was performed according to the Polish Respiratory Society Guidelines [23] with a computer-based spirometer (Jaeger, Dortmund, Germany). Forced vital capacity (FVC) and forced expiratory volume in 1 s (FEV1) were measured, and the Tiffenau index (FEV1/FVC) was calculated. Lung diffusion capacity for carbon monoxide was measured in patients with lung parenchymal disease only (stage II-IV) by the single breath method, using a Lungtest 1000 SB (MES, Cracow, Poland) according to ATS/ERS standards [24]. The values were corrected for the hemoglobin concentration (DLCOc). All data (except the Tiffeneau index) was expressed as % of predicted value.

Peripheral blood sample collection

Blood was collected into 5 ml tubes containing EDTA. For lymphocyte separation, Histopaque-1077 density gradient cell separation medium (Sigma-Aldrich, Poznan, Poland) was used according to the manufacturers protocol. The supernatant was discarded and cells were resuspended in 350 μl RNAlater RNA Stabilization Reagent before being frozen at -80 °C.

miRNA expression analysis

RNA isolation was performed using a mirVana™ miRNA Isolation Kit (Life Technologies, Carlsbad, CA, USA), according to the manufacturer’s protocol. The quality and quantity of isolated RNA was assessed spectrophotometrically (Eppendorf BioPhotometrTM Plus, Eppendorf, Hamburg, Germany). The purity of total RNA was determined using RNA Nano Chips LabChipplates (ratio of 16S to 18S fraction) and miRNA was evaluated by automated electrophoresis using RNA Pico Chips LabChipplates in an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA).

cDNA was transcribed from 100 ng of total RNA using a TaqMan® MicroRNA Reverse Transcription Kit (Applied Biosystems, Carlsbad, CA) in a total volume of 12 μl according to the manufacturer’s protocol.

The relative expression of miRNAs was assessed by qPCR reactions using Custom TaqMan® Array MicroRNA Cards (Applied Biosystems, Carlsbad, CA) according to the manufacturer’s protocol. The following microRNA probes were used for the study: hsa-let-7f (UGAGGUAGUAGAUUGUAUAGUU), hsa-miR-15b (UAGCAGCACA UCAUGGUUUACA), hsa-miR-16 (UAGCAGCACGUAAAUAUUGGCG), hsa-miR-20a (ACUGCAUUAUGAGCACUUAAAG), hsa-miR-27b (AGAGCUUAGCUGAUUGGUGA AC), hsa-miR-128 (UCACAGUGAACCGGUCUCUUU), hsa-miR-130a (CAGUGCAAUG UUAAAAGGGCAU), hsa-miR-192 (CUGCCAAUUCCAUAGGUCACAG), hsa-miR-221 (ACCUGGCAUACAAUGUAGAUUU), hsa-miR-22(CUCAGUAGCCAGUGUAGAUCCU).

The PCR mixture contained 0.90 μl RT Product, 56.25 μl TaqMan® Universal Master Mix II, No AmpErase® UNG (2X) and 55.35 μl Nuclease-free water (Applied Biosystems, Carlsbad, CA). TaqMan Array card was centrifuged according to manufactured protocol and placed in the Applied Biosystems 7900HT Fast Real-Time PCR System (Applied Biosystems, Carlsbad, CA) according to manufactured protocol. The expression levels (RQ values) of the studied miRNA were calculated using the delta delta CT method, with the adjustment to the level of 001973 U6 snRNA expression (endogenous control) and to the expression of the calibrator (Human Lung Total RNA Ambion®), for which RQ value was equal to 1.

Statistical analysis

The Kruskal–Wallis test, Mann–Whitney U-test, Neuman–Keuls’ multiple comparison test, and Spearman’s rank correlation were used to assess the correlations between relative miRNA expression and sarcoidosis in the patient groups classified based on radiological examination (stages I vs. II-IV), acute vs. chronic disease onset, spirometric parameters, DLCO, serum Ca2+ concentration, Ca2+ loss in urine over 24 h, the percentage of lymphocytes in BAL (BALF%), the phenotype of BAL immune cells (CD4+/CD8+), age and sex of patients. P = 0.05 was regarded as the level of statistical significance (StatSoft, Cracow, Poland).

Results

Relative expression levels of the studied miRNAs in sarcoidosis patients vs. controls

BALF cells

The Mann-Whitney U-test revealed significantly greater miRNA expression in sarcoidosis patients miR-192 (P = 0.00004, Mann-Whitney U- test) and for miR-221 (P = 0.0002, Mann-Whitney U- test) (see Fig. 2). Conversely, significantly greater miR-15b expression was noted in the control group (P = 0.04, Mann-Whitney U- test).

Fig. 2
figure 2

Box and whisker plots, presenting a) miR-192 b) miR-221 c) miR-15b expression mean values in BALF cells of patients and controls (P = 0.00004, Mann-Whitney U- test)

Full size image

No statistically significant differences were observed between patients and controls for the remaining analyzed miRNAs (P > 0.05, Mann-Whitney U- test) (see Additional file 1: Figure S1).

Peripheral blood (PB) lymphocytes

Statistically significant differences (Mann-Whitney U- test) between patients and controls were observed for miR-15b (P = 0.00009, Mann-Whitney U- test), miR-27b (P = 0.01, Mann-Whitney U- test), miR-192 (P = 0.005, Mann-Whitney U- test), miR-221 (P, Mann-Whitney U- test), and miR-222 (P = 0.0061, Mann-Whitney U- test) with higher miRNA expression levels in sarcoidosis patients. Statistically significant differences (Mann-Whitney U- test) were also found for miR-let7f (P = 0.00004, Mann-Whitney U- test), and miR-130a (P = 0.02, Mann-Whitney U- test), with higher miRNA expression levels in the control group (see Fig. 3).

Fig. 3
figure 3

Box and whisker plots, presenting a) miR-15b b)miR-27b c) miR-192 d) miR-221 e) miR-222 f) miR-letf7 g) miR-130a expression (mean values) in the PB lymphocytes of patients and controls (P = 0.0002, Mann-Whitney U- test)

Full size image

No statistically significant differences were found between patients and controls for the rest of analyzed miRNAs (P > 0.05, Mann-Whitney U- test) (see Additional file 2: Figure S2).

Relative expression of the studied miRNAs in sarcoidosis patients according to clinical and radiological classification

BALF cells

In patients with radiological stage I, miR-192 demonstrated the greatest mean expression of the evaluated miRNAs (mean RQ = 12.03) and miR-130a the lowest (mean RQ = 0.907). In patients with radiological stages II-IV, the highest mean expression was associated with miR-let7f (mean RQ = 7.92) and the lowest with miR-130a (mean RQ = 0.343). In the acute form of disease, as with radiological stage I, miR-192 demonstrated the greatest mean expression (mean RQ = 12.099) and miR-130a the lowest (mean RQ = 0.637). Finally, in the insidious form of the disease, miR-222 demonstrated the highest mean expression (mean RQ = 8.600) and miR-130a the lowest (mean RQ = 0.615) (see Fig. 4).

Fig. 4
figure 4

Relative expression (mean RQ values) of the studied miRNAs in sarcoidosis patients according to the clinical and radiological classification in BALF cells

Full size image

The expression of miR-130a was decreased in the majority of studied samples (86–96 %), and the expression of miR-192 was increased in the majority of studied samples (74-91 %), regardless the of the clinical disease classification of the group.

Significantly greater miRNA expression was observed in patients without parenchymal involvement (stages I) than those with involvement (Stages II-IV) for miR-27b (P = 0.04), miR-192 (P = 0.02, Mann-Whitney U- test) and for miR-221 (P = 0.03, Mann-Whitney U- test) (see Additional file 3: Figure S3).

Additionally, significantly higher miRNA expression was found in patients with the acute form of disease for miR-27b (P = 0.002, Mann-Whitney U- test), miR-192 (P = 0.01, Mann-Whitney U- test), and miR-221 (P = 0.02, Mann-Whitney U- test) (see Additional file 4: Figure S4).

No statistically significant differences were observed in the expression of the studied miRNAs in BALF cells for patients with normal and abnormal spirometry results (P > 0.05, Mann-Whitney U- test) (see Additional file 5: Figure S5). However, the expression of miRNA-16 and miRNA-20a in BALF cells was found to be significantly elevated in the subgroup of patients with a restrictive spirometric pattern (P = 0.009 and P = 0.04, respectively, Mann-Whitney U- test, see Fig. 5 and Additional file 6: Figure S6).

Fig. 5
figure 5

Box-and-whisker plots representing the expression of a) miR-16 in BALF cells (P = 0.009; Mann-Whitney U- test) in patients with normal spirometry and with lung volume restriction, b) miR-20a in BALF cells (P = 0.04; Mann-Whitney U- test) in patients with normal spirometry and in patients with lung volume restriction

Full size image

PB lymphocytes

In patients with radiological stage I, miR-15b was found to have the highest mean expression of the observed miRNAs (mean RQ = 28.52) and miR-20a the lowest (mean RQ = 1.057). In patients with radiological stages II-IV, miR-221 demonstrated the highest mean expression (mean RQ = 11.5840) and miR-130a the lowest (mean RQ = 1.081).

In the acute onset group, miR-15b was the highest (mean RQ = 26.61) and miR-27b the lowest (mean RQ = 1.478), while in the insidious onset group, miR-221 demonstrated the greatest expression (mean RQ = 14.267) and miR-130a the lowest (mean RQ =0.670) (see Fig. 6).

Fig. 6
figure 6

Relative expression (mean RQ values) of the studied miRNAs in sarcoidosis patients according to the clinical and radiological classification in PB lymphocytes

Full size image

The expression of miR-221 was elevated in the majority of studied samples (89–100 %), while that of miR-130a was decreased (50-78 %), regardless of the clinical disease classification.

Significantly higher expression of miRNA-16 (P = 0.05, Mann-Whitney U- test) and miR-15b (P = 0.001, Mann-Whitney U- test) was observed in patients without parenchymal involvement (see Additional file 7: Figure S7), while miRNA-130a (P = 0.04, Mann-Whitney U – test) and miR-15b (P = 0.004, Mann-Whitney U- test) demonstrated higher expression in patients with the acute form of disease (see Additional file 8: Figure S8). Additionally, significantly increased expression of miRNA-let7f was found in the subgroup of patients with obturation (P = 0.001, Mann-Whitney U- test) (see Additional file 9: Figure S9).

Relative expression levels of the studied miRNAs in relation to patient gender, age, selected biochemical markers, immunological phenotype and lung function test parameters

BALF cells

Several correlations were found between miRNA expression levels and lung function test results, as well as with selected laboratory markers in sarcoidosis patients (see Table 2). No statistically significant correlations were found between miRNA expression and gender (P > 0.05, Mann-Whitney U-test).

Table 2 Significant correlations found in the BALF cells and PB lymphocytes of sarcoidosis patients
Full size table

PB lymphocytes

Significantly higher expression of miRNA-let7f was observed in men with chronic form of disease (P = 0.03, U Mann-Whitney test). Correlations between miRNA expression level and age, immunological phenotype, lung function test parameters, and selected laboratory markers in sarcoidosis patients are shown in Table 2.

Relative expression levels of the studied miRNAs: BALF cells vs PB lymphocytes in sarcoidosis patients

The majority of studied BALF cell samples demonstrated higher expression of miR-let-7f, miR-16, miR-192, miR-222, but lower expression of miR-15b, miR-20a, miR-27b, miR-128, miR-130, miR-221 compared to PB lymphocytes, regardless of the clinical disease classification groups.

Only two of the studied miRNAs, miR-16 and miR-128a, demonstrated overlapping patterns in BALF cells and PB lymphocytes (see Fig. 7)

Fig. 7
figure 7

The expression pattern of miRNAs (mean RQ value) in BALF cells and PB lymphocytes in the sarcoidosis patient group

Full size image

However, statistically significant differences between BALF cells and PB lymphocytes were only observed for three of the non-overlapping miRNAs: for miRNA-15b (P = 0.00001, Mann-Whitney U- test) and miRNA-130a (P = 0.000012, Mann-Whitney U- test), with higher expression in PB lymphocytes, and for miRNA-let7f, with higher expression levels in BALF cells (P = 0.000003, Mann-Whitney U- test).

Discussion

Recent studies have found non-coding regulatory RNAs to play significant roles in the direct regulation of many genes and signaling pathways, as well as other miRNAs [5, 9, 24–26]. However, their biological role is not fully understood. Previous studies regarding miRNA interactions in the lung have mainly focused on the control of innate and adaptive lung inflammation [17, 27, 28], pulmonary fibrosis (IPF, idiopathic pulmonary fibrosis) and the regulation of extracellular matrix (ECM) remodeling [16, 29, 30]. Few studies examine miRNA expression in the course of sarcoidosis [17, 31], and the results which have been obtained are controversial. Therefore, our study offers an innovative and valuable insight into the area, not only because it addresses miRNA expression levels in the peripheral circulation, but it also attempts to determine their value in the clinical classification and monitoring of sarcoidosis patients.

Our findings are difficult to compare with those of other studies, as they use a range of biological materials to investigate the studied miRNAs. In addition, the examined miRNAs may vary with regard to their origin (free circulating miRNA or exosomal miRNA), profile and concentration. Studies on pulmonary sarcoidosis have used different classes of miRNAs and diverse biological material for miRNA isolation (BALF cells, peripheral blood mononuclear cells and serum) [17, 31]. The most recent studies, by Qazi et al [32] and Ho et al. [33], suggest that exosomes may be suitable for miRNA assessment in sarcoidosis. They note the presence of an elevated level of exosomes in BALF from sarcoidosis patients and suggest that exosomes may contribute to the inflammatory state of sarcoidosis via stimulation of cytokine production from peripheral mononuclear cells. On the other hand, exosome production is common in various cells and the mechanisms by which they control the inflammatory process are not fully understood, as well as the biological functions of miRNAs derived from the exosomes.

As miRNAs play a role in cellular communication, the most beneficial approach would be to assess their role in a particular signaling pathway which is important in sarcoidosis, irrespective of their varied origins. It should be emphasized that microarray-based studies have confirmed that alterations in miRNA expression influence the pathogenesis of sarcoidosis, as well as its relationship with the TGF-β and WNT/β-catenin signaling pathways [31]. Therefore, for the purposes of the present study, several novel miRNAs were chosen which were highly complementary to genes involved in TGF-β/Smad (regulation of fibrosis) and HIF-1A/VEGF/ING-4 (angiogenesis/angiostasis homeostasis) signaling. This choice was based on their being awarded high target scores according to the miRDB database.

The present study focuses on these two pathways because an important pathogenic component of sarcoidosis is pulmonary fibrosis, a condition associated with chronic inflammation, whose symptoms are known to be impaired wound healing and angiogenesis. Our earlier study assesses the expression of genes involved in these two pathways, i.e. the regulation of fibrosis and angiogenetic/angiostatic homeostasis, in patients with sarcoidosis [6, 34]. Our findings reveal statistically significant correlations between the expression of some of the studied miRNAs and the expression of genes which had been studied earlier [6, 34].

The present study compares the expression of selected miRNAs in PB lymphocytes and BALF cells to analyze their usefulness as specific systemic and/or local sarcoidosis markers. Although significant differences were found between patients and controls regarding the expression of several examined miRNAs in PB lymphocytes (miR-15b, miR-27b, miR-192, miR-221, miR-222, miR-130a, miR-let-7f), only three miRNAs were found to differ in BALF cells (miR-15b, miR-192, and miR-221). Our findings suggest that the above-mentioned miRNAs may be vital molecular markers for the disease, and merit further examination. It should be pointed that the biomarkers found in BALF cells appear to have high diagnostic and/or prognostic value, and may also confirm the specificity of changes in miRNA expression occurring in response to the lung microenvironment modifications which take place during sarcoidosis. We believe that of the two examined cell types, miRNA analysis may have the highest diagnostic value when performed in BALF cells. Moreover, it is interesting that the expression patterns of two studied miRNAs, miR-16 and miR-128a, overlapped in the BALF cells and PB lymphocytes of sarcoidosis patients, and their expression was higher in patients than in controls (see Fig. 5).

It is commonly known that increases in the numbers of T-cells induced by chemokines secreted by macrophages and epithelial cells, as well as shifts in the Th1/Th2 balance, play crucial roles in the course of sarcoidosis. A recent report by Barczyk et al. [35] demonstrates that CCL4 is highly expressed in T-cells (especially in CD8+, less so in CD4+) and may serve to recruit new T cells to the site of inflammation. As a prolonged state of macrophage activation and CD4+ T-lymphocyte accumulation in the presence of chronic antigenic stimulation is known to lead to inflammation in the early stage of sarcoidosis, “early sarcoid reaction” [36] and granuloma formation, an earlier study by Jazwa et al evaluated the expression of miR-15b and miR-16 as important regulators of CD4+ T cells [17]. Based on clinical and radiological classifications, significant overexpression of miR-15b and miR-16 in PB lymphocytes was observed in patients with stage I disease and in patients with the acute form of the disease. A negative correlation was also found between miR-16 expression and a lung function parameter (FEV1/FVC), and significantly higher miR-16 expression was observed in patients with lung volume restriction. The results suggest that miR-15b and miR-16 may be potentially useful prognostic markers in relation to disease progression and worse prognosis. Therefore, the influence of both miRNAs on VEGF, a target gene in angiogenesis and an inflammatory cascade regulator in sarcoidosis, seems to be pivotal in sarcoidosis pathogenesis, regarding the hypoxia-angiogenesis process and chronic inflammation in early stages of the disease, as emphasized in previous studies [37, 38]. Our previous study also confirmed the significant role played by the up-regulation of VEGF mRNA in sarcoidosis [34].

However, our results regarding miR-15b and miR-16 are not in line with those of Jazwa et al. [17], who did not report any statistically significant differences between patients and healthy controls with regard to the expression of these miRNAs.

Our present study also evaluates the expression of miR-let7f. This miRNA belongs to the class that targets TGF-β [39], a key cytokine responsible for lung fibrosis and epithelial-mesenchymal transition (EMT) [40–42]. It has been confirmed that in the lung tissue, miR-let7 influences the epithelial expression of mesenchymal markers and deposition of collagen without additional injury, which may suggest that the dysregulation of miR-let7 can lead to the overexpression of crucial pro-fibrotic factors and may be involved in the pathogenesis of pulmonary fibrosis [13, 43, 44]. Our findings reveal increased miR-let7f expression in BALF cells in patients, and a positive correlation between miR-let7f expression and Ca2 + U in the PB lymphocytes of the insidious onset group of patients. According to RTG staging (I, II-IV), the greatest mean expression of miR-let7f was observed during RTG stages II-IV. As calcium metabolism dysregulation occurs in the course of sarcoidosis and the presence of hypercalciuria and hypercalcemia in 20-50 % of patients [45], miR-let7f up-regulation may be a valuable marker in the prediction of calcium metabolism dysregulation during sarcoidosis.

To the best of our knowledge, this is the first study to examine the expression of miR-20a and miR-27b in sarcoidosis. Nevertheless, it has already been recognized that the miR-17 ~ 92 cluster, encoding six different miRNAs, is decreased in patients with IPF [46]. In the case of miR-20a and miR-27b, both are recognized as important regulators of the TGF-β/Smad pathway and fibrosis in the lung, with miR-20a additionally targeting VEGF [47]. In the present study, miR-27b expression was found to be increased in the PB lymphocytes of sarcoidosis patients, but this level was significantly lower in the advanced stages of the disease (RTG stage II-IV) than in RTG stage I. Our findings may suggest that miRNA-27b inhibition takes place during disease progression. This miRNA has been shown to play a regulatory role in pro-fibrotic changes in pulmonary cells in response to TGF-β signaling [48]. In addition, reductions in miR-27b expression were also found to be associated with TGF-β stimulation [48].

In the present study, interesting results were obtained in the case of miR 192, which showed significant up-regulation in both BALF cells and PB lymphocytes. However, its expression was found to be significantly reduced during disease progression (in stage II-IV), which is consistent with the results obtained by Hong et al. [49], who confirmed the reduction of miR 192 expression via TGF-β/Smad3 interaction in a study on chronic renal fibrosis [49].

In conclusion, the data presented in this study indicate that selected miRNA may enable patients with sarcoidosis to be distinguished from potentially healthy controls, based on the different expression profiles of the two groups (see Fig. 8).

Fig. 8
figure 8

The differences in miRNA expression profiles (in BALF cells and PB lymphocytes) between controls and patients with sarcoidosis

Full size image

Conclusion

The results of this study indicate that the expression profiles of selected miRNA allow patients with sarcoidosis to be distinguished from potentially healthy controls. Moreover, it should be emphasized that some of these miRNAs may be themselves recognized as markers of disease progression, due to their elevated expression in patients with lung volume restriction or obturation. The deregulated miRNA expression observed in this study may play a vital role in the pathogenesis of sarcoidosis. In particular, the observed significant up-regulation of the miRNAs involved in angiogenesis (miR-192, miR-221, miR-222) or EMT control (miR-192, miR-221) may indicate the presence of a significant relationship between the two processes, and may suggest that they may play a role in the pathogenesis of the disease. The obtained results highlight distinct patterns of miRNA expression in BALF cells and PB lymphocytes in sarcoidosis patients, indicating the importance of the lung microenvironment in the modulation of their expression. However, further studies are needed on larger groups of patients to develop our findings concerning these selected miRNAs and gain a greater understanding of their importance as clinical markers for disease monitoring and making therapeutic decisions.

Abbreviations

BAL:

Bronchoalveolar lavage

BALF:

BAL fluid

DLCO:

Diffusion capacity for carbon monoxide

EBUS-TBNA:

Endobronchial ultrasonography guided transbronchial needle aspiration

FEV1:

Forced expiratory volume in 1st second of expiration

FVC:

Forced vital capacity

LFT:

Lung function test

miRNA:

microRNA

References

  1. Ianuzzi MC, Rybicki BA, Teirstein AS. Sarcoidosis. N Engl J Med. 2007;357:2153–65.

    Article  Google Scholar 

  2. Lazarus A. Sarcoidosis: epidemiology, etiology, pathogenesis, and genetics. Dis Mon. 2009;55:649–60.

    Article  PubMed  Google Scholar 

  3. Spagnolo P, Grunewald J. Recent advances in the genetics of sarcoidosis. J Med Genet. 2013;50:290–7.

    Article  CAS  PubMed  Google Scholar 

  4. Tzouvelekis A, Ntolios P, Karameris A, Koutsopoulos A, Boglou P, Koulelidis A, et al. Expression of hypoxia-inducible factor (HIF)-1a-vascular endothelial growth factor (VEGF)-inhibitory growth factor (ING)-4- axis in sarcoidosis patients. BMC Res Notes. 2013;5:654.

    Article  Google Scholar 

  5. Wu F, Yang Z, Li G. Role of specific microRNAs for endothelial function and angiogenesis. Biochem Biophys Res Commun. 2009;386:549–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Piotrowski WJ, Kiszałkiewicz J, Pastuszak-Lewandoska D, Antczak A, Górski P, Migdalska-Sęk M, et al. TGF-β and SMADs mRNA Expression Levels in Pulmonary Sarcoidosis. Adv Exp Med Biol. 2015;852:59–69.

    Article  CAS  PubMed  Google Scholar 

  7. Salez F, Gosset P, Copin MC, Lamblin Degros C, Tonnel A.B, Wallaert B. Transforming growth factor-β1 in sarcoidosis. Eur Respir J. 1998;12:913–9.

    Article  CAS  PubMed  Google Scholar 

  8. Szlubowski A, Soja J, Grzanka P, Tomaszewska R, Papla B, Kuzdzał J, et al. TGF-beta1 in bronchoalveolar lavage fluid in diffuse parenchymal lung diseases and high-resolution computed tomography score. Pol Arch Med Wewn. 2010;120:270–5.

    CAS  PubMed  Google Scholar 

  9. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;23:281–97.

    Article  Google Scholar 

  10. He L, Hannon GJ. MicroRNAs: small RNAs with a big role in gene regulation. Nat Rev Genet. 2004;5:522–31.

    Article  CAS  PubMed  Google Scholar 

  11. Butz H, Rácz K, Hunyady L, Patócs A. Crosstalk between TGF-β signaling and the microRNA machinery. Trends Pharmacol Sci. 2012;33:382–93.

    Article  CAS  PubMed  Google Scholar 

  12. Ming-Kuen Tang P, Lan HY. MicroRNAs in TGF-b/Smad-mediated Tissue Fibrosis. Curr Pathobiol Rep. 2014;2:235–43.

    Article  Google Scholar 

  13. Pandit KV, Corcoran D, Yousef H, Yarlagadda M, Tzouvelekis A, Gibson KF, et al. Inhibition and role of let-7d in idiopathic pulmonary fibrosis. Am J Respir Crit Care Med. 2010;15:220–9.

    Article  Google Scholar 

  14. Tomankova T, Petrek M, Kriegova E. Involvement of microRNAs in physiological and pathological processes in the lung. Respir Res. 2010;23:159.

    Article  Google Scholar 

  15. Maertzdorf J, Weiner J, Mollenkopf, TBornotTB Network, Bauer T, Prasse A, et al. Common patterns and disease-related signatures in tuberculosis and sarcoidosis. PNAS. 2012; 20: 7853-7858.

  16. Cox TR, Erler JT. Remodeling and homeostasis of the extracellular matrix: implications for fibrotic diseases and cancer. Dis Model Mech. 2011;4:165–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Jazwa A, Kasper L, Bak M, Sobczak M, Szade K, Jozkowicz A, et al. Differential Inflammatory MicroRNA and Cytokine Expression in Pulmonary Sarcoidosis. Arch Immunol Ther Exp. 2015;63:139–46.

    Article  CAS  Google Scholar 

  18. Antoniou KM, Tzouvelekis A, Alexandrakis MG, Sfiridaki K, Tsiligianni I, Rachiotis G, et al. Different angiogenic activity in pulmonary sarcoidosisand idiopathic pulmonary fibrosis. Chest. 2006;130:982–8.

    Article  PubMed  Google Scholar 

  19. Fish JE, Santoro MM, Morton SU, Yu S, Yeh RF, Wythe JD, et al. miR-126 regulates angiogenic signaling and vascular integrity. Dev Cell. 2008;15:272–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Poliseno L, Tuccoli A, Mariani L, Evangelista M, Citti L, Woods K, et al. MicroRNAs modulate the angiogenic properties of HUVECs. Blood. 2006;108:3068–71.

    Article  CAS  PubMed  Google Scholar 

  21. Tzouvelekis A, Anevlavis S, Bouros D. Angiogenesis in interstitial lung diseases: a pathogenetic hallmark or a bystander? Respir Res. 2006;25:82.

    Article  Google Scholar 

  22. American Thoracic Society: European Respiratory Society. World association of sarcoidosis and other granulomatous disorders. Statement on sarcoidosis. Am J Respir Crit Care Med. 1999;160:736–55.

    Article  Google Scholar 

  23. Polish Society of Resiratory Diseases. The recommendations of the Polish Society Pneumological on spirometry. Pneumonol Alergol Pol. 2006;74:(Suppl 1).

  24. European Respiratory Society. Standardized lung function testing. Official statement of the European Respiratory Society. Eur Respir J Suppl. 1993;16:1–100.

    Google Scholar 

  25. Asirvatham AJ, Gregorie CJ, Hu Z, Magner WJ, Tomasi TB. MicroRNA targets in immune genes and the Dicer/Argonaute and ARE machinery components. Mol Immunol. 2008;45:1995–2006.

    Article  CAS  PubMed  Google Scholar 

  26. Xu P, Guo M, Hay BA. MicroRNAs and the regulation of cell death. Trends Genet. 2004;20:617–24.

    Article  CAS  PubMed  Google Scholar 

  27. Ha TY. The role of microRNAs in regulatory T cells and in the immune response. Immune Netw. 2011;11:11–41.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Perry MM, Moschos SA, Williams AE, Shepherd NJ, Larner-Svensson HM, Lindsay MA. Rapid changes in microRNA-146a expression negatively regulate the IL-1beta-induced inflammatory response in human lung alveolar epithelial cells. J Immunol. 2008;180:5689–98.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. 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;2:1589–97.

    Article  Google Scholar 

  30. Williams AE, Moschos SA, Perry MM, Barnes PJ, Lindsay MA. Maternally imprinted microRNAs are differentially expressed during mouse and human lung development. Dev Dyn. 2007;2362:572–80.

    Article  Google Scholar 

  31. Crouser ED, Julian MW, Crawford M, Shao G, Yu L, Planck SR, et al. Differential expression of microRNA and predicted targets in pulmonary sarcoidosis. Biochem Biophys Res Commun. 2012;13:886–91.

    Article  Google Scholar 

  32. Qazi K, Torregrosa Parades P, Dahlberg B, Grunewald J, Eklund A, Gabrielsson S. Proinflammatory exosomes in bronchoalveolar lavage fluid of patients with sarcoidosis. Thorax. 2010;11:1016–24.

    Article  Google Scholar 

  33. Ho LP. Exosomes in lungs of patients with sarcoidosis: a contributor to immune pathogenesis or just another by-product of heightened immune activity? Thorax. 2010;11:947–8.

    Article  Google Scholar 

  34. Piotrowski WJ, Kiszałkiewicz J, Pastuszak-Lewandoska D, Górski P, Antczak A, Migdalska-Sęk, M et al. Expression of HIF-1A/VEGF/ING-4 axis in pulmonary sarcoidosis. Adv Exp Med Biol. 2015;866:61–9.

    Article  CAS  PubMed  Google Scholar 

  35. Barczyk A, Pierzchała E, Caramori G, Sozańska E. Increased expression of CCL4/MIP-1β in CD8+ cells and CD4+ cells in sarcoidosis. Pharmacol. 2014;27:185–93.

    CAS  Google Scholar 

  36. Agostini C, Meneghin A, Semenzato G. T-lymphocytes and cytokines in sarcoidosis. Curr Opin Pulm Med. 2002;8:435–40.

    Article  PubMed  Google Scholar 

  37. Chamorro-Jorganes A, Araldi E, Penalva LO, Sandhu D, Fernández-Hernando C, Suárez Y. MicroRNA-16 and microRNA-424 regulate cell-autonomous angiogenic functions in endothelial cells via targeting vascular endothelial growth factor receptor-2 and fibroblast growth factor receptor-1. Arterioscler Thromb Vasc Biol. 2011;3:2595–606.

    Article  Google Scholar 

  38. Tolnay E, Kuhnen C, Voss B, Wiethege T, Müller KM. Expression and localization of vascular endothelial growth factor and its receptor flt in pulmonary sarcoidosis. Virchows Arch. 1998;432:61–5.

    Article  CAS  PubMed  Google Scholar 

  39. Xie T, Liang J, Guo R, Liu N, Noble PW, Jiang D. Comprehensive microRNA analysis in bleomycin-induced pulmonary fibrosis identifies multiple sites of molecular regulation. Physiol Genomics. 2011;139:479–87.

    Article  Google Scholar 

  40. Khalil N, O’Connor RN, Flanders KC, Unruh H. TGF-beta 1, but not TGF-beta 2 or TGF-beta 3, is differentially present in epithelial cells of advanced pulmonary fibrosis: an immunohistochemical study. Am J Respir Cell Mol Biol. 1996;14:131–8.

    Article  CAS  PubMed  Google Scholar 

  41. Patterson KC, Hogarth K, Husain AN, Sperling AI, Niewold TB. The clinical and immunologic features of pulmonary fibrosis in sarcoidosis. Transl Res. 2012;160:321–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Shi K, Jiang J, Ma T, Xie J, Duan L, Chen R, et al. Pathogenesis pathways of idiopathic pulmonary fibrosis in bleomycin-induced lung injury model in mice. Respir Physiol Neurobiol. 2014;1:113–7.

    Article  Google Scholar 

  43. Levanen B, Wheelock AM, Eklund A, Grunewald J, Nord M. Increased pulmonary Wnt (wingless/integrated)-signaling in patients with sarcoidosis. Respir Med. 2011;105:282–9.

    Article  PubMed  Google Scholar 

  44. Pandit KV, Milosevic J, Kaminski N. MicroRNAs in idiopathic pulmonary fibrosis. Transl Res. 2011;57:191–9.

    Article  Google Scholar 

  45. Kempisty A, Kuś J. Hypercalcemia and renal failure in the course of sarcoidosis-case report. Pneumonol Alergol Pol. 2012;80:570–5.

    PubMed  Google Scholar 

  46. Dakhlallah D, Batte K, Wang Y, Cantemir-Stone CZ, Yan P, Nuovo G, et al. Epigenetic regulation of miR-17 ~ 92 contributes to the pathogenesis of pulmonary fibrosis. Am J Respir Crit Care Med. 2013;187:397–405.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Suárez Y, Sessa WC. MicroRNAs as novel regulators of angiogenesis. Circ Res. 2009;27:442–54.

    Article  Google Scholar 

  48. Graham JR, Williams CM, Yang Z. MicroRNA-27b targets gremlin 1 to modulate fibrotic responses in pulmonary cells. J Cell Biochem. 2014;115:1539–48.

    Article  CAS  PubMed  Google Scholar 

  49. Hong JP, Li XM, Li MX, Zheng FL. VEGF suppresses epithelial-mesenchymal transition by inhibiting the expression of Smad3 and miR-192, a Smad3-dependent microRNA. Int J Mol Med. 2013;31:1436–42.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgments

This work was funded by a grant from the National Science Centre POLAND (grant no. 2011/01/B/NZ5/04239).

Author information

Authors and Affiliations

  1. Department of Molecular Bases of Medicine, 1st Chair of Internal Diseases, Medical University of Lodz, Pomorska 251, 92-213, Lodz, Poland

    Justyna Kiszałkiewicz, Dorota Pastuszak-Lewandoska, Daria Domańska-Senderowska, Monika Migdalska-Sęk, Karolina H. Czarnecka, Ewa Nawrot & Ewa Brzeziańska-Lasota

  2. Department of Pneumonology and Allergy, 1st Chair of Internal Diseases, Medical University of Lodz, Kopcińskiego 22, 90-153, Lodz, Poland

    Wojciech J. Piotrowski, Paweł Górski & Witold Górski

  3. Department of General and Oncological Pulmonology, 1st Chair of Internal Diseases, Medical University of Lodz, Kopcińskiego 22, 90-153, Lodz, Poland

    Adam Antczak

Authors
  1. Justyna Kiszałkiewicz
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Wojciech J. Piotrowski
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Dorota Pastuszak-Lewandoska
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Paweł Górski
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Adam Antczak
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Witold Górski
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Daria Domańska-Senderowska
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Monika Migdalska-Sęk
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Karolina H. Czarnecka
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Ewa Nawrot
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Ewa Brzeziańska-Lasota
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Ewa Brzeziańska-Lasota.

Additional information

Competing interests

Authors declare no conflicts of interest in relation to this article.

Authors’ contributions

EB-L conceived the idea for the study and prepared the final version of the article. JK performed the experiment and contributed the design of the research, DP-L, MM-S, KHC, EN completed the experimental work, D.D performed the statistical analysis, WJP, PG, AA, WG coordinated the collection and selection of patients. All authors were involved in data collection. All authors read and approved the final manuscript.

Additional files

Additional file 1:Figure S1.

Box and whisker plots, presenting a) miR-16 b) miR-20a c) miR-27b d) miR-128a e) miR-130a f) miR-222 g) miR-let7f expression mean values in BALF cells of patients and controls. (JPEG 3074 kb)

Additional file 2:Figure S2.

Box and whisker plots, presenting a) miR-16 b) miR-20a c) miR-128a expression mean values in PB lymphocytes of patients and controls. (JPEG 1357 kb)

Additional file 3:Figure S3.

Box and whisker plots representing the expression of a) miR-15b b) miR-16 c) d) miR-27b, miR-20a e) miR128 f) miR-130a g) miR-192 h) miR-221 i) miR-222 j) miR-let7f  in BALF cells in patients in radiological stags I vs. II-IV (JPEG 4235 kb)

Additional file 4:Figure S4.

Box and whisker plots representing the expression of a) miR-15b b) miR-16 c) miR-20a d) miR-27b e) miR128 f) miR-130a g) miR-192 h) miR-221 i) miR-222 j) miR-let7f  in BALF cells in patients with the acute vs. insidious form of disease. (JPEG 4285 kb)

Additional file 5:Figure S5.

Box and whisker plots representing the expression of a) miR-15b b) miR-16 c) miR-20a d) miR-27b e) miR128 f) miR-130a  g) miR-192 h) miR-221 i) miR-221 j) miR-let7f  in BALF cells in patients with normal vs.abnormal spirometry (JPEG 4306 kb)

Additional file 6:Figure S6.

Box and whisker plots representing the expression of a) miR-15b b) miR-27b c) miR-128a  d) miR-130 e) miR-192 f) miR-221 g) miR-222 h) miR-let7f  in BALF cells in subgroup of patients with a restrictive spirometric pattern vs. normal spirometry (JPEG 3712 kb)

Additional file 7:Figure S7.

Box and whisker plots representing the expression of a) miR-15b b) miR-16  c) miR-20a d) miR-27b e) miR-128a f) miR-130a g) miR-192 h) miR-221 i) miR-222 j) miR-let7f in PB lymphocytes in patients in radiological stages I vs. II-IV (JPEG 4059 kb)

Additional file 8:Figure S8.

Box and whisker plots representing the expression of a) miR-15b b) miR-16 c) miR-20a d) miR-27b e) miR-128a f) miR-130a  g) miR-192 h) miR-221 i) miR-222 j) miR-let7f in PB lymphocytes in patients with acute vs. insidious form of disease. (JPEG 4053 kb)

Additional file 9:Figure S9.

Box and whisker plots representing the expression of a) miR-15b b) miR-16 c) miR-20a d) miR-27b e) miR128 f) miR-130a  g) miR-192 h) miR-221 i miR-222 j) miR-let7f  in PB lymphocytes in patients with obturation vs. abnormal spirometry. (JPEG 4159 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kiszałkiewicz, J., Piotrowski, W.J., Pastuszak-Lewandoska, D. et al. Altered miRNA expression in pulmonary sarcoidosis. BMC Med Genet 17, 2 (2016). https://0-doi-org.brum.beds.ac.uk/10.1186/s12881-016-0266-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12881-016-0266-6

Keywords

BMC Medical Genetics

ISSN: 1471-2350

Contact us