Introduction
Systemic lupus erythematosus (SLE) is an obscure autoim mune disease characterized by its multisystem affection,1,2 with several factors including genetic determinants contribut ing to the disease development. (2 Epigenetic changes such as DNA methylation, histone modifications, and noncoding RNAs play an important role in lupus pathogenesis. (3 SLE patients have distinct miRNA signatures in peripheral blood mononuclear cells (PBMCs), plasma and different tissues com pared to healthy individuals. MiRNAs' dysregulation in SLE is associated with disease activity, organ system involvement, and autoantibody profiles, (4,5 hence miRNAs are considered as potential diagnostic and prognostic disease biomarkers. Among the investigated miRNAs in SLE is the miR-30e member of the miR-30 family. (6 The miR-30 family consists of six mature miRNA molecules (miR-30a, -30b, -30c-1, -30c-2, -30d, -30e). The encoding genes are located on chromosomes 1, 6, and 8. The seed sequence located near the 5'end is shared by these miRNAs but they have different compensatory sequences located near the 3'end. This enables miR-30 fam ily members to target different genes and pathways. (7 Previous studies have revealed that the miRNA-30 family has an impor tant regulatory role in cell differentiation, cellular senescence and apoptosis, (8-15 and are linked to several nervous, circula tory, gastrointestinal and respiratory disorders. (16-19 Moreover, the miR-30 family is involved in immune system aberrations including tumorigenesis, (16,20-23 and inflammatory disorders as Familial Mediterranean Fever (FMF). (24 MiR-30e in partic ular is a multifunctional miRNA25 with reports linking it to cancers, (26,27 and autoimmune diseases. (6,28-30 Interestingly, the mature products from both arms of the miR-30e precursor [miR-30e-5p and miR-30e-3p] were previously verified, (31 yet the 5' product was later shown to be the predominant one. (32 Previous reports demonstrated altered miR-30e-5p expression levels per se in several autoimmune diseases including SLE6,28 and myasthenia gravis. (29,30
The aim of this study was to evaluate the miR-30e-5p expression pattern in serum samples in an Egyptian SLE cohort and determine its potential association with various disease parameters as well as disease activity.
Subjects and methods
Study participants
Forty Egyptian SLE patients attending the Rheumatology and Rehabilitation Outpatient Clinic at Kasr Al Ainy Hospitals, Cairo University, from June 2018 to January 2019 were included in this study. Patients fulfilled the Systemic Lupus Interna tional Collaborating Clinics (SLICC) 2012 criteria. (33 Patients with other co-existing connective tissue, systemic autoim mune diseases, malignancy or infections at the time of sampling were excluded from the study. Demographic and cumulative clinical characteristics were recorded. Disease activity was evaluated using the Systemic Lupus Erythematosus Disease Activity Index-2K (SLEDAI-2K) at the time of assessment. (34 Thirty seven age-and gender-matched healthy volunteers were recruited as a control group. Prior to enroll ment in the study, an informed consent was obtained from all participants. The study adhered to the tenets of the Helsinki Declaration.
Serum separation and storage
Blood samples were dispensed into a sterile serum tube. Tubes were left to clot at room temperature. Samples were centrifuged within 2 h of collection for 10 min at 3000 rpm at 4 °C. The upper serum was placed in another tube and centrifuged for 10 min at 16,000 x g at 4°C. The supernatant was then transferred to another tube.
RNA extraction
Total RNA was extracted from serum using MiRNeasy Mini kit (Qiagen, Germany) following the manufacturer's instructions. The NanoDrop 1000 Spectrophotometer (NanoDrop Technolo gies) was used to determine the concentration & purity of extracted RNA.
Reverse transcription (RT)
TaqMan MicroRNA RT Kit (Applied Biosystems) was used for miR-30e-5p complementary DNA (cDNA) synthesis following the manufacturer's protocol. 5 µl of extracted RNA was used in a 15 µl reaction mix to synthesize the cDNA. The conditions for RT were as follows: 16°C for 30min at, 42°C for 30min, 85 °C for 5 min and 4°C on hold. The cDNA was stored at -20 ° C.
Quantitative Real Time Polymerase Chain Reaction (qRT-PCR) analysis
The expression level of miR-30e-5p was determined by TaqMan qRT-PCR on the step one real time PCR system using TaqMan MicroRNA Assay kit as well as TaqMan Universal Master Mix from Applied Biosystems following the manufacturer's protocol. The samples were normalized to miR-16 as previously reported. (6,35 qPCR was performed under the following settings: 10 min at 95 °C, followed by 40 cycles at 95 °C for 15 s and at 60 °C for 1min. Rela tive miR-30e-5p expression levels were calculated by 2-ΔΔCT as defined by Schmittgen et al., (36 where ACT of both patient and control groups = CT miR-30e gene-CT miR-16 gene. ΔΔCT values were obtained by subtraction of control CT values from CT values. Results were expressed as fold change.
Statistical methods
All statistical calculations were using the computer pro gram IBM SPSS (Statistical Package for the Social Science; IBM Corp, Armonk, NY, USA) release 22 for Microsoft Win dows. Numerical data were expressed in terms of mean and standard deviation or median and range. Qualitative data were expressed in terms of frequencies and percentages. Pearson's Chi-square test or Fisher's exact test was used to deter mine the relation between qualitative variables. Comparison between quantitative variables was done using the Student t-test for normally distributed data or the non-parametric Mann-Whitney test for non-normally distributed data. Correlations between quantitative variables were done using Spearman-rho method. Receiver operator characteristic (ROC) analysis was used to determine the optimum cut off value for the studied diagnostic marker. Accuracy was represented using the terms sensitivity and specificity. All tests were two-tailed. P-values less than 0.05 were considered as statistically significant.
Results
Demographic, clinical characteristics and laboratory parameters
This cross sectional study included a total of 77 female sub jects (40 SLE patients and 37 normal controls). The mean age of the patients was 28.7± 7.9 years and was comparable to the mean age of the controls 31.6 ±5.9 years (P = 0.07). The SLEDAI-2K at the last visit ranged from 0 to 32, with a median score of 6. The patients' cumulative clinical and serologic characteristics, and treatment received at the time of sampling are presented in Table 1.
Association of miR 30e-5p expression with demographic features, disease characteristics and activity
Among the demographic features studied, there was no significant correlation between the miR-30e-5p levels and patients' age (r=0.085, P = 0.593) or the age at onset (r = 0.19, P=0.228). The median miR-30e-5p expression level was 1.748 (0.223-20.485) among SLE patients as opposed to 0.877 (0.058-3.522) in the control group (P=0.02). Patients with serositis expressed higher levels of miR-30e-5p with a median value of 2.214 (0.788-5.682) compared to those without serositis 1.045 (0.223-20.485) (P = 0.013). There was no association between miR 30e-5p levels and any of the other included clin ical manifestations (Table 2), or disease the SLEDAI-2K score (r = -0.045, P = 0.778).
Sensitivity and specificity
Receiver operator characteristic (ROC) analysis was con structed to evaluate the efficiency of serum miR 30e-5p as a potential diagnostic biomarker for lupus patients. At a cutoff value ≥1.06, serum miR-30e-5p showed a sensitivity of 64.3% and a specificity of 61.5% with AUC 0.676 (95% CI: 0.559-0.794, P = 0.02) (Fig. 1).
Discussion
The obscure interplay of genetic, infectious and immunologic factors contributing to the pathogenesis and presentation of SLE further adds to its complexity and hampers targeted and personalized therapy. (37 Among the reported genetic deter minants implicated in the pathogenesis of SLE, aberrant miRNA-mediated regulation could be considered as a key player in the development of the disease and has been asso ciated with distinct disease features. (38 The aim of this study was to investigate the expression of miR-30e, a member of the miRNA 30 family characterized by its versatile functions25 including its role in autoimmune diseases. (6,28-30 To the best of our knowledge, data about miR-30e in SLE is limited to the studies in the Korean6 and Indian populations28 therefore, we evaluated miR-30e-5p expression in serum samples in an Egyptian SLE cohort and investigated its potential association with the disease characteristics. Our study demonstrated sig nificantly higher levels of miR-30e-5p expression among our patients as compared to the control group (P = 0.02). A finding that is similar to that detected by Kim et al.,6 who reported elevated plasma levels of miR-30e-5p in a Korean SLE cohort. Similarly Mishra et al. (28 demonstrated significantly increased miR-30e-5p expression (7 fold) in PBMCs of Indian SLE patients and serum of SLE mouse model with elevated expression of pro-inflammatory cytokines, type I interferons or type I interferons-inducible genes. MiR-30e-5p targets key negative regulators in innate immune responses including Toll-like receptor, Nod-like receptor and type I interferon signaling pathways. (39,40 Mishra et al. (28 reported a significant reduc tion of some targets of miR-30e-5p, such as TANK, TRIM38, SOCS1 and SOCS3 in SLE patients and socs1, socs3, atg5 and atg12 in SLE mice models, proposing that dysregulation of miR-30e expression may elevate innate immune responses, thus may be one of the underlying factors in autoimmune diseases including SLE. Moreover, the authors demonstrated that inhibition of miRNA-30e expression by introduction of stable inhibitor analogs including antagomir and lock nucleic acid inhibitor into SLE patients and SLE mouse model, respectively, significantly decreased SLE associated molecular signatures including type I interferon and proinflammatory cytokines. (28 In view of the role of other members of miR30 family in SLE per se, Liu et al. demonstrated that in B cell lines, miR-30a specifi cally binds the 3'-UTR of Lyn, a key negative regulator of B cell activation. In addition, miR-30a overexpression reduced levels of Lyn in B cells from SLE patients promoting B cell hyper-activity, proliferation and IgG production. (41 Additionally, Dai et al. demonstrated miR-30a up-regulation in renal biopsies of lupus nephritis patients. (42 It is of note that miR-30s tar get some critical components of calcium/calcineurin signaling essential for normal cell physiology including TRPC6, PPP3CA, PPP3CB, PPP3R1, and NFATC3. (43 Remarkably, studies suggest that abnormalities in the calcium-calcineurin pathway in SLE lead to breakage of B cell tolerance and induces differentiation of T cells toward the pathogenic Th17 subtype. (44
In the current study, there was no association between miRNA-30e-5p with any of the investigated demographic fea tures. Apart from serositis (P = 0.013), there was no association with any of the other studied clinical features or with dis ease activity. On the other hand, Kim et al. (6 reported that the expression levels of miR-30e-5p in plasma didn't corre late with any disease manifestations. Interestingly, the role of miR-30 family in serosal inflammation and fibrosis has been demonstrated previously, including the significant downregulation of miR-30e-5p as well as miR-30a-5p, miR-30d-5p in FMF patients than healthy controls, (24 and its association with peritoneal fibrosis following peritoneal hemodialysis. (45 Moreover, miRNAs have been shown to affect mesothelial cells and are potential risk factors of the development of mesothelioma46; with miRNA-30e-5p in particular demon strating distinct histopathologic mesothelial aberrations. (47
To the best of our knowledge, our study is the first to eval uate the expression pattern of miR-30e-5p in serum samples of Egyptian SLE patients. The major limitation of the current work is the rather small number of enrolled participants, nev ertheless, we found that the level of serum miR-30e-5p from lupus patients was significantly higher than that of controls. Significantly high levels of miR-30e-5p were demonstrated in patients with serositis. ROC curve analysis revealed the low accuracy of serum miR-30e-5p in differentiating our lupus patients from controls, thus, additional studies are required to evaluate the applicability of miR-30e-5p as a diagnostic biomarker. Given the fact that miRNA expression is influenced by genetic heterogeneity and exogenous influences as nutrition, exposure to certain environmental factors or medication, (48,49 larger studies including patients with a dif ferent genetic background will be required to assess these issues. Measuring the concentration of a particular miRNA at a specific disease state may miss the dynamically altered miRNA therefore, immunosuppression naïve patients need to be recruited to demonstrate the expression of miRNA-30e-5p at the earliest stages of the disease and then throughout the disease course and in correlation with disease severity. Much remains to be established regarding the role of miR-30s includ ing miRNA-30e-5p as diagnostic biomarkers as well as their role in SLE pathogenesis and establishing organ compromise. Future research addressing these issues might contribute to better understanding of the disease process and may pave the way toward identifying novel therapeutic targets.