Plasma-Derived Exosomal hsa-miR-4488 and hsa-miR-1228-5p: Novel Biomarkers for Dermatomyositis-Associated Interstitial Lung Disease with Anti-Melanoma Differentiation-Associated Protein 5 Antibody-Positive Subset

Zhong, Danli; Wu, Chanyuan; Xu, Dong; Bai, Jingjing; Wang, Qian; Zeng, Xiaofeng

doi:https://doi.org/10.1155/2021/6676107

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Data Availability Ethical Approval Consent Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 6676107 | https://doi.org/10.1155/2021/6676107

Plasma-Derived Exosomal hsa-miR-4488 and hsa-miR-1228-5p: Novel Biomarkers for Dermatomyositis-Associated Interstitial Lung Disease with Anti-Melanoma Differentiation-Associated Protein 5 Antibody-Positive Subset

Danli Zhong,¹Chanyuan Wu,¹Dong Xu,¹Jingjing Bai,¹Qian Wang,¹and Xiaofeng Zeng¹

Academic Editor: Emiliano Antiga

Received03 Oct 2020

Accepted27 Jun 2021

Published29 Jul 2021

Abstract

The present study is aimed at profiling circulating exosome-derived microRNAs (miRNAs/miRs) from patients with dermatomyositis (DM), in particular those complicated with interstitial lung disease (ILD) with anti-melanoma differentiation-associated protein 5 (MDA5) antibody-positive. Fifteen participants were enrolled, including five patients with DM complicated with ILDs prior to treatment with circulating anti-MDA5 antibody-positive status [DM-ILD-MDA5 Ab(+)], five DM patients without ILDs who were negative for 16 detectable myositis-specific antibodies [DM-nonILD-MSA16(-)], and five age- and gender-matched healthy donor controls (HCs). The characteristics of the exosomes extracted by Ribo™ Exosome Isolation Reagent were identified using transmission electron microscopy (TEM), nanoparticle tracking analysis (NTA), and flow cytometry. Differentially expressed miRNAs, determined by next-generation deep sequencing, were identified through the criteria of and . A total of 38 miRNAs were significantly upregulated in exosomes from patients with DM-ILD-MDA5 Ab(+) compared to those from HC, while 21 miRNAs were significantly downregulated. Compared to exosomes derived from patients with DM-nonILD-MSA16(-), 51 miRNAs were significantly upregulated and 33 miRNAs were significantly downregulated from patients with DM-ILD-MDA5 Ab(+). A total of 73 exosomal miRNAs were significantly differentially expressed between DM-nonILD-MSA16(-) and HC. In particular, two miRNAs, Homo sapiens- (hsa-) miR-4488 and hsa-miR-1228-5p, were common differentially expressed miRNAs among three comparisons. GO and KEGG analyses suggested that several pathways may contribute the pathogenesis of DM-ILD-MDA5 Ab(+) and DM-nonILD-MSA16(-), while PPI network analysis of hsa-miR-4488 and hsa-miR-1228-5p indicated that their predicted target genes, DExD-box helicase 39B and MDM2, may be involved in the mechanisms of DM-ILD-MDA5 Ab(+).

1. Introduction

Exosomes, small vesicles (30-150 nm in diameter) of endocytic origin, are formed in endosomal multivesicular compartments and are secreted into the extracellular space when these compartments fuse with the plasma membrane. Exosomes may be purified from cell cultures and body fluids, including plasma/serum, saliva, urine, and cerebrospinal fluid. Exosomes contain various functional small RNAs, mRNAs, and functional proteins, thereby mediating long distant cell-cell communication. MicroRNAs (miRNAs/miRs) are small noncoding RNAs that posttranscriptionally regulate target genes by base-pairing with 3-untranslated regions of target mRNA [1]. It has been reported that exosomal miRNAs are extra stable in various storage conditions [2] and are protected from RNase treatment, thereby providing a protected and enriched source of miRNA [3]. These results suggest that exosome-derived miRNA (exRNA) has the potential for diagnosis, treatment, and prognosis prediction as a stable biomarker, and it overlaps with the even newer field of exRNA-mediated communication.

Dermatomyositis (DM) encompasses a spectrum of rare, heterogeneous autoimmune diseases characterized by proximal symmetric muscle weakness, skin rash, multiple organ involvement, and paraneoplastic syndromes [4]. The etiology of DM is unknown; however, cross-talk between genetic and environmental factors contributes to the pathology [5]. Over the past decades, new subsets of DM have become increasingly recognized. Patients with DM may be subclassified on the basis of clinical and laboratory features, among which myositis-specific autoantibodies (MSAs) have a crucial role in serologically defined subsets of DM [6]. When autoantibodies to a 140 kDa protein were first discovered in 2005 in patients with clinically amyopathic DM [7] who predominately displayed cutaneous manifestations for abnormally long periods of time (≥6 months) without accompanying muscle weakness [8], subsequent evidence demonstrated that melanoma differentiation-associated gene-5 (MDA5) was the target autoantigen [9]. MDA5, a crucial member of the retinoic acid-inducible gene-like receptor family, recognizes viral double-strand RNA with the production of type I interferon, thereby protecting the host against invasion of virus [10]. Patients with anti-MDA5 autoantibodies (anti-MDA5 Ab) are at risk of an acute, rapidly progressive, and potentially fatal form of interstitial lung disease (ILD) that requires aggressive systemic immunosuppressive therapy and is associated with high morbidity and mortality rates [11, 12]. It has been reported that the prevalence of ILD in DM with anti-MDA5 Ab-positive status was higher in Asian populations than that in Europeans [13–15]. The absence of anti-MDA5 Ab in other connective tissue disorders and the presence of anti-MDA5 Ab being mutually exclusive of other MSAs in patients with DM suggested that a specific immune response may govern this subtype of DM. However, the underlying mechanisms of DM-associated ILD with anti-MDA5 Ab remain elusive. As the lack of homogeneous grouping based on the antibody profiles perturbs the diagnostic, therapeutic, and prognostic stratification of these patients, exosome-derived miRNAs of circulating plasma were profiled in the present study using next-generation sequencing (NGS) from two subsets of DM patients: patients with DM complicated with ILDs prior to treatment with circulating anti-MDA5 Ab-positive status [DM-ILD-MDA5 Ab(+)] and patients without ILDs who were negative for 16 detectable MSA [DM-nonILD-MSA16(-)]. Patients were considered to be DM-nonILD-MSA16(-) when 16 MSAs were undetected in serum, including anti-aminoacyl-transfer (t)RNA synthetase autoantibodies (anti-Jo-1, anti-PL-7, anti-PL-12, anti-EJ, and anti-OJ), anti-Ro-52, anti-signal recognition particle antibodies (anti-SRP), anti-PM-SCL75, anti-PM-SCL100, anti-KU, anti-small ubiquitin-like modifier activating enzyme antibodies (anti-SAE), anti-nuclear matrix protein 2 antibodies (anti-NXP2), anti-MDA5, transcriptional intermediary factor-1γ autoantibodies (anti-TIF1γ), anti-MI-2α, and anti-MI-2β antibodies.

2. Materials and Methods

2.1. Characteristics of Eligible Patients

A total of 10 patients diagnosed with DM and five healthy volunteer controls (HCs) were included in the present study from Peking Union Medical College Hospital (PUMCH). The individuals in control group (age- and sex-matched) did not suffer autoimmune or inflammatory diseases. The consecutive patients were enrolled from the Department of Rheumatology, Peking Union Medical College Hospital affiliated to Peking Union Medical College & Chinese Academy of Medical Sciences (Beijing, China) between November 1^st, 2017 and January 1^st, 2018. All patients and controls provided written informed consent according to the Declaration of Helsinki. This study was approved by the Ethics Committee of Peking Union Medical College Hospital (Beijing, China; no. S-478).

The diagnosis of DM was determined according to the Bohan and Peter criteria [16]. ILD was considered if patients fulfilled the following criteria: (i) radiological signs on high-resolution computerized tomography and (ii) impaired pulmonary function (total lung capacity, % of predicted and % of predicted).

Serum MSA profiles including anti-aminoacyl-tRNA synthetase autoantibodies (anti-Jo-1, anti-PL-7, anti-PL-12, anti-EJ, and anti-OJ), anti-RO-52, anti-SRP, anti-PM-SCL75, anti-PM-SCL100, anti-KU, anti-SAE, anti-NXP2, anti-MDA5, anti-TIF1γ, anti-MI-2β, and anti-MI-2α were measured using commercial kits by line blot assays (Myositis Profile Euroline kit; Euroimmun AG) [17]. Patients with DM were divided into two groups according to lung involvement and the presence of MDA5 autoantibodies. A total of 5 patients diagnosed as DM-ILD-MDA5 Ab(+) were eligible and naive to treatment with prednisone or immunosuppressant drugs prior to sample collection. None of the anti-MDA5 Ab(+) DM patients had MSAs other than anti-MDA5 antibody. The other five patients diagnosed with DM without ILD were negative for serum MSAs. The characteristics of the patients included are summarized in Table 1. Patients with malignancy, pulmonary infection, or overlapping connective tissue diseases were excluded.

2.2. Blood Collection and Isolation of Plasma Exosomes

A total of 20 ml of whole blood were collected from the cubital vein of patients using blood collection tubes containing K₂EDTA. Blood samples were centrifuged at 4,000 × g at 4°C for 5 min using a swinging bucket rotor within 4 h of the blood draw. The upper supernatant (plasma) was carefully transferred to a new tube with a conical bottom without disturbing the intermediate buffy coat layer and was stored at -80°C prior to the isolation of plasma exosomes.

Plasma exosomes were extracted using a commercial kit (Ribo™ Exosome Isolation Reagent; RiboBio) following the manufacturer’s protocol. In brief, plasma samples were thawed from -80°C and were centrifuged at 2,000 × g for 20 min at 4°C to remove cells and cell debris. Subsequently, the supernatants were transferred to sterile polyethylene tubes with addition of 1/3 volume of Ribo™ Exosome Isolation Reagent (RiboBio). The samples were incubated at 4°C for 30 min. Subsequently, the samples were centrifuged at 15,000 × g at 4°C for 2 min to obtain white pellets (exosomes) following removal of the upper supernatant. The isolated exosome samples were stored at -80°C for the subsequent assays. The study was designed to be consistent with the minimal experimental requirement for extracellular vesicles from the statement from the International Society for Extracellular Vesicles [18].

2.3. Identifications of Exosomes Using TEM

Biological transmission electron microscopy (TEM) was performed to determine the characteristics of exosomes. Exosomes were deposited onto a copper grid for 5 min at room temperature. Subsequently, the grid was negatively stained with 2% uranyl acetate for 5 min. The grid was then imaged under a biological TEM (Tecnai™ G² spirit BioTWIN; FEI Co.) with an accelerating voltage set to 120 kV.

2.4. Nanoparticle Tracking Analysis (NTA)

Quantification and size determination of plasma exosomes were performed using a NanoSight NS300 instrument (Malvern Instruments, Inc., UK) equipped with the NTA 3.0 software (Malvern Instruments, Inc.). In brief, exosomes were diluted with deionized water by 1 : 10,000, which provided particle concentrations of /ml prior to examination using the NTA system. Exosome samples were then injected into a NanoSight sample cubicle (Malvern Instruments, Inc., UK) by a syringe. A total of 5 recordings of 30 sec each were performed, and data from at least 5,000 individual particle tracks per sample were analyzed. The data of exosomal size distribution and mean number of particles per milliliter were expressed as the .

2.5. Flow Cytometry

Plasma exosomes isolated with the commercial kits were analyzed using a BD Accuri C6 flow cytometer (BD Biosciences) with the FlowJo software 10.0.7 (FlowJo LLC). Exosomes were stained by CD63-FITC (BD Biosciences; cat. no. 557288) and CD81-FITC (BD Biosciences; cat. no. 551108) for 20 min at 4°C, both at a dilution of 1 : 500. In brief, exosomes were resuspended in 100 μl PBS that had been filtered using 0.22 μm filters. The samples were incubated with FITC-conjugated CD63 or CD81 antibody in PBS for 15-20 min at 4°C with gentle agitation. After washing in PBS, the samples were resuspended in 100 μl PBS and analyzed using flow cytometry.

2.6. Exosomal RNA Extraction, Library Preparation, and NGS

Total RNA was extracted from isolated exosomes using TRIzol reagent (Invitrogen; Thermo Fisher Scientific, Inc.) according to the manufacturer’s protocol. The purity and concentration of RNA was assessed using the Nanodrop-1000 (Thermo Fisher Scientific, Inc.). Each RNA sample had an absorbance at and an A260/A230 ratio of >2.0. RNA integrity was evaluated using the Agilent 2200 TapeStation (Agilent Technologies, Inc.), and each sample had an . In brief, RNAs from exosomes were ligated with 3RNA adapter and followed by 5adapter ligation. Subsequently, the adapter-ligated RNAs were subjected to reverse transcription PCR and amplified with a low-cycle as previously described [19]. Then, PCR products were size-selected using PAGE according to the instructions of the NEBNext® Multiplex Small RNA Library Prep Set for Illumina® (New England Biolabs, Inc.). The purified library products were evaluated using the Agilent 2200 TapeStation and were diluted to 10 pM for cluster generation in situ on the HiSeq2500 single-end flow cell followed by sequencing ( bp) on the HiSeq 2500 system.

2.7. Quality Control and Mapping to the Reference Genome

The raw reads were preprocessed to trim 3 and 5 adapters, and the low-quality reads were removed to obtain clean reads. The GC content, Q20, and Q30 of the clean data were also calculated. miRBase (version 21.0; http://www.mirbase.org) was used to build the index of the reference genome for the components and expression profiles of all types of small RNAs [20]. The miRNA expression levels were estimated as the number of reads per million (RPM) using the following formula:

2.8. Differentially Expressed (DE) miRNA Analysis

The edgeR package in Bioconductor was used to analyze DE miRNAs between two groups [21]. DE miRNAs were defined as those fulfilling the criteria of and . Due to the low abundance of small noncoding RNAs in exosomes, the number of DE miRNAs based on adjusted values was too low to perform any further analysis. In fact, thresholds in high-throughput sequencing may be set at unadjusted values to determine the DE miRNAs, which has been reported previously [22]. If the expression level was too low to detect (displayed as zero), the expression level was set as 0.001 in order to calculate the fold change.

The R package multiMiR (http://multimir.ucdenver.edu/) was used to predict candidate target genes (CTGs) of miRNA. This database compiles 14 existing miRNA-target prediction programs from experimental, computational, drug, and disease resources, including DIANA-microT, EIMMo, MicroCosm, miRanda, miRDB, PicTar, PITA, and Targetscan [23].

2.9. Functional Enrichment Analysis

To explore the potential function of miRNAs, Gene Ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG; http://www.genome.jp/) enrichment analyses were performed on predicted target genes of differentially expressed miRNA using the clusterProfiler package [24]. GO functional analysis consists of the category biological process (BP), cellular component (CC), and molecular function (MF). The results of the analysis were extracted under the threshold of an adjusted value <0.05.

2.10. Construction of miRNA-mRNA Regulatory Network and Protein-Protein Interaction (PPI) Network Analysis

To establish the miRNA-mRNA network, the analysis was focused on the DE miRNAs among three comparisons. Cytoscape was used to draw the network for visualization. Subsequently, the DE miRNAs were further analyzed using the Search Tool for the Retrieval of Interacting Genes and Proteins (STRING) database (http://string-db.org) to determine the interactive relationships of CTGs. Based on the results obtained with STRING, the PPI network was visualized utilizing the Cytoscape software (version 3.7.1) [25].

2.11. Statistical Analysis

Demographic data were expressed as the , while differences between two groups were assessed by the unpaired -test as appropriate. values for differential analysis were two-sided, and was considered to indicate statistical significance. All statistical analyses were performed using R packages for R 3.5.1 or GraphPad Prism 7.0 (GraphPad Software, Inc.).

3. Results

3.1. Characterization of Plasma Exosomes

A total of ten patients and five healthy controls were recruited in the present study. The basic demographic characteristics of enrolled individuals are presented in Table 1. There were no significant differences in the sex and age distributions among these three groups. Plasma-derived exosomes were characterized by TEM, NTA, and flow cytometry. The morphology and ultrastructure of the exosomes was assessed using TEM technology, and the characteristic round- or oval-shaped bodies were visualized. The NTA results indicated that exosomes obtained with the commercial kit were nm in size with a concentration of particles/ml, which was expected according to the kit manual. Surface tetraspanins on exosomes, CD63 and CD81, were confirmed using flow cytometry (Figure 1). These results suggested that the commercial kits for exosome separation are reliable and generated acceptable yields as expected.

(a)

(b)

(c)

3.2. Classification of Circulating Exosomal Small RNAs

Exosomes have been demonstrated to encompass various types of small RNAs, including miRNA, tRNA, ribosomal (r)RNA, small nuclear (sn)RNA, small nucleolar (sno)RNA, Piwi-interacting (pi)RNA, Y_RNA, and unannotated RNAs [26, 27]. Therefore, the exosomal small RNAs isolated from patients of the three groups were classified by mapping the clean reads to a noncoding RNA database, including miRBase version 21.0 (miRNA database; http://www.mirbase.org), Rfam 12.1 (tRNA, rRNA, snRNA, and snoRNA databases; http://rfam.xfm.org), and piRNAbank (piRNA database; http://pirnabank.ibab.ac.in). The proportion of miRNA in the total small RNA isolated from exosomes in DM-ILD-MDA5 Ab(+), DM-nonILD-MSA16(-), and HC was 75.12, 72.10, and 53.85%, respectively (Figure 2). Using miRBase version 21.0 as a reference, 912, 895, and 876 known miRNAs were identified in DM-ILD-MDA5 Ab(+), DM-nonILD-MSA16(-), and HC plasma exosomes, respectively.

(a)

(b)

(c)

Figure 2

Pie charts displaying the distribution of different classes of small RNAs derived from exosomes in (a) HC, (b) DM-nonILD-MSA16(-), and (c) DM-ILD-MDA5 Ab(+). HCs: healthy controls; DM-ILD-MDA5 Ab(+): DM patients with ILDs and anti-melanoma differentiation-associated protein 5 antibody-positive status; DM-nonILD-MSA16(-): DM patients without ILDs who were negative for 16 detectable myositis-specific antibodies; ILD: interstitial lung disease; DM: dermatomyositis; miRNA: microRNA; rRNA: ribosomal RNA; snRNA: small nuclear RNA; snoRNA: small nucleolar RNA; piRNA: Piwi-interacting RNA; tRNA: transfer RNA.

3.3. DE miRNA

After filtering out miRNAs with low expression ( in circulating exosomes), using and as a cut-off threshold, numerous miRNAs were indicated to be significantly differentially expressed between two groups. Among the DE miRNAs, 38 miRNAs were upregulated and 21 were downregulated in DM-ILD-MDA5 Ab(+) exosomes as compared to HC exosomes. Compared to patients with DM-nonILD-MSA16(-), 51 miRNAs were significantly upregulated and 33 miRNAs were significantly downregulated in patients with DM-ILD-MDA5 Ab(+). A total of 73 miRNAs were significantly differentially expressed between DM-nonILD-MSA16(-) and HC, including 50 upregulated and 23 downregulated miRNAs in patients with DM-nonILD-MSA16(-) vs. HC. A Venn diagram demonstrated that exosomal miRNA profiles from DM-ILD-MDA5 Ab(+), DM-nonILD-MSA16(-), and HC had two common miRNAs [Homo sapiens- (hsa-) miR-4488 and hsa-miR-1228-5p; Figure 3]. In the present study, hsa-miR-1228-5p was upregulated in the DM-ILD-MDA5 Ab(+) group (, ) and downregulated in the DM-nonILD-MSA16(-) group (, ) when compared with the HC group, while hsa-miR-4488 was upregulated in the DM-ILD-MDA5 Ab(+) group as compared with the DM-nonILD-MSA16(-) group (, ) and the HC group (, ). The DE miRNAs are displayed in Tables SI-III. Heatmaps revealed that the top 50 DE miRNAs had high intragroup similarity and low intergroup similarity (Figure 4).

Figure 3

Comparison between identified miRNAs released from exosomes of DM-ILD-MDA5 Ab(+), DM-nonILD-MSA16(-), and HC. Venn diagram demonstrating that the homologous exosomes from three comparisons have two common miRNAs (hsa-miR-1228-5p and hsa-miR-4488). HC: healthy controls; DM-ILD-MDA5 Ab(+): DM patients with ILDs and anti-melanoma differentiation-associated protein 5 antibody-positive status; DM-nonILD-MSA16(-): DM patients without ILDs who were negative for 16 detectable myositis-specific antibodies; ILD: interstitial lung disease; DM: dermatomyositis; miRNA/miR: microRNA; hsa: Homo sapiens.

(a)

(b)

(c)

Figure 4

Heatmaps of the top 50 significantly differentially expressed miRNAs between pairs of groups. (a) HC vs. DM-ILD-MDA5 Ab (+), (b) DM-nonILD-MSA16 (-) vs. DM-ILD-MDA5 Ab(+), and (c) HC vs. DM-nonILD-MSA16(-). Gene expression data were obtained using next-generation sequencing on the Illumina HiSeq 2500 platform. The miRNAs are displayed on the vertical axis and the different numbered subjects on the horizontal axis. HCs: healthy controls; DM-ILD-MDA5 Ab(+): DM patients with ILDs and anti-melanoma differentiation-associated protein 5 antibody-positive status; DM-nonILD-MSA16(-): DM patients without ILDs who were negative for 16 detectable myositis-specific antibodies; ILD: interstitial lung disease; DM: dermatomyositis; miRNA/miR: microRNA; hsa: Homo sapiens.

3.4. Functional Enrichment of Top 10 DE miRNA-Targeted Genes

Figure 5 presents the top 10 GO terms in the categories of BP (Figure 5(a)), CC (Figure 5(b)), and MF (Figure 5(c)) for CTGs of DE miRNAs, indicating certain similarities and differences among the various subsets when considering the different control groups. In terms of MF, CTGs of DE miRNAs in DM-ILD-MDA5 Ab(+) were substantially associated with protein binding, DNA binding, catalytic activity, organic cyclic compound binding, ion binding, cation binding, metal ion binding, nucleic acid binding, and heterocyclic compound binding when compared to HC and DM-nonILD-MSA16(-). In comparison to HC, CTGs of DE miRNAs in DM-nonILD-MSA16(-) were enriched in enzyme binding rather than DNA binding. Regarding the category of CC, the top 10 terms in DM-nonILD-MSA16(-) vs. DM-ILD-MDA5 Ab(+) and HC vs. DM-nonILD-MSA16(-) were the same in the cell part, intracellular part, intracellular organelle, cytoplasm, organelle part, membrane, cytoplasmic part, and nucleus; membrane-bounded organelle and intracellular membrane-bounded organelle correlated with DM-ILD-MDA5 Ab(+) exosomes when referring to CC in comparison to HC plasma exosomes. In the category of BP, cellular process, single-organism cellular process, regulation of cellular process, metabolism process, developmental processes, and biological regulation were associated with DM-ILD-MDA5 Ab(+) when compared to HC, while biosynthetic process was relevant in the other two comparisons. KEGG analysis revealed that metabolic pathway, viral infection, and oncogenesis have crucial roles in the pathogenesis of DM, while Wnt signaling and mTOR signaling were indicated to contribute to the pathogenesis in the ILD subset (Figure 6).

(a)

(b)

(c)

Figure 5

Top 10 Gene Ontology terms in the categories of (a) biological process, (b) cellular component, and (c) molecular function regulated by candidate target genes of significantly differentially expressed microRNAs in plasma exosomes. HC: healthy controls; DM-ILD-MDA5 Ab(+): DM patients with ILDs and anti-melanoma differentiation-associated protein 5 antibody-positive status; DM-nonILD-MSA16(-): DM patients without ILDs who were negative for 16 detectable myositis-specific antibodies; ILD: interstitial lung disease; DM: dermatomyositis.

(a)

(b)

(c)

Figure 6

Top Kyoto Encyclopedia of Genes and Genomes pathways regulated by candidate target genes of significantly differentially expressed microRNAs. (a) HC vs. DM-ILD-MDA5 Ab(+). (b) DM-nonILD-MSA16(-) vs. DM-ILD-MDA5 Ab(+). (c) HC vs. DM-nonILD-MSA16(-). HCs: healthy controls; DM-ILD-MDA5 Ab(+): DM patients with ILDs and anti-melanoma differentiation-associated protein 5 antibody-positive status; DM-nonILD-MSA16(-): DM patients without ILDs who were negative for 16 detectable myositis-specific antibodies; ILD: interstitial lung disease; DM: dermatomyositis; PKG: protein kinase G; Fox: forkhead box; cGMP: cyclic guanosine monophosphate; cAMP: cyclic adenosine monophosphate; HTLV: human T-lymphotropic virus.

3.5. Construction of miRNA-mRNA Regulatory Network and PPI Network Analysis

Based on the theory of miRNA function, regulatory loops with miRNA-mRNA regulation pairs were constructed. Among the three comparisons, hsa-miR-4488 and hsa-miR-1228-5p were indicated to be the common DE miRNAs. Therefore, the intersecting miRNAs were selected to establish miRNA-mRNA regulatory networks (Figure 7) and PPI networks for CTGs (Figure 8). For the miR-4488-mRNA regulation circuit, the interaction between DExD-box helicase 39B (DDX39B) and hsa-miR-4488 indicated the strongest evidence from six experiments, followed by retinoid X receptor β, RNA polymerase II subunit F, menin 1, cyclin-dependent kinase inhibitor 1A, G protein subunit β 2, G protein subunit α i2, and spalt-like transcription factor 4. In addition, the PPI network based on the CTGs of miR-4488 indicated that H2A histone family member X (H2AFX) was a hub gene because it had the most nodes. In the miR-1228-5p-mRNA regulation circuit, zinc finger and BTB domain containing 22 (ZBTB22) had the strongest evidence for the interaction with miR-1228-5p, followed by solute carrier family 9 member A3, occludin, neurofibromin 2, and macrophage migration inhibitory factor, which were all evidenced twice according to multiMiR database. The PPI network for CTGs of miR-1228-5p indicated that MDM2 was a key gene because it had the most nodes.

(a)

(b)

(a)

(b)

4. Discussion

To the best of our knowledge, the present study was the first to report on circulating exosomes from patients with DM using a small RNA sequencing strategy. It was revealed that circulating exosome miRNA profiles had a distinct pattern between serologically defined DM subsets and HCs. Based on the circulating exosome-derived small RNA classification, eight classes of exosomal small RNA were identified: miRNA, rRNA, tRNA, snRNA, snoRNA, piRNA, Y_RNA, and unannotated small RNA. Furthermore, the DE miRNAs of plasma exosomes and their CTGs were identified. After filtering out miRNAs with low expression ( reads per million in circulating exosomes), using and as the cut-off threshold, several miRNAs were identified to be significantly differentially expressed between two groups. The intersecting miRNAs, hsa-miR-4488 and hsa-miR-1228-5p, were used to construct miRNA-mRNA regulatory networks.

DM is a heterogeneous group of autoimmune disorders characterized by muscle weakness, cutaneous manifestations, internal organ involvement, and paraneoplastic syndrome. Based on MSAs, occurring in idiopathic inflammatory myositis exclusively, DM may be subclassified according to their strong association with different clinical manifestations. Patients who are serologically positive for anti-MDA5 antibody usually have the acute, rapidly progressive, and potentially fatal form of ILD, which was reported to be much worse than that of antisynthetase syndromes [28]. Chest biopsies from patients with anti-MDA5 antibody indicated similar pathological findings to those for acute respiratory distress syndrome, including diffuse alveolar damage and prominent hyaline membrane deposition [29]. However, little is known about the pathogenesis of DM-ILD-MDA5 Ab(+).

Exosomes originating from various subcellular compartments are regulators of transfer of biological information locally and remotely. Exosomes are small membrane vesicles (30-150 nm in diameter) secreted by most nucleated cells [4–6]. In the present study, exosomes enriched with the reagent kit from DM patients and HCs were nm in diameter with concentrations of particles/ml. The exosomal markers CD63 and CD81 were present in those particles, which was consistent with the known characteristics of exosomes.

Exosomal RNAs, which are more stable than cellular RNA, have been reported as novel biomarkers or therapeutic targets in various diseases [30, 31]. The present study demonstrated that exosome-derived miRNAs identified through NGS technology may become potential biomarkers. Among the DE miRNAs, compared to HC exosomes and DM-nonILD-MSA16(-) exosomes, 38 and 51 miRNAs were significantly upregulated, while 21 and 33 miRNAs were significantly downregulated in DM-ILD-MDA5 Ab(+) exosomes, respectively. PPI network analysis revealed that hsa-miR-4488 and hsa-miR-1228-5p may have crucial roles in the pathogenesis of DM-ILD-MDA5 Ab(+) through 112 and 60 candidate target genes, respectively. It was identified that DDX39B may possibly be a target gene of hsa-miR-4488 according to the multiMiR database. It has been reported that DDX39B interacts with the pattern recognition receptor pathway to inhibit NF-κB signaling [32]. In the present study, hsa-miR-4488 was significantly upregulated in DM-ILD-MDA5 Ab(+) when compared to DM-nonILD-MSA16(-) and HC, suggesting that miR-4488 may contribute to systemic inflammation in DM-ILD-MDA5 Ab(+) by upregulating NF-κB signaling through repression of DDX39B. For its part, DDX39B promotes global translation and cell proliferation through upregulating preribosomal RNA levels, thereby contributing to oncogenesis [33], similar to H2AFX as a protooncogene [34]. hsa-miR-1228-5p was upregulated in DM-ILD-MDA5 Ab(+) and was downregulated in DM-nonILD-MSA16(-) when compared to HC. Furthermore, ZBTB22 and MDM2 may be important target genes in the regulation of hsa-miR-1228-5p. MDM2, as a well-known oncogene, may also promote type 17 T-helper cell differentiation by suppressing p53-dependent transcription factor STAT5 signaling cascades, leading to activation of rheumatoid arthritis [35, 36]. However, only scarce information is available for ZBTB22. Furthermore, it was indicated that hsa-miR-4508 was a significant highly expressed miRNA in exosomes. It has been reported that circulating hsa-miR-4508 was significantly upregulated in acute myocardial infarction [37] and patients receiving dexmedetomidine [38]. However, its functions remain elusive. It was hypothesized that hsa-miR-4508 may participate in acute stress.

GO and KEGG pathway enrichment was performed on CTGs of DE exosomal miRNAs to explore the role of miRNAs packaged into exosomes from patients with DM, particularly in the ILD subset. GO analysis indicated that metabolism had a crucial role in the categories BP and CC. Several differentially regulated KEGG pathways were identified. The Wnt signaling pathway, which was indicated to be upregulated in patients with DM-associated ILD [39], and the cyclic adenosine monophosphate signaling pathway, which has been known to antagonize pulmonary fibrosis activated by phosphodiesterase 4 [40], were involved in DM-ILD-MDA5 Ab(+) exosomes when compared to both HC and DM-nonILD-MSA16(-) exosomes, suggesting that these signaling pathways may be unique to ILD, particularly in the DM-ILD-MDA5 Ab(+) subset. Transforming growth factor-β signaling, which was reported to be upregulated in muscle samples with sporadic inclusion body myositis [41], and mTOR signaling may be potential pathways of autoimmune diseases [42] and were included in the CTGs of DE miRNAs in the DM-nonILD-MSA16(-) exosomes, suggesting that these pathways may contribute to the pathogenesis of DM.

In conclusion, in the present study, exosomal miRNAs from plasma were profiled in two DM subsets [DM-ILD-MDA5 Ab(+) and DM-nonILD-MSA16(-)]. Selective packaging of miRNAs into exosomes, particularly hsa-miR-1228-5p and hsa-miR-4488, was observed. The results may provide promising exosomal miRNA targets for the diagnosis, prevention, and therapy of DM, particularly for the DM-ILD-MDA5 Ab(+) subset.

5. Limitation

This study has certain limitations. Firstly, no external data including quantitative PCR and expression matrix involved these DE miRNAs, thereby further work is warranted for the validation of hsa-miR-1228-5p and hsa-miR-4488 in DM-ILD-MDA5 Ab(+) subset. Secondly, the enrolled patients mainly came from China with little representation of the epigenetic modification in DM. In future, further work is warranted to explore the role of exosomal miRNAs in the pathogenesis of DM.

Data Availability

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Ethical Approval

All study procedures in this study were approved by the Ethics Committee of Peking Union Medical College Hospital affiliated to Peking Union Medical College & Chinese Academy of Medical Sciences (Beijing, China).

All patients provided written informed consent form prior to participation in the study.

Conflicts of Interest

The authors declare they have no competing interests.

Authors’ Contributions

DZ collected samples, performed experiments, wrote the manuscript, and revised the manuscript. CW and DX participated in the bioinformatics analysis and produced figures. JB collected clinical information. QW enrolled the patients, designed the study, interpreted data, and revised the manuscript. XZ interpreted data and revised the manuscript. All authors read and approved the final manuscript.

Acknowledgments

The study received grants from the Chinese National Key Technology R&D Program of the Ministry of Science and Technology (grant no. 2017YFC0907604), the National Science and Technology Major Project of the Ministry of Science and Technology of China (grant no. 2019ZX09734001-002-004), the Medical and Health Science and Technology Innovation Project of the Chinese Academy of Medical Sciences (grant no. 2019-I2M-2-008), and the National Natural Science Foundation of China (grant nos. 81601430 and 81471615).

Supplementary Materials

Supplementary Materials Table S1. All differentially expressed miRNAs in HC versus DM-ILD-MDA5 Ab (+). Table S2. All differentially expressed miRNAs in DM-nonILD-MSA16 (-) versus DM-ILD-MDA5 Ab (+). Table S3. All differentially expressed miRNAs in HC versus DM-nonILD-MSA16 (-). (Supplementary Materials)

References

D. P. Bartel, “MicroRNAs: genomics, biogenesis, mechanism, and function,” Cell, vol. 116, no. 2, pp. 281–297, 2004.
View at: Publisher Site | Google Scholar
Q. Ge, Y. Zhou, J. Lu, Y. Bai, X. Xie, and Z. Lu, “miRNA in plasma exosome is stable under different storage conditions,” Molecules, vol. 19, no. 2, pp. 1568–1575, 2014.
View at: Publisher Site | Google Scholar
L. Cheng, R. A. Sharples, B. J. Scicluna, and A. F. Hill, “Exosomes provide a protective and enriched source of miRNA for biomarker profiling compared to intracellular and cell-free blood,” J Extracell Vesicles, vol. 3, no. 1, 2014.
View at: Publisher Site | Google Scholar
I. E. Lundberg, M. de Visser, and V. P. Werth, “Classification of myositis,” Nature Reviews Rheumatology, vol. 14, no. 5, pp. 269–278, 2018.
View at: Publisher Site | Google Scholar
F. W. Miller, J. A. Lamb, J. Schmidt, and K. Nagaraju, “Risk factors and disease mechanisms in myositis,” Nature Reviews Rheumatology, vol. 14, no. 5, pp. 255–268, 2018.
View at: Publisher Site | Google Scholar
N. J. McHugh and S. L. Tansley, “Autoantibodies in myositis,” Nature Reviews Rheumatology, vol. 14, no. 5, pp. 290–302, 2018.
View at: Publisher Site | Google Scholar
S. Sato, M. Hirakata, M. Kuwana et al., “Autoantibodies to a 140-kd polypeptide, CADM-140, in Japanese patients with clinically amyopathic dermatomyositis,” Arthritis and Rheumatism, vol. 52, no. 5, pp. 1571–1576, 2005.
View at: Publisher Site | Google Scholar
S. Sato, M. Hirakata, M. Kuwana et al., “Autoantibodies to a 140kd polypeptide, CADM-140, in Japanese patients with clinically amyopathic dermatomyositis,” Arthritis & Rheumatism, vol. 52, no. 5, pp. 1571–1576, 2005.
View at: Publisher Site | Google Scholar
R. Nakashima, Y. Imura, S. Kobayashi et al., “The RIG-I-like receptor IFIH1/MDA5 is a dermatomyositis-specific autoantigen identified by the anti-CADM-140 antibody,” Rheumatology (Oxford), vol. 49, no. 3, pp. 433–440, 2010.
View at: Publisher Site | Google Scholar
M. Yoneyama, K. Onomoto, M. Jogi, T. Akaboshi, and T. Fujita, “Viral RNA detection by RIG-I-like receptors,” Current Opinion in Immunology, vol. 32, pp. 48–53, 2015.
View at: Publisher Site | Google Scholar
D. J. B. Kurtzman and R. A. Vleugels, “Anti-melanoma differentiation-associated gene 5 (MDA5) dermatomyositis: a concise review with an emphasis on distinctive clinical features,” Journal of the American Academy of Dermatology, vol. 78, no. 4, pp. 776–785, 2018.
View at: Publisher Site | Google Scholar
F. Chen, D. Wang, X. Shu, R. Nakashima, and G. Wang, “Anti-MDA5 antibody is associated with A/SIP and decreased T cells in peripheral blood and predicts poor prognosis of ILD in Chinese patients with dermatomyositis,” Rheumatology International, vol. 32, no. 12, pp. 3909–3915, 2012.
View at: Publisher Site | Google Scholar
S. Moghadam-Kia, C. V. Oddis, S. Sato, M. Kuwana, and R. Aggarwal, “Antimelanoma differentiation-associated gene 5 antibody: expanding the clinical spectrum in north American patients with dermatomyositis,” The Journal of Rheumatology, vol. 44, no. 3, pp. 319–325, 2017.
View at: Publisher Site | Google Scholar
A. Ceribelli, M. Fredi, M. Taraborelli et al., “Prevalence and clinical significance of anti-MDA5 antibodies in European patients with polymyositis/dermatomyositis,” Clinical and Experimental Rheumatology, vol. 32, no. 6, pp. 891–897, 2014.
View at: Google Scholar
Z. Chen, W. Hu, Y. Wang, Z. Guo, L. Sun, and M. Kuwana, “Distinct profiles of myositis-specific autoantibodies in Chinese and Japanese patients with polymyositis/dermatomyositis,” Clinical Rheumatology, vol. 34, no. 9, pp. 1627–1631, 2015.
View at: Publisher Site | Google Scholar
A. Bohan and J. B. Peter, “Polymyositis and dermatomyositis (first of two parts),” The New England Journal of Medicine, vol. 292, no. 7, pp. 344–347, 1975.
View at: Publisher Site | Google Scholar
J. Rönnelid, S. Barbasso Helmers, H. Storfors et al., “Use of a commercial line blot assay as a screening test for autoantibodies in inflammatory myopathies,” vol. 9, no. 1, pp. 58–61, 2009.
View at: Publisher Site | Google Scholar
J. Lötvall, A. F. Hill, F. Hochberg et al., “Minimal experimental requirements for definition of extracellular vesicles and their functions: a position statement from the International Society for Extracellular Vesicles,” Journal of Extracellular Vesicles, vol. 3, no. 1, p. 26913, 2014.
View at: Publisher Site | Google Scholar
P. Landgraf, M. Rusu, R. Sheridan et al., “A mammalian microRNA expression atlas based on small RNA library sequencing,” Cell, vol. 129, no. 7, pp. 1401–1414, 2007.
View at: Publisher Site | Google Scholar
C. J. Creighton, J. G. Reid, and P. H. Gunaratne, “Expression profiling of microRNAs by deep sequencing,” Briefings in Bioinformatics, vol. 10, no. 5, pp. 490–497, 2009.
View at: Publisher Site | Google Scholar
M. D. Robinson, D. J. McCarthy, and G. K. Smyth, “edgeR: a Bioconductor package for differential expression analysis of digital gene expression data,” Bioinformatics, vol. 26, no. 1, pp. 139-140, 2010.
View at: Publisher Site | Google Scholar
X. Lv, Y. Li, H. Li et al., “Molecular function predictions and diagnostic value analysis of plasma exosomal miRNAs in Hirschsprung's disease,” Epigenomics, vol. 12, no. 5, pp. 409–422, 2020.
View at: Publisher Site | Google Scholar
Y. Ru, K. J. Kechris, B. Tabakoff et al., “The multiMiR R package and database: integration of microRNA-target interactions along with their disease and drug associations,” Nucleic Acids Research, vol. 42, no. 17, 2014.
View at: Publisher Site | Google Scholar
G. Yu, L.-G. Wang, Y. Han, and Q.-Y. He, “clusterProfiler: an R package for comparing biological themes among gene clusters,” OMICS: A Journal of Integrative Biology, vol. 16, no. 5, pp. 284–287, 2012.
View at: Publisher Site | Google Scholar
P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003.
View at: Publisher Site | Google Scholar
H. Valadi, K. Ekström, A. Bossios, M. Sjöstrand, J. J. Lee, and J. O. Lötvall, “Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells,” Nature Cell Biology, vol. 9, no. 6, pp. 654–659, 2007.
View at: Publisher Site | Google Scholar
J. Wahlgren, T. D. L. Karlson, M. Brisslert et al., “Plasma exosomes can deliver exogenous short interfering RNA to monocytes and lymphocytes,” Nucleic Acids Research, vol. 40, no. 17, 2012.
View at: Publisher Site | Google Scholar
K. Isoda, T. Kotani, T. Takeuchi et al., “Comparison of long-term prognosis and relapse of dermatomyositis complicated with interstitial pneumonia according to autoantibodies: anti-aminoacyl tRNA synthetase antibodies versus anti-melanoma differentiation-associated gene 5 antibody,” Rheumatology International, vol. 37, no. 8, pp. 1335–1340, 2017.
View at: Publisher Site | Google Scholar
M. Lake, G. George, and R. Summer, “Time to personalize the treatment of anti-MDA-5 associated lung disease,” Annals of the Rheumatic Diseases, vol. 78, no. 6, p. e52, 2019.
View at: Publisher Site | Google Scholar
Y. Wang, F. Zheng, G. Gao et al., “MiR-548a-3p regulates inflammatory response via TLR4/NF‐κB signaling pathway in rheumatoid arthritis,” Journal of Cellular Biochemistry, vol. 120, no. 2, pp. 1133–1140, 2019.
View at: Publisher Site | Google Scholar
F. Wang, L. Li, K. Piontek, M. Sakaguchi, and F. M. Selaru, “Exosome miR-335 as a novel therapeutic strategy in hepatocellular carcinoma,” Hepatology, vol. 67, no. 3, pp. 940–954, 2018.
View at: Publisher Site | Google Scholar
S. J. Szymura, G. M. Bernal, L. Wu et al., “DDX39B interacts with the pattern recognition receptor pathway to inhibit NF-κB and sensitize to alkylating chemotherapy,” BMC Biology, vol. 18, no. 1, 2020.
View at: Publisher Site | Google Scholar
S. Awasthi, B. Chakrapani, A. Mahesh, P. L. Chavali, S. Chavali, and A. Dhayalan, “DDX39B promotes translation through regulation of pre-ribosomal RNA levels,” RNA Biology, vol. 15, no. 9, pp. 1157–1166, 2018.
View at: Publisher Site | Google Scholar
M. Podralska, I. Ziółkowska-Suchanek, M. Żurawek et al., “Genetic variants in ATM, H2AFX and MRE11 genes and susceptibility to breast cancer in the polish population,” BMC Cancer, vol. 18, no. 1, p. 452, 2018.
View at: Publisher Site | Google Scholar
W. Wang, J. J. Qin, M. Rajaei et al., “Targeting MDM2 for novel molecular therapy: beyond oncology,” Medicinal Research Reviews, vol. 40, no. 3, pp. 856–880, 2020.
View at: Publisher Site | Google Scholar
J. S. Park, M. A. Lim, M. L. Cho et al., “p53 controls autoimmune arthritis via STAT-mediated regulation of the Th17 cell/Treg cell balance in mice,” Arthritis & Rheumatism, vol. 65, no. 4, pp. 949–959, 2013.
View at: Publisher Site | Google Scholar
Z. Zhong, H. Wu, W. Zhong, Q. Zhang, and Z. Yu, “Expression profiling and bioinformatics analysis of circulating microRNAs in patients with acute myocardial infarction,” Journal of Clinical Laboratory Analysis, vol. 34, no. 3, article e23099, 2020.
View at: Publisher Site | Google Scholar
X. Yang, H. Chen, Y. Chen et al., “Circulating miRNA expression profiling and target prediction in patients receiving dexmedetomidine,” Cellular Physiology and Biochemistry, vol. 50, no. 2, pp. 552–568, 2018.
View at: Publisher Site | Google Scholar
F. Liu, Z. Liang, J. Xu et al., “Activation of the wnt/beta-catenin signaling pathway in polymyositis, dermatomyositis and duchenne muscular dystrophy,” Journal of Clinical Neurology, vol. 12, no. 3, pp. 351–360, 2016.
View at: Publisher Site | Google Scholar
T. H. Sisson, P. J. Christensen, Y. Muraki et al., “Phosphodiesterase 4 inhibition reduces lung fibrosis following targeted type II alveolar epithelial cell injury,” Physiological Reports, vol. 6, no. 12, article e13753, 2018.
View at: Publisher Site | Google Scholar
S. Noda, H. Koike, S. Maeshima et al., “Transforming growth factor-beta signaling is upregulated in sporadic inclusion body myositis,” Muscle & Nerve, vol. 55, no. 5, pp. 741–747, 2017.
View at: Publisher Site | Google Scholar
Q. B. Xie, Y. Liang, M. Yang, Y. Yang, X. M. Cen, and G. Yin, “DEPTOR-mTOR signaling is critical for lipid metabolism and inflammation homeostasis of lymphocytes in human PBMC culture,” Journal of Immunology Research, vol. 2017, Article ID 5252840, 9 pages, 2017.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Danli Zhong et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

558

Downloads

1163

Citations