Bioinformatics identification and integrative analysis of ferroptosis-related key lncRNAs in patients with osteoarthritis

Background: Ferroptosis and dysregulation of long non-coding RNA (lncRNA) have been described to be strictly relevant to the pathogenesis of osteoarthritis (OA). However, the connection between ferroptosis and lncRNA in OA is poorly appreciated. Herein, we investigated the functional contribution of lncRNA markers correlated with the progression of human OA by comprehensive bioinformatics analysis of a panoramic network of competing endogenous RNA (ceRNA) based on ferroptosis-related genes (FRGs). Methods: FRGs-related competing endogenous RNA (ceRNA) networks were generated using differentially expressed genes based on OA-related whole transcriptome data from the Gene Expression Omnibus (GEO) database via starBase, miRTarBase, and miRWalk databases. The pivotal lncRNAs were ascertained by topological features (degree, betweenness, and closeness) and subceRNA networks were re-visualized. The expression difference of pivotal lncRNAs was verified by quantitative real-time polymerase chain reaction (qRT-PCR). The latent molecular mechanisms of the global ceRNA and subceRNA networks were uncovered by the R package clusterProfiler-based enrichment analysis. Results: A total of 98 dysregulated lncRNA-miRNA-mRNA regulatory relationships were attained in the FRGs-related panoramic ceRNA network of OA, covering 26 mRNAs, 20 miRNAs, and 20 lncRNAs. Three lncRNAs (AC011511.5, AL358072.1, and C9orf139) were ascertained as the central lncRNAs in the panoramic ceRNA network. Functional ensemble analysis illustrated that both the panoramic ceRNA network and the subceRNA network were integrally affiliated with the immune-inflammatory response, oxygen homeostasis, and cell death (apoptosis, autophagy, and ferroptosis). Conclusion: Comprehensive bioinformatics analysis of the FRGs-related ceRNA network determined three molecular biomarkers of lncRNAs that might be affiliated with OA progression.