Integrative analysis of long noncoding RNAs dysregulation and synapse-associated ceRNA regulatory axes in autism

Abstract

Autism spectrum disorder (ASD) is a complex disorder of neurodevelopment, the function of long noncoding RNA (lncRNA) in ASD remains essentially unknown. In the present study, gene networks were used to explore the ASD disease mechanisms integrating multiple data types (for example, RNA expression, whole-exome sequencing signals, weighted gene co-expression network analysis, and protein-protein interaction) and datasets (five human postmortem datasets). A total of 388 lncRNAs and five co-expression modules were found to be altered in ASD. The downregulated co-expression M4 module was significantly correlated with ASD, enriched with autism susceptibility genes and synaptic signaling. Integrating lncRNAs from the M4 module and microRNA (miRNA) dysregulation data from the literature identified competing endogenous RNA (ceRNA) network. We identified the downregulated mRNAs that interact with miRNAs by the miRTarBase, miRDB, and TargetScan databases. Our analysis reveals that MIR600HG was downregulated in multiple brain tissue datasets and was closely associated with 9 autism-susceptible miRNAs in the ceRNA network. MIR600HG and target mRNAs (EPHA4, MOAP1, MAP3K9, STXBP1, PRKCE, and SCAMP5) were downregulated in the peripheral blood by quantitative reverse transcription polymerase chain reaction analysis (false discovery rate <0.05). Subsequently, we assessed the role of lncRNA dysregulation in altered mRNA levels. Experimental verification showed that some synapse-associated mRNAs were downregulated after the MIR600HG knockdown. BrainSpan project showed that the expression patterns of MIR600HG (primate-specific lncRNA) and synapse-associated mRNA were similar in different human brain regions and at different stages of development. A combination of support vector machine and random forest machine learning algorithms retrieved the marker gene for ASD in the ceRNA network, and the area under the curve of the diagnostic nomogram was 0.851. In conclusion, dysregulation of MIR600HG, a novel specific lncRNA associated with ASD, is responsible for the ASD-associated miRNA-mRNA axes, thereby potentially regulating synaptogenesis.