Comprehensive bioinformatics analysis and molecular validation of lncRNAs-mediated ceRNAs network in schizophrenia

The present study aimed to investigate how Schizophrenia (SCZ)-specific long non-coding RNAs (lncRNAs) served as competing endogenous RNAs (ceRNAs) to modulate the biological functions and pathways involved in the pathogenesis of SCZ.

Microarray dataset (GSE54913) was obtained from Gene Expression Omnibus (GEO) database. Differently expressed (DE) lncRNAs and mRNAs were identified by "limma" package. The binding miRNAs of lncRNAs and target mRNAs of shared miRNAs were predicted by miRcode, miRDB, miRTarbase and targetscan databases. Following the ceRNAs theory, interaction network was established and visualized with the cytoscape. Functional enrichment analysis uncovered the concentrated functions and signaling pathways that may be associated with SCZ progression. Protein-protein interaction (PPI) analysis was utilized to determine hub genes. Quantitative real-time PCR (qRT-PCR) and receiver operating characteristic curve (ROC) were performed to evaluate the expression and diagnostic value of ceRNAs members, respectively.

DElncRNAs and DEmRNAs were initially screened from GSE54913 to construct the SCZ-related ceRNAs network with 42 nodes and 53 edges. Functional enrichment analysis revealed that ceRNAs members appeared to be highly correlated with transcription factor activation, cell replication and tumor-related pathways. Once validated, a significant ceRNAs subnetwork was proposed as being implicated in the pathogenesis of SCZ. ROC analysis indicated that SCZ-related ceRNAs members may be sensitive diagnostic biomarkers for SCZ.

The significant SCZ-related ceRNAs subnetworks (lncRNA-C2orf48A/hsa-miR-20b-5p,-17-5p/KIF23, FOXJ2) may represent promising predictive and diagnostic biomarkers and provide novel insights into the mechanism by which lncRNAs act as microRNA sponges and contribute to the pathogenesis of SCZ.