Front. Med.Frontiers in MedicineFront. Med.2296-858XFrontiers Media S.A.10.3389/fmed.2024.1383186MedicineOriginal ResearchConstruction of competing endogenous RNA networks in systemic lupus erythematosus by integrated analysisHeJuanjuan12†DaiYunfeng12†LiuJianwen12LinHe12GaoFei12ChenZhihan12*WuYanfang12*1Fujian Medical University Shengli Clinical Medical College, Fuzhou, China2Department of Rheumatology, Fujian Provincial Hospital, Fuzhou, China
Edited by: Zhiming Lin, Third Affiliated Hospital of Sun Yat-sen University, China
Reviewed by: Runyue Huang, Guangzhou University of Chinese Medicine, China
Chu Yong Liang, Victoria University, Australia
Hong-yan Wen, Second Hospital of Shanxi Medical University, China
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Objective
Systemic lupus erythematosus (SLE) is a disease characterised by immune inflammation and damage to multiple organs. Recent investigations have linked competing endogenous RNAs (ceRNAs) to lupus. However, the exact mechanism through which the ceRNAs network affects SLE is still unclear. This study aims to investigate the regulatory functions of the ceRNAs network, which are important pathways that control the pathophysiological processes of SLE.
Methods
CircRNA microarray for our tested assays were derived from bone marrow samples from three healthy individuals and three SLE patients in our hospital. The other sequencing data of circRNA, miRNA and mRNA were obtained from Gene Expression Omnibus (GEO) datasets. Using the limma package of R program, the differential expression of mRNA and miRNA in the GEO database was discovered. Then predicted miRNA-mRNA and circRNA-miRNA were established using miRMap, miRanda, miRDB, TargetScan, and miTarBase. CircRNA-miRNA-mRNA ceRNA network was constructed using Cytoscape, and hub genes were screened using a protein-protein interaction network. Immune infiltration analysis of the hub gene was also performed by CIBERSORT and GSEA.
Results
230 overlapped circRNAs, 86 DEmiRNAs and 2083 DEmRNAs were identified in SLE patients as compared to healthy controls. We constructed a circRNA–miRNA–mRNA ceRNAs network contained 11 overlapped circRNAs, 9 miRNAs and 51 mRNAs. ESR1 and SIRT1 were the most frequently associated protein-protein interactions in the PPI network. KEGG analysis showed that DEGs was enriched in FoxO signaling pathway as well as lipids and atherosclerosis. We constructed a novel circRNA-miRNA-mRNA ceRNA network (HSA circ 0000345- HSA miR-22-3-P-ESR1/SIRT1) that may have a major impact on SLE.
Conclusion
Through this bioinformatics and integrated analysis, we suggest a regulatory role for ceRNA network in the pathogenesis and treatment of SLE.
systemic lupus erythematosusdifferentially expressed genescompeting endogenous RNAsbioinformatics analysisenrichment analysis2023J0111992023Y9305Natural Science Foundation of FujianJoint Funds for the Innovation of Science and Technology, Fujian provincesection-at-acceptanceRheumatologyIntroduction
Systemic lupus erythematosus (SLE) is a complex autoimmune connective tissue disease that can lead to functional impairment of multiple organs and increased morbidity and mortality (1, 2). Many studies have confirmed the involvement of non-coding RNAs in the pathogenesis of lupus (3), while the exact pathogenesis is not fully understood. Thus, exploring the molecular features and mechanisms of non-coding RNA in SLE is crucial for better insight and treatment of the disease.
Competing endogenous RNAs (ceRNAs) are molecules of both coding and non-coding RNAs that can be targeted by the same miRNAs in the right context and can indirectly regulate each other by competing for them, leading to additional post-transcriptional regulatory layers where non-coding RNAs can find new meaning. miRNA-mediated interactions between different types of RNA molecules have been observed in many different contexts (4). Circular RNAs (circRNAs) lack 5′- and 3′-terminal poly (A) tails, making them less susceptible to nucleic acid exonuclease-driven degradation, resulting in a longer median half-life (5). Therefore, circRNAs are more likely to function as sponges and participate in the competing endogenous RNA (ceRNA) network (6). Researchers have experimentally tested instances of ceRNA crosstalk in a large number of situations and observed in normal and pathological contexts, demonstrating the widespread involvement of this mechanism in brain regenerative mechanisms (7), neuronal and muscle developmental processes (8) and cellular differentiation (9), as well as in diseases which highly complex gene regulatory circuits are most affected by perturbations. Researchers strive to model these regulatory networks to predict and understand how modifications alter the dynamic balance between molecules, leading to the onset of cancer and its progression, cardiovascular problems and neurodegenerative diseases, as well as other pathologies, such as those related to immune and autoimmune responses and degenerative physical conditions. Recently, several studies reported an association between ceRNA network and lupus (10, 11). The development of new computational tools that provide hypothetical predictions of ceRNA interactions facilitates insights into lupus pathogenesis and the discovery of novel therapeutic targets.
CircRNA, miRNA, and mRNA data can be obtained from the Gene Expression Omnibus (GEO) database (12), making it a valuable tool for biological discovery and data mining. Integrative ceRNA regulation circuits can be built to investigate more precise prognostic markers using these publicly available databases. In our study, we constructed ceRNAs networks of lupus by taking the intersection of our self tested circRNA microarray data with datasets from the GEO database. So far, such a study is unique.
In this research, we performed differential expression analysis using circRNA microarray of bone marrow samples from three healthy participants and three patients in our department and circRNA datasets in the GEO dataset. We then used online tools to identify the relationship between miRNA-circRNA and miRNA-mRNA to construct ceRNA network. By using enrichment function analysis, the biological significance of mRNA differential expression was determined. The protein-protein interaction (PPI) network (13) was used to search hub genes. Immune infiltration analysis of the hub gene was also performed by CIBERSORT and GSEA.
MethodsBone marrow samples collection
Bone marrow samples were obtained from 3 healthy participants and 3 patients with systemic lupus erythematosus (SLE) at Fujian Provincial Hospital. All participants, both healthy and those with SLE, were of Han Chinese origin. The inclusion criteria required patients with untreated first-onset SLE and healthy subjects. SLE diagnosed according to the 1997 updated American College of Rheumatology (ACR) classification criteria (14). This research received approval from the Ethics Committee of the Fujian Provincial Committee. Informed consent was obtained from all participants.
Data collection
We searched for SLE-associated aberrantly expressed RNA microarray data from the GEO database.1 The selection criteria were as follows: (a) The dataset included lupus patients and healthy controls. (b) Arrays contained at least 3 samples from patients and 3 samples from the normal group; (c) Uploaded data were available for analysis. Finally, GSE84655 (circRNA, 6 SLE vs. 3 control), GSE175840 (miRNA, 5 SLE vs. 5 control) and GSE175839 (mRNA, 5 SLE vs. 5 control) were included in the study. The downloaded files were calibrated, standardized, and log2 transformed using R software (15).
Identification of the differentially expressed circRNAs, miRNAs and mRNAs
Batch effects were normalised using the sva package. Differential analysis was performed using the limma package. DEcircRNAs were identified using a threshold of |log2FC| over 1.2 and an adjusted p-value less than 0.05. Also, the ballgown package was employed to identify DEmiRNAs and DEmRNAs. The thresholds set were an adjusted p-value less than 0.05 and |log2FC| above 1.2 for DEmiRNAs and DEmRNAs.
ceRNA network
The miRNAs corresponding to circRNAs were predicted by StarBase (16). To obtain the paired miRNA_mRNAs, miRMap (17), miRanda (18), miRDB (19), TargetScan (20), and miTarBase (21) were used to find possible targeted mRNAs of miRNA. Subsequently, the overlapping predictions between the two programs were considered effective target pairs. ceRNA networks were constructed in Cytoscape (version 3.7) (22).
Gene function enrichment analysis
The database for annotation, visualization, and integrated discovery (DAVID) (23) was employed to determine the biological functions of DEGs. The biological process (BP), molecular function (MF), and cellular component (CC) were applied in GO enrichment analyses.2 The Kyoto Encyclopedia of Genes and Genomes (KEGG) serves as a database resource for comprehending biological systems and overarching functions (24). The cut-off criteria were the false discovery rate < 0.1 and adjusted p-value < 0.05.
Formation of protein-protein interaction network and identification of hub genes
A PPI network of the DEmRNAs was established using the Search Tool for the Retrieval of Interacting Genes (STRING).3 Required Confidence (combined score) > 0.4 was selected as the threshold for PPI. The PPI network was visualized using Cytoscape. Molecular Complex Detection (MCODE) app of Cytoscape was used to determine hub genes.
Immune infiltration analysis of hub-gene
CIBERSORT, an universal deconvolution algorithm, was employed to examine the immune cell subset proportions using RNA expression profiles (25). Herein, we obtained a matrix of 22 immune cell subsets in GSE50772 via the R package of CIBERSORT. A bar plot displayed the percent of each immune cell and compare infiltrating levels between SLE patients and controls using the “ggplot2 (v3.3.0)” package, and the “corrplot (v0.90)” package was utilized to depict the relationship of immune cell subsets. p < 0.05 was determined to be statistically significant.
Gene set enrichment analysis
GSEA software (version 3.0) was obtained from the GSEA website.4 The Molecular Signatures Database5 was used to assess the pathways and molecular mechanisms involved in the hub genes. p < 0.05 was determined to be statistically significant.
ResultIdentification of DEcircRNAs in SLE
First, we analyzed 2 profiles of circRNAs including our self data and GSE84655, basic information was shown in Figure 1A. We detected 2,668 different expressed circRNAs (1,602 up-regulated, 1,066 down-regulated) in lupus samples using our tested bone marrow circRNA microarray (Figure 1B). Figure 1C showed the heatmap of top 20 differentially expressed circRNAs. To refine the identification further, a search for DEcircRNAs was conducted in GSE84655 and 1739 DEcircRNAs (847 up-regulated, 892 down-regulated) were found (Figure 1D). Subsequently, the top 20 circRNAs displaying the most significant variations were selected and are exhibited in Figure 1E. Ultimately, we integrated 230 overlapped circRNAs in SLE patients as compared to healthy controls (66 down regulated circRNAs Figure 1F and 164 up regulated Figure 1G).
Identification of DEcircRNAs in SLE (A) Basic information of circRNA datasets. (B) Volcano plots for different expression of circRNAs of bone-marrow in SLE. (C) Heatmap plot showing the top 20 DEcircRNAs in our dataset. (D) Volcano plots for DEcircRNAs of GSE84655. (E) Heatmap plot showing the top 20 DEcircRNAs of GSE84655. (F,G) The Venn plots of DEcircRNAs of two datasets.
Determination of lupus-associated DEmRNAs and DEmiRNA
To explore the potential functions of mRNAs/miRNAs/circRNAs in SLE, The mRNA (GSE175839) and miRNA (GSE175840) datasets were further analyzed. In SLE group, 86 DEmiRNAs (15 up-regulated, 71 down-regulated), and 2083 DEmRNAs (167 up-regulated, 1916 down-regulated) were discovered based on predetermined thresholds (Figure 2A). The top 20 DEmiRNA and DEmRNAs revealed in Figures 2B,C.
Determination of lupus-associated DEmRNAs and DEmiRNA. (A) Basic information of miRNAs and mRNA datasets. (B,C) Hotmap of the 20 top DEmiRNAs and DEmRNAs. (D,E) The Venn plots of predicted miRNAs based on DEmRNAs of GSE175839.
In order to find out the targeted miRNAs corresponding to GSE175839, DEmRNAs were searched in the miRMap, miRanda, miRDB, TargetScan and miTarBase. Based on the principle that paired miRNA-mRNA presented at least four databases, we confirmed 34 pairs of down-regulated miRNA and up-regulated interacted mRNA (Figure 2D) and 115 pairs of up-regulated miRNA and down-interacted mRNA (Figure 2E).
ceRNA network in SLE
Using the overlapped DEcircRNAs discovered in Figure 1, the 156 pairs of miRNAs that were down-regulated and the 1,180 pairs of miRNAs that were up-regulated were extracted from Starbase.
As shown in the flowchart (Figure 3), overlapped DEcircRNAs of (Figure 1) predicted cirRNA-miRNA, mRNA-miRNAs predicted by DEmRNAs of GSE175839 and the DEmiRNAs of GSE175840 were taken interacted, and identified the experimentally strongly supported SLE specific differentially expressed of the circRNA-miRNA-mRNA network. This network contained 11 circRNAs, 9 miRNAs and 51 mRNAs (Table 1). Then, we visualized these candidate genes and constructed the ceRNA network using cytoscapeV3.7 (Figure 4).
A flowchart of DEGs analysis.
Different expression of circRNAs-miRNAs-mRNAs in SLE.
Construction of CeRNAs networks (A) CircRNA(down)-miRNA(up)-mRNA(down) regulatory network (B) CircRNA(up)-miRNA(down)-mRNA(up) regulatory network.
Functional and pathway enrichment analysis for the targeted genes in the ceRNA network
We used DAVID to conduct functional and pathway enrichment analyses to gain more insight into the biological roles of the 51 DEmRNAs implicated in SLE. The top 15 KEGG pathways and enriched GO terms were exhibited in Figure 5. The GO analysis revealed that DEGs were primarily enriched in the following categories: “chromatin organization,” “chromosomal organization,” and “the negative regulation of fat cells” (BPs); “intracellular membrane-bounding organelle,” “nuclear lumen,” and “nucleotic” (CCs); and “kinase binding,” “protein kinase binding,” “RNA polymerase II-specific DNA binding transcription activator activity” (MFs). KEGG pathway analysis presented that DEGs of ceRNA network of SLE were enriched in malignancies, infection, Lipid and atherosclerosis, apoptosis, endocrine resistance, FoxO signaling, HIF-1 signaling, Estrogen signaling and thyroid hormone signaling (Figure 5D).
GO term and KEGG pathways enrichment analysis (A) Biological process (BP), (B) cellular component (CC), (C) molecular function (MF) of the GO annotation (D) KEGG pathways of DEmRNA.
Identification of hub genes
After obtaining the target mRNA of ceRNA network in SLE, 51 DEmRNAs were imported into the PPI network composed of 51nodes and 22edges (Figure 6A). Following the identification of the vital functions of hub genes in the network, we used MOCDE plugin of cytoscape for further analysis, and 9 hub genes were explored (HMOX1, SIRT1, ESR1, AKT3, CDKN18, KDM3A, CHD9, NCOA1, ZNF609).
PPI network and hubgenes of targeted DEmRNAs in SLE (A) PPI network of DEmRNAs (B) The top 9 hubgenes analysed by Cytoscape.
Immune infiltration characteristics of hub gene SIRT1and ESR1
When screening hub genes (Figure 6B), ESR1 and SIRT1 were the most frequently associated protein-protein interactions in the PPI network. To further investigate the mechanism of SIRT1 and ESR1 involvement in lupus pathogenesis, we analyzed SIRT1 and ESR1 using GSEA and also found that SIRT1 and ESR1 were involved in many biological processes of immune cells (Figures 7A,B), which had been demonstrated in SLE. Immune cell infiltration also confirmed the presence of statistically significant differences in CD4+ memory T cells, naive B cells and M0 Macrophages between SLE patients and controls (Figure 7C). Spearman correlation analysis was subsequently performed to demonstrate the relationship between DEOSRGs and various immune cell subsets. As shown in Figure 7D, SIRT1 was positively correlated with CD4+ memory T cells (r = 0.733, p = 0.021) and naive B cells (r = 0.672, p = 0.039), while it was negatively correlated with M0 Macrophages (r = −0.661, p = 0.044). Meanwhile, the immune infiltration analysis (Figure 7E) indicated that ESR1 had a positive correlation with CD4+ memory T cells (r = 0.87, p = 0.0027) and NK cells (r = 0.77, p = 0.014), but a negative correlation with regulatory T cells (r = −0.79, p = 0.0098), eosinophils (r = −0.83, p = 0.0056), and M0 macrophages (r = −0.89, p = 0.0013). Using the CIBERSORT method, we evaluated the immune infiltration features between SLE patients and controls. The bar plot depicted the relative percent of diverse immune cell subsets in each sample (Figure 7F). These results revealed that SIRT1 and ESR1 may be crucial to the progression in SLE.
Immune infiltration characteristics between SLE patients and controls. (A,B) GSEA analysis for SIRT1 and ESR1 in SLE group. (C) Infiltrating difference of immune cells between SLE patients and controls in the box plot. (D,E) Spearman correlation between immune cell subsets for SIRT1 and ESR1. The color of dots denotes the P value. The size of dots denotes the strength of correlation. (F) The relative percent of immune cell subsets in a bar plot. Statistical significance at the level of * <0.05, ** <0.01, and *** <0.001. SLE: Systemic lupus erythematosus; CON: Controls.
Discussion
The expression profiles of non-coding RNAs (ncRNAs), which include miRNAs, long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs), have attracted increasing attention in various diseases due to the rapid growth of bioinformatics (26, 27). Compared with healthy individuals, most circRNAs were found to be abundant in peripheral blood, serum, T cells, peripheral blood mononuclear cells and kidney tissue, indicating their possible roles. In summary, circRNA dysregulation may have implications for key molecular processes involved in the pathogenesis of systemic lupus erythematosus (28). Therefore, a more comprehensive understanding of the importance of circRNAs in SLE is required. We used microarray data analysis to construct the circRNA-miRNA-mRNA network. The mRNAs in the network were used to identify significantly altered GO terms and the KEGG signaling pathway. The hub genes were then identified using the PPI network. In order to open up a new perspective on the therapy of SLE, certain effective substances have also been discovered.
Based on the results of functional enrichment analysis, we hypothesized that in SLE, the aberrant circRNAs utilize protein and RNA binding to regulate chromosome organization, chromatin organization, and negatively regulate adipocyte differentiation, as well as positively regulate colorectal migration of blood vessels. According to KEGG signaling pathway analysis, it is linked to cancer-related signaling pathways. Since the immunological and genetic pathways of SLE pathogenesis and the use of immunosuppressive drugs (ISDs) in treatment may increase the likelihood of such changes, the relevant literature has shown (29) that the risk of malignancy in SLE is considerable. In KEGG pathway analysis, it can be found that ESR1 is mainly related to endocrine resistance, the cancer signaling pathway, the estrogen signaling pathway, and Epstein-Barr virus infection, while SIRT1 is mainly related to the FoxO signaling pathway and Salmonella infection. Salmonella infection disrupts SIRT1/AMPK checkpoint control of mTOR to impair autophagy (30), and alterations of autophagy contribute to the progression of various autoimmune diseases, including systemic lupus erythematosus (SLE) (31). In the FoxO signaling pathway, such observations with Foxo3a in helper T cells imply disease relevance in SLE (32). However, the SIRT1-regulated signaling pathway of these hub genes (ESR1) in autoimmune diseases has not been reported. Whether it is related to the pathogenesis of SLE needs to be confirmed by further experiments.
SLE is more prevalent in women of reproductive age between the ages of 15 and 40, with a male-to-female incidence ratio of 1: 9. As a result, it has been suggested that estrogen may play a role in the development of SLE. In lymphoid nuclei, estrogen primarily interacts with estrogen receptor types 1 and 2 (ESR1 and ESR2, respectively) to carry out biological actions (33). Further evidence that ESR1, not ESR2, produces the lupus phenotype can be found in male mice lacking functioning ESR1. These mice are resistant to developing the lupus phenotype (34). Recent studies demonstrating that ESR1 promotes SLE in F1-generation female mice of the lupus mouse model (35) lend weight to this idea. In a previous investigation, Lee et al. (36) found that ESR1 gene polymorphisms were related to the age at which SLE first manifested in Korean patients; however, this study only included female participants. ESR1 rs2234693 and rs9340799 were revealed to be strongly linked with SLE susceptibility by Wang et al. (37). When Zhou et al. (38) compared the frequencies of C allele rs2234693 and G allele rs9340799 in the ESR1 gene between SLE patients and non-SLE controls, they found that the frequency of the C allele was considerably higher. Faslodex, an ESR1 antagonist, significantly decreased SLEDAI scores in SLE patients in a limited clinical trial (39). These findings support the idea that estrogen plays a part in the development of SLE.
SIRT1 is an NAD+-dependent monomeric protein that plays a significant role in important cellular activities, including cell differentiation, apoptosis, metabolism, aging, and immunological response. In Consiglio’s research, the SIRT1 promoter variant rs3758391 was found to be connected with SLE incidence, and the rs3758391 T allele may be linked to higher SLEDAI and lupus nephritis scores (40). Anti-dsDNA antibodies are well documented to fluctuate with SLE activity and to be closely associated with the symptoms of severe lupus (41, 42). Olivares et al. found that serum anti-dsDNA antibody levels in LN patients were significantly correlated with urinary SIRT1 mRNA levels, which is a valuable marker of renal injury (43). Systemic lupus erythematosus (SLE) patients had higher levels of the plasma protein SIRT1 (SIRT1) compared to healthy controls, and there was a significant correlation between the plasma SIRT1 concentration and disease activity, according to research by Yang et al. (44).
The final circRNA-miRNA-mRNA network contained HSA circ 0000345 and HSA miR-22-3p, which we identified. Only pertinent literature has demonstrated that there is currently no relationship between this circRNA and autoimmune disorders (45, 46). By suppressing HSA circ 0000345, HSA circ 0000345 plays a critical role in the treatment of atherosclerosis. The involvement of HSA circ 0000345 in endothelial cell damage is strongly related to the onset and progression of atherosclerosis. Down-regulation of HSA circ 0000345 has been demonstrated in the literature (47) to have an anti-tumor impact. However, there are no studies indicating that abnormal expression of HSA circ 0000345 contributes to the development of SLE. Future studies should verify the expression status of HSA circ 0000345 in SLE patients and its function. Additionally, Previous studies have demonstrated that miR-22-3p expression levels are increased in the peripheral blood of patients with lupus (48, 49), which is consistent with our findings. In contrast, compared to controls, SLE with lupus nephritis (LN) patients exhibited a downregulation of miR-22-3p (50). The disparate results may be attributed to the varying cell subtypes examined. Further experimentation is required to confirm its status and function.
Nevertheless, it is necessary to acknowledge the limitations of this work. First, the sample size of bone marrow for high-throughput sequencing was limited due to the difficulty in accessing bone marrow specimens from SLE patients and limited research funding. Additionally, the candidate genes screened were not validated with a large enough sample size, which resulted in a lack of rigor and precision. To address these issues, we plan to expand the sample size and validate the function and mechanism of action of the hub genes in SLE in further experiments.
Conclusion
In summary, a novel circRNA-miRNA-mRNA ceRNA network (circ 0000345- miR-22-3-P-ESR1/SIRT1) was identified in SLE by bioinformatics and comprehensive analysis, which could be served as a possible biomarker or therapeutic target for SLE.
Data availability statement
The original contributions presented in the study are publicly available. This data can be found here: GEO, accession GSE266651, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE266651.
Ethics statement
The studies involving humans were approved by the Ethics Committee of Fujian Provincial Hospital (Approval ID: K-2021-11-019). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by Natural Science Foundation of Fujian, grant no. 2023J011199 and Joint Funds for the Innovation of Science and Technology, Fujian province [grant no. 2023Y9305].
We wish to extend our thanks to the patients and their families.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
ReferencesZucchiDElefanteESchiliròDSignoriniVTrentinFBortoluzziA. One year in review 2022: systemic lupus erythematosus. Clin Exp Rheumatol. (2022) 40:4–14. doi: 10.55563/clinexprheumatol/nolysyLerangKGilboeIMSteinarTDGranJT. Mortality and years of potential life loss in systemic lupus erythematosus: a population-based cohort study. Lupus. (2014) 23:1546–52. doi: 10.1177/0961203314551083, PMID: 25209070CatalinaMDOwenKALabonteACGrammerACLipskyPE. The pathogenesis of systemic lupus erythematosus: harnessing big data to understand the molecular basis of lupus. J Autoimmun. (2020) 110:102359. doi: 10.1016/j.jaut.2019.102359, PMID: 31806421AlaU. Competing endogenous RNAs, non-coding RNAs and diseases: an intertwined story. Cells. (2020) 9:1574. doi: 10.3390/cells9071574, PMID: 32605220EnukaYLauriolaMFeldmanMESas-ChenAUlitskyIYardenY. Circular RNAs are long-lived and display only minimal early alterations in response to a growth factor. Nucleic Acids Res. (2016) 44:1370–83. doi: 10.1093/nar/gkv1367, PMID: 26657629WiluszJESharpPA. Molecular biology. A circuitous route to noncoding RNA. Science. (2013) 340:440–1. doi: 10.1126/science.1238522, PMID: 23620042WangGGuoXChengLChuPChenMChenY. An integrated analysis of the circRNA-miRNA-mRNA network reveals novel insights into potential mechanisms of cell proliferation during liver regeneration. Artif Cells Nanomed Biotechnol. (2019) 47:3873–84. doi: 10.1080/21691401.2019.1669623CesanaMCacchiarelliDLegniniISantiniTSthandierOChinappiM. A long noncoding RNA controls muscle differentiation by functioning as a competing endogenous RNA. Cell. (2011) 147:358–69. doi: 10.1016/j.cell.2011.09.028, PMID: 22000014YuYChenYZhangXLuXHongJGuoX. Knockdown of lncRNA KCNQ1OT1 suppresses the adipogenic and osteogenic differentiation of tendon stem cell via downregulating miR-138 target genes PPARgamma and RUNX2. Cell Cycle. (2018) 17:2374–85. doi: 10.1080/15384101.2018.1534510, PMID: 30321077WeiQHeJWuCMakotaPGaoFLinH. Identification of a circRNA-miRNA-mRNA network to explore the effects of circRNAs on pathogenesis and treatment of systemic lupus erythematosus. Clin Exp Rheumatol. (2023) 41:1068–76. doi: 10.55563/clinexprheumatol/kgm1rn, PMID: 36226610YeHWangXWangLChuXHuXSunL. Full high-throughput sequencing analysis of differences in expression profiles of long noncoding RNAs and their mechanisms of action in systemic lupus erythematosus. Arthritis Res Ther. (2019) 21:70. doi: 10.1186/s13075-019-1853-7, PMID: 30836987BarrettTTroupDBWilhiteSELedouxPRudnevDEvangelistaC. NCBI GEO: mining tens of millions of expression profiles--database and tools update. Nucleic Acids Res. (2007) 35:D760–5. doi: 10.1093/nar/gkl887, PMID: 17099226DennisGJrShermanBTHosackDAYangJGaoWLaneHC. DAVID: database for annotation, visualization, and integrated discovery. Genome Biol. (2003) 4:P3. doi: 10.1186/gb-2003-4-5-p3HochbergMC. Updating the American College of Rheumatology revised criteria for the classification of systemic lupus erythematosus. Arthritis Rheum. (1997) 40:1725. doi: 10.1002/art.1780400928RitchieMEPhipsonBWuDHuYLawCWShiW. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. (2015) 43:e47. doi: 10.1093/nar/gkv007, PMID: 25605792LiJHLiuSZhouHQuLHYangJH. StarBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. (2014) 42:D92–7. doi: 10.1093/nar/gkt1248, PMID: 24297251VejnarCEBlumMZdobnovEM. miRmap web: comprehensive microRNA target prediction online. Nucleic Acids Res. (2013) 41:W165–8. doi: 10.1093/nar/gkt430, PMID: 23716633JohnBEnrightAJAravinATuschlTSanderCMarksDS. Human MicroRNA targets. PLoS Biol. (2004) 2:e363. doi: 10.1371/journal.pbio.0020363, PMID: 15502875WongNWangX. miRDB: an online resource for microRNA target prediction and functional annotations. Nucleic Acids Res. (2015) 43:D146–52. doi: 10.1093/nar/gku1104, PMID: 25378301AgarwalVBellGWNamJWBartelDP. Predicting effective microRNA target sites in mammalian mRNAs. eLife. (2015) 4:e05005. doi: 10.7554/eLife.05005, PMID: 26267216HsuSDLinFMWuWYLiangCHuangWCChanWL. miRTarBase: a database curates experimentally validated microRNA-target interactions. Nucleic Acids Res. (2011) 39:D163–9. doi: 10.1093/nar/gkq1107, PMID: 21071411ShannonPMarkielAOzierOBaligaNSWangJTRamageD. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. (2003) 13:2498–504. doi: 10.1101/gr.1239303, PMID: 14597658ShermanBTHaoMQiuJJiaoXBaselerMWLaneHC. DAVID: a web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res. (2022) 50:W216–21. doi: 10.1093/nar/gkac194OgataHGotoSSatoKFujibuchiWBonoHKanehisaM. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. (1999) 27:29–34. doi: 10.1093/nar/27.1.29, PMID: 9847135NewmanAMLiuCLGreenMRGentlesAJFengWXuY. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. (2015) 12:453–7. doi: 10.1038/nmeth.3337, PMID: 25822800ZhouWYCaiZRLiuJWangDSJuHQXuRH. Circular RNA: metabolism, functions and interactions with proteins. Mol Cancer. (2020) 19:172. doi: 10.1186/s12943-020-01286-3, PMID: 33317550HsiaoKYSunHSTsaiSJ. Circular RNA – new member of noncoding RNA with novel functions. Exp Biol Med (Maywood). (2017) 242:1136–41. doi: 10.1177/1535370217708978LoddeVMurgiaGSimulaERSteriMFlorisMIddaML. Long noncoding RNAs and circular RNAs in autoimmune diseases. Biomol Ther. (2020) 10:1044. doi: 10.3390/biom10071044, PMID: 32674342ChoiMYFloodKBernatskySRamsey-GoldmanRClarkeAE. A review on SLE and malignancy. Best Pract Res Clin Rheumatol. (2017) 31:373–96. doi: 10.1016/j.berh.2017.09.013, PMID: 29224679GanesanRHosNJGutierrezSFischerJStepekJMDagliduE. Salmonella typhimurium disrupts Sirt1/AMPK checkpoint control of mTOR to impair autophagy. PLoS Pathog. (2017) 13:e1006227. doi: 10.1371/journal.ppat.1006227, PMID: 28192515JinMZhangY. Autophagy and autoimmune diseases. Adv Exp Med Biol. (2020) 1207:405–8. doi: 10.1007/978-981-15-4272-5_28SelaUDayanMHershkovizRCahalonLLiderOMozesE. The negative regulators Foxj1 and Foxo3a are up-regulated by a peptide that inhibits systemic lupus erythematosus-associated T cell responses. Eur J Immunol. (2006) 36:2971–80. doi: 10.1002/eji.200636137, PMID: 17051618HallJMCouseJFKorachKS. The multifaceted mechanisms of estradiol and estrogen receptor signaling. J Biol Chem. (2001) 276:36869–72. doi: 10.1074/jbc.R100029200, PMID: 11459850FengFNylanJBanyaiMTatumASilverstonAEGavalchinJ. The induction of the lupus phenotype by estrogen is via an estrogen receptor-αdependent pathway. Clin Immunol. (2012) 134:226–36. doi: 10.1016/j.clim.2009.10.004, PMID: 19926524TaborDEGouldKA. Estrogen receptor alpha promotes lupus in (NZB×NZW) F1 mice in a B cell intrinsic manner. Clin Immunol. (2017) 174:41–52. doi: 10.1016/j.clim.2016.10.011, PMID: 27989899LeeYJShinKSKangSWLeeCKYooBChaHS. Association of the oestrogen receptor a gene polymorphisms with disease onset in systemic lupus erythematosus. Ann Rheum Dis. (2004) 63:1244–9. doi: 10.1136/ard.2003.012583WangJNuiteMMcAlindonTE. Association of estrogen and aromatase gene polymorphisms with systemic lupus erythematosus. Lupus. (2010) 19:734–40. doi: 10.1177/0961203309359517, PMID: 20305046ZhouALiuXXiaTLiFWangJLiJ. Estrogen receptor alpha gene (ESR1) polymorphism and its interaction with smoking and drinking contribute to susceptibility of systemic lupus erythematosus. Immunol Res. (2017) 65:951–6. doi: 10.1007/s12026-017-8935-x, PMID: 28707126AbdouNIRiderVGreenwellCLiXKimlerBF. Fulvestrant (Faslodex), an estrogen selective receptor downregulator, in therapy of women with systemic lupus erythematosus. Clinical, serologic, bone density, and T cell activation marker studies: a double-blind placebo-controlled trial. J Rheumatol. (2008) 35:797–803. PMID: 18381791ConsiglioCRJuliana da SilveiraSMonticieloOAXavierRMBrenolJCTChiesJAB. SIRT1 promoter polymorphisms as clinical modifiers on systemic lupus erythematosus. Mol Biol Rep. (2014) 41:4233–9. doi: 10.1007/s11033-014-3294-3ArbuckleMRMcClainMTRubertoneMVScofieldRHDennisGJJamesJA. Development of autoantibodies before the clinical onset of systemic lupus erythematosus. N Engl J Med. (2003) 349:1526–33. doi: 10.1056/NEJMoa021933GilesBMBoackleSA. Linking complement and anti-dsDNA antibodies in the pathogenesis of systemic lupus erythematosus. Immunol Res. (2013) 55:10–21. doi: 10.1007/s12026-012-8345-z, PMID: 22941560OlivaresDPerez-HernandezJFornerMJPerez-SorianoCTormosMCSaezGT. Urinary levels of sirtuin-1 associated with disease activity in lupus nephritis. Clin Sci. (2018) 132:569–79. doi: 10.1042/CS20171410, PMID: 29440621YangCLiRXuWDHuangAF. Increased levels of sirtuin-1 in systemic lupus erythematosus. Int J Rheum Dis. (2022) 25:869–76. doi: 10.1111/1756-185X.14360LiuHMaXWangXMaXMaoZZhuJ. Hsa_circ_0000345 regulates the cellular development of ASMCs in response to xygenized low -density lipoprotein. J Cell Mol Med. (2020) 24:11849–57. doi: 10.1111/jcmm.15801, PMID: 32865338WeiZRanHYangC. CircRSF1 contributes to endothelial cell growth, migration and tube formation under ox-LDL stress through regulating miR-758/CCND2 axis. Life Sci. (2020) 259:118241. doi: 10.1016/j.lfs.2020.118241JiangCLiHLiuFShiLLiuJLiY. Hsa_circ_0000345 inhibits cell proliferation, migration and invasion of nasopharyngeal carcinoma cells via miR-513a-3p/PTEN axis. J Physiol Sci. (2022) 72:10. doi: 10.1186/s12576-022-00834-4, PMID: 35545766WuXNYeYXNiuJWLiYLiXYouX. Defective PTEN regulation contributes to B cell hyperresponsiveness in systemic lupus erythematosus. Sci Transl Med. (2014) 6:246ra99. doi: 10.1126/scitranslmed.3009131, PMID: 25101889MeiHYLiuJShenXPWuR. A novel circRNA, circRACGAP1, hampers the progression of systemic lupus erythematosus via miR-22-3p-mediated AKT signalling. Autoimmunity. (2022) 55:360–70. doi: 10.1080/08916934.2022.2073590, PMID: 35543435WangWMouSWangLZhangMShaoXFangW. Up-regulation of serum MiR130b-3p level is associated with renal damage in early lupus nephritis. Sci Rep. (2015) 5:12644. doi: 10.1038/srep12644, PMID: 26316103