Integrative analysis of transcriptome-wide association study and mRNA expression profiles identifies candidate genes associated with autism spectrum disorders

Autism spectrum disorders (ASD) are a group of highly heritable psychiatric syndromes with high prevalence. The genetic mechanism of ASD remains elusive now. Here we conducted a transcriptome-wide association study (TWAS) of ASD. The GWAS summary data of ASD was driven from the Psychiatric Genomics Consortium (PGC) portal, totally involving 5,305 ASD cases and 5,305 controls. FUSION software was applied to the GWAS summary data for tissue-related TWAS of ASD considering brain and blood. The ASD associated genes identified by TWAS were further validated by mRNA expression profiling of ASD and the Simons Foundation for Autism Research (SFARI) Gene tool. DAVID 6.8 was used to perform gene ontology (GO) enrichment analysis of ASD associated genes identified by TWAS. TWAS identified 85 genes with TWAS P value <0.05 for ASD. Further comparing the 85 genes with the differentially expressed genes identified by mRNA expression profiling of ASD patients found 5 overlapped genes, including MUTYH (P _TWAS = 0.0460, P _mRNA = 0.0040), ARHGAP27 (P _TWAS = 0.0100, P _mRNA = 0.0016), GCA (P _TWAS = 0.0480, P _mRNA = 0.0063), CCDC14 (P _TWAS = 0.0067, P _mRNA = 0.0035), and MED15 (P _TWAS = 0.0324, P _mRNA = 0.0092). Gene Ontology (GO) enrichment analysis of the genes identified by TWAS detected 10 significant GO terms, such as mitochondrion (P = 0.0051), NAD or NADH binding (P = 0.0169), mitochondrial part (P = 0.0386) and 2-oxoglutarate metabolic process (P = 0.0399). In conclusion, this study identified multiple ASD associated genes and gene sets, providing novel clues for revealing the pathogenesis of ASD. Autism Research 2019, 12: 33–38.