Bioinformatics analysis of genomic and immune infiltration patterns in autism spectrum disorder

Ru-Qiong Wei; Wen-Liang Guo; Yin-Teng Wu; Raquel Alarcòn Rodrìguez; Marìa Del Mar Requena Mullor; Yu-Chang Gui; Jian-Wen Xu

doi:10.21037/atm-22-4108

Bioinformatics analysis of genomic and immune infiltration patterns in autism spectrum disorder

Ann Transl Med. 2022 Sep;10(18):1013. doi: 10.21037/atm-22-4108.

Authors

Ru-Qiong Wei^#¹, Wen-Liang Guo^#², Yin-Teng Wu², Raquel Alarcòn Rodrìguez³, Marìa Del Mar Requena Mullor³, Yu-Chang Gui¹, Jian-Wen Xu¹

Affiliations

¹ Department of Rehabilitation Medicine, The First Affiliated Hospital of Guangxi Medical University, Nanning, China.
² Guangxi Medical University, Nanning, China.
³ Faculty of Health Sciences, University of Almerìa, Almeria, Spain.

^# Contributed equally.

Abstract

Background: Autism spectrum disorder (ASD) is a specific type of pervasive developmental disorder, and most studies suggest that the onset of autism may be related to genetic and immune factors. The etiology of autism and the underlying molecular events need to be further addressed.

Methods: The ASD-related dataset GSE18123 was downloaded from the Gene Expression Omnibus (GEO) database. Gene set enrichment analysis (GSEA) was used to screen for Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways that may be associated with autism. The top 5,000 genes with an absolute median difference were obtained, and a co-expression network was constructed using weighted correlation network analysis (WGCNA). In addition, GO and KEGG enrichment analyses were performed for genes in the modules most closely related to ASD. Hub genes were found in the significant modules, and the expression values and receiver operating characteristic (ROC) curves of the hub genes were analyzed and validated. Immune cell infiltration in ASD was calculated using the CIBERSORT algorithm, and the relationship between hub genes and immune cells was analyzed. Finally, GSEA was used to explore the potential pathways of hub genes affecting ASD.

Results: The 5,000 DEGs were divided into eight significant modules by WGCNA. The green module was most significantly associated with ASD, and two hub genes [fatty acid-binding protein 2 (FABP2) and Janus kinase 2 (JAK2)] were found. Immune cell infiltration showed that resting dendritic cells and monocytes differed significantly in ASD and healthy individuals. FABP2 was significantly associated with memory B cells and CD8 T cells. JAK2 was significantly associated with monocytes, CD4 activated memory T cells, CD4 resting memory T cells, activated dendritic cells, gamma delta T cells, regulatory T cells (Tregs), CD8 T cells, and naïve CD4 T cells. FABP2 and JAK2 were found to affect multiple pathways of immunity.

Conclusions: FABP2 and JAK2 may influence the immune microenvironment of ASD by regulating immune cells and immune-related pathways and are candidate molecular markers for the development of ASD.

Keywords: Autism spectrum disorder (ASD); CIBERSORT; immune cell infiltration.