Integrative analysis of long non-coding RNAs and messenger RNA expression profiles in systemic lupus erythematosus

Prediction model for developing neuropsychiatric systemic lupus erythematosus in lupus patients

OBJECTIVE

This study aimed to construct a predictive model for assessing the risk of development of neuropsychiatric systemic lupus erythematosus (NPSLE) among patients with SLE based on clinical, laboratory, and meteorological data.

METHODS

A total of 2232 SLE patients were included and were randomly assigned into training and validation sets. Variables such as clinical and laboratory data and local meteorological data were screened by univariate and least absolute shrinkage and selection operator (LASSO) logistic regression modelling. After 10-fold cross-validation, the predictive model was built by multivariate logistic regression, and a nomogram was constructed to visualize the risk of NPSLE. The efficacy and accuracy of the model were assessed by receiver operating characteristic (ROC) curve and calibration curve analysis. Net clinical benefit was assessed by decision curve analysis.

RESULTS

Variables that were included in the predictive model were anti-dsDNA, anti-SSA, lymphocyte count, hematocrit, erythrocyte sedimentation rate, pre-albumin, retinol binding protein, creatine kinase isoenzyme MB, Nterminal brain natriuretic peptide precursor, creatinine, indirect bilirubin, fibrinogen, hypersensitive C-reactive protein, CO, and mild contamination. The nomogram showed a broad prediction spectrum; the area under the curve (AUC) was 0.895 (0.858-0.931) for the training set and 0.849 (0.783-0.916) for the validation set.

CONCLUSION

The model exhibits good predictive performance and will confer clinical benefit in NPSLE risk calculation. Key Points • Clinical, laboratory, and meteorological data were incorporated into a predictive model for neuropsychiatric systemic lupus erythematosus (NPSLE) in SLE patients. • Anti-dsDNA, anti-SSA, LYM, HCT, ESR, hsCRP, IBIL, PA, RBP, CO, Fib, NT-proBNP, Crea, CO, and mild contamination are predictors of the development of NPSLE and may have potential for research. • The nomogram has good predictive performance and clinical value and can be used to guide clinical diagnosis and treatment.

Long noncoding RNAs (CTC-471J1.2, NeST) as epigenetic risk factors of active juvenile lupus nephritis: a case-control study

BACKGROUND

Measurement of the circulating levels of long-non-coding RNAs (lncRNAs) in lupus nephritis (LN) patients could dramatically explore more insights about the disease pathogenesis. Hence, we aimed to quantify the level of expression of CTC-471J1.2 and NeST in LN patients and to correlate it with the disease activity.

METHOD

This case-control study was conducted on a group of children with juvenile LN attending to Mansoura University Children's Hospital (MUCH). Demographics, clinical, and laboratory findings were collected besides the measurement of lncRNAs by quantitative real-time PCR.

RESULTS

The expression level of lncRNAs-CTC-471J1.2 was significantly down-regulated in children with active LN versus inactive cases or controls. In contrast, the NeST was significantly up-regulated in active LN cases. A significant correlation was found between CTC-471J1.2 expression and LN activity parameters. Additionally, both lncRNAs showed a reasonable sensitivity and specificity in differentiation of active LN. A regression analysis model revealed that CTC-471J1.2 and NeST were independent predictors of active nephritis.

CONCLUSION

The expression level of circulatory lncRNAs-CTC-471J1.2 and NeST can be used as sensitive and specific biomarkers for active LN. Furthermore, both could serve as predictors for nephritis activity.

Higher CD4+CD40+ T cells (Th40 cells) associate with systemic lupus erythematosus activity

The aim of this study was to investigate the characteristics of CD4⁺CD40⁺ T cells (Th40 cells) in Chinese systemic lupus erythematosus (SLE) patients. Flow cytometry was used to identify the percentage of Th40 cells in peripheral blood from 24 SLE patients and 24 healthy individuals and the level of IL-2, IL-4, IL-6, IL-10, IFN-r, and TNF-α in serum (22 cases) from the SLE patients. Systemic Lupus Erythematosus Disease Activity Index 2000 (SLEDAI-2000) was used to assess the SLE disease active state. The percentage of Th40 cells in T cells from SLE patients (19.37 ± 17.43) (%) was significantly higher than that from healthy individuals (4.52 ± 3.16) (%) (P < 0.001). The percentage of Th40 cells was also positively associated with SLEDAI-2000 (P = 0.001) and negatively associated with complement C3 (P = 0.007). The Th40 cell percentage was different in SLE patients with different organs involved. The Th40 cell percentage in SLE patients with lupus serositis (29.29 ± 22.19) was significantly higher than that in patients without serositis (13.41 ± 10.79; P = 0.040), and the percentage in SLE patients with lupus pneumonia involvement (29.11 ± 11.88) was significantly higher than that in patients without lupus pneumonia (16.80 ± 17.99; P = 0.043). After 4 weeks treatment, the Th40 cell percentage decreased significantly (P = 0.005). However, Th40 cell expression was not related to cytokines (IL-2, IL-4, IL-6, IL-10, IFN-r, and TNF-α; P > 0.05). A significantly higher percentage of Th40 cells was found in SLE patients, and the Th40 cell percentage was associated with SLE activity. Thus, Th40 cells may be used as a predictor for SLE disease activity and severity and therapeutic efficacy.

PRDX6AS1 gene polymorphisms and SLE susceptibility in Chinese populations

Background

Systemic lupus erythematosus (SLE) is a complex, multisystem autoimmune disease that is characterized by the production of autoantibodies. Although accumulated evidence suggests that the dysregulation of long non-coding RNAs (lncRNAs) is involved in the pathogenesis of SLE, the genetic contributions of lncRNA coding genes to SLE susceptibility remain largely unknown. Here, we aimed to provide more evidence for the role of lncRNA coding genes to SLE susceptibility.

Methods

The genetic association analysis was first adopted from the previous genome-wide association studies (GWAS) and was then validated in an independent cohort. PRDX6-AS1 is located at chr1:173204199-173446294. It spans a region of approximately 240 kb, and 297 single nucleotide polymorphisms (SNPs) were covered by the previous GWAS. Differential expression at the mRNA level was analyzed based on the ArrayExpress Archive database.

Results

A total of 33 SNPs were associated with SLE susceptibility, with a P<1.68×10^-4. The strongest association signal was detected at rs844649 (P=2.12×10^-6), according to the previous GWAS. Combining the results from the GWAS Chinese cohort and our replication cohort, we pursued a meta-analysis approach and found a pronounced genetic association between PRDX6-AS1 rs844649 and SLE susceptibility (p_meta=1.24×10^-13, OR 1.50, 95% CI: 1.34–1.67). The mRNA expression of PRDX6 was elevated in peripheral blood cells, peripheral blood mononuclear cells (PBMCs), and multiple cell subpopulations, such as B cells, CD4⁺ T cells, CD3⁺ cells, and monocytes in patients with SLE. The PRDX6 protein expression level was also increased in patients with SLE compared with healthy donors.

Conclusion

Our study provides new evidence that variants located in lncRNA coding genes are associated with SLE susceptibility.

Systemic lupus erythematosus dysregulates the expression of long noncoding RNAs in placentas

Background

Systemic lupus erythematosus (SLE) can cause placental dysfunctions, which may result in pregnancy complications. Long noncoding RNAs (lncRNAs) are actively involved in the regulation of immune responses during pregnancy. The present study aimed to determine the lncRNA expression profiles in placentas from women with SLE to gain new insights into the underlying molecular mechanisms in SLE pregnancies.

Methods

RNA sequencing (RNA-seq) analysis was performed to identify SLE-dysregulated lncRNAs and mRNAs in placentas from women with SLE and normal full-term (NT) pregnancies. Bioinformatics analysis was conducted to predict the biological functions of these SLE-dysregulated lncRNAs and mRNAs.

Results

RNA-seq analysis identified 52 dysregulated lncRNAs in SLE placentas, including 37 that were upregulated and 15 downregulated. Additional 130 SLE-dysregulated mRNAs were discovered, including 122 upregulated and 8 downregulated. Bioinformatics analysis revealed that SLE-dysregulated genes were associated with biological functions and gene networks, such as regulation of type I interferon-mediated signaling pathway, response to hypoxia, regulation of MAPK (mitogen-activated protein kinase) cascade, response to steroid hormone, complement and coagulation cascades, and Th1 and Th2 cell differentiation.

Conclusions

This is the first report of the lncRNA profiles in placentas from SLE pregnancies. These results suggest that the aberrant expression and the potential regulatory function of lncRNAs in placentas may play comprehensive roles in the pathogenesis of SLE pregnancies. SLE-dysregulated lncRNAs may potentially serve as biomarkers for SLE.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13075-022-02825-7.

Long Non-coding RNAs in Rheumatology

The last decade has seen an enormous increase in long non-coding RNA (lncRNA) research within rheumatology. LncRNAs are arbitrarily classed as non-protein encoding RNA transcripts that exceed 200 nucleotides in length. These transcripts have tissue and cell specific patterns of expression and are implicated in a variety of biological processes. Unsurprisingly, numerous lncRNAs are dysregulated in rheumatoid conditions, correlating with disease activity and cited as potential biomarkers and targets for therapeutic intervention. In this chapter, following an introduction into each condition, we discuss the lncRNAs involved in rheumatoid arthritis, osteoarthritis and systemic lupus erythematosus. These inflammatory joint conditions share several inflammatory signalling pathways and therefore not surprisingly many commonly dysregulated lncRNAs are shared across these conditions. In the interest of translational research only those lncRNAs which are strongly conserved have been addressed. The lncRNAs discussed here have diverse roles in regulating inflammation, proliferation, migration, invasion and apoptosis. Understanding the molecular basis of lncRNA function in rheumatology will be crucial in fully determining the inflammatory mechanisms that drive these conditions.

Linc8986 and linc0597 in plasma are novel biomarkers for systemic lupus erythematosus

Despite increasing evidence that large intergenic non-coding RNAs (lincRNAs) are widely involved in human disease, the role of lincRNAs in the development of systemic lupus erythematosus (SLE) has remained largely elusive. The purpose of the present study was to investigate the expression of three lincRNAs (linc0597, linc8986 and linc7190) in the plasma of patients with SLE and their potential use as biomarkers for the diagnosis and treatment of SLE. Plasma samples were obtained from 54 patients with SLE, 24 patients with rheumatoid arthritis (RA), 24 patients with Sjogren's syndrome (SS) and 22 healthy controls. LincRNA expression levels were measured by reverse transcription-quantitative PCR. Compared with those in the healthy controls, the plasma levels of linc0597 and linc8986 were significantly increased in the patients with SLE (P<0.001), while the difference in the level of linc7190 was not significant (P=0.052). In addition, there was no significant difference in the levels of linc0597 and linc8986 among patients with RA, patients with SS and the healthy controls (P>0.05). Compared with patients with SLE without lupus nephritis (LN), the levels of linc0597 were significantly higher in patients with LN (P=0.044). For linc7190 and linc8986, there was no significant difference between patients with and without LN (P>0.05). Furthermore, complement component 3 (C3) levels were used to evaluate whether the expression of linc8986 and linc0597 is related to the activity of SLE. The results indicated that the levels of linc8986 and linc0597 were negatively correlated with the level of C3 (P<0.001 and P=0.004, respectively). Further analysis suggested that linc0597 and linc8986 were able to specifically identify patients with SLE and that a combination of linc0597 and linc8986 may improve the diagnostic accuracy. Therefore, the plasma levels of linc0597 and linc8986 may be suitable biomarkers for diagnosing SLE.

Integrated Transcriptome Profiling Revealed That Elevated Long Non-Coding RNA-AC007278.2 Expression Repressed CCR7 Transcription in Systemic Lupus Erythematosus

Purpose

Systemic lupus erythematosus (SLE) is a serious autoimmune disease. Its molecular pathogenesis, especially the long non-coding RNA (lncRNA) function, remains unclear. We want to investigate the lncRNA dysregulation profile and their molecular mechanisms in SLE.

Methods

In this study, we analyzed the transcriptome profiles (RNA-seq) of peripheral blood mononuclear cells (PBMCs) from SLE patients and two published transcriptome datasets to explore lncRNA profiles. The differentially expressed lncRNAs were confirmed by quantitative real-time PCR in another set of female patients. We constructed the lncRNA-mRNA regulatory networks by performing weighted gene co-expression network analysis (WGCNA). Dysregulated lncRNA AC007278.2 was repressed by short hairpin RNA (shRNA) in Jurkat cells. Dual-luciferase reporter gene assay was performed to investigate the regulatory mechanism of AC007278.2 on target gene CCR7.

Results

We observed dominant up-regulation of transcripts, including mRNAs and lncRNAs, in SLE patients. By WGCNA method, we identified three modules that were highly related to SLE. We then focused on one lncRNA, AC007278.2, with a T-helper 1 lineage-specific expression pattern. We observed consistently higher AC007278.2 expression in SLE patients. Co-expression network revealed that AC007278.2 participated in the innate immune response and inflammatory bowel disease pathways. By knocking down AC007278.2 expression, we found that AC007278.2 could regulate the expression of inflammatory and cytokine stimulus response-related genes, including CCR7, AZU1, and TNIP3. AC007278.2 inhibits the functional CCR7 promoter to repress its transcription, thereby regulating autoimmunity and follicular T-helper cell differentiation.

Conclusion

In summary, our study indicated the important regulatory role of lncRNAs in SLE. AC007278.2 may be treated as a novel biomarker for SLE diagnosis and treatment.

m6A RNA Methylation in Systemic Autoimmune Diseases-A New Target for Epigenetic-Based Therapy?

The general background of autoimmune diseases is a combination of genetic, epigenetic and environmental factors, that lead to defective immune reactions. This erroneous immune cell activation results in an excessive production of autoantibodies and prolonged inflammation. During recent years epigenetic mechanisms have been extensively studied as potential culprits of autoreactivity. Alike DNA and proteins, also RNA molecules are subjected to an extensive repertoire of chemical modifications. N6-methyladenosine is the most prevalent form of internal mRNA modification in eukaryotic cells and attracts increasing attention due to its contribution to human health and disease. Even though m6A is confirmed as an essential player in immune response, little is known about its role in autoimmunity. Only few data have been published up to date in the field of RNA methylome. Moreover, only selected autoimmune diseases have been studied in respect of m6A role in their pathogenesis. In this review, I attempt to present all available research data regarding m6A alterations in autoimmune disorders and appraise its role as a potential target for epigenetic-based therapies.

An integrated analysis of lncRNA and mRNA expression profiles in the kidneys of mice with lupus nephritis

Long noncoding RNAs (lncRNAs) are persistently expressed and have been described as potential biomarkers and therapeutic targets in various diseases. However, there is limited information regarding lncRNA expression in the tissue of kidney exhibiting lupus nephritis (LN)a serious complication of systemic lupus erythematosus (SLE). In this study, RNA sequencing (RNA-seq) was performed to characterize the lncRNA and mRNA expression in kidney tissues from LN (MRL/lpr) and control mice. We identified 12,979 novel lncRNAs in mouse. The expression profiles of both mRNAs and lncRNAs were differed significantly between LN and control mice. In particular, there were more upregulated lncRNAs and mRNAs than downregulated ones in the kidney tissues of LN mice. However, GO analysis showed that more downregulated genes were enriched in immune and inflammatory response-associated pathways. KEGG analysis showed that both downregulated and upregulated genes were enriched in a number of pathways, including the SLE pathway, and approximately half of these SLE-associated genes encoded inflammatory factors. Moreover, we observed that 2,181 DElncRNAs may have targeted and regulated the expression of 778 mRNAs in LN kidney tissues. The results of this study showed that 11 DElncRNAs targeted and were co-expressed with six immune and SLE-associated genes. qPCR analysis confirmed that lncRNA Gm20513 positively regulated the expression of the SLE-associated gene H2-Aa. In conclusion, the results of our study demonstrates that lncRNAs influence the progression of LN and provide some cues for further study of lncRNAs in LN. These results regarding the lncRNA-mRNAregulatory network may have important value in LN diagnosis and therapy.

Expression profile and diagnostic value of circRNAs in peripheral blood from patients with systemic lupus erythematosus

Circular RNAs (circRNAs) have gained attention due to their performance in disease diagnosis. However, the characteristics of circRNAs in peripheral blood from patients with systemic lupus erythematosus (SLE) remain unknown. Therefore, the aim of the present study was to determine the expression profile and diagnostic potential of circRNAs in peripheral blood from patients with SLE. The global circRNA expression in the peripheral blood of patients with SLE and healthy controls (HCs) was detected using a circRNA microarray. Then, the expression levels of three upregulated circRNAs were selected for further validation by reverse transcription‑quantitative PCR (RT‑qPCR) in a training set. Moreover, the diagnostic value of these circRNAs was assessed by constructing a receiver operating characteristic curve, and then verified in a blind testing set. In total, 1,566 circRNAs were identified to be dysregulated between patients with SLE and HCs (≥2 fold change, P<0.05). Furthermore, the RT‑qPCR results were consistent with the microarray data, in that all three selected circRNAs, hsa_circ_0082688, hsa_circ_0082689 and hsa_circ_0008675, were significantly upregulated in patients with SLE (P<0.05). Results from the training set demonstrated that the combination of hsa_circ_0082688‑hsa_circ_0082689 may provide the most beneficial diagnostic potential. Moreover, the blind test results indicated that the combination model of hsa_circ_0082688‑hsa_circ_0082689 could effectively discriminate between patients with SLE from patients with rheumatoid arthritis and HCs, with a sensitivity of 91.30%, a specificity of 78.57% and an accuracy of 82.28%. Moreover, the combination model of hsa_circ_0082688‑hsa_circ_0082689 + anti‑dsDNA could more effectively discriminated the SLE group from the control groups, with a sensitivity of 95.65%, a specificity of 100.00% and an accuracy of 98.73%. In addition, correlation analysis results suggested that all three circRNAs in patients with SLE did not correlate with the SLE disease activity index. In conclusion, the expression levels of hsa_circ_0082688‑hsa_circ_0082689 may serve as potential biomarkers for SLE diagnosis.

Exploring the Role of Non-Coding RNAs in the Pathophysiology of Systemic Lupus Erythematosus

Systemic lupus erythematosus (SLE) is a chronic immune-related disorder designated by a lack of tolerance to self-antigens and the over-secretion of autoantibodies against several cellular compartments. Although the exact pathophysiology of SLE has not been clarified yet, this disorder has a strong genetic component based on the results of familial aggregation and twin studies. Variation in the expression of non-coding RNAs has been shown to influence both susceptibility to SLE and the clinical course of this disorder. Several long non-coding RNAs (lncRNAs) such as GAS5, MALAT1 and NEAT1 are dysregulated in SLE patients. Moreover, genetic variants within lncRNAs such as SLEAR and linc00513 have been associated with risk of this disorder. The dysregulation of a number of lncRNAs in the peripheral blood of SLE patients has potentiated them as biomarkers for diagnosis, disease activity and therapeutic response. MicroRNAs (miRNAs) have also been shown to affect apoptosis and the function of immune cells. Taken together, there is a compelling rationale for the better understanding of the involvement of these two classes of non-coding RNAs in the pathogenesis of SLE. Clarification of the function of these transcripts has the potential to elucidate the molecular pathophysiology of SLE and provide new opportunities for the development of targeted therapies for this disorder.

Lnc-ITSN1-2, Derived From RNA Sequencing, Correlates With Increased Disease Risk, Activity and Promotes CD4+ T Cell Activation, Proliferation and Th1/Th17 Cell Differentiation by Serving as a ceRNA for IL-23R via Sponging miR-125a in Inflammatory Bowel Disease

Background: This study aimed to investigate long-non-coding RNA (lncRNA) expression profiles and the correlation of lnc-ITSN1-2 expression with disease risk, activity and inflammation, and its influence on CD4⁺ T cell activation, proliferation, and differentiation of inflammatory bowel disease (IBD).

Methods: LncRNA expression profiles were detected in intestinal mucosa samples from six IBD patients and six healthy controls (HCs). Intestinal mucosa and PBMC lnc-ITSN1-2, IL-23R, and inflammatory cytokines were measured in 120 IBD patients [60 Crohn's disease (CD) and 60 ulcerative colitis (UC)] and 30 HCs. Effect of lnc-ITSN1-2 on IBD CD4⁺ T cell activation, proliferation, and differentiation was determined and its regulatory interaction with miR-125a and IL-23R was detected.

Results: Three-hundred-and-nine upregulated and 310 downregulated lncRNAs were identified in IBD patients by RNA-Sequencing, which were enriched in regulating immune and inflammation related pathways. Large-sample qPCR validation disclosed that both intestinal mucosa and PBMC lnc-ITSN1-2 expressions were increased in IBD patients compared to HCs, and presented with good predictive values for IBD risk, especially for active disease conditions, and they positively correlated with disease activity, inflammation cytokines, and IL-23R in IBD patients. Lnc-ITSN1-2 was decreased after infliximab treatment in active-CD patients. Furthermore, lnc-ITSN1-2 promoted IBD CD4⁺ T cell activation and proliferation, and stimulated Th1/Th17 cell differentiation. Multiple rescue experiments disclosed that lnc-ITSN1-2 functioned in IBD CD4⁺ T cells via targeting miR-125a, then positively regulating IL-23R. Luciferase Reporter assay observed that lnc-ITSN1-2 bound miR-125a, and miR-125a bound IL-23R.

Conclusion: Lnc-ITSN1-2 correlates with increased disease risk, activity, and inflammatory cytokines of IBD, and promotes IBD CD4⁺ T cell activation, proliferation, and Th1/Th17 cell differentiation by serving as a competing endogenous RNA for IL-23R via sponging miR-125a.

Comprehensive expression profile of long non-coding RNAs in Peripheral blood mononuclear cells from patients with neuropsychiatric systemic lupus erythematosus

Background

Long noncoding RNAs (lncRNAs) are emerging as critical regulators of gene expression in the immune system, but their impact on neuropsychiatric systemic lupus erythematosus (NPSLE) remains unknown.

Methods

RNA sequencing analysis was used to screen the comprehensive expression profile of lncRNAs and messenger RNAs (mRNAs) in peripheral blood mononuclear cells (PBMCs) from NPSLE patients, active SLE patients who had never experienced neuropsychiatric manifestations (Non-NPSLE) and healthy controls. Differentially expressed (DE) lncRNA levels were validated by qRT-PCR in 26 NPSLE patients, 31 Non-NPSLE patients and 30 healthy controls. Further, correlations of DE lncRNAs with clinical manifestations of NPSLE patients were analyzed. Finally, a bioinformatic analysis was performed to investigate the potential functions of DE genes.

Results

Four hundred and fifty-one lncRNAs and 272 mRNAs were DE between the NPSLE patients and Non-NPSLE patients, among which, significantly upregulated expression levels of NONHSAT208182.1, NONHSAT182114.1, NONHSAT106801.2, NONHSAT039491.2, ENST00000356215, NONHSAT087499.2 and NONHSAT207026.1 while downregulated expression levels of NONHSAT001281.2 and NONHSAT024353.2 were further validated in PBMCs from NPSLE patients by qRT-PCR. Bioinformatic analysis suggested several gene ontology (GO) terms and signal pathways may play important roles in NPSLE development. Co-expression networks analysis indicated that 170 lncRNAs and 46 mRNAs were included in the co-expression network. The expression level of NONHSAT039491.2 was associated with the activity of SLE and the presence of anti-dsDNA, anti-RNP antibody, dizziness and headache. NONHSAT087499.2 level correlated with anti-RNA antibody, ENST00000356215 level correlated with olfactory threshold and oral ulcer. NONHSAT208182.1 level correlated with the presence of fever, unstable walking and urinary red blood cells. NONHSAT106801.2 correlated with frequency of B cells and the presence of fever. NONHSAT024353.2 level was associated with serum IgG levels and the presence of anti-SSA and disorder of consciousness.

Conclusions

Our data provided comprehensive evidence regarding the differential expression of lncRNAs in PBMCs from NPSLE patients, indicating that these DE lncRNAs may play roles in NPSLE development. Our finding shed light on the understanding of the molecular mechanisms of lncRNAs in the pathogenesis of NPSLE.

The epigenetic face of lupus: Focus on antigen-presenting cells

In recent years, epigenetic mechanisms became widely known due to their ability to regulate and maintain physiological processes such as cell growth, development, differentiation and genomic stability. When dysregulated, epigenetic mechanisms, may introduce gene expression changes and disturbance in immune homeostasis leading to autoimmune diseases. Systemic lupus erythematosus (SLE), the most extensively studied autoimmune disorder, has already been correlated with epigenetic modifications, especially in T cells. Since these cell rely on antigen presentation, it may be assumed that erroneous activity of antigen-presenting cells (APCs), culminates in T cell abnormalities. In this review we summarize and discuss the epigenetic modifications in SLE affected APCs, with the focus on dendritic cells (DCs), B cells and monocytes. Unravelling this aspect of SLE pathogenesis, might result in identification of new disease biomarkers and putative therapeutic approaches.

Analysis of lncRNA expression profiles by sequencing reveals that lnc-AL928768.3 and lnc-AC091493.1 are novel biomarkers for disease risk and activity of rheumatoid arthritis

BACKGROUND

This study aimed to analyze long non-coding RNA (lncRNA) expression profiles in synovium tissue of patients with RA using RNA sequencing, and to further assess the clinical values of dysregulated lncRNAs in RA diagnosis and monitoring.

METHODS

Thirty patients with RA who underwent knee arthroscopy and 30 controls with knee trauma who underwent surgery were consecutively enrolled and synovium tissue samples of both groups were obtained during surgery. In the exploration stage, lncRNA and mRNA expression profiles in three RA samples and three control samples were detected by RNA sequencing and bioinformatic analyses were then performed. In the validation stage, quantitative polymerase chain reaction (qPCR) was subsequently used to detect expression of five candidate lncRNAs in 30 patients with RA and 30 control patients.

RESULTS

A total of 349 lncRNAs and 1582 mRNAs were upregulated and 806 lncRNAs and 1295 mRNAs were downregulated in patients with RA compared with controls. Enrichment analyses revealed that these dysregulated lncRNAs and mRNAs were mainly involved in regulating immune response, leukocyte migration, complement activation, and B cell receptor signaling pathway. Subsequent qPCR validation discovered that lnc-AL928768.3 (P < 0.001) and lnc-AC091493.1 (P < 0.001) were elevated in patients with RA compared with controls and afford good predictive values for RA risk by receiver operating characteristic (ROC) curve analysis. Additionally, the two lncRNAs were positively associated with C-reactive protein level and disease activity score in 28 joints (ESR) (all P < 0.05).

CONCLUSION

Analysis of lncRNA expression profiles by sequencing reveals that lnc-AL928768.3 and lnc-AC091493.1 are novel biomarkers for RA risk and activity.

Cross-Talk between Mitochondrial Dysfunction-Provoked Oxidative Stress and Aberrant Noncoding RNA Expression in the Pathogenesis and Pathophysiology of SLE

Systemic lupus erythematosus (SLE) is a prototype of systemic autoimmune disease involving almost every organ. Polygenic predisposition and complicated epigenetic regulations are the upstream factors to elicit its development. Mitochondrial dysfunction-provoked oxidative stress may also play a crucial role in it. Classical epigenetic regulations of gene expression may include DNA methylation/acetylation and histone modification. Recent investigations have revealed that intracellular and extracellular (exosomal) noncoding RNAs (ncRNAs), including microRNAs (miRs), and long noncoding RNAs (lncRNAs), are the key molecules for post-transcriptional regulation of messenger (m)RNA expression. Oxidative and nitrosative stresses originating from mitochondrial dysfunctions could become the pathological biosignatures for increased cell apoptosis/necrosis, nonhyperglycemic metabolic syndrome, multiple neoantigen formation, and immune dysregulation in patients with SLE. Recently, many authors noted that the cross-talk between oxidative stress and ncRNAs can trigger and perpetuate autoimmune reactions in patients with SLE. Intracellular interactions between miR and lncRNAs as well as extracellular exosomal ncRNA communication to and fro between remote cells/tissues via plasma or other body fluids also occur in the body. The urinary exosomal ncRNAs can now represent biosignatures for lupus nephritis. Herein, we’ll briefly review and discuss the cross-talk between excessive oxidative/nitrosative stress induced by mitochondrial dysfunction in tissues/cells and ncRNAs, as well as the prospect of antioxidant therapy in patients with SLE.

Aberrant expression of long non-coding RNAs in peripheral blood mononuclear cells isolated from patients with gouty arthritis

Gouty arthritis (GA) is the most common inflammatory and immune-associated disease, and its prevalence and incidence exhibit yearly increases. The aim of the present study was to analyse the expression profile variation of long non-coding RNAs (lncRNAs) in GA patients and to explore the role of lncRNAs in the pathogenesis of GA. The peripheral blood mononuclear cells of GA patients and of healthy controls (HCs) were used to detect for the differentially expressed lncRNAs by microarray. The functional annotations and classifications of the differentially expressed transcripts were predicted using Gene Ontology (GO) and pathway analysis. The results were then verified by reverse transcription-quantitative (RT-q)PCR. A total of 1,815 lncRNAs and 971 mRNAs with a >2-fold difference in the levels of expression in the GA patients compared with those in the HCs were identified. According to the GO functional enrichment analysis, the differentially expressed lncRNAs were accumulated in terms including protein binding, catalytic activity and molecular transducer activity. The pathways predicted to be involved were the tumor necrosis factor signaling pathway, osteoclast differentiation, NOD-like receptor signaling pathway and NF-κB signaling pathway. The expression of six lncRNAs was measured by RT-qPCR and the results were consistent with those of the microarrays. Among these lncRNAs, AJ227913 was the most differentially expressed lncRNA in GA patients vs. HCs. The expression of several lncRNAs was significantly changed in GA patients compared with that in HCs, which suggests that these lncRNAs with differential expression levels may have an important role in the development and progression of GA.

Effect of astragalosides on long non-coding RNA expression profiles in rats with adjuvant-induced arthritis

Rheumatoid arthritis (RA) is an autoimmune disease of unknown etiology, which occurs in ~1.0% of the general population. Increasing studies have suggested that long non-coding RNAs (lncRNAs) may serve important roles in various biological processes and may be associated with the pathogenesis of different types of disease, including RA. Astragalosides (AST) has been used as a traditional Chinese medicine for the treatment of RA. However, the mechanism underlying its therapeutic effect has remained unclear to date. Thus, there is an urgent need to elucidate the possible mechanism of AST in the treatment of RA from the perspective of lncRNAs. In the present study, the lncRNAs and mRNAs of a vehicle group, animal model group and AST treatment (control) group were determined by Arraystar Rat lncRNA/mRNA microarray. The differentially expressed genes with a fold change >1.5 and P<0.05 were selected and analyzed. Gene Ontology (GO) and pathway analysis was performed using the Database for Annotation, Visualization and Integration Discovery, and the coding-non-coding gene co-expression network was drawn based on the correlation analysis between the differentially expressed lncRNAs and mRNAs. Based on node degree and the correlation between bioinformatics analysis and RA, the critical differentially expressed lncRNAs were selected, analyzed and verified by reverse transcription-quantitative PCR (RT-qPCR) analysis. The results showed that, following AST treatment, up to 75 lncRNAs and 247 mRNAs were found to be differentially expressed among the three groups. GO and pathway analysis manifested that 135 GO terms and 17 pathways were enriched by differentially expressed genes. Four lncRNAs (MRAK012530, MRAK132628, MRAK003448 and XR_006457) were selected as the critical lncRNAs and their trend in expression showed consistency between the RT-qPCR and microarray data. In conclusion, AST had a regulatory effect on differentially expressed lncRNAs during the development of RA, and four lncRNAs could be selected as critical therapeutic targets of AST.

Identification of Renal Long Non-coding RNA RP11-2B6.2 as a Positive Regulator of Type I Interferon Signaling Pathway in Lupus Nephritis

Objective: Lupus nephritis (LN) is one of the most serious complications of systemic lupus erythematosus (SLE). Type I interferon (IFN-I) is associated with the pathogenesis of LN. Long non-coding RNAs (lncRNAs) have been implicated in the pathogenesis of SLE, however, the roles of lncRNAs in LN are still poorly understood. Here, we identified and investigated the function of LN-associated lncRNA RP11-2B6.2 in regulating IFN-I signaling pathway.

Methods: RNA sequencing was used to analyze the expression of lncRNAs in kidney biopsies from LN patients and controls. Antisense oligonucleotides and CRISPRi system or overexpression plasmids and CRISPRa system were used to perform loss or gain of function experiments. In situ hybridization, imaging flow cytometry, dual-luciferase reporter assay, and ATAC sequencing were used to study the functions of lncRNA RP11-2B6.2. RT-qPCR, ELISA, and western blotting were done to detect RNA and protein levels of specific genes.

Results: Elevated lncRNA RP11-2B6.2 was observed in kidney biopsies from LN patients and positively correlated with disease activity and IFN scores. Knockdown of lncRNA RP11-2B6.2 in renal cells inhibited the expression of IFN stimulated genes (ISGs), while overexpression of lncRNA RP11-2B6.2 enhanced ISG expression. Knockdown of LncRNA RP11-2B6.2 inhibited the phosphorylation of JAK1, TYK2, and STAT1 in IFN-I pathway, while promoted the chromatin accessibility and the transcription of SOCS1.

Conclusion: The expression of lncRNAs is abnormal in the kidney of LN. LncRNA RP11-2B6.2 is a novel positive regulator of IFN-I pathway through epigenetic inhibition of SOCS1, which provides a new therapeutic target to alleviate over-activated IFN-I signaling in LN.

Full high-throughput sequencing analysis of differences in expression profiles of long noncoding RNAs and their mechanisms of action in systemic lupus erythematosus

Background

The specific function of long noncoding RNAs (lncRNAs) in systemic lupus erythematosus (SLE) and the mechanism of their involvement in related pathological changes remain to be elucidated, so, in this study, we analyzed the differences in the expression profiles of lncRNAs and their mechanisms of action in SLE using full high-throughput sequencing, bioinformatics, etc. methods.

Methods

We used high-throughput sequencing to detect differences in the expression profiles of lncRNAs, miRNAs, and mRNAs in PBMCs from patients with SLE at the genome-wide level. Next, we predicted target genes of 30 lincRNAs (long intergenic noncoding RNAs) by constructing a coexpression network of differential lincRNAs and mRNAs and identified the role of lincRNAs. Then, we analyzed the coexpression network of 23 optimized lincRNAs and their corresponding 353 miRNAs, evaluated the cis- and trans-effects of these lincRNAs, and performed GO and KEGG analyses of target genes. We also selected 8 lincRNAs and 2 newly discovered lncRNAs for q-PCR validation and lncRNA–miRNA–mRNA analysis. Finally, we also analyzed respectively the relation between lncRNAs and gender bias in SLE patients using RT-qPCR, the relation between Systemic Lupus Erythematosus Disease Activity Index score and the “IFN signature” using ELISA, and the relation between the differential expression of lncRNAs and a change in the number of a cell type of PBMCs in SLE patients using RT-qPCR.

Results

The profiles of 1087 lncRNAs, 102 miRNAs, and 4101 mRNAs in PBMCs significantly differed between patients with SLE and healthy controls. The coexpression network analysis showed that the network contained 23 lincRNAs and 353 mRNAs. The evaluation of the cis- and trans-effects showed that the 23 lincRNAs acted on 704 target genes. GO and KEGG analyses of the target genes predicted the biological functions of the 23 lincRNAs. q-PCR validation showed 7 lincRNAs and 2 novel lncRNAs were identical to the sequencing results. The ceRNA network contained 7 validated lincRNAs, 15 miRNAs, and 155 mRNAs. In addition, the differential expression of lncRNAs may be gender dependent in SLE patients, SLE patients also exhibit a robust “IFN signature,” and PBMCs exhibiting differential expression of lncRNAs may be due to a change in the number of a cell type.

Conclusion

This work determined specific lncRNAs that play important biological functions in the pathogenesis of lupus and provided a new direction for diagnosis and treatment of disease.

Electronic supplementary material

The online version of this article (10.1186/s13075-019-1853-7) contains supplementary material, which is available to authorized users.

Long noncoding RNA expression profile and association with SLEDAI score in monocyte-derived dendritic cells from patients with systematic lupus erythematosus

BACKGROUND

Monocyte-derived dendritic cells (moDCs) play important roles in the pathogenesis of systemic lupus erythematosus (SLE). Aberrant expression of long noncoding RNAs (lncRNAs) could affect the function of moDCs. The aim of this study was to explore the lncRNA expression profile in moDCs of SLE patients to provide new insights into SLE.

METHODS

LncRNA and mRNA microarrays were performed to identify differentially expressed lncRNAs and mRNAs in moDCs of SLE patients compared with normal controls. Bioinformatics analysis was also performed. Quantitative polymerase chain reaction (qPCR) was used to validate the results, and correlation analysis was used to analyze the relationship between these aberrantly expressed lncRNAs and SLE disease activity index (SLEDAI) scores.

RESULTS

According to the gene expression profiles, 163 lncRNAs were differentially expressed between SLE and normal controls, including 118 that were upregulated and 45 that were downregulated. A total of 137 mRNAs were differentially expressed in moDCs of patients with SLE, including 83 that were upregulated and 54 that were downregulated. Furthermore, qPCR data showed that lncRNA ENST00000604411.1 (18.23-fold, P < 0.001) and ENST00000501122.2 (1.96-fold, P < 0.001) were upregulated and the other two lncRNAs, lnc-HSFY2-3:3 (0.42-fold, P < 0.001) and lnc-SERPINB9-1:2 (0.50-fold, P = 0.040), were downregulated in moDCs of SLE patients. The expression levels of ENST00000604411.1 (r = 0.593, P = 0.020) and ENST00000501122.2 (r = 0.539, P = 0.038) were positively correlated with the SLEDAI score, respectively.

CONCLUSIONS

The results indicate that the abnormal expression of lncRNAs in moDCs may be involved in the pathological processes of SLE. The expression level of ENST00000604411.1 and ENST00000501122.2 may have potential value for the assessment of disease activity in SLE.

The Role of Long Non-coding RNAs in the Pathogenesis of RA, SLE, and SS

Rheumatoid diseases are a group of systemic autoimmune diseases which affect multiple organs with largely unknown etiology. In the past decade, long non-coding RNAs (lncRNAs) have emerged as important regulators of biological processes and contribute deeply to immune cell development and immune responses. Substantial evidences have been accumulated showing that LncRNAs involved in the pathogenesis of the rheumatoid diseases, including rheumatoid arthritis (RA), systemic lupus erythematosus (SLE) and Sjögren's syndrome (SS). In this review, we summarize literature combined with bioinformatics methods to analyze the unique and common lncRNAs patterns in rheumatoid diseases and try to reveal the important function of lncRNAs in RA, SLE and SS.