Abstract

Background. Bullous pemphigoid (BP) is a common T helper 2- (Th2-) dominated autoimmune blistering skin disease with significant mortality. MicroRNAs (miRNAs), which are endogenous noncoding RNA molecules, have been reported to be potential biomarkers for some autoimmune diseases; however, to date, there exist no reports on serum expression profiles of miRNAs in BP patients. Methods. A RNA quantitative PCR- (qPCR-) based array was conducted on sera from 20 active BP patients and 20 healthy controls for screening of miRNAs. Significantly dysregulated miRNAs were validated with use of qPCR as performed on sera samples of 45 active BP patients and 60 healthy controls. Serum CCL17, anti-BP180, and anti-BP230 levels were measured with use of ELISA. Results. Relative baseline expression levels of serum miR-1291 were significantly upregulated in the 45 BP patients as compared with the 60 healthy controls () and significantly decreased in the disease control stage (, ). In addition, these baseline miR-1291 levels showed a significant positive correlation with the baseline levels of serum CCL17 () and anti-BP180 (, ). Like that observed for miR-1291, baseline levels of serum CCL17 were also significantly elevated in the 45 BP patients compared with the 60 healthy controls () and significantly decreased in the disease control stage (, ). However, for anti-BP180, baseline serum levels were significantly elevated in only 38 of the 45 BP patients and significantly decreased in the disease control stage (, ). Conclusions. Relative expression levels of serum miR-1291 can reflect disease activity of BP. miR-1291 may function as an important new serum biomarker for BP.

1. Introduction

Bullous pemphigoid (BP) is a T helper 2- (Th2-) dominated autoimmune subepidermal blistering skin disease resulting from circulating autoantibodies directed against structural components of hemidesmosomes. These autoantibodies mainly consist of bullous pemphigoid antigen 2 (BPAG2, BP180) and bullous pemphigoid antigen 1 (BP230, BPAG1) [1]. Serum levels of anti-BP180 antibody have been reported to be elevated in 72-93% of BP patients and are correlated with disease activity of BP, while anti-BP230 antibody levels are detected in only 57-63% of BP patients and do not reflect BP disease activity [2]. Th2-type inflammation plays an important role in the pathogenesis of BP, and CCL17 is a well-recognized Th2 chemokine. Serum CCL17 levels have been reported to be significantly elevated in active BP patients and reflect disease activity of this condition [35]. In addition to BP, serum CCL17 levels are also significantly elevated in patients with other Th2-dominated diseases including atopic dermatitis (AD) and asthma and can also serve as an index of disease activity for these two conditions [6, 7]. In spite of these relationships, use of anti-BP180 and CCL17 as serum biomarkers for BP is relatively restricted due to their limited sensitivity and specificity. Accordingly, new types of biomarkers are needed for use in the diagnosis of BP.

MicroRNAs (miRNAs) are endogenous noncoding RNAs that play important roles in gene transcription and expression. miRNAs regulate the expression of at least half of the human transcriptome by either repressing the translation of or degrading multiple-target mRNAs. This silencing mechanism can be determined by assessing the extent of base pairings between the miRNA and the target mRNA, when the mRNA binds to the complimentary target sites located in the 3 untranslated region of the target mRNA [8]. miRNAs are secreted by cells through exosomes and extracellular vesicles, and these secreted miRNAs can remain stable in body fluids with a half-life of ≥5 days [9]. Recently, miRNAs have been shown to be key factors in autoimmunity [10, 11] and, in this way, can be used as diagnostic biomarkers for some autoimmune diseases including pemphigus, systemic lupus erythematosus, and systemic sclerosis [1214]. However, the expression profiles of miRNAs, which could serve as biomarkers in sera of BP patients, have not been reported in the literature. The aim of this study was to investigate differentially expressed miRNAs in sera of BP patients and assess their roles as serum biomarkers for BP by comparing levels of these differentially expressed miRNAs with those of serum CCL17 and anti-BP180 levels.

2. Materials and Methods

2.1. Participants

A total of 65 ethnic Han Chinese active BP inpatients were assessed over the period from June 2016 to June 2018. The diagnostic criteria for BP were in accordance with the established European guidelines [15]. BP patients treated with systemic corticosteroids or immunosuppressants, not in the active stage of the disease, showing malignant tumors and/or other autoimmune diseases were excluded from the study. Eighty healthy ethnic Han Chinese served as controls. The study was approved by the ethics committee of the First Hospital of China Medical University (2018-307-3).

2.2. Clinical Data

Patient-level data were obtained from the electronic medical records. The records included age, sex, ethnicity, disease duration, previous history, IIF titer, and treatment regimens of BP.

2.3. Assessment of BP Disease Activity

Baseline was defined as the time at which a physician started treatment for BP. Disease control was defined as the time at which formation of new lesions or pruritic symptoms ceased and established lesions began to heal [15].

2.4. Quantitative Polymerase Chain Reaction- (qPCR-) Based Array

For screening, a qPCR-based array was performed as previously described [16]. Briefly, total RNA was isolated from serum samples of the 20 active BP patients and 20 healthy controls with a miRNeasy Serum/Plasma Kit (Qiagen, Hilden, Germany). The yield of total RNA was 7-15 ng/μL. The cDNA was obtained from 3 μL of RNA with use of the TaqMan® MicroRNA Reverse Transcription Kit (Applied Biosystems, Beverly, USA). The TaqMan® PreAmp Kit (Applied Biosystems) was used for preamplifying with 2.5 μL of cDNA product per specimen. The qPCR-based array analysis was performed with use of the TaqMan Low Density Array Human MicroRNA Panel (CapitalBio, Beijing, China). A total of 768 known miRNAs were quantified with use of the TaqMan® Human MicroRNA Array A and B (Applied Biosystems). Raw cycle threshold (Ct) values were calculated by SDS 2.4 and RQ manager 1.2 software (Applied Biosystems). A Ct value of ≥35.0 was used as a cut-off.

2.5. qPCR Validation Study

C. elegans miR-39 at (Qiagen) was added to the sera prior to RNA isolation. Total RNA was isolated from serum samples of the 45 active BP patients and 60 healthy controls. Extracted RNA was diluted to a uniform concentration of 2 ng/μL, and the reverse transcription reaction was performed with use of a Reverse Transcription Kit (Applied Biosystems). miRNAs were quantified with use of TaqMan™ MicroRNA assays (Applied Biosystems). The qPCR was performed with a 7900HT Fast Real-Time PCR System (Applied Biosystems). Each qPCR test was performed in duplicate. Ct values were obtained with use of SDS 2.4 and RQ manager 1.2 software (Applied Biosystems). The Ct values () of C. elegans miR-39 in the BP patients and the healthy controls in this study were and , respectively. There were no statistically significant difference () and only a very small variation between the two groups. We chose C. elegans miR-39 to normalize the relative expression levels of miRNAs [16, 17].

2.6. Enzyme-Linked Immunosorbent Assay (ELISA)

Serum levels of CCL17 were determined with use of ELISA kits (R&D Systems, USA). The minimal detectable concentration of CCL17 was 7 pg/mL. Serum levels of anti-BP180 and anti-BP230 were also determined using ELISA kits (MBL Systems, Nagoya, Japan). The normal ranges of anti-BP180 and anti-BP230 were both 0-9 U/mL.

2.7. Statistical Methods

For analysis of qPCR-based arrays, data were standardized as based on the global normalization method. The ΔCt value for each individual miRNA was equal to the Ct value, as obtained by subtracting the mean Ct value from all detected miRNAs in each specimen [16]. qPCR-based array data were analyzed using software R (version 3.2.3). Differentially expressed miRNAs were defined as satisfying a criterion of both a minimum of a twofold change in 2∣-ΔΔCt∣ and a value < 0.05 between BP patients and healthy controls. Relative expression levels of serum miRNAs were normalized to cel-miR-39 (ΔCt) and presented as 2-ΔCt in the qPCR validation study.

Results for nonnormally distributed continuous variables were presented as medians (and interquartile ranges (IQR)), and Mann-Whitney tests were used to compare differences between BP patients and healthy controls. Data from qualitative variables were analyzed with use of Fisher’s exact test. Wilcoxon signed-rank tests were performed for paired data. Receiver operating characteristic (ROC) curves and areas under the curve (AUC) were performed to evaluate the diagnostic performance of miRNAs and CCL17. The maximal Youden’s index was calculated with the ROC curve, with . The optimal cut-off value was set as the threshold with the maximal Youden’s index [16]. Correlations between significantly dysregulated miRNAs and significantly elevated levels of CCL17 and anti-BP180 were analyzed with use of Spearman’s rank correlation coefficient. values < 0.05 were required for results to be considered statistically significant. Statistical analysis was performed with use of the GraphPad Prism 7.0 program.

2.8. Bioinformatics Methods

Targeted genes of significantly dysregulated miRNAs as predicted by computer-aided algorithms were obtained from the TargetScan miRNA database, miRDB, and TarBase. As miRNA databases predict targets based on sequences alone, there is some element of chance involved. To avoid this possibility, we considered only those potential targets reported in the literature as being related to BP. The KEGG pathway database [18] and PathCards database [19] were used to assess pathways of the targeted genes.

3. Results

3.1. Characteristics of BP Patients and Healthy Controls

The demographic, clinical, and laboratory data of the participants in the qPCR-based array analysis and the qPCR validation study are presented in Table 1. The initial therapeutic regimens of the 45 BP patients in the qPCR validation study are listed in Table 2.

3.2. Significantly Dysregulated miRNAs in the qPCR-Based Array Analysis

Totally, 327 miRNAs were detected, of which 15 miRNAs were significantly upregulated and 13 downregulated in the sera of the 20 active BP patients as compared with those obtained in the 20 healthy controls (Figure 1).

3.3. Significantly Dysregulated miRNAs in the qPCR Validation Study

Among the 28 significantly dysregulated miRNAs in the qPCR-based array study, relative expression levels (median and IQR) were all significantly upregulated in baseline sera of the 45 active BP patients in comparison with the 60 healthy controls for serum miR-1291 (0.31 (0.13-0.89) vs. 0.06 (0.05-0.12), , Figure 2(a)), miR-27a-5p (0.19 (0.09-0.33) vs. 0.08 (0.04-0.14), , Figure 2(b)), and miR-423-5p (3.09 (1.74-7.74) vs. 1.73 (0.79-3.08), , Figure 2(c)).

3.4. Baseline Serum Levels of CCL17, Anti-BP180, and Anti-BP230 in the 45 Active BP Patients

Baseline levels (median and IQR) of serum CCL17 in the 45 active BP patients were significantly increased as compared with those of the 60 healthy controls (1368 (410-2088) pg/mL vs. 248.5 (143.4-376.7) pg/mL, , Figure 2(d)). Of the 45 active BP patients, 38 (84.44%) had serum levels of (71.88 (41.6-123.6) U/mL) and 25 (55.56%) had serum levels of (82.48 (50.51-106.6) U/mL).

3.5. Relative Expression Levels of Serum miR-1291, CCL17, Anti-BP180, and Anti-BP230 Are Significantly Decreased in the Disease Control Stage

An additional sample of serum was obtained from the 13 BP patients in the disease control stage. Characteristics of these 13 patients are contained in Table 3. Elevated baseline serum levels were obtained in 10/13 BP patients for anti-BP180 and 9/13 for anti-BP230. In comparison with their baseline levels, relative levels (median and IQR) of serum miR-1291 expression in the 13 BP patients were significantly decreased in the disease control stage (0.74 (0.28-2.66) vs. 0.19 (0.08-0.50), , Figure 3(a)). In contrast, no significant reductions were observed in the disease control stage for relative expression levels of miR-27a-5p (0.34 (0.20-0.50) vs. 0.26 (0.11-0.54), ) or miR-423-5p (4.63 (3.11-6.93) vs. 3.123 (1.55-4.3), ). Accordingly, only miR-1291 appears to provide a reflection of BP disease activity. In addition, as compared with the baseline values, levels of serum CCL17 (1557 (387-5906) pg/mL vs. 246.3 (86.4-514.2) pg/mL, , , Figure 3(b)), anti-BP180 (103.5 (66.08-168.2) U/mL vs. 18.09 (4.80-48.58) U/mL, , , Figure 3(c)), and anti-BP230 (80.87 (27.32-102.9) U/mL vs. 21.61 (2.65-59.63) U/mL, , ) were also all significantly decreased in the disease control stage.

3.6. Areas under the Curve (AUC), Sensitivity, and Specificity of miR-1291 and CCL17

The AUC value for miR-1291 was 0.83 (95% CI, 0.75-0.92, , Figure 4(a)). The sensitivity and specificity of miR-1291 for BP (2-ΔCt at cut-off level of >0.13) were 75.56% and 81.03%, respectively. The AUC value for CCL17 was 0.86 (95% CI, 0.77-0.94, , Figure 4(b)). The sensitivity and specificity of CCL17 for BP (at cut-off level of >549.0 pg/mL) were 73.33% and 98.33%, respectively.

3.7. Correlations between miR-1291 and CCL17 and Anti-BP180

Relative levels of miR-1291 expression showed a statistically significant positive correlation with serum CCL17 levels (, , Figure 5(a)) and anti-BP180 (, , Figure 5(b)). No statistically significant correlations were obtained between miR-1291 and anti-BP230.

3.8. Relative Expression of miR-1291 Levels in Different Subgroups of BP Patients

The 45 active BP patients were divided into different subgroups by sex, with/without hypertension, with/without cerebrovascular disease, with/without mucosal involvement, serum CCL17 levels > or ≤549.0 pg/mL, positive/negative anti-BP180, and positive/negative anti-BP230. Relative expression of serum miR-1291 levels (median and IQR) was significantly upregulated in the 33 patients with as compared with the 12 patients showing CCL17 levels of ≤549.0 pg/mL (0.52 (0.16-1.03) vs. 0.13 (0.08-0.24), ). No other statistically significant differences in relative expression of serum miR-1291 levels were obtained among the other subgroups (Table 4). Of the 34 BP patients with serum levels of (cut-off value), 28 (82.35%) had a serum , while in the 11 BP patients with serum levels of , 5 (45.45%) had serum (). The 45 BP patients were divided into two groups according to the median of miR-1291 (0.31), and then the number of BP patients whose serum or was compared. Compared with the patient group with serum , the percent values of patients with or were both significantly higher in the patient group with serum (Table 5).

3.9. BP-Associated Targeted Genes of miR-1291

We found 1575, 469, and 55 predicted targeted genes of miR-1291 in the TargetScan miRNA database, miRDB, and TarBase, respectively. Based on the criteria described in the methods, signal transducer and activator of transcription 6 (STAT6) and interleukin-13 (IL-13) were selected as potential BP-associated targeted genes of miR-1291 in the TargetScan miRNA database, while IL-13 was selected as a potential BP-associated targeted gene of miR-1291 with miRDB [2, 20, 21]. No BP-associated targeted genes were found in the TarBase. Therefore, we chose STAT6 and IL-13 as possible BP-associated target genes of miR-1291. According to the KEGG and the PathCards database, the signal pathways associated with STAT6 and IL-13 are listed in Table 6. We found that the Janus kinase and signal transducer and activator of transcription (JAK-STAT) signaling pathway and Th2 cell differentiation pathway were shared signal pathways in the KEGG and the PathCards database.

4. Discussion

To the best of our knowledge, this study represents the first report on serum expression profiles of miRNAs in BP patients. Our results demonstrate that the relative expression of miR-1291 was significantly upregulated in the sera of active BP patients and significantly decreased after effective treatment. These miR-1291 responses were associated with a 75.56% sensitivity and 81.03% specificity. Such results suggest that miR-1291 can reflect disease activity in BP patients as well as their responses to treatment.

Significantly elevated serum levels of the Th2 chemokine, CCL17, in active BP patients have been reported previously in the literature [35]. In this study, we show that the relative expression levels of serum miR-1291 were significantly upregulated in the 33 active BP patients with as compared with the 12 active BP patients with . We also found that relative baseline expression levels of serum miR-1291 showed a significant positive correlation with baseline levels of serum CCL17 and anti-BP180 in active BP patients. Therefore, like CCL17 and anti-BP180, miR-1291 may also prove to be a good, if not better, marker for assessing disease activity of BP. The percent values of patients with were significantly higher in the patients with serum than in the patients with serum . Compared with the patients with serum , the percent values of patients with or were both significantly higher in the patients with serum . Therefore, we think the high variability of serum miR-1291 levels in the BP patients reflects the various disease activity of BP.

In this study, STAT6 and IL-13 were chosen as potential BP-associated target genes of miR-1291. The expression of STAT6 has been reported to be significantly elevated in BP skin lesions, but not in BP perilesional skin and normal skin [20]. Dysfunction in regulatory T cells (Tregs) triggers STAT6-mediated activation of autoreactive CD4+ Th2 cells and follicular helper T cells, resulting in the production of autoantibodies to BP antigens [22]. With regard to IL-13, it has been reported that eosinophils are able to secrete IL-13 which promotes Th2 polarization [23], and these IL-13-producing cells are often observed in BP skin lesions. IL-13 contributes to the activation of type 2 macrophages, leading to a positive feedback for the induction of Th2 cells, and maintains the Th2-polarized immunological microenvironment in BP [21]. Given this background information, it was considered important to conduct further experiments to verify the regulatory relationships of these two target genes. We found that STAT6 and IL-13 are involved in the JAK-STAT signaling pathway and Th2 cell differentiation pathway. The JAK-STAT pathway is essential for normal functioning of the immune system. Interestingly, JAK-STAT signaling is dysregulated in BP, which may be related to activation of distinct cytokines and neutrophilic and/or eosinophilic infiltration [20]. The signature cytokines of Th2 cells are IL-4 and IL-13. IL-4 can promote the differentiation of Th2 cells. These Th2 cells can then express the chemokine receptors, CCR4 and STAT6, which serve as the major signaling pathway of IL-4-mediated Th2 differentiation. For IL-13, it has been reported that STAT6 induces GATA-3 expression and GATA-3 binds to the promoters of IL-13 and the enhancer of IL-4, thereby promoting their transcription [24].

There were some limitations in our study. Bullous Pemphigoid Disease Area Index (BPDAI) scoring was not performed in all patients. As many patients were from diverse locations, sera collected in the disease control stage were not obtained in all patients but only in some patients who were followed up regularly. As BP occurs mainly in the elderly, we were unable to recruit enough age-matched healthy controls for this study. In addition, no investigation was performed regarding in vivo functions of miR-1291 and the genes targeted by miR-1291. It is not clear what kind of cells or tissues release miR-1291 into serum, and there is still a long way to go between predicting the target genes and understanding the pathogenesis. This is an exploratory study to find the dysregulated microRNAs in sera of BP, making a preliminary discussion of the putative target genes, and give some hints to other researchers. Future studies of target genes and related pathways will be carried out in immune-related cells and lesional tissues of BP patients. While remaining aware of these limitations, we feel that the data presented provide important new information on initial investigations into the expression of miRNAs within the serum of BP patients.

5. Conclusions

Relative expression levels of serum miR-1291 reflect disease activity and treatment responses within BP patients. miR-1291 may function as a serum biomarker for BP.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was funded by the National Natural Science Fund of China (81470142) and Liaoning Provincial Natural Science Foundation of China (2015020509). This manuscript was revised for submission in English by the ED-IT Editorial Service.