Identification of methylation markers for diagnosis of autism spectrum disorder

Autism spectrum disorder (ASD) is a hereditary heterogeneous neurodevelopmental disorder characterized by social and speech dysplasia. We collected the expression profiles of ASD in GSE26415, GSE42133 and GSE123302 from the gene expression omnibus (GEO) database, as well as methylation data of GSE109905. Differentially expressed genes (DEGs) between ASD and controls were obtained by differential expression analysis. Enrichment analysis identified the biological functions and signaling pathways involved by common genes in three groups of DEGs. Protein-protein interaction (PPI) networks were used to identify genes with the highest connectivity as key genes. In addition, we identified methylation markers by associating differentially methylated positions. Key methylation markers were identified using the least absolute shrink and selection operator (LASSO) model. Receiver operating characteristic curves and nomograms were used to identify the diagnostic role of key methylation markers for ASD. A total of 57 common genes were identified in the three groups of DEGs. These genes were mainly enriched in Sphingolipid metabolism and PPAR signaling pathway. In the PPI network, we identified seven key genes with higher connectivity, and used qRT-PCR experiments to verify the expressions. In addition, we identified 31 methylation markers and screened 3 key methylation markers (RUNX2, IMMP2L and MDM2) by LASSO model. Their methylation levels were closely related to the diagnostic effects of ASD. Our analysis identified RUNX2, IMMP2L and MDM2 as possible diagnostic markers for ASD. Identifying different biomarkers and risk genes will contribute to the diagnosis of ASD and the development of new clinical and drug treatments.