- Open Access
NovelmiRNA-25 inhibits AMPD2 in peripheral blood mononuclear cells of patients with systemic lupus erythematosus and represents a promising novel biomarker
Journal of Translational Medicinevolume 16, Article number: 370 (2018)
Systemic lupus erythematosus (SLE) is a multisystemic autoimmune disease with various clinical manifestations. MicroRNAs (miRNAs) and immunometabolism are recognized as key elements in SLE pathogenesis; however, the relationship between miRNAs in peripheral blood mononuclear cells (PBMCs) and metabolism in SLE remains unclear.
We detected PBMC miRNA and mRNA profiles from 3 pooled SLE patients and 3 healthy controls (HCs) using next-generation sequencing, predicted miRNA targets in dysregulated mRNAs, predicted functions and interactions of differentially expressed genes using bioinformatics analysis, validated candidate miRNAs using qRT-PCR, and investigated the association between the expression of candidate miRNAs and SLE clinical characteristics. Moreover, we validated the direct and transcriptional regulatory effect of NovelmiRNA-25 on adenosine monophosphate deaminase 2 (AMPD2) using a dual-luciferase reporter assay and western blot and confirmed AMPD2 mRNA and protein expression in PBMCs using qRT-PCR and western blot, respectively.
Multilayer integrative analysis of microRNA and mRNA regulation showed that 10 miRNAs were down-regulated and 19 miRNAs were up-regulated in SLE patient PBMCs compared with HCs. Bioinformatics analysis of regulatory networks between miRNAs and mRNAs showed that 19 miRNAs were related to metabolic processes. Two candidate miRNAs, NovelmiRNA-25 and miR-1273h-5p, which were significantly increased in the PBMCs of SLE patients (P < 0.05), represented diagnostic biomarkers with sensitivities of 94.74% and 89.47%, respectively (area under the curve = 0.574 and 0.788, respectively). NovelmiRNA-25 expression in PBMCs was associated with disease activity in SLE patients, in both active and stable groups (P < 0.05). NovelmiRNA-25 overexpression downregulated AMPD2 expression in HEK293T cells through direct targeting of the AMPD2 3ʹUTR (P < 0.01), while inhibition of NovelmiRNA-25 activity led to increased AMPD2 expression (P < 0.01). NovelmiRNA-25 overexpression also downregulated AMPD2 protein expression in HEK293T cells; AMPD2 protein expression in SLE patient PBMCs was decreased. Our results show that differentially expressed miRNAs play an important role in SLE.
Our data demonstrate a novel mechanism in SLE development that involves the targeting of AMPD2 expression by NovelmiRNA-25. miRNAs may serve as novel biomarkers for the diagnosis and evaluation of disease activity of SLE and represent potential therapeutic targets for this disease.
Systemic lupus erythematosus (SLE), an autoimmune disease mediated by pathogenic autoantibodies, presents with severe clinical manifestations such as nephritis, multisystem organ failure, or central nervous system disease . Prevalence rates of SLE worldwide vary from approximately 20 to 70 per 100,000 person-years, and incidence rates for women are approximately ten times higher than those for men . As the pathogenesis of SLE remains unclear, treatment, which mainly consists of glucocorticoid-antimalarial drugs and non-steroidal anti-inflammatory drugs (NSAIDs), often focuses on immune suppression to prevent and control disease progression, but flare still occurs [3, 4]. Therefore, the mortality rate of SLE remains high owing to the severity of the disease and the resulting organ damage .
Accumulating evidence has indicated that metabolic syndrome is a major co-morbidity of SLE . Of note, there is a growing understanding that cellular metabolic processes result in changes in inflammatory and immune responses in various immune cell populations , including macrophages, dendritic cells, neutrophils, and peripheral blood mononuclear cells (PBMCs) such as CD4+ T cells and B cells [7,8,9]. Specific metabolic processes are critical checkpoints of effector functions in the immune system, with common as well as cell-specific programs . Lupus-prone mice and SLE patients have been found to exhibit activated metabolism of CD4+ T cells; the use of metabolic inhibitors was found to effectively normalize these features and was associated with therapeutic effects . Furthermore, a proof-of-concept trial found a reduction in clinical flare-ups following metformin add-on treatment in mild or moderate SLE , thereby demonstrating that drugs involved in metabolic regulation may be used as add-on treatment in autoimmune diseases.
MicroRNAs (miRNAs) are small, noncoding RNAs that regulate numerous immunologic, inflammatory, and oncogenic pathways by modulating protein translation [12, 13]. miRNAs have been shown to play crucial roles in SLE pathogenesis , particularly lupus susceptibility, through their fundamental mechanism of post-transcriptional regulation of gene expression, suggesting their potential utility as diagnostic markers or therapeutic targets in this disease . Dysregulation of miRNAs has been associated with disease activity and major organ involvement in patients with SLE . For example, compared with that in healthy individuals, the expression of miRNA miR-21 is elevated in CD4+ T cells from lupus patients; overexpression of miR-21 suppressed the expression of PDCD4, a selective protein translation inhibitor, and regulated aberrant T cell responses . However, data on the relationship between miRNAs in PBMCs and metabolism in SLE has not yet been reported.
In this study, we first evaluated mRNA and miRNA expression profiles in PBMCs from Chinese SLE patients using next-generation sequencing (NGS). Secondly, integrating the results of our bioinformatic analysis of immunometabolism, we sought to identify potential miRNAs associated with the regulation of SLE progression and experimentally validated their differential expression between SLE patients and healthy controls (HCs). Lastly, through dual-luciferase reporter assays and western blotting, we identified the potential metabolism-associated targets that are regulated by miRNAs. The purpose of this study was to elucidate the mechanisms underlying SLE progression and investigate the potential applications of miRNAs in the treatment of this disease.
Subjects and study design
All 28 SLE patients admitted to the Department of Rheumatology and Nephrology, First and Second Affiliated Hospital of Wenzhou Medical University, between June 2016 and January 2018, were enrolled in this study. They were divided into two main independent SLE cohorts. The primary exploratory cohort consisted of three female patients, and the validation cohort comprised 25 patients. Moreover, PBMCs from 5 SLE patients and 5 HCs were selected randomly to validate the protein expression of miRNA target genes. All patients fulfilled the American College of Rheumatology (ACR) 1997 criteria for SLE . Disease activity was assessed according to the Systemic Lupus Erythematosus Disease Activity Index (SLEDAI)  at the time of blood collection. SLE patients with SLEDAI ≥ 5 were defined as having active disease, and those with SLEDAI < 5 were defined as having stable disease. The characteristics of all subjects are summarized in Additional file 3: Table S1. Twenty-eight age- and sex-matched healthy controls (HCs) without arthralgia, heart failure, renal failure, or autoimmune disease, and free from other inflammatory conditions, were recruited based on records available from the Wenzhou local blood bank. The research protocol was approved by the Medical Ethical Committees of the First and Second Affiliated Hospital of Wenzhou Medical University. All subjects who participated in this research provided written informed consent.
PBMC and RNA isolation
PBMCs were isolated from SLE patients and HCs using human peripheral blood lymphocyte separation medium (Tianjin Hao Yang Biological Manufacture, Tianjin, China) within 4 h of collection of the samples. Total RNA was extracted from each sample using TriZol Reagent (Invitrogen Life Technologies®, Grand Island, NY, USA). The isolated RNAs were digested by Dnase I (Invitrogen™, Waltham, MA, USA) to remove the residual DNA, and were then collected in 25 μL of DNase/RNase-free water. The concentration of RNA was quantified and qualified using a NanoDrop instrument (Thermo Fisher Scientific Inc., Waltham, MA, USA) and 1% agarose gel electrophoresis. Isolated RNA was stored at − 80 °C for use. The samples from the exploratory cohort were used for NGS analysis. Samples from the validation cohort were used for validation by quantitative reverse-transcription (qRT)-PCR.
All protocols, including NGS library preparation, were conducted according to the manufacturer’s instructions l (NEBNext® Ultra™ Directional RNA Library Prep Kit for Illumina®). Genewiz Ltd (https://www.genewiz.com.cn/) was used for all steps, from RNA library preparation to the generation of alignment files .
For small RNA sequencing, 3′ SR Adaptors for Illumina were ligated to small RNAs using the 3′ Ligation Enzyme. To prevent adaptor-dimer formation, excess 3′ SR Adaptors were hybridized with SR RT Primer for Illumina. Then, 5′ SR Adaptors for Illumina were ligated to the small RNAs using 5′ Ligation Enzyme, and first-strand cDNA was synthesized using ProtoScript II Reverse Transcriptase. Each sample was then amplified by PCR for 12 cycles using P5 and P7 primers, whose sequences annealed to flow cells to allow bridge PCR to be performed. Further, the P7 primer carried a six-base index allowing for multiplexing. PCR products of ~ 140 bp were recovered and cleaned up using PAGE, validated using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA), and quantified using a Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA, USA). Then, libraries with different indices were multiplexed and loaded on an Illumina HiSeq instrument according to the manufacturer’s instructions (Illumina, San Diego, CA, USA). Sequencing was carried out using a 1 × 50 bp single-end (SE) configuration; image analysis and base calling were conducted with HiSeq Control Software (HCS) + OLB + GAPipeline-1.6 (Illumina) on the HiSeq instrument.
For transcriptome sequencing, the poly(A) mRNA isolation was performed using NEBNext Poly(A) mRNA Magnetic Isolation Module (NEB). mRNA fragmentation and priming were performed using NEBNext First Strand Synthesis Reaction Buffer and NEBNext Random Primers. cDNAs were synthesized, purified, treated, amplified, validated, and quantified, respectively. Then, libraries with different indices were multiplexed and loaded on an Illumina HiSeq 2000 instrument according to the manufacturer’s instructions. Sequencing was carried out using a 2 × 150 bp paired-end (PE) configuration; image analysis and base calling were conducted with HiSeq Control Software (HCS) + OLB + GAPipeline-1.6 (Illumina) on the HiSeq instrument.
RNA-seq results were converted into read count values using Bioconductor Software DESeq (v.1.6.3) or edgeR (v.3.4.6) . Transcript expression levels were calculated using fragments per kilobase of transcript per million fragments mapped (FPKM) values. Fold changes in RNA expression between SLE patients and HCs were calculated from the signal values. Differential expression of miRNAs was considered significant at a fold change > 2 and P-value ≤ 0.05. Candidate circulating miRNAs were chosen for subsequent analysis based on a fold change > 2 and potential for association with target mRNAs. Similar to miRNAs, differential mRNA expression was considered significant at a fold change > 2 and false discovery rate (FDR) ≤ 0.05.
GO and KEGG pathway analysis
The functions of all differentially expressed miRNAs and mRNAs were investigated using Gene Ontology (GO) annotations and Kyoto Encyclopedia of Genes and Genome (KEGG) pathway analysis. Hierarchical clustering of the differentially expressed genes (DEGs) according to the biological process (BP), cellular component (CC), and molecular function (MF) categories was performed by GO analysis to elucidate genetic regulatory networks (http://www.geneontology.org). Pathway analysis using graphical diagrams was performed to explore DEG pathways using the KEGG database (http://www.genome.jp/kegg/). Significance was determined by P-value and FDR.
Target gene prediction and bioinformatics analysis of miRNAs and mRNAs
Based on the miRNA sequences and mRNA expression profiles in the exploratory cohort, potential targets sites for candidate miRNAs were predicted using the miRanda software (https://www.genewiz.com.cn/). Then, miRNAs were further filtered by combining prediction of miRNA-target pairs with mRNA profiling data. Subsequently, the functions of potential miRNAs were investigated by GO analysis using PANTHER software (http://www.pantherdb.org). By integrating the miRNA-mRNA regulatory network, a miRNA-gene regulatory network was constructed. To further explore the biological functions and build a network of the functions of miRNA target genes in SLE, the Thomson Reuters database (https://portal.genego.com)  was used. A coexpression network was constructed for the 29 candidate miRNAs and target genes using Cytoscape software. In addition, RNAhybrid software  was used to predict the resulting secondary structures of mRNA and miRNA interactions.
To identify miRNAs that play important roles in SLE development, the following criteria were applied (1) the key miRNAs and mRNAs should show significant expression changes between SLE patients and HCs; (2) they should exhibit a detectable negative relationship between the expression of the miRNA and the mRNA target genes; (3) the miRNA target genes should be associated with metabolic process in the BP category; (4) the binding energy between the miRNA and mRNA should be < − 31 kcal/mol (Table 1, Additional file 3: Table S7); (5) the 2–8 bp of the miRNA seed region should be strictly matched with target genes; (6) the expression level of target genes of miRNA should be log2foldchange > |2.9|.
qRT-PCR validation of identified miRNAs and miRNA target gene
From the 19 potential miRNAs identified in our comprehensive analysis, three representative miRNAs (miR-874-5p, miR-1273h-5p, and NovelmiRNA-25) were selected and validated by qRT-PCR. The miRNA expression levels in the PBMCs from 50 members of the validation cohort were detected using stem-loop qRT-PCR and the mirVana qRT-PCR miRNA Detection Kit (Ambion, Carlsbad, CA, USA). The miRNA target gene expression levels in the PBMCs from SLE patients of the validation cohort were detected using qRT-PCR. All qRT-PCR reactions were carried out on an Applied BioSystems 7500 Real-Time PCR system (Life Technologies). For each reaction, 1 μL of diluted cDNA was mixed with 5 μL of 2× SYBR Green Reaction Mix (DBI, Ludwigshafen, Germany). A final volume of 10 μL was achieved by the addition of 200 nM of forward and reverse primers. The conditions for PCR amplification were as follows: 95 °C for 3 min, followed by 39 cycles of 95 °C for 10 s and 58 °C for 30 s. The fluorescence signal was measured once every 1 °C. The specificity of the primer amplicons was tested by melting curve analysis. All samples were tested in triplicate. The data were analyzed using the comparative threshold cycle (Ct) method. U6/18S ribosomal RNA was used as a control, and the relative quantification of miRNAs in PBMCs was calculated using the following equation: Amount of target = 2−ΔCt, where ΔCt = Ct PBMC miRNAs − Ct U6/18S ribosomal RNA. The gene-specific primers used are listed in Additional file 3: Table S2.
Cell culture and transfection
HEK293T cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% fetal bovine serum (FBS) at 37 °C in a humidified atmosphere containing 5% CO2 (SANYO, MCO-175, Osaka prefecture, Japan). For cell transfection, HEK293T cells were seeded into 24-well plates at a density of 2 × 104 cells/well. Additionally, 60 nM of NovelmiRNA-25 mimic (5′-AGCAGGACGGUGGCCAUGG-3′; Guangzhou RiboBio Co., Ltd., Guangzhou, China) or a negative control (NC) mimic (Guangzhou RiboBio Co., Ltd.) was mixed with 1 µg of pmirGLO, pmirGLO-h_AMPD2 3′UTR, or pmirGLO-h_AMPD2 3′UTR-mut (GenBank: NM_001257360; Shanghai Bio-Link Co., Ltd.) using X-tremegene HP (Roche, Basel, Switzerland; No. 06366236001) according to the manufacturer’s protocols. Following this, the mixture was added to the 24-well plates to obtain a final concentration of 20 µM of NovelmiRNA-25 mimic or NC mimic. Plates were subsequently incubated at 37 °C for 48 h. Twenty-four wells were tested for each treatment to obtain the mean luciferase activity of each experiment. All experiments were performed in triplicate.
Luciferase reporter assay
Luciferase activity was measured 48 h after transfection using the Dual-Luciferase Reporter Assay System (Promega, Madison, WI, USA; E1910) in a Luminoskan Ascent (BioTek, Synergy HT, Vermont, USA) according to the manufacturer’s protocols. Firefly luciferase activity was normalized to Renilla luciferase activity to account for variations in transfection efficiency between experiments.
Forty-eight hours after NovelmiRNA-25 mimic transfection, HEK293T cells and PBMCs from 5 SLE patients and 5 HCs were lysed using protein lysis buffer (Beyotime Institute of Biotechnology, Beijing, China) supplemented with protease inhibitor cocktail (Pierce, Rockford, IL, USA) at 4 °C for 20 min. Protein samples were separated using 10% sodium dodecyl sulfate (SDS)-polyacrylamide gel electrophoresis and then electrophoretically transferred to polyvinylidene difluoride membranes (Millipore, Billerica, MA, USA). Antibodies against AMPD2 (ab31537, MULTI SCIENCES, Hangzhou, China) were diluted with primary antibody dilution buffer (Beyotime Institute of Biotechnology) at a 1:400 dilution, while those against GAPDH (GOOD HERE, Hangzhou, China) were diluted 1:1000. The membranes were then washed with TBST buffer five times for 5 min each and incubated with horseradish peroxidase (HRP)-conjugated goat anti-rabbit IgG secondary antibody (1:5000 dilution) (MULTI SCIENCES) for 1.5 h at 37 °C. Bands were detected using enhanced chemiluminescence and visualized with a Gel Doc 2000 (BioRad, Hercules, CA, USA).
Statistical analysis was performed with SPSS 22.0 software (SPSS, Inc., Chicago, IL, USA). Data are presented as mean ± standard deviation. Statistical significance between groups was determined by Student’s t test or Chi square test. A P-value < 0.05 was considered to represent statistically significant differences. Receiver operating characteristic (ROC) curve analysis, plotting the true positive rate (sensitivity) versus the false positive rate (1 − specificity) at various threshold settings, was performed for PBMC miRNAs, and the areas under curve (AUC) were calculated using Medcalc 15.2.2. The maximum of the sum of true positive rate and false positive rate was calculated, and cutoff value with higher specificity was selected. Expression graphs and ΔCt values were analyzed using GraphPad Prism version 5.04 software (GraphPad Software, La Jolla, CA, USA).
Overview of RNA-seq data for SLE patients and healthy controls
In order to investigate miRNA and mRNA interactions, the expression profiles of miRNAs and mRNAs in PBMCs from three SLE patients and three HCs (clinical characteristics of patients are shown in Additional file 3: Table S1) were analyzed. Using the Illumina HiSeq 2000 platform, an RNA-seq dataset of 616 million raw reads was obtained from the six samples (ranging from 92 million to 117 million). Using millions of short reads mapped to the human genome, the expression levels of known genes in each sample were quantified using the conventional FPKM parameter . Among annotated RefSeq coding genes, on average, 26,900 genes per sample exhibited detectable expression (RPKM > 0.05, roughly one mapped read for one 1-kb-exon gene; Additional file 3: Table S3). To evaluate the overall quality of our RNA-seq data, the RNA-seq data from multiple samples were compared using Pearson’s correlation and the Spearman and Kendall-tau rank correlation and found significant coefficients. On average, 91.8% (94.2 million per sample) of the short reads uniquely mapped to the human reference genome (hg19) or exon-junction sequences (Additional file 3: Table S3). No significant differences in the percentage of aligned reads between samples (SLE vs. HC: 90.0% vs. 93.2%) were observed, indicating the reliability of the sequencing results.
The application of the small RNA-seq protocol to these six samples (three SLE and three HCs) resulted in a small RNA-seq dataset of 1.88 million 20–25-nt short reads (mean of 313,000; Additional file 3: Table S3). On average, 48.7% (152,700 per sample) mapped to known miRNAs in the miRBase database  (Additional file 3: Table S3). Evaluation of repetitive reads was based on identical read sequences and compared to the corresponding genome location. Further data analysis revealed uniform distribution of the quality of the RNA-seq.
Analysis of SLE-related differentially expressed miRNAs and mRNAs
To identify SLE-related DEGs, we analyzed the log-scale-transformed RPKM data of protein-coding genes across two cohorts. For mRNA expression profiling, a P-value ≤ 0.05 and FDR ≤ 0.05  were set as threshold values, and 1781 genes were determined to exhibit at least a two-fold difference in expression between SLE-PBMCs and HC-PBMCs. Among these, 713 and 1068 genes were upregulated and downregulated, respectively, in SLE-PBMCs compared with levels in HC-PBMCs. A volcano plot revealed the variance in DEGs (Fig. 1a), and differentially expressed mRNAs were visualized by heat map (Fig. 1b).
Based on RPKM values, at a P ≤ 0.05 and fold change > 2 (FDR≤ 0.05), 66 miRNA genes showed significant differential expression: 34 were upregulated and 32 were downregulated in SLE-PBMCs (Fig. 2a, b and Additional file 3: Table S4). Differential expression between miRNAs in the PBMCs of SLE patients and those in HCs was visualized by volcano plot (Fig. 2a) and heat map (Fig. 2b). The fold changes of the significantly differentially expressed miRNAs are listed in Additional file 3: Table S4.
GO and KEGG pathway analysis
GO analysis showed that differentially expressed mRNAs and miRNAs were enriched in some common biological functions. In the MF category, enriched terms included binding, catalytic activity, and molecular transduction activity. In the BP category, enriched terms included biological regulation, cellular process, and metabolic process, and in the CC category, enriched terms included cell part, organelle part, and macromolecular complex. The GO terms of the differentially expressed mRNAs and miRNAs are shown in Figs. 1c and 2c, respectively.
The cellular processes of miRNA targets were further analyzed using the Thomson Reuters database. The results show that differentially expressed miRNAs were enriched for regulatory roles in various cellular and metabolic processes, as well as cell development and differentiation (Additional file 2A and Additional file 3: Table S5). Most relevant networks of differentially expressed miRNAs involved the movement of cell or subcellular component (91.5%) or cell morphogenesis (76.6%) as the top-scored network and cell adhesion (42.0%) or biological adhesion (42.0%) as the second-scored network (Additional file 2B and Additional file 3: Table S6). These results indicate that the differentially expressed mRNAs and target genes of differentially expressed miRNAs are involved in similar biological functions.
Correspondingly, the molecular function pathways were analyzed using KEGG pathway analysis of differentially expressed mRNAs and target genes of differentially expressed miRNAs (Figs. 1d and 2d). Differentially expressed mRNAs and target genes of differentially expressed miRNAs were both enriched in pathways such as the MAPK signaling pathway, NK cell-mediated cytotoxicity, and other receptor and chemokine signaling pathways. The molecular function of the associated MAPK signaling pathway involves regulation of cell differentiation and cellular metabolism . These functional analyses therefore identified relevant metabolic and cellular processes that are important in the development of autoimmune diseases.
Interactions of key differentially expressed mRNAs and miRNAs related to SLE
Accordingly, a three-step analysis was conducted. First, by integrating the RNA-seq data, the small RNA-seq data, and the prediction of miRNA-target mRNA interactions across the six samples, 29 differentially expressed miRNAs were identified, including 10 that were downregulated and 19 that were upregulated in SLE, that showed a significant negative correlation with the expression of the target mRNA (Fig. 3a and Table 1). Next, considering that miRNAs play an important role in metabolic processes involved in SLE occurrence and development, PANTHER enrichment analysis was used to further select, among the 29 potential miRNAs, 19 candidate differentially expressed miRNAs that were associated with metabolic process: miR-1273h-5p, miR-874-5p, miR-330-5p, miR-342-5p, miR-361-5p, miR-4470, miR-6716-5p, miR-7854-3p, miR-92b-5p, NovelmiRNA-25, NovelmiRNA-296, NovelmiRNA-30, NovelmiRNA-376, NovelmiRNA-426, NovelmiRNA-489, NovelmiRNA-72, NovelmiRNA-815, NovelmiRNA-862, and NovelmiRNA-92 (Fig. 3b and Additional file 3: Table S7). Cytoscape software was used to construct the complex regulatory network of the 29 candidate miRNAs. As shown in Additional file 1A, miRNAs and target genes interacted in the miRNA-target mRNA regulatory network. Third, after applying the criteria of binding energy < − 31 kcal/mol, strict matching of the 2–8 bp of the miRNA seed region with the 3′UTR region of the target gene, and miRNA target gene expression level of > |2.9| (log2foldchange), three miRNAs that were upregulated (NovelmiRNA-25, all_hsa-miR-92b-5p and miR-1273h-5p) and one miRNA that was downregulated (miR-874-5p) in SLE patients were selected. The expression of three candidate miRNAs of them was validated in the validation cohort using qRT-PCR. All three miRNAs showed similar expression patterns between the qRT-PCR and NGS analyses. However, only NovelmiRNA-25 (P = 0.018) and miR-1273h-5p (P = 0.001) showed significant differential expression between the PBMCs of patients with SLE and those OF healthy controls (Fig. 3c). Therefore, these two miRNAs were considered candidates for further analysis.
Correlation between candidate miRNA expression and SLE clinical characteristics
Using ROC curve analysis, the potential utility of NovelmiRNA-25 and miR-1273h-5p as new diagnostic biomarkers of SLE was explored further. As shown in Fig. 4a, the AUCs for NovelmiRNA-25 and miR-1273h-5p when distinguishing SLE patients from healthy controls were 0.707 and 0.788, respectively. The diagnostic sensitivities of NovelmiRNA-25 and miR-1273h-5p were 0.91 and 0.89, respectively, and the specificities were 0.60 and 0.56, respectively.
The relationship between candidate miRNA expression levels in PBMCs and the clinical characteristic of patients with SLE was investigated (Table 2). The expression of NovelmiRNA-25 in SLE patients was significantly correlated with leukocyte count (P = 0.045), neutrophil count (P = 0.002), albumin levels (P = 0.034), and low C3 (P = 0.019). Moreover, NovelmiRNA-25 expression in PBMCs was significantly associated with disease activity in patients with SLE, in both the active and stable groups (P < 0.05) (Fig. 4b).
Similarly, low expression of hsa-miR-1273h-5p in SLE patients was correlated with low C4 (P = 0.031) and lupus anticoagulant (P = 0.034). However, no significant association between hsa-miR-1273h-5p expression in PBMCs and disease activity in patients with SLE was observed, in either the active or the stable groups (P = 0.175) (Fig. 4b).
NovelmiRNA-25 directly targets AMPD2
Sequencing analysis showed that NovelmiRNA-25 is located on chromosome 16:72089401–72089448. Next, the potential molecular targets of NovelmiRNA-25 were predicted. The coexpression network of NovelmiRNA-25 revealed an association with AMPD2 (Fig. 5a). The binding energy of NovelmiRNA-25 with AMPD2 mRNA was − 30.1 kcal/mol. AMPD2 is involved in metabolic processes, including purine nucleobase metabolism (GO:0006144) and purine nucleotide metabolism (GO:0006163) (Additional file 3: Table S5). Moreover, the Thomson Reuters database showed that AMPD2 participates in various GO cellular processes, including cellular developmental process and cell differentiation (Additional file 2A and Additional file 3: Table S5). AMPD2 also plays a key role in the second-scored networks, including cell adhesion (42.0%), biological adhesion (42.0%), and the integrin–mediated signaling pathway (20.0%) (Additional file 2B and Additional file 3: Table S6). Therefore, we further explored the binding of NovelmiRNA-25 to the AMPD2 3′UTR. The sequences of both NovelmiRNA-25 and the AMPD2 3′UTR seeding region are highly conserved among primates and mammals (Fig. 5b).
To confirm further whether AMPD2 is a direct downstream target of NovelmiRNA-25, a dual-luciferase reporter assay was performed, in which the direct binding of NovelmiRNA-25 to the pmirGLO AMPD2 3′UTR transcript repressed luciferase activity. HEK293T cells, which do not express NovelmiRNA-25, were used, and the expression of barely detectable levels of NovelmiRNA-25 was confirmed (data not shown). We observed that cotransfection of the NovelmiRNA-25 mimic suppressed the luciferase activity of pmirGLO with the AMPD2 3′UTR by 49% (P < 0.01); moreover, mutation of the NovelmiRNA-25-binding region within the AMPD2 3′UTR abrogated the reduction in the luciferase activity of the vector by 32% (P < 0.01) (Fig. 5c), indicating specificity of the target site of AMPD2.
Next, western blotting was performed to determine whether NovelmiRNA-25 downregulates the protein expression of AMPD2. At 48 h post-transfection with the Novel-miRNA-25 mimic, AMPD2 protein was significantly decreased in HEK293T cells transfected with Novel-miRNA-25 mimic compared with levels in cells transfected with the internal control (P < 0.05, Fig. 5d), and the expression of AMPD2 protein was significantly decreased in PBMCs of SLE patients compared with that in PBMCs from HCs (Fig. 5e). This finding indicates that overexpression of NovelmiRNA-25 regulates AMPD2 expression at the protein level, further confirming that NovelmiRNA-25 directly targets AMPD2.
In the current study, we conducted a multi-tiered integrative analysis using miRNA and mRNA sequencing in PBMCs from subjects with SLE and HCs for the first time. We identified 29 differentially expressed miRNAs, of which 19 were upregulated and 10 downregulated in SLE, with detectable negative regulatory effects on the expression of their target genes in PBMCs from SLE patients. Unlike in previous individual microarray studies [26,27,28,29,30,31,32], we performed NGS analysis of miRNAs and mRNAs from PBMCs of SLE patients simultaneously. Consequently, we identified new differentially expressed miRNAs that were not previously reported in SLE patients.
Metabolomic studies have identified the potential contributions of certain metabolites, including lipids, amino acids, nucleic acids, and carbohydrates, to inflammation; these, therefore, represent promising targets for the diagnosis and disease monitoring of SLE [27,28,29]. Therefore, we selected miRNAs involved in metabolic processes to explore potential links between metabolism and SLE. Among the three selected miRNAs associated with metabolic processes, NovelmiRNA-25 and miR-1273h-5p, which were upregulated in SLE, were validated by qRT-PCR in the validation cohort.
miRNAs may be differentially expressed in the PBMCs of SLE patients and healthy individuals. However, various studies have reported conflicting data (studies are listed in Additional file 3: Table S8). Previous studies have found that miR-19b, miR-20a , miR-126 , miR-125b , miR-27a , miR-155, miR-181b, miR-17 , miR-7 , and miR-210  are dysregulated in the PBMCs of SLE patients resulting in the aberrant expression of essential target genes related to the pathogenesis and development of SLE, such as PTEN  and HIF-1α . In contrast to these previous studies, we identified two novel miRNAs, NovelmiRNA-25 and miR-1273h-5p, that may be involved in the development of SLE.
According to the known miRNA targets published in several databases, we found that the metabolism-associated target of NovelmiRNA-25 was AMPD2, which is involved in purine nucleotide catabolism [37, 38]. The miRNA-target gene interaction results revealed that AMPD2 participates in the top-scored networks, including essential biological and metabolomic processes. Among metabolic processes, AMPD2, involved in purine biosynthesis and metabolism, is essential for cellular energy homeostasis and nucleic acids synthesis . Based on the GeneCard database , we found that AMPD2 was specifically expressed in PBMCs, including CD8+ T cells, CD4+ T cells, B cells, and monocytes. Notably, these clusters of cells participate in the pathogenesis of SLE [41, 42]. AMPD2, which is mainly produced in the liver and spleen , is an enzyme that converts AMP to IMP. As adenosine and ATP are both synthesized from AMP, AMPD2 was involved in several different metabolic pathways by mediating changes in adenosine and ATP levels. Downregulated AMPD2 contributed to AMP accumulation, which mediated higher adenosine and ATP levels. On the one hand, adenosine is known to possess anti-inflammatory and analgesic activities . Koizumi et al.  reported that tofacitinib may increase the levels of adenosine in chondrocytes through the downregulation of AMPD2, indicating its anti-inflammatory activity. On the other hand, ATP is known to have pro-inflammatory effects and activates the inflammasome . In patients with RA, high concentrations of ATP were found in the synovial fluid , implying that ATP was highly involved in RA. Moreover, some studies have revealed the pathologic roles of increased ATP in the development of inflammatory intestinal disorders . Activated SLE increases the level of ATP, which possess pro-inflammation and immune affect, by regulating AMPD2 expression. Therefore, our results indicated that AMPD2 in SLE was more likely to be involved in the pro-inflammatory pathway by increasing ATP level. As AMPD2 is a metabolism-associated target gene of NovelmiRNA-25, these results indicate that NovelmiRNA-25 is a potential functional target in SLE.
Although upregulated miRNA miR-1273h-5p did not appear to be involved in disease flare-ups or organ injury in SLE, our findings did implicate miR-1273h-5p in the pathogenesis of SLE. miR-1273h-5p, with a length of 21 nt, is encoded by the intergenic sequence on chromosome 16 . miR-1273h-5p was predicted to target five genes associated with metabolic processes; these genes encoded a peroxidase, glycosyltransferase, ATP synthase and hydrolase, carbohydrate phosphatase, and exoribonuclease. Additionally, miR-1273h-5p has been found to be dysregulated in cancers and other autoimmune diseases. For example, the expression of miR-1273 was shown to be increased in the pancreas of a mouse model of pancreatic cancer .
Numerous clinical manifestations of SLE are similar to those of other illnesses, making the disease difficult to diagnose or distinguish from other autoimmune diseases. Furthermore, there is no single laboratory indication that can be used for a definitive diagnosis of SLE . In the present study, we found that miRNAs may be promising biomarkers for SLE diagnosis. An ROC curve was used to evaluate the expression of the two candidate miRNAs in the PBMCs of SLE patients, demonstrating that these miRNAs may possess clinical diagnostic value for distinguishing SLE patients from healthy individuals. Moreover, we explored the associations between these miRNAs and the clinical characteristics of SLE patients, including autoantibodies, complement, and other relevant features. The results showed that the expression of the two miRNAs was correlated with complement levels. Further analysis of miRNA levels between patients with active and stable SLE showed that the expression of NovelmiRNA-25 was positively correlated with the degree of SLE activity. This positive correlation suggests that NovelmiRNA-25 may be involved in progression and flare-ups as well as organ damage during SLE and indicates that it is a potential functional target in SLE. The results of a dual-luciferase reporter assay revealed that the overexpression of NovelmiRNA-25 resulted in the downregulation of AMPD2 gene expression in PBMCs through direct targeting of the AMPD2 3′UTR. In addition, the results of western blotting showed that the overexpression of NovelmiRNA-25 downregulated endogenous AMPD2 protein levels in a human cell line. Indicating the function of NovelmiRNA-25 in regulating AMPD2. Moreover, AMPD2 expression in the active SLE group was lower than that in the stable SLE group (Additional file 1B), and AMPD2 protein levels in SLE patients were decreased compared to those in HCs, indicating that AMPD2 was associated with SLE. This evidence further supports the theory that NovelmiRNA-25 is a promising target in the development of therapeutic strategies for SLE.
In summary, in the present work, miRNA and mRNA expression profiling in PBMCs from SLE patients and HCs identified miRNAs that are dysregulated in PBMCs from SLE patients. We further analyzed the regulatory networks of differentially expressed miRNAs and mRNAs and identified NovelmiRNA-25 and miR-1273h-5p as differentially expressed in SLE. The identification of these miRNAs may provide useful insights into the pathophysiology of SLE. Further experimental results point to a novel mechanism underlying SLE pathogenesis involving NovelmiRNA-25 and demonstrated the upregulation of this miRNA in SLE. Our study suggests that miRNAs represent novel diagnostic biomarkers, disease activity markers, and potential therapeutic targets for SLE. Further investigation is needed to clarify the roles of the identified miRNAs in the pathogenesis of SLE, with the aim of developing new methods for the diagnosis and treatment of SLE.
systemic lupus erythematosus
peripheral blood mononuclear cells
non-steroidal anti-inflammatory drugs
American College of Rheumatology
Systemic Lupus Erythematosus Disease Activity Index
Kyoto Encyclopedia of Genes and Genome
differentially expressed genes
receiver operating characteristic
areas under curve
Trager J, Ward MM. Mortality and causes of death in systemic lupus erythematosus. Curr Opin Rheumatol. 2001;13(5):345–51.
Pons-Estel GJ, Alarcon GS, Scofield L, Reinlib L, Cooper GS. Understanding the epidemiology and progression of systemic lupus erythematosus. Semin Arthritis Rheum. 2010;39(4):257–68.
Chen L, Morris DL, Vyse TJ. Genetic advances in systemic lupus erythematosus: an update. Curr Opin Rheumatol. 2017;29(5):423–33.
Dema B, Charles N. Advances in mechanisms of systemic lupus erythematosus. Discov Med. 2014;17(95):247–55.
Jang YN, Baik EJ. JAK-STAT pathway and myogenic differentiation. Jak-stat. 2013;2(2):e23282.
Cameron AR, Morrison VL, Levin D, Mohan M, Forteath C, et al. Anti-inflammatory effects of metformin irrespective of diabetes status. Circ Res. 2016;119(5):652–65.
Morel L. Immunometabolism in systemic lupus erythematosus. Nat Rev Rheumatol. 2017;13(5):280–90.
Pan W, Zhu S, Yuan M, Cui H, Wang L, et al. MicroRNA-21 and microRNA-148a contribute to DNA hypomethylation in lupus CD4+ T cells by directly and indirectly targeting DNA methyltransferase 1. J Immunol. 2010;184(12):6773–81.
Zhao X, Tang Y, Qu B, Cui H, Wang S, et al. MicroRNA-125a contributes to elevated inflammatory chemokine RANTES levels via targeting KLF13 in systemic lupus erythematosus. Arthritis Rheum. 2010;62(11):3425–35.
Yin Y, Choi SC, Xu Z, Perry DJ, Seay H, et al. Normalization of CD4+ T cell metabolism reverses lupus. Sci Transl Med. 2015;7(274):274218.
Wang H, Li T, Chen S, Gu Y, Ye S. Neutrophil extracellular trap mitochondrial DNA and its autoantibody in systemic lupus erythematosus and a proof-of-concept trial of metformin. Arthritis Rheumatol. 2015;67(12):3190–200.
Kuhn AN, Diken M, Kreiter S, Vallazza B, Tureci O, et al. Determinants of intracellular RNA pharmacokinetics: implications for RNA-based immunotherapeutics. RNA Biol. 2011;8(1):35–43.
Hamilton CE, Papavasiliou FN, Rosenberg BR. Diverse functions for DNA and RNA editing in the immune system. RNA Biol. 2010;7(2):220–8.
Zhang H, Huang X, Ye L, Guo G, Li X, et al. B Cell-related circulating MicroRNAs With the potential value of biomarkers in the differential diagnosis, and distinguishment between the disease activity and lupus nephritis for systemic lupus erythematosus. Front Immunol. 2018;9:1473.
Stypinska B, Paradowska-Gorycka A. Cytokines and MicroRNAs as candidate biomarkers for systemic lupus erythematosus. Int J Mol Sci. 2015;16(10):24194–218.
Carlsen AL, Schetter AJ, Nielsen CT, Lood C, Knudsen S, et al. Circulating microRNA expression profiles associated with systemic lupus erythematosus. Arthritis Rheum. 2013;65(5):1324–34.
Stagakis E, Bertsias G, Verginis P, Nakou M, Hatziapostolou M, et al. Identification of novel microRNA signatures linked to human lupus disease activity and pathogenesis: miR-21 regulates aberrant T cell responses through regulation of PDCD4 expression. Ann Rheum Dis. 2011;70(8):1496–506.
Hochberg MC. Updating the American college of rheumatology revised criteria for the classification of systemic lupus erythematosus. Arthritis Rheum. 1997;40(9):1725.
Bombardier C, Gladman DD, Urowitz MB, Caron D, Chang CH. Derivation of the SLEDAI. A disease activity index for lupus patients. The Committee on Prognosis Studies in SLE. Arthritis Rheum. 1992;35(6):630–40.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7(3):562–78.
Krishnan P, Ghosh S, Wang B, Li D, Narasimhan A, et al. Next generation sequencing profiling identifies miR-574-3p and miR-660-5p as potential novel prognostic markers for breast cancer. BMC Genomics. 2015;16:735.
Schuierer S, Tranchevent LC, Dengler U, Moreau Y. Large-scale benchmark of endeavour using MetaCore maps. Bioinformatics. 2010;26(15):1922–3.
Rehmsmeier M, Steffen P, Hochsmann M, Giegerich R. Fast and effective prediction of microRNA/target duplexes. RNA. 2004;10(10):1507–17.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8.
Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ. miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006;34(Database issue):D140–4.
Sampson JN, Boca SM, Moore SC, Heller R, Wren J. FWER and FDR control when testing multiple mediators. Bioinformatics. 2018. https://doi.org/10.1093/bioinformatics/bty064.
Arriens C, Mohan C. Systemic lupus erythematosus diagnostics in the ‘omics’ era. Int J Clin Rheumatol. 2013;8(6):671–87.
Menni C, Zierer J, Valdes AM, Spector TD. Mixing omics: combining genetics and metabolomics to study rheumatic diseases. Nat Rev Rheumatol. 2017;13(3):174–81.
Guma M, Tiziani S, Firestein GS. Metabolomics in rheumatic diseases: desperately seeking biomarkers. Nat Rev Rheumatol. 2016;12(5):269–81.
Teruel R, Perez-Sanchez C, Corral J, Herranz MT, Perez-Andreu V, et al. Identification of miRNAs as potential modulators of tissue factor expression in patients with systemic lupus erythematosus and antiphospholipid syndrome. J Thromb Haemost: JTH. 2011;9(10):1985–92.
Zhao S, Wang Y, Liang Y, Zhao M, Long H, et al. MicroRNA-126 regulates DNA methylation in CD4+ T cells and contributes to systemic lupus erythematosus by targeting DNA methyltransferase 1. Arthritis Rheum. 2011;63(5):1376–86.
Luo X, Zhang L, Li M, Zhang W, Leng X, et al. The role of miR-125b in T lymphocytes in the pathogenesis of systemic lupus erythematosus. Clin Exp Rheumatol. 2013;31(2):263–71.
Sourour SK, Aboelenein HR, Elemam NM, Abdelhamid AK, Salah S, et al. Unraveling the expression of microRNA-27a* & NKG2D in peripheral blood mononuclear cells and natural killer cells of pediatric systemic lupus erythematosus patients. Int J Rheum Dis. 2017;20(9):1237–46.
Kaga H, Komatsuda A, Omokawa A, Ito M, Teshima K, et al. Downregulated expression of miR-155, miR-17, and miR-181b, and upregulated expression of activation-induced cytidine deaminase and interferon-alpha in PBMCs from patients with SLE. Mod Rheumatol. 2015;25(6):865–70.
Li Y, Shen L, Ye Y, Zhang X. The expression and clinical significance of miR-7 in peripheral blood of patients with systemic lupus erythematosus. Basic Clin Med. 2014;34(1):53–7.
Huang Q, Chen SS, Li J, Tao SS, Wang M, et al. miR-210 expression in PBMCs from patients with systemic lupus erythematosus and rheumatoid arthritis. Ir J Med Sci. 2018;187(1):243–9.
Sabina RL, Morisaki T, Clarke P, Eddy R, Shows TB, et al. Characterization of the human and rat myoadenylate deaminase genes. J Biol Chem. 1990;265(16):9423–33.
Mahnke-Zizelman DK, Sabina RL. Cloning of human AMP deaminase isoform E cDNAs. Evidence for a third AMPD gene exhibiting alternatively spliced 5′-exons. J Biol Chem. 1992;267(29):20866–77.
Akizu N, Cantagrel V, Schroth J, Cai N, Vaux K, et al. AMPD2 regulates GTP synthesis and is mutated in a potentially treatable neurodegenerative brainstem disorder. Cell. 2013;154(3):505–17.
Fishilevich S, Nudel R, Rappaport N, Hadar R, Plaschkes I, et al. GeneHancer: genome-wide integration of enhancers and target genes in GeneCards. Database. 2017. https://doi.org/10.1093/database/bax028.
Li Q, Wu H, Liao W, Zhao M, Chan V, et al. A comprehensive review of immune-mediated dermatopathology in systemic lupus erythematosus. J Autoimmun. 2018. https://doi.org/10.1016/j.jaut.2018.07.007.
Suenaga R, Evans MJ, Mitamura K, Rider V, Abdou NI. Peripheral blood T cells and monocytes and B cell lines derived from patients with lupus express estrogen receptor transcripts similar to those of normal cells. J Rheumatol. 1998;25(7):1305–12.
Fredholm BB, Ijzerman AP, Jacobson KA, Linden J, Muller CE. International union of basic and clinical pharmacology. LXXXI. Nomenclature and classification of adenosine receptors—an update. Pharmacol Rev. 2011;63(1):1–34.
Koizumi H, Arito M, Endo W, Kurokawa MS, Okamoto K, et al. Effects of tofacitinib on nucleic acid metabolism in human articular chondrocytes. Mod Rheumatol. 2015;25(4):522–7.
Bours MJ, Swennen EL, Di Virgilio F, Cronstein BN, Dagnelie PC. Adenosine 5′-triphosphate and adenosine as endogenous signaling molecules in immunity and inflammation. Pharmacol Ther. 2006;112(2):358–404.
Ryan LM, Rachow JW, McCarty DJ. Synovial fluid ATP: a potential substrate for the production of inorganic pyrophosphate. J Rheumatol. 1991;18(5):716–20.
Gulbransen BD, Bashashati M, Hirota SA, Gui X, Roberts JA, et al. Activation of neuronal P2X7 receptor-pannexin-1 mediates death of enteric neurons during colitis. Nat Med. 2012;18(4):600–4.
Dang X, Qu X, Wang W, Liao C, Li Y, et al. Bioinformatic analysis of microRNA and mRNA regulation in peripheral blood mononuclear cells of patients with chronic obstructive pulmonary disease. Respir Res. 2017;18(1):4.
Rachagani S, Macha MA, Menning MS, Dey P, Pai P, et al. Changes in microRNA (miRNA) expression during pancreatic cancer development and progression in a genetically engineered KrasG12D;Pdx1-Cre mouse (KC) model. Oncotarget. 2015;6(37):40295–309.
Dai Y, Huang YS, Tang M, Lv TY, Hu CX, et al. Microarray analysis of microRNA expression in peripheral blood cells of systemic lupus erythematosus patients. Lupus. 2007;16(12):939–46.
GG: performed the experiments, analyzed and interpreted the data. GG, WH: drafted the manuscript. SX, YL: performed the experiments and statistical analysis. WK, LK, YS, LB, ZH, LQ, YS: acquired the data and provided material support. CC: analyzed and interpreted the data, revised the manuscript and finally approved the version of the manuscript for publication. XX: contributed to the conception and design of the study, analyzed and interpreted the data, supervised the study, provided the project funding, revised the manuscript and finally approved the version of the manuscript for publication. All authors read and approved the final manuscript.
This study was supported by Zhejiang Science and Technology Bureau (No. 2012C33126 and 2016C33214), the National Key Research and Development Program of China (No. 2016YFC0903902), the Zhejiang Provincial Natural Science Foundation of China (No. LY12H05003 and LY15H100004), and the Wenzhou Municipal Science and Technology Bureau (No. 2014Y0309). We thank the patients and healthy volunteers who participated in this study.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its Additional files.
Consent for publication
Ethics approval and consent to participate
The research protocol was approved by the Medical Ethical Committees of the First and Second Affiliated Hospital of Wenzhou Medical University. All subjects who participated in this research provided written informed consent.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.