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

Thousands of long noncoding RNAs (lncRNAs) have been reported and represent an important subset of pervasive genes associated with a broad range of biological functions. Abnormal expression levels of lncRNAs have been demonstrated in multiple types of human disease. However, the role of lncRNAs in systemic lupus erythematosus (SLE) remains poorly understood. In the present study, the expression patterns of lncRNAs and messenger RNAs (mRNAs) were investigated in peripheral blood mononuclear cells (PBMCs) in SLE using Human lncRNA Array v3.0 (8×60 K; Arraystar, Inc., Rockville, MD, USA). The microarray results indicated that 8,868 lncRNAs (3,657 upregulated and 5,211 downregulated) and 6,876 mRNAs (2,862 upregulated and 4,014 downregulated) were highly differentially expressed in SLE samples compared with the healthy group. Gene ontology (GO) analysis of lncRNA target prediction indicated the presence of 474 matched lncRNA-mRNA pairs for 293 differentially expressed lncRNAs (fold change, ≥3.0) and 381 differentially expressed mRNAs (fold change, ≥3.0). The most enriched pathways were ‘Transcriptional misregulation in cancer’ and ‘Valine, leucine and isoleucine degradation’. Furthermore, reverse transcription-quantitative polymerase chain reaction data verified six abnormal lncRNAs and mRNAs in SLE. The results indicate that the lncRNA expression profile in SLE was significantly changed. In addition, a range of SLE-associated lncRNAs were identified. Thus, the present results provide important insights regarding lncRNAs in the pathogenesis of SLE.