Construction of an immune-related ceRNA network to screen for potential diagnostic markers for autism spectrum disorder

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

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.

Integrated analysis of endoplasmic reticulum stress regulators' expression identifies distinct subtypes of autism spectrum disorder

Endoplasmic reticulum (ER) stress has been demonstrated to play important roles in a variety of human diseases. However, their relevance to autism spectrum disorder (ASD) remains largely unknown. Herein, we aimed to investigate the expression patterns and potential roles of the ER stress regulators in ASD. The ASD expression profiles GSE111176 and GSE77103 were compiled from the Gene Expression Omnibus (GEO) database. ER stress score determined by the single sample gene set enrichment analysis (ssGSEA) was significantly higher in ASD patients. Differential analysis revealed that there were 37 ER stress regulators dysregulated in ASD. Based on their expression profile, the random forest and artificial neuron network techniques were applied to build a classifier that can effectively distinguish ASD from control samples among independent datasets. Weighted gene co-expression network analysis (WGCNA) screened out the turquoise module with 774 genes was closely related to the ER stress score. Through the overlapping results of the turquoise module and differential expression ER stress genes, hub regulators were gathered. The TF/miRNA-hub gene interaction networks were created. Furthermore, the consensus clustering algorithm was performed to cluster the ASD patients, and there were two ASD subclusters. Each subcluster has unique expression profiles, biological functions, and immunological characteristics. In ASD subcluster 1, the FAS pathway was more enriched, while subcluster 2 had a higher level of plasma cell infiltration as well as the BCR signaling pathway and interleukin receptor reaction reactivity. Finally, the Connectivity map (CMap) database was used to find prospective compounds that target various ASD subclusters. A total of 136 compounds were significantly enriched. In addition to some specific drugs which can effectively reverse the differential gene expression of each subcluster, we found that the PKC inhibitor BRD-K09991945 that targets Glycogen synthase kinase 3β (GSK3B) might have a therapeutic effect on both ASD subtypes that worth of the experimental validation. Our finding proved that ER stress plays a crucial role in the diversity and complexity of ASD, which may inform both mechanistic and therapeutic assessments of the disorder.