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.