Identification of a circRNA-mediated comprehensive ceRNA network in spinal cord injury pathogenesis

RNAs are closely associated with human diseases; however, immune-related genes (IRGs) and their potential regulatory networks in relation to spinal cord injury (SCI) are still poorly understood. Here, we investigated the key IRGs as well as the competing endogenous RNA (ceRNA) mechanisms that are associated with SCI pathogenesis based on microarray datasets and the use of a rat SCI model. Specifically, four independent SCI microarray datasets from Gene Expression Omnibus (GEO) database were analyzed and, thereafter, differentially expressed IRGs were annotated via Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. Furthermore, based on the GEO datasets, differentially expressed RNAs (DERNAs), including DEcircRNAs, DEmiRNAs, and DEmRNAs were identified and interactions between them were also predicted using online databases, and to construct a circular RNA (circRNA) mediated ceRNA network, candidate RNAs were also identified. Furthermore, the support vector machine (SVM) and least absolute shrinkage and selection operator (LASSO) methods were used for the identification of critical DERNAs, while differential gene expression was validated using the GSE20907 dataset. Our results were as follows. In the SCI microarray datasets, 32, 58, and 74 DEIRGs, DEcircRNAs, and DEmiRNAs were identified, respectively. In addition, GO and KEGG analyses showed that the DEIRGs were primarily enriched in neutrophil-mediated immunity and nuclear factor-kappa B (NF-κB) and hypoxia-inducible factor-1 (HIF-1) signaling pathways, and based on LASSO and SVM screening, PLXNB2 was identified as a DEIRG, while hsa_circ_0026646 was identified as the key circRNA, showing a higher SCI expression. Furthermore, our results proved that PLXNB2 and hsa_circ_0026646 were upregulated in SCI, whereas miR-331-3p was downregulated, and, interestingly, similar expression profiles were confirmed using the rat SCI model. Furthermore, fluorescent reporter assay indicated that both hsa_circ_0026646 and PLXNB2 have miR-331-3p target sites, and the ceRNA hypothesis suggested the dysregulation of hsa_circ_0026646, miR-331-3p, and PLXNB2 in SCI. Thus, our results suggested that in SCI pathogenesis, hsa_circ_0026646 correlates with PLXNB2 by targeting miR-331-3p.