Recent advances in consortium-scale genome-wide association studies (GWAS) have highlighted the involvement of common genetic variants in autism spectrum disorder (ASD), but our understanding of their etiologic roles, especially the interplay with rare variants, is incomplete. In this work, we introduce an analytical framework to quantify the transmission disequilibrium of genetically regulated gene expression from parents to offsprings. We applied this framework to conduct a transcriptome-wide association study (TWAS) on 7,805 ASD proband-parent trios, and replicated our findings using 35,740 independent samples. We identified 31 associations at the transcriptome-wide significance level. In particular, we identified POU3F2 (p=2.1e-7), a transcription factor (TF) mainly expressed in developmental brain. TF targets regulated by POU3F2 showed a 2.1-fold enrichment for known ASD genes (p=4.6e-5) and a 2.7-fold enrichment for loss-of-function de novo mutations in ASD probands (p=7.1e-5). These results provide a clear example of the connection between ASD genes affected by very rare mutations and an unlinked key regulator affected by common genetic variations.
Results
Transmission disequilibrium of polygenic risk, gene expression, and SNP allelesWe applied multiple analytical approaches to dissect common SNPs' contributions to ASD risk at different scales. First, we performed polygenic transmission disequilibrium test (pTDT) 9 to examine the transmission disequilibrium of ASD polygenic risk in probands. ASD polygenic risk scores (PRS) were constructed using case-control samples from the iPSYCH cohort (N=35,740; Methods). We confirmed a highly significant over-transmission of ASD PRS from parents to probands in multiple datasets (p=1.4e-25 in the meta-analysis), including the SPARK cohort which has not been previously analyzed (p=1.0e-11; Supplementary Figure 1). No significant overtransmission was identified in 3,245 healthy siblings (p=0.88).