Genome-wide, integrative analysis implicates circular RNA dysregulation in autism and the corresponding circular RNA-microRNA-mRNA regulatory axes

Circular RNAs (circRNAs), a class of long non-coding RNAs, are known to be enriched in mammalian brain and neural tissues. While the effects of regulatory genetic variants on gene expression in autism spectrum disorder (ASD) have been widely reported, the role of circRNAs in ASD remains largely unknown. Here, we performed genome-wide circRNA expression profiling in post-mortem brains from individuals with ASD and controls and identified 60 circRNAs and three co-regulated modules that were perturbed in ASD. By integrating circRNA, microRNA, and mRNA dysregulation data derived from the same cortex samples, we identified 8,170 ASD-associated circRNA-microRNA-mRNA interactions. Putative targets of the axes were enriched for ASD risk genes and genes encoding inhibitory postsynaptic density (PSD) proteins, but not for genes implicated in monogenetic forms of other brain disorders or genes encoding excitatory PSD proteins. This result reflects the previous observation that ASD-derived organoids exhibit overproduction of inhibitory neurons. We further confirmed that some ASD risk genes (NLGN1, STAG1, HSD11B1, VIP, and UBA6) were indeed regulated by an upregulated circRNA (circARID1A) via sponging a downregulated microRNA (miR-204-3p) in human neuronal cells. We provided a systems-level view of landscape of circRNA regulatory networks in ASD cortex samples. We also provided multiple lines of evidence for the functional role of ASD for circRNA dysregulation and a rich set of ASD-associated circRNA candidates and the corresponding circRNA-miRNA-mRNA axes, particularly those involving ASD risk genes. Our findings thus support a role for circRNA dysregulation and the corresponding circRNA-microRNA-mRNA axes in ASD pathophysiology.