To capture the full spectrum of genetic risk for autism, we performed a two-stage analysis of rare de novo and inherited coding variants in 42,607 autism cases, including 35,130 new cases recruited online by SPARK. We identified 60 genes with exome-wide significance (P < 2.5 × 10−6), including five new risk genes (NAV3, ITSN1, MARK2, SCAF1 and HNRNPUL2). The association of NAV3 with autism risk is primarily driven by rare inherited loss-of-function (LoF) variants, with an estimated relative risk of 4, consistent with moderate effect. Autistic individuals with LoF variants in the four moderate-risk genes (NAV3, ITSN1, SCAF1 and HNRNPUL2; n = 95) have less cognitive impairment than 129 autistic individuals with LoF variants in highly penetrant genes (CHD8, SCN2A, ADNP, FOXP1 and SHANK3) (59% vs 88%, P = 1.9 × 10−6). Power calculations suggest that much larger numbers of autism cases are needed to identify additional moderate-risk genes.

Despite the known heritable nature of autism spectrum disorder (ASD), studies have primarily identified risk genes with de novo variants (DNVs). To capture the full spectrum of ASD genetic risk, we performed a two-stage analysis of rare de novo and inherited coding variants in 42,607 ASD cases, including 35,130 new cases recruited online by SPARK. In the first stage, we analyzed 19,843 cases with one or both biological parents and found that known ASD or neurodevelopmental disorder (NDD) risk genes explain nearly 70% of the genetic burden conferred by DNVs. In contrast, less than 20% of genetic risk conferred by rare inherited loss-of-function (LoF) variants are explained by known ASD/NDD genes. We selected 404 genes based on the first stage of analysis and performed a meta-analysis with an additional 22,764 cases and 236,000 population controls. We identified 60 genes with exome-wide significance (p < 2.5e-6), including five new risk genes (NAV3, ITSN1, MARK2, SCAF1, and HNRNPUL2). The association of NAV3 with ASD risk is entirely driven by rare inherited LoFs variants, with an average relative risk of 4, consistent with moderate effect. ASD individuals with LoF variants in the four moderate risk genes (NAV3, ITSN1, SCAF1, and HNRNPUL2, n = 95) have less cognitive impairment compared to 129 ASD individuals with LoF variants in well-established, highly penetrant ASD risk genes (CHD8, SCN2A, ADNP, FOXP1, SHANK3) (59% vs. 88%, p= 1.9e-06) . These findings will guide future gene discovery efforts and suggest that much larger numbers of ASD cases and controls are needed to identify additional genes that confer moderate risk of ASD through rare, inherited variants.

This study is the first genetically-informed investigation of avoidant/restrictive food intake disorder (ARFID), an eating disorder that profoundly impacts quality of life for those affected. ARFID is highly comorbid with autism, and we provide the first estimate of its prevalence in a large and phenotypically diverse autism cohort (a subsample of the SPARK study, N = 5,157 probands). This estimate, 21% (at a balanced accuracy 80%), is at the upper end of previous estimates from studies based on clinical samples, suggesting under-diagnosis and potentially lack of awareness among caretakers and clinicians. Although some studies suggest a decrease of disordered eating symptoms by age 6, our estimates indicate that up to 17% (at a balanced accuracy 87%) of parents of autistic children are also at heightened risk for ARFID, suggesting a lifelong risk for disordered eating. We were also able to provide the first estimates of narrow-sense heritability (h2) for ARFID risk, at 0.45. Genome-wide association revealed a single hit near ZSWIM6, a gene previously implicated in neurodevelopmental conditions. While, the current sample was not well-powered for GWAS, effect size and heritability estimates allowed us to project the sample sizes necessary to more robustly discover ARFID-linked loci via common variants. Further genetic analysis using polygenic risk scores (PRS) affirmed genetic links to autism as well as neuroticism and metabolic syndrome.

Most genetic studies consider autism spectrum disorder (ASD) and developmental disorder (DD) separately despite overwhelming comorbidity and shared genetic etiology. Here, we analyzed de novo variants (DNVs) from 15,560 ASD (6,557 from SPARK) and 31,052 DD trios independently and also combined as broader neurodevelopmental disorders (NDDs) using three models. We identify 615 NDD candidate genes (false discovery rate [FDR] < 0.05) supported by ≥1 models, including 138 reaching Bonferroni exome-wide significance ( P < 3.64e–7) in all models. The genes group into five functional networks associating with different brain developmental lineages based on single-cell nuclei transcriptomic data. We find no evidence for ASD-specific genes in contrast to 18 genes significantly enriched for DD. There are 53 genes that show mutational bias, including enrichments for missense ( n = 41) or truncating ( n = 12) DNVs. We also find 10 genes with evidence of male- or female-bias enrichment, including 4 X chromosome genes with significant female burden ( DDX3X , MECP2 , WDR45 , and HDAC8) . This large-scale integrative analysis identifies candidates and functional subsets of NDD genes.

The complexity of autism’s phenotypic spectra is well-known, yet most genetic research uses case-control status as the target trait. It is undetermined if autistic symptom domain severity underlying this heterogeneity is heritable and pleiotropic with other psychiatric and behavior traits in the same manner as autism case-control status. In N = 6064 autistic children in the SPARK cohort, we investigated the common genetic properties of twelve subscales from three clinical autism instruments measuring autistic traits: the Social Communication Questionnaire (SCQ), the Repetitive Behavior Scale-Revised (RBS-R), and the Developmental Coordination Disorder Questionnaire (DCDQ). Educational attainment polygenic scores (PGS) were significantly negatively correlated with eleven subscales, while ADHD and major depression PGS were positively correlated with ten and eight of the autism subscales, respectively. Loneliness and neuroticism PGS were also positively correlated with many subscales. Significant PGS by sex interactions were found—surprisingly, the autism case-control PGS was negatively correlated in females and had no strong correlation in males. SNP-heritability of the DCDQ subscales ranged from 0.04 to 0.08, RBS-R subscales ranged from 0.09 to 0.24, and SCQ subscales ranged from 0 to 0.12. GWAS in SPARK followed by estimation of polygenic scores (PGS) in the typically-developing ABCD cohort (N = 5285), revealed significant associations of RBS-R subscale PGS with autism-related behavioral traits, with several subscale PGS more strongly correlated than the autism case-control PGS. Overall, our analyses suggest that the clinical autism subscale traits show variability in SNP-heritability, PGS associations, and significant PGS by sex interactions, underscoring the heterogeneity in autistic traits at a genetic level. Furthermore, of the three instruments investigated, the RBS-R shows the greatest evidence of genetic signal in both (1) autistic samples (greater heritability) and (2) general population samples (strongest PGS associations).