Screening Key lncRNAs of Ankylosing Spondylitis Using Bioinformatics Analysis

Abstract

Background: Long non-coding RNAs (lncRNAs) are important regulators in ankylosing spondylitis (AS). Few studies have examined the lncRNA-RNA binding protein (RBP) interaction in AS. This study performed bioinformatics analysis and clinical verification to identify key lncRNAs and propose their RBP interaction. Methods: Three GEO datasets of AS were analyzed by differential expression analysis. The differentially expressed lncRNAs between the AS and control groups were screened out, and the intersecting lncRNAs were regarded as target lncRNAs. Functional was performed to identify target lncRNAs by enrichment analysis, co-expressed RNA analysis, and lncRNA-RBP interaction analysis. Finally, this study analyzed the differential expression level and clinical value of lncRNAs between the AS and control groups. Results: Linc00304, linc00926, and MIAT were differentially expressed and upregulated. Enrichment analysis indicated that the key KEGG terms were the T-cell receptor signaling pathway and B-cell receptor signaling pathway. The key molecular function term was protein binding, and the key biological process term was adaptive immune response. In qRT-PCR results, 44 samples were validated. linc00304 expression was positively correlated with bath ankylosing spondylitis disease activity index (BASDAI), bath ankylosing spondylitis functional index (BASFI), erythrocyte sedimentation rate (ESR), and c-reactive protein (CRP). linc00926 expression was only positively correlated with ESR, whereas MIAT expression was positively correlated with BASFI, ESR, and CRP. Logistic regression revealed that linc00304, ESR, and CRP were the independent risk factors for BASDAI activation. The area under the curve (AUC) of serum linc00304 level in the diagnosis of AS was 0.687 (cutoff value: 0.413, specificity: 0.423, sensitivity: 0.900). AUC of linc00926 was 0.664 (cutoff value: 0.299, sensitivity: 0.882, specificity: 0.417). AUC of MIAT was 0.623 (cutoff value: 0.432, specificity: 0.443, sensitivity: 0.890) (all P <0.05). Conclusion: Overall, this study uncovered three novel lncRNAs, which were upregulated in AS, and proposed a new lncRNA-RBP-mRNA interaction that might regulate adaptive immune response.