Hundreds of genes are implicated in autism spectrum disorder (ASD) but the mechanisms through which they contribute to ASD pathophysiology remain elusive. Here, we analyzed leukocyte transcriptomics from 1-4 year-old male toddlers with ASD or typical development from the general population. We discovered a perturbed gene network that includes genes that are highly expressed during fetal brain development and which is dysregulated in hiPSC-derived neuron models of ASD. Highconfidence ASD risk genes emerge as upstream regulators of the network, and many risk genes may impact the network by modulating RAS/ERK, PI3K/AKT, and WNT/β-catenin signaling pathways. We found that the degree of dysregulation in this network correlated with the severity of ASD symptoms in the toddlers.These results demonstrate how the heterogeneous genetics of ASD may dysregulate a core network to influence brain development at prenatal and very early postnatal ages and, thereby, the severity of later ASD symptoms. that the network is mostly preserved at the co-expression level between ASD leukocytes and prenatal brain.Specifically, the direction of correlations (i.e., positive or negative) in the leukocyte transcriptome of subjects with ASD is mostly preserved in prenatal and early postnatal brain ( Fig. 2d). Importantly, this preservation of co-expression was significantly higher in the DE-ASD network than in the DE-TD network (p-value <10 -16 ; Fig. S6).
rASD genes are associated with the DE-ASD networkWe next analyzed the DE-ASD network in the context of other studies to test the relevance of our DE-ASD network to ASD. Parikshak et al. previously reported gene co-expression modules associated with cortical laminae development during prenatal and early postnatal ages 11 . A subset of these modules show enrichment in rASD genes 11 . We examined the overlap of our leukocyte-derived network with all modules from Parikshak et al 11 . The DE-ASD network preferentially overlapped with rASD gene-enriched modules from that study ( Fig. 2e). This suggests that our DE-ASD network is functionally related to rASD genes during neocortical development. Our DE-ASD network also overlapped with the networks of rASD genes reported in other studies 7,9 , indicating the robustness of the results (Fig. 2e). Intriguingly, the prenatal brain co-expression network of high-confidence rASD genes was more similar to that of ASD leukocytes than TD leukocytes ( Fig. 2f), suggesting that neurodevelopmental transcriptional programs related to rASD genes might be more active in the leukocyte transcriptome of toddlers with ASD than in that of TD toddlers.With the observed overlap patterns, we next tested for enrichment of rASD genes in the DE-ASD network. For this analysis, we considered different rASD gene lists of different size and varying confidence levels (Methods). Surprisingly, this analysis demonstrated that rASD genes are not enriched in the DE-ASD network (p-value >0.19).
The DE-ASD network is enriched for regulatory targets of rASD genesMany high confidence rASD...